Roadmap on Electronic Structure Codes in the Exascale Era

Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing.Comment: Submitted as a roadmap article to Modelling and Simulation in Materials Science and Engineering; Address any correspondence to Vikram Gavini ([email protected]) and Danny Perez ([email protected]

Roadmap on Electronic Structure Codes in the Exascale Era

Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing

Roadmap on Electronic Structure Codes in the Exascale Era

Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing

Recent Developments in the General Atomic and Molecular Electronic Structure System

A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods such as the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock, and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as it is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite such as GAMESS is maintained and developed is considered. Future developments are briefly summarized

Large-Scale Atomistic Simulations of Complex and Functional Properties of Ferroic Materials

Ferroelectric (FE) nanostructures have attracted considerable attention as our abilities improve to synthesize them and to predict their properties by theoretical means. Depolarizing field effects at interfaces of FE heterostructures are particularly notable for causing topological defects such as FE vortices and negative dielectric responses in superlattices. In this thesis, I employ two large-scale atomistic techniques, the first-principles-based effective Hamiltonian (HEff) method and the linear-scaling three-dimensional fragment (LS3DF) method. I use these methods to explore optical rotation in FE vortices, electro-optic effects in FE vortices and skyrmions, and voltage amplification via negative capacitance in ferroelectric-paraelectric superlattices. We employ HEff in Monte Carlo and molecular dynamics schemes to maximize spontaneous optical rotation in a BaTiO$_3$/SrTiO$_3$ nanocomposite. For a small bias field, maximal optical rotation was realized at room temperature. The result has acquired greater relevance since Ramesh and coworkers observed ``emergent chirality in FE vortex arrays in PbTiO$_3$/SrTiO$_3$ superlattices. In a similar nanocomposite as above, we use the combined HEff and LS3DF method to study how band gap and band alignment evolves along the path from a polar-toroidal to an electrical skyrmion state. Temperature control of the vortex provides substantially larger range of control of bandgap and band alignment than field control of the skyrmion. Using temperature and electric fields to manipulate polarization and bond angle distortion in both constituent materials provides an additional handle for bandgap engineering in such nanostructures. We then use HEff to study BaTiO$_3$/SrTiO$_3$ superlattices as a platform for negative differential capacitance. We implement an atomistic framework amenable to simulation of negative capacitance in strained superlattices. In these systems, we predict misfit epitaxial strain control allows for broadly extending the operable temperature range for negative. By manipulating this strain, we observed switching of negative capacitance behavior between both constituent materials of the superlattice at low temperature

Complexity Reduction in Density Functional Theory: Locality in Space and Energy

We present recent developments of the NTChem program for performing large scale hybrid Density Functional Theory calculations on the supercomputer Fugaku. We combine these developments with our recently proposed Complexity Reduction Framework to assess the impact of basis set and functional choice on its measures of fragment quality and interaction. We further exploit the all electron representation to study system fragmentation in various energy envelopes. Building off this analysis, we propose two algorithms for computing the orbital energies of the Kohn-Sham Hamiltonian. We demonstrate these algorithms can efficiently be applied to systems composed of thousands of atoms and as an analysis tool that reveals the origin of spectral properties.Comment: Accepted Manuscrip