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

The spin-textures of bound medium-body systems with spin-[Formula: see text] atoms ([Formula: see text]) have been studied. The Hamiltonian is assumed to be dominated by the two-body interaction favoring parallel spins. The system with particle number [Formula: see text] and [Formula: see text] 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 [Formula: see text] and [Formula: see text] 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 [Formula: see text] and [Formula: see text] ([Formula: see text]) particles, respectively. For two bottom-states different in [Formula: see text], the one with a larger [Formula: see text] is higher. For two having the same [Formula: see text], 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 [Formula: see text] if the strength of the [Formula: see text]-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.

Introduction

It is well known that few-body theories and many-body theories are important for understanding the physical systems with the particle number and the systems with a large N (e.g., condensed matter systems), respectively. Basic structures in few-body systems (e.g., some solid pairs and/or clusters)[1-5] and some basic physical processes (e.g., resonance in 2-body scattering and Efimov effect in diffuse 3-body systems) are known to play an important role in many-body systems[6-9]. For example, the -particles (the nucleus of helium atom) model in nuclear theory[10] and the Cooper-pairs in superconductivity[3]. We believe that our understanding of few-body systems would be in general helpful for understanding the complicated phenomena occurring in many-body systems from a more fundamental point. However, the particle number N involved in few-body systems is greatly different from that in many-body systems. The few-body theories aim at finding exact solutions of Schrödinger equation. Each of them contains a complete set of good quantum numbers, and every physical feature of the system can be extracted from the solution. Whereas each many-body theory is skillfully designed and aims at finding approximate solutions for explaining specific physical phenomena (e.g., for explaining the ground state properties). Not all the good quantum numbers are conserved in these approximate solutions. To check the applicability of a many-body solution, an auxiliary way is to study in detail how the related physical features would vary with the increase of N, namely, to clarify the effect of N. This can be (partially) realized via the study of the medium-body systems (e.g., 5). The knowledge from these systems would play the role as a bridge to relate the few- and many-body physics. Thus, the study of medium-body systems is worthy. The present paper is dedicated to these systems which are less studied previously.

An important progress of technology is the realization of extremely low temperature, e.g., pK[11]. For the Bose-Einstein condensates, when T is sufficiently low, the spatial degrees of freedom can be frozen and, for particles with nonzero spin , only the spin-degrees of freedom are active[12]. Since the spin-dependent force is much weaker than the central force, the spin-textures are very fragile and highly sensitive to a change in environment, e.g., a change in particle number and/or an adjustment in the strengths of interaction. Therefore, the spin degrees of freedom dominate the dynamics of such sensitive systems at . Thus, the understanding of the spin-textures is important. Note that there are many literatures dedicated to the study of spin- systems with [13-19] and [18-21] atoms, but very few with high- systems[22-26]. Since the dimension of the spin-space increases greatly with , the physics involved in high- systems would be rich. The present paper is dedicated to the study of the spin-textures in high- systems () which are less studied previously.

Note that the formation of pairs (a coupling of two particles with their momenta and/or spins reversing to each other) is an important phenomenon in many-body systems, e.g., in superconductivity and in the interacting boson model (IBM) of nuclei[2,3]. In a previous paper[27] the role of singlet pairs in high- condensates has been studied. The singlet pairs arise from the part of 2-body interaction which causes the two interacting spins to be anti-parallel. Thus, we think that it would be interesting to see what happens if the interaction leads to parallel spins.

As a continuation of the study in Ref.[27], the present paper is dedicated to the study of high- medium-body systems with the interaction favoring to parallel spins. and 13, and =3 and 4 are chosen as examples for medium-body high- systems. The strength of interaction of the 2-channel is assumed to be much more attractive than those of other channels. The associated Hamiltonian is exactly diagonalized to obtain the spectra and the whole set of eigenstates. A detailed analysis has been made to find out the similarities and differences existing among the eigenstates of a spectrum and among the spectra of different systems (different in N and/or ). The emphasis is placed on finding out some basic spin-textures and the modes of excitation.

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 spin-dependent Hamiltonianwhere is the interaction between particles i and j, is the weighted strength of the -channel (, and 6). A factor , which embodies the effect of spatial profile, has already been included in . is the projector to the -channel. Let denote a pure single spin-state, denote a two-body spin-state in which the spins of i and j are coupled to with magnetic number m. Then . We consider the case that is dominated by the interaction that favors the formation of spin-parallel pair, i.e., 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 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 and , and is denoted as . By using the Fock-states as basis-states, 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 states. The number of 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 with is plotted in Fig. 1, where the excitation energies 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 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 , i.e., a fully polarized core.

Figure 1

The complete spectrum of an 8-body system with spin-3 trapped cold atoms and with the Hamiltonian . The unit for energy is . Totally there are 151 states (levels). The abscissa is for the excitation energies from 0 (the g.s.) to 45.68 (the highest energy state). The ordinate is for the total spin S of each state. The S-group of states are lying along a short horizontal dot-line. of the b-state and the t-state are marked by a red circle and a blue triangle, respectively. The multiplicity of an S-group is marked at the right of the dot-line if it . When multiplicity, only a red circle is marked for the of the state. For an example, the group with contains 5 states and is distributed from 23.42 to 30.75.

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 . We found that , the b-state of each S-group, is mostly composed of two pieces containing and particles respectively (). Let denote a piece with particles, denote a subsystem where the spins are coupled to (in particular, ). We found that, in Zone I marked in Fig. 1 with (where the state is prohibited by symmetry), all the six b-states have and can be written as , where , 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 is an exception). The associated excitation mode is a split of the fully polarized core into a piece, , 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 , all the b-states have and very close to . Thus, a pair of spin-parallel particles (i.e., [2]) have been extracted from the core. In Zone III with , all the b-states are very close to , where a fully polarized triplet [3] has been extracted. In Zone IV with , all the b-states are very close to . 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) and the exact eigen-b-states are listed in Table 1.

Table 1

The overlaps ( in short) and (( in short) between the suggested p-p or p-1-p state and the exact eigen-b-state of the Hamiltonian with and .

S	7+1	(2+1)+5		S	6+2	(6+1)+1	(2+1)+5
24	1			17	0.647	0.949
22	1			16	0.379	0.975
21	1			15	0.975		0.792
20		0.998		14	0.984
19	0.998			13	0.999
				12	1.000

Let j be an index to mark the order of all the eigenstates in up-rising energy. Thus, each state can be equivalently denoted as . We found that the first excited state is just the b-state , where a single spin has been extracted from the core and S is reduced. Let the mode leading to the transition be called the cheapest excitation mode (CExM). Then, for a fully polarized core, the CExM is just to extract a single spin from the core together with a decrease of S from 3N to . Since the 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 that the first excited state always has 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, , where . With this mode, a p-p state becomes a p-1-p state .

From 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 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. For as an example, the weight contributed by alone is , by alone is , while by both is 0.994. For , the weight contributed by alone is , by alone is , 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 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 () contains essentially three species (including and [2]), e.g.,where is a very small component (, and 0.002, respectively, for the above four cases). In these second states, we found again the further excitation of a piece via the CExM.

Furthermore, from Fig. 1 and Table 1 we found that . Thus, the second CExM of a piece is to extract a particle together with a reduction of S by 3, i.e., , where . 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., . With this mode, the p-p state would split further and becomes . We found that there are totally seven having as listed in Table 2. These states with the CExM mode are all higher members of an S-group.

Table 2

All the eigenstates of with and having the overlap are listed in the table.

S	i	c
10	4	0.926
11	3	0.861
14	5	0.891
17	3	0.954
18	3	0.876
19	2	0.958
20	2	0.828

When the energy goes higher, larger pieces with will vanish gradually. Instead, more [2] and pairs will emerge. We are not going to the details of every excited states, but concentrate at the product states of the pairs. Let . There are totally four eigenstates have the overlap . They are , , , and .

Note that, if is most repulsive, the highest energy states would be dominated by the pairs. For , there are three strengths equal to each other, i.e., . In this case, we found that the highest energy state isThus the pair is most important. The very high state iswhere the 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 be the probability of two particles in forming a pair. These probabilities against the index j are plotted in Fig. 2. Figure 2a is for , where is much larger implying the importance of the fully polarized pieces, while is much smaller implying that the appearance of pairs is least probable. In average, . The states with , and 16 are t-states with , and 18 (in Zone I of Fig. 1), respectively. Note that will arrive at a peak at the t-states. Therefore, for , the way to maximize the energy under the conservation of S is to increase the pairs. Figure 2b is for the highest 21 states with , where is remarkably reduced implying the vanish of the polarized pieces, instead, becomes the largest and is also peaked at the t-states. This fact demonstrates again the domination of the pairs in high energy states (note that the highest energy state is completely composed of the pairs). In any cases, remains to be small. Incidentally, the pairs would become dominant in high states when is most repulsive. This can be shown in the reverse of the spectrum given in Ref.[27].

Figure 2

against j for the case of with and . Probabilities, (a) for the states with , (b) for . The symbols in black square, red circle, blue triangle, and olive reverse-triangle are for , and 6, respectively. All the t-states (the highest one of an S-group) are marked with a ’t’ (except those with multiplicity one, which happens when , and 142).

The entire spectrum with

As the next example for medium-body systems, the whole spectrum of with is plotted in Fig. 3.

Figure 3

The spectrum of with and . Totally there are 920 levels, but only the excitation energies of the b- and t-states of every S-group are marked (refer to Fig. 1).

We found that the spectra of and systems are qualitatively similar. The p-p states recovers in the Zone I of Fig. 3, where , ( is prohibited by symmetry), c is ranged from 0.998 to 1.000 (except that has ). In particular, , i.e., it is excited with the CExM, while , i.e., excited with the second CExM. In Zone II of Fig. 3, the b-states , , and c is ranged from 0.949 (when ) to 1. Thus, two spins have been extracted (i.e., ) as before. When S decreases further, more spins would be extracted. For example, in Zone III with , and with c being ranged from 0.977 (when ) to 1.000. In Zone IV, , , and , where c is from 0.998 (when ) to 1.000. In the above cases is so small that the split of the system into two fully polarized pieces is confirmed.

The probabilities are plotted in Fig. 4. The t-states are marked with a ’t’ as before (The , 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, () is more important in Fig. 4a (b), and is mostly peaked at the t-states.

Figure 4

against j for the case of with and (refer to Fig. 2).

The spectra of the Hamiltonians deviated from

To evaluate the effect caused by a deviation from , we consider three cases (i) , (ii) , and (iii) , while remains unchanged and the other two strengths not mentioned in each above case have . The associated three deviated Hamiltonians are denoted as , , and , respectively. Their spectra with are given in Fig. 5a, b, and c to be compared with Fig. 1. These figures are qualitatively similar to each other. For example, of the b-state varies interruptedly when S decreases and crosses over some critical values, these values are the same for and (, and iii). Accordingly, the S-groups can be similarly divided into zones (e.g., the groups with 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 specifying the p-p states are the same.

Figure 5

The spectra of an 8-body system with spin-3 trapped cold atoms. (a), (b), and (c) are for the Hamiltonian , , and respectively.

To reveal the similarity among the b-states with these different Hamiltonians, the overlaps are given in Table 3 as an example, where all the b-states belong to Zone I with (). The data demonstrate that all these states are dominated by the same basis-state , thus they are nearly unchanged under the change .

Table 3

The overlaps (denoted as 7+1 in short) of with and , where is an eigen-b-state of (, ii, and iii) lying in Zone I.

S	H_‖	H_(i)	H_(ii)	H_(iii)
24	1	1	1	1
22	1	1	1	1
21	1	1	1	1
20	0.955	0.955	0.955	0.955
19	0.998	0.998	0.998	0.998
18	1.000	0.997	0.998	0.996

However, quantitatively, of the b-states given in Figs. 1 and 5 are slightly different. Note that, for , the energy E of a spin-state isThus, for the b-states in Zone I, the slight difference arises from the probabilities inherent in the basis-state , which is completely determined by symmetry as shown in Table 4.

Table 4

The set of the p-p states of with and different in S. In all these states, due to the fully polarized 7-body core, is dominant. However, when S decreases more from , the direction of the single spin deviates more from that of the fully polarized core, and therefore leads to a smaller .

S	P_0	P_2	P_4	P_6
24	0	0	0	1
22	0	0	0.153	0.847
21	0	0	0.224	0.776
20	0	0.079	0.163	0.758
19	0	0.171	0.077	0.752
18	0.083	0.147	0.020	0.750

From this table, we know that, among the six states different in S, only the one with contains the pairs. Therefore, only this state would be benefited from the reduction of . This explains why only the red circle for is shifted to the left in Fig. 5a but not in Fig. 1, 5b, and 5c. In Fig.5b, is reduced. From the table we also know that the three states with , and 20 contain the pairs, thus they would be benefited, in particular, the state would be benefited more. This is shown in Fig. 5b. In Fig. 5c is reduced. The energies of all the five states with would be thereby reduced (that of the state would be reduced more) as shown in Fig. 5c.

Examples for the textures of high-lying b-states are shown in Table 5, where the overlaps and are shown.

Table 5

The overlaps (denoted as ) and (), where are the eigen-b-states of the associated Hamiltonians with and , and they are lying in Zone III and IV. Refer to Table 2.

	H_‖	H_‖	H_(i)	H_(i)	H_(ii)	H_(ii)	H_(iii)	H_(iii)
S	5-3	4-4	5-3	4-4	5-3	4-4	5-3	4-4
9	0.995		0.977		0.984		0.971
8	0.997		0.975		0.991		0.968
7	0.999		0.996		0.995		0.995
6	1.000		0.994		0.994		0.992
4		0.999		0.991		0.994		0.989
2		1.000		0.995		0.994		0.993
0		1.000		0.993		0.993		0.990

It is impressive that all the overlaps given in Table 5 are also very close to one. Thus, during , the domination of the textures recovers. The shift of occurring in Zone III and IV can be similarly explained as before.

To go beyond the b-states, we add a superscript in as or to specify the associated Hamiltonian. Let the overlaps . For a given , among the 151 overlaps different in , the largest one is denoted as . We found that, for the whole set of eigenstates of , 125 of them have , thus each of them has an analogue in the set of eigenstates of . The analogy implies that they are qualitatively unchanged during . In particular, out of the 125, 41 of them have implying a high similarity. When , there are 108 eigenstates having (wherein 33 of them have ). When , there are 111 eigenstates having (wherein 38 of them have ). These data demonstrate that, for most eigenstates, the deviation caused by the deviated Hamiltonians is not large. Thus, the picture based on the pieces and the pairs holds qualitatively when is not deviated from seriously.

The case with spin-4 atoms

For the systems with spin-4 atoms, has and . After the numerical diagonalization of , 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 in b-states, the further split of the pieces when energy goes higher, and the emergence of pairs in high-lying states recover. The spectrum with and is plotted in Fig. 6, which is very similar to Fig. 1, e.g., the distribution against of each S-group and the division of the S-groups into zones. Note that in the Zone I of Fig. 1, , where the least value of S is . This is the same in Fig. 6 where the value is . Similarly, in Zone II, , therefore the least S is in Fig. 1 and 16 in Fig. 6. In Zone III, , therefore the least S is in Fig. 1 and 8 in Fig. 6.

Figure 6

The spectrum of with and . Refer to Fig. 1.

Summary and final remarks

Note that, for the high- systems with the pairwise interaction leading to spin anti-parallel pair, all the eigenstates are simply a product-state of the 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 ( and ). Thus the S-groups can be classified according to . For the b-states having the same , the one with a smaller S is higher. For the b-states with different , the one with a larger 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 would gradually vanish, and the eigenstates would be dominated by the pairs. For with positive , the highest energy states are found to be dominated by pairs. Whereas when , it is believed that the pairs would be dominant if alone is most repulsive. This can be seen via the reverse of the spectrum given in Ref.[27] where the pairs are dominant.

We have demonstrated numerically that the above picture holds when the Hamiltonian is not seriously deviated from . Besides, with the same approach of study, we found that the above picture holds also when . 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 , the above picture might hold for arbitrary integer . When is even and is most attractive (repulsive), the triplet with three spins coupled to zero might be dominant in the low-lying (high-lying) states. These assumptions deserve also to be further clarified.

Supplementary Information.