Source author record

Gunnar Bali

Gunnar Bali 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

19works
2topics
4close 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

19 published item(s)

preprint2023arXiv

Sigma terms of the baryon octet in $N_\mathrm{f} = 2+1$ QCD with Wilson quarks

A lot of progress has been made in the direct determination of nucleon sigma terms. Using similar methods, we consider the sigma terms of the other octet baryons as well. These are determined on CLS gauge field ensembles employing the Lüscher-Weisz gluon action and the Sheikholeslami-Wohlert fermion action with $N_\mathrm{f} = 2 + 1$. The ensembles have pion masses ranging from ${410}\,\mathrm{MeV}$ down to the physical value and lattice spacings covering a range between ${0.098}\,\mathrm{fm}$ and ${0.039}\,\mathrm{fm}$. We present some preliminary results for the pion and strange sigma terms and compare to indirect determinations. To do so, we discuss multi-state fits to tackle the well-known problem of excited state contamination comparing the ratio and summation methods also including priors.

preprint2022arXiv

Determination of $m_c$ from $N_f = 2+1$ QCD with Wilson fermions

We present preliminary results for the charm quark mass in the $N_f=4$ RGI scheme. These were obtained using $N_f=2+1$ CLS ensembles with $\mathcal{O}(a)$ non-perturbatively improved Wilson fermions. We employed five different lattice spacings, ranging down to $a\lesssim 0.04$ fm and realized approximately physical pion and kaon masses, with ensembles spread out along three different trajectories in the quark mass plane, enabling a thorough study of the dependence on the lattice spacing and the light and strange sea quark masses. We sketch our analysis strategy and find that the dominant errors at present are due to the renormalization and scale setting uncertainties.

preprint2016arXiv

Nucleon generalized form factors from lattice QCD with nearly physical quark masses

We determine generalized form factors of the nucleon from lattice simulations with $N_f = 2$ mass-degenerate non-perturbatively improved Wilson-Sheikholeslami-Wohlert fermions down to a pion mass of 150 MeV. We also present the resulting isovector quark angular momentum. Possible excited-state contaminations are investigated with correlated simultaneous fits.

preprint2015arXiv

(Approximate) Low-Mode Averaging with a new Multigrid Eigensolver

We present a multigrid based eigensolver for computing low-modes of the Hermitian Wilson Dirac operator. For the non-Hermitian case multigrid methods have already replaced conventional Krylov subspace solvers in many lattice QCD computations. Since the $γ_5$-preserving aggregation based interpolation used in our multigrid method is valid for both, the Hermitian and the non-Hermitian case, inversions of very ill-conditioned shifted systems with the Hermitian operator become feasible. This enables the use of multigrid within shift-and-invert type eigensolvers. We show numerical results from our MPI-C implementation of a Rayleigh quotient iteration with multigrid. For state-of-the-art lattice sizes and moderate numbers of desired low-modes we achieve speed-ups of an order of magnitude and more over PARPACK. We show results and develop strategies how to make use of our eigensolver for calculating disconnected contributions to hadronic quantities that are noisy and still computationally challenging. Here, we explore the possible benefits, using our eigensolver for low-mode averaging and related methods with high and low accuracy eigenvectors. We develop a low-mode averaging type method using only a few of the smallest eigenvectors with low accuracy. This allows us to avoid expensive exact eigensolves, still benefitting from reduced statistical errors.

preprint2015arXiv

Charmed baryon spectroscopy and light flavour symmetry from lattice QCD

We determine the ground state and first excited state masses of singly and doubly charmed spin 1/2 and 3/2 baryons with positive and negative parity. Configurations with $N_f=2+1$ non-perturbatively improved Wilson-clover fermions were employed, with the same quark action also being used for the valence quarks, including the charm. The spectrum is calculated for pion masses in the range $M_π\sim 259-460$ MeV at a lattice spacing $a\sim 0.075$ fm. Finite volume effects are studied comparing lattices with two different linear spatial extents ($1.8\,{\rm fm}$ and $2.4\,{\rm fm}$). The physical point is approached from the SU(3) limit keeping the flavour averaged light quark mass fixed. The baryon masses are extrapolated using expansions in the strange-light quark mass difference. Most particles fall into the expected SU(3) multiplets with well constrained extrapolations, the exceptions having a possibly more complex internal structure. Overall agreement is found with experiment for the masses and splittings of the singly charmed baryons. As part of the calculation an analysis of the lower lying charmonium, $D$ and $D_s$ spectra was performed in order to assess discretisation errors. The gross spectra are reproduced, including the $D^*_{s0}$, $D_{s1}$ and $D_1$ mesons, while at this single lattice spacing hyperfine splittings come out $10-20$ MeV too low.

preprint2015arXiv

Hadronic vacuum polarization and muon g-2 from magnetic susceptibilities on the lattice

We present and test a new method to compute the hadronic vacuum polarization function in lattice simulations. This can then be used, e.g., to determine the leading hadronic contribution to the anomalous magnetic moment of the muon. The method is based on computing susceptibilities with respect to external electromagnetic plane wave fields and allows for a precision determination of both the connected and the disconnected contributions to the vacuum polarization. We demonstrate that the statistical errors obtained with our method are much smaller than those quoted in previous lattice studies, primarily due to a very effective suppression of the errors of the disconnected terms. These turn out to vanish within small errors, enabling us to quote an upper limit. We also comment on the accuracy of the vacuum polarization function determined from present experimental R-ratio data.

preprint2015arXiv

QCD spectroscopy and quark mass renormalisation in external magnetic fields with Wilson fermions

We study the change of the QCD spectrum of low-lying mesons in the presence of an external magnetic field using Wilson fermions in the quenched approximation. Motivated by qualitative differences observed in the spectra of overlap and Wilson fermions for large magnetic fields, we investigate the dependence of the additive quark mass renormalisation on the magnetic field. We provide evidence that the magnetic field changes the critical quark mass both in the free case and on our quenched ensemble. The associated change of the bare quark mass with the magnetic field affects the spectrum and is relevant for the magnetic field dependence of a number of related quantities. We derive Ward identities for lattice and continuum QCD+QED from which we can extract the current quark masses. We also report on a first test of the tuning of the quark masses with the magnetic field using the current quark masses, and show that this tuning resolves the qualitative discrepancy between the Wilson and overlap spectra.

preprint2015arXiv

Topology and glueballs in SU(7) Yang-Mills with open boundary conditions

It is well known that the topology of gauge configurations generated in a Markov Monte-Carlo chain freezes as the continuum limit is approached. The corresponding autocorrelation time increases exponentially with the inverse lattice spacing, affecting the ergodicity of the simulation. In SU(N) gauge theories for large N this problem sets in at much coarser lattice spacings than for N=3. This means that its systematics can be studied on lattices that are smaller in terms of the number of lattice sites. It has been shown that using open boundary conditions in time allows instantons to be created and destroyed, restoring topological mobility and ergodicity. However, with open boundary conditions translational invariance is lost and the influence of spurious states propagating from the boundary into the bulk on physical correlators needs to be carefully evaluated. Moreover, while the total topological charge can be changed, the mobility of instantons across the lattice is still reduced. We consider SU(7) Yang-Mills theory and analyse its topological content in the periodic and open boundary condition cases. We also investigate scalar and pseudo-scalar glueball correlation functions.

preprint2011arXiv

Potentials between pairs of static-light mesons

We give an update on our ongoing investigations of potentials between pairs of static-light mesons in Nf=2 Lattice QCD, in different spin and isospin channels. The question of attraction and repulsion is particularly interesting with respect to the X(3872) charmonium state and charged candidates such as the Z+(4430). We employ the nonperturbatively improved Sheikholeslami-Wohlert fermion and the Wilson gauge actions at two lattice spacings a approx. 0.084 fm and a approx. 0.077 fm with a pseudoscalar mass of mPS approx. 770 MeV and mPS approx. 400 MeV respectively. We use stochastic all-to-all propagator techniques, improved by a hopping parameter expansion. The analysis is based on the variational method, utilizing various source and sink interpolators.

preprint2010arXiv

Hunting the static energy renormalon

We employ Numerical Stochastic Perturbation Theory (NSPT) together with twisted boundary conditions (TBC) to search for the leading renormalon in the perturbative expansion of the static energy. This renormalon is expected to emerge four times faster than the one for the gluon conden- sate in the plaquette. We extract the static energy from Polyakov loop calculations up to 12 loops and present preliminary results, indicating a significant step towards confirming the theoretical expectation.

preprint2010arXiv

Improved Semileptonic Form Factor Calculations in Lattice QCD

We investigate the computational efficiency of two stochastic based alternatives to the Sequential Propagator Method used in Lattice QCD calculations of heavy-light semileptonic form factors. In the first method, we replace the sequential propagator, which couples the calculation of two of the three propagators required for the calculation, with a stochastic propagator so that the calculations of all three propagators are independent. This method is more flexible than the Sequential Propagator Method but introduces stochastic noise. We study the noise to determine when this method becomes competitive with the Sequential Propagator Method, and find that for any practical calculation it is competitive with or superior to the Sequential Propagator Method. We also examine a second stochastic method, the so-called ``one-end trick", concluding it is relatively inefficient in this context. The investigation is carried out on two gauge field ensembles, using the non-perturbatively improved Wilson-Sheikholeslami-Wohlert action with N_f=2 mass-degenerate sea quarks. The two ensembles have similar lattice spacings but different sea quark masses. We use the first stochastic method to extract ${\mathcal O}(a)$-improved, matched lattice results for the semileptonic form factors on the ensemble with lighter sea quarks, extracting f_+(0).

preprint2010arXiv

Meson and baryon masses with low mode averaging

We describe and test a method known in the literature as low mode averaging to improve Euclidean two-point functions in lattice QCD using the low-lying eigenmodes of the Wilson-Dirac operator D. The contribution from the low modes is averaged over all positions of the quark sources while the contribution from high modes is calculated in the traditional way using one source point per lattice. We apply this method to different baryon and meson two-point functions and we compare the improvements using the eigenmodes of the non-hermitian operator D and the eigenmodes of the hermitian operator Q=gamma5 D. The convergence strongly depends on the parity of the states.

preprint2010arXiv

Static-light meson-meson potentials

We investigate potentials between pairs of static-light mesons in Nf=2 Lattice QCD, in different spin channels. The question of attraction and repulsion is particularly interesting with respect to the X(3872) charmonium state and charged candidates such as the Z+(4430). We employ the nonperturbatively improved Sheikholeslami-Wohlert fermion and the Wilson gauge actions at a lattice spacing a approx. 0.084 fm and a pseudoscalar mass mPS approx. 760 MeV. We use stochastic all-to-all propagator techniques, improved by a hopping parameter expansion. The analysis is based on the variational method, utilizing various source and sink interpolators.

preprint2009arXiv

$η'$-$η_c$-mixing with improved stochastic estimators

Charmonia are flavour singlet mesons and thus in principle contributions from disconnected quark line diagrams might affect their masses, either directly or via mixing with other flavour singlet channels. We present a first study that takes both effects into account. We employ improved stochastic all-to-all propagator techniques (including new methods) to calculate the diagrams that appear within the mixing matrix between the $η'$ and the $η_c$. The runs are initially performed on $N_f=2$ $16^3\times 32$ configurations with the non-perturbatively improved Sheikholeslami-Wilson action, both for valence and sea quarks.

preprint2009arXiv

A Stochastic Method for Semileptonic Form Factor Calculations on the Lattice

We investigate an alternative to the Sequential Propagator Method used in Lattice QCD calculations of semileptonic form factors. We replace the sequential propagator with a stochastic propagator so that, in principle, all momentum and sink smearing combinations are available with only a single spin-color inversion. Practically, the stochastic noise is significant and must be reduced at the cost of more inversions. We study the behavior of the stochastic noise and compare the computational costs of this stochastic technique and the Sequential Propagator Method. We also present preliminary semileptonic form factor results using the stochastic technique on N_f=2 configurations with a non-perturbatively improved Sheikoleslami-Wohlert action generated by the QCDSF collaboration. At a fixed cost, measured in terms of the number of heavy-quark inversions, the method provides more correlators for the extraction of the form factors at various q^2's than the Sequential Propagator Method. These additional correlators reduce the total statistical errors of certain kinematic points, although the stochastic error is still comparable to the gauge error at other points.

preprint2009arXiv

Mixing of S-Wave Charmonia with $\mathrm{D}\overline{\mathrm{D}}$ Molecule States

Charmonium states can decay into pairs of $D$ and $\overline{D}$ mesons if their masses are above the allowed decay thresholds. In general $c\bar{c}$ states near threshold will also undergo mixing with $D\overline{D}$ molecular (or tetraquark) states, by creation and annihilation of light quark-antiquark pairs. The investigation of such effects sheds light on the higher Fock state contributions to charmonium wavefunctions and on mass shifts, relative to a scenario where such mixing effects are neglected. A variational approach is applied to a mixing matrix between operators of both sectors, of $c\bar{c}$ and of $D\overline{D}$ molecular type. The efficient calculation of several diagrams appearing in this matrix requires all-to-all propagators, which are realized by sophisticated stochastic estimator techniques. The runs are performed on $n_F=2$ $24^3\times 48$ lattice volumes with $m_π \approx 280$ MeV, using the non-perturbatively improved clover Wilson action, both for valence and for sea quarks.

preprint2009arXiv

Static-static-light baryonic potentials

We determine doubly heavy baryonic potentials as a function of the distance between the two static sources, coupled to a light relativistic quark, for different quantum numbers. We use the variational method to compute the ground state and the first two excitations. These can be used as an input to nonrelativistic models or to NRQCD calculations of properties of doubly heavy baryons. We compare our findings with a factorization model. We employ all-to-all propagator methods, improved by an additional hopping parameter expansion and Wuppertal smearing on QCDSF configurations with two sea quark flavors.

preprint2009arXiv

Strangeness and charm content of the nucleon

We present results on the scalar strangeness and charm contents of the nucleon and of the disconnected contributions to the nucleon spin. These are obtained on two flavour non-perturbatively improved Sheikholeslami-Wilson configurations at a pseudoscalar mass of 290 MeV. We quote f_{T_s} = m_s<N|ss|N>/m_N< = 0.070(22) as our preliminary value for the strange quark fraction of the nucleon mass and a modified minimal subtraction scheme value Delta s = -0.015(10), with as yet unknown systematics, for the strangeness contribution to the spin.