Stark Ionization of Atoms and Molecules within Density Functional Resonance Theory

We show that the energetics and lifetimes of resonances of finite systems under an external electric field can be captured by Kohn--Sham density functional theory (DFT) within the formalism of uniform complex scaling. Properties of resonances are calculated self-consistently in terms of complex densities, potentials and wavefunctions using adapted versions of the known algorithms from DFT. We illustrate this new formalism by calculating ionization rates using the complex-scaled local density approximation and exact exchange. We consider a variety of atoms (H, He, Li and Be) as well as the hydrogen molecule. Extensions are briefly discussed.Comment: 5 pages, 5 figures. This document is the unedited Author's version of a Submitted Work that was subsequently accepted for publication in J.Phys.Chem.Lett., copyright (c) American Chemical Society after peer review. To access the final edited and published work see http://pubs.acs.org/doi/abs/10.1021/jz401110

Localized atomic basis set in the projector augmented wave method

We present an implementation of localized atomic orbital basis sets in the projector augmented wave (PAW) formalism within the density functional theory (DFT). The implementation in the real-space GPAW code provides a complementary basis set to the accurate but computationally more demanding grid representation. The possibility to switch seamlessly between the two representations implies that simulations employing the local basis can be fine tuned at the end of the calculation by switching to the grid, thereby combining the strength of the two representations for optimal performance. The implementation is tested by calculating atomization energies and equilibrium bulk properties of a variety of molecules and solids, comparing to the grid results. Finally, it is demonstrated how a grid-quality structure optimization can be performed with significantly reduced computational effort by switching between the grid and basis representations.Comment: 10 pages, 5 figures. http://prb.aps.org.globalproxy.cvt.dk/abstract/PRB/v80/i19/e19511

Optical Excitations of Chlorophyll a and b Monomers and Dimers

11 pages, 9 figuresPreprin

Recent Progress of the Computational 2D Materials Database (C2DB)

The C2DB is a highly curated open database organizing a wealth of computed properties for more than 4000 atomically thin two-dimensional (2D) materials. Here we report on new materials and properties that were added to the database since its first release in 2018. The set of new materials comprise several hundred monolayers exfoliated from experimentally known layered bulk materials, (homo)bilayers in various stacking configurations, native point defects in semiconducting monolayers, and chalcogen/halogen Janus monolayers. The new properties include exfoliation energies, Bader charges, spontaneous polarisations, Born charges, infrared polarisabilities, piezoelectric tensors, band topology invariants, exchange couplings, Raman- and second harmonic generation spectra. We also describe refinements of the employed material classification schemes, upgrades of the computational methodologies used for property evaluations, as well as significant enhancements of the data documentation and provenance. Finally, we explore the performance of Gaussian process-based regression for efficient prediction of mechanical and electronic materials properties. The combination of open access, detailed documentation, and extremely rich materials property data sets make the C2DB a unique resource that will advance the science of atomically thin materials.Comment: 30 pages, 26 figure

GPAW: open Python package for electronic-structure calculations

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

Roadmap on methods and software for electronic structure based simulations in chemistry and materials

Abstract: This Roadmap article provides a succinct, comprehensive overview of the state of electronic structure methods and software for molecular and materials simulations. Seventeen distinct sections collect insights by 51 leading scientists in the field. Each contribution addresses the status of a particular area, as well as current challenges and anticipated future advances, with a particular eye towards software related aspects and providing key references for further reading. Foundational sections cover density functional theory and its implementation in real-world simulation frameworks, Green's function based many-body perturbation theory, wave-function based and stochastic electronic structure approaches, relativistic effects and semiempirical electronic structure theory approaches. Subsequent sections cover nuclear quantum effects, real-time propagation of the electronic structure, challenges for computational spectroscopy simulations, and exploration of complex potential energy surfaces. The final sections summarize practical aspects, including computational workflows for complex simulation tasks, the impact of current and future high-performance computing architectures, software engineering practices, education and training to maintain and broaden the community, as well as the status of and needs for electronic structure based modeling from the vantage point of industry environments. Overall, the field of electronic structure software and method development continues to unlock immense opportunities for future scientific discovery, based on the growing ability of computations to reveal complex phenomena, processes and properties that are determined by the make-up of matter at the atomic scale, with high precision

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light–matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).This work was supported by the European Research Council (Grant No. ERC-2015-AdG694097), the Cluster of Excellence “Advanced Imaging of Matter” (AIM), Grupos Consolidados (IT1249-19), and SFB925. The Flatiron Institute is a division of the Simons Foundation. X.A., A.W., and A.C. acknowledge that part of this work was performed under the auspices of the U.S. Department of Energy at Lawrence Livermore National Laboratory under Contract No. DE-AC52-07A27344. J.J.-S. gratefully acknowledges the funding from the European Union Horizon 2020 Research and Innovation Program under the Marie Sklodowska-Curie Grant Agreement No. 795246-StrongLights. J.F. acknowledges financial support from the Deutsche Forschungsgemeinschaft (DFG Forschungsstipendium FL 997/1-1). D.A.S. acknowledges University of California, Merced start-up funding.Peer reviewe

Density functional theory based screening of ternary alkali-transition metal borohydrides: A computational material design project

The dissociation of molecules, even the most simple hydrogen molecule, cannot be described accurately within density functional theory because none of the currently available functionals accounts for strong on-site correlation. This problem led to a discussion of properties that the local Kohn-Sham potential has to satisfy in order to correctly describe strongly correlated systems. We derive an analytic expression for the nontrivial form of the Kohn-Sham potential in between the two fragments for the dissociation of a single bond. We show that the numerical calculations for a one-dimensional two-electron model system indeed approach and reach this limit. It is shown that the functional form of the potential is universal, i.e., independent of the details of the two fragments.We acknowledge funding by the Spanish MEC (Grant No. FIS2007-65702-C02-01), “Grupos Consolidados UPV/EHU del Gobierno Vasco” (Grant No. IT-319-07), and the European Community through e-I3 ETSF project (Grant Agreement No. 211956).Peer reviewe