$LoadFeynArts=True;
<<FeynCalc`
$FAVerbose=False;

top=CreateTopologies[1,1->2,ExcludeTopologies->{WFCorrections,SelfEnergies,Tadpoles}];
diags=InsertFields[top,{V[5]}->{V[5],V[5]},Model->"SMQCD",InsertionLevel -> {Particles},ExcludeParticles->{F[__],S[_],V[_],U[1|2|3|4]}];
Paint[diags, ColumnsXRows -> {2, 1},SheetHeader -> False,SheetHeader->None,Numbering -> None,ImageSize->{512,256}];

amps=FCFAConvert[CreateFeynAmp[diags, Truncated -> True,GaugeRules->{},PreFactor->1],IncomingMomenta->{p1},
OutgoingMomenta->{p2,p3},LoopMomenta->{q},DropSumOver->True,UndoChiralSplittings->True,ChangeDimension->D,List->True,SMP->True,LorentzIndexNames->{λ,μ,ν},
SUNIndexNames->{a,b,c}]//Contract//SUNSimplify

amps2=amps[[2]]/.p3->-p2/.p2->p

amps1=amps[[1]]/.p3->-p2/.p2->p

aaa=amps1+amps2;
fcaa=OneLoopSimplify[aaa,q];
AA=OneLoop[q,fcaa]

Simplify[AA//.{FCI[C0[0,SP[p,p],SP[p,p],0,0,0]]->-(((-3+D) B0[FCI[SP[p,p]],0,0])/SP[p,p]),FCI[B0[FCI[SP[p,p]],0,0]]->1/(16*Epsilon*Pi^4)}];
Series[FCI[%]/.D->4-2 Epsilon,{Epsilon,0,0}]//Normal//Collect2[#,Epsilon]&
SelectNotFree[%,1/Epsilon]

ints=FCI[{FCI[B0[SP[p,p],0,0]]->1/(16*Epsilon*Pi^4),FCI[C0[0,SP[p,p],SP[p,p],0,0,0]]->-1/(16*Epsilon*Pi^4*SP[p,p])}]
Series[AA/.ints/.D->4-2 Epsilon,{Epsilon,0,0}]//Normal//Collect2[#,Epsilon]&
SelectNotFree[%,1/Epsilon]