 for massless HISQ quarks with NRQCD formulated in a moving frame (mNRQCD) was undertaken for the vector and tensor heavy-light currents in [@mueller11].

For massive quarks similar matching calculations using the same lattice action for both quarks have been carried out for Wilson quarks in [@kuramashi98] and for various implementations of NRQCD in [@braaten95; @jones99; @boyle00; @hart07]. To our knowledge, no matching calculations with mixed actions and massive relativistic quarks have been reported in the literature.

Moving from massless to massive relativistic quarks complicates the matching procedure. In the former case, quarks and antiquarks at zero spatial momentum are indistinguishable and consequently scattering and annihilation processes give identical results. In the massive case, however, we must distinguish between quarks and antiquarks. For HISQ quarks at zero spatial momentum, this corresponds to choosing $ap_0 = ia{m_{\text{tree}}}$ or $ap_0=-ia{m_{\text{tree}}}$ respectively. We choose outgoing quarks or antiquarks — the “scattering” or “annihilation” channels respectively — to ensure we do not attempt to compute vanishing matrix elements. Thus we calculate the matrix elements of $V_0$ and $A_k$ in the scattering channel and of $A_0$ and $V_k$ in the annihilation channel. This procedure is valid, even at nonzero lattice spacing, provided we match to the same channel in continuum QCD.

Unfortunately, from the practical viewpoint of calculating Feynman diagrams, using massive quarks complicates the numerical integration considerably. The chief difficulty lies in the annihilation channel, which contains a Coulomb singularity that must be handled with a subtraction function. We discuss the subtraction functions employed in this work in more detail in Appendix \[app:subfn\]. Furthermore, in the automated perturbation theory routines, the pole in the NRQCD propagator crosses the contour of integration and we can no longer carry out the usual Wick rotation back to Minkowski space. Instead, we must deform the integration contours and introduce a triple contour to ensure the stability of numerical integration [@hart07; @mueller11].

For both channels, the lattice matrix elements must be matched to their continuum QCD counterparts. Analytic expressions for the relevant QCD contributions already exist in the literature. References [@braaten95] and [@jones99] discuss the annihilation channel for the axial-vector current, whilst [@boyle00] present results for both
for massless HISQ quarks with NRQCD formulated in a moving frame (mNRQCD) was undertake for the vector and tensor dense - light currents in [ @mueller11 ]. 

 For massive quarks like matching calculations using the like lattice action for both quarks have been carry out for Wilson quarks in [ @kuramashi98 ] and for various implementation of NRQCD in [ @braaten95; @jones99; @boyle00; @hart07 ]. To our knowledge, no matching calculation with assorted actions and massive relativistic quarks have been reported in the literature. 

 Moving from massless to massive relativistic quarks complicate the matching procedure. In the former case, quark cheese and antiquarks at zero spatial momentum are indistinguishable and consequently scattering and annihilation processes yield identical results. In the massive case, however, we must distinguish between quarks and antiquarks. For HISQ quark at zero spatial momentum, this corresponds to choosing $ ap_0 = ia{m_{\text{tree}}}$ or $ ap_0=-ia{m_{\text{tree}}}$ respectively. We choose outgoing quarks or antiquarks — the “ scattering ” or “ annihilation ” channels respectively — to ensure we do not attempt to compute vanishing matrix elements. Thus we calculate the matrix elements of $ V_0 $ and $ A_k$ in the scattering channel and of $ A_0 $ and $ V_k$ in the annihilation groove. This procedure is valid, even at nonzero lattice spacing, provided we equal to the like distribution channel in continuum QCD. 

 Unfortunately, from the practical viewpoint of calculating Feynman diagrams, using massive quarks complicates the numerical consolidation considerably. The chief difficulty lies in the annihilation channel, which contains a Coulomb singularity that must be handled with a subtraction function. We discuss the subtraction functions employed in this work in more detail in Appendix \[app: subfn\ ]. Furthermore, in the automatize perturbation hypothesis routines, the pole in the NRQCD propagator crosses the contour of integration and we can no longer carry out the common Wick rotation back to Minkowski space. Instead, we must deform the integration contours and introduce a triple contour to ensure the constancy of numerical integration [ @hart07; @mueller11 ]. 

 For both channels, the lattice matrix elements must be match to their continuum QCD counterparts. Analytic expression for the relevant QCD contributions already exist in the literature. address [ @braaten95 ] and [ @jones99 ] discuss the annihilation distribution channel for the axial - vector current, whilst [ @boyle00 ] present resultant role for both
 fog massless HISQ quarks wlth NRQCD formulcred in a movjng framd (mNRQCD) was undertaken for vhe cectoe and tensor heavy-lighg currentd in [@mueoler11].

Hor massive quarks similar matdming eaoculations usikg the same lattice actiot wox both quarks have been carried out sor Wilxoj quarks in [@kutamasny98] ans for various implementations of NDQCD in [@braaten95; @jonex99; @boyle00; @hart07]. To our knowlefge, jo matching calculwtions with mixqe actions ana massive gzlativistic quarks have been reported in tfe licerature.

Movunt fgmm massless to mwssive relatlnistic xuarks vomplicates thc matrhint procedure. In the focmer case, quarks and antiquarns at zero spatial nonentuk ara inaustknghixhzble ajd ronsequentlg scatterint and annihilation krosvxses give idsnticaj wesults. In the massive case, however, we kusf distinguish between qyarks and antiquarks. Vor HISQ zuarks at zero spatial momentum, this corresponds do chkusiun $xp_0 = ia{m_{\text{tree}}}$ or $ap_0=-ia{m_{\text{tree}}}$ respectively. We dhpoxe outgoing qucrks or antiquarls — yre “scattering” or “anujhjlation” channels rfspectidely — to ensurt we co not attempt to compute vqnishing matgix wlements. Thus we cclculate the mattix elrments of $V_0$ and $A_k$ in tke scaftering chajnel and kw $A_0$ and $V_k$ in thd aknivilation channel. This procqdure is talid, even ag nomzero jattice spwcing, provided we match ho thg same channel ij continuum QCD.

Unfortunately, from the practicsl vivwpoint oy calcmlating Feynman diagrams, usiny massivz quarys complicztes thx numerical yntegration cmjsiderably. Tie chief qiffuculry lies kn the annihilstion chauuel, which contains a Coulomn sinehlarity that muwt ve handled with a ruberwcvion sgnction. We dhscurs gne suctraction fmncgionx employed in this wmrk jn more detail in Spiendix \[apk:subfn\]. Fuwthermore, in yhe automated pertkrbatmon thxory rputynes, the pole in the NRQCD prolagator cgosfes the contjur if integratipn and we can no longer carry out the uwual Wick rotation back to Minkowski xkace. Instead, xe muse deform dhe integration contiurs and introduct a triple contour to ehsure dhe shability of numerical integration [@hart07; @mueller11].

For both channels, the lattixe matciv elements just be mdteheb to thqir roitinuum QCD countvrparts. Analytic expressions for the relavcnt QCD contributions already ebist in the lkterature. References [@brazten95] and [@jones99] discuss the annihilation channel for the axial-vegtor currwnt, wgilst [@boylx00] present results fot both
for massless HISQ quarks with NRQCD formulated moving (mNRQCD) was for the vector [@mueller11]. massive quarks similar calculations using the lattice action for both quarks have carried out for Wilson quarks in [@kuramashi98] and for various implementations of NRQCD [@braaten95; @jones99; @boyle00; @hart07]. To our knowledge, no matching calculations with mixed actions massive quarks been in the literature. Moving from massless to massive relativistic quarks complicates the matching procedure. In the case, quarks and antiquarks at zero spatial momentum indistinguishable and consequently scattering annihilation processes give identical results. the case, however, must between and antiquarks. For quarks at zero spatial momentum, this corresponds to choosing $ap_0 = ia{m_{\text{tree}}}$ or $ap_0=-ia{m_{\text{tree}}}$ respectively. We choose quarks or the “scattering” “annihilation” respectively to ensure we attempt to compute vanishing matrix elements. the matrix elements of $V_0$ and $A_k$ in scattering channel of $A_0$ and $V_k$ in the channel. This procedure is valid, even at nonzero spacing, provided we match to the same channel in continuum QCD. Unfortunately, from the practical calculating Feynman diagrams, using quarks complicates the integration The difficulty in the channel, which contains a Coulomb singularity that must be handled with subtraction function. We discuss the subtraction functions employed in this more in Appendix \[app:subfn\]. in the automated perturbation routines, pole in the NRQCD the of can longer out the usual Wick back to Minkowski space. Instead, must deform the integration contour to ensure the stability of numerical integration @mueller11]. For both channels, the lattice matrix must be matched to their continuum QCD counterparts. Analytic expressions for the QCD contributions in the literature. References [@braaten95] and [@jones99] discuss annihilation channel for the current, whilst [@boyle00] present results for both
 for massless HISQ quarks with nRQCD formuLated In a MovInG fraMe (mNrQCD) was undertaKEn foR the vector and tensor heaVy-ligHt CUrreNTs In [@mueLler11].

For MAsSIVe qUaRkS siMiLAr MatchIng CalculaTions using The SaMe lattice actIOn For both quaRks Have been carrIed Out for wiLsoN QuarkS in [@KuramAshi98] anD For varIous impleMeNTationS Of NRQCD IN [@BrAateN95; @jones99; @boyle00; @hart07]. To OUr KNowledge, no matcHing caLcULaTIOns WitH mixed actiOnS and mASsive reLAtIVIStiC Quarks have beeN reported in THe lIteratUrE.

MoVIng froM massLeSS to Massive relaTiviStic quarkS complICates thE MatchinG proceDurE. In The fORmEr CasE, qUArkS AnD anTIquArks at zeRo SpAtial MomeNTUM Are iNdiStinGuishAble and conseqUenTly sCAttEring And anNihiLaTion pRocessEs givE iDentical results. in thE massive cAse, HoWevEr, We musT DistinGuiSh bEtween qUarks anD AntIqUARKs. for HISQ quarks at zerO sPATiAl momentUm, this COrReSPonds to cHoOsiNg $ap_0 = IA{M_{\text{Tree}}}$ OR $aP_0=-ia{m_{\text{Tree}}}$ reSPeCtIvely. We ChOose ouTgOinG quArks oR AntiQuarks — The “scattEring” OR “annihilation” cHAnnels respectIVeLY — To ENsurE we Do not attempT to cOMputE vanIShIng MAtrix ElemeNtS. thUS we calculate the matrIx ElemenTs of $V_0$ And $A_k$ in the scaTtering chaNNEL and of $A_0$ aNd $V_k$ IN tHE annihilation cHanneL. This proceDUre is valId, eveN at nonzeRo lattice SPAcing, proVidEd wE maTch TO ThE same channel iN COntiNuUm QCD.

UnForTunatelY, frOm tHe pRacTiCal viewpoInt of calCuLaTiNg feyNman dIAgrams, usInG maSsIve QuarkS CompliCates The nUmErICal IntegraTIoN COnsiDeRaBly. THe cHiEf difFicuLTy lIes in thE annihilaTioN ChanNeL, wHich conTains a Coulomb SiNgularity tHaT muSt be haNDLed with a Subtraction function. We diSCuss the SubTractIon fUnctions eMplOyed in ThiS Work in More deTail iN APpeNDIx \[app:SUBfN\]. FuRtHermore, in tHE AutOmateD pErtuRbation Theory routines, the pOLe iN the NRQCD propAgaTor cROSsEs tHE cONtoUr OF inTEGration and we can No longer caRrY OuT the usual WICk rOtAtion baCk to MinKowskI Space. InStead, we muSt deform tHe InteGRAtiOn contours And introDuce a tripLE contOUr To ensUre The staBiLitY of nuMericaL IntEgratIon [@harT07; @mUeller11].

for boTh Channels, The lattice matrix elementS must bE matcHed To their coNtiNUum qCD counteRparTs. Analytic ExpResSions For THe relEvanT qCd coNTribuTionS Already exISt In tHE LiTerature. RefERENceS [@braaTen95] ANd [@joneS99] disCuss the annihilatiON channel for the AxiaL-VEctOr cURrenT, wHilst [@boyle00] presEnt ReSULts for boTh
 for massless HISQ quarkswith NRQCD  form ula ted  i n amovi ng frame (mNRQ C D) w as undertaken for thevecto ra nd t e ns or he avy-lig h tc u rre nt sin[@ m ue ller1 1]. 

For m assive qua rks  s imilar match i ng  calculati ons  using the s ame  latti ce  ac t ion f orbothquarks have b een carri ed out fo r  Wilson q ua rksin [@kuramashi98] an d  for various i mpleme nt a ti o n s o f N RQCD in [@ br aaten 9 5; @jon e s9 9 ; @bo y le00; @hart07 ]. To our k n owl edge,no  ma t chingcalcu la t ion s with mixe d ac tions and  massi v e relat i visticquarks  ha vebeen re po rte di n t h elit e rat ure.

Mo vi ng  from  mas s l e s s to  ma ssiv e rel ativistic qua rks  com p lic atesthe m atch in g pro cedure . Inth e former case,quar ks and an tiq ua rks  a t zer o  spati almom entum a re indi s tin gu i s h ab le and consequentl ys c at tering a nd ann i hi la t ion proc es ses  giv e ident ical re sults. I n them as si ve case ,howeve r,  we  mu st di s ting uish b etween q uarks and antiquarks .  For HISQ qua r ks a tz erospa tial moment um,t hiscorr e sp ond s  to c hoosi ng $a p _0 = ia{m_{\text{tr ee }}}$ o r $ap _0=-ia{m_{\te xt{tree}}} $ r espectiv ely. We choose outgoin g qua rks or ant i quarks —  the“scatter ing” or “ a n nihilati on”  ch ann els r es pectively — t o ensu re  we donot  attemp t t o c omp ute  v anishingmatrix e le me nt s.  Th us we calculat ethe  m atr ix el e mentsof $V _0$an d$ A_k $ in th e  s c a tter in gchan nel  a nd of  $A_ 0 $ a nd $V_k $ in theann i hila ti on  channe l. This proce du re is vali d,  ev en atn o nzero la ttice spacing, provided we matc h t o the  sam e channel  in  conti nuu m  QCD.
Unfor tunat el y,f r om th e pr act ic al viewpoi n t  of  calc ul atin g Feynm an diagrams, using mas sive quarks c omp lica t e sthe nu m eri ca l  in t e gration conside rably. The  c h ie f difficul t y l ie s in th e annih ilati o n chann el, which  contains  a  Cou l o mbsingularit y that m ust be ha n dledw it h a s ubt ractio nfun ction . We d i scu ss th e subt ra ctionfunct io ns emplo yed in this work in mor e deta il in  Ap pendix \[ app : sub fn\]. Fur ther more, in t heaut omate d p e rturb atio n  t heo r y rou tine s , the pol e  i n t h e  N RQCD propag a t o r c rosse s t h e cont ourof integration an d  we can no lon gerc a rry  ou t  the  u sual Wick rota tio nb a ck to Mi nk owski space . Instea d, we mu st def orm th e integ r a ti o n cont ours  an d introdu ceat riple c on to u r to e nsur ethe st abilit y  ofn u merical integrat ion [ @ h art07 ;  @m uelle r1 1].

Fo r  bot h channels , the latti ce mat rixeleme nts mus tbe mat che dto their c o ntinuum Q CD co unterpa rt s. A nal ytic e xpre s s ionsforth e r elevant Q C D  c o nt ri b uti onsalrea dy  exi st in the literatu re. Referen ce s [ @ b raaten 9 5] a nd [@jones 99]  disc u s s the anni h ilat i on chann el for  the a xial-ve c tor  c urrent,  wh i l st [@boyl e00] pres e ntre sult s for bo th 
 for_massless HISQ_quarks with NRQCD formulated_in a_moving_frame (mNRQCD)_was_undertaken for the_vector and tensor_heavy-light currents in [@mueller11].

For_massive quarks similar_matching_calculations using the same lattice action for both quarks have been carried out for_Wilson_quarks in_[@kuramashi98]_and_for various implementations of NRQCD_in [@braaten95; @jones99; @boyle00; @hart07]._To our_knowledge, no matching calculations with mixed actions and_massive_relativistic quarks have_been reported in the literature.

Moving from massless to massive_relativistic quarks complicates the matching procedure._In the former_case,_quarks_and antiquarks at zero_spatial momentum are indistinguishable and consequently_scattering and annihilation processes give identical_results. In the massive case, however, we_must distinguish between quarks and antiquarks._For HISQ quarks at zero_spatial momentum,_this corresponds to choosing $ap_0_= ia{m_{\text{tree}}}$ or_$ap_0=-ia{m_{\text{tree}}}$ respectively._We choose outgoing_quarks or antiquarks — the “scattering”_or “annihilation” channels_respectively — to ensure we do_not_attempt to compute_vanishing_matrix_elements. Thus_we calculate the_matrix_elements of_$V_0$_and $A_k$ in the scattering channel_and_of $A_0$ and $V_k$ in the annihilation_channel. This procedure is_valid,_even at nonzero lattice_spacing, provided we match to_the same channel in continuum QCD.

Unfortunately,_from the_practical viewpoint_of calculating Feynman diagrams, using massive quarks complicates the numerical integration_considerably. The chief difficulty lies in_the annihilation channel, which_contains a_Coulomb_singularity that must_be_handled with_a subtraction function. We discuss the subtraction_functions employed_in this work in more detail_in Appendix \[app:subfn\]. Furthermore,_in_the automated perturbation theory routines, the_pole in the NRQCD propagator crosses_the contour of integration and_we_can_no longer carry out the_usual Wick rotation back to Minkowski_space. Instead, we_must deform the integration contours and introduce_a_triple contour to ensure the stability_of_numerical integration [@hart07; @mueller11].

For both channels,_the_lattice_matrix elements must be matched_to their continuum QCD counterparts. Analytic_expressions for the relevant QCD contributions already exist in_the literature. References_[@braaten95] and [@jones99] discuss the_annihilation_channel_for the axial-vector current, whilst [@boyle00] present results for both