Towards a Heavy Diquark Effective Theory for Weak Decays of Doubly Heavy Baryons

Yu-Ji Shi 1,2 ∗, Wei Wang 2 †, Zhen-Xing Zhao 2,3 ‡, Ulf-G. Meißner 1,4,5 § 1 Helmholtz-Institut für Strahlenund Kernphysik and Bethe Center for Theoretical Physics, Universität Bonn, 53115 Bonn, Germany 2 INPAC, SKLPPC, MOE KLPPC, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China 3 School of Physical Science and Technology, Inner Mongolia University, Hohhot 010021, China 4 Institute for Advanced Simulation, Institut für Kernphysik and Jülich Center for Hadron Physics, Forschungszentrum Jülich, D-52425 Jülich, Germany 5 Tbilisi State University, 0186 Tbilisi, Georgia


I. INTRODUCTION
In the past, the conventionl quark model has successfully explained structures of numerous hadronic states observed in a large number of experiments. However, not all predicted particles by the quark model have been experimentally established. In particular, doubly heavy baryons, that is baryonic states made of two heavy quarks, are of this type. After pursuing the Ξ cc for many years, the LHCb collaboration finally announced in 2017 the observation of Ξ ++ cc , a lowest-lying doubly-charmed baryon whose mass is give as [1] m Ξ ++ cc = (3621.40 ± 0.72 ± 0.27 ± 0.14) MeV.
This inspiring observation follows an earlier prediction Ref. [2], where the Ξ ++ cc is expected to be reconstructed from the decay channel Ξ ++ cc → Λ + c K − π + π + . One year later, LHCb has also successfully measured the Ξ ++ cc 's lifetime [3], and reconstructed this resonance from the Ξ + c π + final state [4]. Thus, the existence of the Ξ ++ cc is unambiguously established. We believe that through continuous experimental efforts [5][6][7], other heavier doubly heavy baryons could be discovered in the future. In addition, there have been numerous theoretical studies aiming to understand the dynamical and spectroscopical properties of the doubly-heavy baryon states, see e.g. Refs. .
However, a comprehensive description of the decay mechanism of doubly heavy baryons is not established yet.
Generally, an ideal platform for studying hadrons is through semi-leptonic weak decays. The main advantage of a semi-leptonic process is its naturalness in separating the QCD relevant and the QCD irrelevant dynamics in the weak decays. All the QCD dynamics is encapsulated in the hadron transition matrix element, which is independent from the leptonic part and can be parametrized by several form factors. However, as a three-body system, a doubly heavy baryon posseses a much more complicated dynamics than a heavy meson.
A straightforward way to consider this problem is to reduce a doubly heavy baryon into a two-body system, where two of the three quarks are treated as a point-like diquarks. Generally, each two quarks in a baryon form a color antitriplet so that they might be bound by an attractive potential. However, for a doubly heavy baryon, it is more reasonable to treat the two heavy quarks behave as static color sources and thus as a diquark, see e.g. [35,36]. The effective distance between the two heavy quarks can be estimated to be r QQ ∼ 1/m Q v, where v is the four-velocity of the baryon. Further, the distance between one of the heavy quarks and the light quark is approximately r Qq ∼ 1/Λ QCD . Since r QQ /r Qq ∼ 1/m Q ≪ 1, the two heavy quarks can be combined to be a pointlike diquark. In the heavy diquark limit, the heavy diquark system can be treated as a static color source in the3 representation, just like a heavy anti-quark. Some earlier papers [38][39][40][41] have used the heavy quark-diquark symmetry to simplify the transition form factors.
In this work, we will try to develop a heavy diquark effective theory (HDiET), whose Lagrangian is expanded in powers of r QQ /r Qq . At leading-order, the diquark appears as a point-like scalar or axial-vector particle described by a scalar or axial-vector field in the color3 representation.
The scalar HDiET has been developed in [42], where the leading order (LO) effective Lagrangian coupling two scalar diquarks and two light quarks was obtained. In this work, we will first construct HDiET for both scalar and axial-vector diquarks. For the transition form factors we will assume the applicability of HDiET, and by assuming the diquark to be a point-like particle, we can construct the weak and electromagnetic transition currents of the diquarks according to the SU(2) heavy flavor symmetry and U(1) symmetry. On the other hand, in the large recoil region, the diquark currents will be derived through the matching between QCD and HDiET at tree level. We then show that the six transition form factors of doubly heavy baryon semi-leptonic decay can be reduced into only one soft function. Furthermore, it will be shown that this soft function is an universal quantity which is nothing but the well known Isgur-Wise function in HQET for heavy meson decays. These results can be used in the phenomenology studies.
This article is organized as follows: In section II, we construct the LO diquark effective theory (DiET) Lagrangian including the kinetic part as well as the terms coupling with weak and electromagnetic fields. The DiET is also transformed to HDiET in the heavy diquark limit. In section III, we derive the diquark transition currents both from symmetry and tree level matching. Section IV focuses on the semi-leptonic decays of doubly heavy baryons. We perform a reduction of the transition matrix element, where a universal soft function is factorized out and the q 2 distributions of all the six form factors are completely determined from it. The resulting form factors are used to predict the semi-leptonic decay widths. Section V contains our conclusions.

A. Effective Lagrangian for Scalar and Axial-vector Diquark
In this section we will construct the DiET at leading order. The first step is to write down the diquark effective Lagrangian. We denote the scalar and axial-vector diquark field as S i and X i µ , where i is the3 color index. The free scalar diquark Lagrangian is simply Here we have assumed that both the scalar and axial-vector diquark have the same mass m X . On the other hand, to construct the axial-vector diquark Lagrangian, one should be aware of that X i µ is a matter field in the color fundamental representation3, instead of the adjoint representation which belongs to the standard gauge fields. Therefore, the axial-vector diquark field is not required to couple with any conserved current, and it seems not necessary to construct the effective Lagrangian with the building blocks of the strength tensor F i µν = ∂ µ X i ν −∂ ν X i µ as is done for Yang-Mills theory. Instead, one can write down a general form However, note that X µ has four components while a spin-1 particle has only three physical degrees of freedom. According to the canonical theory, one needs to introduce two second-class constraints for the Hamiltonian to remove one redundant canonical variable as well as its conjugate momentum.
As a result, one still arrives at a gauge-field-like Lagrangian with an on-shell constraint condition ∂ µ X iµ = 0.
Since the diquark is composed of two flavored heavy quarks, it is natural to dress the diquark fields with certain representation in the flavor space. Notice that in QCD, heavy quarks include bottom and charm. If we approximately assume m b ∼ m c → ∞, the mass matrix for (b, c) T is almost diagonal so that there exists a flavor SU(2) symmetry for the heavy quark sector of the QCD Lagrangian. Furthermore, in HQET, the leading power LagrangianQ v iv · DQ v is exactly invariant under the flavor SU(2) transformation. Such a transformation on a multiplet Q = (b, c) T is denoted as Q = (b, c) T , Q → U Q, U ∈ SU(2). Besides the SU(2) flavor symmetry, there is also a U(1) symmetry which corresponds to the electromagnetic (EM) interaction, As an effective theory of QCD, DiET should also reflect the SU(2) × U(1) symmetry. In the flavor space, a diquark field can be considered to have the structure q i q j , where i, j = b or c are flavor indexes. Thus a diquark field should be represented by a 2 × 2 matrix Note that the representation for a scalar diquark is anti-symmetric while the representation for an axial-vector diquark is symmetric. Under SU(2) × U(1), they transform as With these matrixes as basic building blocks, one can construct a SU(2) × U(1) invariant diquark Lagrangian. An efficient way to realize these symmetries is to apply the spinor representation for the diquark fields. Following Ref. [43], one firstly combines the spin-1 and spin-0 diquarks to be a multiplet, which is described by a bilinear spinor field Σ Σ = (X µ γ µ + Sγ 5 )C, where C is the charge conjugating matrix. A reason to choose such form is due to the Lorentz covariance. Under a general Lorentz transformation X µ → Λ µ ν X ν , one can show that Σ does transform in the expected manner, Σ → Λ 1/2 Σ Λ T 1/2 . In addition, in momentum space the equation of motion of the two constituent heavy quarks yakes the form Note that since the diquark is treated as a point-like particle, both the two constituent heavy quarks and the diquark itself share a common velocity v d , so that it is reasonable to operate with the same slash / v d on the both sides of Σ. Therefore, we can define Σ ′ To obtain the second equality we have used the on-shell constraint v d ·X(v d ) = 0. After transforming Σ ′ (v d ) into coordinates space, we can define a multiplet field K(x) as whereΣ(x) = γ 0 Σ † (x)γ 0 . According to Eq. (7), under SU(2) × U(1) transformation, K andK transform in the same manner as S, X µ and S † , X µ † . Therefore the kinematic Lagrangian of DiET is just the simplest globally SU(2) × U(1) invariant Lagrangian constructed by K,K, m X and one derivative operator where the trace acts in both flavor and spinor spaces. After expressing this equation in terms of X µ and S, the kinematic Lagrangian takes the form of a combination of a spin-1 part and a spin-0 where Tr f only acts in the flavor space. Compared with Eq. (4), this equation has no ∂ µ X † ν ∂ ν X µ term. The reason is that in the heavy quark limit, the diquark field is a very massive field, which is approximately on shell and satisfies the constraint ∂ µ X µ = 0.
Next, let us consider how the diquark field couples to external sources. At the quark level, the weak and the EM coupling come from the coupling terms in QCD Here, λ 1 and λ 2 are two independent coupling constants. After being expressed in terms of X µ and S, the coupling Lagrangians of the X-J-X, S-J-X, X-J-S and S-J-S types are given by where Tr f only acts in the flavor space and J µ = V µ + A µ .F µν = 1 2 ǫ µναβ F αβ is the dual field strength tensor. We have also defined two kinds of commutators in the flavor space

B. Heavy Diquark Effective Theory (HDiET)
A diquark in the color3 representation interacts with gluons in a similar way as a an antiquark. Replacing the ordinary derivatives in Eq. (2) and Eq. (4) with covariant derivatives, one can introduce the coupling of a diquark and a gluon where a , g d is the effective coupling constant between the diquark and the gluon.
In the heavy diquark limit, to expand the Lagrangian in power of 1/m 2 X , one has to separate the diquark field into a static part and a residual part as is done in with the heavy quark in HQET.
For the case of scalar diquark, the 1/m 2 X expansion is trivial. By factorizing out an exponential Note that each covariant derivative scales as Λ QCD . Thus in the heavy diquark limit, the second term in Eq. (23) is suppressed by Λ QCD /m X compared with the first term. Furthermore, at the leading order, S v is massless and its propagator is simply In case of an axial-vector diquark, just factorizing out an exponential phase is not enough. In the heavy diquark limit, one has to separate X µ into a static part exp Also note that both X µ v and Y µ v are dominated by the small momentum k ∼ Λ QCD . Let us introduce two projection operators P µ ν and T µ ν , Using the projection operators, one can project out the static part X µ v and the residual part Y µ v of the heavy axial-vector diquark field X µ which satisfy v · X v = 0 and (· · · X µ v ) † (· · · Y vµ ) = (· · · Y µ v ) † (· · · X vµ ) = 0, where the dots represent any possible insertion of covariant derivatives. Then the full diquark field can be separated as Inserting Eq. (27) into Eq. (22), and using integration by part ← − D = −D to make all the covariant derivatives act on the X, Y fields instead of the X † , Y † fields, one finally arrives at From the Lagrangian Eq. (28), one finds that X µ v is a massless field, while Y µ v is massive due to the non-diagonal mass term −(m 2 To obtain an effective theory containing only the massless field X µ v , one needs to integrate out the heavy degree of freedom Y µ v . One way to realize this is to use the saddle point approximation, where one first solves the equation of motion of the It is not simple to solve this matrix equation directly. To simplify it, we can multiply with v µ on both sides of the equation and introduce a power counting scheme to solve this equation perturbatively. Note that each covariant derivative D scales as Λ QCD which is small compared to m X . So by counting the number After inserting this solution of Y µ v back to Eq. (28), one finally obtains the effective Lagrangian in the form of a power expansion whereḠ µν = G a µνt a is the gluon tensor. In the Eq. (33), the second term represents the heavy diquark kinetic energy while the third term corresponds to the chromomagnetic coupling. These two terms are consistent with those given in Ref. [44][45][46], where a non-relativistic approach is used.
The propagator of the massless heavy axial-vector diquark is The heavy diquark can only couple to soft gluons. Through the following field redifinition, one can decouple the diquark field from gluon field: Using the decoupling transformation, one can replace all the covariant derivatives in Eq. (33) by ordinary derivatives, while the X field should be replaced by the dressed fieldX.

A. Diquark Transition Currents from Symmetry
When using DiET to study doubly heavy baryon decays B bQ → B cQ lν, for instance when the bb diquark turns into the bc diquark through the V − A currentcγ µ (1 − γ 5 )b, or electromagnetic if we approximate the diquark as a point like particle, we require the four most general kinds of diquark currents which correspond to pure axial-vector, axial-vector to scalar, scalar to axial-vector and pure scalar transitions. Note that Γ αβ µ , Γ β µ and Γ µ depend on the momentum of the initial and final diquarks. In the heavy diquark limit we can simply replace the ← − ∂ , ∂ with the four-velocities of the final and initial baryons iv 2 , −iv 1 , with w = v 1 · v 2 close to 1 for the low recoil region.
Consider first the case of V − A weak currentcγ µ (1 − γ 5 )b. According to Eq. (13), it is just a current coupling to the external source V 1 µ + iV 2 µ , which can be found from the expansion Straightforwardly, one can conclude that thecγ µ (1 − γ 5 )b current can be produced by operating with a derivative on the part of the Lagrangian of QCD that contains the couplings to the external On the other hand, on the diquark level, if one performs the same derivative operation on the DiET Lagrangian Eq. (15)(16)(17)(18)(19), one arrives at the V − A currents in the DiET form Explicitly for X → X, X → S and S → S transitions, one has Note that the antisymmetric S has only one non-vanishing component S bc , for flavor changing processes b → c there is no S → S transition. Similarly, the electromagnetic currents I Transition µ can be derived by acting with a derivative on A µ em , where C X is the total electric charge of X. It should be mentioned that all the currents in Eqs. (41)(42)(43)(45)(46)(47)(48) are expressed by the full diquark fields. These expressions are simpler in the heavy diquark limit. According to Eq. (27) and Eq. (32), the full diquark fields X, S are related with the effective ones X v , S v in HDiET as Inserting Eq. (49) into Eq. (41)(42)(43)(45)(46)(47)(48), at leading order, all the derivative operators are simply replaced by the corresponding four velocities where Λ = (λ 1 + λ 2 )m X . Similarly one can obtain the currents at next-to-leading order if the second expansion term of X µ in Eq. (49) is used, but the results will not be shown explicitly here.
It should be mentioned that like the chromomagnetic coupling in the Eq. (33), one can also introduce the magnetic couplings of the axial-vector diquark as those given in Ref. [46] by NRQCD.
Such a term will contribute an extra EM current suppressed by 1/m X in Eq. (53-56).

B. Diquark Transition Currents from Matching
When the recoil is small, to derive the diquark transition currents from symmetries we can assume the diquark as a point-like particle without any internal structure. Therefore, the currents we get in Eq. (41)(42)(43)(45)(46)(47)(48) are only proportional to the constant couplings λ 1 , λ 2 . On the other hand, if the recoil is large, we should consider finite sized diquarks where the transition is dominated by hard internal gluon exchange which can be factorized into short distance coefficients. One way to obtain these short distance coefficients is to perform a matching between DiET and QCD in the large recoil region, where at the quark level one may factorize out a hard kernel, with its tree level form shown in Fig. 1. A hard gluon is exchanged between the two heavy quarks so that the recoil is large, q 2 close to zero, and V µ = γ µ or γ µ γ 5 is the current vertex.
The calculation of the two diagrams in Fig. 1 is straightforward. However, although at tree level we can set the initial and final quarks to be free, the two quark spins are coupled so that the total spin should match with the corresponding diquark spin. Particularly, to match with a scalar or axial-vector diquark, the spinor indexes of the two quarks should be symmetrical or anti-symmetrical. Consider first the X → X transition. By equating the velocities of the initial and final two quarks to be v 1 and v 2 respectively, the amplitude of the two diagrams in Fig. 1 reads where a, b, c, d are spinor indices, and i, j, m, n are color indices. Further, ξ 1 = m Q /(m Q + m b ) and . For the finite-sized diquark, the corresponding weak transition amplitude is Here, X † (v 2 ), X(v 1 ) should be treated as the polarization vectors of the final and initial diquarks, represents the hard kernel. Explicitly, the diquark wave function can be composed of two heavy quark spinors as where i, j, k and β, γ are color and spinor indices, respectively, and N S , N X are normalization factors. Inseting Eq. (59) into Eq. (58) and factorizing an independent color factor C δ k l , one arrives at The tree level matching demands the equivalence of the amplitudes at the quark and the diquark level M QCD = M DiET , thus we can determine the hard kernel as and the color factor is C = −1/3. Similarly, for X → S and S → X transitions, we have Particularly, for the V − A currents, where V µ = γ µ or γ µ γ 5 , the hard kernels are For the EM currents, the X → X, X → S and S → X currents have the same hard kernel as those implies that they are only valid in the large recoil region w → w max .

A. Interpolating Fields
In this section we will focus on semi-leptonic decays of doubly heavy baryons, B bQ → B cQ ℓν.
The transition matrix element of the doubly heavy baryon can be calculated by the reduction where J µ is the current inducing the weak decay. L(P b , P c ) is the operator to pick out the initial and final mass pole residues Φ cQ (x) and Φ bQ (x) are the interpolating fields of the final and initial baryon. Eq. (69) can be expressed both at the quark level and the diquark level. At the quark level, J µ =cγ µ (1 − γ 5 )b, and where χ are the Bargmann-Wigner wave functions [47], where the total spin contributed by the two heavy quarks is j. For a spin-1/2 doubly heavy baryon with j = 0 or j = 1, and a spin-3/2 baryon with j = 1, they are The symmetry indices β, γ project out the spin-1 configuration of the two heavy quarks. The conjugate forms are defined asχ On the other hand, we can equivalently express Eq. (69) at diquark level, with the assumption that the spin-0 and spin-1 heavy diquark field is composed of two heavy quark fields Thus the intepolating field of a doubly heavy baryon can be expressed by the combination of a diquark field and a light quark field Φ 1/2(0) In fact, these normalization factors are related by the heavy flavor symmetry, which leads to However, the relation between N X and N S as well as the relation among N 1/2 , N 1/2(0) and N 3/2 are not obvious. According to Eq. (74) and Eq. (75), we can write the spinor structure of the scalar and axial-vector diquarks in momentum space as Here, we have omitted the color indices. s, s ′ denote the helicity of the spinors u 1 , u 2 , in order.
Since X ss ′ µ has three independent degrees of freedom, while S ss ′ has only one degree of freedom, we can derive the following relation where the sum of all the helicity indices is equivalent to counting the total degrees of freedom. The relations among N 1/2 , N 1/2(0) and N 3/2 can be determined by a similar approach. We transform Eq. (76-78) into the spinor structure in momentum space where r, l, s, s ′ denote the helicities. Since a spin-1/2 particle has two degrees of freedom while a spin-3/2 particle has four, we require the following relations Finally, according to Eq. (81) and Eq. (85) we arrive at where the following properties have been used It should be mentioned that the spinors used here are rescaled from the standard ones as √ m Q u = u QCD . However, as long as we also choose rescaled states as √ m Q | · · · = | · · · QCD , this will never affect our calculations.

B. Transition Matrix Element
With DiET, the transition matrix element defined in Eq. (69) can be calculated at the diquark level. Further, in the heavy diquark limit, utilizing the technique given in Ref. [48], we can reduce the transition matrix element so that it will depend on less unknown form factors. Consider first . The flavor changing current is where j, k are color indices. [Γ ρσ µ ] j k can be factorized as Γ ρσ µ × Cδ j k , and C = −1/3 is given in the last section from matching. To leading power of 1/m 2 X , one can approximate the X µ field as X vµ , so that the ← − ∂ , ∂ in Eq. (88) can be replaced with iv 2 , −iv 1 . According to the reduction formula Eq. (69), the transition matrix element in DiET is where a, b are Dirac indices, while i, j, k, l are color indices. v 1 and v 2 are is four-velocity of the initial and the final baryon, respectively. Using the decoupling transformation defined in Eq. (35), and noting that theX vµ fields are totally decoupled from the soft gluons and also the light quarks, one can factorize the time-ordered matrix element in Eq. (89) to be The last two matrix elements in Eq. (90) can be calculated directly from the free diquark propagator Eq. (34). Using the fact thatχ The dynamics of the light degrees of freedom is completely encapsulated in the following Fourier transformed soft function Next, we need to extract the residues of the mass poles by applying the operator L(P b , P c ) on the correlation function. Near the mass shell, the external momenta P Q can be parameterized as Although the decoupling transformation Eq. In other words, the light particles in the baryon always "know" that they are bound with a heavy diquark. To reflect the confinement, M (k, q; v 2 , v 1 ) is assumed to peak at Operating with L(P Q ) on the denominators, taking the limit ǫ Q , ǫ ⊥ → 0, and noting that there are no poles of 1/ǫ 2 ⊥ , one gets On the other hand, the soft function can be generally parametrized as However, the B(w), C(w), D(w) form factors can be totally absorbed into the the form factor cQ the parametrization takes the form However, if we treat such process by HDiET considering also the heavy flavor symmetry, the number of independent form factors can be greatly reduced. Especially, by combining Eqs. (41)(42)(43) and Eqs. (96-99), one arrives at where only one form factor η(w) is left. This is shared by all the six matrix elements and η(w) is proportional to the soft function A ′ (w) The vector transition shown in Eq. (102) is exactly the same as that given in [39], where the transition matrix element was derived based on heavy quark-diquark symmetry. However, Ref. [39] did not give the result for the axial-current transition. In terms of the complicated factors in Eq. (110), this is determined through the normalization at the zero-recoil point w = 1. From Eq. (41), one can find that the vector current J X→X µ(V ) is conserved ∂ µ J X→X µ(V ) = 0. This implies the conservation of diquark number. Thus we can conclude that where 1 means the diquark number is one. On the other hand, using Eq. (102), and choosing the rest-frame of B bQ (v), v = (1, 0), the same matrix element becomes where we have used γ 0 = / v and / vu = u. Comparing the above two equations, one can conclude that η(1) = 1/6. At the end of last subsection, we have argued that A ′ (w) ∝ ξ(w). Since ξ(1) = 1, it thus follows that η(1) = (1/6)ξ(1).
However, it is necessary to point out that the reduced matrix elements Eqs. (102-109) are only applicable in the region w ∼ 1 or equivalently In the smaller-q 2 region, the large recoil may invalidate the static dynamics of HDiET. As a result, one cannot argue that for any w we have η(w) = (1/6)ξ(w), and an appropriate extension of the form factors from q 2 = q 2 max to q 2 = 0 is necessary. Since the transition matrix elements Eqs. (102-109) are expected to have a lowest-q 2 pole at the mass of B c meson, it is appropriate to multiply η(w) with single pole function B(w) with a suitable normalization B(1) = 1, Finally, we arrive at an explicit expression of the η function Note that for the practical calculation we have to distinguish bewteen the different masses M b , M c .
The Isgur-Wise function was calculated e.g. in Ref. [49], which has the expression GeV 2 is the effective threshold, while M 2 = 3 − 6 GeV 2 is the Borel parameter. In this work, we simply use its center value M 2 = 4.5 ± are the B meson light-cone distribution amplitudes, which have the form where ω 0 = (2/3)Λ [50].    with q µ = P µ b − P µ c the transferred momentum. The f i and g i are related to the F i and G i as form factors, where f ′ 2 (q 2 ) = g ′ 1 (q 2 ) = g ′ 2 (q 2 ) = 0.

D. Semi-Leptonic Decay Widths
Next, using the form factors given in the last section, we will calculate the semi-leptonic decay widths of B , the formula of the differential decay width is given in [8,33] dΓ L dq 2 = G 2 F |V CKM | 2 q 2 p (1 −m 2 l ) 2 384π 3 M 2 The masses of e, µ are neglected here and m τ = 1.78 GeV. Tab. I gives the resulting decay widths and also a comparison with those derived in Ref. [8] within light-front quark model (LFQM). It appears that the two sets of decay width results are consistent.

V. CONCLUSIONS
In summary, we have constructed a heavy diquark effective theory (HDiET), which satisfies the gloabal heavy quark flavor SU(2) symmetry and electromagnetic U(1) symmetry. Imposing these symmetries, we constructed the coupling terms where the diquark fields interact with the external weak and electromagnetic sources. Such coupling terms enable us to obtain the effective diquark transition currents in the small recoil region. On the other hand, for large recoil, the diquark transition currents are derived from the matching between QCD and DiET at tree level.
Furthermore, we simpilfied DiET as HDiET in the heavy diquark limit, from which we reduced the form factors of the doubly heavy baryon transition to only one function η(w). The reduced vector matrix element is the same as those derived by heavy quark-diquark symmetry in earlier works. In addition, we pointed out that η(w) is related with the universal soft function which is proportional to the Isgur-Wise function of heavy meson decays. Thus we obtained the q 2 -dependence of η(q 2 ) by assuming a monopole structure. Finally, the obtained form factors are used to predict the semileptonic decay widths of doubly heavy baryons, and the results are consistent with those derived by LFQM in the earlier works.