Source author record

S. D. Katz

S. D. Katz 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

46works
10topics
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

46 published item(s)

preprint2022arXiv

QED and strong isospin corrections in the hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon

Recently, the Budapest-Marseille-Wuppertal collaboration achieved sub-percent precision in the evaluation of the lowest-order hadronic vacuum polarization contribution to the muon $g_μ-2$ (arXiv:hep-lat/2002.12347v3). At this level of precision, isospin-symmetric QCD is not sufficient. In this contribution we review how QED and strong-isospin-breaking effects have been included in our work. Isospin breaking is implemented by expanding the relevant correlation functions to second order in the electric charge $e$ and to first order in $m_u-m_d$. The correction terms are then computed using isospin-symmetric configurations. The choice of this approach allows us to better distribute the available computing resources among the various contributions.

preprint2020arXiv

Why integral equations should be used instead of differential equations to describe the dynamics of epidemics

It is of vital importance to understand and track the dynamics of rapidly unfolding epidemics. The health and economic consequences of the current COVID-19 pandemic provide a poignant case. Here we point out that since they are based on differential equations, the most widely used models of epidemic spread are plagued by an approximation that is not justified in the case of the current COVID-19 pandemic. Taking the example of data from New York City, we show that currently used models significantly underestimate the initial basic reproduction number ($R_0$). The correct description, based on integral equations, can be implemented in most of the reported models and it much more accurately accounts for the dynamics of the epidemic after sharp changes in $R_0$ due to restrictive public congregation measures. It also provides a novel way to determine the incubation period, and most importantly, as we demonstrate for several countries, this method allows an accurate monitoring of $R_0$ and thus a fine-tuning of any restrictive measures. Integral equation based models do not only provide the conceptually correct description, they also have more predictive power than differential equation based models, therefore we do not see any reason for using the latter.

preprint2016arXiv

Lattice computation of the nucleon scalar quark contents at the physical point

We present a QCD calculation of the $u$, $d$ and $s$ scalar quark contents of nucleons based on $47$ lattice ensembles with $N_f = 2+1$ dynamical sea quarks, $5$ lattice spacings down to $0.054\,\text{fm}$, lattice sizes up to $6\,\text{fm}$ and pion masses down to $120\,\text{MeV}$. Using the Feynman-Hellmann theorem, we obtain $f^N_{ud} = 0.0405(40)(35)$ and $f^N_s = 0.113(45)(40)$, which translates into $σ_{πN}=38(3)(3)\,\text{MeV}$, $σ_{sN}=105(41)(37)\,\text{MeV}$ and $y_N=0.20(8)(8)$ for the sigma terms and the related ratio, where the first errors are statistical and the second are systematic. Using isospin relations, we also compute the individual up and down quark contents of the proton and neutron (results in the main text).

preprint2016arXiv

Lattice QCD for Cosmology

We present a full result for the equation of state (EoS) in 2+1+1 (up/down, strange and charm quarks are present) flavour lattice QCD. We extend this analysis and give the equation of state in 2+1+1+1 flavour QCD. In order to describe the evolution of the universe from temperatures several hundreds of GeV to several tens of MeV we also include the known effects of the electroweak theory and give the effective degree of freedoms. As another application of lattice QCD we calculate the topological susceptibility (chi) up to the few GeV temperature region. These two results, EoS and chi, can be used to predict the dark matter axion's mass in the post-inflation scenario and/or give the relationship between the axion's mass and the universal axionic angle, which acts as a initial condition of our universe.

preprint2016arXiv

Sigma term and strangeness content of octet baryons

By using lattice QCD computations we determine the sigma terms and strangeness content of all octet baryons by means of an application of the Hellmann-Feynman theorem. In addition to polynomial and rational expressions for the quark mass dependence of octet members, we use SU(3) covariant baryon chiral perturbation theory to perform the extrapolation to the physical up and down quark masses. Our N_f=2+1 lattice ensembles include pion masses down to about 190 MeV in large volumes (M_πL > 4), and three values of the lattice spacing. Our main results are the nucleon sigma term σ_{πN} = 39(4)(^{+18}_{-7}) and the strangeness content y_{N} = 0.20(7)(^{+13}_{-17}). Under the assumption of validity of covariant baryon χPT in our range of masses one finds y_{N} = 0.276(77)(^{+90}_{-62}).

preprint2016arXiv

Towards the QCD phase diagram from analytical continuation

We calculate the QCD cross-over temperature, the equation of state and fluctuations of conserved charges at finite density by analytical continuation from imaginary to real chemical potentials. Our calculations are based on new continuum extrapolated lattice simulations using the 4stout staggered actions with a lattice resolution up to $N_t=16$. The simulation parameters are tuned such that the strangeness neutrality is maintained, as it is in heavy ion collisions.

preprint2015arXiv

Ab initio calculation of the neutron-proton mass difference

The existence and stability of atoms rely on the fact that neutrons are more massive than protons. The measured mass difference is only 0.14\% of the average of the two masses. A slightly smaller or larger value would have led to a dramatically different universe. Here, we show that this difference results from the competition between electromagnetic and mass isospin breaking effects. We performed lattice quantum-chromodynamics and quantum-electrodynamics computations with four nondegenerate Wilson fermion flavors and computed the neutron-proton mass-splitting with an accuracy of $300$ kilo-electron volts, which is greater than $0$ by $5$ standard deviations. We also determine the splittings in the $Σ$, $Ξ$, $D$ and $Ξ_{cc}$ isospin multiplets, exceeding in some cases the precision of experimental measurements.

preprint2015arXiv

Axion cosmology, lattice QCD and the dilute instanton gas

Axions are one of the most attractive dark matter candidates. The evolution of their number density in the early universe can be determined by calculating the topological susceptibility $χ(T)$ of QCD as a function of the temperature. Lattice QCD provides an ab initio technique to carry out such a calculation. A full result needs two ingredients: physical quark masses and a controlled continuum extrapolation from non-vanishing to zero lattice spacings. We determine $χ(T)$ in the quenched framework (infinitely large quark masses) and extrapolate its values to the continuum limit. The results are compared with the prediction of the dilute instanton gas approximation (DIGA). A nice agreement is found for the temperature dependence, whereas the overall normalization of the DIGA result still differs from the non-perturbative continuum extrapolated lattice results by a factor of order ten. We discuss the consequences of our findings for the prediction of the amount of axion dark matter.

preprint2015arXiv

Complex Langevin dynamics for dynamical QCD at nonzero chemical potential: a comparison with multi-parameter reweighting

We study lattice QCD at non-vanishing chemical potential using the complex Langevin equation. We compare the results with multi-parameter reweighting both from $μ=0$ and phase quenched ensembles. We find a good agreement for lattice spacings below $\approx$0.15 fm. On coarser lattices the complex Langevin approach breaks down. Four flavors of staggered fermions are used on $N_t=4, 6$ and 8 lattices. For one ensemble we also use two flavors to investigate the effects of rooting.

preprint2015arXiv

Fluctuations and correlations in high temperature QCD

We calculate second- and fourth-order cumulants of conserved charges in a temperature range stretching from the QCD transition region towards the realm of (resummed) perturbation theory. We perform lattice simulations with staggered quarks; the continuum extrapolation is based on $N_t=10\dots24$ in the crossover-region and $N_t=8\dots16$ at higher temperatures. We find that the Hadron Resonance Gas model predictions describe the lattice data rather well in the confined phase. At high temperatures (above $\sim$250 MeV) we find agreement with the three-loop Hard Thermal Loop results.

preprint2015arXiv

Quantum electrodynamics in finite volume and nonrelativistic effective field theories

Electromagnetic effects are increasingly being accounted for in lattice quantum chromodynamics computations. Because of their long-range nature, they lead to large finite-size effects over which it is important to gain analytical control. Nonrelativistic effective field theories provide an efficient tool to describe these effects. Here we argue that some care has to be taken when applying these methods to quantum electrodynamics in a finite volume.

preprint2015arXiv

The QCD phase diagram from analytic continuation

We present the crossover line between the quark gluon plasma and the hadron gas phases for small real chemical potentials. First we determine the effect of imaginary values of the chemical potential on the transition temperature using lattice QCD simulations. Then we use various formulas to perform an analytic continuation to real values of the baryo-chemical potential. Our data set maintains strangeness neutrality to match the conditions of heavy ion physics. The systematic errors are under control up to $μ_B\approx 300$ MeV. For the curvature of the transition line we find that there is an approximate agreement between values from three different observables: the chiral susceptibility, chiral condensate and strange quark susceptibility. The continuum extrapolation is based on $N_t=$ 10, 12 and 16 lattices. By combining the analysis for these three observables we find, for the curvature, the value $κ= 0.0149 \pm 0.0021$.

preprint2014arXiv

Freeze-out parameters from electric charge and baryon number fluctuations: is there consistency?

Recent results for moments of multiplicity distributions of net-protons and net-electric charge from the STAR collaboration are compared to lattice QCD results for higher order fluctuations of baryon number and electric charge by the Wuppertal-Budapest collaboration, with the purpose of extracting the freeze-out temperature and chemical potential. All lattice simulations are performed for a system of 2+1 dynamical quark flavors, at the physical mass for light and strange quarks; all results are continuum extrapolated. We show that it is possible to extract an upper value for the freeze-out temperature, as well as precise baryo-chemical potential values corresponding to the four highest collision energies of the experimental beam energy scan. Consistency between the freeze-out parameters obtained from baryon number and electric charge fluctuations is found. The freeze-out chemical potentials are now in agreement with the statistical hadronization model.

preprint2014arXiv

Isospin splittings in the light baryon octet from lattice QCD and QED

While electromagnetic and up-down quark mass difference effects on octet baryon masses are very small, they have important consequences. The stability of the hydrogen atom against beta decay is a prominent example. Here we include these effects by adding them to valence quarks in a lattice QCD calculation based on $N_f=2+1$ simulations with 5 lattice spacings down to 0.054 fm, lattice sizes up to 6 fm and average up-down quark masses all the way down to their physical value. This allows us to gain control over all systematic errors, except for the one associated with neglecting electromagnetism in the sea. We compute the octet baryon isomultiplet mass splittings, as well as the individual contributions from electromagnetism and the up-down quark mass difference. Our results for the total splittings are in good agreement with experiment.

preprint2014arXiv

Local CP-violation and electric charge separation by magnetic fields from lattice QCD

We study local CP-violation on the lattice by measuring the local correlation between the topological charge density and the electric dipole moment of quarks, induced by a constant external magnetic field. This correlator is found to increase linearly with the external field, with the coefficient of proportionality depending only weakly on temperature. Results are obtained on lattices with various spacings, and are extrapolated to the continuum limit after the renormalization of the observables is carried out. This renormalization utilizes the gradient flow for the quark and gluon fields. Our findings suggest that the strength of local CP-violation in QCD with physical quark masses is about an order of magnitude smaller than a model prediction based on nearly massless quarks in domains of constant gluon backgrounds with topological charge. We also show numerical evidence that the observed local CP-violation correlates with spatially extended electric dipole structures in the QCD vacuum.

preprint2014arXiv

The QCD equation of state in background magnetic fields

We determine the equation of state of 2+1-flavor QCD with physical quark masses, in the presence of a constant (electro)magnetic background field on the lattice. To determine the free energy at nonzero magnetic fields we develop a new method, which is based on an integral over the quark masses up to asymptotically large values where the effect of the magnetic field can be neglected. The method is compared to other approaches in the literature and found to be advantageous for the determination of the equation of state up to large magnetic fields. Thermodynamic observables including the longitudinal and transverse pressure, magnetization, energy density, entropy density and interaction measure are presented for a wide range of temperatures and magnetic fields, and provided in ancillary files. The behavior of these observables confirms our previous result that the transition temperature is reduced by the magnetic field. We calculate the magnetic susceptibility and permeability, verifying that the thermal QCD medium is paramagnetic around and above the transition temperature, while we also find evidence for weak diamagnetism at low temperatures.

preprint2013arXiv

Freeze-out parameters from continuum extrapolated lattice data

We present continuum extrapolated lattice results for the higher order fluctuations of conserved charges in high temperature Quantum Chromodynamics. Through the matching of the grand canonical ensemble on the lattice to the net charge and net baryon distribution realized in heavy ion experiments the temperature and the chemical potential may be estimated at the time of chemical freeze-out

preprint2013arXiv

Freeze-out parameters: lattice meets experiment

We present our results for ratios of higher order fluctuations of electric charge as functions of the temperature. These results are obtained in a system of 2+1 quark flavors at physical quark masses and continuum extrapolated. We compare them to preliminary data on higher order moments of the net electric charge distribution from the STAR collaboration. This allows us to determine the freeze-out temperature and chemical potential from first principles. We also show continuum-extrapolated results for ratios of higher order fluctuations of baryon number. These will allow to test the consistency of the approach, by comparing them to the corresponding experimental data (once they become available) and thus extracting the freeze-out parameters in an independent way.

preprint2013arXiv

Perturbative renormalization of staggered fermion operators with stout improvement: Application to the magnetic susceptibility of QCD

We calculate the fermion propagator and the quark-antiquark Green's functions for a complete set of ultralocal fermion bilinears, ${{\cal O}_Γ}$ [$Γ$: scalar (S), pseudoscalar (P), vector (V), axial (A) and tensor (T)], using perturbation theory up to one-loop and to lowest order in the lattice spacing. We employ the staggered action for fermions and the Symanzik Improved action for gluons. From our calculations we determine the renormalization functions for the quark field and for all ultralocal taste-singlet bilinear operators. The novel aspect of our calculations is that the gluon links which appear both in the fermion action and in the definition of the bilinears have been improved by applying a stout smearing procedure up to two times, iteratively. Compared to most other improved formulations of staggered fermions, the above action, as well as the HISQ action, lead to smaller taste violating effects. The renormalization functions are presented in the RI$'$ scheme; the dependence on all stout parameters, as well as on the coupling constant, the number of colors, the lattice spacing, the gauge fixing parameter and the renormalization scale, is shown explicitly. We apply our results to a nonperturbative study of the magnetic susceptibility of QCD at zero and finite temperature. In particular, we evaluate the "tensor coefficient", $τ$, which is relevant to the anomalous magnetic moment of the muon.

preprint2013arXiv

Thermodynamic properties of QCD in external magnetic fields

We consider the effect of strong external electromagnetic fields on thermodynamic observables in QCD, through lattice simulations with 1+1+1 flavors of staggered quarks at physical quark masses. Continuum extrapolated results are presented for the light quark condensates and for their tensor polarizations, as functions of the temperature and the magnetic field. We find the light condensates to undergo inverse magnetic catalysis in the transition region, in a manner that the transition temperature decreases with growing magnetic field. We also compare the results to other approaches and lattice simulations. Furthermore, we relate the tensor polarization to the spin part of the magnetic susceptibility of the QCD vacuum, and show that this contribution is diamagnetic.

preprint2012arXiv

Anisotropy tuning with the Wilson flow

We use the Wilson flow to define the gauge anisotropy at a given physical scale. We demonstrate the use of the anisotropic flow by performing the tuning of the bare gauge anisotropy in the tree-level Symanzik action for several lattice spacings and target anisotropies. We use this method to tune the anisotropy parameters in full QCD, where we also exploit the diminishing effect of a well chosen smearing on the renormalization of the fermion anisotropy.

preprint2012arXiv

Building a frame and gauge free formulation of quantum mechanics

The wave function of quantum mechanics is not a boost invariant and gauge invariant quantity. Correspondingly, reference frame dependence and gauge dependence are inherited to most of the elements of the usual formulation of quantum mechanics (including operators, states and events). If a frame dependent and gauge dependent formalism is called, in short, a relative formalism, then the aim of the paper is to establish an absolute, i.e., frame and gauge free, reformulation of quantum mechanics. To fulfil this aim, we develop absolute quantities and the corresponding equations instead of the wave function and the Schrödinger equation. The absolute quantities have a more direct physical interpretation than the wave function has, and the corresponding equations express explicitly the independent physical aspects of the system which are contained in the Schrödinger equation in a mixed and more hidden form. Based on the absolute quantities and equations, events, states and physical quantities are introduced also in an absolute way. The formalism makes it possible to obtain some sharper versions of the uncertainty relation and to extend the validity of Ehrenfest's theorem. The absolute formulation allows wide extensions of quantum mechanics. To give examples, we discuss two known nonlinear extensions and, in close detail, a dissipative system. An argument is provided that the absolute formalism may lead to an explanation of the Aharonov-Bohm effect purely in terms of the electromagnetic field strength tensor. At last, on special relativistic and curved spacetimes, absolute quantities and equations instead of the Klein-Gordon wave function and equation are given, and their nonrelativistic limit is derived.

preprint2012arXiv

High-precision scale setting in lattice QCD

Scale setting is of central importance in lattice QCD. It is required to predict dimensional quantities in physical units. Moreover, it determines the relative lattice spacings of computations performed at different values of the bare coupling, and this is needed for extrapolating results into the continuum. Thus, we calculate a new quantity, $w_0$, for setting the scale in lattice QCD, which is based on the Wilson flow like the scale $t_0$ (M. Luscher, JHEP 1008 (2010) 071). It is cheap and straightforward to implement and compute. In particular, it does not involve the delicate fitting of correlation functions at asymptotic times. It typically can be determined on the few per-mil level. We compute its continuum extrapolated value in 2+1-flavor QCD for physical and non-physical pion and kaon masses, to allow for mass-independent scale setting even away from the physical mass point. We demonstrate its robustness by computing it with two very different actions (one of them with staggered, the other with Wilson fermions) and by showing that the results agree for physical quark masses in the continuum limit.

preprint2012arXiv

Magnetic susceptibility of QCD at zero and at finite temperature from the lattice

The response of the QCD vacuum to a constant external (electro)magnetic field is studied through the tensor polarization of the chiral condensate and the magnetic susceptibility at zero and at finite temperature. We determine these quantities using lattice configurations generated with the tree-level Symanzik improved gauge action and N_f=1+1+1 flavors of stout smeared staggered quarks with physical masses. We carry out the renormalization of the observables under study and perform the continuum limit both at T>0 and at T=0, using different lattice spacings. Finite size effects are studied by using various spatial lattice volumes. The magnetic susceptibilities χ_f reveal a diamagnetic behavior; we obtain at zero temperature χ_u=-(2.08 +/- 0.08) 1/GeV^2, χ_d=-(2.02 +/- 0.09) 1/GeV^2 and χ_s=-(3.4 +/- 1.4) 1/GeV^2 for the up, down and strange quarks, respectively, in the MSBar scheme at a renormalization scale of 2 GeV. We also find the polarization to change smoothly with the temperature in the confinement phase and then to drastically reduce around the transition region.

preprint2012arXiv

Precision SU(3) lattice thermodynamics for a large temperature range

We present the equation of state (pressure, trace anomaly, energy density and entropy density) of the SU(3) gauge theory from lattice field theory in an unprecedented precision and temperature range. We control both finite size and cut-off effects. The studied temperature window (0.7...1000 T_c) stretches from the glueball dominated system into the perturbative regime, which allows us to discuss the range of validity of these approaches. We also determine the preferred renormalization scale of the Hard Thermal Loop scheme and we fit the unknown g^6 order perturbative coefficient at extreme high temperatures T>100 T_c. We furthermore quantify the nonperturbative contribution to the trace anomaly using a simple functional form. Our high precision data allows one to have a complete theoretical description of the equation of state from T=0 all the way to the phase transition, through the transition region into the perturbative regime up to the Stefan-Boltzmann limit. We will discuss this description, too.

preprint2012arXiv

QCD equation of state at nonzero chemical potential: continuum results with physical quark masses at order mu^2

We determine the equation of state of QCD for nonzero chemical potentials via a Taylor expansion of the pressure. The results are obtained for N_f=2+1 flavors of quarks with physical masses, on various lattice spacings. We present results for the pressure, interaction measure, energy density, entropy density, and the speed of sound for small chemical potentials. At low temperatures we compare our results with the Hadron Resonance Gas model. We also express our observables along trajectories of constant entropy over particle number. A simple parameterization is given (the Matlab/Octave script parameterization.m, submitted to the arXiv along with the paper), which can be used to reconstruct the observables as functions of T and mu, or as functions of T and S/N.

preprint2012arXiv

QCD quark condensate in external magnetic fields

We present a comprehensive analysis of the light condensates in QCD with 1+1+1 sea quark flavors (with mass-degenerate light quarks of different electric charges) at zero and nonzero temperatures of up to 190 MeV and external magnetic fields B<1 GeV^2/e. We employ stout smeared staggered fermions with physical quark masses and extrapolate the results to the continuum limit. At low temperatures we confirm the magnetic catalysis scenario predicted by many model calculations, while around the crossover the condensate develops a complex dependence on the external magnetic field, resulting in a decrease of the transition temperature.

preprint2012arXiv

QCD thermodynamics with continuum extrapolated Wilson fermions I

QCD thermodynamics is considered using Wilson fermions in the fixed scale approach. The temperature dependence of the renormalized chiral condensate, quark number susceptibility and Polyakov loop is measured at four lattice spacings allowing for a controlled continuum limit. The light quark masses are fixed to heavier than physical values in this first study. Finite volume effects are ensured to be negligible by using approriately large box sizes. The final continuum results are compared with staggered fermion simulations performed in the fixed N_t approach. The same continuum renormalization conditions are used in both approaches and the final results agree perfectly.

preprint2012arXiv

The QCD phase diagram for external magnetic fields

The effect of an external (electro)magnetic field on the finite temperature transition of QCD is studied. We generate configurations at various values of the quantized magnetic flux with $N_f=2+1$ flavors of stout smeared staggered quarks, with physical masses. Thermodynamic observables including the chiral condensate and susceptibility, and the strange quark number susceptibility are measured as functions of the field strength. We perform the renormalization of the studied observables and extrapolate the results to the continuum limit using $N_t=6,8$ and 10 lattices. We also check for finite volume effects using various lattice volumes. We find from all of our observables that the transition temperature $T_c$ significantly decreases with increasing magnetic field. This is in conflict with various model calculations that predict an increasing $T_c(B)$. From a finite volume scaling analysis we find that the analytic crossover that is present at B=0 persists up to our largest magnetic fields $eB \approx 1 \textmd{GeV}^2$, and that the transition strength increases mildly up to this $eB\approx1 \textmd{GeV}^2$.

preprint2011arXiv

FK/Fpi from the Budapest-Marseille-Wuppertal Collaboration

Based on a series of lattice calculations we determine the ratio FK/Fpi in QCD. With experimental data from kaon decay and nuclear double beta decay, we obtain a precise determination of |Vus|. Our simulation includes 2+1 flavours of sea quarks, with three lattice spacings, large volumes and a simulated pion mass reaching down to about 190 MeV for a full control over the systematic uncertainties.

preprint2011arXiv

Lattice QCD at the physical point: light quark masses

Ordinary matter is described by six fundamental parameters: three couplings (gravitational, electromagnetic and strong) and three masses: the electron's (m_e) and those of the up (m_u) and down (m_d) quarks. An additional mass enters through quantum fluctuations: the strange quark mass (m_s). The three couplings and m_e are known with an accuracy of better than a few per mil. Despite their importance, $m_u$, $m_d$ (their average m_{ud}) and m_s are relatively poorly known: e.g. the Particle Data Group quotes them with conservative errors close to 25%. Here we determine these quantities with a precision below 2% by performing ab initio lattice quantum chromodynamics (QCD) calculations, in which all systematics are controlled. We use pion and quark masses down to (and even below) their physical values, lattice sizes of up to 6 fm, and five lattice spacings to extrapolate to continuum spacetime. All necessary renormalizations are performed nonperturbatively.

preprint2011arXiv

Lattice SU(3) thermodynamics and the onset of perturbative behaviour

We present the equation of state (pressure, trace anomaly, energy density and entropy density) of the SU(3) gauge theory from lattice field theory in an unprecedented precision and temperature range. We control both finite size and cut-off effects. The studied temperature window ($0.7... 1000 T_c$) stretches from the glueball dominated system into the perturbative regime, which allows us to discuss the range of validity of these approaches. From the critical couplings on fine lattices we get $T_c/\Lambdamsbar=1.26(7)$ and use this ratio to express the perturbative free energy in $T_c$ units. We also determine the preferred renormalization scale of the Hard Thermal Loop scheme and we fit the unknown $g^6$ order perturbative coefficient at extreme high temperatures $T>100T_c$. We furthermore quantify the nonperturbative contribution to the trace anomaly using two simple functional forms.

preprint2011arXiv

Precision computation of the kaon bag parameter

Indirect CP violation in K \rightarrow ππ decays plays a central role in constraining the flavor structure of the Standard Model (SM) and in the search for new physics. For many years the leading uncertainty in the SM prediction of this phenomenon was the one associated with the nonperturbative strong interaction dynamics in this process. Here we present a fully controlled lattice QCD calculation of these effects, which are described by the neutral kaon mixing parameter B_K . We use a two step HEX smeared clover-improved Wilson action, with four lattice spacings from a\approx0.054 fm to a\approx0.093 fm and pion masses at and even below the physical value. Nonperturbative renormalization is performed in the RI-MOM scheme, where we find that operator mixing induced by chiral symmetry breaking is very small. Using fully nonperturbative continuum running, we obtain our main result B_K^{RI}(3.5GeV)=0.531(6)_{stat}(2)_{sys}. A perturbative 2-loop conversion yields B_K^{MSbar-NDR}(2GeV)=0.564(6)_{stat}(3)_{sys}(6)_{PT}, which is in good agreement with current results from fits to experimental data.

preprint2011arXiv

Spectra of heavy-light and heavy-heavy mesons containing charm quarks, including higher spin states for $N_f=2+ 1$

We study the spectra of heavy-light and heavy-heavy mesons containing charm quarks, including higher spin states. We use two sets of $N_f = 2 + 1$ gauge configurations, one set from QCDSF using the SLiNC action, and the other configurations from the Budapest-Marseille-Wuppertal collaboration, using the HEX smeared clover action. To extract information about the excited states, we choose a suitable basis of operators to implement the variational method.

preprint2011arXiv

The finite temperature QCD transition in external magnetic fields

The effect of an external magnetic field on the finite temperature transition of QCD is studied. We measure thermodynamic observables including the quark condensates and susceptibilities and the strange quark number susceptibility. We generate configurations at various values of the quantized magnetic flux with N_f=2+1 flavors of stout smeared staggered quarks at physical quark masses. We perform the renormalization of our observables and approach the continuum limit with N_t=6,8 and 10 lattices. We also check finite volume effects using various lattice volumes. Our main result is that the transition temperature significantly decreases with growing magnetic field, and that the transition remains an analytic crossover up to our largest external field eB ~ 1 GeV^2.

preprint2011arXiv

The QCD phase diagram at nonzero quark density

We determine the phase diagram of QCD on the μ-T plane for small to moderate chemical potentials. Two transition lines are defined with two quantities, the chiral condensate and the strange quark number susceptibility. The calculations are carried out on N_t =6,8 and 10 lattices generated with a Symanzik improved gauge and stout-link improved 2+1 flavor staggered fermion action using physical quark masses. After carrying out the continuum extrapolation we find that both quantities result in a similar curvature of the transition line. Furthermore, our results indicate that in leading order the width of the transition region remains essentially the same as the chemical potential is increased.

preprint2010arXiv

Lattice QCD at the physical point: Simulation and analysis details

We give details of our precise determination of the light quark masses m_{ud}=(m_u+m_d)/2 and m_s in 2+1 flavor QCD, with simulated pion masses down to 120 MeV, at five lattice spacings, and in large volumes. The details concern the action and algorithm employed, the HMC force with HEX smeared clover fermions, the choice of the scale setting procedure and of the input masses. After an overview of the simulation parameters, extensive checks of algorithmic stability, autocorrelation and (practical) ergodicity are reported. To corroborate the good scaling properties of our action, explicit tests of the scaling of hadron masses in N_f=3 QCD are carried out. Details of how we control finite volume effects through dedicated finite volume scaling runs are reported. To check consistency with SU(2) Chiral Perturbation Theory the behavior of M_π^2/m_{ud} and F_πas a function of m_{ud} is investigated. Details of how we use the RI/MOM procedure with a separate continuum limit of the running of the scalar density R_S(μ,μ') are given. This procedure is shown to reproduce the known value of r_0m_s in quenched QCD. Input from dispersion theory is used to split our value of m_{ud} into separate values of m_u and m_d. Finally, our procedure to quantify both systematic and statistical uncertainties is discussed.

preprint2010arXiv

Sigma term and strangeness content of the nucleon

A status report is given for a joint project of the Budapest-Marseille-Wuppertal collaboration and the Regensburg group to study the quark mass-dependence of octet baryons in SU(3) Baryon XPT. This formulation is expected to extend to larger masses than Heavy-Baryon XPT. Its applicability is tested with 2+1 flavor data which cover three lattice spacings and pion masses down to about 190 MeV, in large volumes. Also polynomial and rational interpolations in M_π^2 and M_K^2 are used to assess the uncertainty due to the ansatz. Both frameworks are combined to explore the precision to be expected in a controlled determination of the nucleon sigma term and strangeness content.

preprint2010arXiv

The ratio FK/Fpi in QCD

We determine the ratio FK/Fpi in QCD with Nf=2+1 flavors of sea quarks, based on a series of lattice calculations with three different lattice spacings, large volumes and a simulated pion mass reaching down to about 190 MeV. We obtain FK/Fpi=1.192 +/- 0.007(stat) +/- 0.006(syst). This result is then used to give an updated value of the CKM matrix element |Vus|. The unitarity relation for the first row of this matrix is found to be well observed.

preprint2009arXiv

The curvature of the QCD phase transition line

We determine the curvature of the phase transition line in the mu-T plane through an analysis of various observables, including the Polyakov loop, the quark number susceptibilities and the susceptibility of the chiral condensate. The second derivative of these quantities with respect to mu was calculated. The measurements were carried out on N_T = 4,6,8 and 10 lattices generated with a Symanzik improved gauge and stout-link improved 2+1 flavour staggered fermion action using physical quark masses.

preprint2004arXiv

Pentaquark hadrons from lattice QCD

We study spin 1/2 isoscalar and isovector candidates in both parity channels for the recently discovered Θ^+(1540) pentaquark particle in quenched lattice QCD. Our analysis takes into account all possible uncertainties, such as statistical, finite size and quenching errors when performing the chiral and continuum extrapolations and we have indications that our signal is separated from scattering states. The lowest mass that we find in the I^P=0^- channel is in complete agreement with the experimental value of the Θ^+ mass. On the other hand, the lowest mass state in the opposite parity I^P=0^+ channel is much higher. Our findings suggests that the parity of the Θ^+ is negative.

preprint2004arXiv

Strong Neutrino-Nucleon Interactions at Ultrahigh Energies as a Solution to the GZK Puzzle

After a short review of the ultrahigh energy cosmic ray puzzle - the apparent observation of cosmic rays originating from cosmological distances with energies above the expected Greisen-Zatsepin-Kuzmin cutoff 4x10^{19} eV - we consider strongly interacting neutrino scenarios as an especially interesting solution. We show that all features of the ultrahigh energy cosmic ray spectrum from 10^{17} eV to 10^{21} eV can be described to originate from a simple power-like injection spectrum of protons, under the assumption that the neutrino-nucleon cross-section is significantly enhanced at center of mass energies above ~ 100 TeV. In such a scenario, the cosmogenic neutrinos produced during the propagation of protons through the cosmic microwave background initiate air showers in the atmosphere, just as the protons. The total air shower spectrum induced by protons and neutrinos shows excellent agreement with the observations. We shortly discuss TeV-scale extensions of the Standard Model which may lead to a realization of a strongly interacting neutrino scenario. We emphasize, however, that such a scenario may even be realized within the standard electroweak model: electroweak instanton/sphaleron induced processes may get strong at ultrahigh energies. Possible tests of strongly interacting neutrino scenarios range from observations at cosmic ray facilities and neutrino telescopes to searches at lepton nucleon scattering experiments.