We describe a complete set of algorithms for ab initio molecular simulations based on numerically tabulated atom-centered orbitals (NAOs) to capture a wide range of molecular and materials properties from quantum-mechanical first principles. The full algorithmic framework described here is embodied in the Fritz Haber Institute "ab initio molecular simulations" (FHI-aims) computer program package. Its comprehensive description should be relevant to any other first-principles implementation based on NAOs. The focus here is on density-functional theory (DFT) in the local and semilocal (generalized gradient) approximations, but an extension to hybrid functionals, Hartree–Fock theory, and MP2/GW electron self-energies for total energies and excited states is possible within the same underlying algorithms. An all-electron/full-potential treatment that is both computationally efficient and accurate is achieved for periodic and cluster geometries on equal footing, including relaxation and ab initio molecular dynamics. We demonstrate the construction of transferable, hierarchical basis sets, allowing the calculation to range from qualitative tight-binding like accuracy to meV-level total energy convergence with the basis set. Since all basis functions are strictly localized, the otherwise computationally dominant grid-based operations scale as O(N) with system size N. Together with a scalar-relativistic treatment, the basis sets provide access to all elements from light to heavy. Both low-communication parallelization of all real-space grid based algorithms and a ScaLapack-based, customized handling of the linear algebra for all matrix operations are possible, guaranteeing efficient scaling (CPU time and memory) up to massively parallel computer systems with thousands of CPUs

No abstract

SummaryX-ray photoelectron spectroscopy (XPS) is a widely used tool for studying the chemical composition of materials and it is a standard technique in surface science and technology. XPS is particularly useful for characterizing nanostructures such as carbon nanomaterials due to their reduced dimensionality. In order to assign the measured binding energies to specific bonding environments, reference energy values need to be known. Experimental measurements of the core level signals of the elements present in novel materials such as graphene have often been compared to values measured for molecules, or calculated for finite clusters. Here we have calculated core level binding energies for variously functionalized or defected graphene by delta Kohn–Sham total energy differences in the real-space grid-based projector-augmented wave density functional theory code (GPAW). To accurately model extended systems, we applied periodic boundary conditions in large unit cells to avoid computational artifacts. In select cases, we compared the results to all-electron calculations using an ab initio molecular simulations (FHI-aims) code. We calculated the carbon and oxygen 1s core level binding energies for oxygen and hydrogen functionalities such as graphane-like hydrogenation, and epoxide, hydroxide and carboxylic functional groups. In all cases, we considered binding energy contributions arising from carbon atoms up to the third nearest neighbor from the functional group, and plotted C 1s line shapes by using experimentally realistic broadenings. Furthermore, we simulated the simplest atomic defects, namely single and double vacancies and the Stone–Thrower–Wales defect. Finally, we studied modifications of a reactive single vacancy with O and H functionalities, and compared the calculated values to data found in the literature.

Published online zzz PACS 71.15.Ap, 71.15.Dx A characteristic feature of the state-of-the-art of real-space methods in electronic structure calculations is the diversity of the techniques used in the discretization of the relevant partial differential equations. In this context, the main approaches include finite-difference methods, various types of finite-elements and wavelets. This paper reports on the results of several code development projects that approach problems related to the electronic structure using these three different discretization methods. We review the ideas behind these methods, give examples of their applications, and discuss their similarities and differences.

We have studied quantum wires using the Green's function technique and the density-functional theory, calculating the electronic structure and the conductance. All the numerics are implemented using the finite-element method with a high-order polynomial basis. For short wires, i.e. quantum point contacts, the zero-bias conductance shows, as a function of the gate voltage and at a finite temperature, a plateau at around 0.7G0. (G0 = 2e 2 /h is the quantum conductance). The behavior, which is caused in our mean-field model by spontaneous spin polarization in the constriction, is reminiscent of the so-called 0.7-anomaly observed in experiments. In our model the temperature and the wire length affect the conductance-gate voltage curves in the same way as in the measured data.

We examine the conductance properties of a chain of Na atoms between two metallic leads in the limit of low bias. Resonant states corresponding to the conductance channel and the local charge neutrality condition cause conductance oscillations as a function of the number of atoms in the chain. Moreover, the geometrical shape of the contact leads influences the conductivity by giving rise to additional oscillations as a function of the lead opening angle.

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, specifically for the large-scale quasihexagonally 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 Χ 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

Lennard-Jones (LJ) parameters are derived for classical nonpolarizable force fields for carbon nanotubes (CNTs) and for CNT-water interaction from van der Waals (vdW) enhanced density functional calculations. The new LJ parameters for carbon-carbon interactions are of the same order as those previously used in the literature but differ significantly for CNT-water interactions. This may partially originate from the fact that in addition to pure vdW interactions the polarization and other quantum mechanics effects are embedded into the LJ-potential.