# Fundamental Composite Higgs Dynamics on the Lattice:

SU(2) with Two Flavors

###### Abstract

In reference Cacciapaglia:2014uja () a unified description, both at the effective and fundamental Lagrangian level, of models of composite Higgs dynamics was proposed. In the unified framework the Higgs itself can emerge, depending on the way the electroweak symmetry is embedded, either as a pseudo-Goldstone boson or as a massive excitation of the condensate.
The most minimal fundamental description consists of an SU(2) gauge theory with two Dirac fermions transforming according to the defining representation of the gauge group. We therefore provide first principle lattice results for the massive spectrum of this theory. We confirm the chiral symmetry breaking phenomenon and determine the lightest spin-one axial and vector masses. The knowledge of the energy scale at which new states will appear at the Large Hadron Collider is of the utmost relevance to guide experimental searches of new physics.

Preprint: CP3-Origins-2014-012 DNRF90 & DIAS-2014-12

## I Introduction

The Standard Model (SM) of particle interactions successfully describes Nature. However, the SM is unappealing. For example the SM Higgs sector simply models spontaneous symmetry breaking, it does not explain it. Furthermore, there is no consistent way to protect the electroweak scale from higher scales, leading to the SM naturalness problem. We refer to Antipin:2013exa () for a mathematical classification of different degrees of naturality.

It is well known that by replacing the SM Higgs sector with a fundamental gauge dynamics featuring fermionic matter fields renders the SM Higgs sector natural. Technicolor Weinberg:1975gm (); Susskind:1978ms () is a time-honored incarnation of this idea. Other ways to use fundamental dynamics to replace the SM Higgs sector appeared later in Kaplan:1983fs (); Kaplan:1983sm (). The Technicolor Higgs Sannino:1999qe (); Hong:2004td (); Dietrich:2005jn (); Dietrich:2005wk (); Sannino:2009za () is the lightest scalar excitation of the fermion condensate responsible for electroweak symmetry breaking. The interplay between the gauge sector and the SM fermion mass sector is relevant because it can reduce the physical mass of the Technicolor Higgs Foadi:2012bb (). If the underlying dynamics has a larger global symmetry group than the one strictly needed to break the electroweak symmetry successfully, one may be able to choose an electroweak embedding in a way that the electroweak symmetry remains intact. Differently from the Technicolor case, here the Higgs state could be identified with one of the Goldstone Bosons (GB) of the theory. In this case the challenges are not only to provide masses to the SM fermions but also to break the electroweak symmetry by means of yet another sector which can also contribute to give mass to the would–be pseudo-GB Higgs. Real progress with respect to the SM Higgs sector shortcomings is achieved, however, only if a more fundamental description exists.

In reference Cacciapaglia:2014uja () a first unified description of models of electroweak composite dynamics was put forward. The description clarified the main similarities, interplay, and shortcomings of the different approaches. In addition a specific underlying realization in terms of fundamental strongly coupled gauge theories was investigated with a clear link to first principle lattice simulations. It was also shown that for a generic electroweak vacuum alignment, the observed Higgs is neither a purely pGB state nor the Technicolor Higgs, but a mixed state. This result has relevant implications for its physical properties and associated phenomenology.

Given a possible underlying gauge theory featuring fermionic matter one can imagine distinct patterns of chiral symmetry breaking Peskin:1980gc (); Preskill:1980mz (); Kosower:1984aw (); Sannino:2004qp (); Dietrich:2006cm (); Sannino:2009aw (); Mojaza:2012zd (). First principle lattice simulations are now in a position to answer these questions Catterall:2007yx (); Catterall:2008qk (); DelDebbio:2008wb (); DelDebbio:2008zf (); Catterall:2009sb (); Hietanen:2009az (); DelDebbio:2009fd (); Kogut:2010cz (); Karavirta:2011zg (); Lewis:2011zb (); Hietanen:2012qd (); Hietanen:2012sz (); Hietanen:2013fya (); Hietanen:2013gva ().

The classification of underlying gauge theories relevant for Technicolor models appeared in Dietrich:2006cm (), while for composite models of the Higgs as a pGB can be found in Mrazek:2011iu (); Bellazzini:2014yua (). In reference Cacciapaglia:2014uja () it was concluded that from the point of view of a fundamental theory with fermionic matter, the minimal scenario to investigate is SU(4)Sp(4) (locally isomorphic to SO(5)), for both a minimal Technicolor as well as composite GB Higgs scenario. The difference being in the way one embeds the electroweak theory within the global flavor symmetry. This pattern of chiral symmetry breaking can be achieved dynamically via an underlying SU()=Sp() gauge theory with 2 Dirac flavors (i.e. four Weyl fermions) transforming according to the fundamental representation of the gauge group.

We will provide here the state-of-the-art lattice results confirming the breaking of the global SU(4) symmetry to Sp(4) (locally isomorphic to SO(5)), first observed in Lewis:2011zb (), via the formation of a non-perturbative fermion condensate, and in addition we will further determine the spin-one spectrum.

## Ii The Lattice Method

In the continuum, the Lagrangian for our technicolor template is

(1) |

which can be discretized in the familiar way to arrive at a Wilson action,

(2) | |||||

where is the gauge field and the gauge coupling in conventional lattice notation. is the doublet of and fermions, and is the 22 diagonal mass matrix.

Mesons will couple to local operators of the form

(3) | |||||

(4) | |||||

(5) |

where denotes any product of Dirac matrices. Baryons (which are diquarks in this two-color theory) will couple to local operators of the form

(6) | |||||

(7) | |||||

(8) |

where the Pauli structure acts on color indices while the charge conjugation operator acts on Dirac indices.

We extract the meson masses from the two-point correlation functions

(9) |

where . The quantities of interest are pseudoscalar , vector (), and axial vector mesons. As a source vector we use single time slice stochastic sources Boyle:2008rh ().

In addition to the meson spectrum we are interested in two other quantities, the quark mass and the Goldstone boson decay constant . We define the quark mass through the Partially Conserved Axial Current (PCAC) relation:

(10) |

where

(11) |

The Goldstone boson decay constant can be calculated as:

(12) |

where is obtained from the asymptotic form of at large :

(13) |

To convert the lattice quantities to physical units, one should determine the lattice spacing for our simulations and the appropriate (mass-independent) renormalization constant. The lattice spacing, in a Technicolor model, is fixed by the requirement that the (renormalized) Goldstone boson decay constant has the value of 246 GeV, giving the correct mass to the electroweak gauge bosons. For the general composite Higgs scenario the electroweak decay constant becomes with the depending on the specific electroweak embedding. The actual value of depends on the electroweak quantum corrections, the top corrections as well as the effects of other possible sources of explicit breaking of the initial SU() symmetry. The Technicolor limit is recovered for while the composite pGB Higgs case corresponds to small, but non-vanishing, . Any other value of the is allowed and corresponds to a combination of these two limits. For the details we refer to Cacciapaglia:2014uja (). For definitiveness we present the results for but we will reinstate the dependence on for the spectrum.

The relevant renormalization constant for , commonly denoted in the literature by , has not been computed non-perturbatively for our simulations. In this work we use the perturbative value which has been calculated in DelDebbio:2008wb (). For fermions in the fundamental representation we have:

(14) |

## Iii The Lattice Results

The lattice simulations used in this work extend the results already published in Ref. Lewis:2011zb (). In particular, we have used larger volumes for the set of parameters closest to the chiral limit. The bare parameters used for our simulations are listed in Table 1, where we also report the number of thermalized trajectories, of length one, used in our analysis below. Thermalizations are estimated by monitoring the average plaquette expectation value and the value of the two-point correlation function in the pseudoscalar channel at a time separation of twelve time slices. These two quantities are shown in Fig. 1 for two representative light quark masses on the finest lattices used in this work.

Volume | Therm. | Conf. | ||
---|---|---|---|---|

2.0 | -0.85, -0.9, -0.94, -0.945, -0.947, -0.949 | 320 | 680 | |

2.0 | -0.947 | 500 | 680 | |

2.2 | -0.60, -0.65, -0.68, -0.70, -0.72, -0.75 | 320 | 680 | |

2.2 | -0.75 | 500 | 2000 | |

2.2 | -0.72,-0.735, -0.75 | 500 | 2000 |

Volume | |||||||
---|---|---|---|---|---|---|---|

2.0 | -0.85 | 0.1919(6) | 0.9163(18) | 1.008(20) | 1.64(3) | 0.1544(4) | |

2.0 | -0.9 | 0.1134(6) | 0.708(3) | 0.821(3) | 1.47(3) | 0.1248(5) | |

2.0 | -0.94 | 0.0476(8) | 0.451(7) | 0.636(5) | 1.15(4) | 0.089(9) | |

2.0 | -0.945 | 0.038(7) | 0.407(5) | 0.57(6) | 1.08(3) | 0.0799(7) | |

2.0 | -0.947 | 0.0327(7) | 0.377(6) | 0.546(7) | 1.02(4) | 0.0754(8) | |

2.0 | -0.949 | 0.0307(8) | 0.374(6) | 0.546(9) | 1.02(3) | 0.075(12) | |

2.0 | 0.947 | 0.0309(3) | 0.3739(14) | 0.536(4) | 1.0(5) | 0.0766(6) | |

2.2 | -0.6 | 0.2296(7) | 0.886(3) | 0.93(3) | 1.371(12) | 0.1119(5) | |

2.2 | -0.65 | 0.1637(7) | 0.755(3) | 0.792(3) | 1.229(9) | 0.0992(5) | |

2.2 | -0.68 | 0.1205(7) | 0.612(4) | 0.671(5) | 1.068(17) | 0.0868(5) | |

2.2 | -0.7 | 0.0968(7) | 0.548(5) | 0.615(6) | 1.018(12) | 0.0793(5) | |

2.2 | -0.72 | 0.0686(6) | 0.455(5) | 0.531(5) | 0.884(19) | 0.0684(5) | |

2.2 | -0.75 | 0.0264(8) | 0.324(8) | 0.445(9) | 0.76(3) | 0.0405(12) | |

2.2 | -0.75 | 0.024(5) | 0.258(4) | 0.359(8) | 0.62(4) | 0.0433(9) | |

2.2 | -0.72 | 0.0661(5) | 0.4475(8) | 0.522(14) | 0.81(3) | 0.0664(5) | |

2.2 | -0.735 | 0.0456(3) | 0.3612(18) | 0.446(4) | 0.75(3) | 0.0568(5) | |

2.2 | -0.75 | 0.0257(5) | 0.2649(16) | 0.363(5) | 0.59(6) | 0.0457(7) |

All the ensembles of gauge configurations were created using the GPU version of the HiRep code DelDebbio:2008zf (). The lattice action used is the plaquette-action SU(2) gauge theory with two flavors ( and ) of mass-degenerate Wilson fermions. The Hybrid Monte Carlo trajectory length was chosen to be one. The autocorrelation times for plaquette expectation values and meson correlators were estimated to be about 10. The errors for all quantities extracted in this work were obtained using a bootstrap procedure.

In a previous work by some of the authors Lewis:2011zb (), a first estimate of the Goldstone spectrum was already obtained. However large finite volume effects were observed – see Fig. 5(b) of Lewis:2011zb () – for the lightest value of the quark mass on the finest lattice used in that work, corresponding to the bare parameter couplings .

Here we perform a more systematic analysis to control finite volume effects, using simulations on three different lattice volumes , , and at the lightest quark mass on the finest lattice. At this quark mass, the results for the smallest volume suffer clearly from finite volume effects, whereas observables measured on the two largest lattices agree within statistical errors, as shown in Fig. 2. From the size of the statistical errors on the two largest lattices, we estimate that the residual finite volume effects at our lightest quark mass are below 2% for , and and below 10% for .

To confirm spontaneous chiral symmetry breaking one should, in principle, reach the chiral regime of the theory, pushing the quark masses light enough that chiral perturbation theory (PT), or the appropriate lattice extension of it, could be used, while keeping under control all other systematic sources of error which are present on the lattice. It is well known, by studies of QCD, that this chiral regime is extremely difficult to reach as lattice artifacts and residual finite volume effects tend to make the predictions of PT difficult to test.

Keeping this in mind, we analyze our data for signs of spontaneous chiral symmetry breaking and check how well predictions from PT fit the measured Goldstone spectrum. The Goldstone boson mass and decay constant are studied as a function of quark mass, defined through the PCAC relation and compared with the expectations from (continuum) chiral perturbation theory at next to leading order (NLO):

(15) |

and

(16) |

where , , and are (unknown) low-energy constants of the theory, and and are known constants. For our theory and . Additional terms in the chiral expansion can also be computed for the NNLO approximations. The relevant expressions can be found in Bijnens:2009qm (), from which the quoted values for and were taken.

Correction to the continuum chiral expansion arise due to the lattice discretisation at non-zero lattice spacing. For Wilson fermions at NLO the functional form of Eqs. (15) and (16) remains unchanged, but the coefficients, and in particular and , depend on the lattice spacing . In the limit one should recover the continuum values for and , however for a fixed lattice spacing and are two additional free parameters of the expansion.

We show in Figs. 3 and 4 our results for and for the two values of the lattice spacing used in this work. It is possible to use NLO Wilson chiral perturbation theory to fit our data at small quark masses for both and . We report in Table 3 the quark mass ranges, and the coefficients and for the fits for the two different values of the lattice spacing used in this work. The relative errors on the fitting parameters are large especially for the coefficients and of the terms of the chiral expansion which suffer from very large uncertainties , and are thus compatible with zero. This can be explained by our data not being yet in a regime where the terms can be clearly distinguished from the polynomial terms in the expansion, even at the lightest quark masses available.

range | /dof | dof | LO coefficient | ||
---|---|---|---|---|---|

2.0 | 0.48 | 2 | |||

2.0 | 2.28 | 2 | |||

2.2 | 0.32 | 1 | |||

2.2 | 2.01 | 1 |

Given that the chiral logs are subdominant, a simple polynomial fit to the data is expected to be an adequate description of and . We therefore fit our data setting . The values for and thus obtained are given in Table 4. The central values of and are compatible, within statistical errors, with the ones obtained using NLO Wilson chiral perturbation theory. In the analysis below we use these values of and as our best estimate, and consider the errors from different fitting procedures as systematic errors. In the last column of Table 4, we use the perturbative value of from Eq. (14) to obtain the renormalized .

We also note that our data are not well described by NLO or NNLO continuum chiral perturbation theory, i.e. when the coefficients and of the logarithmic terms are fixed. In this case we can perform a simultaneous fit of both and , and the fit is thus much more constrained. The resulting shows that this is not a good description of our data.

2.0 | 2.52(12) | 0.052(3) | 0.85 | 0.0439(18) |
---|---|---|---|---|

2.2 | 1.26(03) | 0.033(1) | 0.86 | 0.0285(08) |

The values of vector and axial vector meson masses, as measured from our emsemble of configurations, are plotted in Fig. 5. A linear function represents the data well at small quark masses .

Given that , and are well described by a linear function at small quark masses, one can expect that the two ratios and are also linear functions of the quark mass close to the chiral limit. These ratios are shown in Fig. 6, together with a linear extrapolation to zero quark mass. We will refer to this method of chiral extrapolation as “method 2”, whereas the first method will be named “method 1” in the following.

We use these two different methods as a crosscheck of the chiral extrapolation and to try to quantify the systematic errors due to the choice of extrapolation function. We compare in Table 5 the results in lattice units for the chiral extrapolations. The methods are clearly consistent with each other and method 2 leads to an overall smaller statistical error.

Combining the data from the two lattice spacings available in this study, we can perform a first, crude continuum extrapolation for the masses of the vector and axial vector mesons. As explained above, the lattice spacing is fixed by the requirement that the value of the renormalized Goldstone decay constant satisfies , as required to give the correct masses to the electroweak gauge bosons. For concreteness here we assume , but the dependence on can easily be reinstated when required as done below. The results of the linear extrapolations of and to the continuum limit are reported in Table 6.

Our continuum extrapolation is subject to two major sources of systematic errors. First our simulations are performed only at two lattice spacings, and therefore we do not have a good measure of how well our linear extrapolation describes the data. To take this into account, as a systematic error we quote, quite conservatively, the difference between the value of the continuum extrapolated value and the data point of the finer lattice. The second systematic error stems from the renormalization constant which we do not measure non-perturbatively. As a systematic error we then use the difference between the perturbative value of and . In Table 6 we list the continuum extrapolated values of vector and axial vector mesons for both methods 1 and 2 described above for chiral extrapolation. The vector case is plotted in Fig. 7.

The results produced by both methods are comparable and well inside each other’s error bars. As final results for the meson masses we quote the one obtained by method 2. Square summing the errors, the vector meson reads and the axial vector meson where we have reinstated the dependence on the angle defining the specific electroweak embedding.

Method 1 | Method 2 | ||
---|---|---|---|

2.0 | 8.1(5) | 8.65(8) | |

2.2 | 9.3(4) | 9.22(11) | |

2.0 | 18(2) | 16.6(4) | |

2.2 | 17(3) | 15.5(6) |

Method 1 (GeV) | Method 2 (GeV) | ||
---|---|---|---|

2840(330)(560)(360) | 2520(100)(240)(310) | ||

4000(1800)(200)(430) | 3300(400)(510)(340) |

## Iv Conclusions

The SU(2)-gauge theory with two fundamental fermions unifies both Technicolor and composite pGB Higgs models of electroweak symmetry breaking. In this work, we have calculated the masses of the two lightest non-singlet mesons using the Goldstone boson decay constant to set the scale. We performed the calculations with two different lattice spacings. With conservative error estimates the mass of the lightest vector meson is . This value is clearly above one TeV and outside the current exclusion limits set by the LHC Aad:2012hf ().

To increase the precision of our results for the spectrum, at least one additional lattice spacing is required alongside a nonperturbative determination of the renormalization constant . This would require a significant increase of computational resources. Furthermore we are eager to investigate the scalar sector and the vector decay constants.

###### Acknowledgements.

This work was supported by the Danish National Research Foundation DNRF:90 grant, by a Lundbeck Foundation Fellowship grant, and by NSERC of Canada. The computing facilities were provided by the Danish Centre for Scientific Computing and Canada’s Shared Hierarchical Academic Research Computing Network (SHARCNET: http://www.sharcnet.ca).## References

- (1) G. Cacciapaglia and F. Sannino, arXiv:1402.0233 [hep-ph]. To appear in JHEP.
- (2) O. Antipin, M. Mojaza and F. Sannino, arXiv:1310.0957 [hep-ph]. To appear in Phys. Rev. D.
- (3) S. Weinberg, Phys. Rev. D 13, 974 (1976).
- (4) L. Susskind, Phys. Rev. D 20, 2619 (1979).
- (5) D. B. Kaplan and H. Georgi, Phys. Lett. B 136, 183 (1984).
- (6) D. B. Kaplan, H. Georgi and S. Dimopoulos, Phys. Lett. B 136, 187 (1984).
- (7) F. Sannino and J. Schechter, Phys. Rev. D 60, 056004 (1999) [hep-ph/9903359].
- (8) D. K. Hong, S. D. H. Hsu and F. Sannino, Phys. Lett. B 597, 89 (2004) [hep-ph/0406200].
- (9) D. D. Dietrich, F. Sannino and K. Tuominen, Phys. Rev. D 72, 055001 (2005) [hep-ph/0505059].
- (10) D. D. Dietrich, F. Sannino and K. Tuominen, Phys. Rev. D 73, 037701 (2006) [hep-ph/0510217].
- (11) F. Sannino, Acta Phys. Polon. B 40, 3533 (2009) [arXiv:0911.0931 [hep-ph]].
- (12) R. Foadi, M. T. Frandsen and F. Sannino, Phys. Rev. D 87, 095001 (2013) [arXiv:1211.1083 [hep-ph]].
- (13) M. E. Peskin, Nucl. Phys. B 175, 197 (1980).
- (14) J. Preskill, Nucl. Phys. B 177, 21 (1981).
- (15) D. A. Kosower, Phys. Lett. B 144, 215 (1984).
- (16) F. Sannino and K. Tuominen, Phys. Rev. D 71, 051901 (2005) [hep-ph/0405209].
- (17) D. D. Dietrich and F. Sannino, Phys. Rev. D 75, 085018 (2007) [hep-ph/0611341].
- (18) F. Sannino, Phys. Rev. D 79, 096007 (2009) [arXiv:0902.3494 [hep-ph]].
- (19) M. Mojaza, C. Pica, T. A. Ryttov and F. Sannino, Phys. Rev. D 86, 076012 (2012) [arXiv:1206.2652 [hep-ph]].
- (20) S. Catterall and F. Sannino, Phys. Rev. D 76, 034504 (2007) [arXiv:0705.1664 [hep-lat]].
- (21) S. Catterall, J. Giedt, F. Sannino and J. Schneible, JHEP 0811, 009 (2008) [arXiv:0807.0792 [hep-lat]].
- (22) L. Del Debbio, M. T. Frandsen, H. Panagopoulos and F. Sannino, JHEP 0806, 007 (2008) [arXiv:0802.0891 [hep-lat]].
- (23) L. Del Debbio, A. Patella and C. Pica, Phys. Rev. D 81, 094503 (2010) [arXiv:0805.2058 [hep-lat]].
- (24) S. Catterall, J. Giedt, F. Sannino and J. Schneible, arXiv:0910.4387 [hep-lat].
- (25) A. J. Hietanen, K. Rummukainen and K. Tuominen, Phys. Rev. D 80, 094504 (2009) [arXiv:0904.0864 [hep-lat]].
- (26) L. Del Debbio, B. Lucini, A. Patella, C. Pica and A. Rago, Phys. Rev. D 80, 074507 (2009) [arXiv:0907.3896 [hep-lat]].
- (27) J. B. Kogut and D. K. Sinclair, Phys. Rev. D 81, 114507 (2010) [arXiv:1002.2988 [hep-lat]].
- (28) T. Karavirta, J. Rantaharju, K. Rummukainen and K. Tuominen, JHEP 1205, 003 (2012) [arXiv:1111.4104 [hep-lat]].
- (29) R. Lewis, C. Pica and F. Sannino, Phys. Rev. D 85, 014504 (2012) [arXiv:1109.3513 [hep-ph]].
- (30) A. Hietanen, C. Pica, F. Sannino and U. I. Sondergaard, PoS LATTICE 2012, 065 (2012) [arXiv:1211.0142 [hep-lat]].
- (31) A. Hietanen, C. Pica, F. Sannino and U. I. Sondergaard, Phys. Rev. D 87, no. 3, 034508 (2013) [arXiv:1211.5021 [hep-lat]].
- (32) A. Hietanen, R. Lewis, C. Pica and F. Sannino, arXiv:1308.4130 [hep-ph]. Extended version to appear.
- (33) A. Hietanen, C. Pica, F. Sannino and U. Søndergaard, arXiv:1311.3841 [hep-lat].
- (34) J. Mrazek, A. Pomarol, R. Rattazzi, M. Redi, J. Serra and A. Wulzer, Nucl. Phys. B 853 (2011) 1 [arXiv:1105.5403 [hep-ph]].
- (35) B. Bellazzini, C. Csáki and J. Serra, arXiv:1401.2457 [hep-ph].
- (36) P. A. Boyle, A. Juttner, C. Kelly and R. D. Kenway, JHEP 0808, 086 (2008) [arXiv:0804.1501 [hep-lat]].
- (37) J. Bijnens and J. Lu, JHEP 0911, 116 (2009) [arXiv:0910.5424 [hep-ph]].
- (38) J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 465 (1985).
- (39) G. Aad et al. [ATLAS Collaboration], JHEP 1211, 138 (2012) [arXiv:1209.2535 [hep-ex]].