BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Patrick Rinke

Patrick Rinke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.chem-phcond-mat.otherphysics.comp-phcond-mat.str-elMachine Learningcond-mat.mes-hallphysics.atm-clusphysics.optics

Catalog footprint

What is connected

35works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2023arXiv

Single-atom dopants in plasmonic nanocatalysts

Bimetallic nanostructures combining plasmonic and catalytic metals are promising for tailoring and enhancing plasmonic hot-carrier generation utilized in plasmonic catalysis. In this work, we study the plasmonic hot-carrier generation in noble metal nanoparticles (Ag, Au, Cu) with single-atom dopants (Ag, Au, Cu, Pd, Pt) with first-principles time-dependent density-functional theory calculations. Our results show that the local hot-carrier generation at the dopant atom is significantly altered by the dopant element while the plasmonic response of the nanoparticle as a whole is not significantly affected. In particular, hot holes at the dopant atom originate from the discrete d-electron states of the dopant. The energies of these d-electron states, and hence those of the hot holes, depend on the dopant element, which opens up the possibility to tune hot-carrier generation with suitable dopants.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Accurate computational prediction of core-electron binding energies in carbon-based materials: A machine-learning model combining density-functional theory and $\boldsymbol{GW}$

We present a quantitatively accurate machine-learning (ML) model for the computational prediction of core-electron binding energies, from which x-ray photoelectron spectroscopy (XPS) spectra can be readily obtained. Our model combines density functional theory (DFT) with $GW$ and uses kernel ridge regression for the ML predictions. We apply the new approach to materials and molecules containing carbon, hydrogen and oxygen, and obtain qualitative and quantitative agreement with experiment, resolving spectral features within 0.1 eV of reference experimental spectra. The method only requires the user to provide a structural model for the material under study to obtain an XPS prediction within seconds. Our new tool is freely available online through the XPS Prediction Server.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A balanced treatment of static and dynamic correlation in free- and Mg-porphyrin

We present an ab-initio dynamical configuration interaction (DCI) study of free- and Mg-porphyrin. DCI is a recently developed active space theory based on the Löwdin downfolding technique. In the active space, static correlation is described exactly with full configuration interaction. In the high energy, dynamically correlated subspace, we treat correlation at the quasiparticle level in the $GW$ approximation of Green's function theory. The final theory combines wave function and Green's function methods to give a balanced description of static and dynamic correlation. The theory and algorithm give a multireference treatment of ground and excited states for low computational cost. The four orbital Gouterman model of porphyrin offers an ideal active space in a large, correlated system to test the cost and accuracy of the embedding for large systems. Our parameter free, fully ab-initio DCI calculations in the minimal four-level active space agree well with both experiment and more expensive benchmark theories for the $Q_x$ and $Q_y$ transitions of free- and Mg-porphyrin. Studying the convergence of the excitation energies suggests that DCI correctly captures size extensive correlation effects, making it a promising active space theory for large, strongly-correlated systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Atomic structures and orbital energies of 61,489 crystal-forming organic molecules

Data science and machine learning in materials science require large datasets of technologically relevant molecules or materials. Currently, publicly available molecular datasets with realistic molecular geometries and spectral properties are rare. We here supply a diverse benchmark spectroscopy dataset of 61,489 molecules extracted from organic crystals in the Cambridge Structural Database (CSD), denoted OE62. Molecular equilibrium geometries are reported at the Perdew-Burke-Ernzerhof (PBE) level of density functional theory (DFT) including van der Waals corrections for all 62k molecules. For these geometries, OE62 supplies total energies and orbital eigenvalues at the PBE and the PBE hybrid (PBE0) functional level of DFT for all 62k molecules in vacuum as well as at the PBE0 level for a subset of 30,876 molecules in (implicit) water. For 5,239 molecules in vacuum, the dataset provides quasiparticle energies computed with many-body perturbation theory in the $G_0W_0$ approximation with a PBE0 starting point (denoted GW5000 in analogy to the GW100 benchmark set (M. van Setten et al. J. Chem. Theory Comput. 12, 5076 (2016))).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Boron doping in gallium oxide from first principles

We study the feasibility of boron doping in gallium oxide ($\text{Ga}_2\text{O}_3$) for neutron detection. $\text{Ga}_2\text{O}_3$ is a wide band-gap, radiation hard material which has potential for neutron detection if it can be doped with a neutron active element. We investigate the boron-10 isotope as a possible neutron active dopant. Intrinsic and boron induced defects in $\text{Ga}_2\text{O}_3$ are studied with semi-local and hybrid density-functional-theory calculations. We find that in growth conditions favourable for boron, boron substitutional defects are likely to form making boron doping of $\text{Ga}_2\text{O}_3$ feasible.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient Cysteine Conformer Search with Bayesian Optimization

Finding low-energy molecular conformers is challenging due to the high dimensionality of the search space and the computational cost of accurate quantum chemical methods for determining conformer structures and energies. Here, we combine active-learning Bayesian optimization (BO) algorithms with quantum chemistry methods to address this challenge. Using cysteine as an example, we show that our procedure is both efficient and accurate. After only one thousand single-point calculations and approximately thirty structure relaxations, which is less than 10% computational cost of the current fastest method, we have found the low-energy conformers in good agreement with experimental measurements and reference calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient hyperparameter tuning for kernel ridge regression with Bayesian optimization

Machine learning methods usually depend on internal parameters -- so called hyperparameters -- that need to be optimized for best performance. Such optimization poses a burden on machine learning practitioners, requiring expert knowledge, intuition or computationally demanding brute-force parameter searches. We here address the need for more efficient, automated hyperparameter selection with Bayesian optimization. We apply this technique to the kernel ridge regression machine learning method for two different descriptors for the atomic structure of organic molecules, one of which introduces its own set of hyperparameters to the method. We identify optimal hyperparameter configurations and infer entire prediction error landscapes in hyperparameter space, that serve as visual guides for the hyperparameter dependence. We further demonstrate that for an increasing number of hyperparameters, Bayesian optimization becomes significantly more efficient in computational time than an exhaustive grid search -- the current default standard hyperparameter search method -- while delivering an equivalent or even better accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electron-phonon coupling in d-electron solids: A temperature dependent study of rutile TiO2 by first-principles theory and two-photon photoemission

Rutile TiO2 is a paradigmatic transition metal oxide with applications in optics, electronics, photocatalysis, etc., that are subject to pervasive electron-phonon interaction. To understand how energies of its electronic bands, and in general semiconductors or metals where the frontier orbitals have a strong d-band character, depend on temperature, we perform a comprehensive theoretical and experimental study of the effects of electron-phonon (e-p) interactions. In a two-photon photoemission (2PP) spectroscopy study we observe an unusual temperature dependence of electronic band energies within the conduction band of reduced rutile TiO2, which is contrary to the well understood sp-band semiconductors and points to a so far unexplained dichotomy in how the e-p interactions affect differently the materials where the frontier orbitals are derived from the sp- and d-orbitals. To develop a broadly applicable model, we employ state-of-the-art first-principles calculations that explain how phonons promote interactions between the Ti-3d orbitals of the conduction band within the octahedral crystal field. The characteristics differences in e-p interactions experienced by the Ti 3d-orbitals of rutile TiO2 crystal lattice are contrasted with the more familiar behavior of the Si 2s-orbitals of stishovite SiO2 polymorph, in which the frontier 2s-orbital experience a similar crystal field with the opposite effect...

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Projective Preferential Bayesian Optimization

Bayesian optimization is an effective method for finding extrema of a black-box function. We propose a new type of Bayesian optimization for learning user preferences in high-dimensional spaces. The central assumption is that the underlying objective function cannot be evaluated directly, but instead a minimizer along a projection can be queried, which we call a projective preferential query. The form of the query allows for feedback that is natural for a human to give, and which enables interaction. This is demonstrated in a user experiment in which the user feedback comes in the form of optimal position and orientation of a molecule adsorbing to a surface. We demonstrate that our framework is able to find a global minimum of a high-dimensional black-box function, which is an infeasible task for existing preferential Bayesian optimization frameworks that are based on pairwise comparisons.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Relativistic correction scheme for core-level binding energies from $GW$

We present a relativistic correction scheme to improve the accuracy of 1s core-level binding energies calculated from Green's function theory in the $GW$ approximation, which does not add computational overhead. An element-specific corrective term is derived as the difference between the 1s eigenvalues obtained from the self-consistent solutions to the non- or scalar-relativistic Kohn-Sham equations and the four-component Dirac-Kohn-Sham equations for a free neutral atom. We examine the dependence of this corrective term on the molecular environment and on the amount of exact exchange in hybrid exchange-correlation functionals. This corrective term is then added as a perturbation to the quasiparticle energies from partially self-consistent and single-shot $GW$ calculations. We show that this element-specific relativistic correction, when applied to a previously reported benchmark set of 65 core-state excitations [J. Phys. Chem. Lett. 11, 1840 (2020)], reduces the mean absolute error (MAE) with respect to experiment from 0.55 to 0.30 eV and eliminates the species dependence of the MAE, which otherwise increases with the atomic number. The relativistic corrections also reduce the species dependence for the optimal amount of exact exchange in the hybrid functional used as starting point for the single-shot $G_0W_0$ calculations. Our correction scheme can be transferred to other methods, which we demonstrate for the Delta self-consistent field ($Δ$SCF) approach based on density functional theory.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Accurate absolute and relative core-level binding energies from $GW$

We present an accurate approach to compute X-ray photoelectron spectra based on the $GW$ Green's function method, that overcomes shortcomings of common density functional theory approaches. $GW$ has become a popular tool to compute valence excitations for a wide range of materials. However, core-level spectroscopy is thus far almost uncharted in $GW$. We show that single-shot perturbation calculations in the $G_0W_0$ approximation, which are routinely used for valence states, cannot be applied for core levels and suffer from an extreme, erroneous transfer of spectral weight to the satellite spectrum. The correct behavior can be restored by partial self-consistent $GW$ schemes or by using hybrid functionals with almost 50% of exact exchange as starting point for $G_0W_0$. We include also relativistic corrections and present a benchmark study for 65 molecular 1s excitations. Our absolute and relative $GW$ core-level binding energies agree within 0.3 and 0.2 eV with experiment, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

DScribe: Library of Descriptors for Machine Learning in Materials Science

DScribe is a software package for machine learning that provides popular feature transformations ("descriptors") for atomistic materials simulations. DScribe accelerates the application of machine learning for atomistic property prediction by providing user-friendly, off-the-shelf descriptor implementations. The package currently contains implementations for Coulomb matrix, Ewald sum matrix, sine matrix, Many-body Tensor Representation (MBTR), Atom-centered Symmetry Function (ACSF) and Smooth Overlap of Atomic Positions (SOAP). Usage of the package is illustrated for two different applications: formation energy prediction for solids and ionic charge prediction for atoms in organic molecules. The package is freely available under the open-source Apache License 2.0.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Atomic structure of metal-halide perovskites from first principles: The chicken-and-egg paradox of the organic-inorganic interaction

We have studied the prototype hybrid organic-inorganic perovskite CH3NH3PbI3 and its three close relatives, CH3NH3SnI3, CH3NH3PbCl3 and CsPbI3, using relativistic density function theory. The long-range van der Waals (vdW) interactions were incorporated into the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional using the Tkatchenko-Scheffler pairwise scheme. Our results reveal that hydrogen bonding, which is well described by the PBE functional, plays a decisive role for the structural parameters of these systems, including the position and orientation of the organic cation as well as the deformation of the inorganic framework. The magnitude of the inorganic-framework deformation depends sensitively on the orientation of the organic cation, and directly influences the stability of the hybrid perovskites. Our results suggest that the organic and the inorganic components complement each other: The low symmetry of the organic cation is the origin of the inorganic-framework deformation, which then aids the overall stabilization of the hybrid perovskite structure. This stabilization is indirectly affected by vdW interactions, which lead to smaller unit-cell volumes than in PBE and therefore modulate the interaction between the organic cation and the inorganic framework. The vdW-induced lattice-constant corrections are system dependent and lead to PBE+vdW lattice constants in good agreement with experiment. Further insight is gained by analysing the vdW contributions. In all iodide-based hybrid perovskites the interaction between the organic cation and the iodide anions provides the largest lattice-constant change, followed by iodine-iodine and the organic cation - heavy-metal cation interaction. These corrections follow an almost linear dependence on the lattice constant within the range considered in our study, and are therefore approximately additive.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

First-principles calculations of indirect Auger recombination in nitride semiconductors

Auger recombination is an important non-radiative carrier recombination mechanism in many classes of optoelectronic devices. The microscopic Auger processes can be either direct or indirect, mediated by an additional scattering mechanism such as the electron-phonon interaction. Phonon-assisted Auger recombination is particularly strong in nitride materials and affects the efficiency of nitride optoelectronic devices at high powers. Here we present a first-principles computational formalism for the study of direct and indirect Auger recombination in direct-band-gap semiconductors and apply it to the case of nitride materials. We show that direct Auger recombination is weak in the nitrides and cannot account for experimental measurements. On the other hand, carrier scattering by phonons and alloy disorder enables indirect Auger processes that can explain the observed loss in devices. We analyze the dominant phonon contributions to the Auger recom- bination rate and the influence of temperature and strain on the values of the Auger coefficients. Auger processes assisted by charged-defect scattering are much weaker than the phonon-assisted ones for realistic defect densities and not important for the device performance. The computational formalism is general and can be applied to the calculation of the Auger coefficient in other classes of optoelectronic materials.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Length Dependence of Ionization Potentials of Trans-Acetylenes: Internally-Consistent DFT/GW Approach

We follow the evolution of the Ionization Potential (IP) for the paradigmatic quasi-one-dimensional trans-acetylene family of conjugated molecules, from short to long oligomers and to the infinite polymer trans-poly-acetylene (TPA). Our results for short oligomers are very close to experimental available data. We find that the IP varies with oligomer length and converges to the given value for TPA with a smooth, coupled inverse-length-exponential behavior. Our prediction is based on an "internally-consistent" scheme to adjust the exchange mixing parameter $α$ of the PBEh hybrid density functional, so as to obtain a description of the electronic structure consistent with the quasiparticle approximation for the IP. This is achieved by demanding that the corresponding quasiparticle correction, in the GW@PBEh approximation, vanishes for the IP when evaluated at PBEh($α^{ic}$). We find that $α^{ic}$ is also system-dependent and converges with increasing oligomer length, allowing to capture the dependence of IP and other electronic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Piecewise linearity in the $GW$ approximation for accurate quasiparticle energy predictions

We identify the deviation from the straight line error (DSLE) -- i.e., the spurious non-linearity of the total energy as a function of fractional particle number -- as the main source for the discrepancy between experimental vertical ionization energies and theoretical quasiparticle energies, as obtained from the $GW$ and $GW$+SOSEX approximations to many-body perturbation theory (MBPT). For self-consistent calculations, we show that $GW$ suffers from a small DSLE. Conversely, for perturbative $G_0W_0$ and $G_0W_0$+SOSEX calculations the DSLE depends on the starting point. We exploit this starting-point dependence to reduce (or completely eliminate) the DSLE. We find that the agreement with experiment increases as the DSLE reduces. DSLE-minimized schemes, thus, emerge as promising avenues for future developments in MBPT.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self-Consistent Green Function Embedding for Advanced Electronic Structure Methods Based on a Dynamical Mean-Field Concept

We present an embedding scheme for periodic systems that facilitates the treatment of the physically important part (here the unit cell) with advanced electronic-structure methods, that are computationally too expensive for periodic systems. The rest of the periodic system is treated with computationally less demanding approaches, e.g., Kohn-Sham density-functional theory, in a self- consistent manner. Our scheme is based on the concept of dynamical mean-field theory (DMFT) formulated in terms of Green functions. In contrast to the original DMFT formulation for correlated model Hamiltonians, we here consider the unit cell as local embedded cluster in a first-principles way, that includes all electronic degrees of freedom. Our real-space dynamical mean-field embedding (RDMFE) scheme features two nested Dyson equations, one for the embedded cluster and another for the periodic surrounding. The total energy is computed from the resulting Green functions. The performance of our scheme is demonstrated by treating the embedded region with hybrid functionals and many-body perturbation theory in the GW approach for simple bulk systems. The total energy and the density of states converge rapidly with respect to the computational parameters and approach their bulk limit with increasing cluster (i.e., unit cell) size.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Static correlation and electron localization in molecular dimers from the self-consistent RPA and GW approximation

We investigate static correlation and delocalization errors in the self-consistent GW and random-phase approximation (RPA) by studying molecular dissociation of the H_2 and LiH molecules. Although both approximations contain topologically identical diagrams, the non-locality and frequency dependence of the GW self-energy crucially influence the different energy contributions to the total energy as compared to the use of a static local potential in the RPA. The latter leads to significantly larger correlation energies which allow for a better description of static correlation at intermediate bond distances. The substantial error found in GW is further analyzed by comparing spin-restricted and spin-unrestricted calculations. At large but finite nuclear separation their difference gives an estimate of the so-called fractional spin error normally determined only in the dissociation limit. Furthermore, a calculation of the dipole moment of the LiH molecule at dissociation reveals a large delocalization error in GW making the fractional charge error comparable to the RPA. The analyses are supplemented by explicit formulae for the GW Green's function and total energy of a simplified two-level model providing additional insights into the dissociation limit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Why graphene growth is very different on the C face than on the Si face of SiC: Insights from surface equilibria and the (3$\times$3)-3C-SiC($\bar{\text{1}}\bar{\text{1}}\bar{\text{1}}$) reconstruction

We address the stability of the surface phases that occur on the C-side of 3C-SiC($\bar{1} \bar{1} \bar{1}$) at the onset of graphene formation. In this growth range, experimental reports reveal a coexistence of several surface phases. This coexistence can be explained by a Si-rich model for the unknown (3$\times$3) reconstruction, the known (2$\times$2)$_{C}$ adatom phase, and the graphene covered (2$\times$2)$_{C}$ phase. By constructing an $ab$ $initio$ surface phase diagram using a van der Waals corrected density functional, we show that the formation of a well defined interface structure like the "buffer-layer" on the Si side is blocked by Si-rich surface reconstructions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A multi-scale approach to the electronic structure of doped semiconductor surfaces

The inclusion of the global effects of semiconductor doping poses a unique challenge for first-principles simulations, because the typically low concentration of dopants renders an explicit treatment intractable. In systems which do not display long-range band bending, a satisfactory remedy is offered by the use of "pseudoatoms", with a fractional nuclear charge matching the bulk doping concentration. However, this alone is not always sufficient for charged surfaces, where the width of the space-charge region (SCR) often exceeds realistic supercell dimensions. One generalization of the pseudoatom approach which overcomes this difficulty relies on the introduction of an artificially high doping level within a slab calculation, in conjunction with a multi-scale electrostatic energy correction. Here, we present an alternative technique that naturally extends the pseudoatom approach while bypassing the need for calculations with an unrealistically high doping level. It is based on the introduction of excess charge, mimicking free charge carriers from the SCR, along with a fixed sheet of counter-charge mimicking the SCR-related field. Self-consistency is obtained by imposing charge conservation and Fermi level equilibration between the bulk, treated semi-classically, and the electronic states of the slab/surface, which are treated quantum-mechanically. The method, which we call CREST - the Charge-Reservoir Electrostatic Sheet Technique - can be used with standard electronic structure codes. We validate the approach using a simple tight-binding model, which allows for comparison of its results with calculations encompassing the full SCR explicitly. We then employ it with density functional theory, where it is used to obtain insights into the electronic structures of the "clean-cleaved" Si(111) surface and its buckled (2x1) reconstruction, at various doping densities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First-Principles Description of Charge Transfer in Donor-Acceptor Compounds from Self-Consistent Many-Body Perturbation Theory

We investigate charge transfer in prototypical molecular donor-acceptor compounds using hybrid density functional theory (DFT) and the GW approximation at the perturbative level (G0W0) and at full self-consistency (sc-GW). For the systems considered here, no charge transfer should be expected at large intermolecular separation according to photoemission experiment and accurate quantum-chemistry calculations. The capability of hybrid exchange-correlation functionals of reproducing this feature depends critically on the fraction of exact exchange $α$, as for small values of $α$ spurious fractional charge transfer is observed between the donor and the acceptor. G0W0 based on hybrid DFT yields the correct alignment of the frontier orbitals for all values of $α$. However, G0W0 has no capacity to alter the ground-state properties of the system, because of its perturbative nature. The electron density in donor-acceptor compounds thus remains incorrect for small $α$ values. In sc-GW, where the Green's function is obtained from the iterative solution of the Dyson equation, the electron density is updated and reflects the correct description of the level alignment at the GW level, demonstrating the importance of self-consistent many-body approaches for the description of ground- and excited-state properties in donor-acceptor systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bond Breaking and Bond Formation: How Electron Correlation is Captured in Many-Body Perturbation Theory and Density-Functional Theory

For the paradigmatic case of H2-dissociation we compare state-of-the-art many-body perturbation theory (MBPT) in the GW approximation and density-functional theory (DFT) in the exact-exchange plus random-phase approximation for the correlation energy (EX+cRPA). For an unbiased comparison and to prevent spurious starting point effects both approaches are iterated to full self-consistency (i.e. sc-RPA and sc-GW). The exchange-correlation diagrams in both approaches are topologically identical, but in sc-RPA they are evaluated with non-interacting and in sc-GW with interacting Green functions. This has a profound consequence for the dissociation region, where sc-RPA is superior to sc-GW. We argue that for a given diagrammatic expansion, the DFT framework outperforms the many-body framework when it comes to bond-breaking. We attribute this to the difference in the correlation energy rather than the treatment of the kinetic energy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interface dipoles of organic molecules on Ag(111) in hybrid density-functional theory

We investigate the molecular acceptors 3,4,9,10-perylene-tetracarboxylic acid dianhydride (PTCDA), 2,3,5,6-tetra uoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ), and 4,5,9,10-pyrenetetraone (PYTON) on Ag(111) using densityfunctional theory. For two groups of the HSE(α, ω) family of exchange-correlation functionals (ω= 0 and ω= 0.2Å) we study the isolated components as well as the combined systems as a function of the amount of exact-exchange (α). We find that hybrid functionals favour electron transfer to the adsorbate. Comparing to experimental work-function data, we report for (α) ca. 0.25 a notable but small improvement over (semi)local functionals for the interface dipole. Although Kohn-Sham eigenvalues are only approximate representations of ionization energies, incidentally, at this value also the density of states agrees well with the photoelectron spectra. However, increasing (α) to values for which the energy of the lowest unoccupied molecular orbital matches the experimental electron affinity in the gas phase worsens both the interface dipole and the density of states. Our results imply that semi-local DFT calculations may often be adequate for conjugated organic molecules on metal surfaces and that the much more computationally demanding hybrid functionals yield only small improvements.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Renormalized Second-order Perturbation Theory for The Electron Correlation Energy: Concept, Implementation, and Benchmarks

We present a renormalized second-order perturbation theory (rPT2), based on a Kohn-Sham (KS) reference state, for the electron correlation energy that includes the random-phase approximation (RPA), second-order screened exchange (SOSEX), and renormalized single excitations (rSE). These three terms all involve a summation of certain types of diagrams to infinite order, and can be viewed as "renormalization" of the 2nd-order direct, exchange, and single excitation (SE) terms of Rayleigh-Schrödinger perturbation theory based on an KS reference. In this work we establish the concept of rPT2 and present the numerical details of our SOSEX and rSE implementations. A preliminary version of rPT2, in which the renormalized SE (rSE) contribution was treated approximately, has already been benchmarked for molecular atomization energies and chemical reaction barrier heights and shows a well balanced performance [Paier et al, New J. Phys. 14, 043002 (2012)]. In this work, we present a refined version of rPT2, in which we evaluate the rSE series of diagrams rigorously. We then extend the benchmark studies to non-covalent interactions, including the rare-gas dimers, and the S22 and S66 test sets. Despite some remaining shortcomings, we conclude that rPT2 gives an overall satisfactory performance across different chemical environments, and is a promising step towards a generally applicable electronic structure approach.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Self-consistent GW: All-electron implementation with localized basis functions

This paper describes an all-electron implementation of the self-consistent GW (sc-GW) approach -- i.e. based on the solution of the Dyson equation -- in an all-electron numeric atom-centered orbital (NAO) basis set. We cast Hedin's equations into a matrix form that is suitable for numerical calculations by means of i) the resolution of identity technique to handle 4-center integrals; and ii) a basis representation for the imaginary-frequency dependence of dynamical operators. In contrast to perturbative G0W0, sc-GW provides a consistent framework for ground- and excited-state properties and facilitates an unbiased assessment of the GW approximation. For excited-states, we benchmark sc-GW for five molecules relevant for organic photovoltaic applications: thiophene, benzothiazole, 1,2,5-thiadiazole, naphthalene, and tetrathiafulvalene. At self-consistency, the quasi-particle energies are found to be in good agreement with experiment and, on average, more accurate than G0W0 based on Hartree-Fock (HF) or density-functional theory with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. Based on the Galitskii-Migdal total energy, structural properties are investigated for a set of diatomic molecules. For binding energies, bond lengths, and vibrational frequencies sc-GW and G0W0 achieve a comparable performance, which is, however, not as good as that of exact-exchange plus correlation in the random-phase approximation (EX+cRPA) and its advancement to renormalized second-order perturbation theory (rPT2). Finally, the improved description of dipole moments for a small set of diatomic molecules demonstrates the quality of the sc-GW ground state density.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermodynamic equilibrium conditions of graphene films on SiC

First-principles surface phase diagrams reveal that epitaxial monolayer graphene films on the Si side of 3C-SiC(111) can exist as thermodynamically stable phases in a narrow range of experimentally controllable conditions, defining a path to the highest-quality graphene films. Our calculations are based on a van der Waals corrected density functional. The full, experimentally observed (6 sqrt(3)x 6 sqrt(3))-R30 supercells for zero- to trilayer graphene are essential to describe the correct interface geometries and the relative stability of surface phases and possible defects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Benchmark of GW Methods for Azabenzenes: Is the GW Approximation Good Enough?

Many-body perturbation theory in the GW approximation is a useful method for describing electronic properties associated with charged excitations. A hierarchy of GW methods exists, starting from non-self-consistent G0W0, through partial self-consistency in the eigenvalues (ev-scGW) and in the Green function (scGW0), to fully self-consistent GW (scGW). Here, we assess the performance of these methods for benzene, pyridine, and the diazines. The quasiparticle spectra are compared to photoemission spectroscopy (PES) experiments with respect to all measured particle removal energies and the ordering of the frontier orbitals. We find that the accuracy of the calculated spectra does not match the expectations based on their level of self-consistency. In particular, for certain starting points G0W0 and scGW0 provide spectra in better agreement with the PES than scGW.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Assessment of correlation energies based on the random-phase approximation

The random-phase approximation to the ground state correlation energy (RPA) in combination with exact exchange (EX) has brought Kohn-Sham (KS) density functional theory one step closer towards a universal, "general purpose first principles method". In an effort to systematically assess the influence of several correlation energy contributions beyond RPA, this work presents dissociation energies of small molecules and solids, activation energies for hydrogen transfer and non-hydrogen transfer reactions, as well as reaction energies for a number of common test sets. We benchmark EX+RPA and several flavors of energy functionals going beyond it: second-order screened exchange (SOSEX), single excitation (SE) corrections, renormalized single excitation (rSE) corrections, as well as their combinations. Both the single excitation correction as well as the SOSEX contribution to the correlation energy significantly improve upon the notorious tendency of EX+RPA to underbind. Surprisingly, activation energies obtained using EX+RPA based on a KS reference alone are remarkably accurate. RPA+SOSEX+rSE provides an equal level of accuracy for reaction as well as activation energies and overall gives the most balanced performance, which makes it applicable to a wide range of systems and chemical reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Density-functional Theory for f electron Systems: the α-γ Phase Transition in Cerium

The isostructural α-γ phase transition in cerium is analyzed using density-functional theory with different exchange-correlation functionals, in particular the PBE0 hybrid functional and the exact- exchange plus correlation in the random-phase approximation [(EX+cRPA)@PBE0] approach. We show that the Hartree-Fock exchange part of the hybrid functional actuates two distinct solutions at zero temperature that can be associated with the α and γ phases of cerium. However, despite the relatively good structural and magnetic properties, PBE0 predicts the γ phase to be the stable phase at ambient pressure and zero temperature, in contradiction with low temperature experiments. EX+cRPA reverses the energetic ordering, which emphasizes the importance of correlation for rare- earth systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions

Efficient implementations of electronic structure methods are essential for first-principles modeling of molecules and solids. We here present a particularly efficient common framework for methods beyond semilocal density-functional theory, including Hartree-Fock (HF), hybrid density functionals, random-phase approximation (RPA), second-order Møller-Plesset perturbation theory (MP2), and the $GW$ method. This computational framework allows us to use compact and accurate numeric atom-centered orbitals (popular in many implementations of semilocal density-functional theory) as basis functions. The essence of our framework is to employ the "resolution of identity (RI)" technique to facilitate the treatment of both the two-electron Coulomb repulsion integrals (required in all these approaches) as well as the linear density-response function (required for RPA and $GW$). This is possible because these quantities can be expressed in terms of products of single-particle basis functions, which can in turn be expanded in a set of auxiliary basis functions (ABFs). The construction of ABFs lies at the heart of the RI technique, and here we propose a simple prescription for constructing the ABFs which can be applied regardless of whether the underlying radial functions have a specific analytical shape (e.g., Gaussian) or are numerically tabulated. We demonstrate the accuracy of our RI implementation for Gaussian and NAO basis functions, as well as the convergence behavior of our NAO basis sets for the above-mentioned methods. Benchmark results are presented for the ionization energies of 50 selected atoms and molecules from the G2 ion test set as obtained with $GW$ and MP2 self-energy methods, and the G2-I atomization energies as well as the S22 molecular interaction energies as obtained with the RPA method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Unified description of ground and excited states of finite systems: the self-consistent GW approach

GW calculations with fully self-consistent G and W -- based on the iterative solution of the Dyson equation -- provide an approach for consistently describing ground and excited states on the same quantum mechanical level. We show that for the systems considered here self-consistent GW reaches the same final Green function regardless of the initial reference state. Self-consistency systematically improves ionization energies and total energies of closed shell systems compared to G_0W_0 based on Hartree-Fock and (semi)local density-functional theory. These improvements also translate to the electron density as exemplified by an improved description of dipole moments and permit us to assess the quality of ground state properties such as bond lengths and vibrational frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Beyond the Random Phase Approximation for the Electron Correlation Energy: The Importance of Single Excitations

The random phase approximation (RPA) for the electron correlation energy, combined with the exact-exchange energy, represents the state-of-the-art exchange-correlation functional within density-functional theory (DFT). However, the standard RPA practice -- evaluating both the exact-exchange and the RPA correlation energy using local or semilocal Kohn-Sham (KS) orbitals -- leads to a systematic underbinding of molecules and solids. Here we demonstrate that this behavior is largely corrected by adding a "single excitation" (SE) contribution, so far not included in the standard RPA scheme. A similar improvement can also be achieved by replacing the non-self-consistent exact-exchange total energy by the corresponding self-consistent Hartree-Fock total energy, while retaining the RPA correlation energy evaluated using Kohn-Sham orbitals. Both schemes achieve chemical accuracy for a standard benchmark set of non-covalent intermolecular interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Large-scale surface reconstruction energetics of Pt(100) and Au(100) by all-electron DFT

The low-index surfaces of Au and Pt all tend to reconstruct, a fact that is of key importance in many nanostructure, catalytic, and electrochemical applications. Remarkably, some significant questions regarding their structural energies remain even today, in particular for the large-scale quasihexagonal reconstructed (100) surfaces: Rather dissimilar reconstruction energies for Au and Pt in available experiments, and experiment and theory do not match for Pt. We here show by all-electron density-functional theory that only large enough "(5 x N)" approximant supercells capture the qualitative reconstruction energy trend between Au(100) and Pt(100), in contrast to what is often done in the theoretical literature. Their magnitudes are then in fact similar, and closer to the measured value for Pt(100); our calculations achieve excellent agreement with known geometric characteristics and provide direct evidence for the electronic reconstruction driving force.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Auger recombination rates in nitrides from first principles

We report Auger recombination rates for wurtzite InGaN calculated from first principles density-functional and many-body-perturbation theory. Two different mechanisms are examined -- inter- and intra-band recombination -- that affect different parts of the emission spectrum. In the blue to green spectral region and at room temperature the Auger coefficient can be as large as 2x10^-30cm^6s^-1; in the infrared even larger. Since Auger recombination scales with the cubic power of the free-carrier concentration it becomes an important non-radiative loss mechanism at high current densities. Our results indicate that Auger recombination may be responsible for the loss of quantum efficiency that affects InGaN-based light emitters.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Vertex corrections in localized and extended systems

Within many-body perturbation theory we apply vertex corrections to various closed-shell atoms and to jellium, using a local approximation for the vertex consistent with starting the many-body perturbation theory from a DFT-LDA Green's function. The vertex appears in two places -- in the screened Coulomb interaction, W, and in the self-energy, Σ-- and we obtain a systematic discrimination of these two effects by turning the vertex in Σon and off. We also make comparisons to standard GW results within the usual random-phase approximation (RPA), which omits the vertex from both. When a vertex is included for closed-shell atoms, both ground-state and excited-state properties demonstrate only limited improvements over standard GW. For jellium we observe marked improvement in the quasiparticle band width when the vertex is included only in W, whereas turning on the vertex in Σleads to an unphysical quasiparticle dispersion and work function. A simple analysis suggests why implementation of the vertex only in W is a valid way to improve quasiparticle energy calculations, while the vertex in Σis unphysical, and points the way to development of improved vertices for ab initio electronic structure calculations.

0 citations0 reviewsNo code linkedOpen access