Spin-textures of medium-body boson systems with trapped spin-f cold atoms

The spin-textures of bound medium-body systems with spin-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {f}$$\end{document}f atoms (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {f}\ge 3$$\end{document}f≥3) have been studied. The Hamiltonian is assumed to be dominated by the two-body interaction favoring parallel spins. The system with particle number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N=8$$\end{document}N=8 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {f}=3$$\end{document}f=3 is first chosen, and the Hamiltonian is exactly diagonalized by using Fock-states as basis-states, thereby all the eigenenergies and eigenstates are obtained and a detailed analysis is made. Then the cases with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N=13$$\end{document}N=13 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {f}=4$$\end{document}f=4 are further studied. Since the total spin S is conserved, the eigenstates having the same S form an S-group. Let the lowest (highest) energy state of an S-group be called a bottom-state (top-state). We found that all the bottom-states are bipartite product states with constituent states describing fully polarized subsystems containing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_1$$\end{document}N1 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_2$$\end{document}N2 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\le N_1$$\end{document}≤N1) particles, respectively. For two bottom-states different in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_2$$\end{document}N2, the one with a larger \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_2$$\end{document}N2 is higher. For two having the same \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_2$$\end{document}N2, the one with a smaller S is higher. Whereas all the top-states are found to be essentially a product state of the pairs, in each pair the two spins are coupled to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document}λ if the strength of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document}λ-channel is more repulsive than the others. For the states belonging to an S-group, the higher one would contain more pieces. As the energy goes up, larger pieces (those containing more than two particles) will disappear.

The entire spectrum of a bound cold system with eight spin-3 atoms and with interaction leading to spin-parallel pair We first choose a trapped 8-body system with spin-3 cold atoms as the first example for medium-body systems.
Since we aim at the spin-textures, the temperature is assumed to be sufficiently low that all the spatial degrees of freedom are frozen and all particles fall into the same spatial state φ which is most favorable to binding. Note that, when the spin-orbit interaction appears, the lowest spatial state would in general mix up with higher excited spatial states. How strong the mixing would be, depends on the energy gap. Our assumption on the freezing is equivalent to having an infinite gap. Thus, the mixing is impossible. Therefore, the effect of spin-orbit coupling has been completely suppressed and can be neglected. Then the spin-3 boson system is governed by the spindependent Hamiltonian where V ij is the interaction between particles i and j, g is the weighted strength of the -channel ( = 0, 2, 4 , and 6). A factor φ 4 (r)dr , which embodies the effect of spatial profile, has already been included in g . P i,j is the projector to the -channel. Let χ denote a pure single spin-state, (χ(i)χ(j)) ,m denote a two-body spin-state in which the spins of i and j are coupled to with magnetic number m. Then P i,j = m |(χ(i)χ(j)) ,m ��(χ(i)χ(j)) ,m | . We consider the case that H spin is dominated by the interaction that favors the formation of spin-parallel pair, i.e., g 6 is considered to be much more negative than the other three strengths so that the two interacting spins favor being parallel. The other three are first assumed to be equal to each other. Then, the case that they are unequal will be further considered. Note that the spin-textures would not be changed if {g } as a whole is shifted with a common value and/or if the unit of energy is altered. Thus, for the first case, the Hamiltonian can be in general given as g 6 = −1 and g 0 = g 2 = g 4 = 1 , and is denoted as H . By using the Fock-states as basis-states, H can be diagonalized to obtain the eigenenergies and eigenstates (refer to the Supplementary file Appendix). The total spin S and its Z-component M are good quantum numbers. A group of eigenstates having the same S form an S-group. Since magnetic field is not involved, it is sufficient to consider only M = 0 states. The number of M = 0 states included in an S-group is just the multiplicity of S. Let the lowest energy state of an S-group be called a bottom-state (b-state), while the highest a top-state (t-state). The spectrum of H with N = 8 is plotted in Fig. 1, where the excitation energies E X of the b-and t-states of every S-group are marked. Thus, we can see how an S-group is distributed in the spectrum. When a state has a larger S, the spins are more likely to be parallel to each other. Therefore, when a negative g 6 is dominant, the S-group as a whole would be lower when S is larger as shown in Fig. 1. In particular, the ground state (g.s.) would have the largest S = 3N , i.e., a fully polarized core.
Let a fully polarized subsystem be called a piece (p). Let the i-th eigenstate of an S-group ordered in uprising energy be denoted as S,i . We found that S,1 , the b-state of each S-group, is mostly composed of two pieces containing N 1 and N 2 particles respectively ( . We found that, in Zone I marked in Fig. 1 with 24 ≥ S ≥ 18 (where the S = 23 state is prohibited by symmetry), all the six b-states have N 2 = 1 and can be written as χ) S represents the state in which the spin of a 7-body fully polarized subsystem (denoted by [7]) and the spin of a single atom are coupled to S, and η is an N-body spin-state with an extremely small norm as shown in Table 1 (the case with S = 20 is an exception). The associated excitation mode is a split of the fully polarized core into a [N 1 ] piece, N 1 = N − 1 , together with a single particle. Then the latter two, are coupled in various ways leading to the b-states with various S. For a state with a smaller S, the direction of the single spin would deviate more from that of the piece, therefore the probability of forming spin-parallel pairs is reduced. Accordingly, the state is higher in energy. In Zone II with 12 ≤ S ≤ 15 , all the b-states have N 2 = 2 and very close to ([6][2]) S . Thus, a pair of spin-parallel particles (i.e., [2]) have been extracted from the core. In Zone III with 6 ≤ S ≤ 9 , all the b-states are very close to ([5] [3]) S , where a fully polarized triplet [3] has been extracted. In Zone IV with 0 ≤ S ≤ 4 , all the b-states are very close to ([4][4]) S . Thus, a higher excitation energy is needed to extract more particles from the core. The overlaps of the above suggested piece-piece states (p-p states) ([N 1 ][N 2 ]) S and the exact eigen-b-states S,1 are listed in Table 1.
Let j be an index to mark the order of all the eigenstates in up-rising energy. Thus, each state S,i can be equivalently denoted as j . We found that the first excited state j=2 is just the b-state � 22,1 ≃ ([7]χ) 22 , where a single spin has been extracted from the core and S is reduced. Let the mode leading to the transition 1 → 2 be called the cheapest excitation mode (CExM). Then, for a fully polarized core, the CExM is just to extract a single www.nature.com/scientificreports/ spin from the core together with a decrease of S from 3N to 3N − 2 . Since the S = 3N − 1 state is prohibited by symmetry, the decrease in S implies that the single spin is no more lying along its previous direction but deviates from it as least as possible. We found from numerical diagonalization of H that the first excited state always has S = 3N − 2 for all N. Therefore, the idea of CExM can be generalized to any piece. Thus, when a p-p state is further excited, the lowest way is to excite a piece with CExM, namely, Table 1 we found that, while most b-states are a p-p state, a few of them (those lying between two zones, e.g., the group with S = 10 is lying between Zone II and III) are a mixture of p-p and p-1-p as two basis-states. These basis-states are not orthogonal to each other. From the data given in the table together with the overlaps among the basis-states, the combined weight contributed by both basis-states can be known.  17 alone is (0.949) 2 , while by both is 0.992. Thus, we conclude that all the b-states are composed of the p-p and p-1-p basis-states. An exception is the S = 1 state. This state has multiplicity one, therefore its spin-texture is completely determined by symmetry and irrelevant to dynamics.
The second state of each S-group ( S,2 ) contains essentially three species (including χ and [2]), e.g., short) between the suggested p-p or p-1-p state and the exact eigen-b-state S,1 of the Hamiltonian H with N = 8 and f = 3.  21 . Thus, the second CExM of a piece is to extract a particle together with a reduction of S by 3, i.e., In general, when the core splits up into more pieces and the relative orientations of the piece-spins are more deviated from each other, the associated excitation energy would be higher. As an example, we consider the split of [6] via CExM, i.e., [6] 15 . With this mode, the p-p state ([6][2]) S would split further and becomes We found that there are totally seven S,i having |�� 5−1−2,S |� S,i �| > 0.8 as listed in Table 2. These states with the CExM mode are all higher members of an S-group.
When the energy goes higher, larger pieces [N 1 ] with N 1 ≥ 3 will vanish gradually. Instead, more [2] and (2) =6 pairs will emerge. We are not going to the details of every excited states, but concentrate at the product states of the pairs. Let Note that, if g X is most repulsive, the highest energy states would be dominated by the (2) X pairs. For H , there are three strengths equal to each other, i.e., g 0 = g 2 = g 4 . In this case, we found that the highest energy state is Thus the (2) 4 pair is most important. The very high S = 1 state is where the (2) 2 pair emerges. Incidentally, the study of the highest energy states is meaningful because they would become the lowest energy states when the strengths reverse their signs.
Let P ≡ � |P 1,2 | � be the probability of two particles in forming a (2) pair. These probabilities against the index j are plotted in Fig. 2. Figure 2a is for j = 1 → 20 , where P 6 is much larger implying the importance of the fully polarized pieces, while P 0 is much smaller implying that the appearance of (2) 0 pairs is least probable. In average, P 4 > P 2 . The states with j = 7, 9 , and 16 are t-states with S = 20, 19 , and 18 (in Zone I of Fig. 1), respectively. Note that P 4 will arrive at a peak at the t-states. Therefore, for H , the way to maximize the energy under the conservation of S is to increase the (2) 4 pairs. Figure 2b is for the highest 21 states with j = 130 → 151 , where P 6 is remarkably reduced implying the vanish of the polarized pieces, instead, P 4 becomes the largest and is also peaked at the t-states. This fact demonstrates again the domination of the (2) 4 pairs in high energy states (note that the highest energy state j=151 is completely composed of the (2) 4 pairs). In any cases, P 0 remains to be small. Incidentally, the (2) 0 pairs would become dominant in high states when g 0 is most repulsive. This can be shown in the reverse of the spectrum given in Ref. 27 .

The entire spectrum with N = 13
As the next example for medium-body systems, the whole spectrum of H with N = 13 is plotted in Fig. 3.
We      www.nature.com/scientificreports/ is from 0.998 (when S = 8 ) to 1.000. In the above cases η = (1 − c 2 ) 1/2 is so small that the split of the system into two fully polarized pieces is confirmed. The probabilities {P } are plotted in Fig. 4. The t-states are marked with a 't' as before (The j = 1, 2 , and 3 states are also t-states. Their multiplicity is one, therefore they are not marked). We found that Fig. 4 is qualitatively similar to Fig. 2. For example, P 6 ( P 4 ) is more important in Fig. 4a (b), and P 4 is mostly peaked at the t-states.

The spectra of the Hamiltonians deviated from H
To evaluate the effect caused by a deviation from H , we consider three cases (i) g 0 = 0.5 , (ii) g 2 = 0.5 , and (iii) g 4 = 0.5 , while g 6 = −1 remains unchanged and the other two strengths not mentioned in each above case have g = g ′ = 1 . The associated three deviated Hamiltonians are denoted as H (i) , H (ii) , and H (iii) , respectively. Their spectra with N = 8 are given in Fig. 5a, b, and c to be compared with Fig. 1. These figures are qualitatively similar to each other. For example, E X of the b-state varies interruptedly when S decreases and crosses over some critical values, these values are the same for H and {H (z) } ( z = i, ii , and iii). Accordingly, the S-groups can be similarly divided into zones (e.g., the groups with S = 12 to 15 are contained in Zone II in Fig. 1 and also in Fig. 5a, b, and  c). The interruption implies a remarkable difference in spin-texture during the cross-over. For example, in the same zone of the four figures, the values N 2 specifying the p-p states To reveal the similarity among the b-states with these different Hamiltonians, the overlaps �([7]χ) S |� S,i � are given in Table 3 as an example, where all the b-states belong to Zone I with 24 ≥ S ≥ 18 ( S = 23 ). The data demonstrate that all these states are dominated by the same basis-state ([7]χ) S , thus they are nearly unchanged under the change H � → {H (z) }.
However, quantitatively, E X of the b-states given in Figs. 1 and 5 are slightly different. Note that, for H spin , the energy E of a spin-state is Thus, for the b-states in Zone I, the slight difference arises from the probabilities {P } inherent in the basis-state ([7]χ) S , which is completely determined by symmetry as shown in Table 4.
From this table, we know that, among the six ([7]χ) S states different in S, only the one with S = 18 contains the (2) 0 pairs. Therefore, only this state would be benefited from the reduction of g 0 . This explains why only the red circle for 18,1 is shifted to the left in Fig. 5a but not in Fig. 1, 5b, and 5c. In Fig.5b, g 2 is reduced. From the In all these states, due to the fully polarized 7-body core, P 6 is dominant. However, when S decreases more from 3 × N = 24 , the direction of the single spin deviates more from that of the fully polarized core, and therefore leads to a smaller P 6 . www.nature.com/scientificreports/ table we also know that the three states with S = 18, 19 , and 20 contain the (2) 2 pairs, thus they would be benefited, in particular, the S = 19 state would be benefited more. This is shown in Fig. 5b. In Fig. 5c g 4 is reduced.
The energies of all the five states with 18 ≤ S ≤ 22 would be thereby reduced (that of the S = 21 state would be reduced more) as shown in Fig. 5c. Examples for the textures of high-lying b-states are shown in It is impressive that all the overlaps given in Table 5 are also very close to one. Thus, during H � → H (z) , the domination of the ([N − N 2 ][N 2 ]) S textures recovers. The shift of E X occurring in Zone III and IV can be similarly explained as before.
To go beyond the b-states, we add a superscript in j as

The case with spin-4 atoms
For the systems with spin-4 atoms, H has g 8 = −1 and g 0 = g 2 = g 4 = g 6 = 1 . After the numerical diagonalization of H , the qualitative features of the eigenstates are found to be very similar to those with spin-3 atoms. In particular, the domination of the p-p texture ([N 1 ][N 2 ]) S in b-states, the further split of the pieces when energy goes higher, and the emergence of (2) pairs in high-lying states recover. The spectrum with f = 4 and N = 8 is plotted in Fig. 6, which is very similar to Fig. 1, e.g., the distribution against E X of each S-group and the division of the S-groups into zones. Note that in the Zone I of Fig. 1 Table 2.  www.nature.com/scientificreports/ |7f − f| = 18 . This is the same in Fig. 6 where the value is |7f − f| = 24 . Similarly, in Zone II, � S,1 ≃ ([6][2]) S , therefore the least S is |6f − 2f| = 12 in Fig. 1 and 16 in Fig. 6. In Zone III, � S,1 ≃ ([5] [3]) S , therefore the least S is |5f − 3f| = 6 in Fig. 1 and 8 in Fig. 6.

Summary and final remarks
Note that, for the high-f systems with the pairwise interaction leading to spin anti-parallel pair, all the eigenstates are simply a product-state of the (2) 0 pairs together with certain number of unpaired particles 20,21,27 . Whereas in the present case with the interaction leading to spin parallel pair, we found that the g.s. is just a fully polarized core while the excitation mode is a split of the core, accordingly, the excited states are essentially composed of various pieces (fully polarized subsystems). The textures with more pieces and/or with the directions of the spins of the pieces deviated more from each other would lead to a higher energy. All the eigenstates are first classified into groups according to S. All the b-states of the S-groups have the p-p texture ([N 1 ][N 2 ]) S ( N = N 1 + N 2 and N 1 ≥ N 2 ). Thus the S-groups can be classified according to N 2 . For the b-states having the same N 2 , the one with a smaller S is higher. For the b-states with different N 2 , the one with a larger N 2 is higher. When the excitation energy goes up, the pieces may split further, therefore the number of the pieces will increase and smaller pieces will emerge. In particular, those with particle number N p ≥ 3 would gradually vanish, and the eigenstates would be dominated by the pairs. For H with positive g 0 = g 2 = g 4 , the highest energy states are found to be dominated by (2) 4 pairs. Whereas when g 0 = g 2 = g 4 , it is believed that the (2) o pairs would be dominant if g o alone is most repulsive. This can be seen via the reverse of the spectrum given in Ref. 27 where the (2) 0 pairs are dominant.
We have demonstrated numerically that the above picture holds when the Hamiltonian is not seriously deviated from H . Besides, with the same approach of study, we found that the above picture holds also when N ≤ 7 . It implies that some physical findings found in few-body systems could recover in medium-body systems. However, for many-body systems, a different approach of study is needed. How the above picture would be spoiled by increasing N remains to be further clarified. Although only the cases with spin-3 and spin-4 atoms are considered, we believe that, for H , the above picture might hold for arbitrary integer f ≥ 1 . When f is even and g f is most attractive (repulsive), the triplet with three spins coupled to zero might be dominant in the lowlying (high-lying) states. These assumptions deserve also to be further clarified.

Data availability
All data generated or analysed during this study are included in this published article (and its Supplementary Information files). www.nature.com/scientificreports/