Source author record

Luka Leskovec

Luka Leskovec 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

21works
4topics
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

21 published item(s)

preprint2023arXiv

Doubly charmed tetraquark: isospin channels and diquark-antidiquark interpolators

We perform a lattice simulation to investigate the doubly charmed tetraquark $T^+_{cc}$ observed by the LHCb collaboration, slightly below the $D^{*+}D^0$ threshold, with flavor content $cc\bar{u}\bar{d}$ and isospin-$0$. Two-meson interpolators are implemented to explore the isospin quantum numbers $I=0$ and $I=1$. We observe attraction near the $DD^*$ threshold for $I=0$ and repulsion for $I=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 $T^+_{cc}$. Finally, with a higher mass (close to the bottom quark), the ground state is shifted down significantly. The simulation is performed on $N_f=2+1$ CLS ensembles with $m_π\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 $P$-wave, $I=1$ two-pion final state, as a function of the photon virtuality and $ππ$ invariant mass. The calculation was performed with $2+1$ flavors of clover fermions at a pion mass of approximately $320$ MeV, on a $32^3 \times 96$ lattice with $L\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 $|\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 $q^2$ and $s$ 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)$.

preprint2021arXiv

$I=1/2$ $S$-wave and $P$-wave $Kπ$ scattering and the $κ$ and $K^*$ resonances from lattice QCD

We present a lattice-QCD determination of the elastic isospin-$1/2$ $S$-wave and $P$-wave $Kπ$ scattering amplitudes as a function of the center-of-mass energy using Lüscher's method. We perform global fits of $K$-matrix parametrizations to the finite-volume energy spectra for all irreducible representations with total momenta up to $\sqrt{3}\frac{2π}{L}$; this includes irreps that mix the $S$- and $P$-waves. Several different parametrizations for the energy dependence of the $K$-matrix are considered. We also determine the positions of the nearest poles in the scattering amplitudes, which correspond to the broad $κ$ resonance in the $S$-wave and the narrow $K^*(892)$ resonance in the $P$-wave. Our calculations are performed with $2+1$ dynamical clover fermions for two different pion masses of $317.2(2.2)$ and $175.9(1.8)$ MeV. Our preferred $S$-wave parametrization is based on a conformal map and includes an Adler zero; for the $P$-wave we use a standard pole parametrization including Blatt-Weisskopf barrier factors. The $S$-wave $κ$-resonance pole positions are found to be $\left[0.86(12) - 0.309(50)\,i\right]\:{\rm GeV}$ at the heavier pion mass and $\left[0.499(55)- 0.379(66)\,i\right]\:{\rm GeV}$ at the lighter pion mass. The $P$-wave $K^*$-resonance pole positions are found to be $\left[ 0.8951(64) - 0.00250(21)\,i \right]\:{\rm GeV}$ at the heavier pion mass and $\left[0.8718(82) - 0.0130(11)\,i\right]\:{\rm GeV}$ at the lighter pion mass, which corresponds to couplings of $g_{K^* Kπ}=5.02(26)$ and $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 $ \bar{b}\bar{b}ud $, $ \bar{b}\bar{b}us $ and $ \bar{b}\bar{c}ud $ using lattice QCD. We carry out calculations on several gauge link ensembles with $ 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\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+1$ flavors of clover fermions at a pion mass of approximately $320$ MeV and a lattice size of approximately $3.6$ fm. The required $2$-point and $3$-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\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 $ud\bar b \bar b$ with bottom quarks of finite mass which are treated in the framework of NRQCD. We focus on $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 $u \bar d b \bar b$ four-quark system with quantum numbers $I(J^P)=1(1^+)$ again using static bottom quarks.

preprint2016arXiv

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_π$ = 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ü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 $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.

preprint2015arXiv

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

Recently experimentalists have discovered several charged charmonium-like hadrons $Z_c^+$ with unconventional quark content $\bar cc\bar d u$. We perform a search for $Z_c^+$ with mass below $4.2~$GeV in the channel $I^G(J^{PC})=1^+(1^{+-})$ using lattice QCD. The major challenge is presented by the two-meson states $J/ψ\, π$, $ψ_{2S}π$, $ψ_{1D}π$, $D\bar D^*$, $D^*\bar D^*$, $η_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 $Z_c^+$ below $4.2~$GeV. Possible reasons for not seeing an additional eigenstate related to $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 $Z_c^+$ candidate, which is however not robust after further two-meson states around $4.2~$GeV are implemented.

preprint2015arXiv

Vector and scalar charmonium resonances with lattice QCD

We perform an exploratory lattice QCD simulation of $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 $D \bar D$ 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 $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)$ 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

$D_{s}$ mesons with $DK$ and $D^{*}K$ scattering near threshold

$D_s$ mesons are studied in three quantum channels ($J^P=0^+$, $1^+$ and $2^+$), where experiments have identified the very narrow $D_{s0}^*(2317)$, $D_{s1}(2460)$ and narrow $D_{s1}(2536)$, $D_{s2}^*(2573)$. We explore the effect of nearby $DK$ and $D^*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 $16^3\times 32$ or $32^3\times 64$). In addition to $\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_π\simeq 156~$MeV renders the $D_{s1}(2460)$ as a strong interaction bound state 44(10)MeV below $D^*K$ threshold, which is in agreement with the experiment. The $D_{s0}^*(2317)$ is found 37(17)MeV below $DK$ threshold, close to experiment value of 45MeV. The narrow resonances $D_{s1}(2536)$ and $D_{s2}^*(2573)$ are also found close to the experimental masses.

preprint2014arXiv

$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_π\approx 156$ MeV) a stochastic distillation variant is employed and for another ensemble ($16^3\times 32$, $m_π\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.

preprint2014arXiv

Axial resonances a1(1260), b1(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 $ρπ$ and $ωπ$. Interpolating fields $\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_π\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 $Γ_{a1}(s)=g^2 p/s$ is parametrized in terms of the coupling and we obtain $g_{a_1ρπ}=1.71(39)$ GeV compared to $g_{a_1ρπ}^{exp}=1.35(30)$ GeV derived from $Γ_{a1}^{exp}=425(175)$ MeV. In the $b_1$ channel, we find energy levels related to $π(0)ω(0)$ and $b_1(1235)$, and the lowest level is found at $E_1 \gtrsim m_ω+m_π$ but is within uncertainty also compatible with an attractive interaction. Assuming the coupling $g_{b_1ωπ}$ extracted from the experimental width we estimate $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 $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/ψπ$, $ψ_{2S}π$, $ψ_{1D}π$, $D \bar{D}^*$, $D^* \bar{D}^*$, $η_c ρ$ 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$.

preprint2013arXiv

$D_{s0}^*(2317)$ Meson and $D$-Meson-Kaon Scattering from Lattice QCD

The scalar meson $D_{s0}^*(2317)$ is found 37(17)MeV below DK threshold in a lattice simulation of the $J^P=0^+$ channel using, for the first time, both DK as well as $\bar sc$ interpolating fields. The simulation is done on $N_f=2+1$ gauge configurations with $m_π\simeq 156 $MeV, and the resulting $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 $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/ψω$ 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/ψπ$ 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/ψω$ 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π$ system to moving frames and present an exploratory extraction of the masses and widths for the $K^*$ resonances by simulating $Kπ$ scattering in p-wave with $I=1/2$ on the lattice. Using $Kπ$ systems with non-vanishing total momenta allows the extraction of phase shifts at several values of $Kπ$ relative momenta. A Breit-Wigner fit of the phase renders a $K^*(892)$ resonance mass and $K^*\to K π$ coupling compatible with the experimental numbers. We also determine the $K^*(1410)$ mass assuming the experimental $K^*(1410)$ width. 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.

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 $Z_c^+$ may be that its quantum numbers are not $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.