Position-space sampling for local multiquark operators in lattice QCD using distillation and the importance of tetraquark operators for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:msub> <mml:mi>T</mml:mi> <mml:mrow

Obtaining hadronic two-point functions is a central step in spectroscopy calculations in lattice QCD. This requires solving the Dirac equation repeatedly, which is computationally demanding. The distillation method addresses this difficulty by using the lowest eigenvectors of the spatial Laplacian to construct a subspace in which the Dirac operator can be fully inverted. This approach is efficient for nonlocal operators such as meson-meson and baryon-baryon operators. However, local multiquark operators with four or more (anti)quarks are computationally expensive in this framework: the cost of contractions scales with a high power of the number of Laplacian eigenvectors. We present a position-space sampling method within distillation that reduces this cost scaling by performing the momentum projection only over sparse grids rather than the full spatial lattice. We demonstrate the efficiency of this unbiased estimator for single-meson, single-baryon and local tetraquark operators. Using Wilson-clover fermions at the <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"> <a:mi>S</a:mi> <a:mi>U</a:mi> <a:mo stretchy="false">(</a:mo> <a:mn>3</a:mn> <a:mo stretchy="false">)</a:mo> </a:math> -flavor-symmetric point, we apply this method to study the importance of local tetraquark operators for extracting the finite-volume <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"> <e:msub> <e:mi>T</e:mi> <e:mrow> <e:mi>c</e:mi> <e:mi>c</e:mi> </e:mrow> </e:msub> <e:mo stretchy="false">(</e:mo> <e:mn>3875</e:mn> <e:msup> <e:mo stretchy="false">)</e:mo> <e:mo>+</e:mo> </e:msup> </e:math> spectrum. To this end, we extend a large basis of bilocal <i:math xmlns:i="http://www.w3.org/1998/Math/MathML" display="inline"> <i:mi>D</i:mi> <i:msup> <i:mi>D</i:mi> <i:mo>*</i:mo> </i:msup> </i:math> and <k:math xmlns:k="http://www.w3.org/1998/Math/MathML" display="inline"> <k:msup> <k:mi>D</k:mi> <k:mo>*</k:mo> </k:msup> <k:msup> <k:mi>D</k:mi> <k:mo>*</k:mo> </k:msup> </k:math> scattering operators by including local tetraquark operators. The inclusion of local operators leads to significant shifts in the estimates of several energy levels. Finally, we show the effect of these shifts on the <m:math xmlns:m="http://www.w3.org/1998/Math/MathML" display="inline"> <m:mi>D</m:mi> <m:msup> <m:mi>D</m:mi> <m:mo>*</m:mo> </m:msup> </m:math> scattering phase shift from a single-channel <o:math xmlns:o="http://www.w3.org/1998/Math/MathML" display="inline"> <o:mi>s</o:mi> </o:math> -wave Lüscher analysis.