89 research outputs found
Simple Impurity Embedded in a Spherical Jellium: Approximations of Density Functional Theory compared to Quantum Monte Carlo Benchmarks
We study the electronic structure of a spherical jellium in the presence of a
central Gaussian impurity. We test how well the resulting inhomogeneity effects
beyond spherical jellium are reproduced by several approximations of density
functional theory (DFT). Four rungs of Perdew's ladder of DFT functionals,
namely local density approximation (LDA), generalized gradient approximation
(GGA), meta-GGA and orbital-dependent hybrid functionals are compared against
our quantum Monte Carlo (QMC) benchmarks. We identify several distinct
transitions in the ground state of the system as the electronic occupation
changes between delocalized and localized states. We examine the parameter
space of realistic densities () and moderate depths of the
Gaussian impurity (). The selected 18 electron system (with closed-shell
ground state) presents transitions while the 30 electron system
(with open-shell ground state) exhibits transitions. For the former
system, the accuracy for the transitions is clearly improving with increasing
sophistication of functionals with meta-GGA and hybrid functionals having only
small deviations from QMC. However, for the latter system, we find much larger
differences for the underlying transitions between our pool of DFT functionals
and QMC. We attribute this failure to treatment of the exact exchange within
these functionals. Additionally, we amplify the inhomogeneity effects by
creating the system with spherical shell which leads to even larger errors in
DFT approximations.Comment: 8 pages, 4 figures, submitted to PRB as a regular article revisited
version after revie
Pfaffian pairing wave functions in electronic structure quantum Monte Carlo
We investigate the accuracy of trial wave function for quantum Monte Carlo
based on pfaffian functional form with singlet and triplet pairing. Using a set
of first row atoms and molecules we find that these wave functions provide very
consistent and systematic behavior in recovering the correlation energies on
the level of 95%. In order to get beyond this limit we explore the
possibilities of multi-pfaffian pairing wave functions. We show that a small
number of pfaffians recovers another large fraction of the missing correlation
energy comparable to the larger-scale configuration interaction wave functions.
We also find that pfaffians lead to substantial improvements in fermion nodes
when compared to Hartree-Fock wave functions.Comment: 4 pages, 2 figures, 2 tables, submitted to PR
Approximate and exact nodes of fermionic wavefunctions: coordinate transformations and topologies
A study of fermion nodes for spin-polarized states of a few-electron ions and
molecules with one-particle orbitals is presented. We find exact nodes
for some cases of two electron atomic and molecular states and also the first
exact node for the three-electron atomic system in state using
appropriate coordinate maps and wavefunction symmetries. We analyze the cases
of nodes for larger number of electrons in the Hartree-Fock approximation and
for some cases we find transformations for projecting the high-dimensional node
manifolds into 3D space. The node topologies and other properties are studied
using these projections. We also propose a general coordinate transformation as
an extension of Feynman-Cohen backflow coordinates to both simplify the nodal
description and as a new variational freedom for quantum Monte Carlo trial
wavefunctions.Comment: 7 pages, 7 figure
Coherent "metallic" resistance and medium localisation in a disordered 1D insulator
It is believed, that a disordered one-dimensional (1D) wire with coherent
electronic conduction is an insulator with the mean resistance \simeq
e^{2L/\xi} and resistance dispersion \Delta_{\rho} \simeq e^{L/\xi}, where L is
the wire length and \xi is the electron localisation length. Here we show that
this 1D insulator undergoes at full coherence the crossover to a 1D "metal",
caused by thermal smearing and resonant tunnelling. As a result, \Delta_{\rho}
is smaller than unity and tends to be L/\xi - independent, while grows
with L/\xi first nearly linearly and then polynomially, manifesting the
so-called medium localisation.Comment: 4 pages, 4 figures, RevTeX
Noncovalent Interactions by QMC: Speedup by One-Particle Basis-Set Size Reduction
While it is empirically accepted that the fixed-node diffusion Monte-Carlo
(FN-DMC) depends only weakly on the size of the one-particle basis sets used to
expand its guiding functions, limits of this observation are not settled yet.
Our recent work indicates that under the FN error cancellation conditions,
augmented triple zeta basis sets are sufficient to achieve a benchmark level of
0.1 kcal/mol in a number of small noncovalent complexes. Here we report on a
possibility of truncation of the one-particle basis sets used in FN-DMC guiding
functions that has no visible effect on the accuracy of the production FN-DMC
energy differences. The proposed scheme leads to no significant increase in the
local energy variance, indicating that the total CPU cost of large-scale
benchmark noncovalent interaction energy FN-DMC calculations may be reduced.Comment: ACS book chapter, accepte
Understanding the apparent fractional charge of protons in the aqueous electrochemical double layer
A detailed atomic-scale description of the electrochemical interface is essential to the understanding of electrochemical energy transformations. In this work, we investigate the charge of solvated protons at the Pt(111) | H_2O and Al(111) | H_2O interfaces. Using semi-local density-functional theory as well as hybrid functionals and embedded correlated wavefunction methods as higher-level benchmarks, we show that the effective charge of a solvated proton in the electrochemical double layer or outer Helmholtz plane at all levels of theory is fractional, when the solvated proton and solvent band edges are aligned correctly with the Fermi level of the metal (E_F). The observed fractional charge in the absence of frontier band misalignment arises from a significant overlap between the proton and the electron density from the metal surface, and results in an energetic difference between protons in bulk solution and those in the outer Helmholtz plane
Understanding the apparent fractional charge of protons in the aqueous electrochemical double layer
A detailed atomic-scale description of the electrochemical interface is essential to the understanding of electrochemical energy transformations. In this work, we investigate the charge of solvated protons at the Pt(111) | H_2O and Al(111) | H_2O interfaces. Using semi-local density-functional theory as well as hybrid functionals and embedded correlated wavefunction methods as higher-level benchmarks, we show that the effective charge of a solvated proton in the electrochemical double layer or outer Helmholtz plane at all levels of theory is fractional, when the solvated proton and solvent band edges are aligned correctly with the Fermi level of the metal (E_F). The observed fractional charge in the absence of frontier band misalignment arises from a significant overlap between the proton and the electron density from the metal surface, and results in an energetic difference between protons in bulk solution and those in the outer Helmholtz plane
Acidic Oxygen Evolution Reaction Activity–Stability Relationships in Ru-Based Pyrochlores
Ru-based oxygen evolution reaction (OER) catalysts show significant promise for efficient water electrolysis, but rapid degradation poses a major challenge for commercial applications. In this work, we explore several Ru-based pyrochlores (A2Ru2O7, A = Y, Nd, Gd, Bi) as OER catalysts and demonstrate improved activity and stability of catalytic Ru sites relative to RuO2. Furthermore, we combine complementary experimental and theoretical analysis to understand how the A-site element impacts activity and stability under acidic OER conditions. Among the A2Ru2O7 studied herein, we find that a longer Ru−O bond and a weaker interaction of the Ru 4d and O 2p orbitals compared with RuO2 results in enhanced initial activity. We observe that the OER activity of the catalysts changes over time and is accompanied by both A-site and Ru dissolution at different relative rates depending on the identity of the A-site. Pourbaix diagrams constructed using density functional theory (DFT) calculations reveal a driving force for this experimentally observed dissolution, indicating that all compositions studied herein are thermodynamically unstable in acidic OER conditions. Theoretical activity predictions show consistent trends between A-site cation leaching and OER activity. These trends coupled with Bader charge analysis suggest that dissolution exposes highly oxidized Ru sites that exhibit enhanced activity. Overall, using the stability number (molO2 evolved/molRu dissolved) as a comparative metric, the A2Ru2O7 materials studied in this work show substantially greater stability than a standard RuO2 and commensurate stability to some Ir mixed metal oxides. The insights described herein provide a pathway to enhanced Ru catalyst activity and durability, ultimately improving the efficiency of water electrolyzers
Screening highly active perovskites for hydrogen-evolving reaction via unifying ionic electronegativity descriptor
[[abstract]]Facile and reliable screening of cost-effective, high-performance and scalable electrocatalysts is key for energy conversion technologies such as water splitting. ABO3-δ perovskites, with rich constitutions and structures, have never been designed via activity descriptors for critical hydrogen evolution reaction (HER). Here, we apply coordination rationales to introduce A-site ionic electronegativity (AIE) as an efficient unifying descriptor to predict the HER activities of 13 cobalt-based perovskites. Compared with A-site structural or thermodynamic parameter, AIE endows the HER activity with the best volcano trend. (Gd0.5La0.5)BaCo2O5.5+δ predicted from an AIE value of ~2.33 exceeds the state-of-the-art Pt/C catalyst in electrode activity and stability. X-ray absorption and computational studies reveal that the peak HER activities at a moderate AIE value of ~2.33 can be associated with the optimal electronic states of active B-sites via inductive effect in perovskite structure (~200 nm depth), including Co valence, Co-O bond covalency, band gap and O 2p-band position.[[notice]]補正完
- …