Beyond basis invariants

TL;DR

This paper proposes using basis-covariant objects instead of basis invariants to analyze physical properties of models, simplifying the identification of key features like CP conservation in multi-Higgs-doublet models.

Contribution

It introduces an alternative approach using basis-covariant objects for analyzing basis invariance, improving clarity and efficiency in complex models.

Findings

Basis-covariant approach simplifies analysis of physical properties.

New necessary and sufficient conditions for CP conservation.

Enhanced understanding of basis invariance in multi-Higgs models.

Abstract

Physical observables cannot depend on the basis one chooses to describe fields. Therefore, all physically relevant properties of a model are, in principle, expressible in terms of basis-invariant combinations of the parameters. However, in many cases it becomes prohibitively difficult to establish key physical features exclusively in terms of basis invariants. Here, we advocate an alternative route in such cases: the formulation of basis-invariant statements in terms of basis-covariant objects. We give several examples where the basis-covariant path is superior to the traditional approach in terms of basis invariants. In particular, this includes the formulation of necessary and sufficient basis-invariant conditions for various physically distinct forms of CP conservation in two- and three-Higgs-doublet models.

Equations

$${\varphi }_a(\vec{r},t)\to X_{ab}{\varphi }_b^{\ast }(-\vec{r},t),X\in U(N).$$

$$V=Y_{ab}({\varphi }_a^{†}{\varphi }_b)+Z_{abcd}({\varphi }_a^{†}{\varphi }_b)({\varphi }_c^{†}{\varphi }_d).$$

$$r_0={\varphi }_a^{†}{\varphi }_a,r_i={\varphi }_a^{†}({\sigma }^i{)}_{ab}{\varphi }_b,i=1,2,3,$$

$$r_0=\frac{1}{\sqrt{3}}{\varphi }_a^{†}{\varphi }_a,r_i={\varphi }_a^{†}(t^i{)}_{ab}{\varphi }_b,i=1,\dots ,8,a=1,2,3.$$

$$[t_i,t_j]=i{f}_{ijk}{t}_k,\{t_i,t_j\}=\frac{1}{3}{\delta }_{ij}{1}_3+d_{ijk}{t}_k,$$

$$d_{ijk}{r}_i{r}_j{r}_k+\frac{1}{2\sqrt{3}}{r}_0(r_0^2-3r_i^2)=0.$$

$$V=M_0{r}_0+M_i{r}_i+{\Lambda }_0{r}_0^2+L_i{r}_0{r}_i+{\Lambda }_{ij}{r}_i{r}_j.$$

$$F_i:=f_{ijk}{a}_j{b}_k,D_i:=\sqrt{3}{d}_{ijk}{a}_j{b}_k.$$

$$\left(\!\begin{array}[]{ccc}c_{\alpha}&s_{\alpha}\\ -s_{\alpha}&c_{\alpha}\end{array}\!\right)\quad\mbox{as in Ref.~{}\cite[cite]{\@@bibref{Authors Phrase1YearPhrase2}{Ecker:1987qp}{\@@citephrase{(}}{\@@citephrase{)}}},}\quad\mbox{or}\quad\left(\!\begin{array}[]{ccc}0&e^{i\alpha}\\ e^{-i\alpha}&0\end{array}\!\right)\quad\mbox{as in Ref.~{}\cite[cite]{\@@bibref{Authors Phrase1YearPhrase2}{Weinberg:1995mt}{\@@citephrase{(}}{\@@citephrase{)}}}.}$$

$${\mathcal{I}}_1=(M,M^{(1)},M^{(2)}),{\mathcal{I}}_2=(L,L^{(1)},L^{(2)}),{\mathcal{I}}_3=(M,L,M^{(1)}),{\mathcal{I}}_4=(M,L,L^{(1)}).$$

$$\mathrm{Tr}\Lambda ,\mathrm{Tr}{\Lambda }^2,\mathrm{Tr}{\Lambda }^3,m_0\equiv M_i{M}_i,m_1=M_i{\Lambda }_{ij}{M}_j,m_2=M_i({\Lambda }^2{)}_{ij}{M}_j.$$

$$m_3$$

$$m_4$$

$$\mathcal{I}=I_{3Y3Z}={ϵ}_{ijk}{M}_i(\Lambda M{)}_j({\Lambda }^{2}M{)}_k.$$

$${\mathcal{I}}^2=m_0{m}_2{m}_4+2m_1{m}_2{m}_3-m_0{m}_3^2-m_2^3-m_4{m}_1^2.$$

$$m_k={\lambda }_1^k{M}_1^2+{\lambda }_3^k{M}_3^2,k=0,1,\dots .$$

$$CP(x_2\to -x_2):m_2-m_1({\lambda }_1+{\lambda }_3)+m_0{\lambda }_1{\lambda }_3=0.$$

$$CP(x_1\to -x_1):$$

$$CP(x_3\to -x_3):$$

$${\left[M_1{M}_2{M}_3({\lambda }_1-{\lambda }_2)({\lambda }_2-{\lambda }_3)({\lambda }_3-{\lambda }_1)\right]}^2=0,$$

$$m_1^2=m_0{m}_2.$$

$$I_N=2|f_{ijk}{e}_i{e}_j^{\prime }{e}_k^{\prime \prime }|.$$

$$X=\left(\!\begin{array}[]{ccc}0&1&0\\ -1&0&0\\ 0&0&1\\ \end{array}\!\right)\,.$$

$$x_8\to x_8,(x_1,x_2,x_3)\to -(x_1,x_2,x_3)$$

$$x_4\to x_6,x_6\to -x_4,x_5\to -x_7,x_7\to x_5.$$

$$\Lambda =\left(\begin{array}{ccc}\text{framebox}{}_{3\times 3}& 0& 0\\ 0& \text{framebox}{}_{4\times 4}& 0\\ 0& 0& {\Lambda }_{88}\end{array}\right),$$

$$f_{123}=1,f_{147}=-f_{156}=f_{246}=f_{257}=f_{345}=-f_{367}=\frac{1}{2},f_{458}=f_{678}=\frac{\sqrt{3}}{2},$$

$$d_{146}=d_{157}=-d_{247}=d_{256}=\frac{1}{2},d_{344}=d_{355}=-d_{366}=-d_{377}=\frac{1}{2},$$

$$d_{118}=d_{228}=d_{338}=-d_{888}=\frac{1}{\sqrt{3}},d_{448}=d_{558}=d_{668}=d_{778}=-\frac{1}{2\sqrt{3}}.$$

$$r_1+i{r}_2={\varphi }_1^{†}{\varphi }_2,r_4+i{r}_5={\varphi }_1^{†}{\varphi }_3,r_6+i{r}_7={\varphi }_2^{†}{\varphi }_3,$$

$$r_3=\frac{1}{2}({\varphi }_1^{†}{\varphi }_1-{\varphi }_2^{†}{\varphi }_2),r_8=\frac{1}{2\sqrt{3}}({\varphi }_1^{†}{\varphi }_1+{\varphi }_2^{†}{\varphi }_2-2{\varphi }_3^{†}{\varphi }_3).$$

Full text

I Introduction

When describing the Standard Model (SM) or building models beyond the SM, one always faces the notorious freedom of basis-choices. Complex fields can be rephased, yet the physics emerging from the model must be invariant under these rephasings. In models with several fields with identical quantum numbers, the freedom of basis choices is even larger and includes arbitrary rotations in the space of these fields. One may fix a basis for the initial fields, then arrive at the physical (mass eigenstate) fields and explore their phenomenology. Or one can switch to a different basis and explore the phenomenology there. Although the Lagrangian and the intermediate calculations may look vastly different in different bases, the observables must be the same.

This utterly obvious statement may look less obvious when one actually gets down to practical calculations. Parameters of the Lagrangian depend on the basis choice, and attributing physical importance to them can only be done when the basis choice is also specified. For example, in the two-Higgs-doublet model (2HDM) Lee:1973iz ; Branco:2011iw , the two doublets can acquire non-zero vacuum expectation values $v_{1}$, $v_{2}$, whose ratio is customarily denoted as $\tan\beta=v_{2}/v_{1}$. Although the vast majority of papers on 2HDM phenomenology describe measurable quantities in terms of $\tan\beta$, this parameter is by itself basis-dependent and not an observable. It can become an observable in bases fixed by additional requirements, such as when the $\mathbb{Z}_{2}$ symmetry responsible for natural flavor conservation is manifest Haber:2006ue .

Another vivid example, still within the class of $N$-Higgs-doublet models (NHDM), are the conditions for explicit $CP$-conservation book . A common approach is to simply define a $CP$ transformation in the space of complex scalar fields $\phi_{a}$ $(a=1,\dots,N)$, via $\phi_{a}(\vec{r},t)\to\phi_{a}^{*}(-\vec{r},t)$. Using such a definition, one would naively think that the model explicitly violates $CP$ symmetry, if the scalar potential contains complex coefficients. This is, however, not true in general. This standard definition of how $CP$ acts on scalar fields is basis-dependent. One can define the general $CP$ transformation via Ecker:1981wv ; Ecker:1983hz ; Neufeld:1987wa ; Ecker:1987qp

[TABLE]

If a model is invariant under such a transformation with any matrix $X$, then it is explicitly $CP$ conserving, regardless of whether the potential has complex coefficients Grimus:1989qn . Although the matrix $X$ does depend on the basis choice, the presence of such a symmetry certainly is a basis-independent fact and has observable consequences.

These and other simple examples have led the model building community to appreciate basis-invariant combinations of the parameters of the Lagrangian, or more simply basis invariants. The general procedure for construction of such quantities was presented in Botella:1994cs : recognize transformation properties of the parameters under general basis changes, rewrite them as tensors, and fully contract various tensors to obtain basis invariants of the model. A nice illustration of this strategy is given by the NHDM scalar sector and, in particular, by the issue of $CP$ conservation in the 2HDM Botella:1994cs ; book ; Branco:2005em ; Davidson:2005cw ; Gunion:2005ja ; Varzielas:2016zjc ; Trautner:2018ipq .

There is no doubt that, in any model, all physical observables must be expressible in terms of basis invariants. A major problem is that, beyond the simplest cases, these expressions become exceedingly or even prohibitively complicated. The key message of our paper is: it is not obligatory to formulate physically relevant, basis-invariant statements exclusively in terms of basis invariants. Rather, it is also possible to formulate them in terms of basis-covariant objects. These objects do transform under basis changes but, loosely speaking, their relative properties are basis invariant.

Perhaps surprising, formulating basis-invariant statement in terms of basis-covariant objects sometimes leads to dramatic simplifications as compared to equivalent statements formulated in terms of basis invariants directly. Below, focusing on $CP$ conservation in 2HDMs and 3HDMs, we will collect a few remarkable illustrations. In particular, this includes cases where results in terms of basis invariants are not yet known due to their exceeding complexity. This gives convincing arguments that in sufficiently sophisticated models, working with basis-covariant objects is the method of choice.

The paper is organized as follows. In the next Section we will review the bilinear formalism, which is particularly useful in order to find basis-covariant in the NHDM. Then we describe the problem of explicit $CP$ conservation and various possibilities which exist in the NHDM. Next, we describe how this issue was solved in the 2HDM and 3HDM, and finally summarize our findings.

II Bilinear formalism

II.1 Bilinears

Let us start with a brief review of the bilinear formalism of the 2HDM Nagel:2004sw ; Ivanov:2005hg ; Nishi:2006tg ; Maniatis:2006fs ; Maniatis:2007vn ; Ivanov:2006yq ; Ivanov:2007de ; Nishi:2007dv . The most general renormalizable 2HDM Higgs potential constructed from two Higgs doublets $\phi_{a}$, $a=1,2$ can be compactly written as

[TABLE]

It depends on the Higgs fields via gauge-invariant combinations $\phi_{a}^{\dagger}\phi_{b}$, which can be arranged into components of a real-valued bilinears

[TABLE]

where $\sigma^{i}$ are the familiar Pauli matrices. Each $(r_{0},r_{i})$ in (3) is in one-to-one correspondence with an electroweak gauge orbit in the space of doublets $\phi_{a}$. The map (3) from doublets $\phi_{a}$ to $(r_{0},r_{i})$ does not cover the entire $1+3$-dimensional space but only the region defined by inequalities $r_{0}\geq 0$ and $r_{0}^{2}-r_{i}^{2}\geq 0$. A basis change transformation $\phi_{a}\to\phi^{\prime}_{a}=U_{ab}\phi_{b}$ with $U\in U(2)$ leaves $r_{0}$ invariant and induces an $SO(3)$ rotation of the vector $r_{i}$. Since the map $SU(2)\to SO(3)$ is surjective, any $SO(3)$ rotation in the bilinear space can be realized as a basis change in the space of two doublets.

Beyond two Higgs doublets, the approach remains the same but the complexity of the problem skyrockets. In the 3HDM Nishi:2006tg ; Ivanov:2010ww ; Maniatis:2014oza , we define $1+8$ gauge-invariant bilinear combinations $(r_{0},r_{i})$:

[TABLE]

Here, $t_{i}=\lambda_{i}/2$ are the generators of the $SU(3)$ algebra satisfying

[TABLE]

with the $SU(3)$ structure constants $f_{ijk}$ and the fully symmetric $SU(3)$ invariant tensor $d_{ijk}$. Explicit expressions for the components of $r_{i}$ and a list of non-zero components of the $SU(3)$ invariant tensors are given in the appendix. Group-theoretically, $r_{0}$ is an $SU(3)$ singlet while $r_{i}$ transforms in the adjoint representation of $SU(3)$.

Unlike in the 2HDM, where the bilinears only had to satisfy $r_{0}\geq 0$ and $r_{0}^{2}-r_{i}^{2}\geq 0$, in the 3HDM they must satisfy an additional constraint Ivanov:2010ww :

[TABLE]

Under a basis change in the space of Higgs doublets, $\phi_{a}\to\phi^{\prime}_{a}=U_{ab}\phi_{b}$ with $U\in SU(3)$, $r_{0}$ is invariant while $r_{i}$ rotates as a vector of $SO(8)$. However, not all $SO(8)$ rotations in the adjoint space can be obtained in this way; they must in addition obey the constraint (6) and, therefore, conserve the contraction $d_{ijk}r_{i}r_{j}r_{k}$.

II.2 Constructions in the adjoint space

When passing from Higgs doublets to bilinears, the space we work in becomes more complicated but the objects we study get simpler. The potential $V$ becomes a quadratic, rather than quartic function of variables,

[TABLE]

The bilinear approach and the generic expression for $V$ above hold for any NHDM. All components of the tensors $Y_{ab}$ and $Z_{abcd}$ in (2) fill $M_{0}$, $\Lambda_{0}$, the entries of the real vectors $M$ and $L$ (lying in the adjoint space111Here, we distinguish between the adjoint space of arbitrary real vectors $x\in\mathbb{R}^{N^{2}-1}$ transforming under the adjoint representation of $SU(N)$ and the orbit space which is spanned by those vectors $r\in\mathbb{R}^{N^{2}-1}$ which can be constructed from the doublets and satisfy (6) and the inequalities above. $\mathbb{R}^{N^{2}-1}$), as well as the $(N^{2}-1)\times(N^{2}-1)$ real symmetric matrix $\Lambda$.

In the 2HDM, any $SO(3)$ rotation can be induced by a basis change. Therefore, the matrix $\Lambda$ can always be diagonalized and its eigenvectors can always be aligned with the axes of the adjoint space $(x_{1},x_{2},x_{3})$. These eigenvectors as well as the vectors $M$, $L$ are covariant objects and transform in the same way under basis changes. Using $SO(3)$ invariant tensors $\delta_{ij}$ and $\epsilon_{ijk}$, one can contract these vectors and obtain basis invariants.

For the 3HDM, the potential (7) contains two 8D vectors $M$ and $L$ and the $8\times 8$ real symmetry matrix $\Lambda$. The lack of complete $SO(8)$ rotational freedom implies that it is not guaranteed anymore that $\Lambda$ can be diagonalized by a Higgs-basis change. Nevertheless, $\Lambda$ can always be expanded over its eigensystem, and eigenvalues and eigenvectors can be found numerically.

The fact that $SU(3)$ basis changes do not offer the full $SO(8)$ rotational freedom in the adjoint space certainly feels like a major nuisance factor. However it also offers at our disposal two additional invariant tensors $f_{ijk}$ and $d_{ijk}$. One can use them to define $f$- and $d$-products of any pair of vectors $a$ and $b$ from the adjoint space:

[TABLE]

These products respect group covariance: vectors $F$ and $D$ transform as the adjoint $SU(3)$ representations. These new products are at the heart of the basis-invariant algorithms for detection of various $CP$ symmetries in the 3HDM Nishi:2006tg ; Ivanov:2018ime .

III Explicit $CP$-conservation in the 2HDM

III.1 Different forms of $CP$ symmetry in the 2HDM

Consider the scalar sector of the 2HDM and suppose that it explicitly conserves $CP$. Does this statement unambiguously specify the model (up to basis changes)? The answer is no. There exist several distinct forms of $CP$ symmetry, which cannot be mapped one to another by any basis change. Depending on what kind of $CP$ symmetry one imposes, one obtains physically distinct models. This fact is known since long ago Ecker:1987qp ; book ; Weinberg:1995mt ; Grimus:1995zi ; its application to the 2HDM was discussed at length, for example, in Ferreira:2009wh ; Ferreira:2010yh . Here, we briefly repeat the classification to set up the notation.

The general $CP$ transformation defined in (1) depends on the matrix $X$, whose form is basis-dependent. However any $CP$ transformation possesses a basis-invariant feature: its order, that is, how many times one must apply it to obtain the identity transformation. Starting from an arbitrary unitary $X$, one can bring it to a block-diagonal form Ecker:1987qp ; Weinberg:1995mt , which has on the diagonal either pure phase factors or $2\times 2$ matrices of the following type:

[TABLE]

Applying the $CP$ transformation twice results in a Higgs family transformation with matrix $XX^{*}$. If it happens that $XX^{*}={\boldsymbol{1}}$, which takes place at $\alpha=0$ or $\pi$, the $CP$ transformation is of order 2, which we will generically denote by CP2. If $XX^{*}\not={\boldsymbol{1}}$ but $(XX^{*})^{k}={\boldsymbol{1}}$, which requires $\alpha$ to be a multiple of $\pi/k$, we get a $CP$ transformation of order $2k$ denoted as CP$2k$. If no finite $k$ exists such that $(XX^{*})^{k}={\boldsymbol{1}}$, that is, if $\alpha/\pi$ is not a rational number, we say that the $CP$ transformation is of infinite order, which we will denote as CP$\infty$.

Let us now list the options available in the 2HDM.

CP2. For a CP2 transformation, there always exists a basis in which it takes the standard form $\phi_{a}\to\phi_{a}^{*}$, that is, $X={\boldsymbol{1}}$. This model is referred to as the $CP$-conserving 2HDM; in the classification of Ferreira:2009wh ; Ferreira:2010yh it was denoted as $\mathfrak{CP}1$.222Reference Ferreira:2009wh introduced the notation, also used in Ferreira:2010yh , of $\mathfrak{CP}1$, $\mathfrak{CP}2$, and $\mathfrak{CP}3$, for models obtained in the 2HDM by a GCP symmetry where $XX^{*}=1$, $XX^{*}=-1$, and neither 1 nor $-1$, respectively. We use here the fraktur symbol because our CPn refers to a symmetry where $(CP)^{n}=1$. Thus, in the 2HDM, CP2 $\mathfrak{\mathrel{\widehat{=}}CP}1$, CP4 $\mathrel{\widehat{=}}\mathfrak{CP}2$, and CPn with $n>4$ or CP$\infty$ $\mathrel{\widehat{=}}\mathfrak{CP}3$. In the adjoint space, the standard $CP$ transformation corresponds to the mirror reflection: $x_{1,3}\to x_{1,3}$, $x_{2}\to-x_{2}$. In a different basis, this transformation is still a mirror reflection but with respect to a different axis in the bilinear space. Due to the full $SO(3)$ rotational freedom, any mirror reflection with respect to an arbitrary direction in the bilinear space can be transformed via a basis change to the reflection with respect to $x_{2}$.

The necessary and sufficient condition for the 2HDM potential to possess a CP2 symmetry is the existence of the real basis, that is, a basis in which all coefficients are real Gunion:2005ja . The challenge is how to detect the existence of the real basis in a basis-invariant way. Below, we will list two approaches to solve this problem. 2. 2.

CP4. This transformation implies $XX^{*}\not={\boldsymbol{1}}$, but $(XX^{*})^{2}={\boldsymbol{1}}$, which requires $\alpha=\pi/2$ in (9). This transformation has a remarkably simple geometric interpretation in the adjoint space: it is the point reflection $x_{i}\to-x_{i}$. In the classification of Ferreira:2009wh ; Ferreira:2010yh it was denoted as $\mathfrak{CP}2$. The 2HDM incorporating $\mathfrak{CP}2$ was dubbed in Maniatis:2007de the maximally $CP$-symmetric model.

This geometric picture clearly shows that imposing CP4 on the 2HDM scalar potential is equivalent to simultaneously imposing three CP2s, each performing a mirror reflection with respect to axes $x_{1}$, $x_{2}$, and $x_{3}$. Thus, the maximally $CP$-symmetric model is certainly distinct from the usual $CP$-conserving 2HDM. 3. 3.

Higher-order $CP$. If $\alpha$ in Eq. (9) is not a multiple of $\pi/2$, there is no further simplification possible, and the transformation manifests itself in the bilinear space as a generic rotary reflection (improper rotation). It was dubbed as $\mathfrak{CP}3$ in the classification of Ferreira:2009wh ; Ferreira:2010yh . Group-theoretically, one can define transformations of finite or infinite order, but all of them have the same effect on the scalar potential of the 2HDM: the potential will be invariant under a $O(2)\times\mathbb{Z}_{2}$ symmetry group in the bilinear space. 4. 4.

Combining two CP2s. Finally, one can construct a 2HDM by imposing two different CP2 symmetries at once. Depending on the choice of these symmetry transformations and on their commutation properties, one can end up with different models. In particular, one can arrive in this way at a model which cannot be obtained just by imposing any single GCP. For example, if one CP2 is the standard $CP$ (mirror reflection with respect to $x_{2}$) and the other CP2 is based on $X=\mathrm{diag}(1,\,-1)$ (reflection with respect to $x_{1}$), then the two GCPs commute, and the resulting model is known as the $\mathbb{Z}_{2}$-symmetric 2HDM as it can be enforced by a single sign flip of one of the Higgs doublets.333Group-theoretically, imposing two different GCPs, in this case, is equivalent to imposing one GCP and a group of family symmetries. Thus, in this way one does not arrive at a completely new, previously overlooked 2HDM.

It turns out that for all of the cases listed above, the total symmetry group can always be factorized as $\mathrm{CP2}\times H$, where $H$ is a family symmetry not involving any $CP$ transformation. In other words, the more exotic $CP$ transformations of the cases 2–4 above are related to CP2 by adjoining it with elements from $H$, which act in the adjoint space as

$180^{\circ}$ rotation in the $(x_{1},x_{3})$ plane; 2. 3.

rotation by any angle in the $(x_{1},x_{3})$ plane and $180^{\circ}$ rotation in the $(x_{1},x_{2})$ plane; 3. 4.

$180^{\circ}$ rotation in the $(x_{1},x_{2})$ plane.

Consequently, for the 2HDM there are two ways to detect the enlarged symmetry groups: either by detecting the presence of multiple $CP$ symmetries (corresponding to the composition of usual $CP$ with unitary symmetries), or by directly detecting the unavoidable presence of the additional unitary symmetries. Crucially, this does not apply for the 3HDM where there are genuinely distinct $CP$ symmetries that cannot be factored out, see Sec. IV. Here, for the 2HDM, we focus on the first strategy and show how to detect the presence of multiple $CP$ symmetries.

III.2 Explicit $CP$-conservation via $CP$-odd basis invariants

All the different forms of $CP$ symmetry have a common consequence: all $CP$-odd physical observables are zero. Thus, to detect physical $CP$ invariance with respect to any form of $CP$ transformation, one has to make sure that all $CP$-odd basis invariants are zero.

Constructing a $CP$-odd basis invariant out of the couplings $Y_{ab}$ and $Z_{abcd}$ is a rather straightforward exercise Botella:1994cs . The challenge is to find the minimal number of $CP$-odd invariants such that setting them to zero implies that all other $CP$-odd invariants are zero as well. This problem was first solved in 1994 for the 2HDM after electroweak symmetry breaking Lavoura:1994fv , while the solution for the 2HDM before symmetry breaking was discovered in the mid-2000’s Branco:2005em ; Davidson:2005cw ; Gunion:2005ja . Four $CP$-odd basis-invariant combinations were constructed, labeled $I_{Y3Z}$, $I_{2Y2Z}$, $I_{3Y3Z}$, $I_{6Z}$ in Gunion:2005ja , according to the powers of tensors $Y_{ab}$ and $Z_{abcd}$ used. Setting these invariants to zero ($I_{i}=0$) implies that all other $CP$-odd invariants vanish, too. This gives the necessary and sufficient condition for the 2HDM scalar sector to be explicitly $CP$-conserving, which can be checked in any basis. Very recently, a powerful method based on Hilbert series and plethystic logarithm was proposed in Trautner:2018ipq which allows one to efficiently construct the full ring of $CP$-even and $CP$-odd invariants. This offers a shortcut to find these four $CP$-odd basis invariants, and allows for a concise proof that the vanishing of these four invariants is indeed sufficient for explicit $CP$ conservation.

$CP$-odd basis invariants can also be constructed using the bilinear formalism Ivanov:2005hg ; Nishi:2006tg ; Maniatis:2007vn . Since $CP$-transformations always correspond to reflections in the adjoint space, $CP$-odd invariants can be constructed as triple products of the vectors constructed from $L$, $M$, and $\Lambda$. Defining $L^{(p)}_{i}:=(\Lambda^{p})_{ij}L_{j}$ and $M^{(p)}_{i}:=(\Lambda^{p})_{ij}M_{j}$, and denoting the triple product as $(A,B,C):=\epsilon_{ijk}A_{i}B_{j}C_{k}$, one can construct the four $CP$-odd invariants

[TABLE]

Also, one can prove that the model is explicitly $CP$-conserving if and only if all four invariants ${\cal I}_{i}=0$. The relation between these ${\cal I}_{i}$ and $I_{Y3Z}$, $I_{2Y2Z}$, $I_{3Y3Z}$, and $I_{6Z}$ of Gunion:2005ja was established in Nishi:2006tg .

The conditions $I_{i}=0$ or ${\cal I}_{i}=0$ indicate that the 2HDM possesses a $CP$ symmetry. However, one cannot distinguish whether it is just a single CP2, or a higher-order $CP$, or a combination of several $CP$ symmetries imposed simultaneously. These different forms of $CP$ invariance do lead to physically distinct $CP$-conserving 2HDMs. But $CP$-odd invariants, by construction, cannot tell the difference between the physically distinct models as they are not sensitive to the matrix $X$ in (1). Thus, it is mandatory to go beyond $CP$-odd invariants in order to recover this information.

III.3 Explicit $CP$-conservation via $CP$-even basis invariants

It is known that the presence of $CP$ violation can be detected exclusively via $CP$-even invariants. For example, in the Standard Model, precise quark sector measurements of $|V_{ud}V_{us}|$, $|V_{cd}V_{cs}|$, and $|V_{td}V_{ts}|$ imply that the unitarity triangle has a non-vanishing area, which is a measure of $CP$ violation.444We thank João Silva for reminding us of this example. What is more important is that $CP$-even invariants, being non-zero for $CP$ conserving models, can reveal which form of $CP$ symmetry is imposed. The distinction appears in the form of additional relations among these invariants.

Let us illustrate this statement using the bilinear formalism. For the sake of the argument it suffices to treat the simplified case where $L=0$ by assumption. We are then left with one real 3D vector $M$ and the real symmetric $3\times 3$ matrix $\Lambda$. In total, they have $6+3=9$ components. The basis-change freedom is characterized in the bilinear space by the group $SO(3)$. Thus, all inequivalent models can be characterized by $9-3=6$ basis-invariant parameters. A possible choice is555As stressed in Trautner:2018ipq , the trace basis for invariants of $\Lambda$ may not be the most convenient choice for many applications. Nevertheless, it suffices for our argument here.

[TABLE]

All higher-order invariants are then expressible in terms of these six invariants, for example,

[TABLE]

Suppose, the model had a CP2 symmetry. Since we have set $L=0$, the only remaining non-trivial $CP$-odd invariant in Eq. (10) is

[TABLE]

The condition ${\cal I}=0$ is compact and basis invariant, but it does not distinguish what particular CP2 we have imposed: The usual $CP$, which amounts to $x_{2}\to-x_{2}$, a different mirror reflection $x_{1}\to-x_{1}$, or both of them simultaneously. One can also square ${\cal I}$ in order to express it exclusively via $CP$-even invariants,

[TABLE]

Using (12) and (13), this can further be reduced to the invariants of Eq. (11). Setting this expression to zero represents the basis-invariant condition for explicit $CP$ conservation written in terms of $CP$-even invariants. However, this single relation is still equivalent to ${\cal I}=0$ and cannot by itself distinguish which, or how many reflections are imposed in addition. We need an additional relation among $CP$-even invariants to settle the issue.

To derive it, let us first fix the basis. The intermediate relations will rely explicitly on this basis choice, but the final result will be basis-independent.666It is crucial that basis-independent relations are sought in the end. The specific basis choice is just an auxiliary tool for constructing basis-independent relations. In general, it is not recommended to start the classification of symmetries in a basis where $\Lambda$ is diagonal. One may run into conditions which are not renormalization group invariant. That is, they would not correspond to actual symmetry classes. A specific example is given in Eqs. (129)-(136) of Ferreira:2010yh .

We choose the basis in which $\Lambda$ is diagonal, with eigenvalues $\lambda_{1,2,3}$, which are, generically, non-zero and distinct, while $M$ has initially three non-zero components. In this basis, invariance under the usual $CP$ ($\mathfrak{CP}1$) implies $M_{2}=0$, which allows us to write the invariants $m_{k}$ as

[TABLE]

From here we deduce an extra relation among invariants $m_{k}$:

[TABLE]

Although the set of eigenvalues is basis-independent, their ordering depends on the basis choice, and this relation explicitly distinguishes $\lambda_{2}$ from $\lambda_{1,3}$.

Imposing a different mirror reflection would produce the analogous relations,

[TABLE]

Thus, after fixing the basis, we do distinguish among different mirror-reflection symmetries. As a cross check, to return to the basis-independent formulation of the same condition, one can multiply all three expressions and set the product to zero. After some algebra we get

[TABLE]

which is exactly ${\cal I}^{2}=0$.

Let us now return to the fixed basis and impose (17) and (18) simultaneously. Assuming that the $\lambda_{i}$ are different, we first get two simplified basis-dependent relations, $m_{1}=m_{0}\lambda_{3}$ and $m_{2}=m_{1}\lambda_{3}$, and in general $m_{k+1}=m_{k}\lambda_{3}$, from which we deduce a new basis-independent relation:

[TABLE]

Thus, imposing condition (21) not only guarantees that the model is $CP$-conserving (which follows from direct substitution in Eq. (15), showing that ${\cal I}^{2}$ vanishes), but also fixes the specific $CP$-symmetry to be the one of case 4: namely, simultaneously imposing two commuting mirror reflections. This relation also signals the presence of a unitary $\mathbb{Z}_{2}$ symmetry, distinct of $CP$, with a clear geometrical interpretation, see Sec. III.4. This demonstrates that relations between $CP$-even invariants can indicate the presence of symmetries other than $CP$.

If one imposes the three commuting CP2 symmetries simultaneously, which in the 2HDM is equivalent to imposing CP4, one notices, in the $\Lambda$-diagonal basis, that each component of the vector $M$ must be zero, and therefore $m_{0}=0$. This is the $CP$-even basis-invariant condition for existence of the CP4 symmetry.

Notice that we needed to fix a basis here and perform the intermediate calculations in a basis-dependent manner in order to arrive at a new relation among $CP$-even invariants. This is not necessary in general and a universal approach to derive such relations, even in the absence of the simplifying assumption $L=0$, has recently been introduced in Trautner:2018ipq . Because the final relation occurs between basis-invariant quantities, and since it does correspond to a legitimate symmetry, it is also renormalization group invariant.

III.4 Explicit $CP$-conservation via basis-covariant objects

The bilinear space approach outlined in section II offers a more direct insight into the structural properties of the scalar potential. Due to the full $SO(3)$ rotational freedom in the adjoint space, the scalar sector of the 2HDM is fully specified by the eigenvalues of the matrix $\Lambda$ and by the orientation of the two real vectors $L$ and $M$ with respect to the eigenvectors of $\Lambda$. We stress that although all of these vectors are basis covariant objects, their relative orientation can be specified in basis-invariant terms. It is this relative orientation that encodes various forms of $CP$ symmetry in the 2HDM Ferreira:2010yh .

A model is explicitly $CP$-conserving if and only if there exists an eigenvector of $\Lambda$ which is orthogonal to both $M$ and $L$ Ivanov:2005hg ; Nishi:2006tg ; Maniatis:2007vn ; Ferreira:2010yh . This geometric criterion is basis-invariant. Using linear algebra, one can rewrite this geometric observation in the form of four $CP$-odd invariants being equal to zero, cf. (10). In fact, this is how these invariants were found in Nishi:2006tg . 2. 2.

Since CP4 implies the point reflection in the bilinear space, $x_{i}\to-x_{i}$, imposing it on the 2HDM implies $M=0$ and $L=0$, with no restriction on $\Lambda$. Each individual mirror reflection $x_{1}\to-x_{1}$, $x_{2}\to-x_{2}$, and $x_{3}\to-x_{3}$ along each eigenvector of $\Lambda$ is then a valid $CP$ symmetry. 3. 3.

Imposing symmetry under a higher-order $CP$ transformation implies that, in addition to $M=L=0$, the matrix $\Lambda$ has a pair of degenerate eigenvalues. This implies infinitely many mirror reflections as symmetries. 4. 4.

Finally, imposing two commuting CP2 symmetries means that there exist two eigenvectors of $\Lambda$ orthogonal to both $M$ and $L$. Within the 3D adjoint space, this is equivalent to the statement that there exists an eigenvector of $\Lambda$ parallel to both $M$ and $L$. The algebraic expression (21) encodes precisely this information.

There are two lessons that we learn from this list. Firstly, formulating physically relevant features of the model solely in terms of basis invariants is not the only option. Basis-invariant statements can also be expressed in terms of basis-covariant objects. Secondly, and perhaps more importantly, we see that basis-invariant statement in terms of basis-covariant objects can often be formulated more concisely and derived more directly than via the brute force contraction of tensors. This applies to both, $CP$-even and $CP$-odd invariants.

Having seen the three approaches to establishing various forms of explicit $CP$ conservations in the 2HDM and having established a “dictionary” between some of them, one may be tempted to think that it is just a matter of taste which approach to use. However, in the next section we will consider the problem of explicit $CP$ conservation in the 3HDM and demonstrate that working with basis-covariant objects substantially simplifies the analysis as compared to working only with $CP$-even or $CP$-odd invariants.

IV Explicit $CP$ conservation in the 3HDM

Let us first remark that in the 3HDM, necessary and sufficient conditions for explicit $CP$ conservation in terms of basis invariants are not known. Of course, one could follow the standard procedure and construct an arbitrary number of $CP$-odd invariants out of $Y_{ab}$ and $Z_{abcd}$ and require them to vanish, Botella:1994cs ; book ; Branco:2005em ; Davidson:2005cw ; Gunion:2005ja ; Varzielas:2016zjc . However, it is not known when one can stop this routine such that it is guaranteed that all higher-power $CP$-odd invariants also vanish. Also, constructing invariants and checking their algebraic independence is a cumbersome task, which is unavoidably delegated to a computer. In addition, as in the 2HDM above, $CP$-odd invariants by themselves do not distinguish physically distinct forms of $CP$ symmetry. Therefore, one is again forced to involve $CP$-even invariants to distinguish different cases. The method proposed in Trautner:2018ipq may help overcome these difficulties, but it must first be extended to $SU(3)$.

It turns out that all of these conditions have been derived using basis-covariant objects, derived in the bilinear formalism. The necessary and sufficient conditions for explicit CP2 conservation in the 3HDM were formulated already in 2006 in Nishi:2006tg . The conditions for CP4 symmetry, as well as for simultaneously occurring CP2 and CP4 symmetries were derived recently in Ivanov:2018ime . All these results are based on the (relative orientation between) vectors $M$ and $L$ and on the eigenvectors and eigenvalues of the $8\times 8$ real symmetric matrix $\Lambda$. We will briefly recapitulate these results in this section.

IV.1 CP2 conservation in the 3HDM

For any CP2 transformation, there exists a basis in which $X$ is the unit matrix and the $CP$ transformation takes the standard form: $\phi_{a}\xrightarrow{CP}\phi_{a}^{*}$, $a=1,2,3$. A necessary and sufficient condition for the potential (2) to be explicitly CP2-conserving is that in this basis all coupling coefficients are real.

In the adjoint space, in the basis where $X$ is the unit matrix, the CP2 transformation leaves invariant all vectors in the 5D subspace $V_{+}=(x_{3},\,x_{8},\,x_{1},\,x_{4},\,x_{6})$, while it flips the sign of all vectors in the 3D subspace $V_{-}=(x_{2},\,x_{5},\,x_{7})$. Therefore, the 3HDM potential is explicitly CP2-invariant if and only if there exists a basis in which the vectors $M,L\in V_{+}$ and $\Lambda$ is block-diagonal with a $5\times 5$ block in $V_{+}$ and a $3\times 3$ block in $V_{-}$.

The challenge then is to formulate this splitting in a basis-invariant form. This problem was solved in Nishi:2006tg with the aid of the eigenvectors of the matrix $\Lambda$. First, consider the $3\times 3$ block of the matrix $\Lambda$ in the 3D subspace $V_{-}$. Using the explicit expressions for the tensor $f_{ijk}$ given in (26), one can verify that vectors of this subspace are closed under the $f$-product defined in (8): if $a,b\in V_{-}$, then $F_{i}=f_{ijk}a_{j}b_{k}\in V_{-}$. Thus, CP2 invariance implies that there exist three mutually orthogonal eigenvectors of $\Lambda$, which we denote $e$, $e^{\prime}$, and $e^{\prime\prime}$, which are closed under the $f$-product.

Next, associating a vector $a$ in the adjoint space to a traceless Hermitian $3\times 3$ matrix $A$ à la $A:=2a_{i}t_{i}$, one can recognize that the $f$-product of $a$ and $b$ corresponds to the commutator of $A$ and $B$ Ivanov:2018ime . Thus, CP2 invariance implies that the three traceless Hermitian matrices $E$, $E^{\prime}$, and $E^{\prime\prime}$, corresponding to the three eigenvectors $e$, $e^{\prime}$, and $e^{\prime\prime}$, are closed under taking commutators. In short, $E$, $E^{\prime}$, and $E^{\prime\prime}$ form a 3D subalgebra of $su(3)$.

There exist two options for 3D subalgebras of $su(3)$: $su(2)$ or $so(3)$. In the adjoint space, the former corresponds, in a suitable basis, to the subspace $(x_{1},x_{2},x_{3})$ while the latter corresponds to $V_{-}$. The difference between them is numerical: $f_{123}=1$ while $f_{257}=1/2$. Therefore, if one finds that three orthonormal eigenvectors of $\Lambda$ are closed under the $f$-product, in order to decide whether one has found the correct subalgebra, one computes the invariant

[TABLE]

If $I_{N}=1$, we have found three eigenvectors which, in a suitable basis, span $V_{-}$. The last step is to check whether or not the vectors $M$ and $L$ have components in $V_{-}$. If not, $M_{i}e_{i}=M_{i}e^{\prime}_{i}=M_{i}e^{\prime\prime}_{i}=0$ and $L_{i}e_{i}=L_{i}e^{\prime}_{i}=L_{i}e^{\prime\prime}_{i}=0$. If and only if all these conditions are satisfied, we have an explicitly CP2 conserving 3HDM Nishi:2006tg .

This sequence of checks represents the necessary and sufficient conditions for a CP2 symmetry in the 3HDM. Despite involving basis-covariant objects, the ultimate conditions are basis-invariant and, therefore, can be checked in any basis. It is clear that there must exist a formulation of these conditions in terms of basis invariants, either $CP$-even or $CP$-odd. However, such a formulation is likely to be extremely complicated as it has not been found yet.

IV.2 CP4 conservation in the 3HDM

Using three Higgs doublets, it is possible to construct a model whose only symmetry is a $CP$-symmetry of order 4 (CP4) Ivanov:2011ae . This model, proposed in Ivanov:2015mwl and denoted CP4 3HDM, has a peculiar property: despite explicit $CP$ conservation, it contains irremovable complex coefficients in the scalar potential. Unlike in the 2HDM, imposing a CP4 symmetry on the 3HDM does not produce any accidental CP2 symmetry. Thus, CP4 opens the path to a new model, physically distinct from any CP2-conserving situation and bearing its own interesting phenomenology Aranda:2016qmp ; Ferreira:2017tvy ; Ivanov:2017bdx ; Haber:2018iwr ; Cherchiglia:2019gll . The basis-invariant necessary and sufficient conditions for the CP4 invariance in the 3HDM are not known in terms of basis invariants. However, they are known in terms of basis-covariant objects Ivanov:2018ime , and we will briefly review them here.

A $CP$ transformation of order 4 is a transformation $\phi_{a}\xrightarrow{CP}X_{ab}\phi_{b}^{*}$, whose matrix $X$, in an appropriate basis, takes the form

[TABLE]

In this basis, CP4 acts on the adjoint space as

[TABLE]

For the potential to be CP4 invariant, $M$ and $L$ must be aligned with $x_{8}$, the only 1D subspace invariant under CP4, while $\Lambda$ must have the block diagonal form

[TABLE]

with an arbitrary $3\times 3$ block in the subspace $(x_{1},x_{2},x_{3})$ and very specific correlation patterns in the $4\times 4$ block of the $(x_{4},x_{5},x_{6},x_{7})$ subspace, cf. Ivanov:2018ime .

The block-diagonal form of $\Lambda$ in Eq. (25) has two key features: the stand-alone direction $x_{8}$ and the 3D block in the subspace $(x_{1},x_{2},x_{3})$. Both features can be detected in a basis-invariant way.

First, one can take the $d$-product (8) of any vector $a$ with itself: $a_{i}\mapsto\sqrt{3}d_{ijk}a_{j}a_{k}$. If the resulting vector is parallel to $a$, we call it self-aligned. If a vector is self-aligned, then there exists a basis in which it lies along the $x_{8}$ direction Ivanov:2018ime . Therefore, the basis-invariant criterion for splitting the $\Lambda_{88}$ entry from the rest is that there exists an eigenvector of $\Lambda$ which is self-aligned. Let us denote this eigenvector as $e^{(8)}$.

Next, the $3\times 3$ block in the subspace $(x_{1},x_{2},x_{3})$ implies that there are three eigenvectors lying in this subspace. Vectors in this subspace have a remarkable property: they are $f$-orthogonal to the previously identified $e^{(8)}$. The converse also holds: if there is a vector $a$ which is both orthogonal to $e^{(8)}$ (meaning $a_{i}e^{(8)}_{i}=0$) and $f$-orthogonal to it (meaning the vector $f_{ijk}a_{j}e^{(8)}_{k}=0$), then it must lie in the $(x_{1},x_{2},x_{3})$ subspace and nowhere else. Thus, we arrive at the basis-invariant condition for the $3\times 3$ block to split from the rest: there must exist three mutually orthogonal eigenvectors of $\Lambda$ which are both orthogonal and $f$-orthogonal to the previously identified $e^{(8)}$.

Once the block-diagonal form of $\Lambda$ is established, what remains to be checked for the CP4 invariance is that all vectors are aligned with $e^{(8)}$. This concerns not only $M$ and $L$ but also vectors $K_{i}=d_{ijk}\Lambda_{jk}$ and $K^{(2)}_{i}=d_{ijk}(\Lambda^{2})_{jk}$. A byproduct of these checks is that the eigenvalues of $\Lambda$ in the $(x_{4},x_{5},x_{6},x_{7})$ subspace are, at least, pairwise degenerate. Having established the block-diagonal form of $\Lambda$ and having checked that $M$, $L$, $K$, and $K^{(2)}$ are aligned with $e^{(8)}$, one concludes that the model respects a CP4 symmetry, see details in Ivanov:2018ime .

IV.3 Combining CP4 and CP2

From the classification of all discrete symmetry-based 3HDMs presented in Ivanov:2012ry ; Ivanov:2012fp , it follows that if CP4 is combined with any other symmetry, the resulting symmetry group will unavoidably include a CP2 transformation. Thus, if one wishes to detect the presence of CP4-symmetry in the 3HDM while at the same time excluding the presence of any other accidental symmetry, it is sufficient to check the absence of any CP2 symmetry.

In principle, one can check the conditions for CP2 independently from CP4. However, using the algorithm explained above it turns out that it is actually shorter to check CP2 in addition to CP4 Ivanov:2018ime . For this, one simply checks whether the $f$-product of eigenvectors of the previously identified $(x_{4},x_{5},x_{6},x_{7})$-subspace is itself an eigenvector of $\Lambda$. If this is not the case, one has found a pure CP4 3HDM model without any accidental symmetry.

The problem of the basis-invariant recognition of an additional CP2 symmetry in a CP4 symmetric 3HDM was also studied recently in Haber:2018iwr . Starting with a CP4 symmetric 3HDM, the authors discovered a $CP$-even basis-invariant, called $\mathcal{N}$, which is zero if and only if the model possesses an additional CP2 symmetry that commutes with CP4. However, the constructed invariant is quite involved, being a high-degree polynomial of the quartic coefficients of the potential and the vacuum expectation values of the doublets. In contrast, the algorithm presented in Ivanov:2018ime is short, transparent, covers both commuting and non-commuting cases, and does not rely on vacuum expectation values. This nicely illustrates the power of working with eigenvectors of $\Lambda$, i.e. basis covariant objects, as compared to manipulations directly with invariants. Since in the end one also wishes to link the basis invariant statements to physical observables, as partly achieved in Haber:2018iwr , it seems likely that a combination of the different techniques could be fruitful.

V Conclusions

Physical observables must be independent of arbitrary choices of basis and, therefore, expressible in terms of basis-invariant combinations of the Lagrangian parameters. In addition, the potentially vast number of basis choices can obscure physically relevant features of a model. Together this justifies efforts to formulate various physically important features of New Physics models, such as the issue of detecting explicit $CP$ conservation in multi-Higgs models, in terms of basis invariants.

However, in many cases it turns out to be prohibitively difficult to corner a physical feature of a model directly in terms of basis invariants. The main message of the present paper is that in such cases it is still possible to establish basis-invariant statements as relations between basis-covariant objects.

Using relations between basis-covariant objects often offers a more direct and transparent way to basis-invariant statements. This is firmly supported by the fact that a recently proposed systematic method for construction of basis invariants uses basis-covariant objects as building blocks Trautner:2018ipq . Resolving the substructure of basis invariants in terms of basis-covariant objects is essential there, also to derive relations between different basis invariants. In the present paper, we have supported this message with many examples of conditions for $CP$ conservation in the two- and three-Higgs-doublet models.

The 2HDM here serves as a warm-up exercise because the conditions are known both in terms of invariants, and in term of basis-covariant objects. Therefore, a “dictionary” could be established. We have shown that expressing various forms of $CP$ conservation via vectors in the adjoint space, including the eigenvectors of the matrix $\Lambda$, leads directly to the basis-invariant conditions.

In the 3HDM, the power of using the (basis-covariant) eigenvectors of $\Lambda$ becomes evident. All cases of necessary and sufficient conditions for the 3HDM scalar potential to be invariant under CP2, CP4, or CP2 and CP4 simultaneously, were derived in this approach. A reformulation of these conditions directly in terms of basis invariants is not yet known.

All these examples support our point that, when dealing with sophisticated models with a large number of basis-choices, it can be much more efficient to work with basis-covariant building blocks rather than with basis invariants directly. These building blocks do not have to be derived from the bilinear approach; alternative paths are conceivable Trautner:2018ipq .

Finally, we stress that the entire approach is in no way limited to our illustrative example of $CP$ violation in multi-Higgs models. Analogous methods can be applied to achieve the basis-invariant detection of other symmetries in other models.

Acknowledgments

We thank João P. Silva for many useful discussions and numerous comments on the paper. I.P.I. acknowledges funding from the Portuguese Fundação para a Ciência e a Tecnologia (FCT) through the FCT Investigator contract IF/00989/2014/CP1214/CT0004 under the IF2014 program, and through the contracts UID/FIS/00777/2013, CERN/FIS-NUC/0010/2015, and PTDC/FIS-PAR/29436/2017, which are partially funded through POCI (FEDER), COMPETE, QREN, and the EU. I.P.I. also acknowledges the support from National Science Center, Poland, via the project Harmonia (UMO-2015/18/M/ST2/00518). C.C.N. acknowledges partial support by Brazilian funding agencies Fapesp through grant 2014/19164-6 and CNPq through grant 308578/2016-3. The work of A.T. is partly supported by a postdoc fellowship of the German Academic Exchange Service (DAAD). A.T. is grateful to Stuart Raby and the Physics Department of Ohio State University for hospitality during the completion of this work.

Appendix A Invariant tensors of $SU(3)$

In this appendix we give explicit expressions for the $SU(3)$ invariant tensors $f_{ijk}$ and $d_{ijk}$ of equation (5), as well as for the $1+8$ gauge-invariant bilinear combinations $(r_{0},r_{i})$ defined in (4). With the usual choice of basis for the Gell-Mann matrices $\lambda_{i}$, the totally anti-symmetric $f_{ijk}$ and symmetric $d_{ijk}$ have the non-zero components

[TABLE]

as well as

[TABLE]

The coefficient in the definition of the $SU(3)$ singlet $r_{0}$ is not fixed by the bilinear construction, but the exact normalization is inessential here. We use the definition adopted from Ivanov:2010ww but alternative normalization factors are possible, see e.g. Maniatis:2014oza . In the Gell-Mann basis, the bilinears $r_{i}$ read

[TABLE]