Ion-Specific Effects on Ion and Polyelectrolyte Solvation

Ion-specific effects on aqueous solvation of monovalent counter ions, Na+, K+, Cl-, and Br-, and two model polyeletrolytes (PEs), poly(styrene sulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) were here studied with ab initio molecular dynamics (AIMD) and classical molecular dynamics (MD) simulations based on the OPLS-aa force-field which is an empirical fixed point-charge force-field. Ion-specific binding to the PE charge groups was also characterized. Both computational methods predict similar response for the solvation of the PEs but differ notably in description of ion solvation. Notably, AIMD captures the experimentally observed differences in Cl- and Br- anion solvation and binding with the PEs, while the classical MD simulations fail to differentiate the ion species response. Furthermore, the findings show that combining AIMD with the computationally less costly classical MD simulations allows benefiting from both the increased accuracy and statistics reach.

Exploration of Vitamin B6-Based Redox-Active Pyridinium Salts towards the Application in Aqueous Organic Flow Batteries

Pyridoxal hydrochloride, a vitamin B6 vitamer, was synthetically converted to a series of diverse redox-active benzoyl pyridinium salts. Cyclic voltammetry studies demonstrated redox reversibility under basic conditions, and two of the most promising salts were subjected to laboratory-scale redox flow battery tests involving galvanostatic cycling at 10 mM in 0.1 M NaOH. In these tests, the battery was charged completely, corresponding to the transfer of two electrons to the electrolyte, but no discharge was observed. Both CV analysis and electrochemical simulations confirmed that the redox wave observed in the experimental voltammograms corresponds to a two-electron process. To explain the irreversibility in the battery tests, we conducted bulk electrolysis with the benzoyl pyridinium salts, affording the corresponding benzylic secondary alcohols. Computational studies suggest that the reduction proceeds in three consecutive steps: first electron transfer (ET), then proton-coupled electron transfer (PCET) and finally proton transfer (PT) to give the secondary alcohol. 1H NMR deuterium exchange studies indicated that the last PT step is not reversible in 0.1 M NaOH, rendering the entire redox process irreversible. The apparent reversibility observed in CV at the basic media likely arises from the slow rate of the PT step at the timescale of the measurement.

DFT Mechanistic Investigation into Ni(II)-Catalyzed Hydroxylation of Benzene to Phenol by H2O2

Introduction of oxygen into aromatic C-H bonds is intriguing from both fundamental and practical perspectives. Although the 3d metal-catalyzed hydroxylation of arenes by H2O2 has been developed by several prominent researchers, a definitive mechanism for these crucial transformations remains elusive. Herein, density functional theory calculations were used to shed light on the mechanism of the established hydroxylation reaction of benzene with H2O2, catalyzed by [NiII(tepa)]2+ (tepa = tris[2-(pyridin-2-yl)ethyl]amine). Dinickel(III) bis(μ-oxo) species have been proposed as the key intermediate responsible for the benzene hydroxylation reaction. Our findings indicate that while the dinickel dioxygen species can be generated as a stable structure, it cannot serve as an active catalyst in this transformation. The calculations allowed us to unveil an unprecedented mechanism composed of six main steps as follows: (i) deprotonation of coordinated H2O2, (ii) oxidative addition, (iii) water elimination, (iv) benzene addition, (v) ketone generation, and (vi) tautomerization and regeneration of the active catalyst. Addition of benzene to oxygen, which occurs via a radical mechanism, turns out to be the rate-determining step in the overall reaction. This study demonstrates the critical role of Ni-oxyl species in such transformations, highlighting how the unpaired spin density value on oxygen and positive charges on the Ni-O• complex affect the activation barrier for benzene addition.

Kinetics of N2 Release from Diazo Compounds: A Combined Machine Learning-Density Functional Theory Study

Diazo compounds are commonly employed as carbene precursors in carbene transfer reactions during a variety of functionalization procedures. Release of N2 gas from diazo compounds may lead to carbene formation, and the ease of this process is highly dependent on the characteristics of the substituents located in the vicinity of the diazo moiety. A quantum mechanical density functional theory assisted by machine learning was used to investigate the relationship between the chemical features of diazo compounds and the activation energy required for N2 elimination. Our results suggest that diazo molecules, possessing a higher positive partial charge on the carbene carbon and more negative charge on the terminal nitrogen, encounter a lower energy barrier. A more positive C charge decreases the π-donor ability of the carbene lone pair to the π* orbital of N2, while the more negative N charge is a result of a weak interaction between N2 lone pair and vacant p orbital of the carbene. The findings of this study can pave the way for molecular engineering for the purpose of carbene generation, which serves as a crucial intermediate for many chemical transformations in synthetic chemistry.

Density functional theory and machine learning for electrochemical square-scheme prediction: an application to quinone-type molecules relevant to redox flow batteries

Proton-electron transfer (PET) reactions are rather common in chemistry and crucial in energy storage applications. How electrons and protons are involved or which mechanism dominates is strongly molecule and pH dependent. Quantum chemical methods can be used to assess redox potential (Ered.) and acidity constant (pKa) values but the computations are rather time consuming. In this work, supervised machine learning (ML) models are used to predict PET reactions and analyze molecular space. The data for ML have been created by density functional theory (DFT) calculations. Random forest regression models are trained and tested on a dataset that we created. The dataset contains more than 8200 quinone-type organic molecules that each underwent two proton and two electron transfer reactions. Both structural and chemical descriptors are used. The HOMO of the reactant and LUMO of the product participating in the oxidation reaction appeared to be strongly associated with Ered.. Trained models using a SMILES-based structural descriptor can efficiently predict the pKa and Ered. with a mean absolute error of less than 1 and 66 mV, respectively. Good prediction accuracy of R2 > 0.76 and >0.90 was also obtained on the external test set for Ered. and pKa, respectively. This hybrid DFT-ML study can be applied to speed up the screening of quinone-type molecules for energy storage and other applications.

Understanding Electron Transfer Reactions Using Constrained Density Functional Theory: Complications Due to Surface Interactions

The kinetic rates of electrochemical reactions depend on electrodes and molecules in question. In a flow battery, where the electrolyte molecules are charged and discharged on the electrodes, the efficiency of the electron transfer is of crucial importance for the performance of the device. The purpose of this work is to present a systematic atomic-level computational protocol for studying electron transfer between electrolyte and electrode. The computations are done by using constrained density functional theory (CDFT) to ensure that the electron is either on the electrode or in the electrolyte. The ab initio molecular dynamics (AIMD) is used to simulate the movement of the atoms. We use the Marcus theory to predict electron transfer rates and the combined CDFT-AIMD approach to compute the parameters for the Marcus theory where it is needed. We model the electrode with a single layer of graphene and methylviologen, 4,4'-dimethyldiquat, desalted basic red 5, 2-hydroxy-1,4-naphthaquinone, and 1,1-di(2-ethanol)-4,4-bipyridinium were selected for the electrolyte molecules. All of these molecules undergo consecutive electrochemical reactions with one electron being transferred at each stage. Because of significant electrode-molecule interactions, it is not possible to evaluate outer-sphere ET. This theoretical study contributes toward the development of a realistic-level prediction of electron transfer kinetics suitable for energy storage applications.

Temperature-Controlled Syngas Production via Electrochemical CO2 Reduction on a CoTPP/MWCNT Composite in a Flow Cell

The mixture of CO and H₂, known as syngas, is a building block for many substantial chemicals and fuels. Electrochemical reduction of CO₂ and H₂O to syngas would be a promising alternative approach for its synthesis due to negative carbon emission footprint when using renewable energy to power the reaction. Herein, we present temperature-controlled syngas production by electrochemical CO₂ and H₂O reduction on a cobalt tetraphenylporphyrin/multiwalled carbon nanotube (CoTPP/MWCNT) composite in a flow cell in the temperature range of 20-50 °C. The experimental results show that for all the applied potentials the ratio of H₂/CO increases with increasing temperature. Interestingly, at -0.6 V _RHE and 40 °C, the H₂/CO ratio reaches a value of 1.2 which is essential for the synthesis of oxo-alcohols. In addition, at -1.0 V _RHE and 20 °C, the composite shows very high selectivity toward CO formation, reaching a Faradaic efficiency of ca. 98%. This high selectivity of CO formation is investigated by density functional theory modeling which underlines that the potential-induced oxidation states of the CoTPP catalyst play a vital role in the high selectivity of CO production. Furthermore, the stability of the formed intermediate species is evaluated in terms of the pK_a value for further reactions. These experimental and theoretical findings would provide an alternative way for syngas production and help us to understand the mechanism of molecular catalysts in dynamic conditions.

Universal Trends between Acid Dissociation Constants in Protic and Aprotic Solvents

pKa values in non-aqueous solvents are of critical importance in many areas of chemistry. Our knowledge is, despite their relevance, still limited to the most fundamental properties and few pKa values in the most common solvents. Taking advantage of a recently introduced computationally efficient procedure we computed the pKa values of 182 compounds in 21 solvents. This data set is used to establish for the first time universal trends between all solvents. Our computations indicate, that the total charge of the molecule and the charge of the acidic group combined with the Kamlet-Taft solvatochromic parameters are sufficient to predict pKa values with at least semi- quantitative accuracy. We find, that neutral acids such as alcohols are strongly affected by the solvent properties. This is contrasted by cationic acids like ammonium ions whose pKa is often almost completely independent from the choice of solvent.

Reassignment of magic numbers for icosahedral Au clusters: 310, 564, 928 and 1426

Icosahedral Au clusters with three and four shells of atoms are found to deviate significantly from the commonly assumed Mackay structures. By introducing additional atoms in the surface shell and creating a vacancy in the center of the cluster, the calculated energy per atom can be lowered significantly, according to several different descriptions of the interatomic interaction. Analogous icosahedral structures with five and six shells of atoms are generated using the same structural motifs and are similarly found to be more stable than Mackay icosahedra. The lowest energy per atom obtained here is for clusters containing 310, 564, 928 and 1426 atoms, as compared with the commonly assumed magic numbers of 309, 561, 923 and 1415. Some of the vertices in the optimized clusters have a hexagonal ring of atoms, rather than a pentagon, with the vertex atom missing. An inner shell atom in some cases moves outwards by more than an Ångström into the surface shell at such a vertex site. This feature, as well as the wide distribution of nearest-neighbor distances in the surface layer, can strongly influence the properties of icosahedral clusters, for example catalytic activity. The structural optimization is initially carried out using the GOUST method with atomic forces estimated with the EMT empirical potential function, but the atomic coordinates are then refined by minimization using electron density functional theory (DFT) or Gaussian approximation potential (GAP). A single energy barrier is found to separate the Mackay icosahedron from a lower energy structure where a string of atoms moves outwards in a concerted manner from the center so as to create a central vacancy while placing an additional atom in the surface shell.

How to Predict the pK a of Any Compound in Any Solvent

Acid-base properties of molecules in nonaqueous solvents are of critical importance for almost all areas of chemistry. Despite this very high relevance, our knowledge is still mostly limited to the pK _a of rather few compounds in the most common solvents, and a simple yet truly general computational procedure to predict pK _a's of any compound in any solvent is still missing. In this contribution, we describe such a procedure. Our method requires only the experimental pK _a of a reference compound in water and a few standard quantum-chemical calculations. This method is tested through computing the proton solvation energy in 39 solvents and by comparing the pK _a of 142 simple compounds in 12 solvents. Our computations indicate that the method to compute the proton solvation energy is robust with respect to the detailed computational setup and the construction of the solvation model. The unscaled pK _a's computed using an implicit solvation model on the other hand differ significantly from the experimental data. These differences are partly associated with the poor quality of the experimental data and the well-known shortcomings of implicit solvation models. General linear scaling relationships to correct this error are suggested for protic and aprotic media. Using these relationships, the deviations between experiment and computations drop to a level comparable to that observed in water, which highlights the efficiency of our method.

High-Quality Graphene Using Boudouard Reaction

Following the game-changing high-pressure CO (HiPco) process that established the first facile route toward large-scale production of single-walled carbon nanotubes, CO synthesis of cm-sized graphene crystals of ultra-high purity grown during tens of minutes is proposed. The Boudouard reaction serves for the first time to produce individual monolayer structures on the surface of a metal catalyst, thereby providing a chemical vapor deposition technique free from molecular and atomic hydrogen as well as vacuum conditions. This approach facilitates inhibition of the graphene nucleation from the CO/CO2 mixture and maintains a high growth rate of graphene seeds reaching large-scale monocrystals. Unique features of the Boudouard reaction coupled with CO-driven catalyst engineering ensure not only suppression of the second layer growth but also provide a simple and reliable technique for surface cleaning. Aside from being a novel carbon source, carbon monoxide ensures peculiar modification of catalyst and in general opens avenues for breakthrough graphene-catalyst composite production.

Dynamics and Surface Propensity of H+ and OH- within Rigid Interfacial Water: Implications for Electrocatalysis

Facile solvent reorganization promoting ion transfer across the solid-liquid interface is considered a prerequisite for efficient electrocatalysis. We provide first-principles insight into this notion by examining water self-ion dynamics at a highly rigid NaCl(100)-water interface. Through extensive density functional theory molecular dynamics simulations, we demonstrate for both acidic and alkaline solutions that Grotthuss dynamics is not impeded by a rigid water structure. Conversely, decreased proton transfer barriers and a striking propensity of H₃O⁺ and OH^- for stationary interfacial water are found. Differences in the ideal hydration structure of the ions, however, distinguish their behavior at the water contact layer. While hydronium can maintain its optimal solvation, the preferentially hypercoordinated hydroxide is repelled from the immediate vicinity of the surface due to interfacial coordination reduction. This has implications for alkaline hydrogen electrosorption in which the formation of undercoordinated OH^- at the surface is proposed to contribute to the observed sluggish kinetics.

From absolute potentials to a generalized computational standard hydrogen electrode for aqueous and non-aqueous solvents

We describe a simple and efficient procedure to compute a conversion factor for the absolute potential of the standard hydrogen electrode in water to any other solvent. In contrast to earlier methods our procedure only requires the pKa of an arbitrary acid in water and few simple quantum chemical calculations as input. Thus, it is not affected adversely by experimental shortcomings related to measurements in non-aqueous solvents. By combining this conversion factor with the absolute potential in water, the absolute potential in the solvent of interest is obtained. Based on this procedure a new generalized computational standard hydrogen electrode for the computation of electron transfer and proton-coupled electron transfer potentials in non-aqueous solvents and ionic liquids is developed. This enables for the first time the reliable prediction of redox potentials in any solvent. The method is tested through calculation of absolute potentials in 36 solvents. Using the Kamlet-Taft linear solvation energy model we find that the relative absolute potentials consistently increase with decreasing polarisability and decreasing hydrogen bonding ability. For protic solvents good agreement with literature is observed while significant deviations are found for aprotic solvents. The obtained conversion factors are independent of the quantum chemical method, while minor differences are observed between solvation models. This does, however, not affect the global trends.

Hydrogen adsorption trends on two metal-doped Ni2P surfaces for optimal catalyst design

In this study, we looked at the hydrogen evolution reaction on the doubly doped Ni3P2 terminated Ni2P surface. Two Ni atoms in the first three layers of the Ni2P surface model were exchanged with two transition metal atoms. We limited our investigation to combinations of Al, Co, and Fe based on their individual effectiveness as Ni2P dopants in our previous computational studies. The DFT calculated hydrogen adsorption free energy was employed as a predictor of the materials' catalytic HER activity. Our results indicate that the combination of Co and Fe dopants most improves the catalytic activity of the surface through the creation of multiple novel and active catalytic sites.

Method for the accurate prediction of electron transfer potentials using an effective absolute potential

A protocol for the accurate computation of electron transfer (ET) potentials from ab initio and density functional theory (DFT) calculations is described. The method relies on experimental pKa values, which can be measured accurately, to compute a computational setup dependent effective absolute potential. The effective absolute potentials calculated using this protocol display strong variations between the different computational setups and deviate in several cases significantly from the "generally accepted" value of 4.28 V. The most accurate estimate, obtained from CCSD(T)/aug-ccpvqz, indicates an absolute potential of 4.14 V for the normal hydrogen electrode (nhe) in water. Using the effective absolute potential in combination with CCSD(T) and a moderately sized basis, we are able to predict ET potentials accurately for a test set of small organic molecules (σ = 0.13 V). Similarly we find the effective absolute potential method to perform equally good or better for all considered DFT functionals compared to using one of the literature values for the absolute potential. For, M06-2X, which comprises the most accurate DFT method, standard deviation of 0.18 V is obtained. This improved performance is a result of using the most appropriate effective absolute potential for a given method.

Water Oxidation at Neutral pH using a Highly Active Copper-Based Electrocatalyst

The sluggish kinetics of the oxygen evolution reaction (OER) at the anode severely limit hydrogen production at the cathode in water splitting systems. Although electrocatalytic systems based on cheap and earth-abundant copper catalysts have shown promise for water oxidation under basic conditions, only very few examples with high overpotential can be operated under acidic or neutral conditions, even though hydrogen evolution in the latter case is much easier. This work presents an efficient and robust Cu-based molecular catalyst, which self-assembles as a periodic film from its precursors under aqueous conditions on the surface of a glassy carbon electrode. This film catalyzes the OER under neutral conditions with impressively low overpotential. In controlled potential electrolysis, a stable catalytic current of 1.0 mA cm^-2 can be achieved at only 2.0 V (vs. RHE) and no significant decrease in the catalytic current is observed even after prolonged bulk electrolysis. The catalyst displays first-order kinetics and a single site mechanism for water oxidation with a TOF (k_cat ) of 0.6 s^-1 . DFT calculations on of the periodic Cu(TCA)₂ (HTCA=1-mesityl-1H-1,2,3-triazole-4-carboxylic acid) film reveal that TCA defects within the film create Cu^I active sites that provide a low overpotential route for OER, which involves Cu^I , Cu^II -OH, Cu^III =O and Cu^II -OOH intermediates and is enabled at a potential of 1.54 V (vs. RHE), requiring an overpotential of 0.31 V. This corresponds well with an overpotential of approximately 0.29 V obtained experimentally for the grown catalytic film after 100 CV cycles at pH 6. However, to reach a higher current density of 1 mA cm^-2 , an overpotential of 0.72 V is required.

Revisiting the Volmer-Heyrovský mechanism of hydrogen evolution on a nitrogen doped carbon nanotube: constrained molecular dynamics versus the nudged elastic band method

Density functional theory (DFT) based computational electrochemistry has the potential to serve as a tool with predictive power in the rational development and screening of electrocatalysts for renewable energy technologies. It is, however, of paramount importance that simulations are conducted rigorously at a level of theory that is sufficiently accurate in order to obtain physicochemically sensible results. Herein, we present a comparative study of the performance of the static climbing image nudged elastic band method (CI-NEB) vs. DFT based constrained molecular dynamics simulations with thermodynamic integration in estimating activation and reaction (free) energies of the Volmer-Heyrovský mechanism on a nitrogen doped carbon nanotube. Due to cancellation of errors within the CI-NEB calculations, static and dynamic activation barriers are observed to be surprisingly similar, while a substantial decrease in reaction energies is seen upon incorporation of solvent dynamics. This finding is attributed to two competing effects; (1) solvent reorganization that stabilizes the transition and, in particular, the product states with respect to the reactant state and (2) destabilizing entropic contributions due to solvent fluctuations. Our results highlight the importance of explicitly sampling the interfacial solvent dynamics when studying hydrogen evolution at solid-liquid interfaces.

Active site manipulation in MoS2 cluster electrocatalysts by transition metal doping

The development of non-platinum group metal catalysts for the hydrogen evolution reaction (HER) in water electrolyser devices is essential for their widespread and sustainable deployment. In recent years, molybdenum disulfide (MoS2) catalysts have received significant attention as they not only exhibit good electrocatalytic HER activity but also, crucially, acid-stability. However, further performance enhancement is required for these materials to be competitive with Pt and to that end transition metal doping of MoS2 has been explored as a route to further increasing its catalytic activity. In this work, cluster beam deposition was employed to produce controlled cobalt-doped MoS2 clusters (MoS2-Co). We demonstrate that, in contrast to previous observations of performance enhancement in MoS2 resulting from nickel doping (MoS2-Ni), the introduction of Co has a detrimental effect on HER activity. The contrasting behaviours of Ni and Co doping are rationalized by density functional theory (DFT) calculations, which suggest that HER-active surface vacancies are deactivated by combination with Co dopant atoms, whilst their activity is retained, or even partially enhanced, by combination with Ni dopant atoms. Furthermore, the adatom dopant-vacancy combination kinetics appear to be more than three orders of magnitude faster in MoS2-Co than for MoS2-Ni. These findings highlight a fundamental difference in the influence of transition metal dopants on the HER performance of MoS2 electrocatalysts and stress the importance of considering surface atomic defects when predicting their behaviour.

Hydrogen adsorption trends on various metal-doped Ni2P surfaces for optimal catalyst design

This article shows that the presence of dopants, especially Fe and Co, can enhance the HER catalytic activity of Ni₂P.

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.

Diabatic model for electrochemical hydrogen evolution based on constrained DFT configuration interaction

The accuracy of density functional theory (DFT) based kinetic models for electrocatalysis is diminished by spurious electron delocalization effects, which manifest as uncertainties in the predicted values of reaction and activation energies. In this work, we present a constrained DFT (CDFT) approach to alleviate overdelocalization effects in the Volmer-Heyrovsky mechanism of the hydrogen evolution reaction (HER). This method is applied a posteriori to configurations sampled along a reaction path to correct their relative stabilities. Concretely, the first step of this approach involves describing the reaction in terms of a set of diabatic states that are constructed by imposing suitable density constraints on the system. Refined reaction energy profiles are then recovered by performing a configuration interaction (CDFT-CI) calculation within the basis spanned by the diabatic states. After a careful validation of the proposed method, we examined HER catalysis on open-ended carbon nanotubes and discovered that CDFT-CI increased activation energies and decreased reaction energies relative to DFT predictions. We believe that a similar approach could also be adopted to treat overdelocalization effects in other electrocatalytic proton-coupled electron transfer reactions, e.g., in the oxygen reduction reaction.

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 H₂O 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.

Atomic Layer Deposition of Zinc Oxide: Study on the Water Pulse Reactions from First-Principles

Atomic layer deposition (ALD) of zinc oxide thin films has been under intense research in the past few years. The most common precursors used in this process are diethyl zinc (DEZ) and water. The surface chemistry related to the growth of a zinc oxide thin film via atomic layer deposition is not entirely clear, and the ideal model of the process has been contradicted by experimental data, e.g., the incomplete elimination of the ethyl ligands from the surface and the non-negative mass change during the water pulse. In this work we investigate the surface reactions of water during the atomic layer deposition of zinc oxide. The adsorption and ligand-exchange reactions of water are studied on ethyl-saturated surface structures to grasp the relevant surface chemistry contributing to the deposition process. The complex ethyl-saturated surface structures are adopted from a previous publication on the DEZ/H2O-process, and different configurations are sampled using ab initio molecular dynamics in order to find a suitable minimum structure. Water molecules are found to adsorb exothermically onto the ethyl-covered surface at all the ethyl concentrations considered. We do not observe an adsorption barrier for water at 0 K; however, the adsorption energy for any additional water molecules decreases rapidly at high ethyl concentrations. Ligand-exchange reactions are studied at various surface ethyl coverages. The water pulse ligand-exchange reactions have overall larger activation energies than surface reactions for diethyl zinc pulse. For some of the configurations considered, the reaction barriers may be inaccessible at the process conditions, suggesting that some ligands may be inert toward ligand-exchange with water. The activation energies for the surface reactions show only a weak dependence on the surface ethyl concentration. The sensitivity of the adsorption of water at high ethyl coverages suggests that at high ligand-coverages the kinetics may be somewhat hindered due to steric effects. Calculations on the ethyl-covered surfaces are compared to a simple model containing a single monoethyl zinc group. The calculated activation energy for this model is in line with calculations done on the complex model, but the adsorption of water is poorly described. The weak adsorption bond onto a single monoethyl zinc is probably due to a cooperative effect between the surface zinc atoms. A cooperative effect between water molecules is also observed; however, the effect on the activation energies is not as significant as has been reported for other ALD processes.

Hydrogen adsorption trends on Al-doped Ni2P surfaces for optimal catalyst design

Nanoparticles of nickel phosphide are promising materials to replace the currently used rare Pt-group metals at cathode-side electrodes in devices for electrochemical hydrogen production.

Computational exploration of Fe55@C240-catalyzed Fischer–Tropsch synthesis

DFT calculations showed possible hydrocarbon chain growth on Fe55@C240 preferentially via a CO insertion mechanism.

Efficient Constrained Density Functional Theory Implementation for Simulation of Condensed Phase Electron Transfer Reactions

Constrained density functional theory (CDFT) is a versatile tool for probing the kinetics of electron transfer (ET) reactions. In this work, we present a well-scaling parallel CDFT implementation relying on a mixed basis set of Gaussian functions and plane waves, which has been specifically tailored to investigate condensed phase ET reactions using an explicit, quantum chemical representation of the solvent. The accuracy of our implementation is validated against previous theoretical results for predicting electronic couplings and charge transfer energies. Subsequently, we demonstrate the efficiency of our method by studying the intramolecular ET reaction of an organic mixed-valence compound in water using a CDFT based molecular dynamics simulation.

Hydrogen adsorption on MoS2-surfaces: a DFT study on preferential sites and the effect of sulfur and hydrogen coverage

H-Adsorption on MoS₂-surfaces is studied as a function of structural parameters and an assessment of the intricate structure–property relations is conducted.

Isomerism of Trimeric Aluminum Complexes in Aqueous Environments: Exploration via DFT-Based Metadynamics Simulation

The chemistry of aluminum or oxo-aluminum in water is still relatively unknown, although it is the basis for many chemical and industrial processes, including flocculation in water treatment plants. Trimeric species have a predominant role in the formation of the Keggin cations, which are the basic building blocks of aluminum-based chemicals. Despite this, details of the structural evolution of these small solvated clusters and how this is related to the processes leading to the formation of larger aggregates are still an open issue. To address these questions, here, we have applied the metadynamics (MTD) simulation technique [ Barducci , A. ; Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2010 , 1 , 826 - 843 ] with density functional theory-based molecular dynamics to disclose the dynamics and structural conversions of trimeric aluminum complexes in an aqueous environment. The existence of a variety of competing metastable conformations, for example, book-like, cyclic boat, and linear shape conformations, is revealed in the MTD simulation. Furthermore, equilibrium simulations of the various intermediate states encountered along the MTD trajectory are used to assess their (meta)stability, determine the rearrangement of the OH ligands, and discuss the role of the solvating water.

DFT simulations and microkinetic modelling of 1-pentyne hydrogenation on Cu20 model catalysts

Adsorption and dissociation of H2 and hydrogenation of 1-pentyne on neutral and anionic Cu20 clusters have been investigated using the density functional theory and microkinetic modelling. Molecular adsorption of H2 is found to occur strictly at atop sites. The H2 dimer is activated upon adsorption, and the dissociation occurs with moderate energy barriers. The dissociated H atoms reside preferentially on 3-fold face and 2-fold edge sites. Based on these results, the reaction paths leading to the partial and total hydrogenation of 1-pentyne have been studied step-by-step. The results suggest that copper clusters can display selective activity on the hydrogenation of alkyne and alkene molecules. The hydrogenated products are more stable than the corresponding initial reactants following an energetic staircase with the number of added H atoms. Stable semi-hydrogenated intermediates are formed before the partial (1-pentene) and total (pentane) hydrogenation stages of 1-pentyne. The microkinetic model analysis shows that C5H10 is the dominant product. Increasing the reactants (C5H8/H2) ratio enhances the formation of products (C5H10 and C5H12).

Dissociative adsorption of O2 on negatively charged nitrogen-doped single-walled carbon nanotubes: first-principles calculations

Spin unrestricted DFT calculations have been used to study the molecular and dissociative adsorption of O₂ on achiral substitutional nitrogen-doped single-walled carbon nanotubes with and without additional charges.

Charge distribution and Fermi level in bimetallic nanoparticles

Fermi level equilibration driven charge redistribution and electric dipole formation was quantified using a simple nanocapacitor model for bimetallic nanoparticles.

Theoretical Insight into the Hydrogen Evolution Activity of Open-Ended Carbon Nanotubes

Carbon nanotubes (CNTs), while inactive by themselves, are often used as a platform in the search of new catalysts for the hydrogen evolution reaction (HER) by introducing metal nanoparticles or other dopants. Here, we examine the HER activity of pristine open-ended CNTs considering both the effects of chirality and hydrogen coverage using electronic structure calculations. The results indicate that the formation of different 5-ring structures at the end of the CNT introduces surface sites that are highly active toward HER, whereas the activity of traditional 6-ring sites is not greatly altered by tube termination. At fixed hydrogen coverage, the enhanced activity of these sites was attributed to valence orbitals residing close to the highest occupied molecular level facilitating electron transfer to protons.

First principles study of the atomic layer deposition of alumina by TMA-H2O-process

Atomic layer deposition (ALD) is a coating technology used to produce highly uniform thin films. Aluminiumoxide, Al2O3, is mainly deposited using trimethylaluminium (TMA) and water as precursors and is the most studied ALD-process to date. However, only few theoretical studies have been reported in the literature. The surface reaction mechanisms and energetics previously reported focus on a gibbsite-like surface model but a more realistic description of the surface can be achieved when the hydroxylation of the surface is taken into account using dissociatively adsorbed water molecules. The adsorbed water changes the structure of the surface and reaction energetics change considerably when compared to previously studied surface model. Here we have studied the TMA-H2O process using density functional theory on a hydroxylated alumina surface and reproduced the previous results for comparison. Mechanisms and energetics during both the TMA and the subsequent water pulse are presented. TMA is found to adsorb exothermically onto the surface. The reaction barriers for the ligand-exchange reactions between the TMA and the surface hydroxyl groups were found to be much lower compared to previously presented results. TMA dissociation on the surface is predicted to saturate at monomethylaluminium. Barriers for proton diffusion between surface sites are observed to be low. TMA adsorption was also found to be cooperative with the formation of methyl bridges between the adsorbants. The water pulse was studied using single water molecules reacting with the DMA and MMA surface species. Barriers for these reactions were found to reasonable in the process conditions. However, stabilizing interactions amongst water molecules were found to lower the reaction barriers and the dynamical nature of water is predicted to be of importance. It is expected that these calculations can only set an upper limit for the barriers during the water pulse.

Impact of Ga-V Codoping on Interfacial Electron Transfer in Dye-Sensitized TiO2

The improvement of charge transfer between an organic molecule and a semiconductor is an important and challenging goal in the fields of photovoltaics and photocatalysis. In this work, we present a time-dependent density functional theory investigation of the impact of Ga-V codoping of TiO2 on the excited-state electron injection from perylene-3-carboxylic acid. The doping is shown to raise the charge-transfer efficiency for the highest possible surface dye uptake by ∼16%. The strength of the effect depends on the dopant-pair-dye separation, dopant concentration, and distribution of Ga, V atoms in TiO2. The doping of the superficial level turns out to be more favorable than those in the bulk. The changes in electron injection dynamics are attributed to the modification of accepting semiconductor levels and hybridization profile between molecular and semiconductor states.

CO oxidation catalyzed by neutral and anionic Cu20 clusters: relationship between charge and activity

Reactions of CO and O2 on neutral and anionic Cu20 clusters have been investigated by spin-polarized density functional theory. Three reaction mechanisms of CO oxidation are explored: reactions with atomic oxygen (dissociated O2) as well as reactions with molecular oxygen, including Langmuir-Hinshelwood (LH) and Eley-Rideal (ER) mechanisms. The adsorption energies, reaction pathways, and reaction barriers for CO oxidation are calculated systematically. The anionic Cu20(-) cluster can adsorb CO and O2 more strongly than the neutral counterpart due to the superatomic shell closing of 20 valence electrons which leaves one electron above the band gap. The activation of O2 molecule upon adsorption is crucial to determine the rate of CO oxidation. The CO oxidation proceeds efficiently on both Cu20 and Cu20(-) clusters, when O2 is pre-adsorbed dissociatively. The ER mechanism has a lower reaction barrier than the LH mechanism on the neutral Cu20. In general, CO oxidation occurs more readily on the anionic Cu20(-) (effective reaction barriers 0.1-0.3 eV) than on the neutral Cu20 cluster (0.3-0.5 eV). Moreover, Cu20(-) exhibits enhanced binding for CO2. From the analysis of the reverse direction of CO oxidation, it is observed that the transition of CO2 to CO + O can occur on the Cu20(-) cluster, which demonstrates that Cu clusters may serve as good catalyst for CO2 chemistry.

Ab initio Kinetic Monte Carlo simulations of dissolution at the NaCl-water interface

We have used ab initio molecular dynamics (AIMD) simulations to study the interaction of water with the NaCl surface. As expected, we find that water forms several ordered hydration layers, with the first hydration layer having water molecules aligned so that oxygen atoms are on average situated above Na sites. In an attempt to understand the dissolution of NaCl in water, we have then combined AIMD with constrained barrier searches, to calculate the dissolution energetics of Na(+) and Cl(-) ions from terraces, steps, corners and kinks of the (100) surface. We find that the barrier heights show a systematic reduction from the most stable flat terrace sites, through steps to the smallest barriers for corner and kink sites. Generally, the barriers for removal of Na(+) ions are slightly lower than for Cl(-) ions. Finally, we use our calculated barriers in a Kinetic Monte Carlo as a first order model of the dissolution process.

The molecular and magnetic structure of carbon-enclosed and partially covered Fe55 particles

The structure and magnetic moment distribution are studied for an iron nanoparticle with varying degree of carbon adatom coverage. The limiting models of the study are the clean icosahedral Fe55 particle and the iron particle completely enclosed in carbon cages. Between the two extrema, partially covered particles are considered. The iron cluster with partial coverage of carbon adatoms represents a model of active catalysts in the chemical vapor deposition synthesis of carbon nanotubes. The investigated structures are the bare Fe55 cluster, Fe55N4C(x) (x = 27, 37, 47, 54, 65), and Fe55 encapsulated inside C180 and C240. The two latter are extreme examples of an iron particle completely enclosed in a carbon network. Fe55@C180 and Fe55@C240 present novel structures resembling the endohedral metallofullerenes. Two structural isomers of the Fe55@C180 are considered. Enclosing the Fe55 cluster inside C180 and C240 fullerenes gives rise to changes in the Fe-Fe bond lengths. This alters the magnetic structure of the iron cluster considerably. The interaction between the fullerenes and the enclosed iron cluster is reflected in a charge transfer of 8-13 electrons in the considered endohedral complexes. The localization of atomic charges on the C180 and C240 cages suggests site-selective reactivity of the endohedral complexes. The total magnetic moments of the Fe55N4C(x) nanoparticles vary with the degree of adatom coverage. The magnetic moments of individual Fe atoms depend strongly on the element of the nearest-neighbor atoms and on the coordination number and carry therefore information about the local chemistry.

Dissolution of NaCl nanocrystals: an ab initio molecular dynamics study

NaCl nanocrystal dissolution was investigated in atomistic detail revealing a difference in the solvation of two different ionic species.