\right) \,.
\end{split}$$ This expression can be written as $I=\sum_r^d f(a,b,c,d,r)$, where $f$ is a well-behaved function. Furthermore, we can write this result as $$I_{abcd}=\sum_{r=0}^d \exp{\log{[f(a,b,c,d,r)]}}\,,$$ where $f$ must be positive.\
\
The logarithm $\log{(f)}$ is given by $$\begin{split}
    \log{(f)}=&-\frac{1}{2}\log{(2)}-2\log{(\pi)}\\ &+\frac{1}{2}\left(\log{(c!)}+\log{(d!)}
    -\log{(a!)}-\log{(b!)}\right)\\
    &-\log{\left((d-r)!\right)}-\log{\left((c-d+r)!\right)}-\log{(r!)}\\
    &+\log{\left(\Gamma\left(\frac{a+b-c+d+1}{2}-r\right)\right)}\\
    &+\log{\left(\Gamma\left(\frac{a-b+c-d+1}{2}+r\right)\right)}\\
    &+\log{\left(\Gamma\left(\frac{-a+b+c-d+1}{2}+r\right)\right)}\,.
\end{split}$$ Notice that due to the presence of the gamma functions, the previous expression could require the evaluation of a logarithm of a negative quantity. For this reason, we calculate instead $\log{(\abs{\Gamma(n)})}$, and compute the associated phase separately. This phase is introduced in the final expression after the exponentiation. Thus, the final expression of the integral taking into account the possible negative values of gamma functions reads $$I_{abcd}=\sum_{r=0}^d p\exp{\log{[f(a,b,c,d,r)]}}\,,$$ where $p$ is the phase, given by $p=\prod_{n=1}^3 p_n$, with $p_n$ being the phase generated by each gamma function.

![Energy spectrum of the two-particle system as a function of the interaction strength. The calculations are performed using 100 single-particle modes (solid green line). The analytical spectrum (
\right) \, . 
 \end{split}$$ This expression can be written as $ I=\sum_r^d f(a, b, c, d, r)$, where $ f$ is a well - behaved function. Furthermore, we can spell this resultant role as $ $ I_{abcd}=\sum_{r=0}^d \exp{\log{[f(a, b, c, d, r)]}}\,,$$ where $ f$ must be positive.\ 
 \ 
 The logarithm $ \log{(f)}$ is given by $ $ \begin{split } 
     \log{(f)}=&-\frac{1}{2}\log{(2)}-2\log{(\pi)}\\ & + \frac{1}{2}\left(\log{(c!)}+\log{(d !) } 
     -\log{(a!)}-\log{(b!)}\right)\\ 
     & -\log{\left((d - r)!\right)}-\log{\left((c - d+r)!\right)}-\log{(r!)}\\ 
     & + \log{\left(\Gamma\left(\frac{a+b - c+d+1}{2}-r\right)\right)}\\ 
     & + \log{\left(\Gamma\left(\frac{a - b+c - d+1}{2}+r\right)\right)}\\ 
     & + \log{\left(\Gamma\left(\frac{-a+b+c - d+1}{2}+r\right)\right)}\, . 
 \end{split}$$ Notice that due to the presence of the gamma function, the previous expression could want the evaluation of a logarithm of a minus quantity. For this rationality, we calculate instead $ \log{(\abs{\Gamma(n)})}$, and calculate the associated phase separately. This phase is introduced in the final formula after the exponentiation. Thus, the final expression of the integral taking into history the possible negative values of da gamma functions reads $ $ I_{abcd}=\sum_{r=0}^d p\exp{\log{[f(a, b, c, d, r)]}}\,,$$ where $ p$ is the phase, give by $ p=\prod_{n=1}^3 p_n$, with $ p_n$ being the phase render by each gamma routine. 

! [ Energy spectrum of the two - particle system as a affair of the interaction strength. The calculations are performed using 100 single - particle modes (solid green line). The analytical spectrum (
\rigjt) \,.
\end{split}$$ This expresslon can be writtgn as $I=\snm_r^d f(a,g,c,d,r)$, whefe $f$ is a well-behaved functiln. Furtyermore, we can write tfis resuln as $$I_{abce}=\sum_{c=0}^d \exp{\log{[f(a,b,c,d,r)]}}\,,$$ xgere $f$ must bs posntmve.\
\
The logarithk $\log{(f)}$ is civen by $$\begin{vpuic}
    \log{(f)}=&-\frac{1}{2}\log{(2)}-2\log{(\pi)}\\ &+\frac{1}{2}\left(\log{(c!)}+\lod{(d!)}
    -\lpg{(w!)}-\log{(b!)}\right)\\
    &-\jog{\ltft((q-r)!\rifht)}-\log{\left((c-d+r)!\right)}-\log{(r!)}\\
    &+\log{\left(\Gzmma\lefu(\frac{a+b-c+d+1}{2}-r\right)\rignt)}\\
    &+\log{\left(\Gamma\left(\frac{w-b+c-d+1}{2}+g\right)\right)}\\
    &+\log{\peft(\Gamma\ledt(\frwx{-a+b+c-d+1}{2}+r\right)\rkght)}\,.
\end{splpc}$$ Notice thzt due to the presence of the gxmma yunctions, tye prfeious exprewsion could requivv the eealuatipn of a logaribhm oh a begative quantity. For this reason, we cajculate itscead $\log{(\abs{\Gamma(n)})}$, and cimputg the assuxiaged piass sepagatxly. This phzse is intriduced in the final evitession afted the qx[onentiation. Thus, the final expression mf fhe integral taking inti account the possiblg negative values of gamma functions reads $$I_{abcd}=\sum_{r=0}^d p\exp{\lmg{[f(a,b,r,d,f)]}}\,,$$ wkcve $p$ us the phase, given by $p=\prod_{n=1}^3 p_n$, with $p_n$ being trs khaxe generated bn each gamma functooj.

![Emgrgy spectrum uf the twk-particle system ad a funstion of the igtersction strength. The calcularions are pegforned using 100 single-pcrticle moder (sokid gteen line). The analyticau spsctrum (
\right) \,. \end{split}$$ This expression can be $I=\sum_r^d where $f$ a well-behaved function. result $$I_{abcd}=\sum_{r=0}^d \exp{\log{[f(a,b,c,d,r)]}}\,,$$ where must be positive.\ The logarithm $\log{(f)}$ is given by \log{(f)}=&-\frac{1}{2}\log{(2)}-2\log{(\pi)}\\ &+\frac{1}{2}\left(\log{(c!)}+\log{(d!)} -\log{(a!)}-\log{(b!)}\right)\\ &-\log{\left((d-r)!\right)}-\log{\left((c-d+r)!\right)}-\log{(r!)}\\ &+\log{\left(\Gamma\left(\frac{a+b-c+d+1}{2}-r\right)\right)}\\ &+\log{\left(\Gamma\left(\frac{a-b+c-d+1}{2}+r\right)\right)}\\ &+\log{\left(\Gamma\left(\frac{-a+b+c-d+1}{2}+r\right)\right)}\,. \end{split}$$ Notice that due to the of the gamma functions, the previous expression could require the evaluation of a of negative For reason, we calculate instead $\log{(\abs{\Gamma(n)})}$, and compute the associated phase separately. This phase is introduced in final expression after the exponentiation. Thus, the final of the integral taking account the possible negative values gamma reads $$I_{abcd}=\sum_{r=0}^d where is phase, given by p_n$, with $p_n$ being the phase generated by each gamma function. ![Energy spectrum of the two-particle system a function interaction strength. calculations performed 100 single-particle modes line). The analytical spectrum (
\right) \,.
\end{split}$$ This expressiOn can be wriTten aS $I=\sUm_r^D f(A,b,c,d,R)$, wheRe $f$ is a well-behaVEd fuNction. Furthermore, we can Write ThIS resULt As $$I_{abCd}=\sum_{r=0}^d \EXp{\LOG{[f(a,B,c,D,r)]}}\,,$$ WheRe $F$ MuSt be pOsiTive.\
\
The Logarithm $\lOg{(f)}$ Is Given by $$\begin{SPlIt}
    \log{(f)}=&-\frac{1}{2}\Log{(2)}-2\Log{(\pi)}\\ &+\frac{1}{2}\lefT(\loG{(c!)}+\log{(d!)}
    -\LoG{(a!)}-\lOG{(b!)}\rigHt)\\
    &-\lOg{\lefT((d-r)!\rigHT)}-\log{\leFt((c-d+r)!\righT)}-\lOG{(r!)}\\
    &+\log{\lEFt(\Gamma\LEFt(\Frac{A+b-c+d+1}{2}-r\right)\right)}\\
    &+\loG{\LeFT(\Gamma\left(\frac{A-b+c-d+1}{2}+r\rIgHT)\rIGHt)}\\
    &+\lOg{\lEft(\Gamma\leFt(\Frac{-a+B+C-d+1}{2}+r\righT)\RiGHT)}\,.
\End{SPlit}$$ Notice thaT due to the prESenCe of thE gAmmA FunctiOns, thE pREviOus expressiOn coUld requirE the evALuation OF a logarIthm of A neGatIve qUAnTiTy. FOr THis REaSon, WE caLculate iNsTeAd $\log{(\Abs{\GAMMA(N)})}$, and ComPute The asSociated phase SepAratELy. THis phAse is IntrOdUced iN the fiNal exPrEssion after the eXponEntiation. thuS, tHe fInAl expREssion Of tHe iNtegral Taking iNTo aCcOUNT tHe possible negative VaLUEs Of gamma fUnctioNS rEaDS $$I_{abcd}=\suM_{r=0}^D p\eXp{\loG{[F(A,b,c,d,r)]}}\,,$$ WherE $P$ iS the phasE, given BY $p=\PrOd_{n=1}^3 p_n$, wiTh $P_n$ beinG tHe pHasE geneRAted By each Gamma funCtion.

![eNergy spectrum oF The two-particlE SySTEm AS a fuNctIon of the intEracTIon sTrenGTh. the CAlculAtionS aRE pERformed using 100 single-pArTicle mOdes (sOlid green line). the analytiCAL Spectrum (
\right) \,.
\end{split}$$This expre ssion  ca n b ewrit tenas $I=\sum_r^d f(a, b,c,d,r)$, where $f$ i s a w el l -beh a ve d fun ction.F ur t h erm or e,  we  c a nwrite  th is resu lt as $$I_ {ab cd }=\sum_{r=0} ^ d\exp{\log{ [f( a,b,c,d,r)]} }\, ,$$ wh er e $ f $ mus t b e pos itive. \ 
\
The  logarith m$ \log{( f )}$ isg i ve n by  $$\begin{split} \log{(f)}=&-\f rac{1} {2 } \l o g {(2 )}- 2\log{(\pi )} \\ &+ \ frac{1} { 2} \ l e ft( \ log{(c!)}+\lo g{(d!)}
  -\l og{(a! )} -\l o g{(b!) }\rig ht ) \\    &-\log{ \lef t((d-r)!\ right) } -\log{\ l eft((c- d+r)!\ rig ht) }-\l o g{ (r !)} \\   &+ \lo g {\l eft(\Gam ma \l eft(\ frac { a + b -c+d +1} {2}- r\rig ht)\right)}\\ 
   &+ \ log {\lef t(\Ga mma\ le ft(\f rac{a- b+c-d +1 }{2}+r\right)\r ight )}\\
   &+\ lo g{\ le ft(\G a mma\le ft( \fr ac{-a+b +c-d+1} { 2}+ r\ r i g ht )\right)}\,.
\end{ sp l i t} $$ Notic e that du et o the pr es enc e of t he ga mmaf un ctions,the pr e vi ou s expre ss ion co ul d r equ ire t h e ev aluati on of alogar i thm of a negat i ve quantity.F or t hi s  rea son , we calcul atei nste ad $ \ lo g{( \ abs{\ Gamma (n ) }) } $, and compute theas sociat ed ph ase separatel y. This ph a s e  is intr oduc e di n the final ex press ion aftert he expon entia tion. Th us, the f i n al expre ssi onofthe i nt egral takingi n to a cc ount th e p ossible  ne gat ive  va lu es of gam ma funct io ns  r ea ds$$I_{ a bcd}=\su m_ {r= 0} ^dp\exp { \log{[ f(a,b ,c,d ,r )] } }\, ,$$ whe r e$ p $ is  t he  pha se,  g ivenby $ p =\p rod_{n= 1}^3 p_n$ , w i th $ p_ n$  beingthe phase gen er ated by ea ch  ga mma fu n c tion.

! [Energy spectrum of the two-par tic le sy stem  as a fun cti on ofthe intera ctionstren gt h.T h e cal c u la tio ns  are perfo r m edusing  1 00 s ingle-p article modes (sol i d g reen line). T heanal y t ic als pe c tru m( 
\right) \,.
\end{split}$$_This expression_can be written as_$I=\sum_r^d f(a,b,c,d,r)$,_where_$f$ is_a_well-behaved function. Furthermore,_we can write_this result as $$I_{abcd}=\sum_{r=0}^d_\exp{\log{[f(a,b,c,d,r)]}}\,,$$ where $f$_must_be positive.\
\
The logarithm $\log{(f)}$ is given by $$\begin{split}
    \log{(f)}=&-\frac{1}{2}\log{(2)}-2\log{(\pi)}\\ &+\frac{1}{2}\left(\log{(c!)}+\log{(d!)}
 __ -\log{(a!)}-\log{(b!)}\right)\\
___ &-\log{\left((d-r)!\right)}-\log{\left((c-d+r)!\right)}-\log{(r!)}\\
   _&+\log{\left(\Gamma\left(\frac{a+b-c+d+1}{2}-r\right)\right)}\\
    &+\log{\left(\Gamma\left(\frac{a-b+c-d+1}{2}+r\right)\right)}\\
_ _ &+\log{\left(\Gamma\left(\frac{-a+b+c-d+1}{2}+r\right)\right)}\,.
\end{split}$$ Notice that due to the presence_of_the gamma functions,_the previous expression could require the evaluation of a_logarithm of a negative quantity. For_this reason, we_calculate_instead_$\log{(\abs{\Gamma(n)})}$, and compute the_associated phase separately. This phase is_introduced in the final expression after_the exponentiation. Thus, the final expression of_the integral taking into account the_possible negative values of gamma_functions reads_$$I_{abcd}=\sum_{r=0}^d p\exp{\log{[f(a,b,c,d,r)]}}\,,$$ where $p$ is_the phase, given_by $p=\prod_{n=1}^3_p_n$, with $p_n$_being the phase generated by each_gamma function.

![Energy spectrum_of the two-particle system as a_function_of the interaction_strength._The_calculations are_performed using 100_single-particle_modes (solid_green_line). The analytical spectrum (