Researcher profile

Jacob Finkenrath

Jacob Finkenrath contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

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

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.

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.

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.

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

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.

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.

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

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.

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.

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.