recently , there has been a wide interest in relativistic potentials involving mixtures of vector and scalar potentials with opposite signs .
the interest lies on attempts to explain the pseudospin symmetry in nuclear physics .
@xcite , using a dirac hamiltonian with scalar @xmath1 and vector @xmath2 potentials quadratic in space coordinates , found a harmonic - oscillator - like second order equation which can be solved analytically for @xmath3 , as considered before by kukulin _
@xcite , and also for @xmath4 .
very recently , ginocchio solved the triaxial , axial , and spherical harmonic oscillators for the case @xmath5 and applied it to the study of antinucleons embedded in nuclei @xcite .
the case @xmath6 is particularly relevant in nuclear physics , since it is usually pointed out as a necessary condition for occurrence of pseudospin symmetry in nuclei @xcite-@xcite . also the observed isospin dependence of pseudospin
can be explained by the effect of the @xmath7 potential on the shape of the @xmath8 potential @xcite-@xcite . as shown by bell and ruegg @xcite
, the @xmath9 and @xmath6 cases correspond to a @xmath10 symmetry of the dirac hamiltonian with only scalar and vector potentials which is independent of the particular shapes of these potentials .    a generalized harmonic oscillator ( ho ) for spin 1/2 particles that includes not only those combinations of scalar and vector potentials , but also a linear tensor potential , giving rise to the so - called dirac oscillator , has also been considered in @xcite .
the research of this kind of interaction was started by it @xcite and has been revived by moshinsky and szczepaniak @xcite ( see @xcite for a complete reference list ) .    in ref .
@xcite a special attention was paid to the case when the scalar and vector potentials are such that @xmath11 , which , as stated before , has been related to the nuclear pseudospin symmetry .
there it has been concluded that only when the tensor potential is turned on one gets negative energy bound - state solutions along with the usual positive energy solutions .
those negative energy solutions are important to describe the possible existence of antinucleons in nuclei @xcite . when the tensor potential is absent the special conditions relating the scalar and vector potentials ( @xmath12 ) are needed to have a standard positive ho potential which is only able to bind fermions , and therefore exclude the negative bound - state solutions from the spectra
actually , as was discussed in another context in @xcite , with an appropriate choice of signs , potentials fulfilling the relations @xmath12 are able to bind either fermions or antifermions . as an example
, we will see later that a ( negative ) harmonic potential , with no tensor potential , binds antifermions and does not allow positive bound - state solutions .
closely related to this is the fact that , in the nucleus , the charge - conjugation transformation relates the spin symmetry of the negative bound - state solutions ( antinucleons ) to the pseudospin symmetry of the positive bound - state solutions ( nucleons ) @xcite .
this has also been discussed recently in @xcite , analyzing the ho for antinucleons with spin symmetry ( @xmath5 ) and without the tensor interaction .
again , the charge conjugation connects this problem with the ho for nucleons with pseudospin symmetry ( @xmath13 ) , but in this case the positive harmonic binding potential is @xmath14 .
note that , because the harmonic potential does not vanish for large distances but is a confining potential , it is possible to obtain particle bound - state solutions in the exact pseudospin limit @xcite .
therefore , we believe that this connection between spin and pseudospin symmetry obtained by charge conjugation deserves to be more explored . in addition , after the pseudospin symmetry has been revealed as a relativistic symmetry by ginocchio @xcite , a link of it to chiral symmetry has been suggested based on sum rules calculations in relativistic chiral models that support a small value for @xmath8 @xcite . however , until now , a clear connection of pseudospin and chiral symmetries has not been established .
the main motivation of this paper is to shed some light in the relation between spin and pseudospin symmetries by means of charge - conjugation and @xmath0 chiral transformations . in order to do so
, we solve the mixed scalar - vector - pseudoscalar ho potential in 1 + 1 dimensions , taking advantage of the simplicity of the lowest dimensionality of the space - time .
this approach is equivalent to consider fermions in 3 + 1 dimensions which are restricted to move in one direction @xcite .
we explore the spectra when it is possible to obtain analytical solutions , _
i.e. _ , in the particular cases when @xmath5 and @xmath6 . as referred before , if the pseudocalar term is turned off , these cases correspond , respectively , to have exact spin and pseudospin symmetry in 3 + 1 dimensions .
we explore all the possible signs of the quadratic ( @xmath15 or @xmath16 ) and linear ( pseudoscalar ) potentials , thus paying attention to bound states of fermions and antifermions as well .
we compare both cases @xmath9 and @xmath6 to establish the charge - conjugation connection discussed above in the presence of the pseudoscalar term .
we also consider the case of zero mass to look for the connection between spin and pseudospin symmetries by means of the chiral transformation .
an added motivation for this paper is given by the recent demonstration that the relativistic 3 + 1 ho with scalar and vector potentials can describe successfully the heavy nucleus spectrum @xcite .
the parameters of the ho are determined by fitting the scalar and vector potentials derived from relativistic mean field nuclear calculations ( rmf ) .
the @xmath17pb neutron single - particle levels obtained in the ho are very similar to the rmf results .
the pseudospin symmetry is shown to be almost satisfied for a heavy nucleus in the ho description .
the possibility of breaking perturbatively this symmetry in the 3 + 1 ho by a tensor interaction ( in our 1 + 1 case a pseudoscalar coupling ) has been discussed in @xcite and can be included in a calculation like the one of ref .
@xcite in order to improve the results .
thus , the study of all the possible eigenenergies of the 1 + 1 ho , presented in this paper , considering not only the positive energy solutions already obtained for the 3 + 1 case in @xcite but also the negative ones , can be applied to analyze the existence of antinucleons in nuclei and their spectrum .
the particle and antiparticle spectra we obtain for the cases @xmath6 and @xmath9 demonstrate the charge - conjugation connection referred above .
they also show that the spectra of massless @xmath5 states are degenerate with massless @xmath11 states . further imposing chiral symmetry to these states , such that they would become eigenstates of @xmath0 ,
would imply that all potentials would be zero and thus one would have just massless free fermions .
our conclusions are valid in the case of 3 + 1 dimensions , because , as will be shown later , one basically has to add the angular momentum quantum numbers @xmath18 and @xmath19 to the coefficients of the functions appearing in the 1 + 1 eigenvalue equations in order to get the corresponding 3 + 1 equations .
this paper is organized as follows . in sec .  [
sec : dirac_1 + 1 ] we present the general dirac equation in 1 + 1 dimensions with a potential with a completely general lorentz structure .
the effect of charge - conjugation and the @xmath0 chiral transformations upon the general 1 + 1 dirac hamiltonian is discussed in secs .
[ subsec : charge_conj ] and [ subsec : chiral ] respectively . in secs .
[ subsec : eq_motion ] and [ subsec : sturm - liouville ] we present the equations of motion and discuss the isolated solutions , _
i.e. _ , the solutions out of the sturm - liouville problem , for @xmath9 and @xmath6 . the nonrelativistic limits of the @xmath9 and @xmath6 cases are discussed in sec .
[ subsec : nonrel_limit ] , where we show that when @xmath6 there is no nonrelativistic limit when the pseudoscalar potential vanishes or is very small . in sec .
[ sec : rel_h_oscillator ] we introduce the general relativistic harmonic oscillator in 1 + 1 dimensions for @xmath9 by assigning a quadratic potential @xmath20 for @xmath8 and a linear potential @xmath21 for the pseudoscalar potential @xmath22 .
there we also present the eigenvalue equation and the wave functions for bound states , discuss the dirac oscillator case ( @xmath23 and @xmath24 ) and show how one can use the charge - conjugation and chiral transformations to get the eigenenergies for the @xmath6 case from the @xmath9 case .
the following section is devoted to a detailed analysis of the solutions of the @xmath9 relativistic ho for all signs of the parameters @xmath25 and @xmath26 .
finally in sec .
[ sec : sol_sym_arg ] we discuss the spectrum for @xmath6 and massless fermions and in sec .  [ sec : conclusions ] we draw the conclusions .
the two - dimensional dirac equation can be obtained from the four - dimensional one with the mixture of spherically symmetric scalar , vector and anomalous magnetic - like ( tensor ) interactions . if we limit the fermion motion to the @xmath27-direction ( @xmath28 ) the four - dimensional dirac equation decomposes into two equivalent two - dimensional equations with 2-component spinors and 2@xmath292 matrices @xcite
. then , there results that the scalar and vector interactions preserve their lorentz structures whereas the anomalous magnetic - like interaction turns out to be a pseudoscalar interaction .
furthermore , in the 1 + 1 world there is no angular momentum so that the spin is absent .
therefore , the 1 + 1 dimensional dirac equation allows us to explore the physical consequences of the negative - energy states in a mathematically simpler and more physically transparent way . in this spirit
the two - dimensional version of the anomalous magnetic - like interaction linear in the space coordinate has also received attention @xcite-@xcite .
later this system was shown to be a particular case of a more general class of exactly solvable problems asc2 .
let us begin by presenting the dirac equation in 1 + 1 dimensions .
the dirac equation for a particle with mass @xmath30 and a general potential @xmath31 , to be defined later , is formally the same as its 3 + 1 counterpart : @xmath32here @xmath33 and the @xmath34 matrices @xmath35 obey the usual relations @xmath36 , where @xmath37 is the @xmath34 identity matrix .
the metric tensor and the covariant derivative operator are , respectively , @xmath38 and @xmath39 .
the hamiltonian is defined in the usual way , such that eq .
( eq_dirac_gen ) is written as @xmath40where @xmath41 .
the positive definite function @xmath42 , satisfying a continuity equation , is interpreted as a position probability density and its norm is a constant of motion .
this interpretation is completely satisfactory for single - particle states @xcite .
the traceless matrices @xmath43 and @xmath44 are defined by @xmath45 and @xmath46 , and obey the relations @xmath47 , @xmath48 .
we set @xmath31 to be @xmath49where @xmath50 .
this is the most general combination of lorentz structures because there are only four linearly independent 2@xmath292 matrices .
the subscripts for the terms of the potential denote their properties under a lorentz transformation : @xmath51 and @xmath52 for the time and space components of the 2-vector potential , @xmath53 and @xmath54 for the scalar and pseudoscalar terms , respectively .
if the terms in the potential @xmath31 are time - independent , the dirac equation ( [ def_hamiltoniano_d ] ) becomes @xmath55 where @xmath56 and @xmath57 is the fermion energy . to have an explicit expression for the @xmath58 and @xmath59 matrices one can choose 2@xmath292 pauli matrices which satisfy the same algebra .
we use @xmath60 , @xmath61 and thus @xmath62 .
the hamiltonian ( [ time_ind_h ] ) is invariant under the parity operation , _
i.e. _ , when @xmath63 , @xmath64 and @xmath65 change sign , whereas @xmath66 and @xmath67 remain the same .
this is because the parity operator is @xmath68 , where @xmath69 is a constant phase and @xmath70 changes @xmath27 into @xmath71 . since this unitary operator anticommutes with @xmath58 and @xmath72 , they change sign under a parity transformation , whereas @xmath37 and @xmath59 , which commute with @xmath73 , remain the same .
when one writes down the explicit equations of motion in terms of the components of the spinor @xmath74 , the combinations @xmath75 and @xmath76 of the vector and scalar components arise naturally . therefore , it is convenient to rewrite the hamiltonian ( [ time_ind_h ] ) in terms of these potentials .
we have @xmath77      the charge - conjugation operation changes the sign of the electromagnetic interaction , _
i.e. _ , changes the sign of the time and space potentials in ( [ def_v ] ) .
this is accomplished by the transformation ( see , for example , itzykson and zuber @xcite ) @xmath78where @xmath79 denotes matrix transposition and @xmath80 is a matrix such that @xmath81 . in @xmath82 dimensions one matrix that satisfies
this relation is @xmath83if we choose the phase factor @xmath84 equal to 1 , we have @xmath85after applying this charge - conjugation operation to the dirac equation ( def_hamiltoniano_d ) , the time - independent dirac equation becomes @xmath86where @xmath87 and @xmath88 is given by @xmath89 in terms of the potentials @xmath15 and
@xmath16 , this hamiltonian reads @xmath90one sees that the charge - conjugation operation changes the sign of the energy and of the potentials @xmath91 , @xmath92 and @xmath93 . in turn , this means that @xmath16 turns into @xmath94 and @xmath15 into @xmath95 .
therefore , to be invariant under charge conjugation , the hamiltonian must contain only a scalar potential .
the chiral operator for a dirac spinor is the matrix @xmath0 , and we will call `` chiral transformation '' the transformation associated with it .
thus , the transformed spinor is given by @xmath96 and the transformed hamiltonian @xmath97 .
since @xmath98 anticommutes with @xmath59 , the time - independent chiral - transformed dirac equation is @xmath99 where @xmath100 is given by @xmath101 or @xmath102 in terms of @xmath8 and @xmath14 .
this means that the chiral transformation changes the sign of the mass and of the scalar and pseudoscalar potentials , thus turning @xmath8 into @xmath14 and vice - versa .
a chiral - invariant hamiltonian needs to have zero mass and @xmath103 and @xmath22 zero everywhere .
the space component of the 2-vector potential in ( [ def_v ] ) can be absorbed into the wave function by defining a new spinor @xmath104 such that @xmath105 in which @xmath106 , since we have @xmath107 . from this point on
, we will refer to @xmath91 as simply a vector potential , following the common usage of this term ( usually denoted by @xmath2 ) .
if we now write the spinor @xmath104 in terms of its components , @xmath108the dirac equation gives rise to two coupled first - order equations for @xmath109 and @xmath110 : @xmath111where the prime denotes differentiation with respect to @xmath27 . in terms of @xmath109 and @xmath110 the spinor is normalized as @xmath112 , so that @xmath109 and @xmath110 are square integrable functions .
it is clear from the pair of coupled first - order differential equations ( eq_mot_1 ) and ( [ eq_mot_2 ] ) that @xmath109 and @xmath110 have opposite parities if the dirac equation is covariant under @xmath63 .    under the charge - conjugation and chiral transformations the spinor ( def_psi_comp ) becomes @xmath113 and @xmath114 respectively .      using the expression for @xmath115 obtained from ( [ eq_mot_2 ] ) with @xmath116 , _ viz .
_ @xmath117 and inserting it in ( [ eq_mot_1 ] ) one arrives at the following second - order differential equation for @xmath118 : @xmath119 \psi_{+ } = 0\ . \label{2_order_gen_psi+}\ ] ]    in a similar way , using the expression for @xmath120 obtained from ( [ eq_mot_1 ] ) with @xmath121 , _ viz .
_ @xmath122 and inserting it in ( [ eq_mot_2 ] ) one arrives at the following second - order differential equation for @xmath123 : @xmath124 \psi_{- } = 0   \label{2_order_gen_psi-}\ ] ] for @xmath5 with @xmath125 , ( [ psi-_gen ] ) and ( 2_order_gen_psi+ ) reduce to @xmath126 @xmath127 \psi_{+ } = \frac{e^{2}-m^{2}c^{4}}{2mc^{2}}\psi_{+}\ , \label{2_order_delta=0_psi+}\ ] ] and for @xmath11 with @xmath128 , ( [ psi+_gen ] ) and ( [ 2_order_gen_psi- ] ) reduce to @xmath129 @xmath130 \psi_{- } = \frac{e^{2}-m^{2}c^{4}}{2mc^{2}}\psi_{-}\ . \label{2_order_sigma=0_psi-}\ ] ] either for @xmath5 with @xmath125 or @xmath11 with @xmath131 the solution of the relativistic problem is mapped into a sturm - liouville problem in such a way that solution can be found by solving a schrdinger - like problem . if one considers potentials @xmath91 and @xmath1 quadratic in @xmath27 with @xmath132 and a potential @xmath93 linear in @xmath27 one obtains schr'odinger - like equations for the harmonic oscillator potential .
the solutions for @xmath5 with @xmath133 and @xmath11 with @xmath134 , excluded from the sturm - liouville problem , can be obtained directly from the original first - order equations ( [ eq_mot_1 ] ) and ( eq_mot_2 ) .
they are @xmath135 \\ \\
\displaystyle\psi_{- } ^{\prime } + \frac{v_{p}}{\hbar c}\,\psi_{- } = -\frac{i}{% \hbar c}\left ( \sigma + 2mc^{2}\right ) \psi_{+}% \end{array } \right\ } , \mathrm{\;for\;}\delta = 0\mathrm{\;with\;}e =- mc^{2 } \label{isol_sol_delta=0}\ ] ] and @xmath136 \\ \\
\displaystyle\psi_{+ } ^{\prime } -\frac{v_{p}}{\hbar c}\,\psi_{+ } = -\frac{i}{% \hbar c}\left ( \delta -2mc^{2}\right ) \psi_{-}% \end{array } \right\ } , \mathrm{\;for\;}\sigma = 0\mathrm{\;with\;}e = mc^{2 } \label{isol_sol_sigma=0}\ ] ] where @xmath137 and @xmath138 are normalization constants .
since , by parity conservation , @xmath22 must be an odd function of @xmath139 , its integral must be an even function of @xmath27 .
thus , if the integral of @xmath22 has a definite sign , its overall sign can determine whether the isolated solutions ( [ isol_sol_delta=0 ] ) and ( [ isol_sol_sigma=0 ] ) are bound states .
in particular , if @xmath22 and @xmath140 or @xmath14 are powers of @xmath27 , one sees from each pair of equations ( [ isol_sol_delta=0 ] ) and ( [ isol_sol_sigma=0 ] ) that , in general , one of the spinor components has to be zero , either @xmath118 in the first case ( sign of the @xmath22 coefficient positive ) or @xmath115 in the last case ( sign of the @xmath22 coefficient negative ) .
this is the case of the harmonic oscillator potentials , as we will see in the next section .
however , it is possible find solutions with the @xmath22 coefficient negative in the first case ( @xmath5 , @xmath133 ) and with the @xmath22 coefficient positive in the second case ( @xmath11 , @xmath141 ) , provided a certain relation between the coefficients of all the potentials is satisfied , as we will see later .      in the nonrelativistic limit , @xmath142 , eq .
( [ 2_order_gen_psi+ ] ) becomes @xmath143+\sigma\psi_{+}=\mathcal{e}\psi_{+ } \ , \ ] ] where @xmath144 .
since , in the nonrelativistic regime , for the range of values of @xmath139 in which the wave function is not negligible , @xmath145 , from eq .
( [ psi-_gen ] ) we see that @xmath123 is of order @xmath146 relative to @xmath118 , provided , of course , one has also @xmath147 . on the other hand , since in these conditions @xmath148 and the @xmath149 term in eq .
( [ eq_mot_nr2a ] ) is suppressed relative to the other terms , we also see that @xmath118 obeys a schrdinger equation with binding energy equal to @xmath150 .
we shall consider now the nonrelativistic limit in the special cases @xmath151 and @xmath6 , which we will also consider later in the paper when we solve the dirac equation for harmonic oscillator potentials .
when @xmath9 , eq .
( [ eq_mot_nr2a ] ) becomes @xmath152    here we clearly see that @xmath8 plays the role of a binding potential in the nonrelativistic limit and that @xmath22 gives rise to effective binding potentials proportional to @xmath153 and @xmath154 .
this means that even a pseudoscalar potential unbounded from below could be a confining potential .
note that if @xmath155 is of the same order that @xmath156 in a @xmath157 expansion within the classically accessible region ( positive kinetic energy ) where the wave function is not negligible , and if its variation over this region is such that @xmath158 ( @xmath22 changes very little over a distance of the order of a compton wavelength ) , @xmath22 gets suppressed relative to @xmath140 .    for the case @xmath6 and if @xmath159 is very small in the classically accessible region , eq .
( [ eq_mot_nr2a ] ) becomes @xmath160 if @xmath161 or if , as explained before , the @xmath154 and @xmath162 terms are higher - order terms in a @xmath157 expansion within the classically accessible region , eq .
( [ eq_mot_nr2a_sigma=0 ] ) becomes a free particle schrdinger equation .
this is in agreement with what was found in ref .
@xcite , namely that when @xmath163 and @xmath14 is a three - dimensional harmonic oscillator potential , there is no nonrelativistic limit . in this case
we were able to show that this is true in a more general framework . a model for a relativistic particle in which the vector and scalar potentials are such that @xmath6 , and where the pseudoscalar potential is suppressed , is intrinsically relativistic , as far as bound states are concerned .
this last result deserves some more comments .
the fact that in a dirac equation one can have lorentz scalar and pseudoscalar potentials leads to remarkable results that are at odds with what is known from nonrelativistic quantum mechanics , where the lorentz structure plays no role , since one has only potentials which couple to the energy , _
i.e. _ , behave as time components of a relativistic vector . among these results
is the fact that the dirac equation is not invariant under a simultaneous shift of the energy , the scalar potential , and the pseudoscalar potential .
it has already been verified that a constant added to the pseudoscalar screened coulomb potential @xcite or to the pseudoscalar inversely linear potential asc100 is physically relevant and plays a crucial role in ensuring the existence of bound states . as an example of another notable result
, it is well known that a confining potential for a nonrelativistic fermion can not confine it in the relativistic regime when the potential is considered as a lorentz vector ( see , _ e.g. _ , @xcite ) , simply because there is pair creation and the single - particle picture does not hold anymore . it is surprising that the converse is also true , _ i.e. _ , that relativistic binding potentials may not bind in the nonrelativistic limit .
this is the case of @xmath164 ( @xmath11 ) referred before , and is basically due to the fact that vector potentials couple in a different way than scalar potentials in the dirac equation , whereas there is no such distinction in the schrdinger equation . actually , as we have seen before , @xmath8 , which is the sum of a vector and a scalar potential , plays the role of a binding potential in the nonrelativistic limit , whereas a weak enough @xmath14 can not bind a fermion with nonrelativistic energy if @xmath6 .
note that for potentials going to zero at infinity , like mean - field nuclear potentials , there are no bound states when @xmath6 , even if @xmath14 is not small compared with fermion rest energy @xcite-@xcite .    regarding antiparticles , an attractive vector potential for a particle is , of course , repulsive for its corresponding antiparticle , and vice versa .
however , the attractive ( repulsive ) scalar potential for particles is also attractive ( repulsive ) for antiparticles .
this is expressed in the change @xmath165 and @xmath166 under charge conjugation , as referred in sec .
[ subsec : charge_conj ] .
this means that an antiparticle with nonrelativistic energies ( @xmath167 ) will not be bound when @xmath9 , provided @xmath22 is small enough .
let us consider    @xmath168    as we have seen in the last section , the chiral transformation performs the changes @xmath169 , @xmath170 , @xmath171 and @xmath172 .
morever , since @xmath0 interchanges the upper and lower components ( see ( [ psi_comp_chiral ] ) ) , the resulting pair of transformed equations of motion are formally the same , so that their solutions have the same energy eigenvalues .
this symmetry can be clearly seen from the two equation pairs ( [ psi-_gen])-([2_order_gen_psi+ ] ) and ( [ psi+_gen])-([2_order_gen_psi- ] ) as well as from the isolated solutions ( [ isol_sol_delta=0 ] ) and ( [ isol_sol_sigma=0 ] ) which are converted into each other by this kind of transformation .
this means that one can restrict the discussion to the @xmath8 case ( @xmath5 ) .
the results for the case when @xmath173 can be obtained immediately by just changing the sign of @xmath30 and of @xmath26 in the relevant expressions .
the dirac spinor corresponding to the isolated solution with @xmath133 is obtained from ( [ isol_sol_delta=0 ] ) . only for @xmath174
there is a normalizable dirac spinor and , as commented before , the upper component vanishes whereas the lower component is given by @xmath175 , regardless the value of @xmath26 . for @xmath125 , eq .
( [ 2_order_delta=0_psi+ ] ) takes the form @xmath176 where @xmath177    the well - behaved solution for ( [ 2_order_delta=0_oh_psi+ ] ) , with @xmath178 necessarily real and positive , is the well - known solution of the schrdinger equation for the nonrelativistic harmonic oscillator ( see , _
e.g. _ , @xcite ) : @xmath179 where @xmath180 , @xmath181 is a normalization constant , @xmath182 is a @xmath183-th degree hermite polynomial and @xmath184    the lower component of the dirac spinor is obtained from ( psi-_delta=0 ) .
when one uses the recursion relations for the hermite polynomials ( see , e.g. , @xcite ) one gets @xmath185
e^{-\lambda ^{2}x^{2}/2 } \label{psi-_delta=0_oh}\ ] ] since the hermite polynomial @xmath186 has @xmath187 distinct zeros one may conclude that @xmath120 has @xmath183 nodes , and the expression for @xmath123 suggests that it has @xmath188 nodes , depending on the sign @xmath26 . in fact ,
if @xmath189 , this last expression reduces to @xmath190 e^{-\lambda ^{2}x^{2}/2 } \ , \label{psi-_delta=0_k1=0}\ ] ] showing that @xmath115 is proportional to @xmath191 if @xmath26 is positive . if @xmath192 , the recursion relations of hermite polynomials imply that @xmath115 is proportional to @xmath193 ( see @xcite ) .    from eqs .
( [ def_k ] ) , ( [ def_ee ] ) and ( [ delta=0_oh_ee ] ) we can get the quantization condition for the dirac eigenenergies : @xmath194 the nonrelativistic limit is reached when @xmath195 and @xmath196 . in this case
@xmath197 for small quantum numbers , and eq .
( [ eigen_eq_delta=0 ] ) becomes @xmath198 so that @xmath199 is restricted to @xmath200 .
in general there is no requirements on the signs of @xmath199 and @xmath201 , except that , as stated before , @xmath202 , and therefore , from ( [ def_k ] ) , the dirac eigenenergies corresponding to the bound - state solutions must be within the limits @xmath203 for @xmath204 , when there is a pure pseudoscalar potential linear in @xmath27 ( the two - dimensional dirac oscillator ) , eq .
( [ eigen_eq_delta=0 ] ) reduces to @xmath205\ . \label{energy_eigen_k1=0}\ ] ] note that @xmath206 for @xmath207 and @xmath208 for @xmath192 , because for the latter case the lower component is proportional to a hermite polynomial of degree @xmath209 .
this result , together with eq .
( [ energy_eigen_k1=0 ] ) , allows us to conclude that the dirac eigenvalues are given by @xmath210 , where @xmath211 , @xmath212 , irrespective of the sign of @xmath201 .
this means that _ the spectrum is independent of the sign of @xmath26_. the solutions @xmath141 and @xmath133 correspond to isolated solutions for @xmath192 and @xmath207 , respectively ( see eqs .
( isol_sol_delta=0 ) and ( [ isol_sol_sigma=0 ] ) ) . however , the spinors do depend on the sign of @xmath26 , because different signs induce different node structures , as can be seen from ( [ psi-_delta=0_k1=0 ] ) .
since the charge - conjugation transformation changes the sign of @xmath26 , it connects both kinds of solutions , but each can have both positive and negative energies .
thus a complete set of states for a certain @xmath26 includes both particle ( positive energy ) and antiparticle ( negative energy ) states but they are not related to each other by a charge - conjugation transformation .
the spectral gap between these solutions is @xmath213 and does not disappear when @xmath214 .
therefore , there is no room for transitions from positive- to negative - energy solutions and klein s paradox never comes into play .
there are no scattering states when @xmath204 , but for @xmath215 the continuum states are those ones which have energies beyond the limits expressed by ( limit_values_k1_p)-([limit_values_k1_n ] ) , for which @xmath216 .
as far as the chiral transformation is concerned , one can note that @xmath0 acts much in the same way as the charge - conjugation transformation , except that it does not change the sign of the energy and changes the mass sign .
therefore , the chiral transformation connects solutions with the same energy sign but that have opposite @xmath26 signs .
when @xmath217 , _
i.e. _ , no pseudoscalar potential , there is no isolated bounded solutions and the quantization condition expressed by eq .
( [ eigen_eq_delta=0 ] ) takes the form @xmath218 where @xmath219 stands for the sign function .
this means that when @xmath220 there are only bound states for fermions with @xmath221 and those bounded solutions are kept apart from the negative - energy continuum ( @xmath222 ) by an energy gap greater than @xmath223 . on the other hand , for @xmath224
there are only bound states for antifermions with @xmath222 and those bounded solutions also do not mix with continuum states with @xmath225 .
therefore , the positive - energy bound states of fermions never sink into the dirac sea of negative energies , and the negative - energy bound states of antifermions never reach the upper continuum , so again there is no danger of reaching the conditions for klein s paradox .
it is important to note that the solutions with @xmath226 _ are not _ the antiparticles of the solutions with @xmath227 because charge conjugation changes @xmath8 into @xmath228 and vice - versa .
this means that the antiparticles of the states with @xmath229 and @xmath227 are the negative - energy solutions of the relativistic system with @xmath230 , with @xmath25 positive .    to conclude ,
let us present , as an example of the rule stated at the beginning of the section , the eigenvalue equation for the relativistic harmonic oscillator when @xmath231 : @xmath232here @xmath233 is the number of nodes ( the degree of the hermite polynomial ) for the lower component , since the chiral transformation interchanges the upper and lower components , so that the second - order equation to solve in this case is the one for the lower component , eq .
( [ 2_order_sigma=0_psi- ] ) . the upper component would be obtained from the corresponding first - order equation eq .
( [ psi+_sigma=0 ] ) . as a final remark
, one may note that for massless particles , the @xmath5 and @xmath234 oscillators have the same spectrum with the sign of @xmath201 reversed .    in 3 + 1 dimensions the corresponding eigenvalue equations to eqs .
( eigen_eq_delta=0 ) and ( [ eigen_eq_sigma=0 ] ) have a very similar structure .
from @xcite we see that , for @xmath9 , one just adds the dirac equation quantum number @xmath19 to the coefficient of @xmath26 and modifies the coefficient of the square root to include the orbital angular momentum @xmath235 of the upper component @xmath236 while for @xmath6 one uses the orbital angular momentum @xmath237 of the lower component @xmath238 where @xmath239 ( see @xcite ) .
if we compare these equations with ( [ eigen_eq_delta=0 ] ) and ( [ eigen_eq_sigma=0 ] ) , we see that the main features are identical , and thus all the main results obtained in the next section remain valid for the relativistic harmonic oscillator in three spatial dimensions .
we turn now to the solutions of eq .
( [ eigen_eq_delta=0 ] ) . in order to obtain a deeper insight on those solutions for arbitrary
@xmath199 and @xmath201 one has to seek a convenient method for solving this transcendental equation . if we square eq .
( [ eigen_eq_delta=0 ] ) the resulting equation is in general a quartic algebraic equation in @xmath57 , which in principle can be solved analytically .
the price paid is that some spurious roots can appear in this process , although , of course , these can be eliminated by checking whether they satisfy the original equation .
a more instructive procedure is follow a graphical method , by which one seeks the intersection points of the functions of energy in eq .
( [ eigen_eq_delta=0 ] ) : a parabola on the left - hand side @xmath240 and a square - root function on the right - hand side @xmath241 in this way , it is easy to see that there are at most two acceptable solutions of ( [ eigen_eq_delta=0 ] ) , since the parabola and the square - root function can have at most two intersection points .
the construction of the graph of @xmath242 is split into three distinct classes depending on the @xmath199 values .
below we discuss these three classes of problems .
it should not be forgotten , though , that the solutions with @xmath133 are not included in this discussion , since they correspond to isolated solutions for @xmath174 .
this kind of analysis and the conclusions that will be drawn from it in the following subsections can be applied with appropriate modifications to the 3 + 1 case , and therefore to nucleon and antinucleon single - particle spectra .
this is especially true as far the relations between the spectra of systems with @xmath9 and @xmath6 are concerned , since these are obtained from the general properties of charge - conjugation and @xmath0 transformations which hold in 3 + 1 dimensions .
although we have already solved this case , it is presented here for the sake of completeness . for @xmath204
, @xmath242 is just a nonnegative constant function as can be seen from eq .
( [ energy_eigen_k1=0 ] ) , having two symmetric intersections with the parabola .
[ fig1 ] and fig .
[ fig2 ] present the results for the first four quantum numbers . in these and the following plots we use the scaled quantities @xmath243 , @xmath244 and @xmath245 .
the first shows how the values of the energy are obtained from the intersection of the two curves and the second how the eigenvalues vary with @xmath26 .
the energies for the zero mass case can also be computed from the intersections with the dashed parabola in fig .
[ fig1 ] . from these figures
it is clear the symmetry of the energy levels .
note again that the isolated solution is not plotted .
for the first four levels when @xmath246 and @xmath247 .
the black dots represent the values of energies coming from the intersection of the parabola and the horizontal lines correspond to the scaled function @xmath248 for @xmath249 .
the dashed parabola is the function @xmath250 for @xmath251 . in this case , the horizontal lines are just @xmath242 for @xmath249 , @xmath246 in appropriate units and the horizontal scale is in the units of the energy @xmath57 .
, width=377 ]     for the first four levels when @xmath252 as a function of @xmath253 .
for @xmath251 the plots are qualitatively the same , except that the lines for the different levels begin at the origin and the vertical scale would be in units of the energy @xmath57.,width=377 ]      when @xmath220 , @xmath242 is a monotonically increasing function of @xmath57 . for each @xmath183
, it can intersect the parabola at one or two points , depending on the particular values of @xmath25 and @xmath26 . to determine the conditions for this to occur
, we have first to determine the minimum of @xmath242 , which occurs when the expression under the square root in ( [ def_f_s ] ) is zero : @xmath254 this value sets a lower bound for the energies of all the solutions of eq .
( [ eigen_eq_delta=0 ] ) . for @xmath255 ,
the function @xmath242 takes the value @xmath256 and @xmath257 takes the value @xmath258 when @xmath227 , the minimum of @xmath259 occurs for negative values of the energy , that is , on the left ( and decreasing ) side of the parabola .
if @xmath260 there is only one solution of eq .
( [ eigen_eq_delta=0 ] ) and it has positive energy . from ( [ f_s_min ] ) , this can happen only for @xmath207 , since @xmath261 .
if @xmath262 there are _ always _ two solutions , one negative and one positive , because @xmath263 is always positive , no matter how large is @xmath264 , for a positive @xmath25 .
furthermore , the positive and negative solutions are such that @xmath265 , since , when @xmath266 and @xmath26 negative there is a solution with @xmath267 , but this is an isolated solution out of the sturm - liouville problem we are considering .    the considerations of the previous paragraph can be further elaborated by analyzing the behavior of the function @xmath268 @xmath269 according to what was stated in the previous paragraph , if @xmath270 one has two solutions for each @xmath183 , and if @xmath271 there is just one solution for each @xmath183 . to determine the signs of @xmath268 for each value of @xmath26 , we will need to know the values of @xmath26 which are roots of @xmath272 , or of the equation @xmath273 excluding the solution @xmath274 . applying descartes rule of signs to this last equation , we see that it can only have a positive real solution .
the rule states that the number of positive real roots @xmath275 of an algebraic equation with real coefficients @xmath276 is never greater than the number of changes of signs @xmath277 in the sequence @xmath278 ( not counting the null coefficients ) and , if less , then differs from @xmath277 by an even number .
this means that @xmath275 can take the values @xmath279 , @xmath212 . in this case , for positive @xmath25 , there is only one such sign change and therefore there is only a positive and unique solution , which we shall call @xmath280 .
the function @xmath268 changes sign at @xmath281 and @xmath282 , being negative in the interval @xmath283 and positive outside this interval .
therefore , according to what we have said before , there is only one positive - energy solution when @xmath284 and two solutions ( one positive and one negative ) when @xmath192 or @xmath285 . if @xmath274 , the two functions intersect at @xmath286 , but again , since this solution was excluded from the sturm - liouville problem leading to the eigenvalue equation ( [ eigen_eq_delta=0 ] ) , we consider only the harmonic oscillator solutions with energy greater than @xmath287 .    ) for @xmath288 and three values of @xmath253 .
the circle , square and diamond symbols mark the energy solutions for @xmath289 and @xmath290 , respectively . in the case of @xmath291 ,
the square root function for @xmath292 is also represented and the respective energy solutions are marked by the triangle symbols.,width=377 ]    the graphical representations of @xmath257 and @xmath293 are shown in fig .
fig3 for @xmath294 .
we plot the square - root function for three values of @xmath295 ( the scaled @xmath26 ) , for a positive @xmath25 .
we plot @xmath259 with @xmath296 for @xmath297 and with @xmath289 and @xmath292 for @xmath298 .
the value of @xmath299 was chosen such that @xmath268 has a positive value for @xmath300 .
there is always an infinite and unbounded set of positive - energy states with @xmath301 for any value of @xmath26 , and there are always negative states when @xmath26 is negative .
note that increasing the value of @xmath183 does not change qualitatively these features , as can be seen from the plots of @xmath259 when @xmath26 is negative .    in fig .
[ fig4 ] it is clearly seen that for values of @xmath26 between zero and @xmath302 there are only positive - energy solutions .
moreover , one can see from eq .
( [ fs_min = fp_min ] ) that the value of @xmath302 increases with @xmath25 , meaning that when @xmath303 one reaches the limit of the harmonic oscillator ( @xmath304 ) where there are only positive bound - state solutions with @xmath305 .
thus , the existence of the pseudoscalar potential gives rise to negative - energy solutions coming from the dirac sea .
for @xmath306 these kind of solutions always exist .
we excluded the isolated level with @xmath289 for negative @xmath26 .
note also that for decreasing negative @xmath307 the positive - energy solution for @xmath289 goes to @xmath308 , which was the case we discarded in the discussion of the dirac oscillator in the previous section .
when @xmath309 there is an infinite set of degenerate solutions with @xmath255 , because in this case the intersection point of @xmath259 with the negative - energy branch of @xmath250 is the initial ( minimum ) point of the square - root function , which of course does not change with @xmath183 ( see fig .
fig3 ) . however , when @xmath255 , the value of @xmath178 defined in ( [ def_k ] ) is zero , and thus , from ( [ def_lambda ] ) , @xmath310 and therefore these solutions do not correspond to bound states ( see eqs .
( [ psi+_delta=0_oh ] ) and ( [ psi-_delta=0_oh ] ) ) . for @xmath311
the degeneracy is removed and we have a very high density of very delocalized states ( since @xmath312 ) which then decreases with increasing @xmath313 . in fig .
[ fig4 ] we also plotted the parabola @xmath314 ( in scaled units @xmath315 ) which corresponds to the asymptotic values of the negative - energy states when @xmath316 and thus represents the set of lower energy bounds of the negative - energy solutions for each pair of @xmath25 and @xmath26 values .
this parabola can be also seen as delimiting the region in the energy-@xmath26 plane in which @xmath216 , _
i.e. _ , where there are no bound - state solutions .
this is related to the discussion in the previous section leading to eqs .
( [ limit_values_k1_p ] ) and ( [ limit_values_k1_n ] ) .    as referred before , for a given energy @xmath317 and positive @xmath25 there
are not bound states for values of @xmath26 in the interval @xmath318 .
in fact , from what was said above in this section , one can be more precise and state that there are no negative bound states for @xmath319 . since @xmath302 is zero
when @xmath189 and increases with increasing @xmath25 , this means that when @xmath189 ( the dirac oscillator ) there are no continuum states and when @xmath320 ( or when @xmath321 ) , all the negative - energy states are continuum states .     as a function of @xmath253 .
the @xmath289 level for negative @xmath26 and negative energy was not plotted because it corresponds to the isolated solution with @xmath322 ( @xmath267 ) .
the value of @xmath323 is marked with a cross in the energy axis .
the hollow circles mark the points ( @xmath274 , @xmath322 ) and ( @xmath282 , @xmath324 for which there are no bounded solutions . the dashed parabola represents the function @xmath325 , whose meaning is given in the text.,width=377 ]      for @xmath224 the function @xmath242 decreases monotonically with the energy , reaching again its minimum when the expression under the square root in ( [ def_f_s ] ) is zero , setting now a maximum value for the energy of the solutions of eq .
( [ eigen_eq_delta=0 ] ) .
the expressions for @xmath326 , @xmath327 and @xmath328 are the same as corresponding ones in the previous subsection , respectively eqs .
( [ e_min ] ) , ( [ f_s_min ] ) and ( [ f_p_min ] ) , except that , of course , now @xmath199 is negative .    noting that @xmath293 at @xmath133 has a nonnegative value given by @xmath329
, we see that there is always an infinite and unbounded set of negative - energy solutions ( one for each @xmath183 ) with @xmath222 . for @xmath289 and @xmath330 ,
there is a solution with @xmath331 which was already discussed in the previous section .
however , when @xmath224 , one can solve eq .
( [ isol_sol_delta=0 ] ) when @xmath133 and obtain a bound - state solution with both upper and lower components different from zero , if and only if @xmath332 .
one can show that this corresponds to a @xmath289 solution given by the functions ( psi+_delta=0_oh ) and ( [ psi-_delta=0_oh ] ) taking the limit @xmath333 .
in addition to that set of solutions , there can be another set of solutions only if the function @xmath334 ( see eq .
[ def_f(k2 ) ] ) is zero or positive , as explained before . to determine the values of @xmath201 for which this happens
, one has to compute the roots @xmath335 of @xmath336 .
when @xmath224 , descartes rule of signs applied to @xmath336 states that one can have a positive root and two or zero negative roots ( this is done by changing @xmath337 in ( [ descartes_rule ] ) ) , in addition to , of course , the root @xmath338 .
the two negative roots arise for a value of @xmath339 such that @xmath340 .
in scaled units , this value is @xmath341 , corresponding to the double negative root @xmath342 . for @xmath343
there is only a positive root , while for @xmath344 there are two negative roots and a positive root .
because @xmath334 changes sign at @xmath345 ( except when @xmath346 is a double root ) and @xmath334 is always positive when @xmath347 , this means that when @xmath343 there is only one range of values of @xmath201 where @xmath348 whereas when @xmath349 there are two ranges of @xmath201 values where this happens . within these ranges ,
only negative - energy solutions exist .
when @xmath350 , the value @xmath351 corresponds to an infinite set of degenerate negative - energy solutions with scaled energy @xmath352 , which , as mentioned before , are not bound - state solutions .
these features of the spectrum can be seen in figs .
[ fig5 ] and [ fig6 ] .     as a function of @xmath253 . as in fig .
fig4 the three non - zero values of the roots of @xmath268 , @xmath353 , are marked with a cross in the energy axis .
as in that figure , the hollow circles mark the point ( @xmath354 , @xmath322 ) and the three points where ( @xmath282 , @xmath355 ) .
the dashed parabola represents the function @xmath356 ( see text).,width=377 ]     as a function of @xmath253 .
the conventions and meaning of the symbols are the same as in figs .
[ fig4 ] and [ fig5 ] .
the dashed parabola represents the same function as in fig .
[ fig5].,width=377 ]    in fig .
[ fig5 ] the solutions with @xmath357 with @xmath358 are basically degenerate ( they follow closely the parabola @xmath359 ) and the @xmath289 solution is very close to them . since ,
as remarked before , on the parabola the value of @xmath178 is zero , these are very weakly bounded solutions .
note that the @xmath289 positive - energy solution with @xmath360 goes to @xmath287 when @xmath361 , thus becoming the isolated solution of the dirac oscillator ( @xmath189 ) .
this can also be seen in fig .
[ fig6 ] , where one may also remark that there is a value of @xmath362 around @xmath363 for which the solutions for distinct values of @xmath183 and @xmath357 are almost degenerate .
this point corresponds to the local minimum of @xmath268 , which , for @xmath364 , is very close to the double root @xmath365 .
similarly as in the case of @xmath227 , the solutions with @xmath366 have infinite degeneracy , but they do not correspond to bounded solutions . again , for @xmath26 near @xmath302 there is a high density of very delocalized states .
the parabola @xmath359 sets now an upper bound ( @xmath316 ) for the positive - energy states .    the limits @xmath367 and @xmath347 correspond , respectively , to the harmonic oscillator solutions for negative - energy ( antiparticle ) states and to the dirac oscillator .
in the first case there are no bound states for particles ( positive - energy states ) . as @xmath368 increases , positive bound states start to appear .      to complete the study of solutions of equation ( [ eigen_eq_delta=0 ] ) we present now a plot of its solutions
when @xmath26 is kept fixed and let @xmath25 vary .
note that the existence of one or two solutions is governed by the sign of the function @xmath369 identical with @xmath268 , except that now @xmath26 acts as a parameter and @xmath25 is the variable .
the critical value is now @xmath370 . in fig .
fig7 we present a plot of the solutions for @xmath207 , when there are two possible values for @xmath371 , the roots of @xmath372 .
apart the factor @xmath373 , this is just a quadratic equation in @xmath25 , so that it is easily shown that when @xmath374 it has no real solutions .
this means that , since in this case @xmath375 , there are always two sets of solutions . note that now the points in the energy-@xmath25 plane for which @xmath376 are on the hyperbole @xmath377 , such that only in the region between the hyperbole branches , where @xmath178 is positive , bound states may exist .     as a function of @xmath378 .
as before , the hollow circles mark the energy solutions for @xmath379 .
the dashed hyperbole represents the function @xmath356 ( see text).,width=377 ]    as before , we can recognize the dirac oscillator and the harmonic oscillator limits , as one sets @xmath189 and @xmath380 , respectively . in this last case , one has the harmonic oscillator solutions for particles or antiparticles depending whether @xmath25 goes to @xmath381 or @xmath382 .
one can again see that the energy levels of the dirac oscillator are symmetric about @xmath383 and that there are no scattering states . here
this follows simply from the fact that @xmath189 is an asymptote of the hyperbole branches which separate the continuum region from the bound - state region .
as referred before , the case @xmath6 can be obtained from the @xmath9 case by applying the charge - conjugation transformation .
we recall that this transformation performs the changes @xmath384 , @xmath385 and @xmath386 , or , for harmonic oscillator potentials , the changes @xmath387 and @xmath388 .
this can be seen if we solve the eigenvalue equation for @xmath6 , eq .
( [ eigen_eq_sigma=0 ] ) , for @xmath389 and plot the solutions for several values of @xmath362 .
the result is shown in fig .  [ fig8 ] .     as a function of @xmath253 and @xmath6 .
the conventions are the same as in fig .
the dashed parabola represents now the function @xmath390 .
@xmath391 stands for the quantum number of the solutions of eq .
( [ eigen_eq_sigma=0]).,width=377 ]    if we compare this figure with fig .
[ fig5 ] ( @xmath392 ) we see that the plots are identical if we reverse both the vertical and horizontal axes , _ i.e _ , if we reverse the sign of the energy and of @xmath26 , respectively
. of course , the identification of the particle and antiparticle levels is also reversed , as it should be , since we are applying the charge - conjugation transformation .    the existence of both positive- and negative - energy solutions in a system with @xmath6 can be relevant for nuclei .
in the 3 + 1 dirac equation of relativistic nuclear mean - field theories there is a connection between the ( isoscalar ) vector and tensor potential .
the latter is proportional to the derivative of the former @xcite .
for the harmonic oscillator potentials used in this paper , this would give , using the notation of ref .
@xcite , the relation @xmath393the constant @xmath394 can take values up to around 1.3 in nuclei using relativistic mean - field potentials of woods - saxon type @xcite . for harmonic oscillator mean - field potentials , there should be also an upper value for @xmath395 , and therefore the relation ( [ k_2_k_1_f_v ] ) sets a maximum for @xmath396 in nuclei .
as was seen above , this is relevant to know whether there can be simultaneously nucleon and antinucleon bound states in nuclei .
this would depend very much on the strength of the @xmath397 potential , which is given by @xmath25 . as is apparent from fig .
[ fig8 ] , for small values of the scaled @xmath25 , @xmath398 , the maximum value of @xmath399 allowed by eq .
( [ k_2_k_1_f_v ] ) could be such that there would exist only positive - energy states .
on the other hand , for stronger @xmath14 potentials , the relation ( [ k_2_k_1_f_v ] ) could allow both positive- and negative - energy bound states , as can be seen by applying the axis reversal described above to the @xmath9 spectrum shown in fig .
[ fig6 ] .    in sec .
[ sec : rel_h_oscillator ] we have seen that the @xmath0 chiral transformation can also be used to get formally @xmath6 from @xmath151 solutions , involving , however , the unphysical transformation @xmath400 .
this objection is overcome when we deal with ultrarelativistic particles , for which we may consider @xmath401 . in this case
, the chiral transformation relates massless physical states with @xmath402 and a pseudoscalar potential @xmath22 to massless states with @xmath403 and @xmath404 , where @xmath405 can be an arbitrary potential . for harmonic oscillator potentials , one has @xmath406 and @xmath388 .
since , as explained in sec .
[ sec : rel_h_oscillator ] , the chiral transformation does not change the equations of motion , the energy of their solutions is not changed .
this can be seen in figs .
[ fig9 ] and [ fig10 ] , where are plotted , respectively , the solutions of the @xmath9 and @xmath407 eigenvalue equations , eqs .
( [ eigen_eq_delta=0 ] ) and ( eigen_eq_sigma=0 ) , when @xmath251 , for @xmath408 and several @xmath26 values .
the basic features of these plots were already discussed in the previous section .
as before , the @xmath383 isolated solution was excluded .     as a function of @xmath253 for @xmath9 and @xmath251 .
the conventions are the same as in fig .
the dashed parabola represents the function @xmath409.,width=377 ]     as a function of @xmath253 for @xmath6 and @xmath251 .
the conventions are the same as in fig .
the dashed parabola represents the function @xmath409.,width=377 ]    we see that the spectra are the same apart from the @xmath388 change , which shows up as a reversal of the horizontal axis . for @xmath274 , the pure harmonic oscillator solution , both spectra are _ identical
_ , with no negative - energy solutions and , of course , no nonrelativistic limit .
note , however , that the @xmath9 and @xmath6 energy eigenstates _ are not the same _ , since @xmath0 interchanges the upper and lower components .
this amounts to say that these states are not eigenstates of @xmath0 , which of course they could not be , since , because of @xmath1 , it does not commute with the hamiltonian ( [ hamiltonian_d_sigma_delta ] ) even with @xmath410 and @xmath161 . if @xmath103 would be zero , then , because one has either @xmath411 or @xmath6 , @xmath412 would have also to be zero , and we would have a free ultrarelativistic particle , which has a good chirality , _
i.e. _ , is an eigenstate of @xmath0 .
in this work we studied the solutions of 1 + 1 dirac equation with a potential with the most general lorentz structure which leads to harmonic oscillator potentials in the second - order differential equation for the spinor components .
this is achieved when the vector @xmath412 and scalar @xmath103 potentials are quadratic and satisfy either @xmath3 or @xmath4 , and the pseudoscalar potential is linear .
we showed that there are no bound - state solutions in the nonrelativistic limit for @xmath11 and @xmath22 very small , meaning in this case we have an intrinsic relativistic problem .
we analyzed in great detail , by use of a graphical procedure , the solutions of the transcendental eigenvalue equation in the case of @xmath9 , for all signs of the potentials and for both particle and antiparticle states .
as referred in the introduction , there is an increasing interest in studying the behavior of antinucleon in nuclei , and in particular how well pseudospin and spin symmetries hold for nucleons and antinucleons respectively .
we discussed in detail the isolated solutions , _
i.e. _ , solutions out of the sturm - liouville problem , as well as the critical points in the parameter space of @xmath25 and @xmath26 , the coefficients of the quadratic and linear potentials . near these points
there is a high density of bounded states which are very extended in space .
these critical points are also relevant to determine which values of @xmath25 and @xmath26 allow to have a spectrum of both fermion and antifermion bounded solutions simultaneously or just one of these kind of solutions .
we established a connection to 3 + 1 dimensional relativistic oscillator harmonic potentials applied to nuclei , where there is a relation between those coefficients and concluded that , for a sufficiently strong @xmath14 potential , both kinds of solutions ( _ i.e _ , for both nucleons and antinucleons ) can coexist .
we showed how to obtain the solutions for @xmath6 from the @xmath9 case , using the charge - conjugation and chiral transformations .
in particular , we used the chiral transformation to show that , for massless fermions and zero pseudoscalar potential , the @xmath5 and @xmath11 solutions have the same spectrum .
these conclusions are quite general and do not depend on the number of dimensions so that they remain true in 3 + 1 dimensions .
also , most of the features of the 1 + 1 spectra presented in this work are also present in the 3 + 1 relativistic harmonic oscillator . indeed , since the 1 + 1 eigenvalue equations for @xmath9 and @xmath6 are very similar to the corresponding 3 + 1 equations , we can apply the same graphical methods to the 3 + 1 case , and therefore most of the analysis of the 1 + 1 case remains true with appropriate changes .
for instance , the shape of the parabola in the energy-@xmath26 plane related to @xmath413 would depend now on the value of the spin - orbit quantum number @xmath19 .
the function @xmath268 would also depend on @xmath19 , and therefore the roots @xmath302 and critical value @xmath414 would also depend on @xmath19 .
for @xmath274 there would be the known harmonic oscillator @xmath415 degeneracy , produced by the u(3 ) symmetry studied very recently by ginocchio @xcite in connection to spin and isospin symmetry .
finally , the chiral transformation enables us to switch between the so - called pseudospin ( @xmath6 and @xmath161 ) and spin symmetries ( @xmath5 and @xmath161 ) @xcite _ for any potential_. in the context of a study of the symmetries of a 3-dimensional relativistic harmonic oscillator , ginocchio already showed that one can transform the generators of the su(2 ) spin symmetry into the pseudospin generators by a @xmath0 transformation @xcite .
we showed that this transformation allows us to relate the @xmath402 and @xmath403 problems for any potential @xmath405 , and , for massless fermions and zero pseudoscalar potential , also to conclude that the two energy spectra are the same .