\section{Introduction}
The determination of the compositeness, or molecular component of physical hadronic states has been the subject of multiple discussions starting from the pioneer work of Weinberg~\cite{PhysRev.137.B672}. In one of the most popular variants of the relationship of low energy scattering observables to the compositeness $X_{W}$, or probability to have a molecular state of two hadrons, one finds  
\begin{align}\label{1_1}
    a=R\bigg[\frac{2X_W}{1+X_W}+O(\frac{R_{\mathrm{typ}}}{R})\bigg]
\end{align}
\begin{align}\label{1_2}
    r_0=R\bigg[-\frac{1-X_W}{X_W}+O(\frac{R_{\mathrm{typ}}}{R})\bigg]
\end{align}
$R=1/\sqrt{2\mu B}$, $\mu=\frac{m_1m_2}{m1+m_2}$, $B=E_{\mathrm{th}}-E_0$, $E_0$ and $E_{\mathrm{th}}$ the energy of the bound state and the threshold energy of the two particles with mass $m_1$and $m_2$, with $a$, $r_0$ the scattering length and effective range in the approximation of the scattering matrix,
\begin{align}\label{1_3}
    f=\frac{1}{k\cot\delta-ik}\approx \frac{1}{-\frac{1}{a}+\frac{1}{2}r_0k^2-ik}
\end{align}

Eqs.~(\ref{1_1}),~(\ref{1_2}) hold in the limit of small binding energy $B$ and $O(\frac{R_{\mathrm{typ}}}{R})$ are corrections to the formula stemming from a typical scale $R_{\mathrm{typ}}$ related to the range of the interaction.

If we invert Eqs.~(\ref{1_1}),~(\ref{1_2}), we find 
\begin{align}\label{2_1}
    X_W=\frac{a}{2R-a}+O(\frac{R_{\mathrm{typ}}}{R})
\end{align}
\begin{align}\label{2_2}
    X_W=\frac{R}{R-r_0}+O(\frac{R_{\mathrm{typ}}}{R})
\end{align}
and one also defines $Z=1-X$ as the amount of non molecular component of the bound state.

If we apply these rules to the deuteron, were for $I=0$, $J=1$, we have
\begin{align}\label{2_3}
    a=5.419(7)\mathrm{fm},\qquad~r_0=1.766(8)\mathrm{fm},\qquad~B=2.224575(9)\mathrm{MeV}
\end{align}
one finds the well known surprise~\cite{Guo:2017jvc,Esposito:2021vhu,Li:2021cue} that 
\begin{align}\label{2_4}
    X=1.68\quad \mathrm{from\quad Eq.~(\ref{2_1})},\qquad~X=1.69 \quad \mathrm{from\quad Eq.~(\ref{2_2})}
\end{align}
with an unacceptable result since $X\leq1$. Eq.~(\ref{1_2}) also yields an unacceptable result, which forces $r_0$ to be negative if $X_W$ has an acceptable value around $1$, Yet, the test is accepted as a good guess for a molecular $np$ structure of the deuteron because the alternative is even worse: For $X=0.1$, Eq.~(\ref{1_2}) yields $r_0=-38.86\mathrm{fm}$, very far from the relative small number of Eq.~(\ref{2_3}). Indeed, this is the argument of Weinberg quoting ``The true token that the deuteron is composite is that $r_0$ is small and positive rather than large and negative''. Since for different reasons the molecular picture of the deuteron is acceptable, one can only conclude that the corrections $O(\frac{R_{\mathrm{typ}}}{R})$ in Eqs.~(\ref{1_1}),~(\ref{1_2}) are very large. However, they are difficult to quantity in terms of low energy variables.

Lots of energies have been devoted to improve the Weinberg conditions, Eqs.~(\ref{1_1}),~(\ref{1_2}), and understand the meaning of the compositeness in the presence of coupled channels or energy dependent potentials~\cite{BARU200453,Gamermann:2009uq,Hyodo:2011qc,Baru:2010ww,Hanhart:2011jz,Aceti:2012dd,Hyodo:2013iga,Sekihara:2014kya,Hanhart:2014ssa,Guo:2015daa,Sekihara:2015gvw,Kamiya:2015aea,Sekihara:2016xnq,Kamiya:2016oao,Matuschek:2020gqe,Guo:2017jvc}. Very recently the role of the range of the interaction has been revisited in~\cite{Kinugawa:2021ykv,Li:2021cue}. In~\cite{Kinugawa:2021ykv} a simple field theoretical approach is done containing a point like interaction and a derivative coupling term. A different version is also studied incorporating an additional field that couples to the two body scattering state. In these cases it is shown that the effective range $r_0$ originates from the derivative coupling interaction or from the channel coupling to the bare state. It is suggested that $R_{\mathrm{typ}}$ in Eqs.~(\ref{1_1}),~(\ref{1_2}) should be $R_{\mathrm{typ}}=\mathrm{max}[R_{\mathrm{int}},R_{\mathrm{eff}}]$, where $R_{\mathrm{int}}$ is $1/\Lambda$($\Lambda$, cutoff or range of the interaction) and $R_{\mathrm{eff}}$ is a length scale in the effective range expansion of Eq.~(~\ref{1_3}). The developments are formal and some examples are worked out in~\cite{Kinugawa:2021ybb}.

A different approach is followed in~\cite{Li:2021cue}. The Low equation~\cite{PhysRev.137.B672} for the scattering matrix is used and the non pole term is neglected. Some form factors are introduced and an expression for $X$ is obtained which depends on phase shifts obtained with some approximations, which has the virtue of being smaller than $1$. A discussion is done on how the results can depend on the form factors assumed.

The purpose of the present work is different. We share with the former works the idea that the range of the interaction is relevant in the determination of $X$. The aim, however, is to see how much information can we obtain from the combined knowledge of $a$, $r_0$ and the binding, not knowing exactly which is the range of the interaction. For this purpose we start with a formalism that takes into account the range of the interaction and incorporates a possible energy dependence of the potential. With the information on the binding, $a$ and $r_0$ we investigate the range of values of $X$ that the formalism provides. As we shall see, we can obtain a more accurate prediction for $X$ for the deuteron, and we can obtain a positive $r_0$, while at the same time we obtain a qualitative information on the range of the interaction. We also discuss two more states, the $D^{*}_{s0}(2317)$ and $D^{*}_{s1}(2460)$, which were studied in detail from a lattice QCD perspective in~\cite{MartinezTorres:2014kpc}, and the value of $X$, $a$, $r_0$ were determined. The lattice data of~\cite{MartinezTorres:2014kpc} contain more information than just $a$, $r_0$ and the binding, which allowed one to determine $X$ with relative precision. In the present work we shall exploit how much information we can get on $X$ from the knowledge of $a$, $r_0$ and $B$ alone.

\section{formalism}
\subsection{Formalism for the meson meson interaction}
Let us start with a meson meson interaction and the formalism employed in the studies of the chiral unitary approach. To put the problem in perspective let us start with a two channel problem, in which one channel is more important than the other and can lead to a bound state with this interaction. Let the interaction be given by the matrix 
\begin{align}\label{5_1}
    V=
    \left(
  \begin{array}{cc}
    V_{11} & V_{12} \\ 
    V_{12} & 0 \\
  \end{array}
\right),
\end{align}
where for simplicity we have made $V_{22}=0$. The formalism can be generalized to more coupled channels keeping all $V_{ij}$ terms~\cite{Hyodo:2013nka}. We assume $V_{ij}$ in Eq.~(\ref{5_1}) to be energy independent. The $T$ matrix is then given by ~\cite{Oller:1997ti} 
\begin{align}\label{5_2}
    T=[1-VG]^{-1}V
\end{align}
where $G=\mathrm{diag}(G_l)$ is the loop function of two mesons, which we regularize with the cut off method~\cite{Oller:1997ti}
\begin{align}\label{5_3}
   G_l=\int_{|q|< q_{\mathrm{max}}}\frac{d^3q}{(2\pi)^3}\frac{w_1(q)+w_2(q)}{2w_1(q)w_2(q)}\frac{1}{s-(w_1(q)+w_2(q))^2+i\epsilon}
\end{align}
where $w_i=\sqrt{\Vec{q}^2+m^2_i}$ and $s$ is the square of the incoming energy in the rest frame of the two particles {\footnote{In Ref.~\cite{Gamermann:2009uq} it is shown that Eqs.~(\ref{5_1}),~(\ref{5_2}),~(\ref{5_3}) can be obtained starting with a potential $V_{ij}(\textbf{q},\textbf{q}')=V_{ij}\theta({q_{\mathrm{max}}-|\bar{\vec{q}}|})\theta({q_{\mathrm{max}}-|\bar{\vec{q}}'|})$. Hence, one has to see $q_{\mathrm{max}}$ as the range of the interaction in momentum space.}}. 

Let us assume that Eq.~(\ref{5_2}) has a pole at $s_R$ corresponding to a bound state in channels $1$ and $2$. We define the couplings of the state to the channels $1$ and $2$ as
\begin{align*}
  g_1^2=\lim\limits_{s-s_R}(s-s_R)t_{11},\qquad g_2^2=\lim\limits_{s-s_R}(s-s_R)t_{22}
\end{align*}
and we can explicitly prove that (see section 5 of ~\cite{Aceti:2014ala})
\begin{align}\label{6_1}
    -g_1^2\frac{\partial G}{\partial s}|_{s=s_R}-g_2^2\frac{\partial G}{\partial s}|_{s=s_R}=P_1+P_2=1
\end{align}
where $P_1$, $P_2$ represent the probability to find the state in the channel $1$ and $2$. A general proof for many channels can be seen in ~\cite{Gamermann:2009uq,Hyodo:2013nka}. As shown in section 6 of \cite{Aceti:2014ala}, one can construct an effective potential for channel $1$
\begin{align}\label{6_2}
    V_{\mathrm{eff}}=V_{11}+V_{12}^2G_2
\end{align}
such that $t_{11}$ and $g_1^2$ are the same in channel $1$ as a single channel with $V_{\mathrm{eff}}$  than with the two body problem, and hence $-g_1^2\frac{\partial G_1}{\partial s}$ is still the probability to the find channel $1$ in the molecular state. Note that the price we payed to eliminate channel $2$ is that the effective potential in channel $1$, Eq.~(\ref{6_2}) is now energy dependent because $G_2$ depends on the energy. Now we have with one channel 
\begin{align}
    T_{11}=[1-V_{\mathrm{eff}}G_1]^{-1}V_{\mathrm{eff}}=\frac{1}{V_{\mathrm{eff}}^{-1}-G_1}
\end{align}
\begin{align*}
  g_1^2=\lim\limits_{s-s_R}(s-s_R)t_{11}=\frac{1}{\frac{\partial{V_{\mathrm{eff}}^{-1}}}{\partial s}|_{s_R}-\frac{\partial G}{\partial s}|_{s_R}}  
\end{align*}
and hence
\begin{align}\label{6_3}
    -g_1^2\frac{\partial G_1}{\partial s}|_{s_R}=P_1,\qquad
    g_1^2\frac{\partial{V_{\mathrm{eff}}^{-1}}}{\partial s}|_{s_R}=--g_1^2\frac{1}{V_{\mathrm{eff}}^{2}}\frac{\partial{V_{\mathrm{eff}}}}{\partial s}=1-P_1\equiv P_2\equiv Z
\end{align}
To have a feeling of the energy dependence of $V_{\mathrm{eff}}$, we particularize to a case that will be studied later, the $D^{*}_{s0}(2317)$, which qualifies as a molecular state of $DK$ and $\eta D_s$, mostly $DK$~\cite{Kolomeitsev:2003ac,Guo:2006fu,Gamermann:2006nm,Guo:2009ct,Wang:2012bu,Liu:2012zya,Altenbuchinger:2013gaa}. In Fig.~\ref{G_s} we show the $G(s)$ function for the $DK$ and $\eta D_s$ channels. We can see that around the energy of the $D^{*}_{s0}(2317)$, between $\sqrt{s}=2.236$~MeV and the $DK$ threshold, the $G_{\eta D_s}(s)$ function is very well represented by a linear function in $s$. 
\begin{figure}[H]
  \centering
 
  \includegraphics[width=0.45\textwidth]{G_fun.eps}
  \caption{The $G(s)$ function for the $DK$ and $\eta D_s$ channels at $q_{\mathrm{max}}=650$~MeV. The black solid line and red dashed line present the $DK$ and $\eta D_s$ channels respectively, with the black dotted line and red dashed dotted line corresponding to the threshold of the $KD$ and $\eta D_s$ channels.}\label{G_s}
\end{figure}
We shall then assume as a general rule that the $G$ function of possible missing channels, when studying a given one, can be approximated by a linear function of $s$ in the neighborhood of the pole of the studied channel, such that we can well approximate $V_{\mathrm{eff}}$ by
\begin{align}\label{7_1}
   V_{\mathrm{eff}}=V_{0}+\beta(s-s_{0})
\end{align}
with $\beta$ negative, as seen in Fig.~\ref{G_s}. If we eliminated some channel with a threshold below the threshold of the studied channel, $\beta$ would be positive. It is possible to tackle also this case, but it is convenient to consider explicitly the possible decay channels. Hence, we stick to Eq.~(\ref{7_1}) with $\beta\leq 0$.

The scattering matrix with potential $V_{\mathrm{eff}}$ of Eq.~(\ref{7_1}) is given for the one channel system that we wish to investigate by
\begin{align}\label{8_1}
    T(s)=\frac{1}{[V_{0}+\beta(s-s_{0})]^{-1} -G(s)}
\end{align}
and we impose that it has a pole at $s_0$, giving rise to a state whose nature we want to investigate. We will have 
\begin{align}\label{8_11}
    V_{0}^{-1} -G(s_0)=0,\qquad V_{0}=\frac{1}{G(s_0)}
\end{align}
Thus,
\begin{align}\label{8_2}                          T(s)    =&\frac{1}{[\frac{1}{G(s_0)}+\beta(s-s_0)]^{-1} -G(s)}
\end{align}
where we have eliminated the unknown $V_0$.

It is easy to establish the connection of this expression with the current amplitude of Quantum Mechanics by restoring the normalization of $T(s)$, governed by the normalization of the $G(s)$ function of Eq.~(\ref{5_3}). We have 
\begin{align}\label{8_3}
    \frac{1}{C\bigg\{[\frac{1}{G(s_0)}+\beta(s-s_0)]^{-1} -G(s)\bigg\}}\approx\frac{1}{-\frac{1}{a}+\frac{1}{2}r_0k^2-ik}
\end{align}
Since 
\begin{align}
   \mathrm{Im}G=&
   -\frac{1}{8\pi}\frac{k}{\sqrt{s}}
\end{align}
we can establish the connection
\begin{align}
    8\pi\sqrt{s}\bigg\{[\frac{1}{G(s_0)}+\beta(s-s_0)]^{-1} -G(s)\bigg\}+ik
        \approx &\frac{1}{a}-\frac{1}{2}r_0k^2+ik
\end{align}
or
\begin{align}\label{8_4}
    8\pi\sqrt{s}\bigg\{[\frac{1}{G(s_0)}+\beta(s-s_0)]^{-1} -\mathrm{Re}G(s)\bigg\}
        \approx &\frac{1}{a}-\frac{1}{2}r_0k^2
\end{align}

From Eq.~(\ref{8_4}), we obtain two equations evaluating the expression on the left at threshold and its derivative with respect to $k^2$ at threshold for positive increases of $s$. We have 
\begin{align}\label{9_1}
    8\pi\sqrt{s_\mathrm{th}}\bigg\{[\frac{1}{G(s_0)}+\beta(s_{\mathrm{th}}-s_0)]^{-1} -\mathrm{Re} G(s_{\mathrm{th}})\bigg\}
        =&\frac{1}{a}
\end{align}
and talking with account that
  $s=(w_1(k)+w_2(k))^2$, $w_1(k)=\sqrt{m_1^2+k^2}$, $w_2(k)=\sqrt{m_2^2+k^2}$,
\begin{align}
    \frac{\partial }{\partial k^2}=\frac{\partial }{\partial s}\frac{\partial s}{\partial k^2}=2(w_1(k)+w_2(k))(\frac{1}{2w_1(k)}+\frac{1}{2w_2(k)})\frac{\partial }{\partial s}=\frac{(w_1(k)+w_2(k))^2}{w_1(k)w_2(k)}\frac{\partial }{\partial s}.
\end{align}
we obtain
\begin{align}\label{9_2}
    &\frac{1}{2\sqrt{s_\mathrm{th}}}8\pi\bigg[[\frac{1}{G(s_0)}+\beta(s_\mathrm{th}-s_0)]^{-1} -\mathrm{Re}G(s)_\mathrm{th}\bigg]\frac{s}{w_1(k)w_2(k)}|_{s_{\mathrm{th}}}\\\nonumber
    &+
    8\pi\sqrt{s_\mathrm{th}}\bigg[-\beta[\frac{1}{G(s_0)}+\beta(s_\mathrm{th}-s_0)]^{-2}-\frac{\partial{\mathrm{Re}[G(s)]}}{\partial s}|_{s_{\mathrm{th}}}\bigg]\frac{s}{w_1(k)w_2(k)}|_{s_{\mathrm{th}}}\\\nonumber
    =&    
    -\frac{1}{2}r_0
\end{align}
From Eq.~(\ref{9_1}) we get the value of $\beta$ from $a$, as
\textbf{\begin{align}\label{9_3}
    \beta=&\frac{1}{s_{\mathrm{th}}-s_0}\bigg\{[\frac{1}{a}\frac{1}{8\pi}\frac{1}{\sqrt{s_{\mathrm{th}}}}+\mathrm{Re}G(s_{\mathrm{th}})]^{-1}-\frac{1}{G(s_0)}\bigg\}
\end{align}}
Note that using Eqs.~(\ref{6_3}),~(\ref{7_1}),~(\ref{8_11}) we have 
\begin{align}\label{9_4}
    P_2=1-P_1=Z=-g^2G(s_0)^2\beta
\end{align}
with 
\begin{align}\label{9_5}
        g^2&=\lim\limits_{s-s_0}(s-s_0)T(s)=\lim\limits_{s-s_0}\frac{s-s_0}{[\frac{1}{G(s_0)}+\beta(s-s_0)]^{-1} -G(s)}\\\nonumber
        &=\frac{1}{-(\frac{1}{G(s_0)})^{-2}\beta-\frac{\partial G}{\partial s}|_{s_0}}=\frac{1}{-G(s_0)^2\beta-\frac{\partial G}{\partial s}|_{s_0}}
\end{align}
where we have applied L'Hospital's rule to calculate the limit. Eq.~(\ref{9_4}) with $\beta \leq 0$ will always guarantee that $Z$ is positive, as it should be (note that $\frac{\partial G}{\partial s}|_{s_0}< 0$).

Eq.~(\ref{9_3}), with the knowledge of $a$ and $q_{\mathrm{max}}$ and the position of the pole will give the value of $\beta$  and by means of Eqs.~(\ref{9_4}),~(\ref{9_5}), we will get the value of $Z$, the complement of the compositeness of the state that we study. With the value of $\beta$ obtained we can go to Eq.~(\ref{9_2}) and establish that 
\begin{align}\label{10_1}
   R_0\equiv &-\frac{1}{2\sqrt{s_\mathrm{th}}}16\pi\bigg[[\frac{1}{G(s_0)}+\beta(s_\mathrm{th}-s_0)]^{-1} -\mathrm{Re}G(s_\mathrm{th})\bigg]\frac{s}{w_1(k)w_2(k)}|_{s_{\mathrm{th}}}\\\nonumber
    &+
    16\pi\sqrt{s_\mathrm{th}}\bigg[\beta[\frac{1}{G(s_0)}+\beta(s_\mathrm{th}-s_0)]^{-2}+\frac{\partial{\mathrm{Re}[G(s)]}}{\partial s}|_{s_{\mathrm{th}}}\bigg]\frac{s}{w_1(k)w_2(k)}|_{s_{\mathrm{th}}}\\\nonumber
    =&    
    r_0
\end{align}
Given $a$ and $s_0$, $R_0$ in Eq.~(\ref{10_1}) is now a function of $q_{\mathrm{max}}$ and we can see which value of $q_{\mathrm{max}}$ we need to satisfy the equation or how close or far we are from satisfying it.

Note that we have $\frac{\partial G}{\partial s}$ which we evaluate numerically. For this we evaluate Eq.~(\ref{5_3}) analytically, using the formula given in Ref.~\cite{Oller:1998hw}.

\subsection{Formalism for the nucleon nucleon interaction}
We follow the steps of the former subsection but we use the variable $E$ instead of $s$ ($E^2=s$). The $G$ function is now defined as 
\begin{align}\label{11_1}
   G(E)=\int_{|q|< q_{\mathrm{max}}}\frac{d^3q}{(2\pi)^3}\frac{m_1(q)m_2(q)}{E_1(q)E_2(q)}\frac{1}{\sqrt{s}-E_1(q)-E_2(q)+i\epsilon}
\end{align}
It corresponds to Eq.~(\ref{5_3}) multiplied by $2m_12m_2$ for reasons of normalization of the fields and neglecting the negative energy parts of the relativistic propagator, as appropriate for heavy particles as the nucleons. In practice, we use the same formula of Ref.~\cite{Oller:1998hw} multiplied by $4m_1m_2$.

The potential is now
\begin{align}
   V=V_{0}+\beta(E-E_0) 
\end{align}
Then,
\begin{align}\label{T_baryoncheck}
    T=\frac{1}{[V_{0}+\beta(E-E_0)]^{-1} -G(E)}
\end{align}
The pole at $E_0$ implies
\begin{align}\label{11_1}
    V_{0}^{-1} -G(E_0)=0,\qquad V_{0}=\frac{1}{G(E_0)}
\end{align}
hence,
\begin{align}\label{11_2}
    T       =&\frac{1}{[\frac{1}{G(E_0)}+\beta(E-E_0)]^{-1} -G(E)}\\\nonumber
\end{align}
but now
\begin{align}
   \mathrm{Im}~G=-\frac{1}{8\pi}k\frac{2m_12m_2}{\sqrt{s}}
   =-\frac{1}{2\pi}k\frac{m_1m_2}{E}
\end{align}
and equivalently to Eq.~(\ref{8_4}) we have now
\begin{align}\label{11_3}
    \frac{2\pi E}{m_1m_2}\bigg\{[\frac{1}{G(E_0)}+\beta(E-E_0)]^{-1} -\mathrm{Re}G(E)\bigg\}
        \approx &\frac{1}{a}-\frac{1}{2}r_0k^2
\end{align}
Evaluating Eq.~(\ref{11_3}) at threshold we get 
\begin{align}\label{12_1}
    \frac{2\pi E_{\mathrm{th}}}{m_1m_2}\bigg\{[\frac{1}{G(E_0)}+\beta(E-E_0)]^{-1} -\mathrm{Re}G(E_{\mathrm{th}})\bigg\}
        =&\frac{1}{a}
\end{align} 
and using $E=E_1(k)+E_2(k)$, $E_i(k)=\sqrt{m_i^2+k^2}$,
\begin{align}
    \frac{\partial }{\partial k^2}=\frac{\partial}{\partial E}\frac{\partial E}{\partial k^2}=(\frac{1}{2E_1(k)}+\frac{1}{2E_2(k)})\frac{\partial}{\partial E}=\frac{1}{2}\frac{E_1+E_2}{E_1(k)E_2(k)}\frac{\partial}{\partial E},
\end{align}
from the derivative of Eq.~(\ref{11_3}) at threshold we get
\begin{align}\label{12_2}
   &\frac{2\pi}{m_1m_2}\bigg[[\frac{1}{G(E_0)}+\beta(E_\mathrm{th}-E_0)]^{-1} -\mathrm{Re}G(E)_\mathrm{th}\bigg]\frac{E}{2E_1(k)E_2(k)}|_{E_{\mathrm{th}}}\\\nonumber
    &+
    \frac{2\pi E_{\mathrm{th}}}{m_1m_2}\bigg[-\beta[\frac{1}{G(E_0)}+\beta(E_\mathrm{th}-E_0)]^{-2}-\frac{\partial{\mathrm{Re}[G(E)]}}{\partial E}|_{E_{\mathrm{th}}}\bigg]\frac{E}{2E_1(k)E_2(k)}|_{E_{\mathrm{th}}}\\\nonumber
    =&    
    -\frac{1}{2}r_0
\end{align}
or
\begin{align}\label{12_2}
   R_0\equiv &-\frac{2\pi}{m_1m_2}\bigg[[\frac{1}{G(E_0)}+\beta(E_\mathrm{th}-E_0)]^{-1} -\mathrm{Re}G(E_\mathrm{th})\bigg]\frac{m_1+m_2}{m_1m_2}\\\nonumber
    &+
    \frac{2\pi (m_1+m_2}{m_1m_2})\bigg[\beta[\frac{1}{G(E_0)}+\beta(E_\mathrm{th}-E_0)]^{-2}+\frac{\partial{\mathrm{Re}[G(E)]}}{\partial E}|_{E_{\mathrm{th}}}\bigg]\frac{m_1+m_2}{m_1m_2}\\\nonumber
    =&
    r_0
\end{align}
Once again we will compare $R_0$ versus $r_0$ as a function of $q_{\mathrm{max}}$,. The couplings are now defined as  
\begin{align}
        g^2&=\lim\limits_{E-E_0}(E-E_0)T=\lim\limits_{E-E_0}\frac{E-E_0}{[\frac{1}{G(E_0)}+\beta(E-E_0)]^{-1} -G(E)}\\\nonumber
        &=\frac{1}{-[\frac{1}{G(E_0)}]^{-2}\beta-\frac{\partial G(E)}{\partial E}|_{E_0}}=\frac{1}{-G(E_0)^2\beta-\frac{\partial G(E)}{\partial E}|_{E_0}}
\end{align}
and 
\begin{align}\label{13_0}
  Z=1-X=-g^2G(E_0)^2\beta, \qquad   X=-g^2\frac{\partial G(E)}{\partial E}|_{E_0}
\end{align}

The magnitude $\beta$ is now evaluated from Eq.~(\ref{12_1})
\begin{align}\label{13_1}
    \beta=&\frac{1}{E_{\mathrm{th}}-E_0}\bigg\{[\frac{1}{a}\frac{1}{2\pi}\frac{m_1m_2}{m_1+m_2}+\mathrm{Re}G(E_{\mathrm{th}})]^{-1}-\frac{1}{G(E_0)}\bigg\}
\end{align}
which determines $\beta$ from $a$, $E_0$ and $q_{\mathrm{max}}$. Substituted in Eq.~(\ref{12_2}), it allows to check $R_0$ versus $r_0$ as a function of $q_{\mathrm{max}}$.

\section{results}
\subsection{The deuteron case}
We have the data of Eqs.~(\ref{2_3}). We determine $\beta$ in terms of this input and Eq.~(\ref{13_1}) and then the value of $Z$ of Eq.~(\ref{13_0}) as a function of $q_{\mathrm{max}}$ and plot $Z$ in Fig.~\ref{NNZ}
\begin{figure}[H]
  \centering
 
  \includegraphics[width=0.45\textwidth]{NN.eps}
  \caption{Z for the deuteron case as a function of $q_{\mathrm{max}}$.}\label{NNZ}
\end{figure}

We can see in Fig.~\ref{NNZ}, that starting from $q_{\mathrm{max}}=100$~MeV, $Z$ is a small number, smaller than $0.25$ indicating a strong molecular $pn$ component. If we go beyond $q_{\mathrm{max}}=140$~MeV then $Z$ becomes negative and we should discard this situation. The comparison of the theoretical $R_0$ versus the experimental $r_0$ value is shown in Fig.~\ref{deotron_r}. There we can see that $R_0$ is close to $r_0$ for values of $q_{\mathrm{max}}\geq140$~MeV, but the disagreement becomes noticeable below this value. Altogether, there is a optimal situation around $q_{\mathrm{max}}\approx140$~MeV, where $Z\sim0$, indicating that the deuteron is a molecular state and the value of $R_0$ is very close to $r_0$. This situation is realistic indicating that the range of the $NN$ interaction in the deuteron in $r$-space is rather large, and hence $q_{\mathrm{max}}$ is small. If we look at the deuteron wave function in momentum space from the Bonn potential~\cite{Machleidt:2000ge}, we see that at $q=140$~MeV, the wave function is $7$\%, the value at the origin and $\Phi^2(q=140~\mathrm{MeV})\approx0.006\Phi^2(0)$. We find then that the optimal agreement of the theory with the $a$ and $r_0$ data is obtained with a value of $q_{\mathrm{max}}$ which reflects realistically the actual deuteron wave function and tells us that $Z\sim0$, hence the deuteron is mostly a $pn$ molecule. There is no point demanding more than this information knowing that apart from the $S-$wave potential, there is also an important tensor part in the deuteron and a non negligible $D$-wave part in the wave function.
As we can see, the range of the interaction has been essential to obtain this acceptable picture. From Fig.~\ref{NNZ} we can see that if $q_{\mathrm{max}}$ increases, indicating short range interaction in $r$-space, $Z$ becomes negative and $X\approx1.6$ similar to what one obtains with the Weinberg formula of Eq.~(\ref{2_4}), not surprising since one implicitly is making  this assumption in the derivation of Weinberg formulas of Eqs.~(\ref{1_1}),~(\ref{1_2}) (see Ref.~\cite{Kinugawa:2021ykv}). Note that the consideration of the range of the interaction has also allowed us to obtain positive values of $R_0$ and close to the experimental one of $r_0$. Note also that in Eq.~(\ref{12_2}) the first term of $R_0$ is small and most of the contribution comes from the second term. If $\frac{\partial{\mathrm{Re}[G(E)]}}{\partial E}|_{E_{\mathrm{th}}}$ were zero then $R_0$ would be negative and proportional to $Z$, as found in Eq.~(\ref{1_2}). Actually this is the situation with non relativistic kinematics when the range $q_{\mathrm{max}}$ is set to infinity (see non relativistic dimensional regularization results for $G$ in Ref~\cite{Kaplan:1996xu} where $\mathrm{Re}G=\mathrm{const}$ above threshold). In summary, the consideration of the range of the interaction has rendered us a picture of the deuteron far closer to the actual molecular nature than Weinberg`s equations, and in return has shown that the interaction has to be of long range in $r$-space, with a realistic value of the range when compared to the actual deuteron wave function.
\begin{figure}[H]
  \centering
 
  \includegraphics[width=0.45\textwidth]{deotron_r.eps}
  \caption{Comparison of $R_0$ and $r_0$ for the deuteron case as a function of $q_{\mathrm{max}}$.}\label{deotron_r}
\end{figure}
\section{The $D^*_{s0}(2317)$}
Next we pay attention to the $D^*_{s0}(2317)$ state considered as a $KD$ molecule. We will take the values of $a$, $r_0$ and the binding from the QCD lattice analysis of the finite volume levels of~\cite{MartinezTorres:2014kpc} (note $\frac{1}{a}$ in~\cite{MartinezTorres:2014kpc} versus $-\frac{1}{a}$ here in Eq.~(\ref{1_3}))
\begin{align}
    a(KD)&=+1.3\pm0.5\pm0.1~\mathrm{fm}\\\nonumber
    r_0(KD)&=-0.1\pm0.3\pm0.1~\mathrm{fm}    
\end{align}
and we take the nominal mass $2317$~MeV for the mass of the state.

In Fig.~\ref{Z_2317} we see that $Z$ takes unrealistically large values for $q_{\mathrm{max}}< 300$~MeV. On the other hand, in Fig.~\ref{2317_r} we plot $R_0$ versus $r_0$ and we see that for $q_{\mathrm{max}}< 300$~MeV. The deviation of $R_0$ and $r_0$ becomes gradually large and unacceptable. We could say that for values of $q_{\mathrm{max}}> 400$~MeV we already obtain an acceptable agreement of $R_0$ versus $r_0$. We see then from Fig.~\ref{Z_2317} that in this case $Z< 0.4$, indicating a $DK$ molecular component with probability larger than $60$\%. This would be in agreement with the findings in~\cite{MartinezTorres:2014kpc} where it was found that 
\begin{align}
    P(DK)=(72\pm13\pm5)\%
\end{align}
\begin{figure}[H]
  \centering
 
  \includegraphics[width=0.45\textwidth]{Z_2317.eps}
  \caption{The value of $Z$ for $KD$ scattering forming the $D^{*}_{s0}(2317)$.}\label{Z_2317}
\end{figure}
\begin{figure}[H]
  \centering
 
  \includegraphics[width=0.45\textwidth]{R_0_2317.eps}
  \caption{The comparison of $R_0$ and $r_0$ for the $D^{*}_{s0}(2317)$ case as a function of $q_{\mathrm{max}}$.}\label{2317_r}
\end{figure}


We also see that the range, with $q_{\mathrm{max}}\ge 400$~MeV, corresponds to a shorter range in $r$-space then in the deuteron case. In the dynamical picture of~\cite{Gamermann:2006nm} the interaction is driven by vector meson exchange, contrary to the $np$ interaction where pion exchange plays a dominant role. 

Note also from Fig.~\ref{Z_2317} that around $q_{\mathrm{max}}= 725$~MeV $Z$ becomes zero and negative from there on. This also gives us an idea of the range demanded by the data, which reflects realistically what one can expect from light vector exchange. From comparison recall that $q_{\mathrm{max}}=630$~MeV was used in the study of the $\Bar{K}N$ interaction in~\cite{Oset:1997it}.

\section{The $D^*_{s1}(2460)$}
Now we analyse the $D^*_{s1}(2460)$ state shown to be mostly molecular in the $KD^{*}$ channel from the analysis of the lattice QCD  levels in~\cite{MartinezTorres:2014kpc}. There it was found
\begin{align}
    a(KD^*)&=+1.1\pm0.5\pm0.2~\mathrm{fm}\\\nonumber
    r_0(KD^*)&=-0.2\pm0.3\pm0.1~\mathrm{fm}    
\end{align} 
and we take the nominal mass $2460$~MeV for the state. In Fig.~\ref{Z_2460} we show the value of $Z$ as a function of $q_{\mathrm{max}}$. We have to go to values of $q_{\mathrm{max}}$ bigger than $280$~MeV to have $Z$ smaller than $1$. We can complement this information from Fig.~\ref{2460_r} which shows a big discrepancy of $R_0$ with $r_0$ for values of $q_{\mathrm{max}}$ smaller than $400$~MeV.

We can see now that $Z$ never becomes zero, independent of $q_{\mathrm{max}}$, reaching a value of $0.2$ for large $q_{\mathrm{max}}$. If we take a range of $q_{\mathrm{max}}$ like in the former case $400~\mathrm{MeV}<q_{\mathrm{max}}<750~\mathrm{MeV}$, the $Z$ range becomes $0.3< Z< 0.6$. What we can see is that $Z$ is now bigger than in the case of the $D^*_{s0}(2317)$. The compositeness is now around or bigger than $40$\% ($0.4< X< 0.7$) in agreement with the findings of~\cite{MartinezTorres:2014kpc}, where one obtains 
\begin{align}
    P(KD^*)=(57\pm21\pm6)\%
\end{align}
In this case the $\eta D^*_{s}$ channel is mostly responsible for the remaining probability. It is interesting to see that the analysis done here renders values of $Z$ which are in good agreement with those found in~\cite{MartinezTorres:2014kpc}.

What we have found in the three cases studied is that the range of the interaction is very important to consider, even in cases little bound like the deuteron. We could see that the combined information of $a$, $r_0$ and the binding could provide a fair information on the $D^*_{s0}(2317)$ and $D^*_{s1}(2460)$ which are bound by about $40-45$~MeV. At the same time the analysis gives us some idea about the range of the interaction, with the $NN$ interaction for the deuteron being of longer range in $r$-space than the $KD$ and $KD^*$ in the cases of the $D^*_{s0}(2317)$ and $D^*_{s1}(2460)$ states. The information obtained is more accurate and richer than that obtained by the limiting equations of (\ref{1_1}), (\ref{1_2}). The algorithms resulting from our study are also easy to implement and can be used to study other cases. 

\section{conclusion}
We have used a formalism for the scattering of hadrons following the steps used in the chiral unitary approach. When studying the scattering of the two hadrons sometimes a bound state appears below threshold and it is logical to ask oneself whether that state corresponds to a molecular state of these hadronic components, stemming from this interaction, or it corresponds to a different structure. Sometimes the state could have a mixture of another pair of hadrons, or even have a component of a compact quark cluster. To account for all these cases we have assumed an energy dependent potential, and our formalism allows one to determine the scattering length and the variable $q_{\mathrm{max}}$, the regulator in the loop functions, which gives the range of the interaction in momentum space. The formalism also allows one to obtain the effective range $r_0$, and not always is it possible to get agreement with experiment, indicating that more information beyond the effective range expansion would be necessary for a better analysis of the data. Yet, the comparison of the theoretical value $R_0$ and $r_0$  shows that there are region of  $q_{\mathrm{max}}$ where the disagreement is too big and unacceptable. Looking at the regions of  $q_{\mathrm{max}}$ where $R_0$ is closer to $r_0$ one can find a double information, which is the range of the interaction and the value of the $Z$, or the molecular compositeness $X=1-Z$. The combined analysis using the information of the binding, $a$ and $r_0$ renders us reasonable values of $Z$ and the range of the interaction for the three systems studied, with very different binding energy: the deuteron with $2.22$~MeV binding and the $D^*_{s0}(2317)$ and $D^*_{s1}(2460)$ states with binding around $40-45$~MeV. The information obtained with this method is more accurate and richer than that obtained from the formalism of Weinberg, derived in the limit of very small binding and zero range of the interaction in $r$-space.

\begin{figure}[H]
  \centering
 
  \includegraphics[width=0.45\textwidth]{Z_2460.eps}
  \caption{Z for the $D^{*}_{s1}(2460)$ case as a function of $q_{\mathrm{max}}$.}\label{Z_2460}
\end{figure}

\begin{figure}[H]
  \centering
 
  \includegraphics[width=0.45\textwidth]{R_0_2460.eps}
  \caption{The comparison of $R_0$ and $r_0$ for the $D^{*}_{s1}(2460)$ case as a function of $q_{\mathrm{max}}$.}\label{2460_r}
\end{figure}


\section{Acknowledgements}
We would like to thank Juan Nieves and Miguel Albaladejo for useful discussions. This work is partly supported by the National Natural Science Foundation of China under Grants Nos. 11975009, 12175066 and 12147219. This work is also partly supported by the Spanish Ministerio de Economia y Competitividad (MINECO) and European FEDER funds under Contracts No. FIS2017-84038-C2-1-PB, PID2020-112777GB-I00, and by Generalitat Valenciana under contract PROMETEO/2020/023. This project has received funding from the European Union Horizon 2020 research and innovation programme under the program H2020-INFRAIA-2018-1, grant agreement No. 824093 of the STRONG-2020 project. One of us, Jing Song wishes to acknowledge support from China Scholarship Council.






















































































































































































