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
Psigma1
σ1
You can see the matrix contents of σ1 in different ways, for example:
Psigma1,matrix
σ1=0110
or for generic expressions involving tensors that represent matrices use
Library:-RewriteInMatrixForm⁡
0110
Psigma1⁢Psigma2+Psigma2⁢Psigma1
σ1⁢σ2+σ2⁢σ1
0110·0−II0+0−II0·0110
Besides Library:-RewriteInMatrixForm, which shows the matrix contents of a tensorial expression, to perform those matrix operations you can use
Library:-PerformMatrixOperations⁡
0000
Among the basic properties of Pauli matrices, there are
Psigma1⁢Psigma1
σ12
Trace⁡
2
Psigma1⁢Psigma2
σ1⁢σ2
0
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μ,ν
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
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
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
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
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=+ - - -
Psigma1=Psigma0
σ1=σ1
σ1=1001
You can still refer to σx indexing with the letter x
Psigmax=Library:-RewriteInMatrixForm⁡Psigmax
σ2=0110
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
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,εα,β,μ,ν
Now we have σ1 = σx
Note that in this case the algebra is expressed in terms of the 3D metric, gamma_, displayed as γi,j.
σi,σj−=2⁢I⁢εi,j,k⁢σk,σi,σj+=2⁢γi,j
σ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
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
Note that the redefinition requires passing Psigma indexed with a su2 index
Define⁡redo,Psigmaa
γμ,σa,∂μ,gμ,ν,γi,j,εα,β,μ,ν
Now, again, we have σ1 = σx
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
σa,σb−=2⁢I⁢εa,b,c⁢σc,σa,σb+=2⁢δa,b
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
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
Compatibility
The Physics[Psigma] command was updated in Maple 2020.
Download Help Document