The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C2 singlet ground state is known from correlated wave function theory to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B2, C2, O2, and F2 in the local spin density approximation (LSDA), the Perdew-Burke-Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C2. This work shows that spin symmetry breaking in singlet C2, which involves the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that underbinding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking and the insights into strong correlation that symmetry breaking can bring. The normally correlated low-lying triplet excited state has the right vertical excitation energy in SCAN but not in LSDA or PBE, where the triplet is a false ground state. Fractional occupation numbers are found only for the symmetry-unbroken singlet and only in LSDA and PBE GGA.
Polymer additives [like polyethylene oxide (PEO)] are widely used for smooth electrode deposition in aqueous zinc and a number of other battery systems currently investigated for energy storage applications. However, the precise mechanism by which they regulate morphology and suppress dendrite formation remains unclear. In this study, we address this knowledge gap by using in-situ electrochemical atomic force microscopy (EC-AFM) to directly observe the interfacial evolution during Zn electrodeposition and polymer adsorption on copper (Cu) substrates in the presence of varying concentrations of ZnSO4 and PEO. Contrary to previous literature assumptions which emphasize the binding to the growing Zn crystal surfaces or Zn2+ ions, our results demonstrate that PEO smooths Zn films by promoting nucleation of (002)-oriented Zn platelets through interactions with the Cu substrate. Density functional theory (DFT) simulations support this finding by showing that PEO adsorption on Cu modifies the interfacial energy of Zn/Cu/electrolyte interfaces, favoring the stabilization of Zn (002) on the Cu substrate, as well as confines Zn electrodeposition to a narrow near-surface region. These findings elucidate a novel design principle for electrode smoothing, emphasizing the importance of substrate selection paired with polymer additives that exhibit an attractive interaction with the substrate, but minimal interaction with growing crystals, offering a mechanistic perspective for improved battery performance.
Understanding the structure and composition of aluminate complexes in extremely alkaline systems such as Bayer liquors has received enormous attention due to their fundamental and industrial importance. However, obtaining direct molecular information of the underlying ion–ion interactions using traditional approaches such as NMR spectroscopy or Raman spectroscopy is challenging due to the weakness of these interactions and/or their complex overlapping spectral signatures. Here, we exploit in situ liquid secondary-ion mass spectrometry (SIMS) as a new approach and show how it enables new insights. In contrast with traditional techniques, using SIMS we succeeded in acquiring information on dominant ion clusters in these alkaline systems. In Na+/K+ mixed alkaline aluminate solutions, we clearly observe preferential formation of Na+-anion clusters over K+-anion clusters. Evaluation of these clusters by density functional theory (DFT) calculations shows that these structures are stable and that their relative bond energies are consistent with their observed SIMS signal intensity differences. This demonstrates a key advantage of in situ liquid SIMS for overcoming ambiguities obscuring important information in these systems on constituent molecular clusters defined by relatively weak ion-pair competition and ion–solvent interactions.
Abstract A procedure for defining virtual spaces, and the periodic one-electron and two-electron integrals, for plane-wave second quantized Hamiltonians has been developed and demonstrated using full configuration interaction (FCI) simulations and variational quantum eigensolver (VQE) circuits on Quantinuum's ion trap quantum computers accessed through Microsoft's Azure Quantum service. This work is an extension to periodic systems of a new class of algorithms in which the virtual spaces were generated by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. In this extension, the integration of the first Brillouin zone is automatically incorporated into the two-electron integrals. With these procedures we have been able to derive virtual spaces, containing only a few orbitals, that were able to capture a significant amount of correlation. The focus in this manuscript is on comparing the simulations of small molecules calculated with plane-wave basis sets with large periodic unit cells at the Γ-point, including images, to results for plane-wave basis sets with aperiodic unit cells. The results for this approach were promising as we were able to obtain good agreement between periodic and aperiodic results for an LiH molecule. Simulations performed on the Quantinuum H1-1 quantum computer were able to produce surprisingly good energies, reproducing the FCI values for the 1 COVO Hamiltonian to within 11 milliHartree (6.9 kcal/mol), when corrected for noise.
Solar photoexcitation of chromophoric groups in dissolved organic matter (DOM), when coupled to photoreduction of ubiquitous Fe(III)-oxide nanoparticles, can significantly accelerate DOM degradation in near-surface terrestrial systems, but the mechanisms of these reactions remain elusive. We examined the photolysis of chromophoric soil DOM coated onto hematite nanoplatelets featuring (001) exposed facets using a combination of molecular spectroscopies and density functional theory (DFT) computations. Reactive oxygen species (ROS) probed by electron paramagnetic resonance (EPR) spectroscopy revealed that both singlet oxygen and superoxide are the predominant ROS responsible for DOM degradation. DFT calculations confirmed that Fe(II) on the hematite (001) surface, created by interfacial electron transfer from photoexcited chromophores in DOM, can reduce dioxygen molecules to superoxide radicals (•O2–) through a one-electron transfer process. 1H nuclear magnetic resonance (NMR) and electrospray ionization Fourier-transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS) spectroscopies show that the association of DOM with hematite enhances the cleavage of aromatic groups during photodegradation. The findings point to a pivotal role for organic matter at the interface that guides specific ROS generation and the subsequent photodegradation process, as well as the prospect of using ROS signatures as a forensic tool to help interpret more complicated field-relevant systems.
Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.
We study symmetry breaking in the mean field solutions to the 2 electron hydrogen molecule within Kohn Sham (KS) local spin density function theory with Dirac exchange (the XLDA model). This simplified model shows behavior related to that of the (KS) spin density functional theory (SDFT) predictions in condensed and molecular systems. The Kohn Sham solutions to the constrained SDFT variation problem undergo spontaneous symmetry breaking as the relative strength of the non-convex exchange term increases. This results in the change of the molecular ground state from a paramagnetic state to an antiferromagnetic ground states and a stationary symmetric delocalized 1st excited state. We further characterize the limiting behavior of the minimizer when the strength of the exchange term goes to infinity. This leads to further bifurcations and highly localized states with varying character. The stability of the various solution classes is demonstrated by Hessian analysis. Finite element numerical results provide support for the formal conjectures.
For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. The use of virtual spaces derived from the one-electron operators has also been tried, and while some correlations are captured, the amount is quite low. To overcome these limitations, we have been developing new classes of algorithms to define virtual spaces by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, which are able to capture a significant amount of correlation. The focus in this manuscript is on using these derived basis sets to target full CI (FCI) quality results for H 2 on near-term quantum computers. However, the initial results for this approach were promising. We were able to obtain good agreement with FCI/cc-pVTZ results for this system with just 4 virtual orbitals, using both FCI and quantum simulations. The quality of the results using COVOs suggests that it may be possible to use them in other many-body approaches, including coupled cluster and Møller–Plesset perturbation theories, and open up the door to many-body calculations for pseudopotential plane-wave basis set methods.
