Psigma - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Mozilla Firefox.

Online Help

All Products    Maple    MapleSim


Physics[Psigma] - the Pauli's 2 x 2 sigma matrices

Calling Sequence

Psigma[n]

Parameters

n

-

an integer between 0 and 4, or an algebraic expression representing it, identifying a Pauli matrix

Description

• 

The Psigma[n] command, where n ranges from 1 to 3, represents the three Pauli matrices, displayed on the screen as σn; these are the Hermitian and unitary matrices

σ1=0110,σ2=0−II0,σ3=100−1

  

where I is the imaginary unit (to represent it with a lowercase i, see interface,imaginaryunit). Psigma[...] can be also be indexed with the letters x, y and z, with the standard correspondence σx=σ1 σy=σ2 and σz=σ3, and also with `+` and `-` (include the ``), representing the ladder operators

σ+≡σ1+⁢I⁢σ2=0200

σ-≡σ1−I⁢σ2⁢=0020

  

To see the matrix form of any of these formulas use Physics:-Library:-RewriteInMatrixForm.

• 

The Pauli matrices satisfy the commutation relations

σa,σb−=2⁢I⁢εa,b,c⁢σc

σa,σb+=2⁢δa,b

  

where δa,b and εa,b,c are respectively the KroneckerDelta and LeviCivita symbols, and a,b,c range from 1 to 3. The Pauli matrices satisfy Det⁡σa=−1,Trace⁡σa=0, and σa2=1 (the 2 x 2 identity matrix), where Det represents the determinant, and Trace represents/computes the trace.

• 

When Physics is loaded, the three Pauli matrices Psigma[a], with a ranging from 1 to 3, together with Psigma[4] representing the 2 x 2 identity, are the components of a 4-vector in spacetime, Psigma[mu], with mu ranging from 1 to 4 and displayed as σμ. As with all spacetime tensors, you can use the value 0 of an index to refer to the position of the timelike component (by default equal to 4). The defining algebra for the 3D Pauli matrices σa shown above is extended to the 4D σμ as follows

σμ,σν−=2⁢I⁢ε4,μ,να4,μ,να⁢σα

σμ,σν+=2⁢σμ⁢δν4ν4+2⁢σν⁢δμ4μ4−2⁢gμ,ν

• 

Note: The default metric is of Minkowski type with signature (---+), so the space components of the contravariant σμ change in sign with respect to the definition of the Pauli matrices shown above.

• 

You can use Setup to change the metric or set a different signature, for example with the timelike component in position 1, as in (+---) or (-+++). That, however, makes the values 1, 2 and 3 of the index μ respectively refer to the identity matrix and the Pauli matrices 1, 2. To work with (+---) or (-+++) and avoid the inconvenience of having Psigma[1] referring to the identity matrix instead of σx, set spaceindices or su2indices, for example via Setup⁡su2indices=lowercaselatin_ah, and use Define with its redo option to redefine the type of tensor of Psigma, for example entering Define⁡redo,Psigmaa. In that way Psigma[a] becomes a 3D tensor, free of the issue mentioned about the numerical value of the index 1 in Psigma[mu].

Examples

> 

with⁡Physics:

> 

Setup⁡mathematicalnotation=true

mathematicalnotation=true

(1)
> 

Psigma1

σ1

(2)

You can see the matrix contents of σ1 in different ways, for example:

> 

Psigma1,matrix

σ1=0110

(3)

or for generic expressions involving tensors that represent matrices use

> 

Library:-RewriteInMatrixForm⁡

0110

(4)
> 

Psigma1⁢Psigma2+Psigma2⁢Psigma1

σ1⁢σ2+σ2⁢σ1

(5)
> 

Library:-RewriteInMatrixForm⁡

0110·0−II0+0−II0·0110

(6)

Besides Library:-RewriteInMatrixForm, which shows the matrix contents of a tensorial expression, to perform those matrix operations you can use

> 

Library:-PerformMatrixOperations⁡

0000

(7)

Among the basic properties of Pauli matrices, there are

> 

Psigma1⁢Psigma1

σ12

(8)
> 

Trace⁡

2

(9)
> 

Psigma1⁢Psigma2

σ1⁢σ2

(10)
> 

Trace⁡

0

(11)

To see the algebra satisfied by the Pauli matrices at any moment use Library:-DefaultAlgebraRules

> 

Library:-DefaultAlgebraRules⁡Psigma

σμ,σν−=2⁢I⁢ε4,μ,να4,μ,να⁢σα,σμ,σν+=2⁢σμ⁢δν4ν4+2⁢σν⁢δμ4μ4−2⁢gμ,ν

(12)

These equations can be verified in different ways. For example, construct an array with their components, then use its simplifier option to evaluate the commutators and anticommutators

> 

TensorArray⁡

σ1,σ1−=0σ1,σ2−=2⁢I⁢σ3σ1,σ3−=−2⁢I⁢σ2σ1,σ4−=0σ2,σ1−=−2⁢I⁢σ3σ2,σ2−=0σ2,σ3−=2⁢I⁢σ1σ2,σ4−=0σ3,σ1−=2⁢I⁢σ2σ3,σ2−=−2⁢I⁢σ1σ3,σ3−=0σ3,σ4−=0σ4,σ1−=0σ4,σ2−=0σ4,σ3−=0σ4,σ4−=0,σ1,σ1+=2σ1,σ2+=0σ1,σ3+=0σ1,σ4+=2⁢σ1σ2,σ1+=0σ2,σ2+=2σ2,σ3+=0σ2,σ4+=2⁢σ2σ3,σ1+=0σ3,σ2+=0σ3,σ3+=2σ3,σ4+=2⁢σ3σ4,σ1+=2⁢σ1σ4,σ2+=2⁢σ2σ4,σ3+=2⁢σ3σ4,σ4+=4⁢σ4−2

(13)

Note that in the lines above the matricial operations are performed abstractly, with the 2x2 matrices 0 and 1 (identity) omitted. To represent the algebra of the Pauli matrices with those two matrices not omitted, see the approach used in the MaplePrimes post Algebra of the Dirac matrices with an identity matrix on the right-hand side.

> 

TensorArray⁡,simplifier=value

0=02⁢I⁢ε1,2,4β1,2,4β⁢σβ=2⁢I⁢σ32⁢I⁢ε1,3,4β1,3,4β⁢σβ=−2⁢I⁢σ20=0−2⁢I⁢ε1,2,4β1,2,4β⁢σβ=−2⁢I⁢σ30=02⁢I⁢ε2,3,4β2,3,4β⁢σβ=2⁢I⁢σ10=0−2⁢I⁢ε1,3,4β1,3,4β⁢σβ=2⁢I⁢σ2−2⁢I⁢ε2,3,4β2,3,4β⁢σβ=−2⁢I⁢σ10=00=00=00=00=00=0,2=20=00=02⁢σ1=2⁢σ10=02=20=02⁢σ2=2⁢σ20=00=02=22⁢σ3=2⁢σ32⁢σ1=2⁢σ12⁢σ2=2⁢σ22⁢σ3=2⁢σ34⁢σ4−2=4⁢σ4−2

(14)

Alternatively, for instance, rewrite in matrix form the equations before computing the commutators, then activate the inert commutators using value

> 

Library:-RewriteInMatrixForm⁡1

0110,0110−=00000110,0−II0−=2⁢I⁢100−10110,100−1−=−2⁢I⁢0−II00110,1001−=00000−II0,0110−=−2⁢I⁢100−10−II0,0−II0−=00000−II0,100−1−=2⁢I⁢01100−II0,1001−=0000100−1,0110−=2⁢I⁢0−II0100−1,0−II0−=−2⁢I⁢0110100−1,100−1−=0000100−1,1001−=00001001,0110−=00001001,0−II0−=00001001,100−1−=00001001,1001−=0000

(15)
> 

map⁡expand,value⁡

0000=00002⁢I00−2⁢I=2⁢I00−2⁢I0−220=0−2200000=0000−2⁢I002⁢I=−2⁢I002⁢I0000=000002⁢I2⁢I0=02⁢I2⁢I00000=000002−20=02−200−2⁢I−2⁢I0=0−2⁢I−2⁢I00000=00000000=00000000=00000000=00000000=00000000=0000

(16)

A notational issue, correct but that could be seen as an inconvenience, happens when you set the signature with the timelike component in position 1, as in (+---) or (-+++), in that Psigma[1] points to Psigma[0], the identity 2x2 matrix instead of to Psigma[x]

> 

Setup⁡signature=`+---`

signature=+ - - -

(17)
> 

Psigma1=Psigma0

σ1=σ1

(18)
> 

Psigma1,matrix

σ1=1001

(19)

You can still refer to σx indexing with the letter x

> 

Psigmax=Library:-RewriteInMatrixForm⁡Psigmax

σ2=0110

(20)

To avoid this potential inconvenience you could set Psigma to be a 3D tensor. One way of doing that is to set spaceindices and redefine Psigma using Define with its redo option (necessary)

> 

Setup⁡spaceindices=lowercaselatin_is

spaceindices=lowercaselatin_is

(21)

Note that the redefinition requires passing Psigma indexed with a space index

> 

Define⁡redo,Psigmaj

Defined Pauli sigma matrices (Psigma): σ1 , σ2 , σ3

__________________________________________________

Defined objects with tensor properties

γμ,σj,∂μ,gμ,ν,γi,j,εα,β,μ,ν

(22)

Now we have σ1 = σx

> 

Psigma1,matrix

σ1=0110

(23)

Note that in this case the algebra is expressed in terms of the 3D metric, gamma_, displayed as γi,j.

> 

Library:-DefaultAlgebraRules⁡Psigma

σi,σj−=2⁢I⁢εi,j,k⁢σk,σi,σj+=2⁢γi,j

(24)
> 

TensorArray⁡

σ1,σ1−=0σ1,σ2−=2⁢I⁢σ3σ1,σ3−=−2⁢I⁢σ2σ2,σ1−=−2⁢I⁢σ3σ2,σ2−=0σ2,σ3−=2⁢I⁢σ1σ3,σ1−=2⁢I⁢σ2σ3,σ2−=−2⁢I⁢σ1σ3,σ3−=0,σ1,σ1+=2σ1,σ2+=0σ1,σ3+=0σ2,σ1+=0σ2,σ2+=2σ2,σ3+=0σ3,σ1+=0σ3,σ2+=0σ3,σ3+=2

(25)

Alternatively, to entirely detach the definition of the Pauli matrices from the details of the spacetime or space metric and signatures you can set Psigma as a tensor of a generic SU(2) space setting su2indices and redefining Psigma in the same way

> 

Setup⁡su2indices=lowercaselatin_ah

su2indices=lowercaselatin_ah

(26)

Note that the redefinition requires passing Psigma indexed with a su2 index

> 

Define⁡redo,Psigmaa

Defined objects with tensor properties

γμ,σa,∂μ,gμ,ν,γi,j,εα,β,μ,ν

(27)

Now, again, we have σ1 = σx

> 

Psigma1,matrix

σ1=0110

(28)

This time the algebra is expressed using δa,b, the KroneckerDelta, used to represent the metric in the SU(2) space, and the components of these tensorial equations are the same as those computed for Psigma as a 3D space tensor lines above

> 

Library:-DefaultAlgebraRules⁡Psigma

σa,σb−=2⁢I⁢εa,b,c⁢σc,σa,σb+=2⁢δa,b

(29)
> 

TensorArray⁡

σ1,σ1−=0σ1,σ2−=2⁢I⁢σ3σ1,σ3−=−2⁢I⁢σ2σ2,σ1−=−2⁢I⁢σ3σ2,σ2−=0σ2,σ3−=2⁢I⁢σ1σ3,σ1−=2⁢I⁢σ2σ3,σ2−=−2⁢I⁢σ1σ3,σ3−=0,σ1,σ1+=2σ1,σ2+=0σ1,σ3+=0σ2,σ1+=0σ2,σ2+=2σ2,σ3+=0σ3,σ1+=0σ3,σ2+=0σ3,σ3+=2

(30)

To activate the inert commutators and anticommutators use value

> 

value⁡

0=02⁢I⁢σ3=2⁢I⁢σ3−2⁢I⁢σ2=−2⁢I⁢σ2−2⁢I⁢σ3=−2⁢I⁢σ30=02⁢I⁢σ1=2⁢I⁢σ12⁢I⁢σ2=2⁢I⁢σ2−2⁢I⁢σ1=−2⁢I⁢σ10=0,2=20=00=00=02=20=00=00=02=2

(31)

See Also

Physics, Physics conventions, Physics examples, Physics Updates, Tensors - a complete guide, Mini-Course Computer Algebra for Physicists, Physics[*], Physics[Dgamma], Physics[Library], TensorArray, Trace, value