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

Jacob Finkenrath

Jacob Finkenrath 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.NAquant-ph

Catalog footprint

What is connected

24works
5topics
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

24 published item(s)

preprint2024arXiv

Low-lying baryon masses using twisted mass fermions ensembles at the physical pion mass

We investigate the low-lying baryon spectrum using three $N_f$=2+1+1 ensembles simulated with physical values of the quark masses and lattice spacings of 0.080, 0.069, and 0.057 fm. The ensembles are generated using twisted mass clover-improved fermions and the Iwasaki gauge action. The spatial length is kept approximately the same at about 5.1 fm to 5.5 fm fulfilling the condition $m_πL$> 3.6. We investigate isospin splitting within isospin multiples and verify that for most cases the isospin splitting for these lattice spacing is consistent with zero. In the couple of cases, for which there is a non-zero value, in the continuum limit, the mass splitting goes to zero. The baryon masses are extrapolated to the continuum limit using the three $N_f$=2+1+1 ensembles and are compared to other recent lattice QCD results. For the strange and charm quark masses, we find, respectively, $m_s$(2 GeV)=99.2(2.7) MeV and $m_c$(3 GeV)=1.015(39) GeV. The values predicted for the masses of the doubly charmed $Ξ_{cc}^\star$, $Ω_{cc}$ and $Ω_{cc}^\star$ baryons are 3.676(55) GeV, 3.703(51) GeV and 3.803(50) GeV, respectively, and for the triply charmed $Ω_{ccc}$ baryon is 4.785(71) GeV.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Implementation of Simultaneous Inversion of a Multi-shifted Dirac Matrix for Twisted-Mass Fermions within DD-αAMG

At physical light quark masses, efficient linear solvers are crucial for carrying out the millions of inversions of the Dirac matrix required for obtaining high statistics in quark correlation functions. Adaptive algebraic multi-grid methods have proven to be very efficient in such cases, exhibiting mild critical slowing down towards very light quark masses and outperforming traditional solver methods, such as the conjugate gradient method, at the physical point. We will discuss our implementations of simultaneous inversion of a (degenerate) Dirac matrix for twisted-mass fermions for multiple right-hand-sides (rhs) with multi-shifts and block-Krylov solvers. The implementation is carried out within the community library DD$α$AMG, which implements aggregation-based Domain Decomposition adaptive algebraic multi-grid methods. The block-Krylov solvers are provided via the Fast Accurate Block Linear krylOv Solver (Fabulous) library and can be used at coarser levels. Our code inverts Dirac matrices with different twisted-mass terms and for multiple rhs simultaneously and is thus also suitable for components within a typical lattice QCD simulation workflow, such as the rational approximation. We show preliminary results on scalability and compare the performance of our implementation when using different Block-Krylov solver techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lyncs-API: a Python API for Lattice QCD applications

We present Lyncs-API, a Python API for Lattice QCD applications currently under development. Lyncs aims to bring several widely used libraries for Lattice QCD under a common framework. Lyncs flexibly links to libraries for CPUs and GPUs in a way that can accommodate additional computing architectures as these arise, achieving performance-portability for the calculations while maintaining the same high-level workflow. Lyncs distributes calculations using Dask and mpi4py, with bindings to the libraries automatically generated by cppyy. While Lyncs is designed to allow linking to multiple libraries, we focus on a set of targeted packages that include DDalphaAMG, tmLQCD, QUDA and c-lime. More libraries will be added in the future. We also develop generic-purpose tools for facilitating the usage of Python in Lattice QCD and HPC in general. The project is open-source, community-oriented and available on Github.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Pion-pole contribution to HLbL from twisted mass lattice QCD at the physical point

We report on our computation of the pion transition form factor ${\cal F}_{P\rightarrow γ^*γ^*}$ from twisted mass lattice QCD in order to determine the numerically dominant light pseudoscalar pole contribution in the hadronic light-by-light scattering contribution to the anomalous magnetic moment of the muon $a_μ=(g-2)_μ$. The pion transition form factor is computed directly at the physical point. We present first results for our estimate of the pion-pole contribution with kinematic setup for the pion at rest.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Tackling critical slowing down using global correction steps with equivariant flows: the case of the Schwinger model

We propose a new method for simulating lattice gauge theories in the presence of fermions. The method combines flow-based generative models for local gauge field updates and hierarchical updates of the factorized fermion determinant. The flow-based generative models are restricted to proposing updates to gauge-fields within subdomains, thus keeping training times moderate while increasing the global volume. We apply our method performs to the 2-dimensional (2D) Schwinger model with $N_f=2$ Wilson Dirac fermions and show that no critical slowing down is observed in the sampling of topological sectors up to $β=8.45$. Furthermore, we show that fluctuations can be suppressed exponentially with the distance between active subdomains, allowing us to achieve acceptance rates of up to $99\%$ for the outer-most accept/reject step on lattices volumes of up to $V=128\times128$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Twisted mass gauge ensembles at physical values of the light, strange and charm quark masses

Lattice QCD simulations directly at physical masses of dynamical light, strange and charm quarks are highly desirable especially to remove systematic errors due to chiral extrapolations. However such simulations are still challenging. We discuss the adaption of efficient algorithms, like multi-grid methods or higher order integrators, within the molecular dynamic steps of the Hybrid Monte Carlo algorithm, that are enabling simulations of a new set of gauge ensembles by the Extended Twisted Mass collaboration (ETMC). We present the status of the on-going ETMC simulation effort that aims to enabling studies of finite size and discretization effects. We work within the twisted mass discretization which is free of odd-discretization effects at maximal twist and present our tuning procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Ruling Out the Massless Up-Quark Solution to the Strong CP Problem by Computing the Topological Mass Contribution with Lattice QCD

The infamous strong CP problem in particle physics can in principle be solved by a massless up quark. In particular, it was hypothesized that topological effects could substantially contribute to the observed nonzero up-quark mass without reintroducing CP violation. Alternatively to previous work using fits to chiral perturbation theory, in this Letter, we bound the strength of the topological mass contribution with direct lattice QCD simulations, by computing the dependence of the pion mass on the dynamical strange-quark mass. We find that the size of the topological mass contribution is inconsistent with the massless up-quark solution to the strong CP problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quark masses and decay constants in $N_f=2+1+1$ isoQCD with Wilson clover twisted mass fermions

We present a preliminary study of the pion, kaon and D-meson masses and decay constants in isosymmetric QCD, as well as a preliminary result for the light-quark renormalized mass. The analysis is based on the gauge ensembles produced by ETMC with $N_f=2+1+1$ flavours of Wilson-clover twisted mass quarks, spanning a range of lattice spacings from $\sim0.10$ to $0.07$ fm and include configurations at the physical pion point on lattices with linear size up to $L~\sim~5.6$~fm

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tetraquark interpolating fields in a lattice QCD investigation of the $D_{s0}^\ast(2317)$ meson

We investigate the $D_{s0}^\ast(2317)$ meson using lattice QCD and considering correlation functions of several $\bar{c} s$ two-quark and $\bar{c} s (\bar{u} u + \bar{d} d)$ four-quark interpolating fields. These interpolating fields generate different structures in color, spin and position space including quark-antiquark pairs, tetraquarks and two-meson scattering states. For our computation we use an ensemble simulated with pion mass $m_π\approx 0.296 \, \textrm{GeV}$ and spatial volume of extent $2.90 \, \textrm{fm}$. We find in addition to the expected spectrum of two-meson scattering states another state around $60 \, \textrm{MeV}$ below the $D K$ threshold, which we interpret as the $D_{s0}^\ast(2317)$ meson. This state couples predominantly to a quark-antiquark interpolating field and only weakly to a $D K$ two-meson interpolating field. The coupling to the tetraquark interpolating fields is essentially zero, rendering a tetraquark interpretation of the $D_{s0}^\ast(2317)$ meson rather unlikely. Moreover, we perform a scattering analysis using Lüscher's method and the effective range approximation to determine the $D_{s0}^\ast(2317)$ mass for infinite spatial volume. We find this mass $51 \, \textrm{MeV}$ below the $D K$ threshold, rather close to both our finite volume result and the experimentally observed value.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adapted nested force-gradient integrators: the Schwinger model case

We study a novel class of numerical integrators, the adapted nested force-gradient schemes, used within the molecular dynamics step of the Hybrid Monte Carlo (HMC) algorithm. We test these methods in the Schwinger model on the lattice, a well known benchmark problem. We derive the analytical basis of nested force-gradient type methods and demonstrate the advantage of the proposed approach, namely reduced computational costs compared with other numerical integration schemes in HMC.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adaptive Aggregation-based Domain Decomposition Multigrid for Twisted Mass Fermions

The Adaptive Aggregation-based Domain Decomposition Multigrid method (arXiv:1303.1377) is extended for two degenerate flavors of twisted mass fermions. By fine-tuning the parameters we achieve a speed-up of the order of hundred times compared to the conjugate gradient algorithm for the physical value of the pion mass. A thorough analysis of the aggregation parameters is presented, which provides a novel insight into multigrid methods for lattice QCD independently of the fermion discretization.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

DDalphaAMG for Twisted Mass Fermions

We present the Adaptive Aggregation-based Domain Decomposition Multigrid method extended to the twisted mass fermion discretization action. We show comparisons of results as a function of tuning the parameters that enter the twisted mass version of the DDalphaAMG library (https://github.com/sbacchio/DDalphaAMG). Moreover, we linked the DDalphaAMG library to the tmLQCD software package and give details on the performance of the multigrid solver during HMC simulations at the physical point.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Disconnected diagrams with twisted-mass fermions

The latest results from the Twisted-Mass collaboration on disconnected diagrams at the physical value of the pion mass are presented. In particular, we focus on the sigma terms, the axial charges and the momentum fraction, all of them for the nucleon. A detailed error analysis for each observable follows, showing the strengths and weaknesses of the one-end trick. Alternatives are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Importance of closed quark loops for lattice QCD studies of tetraquarks

To investigate the light scalar tetraquark candidate $a_0(980)$ (quantum numbers $J^P =0^+$), a correlation matrix including a variety of two- and four-quark interpolating operators has to be computed. We discuss efficient techniques to compute the elements of this correlation matrix, in particular diagrams with closed quark loops. Furthermore, we present evidence that such diagrams are not negligible given our precision, and their contribution is essential to obtain physically meaningful results. In particular, we find indications of the existence of an "additional" state around the two-particle thresholds of $K+ \bar K$ and $η+π$, which could correspond to the $a_0(980)$ meson.

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

On the effects of heavy sea quarks at low energies

We present a factorisation formula for the dependence of light hadron masses and low energy hadronic scales on the mass $M$ of a heavy quark: apart from an overall factor $Q$, ratios such as $r_0(M)/r_0(0)$ are computable in perturbation theory at large $M$. The mass-independent factor $Q$ is obtained from the theory in the limit $M\to0$ and the decoupled theory with the heavy quark removed. The perturbation theory part is stable concerning different loop orders and our non-perturbative results match on quantitatively to the perturbative prediction. Upon taking ratios of different hadronic scales at the same mass, the perturbative function drops out and the ratios are given by the decoupled theory up to $M^{-2}$ corrections. Our present numerical results are obtained in a model calculation where there are no light quarks and a heavy doublet of quarks is decoupled. They are limited to masses a factor two below the charm. This is not large enough to see the $M^{-2}$ scaling predicted by the theory, but it is sufficient to verify - in the continuum limit - that the sea quark effects of quarks with masses around the charm mass are very small.

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