Max Planck Institute for Chemical Energy Conversion

This short update provides an overview of the capabilities that have been added to the ORCA electronic structure package (version 4.0) since publication of the first article in 2012. WIREs Comput Mol Sci 2018, 8:e1327. doi: 10.1002/wcms.1327 This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Electronic Structure Theory > Density Functional Theory Software > Quantum Chemistry

CO2 conversion covers a wide range of possible application areas from fuels to bulk and commodity chemicals and even to specialty products with biological activity such as pharmaceuticals. In the present review, we discuss selected examples in these areas in a combined analysis of the state-of-the-art of synthetic methodologies and processes with their life cycle assessment. Thereby, we attempted to assess the potential to reduce the environmental footprint in these application fields relative to the current petrochemical value chain. This analysis and discussion differs significantly from a viewpoint on CO2 utilization as a measure for global CO2 mitigation. Whereas the latter focuses on reducing the end-of-pipe problem “CO2 emissions” from todays’ industries, the approach taken here tries to identify opportunities by exploiting a novel feedstock that avoids the utilization of fossil resource in transition toward more sustainable future production. Thus, the motivation to develop CO2-based chemistry does not depend primarily on the absolute amount of CO2 emissions that can be remediated by a single technology. Rather, CO2-based chemistry is stimulated by the significance of the relative improvement in carbon balance and other critical factors defining the environmental impact of chemical production in all relevant sectors in accord with the principles of green chemistry.

ADVERTISEMENT RETURN TO ISSUEPREVReviewNEXTHydrogenasesWolfgang Lubitz*, Hideaki Ogata, Olaf Rüdiger, and Edward ReijerseView Author Information Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany*E-mail: [email protected]Cite this: Chem. Rev. 2014, 114, 8, 4081–4148Publication Date (Web):March 21, 2014Publication History Received16 October 2013Published online21 March 2014Published inissue 23 April 2014https://pubs.acs.org/doi/10.1021/cr4005814https://doi.org/10.1021/cr4005814review-articleACS PublicationsCopyright © 2014 American Chemical SocietyRequest reuse permissionsArticle Views26977Altmetric-Citations1637LEARN ABOUT THESE METRICSArticle Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated. Share Add toView InAdd Full Text with ReferenceAdd Description ExportRISCitationCitation and abstractCitation and referencesMore Options Share onFacebookTwitterWechatLinked InRedditEmail Other access optionsGet e-Alertsclose SUBJECTS:Cluster chemistry,Ligands,Monomers,Peptides and proteins,Redox reactions Get e-Alerts

In previous publications, it was shown that an efficient local coupled cluster method with single- and double excitations can be based on the concept of pair natural orbitals (PNOs) [F. Neese, A. Hansen, and D. G. Liakos, J. Chem. Phys. 131, 064103 (2009)]. The resulting local pair natural orbital-coupled-cluster single double (LPNO-CCSD) method has since been proven to be highly reliable and efficient. For large molecules, the number of amplitudes to be determined is reduced by a factor of 10(5)-10(6) relative to a canonical CCSD calculation on the same system with the same basis set. In the original method, the PNOs were expanded in the set of canonical virtual orbitals and single excitations were not truncated. This led to a number of fifth order scaling steps that eventually rendered the method computationally expensive for large molecules (e.g., >100 atoms). In the present work, these limitations are overcome by a complete redesign of the LPNO-CCSD method. The new method is based on the combination of the concepts of PNOs and projected atomic orbitals (PAOs). Thus, each PNO is expanded in a set of PAOs that in turn belong to a given electron pair specific domain. In this way, it is possible to fully exploit locality while maintaining the extremely high compactness of the original LPNO-CCSD wavefunction. No terms are dropped from the CCSD equations and domains are chosen conservatively. The correlation energy loss due to the domains remains below <0.05%, which implies typically 15-20 but occasionally up to 30 atoms per domain on average. The new method has been given the acronym DLPNO-CCSD ("domain based LPNO-CCSD"). The method is nearly linear scaling with respect to system size. The original LPNO-CCSD method had three adjustable truncation thresholds that were chosen conservatively and do not need to be changed for actual applications. In the present treatment, no additional truncation parameters have been introduced. Any additional truncation is performed on the basis of the three original thresholds. There are no real-space cutoffs. Single excitations are truncated using singles-specific natural orbitals. Pairs are prescreened according to a multipole expansion of a pair correlation energy estimate based on local orbital specific virtual orbitals (LOSVs). Like its LPNO-CCSD predecessor, the method is completely of black box character and does not require any user adjustments. It is shown here that DLPNO-CCSD is as accurate as LPNO-CCSD while leading to computational savings exceeding one order of magnitude for larger systems. The largest calculations reported here featured >8800 basis functions and >450 atoms. In all larger test calculations done so far, the LPNO-CCSD step took less time than the preceding Hartree-Fock calculation, provided no approximations have been introduced in the latter. Thus, based on the present development reliable CCSD calculations on large molecules with unprecedented efficiency and accuracy are realized.

In this work, the extension of the previously developed domain based local pair-natural orbital (DLPNO) based singles- and doubles coupled cluster (DLPNO-CCSD) method to perturbatively include connected triple excitations is reported. The development is based on the concept of triples-natural orbitals that span the joint space of the three pair natural orbital (PNO) spaces of the three electron pairs that are involved in the calculation of a given triple-excitation contribution. The truncation error is very smooth and can be significantly reduced through extrapolation to the zero threshold. However, the extrapolation procedure does not improve relative energies. The overall computational effort of the method is asymptotically linear with the system size O(N). Actual linear scaling has been confirmed in test calculations on alkane chains. The accuracy of the DLPNO-CCSD(T) approximation relative to semicanonical CCSD(T0) is comparable to the previously developed DLPNO-CCSD method relative to canonical CCSD. Relative energies are predicted with an average error of approximately 0.5 kcal∕mol for a challenging test set of medium sized organic molecules. The triples correction typically adds 30%-50% to the overall computation time. Thus, very large systems can be treated on the basis of the current implementation. In addition to the linear C150H302 (452 atoms, >8800 basis functions) we demonstrate the first CCSD(T) level calculation on an entire protein, Crambin with 644 atoms, and more than 6400 basis functions.

The methodology presented within this work is a result of years of interactions between many junior and senior X-ray Photoelectron Spectroscopy (XPS) users operating within the CasaXPS spectral processing and interpretation program framework. In particular, discussions arising from a series of workshops have been a significant source for developing the overall XPS data processing concept and are the motivation for creating this work. These workshops organized by the Institut des Matériaux Jean Rouxel (IMN), Nantes gather both experienced and novice users of XPS for a week of discourse in conceptual experiment design and the resulting data processing. However, the framework constructed and utilized within these workshops encouraged the dissemination of knowledge beyond XPS data analysis and emphasized the importance of a multi-disciplinary collaborative approach to surface analysis problem-solving. The material presented here embodies data treatment originating from data made available to the first CNRS Thematic Workshop presented at Roscoff 2013. The methodology described here has evolved over the subsequent workshops in 2016 and 2019 and currently represents the philosophy used in CasaXPS spectral data processing paradigm.

The material basis of a sustainable society will depend on chemical products and processes that are designed following principles that make them conducive to life. Important inherent properties of molecules need to be considered from the earliest stage-the design stage-to address whether compounds and processes are depleting versus renewable, toxic versus benign, and persistent versus readily degradable. Products, feedstocks, and manufacturing processes will need to integrate the principles of green chemistry and green engineering under an expanded definition of performance that includes sustainability considerations. This transformation will require the best of the traditions of science and innovation coupled with new emerging systems thinking and systems design that begins at the molecular level and results in a positive impact on the global scale.

NiFe and CoFe (MFe) layered double hydroxides (LDHs) are among the most active electrocatalysts for the alkaline oxygen evolution reaction (OER). Herein, we combine electrochemical measurements, operando X-ray scattering and absorption spectroscopy, and density functional theory (DFT) calculations to elucidate the catalytically active phase, reaction center and the OER mechanism. We provide the first direct atomic-scale evidence that, under applied anodic potentials, MFe LDHs oxidize from as-prepared α-phases to activated γ-phases. The OER-active γ-phases are characterized by about 8% contraction of the lattice spacing and switching of the intercalated ions. DFT calculations reveal that the OER proceeds via a Mars van Krevelen mechanism. The flexible electronic structure of the surface Fe sites, and their synergy with nearest-neighbor M sites through formation of O-bridged Fe-M reaction centers, stabilize OER intermediates that are unfavorable on pure M-M centers and single Fe sites, fundamentally accounting for the high catalytic activity of MFe LDHs.

Domain based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations (DLPNO-CCSD(T)) is a highly efficient local correlation method. It is known to be accurate and robust and can be used in a black box fashion in order to obtain coupled cluster quality total energies for large molecules with several hundred atoms. While previous implementations showed near linear scaling up to a few hundred atoms, several nonlinear scaling steps limited the applicability of the method for very large systems. In this work, these limitations are overcome and a linear scaling DLPNO-CCSD(T) method for closed shell systems is reported. The new implementation is based on the concept of sparse maps that was introduced in Part I of this series [P. Pinski, C. Riplinger, E. F. Valeev, and F. Neese, J. Chem. Phys. 143, 034108 (2015)]. Using the sparse map infrastructure, all essential computational steps (integral transformation and storage, initial guess, pair natural orbital construction, amplitude iterations, triples correction) are achieved in a linear scaling fashion. In addition, a number of additional algorithmic improvements are reported that lead to significant speedups of the method. The new, linear-scaling DLPNO-CCSD(T) implementation typically is 7 times faster than the previous implementation and consumes 4 times less disk space for large three-dimensional systems. For linear systems, the performance gains and memory savings are substantially larger. Calculations with more than 20 000 basis functions and 1000 atoms are reported in this work. In all cases, the time required for the coupled cluster step is comparable to or lower than for the preceding Hartree-Fock calculation, even if this is carried out with the efficient resolution-of-the-identity and chain-of-spheres approximations. The new implementation even reduces the error in absolute correlation energies by about a factor of two, compared to the already accurate previous implementation.

Efficient reversible oxygen electrodes for both the oxygen reduction reaction (ORR) and the oxygen evolution reaction (OER) are vitally important for various energy conversion devices, such as regenerative fuel cells and metal-air batteries. However, realization of such electrodes is impeded by insufficient activity and instability of electrocatalysts for both water splitting and oxygen reduction. We report highly active bifunctional electrocatalysts for oxygen electrodes comprising core-shell Co@Co3O4 nanoparticles embedded in CNT-grafted N-doped carbon-polyhedra obtained by the pyrolysis of cobalt metal-organic framework (ZIF-67) in a reductive H2 atmosphere and subsequent controlled oxidative calcination. The catalysts afford 0.85 V reversible overvoltage in 0.1 m KOH, surpassing Pt/C, IrO2 , and RuO2 and thus ranking them among one of the best non-precious-metal electrocatalysts for reversible oxygen electrodes.

A heterogeneous catalyst is a functional material that continually creates active sites with its reactants under reaction conditions. These sites change the rates of chemical reactions of the reactants localized on them without changing the thermodynamic equilibrium between the materials.

In this communication, an improved perturbative triples correction (T) algorithm for domain based local pair-natural orbital singles and doubles coupled cluster (DLPNO-CCSD) theory is reported. In our previous implementation, the semi-canonical approximation was used and linear scaling was achieved for both the DLPNO-CCSD and (T) parts of the calculation. In this work, we refer to this previous method as DLPNO-CCSD(T0) to emphasize the semi-canonical approximation. It is well-established that the DLPNO-CCSD method can predict very accurate absolute and relative energies with respect to the parent canonical CCSD method. However, the (T0) approximation may introduce significant errors in absolute energies as the triples correction grows up in magnitude. In the majority of cases, the relative energies from (T0) are as accurate as the canonical (T) results of themselves. Unfortunately, in rare cases and in particular for small gap systems, the (T0) approximation breaks down and relative energies show large deviations from the parent canonical CCSD(T) results. To address this problem, an iterative (T) algorithm based on the previous DLPNO-CCSD(T0) algorithm has been implemented [abbreviated here as DLPNO-CCSD(T)]. Using triples natural orbitals to represent the virtual spaces for triples amplitudes, storage bottlenecks are avoided. Various carefully designed approximations ease the computational burden such that overall, the increase in the DLPNO-(T) calculation time over DLPNO-(T0) only amounts to a factor of about two (depending on the basis set). Benchmark calculations for the GMTKN30 database show that compared to DLPNO-CCSD(T0), the errors in absolute energies are greatly reduced and relative energies are moderately improved. The particularly problematic case of cumulene chains of increasing lengths is also successfully addressed by DLPNO-CCSD(T).

The domain based local pair natural orbital coupled cluster method with single-, double-, and perturbative triple excitations (DLPNO–CCSD(T)) is an efficient quantum chemical method that allows for coupled cluster calculations on molecules with hundreds of atoms. Because coupled-cluster theory is the method of choice if high-accuracy is needed, DLPNO–CCSD(T) is very promising for large-scale chemical application. However, the various approximations that have to be introduced in order to reach near linear scaling also introduce limited deviations from the canonical results. In the present work, we investigate how far the accuracy of the DLPNO–CCSD(T) method can be pushed for chemical applications. We also address the question at which additional computational cost improvements, relative to the previously established default scheme, come. To answer these questions, a series of benchmark sets covering a broad range of quantum chemical applications including reaction energies, hydrogen bonds, and other noncovalent interactions, conformer energies, and a prototype organometallic problem were selected. An accuracy of 1 kcal/mol or better can readily be obtained for all data sets using the default truncation scheme, which corresponds to the stated goal of the original implementation. Tightening of the three thresholds that control DLPNO leads to mean absolute errors and standard deviations from the canonical results of less than 0.25 kcal/mol (<1 kJ/mol). The price one has then to pay is an increased computational time by a factor close to 3. The applicability of the method is shown to be independent of the nature of the reaction. On the basis of the careful analysis of the results, three different sets of truncation thresholds (termed “LoosePNO”, “NormalPNO”, and “TightPNO”) have been chosen for “black box” use of DLPNO–CCSD(T). This will allow users of the method to optimally balance performance and accuracy.

A procedure was developed to automatically generate auxiliary basis sets (ABSs) for use with the resolution of the identity (RI) approximation, starting from a given orbital basis set (OBS). The goal is to provide an accurate and universal solution for cases where no optimized ABSs are available. In this context, “universal” is understood as the ability of the ABS to be used for Coulomb, exchange, and correlation energy fitting. The generation scheme (denoted AutoAux) works by spanning the product space of the OBS using an even-tempered expansion for each atom in the system. The performance of AutoAux in conjunction with different OBSs [def2-SVP, def2-TZVP, def2-QZVPP, and cc-pwCV n Z ( n = D, T, Q, 5)] has been evaluated for elements from H to Rn and compared to existing predefined ABSs. Due to the requirements of simplicity and universality, the generated ABSs are larger than the optimized ones but lead to similar errors in MP2 total energies (on the order of 10 –5 to 10 –4 E h /atom).

In industrially relevant Cu/ZnO/Al2 O3 catalysts for methanol synthesis, the strong metal support interaction between Cu and ZnO is known to play a key role. Here we report a detailed chemical transmission electron microscopy study on the nanostructural consequences of the strong metal support interaction in an activated high-performance catalyst. For the first time, clear evidence for the formation of metastable "graphite-like" ZnO layers during reductive activation is provided. The description of this metastable layer might contribute to the understanding of synergistic effects between the components of the Cu/ZnO/Al2 O3 catalysts.

Producing propene by the oxidative dehydrogenation of propane (ODH) has become an attractive and feasible route for bridging the propene production-demand gap, either as a complementary route of the existing oil-based processes or as a new alternative from propane separated from natural gas. The industrial application of propane ODH has not succeeded so far due to low propene yields. Therefore, propane ODH has been extensively investigated in recent decades using different catalysts and reaction conditions. Although several important aspects have been discussed in previous reviews (e.g., supported vanadium oxide catalysts, bulk catalysts, productivity toward propene, etc.), other relevant aspects have not been addressed (e.g., support effects, loading effects, vanadia precursor or catalyst synthesis methods, surface impurities, structure–reactivity relationships, etc.). In this review, we endeavor to cover the majority of the publications with an emphasis on the following: (1) catalyst synthesis: to focus on the influence of synthesis methods on the final vanadium oxide surface species, (2) catalyst characterization: to identify the molecular structures of the supported vanadium oxide species as well as the oxide support surface physical and chemical characteristics, (3) kinetics: to understand how reaction rates depend on variables such as concentration of gas-phase reactants and temperature, (4) structure–activity relationship: to examine the influence of the concentration as well as molecular structures of the surface vanadium oxide species on the reaction kinetics, and (5) reaction mechanism: to use the structure–activity relationships as well as kinetic studies plus theoretical calculations to corroborate and/or propose reaction pathways that account for the overall ODP reaction mechanism.

Atrazine, propazine, and terbuthylazine are chlorotriazine herbicides that have been frequently used in agriculture and thus are potential drinking water contaminants. Hydroxyl radicals produced by advanced oxidation processes can degrade these persistent compounds. These herbicides are also very reactive with sulfate radicals (2.2-3.5 × 10(9) M(-1) s(-1)). However, the dealkylated products of chlorotriazine pesticides are less reactive toward sulfate radicals (e.g., desethyl-desisopropyl-atrazine (DEDIA; 1.5 × 10(8) M(-1) s(-1))). The high reactivity of the herbicides is largely due to the ethyl or isopropyl group. For example, desisopropyl-atrazine (DIA) reacts quickly (k = 2 × 10(9) M(-1) s(-1)), whereas desethyl-atrazine (DEA) reacts more slowly (k = 9.6 × 10(8) M(-1) s(-1)). The tert-butyl group does not have a strong effect on reaction rate, as shown by the similar second order reaction rates between desethyl-terbuthylazine (DET; k = 3.6 × 10(8) M(-1) s(-1)) and DEDIA. Sulfate radicals degrade a significant proportion of atrazine (63%) via dealkylation, in which deethylation significantly dominates over deisopropylation (10:1). Sulfate and hydroxyl radicals react at an equally fast rate with atrazine (k (hydroxyl radical + atrazine) = 3 × 10(9) M(-1) s(-1)). However, sulfate and hydroxyl radicals differ considerably in their reaction rates with humic acids (k (sulfate radical + humic acids) = 6.8 × 10(3) L mgC(-1) s(-1) (mgC = mg carbon); k (hydroxyl radical + humic acids) = 1.4 × 10(4) L mgC(-1) s(-1)). Thus, in the presence of humic acids, atrazine is degraded more efficiently by sulfate radicals than by hydroxyl radicals.

A series of two-coordinate complexes of iron(II) were prepared and studied for single-molecule magnet behavior. Five of the compounds, Fe[N(SiMe3)(Dipp)]2 (1), Fe[C(SiMe3)3]2 (2), Fe[N(H)Ar′]2 (3), Fe[N(H)Ar*]2 (4), and Fe(OAr′)2 (5) feature a linear geometry at the FeII center, while the sixth compound, Fe[N(H)Ar#]2 (6), is bent with an N–Fe–N angle of 140.9(2)° (Dipp = C6H3-2,6-Pri2; Ar′ = C6H3-2,6-(C6H3-2,6-Pri2)2; Ar* = C6H3-2,6-(C6H2-2,4,6-Pri2)2; Ar# = C6H3-2,6-(C6H2-2,4,6-Me3)2). Ac magnetic susceptibility data for all compounds revealed slow magnetic relaxation under an applied dc field, with the magnetic relaxation times following a general trend of 1 > 2 > 3 > 4 > 5 ≫ 6. Arrhenius plots created for the linear complexes were fit by employing a sum of tunneling, direct, Raman, and Orbach relaxation processes, resulting in spin reversal barriers of Ueff = 181, 146, 109, 104, and 43 cm−1 for 1–5, respectively. CASSCF/NEVPT2 calculations on the crystal structures were performed to explore the influence of deviations from rigorous D∞h geometry on the d-orbital splittings and the electronic state energies. Asymmetry in the ligand fields quenches the orbital angular momentum of 1–6, but ultimately spin–orbit coupling is strong enough to compensate and regenerate the orbital moment. The lack of simple Arrhenius behavior in 1–5 can be attributed to a combination of the asymmetric ligand field and the influence of vibronic coupling, with the latter possibility being suggested by thermal ellipsoid models to the diffraction data.

Abstract Methanol, an important chemical, fuel additive, and precursor for clean fuels, is produced by hydrogenation of carbon oxides over Cu‐based catalysts. Despite the technological maturity of this process, the understanding of this apparently simple reaction is still incomplete with regard to the reaction mechanism and the active sites. Regarding the latter, recent progress has shown that stepped and ZnO x ‐decorated Cu surfaces are crucial for the performance of industrial catalysts. Herein, we integrate this insight with additional experiments into a full microkinetic description of methanol synthesis. In particular, we show how the presence or absence of the Zn promoter dramatically changes not only the activity, but unexpectedly the reaction mechanism itself. The Janus‐faced character of Cu with two different sites for methanol synthesis, Zn‐promoted and unpromoted, resolves the long‐standing controversy regarding the Cu/Zn synergy and adds methanol synthesis to the few major industrial catalytic processes that are described on an atomic level.

Hydrogen plays a key role in many industrial applications and is currently seen as one of the most promising energy vectors. Many efforts are being made to produce hydrogen with zero CO2 footprint via water electrolysis powered by renewable energies. Nevertheless, the use of fossil fuels is essential in the short term. The conventional coal gasification and steam methane reforming processes for hydrogen production are undesirable due to the huge CO2 emissions. A cleaner technology based on natural gas that has received special attention in recent years is methane pyrolysis. The thermal decomposition of methane gives rise to hydrogen and solid carbon, and thus, the release of greenhouse gases is prevented. Therefore, methane pyrolysis is a CO2-free technology that can serve as a bridge from fossil fuels to renewable energies.