Source author record

Karl Jansen

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

128works
21topics
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

128 published item(s)

preprint2026arXiv

Computing Statistical Properties of Velocity Fields on Current Quantum Hardware

Quantum algorithms are gaining attention in Computational Fluid Dynamics (CFD) for their favorable scaling, as encoding physical fields into quantum probability amplitudes enables representation of two to the power of n spatial points with only n qubits. A key challenge in Quantum CFD is the efficient readout of simulation results, a topic that has received limited attention in literature. This work presents methods to extract statistical properties of spatial velocity fields, such as central moments and structure functions, directly from parameterized ansatz circuits, avoiding full quantum state tomography. As a proof of concept, we implement our approach for 1D velocity fields, encoding 16 spatial points with 4 qubits, and analyze both a sine wave signal and four snapshots from Burgers' equation evolution. Using Qedma's error mitigation software QESEM, we demonstrate that such computations achieve high accuracy on current quantum devices, specifically IBMQ's Heron2 system ibm_fez.

preprint2025arXiv

Processing through encoding: Quantum circuit approaches for point-wise multiplication and convolution

This paper introduces quantum circuit methodologies for pointwise multiplication and convolution of complex functions, conceptualized as "processing through encoding". Leveraging known techniques, we describe an approach where multiple complex functions are encoded onto auxiliary qubits. Applying the proposed scheme for two functions $f$ and $g$, their pointwise product $f(x)g(x)$ is shown to naturally form as the coefficients of part of the resulting quantum state. Adhering to the convolution theorem, we then demonstrate how the convolution $f*g$ can be constructed. Similarly to related work, this involves the encoding of the Fourier coefficients $\mathcal{F}[f]$ and $\mathcal{F}[g]$, which facilitates their pointwise multiplication, followed by the inverse Quantum Fourier Transform. We discuss the simulation of these techniques, their integration into an extended \verb|quantumaudio| package for audio signal processing, and present initial experimental validations. This work offers a promising avenue for quantum signal processing, with potential applications in areas such as quantum-enhanced audio manipulation and synthesis.

preprint2022arXiv

Digitising SU(2) Gauge Fields and the Freezing Transition

Efficient discretisations of gauge groups are crucial with the long term perspective of using tensor networks or quantum computers for lattice gauge theory simulations. For any Lie group other than U$(1)$, however, there is no class of asymptotically dense discrete subgroups. Therefore, discretisations limited to subgroups are bound to lead to a freezing of Monte Carlo simulations at weak couplings, necessitating alternative partitionings without a group structure. In this work we provide a comprehensive analysis of this freezing for all discrete subgroups of SU$(2)$ and different classes of asymptotically dense subsets. We find that an appropriate choice of the subset allows unfrozen simulations for arbitrary couplings, though one has to be careful with varying weights of unevenly distributed points. A generalised version of the Fibonacci spiral appears to be particularly efficient and close to optimal.

preprint2022arXiv

Lattice QCD Study of Transverse-Momentum Dependent Soft Function

In this work, we perform a lattice QCD study of the intrinsic, rapidity-independent soft function within the framework of large momentum effective theory. The computation is carried out using a gauge ensemble of $N_f=2+1+1$ clover-improved twisted mass fermion. After applying an appropriate renormalization procedure and the removal of significant higher-twist contamination, we obtain the intrinsic soft function that is comparable to the one-loop perturbative result at large external momentum. The determination of the nonperturbative soft function from first principles is crucial to sharpen our understanding of the processes with small transverse momentum such as the Drell-Yan production and the semi-inclusive deep inelastic scattering. Additionally, we calculate the Collins-Soper evolution kernel using the quasi-transverse-momentum-dependent wave function as input.

preprint2022arXiv

Measurement Error Mitigation in Quantum Computers Through Classical Bit-Flip Correction

We develop a classical bit-flip correction method to mitigate measurement errors on quantum computers. This method can be applied to any operator, any number of qubits, and any realistic bit-flip probability. We first demonstrate the successful performance of this method by correcting the noisy measurements of the ground-state energy of the longitudinal Ising model. We then generalize our results to arbitrary operators and test our method both numerically and experimentally on IBM quantum hardware. As a result, our correction method reduces the measurement error on the quantum hardware by up to one order of magnitude. We finally discuss how to pre-process the method and extend it to other errors sources beyond measurement errors. For local Hamiltonians, the overhead costs are polynomial in the number of qubits, even if multi-qubit correlations are included.

preprint2022arXiv

New Directions in Quantum Music: concepts for a quantum keyboard and the sound of the Ising model

We explore ideas for generating sounds and eventually music by using quantum devices in the NISQ era using quantum circuits. In particular, we first consider a concept for a "qeyboard", i.e. a quantum keyboard, where the real-time behaviour of expectation values using a time evolving quantum circuit can be associated to sound features like intensity, frequency and tone. Then, we examine how these properties can be extracted from physical quantum systems, taking the Ising model as an example. This can be realized by measuring physical quantities of the quantum states of the system, e.g. the energies and the magnetization obtained via variational quantum simulation techniques.

preprint2022arXiv

Phase structure of the CP(1) model in the presence of a topological $θ$-term

We numerically study the phase structure of the CP(1) model in the presence of a topological $θ$-term, a regime afflicted by the sign problem for conventional lattice Monte Carlo simulations. Using a bond-weighted tensor renormalization group method, we compute the free energy for inverse couplings ranging from $0\leq β\leq 1.1$ and find a CP-violating, first-order phase transition at $θ=π$. In contrast to previous findings, our numerical results provide no evidence for a critical coupling $β_c<1.1$ above which a second-order phase transition emerges at $θ=π$ and/or the first-order transition line bifurcates at $θ\neqπ$. If such a critical coupling exists, as suggested by Haldane's conjecture, our study indicates that is larger than $β_c>1.1$.

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

Proton generalized parton distributions from lattice QCD

Momentum and spatial distributions of quarks and gluons inside hadrons are typically encoded in the so-called generalized parton distributions (GPDs). GPDs are multi-dimensional quantities that are very challenging to extract, both experimentally and within lattice QCD. We present the first lattice results on the $x$-dependence of isovector unpolarized, helicity and transversity GPDs of the proton, obtained from lattice QCD using an ensemble of $N_f=2+1+1$ maximally twisted mass fermions, with pion mass $M_π=260$ MeV and lattice spacing $a\simeq 0.093$ fm. Our calculations use the quasi-distribution formalism and the final distributions are presented in the MS-bar scheme at a renormalization scale of 2 GeV.

preprint2022arXiv

Snowmass White Paper: Quantum Computing Systems and Software for High-energy Physics Research

Quantum computing offers a new paradigm for advancing high-energy physics research by enabling novel methods for representing and reasoning about fundamental quantum mechanical phenomena. Realizing these ideals will require the development of novel computational tools for modeling and simulation, detection and classification, data analysis, and forecasting of high-energy physics (HEP) experiments. While the emerging hardware, software, and applications of quantum computing are exciting opportunities, significant gaps remain in integrating such techniques into the HEP community research programs. Here we identify both the challenges and opportunities for developing quantum computing systems and software to advance HEP discovery science. We describe opportunities for the focused development of algorithms, applications, software, hardware, and infrastructure to support both practical and theoretical applications of quantum computing to HEP problems within the next 10 years.

preprint2022arXiv

Strategies for the Determination of the Running Coupling of $(2+1)$-dimensional QED with Quantum Computing

We propose to utilize NISQ-era quantum devices to compute short distance quantities in $(2+1)$-dimensional QED and to combine them with large volume Monte Carlo simulations and perturbation theory. On the quantum computing side, we perform a calculation of the mass gap in the small and intermediate regime, demonstrating, in the latter case, that it can be resolved reliably. The so obtained mass gap can be used to match corresponding results from Monte Carlo simulations, which can be used eventually to set the physical scale. In this paper we provide the setup for the quantum computation and show results for the mass gap and the plaquette expectation value. In addition, we discuss some ideas that can be applied to the computation of the running coupling. Since the theory is asymptotically free, it would serve as a training ground for future studies of QCD in $(3+1)$-dimensions on quantum computers.

preprint2022arXiv

Transversity GPDs of the proton from lattice QCD

We present the first calculation of the $x$-dependence of the isovector transversity generalized parton distributions (GPDs) for the proton within lattice QCD. We compute the matrix elements with non-local operators containing a Wilson line. The calculation implements the Breit symmetric frame. The proton momenta are chosen as $0.83,\,1.25,\,1.67$ GeV, and the values of the momentum transfer squared are $0.69,\,1.02$ GeV$^2$. These combinations include cases with zero and nonzero skewness. The calculation is performed using one ensemble of two degenerate-mass light, a strange and a charm quark of maximally twisted mass fermions with a clover term. The lattice results are renormalized non-perturbatively and finally matched to the light-cone GPDs using one-loop perturbation theory within the framework of large momentum effective theory. The final GPDs are given in the $\overline{\rm MS}$ scheme at a scale of 2 GeV. In addition to the individual GPDs, we form the combination of the transversity GPDs that is related to the transverse spin structure of the proton. Finally, we extract the lowest two moments of GPDs and draw a number of important qualitative conclusions.

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.

preprint2022arXiv

Using classical bit-flip correction for error mitigation including 2-qubit correlations

We present an error mitigation scheme which corrects readout errors on Noisy Intermediate-Scale Quantum (NISQ) computers [1,2]. After a short review of applying the method to one qubit, we proceed to discuss the case when correlations between different qubits occur. We demonstrate how the readout error can be mitigated in this case. By performing experiments on IBMQ hardware, we show that such correlations do not have a strong effect on the results, justifying to neglect them.

preprint2022arXiv

Zeta-regularized Lattice Field Theory with Lorentzian background metrics

Lattice field theory is a very powerful tool to study Feynman's path integral non-perturbatively. However, it usually requires Euclidean background metrics to be well-defined. On the other hand, a recently developed regularization scheme based on Fourier integral operator $ζ$-functions can treat Feynman's path integral non-pertubatively in Lorentzian background metrics. In this article, we formally $ζ$-regularize lattice theories with Lorentzian backgrounds and identify conditions for the Fourier integral operator $ζ$-function regularization to be applicable. Furthermore, we show that the classical limit of the $ζ$-regularized theory is independent of the regularization. Finally, we consider the harmonic oscillator as an explicit example. We discuss multiple options for the regularization and analytically show that they all reproduce the correct ground state energy on the lattice and in the continuum limit. Additionally, we solve the harmonic oscillator on the lattice in Minkowski background numerically.

preprint2021arXiv

A resource efficient approach for quantum and classical simulations of gauge theories in particle physics

Gauge theories establish the standard model of particle physics, and lattice gauge theory (LGT) calculations employing Markov Chain Monte Carlo (MCMC) methods have been pivotal in our understanding of fundamental interactions. The present limitations of MCMC techniques may be overcome by Hamiltonian-based simulations on classical or quantum devices, which further provide the potential to address questions that lay beyond the capabilities of the current approaches. However, for continuous gauge groups, Hamiltonian-based formulations involve infinite-dimensional gauge degrees of freedom that can solely be handled by truncation. Current truncation schemes require dramatically increasing computational resources at small values of the bare couplings, where magnetic field effects become important. Such limitation precludes one from `taking the continuous limit' while working with finite resources. To overcome this limitation, we provide a resource-efficient protocol to simulate LGTs with continuous gauge groups in the Hamiltonian formulation. Our new method allows for calculations at arbitrary values of the bare coupling and lattice spacing. The approach consists of the combination of a Hilbert space truncation with a regularization of the gauge group, which permits an efficient description of the magnetically-dominated regime. We focus here on Abelian gauge theories and use $2+1$ dimensional quantum electrodynamics as a benchmark example to demonstrate this efficient framework to achieve the continuum limit in LGTs. This possibility is a key requirement to make quantitative predictions at the field theory level and offers the long-term perspective to utilise quantum simulations to compute physically meaningful quantities in regimes that are precluded to quantum Monte Carlo.

preprint2021arXiv

Estimation of Thermodynamic Observables in Lattice Field Theories with Deep Generative Models

In this work, we demonstrate that applying deep generative machine learning models for lattice field theory is a promising route for solving problems where Markov Chain Monte Carlo (MCMC) methods are problematic. More specifically, we show that generative models can be used to estimate the absolute value of the free energy, which is in contrast to existing MCMC-based methods which are limited to only estimate free energy differences. We demonstrate the effectiveness of the proposed method for two-dimensional $ϕ^4$ theory and compare it to MCMC-based methods in detailed numerical experiments.

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.

preprint2021arXiv

Towards Quantum Simulations in Particle Physics and Beyond on Noisy Intermediate-Scale Quantum Devices

We review two algorithmic advances that bring us closer to reliable quantum simulations of model systems in high energy physics and beyond on noisy intermediate-scale quantum (NISQ) devices. The first method is the dimensional expressivity analysis of quantum circuits, which allows for constructing minimal but maximally expressive quantum circuits. The second method is an efficient mitigation of readout errors on quantum devices. Both methods can lead to significant improvements in quantum simulations, e.g., when variational quantum eigensolvers are used.

preprint2020arXiv

Avoiding the sign-problem in lattice field theory

In lattice field theory, the interactions of elementary particles can be computed via high-dimensional integrals. Markov-chain Monte Carlo (MCMC) methods based on importance sampling are normally efficient to solve most of these integrals. But these methods give large errors for oscillatory integrands, exhibiting the so-called sign-problem. We developed new quadrature rules using the symmetry of the considered systems to avoid the sign-problem in physical one-dimensional models for the resulting high-dimensional integrals. This article gives a short introduction to integrals used in lattice QCD where the interactions of gluon and quark elementary particles are investigated, explains the alternative integration methods we developed and shows results of applying them to models with one physical dimension. The new quadrature rules avoid the sign-problem and can therefore be used to perform simulations at until now not reachable regions in parameter space, where the MCMC errors are too big for affordable sample sizes. However, it is still a challenge to develop these techniques further for applications with physical higher-dimensional systems.

preprint2020arXiv

Comparison of topological charge definitions in Lattice QCD

In this paper, we show a comparison of different definitions of the topological charge on the lattice. We concentrate on one small-volume ensemble with 2 flavours of dynamical, maximally twisted mass fermions and use three more ensembles to analyze the approach to the continuum limit. We investigate several fermionic and gluonic definitions. The former include the index of the overlap Dirac operator, the spectral flow of the Wilson--Dirac operator and the spectral projectors. For the latter, we take into account different discretizations of the topological charge operator and various smoothing schemes to filter out ultraviolet fluctuations: the gradient flow, stout smearing, APE smearing, HYP smearing and cooling. We show that it is possible to perturbatively match different smoothing schemes and provide a well-defined smoothing scale. We relate the smoothing parameters for cooling, stout and APE smearing to the gradient flow time $τ$. In the case of hypercubic smearing the matching is performed numerically. We investigate which conditions have to be met to obtain a valid definition of the topological charge and susceptibility and we argue that all valid definitions are highly correlated and allow good control over topology on the lattice.

preprint2020arXiv

From spin chains to real-time thermal field theory using tensor networks

One of the most interesting directions in theoretical high-energy and condensed-matter physics is understanding dynamical properties of collective states of quantum field theories. The most elementary tool in this quest is retarded equilibrium correlators governing the linear response theory. In this article we examine tensor networks as a way of determining them in a fully ab initio way in a class of (1+1)-dimensional quantum field theories arising as infrared descriptions of quantum Ising chains.We show that, complemented with signal analysis using the Prony method, tensor network calculations for intermediate times provide a powerful way to explore the structure of singularities of the correlator in the complex frequency plane and to make predictions about the thermal response to perturbations in a class of nonintegrable interacting quantum field theories.

preprint2020arXiv

Improvement, generalization, and scheme conversion of Wilson-line operators on the lattice in the auxiliary field approach

Nonlocal quark bilinear operators connected by link paths are used for studying parton distribution functions (PDFs) and transverse momentum-dependent PDFs of hadrons using lattice QCD. The nonlocality makes it difficult to understand the renormalization and improvement of these operators using standard methods. In previous work, we showed that by introducing an auxiliary field on the lattice, one can understand an on-axis Wilson-line operator as the product of two local operators in an extended theory. In this paper, we provide details about the calculation in perturbation theory of the factor for conversion from our lattice-suitable renormalization scheme to the MS-bar scheme. Extending our work, we study Symanzik improvement of the extended theory to understand the pattern of discretization effects linear in the lattice spacing, $a$, which are present even if the lattice fermion action exactly preserves chiral symmetry. This provides a prospect for an eventual $O(a)$ improvement of lattice calculations of PDFs. We also generalize our approach to apply to Wilson lines along lattice diagonals and to piecewise-straight link paths.

preprint2020arXiv

Multilevel Monte Carlo for quantum mechanics on a lattice

Monte Carlo simulations of quantum field theories on a lattice become increasingly expensive as the continuum limit is approached since the cost per independent sample grows with a high power of the inverse lattice spacing. Simulations on fine lattices suffer from critical slowdown, the rapid growth of autocorrelations in the Markov chain. This causes a strong increase in the number of lattice configurations that have to be generated to obtain statistically significant results. This paper discusses hierarchical sampling methods to tame the growth in autocorrelations. Combined with multilevel variance reduction, this significantly reduces the computational cost of simulations for given tolerances $ε_{\text{disc}}$ on the discretisation error and $ε_{\text{stat}}$ on the statistical error. For observables with lattice errors of order $α$ and integrated autocorrelation times that grow like $τ_{\mathrm{int}}\propto a^{-z}$, multilevel Monte Carlo (MLMC) reduces the cost from $\mathcal{O}(ε_{\text{stat}}^{-2}ε_{\text{disc}}^{-(1+z)/α})$ to $\mathcal{O}(ε_{\text{stat}}^{-2}\vert\log ε_{\text{disc}} \vert^2+ε_{\text{disc}}^{-1/α})$ or $\mathcal{O}(ε_{\text{stat}}^{-2}+ε_{\text{disc}}^{-1/α})$. Higher gains are expected for simulations of quantum field theories in $D$ dimensions. The efficiency of the approach is demonstrated on two model systems, including a topological oscillator that is badly affected by critical slowdown from topological charge freezing. On fine lattices, the new methods are orders of magnitude faster than standard Hybrid Monte Carlo sampling. For high resolutions, MLMC can be used to accelerate even the cluster algorithm for the topological oscillator. Performance is further improved through perturbative matching which guarantees efficient coupling of theories on the multilevel hierarchy.

preprint2020arXiv

Parton distribution functions of $Δ^+$ on the lattice

We perform a first calculation for the unpolarized parton distribution function of the $Δ^+$ baryon using lattice QCD simulations within the framework of Large Momentum Effective Theory. Two ensembles of $N_f=2+1+1$ twisted mass fermions are utilized with a pion mass of 270 MeV and 360 MeV, respectively. The baryon, which is treated as a stable single-particle state, is boosted with momentum $P_3$ with values $\{0.42,0.83,1.25\}$ GeV, and we utilize momentum smearing to improve the signal. The unpolarized parton distribution function of $Δ^+$ is obtained using a non-perturbative renormalization and a one-loop formula for the matching, with encouraging precision. In particular, we compute the $\overline{d}(x)-\overline{u}(x)$ asymmetry and compare it with the same quantity in the nucleon, in a first attempt towards resolving the physical mechanism responsible for generating such asymmetry.

preprint2020arXiv

Topological vacuum structure of the Schwinger model with matrix product states

We numerically study the single-flavor Schwinger model with a topological $θ$-term, which is practically inaccessible by standard lattice Monte Carlo simulations due to the sign problem. By using numerical methods based on tensor networks, especially the one-dimensional matrix product states, we explore the non-trivial $θ$-dependence of several lattice and continuum quantities in the Hamiltonian formulation. In particular, we compute the ground-state energy, the electric field, the chiral fermion condensate, and the topological vacuum susceptibility for positive, zero, and even negative fermion mass. In the chiral limit, we demonstrate that the continuum model becomes independent of the vacuum angle $θ$, thus respecting CP invariance, while lattice artifacts still depend on $θ$. We also confirm that negative masses can be mapped to positive masses by shifting $θ\rightarrow θ+π$ due to the axial anomaly in the continuum, while lattice artifacts non-trivially distort this mapping. This mass regime is particularly interesting for the (3+1)-dimensional QCD analog of the Schwinger model, the sign problem of which requires the development and testing of new numerical techniques beyond the conventional Monte Carlo approach.

preprint2020arXiv

Unpolarized and helicity generalized parton distributions of the proton within lattice QCD

We present the first calculation of the $x$-dependence of the proton generalized parton distributions (GPDs) within lattice QCD. Results are obtained for the isovector unpolarized and helicity GPDs. We compute the appropriate matrix elements of fast-moving protons coupled to non-local operators containing a Wilson line. We present results for proton momenta $0.83,\,1.25,\,1.67$ GeV, and momentum transfer squared $0.69,\,1.38$ GeV$^2$. These combinations include cases with zero and nonzero skewness. The calculation is performed using one ensemble of two degenerate mass light, a strange and a charm quark of maximally twisted mass fermions with a clover term. The lattice results are matched to the light-cone GPDs using one-loop perturbation theory within the framework of large momentum effective theory. The final GPDs are given in the $\overline{\rm MS}$ scheme at a scale of 2 GeV.

preprint2016arXiv

A Lattice Calculation of Parton Distributions

We present results for the $x$ dependence of the unpolarized, helicity, and transversity isovector quark distributions in the proton using lattice QCD, employing the method of quasi-distributions proposed by Ji in 2013. Compared to a previous calculation by us, the errors are reduced by a factor of about 2.5. Moreover, we present our first results for the polarized sector of the proton, which indicate an asymmetry in the proton sea in favor of the $u$ antiquarks for the case of helicity distributions, and an asymmetry in favor of the $d$ antiquarks for the case of transversity distributions.

preprint2016arXiv

Applying recursive numerical integration techniques for solving high dimensional integrals

The error scaling for Markov-Chain Monte Carlo techniques (MCMC) with $N$ samples behaves like $1/\sqrt{N}$. This scaling makes it often very time intensive to reduce the error of computed observables, in particular for applications in lattice QCD. It is therefore highly desirable to have alternative methods at hand which show an improved error scaling. One candidate for such an alternative integration technique is the method of recursive numerical integration (RNI). The basic idea of this method is to use an efficient low-dimensional quadrature rule (usually of Gaussian type) and apply it iteratively to integrate over high-dimensional observables and Boltzmann weights. We present the application of such an algorithm to the topological rotor and the anharmonic oscillator and compare the error scaling to MCMC results. In particular, we demonstrate that the RNI technique shows an error scaling in the number of integration points $m$ that is at least exponential.

preprint2016arXiv

Chiral condensate in the Schwinger model with Matrix Product Operators

Tensor network (TN) methods, in particular the Matrix Product States (MPS) ansatz, have proven to be a useful tool in analyzing the properties of lattice gauge theories. They allow for a very good precision, much better than standard Monte Carlo (MC) techniques for the models that have been studied so far, due to the possibility of reaching much smaller lattice spacings. The real reason for the interest in the TN approach, however, is its ability, shown so far in several condensed matter models, to deal with theories which exhibit the notorious sign problem in MC simulations. This makes it prospective for dealing with the non-zero chemical potential in QCD and other lattice gauge theories, as well as with real-time simulations. In this paper, using matrix product operators, we extend our analysis of the Schwinger model at zero temperature to show the feasibility of this approach also at finite temperature. This is an important step on the way to deal with the sign problem of QCD. We analyze in detail the chiral symmetry breaking in the massless and massive cases and show that the method works very well and gives good control over a broad range of temperatures, essentially from zero to infinite temperature.

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.

preprint2016arXiv

Finite Size Scaling of the Higgs-Yukawa Model near the Gaussian Fixed Point

We study the scaling properties of Higgs-Yukawa models. Using the technique of Finite-Size Scaling, we are able to derive scaling functions that describe the observables of the model in the vicinity of a Gaussian fixed point. A feasibility study of our strategy is performed for the pure scalar theory in the weak-coupling regime. Choosing the on-shell renormalisation scheme gives us an advantage to fit the scaling functions against lattice data with only a small number of fit parameters. These formulae can be used to determine the universality of the observed phase transitions, and thus play an essential role in future investigations of Higgs-Yukawa models, in particular in the strong Yukawa coupling region.

preprint2016arXiv

New polynomially exact integration rules on U(N) and SU(N)

In lattice Quantum Field Theory, we are often presented with integrals over polynomials of coefficients of matrices in U(N) or SU(N) with respect to the Haar measure. In some physical situations, e.g., in presence of a chemical potential, these integrals are numerically very difficult since their integrands are highly oscillatory which manifests itself in form of the sign problem. In these cases, Monte Carlo methods often fail to be adequate, rendering such computations practically impossible. We propose a new class of integration rules on U(N) and SU(N) which are derived from polynomially exact rules on spheres. We will examine these quadrature rules and their efficiency at the example of a 0+1 dimensional QCD for a non-zero quark mass and chemical potential. In particular, we will demonstrate the failure of Monte Carlo methods in such applications and that we can obtain polynomially exact, arbitrary precision results using the new integration rules.

preprint2016arXiv

Non-perturbative running of renormalization constants from correlators in coordinate space using step scaling

Working in a quenched setup with Wilson twisted mass valence fermions, we explore the possibility to compute non-perturbatively the step scaling function using the coordinate (X-space) renormalization scheme. This scheme has the advantage of being on-shell and gauge invariant. The step scaling method allows us to calculate the running of the renormalization constants of quark bilinear operators. We describe here the details of this calculation. The aim of this exploratory study is to identify the feasibility of the X-space scheme when used in small volume simulations required by the step scaling technique. Eventually, we translate our final results to the continuum MSbar scheme and compare against four-loop analytic formulae finding satisfactory agreement.

preprint2016arXiv

Recent results for the proton spin decomposition from lattice QCD

The exact decomposition of the proton spin has been a much debated topic, on the experimental as well as the theoretical side. In this talk we would like to report on recent non-perturbative results and ongoing efforts to explore the proton spin from lattice QCD. We present results for the relevant generalized form factors from gauge field ensembles that feature a physical value of the pion mass. These generalized form factors can be used to determine the total spin and angular momentum carried by the quarks. In addition we present first results for our ongoing effort to compute the angular momentum of the gluons in the proton.

preprint2016arXiv

Step scaling in coordinate space: running of the quark mass

We perform a benchmark study of the step scaling procedure for the ratios of renormalization constants extracted from position space correlation functions. We work in the quenched approximation and consider the pseudoscalar, scalar, vector and axial vector bilinears. The pseudoscalar/scalar cases allow us to obtain the non-perturbative running of the quark mass over a wide range of energy scales - from around 17 GeV to below 1.5 GeV - which agrees well with the 4-loop prediction of continuum perturbation theory. We find that step scaling is feasible in X-space and we discuss its advantages and potential problems.

preprint2016arXiv

The multi-flavor Schwinger model with chemical potential - Overcoming the sign problem with Matrix Product States

During recent years there has been an increasing interest in the application of matrix product states, and more generally tensor networks, to lattice gauge theories. This non-perturbative method is sign problem free and has already been successfully used to compute mass spectra, thermal states and phase diagrams, as well as real-time dynamics for Abelian and non-Abelian gauge models. In previous work we showed the suitability of the method to explore the zero-temperature phase structure of the multi-flavor Schwinger model at non-zero chemical potential, a regime where the conventional Monte Carlo approach suffers from the sign problem. Here we extend our numerical study by looking at the spatially resolved chiral condensate in the massless case. We recover spatial oscillations, similar to the theoretical predictions for the single-flavor case, with a chemical potential dependent frequency and an amplitude approximately given by the homogeneous zero density condensate value.

preprint2015arXiv

A Lattice Calculation of Parton Distributions

We report on our exploratory study for the direct evaluation of the parton distribution functions from lattice QCD, based on a recently proposed new approach. We present encouraging results using Nf = 2 + 1 + 1 twisted mass fermions with a pion mass of about 370 MeV. The focus of this work is a detailed description of the computation, including the lattice calculation, the matching to an infinite momentum and the nucleon mass correction. In addition, we test the effect of gauge link smearing in the operator to estimate the influence of the Wilson line renormalization, which is yet to be done.

preprint2015arXiv

A lattice study of a chirally invariant Higgs-Yukawa model including a higher dimensional $Φ^6$-term

We discuss the non-thermal phase structure of a chirally invariant Higgs-Yukawa model on the lattice in the presence of a higher dimensional $Φ^6$-term. For the exploration of the phase diagram we use analytical, lattice perturbative calculations of the constraint effectice potential as well as numerical simulations. We also present first results of the effects of the $Φ^6$-term on the lower Higgs boson mass bounds.

preprint2015arXiv

Disconnected quark loop contributions to nucleon observables using $N_f=2$ twisted clover fermions at the physical value of the light quark mass

We compute the disconnected quark loops contributions entering the determination of nucleon observables, by using a $N_f = 2$ ensemble of twisted mass fermions with a clover term at a pion mass $m_π= 133$ MeV. We employ exact deflation and implement all calculations in GPUs, enabling us to achieve large statistics and a good signal.

preprint2015arXiv

Lattice study of the Higgs-Yukawa model in and beyond the Standard Model

We derive finite-size scaling formulae for four-dimensional Higgs-Yukawa models near the Gaussian fixed point. These formulae will play an essential role in future, detailed investigation of such models. In particular, they can be used to determine the nature of the observed phase transitions, and confirm or rule out the possibility of having non-trivial fixed points in the Higgs-Yukawa models. Our scaling formula for Binder's cumulant is tested against lattice simulations carried out at weak couplings, and good agreement is found. As a separate project, we also present preliminary results from our study of a chirally-invariant Higgs-Yukawa model including a dimension-six operator at finite temperature. Our work provides first indications of first-order temperature-induced phase transitions near the infinite cutoff limit in this model.

preprint2015arXiv

Leading hadronic contributions to the running of the electroweak coupling constants from lattice QCD

The quark-connected leading-order hadronic contributions to the running of the electromagnetic fine structure constant, $α_{\rm QED}$, and the weak mixing angle, $θ_W$, are determined by a four-flavour lattice QCD computation with twisted mass fermions. Full agreement of the results with a phenomenological analysis is observed with an even comparable statistical uncertainty. We show that the uncertainty of the lattice calculation is dominated by systematic effects which then leads to significantly larger errors than obtained by the phenomenological analysis.

preprint2015arXiv

Leading-order hadronic contributions to the lepton anomalous magnetic moments from the lattice

The hadronic leading-order (hlo) contribution to the lepton anomalous magnetic moments $a_l^\mathrm{hlo}$ of the Standard Model leptons still accounts for the dominant source of the uncertainty of the Standard Model estimates. We present the results of an investigation of the hadronic leading order anomalous magnetic moments of the electron, muon and tau lepton from first principles in twisted mass lattice QCD. With lattice data for multiple pion masses in the range $230 \mathrm{~MeV} \lesssim m_{PS} \lesssim 490 \mathrm{~MeV}$, multiple lattice volumes and three lattice spacings we perform the extrapolation to the continuum and to the physical pion mass and check for all systematic uncertainties in the lattice calculation. As a result we calculate $a_{l}^\mathrm{hlo}$ for the three Standard Model leptons with controlled statistical and systematic error in agreement with phenomenological determinations using dispersion relations and experimental data. In addition, we also give a first estimate of the hadronic leading order anomalous magnetic moments from simulations directly at the physical value of the pion mass.

preprint2015arXiv

Non-perturbative Test of the Witten-Veneziano Formula from Lattice QCD

We compute both sides of the Witten-Veneziano formula using lattice techniques. For the one side we perform dedicated quenched simulations and use the spectral projector method to determine the topological susceptibility in the pure Yang-Mills theory. The other side we determine in lattice QCD with $N_f=2+1+1$ dynamical Wilson twisted mass fermions including for the first time also the flavour singlet decay constant. The Witten-Veneziano formula represents a leading order expression in the framework of chiral perturbation theory and we also employ leading order chiral perturbation theory to relate the flavor singlet decay constant to the relevant decay constant parameters in the quark flavor basis and flavor non-singlet decay constants. After taking the continuum and the SU$(2)$ chiral limits we compare both sides and find good agreement within uncertainties.

preprint2015arXiv

Nucleon electromagnetic form factors from twisted mass lattice QCD

Results on the electromagnetic form factors of the nucleon using twisted mass fermion configurations are presented. These include a gauge field ensemble simulated with two degenerate light quarks yielding a pion mass of around 130 MeV, as well as two ensembles that include strange and charm quarks in the sea yielding pion masses of 210 MeV and 373 MeV. Details of the methods used and systematic errors are discussed, such as noise reduction techniques and the effect of excited states contamination.

preprint2015arXiv

Phase structure and Higgs boson mass in a Higgs-Yukawa model with a dimension-6 operator

We investigate the impact of a $λ_6 φ^6$ term included in a chirally invariant lattice Higgs-Yukawa model. Such a term could emerge from BSM physics at some larger energy scale. We map out the phase structure of the Higgs-Yukawa model with positive $λ_6$ and negative quartic self coupling of the scalar fields. To this end, we evaluate the constraint effective potential in lattice perturbation theory and also determine the magnetization of the model via numerical simulations which allow us to reach also non-perturbative values of the couplings. As a result, we find a complex phase structure with first and second order phase transitions identified through the magnetization. Further we analyze the effect of such a $φ^6$ term on the lower Higgs boson mass bound to see, whether the standard model lower mass bound can be altered.

preprint2015arXiv

Short distance singularities and automatic O($a$) improvement: the cases of the chiral condensate and the topological susceptibility

Short-distance singularities in lattice correlators can modify their Symanzik expansion by leading to additional O($a$) lattice artifacts. At the example of the chiral condensate and the topological susceptibility, we show how to account for these lattice artifacts for Wilson twisted mass fermions and show that the property of automatic O($a$) improvement is preserved at maximal twist.

preprint2015arXiv

Strangeness of the nucleon from Lattice Quantum Chromodynamics

We present a non-perturbative calculation of the strangeness of the nucleon $y_N$ within the framework of lattice QCD. This observable is known to be an important cornerstone to interpret results from direct dark matter detection experiments. We perform a lattice computation for $y_N$ with an analysis of systematic effects originating from discretization, finite size, chiral extrapolation and excited state effects leading to the value of $y_N= 0.173(50)$. The rather large uncertainty of this value of $y_N$ is dominated by systematic uncertainties which we were able to quantify in this work.

preprint2015arXiv

Testing the Witten-Veneziano Formula on the Lattice

In this proceeding contribution we report on a test of the famous Witten-Veneziano formula using lattice techniques. We perform dedicated quenched simulations and apply the spectral projector method to determine the topological susceptibility in pure Yang-Mills theory. In order to compute the relevant meson masses and the flavor singlet decay constant we employ lattice QCD with $N_f=2+1+1$ dynamical Wilson twisted mass fermions. Taking the continuum and the $SU(2)$ chiral limits we find good agreement within uncertainties for both sides of the formula.

preprint2015arXiv

The electric dipole moment of the neutron from $N_f=2+1+1$ twisted mass fermions

We extract the neutron electric dipole moment (nEDM) $ \vert \vec{d}_n \vert $ on configurations produced with $N_f=2+1+1$ twisted mass fermions with lattice spacing of $a \simeq 0.082$ fm and a light quark mass that corresponds to $M_π \simeq 373$ MeV. We do so by evaluating the $CP$-odd form factor $F_3$ for small values of the $CP$-violation parameter $θ$ in the limit of zero momentum transfer. This limit is extracted using the usual parametrization but in addition position space methods. The topological charge is computed via cooling and gradient flow using the Wilson, Symanzik tree-level improved and Iwasaki actions for smoothing. We obtain consistent results for all choices of smoothing procedures and methods to extract $F_3$ at zero momentum transfer. For the ensemble analyzed we find a value of nEDM of $\vert \vec{d}_n \vert / θ= -0.045(6)(1) {\rm e} \cdot {\rm fm}$.

preprint2015arXiv

The hadronic vacuum polarization and automatic O(a) improvement for twisted mass fermions

The vacuum polarization tensor and the corresponding vacuum polarization function are the basis for calculations of numerous observables in lattice QCD. Examples are the hadronic contributions to lepton anomalous magnetic moments, the running of the electroweak and strong couplings and quark masses. Quantities which are derived from the vacuum polarization tensor often involve a summation of current correlators over all distances in position space leading thus to the appearance of short-distance terms. The mechanism of O(a) improvement in the presence of such short-distance terms is not directly covered by the usual arguments of on-shell improvement of the action and the operators for a given quantity. If such short-distance contributions appear, the property of O(a) improvement needs to be reconsidered. We discuss the effects of these short-distance terms on the vacuum polarization function for twisted mass lattice QCD and find that even in the presence of such terms automatic O(a) improvement is retained if the theory is tuned to maximal twist.

preprint2014arXiv

A Stochastic Method for Computing Hadronic Matrix Elements

We present a stochastic method for the calculation of baryon three-point functions that is more versatile compared to the typically used sequential method. We analyze the scaling of the error of the stochastically evaluated three-point function with the lattice volume and find a favorable signal-to-noise ratio suggesting that our stochastic method can be used efficiently at large volumes to compute hadronic matrix elements.

preprint2014arXiv

Computation of disconnected contributions to nucleon observables

We compare several methods for computing disconnected fermion loops contributing to nucleon three-point functions. The comparison is carried out using one ensemble of $N_f=2+1+1$ twisted mass fermions with pion mass of 373 MeV. The complete set of operators up to one-derivative are examined by developing optimized code for mutli-GPUs. Simple guidelines are given as to the preferable method for each class of operators.

preprint2014arXiv

First results with twisted mass fermions towards the computation of parton distribution functions on the lattice

We report on our exploratory study for the evaluation of the parton distribution functions from lattice QCD, based on a new method proposed in Ref.~arXiv:1305.1539. Using the example of the nucleon, we compare two different methods to compute the matrix elements needed, and investigate the application of gauge link smearing. We also present first results from a large production ensemble and discuss the future challenges related to this method.

preprint2014arXiv

Four-Flavour Leading-Order Hadronic Contribution To The Muon Anomalous Magnetic Moment

We present a four-flavour lattice calculation of the leading-order hadronic vacuum polarisation contribution to the anomalous magnetic moment of the muon, $a_\mathrmμ^{\rm hvp}$, arising from quark-connected Feynman graphs. It is based on ensembles featuring $N_f=2+1+1$ dynamical twisted mass fermions generated by the European Twisted Mass Collaboration (ETMC). Several light quark masses are used in order to yield a controlled extrapolation to the physical pion mass. We employ three lattice spacings to examine lattice artefacts and several different volumes to check for finite-size effects. Incorporating the complete first two generations of quarks allows for a direct comparison with phenomenological determinations of $a_\mathrmμ^{\rm hvp}$. Our final result including an estimate of the systematic uncertainty $$a_{\mathrmμ}^{\rm hvp} = 6.74(21)(18) \cdot 10^{-8}$$ shows a good overall agreement with these computations.

preprint2014arXiv

Hadronic contributions to the muon anomalous magnetic moment Workshop. $(g-2)_μ$: Quo vadis? Workshop. Mini proceedings

We present the mini-proceedings of the workshops Hadronic contributions to the muon anomalous magnetic moment: strategies for improvements of the accuracy of the theoretical prediction and $(g-2)_μ$: Quo vadis?, both held in Mainz from April 1$^{\rm rst}$ to 5$^{\rm th}$ and from April 7$^{\rm th}$ to 10$^{\rm th}$, 2014, respectively.

preprint2014arXiv

Lepton anomalous magnetic moments from twisted mass fermions

We present our results for the leading-order hadronic quark-connected contributions to the electron, the muon, and the tau anomalous magnetic moments obtained with four dynamical quarks. Performing the continuum limit and an analysis of systematic effects, full agreement with phenomenological results is found. To estimate the impact of omitting the quark-disconnected contributions to the hadronic vacuum polarisation we investigate them on one of the four-flavour ensembles. Additionally, the light quark contributions on the four-flavour sea are compared to the values obtained for $N_f=2$ physically light quarks. In the latter case different methods to fit the hadronic vacuum polarisation function are tested.

preprint2014arXiv

Mesoscopic behavior from microscopic Markov dynamics and its application to calcium release channels

A major challenge in biology is to understand how molecular processes determine phenotypic features. We address this fundamental problem in a class of model systems by developing a general mathematical framework that allows the calculation of mesoscopic properties from the knowledge of microscopic Markovian transition probabilities. We show how exact analytic formulae for the first and second moments of resident time distributions in mesostates can be derived from microscopic resident times and transition probabilities even for systems with a large number of microstates. We apply our formalism to models of the inositol trisphosphate receptor, which plays a key role in generating calcium signals triggering a wide variety of cellular responses. We demonstrate how experimentally accessible quantities, such as opening and closing times and the coefficient of variation of inter-spike intervals, and other, more elaborated, quantities can be analytically calculated from the underlying microscopic Markovian dynamics. A virtue of our approach is that we do not need to follow the detailed time evolution of the whole system, as we derive the relevant properties of its steady state without having to take into account the often extremely complicated transient features. We emphasize that our formulae fully agree with results obtained by stochastic simulations and approaches based on a full determination of the microscopic system's time evolution. We also illustrate how experiments can be devised to discriminate between alternative molecular models of the inositol trisphosphate receptor. The developed approach is applicable to any system described by a Markov process and, owing to the analytic nature of the resulting formulae, provides an easy way to characterize also rare events that are of particular importance to understand the intermittency properties of complex dynamic systems.

preprint2014arXiv

Nucleon observables and axial charges of other baryons using twisted mass fermions

We present results on the nucleon scalar, axial and tensor charges, as well as, on the first moments of the unpolarized, polarized and transversity parton distributions using $N_f=2$ and $N_f=2+1+1$ twisted mass fermions. These include an ensemble that yields the physical value of the ratio of the nucleon to the pion mass. Results on the axial charges of hyperons and charmed baryons are also presented for a range of pion masses including the physical one.

preprint2014arXiv

Results on the disconnected contributions for hadron structure

We present results on the disconnected contributions to three point functions entering in studies of hadron structure. We use $N_F = 2+1+1$ twisted mass fermions and give a detailed description on the results of the nucleon σ-terms, isoscalar axial charge and first moments of bare parton distributions for a range of pions masses. In addition we give the σ-terms and the computations are performed using QUDA code implemented on GPUs.

preprint2014arXiv

The temperature dependence of the chiral condensate in the Schwinger model with Matrix Product States

We present our recent results for the tensor network (TN) approach to lattice gauge theories. TN methods provide an efficient approximation for quantum many-body states. We employ TN for one dimensional systems, Matrix Product States, to investigate the 1-flavour Schwinger model. In this study, we compute the chiral condensate at finite temperature. From the continuum extrapolation, we obtain the chiral condensate in the high temperature region consistent with the analytical calculation by Sachs and Wipf.

preprint2014arXiv

Topological susceptibility from the twisted mass Dirac operator spectrum

We present results of our computation of the topological susceptibility with $N_f=2$ and $N_f=2+1+1$ flavours of maximally twisted mass fermions, using the method of spectral projectors. We perform a detailed study of the quark mass dependence and discretization effects. We make an attempt to confront our data with chiral perturbation theory and extract the chiral condensate from the quark mass dependence of the topological susceptibility. We compare the value with the results of our direct computation from the slope of the mode number. We emphasize the role of autocorrelations and the necessity of long Monte Carlo runs to obtain results with good precision. We also show our results for the spectral projector computation of the ratio of renormalization constants $Z_P/Z_S$.

preprint2013arXiv

Applicability of Quasi-Monte Carlo for lattice systems

This project investigates the applicability of quasi-Monte Carlo methods to Euclidean lattice systems in order to improve the asymptotic error scaling of observables for such theories. The error of an observable calculated by averaging over random observations generated from ordinary Monte Carlo simulations scales like $N^{-1/2}$, where $N$ is the number of observations. By means of quasi-Monte Carlo methods it is possible to improve this scaling for certain problems to $N^{-1}$, or even further if the problems are regular enough. We adapted and applied this approach to simple systems like the quantum harmonic and anharmonic oscillator and verified an improved error scaling of all investigated observables in both cases.

preprint2013arXiv

Chiral condensate from the twisted mass Dirac operator spectrum

We present the results of our computation of the chiral condensate with $N_f=2$ and $N_f=2+1+1$ flavours of maximally twisted mass fermions. The condensate is determined from the Dirac operator spectrum, applying the spectral projector method proposed by Giusti and Luscher. We use 3 lattice spacings and several quark masses at each lattice spacing to perform the chiral and continuum extrapolations. We study the effect of the dynamical strange and charm quarks by comparing our results for $N_f=2$ and $N_f=2+1+1$ dynamical flavours.

preprint2013arXiv

Computation of the chiral condensate using $N_f=2$ and $N_f=2+1+1$ dynamical flavors of twisted mass fermions

We apply the spectral projector method, recently introduced by Giusti and Lüscher, to compute the chiral condensate using $N_f=2$ and $N_f=2+1+1$ dynamical flavors of maximally twisted mass fermions. We present our results for several quark masses at three different lattice spacings which allows us to perform the chiral and continuum extrapolations. In addition we report our analysis on the $O(a)$ improvement of the chiral condensate for twisted mass fermions. We also study the effect of the dynamical strange and charm quarks by comparing our results for $N_f=2$ and $N_f=2+1+1$ dynamical flavors.

preprint2013arXiv

Computing the hadronic vacuum polarization function by analytic continuation

We propose a method to compute the hadronic vacuum polarization function on the lattice at continuous values of photon momenta bridging between the spacelike and timelike regions. We provide two independent demonstrations to show that this method leads to the desired hadronic vacuum polarization function in Minkowski spacetime. We show with the example of the leading-order QCD correction to the muon anomalous magnetic moment that this approach can provide a valuable alternative method for calculations of physical quantities where the hadronic vacuum polarization function enters.

preprint2013arXiv

Constraining a fourth generation of quarks: non-perturbative Higgs boson mass bounds

We present a non-perturbative determination of the upper and lower Higgs boson mass bounds with a heavy fourth generation of quarks from numerical lattice computations in a chirally symmetric Higgs-Yukawa model. We find that the upper bound only moderately rises with the quark mass while the lower bound increases significantly, providing additional constraints on the existence of a straight-forward fourth quark generation. We examine the stability of the lower bound under the addition of a higher dimensional operator to the scalar field potential using perturbation theory, demonstrating that it is not significantly altered for small values of the coupling of this operator. For a Higgs boson mass of $\sim125\mathrm{GeV}$ we find that the maximum value of the fourth generation quark mass is $\sim300\mathrm{GeV}$, which is already in conflict with bounds from direct searches.

preprint2013arXiv

Higgs boson mass bounds in the presence of a heavy fourth quark family

We present Higgs boson mass bounds in a lattice regularization allowing thus for non-perturbative investigations. In particular, we employ a lattice modified chiral invariant Higgs-Yukawa model using the overlap operator. We show results for the upper and lower Higgs boson mass bounds in the presence of a heavy mass-degenerate quark doublet with masses ranging up to 700 GeV. We perform infinite volume extrapolations in most cases, and examine several values of the lattice cutoff. Furthermore, we argue that the lower Higgs boson mass bound is stable with respect to the addition of higher dimensional operators to the scalar field potential. Our results have severe consequences for the phenomenology of a fourth generation of quarks if a light Higgs boson is discovered at the LHC.

preprint2013arXiv

Higgs-Yukawa model in chirally-invariant lattice field theory

Non-perturbative numerical lattice studies of the Higgs-Yukawa sector of the standard model with exact chiral symmetry are reviewed. In particular, we discuss bounds on the Higgs boson mass at the standard model top quark mass, and in the presence of heavy fermions. We present a comprehensive study of the phase structure of the theory at weak and very strong values of the Yukawa coupling as well as at non-zero temperature.

preprint2013arXiv

Lattice Monte Carlo methods for systems far from equilibrium

We present a new numerical Monte Carlo approach to determine the scaling behavior of lattice field theories far from equilibrium. The presented methods are generally applicable to systems where classical-statistical fluctuations dominate the dynamics. As an example, these methods are applied to the random-force-driven one-dimensional Burgers' equation - a model for hydrodynamic turbulence. For a self-similar forcing acting on all scales the system is driven to a nonequilibrium steady state characterized by a Kolmogorov energy spectrum. We extract correlation functions of single- and multi-point quantities and determine their scaling spectrum displaying anomalous scaling for high-order moments. Varying the external forcing we are able to tune the system continuously from equilibrium, where the fluctuations are short-range correlated, to the case where the system is strongly driven in the infrared. In the latter case the nonequilibrium scaling of small-scale fluctuations are shown to be universal.

preprint2013arXiv

Leading-order hadronic contribution to the anomalous magnetic moment of the muon from N_f=2+1+1 twisted mass fermions

We present results for the leading order QCD correction to the anomalous magnetic moment of the muon including the first two generations of quarks as dynamical degrees of freedom. Several light quark masses are examined in order to yield a controlled extrapolation to the physical pion mass. We analyse ensembles for three different lattice spacings and several volumes in order to investigate lattice artefacts and finite-size effects, respectively. We also provide preliminary results for this quantity for two flavours of mass-degenerate quarks at the physical value of the pion mass.

preprint2013arXiv

Looking at the gluon moment of the nucleon with dynamical twisted mass fermions

To understand the structure of hadrons it is important to know the PDF of their constituents, the quarks and gluons. In our work we aim to compute the first moment of the gluon PDF $\langle x \rangle_g$ for the nucleon. We follow two possible approaches in order to extract the gluon moment: the Feynman-Hellmann theorem and a direct method with smearing of the gluon operator. We present preliminary results computed on $24^3 \times 48$ lattices for the case where the Feynman-Hellman theorem is used and $32^3 \times 64$ lattices for the direct method, employing $N_f=2+1+1$ maximally twisted mass fermions.

preprint2013arXiv

Matrix Product States for Lattice Field Theories

The term Tensor Network States (TNS) refers to a number of families of states that represent different ansätze for the efficient description of the state of a quantum many-body system. Matrix Product States (MPS) are one particular case of TNS, and have become the most precise tool for the numerical study of one dimensional quantum many-body systems, as the basis of the Density Matrix Renormalization Group method. Lattice Gauge Theories (LGT), in their Hamiltonian version, offer a challenging scenario for these techniques. While the dimensions and sizes of the systems amenable to TNS studies are still far from those achievable by 4-dimensional LGT tools, Tensor Networks can be readily used for problems which more standard techniques, such as Markov chain Monte Carlo simulations, cannot easily tackle. Examples of such problems are the presence of a chemical potential or out-of-equilibrium dynamics. We have explored the performance of Matrix Product States in the case of the Schwinger model, as a widely used testbench for lattice techniques. Using finite-size, open boundary MPS, we are able to determine the low energy states of the model in a fully non-perturbative manner and to accurately extract the mass spectrum. In this work we extend the analysis to the determination of the chiral condensate, both for massless and massive fermions. The method allows for accurate finite size and continuum limit extrapolations and produces remarkably precise results, thus showing the feasibility of these techniques for gauge theory problems.

preprint2013arXiv

Nucleon transversity generalized form factors with twisted mass fermions

We present results on the nucleon tensor form factors and first moment of the transversity distribution using maximally twisted mass fermions. We analyze two Nf=2+1+1 ensembles having pion masses of 213 MeV and 373 MeV with lattice spacing a=0.064 fm and a=0.082 fm, respectively. First results using an Nf=2 ensemble of twisted mass fermions with a clover term at a physical pion mass and lattice spacing a=0.094 fm are also presented. The renormalization function for the local tensor form factors is evaluated non-perturbatively with a perturbative subtraction of O(a^2)-terms, while for the first moment of the transversity we use a perturbative estimate. Results are given in the MS-bar scheme at a scale of mu=2 GeV, and are compared with recent results obtained using different discretization schemes.

preprint2013arXiv

Properties of pseudoscalar flavour-singlet mesons from 2+1+1 twisted mass lattice QCD

We study properties of pseudoscalar flavour-singlet mesons from Wilson twisted mass lattice QCD with Nf=2+1+1 dynamical quark flavors. Results for masses are presented at three values of the lattice spacing and light quark masses corresponding to values of the pion mass from 230 MeV to 500 MeV. We briefly discuss scaling effects and the light and strange quark mass dependence of the eta mass. In addition we present an exploratory study using Osterwalder-Seiler type strange and charm valence quarks. This approach avoids some of the complications of the twisted mass heavy doublet. We present first results for matching valence and unitary actions and a comparison of statistical uncertainties.

preprint2013arXiv

Stabilizing the electroweak vacuum by higher dimensional operators in a Higgs-Yukawa model

The Higgs boson discovery at the LHC with a mass of approximately 126 GeV suggests, that the electroweak vacuum of the standard model may be metastable at very high energies. However, any new physics beyond the standard model can change this picture. We want to address this important question within a lattice Higgs-Yukawa model as the limit of the standard model (SM). In this framework we will probe the effect of a higher dimensional operator for which we take a $(ϕ^{\dagger}ϕ)^3$-term. Such a term could easily originate as a remnant of physics beyond the SM at very large scales. As a first step we investigate the phase diagram of the model including such a $(ϕ^{\dagger}ϕ)^3$ operator. Exploratory results suggest the existence of regions in parameter space where first order transitions turn to second order ones, indicating the existence of a tri-critical line. We will explore the phase structure and the consequences for the stability of the SM, both analytically by investigating the constraint effective potential in lattice perturbation theory, and by studying the system non-perturbatively using lattice simulations.

preprint2013arXiv

The phase structure of a chirally-invariant Higgs-Yukawa model

We present new results of our ongoing project on the investigation of the phase structure of the Higgs-Yukawa model at small and large bare Yukawa couplings. The critical exponents of the second order bulk phase transitions of this model are determined from finite-size analyses and compared to the pure O(4)-model to test for triviality and the possibility of having a non-Gaussian fixed point. In addition, we will present a first study of Higgs boson masses and fermion correlation functions.

preprint2013arXiv

Topological susceptibility from twisted mass fermions using spectral projectors

We discuss the computation of the topological susceptibility using the method of spectral projectors and dynamical twisted mass fermions. We present our analysis concerning the O(a)-improvement of the topological susceptibility and we show numerical results for Nf=2 and Nf=2+1+1 flavours, performing a study of the quark mass dependence in terms of leading order chiral perturbation theory.

preprint2013arXiv

Using analytic continuation for the hadronic vacuum polarization computation

We present two examples of applications of the analytic continuation method for computing the hadronic vacuum polarization function in space- and time-like momentum regions. These examples are the Adler function and the leading order hadronic contribution to the muon anomalous magnetic moment. We comment on the feasibility of the analytic continuation method and provide an outlook for possible further applications.

preprint2012arXiv

Investigation of the phase structure of a chirally-invariant Higgs-Yukawa model

We present new data on our ongoing project on the investigation of the phase structure of the Higgs-Yukawa model at large bare Yukawa couplings. The data presented last year are extended in terms of statistics, the number of bare Yukawa couplings at existing, and new larger volumes. In addition, this study is extended by a finite temperature project at the physical top quark mass m_t =175 GeV and a hypothetical fourth generation top quark with a mass of m_t' =700 GeV .

preprint2012arXiv

Leading-order hadronic contributions to a_μ and α_{QED} from N_f=2+1+1 twisted mass fermions

We present the first four-flavour lattice calculation of the leading-order hadronic vacuum-polarisation contribution to the anomalous magnetic moment of the muon, a_μ^{hvp}, and the hadronic running of the QED coupling constant, Δα_{QED}^{hvp} (Q^2). In the heavy sector a mixed-action setup is employed. The bare quark masses are determined from matching the K- and D-meson masses to their physical values. Several light quark masses are used in order to yield a controlled extrapolation to the physical pion mass by utilising a recently proposed improved method. We demonstrate that this method also works in the four-flavour case.

preprint2012arXiv

Non-perturbative renormalization in coordinate space for $N_f=2$ maximally twisted mass fermions with tree-level Symanzik improved gauge action

We present results of a lattice QCD application of a coordinate space renormalization scheme for the extraction of renormalization constants for flavour non-singlet bilinear quark operators. The method consists in the analysis of the small-distance behaviour of correlation functions in Euclidean space and has several theoretical and practical advantages, in particular: it is gauge invariant, easy to implement and has relatively low computational cost. The values of renormalization constants in the X-space scheme can be converted to the MSbar scheme via 4-loop continuum perturbative formulae. Our results for $N_f=2$ maximally twisted mass fermions with tree-level Symanzik improved gauge action are compared to the ones from the RI-MOM scheme and show full agreement with this method.

preprint2012arXiv

Nonperturbative QCD corrections to electroweak observables

Nonperturbative QCD corrections are important to many low-energy electroweak observables, for example the muon magnetic moment. However, hadronic corrections also play a significant role at much higher energies due to their impact on the running of standard model parameters, such as the electromagnetic coupling. Currently, these hadronic contributions are accounted for by a combination of experimental measurements, effective field theory techniques and phenomenological modeling but ideally should be calculated from first principles. Recent developments indicate that many of the most important hadronic corrections may be feasibly calculated using lattice QCD methods. To illustrate this, we will examine the lattice computation of the leading-order QCD corrections to the muon magnetic moment, paying particular attention to a recently developed method but also reviewing the results from other calculations. We will then continue with several examples that demonstrate the potential impact of the new approach: the leading-order corrections to the electron and tau magnetic moments, the running of the electromagnetic coupling, and a class of the next-to-leading-order corrections for the muon magnetic moment. Along the way, we will mention applications to the Adler function, which can be used to determine the strong coupling constant, and QCD corrections to muonic-hydrogen.

preprint2012arXiv

Overlap valence quarks on a twisted mass sea: a case study for mixed action Lattice QCD

We discuss a Lattice QCD mixed action investigation employing Wilson maximally twisted mass sea and overlap valence fermions. Using four values of the lattice spacing, we demonstrate that the overlap Dirac operator assumes a point-like locality in the continuum limit. We also show that by adopting suitable matching conditions for the sea and valence theories a consistent continuum limit for the pion decay constant and light baryon masses can be obtained. Finally, we confront results for sea-valence mixed meson masses and the valence scalar correlator with corresponding expressions of chiral perturbation theory. This allows us to extract low energy constants of mixed action chiral perturbation which characterize the strength of unitarity violations in our mixed action setup.

preprint2012arXiv

Sigma terms and strangeness content of the nucleon with $N_f=2+1+1$ twisted mass fermions

We study the nucleon matrix elements of the quark scalar-density operator using maximally twisted mass fermions with dynamical light ($u$,$d$), strange and charm degrees of freedom. We demonstrate that in this setup the nucleon matrix elements of the light and strange quark densities can be obtained with good statistical accuracy, while for the charm quark counterpart only a bound can be provided. The present calculation which is performed at only one value of the lattice spacing and pion mass serves as a benchmark for a future more systematic computation of the scalar quark content of the nucleon.

preprint2012arXiv

The Higgs boson resonance width from a chiral Higgs-Yukawa model on the lattice

The Higgs boson is a central part of the electroweak theory and is crucial to generate masses for quarks, leptons and the weak gauge bosons. We use a 4-dimensional Euclidean lattice formulation of the Higgs-Yukawa sector of the electroweak model to compute physical quantities in the path integral approach which is evaluated by means of Monte Carlo simulations thus allowing for fully non perturbative calculations. The chiral symmetry of the model is incorporated by using the Neuberger overlap Dirac operator. The here considered Higgs-Yukawa model does not involve the weak gauge bosons and furthermore, only a degenerate doublet of top- and bottom quarks are incorporated. The goal of this work is to study the resonance properties of the Higgs boson and its sensitivity to the strength of the quartic self coupling.

preprint2011arXiv

Anomalous scaling in the random-force-driven Burgers equation: A Monte Carlo study

We present a new approach to determine the small-scale statistical behavior of hydrodynamic turbulence by means of lattice simulations. Using the functional integral representation of the random-force-driven Burgers equation we show that high-order moments of velocity differences satisfy anomalous scaling. The general applicability of Monte Carlo methods provides the opportunity to study also other systems of interest within this framework.

preprint2011arXiv

Anomalous scaling in the random-force-driven Burgers equation: A Monte Carlo study

We present a new approach to determine numerically the statistical behavior of small-scale structures in hydrodynamic turbulence. Starting from the functional integral representation of the random-force-driven Burgers equation we show that Monte Carlo simulations allow us to determine the anomalous scaling of high-order moments of velocity differences. Given the general applicability of Monte Carlo methods, this opens up the possibility to address also other systems relevant to turbulence within this framework.

preprint2011arXiv

Charm Current-Current Correlators in Twisted Mass Lattice QCD

The current correlator method has been shown to be a practical tool to extract the charm quark mass and strong coupling constant from Lattice QCD data as an alternative to the sum rule approach using experimental electron-positron annihilation cross section data. We report on the progress of an investigation of charm current-current correlators in N_f=2 Twisted Mass Lattice QCD. Upon determining the temporal moments of the current correlators we compare to the low-energy expansion of the moments in perturbative QCD and calculate the charm quark mass and strong coupling constant and in case of the vector current correlator directly compare both methods.

preprint2011arXiv

Cosmological Phase transitions from Lattice Field Theory

In this proceedings contribution we discuss the fate of the electroweak and the quantum chromodynamics phase transitions relevant for the early stage of the universe at non-zero temperature. These phase transitions are related to the Higgs mechanism and the breaking of chiral symmetry, respectively. We will review that non-perturbative lattice field theory simulations show that these phase transitions actually do not occur in nature and that physical observables show a completely smooth behaviour as a function of the temperature.

preprint2011arXiv

Dark matter Search and the Scalar Quark Contents of the Nucleon

We present lattice QCD simulation results from the European Twisted Mass Collaboration (ETMC) for the light, strange and charm quark contents of the nucleon. These quantities are important ingredients to estimate the cross-section for the detection of WIMPs as Dark Matter candidates. By employing a particular lattice QCD formulation, i.e. twisted mass fermions, accurate results of the light and strange scalar contents of the nucleon can be obtained. In addition, we provide a bound for the charm quark content of the nucleon.

preprint2011arXiv

Excited State Effects in Nucleon Matrix Element Calculations

We perform a high-statistics precision calculation of nucleon matrix elements using an open sink method allowing us to explore a wide range of sink-source time separations. In this way the influence of excited states of nucleon matrix elements can be studied. As particular examples we present results for the nucleon axial charge $g_A$ and for the first moment of the isovector unpolarized parton distribution $<x>_{u-d}$. In addition, we report on preliminary results using the generalized eigenvalue method for nucleon matrix elements. All calculations are performed using $N_f=2+1+1$ maximally twisted mass Wilson fermions.

preprint2011arXiv

Higgs boson resonance parameters and the finite temperature phase transition in a chirally invariant Higgs-Yukawa model

We study a chirally invariant Higgs-Yukawa model regulated on a space-time lattice. We calculate Higgs boson resonance parameters and mass bounds for various values of the mass of the degenerate fermion doublet. Also, first results on the phase transition temperature are presented. In general, this model may be relevant for BSM scenarios with a heavy fourth generation of quarks.

preprint2011arXiv

Lambda_msbar from the static potential for QCD with n_f=2 dynamical quark flavors

We determine Lambda_msbar for QCD with n_f=2 dynamical quark flavors by fitting the Q-Q-bar static potential known analytically in the perturbative regime up to terms of O(alpha_s^4) and ~alpha_s^4 ln(alpha_s) to corresponding results obtained from lattice simulations. This has become possible, due to recent advances in both perturbative calculations, namely the determination and publication of the last missing contribution to the Q-Q-bar static potential at O(alpha_s^4), and lattice simulations with n_f=2 dynamical quark flavors performed at the rather fine lattice spacing of a~0.042 fm. Imposing conservative error estimates we obtain Lambda_msbar=315(30) MeV.

preprint2011arXiv

Leading-order hadronic contribution to g-2 from lattice QCD

We calculate the leading-order hadronic correction to the anomalous magnetic moments of each of the three charged leptons in the Standard Model: the electron, muon and tau. Working in two-flavor lattice QCD, we address essentially all sources of systematic error: lattice artifacts, finite-size effects, quark-mass extrapolation, momentum extrapolation and disconnected diagrams. The most significant remaining systematic error, the exclusion of the strange and charm quark contributions, will be addressed in our four-flavor calculation. We achieve a statistical accuracy of 2% or better for the physical values for each of the three leptons and the systematic errors are at most comparable.

preprint2011arXiv

Nucleon matrix elements with $N_f=2+1+1$ maximally twisted fermions

We present the first lattice calculation of nucleon matrix elements using four dynamical flavors. We use the $N_f=2+1+1$ maximally twisted mass formulation. The renormalization is performed non-perturbatively in the RI$^\prime$-MOM scheme and results are given for the vector and axial vector operators with up to one-derivative. Our calculation of the average momentum of the unpolarized non-singlet parton distribution is presented and compared to our previous results obtained from the $N_f=2$ case.

preprint2011arXiv

Precision Study of Excited State Effects in Nucleon Matrix Elements

We present a dedicated analysis of the influence of excited states on the calculation of nucleon matrix elements. This calculation is performed at a fixed value of the lattice spacing, volume and pion mass that are typical of contemporary lattice computations. We focus on the nucleon axial charge, g_A, for which we use about 7,500 measurements, and on the average momentum of the unpolarized isovector parton distribution, <x>_{u-d}, for which we use about 23,000 measurements. All computations are done employing N_f=2+1+1 maximally-twisted-mass Wilson fermions and using non-perturbatively calculated renormalization factors. Excited state effects are shown to be negligible for g_A whereas they lead to an O(10%) downward shift for <x>_{u-d}.

preprint2011arXiv

Resonance Parameters of the rho-Meson from Lattice QCD

We perform a non-perturbative lattice calculation of the P-wave pion-pion scattering phase in the rho-meson decay channel using two flavors of maximally twisted mass fermions at pion masses ranging from 480 MeV to 290 MeV. Making use of finite-size methods, we evaluate the pion-pion scattering phase in the center-of-mass frame and two moving frames. Applying an effective range formula, we find a good description of our results for the scattering phase as a function of the energy covering the resonance region. This allows us to extract the rho-meson mass and decay width and to study their quark mass dependence.

preprint2011arXiv

Study of the Higgs-Yukawa theory in the strong-Yukawa coupling regime

In this article, we present an ongoing lattice study of the Higgs-Yukawa model, in the regime of strong-Yukawa coupling, using overlap fermions. We investigated the phase structure in this regime by computing the Higgs vacuum expectation value, and by exploring the finite-size scaling behaviour of the susceptibility corresponding to the magnetisation. Our preliminary results indicate the existence of a second-order phase transition when the Yukawa coupling becomes large enough, at which the Higgs vacuum expectation value vanishes and the susceptibility diverges.

preprint2011arXiv

Two-flavor QCD correction to lepton magnetic moments at leading-order in the electromagnetic coupling

We present a reliable nonperturbative calculation of the QCD correction, at leading order in the electromagnetic coupling, to the anomalous magnetic moment of the electron, muon, and tau leptons using two-flavor lattice QCD. We use multiple lattice spacings, multiple volumes, and a broad range of quark masses to control the continuum, infinite-volume, and chiral limits. We examine the impact of the commonly ignored disconnected diagrams and introduce a modification to the previously used method that results in a well-controlled lattice calculation. We obtain 1.513(43) 10^(-12), 5.72(16) 10^(-8), and 2.650(54) 10^(-6) for the leading-order two-flavor QCD correction to the anomalous magnetic moment of the electron, muon, and tau, respectively, each accurate to better than 3%.

preprint2010arXiv

Comparing topological charge definitions using topology fixing actions

We investigate both the hyperbolic action and the determinant ratio action designed to fix the topological charge on the lattice. We show to what extent topology is fixed depending on the parameters of these actions, keeping the physical situation fixed. At the same time the agreement between different definitions of topological charge - the field theoretic and the index definition - is directly correlated to the degree topology is fixed. Moreover, it turns out that the two definitions agree very well. We also study finite volume effects arising in the static potential and related quantities due to topology fixing.

preprint2010arXiv

Computing K and D meson masses with N_f = 2+1+1 twisted mass lattice QCD

We discuss the computation of the mass of the K and D mesons within the framework of N_f = 2+1+1 twisted mass lattice QCD from a technical point of view. These quantities are essential, already at the level of generating gauge configurations, being obvious candidates to tune the strange and charm quark masses to their physical values. In particular, we address the problems related to the twisted mass flavor and parity symmetry breaking, which arise when considering a non-degenerate (c,s) doublet. We propose and verify the consistency of three methods to extract the K and D meson masses in this framework.

preprint2010arXiv

Continuum Limit of Overlap Valence Quarks on a Twisted Mass Sea

We study a lattice QCD mixed action with overlap valence quarks on two flavours of Wilson maximally twisted mass sea quarks. Employing three different matching conditions to relate both actions to each other, we investigate the continuum limit by using three values of the lattice spacing ranging from 0.05 fm to 0.08 fm. A particular emphasis is put on the effect on physical observables of the topological zero modes appearing in the valence overlap operator. We estimate the region of parameter space where the contribution from these zero modes is sufficiently small such that their effects can be safely controlled and a restoration of unitarity of the mixed action in the continuum limit is reached.

preprint2010arXiv

Effects of a potential fourth fermion generation on the upper and lower Higgs boson mass bounds

We study the effect of a potential fourth fermion generation on the upper and lower Higgs boson mass bounds. This investigation is based on the numerical evaluation of a chirally invariant lattice Higgs-Yukawa model emulating the same Higgs-fermion coupling structure as in the Higgs sector of the electroweak Standard Model. In particular, the considered model obeys a Ginsparg-Wilson version of the underlying ${SU}(2)_L\times {U}(1)_Y$ symmetry, being a global symmetry here due to the neglection of gauge fields in this model. We present our results on the modification of the upper and lower Higgs boson mass bounds induced by the presence of a hypothetical very heavy fourth quark doublet. Finally, we compare these findings to the standard scenario of three fermion generations.

preprint2010arXiv

Leading order hadronic contribution to g-2 from twisted mass QCD

We calculate the leading order hadronic contribution to the muon anomalous magnetic moment using twisted mass lattice QCD. The pion masses range from 330 MeV to 650 MeV. We use two lattice spacings, a=0.079 fm and 0.063 fm, to study lattice artifacts. Finite-size effects are studied for two values of the pion mass, and we calculate the disconnected contributions for four ensembles. Particular attention is paid to the dominant contributions of the vector mesons, both phenomenologically and from our lattice calculation.

preprint2010arXiv

Light hadrons from lattice QCD with light (u,d), strange and charm dynamical quarks

We present results of lattice QCD simulations with mass-degenerate up and down and mass-split strange and charm (N_f = 2+1+1) dynamical quarks using Wilson twisted mass fermions at maximal twist. The tuning of the strange and charm quark masses is performed at two values of the lattice spacing a~0.078 fm and a~0.086 fm with lattice sizes ranging from L~1.9 fm to L~2.8 fm. We measure with high statistical precision the light pseudoscalar mass m_PS and decay constant f_PS in a range 270 < m_PS < 510 MeV and determine the low energy parameters f_0, l_3 and l_4 of SU(2) chiral perturbation theory. We use the two values of the lattice spacing, several lattice sizes as well as different values of the light, strange and charm quark masses to explore the systematic effects. A first study of discretisation effects in light-quark observables and a comparison to N_f=2 results are performed.

preprint2010arXiv

Low lying baryon spectrum with $N_f = 2+1+1$ dynamical twisted quarks

We present first results on the octet and decuplet strange baryon spectrum with $N_f=2+1+1$ twisted mass quarks. We use an Osterwalder Seiler valence strange quark with a mass tuned to the kaon and compare the results with those obtained in the unitary setup. This comparison allows to perform a first study of the lattice artefacts introduced by the mixed action approach. We investigate the effect of the strange and charm quarks in the sea by using two lattice spacings and comparing with preceding $N_f = 2$ twisted mass fermion calculations.

preprint2010arXiv

Overlap Valence Quarks on a Twisted Mass Sea

We present the results of an investigation of a mixed action approach of overlap valence and maximally twisted mass sea quarks. Employing a particular matching condition on the pion mass, we analyze the continuum limit scaling of the pion decay constant and the role of chiral zero modes of the overlap operator in this process. We employ gauge field configurations generated by the European Twisted Mass Collaboration with linear lattice size $L$ ranging from 1.3 to 1.9 fm. The continuum limit is taken at a fixed value of $L=1.3$ fm, employing three values of the lattice spacing and two values of the pion mass constructed from sea quarks only.

preprint2010arXiv

Pseudoscalar decay constants from N_f=2+1+1 twisted mass lattice QCD

We present first results for the pseudoscalar decay constants $f_K$, $f_D$ and $f_{D_s}$ from lattice QCD with N_f=2+1+1 flavours of dynamical quarks. The lattice simulations have been performed by the European Twisted Mass collaboration (ETMC) using maximally twisted mass quarks. For the pseudoscalar decay constants we follow a mixed action approach by using so called Osterwalder-Seiler fermions in the valence sector for strange and charm quarks. The data for two values of the lattice spacing and several values of the up/down quark mass is analysed using chiral perturbation theory.

preprint2010arXiv

The pi+ pi+ scattering length from maximally twisted mass lattice QCD

We calculate the s-wave pion-pion scattering length in the isospin I=2 channel in lattice QCD for pion masses ranging from 270 Mev to 485 Mev using two flavors of maximally twisted mass fermions at a lattice spacing of 0.086 fm. Additionally, we check for lattice artifacts with one calculation at a finer lattice spacing of 0.067 fm. We use chiral perturbation theory at next-to-leading order to extrapolate our results. At the physical pion mass, we find m_pi a_pipi(I=2)=-0.04385(28)(38) for the scattering length, where the first error is statistical and the second is our estimate of several systematic effects.

preprint2009arXiv

Chirally rotated Schroedinger functional: non-perturbative tuning in the quenched approximation

The use of chirally rotated boundary conditions provides a formulation of the Schroedinger functional that is compatible with automatic O(a)-improvement of Wilson fermions in the bulk. The elimination of bulk O(a) terms requires the non-perturbative tuning of the critical mass and one additional boundary counterterm. We present the results of such a tuning in the quenched approximation at three values of the renormalised gauge coupling and for a range of lattice spacings.

preprint2009arXiv

Continuum-limit scaling of overlap fermions as valence quarks

We present the results of a mixed action approach, employing dynamical twisted mass fermions in the sea sector and overlap valence fermions, with the aim of testing the continuum limit scaling behaviour of physical quantities, taking the pion decay constant as an example. To render the computations practical, we impose for this purpose a fixed finite volume with lattice size $L\approx1.3$ fm. We also briefly review the techniques we have used to deal with overlap fermions.

preprint2009arXiv

Scattering from finite size methods in lattice QCD

Using two flavors of maximally twisted mass fermions, we calculate the S-wave pion-pion scattering length in the isospin I=2 channel and the P-wave pion-pion scattering phase in the isospin I=1 channel. In the former channel, the lattice calculations are performed at pion masses ranging from 270 MeV to 485 MeV. We use chiral perturbation theory at next-to-leading order to extrapolate our results. At the physical pion mass, we find m_pi a_pipi(I=2)=-0.04385(28)(38) for the scattering length. In the latter channel, the calculation is currently performed at a single pion mass of 391 MeV. Making use of finite size methods, we evaluate the scattering phase in both the center of mass frame and the moving frame. The effective range formula is employed to fit our results, from which the rho resonance mass and decay width are evaluated.

preprint2008arXiv

Higgs mass bounds from a chirally invariant lattice Higgs-Yukawa model with overlap fermions

We study the parameter dependence of the Higgs mass in a chirally invariant lattice Higgs-Yukawa model emulating the same Higgs-fermion coupling structure as in the Higgs sector of the electroweak Standard Model. Eventually, the aim is to establish upper and lower Higgs mass bounds. Here we present our preliminary results on the lower Higgs mass bound at several selected values for the cutoff and give a brief outlook towards the upper Higgs mass bound.

preprint2006arXiv

Status Report on ILDG activities

A status report about the International Lattice Data Grid (ILDG) is given. Different countries participating in the ILDG have created regional lattice data grid solutions that are implemented, working and used. The remaining task and the focus of present activities is the development of the interoperability of these regional grids. A first, successful step in this direction is a metadata catalogue service which is already working interoperable.

preprint2001arXiv

Non-perturbative renormalization of the quark condensate in Ginsparg-Wilson regularizations

We present a method to compute non-perturbatively the renormalization constant of the scalar density for Ginsparg-Wilson fermions. It relies on chiral symmetry and is based on a matching of renormalization group invariant masses at fixed pseudoscalar meson mass, making use of results previously obtained by the ALPHA Collaboration for O(a)-improved Wilson fermions. Our approach is quite general and enables the renormalization of scalar and pseudoscalar densities in lattice regularizations that preserve chiral symmetry and of fermion masses in any regularization. As an application we compute the non-perturbative factor which relates the renormalization group invariant quark condensate to its bare counterpart, obtained with overlap fermions at beta=5.85 in the quenched approximation.

preprint2001arXiv

Scalar condensate and light quark masses from overlap fermions

We have studied pseudoscalar correlation functions computed using the overlap operator. Within the accuracy of our calculation we find that the quark mass dependence agrees with the prediction of lowest-order Chiral Perturbation Theory (ChPT) for quark masses in the range of m ~ m_s/2-2m_s. We present the results of an analysis which assumes lowest-order ChPT to be valid to extract the low-energy constants Sigma and f_P, as well as the strange quark mass. Non-perturbative renormalization is implemented via a matching procedure with data obtained using Wilson fermions in the Schroedinger functional set-up. We find that the scalar condensate computed here agrees with the one obtained previously through a finite-size scaling analysis.

preprint1997arXiv

Liapunov Exponents and the Reversibility of Molecular Dynamics Algorithms

We study the phenomenon of lack of reversibility in molecular dynamics algorithms for the case of Wilson's lattice QCD. We demonstrate that the classical equations of motion that are employed in these algorithms are chaotic in nature. The leading Liapunov exponent is determined in a range of coupling parameters. We give a quantitative estimate of the consequences of the breakdown of reversibility due to round-off errors.

preprint1994arXiv

Domain Wall Fermions and Chiral Gauge Theories

We review the status of the domain wall fermion approach to construct chiral gauge theories on the lattice. In this model an extra, fifth dimension is added and our 4-dimensional world lives on a domainwall induced by a soliton shaped mass defect that depends on the extra dimension only. We demonstrate that the domain wall model gives the correct anomaly structure when external gauge fields are used. We discuss two ways of adding dynamical gauge fields aiming at a lattice regularized chiral gauge theory. An approach is presented to keep the lattice infinite by regarding the fifth direction as the time-axis of a 4-dimensional Hamiltonian. Finally, a prospect to use domain wall fermions for simulating QCD is given.

preprint1993arXiv

The Higgs Model with a Complex Ghost Pair

A higher derivative term is introduced in the kinetic energy of the Higgs Lagrangian in the minimal Standard Model. A logically consistent and {\it finite} field theory is obtained when some excitations of the Higgs field are quantized with indefinite metric in the Hilbert space. The Landau ghost phenomenon of the conventional triviality problem is replaced by the state vectors of a complex ghost pair at a finite mass scale with observable physical consequences. It is shown that the ghost states exhibit unusual resonance properties and correspond to a complex conjugate pair of Pauli-Villars regulator masses in the euclidean path integral formulation of the theory. An argument is given that microscopic acausality effects associated with the ghost pair remain undetectable in scattering processes with realistic wave packects, and the S-matrix should exhibit unitarity in the observable sector of the Hilbert space. Part One of Extended UCSD-PTH 92-40