the schrdinger equation of an interacting many - body system can be solved analytically only for specific interaction potentials and for a very restricted number of particles .
a wide range of approximation methods has been developed to determine the ground and excited states of many - body systems found in nuclear , atomic , and condensed matter systems . here
, we focus on quasi two - dimensional electronic systems , which are experimentally realized in quantum - dots . quantum dots are often embedded in layered semiconductor structures .
the strong confinement along the vertical direction perpendicular to the layer suggests an effectively two - dimensional description along the remaining two lateral dimensions . besides the electron - electron interactions ,
an additional external confinement potential along the lateral directions is created by etching and gating of the semiconductor device .
quantum dots can be regarded as _ artificial atoms _ , albeit with a central potential different from the nuclear coulomb attraction .
their excitation spectra are probed by electronic and optical measurements , also under the influence of additional magnetic fields , which gives rise to a landau - level structure .    with the discovery of the fractional quantum hall effect ( fqhe ) in systems of interacting electrons in quasi two - dimensional systems , many theoretical approaches for studying interacting electrons
have been proposed , including the celebrated laughlin wave - function for electrons in the lowest landau level . in the
experimentally realized hall devices translational symmetry is broken by the current source and drain contacts , which lead to the formation of hot - spots with high electric field values @xcite .
a comparison of the various theoretical approaches with numerical methods is often performed for few - electron quantum dots .
analytic @xcite and numerical @xcite solutions for the two electron case show the interplay of coulomb interaction and confinement potential , leading to alternating spin polarization of the ground state in a quantum - dot as function of magnetic field @xcite . by separation of the two - electron case into the center - of - mass motion and the relative part ,
semiclassical solutions have been constructed @xcite and compared to the exact solutions @xcite . for three electrons , the relative - coordinate basis set and
expressions for the coulomb matrix elements are given in @xcite .
in the popular _ configuration interaction _ ( ci ) approach a diagonalization of the hamiltonian including the coulomb - interaction is done with respect to a product basis of single - particle slater determinants .
convergence is generally checked by systematically enlarging the size of the basis set , for instance by increasing the value of the cut - off energy up to which all states are included for the non - interacting system @xcite .    with increasing number of basis states , the computation of all interacting matrix - elements and the final diagonalization of the hamiltonian quickly exhausts computer resources .
thus a reduction of the basis set while maintaining accuracy is required for studying the excitation spectra of few body systems .
if the reduction is based on underlying group - theoretical properties , an exact decomposition of the basis set into irreducible blocks is possible .
a diagonalization of the smaller blocks is more easily achievable and the decomposition provides physical insight into the system .
one example is the separation of the spatial part of the wavefunction from the spinor states for electrons .
if we require that the product of spatial and spin - part represents electrons ( antisymmetric states ) , the basis set is decomposed into separate spin states .
this results in spin - adapted ( sa ) configuration methods , for instance implemented in the saci program code @xcite .    at this point saci methods used in condensed matter theory for
more than 3 electrons stop and do not exploit additional properties of the basis set , such as a decomposition into relative and ( interaction - free ) center - of - mass ( _ cm _ ) states .
the separation of the basis set into a product state of relative and center - of - mass part is of particular use if the hamiltonian of the @xmath0-electron system is given as a sum of operators acting on center - of - mass and relative coordinates respectively . for electrons confined to a two - dimensional plane perpendicular to a uniform magnetic field
this separation of the hamiltonian is possible as long as the confinement potential in the two - dimensional plane does not lead to an admixing of relative and center - of - mass coordinates @xcite . for molecular or atomic systems moving in three dimensions across uniform magnetic fields , additional couplings between relative and center - of - mass coordinates occur , i.e.  the motional stark effect @xcite .
the importance of a systematic separation of relative coordinates from the center - of - mass excitation for studying correlated two - dimensional electron systems for three electrons was pointed out early by laughlin @xcite and noticed recently again for more than three electrons @xcite .
laughlin s trial wave function by construction is expressed only in relative coordinates and can be interpreted and numerically handled similar to a classical plasma of interacting particles @xcite .
we note that both in many - electron and nuclear physics , the interactions between the particles are exclusively functions of relative coordinates .
then the same shift vector applied to any single particle coordinate does not change the hamiltonian . by the well - known noether theory of symmetry it then follows that the total momentum @xmath1 of the system is conserved . for an unusual but powerful example of a conserved quantity we refer to @xcite .
relative coordinates , in contrast to the single particle coordinates , are free of the @xmath2 coordinate . in quantum mechanics ,
the components of @xmath1 are proportional to the partial derivatives of the @xmath2 coordinate .
therefore , if the states depend on relative coordinates only , we have @xmath3 . another possible option which we shall adopt
is to demand the @xmath2 state to have zero oscillator excitation .
the paradigms of nuclear theory show @xcite that relative coordinates are the natural tools to understand and compute nuclear scattering and reactions . in nuclear structure ,
such coordinates may be adapted and serve to describe few - particle systems @xcite , clustering , and decay processes .
the construction of orbital states can be achieved by use of group and representation theory @xcite .
the orbital states are coupled with spin states to antisymmetric states .
the adaption to orbital permutational symmetry for composite systems in relative coordinates becomes more involved , but was extensively studied in @xcite .
we report here work on 2d electrons in an external perpendicular magnetic field , so that the 2d total momentum is still conserved @xcite .
we use harmonic oscillator states and the methods of bargmann hilbert space .
the electron - electron coulomb repulsion is computed in matrix form and included into the diagonalization . for @xmath4 electrons
we use the relative tetrahedral coordinates from @xcite with simple properties under permutations .
the spin states for total spin @xmath5 are explicitly determined by projection operators .
+ single - particle coordinates are used to construct antisymmetric states as slater determinants of single - particle states .
if these states are taken as harmonic oscillator states , slater determinants can not avoid the intrusion of @xmath2 excitations into state space which give raise , in addition to the proper states , to so - called spurious states of spurious angular momentum .
the true role of spurious states can be seen only by the parallel use of relative and @xmath2 coordinates .
we shall present in section 7 in particular for total spin @xmath6 , a complete analysis of these spurious states .
their intrusion of state space is demonstrated in the spectrum of determinantal eigenstates of the standard saci approach @xcite . moreover , if determinantal eigenstates are variationally computed in a subspace defined by a fixed maximal oscillator excitation , we find that the admixture of spurious states distorts this subspace and the variational results . in our proper treatment ,
any eigenstate in relative coordinates is the basis to a set of equidistant spurious states of increasing @xmath2 excitation .
the total spin of @xmath0 electrons ranges from the maximum @xmath7 to the minimum value , @xmath8 for @xmath0 odd and @xmath6 for @xmath0 even . the permutational symmetry of the spin state is characterized by a partition @xmath9,\ : f_1^s+f_2^s = n,\ : f_1^s \geq f_2^s$ ] of two rows @xmath10 $ ] only , because there are for any electron only the two spin ( up , down ) states .
the total spin @xmath11 is related one - to - one to the spin partition @xmath12 by @xmath13 since the electrons are fermions , the combined orbital plus spin state must be antisymmetric with partition @xmath14 $ ] . then is follows that the only allowed orbital partition @xmath15 is associate to the spin partition @xmath12 , that is , its young diagram consists of @xmath0 square boxes , arranged in @xmath16 rows of length 2 , and @xmath17 rows of length 1 . to the orbital partition @xmath15 there
belong orthogonal partner functions whose number we denote by @xmath18
. they may be labelled by young tableaus @xmath19 .
any young tableau is a filling of the boxes of the young diagram @xmath15 by the numbers @xmath20 , such that the numbers increase both to the right and downwards .
the application of permutations to an orbital state of @xmath0 particles , denoted as @xmath21 , is completely determined by the tableau @xmath22 .
the antisymmetric fermion state , with respect to permutations , takes the form of a sum of products of orbital and spin states , @xmath23\alpha^n   s \rangle= |f|^{-\frac{1}{2 } } \sum_r |\alpha^n f r\rangle |s f^s r^s\rangle,\ ] ] where @xmath24 is the young tableau of the spin partition @xmath12 , determined from the orbital tableau @xmath22 by reflection in the main diagonal of its orbital young diagram .
the operators @xmath25 for the external magnetic field , the central field and the interactions are all invariant under permutations and independent of the spin . since the spin basis states @xmath26
are orthogonal , the matrix elements of all spin - independent operators @xmath25 between antisymmetric states reduce to the spin - free sum of orbital matrix elements @xmath27 here the total spin @xmath11 is fixed from eq .  [ es1 ] by half the difference of the column length of the orbital partition @xmath15 . from eq .
[ es3 ] it follows that we never need the spin states but must compute @xmath18 orbital diagonal matrix elements .
_ example _ : for @xmath4 particles , the possible values of the total spin are @xmath5
. the corresponding spin partitions are @xmath28,[31],[22]$ ] , and the orbital partitions are @xmath29 ,   [ 211],[22 ] $ ] of dimensions @xmath18 = 1 , 3 , 2 respectively .
so we need the construction of @xmath18 orbital states with an orbital partition @xmath15 , fixed by the total spin @xmath11 , and by its young tableaus @xmath22 .
to project states of orbital symmetry @xmath15 with the young tableau @xmath30 , we use the matrix elements of the representation @xmath31 which are determined by the tableaus .
a young diagram is an arrangement of square boxes in rows of length @xmath32 . in a young tableau @xmath22
, we label these boxes by @xmath20 such that the numbers increase to the right and downwards .
there are exactly @xmath18 young diagrams for the partition @xmath15 .
the row and column of the box @xmath33 are then @xmath34 . for a pair of boxes @xmath35
we define the axial distance @xcite p.28 by @xmath36 we first determine the matrices of the representation for the @xmath37 transpositions @xmath38 of @xmath39 .
these transposition by multiplication generate all elements of @xmath39 , and so all representation matrices can be built by matrix multiplication from the representations of these transpositions .
for the permutation which interchanges @xmath33 with @xmath40 and the young tableau @xmath22 , there may be a second tableau @xmath41 with boxes @xmath35 interchanged .
if there is no such second tableau , the two boxes must appear in succession in the same row or the same column of the young tableau .
in these two cases the matrix element of the representation are given by @xmath42 otherwise , the non - trivial part of the representation is given by the @xmath43 matrix @xmath44.\ ] ]    the operator of the transposition @xmath45 interchanges the coordinates of the single particle indexed by the number @xmath46 with those of the particle indexed by @xmath47 .
following the wigner @xcite pp.102 - 6 prescription , to assure a homomorphic action we must associate to any permutation @xmath48 the inverse action @xmath49 of the permutation on the particle indices .
the young operators @xcite pp .
25 - 26 are given by @xmath50^{\dagger}= c^f_{sr}.\ ] ] these operators must be applied to the highest polynomials eq .
[ es21 ] of fixed pseudo angular momentum @xmath51 ( see sect .  [
sec : pseudo ] ) , provided that from @xcite we know that @xmath51 can yield the orbital partition @xmath15 . if a young operator gives zero when applied to the polynomial , we first apply a lowering operator eq .
[ es23 ] and then again try the young operator . exhausting all values of the pseudo angular component we must find all the states of orbital symmetry .    for an antisymmetric state eq .
[ es2 ] we use orbital transposition operators @xmath52 to write @xmath53s|\sum_{s < t } v(s , t)|\alpha^n[1^n]s\rangle   \\ \nonumber & & = \frac{1}{|f|}\sum_{s < t } \sum_r \langle f r|t(s,1 ) t(t,2 ) v(1,2)t(t,2)t(s,1)|f r\rangle \\
\nonumber   & & = \frac{1}{|f| } \sum_{s
< t } \sum_{r}\sum_{r',r '' }   \langle fr| t(s,1 ) t(t,2)|fr'\rangle   \langle f r'|v(1,2)|f r''\rangle \langle fr''|t(t,2)t(s,1)|fr\rangle \\ \nonumber   & & = \frac{n(n-1)}{2 } \frac{1}{|f|}\sum_r \langle f r|v(1,2)|f r\rangle . \end{aligned}\ ] ] the third line of eq .
[ es13e ] follows from the orthogonality of the spin states in eq .
[ es2 ] . in the last line we use the involutive property of transpositions and the orthogonality of young tableaus ,    @xmath54 ^ 2=[t(t,2)]^2=1 ,    \\ \nonumber & & \langle f r''|fr'\rangle   = \delta_{r'',r'}.\end{aligned}\ ] ]    * prop * : the matrix element of orbital twobody interactions , invariant under permutations , between antisymmetric states of fixed total spin @xmath11 , can be reduced to a sum over diagonal orbital matrix elements of the interaction of the first pair between states characterized by young tableaus , multiplied by the number of pairs .    to construct states with other tableaus of the same partition
, we provide a ladder procedure by employing the properties of the representation .
given an initial basis function @xmath55 , we know from eq .
[ es12 ] that the permutation @xmath56 connects it at most to a second basis function @xmath57 whose young tableau differs by an interchange of the numbers @xmath56 .
we rewrite eq .
[ es12 ] in the form @xmath58 therefore once we know how to apply @xmath59 to the state @xmath55 , we can construct @xmath57 as the linear combination eq .
[ es14 ] , second line . with the operator eq .
[ es13 ] we have constructed the state with the highest young tableau .    * prop * : with steps as in eq .
[ es14 ] we can generate all other young tableaus from an initial one .
examples are given in the following subsections .
in an orbital young operator eq .
[ es13 ] , the sum runs over all permutations of @xmath39 .
this sum may be broken into classes @xmath60 , characterized by @xmath61 independent and commuting cyclic permutations of length @xmath62 where @xmath63 any cyclic permutations in turn can be factorized into transpositions as @xmath64 and any transposition is conjugate to one of the @xmath0 generators @xmath65 of @xmath39 , see eq .
[ es28 ] . by the representation @xmath66 ,
all the factorizations of permutations are converted into matrix factorizations .
a systematic enumeration of all @xmath67 permutations can go in terms of classes @xmath60 of conjugate elements .
the number of permutations in the class @xmath60 corresponding to eq .  [ es14a ] according to weyl @xcite p.329 is @xmath68
for a single electron in 2d , we label its state by the oscillator excitation @xmath69 , by the numbers @xmath70 of right and left circular oscillator quanta , and by the angular momentum @xmath71 .
note the following relations for coordinates and derivatives of a single particle in 2d , passed to a circular setting : @xmath72 and the quadratic relations @xmath73 the last relations are easily extended to scalar products of several vectors , labelled and contracted wrt lower indices .
these yield generators of the symplectic group @xmath74
.    the states of @xmath0 electrons are then labelled by the summed up number @xmath75 of oscillator quanta and by the summed up angular momentum @xmath76 it is desirable to find other ( integer ) orbital quantum numbers which if possible can label at least in part an orthogonal set of orbital states .
for this purpose we use the scheme of @xcite whose quantum numbers arise from group representations .
the unitary group of all orbital degrees of freedom is the unitary group @xmath77 , its commuting subgroups acting on the 2 orbital or the @xmath37 relative motion degrees of freedom are @xmath78 .
their representations are both given by the same two - component integer partitions @xmath79
, h_1+h_2=n,\ : h_1\geq h_2\geq 0 $ ] .
the representation of the subgroup @xmath80 is characterized by the angular momentum @xmath81 . for given partition @xmath79 $ ] the angular momentum ranges from its maximum @xmath82 to its minimum value @xmath83 or @xmath84 for @xmath85 odd or even respectively .      in the 2d magnetic field ,
the lowest landau level @xmath86 arises if the summed up angular momentum @xmath81 equals the summed up excitation @xmath75 . in terms of the partition
@xmath79 $ ] for the unitary group @xmath87 , we have @xmath88 .
the weights , i.e. the degrees @xmath89 of the vectors indexed @xmath90 , are restricted to the range @xmath91 , and the maximum angular momentum is @xmath92 .
the lowest landau level therefore enforces @xmath93 .
this condition selects among the partitions of @xmath75 the single case @xmath79=[n0 ] , l = n$ ] . only for this partition
can the excitation @xmath75 be carried by a single vector . for higher landau levels
this condition is relaxed and new unitary partitions are allowed .
the group @xmath94 has an orthogonal subgroup @xmath95 whose representation label we term the pseudo angular momentum @xmath51 .
the possible values of @xmath51 were given in @xcite . here
the authors use the subgroup chain @xmath96 with branching of representation labels @xmath97 > ( \lambda',\mu ' ) > \lambda,\ : ( \lambda',\mu')=(h_1-h_2 , h_2).\ ] ] the range of these labels is determined by specific branchings from the representation in the next higher subgroup which we shall discuss in more detail .
we treat the oscillator states by the methods of bargmann hilbert space of analytic functions @xcite .
we specifically note that in what follows we must _ interchange the role of the groups and subgroups _
@xmath98 compared to @xcite , and _ convert angular into pseudo angular momentum .
_ this allows us to take advantage of the group theory given in @xcite , but in a new and different interpretation .    by _ upper indices _
@xmath99 ( or upper @xmath100 ) , see below , we denote the cartesian vector components in 2d , and by _ lower cartesian indices _
@xmath101 or spherical indices @xmath102 three ( or more ) relative vectors of the four ( or more ) particles .
the normalized states of angular momentum @xmath103 of a single particle with oscillator excitation @xmath75 in the bargmann representation are @xmath104^{\frac{1}{2}}.\ ] ] with @xmath105 a solid spherical harmonic @xcite p.69 .
the coordinates may then be taken as the six linear combinations of cartesian coordinates @xmath106 ,   \\
\nonumber   & & z^+_-= \sqrt{\frac{1}{2}}[(z_1 ^ 1-iz_2 ^ 1)+ i(z_1 ^ 2-iz_2 ^ 2)],\ : z^+_0=z^+_3 ,   \\
\nonumber   & & z^-_+= \sqrt{\frac{1}{2}}[(z_1 ^ 1+iz_2 ^ 1)- i(z_1 ^ 2+iz_2 ^ 2 ) ] ,   \\ \nonumber   & & z^-_-= \sqrt{\frac{1}{2}}[(z_1 ^ 1-iz_2 ^ 1)- i(z_1 ^ 2-iz_2 ^ 2 ) ] , \ : z^-_0=z^-_3 .   \\
\nonumber \end{aligned}\ ] ] the lower cartesian indices @xmath107 denote three relative coordinates like the jacobi or tetrahedral ones .
the permutations of particles act on these coordinates and introduce linear combinations wrt the lower indices .
next we introduce from @xcite p.187 for @xmath4 the elementary polynomials @xmath108_+ \\ \nonumber & s:=z_+^1z_-^1 + z_0 ^ 1z_0 ^ 1 , \\
\nonumber   & t:= ( z^1\cdot z^1)(z^2\cdot z^2)-(z^1\cdot z^2)^2 , \\ \nonumber \end{aligned}\ ] ] in @xcite p. 187
we find in the scheme eq .
[ es16 ] the linearly independent polynomials @xmath109 here all numbers in exponentials must be non - negative , and @xmath110 is an integer which labels the multiplicity of a repeated pseudo angular momentum @xmath51 for a given partition @xmath79 $ ] of @xmath94 .
these multiple states of equal pseudo angular momentum are indexed by @xmath110 and linearly independent , but not yet orthogonal .
orthogonality must be achieved by standard matrix methods .
the permutation group @xmath111 acting on relative coordinates is a subgroup of the pseudo rotation group @xmath95 .
the multiplicity of its representations , see table  [ table1 ] , was given for all orbital partitions @xmath15 of @xmath112 particles in @xcite .
the polynomials eq .
( [ es21 ] ) have the maximum component @xmath51 of the pseudo angular momentum .
lowering of the component is achieved by applying the first order lowering differential operator @xmath113    finally we can also lower the angular momentum @xmath81 from its maximum value @xmath114 in steps of @xmath115 by application of the angular momentum operators @xmath116      @xmath117 & [ 211 ] & [ 22 ]   & [ 31 ] & [ 4 ]   & \lambda^{\pi}[1 ^ 4]&\lambda^{\pi}[4]\\ \hline   \lambda^+ & [ 4 ]    & [ 31 ] & [ 22 ]    & [ 211]&[1 ^ 4 ] &   & \\ \hline    0         & 1    &    &    &    &    & 0 ^ -&0^+\\   1         &      &    &    & 1   &    &    & \\   2         &      & 1   & 1   &    &    &    & \\   3         &      & 1   &    & 1   & 1   & 3^+&3 ^ -\\   4         & 1    & 1   & 1   & 1   &    & 4 ^ -&4^+\\   5         &      & 1   & 1   & 2   &    &    &    \\   6         & 1    & 2   & 1   & 1   & 1   & 6^+6 ^ -&6 ^ -6^+\\   7         &      & 2   & 1   & 2   & 1   & 7^+&7 ^ -\\   8         & 1    & 2   & 2   & 2   &    & 8 ^ -&8^+   \\   9         & 1     & 2   & 1   & 3   & 1   & 9^+9 ^ -&9 ^ -9^+\\   10        & 1     & 3   & 2   & 2   & 1   & 10 ^ -10^+&10^+10 ^ -\\   11        &      & 3   & 2   & 3   & 1   & 11^+&11 ^ - \\   12        & 2     & 1   & 1   & 1   & 1   & ( 12 ^ -)^212^+&(12^+)^212 ^ -\\   13        & 1     & 1   & 1   & 2   & 1   & 13 ^ -13^+&13^+13 ^ -\\   14        & 1     & 2   & 2   & 1   & 1   & 14 ^ -14^+&14^+14 ^ -\\   15        & 1     & 2   & 1   & 2   & 2   & 15 ^ -(15^+)^2&15^+(15-)^2\\   16        & 1     & 2   & 2   & 2   & 1   & 16 ^ -16^+&16^+16 ^ -\\   17        & 1     & 2   & 2   & 3   & 1   & 17 ^ -17^+&17^+17 ^ -\\   18        & 2     & 3   & 2   & 2   & 2   & ( 18 ^ -)^2(18^+)^2&(18^+)^2(18 ^ -)^2\\   19        & 1     & 3   & 2   & 3   & 2   & 19 ^ -(19^+)^2&19^+(19 ^ -)^2\\   20        & 2     & 3   & 3   & 3   & 1   & ( 20 ^ -)^220^+&(20^+)^2 20 ^ -\\   21        & 2     & 3   & 2   & 4   & 2   & ( 21 ^ -)^2(21^+)^2&(21^+)^2(21 ^ -)^2\\   22        & 2     & 4   & 3   & 3   & 2   &
( 22^+)^2(22 ^ -)^2&(22 ^ -)^2(22^+)^2\\   23        & 1     & 4   & 3   & 4   & 2   & 23 ^ -(23^+)^2&23^+(23 ^ -)^2\\   24        & 3     & 2   & 2   & 2   & 2   & ( 24 ^ -)^3(24^+)^3&(24^+)^2(24 ^ -)^2\\ \hline   \end{array } $ ]    the orbital partitions @xmath118,[4]$ ] do not apply for electrons with spin @xmath119 .
the entries for higher values of @xmath51 than given in table [ table1 ] are found by the recursive methods from @xcite eq .
a full orbital 4-particle state is characterized by the quantum numbers @xmath120 , \lambda , l ,   f r\rangle,\ ] ] where @xmath22 ranges over the orbital young tableaus of the diagram @xmath15 .
the total spin @xmath11 is determined by antisymmetry and by the orbital partition @xmath15 , see section 2 .
the orbital state has partition @xmath121 $ ] of dimension @xmath122|=1 $ ] .
the operator eq .
[ es13 ] applied to the initial state @xmath123 becomes @xmath124 is the sign of the permutation @xmath48 , @xmath125 for even and @xmath126 for odd permutations .
the orbital partition is @xmath127 $ ] of dimension @xmath128|=3 $ ] .
the three young diagrams ( written for simplicity in rectangular frames ) read @xmath129 , & r_2=\left [ \begin{array}{ll } 1 & 3\\ 2&\\ 4&\\ \end{array } \right ] , & r_3=\left [ \begin{array}{ll } 1 & 4\\ 2&\\ 3&\\ \end{array } \right]\\ \end{array}.\ ] ] assume that the initial state @xmath130 r_1\rangle$ ] has been generated by application of a young operator eq.[es13 ] .
we find the two partner states of @xmath130r_1\rangle$ ] by the operator relations eq [ es14 ] and the axial distances , @xmath131 r_2\rangle = |\alpha^4[211]r_1\rangle   \sqrt{\frac{1}{3 } }            + t(2,3)|\alpha^4[211]r_1\rangle \sqrt{\frac{4}{3}},\:\ : \tau_{23}(r_1)=-2,\\ \nonumber             & |\alpha^4[211 ] r_3\rangle = |\alpha^4[211]r_2\rangle \sqrt{\frac{1}{8 } }            + t(3,4)|\alpha^4[211]r_2\rangle \sqrt{\frac{9}{8}},\:\ : \tau_{34}(r_2)= -3.\end{aligned}\ ] ]      the orbital partition is @xmath132 $ ] of dimension @xmath133|=2 $ ] .
the two young diagrams read @xmath134 ,   r_2=\left [ \begin{array}{ll } 1&3\\ 2&4\\ \end{array } \right].\ ] ] again we generate the initial orbital state @xmath135r_1\rangle$ ] by application of a young operator eq .
we find the second state by the operator relation and the axial distance @xmath136r_1\rangle \sqrt{\frac{4}{3}},\:\ : \tau_{23}(r_1)=-2.\ ] ]
since we are dealing with particle - particle interactions , the _ cm _ vector is not affected , the total momentum is conserved , and the system is amenable to a description in relative coordinates .
a standard choice are the jacobi relative coordinates @xmath137    in @xcite there was used a particular tetrahedral set of three doubledot relative coordinates .
this set has the virtue that under the application of any permutation these relative coordinates are at most permuted , eq .
the transformation properties of the doubledot coordinates are summarized in @xcite table 3 .
the relative tetrahedral coordinates @xmath138 from @xcite , but now enumerated in the order of @xcite for @xmath4 , are    @xmath139    it is found from @xcite , table  3 , that the tetrahedral coordinates eq .
[ es30 ] have extremely simple matrix transformations under the generating transpositions @xmath140 : @xmath141 = \left[\begin{array}{rrr }    1&0&0\\    0&0&-1\\    0&-1&0\\    \end{array } \right ] \left[\begin{array}{l }        \eta_1\\        \eta_2\\        \eta_3\\        \end{array}\right ] , \\
\nonumber   & & t^{\eta}(2,3)\left[\begin{array}{l }        \eta_1\\        \eta_2\\        \eta_3\\        \end{array } \right ] = \left[\begin{array}{lll }    0&1&0\\    1&0&0\\    0&0&1\\    \end{array } \right ] \left[\begin{array}{l }        \eta_1\\        \eta_2\\        \eta_3\\        \end{array}\right ] , \\
\nonumber   & & t^{\eta}(3,4)\left[\begin{array}{l }        \eta_1\\        \eta_2\\        \eta_3\\        \end{array } \right ] = \left[\begin{array}{lll }    1&0&0\\    0&0&1\\    0&1&0\\    \end{array } \right ] \left[\begin{array}{l }        \eta_1\\        \eta_2\\        \eta_3\\        \end{array}\right].\end{aligned}\ ] ]    the relative coordinates in gusev @xcite extend the tetrahedral coordinates eq .
[ es30 ] to @xmath142 particles
. see appendix b for their transforms under permutations .
we give now a first example of the polynomial state construction in the tetrahedral coordinates eq .
[ es30 ] for the lowest excitation @xmath143=[30 ] ,   \lambda=3 , \ : f=[211]$ ] with spin @xmath144 and the orbital young tableaus @xmath145 of eq .
[ es16 ] : @xmath146 from matrix elements of the coulomb interaction to the energy level structure .
------------------------------------------------------------------------------    with the polynomial states for a specific set of quantum numbers at hand , we proceed to evaluate the various parts of the four - electron hamiltonian in the presence of a magnetic field @xmath147 perpendicular to the @xmath148 plane , @xmath149 with the larmor frequency @xmath150 , the effective mass @xmath151 and the combined frequency @xmath152 due to an external confinement potential @xmath153 . in all following expressions we work within the symmetric gauge , which allows one to express the uniform magnetic field in terms of two - dimensional oscillator excitations @xcite . the resulting four particle hamiltonian is readily decomposed into _
cm _ and relative parts in the plane perpendicular to the magnetic field direction . without the coulomb part
it resembles the addition of 2d oscillator systems with frequencies @xmath154 .
the magnetic field induces an angular - momentum dependent energy shift and the non - interacting eigenenergies read @xmath155 where we followed fock s notation . in terms of the excitations
@xmath156 these quantum - numbers read @xmath157 the total energy is given by the addition of _ cm _ and relative parts @xmath158 for the evaluation of the coulomb interaction we express the states as polynomials in the jacobi relative vectors .
there is , independent of the @xmath2 coordinate , a direct linear matrix transformation from relative tetrahedral to jacobi vectors obtained from eqs .
[ es29b ] , [ es30 ] , which for short we denote as @xmath159 the first jacobi vector from eq .  [ es29b ] is the relative vector between particles 1,2 and according to eq .
[ es13e ] it suffices to compute the full coulomb interaction .    in the relative - coordinate polynomials
the dimensionless coulomb interaction reads @xmath160 where the factor @xmath161 arises from the number of pairs for the two - body operator in eq .
[ es13e ] .
the coulomb matrix element between two states expressed in the three jacobi vectors becomes @xmath162 , { \rm min}[m^+_1 , m^-_1 ] , |n^-_1 - n^+_1|),\end{aligned}\ ] ] with @xmath163    we construct for a given set of relative - coordinate excitations all relative - coordinate states @xmath164 and express each state @xmath165 as sum of polynomials in jacobi coordinates with exponents @xmath166 @xmath167 then the interaction matrix is given by @xmath168    the coulomb interaction does not connect states with different angular momenta and so the blocks of states @xmath169 fall into subblocks labelled by @xmath170 .
however , different excitations of the relative oscillator are coupled by the coulomb interaction .
the addition of the relative oscillators , including the coulomb interaction is given by diagonalization of the matrix @xmath171 if we augment the relative eigenvalues with the analytic results for the _ cm _ energies , we are in a position to compare the results with the configuration - interaction methods using slater determinantes of single - particle states .
.list of 34 saci coulomb eigenvalues ( @xmath173 ) in units of @xmath174 for @xmath172 obtained by diagonalization of the slater determinants of the 34 single - particle product states @xmath175 yielding @xmath176 within the lowest landau level .
the entries in boldface indicate the seven states with the _ cm _ oscillator in the ground state obtained by the polynomial method described in the text for the lowest landau level ( @xmath177 ) , while the remaining 27 eigenvalues correspond to spurious states . [
cols="^,^,^,^ " , ]
the construction of the basis set of relative coordinate polynomials for @xmath0 electrons entailed the following steps : + ( 1 ) spin , partition : given the number @xmath0 of electrons , we determine the possible spins @xmath11 , orbital partitions @xmath15 and angular momenta @xmath81 .
the state space falls into independent and orthogonal subspaces labeled by @xmath178 .
steps ( 1 - 7 ) work for n=4 . for @xmath142 and steps ( 4 - 7 )
an extended analysis needs to be developed , some results are given in the appendices .
\(2 ) unitary quantum numbers : we provide the orbital quantum numbers of the system . for a maximal relative oscillator excitation
@xmath75 we list the unitary partitions @xmath79 , h_1+h_2=n$ ] for the groups @xmath179 .
a restriction to the lowest landau level @xmath86 with @xmath79 \rightarrow [ n0]$ ] greatly simplifies the analysis .
\(3 ) angular momenta : the unitary partition @xmath180 $ ] from @xmath181 fixes the ladder of angular momenta @xmath81 from the maximum value @xmath182 in steps of @xmath115 down to @xmath84 or @xmath83 and provides the final blocks @xmath178 .
\(4 ) pseudo angular momentum and partition : the pseudo angular momentum @xmath51 is obtained from @xmath180 $ ] by the representation theory of the rotation subgroup in @xmath183 . for @xmath4 , any fixed pseudo angular momentum @xmath51 yields from @xcite p. 260
a precise list of possible orbital partitions @xmath15 .
multiple values of pseudo angular momentum yield linearly independent polynomials .
\(5 ) state polynomials : given this list of quantum numbers , we construct state polynomials in relative coordinates .
the initial polynomial is known from @xcite .
the component of pseudo angular momentum and of total angular momentum are lowered independently .
\(6 ) young tableau : for the known orbital partitions @xmath15 , by applying permutation operators we project states with definite young tableau .
\(7 ) coulomb interaction and energies : for state polynomials in a block @xmath178 , we isolate the first jacobi vector in bra and a ket states , evaluate the matrix elements of the coulomb interaction , determine the hamiltonian matrix , diagonalize it and find the eigenstates of relative motion .
the study of the four - electron quantum - dot with coulomb interactions in the relative - coordinate basis elucidates and overcomes some of the problems inherent in the commonly used configuration interaction methods such as saci @xcite : if the @xmath2 excitations , as in determinantal state space , can not be separated , a large cloud of spurious @xmath2 excited states , augmented by spurious angular momenta , intrudes the variational state space and obscures the variational energy spectrum reflecting the coulomb interaction . within saci , higher
eigenenergies are distorted by the coulomb interaction and a loss in accuracy is seen , which can be avoided by the separation of the relative and @xmath2 coordinates .
this allows one to construct polynomial states with labels from representations of unitary and rotation groups and to diagonalize the coulomb - interaction part in smaller subblocks .
our approach opens a clear view on the energy spectrum as shown in figs . 1 and 2 and preserves the exact ladder spacing @xmath184 of _ cm _ copies of the relative coordinates variational subspace .
by increasing the basis size with respect to the relative excitations , the eigenenergies and eigenvectors obtained by diagonalization converge toward the exact eigenstates of the system .
tk has been supported by the heisenberg program of the dfg ( kr 2889/5 ) .
an independent view on orbital basis states for @xmath112 electons under permutations arises by use of the subgroup @xmath186 from @xcite .
this method can be extended to higher excitations .
it allows in tables  [ table8],[table9 ] to identify states of permutational symmetry for any pseudo angular momentum almost without projection .
we use @xmath187 to denote a standard state of pseudo angular momentum @xmath188 with component @xmath189 , @xmath190 for parity , and @xmath191 with the states defined in eq .
as mentioned in @xcite p. 263
, the group @xmath192 is the largest subgroup of @xmath111 that leaves invariant the double - dot @xmath193 axis up to a reflection . with the frobenius multiplicity expression we assign from the characters of @xmath192 , table  [ table6 ] , in table  [ table7 ] the reduction of ir f of @xmath111 to those of @xmath192 .
the numbers @xmath194 are taken from @xcite and used in table  [ table6 ] .
they determine the irreducible bases of representation of @xmath192 but from table  [ table7 ] also basis states of permutational symmetry .
the basis states @xmath118_3 $ ] and @xmath195_3 $ ] are unique , and so their partner states may be derived by the ladder procedure [ es14 ] .
the basis functions for @xmath196,[211]$ ] here transform like the doubledot representation . by @xcite eqs .
( 6.1 , 6.9 ) we can always return from the doubledot basis to the jacobi basis .
the basis states of @xmath132 $ ] can be separated from the bases of @xmath197 $ ] and @xmath198 $ ] by projection .
we fix for the moment the labels @xmath199 and introduce @xmath191 in eq .
motivated by @xcite , but independent of it , we redefine two new linear combinations of spherical harmonics @xcite ( 2.5.7 ) by @xmath200,\ : \mu \geq 0,\\   \nonumber   & & \rho=-1 : |-1(\lambda,\mu)\rangle = i^{-1}g_{\mu}[|(\lambda,\mu)\rangle-|(\lambda,-\mu)\rangle ( -)^{\lambda+\mu } ] , \ :   \mu>0 . \end{aligned}\ ] ] these states by consulting @xcite table 3 may be shown to be orthogonal eigenstates of the permutation operator @xmath201 with eigenvalue @xmath202 , which from @xcite fig 1 reverses the axis @xmath193 .
states characterized by @xmath203 we show to belong to irreps of the subgroup @xmath192 . for the proof
we employ the class operators @xcite , which are group operators @xmath204 summed over a class of equivalent group elements , and commute with any group operator . the five class operators @xmath205 @xcite of the group @xmath192 , written as sums of permutation operators from @xcite , table  4 are @xmath206 by use of the representation matrices of permutations we express the class operators eq .
[ a2 ] in matrix form , table [ table12 ] .
the matrix forms depend on the chosen representation , and we make convenient choices : for the representations @xmath207}$ ] and @xmath208 } , d^{[211]}(i , j)= i d^{[31]}((i , j)$ ] we use the basis eq .
[ es30 ] of @xmath209 , for the representation @xmath210}$ ] the young representation @xcite p. 226 . by virtue of these choices
we find : all class operators eq .  [ a2 ] have diagonal representations , see table [ table12 ] .
this proves : + * prop * : all basis states of the irreps of the subgroup @xmath192 in the chosen representations are basis states of irreducible representations of the bigger group @xmath111 .
+ this remarkable result allows in table  [ table7 ] , up to certain ambiguities , to almost avoid the use and projection with young operators for the bigger group @xmath111 . in tables
[ table8],[table9 ] we use it in relation with the full scheme of groups including @xmath211 and subgroups to assign orbital symmetry to the oscillator states .
@xmath212    @xmath213}(ki ) & d^{[31]}(ki ) & d^{[211]}(ki ) & d^{[22]}(ki ) & d^{[1 ^ 4]}(ki)\\ \hline k1      & \begin{array}{|l| } \hline      1\\ \hline    \end{array }    & \begin{array}{|l|l|l| }   \hline             1 &   &   \\   \hline                          & 1 &
\\    \hline             &   & 1\\    \hline          \end{array } &   \begin{array}{|l|l|l| }   \hline                   1 &   &   \\    \hline                        & 1 & \\     \hline             &   & 1\\     \hline      \end{array } & \begin{array}{|l|l| } \hline             1 &    \\ \hline             & 1\\ \hline      \end{array }   & \begin{array}{|l| } \hline      1\\   \hline    \end{array }   \\ \hline     k2      & \begin{array}{|l| } \hline      1\\ \hline    \end{array }    & \begin{array}{|l|l|l| }   \hline            -1 &   &   \\   \hline                          & -1 &
\\    \hline             &   & 1\\    \hline          \end{array } &   \begin{array}{|l|l|l| }   \hline                  -1 &   &   \\    \hline                        & -1 &
\\     \hline             &   & 1\\     \hline      \end{array } & \begin{array}{|l|l| } \hline             1 &    \\
\hline             & 1 \\ \hline      \end{array } & \begin{array}{|l| } \hline      1\\ \hline    \end{array}\\ \hline    k3      & \begin{array}{|l| } \hline      2\\ \hline    \end{array }    & \begin{array}{|l|l|l| }   \hline              &   &   \\   \hline                          &   & \\    \hline             &   & -2\\    \hline          \end{array } &   \begin{array}{|l|l|l| }   \hline                    &   &   \\    \hline                        &   & \\     \hline             &   & 2\\     \hline      \end{array } & \begin{array}{|l|l| } \hline             1 &    \\ \hline             & -1\\ \hline      \end{array }    &
\begin{array}{|l| } \hline      -2\\ \hline    \end{array}\\ \hline      k4 & \begin{array}{|l| } \hline      2\\ \hline    \end{array }    & \begin{array}{|l|l|l| }   \hline              &   &   \\   \hline                          &   & \\    \hline             &   & -2\\    \hline          \end{array } &   \begin{array}{|l|l|l| }   \hline                    &   &
\\    \hline                        &   & \\     \hline             &   & -2\\     \hline      \end{array } & \begin{array}{|l|l| } \hline             1 &
\\ \hline             & 1\\ \hline      \end{array }    & \begin{array}{|l| } \hline      2\\ \hline    \end{array}\\ \hline    k5 & \begin{array}{|l| } \hline      2\\ \hline    \end{array }    & \begin{array}{|l|l|l| }   \hline              &   &   \\   \hline                          &   & \\
\hline             &   & 2\\    \hline          \end{array } &   \begin{array}{|l|l|l| }   \hline                    &   &   \\    \hline                        &   & \\     \hline             &   & -2\\     \hline      \end{array } & \begin{array}{|l|l| } \hline             2 &
\\ \hline             & -2\\ \hline      \end{array }   & \begin{array}{|l| } \hline      -2\\ \hline    \end{array}\\    \hline   \end{array }    $ ]    @xmath214_2 | [ 4 ]     & [ 211]_3           \\ \hline   0       & -1             & [ 22]_1   |[1 ^ 4 ]   & [ 31]_3            \\
\hline   0       & 0              & [ 31]_1 | [ 211]_1&[31]_2 | [ 211]_2    \\ \hline 1        & + 1             & [ 22]_1 | [ 1 ^ 4 ]   & [ 31]_3            \\ \hline 1        & -1             & [ 22]_2 | [ 4 ]     & [ 211]_3           \\ \hline 1        & 0              & [ 211]_1| [ 31]_1&[211]_2 | [ 31]_2     \\
\hline \end{array } $ ]    of two states separated by @xmath215 , one and only one can belong to the listed tableau .
the states @xmath216,[1 ^ 4]$ ] are identified as eigenstates under the transposition @xmath217 with eigenvalue @xmath218 respectively .
if a state is not reproduced under @xmath217 , it necessarily belongs to @xmath132 $ ] and spin @xmath6 .
we conclude that the states eq .  [ a1 ]
yield all the bases of the orbital young tableaus .
@xmath219 +    & -           & ( 0,0 )            & 0^+           & [ 4]\\ \hline 1 ^ -&1 ^ -            & 1     & 1         & 0        & [ 31]1 +   & [ 31]2 +       & ( 1,0 )            & 1 ^ -            & [ 31]\\ \hline     &               & 0     & 1         & 1        & -        & [ 31]3 +      &                 &               & \\    \hline    2^+&2^+            & 2     & 0         & -1       & [ 22]1 +   & [ 31]3 +       & ( 2,0 )            & 2^+            & [ 31][22]\\ \hline     &               & 1     & 0         & 0        & [ 31]1 +   & [ 31]2 +       &                 &               &        \\
\hline     &               & 0     & 0         & 1        & [ 22]2 +   & -            &                 &               &         \\   \hline        & 1^+            & 1     & 0         & 0        & [ 211]1 + & [ 211]2 +      & ( 0,1 )            & 1^+            & [ 211 ]   \\
\hline     &               & 0     & 0         & 1        & -            & [ 211]3 + &                 &               &   \\ \hline     & 0^+            & 0     & 0         & 1        & [ 4]+     & -           & ( 2,0 )            & 0^+            & [ 4 ]   \\ \hline
3 ^ -&3 ^ -            & 3     & 1         & 0        & [ 211]1|[31]1*&[211]2|[31]2*&(3,0 )          & 3 ^ -            & [ 211][31][4]\\ \hline      &               & 2     & 1         & -1       & [ 4]+     & [ 211]3 +      &                 &               & \\
\hline     &               & 1     & 1         & 0        & [ 211]1|[31]1*&[211]2|[31]2 * &               &               & \\ \hline     &               & 0     & 1         & 1        & -            & [ 31]3 +   &                 &               &   \\ \hline     & 2 ^ -            & 2     & 1         & -1       & [ 22]2 +   & [ 211]3 +      & ( 1,1 )            & 2 ^ -            & [ 211][22 ] \\
\hline     &               & 1     & 1         & 0        & [ 211]1 + & [ 211]2 +      &                 &               &
\\ \hline     &               & 0     & 1         & 1        & [ 22]1 +   & -            &                 &               &   \\ \hline     & 1 ^ -            & 1     & 1         & 0        & [ 31]1 +   & [ 31]2 +       & ( 3,0),(1,1 )      & ( 1 ^ -)^2        & [ 31]^2
\\ \hline     &               & 0     & 1         & 1        & -        & [ 31]3 +       &                 &               &    \\
\hline 4^+&4^+            & 4     & 0         & 1        & [ 22]2|[4]*&[211]3 *      & ( 4,0 )            & 4^+            & [ 4][31][22][211]\\ \hline     &               & 3     & 0         & 0        & [ 31]1|[211]1*&[31]|2[211]2 * &                 &               &    \\
\hline     &               & 2     & 0         & -1       & [ 22]1 +   & [ 31]3 + &                 &               &    \\
\hline     &               & 1     & 0         & 0        & [ 31]1|[211]1*&[31]2|[211]2 * &                 &               &    \\ \hline     &               & 0     & 0         & 1        & [ 22]2|[4 ] *    & [ 211]3 *     &                 &               &
\\ \hline     & 3^+            & 3     & 0         & 0        & [ 31]1|[211]1*&[31]2|[211]2*&(2,1 )            & 3^+            & [ 31][211][1 ^ 4]\\ \hline     &               & 2     & 0         & -1       & [ 1 ^ 4 ]        & [ 31]3 +   &                 &               &      \\
\hline     &               & 1     & 0         & 0        & [ 31]1|[211]1*&[31]2|[211]2 * &                 &               & \\ \hline     &               & 0     & 0         & 1        & -           & [ 211]3 + &                 &               & \\
\hline     & 2^+            & 2     & 0         & -1       & [ 22]1 +   & [ 31]3 +       & ( 4,0),(2,1),(0,2)&(2^+)^3           & [ 31]^3[22]^3 \\
\hline     &               & 1     & 0         & 0        & [ 31]1 +   & [ 31]2 +       &                 &               &
\\ \hline     &               & 0     & 0         & 1        & [ 22]2 +   & -            &                 &               & \\ \hline     & 1^+            & 1     & 0         & 0        & [ 211]1+&[211]2 +       & ( 2,1 )           & 1^+            & [ 211]\\   \hline     &               & 0     & 0         & 1        & -        & [ 211]3 +      &                 &               & \\ \hline     & 0^+            & 0     & 0         & 1        & [ 4]+     & -            & ( 4,0),(0,2 )      & ( 0^+)^2            & [ 4]^2\\ \hline   \end{array } $ ]    @xmath2201|[31]1*&[211]2|[31]2 * & ( 5,0 )            & 5 ^ -            & [ 211][22][31]^2\\ \hline           &                   & 4           & 1         & 1        & [ 22]1 +        & [ 31]3 +      &                  &               &
\\ \hline     &                   & 3           & 1         & 0        & [ 211]1|[31]1*&[211]2|[31]2 * &                &               & \\ \hline     &                   & 2           & 1         & -1       & [ 22]2 +        & [ 211]3 +     &                  &               &
\\ \hline     &                   & 1           & 1         & 0        & [ 211]1|[31]1*&[211]2|[31]2 * &                &               &   \\ \hline     &                   & 0           & 1         & 1        & [ 22]1 +        & [ 31]3 +      &                  &               &
\\ \hline     & 4 ^ -                & 4           & 1         & 1        & [ 1 ^ 4]|[22]1 * & [ 31]3 +      & ( 3,1 )             & 4 ^ -            & [ 1 ^ 4][211][22][31]\\ \hline     &                   & 3           & 1         & 0        & [ 211]1|[31]1*&[211]2|[31]2 * &                &               &   \\
\hline     &                   & 2           & 1         & -1       & [ 22]2 +        & [ 211]3 +      &                 &               &   \\
\hline     &                   & 1           & 1         & 0        & [ 211]1|[31]1*&[211]2|[31]2 * &                &               &   \\ \hline     &                   & 0           & 1         & 1        & [ 22]1|[1 ^ 4 ] * & [ 31]3 +       &                 &               &   \\
\hline     & 3 ^ -                & 3           & 1         & 0        & [ 211]1|[31]1*&[211]2|[31]2*&(5,0)(3,1)(1,2)&(3 ^ -)^3            & [ 211]^3[31]^3[4]^3 \\ \hline       &                   & 2           & 1         & -1       & [ 4]+          & [ 211]3 +     &                  &               &   \\
\hline     &                   & 1           & 1         & 0        & [ 211]1|[31]1*&[211]2|[31]2 * &                &               &    \\ \hline     &                   & 0           & 1         & 1        & -             & [ 31]3 +      &                  &               &    \\
\hline     & 2 ^ -                & 2           & 1         & -1       & [ 22]2 +        & [ 211]3 +     & ( 3,1)(1,2 )        & ( 2 ^ -)^2            & [ 211]^2[22]^2   \\ \hline     &                   & 1           & 1         & 0        & [ 211]1 +       & [ 211]2 +     &                  &               &
\\ \hline     &                   & 0           & 1         & 1        & [ 22]1 +        & -          &                  &               &     \\
\hline     & 1 ^ -                & 1           & 1         & 0        & [ 31]1 +        &   [ 31]2 +    & ( 5,0)(3,1)(1,2 )   & ( 1 ^ -)^3        & [ 31]^3 \\
\hline       &                   & 0           & 1         & 1        & -            & [ 31]3 +       &                  &               &   \\ \hline \end{array }     $ ]
the efficiency of the tetrahedral coordinates raises the question if similar relative coordinates exist for @xmath142 .
as a generalization of the tetrahedral coordinates from @xcite , new symmetrized coordinates for @xmath0 particles were proposed by gusev et al @xcite .
the matrix that gives the @xmath0 new coordinates @xmath221 in terms of the old ones @xmath222 reads @xmath223 ,   \\
\nonumber   & & c = c^t,\ : c^{-1}=c,\ : c^2=i ,   \\
\nonumber   & & a=[1-\sqrt{n}]^{-1},\ : b = a+ \sqrt{n}.   \end{aligned}\ ] ] the _ cm _ coordinate is included as @xmath224 .
we shall explore the properties of these coordinates under the action of permutations .
all permutations are generated by products of the @xmath225 transpositions @xmath226 .
we need the action of these transpositions on the new relative coordinates .
we shall see in eq .
[ b4 ] that under permutations @xmath227 acting on particles @xmath228 the symmetrized coordinates transform like single particle coordinates .
so we focus on the remaining action @xmath229 of the first transposition @xmath230 with respect to the new coordinates .      for @xmath231
we obtain for the linear action of the transposition @xmath230 on the new coordinates @xmath232   \left [   \begin{array}{lllll }   0&1&0&0&0\\   1&0&0&0&0\\   0&0&1&0&0\\   0&0&0&1&0\\   0&0&0&0&1\\   \end{array }   \right ]    \left [   \begin{array}{lllll }   1&1&1&1&1\\   1&b&a&a&a\\   1&a&b&a&a\\
1&a&a&b&a\\   1&a&a&a&b\\   \end{array }   \right ]    \\ \nonumber & & = \frac{1}{5 } \left [   \begin{array}{l|l|l|l|l } \hline   5      & 1 + 3a+b       & 1 + 3a+b      & 1 + 3a+b      & 1 + 3a+b\\ \hline   1 + 3a+b&2b+3a^2      & a+b+2a^2+ab&a+b+2a^2+ab&a+b+2a^2+b^2\\ \hline   1 + 3a+b&a+b+2a^2+ab & 2a+2a^2+b^2&2a+a^2 + 2ab &
2a+a^2 + 2ab\\ \hline   1 + 3a+b&a+b+2a^2+ab & 2a+a^2 + 2ab & 2a+2a^2+b^2&2a+a^2 + 2ab\\ \hline   1 + 3a+b&a+b+2a^2+b^2&2a+a^2 + 2ab & 2a+a^2 + 2ab & 2a+2a^2+b^2\\ \hline   \end{array }   \right ] \end{aligned}\ ] ] this matrix is symmetric and orthogonal by construction , its square is the unit matrix .
we compute from eqs.[es21],[es26],[es30 ] the numbers @xmath233,\ : b= -\frac{1}{4}[1 - 3\sqrt{5 } ] , \\
\nonumber & & ( 002300)=2b+3a^2=\frac{1}{8}[5-\sqrt{5 } ] , \\
\nonumber & & ( 011201)=a+b+2a^2+ab=\frac{5}{8}[-1+\sqrt{5 } ] , \\
\nonumber & & ( 020210)=2a+2a^2+b^2=\frac{5}{8}[5-\sqrt{5 } ] , \\
\nonumber   & & ( 020102)=2a+a^2 + 2ab=-\frac{5}{8}[3+\sqrt{5}].\end{aligned}\ ] ] for @xmath231 , all the numbers and matrices can be decomposed into a rational part and a second part rational but proportional to @xmath234 .
all these numbers form a module that closes under addition , multiplication and division , compare the example eq .
we give this decomposition in the following examples .
a similar module decomposition applies for any @xmath0 whose square root is not an integer .
@xmath235 , b=\frac{1}{4\cdot 5 } \left [ \begin{array}{rrrrr } 4      & 4 & 4 & 4 & 4\\ 4      & -1&-1&-1&-1\\ 4      & -1&-1&-1&-1\\ 4      & -1&-1&-1&-1\\ 4      & -1 & -1&-1 & -1\\ \end{array } \right ] , \\
\nonumber & & d^{\eta}(1,2)= a^{\eta}+ \sqrt{5}b^{\eta } ,   \\ \nonumber & & a^{\eta}= \frac{1}{5\cdot8 } \left [ \begin{array}{rrrrr } 5\cdot 8&0 & 0 & 0 & 0\\ 0        & 5      & -5&-5 & -5\\ 0        & -5&25&-15&-15\\ 0        & -5&-15&25&-15\\ 0        & -5&-15&-15&25\\ \end{array } \right ] , \ : b^{\eta}=   \frac{1}{5\cdot 8 } \left [ \begin{array}{rrrrr } 0       & 0    & 0    & 0    & 0\\ 0       & -1    & 5 & 5   & 5\\ 0       & 5    & -5   & -5   & -5\\ 0       & 5    & -5   & -5   & -5\\ 0       & 5    & -5   & -5   & -5\\ \end{array } \right].\end{aligned}\ ] ]    we find with eq.[b3 ] from the values of @xmath236 for @xmath231 , @xmath237 and so the product matrix eq.[b2 ] reduces to its diagonal @xmath238 submatrix that involves only the relative coordinates @xmath239 .
this is always necessary , since the transposition @xmath230 can not affect the _ cm _ coordinate @xmath224 . for the other transpositions
@xmath240 we find the result @xmath241 so all transpositions @xmath56 in single particle coordinates @xmath242 act on the relative coordinates @xmath243 like transpositions of relative coordinates which we denote as @xmath244 .
we summarize eq .
[ b4 ] in symbolic form as @xmath245 the transposition @xmath246 , missing in eq .
[ b7 ] , and acting on the symmetrized coordinates , is given by the matrix @xmath229 eq .
[ b10 ] .
[ [ permutations - of - symmetrized - relative - coordinates - for - n - electrons . ] ] permutations of symmetrized relative coordinates for @xmath0 electrons .
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~    the results of the last subsection , in particular eq .
[ b7 ] generalize to any @xmath0 . since the transpositions eq .  [ b4 ] generate the subgroup @xmath247 , all permutations of particles @xmath248 are mapped by eq .
[ b7 ] into permutations of the relative coordinates @xmath249 . the correspondence eq .
[ b7 ] must be kept in mind when working with the new relative coordinates .
we could easily antisymmetrize a polynomial state wrt to the ( n-1 ) particles @xmath250 , but must augment the antisymmetrizer @xmath251 to include particle @xmath83 , see the next section . in terms of the symmetric group , we have the subgroup @xmath247 acting on particles @xmath250 . all elements @xmath252 can be written as @xmath253 with @xmath254 , multiplied by the @xmath0 coset generators which may be chosen as transpositions @xmath255 .
these coset generators can be rewritten as @xmath256 and so their matrices @xmath257 acting on relative coordinates reduce to the element @xmath229 whose matrix we compute in eq .
[ b7 ] . the representations @xmath258 appearing in eq .
[ b5 ] simply mean the interchange of two rows or two columns .
so our main task is to give the generalization of eq .
[ b2 ] for any @xmath0 .
we denote the matrix elements of @xmath229 for short by six integers @xmath259^{-1 } , b = a+ \sqrt{n},\end{aligned}\ ] ] and obtain for the matrix representation in @xmath260-coordinates and general @xmath0 the result @xmath261   \left [   \begin{array}{llllll }   0&1&0&0&0&.\\   1&0&0&0&0&.\\   0&0&1&0&0&.\\   0&0&0&1&0&.\\   0&0&0&0&1&.\\   .&.&.&.&.&.\\
\end{array }   \right ]    \left [   \begin{array}{llllll }   1&1&1&1&1&.\\   1&b&a&a&a&.\\   1&a&b&a&a&.\\   1&a&a&b&a&.\\   1&a&a&a&b&.\\   .&.&.&.&.&.\\
\end{array }   \right ]    \\ \nonumber & & = \frac{1}{n } \left [   \begin{array}{l|l|l|l|l|l } \hline   ( n00000 ) & ( 1n\!-\!21000 ) & ( 1n\!-\!21000 ) & ( 1n\!-\!21000 ) & ( 1n\!-\!21000 )   & .\\
\hline   ( 1n\!-\!21000)&(002n\!-\!200)&(011n\!-\!301)&(011n\!-\!301)&(011n\!-\!301)&.\\ \hline   ( 1n\!-\!21000)&(011n\!-\!301)&(020n\!-\!310)&(020n\!-\!402)&(020n\!-\!402)&.\\ \hline   ( 1n\!-\!21000)&(011n\!-\!301)&(020n\!-\!402)&(020n\!-\!310)&(020n\!-\!402)&.\\ \hline   ( 1n\!-\!21000)&(011n\!-\!301)&(020n\!-\!402)&(020n\!-\!402)&(020n\!-\!310)&.\\ \hline   ..              & ..             & ..             & ..             & ..              & .\\
\hline   \end{array } \right]\end{aligned}\ ] ] the matrix elements are given in the notation eq .
[ b9 ] . from eq .
[ b10 ] the full matrix @xmath229 is constructed for any @xmath0 : the numbers @xmath236 are inserted from eq .
[ b1 ] as functions of @xmath0 .
the diagonal elements indexed by @xmath262 are repeated along the diagonal , the elements in the lines on top , and columns below the diagonal , are repeated in rows to the right and in columns downwards respectively by the same functions eq .
[ b9 ] of @xmath0 .
again @xmath229 is orthogonal and symmetric , its square is the unit matrix .
we find in general from eq .
[ b1 ] for the entries of the first row and column of @xmath229 @xmath263 and so the @xmath264 matrix @xmath229 reduces to its @xmath265 diagonal submatrix acting exclusively on the relative coordinates @xmath243 . in building orbital states from the symmetrized gusev coordinates we must take care in the application of permutation operators .
99 aguilera - navarro v c , kramer p , moshinsky m + _ harmonic oscillator states and the @xmath266-particle ii .
configuration - space states of arbitrary symmetry .
_ ann of physics ( n y ) * 54 * ( 1969 ) 379393 avron j e , herbst i w , and simon b + _ separation of center of mass in homogeneous magnetic fields _ annals of physics * 114 * ( 1978 ) 431451 bargmann v and moshinsky m + _ group theory of harmonic oscillators ( i ) . _ + nuclear physics * 18 * ( 1960 ) 697712 bargmann v and moshinsky m + _ group theory of harmonic oscillators ( ii ) . _ + nuclear physics * 23 * ( 1961 ) 177199 barut a o and raczka r + _ theory of group representations and applications . _ + pwn polish scientific publishers , warsaw ( 1977 ) edmonds a r + _ angular momentum in quantum mechanics . _ + princeton university press .
princeton ( 1957 ) escobar - ruiz m a and turbiner a v + _ two charges on plane in a magnetic field i. `` quasi - equal '' charges and neutral quantum system at rest cases . _ + ann phys * 340 * ( 2014 ) 3759 grossmann f and kramer t + _ spectra of harmonium in a magnetic field using an initial value representation of the semiclassical propagator .
_ + j phys a * 44 * ( 2011 ) 445309 gusev a et al .
+ _ symbolic - numerical algorithm for generating cluster eigenfunctions : identical particles with pair oscillator interactions .
_ in : computer algebra in scientific computing , springer , berlin ( 2013 ) 155168 hamermesh m + _ group theory and its application to physical problems . _ + addison - wesley , reading ( 1962 ) 182278 johnson b r , hirschfelder j o , and yang k - h + _ interaction of atoms , molecules , and ions with constant electric and magnetic fields .
_ rev mod phys * 55 * ( 1983 ) 109153 klama s and mishchenko e + _ two electrons in a quantum dot : a semiclassical approach . _ + j.  phys . : condens .
matter * 10 * ( 1998 ) 34113416 kohn w + _ cyclotron resonance and de haas - van alphen oscillations of an interacting electron gas . _ + phys rev * 123 * ( 1961 ) 12421244 kramer p and moshinsky m + _ group theory of harmonic oscillators ( iii ) : states with permutational symmetry .
_ nuclear physics * 82 * ( 1966 ) 241273 kramer p , john g and schenzle d + _ group theory and the interaction of composite nucleon systems . _ + vieweg , braunschweig ( 1971 ) kramer p and moshinsky m + _ group theory of harmonic oscillators ( iv ) : dynamics of the 2,3 , and 4 particle systems .
_ nuclear physics * a 107 * ( 1968 ) 481522 kramer p and moshinsky m + _ group theory of harmonic oscillators and nuclear structure .
_ in : group theory and its application , + ed e m loebl , academic press , new york ( 1968 ) 340368 kramer p + _ permutation group in light nuclei . in : symmetry properties of nuclei , proc 15th solvay conf
. on physics 1970
_ , gordon and breach , new york ( 1974 ) 211232 kramer p , kramer t , and manko v i + _ propagation in time of some classical / quantum systems and nielsen s conserved quantity . _ + phys scr * 79 * ( 2009 ) 055006 kramer t + _ two interacting electrons in a magnetic field : comparison of semiclassical , quantum , and variational solutions . _
+ aip conference proceedings * 1323 * ( 2010 ) 178 kramer t , krueckl v , heller e j and parrott r e , + _ self - consistent calculation of electric potentials in hall devices . _
+ phys rev b * 81 * ( 2010 ) 205306 kramer t , kreisbeck c , krueckl v , heller e j , parrott r e and liang c - t + _ theory of the quantum hall effect in finite graphene devices . _ + phys rev b * 81 * ( 2010 ) 081410(r ) kramer t + _ interacting electrons in a magnetic field : mapping quantum mechanics to a classical ersatz - system .
_ + aip conference proceedings * 1488 * ( 2012 ) 394 laughlin r + _ anomalous quantum hall effect : an incompressible quantum fluid with fractionally charged excitations . _
+ phys rev lett * 50 * ( 1983 ) 13951398 moshinsky m and quesne c + _ oscillator systems . in : symmetry properties of nuclei ,
proc 15th solvay conf . on physics 1970 .
_ gordon and breach , new york ( 1974 ) 233300 moshinsky m and smirnov y + _ the harmonic oscillator in modern physics .
_ + harwood academic publishers , amsterdam ( 1996 ) muhandiramge r and wang j + _ electronic structure of multi - electron quantum dots . _ + mathematica journal * 10:2 * ( 2006 ) 396414 muhandiramge r d and wang j + _ spinadapted states : a basis for quantum dot structure calculation . _ international journal of quantum chemistry * 106.1 * ( 2006 ) 2753 su b , goldman v j , and cunningham j e + _ evidence for spin singlet - triplet transitions of two electrons in a quantum dot observed via single - electron tunneling _ surf sci * 305 * ( 1994 ) 566 - 570 taut m + _ two electrons in a homogeneous magnetic field : particular analytical solutions . _
+ j phys a : math gen * 27 * ( 1994 ) 47234724 tavernier m b et al .
+ _ four - electron quantum dot in a magnetic field . _ + phys .
rev b * 68 * ( 2003 ) 205305 simonovi n s + _ three electrons in a two - dimensional parabolic trap : the relative motion solution .
_ + few - body systems * 38 * ( 2006 ) 139145 verin a + _ two anyons in a static , uniform magnetic field
. exact solution . _ + physics letters b * 260 * ( 1991 ) , 120124 wagner m , merkt u , and chaplik a v + _ spin - singlet spin - triplet oscillations in quantum dots . _
+ phys rev b * 45 * ( 1992 ) 19511954 weyl h + _ the theory of groups and quantum mechanics .
_ dover , new york ( 1930 ) wigner e p + _ group theory and its application to the quantum mechanics of atomic spectra . _ academic press , new york ( 1959 ) yannouleas c and landman u + _ quantal molecular description and universal aspects of the spectra of bosons and fermions in the lowest landau level . _
a * 81 * ( 2010 ) 023609 yukich j n , kramer t , and bracher c + _ observed photodetachment in parallel electric and magnetic fields .
a * 68 * ( 2003 ) 033412