Message Boards

GROUPS:

Wolfram Science Physics Group Theory Wolfram Language Wolfram Summer School Wolfram Function Repository

Chang Wu

[WWS21] Exploring the gauge symmetries in Wolfram Model

Chang Wu, Department of Physics, Technion, Israel Institute of Technology

Posted 7 months ago

829 Views

0 Replies

2 Total Likes

Follow this post

Exploring the gauge symmetries in Wolfram Model
Chang Wu
Institut für Theoretische Physik, Georg-August-Universität Göttingen

According to the previous project “Full Discretization of Local Gauge Invariance”, the discrete gauge invariance techniques were previously developed [1]. Therefore, continue exploring the symmetry and related to some known physics in QFT and QM is interesting.

In order to study the gauge symmetries with the wolfram model, we need to reproduce a function like “QuantumToMultiwaySystem”, but instead of each state vertex being a pure state, each state vertex would now be a configuration of a field over one-dimensional lattice, with the evolution rules being some discrete approximation to the equation of motion, such as Klein-Gordon equation for the scalar field.

As a first step, we use the “QuantumToMultiwaySystem” [2] to study a few well known operators, such as the Pauli and Gell-Mann matrices. And apply it to calculate the colour factor in QCD.

A First Example

As a first example we study a simple operator which create a superposition of basis state. It has the matrix:

H

1	1
1	-1

,

This operator is called as Hadamard gate in quantum computing [3], apply it to one of the basis gives

H0〉0〉+1〉.H1〉0〉-1〉.

Therefore, we could apply the “QuantumToMultiwaySystem”, which gives the evolution graph as

In[]:=

ResourceFunction["QuantumToMultiwaySystem"][<|"Operator"{{1,1},{1,-1}},"Basis"{{1,0},{0,1}}|>,{1,1},2,"EvolutionGraphFull","IncludeStatePathWeights"True,VertexLabels"VertexWeight",AspectRatio0.5]

Out[]=

In[]:=

ResourceFunction["QuantumToMultiwaySystem"][<|"Operator"{{1,1},{1,-1}},"Basis"{{1,0},{0,1}}|>,{1,1},2,"EvolutionCausalGraph"]

Out[]=

This simple operator can be decomposed into a linear combination of Pauli matrices

H

where the Pauli matrices are usually taken in the form



0	1
1	0



0	-i
i	0



1	0
0	-1

“ 4-dimensional” Pauli Matrices

It is also possible to construct the Hilbert space for systems containing 2 spin systems, for example spinors in relativistic quantum mechanics. Using the tensor product of each vector, we have

{0〉,1〉}×{0〉,1〉}{00〉,01〉,10〉,11〉},

which the matrix representation are

Similarly, we could construct the “ 4-dimensional” Pauli Matrices



for simplicity, we only collect the result below:



0	1	1	0
1	0	0	1
1	0	0	1
0	1	1	0



0	- i	- i	0
i	0	0	- i
i	0	0	- i
0	i	i	0



2	0	0	0
0	0	0	0
0	0	0	0
0	0	0	-2

which is easy to varify the commutation relations, for example:

In[]:=

sigmax={{0,1,1,0},{1,0,0,1},{1,0,0,1},{0,1,1,0}};sigmay={{0,-I,-I,0},{I,0,0,-I},{I,0,0,-I},{0,I,I,0}};sigmax.sigmay-sigmay.sigmax//MatrixForm

Out[]//MatrixForm=

4	0	0	0
0	0	0	0
0	0	0	0
0	0	0	-4

Thus, we can apply the “QuantumToMultiwaySystem” to study

act on the basis states

In[]:=

ResourceFunction["QuantumToMultiwaySystem"][<|"Operator"sigmax,"Basis"{{1,0,0,0},{0,1,0,0},{0,0,1,0},{0,0,0,1}}|>,{1,0,0,0},3,"EvolutionGraphFull","IncludeStatePathWeights"True,VertexLabels"VertexWeight",AspectRatio0.8]

Out[]=

In[]:=

ResourceFunction["QuantumToMultiwaySystem"][<|"Operator"sigmax,"Basis"{{1,0,0,0},{0,1,0,0},{0,0,1,0},{0,0,0,1}}|>,{1,0,0,0},4,"EvolutionCausalGraph"]

Out[]=

The “total spin” S can be identified from the eigenvalue of the operator, where in Multiway system, the eigenvalue can be obtained from the vertex weight



0	1	1	0
1	0	0	1
1	0	0	1
0	1	1	0

0	- i	- i	0
i	0	0	- i
i	0	0	- i
0	i	i	0

2	0	0	0
0	0	0	0
0	0	0	0
0	0	0	-2

In[]:=

ResourceFunction["QuantumToMultiwaySystem"][<|"Operator"{{8,0,0,0},{0,4,4,0},{0,4,4,0},{0,0,0,8}},"Basis"{{1,0,0,0},{0,1,0,0},{0,0,1,0},{0,0,0,1}}|>,{1,0,0,0},1,"EvolutionGraphFull","IncludeStatePathWeights"True,VertexLabels"VertexWeight"]

Out[]=

The

4×4

matrix operator

can also formed in other 4-dim representation, such as



σ	0
0	σ

In[]:=

sigma2x={{0,1,0,0},{1,0,0,0},{0,0,0,1},{0,0,1,0}};sigma2y={{0,-I,0,0},{I,0,0,0},{0,0,0,-I},{0,0,I,0}};sigma2z={{1,0,0,0},{0,-1,0,0},{0,0,1,0},{0,0,0,-1}};sigma2x.sigma2y-sigma2y.sigma2x//MatrixForm

Out[]//MatrixForm=

2	0	0	0
0	-2	0	0
0	0	2	0
0	0	0	-2

In[]:=

ResourceFunction["QuantumToMultiwaySystem"][<|"Operator"sigma2x,"Basis"{{1,0,0,0},{0,1,0,0},{0,0,1,0},{0,0,0,1}}|>,{1,1,1,1},4,"EvolutionCausalGraph"]

Out[]=

Moreover, the Dirac–Pauli representation is useful to express the Dirac equation in covariant form

i

∂

-mψ0,

where the

matrices are



I	0
0	-I



0	σ
-σ	0

In[]:=

sigma3={{1,0,0,0},{0,1,0,0},{0,0,-1,0},{0,0,0,-1}};sigma3x={{0,0,0,1},{0,0,1,0},{0,-1,0,0},{-1,0,0,0}};sigma3y={{0,0,0,-I},{0,0,I,0},{0,I,0,0},{-I,0,0,0}};sigma3z={{0,0,1,0},{0,0,0,-1},{-1,0,0,0},{0,1,0,0}};sigma3x.sigma3z-sigma3z.sigma3x//MatrixForm

In[]:=

ResourceFunction["QuantumToMultiwaySystem"][<|"Operator"sigma3x,"Basis"{{1,0,0,0},{0,1,0,0},{0,0,1,0},{0,0,0,1}}|>,{1,1,1,1},4,"EvolutionCausalGraph"]

Out[]=

Gell-Mann Matrices

Start from the QCD Lagrangian, where the color is described by the local gauge transformations with respect to color SU(3) group

ψ(x)

(x)

ψ(x),

and similar to QED case, we could introduce the covariant derivative



∂

-ig

(x)

where

(x)

is the gauge (gluon) fields, transfer as

(x)

(x)-

∂

(x)-

abc

(x)

One the other hand, one can also start by describe a quark as Dirac spinor

u(p)

plus a colour D.O.F. described by a column colour vector:

r

,b

,r

And mathematically, quark colour transforms under the fundamental representation, while gluon responsible for exchanging momentum and colour between quarks, are transfer under the adjoint representation of SU(3), which are described by the Gell-Mann matrices

In[]:=

lambda1={{0,1,0},{1,0,0},{0,0,0}};(*lambda1//MatrixForm*)lambda2={{0,-I,0},{I,0,0},{0,0,0}};(*lambda2//MatrixForm*)lambda3={{1,0,0},{0,-1,0},{0,0,0}};(*lambda3//MatrixForm*)lambda4={{0,0,1},{0,0,0},{1,0,0}};(*lambda4//MatrixForm*)lambda5={{0,0,-I},{0,0,0},{I,0,0}};(*lambda5//MatrixForm*)lambda6={{0,0,0},{0,0,1},{0,1,0}};(*lambda6//MatrixForm*)lambda7={{0,0,0},{0,0,-I},{0,I,0}};(*lambda7//MatrixForm*)lambda8={{1,0,0},{0,1,0},{0,0,-2}}/Sqrt[3];(*lambda8//MatrixForm*)TabL={lambda1,lambda2,lambda3,lambda4,lambda5,lambda6,lambda7,Sqrt[3]lambda8};

There are 8 independent color states, and one colour singlet state, as

3×38+1

for SU(3). For example:



0	1	0
1	0	0
0	0	0

r

The Gell-Mann matrices are one possible representation of the infinitesimal generators of the special unitary group called SU(3), and we can use the Multiway system to study it’s eigenstate. For example

In[]:=

ResourceFunction["QuantumToMultiwaySystem"][<|"Operator"TabL[[8]],"Basis"{{1,0,0},{0,1,0},{0,0,1}}|>,{1,1,1},2,"EvolutionGraphFull","IncludeStatePathWeights"True,VertexLabels"VertexWeight",EdgeLabels"EdgeWeight",AspectRatio0.5]

Out[]=

Colour Factor

The Gell-Mann matrices is useful in scattering theory. For example: at tree level, the amplitude can be decomposed as [4]

M(ng)

∑

perm(2,⋯,n)

Tr(

⋯λ

)A(1⋯n)

Here, we use the

q

scattering process to study this colour decomposition. To write down the matrix element, we need to using the Feynman rule and follow the fermions arrow backwards

M

[

][

]

This matrix element is similar to electromagnetic scattering except the strong coupling constant, i.e.

e

, and the addition term

called as colour factor, where



†

and

is the colour vector, i.e. r, g, b. Moreover, we could define the colour factor as

f

. Then, in the lowest order approximation, the dynamics of q\bar{q} is the same as

scattering, where the Coloumb-like potential is

-

In order to calculate the colour factor, we choose colours for

and

. And the theory is invariant under rotations in colour space, thus any choice of colours will give the same answer. For this process, we have three colour options, i.e.

r

and

b

. Thus we have



∑



∑



∑

Now, in order to obtain the matrices with non-zero element, we can apply the Multiway system to the Gell-Mann matrices:

In[]:=

Table[ResourceFunction["QuantumToMultiwaySystem"][<|"Operator"TabL[[n]],"Basis"{{1,0,0},{0,1,0},{0,0,1}}|>,{1,1,1},1,"EvolutionGraphFull",VertexLabels"VertexWeight",AspectRatio1],{n,1,8}]

Out[]=





In[]:=

Table[ResourceFunction["QuantumToMultiwaySystem"][<|"Operator"TabL[[n]],"Basis"{{1,0,0},{0,1,0},{0,0,1}}|>,{1,1,1},1,"EvolutionCausalGraph",VertexLabels"VertexWeight",AspectRatio1],{n,1,8}]

Out[]=





In[]:=

Table[VertexList[ResourceFunction["QuantumToMultiwaySystem"][<|"Operator"TabL[[n]],"Basis"{{1,0,0},{0,1,0},{0,0,1}}|>,{1,1,0},1,"EvolutionGraphFull"]][[2]][[2]],{n,1,8}]

Out[]=

{{0, 1, 0},{0, I, 0},{0, -1, 0},{0, 0, 1},{0, 0, I},{0, 0, 1},{0, 0, I},{0, 1, 0}}

Where the only non-zero result for

and

, i.e.

and

are

and

. And we can get the eigenvalue from the vertex weight. Thus we have



(1)(1)+



and



(1)(-1)+

-

Similarly, the only matrices with non-zero element for

are

and

.Therefore,theresultreadas



[(-i)(i)+(1)(1)]

Again, we can compare to the

potential, where the minus sign indicate it is repulsive. Furthermore, we could also apply to the colour singlet case for both out-going and in-coming q and

r



The colour factor can be obtained as

f

∑

Tr(

)

= 4/3, which is known as SU(3) quadratic Casimir.

Next Step

In this post, we have applied the Multiway system to study the QCD colour factor, where the higher-dimensional colour representation [5] is important to study the sub-leading colour effects. Therefore, we could continue developing this package and applying it to some processes with more complicated colour structure.

For a more interesting direction, we could develop a new function similar to the “QuantumToMultiwaySystem”, but each state vertex is replaced by a real scalar field over one-dimensional lattice structure, for example, the 1-dim Ising model. The Ising Model is a model of spins on a lattice with nearest-neighbor interactions. [6]

From the point of view in Wilsonian effective field theory, we can zoom out by ‘integrating out’ half of the spins on the lattice, which leaves a new effective theory for the remainder. Thus, at long distance the model is described by the QFT with action

S

(∂ϕ)

λϕ

In Multiway system, the Ising model have been studied with the string substitution systems [7]. For example:

In[]:=

ResourceFunction["MultiwaySystem"][{"A""AB","B""BA"},"A",5,"BranchialGraph"]

Out[]=

Moreover, the ultimate goal would be looking an alternate way for spontaneous symmetry-breaking that arise from the discreteness of the lattice over which the scalar field is defined, where the EOM is the discretized Klein-Gordon equation with wavefunctions replaced by operators acting on states, analogous to the generation of Goldstone bosons.

Keywords

◼

Multiway System

◼

Gauge Symmetries

◼

Pauli and Gell-Mann Matrices

◼

SU(3) quadratic Casimir

Acknowledgment

Mentor: Jon Lederman

The author would like to thank Jon Lederman, Jonathan Gorard, Xerxes D. Arsiwalla, and Stephan Wolfram for the useful discussions, and especially thank Xerxes D. Arsiwalla for pointed a possible new direction for this project.

References

◼

[1] https://community.wolfram.com/groups/-/m/t/2030337

◼

[2] https://resources.wolframcloud.com/FunctionRepository/resources/QuantumToMultiwaySystem

◼

[3] https://en.wikipedia.org/wiki/Quantum_logic_gate#Hadamard_(H)_gate

◼

[4] M. L. Mangano, S. J. Parke, and Z. Xu, “Duality and Multi - Gluon Scattering,” Nucl. Phys. B, vol. 298, pp. 653–672, 1988.

◼

[5] Y. Bao and A. J. Larkoski, Calculating Pull for Non-Singlet Jets, arXiv:1910.02085.

◼

[6] Jared Kaplan - QFT Lectures Notes

◼

[7] https://community.wolfram.com/groups/-/m/t/2029608

POSTED BY: Chang Wu

Answer

Reply to this discussion

Reply Preview

Attachments

Remove Add a file to this post

Follow this discussion

or Discard

Group Abstract

Feedback