excitonic molecules - or biexcitons - are known and well studied since several decades in conventional semiconductors .
their formations due to the coupling of a pair of exciton and their photophysics properties have been well established both theoretically and experimentally @xcite .
the outstanding facts are related to the giant oscillator strength associated with the one exciton - biexciton transition ; we may cite , for instance , the giant two photon absorption , shown theoretically by hanamura @xcite and observed first on cucl by gale and mysyrowicz @xcite .
there , five or six orders of magnitude more than for typical two - photon interband transitions are observed .    in organic compounds ,
the interest for excitonic bound states - or excitonic complexes or last , to use a terminology coming from the bethe ansatz technique , excitonic strings - is rather new since the main studies started about five or six years ago .
they concern three types of compounds : some organic charge transfer solids ( ct ) , _ j _ or _ h_-aggregates ( _ j - h_a ) and conjugated polymers ( cp ) .
the organic ct solids considered here , consist of planar aromatic donor ( d ) and acceptor ( a ) molecules alternately arranged along a one - dimensional stack . a large distance between molecules of the same stack ( more than 3 @xmath0 )
gives a rather small hopping integral compared to the characteristic coulomb repulsion .
moreover , a very large distance between neighbouring chains ensures the quasi - one dimensionality of these systems .
a prototype is given by the anthracene - pmda ( pyromellitic acid dianhydride ) ; for this compound , by comparison between theoretical calculations on small clusters and experimental results of differential transmission spectroscopies , the existence of biexcitons - and even triexcitons - has been demonstrated recently @xcite .
this result stays , for the moment , the only clear evidence for excitonic strings in organic compounds .
_ j _ and _ h_-aggregates are also , in many cases , one dimensional stacks of organic molecules .
but , on the contrary to the previous class of compounds , they have no ionic character and the attractive exciton - exciton interaction is smaller .
they are typical organic solids in the sense that their low - lying excitations can be understood in terms of frenkel excitons .
these compounds are claimed to be the ideal ones to get biexcitons but , without convincing demonstration for the moment - possibly be due to a too small binding energy @xcite .
conjugated polymers are very complex compounds .
as far as we are concerned with low - lying excitations , we may consider them , at first approximation , as an one dimensional interacting electron gas , which strongly interacts with the lattice too @xcite .
moreover , for a more realistic description , in many cases , three dimensional effects and inter- and intra - chain disorder should be considered .
because of these difficulties , despite a large amount of work , even the nature of the primary excitations remains subject of intense controversies @xcite .
anyway , some photoinduced absorption ( pa ) experiments @xcite and two photon absorption ( tpa ) experiments @xcite show some features which are possibly due to biexcitons .
only a few theoretical investigations of excitonic strings in organic compounds have been done @xcite .
most of them are based on exact diagonalizations of small clusters using different kinds of model hamiltonians depending on the system under consideration @xcite . from these results and based on physical reasoning , it has been argued that , when biexcitonic states exist , the intensity of the optical one exciton - biexciton transition may decrease when the binding energy of the biexciton increases ; moreover , this intensity is shown to be less important than the one of the one exciton - two excitons ( two free excitons ) transitions and , even , is supposed to be independent of the system size for sufficiently large binding energies @xcite . therefore , if the tendencies observed on finite clusters remain the same for the infinite system , on the contrary to the situation observed in some inorganic semiconductors @xcite , it may be difficult to create biexcitons by tpa
. a more appropriated experimental technique to observe excitonic strings would then be the differential transmission spectroscopy used by the way in @xcite . throughout this manuscript
, we often refer to these numerical works for comparisons , especially concerning the intensity of the transition from one exciton to biexciton states .    in this paper
, we propose a model for studying excitonic complexes in organic compounds . based on physical arguments ,
we think , this model is relevant for the three classes of compounds enumerated above . for organic crystals ,
as the @xmath1 and @xmath2 aggregates , our model is equivalent , in the low energy sector , to an effective _ xxz _ spin one - half heisenberg model already used to study this problem @xcite . for the ct solids and cp
, some extensions are necessary because the proper excitons are not of frenkel type .
they are realized by introducing effective hard core boson particles with a finite extension on a lattice , which mimics the exciton extension .
next , we propose an effective hamiltonian in terms of these particles which captures the essential physical ingredients to describe excitonic strings in organic compounds .
this model is considerably easier to study than the traditionally employed models and therefore some analytical investigations are possible . in this paper , we study the proposed model for two and three particles in a finite and infinite lattice .
our main interests concern the behaviours of the binding energies of the excitonic strings with the parameters of the model and of the one exciton - biexciton transition oscillator strengths ; are they giant as in inorganic semiconductors or , on the contrary , small as suggested by exact diagonalizations on small clusters @xcite ?
this question is crucial in order to clarify the role of hypothetical biexcitons - or even more extended excitonic strings - in the photophysics of organic compounds .
the paper is organized as follow . in section
ii , we introduce our hard core boson hamiltonian . in section iii ,
iv and v , we study this hamiltonian for one , two and three particles respectively .
first , for each class of systems mentioned in this study , we emphasize the main characteristics which are important for our purpose .
_ j _ and _ h _ aggregates are typical organic crystals . the ground state of these systems is the tensorial product of the molecular ground states .
the low - lying excitations are frenkel - type excitons where the hole and the electron of a monoexcitation are located on the same molecular site @xcite .
we can draw a picture of these excitations in the following way : if we represent symbolically the ground state of such a quasi - one dimensional system as @xmath3 , where @xmath4 stands for a molecule in its ground state , the frenkel excitons are the bloch states of local configurations such as @xmath5 where the superscripts @xmath6 and @xmath7 are for the electron and the hole , respectively ; @xmath8 is then a monoexcited molecule .
organic ct solids are stacks of donor ( d ) and acceptor ( a ) molecules .
we already mentioned that they are , with a good approximation , quasi - one dimensional systems with a nearest - neighbour hopping term small compared to the coulombic term .
the strong coupling regime is then a good starting point to consider such systems .
moreover , if one consider only a short range coulombic term , their ground state would be a charge density wave depicted roughly by the unique configuration @xmath9 where d molecules get two electrons and a molecules none .
the low - lying excitations would then be the excitons where an electron is transfered from a d to an a molecule ; they would be bloch states of excited configurations such as @xmath10 ( more details could be found in ref @xcite ) .
we get a so - called charge transfer exciton ( cte ) . for a more realistic coulomb potential ,
the picture proposed above will remain reasonable , at first approximation .
conjugated polymers are quite different from the other systems just considered since , obviously , they are not molecular crystals . on the contrary ,
they are usually described as an one dimensional electron gas governed by the pariser - parr - pople ( ppp ) hamiltonian which includes both long range coulomb interaction and , semi - classically , electron - phonon interaction @xcite .
numerical studies of this hamiltonian give , for a proper choice of parameters , excitonic states of small radius , typically a few monomeric units @xcite .
on the other hand , the ground state may be analyzed in terms of basis sets completely localized on the monomers .
then , the ground state of finite clusters is shown to contain mainly intermonomer charge fluctuations of very short range - especially the nearest - neighbours ones @xcite . according to these specificities ,
an effective model has been proposed very recently @xcite .
starting from the ppp hamiltonian , this model reaches a kind of molecular - crystal description for the low - lying excitations in conjugated polymers . in this new approach the vacuum ( or ground state ) is still not the simple molecular - type of ground state , proper to the simpson related models @xcite , but a matrix - product like state composed by local configurations extended at most , over two monomers .
excitons are then composed mainly by local configurations where either the electron and the hole are localized on the same monomer or distant by only one monomer ; the weight of the configurations where the electron and the hole are more separated decreases exponentially in the wave functions . therefore , the radius of the exciton , @xmath11 , defined as an average distance between the hole and the electron is rather small .
if we assume @xmath12 , the following pictorial representation can be done : if we represent the ground state as @xmath3 , where @xmath4 is a monomer , the excitons will be the bloch states of such configurations @xmath13 .
because of the electron - hole symmetry , the charge transfer of the electron on the right side and on the left side must be combined either symmetrically or antisymmetrically .
in contrast to the two previous cases , excitons have here , a structure .    in conclusion , in these three cases , excitons are bloch states of very localized excitations .
the vacuum and the nature of these local excitations vary , of course , following the system under consideration : crystal molecular ground state and frenkel excitons for @xmath1 and @xmath2 aggregates , cdw and cte for ct solids and finally , matrix - product state and some kind of cte but including the two possible arrangements of the charges , for cp .
however , all these very localized excitations can be thought in terms of effective hard core bosons ( hcb ) extended over a few bonds - with an extension of the order of the exciton s radius .
we adopt such a view in this work , and describe the dynamics of these excitations within the same model which we introduce now .    in this study
, we consider only one dimensional lattices with n sites , a site being a molecular - site in the case of molecular crystals@xcite and a monomer - site in the case of cp @xcite .
as we already mentioned , this is a good approximation for a large class of compounds studied here .
however , the generalization to higher dimensions should be possible and the subject of further studies .
excitons are extended over @xmath11 sites , the radius of the excitons , which varies from zero to a few units in the cases considered here . in this study , to mimic the excitons , we introduce hcb particles extended over one bond ( @xmath12 ) .
this minimal extension is adopted more by a pictorial preoccupation but , also , in accordance with the work developed in ref @xcite .
indeed , at the thermodynamic limit where we will derive all the analytical formulae , this extension plays obviously no role , since it remains small compared to the system size .
moreover , for finite size cluster calculations , the qualitative behaviours shown in this paper could also be obtained by choosing @xmath14 , @xmath15 or even other small @xmath11 . the restriction to @xmath12 does not affect the conclusions but , if necessary , more ( or less ) extended particles could also be considered without any difficulties . working with extended particles , instead of particles localized on a site ,
does not imply technical complications .
the coulomb interaction between excitons must be included in order to describe excitonic strings .
for that purpose , the relative position between the electron and the hole , the two elementary constituents of an exciton , becomes important .
we introduce then two species of hcb extended over one bond :    * the so - called right - bosons ( @xmath16 ) which mimic nearest neighbour electron - hole pairs with the electron on the right ( see figure [ rlbosons].a ) . * the so - called left - bosons ( @xmath17 ) which also mimic nearest neighbour electron - hole pairs but with the electron on the left ( see figure [ rlbosons].b ) .    the @xmath16 and @xmath17 particles obey the following hamiltonian    @xmath18    @xmath19 and @xmath20 are the destruction ( creation ) operators of the @xmath16 and @xmath17 bosons at site @xmath21 , @xmath22 , the operators number of hcb @xmath16 ( @xmath17 ) at site @xmath21 .
@xmath23 is the excitonic energy ( @xmath24 ) .
the @xmath1 term is the hopping term of the hcb ( @xmath25 ) ; it is directly related to the exciton band width given by @xmath26 .
@xmath27 is an effective term for virtual interactions which couple locally @xmath16 and @xmath17 bosons ( @xmath28 ) .
the @xmath16 and @xmath17 particles are effective representations of more complex particles , invoking a more sophisticated level of description ; these `` high level '' particles possess a structure which can be very complicated as it is the case for most of the cp .
the interaction @xmath27 reflects the existence of virtual processes exchanging the relative position of the electron and the hole of an exciton ; this can be due , for instance , to kinetic terms .
indeed , by a second order process involving a n.n . hopping term , it is readily possible to exchange the relative position of the two building particles of an exciton .
the last term of ( [ hamiltonian ] ) , @xmath29 ( @xmath30 ) , is the interaction between excitons restricted to nearest - neighbours .
we have to distinguish between two situations : for two excitons of the same species , @xmath31 or @xmath32 , the nearest charges are of opposite sign so that the resulting interaction must be attractive ; the situation is reversed for two excitons of opposite species , @xmath33 or @xmath34 , where the resulting interaction is repulsive .    for structureless excitons , frenkel or ct - excitons , the definition of two species of degenerate hcb makes no sense . to study strings of such excitons with the general hamiltonian ( [ hamiltonian ] ) , it is sufficient to : ( i ) identify the @xmath16 and @xmath17 particles and ( ii ) to take the limit @xmath35 . to be more precise
, @xmath16 and @xmath17 particles may exist in the case of ct solids , but with a large energy difference which results , in our model , in a very small @xmath27 ; therefore , it is a good approximation to simply neglect this term in this case .    in the following , we will study the hamiltonian ( [ hamiltonian ] ) for one , two and three particles , which are the cases of experimental interest .
excitons appear as bloch states of symmetric or antisymmetric combinations of @xmath16 and @xmath17 bosons , corresponding to two different symmetry classes    @xmath36    where @xmath37 is the state without any boson .
the separation in energy between this two states is given by @xmath38 , the @xmath39 state being lower .
in the case of _ j _ ( _ h _ ) aggregates and ct solids , we identify @xmath16 and @xmath17 bosons and get only one excitonic state @xmath40    in all cases , exciton states result in a very intense peak in the linear absorption spectrum for @xmath41 , while the electron - hole continuum is almost not visible .
this very unusual feature , if we think about conventional semiconductors , is an important characteristic of these one dimensional compounds @xcite .
this is easily understandable since the radius of the exciton is very small and the ground state contains mainly short range charge fluctuations @xcite .
now , we write down the elementary transition moments due to a localized hcb , @xmath42 , in the case of cp and @xmath43 , in the case of _ j _ ( _ h _ ) aggregates and ct solids ( @xmath44 is for organic crystal ) ; these quantities will serve in the following to express the transition moments between one exciton and two exciton states    @xmath45    where @xmath46 is the charge of the electron and @xmath47 the position operator .
the intensity of the ground state - one exciton transition at @xmath41 is proportional to @xmath48 or @xmath49 , depending on the system under consideration @xcite . in the case of cp
, only the excitonic state @xmath50 is observable in linear absorption .
however , the state @xmath51 could become important , as it is the case in electroabsorption experiments for instance .
the wave functions with two excitons are written as @xmath52where @xmath53 and @xmath54 are @xmath16 or @xmath17 bosons , @xmath55 the ket relative to the situation where the particle @xmath53 is in @xmath56 and @xmath54 in @xmath57 with @xmath58 , @xmath59 is the corresponding amplitude .    as usual
, we treat separately the movement of the center of mass and of the relative position of the two particles .
moreover , we have to specify the nature , @xmath16 or @xmath17 , of the particles located on @xmath56 and @xmath57 , it is then convenient to use matrix formalism .
we write @xmath60_{n}\mid n_{1},n_{2}>\]]with @xmath61 ( @xmath62 ) and    @xmath63_{n}=\left ( \begin{array}{l}\varphi_{rr}(n)\\ \varphi_{rl}(n)\\ \varphi_{ll}(n)\\ \varphi_{lr}(n ) \end{array}\right)\ ] ]    where the small left ( right ) indices refer to the nature of the left ( right ) particle , @xmath16 or @xmath17 bosons .
we apply periodic boundary conditions ( pbc ) to the wave functions which give the following constraints    @xmath64_{-n+n}=[\varphi]_{n } \quad \mbox{where } \quad n = n_{2}-n_{1 } \end{array}\right.\ ] ]    moreover , we work with hard core bosons extended over one link so that we must preserve @xmath63_{0}=[\varphi]_{1}=[0]\]]where @xmath65 $ ] is the zeroth four component vector .    with the wave function ( [ wf ] ) ,
the eigenvalue equation is written @xmath66_{n}=-2j\cos q[i]\{[\varphi]_{n-1}+[\varphi]_{n+1}\}-\alpha[\sigma][\varphi]_{n}-\delta_{n,2}v[\lambda][\varphi]_{n}\]]where @xmath67 $ ] is the 4 by 4 identity matrix , and @xmath68=\left ( \begin{array}{c }   0 \quad 1 \quad 0 \quad
1 \\ 1 \quad 0 \quad 1 \quad 0 \\ 0 \quad 1 \quad 0 \quad 1 \\ 1 \quad 0 \quad 1 \quad 0 \end{array}\right ) \quad , \quad   [ \lambda]=\left ( \begin{array}{c } 1 \quad \quad \quad \quad \quad \quad \\ \quad \quad -1 \quad \quad \quad \quad\\ \quad \quad \quad \quad 1 \quad \quad \\
\quad \quad \quad \quad \quad \quad -1 \end{array}\right)\ ] ]    with the equation([vp ] ) , we have reached an impurity like problem where free particles , belonging to four different continuum , interact with impurities pinned at @xmath69 .
first , by applying a local unitary transformation , we solve ( [ vp ] ) away from the impurities .
@xmath70_{n}=[u][\varphi]_{n}\\ \quad \\
\left(\begin{array}{c}\varphi_{0}(n)\\ \varphi_{\bar{0}}(n)\\ \varphi_{+}(n)\\ \varphi_{-}(n)\end{array}\right ) =   \left(\begin{array}{crlrlrl }   \frac{1}{\sqrt{2 } } &     & 0                  & -   & \frac{1}{\sqrt{2 } } & & 0 \\   0                   &     & \frac{1}{\sqrt{2 } } &    & 0                   & -   & \frac{1}{\sqrt{2 } } \\
\frac{1}{2 }         &     & \frac{1}{2 }        &    & \frac{1}{2 }         & & \frac{1}{2 } \\   \frac{1}{2 }         & -   &   \frac{1}{2 }       &    & \frac{1}{2 }         & - & \frac{1}{2 }   \end{array}\right )    \left(\begin{array}{c}\varphi_{rr}(n)\\ \varphi_{rl}(n)\\ \varphi_{ll}(n)\\ \varphi_{lr}(n ) \end{array}\right ) \end{array}\ ] ]    we can now rewrite equation ( [ vp ] ) in terms of @xmath71_{n}$ ] ; we find @xmath72_{n}=-2j\cos q[i]\{[\bar{\varphi}]_{n-1}+[\bar{\varphi}]_{n+1}\}-2\alpha[\sigma_{2}][\bar{\varphi}]_{n}-\delta_{n,2}v\{[\sigma_{1}][\bar{\varphi}]_{n}+[\sigma_{3}][\bar{\varphi}]_{n}\}\]]where @xmath73=\left(\begin{array}{c}1 \quad \quad \quad \quad \quad \quad \\
\quad \quad -1
\quad \quad \quad \quad \\ \quad \quad \quad \quad   0 \quad \quad \\
\quad \quad \quad \quad \quad \quad 0\end{array}\right ) , \quad [ \sigma_{2}]=\left(\begin{array}{c } 0 \quad \quad \quad \quad \quad \quad \\
\quad \quad 0 \quad \quad \quad \quad \\ \quad \quad \quad \quad 1 \quad \quad \\
\quad \quad   \quad \quad \quad \quad -1\end{array}\right ) , \quad [ \sigma_{3}]=\left(\begin{array}{c}0 \quad \quad \quad \quad \quad \quad \\
\quad \quad 0 \quad \quad \quad \quad \\
\quad \quad   \quad \quad 0 \quad 1 \\
\quad \quad \quad \quad 1 \quad 0\end{array}\right)\ ] ]    without the impurity like states at @xmath69 , there is no mixing between @xmath74 , @xmath75 , @xmath76 and @xmath77 components ; with the scattering of particles at @xmath69 , only the @xmath76 and @xmath77 components are mixed .    from equations ( [ moment ] ) and ( [ unitary ] ) ,
it is clear that only the @xmath78 and the @xmath79 components of @xmath71_{n}$ ] are relevant to study the optical responses of two exciton states .
since we are mainly interested in spectroscopic properties of biexcitons , we consider only these two components in the following .
the @xmath80 and @xmath81 components have no physical meaning for structureless excitons , but could be important for cp ; their characteristics could be calculated , following exactly the same way as the one presented here .
we write , following h. bethe @xcite , @xmath82which give the following expressions for the energy where @xmath83 and where we choose @xmath84 as reference @xmath85    @xmath86 , @xmath87 , @xmath88 and @xmath89 , @xmath90 and @xmath91 are the constants to be determined . note that , if in ( [ fonctions ] ) @xmath90 and @xmath91 are purely imaginary numbers , we get free excitons .
on the contrary , if they are purely real numbers , we get bound - states , the biexcitons .
the ansatz ( [ fonctions ] ) gives obviously the exact solution of ( [ vp2 ] ) without impurities ( @xmath92 ) ; with impurities , ( [ fonctions ] ) must fulfill the following constraints in order to remain solution of ( [ vp2 ] ) @xcite @xmath93    moreover , the pbc imply @xmath94    in the following , we give the solutions of the above mentioned problem for different choices of parameters .      in this trivial case ,
the two excitons are free ; we get two degenerated continuum for @xmath78 and @xmath79 within the energy range @xmath95 $ ] .
it is easy to show that a very intense linear absorption takes place for @xmath96 and @xmath97 only ; then the transition moment of the one exciton - two free exciton transition at the edge of the continuum behaves as    @xmath98    plus corrections in @xmath99 . for @xmath100 and/or @xmath101 ,
this strong absorption remains but slightly shifted toward the high energies .
this transition is similar to the ground state - one exciton transition ; indeed , the excitation energies and the behaviour of the transition moments are the same in these two cases .
this characteristic was already used in @xcite to explain the lack of bleaching signal in pa experiments for charge transfer and frenkel exciton systems .
this case concerns ct solids and _ j _ or _
h _ aggregates for instance . without the @xmath27 term ,
our model is equivalent to a standard interacting frenkel exciton hamiltonian already studied numerically @xcite and analytically @xcite .
our results , obtained with a different technique , are in accordance with the conclusions reach by these various studies .    since @xmath102
, bound states - or biexcitons - may exist ; in this subsection we consider the biexcitonic states only and then @xmath90 and @xmath91 of equations ( [ fonctions ] ) , as real quantities .    when @xmath35 , equations ( [ energie ] ) imply that @xmath103 .
from equations ( [ impurities ] ) together with the pbc ( [ pbc ] ) , we get at the thermodynamic limit , the following set of equations @xmath104where @xmath105 .    additionally to the continuum ,
if @xmath106 , i.e @xmath107 , we get with the second equation of ( [ set ] ) and equations ( [ energie ] ) , two bound states with energies given by the following simple expressions @xmath108 by using a two - particle green function approach @xcite , these equations have already been obtained for @xmath96 in @xcite and for every @xmath109 in @xcite together with the biexciton condition @xmath110 .
the first equation gives a bound state below the continuum of two free exciton states , this is the biexciton of interest for us ; the second equation gives another bound state on top of the continuum . in the following ,
we consider only the first state and calculate its wave function and the transition moment associated with the one exciton - biexciton transition .
first , we write @xmath111 and @xmath112 . with the pbc ( [ pbc ] ) and equations ( [ set ] ) , we get the expression for @xmath113 @xmath114then , with equations ( [ fonctions ] ) and ( [ set ] ) , we get @xmath115)\]]valid in the interval @xmath116 .
@xmath117 is the normalization constant given by @xmath118^{-\frac{1}{2}}\ ] ]    the first equation of ( [ set ] ) shows that the wave function of the lowest bound state is a symmetric mixture of the @xmath76 and @xmath77 components resulting in configurations with two like particles @xmath31 or @xmath32 . on the contrary
, the highest bound state contains only configurations with two unlike particles @xmath33 or @xmath34 .
we recall that for this case @xmath16 and @xmath17 particles are identified so that the distinction between the case with @xmath31 or @xmath32 bosons and @xmath33 or @xmath34 bosons is then a way to treat the case of attractive and repulsive interaction between bosons at the same time . in case of attractive interaction ( two bosons of same species ) , we get a bound state below the two exciton continuum ; the situation is reversed in the case of repulsive interaction ( two bosons of opposite species ) , where a bound state appears above the two exciton continuum @xcite .
for large systems , the transition moment of the one exciton - biexciton transition is non - zero only for @xmath96 and @xmath41 ( @xmath109 , the center of mass of the biexciton , @xmath119 , the exciton momentum ) ; its expression is then given with the help of the elementary transition moment ( [ moment ] ) by ( more precisely , the square of the transition moment is shown here ) @xmath120going to the thermodynamic limit , we found the simple result @xmath121which is exactly the result found in @xcite by using a two particle green function approach .
we recall that the bound states ( biexcitons ) exist only if @xmath122 and that the biexciton binding energy increases when @xmath123 approaches @xmath124 . with this formula
, the behaviour of the absorption from the one exciton to biexciton state appears very clearly .
first les us consider its behavior when @xmath123 approaches @xmath124 or @xmath125 , at the thermodynamic limit ( equation ( [ momentthermo ] ) ) : for @xmath126 , the transition moment reaches saturation to @xmath127 which is easily understandable since , in this case , the two excitons are strongly bounded together in two n.n .
sites ; for @xmath128 , the transition moment diverges following the asymptotic behaviour in @xmath129 .
next , we can make some conclusions by studying the equation ( [ momentn ] ) .    * at large and fixed @xmath130 , with decreasing @xmath123 ( which means , with increasing the biexciton binding energy ) , the intensity of the transition decreases .
* at fixed @xmath123 and sufficiently large @xmath130 , @xmath131 is almost independent of the system size . on the contrary
, the intensity of the transition between one exciton and two free excitons increases linearly with the system size ( cf .
. [ momentfree ] ) .
however , as expected , when @xmath123 approaches @xmath125 , @xmath132 recovers progressively a linear dependence in @xmath130 .
both of these behaviours depict inverse mechanisms than the ones theoretically predicted for the ground state - one exciton transition @xcite , where the intensity of the excitonic peak is proportional to the size of the system and increases with the binding energy .
these behaviours are illustrated in figure [ intensity ] where the equation ( [ momentn ] ) is evaluated for @xmath133 .
there , we see that the intensity of the one exciton - biexciton transition decreases dramatically when @xmath123 decreases ( which means , when the binding energy increases ) reaching saturation for small @xmath123 . then , this intensity is independent of @xmath130 which gives in pa and tpa spectroscopy an intensity proportional to @xmath134 instead of @xmath130 for the transition to the continuum of two free excitons ; this last behaviour can be seen by comparison with equation ( [ momentthermo ] ) evaluated at the large @xmath130 limit ( dashed line ) .
consequently , at the thermodynamic limit , biexcitons may be observable by spectroscopy experiments , but with the need of very clean compounds , otherwise we may expect at first analyse , that the intense transition to two free excitons would not permit to detect possible biexcitons . within our model , no two - photon giant resonance @xcite can be expected for @xmath1 , @xmath2 aggregates and organic charge transfer solids .
the very same conclusions were already pointed out in ref @xcite from numerical calculations and in @xcite from analytical studies .
this case concerns cp . as in the previous subsection , we look for bound states only so that @xmath90 and @xmath91 are again assumed to be real quantities .
the question we want to address more specifically here is how the @xmath27 term will affect the binding energy , the wave function and the transition moment related to the biexcitons .
for @xmath136 , from equations ( [ impurities ] ) and ( [ pbc ] ) , we get at the thermodynamic limit the following set of equations    @xmath137    from equations ( [ energie ] ) and the second equation of ( [ set2 ] ) , we get @xmath138where @xmath105 and @xmath139 .
this expression together with equations ( [ energie ] ) , gives the energy of the lowest biexciton .
the binding energy of the biexciton , @xmath140 , is then given by the following equation @xmath141\]]the study of this expression shows that the biexciton binding energy decreases when @xmath27 increases ( see figure [ bindingenergy ] ) .
this behaviour is not surprising since the @xmath27 term mixed together @xmath16 and @xmath17 particles ; by doing so , configurations with two unlike particles , @xmath33 or @xmath34 , appear in the wave function which increase its energy and then reduce the binding energy with respect to the case for @xmath142 .
the critical value @xmath143 for which the binding energy becomes zero is given by the solution of the equation @xmath144 ,    @xmath145    the critical value , @xmath143 , increases when @xmath123 decreases , to diverge for @xmath146 : for a very large binding energy , @xmath27 has no effect on the bound states anymore .
another interesting quantity along this critical line , is the critical value @xmath147
@xmath148    from this equation , we can see without any surprise , than @xmath147 - the critical value above which one gets biexcitonic states - continuously increases when @xmath149 increases .
the attraction between excitons must be stronger when the mixing between @xmath16 and @xmath17 bosons becomes more important .
we analyze now the wave function and the resulting transition moment associated with the one exciton - biexciton transition at @xmath96 .
we proceed as in the previous case ; we write @xmath150 , @xmath151 and @xmath152 , @xmath153 and get with equations ( [ pbc ] ) and ( [ set2 ] )    @xmath154)+(-1)^{\delta_{\sigma_{1},\sigma_{2}}}e^{(\mu+\ln\lambda)(n-3)}\cosh([\mu+2\ln\lambda][n-\frac{n}{2}])\right)\ ] ]    where again @xmath155 .
@xmath156 is the normalization constant which could be easily calculated if necessary .
the first term of ( [ wf2 ] ) comes from the @xmath76 component , the second term from the @xmath77 component . for @xmath35 , @xmath157 and
the previous result is recovered . for @xmath158 ,
the weight of the @xmath77 component decreases smoothly in the wave function .
consequently , the a priori unfavourable configurations @xmath33 or @xmath34 appear in the biexciton wave function , on the contrary to the previous case without the @xmath27 interaction .
the transition moment of the one exciton - biexciton transition at @xmath96 can be easily calculated with the expression of the wave function ( [ wf2 ] ) , but contains no new information compared to what was already pointed out for @xmath35 .
the results obtained in the previous subsection remain valid in this case .
the only difference comes from the fact that the @xmath77 component is non active in linear absorption .
consequently , for the same binding energy , if we assume @xmath159 , the biexcitonic peak would be a bit more intense for @xmath1 ( @xmath2 ) aggregates or ct solids than for cp .
for instance , for @xmath35 , we get @xmath160 ; for @xmath158 , the @xmath76 becomes more important than the @xmath77 component in the wave function : the intensity of the one exciton - biexciton transition then increases .
now , let us take reasonable parameters for conjugated polymers .
extracting from experimental data ( we take roughly the binding energy of the @xmath161 exciton state observed in polydiacetylene compounds ) , we get @xmath162 @xcite . by calculating the coulomb interaction between two hcb using the mataga potential ( for instance ) @xcite and for a reasonable lattice constant ( @xmath163 ) ,
we get @xmath164 the @xmath27 term may be evaluated from the energy difference between the @xmath161 exciton and the @xmath165 `` ionic '' exciton .
these states are represented here , in an effective way , by @xmath50 and @xmath51 respectively ( see eq .
[ 1exciton ] ) .
this quantity could be estimated , for instance , from exact resolutions of the ppp hamiltonian for small clusters or , in a better way , from electroabsorption experiments . there
, the state @xmath166 shows a red - shift with no change of the spectral line and a quadratic dependence on the applied field .
this quadratic stark effect of the exciton can be reasonably well reproduced with the state @xmath51 only , neglecting , in particular , the continuum effects @xcite . within this assumption , it is then very easy to estimate @xmath27 .
we take here a reasonable value , @xmath167 . with this choice of parameters ,
we get a binding energy of @xmath168 for the biexcitons which is in accordance with experimental results @xcite . by the way , it is easy within our model to get such a value for the binding energy for other values of @xmath27 or / and @xmath1 by changing @xmath29 , which is more flexible .    in figure [ finitesize ] , a calculation for a finite cluster of ten sites is shown .
the parameters are @xmath169 and @xmath170 which give a typical spectrum .
be is for the peak associated with the biexciton and , fe , for the peak associated with the more intense two free exciton state ( even if such concepts are not well defined for finite clusters ) .
the be peak is a bit smaller than the fe one , which is in accordance with the results of ref .
@xcite from exact calculations using the extended hubbard hamiltonian , a short version of the ppp hamiltonian . at the thermodynamics limit
, only the fe peaks will survive , slightly shifted toward the low energies .    to conclude with the two particle case
, we may summarize our main results . for @xmath35
, biexcitons exist if @xmath122 ; the transition moment of the one exciton - biexciton transition shows inverse behaviour compared to the ground - state - one exciton transition : it decreases when the biexciton binding energy increases and , for sufficiently large @xmath130 , is independent of the system size . for @xmath171 , the conclusions aforementioned remain valid at the condition that @xmath172 for @xmath1 and @xmath29 fixed . above this critical value
, biexcitons fall in the continuum .
finite size cluster calculations give qualitatively similar results as other calculations performed with more complicated hamiltonians .
for more than two particles , the question about the integrability of the hamiltonian ( [ hamiltonian ] ) becomes important . to answer this question ,
it is sufficient to analyze carefully the case with only two particles @xcite . for that purpose ,
let us first rewrite the hamiltonian ( [ hamiltonian ] ) in a more convenient way by introducing two new operators    @xmath173    then , the hamiltonian becomes diagonal in @xmath27 @xmath174to illustrate this new picture , note that the states with one exciton ( [ 1exciton ] ) are the bloch functions of one @xmath175 or one @xmath176 particle .
let us now consider the case with two like particles ( @xmath177 or @xmath178 ) located in @xmath56 and @xmath57 ( @xmath58 ) . because the last term of ( [ hamiltonian2 ] ) exchanges the color ( @xmath175 or @xmath176 ) of the two particles , the corresponding amplitude in the wave function
is written as @xmath179the constants @xmath180 being complex numbers determined in the previous section by solving the schrdinger equation .
the system is invariant by translation , hence the total momentum is conserved : @xmath181 .
we then write @xmath182 , @xmath183 and @xmath184 , @xmath185 .
now , let us do a scattering experiment. starting , for instance , with two @xmath175 particles , after the scattering processes we get two @xmath176 particles .
the collision being elastic , the energy ( [ energie ] ) must be conserved @xmath186    depending on the @xmath27 value , we can make some conclusions .    * if @xmath100 , @xmath187 ( see eq .
[ wf2 ] for instance ) .
hence , there is a change in the momenta during the scattering processes . in other words ,
diffractive processes are involved during the scattering which is sufficient to conclude that the model defined by ( [ hamiltonian ] ) ( or ( [ hamiltonian2 ] ) ) is non - integrable @xcite . * if @xmath35 , @xmath188 .
there is no diffractive processes anymore ; then , as we look for the lowest bound - states , we will see that the model ( 1 ) may be expressed , in that case , as an effective spin-1/2 _ xxz _ hamiltonian which is well known to be integrable @xcite . rigorously , to conclude about the integrability of this model in the general case , the study of the two particles @xmath189 matrix has to be done @xcite .
we will not attempt such study here .    in the following two subsections
, we will study both cases , starting with the simplest one where @xmath35 .      for @xmath35 , by interacting with each other , particles
do not exchange their color ; the interaction is reduced to an exchange of momenta alone . in such case , it is better to consider the hamiltonian ( [ hamiltonian ] ) expressed in terms of @xmath16 and @xmath17 particles .    since @xmath16 and @xmath17 particles are hcb , an _
n_-electron configuration is subdivided into _
n _ separated regions .
each particle stays in its own area of the chain ; within optical terminology , only reflections are allowed .
the hamiltonian is studied independently for each configuration of color which corresponds to a configuration of pair - interactions , attractive or repulsive , depending on the nature of the neighbouring particles . in the general case ,
the problem is very complex .
however , the interesting case for us is given by the configurations where all the particles have the same color . in this case
, the model can be written in the form of a _ xxz _ ferromagnetic hamiltonian for spin-1/2 which has been solved exactly by using the bethe - ansatz @xcite .
first , our goal is to show the equivalence of our model ( [ hamiltonian ] ) in the low energy sector and the _ xxz _ model . for simplicity
, we consider here particles without extension ( @xmath14 ) .
of course , it is possible to work with our defined @xmath16 and @xmath17 particles , but the comparison with the _ xxz _ model would then claimed some additional notations unnecessary for our purpose . by the way , in the second subsection of this section , we will show more formally the equivalence between the two models working with the extended particles .
second , since we consider the simplest case with _
n_-particles like - for instance , the states with _
n _ _ r _ particles - the hamiltonian of the system may be expressed in terms of spin - one - half operators : @xmath190 being for a site @xmath21 without any particle and @xmath191 , for a site occupied by a @xmath16 particle ( an exciton ) .
the model becomes then equivalent to the spin-1/2 @xmath192 model with an attractive anisotropic term @xcite    @xmath193    the ground state of ( [ xxz ] ) is the state with only up - spins , @xmath194 .
@xmath195 creates a down - spin ( a @xmath16 particle or exciton ) localized at site @xmath21 , @xmath196 .
we introduce here , as usual , the parameter of anisotropy @xmath197 @xcite .
this model has been subject of intense studies since several decades , and many results are known about it @xcite . concerning the bound states , or more precisely , the so - called `` string - states '' , their properties depend on the anisotropic parameter @xmath198 .
within the `` string hypothesis '' , the following results were found@xcite .    * if @xmath199 , the bound states of @xmath21-spins , here @xmath21-excitons , exist without any restriction on @xmath21 @xcite .
the energy of a @xmath21-excitonic string is then given by the following formula first appeared in @xcite @xmath200with @xmath201 and @xmath109 being the wave number of the n - string center of mass .
for @xmath69 , the energy ( [ energiebiexciton ] ) is of course recovered . *
if @xmath202 , the @xmath21-excitonic strings may exist but with a strong restriction on @xmath21 @xmath203with @xmath204 , @xmath205 .
moreover , these excitonic strings fall in the continuum spectrum .
we have seen in the previous section , that biexcitons exist if @xmath199 . in this case , from the above conclusions about the _
xxz_-model , we can say that @xmath21-excitonic strings will also exist and that their binding energies will increase with @xmath21 ( [ strings ] ) ; this result is in agreement with numerical studies on small clusters @xcite .
however , as we have shown in detail for the biexcitonic case , it is not easy to observe such states experimentally , since the particles are tightly bounded inside the complex , resulting in a lack of oscillator strength . for @xmath206
, excitonic strings may still exist but with an energy falling in the continuum . in this limit ,
the excitons are less bounded and , therefore , it could be interesting to study the signatures of these states under the influence of an electric field .
we leave this question to further work .      for @xmath100 ,
the interactions between particles involve both exchange of momenta and exchange of color . by a simple argument , showing the appearance of diffractive terms in the two particles scattering case ,
the non - integrability of the model ( [ hamiltonian2 ] ) has been shown .
therefore , to study _ n_-excitonic strings with @xmath207 one needs to do some approximations . in the following ,
we propose a trial wave function for the three particle case build from the bethe ansatz solution of the _ xxz_-heisenberg model within the string hypothesis .
let us consider the expression ( [ hamiltonian2 ] ) of our model expressed in terms of @xmath176 and @xmath175 particles .
for the three particle case , eight configurations have to be distinguished depending on the color of the particles and on their relative positions ; we denote these configurations as @xmath208 , where @xmath209 is for the color ( @xmath176 or @xmath175 ) .
these configurations are separated into two disconnected channels : on one hand , we have @xmath210 , @xmath211 , @xmath212 and @xmath213 ; on the other hand , @xmath214 , @xmath215 , @xmath216 and @xmath217 .
we consider only the latter channel which is the lowest in energy .
the wave function for the triexciton is expressed as a linear combination @xmath218where @xmath219 stands for the @xmath214 , @xmath215 , @xmath216 and @xmath217 states respectively ( @xmath220 for _ right _ , _ center _ and _ left _ refer to the position of the @xmath175 particle in the three particles complex ) .
@xmath221 are the coefficients to determine . every component @xmath222 is assumed to have the bethe ansatz structure @xmath223where , more precisely , the particles are again separated by one bond , at least , and with @xmath224the @xmath225 are the quasi - momenta and the @xmath226 the phase terms @xcite .
the summation runs over the @xmath227 permutations of the indices .
we apply periodic boundary conditions to the wave functions , @xmath228 , and we get the well known equations @xcite    @xmath229    since we are looking for excitonic complexes , we are interested in complex solutions of these equations . following ovchinnikov @xcite , we assume that the only relevant phase terms are @xmath230 and @xmath231 , then @xmath232where @xmath233 are integers between @xmath124 and @xmath234 .
it follows @xmath235    next , we assume @xmath236 and @xmath237 without loss of generality ; then , at the thermodynamic limit , among the summation ( [ ba ] ) only one term remains with an exponential accuracy in @xmath130    @xmath238    we may notice , that the used assumption about the phase terms @xmath226 , is rigorously correct in the case without @xmath27 @xcite .
the wave function ( [ wf3 ] ) with the simplification ( [ wfstring ] ) , must satisfy the following constraints due to the @xmath29 term , in order to be an eigenfunction of ( [ hamiltonian2 ] ) :    * for @xmath239 , ( @xmath240 ) @xmath241 * for @xmath242 , ( @xmath243 ) @xmath244    moreover , by symmetry , we have @xmath245 , @xmath246 and @xmath247 .
then , the previous set of equations together with the expression ( [ wfstring ] ) gives @xmath248-ve^{ik_{2}^{0}}\frac{v}{j}\frac{e^{ik_{2}^{l}}}{e^{i(k_{1}^{l}+k_{2}^{l})}+1}=0\\ j[e^{i(k_{2}^{0}+k_{3}^{0})}+1]-ve^{ik_{3}^{0}}\frac{v}{j}\frac{e^{ik_{3}^{l}}}{e^{i(k_{2}^{l}+k_{3}^{l})}+1}=0\\ j[e^{i(k_{1}^{l}+k_{2}^{l})}+1]-ve^{ik_{2}^{l}}\frac{v}{j}\frac{e^{ik_{2}^{c}}}{e^{i(k_{1}^{c}+k_{2}^{c})}+1}=0\\ j[e^{i(k_{2}^{l}+k_{3}^{l})}+1]-ve^{ik_{3}^{l}}\frac{v}{j}\frac{e^{ik_{3}^{c}}}{e^{i(k_{2}^{c}+k_{3}^{c})}+1}=0 \end{array } \right.\]]to which we add the momentum conservation law @xmath249@xmath109 being the momentum associated with the motion of the center of mass of the excitonic complex .
this system of equations is obviously non - soluble in the general case .
hence , we add an additional simplification to the trial wave function .    if @xmath35 , @xmath250 , and therefore @xmath251 ; it is easy to verify that the system ( [ system ] ) becomes equivalent , in that case , to the corresponding equations for the _ xxz _ heisenberg model @xcite as we already pointed out in the previous subsection .
then the system ( [ system ] ) leads to the energy ( [ strings ] ) , for @xmath199 . within the string hypothesis , @xmath252 and @xmath253
are complex conjugate , @xmath254 , a purely real quantity .
enlightened by the exact results for @xmath35 @xcite , we assume for @xmath100 @xmath255where @xmath256 , a real number , is given by the following equation obtained from the energy conservation law @xmath257the left - hand - side is the energy of the @xmath214 component , the right - hand - side , the energy of the three other components .    with this assumption on the momenta @xmath225 ,
the solutions stay very close to the bethe ansatz form ; the binding between the two bordered particles of the complex ( with momenta @xmath258 and @xmath259 ) , ensures the cohesion of the three particle bound states .
the set of equations composed by the two first equations of ( [ system ] ) and equations ( [ q ] ) , ( [ ansatz ] ) and ( [ epsilon ] ) are then solved numerically .
it is important to realize that our approximation is variational .
this important statement could be checked in the following way .
one starts with the trial wave function given by ( [ wf3 ] ) where the components @xmath222 are assumed to be of the form defined by the equations ( [ amplitude ] ) , ( [ wfstring ] ) and ( [ ansatz ] ) together with the restrictions ( [ q ] ) and ( [ epsilon ] ) . the coefficients @xmath221 of ( [ wf3 ] ) , and their complex conjugates , are determined variationally by minimizing the functional @xmath260where @xmath261 is the lagrange parameter for the normalization constraint of the wave function .
we get @xmath262and similar equations for all the other minimizations . at the end of the day
, we obtain nothing else than the usual secular equations and one way to solve them , following a. bethe @xcite , is to impose the constraints given by the sets of equations ( [ constraints1 ] ) and ( [ constraints2 ] ) .
our procedure is then variational and the results found are an upper bound for the problem .
typical results are shown in figure ( [ bindingenergy ] ) .
the binding energy of the biexciton got in the previous section and the one of the triexciton , obtained with our variational calculation , are represented for @xmath96 and @xmath169 as a function of @xmath263 .
the triexciton binding energy is defined with respect to the continuum of one biexciton plus one free exciton which is the lowest in energy .
the following remarks can be made .
( i ) the triexciton binding energy for @xmath35 is the exact one .
( ii ) the binding energy of the triexciton is larger than the one of the biexciton .
( iii ) a critical value of @xmath27 exists also for the triexciton ; it is larger than @xmath143 ( [ alphac ] ) for the biexciton .
it follows that triexcitons may exist without biexcitons ( in the sense that @xmath264 ) .
we may suspect the same behaviour for the binding energy of n - string for any @xmath21 : the binding energy and the critical @xmath27 may be larger for a n - string than for a ( n-1)-string in accordance with the results for the _ xxz _ model .
( iv ) since our result is variational , we get an upper bound for the triexciton binding energy so that our conclusions are qualitatively correct .    to illustrate the behaviour of the binding energy of a
_ n_-string at strong @xmath27 , we study this limit perturbatively . for @xmath265 , the hilbert space of ( [ hamiltonian2 ] )
is naturally separated in two subspaces : the first , lower in energy , with configurations without @xmath176 particles and , the second , higher in energy , with configurations with at least one @xmath176 particle .
the first one is associated with the projector @xmath266 , the second one with the projector @xmath267 ( @xmath268 ) .
the hamiltonian may be written as @xmath269    by a schrieffer - wolf type of canonical transformation @xcite , we want to derive an effective hamiltonian @xmath270 , with @xmath271 , @xmath272 being the effective term to be determined .
the study of the n - excitonic strings will then be performed by diagonalizing the projection of this model into the space without @xmath176 particle , @xmath273 . introducing the unitary transformation @xmath274 ( with @xmath275 ) , and following the method described in detail in @xcite
, we found the simple effective hamiltonian valid in the strong coupling limit where @xmath265 , at second order in @xmath276    @xmath277    we have recovered , once again , the spin-1/2 _ xxz _ hamiltonian ( [ xxz ] ) .
however , the anisotropic parameter is now given by @xmath278 , instead of @xmath279 , and the interpretation of the spin operators is different : @xmath190 is still for a site without any particle but , @xmath191 , is now for a site occupied by a @xmath175 particle instead of a @xmath16 particle .    in the effective model ,
the attractive interaction is reduced by the @xmath27 term and , for the strong limit ( @xmath280 ) where this derivation is valid , one reaches a _
xx_-heisenberg hamiltonian which is well known to be equivalent to a spinless free fermion system @xcite .
hence , in this limit , the excitonic strings do not exist which is consistent with the results of our variational calculation for @xmath281 .
the states with free @xmath175 particles are then preferred .
in this paper , we propose a simplified model to describe excitonic strings in quasi - one dimensional organic compounds including organic charge transfer solids , _
j_- and _ h_-aggregates and conjugated polymers .
for all these compounds , the excitonic states are characterized by small radius .
hence , our model starts with the definition of hard core boson particles extended over one bond which describe , in an effective way , the excitons in one dimensional organic compounds .
there are two possibilities in one dimension to place the electron and the hole , the two elementary constituents of an exciton .
either the electron is on the right or on the left of the hole ; we introduce two kinds of hcb which illustrate these two situations : the _ right_-bosons and the _ left_-bosons .
the proposed model in terms of these effective particles contains four terms : the excitonic energy ( @xmath23 ) , which separate in energy the space with @xmath21 particles from the space with @xmath282 particles ; the kinetic energy of the hcb ( @xmath1 ) ; the local interaction between the two species of bosons ( @xmath27 ) , which illustrates the fact that there exist effective interactions which can exchange the relative position of the hole and the electron of an exciton ; and last , the interaction between two excitons which are attractive between particles of the same species and repulsive between particles of opposite species ( interaction of intensity @xmath29 ) . the kinetic energy and the @xmath27-interaction combined with the repulsive two body interaction act against the attractive two body interaction - they tend to delocalize the excitons on the contrary to the two particle attraction which tends to create bound states .
the cases with two and three particles are mainly studied in this work .
we summarize briefly our results considering the case where @xmath109 , the momentum of the center of mass of the excitonic complexes , is zero .
the case with two particles is solved exactly . for @xmath35 ,
biexcitonic states exist if @xmath283 .
then , the intensity of the transition from the one exciton state decreases dramatically with the increase of the binding energy , to saturate and become independent of the system size .
this behaviour makes the biexcitonic state difficult to observe at the thermodynamic limit . for @xmath100 , a critical value @xmath143 is found , above which biexcitonic states do not exist
otherwise , the behaviour of the oscillator strength remains qualitatively unchanged .    for @xmath35 , the @xmath21 particle case can be solved in the low energy sector .
then , the model is equivalent to the spin one half _ xxz _ heisenberg model which is exactly solvable in one dimension . for @xmath283
, it is well known that @xmath21-excitonic strings exist with a binding energy which increases with @xmath21 @xcite . for @xmath284 ,
the model is non - integrable .
the three particle case is then studied using a variational ansatz , which may be extended to a more general case with @xmath21 particles .
again , a 3-string is found with a binding energy larger than the one for the biexciton .
we also found a critical value for @xmath27 , larger than the @xmath143 valid for the 2-strings .    from numerical studies on small clusters @xcite
, it has been argued that _
n_-strings may exist for any value of _ n _ in organic compounds .
our work gives some confidence to this statement .
indeed , with our hcb hamiltonian , if biexcitons exist , triexcitons will exist with a larger binding energy ; moreover , if @xmath35 , _ n_-strings will also exist with binding energies which increase with _
n _ and , we believe , from our present results , that it will also be the case for @xmath100
.    the next important question concerning the _ n_-string peaks in pa and tpa experiments , is fully clarified for the biexciton ( the most important case ) and confirms some conclusions of previous works @xcite .
the intensity of the transition from the one exciton to the biexciton decreases with the binding energy and , for sufficient binding energy , becomes independent of the system size .
hence , the pa intensity to the biexciton is proportional to @xmath134 instead of @xmath130 for the two free exciton case .
this property comes from the fact that the two excitons are tightly bounded within the complex , rendering the 2-string states difficult to observe for infinite systems ( clean samples are needed ) and more accessible experimentally for small oligomers @xcite .    in the case of conjugated polymers ,
we think some additional studies are needed going in two different directions .
first , in the context of pa experiments , excitations extended over two molecules have been invoked @xcite ; extension of our work to two dimensions or , at least , for two coupled polymers would be suitable .
second , conjugated polymers are disordered systems with several possible sources of disorder which result , in practice , in the vague definition of the so - called conjugation length @xcite .
as we have already mentioned , since the transitions to biexcitons are much less intense than the ones involving two free excitons , we may expect at first analyse - and it was implicitly our point of view all along this manuscript - that a sufficient disorder will render the biexcitons not observable .
however , some relatively recent results about the two interacting particle problem in a disordered medium @xcite show very unexpected behaviours which may have dramatic consequences for our problem , and possibly in a way reverse to our intuition . in the light of the intriguing results of @xcite ,
the effects of disorder on transition moments should be clarified for a correct interpretation of pa and tpa results .
we leave these considerations to further works .
i would like to thank prof .
ovchinnikov for his support and illuminating discussions .
this work would certainly not be developed in the same way without his help .
this work has been supported by the european commission through the tmr network contract erbfnrx - ct96 - 0079 ( qucex ) .
m. ueta , h. kanzaki , h. kobayashi , y. toyozawa and e. hanamura , excitonic processes in solids ( springer , berlin , 1986 ) .
e. hanamura , solid state comm .
* 12 * , 951 ( 1973 ) .
g.m . gale and a. mysyrowicz ,
lett . * 32 * , 727 ( 1974 ) .
s. mazumdar , f. guo , k. meissner , b. fluegel , n. peyghambarian , m. kuwata - gonokami , y. sato , k. ema , r. shimano , t. tokihiro , h. ezaki and e. hanamura , j. chem .
phys * 104 * , 9283 ( 1996 ) .
h. ezaki , t. tokihiro and e. hanamura , phys . rev .
* b 50 * , 10 506 ( 1994 ) .
a. chakrabarti , a. schmidt , v. valencia , b. fluegel , s. mazumdar , n. armstrong and n. peyghambarian , phys . rev . * b 57 * , r4206 , ( 1998 ) .
d. baeriswyl , d.k .
campbell and s. mazumdar , in _ conjugated conducting polymers _ , edited by h. kiess ( springer - verlag , heidelberg , 1992 ) , pp7 - 133 . ,
edited by n.s .
sariciftci ( world scientific publishing , singapore , 1997 ) .
sinclair , d. mcbranch , t.w .
hagler and a.j .
heeger , synth . met .
* 49 - 50 * , 593 ( 1992 ) .
klimov , d.w .
mcbranch , n. barashkov and j. ferraris , phys .
rev . * 58 * , 7654 ( 1998 ) .
townsend , w.s .
fann , s. etemad , g.l .
baker , z.g .
soos and p.c.m .
mcwilliams , chem .
180 * , 485 ( 1991 )
. v.a . shakin and s. abe , phys .
* b 50 * , 4306 ( 1994 ) .
f. guo , m. chandross and s. mazumdar , phys .
lett . * 74 * , 2086 ( 1995 ) .
g. vektaris , j. chem . phys .
* 101 * , 3031 ( 1994 ) .
f.c . spano , chem .
* 234 * , 29 ( 1995 ) .
spano and e. manas , j. chem . phys . * 103 * , 5939 ( 1995 ) .
wu , x. sun and k. nasu , phys . rev . * b52 * , 4849 ( 1995 ) .
f.b . gallagher and f.c .
spano , phys .
b 53 * , 3790 ( 1996 ) .
davydov , theory of molecular excitons , new york , 1971 .
m. chandross , y. shimoi and s. mazumdar , phys .
b 59 * , 4822 ( 1999 ) .
s. pleutin and j.l .
fave , j. phys .
* 10 * , 3941 ( 1998 ) . s. pleutin , in preparation .
simpson , j. am .
soc . * 77 * , 6164 ( 1955 ) .
m.j . rice and y.n .
gartstein , phys .
lett . * 73 * , 2504 ( 1994 ) .
a. horvath , g. weiser , c. lapersonne - meyer , m. schott and s. spagnoli , phys
b 53 * , 13507 ( 1996 ) .
h. bethe , z. physik * 71 * , 205 ( 1931 ) .
b. sutherland , synth . met .
* 84 * , 11 ( 1997 ) . y.a
izyumov and y.n .
skryabin , `` statistical mechanics of magnetically ordered systems '' , consultants bureau , new - york and london ( 1988 ) .
m. takahashi , in `` conformal field theories and integrable models '' , springer - verlag , 1997 ( berlin ) .
m. gaudin , `` la fonction donde de bethe '' , masson ( 1983 ) .
ovchinnikov , jetp lett . * 5 * , 38 ( 1967 ) .
schrieffer and p.a .
wolf , phys .
rev . * 149 * , 491 ( 1966 ) k.a .
chao , j. spalek and a.m. oles , j. phys c : solid state phys * 10 * , l271 ( 1977 ) .
e. lieb , t. schultz and d. mattis , annals of phys .
* 16 * , 407 ( 1961 ) .
hsu , m. yan , t.m .
jedju , l.j .
rothberg and b.r .
hsieh , phys .
* b 49 * , 712 ( 1994 ) .
dorokhov , jetp * 71 * , 360 ( 1990 ) ; d.l .
shepelyansky , phys .
lett . * 73 * , 2607 ( 1994 ) .