Chromomagnetic and chromoelectric dipole moments of quarks in the reduced 331 model

The one-loop contributions to the chromomagnetic dipole moment and electric dipole moment of the top quark are calculated within the reduced 331 model (RM331) for non-zero . It is argued that the results are gauge independent and thus represent valid observable quantities. In the RM331, receives new contributions from two heavy gauge bosons, namely and , and one neutral scalar boson , along with a new contribution from the standard model's Higgs boson via flavor changing neutral currents. The latter, which is also mediated by the gauge boson and the scalar boson , can provide a non-vanishing if there is a -violating phase. The analytical results are presented in terms of both Feynman parameter integrals and Passarino-Veltman scalar functions, which are useful to cross-check the numerical results. Both and are numerically evaluated for parameter values still allowed by the constraints from experimental data. It is found that the new one-loop contributions of the RM331 to the real (imaginary) part of are of the order of ( ), which means at least three orders of magnitude smaller than the standard model prediction but larger than the predictions of other models of new physics. In the RM331, the dominant contribution arises from the gauge boson for in the 30-1000 GeV interval and a mass of the order of a few hundreds of GeV. As for , it receives its largest contribution from exchange and can reach values of the order of , i.e., smaller than the contributions predicted by other standard model extensions.


Puebla, Méxicoμ
Abstract: The one-loop contributions to the chromomagnetic dipole moment and electric dipole moment of the top quark are calculated within the reduced 331 model (RM331) for non-zero . It is argued that the results are gauge independent and thus represent valid observable quantities. In the RM331, receives new contributions from two heavy gauge bosons, namely and , and one neutral scalar boson , along with a new contribution from the standard model's Higgs boson via flavor changing neutral currents. The latter, which is also mediated by the gauge boson and the scalar boson , can provide a non-vanishing if there is a -violating phase. The analytical results are presented in terms of both Feynman parameter integrals and Passarino-Veltman scalar functions, which are useful to cross-check the numerical results. Both and are numerically evaluated for parameter values still allowed by the constraints from experimental data. It is found that the new one-loop contributions of the RM331 to the real (imaginary) part of are of the order of ( ), which means at least three orders of magnitude smaller than the standard model prediction but larger than the predictions of other models of new physics. In the RM331, the dominant contribution arises from the gauge boson for in the 30-1000 GeV interval and a mass of the order of a few hundreds of GeV. As for , it receives its largest contribution from exchange and can reach values of the order of , i.e., smaller than the contributions predicted by other standard model extensions. The anomalous magnetic dipole moment (MDM) and the electric dipole moment are among the lepton properties that have attracted more interest in the experimental and theoretical areas. Currently, there is a discrepancy between the theoretical standard model (SM) prediction of the muon anomalous MDM and its experimental measurement, which might be a hint of new physics [1]. Moreover, any experimental evidence of an electric dipole moment would give a clear signal of new sources of violation as the SM contributions are negligibly small. With the advent of the LHC, anomalous contributions to the coupling have also become a focus of interest. In analogy with the lepton electromagnetic vertex , the anomalous coupling can be written as µ qd q G µυ a T a S U (3) where is the quark chromomagnetic dipole moment (CMDM) and is the quark chromoelectric dipole moment (CEDM), whereas is the gluon field tensor and are the color generators. It is also customary to define the CMDM and CEDM in their dimensionless forms [2] ttg On the experimental side, the search for evidences of the anomalous top quark coupling is underway at the μ tdt LHC [2][3][4]. The most recent bounds on the top quark CM-DM and CEDM were obtained by the CMS collaboration [4,5], which managed to improve the previous bounds [2] by one order of magnitude. Thus, one would expect that tighter constraints on and could be set in the near future.μ t Z ′ As far as the theoretical predictions are concerned, in the SM the CMDM is induced at the one-loop level or higher orders via electroweak (EW) and QCD contributions, whereas the CEDM can only arise up to the threeloop level [6][7][8]. The SM contributions to the on-shell have already been studied in [9][10][11], and more recently the scenario with an off-shell gluon was studied in [12,13] to address some ambiguities of previous calculations, particularly about the on-shell CMDM, which is divergent and meaningless in perturbative QCD. Given that both the top quark CMDM and CEDM could receive a considerable enhancement from new physics contributions, several calculations have been reported in the literature within the framework of extension theories such as the two-Higgs doublet models (THDMs) [14], the fourgeneration THDM [15], models with a heavy gauge boson [11], little Higgs models [16,17], the minimal supersymmetric standard model (MSSM) [18], unparticle models [19], and vector like multiplet models [20]. In this study, we are interested in the contributions to the top quark CMDM and CEDM in the reduced 331 model [21].
The study of elementary particle models based on the gauge symmetry dates back to the 1970s, when it was not clear that Weinberg's model was the right theory of electroweak interactions [22]. After the discovery of the Z and W gauge bosons, given that the electroweak gauge group is embedded into , the so called 331 models [23,24] became serious candidates to extend the SM and explain some issues with no answer, such as the flavor problem and the large splitting between the mass of the top quark and those of the remaining fermions. Several realizations of the 331 model have been proposed in the literature, which predict new fermions, gauge bosons, and scalar bosons. Their phenomenologies have been considerably studied [25][26][27][28][29][30][31][32].
The minimal 331 model [23,24] requires a very large scalar sector that introduces three scalar triplets to give masses to the new heavy gauge bosons and one scalar sextet to endow the leptons with small masses. The complexity of this model has led to the appearance of alternative 331 models aimed to economize the scalar sector. In particular, the reduced 331 model (RM331) [21] only requires two scalar triplets, thereby being considerably simpler than the minimal version [33,34]. In the RM331, the physical scalar states obtained after the symmetry breaking are two neutral scalar bosons only, with the lightest one being identified with the SM Higgs boson [35], and a doubly charged one. Unlike other 331 models, no singly N charged scalar boson arises in the RM331 [36][37][38]. In the  gauge sector, there is one new neutral gauge boson , a  new pair of singly charged gauge bosons  , and a pair  of doubly charged gauge bosons . Similar to other 331 models, the RM331 also predicts three new exotic quarks. The original RM331 is strongly disfavored by experimental data [39], though it would still be allowed as long as left-handed quarks are introduced via a particular representation [40,41], which in fact would give rise to flavor changing neutral current (FCNC) effects.
The contributions to the electron and muon anomalous MDM have been already studied in the RM331 [26] within another 331 realization [42]. As for the CMDM of quarks, there is only a previous calculation in the context of an old version of the 331 model [9], though such a calculation is limited to the on-shell case. However, given that the on-shell CMDM is infrared divergent in the SM [12], a calculation of the off-shell CMDM is mandatory. To the best of our knowledge, there is no calculation of the off-shell CMDM of quarks, let alone their off-shell CEDM, in 331 models. Furthermore, in the model studied in [9], the new contributions only arise in the gauge sector, whereas in the RM331 there are additional contributions from the neutral scalar bosons, which are absent in other 331 models.t tg In this paper, we present a study on the contributions of the RM331 to the off-shell CMDM and CEDM of the top quark. The manuscript is organized as follows. In Section II, we present a brief description of the RM331, with the Feynman rules necessary for our calculations presented in Appendix A. The analytical calculations of the new contributions to the dipole form factors of the vertex are presented in Section III; our results in terms of Feynman parameter integrals and Passarino-Veltman scalar functions are presented in Appendix B. Section IV is devoted to a review of the current constraints on the parameter space of the model and the numerical analysis of the off-shell CMDM and CEDM of the top quark. Finally, in Section V, the conclusions and outlook are presented.

II. BRIEF OUTLINE OF THE RM331
We briefly describe the main features of each sector of the RM331, with a focus only on those details relevant to our calculations.

A. Scalar and gauge boson eigenstates
As far as the scalar sector is concerned, the scalar potential is given by A. I. Hernández-Juárez, G. Tavares-Velasco, A. Moyotl Chin. Phys. C 45, 113101 (2021) where the scalar triplets transform as and . To induce the spontaneous symmetry breaking (SSB), the neutral scalar bosons and develop non-zero vacuum expectation values (VEVs) under the shifting of the fields as which leads to the following constraints The breaks down into the SM gauge group following the pattern where can be identified with the SM Higgs VEV . The left-over of SSB are two neutral scalar bosons and a pair of doubly charged ones , as explained below.
The mass matrix of the neutral scalar bosons in the basis is where . After diagonalization, the mass eigenstates in the limit are with masses where , and . The SM Higgs boson h can be recovered in the limit, thus must be identified with the Higgs boson discovered at the LHC. Given that GeV, from Eq. (8), we obtain the relation [1]. In the case , , and , we obtain , which recovers the SM case and thus .
In the gauge sector, there are two new singly charged gauge bosons , two doubly charged gauge bosons , and a neutral gauge boson . They acquire their masses as follows. The would-be Goldstone bosons are eaten by the singly charged gauge bosons, whereas a linear combination of the doubly charged would-be Goldstone bosons and are absorbed by the doubly charged gauge boson . Also, the orthogonal combination of and gives rise to a physical doubly charged scalar boson pair . Finally, the would-be Goldstone boson becomes the longitudinal components of the gauge boson. Thus, the masses of the new gauge bosons at leading order at are [43] As far as the SM gauge bosons are concerned, the would-be Goldstone bosons and endow the Z and gauge bosons with masses, respectively.

B. Gauge and scalar boson couplings to the top quark
The number of new fermions necessary to fill out the multiplets as well as their quantum numbers depend on the particular 331 model version. There are no new leptons in the RM331, but a new quark is required for each quark triplet. They transform as (1) with the numbers between parentheses representing the field transformations under the gauge group, whereas , , and are the new exotic quarks with electric charges and . Under this representation, the theory is anomaly free [40].

Charged currents
In the quark sector, the charged currents relevant for our calculations are given by the following Lagrangian Chromomagnetic and chromoelectric dipole moments of quarks in the reduced 331 model Chin. Phys. C 45, 113101 (2021) 113101-3 where the family index a runs over 1, 2, and 3, whereas i and l run over 1 and 2. Also stands for the Cabibbo-Kobayashi-Maskawa matrix, with the mixing matrices ( ) transforming the left-handed up (down) quarks flavor eigenstates into their mass eigenstates. It is assumed that the new quarks are given in their diagonal basis. Note that the doubly charged gauge boson does not couple to the top quark.

FCNC currents
Given that the gauge boson couplings to the quarks are non-universal, flavor changing neutral currents (FCNCs) are induced at the tree level. The corresponding Lagrangian for the up quark sector reads where the up quarks are in the flavor basis. It is evident that the above Lagrangian induces FCNC at the tree level after the rotation to the mass eigenstate basis. In addition, the interactions between up quarks and the neutral scalar bosons arise from the Lagrangian where is an up quark triplet and being the quark mass matrix in the flavor basis [40]. After rotating to the mass eigenstate basis, only the terms proportional to are diagonalized, whereas the remaining term gives rise to FCNC couplings, which can be written as where Through the parametrization reported in [44] for the mixing matrices, it is possible to obtain numerical values for the entries of the matrix. Under this framework, , , and .

THE RM331
Apart from the pure SM contributions, at the one-loop level, there are new contributions to the CMDM of the top quark arising in both the gauge and scalar sectors of the RM331. The corresponding Feynman diagrams are depicted in Fig. 1. In the gauge sector, the new contributions arise from the neutral gauge boson and are induced by both diagonal and non-diagonal couplings. There is also a new contribution from the singly-charged gauge boson , which is accompanied by the new exotic quark . As already mentioned, the doubly-charged gauge boson does not couple to the top quark, thus there is no contribution from this gauge boson to the top quark CMDM and CEDM. As for the scalar sector, there are new contributions from the neutral scalar bosons and , which in fact are the novel contributions from the RM331 as they are absent in other 331 model versions. The SM-like Higgs boson yields new contributions arising from its FCNC couplings, which are induced at the tree-level and also from its diagonal coupling, which has a small deviation from its SM value. As for the new Higgs boson , it also contributes via both diagonal and non-diagonal couplings. We would like to point out that such scalar contributions are absent in the 331 model studied in Ref. [9], where the on-shell CMDM of the top quark was calculated. Furthermore, as long as complex FCNC couplings are considered, there are non-vanishing contributions to the CEDM. This class of contributions has not been studied either in the context of 331 models.
We are interested in the off-shell CMDM and CEDM of the top quark. Given that off-shell Green functions are not associated with an S-matrix element, they can suffer from issues such as being gauge non-invariant, gauge dependent, or ultraviolet divergent. In this context, the A. I. Hernández-Juárez, G. Tavares-Velasco, A. Moyotl Chin. Phys. C 45, 113101 (2021) 113101-4 pinch technique (PT) was meant to provide a systematic approach to construct well-behaved Green functions [45], from which valid observable quantities can be extracted. It was later found that there is an equivalence at least at the one-loop level between the results found via the PT and those obtained through the background field method (BFM) via the Feynman gauge [46]. This provides a straightforward computational method to obtain gauge independent Green functions. It is thus necessary to verify whether the RM331 contributions to the CMDM and CEDM of quarks are gauge independent for . Nevertheless, from the Feynman diagrams in Fig. 1, note that the gauge parameter only enters the amplitudes of the Feynman diagrams (a) and (b) via the propagators of the gauge bosons and their associated would-be Goldstone bosons. Those types of diagrams have an amplitude that shares the same structure to those mediated by the electroweak gauge bosons Z and W in the SM, which are known to yield a gauge independent contribution to the CMDM for an off-shell gluon when the contribution of their associated would-be Goldstone bosons are added up. See for instance Ref. [12], where we calculated the electroweak contribution to the CMDM of quarks in the conventional linear gauge and verified that the gauge parameter drops out. Furthermore, the dipole form factors cannot receive contributions from self-energy diagrams, which are required to cancel gauge dependent terms appearing in the monopolar terms via the PT approach. Thus, both the CMDM and CEDM must be gauge independent for an off-shell gluon and hence valid observable quantities.
Below, we present the analytical results of our calculations in a model-independent way, from which the results for the RM331 and other SM extensions would follow easily. The corresponding coupling constants for the RM331 are presented in Appendix A. For loop integration, we used the Passarino-Veltman reduction method, and for completeness, our calculations were also performed through Feynman parametrization via the unitary gauge, which provides alternative expressions to crosscheck the numerical results. The Dirac algebra and the Passarino-Veltman reduction were done in Mathematica with the help of Feyncalc [47] and Package-X [48].

A. New gauge boson contributions
We first consider the generic contribution of a new gauge boson V with the following interaction to the quarks where the coupling constants are taken in general as complex quantities. By hermicity they should obey . The above interaction gives rise to a new contribution to the quark CMDM and CEDM via a Feynman diagram similar to that of Fig. 1(a). The corresponding contribution to the quark CMDM can be written as where we introduced the auxiliary variable , and the function is presented in Appendix B in terms of Feynman parameter integrals and Passarino-Veltman scalar functions. The second term of the righthand side stands for the first term with the indicated replacements. As for the contribution to the quark CEDM, it can arise as long as there are flavor changing complex couplings and is given by where again the function is presented in Appendix B.
From Eqs. (17) and (18), we can straightforwardly obtain the contributions to the quark CMDM and CEDM of the neutral gauge boson and the singly charged gauge boson after replacing the coupling constants and the gauge boson masses. B. New scalar boson contributions Following the same approach described above, we next present the generic contribution to the quark CM-DM and CEDM arising from FCNC mediated by a new scalar boson S, which arises from the Feynman diagram in Fig 1(c). We consider an interaction of the form The above scalar interaction leads to the following contribution to the quark CMDM whereas the corresponding contribution to the quark CEDM is given by where the and functions are presented in Appendix B.
From the above expression we can obtain the contribution of the new scalar Higgs boson of the RM331 as well as the contribution of the SM Higgs boson, which has tree-level FCNC couplings in the RM331.

IV. NUMERICAL ANALYSIS AND DISCUSSION
We next address the numerical analysis. The coupling constants that enter the Feynman rules and are necessary to evaluate the CMDM and CEDM of the top quark [c.f. Eqs. (16) through (21)] are presented in Tables A1 and A2 of Appendix A. Note that these couplings depend on several free parameters, such as the mass parameter , the VEV , the parameters of the scalar potential and , as well as the entries of the matrices , , and . To obtain an estimate of the contributions of the RM331 to the CMDM and CEDM of the top quark we need to discuss the most up-to-date constraints on these parameters from current experimental data.
A. Constraints on the parameter space

Heavy particle masses
As already mentioned, the mass parameter can be identified with the top quark mass [40], whereas the VEV determines the masses of the heavy gauge bosons and the heavy quark . As for the mass of the new scalar bo- , it is determined by the parameters and , along with the VEV , which also determines the mixing angle .
We will first discuss the current indirect constraints on the heavy neutral gauge boson masses. From the muon discrepancy, the constraint TeV [41] was obtained, from which bounds on the heavy gauge boson masses follow. Nevertheless, there are also indirect constraints obtained through the experimental data on oscillations. The RM331 contribution to arises from FCNC couplings mediated by the gauge boson and the and scalar bosons [40,43]. Then, using the parametrization reported in [44], the experimental limit on leads to the following bounds: TeV, TeV and TeV [40]. Similar limits have been imposed using the mass difference in the and systems [43]. In addition, the current experimental bounds on the masses of new neutral and charged heavy gauge bosons from collider searches are model dependent [1]. At the LHC, the ATLAS and CMS Collaborations have searched for an extra charged gauge boson at TeV via the decay modes [49,50] and . The most stringent bounds are obtained for a gauge boson with SM couplings (sequential SM). The respective lower bounds on are 6.0 TeV (5.1 TeV) for the ( ) decay channel, whereas for the decay , the corresponding bound is less stringent, of the order of TeV [51,52]. As far as an extra neutral gauge boson is concerned, the search at the LHC at TeV via its decays into a lepton pair has been useful to impose the lower limit TeV for a gauge boson model arising in the sequential SM and in an -motivated Grand Unification model [53,54]. In this context, it has been pointed out recently that the LHC might be able to constrain the mass of the heavy boson up to the TeV level in several 331 models [55][56][57]. Although these bounds are model dependent and rely on several assumptions, if we consider the conservative value of 5 TeV for the gauge boson masses, we obtain a lower constraint on of the order of 10 TeV. Thus, we will use this value in our analysis to be consistent with experimental constraints and limits from FCNC couplings. As far as direct constraints on the mass of exotic quarks are concerned, the ATLAS and CMS Collaborations have used the TeV data to search for vector-like quarks with electric charge of via their decay into a top quark and a W gauge boson, with the final state consisting of a single charged lepton (muon or electron), missing transverse momentum, and several jets. A mass exclusion limit up to TeV is obtained depending on the properties of the vector-like quark [58][59][60]. We will thus use TeV to be consistent with the experimental bound. According to Eq. (8), the mass of the SM-like Higgs boson receives new corrections through the and parameters. As discussed above, the SM case is recovered when and . Thus, the new corrections to must lie within the experimental error of the SM Higgs boson mass GeV [1]. This allows constraining the and parameters, which in turn translates into constraints on and once the value is fixed. Again, we follow a conservative approach and only consider the experimental uncertainty in the Higgs boson mass, whereas theoretical uncertainties from higher order corrections are not taken into account. Note in Fig. 2 that the allowed regions in the planes vs. and vs. are consistent with the experimental error of the Higgs boson mass at 95% C.L. Note also that for a given , must be approximately one order of magnitude below. In our calculations, we used and , though there is no significant sensitivity of the top quark CMDM and CEDM to mild changes in the values of these parameters. In addition, we found that values ranging from to are allowed for provided that TeV and GeV, which is consistent with recent searches for new neutral scalar bosons at the LHC [1].

Mixing matrices
As for the mixing matrices, we can obtain the absolute values for the entries of the matrices , , and . The entries of the last matrix are given in terms of , , and the matrix elements, and their values are obtained following the parametrization used in [44]. In general, and are expressed in terms of the entries of and , i.e., the complex matrices that diagonalize the mass matrices of up quarks. These matrices can be assumed to be triangular. Then, using the experimental data on quark masses and the mixing angles, it is possible to CP 10 −3 obtain values of their entries [61]. It is also assumed that the only non-negligible mixing is that arising between the third and second fermion families. Furthermore, given that the violation phases are expected to be very small, we follow a conservative approach and assume complex phases of the order of . We present a summary of the values we used in our numerical evaluation in Table 1 .

B. Top quark CMDM
As already mentioned, in the RM331 there are new contributions to the off-shell top quark CMDM arising from the heavy gauge bosons and as well as the neutral scalar bosons and . Below we will use the notation for the contribution of particle A due to Table 1. Values of the parameters used in our evaluation of the CMDM and CEDM of the top quark in the RM331. For the entries of the matrices , , and , we used the values obtained in [40] using the parametrization of [44], where the mass parameter is identified with the top quark mass. We also used and values allowed by the experimental error in the Higgs boson mass and assumed that the only non-negligible mixing is that arising between the third and second fermion families. GeV [1] at 95% C.L. We consider and , which yield the SM limit.
Chromomagnetic and chromoelectric dipole moments of quarks in the reduced 331 model Chin. Phys. C 45, 113101 (2021) coupling. Thus, for instance, will denote the contribution of the loop with the gauge boson due to the coupling. Given that we would like to assess the magnitude of the new physics contributions to , we will extract the pure SM contributions from our calculations. Thus, apart from the contribution due to the treelevel FCNCs of the SM-like Higgs boson , we only consider the contribution arising from the small deviation of the diagonal coupling from the SM coupling. This contribution will be denoted by .
We will examine the behavior of the CMDM of the top quark as a function of , where q is the gluon four-momentum. In the left plot of Fig. 3 we show the real part of the partial contributions to as a function of for the parameter values in Table 1, whereas the real and imaginary parts of the total contribution are shown in the right plot. In general, there is little dependence of on , except for the , , and contributions, which have a change sign. Note also that the contribution is the largest one, whereas the remaining contributions are negligible, with the contribution being the smallest one. Thus, the curve for the real part of the total contribution seems to overlap with that of the contribution, though the former shows a small peak at . This can be explained by the peak appearing in the contribution, which can be as large as the contribution for . We conclude that can have a real part of the order of .μ Concerning the imaginary parts of the partial contributions to , they are several orders of magnitude smaller than the corresponding real parts. As observed in the right plot of Fig. 3, the imaginary part of the total contribution is negligible for but increases up to approximately around GeV, where it starts to decrease up to one order of magnitude as increases up to 1 TeV.
Plots similar to those in Fig. 3 but depicting the behavior of as a function of for GeV and the parameter values in Table 1 are shown in Fig. 4. In this case, note that the real parts of the partial contributions to show a variation of approximately one order of magnitude when increases from 10 TeV to 20 TeV. As already mentioned, the contribution yields the bulk of the total contribution to , whose imaginary part is slightly larger than its real part. Therefore, both real and imaginary contributions of the RM331 to the top quark CMDM can be as large as .
In summary, for TeV, the real part of the new contribution of the RM331 to would be three orders of magnitude smaller than the real part of the SM electroweak contribution [12], whereas its imaginary part can be as large as its real part. In general, there is no appreciable variation in the magnitude of for mild changes in the parameters listed in Table 1 Table 1. The real and imaginary parts of the total contribution are shown in the right plot. Finally, note that the RM331 can render a contribution larger than the ones predicted by other extension models, where a new neutral Z gauge boson is predicted [11]. The real and imaginary parts of the top quark CMDM are of order and respectively in such models.

C. Top quark CEDM
CP η u A potential new source of violation can arise in the RM331 through the FCNC couplings mediated by neutral scalar bosons, which are proportional to the entries of the non-symmetric complex mixing matrix [40], thereby allowing the presence of a non-zero CEDM, which is absent in other 331 models. Thus, it is a novel prediction of the RM331.
There are only two partial contributions to the top quark CEDM in the RM331. Thus, we only analyzed the behavior of the total contribution. Fig. 5 shows the contour lines of the real part (left plot) and the imaginary part (right plot) of in the vs. plane for the parameter values listed in Table 1. We found that the new scalar boson yields the dominant contribution to , whose real (imaginary) part can be as large as ( ), whereas the contribution from the scalar boson is three or more orders of magnitude below. Note also that the real part of decreases as and increase, while the imaginary part remains almost constant. For GeV, the RM331 contribution to the CEDM of the top quark is expected to be below the level, which seems to be much smaller than the values predicted in other extension models [11], where the real and imaginary parts are of order and , respectively. In the range 2 TeV 10 TeV, our results for are enhanced by one order of magnitude, but as already mentioned, this interval is disfavored by current constraints.
For comparison, a compilation of the predictions of several extension models of the top quark CMDM and CEDM for is presented in Table 2. We would like to stress that, to the best of our knowledge, there is no previous estimate of the top quark CEDM in 331 models. Note also that, although these values seem to be much larger than the results obtained for in the RM331, the dipole form factors are expected to decrease as increases. Such a behavior is indeed observed in the SM case [12], where the magnitude of decreases as increases.

V. CONCLUSIONŜ
We present a calculation of the one-loop contributions to the CMDM and CEDM, i.e., and , of the top quark in the framework of the RM331, which is an economic version of the so-called 331 models with a scalar sector composed of two scalar triplets only. We  . We used the parameter values in Table 1.
Chromomagnetic and chromoelectric dipole moments of quarks in the reduced 331 model Chin. Phys. C 45, 113101 (2021) 113101-9 μ t (q 2 ) consider the general case of an off-shell gluon, as it has been pointed out before, so that the QCD contribution to is infrared divergent and the CMDM has no physical meaning for . We argue that the results are gauge independent for and represent valid observable quantities given that the structure of the gauge boson contributions is similar to those arising in the SM. To the best of our knowledge, no previous calculations of the off-shell CMDM and CEDM of the top quark have been presented before in the context of 331 models.

CP
Apart from the usual SM contributions, in the RM331, the CMDM of the top quark receives new contributions from two new heavy gauge bosons and as well as one new neutral scalar boson , along with a new contribution from the neutral scalar boson , which must be identified with the 125 GeV scalar boson detected at the LHC. This model also predicts tree-level FCNCs mediated by the gauge boson and the two neutral scalar bosons and , which at the one-loop level can also give rise to a non-vanishing CEDM provided that there is a -violating phase. The analytical results are presented in terms of both Feynman parameter integrals and Passarino-Veltman scalar functions, which are useful to cross-check the numerical results. We present an analysis of the region of the parameter space of the model consistent with experimental data and evaluate the CMDM and CEDM of the top quark for parameter values still allowed. It was found that the new one-loop contributions of the RM331 to the real (imaginary) part of are of the order of ( ), which are larger than the predictions of other SM extensions [11], with the dominant contribution arising from the gauge boson, whereas the remaining contributions are considerably smaller. It was also found that there is little dependence of on in the 30-1000 GeV interval for a mass of the order of a few hundreds of GeV. As far as the CEDM of the top quark is concerned, it is mainly induced by the loop with exchange and can reach values of the order of for realistic values of the -violating phases. Such a contribution is smaller than the ones predicted by other SM extensions [11].

ACKNOWLEDGMENTS
We acknowledge support from Consejo Nacional de Ciencia y Tecnología and Sistema Nacional de Investigadores. Partial support from Vicerrectoría de Investigación y Estudios de Posgrado de la Benemérita Universidad Autónoma de Puebla is also acknowledged.

APPENDIX A: FEYNMAN RULES
We next present Tables A1 and A2 containing the coupling constants that enter the Feynman rules [26,40,43] that follow from Eqs. (16) and (19). These are neces-sary for the evaluation of the CMDM and CEDM of the top quark in the RM331.

LOOP INTEGRALS
In this appendix, we present the loop integrals appearing in Eqs. (17), (18), (20), and (21) in terms of Feynman parameter integrals and Passarino-Veltman scalar functions both for non-zero and zero . We have verified that all the ultraviolet divergences cancel out. Furthermore, unlike the QCD contribution, all the contributions of the RM331 are finite for .

Feynman parameter integrals
The function in Eq. (17) can be written as Table A1. Coupling constants for the interactions between gauge bosons and quarks in the RM331. We follow the notation of Lagrangian (16). Here are entries of the complex mixing matrix , where the subscript q runs over u and c. This matrix is given in terms of the unitary complex matrix that diagonalizes the mass matrix of up quarks and can be written as . Here, . Table A2. Coupling constants for the interactions between scalar bosons and quarks necessary for the evaluation of the one-loop contributions to the CMDM and CEDM in the RM331. We follow the notation of Lagrangian (19). Here are entries of the complex mixing matrix , where the subscript q runs over u and c. This matrix is given in terms of the unitary complex matrices and that diagonalize the mass matrix of up quarks and can be written as and given that the matrix is not symmetric.
, with and . ] . (B2) As far as the function in Eq. (18) is concerned, it is given by which leads to The function in Eq. (20) is which for simplifies to Finally, the loop function in Eq. (21) reads which yields

Passarino-Veltman results
We next present the results for the loop functions in terms of Passarino-Veltman scalar functions, which can be numerically evaluated by either LoopTools [64] or Collier [65], thereby enabling cross-checking the results. We introduce the following notation for two-and threepoint scalar functions in the customary notation used in the literature: for and . We also define and .
For non-zero , the loop functions in Eqs. (17) and (18) are given by As far as the results for are concerned, they read and The loop functions in Eqs. (20) and (21) are given by and For , we obtain and

Two-point scalar functions
To conclude the manuscript, we present the closed form solutions for the two-point Passarino-Veltman scalar functions appearing in the above calculations. The three-point scalar functions are too lengthy to be shown here.