LRb}/N_c;&M^{d4}_{LL}&=&A^F_{SPb}.
\end{array}\end{aligned}$$ where $T^F_{Xa}$ and $T^F_{Xb}$ ($X=LL,LR,SP$) represent the factorizable emission diagram contributions, $T^N_{Xa}$ and $T^N_{Xb}$ ($X=LL,LR,SP$) are the non-factorizable emission diagram contributions. $A^F_{Xa}$, $A^F_{Xb}$ and $A^N_{Xa}$, $A^N_{Xb}$ ($X=LL,LR,SP$) denote the so-called factorizable and non-factorizable annihilation diagram contributions respectively. Their detailed definitions and general formalisms are presented in the Appendix.

Treatment of Singularities  {#sec:TOD}
===========================

In the evaluation of hadronic matrix elements, there are two kinds of singularities, one is caused by the infrared divergence of gluon exchanging interaction, and the other arises from the on-mass shell divergence of internal quark propagator. As the quark propagator singularity is a physical-region singularity, one can simply add $i\epsilon$ to the denominator of quark propagator and apply the Cutkosky rule [@cutkosky] to avoid such a singularity. It then allows us to obtain the virtual part of amplitudes as the Cutkosky rule gives a compact expression for the discontinuity across the cut arising from a physical-region singularity. In general, a Feynman diagram will yield an imaginary part for the decay amplitudes when the virtual particles in the diagram become on mass-shell, and the resulting diagram can be considered as a genuine physical process. It is well-known that when applying the Cutkosky rule to deal with a physical-region singularity of all propagators, the following formula holds $$\begin{aligned}
\frac{1}{p^2-m_b^2+i\epsilon}&=&P\biggl[\frac{1}{p^2-m_b^2}
\biggl]-i\pi\delta[p^2-m_b^2],\nonumber\\
\frac{1}{p^2-m_q^2+i\epsilon}&=&P\biggl[\frac{1}{p^2-m_q^2}
\biggl]-i\pi\delta[p^2
LRb}/N_c;&M^{d4}_{LL}&=&A^F_{SPb }. 
 \end{array}\end{aligned}$$ where $ T^F_{Xa}$ and $ T^F_{Xb}$ ($ X = LL, LR, SP$) represent the factorizable emission diagram contributions, $ T^N_{Xa}$ and $ T^N_{Xb}$ ($ ten = LL, LR, SP$) are the non - factorizable discharge diagram contributions. $ A^F_{Xa}$, $ A^F_{Xb}$ and $ A^N_{Xa}$, $ A^N_{Xb}$ ($ X = LL, LR, SP$) denote the then - called factorizable and non - factorizable annihilation diagram contributions respectively. Their detailed definition and cosmopolitan formalisms are presented in the Appendix. 

 Treatment of Singularities   { # sec: TOD } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In the evaluation of hadronic matrix elements, there are two kinds of singularity, one is caused by the infrared divergence of gluon commute interaction, and the early arises from the on - mass shell deviation of internal quark propagator. As the quark propagator singularity is a physical - area singularity, one can simply add $ i\epsilon$ to the denominator of quark propagator and apply the Cutkosky principle   [ @cutkosky ] to avoid such a singularity. It then allows us to receive the virtual part of amplitude as the Cutkosky rule gives a compact expression for the discontinuity across the cut get up from a physical - region singularity. In general, a Feynman diagram will yield an imaginary part for the decay amplitudes when the virtual particles in the diagram become on mass - shell, and the resulting diagram can be considered as a genuine physical process. It is well - known that when applying the Cutkosky convention to deal with a physical - area singularity of all propagator, the follow formula holds $ $ \begin{aligned } 
 \frac{1}{p^2 - m_b^2+i\epsilon}&=&P\biggl[\frac{1}{p^2 - m_b^2 } 
 \biggl]-i\pi\delta[p^2 - m_b^2],\nonumber\\ 
 \frac{1}{p^2 - m_q^2+i\epsilon}&=&P\biggl[\frac{1}{p^2 - m_q^2 } 
 \biggl]-i\pi\delta[p^2
LRb}/J_c;&M^{d4}_{LL}&=&A^F_{SPb}.
\end{array}\end{allgned}$$ where $T^F_{Xa}$ and $T^F_{Eb}$ ($X=LL,LD,SP$) reprdsent the factorizable emissmon eiagrqm contributions, $T^N_{Xa}$ xnd $T^N_{Xb}$ ($Q=LL,LR,SP$) aee tie non-factorizable emissljn djwgrak contributions. $A^F_{Xa}$, $A^F_{Xb}$ and $A^N_{Xa}$, $A^N_{Xb}$ ($X=UL,PR,SP$) denote the so-called factorizabje and moj-factorizable wnninylatjon diagram contributions respectibely. Thtir detailed definotions and general formalidms wre presented in tje Appendix.

Jdeaenent of Singjlarities  {#sec:TOD}
===========================

In thg evaluation of hadronic matrix dlemeuts, there atz twl kinds of smngulagities, one is caused ty the onfrared divernence of gluon exchanging intxraction, and the othgr arises xrkm the on-mass sheol divetgenca of untdrnzl qhark pgopegator. As tge quark pripagator singularitu yw a physical-rsgion figgularity, one can simply add $i\epsilon$ tm tge denominator of quark propagator and apply the Cutkjsky rule [@cutkosky] to avoid such a singularity. It dhen eluoww ms tu ohtain the virtual part of amplitudes as the Ctfkpsly rule gives c compact expresxiln sor the discovtinuicg zcross the cut ariding frjm a physical-rtgion singularity. In general, a Fwynman diagrcm qill yield an imagnnary part fur tne devay amplitudes when the virfual particpes in ths diagram become un kavs-shell, and the resulting qiagram cen be considdred as a denuine phjsical process. It is welp-knowu thad when apppying the Cutkosky rule to deal with a physicsl-segpon singujaritn of all propagwtors, the follpwing fjrmulx holds $$\befin{aligied}
\frac{1}{p^2-m_b^2+i\e[silon}&=&P\biggl[\fswc{1}{p^2-m_b^2}
\biggl]-i\pm\delta[p^2-m_b^2],\gonunber\\
\drac{1}{p^2-m_q^2+k\dpsilon}&=&P\biggl[\ftac{1}{p^2-m_q^2}
\bigyj]-i\pi\eelta[p^2
LRb}/N_c;&M^{d4}_{LL}&=&A^F_{SPb}. \end{array}\end{aligned}$$ where $T^F_{Xa}$ and $T^F_{Xb}$ ($X=LL,LR,SP$) factorizable diagram contributions, and $T^N_{Xb}$ ($X=LL,LR,SP$) contributions. $A^F_{Xb}$ and $A^N_{Xa}$, ($X=LL,LR,SP$) denote the factorizable and non-factorizable annihilation diagram contributions Their detailed definitions and general formalisms are presented in the Appendix. Treatment of {#sec:TOD} =========================== In the evaluation of hadronic matrix elements, there are two kinds singularities, is by infrared divergence of gluon exchanging interaction, and the other arises from the on-mass shell divergence of quark propagator. As the quark propagator singularity is physical-region singularity, one can add $i\epsilon$ to the denominator quark and apply Cutkosky [@cutkosky] avoid such a It then allows us to obtain the virtual part of amplitudes as the Cutkosky rule gives a expression for across the arising a singularity. In general, diagram will yield an imaginary part amplitudes when the virtual particles in the diagram on mass-shell, the resulting diagram can be considered a genuine physical process. It is well-known that applying the Cutkosky rule to deal with a physical-region singularity of all propagators, the following $$\begin{aligned} \frac{1}{p^2-m_b^2+i\epsilon}&=&P\biggl[\frac{1}{p^2-m_b^2} \biggl]-i\pi\delta[p^2-m_b^2],\nonumber\\ \frac{1}{p^2-m_q^2+i\epsilon}&=&P\biggl[\frac{1}{p^2-m_q^2}
LRb}/N_c;&M^{d4}_{LL}&=&A^F_{SPb}.
\end{array}\end{Aligned}$$ wheRe $T^F_{XA}$ anD $T^F_{xb}$ ($x=LL,Lr,SP$) rEpresent the facTOrizAble emission diagram conTribuTiONs, $T^N_{xA}$ aNd $T^N_{XB}$ ($X=LL,LR,Sp$) ArE THe nOn-FaCtoRiZAbLe emiSsiOn diagrAm contribuTioNs. $a^F_{Xa}$, $A^F_{Xb}$ and $A^n_{xa}$, $a^N_{Xb}$ ($X=LL,LR,Sp$) deNote the so-calLed FactorIzAblE And noN-faCtoriZable aNNihilaTion diagrAm COntribUTions reSPEcTiveLy. Their detailed deFInITions and generaL formaLiSMs ARE prEseNted in the APpEndix.

tReatmenT Of sINGulARities  {#sec:TOD}
===========================

IN the evaluatIOn oF hadroNiC maTRix eleMents, ThERe aRe two kinds oF sinGularitieS, one is CAused by THe infraRed divErgEncE of gLUoN eXchAnGIng INtEraCTioN, and the oThEr AriseS froM THE On-maSs sHell DiverGence of internAl qUark PRopAgatoR. As thE quaRk PropaGator sIngulArIty is a physical-rEgioN singularIty, OnE caN sImply ADd $i\epsIloN$ to The denoMinator OF quArK PROpAgator and apply the CUtKOSkY rule [@cutKosky] tO AvOiD Such a sinGuLarIty. IT THen alLows US tO obtain tHe virtUAl PaRt of ampLiTudes aS tHe CUtkOsky rULe giVes a coMpact expRessiON for the discontINuity across thE CuT ARiSIng fRom A physical-reGion SInguLariTY. IN geNEral, a feynmAn DIaGRam will yield an imagiNaRy part For thE decay amplituDes when the VIRTual partIcleS In THe diagram becomE on maSs-shell, and THe resultIng diAgram can Be consideRED as a genuIne PhySicAl pROCeSs. It is well-knoWN That WhEn applyIng The CutkOskY ruLe tO deAl With a physIcal-regiOn SiNgUlAriTy of aLL propagaToRs, tHe FolLowinG FormulA holdS $$\begIn{AlIGneD}
\frac{1}{p^2-m_B^2+I\ePSIlon}&=&p\bIgGl[\frAc{1}{p^2-M_b^2}
\Biggl]-I\pi\dELta[P^2-m_b^2],\nonuMber\\
\frac{1}{p^2-M_q^2+i\EPsilOn}&=&p\bIggl[\fraC{1}{p^2-m_q^2}
\biggl]-i\pi\dElTa[p^2
LRb}/N_c;&M^{d4}_{LL}&=&A^ F_{SPb}.
\ end{a rra y}\ en d{al igne d}$$ where $T^ F _{Xa }$ and $T^F_{Xb}$ ($X= LL,LR ,S P $) r e pr esent  the fa c to r i zab le  e mis si o ndiagr amcontrib utions, $T ^N_ {X a}$ and $T^N _ {X b}$ ($X=LL ,LR ,SP$) are th e n on-fac to riz a ble e mis siondiagra m  contr ibutions.  $ A ^F_{Xa } $, $A^F _ { Xb }$ a nd $A^N_{Xa}$, $A ^ N_ { Xb}$ ($X=LL,LR ,SP$)de n ot e the  so -called fa ct oriza b le andn on - f a cto r izable annihi lation diag r amcontri bu tio n s resp ectiv el y . T heir detail ed d efinition s andg eneralf ormalis ms are  pr ese nted in  t heAp p end i x. 

T r eat ment ofSi ng ulari ties { # sec: TOD }
== ===== ============= === ==== 
In  theevalu atio nof ha dronic  matr ix  elements, ther e ar e two kin dsof  si ng ulari t ies, o neiscausedby thei nfr ar e d di vergence of gluonex c h an ging int eracti o n,  a n d the ot he r a rise s fromtheo n- mass she ll div e rg en ce of i nt ernalqu ark  pr opaga t or.As the  quark p ropag a tor singularit y  is a physica l -r e g io n  sin gul arity, onecans impl y ad d  $ i\e p silon $ toth e  d e nominator of quarkpr opagat or an d apply the C utkosky ru l e [@cutkos ky]t oa void such a si ngula rity. It t h en allow s usto obtai n the vir t u al partofamp lit ude s as  the Cutkosky r ulegi ves a c omp act exp res sio n f orth e discont inuity a cr os sth e c ut ar i sing fro ma p hy sic al-re g ion si ngula rity .In gen eral, a Fe y n mandi ag ramwil lyield  ani mag inary p art for t hed ecay  a mp litudes  when the vir tu al particl es  in  the d i a gram bec ome on mass-shell, andt he resu lti ng di agra m can becon sidere d a s  a gen uine p hysic al  pr o c ess.I t  i s w el l-known th a t  wh en ap pl ying  the Cu tkosky rule to dea l  wi th a physical -re gion s in gul a ri t y o fa llp r opagators, thefollowingfo r mu la holds $ $ \be gi n{align ed}
\fr ac{1} { p^2-m_b ^2+i\epsi lon}&=&P\ bi ggl[ \ f rac {1}{p^2-m_ b^2}
\bi ggl]-i\pi \ delta [ p^ 2-m_b ^2] ,\nonu mb er\ \
\fr ac{1}{ p ^2- m_q^2 +i\eps il on}&=& P\big gl [\frac{1 }{p^2-m_q^2}
\biggl]-i\ pi\del ta[p^ 2
LRb}/N_c;&M^{d4}_{LL}&=&A^F_{SPb}.
\end{array}\end{aligned}$$ where_$T^F_{Xa}$ and_$T^F_{Xb}$ ($X=LL,LR,SP$) represent the_factorizable emission_diagram_contributions, $T^N_{Xa}$_and_$T^N_{Xb}$ ($X=LL,LR,SP$) are_the non-factorizable emission_diagram contributions. $A^F_{Xa}$, $A^F_{Xb}$_and $A^N_{Xa}$, $A^N_{Xb}$_($X=LL,LR,SP$)_denote the so-called factorizable and non-factorizable annihilation diagram contributions respectively. Their detailed definitions and_general_formalisms are_presented_in_the Appendix.

Treatment of Singularities _{#sec:TOD}
===========================

In the evaluation of hadronic_matrix elements,_there are two kinds of singularities, one is_caused_by the infrared_divergence of gluon exchanging interaction, and the other arises_from the on-mass shell divergence of_internal quark propagator._As_the_quark propagator singularity is_a physical-region singularity, one can simply_add $i\epsilon$ to the denominator of_quark propagator and apply the Cutkosky rule [@cutkosky]_to avoid such a singularity. It_then allows us to obtain_the virtual_part of amplitudes as the_Cutkosky rule gives_a compact_expression for the_discontinuity across the cut arising from_a physical-region singularity._In general, a Feynman diagram will_yield_an imaginary part_for_the_decay amplitudes_when the virtual_particles_in the_diagram_become on mass-shell, and the resulting_diagram_can be considered as a genuine physical_process. It is well-known_that_when applying the Cutkosky_rule to deal with a_physical-region singularity of all propagators, the_following formula_holds $$\begin{aligned}
\frac{1}{p^2-m_b^2+i\epsilon}&=&P\biggl[\frac{1}{p^2-m_b^2}
\biggl]-i\pi\delta[p^2-m_b^2],\nonumber\\
\frac{1}{p^2-m_q^2+i\epsilon}&=&P\biggl[\frac{1}{p^2-m_q^2}
\biggl]-i\pi\delta[p^2