quantum information theory is an area of scientific investigation that has witnessed an enormous progress in the last decade @xcite . in the framework of quantum information science , paradigmatic systems of condensed matter physics such as , e.g. ,
spin chains and harmonic lattices , are analyzed from the point of view of their information content and ability to process and transfer information .
fundamental concepts of statistical mechanics and probability , such as the shannon - von neumann entropy @xcite , play a central role in the quantification of bipartite quantum entanglement @xcite . on the other hand , exciting advances in quantum information research
are in a sense paying back the debt : the mathematical and theoretical toolkit developed for the characterization and quantification of quantum entanglement has proven useful to tackle questions and improve our understanding in the investigation of strongly correlated systems and quantum phase transitions @xcite .
perhaps the most interesting development stemming from the application of the tools of entanglement theory concerns the study of correlation scaling in diversely connected systems , with its relations to conformal field theory and the proof of longstanding conjectures on entropic area laws in spin and harmonic lattice systems @xcite .
these investigations and the associated results are of fundamental interest because they allow to enhance the control and establish firmer bounds and limits of applicability on many important tools of simulation and numerical analysis of complex many - body systems , a fundamental task in condensed matter , such as the density matrix renormalization group ( dmrg ) algorithms @xcite , matrix product representations @xcite , the entanglement renormalization ansatz @xcite , and weighted graph states methods @xcite .
the pioneering application of quantum information concepts to condensed matter was the observation that two - body entanglement , as quantified by the concurrence , in the ground state of a cooperative system , exhibits peculiar scaling features approaching a quantum critical point @xcite .
these seminal studies helped to clarify that at quantum phase transitions @xcite , the dramatic change in the ground state of a many - body system is associated to , or reflected by , a change in the way that quantum fluctuations are correlated , that is , in the way entanglement is distributed among the elementary constituents .
many recent efforts have thus been aimed at understanding the behavior of different measures of entanglement at quantum phase transitions , and assessing the corresponding enhancement of properties useful for applications in quantum technology and quantum engineering @xcite . however , since most entanglement measures can be often rewritten in terms of the conventional n - point correlation functions , the presence of a divergence in the former at quantum criticality can be directly traced back to a divergence in the latter @xcite , although in some particular instances entanglement - based studies allowed to discover novel types of phase transitions @xcite .    finally , a more technologically oriented product of the interaction between these two areas of research is the study of quantum spin chains as natural information carriers and distributors @xcite , able to realize tasks such as quantum state transfer @xcite and storage @xcite , possibly with maximal fidelity in chains with suitably engineered couplings and dynamics .
quantum engineering and quantum technology are being indeed developed at a fast pace .
quantum devices are vigorously pursued for applications ranging from nano - sciences to quantum computation and entanglement - enhanced metrology @xcite . despite a large variety of possible implementations involving different physical systems ,
many relevant properties of such devices can be in fact investigated in a unified setting by appropriate mappings to quantum spin models @xcite .
control of separability and entanglement in ground states of quantum spin models plays an important role in entanglement - enabled quantum technology applications , and needs to be characterized and quantified in different physical regimes @xcite . from a more fundamental perspective
, the determination of exact solutions endowed with precisely known properties of separability or entanglement , can be of great relevance in the study of advanced models of condensed matter and cooperative systems that are in general not exactly solvable .    in this paper
we exploit quantum information tools to tackle another important problem in condensed matter physics : the determination and characterization of exact ground states of cooperative systems characterized by the property of being in the form of a fully factorized , tensor product of single - particle pure states , with no quantum correlations between the individual constituents .
historically , the occurrence of totally factorized ( unentangled ) ground states of quantum many - body systems was first discovered in the one - dimensional anisotropic heisenberg model with nearest - neighbor interactions by kurmann , thomas and mller @xcite , by adopting a direct method in terms of product state ansatz .
this result was later re - derived and extended to two dimensions using quantum monte carlo numerical methods @xcite .
the direct product - state ansatz method of kurmann , thomas and mller has been later extended to anisotropic rings with interactions of arbitrary range , by exploiting the observation that factorized ground states break the parity symmetry @xcite .
complex quantum systems exhibiting cooperative behaviors , whose ground states are typically entangled @xcite , may thus admit , for some non trivial values of the hamiltonian parameters , a ground state which is completely separable .
the occurrence of factorization at finite or even strong values of the couplings is thus an effect of a delicate balancing between interactions and external fields .
the phenomenon of ground state factorization is particularly intriguing as it appears to be associated with the presence of an `` entanglement phase transition '' with no classical counterpart @xcite .
furthermore , for the purposes of quantum engineering applications , that employ distributed entanglement in order to manipulate and transfer information @xcite , factorization points need to be exactly identified and either avoided , in order to guarantee the reliable implementation of quantum devices , or properly exploited for the dynamical creation of strongly entangled multipartite states of large assemblies of microsystems ( graph and cluster states ) . finally , and also quite importantly , for models not admitting exact general solutions , achieving knowledge of the exact ground state , even if only for the restricted nontrivial set of parameters associated to factorization , would allow ( i ) to prove the existence of an ordered phase and characterize it ; ( ii ) to build variational or perturbative approximation schemes around the exact factorized solution , that may then be used as test benchmarks for the validity and the precision of numerical algorithms and simulations .
unfortunately , the direct method based on the ansatz of a factorized solution is neither sufficiently general nor mathematically tractable , apart from the notable exceptions of one - dimensional spin chains with short - range @xcite , infinite - range @xcite , and some types of long - range exchange interactions @xcite .    in a recent work
@xcite we introduced an analytic method that allows to determine exactly the existence of factorized ground states and to characterize their properties in quantum spin models defined on regular lattices of any size ( finite as well as at the thermodynamic limit ) , in any spatial dimension , and with spin - spin interactions of arbitrary range . in correspondence to rigorously established ground state factorizability , the method also allows to determine novel sets of exact solutions in generally non exactly solvable models .
all the previous particular findings can be rigorously re - derived and extended within a unified framework inspired by concepts of quantum information science .    in the present work
we extend the original method introduced in ref .
@xcite and we develop a general self - contained theory of ground state factorization in ( frustration - free ) quantum spin systems : with no _ a priori _ assumption on the nature of the magnetic ordering , we derive a complete and closed set of conditions that have to be satisfied by the hamiltonian interaction parameters in order for the ground state to be factorized , at precisely determined values of the external magnetic field , with well defined values of the magnetization and of the energy .
besides the completely rigorous re - derivation of the few previously known results on quantum factorization points @xcite , our general analytic method allows to determine novel , exact factorization points and factorized ground states in , generally non exactly solvable , spin-@xmath0 models as well as in the corresponding hamiltonian models of any higher spin .
these results hold true regardless of spatial dimensionality and interaction range , and can be easily extended even to non - translationally invariant systems .
we illustrate the versatility of the method with a series of different applications , including models with nearest - neighbor interactions ( for which we provide the full factorization diagram , generalizing the results of refs .
@xcite and of our recent work @xcite ) , models with long - range interactions ( notably including , among others , the fully connected infinite - range lipkin - meshkov - glick ( lmg ) model @xcite and systems with unbalanced next - to - nearest - neighbor interactions on cubic lattices ) , and models with spatial anisotropies .
the paper is organized as follows .
an excursus on the problem of ground state factorization is presented in sec .
[ sec.history ] . in sec .
[ sec.squoexe ] we recall the quantum information toolkit upon which our method is based , in particular the formalism of single - spin , or single - qubit , unitary operations and associated entanglement excitation energies , previously introduced for the characterization and quantification of bipartite entanglement from an abstract geometrical perspective @xcite . the general method for the rigorous determination of factorized ground states , that generalizes the scheme introduced in ref .
@xcite , is discussed in full detail in sec .
[ sec.method ] .
relevant applications and examples , including models defined on lattices of different spatial dimensions and with different interaction ranges as well as extensions to models with spatial anisotropies and with dimerized interactions are discussed in sec .
[ sec.short ] , sec .
[ sec.long ] , sec .
[ sec.lmg ] , and sec .
[ sec.anisotropy ] .
concluding remarks and outlooks on future perspectives are presented in sec .
[ sec.concl ] .
quantum fluctuations in the ground states of cooperative many - body systems are typically highly correlated , that is , those states contain in the most general instance a strong degree of entanglement distributed among the individual components .
this is one of the reasons why their exact characterization is often hard to be accomplished by analytical or even numerical methods .
however , there may exist particular solutions corresponding to special values of the hamiltonian parameters such that the ground state becomes `` classical - like '' , in the sense of being a full product state of single - site factors ( completely separable , fully factorized state ) .
it is quite surprising that despite the inherent appeal of such exact , special solutions , and the renewed interest stemming from possible applications to quantum technology , very few analytical advances were obtained in the last quarter of a century .      the first systematic work on ground state factorization in quantum spin models was completed by kurmann , thomas , and mller in 1982 @xcite .
they focused on the one - dimensional nearest - neighbor anisotropic heisenberg - like antiferromagnetic spin-@xmath0 model in transverse field , with the assumption that all the spin - spin couplings take nonnegative values .
et al . _ identified a precise value of the external magnetic field , the `` factorizing field '' , corresponding to which the system admits a fully factorized ground state .
they employed a very direct ,
brute force " method , i.e. formulating a product state ansatz and verifying if and for what values of the hamiltonian parameters it satisfies the stationary schrdinger problem with lowest energy . even though this approach appears obvious and extremely easy to pursue , the positive result obtained by kurmann _
_ in the case of the 1-d heisenberg model was enabled by some special properties of the case they analyzed .
in fact , in order to verify the eigenvalue equation for the product state ansatz , one needs to decompose the total hamiltonian of the system into a series of pairwise terms , which could be done easily for the model studied by kurmann _
_ since all the `` pairwise components '' of the factorizing magnetic field had the same weight in this particular case . however , this direct method of determining factorization is in fact significantly less obvious , useful , and comprehensive , when one tries to apply it to more general cases .
for instance , considering models with interaction terms of growing spatial range , and/or increasing the spatial dimensionality of the lattice , and/or treating extensions to spatially anisotropic or non - translationally invariant models , results in increasingly hard `` guesses '' concerning the partitioning of the external field into pairwise components and in increasingly nontrivial verification steps . depending on the structure of the model and on the form of the product state ansatz , the direct method may either fail to detect factorization points at all , or it may detect only a subset of all physically realizable factorized solutions , or finally it may lead to intractable sets of conditions . these limitations of the direct method become evident , as we will see in the following , even in the simplest 1-d @xmath1 model with nearest - neighbor interactions analyzed by kurmann and coworkers . in their original work , due to the insufficient generality of their proposed product ansatz
, they failed to identify a wide range of situations that allow the existence of factorized ground states . as we will see when
discussing models with competing interactions of different spatial ranges , a further serious fault of the direct method , that becomes incurable in the presence of frustration , is its incorrect assessment of factorized ground states that are in fact excited energy eigenstates of the system .    following the seminal analysis by kurmann _
_ , the direct method has been applied to the fully connected infinite - range lipkin - meshkov - glick model by dusuel and vidal @xcite , and to the anisotropic heisenberg chain with the same interactions of arbitrary range along the three spatial directions by hoeger et al .
@xcite and by rossignoli _
we would like to remark that these works , that include the results of refs .
@xcite as special cases , provide an interesting characterization of factorization points in finite ( one - dimensional ) translationally invariant lattice spin models with periodic boundary conditions as those points where the parity symmetry is broken .      a general and rigorous approach to the problem of ground state separability , completely different from the direct method based on the product state ansatz , has been introduced by us in ref .
@xcite and is fully developed in the present work . this analytic method allows _
inter alia _ to complete the factorization diagram of the anisotropic spin-@xmath0 models without restrictions on the sign of the couplings , and to go significantly beyond the limitations of the previously mentioned studies .
our approach is based on an additional characterization of product states ( including product ground states ) that employs concepts and tools of quantum information and entanglement theory .
the strong point of the quantum informatic approach to ground state factorization lies in its generality and its ensuing independence from the types and ranges of interactions present in the hamiltonian , the lattice size , and the spatial dimension , which makes it applicable in the most general cases and therefore especially useful in all the situations that are either intractable or incomplete if one resorts to the direct method or to numerical techniques .
in fact , in the original formulation of the quantum informatic method @xcite there was still trace of a residual `` ansatz '' in the initial stage , since the magnetic ordering of the candidate factorized ground state had to be in some sense guessed , and imposed _
a priori_. moreover , in the final stage , we were not yet able to derive a completely general set of conditions in order to determine in each situation whether a candidate factorized state that turned out to belong to the spectrum of the hamiltonian was actually a ground state or belonged to the set of excited states .
this had then to be decided case by case in each specific instance .    in the present work
we proceed further and develop the complete method to determine and characterize factorized ground states in quantum cooperative systems .
we provide the systematics that allows to determine the occurrence of full ground state separability given a generic quantum spin hamiltonian .
the magnetic ordering is derived as a consequence of the theory and not imposed _ a priori_. a complete set of general conditions is established for the verification of ground state factorization , and the elegance of the method does not hinder practical usefulness in the application to specific models and situations .    we would like to recall that numerical advances on the problem of ground state factorization have been realized by the monte carlo study of ground state factorization in two - dimensional anisotropic nearest - neighbor spin models @xcite .
unfortunately , beyond this important result the numerical approach appears to have a very restricted potential for extensions , since the required computational power increases rapidly with the spatial dimension and/or with the range of interactions .
we will show that the difficulties inherent to the direct analytic approach and to the numerical method are intrinsically absent in our formalism .
in the present section we provide a brief overview of some recent results on the geometrical interpretation of the entanglement measure known as `` tangle '' ( linear entropy ) @xcite and the consequences it entails for the relationship between energy and pairwise entanglement @xcite .
such results embody the premises for our all - analytic characterization of ground state separability in quantum spin systems .    here and in the following we consider generic systems made of @xmath2 spin-@xmath0 elementary constituents , or , in the language of quantum information , @xmath2 qubits .
for such a system , we introduce the set of single - qubit unitary operations ( _ squos _ ) @xmath3 @xcite , defined as the unitary transformations that act separately as the identity on any spin of the system except the @xmath4-th one , on which they act instead as hermitian , unitary , and traceless operators : @xmath5 in eq .
( [ e.uk ] ) @xmath6 stands for the identity operator on the @xmath7-th spin of the system while the generic hermitian , unitary , and traceless operator @xmath8 can be expressed as a linear combination of the standard spin operators ( @xmath9 ) defined on the @xmath4-th spin : @xmath10 in the previous definition the parameters @xmath11 and @xmath12 take values respectively in the ranges @xmath13 $ ] and @xmath14 .
physically , @xmath15 corresponds to a rotation in the spin space , and the traceless condition imposes the rotation operation to be orthogonal to the identity , so basically a combination of spin flip and phase flip .
a squo transforms every pure state @xmath16 into a new state @xmath17 , that in general differs from @xmath18 .
we may quantify the action of the squo in terms of the trace distance between the original and the transformed state , defined as @xmath19 .
the distance @xmath20 varies in the interval @xmath21 $ ] and vanishes if and only if the two states coincide , meaning that the considered squo leaves the original state unchanged . for a given initial state and for any site @xmath4
, one can determine the _ extremal - squo _
( e - squo ) that minimizes such a distance .
one can prove that the e - squo is uniquely determined by the following conditions on the angular variables @xcite : @xmath22 where @xmath23 denotes the expectation value of the spin operator @xmath24 on the unperturbed state @xmath18 ( @xmath25 ) .
the e - squo plays a crucial role in our analysis : for any pure state defined on a general system of interacting qubits one can prove @xcite that the square of the distance associated to the action of the e - squo coincides with the linear entropy @xmath26 $ ] of the @xmath4-th qubit ( where @xmath27 denotes the reduced density matrix of qubit @xmath4 ) .
this quantity , also known as the _ tangle _ @xmath28 in the literature , is an entanglement monotone for qubit systems and satisfies important monogamy constraints @xcite .
being a monotonic function of the von neumann entropy of entanglement @xcite , it quantifies the bipartite block entanglement present in the state @xmath18 between spin @xmath4 and the remaining @xmath29 constituents , and is thus a proper measure of single - site entanglement @xcite .
one has the following relation @xcite :    @xmath30 \ ; .\ ] ]    this geometric - operational characterization of the tangle allows to develop a quantum informatic approach to the problem of characterizing the correlation properties of the ground state in many - body quantum systems @xcite .
let us consider a collection of spin-@xmath0 systems defined on a regular lattice , with hamiltonian @xmath31 , and let @xmath18 denote from now on the ground state of @xmath31 . if the system is in the ground state , the application of a squo on a generic spin of the lattice is necessarily associated to an energy transfer that one must provide to the system in order to realize the selected squo . in other words ,
a squo perturbs the system , and this perturbation results in an increase of the average energy ( obviously , in general the transformed state needs not to be an eigenstate of the hamiltonian ) .
such energy deviation from the ground - state value can be quantitatively defined as @xmath32 because @xmath18 is the ground state of @xmath31 , @xmath33 is a non negative defined quantity . considering the special relation between the e - squo and the tangle
we may expect that the associated amount of energy difference @xmath34 , appropriately named ` entanglement excitation energy ' ( exe ) , will have as well a strong , direct connection with the tangle and in general with the single - site entanglement in the ground state .
in fact , starting from the above definition of eq .
( [ e.ees ] ) and taking into account the property of the e - squo of leaving unchanged a fully disentangled state , it has been established in full generality @xcite that if the system is invariant under spatial translation and the hamiltonian of the system does not commute with any possible squo , then the vanishing of the exe is a _ necessary and sufficient _ condition to admit a fully factorized ground state @xmath18 @xcite .
hence , if the system that we want to analyze satisfies all the former hypotheses , the vanishing of the exe provides an additional condition , besides the vanishing of the tangle @xmath28 , that one needs to include in establishing a general analytic approach to the problem of ground state factorization .
the general method that was first partially developed in ref .
@xcite for the determination of existence , location , and exact form of factorized ground states rests on the following two main observations :    the ground state @xmath18 of a spin-@xmath0 hamiltonian @xmath31 is factorized if and only if the single - site entanglement ( tangle ) @xmath28 vanishes for all spins @xmath4 in the lattice ;    provided that @xmath35\neq 0 $ ] for every possible squo @xmath36 of the form eq .
( [ e.uk ] ) , the ground state @xmath18 of a spin-@xmath0 hamiltonian @xmath31 is factorized if and only if the exe @xmath37 vanishes @xmath38 .    in ref .
@xcite , the strategy to the determination and understanding of factorization for a given hamiltonian bore on the above two facts and on the requirement to assume as working point a phase endowed with some kind of magnetic order . by imposing the vanishing of the exe and of the tangle
, we showed how to determine uniquely the form of the completely disentangled state compatible with the assumed magnetic order , candidate to be the ground state of the system at a particular value of the external magnetic field ( labeled as `` factorizing field '' ) .
the final steps of the procedure concern the determination of exact analytic conditions for the candidate factorized state to be an eigenstate of the hamiltonian with lowest energy . in ref .
@xcite , we derived a closed set of eigenstate conditions , but establishing whether a factorized eigenstate is actually a ground state had to be left to a case by case analysis .    in the present paper
we generalize the method , extend its scope , and apply it to different classes of quantum spin models .
firstly , we will show that there is no need for _ a priori _ assumptions on the magnetic order .
in fact , we will prove that the existence of a factorized ground state is a necessary and sufficient condition for the existence of quantum phase transitions and ordered phases in translationally - invariant quantum spin models with exchange interactions and in the presence of external fields .
moreover , the actual structure of the factorized ground state determines automatically the kind of ordered phase .
further , we derive a complete set of rigorous conditions for the candidate disentangled state to be an eigenstate and a ground state . in particular , for systems with no frustration
we prove that if the factorized state is an eigenstate , then it is always a ground state of the system .
we illustrate the method in detail by considering its application to the general , translationally invariant , exchange hamiltonian @xmath31 for spin-@xmath0 systems on a @xmath20-dimensional regular lattices , with spin - spin interactions of arbitrary range and arbitrary anisotropic couplings : @xmath39 here @xmath40 ( and similarly @xmath41 ) is a @xmath20-dimensional index vector identifying a site in the lattice , @xmath42 @xmath43 stands for the spin-@xmath0 operator on site @xmath40 , @xmath44 is external field directed along the @xmath45 direction , @xmath46 is the distance between two lattice sites , and @xmath47 is the spin - spin coupling along the @xmath48 direction .
translational invariance is ensured by the fact that all the couplings depend only on the distance @xmath49 between the spins .
this type of hamiltonian encompasses a large variety of models and spans over several universality classes including , among others , the ising , xy , heisenberg , and xyz symmetry classes .
besides their importance in quantum statistical mechanics and in the theory of quantum critical phenomena , these models play an important role in the study of various schemes of quantum information and quantum communication tasks with quantum many - body systems @xcite .
it is a rather straightforward exercise to verify that the hamiltonian eq .
( [ e.hamiltonian ] ) never commutes with the squos @xmath50 .
translational invariance and non - commutativity of the hamiltonian with the squos guarantee that the vanishing of the exe is a necessary and sufficient condition for the full separability of the ground state [ statement ( ii ) above ] .
since the external field @xmath44 is uniform on the lattice , it forces the existence of a site - independent , nonvanishing magnetization @xmath51 along the @xmath45 axis : @xmath52 . on the other hand , at a factorization point , statement ( i ) above imposes the vanishing of the single - site entanglement on all sites of the lattice .
therefore , if a factorization point exists , the magnetizations along the directions orthogonal to the external field must assume the following form : @xmath53 where @xmath54 is the modulus of the projection of the magnetization on the @xmath55 plane and the local orientations @xmath56 , and hence the type of magnetic ordering in the @xmath55 plane , remain undetermined .
since @xmath51 is site - independent , then @xmath57 is site - independent as well .
taking into account eq .
( [ e.minimumpoint ] ) of sec .
[ sec.squoexe ] one sees that , at least at the factorization point , @xmath58 , i.e. one of the two angles that fix the orientation of the e - squo , does not depend on the site index
. therefore @xmath59 @xmath60 .    collecting the former results and considering that in the presence of a factorized ground state , if it exists , all the correlation functions are products of single - site expectations values
, we obtain that the condition on the vanishing of the exe associated to the generic site @xmath61 reads :    @xmath62    in the former equation @xmath63 is the _ factorizing field _ , i.e. the value , if it exists , of the external field for which ground state factorization occurs .
since we are looking for a fully factorized ground state , eq .
( [ e.exe ] ) must be satisfied simultaneously for all @xmath61 and , henceforth , eq .  ( [ e.exe ] ) must be site - independent .
this fact in turn implies that the following relation must be satisfied by every angle @xmath56 : @xmath64 eq .
( [ e.exe ] ) must hold independently of the site on which the e - squo acts and the value of the constant @xmath65 must be identified by determining the expression of the energy density , i.e. the energy per site , associated to the factorized ground state .
given the total ground - state energy @xmath66 and the total number of sites @xmath2 , the latter quantity can be written as @xmath67 where @xmath68 is the coordination number , that is the number of sites at distance @xmath49 from a given site . in order to determine the magnetic order in the candidate factorized ground state
, one can now exploit the fact that , by definition , it must minimize the energy of the system .
therefore such state will be simply characterized by the magnetic order that minimizes the per - site contribution @xmath65 , depending on the hamiltonian parameters and the geometric structure of the lattice . thus the general conditions for full separability of the ground state lead to the uncovering of the magnetic order in the @xmath55 plane . in the original formulation of the quantum - informatic method @xcite one had to assume from the beginning the type of ordered phase , and this resulted in constraints on the magnetization and limitations in the range of the coupling coefficients .
the outcome of the minimization of the energy density depends on the presence or the absence of frustration in the system , either due to the lattice geometry , or to the structure of the hamiltonian , or to both causes .
for hamiltonians belonging to the class of eq .
( [ e.hamiltonian ] ) , frustration can be due both to the geometry of the lattice , e.g. for systems with nearest - neighbor antiferromagnetic interactions defined on a lattice that admits closed loops with an odd number of spins , and to the competition of interactions of different spatial orders , e.g. in systems defined on linear chains with nearest - neighbor ferromagnetic and next - nearest - neighbor antiferromagnetic interactions .
frustration effects in systems with interactions of different spatial ranges arise whenever it happens that the minimization of the energy associated to interactions over a certain spatial scale precludes the possibility of minimizing the energy associated to interactions over a different spatial scale . in the absence of frustration , i.e. when the lattice geometry and the interactions between the spins are such to allow all energy minimizations simultaneously , it is straightforward to prove that every possible factorized ground state must be characterized by one of only four magnetic orders out of the many possible ones .
the complete proof of this statement is reported in appendix [ appendixfourcases ] .
this restricted set includes ferromagnetic ordering along the @xmath69 direction ( @xmath70 ) or along the @xmath71 direction ( @xmath72 ) , and antiferromagnetic ordering along the @xmath69 direction ( @xmath73 ) or along the @xmath71 direction ( @xmath74 ) , the type of order effectively present in the ground state can be determined by comparing the values of the different `` net interactions '' along the @xmath69 and @xmath71 directions , where the net interactions are defined as follows : @xmath75 with @xmath76 .
the net interaction along the @xmath45 direction is defined equivalently as @xmath77 as a function of the net interactions , the magnetic orders in the system are determined by the values of @xmath78 : @xmath79 it is interesting to observe that , in terms of the net interactions , one can immediately establish the presence or absence of frustration in the system : any frustration arising from the presence of competing interactions would imply that @xmath80 , with @xmath76 depending on which of the two axes is characterized by a non vanishing value of the magnetization .
we remark that here and in the following we are not considering the particular situation in which a _ saturation _ occurs rather than a true ground state factorization : this instance will be discussed separately in subsec .
[ s.sec.saturation ] .    from now on , without loss of generality , we will specialize our analysis to the the case of an antiferromagnetic ordering along the @xmath69 direction . with trivial modifications , all the steps and results that will be obtained in the following hold as well in the remaining three cases of eq .
( [ fourcases ] ) . at the end of the procedure we will present a summary table collecting the main results for all the four types of magnetic orders compatible with factorization .
given an antiferromagnetic ordering along the @xmath69 direction , and keeping in mind that in such case @xmath73 , eq .
( [ e.magnxy ] ) implies that @xmath81 while @xmath82 with @xmath83 .
the @xmath84 sign reflects the fact that the antiferromagnetic order discriminates between two sublattices , each characterized by the opposite sign of the staggered magnetization along @xmath69 . consequently , the general condition eq .
( [ e.exe ] ) takes the form @xmath85 \ , , \end{aligned}\ ] ] that admits two solutions for @xmath86 , @xmath87 however , since the e - squo must be unique @xcite , the two solutions must coincide .
factorization thus requires that @xmath88 moreover , the requirement of a vanishing tangle ( single - site entanglement ) imposes @xmath89 combining eqs.([e.mz][e.mx ] ) via eq .
( [ e.minimumpoint ] ) yields a closed expression for the phase @xmath86 as a function of the hamiltonian parameters and of the factorizing field @xmath63 : @xmath90 therefore , eq .
( [ e.theta ] ) , together with the knowledge of @xmath91 , once inserted into eq .
( [ e.ok ] ) , determines , independently of the value of the magnetizations , the form of the e - squo at factorization : @xmath92 the form of the candidate factorized state @xmath93 is then readily determined by imposing that the action of the e - squo on all lattice sites leaves the ground state invariant : @xmath94    > from eqs.([e.inalterate][e.oz ] ) we find then that the candidate factorized ground state @xmath93 has to be a tensor product of single - site states that are the exact eigenstates of the operators @xmath95 with eigenvalue @xmath0 : @xmath96    equation ( [ e.candidate ] ) expresses the general form that a factorized ground state must assume given an antiferromagnetic order along @xmath69 .
this result is independent of the actual values assumed by the staggered magnetizations @xmath97 s .
we remark that in order to establish the general form of the candidate ground state it is essential to impose the vanishing of the exe .
in fact , without requiring it , the phase @xmath86 and the candidate factorized ground state become explicitly dependent on the magnetizations . then , imposing the conditions for eq .
( [ e.candidate ] ) to be an eigenstate of the hamiltonian , one merely obtains a relation between the factorizing field @xmath63 and the magnetizations @xmath98 s . if the analytic expression of at least one of the @xmath98 s is known , it is still possible to obtain an expression of the field @xmath63 in terms of the hamiltonian parameters .
however , this actually occurs in few special , exactly solvable , cases ( for instance , the xy model @xcite ) .
otherwise , in general one needs to resort to numerical evaluations of the magnetizations or introduce other approximations in the analysis .
this problem is completely eliminated in the general analytic framework based on the invariance of the factorized ground states under the action of e - squos and on the vanishing of the exe at factorization .    before proceeding
further let us clarify a relevant point .
because of the even - odd symmetry of the energy spectra of the system , as elucidated also in refs .
@xcite , the assignment of one of the two different staggered magnetizations to each sublattice is arbitrary . according to this observation
, one sees that , inverting the sign of the magnetization along the @xmath69 direction on each site , all the conditions for factorization are still satisfied but the candidate ground state takes the form @xmath99 obviously the two states in eq .
( [ e.candidate ] ) and in eq .
( [ e.candidate-1 ] ) are distinguishable and both @xmath93 and @xmath100 are legitimate candidate ground states of the system , as well as general linear combinations of the two ( the latter will be in general highly entangled ) . here and in the following we
assume to be working in a situation of broken symmetry , i.e. after that a small perturbing external field along the @xmath69 direction has lifted the degeneracy between @xmath93 and @xmath100 , e.g. lowering the energy of the former and rising that of the latter .
after a time long enough to ensure convergence to equilibrium , i.e. relaxation to the state @xmath93 , the perturbation is switched off : in the thermodynamic limit , this ensures that the system will remain indefinitely in the state @xmath93 .
having determined the exact form of the candidate factorized ground state , the subsequent step concerns the determination of the conditions for its occurrence , i.e. the conditions under which a state of the form eq .
( [ e.candidate ] ) is indeed the eigenstate of @xmath31 , eq .
( [ e.hamiltonian ] ) , with the lowest energy eigenvalue . to this aim ,
taking the value of the external field at the factorization point , @xmath101 , it is useful to decompose the total hamiltonian as a sum of pairwise hamiltonian terms @xmath102 , so that @xmath103 , where each term @xmath102 involves only the degrees of freedom of a single pair of spins : @xmath104 here the quantity @xmath105 plays the role of the `` component '' of the factorizing field that acts on the selected pair of spins and obeys the following relation : @xmath106 the above ensures that @xmath107 satisfies eq .  ( [ e.theta ] ) .
proving that the candidate factorized ground state @xmath93 is an eigenstate of the total hamiltonian @xmath31 is equivalent to prove that it is a simultaneous eigenstate of all pair hamiltonians @xmath108 , or , more precisely , that the projection of @xmath93 onto the subspace of two given spins @xmath40 and @xmath109  which is still a pure state since @xmath93 is a tensor product of single - spin pure states  is an eigenstate of @xmath108 for every pair @xmath110 .
to proceed , we associate to each lattice site a set of orthogonal spin operators defined as follows : @xmath111 it is immediate to observe that @xmath112 . by inverting eqs .
( [ e.axayaz ] ) , we can conveniently re - express the standard spin operators as functions of the new set of operators @xmath113 : @xmath114 inserting eq .
( [ e.sxsysz ] ) in eq .
( [ e.couple ] ) we obtain the expression of the pair hamiltonian @xmath115 as a function of the sets of operators @xmath113 and @xmath116 :    @xmath117 -2 h_f^r \sin \theta \left [ { a_{\underline{i}}^x \left ( a_{\underline{j}}^z -\frac{1}{2}\right ) + \left ( a_{\underline{i}}^z-\frac{1}{2}\right )
a_{\underline{j}}^x } \right ] \nonumber \\ & & + \ a_{\underline{i}}^x a_{\underline{j}}^x \left ( \sin^2 \theta j_z^r+\cos^2\theta ( -1)^r j_x^r \right)+ a_{\underline{i}}^y a_{\underline{j}}^y j_y^r ( -1)^r - h_f^r \cos \theta \left ( { a_{\underline{i}}^z + a_{\underline{j}}^z } \right ) \ , , \end{aligned}\ ] ]    where we recall that @xmath49 is the distance between sites @xmath40 and @xmath109 and that , in the presence of an antiferromagnetic order along the @xmath69 direction , @xmath118 .
the spin - pair projection @xmath119 of @xmath93 is an eigenstate of the pair hamiltonian @xmath120 either if all terms in eq .
( [ e.pairhij ] ) admit @xmath119 as an eigenstate or if they annihilate it . by construction , recalling that @xmath121 and that @xmath122 is the eigenstate of @xmath123 with eigenvalue equal to @xmath0 , the spin - pair projection of @xmath93 is an eigenstate of the pair hamiltonian if the following condition holds : @xmath124 for general models of the type eq .
( [ e.hamiltonian ] ) with a given , arbitrary , maximum spatial range of interaction @xmath125 , eq .
( [ e.condition ] ) admits the following solution : @xmath126 for all distances @xmath127 , i.e. for all distances associated with nonvanishing couplings ( we recall that @xmath128 and @xmath129 ) . the set of equations ( [ e.solcondition ] ) determines the conditions that must be satisfied simultaneously by the phase @xmath86 , which is unique , in order for @xmath93 to be an eigenstate of the total hamiltonian eq .
( [ e.hamiltonian ] ) .
( [ e.solcondition ] ) discriminate quite clearly between short - range models , i.e. models with only one nonvanishing coupling at @xmath130 ( models with nearest - neighbor interactions ) and all other possible models containing finite- or infinite - range interaction terms .
specifically , models with only nearest - neighbor interactions are characterized by the fact that eqs .
( [ e.solcondition ] ) reduce to a single condition and hence if the value of the r.h.s . of eq .
( [ e.solcondition ] ) falls in the interval @xmath21 $ ] ( associated to the permitted values of @xmath131 in terms of the interaction parameters ) , one can immediately conclude that the models under investigation admit a factorized eigenstate . on the other hand , for models containing interaction terms of longer range  so that eqs .
( [ e.solcondition ] ) include two or more conditions  factorized energy eigenstates are allowed if and only if the r.h.s . of all conditions in the set ( [ e.solcondition ] )
take values in the interval @xmath21 $ ] and , moreover , they all coincide .
summing term by term over the index @xmath49 all the relations in eq .
( [ e.condition ] ) , taking into account eq .
( [ e.theta ] ) , and solving for @xmath63 , one eventually obtains the exact expression of the factorizing field as a function of the net interactions : @xmath132 we remark that eq .
( [ e.factorizingfield ] ) , which was derived in ref .
@xcite using some unnecessary auxiliary assumptions , is completely general and holds for lattices of arbitrary spatial dimension and for spin - spin interactions of arbitrary range .
accordingly , the angle @xmath86 that determines the direction of the e - squo can be also expressed as a function of the net interactions as follows @xmath133 in the case of a maximum range of interaction @xmath134 , relation eq .
( [ e.theta2 ] ) for the net interactions is a necessary but not sufficient condition for factorization , and further use of eqs .
( [ e.condition ] ) and eqs .
( [ e.solcondition ] ) is needed , as we will show with explicit examples in sec .
[ sec.long ] and sec .
[ sec.anisotropy ] .
we are finally left with the problem of establishing conditions for a factorized energy eigenstate @xmath93 to be indeed a ground state of the system .
a very simple sufficient but not necessary condition is that the spin - pair projection of @xmath93 onto the subspace of a pair of spins be the simultaneous ground state of every pair hamiltonian @xmath108 .
in fact , given the relation existing between all the pair hamiltonians @xmath108 in eq .
( [ e.couple ] ) and the total hamiltonian @xmath31 eq .  ( [ e.hamiltonian ] ) , it follows that the simultaneous ground state of all pair hamiltonians @xmath108 is also ground state of @xmath31 . to proceed ,
let us re - express the pair hamiltonians eq .
( [ e.pairhij ] ) in matrix form in the basis spanned by the eigenstates of the operators @xmath135 .
one has @xmath136 where @xmath137 and @xmath138 . making use of eq .
( [ e.theta2 ] ) and eq .
( [ e.hfij ] ) the eigenvalues of the matrix eq .
( [ e.matrix ] ) can be put in the form @xmath139 \nonumber \\
\varepsilon_1 & = & \frac{1}{8 } \left [ 3 ( -1)^r j_x^r   - j_z^r   -    \left ( { j_z^r -(-1)^r j_x^r } \right)\cos \left ( { 2\theta } \right )    \right ] \nonumber \\
\varepsilon_2 & = & \frac{1}{8 } \left [ - 3\left ( { ( -1)^rj_x^r   + j_z^r } \right )     + \left ( { j_z^r-(-1)^r j_x^r   } \right)\cos \left ( { 2\theta } \right )    \right]\nonumber \\
\varepsilon_3 & = & \frac{1}{8 } \left [ { - 3(-1)^rj_x^r   + 5j_z^r   +    \left ( { j_z^r- ( -1)^r j_x^r } \right)\cos \left ( { 2\theta } \right ) }   \right ] \nonumber \\   & & \end{aligned}\ ] ]    the eigenvalue @xmath140 is the one associated to the projection @xmath119 of the factorized eigenstate @xmath93 onto the four - dimensional hilbert space associated to the pair of spins @xmath141 .
therefore , the spin - pair projection is the ground state of the pair hamiltonian if @xmath142 . imposing this condition yields the following inequalities : @xmath143 & \ge & 0 \ ; ; \label{e.suffcondition.1 } \\ ( -1)^r j_x^r - j_z^r & \le & 0 \ ; . \label{e.suffcondition.2}\end{aligned}\ ] ] notably , these inequalities are automatically satisfied for any interacting spin system that is not frustrated and that verifies all of the eqs .
( [ e.solcondition ] ) .
namely , for every even value of the distance @xmath49 in the interaction range ( @xmath127 ) eq .
( [ e.suffcondition.2 ] ) implies that the denominator in the corresponding eq .
( [ e.solcondition ] ) is non - positive .
therefore , an acceptable value of @xmath131 can be obtained only if also the numerator is non - positive and exceeds or equals the denominator . these conditions in turn imply the order relation , @xmath144 . on the other hand , comparing eq .
( [ e.suffcondition.2 ] ) and eq .
( [ e.suffcondition.1 ] ) , one has that @xmath145 . therefore , comparing the two order relations yields finally @xmath146 ( @xmath49 even ) .
this condition is always verified in systems with an antiferromagnetic order along the @xmath69 direction and in absence of frustration and therefore the projection of the factorized energy eigenstate on the hilbert space of every pair of spins is indeed the ground state of every pair hamiltonian .
on the other hand , in the case of odd @xmath49 , eq .
( [ e.suffcondition.2 ] ) implies that the denominator of eq .
( [ e.solcondition ] ) is non - negative and , similarly to the former case , the numerator must be non - negative and not exceeding the denominator .
hence we obtain @xmath147 while from eq .
( [ e.suffcondition.1 ] ) we recover @xmath148 that , in turn , implies @xmath149 .
again , such inequalities are always verified in the presence of an antiferromagnetic order and in the absence of frustration .    [ cols="^,^,^,^ " , ]     collecting all these results , we have proved that the following holds :    * theorem*. _ for any cooperative system of spin-@xmath0 particles described by hamiltonians @xmath31 of the type [ eq .
( [ e.hamiltonian ] ) ] , characterized by an anti - ferromagnetic order along the @xmath69 axis as emerging from eq .
( [ fourcases ] ) and in the absence of frustration , the simultaneous verification of all eqs .
( [ e.condition ] ) is necessary and sufficient for the fully factorized state eq .
( [ e.candidate ] ) to be the exact ground state of @xmath31 when the external magnetic field takes the value @xmath101 determined by eq .
( [ e.factorizingfield])_.    this central result holds as well , with obvious modifications , when one considers the other different ordered phases that form the set identified in eq .
( [ fourcases ] ) . in tab .
[ t.summary ] we summarize and compare results for each of the four different possibilities .      before moving to apply the general method to specific models and examples of conceptual and physical relevance
, we discuss briefly the meaning of factorization and its marked differences with the phenomenon of  saturation " ( see below ) . as we have seen , and as it will appear even clearer in the discussion of the examples in the following sections , the physical mechanism of ground state factorization in quantum spin systems is due to a delicate kind of  balancing " between the coupling strengths that regulate the intensity ( and range ) of the interactions and the aligning effect of the external field , that tends to orient all spins along a given direction and to destroy all quantum correlations .
the remarkable aspect of factorization is then that it occurs , according to well defined conditions and constraints , at precisely defined , finite values of the couplings and of the fields in an ordered ( or symmetry - broken ) phase .
this means that factorization occurs when the system is relatively strongly interacting and the external aligning field is relatively weak .
moreover , we will see that the factorization point is always a precursor from below ( from the ordered phase ) of the quantum critical point , as first observed for the @xmath1 model with nearest - neighbor interactions @xcite .
saturation is the phenomenon of ground state factorization that occurs trivially when the value of the external field grows unboundedly compared to all other hamiltonian parameters .
therefore its aligning effect on the spins prevails against all other effects and quantum correlations are suppressed .
the main physical difference between proper factorization and saturation is then clear and is reflected at the level of the hilbert space of states .
namely , suppose that the system under study admits a fully factorized state that is the true ground state of the system at @xmath101 . then , when one moves away from the factorizing field , i.e. when @xmath150 , the factorized state ceases to be an eigenstate of the hamiltonian . on the contrary , in the presence of saturation ,
the trivially ensuing factorized state is always an eigenstate of the hamiltonian and remains the ground state of the system when perturbing the value of the external field @xmath44 .
let us formalize the above discussion .
consider a system that admits a fully factorized ground state . in this case the pair hamiltonian in eq .
( [ e.pairhij ] ) admits as eigenstate a factorized spin - pair state .
however , eq .  ( [ e.pairhij ] ) is the expression of the pair hamiltonian only at @xmath101 . if we consider a value of the external field different from the factorizing field we must include in eq .
( [ e.pairhij ] ) a term proportional to @xmath151 . obviously , for nonvanishing values of @xmath152 the presence of such term prevents the factorized state from being an eigenstate of the hamiltonian of the system .
on the contrary , if @xmath153 and hence if @xmath154 the extra hamiltonian term vanishes and the factorized state @xmath155 remains an energy eigenstate of the system for all values of the external field .
this result can be understood from a different point of view .
when @xmath156 , from eqs.([e.solcondition ] ) we immediately obtain that @xmath157 , which implies that the total magnetization along the @xmath45 axis is preserved , @xmath158=0 $ ] . therefore the hamiltonian of the system and the total magnetization operator admit a complete set of common eigenstates , discriminated by the eigenvalue of the total magnetization : the factorized state is the one ( and only one ) characterized by the maximum eigenvalue .
we now move on to show that the factorized eigenstate , in the case of saturation , is the system s ground state in an infinite range of values of the external field . here
and in the following we recall that since the external field is directed along the conventional @xmath45 direction , only saturation with spin alignment along the @xmath45-axis is allowed . at very low values of @xmath44
the ground state of the system will possess a small total magnetization along the field direction .
however , as @xmath44 increases , also the total ground - state magnetization along the @xmath45 direction will increase until all spins align with the field and saturation occurs . by increasing @xmath44
further , it is not possible for the system to evolve in a state with higher values of the magnetization along the field direction .
hence , the factorized state remains the ground state of the system for all values of the external field @xmath159 , where by @xmath160 we denote the value of the field for which the system reaches saturation ( `` saturating field '' ) .
let us consider the matrix representation , similarly to eq .
( [ e.matrix ] ) , of the pair hamiltonian at arbitrary values of @xmath44 . taking into account the condition @xmath153 we have @xmath161 it is immediate to verify that if @xmath162 , the projection of the factorized state onto the considered pair of spins is associated to the lowest energy .
re - summing over all values of the spatial range index @xmath49 and taking into account the coordination number @xmath68 of the lattice , we have that the factorized state is the ground state of the total system if @xmath163 eq .
( [ e.hfsaturation ] ) is the rigorous condition for saturation .
we notice that it coincides with eq .
( [ e.factorizingfield ] ) in the limit @xmath164 , meaning that in this simple instance no ground - state factorization occurs other than saturation .    in the following sections we apply the general theory of factorization to different quantum spin models .
we will derive a series of exact results concerning the occurrence of factorized ground states in systems with various types of two - body interactions of different spatial ranges defined on regular lattices of arbitrary spatial dimensionality .
the general method enables to determine many novel exact factorized ground state solutions as well as the existence and the nature of ordered phases in models that are , in general , non exactly solvable .
moreover , it allows to recover quite straightforwardly the existing analytical and numerical results . in the process
, we will come to appreciate that ground state factorization is a phenomenon more common than previously believed . determining the existence of exact factorized solutions in non exactly solvable models should in turn allow to envisage controlled approximation schemes , e.g. perturbative or variational , to investigate the physics of non exactly solvable quantum spin models in the vicinity of quantum factorization points .
in this and in the following sections we will consider quantum spin models of the @xmath1 type eq .
( [ e.hamiltonian ] ) defined on regular lattices of arbitrary dimension .
let us first consider models with short - range interactions including only nearest - neighbor couplings : @xmath165 as in sec .
[ sec.method ] we restrict ourselves to the case in which such models are defined on regular lattices of a generic dimension , whose geometry does not induce any frustration , i.e. in which any closed loop is characterized by an even number of spins .
the absence of geometrical frustration together with the fact that all the couplings beyond the nearest neighbor spins vanish , ensure that it is always possible for the system to satisfy all the conditions on the net magnetic interactions .
therefore , as shown in sec .
[ sec.method ] , if the r.h.s . of eq .
( [ e.solcondition ] )  or , better , if the r.h.s . of the equation relative to each possible magnetic order ( see tab .
[ t.summary ] ) takes an acceptable value ( @xmath166 )  then the system admits a fully factorized ground state . in the case of antiferromagnetic
ordering along the @xmath69 direction , which is obtained if @xmath167 , from eq .
( [ e.solcondition ] ) it follows that a factorized ground state exists if @xmath168 .
analogous conditions are obtained for the other types of magnetic orders compatible with full ground state separability .
the results are summarized in fig .
[ f.factor ] , where @xmath169 is set assume positive values for reference . a very similar diagram of factorization holds in the case @xmath170 .     and
the coupling @xmath169 is set to be positive for reference .
analogous results hold in the case @xmath170 .
blank regions : no factorization allowed .
rightmost ( yellow online ) shaded region : presence of a fully factorized ground state supporting an antiferromagnetic order along the @xmath69 direction .
leftmost ( light blue online ) shaded region : fully factorized ground state with a ferromagnetic order along @xmath69 .
bottommost ( light red online ) shaded region : fully factorized ground state with a ferromagnetic order along the @xmath71 direction .
topmost ( light green online ) shaded region : fully factorized ground state with an antiferromagnetic order along @xmath71 .
the solid line @xmath172 corresponds to models that exhibit saturation rather than true factorization ; the dashed line @xmath173 accommodates for models with unresolved superpositions of two fully factorized ground states each one supporting a different magnetic order .
either one of the two possible ground states is then `` chosen '' by the system depending on the way the symmetry is broken , see text for further details .
all plotted quantities are dimensionless.,width=321 ]    as one can see from fig .
[ f.factor ] , the type of magnetic order corresponding to the factorized ground state determines four different regions in the space of the hamiltonian parameters ( one for each different ordering ) associated to the existence of a factorized ground state .
the four different domains of the interaction parameters are separated by two bisectrices that individuate special cases .
the line @xmath172 is constituted by the set of points in which saturation occurs in place of genuine factorization ( see subsec .  [ s.sec.saturation ] ) . on the other hand ,
the line @xmath174 , corresponds to models that , for a given value of the factorizing field @xmath63 and of the e - squo orientation @xmath86 , admit two degenerate , but distinct factorized states .
the choice of the ground state is then determined by the usual mechanisms of symmetry breaking : if an infinitesimal magnetic field is applied in the @xmath55 plane , the axis of magnetic alignment in the ground state is determined by the order in which the two limits @xmath175 and @xmath176 are realized .
it is important to observe that all the shaded regions of fig .
[ f.factor ] correspond to exact analytical ground - state solutions of the investigated short - range models , which are , in general , non - exactly solvable .
the exact , fully factorized ground states are tensor products of the form reported in the fourth column of tab .
[ t.summary ] , in which the value of the orientation parameter @xmath86 of the e - squo is retrieved by solving the equation in the second column of tab .
[ t.summary ] . obviously , such an exact solution holds only at the factorizing field @xmath101 . in the case of antiferromagnetic
order either along the @xmath69 or the @xmath71 direction , the latter is given by ( see tab . [ t.summary ] ) @xmath177 while in the case of a ferromagnetic phase , again unrespectfully if along @xmath69 or along @xmath71 , it reads @xmath178 given that the ground state is known exactly at the factorization point , all the relevant physical observables can be evaluated straightforwardly .
for instance , the ground state energy density at factorization reads @xmath179 in the presence of an antiferromagnetic ordering , while in the ferromagnetic case it is @xmath180    both the factorizing field and the ground state energy , unlike the e - squo orientation @xmath86 and the expression of the factorized ground states , are functions of the number of nearest neighbors per site ( the coordination number @xmath181 , where the index @xmath182 denotes the nearest - neighbor range of the interactions ) .
hence , they depend on the geometry and the spatial dimension of the lattice . to begin with ,
let us consider one - dimensional lattices ( spin chains , @xmath183 ) , for which the coordination number @xmath184 . substituting this value in the previous equations
, we immediately recover , among others , the value of the factorizing field and the ground state energy obtained analytically by kurmann _
et al . _ for the @xmath1 spin chain using the direct method @xcite .
et al . _ had to restrict their investigation to coupling constants @xmath185 taking values in the interval @xmath21 $ ] , thus establishing ground state factorization only in a restricted set of conditions .
instead , the general method allows to assess the occurrence of factorized ground states associated to different types of magnetic orderings in the entire region of the hamiltonian parameters .
one of the most important advantages of the general method is that it is easily applicable irrespective of the spatial dimensionality , while the direct method becomes increasingly difficult to apply with growing spatial dimensions of the lattice .
numerical techniques can overcome some of the limits of the direct method , although they fall short of the power and extension of the general analytic method .
the most important numerical result has been obtained by applying quantum monte carlo techniques to the study of the @xmath1 model in @xmath186 space dimensions ( square lattice ) @xcite . applying the present method to the @xmath1 model on the square lattice
allows not only to immediately reobtain analytically the previous numerical findings , but also , again , to extend the investigation to the entire space of the hamiltonian parameters by simply taking into account that for the square lattice @xmath187 .
moreover , the flexibility of the general method allows to include in the analysis different planar models .
for instance , considering models defined on planar hexagonal lattices , the previous analysis needs to be modified only in that the value of the coordination number , in this case is @xmath188 .
let us now move to @xmath189 and investigate ground state factorization in nearest - neighbor @xmath1 models on three - dimensional regular lattices , for which there were so far no results ( apart the preliminary ones obtained in our previous work @xcite ) either analytic or numerical ( such a task would be extremely demanding if approached either via the direct method or via numerical techniques ) . resorting to the general method , we can instead fully characterize the factorized ground state and the associated physical quantities just by entering the correct value of the coordination number , e.g. @xmath190 for a simple cubic lattice .
obviously , the results extend straightforwardly also for models defined on lattices of arbitrary higher dimension @xmath191 .
the analytic method can be applied with minor complications to models with interactions of longer spatial range . as already mentioned in sec .
[ sec.method ] , concerning ground state factorization , the main difference between models with short- and finite - range interactions is that for the latter there exist more than one equation in the set ( [ e.solcondition ] ) .
each of these equations is associated to a different value of the range @xmath127 , where @xmath125 is the maximum interaction range , that is the maximum distance between two directly interacting spins : all the couplings vanish for pairs of spins with inter - spin distance greater or equal to @xmath192 .    clearly , the fact that all the equations have to be solved simultaneously yields tighter conditions on factorizability compared to the case of models with nearest - neighbor interactions .
namely , it is not sufficient that the r.h.s . of each eq .
( [ e.solcondition ] ) takes values in the interval @xmath21 $ ] .
it is also necessary to require that all these values coincide .
this constraint yields a set of further conditions that must be satisfied by all the nonvanishing couplings .
for instance , let us suppose that the system is in an antiferromagnetic ordered phase along the @xmath69 direction . labeling by @xmath193 the value of @xmath131 and imposing that for every @xmath49 the r.h.s . of eqs.([e.solcondition
] ) must equate @xmath193 , we find that @xmath194 where , to ensure the absence of frustration , @xmath195 must obey the following inequality : @xmath196 .
analogously , for an antiferromagnetic order along @xmath71 : @xmath197 and to exclude frustration we must also have that @xmath198 . in the presence of a ferromagnetic phase ,
respectively for an ordering along @xmath69 and along @xmath71 , the constraints and the no - frustration condition become @xmath199 eqs .
( [ e.exe.cantix][e.exe.cferroy ] ) generalize the instance analyzed in ref .
@xcite , in which we limited the analysis to the case of vanishing @xmath169 that , according to eqs .
( [ e.exe.cantix][e.exe.cferroy ] ) , implies @xmath200 for all @xmath127 .    especially in the presence of a ferromagnetic order in the lattice , eqs.([e.exe.cferrox][e.exe.cferroy ] )
allow us to simplify the expressions of both the factorizing field and the energy per site .
namely , taking into account the definition of the ferromagnetic net interactions eq .
( [ e.netinteractions ] ) and the expression of the factorizing field in the presence of a ferromagnetic order ( see tab .
[ t.summary ] ) , we have @xmath201 respectively for an order along the @xmath69 and the @xmath71 directions .
similarly , for the energy density we obtain @xmath202    we now discuss a specific example ( out of many possible ones ) which demonstrates the power of the analytic method in conditions where the more conventional direct approach based on the factorized ansatz of refs .
@xcite is completely ineffective without the supplement of a remarkable intuition in guessing the correct forms of the possible factorized ground states .
let us consider a @xmath1-type spin model defined on a spatially isotropic cubic lattice of arbitrary size , possessing nearest - neighbor and next - nearest - neighbor interactions , and admitting arbitrarily different couplings along the @xmath69 , @xmath71 , and @xmath45 directions .
the coordination numbers associated to the two types of interaction can be straightforwardly computed from the geometry of the problem and are @xmath190 and @xmath203 .
suppose the model is deduced from an explicit experimental realization of a spin lattice and we are provided with explicit values of the various interaction strengths .
for example , @xmath204 .
the application of our `` on - demand '' ground - state factorization search engine is now immediate . from eq .
( [ fourcases ] ) we have that the candidate factorized state must possess antiferromagnetic order along the @xmath69 direction .
therefore the general conditions for the existence of a factorization point , as determined by eq .
( [ e.solcondition ] ) , can be summarized as follows @xmath205 this condition must hold together with the frustration - free constraint .
it is readily verified that @xmath206 provides a solution to the combined set of necessary and sufficient conditions for ground state factorization .
therefore if the experimentalist tunes the external field @xmath44 at the factorizing value @xmath207 , obtained from eq .
( [ e.factorizingfield ] ) , the engineered cubic spin lattice system will relax in a completely factorized ground state of the form given in table [ t.summary ] ( second row , last column ) with @xmath208 .
many other applications can be considered to models with arbitrary spin - spin interactions of longer range , and they can be treated with the same simplicity within the framework of the general analytic method .
for instance , a very interesting example of a model with infinite - range interactions will be treated in the next section . however , before moving to the study of this case , we should comment on the important issue of the interplay between factorization and frustration in quantum spin systems with multiple spatial scales of interaction .
the investigation of general models including many interaction terms of different spatial ranges needs to be carried out with particular care , because the interplay between the different interactions can lead in general to important effects of frustration that tend to suppress the occurrence of ground state factorization .
the very delicate and tricky nature of the problem is at the origin of a recent incorrect prediction about the existence of factorized ground states in one - dimensional @xmath1 models on finite rings with long - range interactions , obtained exploiting the direct method @xcite . in ref .
@xcite the author , resorting to the direct method , concludes that all the factorized energy eigenstates in one - dimensional long - range @xmath1 models , with interaction terms of spatial ranges that are in integer ratio with the total number of sites , are also the states of lowest energy .
however , the proof of this statement obtained using only the direct method requires that there must be no frustration in the system @xcite .
unfortunately , the interplay between the interactions of different spatial range can in general introduce frustration effects ( e.g. in the case of antiferromagnetic nearest - neighbor and next - to - nearest neighbor interactions , and so on ) .
the existence of these effects nullifies all proofs , based on the direct method , of the existence of fully separable ground states .
in fact , immediate counterexamples can be given of fully factorized energy eigenstates that are not ground states .
let us consider for instance an @xmath209 model on a finite ring of six spins with nearest - neighbor and next - to - nearest - neighbor antiferromagnetic interactions , and @xmath210 .
following ref .
@xcite , this model should admit a factorized ground state .
in fact , a factorized energy eigenstate with antiferromagnetic order does indeed exist , with an energy @xmath211 .
however , solving the model by exact diagonalization yields that the lowest energy is @xmath212 and that the associated ground state is entangled .
essentially the same type of counterarguments nullifies a recent claim that the factorized energy eigenstates are the factorized ground states of mixed - spin models of ferrimagnetism of the @xmath1 type @xcite .
a general analysis of ground state factorization in ferrimagnetic models with mixed spin-@xmath0-spin-@xmath182 interactions will be presented elsewhere @xcite , based on a proper extension of the formalism of squos to include spin-@xmath182 systems , as sketched in the final part of ref .
@xcite .    in conclusion
, the study of factorization in quantum spin models with many different finite- and long - range interaction terms requires a very careful analysis and the use of the general analytic method whenever frustration effects are present
. a rigorous and systematic study of the crucial ( and subtle ) interplay between frustration and factorization in quantum spin systems will be the subject of a forthcoming paper @xcite .
a very interesting limit of the cases analyzed in subsec .
[ sec.long ] is given by models with infinite range interactions , such as the fully connected or lipkin - meshkov - glick ( lmg ) model @xcite , extensively studied in condensed matter physics @xcite , which is described by the following hamiltonian @xmath213 with @xmath214 .
the pre - factor @xmath215 ( with @xmath2 denoting the total number of spins ) ensures the linear divergence of the ground state energy as @xmath2 increases .
let us remark that normalizing the interaction part of the lipkin - meshkov - glick hamiltonian requires some care : being a connected interaction , it scales as @xmath216 , i.e. as the product of the total number of spins and the number of pairs of spins .
therefore , to ensure that the interaction energy is extensive , i.e. scales with @xmath2 , the interaction part of the hamiltonian must be normalized by @xmath29 .
usually , for calculational ease , the interaction part is normalized by @xmath2 , so that it scales as @xmath29 . using this convention is strictly speaking incorrect , but is clearly harmless for large @xmath2 and in the thermodynamic limit .
for instance , concerning the separability of the ground state , it introduces a slight , spurious @xmath2-dependence of the factorization point through the multiplicative factor @xmath217 , and the error vanishes in the thermodynamic limit .
however , on fundamental physical grounds , the factorization points of translationally invariant systems , if they exist , must be independent of the size of the system , and must take the same value both at finite @xmath2 and in the thermodynamic limit . applying the kurmann - thomas - mller factorized ansatz @xcite and the @xmath218 normalization , dusuel and vidal obtained an analytic expression of the factorizing field for the lmg model that is @xmath2-dependent through the multiplicative factor @xmath217 @xcite .
although defined according to an anomalous convention , this result is essentially exact in the limit of large @xmath2 .
we first observe that the lmg model in eq .
( [ e.exe.hlmg ] ) can be obtained as the limiting case of a generic model described by eq .
( [ e.hamiltonian ] ) assuming that the range of interactions diverges ( i.e @xmath219 ) while the individual strengths become independent of the distance ( @xmath220 ) , scale with @xmath2 ( @xmath221 ) , and vanish along the @xmath45-direction ( @xmath222 ) . therefore , assuming @xmath223 , @xmath224 and @xmath222 in eq .
( [ e.hamiltonian ] ) we recover eq .
( [ e.exe.hlmg ] ) .
because of the sign of the interactions and of the range of values that can be assumed by @xmath225 , we immediately obtain that the energy per site ( [ e.energy ] ) is minimized by a state possessing ferromagnetic order along the @xmath69 direction .
moreover , taking into account that the coupling strengths are proportional to the ratio @xmath218 , we have that the net interactions converge to finite values in the thermodynamic limit : @xmath226 , @xmath227 .
it is then rather straightforward to apply the analytic method to the lmg model and hence prove , exactly , that the ground state is fully factorized at the factorizing field @xmath228 , which is , correctly , independent of @xmath2 .
the form of the fully separable ground state is the tensor product of local states of the form @xmath229 with @xmath230 and the energy density reads @xmath231 .
these results are in agreement with those obtained , using the direct method , by dusuel and vidal @xcite , modulo a proper normalization factor .
the versatility and generality of the analytic method can be demonstrated further by considering extensions of the translationally - invariant models discussed in the previous sections .
in fact , till now , we have always considered systems in which the coupling amplitudes depend only on the distance between the spins involved . in the present subsection
we shall focus instead on models which contemplate the possibility of spatial anisotropies , i.e. given the spin @xmath232 associated to site @xmath40 , the interaction coupling with spin @xmath233 associated to site @xmath109 depends not only on the distance @xmath234 but also on the location of site @xmath109 relative to site @xmath40 .    to simplify the notations , the physical analysis , and the number of possible cases that need to be considered , we limit ourselves to models with short - range interactions , i.e. with all the couplings vanishing for @xmath235
. we will comment briefly on general models with spatial anisotropies and interactions of arbitrary range at the end of this section .
let us assume that for each site @xmath40 the total number of nearest neighbors is @xmath236 , where @xmath237 is the number of nearest neighbors whose coupling strength with @xmath232 takes a certain value @xmath238 , and @xmath239 is the number of nearest neighbors whose coupling with @xmath232 is given by a different value @xmath240 .
obviously , we could consider situations more complicated at will , with an arbitrary number of @xmath241 different types of nearest neighbors .
for the given example , the system hamiltonian reads @xmath242 , \nonumber \\\end{aligned}\ ] ] where the index @xmath243 runs over the first @xmath237 nearest neighbors , and the index @xmath244 runs over the remaining @xmath239 nearest neighbors for each spin @xmath232 . according to the different values of @xmath237 and
@xmath239 this hamiltonian describes various possible models .
for instance , if @xmath245 and @xmath246 eq .
( [ e.exe.hamiltonian.anis ] ) describes @xmath1 models on ladders constituted by two coupled linear chains ; the case @xmath247 and @xmath246 corresponds to models defined on two coupled square lattices . in both cases ,
every spin in each chain ( plane ) , besides the usual interactions within a single chain ( plane ) involving @xmath237 nearest neighbors , is also coupled to @xmath239 companion spins located on the second chain ( plane ) .
the instance @xmath248 corresponds to a @xmath1-type model defined on a linear chain with nearest - neighbor interactions of alternating strengths .
the factorization properties of this latter model have been studied by g. l. giorgi in the limit @xmath249 that , belonging to the @xmath209 symmetry class , can be solved exactly @xcite .    as we have seen in sec .
[ sec.method ] , the first step of the method is to single out the magnetic order that minimizes the energy associated to the candidate factorized ground state , at fixed magnetization along the @xmath45 direction ( the direction of the external field ) .
provided there are no effects of frustration , in the present case we obtain that every possible hypothetical factorized ground state can assume one of eight different magnetic orders .
these eight possibilities stem from the four possible orders of eq .
( [ fourcases ] ) , that the system can support when the effects of @xmath240 can be neglected , duplicated according to the two possible arrangements , parallel or antiparallel , induced by the interactions associated to @xmath240 .
the four possible orderings along a single direction are represented in fig .
[ f.disegno ] for a model of two transversally coupled one - dimensional spin chains .     or
@xmath71 ) characterizing the candidate factorized ground states of models with spatially anisotropic interactions .
the case considered corresponds to a ladder of two coupled one - dimensional spin chains with coordination numbers @xmath245 and @xmath246 .
topmost row : parallel ( a ) and antiparallel ( b ) ferromagnetic order . bottommost row : parallel ( c ) and antiparallel ( d ) antiferromagnetic order.,width=321 ]    by imposing the minimization of the energy per site , eq .
( [ e.energy ] ) , the kind of order actually present in the ground state can be determined comparing the antiferromagnetic and ferromagnetic net interactions along the @xmath69 and @xmath71 axes , defined as @xmath250 where the @xmath84 sign discriminates between parallel / antiparallel ordering with respect to @xmath251 . as a function of the net interactions ,
the type of magnetic order present in the system is determined by the value of the minimum , i.e. of @xmath252 .
accordingly , we have the following correspondences : @xmath253 obviously , also in the presence of spatial anisotropies we may observe that , in terms of the net interactions , the fact that in the system there is no frustration ( and especially no frustration arising from the spatial anisotropies ) implies that latexmath:[$\mu = - z_1^{(1 ) }    according to which axis is characterized by the given magnetic order .    without loss of generality ,
let us consider the situation in which the system realizes a parallel antiferromagnetic order along the @xmath69 direction , see fig .
[ f.disegno](c ) .
all other cases can be treated analogously .
following once again the steps described in sec .
[ sec.method ] , and in complete analogy with eq .
( [ e.condition ] ) , it is possible to derive the set of conditions that the orientation @xmath86 of the e - squo must satisfy simultaneously in order for the system to admit a factorized energy eigenstate : @xmath255 for the two conditions to be satisfied simultaneously the two r.h.s . must coincide and must assume values in the interval @xmath21 $ ] .
hence , by equating them it is possible to derive an equation that plays the same role of eqs .
( [ e.exe.cantix][e.exe.cferroy ] ) . in the hypothesis that the model under investigation satisfies both eqs .
( [ e.exe.theta.anis.1][e.exe.theta.anis.2 ] ) , what is left to prove is that the associated factorized state be indeed the ground state . following the route illustrated in sec .
[ sec.method ] it is not difficult to verify that the following inequalities must hold : @xmath256 it is straightforward to establish that they are both satisfied , provided there are no frustration effects .
namely , for spins interacting via couplings of the first kind @xmath238 ( referring to fig .
[ f.disegno ] , these are the intra - chain interactions ) , the second of eqs .
( [ e.exe.suffcondition.anis.1 ] ) yields that the denominator in the r.h.s . of eq .
( [ e.exe.theta.anis.1 ] ) is non - negative and , therefore , the corresponding numerator must be non - negative and not exceeding the denominator .
hence , we obtain @xmath257 . on the other hand , from the first of eqs .
( [ e.exe.suffcondition.anis.1 ] ) it follows that @xmath258 which , in turn , implies @xmath259 .
the latter relation is in agreement with the hypothesis of absence of frustration .
furthermore , for pairs of neighboring spins interacting with a coupling of the second kind @xmath240 ( referring to fig .  [ f.disegno ] , these are the inter - chain interactions ) , the second of eqs .
( [ e.exe.suffcondition.anis.2 ] ) yields that the denominator in the r.h.s . of eq .
( [ e.exe.theta.anis.2 ] ) must be non - positive .
hence , also the corresponding numerator must be non - positive and therefore @xmath260 . simultaneously
, the first of eqs . ( [ e.exe.suffcondition.anis.2 ] ) implies @xmath261 and therefore @xmath262 , again in complete agreement with the hypothesis of absence of frustration . collecting all these results
we can conclude that in the absence of frustration , if a system described by the model hamiltonian eq .
( [ e.exe.hamiltonian.anis ] ) satisfies simultaneously all conditions in eqs.([e.exe.theta.anis.1][e.exe.theta.anis.2 ] ) , there exist a factorizing field @xmath63 at which the ground state is fully factorized and characterized by a parallel antiferromagnetic order along @xmath69 .
the factorizing field is determined according to the procedure described in sec .
[ sec.method ] , and it is not difficult to verify that it reads : @xmath263 the exact form of the factorized ground state reads @xmath264 in eq .
( [ e.antiferroparrallelolungox ] ) , the superscript `` @xmath265 '' ( `` @xmath266 '' ) denotes the upper ( lower ) chain of the ladder model considered in fig .
[ f.disegno](c ) .
we observe that @xmath267 does not depend on the choice of the chain , implying that the two generic @xmath4-th spins in both legs of the ladder are in the same state equipped with an antiferromagnetic order with respect to the nearest neighbor sites along the corresponding chain .
analogous results hold , with the appropriate trivial modifications , when considering the other possible magnetic orders compatible with factorization ( see eq .
( [ e.exe.eightcases ] ) ) .
the present work has been motivated both by the need for exact solutions to complex many - body quantum spin models , in particular by the question `` when is a mean - field solution exact ? '' , and by the necessity to acquire improved knowledge and control on the structure of quantum correlations for potential technological applications of spin systems .
we have introduced a simple and rather powerful all - analytic general method to determine the conditions for the existence and the properties of fully factorized ( fully separable ) ground states in translationally invariant quantum spin models with general two - body exchange interactions and subject to external fields .
the theory developed in the present paper builds on and extends the scheme originally introduced in ref .
it provides a readily useful procedure to establish , in terms of the hamiltonian parameters , whether the ground state is entangled or completely factorized and , in the latter case , what is the exact form of the state , the expression of the ground - state energy and of the magnetic observables , and the type of ordered phase ( magnetic order ) compatible with factorization .
the task is achieved by exploiting tools imported from quantum information theory such as the formalism of single - qubit unitary operations introduced in refs . @xcite and the corresponding entanglement excitation energies , whose property of vanishing if and only if a quantum ground state is totally factorized plays a central role in the analysis of ground - state separability . using these tools
we have completed the program initiated in ref .
@xcite by determining a general set of exact relations and inequalities that , for frustration - free systems , allow to establish rigorously the occurrence of all types of possible factorized ground states .
we have applied the analytic method to spin-@xmath0 models with general heisenberg - like exchange interactions of arbitrary range ( short , finite , long ) , either isotropic or anisotropic , defined on regular lattices of arbitrary dimensions .
the method is insensitive to the size of the system and applies rigorously in the thermodynamic limit as well as for finite lattices .
the method allows to establish novel exact , fully factorized ground state solutions of generally non exactly solvable models , corresponding to nontrivial sets of values of the hamiltonian parameters .
key results include , for instance , the rather straightforward determination of the existence and form of factorized ground states in @xmath1-type models with different types of short - range and long - range interactions on cubic lattices , a task that would be of formidable complexity if tackled with numerical techniques or resorting to the simple - minded direct method based on the product ansatz .    the complete factorization diagram for nearest - neighbor heisenberg - like anisotropic models might be particularly useful for those technological implementations which employ spin systems as information processors exploiting their ground state entanglement @xcite .
we have provided in fig .
[ f.factor ] a `` minefield '' map of the unentangled working points that need to be avoided , when engineering the couplings and tuning the external fields , in order to achieve satisfactory quantum transmission performances . on the other hand , there are also alternative tasks in quantum information which instead require as initial working points factorized states , e.g. for engineering , via suitable dynamics , specific classes of long - distance entangled states useful for quantum state transfer and quantum teleportation @xcite or the instantaneous creation of strongly entangled graphs or cluster states involving a mesoscopic or macroscopic number of spins for one - way quantum computation @xcite : in these instances the minefield may turn into a treasure map . the analysis of ground - state factorization would be particularly useful in the study of quantum spin models on open - end lattices .
such systems are not translationally invariant , and therefore the analytic method requires to be suitably extended to this type of instances .
this generalization is under way and should be in reach in the near future .
notably , according to a general theorem by kurmann , thomas , and mller on factorization @xcite , given any hamiltonian of the form eq .
( [ e.hamiltonian ] ) involving higher spins @xmath268 , the ground state of that spin-@xmath269 system is fully factorized at the same value of the external field @xmath101 eq .
( [ e.factorizingfield ] ) , at which factorization occurs in the corresponding spin-@xmath0 model .
therefore , the method and the results derived in the present paper have a much wider scope of application and can be straightforwardly generalized to interacting systems with arbitrarily high values of the spin , provided they are of the same hamiltonian structure as in the spin-@xmath0 case .
absence of frustration effects , either of geometric or dynamical nature , has been essential to the theoretical scheme derived in the present work .
for instance , we have determined the factorization diagram for specific examples of frustration - free systems endowed with isotropic interaction terms of different spatial ranges , and for systems with nearest - neighbor anisotropic interactions .
however , it is clear that there are many extremely interesting instances of more complex systems that would in general be subject to frustration effects .
continuing a step - by - step strategy of applying the general method to models of increasing complexity , natural further stages of investigation would include models subject to frustration , e.g. , models on linear chains with competing isotropic interactions of different range such as a 1-d @xmath1-type model with nearest - neighbor ferromagnetic interactions and next - to - nearest neighbor antiferromagnetic interactions . increasing further the degree of complexity , one could then consider models with competing anisotropic interactions of different ranges , that according to the type of anisotropy could be subject to both geometrically- and dynamically - induced frustrations .
when considering frustrated models , except for very special cases , even when the considered system admits a fully factorized state as an energy eigenstate , the set of conditions derived in the present work are not sufficient to establish that such an eigenstate is indeed the ground state of the system .
in fact , a competition arises between factorization , which requires a magnetic ordering , and frustration , which tends to destroy magnetic orders , ultimately leading to a chaotic correlating behavior .
the latter unavoidably tends to enhance ground - state entanglement , thus suppressing the occurrence of situations where the most favorable state , in terms of energy content , is completely uncorrelated .
in such cases , a naive application of the direct method based on the factorized ansatz , for instance to models including anisotropic interaction terms of different magnetic types on many different spatial scales ( or even on all scales ) , leads unavoidably to overestimate the factorizing effect of the balancing between interactions and external fields and thus to badly incorrect predictions on the occurrence of fully factorized ground states in these classes of models , as done , unfortunately , in the second part of ref .
@xcite . in a forthcoming work
@xcite we will present a systematic and thorough study of some important classes of frustrated quantum spin models in order to establish the conditions for the occurrence of true ground state factorization and to characterize it in terms of compatibility thresholds with frustration .
an intriguing issue that arises in connection with frustration is the relationship between factorization and the existence of ordered phases . for frustration - free hamiltonian systems that conserve the even - odd parity ,
spontaneous symmetry breaking and the existence of an ordered phase are necessary conditions for ground state factorization , because factorized states have no definite parity .
however , in the presence of frustration there may exist factorized energy eigenstates without ground state factorization
. it would then be very interesting to understand whether in these situations there is still a causal relation between factorization and the existence of phase transitions in the system .
factorization might then be used as an heuristic tool to gain insight in the phase diagram of frustrated quantum systems .
it would also be important to understand whether chiral interactions can enhance factorization beyond the limits imposed by the presence of frustration effects .
in fact , since factorization can be seen as `` mean field becoming exact '' , it seems to require quite naturally a balancing between interactions and external fields as the only possibility for its occurrence .
however , one can not exclude , in principle , that other mechanisms might lead to the same effect without requiring the presence of external driving fields .
future investigations will be concerned with the application of the analytic method to the determination of the spectrum of factorization , that is the study of the occurrence of factorized excited states in different classes of quantum spin hamiltonians , and the implications of factorization diagrams on the structural and informational properties of quantum spin models for potential applications in quantum technology .
we acknowledge cnr - infm research and development center `` coherentia '' , infn , and isi foundation for financial support .
one of us ( f.i . ) , would like to thank ylenia dautilia for invaluable exchanges .
in sec .  [ sec.method ] we have determined the general conditions for ground state factorization by selecting as candidate factorized ground states only those characterized by a ferromagnetic or antiferromagnetic order along the @xmath69 or @xmath71 directions depending on the value of the minimum value @xmath270 in the set of the net interactions [ see eq .
( [ fourcases ] ) ] , and neglecting all other possible magnetic arrangements . in the following
, we will prove that indeed , given any spin hamiltonian of the type eq .
( [ e.hamiltonian ] ) , if a factorized ground state exists , then it is characterized by one of the four above - mentioned magnetic orders .
let us consider the most general factorized state compatible with the constraint of a site - independent magnetization @xmath51 : @xmath271 where the local spin state @xmath272 is eigenstate of @xmath273 with eigenvalue equal to @xmath0 .
let us suppose that there exists an hamiltonian of the type eq .
( [ e.hamiltonian ] ) that admits , for some particular value of the external field @xmath101 , the state in eq .
( [ e.candidate.gen ] ) as an eigenstate . under these hypotheses we must have that the projection of the factorized eigenstate onto the four - dimensional hilbert space of a pair of spins in the lattice must be an eigenstate of the pair - hamiltonian term @xmath102 defined as in eq .
( [ e.couple ] ) .
following the same strategy as in sec .
[ sec.method ] let us introduce for each site @xmath274 a set of three auxiliary mutually orthogonal spin operators among which one , that we name @xmath123 , coincides with @xmath275 : @xmath276 by inverting eqs .
( [ e.axayaz.gen ] ) , we can conveniently re - express the conventional single - spin operators as functions of the new set of operators @xmath113 : @xmath277 substituting eq .
( [ e.sxsysz.gen ] ) into the pair hamiltonian we obtain the expression for @xmath115 at factorization as a function of the sets of operators @xmath113 and @xmath116 , which reads    @xmath278 -2 h_f^r \sin \theta \left [ { a_{\underline{i}}^x \left ( a_{\underline{j}}^z -\frac{1}{2}\right ) + \left ( a_{\underline{i}}^z-\frac{1}{2}\right ) a_{\underline{j}}^x } \right ] \nonumber \\ & & + a_{\underline{i}}^x a_{\underline{j}}^x \left ( j_z^r-2 h_f^r \cos \theta\right)+ a_{\underline{i}}^y a_{\underline{j}}^y \left ( j_x^r \sin\varphi_{\underline{i } } \sin \varphi_{\underline{j } } + j_y^r \cos \varphi_{\underline{i } } \cos \varphi_{\underline{j } } \right ) - h_f^r \cos \theta \left ( { a_{\underline{i}}^z + a_{\underline{j}}^z } \right ) \nonumber \\ & & + \left ( a_{\underline{i}}^z a_{\underline{j}}^y \cos \theta + a_{\underline{i}}^x a_{\underline{j}}^y \sin \theta \right)\left(-j_x^r \cos \varphi_{\underline{i } } \sin \varphi_{\underline{j}}+ j_y^r \sin \varphi_{\underline{i } } \cos \varphi_{\underline{j } } \right ) \\ & & + \left ( a_{\underline{i}}^y a_{\underline{j}}^z \cos \theta + a_{\underline{i}}^y a_{\underline{j}}^x \sin \theta \right)\left ( -j_x^r \sin \varphi_{\underline{i } } \cos \varphi_{\underline{j}}+ j_y^r \cos \varphi_{\underline{i } } \sin \varphi_{\underline{j } } \right ) \ , , \nonumber\end{aligned}\ ] ]        > from eq .  ( [ e.pairhij.gen ] )
, one has that to ensure that the projection of the factorized state onto the four - dimensional hilbert space of spins @xmath232 and @xmath233 is an eigenstate of the pair hamiltonian , it is necessary to require that the terms associated to the operators @xmath281 , @xmath282 , @xmath283 and @xmath284 vanish .
this condition yields that @xmath285 solving eq .
( [ e.condition2.1.gen ] ) and eq .
( [ e.condition2.2.gen ] ) we obtain that if @xmath286 , then the phases @xmath287 obey to the following relations : @xmath288 with @xmath241 integer .
it is immediate to verify that these relations are consistent only with a ferromagnetic or with an antiferromagnetic order along the @xmath69 or the @xmath71 direction .
this concludes the proof .
g. vidal , j. i. latorre , e. rico , and a. kitaev , phys .
* 90 * , 227902 ( 2003 ) ; v. e. korepin , phys .
lett . * 92 * , 096402 ( 2004 ) ; g. c. levine , phys . rev . lett . * 93 * , 266402 ( 2004 ) ; g. refael and j. e. moore , phys .
93 * , 260602 ( 2004 ) ; p. calabrese and j. cardy , j. stat . mech .
* p06002 * ( 2004 ) ; m.  b. plenio , j.  eisert , j.  dreiig , and m.  cramer , phys . rev . lett . * 94 * , 060503 ( 2005 ) ; m.  cramer , j.  eisert , and m.  b. plenio , phys . rev
. lett . * 98 * , 220603 ( 2007 ) ; j. eisert , m. cramer , and m. b. plenio , arxiv:0808.3773 .
g. vidal , phys .
* 99 * , 220405 ( 2007 ) ; m. rizzi , s. montangero , and g. vidal , phys . rev .
a * 77 * , 052328 ( 2008 ) ; v. giovannetti , s. montangero , and r. fazio , phys . rev . lett .
* 101 * , 180503 ( 2008 ) .
s. bose , phys .
lett . * 91 * , 207901 ( 2003 ) ; m. christandl , n. datta , a. ekert , and a. j. landahl , phys .
lett . * 92 * , 187902 ( 2004 ) ;
d. burgarth and s. bose , phys . rev .
a * 71 * , 052315 ( 2005 ) ; a. kay , phys .
rev . a 73 , 032306 ( 2006 ) .
y. makhlin , g. schn , and a. shnirman , rev .
* 73 * , 357 ( 2001 ) ; v. giovannetti , s. lloyd , and l. maccone , science * 306 * , 1330 ( 2004 ) ; f. dellanno , s. de siena , and f. illuminati , phys
. rep . * 428 * , 53 ( 2006 ) .
m. j. hartmann , f. g. s. l. brando , and m. b. plenio , phys .
. lett . * 99 * , 160501 ( 2007 ) ; d. g. angelakis , m. f. santos , and s. bose , phys .
a 76 , 031805(r ) ( 2007 ) ; a. kay and d. g. angelakis , europhys . lett . * 84 * , 20001 ( 2008 ) .    a. srensen and k. mlmer , phys .
83 * , 2274 ( 1999 ) ; a. b. kuklov and b. v. svistunov , phys .
lett . * 90 * , 100401 ( 2003 ) ; l .- m .
duan , e. demler , and m. d. lukin , phys .
* 91 * , 090402 ( 2003 ) .