Two Allogons of an O2 -activating Bis(disiloxido)ferrate(II) Accessible Selectively just by Variation of the Crystallization Temperature.

Deprotonation of O(iPr2 SiOH)2 (iPr LH2 ) with LiOtBu followed by reaction with FeCl2 in THF led to the complex [iPr L2 Fe][Li(THF)2 ]2 , 2, which represents a structural and spectroscopic model of the α-Fe sites of Fe/ZSM-5. Reaction with O2 in THF solution proceeds rather fast and is complete within 200 ms; an intermediate O2 adduct could not be identified by stopped-flow methods. Cooling blue solutions of 2 to -80 °C led to the growth of blue crystals of 2⋅THF, the analysis of which by XRD revealed a FeO4 core that is somewhat distorted from planarity towards a tetrahedral structure. By contrast, cooling such solutions to -30 °C led to pink crystals of an allogon featuring a perfectly square planar FeO4 entity. Hence, 2 represents a unique case where two different structural isomers (allogons) can be crystallized from the same solvent selectively, controlled by the temperature. DFT calculations were performed to understand this finding.

Reaction barriers on non-conducting surfaces beyond periodic local MP2: Diffusion of hydrogen on α-Al2O3(0001) as a test case.

The quest for "chemical accuracy" is becoming more and more demanded in the field of structure and kinetics of molecules at solid surfaces. In this paper, as an example, we focus on the barrier for hydrogen diffusion on a α-Al2O3(0001) surface, aiming for a couple cluster singles, doubles, and perturbative triples [CCSD(T)]-level benchmark. We employ the density functional theory (DFT) optimized minimum and transition state structures reported by Heiden, Usvyat, and Saalfrank [J. Phys. Chem. C 123, 6675 (2019)]. The barrier is first evaluated at the periodic Hartree-Fock and local Møller-Plesset second-order perturbation (MP2) level of theory. The possible sources of errors are then analyzed, which includes basis set incompleteness error, frozen core, density fitting, local approximation errors, as well as the MP2 method error. Using periodic and embedded fragment models, corrections to these errors are evaluated. In particular, two corrections are found to be non-negligible (both from the chemical accuracy perspective and at the scale of the barrier value of 0.72 eV): the correction to the frozen core-approximation of 0.06 eV and the CCSD(T) correction of 0.07 eV. Our correlated wave function results are compared to barriers obtained from DFT. Among the tested DFT functionals, the best performing for this barrier is B3LYP-D3.

Full configuration interaction quantum Monte Carlo treatment of fragments embedded in a periodic mean field.

We present an embedded fragment approach for high-level quantum chemical calculations on local features in periodic systems. The fragment is defined as a set of localized orbitals (occupied and virtual) corresponding to a converged periodic Hartree-Fock solution. These orbitals serve as the basis for the in-fragment post-Hartree-Fock treatment. The embedding field for the fragment, consisting of the Coulomb and exchange potential from the rest of the crystal, is included in the fragment's one-electron Hamiltonian. As an application of the embedded fragment approach, we investigate the performance of full configuration interaction quantum Monte Carlo (FCIQMC) with the adaptive shift. As the orbital choice, we use the natural orbitals from the distinguishable cluster method with singles and doubles. FCIQMC is a stochastic approximation to the full CI method and can be routinely applied to much larger active spaces than the latter. This makes this method especially attractive in the context of open shell defects in crystals, where fragments of adequate size can be rather large. As a test case, we consider dissociation of a fluorine atom from a fluorographane surface. This process poses a challenge for high-level electronic structure models as both the static and dynamic correlations are essential here. Furthermore, the active space for an adequate fragment (32 electrons in 173 orbitals) is already quite large even for FCIQMC. Despite this, FCIQMC delivers accurate dissociation and total energies.

Insights into Reaction Kinetics in Confined Space: Real Time Observation of Water Formation under a Silica Cover.

We offer a comprehensive approach to determine how physical confinement can affect the water formation reaction. By using free-standing crystalline SiO2 bilayer supported on Ru(0001) as a model system, we studied the water formation reaction under confinement in situ and in real time. Low-energy electron microscopy reveals that the reaction proceeds via the formation of reaction fronts propagating across the Ru(0001) surface. The Arrhenius analyses of the front velocity yield apparent activation energies (Eaapp) of 0.32 eV for the confined and 0.59 eV for the nonconfined reaction. DFT simulations indicate that the rate-determining step remains unchanged upon confinement, therefore ruling out the widely accepted transition state effect. Additionally, H2O accumulation cannot explain the change in Eaapp for the confined cases studied because its concentration remains low. Instead, numerical simulations of the proposed kinetic model suggest that the H2 adsorption process plays a decisive role in reproducing the Arrhenius plots.

Growth and Atomic-Scale Characterization of Ultrathin Silica and Germania Films: The Crucial Role of the Metal Support.

The present review reports on the preparation and atomic-scale characterization of the thinnest possible films of the glass-forming materials silica and germania. To this end state-of-the-art surface science techniques, in particular scanning probe microscopy, and density functional theory calculations have been employed. The investigated films range from monolayer to bilayer coverage where both, the crystalline and the amorphous films, contain characteristic XO4 (X=Si,Ge) building blocks. A side-by-side comparison of silica and germania monolayer, zigzag phase and bilayer films supported on Mo(112), Ru(0001), Pt(111), and Au(111) leads to a more general comprehension of the network structure of glass former materials. This allows us to understand the crucial role of the metal support for the pathway from crystalline to amorphous ultrathin film growth.

The color center singlet state of oxygen vacancies in TiO2.

Oxygen vacancies are ubiquitous in TiO2 and play key roles in catalysis and magnetism applications. Despite being extensively investigated, the electronic structure of oxygen vacancies in TiO2 remains controversial both experimentally and theoretically. Here, we report a study of a neutral oxygen vacancy in TiO2 using state-of-the-art quantum chemical electronic structure methods. We find that the ground state is a color center singlet state in both the rutile and the anatase phases of TiO2. Specifically, embedded coupled cluster with singles, doubles, and perturbative triples calculations find, for an oxygen vacancy in rutile, that the lowest triplet state energy is 0.6 eV above the singlet state, and in anatase, the triplet state energy is higher by 1.4 eV. Our study provides fresh insights into the electronic structure of the oxygen vacancy in TiO2, clarifying earlier controversies and potentially inspiring future studies of defects with correlated wave function theories.

Evolutionary Algorithm-Based Crystal Structure Prediction for Copper(I) Fluoride.

Despite numerous experimental studies since 1824, the binary copper(I) fluoride remains unknown. A crystal structure prediction has been carried out for CuF using the USPEX evolutionary algorithm and a dispersion-corrected hybrid density functional method. In total about 5000 hypothetical structures were investigated. The energetics of the predicted structures were also counter-checked with local second-order Møller-Plesset perturbation theory. Herein 39 new hypothetical copper(I) fluoride structures are reported that are lower in energy compared to the previously predicted cinnabar-type structure. Cuprophilic Cu-Cu interactions are present in all the low-energy structures, leading to ordered Cu substructures such as helical or zig-zag-type Cu-Cu motifs. The lowest-energy structure adopts a trigonal crystal structure with space group P31 21. From an electronic point of view, the predicted CuF modification is a semiconductor with an indirect band gap of 2.3 eV.

A Two-Dimensional 'Zigzag' Silica Polymorph on a Metal Support.

We present a new polymorph of the two-dimensional (2D) silica film with a characteristic 'zigzag' line structure and a rectangular unit cell which forms on a Ru(0001) metal substrate. This new silica polymorph may allow for important insights into growth modes and transformations of 2D silica films as a model system for the study of glass transitions. Based on scanning tunneling microscopy, low energy electron diffraction, infrared reflection absorption spectroscopy, and X-ray photoelectron spectroscopy measurements on the one hand, and density functional theory calculations on the other, a structural model for the 'zigzag' polymorph is proposed. In comparison to established monolayer and bilayer silica, this 'zigzag' structure system has intermediate characteristics in terms of coupling to the substrate and stoichiometry. The silica 'zigzag' phase is transformed upon reoxidation at higher annealing temperature into a SiO2 silica bilayer film which is chemically decoupled from the substrate.

Dispersion interactions in silicon allotropes.

van der Waals interactions are known to play a key role in the formation of weakly bound solids, such as molecular or layered crystals. Here we show that the correct quantum-chemical description of van der Waals dispersion is also essential for a correct description of the relative stability between purely covalently-bound solids like silicon allotropes. To this end, we apply periodic local MP2 and DFT with Grimme's empirical -D3 correction to 11 experimentally determined or yet hypothetical crystalline silicon structures, including the most recently discovered silicon allotropes. Both methods provide similar energy ordering of the polymorphs, which, at the same time, noticeably deviate from the order predicted by standard DFT without an appropriate description of the van der Waals dispersion.

A comparison between quantum chemistry and quantum Monte Carlo techniques for the adsorption of water on the (001) LiH surface.

We present a comprehensive benchmark study of the adsorption energy of a single water molecule on the (001) LiH surface using periodic coupled cluster and quantum Monte Carlo theories. We benchmark and compare different implementations of quantum chemical wave function based theories in order to verify the reliability of the predicted adsorption energies and the employed approximations. Furthermore we compare the predicted adsorption energies to those obtained employing widely used van der Waals density-functionals. Our findings show that quantum chemical approaches are becoming a robust and reliable tool for condensed phase electronic structure calculations, providing an additional tool that can also help in potentially improving currently available van der Waals density-functionals.

A hierarchy of local coupled cluster singles and doubles response methods for ionization potentials.

We present a hierarchy of local coupled cluster (CC) linear response (LR) methods to calculate ionization potentials (IPs), i.e., excited states with one electron annihilated relative to a ground state reference. The time-dependent perturbation operator V(t), as well as the operators related to the first-order (with respect to V(t)) amplitudes and multipliers, thus are not number conserving and have half-integer particle rank m. Apart from calculating IPs of neutral molecules, the method offers also the possibility to study ground and excited states of neutral radicals as ionized states of closed-shell anions. It turns out that for comparable accuracy IPs require a higher-order treatment than excitation energies; an IP-CC LR method corresponding to CC2 LR or the algebraic diagrammatic construction scheme through second order performs rather poorly. We therefore systematically extended the order with respect to the fluctuation potential of the IP-CC2 LR Jacobian up to IP-CCSD LR, keeping the excitation space of the first-order (with respect to V(t)) cluster operator restricted to the m=½âŠ•3/2 subspace and the accuracy of the zero-order (ground-state) amplitudes at the level of CC2 or MP2. For the more expensive diagrams beyond the IP-CC2 LR Jacobian, we employ local approximations. The implemented methods are capable of treating large molecular systems with hundred atoms or more.

High precision quantum-chemical treatment of adsorption: Benchmarking physisorption of molecular hydrogen on graphane.

A multilevel hierarchical ab initio protocol for calculating adsorption on non-conducting surfaces is presented. It employs fully periodic treatment, which reaches local Møller-Plesset perturbation theory of second order (MP2) with correction for the basis set incompleteness via the local F12 technique. Post-MP2 corrections are calculated using finite clusters. That includes the coupled cluster treatment in the local and canonical frameworks (up to perturbative quadruples) and correlated core (with MP2). Using this protocol, the potential surface of hydrogen molecules adsorbed on graphane was computed. According to the calculations, hydrogen molecules are adsorbed on graphane in a perpendicular to the surface orientation with the minimum of the potential surface of around -3.6 kJ/mol located at the distance of 3.85 Å between the bond center of the hydrogen molecule and the mid-plane of graphane. The adsorption sites along the path from the downward-pointing carbon to the ring center of the graphane are energetically virtually equally preferable, which can enable nearly free translations of hydrogen molecules along these paths. Consequently, the hydrogen molecules on graphane most likely form a non-commensurate monolayer. The analysis of the remaining errors reveals a very high accuracy of the computed potential surface with an error bar of a few tenths of a kJ/mol. The obtained results are a high-precision benchmark for further theoretical and experimental studies of hydrogen molecules interacting with graphane.

Periodic local MP2 method employing orbital specific virtuals.

We introduce orbital specific virtuals (OSVs) to represent the truncated pair-specific virtual space in periodic local Møller-Plesset perturbation theory of second order (LMP2). The OSVs are constructed by diagonalization of the LMP2 amplitude matrices which correspond to diagonal Wannier-function (WF) pairs. Only a subset of these OSVs is adopted for the subsequent OSV-LMP2 calculation, namely, those with largest contribution to the diagonal pair correlation energy and with the accumulated value of these contributions reaching a certain accuracy. The virtual space for a general (non diagonal) pair is spanned by the union of the two OSV sets related to the individual WFs of the pair. In the periodic LMP2 method, the diagonal LMP2 amplitude matrices needed for the construction of the OSVs are calculated in the basis of projected atomic orbitals (PAOs), employing very large PAO domains. It turns out that the OSVs are excellent to describe short range correlation, yet less appropriate for long range van der Waals correlation. In order to compensate for this bias towards short range correlation, we augment the virtual space spanned by the OSVs by the most diffuse PAOs of the corresponding minimal PAO domain. The Fock and overlap matrices in OSV basis are constructed in the reciprocal space. The 4-index electron repulsion integrals are calculated by local density fitting and, for distant pairs, via multipole approximation. New procedures for determining the fit-domains and the distant-pair lists, leading to higher efficiency in the 4-index integral evaluation, have been implemented. Generally, and in contrast to our previous PAO based periodic LMP2 method, the OSV-LMP2 method does not require anymore great care in the specification of the individual domains (to get a balanced description when calculating energy differences) and is in that sense a black box procedure. Discontinuities in potential energy surfaces, which may occur for PAO-based calculations if one is not careful, virtually disappear for OSV-LMP2. Moreover, due to much increased compactness of the pair-specific virtual spaces, the OSV-LMP2 calculations are faster and require much less memory than PAO-LMP2 calculations, despite the noticeable overhead of the initial OSV construction procedure.

Communication: Improved pair approximations in local coupled-cluster methods.

In local coupled cluster treatments the electron pairs can be classified according to the magnitude of their energy contributions or distances into strong, close, weak, and distant pairs. Different approximations are introduced for the latter three classes. In this communication, an improved simplified treatment of close and weak pairs is proposed, which is based on long-range cancellations of individually slowly decaying contributions in the amplitude equations. Benchmark calculations for correlation, reaction, and activation energies demonstrate that these approximations work extremely well, while pair approximations based on local second-order Møller-Plesset theory can lead to errors that are 1-2 orders of magnitude larger.

Range-separated double-hybrid density-functional theory applied to periodic systems.

Quantum chemistry methods exploiting density-functional approximations for short-range electron-electron interactions and second-order Møller-Plesset (MP2) perturbation theory for long-range electron-electron interactions have been implemented for periodic systems using Gaussian-type basis functions and the local correlation framework. The performance of these range-separated double hybrids has been benchmarked on a significant set of systems including rare-gas, molecular, ionic, and covalent crystals. The use of spin-component-scaled MP2 for the long-range part has been tested as well. The results show that the value of µ = 0.5 bohr(-1) for the range-separation parameter usually used for molecular systems is also a reasonable choice for solids. Overall, these range-separated double hybrids provide a good accuracy for binding energies using basis sets of moderate sizes such as cc-pVDZ and aug-cc-pVDZ.

Diffraction of helium on MgO(100) surface calculated from first-principles.

In this work we simulate the diffraction peak intensities of He beams scattered on the MgO(100) surface from first principles. It turns out that diffraction peak intensities are extremely sensitive to the quality of the potential describing the He-MgO surface interaction. Achieving the required accuracy in first principles calculations is very challenging indeed. The present work describes a first principles protocol able to achieve very high accuracy for reasonable computational cost. This method is based on periodic local second-order Møller-Plesset perturbation theory where systematic corrections for basis set truncation and for high-order electronic correlation are introduced using coupled cluster calculations on finite model systems mimicking the target system. For the He-MgO system the requirements with respect to the level of theory are very high; it turns out that contributions from connected quadruple excitations are non-negligible. Here we demonstrate that using this protocol, it is possible to reach the accuracy in the He-MgO potential that is required to predict the observed He diffraction peak intensities.

Efficient and accurate treatment of weak pairs in local CCSD(T) calculations. II. Beyond the ring approximation.

In order to arrive at linear scaling of the computational cost with molecular size, local coupled cluster methods discriminate pairs of local molecular orbitals according to the spatial separation R of the latter. Only strong pairs are treated at the full coupled cluster level, whereas for weak pairs a lower level of theory (usually Møller-Plesset perturbation theory of second order, MP2) is used. Yet an MP2 treatment of weak pairs is inadequate in certain situations (for example, for describing π-stacking), which calls for an improved but still inexpensive method for dealing with the weak pairs. In a previous contribution, we proposed as a substituent for MP2 the LrCCD3 method, which is based on ring coupled cluster doubles (ring-CCD) and includes all third-order diagrams with energy contributions decaying not quicker than R(-6). In the present work, we explore a still more accurate method, which is based on the same principles. It turned out to be essential to abandon the restriction to ring-CCD, i.e., to include further CCD diagrams beyond the ring approximation. The occurring intermediates turn out to be formally very similar to LMP2 density matrices, such that an efficient evaluation of these non-ring CCD diagrams is possible. Furthermore, a computationally cheap a posteriori estimate for the fourth-order singles contribution to the weak pair energy, which also exhibits a decay behavior of R(-6), is introduced. The resulting method, denoted as LCCD[S]-R(-6), indeed provides a substantial improvement in accuracy over the previous LrCCD3 method at a relatively modest additional computational cost.

Hybrid RPA:DFT Approach for Adsorption on Transition Metal Surfaces: Methane and Ethane on Platinum (111).

The hybrid QM:QM approach is extended to adsorption on transition metal surfaces. The random phase approximation (RPA) as the high-level method is applied to cluster models and, using the subtractive scheme, embedded in periodic models which are treated with density functional theory (DFT) that is the low-level method. The PBE functional, both without dispersion and augmented with the many-body dispersion (MBD), is employed. Adsorption of methane and ethane on the Pt(111) surface is studied. For methane in a 2 × 2 surface cell, the hybrid RPA:PBE and RPA:PBE+MBD results, -14.3 and -16.0 kJ mol-1, respectively, are in close agreement with the periodic RPA value of -13.8 kJ mol-1 at significantly reduced computational cost (factor of ∼50). For methane and ethane, the RPA:PBE results (-14.3 and -17.8 kJ mol-1, respectively) indicate underbinding relative to energies derived from experimental desorption barriers for relevant loadings (-15.6 ± 1.6 and -27.2 ± 2.9 kJ mol-1, respectively), whereas the hybrid RPA:PBE+MBD results (-16.0 and -24.9 kJ mol-1, respectively) agree with the experiment well within experimental uncertainty limits (deviation of -0.4 ± 1.5 and +2.3 ± 2.9 kJ mol-1, respectively). Finding a cluster that adequately and robustly represents the adsorbate at the bulk surface is important for the success of the RPA-based QM:QM scheme for metals.

Linear-scaling explicitly correlated treatment of solids: periodic local MP2-F12 method.

Theory and implementation of the periodic local MP2-F12 method in the 3*A fixed-amplitude ansatz is presented. The method is formulated in the direct space, employing local representation for the occupied, virtual, and auxiliary orbitals in the form of Wannier functions (WFs), projected atomic orbitals (PAOs), and atom-centered Gaussian-type orbitals, respectively. Local approximations are introduced, restricting the list of the explicitly correlated pairs, as well as occupied, virtual, and auxiliary spaces in the strong orthogonality projector to the pair-specific domains on the basis of spatial proximity of respective orbitals. The 4-index two-electron integrals appearing in the formalism are approximated via the direct-space density fitting technique. In this procedure, the fitting orbital spaces are also restricted to local fit-domains surrounding the fitted densities. The formulation of the method and its implementation exploits the translational symmetry and the site-group symmetries of the WFs. Test calculations are performed on LiH crystal. The results show that the periodic LMP2-F12 method substantially accelerates basis set convergence of the total correlation energy, and even more so the correlation energy differences. The resulting energies are quite insensitive to the resolution-of-the-identity domain sizes and the quality of the auxiliary basis sets. The convergence with the orbital domain size is somewhat slower, but still acceptable. Moreover, inclusion of slightly more diffuse functions, than those usually used in the periodic calculations, improves the convergence of the LMP2-F12 correlation energy with respect to both the size of the PAO-domains and the quality of the orbital basis set. At the same time, the essentially diffuse atomic orbitals from standard molecular basis sets, commonly utilized in molecular MP2-F12 calculations, but problematic in the periodic context, are not necessary for LMP2-F12 treatment of crystals.

Efficient and accurate treatment of weak pairs in local CCSD(T) calculations.

Local coupled cluster theory is based on (i) a restriction of the list of pairs (or triples) of occupied molecular orbitals, and (ii) a truncation of the virtual space to orbital pair (or triple) specific subspaces. The latter is motivated by an exponential decay of the contributions to the pair energy with respect to the distance between related local occupied and virtual orbitals; the former only by a polynomial R(-6) decay with respect to the distance R between the two occupied orbitals of the pair. Consequently, the restriction of the pair list is more critical, and contributions of pairs should not be neglected unless the corresponding interorbital distance is really large. In local coupled cluster theory pairs are usually discriminated on the basis of the interorbital distance, or the size of the 2nd order Møller-Plesset perturbation theory (MP2) estimate of the pair energy. Only strong pairs are treated at the full coupled cluster level, while weak pairs are treated just at the level of MP2. Yet MP2 might be problematic in certain cases, for example, π-stacking is badly described by MP2, etc. We propose to substitute the MP2 treatment of weak pairs by an approach based on ring-CCD by including third-order diagrams with R(-6) decay behavior. Such an approach is clearly superior; it provides higher accuracy, while the computational cost is not significantly higher than that of a MP2 based treatment of weak pairs.