Comparison between XY Spin Chains with Spin 1/2 or 1 Interacting with Quantized Electromagnetic Field by One and Two Photon Jaynes-Cummings Model

: This paper describes two cases of interaction between a quantized electromagnetic ﬁeld and two different XY spin molecules; one with spins 12 , and the other with spins 1. Both interact with a quantized electromagnetic ﬁeld, with one of the spins in the chain interacting with the electromagnetic ﬁeld. The interaction between the ﬁeld mode and the spin chain with spins 1 is described by the one- and two-photon Jaynes-Cummings model (JC model). On the other hand, the interaction between the spins 1 2 and the electromagnetic ﬁeld is described only by the one-photon Jaynes-Cummings model. Analytical and numerical calculations were made for the case of a different number of photons in the ﬁeld mode, a different number of spins, and a different position of spin, interacting with the electromagnetic ﬁeld. The invariant and block structures of such a chain are shown with a comparison made between the evolution of the magnetic moment and the number of photons in both cases.


Introduction
The main task in quantum optics is the interaction between a quantized electromagnetic field and the energy levels of the atom or ion. Such an interaction is enshrined in the be Jaynes-Cummings quantum model [1][2][3]. The Jaynes-Cummings (JC) model describes the dipole transition between two levels in a quantum system. The model describes the transition between electronic levels of an atom due to absorption or emission of an electromagnetic quantum. This model is used to describe the transition between states of electron in an ion situated in an ion trap [4][5][6]; the manipulation of the state of a quantum dot [7,8].
The two-photon JC model with three-level atom has also been extensively studied [19]. The transitions between levels, from third to second and from second to first, are made by the one-photon JC model, with one photon emission. The two-photon JC model also includes two photon absorption/emission, and has been used for transferring the population between the first and third levels. The two-photon JC model was an important step in developing advanced methods for better population transfer, such as STIRAP (stimulated-Raman adiabatic passage) [20]. Nowadays, STIRAP has applications in many fields [21], such as the formation of ultracold molecules [21], the creation of gates for quantum information in different schemes [21], the control of nitrogen vacancy centers and semiconductor quantum dots [21], and the quantum optical analogues that are used in waveguide optics [21] and frequency conversion [21]. The two photon JC model has also been used to study the effects of the spin orbital coupling [22,23] of ultracold atoms located in the ring optical cavity. Such atoms play the role of pseudo-spins. Spin orbital coupling

Spin Chain with Spin 1/2 Interacting with the Electromagnetic Field
The system under consideration here is a linear spin chain consisting of N spins, located in the cavity. A laser field is applied to the cavity. The laser field is focused on the k-th spin. The cavity leads to quantization of the field, so let there be one mode with frequency c 2 . The electromagnetic field interacts with the z-th component of the spin. The difference between the two levels of the k-th spin is c 1 . An external magnetic field is applied to the cavity coinciding with the Z-detection of the spins. The magnitude of the external magnetic field is c 0 . All spins have spin 1 2 , and interact with their nearest neighbors through the X and Y components of their spin (described by the XY model).

Hamiltonian, Invariant, and Block Structure of the Hamiltonian
Let us first consider the case of a linear chain consisting of spins 1/2 interacting with a single mode of quantized electromagnetic field. The spin chain is open, and there is only first neighbor interaction. The interaction between the spins is in the XY model. This case has been studied in detail in the articles [42,43], and the system is shown in Figure 1.
Magnetochemistry 2021, 7, x FOR PEER REVIEW 3 of 18 spin chain. The last section, Section 3, is the conclusion, comparing the two spin chain cases.

Spin Chain with Spin 1/2 Interacting with the Electromagnetic Field
The system under consideration here is a linear spin chain consisting of spins, located in the cavity. A laser field is applied to the cavity. The laser field is focused on the k-th spin. The cavity leads to quantization of the field, so let there be one mode with frequency . The electromagnetic field interacts with the z-th component of the spin. The difference between the two levels of the k-th spin is . An external magnetic field is applied to the cavity coinciding with the Z-detection of the spins. The magnitude of the external magnetic field is . All spins have spin ½, and interact with their nearest neighbors through the X and Y components of their spin (described by the XY model).

Hamiltonian, Invariant, and Block Structure of the Hamiltonian
Let us first consider the case of a linear chain consisting of spins 1/2 interacting with a single mode of quantized electromagnetic field. The spin chain is open, and there is only first neighbor interaction. The interaction between the spins is in the XY model. This case has been studied in detail in the articles [42,43], and the system is shown in Figure 1. ( With N denoting the number of spins in the chain. is the interaction constant between spin and field mode. Exchange coupling constant of spin-spin interaction between the nearest neighboring spins (at position and 1) is J. Both constants and J can be positive or negative. The Hamiltonian can be block diagonalized using the invariant, as shown in [42,43].
The operators that change the number of photons in the mod of the field are creation operators of the mode and annihilation operators of the mode : In wave function | , … ⟩, first is written the number of photons in the field (denoted with n), second are written the states of spins in the chain, (can be ↑ for spin up and ↓ for spin down).
The state of the k-th spin (denoted by | , … … ⟩ ≡ | ⟩ ) is changed by the following operators: The Hamiltonian of this system is: With N denoting the number of spins in the chain. G 1 is the interaction constant between spin and field mode. Exchange coupling constant of spin-spin interaction between the nearest neighboring spins (at position i and i + 1) is J. Both constants G 1 and J can be positive or negative. The Hamiltonian can be block diagonalized using the invariant, as shown in [42,43].
The operators that change the number of photons in the mod of the field are creation operators of the modeâ † and annihilation operators of the modeâ: In wave function |n, s 1 . . . s N , first is written the number of photons in the field (denoted with n), second are written the states of spins in the chain, s i (can be ↑ for spin up and ↓ for spin down).
The state of the k-th spin (denoted by |n, s 1 . . . s k . . . s N ≡ |s k ) is changed by the following operators: For each position in the chain i ∈ [1, N − 1], the operatorsŜ + iŜ − i+1 andŜ − iŜ + i+1 turn down the state of one spin, and up the state of the other spin. The total number of spins oriented upwards in the system remains constant. The operatorsâ †Ŝ− k andâŜ + k turn the spin state down, and one further photon is emitted, or the number of photons is reduced when one spin is turned up. I is an imaginary unit.
The parts of the Hamiltonian: do not change the state of the spin or the number of photons. The invariant is defined by Ĥ ,Înv = 0. Then the invariant in case of spin 1/2 of the system isÎnv 1/2 [54]:Î whereŜ z i is the Pauli matrix (of the i-th spin) in the Z direction, andn =â †â is the operator for number of photons in the field's mode.
When the operator,Înv 1/2 , is applied to one of the Hamiltonian eigenstates, | ψ(t i ) , the corresponding eigenvalue, µ k , is whole-or half-integer. The eigenvalue, µ k , makes sense of the total number of excitations of the state [42].
where | ψ(t i ) are the states of the system at time t i . The application of Hamiltonian to each state results in another state with the same invariant value,Înv 1/2 . By reordering the basis vectors of the matrix, a block diagonal matrix can be obtained with each block corresponding to a different value of the invariant whereĤ(µ w ) are block matrices and each block corresponds to a different value of the invariant µ w . Index w is the integer and corresponds to the block number. The block size for different invariant values is: The maximum value of Dim Ĥ (µ w ) is obtained when k = N. Then, as the value of w increases, the block size remains constant. The matrix elements of these block matrices can be calculated as follows: where | n 1 , s 1 s 2 . . . s i . . . s N and n 2 , z 1 z 2 . . . z j . . . z N are two states of the system. The states have n 1 and n 2 photons, respectively. For vectors s i and z j , the states of the i-th spin of the system can be ↑ or ↓.

Evolution of Magnetic Moment and Number of Photons
One way to calculate the evolution is by using a full Hamiltonian to obtain an evolution operator,Ê. The task can be simplified by taking into account the fact that during the state evolution, the total number of excitations remains constant. Therefore, the evolution of each state can be found using the block matrix with the same number of excitations as the initial state: Diagonalizing the Hamiltonian makes it easy to find the operator of evolution. The diagonal elements of the Hamiltonian correspond to the energy spectrum. Diagonalization also helps to find the eigenstates of the system, and shows degeneracy.
Using the evolution operator, the state of the system can be calculated from the formula: where | ψ(t i ) is the current state of the system in the time t i .
whereM z is the magnetic moment operator in the Z direction. Mz(t i ) is the magnetic moment of the chain in the Z direction, when the system is in state | ψ(t i ) . The number of photons in the mode of the field, in the state of the system | ψ(t i ) is n(t i ). Figure 2 shows a numerical simulation of the evolution of the magnetic moment and the number of photons of the spin chain interacting with a quantized electromagnetic field. The spin chain consists of 1 spin in Figure 2a The Hamiltonian connects all states that have the same number of excitations. In the case of more than one spin, there is an interaction between spins. The spin flipping part by itself does not change the number of photons in the mode (or magnetization), but when the spin chain consists of more than one spin, the excitations can pass between spins in the chain. It can be seen that increasing the number of spins increases the number of evolution-connected states. The interferences that occur lead to a difference in the evolution of the number of photons in the mode for more spins than the evolution of the number of photons in the mode in the case of only one spin. During evolution, the system can be in superposition of the states. The magnitude of the spin-spin interaction determines how much each spin in the chain contributes to magnetization, depending on the distance of the spin interacting with the field. These are the reasons why magnetization has continuous spectrum.  The Hamiltonian connects all states that have the same number of excitations. In the case of more than one spin, there is an interaction between spins. The spin flipping part of the Hamiltonian ∑ by itself does not change the number of photons in the mode (or magnetization), but when the spin chain consists of more than one spin, the excitations can pass between spins in the chain. It can be seen that increasing the number of spins increases the number of evolution-connected states. The interferences that occur lead to a difference in the evolution of the number of photons in the mode for more spins than the evolution of the number of photons in the mode in the case of only one spin. During evolution, the system can be in superposition of the states. The magnitude of the spin-spin interaction determines how much each spin in the chain contributes to magnetization, depending on the distance of the spin interacting with the field. These are the reasons why magnetization has continuous spectrum.
In the case of one spin, the Hamiltonian connects the following states, | , ↑⟩ and | 1, ↓⟩. There is only one excitation that can be exchanged between the spin "chain" and the mode. The invariant value is 0.5. This is the typical Rabi oscillation, well known from quantum optics.  In the case of one spin, the Hamiltonian connects the following states, | n, ↑ and | n + 1, ↓ . There is only one excitation that can be exchanged between the spin "chain" and the mode. The invariant value is 0.5. This is the typical Rabi oscillation, well known from quantum optics. Figure 2c,d shows two spins and two excitations. In this case the invariant value is 1. The eigenvectors connected during evolution are | 2, ↓↓ , | 1, ↓↑ , | 1, ↑↓ and | 0, ↑↑ . The states with the same number of photons | 1, ↓↑ and | 1, ↑↓ have the same value of Mz equal to zero. The states | 0, ↑↑ and | 2, ↓↓ have magnetization values 1 and −1, respectively.

Spin Chain with Spin 1 Interacting with the Electromagnetic Field
The system of our study is a linear spin chain consisting of N spins, all of them are with spin 1. All spins are located in a cavity, and the laser field is focused on the k-th spin, when the field interacts with the z-th component of the spin. The electromagnetic field is quantized in the cavity, and has a frequency of c 2 . An external magnetic field is applied to the cavity, coinciding with the Z detection of the spins. The magnitude of the external magnetic field is c 0 . The k-th spin between neighboring levels has the same transition frequency, c 1 . Spins interact with their closest neighbors through the X and Y components of their spin.

Hamiltonian
A spin chain consisting of spin 1 particles interacts with one or two photons, and one of the spins interacts with an electromagnetic field mode. The Hamiltonian of such a system is: whereĤ XY is the Hamiltonian of spin-spin interaction through the XY model. The spin chain is open, and each spin only interacts with its first neighbors. The Hamiltonian of the two-photon Jaynes-Cummings model,Ĥ JC2 , in the system Hamiltonian,Ĥ, describes the interaction between one spin in the chain and two photons in one mode from a quantized electromagnetic field: Here,Ĥ 1P andĤ 2P correspond to interaction between the k-th spin and the field mode, whileĤ 1P corresponds to the one-photon absorption, andĤ 2P corresponds to the twophoton absorption effect. In this case, the spin has to absorb two photons from the mode to change its state to two levels higher, and accordingly to emit the corresponding number of photons to go to a lower level. G 1 and G 2 are the interaction constants between the k-th spin and field mode. Both constants, G 1 and G 2 , can be positive or negative. The constant c 1 is the energy difference between the two states of the k-th spin, and c 2 is the energy carried from the photons in the mode of the field.â † andâ are creation and annihilation operators of one photon in the mode, correspondingly, and the creation and annihilation operators for k-th spin areς + k andς − k . The Pauli matrix analogue,ς Z k , corresponds with the spin Z component. Concerning the Z axis, the other two axes, X and Y, are equivalent. From a mathematical point of view, operatorsς + k andς − k are a linear combination of the Pauli matrix analogues:ς A Hamiltonian part describing interaction between spins, according to the XY model, isĤ XY . Another Hamiltonian part describes an external magnetic field,Ĥ EF .
where k is the position of the spin interacting with the quantized electromagnetic field; constant c 0 corresponds to the external magnetic field, andς Z i is the Pauli matrix analogue that corresponds to the third component of the spin. Here, J is the exchange coupling constant of the spin-spin interaction between the i-th and (i + 1)-th spin. Constant J can be positive or negative. The spin chain consists of N spins. Illustrative example for such a system is shown in Figure 1.

Invariant and Block Structure
An important feature of the model describing the chain of particles with spin 1 interacting with the quantized electromagnetic field is the invariant. By definition, the invariant is: (27) and The Hamiltonian (Equation (21)) retains the total number of excitations during evolution, and has the following invariant: i is the magnetic moment in the Z direction, andn =â †â is number of photons in the mode of the field. WhenÎnv 1 is applied to one of the Hamiltonian eigenstates, the corresponding eigenvalue, µ i , is an integer. The operators,ς Z i , are the Z-th Pauli matrix analogue for spin 1. In this form the invariant does not match the number of excitations in the system. The number of excitations can be found by applying the following operator to the state: Ex 1 gives the number of excitations of the system, and is also an invariant of the system. The analogs of the Pauli matrix for spin 1 are: Here, the spin projection in direction Z is considered. The eigenvalues of the Pauli matrix for spin 1 are −1, 0, and 1. The states will be denoted as ↓, →, and ↑, respectively. By using the invariant, the Hamiltonian can be block diagonalized by analogy with the 1/2 spin systems. The maximal size of the block is N 2 . After reaching this maximum amount, the block matrix size remains constant when increasing the number of photons in the mode. The formula for the block matrix size of a spin chain, with spin 1 interacting with electromagnetic field, can be obtained analytically. The expression for the block matrix size is the following: In this formula n is the number of the spins, µ w is the number if the excitations, and m 1 = |µ w − j 0 − 2j 2 | is the number of the spins in the state zero (one excitation). The number of spins with zero excitations, m 0 (spin state is −1), and the number of spins with two excitations is m 2 (spin state is 1). The brackets and round the number in them down. The condition for the number of spins is obtained by m 1 + 2m 2 ≤ N.
The matrix elements of the Hamiltonian can be found as follows: Similarly, to the spin 1 2 case, block matrices have the same size and the corresponding positions of the matrix elements, but the elements can have different values, no matter which spin interacts with the photon.
The Equations (35)-(37) describe examples of block matrices of the open two-spin chain with spin 1, where one of the spins interacts with the quantized electromagnetic field. The interaction is only between the closest neighbors J i,i+1 = J and J i,j =i+1 = 0. The number of excitations is 1, 2, 3, and 4, with blocks sizes 1, 3, and 6, respectively: The block matrices reach their maximum size when the invariant value is equal to 2N. The maximum size can be calculated using the number of spin states (in this case equal to three) and the number of spins in the chain, 3 N .
In the same way, as is shown in Section 2.1.2., evolution can be calculated using the full Hamiltonian or using only its block corresponding to the relevant number of excitations and equal to that of the initial state of the system: System states can be calculated using the operator of an evolution formula: where |ψ(t i ) are the system states at time t i ;M z (t i ) is the magnetic moment the operator in the Z direction, and the same moment in time t i andn is the number of photons in the field mode. Each state of the system has its magnetic moment in the Z direction, Mz(t i ), and number of photons, n(t i ). They can be computed using equations:

Numerical Simulation for Energy Spectrum
The eigenvalues of the Hamiltonian can be used for facilitated computation of the evolution operator. They also show connections between eigenvalues with changing parameters. The following parameterizations, G 1 = cos (ϕ), G 2 = cos 2 (ϕ), and J i,i+1 = sin(ϕ), can be used. This description of these constants is general, since arbitrary values of G and J can be written as G = A.cos (ϕ) and J = A.sin (ϕ). Here the amplitude is The Hamiltonian in the case of single-photon absorption can be obtained from Equations (19)-(26), with G 1 = 0 and G 2 = 0. Similarly, in the case of two-photon absorption, they can be obtained from the same equations, when G 1 = 0 and G 2 = 0. In both cases the invariant is the same and is equal to the one, where G 1 = 0 and G 2 = 0. Important limiting cases are when there is only one spin or when there is no spin-spin interaction between the spins. This case can be reduced to a Jaynes-Cummings model of a two or three-level system, similar to a two or three-level atom. In the case of singlephoton absorption, the energy for the transition between levels 1→2 and 2→3 is equal. The two-photon simultaneous absorption case without spin-spin interaction can be effectively reduced to one-photon absorption for spin 1 2 . Figures 3-6 show the energy spectrum for spin chains consisting of 4 spins and a second spin interacting with a quantized electromagnetic field. Each figure corresponds to a different number of excitations: Figure 3 for one excitation, Figure 4 for two, Figure 5 for three, and Figure 6 for four. In the figures, Figures 3a, 4a, 5a and 6a corresponds to a spin chain consisting of spins 1 2 . Figures 3a, 4b, 5b and 6b correspond to a spin chain consisting of spins 1 interacting when only one photon emission/absorption is present. Figures 3b, 4c ,

Numerical Simulation for Energy Spectrum
The eigenvalues of the Hamiltonian can be used for facilitated computation of the evolution operator. They also show connections between eigenvalues with changing parameters. The following parameterizations, cos φ , cos φ , and , sin φ , can be used. This description of these constants is general, since arbitrary values of G and J can be written as . cos φ and . sin φ . Here the amplitude is . The Hamiltonian in the case of single-photon absorption can be obtained from Equations (19)-(26), with 0 and 0. Similarly, in the case of two-photon absorption, they can be obtained from the same equations, when 0 and 0. In both cases the invariant is the same and is equal to the one, where 0 and 0. Important limiting cases are when there is only one spin or when there is no spinspin interaction between the spins. This case can be reduced to a Jaynes-Cummings model of a two or three-level system, similar to a two or three-level atom. In the case of singlephoton absorption, the energy for the transition between levels 1→2 and 2→3 is equal. The two-photon simultaneous absorption case without spin-spin interaction can be effectively reduced to one-photon absorption for spin ½. Figures 3-6 show the energy spectrum for spin chains consisting of 4 spins and a second spin interacting with a quantized electromagnetic field. Each figure corresponds to a different number of excitations: Figure 3 for one excitation, Figure 4 for two, Figure 5 for three, and Figure 6 for four. In the figures, Figures 3a, 4a, 5a and 6a corresponds to a spin chain consisting of spins ½. Figures 3a, 4b, 5b and 6b correspond to a spin chain consisting of spins 1 interacting when only one photon emission/absorption is present. Figures 3b, 4c, 5c and 6c    The number of excitations is equal to one. (a) Results for the spectrum are the same for the case of spin chain consist of spins 1/2; the case of spin chain consists of spin 1, G 1 = cos (ϕ), and G 2 = 0; the case of spin chain consists of spin 1, G 1 = cos (ϕ) and G 2 = cos 2 (ϕ) (b) spin chain consists of spin 1, G 1 = 0 and G 2 = cos 2 (ϕ).
In the case of one excitation, the energy spectrum of both cases of spin 1 2 and spin 1 (with only one photon and with one and two photon interaction) is the same, because there is not enough excitation for spin to go to state 1. When there is only one spin interacting with the electromagnetic field, and there are no other spins interacting with the first one, the system only has two states, |1, ↓ and |0, → .
In the case of more excitations, when J = G 1 = 0 and G 2 = 0 the eigenvalues are close to those with spin 1 2 chain, with two times more excitations and changed values of parameters of spin 1 2 Hamiltonian-c 1 , c 2 , and G 1 . When J = 0 and G 1 = G 2 = 0 the energy spectrum is the same as that of the spin 1 2 case, with the same number of excitations. With a change of ϕ between ±π/2 and 0, these states "connect" with each other. Differences in energy spectrum occur due to different numbers of excitations transported by different parts of the Hamiltonian. The XY interaction between spins in the chain transports only one excitation between the spins at a time. However, the spin-light interaction transports "quanta" of two excitations at once.     In the case of one excitation, the energy spectrum of both cases of spin ½ and spin 1 (with only one photon and with one and two photon interaction) is the same, because there is not enough excitation for spin to go to state 1. When there is only one spin interacting with the electromagnetic field, and there are no other spins interacting with the first one, the system only has two states, |1, ↓⟩ and |0, →⟩.
In the case of more excitations, when 0 and 0 the eigenvalues are close to those with spin ½ chain, with two times more excitations and changed values of parameters of spin ½ Hamiltonian-, , and . When 0 and 0 the energy spectrum is the same as that of the spin ½ case, with the same number of excitations. With a change of φ between π/2 and 0, these states "connect" with each other. Differences in energy spectrum occur due to different numbers of excitations transported by different parts of the Hamiltonian. The XY interaction between spins in the chain transports only one excitation between the spins at a time. However, the spin-light interaction transports "quanta" of two excitations at once. From a mathematical point of view cos φ and sin φ have the same value √2/2 when φ π/4, so that both cos φ and sin φ have a similar contribution when only one photon interaction is present (as in the case of spin 1 and spin ½ chain). When there are only two photon absorptions, cos φ is much smaller than cos φ , but the contribution of the spin-spin interaction is the same as in the case when only two photon absorptions are present.

Numerical Simulation of Magnetic Moment and the Number of Photons
The invariant value and the number of excitations do not change during evolution. They are simply offset from each other, and are spaced equal to lnv=|n-N|. This is obtained from Equations (30) and (31) for the invariant of the system.
The difference between the magnetic moment behavior and the number of photons for spins ½ and 1 comes from the different invariant equations (6 and 30, respectively). The absorption of a photon by a particle chain with spins 1 changes the magnetic moment of the spin chain by half the difference from the magnetic moments at the two opposite spin states (Figure 7). In the case of a spin chain with spins ½, upon absorption of a photon, Figure 6. Energy spectrum of four spins in a chain and a second spin interacting with the quantized electromagnetic field. The number of excitations is equal to four. (a) Spin chain consist of spins 1 2 ; (b) Spin chain consists of spin 1, G 1 = cos (ϕ) and G 2 = 0; (c) Spin chain consists of spin 1, G 1 = 0 and G 2 = cos 2 (ϕ); (d) Spin chain consists of spin 1, G 1 = cos (ϕ) and G 2 = cos 2 (ϕ).
From a mathematical point of view cos (ϕ) and sin (ϕ) have the same value √ 2/2 when ϕ = π/4, so that both G 1 = cos (ϕ) and J = sin (ϕ) have a similar contribution when only one photon interaction is present (as in the case of spin 1 and spin 1 2 chain). When there are only two photon absorptions, G 2 = cos 2 (ϕ) is much smaller than G 1 = cos (ϕ), but the contribution of the spin-spin interaction is the same as in the case when only two photon absorptions are present.

Numerical Simulation of Magnetic Moment and the Number of Photons
The invariant value ψ(t i ) Î nv ψ(t i ) and the number of excitations do not change during evolution. They are simply offset from each other, and are spaced equal to lnv=|n-N|. This is obtained from Equations (30) and (31) for the invariant of the system.
The difference between the magnetic moment behavior and the number of photons for spins 1 2 and 1 comes from the different invariant equations (6 and 30, respectively). The absorption of a photon by a particle chain with spins 1 changes the magnetic moment of the spin chain by half the difference from the magnetic moments at the two opposite spin states (Figure 7). In the case of a spin chain with spins 1 2 , upon absorption of a photon, the magnetic moment changes with the difference from the magnetic moments of the two opposite states (Figure 2). Figure 7 shows a numerical simulation of the evolution of the number of photons (a, c, e) and the magnetic (b, d, f) moment of a spin chain, consisting of different numbers of spins, interacting with an electromagnetic field. In the case of more than two spins, the electromagnetic field interacts with the second spin. Figure 7a During evolution the Hamiltonian connects all states that have the same number of excitations. As in the case of spin ½, the interaction between spins 1 does not change the number of photons or the magnetization. Only spin-mode interaction in the JC model changes the number of photons. However, the chain interferences occur again when the excitations "go through". Increasing the number of spins increases the number of states connected through the evolution. The system can be in a superposition of states, and the magnitude of the spin-spin interaction gives a continuous spectrum magnetization and number of photons. show the magnetic moment, Mz, along the Z axis of the system at time t during the evolution of the system. The red line corresponds to the number of photons and the green to the magnetic moment in the case of single-photon and two-photon absorption simultaneously. The purple and teal lines correspond to the case where there is only two-photon absorption. The brown and blue lines correspond to the case where there is only a one photon absorption. The first two (a,b) correspond to an initial state |2, ↓↓ , the next two (c,d) to an initial state |3, ↓↓ . The remaining two correspond to an initial state |3, ↓↓↓↓ .
During evolution the Hamiltonian connects all states that have the same number of excitations. As in the case of spin 1 2 , the interaction between spins 1 does not change the number of photons or the magnetization. Only spin-mode interaction in the JC model changes the number of photons. However, the chain interferences occur again when the excitations "go through". Increasing the number of spins increases the number of states connected through the evolution. The system can be in a superposition of states, and the magnitude of the spin-spin interaction gives a continuous spectrum magnetization and number of photons.
In the case when only a two-photon spin-light interaction is included from the JC model, the excitation transfer between the mode and the field very much depends on the number of photons in the mode. The amplitude of the two photon interaction with the mode is much smaller compared to one photon interaction, because G 2 = G 2 1 < G 1 . When the number of photons in the field is one, it plays no role in evolution. Since these terms change the number of photons only by two, there is no state connected with the state with one excitation in the mode. When the number of the excitations in the system is two or more, both spin-light interaction parts of the Hamiltonian govern the evolution.
As the number of photons in the field increases, the excitations transfer coming from this interaction also increases, as shown in Figure 8. The Hamiltonian part of the two photon JC model does not play any role in evolution when there is only one photon in the mode. That is the reason why, when the only spin-light interaction is through a one photon JC model and when both a one and two photon JC model interaction are present, they have the same evolution. This is shown in Figure 8a for a spin chain consisting of 5 spins. Evolution of the system when the both spin-light terms of the Hamiltonian are included, at the beginning is close to the evolution with only one photon interaction. With a longer evolution, however, the difference between them becomes significantly larger. As the number of photons in the mode increases, the two photon interaction part of the JC Hamiltonian has a bigger impact on the overall evolution of the system. This can be seen in Figure 8b for a spin chain consisting of 5 spins. The same is also true for a changing magnetization due to the interaction of the field with the spin chain. Figure 8 shows the evolution of the number of photons for a 5-spin system. The first spin interacts with the photons. The other parameters , , , , , and , are the same as in the previous example. The initial state is when all the spins are down, and the number of excitations is equal to the number of photons. It can be seen that when there are a small number of photons, the single-photon absorption has a major influence on evolution. As the number of photons increases, the effect of two-photon absorption also increases.  Figure 8. The number of photons, n ,in the system at time t during the evolution of the system is shown. The first spin of the system interacts with the electromagnetic field. Purple, brown, and red correspond to cases where there is only twophoton absorption, only one-photon absorption, and absorption in both cases, respectively. (a) corresponds to 1 photon, (b) corresponds to 2 photon, (c) corresponds to 3 photon, (d) corresponds to 4, (e-l) corresponds to 5-12 photon.

Conclusions
This paper compares XY spin chains with a different spin; ½ and 1 interacting with a quantized electromagnetic field. In both cases, there is an invariant corresponding to the number of excitations. In the case of spin 1 three cases of Hamiltonian are shown: First, with only a one photon spin-light interaction between the mode and the chain. Second, only a two photon spin-light interaction is present. Third, both a one photon and two spin Figure 8. The number of photons, n, in the system at time t during the evolution of the system is shown. The first spin of the system interacts with the electromagnetic field. Purple, brown, and red correspond to cases where there is only two-photon absorption, only one-photon absorption, and absorption in both cases, respectively. (a) corresponds to 1 photon, (b) corresponds to 2 photon, (c) corresponds to 3 photon, (d) corresponds to 4, (e-l) corresponds to 5-12 photon. Evolution of the system when the both spin-light terms of the Hamiltonian are included, at the beginning is close to the evolution with only one photon interaction. With a longer evolution, however, the difference between them becomes significantly larger. As the number of photons in the mode increases, the two photon interaction part of the JC Hamiltonian has a bigger impact on the overall evolution of the system. This can be seen in Figure 8b for a spin chain consisting of 5 spins. The same is also true for a changing magnetization due to the interaction of the field with the spin chain. Figure 8 shows the evolution of the number of photons for a 5-spin system. The first spin interacts with the photons. The other parameters c 0 , c 1 , c 2 , J, G 1 , and G 2 , are the same as in the previous example. The initial state is when all the spins are down, and the number of excitations is equal to the number of photons. It can be seen that when there are a small number of photons, the single-photon absorption has a major influence on evolution. As the number of photons increases, the effect of two-photon absorption also increases.

Conclusions
This paper compares XY spin chains with a different spin; 1 2 and 1 interacting with a quantized electromagnetic field. In both cases, there is an invariant corresponding to the number of excitations. In the case of spin 1 three cases of Hamiltonian are shown: First, with only a one photon spin-light interaction between the mode and the chain. Second, only a two photon spin-light interaction is present. Third, both a one photon and two spin light interaction are present. The invariant is defined by the number of excitations, and can be obtained by using the number of photons in the mode and the number of excited spins, and it is the same in all three cases.
Using the invariant, the Hamiltonian can be diagonalized in blocks. The block structure is independent of the position of the interacting spin, and can be used to calculate the evolution of the state with the corresponding number of excitations. The different invariant of the chains with spin 1 2 and 1 leads to a difference in the evolution, and a corresponding difference in photon absorption and different magnetic moment.
The energy spectrum is obtained for a Hamiltonian of spin 1 chain, were energy is found for different values of spin-spin and spin-light interaction terms. By using the energy spectrum, it can be seen that when there is only one photon in the mode, the two photon interaction part of the Hamiltonian does not play any role. The contribution of the two photon spin-light interaction term increases with increasing the number of photons. The spin-light interaction in the Hamiltonian can be effectively reduced to single-photon absorption, such as in the case of a spin 1 2 chain; with different values of other parameters when there is only two-photon absorption and spin-spin interaction is not present. The addition of a spin-spin interaction leads to a different evolution, as the spins exchange only one excitation, but the spin-light interaction exchanges two.
The interaction of two photons with spin-light plays a much smaller role at the beginning of the evolution, when the amplitude of the two photon spin-light interaction is smaller than the amplitude of the one photon spin-light interaction and the number of photons in the mode is small. The role of the two photon spin-light interaction increases when the number of photons in the mode increases.