k,\lambda)$ at $k^*$. This makes the search for the zeros of $\mathbold{S}(k,\lambda)^{-1}$ a numerically hard problem, particularly when $k \ll 1$. In that case $\mathbold{S}(k,\lambda)$ has zeros and poles close to each other and it does not satisfy the smoothness conditions required to ensure rapid convergence if, for instance, Newton iterations are used to find the zeros.

For single-channel radial problems, as in section \[sphericalsymmetry\], it was proposed in [@Broeckhove2009] to use a regularized function rather than the $S$-matrix itself to converge and to track down bound states and resonances. This function was defined as, per $l$-channel, $$\label{eq:one_channel_regularized}
    F_{l}(k,\lambda) = \frac{k^{2l+1}}{S_l(k,\lambda)-1},$$ with $S_{l}$ as in. It was shown that for $k\neq 0 $ this function has a zero iff $S_l(k,\lambda)$ has a pole. In addition this function is bounded as long as $\int_0^\infty dr\, \hat{j}_l(kr) V(r,\lambda) \psi_l(k,r) \neq 0$, which does not hold only for very specific potentials.

In this paper we extend this result to the more general class of coupled channel problems described in the previous section. We introduce, in extension of, the $(l_\text{max} + 1)$ by $(l_\text{max} + 1)$ square matrix function $$\label{eq:multi_channel_regularized}
    \mathbold{F}(k,\lambda) = \mathcal{K}
    \left( \mathbold{S}(k,\lambda) - \mathbold{I} \right)^{-1},$$ where $\mathcal{K}$ is a diagonal matrix with elements $(\mathcal{K})_{ii} = k^{2l_{i}+1}$ on the diagonal and where we extract the $S$-matrix using from the numerical solution of, and. The function used for continuation then is $$\label{eq:det_multi_channel_regularized}
    \det(\mathbold{F}(k,\lambda)) = \left( \prod_{l=0}^{l_{\text{max}}} k^{
k,\lambda)$ at $ k^*$. This makes the search for the zeros of $ \mathbold{S}(k,\lambda)^{-1}$ a numerically hard problem, peculiarly when $ thousand \ll 1$. In that case $ \mathbold{S}(k,\lambda)$ has zero and pole close to each other and it does not meet the smoothness conditions required to ensure rapid overlap if, for case, Newton iterations are use to find the zeros. 

 For unmarried - distribution channel radial problems, as in section \[sphericalsymmetry\ ], it was proposed in [ @Broeckhove2009 ] to use a regularized affair rather than the $ S$-matrix itself to converge and to track down bound state and resonances. This function was defined as, per $ l$-channel, $ $ \label{eq: one_channel_regularized } 
     F_{l}(k,\lambda) = \frac{k^{2l+1}}{S_l(k,\lambda)-1},$$ with $ S_{l}$ as in. It was shown that for $ k\neq 0 $ this affair has a zero iff $ S_l(k,\lambda)$ has a pole. In addition this function is restrict as long as $ \int_0^\infty dr\, \hat{j}_l(kr) V(r,\lambda) \psi_l(k, r) \neq 0 $, which does not hold only for very specific potentials. 

 In this paper we extend this result to the more general class of coupled channel problems described in the previous section. We introduce, in propagation of, the $ (l_\text{max } + 1)$ by $ (l_\text{max } + 1)$ square matrix function $ $ \label{eq: multi_channel_regularized } 
     \mathbold{F}(k,\lambda) = \mathcal{K } 
     \left (\mathbold{S}(k,\lambda) - \mathbold{I } \right)^{-1},$$ where $ \mathcal{K}$ is a aslant matrix with chemical element $ (\mathcal{K})_{ii } = k^{2l_{i}+1}$ on the aslant and where we extract the $ S$-matrix using from the numerical solution of, and. The function used for continuation then is $ $ \label{eq: det_multi_channel_regularized } 
     \det(\mathbold{F}(k,\lambda) ) = \left (\prod_{l=0}^{l_{\text{max } } } k^ {
k,\lalbda)$ at $k^*$. This makes the search for the zeros mf $\matgbold{S}(k,\lxmbda)^{-1}$ a numerically hard proulem, partucularly when $k \ll 1$. In that casv $\mathbole{S}(k,\lembda)$ has zeros ehd poles closs to zarh other and it does not vatisfy the smmoghuess conditions required to ensure rwpid comvfrgence if, for insuanse, Nswton iterations are used to find fhe zerms.

For single-cnannel radial problems, as ln sfction \[sphericalsylmetry\], it wqs pwiposed in [@Brueckhove2009] to use a regumarized function rather than thd $S$-macrix itself ti cltverge and vo trabk down bound states dnd respnances. This fmnctimn qas defined as, per $l$-ciannel, $$\label{eq:one_chagnel_reguldrnzed}
    F_{l}(k,\lambda) = \frqc{j^{2l+1}}{S_l(k,\nambga)-1},$$ wkrh $R_{l}$ zs ih. It wws ahown that for $k\neq 0 $ this function has s svto iff $S_l(k,\lajbda)$ hws a pole. In addition this function is bolndes as long as $\int_0^\infty de\, \hat{j}_l(kr) V(r,\lambda) \psl_l(k,r) \neq 0$, which does not hold only for very specific potendials.

Mn thnw papde ae extend this result to the more general clafa pf coupled channcl problems descrinef og the previour sectnkn. We introduce, in edtensiog of, rhe $(l_\text{iax} + 1)$ by $(l_\text{max} + 1)$ square matrux function $$\jqbel{eq:multi_channel_xegularized}
    \msthbokd{F}(k,\lambda) = \mathcal{K}
    \leff( \mathbold{S}(n,\lambda) - \jxthbold{I} \right)^{-1},$$ wferv $\madhcal{K}$ is a diagonal matriv with elxmentx $(\mathcxl{K})_{io} = k^{2l_{i}+1}$ on the diwgonal and where we extrwct tke $S$-mdtrix usinh from the numerical solution oh, and. The funvthon used fox contlnuation then if $$\label{eq:det_mukti_chanuel_regjlarized}
    \det(\mavhbold{F}(k,\lambqa)) = \left( \prod_{n=0}^{p_{\text{max}}} k^{
k,\lambda)$ at $k^*$. This makes the search zeros $\mathbold{S}(k,\lambda)^{-1}$ a hard problem, particularly that $\mathbold{S}(k,\lambda)$ has zeros poles close to other and it does not satisfy smoothness conditions required to ensure rapid convergence if, for instance, Newton iterations are to find the zeros. For single-channel radial problems, as in section \[sphericalsymmetry\], it proposed [@Broeckhove2009] use regularized function rather than the $S$-matrix itself to converge and to track down bound states and This function was defined as, per $l$-channel, $$\label{eq:one_channel_regularized} = \frac{k^{2l+1}}{S_l(k,\lambda)-1},$$ with $S_{l}$ in. It was shown that $k\neq $ this has zero $S_l(k,\lambda)$ has a In addition this function is bounded as long as $\int_0^\infty dr\, \hat{j}_l(kr) V(r,\lambda) \psi_l(k,r) \neq 0$, which not hold very specific In paper extend this result more general class of coupled channel the previous section. We introduce, in extension of, $(l_\text{max} + by $(l_\text{max} + 1)$ square matrix $$\label{eq:multi_channel_regularized} \mathbold{F}(k,\lambda) = \mathcal{K} \left( \mathbold{S}(k,\lambda) - \mathbold{I} where $\mathcal{K}$ is a diagonal matrix with elements $(\mathcal{K})_{ii} = k^{2l_{i}+1}$ on the diagonal and extract the $S$-matrix using the numerical solution and. function for then is \det(\mathbold{F}(k,\lambda)) = \left( \prod_{l=0}^{l_{\text{max}}} k^{
k,\lambda)$ at $k^*$. This makes the seaRch for the zEros oF $\maThbOlD{S}(k,\lAmbdA)^{-1}$ a numerically hARd prOblem, particularly when $k \Ll 1$. In tHaT Case $\MAtHbold{s}(k,\lambdA)$ HaS ZEroS aNd PolEs CLoSe to eAch Other anD it does not SatIsFy the smoothnESs Conditions ReqUired to ensurE raPid conVeRgeNCe if, fOr iNstanCe, NewtON iteraTions are uSeD To find THe zeros.

fOR sInglE-channel radial proBLeMS, as in section \[spHericaLsYMmETRy\], iT waS proposed iN [@BRoeckHOve2009] to usE A rEGULarIZed function raTher than the $s$-MatRix itsElF to COnvergE and tO tRAck Down bound stAtes And resonaNces. ThIS functiON was defIned as, Per $L$-chAnneL, $$\LaBeL{eq:OnE_ChaNNeL_reGUlaRized}
    F_{l}(k,\LaMbDa) = \fraC{k^{2l+1}}{S_L(K,\LAMbda)-1},$$ WitH $S_{l}$ aS in. It Was shown that fOr $k\Neq 0 $ tHIs fUnctiOn has A zerO iFf $S_l(k,\Lambda)$ Has a pOlE. In addition this FuncTion is bouNdeD aS loNg As $\int_0^\INfty dr\, \Hat{J}_l(kR) V(r,\lambDa) \psi_l(k,R) \Neq 0$, WhICH DoEs not hold only for veRy SPEcIfic poteNtials.

iN tHiS Paper we eXtEnd This RESult tO the MOrE general Class oF CoUpLed chanNeL problEmS deScrIbed iN The pReviouS section. we intROduce, in extensiON of, the $(l_\text{maX} + 1)$ By $(L_\TExT{Max} + 1)$ sQuaRe matrix funCtioN $$\LabeL{eq:mULtI_chANnel_rEgulaRiZEd}
    \MAthbold{F}(k,\lambda) = \mathCaL{K}
    \left( \MathbOld{S}(k,\lambda) - \maThbold{I} \rigHT)^{-1},$$ WHere $\mathCal{K}$ IS a DIagonal matrix wIth elEments $(\mathCAl{K})_{ii} = k^{2l_{i}+1}$ On the Diagonal And where wE EXtract thE $S$-mAtrIx uSinG FRoM the numerical SOLutiOn Of, and. ThE fuNction uSed For ConTinUaTion then iS $$\label{eq:DeT_mUlTi_ChaNnel_rEGularizeD}
    \dEt(\mAtHboLd{F}(k,\lAMbda)) = \leFt( \proD_{l=0}^{l_{\tExT{mAX}}} k^{
k,\lambda)$ at $k^*$. This  makes the  sear chfor  t he z eros  of $\mathbold { S}(k ,\lambda)^{-1}$ a nume rical ly hard pr oblem , parti c ul a r lywh en  $k  \ l l1$. I n t hat cas e $\mathbo ld{ S} (k,\lambda)$ ha s zeros an d p oles close t o e ach ot he r a n d itdoe s not  satis f y thesmoothnes sc onditi o ns requ i r ed  toensure rapid conv e rg e nce if, for in stance ,N ew t o n i ter ations are  u sed t o  find t h ez e r os. 
For single-c hannel radi a l p roblem s,  as in sec tion\[ s phe ricalsymmet ry\] , it waspropos e d in [@ B roeckho ve2009 ] t o u se a re gu lar iz e d f u nc tio n  ra ther tha nth e $S$ -mat r i x itse lfto c onver ge and to tra ckdown bou nd st atesandre sonan ces. T his f un ction was defin ed a s, per $l $-c ha nne l,  $$\l a bel{eq :on e_c hannel_ regular i zed } F _{l}(k,\lambda) =\f r a c{ k^{2l+1} }{S_l( k ,\ la m bda)-1}, $$  wi th $ S _ {l}$as i n .It was s hown t h at  f or $k\n eq  0 $ t hi s f unc tionh as a  zeroiff $S_l (k,\l a mbda)$ has a p o le. In additi o nt h is func tio n is bounde d as long  as$ \i nt_ 0 ^\inf ty dr \, \h a t{j}_l(kr) V(r,\lam bd a) \ps i_l(k ,r) \neq 0$,which does n o t hold o nlyf or very specificpoten tials.

In this pap er we  extendthis resu l t  to themor e g ene ral c la ss of coupled c hann el  proble msdescrib edinthe  pr ev ious sect ion. Wein tr od uc e,in ex t ension o f,  th e$(l _\tex t {max}+ 1)$  by$( l_ \ tex t{max}+  1 ) $  squ ar ematr ixfu nctio n $$ \ lab el{eq:m ulti_chan nel _ regu la ri zed}
   \mathbold{F }( k,\lambda)  =  \m athcal { K }
    \l eft( \mathbold{S}(k,\la m bda) -\ma thbol d{I}  \right)^ {-1 },$$ w her e  $\mat hcal{K }$ is  a  di a g onalm a tr ixwi th element s $(\ mathc al {K}) _{ii} =  k^{2l_{i}+1}$ ont hediagonal andwhe re w e ex tra c tt he$S $ -ma t r ix using from t he numeric al so lution of, and .The fun ction u sed f o r conti nuation t hen is $$ \l abel { e q:d et_multi_c hannel_r egularize d }
  \d et(\m ath bold{F }( k,\ lambd a)) =\ lef t( \p rod_{l =0 }^{l_{ \text {m ax}}} k^ {
k,\lambda)$ at_$k^*$. This_makes the search for_the zeros_of_$\mathbold{S}(k,\lambda)^{-1}$ a_numerically_hard problem, particularly_when $k \ll_1$. In that case_$\mathbold{S}(k,\lambda)$ has zeros_and_poles close to each other and it does not satisfy the smoothness conditions required_to_ensure rapid_convergence_if,_for instance, Newton iterations are_used to find the zeros.

For_single-channel radial_problems, as in section \[sphericalsymmetry\], it was proposed_in_[@Broeckhove2009] to use_a regularized function rather than the $S$-matrix itself to_converge and to track down bound_states and resonances._This_function_was defined as, per_$l$-channel, $$\label{eq:one_channel_regularized}
    F_{l}(k,\lambda)_= \frac{k^{2l+1}}{S_l(k,\lambda)-1},$$ with $S_{l}$ as in._It was shown that for $k\neq 0_$ this function has a zero_iff $S_l(k,\lambda)$ has a pole._In addition_this function is bounded as_long as $\int_0^\infty_dr\, \hat{j}_l(kr)_V(r,\lambda) \psi_l(k,r) \neq_0$, which does not hold only_for very specific_potentials.

In this paper we extend this_result_to the more_general_class_of coupled_channel problems described_in_the previous_section._We introduce, in extension of, the_$(l_\text{max}_+ 1)$ by $(l_\text{max} + 1)$ square_matrix function $$\label{eq:multi_channel_regularized}
 __ \mathbold{F}(k,\lambda) = \mathcal{K}
_   \left( \mathbold{S}(k,\lambda)_- \mathbold{I} \right)^{-1},$$ where $\mathcal{K}$ is_a diagonal_matrix with_elements $(\mathcal{K})_{ii} = k^{2l_{i}+1}$ on the diagonal and where we extract_the $S$-matrix using from the numerical_solution of, and. The_function used_for_continuation then is_$$\label{eq:det_multi_channel_regularized}
_ _ \det(\mathbold{F}(k,\lambda)) = \left( \prod_{l=0}^{l_{\text{max}}} k^{