Source author record

Z. Fodor

Z. Fodor 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

96works
11topics
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

96 published item(s)

preprint2024arXiv

Search for a critical point of strongly-interacting matter in central $^{40}$Ar +$^{45}$Sc collisions at 13$A$-75$A$ GeV/$c$ beam momentum

The critical point of strongly interacting matter is searched for at the CERN SPS by the NA61/SHINE experiment in central $^{40}$Ar +$^{45}$Sc collisions at 13$A$, 19$A$, 30$A$, 40$A$, and 75$A$ GeV/$c$. The dependence of the second-order scaled factorial moments of proton multiplicity distributions on the number of subdivisions in transverse momentum space is measured. The intermittency analysis uses statistically independent data sets for every subdivision in transverse and cumulative-transverse momentum variables. The results obtained do not indicate the searched intermittent pattern. An upper limit on the fraction of correlated protons and the intermittency index is obtained based on a comparison with the Power-law Model.

preprint2024arXiv

Thermal QCD phase transition with dynamical chiral fermions

We discuss properties of Quantum Chromodynamics at finite temperature obtained by means of lattice simulations with overlap fermions. This fermion discretization preserves chiral symmetry even at finite lattice spacing. We present details of the lattice formulation, first results for the chiral observables and discuss the behaviour of the system near the chiral thermal phase transition.

preprint2022arXiv

$K^{*}(892)^0$ meson production in inelastic $p+p$ interactions at 40 and 80 GeV/$c$ beam momenta measured by NA61/SHINE at the CERN SPS

Measurements of $K^{*}(892)^0$ resonance production via its $K^{+}π^{-}$ decay mode in inelastic $p+p$ collisions at beam momenta 40 and 80 GeV/$c$ ($\sqrt{s_{NN}}=8.8$ and 12.3 GeV) are presented. The data were recorded by the NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. The \textit{template} method was used to extract the $K^{*}(892)^0$ signal. Transverse momentum and rapidity spectra were obtained. The mean multiplicities of $K^{*}(892)^0$ mesons were found to be $(35.1 \pm 1.3 \mathrm{(stat)} \pm 3.6 \mathrm{(sys)) \cdot 10^{-3}}$ at 40 GeV/$c$ and $(58.3 \pm 1.9 \mathrm{(stat)} \pm 4.9 \mathrm{(sys)) \cdot 10^{-3}}$ at 80 GeV/$c$. The NA61/SHINE results are compared with the EPOS1.99 and Hadron Resonance Gas models as well as with world data. The transverse mass spectra of $K^{*}(892)^0$ mesons and other particles previously reported by NA61/SHINE were fitted within the Blast-Wave model. The transverse flow velocities are close to 0.1--0.2 of the speed of light and are significantly smaller than the ones determined in heavy nucleus-nucleus interactions at the same beam momenta.

preprint2022arXiv

$K^{0}_{S}$ meson production in inelastic $\textit{p+p}$ interactions at 158 GeV/c beam momentum measured by NA61/SHINE at the CERN SPS

The production of $K^{0}_{S}$ mesons in inelastic $\textit{p+p}$ collisions at beam momentum 158 GeV/c ($\sqrt{s_{NN}}=17.3$ GeV) was measured with the NA61/SHINE spectrometer at the CERN Super Proton Synchrotron. Double-differential distributions were obtained in transverse momentum and rapidity. The mean multiplicity of $K^{0}_{S}$ was determined to be $0.162 \pm 0.001 (stat.) \pm 0.011 (sys.)$. The results on $K^{0}_{S}$ production are compared with model predictions (EPOS 1.99, SMASH 2.0, PHSD and UrQMD 3.4 models) as well as with published world data.

preprint2022arXiv

Measurements of $Ξ^{-}$ and $\overlineΞ^{+}$ production in proton-proton interactions at $\sqrt{s_{NN}}$ = 17.3 GeV in the NA61/SHINE experiment

The production of $Ξ(1321)^{-}$ and $\overlineΞ(1321)^{+}$ hyperons in inelastic p+p interactions is studied in a fixed target experiment at a beam momentum of 158 GeV/textit{c}. Double differential distributions in rapidity y and transverse momentum $p_{T}$ are obtained from a sample of 33M inelastic events. They allow to extrapolate the spectra to full phase space and to determine the mean multiplicity of both $Ξ^{-}$ and $\overlineΞ^{+}$. The rapidity and transverse momentum spectra are compared to transport model predictions. The $Ξ^{-}$ mean multiplicity in inelastic p+p interactions at 158~\GeVc is used to quantify the strangeness enhancement in A+A collisions at the same centre-of-mass energy per nucleon pair.

preprint2022arXiv

Precision study of the continuum SU(3) Yang-Mills theory: how to use parallel tempering to improve on supercritical slowing down for first order phase transitions

We perform large scale simulations to characterize the transition in quenched QCD. It is shown by a rigorous finite size scaling that the transition is of first order. After this qualitative feature quantitative results are obtained with unprecedented precision: we calculate the transition temperature $w_0T_c$=0.25384(23), -- which is the first per-mill accurate result in QCD thermodynamics -- and the latent heat $ΔE/T_c^4$=1.025(21)(27) in both cases carrying out controlled continuum and infinite volume extrapolations. As it is well known the cost of lattice simulations explodes in the vicinity of phase transitions, a phenomenon called critical slowing down for second order phase transitions and supercritical slowing down for first order phase transitions. We show that a generalization of the parallel tempering algorithm of Marinari and Parisi [Europhys. Lett. 19, 451 (1992)] originally for spin systems can efficiently overcome these difficulties even if the transition is of first order, like in the case of QCD without quarks, or with very heavy quarks. We also report on our investigations on the autocorrelation times and other details.

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.

preprint2021arXiv

Measurement of the production cross section of 31 GeV/$c$ protons on carbon via beam attenuation in a 90-cm-long target

The production cross section of 30.92 GeV/$c$ protons on carbon is measured by the NA61/SHINE spectrometer at the CERN SPS by means of beam attenuation in a copy (replica) of the 90-cm-long target of the T2K neutrino oscillation experiment. The employed method for direct production cross-section estimation minimizes model corrections for elastic and quasi-elastic interactions. The obtained production cross section is $σ_\mathrm{prod}~=~227.6~\pm~0.8\mathrm{(stat)}~_{-~3.2}^{+~1.9}\mathrm{(sys)}~{-~0.8}\mathrm{(mod)}$ mb. It is in agreement with previous NA61/SHINE results obtained with a thin carbon target, while providing improved precision with a total fractional uncertainty of less than 2$\%$. This direct measurement is performed to reduce the uncertainty on the T2K neutrino flux prediction associated with the re-weighting of the interaction rate of neutrino-yielding hadrons.

preprint2021arXiv

Measurements of $π^-$ production in $^7$Be+$^9$Be collisions at beam momenta from 19$A$ to 150$A$GeV/$c$ in the NA61/SHINE experiment at the CERN SPS

The NA61/SHINE collaboration studies at the CERN Super Proton Synchrotron (SPS) the onset of deconfinement in hadronic matter by the measurement of particle production in collisions of nuclei with various sizes at a set of energies covering the SPS energy range. This paper presents results on inclusive double-differential spectra and mean multiplicities of $π^{-}$ mesons produced in the 5\% most \textit{central} $^7$Be+$^9$Be collisions at beam momenta of 19$A$, 30$A$, 40$A$, 75$A$ and 150$A$ GeV/$c$ obtained by the so-called $h^-$ method which does not require any particle identification. The shape of the transverse mass spectra differs from the shapes measured in central Pb+Pb collisions and inelastic p+p interactions. The normalized width of the rapidity distribution decreases with increasing collision energy and is in between the results for inelastic nucleon-nucleon and central Pb+Pb collisions. The mean multiplicity of pions per wounded nucleon in \textit{central} $^7$Be+$^9$Be collisions is close to that in central Pb+Pb collisions up to 75$A$GeV/$c$. However, at the top SPS energy the result is close to the one for inelastic nucleon-nucleon interactions. The results are discussed in the context of predictions for the onset of deconfinement at the CERN SPS collision energies.

preprint2021arXiv

Spectra and mean multiplicities of $π^{-}$ in $central$ ${}^{40}$Ar+${}^{45}$Sc collisions at 13$A$, 19$A$, 30$A$, 40$A$, 75$A$ and 150$A$ GeV/$c$ beam momenta measured by the NA61/SHINE spectrometer at the CERN SPS

The physics goal of the strong interaction program of the NA61/SHINE experiment at the CERN Super Proton Synchrotron (SPS) is to study the phase diagram of hadronic matter by a scan of particle production in collisions of nuclei with various sizes at a set of energies covering the SPS energy range. This paper presents differential inclusive spectra of transverse momentum, transverse mass and rapidity of $π^{-}$ mesons produced in $central$ ${}^{40}$Ar+${}^{45}$Sc collisions at beam momenta of 13$A$, 19$A$, 30$A$, 40$A$, 75$A$ and 150$A$ GeV/$c$. Energy and system size dependence of parameters of these distributions -- mean transverse mass, the inverse slope parameter of transverse mass spectra, width of the rapidity distribution and mean multiplicity -- are presented and discussed. Furthermore, the dependence of the ratio of the mean number of produced pions to the mean number of wounded nucleons on the collision energy was derived. The results are compared to predictions of several models.

preprint2021arXiv

Two-particle correlations in azimuthal angle and pseudorapidity in central $^7$Be+$^9$Be collisions at the CERN Super Proton Synchrotron

A measurement of charged hadron pair correlations in two-dimensional $ΔηΔϕ$ space is presented. The analysis is based on total 30 million central Be+Be collisions observed in the NA61/SHINE detector at the CERN SPS for incident beam momenta of 19$A$, 30$A$, 40$A$, 75$A$, and 150$A$ GeV/$c$. Measurements were carried out for unlike-sign and like-sign charge hadron pairs independently. The $C(Δη,Δϕ)$ correlation functions were compared with results from a similar analysis on p+p interactions at similar beam momenta per nucleon. General trends of the back-to-back correlations are similar in central Be+Be collisions and p+p interactions, but are suppressed in magnitude due to the increased combinatorial background. Predictions from the EPOS and UrQMD models are compared to the measurements. Evolution of an enhancement around $(Δη,Δϕ) = (0,0)$ with incident energy is observed in central Be+Be collisions. It is not predicted by both models and almost non-existing in proton-proton collisions at the same momentum per nucleon.

preprint2020arXiv

$K^{*}(892)^0$ meson production in inelastic p+p interactions at 158 GeV/$c$ beam momentum measured by NA61/SHINE at the CERN SPS

The measurement of $K^{*}(892)^0$ resonance production via its $K^{+}π^{-}$ decay mode in inelastic p+p collisions at beam momentum 158 GeV/$c$ ($\sqrt{s_{NN}}=17.3$ GeV) is presented. The data were recorded by the NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. The \textit{template} method was used to extract the $K^{*}(892)^0$ signal and double differential transverse momentum and rapidity spectra were obtained. The full phase-space mean multiplicity of $K^{*}(892)^0$ mesons was found to be $(78.44 \pm 0.38 \mathrm{(stat)} \pm 6.0 \mathrm{(sys)) \cdot 10^{-3}}$. The NA61/SHINE results are compared with the EPOS1.99 and Hadron Resonance Gas models as well as with world data from p+p and nucleus-nucleus collisions.

preprint2020arXiv

Ab-initio calculation of the proton and the neutron's scalar couplings for new physics searches

Many low-energy, particle-physics experiments seek to reveal new fundamental physics by searching for very rare scattering events on atomic nuclei. The interpretation of their results requires quantifying the non-linear effects of the strong interaction on the spin-independent couplings of this new physics to protons and neutrons. Here we present a fully-controlled, ab-initio calculation of these couplings to the quarks within those constituents of nuclei. We use lattice quantum chromodynamics computations for the four lightest species of quarks and heavy-quark expansions for the remaining two. We determine each of the six quark contributions with an accuracy better than 15%. Our results are especially important for guiding and interpreting experimental searches for our universe's dark matter.

preprint2020arXiv

Measurements of hadron production in $π^{+}$ + C and $π^{+}$ + Be interactions at 60 GeV/$c$

Precise knowledge of hadron production rates in the generation of neutrino beams is necessary for accelerator-based neutrino experiments to achieve their physics goals. NA61/SHINE, a large-acceptance hadron spectrometer, has recorded hadron+nucleus interactions relevant to ongoing and future long-baseline neutrino experiments at Fermi National Accelerator Laboratory. This paper presents three analyses of interactions of 60 GeV/$c$ $π^+$ with thin, fixed carbon and beryllium targets. Integrated production and inelastic cross sections were measured for both of these reactions. In an analysis of strange, neutral hadron production, differential production multiplicities of $K^0_{S}$, $Λ$ and anti-$Λ$ were measured. Lastly, in an analysis of charged hadron production, differential production multiplicities of $π^+$, $π^-$, $K^+$, $K^-$ and protons were measured. These measurements will enable long-baseline neutrino experiments to better constrain predictions of their neutrino flux in order to achieve better precision on their neutrino cross section and oscillation measurements.

preprint2020arXiv

Proton-proton interactions and onset of deconfinement

The NA61/SHINE experiment at the CERN SPS is performing a uniqe study of the phase diagram of strongly interacting matter by varying collision energy and nuclear mass number of colliding nuclei. In central Pb+Pb collisions the NA49 experiment found structures in the energy dependence of several observables in the CERN SPS energy range that had been predicted for the transition to a deconfined phase. New measurements of NA61/SHINE find intriguing similarities in p+p interactions for which no deconfinement transition is expected at SPS energies. Possible implications will be discussed.

preprint2020arXiv

Search for an exotic S=-2, Q=-2 baryon resonance in proton-proton interactions at $\sqrt{s_{NN}}$ = 17.3 GeV

Pentaquark states have been extensively investigated theoretically in the context of the constituent quark model. In this paper results of an experimental search for pentaquarks in the \Xim\pim, \Xim\pip, \Xip\pim and \Xip\pip invariant mass spectra in proton-proton interactions at $\sqrt{s}$=17.3~\GeV are presented. Previous possible evidence from the NA49 collaboration of the existence of a narrow \Xim\pim baryon resonance in p+p interactions is not confirmed with almost 10 times greater event statistics. The search was performed using the \NASixtyOne detector which reuses the main components of the NA49 apparatus. No signal was observed with either the selection cuts of NA49 or newly optimised cuts.

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

Centrality dependence of subthreshold $ϕ$ meson production in Ni+Ni collisions at 1.9A GeV

We analysed the $ϕ$ meson production in central Ni+Ni collisions at the beam kinetic energy of 1.93A GeV with the FOPI spectrometer and found the production probability per event of $[8.6 ~\pm~ 1.6 ~(\text{stat}) \pm 1.5 ~(\text{syst})] \times 10^{-4}$. This new data point allows for the first time to inspect the centrality dependence of the subthreshold $ϕ$ meson production in heavy-ion collisions. The rise of $ϕ$ meson multiplicity per event with mean number of participants can be parameterized by the power function with exponent $α= 1.8 \pm 0.6$. The ratio of $ϕ$ to $\text{K}^-$ production yields seems not to depend within the experimental uncertainties on the collision centrality, and the average of measured values was found to be $0.36 \pm 0.05$.

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

Measurements of $π^{\pm}$ differential yields from the surface of the T2K replica target for incoming 31 GeV/c protons with the NA61/SHINE spectrometer at the CERN SPS

Measurements of particle emission from a replica of the T2K 90 cm-long carbon target were performed in the NA61/SHINE experiment at CERN SPS, using data collected during a high-statistics run in 2009. An efficient use of the long-target measurements for neutrino flux predictions in T2K requires dedicated reconstruction and analysis techniques. Fully-corrected differential yields of $π^\pm$-mesons from the surface of the T2K replica target for incoming 31 GeV/c protons are presented. A possible strategy to implement these results into the T2K neutrino beam predictions is discussed and the propagation of the uncertainties of these results to the final neutrino flux is performed.

preprint2016arXiv

Measurements of $π^\pm$, $K^\pm$, $K^0_S$, $Λ$ and proton production in proton-carbon interactions at 31 GeV/$c$ with the NA61/SHINE spectrometer at the CERN SPS

Measurements of hadron production in p+C interactions at 31 GeV/c are performed using the NA61/ SHINE spectrometer at the CERN SPS. The analysis is based on the full set of data collected in 2009 using a graphite target with a thickness of 4% of a nuclear interaction length. Inelastic and production cross sections as well as spectra of $π^\pm$, $K^\pm$, p, $K^0_S$ and $Λ$ are measured with high precision. These measurements are essential for improved calculations of the initial neutrino fluxes in the T2K long-baseline neutrino oscillation experiment in Japan. A comparison of the NA61/SHINE measurements with predictions of several hadroproduction models is presented.

preprint2016arXiv

Multiplicity and transverse momentum fluctuations in inelastic proton-proton interactions at the CERN Super Proton Synchrotron

Measurements of multiplicity and transverse momentum fluctuations of charged particles were performed in inelastic p+p interactions at 20, 31, 40, 80 and 158 GeV/c beam momentum. Results for the scaled variance of the multiplicity distribution and for three strongly intensive measures of multiplicity and transverse momentum fluctuations \$Δ[P_{T},N]\$, \$Σ[P_{T},N]\$ and \$Φ_{p_T}\$ are presented. For the first time the results on fluctuations are fully corrected for experimental biases. The results on multiplicity and transverse momentum fluctuations significantly deviate from expectations for the independent particle production. They also depend on charges of selected hadrons. The string-resonance Monte Carlo models EPOS and UrQMD do not describe the data. The scaled variance of multiplicity fluctuations is significantly higher in inelastic p+p interactions than in central Pb+Pb collisions measured by NA49 at the same energy per nucleon. This is in qualitative disagreement with the predictions of the Wounded Nucleon Model. Within the statistical framework the enhanced multiplicity fluctuations in inelastic p+p interactions can be interpreted as due to event-by-event fluctuations of the fireball energy and/or volume.

preprint2016arXiv

Production of $Λ$ hyperons in inelastic p+p interactions at 158 GeV/$c$

Inclusive production of $Λ$-hyperons was measured with the large acceptance NA61/SHINE spectrometer at the CERN SPS in inelastic p+p interactions at beam momentum of 158~\GeVc. Spectra of transverse momentum and transverse mass as well as distributions of rapidity and x$_{_F}$ are presented. The mean multiplicity was estimated to be $0.120\,\pm0.006\;(stat.)\,\pm 0.010\;(sys.)$. The results are compared with previous measurements and predictions of the EPOS, UrQMD and FRITIOF models.

preprint2016arXiv

Production of deuterium, tritium, and $^3$He in central Pb+Pb collisions at 20A, 30A, 40A, 80A, and 158A GeV at the CERN SPS

Production of $d$, $t$, and $^3$He nuclei in central Pb+Pb interactions was studied at five collision energies ($\sqrt{s_{NN}}=$ 6.3, 7.6, 8.8, 12.3, and 17.3 GeV) with the NA49 detector at the CERN SPS. Transverse momentum spectra, rapidity distributions, and particle ratios were measured. Yields are compared to predictions of statistical models. Phase-space distributions of light nuclei are discussed and compared to those of protons in the context of a coalescence approach. The coalescence parameters $B_2$ and $B_3$, as well as coalescence radii for $d$ and $^3$He were determined as a function of transverse mass at all energies.

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

Strange meson production in Al+Al collisions at 1.9A GeV

The production of K$^+$, K$^-$ and $φ$(1020) mesons is studied in Al+Al collisions at a beam energy of 1.9A GeV which is close or below the production threshold in NN reactions. Inverse slopes, anisotropy parameters, and total emission yields of K$^{\pm}$ mesons are obtained. A comparison of the ratio of kinetic energy distributions of K$^-$ and K$^+$ mesons to the HSD transport model calculations suggests that the inclusion of the in-medium modifications of kaon properties is necessary to reproduce the ratio. The inverse slope and total yield of $ϕ$ mesons are deduced. The contribution to K$^-$ production from $ϕ$ meson decays is found to be [17 $\pm$ 3 (stat) $^{+2}_{-7}$ (syst)] %. The results are in line with previous K$^{\pm}$ and $ϕ$ data obtained for different colliding systems at similar incident beam energies.

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.

preprint2016arXiv

Up and down quark masses and corrections to Dashen's theorem from lattice QCD and quenched QED

In a previous letter (arXiv:1306.2287) we determined the isospin mass splittings of the baryon octet from a lattice calculation based on quenched QED and $N_f{=}2{+}1$ QCD simulations with 5 lattice spacings down to $0.054~\mathrm{fm}$, lattice sizes up to $6~\mathrm{fm}$ and average up-down quark masses all the way down to their physical value. Using the same data we determine here the corrections to Dashen's theorem and the individual up and down quark masses. For the parameter which quantifies violations to Dashens's theorem, we obtain $ε=0.73(2)(5)(17)$, where the first error is statistical, the second is systematic, and the third is an estimate of the QED quenching error. For the light quark masses we obtain, $m_u=2.27(6)(5)(4)~\mathrm{MeV}$ and $m_d=4.67(6)(5)(4)~\mathrm{MeV}$ in the $\bar{\mathrm{MS}}$ scheme at $2~\mathrm{GeV}$ and the isospin breaking ratios $m_u/m_d=0.485(11)(8)(14)$, $R=38.2(1.1)(0.8)(1.4)$ and $Q=23.4(0.4)(0.3)(0.4)$. Our results exclude the $m_u=0$ solution to the strong CP problem by more than $24$ standard deviations.

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

Critical fluctuations of the proton density in A+A collisions at $158A$ GeV

We look for fluctuations expected for the QCD critical point using an intermittency analysis in the transverse momentum phase space of protons produced around midrapidity in the 12.5\% most central C+C, Si+Si and Pb+Pb collisions at the maximum SPS energy of 158$A$~GeV. We find evidence of power-law fluctuations for the Si+Si data. The fitted power-law exponent $ϕ_{2} = 0.96^{+0.38}_{-0.25}\text{ (stat.)}$ $\pm 0.16\text{ (syst.)}$ is consistent with the va\-lue expected for critical fluctuations. Power-law fluctuations had previously also been observed in low-mass $π^+ π^-$ pairs in the same Si+Si collisions.

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

Measurement of event-by-event transverse momentum and multiplicity fluctuations using strongly intensive measures $Δ[P_T, N]$ and $Σ[P_T, N]$ in nucleus-nucleus collisions at the CERN Super Proton Synchrotron

Results from the NA49 experiment at the CERN SPS are presented on event-by-event transverse momentum and multiplicity fluctuations of charged particles, produced at forward rapidities in central Pb+Pb interactions at beam momenta 20$A$, 30$A$, 40$A$, 80$A$, and 158$A$ GeV/c, as well as in systems of different size ($p+p$, C+C, Si+Si, and Pb+Pb) at 158$A$ GeV/c. This publication extends the previous NA49 measurements of the strongly intensive measure $Φ_{p_T}$ by a study of the recently proposed strongly intensive measures of fluctuations $Δ[P_T, N]$ and $Σ[P_T, N]$. In the explored kinematic region transverse momentum and multiplicity fluctuations show no significant energy dependence in the SPS energy range. However, a remarkable system size dependence is observed for both $Δ[P_T, N]$ and $Σ[P_T, N]$, with the largest values measured in peripheral Pb+Pb interactions. The results are compared with NA61/SHINE measurements in $p+p$ collisions, as well as with predictions of the UrQMD and EPOS models.

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

Azimuthal Emission Patterns of $K^{+}$ and of $ K^{-} $ Mesons in Ni + Ni Collisions near the Strangeness Production Threshold

Azimuthal emission patterns of $K^\pm$ mesons have been measured in Ni + Ni collisions with the FOPI spectrometer at a beam kinetic energy of 1.91 A GeV. The transverse momentum $p_{T}$ integrated directed and elliptic flow of $K^{+}$ and $K^{-}$ mesons as well as the centrality dependence of $p_{T}$ - differential directed flow of $K^{+}$ mesons are compared to the predictions of HSD and IQMD transport models. The data exhibits different propagation patterns of $K^{+}$ and $K^{-}$ mesons in the compressed and heated nuclear medium and favor the existence of a kaon-nucleon in-medium potential, repulsive for $K^{+}$ mesons and attractive for $K^{-}$ mesons.

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

Influence of $ϕ$ mesons on negative kaons in Ni+Ni collisions at 1.91A GeV beam energy

$ϕ$ and K$^-$ mesons from Ni+Ni collisions at the beam energy of 1.91A GeV have been measured by the FOPI spectrometer, with a trigger selecting central and semi-central events amounting to 51% of the total cross section. The phase space distributions, and the total yield of K$^-$, as well as the kinetic energy distribution and the total yield of $ϕ$ mesons are presented. The $ϕ$\K$^-$ ratio is found to be $0.44 \pm 0.07(\text{stat}) ^{+0.18}_{-0.12} (\text{syst})$, meaning that about 22% of K$^-$ mesons originate from the decays of $ϕ$ mesons, occurring mostly in vacuum. The inverse slopes of direct kaons are up to about 15 MeV larger than the ones extracted within the one-source model, signalling that a considerable share of gap between the slopes of K$^+$ and K$^-$ could be explained by the contribution of $ϕ$ mesons to negative kaons.

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

Measurement of negatively charged pion spectra in inelastic p+p interactions at $p_{lab}$ = 20, 31, 40, 80 and 158 GeV/c

We present experimental results on inclusive spectra and mean multiplicities of negatively charged pions produced in inelastic p+p interactions at incident projectile momenta of 20, 31, 40, 80 and 158 GeV/c ($\sqrt{s} = $ 6.3, 7.7, 8.8, 12.3 and 17.3 GeV, respectively). The measurements were performed using the large acceptance NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. Two-dimensional spectra are determined in terms of rapidity and transverse momentum. Their properties such as the width of rapidity distributions and the inverse slope parameter of transverse mass spectra are extracted and their collision energy dependences are presented. The results on inelastic p+p interactions are compared with the corresponding data on central Pb+Pb collisions measured by the NA49 experiment at the CERN SPS. The results presented in this paper are part of the NA61/SHINE ion program devoted to the study of the properties of the onset of deconfinement and search for the critical point of strongly interacting matter. They are required for interpretation of results on nucleus-nucleus and proton-nucleus collisions.

preprint2014arXiv

NA61/SHINE facility at the CERN SPS: beams and detector system

NA61/SHINE (SPS Heavy Ion and Neutrino Experiment) is a multi-purpose experimental facility to study hadron production in hadron-proton, hadron-nucleus and nucleus-nucleus collisions at the CERN Super Proton Synchrotron. It recorded the first physics data with hadron beams in 2009 and with ion beams (secondary 7Be beams) in 2011. NA61/SHINE has greatly profited from the long development of the CERN proton and ion sources and the accelerator chain as well as the H2 beamline of the CERN North Area. The latter has recently been modified to also serve as a fragment separator as needed to produce the Be beams for NA61/SHINE. Numerous components of the NA61/SHINE set-up were inherited from its predecessors, in particular, the last one, the NA49 experiment. Important new detectors and upgrades of the legacy equipment were introduced by the NA61/SHINE Collaboration. This paper describes the state of the NA61/SHINE facility - the beams and the detector system - before the CERN Long Shutdown I, which started in March 2013.

preprint2014arXiv

QCD and strongly coupled gauge theories: challenges and perspectives

We highlight the progress, current status, and open challenges of QCD-driven physics, in theory and in experiment. We discuss how the strong interaction is intimately connected to a broad sweep of physical problems, in settings ranging from astrophysics and cosmology to strongly-coupled, complex systems in particle and condensed-matter physics, as well as to searches for physics beyond the Standard Model. We also discuss how success in describing the strong interaction impacts other fields, and, in turn, how such subjects can impact studies of the strong interaction. In the course of the work we offer a perspective on the many research streams which flow into and out of QCD, as well as a vision for future developments.

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

Inclusive production of protons, anti-protons, neutrons, deuterons and tritons in p+C collisions at 158 GeV/c beam momentum

The production of protons, anti-protons, neutrons, deuterons and tritons in minimum bias p+C interactions is studied using a sample of 385 734 inelastic events obtained with the NA49 detector at the CERN SPS at 158 GeV/c beam momentum. The data cover a phase space area ranging from 0 to 1.9 GeV/c in transverse momentum and in Feynman x from -0.80 to 0.95 for protons, from -0.2 to 0.4 for anti-protons and from 0.2 to 0.95 for neutrons. Existing data in the far backward hemisphere are used to extend the coverage for protons and light nuclear fragments into the region of intranuclear cascading. The use of corresponding data sets obtained in hadron-proton collisions with the same detector allows for the detailed analysis and model-independent separation of the three principle components of hadronization in p+C interactions, namely projectile fragmentation, target fragmentation of participant nucleons and intranuclear cascading.

preprint2013arXiv

Measurements of Production Properties of K0S mesons and Lambda hyperons in Proton-Carbon Interactions at 31 GeV/c

Spectra of K0S mesons and Lambda hyperons were measured in p+C interactions at 31 GeV/c with the large acceptance NA61/SHINE spectrometer at the CERN SPS. The data were collected with an isotropic graphite target with a thickness of 4% of a nuclear interaction length. Interaction cross sections, charged pion spectra, and charged kaon spectra were previously measured using the same data set. Results on K0S and Lambda production in p+C interactions serve as reference for the understanding of the enhancement of strangeness production in nucleus-nucleus collisions. Moreover, they provide important input for the improvement of neutrino flux predictions for the T2K long baseline neutrino oscillation experiment in Japan. Inclusive production cross sections for K0S and Lambda are presented as a function of laboratory momentum in intervals of the laboratory polar angle covering the range from 0 up to 240 mrad. The results are compared with predictions of several hadron production models. The K0S mean multiplicity in production processes <n_K0S> and the inclusive cross section for K0S production were measured and amount to 0.127 +- 0.005 (stat) +- 0.022 (sys) and 29.0 +- 1.6 (stat) +- 5.0 (sys) mb, respectively.

preprint2013arXiv

Phase-space dependence of particle-ratio fluctuations in Pb+Pb collisions from 20A to 158A GeV beam energy

A novel approach, the identity method, was used for particle identification and the study of fluctuations of particle yield ratios in Pb+Pb collisions at the CERN Super Proton Synchrotron (SPS). This procedure allows to unfold the moments of the unknown multiplicity distributions of protons (p), kaons (K), pions ($π$) and electrons (e). Using these moments the excitation function of the fluctuation measure $ν_{\text{\text{dyn}}}$[A,B] was measured, with A and B denoting different particle types. The obtained energy dependence of $ν_{\text{dyn}}$ agrees with previously published NA49 results on the related measure $σ_{\text{dyn}}$. Moreover, $ν_{\text{dyn}}$ was found to depend on the phase space coverage for [K,p] and [K,$π$] pairs. This feature most likely explains the reported differences between measurements of NA49 and those of STAR in central Au+Au collisions.

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

Hadron calorimeter with MAPD readout in the NA61/SHINE experiment

The modular hadron calorimeter with micro-pixel avalanche photodiodes readout for the NA61/SHINE experiment at the CERN SPS is presented. The calorimeter consists of 44 independent modules with lead-scintillator sandwich structure. The light from the scintillator tiles is captured by and transported with WLS-fibers embedded in scintillator grooves. The construction provides a longitudinal segmentation of the module in 10 sections with independent MAPD readout. MAPDs with pixel density of $~10^{4}$/mm$^2$ ensure good linearity of calorimeter response in a wide dynamical range. The performance of the calorimeter prototype in a beam test is reported.

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

Pion emission from the T2K replica target: method, results and application

The T2K long-baseline neutrino oscillation experiment in Japan needs precise predictions of the initial neutrino flux. The highest precision can be reached based on detailed measurements of hadron emission from the same target as used by T2K exposed to a proton beam of the same kinetic energy of 30 GeV. The corresponding data were recorded in 2007-2010 by the NA61/SHINE experiment at the CERN SPS using a replica of the T2K graphite target. In this paper details of the experiment, data taking, data analysis method and results from the 2007 pilot run are presented. Furthermore, the application of the NA61/SHINE measurements to the predictions of the T2K initial neutrino flux is described and discussed.

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

Energy dependence of kaon-to-proton ratio fluctuations in central Pb+Pb collisions from $\sqrt{s_{NN}}$ = 6.3 to 17.3 GeV

Kaons and protons carry large parts of two conserved quantities, strangeness and baryon number. It is argued that their correlation and thus also fluctuations are sensitive to conditions prevailing at the anticipated parton-hadron phase boundary. Fluctuations of the $(\mathrm{K}^+ + \mathrm{K}^-)/(\mathrm{p}+\bar{\mathrm{p}})$ and $\mathrm{K}^+/\mathrm{p}$ ratios have been measured for the first time by NA49 in central Pb+Pb collisions at 5 SPS energies between $\sqrt{s_{NN}}$= 6.3 GeV and 17.3 GeV. Both ratios exhibit a change of sign in $σ_{\mathrm{dyn}}$, a measure of non-statistical fluctuations, around $\sqrt{s_{NN}}$ = 8 GeV. Below this energy, $σ_{\mathrm{dyn}}$ is positive, indicating higher fluctuation compared to a mixed event background sample, while for higher energies, $σ_{\mathrm{dyn}}$ is negative, indicating correlated emission of kaons and protons. The results are compared to UrQMD calculations which which give a good description at the higher SPS energies, but fail to reproduce the transition to positive values.

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

Measurement of Production Properties of Positively Charged Kaons in Proton-Carbon Interactions at 31 GeV/c

Spectra of positively charged kaons in p+C interactions at 31 GeV/c were measured with the NA61/SHINE spectrometer at the CERN SPS. The analysis is based on the full set of data collected in 2007 with a graphite target with a thickness of 4% of a nuclear interaction length. Interaction cross sections and charged pion spectra were already measured using the same set of data. These new measurements in combination with the published ones are required to improve predictions of the neutrino flux for the T2K long baseline neutrino oscillation experiment in Japan. In particular, the knowledge of kaon production is crucial for precisely predicting the intrinsic electron neutrino component and the high energy tail of the T2K beam. The results are presented as a function of laboratory momentum in 2 intervals of the laboratory polar angle covering the range from 20 up to 240 mrad. The kaon spectra are compared with predictions of several hadron production models. Using the published pion results and the new kaon data, the K+/π+ ratios are computed.

preprint2011arXiv

Measurements of Cross Sections and Charged Pion Spectra in Proton-Carbon Interactions at 31 GeV/c

Interaction cross sections and charged pion spectra in p+C interactions at 31 GeV/c were measured with the large acceptance NA61/SHINE spectrometer at the CERN SPS. These data are required to improve predictions of the neutrino flux for the T2K long baseline neutrino oscillation experiment in Japan. A set of data collected during the first NA61/SHINE run in 2007 with an isotropic graphite target with a thickness of 4% of a nuclear interaction length was used for the analysis. The measured p+C inelastic and production cross sections are 257.2 +- 1.9 +- 8.9 mb and 229.3 +- 1.9 +- 9.0 mb, respectively. Inclusive production cross sections for negatively and positively charged pions are presented as a function of laboratory momentum in 10 intervals of the laboratory polar angle covering the range from 0 up to 420 mrad. The spectra are compared with predictions of several hadron production models.

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

Proton -- Lambda Correlations in Central Pb+Pb Collisions at sqrt(s_{NN}) = 17.3 GeV

The momentum correlation between protons and lambda particles emitted from central Pb+Pb collisions at sqrt(s_{NN}) = 17.3 GeV was studied by the NA49 experiment at the CERN SPS. A clear enhancement is observed for small relative momenta (q_{inv} < 0.2 GeV). By fitting a theoretical model, which uses the strong interaction between the proton and the lambda in a given pair, to the measured data a value for the effective source size is deduced. Assuming a static Gaussian source distribution we derive an effective radius parameter of R_G = 3.02 \pm 0.20$(stat.)^{+0.44}_{-0.16}(syst.) fm.

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

Systematics of azimuthal asymmetries in heavy ion collisions in the 1 A GeV regime

Using the large acceptance apparatus FOPI, we study central and semi-central collisions in the reactions (energies in A GeV are given in parentheses): 40Ca+40Ca (0.4, 0.6, 0.8, 1.0, 1.5, 1.93), 58Ni+58Ni (0.15, 0.25, 0.4), 96Ru+96Ru (0.4, 1.0, 1.5), 96Zr+96Zr (0.4, 1.0, 1.5), 129Xe+CsI (0.15, 0.25, 0.4), 197Au+197Au (0.09, 0.12, 0.15, 0.25, 0.4, 0.6, 0.8, 1.0, 1.2, 1.5). The observables include directed and elliptic flow. The data are compared to earlier data where possible and to transport model simulations. A stiff nuclear equation of state is found to be incompatible with the data. Evidence for extra-repulsion of neutrons in compressed asymmetric matter is found.

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

Centrality dependence of proton and antiproton spectra in Pb+Pb collisions at 40A GeV and 158A GeV measured at the CERN SPS

The yields of (anti-)protons were measured by the NA49 Collaboration in centrality selected Pb+Pb collisions at 40A GeV and 158A GeV. Particle identification was obtained in the laboratory momentum range from 5 to 63 GeV/c by the measurement of the energy loss dE/dx in the TPC detector gas. The corresponding rapidity coverage extends 1.6 units from mid-rapidity into the forward hemisphere. Transverse mass spectra, the rapidity dependences of the average transverse mass, and rapidity density distributions were studied as a function of collision centrality. The values of the average transverse mass as well as the midrapidity yields of protons when normalized to the number of wounded nucleons show only modest centrality dependences. In contrast, the shape of the rapidity distribution changes significantly with collision centrality, especially at 40A GeV. The experimental results are compared to calculations of the HSD and UrQMD transport models.

preprint2010arXiv

Inclusive production of charged kaons in p+p collisions at 158 GeV/c beam momentum and a new evaluation of the energy dependence of kaon production up to collider energies

New data on the production of charged kaons in p+p interactions are presented. The data come from a sample of 4.8 million inelastic events obtained with the NA49 detector at the CERN SPS at 158 GeV/c beam momentum. The kaons are identified by energy loss in a large TPC tracking system. Inclusive invariant cross sections are obtained in intervals from 0 to 1.7 GeV/c in transverse momentum and from 0 to 0.5 in Feynman x. Using these data as a reference, a new evaluation of the energy dependence of kaon production, including neutral kaons, is conducted over a range from 3 GeV to p+anti-p collider energies.

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

Systematics of central heavy ion collisions in the 1A GeV regime

Using the large acceptance apparatus FOPI, we study central collisions in the reactions (energies in A GeV are given in parentheses): 40Ca+40Ca (0.4, 0.6, 0.8, 1.0, 1.5, 1.93), 58Ni+58Ni (0.15, 0.25, 0.4), 96Ru+96Ru (0.4, 1.0, 1.5), 96Zr+96Zr (0.4, 1.0, 1.5), 129Xe+CsI (0.15, 0.25, 0.4), 197Au+197Au (0.09, 0.12, 0.15, 0.25, 0.4, 0.6, 0.8, 1.0, 1.2, 1.5). The observables include cluster multiplicities, longitudinal and transverse rapidity distributions and stopping, and radial flow. The data are compared to earlier data where possible and to transport model simulations.

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

Inclusive production of protons, anti-protons and neutrons in p+p collisions at 158 GeV/c beam momentum

New data on the production of protons, anti-protons and neutrons in p+p interactions are presented. The data come from a sample of 4.8 million inelastic events obtained with the NA49 detector at the CERN SPS at 158 GeV/c beam momentum. The charged baryons are identified by energy loss measurement in a large TPC tracking system. Neutrons are detected in a forward hadronic calorimeter. Inclusive invariant cross sections are obtained in intervals from 0 to 1.9 GeV/c (0 to 1.5 GeV/c) in transverse momentum and from -0.05 to 0.95 (-0.05 to 0.4) in Feynman x for protons (anti-protons), respectively. pT integrated neutron cross sections are given in the interval from 0.1 to 0.9 in Feynman x. The data are compared to a wide sample of existing results in the SPS and ISR energy ranges as well as to proton and neutron measurements from HERA and RHIC.

preprint2009arXiv

QPACE -- a QCD parallel computer based on Cell processors

QPACE is a novel parallel computer which has been developed to be primarily used for lattice QCD simulations. The compute power is provided by the IBM PowerXCell 8i processor, an enhanced version of the Cell processor that is used in the Playstation 3. The QPACE nodes are interconnected by a custom, application optimized 3-dimensional torus network implemented on an FPGA. To achieve the very high packaging density of 26 TFlops per rack a new water cooling concept has been developed and successfully realized. In this paper we give an overview of the architecture and highlight some important technical details of the system. Furthermore, we provide initial performance results and report on the installation of 8 QPACE racks providing an aggregate peak performance of 200 TFlops.

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.

preprint2008arXiv

Status of the QPACE Project

We give an overview of the QPACE project, which is pursuing the development of a massively parallel, scalable supercomputer for LQCD. The machine is a three-dimensional torus of identical processing nodes, based on the PowerXCell 8i processor. The nodes are connected by an FPGA-based, application-optimized network processor attached to the PowerXCell 8i processor. We present a performance analysis of lattice QCD codes on QPACE and corresponding hardware benchmarks.

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.

preprint2001arXiv

Ultra high energy cosmic rays: clustering, GUT scale and neutrino masses

The clustering of ultra high energy (above 5\cdot 10^{19} eV) cosmic rays (UHECR) suggests that they might be emitted by compact sources. We present a statistical analysis on the source density based on the multiplicities. The propagation of UHECR protons is studied in detail. The UHECR spectrum is consistent with the decay of GUT scale particles and/or with the Z-burst. The predicted GUT mass is m_X=10^b GeV, where b=14.6_{-1.7}^{+1.6}. Our neutrino mass prediction depends on the origin of the power part of the spectrum: m_ν=2.75^{+1.28}_{-0.97} eV for halo and 0.26^{+0.20}_{-0.14} eV for extragalactic (EG) origin.

preprint1999arXiv

Electroweak Phase Transitions

Recent developments on the four dimensional (4d) lattice studies of the finite temperature electroweak phase transition (EWPT) are summarized. The phase diagram is given in the continuum limit. The finite temperature SU(2)-Higgs phase transition is of first order for Higgs-boson masses m_H<66.5 +/- 1.4 GeV. Above this endpoint only a rapid cross-over can be seen. The full 4d result agrees completely with that of the dimensional reduction approximation. The Higgs-boson endpoint mass in the Standard Model (SM) would be 72.1 +/- 1.4 GeV. Taking into account the LEP Higgs-boson mass lower bound excludes any EWPT in the SM. A one-loop calculation of the static potential in the SU(2)-Higgs model enables a precise comparison between lattice simulations and perturbative results. The most popular extension of the SM, the Minimal Supersymmetric SM (MSSM) is also studied on 4d lattices.

preprint1998arXiv

Perturbative and non-perturbative studies of the SU(2)-Higgs model on lattices with asymmetric lattice spacings

We present a calculation of the O(g^2,λ) perturbative corrections to the coupling anisotropies of the SU(2)-Higgs model on lattices with asymmetric lattice spacings. These corrections are obtained by a one-loop calculation requiring the rotational invariance of the gauge and Higgs boson propagators in the continuum limit. The coupling anisotropies are also determined from numerical simulations of the model on appropriate lattices. The one-loop perturbation theory and the simulation results agree with high accuracy. It is demonstrated that rotational invariance is also restored for the static potential determined from space-space and space-time Wilson loops.

preprint1995arXiv

PERTURBATIVE STUDY OF THE ELECTROWEAK PHASE TRANSITION

The electroweak phase transition is studied at finite temperature. The effective action is given to higher orders, including wave function correction factors and the full $g^4,λ^2$ effective potential. An upper bound for the Higgs mass $m_{H} \approx 70\ GeV$ is concluded for the reliability of the perturbative approach. A gauge invariant treatment of the phase transition is presented.

preprint1995arXiv

THERMODYNAMICS OF THE ELECTROWEAK PHASE TRANSITION

We discuss several general aspects of the free energy of the standard model at high temperatures. In particular the Clausius-Clapeyron equation is shown to yield a relation between the latent heat and the jump in the order parameter. The free energy is calculated as function of temperature in resummed perturbation theory to two-loop order. A new resummation procedure is proposed in which the symmetric phase and the Higgs phase are treated differently. A quantitative description of the phase transition is achieved for Higgs masses below $\sim 70$ GeV. The results are found to be in agreement with recent numerical simulations on large lattices. The phase transition provides no evidence for strong non-perturbative effects in the symmetric phase.