BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Björn Leder

Björn Leder appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-phhep-thmath.NA

Catalog footprint

What is connected

19works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

19 published item(s)

preprint2020arXiv

New linearization and reweighting for simulations of string sigma-model on the lattice

We study the discretized worldsheet of Type IIB strings in the Gubser-Klebanov-Polyakov background in a new setup, which eliminates a complex phase previously detected in the fermionic determinant. A sign ambiguity remains, which a study of the fermionic spectrum shows to be related to Yukawa-like terms, including those present in the original Lagrangian before the linearization standard in a lattice QFT approach. Monte Carlo simulations are performed in a large region of the parameter space, where the sign problem starts becoming severe and instabilities appear due to the zero eigenvalues of the fermionic operator. To face these problems, simulations are conducted using the absolute value of a fermionic Pfaffian obtained introducing a small twisted-mass term, acting as an infrared regulator, into the action. The sign of the Pfaffian and the low modes of the quadratic fermionic operator are then taken into account by a reweighting procedure of which we discuss the impact on the measurement of the observables. In this setup we study bosonic and fermionic correlators and observe a divergence in the latter, which we argue - also via a one-loop analysis in lattice perturbation theory - to originate from the U(1)-breaking of our Wilson-like discretization for the fermionic sector.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Green-Schwarz superstring on the lattice

We consider possible discretizations for a gauge-fixed Green-Schwarz action of Type IIB superstring. We use them for measuring the action, from which we extract the cusp anomalous dimension of planar $\mathcal{N}=4$ SYM as derived from AdS/CFT, as well as the mass of the two $AdS$ excitations transverse to the relevant null cusp classical string solution. We perform lattice simulations employing a Rational Hybrid Monte Carlo (RHMC) algorithm and two Wilson-like fermion discretizations, one of which preserves the global $SO(6)$ symmetry of the model. We compare our results with the expected behavior at various values of $g=\frac{\sqrtλ}{4π}$. For both the observables, we find a good agreement for large $g$, which is the perturbative regime of the sigma-model. For smaller values of $g$, the expectation value of the action exhibits a deviation compatible with the presence of quadratic divergences. After their non-perturbative subtraction the continuum limit can be taken, and suggests a qualitative agreement with the non-perturbative expectation from AdS/CFT. Furthermore, we detect a phase in the fermion determinant, whose origin we explain, that for very small $g$ leads to a sign problem not treatable via standard reweigthing. The continuum extrapolations of the observables in the two different discretizations agree within errors, which is strongly suggesting that they lead to the same continuum limit. Part of the results discussed here were presented earlier in arXiv:1601.04670.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Lattice and string worldsheet in AdS/CFT: a numerical study

We consider a possible discretization for the gauge-fixed Green-Schwarz (two-dimensional) sigma-model action for the Type IIB superstring and use it for measuring the cusp anomalous dimension of planar $\mathcal{N}=4$ SYM as derived from string theory. We perform lattice simulations employing a Rational Hybrid Monte Carlo (RHMC) algorithm and a Wilson-like fermion discretization. In this preliminary study, we compare our results with the expected behavior for very large values of $g=\frac{\sqrtλ}{4π}$, which is the perturbative regime of the sigma-model, and find a qualitative agreement at finite lattice spacing. For smaller $g$ the continuum limit is obstructed by a divergence. We also detect a phase in the fermion determinant, whose origin we explain, which for small $g$ leads to a sign problem not treatable via standard reweigthing. Results presented here are discussed thoroughly in~\cite{toappear}.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Isospin Effects by Mass Reweighting

Most of today's lattice simulations are performed in the isospin symmetric limit of the light quark sector. Mass reweighting is a technique to include effects of isospin breaking in the sea quarks at moderate numerical cost. We will give a summary of our recent results on fine lattices with light quark masses and will show how light quark masses can be extracted by introducing suitable tuning conditions for the bare mass parameters. In general the reweighting factor introduces additional fluctuations and thus increases the statistical uncertainties. In the case of isospin reweighting this factor is a ratio of fermion determinants. The stochastic evaluation of the determinants potentially leads to stochastic noise in observables. We show the quark mass and the volume dependence of these fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Perturbative versus non-perturbative decoupling of heavy quarks

We simulate a theory with $N_f=2$ heavy quarks of mass $M$. At energies much smaller than $M$ the heavy quarks decouple and the theory can be described by an effective theory which is a pure gauge theory to leading order in $1/M$. We present results for the mass dependence of ratios such as $t_0(M)/t_0(0)$. We compute these ratios from simulations and compare them to the perturbative prediction. The latter relies on a factorisation formula for the ratios which is valid to leading order in $1/M$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tuning of the strange quark mass with optimal reweighting

Quark mass reweighting can be used to tune the mass of dynamical quarks. The basic idea is to use gauge field ensembles generated at some bare mass parameters to evaluate observables at different bare sea quark masses. This involves the computation of so called reweighing factors which are given as ratios of fermion determinants. In the case of simulations including the strange quark, reweighting can be used to improve the approach towards physical quark masses. Optimal reweighting strategies combine a change of the strange quark mass with a change of the light quark masses in order to minimize the fluctuations of the reweighting factor. We present numerical test of such strategies for recent CLS2 simulations and a software package for mass reweighting based on openQCD.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Adaptive Aggregation Based Domain Decomposition Multigrid for the Lattice Wilson Dirac Operator

In lattice QCD computations a substantial amount of work is spent in solving discretized versions of the Dirac equation. Conventional Krylov solvers show critical slowing down for large system sizes and physically interesting parameter regions. We present a domain decomposition adaptive algebraic multigrid method used as a precondtioner to solve the "clover improved" Wilson discretization of the Dirac equation. This approach combines and improves two approaches, namely domain decomposition and adaptive algebraic multigrid, that have been used seperately in lattice QCD before. We show in extensive numerical test conducted with a parallel production code implementation that considerable speed-up over conventional Krylov subspace methods, domain decomposition methods and other hierarchical approaches for realistic system sizes can be achieved.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multigrid Preconditioning for the Overlap Operator in Lattice QCD

The overlap operator is a lattice discretization of the Dirac operator of quantum chromodynamics, the fundamental physical theory of the strong interaction between the quarks. As opposed to other discretizations it preserves the important physical property of chiral symmetry, at the expense of requiring much more effort when solving systems with this operator. We present a preconditioning technique based on another lattice discretization, the Wilson-Dirac operator. The mathematical analysis precisely describes the effect of this preconditioning in the case that the Wilson-Dirac operator is normal. Although this is not exactly the case in realistic settings, we show that current smearing techniques indeed drive the Wilson-Dirac operator towards normality, thus providing a motivation why our preconditioner works well in computational practice. Results of numerical experiments in physically relevant settings show that our preconditioning yields accelerations of up to one order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

One flavor mass reweighting in lattice QCD

One flavor mass reweighting can be used in lattice QCD computations to fine tune the quark masses to their physical values. We present a new method that utilizes an unbiased stochastic estimation of the one flavor determinant. The stochastic estimation is based on the integral representation of the determinant of a complex matrix, which we prove. In contrast to other methods it can also be applied in situations where the determinant has a complex phase. The stochastic error is controlled by determinant factorizations based on mass interpolation and Schur decomposition. As an example of an application we demonstrate how the method can be used to tune the up-down quark mass difference.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

One flavor mass reweighting: foundations

Reweighting is not a new method in lattice QCD, but a comprehensive analysis is missing in the literature. We close this gap by presenting: (i) a proof of an integral representation of the complex determinant of a complex matrix, (ii) a method to control the stochastic error of its Monte Carlo estimation, (iii) expansions of the stochastic error and the ensemble fluctuations of the one flavor reweighting factor. Based on (iii) we present a detailed scaling analysis and optimized reweighting strategies. As an application we analyze the ensemble fluctuations of the reweighting factor corresponding to the sea contribution to isospin splitting and predict at physical quark masses a standard deviation of +-20%.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Physical and cut-off effects of heavy sea quarks

We simulate a theory with two dynamical O($a$) improved Wilson quarks whose mass $M$ ranges from a factor eight up to a factor two below the charm quark mass and at three values of the lattice spacing ranging from 0.066 to 0.034 fm. This theory is a prototype to study the decoupling of heavy quarks. We measure the mass and cut-off dependence of ratios of gluonic observables defined from the Wilson flow or the static potential. The size of the 1/$M$ corrections can be determined and disentangled from the lattice artifacts. The difference with the pure gauge theory is at the percent level when two quarks with a mass of the charm quark are present.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Application of Domain Decomposition to the Evaluation of Fermion Determinant Ratios

We analyze the fluctuations in the case of mass reweighting for N_f=2 Wilson fermions. We use a domain decomposition factorization of the fermion determinant. Ratios of determinants are estimated stochastically. We study the stochastic and the ensemble fluctuations as a function of the volume V and the mass shift Delta m. With our result it is possible to estimate the cost and the effectiveness of mass reweighting. In addition we introduce a stochastic estimation for the one flavor case without using the square root.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fermions as Global Correction: the QCD Case

It is widely believed that the fermion determinant cannot be treated in global acceptance-rejection steps of gauge link configurations that differ in a large fraction of the links. However, for exact factorizations of the determinant that separate the ultraviolet from the infrared modes of the Dirac operator it is known that the latter show less variation under changes of the gauge field compared to the former. Using a factorization based on recursive domain decomposition allows for a hierarchical algorithm that starts with pure gauge updates of the links within the domains and ends after a number of filters with a global acceptance-rejection step. Ratios of determinants have to be treated stochastically and we construct techniques to reduce the noise. We find that the global acceptance rate is high on moderate lattice sizes and demonstrate the effectiveness of the hierarchical filter.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Aggregation-based Multilevel Methods for Lattice QCD

In Lattice QCD computations a substantial amount of work is spent in solving the Dirac equation. In the recent past it has been observed that conventional Krylov solvers tend to critically slow down for large lattices and small quark masses. We present a Schwarz alternating procedure (SAP) multilevel method as a solver for the Clover improved Wilson discretization of the Dirac equation. This approach combines two components (SAP and algebraic multigrid) that have separately been used in lattice QCD before. In combination with a bootstrap setup procedure we show that considerable speed-up over conventional Krylov subspace methods for realistic configurations can be achieved.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The shape of the static potential with dynamical fermions

We present the analysis of the static potential extracted from Wilson loops measured on CLS ensembles generated with Wilson gauge action and Nf=2 flavors of O(a) improved Wilson quarks at three different lattice spacings and a range of quark masses. The shape of the static potential at distances well below the string breaking region is studied in terms of renormalized couplings derived from the static force and its derivative. We comment on the (im)possibility of extracting the Lambda parameter at our smallest lattice spacing a=0.05 fm. Finally we give an update on the scale determination through r0.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The strange quark mass and Lambda parameter of two flavor QCD

We complete the non-perturbative calculations of the strange quark mass and the Lambda parameter in two flavor QCD by the ALPHA collaboration. The missing lattice scale is determined via the kaon decay constant, for whose chiral extrapolation complementary strategies are compared. We also give a value for the scale r_0 in physical units as well as an improved determination of the renormalization constant Z_A.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Determination of the Static Potential with Dynamical Fermions

We present in detail a technique to extract the potential between a static quark and anti-quark pair from Wilson loops measured on dynamical configurations. This technique is based on HYP smearing and leads to an exponential improvement of the noise-to-signal ratio of Wilson loops. We explain why the correct continuum potential is obtained and show numerical evidence that the cut-off effects are small. We present precise results for the non-perturbative potential. As applications, we determine the scale r_0/a and study the shape of the static potential in the range of distances around r_0, where it can be compared with phenomenological potential models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fermions as global correction in lattice QCD

The fermion determinant is a highly non-local object and its logarithm is an extensive quantity. For these reasons it is widely believed that the determinant cannot be treated in acceptance steps of gauge link configurations that differ in a large fraction of the links. However, for exact factorisations of the determinant that separate the ultraviolet from the infra-red modes of the Dirac operator it is known that the latter show less variation under changes of the gauge field compared to the former. Using a factorisation based on recursive domain decomposition allows for a hierarchical algorithm that starts with pure gauge updates of the links within the domains and ends after a number of filters with a global acceptance step. We find that the global acceptance rate is high on moderate lattice sizes. Whether this type of algorithm can help in curing the problem of critical slowing down is presently under study.

0 citations0 reviewsNo code linkedOpen access