Source author record

Luka Leskovec

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

21works
4topics
4close collaborators

Published work

21 publications

preprint2023arXiv

Doubly charmed tetraquark: isospin channels and diquark-antidiquark interpolators

We perform a lattice simulation to investigate the doubly charmed tetraquark Tcc+T^+_{cc} observed by the LHCb collaboration, slightly below the D∗+D0D^{*+}D^0 threshold, with flavor content ccuˉdˉcc\bar{u}\bar{d} and isospin-00. Two-meson interpolators are implemented to explore the isospin quantum numbers I=0I=0 and I=1I=1. We observe attraction near the DD∗DD^* threshold for I=0I=0 and repulsion for I=1I=1. Moreover, we also include diquark-antidiquark interpolators to study their effect on the energy spectrum. There is no significant shift in the ground state energy when adding diquark-antidiquark interpolators to the interpolator basis when the heavy quark mass is close to the physical charm quark mass. However, we observe a non-negligible shift in the second energy level. This effect has to be taken into account to extract the scattering amplitude of the Tcc+T^+_{cc}. Finally, with a higher mass (close to the bottom quark), the ground state is shifted down significantly. The simulation is performed on Nf=2+1N_f=2+1 CLS ensembles with mπ≃280m_π\simeq 280 MeV.

preprint2022arXiv

The πγ→πππγ\to ππ transition and the ρρ radiative decay width from lattice QCD

We report a lattice QCD determination of the πγ→πππγ\to ππ transition amplitude for the PP-wave, I=1I=1 two-pion final state, as a function of the photon virtuality and ππππ invariant mass. The calculation was performed with 2+12+1 flavors of clover fermions at a pion mass of approximately 320320 MeV, on a 323×9632^3 \times 96 lattice with L≈3.6L\approx 3.6 fm. We construct the necessary correlation functions using a combination of smeared forward, sequential and stochastic propagators, and determine the finite-volume matrix elements for all ππππ momenta up to ∣P⃗∣=32πL|\vec{P}|= \sqrt{3} \frac{2π}{L} and all associated irreducible representations. In the mapping of the finite-volume to infinite-volume matrix elements using the Lellouch-Lüscher factor, we consider two different parametrizations of the ππππ scattering phase shift. We fit the q2q^2 and ss dependence of the infinite-volume transition amplitude in a model-independent way using series expansions, and compare multiple different truncations of this series. Through analytic continuation to the ρρ resonance pole, we also determine the πγ→ρπγ\to ρ resonant transition form factor and the ρρ meson photocoupling, and obtain ∣Gρπγ∣=0.0802(32)(20)|G_{ρπγ}| = 0.0802(32)(20).

preprint2021arXiv

I=1/2I=1/2 SS-wave and PP-wave KπKπ scattering and the κκ and K∗K^* resonances from lattice QCD

We present a lattice-QCD determination of the elastic isospin-1/21/2 SS-wave and PP-wave KπKπ scattering amplitudes as a function of the center-of-mass energy using Lüscher's method. We perform global fits of KK-matrix parametrizations to the finite-volume energy spectra for all irreducible representations with total momenta up to 32πL\sqrt{3}\frac{2π}{L}; this includes irreps that mix the SS- and PP-waves. Several different parametrizations for the energy dependence of the KK-matrix are considered. We also determine the positions of the nearest poles in the scattering amplitudes, which correspond to the broad κκ resonance in the SS-wave and the narrow K∗(892)K^*(892) resonance in the PP-wave. Our calculations are performed with 2+12+1 dynamical clover fermions for two different pion masses of 317.2(2.2)317.2(2.2) and 175.9(1.8)175.9(1.8) MeV. Our preferred SS-wave parametrization is based on a conformal map and includes an Adler zero; for the PP-wave we use a standard pole parametrization including Blatt-Weisskopf barrier factors. The SS-wave κκ-resonance pole positions are found to be [0.86(12)−0.309(50) i] GeV\left[0.86(12) - 0.309(50)\,i\right]\:{\rm GeV} at the heavier pion mass and [0.499(55)−0.379(66) i] GeV\left[0.499(55)- 0.379(66)\,i\right]\:{\rm GeV} at the lighter pion mass. The PP-wave K∗K^*-resonance pole positions are found to be [0.8951(64)−0.00250(21) i] GeV\left[ 0.8951(64) - 0.00250(21)\,i \right]\:{\rm GeV} at the heavier pion mass and [0.8718(82)−0.0130(11) i] GeV\left[0.8718(82) - 0.0130(11)\,i\right]\:{\rm GeV} at the lighter pion mass, which corresponds to couplings of gK∗Kπ=5.02(26)g_{K^* Kπ}=5.02(26) and gK∗Kπ=4.99(22)g_{K^* Kπ}=4.99(22), respectively.

preprint2020arXiv

Investigation of Doubly Heavy Tetraquark Systems using Lattice QCD

We search for possibly existent bound states in the heavy-light tetraquark channels with quark content bˉbˉud \bar{b}\bar{b}ud , bˉbˉus \bar{b}\bar{b}us and bˉcˉud \bar{b}\bar{c}ud using lattice QCD. We carry out calculations on several gauge link ensembles with Nf=2+1 N_f=2+1 flavours of domain-wall fermions and consider a basis of local and non-local interpolators. Besides extracting the energy spectrum from the correlation matrices, we also perform a Lüscher analysis to extrapolate our results to infinite volume.

preprint2016arXiv

A study of the radiative transition ππ→πγ∗ππ\to πγ^{*} with lattice QCD

Lattice QCD calculations of radiative transitions between hadrons have in the past been limited to processes of hadrons stable under the strong interaction. Recently developed methods for 1→21\to2 transition matrix elements in a finite volume now enable the determination of radiative decay rates of strongly unstable particles. Our lattice QCD study focuses on the process ππ→πγ∗ππ\to πγ^{*}, where the ρρ meson is present as an enhancement in the cross-section. We use 2+12+1 flavors of clover fermions at a pion mass of approximately 320320 MeV and a lattice size of approximately 3.63.6 fm. The required 22-point and 33-point correlation functions are constructed from a set of forward, sequential and stochastic light quark propagators. In addition to determining the ρρ meson resonance parameters via the Lüscher method, the scattering phase shift is used in conjunction with the 1→21\to2 transition matrix element formalism of Briceño, Hansen and Walker-Loud to compute the ππ→πγ∗ππ\toπγ^{*} amplitude at several values of the momentum transfer and ππππ invariant mass.

preprint2016arXiv

Lattice QCD study of heavy-heavy-light-light tetraquark candidates

We investigate heavy-light four-quark systems udbˉbˉud\bar b \bar b with bottom quarks of finite mass which are treated in the framework of NRQCD. We focus on I(JP)=0(1+)I(J^P)=0(1^+), where we recently found evidence for the existence of a tetraquark state using static bottom quarks. Furthermore, we report on an investigation of the udˉbbˉu \bar d b \bar b four-quark system with quantum numbers I(JP)=1(1+)I(J^P)=1(1^+) again using static bottom quarks.

preprint2016arXiv

Positive parity DsD_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 Ds0∗(2317)D_{s0}^*(2317) and Ds1(2460)D_{s1}(2460) have challenged theory for quite some time. The dynamical lattice QCD simulations are performed at two distinct pion masses, mπm_π = 266 MeV and 156 MeV, using both cˉs\bar{c}s as well as DKDK and D∗KD^*K scattering operators in the construction of the correlation matrix in order to take into the account threshold effects. While the JP=0+J^P = 0^+ channel benefited most from the inclusion of scattering operators, it was also crucial for the case of the Ds1(2460)D_{s1}(2460). Using the Lüscher 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 Ds0∗(2317)D_{s0}^*(2317) and Ds1(2460)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.

preprint2015arXiv

Study of the Zc+Z_c^+ channel using lattice QCD

Recently experimentalists have discovered several charged charmonium-like hadrons Zc+Z_c^+ with unconventional quark content cˉcdˉu\bar cc\bar d u. We perform a search for Zc+Z_c^+ with mass below 4.2 4.2~GeV in the channel IG(JPC)=1+(1+−)I^G(J^{PC})=1^+(1^{+-}) using lattice QCD. The major challenge is presented by the two-meson states J/ψ πJ/ψ\, π, ψ2Sπψ_{2S}π, ψ1Dπψ_{1D}π, DDˉ∗D\bar D^*, D∗Dˉ∗D^*\bar D^*, ηcρη_cρ that are inevitably present in this channel. The spectrum of eigenstates is extracted using a number of meson-meson and diquark-antidiquark interpolating fields. For our pion mass of 266~MeV we find all the expected two-meson states but no additional candidate for Zc+Z_c^+ below 4.2 4.2~GeV. Possible reasons for not seeing an additional eigenstate related to Zc+Z_c^+ are discussed. We also illustrate how a simulation incorporating interpolators with a structure resembling low-lying two-mesons states seems to render a Zc+Z_c^+ candidate, which is however not robust after further two-meson states around 4.2 4.2~GeV are implemented.

preprint2015arXiv

Vector and scalar charmonium resonances with lattice QCD

We perform an exploratory lattice QCD simulation of DDˉD \bar D 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ˉD \bar D scattering in p-wave yields the well-known vector resonance ψ(3770)ψ(3770). For mπ=156m_π= 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)χ_{c0}(1P) is well understood, while there is no commonly accepted candidate for its first excitation. We simulate DDˉD \bar D 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)χ_{c0}(1P) agrees with the energy-dependence of our phase shift. Further lattice QCD simulations and experimental efforts are needed to resolve the puzzle of the excited scalar charmonia.

preprint2014arXiv

DsD_{s} mesons with DKDK and D∗KD^{*}K scattering near threshold

DsD_s mesons are studied in three quantum channels (JP=0+J^P=0^+, 1+1^+ and 2+2^+), where experiments have identified the very narrow Ds0∗(2317)D_{s0}^*(2317), Ds1(2460)D_{s1}(2460) and narrow Ds1(2536)D_{s1}(2536), Ds2∗(2573)D_{s2}^*(2573). We explore the effect of nearby DKDK and D∗KD^*K thresholds on the subthreshold states using lattice QCD. Our simulation is done on two very different ensembles of gauge configurations (2 or 2+1 dynamical quarks, Pion mass of 266 or 156 MeV, lattice size 163×3216^3\times 32 or 323×6432^3\times 64). In addition to qˉq\bar{q}q operators we also include meson-meson interpolators in the correlation functions. This clarifies the identification of the states above and below the scattering thresholds. The ensemble with mπ≃156 m_π\simeq 156~MeV renders the Ds1(2460)D_{s1}(2460) as a strong interaction bound state 44(10)MeV below D∗KD^*K threshold, which is in agreement with the experiment. The Ds0∗(2317)D_{s0}^*(2317) is found 37(17)MeV below DKDK threshold, close to experiment value of 45MeV. The narrow resonances Ds1(2536)D_{s1}(2536) and Ds2∗(2573)D_{s2}^*(2573) are also found close to the experimental masses.

preprint2014arXiv

DKDK and D∗KD^* K scattering near threshold

We study the three DsD_s quantum channels JP=0+J^P = 0^+, 1+1^+ and 2+2^+ where experiments have identified the charm-strange states Ds0∗(2317)D^*_{s0} (2317), Ds1(2460)D_{s1}(2460), Ds1(2536)D_{s1}(2536) near the DKDK and D∗KD^*K thresholds, and Ds2∗(2573)D^*_{s2}(2573). We consider correlation functions for sets of q‾q\overline q q operators and, for JP=0+J^P = 0^+, 1+1^+, also the DKDK and D∗KD^*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 (323×6432^3\times 64, mπ≈156m_π\approx 156 MeV) a stochastic distillation variant is employed and for another ensemble (163×3216^3\times 32, mπ≈266m_π\approx 266 MeV) we use the full distillation method. Both, Ds0∗(2317)D^*_{s0} (2317) and Ds1(2460)D_{s1}(2460), are found as bound states below threshold, whereas Ds1(2536)D_{s1}(2536), and Ds2∗(2573)D^*_{s2}(2573) are identified as narrow resonances close to the experimental masses.

preprint2014arXiv

Axial resonances a1(1260), b1(1235) and their decays from the lattice

The light axial-vector resonances a1(1260)a_1(1260) and b1(1235)b_1(1235) are explored in Nf=2 lattice QCD by simulating the corresponding scattering channels ρπρπ and ωπωπ. Interpolating fields qˉq\bar{q} q and ρπρπ or ωπωπ are used to extract the s-wave phase shifts for the first time. The ρρ and ωω are treated as stable and we argue that this is justified in the considered energy range and for our parameters mπ≃266 m_π\simeq 266~MeV and L≃2 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 a1(1260)a_1(1260) resonance mass ma1res=1.435(53)(−109+0)m_{a1}^{res}=1.435(53)(^{+0}_{-109}) GeV compared to ma1exp=1.230(40)m_{a1}^{exp}=1.230(40) GeV. The a1a_1 width Γa1(s)=g2p/sΓ_{a1}(s)=g^2 p/s is parametrized in terms of the coupling and we obtain ga1ρπ=1.71(39)g_{a_1ρπ}=1.71(39) GeV compared to ga1ρπexp=1.35(30)g_{a_1ρπ}^{exp}=1.35(30) GeV derived from Γa1exp=425(175)Γ_{a1}^{exp}=425(175) MeV. In the b1b_1 channel, we find energy levels related to π(0)ω(0)π(0)ω(0) and b1(1235)b_1(1235), and the lowest level is found at E1≳mω+mπE_1 \gtrsim m_ω+m_π but is within uncertainty also compatible with an attractive interaction. Assuming the coupling gb1ωπg_{b_1ωπ} extracted from the experimental width we estimate mb1res=1.414(36)(−83+0)m_{b_1}^{res}=1.414(36)(^{+0}_{-83}).

preprint2014arXiv

K pi scattering and the K* decay width from lattice QCD

K* mesons and in particular the K*(892) were frequently addressed in lattice simulations, but always while ignoring that the K*(892) decays strongly. We present an exploratory extraction of the masses and widths for the K* resonances by simulating K pi scattering in p-wave with I = 1/2 on the lattice. The K pi system with total momenta P = 2*pi/L e_z, 2*pi/L (e_x + e_y) and 0, that allows the extraction of phase shifts at several values of K pi relative momenta, is studied. A Breit-Wigner fit of the phase renders a K*(892) resonance mass m^{lat} = 891 +/- 14 MeV and the K*(892) -> K pi coupling g^{lat} = 5.7 +/- 1.6 compared to the experimental values m^{exp} = 892 MeV and g^{exp} = 5.72 +/- 0.06, where g parametrizes the K* -> K pi width. When extracting the phase shift around the K*(1410) and K2*(1430) resonances we take into account the mixing of p-wave with d-wave and assume that the scattering is elastic in our simulation. This gives us an estimate of the K*(1410) resonance mass m^{lat} = 1.33 +/- 0.02 GeV compared to m^{exp} = 1.414 +/- 0.0015 GeV assuming the experimental K*(1410) -> K pi coupling. We contrast the resonant I = 1/2 channel with the repulsive non-resonant I = 3/2 channel, where the phase is found to be negative and small, in agreement with experiment.

preprint2014arXiv

Study of the Zc+ channel in lattice QCD

Several charged charmonium-like hadrons called ZcZ_c have been recently discovered by different experiments. In contrast to conventional hadrons these contain at least two valence quarks and antiquarks (cˉcdˉu\bar{c}c\bar{d}u). We perform a lattice QCD simulation of the IG(JPC)=1+(1+−)I^G(J^{PC})=1^+(1^{+-}) channel including all relevant two-meson operators under 4.3 GeV: J/ψπJ/ψπ, ψ2Sπψ_{2S}π, ψ1Dπψ_{1D}π, DDˉ∗D \bar{D}^*, D∗Dˉ∗D^* \bar{D}^*, ηcρη_c ρ as well as additional diquark anti-diquark operators. In our Nf=2N_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 ZcZ_c.

preprint2013arXiv

Ds0∗(2317)D_{s0}^*(2317) Meson and DD-Meson-Kaon Scattering from Lattice QCD

The scalar meson Ds0∗(2317)D_{s0}^*(2317) is found 37(17)MeV below DK threshold in a lattice simulation of the JP=0+J^P=0^+ channel using, for the first time, both DK as well as sˉc\bar sc interpolating fields. The simulation is done on Nf=2+1N_f=2+1 gauge configurations with mπ≃156m_π\simeq 156 MeV, and the resulting MDs0∗−14(MDs+3MDs∗)=266(16)M_{D_{s0}^*}-\tfrac{1}{4}(M_{D_s}+3M_{D_s^*})=266(16) MeV is close to the experimental value 241.5(0.8)MeV. The energy level related to the scalar meson is accompanied by additional discrete levels due to DK scattering states. The levels near threshold lead to the negative DK scattering length a0=−1.33(20)a_0=-1.33(20) fm that indicates the presence of a state below threshold.

preprint2013arXiv

Charmonium-like states from scattering on the lattice

Three charmonium-like states were studied using lattice QCD. The candidate for X(3872) was found slightly below DD* threshold in the channel with J^PC=1^++ and I=0, where cc as well as DD* and J/ψωJ/ψω interpolating operators were used. A charmonium-like channel with J^PC=1^+- and I=1 was also studied, as the recently discovered Z_c^+(3900) might reside there. Here J/ψπJ/ψπ and DD* scattering states were found, but no candidate for the Z_c^+(3900). We present also preliminary results for the J^PC=0^++ charmonium channel, where cc, DD and J/ψωJ/ψω interpolating operators were used. A candidate for a resonance, chi_c0^', that couples to DD in J^PC=0^++ was found.

preprint2013arXiv

Evidence for X(3872) from DD* scattering on the lattice

A candidate for the charmonium(like) state X(3872) is found 11 +/- 7 MeV below the DD* threshold using dynamical Nf=2 lattice simulation with J^PC=1^++ and I=0. This is the first lattice simulation that establishes a candidate for X(3872) in addition to the nearby scattering states D D* and J/psi omega, which inevitably have to be present in dynamical QCD. We extract large and negative DD* scattering length a_0=-1.7 +/- 0.4 fm and the effective range r_0=0.5 +/- 0.1 fm, but their reliable determination will have to wait for a simulation on a larger volume. In I=1 channel, only the D D* and J/psi rho scattering states are found and no candidate for X(3872). This is in agreement with the interpretation that X(3872) is dominantly I=0, while its small I=1 component arises solely from the isospin breaking and is therefore absent in our simulation with m_u=m_d.

preprint2013arXiv

K pi scattering in moving frames

We extend our study of the KπKπ system to moving frames and present an exploratory extraction of the masses and widths for the K∗K^* resonances by simulating KπKπ scattering in p-wave with I=1/2I=1/2 on the lattice. Using KπKπ systems with non-vanishing total momenta allows the extraction of phase shifts at several values of KπKπ relative momenta. A Breit-Wigner fit of the phase renders a K∗(892)K^*(892) resonance mass and K∗→KπK^*\to K π coupling compatible with the experimental numbers. We also determine the K∗(1410)K^*(1410) mass assuming the experimental K∗(1410)K^*(1410) width. We contrast the resonant I=1/2I=1/2 channel with the repulsive non-resonant I=3/2I=3/2 channel, where the phase is found to be negative and small, in agreement with experiment.

preprint2013arXiv

Search for Z_c^+(3900) in the 1^+- channel on the lattice

Recently three experiments reported a discovery of manifestly exotic Z_c^+(3900) in the decay to J/psi pi+, while J and P are experimentally unknown. We search for this state on the lattice by simulating the channel with J^PC=1^+- and I=1, and we do not find a candidate for Z_c^+(3900). Instead, we only find discrete scattering states D D* and J/psi pi, which inevitably have to be present in a dynamical QCD. The possible reasons for not finding Zc+Z_c^+ may be that its quantum numbers are not 1+−1^{+-} or that the employed interpolating fields are not diverse enough. Simulations with additional types of interpolators will be needed to reach a more definite conclusion.

preprint2012arXiv

K pi scattering for isospin 1/2 and 3/2 in lattice QCD

We simulate K pi scattering in s-wave and p-wave for both isospins I=1/2, 3/2 using quark-antiquark and meson-meson interpolating fields. We extract the elastic phase shifts delta at several values of the K-pi relative momenta. The resulting phases exhibit qualitative agreement with the experimental phases in all four channels. We express the s-wave phase shifts near threshold in terms of the scattering length and the effective range. Our K pi system has zero total momentum and is simulated on a single ensemble with two dynamical quarks, so results apply for mpi=266 MeV and mK=552 MeV in our simulation. The backtracking contractions in both I=1/2 channels are handled by the use of Laplacian-Heavyside smeared quarks within the distillation method. Elastic phases are extracted from the energy levels using Luscher's relations. In all four channels we observe the expected K(n)pi(-n) scattering states, which are shifted due to the interaction. In both attractive I=1/2 channels we observe additional states that are related to resonances; we attribute them to K_0^*(1430) in s-wave and K*(892), K*(1410) and K*(1680) in p-wave.

preprint2012arXiv

Scattering phase shifts for two particles of different mass and non-zero total momentum in lattice QCD

We derive the relation between the scattering phase shift and the two-particle energy in the finite box, which is relevant for extracting the strong phase shifts in lattice QCD. We consider elastic scattering of two particles with different mass and with non-zero total momentum in the lattice frame. This is a generalization of the Luscher formula, which considers zero total momentum, and the generalization of Rummukainen-Gottlieb's formula, which considers degenerate particles with non-zero total momentum. We focus on the most relevant total momenta in practice, i.e. P=(2π/L) e_z and P=(2π/L) (e_x+e_y) including their multiples and permutations. We find that the P-wave phase shift can be reliably extracted from the two-particle energy if the phase shifts for l>=2 can be neglected, and we present the corresponding relations. The reliable extraction of S-wave phase shift is much more challenging since delta(l=0) is always accompanied by delta(l=1) in the phase shift relations, and we propose strategies for estimating delta(l=0). We also propose the quark-antiquark and meson-meson interpolators that transform according the considered irreducible representations.

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.