Kink-Controlled Gold Nanoparticles for Electrochemical Glucose Oxidation.

Enzymes in nature efficiently catalyze chiral organic molecules by elaborately tuning the geometrical arrangement of atoms in the active site. However, enantioselective oxidation of organic molecules by heterogeneous electrocatalysts is challenging because of the difficulty in controlling the asymmetric structures of the active sites on the electrodes. Here, we show that the distribution of chiral kink atoms on high-index facets can be precisely manipulated even on single gold nanoparticles; and this enabled stereoselective oxidation of hydroxyl groups on various sugar molecules. We characterized the crystallographic orientation and the density of kink atoms and investigated their specific interactions with the glucose molecule due to the geometrical structure and surface electrostatic potential.

Anomalous π-backbonding in complexes between B(SiR3)3 and N2: catalytic activation and breaking of scaling relations.

Chemical transformations of molecular nitrogen (N2), including the nitrogen reduction reaction (NRR), are difficult to catalyze because of the weak Lewis basicity of N2. In this study, it is shown that Lewis acids of the types B(SiR3)3 and B(GeR3)3 bind N2 and CO with anomalously short and strong B-N or B-C bonds. B(SiH3)3·N2 has a B-N bond length of 1.48 Å and a complexation enthalpy of -15.9 kcal mol-1 at the M06-2X/jun-cc-pVTZ level. The selective binding enhancement of N2 and CO is due to π-backbonding from Lewis acid to Lewis base, as demonstrated by orbital analysis and density difference plots. The π-backbonding is found to be a consequence of constructive orbital interactions between the diffuse and highly polarizable B-Si and B-Ge bond regions and the π and π* orbitals of N2. This interaction is strengthened by electron donating substituents on Si or Ge. The π-backbonding interaction is predicted to activate N2 for chemical transformation and reduction, as it decreases the electron density and increases the length of the N-N bond. The binding of N2 and CO by the B(SiR3)3 and B(GeR3)3 types of Lewis acids also has a strong σ-bonding contribution. The relatively high σ-bond strength is connected to the highly positive surface electrostatic potential [VS(r)] above the B atom in the tetragonal binding conformation, but the σ-bonding also has a significant coordinate covalent (dative) contribution. Electron withdrawing substituents increase the potential and the σ-bond strength, but favor the binding of regular Lewis acids, such as NH3 and F-, more strongly than binding of N2 and CO. Molecules of the types B(SiR3)3 and B(GeR3)3 are chemically labile and difficult to synthesize. Heterogenous catalysts with the wanted B(Si-)3 or B(Ge-)3 bonding motif may be prepared by boron doping of nanostructured silicon or germanium compounds. B-doped and hydrogenated silicene is found to have promising properties as catalyst for the electrochemical NRR.

Asymmetric Synthesis of Adjacent Tri- and Tetrasubstituted Carbon Stereocenters: Organocatalytic Aldol Reaction of an Hydantoin Surrogate with Azaarene 2-Carbaldehydes.

A bifunctional amine/squaramide catalyst promoted direct aldol addition of an hydantoin surrogate to pyridine 2-carbaldehyde N-oxides to afford adducts bearing two vicinal tertiary/quaternary carbons in high diastereo- and enantioselectivity (d.r. up to >20:1; ee up to 98 %) is reported. Acid hydrolysis of adducts followed by reduction of the N-oxide group yields enantiopure carbinol-tethered quaternary hydantoin-azaarene conjugates with densely functionalized skeletons. DFT studies of the potential energy surface (B3LYP/6-31+G(d)+CPCM (dichloromethane)) of the reaction correlate the activity of different catalysts and support an intramolecular hydrogen-bond-assisted activation of the squaramide moiety in the transition state of the catalytic reaction.

The local electron attachment energy and the electrostatic potential as descriptors of surface-adsorbate interactions.

Two local reactivity descriptors computed by Kohn-Sham density functional theory (DFT) are used to predict and rationalize interactions of nucleophilic molecules (exemplified by CO and H2O) with transition metal (TM) and oxide surfaces. The descriptors are the electrostatic potential, VS(r), and the local electron attachment energy, ES(r), evaluated on surfaces defined by the 0.001 e Bohr-3 isodensity contour. These descriptors have previously shown excellent abilities to predict regioselectivity and rank molecular as well as nanoparticle reactivities and interaction affinities. In this study, we generalize the descriptors to fit into the framework of periodic DFT computations. We also demonstrate their capabilities to predict local surface propensity for interaction with Lewis bases. It is shown that ES(r) and VS(r) can rationalize the interaction behavior of TM oxides and of fcc TM surfaces, including low-index, stepped and kinked surfaces spanning a wide range of interaction sites with varied coordination environments. Broad future applicability in surface science is envisaged for the descriptors, including heterogeneous catalysis and electrochemistry.

σ-Holes and σ-lumps direct the Lewis basic and acidic interactions of noble metal nanoparticles: introducing regium bonds.

Using local DFT-based probes for electrostatic as well as charge transfer/polarization interactions, we are able to characterize Lewis basic and acidic sites on copper, silver and gold nanoparticles. The predictions obtained using the DFT-probes are compared to the interaction energies of the electron donating (CO, H2O, NH3 and H2S) and the electron accepting (BH3, BF3, HCl [H-down] and Na+) compounds. The probes include the local electron attachment energy [E(r)], the average local ionization energy [I(r)], and the electrostatic potential [V(r)] and are evaluated on isodensity surfaces located at distances corresponding to typical interaction distances. These probes have previously been successful in characterizing molecular interactions. Good correlations are found between Lewis acidity and maxima in V(r), appearing as a consequence of σ-holes, as well as minima in E(r), of the noble metal nanoparticles. Similarly are Lewis basic sites successfully described by surface minima in V(r) and I(r); the former are indicative of σ-lumps, i.e. regions of enhanced σ-density. The investigated probes are anticipated to function as reliable tools in nanoparticle reactivity and interaction characterization, and may act as suitable descriptors in large-scale screenings for materials of specific properties, e.g. in heterogeneous catalysis. Because of the similarity between the noble metal nanoparticle's interactions with Lewis bases and the concepts of halogen and hydrogen bonding, a new class of bonds is introduced - regium bonds - taking place between a σ-hole of a Cu, Ag or Au compound and an electron donor.

Theoretical Investigation into Rate-Determining Factors in Electrophilic Aromatic Halogenation.

The halogenation of monosubstituted benzenes in aqueous solvent was studied using density functional theory at the PCM-M06-2 X/6-311G(d,p) level. The reaction with Cl2 begins with the formation of C atom coordinated π-complex and is followed by the formation of the σ-complex, which is rate-determining. The final part proceeds via the abstraction of the proton by a water molecule or a weak base. We evaluated the use of the σ-complex as a model for the rate-determining transition state (TS) and found that this model is more accurate the later the TS comes along the reaction coordinate. This explains the higher accuracy of the model for halogenations (late TS) compared to nitrations (early TS); that is, the more deactivated the substrate the later the TS. The halogenation with Br2 proceeds with a similar mechanism as the corresponding chlorination, but the bromination has a very late rate-determining TS that is similar to the σ-complex in energy. The iodination with ICl follows a different mechanism than chlorination and bromination. After the formation of the π-complex, the reaction proceeds in a concerted manner without a σ-complex. This reaction has a large primary hydrogen kinetic isotope effect in agreement with experimental observations.

Extending the σ-Hole Concept to Metals: An Electrostatic Interpretation of the Effects of Nanostructure in Gold and Platinum Catalysis.

Crystalline surfaces of gold are chemically inert, whereas nanoparticles of gold are excellent catalysts for many reactions. The catalytic properties of nanostructured gold have been connected to increased binding affinities of reactant molecules for low-coordinated Au atoms. Here we show that the high reactivity at these sites is a consequence of the formation of σ-holes, i.e., maxima in the surface electrostatic potential (VS,max), due to the overlap of mainly the valence s-orbitals when forming the bonding σ-orbitals. The σ-holes are binding sites for Lewis bases, and binding energies correlate with the magnitudes of the VS,max. For symmetrical Au clusters, of varying sizes, the most positive VS,max values are found at the corners, edges, and surfaces (facets), decreasing in that order. This is in agreement with the experimentally and theoretically observed dependence of catalytic activity on local structure. The density of σ-holes can explain the increasing catalytic activity with decreasing particle size for other transition metal catalysts also, such as platinum.

Nucleophilic Aromatic Substitution Reactions Described by the Local Electron Attachment Energy.

A local multiorbital electrophilicity descriptor, the local electron attachment energy [E(r)], is used to study the nucleophilic aromatic substitution reactions of SNAr and VNS (vicarious nucleophilic substitution). E(r) considers all virtual orbitals below the free electron limit and is determined on the molecular isodensity contour of 0.004 atomic units. Good (R2 = 0.83) to excellent (R2 = 0.98) correlations are found between descriptor values and experimental reactivity data for six series of electron deficient arenes. These include homo- and heteroarenes, rings of five to six atoms, and a variety of fluorine, bromine, and hydride leaving groups. The solvent, temperature, and nucleophile are in addition varied across the series. The surface E(r) [ES(r)] is shown to provide reactivity predictions better than those of transition-state calculations for a concerted SNAr reaction with a bromine nucleofug, gives correlations substantially stronger than those of LUMO energies, and is overall more reliable than the molecular electrostatic potential. With the use of ES(r), one can identify the various electrophilic sites within a molecule and correctly predict isomeric distributions. Since the calculations of ES(r) are computationally inexpensive, the descriptor offers fast but accurate reactivity predictions for the important nucleophilic aromatic substitution class of reactions. Applications in, e.g., drug discovery, synthesis, and toxicology studies are envisaged.

Dehydrogenation of methanol on Cu2O(100) and (111).

Adsorption and desorption of methanol on the (111) and (100) surfaces of Cu2O have been studied using high-resolution photoelectron spectroscopy in the temperature range 120-620 K, in combination with density functional theory calculations and sum frequency generation spectroscopy. The bare (100) surface exhibits a (3,0; 1,1) reconstruction but restructures during the adsorption process into a Cu-dimer geometry stabilized by methoxy and hydrogen binding in Cu-bridge sites. During the restructuring process, oxygen atoms from the bulk that can host hydrogen appear on the surface. Heating transforms methoxy to formaldehyde, but further dehydrogenation is limited by the stability of the surface and the limited access to surface oxygen. The (√3 × âˆš3)R30°-reconstructed (111) surface is based on ordered surface oxygen and copper ions and vacancies, which offers a palette of adsorption and reaction sites. Already at 140 K, a mixed layer of methoxy, formaldehyde, and CHxOy is formed. Heating to room temperature leaves OCH and CHx. Thus both CH-bond breaking and CO-scission are active on this surface at low temperature. The higher ability to dehydrogenate methanol on (111) compared to (100) is explained by the multitude of adsorption sites and, in particular, the availability of surface oxygen.

Dual Lewis Acid/Lewis Base Catalyzed Acylcyanation of Aldehydes: A Mechanistic Study.

A mechanistic investigation, which included a Hammett correlation analysis, evaluation of the effect of variation of catalyst composition, and low-temperature NMR spectroscopy studies, of the Lewis acid-Lewis base catalyzed addition of acetyl cyanide to prochiral aldehydes provides support for a reaction route that involves Lewis base activation of the acyl cyanide with formation of a potent acylating agent and cyanide ion. The cyanide ion adds to the carbonyl group of the Lewis acid activated aldehyde. O-Acylation by the acylated Lewis base to form the final cyanohydrin ester occurs prior to decomplexation from titanium. For less reactive aldehydes, the addition of cyanide is the rate-determining step, whereas, for more reactive, electron-deficient aldehydes, cyanide addition is rapid and reversible and is followed by rate-limiting acylation. The resting state of the catalyst lies outside the catalytic cycle and is believed to be a monomeric titanium complex with two alcoholate ligands, which only slowly converts into the product.

Local Electron Attachment Energy and Its Use for Predicting Nucleophilic Reactions and Halogen Bonding.

A new local property, the local electron attachment energy [E(r)], is introduced and is demonstrated to be a useful guide to predict intermolecular interactions and chemical reactivity. The E(r) is analogous to the average local ionization energy but indicates susceptibility toward interactions with nucleophiles rather than electrophiles. The functional form E(r) is motivated based on Janak's theorem and the piecewise linear energy dependence of electron addition to atomic and molecular systems. Within the generalized Kohn-Sham method (GKS-DFT), only the virtual orbitals with negative eigenvalues contribute to E(r). In the present study, E(r) has been computed from orbitals obtained from GKS-DFT computations with a hybrid exchange-correlation functional. It is shown that E(r) computed on a molecular isodensity surface, ES(r), reflects the regioselectivity and relative reactivity for nucleophilic aromatic substitution, nucleophilic addition to activated double bonds, and formation of halogen bonds. Good to excellent correlations between experimental or theoretical measures of interaction strengths and minima in ES(r) (ES,min) are demonstrated.

Searching for the thermodynamic limit--a DFT study of the step-wise water oxidation of the bipyramidal Cu7 cluster.

Oxidative degradation of copper in aqueous environments is a major concern in areas such as catalysis, electronics and construction engineering. A particular challenge is to systematically investigate the details of this process for non-ideal copper surfaces and particles under the conditions found in most real applications. To this end, we have used hybrid density functional theory to study the oxidation of a Cu7 cluster in water solution. Especially, the role of a large water coverage is explored. This has resulted in the conclusion that, under atmospheric H2 pressures, the thermodynamically most favored state of degradation is achieved upon the generation of four H2 molecules (i.e. Cu7 + 8H2O → Cu7(OH)8 + 4H2) in both condensed and gas phases. This state corresponds to an average oxidation state below Cu(I). The calculations suggest that the oxidation reaction is slow at ambient temperatures with the water dissociation as the rate-limiting step. Our findings are expected to have implication for, among other areas, the copper catalyzed water-gas shift reaction, and for the general understanding of copper corrosion in aqueous environments.

Interrupted SNAr-Alkylation Dearomatization.

Dearomatizations provide powerful synthetic routes to rapidly assemble substituted carbocycles and heterocycles found in a plethora of bioactive molecules. Harnessing the advantages of dearomatization typically requires vigorous reagents because of the difficulty in disrupting the stable aromatic core. A relatively mild dearomatization strategy is described that employs lithiated nitriles or isocyanides in a simple SNAr-type addition to form σ-complexes that are trapped by alkylation. The dearomatizations are diastereoselective and efficient and rapidly install two new carbon-carbon bonds, one of which is a quaternary center, as well as nitrile, isocyanide, and cyclohexadiene functionalities.

Control of the ambident reactivity of the nitrite ion.

In previous studies, it was reported that a neighbouring equatorial ester group is essential for a good yield of nitrite-mediated triflate inversion, whereas with neighbouring benzyl ether groups or axial ester groups, mixtures are generally produced. In the present study, the origin of this difference was addressed. The ambident reactivity of the nitrite ion has been found to be the cause of the complex product formation observed, which can be controlled by a neighbouring equatorial ester group. Both N-attack and O-attack occur in the absence of the ester group, whereas O-attack is favoured in its presence. A neighbouring group assistance mechanism is proposed, in addition to steric effects, based on secondary interactions between the neighbouring ester group and the incoming nucleophile. High-level quantum mechanical calculations were carried out in order to delineate this effect. The theoretical results are in excellent agreement with experiments, and suggest a catalytic role for the neighbouring equatorial ester group.

On the method-dependence of transition state asynchronicity in Diels-Alder reactions.

This work discusses the dependence of transition state geometries on the choice of quantum chemical optimization method for the extensively studied Diels-Alder reaction. Rather significant differences are observed between post-Hartree-Fock methods and (hybrid) density functional theory, where the latter predicts larger asynchronicities. The results show that the low MP2 asynchronicity observed is likely artificial. Still, there are significant discrepancies between hybrid and pure density functionals. The role of the exchange functional seems to be most prominent in less activated reacting systems, while the importance of the correlation functional seems to increase as they become more activated by, e.g., an electron-donating group on the diene. To correct the dubious MP2 geometries, we employed the SCS-MP2 protocol for transition state optimization, which leads to significantly better results with respect to CCSD/6-31+G(d) level calculations. We conclude that in order for hybrid functionals to give descriptions consistent with the sample post-Hartree-Fock methods, a balanced combination of both Hartree-Fock exchange (with a couple of exceptions) and a well-behaved correlation functional is required. Given that the benchmark CCSD/6-31+G(d) geometries are sufficient representations, the best geometries were obtained using ωB97X(D), B2PLYP(D) and M06-2X.

Reactivity of metal oxide clusters with hydrogen peroxide and water--a DFT study evaluating the performance of different exchange-correlation functionals.

We have performed a density functional theory (DFT) investigation of the interactions of H2O2, H2O and HO radicals with clusters of ZrO2, TiO2 and Y2O3. Different modes of H2O adsorption onto the clusters were studied. In almost all the cases the dissociative adsorption is more exothermic than molecular adsorption. At the surfaces where H2O has undergone dissociative adsorption, the adsorption of H2O2 and the transition state for its decomposition are mediated by hydrogen bonding with the surface HO groups. Using the functionals B3LYP, B3LYP-D and M06 with clusters of 26 and 8 units of ZrO2, the M06 functional performed better than B3LYP in describing the reaction of decomposition of H2O2 and the adsorption of H2O. Additionally, we investigated clusters of the type (ZrO2)2, (TiO2)2 and (Y2O3) and the performance of the functionals B3LYP, B3LYP-D, B3LYP*, M06, M06-L, PBE0, PBE and PWPW91 in describing H2O2, H2O and HO˙ adsorption and the energy barrier for decomposition of H2O2. The trends obtained for HO˙ adsorption onto the clusters are discussed in terms of the ionization energy of the metal cation present in the oxide. In order to correctly account for the existence of an energy barrier for the decomposition of H2O2, the functional used must include Hartree-Fock exchange. Using minimal cluster models, the best performance in describing the energy barrier for H2O2 decomposition was obtained with the M06 and PBE0 functionals - the average absolute deviations from experiments are 6 kJ mol(-1) and 5 kJ mol(-1) respectively. With the M06 functional and a larger monoclinic (ZrO2)8 cluster model, the performance is in excellent agreement with experimental data. For the different oxides, PBE0 was found to be the most effective functional in terms of performance and computational time cost.

Utilizing the σ-complex stability for quantifying reactivity in nucleophilic substitution of aromatic fluorides.

A computational approach using density functional theory to compute the energies of the possible σ-complex reaction intermediates, the "σ-complex approach", has been shown to be very useful in predicting regioselectivity, in electrophilic as well as nucleophilic aromatic substitution. In this article we give a short overview of the background for these investigations and the general requirements for predictive reactivity models for the pharmaceutical industry. We also present new results regarding the reaction rates and regioselectivities in nucleophilic substitution of fluorinated aromatics. They were rationalized by investigating linear correlations between experimental rate constants (k) from the literature with a theoretical quantity, which we call the sigma stability (SS). The SS is the energy change associated with formation of the intermediate σ-complex by attachment of the nucleophile to the aromatic ring. The correlations, which include both neutral (NH3) and anionic (MeO(-)) nucleophiles are quite satisfactory (r = 0.93 to r = 0.99), and SS is thus useful for quantifying both global (substrate) and local (positional) reactivity in SNAr reactions of fluorinated aromatic substrates. A mechanistic analysis shows that the geometric structure of the σ-complex resembles the rate-limiting transition state and that this provides a rationale for the observed correlations between the SS and the reaction rate.

Stepwise Diels-Alder: more than just an oddity? A computational mechanistic study.

We have employed hybrid DFT and SCS-MP2 calculations at the SMD-PCM-6-311++G(2d,2p)//6-31+G(d) level to investigate the relationship between three possible channels for forming a Diels-Alder adduct from a highly nucleophilic diene and moderately to highly electrophilic dienophiles. We discuss geometries optimized using the B3LYP and M06-2X functionals with the 6-31+(d) basis set. The transition states and intermediates are characterized on the basis of geometric and electronic properties, and we also address the possibility of predicting detectability of a zwitterionic intermediate based on its relative stability. Our results show that a conventional Diels-Alder transition state conformation yields intermediates in all four investigated cases, but that these are too short-lived to be detected experimentally for the less activated reactants. The stepwise trans pathway, beginning with a conjugate addition-like transition state, becomes increasingly competitive with more activated reactants and is indeed favored for the most electrophilic dienophiles. Addition of a trans diene leads to a dead-end as the trans intermediates have insurmountable rotation barriers that prohibit formation of the second bond, unless another, heterocyclic intermediate is formed. We also show that introduction of a hydrogen bond donating catalyst favors a stepwise pathway even for less activated dienophiles.

Predicting regioselectivity in nucleophilic aromatic substitution.

We have investigated practical and computationally efficient methods for the quantitative prediction of regioisomer distribution in kinetically controlled nucleophilic aromatic substitution reactions. One of the methods is based on calculating the relative stabilities of the isomeric σ-complex intermediates using DFT. We show that predictions from this method can be used quantitatively both for anionic nucleophiles with F(-) as leaving group, as well as for neutral nucleophiles with HF as leaving group. The σ-complex approach failed when the leaving group was Cl/HCl or Br/HBr, both for anionic and neutral nucleophiles, because of difficulties in finding relevant σ-complex structures. An approach where we assumed a concerted substitution step and used such transition state structures gave quantitatively useful results. Our results are consistent with other theoretical works, where a stable σ-complex has been identified in some cases, whereas others have been indicated to proceed via a concerted substitution step.

Computational design of a Diels-Alderase from a thermophilic esterase: the importance of dynamics.

A novel computational Diels-Alderase design, based on a relatively rare form of carboxylesterase from Geobacillus stearothermophilus, is presented and theoretically evaluated. The structure was found by mining the PDB for a suitable oxyanion hole-containing structure, followed by a combinatorial approach to find suitable substrates and rational mutations. Four lead designs were selected and thoroughly modeled to obtain realistic estimates of substrate binding and prearrangement. Molecular dynamics simulations and DFT calculations were used to optimize and estimate binding affinity and activation energies. A large quantum chemical model was used to capture the salient interactions in the crucial transition state (TS). Our quantitative estimation of kinetic parameters was validated against four experimentally characterized Diels-Alderases with good results. The final designs in this work are predicted to have rate enhancements of ≈ 10(3)-10(6) and high predicted proficiencies. This work emphasizes the importance of considering protein dynamics in the design approach, and provides a quantitative estimate of the how the TS stabilization observed in most de novo and redesigned enzymes is decreased compared to a minimal, 'ideal' model. The presented design is highly interesting for further optimization and applications since it is based on a thermophilic enzyme (T (opt) = 70 °C).