The application of density functional theory (DFT) methods in catalysis has been growing fast in the last few decades thanks to both the availability of more powerful high computing resources and the development of new efficient approximations and approaches. DFT calculations allow for the understanding of crucial catalytic aspects that are difficult or even impossible to access by experiments, thus contributing to faster development of more efficient and selective catalysts. Depending on the catalytic system and properties under investigation, different approaches should be used. Moreover, the reliability of the obtained results deeply depends on the approximations involved in both the selected method and model. This review addresses chemists, physicists and materials scientists whose interest deals with the application of DFT-based computational tools in both homogeneous catalysis and heterogeneous catalysis. First, a brief introduction to DFT is presented. Then, the main approaches based on atomic centered basis sets and plane waves are discussed, underlining the main differences, advantages and limitations. Eventually, guidance towards the selection of the catalytic model is given, with a final focus on the evaluation of the energy barriers, which represents a crucial step in all catalytic processes. Overall, the review represents a rational and practical guide for both beginners and more experienced users involved in the wide field of catalysis.
The chemical shift (CS) values obtained by 1H NMR spectroscopy for the hydrogen atoms of a tetradentate N2O2-substituted Salphen ligand (H2L1) are differently shifted in its complexes of nickel(II), palladium(II), platinum(II), and zinc(II), all bearing the same charge on the metal ions. To rationalize the observed trends, DFT calculations have been performed in the implicit d6-DMSO solvent in terms of the electronic effects induced by the metal ion and of the nature and strength of the metal-N and metal-O bonds. Overall, the results obtained point out that, in the complexes involving group 10 elements, the CS values show the greater shift when considering the two hydrogen atoms at a shorter distance from the coordinated metal center and follow the decreasing metal charge in the order Ni > Pd > Pt. This trend suggests a more covalent character of the ligand-metal bonds with the increase of the metal atomic number. Furthermore, a slightly poorer agreement between experimental and calculated data is observed in the presence of the nickel(II) ion. Such discrepancy is explained by the formation of stacked oligomers, aimed at minimizing the repulsive interactions with the polar DMSO solvent.
Electron-deficient acridones and in situ generated acridinium salts are reported as potent, closed-shell photooxidants that undergo surprising mechanisms. When bridging acyclic triarylamine catalysts with a carbonyl group (acridones), this completely diverts their behavior away from open-shell, radical cationic, 'beyond diffusion' photocatalysis to closed-shell, neutral, diffusion-controlled photocatalysis. Brønsted acid activation of acridones dramatically increases excited state oxidation power (by +0.8â V). Upon reduction of protonated acridones, they transform to electron-deficient acridinium salts as even more potent photooxidants (*E1/2 =+2.56-3.05â V vs SCE). These oxidize even electron-deficient arenes where conventional acridinium salt photooxidants have thusfar been limited to electron-rich arenes. Surprisingly, upon photoexcitation these electron-deficient acridinium salts appear to undergo two electron reductive quenching to form acridinide anions, spectroscopically-detected as their protonated forms. This new behaviour is partly enabled by a catalyst preassembly with the arene, and contrasts to conventional SET reductive quenching of acridinium salts. Critically, this study illustrates how redox active chromophoric molecules initially considered photocatalysts can transform during the reaction to catalytically active species with completely different redox and spectroscopic properties.
We prepared a series of phenothiazine (PTZ)-anthraquinone (AQ) electron donor-acceptor dyads to study the relationship between molecular structures and the possibility of charge transfer (CT) and intersystem crossing (ISC). As compared to the previously reported PTZ-AQ dyad with a direct connection of two units via a C-N single bond, the PTZ and AQ units are connected via a p-phenylene or p-biphenylene linker. Conformation restriction is imposed by attaching ortho-methyl groups on the phenylene linker. UV-vis absorption spectra indicate electronic coupling between the PTZ and AQ units in the dyads without conformation restriction. Different from the previously reported PTZ-AQ, thermally activated delayed fluorescence (TADF) is observed for the dyads containing one phenylene linker (PTZ-Ph-AQ and PTZ-PhMe-AQ). The prompt fluorescence lifetime in cyclohexane is exceptionally long (τPF = 62.0 ns, population ratio: 99.2%) and 245.0 ns (93.5%) for PTZ-Ph-AQ and PTZ-PhMe-AQ, respectively (normally τPF <20 ns); the delayed fluorescence lifetimes for these two dyads were determined as τDF = 2.4 µs (6.5%) and 7.6 µs (0.8%), respectively. For the dyad containing a biphenylene linker (PTZ-Ph2Me-AQ), no TADF was observed. Charge-separated (CS) states were observed for PTZ-Ph-AQ and PTZ-PhMe-AQ, and the lifetimes were determined as 7.0 and 1.3 µs, respectively, indicating the triplet spin multiplicity of the CS state. The 3CS state lifetimes are shortened to 100 ns and 440 ns for the two dyads, respectively, in the polar solvent acetonitrile. For dyads with a longer linker, i.e., PTZ-Ph2Me-AQ, the CS state lifetime is not sensitive to solvent polarity (τCS = 1.8 and 1.3 µs in cyclohexane and acetonitrile, respectively). In reference dyads, where the PTZ unit is oxidized to sulfoxide, no CT absorption band and TADF were observed, which is attributed to the increased CS state energy (>3 eV) becoming higher than that of the AQ triplet (3AQ*) state (ca. 2.7 eV). These experimental evidence show that the presence of 1CS, 3CS, and 3LE (LE: locally excited) states sharing similar energy is essential for the occurrence of TADF. Population of the long-lived 3CS state (with a lifetime of a few µs) does not produce by itself TADF, because the ISC process of 1CSâ3CS is nonsufficient. Femtosecond transient absorption spectra show that charge separation (CS) occurs readily (<5 ps) for most dyads, even in nonpolar solvents. Nanosecond pulsed laser-excited time-resolved electron paramagnetic resonance (TREPR) spectra show that either a spin correlated radical pair (SCRP) is formed, with the electron exchange energy 2J = +2.14 mT, or radical pairs with stronger interaction, |2J| > 6.57 mT. These studies are useful for in-depth understanding of the CS and ISC in compact electron donor-acceptor dyads and for design of efficient TADF emitters.
The utilization of CO2 as building block for the production of cyclic carbonate is a promising route to simultaneously mitigate the global warming issue and obtain valuable commercial chemicals. In this work, the activity of nicotinamidium halide catalysts towards the CO2 conversion into cyclic carbonate has been explored by means of density functional theory (DFT) calculations. DFT calculations support the ability, suggested experimentally, of the pyridium α-C-H proton of the catalysts to activate the epoxide ring via a hydrogen bond. Interestingly, DFT calculations underline the involvement of the n-octyl substituent of the pyridyl ring in the epoxide activation, while the hydrogen atom of the amide group N-H is rather involved in the stabilization of the iodide trough electrostatic interactions. Moreover, the replacement of the pyridium α-C-H proton with the bulkier methyl group leads to a different reaction mechanism. The calculated energy barriers well reproduce the experimental trends of the studied catalysts, and the computed activation barrier of 29.0â kcal/mol, relative to the ring opening step of the most active catalyst, is in line with the experimental working temperature of 80 °C. Those results shed light on the CO2 fixation reaction contributing to the development of more efficient catalytic systems.
Titanium dioxide (TiO2) has been extensively studied as a suitable material for a wide range of fields including catalysis and sensing. For example, TiO2-based nanoparticles are active in the catalytic conversion of glucose into value-added chemicals, while the good biocompatibility of titania allows for its application in innovative biosensing devices for glucose detection. A key process for efficient and selective biosensors and catalysts is the interaction and binding mode between the analyte and the sensor/catalyst surface. The relevant features regard both the molecular recognition event and its effects on the nanoparticle electronic structure. In this work, we address both these features by combining two first-principles methods based on periodic boundary conditions and cluster approaches (CAs). While the former allows for the investigation of extended materials and surfaces, CAs focus only on a local region of the surface but allow for using hybrid functionals with low computational cost, leading to a highly accurate description of electronic properties. Moreover, the CA is suitable for the study of reaction mechanisms and charged systems, which can be cumbersome with PBC. Here, a direct and detailed comparison of the two computational methodologies is applied for the investigation of d-glucose on the TiO2 (100) anatase surface. As an alternative to the commonly used PBC calculations, the CA is successfully exploited to characterize the formation of surface and subsurface oxygen vacancies and to determine their decisive role in d-glucose adsorption. The results of such direct comparison allow for the selection of an efficient, finite-size structural model that is suitable for future investigations of biosensor electrocatalytic processes and biomass conversion catalysis.
We provide a comprehensive DFT investigation of the mechanistic details of CO2 fixation into styrene oxide to form styrene carbonate, catalyzed by potassium iodide-tetraethylene glycol complex. A detailed view on the intermediate steps of the overall reaction clarifies the role of hydroxyl substances as co-catalysts for the alkali halide-catalyzed cycloaddition. The increase of iodide nucleophilicity in presence of tetraethylene glycol is examined and rationalized by NBO and Hirshfeld charge analysis, and bond distances. We explore how different alkali metal salts and glycols affect the catalytic performance. Our results provide important hints on the synthesis of cyclic carbonates from CO2 and epoxides promoted by alkali halides and glycol complexes, allowing the development of more efficient catalysts.
Due to the high oxygen evolution reaction (OER) activity, stability, and abundance of NiOx materials, they are found to be promising catalysts, competitive with expensive metal oxides such as IrO2 and RuO2. From a theoretical point of view, studies reported in the literature so far are mostly based on density functional theory using periodic slab models for the bulk and surface of ß-NiOOH, one of the active NiOx phases. However, cluster models are a valid method to investigate many aspects about structure, charge carrier transport properties, and OER activity of ß-NiOOH. Hence, here we present new cluster models for the surface of ß-NiOOH, where the oxygen atoms are bonded to Mg effective core potentials (ECPs) mimicking neighboring atom cores. This cluster embedding procedure is superior to saturating the cluster with hydrogen atoms, and to using other atomic ECPs for ß-NiOOH. We find that layered materials such as ß-NiOOH are more vulnerable to geometrical rupture and therefore a cluster approach requires additional care in choosing the embedding approach. We evaluated the models by using them to calculate the energy required for water adsorption and deprotonation, which are essential ingredients for OER. Specifically, our results agree with previous slab models that the first deprotonation reaction step requires a large amount of energy. In addition, we find that water and hydroxyl groups have high adsorption energy and therefore the first deprotonation step is limiting the reaction efficiency.
Mg and Ca ß-diketiminato silylamides [HC{(Me)CN(2,6-iPr2C6H3)}2M(THF)n{N(SiMe3)2}] (M = Mg, n = 0; M = Ca, n = 1) were studied as precatalysts for the dehydrogenation/dehydrocoupling of secondary amine-boranes R2HNBH3 . By reaction with equimolar quantities of amine-boranes, the corresponding amidoborane derivatives are formed, which further react to yield dehydrogenation products such as the cyclic dimer [BH2-NMe2]2. DFT was used here to explore the mechanistic alternatives proposed on the basis of the experimental findings for both Mg and Ca amidoboranes. The influence of the steric demand of amine-boranes on the course of the reaction was examined by performing calculations on the dehydrogenation of dimethylamine-borane (DMAB), pyrrolidine-borane (PB), and diisopropylamine-borane. In spite of the analogies in the catalytic activity of Mg- and Ca-based complexes in the dehydrocoupling of amine-boranes, our theoretical analysis confirmed the experimentally observed lower reactivity of Ca complexes. Differences in catalytic activity of Mg- and Ca-based complexes were examined and rationalized. As a consequence of the increase in ionic radius on going from Mg(2+) to Ca(2+), the dehydrogenation mechanism changes and formation of a key metal hydride intermediate becomes inaccessible. Dimerization is likely to occur off-metal in solution for DMAB and PB, whereas steric hindrance of iPr2NHBH3 hampers formation of the cyclic dimer. The reported results are of particular interest because, although amine-borane dehydrogenation is now well established, mechanistic insight is still lacking for many systems.
A detailed theoretical investigation of the charge transport mechanism in poly(4-vinyl-imidazole) (P4VI), the parent polymer of a series of N-heterocyclic-based membranes used as an electrolyte in proton exchange membrane fuel cells, is presented. In particular, Density Functional Theory (DFT) results obtained for small model systems (protonated imidazole dimers and trimers) suggest that the commonly accepted conduction mechanism, based on a sequential proton transfer between imidazole moieties, could be impeded by the geometrical constraints imposed by the polymeric backbone. Indeed only one kind of proton transfer reaction is energetically allowed between adjacent imidazoles, so that a rotation of the protonated imidazole is required for a second proton transfer. Molecular dynamics simulations on a larger model (15 oligomers with an excess proton) show that the rotation of the imidazole carrying the excess proton is a soft large amplitude motion. These results allow us to propose a new proton conduction mechanism in P4VI, where a frustrated rotation of the protonated imidazole before each proton transfer reaction represents the rate-limiting step. Furthermore, in contrast with the Grotthuss proton transport mechanism in water, our results indicate that here it is the same proton which could be successively transferred. From a chemical point of view, these new insights into the mechanism are relevant for a rational design of modified azole-based systems for Proton Exchange Membrane Fuel Cells.
The recently synthesized rhodium complex [Rh{P(C(5)H(9))(2)(η(2)-C(5)H(7))}(Me(2)HNBH(3))(2)]BAr(F)(4) (2), which incorporates two amine-boranes coordinated to the rhodium center with two different binding modes, namely η(1) and η(2), has been used to probe whether bis(σ-amine-borane) motifs are important in determining the general course of amine-boranes dehydrocoupling reactions. DFT calculations have been carried out to explore mechanistic alternatives that ultimately lead to the formation of the amine-borane cyclic dimer [BH(2)NMe(2)](2) (A) by hydrogen elimination. Sequential concerted, on- or off-metal, intramolecular dehydrogenations provide two coordinated amine-borane molecules. Subsequent dimerization is likely to occur off the metal in solution. In spite of the computationally confirmed presence of a BHâ â â NH hydrogen bond between amine-borane ligands, neither a simple intermolecular route for dehydrocoupling of complex 2 is operating, nor seems [Rh{P(C(5)H(9))(2)(η(2)-C(5)H(7))}B](+) to be important for the whole dehydrocoupling process.
The platination processes of DNA bases with second- and third-generation Pt(II) anticancer drugs have been investigated using density functional theory (DFT) combined with the conductor-like dielectric continuum model (CPCM) approach, in order to describe their binding mechanisms and to obtain detailed data on the reaction energy profiles. Although there is no doubt that a Pt-N7 bond forms during initial attack, the energetic profiles for the formation of the monofunctional adducts are not known. Herein, a direct comparison between the rate of formation of the monofunctional adducts of the second- and third-generation anticancer drugs with guanine (G) and adenine (A) DNA bases has been made in order to spotlight possible common or different behavior. The guanine as target for platination process is confirmed to be preferred over adenine for all the investigated compounds and for both the hydrolyzed forms considered in our investigation. The preference for G purine base is dominated by electronic factors and promoted by a more favorable hydrogen-bonds pattern, confirming the important role played by H-bonds in determining both structural and kinetic control on the purine platination process.
Adenine/metabolism , Antineoplastic Agents/metabolism , DNA Adducts/chemistry , DNA Adducts/metabolism , Guanine/metabolism , Organoplatinum Compounds/metabolism , Carboplatin/metabolism , Hydrogen-Ion Concentration , Models, Molecular , Oxaliplatin , Quantum Theory , Thermodynamics
