text:
# Spectra of supersymmetric Yang-Mills quantum mechanics

## 1 Introduction

Since the original conjecture of Banks, Fischler, Shenker and Susskind of the equivalence between M-theory and the D=10 supersymmetric Yang-Mills quantum mechanics (SYMQM) , a lot of effort has been put to understand, and ultimately solve, the latter -.

The general, not necessarily gauge, supersymmetric quantum mechanical systems have much longer history . Claudson and Halpern considered for the first time the gauge systems also well before the BFSS hypothesis . In particular a complete solution of the D=2, N=2 SYMQM was given there (see also Ref. ) . Later the specific gauge models were studied in more detail and some candidates for the ground state were constructed in the, gauge invariant, Born-Oppenheimer approximation .

Another important development was achieved by de Wit, Lüscher and Nicolai who have shown that the spectrum of supersymmetric Yang-Mills quantum mechanics is continuous  due to the cancellations between fermions and bosons. This is different from the pure Yang-Mills case where the transverse fluctuations across the vacuum valleys do not cancel and effectively block the valley, resulting in the discrete spectrum of the 0-volume glueballs . The continuous spectrum was first regarded as a setback, however it has turned out into a virtue with the advent of the BFSS hypothesis and new interpretation in terms of the scattering states. Still this new connection requires existence of the localized state at the threshold of the spectrum - the supergraviton . This question triggered intense studies of the Witten index of SYMQM for various D and different gauge groups . Powerful techniques were developed to calculate analytically nonabelian integrals  related to the Witten index. They were accompanied by complementary and original numerical methods . The emerging picture is rather satisfactory indeed: for $`D<10`$ there is no threshold bound state while for $`D=10`$ there is one, exactly as required by BFSS. This has been proven for N=2, but the evidence is being accumulated that it holds for higher N, as well as for other gauge groups.

The large N limit of SYMQM, which is relevant to M-theory, was studied in the framework of the mean field approximation in Refs . This provided an interesting realization of the black hole thermodynamics predicted by the M-theory .

In spite of all above results, SYMQM remains unsolved for $`D>2`$. Therefore it seems natural to study this model with lattice methods which proved so successful in treating more complex field theoretical systems. Such a programme has been proposed in Ref. , beginning with the yet simpler (quenched, D=4, N=2) case which was later extended to higher gauge groups $`2<N<9`$ . The asymptotic behavior in N was observed, and indications of an onset of the interesting phase structure was found in agreement with Refs.. Including dynamical fermions is possible for D=4, and may be feasible at D=10 for the first few N. However an arbitrary N case is plagued by the sign problem caused by the generally complex fermionic determinant.

Another interesting approach studied by many authors, follows the Eguchi-Kawai trick to trade entirely the D-dimensional configuration space into a group space . Proceeding in this way one is led to consider the fully reduced (to a single point in the Euclidean space) model which nevertheless possesses a version of the superymmetry . Many analytical and numerical results were obtained in this way (for a review see e.g. ) . Numerical simulations were pushed to quite high N in simplified models . However the sign problem which is also present there limits the Monte Carlo approach. This may be alleviated by the new method proposed recently in Ref..

It is important to remember that the subject has a lot of overlap with the small volume study of gauge theories where the valuable expertise has been accumulating for a long time . Although the final goal there is to increase the space volume and to match ultimately the standard large volume physics, the starting point is the pure Yang-Mills quantum mechanics identical to the bosonic part of our supersymmetric systems. In fact the classical results of Lüscher and Münster are the special case of the solution of the D=4 SYMQM in the gluino-free sector of the latter. This will be shown in Sec.5. Recently van Baal has used the full machinery of the small volume approach to study the supersymmetric vacuum state of the more complicated, compact version of D=4 supersymmetric YM quantum mechanics .

Summarizing, even though above quantum mechanics is much simpler that the M-theory, it remains unsolved and thereby still poses an interesting challenge.

In this paper we propose a new approach and present a series of quantitative results for simpler systems not always solved up to now. To this end we use the standard hamiltonian formulation of quantum mechanics in the continuum111This work was inspired by the discussion with C. M. Bender and the methods he developed in studying various quantum mechanical systems ., construct explicitly the (finite) basis of physical states and calculate algebraically matrix representations of all relevant observables. This done, we proceed to calculate numerically the complete spectrum, the energy eigenstates, and identify (super)symmetry multiplets. This approach is entirely insensitive to the sign problem and is equally well applicable to systems with and without fermions. Similarly to the lattice approach, the method has an intrinsic cut-off: any quantitative results can be obtained only within the finite dimensional subspace of the whole Hilbert space. It turns out, however, that in all cases studied below many important characteristics (in particular the low energy spectrum) can be reliably obtained before the number of basis vectors grows out of control.

In the next section we describe the method and its algebraic computer implementation. In sections 3 and 4 the spectra of Wess-Zumino quantum mechanics and D=2 supersymmetric Yang-Mills SU(2) quantum mechanics are derived. Section 5 contains main results of this paper. It is devoted to D=4 supersymmetric Yang-Mills quantum mechanics with the SU(2) gauge group. The global picture of the whole system, with the quantitative spectra in all fermionic channels and their supersymmetric interrelations, will be presented for the first time. Summary and discussion of the future applications follow in the last two Sections.

## 2 Quantum mechanics in a PC

We begin with the simple observation that the action of any quantum mechanical observable can be efficiently implemented (e.g. in an algebraic program) if we use the discrete eigen basis

$$\{|n>\},|n>=\frac{1}{\sqrt{n!}}(a^{})^n|0>,$$
(1)
of the occupation number operator $`a^{}a`$. For example, the bosonic coordinate and momentum operators can be written as

$$x=\frac{1}{\sqrt{2}}(a+a^{}),p=\frac{1}{i\sqrt{2}}(aa^{}),$$
(2)
where all dimensionfull parameters are set to 1. Since typical quantum observables are relatively simple functions of $`x`$ and $`p`$, they can be represented as the multiple actions of the basic creation and annihilation operators222The method can be also extended to non polynomial potentials.. Including fermionic observables is straightforward and will be done in subsequent Sections. Generalization to more degrees of freedom is also evident and will be carried out separately for individual systems.

The first step to quantify above considerations is to implement the Hilbert space in an algebraic program. Any quantum state is a superposition of arbitrary number, $`n_s`$, of elementary states $`|n>`$

$$|st>=\mathrm{\Sigma }_I^{n_s}a_I|n^{(I)}>,$$
(3)
and will be represented as a Mathematica list 333Of course any other algebraic language can be used. In more advanced and time consuming applications we shall be using much faster, compiler based, languages anyway.

$$st=\{n_s,\{a_1,\mathrm{},a_{n_s}\},\{n^{(1)}\},\{n^{(2)}\},\mathrm{},\{n^{(n_s)}\}\},$$
(4)
with $`n_s+2`$ elements. The first element specifies the number of elementary states entering the linear combination, Eq.(3), the second is the sublist supplying all (in general complex) coefficients $`a_I,I=1,\mathrm{},n_s`$, and the remaining $`n_s`$ sublists specify the occupation numbers of elementary, basis states. In particular, convention (4) implies for an elementary state

$$|n>\{1,\{1\},\{n\}\}.$$
(5)
In the case of more degrees of freedom the single occupation numbers $`n^{(I)}`$ in Eq.(4) will become lists themselves containing all occupation numbers corresponding to the individual degrees of freedom.

In the next step we implement basic operations defined in the Hilbert space: addition of two states, multiplication by a number and the scalar product. All these can be simply programmed as definite operations on Mathematica lists. Table 1 displays corresponding routines which transform lists in accord with the principles of quantum mechanics. This done, it is a simple matter to define the creation and annihilation operators which act as a list-valued functions on above lists. This also defines the action of the position and momentum operators according to Eq.(2). Then we proceed to define any quantum observables of interest: hamiltonian, angular momentum, generators of gauge transformations and supersymmetry generators, to name only the most important examples.

Now our strategy should be obvious: given a particular system, we define the list corresponding to the empty state, then we generate a finite basis of $`N_{cut}`$ vectors and calculate matrix representations of the hamiltonian and other quantum operators using above rules. Thereby the problem is reduced to a simple question in linear algebra, namely the spectrum of the system is given by the eigenvalues and eigenvectors of the hamiltonian matrix, their behaviour under rotations by the angular momentum matrix, etc.

Of course the most important question is how much the ultimate physics at $`N_{cut}=\mathrm{}`$ is distorted by the finite cut-off. This can be answered quantitatively by inspecting the dependence of our results on $`N_{cut}`$ for practically available sizes of the bases. In all systems studied so far (and discussed below) the answer is positive: one can extract the meaningful (i.e. $`N_{cut}=\mathrm{}`$) results before the size of the basis becomes unmanageable. However this can be answered only a posteriori and individually for every system.

The method outlined above turns out to be a rather powerful tool capable to solve quantitatively various quantum mechanical problems with finite but large number of degrees of freedom.

## 3 Wess-Zumino quantum mechanics

This system is obtained by dimensional reduction of the four-dimensional Wess-Zumino model leaving only the time dependence of the dynamical degrees of freedom. Resulting quantum mechanics (WZQM) has one complex bosonic variable $`\varphi (t)=x(t)+iy(t)x_1(t)+ix_2(t)`$ and two complex Grassmann-valued fermions $`\psi _\alpha (t),\alpha =1,2`$ . The hamiltonian reads

$$H=\frac{1}{2}(p_x^2+p_y^2)+|F|^2+((m+2g\varphi )\psi _1\psi _2+h.c.),$$
(6)
where the auxiliary field $`F`$ is determined in terms of $`\varphi `$ by the classical equation of motion,

$$F=(m\varphi +g\varphi ^2),$$
(7)
since it enters only quadratically into the action. Independent variables satisfy the standard commutation relations

$$[x_i,p_k]=i\delta _{ik},\{\psi _\alpha ,\psi _\beta ^{}\}=\delta _{\alpha \beta },$$
(8)
which allow us to parameterize these coordinates in terms of the creation and annihilation operators

$`x={\displaystyle \frac{1}{\sqrt{2}}}(a_x+a_x^{}),p_x={\displaystyle \frac{1}{i\sqrt{2}}}(a_xa_x^{}),`$ (9)
$`y={\displaystyle \frac{1}{\sqrt{2}}}(a_y+a_y^{}),p_y={\displaystyle \frac{1}{i\sqrt{2}}}(a_ya_y^{}),`$ (10)
$`\psi _\alpha =f_\alpha ,\psi _\alpha ^{}=f_\alpha ^{},`$ (11)
which in turn satisfy

$$[a_i,a_k^{}]=\delta _{ik},\{f_\alpha ,f_\beta ^{}\}=\delta _{\alpha \beta }.$$
(12)
Now we implement this quantum mechanical system in the computer as explained in the previous Section. The empty state will be represented by a list

$$|(0,0),(0,0)>\{1,\{1\},\{\{0,0\},\{0,0\}\},$$
(13)
where the first brackets specify occupation numbers of the bosonic, and the second fermionic, states. To define properly the fermionic creation/annihilation operators we follow the original construction of Jordan and Wigner

$$f_\alpha =\mathrm{\Pi }_{\beta <\alpha }(1)^{F_\beta }\sigma _\alpha ^{}f_\alpha ^{}=\mathrm{\Pi }_{\beta <\alpha }(1)^{F_\beta }\sigma _\alpha ^+,$$
(14)
where the fermionic number operator $`F_\alpha =f_\alpha ^{}f_\alpha `$ (no sum) and $`\sigma _\alpha ^\pm `$ are rising and lowering operators, commuting for different $`\alpha `$, with $`(\sigma _\alpha ^\pm )^2=0`$ . The Jordan-Wigner phases ensure uniform anticommutation rules for fermionic operators.

We proceed to construct the finite eigen basis of the fermionic, $`F_\alpha `$, and bosonic, $`B_i`$, number operators. For more than one degree of freedom the organization of a basis and definition of the cut-off, $`N_{cut}`$, is not unique. In this case we have chosen the cut-off to be the maximal number of bosonic quanta. That is, the basis subject to the cut-off $`N_{cut}`$ consists of all orthonormal elementary states with $`B=B_1+B_2N_{cut}`$ and all allowed fermionic quanta, i.e. $`F=F_1+F_22`$. In practice the basis is created by action on the empty state with all elementary independent monomials of bosonic creation operators up to $`N_{cut}`$’th order, followed by the three independent monomials of fermionic creation operators. Hence the Hilbert space, subject to the cut-off $`N_{cut}`$, has $`2(N_{cut}+1)(N_{cut}+2)`$ dimensions.

Given the basis, the matrix representation of the hamiltonian (and any other observable of interest) can be readily calculated with our computer based ”quantum algebra”. The spectrum and energy eigenstates are then obtained by numerical diagonalization of the hamiltonian matrix.

Figure 1 shows the spectrum of low energy states as a function of the cut-off up to $`N_{cut}`$=15 for $`m=g=1`$. Energies of the fourth and fifth level are shifted by 1 to avoid confusion of levels for low $`N_{cut}`$. Clear convergence with $`N_{cut}`$ is seen, well within a capacity of a medium class PC. The convergence is faster for lower states. This general feature of our method can be easily understood. The basis generated by creation operators, Eq.(12), is nothing but the eigen basis of some normalized harmonic oscillator. It is obvious that the ground state can be easier approximated by a series of harmonic oscillator wave functions that the higher states with all their detailed structures, zeroes, etc. Quantitatively, at $`N_{cut}`$=10 the exact supersymmetric ground state energy is reproduced with the 3% accuracy which further improves to below 1% at $`N_{cut}`$=15 444Since the exact value is zero, we take the first excited energy as a reference scale.. Supersymmetric pattern of the spectrum is also evident. First, the lowest bosonic state has no fermionic counterpart, its energy tends to zero, hence we identify it with the supersymmetric vacuum of the model 555The model exhibits additional two-fold degeneracy in bosonic and fermionic sectors which is not related to supersymmetry . Second, all higher states group into bosonic-fermionic multiplets with the same energy. Splittings inside the multiplets decrease with $`N_{cut}`$ as shown in Table 2. We conclude that supersymmetry at low energies is restored at the level of few percent in the cut Hilbert space with states containing up to 15 bosonic quanta.

Let us see how the above convergence to the supersymmetric spectrum shows up in the Witten index, which in this approach can be calculated directly from the definition

$$I_W(T)=\mathrm{\Sigma }_i(1)^{F_i}\mathrm{exp}(TE_i).$$
(15)
Fig.2 shows Witten index for $`4<N_{cut}<15`$. Nice convergence to the exact, time independent, result $`I_W(T)=2`$ is observed. Two things are happening which ensure this behaviour. First, there are two (due to the above degeneracy within bosonic and fermionic sectors) supersymmetric ground states. At the largest cut-off, $`N_{cut}`$=15, their energies are not exactly zero but tiny $`E_{0,0}=0.0046,E_{0,1}=0.0071`$, hence giving rise to the exponential fall-off at large $`T`$ with a very small slope. Second, cancellations between non zero energy states within supersymmetric multiplets are not exact and still leave exponential terms, albeit with smaller and smaller (with increasing $`N_{cut}`$) coefficients. This gives the residual time dependence which is however much weaker than the exponential one.

Another, rather general feature of $`I_W(T)`$, can be also observed here. Exactly at $`T=0`$ Witten index vanishes for every $`N_{cut}`$ just because the bases generated with our procedure automatically satisfy a ”global” supersymmetry requirement. Namely the number of bosonic and fermionic states in a basis is the same for every $`N_{cut}`$. Therefore we conclude that

$$\underset{N_{cut}\mathrm{}}{lim}I_W(0,N_{cut})=0.$$
(16)
On the other hand exact supersymmetry requires that this complete correspondence between fermionic and bosonic states is violated since the ground state at $`E=0`$ does not have its fermionic counterpart while each non zero energy state has. This apparent paradox is reconciled by noting that exact supersymmetry emerges only in the limit of the infinite $`N_{cut}`$. In another words the unbalanced fermionic state is pushed (with increasing $`N_{cut}`$) to higher and higher energies and eventually ”vanishes from the spectrum”. However it leaves a visible effect - Witten index has a discontinuity at $`T=0`$ which can be expressed as the noncommutativity of the $`T0`$ and $`N_{cut}\mathrm{}`$ limits

$$\underset{N_{cut}\mathrm{}}{lim}\underset{T0}{lim}I_W(T,N_{cut})\underset{T0}{lim}\underset{N_{cut}\mathrm{}}{lim}I_W(T,N_{cut}).$$
(17)
The right hand side of the above equation, known as the bulk contribution, may not be integer for the continuous spectrum.

Above results were obtained for $`m=g=1`$, but the method works equally well for any other choices of parameters. The case $`m=0`$ has a rotational symmetry and can be solved by other methods . To our knowledge no quantitative calculation of the spectra and wave functions for $`m0`$ exists. As mentioned earlier this approach provides a complete quantum mechanics of a system. In particular, wave functions in the coordinate representation are simply given by the well defined linear combinations of the harmonic oscillator wave functions. More detailed discussion of the model, and comparison with other authors will be presented elsewhere.

## 4 D=2, supersymmetric, SU(2) Yang-Mills quantum mechanics

This system, although solved analytically, has two new features: gauge invariance and continuous spectrum. Therefore we have chosen it as a next exercise. Reducing from $`D=2`$ to one time dimension one is left with the three real, colored, bosonic variables $`x_a(t)`$ and three complex, fermionic degrees of freedom $`\psi _a(t)`$ also in the adjoint representation of SU(2) , $`a=1,2,3`$.

The hamiltonian reads

$$H=\frac{1}{2}p_ap_a+igϵ_{abc}\psi _a^{}x_b\psi _c,$$
(18)
where the quantum operators $`x,p,\psi ,\psi ^{}`$ satisfy

$$[x_a,p_b]=i\delta _{ab},\{\psi _a,\psi _b^{}\}=\delta _{ab}.$$
(19)
Hence they can be written in terms of the creation and annihilation operators as before

$`x_a={\displaystyle \frac{1}{\sqrt{2}}}(a_a+a_a^{}),`$ $`p_a={\displaystyle \frac{1}{i\sqrt{2}}}(a_aa_a^{}),`$ (20)
$`\psi _a=f_a,`$ $`\psi _a^{}=f_a^{}.`$ (21)
To implement fermionic creation and annihilation operators we again used the Jordan-Wigner construction.

The system has a local gauge invariance with the generators

$$G_a=ϵ_{abc}(x_bp_ci\psi _b^{}\psi _c).$$
(22)
Therefore the physical Hilbert space consists only of the gauge invariant states. This can be easily accommodated in our scheme noting that the gauge generators of the SU(2) are just the angular momentum operators acting in color space. In fact the fermionic part of Eq.(22) can be also interpreted in this way since the momentum canonically conjugate to $`\psi `$: $`\pi _\psi =i\psi ^{}`$. Therefore we construct all possible invariant under SU(2) combinations of the creation operators (referred for short as creators) and use them to generate a complete gauge invariant basis of states. There are four lower order creators:

$`a_a^{}a_a^{}(aa),`$ (23)
$`a_a^{}f_a^{}(af),`$
$`ϵ_{abc}a_a^{}f_b^{}f_c^{}(aff),`$
$`ϵ_{abc}f_a^{}f_b^{}f_c^{}(fff).`$ (24)
Pauli principle implies that $`(ff)=(af)^2=(aff)^2=(fff)^2=0`$, therefore the whole basis can be conveniently organized into the four towers of states, each tower beginning with one of the following states

$$|0_F>=|0>;|1_F>=(af)|0>,|2_F>=(aff)|0>,|3_F>=(fff)|0>,$$
(25)
where we have labeled the states by the gauge invariant fermionic number $`F=f_a^{}f_a`$. The empty state (and its Mathematica representation) reads

$$|(0,0,0),(0,0,0)>\{1,\{1\},\{\{0,0,0\},\{0,0,0\}\}\},$$
(26)
with the obvious assignment of bosonic and fermionic occupation numbers. To obtain the whole basis it is now sufficient to repeatedly act on the four vectors, Eq.(25), with the bosonic creator (aa), since action with other creators either gives zero, due to the Pauli blocking, or produces linearly dependent state from another tower. Basis with the cut-off $`N_{cut}`$ is then obtained by applying $`(aa)`$ up to $`N_{cut}`$ times to each of the four ”ground” states. Of course our cut-off is gauge invariant, since it is defined in terms of the gauge invariant creators. Moreover, since the gauge generators can be implemented into the algebraic program as any other observables, we can directly verify gauge invariance of the basis and any other state in question.

Once the above basis is constructed, we proceed to calculate the matrix representation of the hamiltonian, the spectrum and the eigenstates according to Sec.2 666All numerical results are for $`g=1`$.. Again all computation can be done on a small PC and the longest run, corresponding to 20 bosonic quanta, takes approximately 500 sec. It is worth to mention that the most time consuming are the algebraic operations on the states (c.f. Table 1) in the PC-based abstract Hilbert space. Consequently the program spends most of the time calculating the matrix representations of various operators 777In fact this part can be (and often is) done analytically in the integer arithmetic. while the numerical diagonalization is rather fast.

The hamiltonian, Eq.(18) can be rewritten

$$H=\frac{1}{2}p_ap_a+gx_aG_a,$$
(27)
hence it reduces to that of a free particle in the physical, gauge invariant basis. It follows that it preserves the gauge invariant fermionic number and, consequently, can be diagonalized independently in each sector spanned by the four towers in Eq.(25).

The spectrum is doubly degenerate because of the particle-hole symmetry which relates empty and filled fermionic states ($`|0_F>|3_F>`$) and their 1-particle 1-hole counterparts ($`|1_F>|2_F>`$). This symmetry is not violated by the cut-off and indeed we see it exactly in the spectrum. On the other hand, supersymmetry connects sectors which differ by 1 in the fermionic number e.g. ( $`|0_F>|1_F>`$ ) and in general is restored only at infinite $`N_{cut}`$ . A good measure of the SUSY violation is provided by the energy of the ground state which should be 0. We see that the energy of the lowest (doubly degenerate) state converges to 0 but rather slowly, c.f. Fig.3. This is interpreted as the indication that the spectrum is continuous at infinite cut-off. Indeed it is hard to approximate the non localized state by the harmonic oscillator states (this in fact is our basis) and consequently the convergence of such an approximation must be slow 888This argument can be turned into a proof that the free spectrum converges like 1/$`N_{cut}`$ which is also confirmed here (c.f. Fig.3).. This deficiency can be sometimes turned into an advantage as will be seen in the D=4 case.

Surprisingly however, there exists a particular scheme of increasing the basis which renders exact supersymmetry at every finite cut-off $`N_{cut}`$ in this model. To see that, let us first discuss tests of the SUSY on the operator level. SUSY generators read

$$Q=\psi _ap_a,\overline{Q}_a=\overline{\psi }_ap_a,$$
(28)
with $`\overline{\psi }=\psi ^{}`$, and satisfy the algebra

$$\{Q,\overline{Q}\}=2Hgx_aG_a,Q^2=\overline{Q}^2=0.$$
(29)
Since the matrix elements of the SUSY generators can be easily calculated in our PC-based Hilbert space approach, one can readily check the matrix element version of Eqs.(29). Indeed we find that these relations are satisfied to better and better precision with increasing $`N_{cut}`$. If so, it is natural to ask what is the spectrum of the hamiltonian matrix $`H^{Q\overline{Q}}`$ defined by Eq.(29) at finite cut-off. To this end we diagonalize the matrix

$$H_{N,N^{}}^{Q\overline{Q}}=\frac{1}{2}(<N|Q|M><M|\overline{Q}|N^{}>+<N|\overline{Q}|M><M|Q|N^{}>),$$
(30)
where $`N,N^{}`$ and $`M`$ label vectors of our finite basis999The second term in Eq.(29) does not contribute in the physical basis.. Of course SUSY generators mix the fermionic number, therefore one should combine all four towers of Eq.(25) into one big basis. One more trick is required to achieve exact SUSY at finite cut-off: we have to increase the basis allowing from the beginning for the disparity between the fermionic and bosonic sectors. That is, the size of the basis grows as: $`2+4+4+\mathrm{}`$. Hence dimension of the cut Hilbert space is $`2+4N_{cut}`$ in this scheme. With this choice the spectrum of the $`H^{Q\overline{Q}}`$ has exact superymmetry as shown in Table 3.

As required, the ground state has zero energy and does not have the supersymmetric image while higher states form supersymmetric doublets with the same energy. Additional degeneracy is caused by the particle-hole symmetry as explained earlier. All non zero eigenvalues (with finite ”principal” quantum number) tend to zero with increasing $`N_{cut}`$ and form, at $`N_{cut}=\mathrm{}`$, the continuum spectrum of a free hamiltonian, Eq.(27). To reach the solutions with non zero energy, one should appropriately scale the index of an eigenvalue with $`N_{cut}`$ . The four towers of eigenstates are in the direct correspondence with the four gauge invariant plane waves constructed in Ref..

To our knowledge this is the first case where the cut-off preserves exactly the supersymmetry. On the other hand this success of the particular cut-off scheme may be related to the simple structure of the D=2 SYM quantum mechanics and its complete solubility. Nevertheless, it is important to keep in mind that at finite $`N_{cut}`$ one is free to define the hamiltonian within the ”O(1/$`N_{cut}`$)” tolerance, i.e. as long as various definitions converge to the same limit at infinite $`N_{cut}`$. We decided to use $`H^{Q\overline{Q}}`$ because it is the anticommutator of charges which is used to derive basic properties of the SUSY spectrum. Further study of this exact realization of supersymmetry will be continued elsewhere.

Finally let us shortly discuss the Witten index for this model. Because the particle-hole symmetry interchanges odd and even fermionic numbers, the Witten index vanishes identically. This is also true for any finite cut-off since $`N_{cut}`$ preserves particle-hole symmetry. Nevertheless one can obtain a non trivial and interesting information by defining the index restricted to the one pair of fermion-boson sectors. For example

$$I_W(T)^{(0,1)}=\mathrm{\Sigma }_{i,F_i=0,1}(1)^{F_i}\mathrm{exp}(TE_i),$$
(31)
where the sum is now restricted only to the $`F=0`$ and $`F=1`$ sectors. Since supersymmetry balances fermionic and bosonic states between these sectors (with the usual exception of the vacuum), the restricted index is a good and non trivial measure of the amount of the violation/restoration of SUSY, even when the total Witten index vanishes. Obviously $`I_W^{(2,3)}=I_W^{(0,1)}`$, and $`I_W^{(0,3)}=I_W^{(1,2)}=0`$ . Studying restricted index is particularly interesting in this model since, due to the continuum spectrum, it does not have to be integer and provides some information about the density of states. We have calculated the index from our spectrum of the original hamiltonian, Eq.(18), for a range of cut-offs: $`5<N_{cut}<20`$. Figure 4 shows the result for the $`2+4+4+\mathrm{}`$ scheme, i.e. when the basis is increased by every four vectors allowing from the beginning for the two unbalanced states from the empty and filled sectors each. We see a slow approach towards the time independent constant which seems to be $`1/2`$.

Witten index, as a sum over all states, provides an average measure of the breaking/restoration of SUSY. Indeed even though individual supersymmetric multiplets have not yet clearly formed (at cut-offs used for this calculation), we see definite flattening of the $`T`$ dependence of $`I_W(T)`$ which must then be a result of the average cancellation between many levels. Therefore it is easier to see the onset of the supersymmetric behavior in the Witten index than in the location of the individual levels. As before Fig. 4 shows that the limiting value of the index is discontinuous at $`T=0`$ with the same interpretation as in Sect.3. However, contrary to the discrete spectrum of the Wess-Zumino model, the value of the index is not integer 101010It would be an interesting exercise to calculate the restricted index from the continuous spectrum of Ref... The value of the index at $`T=0`$ is the direct consequence of our scheme ($`2+4+4+\mathrm{}`$) of increasing the basis and, of course, is scheme dependent. As a final remark we note an intriguing existence of a ”universal” point in $`T`$ at which the asymptotic value seems to be attained at all $`N_{cut}`$.

## 5 D=4 supersymmetric SU(2) Yang-Mills quantum mechanics

### 5.1 Preliminaries

Dimensionally reduced (in space) system is described by nine bosonic coordinates $`x_a^i(t)`$, $`i=1,2,3;a=1,2,3`$ and six independent fermionic coordinates contained in the Majorana spinor $`\psi _a^\alpha (t)`$, $`\alpha =1,\mathrm{},4`$. Equivalently (in D=4) one could have chosen to impose the Weyl condition and worked with Weyl spinors. Hamiltonian reads

$`H=H_B+H_F,`$ (32)
$`H_B={\displaystyle \frac{1}{2}}p_a^ip_a^i+{\displaystyle \frac{g^2}{4}}ϵ_{abc}ϵ_{ade}x_b^ix_c^jx_d^ix_e^j,`$
$`H_F={\displaystyle \frac{ig}{2}}ϵ_{abc}\psi _a^T\mathrm{\Gamma }^k\psi _bx_c^k,`$ (33)
where $`\psi ^T`$ is the transpose of the real Majorana spinor, and $`\mathrm{\Gamma }`$ in D=4 are just the standard Dirac $`\alpha `$ matrices. In all explicit calculations we are using Majorana representation of Ref..

Even though the three dimensional space was reduced to a single point, the system still has the internal spin(3) rotational symmetry, inherited from the original theory, and generated by the angular momentum

$$J^i=ϵ^{ijk}\left(x_a^jp_a^k\frac{1}{4}\psi _a^T\mathrm{\Sigma }^{jk}\psi _a\right),$$
(34)
with

$$\mathrm{\Sigma }^{jk}=\frac{i}{4}[\mathrm{\Gamma }^j,\mathrm{\Gamma }^k].$$
(35)
Further, the system has gauge invariance with the generators

$$G_a=ϵ_{abc}\left(x_b^kp_c^k\frac{i}{2}\psi _b^T\psi _c\right),$$
(36)
and is invariant under the supersymmetry transformations generated by

$$Q_\alpha =(\mathrm{\Gamma }^k\psi _a)_\alpha p_a^k+igϵ_{abc}(\mathrm{\Sigma }^{jk}\psi _a)_\alpha x_b^jx_c^k.$$
(37)
The bosonic potential (written now in the vector notation in the color space)

$$V=\frac{g^2}{4}\mathrm{\Sigma }_{jk}(\stackrel{}{x}^j\times \stackrel{}{x}^k)^2,$$
(38)
exhibits the famous flat directions responsible for a rich structure of the spectrum. At first sight one might expect that the spectrum of the purely bosonic (hence non supersymmetric) system does not have localized states because of these flat valleys (c.f. Introduction). However the flat directions are blocked by the energy of the transverse quantum fluctuations since valleys narrow as we move from the origin. As a consequence the spectrum of the model is discrete. Energies of the first lower states, known as zero volume glueballs, were first calculated by Lüscher and Münster . On the other hand, in the supersymmetric system, transverse fluctuations cancel among bosons and fermions, valleys are not blocked, and the spectrum of the model is expected to be continuous . The story has its interesting continuation in D=10 model where the evidence of the threshold bound state was accumulated . Existence of such a state is necessary for the M-theory interpretation of the model where it is considered as a prototype of the graviton multiplet. Therefore detailed study of the low energy spectra of the whole family of Yang-Mills quantum mechanical systems, including identification of the localized and non localized states is an important and fascinating subject. We shall face some of these questions also in the D=4 system.

### 5.2 Creation and annihilation operators

There are many ways to write quantum coordinates in terms of creation and annihilation operators

$$[a_a^i,a_b^k]=\delta _{ab}^{ik},\{f_a^\rho ,f_b^\sigma \}=\delta _{ik}^{\rho \sigma },\rho ,\sigma =1,2.$$
(39)
The only constraint comes from the canonical (anti)commutation relations

$$[x_a^i,p_b^k]=i\delta ^{ik}\delta _{ab},\{\psi _a^\alpha ,\psi _b^\beta \}=\delta _{ab}^{\alpha \beta }.$$
(40)
For bosonic variables we just use the standard extension of Eq.(21) to more degrees of freedom

$$x_a^i=\frac{1}{\sqrt{2}}(a_a^i+a_a^i),p_a^i=\frac{1}{i\sqrt{2}}(a_a^ia_a^i).$$
(41)
For fermionic variables we begin with the classical Majorana fermion in the Weyl representation

$$\psi _W=(\zeta _2^{},\zeta _1^{},\zeta _1,\zeta _2),$$
(42)
replace the classical Grassmann variables $`\zeta ,\zeta ^{}`$ by fermionic creation and annihilation operators $`f,f^{}`$, and transform to the Majorana representation. Final result is a quantum hermitean Majorana spinor

$$\psi _a=\frac{1+i}{2\sqrt{2}}\left(\begin{array}{c}f_a^1if_a^2+if_a^1+f_a^2\\ +if_a^1f_a^2f_a^1+if_a^2\\ f_a^1+if_a^2+if_a^1f_a^2\\ if_a^1f_a^2+f_a^1+if_a^2\end{array}\right),$$
(43)
which satisfies Eq.(40) due to (39). Other choices of fermionic creation and annihilation operators are also possible .

The next step is to define an empty state and its, computer based, algebraic representation

$`|(0,0,0),(0,0,0),(0,0,0),(0,0,0),(0,0,0)>`$ (45)
$`\{1,\{1\},\{\{0,0,0\},\{0,0,0\},\{0,0,0\},\{0,0,0\},\{0,0,0\}\}\},`$
where the first three vectors (in color) specify bosonic, and the last two fermionic, occupation numbers.

The Jordan-Wigner transformation requires additional specification in this case. As before, the action of $`f_a^i`$ and $`f_a^i`$ on elementary states is defined by the spin-like raising and lowering operators corrected by the non local Jordan-Wigner phases

$$f_A=\mathrm{\Pi }_{B<A}(1)^{F_B}\sigma _B^{}f_A^{}=\mathrm{\Pi }_{B<A}(1)^{F_B}\sigma _A^+.$$
(46)
However, since individual states are now labeled by a double index $`A`$, one must define the ordering of the two dimensional indices. We choose the lexicographic order. If $`A=(a,\alpha )`$ and $`B=(b,\beta )`$ than

$$\mathrm{\Pi }_{B<A}=\mathrm{\Pi }_{\beta <\alpha ,b3}\mathrm{\Pi }_{\beta =\alpha ,b<a}.$$
(47)
Any other unambiguous definition of the ordering is admissible.

### 5.3 Constructing the basis

There are many ways to construct bases in higher dimensional systems. Apart from theoretical requirements one must also take into account practical limitations of computer implementation. Of course the number of states at finite cut-off is generally bigger for higher D. However, the theory also has more symmetries and consequently the number of states in a particular channel can be kept manageable.

First important symmetry is the conservation of the number of Majorana fermions

$$[F,H]=0F=f_a^if_a^i.$$
(48)
That is, the vector interaction $`H_F`$ cannot produce Majorana pairs in this model. As a consequence the whole Hilbert space splits into seven sectors of fixed $`F=0,1,\mathrm{},6`$.

Second, the system has the particle-hole symmetry

$$F6F,$$
(49)
therefore the first four sectors $`F=0,1,2,3`$ contain all information about the spectrum.

Further, the gauge invariance requires that our basis is built only from the gauge invariant creators. All this is not different from the earlier D=2 case.

The new element is the SO(3) invariance, with the angular momentum Eq.(34), which can be used to split further the problem into the different channels of fixed $`J`$. However, more we specify the basis more complicated its vectors become, and consequently more computer time goes into the generation of such a basis and subsequent calculation of matrix elements. Similarly one might contemplate using the reduced version of Lorentz invariant composite operators as possible creators of a basis. Again this would generate more complex basis since fields contain both creation and annihilation operators. Also, Lorentz covariance is not that relevant in our fixed frame, Hamiltonian formulation.

Taking all above into account we have decided to produce the simplest basis of gauge invariant vectors using elementary creation operators. To this end we proceed as follows. Consider each fermionic sector separately, i.e. fix the fermionic number $`0F3`$. At given $`F`$ create all independent vectors with fixed number of bosonic quanta $`B=a_a^ia_a^i`$, and define a cut-off as the maximal number of bosonic quanta $`BN_{cut}`$, independently for each $`F`$. Hence $`N_{cut}`$ can depend on $`F`$. To create all independent, gauge invariant states at fixed $`F`$ and $`B`$ consider all possible contractions of color indices in a creator of $`(F,B)`$ order

$$a_{a_1}^{i_1}\mathrm{}a_{a_B}^{i_B}f_{b_1}^{\sigma _1}..f_{b_F}^{\sigma _F},$$
(50)
for all values of the spatial indices $`i`$ and $`\sigma `$. All color contractions fall naturally into different gauge invariant classes. Creators from different classes differ by color contractions between bosonic and fermionic operators. For example

$$a_a^ia_a^ja_b^ka_b^lf_c^\sigma f_c^\rho $$
(51)
and

$$a_a^ia_a^ja_b^ka_c^lf_c^\sigma f_b^\rho $$
(52)
belong to different gauge invariant classes according to our definition. Another example involves odd number of operators where one ”contraction” consists of one triple of color indices coupled with the Levi-Civita symbol 111111For the SU(2) gauge group, this step can be generalized to higher SU(N) groups.. Namely

$$ϵ_{cde}a_a^ia_a^ja_b^ka_b^la_c^mf_d^\sigma f_e^\rho $$
(53)
and

$$ϵ_{cde}a_a^ia_a^ja_b^ka_c^la_d^mf_e^\sigma f_b^\rho $$
(54)
also belong to different gauge invariant classes. Different contractions within bosonic and fermionic family of operators do not lead to independent states hence are considered to be in the same class. Still there will be many linearly dependent states generated by above creators with all values of spatial indices $`i`$ and $`\sigma `$. It is however a simple matter to recognize linearly independent ones given our rules of ”quantum algebra”. We therefore leave it to the computer for the time being. At fixed $`F`$ and $`B`$ the final procedure is as follows: 1) identify all gauge invariant classes of creators, 2) loop over all values of spatial indices and for each $`i_1,\mathrm{},i_B,\sigma _1,\mathrm{},\sigma _F`$ create corresponding state from the empty state, Eq.(45), 3) identify and reject linearly dependent vectors, 4) orthonormalize the remaining set of states. In principle this procedure depends exponentially on $`N_{cut}`$, hence can be further shortened and improved, however it is sufficiently fast for our present purposes. It also has an advantage of generating all channels of angular momentum (available within the cut-off $`N_{cut}`$). This will prove useful in studying supersymmetry.

Table IV shows the sizes of bases we have reached so far in each of the fermionic sectors. Due to the particle-hole symmetry the structure in the $`F=4,5,6`$ sectors is identical with that in $`F=2,1,0`$ respectively.

Calculation of the spectrum of H proceeds exactly as earlier. Naturally we can also readily obtain the spectrum of other observables, e.g. of the angular momentum, Eq.(34), as well as the angular momentum content of the energy eigenstates. For example the last row of Table 4 gives the highest angular momentum which can be constructed in each sector with the bases generated so far.

### 5.4 Results

The spectrum of the theory is shown in Fig.5. A sample of states is labeled with their angular momenta. It is important that our cut-off preserves the SO(3) symmetry. Consequently, the basis described in previous section contains complete representations of the rotation group for any $`N_{cut}`$ . Hence the spectrum displayed in Fig.5 has appropriate degeneracies for each value of the angular momentum quantum number $`j`$.

To maintain some clarity of the figure we have cut arbitrarily the upper part of the spectrum, which extends to about $`E_{max}35`$ with current $`N_{cut}`$. One expects that the individual higher states have a considerable dependence on $`N_{cut}`$. However they also carry some relevant information, e.g. for quantum averages.

Apart from relating the bases in corresponding fermionic sectors, particle-hole symmetry also implies equality of all corresponding energy levels. Since our cut-off respects this symmetry we indeed observe, in sectors with $`F=6,5,4`$, exact repetition of the eigenvalues from $`F=0,1,2`$ subspaces. Therefore only first four sectors are displayed in Fig.5. Evidently the spectrum of D=4 theory is very rich and raises many interesting questions. We shall discuss some of them, beginning with the relation to the already known results.

#### 5.4.1 Correspondence with pure Yang-Mills quantum mechanics

It turns out that the fermionic part of the hamiltonian, Eq.(33), annihilates any state in the $`F=0`$ sector

$$H_F|0_F>=0.$$
(55)
Therefore the supersymmetric hamiltonian, Eq.(33) reduces in that sector to the hamiltonian of the pure Yang-Mills quantum mechanics, well known from the small volume approach . This offers us a chance to compare one special case of the present spectrum with the classical results of Ref.. Figure 6 shows such a comparison for the first three states as a function of $`N_{cut}`$. Nice convergence of our results towards those of Lüscher and Münster is seen. They used 120 states in each angular momentum channel and variationally optimized frequencies of the basis states. Our $`F=0`$ basis contains 238 states which span all angular momenta up to $`j=8`$, and has only 12 states with $`j=0`$. No variational adjustment of the frequencies was attempted. In view of this we consider above agreement as a rather satisfactory confirmation of the present approach. The assignment of the angular momentum, which is done automatically within our method, is also the same. In particular, the characteristic to the 0-volume, ordering of the tensor glueball and the first scalar excitation is reproduced.

#### 5.4.2 The cut-off dependence

The dependence of the spectrum on $`N_{cut}`$ is of course crucial, as it provides the quantitative criterion if the size of the basis is adequate. It is displayed in Fig. 7 for the first few lower levels in all fermionic sectors. It is clear that the cut-off dependence is different in various fermionic sectors. For $`F=0,1`$ (and $`F=6,5`$) the dependence is rather weak, suggesting that the energies of lower levels have already converged (within few percent) to their asymptotic values. On the other hand, in $`F=2,3,4`$ sectors convergence is considerably slower. As explained earlier, one can use the rate of convergence with $`N_{cut}`$ to distinguish between the continuous and discrete spectrum. It was proved that in the case of the free spectrum the convergence is O(1/$`N_{cut}`$) i.e. slow . For the localized bound states we conjecture that the rate of convergence with $`N_{cut}`$ is sensitive to the large $`x`$ asymptotic of the wave function. We have observed this in a simple anharmonic oscillator and other models. This regularity is also clearly confirmed in the Wess-Zumino and D=2 SYMQM models discussed in previous Sections.

Taking this into account we claim that the low energy spectrum of D=4 SYMQM is discrete in $`F=0,1,5,6`$ and continuous in the $`F=2,3,4`$ sectors. This is an interesting quantification of the result of Ref. which was mentioned earlier. Since the fermionic modes are crucial to provide continuous spectrum, it is natural that it does not show up in the sectors where they do not exist at all, or are largely freezed out by Pauli blocking.

Second result, which is evident from Fig.7, concerns the supersymmetric vacuum in this model. Assuming that indeed the eigen energies have approximately converged in $`F=0,1`$ sectors (none of them to zero) it follows that the SUSY vacuum can not be in empty and filled sectors ($`F=1,5`$ is also ruled out by the angular momentum) 121212In early attempts, empty and filled sectors were considered as possible candidates for SUSY vacuum.. The obvious candidates are the lowest states in $`F=2,4`$ sectors with their energy consistent with zero at infinite $`N_{cut}`$ , c.f. Fig. 7. It follows, together with the conjectured correlation between localizability and $`N_{cut}`$ dependence, that the SUSY vacuum in this model is non normalizable. Further evidence is based on the structure of the supersymmetric multiplets, and will be presented in the next section. Very recently van Baal studied $`F=2`$ sector in a more complex, non compact, supersymmetric model of the same family . Although our results cannot be directly compared, due to the specific boundary conditions required in , he finds that the energy of the lowest state in $`F=2`$ sector is indeed very close to zero. The deviations from exact zero are caused by the SUSY violating boundary conditions.

It is important to note that not all states in $`F=2,3,4`$ sectors have to belong to the continuum. Supersymmetry together with the discrete spectrum in $`F=0,1`$ channels implies existence of the normalizable states among the continuum of $`F=2,3,4`$ states. This will be discussed later, here we only give one explicit example. Indeed the energy of the $`F=2,j=1`$ state, shown in Fig. 7 (flatter curve beginning at E=6), has definitely weaker dependence on $`N_{cut}`$ than the others. We interpret this as a signature that the lowest $`|2_F,1_j>`$ state is localized. This situation may be a precursor of the more complex phenomenon expected in the D=10 theory. There, the zero energy, localized bound state of $`D_0`$ branes should exist at the threshold of the continuous spectrum. Present example suggests that one way to distinguish such a state from the continuum may be by the different $`N_{cut}`$ dependence.

#### 5.4.3 Supersymmetry

Operator level. To check the supersymmetry algebra

$$\{Q_\alpha ,Q_\beta \}=2\delta _{\alpha \beta }H+2g\mathrm{\Gamma }_{\alpha \beta }^kx_a^kG_a,$$
(56)
at finite cut-off, we sandwich both sides of this operator relation between our basis , introduce the intermediate complete set of states to saturate the product of $`Q`$’s and test the equivalent relations between matrix elements

$$<N|\{Q_\alpha ,Q_\beta \}|N^{}>=2\delta _{\alpha \beta }<N|H|N^{}>.$$
(57)
The gauge term in Eq.(56) does not contribute since the basis is gauge invariant at any $`N_{cut}`$. There is, however, one limitation. As is evident from Eq.(37)), supersymmetry generators change fermionic number by one, and can change number of bosonic quanta, $`B`$, at most by two

$$<N|Q_\alpha |M>0ifF_N=F_M\pm 1,|B_MB_N|2.$$
(58)
It follows that, if a cut off representation of unity $`|M><M|,B_M<N_{cut}`$ is used to saturate any matrix element of a product $`<N|Q|M><M|Q|N^{}>`$, it will be biased for high states in $`N,N^{}`$ sectors. Quantitatively, the supersymmetric relations, Eq.(57) are valid if

$$B_N,B_N^{}N_{cut}2.$$
(59)
This limits the number of matrix elements available for the test. Still many interesting predictions for lower states can be verified and should be satisfied exactly at finite $`N_{cut}`$.

One more property of supersymmetry generators should be kept in mind. Since they change $`F`$ by 1, they move between different fermionic sectors. Hence the generic matrix equation (57), when restricted to a particular sector, reads ($`\alpha =\beta `$, no sum)

$`<F|H|F^{}>=`$
$`<F|Q_\alpha |F1><F1|Q_\alpha |F^{}>+<F|Q_\alpha |F+1><F+1|Q_\alpha |F^{}>.`$
That is, the hamiltonian matrix receives contribution from a single fermionic sector, only if it is computed in empty or filled sectors. Of course above SUSY relations have many other consequences. Some of them will be used later, when discussing emergence of the supersymmetric multiplets.

Within the above limitations we have checked some consequences of SUSY algebra at finite $`N_{cut}`$. The spectrum of the $`Q^2`$ matrix agrees exactly with that shown in Fig.5 and derived directly form the hamiltonian. Moreover, the off-diagonal commutators indeed vanish in our basis. Many other interesting questions emerge and will be studied elsewhere.

SUSY on the level of states. The next goal is to identify supersymmetric multiplets in the spectrum and to asses the effect of the cut-off on their splittings and mixings. It is clear from Fig.5 that any direct search for approximately degenerate states is difficult and may not be conclusive. However our construction provides another simple way to move within the SUSY multiplets. Namely it is sufficient to use the explicit matrix representation of supersymmetry generators in bases summarized in Table 4. If supersymmetry was exact, acting with $`Q`$ on an eigenstate of $`H`$, would generate states within the same multiplet. At finite $`N_{cut}`$ this is no longer true, however we expect that the resulting state, $`Q|\mathrm{\Psi }>`$ say , will be spread around the appropriate multiplet member (or members) in another fermionic sector. This information will then be correlated with the spectrum of Fig. 5. Important simplification comes from the rotational symmetry, which is exact at every $`N_{cut}`$ in this approach. Supersymmetric generators carry the angular momentum 1/2 which clearly limits the range of possible targets $`Q|\mathrm{\Psi }>`$. Indeed, we confirm that

$$<F,j|Q|F\pm 1,j^{}>=0if|jj^{}|\frac{1}{2}.$$
(61)
To proceed, we denote the $`k`$-th energy eigenvector in the $`m`$-th fermionic sector and the $`n`$-th angular momentum channel by $`|m_F,n_j;k>`$. The subscript $`j`$ will be omitted where evident. A summary of various transitions we have analyzed is shown in Fig.8. For example, the action of $`Q_1`$ on the first state in $`F=0`$ sector gives

$$\frac{(|<1_F,1/2;1|Q_1|0_F,0_j;1>|^2+|<1_F,1/2;2|Q_1|0_F,0_j;1>|^2)}{Q_1|0_F,0_j;1>^2}=95\%,$$
(62)
which provides rather satisfactory confirmation of the first-glance guess from Fig.5. For higher states situation is more complex. The second scalar state in this sector goes under $`Q_1`$ into $`Q_1|0_F,0_j;2>`$ which decomposes in 40% into the $`|1_F,1/2;3,4>`$ doublet and 43% of yet higher, but also close in energy, $`|1_F,1/2;5,6>`$ doublet. This is not surprising, since higher states have not yet converged to their infinite $`N_{cut}`$ limit.

Even more interesting is to analyze SUSY relations between $`F=1`$ and $`F=2`$ sectors. Again the lowest $`|1_F,1/2;1>`$ state transforms predominantly (84%) into the third scalar $`|2_F,0_j;3>`$ which is closest in the energy. Another consequence of the SUSY algebra can be seen here. Naively one would expect that combining above transitions $`|0_F,0_j;1>|1_F,1/2;1>`$ and $`|1_F,1/2;1>|2_F,0_j;3>`$ would lead to the connection of $`F=0`$ and $`F=2`$ sectors and extension of the multiplet. This contradicts Eq.(57) which implies that the $`Q_1^2`$ is diagonal in the eigen basis of the hamiltonian. Solution of this apparent paradox is simple. There are two states with $`j=1/2`$ and while the multiplet beginning in $`F=0`$ sector ends on one linear combination, the second multiplet must start on the orthogonal one. In other words the supersymmetric images of $`|0_F,0_j;1>`$ and $`|2_F,0_j;3>`$ in $`F=1`$ sector must be orthogonal. We confirm this exactly (e.g. to the full precision of Mathematica numerics)131313This is because we use the limited basis Eq.(59), which guarantees Eq.(57) exactly.

$$<2_F,0_j|Q_1^{}Q_1|0_F,0_j>10^{16}.$$
(63)
The next (in the energy) $`|1_F,3/2;1>`$ state goes under $`Q_1`$ in 30% to $`|2_F,1_j;1,2,3>`$ triplet, with the rest shared among higher $`F=2`$ states. In fact, the $`|2_F,1_j;1,2,3>`$ states are those with the fast $`N_{cut}`$ convergence, which we have already pointed earlier as good candidates for localized states. Very little of $`|1_F,3/2;1>`$ goes back to the $`F=0`$ sector (c.f. fat head arrows in Fig.8)

$$\frac{Q_1|1_F,3/2;1>_{F=0}^2}{Q_1|1_F,3/2;1>^2}=6.3\%,$$
(64)
since there is no state with similar energy there 141414All surrounding states have already converged within few percent which suggests that the above small rate is also a O(1/$`N_{cut}`$) effect.. This means that the $`j=3/2`$ states start a new SUSY multiplet which extends into $`F=2`$ sector and not into $`F=0`$ one.

We are using only one of the four SUSY generators just as an example. Similarly we do not consider systematically all members of $`j0`$ multiplets. The full analysis, which in fact would be quite straightforward in the present setup, would involve complete discussion of the extended, $`N=4`$, SUSY algebra. However it is out of the scope of this article and will be done elsewhere.

All above examples have a common feature. Clear SUSY regularities show up if the states have already converged, i.e. if their energy depends weakly on $`N_{cut}`$ . This is not the case for the lowest states in $`F=2,3,4`$ sectors which, we believe, form a continuum at infinite cut-off. Accordingly, and similarly to the D=2 SYM QM case, identification of multiplets and the SUSY vacuum in this sector is not clear at the moment. As one possible signature of the vacuum one might take the average norm of the supersymmetric image of the lowest state $`Q_1|2_F,0_j;1>^2/dim[3_F]`$. Indeed it is smaller than that for other states and falls with $`N_{cut}`$. However it is too early to draw definite conclusions.

Clearly we have only begun. The aim here is just to introduce the new approach and its potential 151515We know of no other method capable to provide such a detailed and quantitative characteristics as discussed in this Section., leaving specific applications for more focused articles.

Witten index. The total number of bosonic and fermionic states is the same for each $`B`$, c.f. last column of Table 4. This is not a direct consequence of supersymmetry, but rather of a combinatorics of binary fermionic systems, which is not spoiled by gauge invariance. Total number of states gives the Witten index at $`T=0`$. Hence we expect

$$\underset{N_{cut}\mathrm{}}{lim}I_W(0,N_{cut})=0.$$
(65)
As was already discussed Witten index is discontinuous at $`T=0`$. In particular $`I_W(0)`$ depends on the regularization, therefore this number is tied to our B - F symmetric scheme of increasing the basis.

Figure 9 shows the (euclidean) time dependence of the Witten index calculated from our spectrum with up to five bosonic quanta. Clearly we are far from the convergence, but some interesting signatures can be already observed. First, one can see how the discontinuity at $`T=0`$ emerges even at this early stage. Second, the exponential fall-off common to all curves at large $`(T>5)`$ times is evident. This is where the lowest state dominates. Since at this $`N_{cut}`$ its energy is not zero, the fall-off is exponential. Finally, and most interestingly, we observe the flattening shoulder appearing at $`T23`$. This is the signal of the supersymmetric cancellations which occur on the average, even though the exact SUSY correspondence between individual states does not yet appear. The behavior with $`N_{cut}`$ is not inconsistent with the exact bulk value 1/4  obtained also from the nonabelian integrals . Clearly higher $`N_{cut}`$ are needed, and, what is important, they are perfectly within the range of present computers.

## 6 Summary

The new approach to quantum mechanical problems is proposed. The Hilbert space of quantum states is algebraically implemented in the computer code. In particular realization, used here, states are represented as Mathematica lists. All basic operations on quantum states are mirrored as definite operations on above lists. Any quantum observable is represented as a well defined function on these lists. This allows for automatic calculation of matrix representation of a hamiltonian and other quantum operators of interest. To this end we use the discrete eigen basis of the operators of the numbers of quanta. The length of above lists is not fixed. Similarly to the length of an arbitrary quantum combination of basic states, it can vary dynamically allowing for any number of quanta.

Of course in any finite computer the maximal length of a combination must be limited. Therefore we impose a cut-off $`N_{cut}`$ which bounds from above the number of allowed quanta. The stringent and quantitative test of this approach is provided by checking the dependence of any physical observable, e.g. the spectrum or the wave functions, on the cut-off. This is similar to studying the infinite volume statistical systems via the finite size scaling. We have applied above technique to the three progressively more complicated systems: Wess-Zumino quantum mechanics, supersymmetric Yang-Mills quantum mechanics in D=1+1 dimensions, and D=4 SYM quantum mechanics, both for the SU(2) gauge group. In distinction from many other approaches (for example lattice simulations) the method is completely insensitive to the sign problem and works equally well for bosonic and fermionic systems.

For Wess-Zumino quantum mechanics we can calculate the discrete spectrum for any values of parameters. Clear restoration of the supersymmetry was observed for the cut-offs well within the capabilities of present computers. Witten index was also obtained and its convergence to the known result $`I_W(T)=2`$ is clearly seen.

The next system, D=2 SYM QM, possesses the gauge invariance which was readily incorporated in our approach. The physical subspace of gauge invariant states was explicitly constructed. Known structure of the solutions in terms of four fermionic sectors was reproduced. Convergence of the method, and emergence of the supersymmetry, was also studied in this more difficult case of the continuum spectrum. Witten index, restricted to one supersymmetric branch of the model, was defined and computed for the first time. Clear, albeit slow, convergence to the time independent fractional value $`I_W^R(T)=\pm 1/2`$ was observed. Moreover, we have found a special scheme of increasing the basis such that the supersymmetry is exact at any finite cut-off. This however, may be related to the exact solubility of the model.

Finally, the method was applied to the unsolved up to now D=4 SYMQM. This much richer system has the SO(3) rotational symmetry inherited from its space extended predecessor. Our approach preserves this symmetry exactly at any finite cut-off. The Hilbert space splits again into seven sectors with fixed fermionic number. We have obtained the complete spectra in all these sectors and studied their cut-off dependence. The spectrum in $`F=0`$ sector agrees with the classical 0-volume calculation for pure Yang-Mills quantum mechanics . An efficient method to distinguish between the discrete and continuous spectrum was proposed. It turns out that the asymptotics of the wave function at large distances determines the convergence of our calculations with the number of allowed quanta $`N_{cut}`$. The continuous spectrum with non localized wave functions converges slowly $`(O(1/N_{cut})`$, while discrete, localized bound states lead to faster (sometimes even exponential) convergence. Accordingly, we have found an evidence that the spectrum of D=4 SYM QM is discrete in $`F=0,1,5,6`$ sectors while it is continuous in the $`F=2,3,4`$ sectors. This provides an explicit realization of the claim of Ref. in particular fermionic sectors. Interestingly, localized states exist also in the sectors with continuous spectrum. This is a simple consequence of supersymmetry whose transformations move between adjacent fermionic sectors. Our method allows to monitor directly the action of SUSY generators and analyse supersymmetric images of any state. In this way we have identified some candidates for lowest supersymmetric multiplets. They do not have the same energy at current values of the cut-off, however the splittings are small and consistent with vanishing at infinite $`N_{cut}`$. In particular, the 0-volume glueballs found in $`F=0`$ sector are relevant to fully supersymmetric theory in that there exist gluino-gluon bound states with the same masses. We have also found a candidate for the supersymmetric vacuum which seems to belong to the continuum. However identification of SUSY multiplets in the continuum part of the spectrum requires higher $`N_{cut}`$. Supersymmetry was also tested on the operator level. In particular we confirmed that the spectrum of the hamiltonian coincides with that of $`Q_1^2`$. As a final application we have calculated the Witten index for this theory. It still depends strongly on $`N_{cut}`$ . Nevertheless an early evidence for some of its asymptotic properties can be already seen.

## 7 Future prospects

The main goal of this work is to asses the feasibility of attacking with the new approach the BFSS model of M-theory. Certainly the method is more efficient for smaller number of degrees of freedom. Application to the Wess-Zumino quantum mechanics and D=2 Yang-Mills quantum mechanics give quite satisfactory and quantitative results including the new intriguing scheme which preserves exact supersymmetry for finite cut-off. For more complex, and correspondingly richer, D=4 SYMQM the method is able to provide new results including detailed information about the supersymmetric structure of the spectrum and the observables. Obviously we also see a room for improvement which is especially needed in the continuum sector. However, and that we find most important, further increase of the size of the Fock space is possible within the available technology. The whole programme can be (and is being) implemented in the standard, compiler based languages which usually improves performance by a factor 10-100. Further, the action of the quantum operators can be optimized taking into account symmetries of the states. Finally, one can go for more powerful computers.

Taking all above into account, one can reasonably expect substantial improvement in the quality of the present D=4, N=2 results. At the same time one should be able to study D=4 systems with higher N. To answer the main question: yes, we think that the quantitative study of the $`D=10`$ theory is feasible, and reaching current quality results for the D=10 is realistic, beginning with the SU(2) gauge group. $`D=10`$ hamiltonian and $`spin(9)`$ generators do not conserve fermionic number . The reason for this complication should be better understood in the first place.

Apart from increasing the number of degrees of freedom, higher gauge groups pose an interesting problem of constructing gauge invariant states. This has already been done in the small volume approach for SU(3) , and should not present any fundamental problem for higher N as well. At the same time the possibility of some large N simplification, specifically within the present approach, should be investigated.

Last, but not least, we would like to mention a host of applications of this method to other quantum mechanical systems. For example one may simply extend the 0-volume calculations to the full non supersymmetric QCD with dynamical quarks in fundamental representation. Apart from known glueballs, this would give us a spectrum of quark-made-hadrons in the ”femtouniverse”.

As another rather different application, we mention that this program is already being used as a routine in solving, to a high precision, quantum mechanics of a simple two dimensional building blocks of a prototype quantum computer . In particular, a complete quantum evolution in time with the time dependent hamiltonian, was straightforward to simulate.

The D=4, N=2 SYMQM studied here has 15 degrees of freedom. There are many unsolved quantum mechanical systems with this or smaller complexity. Present approach should be well applicable in some of these cases.

Obviously there are many routes which can be followed from this point and we are looking forward to explore some of them.

Acknowledgments. I would like to thank C. M. Bender for an instructive discussion which inspired this approach. I also thank P. van Baal, P. Breitenlohner, L. Hadasz, M. Rostworowski and H. Saller for intensive discussions. This work is supported by the Polish Committee for Scientific Research under the grant no. PB 2P03B01917.


format: text

text:
# Local transformation of mixed states of two qubits to Bell diagonal states

## ACKNOWLEDGMENTS

This work was supported in part by the Postdoctoral Science Foundation and the special funds for Major State Basic Research Project No. G001CB3095 of China.