Energy-States of Particles with Representational Spin

In this paper, an algorithm to produce a transition matrix between all states in ideal anti-paramagetic system under the simulated annealing condition that only moves to lower energy states are accepted. The check the accuracy of the transition matrix is confirmed by computer simulation, showing close agreement between model and simulation results.


Introduction
In this paper we compute all transition probabilities of a system of N isolated particles of spin M to move to any lower energy state.In each move, the current state ψ , is compared with a random trial state, t ψ .The move is accepted if it the trial state has a lower energy, and rejected if the energy is equal, or higher.No restrictions are placed upon the trial state: it is totally random, thus the states form a Markov chain [1] which the next state chosen at random, via the Monte Carlo method [2].This is a kind of simulated annealing [3], but with a simulation temperature of zero [4], which continues to be relevant within the context of recent work [5].The such as spin glasses, simulated annealing is used to find the optimal ground state [6].By contrast, in this paper, the system is rather simple and the optimal state is known.Thus we can compute the probability to move from any state to any other state which should lend insight into how this very powerful tool (simulated annealing) operates.We can pause the simulation at any point, and since we know the optimal state, measure how far the system is from optimal, The number of moves to find the best state increases rapidly with the complexity of the system.In the magnetic systems we study in this paper, this complexity can be increased by increasing the number of spins in the system, or by increasing the spin of each particle.Either way greatly increases the number of states for the system to explore.

Plan of the Paper
Beginning with two, three and four spin one half particle states, the transition probabilities between all states are enumerated.These lists of probabilities are formed into matrices.This is repeated for spin one particles in section 2 where the transition probabilities are derived from the degeneracy of the states.These degeneracies employ the generalized binomial coefficients, called the multinomials.In Section 3, the spin 3/2 particles are considered.Then, in Section 4, the entries in the transition matrices are derived from the degeneracy of the states, using a simple summation rule.In Section 5, degeneracy vectors are given for various systems.In Section 6, the results of numerical simulations are provided which confirm the results of the matrix calculations.

Assumptions of Model
Our underlying assumptions are: • 0 T → , thus eliminating the effects of entropy 0 S ∴ → • ψ states gravitate towards more preferred states of lower magnetism number due to our model being an anti-paramagnet • Particles in ψ are indistinguishable.
In any system N-particles with spin we will assign values.Those values will be relative to the system we are working with.We will assign +1 for spin up, and values −1 for spin down in a spin one half system.+1, 0, −1 for a spin one system, and 3 ± , 1 ± for a spin 3/2 system.The magnetism number will be achieved by virtue of the magnetism operator S.
Definition 1.1.For any state ψ in a N-particle system we define S to be the operator such that S ψ θ γ = − where θ is the total spin up, γ Total spin down.This leads us to the following proposition.Proposition 1.6.There can only be one idealized Ground state; that is 1 The only way to get magnetism number N − for an N-particle system, is to have all the particles in that system be spin down.Since particles are indistinguishable by assumption 3, there is only one state that gives magnetism number N − , thus 1 is the set of all states that are preferred to a state with magnetism number .i

Markov Chains
Definition 1.8.A system with a finite or countably infinite number of states with the property that given a present state, past states have no influence on the future is known as a Markov chain Definition 1.9.Markov chains have the property that for a random variable n X , ( ) ( ) are in the state space .This property is known as the Markov property.Definition 1.10.Let , x y be in the state space.The conditional probabilities ( ) x y in the sample space is known as the transition function.
From the transition function ( ) , p x y we construct a square matrix of the form    P is known as the transition matrix for the Markov chain.We will now apply the concepts of a Markov chain directly to our model.

Markov Chains and Our Model
We will use Markov chains to model the transition states for spin 1/2, spin 1, or spin 3/2 particles.For our model of N-spin particles the state space Ψ will be composed of all k Ψ with different magnetism numbers.We begin by defining the probability function for our model.
where N is total number of particles in the system, and for spin 1/2, spin 1, spin 3/2 particles respectively.Since our Markov chain has the property of only transferring to preferred states, our transition probabilities are defined by the following transition function.

( ) ( ) ( )
This transition function shows how transition probabilities are absorbed when we move to a preferred state from a less desired state.This probability absorption allows us to construct upper triangular square transition matrices of the following form.
( ) ( ) We will now look at Markov chains and the spin 1/2 particle.

Two Spin One Half Particles
For a two spin one half particles, the only possibilities of spin configuration of the fermions look like this , , − Ψ Ψ Ψ make up the state space of the two fermion system with our underlying 3 assumptions.Giving us the following probability transition table For example we would read ( ) | 0 P Ψ Ψ =; it is our assumption 2 that makes all of our probability matrices upper triangular.Which for our model denotes that system is probabilistically approaching the idealized ground state.Since read ( ) Ψ is an absorbing state, physically it shows that it is the idealized ground state.

Three Spin One Half Particles
For a system of three spin one half particles, the only possibilities of spin configuration of the fermions look like this make up the state space Ψ of the three fermion system with our underlying 3 assumptions.
Although we can clearly see what the probability of landing in each state is, we can use basic counting arguments to show our results mathematically.Thus • ( ) Giving us the following probability transition table Given the transition matrix for the three spin one half particle system, above, we can calculate the probablity of moving to a preferred state.
• ( ) With our understanding of the three spin one half particle case and counting arguments, we can develop a generalization for any N spin one half particle system.

N-Spin One Half Particle System
Now given our understanding of the probabilities we can build the probability transition matrix for any N-fer- Leading us to the following definition.

Definition 3.2. For any N-fermion system we can calculate ( )
| m m P Ψ Ψ where m Ψ is not the idealized ground state.Let 0 k > ∈  be the number of spin up or spin down particles for magnetism number m.Then ( ) For any N-fermion system we can calculate ( ) Ψ is not the idealized ground state.
( ) With these result we can calculate the probability of our system moving to a better state, given any N-fermions and m Ψ .We now turn our attention to systems of spin one particles.

Spin One Particles
For the spin one particles we will assign { } 1, 1, 0 + − for , , ↑ ↓ → respectively.The assumptions to the model will not change.The goal is examine N-spin one systems, and to see how their transition matrices look.

System of Two Spin One Particles
A system with two spin one particles has the following states: S ψ = − Thus for a system of two spin one particles, the state space Ψ has possible states 2 Ψ , which give us the following transition table

System of Three Spin One Particles
For a system of three spin one particles the possible states are given by the macro-states 3 Ψ .Whose micro-state configuration looks like this.
Notice that listing out all appropriate micro states gets a bit cumbersome, but unlike the spin 1/2 particles, the spin one particles have 3 possible spins.So to calculate their probabilities with the aid multinomial probability formula.
( ) , , , Now with this formula we are able to calculate the probability of the macro states for a 4 spin one particle system.

4 Spin One Particle System
For a spin one system with 4 particles we find macro states n Ψ such that { } 0, 1, 2, 3, 4 n ∈ ± ± ± ± .We do this by calculating their degenerate states via the multinomial formula.Giving us the following results: We now generalize our results with the following formulas Definition 4.1.For any N-spin one particle system and a given macro state m Ψ s.t.m N < the following hold.

even m-odd, or N-odd m-even
And for the zero case we get the following.
( ) Ψ be a macro state, for any N spin one particle system, then the probability of moving to a preferred state is ( ) Now with this result we can calculate the probability of our system moving to a preferred state, given N-spin one particles and m Ψ .We now turn our attention to a system of spin 3/2 particles.

Spin 3/2 Particles
Now we turn our attention to spin 3/2 particles.We will de-note all the possible states a particle can take with the set To give the spin 3/2 particles integer coefficients we simply represent them in the set Once again we will hold the assumptions as the spin 1/2, and spin 1 cases to find the Probability transition matrix between macro states.

Two Spin 3/2 Particle System
A two spin 3/2 particle system has the following states With the following transition matrix

Three Spin 3/2 Particle System
For a system of three spin 3/2 particles we calculate the state space in the same manner as we did in the two spin 3/2 particle system.Giving us the following state space { }

Four Spin 3/2 Particle System
With the aid of the multinomial formula we are able to calculate the probabilities of the state space With D being the degeneracy matrix.Notice that we can construct any transition probability matrix from just understanding degeneracy matrix D. Thus given our systems it suffices to understand degeneracies in order to construct their transition probabilities.

Degeneracy Matrix Algorithm
As we saw in the previous section we can find all the transition probabilities from any j spin-N particle system under assumptions: • 0 T → , thus eliminating the effects of entropy 0 S ∴ → • ψ states gravitate towards more preferred states of lower magnetism number due to our model being an anti-paramagnet • Particles in ψ are indistinguishable.
With Probability transition matrix P To construct D let d be the degeneracy vector, composed from all the degeneracies form all the i Ψ states We will now use the degeneracy vector to construct the degeneracy matrix.

Algorithm for Constructing D
For simplicity we will construct D from our degeneracy vector which we will denote 1 , , , m a a a = d  as the base for constructing D. We can construct the entire degeneracy matrix D from our degeneracy vector 1 d Now with this generalized algorithm we can construct any degeneracy matrix D, which allows us to construct any probability transition matrix P.

Degeneracy Vector (d1) for Spin One Particles for n = 5, 6, 7, 8, 9, 10
From the preceding section, we see that we can construct the probability transition matrix of a system solely by means of the degeneracy vector.Now here are the degeneracy vectors for n = 5, 6, 7, 8, 9, 10.

Results from Numerical Simulations
Two systems are studied: a small system with 10 spin 1/2 particles giving 1024 possible states, and a larger system with 20 spin 1/2 particles, with about a million possible states.The system begins in the least optimal state, all spins pointing up, and at each move, a random trial state is generated.If the trial state is more optimal, meaning it has fewer spins pointing up, it is accepted.It is highly probable the first move will be accepted.After many moves, the likelihood of improvement diminishes.We consider three "times" to illustrate the relative movement slowing at large time: after 9, 81 and 729 moves, so the moves are evenly spaced in log time.The simulation is run one million times for the small system, and one hundred thousand times for the large system.The results are displayed in Figure 1 and Figure 2, and are consistent with the calculation of predicted number at each, as shown in Table 1 and Table 2,   Ten spin one half particles are simulated one million times.The x-axis is the number of spins pointing up, at various times.The y-axis is the number of replicas of the system where this was found, proportional to the probability.Three times are given, in different colors, receding into the page, as: red (9 steps) green (81 time steps) and blue (729 time steps).The distribution is roughly Gaussian, and the distribution moves toward lower numbers of spins pointing up over time, but slowly.Numeric values are given in Table 1, along with computed values from the transition matrices from the formula above.Histogram of number of states having various net spin as a function of time.The x-axis is the number of spins pointing up.The y-axis is the number of replicas of the system which have the given x-axis value of spins up, and so the y-value is proportional to probability.Twenty spin one half particles are simulated over 729 Monte Carlo steps.One hundred thousand replicas of the system are used.Three times are displayed: after nine steps (red), after 81 steps (green) and after 729 steps (blue).The system has just over one million possible states, and so unlike in figure one, above, 729 steps are not nearly enough to get any significant fraction of the system into the optimal (x-axis value zero) state.However, despite the odds of moving to the optimal state being somewhat worse then one in a million, after only 729 moves, the system is pretty close to optimal: most likely the system has three or four spins pointing up out of 20.Beyond this time, progress is slow.

Definition 1 . 2 . 3 . 5 .
For any two states i For any N-particle system.i ψ will be the idealized ground state iff i The cardinality of the set n Ψ , denoted n Ψ is the total number of states with magnetism number n.
These probabilities give us the following transition table. below.

Figure 1 .
Figure 1.Histogram of number of states having various net spin as a function of time.Ten spin one half particles are simulated one million times.The x-axis is the number of spins pointing up, at various times.The y-axis is the number of replicas of the system where this was found, proportional to the probability.Three times are given, in different colors, receding into the page, as: red (9 steps) green (81 time steps) and blue (729 time steps).The distribution is roughly Gaussian, and the distribution moves toward lower numbers of spins pointing up over time, but slowly.Numeric values are given in Table1, along with computed values from the transition matrices from the formula above.

Figure 2 .
Figure 2.Histogram of number of states having various net spin as a function of time.The x-axis is the number of spins pointing up.The y-axis is the number of replicas of the system which have the given x-axis value of spins up, and so the y-value is proportional to probability.Twenty spin one half particles are simulated over 729 Monte Carlo steps.One hundred thousand replicas of the system are used.Three times are displayed: after nine steps (red), after 81 steps (green) and after 729 steps (blue).The system has just over one million possible states, and so unlike in figure one, above, 729 steps are not nearly enough to get any significant fraction of the system into the optimal (x-axis value zero) state.However, despite the odds of moving to the optimal state being somewhat worse then one in a million, after only 729 moves, the system is pretty close to optimal: most likely the system has three or four spins pointing up out of 20.Beyond this time, progress is slow.

Definition 2.1. Let
Our understanding of the transition function and the transition matrix of our model leads to the following 2 propositions .
The two particle system gives us the following states gives us the following degeneracy transition table, with