Search NASA⌕ Search

Engineering topics

Leskovec, Luka

Publications and source records attributed to Leskovec, Luka.

Low-energy pion-nucleon scattering and the Δ resonance in lattice QCD

We report on our investigation of low-energy pion-nucleon scattering from lattice QCD with Wilson type fermions. Our focus is on the isospin I = 3/2 channel, which couples to the Delta resonance. We discuss our calculations aiming at the extraction of the Delta resonance mass, width and effective coupling at pion mass 250 MeV, using the Lüscher finite-volume method to extract the phase shift δ J=3/2,l=1 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Calculating the $\rho$ radiative decay width with lattice QCD

We present the results of our lattice QCD study of the $\pi\gamma\to\pi\pi$ process, where the $\rho$ resonance appears as an enhancement in the transition amplitude. We use $N_f=2+1$ clover fermions on a lattice of $L=3.6$ fm and a pion mass of $320$ MeV. Using a combination of forward, stochastic, and sequential propagators, we calculate the two-point and three-point functions that allow us to determine the $\pi\gamma\to\pi\pi$ matrix elements for several values of the invariant mass $s$ and momentum transfer $q^2$. To fit the $q^2$ and $s$ dependence of the $\pi\gamma\to\pi\pi$ amplitude, we explore a set of general parametrizations based on a Taylor expansion. By analytic continuation to the complex pole corresponding to the $\rho$ resonance, we determine the resonant form factors and calculate the radiative decay width of the $\rho$.

Leskovec, Luka↗

P -wave π π scattering and the ρ resonance from lattice QCD

Here we calculate the parameters describing elastic I = 1, P-wave ππ scattering using lattice QCD with 2 +1 flavors of clover fermions. Our calculation is performed with a pion mass of m π ≈ 320 MeV and a lattice size of L ≈ 3.6 fm. We construct the two-point correlation matrices with both quark-antiquark and two-hadron interpolating fields using a combination of smeared forward, sequential and stochastic propagators. The spectra in all relevant irreducible representations for total momenta |$\vec{P}$| ≤ $\sqrt{3}\frac{2π}{L}$ are extracted with two alternative methods: a variational analysis as well as multiexponential matrix fits. We perform an analysis using Lüscher’s formalism for the energies below the inelastic thresholds, and investigate several phase shift models, including possible nonresonant contributions. We find that our data are well described by the minimal Breit-Wigner form, with no statistically significant nonresonant component. In determining the ρ resonance mass and coupling we compare two different approaches: fitting the individually extracted phase shifts versus fitting the t-matrix model directly to the energy spectrum. We find that both methods give consistent results, and at a pion mass of am π = 0.18295 (36) stat obtain gρππ = 5.69(13) stat (16) sys , am ρ = 0.4609(16) stat (14) sys , and am ρ /am N = 0.7476(38) stat (23) sys , where the first uncertainty is statistical and the second is the systematic uncertainty due to the choice of fit ranges.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Positive Parity $D_s$ Mesons

We study the positive parity charmed strange mesons using lattice QCD, the only reliable ab initio method to study QCD at low energies. Especially the experimentally observed $D_{s0}^*(2317)$ and $D_{s1}(2460)$ have challenged theory for quite some time. The dynamical lattice QCD simulations are performed at two distinct pion masses, $m_{\pi}$ = 266 MeV and 156 MeV, using both $\bar{c}s$ as well as $DK$ and $D^*K$ scattering operators in the construction of the correlation matrix in order to take into the account threshold effects. While the $J^P = 0^+$ channel benefited most from the inclusion of scattering operators, it was also crucial for the case of the $D_{s1}(2460)$. Using the L\"uscher method, which relates the finite volume spectrum to the infinite volume scattering matrix, we were able to determine the near threshold behavior of the scattering amplitude. From it we extracted the binding momenta and the masses of the below-threshold bound states $D_{s0}^*(2317)$ and $D_{s1}(2460)$ by determining the pole positions of the scattering amplitudes. Our results compare well with experiment, resolving a long standing discrepancy between theory and experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Vector and scalar charmonium resonances with lattice QCD

We perform an exploratory lattice QCD simulation of DD¯ scattering, aimed at determining the masses as well as the decay widths of charmonium resonances above open charm threshold. Neglecting coupling to other channels, the resulting phase shift for DD¯ scattering in p-wave yields the well-known vector resonance ψ(3770). For m π = 156 MeV, the extracted resonance mass and the decay width agree with experiment within large statistical uncertainty. The scalar charmonium resonances present a puzzle, since only the ground state Χ c0 (1P) is well understood, while there is no commonly accepted candidate for its first excitation. We simulate DD¯ scattering in s-wave in order to shed light on this puzzle. The resulting phase shift supports the existence of a yet-unobserved narrow resonance with a mass slightly below 4 GeV. A scenario with this narrow resonance and a pole at Χ c0 (1P) agrees with the energy-dependence of our phase shift. In addition, further lattice QCD simulations and experimental efforts are needed to resolve the puzzle of the excited scalar charmonia.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

$DK$ and $D^* K$ scattering near threshold

We study the three $D_s$ quantum channels $J^P = 0^+$, $1^+$ and $2^+$ where experiments have identified the charm-strange states $D^*_{s0} (2317)$, $D_{s1}(2460)$, $D_{s1}(2536)$ near the $DK$ and $D^*K$ thresholds, and $D^*_{s2}(2573)$. We consider correlation functions for sets of $\overline q q$ operators and, for $J^P = 0^+$, $1^+$, also the $DK$ and $D^*K$ meson-meson interpolators and determine for these cases values of the elastic scattering amplitude. Constructing the full set of correlators requires propagators which connect any pair of lattice sites. For one ensemble of gauge configurations ($32^3\times 64$, $m_\pi\approx 156$ MeV) a stochastic distillation variant is employed and for another ensemble ($16^3\times 32$, $m_\pi\approx 266$ MeV) we use the full distillation method. Both, $D^*_{s0} (2317)$ and $D_{s1}(2460)$, are found as bound states below threshold, whereas $D_{s1}(2536)$, and $D^*_{s2}(2573)$ are identified as narrow resonances close to the experimental masses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Study of the $Z_c^+$ channel in lattice QCD

Several charged charmonium-like hadrons called $Z_c$ have been recently discovered by different experiments. In contrast to conventional hadrons these contain at least two valence quarks and antiquarks ($\bar{c}c\bar{d}u$). We perform a lattice QCD simulation of the $I^G(J^{PC})=1^+(1^{+-})$ channel including all relevant two-meson operators under 4.3 GeV: $J/\psi \pi$, $\psi_{2S}\pi$, $\psi_{1D}\pi$, $D \bar{D}^*$, $D^* \bar{D}^*$, $\eta_c \rho$ as well as additional diquark anti-diquark operators. In our $N_f = 2$ simulation with pion mass at 266 MeV we are able to identify all two-meson levels within the energy region of interest. However we find no additional level identifiable as a candidate for $Z_c$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Axial resonances a$_{1}$(1260), b$_{1}$(1235) and their decays from the lattice

The light axial-vector resonances $a_1(1260)$ and $b_1(1235)$ are explored in Nf=2 lattice QCD by simulating the corresponding scattering channels $\rho\pi$ and $\omega\pi$. Interpolating fields $\bar{q} q$ and $\rho\pi$ or $\omega\pi$ are used to extract the s-wave phase shifts for the first time. The $\rho$ and $\omega$ are treated as stable and we argue that this is justified in the considered energy range and for our parameters $m_\pi\simeq 266~$MeV and $L\simeq 2~$fm. We neglect other channels that would be open when using physical masses in continuum. Assuming a resonance interpretation a Breit-Wigner fit to the phase shift gives the $a_1(1260)$ resonance mass $m_{a1}^{res}=1.435(53)(^{+0}_{-109})$ GeV compared to $m_{a1}^{exp}=1.230(40)$ GeV. The $a_1$ width $\Gamma_{a1}(s)=g^2 p/s$ is parametrized in terms of the coupling and we obtain $g_{a_1\rho\pi}=1.71(39)$ GeV compared to $g_{a_1\rho\pi}^{exp}=1.35(30)$ GeV derived from $\Gamma_{a1}^{exp}=425(175)$ MeV. In the $b_1$ channel, we find energy levels related to $\pi(0)\omega(0)$ and $b_1(1235)$, and the lowest level is found at $E_1 \gtrsim m_\omega+m_\pi$ but is within uncertainty also compatible with an attractive interaction. Lastly, assuming the coupling $g_{b_1\omega\pi}$ extracted from the experimental width we estimate $m_{b_1}^{res}=1.414(36)(^{+0}_{-83})$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗