• ### Determination of N* amplitudes from associated strangeness production in p+p collisions(1703.01978)

Sept. 6, 2018 nucl-ex
We present the first determination of the energy-dependent production amplitudes of N$^{*}$ resonances with masses between 1650 MeV/c$^{2}$ and 1900 MeV/c$^{2}$ for an excess energy between $0$ and $600$ MeV. A combined Partial Wave Analysis of seven exclusively reconstructed data samples for the reaction p+p $\rightarrow pK\Lambda$ measured by the COSY-TOF, DISTO, FOPI and HADES collaborations in fixed target experiments at kinetic energies between 2.14 and 3.5 GeV is used to determine the amplitude of the resonant and non-resonant contributions.
• ### Measurements of $\pi^\pm$, K$^\pm$, p and $\bar{\textrm{p}}$ spectra in proton-proton interactions at 20, 31, 40, 80 and 158 GeV/c with the NA61/SHINE spectrometer at the CERN SPS(1705.02467)

Sept. 27, 2017 nucl-ex
Measurements of inclusive spectra and mean multiplicities of $\pi^\pm$, K$^\pm$, p and $\bar{\textrm{p}}$ 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) were performed at the CERN Super Proton Synchrotron using the large acceptance NA61/SHINE hadron spectrometer. Spectra are presented as function of rapidity and transverse momentum and are compared to predictions of current models. The measurements serve as the baseline in the NA61/SHINE study of the properties of the onset of deconfinement and search for the critical point of strongly interacting matter.
• ### Measurement of Meson Resonance Production in $\pi^{-} +$C Interactions at SPS energies(1705.08206)

May 23, 2017 nucl-ex, astro-ph.HE
We present measurements of $\rho^0$, $\omega$ and K$^{*0}$ spectra in $\pi^{-} +$C production interactions at 158 GeV/c and $\rho^0$ spectra at 350 GeV/c using the NA61/SHINE spectrometer at the CERN SPS. Spectra are presented as a function of the Feynman's variable $x_\text{F}$ in the range $0 < x_\text{F} < 1$ and $0 < x_\text{F} < 0.5$ for 158 GeV/c and 350 GeV/c respectively. Furthermore, we show comparisons with previous measurements and predictions of several hadronic interaction models. These measurements are essential for a better understanding of hadronic shower development and for improving the modeling of cosmic ray air showers.
• ### Two-particle correlations in azimuthal angle and pseudorapidity in inelastic p+p interactions at the CERN Super Proton Synchrotron(1610.00482)

Feb. 7, 2017 nucl-ex
Results on two-particle $\Delta\eta\Delta\phi$ correlations in inelastic p+p interactions at 20, 31, 40, 80, and 158~GeV/c are presented. The measurements were performed using the large acceptance NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. The data show structures which can be attributed mainly to effects of resonance decays, momentum conservation, and quantum statistics. The results are compared with the EPOS and UrQMD models.
• ### Workshop on Excited Hyperons in QCD Thermodynamics at Freeze-Out (YSTAR2016) Mini-Proceedings(1701.07346)

Feb. 2, 2017 hep-ph
This Workshop brought top experts, researchers, postdocs, and students from high-energy heavy ion interactions, lattice QCD and hadronic physics communities together. YSTAR2016 discussed the impact of "missing" hyperon resonances on QCD thermodynamics, on freeze-out in heavy ion collisions, on the evolution of early universe, and on the spectroscopy of strange particles. Recent studies that compared lattice QCD predictions of thermodynamic properties of quark-gluon plasma at freeze-out with calculations based on statistical hadron resonance gas models as well as experimentally measured ratios between yields of different hadron species in heavy ion collisions provide indirect evidence for the presence of "missing" resonances in all of these contexts. The aim of the YSTAR2016 Workshop was to sharpen these comparisons and advance our understanding of the formation of strange hadrons from quarks and gluons microseconds after the Big Bang and in todays experiments at LHC and RHIC as well as at future facilities like FAIR, J-PARC and KL at JLab. It was concluded that the new initiative to create a secondary beam of neutral kaons at JLab will make a bridge between the hardron spectroscopy, heavy-ion experiments and lattice QCD studies addressing some major issues related to thermodynamics of the early universe and cosmology in general.
• ### Up and down quark masses and corrections to Dashen's theorem from lattice QCD and quenched QED(1702.00309)

Feb. 1, 2017 hep-lat
We present a determination of the corrections to Dashen's theorem and of the individual up and down quark masses from a lattice calculation based on quenched QED and $N_f=2+1$ QCD simulations with 5 lattice spacings down to 0.054 fm. The simulations feature lattice sizes up to 6 fm and average up-down quark masses all the way down to their physical value. For the parameter which quantifies violations to Dashens's theorem we obtain $\epsilon=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) \, MeV$ and $m_d=4.67(6)(5)(4) \, MeV$ in the $\overline{MS}$ scheme at $2 \, 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.
• ### Slope and curvature of the hadron vacuum polarization at vanishing virtuality from lattice QCD(1612.02364)

Dec. 7, 2016 hep-ph, hep-lat
We compute the slope and curvature, at vanishing four-momentum transfer squared, of the leading order hadron vacuum polarization function, using lattice QCD. Calculations are performed with 2+1+1 flavors of staggered fermions directly at the physical values of the quark masses and in volumes of linear extent larger than 6fm. The continuum limit is carried out using six different lattice spacings. All connected and disconnected contributions are calculated, up to and including those of the charm.
• ### Measurements of $\pi^{\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(1603.06774)

Nov. 29, 2016 hep-ex, physics.ins-det
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 $\pi^\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.
• ### Multiplicity and transverse momentum fluctuations in inelastic proton-proton interactions at the CERN Super Proton Synchrotron(1510.00163)

Aug. 30, 2016 hep-ex, nucl-ex
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 \$\Delta[P_{T},N]\$, \$\Sigma[P_{T},N]\$ and \$\Phi_{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.
• ### The QCD equation of state at finite density from analytical continuation(1607.02493)

July 8, 2016 hep-ph, nucl-th, hep-lat
We determine the equation of state of QCD at finite chemical potential, to order $(\mu_B/T)^6$, for a system of 2+1 quark flavors. The simulations are performed at the physical mass for the light and strange quarks on several lattice spacings; the results are continuum extrapolated using lattices of up to $N_t=16$ temporal resolution. The QCD pressure and interaction measure are calculated along the isentropic trajectories in the $(T,~\mu_B)$ plane corresponding to the RHIC Beam Energy Scan collision energies. Their behavior is determined through analytic continuation from imaginary chemical potentials of the baryonic density. We also determine the Taylor expansion coefficients around $\mu_B=0$ from the simulations at imaginary chemical potentials. Strangeness neutrality and charge conservation are imposed, to match the experimental conditions.
• ### Lattice QCD for Cosmology(1606.07494)

June 27, 2016 hep-ph, hep-lat, astro-ph.CO
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.
• ### Production of deuterium, tritium, and $^3$He in central Pb+Pb collisions at 20A, 30A, 40A, 80A, and 158A GeV at the CERN SPS(1606.04234)

June 15, 2016 nucl-ex
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.
• ### Centrality dependence of subthreshold $\phi$ meson production in Ni+Ni collisions at 1.9A GeV(1602.04378)

June 9, 2016 nucl-ex
We analysed the $\phi$ 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 $\phi$ meson production in heavy-ion collisions. The rise of $\phi$ meson multiplicity per event with mean number of participants can be parameterized by the power function with exponent $\alpha = 1.8 \pm 0.6$. The ratio of $\phi$ 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$.
• ### Strange meson production in Al+Al collisions at 1.9A GeV(1512.06988)

May 17, 2016 nucl-ex
The production of K$^+$, K$^-$ and $\varphi$(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 $\phi$ mesons are deduced. The contribution to K$^-$ production from $\phi$ meson decays is found to be [17 $\pm$ 3 (stat) $^{+2}_{-7}$ (syst)] %. The results are in line with previous K$^{\pm}$ and $\phi$ data obtained for different colliding systems at similar incident beam energies.
• ### Lattice computation of the nucleon scalar quark contents at the physical point(1510.08013)

May 11, 2016 hep-ph, hep-lat
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 $\sigma_{\pi N}=38(3)(3)\,\text{MeV}$, $\sigma_{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).
• ### Up and down quark masses and corrections to Dashen's theorem from lattice QCD and quenched QED(1604.07112)

April 25, 2016 hep-ph, hep-lat
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 $\epsilon=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.
• ### Production of $\Lambda$ hyperons in inelastic p+p interactions at 158 GeV/$c$(1510.03720)

April 25, 2016 hep-ex, nucl-ex
Inclusive production of $\Lambda$-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.
• ### Sigma term and strangeness content of octet baryons(1109.4265)

Feb. 27, 2016 hep-ph, hep-lat
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_\pi L > 4), and three values of the lattice spacing. Our main results are the nucleon sigma term \sigma_{\pi N} = 39(4)(^{+18}_{-7}) and the strangeness content y_{N} = 0.20(7)(^{+13}_{-17}). Under the assumption of validity of covariant baryon \chi PT in our range of masses one finds y_{N} = 0.276(77)(^{+90}_{-62}).
• ### Measurements of $\pi^\pm$, $K^\pm$, $K^0_S$, $\Lambda$ and proton production in proton-carbon interactions at 31 GeV/$c$ with the NA61/SHINE spectrometer at the CERN SPS(1510.02703)

Feb. 24, 2016 hep-ex, nucl-ex
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 $\pi^\pm$, $K^\pm$, p, $K^0_S$ and $\Lambda$ 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.
• ### Towards the QCD phase diagram from analytical continuation(1601.00466)

Jan. 4, 2016 nucl-th, hep-lat
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.
• ### The QCD phase diagram from analytic continuation(1507.07510)

Nov. 12, 2015 hep-ph, nucl-th, hep-lat
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 $\mu_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 $\kappa = 0.0149 \pm 0.0021$.
• ### Critical fluctuations of the proton density in A+A collisions at $158A$ GeV(1208.5292)

Nov. 6, 2015 hep-ph, hep-ex, nucl-ex
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 $\phi_{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 $\pi^+ \pi^-$ pairs in the same Si+Si collisions.
• ### Complex Langevin dynamics for dynamical QCD at nonzero chemical potential: a comparison with multi-parameter reweighting(1508.05260)

Oct. 27, 2015 hep-ph, hep-lat
We study lattice QCD at non-vanishing chemical potential using the complex Langevin equation. We compare the results with multi-parameter reweighting both from $\mu=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.
• ### Measurement of event-by-event transverse momentum and multiplicity fluctuations using strongly intensive measures $\Delta[P_T, N]$ and $\Sigma[P_T, N]$ in nucleus-nucleus collisions at the CERN Super Proton Synchrotron(1509.04633)

Sept. 15, 2015 nucl-ex
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 $\Phi_{p_T}$ by a study of the recently proposed strongly intensive measures of fluctuations $\Delta[P_T, N]$ and $\Sigma[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 $\Delta[P_T, N]$ and $\Sigma[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.
• ### Axion cosmology, lattice QCD and the dilute instanton gas(1508.06917)

Aug. 27, 2015 hep-ph, hep-lat
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 $\chi(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 $\chi(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.