Self-interaction corrected SCAN functional for molecules and solids in the numeric atom-center orbital framework

Semilocal density-functional approximations (DFAs), including the state-of-the-art SCAN functional, are plagued by the self-interaction error (SIE). While this error is explicitly defined only for one-electron systems, it has inspired the self-interaction correction method proposed by Perdew and Zunger (PZ-SIC), which has shown promise in mitigating the many-electron SIE. However, the PZ-SIC method is known for its significant numerical instability. In this study, we introduce a novel constraint that facilitates self-consistent localization of the SIC orbitals in the spirit of Edmiston-Ruedenberg orbitals [Rev. Mod. Phys. 35, 457 (1963)]. Our practical implementation within the all-electron numeric atom-centered orbitals code FHI-aims guarantees efficient and stable convergence of the self-consistent PZ-SIC equations for both molecules and solids. We further demonstrate that our PZ-SIC approach effectively mitigates the SIE in the meta-generalized gradient approximation SCAN functional, significantly improving the accuracy for ionization potentials, charge-transfer energies, and bandgaps for a diverse selection of molecules and solids. However, our PZ-SIC method does have its limitations. It cannot improve the already accurate SCAN results for properties such as cohesive energies, lattice constants, and bulk modulus in our test sets. This highlights the need for new-generation DFAs with more comprehensive applicability.

Developing correlation-consistent numeric atom-centered orbital basis sets for krypton: Applications in RPA-based correlated calculations

Localized atomic orbitals are the preferred basis set choice for large-scale explicit correlated calculations, and high-quality hierarchical correlation-consistent basis sets are a prerequisite for correlated methods to deliver numerically reliable results. At present, numeric atom-centered orbital (NAO) basis sets with valence correlation consistency (VCC), designated as NAO-VCC-nZ, are only available for light elements from hydrogen (H) to argon (Ar) [Zhang et al., New J. Phys. 15, 123033 (2013)]. In this work, we extend this series by developing NAO-VCC-nZ basis sets for krypton (Kr), a prototypical element in the fourth row of the periodic table. We demonstrate that NAO-VCC-nZ basis sets facilitate the convergence of electronic total-energy calculations using the Random Phase Approximation (RPA), which can be used together with a two-point extrapolation scheme to approach the complete basis set limit. Notably, the Basis Set Superposition Error (BSSE) associated with the newly generated NAO basis sets is minimal, making them suitable for applications where BSSE correction is either cumbersome or impractical to do. After confirming the reliability of NAO basis sets for Kr, we proceed to calculate the Helmholtz free energy for Kr crystal at the theoretical level of RPA plus renormalized single excitation correction. From this, we derive the pressure-volume (P-V) diagram, which shows excellent agreement with the latest experimental data. Our work demonstrates the capability of correlation-consistent NAO basis sets for heavy elements, paving the way toward numerically reliable correlated calculations for bulk materials.

Utilization of Fe-Ethylenediamine-N,N'-Disuccinic Acid Complex for Electrochemical Co-Catalytic Activation of Peroxymonosulfate under Neutral Initial pH Conditions

The ethylenediamine-N,N'-disuccinic acid (EDDS) was utilized to form Fe-EDDS complex to activate peroxymonosulfate (PMS) in the electrochemical (EC) co-catalytic system for effective oxidation of naphthenic acids (NAs) under neutral pH conditions. 1-adamantanecarboxylic acid (ACA) was used as a model compound to represent NAs, which are persistent pollutants that are abundantly present in oil and gas field wastewater. The ACA degradation rate was significantly enhanced in the EC/PMS/Fe(III)-EDDS system (96.6%) compared to that of the EC/PMS/Fe(III) system (65.4%). The addition of EDDS led to the formation of a stable complex of Fe-EDDS under neutral pH conditions, which effectively promoted the redox cycle of Fe(III)-EDDS/Fe(II)-EDDS to activate PMS to generate oxidative species for ACA degradation. The results of quenching and chemical probe experiments, as well as electron paramagnetic resonance (EPR) analysis, identified significant contributions of •OH, 1O2, and SO4•- in the removal of ACA. The ACA degradation pathways were revealed based on the results of high resolution mass spectrometry analysis and calculation of the Fukui index. The presence of anions, such as NO3-, Cl-, and HCO3-, as well as humic acids, induced nonsignificant influence on the ACA degradation, indicating the robustness of the current system for applications in authentic scenarios. Overall results indicated the EC/PMS/Fe(III)-EDDS system is a promising strategy for the practical treatment of NAs in oil and gas field wastewater.

Hydrogen sulfide functions as a micro-modulator bound at the copper active site of Cu/Zn-SOD to regulate the catalytic activity of the enzyme

The present study examines whether there is a mechanism beyond the current concept of post-translational modifications to regulate the function of a protein. A small gas molecule, hydrogen sulfide (H₂S), was found to bind at active-site copper of Cu/Zn-SOD using a series of methods including radiolabeled binding assay, X-ray absorption near-edge structure (XANES), and crystallography. Such an H₂S binding enhanced the electrostatic forces to guide the negatively charged substrate superoxide radicals to the catalytic copper ion, changed the geometry and energy of the frontier molecular orbitals of the active site, and subsequently facilitated the transfer of an electron from the superoxide radical to the catalytic copper ion and the breakage of the copper-His61 bridge. The physiological relevance of such an H₂S effect was also examined in both in vitro and in vivo models where the cardioprotective effects of H₂S were dependent on Cu/Zn-SOD.

Pulsed versus direct current electrochemical co-catalytic peroxymonosulfate-based system: Elevated degradation and energy efficiency with enhanced oxidation mechanisms

In this work, the pulsed electrochemical (PE) system was investigated to activate peroxymonosulfate (PMS) with the addition of Fe(III) to achieve efficient degradation of sulfamethoxazole (SMX) with reduced energy consumption, in comparison with the direct current (DC) electrochemical system. The operational conditions of PE/PMS/Fe(III) system were optimized as 4 kHz pulse frequency, 50% duty cycle, and pH 3, at which 67.6% reduction of energy consumption and enhanced degradation performance were achieved compared to the DC/PMS/Fe(III) system. Results of electron paramagnetic resonance spectroscopy analysis and quenching and chemical probe experiment revealed the presence of ^•OH, SO₄^•-, and ¹O₂ in the system, with ^•OH being the dominant role. The concentrations of these active species were averagely 15 ± 1% higher in the PE/PMS/Fe(III) system than those of the DC/PMS/Fe(III) system. Identification of SMX byproducts was achieved based on high resolution mass spectrometry analysis to predict the degradation pathways. The SMX byproducts could eventually be eliminated by the PE/PMS/Fe(III) system with extended treatment time. Overall, the PE/PMS/Fe(III) system was demonstrated with high energy and degradation performance, and is appear to be an robust strategy for practical treatment of wastewater.

Efficient removal of recalcitrant naphthenic acids with electro-cocatalytic activation of peroxymonosulfate by Fe(III)-nitrilotriacetic acid complex under neutral initial pH condition

This work investigated the activation of peroxymonosulfate by electrochemical (EC) system assisted with Fe(III)-nitrilotriacetic acid (NTA) complex for degradation of persistent naphthenic acids (NAs) under neutral initial pH conditions. As NAs are a complicated mixture, 1-adamantanecarboxylic acid (ACA) was selected as the model NA compound for degradation experiment. The addition of NTA is to chelate with Fe(III), gaining stability under neutral pH condition to facilitate the circulation of Fe(II)/Fe(III) by the electrochemical process to activate PMS. The EC/Fe(III)-NTA/PMS system was explored with applicable pH range of 3-9 and an optimized molar ratio 1: 2 for Fe: NTA. Results of quenching and chemical probe experiment together with results of electron paramagnetic resonance (EPR) analysis revealed the main reactive species of the system, including ^•OH, SO₄^•-, ¹O₂ and possibly Fe(IV). With the addition of NTA, the yields of ^•OH, SO₄^•-, ¹O₂ were enhanced. Results of mass spectrometry analysis and DFT calculations indicated the formation of 9 degradation byproducts of ACA via three primary degradation pathways such as hydroxyl substitution, carbonyl substitution, and decarboxylation. Furthermore, the EC/Fe(III)-NTA/PMS system could achieve excellent removal efficiency of ACA with different anions such as Cl^-, HCO₃^-, NO₃^- and H₂PO₄^- in the background. The practical applicability of the system was also verified with the high removal of commercial NAs mixture standard. Overall results have indicated the EC/Fe(III)-NTA/PMS system could be utilized for efficient reclamation of authentic oil and gas industrial wastewater under natural pH conditions.

Assessment of advanced xDH@B3LYP methods in describing various potential energy curves driven by π-π, CH/π, and SH/π non-bonded interactions

Accurate description of potential energy curves driven by non-bonded interactions remains a great challenge for pure density functional approximations (DFAs). It is because the R−6 decay behavior of dispersion cannot be intrinsically captured by the (semi)-local ingredients and the exact-exchange used in the popular hybrid DFAs. Overemphasizing the accuracy on the equilibrium region for the functional construction would likely deteriorate the overall performance on the other regions of potential energy surfaces. In consequence, the empirical dispersion correction becomes the standard component in DFAs to treat the non-bonded interactions. In this Letter, we demonstrate that without the use of empirical dispersion correction, doubly hybrid approximations, in particular two recently proposed revXYG3 and XYG7 functionals, hold the promise to have a balanced description of non-bonded interactions on the whole potential energy curves for several prototypes of π- π, CH/ π, and SH/ π interactions. The error of revXYG3 and XYG7 for non-bonded interactions is around 0.1 kcal/mol, and their potential energy curves almost coincide with the accurate CCSD(T)/CBS curves.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Doubly Hybrid Functionals Close to Chemical Accuracy for Both Finite and Extended Systems: Implementation and Test of XYG3 and XYGJ-OS

While being widely used to understand the chemical reactions in heterogeneous catalysis or other multidisciplinary systems, a great challenge that semilocal and hybrid density functional approximations (DFAs) are facing is to deliver a uniformly accurate description for both finite and extended systems. Herein, we perform reliable and well-converged periodic calculations of two doubly hybrid approximations (DHAs), XYG3 and XYGJ-OS, and demonstrate that the good accuracy of DHAs achieved for molecules is transferable to the semiconductors and insulators. Such an accuracy is not only for energetic properties but also for the first- and second-order response properties, which is general for different kinds of chemical environments, including simple cubic bulks, perovskite-type transition metal oxides like TiO2, and heterogeneous systems like CO adsorption on the NaCl(100) surface. The present finding has strengthened the predictive power of DFT, which not only will inspire the future development of the top-rung DFAs but also will boost their applications in multidisciplinary studies with high accuracy and efficiency.

Exploring the Limits of the XYG3-Type Doubly Hybrid Approximations for the Main-Group Chemistry: The xDH@B3LYP Model

Despite being the most accurate class of density functional approximations for the main-group chemistry, doubly hybrid approximations (DHAs) are generally considered to be incomplete in describing the medium- to long-range dispersive interactions. The existing DHAs are often supplemented with empirical long-range dispersion corrections. By using the extensive and chemically diverse GMTKN55 database, we explore the limits of the XYG3-type DHAs using the B3LYP reference orbitals, namely, xDH@B3LYP, with a gradually relaxed constraint on the mixing parameters of DHAs. Our results demonstrate that the xDH@B3LYP model can provide a balanced description of both covalent and noncovalent interactions with the accuracy and robustness comparable to or even better than the very expensive composite methods in wave function theory. Such an accuracy can be achieved without resorting to the use of any long-range correction scheme, shedding new light on the development of DHAs.

Polyvinylpyrrolidone-Coordinated Single-Site Platinum Catalyst Exhibits High Activity for Hydrogen Evolution Reaction

The essence of developing a Pt-based single-atom catalyst (SAC) for hydrogen evolution reaction (HER) is the preparation of well-defined and stable single Pt sites with desired electrocatalytic efficacy. Herein, we report a facile approach to generate uniformly dispersed Pt sites with outstanding HER performance via a photochemical reduction method using polyvinylpyrrolidone (PVP) molecules as the key additive to significantly simplify the synthesis and enhance the catalytic performance. The as-prepared catalyst displays remarkable kinetic activities (20 times higher current density than the commercially available Pt/C) with excellent stability (76.3 % of its initial activity after 5000 cycles) for HER. EXAFS measurements and DFT calculations demonstrate a synergetic effect, where the PVP ligands and the support together modulate the electronic structure of the Pt atoms, which optimize the hydrogen adsorption energy, resulting in a considerably improved HER activity.

Resolving the chemical identity of H2SO4 derived anions on Pt(111) electrodes: they're sulfate

Understanding how electrolyte composition controls electrocatalytic reactions requires molecular-level insight into electrode/electrolyte interaction. Perhaps the most basic aspect of this interaction, the speciation of the interfacial ion, is often controversial for even relatively simple systems. For example, for Pt(111) in 0.5 M H₂SO₄ it has long been debated whether the adsorbed anion is SO₄^2-, HSO₄^- or an H₃O⁺SO₄^2- ion pair. Here we apply interface-specific vibrational sum frequency (VSF) spectroscopy and theory to this problem and perform an isotope exchange study: we collect VSF spectra of Pt(111) in H₂SO₄(H₂O) and D₂SO₄(D₂O) as a function of bias and show that at all potentials they are identical. This is the most direct spectroscopic evidence to date that SO₄^2- is the dominant adsorbate, despite the fact that at 0.5 M H₂SO₄ bulk solution is dominated by HSO₄^-. This approach is based on the unique selection rule of the VSF spectroscopy and thus offers a new way of accessing general electrode/electrolyte interaction in electrocatalysis.

Simultaneous Attenuation of Both Self-Interaction Error and Nondynamic Correlation Error in Density Functional Theory: A Spin-Pair Distinctive Adiabatic-Connection Approximation

We present a spin-pair distinctive algorithm in the context of adiabatic-connection fluctuation-dissipation (ACFD) theorem, which enables to quantify the self-interaction error (SIE) and the nondynamic/strong correlation error (NCE) in the direct random-phase approximation (dRPA). Using this knowledge, we propose a spin-component scaled dRPA (scsRPA) correlation model with simultaneous attenuation of both the SIE and the NCE. Along with the exact exchange, scsRPA is shown to present a comprehensive improvement over dRPA, as well as the well-established PBE and PBE0 functionals, for bonding energies of pronounced multireference characters and transition-metal complexes of strongly correlated systems, while consistently providing an accurate description for reaction energies, reaction barriers, and noncovalent bond interactions of weakly correlated systems.

Accurate heats of formation of polycyclic saturated hydrocarbons predicted by using the XYG3 type of doubly hybrid functionals

Polycyclic saturated hydrocarbons (PSHs) are attractive candidates as hydrocarbon propellants. To assess their potential values, one of the key factors is to determine their energy contents, such as to calculate their heats of formation (HOF). In this work, we have calculated HOFs for a set of 36 PSHs including exo-Tricyclo[5.2.1.0^(2,6) ] decane, the principal component of the high-energy density hydrocarbon fuel commonly identified as JP-10. The results from B3LYP, B3LYP-D3BJ, M06-2X, B2PLYP, B2PLYP-D3BJ, and the XYG3 type of doubly hybrid (xDH) functionals are presented. It is demonstrated here that the xDH functionals yield accurate HOFs in good agreement with those from experiments or the G4 theory. In particular, XYGJ-OS, a low scaling xDH functional, is shown to hold the promise for accurate prediction of HOFs for PSHs of larger sizes. © 2018 Wiley Periodicals, Inc.

Understanding band gaps of solids in generalized Kohn-Sham theory

The fundamental energy gap of a periodic solid distinguishes insulators from metals and characterizes low-energy single-electron excitations. However, the gap in the band structure of the exact multiplicative Kohn-Sham (KS) potential substantially underestimates the fundamental gap, a major limitation of KS density-functional theory. Here, we give a simple proof of a theorem: In generalized KS theory (GKS), the band gap of an extended system equals the fundamental gap for the approximate functional if the GKS potential operator is continuous and the density change is delocalized when an electron or hole is added. Our theorem explains how GKS band gaps from metageneralized gradient approximations (meta-GGAs) and hybrid functionals can be more realistic than those from GGAs or even from the exact KS potential. The theorem also follows from earlier work. The band edges in the GKS one-electron spectrum are also related to measurable energies. A linear chain of hydrogen molecules, solid aluminum arsenide, and solid argon provide numerical illustrations.

Towards Efficient Orbital-Dependent Density Functionals for Weak and Strong Correlation

We present a new paradigm for the design of exchange-correlation functionals in density-functional theory. Electron pairs are correlated explicitly by means of the recently developed second order Bethe-Goldstone equation (BGE2) approach. Here we propose a screened BGE2 (sBGE2) variant that efficiently regulates the coupling of a given electron pair. sBGE2 correctly dissociates H_{2} and H_{2}^{+}, a problem that has been regarded as a great challenge in density-functional theory for a long time. The sBGE2 functional is then taken as a building block for an orbital-dependent functional, termed ZRPS, which is a natural extension of the PBE0 hybrid functional. While worsening the good performance of sBGE2 in H_{2} and H_{2}^{+}, ZRPS yields a remarkable and consistent improvement over other density functionals across various chemical environments from weak to strong correlation.

Reaching a Uniform Accuracy for Complex Molecular Systems: Long-Range-Corrected XYG3 Doubly Hybrid Density Functional

An unbiased understanding of complex molecular systems from first-principles critically demands theoretical methods with uniform accuracy for diverse interactions with different natures covering short-, medium-, and long-range correlations. Among the state-of-the-art density functional approximations (DFAs), doubly hybrid (DH) DFAs (e.g., XYG3 in this Letter) provide a remarkable improvement over the conventional DFAs (e.g., B3LYP in this Letter). Even though XYG3 works quite well in many cases of noncovalent bonding interactions (NCIs), it is incomplete in describing the pure long-range dispersive interactions. Here, we address such concerns by adding a scaled long-range contribution from the second-order perturbation theory (PT2). The long-range-corrected XYG3 (lrc-XYG3) is proposed without reparameterizing the three parameters in the original XYG3. Due to its overall excellent performance for all testing sets constructed for various purposes, lrc-XYG3 is the recommended method, which is expected to provide a balanced description of diverse interactions in complex molecular systems.

A fast doubly hybrid density functional method close to chemical accuracy using a local opposite spin ansatz

We develop and validate the XYGJ-OS functional, based on the adiabatic connection formalism and Görling-Levy perturbation theory to second order and using the opposite-spin (OS) ansatz combined with locality of electron correlation. XYGJ-OS with local implementation scales as N(3) with an overall accuracy of 1.28 kcal/mol for thermochemistry, bond dissociation energies, reaction barrier heights, and nonbonded interactions, comparable to that of 1.06 kcal/mol for the accurate coupled-cluster based G3 method (scales as N(7)) and much better than many popular density functional theory methods: B3LYP (4.98), PBE0 (4.36), and PBE (12.10).