Source author record

M. Puhr

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
5topics
2close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

6 published item(s)

preprint2016arXiv

A method to compute derivatives of functions of large complex matrices

A recently developed numerical method for the calculation of derivatives of functions of general complex matrices, which can also be combined with implicit matrix function approximations such as Krylov-Ritz type algorithms, is presented. An important use case for the method in the context of lattice gauge theory is the overlap Dirac operator at finite quark chemical potential. Derivatives of the lattice Dirac operator are necessary for the computation of conserved lattice currents or the fermionic force in Hybrid Monte-Carlo and Langevin simulations. To calculate the overlap Dirac operator at finite chemical potential the product of the sign function of a non-Hermitian matrix with a vector has to be computed. For non-Hermitian matrices it is not possible to efficiently approximate the sign function with polynomials or rational functions. Implicit approximation algorithms, like Krylov-Ritz methods, that depend on the source vector have to be used instead. Our method can also provide derivatives of such implicit approximations. We show how a generalised deflation prescription can be used to improve the performance of the method, if some eigenvalues and eigenvectors of the matrix being differentiated are known. To show that the method is efficient and well suited for practical calculations we provide test results for the two-sided Lanczos approximation of the finite-density overlap Dirac operator on SU(3) gauge field configurations on lattices with sizes up to $14 \times 14^3$.

preprint2016arXiv

A numerical method to compute derivatives of functions of large complex matrices and its application to the overlap Dirac operator at finite chemical potential

We present a method for the numerical calculation of derivatives of functions of general complex matrices. The method can be used in combination with any algorithm that evaluates or approximates the desired matrix function, in particular with implicit Krylov-Ritz-type approximations. An important use case for the method is the evaluation of the overlap Dirac operator in lattice Quantum Chromodynamics (QCD) at finite chemical potential, which requires the application of the sign function of a non-Hermitian matrix to some source vector. While the sign function of non-Hermitian matrices in practice cannot be efficiently approximated with source-independent polynomials or rational functions, sufficiently good approximating polynomials can still be constructed for each particular source vector. Our method allows for an efficient calculation of the derivatives of such implicit approximations with respect to the gauge field or other external parameters, which is necessary for the calculation of conserved lattice currents or the fermionic force in Hybrid Monte-Carlo or Langevin simulations. We also give an explicit deflation prescription for the case when one knows several eigenvalues and eigenvectors of the matrix being the argument of the differentiated function. We test the method for the two-sided Lanczos approximation of the finite-density overlap Dirac operator on realistic $SU(3)$ gauge field configurations on lattices with sizes as large as $14\times14^3$ and $6\times18^3$.

preprint2016arXiv

Chiral Magnetic Effect in finite-size samples of parity-breaking Weyl semimetals

We study the static electric current due to the Chiral Magnetic Effect in samples of Weyl semimetals with slab geometry, where the magnetic field is parallel to the boundaries of the slab. We use the Wilson-Dirac Hamiltonian as a simplest model of parity-breaking Weyl semimetal with two-band structure. We find that the CME current is strongly localized at the open boundaries of the slab, where the current density in the direction of the magnetic field approaches the conventional value $j = μ_5 B /2 π^2$ at sufficiently small values of the chiral chemical potential $μ_5$ and magnetic field $B$. On the other hand, very large values of magnetic field tend to suppress the CME response. We observe that the localization width of the current is independent of the slab width and is given by the magnetic length $l_B = 1/\sqrt{B}$ when $μ_5 \ll \sqrt{B}$. In the opposite regime when $μ_5 \gg \sqrt{B}$ the localization width is determined solely by $μ_5$

preprint2015arXiv

Chiral magnetic conductivity in an interacting lattice model of parity-breaking Weyl semimetal

We report on the mean-field study of the Chiral Magnetic Effect (CME) in static magnetic fields within a simple model of a parity-breaking Weyl semimetal given by the lattice Wilson-Dirac Hamiltonian with constant chiral chemical potential. We consider both the mean-field renormalization of the model parameters and nontrivial corrections to the CME originating from re-summed ladder diagrams with arbitrary number of loops. We find that on-site repulsive interactions affect the chiral magnetic conductivity almost exclusively through the enhancement of the renormalized chiral chemical potential. Our results suggest that nontrivial corrections to the chiral magnetic conductivity due to inter-fermion interactions are not relevant in practice, since they only become important when the CME response is strongly suppressed by the large gap in the energy spectrum.

preprint2014arXiv

A lattice mean-field study of the phase diagram of interacting parity-breaking Weyl semimetals

We perform a mean-field study of the phase diagram of interacting Weyl semimetals with broken parity, that is, with different densities of right- and left-handed quasiparticles. As a simple model system, we consider the Wilson-Dirac Hamiltonian with the chiral chemical potential and on-site repulsive interactions. We find that the chiral chemical potential somewhat shrinks the region of the pion condensation (Aoki phase) in the parameter space of the bare mass and the interaction strength, so that the condensation thresholds are at smaller interaction strengths. The renormalized chiral chemical potential monotonously grows with interaction strength everywhere in the phase diagram, and only the growth rate is discontinuous across the phase transition lines. These findings are in full agreement with previous results obtained by one of the authors for the continuum Dirac Hamiltonian, except for the fact that for our lattice model with explicitly broken chiral symmetry the boundaries of the Aoki phase remain sharp second-order phase transitions even at nonzero chiral chemical potential and there are no signatures of Cooper-type instabilities in the weakly interacting regime.

preprint2014arXiv

A Method to Calculate Conserved Currents and Fermionic Force for the Lanczos Approximation to the Overlap Dirac Operator

The overlap Dirac operator obeys the Ginsparg-Wilson equation and offers a possibility to introduce chiral symmetry on the lattice. Evaluating the overlap operator is numerically very expensive and one has to rely on approximation methods. At finite chemical potential the overlap operator can be efficiently computed with the two-sided Lanczos algorithm. To calculate conserved currents on the lattice, or to evaluate the fermionic force in HMC calculations, one needs to compute derivatives of the Dirac operator with respect to gauge fields. In this paper we present a method to simultaneously compute the action of the overlap operator and its derivative on a source vector.