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.

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.

Improving XYG3-type doubly hybrid approximation using self-interaction corrected SCAN density and orbitals via the PZ-SIC framework: The xDH@SCAN(SIC) approach.

XYG3-type doubly hybrid (xDH) approximations have gained widespread recognition for their accuracy in describing a diverse range of chemical and physical interactions. However, a recent study [Song et al., J. Phys. Chem. Lett. 12, 800-807 (2021)] has highlighted the limitation of xDH methods in calculating the dissociation of NaCl molecules. This issue has been related to the density and orbitals used for evaluating the energy in xDH methods, which are obtained from lower-rung hybrid density functional approximations (DFAs) and display substantial density errors in the dissociation limit. In this work, we systematically investigate the influence of density on several challenging datasets and find that xDH methods are less sensitive to density errors compared to semi-local and hybrid DFAs. Furthermore, we demonstrate that the self-interaction corrected SCAN density approach offers superior accuracy compared to the self-consistent SCAN density and Hartree-Fock density approaches, as evidenced by performing charge analysis on the dissociation of heterodimers, such as NaCl and LiF. Building on these insights, we propose a five-parameter xDH method using the SCAN density and orbitals corrected by the PZ-SIC scheme. This new xDH@SCAN(SIC) method provides a balanced and accurate description across a wide range of challenging systems.

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.

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, SO4•-, and 1O2 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.

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 (H2S), 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 H2S 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 H2S effect was also examined in both in vitro and in vivo models where the cardioprotective effects of H2S were dependent on Cu/Zn-SOD.

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, SO4•-, 1O2 and possibly Fe(IV). With the addition of NTA, the yields of •OH, SO4•-, 1O2 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-, HCO3-, NO3- and H2PO4- 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.

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.

Identification of Water Hexamer on Cu(111) Surfaces.

While being considered as the building block of ice on a hydrophobic metal surface, the global minimum of the water hexamer is still elusive, which has impeded our understanding of water/metal interfaces. Herein, we comprehensively investigate water hexamer on Cu(111) theoretically and propose the boat configuration as the new in situ adsorption configuration from the scanning tunneling microscope experiments. All existing experimental measurements can therefore be well reproduced. Calculations in high-level theories reveal that the boat configuration is indeed the global minimum under experimental conditions, solving a long-standing discrepancy.

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 H2SO4 it has long been debated whether the adsorbed anion is SO42-, HSO4- or an H3O+SO42- 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 H2SO4(H2O) and D2SO4(D2O) as a function of bias and show that at all potentials they are identical. This is the most direct spectroscopic evidence to date that SO42- is the dominant adsorbate, despite the fact that at 0.5 M H2SO4 bulk solution is dominated by HSO4-. 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.

Massive-Parallel Implementation of the Resolution-of-Identity Coupled-Cluster Approaches in the Numeric Atom-Centered Orbital Framework for Molecular Systems.

We present a massive-parallel implementation of the resolution of identity (RI) coupled-cluster approach that includes single, double, and perturbatively triple excitations, namely, RI-CCSD(T), in the FHI-aims package for molecular systems. A domain-based distributed-memory algorithm in the MPI/OpenMP hybrid framework has been designed to effectively utilize the memory bandwidth and significantly minimize the interconnect communication, particularly for the tensor contraction in the evaluation of the particle-particle ladder term. Our implementation features a rigorous avoidance of the on-the-fly disk storage and excellent strong scaling of up to 10 000 and more cores. Taking a set of molecules with different sizes, we demonstrate that the parallel performance of our CCSD(T) code is competitive with the CC implementations in state-of-the-art high-performance-computing computational chemistry packages. We also demonstrate that the numerical error due to the use of RI approximation in our RI-CCSD(T) method is negligibly small. Together with the correlation-consistent numeric atom-centered orbital (NAO) basis sets, NAO-VCC-nZ, the method is applied to produce accurate theoretical reference data for 22 bio-oriented weak interactions (S22), 11 conformational energies of gaseous cysteine conformers (CYCONF), and 32 isomerization energies (ISO32).

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 the Nonplanarity in Aromatic Metallabenzenes: A σ-Control Mechanism.

Metallabenzenes, the organometallic counterparts of benzene with one of the C atoms being replaced by a metal atom, expand the family of aromatics and further create prospective candidates for novel applications as functional materials. One intriguing feature of these complexes is that their MC5 rings do not always constrict into a planar configuration as in the C6 ring of benzene. Such a deviation has often been attributed to the unfavorable antibonding interactions between an occupied metal d orbital and the π orbitals of the C5 moiety. We herein scrutinize the frontier orbital interactions in both σ and π spaces in a plethora of metallabenzene complexes using extensive density functional theory calculations. Unexpectedly, the nonplanarity in metallabenzenes is found to be hardly related to the π orbitals. It is the antibonding interaction between an occupied metal d orbital and the σ orbitals of the C5 moiety that dominates the observed distortion. Such a σ-control mechanism not only provides an explanation for the commonly observed nonplanarity in metallabenzenes but also points out a novel direction toward the rational design of functional materials with enhanced metalla-aromaticity.

Experimentally quantifying anion polarizability at the air/water interface.

The adsorption of large, polarizable anions from aqueous solution on the air/water interface controls important atmospheric chemistry and is thought to resemble anion adsorption at hydrophobic interfaces generally. While the favourability of adsorption of such ions is clear, quantifying adsorption thermodynamics has proven challenging because it requires accurate description of the structure of the anion and its solvation shell at the interface. In principle anion polarizability offers a structural window, but to the best of our knowledge there has so far been no experimental technique that allowed its characterization with interfacial specificity. Here, we meet this challenge using interface-specific vibrational spectroscopy of Cl-O vibrations of the [Formula: see text] anion at the air/water interface and report that the interface breaks the symmetry of the anion, the anisotropy of [Formula: see text]'s polarizability tensor is more than two times larger than in bulk water and concentration dependent, and concentration-dependent polarizability changes are consistent with correlated changes in surface tension.

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.

Insight into organic reactions from the direct random phase approximation and its corrections.

The performance of the random phase approximation (RPA) and beyond-RPA approximations for the treatment of electron correlation is benchmarked on three different molecular test sets. The test sets are chosen to represent three typical sources of error which can contribute to the failure of most density functional approximations in chemical reactions. The first test set (atomization and n-homodesmotic reactions) offers a gradually increasing balance of error from the chemical environment. The second test set (Diels-Alder reaction cycloaddition = DARC) reflects more the effect of weak dispersion interactions in chemical reactions. Finally, the third test set (self-interaction error 11 = SIE11) represents reactions which are exposed to noticeable self-interaction errors. This work seeks to answer whether any one of the many-body approximations considered here successfully addresses all these challenges.