
The matter fields interaction terms are:=matterinteraction

 

                                                  2
Matterinteraction=(8*el(-1)*er(1)*hig(1)*c(1,4)*wm

                                                         2
                      + 8*hig(-1)*er(-1)*el(1)*c(-1,4)*wm

                                                       2   2
                      - hig(-1)*hig(1)*hig(-1)*hig(1)*g *hm

                                           2   2
                      + 4*hig(-1)*hig(1)*hm *wm

                                                            2
                      + 8*hig(-1)*q1dr(-1)*q1l(1)*c(-1,3)*wm

                                                               2
                      + 8*higtid(-1)*q1ur(-1)*q1l(1)*c(-1,2)*wm

                                                          2
                      + 8*q1l(-1)*q1dr(1)*hig(1)*c(1,3)*wm

                                                             2       2
                      + 8*q1l(-1)*q1ur(1)*higtid(1)*c(1,2)*wm )/(8*wm

                     )

 

                        2   2
                      hm *wm
The vancum energy is=---------
                          2
                       2*g

 
 ********* The following is equation list ********
 

 0=0

 0=0

 0=0

 0=0

 0=0

 0=0

 0=0

 0=0

 0=0

 0=0

 0=0

 0=0

 0=0

 0=0

 
 ********* The following is mass list ********
 

 gsg=0

        wm
 z=------------
    cos(theta)

 gs=0

 
 ********* The following is the propagtors list ********
 

 

 gs' progaptor := - v.va

 

 w' progaptor :=r1.v*r1.va*dmas(wm,2) - v.va

 

                                     wm
 z' progaptor :=r1.v*r1.va*dmas(------------,2) - v.va
                                 cos(theta)

 

 p' progaptor := - v.va

 

 gsg' progaptor :=1

 

 h0' progaptor :=1

 

 qt' progaptor := - (g(l,r1) + fmt)

 

 qc' progaptor := - (g(l,r1) + fmc)

 

 qu' progaptor := - (g(l,r1) + fmu)

 

 tau' progaptor := - (g(l,r1) + fmtau)

 

 mu' progaptor := - (g(l,r1) + fmmu)

 

 ef' progaptor := - (g(l,r1) + fme)

 

 qb' progaptor := - (g(l,r1) + fmb)

 

 qs' progaptor := - (g(l,r1) + fms)

 

 qd' progaptor := - (g(l,r1) + fmd)

 

 nut' progaptor := - g(l,r1)

 

 numu' progaptor := - g(l,r1)

 

 nue' progaptor := - g(l,r1)

 
 *********** The following is the vertex list ***********
 
 ***** The following is the 4-particle vertices ****
 
   Vertex 47: h0(0,r1)---h0(0,r2)---h0(0,r3)---h0(0,r4)---
                   2   2
              - 3*g *hm *i
    vvv(47)=---------------
                     2
                 4*wm

 
   Vertex 8: z(0,r3,vb)---z(0,r4,vc)---w(-1,r2,va)---w(1,r1,
v)---
                     2  2
    vvv(8)=cos(theta) *g *i

           *( - 2*v.va*vb.vc + v.vb*va.vc + v.vc*va.vb)

 
   Vertex 7: z(0,r3,vb)---p(0,r4,vc)---w(-1,r2,va)---w(1,r1,
v)---
                                  2
    vvv(7)=cos(theta)*sin(theta)*g *i

           *( - 2*v.va*vb.vc + v.vb*va.vc + v.vc*va.vb)

 
   Vertex 6: p(0,r3,vb)---p(0,r4,vc)---w(-1,r2,va)---w(1,r1,
v)---
                     2  2
    vvv(6)=sin(theta) *g *i

           *( - 2*v.va*vb.vc + v.vb*va.vc + v.vc*va.vb)

 
   Vertex 5: w(1,r1,v)---w(1,r2,va)---w(-1,r3,vb)---w(-1,r4,
vc)---
            2
    vvv(5)=g *i*(2*v.va*vb.vc - v.vb*va.vc - v.vc*va.vb)

 
   Vertex 2: z(0,r3,v)---z(0,r4,va)---h0(0,r2)---h0(0,r1)---
                    2
              v.va*g *i
    vvv(2)=---------------
                        2
            2*cos(theta)

 
   Vertex 1: h0(0,r1)---h0(0,r2)---w(1,r3,v)---w(-1,r4,va
)---
                  2
            v.va*g *i
    vvv(1)=-----------
                2
 ***** The following is the 3-particle vertices ****
 
   Vertex 54: qt(-2/3,r1)---qt(2/3,r2)---h0(0,r3)---
             fmt*g*i
    vvv(54)=---------
              2*wm

 
   Vertex 53: qc(-2/3,r1)---qc(2/3,r2)---h0(0,r3)---
             fmc*g*i
    vvv(53)=---------
              2*wm

 
   Vertex 52: qu(-2/3,r1)---qu(2/3,r2)---h0(0,r3)---
             fmu*g*i
    vvv(52)=---------
              2*wm

 
   Vertex 51: qb(1/3,r1)---qb(-1/3,r2)---h0(0,r3)---
             fmb*g*i
    vvv(51)=---------
              2*wm

 
   Vertex 50: qs(1/3,r1)---qs(-1/3,r2)---h0(0,r3)---
             fms*g*i
    vvv(50)=---------
              2*wm

 
   Vertex 49: qd(1/3,r1)---qd(-1/3,r2)---h0(0,r3)---
             fmd*g*i
    vvv(49)=---------
              2*wm

 
   Vertex 48: h0(0,r1)---h0(0,r2)---h0(0,r3)---
                      2
              - 3*g*hm *i
    vvv(48)=--------------
                 2*wm

 
   Vertex 46: tau(1,r1)---tau(-1,r2)---h0(0,r3)---
             fmtau*g*i
    vvv(46)=-----------
               2*wm

 
   Vertex 45: mu(1,r1)---mu(-1,r2)---h0(0,r3)---
             fmmu*g*i
    vvv(45)=----------
               2*wm

 
   Vertex 44: ef(1,r1)---ef(-1,r2)---h0(0,r3)---
             fme*g*i
    vvv(44)=---------
              2*wm

 
   Vertex 43: z(0,r3,v)---qt(2/3,r2)---qt(-2/3,r1)---
    vvv(43)=(g*i

                                               2
             *(3*g(l,a,v) + 8*g(l,v)*cos(theta)  - 5*g(l,v))

             )/(12*cos(theta))

 
   Vertex 42: z(0,r3,v)---qc(2/3,r2)---qc(-2/3,r1)---
    vvv(42)=(g*i

                                               2
             *(3*g(l,a,v) + 8*g(l,v)*cos(theta)  - 5*g(l,v))

             )/(12*cos(theta))

 
   Vertex 41: z(0,r3,v)---qu(2/3,r2)---qu(-2/3,r1)---
    vvv(41)=(g*i

                                               2
             *(3*g(l,a,v) + 8*g(l,v)*cos(theta)  - 5*g(l,v))

             )/(12*cos(theta))

 
   Vertex 40: p(0,r3,v)---qt(2/3,r2)---qt(-2/3,r1)---
             2*sin(theta)*g*i*g(l,v)
    vvv(40)=-------------------------
                        3

 
   Vertex 39: p(0,r3,v)---qc(2/3,r2)---qc(-2/3,r1)---
             2*sin(theta)*g*i*g(l,v)
    vvv(39)=-------------------------
                        3

 
   Vertex 38: p(0,r3,v)---qu(2/3,r2)---qu(-2/3,r1)---
             2*sin(theta)*g*i*g(l,v)
    vvv(38)=-------------------------
                        3

 
   Vertex 37: qu(-2/3,r1)---qd(-1/3,r2)---w(1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(37)=---------------------------------
                            4

 
   Vertex 36: qc(-2/3,r1)---qs(-1/3,r2)---w(1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(36)=---------------------------------
                            4

 
   Vertex 35: qt(-2/3,r1)---qb(-1/3,r2)---w(1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(35)=---------------------------------
                            4

 
   Vertex 34: qd(1/3,r1)---qu(2/3,r2)---w(-1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(34)=---------------------------------
                            4

 
   Vertex 33: qs(1/3,r1)---qc(2/3,r2)---w(-1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(33)=---------------------------------
                            4

 
   Vertex 32: qb(1/3,r1)---qt(2/3,r2)---w(-1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(32)=---------------------------------
                            4

 
   Vertex 31: z(0,r3,v)---qb(-1/3,r2)---qb(1/3,r1)---
                                                     2
    vvv(31)=(g*i*( - 3*g(l,a,v) + 4*g(l,v)*sin(theta)

                 - 3*g(l,v)))/(12*cos(theta))

 
   Vertex 30: z(0,r3,v)---qs(-1/3,r2)---qs(1/3,r1)---
                                                     2
    vvv(30)=(g*i*( - 3*g(l,a,v) + 4*g(l,v)*sin(theta)

                 - 3*g(l,v)))/(12*cos(theta))

 
   Vertex 29: z(0,r3,v)---qd(-1/3,r2)---qd(1/3,r1)---
                                                     2
    vvv(29)=(g*i*( - 3*g(l,a,v) + 4*g(l,v)*sin(theta)

                 - 3*g(l,v)))/(12*cos(theta))

 
   Vertex 28: p(0,r3,v)---qb(-1/3,r2)---qb(1/3,r1)---
              - sin(theta)*g*i*g(l,v)
    vvv(28)=--------------------------
                        3

 
   Vertex 27: p(0,r3,v)---qs(-1/3,r2)---qs(1/3,r1)---
              - sin(theta)*g*i*g(l,v)
    vvv(27)=--------------------------
                        3

 
   Vertex 26: p(0,r3,v)---qd(-1/3,r2)---qd(1/3,r1)---
              - sin(theta)*g*i*g(l,v)
    vvv(26)=--------------------------
                        3

 
   Vertex 25: z(0,r3,v)---w(-1,r2,vb)---w(1,r1,va)---
    vvv(25)=cos(theta)*g*i*( - r1.v*va.vb + r1.vb*v.va

                + r2.v*va.vb - r2.va*v.vb + r3.va*v.vb

                - r3.vb*v.va)

 
   Vertex 24: p(0,r3,v)---w(-1,r2,vb)---w(1,r1,va)---
    vvv(24)=sin(theta)*g*i*( - r1.v*va.vb + r1.vb*v.va

                + r2.v*va.vb - r2.va*v.vb + r3.va*v.vb

                - r3.vb*v.va)

 
   Vertex 23: z(0,r3,v)---nut(0,r2)---nut(0,r1)---
             g*i*(g(l,a,v) + g(l,v))
    vvv(23)=-------------------------
                  4*cos(theta)

 
   Vertex 22: z(0,r3,v)---numu(0,r2)---numu(0,r1)---
             g*i*(g(l,a,v) + g(l,v))
    vvv(22)=-------------------------
                  4*cos(theta)

 
   Vertex 21: z(0,r3,v)---nue(0,r2)---nue(0,r1)---
             g*i*(g(l,a,v) + g(l,v))
    vvv(21)=-------------------------
                  4*cos(theta)

 
   Vertex 20: nue(0,r1)---ef(-1,r2)---w(1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(20)=---------------------------------
                            4

 
   Vertex 19: numu(0,r1)---mu(-1,r2)---w(1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(19)=---------------------------------
                            4

 
   Vertex 18: nut(0,r1)---tau(-1,r2)---w(1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(18)=---------------------------------
                            4

 
   Vertex 17: ef(1,r1)---nue(0,r2)---w(-1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(17)=---------------------------------
                            4

 
   Vertex 16: mu(1,r1)---numu(0,r2)---w(-1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(16)=---------------------------------
                            4

 
   Vertex 15: tau(1,r1)---nut(0,r2)---w(-1,r3,v)---
             sqrt(2)*g*i*(g(l,a,v) + g(l,v))
    vvv(15)=---------------------------------
                            4

 
   Vertex 14: z(0,r3,v)---tau(-1,r2)---tau(1,r1)---
    vvv(14)=(g*i

                                                2
             *( - g(l,a,v) + 4*g(l,v)*sin(theta)  - g(l,v)))

            /(4*cos(theta))

 
   Vertex 13: z(0,r3,v)---mu(-1,r2)---mu(1,r1)---
    vvv(13)=(g*i

                                                2
             *( - g(l,a,v) + 4*g(l,v)*sin(theta)  - g(l,v)))

            /(4*cos(theta))

 
   Vertex 12: z(0,r3,v)---ef(-1,r2)---ef(1,r1)---
    vvv(12)=(g*i

                                                2
             *( - g(l,a,v) + 4*g(l,v)*sin(theta)  - g(l,v)))

            /(4*cos(theta))

 
   Vertex 11: p(0,r3,v)---tau(-1,r2)---tau(1,r1)---
    vvv(11)= - sin(theta)*g*i*g(l,v)

 
   Vertex 10: p(0,r3,v)---mu(-1,r2)---mu(1,r1)---
    vvv(10)= - sin(theta)*g*i*g(l,v)

 
   Vertex 9: p(0,r3,v)---ef(-1,r2)---ef(1,r1)---
    vvv(9)= - sin(theta)*g*i*g(l,v)

 
   Vertex 4: z(0,r2,v)---z(0,r3,va)---h0(0,r1)---
            v.va*g*i*wm
    vvv(4)=-------------
                      2
            cos(theta)

 
   Vertex 3: h0(0,r1)---w(1,r2,v)---w(-1,r3,va)---
    vvv(3)=v.va*g*i*wm
