Search

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE) providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation (BSE), variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support of GPU acceleration has been achieved with minor modifications of the GPAW code thanks to the CuPy library. We end the review with an outlook describing some future plans for GPAW.

How does an electrochemical interface respond to changes in the electrode potential? How does the response affect the key properties of the system - energetics, excess charge, capacitance? Essential questions key to ab-initio simulations of electrochemical systems, which we address in this work on the basis of a rigorous mathematical evaluation of the interfacial energetics at constant applied potential. By explicitly taking into account the configurational and electronic degrees of freedom we derive important statements about stationary points in the electronically grand canonical ensemble. We analyze their geometric response to changes in electrode potential and show that it can be mapped identically onto an additional contribution to the system's capacitance. We draw similar conclusions for the constant charge ensemble which equally allows to assess the respective stationary points. Our analysis of the relation between the canonical and grand canonical energetics reveals, however, one key difference between both ensembles. While the constant potential ensemble yields in general positive capacitances at local minima, the capacitance of local minima in the constant charge ensemble might become negative. We trace back this feature to the possibility of character switching of stationary points when switching between the ensembles causing the differences in the response to perturbations. Our systematical analysis not only provides a detailed qualitative and quantitative understanding of the interplay between electronic and configurational degrees of freedom and their contributions to the energetics of electrified interfaces but also highlights the similarities and subtle dissimilarities between the canonical and grand canonical description of the electronic degrees of freedom, which is crucial for a better understanding of theoretical calculations with and without potentiostat., Comment: 33 pages, 0 figures, 2 tables

We find that ion creation/destruction dominates the behavior of electrochemical reaction barriers, through grand-canonical electronic structure calculations of proton-deposition on transition metal surfaces. We show that barriers respond to potential in a nonlinear manner and trace this to the continuous degree of electron transfer as an ion is created or destroyed. This explains both Marcus-like curvature and Hammond-like shifts. Across materials, we find the barrier energy to be driven primarily by the charge presented on the surface, which in turn is dictated by the native work function, a fundamentally different driving force than non-electrochemical systems.

In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a first crucial step to evaluate the reliability of such computations. We discuss here general recommendations for verification studies aiming at further testing precision and transferability of density-functional-theory computational approaches and codes. We illustrate such recommendations using a greatly expanded protocol covering the whole periodic table from Z=1 to 96 and characterizing 10 prototypical cubic compounds for each element: 4 unaries and 6 oxides, spanning a wide range of coordination numbers and oxidation states. The primary outcome is a reference dataset of 960 equations of state cross-checked between two all-electron codes, then used to verify and improve nine pseudopotential-based approaches. Such effort is facilitated by deploying AiiDA common workflows that perform automatic input parameter selection, provide identical input/output interfaces across codes, and ensure full reproducibility. Finally, we discuss the extent to which the current results for total energies can be reused for different goals (e.g., obtaining formation energies)., Comment: Main text: 23 pages, 4 figures. Supplementary: 68 pages. Nature Review Physics 2023

We present a modified nudged elastic band routine that can reduce the number of force calls by more than 50% for bands with non-uniform convergence. The method, which we call "dyNEB", dynamically and selectively optimizes states based on the perpendicular forces and parallel spring forces acting on that region of the band. The convergence criteria are scaled to focus on the region of interest, i.e., the saddle point, while maintaining continuity of the band and avoiding truncation. We show that this method works well for solid state reaction barriers---non-electrochemical in general and electrochemical in particular---and that the number of force calls can be significantly reduced without loss of resolution at the saddle point.

Platinum is a nearly perfect catalyst for the hydrogen evolution reaction, and its high activity has conventionally been explained by its close-to-thermoneutral hydrogen binding energy (G~0). However, many candidate non-precious metal catalysts bind hydrogen with similar strengths, but exhibit orders-of-magnitude lower activity for this reaction. In this study, we employ electronic structure methods that allow fully potential-dependent reaction barriers to be calculated, in order to develop a complete working picture of hydrogen evolution on platinum. Through the resulting ab initio microkinetic models, we assess the mechanistic origins of Pt's high activity. Surprisingly, we find that the G~0 hydrogen atoms are kinetically inert, and that the kinetically active hydrogen atoms have G's much weaker, similar to that of gold. These on-top hydrogens have particularly low barriers, which we compare to those of gold, explaining the high reaction rates, and the exponential variations in coverages can uniquely explain Pt's strong kinetic response to the applied potential. This explains the unique reactivity of Pt that is missed by conventional Sabatier analyses, and suggests true design criteria for non-precious alternatives.

Azulene (Az) is a non-alternating, aromatic hydrocarbon composed of a five-membered, electron-rich and a seven-membered, electron-poor ring; an electron distribution that provides intrinsic redox activity. By varying the attachment points of the two electrode-bridging substituents to the Az centre, the influence of the redox functionality on charge transport is evaluated. The conductance of the 1,3 Az derivative is at least one order of magnitude lower than those of the 2,6 Az and 4,7 Az derivatives, in agreement with density functional theory (DFT) calculations. In addition, only 1,3 Az exhibits pronounced nonlinear current-voltage characteristics with hysteresis, indicating a bias-dependent conductance switching. DFT identifies the LUMO to be nearest to the Fermi energy of the electrodes, but to be an active transport channel only in the case of the 2,6 and the 4,7 Az derivatives, whereas the 1,3 Az derivative uses the HOMO at low and the LUMO+1 at high bias. In return, the localized, weakly coupled LUMO of 1,3 Az creates a slow electron-hopping channel responsible for the voltage-induced switching due to the occupation of a single MO., Comment: 6 pages, 3 figures

In our theoretical study where we combine a nonequilibrium Green's function (NEGF) approach with density functional theory (DFT) we investigate branched compounds containing ferrocene moieties in both branches which due to their metal centers are designed to allow for asymmetry induced by local charging. In these compounds the ferrocene moieties are connected to pyridyl anchor groups either directly or via acetylenic spacers in a meta-connection where we also compare our results with those obtained for the respective single-branched molecules with both meta- and para-connections between the metal center and the anchors. We find a destructive quantum interference (DQI) feature in the transmission function slightly below the lowest unoccupied molecular orbital (LUMO) which dominates the conductance even for the uncharged branched compound with spacer groups inserted. In an analysis based on mapping the structural characteristics of the range of molecules in our article onto tight-binding models, we identify the structural source of the DQI minimum as the through-space coupling between the pyridyl anchor groups. We also find that local charging on one of the branches only changes the conductance by about one order of magnitude which we explain in terms of the spatial distributions of the relevant molecular orbitals for the branched compounds., Comment: Accepted for publication in Physical Review B, 12 pages in printed version, 10 figures

Charge transport through single molecules can be influenced by the charge and spin states of redox-active metal centres placed in the transport pathway. These molecular intrinsic properties are usually addressed by varying the molecules electrochemical and magnetic environment, a procedure that requires complex setups with multiple terminals. Here we show that oxidation and reduction of organometallic compounds containing either Fe, Ru or Mo centres can solely be triggered by the electric field applied to a two-terminal molecular junction. Whereas all compounds exhibit bias-dependent hysteresis, the Mo-containing compound additionally shows an abrupt voltage-induced conductance switching, yielding high to low current ratios exceeding 1000 at voltage stimuli of less than 1.0 V. DFT calculations identify a localized, redox active molecular orbital that is weakly coupled to the electrodes and closely aligned with the Fermi energy of the leads because of the spin-polarised ground state unique to the Mo centre. This situation opens an additional slow and incoherent hopping channel for transport, triggering a transient charging effect of the entire molecule and a strong hysteresis with unprecedented high low-to-high current ratios., Comment: 9 pages, 4 figures

Besides active, functional molecular building blocks such as diodes or switches, passive components as, e.g., molecular wires, are required to realize molecular-scale electronics. Incorporating metal centers in the molecular backbone enables the molecular energy levels to be tuned in respect to the Fermi energy of the electrodes. Furthermore, by using more than one metal center and sp-bridging ligands, a strongly delocalized electron system is formed between these metallic "dopants", facilitating transport along the molecular backbone. Here, we study the influence of molecule--metal coupling on charge transport of dinuclear X(PP)$_2$FeC$_4$Fe(PP)$_2$X molecular wires (PP = Et$_2$PCH$_2$CH$_2$PEt$_2$); X = CN (1), NCS (2), NCSe (3), C$_4$SnMe$_3$ (4) and C$_2$SnMe$_3$ (5)) under ultra-high vacuum and variable temperature conditions. In contrast to 1 which showed unstable junctions at very low conductance ($8.1\cdot10^{-7}$ G$_0$), 4 formed a Au-C$_4$FeC$_4$FeC$_4$-Au junction 4' after SnMe$_3$ extrusion which revealed a conductance of $8.9\cdot10^{-3}$ G$_0$, three orders of magnitude higher than for 2 ($7.9\cdot10^{-6}$ G$_0$) and two orders of magnitude higher than for 3 ($3.8\cdot10^{-4}$ G$_0$). Density functional theory (DFT) confirmed the experimental trend in the conductance for the various anchoring motifs. The strong hybridization of molecular and metal states found in the C--Au coupling case enables the delocalized electronic system of the organometallic Fe$_2$ backbone to be extended over the molecule-metal interfaces to the metal electrodes to establish high-conductive molecular wires.

For defining the conductance of single molecule junctions with a redox functionality in an electrochemical cell, two conceptually different electron transport mechanisms, namely coherent tunnelling and vibrationally induced hopping compete with each other, where implicit parameters of the setup such as the length of the molecule and the applied gate voltage decide which mechanism is the dominant one. Although coherent tunnelling is most efficiently described within Landauer theory, while the common theoretical treatment of electron hopping is based on Marcus theory, both theories are adequate for the processes they describe without introducing accuracy limiting approximations. For a direct comparison, however, it has to be ensured that the crucial quantities obtained from electronic structure calculations, i.e. the transmission function T(E) in Landauer theory, and the transfer integral V, the reorganisation energy $\lambda$ and the driving force $\Delta G^0$ in Marcus theory, are derived from similar grounds as pointed out by Nitzan and co-workers in a series of publications. In this article our framework is a single particle picture, where we perform density functional theory calculations for the conductance corresponding to both transport mechanisms for junctions with the central molecule containing one, two or three Ruthenium centers, respectively, where we extrapolate our results in order to define the critical length for the transition point of the two regimes which we identify at 5.76 nm for this type of molecular wire. We also discuss trends in dependence on an electrochemically induced gate potential.

There are various quantum chemical approaches for an ab initio description of transfer integrals within the framework of Marcus theory in the context of electron transfer reactions. In our article we aim to calculate transfer integrals in redox-active single molecule junctions, where we focus on the coherent tunneling limit with the metal leads taking the position of donor and acceptor and the molecule acting as a transport mediating bridge. This setup allows us to derive a conductance, which can be directly compared with recent results from a non-equilibrium Green's function approach. Compared with purely molecular systems we face additional challenges due to the metallic nature of the leads, which rules out some of the common techniques, and due to their periodicity, which requires {\bf k} space integration. We present three different methods, all based on density functional theory, for calculating the transfer integral under these constraints, which we benchmark on molecular test systems from the relevant literature. We also discuss manybody effects and apply all three techniques to a junction with a Ruthenium complex in different oxidation states.

For adjusting the charging state of a molecular metal complex in the context of a density functional theory description of coherent electron transport through single molecule junctions, we correct for self interaction effects by fixing the charge on a counterion, which in our calculations mimics the effect of the gate in an electrochemical STM setup, with two competing methods, namely the generalized $\Delta$ SCF technique and screening with solvation shells. One would expect a transmission peak to be pinned at the Fermi energy for a nominal charge of +1 on the molecule in the junction but we find a more complex situation in this multicomponent system defined by the complex, the leads, the counterion and the solvent. In particular equilibrium charge transfer between the molecule and the leads plays an importanty role, which we investigate in dependence on the total external charge in the context of electronegativity theory.