Source author record

Stephan Mohr

Stephan Mohr appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.mtrl-sci physics.chem-ph physics.comp-ph cond-mat.mes-hall cond-mat.other cond-mat.stat-mech Distributed, Parallel, and Cluster Computing

Catalog footprint

What is connected

15works

7topics

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

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.

preprint2022arXiv

Unveiling the Zero-Phonon Line of the Boron Vacancy Center by Cavity Enhanced Emission

Negatively charged boron vacancies ($V_B^-$) in hexagonal boron nitride (hBN) exhibit a broad emission spectrum due to strong electron-phonon coupling and Jahn-Teller mixing of electronic states. As such, the direct measurement of zero-phonon line (ZPL) of $V_B^-$ has remained elusive. Here, we measure the room-temperature ZPL wavelength to be $773\pm2$ nm by coupling the hBN layer to the high-Q nanobeam cavity. As the wavelength of cavity mode is tuned, we observe a pronounced intensity resonance, indicating the coupling to $V_B^-$. Our observations are consistent with the spatial redistribution of $V_B^-$ emission. Spatially resolved measurements show a clear Purcell effect maximum at the midpoint of the nanobeam, in accord with the optical field distribution of the cavity mode. Our results are in good agreement with theoretical calculations, opening the way to using $V_B^-$ as cavity spin-photon interfaces.

preprint2021arXiv

Lifetime effects and satellites in the photoelectron spectrum of tungsten metal

Tungsten is an important and versatile transition metal and has a firm place at the heart of many technologies. A popular experimental technique for the characterisation of tungsten and tungsten-based compounds is X-ray photoelectron spectroscopy (XPS), which enables the assessment of chemical states and electronic structure through the collection of core level and valence band spectra. However, in the case of metallic tungsten, open questions remain regarding the origin, nature, and position of satellite features that are prominent in the photoelectron spectrum. These satellites are a fingerprint of the electronic structure of the material and have not been thoroughly investigated, at times leading to their misinterpretation. The present work combines high-resolution soft and hard X-ray photoelectron spectroscopy (SXPS and HAXPES) with reflection electron energy loss spectroscopy (REELS) and a multi-tiered ab-initio theoretical approach, including density functional theory (DFT) and many-body perturbation theory (G0W0 and GW+C), to disentangle the complex set of experimentally observed satellite features attributed to the generation of plasmons and interband transitions. This combined experiment-theory strategy is able to uncover previously undocumented satellite features, improving our understanding of their direct relationship to tungsten's electronic structure. Furthermore, it lays the groundwork for future studies into tungsten based mixed-metal systems and holds promise for the re-assessment of the photoelectron spectra of other transition and post-transition metals, where similar questions regarding satellite features remain.

preprint2021arXiv

Ranking the efficiency of gas hydrate anti agglomerants through molecular dynamic simulations

Using both computational and experimental methods, the capacity of four different surfactant molecules to inhibit the agglomeration of sII hydrate particles was assessed. The computational simulations were carried out using both steered and non-steered Molecular Dynamics (MD), simulating the coalescence process of a hydrate slab and a water droplet, both covered with surfactant molecules. The experimental work was based on rocking cell measurements, determining the minimum effective dose necessary to inhibit agglomeration. Overall, good agreement was obtained between the performance predicted by the simulations and the experimental measurements. Moreover, the simulations allowed to gain additional insights that are not directly accessible via experiments, such as an analysis of the mass density profiles, the diffusion coefficients, or the orientations of the long tails.

preprint2020arXiv

Complexity Reduction in Density Functional Theory Calculations of Large Systems: System Partitioning and Fragment Embedding

With the development of low order scaling methods for performing Kohn-Sham Density Functional Theory, it is now possible to perform fully quantum mechanical calculations of systems containing tens of thousands of atoms. However, with an increase in the size of system treated comes an increase in complexity, making it challenging to analyze such large systems and determine the cause of emergent properties. To address this issue, in this paper we present a systematic complexity reduction methodology which can break down large systems into their constituent fragments, and quantify inter-fragment interactions. The methodology proposed here requires no a priori information or user interaction, allowing a single workflow to be automatically applied to any system of interest. We apply this approach to a variety of different systems, and show how it allows for the derivation of new system descriptors, the design of QM/MM partitioning schemes, and the novel application of graph metrics to molecules and materials.

preprint2020arXiv

SIESTA: recent developments and applications

A review of the present status, recent enhancements, and applicability of the SIESTA program is presented. Since its debut in the mid-nineties, SIESTA's flexibility, efficiency and free distribution has given advanced materials simulation capabilities to many groups worldwide. The core methodological scheme of SIESTA combines finite-support pseudo-atomic orbitals as basis sets, norm-conserving pseudopotentials, and a real-space grid for the representation of charge density and potentials and the computation of their associated matrix elements. Here we describe the more recent implementations on top of that core scheme, which include: full spin-orbit interaction, non-repeated and multiple-contact ballistic electron transport, DFT+U and hybrid functionals, time-dependent DFT, novel reduced-scaling solvers, density-functional perturbation theory, efficient Van der Waals non-local density functionals, and enhanced molecular-dynamics options. In addition, a substantial effort has been made in enhancing interoperability and interfacing with other codes and utilities, such as Wannier90 and the second-principles modelling it can be used for, an AiiDA plugin for workflow automatization, interface to Lua for steering SIESTA runs, and various postprocessing utilities. SIESTA has also been engaged in the Electronic Structure Library effort from its inception, which has allowed the sharing of various low level libraries, as well as data standards and support for them, in particular the PSML definition and library for transferable pseudopotentials, and the interface to the ELSI library of solvers. Code sharing is made easier by the new open-source licensing model of the program. This review also presents examples of application of the capabilities of the code, as well as a view of on-going and future developments.

preprint2020arXiv

The CECAM Electronic Structure Library and the modular software development paradigm

First-principles electronic structure calculations are very widely used thanks to the many successful software packages available. Their traditional coding paradigm is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, from the compiler up, with the exception of linear-algebra and message-passing libraries. This model has been quite successful for decades. The rapid progress in methodology, however, has resulted in an ever increasing complexity of those programs, which implies a growing amount of replication in coding and in the recurrent re-engineering needed to adapt to evolving hardware architecture. The Electronic Structure Library (\esl) was initiated by CECAM (European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure programs and redesign them as free, open-source libraries. They include "heavy-duty" ones with a high degree of parallelisation, and potential for adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by scientists when implementing new ideas. It is a community effort, undertaken by developers of various successful codes, now facing the challenges arising in the new model. This modular paradigm will improve overall coding efficiency and enable specialists (computer scientists or computational scientists) to use their skills more effectively. It will lead to a more sustainable and dynamic evolution of software as well as lower barriers to entry for new developers.

preprint2016arXiv

Challenges in Large Scale Quantum Mechanical Calculations

During the past decades, quantum mechanical methods have undergone an amazing transition from pioneering investigations of experts into a wide range of practical applications, made by a vast community of researchers. First principles calculations of systems containing up to a few hundred atoms have become a standard in many branches of science. The sizes of the systems which can be simulated have increased even further during recent years, and quantum-mechanical calculations of systems up to many thousands of atoms are nowadays possible. This opens up new appealing possibilities, in particular for interdisciplinary work, bridging together communities of different needs and sensibilities. In this review we will present the current status of this topic, and will also give an outlook on the vast multitude of applications, challenges and opportunities stimulated by electronic structure calculations, making this field an important working tool and bringing together researchers of many different domains.

preprint2015arXiv

Accurate and efficient linear scaling DFT calculations with universal applicability

Density Functional Theory calculations traditionally suffer from an inherent cubic scaling with respect to the size of the system, making big calculations extremely expensive. This cubic scaling can be avoided by the use of so-called linear scaling algorithms, which have been developed during the last few decades. In this way it becomes possible to perform ab-initio calculations for several tens of thousands of atoms or even more within a reasonable time frame. However, even though the use of linear scaling algorithms is physically well justified, their implementation often introduces some small errors. Consequently most implementations offering such a linear complexity either yield only a limited accuracy or, if one wants to go beyond this restriction, require a tedious fine tuning of many parameters. In our linear scaling approach within the BigDFT package, we were able to overcome this restriction. Using an ansatz based on localized support functions expressed in an underlying Daubechies wavelet basis -- which offers ideal properties for accurate linear scaling calculations -- we obtain an amazingly high accuracy and a universal applicability while still keeping the possibility of simulating large systems with only a moderate demand of computing resources.

preprint2015arXiv

Fragment Approach to Constrained Density Functional Theory Calculations using Daubechies Wavelets

In a recent paper we presented a linear scaling Kohn-Sham density functional theory (DFT) code based on Daubechies wavelets, where a minimal set of localized support functions is optimized in situ and therefore adapted to the chemical properties of the molecular system. Thanks to the systematically controllable accuracy of the underlying basis set, this approach is able to provide an optimal contracted basis for a given system: accuracies for ground state energies and atomic forces are of the same quality as an uncontracted, cubic scaling approach. This basis set offers, by construction, a natural subset where the density matrix of the system can be projected. In this paper we demonstrate the flexibility of this minimal basis formalism in providing a basis set that can be reused as-is, i.e. without reoptimization, for charge-constrained DFT calculations within a fragment approach. Support functions, represented in the underlying wavelet grid, of the template fragments are roto-translated with high numerical precision to the required positions and used as projectors for the charge weight function. We demonstrate the interest of this approach to express highly precise and efficient calculations for preparing diabatic states and for the computational setup of systems in complex environments.

preprint2014arXiv

Daubechies Wavelets for Linear Scaling Density Functional Theory

We demonstrate that Daubechies wavelets can be used to construct a minimal set of optimized localized contracted basis functions in which the Kohn-Sham orbitals can be represented with an arbitrarily high, controllable precision. Ground state energies and the forces acting on the ions can be calculated in this basis with the same accuracy as if they were calculated directly in a Daubechies wavelets basis, provided that the amplitude of these contracted basis functions is sufficiently small on the surface of the localization region, which is guaranteed by the optimization procedure described in this work. This approach reduces the computational costs of DFT calculations, and can be combined with sparse matrix algebra to obtain linear scaling with respect to the number of electrons in the system. Calculations on systems of 10,000 atoms or more thus become feasible in a systematic basis set with moderate computational resources. Further computational savings can be achieved by exploiting the similarity of the contracted basis functions for closely related environments, e.g. in geometry optimizations or combined calculations of neutral and charged systems.

preprint2014arXiv

Minima Hopping Guided Path Search: An Efficient Method for Finding Complex Chemical Reaction Pathways

The Minima Hopping global optimization method uses physically realizable molecular dynamics moves in combination with an energy feedback that guarantees the escape from any potential energy funnel. For the purpose of finding reactions pathways, we argue that Minima Hopping is particularly suitable as a guide through the potential energy landscape and as a generator for pairs of minima that can be used as input structures for methods capable of finding transition states between two minima. For Lennard-Jones benchmark systems we compared this Minima Hopping guided path search method to a known approach for the exploration of potential energy landscapes that is based on deterministic mode-following. Although we used a stabilized mode-following technique that reliably allows to follow distinct directions when escaping from a local minimum, we observed that Minima Hopping guided path search is far superior in finding lowest-barrier reaction pathways. We therefore suggest that Minima Hopping guided path search can be used as a simple and efficient way to identify energetically low-lying chemical reaction pathways. Finally we applied the Minima Hopping guided path search approach to 75-atom and 102-atom Lennard Jones systems. For the 75-atom system we found pathways whose highest energies are significantly lower than the highest energy along the previously published lowest-barrier pathway. Furthermore, many of these pathways contain a smaller number of intermediate transition states than the previously publish lowest-barrier pathway. In case of the 102-atom system Minima Hopping guided path search found a previously unknown and energetically low-lying funnel.

preprint2013arXiv

Boron aggregation in the ground states of boron-carbon fullerenes

We present novel structural motifs for boron-carbon nano-cages of the stochiometries B12C48 and B12C50, based on first principle calculations. These configurations are distinct from those proposed so far by the fact that the boron atoms are not isolated and distributed over the entire surface of the cages, but rather aggregate at one location to form a patch. Our putative ground state of B12C48 is 1.8 eV lower in energy than the previously proposed ground state and violates all the suggested empirical rules for constructing low energy fullerenes. The B12C50 configuration is energetically even more favorable than B12C48, showing that structures derived from the C60 Buckminsterfullerene are not necessarily magic sizes for heterofullerenes.

preprint2013arXiv

Metrics for measuring distances in configuration spaces

In order to characterize molecular structures we introduce configurational fingerprint vectors which are counterparts of quantities used experimentally to identify structures. The Euclidean distance between the configurational fingerprint vectors satisfies the properties of a metric and can therefore safely be used to measure dissimilarities between configurations in the high dimensional configuration space. We show that these metrics correlate well with the RMSD between two configurations if this RMSD is obtained from a global minimization over all translations, rotations and permutations of atomic indices. We introduce a Monte Carlo approach to obtain this global minimum of the RMSD between configurations.

preprint2010arXiv

Efficient moves for global geometry optimization methods and their application to binary systems

We show that molecular dynamics based moves in the Minima Hopping (MH) method are more efficient than saddle point crossing moves which select the lowest possible saddle point. For binary systems we incorporate identity exchange moves in a way that allows to avoid the generation of high energy configurations. Using this modified Minima Hopping method we reexamine the binary Lennard Jones (BLJ) benchmark system with up to 100 atoms and we find a large number of new putative global minima structures.

Stephan Mohr

What is connected

Connect this record

See the researcher in context

Building this map preview

15 published item(s)

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

Unveiling the Zero-Phonon Line of the Boron Vacancy Center by Cavity Enhanced Emission

Lifetime effects and satellites in the photoelectron spectrum of tungsten metal

Ranking the efficiency of gas hydrate anti agglomerants through molecular dynamic simulations

Complexity Reduction in Density Functional Theory Calculations of Large Systems: System Partitioning and Fragment Embedding

SIESTA: recent developments and applications

The CECAM Electronic Structure Library and the modular software development paradigm

Challenges in Large Scale Quantum Mechanical Calculations

Accurate and efficient linear scaling DFT calculations with universal applicability

Fragment Approach to Constrained Density Functional Theory Calculations using Daubechies Wavelets

Daubechies Wavelets for Linear Scaling Density Functional Theory

Minima Hopping Guided Path Search: An Efficient Method for Finding Complex Chemical Reaction Pathways

Boron aggregation in the ground states of boron-carbon fullerenes

Metrics for measuring distances in configuration spaces

Efficient moves for global geometry optimization methods and their application to binary systems