Binding of the three-hadron <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mi>D</mml:mi><mml:msup><mml:mi>D</mml:mi><mml:mo>*</mml:mo></mml:msup><mml:mi>K</mml:mi></mml:math> system from the lattice effective field the

We employ the nuclear lattice effective field theory (NLEFT), an efficient tool for nuclear calculations, to solve the asymmetric multihadron systems. We take the <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mi>D</a:mi><a:msup><a:mi>D</a:mi><a:mo>*</a:mo></a:msup><a:mi>K</a:mi></a:math> three-body system as an illustration to demonstrate the capability of the method. Here the two-body chiral interactions between <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"><c:mi>D</c:mi></c:math>, <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"><e:msup><e:mi>D</e:mi><e:mo>*</e:mo></e:msup></e:math>, and <g:math xmlns:g="http://www.w3.org/1998/Math/MathML" display="inline"><g:mi>K</g:mi></g:math> are regulated with a soft lattice regulator and calibrated with the binding energies of the <i:math xmlns:i="http://www.w3.org/1998/Math/MathML" display="inline"><i:msubsup><i:mi>T</i:mi><i:mrow><i:mi>c</i:mi><i:mi>c</i:mi></i:mrow><i:mo>+</i:mo></i:msubsup></i:math>, <k:math xmlns:k="http://www.w3.org/1998/Math/MathML" display="inline"><k:msubsup><k:mi>D</k:mi><k:mrow><k:mi>s</k:mi><k:mn>0</k:mn></k:mrow><k:mo>*</k:mo></k:msubsup><k:mo stretchy="false">(</k:mo><k:mn>2317</k:mn><k:mo stretchy="false">)</k:mo></k:math>, and <o:math xmlns:o="http://www.w3.org/1998/Math/MathML" display="inline"><o:msub><o:mi>D</o:mi><o:mrow><o:mi>s</o:mi><o:mn>1</o:mn></o:mrow></o:msub><o:mo stretchy="false">(</o:mo><o:mn>2460</o:mn><o:mo stretchy="false">)</o:mo></o:math> molecular states. We then calculate the three-body binding energy using the NLEFT and analyze the systematic uncertainties due to the finite volume effects, the sliding cutoff, and the leading-order three-body forces. Even when the three-body interaction is repulsive (even as large as the infinite repulsive interaction), the three-body system has a bound state unambiguously with binding energy no larger than the <s:math xmlns:s="http://www.w3.org/1998/Math/MathML" display="inline"><s:msub><s:mi>D</s:mi><s:mrow><s:mi>s</s:mi><s:mn>1</s:mn></s:mrow></s:msub><s:mo stretchy="false">(</s:mo><s:mn>2460</s:mn><s:mo stretchy="false">)</s:mo><s:mi>D</s:mi></s:math> threshold. To check the renormalization group invariance of our framework, we extract the first excited state. We find that when the ground state is fixed, the first excited states with various cutoffs coincide with each other when the cubic size goes larger. In addition, the standard angular momentum and parity projection technique is implemented for the quantum numbers of the ground and excited states. We find that both of them are <w:math xmlns:w="http://www.w3.org/1998/Math/MathML" display="inline"><w:mi>S</w:mi></w:math>-wave states with quantum number <y:math xmlns:y="http://www.w3.org/1998/Math/MathML" display="inline"><y:msup><y:mi>J</y:mi><y:mi>P</y:mi></y:msup><y:mo>=</y:mo><y:msup><y:mn>1</y:mn><y:mo>−</y:mo></y:msup></y:math>. Because the three-body state contains two charm quarks, it is easier to be detected in the Large Hadron Collider.