Search

Modern software engineering of electronic structure codes has seen a paradigm shift from monolithic workflows towards object-based modularity. Software objectivity allows for greater flexibility in the application of electronic structure calculations, with particular benefits when integrated with approaches for data-driven analysis. Here, we discuss different approaches to create "deep" modular interfaces that connect big-data workflows and electronic structure codes, and explore the diversity of use cases that they can enable. We present two such interface approaches for the semi-empirical electronic structure package, DFTB+. In one case, DFTB+ is applied as a library and provides data to an external workflow; and in another, DFTB+ receives data via external bindings and processes the information subsequently within an internal workflow. We provide a general framework to enable data exchange workflows for embedding new machine-learning-based Hamiltonians within DFTB+, or to enabling deep integration of DFTB+ in multiscale embedding workflows. These modular interfaces demonstrate opportunities in emergent software and workflows to accelerate scientific discovery by harnessing existing software capabilities., Comment: 12 pages, 5 figures, Journal of Chemical Physics special issue "modular and interoperable software for chemical physics"

Accurate electronic bandstructures of solids are indispensable for a wide variety of applications and should provide a sound prediction of phonon-induced band gap renormalization at finite temperatures. We employ our previously introduced formalism of general hybrid functionals within the approximate density functional method, DFTB, to present first insights into the accuracy of temperature dependent band gaps obtained by a dielectric-dependent global hybrid functional. The work targets the prototypical group-IV semiconductors diamond and silicon. Following Zacharias et al. [Phys. Rev. Lett. 115, 177401 (2015)], we sample the nuclear wave function by stochastic Monte-Carlo integration as well as the deterministic one-shot procedure [Phys. Rev. B 94, 075125 (2016)] derived from it. The computational efficiency of DFTB enables us to further compare these approaches, which fully take nuclear quantum effects into account, with classical Born-Oppenheimer molecular dynamic (BOMD) simulations. While the quantum mechanical treatments of Zacharias et al. yield band gaps in good agreement with experiment, calculations based on BOMD snapshots inadequately describe the renormalization effect at low temperatures. We demonstrate the importance of properly incorporating nuclear quantum effects by adapting the stochastic approach to normal amplitudes that arise from the classical equipartition principle. For low temperatures the results thus obtained closely resemble the BOMD predictions, while anharmonic effects become important beyond $500\,\mathrm{K}$. Comparisons between DFTB parametrized from semi-local DFT, and global hybrid DFTB, suggest that Fock-type exchange systematically yields a slightly more pronounced electron-phonon interaction, hence stronger gap renormalization and zero-point corrections.

Screened range-separated hybrid (SRSH) functionals within generalized Kohn-Sham density functional theory (GKS-DFT) have been shown to restore a general $1/(r\varepsilon)$ asymptotic decay of the electrostatic interaction in dielectric environments. Major achievements of SRSH include an improved description of optical properties of solids and correct prediction of polarization-induced fundamental gap renormalization in molecular crystals. The density functional tight-binding method (DFTB) is an approximate DFT that bridges the gap between first principles methods and empirical electronic structure schemes. While purely long-range corrected RSH are already accessible within DFTB for molecular systems, this work generalizes the theoretical foundation to also include screened range-separated hybrids, with conventional pure hybrid functionals as a special case. The presented formulation and implementation is also valid for periodic boundary conditions (PBC) beyond the $\Gamma$-point. To treat periodic Fock exchange and its integrable singularity in reciprocal space, we resort to techniques successfully employed by DFT, in particular a truncated Coulomb operator and the minimum image convention. Starting from the first principles Hartree-Fock operator, we derive suitable expressions for the DFTB method, using standard integral approximations and their efficient implementation in the DFTB+ software package. Convergence behavior is investigated and demonstrated for the polyacene series as well as two- and three-dimensional materials. Benzene and pentacene molecular and crystalline systems show the correct polarization-induced gap renormalization by SRSH-DFTB at heavily reduced computational cost compared to first principles methods.

We measure the relaxation time of a square lattice Ising ferromagnet that is quenched to zero-temperature from supercritical initial conditions. We reveal an anomalous and seemingly overlooked timescale associated with the relaxation to "frozen" two-stripe states. While close to a power law of the form $\sim L^{\nu}$ , we argue this timescale actually grows as $\sim L^{2}\ln L$, with L the linear dimension of the system. We uncover the mechanism behind this scaling form by using a synthetic initial condition that replicates the late time ordering of two-stripe states, and subsequently explain it heuristically.

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.

The increasing need to simulate the dynamics of photoexcited molecular and nanosystems in the sub-picosecond regime demands new efficient tools able to describe the quantum nature of matter at a low computational cost. By combining the power of the approximate DFTB method with the semiclassical Ehrenfest method for nuclear-electron dynamics we have achieved a real-time time-dependent DFTB (TD-DFTB) implementation that fits such requirements. In addition to enabling the study of nuclear motion effects in photoinduced charge transfer processes, our code adds novel features to the realm of static and time-resolved computational spectroscopies. In particular, the optical properties of periodic materials such as graphene nanoribbons or the use of corrections such as the "LDA+U" and "pseudo SIC" methods to improve the optical properties in some systems, can now be handled at the TD-DFTB level. Moreover, the simulation of fully-atomistic time-resolved transient absorption spectra and impulsive vibrational spectra can now be achieved within reasonable computing time, owing to the good performance of the implementation and a parallel simulation protocol. Its application to the study of UV/visible light-induced vibrational coherences in molecules is demonstrated and opens a new door into the mechanisms of non-equilibrium ultrafast phenomena in countless materials with relevant applications., Comment: 51 pages, 7 figures, 1 TOC graphic, submitted to JCTC

We study the time evolution of the two-dimensional kinetic Ising model in finite systems with a non-conserved order parameter, considering nearest-neighbour interactions on the square lattice with periodic and open boundary conditions. Universal data collapse in spin product correlation functions is observed which, when expressed in rescaled units, is valid across the entire time evolution of the system at all length scales, not just within the time regime usually considered in the dynamical scaling hypothesis. Consequently, beyond rapidly decaying finite size effects, the evolution of correlations in small finite systems parallels arbitrarily larger cases, even at large fractions of the size of these finite systems., Comment: 6 pages 6 figures

The nature and interplay of intrinsic point and extended defects in n-type Si-doped InN epilayers with free carrier concentrations up to 6.6x10E20cm-3 are studied using positron annihilation spectroscopy and transmission electron microscopy and compared to results from undoped irradiated films. In as-grown Si-doped samples, V_In-V_N complexes are the dominant III-sublattice related vacancy defects. Enhanced formation of larger V_In-mV_N clusters is observed at the interface, which speaks for high concentrations of additional V_N in the near-interface region and coincides with an increase in the density of screw and edge type dislocations in that area., Comment: 4 pages, 3 figures