quantum correlations have a central role in quantum information processing , in quantum foundations , as well as in the physics of strongly correlated systems  @xcite . on one hand , quantum correlations , and in particular entanglement , are a resource that allows one to go beyond what classically possible in many scenarios , from communication tasks , to ( measurement - based ) quantum computing , to quantum cryptography . on the other hand ,
the non - classicality of quantum correlations  be it in the form of non - locality , steering , entanglement , or discord  is one of the most distinctive traits of quantum mechanics , and challenges our understanding of quantum mechanics itself .
the interplay between local and global properties of quantum states is a key aspect in the study of quantum correlations , both from a fundamental perspective and an applicative one .
for example , we may want to certify the presence of multipartite entanglement in large systems without the  often inaccessible  knowledge of the global state , using instead only the information that comes from reduced states . on the other hand , in condensed - matter physics , because of the typically local
, two - body  interactions , relevant properties are dictated by the interplay between the allowed reduced states and global correlations , giving raise to phenomena like frustration  @xcite .
the general study of the relations between the properties of the reduced states and the properties of the global state is known as the quantum marginal problem , which has seen a growing interest in the past years also for the reasons above  @xcite .    in this work
we study what can be inferred about the quality of the correlations of the global state given information about the two - body reduced states , aiming at answering the question : _ what correlations need to be present globally to explain what we see locally ? _ in  @xcite a characterization of multipartite entanglement in terms of
even just single - party reduced states ( actually , single - party spectra ) was given , but under the assumption of dealing with a pure or quasi - pure global state .
in  @xcite the possibility of dealing with global mixed states is taken into account , and examples are given where two - qubit separable states are only compatible with global entanglement , intended in the sense of lack of total separability ( see section  [ sec : definitions ] for definitions ) . in  @xcite examples
are also given where genuine multipartite entanglement  a much stronger notion of global entanglement  can still be deduced from the properties of the two - body reduced states , but only when these reduced states exhibit bipartite entanglement themselves .
in this work we present several results that complement and generalize those of  @xcite . we offer a brief summary of our findings in table  [ table ] .
firstly , we provide examples of triples of bipartite reduced states  in the simplest case , two - qubit states  that , albeit separable , are only compatible with genuine tripartite entanglement ( lower - right corner of table  [ table ] ) . as far as we know , this is the strongest `` gap '' known between the entanglement properties of the marginals and of the global state .
secondly , we address the issue of relating the general quantumness of correlations  @xcite of reduced states to the quantum correlations of the global state .
we find that strictly classical reduced states may still be compatible only with global entanglement ( upper - left corner of table  [ table ] ) .
thirdly , we find that , at least for qubits , the strict classicality of the two - body correlations makes it impossible to certify genuine tripartite entanglement based on the knowledge of the reductions : strictly classical two - qubit reduced states are always compatible with a global state that is not genuine tripartite entangled ( upper - right corner of table  [ table ] ) .
the rest of the paper is organized as follows . in section  [ sec : definitions ] we define the relevant notions of correlations and classicality , and of compatibility of two - body reduced states in tripartite systems . in section  [ sec : classicalcomp ] we study the relation between the classicality of reductions and their compatibility . in section  [ sec : genuine ] we prove that unentangled reduced states may only be compatible with genuine multipartite entanglement at the level of the global state .
finally , we conclude in section  [ sec : conclusions ] .
l l|>m1.5cm|>m3.5cm| & & + & & entangled &     genuinely multipartite + entangled    &    fully + classical    & & ( qubits ) &    separable    &  @xcite &
we begin by formally defining qualitatively different types of correlations . [ def : separable ] any tripartite mixed state can be written as a mixture of an ensemble of pure states as @xmath0 .
we say that @xmath1 is :    * _ fully separable _ , if we can take each @xmath2 to be fully factorized , e.g. @xmath3 ; * _ biseparable _ , if we can take each @xmath2 to be unentangled in at least one partition , e.g. , @xmath4 , @xmath5 , or @xmath6 ; * _ genuinely multipartite entangled _ , if for any ensemble there is at least one @xmath7 with @xmath8 that is not factorized with respect to any bipartition , i.e. , if @xmath1 is not biseparable ; * _ fully classical _ , if we can take each @xmath2 to be of the form @xmath9 , with @xmath10 orthonormal basis on @xmath11 , @xmath12 , and @xmath13 , respectively , so that , overall , @xmath14 .
bipartite full classicality and separability are defined similarly : @xmath15 is _ fully classical _ if @xmath16 , for @xmath17 orthonormal bases , and _ separable _ if @xmath18 .
a bipartite state is _ entangled _ if it is not separable .
the notions of full separability and biseparability are redundant for bipartite states .
the set of fully classical states , @xmath19 , is a subset of the set of fully separable states , @xmath20 , which in turn is a subset of the set of biseparable states , @xmath21 ; the set of genuinely multipartite entangled states , @xmath22 is the is the complement of @xmath21 in the space of all states , @xmath23 ( see figure  [ fig : sets ] ) .
all the mentioned sets apart from @xmath19 and @xmath22 are convex
. a biseparable state may be either separable or entangled in any bipartition , but it is by definition the convex sum of three states that are each separable in one of the bipartitions @xmath24 , @xmath25 , and @xmath26 .
we now move to define formally the notion of compatibility for reduced states .
[ def : compatible ] given a triple of three two - party states @xmath27 , its _ compatibility set _ is defined as @xmath28
. any compatibility set is a convex set  @xcite , and the property of being part of a given compatibility set defines an equivalence relation .
we find it useful to denote by @xmath29 the compatibility set associated with the reduced states of @xmath30 , i.e. , the set of all states that have the same reductions as @xmath30 .
a triple of two - party states @xmath31 is said to be _ compatible _ ( so that we refer to the triple as _ triple of reductions _ ) if @xmath32 , i.e. , if there is at least one global state with those reductions .
the following definition links the compatibility of reduced states to the correlation properties of global states .
[ def : compatible2 ] we say that the reductions @xmath31 are _ incompatible with a set @xmath33 _ ( or with the defining correlation property of @xmath33 ) if @xmath34 .
we say that ( a compatible ) @xmath31 is only compatible with genuine multipartite entanglement if it is incompatible with @xmath21 .
figure  [ fig : sets ] illustrates the problem of deciding whether certain bipartite reductions necessarily require the presence of global correlations of a certain kind .    we can always expand a generic tripartite state as @xmath35 where @xmath36 is the reduced states of party @xmath37 and @xmath38 is the normalized identity operator on the hilbert space @xmath39 . the bipartite correlation matrices @xmath40
can be defined via @xmath41 , and satisfy @xmath42 = { \tr}_l [ \chi_{kl}]=0 $ ] .
it is worth noting that when a bipartite marginal state @xmath43 is fully classical , then @xmath44 = 0 $ ] , and @xmath45 is also diagonal in the same product basis as @xmath46 .
the tripartite correlation matrix @xmath47 , which for a fixed @xmath1 can be defined via , satisfies @xmath48 = 0 $ ] , for all @xmath49 . for compatible reductions @xmath50 ,
the compatibility set @xmath51 is spanned by choosing the tripartite correlation matrix @xmath47 so that the resulting operator in is a physical state , i.e. , positive semidefinite . on the other hand , to determine whether a triple of bipartite states is compatible , we first check the basic necessary condition that the single party marginals be the same , i.e. , @xmath52 = \tr_k[\rho_{ik}]$ ] for all @xmath53 .
next , we have to search for a tripartite correlation matrix , @xmath47 , such that eq .
is physical .
if no such @xmath47 exists , the given states are not compatible .
consider the marginals from the well - known greenberger - horne - zeilinger ( ghz ) state @xmath54 : @xmath55 .
these are fully classical marginals coming from a genuinely tripartite entangled state
. however , these marginals are also compatible with @xmath56 , which is fully - classical .    in this section
we will first provide an example where fully classical two - body reduced states are not compatible with a global fully classical state , and actually require the presence of entanglement . on the other hand , we will prove that , in the case of three qubits , the fully classical two - body reductions are always compatible with a global states that is not genuine multipartite entangled .      we will now derive conditions to ensure that some fully - classical marginals can not be compatible with any global fully - classical state .
we start with the following lemma .
[ le : commutatorlemma ] suppose three states @xmath57 , @xmath58 , and @xmath59 are fully classical .
consider the commutators @xmath60 = \left [ \rho_{ij } \ox \openone_k , \ ; \rho_{ik } \ox \openone_j \right]$ ] , where the second equality is due to the assumed classicality , i.e. , @xmath61 = 0 $ ] .
then :    a.   all commutators @xmath62 vanish if and only if there are orthonormal basis @xmath63 , @xmath64 , @xmath65 , such that @xmath66 , @xmath67 , and @xmath68 . b.   if some commutator @xmath62 does not vanish , then : a.   at least one @xmath69 is degenerate ( i.e. , at least two eigenvalues of some @xmath69 are the same ) ; b.   there does not exist a tripartite fully classical state that is compatible with @xmath57 , @xmath59 , and @xmath58 .
\(i ) the `` if '' part is trivial .
let us prove the `` only if '' part .
by hypothesis we may assume @xmath70 with the orthonormal basis @xmath71 on @xmath11 , @xmath72 on @xmath12 , @xmath73 on @xmath13 , and @xmath74 , @xmath75 ( similarly for @xmath76 , etc ) .
then @xmath77 implies @xmath78 \label{eq : commute } \notag\\ & \quad=\left [ ( \r_{ij}-\r_i\ox\r_j ) \ox \openone_k , \ ; ( \r_{ik}-\r_i\ox\r_k ) \ox
\openone_j \right ] \notag\\ & \quad=\left [ \r_{ij } \ox \openone_k , \ ; \r_{ik } \ox \openone_j \right]=0,\end{aligned}\ ] ] with @xmath79 . by setting @xmath80 in ,
we have @xmath81=0 $ ] , @xmath82 . thus the states @xmath83 are simultaneously diagonalizable in the orthonormal basis @xmath84 .
so we may replace the basis @xmath63 and @xmath85 in and by @xmath84 .
this replacement may result in the change of @xmath86 and @xmath87 . since there is no confusion ,
we still use them in and .
next by setting @xmath88 in , we have @xmath89=0 $ ] , @xmath82 . thus the states @xmath90 are simultaneously diagonalizable in the orthonormal basis @xmath91 .
so we may replace the basis @xmath64 and @xmath92 in and by @xmath91 .
this replacement may result in the change of @xmath93 . since there is no confusion ,
we still use them in .
third we set @xmath94 in and repeat the above argument to show that the basis @xmath65 and @xmath95 in and can be replaced by the orthonormal basis @xmath96 .
so the assertion follows .
\(ii ) suppose either condition ( a ) or ( b ) is violated .
we have that either @xmath97 , @xmath98 , and @xmath99 are all non - degenerate respectively in the orthonormal basis @xmath63 , @xmath64 , @xmath65 ( violation of ( a ) ) , or that there exists a tripartite fully classical state @xmath100 that is compatible with @xmath57 , @xmath59 , and @xmath58 ( violation of ( b ) ) .
in either case we have @xmath101 , @xmath102 and @xmath67 .
so ( i ) implies that all commutators @xmath103 vanish , and we have a contradiction .
an immediate consequence of lemma  [ le : commutatorlemma ] is the following .
[ thm : commuteimpliesclassical ] let @xmath57 , @xmath59 , and @xmath58 be three compatible bipartite fully classical states , such that    \(i ) they all commute ( all commutators @xmath103 of lemma  [ le : commutatorlemma ] vanish ) ,    or    \(ii ) all three one - body reduced states @xmath97 , @xmath104 and @xmath99 are non - degenerate .
then @xmath57 , @xmath59 , and @xmath58 are compatible with a fully classical tripartite state .
the fact that all three single - system reductions are not degenerate implies , by lemma  [ le : commutatorlemma].(iia ) , that all the commutators @xmath103 defined in lemma  [ le : commutatorlemma ] vanish . by lemma  [ le :
commutatorlemma].(i ) we have that @xmath97 , @xmath98 , and @xmath99 , are diagonal in the orthonormal bases @xmath63 , @xmath105 , and @xmath106 , respectively , in which @xmath57 , @xmath59 , and @xmath58 are explicitly classical .
most importantly , we have @xmath107 with @xmath108 and @xmath109 .
let @xmath1 be any tripartite state with which the three two - body reductions are compatible .
then also the fully classical tripartite state @xmath110 has bipartite reduced density matrices @xmath57 , @xmath59 , and @xmath58 .
given theorem [ thm : commuteimpliesclassical ] , in order to construct an example where @xmath57 , @xmath58 , and @xmath59 are all fully classical but not compatible with any fully - classical state , we have first of all to construct an example where @xmath57 , @xmath58 , and @xmath59 are classical but do not commute with each other . for this
, we will need the following lemma .
we recall that for a bipartite state @xmath111 acting on the hilbert space @xmath112 , the partial transpose computed in the standard orthonormal basis @xmath113 of system @xmath114 , is defined by @xmath115 .
one can similarly define the partial transpose @xmath116 on the system @xmath117 .
[ le : classicalmain ] consider three classical - classical two - qubit states    @xmath118    where @xmath119 , and any one of @xmath63 ,
@xmath64 , @xmath65 is a real and orthonormal basis in @xmath120 .
let @xmath121 then :    a.   if @xmath122 then @xmath1 is separable with respect to to the partition @xmath24 , @xmath25 and @xmath26 .
b.   @xmath57 , @xmath59 and @xmath58 are compatible if and only if they are compatible with the biseparable state @xmath1 in . c.   suppose @xmath57 , @xmath59 and @xmath58 are compatible
they are compatible with a fully separable state if and only if @xmath1 is fully separable .
\(i ) one may directly verify that the state is invariant under partial transposition with respect to any system , i.e. , @xmath123 for @xmath124 .
since @xmath122 , the assertion follows from  ( * ? ? ?
* theorem 2 ) .
\(ii ) the `` if '' part is trivial ; let us prove the `` only if '' part .
suppose the bipartite marginals @xmath57 , @xmath59 and @xmath58 are compatible with a tripartite state @xmath125 .
since @xmath57 , @xmath59 , and @xmath58 are real , they are also compatible with the real state @xmath126 , so we can assume that @xmath127 is real without loss of generality . by eq .
( 1 ) , there is a hermitian matrix @xmath47 such that @xmath128 since in our case both @xmath129 and @xmath130 are real , also @xmath131 is real .
it follows from eqs .
 , and the fact that @xmath63 , @xmath64 , @xmath65 are real and orthonormal bases , that @xmath57 , @xmath59 and @xmath58 are invariant under the local unitary @xmath132 .
so they are compatible with the state @xmath133 where we have used that , from  , @xmath134 we will now argue that , for a real @xmath131 , @xmath135 so that proves that , for a physical state @xmath130 , @xmath129 is also physical , as it corresponds to the convex combination of physical states .
the starting point in proving is to observe that every three - qubit correlations matrix @xmath131 is by definition the linear combination of traceless pauli matrices , i.e. @xmath136 since @xmath131 is hermitian , all coefficients @xmath137 are real .
moreover , for a real ( and hence symmetric ) @xmath131 only terms with an even number of @xmath138 are present in the expansion , because @xmath139 , while @xmath140 and @xmath141 are symmetric . on the other hand , @xmath142 for @xmath143 , while , obviously , @xmath144 .
since each non - zero term in the expansion   of @xmath131 contains an even number of @xmath145 s , it will change sign after conjugation by @xmath146 , i.e. , @xmath147 as argued , this implies @xmath122 , with biseparability following from ( i ) .
\(iii ) the `` if '' part follows from ( ii ) , let us prove the `` only if '' part .
suppose @xmath57 , @xmath59 and @xmath58 are compatible with a fully separable state @xmath125 . from , @xmath1 is the convex sum of a few fully separable states .
so the assertion follows .
this completes the proof .
we are now ready to present our example where @xmath57 , @xmath58 , and @xmath59 are all fully classical but not compatible with any fully - classical state .
[ eg : commutatorex ] consider the three - qubit state @xmath148 where @xmath149 , with @xmath150 are the pauli matrices and @xmath151 ( for @xmath152 outside of this interval the matrix @xmath153 is not positive semi - definite ) , @xmath154 .
it is not hard to see that each of the bipartite marginals states is fully classical but with respect to different bases .
a quick method to verify this assertion is using theorem 1 in  @xcite .
moreover the reductions of @xmath153 do not commute with each other .
part ( ii ) of lemma  [ le : commutatorlemma ] implies that there is no tripartite fully classical state that is compatible with these bipartite marginals : @xmath155 .
surprisingly , we can find a value of @xmath152 for which the fully classical marginals in fact require some global entanglement .
consider the state @xmath156 .
the range of @xmath157 does not contain any fully factorized pure state , hence it can not be fully separable .
now , up to a local unitary , the classical bipartite marginals @xmath158 , @xmath159 , @xmath160 , and the entangled state @xmath157 can be written as  - , respectively .
it follows from lemma [ le : classicalmain].(iii ) that @xmath158 , @xmath159 , and @xmath160 can not be compatible with any fully separable state .    in appendix
[ sec : appuni ] , which focuses on the uniqueness of global states with fixed two - body reductions , we provide also an example of classical reductions compatible with a _
unique _ global state that is not fully classical , although fully separable .
although , as we have just seen , compatible classical marginals may require global quantum correlations or even entanglement , it turns out that for the case of _ three qubits _ they will never require the global state to be genuinely multipartite entangled .    to see this , will need an additional lemma .
[ le : nongenuine ] let @xmath63 , @xmath64 be an orthonormal basis on @xmath11 and @xmath12 , respectively .
suppose the bipartite marginals @xmath57 , @xmath58 and @xmath59 are compatible .
they are compatible to a non - genuinely entangled tripartite state when one of the following conditions is satisfied : ( i ) @xmath161 and @xmath162 .
( ii ) @xmath163 .
suppose @xmath57 , @xmath58 and @xmath59 are compatible with @xmath1 . if hypothesis ( i ) is satisfied , then they are also compatible with the state @xmath164 , which is biseparable . on the other hand
if hypothesis ( ii ) is satisfied , then it follows from  @xcite that there is a quantum channel @xmath165 on @xmath13 such that @xmath166 , for @xmath167 .
so we obtain @xmath168 and @xmath169 .
then @xmath57 , @xmath58 and @xmath59 are compatible with the fully separable state @xmath170 .
this completes the proof .
now we are ready to give the proof of the following .
[ thm : classical ] any three compatible classical - classical two - qubit states @xmath57 , @xmath59 and @xmath58 are compatible with a tripartite biseparable state .
suppose some fully classical bipartite marginals @xmath57 , @xmath59 and @xmath58 are only compatible with genuinely entangled states @xmath1 . then the one - party reduced density operators @xmath171 have to be the maximally mixed states , @xmath172 .
indeed , if , without loss of generality in the argument , @xmath173 is non - degenerate , the two - party reduced states would also be compatible with a global state given by the a locally ( on @xmath114 ) dephased version of @xmath1 , which would be separable in @xmath24 , leading to a contradiction .
thus , up to local unitaries , we may assume @xmath174 where @xmath175 and @xmath176 . since @xmath177 , we have @xmath178 .
so we obtain @xmath179 $ ] . since @xmath1 is genuinely entangled , the cases @xmath180 and @xmath181 are excluded by lemma  [ le : nongenuine].(i ) and  [ le : nongenuine].(ii ) , respectively .
so we obtain @xmath57 as eq .  .
by similar arguments and performing suitable diagonal local unitary gates on systems @xmath182 , the classical - classical two - qubit states @xmath59 and @xmath58 can be simplified to the forms eqs .   and , respectively .
meantime , @xmath57 is unchanged . since @xmath57 , @xmath59 and @xmath58
are compatible , it follows from lemma  [ le : classicalmain].(ii ) that they are compatible with a biseparable state .
it gives us a contradiction .
so there are no compatible bipartite marginals that are only compatible with genuine entangled states .
this completes the proof .
we have seen that the condition of classicality of marginals is strong enough to exclude the need for global genuine multipartite entanglement .
we will now construct non - classical separable marginals that are only compatible with global genuine multipartite entanglement , but first we need to establish some more definitions and notation .
we set @xmath183 , @xmath184 and @xmath185 .
we denote @xmath186 , @xmath187 the rank and range of any square matrix @xmath188 , respectively .
a quantum state is a positive semidefinite linear operator @xmath189 with @xmath190 .
we say @xmath1 is a @xmath191 state by meaning that the reduced density operators satisfy @xmath192 , @xmath193 , @xmath194 .
the ranks of the reduced density operators of @xmath1 are invariant when we perform an invertible local operator ( ilo ) on @xmath1 .
that is , let @xmath195 such that @xmath196 .
then @xmath197 , @xmath198 and @xmath199 where @xmath200 .
we also denote @xmath201 as the vector whose components are the complex conjugate of those of @xmath202 .
so @xmath202 is real when @xmath203 .
evidently @xmath204 , where @xmath205 , we recall , denotes partial transposition .
we call the integer pair @xmath206 the _ birank _ of @xmath111 , and the two integers may be different . for such examples of two - qubit and qubit - qutrit separable states ,
we refer the readers to  ( * ? ? ?
* table i , ii ) .
furthermore , we say that @xmath111 is a ppt [ npt ] state if @xmath207 [ @xmath208 has at least one negative eigenvalue ] .
evidently , a separable state must be ppt .
the converse is true only if @xmath209  @xcite .
we say a bipartite state @xmath57 is _
@xmath114-finite _ when for any subspace @xmath210 , @xmath211 and any state @xmath212 , it holds that @xmath213 . in other words
, @xmath57 is not @xmath114-finite when @xmath214 contains a 2-dimensional subspace spanned by @xmath215 with some linearly independent states @xmath216 .
so if @xmath57 is not @xmath114-finite , there must be infinitely many product states in @xmath214 .
besides these notions and notation , we will need also the following lemma .
[ le : symmetricab ] suppose the bipartite marginals @xmath57 , @xmath59 and @xmath58 are compatible with @xmath1 , and @xmath57 is @xmath114-finite and @xmath117-finite
. then @xmath1 is either separable with respect to the partition @xmath217 , or genuinely multipartite entangled .
suppose @xmath1 is biseparable , so @xmath218 , with @xmath219 separable in the @xmath24 partition ( similarly for @xmath220 and @xmath221 .
we argue that @xmath222 is fully separable , and a similar argument will apply to @xmath223 .
let @xmath224 . since @xmath57 is @xmath117-finite and @xmath225 , any @xmath226 must be a product state .
so @xmath222 is fully separable .
similarly one can show that @xmath223 is also fully separable , so @xmath1 is separable with respect to @xmath217 .
we are now in the position to prove the following .
[ thm : main ] suppose the triple @xmath227 is compatible .
then , @xmath228 if all the following conditions are met : ( i ) for any @xmath229 , the state @xmath230 is @xmath231-finite and @xmath232-finite ; ( ii ) @xmath59 has birank @xmath233 ; ( iii ) @xmath57 has birank @xmath234 , @xmath235 .
suppose @xmath57 , @xmath59 and @xmath58 are compatible with a biseparable state @xmath1 . by hypothesis ( i ) and
lemma  [ le : symmetricab ] , we obtain that @xmath1 is separable with respect to the partition @xmath217 .
let @xmath236 where @xmath8 , the @xmath237 on @xmath13 are pairwise linearly independent .
note that the product subspace @xmath238 , @xmath239 .
this fact and ( i ) imply that @xmath240 , @xmath239 . by similar arguments we have @xmath241 .
we may assume @xmath242 .
hence @xmath243 . by ( ii ) we have @xmath244 where @xmath245 . without loss of generality
, we may assume that the states @xmath246 , @xmath247 span @xmath248 .
we choose a suitable ilo @xmath249 such that @xmath250 , @xmath251 and @xmath252 is real . by performing @xmath249 on the state @xmath1
, we have @xmath253 and @xmath254 for any @xmath231 .
since the operation @xmath249 does not change the rank of quantum states , it follows from ( ii ) that @xmath255 has birank @xmath256 .
recall that the @xmath237 are pairwise linearly independent .
since ( i ) is not changed under ilos , it follows from ( i ) that @xmath252 is not parallel to any state @xmath257otherwise @xmath258 would not be @xmath117-finite . since @xmath252 is real , the two @xmath259-dimensional subspaces @xmath260 and @xmath261 are equal and spanned by @xmath262 and @xmath263 .
so the states @xmath264 , for any @xmath265 .
then ( i ) implies that these @xmath266 are real up to an overall phase .
so @xmath267 .
it implies @xmath268 which is a contradiction with ( iii ) .
therefore @xmath57 , @xmath59 and @xmath58 are not compatible to any non - genuinely entangled state .
this completes the proof .
we will now make use of theorem  [ thm : main ] and offer an example of separable marginals only compatible with genuinely entangled tripartite states .
[ ex : main ]    consider the family of rank-@xmath259 states on @xmath269 given by @xmath270 , with @xmath271 @xmath272 , @xmath273 , and @xmath274 .
it is easy to see that the only biseparable pure state in @xmath275 is @xmath276 .
the bipartite reduced density operators of @xmath277 are @xmath278 with @xmath279 and @xmath280 .
the three two - qubit marginals @xmath281 , @xmath255 and @xmath282 are positive under partial transposition ( ppt ) , so they are separable  @xcite .
hence , @xmath57 , @xmath59 and @xmath58 are separable too ; they also evidently satisfy condition ( i ) of theorem  [ thm : main ] .
furthermore , @xmath59 has birank @xmath283 , while @xmath284 , and @xmath57 has birank @xmath285
. so also conditions ( ii ) and ( iii ) of theorem  [ thm : main ] are satisfied , and we conclude that @xmath57 , @xmath59 and @xmath58 are only compatible with genuinely tripartite entangled states .
the example shows that for any fixed local dimension @xmath286 , there exist triples of two - qudit separable states that are only compatible with genuine multipartite entanglement .
the `` core '' of our construction is the genuine multipartite entangled three - qubit state @xmath287 of eq .  .
it turns out that @xmath287 is actually the only state compatible with its reductions .
the proof of this is given in the appendix .
it is worth comparing this with the results of  @xcite .
there it was proven that for almost all pure entangled states of three qubits @xmath288 it holds @xmath289 , with the exception of states of the generalized ghz form @xmath290 ( up to local unitary transformations ) , which satisfy , e.g. , @xmath291 .
interestingly , the only three - qubit pure states that have separable reduction are of the generalized ghz form  @xcite .
this implies that any three - qubit state @xmath129 such that ( i ) its reductions are separable and ( ii ) @xmath292 , must be mixed . since the state @xmath287 has rank two , we can think of it as the simplest possible example that satisfies ( i ) and ( ii ) , with the additional property of being uniquely determined by its reductions .
we generalize example  [ ex : main ] in several ways , all presented in the appendix .
while we showed that there exist genuine multipartite states whose compatibility set contains only genuine multipartite states , it is natural to ask how common this phenomenon is , i.e. , whether such states have finite volume in the set of all states .
this is important also from the point of view of the potential realization of such states in the lab , which can never be perfect .
we answer this question in the affirmative .
we introduce a parameter of compatibility of a tripartite state @xmath30 with @xmath293 as @xmath294 , where we have used the hilbert - schmidt norm @xmath295   in its arguments and the fact that @xmath296 is positive and vanishes if an only if @xmath297 . ] .
we further define @xmath298 .
we have @xmath299 for any triple @xmath300 such that @xmath301 , even if the triple of reduced states is compatible , as in example  [ ex : main ] .
finally , given a tripartite state @xmath287 , we define @xmath302 and @xmath303 .    now , consider a genuinely entangled multipartite state @xmath304 with separable reductions such that @xmath305 , and the convex combination of @xmath304 with an arbitrary fully separable state @xmath306 : @xmath307 for @xmath308 .
since the set of biseparable states is closed , there exists @xmath309 such that @xmath310 is genuine multipartite entangled for all @xmath306 and all @xmath311 . since @xmath306 is fully separable , so are the two - party reduced states of @xmath312 . furthermore , since @xmath313 is continuous , there exists @xmath314 such that @xmath315 for all @xmath306 and all @xmath316 .    for any local finite dimensions ,
the set of fully separable states has non - zero volume among all states , because there exists a ball of fully separable states around the maximally mixed state  @xcite .
thus , the argument above proves that also the set of tripartite states whose two - party marginals are separable but only compatible with genuine multipartite entanglement has non - zero volume .
we analyzed the relation between the character of correlations of tripartite states and the ones exhibited by their bipartite reductions , i.e. , a version of the quantum marginal problem that focuses on the compatibility of bipartite reductions with certain global properties .
we constructed examples where separable reductions are only compatible with genuine multipartite entanglement . up to our knowledge , this is the `` largest '' known separation between the character of correlations of bipartite reductions and what can be inferred about the quality of correlations of the global state , based only on the knowledge of the reductions . on the other hand , at least for qubits we were able to prove that compatible reductions that are fully classical can always originate from a biseparable global state .
nonetheless , bipartite reductions that are fully classical may still require the presence of some entanglement in the global state .
our results show that the relation between global and `` local '' correlations is far from trivial .
notably , the notion of fully classical correlations is strong enough to `` break '' the need for genuine multipartite entanglement , but not the potential need for global entanglement altogether .
an interesting open question is whether compatible classical - classical marginals in high dimension are always compatible with a biseparable tripartite state .
another question is how to quantitatively bound the certifiable genuine multipartite entanglement in terms of the non - classicality of the two - body reductions , at least in the three - qubit case .
the latter problem is reminiscent of the case of entanglement distribution , where the non - classicality of correlations  rather than the entanglement  present between a quantum carrier and distant labs constitutes a bound on the entanglement that can be generated between the labs by exchanging the carrier  @xcite .
we are grateful to a. acn and o. ghne for useful discussions .
lc started this work when at the institute for quantum computing & department of pure mathematics , university of waterloo , and was mainly supported by mitacs and nserc ; lc s research was funded in part by the singapore national research foundation under nrf award nrf - nrff2013 - 01 .
og is grateful for the support of the austrian science fund ( fwf ) and marie curie actions ( erwin schrdinger stipendium j3312-n27 ) .
km was supported by the john templeton foundation , national research foundation and the ministry of education ( singapore ) during the completion of this work .
mp acknowledges support from nserc , cifar , and ontario centres of excellence .
we first prove that @xmath287 in eq .   in the paper is the only state compatible with its reductions , a fact of interest in its own . [ le : entunique ] for @xmath287 in eq .
it holds @xmath317 .
suppose @xmath318 has the same reductions as @xmath287 , i.e. , @xmath319 .
we can always write its spectral decomposition as @xmath320 , where @xmath321 , and @xmath322 , @xmath239 .
we have @xmath323 .
it follows from eq .
( 5 ) in the paper that @xmath324 . so
any four states of @xmath325 are linearly dependent . using the freedom in the choice of the pure - state ensemble representation of a mixed states  @xcite
, we can choose a suitable linear combination of @xmath226 , @xmath326 , such that it is equal to @xmath327 .
so the state can be written as @xmath328 where @xmath329 . by applying this procedure to another four states
@xmath330 with @xmath331 respectively , we can realize @xmath332 .    by relabeling the states ,
we can write @xmath333 with @xmath334 on @xmath335 . we have @xmath336 and @xmath337 by eq .  .
since @xmath281 is x - finite for @xmath338 , we have @xmath339 . by the similar argument
we can show @xmath340 .
so we have @xmath341 .
let @xmath342 where @xmath343 and @xmath344 . by eq .
we have @xmath345 , @xmath346 . by eq .  , we have @xmath347 and @xmath348 .
these equations imply for @xmath343 , we have    @xmath349    the coefficients of @xmath350 in both @xmath281 and @xmath282 are @xmath351 , so @xmath352 . by replacing @xmath353 by a suitable linear combination of them , we may assume @xmath354 .
so the tripartite state can be rewritten as @xmath355 , where @xmath356    since @xmath324 , we have @xmath357 . by eq .
, we have @xmath358 .
so @xmath359 . by eq .
again , we have @xmath360 , @xmath361 , and @xmath362 .
now we see @xmath363 in eq .  .
this completes the proof .
we further derive ( and later use in appendix [ sec : appgen ] ) the following lemma .    [
le : unique ] suppose @xmath57 , @xmath59 , @xmath58 are only compatible with a tripartite state @xmath1 , and @xmath281 , @xmath255 , @xmath282 are compatible with another tripartite state @xmath277 .
if @xmath364 , then @xmath277 is the only state with which @xmath281 , @xmath255 , @xmath282 are compatible .
suppose @xmath365 are compatible with another state @xmath366 .
since @xmath367 , we may find a small enough @xmath368 and a tripartite state @xmath369 , such that @xmath370 .
so the bipartite reductions @xmath371 are compatible with the state @xmath372 , which is different from @xmath1 .
it gives us a contradiction .
we conclude this section presenting separable marginals that are only compatible with a unique quantum correlated ( unentangled ) state .
[ pp : sepunique ] the separable states @xmath373 , where @xmath374 , are only compatible with the separable state @xmath375 .
we will use the following observation in the proof and it is easy to verify . for any @xmath376 , there are only two product states @xmath377 , they also span the space @xmath378 . that is , any state in @xmath378 is the linear combination of @xmath379 and @xmath380 .
it is clear that @xmath381 are compatible with @xmath1 .
suppose they are compatible with another three - qubit state @xmath382 . by applying the observation to system
@xmath383 we have @xmath384 with some complex numbers @xmath385 . by applying the observation to system @xmath182 we have @xmath386 , and hence @xmath387 .
so we may assume @xmath388 .
as a result , the range of the state @xmath277 is spanned by the product states @xmath389 . by simple algebra
one can see that the only feasible @xmath277 compatible to @xmath390 and @xmath58 is the convex sum of @xmath391 and @xmath392 . by using the condition
@xmath393 we obtain @xmath394 .
this completes the proof .
we provide here some further examples of states with separable reductions that are only compatible with genuine multipartite entanglement , also making use of proposition  [ le : entunique ] and lemma [ le : unique ] .
note that @xmath395 for @xmath287 in eq .  .
it follows from lemma  [ le : unique ] that for any @xmath396 , the separable states @xmath397 , @xmath398 , and @xmath399 are uniquely compatible with the state @xmath400 .
so we have generated a family of separable bipartite marginals which are uniquely compatible with a genuinely entangled state , extending example  [ ex : main ] .
we now generalize example  [ ex : main ] to a different family of states that satisfy the conditions in theorem  [ thm : main ] .
let @xmath277 be as in eq .  , and the product state @xmath401 .
such product state always exists because @xmath402 , and there is a product state in any 2-dimensional two - qubit subspace . for example , we can choose @xmath403 .
we have the following corollary now .
[ cr : mixture ] let @xmath404 , @xmath405 , @xmath273 , @xmath406 . for
@xmath407 suppose the product states @xmath408 where @xmath409 is real and @xmath255 is the reduced density operator of the state @xmath277 in eq .  .
the three reduced density operators of the three - qubit state @xmath410 are only compatible with genuinely entangled states .
it is sufficient to show that the three reduced density operators @xmath411 , @xmath412 , and @xmath413 satisfy the three conditions ( i ) , ( ii ) , and ( iii ) in theorem  [ thm : main ] . recall that @xmath277 satisfies these conditions .
since @xmath414 and @xmath415 , we have @xmath416 for any @xmath376 .
so condition ( i ) is satisfied .
next , the same argument shows that @xmath411 has birank @xmath417 , which is exactly condition ( iii ) .
third , the hypothesis @xmath418 and @xmath409 is real imply that the birank of @xmath413 is @xmath419 .
so condition ( ii ) is also satisfied .
k. binder and a. p. young , rev .
phys . * 58 * , 801 ( 1986 ) ; h. t. diep ( ed . ) , _ frustrated spin systems _ ( world scientific , singapore , 2005 ) ; c. m. dawson and m. a. nielsen , phys .
a * 69 * , 052316 ( 2004 )


the determination of scattering amplitudes between pairs of mesons is an interesting problem in strong qcd .
it is also a complicated problem , because both @xmath14 annihilation to @xmath15-channel resonances and `` nonresonant '' scattering are important effects , and it is often difficult to separate the various contributions .
however by specializing to annihilation - free channels such as i=2 @xmath0 and @xmath16 , i=3/2 k@xmath17 , kn and nn , one may study nonresonant scattering in relative isolation .
the determination of resonance parameters , reaction mechanisms , and many other aspects of hadron physics are complicated by the presence of nonresonant scattering , which is treated as an ( often poorly understood ) initial - state and final - state rescattering effect .
developing an accurate description of nonresonant scattering would help clarify many other aspects of hadron physics .
a further interesting possibility is that sufficiently attractive nonresonant scattering may lead to weakly bound hadron - hadron or multihadron states , as does happen in nuclei and hypernuclei .
we may also find a rich spectrum of meson - meson bound states , the study of which will extend nuclear physics into the largely unexplored field of  mesonic nuclei " or  molecules " @xcite .
an understanding of psps , psv and other meson - meson scattering amplitudes is also important for the interpretation of non - qcd processes such as nonleptonic weak decays , since these show evidence of important hadronic final state interactions .
the @xmath18i=1/2 rule is a well known example .
similarly , a recent study of d and d@xmath19 decays to k@xmath20 @xcite found that the dalitz plots are dominated by two - meson isobars , including @xmath21 , k@xmath22 and k@xmath23 , and complex relative amplitudes are required to describe the d@xmath24 dalitz plot . without final state interactions one would expect relatively real couplings to these final states .
one finds a surprising variety of approaches to strong hadron - hadron scattering in the literature .
there are many studies using effective hadronic lagrangians , such as the  chiral perturbation theory " description of the psps sector .
although this method is convenient because it uses perturbative qft techniques , it is incomplete in that it takes effective lagrangian vertex strengths from the data ; one should be able to calculate these hadronic couplings directly from quark - gluon forces .
second , there are studies that model the low energy hadron - hadron scattering mechanism , which include the apparently dissimilar meson exchange and quark - gluon descriptions of hadronic forces .
meson exchange models are again attractive for their simplicity , since they use perturbative qft techniques to determine scattering amplitudes .
this approach has been elaborated in greatest detail in models of the nn force @xcite , in which a large number of meson exchanges is assumed . with this large parameter space
a good description of this interaction is possible , although there is a concern that one may be parametrizing other scattering mechanisms in addition to @xmath25-channel meson exchange .
alternatively , one may calculate hadron - hadron forces directly from the fundamental quark - gluon interaction , using quark model hadron wavefunctions .
this approach has also seen its most detailed development in studies of the nn interaction @xcite , and is most successful in describing the short - ranged repulsive core .
maltman and isgur @xcite also found a physically reasonable intermediate ranged attraction from a color van der waals effect in the quark - gluon approach , which is _ not _ equivalent to the usual @xmath0 or @xmath26 meson exchange explanation of this force .
the quark description of hadron - hadron interactions is complicated by the combinatorics of matrix elements between quark bound states , but has the advantage that it can easily be extended to a wide range of spin and flavor channels through a simple change of the external hadron wavefunctions .
a third promising approach is to infer hadron scattering amplitudes from lgt .
to date lgt has seen little application to scattering problems because of the difficulty of treating systems that are not in their ground states .
estimates of the i=0 and i=2 @xmath0 scattering lengths have been obtained by exploiting a theoretical relation to finite - size effects @xcite , and more recently very interesting results for nuclear physics potentials in the @xmath27 system were reported @xcite . in future it may be possible to improve hadron scattering models through comparisons with similar `` lgt data '' .
in this paper we are concerned with the derivation of meson - meson scattering amplitudes from quark - gluon forces .
we derive meson - meson scattering amplitudes at lowest order in the quark - quark interaction , which leads to a quark interchange model described by `` quark born diagrams '' @xcite .
since the quark - quark interaction is already well established from hadron spectroscopy , our predictions have little parameter freedom .
in previous work we and others ( usually assuming oge hyperfine dominance ) have shown that this approach gives a reasonably accurate description of s - wave scattering in a wide range of channels without @xmath14 annihilation , including i=2 @xmath0 @xcite , i=3/2 k@xmath17@xcite , i=0,1 kn@xcite , i=0,1 bb@xcite ( compared to lgt data ) , and the nn repulsive core @xcite .
this approach has also been applied to @xmath28 @xcite and other reactions relevant to heavy ion collisions , where the experimental low energy cross sections are as yet unclear .
the principal new contribution of this paper is a detailed analytical derivation of the meson - meson scattering amplitudes that follow from the complete quark - quark interaction , including color coulomb , linear scalar confinement , oge spin - spin , oge spin - orbit , oge tensor and linear spin - orbit forces . as a future application of these results , one might hope to clarify the relationship between meson exchange and quark interchange models by a detailed comparison of the spin dependence of hadron - hadron scattering amplitudes , which we expect to be sensitive to the details of the scattering mechanism .    here
we consider both pseudoscalar - pseudoscalar ( psps ) and pseudoscalar - vector ( psv ) scattering .
the former is a `` standard benchmark '' for meson scattering models , because i=2 @xmath0 low energy scattering has no @xmath15-channel resonances and has been the subject of many experimental phase shift analyses .
although we find reasonable agreement with s - wave i=2 @xmath0 scattering , this channel has no spin degree of freedom , and so can not be used to test the characteristic spin dependences predicted by the quark model s oge and linear scalar confinement forces .
we find in contrast that psv is an excellent theoretical laboratory for the study of spin - dependent forces , as it can accommodate both meson - meson spin - orbit and tensor interactions .
the spin - dependent forces at the meson - meson level are closely related to the corresponding terms in the quark - quark interaction in our approach .
although the study of psv scattering is essentially a theoretical exercise at present , these phase shifts are accessible experimentally , for example through measurement of the relative s and d final state phases in @xmath29 . thus it should be possible to measure psv phase shifts from resonance decays to multiamplitude psv final states .
before we proceed to our detailed results , we note that some work has already appeared on meson - meson scattering in psv systems .
numerical results for many light s - wave psv meson channels were previously reported by swanson @xcite using a similar quark model approach that incorporated oge spin - spin and linear confinement forces .
theoretical results for psv scattering ( @xmath16 in particular ) in a meson exchange model were published by janssen _
_ @xcite and bckmann _ et al . _
@xcite , assuming @xmath17 , vector and @xmath30 exchange .
since the @xmath31 , @xmath32 and @xmath33 vertex strengths are relatively well established , it was possible to evaluate these scattering amplitudes numerically .
these papers did not consider the exotic i=2 channel , so a direct comparison with our quark model psv results is not possible at present .
we approximate the full hadron - hadron scattering amplitude by a single ( born - order ) matrix element of the quark - quark interaction hamiltonian @xmath34 .
since @xmath34 is @xmath35 in color , one must then have quark line rearrangement to have a nonvanishing overlap with two color - singlet mesons in the final state . in @xmath36 scattering
there are four independent born - order diagrams , which we label according to which pair of constituents interacted ; these are  transfer@xmath37 " ( t1 ) ,  transfer@xmath38 " ( t2 ) ,  capture@xmath37 " ( c1 ) and  capture@xmath38 " ( c2 ) , which are shown in fig.1 . in the special case of identical quarks
_ and _ identical antiquarks , which is relevant here , there is a second set of four  symmetrizing " diagrams t1@xmath39 c2@xmath40 , which are identical to t1@xmath41c2 except that the quark lines are interchanged rather than the antiquark lines .
the hadron - hadron t - matrix element @xmath42 for each diagram can conveniently be written as an overlap integral of the meson wavefunctions times the underlying quark @xmath42 .
these overlap integrals ( specializing ref.@xcite to the case of equal quark and antiquark masses ) are    @xmath43    the quark @xmath42 has momentum arguments @xmath44 , which are defined in fig.2 . in this paper
we will evaluate these overlap integrals with standard gaussian quark model wavefunctions and the quark @xmath42 for the complete set of oge color coulomb , linear scalar , oge spin - spin , oge and linear scalar confinement spin - orbit and oge tensor interactions .
these interactions are given in app.a .
we specialize the general problem of psps scattering without @xmath14 annihilation to i=2 @xmath0 because many experiments have published phase shift analyses of this channel .
the other @xmath0 channels have large @xmath15-channel @xmath14 annihilation contributions .
the full i=2 @xmath0 born - order t - matrix element is determined by adding the individual contributions of app.b , with psps spin matrix elements given in app.c , part c2 .
there are also flavor and color factors for each diagram and an overall  signature " phase of @xmath45 , and a second set of  symmetrizing " diagrams for identical quarks and identical antiquarks , as discussed in detail in ref.@xcite .
on summing these contributions we find    @xmath46    where @xmath47 is an abbreviation for the confluent hypergeometric ( kummer ) function @xmath48 .
the three separate expressions above are the oge spin - spin , color coulomb and linear confinement contributions respectively .
the @xmath49 terms come from the transfer diagrams , and the remaining , isotropic , terms come from the capture diagrams .
the spin matrix elements of the spin - orbit and tensor terms vanish identically in the psps channel .    since @xmath50 and @xmath51 , one can equivalently write this amplitude as a function of the cm momentum and scattering angle using @xmath52 , where @xmath53 .
the bose symmetry required for this @xmath0 scattering amplitude is evident .
we may derive the elastic born - order i=2 @xmath0 phase shifts from eq.([pipit ] ) , using the relation between phase shifts and the t - matrix given in app.d , especially eq.(d17 ) , and the integrals in app.g .
the result we find for the s - wave is @xmath54 where we have introduced @xmath55 .
the total born - order s - wave phase shift is the sum of these three contributions .
this s - wave phase shift is shown in fig.3 with our standard quark model parameter set @xmath56 , @xmath57  gev , @xmath58  gev and @xmath59  gev@xmath60 .
we also use @xmath61  gev throughout .
this confirms that the color coulomb and linear confinement interactions make relatively small contributions to the i=2 @xmath0 s - wave at moderate energies .
the weakly repulsive linear confining interaction in i=2 @xmath0 near threshold was previously found numerically by swanson@xcite .
one might be concerned about the approximation of using sho wavefunctions , especially at higher energy scales where there should be strong short - distance components in the pion wavefunction due to the attractive spin - spin hyperfine interaction . to test the sensitivity to sho wavefunctions we evaluated the i=2 @xmath0 scattering amplitudes and phase shifts numerically using coulomb plus linear plus hyperfine @xmath14 wavefunctions and monte carlo integration of the real space integrals corresponding to the t - matrix integrals ( [ t1def]-[c2def ] ) . as usual
this requires a `` smearing '' of the contact hyperfine term , @xmath62 , to allow solution of the schrdinger equation with an attractive delta - function interaction . in the literature the inverse smearing length is typically taken to be @xmath63  gev .
( a calculation of i=2 @xmath0 scattering with this interaction and @xmath64  gev was reported previously by swanson @xcite . ) with our standard light - quark parameter set @xmath56 , @xmath65  gev@xmath60 and @xmath66  gev , we found that fitting the @xmath67 splitting required a value of @xmath68  gev . to illustrate the dependence of the scattering amplitude on this parameter , in fig.4 we show the i=2 @xmath0 s - wave that follows from our standard quark model set @xmath69 , with @xmath70 and @xmath71  gev .
clearly the predicted phase shift is rather similar to the sho result of fig.3 , although the effect of short - distance peaking in the @xmath17 wavefunction is evident above @xmath72  gev .
we also show most of the higher statistics experimental results for the i=2 @xmath0 s - wave phase shift in figs.3 and 4 .
the references shown are colton _ et al._@xcite , durusoy _ et al._@xcite ( ope extrapolation , solid ; ope + dp form factor , open , slightly displaced in @xmath73 for visibiity ) , hoogland _ et al._@xcite ( extrapolation b ) , and losty _ et al._@xcite .
et al._@xcite found a wide range of results from three different off - shell extrapolations , so we simply quote their fitted scattering length below .
clearly there is already reasonable agreement with the experimental s - wave phase shift at lower energies without fitting the quark model parameters .
the model predicts a rather dramatic extremum in this phase shift near @xmath74  gev , which is unfortunately beyond the limiting invariant mass of most of the experimental studies .
there are some measurements of this phase shift at higher invariant mass with lower accuracy due to durusoy _
et al._@xcite , which are also shown in the figure .
the durusoy _
_ results support our predicted extremum near @xmath74  gev ; indeed , their phase shift above @xmath74  gev appears to fall even more rapidly than we predict .
we have investigated optimal parameter fits of the s - wave phase shift formula eq.(6 ) to the data , but we find that these are rather unstable because the color coulomb and linear confinement contributions are small and are qualitatively similar functions . in any case
the durusoy _
data and fig.4 show that the hyperfine smearing distance @xmath26 is an important parameter , and this will not be well determined until accurate phase shift data becomes available at higher invariant mass .
an accurate measurement of i=2 @xmath0 scattering amplitudes near and above @xmath74  gev would clearly be very useful as a test of this and other models of meson - meson scattering .
the i=2 @xmath0 scattering length is defined by @xmath75 .
the results we find from eq.(6 ) are @xmath76 and their numerical values with our standard quark model parameters set are @xmath77 } & { $ s\cdot s$ }    \cr   -0.007 \ { \rm [ fm ] } & { cou . }
\cr   -0.017
\ { \rm [ fm ] } & { lin . }    \cr   -0.109 \
{ \rm [ fm ] } & { \rm total}. } \ ] ] the coulomb and linear contributions were independently checked by monte carlo integration of the corresponding real - space overlap integrals .
the relative sizes of these numerical contributions _ a posteriori _ justify the approximation of neglecting the color coulomb and linear terms in i=2 @xmath0 scattering .
the i=2 @xmath0 scattering length has been calculated previously using many other theoretical approaches .
a summary of some of these predictions is given below .
( we use a current value of @xmath78  mev in weinberg s pcac formula @xmath79 . )
@xmath80 } & { lgt \cite{lgt_pipi } }   \cr   -0.052   \ { \rm [ fm ] } & { meson exchange \cite{mesex , fs } } \cr   -0.053   \ { \rm [ fm ] } & { roy eqs.\cite{ana } }    \cr   -0.063    \ { \rm [ fm ] } & { pcac \cite{wein}}.    } \ ] ]    although @xmath81 and other effective lagrangian parameters are normally taken from experiment , these parameters are of course calculable from quark - gluon forces . as an example
, our result for @xmath82 yields the following expression for @xmath81 , @xmath83 the dominant contribution comes from the @xmath84 oge @xmath85 term .
experimental determinations of the scattering length have yielded results which are larger than theoretical expectations ; @xmath86 } & { losty { \it et al.}\cite{losty } }    \cr   -0.24(2 ) , -0.22{+0.03\atop -0.04 } \
{ \rm [ fm ] } & { prukop { \it et al . } \cite{prukop } }   \ . } \ ] ] we speculate that this discrepancy is due to the use of a simple @xmath87 effective range formula in extrapolation .
the difficulty of extrapolating experimental phase shifts to threshold has been stressed by morgan and pennington @xcite .
we advocate the use of a  generalized specific heat plot " of low energy phase shifts for this purpose @xcite .
this plot takes into account the threshold behavior seen in eq.(6 ) , @xmath88 where @xmath89 is a relatively slowly varying function of @xmath90 .
thus the threshold behavior is approximately proportional to @xmath91 rather than just @xmath92 , and since @xmath93 is quite small this leads to rapid variation near threshold and makes the linear-@xmath92 extrapolation inaccurate .
we suggest removal of all this dependence by displaying @xmath94 versus @xmath1 .
the intercept in this plot is the scattering length , and the slope at intercept implies the effective range .
this generalized specific heat plot is shown in fig.5 for i=2 @xmath0 scattering .
an extrapolation of the moderate-@xmath92 data can now be seen to be much closer to the theoretical scattering lengths .
the small-@xmath92 dependence of @xmath95 was calculated by donoghue @xcite in a chiral effective lagrangian , which gave the weinberg result at @xmath96 and an @xmath97 correction factor of @xmath98 .
this is just the correction due to an overall factor of @xmath99 , so this predicts a zero slope in @xmath1 for @xmath94 at threshold .    the jlich meson - exchange model @xcite , which is dominated by @xmath25-channel @xmath100 exchange in this process , also predicts rapid variation in @xmath95 near threshold .
the prediction of this model for @xmath101 @xcite , shown in fig.5 , is rather similar to our quark model result .
low energy `` phase shift equivalent '' gaussian i=2 @xmath0 potentials , derived using the method of mott and massey @xcite as described in app.e , are given below .
we quote separate gaussians for the transfer and capture contributions from each of the three interactions , spin - spin contact , color coulomb and linear hyperfine .
however their predicted phase shift decays more slowly at large momentum , probably due to the use of power law form factors in their vertices .
@xmath102    in ref.@xcite we derived i=2 @xmath0 potentials for the spin - spin contact interaction using the  locality expansion " method of ref.@xcite .
this gave an identical result for the spin - spin transfer potential , because this amplitude ( before bose symmetrization ) is a function of @xmath25 only .
however for the capture diagrams the mott - massey approach used here gives a different potential , since it is constrained to reproduce the @xmath103 series expansion of the phase shift in eq.(6 ) , but the local approximation is not .
the two capture potentials reproduce the scattering length , but the local approximation gives an incorrect effective range .    the low energy mott - massey i=2 @xmath0 potential is shown in fig.6 for our standard quark model parameters @xmath56 , @xmath65  gev@xmath60 and @xmath66  gev .
the spin - spin hyperfine contribution is dominant over the range shown .
the higher partial waves ( l @xmath105 ) may be evaluated similarly . according to eq.([pipit ] ) , these receive contributions only from the transfer diagrams . the born - order d - wave phase shift with sho wavefunctions is given by    @xmath106    these three expressions are numerically rather small , and their phases are such that they approximately cancel ; at @xmath74  gev they are respectively @xmath107 , @xmath108 and @xmath109 . to see this more clearly , the leading @xmath110 behavior predicted by eq.(14 ) is @xmath111 and the three dimensionless combinations @xmath112 , @xmath113 and @xmath114 are comparable in size .
we have also evaluated this d - wave phase shift using coulomb plus linear plus hyperfine wavefunctions .
the result is shown in fig.4 , and is numerically similar to the sho d - wave , eq.(14 ) .    in comparison the experimental d - waves reported by durusoy _
_ @xcite and hoogland _ et al .
_ @xcite are @xmath115 at @xmath74  gev ( see fig.4 ) .
( losty _ et al _ @xcite report a rather larger but inconsistent low - energy d - wave . )
this is clearly larger than our prediction , although the rather slow variation of the durusoy _ et al . _ and hoogland _
d - waves with @xmath116 may indicate a problem with the measurements ; the expected threshold behavior of @xmath117 is much more rapid than the observed energy dependence . unfortunately the dispersion relations represented by the roy equations have technical difficulties with determining d- and higher waves @xcite .
they do however lead to predictions of a _ positive _ d - wave close to threshold , which is not evident in the data .
the d - wave may well have important meson exchange contributions , since this type of model can accommodate the reported experimental phase shift @xcite .
for simplicitly we will initially quote results only for i=2 @xmath16 .
the other isospin channels are simply related by flavor factors , which we will discuss subsequently .
we assume identical spatial wavefunctions , so only the @xmath100 spin degree of freedom and difference in phase space distinguish this case from @xmath0 .
summing the individual contributions in app.b with the appropriate flavor and color factors and the @xmath45 signature phase , and using the psv spin matrix elements of app.c part c3 , we find for the i=2 @xmath16 t - matrix    @xmath118    @xmath119    @xmath120    @xmath121    @xmath122 \hfill\ ] ]    @xmath123 \hfill\ ] ]    @xmath124   \hfill\ ] ]    @xmath125 + \bigg [ \vec s_\rho \cdot \vec c\;\vec
s_\rho \cdot \vec c     - { 2\over 3}\vec c^{\ , 2 } \bigg ] \bigg ] \hfill \ .\ ] ]    @xmath126    the individual contributions in this result are respectively oge spin - spin ; oge color coulomb ; linear confinement ; oge spin - orbit ; linear scalar confinement spin - orbit ; oge tensor ( transfer diagrams ) ; and oge tensor ( capture diagrams ) . in all these we list transfer followed by capture contributions .
@xmath127 is the @xmath100 meson spin vector , @xmath128 and @xmath129 are the initial and final @xmath17 momenta , @xmath130 and @xmath131 are the initial and final @xmath100 momenta , and @xmath132 as in @xmath0 . since this result was derived in the cm frame , @xmath133 .
this @xmath42 evidently describes @xmath16 spin - orbit and tensor interactions , in addition to spin - independent scattering .
it is interesting that there is a one - to - one mapping between the quark - quark spin - orbit and tensor interactions and these @xmath16 spin - orbit and tensor terms .
this simple result need not be true in general ; a given spin - dependent interaction at the quark level may give rise to a different type of hadron - hadron interaction .
as an example , a mapping of a tensor nucleon - nucleon force into a nucleon - nucleus spin - orbit interaction was discussed by stancu , brink and flocard.@xcite    to evaluate phase shifts and inelasticities it is convenient to calculate the matrix element of our @xmath16 t - matrix eq.([pirhot ] ) between general @xmath134 states , which gives the reduced matrix element    @xmath135    @xmath136    as discussed in app.d .
this is a straightforward exercise , although integrals of special functions and a careful angular analysis of the spin - orbit and tensor terms are required ; the details are discussed in apps.g and h. this matrix element is diagonal in @xmath137 except for the tensor interaction , which has both diagonal and off - diagonal ( transfer ) and fully off - diagonal ( capture ) contributions .
the @xmath137-diagonal results , again showing transfer diagram contributions followed by capture , are    @xmath138    and the off - diagonal tensor matrix elements are 0.5 cm @xmath139^{1/2 } \over   ( 2j+1)\hskip0.3 cm   } \ ; x \bigg ( { \cal f}^{(j-1)}_{{5\over 2},{7\over 2}}(x )   + 2 { \cal f}^{(j)}_{{5\over 2},{7\over 2}}(x )   + { \cal f}^{(j+1)}_{{5\over 2},{7\over 2}}(x )   \bigg)\ ; \ ; e^{-2x }   \\ & & \nonumber \\ & & t^j_{l'\ne l}|^{\it capture}_{oge\ t }   =    { \pi^2 \alpha_s \over m^2 } \ ;   \delta_{j1 } \ ; \big ( \delta_{l2 } \delta_{l'0 }   + \delta_{l0 } \delta_{l'2 }   \big )   \ ; { 2^{19/2 } \over 3^{11/2}\cdot 5}\ ; x\ ;   f_{{5\over 2},{7\over 2}}(2x/3 ) \ ; e^{-2x } \ .\end{aligned}\ ] ] 0.5 cm in these formulas @xmath140 is a modified spherical bessel function , the tensor @xmath141 matrix element between @xmath142 @xmath16 states is @xmath143 and the integral @xmath144 is evaluated in app.g .      in s - wave to s - wave scattering the spin - orbit and tensor @xmath16 t - matrix contributions vanish , and we are left with color coulomb , linear and spin - spin contributions , just as in the i=2 @xmath0 case .
the i=2 @xmath16 s - wave phase shifts that result from these interactions , again using eq.(d17 ) , are @xmath145 where @xmath146 , and again @xmath55 . in fig.6
we show these individual components and the total s - wave phase shift with our standard quark model parameter set and meson masses ( used throughout ) of @xmath147  gev and @xmath148  gev .
the forces considered here evidently lead to strong repulsion in the i=2 @xmath16 channel .
the spin - orbit and tensor terms in eqs.(21 - 23 ) all contribute to @xmath150 @xmath16 scattering , and there is also an odd-@xmath137 , @xmath151-independent term due to the oge spin - spin interaction in eq.(18 ) , which is not symmetric under @xmath152 .
the color coulomb and linear confinement spin - independent terms , eqs.(19,20 ) , contribute only to even @xmath137 .
adding the various diagonal matrix elements of eqs.(18 - 23 ) and using eq.(d17 ) gives phase shifts for each @xmath153l@xmath154 partial wave .
in fig.8 we show results for all p - wave channels and for j = l@xmath155 in d- and f - wave .
note that there is a large , inverted spin - orbit force in the p - wave , so the @xmath153p@xmath156 phase shift is widely separated from @xmath153p@xmath38 , and has an even larger maximum phase shift than the s - wave .
the higher - l channels show decreasing phase shifts with increasing l , as expected for short - ranged quark - gluon forces .
the relative importance of the individual contributions to the spin - dependent force is of considerable interest . in fig.9
we show the various spin - dependent contributions to the i=2 @xmath153p@xmath38 @xmath16 phase shift .
the largest contribution arises from oge spin - orbit , in particular from the transfer diagrams .
the oge and confinement spin - orbit capture diagrams give smaller contributions of the same sign .
finally , the confinement spin - orbit transfer diagrams have a sign opposite to all these and reduces the total spin - orbit force somewhat .
this dominance of the psv spin - orbit by oge is an interesting result , especially since mukhopadhyay and pirner @xcite found the opposite result in kn . in that system
they concluded that confinement , not oge , makes the largest contribution to the spin - orbit force .
the oge tensor in i=2 @xmath153p@xmath38 @xmath16 is weakly repulsive ; it makes a much larger contribution to @xmath153p@xmath37 and @xmath153p@xmath156 , where the tensor matrix element is respectively @xmath157 and @xmath158 times as large .
the oge tensor is evident in fig.8 , in the departure of the ratio ( @xmath153p@xmath159p@xmath37):(@xmath153p@xmath160p@xmath156 ) from the pure spin - orbit value of 2:1 at higher energies .
there is also an off - diagonal coupling due to the oge tensor terms , given by eqs.(24,25 ) , but we have neglected this in calculating phase shifts because we find that it is numerically a small effect . the largest coupling at low energies is @xmath161s@xmath162d@xmath37 , which leads to an inelasticity of only @xmath163 by @xmath164  gev ( calculated using eqs.(d18-d20 ) ) .
we may determine low energy gaussian equivalent @xmath16 potentials from the phase shifts , as discussed in app.e .
the most interesting potential phenomenologically is the spin - orbit one , since the origin of the spin - orbit interaction in the nn system is a long - standing and still poorly understood problem .
in particular we derived gaussian potentials corresponding to the p - wave phase shifts due to the oge and linear scalar confinement spin - orbit interactions , using eq.(e4 ) of app.e .
the results for the transfer and capture contributions to these potentials are    @xmath165    the oge , linear and total spin - orbit potentials for the @xmath153p@xmath38 wave of i=2 @xmath16 are shown in fig.10 for our standard parameter set . the largest contribution to
the @xmath16 spin - orbit force comes from oge transfer diagrams ; the linear confinement spin - orbit from the transfer diagrams is about half as large and opposite in sign , and the two capture diagram contributions are much smaller .
since the confinement capture and transfer diagrams have opposite signs , the net result is dominance of the psv spin - orbit by the oge contribution .      0.5 cm    .
j@xmath3 exotic states in psv . [ cols="^,^,^,^,^,^,^,^,^,^ " , ]     [ table2 ]
inspection of the tables shows that the largest exotic scattering amplitude should be in the i=2 @xmath16 @xmath167 p - wave .
the elastic phase shift in this channel is the @xmath153p@xmath156 curve in fig.8 .
the large negative phase shift shows that this is a strongly repulsive channel ; the maximum phase shift is predicted to be a quite large @xmath168 at @xmath169 gev , which exceeds even the s - wave phase shift maximum .
the largest _ attractive _ exotic phase shift we have found in psv is the i=0 partner , which is @xmath170 of i=2 , giving a maximum phase shift of @xmath171 at the same mass .
we do not find sufficient attraction to form a meson - meson `` molecular '' bound state in any of these nonstrange j@xmath3-exotic psv channels .
the @xmath172 channels are relatively weak because only the @xmath173 part of the @xmath172 contributes to these diagonal scattering amplitudes ; the @xmath174 component leads to open - strange final states ( k@xmath175 for example ) after quark line interchange .    regarding candidate exotic resonances ,
there have been speculations that the determination of the mass and width of the exotic candidate @xmath4 may have been compromised by inelastic rescattering effects @xcite , analogous to the `` deck effect '' proposed as a nonresonant explanation of the @xmath176 .
for example , crossing the @xmath177 threshold at @xmath178  gev in the process @xmath179 might mimic resonant phase motion if this process has a rapidly varying inelasticity .
we can test this and other nonresonant possibilities by calculating the elementary @xmath180 scattering amplitudes using our quark model approach .
some important results follow from simple flavor factors .
note in particular that the nonresonant scattering amplitude @xmath181 vanishes in any i=1 channel , including the @xmath182 exotic one .
this is a general result whenever the quark line diagram of fig.1 dominates ; clearly a pair of oppositely charged , nonstrange @xmath14 mesons @xmath183 can not scatter into another charged pair @xmath184 under quark interchange .
a comparison with isospin matrix elements shows that this implies that scattering of any two @xmath14 isovectors in i=1 vanishes .
this isospin selection rule eliminates two subprocesses discussed by donnachie and page @xcite as deck effect backgrounds that might shift a higher - mass exotic resonance to an apparent @xmath4 , @xmath185 and @xmath186 .
independent of any scattering model , one should note that the coupling @xmath187 is probably small because of the strong ves experimental limit ( reported by v.dorofeev @xcite ) of @xmath188 since @xmath189 is a large mode ( @xmath190 @xcite ) , if @xmath187 rescattering were important we would also expect to observe a large @xmath191 branching fraction .
we also expect the final background process suggested by donnachie and page ( @xmath192 ) to be small , because the direct time ordering @xmath193 vanishes in p - wave in this model
. finally , the rescattering process they propose , @xmath194 , does not vanish in the quark interchange model , although the required @xmath195 and suppressed @xmath172 flavor factor may nonetheless make this a relatively weak amplitude .
a calculation of this and related scattering amplitudes is planned for a future publication .
although there is little experimental information about psv interactions at present , these phase shifts actually are experimentally accessible in existing data , for example as relative fsi phases in the d and s amplitudes in @xmath196 .
these are usually , and incorrectly , taken to be relatively real amplitudes .
the relative phase including the fsi is @xmath197 @xcite , and is observable for example as a reduction in the strength of the sd cross term in the @xmath198 angular distribution by @xmath199 . since
this method requires individual measurements of the s@xmath60 , d@xmath60 and sd cross term in the angular distribution , it should be applicable to cases such as @xmath200 and @xmath201 where s and d are of comparable magnitude @xcite .
the @xmath202 phases we predict at these masses ( which are calculated from + 1/2 times the @xmath153s@xmath37 and @xmath153d@xmath37 i=2 @xmath16 phases in figs.7,8 ) are @xmath203 and @xmath204 .
this proposed technique is similar to that used in k@xmath205 decays @xcite , in which the low energy i=0 @xmath0 s - wave phase shift is actually observed as the difference between the i=0 s - wave and i=1 p - wave @xmath0 fsi phases .
in this paper we have derived meson - meson scattering amplitudes , including spin - dependent forces , from a calculation of the born - order matrix element of the quark - quark interaction between two - meson states . since @xmath14 annihilation
is not included in the model , it describes scattering that does not involve coupling to @xmath15-channel resonances .
this includes for example i=2 and the nonresonant backgrounds in all channels , including exotic j@xmath3 .
we considered the cases of psps and psv scattering , and derived the scattering amplitudes in all @xmath206 channels for these cases .
the parameters of the model were previously fixed by quark model studies of hadron spectroscopy .
where possible we have compared the results to experiment .    in i=2 @xmath0 ( the best studied psps case )
the results were shown to be in reasonable agreement with experiment in s - wave scattering , and an extremum predicted near @xmath74  gev is supported by the data .
rapid variation of @xmath95 is predicted near threshold , which may reconcile theoretical expectations of a small scattering length with larger reported experimental values based on extrapolation in @xmath1 . the experimental d - wave , although quite small , is clearly larger than the model predicts .
the psv system is a convenient theoretical laboratory for studying spin - dependent forces , since it can accommodate both spin - orbit and tensor interactions , and is simpler than kn or nn . we derived analytical results for these spin - dependent psv interactions ( t - matrices and phase shifts ) given sho wavefunctions and the standard spin - dependent quark model forces .
the quark - quark spin - orbit and tensor forces map directly into spin - orbit and tensor psv interactions .
we find that the oge spin - orbit force in the psv system is quite large in p - wave , and so is expected to be large in many other hadron - hadron systems as well .
there is no psv phase shift data at present .
we noted however that psv phase shifts actually can be measured in multiamplitude resonance decays to psv final states , so it should be possible to test theoretical predictions for psv scattering amplitudes in future experimental studies .
our predictions for scattering in j@xmath3-exotic channels are of current interest because the reported exotics might be complicated by large and rapidly varying nonresonant inelasticities .
one speculation is that the @xmath4 parameters might be strongly affected by the opening of inelastic couplings to the @xmath177 channel . in our model
( and in any @xmath14 constituent interchange model ) several of these nonresonant processes can be rejected as significant complications because of vanishing flavor factors .    in future
we plan to extend our calculations to other exotic meson - meson channels , such as s+p , to test whether strong attractive interactions are predicted that might support
 multiquark exotics " such as s+p molecules .
we also plan to apply the current approach to the study of spin - dependent interactions in other hadronic systems , including kn , nn and light hadron + charmonium systems .
we are grateful to our colleagues d.blaschke , n.cason , h.g.dosch , s.krewald , e.klempt , n.isgur , c.michael , m.r.pennington , b.pick , h .- j.pirner , g.rpke , f.sassen , j.speth , fl.stancu and c .-
y.wong for useful discussions of various issues relating to meson - meson scattering .
this research was supported in part by the doe division of nuclear physics , at ornl , managed by ut - battelle , llc , for the us department of energy under contract no .
de - ac05 - 00or22725 .
es acknowledges support from the doe under grant de - fg02 - 96er40944 and doe contract de - ac05 - 84er40150 under which the southeastern universities research association operates the thomas jefferson national accelerator facility .
tb acknowledges additional support from the deutsche forschungsgemeinschaft dfg at the university of bonn and the forschungszentrum jlich under contract bo 56/153 - 1 .
the various contributions to the quark - quark @xmath42 ( with color factors of @xmath207 removed ) are + @xmath208   & { \rm oge spin - spin}\cr   + { 4\pi \alpha_s \over    \vec q^{\ , 2 } } & { \rm oge color coulomb}\cr   + { 6\pi b \over \vec q^{\ , 4 } } & { \rm linear conft . } \cr   + { 4i\pi \alpha_s \over   m^2\vec q^{\ , 2}}\ ; \bigg [ \vec s_1\cdot \big ( \vec q\times ( { \vec p_1\over 2}-\vec p_2)\big ) + \vec s_2\cdot \big(\vec q\times ( \vec p_1-{\vec p_2\over 2 } ) \big )   \bigg ] & { \rm oge spin - orbit}\cr   - { 3i\pi   b \over m^2\vec q^{\ , 4}}\ ; \bigg [ \vec s_1\cdot ( \vec q\times \vec p_1 ) - \vec s_2\cdot ( \vec q\times \vec p_2 ) \bigg ] & { \rm linear spin - orbit } \cr + { 4\pi \alpha_s \over m^2\vec q^{\ ; 2}}\ ; \bigg [ \vec s_1 \cdot \vec q \ ; \vec s_2 \cdot \vec q
-{1\over 3}\ ; \vec q^{\ ; 2}\ , \vec s_1\cdot \vec s_2 \bigg ] & { \rm oge tensor.}\cr   }   \label{quarkt}\ ] ]    the standard @xmath14 quark model gaussian wavefunction is given by @xmath209 where in general @xmath210 and for our special case of equal quark and antiquark masses @xmath211
in this appendix we give the explicit meson - meson t - matrix elements that follow from the overlap integrals eqs.([t1def]-[c2def ] ) with gaussian wavefunctions and the various quark t - matrix elements .
the oge spin - spin hyperfine , color coulomb and linear confinement results were derived previously @xcite . for completeness
we quote the formulas here , as well as giving the new spin - orbit and tensor results .
the multiplicative diagram - dependent color and flavor factors and the signature phase ( which is @xmath45 for these meson - meson scattering diagrams ) are not included in the results given below .
these formulas abbreviate the confluent hypergeometric function @xmath48 as @xmath47 , and @xmath212 .
these simple contact matrix elements were evaluated previously , for example in ref .
@xcite ( in an equivalent form , but incorporating color factors and the signature phase , as eqs.(71 - 73 ) of that reference ) .
the results are @xmath213 \\ & & { t_{fi}}^{\rm ( t2 ) } =   { t_{fi}}^{\rm ( t1)}(\vec c \rightarrow -\vec c )   \\ \rule{0mm}{1 cm } & & { t_{fi}}^{\rm ( c1 ) } =   -\ ; { 2 ^ 6 \over 3^{5/2}}\ ; { \pi\alpha_s \over m^2}\ ; e^{-\vec a^{\ , 2}/3\beta^2 } \bigg [ \vec s_1\cdot \vec s_2 \bigg ] \\ & & { t_{fi}}^{\rm ( c2 ) } =   { t_{fi}}^{\rm ( c1 ) }   \ . \end{aligned}\ ] ]      the contribution of the oge color coulomb interaction to the meson - meson t - matrix follows from the evaluation of the integrals eqs.([t1def]-[c2def ] ) with the second quark - quark @xmath42 in eq.([quarkt ] ) .
the results are    @xmath214      the linear confinement integrals were carried out in coordinate space , since the fourier transform of the linear potential is singular .
the results are @xmath215 one may also obtain these results using the momentum space integrals eqs.([t1def]-[c2def ] ) , but the @xmath216 quark - quark @xmath42 in eq.([quarkt ] ) must include a long - distance regularization in the intermediate stages of the integration . the final result is well defined due to the gaussian damping provided by the wavefunctions .
we have checked both the linear and color coulomb @xmath42 results by comparing the expressions eqs.([b5]-[b12 ] ) with monte carlo evaluations of the corresponding real - space overlap integrals .
the four oge spin - orbit overlap integrals can be evaluated similarly using the fourth quark - quark @xmath42 in eq.([quarkt ] ) , which gives    @xmath217 \\ & & { t_{fi}}^{\rm ( t2 ) } =   { t_{fi}}^{\rm ( t1)}(\vec c \rightarrow -\vec c )   \\ \rule{0mm}{1 cm } & & { t_{fi}}^{\rm ( c1 ) } =   + { 4 \over 3^{5/2 } } \ ; { \pi\alpha_s\over m^2\beta^2}\ ; { \rm f}_{{3\over 2},{5\over 2}}(\vec a^{\ , 2 } /6\beta^2 )   \ ; e^{-\vec a^{\ , 2}/2\beta^2}\ ; \bigg [ ( \vec s_1 - \vec s_2)\cdot i(\vec a \times \vec c ) \bigg ] \\ & & { t_{fi}}^{\rm ( c2 ) } =   -\ ; { t_{fi}}^{\rm ( c1 ) }   \ . \end{aligned}\ ] ]      the matrix elements eqs.([t1def]-[c2def ] ) , of the scalar confinement spin - orbit interaction in eq.([quarkt ] ) are    @xmath218 \\ & & { t_{fi}}^{\rm ( t2 ) } =   { t_{fi}}^{\rm ( t1)}(\vec c \rightarrow -\vec c )   \\ \rule{0mm}{1 cm } & & { t_{fi}}^{\rm ( c1 ) } =   + { 1   \over 2\cdot 3^{1/2 } } \ ; { \pi b\over m^2\beta^4}\ ; { \rm f}_{{1\over 2},{5\over 2}}(\vec a^{\ , 2 } /6\beta^2 )   \ ; e^{-\vec a^{\ , 2}/2\beta^2}\ ; \bigg [ ( \vec s_1 - \vec s_2)\cdot i ( \vec a \times \vec c ) \bigg ] \\ & & { t_{fi}}^{\rm ( c2 ) } =   -\ ; { t_{fi}}^{\rm ( c1 ) }   \ . \end{aligned}\ ] ]      finally , for the oge tensor @xmath42 ( the last entry in eq.([quarkt ] ) we find    @xmath219 \\ & & { t_{fi}}^{\rm ( t2 ) } =   { t_{fi}}^{\rm ( t1)}(\vec c \rightarrow -\vec c )   \\ & & { t_{fi}}^{\rm ( c1 ) } =   + \ ; { 2 ^ 5 \over 3^{5/2}\cdot 5}\ ; { \pi\alpha_s\over m^2\beta^2}\ ; { \rm f}_{{5\over 2},{7\over 2}}(\vec a^{\ , 2 } /6\beta^2)\ ;   e^{-\vec a^{\ , 2}/2\beta^2}\ ; \bigg [ \vec s_1 \cdot \vec a\ ; \vec s_2 \cdot \vec a -{1\over 3}\ , \vec a^{\ , 2}\ ; \vec s_1 \cdot \vec s_2 \bigg ] \label{b23 } \\ & & { t_{fi}}^{\rm ( c2 ) } =   { t_{fi}}^{\rm ( c1 ) }   \label{b24 } \ . \end{aligned}\ ] ]    the tensor matrix elements in the capture diagrams , eqs.([b23],[b24 ] ) , are the only cases in which we have found a post - prior discrepancy in these psps and psv scattering amplitudes ; the post forms of these matrix elements involve a tensor in @xmath129 rather than @xmath128 , @xmath220 \ .\ ] ] these capture tensor terms vanish in the psps channel .
they do make a small , off - diagonal contribution to psv scattering , _ albeit _ only in the @xmath161s@xmath221d@xmath37 amplitude .
in these scattering amplitude calculations the matrix elements of spin - dependent quark interactions ( the spin - spin , spin - orbit and tensor forces ) involve matrix elements of linear and bilinear quark spin operators .
since the quark spins are not directly observed , it is useful to replace them by the spins of the external hadrons .
this appendix gives the ( diagram dependent ) mapping from quark spins to hadron spins in the psps and psv cases considered in this paper .
the spin matrix elements we require are @xmath222 , @xmath223 , @xmath224 and @xmath225 between general initial and final psps and psv spin states .
our convention for the diagrams ( fig.1 ) is that mesons a and c are always ps ( _ e.g. _ @xmath17 ) , and b and d are ps or v ( _ e.g. _ @xmath17 or @xmath100 ) .      first , in psps scattering there are no external meson spins , so the quark spin matrix elements are proportional to geometrical tensors such as @xmath226
. the matrix elements by diagram are @xmath227 note that the spin - orbit and tensor terms vanish identically in psps scattering ; this follows from applying eqs.(c2,c3 ) to eq.(a1 ) .      in psv scattering the vector ( _ e.g. _ @xmath100 )
meson spin @xmath228 provides an additional degree of freedom , and the linear and quadratic quark spin matrix elements can be expressed in terms of the @xmath100 spin matrix elements @xmath229 and @xmath230 .
the mapping of quark to meson spins is    @xmath231
since total angular momentum is conserved , the t - matrix is block diagonal in a total angular momentum basis , and can be written as @xmath232 in the special case of spinless scattering these basis states are eigenstates of orbital angular momentum @xmath233 so the t - matrix is given by    @xmath234 between momentum eigenstates , @xmath235 are related by @xmath236 and hence @xmath237 where @xmath238 .
a more familiar quantum - mechanical result follows from fixing the incident direction @xmath239 in eq.([tp ] ) and integrating over final angles with a @xmath240 weight , which gives @xmath241    since we define @xmath242 by @xmath243 it is related to the lorentz invariant @xmath180 scattering amplitude @xmath244 defined by the pdg@xcite ( their eq.(35.8 ) ) by @xmath245 and hence to the cm scattering amplitude @xmath246 ( their eq.(35.48 ) ) by @xmath247 ( here @xmath248 is the cm momentum of any particle , @xmath249 and the relation between a diagonal partial wave amplitude @xmath250 and the phase shift @xmath251 allow us to determine @xmath252 from @xmath42 . for purely elastic scattering , and assuming small phase shifts so that @xmath253 , eqs.([tp]-[al ] ) give @xmath254 this relation was used previously to determine for example k@xmath17 @xcite and @xmath137-diagonal kn @xcite elastic phase shifts . in the case of elastic scattering of identical bosons , such as i=2 @xmath0
, there is an additional factor of two for identical particles @xcite , so the relation between the born - order phase shift and the t - matrix element becomes @xmath255    since the angular integral in eq.([deltal1 ] ) is proportional to the amplitude @xmath256 , we may also write the born - order elastic phase shift ( for distinguishable particles ) directly in terms of the t - matrix , @xmath257    this formula has a straightforward generalization to the case of external hadrons with spin , which we use to evaluate @xmath16 phase shifts and inelasticities . in the case of an @xmath137- and @xmath15-diagonal interaction
this is @xmath258 this is adequate for diagonal forces such as our @xmath16 spin - orbit interactions .
the tensor force however is not @xmath137-diagonal , so for this case we must introduce a more general parametrization .
since the tensor interaction couples @xmath259 channel pairs which have the same @xmath151 but differ by @xmath260 , it leads to @xmath261 @xmath262-matrices .
these can be parametrized as    @xmath263 \ .\ ] ]    in our calculation , both of the phase shifts and the inelasticity @xmath264 are @xmath265 , so to this order we can relate these linearly to the matrix elements of @xmath34 .
the born order phase shift formula ( d17 ) remains valid for both channels , and the inelasticity is @xmath266 the overall phase of @xmath267 is dependent on the state normalizations , but the familiar @xmath268 is unique .
potentials provide a very useful representation of hadron - hadron interactions .
they have a clear intuitive meaning , and can easily be used in the nonrelativistic schrdinger equation in searches for bound states or in coupled channel problems .
unfortunately , one may define hadron - hadron potentials in many different ways .
ideally they should reproduce phase shifts or t - matrix elements , at least in the low energy limit .
the assumption of a unique , purely local potential is in general overly restrictive , as it leads to a scattering amplitude that is a function of @xmath25 only . in general
we find @xmath180 scattering amplitudes that depend on both @xmath15 and @xmath25 .
one approach to this problem is to introduce nonlocal  gradient " corrections to the potential @xcite , which can be expressed for example as @xmath269 terms ; this approach leads to the familiar breit - fermi hamiltonian for one - photon and one - gluon exchange , and was used to define hadron - hadron potentials in our previous work @xcite .
alternatively one may project the scattering amplitude onto a given angular channel @xmath137 so that only @xmath15 dependence remains , and find a local @xmath137-wave potential that describes the scattering in that channel .
this definition of potentials was discussed by mott and massey @xcite , and is equivalent to the definition we shall use here .
this approach was previously used by swanson @xcite to define meson - meson potentials from scattering amplitudes .    the quantum mechanical relation between the phase shift @xmath270 and the radial wavefunction @xmath271 in potential scattering of a mass @xmath272 particle ( which becomes the reduced mass below ) in first born approximation is @xmath273 since our t - matrix elements implicitly determine the elastic scattering phase shifts , for example the @xmath274 @xmath0 s - wave in eq.(6 ) , we can invert this formula for each @xmath137 to determine the corresponding @xmath137-wave local potential @xmath275 . in practice
we find that our phase shifts are sufficiently  hard " at high energies to require singular potentials ; this is presumably an artifact of our approximations , such as assuming a contact spin - spin interaction .
for this reason we do not completely invert the phase shift relation ( e1 ) , and instead simply fit an assumed gaussian form @xmath276 to our theoretical _ low energy _ phase shift .
the two gaussian parameters are determined from the @xmath277 and @xmath103 terms in the expansion of the phase shift near threshold , which in the s - wave case are equated to @xmath278 the generalization to higher @xmath137 , also using eqs.(e1 ) , is @xmath279    in determining the gaussian potentials that correspond to our derived phase shifts such as eqs.(6,14,28 ) , we set the external factors of @xmath99 and @xmath280 equal to @xmath281 and @xmath282 before expanding in @xmath248 .
this corresponds to using nonrelativistic phase space and nonrelativistic external hadron line normalizations in our t - matrix calculations , which we assume is the appropriate choice for the derivation of a nonrelativistic equivalent potential .
the  post - prior discrepancy " is a familiar problem in rearrangement collisions ; the diagrams of fig.1 treat the scattering as due to an interaction between the initial hadrons a and b (  prior " ) , but we could equally well have written the interaction between the two final hadrons c and d (  post " ) . the post t - matrices may be obtained from the prior ones by exchanging the initial and final mesons and transforming the momenta @xmath283 , @xmath284 and @xmath285 .
thus for example the post c1 t - matrix is @xmath286 @xmath287    one may show that the post and prior results for the scattering amplitude are equal provided that the external wavefunctions are eigenfunctions of the hamiltonian @xcite .
swanson @xcite shows an example of convergence of post and prior results for meson - meson scattering amplitudes derived from quark born diagrams as the external wavefunctions approach exact hamiltonian eigenstates .
of course the gaussian wavefunctions we use to derive our analytical results are only approximations to the eigenfunctions of the full oge plus linear hamiltonian , so in general we can expect to find a post - prior discrepancy . in this study
we actually find such a discrepancy only in part of the tensor interaction in psv scattering , for which we take the mean of the two results , @xmath288
it is useful in the partial wave decomposition of scattering amplitudes to expand functions of the sum and difference momentum transfers @xmath290 ( here @xmath291 ) in spherical harmonics , @xmath292 @xmath293 this expansion may be inverted to determine the coefficient functions @xmath294 , @xmath295 @xmath296    many of the scattering amplitudes derived in this paper are proportional to confluent hypergeometric functions in @xmath297 or @xmath298 , and their partial wave decomposition requires the integral of a legendre polynomial times a shifted confluent hypergeometric function .
this integral is given by ( abbreviating @xmath299f@xmath300 as @xmath301 ) @xmath302 where @xmath303 and the pochhammer symbol of negative index is @xmath304 in numerical evaluations it is often useful to transform the confluent hypergeometric function in eq.([teds_excellent_int ] ) to a more rapidly converging negative - argument form , using @xmath305 .
as noted in app.d , when evaluating phase shifts and inelasticities it is useful to determine t - matrix elements between @xmath134 states . in the psv system there are spin - orbit and tensor contributions to the t - matrix , and determination of the @xmath206-basis matrix elements of these terms is a complicated problem in angular analysis . here
we show how these matrix elements may be evaluated .        where @xmath307 .
( additional dependence on the rotational scalar @xmath308 is a trivial modification of this angular decomposition . ) to proceed , we expand @xmath309 in spherical harmonics , as in eq.(g1 ) ; @xmath310 and introduce spherical components for the spin and momentum vectors , @xmath311 @xmath312 and the usual state vector expansion , @xmath313 with these substitutions one may determine the @xmath16 @xmath314 matrix elements ( analogous to the spinless matrix element @xmath315 of eq.(d7 ) ) for the spin - orbit term ( h1 ) .
the result involves a sum over a product of six clebsch - gordon coefficients , and can be written as the product of two wigner @xmath316 symbols and two @xmath317 symbols ,        the constraints of the @xmath316 and @xmath320 symbols force this matrix element to be diagonal in @xmath259 , and imply that the only radial components of @xmath321 in eq.([fexpand ] ) that contribute are @xmath322 .
substitution of the explicit @xmath316 and @xmath317 symbols gives our final result for @xmath323 matrix elements of spin - orbit ( h1 ) type , @xmath324    \over 2 ( 2l+1 ) } \ , { \vec a\ , } ^2   \bigg ( f_{l-1 } - \ , f_{l+1 } \bigg ) \ .\ ] ] this result has the overall @xmath325 dependence that one would expect from a spin - orbit force .        and
@xmath327 \ .\ ] ] both tensor matrix elements have @xmath328 contributions , unlike the other interactions we have considered .
the general results in terms of wigner @xmath329 and @xmath330 symbols are @xmath331 @xmath332 and                  r.machleidt , k.holinde and c.elster , phys .
149 , 1 ( 1987 ) . a recent variant of this approach derives nn scattering amplitudes from meson exchange at the quark level ; see d.bartz and fl.stancu , hep - ph/0006012 .
there is an extensive literature on the derivation of the nn interaction from quark - gluon forces in the quark model , beginning with d.a.liberman , phys .
d16 , 1542 ( 1977 ) .
most studies have employed resonating group methods ( see for example c.s.warke and r.shanker , phys .
c21 , 2643 ( 1980 ) ; m.oka and k.yazaki , phys .
90b , 41 ( 1980 ) ; j.e.t.ribeiro , z.phys.c5 , 27 ( 1980 ) ) and variational methods ( for example k.maltman and n.isgur , phys .
d29 , 952 ( 1984 ) ) . the born - order quark - gluon approach described here
was applied to nn by t.barnes , s.capstick , m.d.kovarik and e.s.swanson , phys .
c48 , 539 ( 1993 ) , and in a closely related approach by d.hadjimichef , g.krein , s.szpigel and j.s.daveiga , phys .
b367 , 317 ( 1996 ) .
m.lscher , comm .
104 , 177 ( 1986 ) ; 105 , 153 ( 1986 ) ; nucl .
b354 , 531 ( 1991 ) ; s.r.sharpe , r.gupta and g.kilcup , nucl .
b383 , 309 ( 1992 ) ; r.gupta , a.patel and s.r.sharpe , phys .
d48 , 388 ( 1993 ) ; m.fukugita _ et al .
71 , 2387 ( 1993 ) ; phys .
d52 , 3003 ( 1995 ) ; s.aoki _ et al .
_ , hep - lat/9911025 .
( fig.5 shows the scattering length quoted by aoki _
et al._. )      t.barnes and e.s.swanson , phys .
d46 , 131 ( 1992 ) ; see also g .- q.zhao , x.g.jing and j-c.su , phys . rev .
d58 , 117503 - 1 ( 1998 ) ; d.hadjimichef , g.krein , s.szpigel and j.s.daveiga , phys .
b367 , 317 ( 1996 ) ; ann .
268 , 105 ( 1998 ) ; b.masud , j.paton , a.m.green and g.q.liu , nucl .
a528 , 477 ( 1991 ) ; b.masud , phys .
d50 , 6783 ( 1994 ) .
t.barnes , n.black , d.j.dean and e.s.swanson , phys .
c60 , 045202 ( 1999 ) .
c version has an incorrect figure ( a duplicate of fig.4 was used as fig.5 ) , and minus signs were deleted from the lower @xmath337-axis in both figures .
the original , archived version nucl - th/9902068 has the correct figures .
d.blaschke and g.rpke , phys .
b299 , 332 ( 1993 ) ; k.martins , d.blaschke and e.quack , phys .
c51 , 2723 ( 1995 ) ; k.martins , prog .
nucl . phys .
36 , 409 ( 1996 ) ; d.blaschke , g.burau and yu.l.kalinovsky , nucl - th/0006071 .                                            for @xmath338 see v.dorofeev ( ves collaboration ) , in proceedings of whs99 , hep - ex/9905002 ; g.s.adams _ et al .
( e852 collaboration ) , phys .
81 , 5760 ( 1998 ) . for @xmath4
see d.r.thompson _ et al . _
( e852 collaboration ) , phys .
79 , 1630 ( 1997 ) ; a.abele _ et al .
_ ( crystal barrel collaboration ) , phys .
b423 , 175 ( 1998 ) ; s.u.chung _ et al . _
( e852 collaboration ) , phys .
d60 , 092001 ( 1999 ) .
n.isgur , r.kokoski and j.paton , phys .
54 , 869 ( 1985 ) ; t.barnes , f.e.close and e.s.swanson , phys . rev .
d52 , 5242 ( 1995 ) ; c.bernard _ et al . _
( milc collaboration ) , phys .
d56 , 7039 ( 1997 ) ; p.lacock _ et al .
_ ( ukqcd collaboration ) , phys .
b401 , 308 ( 1997 ) .