2014lla] these MI can be expressed in terms of Goncharov polylogarithms. The complexity of the resulting expressions is certainly promising that the project of computing all MI relevant to massless QCD, namely all eight-denominator MI with arbitrary configuration of the external momenta, is feasible. Having such a complete library of two-loop MI, the analog of $A_0, B_0, C_0, D_0$ scalar integrals at one loop, the reduction of an arbitrary two-loop amplitude [à la OPP]{} can pave the road for a NNLO automation in the near future.

As experience shows, there are several issues that will need to attract our attention in order to accomplish our goal. First of all in order to systematize the whole procedure of reducing an arbitrary Feynman Integral in terms of MI in an efficient way, a deepening of our current understanding of IBP identities [@Baikov:1996iu; @Smirnov:2003kc] is necessary [@Ita:2015tya]. Secondly, further standardising the procedure to obtain a canonical form of DE [@Henn:2013pwa], which drastically simplifies the expression of MI in terms of GPs, is certainly a very desirable feature. Thirdly, the inclusion of MI with massive internal propagators, at a first stage with one mass scale corresponding to the heavy top quark, will provide the complete basis for NNLO QCD automated computations. Moreover, the calculation of boundary terms for the DE can benefit from further developments and exploitations of expansion-by-regions techniques, in conjunction with Mellin-Barnes representation of the resulting integrals. Finally, on the numerical side, a more efficient computation of polylogarithms is also necessary.

Acknowledgements {#acknowledgements.unnumbered}
----------------

This research was supported by the Research Funding Program ARISTEIA, HOCTools (co-financed by the European Union (European Social Fund ESF) and Greek national funds through the Operational Program “Education and Lifelong Learning” of the National Strategic Reference Framework (NSRF)). We gratefully acknowledge fruitful discussions with J. Henn, A. V. Smirnov, V. A. Smirnov, M. Czakon, S. Weinzierl, C. Duhr and E. Panzer during various stages of this project. We would like also to thank H. Frell
2014lla ] these MI can be expressed in terms of Goncharov polylogarithms. The complexity of the resulting expression is surely promising that the project of calculate all MI relevant to massless QCD, namely all eight - denominator MI with arbitrary shape of the external momenta, is feasible. Having such a arrant library of two - closed circuit MI, the analogue of $ A_0, B_0, C_0, D_0 $ scalar integral at one loop, the reduction of an arbitrary two - closed circuit amplitude [ à la OPP ] { } can pave the road for a NNLO automation in the near future. 

 As experience shows, there are respective issues that will need to attract our care in order to carry through our finish. First of all in order to systematize the whole routine of reducing an arbitrary Feynman Integral in terms of michigan in an efficient manner, a deepening of our current understanding of IBP identities   [ @Baikov:1996iu; @Smirnov:2003kc ] is necessary   [ @Ita:2015tya ]. Secondly, further standardising the operation to obtain a canonical form of DE   [ @Henn:2013pwa ], which drastically simplifies the expression of MI in terms of GPs, is certainly a very desirable feature. Thirdly, the inclusion of MI with massive internal propagators, at a first stage with one mass plate match to the heavy circus tent quark, will provide the complete basis for NNLO QCD automatize computations. furthermore, the calculation of boundary terms for the DE can benefit from further developments and exploitation of expansion - by - regions techniques, in conjunction with Mellin - Barnes representation of the resulting integrals. Finally, on the numerical side, a more efficient computation of polylogarithms is besides necessary. 

 Acknowledgements { # acknowledgements.unnumbered } 
 ---------------- 

 This inquiry was supported by the Research Funding Program ARISTEIA, HOCTools (co - financed by the European Union (European Social Fund ESF) and Greek national funds through the Operational Program “ Education and Lifelong Learning ” of the National Strategic Reference Framework (NSRF) ). We gratefully acknowledge fruitful discussion with J.   Henn, A.   V.   Smirnov, V.   A.   Smirnov, M.   Czakon, S.   Weinzierl, C.   Duhr and E.   Panzer during various stages of this project. We would wish also to thank H.   Frell
2014lla] these MI can be expresstd in terms of Gouxharov polylkgarithmr. The complexity of the resuptung ezpressions is certainlh promisijg that rhe kroject of computmhg all MI relsyant co massless QCD, kamely all aight-denominatmr MN with arbitrary configuration of thq extermap momenta, is fgasibkq. Habpnn such a complete library of two-moop MI, the analog og $A_0, B_0, C_0, D_0$ scalar integrals at lne loop, the reduchion of an qrbieeary two-loop amplitude [à la OPP]{} can pave the road for a NNLO automxtion in the neqr fuhore.

As experixnce srows, there avv severdl issurs that will nced tm artract our attention mn order to accomplifh our godl. First of all in ireer tm syvtemxrizd tge wgole pgocxdure of resucing an aebitrary Feynman Inuegwql in terms or MI ig wn efficient way, a deepening of our curgent understanding of IBP ieentities [@Baikov:1996iu; @Smitnov:2003kc] is gecessary [@Ita:2015tya]. Secondly, further standardising tha proreaurt bj ibhain a canonical form of DE [@Henn:2013pwa], which drasejcslky simplifies bhe expression of KI im terms of GPs, is cexfajnly a very desirahle feajure. Tyirdly, tht inckusion of MI with massive ibternal propcgarors, at a first stcge with one masx scake corresponding to the heaby top quarn, will prkxide the completd bsshs for NNLO QCD automated somputatiins. Koreovef, thg calcujation of houndary terms for the DF can banefit frol further developments and exploitations of ecpdnspon-by-regijns tcchniques, in cogjunction with Mellin-Yarnes representztion oh the resultyng integrals. Vinally, on tie numerisal wide, a more dfficient compitation oy poltlogarithms is alsp ngcsssary.

Acknowledytmebts {#acknowledgekengs.ugnlmbxred}
----------------

Trhs research fas rupoprted by the Rescarzh Finding Program ARISTAIA, GOCTools (co-financec ny the Eutopean Unyon (European Xocial Fund ESF) anf Grexk natmonal gunqs through the Operational Profram “Educwtijn and Lifeljng Oearning” of che National Strategic Reference Framewock (NSRF)). We gratefully axknowledge fruitful ciscussions with J. Henn, A. V. Vmirnov, V. A. Smirnov, M. Czakon, S. Weinzierl, G. Duhr and E. Panzer durihg varhous dtages of this project. We would like also to thank H. Frell
2014lla] these MI can be expressed in Goncharov The complexity the resulting expressions project computing all MI to massless QCD, all eight-denominator MI with arbitrary configuration the external momenta, is feasible. Having such a complete library of two-loop MI, analog of $A_0, B_0, C_0, D_0$ scalar integrals at one loop, the reduction an two-loop [à OPP]{} can pave the road for a NNLO automation in the near future. As experience shows, are several issues that will need to attract attention in order to our goal. First of all order systematize the procedure reducing arbitrary Feynman Integral terms of MI in an efficient way, a deepening of our current understanding of IBP identities [@Baikov:1996iu; is necessary further standardising procedure obtain canonical form of which drastically simplifies the expression of of GPs, is certainly a very desirable feature. the inclusion MI with massive internal propagators, at first stage with one mass scale corresponding to heavy top quark, will provide the complete basis for NNLO QCD automated computations. Moreover, the boundary terms for the can benefit from developments exploitations expansion-by-regions in conjunction Mellin-Barnes representation of the resulting integrals. Finally, on the numerical side, more efficient computation of polylogarithms is also necessary. Acknowledgements {#acknowledgements.unnumbered} research supported by the Funding Program ARISTEIA, HOCTools by European Union (European Social and national Operational “Education Lifelong Learning” of the Strategic Reference Framework (NSRF)). We acknowledge fruitful discussions with V. A. Smirnov, M. Czakon, S. Weinzierl, C. and E. Panzer during various stages of project. We would like also to thank H. Frell
2014lla] these MI can be expressed iN terms of GoNcharOv pOlyLoGariThms. the complexity oF The rEsulting expressions is cErtaiNlY PromISiNg thaT the proJEcT OF coMpUtIng AlL mI RelevAnt To masslEss QCD, nameLy aLl Eight-denominAToR MI with arbItrAry configuraTioN of the ExTerNAl momEntA, is feAsible. hAving sUch a complEtE LibrarY Of two-loOP mI, The aNalog of $A_0, B_0, C_0, D_0$ scalar INtEGrals at one loop, The redUcTIoN OF an ArbItrary two-lOoP amplITude [à la opP]{} CAN PavE The road for a NNlO automatioN In tHe near FuTurE.

as expeRiencE sHOws, There are sevEral Issues thaT will nEEd to attRAct our aTtentiOn iN orDer tO AcCoMplIsH Our GOaL. FiRSt oF all in orDeR tO systEmatIZE THe whOle ProcEdure Of reducing an aRbiTrarY feyNman INtegrAl in TeRms of mI in an EfficIeNt way, a deepening Of ouR current uNdeRsTanDiNg of Ibp identItiEs [@BAikov:1996iu; @smirnov:2003KC] is NeCESSaRy [@Ita:2015tya]. Secondly, fuRtHER sTandardiSing thE PrOcEDure to obTaIn a CanoNICal foRm of de [@HEnn:2013pwa], whIch draSTiCaLly simpLiFies thE eXprEssIon of mi in tErms of gPs, is cerTainlY A very desirable FEature. Thirdly, THe INClUSion Of Mi with massivE intERnal PropAGaTorS, At a fiRst stAgE WiTH one mass scale corresPoNding tO the hEavy top quark, wIll provide THE Complete BasiS FoR nNLO QCD automatEd comPutations. MOReover, thE calcUlation oF boundary TERms for thE DE Can BenEfiT FRoM further develOPMentS aNd exploItaTions of ExpAnsIon-By-rEgIons technIques, in cOnJuNcTiOn wIth MeLLin-BarneS rEprEsEntAtion OF the reSultiNg inTeGrALs. FInally, oN ThE NUmerIcAl Side, A moRe EfficIent COmpUtation Of polylogAriTHms iS aLsO necessAry.

AcknowledgEmEnts {#acknowLeDgeMents.uNNUmbered}
----------------

THis research was supported BY the ResEarCh FunDing program ARiSTeIA, HOCtooLS (co-finAnced bY the EUrOpeAN union (eURoPeaN SOcial Fund Esf) And greek NaTionAl funds Through the OperatioNAl PRogram “EducatiOn aNd LiFELoNg LEArNIng” Of THe NATIonal Strategic REference FrAmEWoRk (NSRF)). We grATefUlLy acknoWledge fRuitfUL discusSions with j. Henn, A. V. SmIrNov, V. a. sMirNov, M. Czakon, s. WeinzieRl, C. Duhr anD e. PanzER dUring VarIous stAgEs oF this ProjecT. we wOuld lIke alsO tO thank h. FrelL
2014lla] these MI can be e xpressed i n ter msofGo ncha rovpolylogarithms .  The  complexity of the res ultin ge xpre s si ons i s certa i nl y pro mi si ngth a tthe p roj ect ofcomputingall  M I relevant t o  m assless QC D,namely all e igh t-deno mi nat o r MIwit h arb itrary config uration o ft he ext e rnal mo m e nt a, i s feasible. Havin g  s u ch a completelibrar yo ft w o-l oop  MI, the a na log o f  $A_0,B _0 , C _0, D_0$ scalar i ntegrals at one  loop,  t her educti on of  a n  ar bitrary two -loo p amplitu de [àl a OPP]{ }  can pa ve the  ro adfora  N NL O a ut o mat i on  in the  near fu tu re .

As  exp e r i e ncesho ws,there  are severaliss uest hat  will  need  toat tract  our a ttent io n in order to a ccom plish our  go al . F ir st of all in  or der  to sys tematiz e  th ew h o le  procedure of redu ci n g  a n arbitr ary Fe y nm an Integral  i n t erms o f MIin a n  e fficient  way,a  d ee peningof  our c ur ren t u nders t andi ng ofIBP iden titie s  [@Baikov:1996 i u; @Smirnov:2 0 03 k c ]i s ne ces sary [@Ita: 2015 t ya].  Sec o nd ly, furth er st an d ar d ising the procedure  t o obta in acanonical for m of DE [@ H e n n:2013pw a],w hi c h drasticallysimpl ifies thee xpressio n ofMI in te rms of GP s ,  is cert ain lya v ery d es irable featur e .  Thi rd ly, the  in clusion  of  MI  wi thma ssive int ernal pr op ag at or s,at af irst sta ge  wi th  on e mas s  scale  corr espo nd in g  to  the he a vy t op q ua rk , wi llpr ovide  the com plete b asis forNNL O  QCD  a ut omatedcomputations.  M oreover, t he  ca lculat i o n of bou ndary terms for the DEc an bene fit  from  fur ther deve lop mentsand exploi tation s ofex pan s i on-by - r eg ion stechniques , inconju nc tion  with M ellin-Barnes repre s ent ation of theres ulti n g  i nte g ra l s.Fi n all y ,  on the numeric al side, a  m o re  efficient com pu tationof poly logar i thms is  also nec essary.
Ac know l e dge ments {#ac knowledg ements.un n umber e d} 
---- --- ------ -- -
Thisresear c h w as su pporte dby the  Rese ar ch Fundi ng Program ARISTEIA, HO CTools  (co- fin anced bythe Eur opean Uni on ( European S oci alFundESF )  andGree k  n ati o nal f unds through t h eOpe r a ti onal Progra m “ Edu catio n a n d Life long  Learning” of the National Strat egic R efe ren c e Fr am ework (NSRF)).  We  g r a tefullyac knowledge f ruitfuldi s cussi ons wi th J.Henn, A . V. Smirno v, V . A . Smirnov , M .C zakon,S.  W e inzier l, C .Duhr a nd E.P anze r during various s tages o f thi s  pr oject .We woul d  lik e also tothank H. Fr ell
2014lla] these_MI can_be expressed in terms_of Goncharov_polylogarithms._The complexity_of_the resulting expressions_is certainly promising_that the project of_computing all MI_relevant_to massless QCD, namely all eight-denominator MI with arbitrary configuration of the external momenta,_is_feasible. Having_such_a_complete library of two-loop MI,_the analog of $A_0, B_0,_C_0, D_0$_scalar integrals at one loop, the reduction of_an_arbitrary two-loop amplitude_[à la OPP]{} can pave the road for a_NNLO automation in the near future.

As_experience shows, there_are_several_issues that will need_to attract our attention in order_to accomplish our goal. First of_all in order to systematize the whole_procedure of reducing an arbitrary Feynman_Integral in terms of MI_in an_efficient way, a deepening of_our current understanding_of IBP_identities [@Baikov:1996iu; @Smirnov:2003kc] is_necessary [@Ita:2015tya]. Secondly, further standardising the procedure_to obtain a_canonical form of DE [@Henn:2013pwa], which drastically_simplifies_the expression of_MI_in_terms of_GPs, is certainly_a_very desirable_feature._Thirdly, the inclusion of MI with_massive_internal propagators, at a first stage with_one mass scale corresponding_to_the heavy top quark,_will provide the complete basis_for NNLO QCD automated computations. Moreover,_the calculation_of boundary_terms for the DE can benefit from further developments and exploitations_of expansion-by-regions techniques, in conjunction with_Mellin-Barnes representation of the_resulting integrals._Finally,_on the numerical_side,_a more_efficient computation of polylogarithms is also necessary.

Acknowledgements_{#acknowledgements.unnumbered}
----------------

This research_was supported by the Research Funding_Program ARISTEIA, HOCTools (co-financed_by_the European Union (European Social Fund_ESF) and Greek national funds through_the Operational Program “Education and_Lifelong_Learning”_of the National Strategic Reference_Framework (NSRF)). We gratefully acknowledge fruitful_discussions with J. Henn,_A. V. Smirnov, V. A. Smirnov, M. Czakon, S. Weinzierl, C. Duhr and E. Panzer_during_various stages of this project. We_would_like also to thank H. Frell