Source author record

Florian Merz

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

ResearcherUnclaimed source record

physics.plasm-ph cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.stat-mech Data Structures and Algorithms physics.comp-ph

Catalog footprint

What is connected

6works

6topics

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.

preprint2021arXiv

All-electron periodic $G_0W_0$ implementation with numerical atomic orbital basis functions: algorithm and benchmarks

We present an all-electron, periodic {\GnWn} implementation within the numerical atomic orbital (NAO) basis framework. A localized variant of the resolution-of-the-identity (RI) approximation is employed to significantly reduce the computational cost of evaluating and storing the two-electron Coulomb repulsion integrals. We demonstrate that the error arising from localized RI approximation can be reduced to an insignificant level by enhancing the set of auxiliary basis functions, used to expand the products of two single-particle NAOs. An efficient algorithm is introduced to deal with the Coulomb singularity in the Brillouin zone sampling that is suitable for the NAO framework. We perform systematic convergence tests and identify a set of computational parameters, which can serve as the default choice for most practical purposes. Benchmark calculations are carried out for a set of prototypical semiconductors and insulators, and compared to independent reference values obtained from an independent $G_0W_0$ implementation based on linearized augmented plane waves (LAPW) plus high-energy localized orbitals (HLOs) basis set, as well as experimental results. With a moderate (FHI-aims \textit{tier} 2) NAO basis set, our $G_0W_0$ calculations produce band gaps that typically lie in between the standard LAPW and the LAPW+HLO results. Complementing \textit{tier} 2 with highly localized Slater-type orbitals (STOs), we find that the obtained band gaps show an overall convergence towards the LAPW+HLO results. The algorithms and techniques developed in this work pave the way for efficient implementations of correlated methods within the NAO framework.

preprint2020arXiv

NECI: N-Electron Configuration Interaction with emphasis on state-of-the-art stochastic methods

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24000 CPU cores. In this paper, we describe the core functionalities of NECI and recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs and it is licensed under GPL-3.0.

preprint2014arXiv

PReaCH: A Fast Lightweight Reachability Index using Pruning and Contraction Hierarchies

We develop the data structure PReaCH (for Pruned Reachability Contraction Hierarchies) which supports reachability queries in a directed graph, i.e., it supports queries that ask whether two nodes in the graph are connected by a directed path. PReaCH adapts the contraction hierarchy speedup techniques for shortest path queries to the reachability setting. The resulting approach is surprisingly simple and guarantees linear space and near linear preprocessing time. Orthogonally to that, we improve existing pruning techniques for the search by gathering more information from a single DFS-traversal of the graph. PReaCH-indices significantly outperform previous data structures with comparable preprocessing cost. Methods with faster queries need significantly more preprocessing time in particular for the most difficult instances.

preprint2011arXiv

Dynamic Procedure for Filtered Gyrokinetic Simulations

Large Eddy Simulations (LES) of gyrokinetic plasma turbulence are investigated as interesting candidates to decrease the computational cost. A dynamic procedure is implemented in the GENE code, allowing for dynamic optimization of the free parameters of the LES models (setting the amplitudes of dissipative terms). Employing such LES methods, one recovers the free energy and heat flux spectra obtained from highly resolved Direct Numerical Simulations (DNS). Systematic comparisons are performed for different values of the temperature gradient and magnetic shear, parameters which are of prime importance in Ion Temperature Gradient (ITG) driven turbulence. Moreover, the degree of anisotropy of the problem, that can vary with parameters, can be adapted dynamically by the method that shows Gyrokinetic Large Eddy Simulation (GyroLES) to be a serious candidate to reduce numerical cost of gyrokinetic solvers.

preprint2010arXiv

Impurity transport in ITG and TE mode dominated turbulence

In the present study, turbulent impurity transport in Deuterium tokamak plasmas, driven by Ion Temperature Gradient (ITG) and Trapped Electron (TE) modes, has been investigated using fluid and gyrokinetic models. The impurity diffusivity ($D_Z$) and convective velocity ($V_Z$) are calculated, and from these the zero-flux peaking factor ($PF_0$) is derived. This quantity expresses the impurity density gradient at which the convective and diffusive transport of impurities are exactly balanced. The sign of $PF_0$ is of special interest, as it determines whether the impurities are subject to an inward ($PF_0 > 0$) or outward ($PF_0 < 0$) pinch.

preprint2002arXiv

Negative scaling dimensions and conformal invariance at the Nishimori point in the +/-J random-bond Ising model

We reexamine the disorder-dominated multicritical point of the two-dimensional +/-J Ising model, known as the Nishimori point (NP). At the NP we investigate numerically and analytically the behavior of the disorder correlator, familiar from the self-dual description of the pure critical point of the two-dimensional Ising model. We consider the logarithmic average and the q-th moments of this correlator in the ensemble average over randomness, for continuous q in the range 0<q<2.5, and demonstrate their conformal invariance. At the NP we find, in contrast to the self-dual pure critical point, that the disorder correlators exhibit multi-scaling in q which is different from that of spin-spin correlators and that their scaling dimension becomes negative for q>1 and q<0. Using properties on the Nishimori line we show that the first moment (q=1) of the disorder correlator is exactly one for all separations. The spectrum of scaling dimensions at the NP is not parabolic in q.