BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Thomas D. Kühne

Thomas D. Kühne contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.comp-phphysics.chem-phcond-mat.mtrl-scicond-mat.dis-nncond-mat.stat-mechDistributed, Parallel, and Cluster Computingmath.NAquant-phcond-mat.mes-hallcond-mat.softcond-mat.str-elmath.RANumerical AnalysisPerformanceQuantitative Methods

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
15topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 published item(s)

preprint2026arXiv

Hitchhiker's guide to second-generation Car-Parrinello ab-initio molecular dynamics

In a recent letter [T. D. Kühne, M. Krack, F. Mohamed and M. Parrinello, Phys. Rev. Lett. 98, 066401 (2007)], we outlined a new Car-Parrinello-like approach to Born-Oppenheimer molecular dynamics. Here, we provide a guide to performing actual calculations using our method and demonstrate this on liquid water at ambient conditions. We do not go into methodological details beyond those necessary for applying this approach, but focus on practical details pertinent to our particular implementation within the CP2K/Quickstep code [T. D. Kühne et al., J. Chem. Phys. 152, 194103 (2020)].

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

The CP-PAW code package for first-principles calculations from a user's perspective

CP-PAW is a combined electronic structure and ab-initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave method with the Car-Parrinello approach to determine not only the electronic and nuclear ground state of condensed matter, but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Breaking the Exascale Barrier for the Electronic Structure Problem in Ab-Initio Molecular Dynamics

The non-orthogonal local submatrix method applied to electronic-structure based molecular dynamics simulations is shown to exceed 1.1 EFLOP/s in FP16/FP32 mixed floating-point arithmetic when using 4,400 NVIDIA A100 GPUs of the Perlmutter system. This is enabled by a modification of the original method that pushes the sustained fraction of the peak performance to about 80%. Example calculations are performed for SARS-CoV-2 spike proteins with up to 83 million atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

CP2K on the road to exascale

The CP2K program package, which can be considered as the swiss army knife of atomistic simulations, is presented with a special emphasis on ab-initio molecular dynamics using the second-generation Car-Parrinello method. After outlining current and near-term development efforts with regards to massively parallel low-scaling post-Hartree-Fock and eigenvalue solvers, novel approaches on how we plan to take full advantage of future low-precision hardware architectures are introduced. Our focus here is on combining our submatrix method with the approximate computing paradigm to address the immanent exascale era.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Parallel Quantum Chemistry on Noisy Intermediate-Scale Quantum Computers

A novel parallel hybrid quantum-classical algorithm for the solution of the quantum-chemical ground-state energy problem on gate-based quantum computers is presented. This approach is based on the reduced density-matrix functional theory (RDMFT) formulation of the electronic structure problem. For that purpose, the density-matrix functional of the full system is decomposed into an indirectly coupled sum of density-matrix functionals for all its subsystems using the adaptive cluster approximation to RDMFT. The approximations involved in the decomposition and the adaptive cluster approximation itself can be systematically converged to the exact result. The solutions for the density-matrix functionals of the effective subsystems involves a constrained minimization over many-particle states that are approximated by parametrized trial states on the quantum computer similarly to the variational quantum eigensolver. The independence of the density-matrix functionals of the effective subsystems introduces a new level of parallelization and allows for the computational treatment of much larger molecules on a quantum computer with a given qubit count. In addition, for the proposed algorithm techniques are presented to reduce the qubit count, the number of quantum programs, as well as its depth. The new approach is demonstrated for Hubbard-like systems on IBM quantum computers based on superconducting transmon qubits.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Towards Electronic Structure-Based Ab-Initio Molecular Dynamics Simulations with Hundreds of Millions of Atoms

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An automated approach for developing neural network interatomic potentials with FLAME

The performance of machine learning interatomic potentials relies on the quality of the training dataset. In this work, we present an approach for generating diverse and representative training data points which initiates with \it{ab initio} calculations for bulk structures. The data generation and potential construction further proceed side-by-side in a cyclic process of training the neural network and crystal structure prediction based on the developed interatomic potentials. All steps of the data generation and potential development are performed with minimal human intervention. We show the reliability of our approach by assessing the performance of neural network potentials developed for two inorganic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Artificial Neural Networks as Trial Wave Functions for Quantum Monte Carlo

Inspired by the universal approximation theorem and widespread adoption of artificial neural network techniques in a diversity of fields, we propose feed-forward neural networks as a general purpose trial wave function for quantum Monte Carlo simulations of continous many-body systems. Whereas for simple model systems the whole many-body wave function can be represented by a neural network, the antisymmetry condition of non-trivial fermionic systems is incorporated by means of a Slater determinant. To demonstrate the accuracy of our trial wave functions, we have studied an exactly solvable model system of two trapped interacting particles, as well as the hydrogen dimer.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accurate Sampling with Noisy Forces from Approximate Computing

In scientific computing, the acceleration of atomistic computer simulations by means of custom hardware is finding ever growing application. A major limitation, however, is that the high efficiency in terms of performance and low power consumption entails the massive usage of low-precision computing units. Here, based on the approximate computing paradigm, we present an algorithmic method to rigorously compensate for numerical inaccuracies due to low-accuracy arithmetic operations, yet still obtaining exact expectation values using a properly modified Langevin-type equation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Disordered Crystals from First Principles II: Transport Coefficients

This is the second part of a project on the foundations of first-principle calculations of the electron transport in crystals at finite temperatures, aiming at a predictive first-principles platform that combines ab-initio molecular dynamics (AIMD) and a finite-temperature Kubo-formula with dissipation for thermally disordered crystalline phases. The latter are encoded in an ergodic dynamical system $(Ω,\mathbb G,{\rm d}\mathbb P)$, where $Ω$ is the configuration space of the atomic degrees of freedom, $\mathbb G$ is the space group acting on $Ω$ and ${\rm d}\mathbb P$ is the ergodic Gibbs measure relative to the $\mathbb G$-action. We first demonstrate how to pass from the continuum Kohn-Sham theory to a discrete atomic-orbitals based formalism without breaking the covariance of the physical observables w.r.t. $(Ω,\mathbb G,{\rm d}\mathbb P)$. Then we show how to implement the Kubo-formula, investigate its self-averaging property and derive an optimal finite-volume approximation for it. We also describe a numerical innovation that made possible AIMD simulations with longer orbits and elaborate on the details of our simulations. Lastly, we present numerical results on the transport coefficients of crystal silicon at different temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient Ab-Initio Molecular Dynamic Simulations by Offloading Fast Fourier Transformations to FPGAs

A large share of today's HPC workloads is used for Ab-Initio Molecular Dynamics (AIMD) simulations, where the interatomic forces are computed on-the-fly by means of accurate electronic structure calculations. They are computationally intensive and thus constitute an interesting application class for energy-efficient hardware accelerators such as FPGAs. In this paper, we investigate the potential of offloading 3D Fast Fourier Transformations (FFTs) as a critical routine of plane-wave-based electronic structure calculations to FPGA and in conjunction demonstrate the tolerance of these simulations to lower precision computations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Energy Transfer within the Hydrogen Bonding Network of Water Following Resonant Terahertz Excitation

Energy dissipation in water is very fast and more efficient than in many other liquids. This behavior is commonly attributed to the intermolecular interactions associated with hydrogen bonding. Here, we investigate the dynamic energy flow in the hydrogen-bond network of liquid water by a pump-probe experiment. We resonantly excite intermolecular degrees of freedom with ultrashort single-cycle terahertz pulses and monitor its Raman response. By using ultrathin sample-cell windows, a background-free bipolar signal whose tail relaxes mono-exponentially is obtained. The relaxation is attributed to the molecular translational motions, using complementary experiments, force-field and ab initio molecular dynamics simulations. They reveal an initial coupling of the terahertz electric field to the molecular rotational degrees of freedom whose energy is rapidly transferred, within the excitation pulse duration, to the restricted-translational motion of neighboring molecules. This rapid energy transfer may be rationalized by the strong anharmonicity of the intermolecular interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modulation of Nearly Free Electron States in Hydroxyl-Functionalized MXenes: A First-Principles Study

The transition metal carbides (namely MXenes) and their functionalized derivatives exhibit various physical and chemical characteristics and offer many potential applications in electronic devices and sensors. Using density functional theory (DFT), it is revealed that the nearly free electron (NFE) states are near the Fermi levels in hydroxyl (OH) functionalized MXenes. Most of the OH-terminated MXene are metallic, but some of them, e.g. Sc2C(OH)2, are semiconductors and the NFE states are conduction bands. In this paper, to investigate the NFE states in MXenes, an attractive image-potential well model is adopted. Compared the solutions of this model with the DFT calculations, it is found that due to the overlap of spatially extensive wave functions of NFE states and their hybridization between the artificial neighboring layers imposed by the periodical boundary conditions (PBCs), the DFT results represent the properties of multiple layers, intrinsically. Based on the DFT calculations, it is found that the energy gap widths are affected by the interlayer distances. We address that the energetics of the NFE states can be modulated by the external electric fields and it is possible to convert semiconducting MXenes into metals. This band-gap manipulation makes the OH-terminated semiconducting MXenes an excellent candidate for electronic switch applications. Finally, using a set of electron transport calculations, I-V characteristics of Sc2C(OH)2 devices are investigated with the gate voltages.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tumbling with a Limp: Local Asymmetry in Water's Hydrogen Bond Network and its Consequences

Ab initio molecular dynamics simulations of liquid water under equilibrium ambient conditions, together with a novel energy decomposition analysis, have recently shown that a substantial fraction of water molecules exhibit a significant asymmetry between the strengths of the two donor and/or the two acceptor interactions. We refer to this recently unraveled aspect as the "local asymmetry in the hydrogen bond network". We discuss how this novel aspect was first revealed, and provide metrics that can be consistently employed on simulated water trajectories to quantify this local heterogeneity in the hydrogen bond network and its dynamics. We then discuss the static aspects of the asymmetry, pertaining to the frozen geometry of liquid water at any given instant of time and the distribution of hydrogen bond strengths therein, and also its dynamic characteristics pertaining to how fast this asymmetry decays and the kinds of molecular motions responsible for this decay. Following this we discuss the spectroscopic manifestations of this asymmetry, from ultrafast X-ray absorption spectra to infrared spectroscopy and down to the much slower terahertz regime. Finally, we discuss the implications of these findings in a broad context and its relation to the current notions about the structure and dynamics of liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Water structure near the surface of Weyl semimetals as catalysts in photocatalytic proton reduction

In this work, second-generation Car-Parrinello-based QM/MM molecular dynamics simulations of small nanoparticles of NbP, NbAs, TaAs and 1T-TaS$_2$ in water are presented. The first three materials are topological Weyl semimetals, which were recently discovered to be active catalysts in photocatalytic water splitting. The aim of this research was to correlate potential differences in the water structure in the vicinity of the nanoparticle surface with the photocatalytic activity of these materials in light induced proton reduction. The results presented herein allow to explain the catalytic activity of these Weyl semimetals: the most active material, NbP, exhibits a particularly low water coordination near the surface of the nanoparticle, whereas for 1T-TaS$_2$, with the lowest catalytic activity, the water structure at the surface is most ordered. In addition, the photocatalytic activity of several organic and metalorganic photosensitizers in the hydrogen evolution reaction was experimentally investigated with NbP as proton reduction catalyst. Unexpectedly, the charge of the photosensitizer plays a decisive role for the photocatalytic performance.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A General Algorithm to Calculate the Inverse Principal $p$-th Root of Symmetric Positive Definite Matrices

We address the general mathematical problem of computing the inverse $p$-th root of a given matrix in an efficient way. A new method to construct iteration functions that allow calculating arbitrary $p$-th roots and their inverses of symmetric positive definite matrices is presented. We show that the order of convergence is at least quadratic and that adaptively adjusting a parameter $q$ always leads to an even faster convergence. In this way, a better performance than with previously known iteration schemes is achieved. The efficiency of the iterative functions is demonstrated for various matrices with different densities, condition numbers and spectral radii.

0 citations0 reviewsNo code linkedOpen access