Search

We developed a general framework for hybrid quantum-classical computing of molecular and periodic embedding approaches based on an orbital space separation of the fragment and environment degrees of freedom. We demonstrate its potential by presenting a specific implementation of periodic range-separated DFT coupled to a quantum circuit ansatz, whereby the variational quantum eigensolver and the quantum equation-of-motion algorithm are used to obtain the low-lying spectrum of the embedded fragment Hamiltonian. Application of this scheme to study localized electronic states in materials is showcased through the accurate prediction of the optical properties of the neutral oxygen vacancy in magnesium oxide (MgO). Despite some discrepancies in the position of the main absorption band, the method demonstrates competitive performance compared to state-of-the-art ab initio approaches, particularly evidenced by the excellent agreement with the experimental photoluminescence emission peak.

Water absorption and dissociation processes on pristine low-index TiO$_2$ interfaces are important but poorly understood outside the well-studied anatase (101) and rutile (110). To understand these, we construct three sets of machine learning potentials that are simultaneously applicable to various TiO$_2$ surfaces, based on three density-functional-theory approximations. We compute the water dissociation free energies on seven pristine TiO$_2$ surfaces, and predict that anatase (100), anatase (110), rutile (001), and rutile (011) surfaces favor water dissociation, anatase (101) and rutile (100) surfaces have mostly molecular absorption, while the calculated dissociation fraction on rutile (110) surface may depend on the density functional assumed. Moreover, using an automated algorithm, we reveal that these surfaces follow different types of atomistic mechanisms for proton transfer and water dissociation: one-step, two-step, or both. Our finding thus demonstrates that the different pristine TiO$_2$ surfaces react with water in distinct ways, and cannot be represented using just the low-energy anatase (101) and rutile (110) surfaces.

Linear-response time-dependent density functional theory (LR-TDDFT) for core level spectroscopy using standard local functionals suffers from self-interaction error and a lack of orbital relaxation upon creation of the core hole. As a result, LR-TDDFT calculated X-ray absorption near edge structure (XANES) spectra need to be shifted along the energy axis to match experimental data. We propose a correction scheme based on many body perturbation theory to calculate the shift from first principles. The ionization potential of the core donor state is first computed and then substituted for the corresponding Kohn--Sham orbital energy, thus emulating Koopmans' condition. Both self-interaction error and orbital relaxation are taken into account. The method exploits the localized nature of core states for efficiency and integrates seamlessly in our previous implementation of core level LR-TDDFT, yielding corrected spectra in a single calculation. We benchmark the correction scheme on molecules at the K- and L-edges as well as for core binding energies and report accuracies comparable to higher order methods. We also demonstrate applicability in large and extended systems and discuss efficient approximations., Comment: The following article has been submitted to The Journal of Chemical Physics

The GW approximation of many-body perturbation theory is an accurate method for computing electron addition and removal energies of molecules and solids. In a canonical implementation, however, its computational cost is $O(N^4)$ in the system size N, which prohibits its application to many systems of interest. We present a full-frequency GW algorithm in a Gaussian type basis, whose computational cost scales with $N^2$ to $N^3$. The implementation is optimized for massively parallel execution on state-of-the-art supercomputers and is suitable for nanostructures and molecules in the gas, liquid or condensed phase, using either pseudopotentials or all electrons. We validate the accuracy of the algorithm on the GW100 molecular test set, finding mean absolute deviations of 35 meV for ionization potentials and 27 meV for electron affinities. Furthermore, we study the length-dependence of quasiparticle energies in armchair graphene nanoribbons of up to 1734 atoms in size, and compute the local density of states across a nanoscale heterojunction.

We push the boundaries of electronic structure-based \textit{ab-initio} molecular dynamics (AIMD) beyond 100 million atoms. This scale is otherwise barely reachable with classical force-field methods or novel neural network and machine learning potentials. We achieve this breakthrough by combining innovations in linear-scaling AIMD, efficient and approximate sparse linear algebra, low and mixed-precision floating-point computation on GPUs, and a compensation scheme for the errors introduced by numerical approximations. The core of our work is the non-orthogonalized local submatrix method (NOLSM), which scales very favorably to massively parallel computing systems and translates large sparse matrix operations into highly parallel, dense matrix operations that are ideally suited to hardware accelerators. We demonstrate that the NOLSM method, which is at the center point of each AIMD step, is able to achieve a sustained performance of 324 PFLOP/s in mixed FP16/FP32 precision corresponding to an efficiency of 67.7% when running on 1536 NVIDIA A100 GPUs., Comment: 12 pages, 11 figures

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2k to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension., Comment: 51 pages, 5 figures

Advanced algorithms for large-scale electronic structure calculations are mostly based on processing multi-dimensional sparse data. Examples are sparse matrix-matrix multiplications in linear-scaling Kohn-Sham calculations or the efficient determination of the exact exchange energy. When going beyond mean field approaches, e.g. for Moller-Plesset perturbation theory, RPA and Coupled-Cluster methods, or the GW methods, it becomes necessary to manipulate higher-order sparse tensors. Very similar problems are also encountered in other domains, like signal processing, data mining, computer vision, and machine learning. With the idea that the most of the tensor operations can be mapped to matrices, we have implemented sparse tensor algebra functionalities in the frames of the sparse matrix linear algebra library DBCSR (Distributed Block Compressed Sparse Row). DBCSR has been specifically designed to efficiently perform blocked-sparse matrix operations, so it becomes natural to extend its functionality to include tensor operations. We describe the newly developed tensor interface and algorithms. In particular, we introduce the tensor contraction based on a fast rectangular sparse matrix multiplication algorithm., Comment: 10 pages, 3 figures, submitted at the ParCO2019 conference, Prague, Czech Republic, 10-13 September 2019

Most, if not all the modern scientific simulation packages utilize matrix algebra operations. Among the operation of the linear algebra, one of the most important kernels is the multiplication of matrices, dense and sparse. Examples of application of such a kernel are in electronic structure calculations, machine learning, data mining, graph processing, and digital signal processing. Several optimized libraries exist that can achieve high-performance on distributed systems. Only a few of them target distributed GPU-accelerated systems. In most of the cases, these libraries are provided and optimized by system vendors for their specific computer systems. In this paper, we present the DBCSR library (Distributed Block Compressed Sparse Row) for the distributed dense matrix-matrix multiplications. Although the library is specifically designed for block-sparse matrix-matrix multiplications, we optimized it for the dense case on GPU-accelerated systems. We show that the DBCSR outperforms the multiplication of matrices of different sizes and shapes provided by a vendor optimized GPU version of the ScaLAPACK library up to 2.5x (1.4x on average)., Comment: 5 pages, 4 figures, proceeding of the 2019 International Multi-Conference on Engineering, Computer and Information Sciences (SIBIRCON)

Surfaces of metal oxides at working conditions are usually electrified due to the acid-base chemistry. The charged interface compensated with counterions forms the so-called electric double layer. The coupling of surface chemistry and electric double layer is considered to be crucial but poorly understood because of lacking the information at the atomistic scale. Here, we used the latest development in density functional theory based finite-field molecular dynamics simulation to investigate pH-dependence of the Helmholtz capacitance at electrified rutile TiO$_2$ (110)-NaCl electrolyte interfaces. It is found that, due to competing forces from surface adsorption and from electric double layer, water molecules have a stronger structural fluctuation at high pH and this leads to a much larger capacitance. It is also seen that, interfacial proton transfers at low pH increase significantly the capacitance value. These findings elucidate the microscopic origin for the same trend observed in titration experiments.

Multiplication of two sparse matrices is a key operation in the simulation of the electronic structure of systems containing thousands of atoms and electrons. The highly optimized sparse linear algebra library DBCSR (Distributed Block Compressed Sparse Row) has been specifically designed to efficiently perform such sparse matrix-matrix multiplications. This library is the basic building block for linear scaling electronic structure theory and low scaling correlated methods in CP2K. It is parallelized using MPI and OpenMP, and can exploit GPU accelerators by means of CUDA. We describe a performance comparison of DBCSR on systems with Intel Xeon Phi Knights Landing (KNL) processors, with respect to systems with Intel Xeon CPUs (including systems with GPUs). We find that the DBCSR on Cray XC40 KNL-based systems is 11%-14% slower than on a hybrid Cray XC50 with Nvidia P100 cards, at the same number of nodes. When compared to a Cray XC40 system equipped with dual-socket Intel Xeon CPUs, the KNL is up to 24% faster., Comment: Submitted to the ParCo2017 conference, Bologna, Italy 12-15 September 2017

Understanding lithium solvation and transport in ionic liquids is important due to their possible application in electrochemical devices. Using first-principles simulations aided by a metadynamics approach we study the free-energy landscape for lithium ions at infinite dilution in ethylammonium nitrate, a protic ionic liquid. We analyze the local structure of the liquid around the lithium cation and obtain a quantitative picture in agreement with experimental findings. Our simulations show that the lowest two free energy minima correspond to conformations with the lithium ion being solvated either by three or four nitrate ions with a transition barrier between them of 0.2 \eV. Other less probable conformations having different solvation pattern are also investigated.

Matrix-matrix multiplication is a basic operation in linear algebra and an essential building block for a wide range of algorithms in various scientific fields. Theory and implementation for the dense, square matrix case are well-developed. If matrices are sparse, with application-specific sparsity patterns, the optimal implementation remains an open question. Here, we explore the performance of communication reducing 2.5D algorithms and one-sided MPI communication in the context of linear scaling electronic structure theory. In particular, we extend the DBCSR sparse matrix library, which is the basic building block for linear scaling electronic structure theory and low scaling correlated methods in CP2K. The library is specifically designed to efficiently perform block-sparse matrix-matrix multiplication of matrices with a relatively large occupation. Here, we compare the performance of the original implementation based on Cannon's algorithm and MPI point-to-point communication, with an implementation based on MPI one-sided communications (RMA), in both a 2D and a 2.5D approach. The 2.5D approach trades memory and auxiliary operations for reduced communication, which can lead to a speedup if communication is dominant. The 2.5D algorithm is somewhat easier to implement with one-sided communications. A detailed description of the implementation is provided, also for non ideal processor topologies, since this is important for actual applications. Given the importance of the precise sparsity pattern, and even the actual matrix data, which decides the effective fill-in upon multiplication, the tests are performed within the CP2K package with application benchmarks. Results show a substantial boost in performance for the RMA based 2.5D algorithm, up to 1.80x, which is observed to increase with the number of involved processes in the parallelization., Comment: In Proceedings of PASC '17, Lugano, Switzerland, June 26-28, 2017, 10 pages, 4 figures

An integral scheme for the efficient evaluation of two-center integrals over contracted solid harmonic Gaussian functions is presented. Integral expressions are derived for local operators that depend on the position vector of one of the two Gaussian centers. These expressions are then used to derive the formula for three-index overlap integrals where two of the three Gaussians are located at the same center. The efficient evaluation of the latter is essential for local resolution-of-the-identity techniques that employ an overlap metric. We compare the performance of our integral scheme to the widely used Cartesian Gaussian-based method of Obara and Saika (OS). Non-local interaction potentials such as standard Coulomb, modified Coulomb and Gaussian-type operators, that occur in range-separated hybrid functionals, are also included in the performance tests. The speed-up with respect to the OS scheme is up to three orders of magnitude for both, integrals and their derivatives. In particular, our method is increasingly efficient for large angular momenta and highly contracted basis sets., Comment: 18 pages, 2 figures; accepted manuscript. v2: supplementary material included

We determine and compare structural, dynamical, and electronic properties of liquid water at near ambient conditions through density-functional molecular dynamics simulations, when using either plane-wave or atomic-orbital basis sets. In both frameworks, the electronic structure and the atomic forces are self-consistently determined within the same theoretical scheme based on a nonlocal density functional accounting for van der Waals interactions. The overall properties of liquid water achieved within the two frameworks are in excellent agreement with each other. Thus, our study supports that implementations with plane-wave or atomic-orbital basis sets yield equivalent results and can be used indiscriminately in study of liquid water or aqueous solutions.

We present a molecular dynamics simulation study of CH$_3$NH$_3$PbI$_3$ based on forces calculated from density functional theory. The simulation were performed on model systems having 8 and 27 unit cells, and for a total simulation time of 40 ps in each case. Analysis of the finite size effects, in particular the mobility of the organic component, suggests that the smaller system is over correlated through the long range electrostatic interaction. In the larger system this finite size artifact is relaxed producing a more reliable description of the anisotropic rotational behavior of the methyl ammonium molecules. The thermal effects on the optical properties of the system were also analyzed. The HOMO-LUMO energy gap fluctuates around its central value with a standard deviation of approximately 0.1 eV. The projected density of states consistently place the Fermi level on the $p$ orbitals of the I atoms, and the lowest virtual state on $p$ orbitals of the Pb atoms throughout the whole simulation trajectory., Comment: 16 pages, 11 figures

Recently, the nonlocal van der Waals (vdW) density functionals [M. Dion, H. Rydberg, E. Schroeder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)] have attracted considerable attention due to their good performance for systems where weak interactions are important. Since the physics of dispersion is included in these functionals, they are usually more accurate and show less erratic behavior than the semilocal and hybrid methods. In this work, several variants of the vdW functionals have been tested on rare-gas dimers (from He2 to Kr2) and solids (Ne, Ar, and Kr) and their accuracy compared to standard semilocal approximations supplemented or not by an atom-pairwise dispersion correction [S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, J. Chem. Phys. 132, 154104 (2010)]. An analysis of the results in terms of energy decomposition is also provided.

Calculating highly accurate thermochemical properties of condensed matter via wave function-based approaches (such as e.g. Hartree-Fock or hybrid functionals) has recently attracted much interest. We here present two strategies providing accurate Hartree-Fock energies for solid LiH in a large Gaussian basis set and applying periodic boundary conditions. The total energies were obtained using two different approaches, namely a supercell evaluation of Hartree-Fock exchange using a truncated Coulomb operator and an extrapolation toward the full-range Hartree-Fock limit of a Pad\'e fit to a series of short-range screened Hartree-Fock calculations. These two techniques agreed to significant precision. We also present the Hartree-Fock cohesive energy of LiH (converged to within sub-meV) at the experimental equilibrium volume as well as the Hartree-Fock equilibrium lattice constant and bulk modulus., Comment: 7.5 pages, 2 figures, submitted to Phys. Rev. B; v2: typos removed, References added