Source author record

Anders M. N. Niklasson

Anders M. N. Niklasson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.chem-ph cond-mat.mtrl-sci physics.comp-ph math.NA cond-mat.other Numerical Analysis Performance quant-ph

Catalog footprint

What is connected

15works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Quantum perturbation theory using Tensor cores and a deep neural network

Time-independent quantum response calculations are performed using Tensor cores. This is achieved by mapping density matrix perturbation theory onto the computational structure of a deep neural network. The main computational cost of each deep layer is dominated by tensor contractions, i.e. dense matrix-matrix multiplications, in mixed precision arithmetics which achieves close to peak performance. Quantum response calculations are demonstrated and analyzed using self-consistent charge density-functional tight-binding theory as well as coupled-perturbed Hartree-Fock theory. For linear response calculations, a novel parameter-free convergence criterion is presented that is well-suited for numerically noisy low precision floating point operations and we demonstrate a peak performance of almost 200 Tflops using the Tensor cores of two Nvidia A100 GPUs.

preprint2021arXiv

Performance Optimizations of Recursive Electronic Structure Solvers targeting Multi-Core Architectures (LA-UR-20-26665)

As we rapidly approach the frontiers of ultra large computing resources, software optimization is becoming of paramount interest to scientific application developers interested in efficiently leveraging all available on-Node computing capabilities and thereby improving a requisite science per watt metric. The scientific application of interest here is the Basic Math Library (BML) that provides a singular interface for linear algebra operation frequently used in the Quantum Molecular Dynamics (QMD) community. The provisioning of a singular interface indicates the presence of an abstraction layer which in-turn suggests commonalities in the code-base and therefore any optimization or tuning introduced in the core of code-base has the ability to positively affect the performance of the aforementioned library as a whole. With that in mind, we proceed with this investigation by performing a survey of the entirety of the BML code-base, and extract, in form of micro-kernels, common snippets of code. We introduce several optimization strategies into these micro-kernels including 1.) Strength Reduction 2.) Memory Alignment for large arrays 3.) Non Uniform Memory Access (NUMA) aware allocations to enforce data locality and 4.) appropriate thread affinity and bindings to enhance the overall multi-threaded performance. After introducing these optimizations, we benchmark the micro-kernels and compare the run-time before and after optimization for several target architectures. Finally we use the results as a guide to propagating the optimization strategies into the BML code-base. As a demonstration, herein, we test the efficacy of these optimization strategies by comparing the benchmark and optimized versions of the code.

preprint2021arXiv

Shadow Lagrangian dynamics for superfluidity

Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without any additional costs. The findings are illustrated in numerical experiments.

preprint2020arXiv

Density-matrix based Extended Lagrangian Born-Oppenheimer Molecular Dynamics

Extended Lagrangian Born-Oppenheimer molecular dynamics [{\em Phys.\ Rev.\ Lett.\ } {\bf 2008}, {\em 100}, 123004] is presented for Hartree-Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born-Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born-Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, ${\cal K}$, that determines the metric tensor, ${\cal T}\equiv {\cal K}^T{\cal K}$, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. The formulation and algorithms provide a general guide to implement extended Lagrangian Born-Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

preprint2020arXiv

Modeling solid-liquid interface reactions with next generation extended Lagrangian quantum-based molecular dynamics

We demonstrate the applicability of extended Lagrangian Born-Oppenheimer quantum-based molecular dynamics (XL-BOMD) to model electron transfer reactions occurring on solid-liquid interfaces. Specifically, we consider the reduction of O$_2$ as catalyzed at the interface of an N-doped graphene sheet and H$_2$O at fuel cell cathodes. This system is a good testbed for next-generation computational chemistry methods since the electrochemical functionalities strongly depend on atomic-scale quantum mechanics. As opposed to prior iterations of first principles molecular dynamics, XL-BOMD only requires a full self-consistent-charge relaxation during the initial time step. The electronic ground state and total energy are stabilized thereafter through nuclear and electronic equations of motion assisted by an inner-product kernel updated with low-rank approximations. A species charge analysis reveals that the kernel-based XL-BOMD simulation can capture an electron transfer between the PGM-free catalyst and a solvated O$_2$ molecule mediated by H$_2$O, which results in the molecular dissociation of O$_2$.

preprint2016arXiv

Graph-based linear scaling electronic structure theory

We show how graph theory can be combined with quantum theory to calculate the electronic structure of large complex systems. The graph formalism is general and applicable to a broad range of electronic structure methods and materials, including challenging systems such as biomolecules. The methodology combines well-controlled accuracy, low computational cost, and natural low-communication parallelism. This combination addresses substantial shortcomings of linear scaling electronic structure theory, in particular with respect to quantum-based molecular dynamics simulations.

preprint2015arXiv

Canonical density matrix perturbation theory

Density matrix perturbation theory [Niklasson and Challacombe, Phys. Rev. Lett. 92, 193001 (2004)] is generalized to canonical (NVT) free energy ensembles in tight-binding, Hartree-Fock or Kohn-Sham density functional theory. The canonical density matrix perturbation theory can be used to calculate temperature dependent response properties from the coupled perturbed self-consistent field equations as in density functional perturbation theory. The method is well suited to take advantage of sparse matrix algebra to achieve linear scaling complexity in the computational cost as a function of system size for sufficiently large non-metallic materials and metals at high temperatures.

preprint2014arXiv

First principles molecular dynamics without self-consistent field optimization

We present a first principles molecular dynamics approach that is based on time-reversible ex- tended Lagrangian Born-Oppenheimer molecular dynamics [Phys. Rev. Lett. 100, 123004 (2008)] in the limit of vanishing self-consistent field optimization. The optimization-free dynamics keeps the computational cost to a minimum and typically provides molecular trajectories that closely follow the exact Born-Oppenheimer potential energy surface. Only one single diagonalization and Hamiltonian (or Fockian) costruction are required in each integration time step. The proposed dy- namics is derived for a general free-energy potential surface valid at finite electronic temperatures within hybrid density functional theory. Even in the event of irregular functional behavior that may cause a dynamical instability, the optimization-free limit represents an ideal starting guess for force calculations that may require a more elaborate iterative electronic ground state optimization. Our optimization-free dynamics thus represents a flexible theoretical framework for a broad and general class of ab initio molecular dynamics simulations.

preprint2014arXiv

Generalized Extended Lagrangian Born-Oppenheimer Molecular Dynamics

Extended Lagrangian Born-Oppenheimer molecular dynamics based on Kohn-Sham density functional theory is generalized in the limit of vanishing self-consistent field optimization prior to the force evaluations. The equations of motion are derived directly from the extended Lagrangian under the condition of an adiabatic separation between the nuclear and the electronic degrees of freedom. We show how this separation is automatically fulfilled and system independent. The generalized equations of motion require only one diagonalization per time step and are applicable to a broader range of materials with improved accuracy and stability compared to previous formulations.

preprint2013arXiv

Accelerated density matrix expansions for Born-Oppenheimer molecular dynamics

An accelerated polynomial expansion scheme to construct the density matrix in quantum mechanical molecular dynamics simulations is proposed. The scheme is based on recursive density matrix expansions, e.g. [Phys. Rev. B. 66 (2002), p. 155115], which are accelerated by a scale-and-fold technique [J. Chem. Theory Comput. 7 (2011), p. 1233]. The acceleration scheme requires interior eigenvalue estimates, which may be expensive and cumbersome to come by. Here we show how such eigenvalue estimates can be extracted from the recursive expansion by a simple and robust procedure at a negligible computational cost. Our method is illustrated with density functional tight-binding Born-Oppenheimer molecular dynamics simulations, where the computational effort is dominated by the density matrix construction. In our analysis we identify two different phases of the recursive polynomial expansion, the conditioning and purification phases, and we show that the acceleration represents an improvement of the conditioning phase, which typically gives a significant reduction of the computational cost.

preprint2013arXiv

Extended Lagrangian Born-Oppenheimer molecular dynamics in the limit of vanishing self-consistent field optimization

We present an efficient general approach to first principles molecular dynamics simulations based on extended Lagrangian Born-Oppenheimer molecular dynamics in the limit of vanishing self-consistent field optimization. The reduction of the optimization requirement reduces the computational cost to a minimum, but without causing any significant loss of accuracy or longterm energy drift. The optimization-free first principles molecular dynamics requires only one single diagonalization per time step and yields trajectories at the same level of accuracy as "exact", fully converged, Born-Oppenheimer molecular dynamics simulations. The optimization-free limit of extended Lagrangian Born-Oppenheimer molecular dynamics therefore represents an ideal starting point for a robust and efficient formulation of a new generation first principles quantum mechanical molecular dynamics simulation schemes.

preprint2012arXiv

Fast method for quantum mechanical molecular dynamics

With the continuous growth of processing power for scientific computing, first principles Born-Oppenheimer molecular dynamics (MD) simulations are becoming increasingly popular for the study of a wide range of problems in materials science, chemistry and biology. Nevertheless, the computational cost still remains prohibitively large in many cases, particularly in comparison to classical MD simulations using empirical force fields. Here we show how to circumvent the major computational bottleneck in Born-Oppenheimer MD simulations arising from the self-consistent-charge optimization. The optimization-free quantum mechanical MD method is demonstrated for density functional tight-binding theory. The molecular trajectories are almost indistinguishable from an "exact" microcanonical Born-Oppenheimer MD simulation even when linear scaling sparse matrix algebra is used. Our findings drastically reduce the computational gap between classical and quantum mechanical MD simulations.

preprint2011arXiv

Extended Lagrangian free energy molecular dynamics

Extended free energy Lagrangians are proposed for first principles molecular dynamics simulations at finite electronic temperatures for plane-wave pseudopotential and local orbital density matrix based calculations. Thanks to the extended Lagrangian description the electronic degrees of freedom can be integrated by stable geometric schemes that conserve the free energy. For the local orbital representations both the nuclear and electronic forces have simple and numerically efficient expressions that are well suited for reduced complexity calculations. A rapidly converging recursive Fermi operator expansion method that does not require the calculation of eigenvalues and eigenfunctions for the construction of the fractionally occupied density matrix is discussed. An efficient expression for the Pulay force that is valid also for density matrices with fractional occupation occurring at finite electronic temperatures is also demonstrated.

preprint2010arXiv

Wavefunction extended Lagrangian Born-Oppenheimer molecular dynamics

Extended Lagrangian Born-Oppenheimer molecular dynamics [Niklasson, Phys. Rev. Lett. 100 123004 (2008)] has been generalized to the propagation of the electronic wavefunctions. The technique allows highly efficient first principles molecular dynamics simulations using plane wave pseudopotential electronic structure methods that are stable and energy conserving also under incomplete and approximate self-consistency convergence. An implementation of the method within the planewave basis set is presented and the accuracy and efficiency is demonstrated both for semi-conductor and metallic materials.

preprint2008arXiv

Molecular-orbital-free algorithm for excited states in time-dependent perturbation theory

A non-linear conjugate gradient optimization scheme is used to obtain excitation energies within the Random Phase Approximation (RPA). The solutions to the RPA eigenvalue equation are located through a variational characterization using a modified Thouless functional, which is based upon an asymmetric Rayleigh quotient, in an orthogonalized atomic orbital representation. In this way, the computational bottleneck of calculating molecular orbitals is avoided. The variational space is reduced to the physically-relevant transitions by projections. The feasibility of an RPA implementation scaling linearly with system size, N, is investigated by monitoring convergence behavior with respect to the quality of initial guess and sensitivity to noise under thresholding, both for well- and ill-conditioned problems. The molecular- orbital-free algorithm is found to be robust and computationally efficient providing a first step toward a large-scale, reduced complexity calculation of time-dependent optical properties and linear response. The algorithm is extensible to other forms of time-dependent perturbation theory including, but not limited to, time-dependent Density Functional theory.