Exotic Particle Dynamics Using Novel Hermitian Spin Matrices

: In this work, an analogue to the Pauli spin matrices is presented and investigated. The proposed Hermitian spin matrices exhibit four symmetries for spin-1/ n particles. The spin projection operators are derived, and the electrodynamics for hypothetical spin-1/2 fermions are explored using the proposed spin matrices. The fermionic quantum Heisenberg model is constructed using the proposed spin matrices, and comparative studies against simulation results using the Pauli spin matrices are conducted. Further analysis of the key ﬁndings as well as discussions on extending the proposed spin matrix framework to describe hypothetical bosonic systems (spin-1 particles) are provided.


Introduction
In this work, four novel Hermitian spin matrices are introduced as analogues to the Pauli spin matrices.The proposed spin matrices exhibit the following symmetrical relations with their respective basis using the novel spin matrices for particles with spin 1/n, where n > 0: Relation II : Relation III : Relation IV : where σ 2 = −σ 3 , and the complex conjugate relation σ 0 = σ 1 .Four possible symmetrical relations with their respective basis are then ascertained using the proposed spin matrices for n > 0. Although the analysis in this work could be carried out with any of the symmetrical relations given in Equation ( 1), the scope of the present manuscript is limited to Relation I. Thus, this work aims to explore hypothetical fermion dynamics using the proposed Hermitian spin matrices for Relation I in Equation ( 1) at n = 2.
explored the topology of a pseudo-Hermitian Chern insulator defined using the basis of q-deformed Pauli matrices (related to deformed algebras).The key finding of that work was obtained for a completely nonequilibrium case where the quantum evolution after quenching was dictated by a Floquet pseudo-Hermitian Hamiltonian.Similarly, in Fring and Taira (2020) [15] the authors employed a pseudo-Hermitian approach to Goldstone's theorem in non-Abelian and non-Hermitian quantum field theories.In that work, a detailed analysis of a non-Hermitian field theory with two complex scalar fields (with two components) exhibiting SU(2) symmetry was presented.Zhu et al. (2021) [16] demonstrated that several two-and three-dimensional pseudo-Hermitian phases could be constructed using q-deformed matrices.In addition to investigations related to topological bulk states, quantum metrics and non-Abelian tensor Berry connections, an experimental protocol was proposed for the empirical validation (of the proposed models).The following works provide current developments on research efforts in the direction of non-Hermitian and pseudo-Hermitian frameworks: Okuma and Sato, (2023) [17], Ashida et al., (2020) [5], Kunst and Dwivedi (2019) [18] and Feinberg and Riser, (2021) [19].This paper is organized as follows: The second and third sections describe the properties of the proposed spin matrices and their respective projection operators.In the fourth section, the electrodynamics of hypothetical fermions are explored using the proposed spin matrices as operators in the Schrödinger-Pauli and Dirac equations.The fifth section compares spin chain simulation results using the proposed spin matrices and the Pauli spin matrices.The sixth section provides further analysis of bosonic systems.This paper ends with the key conclusions from this research along with some ideas on directions for future research.

Analogue Spin Matrices
Pauli spin matrices serve as quantum operators corresponding to observables for the spin of fermions in each spatial direction: In this work, four spin matrices are proposed as analogues to the Pauli spin matrices, σ i (see Equation ( 1)).The proposed spin matrices, σ i , are constructed from the Pauli matrices as follows: For a particle with spin 1 n , the proposed Hermitian spin matrices follow the complex conjugate transpose relations: The proposed matrices in Relation I exhibit the following symmetry: In addition, the identity matrix is represented as follows: The parameter, n > 0, is real-valued.If the parameter, n, is complex-valued, then the proposed matrices will lose their Hermitian properties.Similar to the Pauli matrices, the following matrix properties apply to arbitrary n: However, unlike the Pauli matrices, the proposed matrices are not involutory.The following relationships exist between the proposed matrices: σ 2 = −σ 3 and σ 0 = σ 1 .This introduces the possibility of simplifying the analysis of three-dimensional observables to two dimensions.The proposed matrices have the following commutative properties: The commutation relations for the proposed matrices, where [a, b] = ab − ba, are as follows: On the other hand, the anti-commutation relations for the proposed matrices, where {a, b} = ab + ba, are as follows: The anti-commutator operation for the squared spin matrices, (σ 0 ) 2 and (σ 1 ) 2 , is as follows: The vector properties of the proposed spin matrices are proved in the Appendix A. Using the proposed matrices, an analogue to the gamma matrices in the Dirac basis is then constructed for Relation I: The proposed time-like gamma matrix is then as follows: As with the Pauli matrices, all the proposed gamma matrices are Hermitian.Hence, these matrices yield real eigenvalues (real quantum energy states).The anticommutation relations for the proposed gamma matrices are as follows: As with the conventional gamma matrices, the matrix properties for n = 2 are as follows: det

Spin Projection Operators
In this section, the spin projection operators are obtained for the exotic fermions using the proposed spin matrices.A formulation for spin projection operators for fermions with spin-1 n at n = 2 is provided.The two-component spinor is then employed to represent the quantum state of a fermion using the proposed spin matrices.The spin projection operator, S 1 , is expressed as follows: where is the reduced Planck's constant.The eigenvalues for S 1 are The eigenspinor representations are then as follows: The spin projection operator, S 2 , is expressed as follows: The eigenvectors for S 2 are v 1 = [0, 1] and v 2 = [1, 0], and the eigenspinor represen- tations are then as follows: The spin projection operator, S 3 , is expressed as follows: The eigenvectors for S 3 are v 1 = [1, 0] and v 2 = [0, 1], and the eigenspinor represen- tations are then as follows: An interesting difference with the Pauli spin projection operators is that for the operator S 1 , the |↓ spin state contains a small imaginary contribution of 3.92523 × 10 −17 i.

Electrodynamics
In this section, the non-relativistic electrodynamics of exotic fermions are explored using the Schrödinger-Pauli equation [20,21].The Schrödinger-Pauli equation is presented using the proposed spin matrices at n = 2: where m is the fermion mass, p is the vector form of the momentum operator, → A is the magnetic vector potential, q is the fermion's electric change, φ is the electric scalar potential, Ĥ0 is the Hamiltonian (using the proposed spin matrices) and |ψ is the quantum state.The spin matrices are represented as the following vector: σ = σ i , σ j , σ k .By implementing the proposed spin matrices, the Schrödinger-Pauli equation yields solutions for the Hamiltonian, Ĥ. Considering that the fermion (at n = 2) is subjected to a constant magnetic field within the Landau gauge [22], where B is the uniform magnetic field.Solving the Schrödinger-Pauli equation for Ĥ0 gives the following Hamiltonian operators: This analysis is extended using the Landau symmetric gauge, where the magnetic vector potential is given as follows: Solving the Schrödinger-Pauli equation within the Landau symmetric gauge using the spin matrix, σ , gives the following Hamiltonian operators: In the context of relativistic quantum electrodynamics, the following Dirac equation is considered to determine the rest energy of the fermion placed in the electric potential, where m is the mass of the fermion, c is the speed of light and γ j is the proposed gamma matrices.As with the conventional gamma matrices in the Dirac basis, it can be shown using Equation ( 25) that the rest energy of the hypothetical fermion at → A = 0 is equivalent to the energy of a particle placed in the electric potential, cqA 0 .

Simulation with Spin Chains
The proposed spin matrices are compared with the Pauli spin matrices by analyzing their respective behaviors in the context of a spin chain using the spin − 1 2 quantum Heisenberg model [23,24].The one-dimensional Heisenberg model is utilized where magnetic interactions take place specifically between adjacent dipoles: where J is the coupling constant, h is the external magnetic field, and the dipoles are described as the quantum operator acting on the Kronecker product of dimensions, 2 N .Considering the coupling constant, J = J x , J y , J z , to be real-valued, the Hamiltonian operator is represented as follows: where σ k j is the k th Pauli matrix on the j th lattice point with periodic boundary conditions.The Pauli matrices, σ k j , is defined as . In this work, the Heisenberg XXX model with J > 0 is employed, where the coupling constants conform to the following equivalence: J = J x = J y = J z .The simulation is carried out using the Python programming language.To simplify the simulation, a qubit system is considered where the number of spins in the chain is limited to two.The following parameters are fixed in the simulation: coupling constant, J = 1 × 10 −19 and B = 1 × 10 −19 T. The simulation is performed at three temperature values: T 1 = 10 K (low), T 2 = 1000 K (medium), T 4 = 10,000 K (high) and T 4 = 100,000 K (very high).The spectrum of the Hamiltonian (i.e., energy states, λ(σ) = H(σ)), is then obtained for each spin state: σ = {(1, 1), (1, −1), (−1, 1), (1, 1)}.The configuration probability for each spin state, P(σ), is determined as follows: where Z is the partition function and k B is the Boltzmann constant: Computational experiments were performed by executing the simulation using Pauli matrices and the proposed spin matrices.The simulation routine for the spin chain used in this work is given in Algorithm 1: Randomly generate initial spin states.
Define Zeeman term to account for an external magnetic field.5.
Compute the energy of the system and determine quantum states (eigenvalues).6.
Compute partition function and determine density matrix (state probabilities).7.
Apply Bures metric to determine the difference in the generated density states.8.
Output: Print and plot state probabilities and energy levels.

END
Figure 1 provides the state probabilities generated from the simulations using the proposed spin matrices and Pauli matrices at different temperatures.In Figure 1, the comparison of the state probabilities obtained using the proposed spin matrices and those generated using the proposed spin matrices for temperatures T = 100,000 K, 10,000 K, 1000 K and 10 K is shown.It can be observed that the overall dynamics of the system are similar to that of spin chains using Pauli matrices where the particles experience losses in magnetic spin orientation at higher temperatures, while the restoration of orientation rises at lower temperatures.To measure the difference between the state probabilities of the proposed matrices and the Pauli matrices, the Bures metric (or the Bures distance) is employed [25].The Bures metric measures the statistical distance between two quantum states represented by their respective density matrices,  .The Bures distance is measured as follows: Tr    (30) In Figure 1, the comparison of the state probabilities obtained using the proposed spin matrices and those generated using the proposed spin matrices for temperatures T = 100,000 K, 10,000 K, 1000 K and 10 K is shown.It can be observed that the overall dynamics of the system are similar to that of spin chains using Pauli matrices where the particles experience losses in magnetic spin orientation at higher temperatures, while the restoration of orientation rises at lower temperatures.To measure the difference between the state probabilities of the proposed matrices and the Pauli matrices, the Bures metric (or the Bures distance) is employed [25].The Bures metric measures the statistical distance between two quantum states represented by their respective density matrices, ρ.The Bures distance is measured as follows: Figure 2 shows the Bures distance for the density states produced using the proposed matrices versus the Pauli matrices in the spin chain simulations using the Heisenberg spin model.It can be seen in Figure 1 that at very high temperatures (T = 100,000 K), simulations undertaken using both sets of matrices generate similar state probabilities.At T = 100,000 K, the results show a somewhat equal distribution of probabilities among the four states.Despite this general similarity, the exact values of the state probabilities generated by the simulations using the Pauli matrices are distinct from those obtained using the proposed spin matrices (see Figure 1).However, as the temperature is lowered, the simulations conducted using the proposed spin matrices favor the state (1, −1).On the other hand, the simulations carried out using the Pauli matrices lean toward equal probabilities between states (1, 1) and (−1, 1) as the temperature is lowered (T → 0 K).This behavior is captured accurately using the Bures metric.The state probabilities (density matrices) generated by the simulations using the proposed spin matrices differ significantly from those obtained using the Pauli spin matrices as the temperature is lowered, where the Bures distance (D) reaches its maximum value.

Analysis
Another key distinctive feature of the proposed spin matrices compared to the Pauli matrices is that they introduce a different type of algebra in relation to their commuting and anti-commuting properties (see Section 2).As seen in Section 3, the proposed spin matrices also produce spin states that contain a constant small imaginary contribution of 3.92523 × 10  .The electrodynamics exploration conducted in Section 4 using the Schrödinger-Pauli equation and the proposed spin matrices yielded different Hamiltonian expressions compared to that using Pauli matrices.The computation of the rest energy of a theoretical fermion using the proposed gamma matrices is consistent with the analysis of the Dirac equation performed using conventional gamma matrices.As per this line of reasoning, it is also possible to construct higher spin systems (e.g., for bosons) using the proposed spin matrices.For instance, a set of Hermitian spin-1 matrices for triplet states using the proposed spin matrices at  = 1 is constructed.The spin projection operator,  , is expressed as follows: It can be seen in Figure 1 that at very high temperatures (T = 100,000 K), simulations undertaken using both sets of matrices generate similar state probabilities.At T = 100,000 K, the results show a somewhat equal distribution of probabilities among the four states.Despite this general similarity, the exact values of the state probabilities generated by the simulations using the Pauli matrices are distinct from those obtained using the proposed spin matrices (see Figure 1).However, as the temperature is lowered, the simulations conducted using the proposed spin matrices favor the state (1, −1).On the other hand, the simulations carried out using the Pauli matrices lean toward equal probabilities between states (1, 1) and (−1, 1) as the temperature is lowered ( T → 0 K ).This behavior is captured accurately using the Bures metric.The state probabilities (density matrices) generated by the simulations using the proposed spin matrices differ significantly from those obtained using the Pauli spin matrices as the temperature is lowered, where the Bures distance (D) reaches its maximum value.

Analysis
Another key distinctive feature of the proposed spin matrices compared to the Pauli matrices is that they introduce a different type of algebra in relation to their commuting and anti-commuting properties (see Section 2).As seen in Section 3, the proposed spin matrices also produce spin states that contain a constant small imaginary contribution of 3.92523 × 10 −17 i.The electrodynamics exploration conducted in Section 4 using the Schrödinger-Pauli equation and the proposed spin matrices yielded different Hamiltonian expressions compared to that using Pauli matrices.The computation of the rest energy of a theoretical fermion using the proposed gamma matrices is consistent with the analysis of the Dirac equation performed using conventional gamma matrices.As per this line of reasoning, it is also possible to construct higher spin systems (e.g., for bosons) using the proposed spin matrices.For instance, a set of Hermitian spin-1 matrices for triplet states using the proposed spin matrices at n = 1 is constructed.The spin projection operator, S 1 , is expressed as follows: The eigenvalues for S 1 are . The eigenspinor representations are then as follows: The spin projection operator, S 2 , is expressed as follows: The eigenvalues for S 2 are λ 1 = − , λ 2 = and λ 3 = 0.The eigenvectors for S 1 are v 1 = [0, 0, 1], v 2 = [1, 0, 0] and v 3 = [0, 1, 0].The eigenspinor representations are then as follows: The spin projection operator, S 3 , is expressed as follows: The eigenvalues for S 3 are λ 1 = − , λ 2 = and λ 3 = 0.The eigenvectors for S 1 are v 1 = [1, 0, 0], v 2 = [0, 0, 1] and v 3 = [0, 1, 0].The eigenspinor representations are then as follows: The proposed spin matrices allow theoretical explorations of particles with arbitrary spins by defining the appropriate n.In Figure 1, it can be observed that the spin states obtained by the simulations of the quantum Heisenberg model using the proposed spin matrices differ from those generated using the Pauli spin matrices.Since the system of proposed spin matrices is completely Hermitian, real eigenstates (energy states) are consistently obtained.As observed in Section 2, the relationships σ 0 = σ 1 and σ 2 = −σ 3 could be established.This introduces the possibility of reducing the conventional threedimensional observable on spin systems to a two-dimensional one.Considering this reduction in the dimension, the proposed spin matrices may have applications in describing two-dimensional quasiparticles-e.g., anyons and the fractional quantum Hall effect [26,27].It is also conjectured that the proposed spin matrices and their underlying symmetry may have applications in the particle physics of more exotic forms of matter-e.g., dark matter.The analysis presented in this work provides another perspective on the magnetic properties (spin) of quantum particles.The proposed gamma matrices generated within the Majorana or Weyl basis may have applications in the design of superconductors [28,29].

Conclusions & Future Work
In this work, a novel basis of four Hermitian spin matrices constructed using Pauli matrices for arbitrary spin, 1/n, is proposed.The primary strength of the proposed spin matrices is their Hermitian property (similar to the Pauli matrices).This property plays a central role in generating physical observables by providing real-valued solutions to the energy eigenstates of quantum systems.This was demonstrated in this work by applying the proposed spin matrices to the Pauli equation and the spin chain simulation (using the Heisenberg model).Another strength of the proposed spin matrices is that it could be used to generate gamma matrices for applications in relativistic quantum mechanics.Additionally, the symmetry property of the proposed spin matrices could be generalized to particles with arbitrary spins.Another key advantage of the proposed spin matrices is their potential to describe two-dimensional observables in three-dimensional spatial settings.In addition to providing a means to perform more simplified analyses on quantum systems, they may also be a gateway towards understanding the behaviors of quasiparticles (anyons).However, a weakness of the proposed spin matrices is that they are not involutory (unlike the Pauli matrices).
Future work could be directed towards exploring other spin state systems (where n = 2) as provided in this section.An in-depth investigation of the key characteristics of all the relations in Equation ( 1) could be performed.Generalizations of the proposed matrices to produce analogues to the Gell-Mann matrices to theoretically investigate particle physics involving strong interactions could be conducted.In addition, the proposed gamma matrices and their formulations in the Weyl and Majorana basis could also be carried out.Finally, the implications of the proposed gamma matrices on quantum interactions via field theory would be an interesting avenue for potential research.

Figure 1 .
Figure 1.State probabilities from simulations using the proposed spin matrices and Pauli matrices at different temperatures.

Figure 1 .
Figure 1.State probabilities from simulations using the proposed spin matrices and Pauli matrices at different temperatures.

Figure 2 .
Figure 2. Bures distance for the density states generated using the proposed matrices versus the Pauli matrices.

Figure 2 .
Figure 2. Bures distance for the density states generated using the proposed matrices versus the Pauli matrices.