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

Volker Blum

Volker Blum 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.comp-phphysics.chem-phBiomoleculesphysics.atm-clusBiological Physicscond-mat.other

Catalog footprint

What is connected

25works
7topics
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

25 published item(s)

preprint2022arXiv

Nuclear-Electronic Orbital Approach to Quantization of Protons in Periodic Electronic Structure Calculations

The nuclear-electronic orbital (NEO) method is a well-established approach for treating nuclei quantum mechanically in molecular systems beyond the usual Born-Oppenheimer approximation. In this work, we present a strategy to implement the NEO method for periodic electronic structure calculations, particularly focused on multicomponent density functional theory (DFT). The NEO-DFT method is implemented in an all-electron electronic structure code, FHI-aims, using a combination of analytical and numerical integration techniques as well as a resolution of the identity scheme to enhance computational efficiency. After validating this implementation, proof-of-concept applications are presented to illustrate the effects of quantized protons on the physical properties of extended systems such as two-dimensional materials and liquid-semiconductor interfaces. Specifically, periodic NEO-DFT calculations are performed for a trans-polyacetylene chain, a hydrogen boride sheet, and a titanium oxide-water interface. The zero-point energy effects of the protons, as well as electron-proton correlation, are shown to noticeably impact the density of states and band structures for these systems. These developments provide a foundation for the application of multicomponent DFT to a wide range of other extended condensed matter systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

All-electron periodic $G_0W_0$ implementation with numerical atomic orbital basis functions: algorithm and benchmarks

We present an all-electron, periodic {\GnWn} implementation within the numerical atomic orbital (NAO) basis framework. A localized variant of the resolution-of-the-identity (RI) approximation is employed to significantly reduce the computational cost of evaluating and storing the two-electron Coulomb repulsion integrals. We demonstrate that the error arising from localized RI approximation can be reduced to an insignificant level by enhancing the set of auxiliary basis functions, used to expand the products of two single-particle NAOs. An efficient algorithm is introduced to deal with the Coulomb singularity in the Brillouin zone sampling that is suitable for the NAO framework. We perform systematic convergence tests and identify a set of computational parameters, which can serve as the default choice for most practical purposes. Benchmark calculations are carried out for a set of prototypical semiconductors and insulators, and compared to independent reference values obtained from an independent $G_0W_0$ implementation based on linearized augmented plane waves (LAPW) plus high-energy localized orbitals (HLOs) basis set, as well as experimental results. With a moderate (FHI-aims \textit{tier} 2) NAO basis set, our $G_0W_0$ calculations produce band gaps that typically lie in between the standard LAPW and the LAPW+HLO results. Complementing \textit{tier} 2 with highly localized Slater-type orbitals (STOs), we find that the obtained band gaps show an overall convergence towards the LAPW+HLO results. The algorithms and techniques developed in this work pave the way for efficient implementations of correlated methods within the NAO framework.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

GPU-Acceleration of the ELPA2 Distributed Eigensolver for Dense Symmetric and Hermitian Eigenproblems

The solution of eigenproblems is often a key computational bottleneck that limits the tractable system size of numerical algorithms, among them electronic structure theory in chemistry and in condensed matter physics. Large eigenproblems can easily exceed the capacity of a single compute node, thus must be solved on distributed-memory parallel computers. We here present GPU-oriented optimizations of the ELPA two-stage tridiagonalization eigensolver (ELPA2). On top of cuBLAS-based GPU offloading, we add a CUDA kernel to speed up the back-transformation of eigenvectors, which can be the computationally most expensive part of the two-stage tridiagonalization algorithm. We benchmark the performance of this GPU-accelerated eigensolver on two hybrid CPU-GPU architectures, namely a compute cluster based on Intel Xeon Gold CPUs and NVIDIA Volta GPUs, and the Summit supercomputer based on IBM POWER9 CPUs and NVIDIA Volta GPUs. Consistent with previous benchmarks on CPU-only architectures, the GPU-accelerated two-stage solver exhibits a parallel performance superior to the one-stage counterpart. Finally, we demonstrate the performance of the GPU-accelerated eigensolver developed in this work for routine semi-local KS-DFT calculations comprising thousands of atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

ELSI -- An Open Infrastructure for Electronic Structure Solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn-Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. Finally, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pentacene and Tetracene Molecules and Films on H/Si(111): Level Alignment from Hybrid Density Functional Theory

The electronic properties of hybrid organic-inorganic semiconductor interfaces depend strongly on the alignment of the electronic carrier levels in the organic/inorganic components. In the present work, we address this energy level alignment from first principles theory for two paradigmatic organic-inorganic semiconductor interfaces, the singlet fission materials tetracene and pentacene on H/Si(111), using all-electron hybrid density functional theory. For isolated tetracene on H/Si(111), a type I-like heterojunction (lowest-energy electron and hole states on Si) is found. For isolated pentacene, the molecular and semiconductor valence band edges are degenerate. For monolayer films, we show how to construct supercell geometries with up to 1,192 atoms, which minimize the strain between the inorganic surface and an organic monolayer film. Based on these models, we predict the formation of type II heterojunctions (electron states on Si, hole-like states on the organic species) for both acenes, indicating that charge separation at the interface between the organic and inorganic components is favored. The paper discusses the steps needed to find appropriate low-energy interface geometries for weakly bonded organic molecules and films on inorganic substrates from first principles, a necessary prerequisite for any computational level alignment prediction.

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
preprint2020arXiv

SIESTA: recent developments and applications

A review of the present status, recent enhancements, and applicability of the SIESTA program is presented. Since its debut in the mid-nineties, SIESTA's flexibility, efficiency and free distribution has given advanced materials simulation capabilities to many groups worldwide. The core methodological scheme of SIESTA combines finite-support pseudo-atomic orbitals as basis sets, norm-conserving pseudopotentials, and a real-space grid for the representation of charge density and potentials and the computation of their associated matrix elements. Here we describe the more recent implementations on top of that core scheme, which include: full spin-orbit interaction, non-repeated and multiple-contact ballistic electron transport, DFT+U and hybrid functionals, time-dependent DFT, novel reduced-scaling solvers, density-functional perturbation theory, efficient Van der Waals non-local density functionals, and enhanced molecular-dynamics options. In addition, a substantial effort has been made in enhancing interoperability and interfacing with other codes and utilities, such as Wannier90 and the second-principles modelling it can be used for, an AiiDA plugin for workflow automatization, interface to Lua for steering SIESTA runs, and various postprocessing utilities. SIESTA has also been engaged in the Electronic Structure Library effort from its inception, which has allowed the sharing of various low level libraries, as well as data standards and support for them, in particular the PSML definition and library for transferable pseudopotentials, and the interface to the ELSI library of solvers. Code sharing is made easier by the new open-source licensing model of the program. This review also presents examples of application of the capabilities of the code, as well as a view of on-going and future developments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The CECAM Electronic Structure Library and the modular software development paradigm

First-principles electronic structure calculations are very widely used thanks to the many successful software packages available. Their traditional coding paradigm is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, from the compiler up, with the exception of linear-algebra and message-passing libraries. This model has been quite successful for decades. The rapid progress in methodology, however, has resulted in an ever increasing complexity of those programs, which implies a growing amount of replication in coding and in the recurrent re-engineering needed to adapt to evolving hardware architecture. The Electronic Structure Library (\esl) was initiated by CECAM (European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure programs and redesign them as free, open-source libraries. They include "heavy-duty" ones with a high degree of parallelisation, and potential for adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by scientists when implementing new ideas. It is a community effort, undertaken by developers of various successful codes, now facing the challenges arising in the new model. This modular paradigm will improve overall coding efficiency and enable specialists (computer scientists or computational scientists) to use their skills more effectively. It will lead to a more sustainable and dynamic evolution of software as well as lower barriers to entry for new developers.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Ab Initio Bethe-Salpeter Equation Approach to Neutral Excitations in Molecules with Numeric Atom-Centered Orbitals

The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as "Thiel's set". Literature reference data based on Gaussian-type orbitals are reproduced to about one meV precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted "tier2+aug2" is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

GPGPU Acceleration of All-Electron Electronic Structure Theory Using Localized Numeric Atom-Centered Basis Functions

We present an implementation of all-electron density-functional theory for massively parallel GPGPU-based platforms, using localized atom-centered basis functions and real-space integration grids. Special attention is paid to domain decomposition of the problem on non-uniform grids, which enables compute- and memory-parallel execution across thousands of nodes for real-space operations, e.g. the update of the electron density, the integration of the real-space Hamiltonian matrix, and calculation of Pulay forces. To assess the performance of our GPGPU implementation, we performed benchmarks on three different architectures using a 103-material test set. We find that operations which rely on dense serial linear algebra show dramatic speedups from GPGPU acceleration: in particular, SCF iterations including force and stress calculations exhibit speedups ranging from 4.5 to 6.6. For the architectures and problem types investigated here, this translates to an expected overall speedup between 3-4 for the entire calculation (including non-GPU accelerated parts), for problems featuring several tens to hundreds of atoms. Additional calculations for a 375-atom Bi$_2$Se$_3$ bilayer show that the present GPGPU strategy scales for large-scale distributed-parallel simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

MatD3: A Database and Online Presentation Package for Research Data Supporting Materials Discovery, Design, and Dissemination

The discovery of new materials as well as the determination of a vast set of materials properties for science and technology is a fast growing field of research, with contributions from many groups worldwide. Materials data from individual research groups is traditionally disseminated by means of loosely interconnected, peer-reviewed publications. MatD3 is an open-source, dedicated database and web application framework designed to store, curate and disseminate experimental and theoretical materials data generated by individual research groups or research consortia. A research group can set up its own instance of MatD3 and publish scientific results or simply use an existing online MatD3 instance. Disseminating research data in this form enables broader access, reproducibility, and repurposing of scientific products. MatD3 is a general purpose database that does not focus on any specific level of theory or experimental method. Instead, the focus is on storing and making accessible the data and making it straightforward to curate them.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Molecular NMR shieldings, J-couplings, and magnetizabilities from numeric atom-centered orbital based density-functional calculations

We describe an accurate and scalable implementation for the computation of molecular nuclear magnetic resonance shieldings, J-couplings, and magnetizabilities within nonrelativistic semilocal density functional theory, based on numeric atom-centered orbital (NAO) basis sets. We compare the convergence to the basis set limit for two established types of NAO basis sets, called NAO-VCC-nZ and FHI-aims-09, to several established Gaussian-type basis sets. The basis set limit is reached faster for the NAO basis sets than for standard correlation consistent Gaussian-type basis sets (cc-pVnZ, aug-cc-pVnZ, cc-pCVnZ, aug-cc-pCVnZ). For shieldings, the convergence properties and accuracy of the NAO-VCC-nZ basis sets are similar to Jensen's polarization consistent (pc) basis sets optimized for shieldings (pcS-n). For J-couplings, we develop a new type of NAO basis set (NAO-J-n) by augmenting the NAO-VCC-nZ basis sets with tight s-functions from Jensen's pcJ-n basis sets, which are optimized for J-couplings. We find the convergence of the NAO-J-n to be similar to the pcJ-n basis sets. Large scale applicability of the implementation is demonstrated for shieldings and J-couplings in a system of over 1,000 atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

First-principles data set of 45,892 isolated and cation-coordinated conformers of 20 proteinogenic amino acids

We present a structural data set of the 20 proteinogenic amino acids and their amino-methylated and acetylated (capped) dipeptides. Different protonation states of the backbone (uncharged and zwitterionic) were considered for the amino acids as well as varied side chain protonation states. Furthermore, we studied amino acids and dipeptides in complex with divalent cations (Ca2+, Ba2+, Sr2+, Cd2+, Pb2+, and Hg2+). The database covers the conformational hierarchies of 280 systems in a wide relative energy range of up to 4 eV (390 kJ/mol), summing up to an overall of 45,892 stationary points on the respective potential-energy surfaces. All systems were calculated on equal first-principles footing, applying density-functional theory in the generalized gradient approximation corrected for long-range van der Waals interactions. We show good agreement to available experimental data for gas-phase ion affinities. Our curated data can be utilized, for example, for a wide comparison across chemical space of the building blocks of life, for the parametrization of protein force fields, and for the calculation of reference spectra for biophysical applications.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Towards rational design of carbon nitride photocatalysts: Identification of cyanamide "defects" as catalytically relevant sites

The heptazine-based polymer melon (also known as graphitic carbon nitride, g-C3N4), is a promising photocatalyst for hydrogen evolution. Nonetheless, attempts to improve its inherently low activity are rarely based on rational approaches due to a lack of fundamental understanding of its mechanistic operation. Here, we employ molecular heptazine-based model catalysts to identify the cyanamide moiety as a photocatalytically relevant "defect". We exploit this knowledge for the rational design of a carbon nitride polymer populated with cyanamide groups, yielding a material with 12- and 16-times the hydrogen evolution rate and apparent quantum efficiency (400 nm), respectively, compared to the benchmark melon. Computational modelling and material characterization suggest this moiety improves co-ordination (and, in turn, charge transfer kinetics) to the platinum co-catalyst and enhances the separation of the photo-generated charge carriers. The demonstrated knowledge transfer for rational catalyst design presented here provides the conceptual framework for engineering high performance heptazine-based photocatalysts.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Trends for isolated amino acids and dipeptides: Conformation, divalent ion binding, and remarkable similarity of binding to calcium and lead

We derive structural and binding energy trends for twenty amino acids, their dipeptides, and their interactions with the divalent cations Ca$^{2+}$, Ba$^{2+}$, Sr$^{2+}$, Cd$^{2+}$, Pb$^{2+}$, and Hg$^{2+}$. The underlying data set consists of 45,892 first-principles predicted conformers with relative energies up to about 4 eV (about 400kJ/mol). We show that only very few distinct backbone structures of isolated amino acids and their dipeptides emerge as lowest-energy conformers. The isolated amino acids predominantly adopt structures that involve an acidic proton shared between the carboxy and amino function. Dipeptides adopt one of two intramolecular-hydrogen bonded conformations C$_5$ or equatorial C$_7$. Upon complexation with a divalent cation, the accessible conformational space shrinks and intramolecular hydrogen bonding is prevented due to strong electrostatic interaction of backbone and side chain functional groups with cations. Clear correlations emerge from the binding energies of the six divalent ions with amino acids and dipeptides. Cd$^{2+}$ and Hg$^{2+}$ show the largest binding energies - a potential correlation with their known high acute toxicities. Ca$^{2+}$ and Pb$^{2+}$ reveal almost identical binding energies across the entire series of amino acids and dipeptides. This observation validates past indications that ion-mimicry of calcium and lead should play an important role in a toxicological context.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

First-principles molecular structure search with a genetic algorithm

The identification of low-energy conformers for a given molecule is a fundamental problem in computational chemistry and cheminformatics. We assess here a conformer search that employs a genetic algorithm for sampling the low-energy segment of the conformation space of molecules. The algorithm is designed to work with first-principles methods, facilitated by the incorporation of local optimization and blacklisting conformers to prevent repeated evaluations of very similar solutions. The aim of the search is not only to find the global minimum, but to predict all conformers within an energy window above the global minimum. The performance of the search strategy is: (i) evaluated for a reference data set extracted from a database with amino acid dipeptide conformers obtained by an extensive combined force field and first-principles search and (ii) compared to the performance of a systematic search and a random conformer generator for the example of a drug-like ligand with 43 atoms, 8 rotatable bonds and 1 cis/trans bond.

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

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

Embedded-Cluster Calculations in a Numeric Atomic Orbital Density-Functional Theory Framework

We integrate the all-electron electronic structure code FHI-aims into the general ChemShell package for solid-state embedding (QM/MM) calculations. A major undertaking in this integration is the implementation of pseudopotential functionality into FHI-aims to describe cations at the QM/MM boundary through effective core potentials and therewith prevent spurious overpolarization of the electronic density. Based on numeric atomic orbital basis sets, FHI-aims offers particularly efficient access to exact exchange and second order perturbation theory, rendering the established QM/MM setup an ideal tool for hybrid and double-hybrid level DFT calculations of solid systems. We illustrate this capability by calculating the reduction potential of Fe in the Fe-substituted ZSM-5 zeolitic framework and the reaction energy profile for (photo-)catalytic water oxidation at TiO2(110).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How Cations Change Peptide Structure

Specific interactions between cations and proteins have a strong impact on peptide and protein structure. We here shed light on the nature of the underlying interactions, especially regarding the effects on the polyamide backbone structure. To do so, we compare the conformational ensembles of model peptides in isolation and in the presence of either Li+ or Na+ cations by state-of-the-art density-functional theory (including van der Waals effects) and gas-phase infrared spectroscopy. These monovalent cations have a drastic effect on the local backbone conformation of turn-forming peptides, by disruption of the H bonding networks and the resulting severe distortion of the backbone conformations. In fact, Li+ and Na+ can even have different conformational effects on the same peptide. We also assess the predictive power of current approximate density functionals for peptide-cation systems and compare to results from established protein force fields as well as to high-level quantum chemistry (CCSD(T)).

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

Impact of vibrational entropy on the stability of unsolvated peptide helices with increasing length

Helices are a key folding motif in protein structure. The question which factors determine helix stability for a given polypeptide or protein is an ongoing challenge. Here we use van der Waals corrected density-functional theory to address a part of this question in a bottom-up approach. We show how intrinsic helical structure is stabilized with length and temperature for a series of experimentally well studied unsolvated alanine based polypeptides, Ac-Alan-LysH+. By exploring extensively the conformational space of these molecules, we find that helices emerge as the preferred structure in the length range n=4-8 not just due to enthalpic factors (hydrogen bonds and their cooperativity, van der Waals dispersion interactions, electrostatics), but importantly also by a vibrational entropic stabilization over competing conformers at room temperature. The stabilization is shown to be due to softer low-frequency vibrational modes in helical conformers than in more compact ones. This observation is corroborated by including anharmonic effects explicitly through \emph{ab initio} molecular dynamics, and generalized by testing different terminations and considering larger helical peptide models.

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
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