Here's how their interactions work. I'll be focusing on elementary fermions and Higgs particles, because that's where most of the free parameters in the Standard Model are. Gauge interactions can be completely predicted from the symmetry group and multiplet structure, to within an overall coupling-constant or "charge" value.
For the
Standard Model with neutrino seesaw masses, one has
y
u,ijH'.Q
i.U*
j + y
d,ijH.Q
i.D*
j + y
n,ijH'.L
i.N*
j + y
e,ijH.L
i.E*
j + m
N,ij.N*
i.N*
j + (Hermitian conjugates) + m
H2*|H|
2 + g
H*|H|
4
H = {H
c,H
n} (charged and neutral Higgs, both complex-valued fields), H' = flipped version = {H
n*,-H
c*}
m
H is the unbroken-electroweak Higgs mass, and it's imaginary, making the Higgs particle tachyonic. That makes it unstable, but the instability is limited by its self-interaction, with parameter g
H. That's what makes electroweak symmetry breaking, the nonzero vacuum value of the Higgs field -- {0,v} -- and Standard Model particles' masses.
In the
MSSM, both the EF's and the Higgs particles are Wess-Zumino multiplets, the simplest kind of 4D supersymmetric field. WZ multiplets are composed of SUSY-related spin-0 and spin-1/2 particles. Their interactions look much like the SM ones:
y
u,ijHu.Q
i.U*
j + y
d,ijHd.Q
i.D*
j + y
n,ijHu.L
i.N*
j + y
e,ijHd.L
i.E*
j + m
HudHu.Hd + m
N,ij.N*
i.N*
j + (Hermitian conjugates) + (oodles of SUSY-breaking terms)
When one expands it out, one gets Higgs masses and self-interactions. The Higgs particles Hu and Hd are separate, unlike in the SM -- that's a consequence of SUSY and their WZ nature.
The y's are 3*3 matrices of dimensionless complex quantities, giving 18 free parameters each, and a total of 72. m
N is a 3*3 symmetric complex mass matrix, giving 12 free parameters. There is only one of m
Hud, and it can be complex. However, redefining the fields can easily absorb a large fraction of these parameters, and one gets quark and lepton masses and mixing angles.
In the Next to MSSM, or NMSSM, m
HudHu.Hd gets replaced with g
HSS.Hu.Hd with a Higgs singlet, S, and coupling constant g
HS. Like the right-handed neutrino N, S is a gauge singlet, and S is also a gauge singlet in SU(5) and SO(10) GUT's.
Let's see what happens in GUT's.
In
SU(5), the interactions become
y
u,ijHu.F(10)
i.F(10)
j + y
d,ijHd.F(10)
i.F(5*)
j + y
n,ijHn.F(1)
i.F(5*)
j + m
HudHu.Hd (or g
HSS.Hu.Hd) + m
N,ijF(1)
i.F(1)
j + (Hermitian conjugates)
where F is the EF's, y
e,ij = y
d,ji, and y
u,ij is symmetric.
One gets bottom-tau Higgs-coupling unification, and thus mass unification, and extrapolation to GUT energies in the MSSM can get close masses. A calculation I've found is
[1206.5909] Updated values of running quark and lepton masses at GUT scale in SM, 2HDM and MSSM, m(tau) ~ (1.2 or 1.3) * m(bottom) at GUT energies. Unification is not nearly as good for the lighter generations.
Still a lot of terms, however.
Now,
SO(10). We get
y
ijH.F
i.F
j + m
HudH.H (or g
HSS.H.H) + (Hermitian conjugates)
where F is the EF's and H contains both MSSM Higgs doublets. One gets Higgs-coupling unification: y is symmetric, and y
u, y
d, y
n, and y
e are all equal to it. Complete with no cross-generation mixing. It must therefore be the result of breaking of SO(10) symmetry.
The right-handed-neutrino mass term drops out. It must also be produced by that symmetry breaking.
Interestingly, these interaction terms are the only terms with F and H that are possible with nonnegative mass dimension. Negative mass dimension gives bad behavior for energy scales above the interaction's mass scale. Historically, that was a problem for the original theory of the weak interactions, Fermi's four-fermion contact interaction. That was eventually solved by electroweak unification.
Turning to
E6, we get even more simplification.
y
ijkX
iX
jX
k + (Hermitian conjugates)
where y is symmetric, and X = the 27 multiplet: (elementary fermions: F) + (Higgs doublets: H) + (Higgs singlet: S)
Here again, these terms are the only ones possible with nonnegative mass dimension.
Going from E6 to SO(10), the possible interactions are H.F.F and S.H.H. No m
HudH.H term, but instead the NMSSM term g
HSS.H.H.
So E6 gets it right.
Going to
E8, it has only one free parameter, the gauge coupling constant, and
all the interactions must result from its self-interactions.
So we get more and more simple as we go SM - MSSM - SU(5) - SO(10) - E6 - E8. The complexity that we find is due to symmetry breaking -- a *lot* of it -- and SUSY and GUT symmetry breaking aren't very well-understood.