Source author record

Savvas Zafeiropoulos

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

23works
9topics
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

23 published item(s)

preprint2022arXiv

Extraction of $B_s \to D^{(*)}_s$ form factors from ${ N_f}=2$ lattice QCD

We report on a two-flavour lattice QCD study of the $B_s \to D_s$ and $B_s \to D^*_s$ transitions parameterized, in the heavy quark limit, by the form factors ${\cal G}$, and $h_{A_1}$, $h_{A_2}$ and $h_{A_3}$, respectively. In the search for New Physics through tests of lepton flavour universality, $B_s$ decay channels are complementary to the $B$ decays widely studied at $B$ factories and LHCb, while on the theory side they can be better controlled than the $B_c$ and $Λ_b$ decays. The purpose of this exploratory two-flavour study is, in preparation for future analyses of lattice QCD simulations with ${N_f>2}$ and physical quark-masses, to gain experience on a suitable method for a lattice extraction of form factors associated with $b \to c$ currents that may yield tighter control over systematic effects like contamination from excited states and cut-off effects. We obtain the zero-recoil values ${\cal G}^{B_s \to D_s}(1)=1.03(14)$ and $h^{B_s \to D^*_s}(1)=0.85(16)$.

preprint2022arXiv

Properties, ensembles and hadron spectra with Stabilised Wilson Fermions

In this joint contribution we announce the formation of the "OPEN LATtice initiative", https://openlat1.gitlab.io, to study Stabilised Wilson Fermions (SWF). They are a new avenue for QCD calculations with Wilson-type fermions and we report results on our continued study of this framework: Tuning the clover improvement coefficient, and extending the reach of lattice spacings to $a=0.12$ fm. We fix the flavor symmetric points $m_π=m_K=412$ MeV at $a=0.055,0.064, 0.077, 0.094, 0.12$ fm and define the trajectories to the physical point by fixing the trace of the quark mass matrix. Currently our pion mass range extends down to $m_π\sim200$ MeV. We outline our tuning goals and strategy as well as our future planned ensembles. First scaling studies are performed on $f_π$ and $m_π$. Additionally results of a preliminary continuum extrapolation of $m_N$ at the flavor symmetric point are presented. Going further a first determination of the light and strange hadron spectrum chiral dependence is shown, which serves to check the quality of the action for precision measurements. We also investigate other quantities such as flowed gauge observables to study how the continuum limit is approached. Taken together we observe the SWF enable us to perform stable lattice simulations across a large range of parameters in mass, volume and lattice spacing. Pooling resources our new initiative has made our reported progress possible and through it we will share generated gauge ensembles under an open science philosophy.

preprint2022arXiv

Reconstructing QCD Spectral Functions with Gaussian Processes

We reconstruct ghost and gluon spectral functions in 2+1 flavor QCD with Gaussian process regression. This framework allows us to largely suppress spurious oscillations and other common reconstruction artifacts by specifying generic magnitude and length scale parameters in the kernel function. The Euclidean propagator data are taken from lattice simulations with domain wall fermions at the physical point. For the infrared and ultraviolet extensions of the lattice propagators as well as the low-frequency asymptotics of the ghost spectral function, we utilize results from functional computations in Yang-Mills theory and QCD. This further reduces the systematic error significantly. Our numerical results are compared against a direct real-time functional computation of the ghost and an earlier reconstruction of the gluon in Yang-Mills theory. The systematic approach presented in this work offers a promising route towards unveiling real-time properties of QCD.

preprint2021arXiv

The transversity parton distribution function of the nucleon using the pseudo-distribution approach

We present a determination of the non-singlet transversity parton distribution function (PDF) of the nucleon, normalized with respect to the tensor charge at $μ^2=2$ GeV$^2$ from lattice quantum chromodynamics. We apply the pseudo-distribution approach, using a gauge ensemble with a lattice spacing of 0.094 fm and the light quark mass tuned to a pion mass of 358 MeV. We extract the transversity PDF from the analysis of the short-distance behavior of the Ioffe-time pseudo-distribution using the leading-twist next-to-leading order (NLO) matching coefficients calculated for transversity. We reconstruct the $x$-dependence of the transversity PDF through an expansion in a basis of Jacobi polynomials in order to reduce the PDF ansatz dependence. Within the limitations imposed by a heavier-than-physical pion mass and a fixed lattice spacing, we present a comparison of our estimate for the valence transversity PDF with the recent global fit results based on single transverse spin asymmetry. We find the intrinsic nucleon sea to be isospin symmetric with respect to transversity.

preprint2020arXiv

Lefschetz thimbles decomposition for the Hubbard model on the hexagonal lattice

We propose a framework to study the properties of the Lefschetz thimbles decomposition for lattice fermion models approaching the thermodynamic limit. The proposed set of algorithms includes the Schur complement solver and the exact computation of the derivatives of the fermion determinant. It allows us to solve the gradient flow (GF) equations taking into account the fermion determinant exactly, with high performance. We can find both real and complex saddle points and describe the structure of the Lefschetz thimbles decomposition for large enough lattices to extrapolate the results to the thermodynamic limit. The algorithms are described for a general lattice fermion model, with emphasis on two types of lattice discretizations for relativistic fermions (staggered and Wilson), as well as on interacting tight-binding models for condensed matter systems. We apply these algorithms to the Hubbard model on a hexagonal lattice, dealing with lattice volumes as large as 12x12 with 256 steps in Euclidean time, in order to capture the properties of the thimbles decomposition as the thermodynamic, low-temperature, and continuum limits are approached. The complexity of the thimbles decomposition appears to be dependent on the form of the Hubbard-Stratonovich (HS) transformation. In particular, the evidence is provided for the existence of an optimal regime for the hexagonal lattice Hubbard model, with a reduced number of thimbles being important in the overall sum. We have performed quantum Monte Carlo (QMC) simulations using GF to deform the integration contour into the complex plane and demonstrated the agreement with exact diagonalization on small volumes (8 sites in space) if optimal setup for HS transformation is used. The residual sign problem is compared with the state-of-the-art BSS-QMC. We show that the average sign can be kept substantially higher using the Lefschetz thimbles approach.

preprint2020arXiv

Parton Distribution Functions from Ioffe time pseudo-distributions

In this paper, we present a detailed study of the unpolarized nucleon parton distribution function (PDF) employing the approach of parton pseudo-distribution functions. We perform a systematic analysis using three lattice ensembles at two volumes, with lattice spacings $a=$ 0.127 fm and $a=$ 0.094 fm, for a pion mass of roughly 400 MeV. With two lattice spacings and two volumes, both continuum limit and infinite volume extrapolation systematic errors of the PDF are estimated. In addition to the $x$ dependence of the PDF, we compute their first two moments and compare them with the pertinent phenomenological determinations.

preprint2016arXiv

A real-time lattice simulation of the thermalization of a gluon plasma: first results

To achieve an understanding of the thermalization of a quark-gluon plasma, starting from QCD without using model assumptions, is a formidable task. We study the early stage dynamics of a relativistic heavy ion collision in the framework of real time simulations of classical Yang-Mills theory in a static box with the color glass condensate as initial condition. Our study generalizes a previous one by Fukushima and Gelis from SU(2) to the realistic case of SU(3). We calculate the chromo-electric and chromo-magnetic energy densities as well as the ratio of longitudinal and transverse pressure as a function of time as probes for thermalization. Our preliminary results on coarse lattices show the occurrence of Weibel instabilities prior to thermalization.

preprint2016arXiv

Magnetic Catalysis in Graphene Effective Field Theory

We report on the first observation of magnetic catalysis at zero temperature in a fully nonperturbative simulation of the graphene effective field theory. Using lattice gauge theory, a nonperturbative analysis of the theory of strongly-interacting, massless, (2+1)-dimensional Dirac fermions in the presence of an external magnetic field is performed. We show that in the zero-temperature limit, a nonzero value for the chiral condensate is obtained which signals the spontaneous breaking of chiral symmetry. This result implies a nonzero value for the dynamical mass of the Dirac quasiparticle. This in turn has been posited to account for the quantum-Hall plateaus that are observed at large magnetic fields.

preprint2015arXiv

Dirac Spectrum of the Wilson Dirac Operator for QCD with Two Colors

We study the lattice artefacts of the Wilson Dirac operator for QCD with two colors and fermions in the fundamental representation from the viewpoint of chiral perturbation theory. These effects are studied with the help of the following spectral observables: the level density of the Hermitian Wilson Dirac operator, the distribution of chirality over the real eigenvalues, and the chiral condensate for the quenched as well as for the unquenched theory. We provide analytical expressions for all these quantities. Moreover we derive constraints for the level density of the real eigenvalues of the non-Hermitian Wilson Dirac operator and the number of additional real modes. The latter is a good measure for the strength of lattice artefacts. All computations are confirmed by Monte Carlo simulations of the corresponding random matrix theory which agrees with chiral perturbation theory of two color QCD with Wilson fermions.

preprint2015arXiv

Discretization effects in $N_c=2$ QCD and Random Matrix Theory

We summarize the analytical solution of the Chiral Perturbation Theory for the Hermitian Wilson Dirac operator of $N_c=2$ QCD with quarks in the fundamental representation. Results have been obtained for the quenched microscopic spectral density, the distribution of the chiralities over the real modes and the chiral condensate. The analytical results are compared with results from a Monte Carlo simulation of the corresponding Random Matrix Theory.

preprint2015arXiv

Magnetic Catalysis in Graphene

One of the most important developments in condensed matter physics in recent years has been the discovery and characterization of graphene. A two-dimensional layer of Carbon arranged in a hexagonal lattice, graphene exhibits many interesting electronic properties, most notably that the low energy excitations behave as massless Dirac fermions. These excitations interact strongly via the Coulomb interaction and thus non-perturbative methods are necessary. Using methods borrowed from lattice QCD, we study the graphene effective theory in the presence of an external magnetic field. Graphene, along with other $(2+1)$-dimensional field theories, has been predicted to undergo spontaneous breaking of flavor symmetry including the formation of a gap as a result of the external magnetic field. This phenomenon is known as magnetic catalysis. Our study investigates magnetic catalysis using a fully non-perturbative approach.

preprint2015arXiv

Phase Diagram of Dynamical Twisted Mass Wilson Fermions at Finite Isospin Chemical Potential

We consider the phase diagram of twisted mass Wilson fermions of two-flavor QCD in the parameter space of the quark mass, the isospin chemical potential, the twist angle and the lattice spacing. This work extends earlier studies in the continuum and those at zero chemical potential. We evaluate the phase diagram as well as the spectrum of the (pseudo-)Goldstone bosons using the chiral Lagrangian for twisted mass Wilson fermions at non-zero isospin chemical potential. The phases are obtained from a mean field analysis. At zero twist angle we find that already an infinitesimal isospin chemical potential destroys the Aoki phase. The reason is that in this phase we have massless Goldstone bosons with a non-zero isospin charge. At finite twist angle only two different phases are present, one phase which is continuously connected to the Bose condensed phase at non-zero chemical potential and another phase which is continuously connected to the normal phase. For either zero or maximal twist the phase diagram is more complicated as the saddle point equations allow for more solutions.

preprint2015arXiv

Renormalization of quark propagator, vertex functions and twist-2 operators from twisted-mass lattice QCD at $N_f$=4

We present a precise non-perturbative determination of the renormalization constants in the mass independent RI'-MOM scheme. The lattice implementation uses the Iwasaki gauge action and four degenerate dynamical twisted mass fermions. The gauge configurations are provided by the ETM Collaboration. Renormalization constants for scalar, pseudo-scalar, vector and axial operators, as well as the quark propagator renormalization, are computed at three different values of the lattice spacing, two volumes and several twisted mass parameters. The method we developed allows for a precise cross-check of the running, thanks to the particular proper treatment of hypercubic artifacts. Results for the twist-2 operator $O_{44}$ are also presented.

preprint2015arXiv

The microscopic Twisted Mass Dirac spectrum and the spectral determination of the LECs of Wilson $χ$-PT

We present the comparison of the analytical microscopic spectral density for lattice QCD with $N_{\rm f}=2$ twisted mass fermions with the one obtained on the lattice utilizing configurations produced by the ETM collaboration. We extract estimates for the chiral condensate as well as the low-energy constant $W_8$ of Wilson $χ$-PT by employing spectral information of the Wilson Dirac operator with fixed index at finite volume.

preprint2014arXiv

Dirac Spectra of 2-dimensional QCD-like theories

We analyze Dirac spectra of two-dimensional QCD like theories both in the continuum and on the lattice and classify them according to random matrix theories sharing the same global symmetries. The classification is different from QCD in four dimensions because the anti-unitary symmetries do not commute with $γ_5$. Therefore in a chiral basis, the number of degrees of freedom per matrix element are not given by the Dyson index. Our predictions are confirmed by Dirac spectra from quenched lattice simulations for QCD with two or three colors with quarks in the fundamental representation as well as in the adjoint representation. The universality class of the spectra depends on the parity of the number of lattice points in each direction. Our results show an agreement with random matrix theory that is qualitatively similar to the agreement found for QCD in four dimensions. We discuss the implications for the Mermin-Wagner-Coleman theorem and put our results in the context of two-dimensional disordered systems.

preprint2014arXiv

Renormalization constants for $N_{\rm f}=2+1+1$ twisted mass QCD

We summarize recent non-perturbative results obtained for the renormalization constants computed in the RI'-MOM scheme for $N_{\rm f}=2+1+1$ twisted mass QCD. Our implementation employs the Iwasaki gauge action and four dynamical degenerate twisted mass fermions. Renormalization constants for scalar, pseudo-scalar, vector and axial operators, as well as the quark propagator renormalization, are computed at three different values of the lattice spacing, two different volumes and several values of the twisted mass. Our method allows for a precise cross-check of the running, because of the particular proper treatment of the hypercubic artifacts. Preliminary results for twist-2 operators are also presented.

preprint2013arXiv

A classification of 2-dim Lattice Theory

A unified classification and analysis is presented of two dimensional Dirac operators of QCD-like theories in the continuum as well as in a naive lattice discretization. Thereby we consider the quenched theory in the strong coupling limit. We do not only consider the case of a lattice which has an even number of lattice sites in both directions and is thus equivalent to the case of staggered fermions. We also study lattices with one or both directions with an odd parity to understand the general mechanism of changing the universality class via a discretization. Furthermore we identify the corresponding random matrix ensembles sharing the global symmetries of these QCD-like theories. Despite the Mermin-Wagner-Coleman theorem we find good agreement of lattice data with our random matrix predictions.

preprint2013arXiv

Discretization Effects in the $ε$ Domain of QCD

At nonzero lattice spacing the QCD partition function with Wilson quarks undergoes either a second order phase transition to the Aoki phase for decreasing quark mass or shows a first order jump when the quark mass changes sign. We discuss these phase transitions in terms of Wilson Dirac spectra and show that the first order scenario can only occur in the presence of dynamical quarks while in the quenched case we can only have a transition to the Aoki phase. The exact microscopic spectral density of the non-Hermitian Wilson Dirac operator with dynamical quarks is discussed as well. We conclude with some remarks on discretization effects for the overlap Dirac operator.

preprint2013arXiv

Random Matrix Models for the Hermitian Wilson-Dirac operator of QCD-like theories

We introduce Random Matrix Models for the Hermitian Wilson-Dirac operator of QCD-like theories. We show that they are equivalent to the $ε$-limit of the chiral Lagrangian for Wilson chiral perturbation theory. Results are obtained for two-color QCD with quarks in the fundamental representation of the color group as well as any-color QCD with quarks in the adjoint representation. For $N_c=2$ we also have obtained the lattice spacing dependence of the quenched average spectral density for a fixed value of the index of the Dirac operator. Comparisons with direct numerical simulations of the random matrix ensemble are shown.

preprint2013arXiv

Spectral Properties of the Wilson Dirac Operator and random matrix theory

Random Matrix Theory has been successfully applied to lattice Quantum Chromodynamics. In particular, a great deal of progress has been made on the understanding, numerically as well as analytically, of the spectral properties of the Wilson Dirac operator. In this paper, we study the infra-red spectrum of the Wilson Dirac operator via Random Matrix Theory including the three leading order $a^2$ correction terms that appear in the corresponding chiral Lagrangian. A derivation of the joint probability density of the eigenvalues is presented. This result is used to calculate the density of the complex eigenvalues, the density of the real eigenvalues and the distribution of the chiralities over the real eigenvalues. A detailed discussion of these quantities shows how each low energy constant affects the spectrum. Especially we consider the limit of small and large (which is almost the mean field limit) lattice spacing. Comparisons with Monte Carlo simulations of the Random Matrix Theory show a perfect agreement with the analytical predictions. Furthermore we present some quantities which can be easily used for comparison of lattice data and the analytical results.

preprint2013arXiv

The Effect of the Low Energy Constants on the Spectral Properties of the Wilson Dirac Operator

We summarize recent analytical results obtained for lattice artifacts of the non-Hermitian Wilson Dirac operator. Hereby we discuss the effect of all three low energy constants. In particular we study the limit of small lattice spacing and also consider the regime of large lattice spacing which is closely related to the mean field limit. Thereby we extract simple relations between measurable quantities like the average number of additional real modes and the low energy constants. These relations may improve the fitting for the low energy constants.

preprint2012arXiv

Eigenvalue Density of the non-Hermitian Wilson Dirac Operator

We find the lattice spacing dependence of the eigenvalue density of the non-Hermitian Wilson Dirac operator in the $ε$-domain. The starting point is the joint probability density of the corresponding random matrix theory. In addition to the density of the complex eigenvalues we also obtain the density of the real eigenvalues separately for positive and negative chiralities as well as an explicit analytical expression for the number of additional real modes.

preprint2012arXiv

Random Matrix Models for Dirac Operators at finite Lattice Spacing

We study discretization effects of the Wilson and staggered Dirac operator with $N_{\rm c}>2$ using chiral random matrix theory (chRMT). We obtain analytical results for the joint probability density of Wilson-chRMT in terms of a determinantal expression over complex pairs of eigenvalues, and real eigenvalues corresponding to eigenvectors of positive or negative chirality as well as for the eigenvalue densities. The explicit dependence on the lattice spacing can be readily read off from our results which are compared to numerical simulations of Wilson-chRMT. For the staggered Dirac operator we have studied random matrices modeling the transition from non-degenerate eigenvalues at non-zero lattice spacing to degenerate ones in the continuum limit.