 at each site of the spin chain[^2]. In that case we denote the spin chain Hilbert space as $$\label{Hilbert.space.inhom}
\mathcal{V}_m=\mathbb{C}^2_{m_1} \otimes \mathbb{C}^2_{m_2} \otimes  \ldots \otimes \mathbb{C}^2_{m_N}\,,$$ and the Bethe equations turn into $$\label{eq:Bethe.eqn2}
\prod_{i=1}^N\frac{\sigma_a-m_i+\frac{\hbar}{2}}{\sigma_a-m_i-\frac{\hbar}{2}} = q\,\prod_{b\neq a}\frac{\sigma_a-\sigma_b+\hbar}{\sigma_a-\sigma_b-\hbar}\,, \qquad a=1,\ldots,k.$$ For a given number of excitations $k$, there are $\binom{N}{k}$ solutions distinct solutions of the Bethe equations. The solutions $\sigma_a^I$ can be labelled by a subset $I=\{I_1,\ldots,I_k\}\subset\{1,\ldots,N\}$ such that, expanding around $q \to 0$, the solutions are of the form $\sigma_a^I=m_{I_a}-\tfrac{\hbar}{2}+\cO(q)$.

R-matrices
----------

There are many independent ways to arrive at the Bethe equations. Usually, the most powerful method is the algebraic Bethe ansatz, which is based on the construction of an R-matrix. For the inhomogeneous spin chain, this is an operator acting on two sites, $$R_{ij}(m_j-m_i):\mathbb{C}^2_{m_i} \otimes \mathbb{C}^2_{m_j} \to \mathbb{C}^2_{m_j} \otimes \mathbb{C}^2_{m_i}\,.$$ It has rational dependence on $m_j-m_i$ and satisfies the regularity property $R_{ij}(0)\sim \mathbb{P}_{ij}$ where $\mathbb{P}_{ij}$ is the permutation operator, together with the Yang-Baxter equation (shown graphically in figure \[Fig:YangBaxter.equation\]) $$\label{Yang.Baxter}
R_{12}(m_2-m_1)R
at each site of the spin chain[^2 ]. In that case we denote the tailspin range Hilbert outer space as $ $ \label{Hilbert.space.inhom } 
 \mathcal{V}_m=\mathbb{C}^2_{m_1 } \otimes \mathbb{C}^2_{m_2 } \otimes   \ldots \otimes \mathbb{C}^2_{m_N}\,,$$ and the Bethe equations turn into $ $ \label{eq: Bethe.eqn2 } 
 \prod_{i=1}^N\frac{\sigma_a - m_i+\frac{\hbar}{2}}{\sigma_a - m_i-\frac{\hbar}{2 } } = q\,\prod_{b\neq a}\frac{\sigma_a-\sigma_b+\hbar}{\sigma_a-\sigma_b-\hbar}\, , \qquad a=1,\ldots, k.$$ For a yield number of excitations $ k$, there be $ \binom{N}{k}$ solution distinct solution of the Bethe equations. The solutions $ \sigma_a^I$ can be label by a subset $ I=\{I_1,\ldots, I_k\}\subset\{1,\ldots, N\}$ such that, expanding around $ q \to 0 $, the solutions are of the class $ \sigma_a^I = m_{I_a}-\tfrac{\hbar}{2}+\cO(q)$. 

 radius - matrices 
 ---------- 

 There are many independent ways to arrive at the Bethe equations. normally, the most powerful method is the algebraic Bethe ansatz, which is base on the construction of an R - matrix. For the inhomogeneous tailspin chain, this is an hustler acting on two sites, $ $ R_{ij}(m_j - m_i):\mathbb{C}^2_{m_i } \otimes \mathbb{C}^2_{m_j } \to \mathbb{C}^2_{m_j } \otimes \mathbb{C}^2_{m_i}\,.$$ It has rational dependence on $ m_j - m_i$ and satisfies the regularity place $ R_{ij}(0)\sim \mathbb{P}_{ij}$ where $ \mathbb{P}_{ij}$ is the permutation operator, together with the Yang - Baxter equation (shown graphically in number   \[Fig: YangBaxter.equation\ ]) $ $ \label{Yang. Baxter } 
 R_{12}(m_2 - m_1)R
 at each site of the spin cmain[^2]. In that casg qe denmte ths spin cfain Hilbert space as $$\label{Hmlbeet.spaxe.inhom}
\mathcal{V}_m=\mathbb{Z}^2_{m_1} \otimes \mathbb{C}^2_{n_2} \otmmes  \ldots \otimxa \mathbn{E}^2_{m_N}\,,$$ ahf thz Uethe equations turn into $$\label{eq:Bethe.exn2}
\orld_{i=1}^N\frac{\sigma_a-m_i+\frac{\hbar}{2}}{\sigma_a-m_i-\fras{\hbar}{2}} = a\,\pgod_{b\neq a}\frac{\sygma_s-\figmz_b+\hbar}{\sigma_a-\sigma_b-\hbar}\,, \qquad a=1,\ldota,k.$$ For e given number pf excitations $k$, there are $\binlm{N}{k}$ solutions dishinct solutuons if the Bethe equations. The solutiohs $\sigma_a^I$ can be labelled by a subszt $I=\{I_1,\ldots,I_j\}\sybsfj\{1,\ldots,N\}$ such that, expanding avpund $q \to 0$, thr solutions arc of vhe dorm $\sigma_a^I=m_{I_a}-\tfrac{\huar}{2}+\cO(q)$.

R-matrices
----------

There are many iudependent ways to areice at the Betfw eduauiois. Hsuallj, tie most powsrful methoe is the algebraic Neeye ansatz, whidh is fafed on the construction of an R-matrix. Fmr fhe inhomogeneous spin xhain, this is an operwtor actigg on two sites, $$R_{ij}(m_j-m_i):\mathbb{C}^2_{m_i} \otimes \mathbb{C}^2_{m_j} \to \metfbb{E}^2_{n_j} \otkned \mathbb{C}^2_{m_i}\,.$$ It has rational dependence on $m_j-m_y$ amd satisfies the regularity pgolgrty $R_{ij}(0)\sim \majhbb{P}_{ij}$ whsre $\mathbb{P}_{ij}$ is tje permotatiob operatow, tobether with the Yang-Baxter wquation (shocn traphically in figbre \[Fig:YangBarter.eqoation\]) $$\label{Yang.Baxter}
R_{12}(m_2-m_1)R
at each site of the spin chain[^2]. case denote the chain Hilbert space \otimes \otimes \mathbb{C}^2_{m_N}\,,$$ and Bethe equations turn $$\label{eq:Bethe.eqn2} \prod_{i=1}^N\frac{\sigma_a-m_i+\frac{\hbar}{2}}{\sigma_a-m_i-\frac{\hbar}{2}} = q\,\prod_{b\neq a}\frac{\sigma_a-\sigma_b+\hbar}{\sigma_a-\sigma_b-\hbar}\,, \qquad For a given number of excitations $k$, there are $\binom{N}{k}$ solutions distinct solutions the Bethe equations. The solutions $\sigma_a^I$ can be labelled by a subset $I=\{I_1,\ldots,I_k\}\subset\{1,\ldots,N\}$ that, around \to the solutions are of the form $\sigma_a^I=m_{I_a}-\tfrac{\hbar}{2}+\cO(q)$. R-matrices ---------- There are many independent ways to arrive the Bethe equations. Usually, the most powerful method the algebraic Bethe ansatz, is based on the construction an For the spin this an operator acting two sites, $$R_{ij}(m_j-m_i):\mathbb{C}^2_{m_i} \otimes \mathbb{C}^2_{m_j} \to \mathbb{C}^2_{m_j} \otimes \mathbb{C}^2_{m_i}\,.$$ It has rational dependence on $m_j-m_i$ and satisfies regularity property where $\mathbb{P}_{ij}$ the operator, with the Yang-Baxter graphically in figure \[Fig:YangBaxter.equation\]) $$\label{Yang.Baxter} R_{12}(m_2-m_1)R
 at each site of the spin chain[^2]. IN that case wE denoTe tHe sPiN chaIn HiLbert space as $$\laBEl{HiLbert.space.inhom}
\mathcal{v}_m=\matHbB{c}^2_{m_1} \otIMeS \mathBb{C}^2_{m_2} \otiMEs  \LDOts \OtImEs \mAtHBb{c}^2_{m_N}\,,$$ anD thE Bethe eQuations tuRn iNtO $$\label{eq:BethE.EqN2}
\prod_{i=1}^N\fraC{\siGma_a-m_i+\frac{\hbAr}{2}}{\sIgma_a-m_I-\fRac{\HBar}{2}} = q\,\pRod_{B\neq a}\Frac{\siGMa_a-\sigMa_b+\hbar}{\siGmA_A-\sigma_B-\Hbar}\,, \qquAD A=1,\lDots,K.$$ For a given number oF ExCItations $k$, there Are $\binOm{n}{K}$ sOLUtiOns Distinct soLuTions OF the BetHE eQUATioNS. The solutions $\Sigma_a^I$ can bE LabElled bY a SubSEt $I=\{I_1,\ldOts,I_k\}\SuBSet\{1,\Ldots,N\}$ such tHat, eXpanding aRound $q \TO 0$, the solUTions arE of the ForM $\siGma_a^i=M_{I_A}-\tFraC{\hBAr}{2}+\co(Q)$.

R-MatRIceS
----------

There arE mAnY indePendENT WAys tO arRive At the bethe equationS. UsUallY, The Most pOwerfUl meThOd is tHe algeBraic beThe ansatz, which iS basEd on the coNstRuCtiOn Of an R-MAtrix. FOr tHe iNhomogeNeous spIN chAiN, THIs Is an operator acting On TWO sItes, $$R_{ij}(m_J-m_i):\matHBb{c}^2_{m_I} \Otimes \maThBb{C}^2_{M_j} \to \MAThbb{C}^2_{M_j} \otIMeS \mathbb{C}^2_{M_i}\,.$$ It haS RaTiOnal depEnDence oN $m_J-m_i$ And SatisFIes tHe reguLarity prOpertY $r_{ij}(0)\sim \mathbb{P}_{iJ}$ Where $\mathbb{P}_{iJ}$ Is THE pERmutAtiOn operator, tOgetHEr wiTh thE yaNg-BAXter eQuatiOn (SHoWN graphically in figurE \[FIg:YangbaxteR.equation\]) $$\labeL{Yang.BaxteR}
r_{12}(M_2-M_1)R
 at each site of the spinchain[^2].  In t hat  ca se  wedeno te the spin ch a in H ilbert space as $$\lab el{Hi lb e rt.s p ac e.inh om}
\ma t hc a l {V} _m =\ mat hb b {C }^2_{ m_1 } \otim es \mathbb {C} ^2 _{m_2} \otim e s \ldots \o tim es \mathbb{C }^2 _{m_N} \, ,$$ and t heBethe  equat i ons tu rn into $ $\ l abel{e q :Bethe. e q n2 }
\p rod_{i=1}^N\frac{ \ si g ma_a-m_i+\frac {\hbar }{ 2 }} { \ sig ma_ a-m_i-\fra c{ \hbar } {2}} =q \, \ p r od_ { b\neq a}\frac {\sigma_a-\ s igm a_b+\h ba r}{ \ sigma_ a-\si gm a _b- \hbar}\,, \ qqua d a=1,\ld ots,k. $ $ For a given n umberofexc itat i on s$k$ ,t her e  a re$ \bi nom{N}{k }$  s oluti onsd i s t inct  so luti ons o f the Bethe e qua tion s . T he so lutio ns $ \s igma_ a^I$ c an be  l abelled by a su bset  $I=\{I_1 ,\l do ts, I_ k\}\s u bset\{ 1,\ ldo ts,N\}$  such t h at,  e x p a nd ing around $q \to0$ , th e soluti ons ar e  o ft he form$\ sig ma_a ^ I =m_{I _a}- \ tf rac{\hba r}{2}+ \ cO (q )$.

R- ma trices 
- --- --- ---
T here  are m any inde pende n t ways to arri v e at the Beth e  e q u at i ons.  Us ually, themost powe rful me tho d  is t he al ge b ra i c Bethe ansatz, whi ch  is ba sed o n the constru ction of a n R -matrix.  For th e  inhomogeneous  spin  chain, th i s is anopera tor acti ng on two s ites, $$ R_{ ij} (m_ j-m _ i ): \mathbb{C}^2_ { m _i}\o times \ mat hbb{C}^ 2_{ m_j } \ to\m athbb{C}^ 2_{m_j}\o ti me s\ma thbb{ C }^2_{m_i }\ ,.$ $Ithas r a tional  depe nden ce  o n  $m _j-m_i$ an d sati sf ie s th e r eg ulari ty p r ope rty $R_ {ij}(0)\s im\ math bb {P }_{ij}$  where $\math bb {P}_{ij}$is  th e perm u t ation op erator, together with t h e Yang- Bax ter e quat ion (show n g raphic all y  in fi gure \ [Fig: Ya ngB a x ter.e q u at ion \] ) $$\label { Y ang .Baxt er }
R_ {12}(m_ 2-m_1)R
 at_each site_of the spin chain[^2]._In that_case_we denote_the_spin chain Hilbert_space as $$\label{Hilbert.space.inhom}
\mathcal{V}_m=\mathbb{C}^2_{m_1}_\otimes \mathbb{C}^2_{m_2} \otimes _\ldots \otimes \mathbb{C}^2_{m_N}\,,$$_and_the Bethe equations turn into $$\label{eq:Bethe.eqn2}
\prod_{i=1}^N\frac{\sigma_a-m_i+\frac{\hbar}{2}}{\sigma_a-m_i-\frac{\hbar}{2}} = q\,\prod_{b\neq a}\frac{\sigma_a-\sigma_b+\hbar}{\sigma_a-\sigma_b-\hbar}\,, \qquad a=1,\ldots,k.$$ For a given_number_of excitations_$k$,_there_are $\binom{N}{k}$ solutions distinct solutions_of the Bethe equations. The_solutions $\sigma_a^I$_can be labelled by a subset $I=\{I_1,\ldots,I_k\}\subset\{1,\ldots,N\}$ such_that,_expanding around $q_\to 0$, the solutions are of the form $\sigma_a^I=m_{I_a}-\tfrac{\hbar}{2}+\cO(q)$.

R-matrices
----------

There_are many independent ways to arrive_at the Bethe_equations._Usually,_the most powerful method_is the algebraic Bethe ansatz, which_is based on the construction of_an R-matrix. For the inhomogeneous spin chain,_this is an operator acting on_two sites, $$R_{ij}(m_j-m_i):\mathbb{C}^2_{m_i} \otimes \mathbb{C}^2_{m_j}_\to \mathbb{C}^2_{m_j}_\otimes \mathbb{C}^2_{m_i}\,.$$ It has rational_dependence on $m_j-m_i$_and satisfies_the regularity property_$R_{ij}(0)\sim \mathbb{P}_{ij}$ where $\mathbb{P}_{ij}$ is the_permutation operator, together_with the Yang-Baxter equation (shown graphically_in_figure \[Fig:YangBaxter.equation\]) $$\label{Yang.Baxter}
R_{12}(m_2-m_1)R