the understanding of the relationship between excitation parameters and oscillation regimes is a classical topic concerning bowed string instruments : for instance shelleng @xcite , or guettler guettler , or demoucron et al @xcite proposed 2d diagrams with respect to either bow force and bow position , or bow position and bow velocity .  for reed instruments ,
this kind of diagrams are less numerous : in ref .
@xcite , dalmont et al proposed a diagram with respect to excitation pressure and reed opening , and recently almeida et al almeida proposed a diagram with respect to blowing pressure and lip force , related to the reed opening .    in the minimum model of reed , clarinet - like instruments , three parameters are considered : i ) the mouth pressure , ii ) the reed opening at rest , iii ) the length of the resonator  assumed to be cylindrical . in ref .
dalmont:05 , a supplementary parameter was added : the loss parameter of the resonator ( using the raman model , that considers resonator losses independent of frequency ) .
this allowed explaining the extinction of sound when the mouth pressure becomes very large .
we notice that the agreement of the theoretical results  with experimental results was qualitatively satisfactory ( see ref .
dalmont:07 ) .
the objective of the present paper is to revisit the paper by dalmont et al@xcite : the focus is the search for a diagram of oscillation regimes of reed instruments with respect to two parameters : the reed opening , and the loss parameter .
the choice of these two parameters is justified by the fact that the third parameter , the blowing pressure , is the easiest to modify for the instrumentalist .
the elements of this diagram were rather complete in ref .
@xcite , but phenomena occurring for strong losses , especially at the emergence of the sound , were not investigated .
the use of simplified models for the prediction of the oscillation regimes is classical for musical instruments producing self - sustained oscillations . for the calculation of the instability thresholds ,
linearization was used ( see wilson and beavers @xcite or silva et al @xcite ) , while for ab initio computation , the iterated map scheme was studied ( see mc intyre et al @xcite , maganza et al @xcite , taillard et al @xcite ) .
the interest of the model chosen in the present paper is that analytical formulas are possible , given clear dependence of the parameters ( for other models , numerical calculations could be possible by using similar basic ideas , using for instance continuation methodskarkar or time - domain methods @xcite ) .
the analytical calculations presented hereafter are slightly different than those of ref .
they also are limited to the limit cycles corresponding to the two - state oscillating regime , but are based upon a generalization of the fact that for this regime , when no losses are present , the flow rate is a constant .
this regime is the most similar to what musicians consider as a normal " sound . in particular the method allows studying the character of the emergence and extinction bifurcations of this regime , which are important properties related to the possibility to play pianissimo or not , and more generally to the ease of playing . in ref .
@xcite it was shown that the extinction bifurcation can be direct or inverse ; here it is shown that this is true also for the emergence bifurcation .    in section [ model ] , the basic model of ref .
@xcite is reminded ( see also @xcite ) , and a treatment of the problem based upon a unique quantity , the pressure difference @xmath0 is presented in section [ recu ] .
this leads to a simple graphical analysis of the two - state regime , explained in section [ exis ] , yielding a proof that it can not exist with reverse flow , and an easy calculation method . then some blowing pressure thresholds ( stability , existence ,
... ) are calculated with respect to the parameters of interest ( loss , reed opening ) .    in section [ ct ] the thresholds related to the instability of the regimes
are calculated . then by making two thresholds equal ,
the mouth pressure can be eliminated and limits of existence and stability of the static and two - state regimes are found in section [ dia ] : this allows drawing the diagram sought .
finally a discussion is proposed concerning the existence of oscillating regimes ( section [ disc ] ) , with an attempt to consider more realistic models and a discussion about musical consequences .
we briefly remind the basic elements of the model , the non - linear characteristic of the exciter , and the origin of the iteration method , thanks to a simplified treatment of the resonator .      in a quasi static regime , the flow @xmath1 entering the resonant cavity
is modeled with the help of an approximation of the bernoulli equation , as discussed e.g. in @xcite .
comparison with experiment can be found in ref .
we note @xmath2 the acoustic pressure inside the mouthpiece , assumed to be equal to the one at the output of the reed channel , @xmath3 the pressure inside the mouth of the player . for small values of the difference : @xmath4the reed remains close to its equilibrium position , and the flow @xmath1 is proportional to @xmath5 ; for larger values of this difference , the reed moves and , when the difference reaches the closure pressure @xmath6 , it completely blocks the flow .
these two effects are included by assuming that if @xmath7 the flow @xmath8 is proportional to @xmath9 $ ] , and if @xmath10 , the flow vanishes.introducing the dimensionless quantities : @xmath11where @xmath12 is the characteristic acoustic impedance of the cylindrical resonator , having the cross section @xmath13 ( @xmath14 is the density of air , @xmath15 the velocity of sound ) , we obtain : @xmath16with@xmath17the parameter @xmath18 characterizes the intensity of the flow and is defined as : @xmath19where @xmath20 is the opening cross section of the reed channel at rest .
@xmath21 is inversely proportional to square root of the reed stiffness , contained in @xmath6 . in real clarinet - like instruments , typical values of the parameters
are @xmath22 $ ] ; @xmath23 $ ] ; values @xmath24 will not be considered here , since they correspond to multi - valued functions to be solved ( see ref .
@xcite ) , and this case does not seem very realistic in practice for clarinet - like instruments .
the function @xmath25 is obviously non - analytic ; it is made of three separate analytic pieces , with a singular point at @xmath26 , and its derivative is discontinuous at @xmath27 .
the resonator of length @xmath28 is assumed to be cylindrical , with zero terminal impedance . using the dalembert decomposition , a change in variables at the entry of the resonator
can be chosen as:@xmath29with the following relationship between incoming wave @xmath30 and outcoming wave @xmath31:@xmath32where @xmath33 is the loss parameter , assumed to be independent of frequency .
this is a strong assumption , necessary to obtain square signals : for certain initial conditions , all quantities remain constant in the time interval @xmath34 .
the approximation is rough for certain characteristics of the signal , such as the spectrum , but it is useful for the study of the existence , stability and amplitude of the produced sound . with this assumption ,
the resonator is characterized by a unique ( recurrence ) relation:@xmath35where @xmath36 is the time unit . as discussed further in section [ rs ] , losses can occur either at the extremity ( radiation ) or during propagation ( e.g. in the boundary layers ) : in the latter case , if @xmath37 is the attenuation constant per unit length , @xmath38 and the dimensionless input impedance at zero frequency is @xmath39while the input impedance at the operating frequency is @xmath40 . in what follows ,
the losses are characterized by the parameter @xmath41 , varying between @xmath42 ( no losses ) @xmath43and @xmath44 ( very strong losses , no wave reflection ) .
this parameter , the reed opening @xmath18 and the mouth pressure @xmath45 are the three parameters of the problem .
several combination parameters will be useful:@xmath46@xmath47 is proportional to the input impedance at zero frequency , while @xmath48 is proportional to the input admittance at the operating frequency .
two pairs of parameters can be used : either @xmath49 or @xmath50 notice that because @xmath18 and @xmath41 are smaller than unity , @xmath51 and @xmath52 , and @xmath53 .
other parameters will be useful and @xmath18 are with dimension , except in the appendix , and are denoted @xmath54 and @xmath55 , respectively .
@xmath47 , @xmath56 and @xmath57 are defined in the same way than in the present paper.]:@xmath58
using eqs . ( [ e1 ] ) , ( [ e3a ] ) and ( [ e49 ] ) , the recurrence relation ( [ e4 ] ) can be rewritten with respect to the quantities @xmath59 and @xmath60 .
the result is:@xmath61because the flow rate @xmath62 is a function of the pressure difference @xmath0 this relation is a recurrence for the quantity @xmath63 , equivalent to the recurrence used in ref .
@xcite for the quantity @xmath64:@xmath65 \text { \ \ with } \label{pa2 } \\
h(x ) & = & x+\zeta f(x)\text { , \ \ } k(x)=\gamma ( 1+\lambda ) -\lambda ( x-\zeta f(x ) ) .
\notag\end{aligned}\]]the inverse of function @xmath66 can be found in ref .
@xcite ( appendix a ) .
( [ pa1 ] ) is interesting in particular for the calculation of the limit cycles . for the static regime ,
@xmath67 is a constant , then@xmath68it is possible to calculate @xmath67 from the value of @xmath45 , or vice - versa .
concerning the stability , if the iteration function is denoted @xmath69 , $ ] the classical stability condition is @xmath70 . because @xmath71 , $ ] @xmath72 , and the condition can be written as is equivalent to @xmath73:@xmath74 ^{2}<1\text { \ \ or \ } \left [ f^{\prime } ( \delta p)+\beta _ { 1}\right ] \left [ 1+\beta f^{\prime } ( \delta p)\right ] > 0\text{\ } \label{e10}\ ] ]      for the two - state regime , because @xmath75 , the following expression is found by eliminating @xmath45 from the equation ( [ pa1 ] ) written for the pairs ( @xmath76 , @xmath77 ) and ( @xmath77 , @xmath78):@xmath79an important property of the two - state regime is the square shape of the signal , which can be decomposed into the sum of a mean value @xmath80 ( zero frequency component ) and an acoustic component @xmath81 ( sum of the odd harmonics of the operating frequency ) , with zero mean value:@xmath82notice that @xmath83 , and that a similar equation can be written for the flow rate .
( [ e889 ] ) generalizes the result obtained when losses are ignored ( @xmath84 , i.e. the constant flow rate , and it is nothing else than the input impedance relation for the acoustic component:@xmath85    it is possible to calculate the values of the two states without knowledge of @xmath45 , starting e.g. with the value of @xmath86 : @xmath87 and @xmath45 are successively deduced from eqs . ( [ e889 ] ) and eq .
( [ pa1 ] ) . adding the two equations ( [ pa1 ] ) for the pairs ( @xmath88 , @xmath77 ) and ( @xmath77 , @xmath78 ) , it is obtained : @xmath89 \text { , or }
\label{e921a } \\
\gamma _ { ij } & = & \frac{1}{2}(\delta p_{i}+\delta p_{j})(1-\beta \beta _ { 1})+\beta h_{1}\text { , }   \label{e921}\end{aligned}\]]with @xmath90 and @xmath91 , and @xmath92 . similarly for regimes with more than two states , it could be possible to start the calculation from a given state , and to deduce the other states , the prior knowledge of @xmath45 being unnecessary .
( [ e921a ] ) is the input impedance relation for the mean value component : @xmath93    for the two - state regime , the stability condition is @xmath94 , and after some algebra , the result of ref .
@xcite is found:@xmath95
the previous results can be applied whatever the shape of the function @xmath96 the present section investigates the existence of the static and two - state regimes for the particular shape of the function given by eqs .
( [ e2 ] ) and ( [ e2a ] ) .      for negative @xmath67 ,
the function @xmath25 is negative too , therefore the static regime does not exist for negative flow ( and positive excitation pressure @xmath45 ) . for @xmath97 , the static regime exists for @xmath98 .
the pressure @xmath99 vanishes : this is obvious because the reed closes the input of the resonator .    otherwise , for a non - beating reed , the study of the function in the right - hand side of eq .
( [ e8 ] ) shows that it increases from @xmath42 to @xmath44 when @xmath67 increases from @xmath42 to @xmath100 therefore a unique solution \right ] ^{2 } \\
\text{with } \kappa & = & \left [ -2 + 9(-1 + 3\gamma ) \beta ^{2}\right ] /\left
[ 2(1 + 3\beta ^{2})^{3/2}\right ] .\end{aligned}\ ] ] ] exists  for @xmath101 .
for the two - state regime , two values of @xmath67 , with the same value of the function @xmath102 are sought .
they do not depend on the impedance at zero frequency , i.e. on the value of @xmath103 the study of @xmath102 leads to the following results ( see fig .
[ fighx ] ) : for negative @xmath104 , the derivative @xmath105 is always positive , while for positive @xmath104 , it is positive up to @xmath106 , where @xmath107 is the positive root of the following equation:@xmath108 ^{2}.   \label{e912}\end{aligned}\]]for @xmath109 ( beating reed ) , the value of the derivative is @xmath56 , and is positive .
two cases need to be distinguished : if @xmath110 , i.e. if @xmath111 , the derivative is always positive and it is impossible to find two values of @xmath104 with the same value @xmath112 on the contrary , for @xmath113the function decreases from @xmath114 to @xmath56 when @xmath104 increases from @xmath115 to @xmath44 , then re - increases . solutions of eq .
( [ e889 ] ) can be found in this case only , and the corresponding value of the function is necessarily larger than @xmath56 .
a consequence is that eq .
( [ e889 ] ) has no solution for @xmath116 , and this is true in particular for negative @xmath102
.  therefore no two - state regime can be found with negative flow .
this conclusion is compatible with the general result obtained for all possible regimes in ref .
@xcite . because of this result , the paper is now focused on positive pressure differences @xmath117    for @xmath118 $ ] , three values of @xmath104 , defined as @xmath119 lead to the same value of the function ( see fig . [ fighx ] ) , and three two - state regimes are possible , which can be either non - beating ( for the pair ( @xmath120 , @xmath121 ) ) or beating ( for the pairs ( @xmath120 , @xmath122 ) and ( @xmath121 , @xmath122 ) ) .
the intervals of the three solutions are as follows : @xmath123 $ ] ; @xmath124 $ ] ; @xmath125 .    [ cols="^,^ " , ]      the radiation losses are small at low frequencies , therefore it is very simple and classical to deduce them by perturbation from the output flow rate , considering a monopole radiation .
two cases have to be distinguished : losses occur at the output , or losses occur during propagation into the tube .
both cases give the same input impedance , but not the same transfer functions between input and output .
the latter case is more realistic , and is considered here . if the output impedance is @xmath126 and losses due to boundary layers only , the output acoustic flow rate is given by the standard transmission lines relationships : @xmath127 , where @xmath128 ( see section [ besr ] ) .
therefore the amplitude relationship is the following:@xmath129the maximum output acoustic flow rate is obtained for the saturation threshold ; at the saturation threshold , because the reed is beating , @xmath130(see eq . ( [ e18h ] ) ) . for @xmath131
the saturation threshold is the beating threshold ( see section [ eb ] ) , and @xmath132 , therefore:@xmath133for a given value of @xmath18 , this is monotonously decreasing function of @xmath134 ; for a given value of @xmath41 , this is a increasing function of @xmath135
the stability conditions ( [ e10 ] ) and ( [ e301 ] ) are calculated using the expression of the derivative @xmath136:@xmath137      the condition ( [ e10 ] ) generalizes the condition @xmath138 , as discussed in ref .
@xcite for the lossless case ( see page 349 ) .
for the ( static ) beating reed case , which exists for @xmath139 ( see previous section ) , @xmath140 , thus the static regime is always stable . for the non - beating reed case ,
the first factor of inequality ( [ e10 ] ) is equal to @xmath141 , thus it vanishes when @xmath63 satisfies eq .
( [ e914 ] ) , i.e. @xmath142 . for this value , @xmath115 , the second factor of inequality ( [ e10 ] )
is @xmath143 , which is positive ) vanishes : this leads to the solution of either @xmath144 when @xmath97 , or @xmath145 , with @xmath146 .
these two equations have no solution . ] . thus ( [ e912 ] ) together with eqs .
( [ e8 ] ) gives the threshold @xmath147 , given by eq .
( [ e923 ] ) .
therefore the condition @xmath148 gives both the limit of existence of the two - state regime ( see eq .
( [ e923 ] ) ) , and the instability threshold of the static regime .
nevertheless the nature ( direct or inverse ) of the bifurcation between the two regimes is not yet solved by this result was assumed to be very small in practice , and eq .
( [ e923 ] ) was simplified in @xmath149 , but the complete eq .
( [ e923 ] ) was given for the threshold of existence for the two - state regime . ] .
the term @xmath150 is the static pressure @xmath151  in the mouthpiece ; the presence of the parameter @xmath47 indicates that the threshold depends of the impedance at zero frequency .        for the beating regime @xmath152 , with @xmath153
thus , using condition ( [ e301 ] ) , @xmath154  for @xmath155 or @xmath156 .
because @xmath157 , and because @xmath158 , the second inequality ( [ e301 ] ) is never valid
. therefore the stability is defined by the condition @xmath159    when @xmath160 , this condition is always satisfied , and the beating two - state regime is stable ( no overcritical threshold @xmath161 exists , as noticed in section [ eb ] ) .    when @xmath162 , eq .
( [ e18f ] ) is used , yielding @xmath163 .
the function @xmath164 being always increasing for positive @xmath165 , the inequality @xmath166 holds if @xmath167 , or @xmath168 this distinguishes in the ( @xmath169 ) plane the two branches separated by the overcritical threshold : the upper one is stable , while the lower one is unstable : because @xmath170 , the solution @xmath171 at the overcritical threshold is larger than @xmath115 , thus it is always larger than @xmath172 instability occurs for the pair ( @xmath173 @xmath122 ) .
this is confirmed by figs .
[ pfigiii ] and [ pfigi ] . ] .
the unstable branch is the branch joining the static regime , because the two regimes can not be stable for the same value of the parameter @xmath45 when they converge to the same point .
this can be explained with mathematical arguments , based upon either a perturbation method ( see ref .
@xcite ) or the topological degree ( see ref .
this can be also studied by using inequalities ( [ e301 ] ) , as done in ref .
@xcite .
it can be noticed that because @xmath161 is larger than unity , the beating , two - state regime is always stable for @xmath174 , whatever the value of @xmath175      \i ) for the non - beating regime ( @xmath155 , @xmath176 ) , an expression of the instability threshold was given in ref . @xcite and it was explained that the threshold is given by @xmath177 .
some errors were done in the derivation , and they are corrected in appendix a of the present paper . when the excitation pressure is larger than this threshold , here denoted @xmath178 , period doubling can occur , then a complex bifurcation scenario ( see ref .
@xcite ) .
\ii ) otherwise , it can be checked that the second condition @xmath179 leads to the subcritical emergence threshold @xmath180 : it separates two branches in the ( @xmath181 plane , the upper one being stable while the lower one is unstable .
this is similar to what happens for the overcritical threshold .
when it exists , the emergence bifurcation is as follows : when playing crescendo , the oscillation starts for @xmath182 , while playing decrescendo , the oscillation stops for @xmath183    when this subcritical threshold can exist ? from eq .
( [ a ] ) and ( [ cz ] ) , it is found that @xmath184 \label{cq}\]]where @xmath185 and @xmath186 , with @xmath187    two cases are possible :    * @xmath188 : the emergence bifurcation is direct : stable solutions exist for @xmath189 . *
@xmath190 : the emergence bifurcation is inverse , and stable solutions exist for @xmath191    when @xmath192 continues to decrease below @xmath193 , the bifurcation remains inverse , but the subcritical threshold @xmath180 becomes the beating threshold @xmath194 .
this happens when @xmath195 , or @xmath196 ( notice that the inequality @xmath197 always holds ) .
the discussion is extended in section [ ol ] .
from the different expressions of the thresholds , it is possible to deduce the limits separating different domains in the plane @xmath199 , as shown in fig .
[ pertes4 ] . above the diagonal ( region 0 ) , no two - state regime can exist .
other regions of the plane are defined by the nature of the emergence and extinction bifurcations : in regions i and iii , the extinction bifurcation is inverse , while in regions ii  and iv , it is direct .
what is new in this paper is the separation between regions i and ii , with direct emergence bifurcation , and regions iii and iv , with inverse emergence bifurcation .
finally , in region v , the two - state regime can be unstable , and can be replaced by more complicated regimes , with period doubling , chaos , intermittences , etc ... ( see ref .
@xcite ) . figures [ pfigi ] , and [ pfigii ] to [ pfigiv ] show examples of bifurcation schemes for the regions i to iv .  the calculation is made using the analytical formulas given hereafter , and
have been verified using the iterated map algorithm ( ref .
@xcite ) : obviously only stable solutions are found with the latter method . the limits corresponding to the beating threshold ( between regions ii and iv , and between i and iii ) are discussed in section [ ol ] .
the four regions i to iv correspond to the different cases presented for the examples shown in fig .
[ pfigi ] , and [ pfigii ] to [ pfigiv ] .
the chosen values for the examples are indicated by a circle . in region v period doubling or other regimes can exist , in region 0 no sound is possible .
the numbers refer to the equation number . above line ( [ e22 ] ) ,
no sound is possible .
curve ( [ e23 ] ) distinguishes the extinction bifurcation ( direct above , inverse below ) .
curve ( [ e24 ] ) distinguishes the emergence bifurcation ( direct above , inverse below ) . below curve ( [ e237 ] ) , complicated regimes can appear by destabilization of the two - state regime .
the two curves close to the curve ( [ e24 ] ) correspond to equations ( [ e25 ] ) , and ( [ e236 ] ) , from to the top to the bottom ( see section [ ol]).,width=377 ]    , @xmath200 ( point ii of fig . [ pertes4 ] ) .
see figure caption of fig .
[ pfigi ] left . both emergence @xmath201 and extinction bifurcations
@xmath202 are direct .
the point @xmath203 indicates the beating threshold @xmath204 of the two - state regime . because @xmath205 , saturation occurs at the beating threshold @xmath206    , @xmath207 ( point iii of fig .
[ pertes4 ] .
see figure caption of fig .
[ pfigi ] left .
both emergence and extinction bifurcations are inverse .
the point @xmath208 indicates the overcritical ( extinction ) threshold @xmath209 of the two - state regime , and the point @xmath202 the closure threshold @xmath210 of the static regime . ]    , @xmath211 ( point iv of fig .
[ pertes4 ] ) .
see figure caption of fig .
[ pfigi ] left . the emergence bifurcation is inverse while the extinction one is direct . ]
* when the instability threshold @xmath212(eq .
( [ e923 ] ) ) of the static regime reaches the closing threshold @xmath213 , the static regime becomes stable whatever the values of all parameters , and no sound can be expected .
this happens has two solutions : @xmath214 , and @xmath215 .
because @xmath47 is necessarily less than unity , the latter solution is larger than unity . ] for @xmath111 , and this confirms the result explained in section [ nsol ] that no two - state regime can exist ( this discussion is extended in the next section ) .
the condition for the existence of sound can be also written as : @xmath216 * when the overcritical threshold @xmath161 ( eq .
( [ e18q ] ) ) reaches the closure threshold of the static regime @xmath213 ( eq . ( e19 ) ) , the extinction bifurcation becomes direct instead of inverse , as explained in the previous section . for @xmath217
the bifurcation is inverse : this is probably the most frequent case for real clarinets and clarinettists , and corresponds to:@xmath218between the two limits ( [ e22 ] ) and ( [ e23 ] ) , the extinction bifurcation is direct ( regions i and iv ) .
* when the subcritical threshold @xmath180 ( eq .
( [ cy ] ) ) reaches the instability threshold of the static regime @xmath147 ( eq . ( e923 ) , the emergence bifurcation becomes direct . from eq .
( [ cq ] ) , this happens when @xmath219 .
the bifurcation is direct in regions i and ii , with the following condition : @xmath220      finally , when the instability threshold @xmath178 of the oscillating regime given by the condition @xmath177 ( see section sb ) reaches the beating threshold @xmath194 , the limit was given in ref .
@xcite ( with a small error ) .
the formula can be obtained from eq .
( [ e301 ] ) , for @xmath221 , and @xmath222 , @xmath223 the limit @xmath224 is given by the following equation:@xmath225this leads to a second order equation in @xmath226:@xmath227a particular value is @xmath228 , @xmath229 , @xmath230 , @xmath231 , @xmath232 . for a given loss coefficient @xmath41 @xmath233 , the two - state regime is always stable when @xmath234 . for sake of simplicity ,
we remark that an excellent approximation of the limit , better than 1% , is based upon the fact that for small losses , the coefficient @xmath235 is small , therefore @xmath236:@xmath237      another limit is reached when the subcritical threshold @xmath180 becomes the beating - reed threshold @xmath194 , thus the stable two - state regime is always beating
. however the unstable branch is non - beating ( see figs . [ pfigiii ] and [ pfigiv ] ) .
( [ e ] ) , it is found to be given by @xmath238 , i.e.@xmath239    a last limit is reached when the instability threshold of the static regime @xmath147 ( [ e923 ] ) reaches the beating threshold of the two - state regime @xmath194 , i.e. when @xmath240 in eq .
( [ a]):@xmath241for both limits given by eqs .
( [ e25 ] ) and ( [ e236 ] ) , the corresponding curves are therefore within regions iii and iv of figure pertes4 .
they are very close to the limit given by eq .
( [ e24 ] ) , corresponding to the change in nature of the emergence bifurcation .
the three limits are even equal for @xmath242 , @xmath243 , @xmath244 for a given @xmath18 , when losses are small ( @xmath41 small ) , the emergence bifurcation is direct . then , when @xmath47 reaches the limit ( e24 ) , the bifurcation becomes inverse , with a non - beating reed .
then , when @xmath47 reaches the limit given by eq .
( [ e25 ] ) , the bifurcation remains inverse , but the reed becomes always beating in the two - state regime .
[ figbdzoom ] shows details of the bifurcation scheme ( @xmath245 near the subcritical threshold for different cases between the non - beating case @xmath246 and the beating case @xmath247 ( for the latter case , the beating threshold is the subcritical threshold , as for the points iii and iv ) .
the corresponding values of the parameter ( @xmath248 ) are extremely close together .
* left up : * case @xmath249 , @xmath250 * right up : * case @xmath251 , @xmath252 ; * left down : * case @xmath253 , @xmath254 * right down : * case @xmath255 , @xmath256 , title="fig:",width=188 ]   * left up : * case @xmath249 , @xmath250 * right up : * case @xmath251 , @xmath252 ; * left down : * case @xmath253 , @xmath254 * right down : * case @xmath255 , @xmath256 , title="fig:",width=188 ]   * left up : * case @xmath249 , @xmath250 * right up : * case @xmath251 , @xmath252 ; * left down : * case @xmath253 , @xmath254 * right down : * case @xmath255 , @xmath256 , title="fig:",width=188 ]   * left up : * case @xmath249 , @xmath250 * right up : * case @xmath251 , @xmath252 ; * left down : * case @xmath253 , @xmath254 * right down : * case @xmath255 , @xmath256 , title="fig:",width=188 ]
when losses tend to infinity ( @xmath41 tends to unity ) , no sound is possible whatever the value of @xmath18 : this is in accordance with the intuition that if no reflection exists at the input of the resonator , no self - sustained oscillation can happen .
nevertheless , we do not prove that other types of regimes can not exist , such as four - state , eight - state , ... in this section , this issue  is discussed together with the influence of the choice of the model . moreover ,
some musical consequences are discussed .
when the static regime is always stable ( @xmath111 ) , it has been proved that the two - state regime can not exist .
it is probable that other types of regimes do not exist , but the general proof is difficult . using the calculation of the successive iterates functions for different values of the initial condition ( see ref .
@xcite ) , we verified that when @xmath257 , the successive iterates converge to the unique point that is the limit cycle of the static regime , thus no other regime can exist .
this can be done for every set of parameter values , @xmath45 , @xmath18 , @xmath258 : if the convergence is always to a unique state , then it is sure that no other regime than the static one can exist .
obviously , this verification is not possible in practice , but the verification has been done for some set of parameter values .    we conclude that no oscillating regime exists for @xmath111 , even if the rigorous sentence should be : no oscillating regime exists by destabilization of the static regime .
a similar discussion could be done for the destabilization of the two - state regime into more complicated regimes ( see section [ lit ] ) .
the method used in the present paper can also be used for any shape of the nonlinear characteristic , at least numerically ( the condition being that a nonlinear characteristic is static ) .
all equations of sections [ model ] and [ recu ] remain valid when modifying the nonlinear function @xmath259 in particular if a smooth beating transition is chosen with no singularity , fig .
[ fighx ] shows that the condition ( [ e913 ] ) , @xmath260 , can be generalized into the following condition : the function @xmath261 goes through both a maximum and a minimum .
the raman model is interesting because all quantities can be determined analytically , but it is not very realistic .
it is based upon two important assumptions : losses do not depend on frequency , and the reed has no dynamics .
a generalization of the present study is out of the scope of the present paper , but it is interesting to note that the condition @xmath262 can be easily generalized when these assumptions are not done , as explained hereafter .
when losses depend on frequency , it is possible to use the characteristic equation obtained by linearizing the nonlinear equation around the pressure of the static regime @xmath151 , and writing the approximation of the first harmonic:@xmath263where @xmath264 is the admittance of the fundamental frequency . the characteristic equation is written as:@xmath265as it is well known ( see ref .
@xcite ) , this gives the condition @xmath266 , thus the playing frequency @xmath267 at the threshold can be deduced .
moreover , if at this frequency , @xmath268 , the pressure threshold is given by : @xmath269because @xmath270 therefore @xmath271 satisfies eq .
( [ e914 ] ) : @xmath272as expected in section [ is ] . as a consequence , when the losses depend on frequency , the value of the threshold is the same as for the raman model , and the limit of existence of the two - state regime @xmath242 is unchanged .
nevertheless the hypothesis has been done that the small oscillations are sinusoidal ( on this subject , see refs .
@xcite , grand or @xcite ) , and this is not true for inverse bifurcation .
when several harmonics interact , the problem becomes much more intricate , especially because of resonance inharmonicity .
moreover taking into account the frequency dependence of the losses leads to a distinction between the threshold of the fundamental regime and the overblown " regimes ( see ref .
@xcite ) : this distinction does not exist with the raman model , which allows a distinction based upon the initial conditions only .    as a summary
, it can be said that if we suppose that the impedance peak of the operating frequency is higher than the other ones , the emergence bifurcation is direct and the limit @xmath273 is valid .
this is true in particular for the first register of a clarinet and a part of the second register .
the large increase of radiation losses at higher frequency due to the open toneholes lattice ( see e.g. ref .
@xcite ) does not affect the highest peak , therefore the main result of the present paper can be extrapolated to a large number of notes of a real clarinet .
when the reed dynamics is taken into account as that of a 1 dof oscillator , the following characteristic equation has been obtained by silva et al silva:@xmath274   \label{e80}\]]where @xmath275 , @xmath276 is the reed - resonance angular frequency , and @xmath277 its quality factor . the threshold pressure and frequency can be deduced from this equation , and were studied in this paper ; here we are interested in the limit for which the static regime becomes always stable , i.e. when @xmath278 if the input impedance of the resonator is considered around a resonance frequency @xmath279 , it is possible to write:@xmath280\]]where @xmath281 is the quality factor of the resonance . for @xmath282 the real part of eq .
( [ e80 ] ) leads to the following result:@xmath283for a lossless reed , @xmath284 : the limit of the losses in the tube allowed for having a sound is increased by the reed dynamics , which favors the sound production .
but the effect of the reed losses is to decrease the limit .
it can be concluded that reed losses and resonator losses act in the same sense concerning the range of parameter allowing sound production ( this conclusion is valid for the direct - bifurcation case ) .
the previous results can be useful in order to understand and teach important aspects of the sound control by the instrumentalist .
such an objective knowledge should largely increase the pedagogical efficiency .
otherwise the approach of the problems remains more subjective and the explanations can be lengthy and less clear .
one of the most useful aspects is about _ pianissimo _ playing .
the bifurcation diagrams show that the players have two possibilities : near the emergence and near the extinction .
the first possibility is used for playing _ dolce _ , with a quasi monochromatic sound , but the sound is noisy and can not be sustained for a long time , due to high value of the airflow @xmath285 ( see fig .
[ pfigi ] ) .
the dynamic is not easy to control because of the steepness of the bifurcation diagram near @xmath147 .
the second possibility , near @xmath161 , conducts to a clean pianissimo , with a sound richer in high harmonics .
this can however only be achieved by crossing the curve ( [ e23 ] ) in fig.pertes4 , in order to reach the region ii where the extinction bifurcation is direct .
this property is usually unknown by the players ( the ability of playing such a magical `` pianissimo is often attributed exclusively to the quality '' of the reed ) .
the beginners reduce the reed opening ( @xmath286 by biting `` the reed and this causes unwanted effects : the pitch rises considerably and such a bending stress can cause a plastic ( irreversible ) deformation of the reed . the skilled player can reach region ii by increasing the damping due to the lip and use a high blowing pressure near the extinction threshold , playing in the reverse way ( decreasing the mouth pressure for playing louder , see fig .
[ pfigii ] ) .
the lip comes very near to the tip of the reed , with a moderate lip pressure .
this effect is probably similar to an increase of @xmath48 , resulting in a displacement of the playing parameters on the @xmath287 plane exactly in the wanted direction , increasing the value of @xmath41 and decreasing the value of the parameter @xmath18 proportionally to @xmath288 .
the decrease in @xmath18 is probably due to an increase of the reed stiffness . decreasing @xmath18 without biting '' too much
could also be achieved by modifying the hydrodynamics of the airflow entering the channel , in order to increase the _ vena contracta _ , but to our knowledge no experimental evidence shows that the player can indeed modify significantly the _ vena contracta_. conversely , it seems that the parameter @xmath47 can not be significantly controlled by the player ( in another way than by modifying the length of the air column ) . in real life
, this parameter of static airflow resistance may not be determined only by the length and the diameter of the bore but certainly also by hydrodynamic effects near the channel , due to the viscosity of the air .
this acts in a similar but possibly stronger way than the static resistance of the bore .
practical tests show that the effects predicted in the zones iii and iv are indeed observed in some pathological situations , despite the fact that our theoretical model would require many meters of a tube of small diameter to obtain such high values of @xmath47 .    to include the musician mouth in the model is obviously rather complicated , even if at low frequencies , the effect of the vocal tract is not important
. therefore the previous analysis requires some conjectures . besides the problems of bifurcation , the analysis of the raman model permits establishing some facts useful for the musician :    * the calculation of the mean flow shows that the most economical blowing pressure is near the beating threshold in region i ( corresponding to normal playing )
this explains that skilled players can sustain the sound significantly longer than beginners . *
the transients are much faster if @xmath18 is large ( see ref .
kergomard:95 ) ; weak reeds help doing this , as well as using a moderate lip pressure .
this simplifies the staccato learning . * the effects of leaks in the air column ( misplacement of a finger , defective pads ) increase the value of @xmath56 , so that regions 0 , and probably iii or iv can be possibly visited ( see section [ fdl ] ) .
almost any control can be destroyed over the dynamics ( or at least rendering the dynamic control more difficult ) , despite of the attempts of the clarinettist to supply more energy to the instrument by opening the embouchure , increasing @xmath18 .
the present paper is focused on limit cycles corresponding to two - state regimes , and is a complement to the paper @xcite , which was focused on transients .
thanks to a formulation focused on the pressure difference between mouth and mouthpiece , the effect of the nonlinear function on the production of the two - state regime can be analyzed , and especially the role of the losses .
the map shown in fig [ pertes4 ] can certainly be improved by using more complex models , but we think that some results are robust .
when the reed opening at rest is very small or when the reed stiffness is very large ( i.e. when the dimensionless parameter @xmath18 is very small ) , losses can be too large and the sound production becomes impossible .
a complement to this conclusion is the following : for @xmath18 larger than @xmath289 , when losses increase , the emergence bifurcation becomes inverse before the sound disappears , and the instrument becomes more difficult to play . for @xmath18
smaller than this value , when losses increase , there is a direct passage from the direct emergence bifurcation to the absence of sound .
this work was supported by the french national agency anr within the sdns - aimv and cagima projects .
we thank also the high school arc - engineering in neuchtel .
finally we wish to thank philippe guillemain and christophe vergez for fruitful discussions .
the instability threshold @xmath178 of the oscillating regime is given by the condition @xmath177 ( see ineqs .
( [ e301 ] ) , it is the limit of unstable solutions toward period - doubling of the two - state regime ) .
this leads to the following equation , if @xmath290 and @xmath291 are defined by eq .
( [ aaa ] ) :      together with eq .
( [ azr ] ) leads to a fourth - order equation in @xmath290 ; from the solution @xmath290 the value the threshold of instability @xmath293 is deduced by using eq .
( [ a ] ) .
another method is to start from certain values of the parameter @xmath56 and of the solution @xmath290 , then to deduce @xmath291 using eq .
( [ azr ] ) , then @xmath47 , which is solution of a second order equation    in ref .
@xcite , eq .
( a23 ) was correct , but a factor 4 was missing in eq .
( a24 ) , the correct equation being the present eq .
( [ e21 ] ) .
( a25 ) needs to be corrected by introducing a factor 4 on the right - hand side , and the factor @xmath294 needs to be replaced by @xmath295 in eq .
( a28 ) and similarly @xmath296 needs to be replaced by @xmath297 in eq .
( a30 ) .    concerning the limit @xmath298 , eqs .
( a32 ) to ( a34 ) of ref .
@xcite are corrected in section [ lit ] of the present paper .
the last equation gives the coefficient @xmath56 as a series expansion of the limit @xmath224:@xmath299   \label{216}\]]this expression corrects eq .
( a34 ) of the previous paper ( the correction is small , because the order 6 in @xmath18 only is concerned ) , but this approximation is much less accurate than the present eq .
( [ e238 ] ) .
a. almeida , j. lemare , m. sheahan , j. judge , r. auvray , kim son dang , s. john , j. geoffroy , j. katupitiya , p. santus , a. skougarevsky , j. smith and j. wolfe , clarinet parameter cartography : automatic mapping of the sound produced as a function of blowing pressure and reed force proceedings of the international symposium on music acoustics , 2010 , sydney , australia .
kergomard , j. , `` elementary considerations on reed - instrument oscillations '' . in _ mechanics of musical instruments _ , vol .
* 335 * ( a.  hirschberg/ j. kergomard/  g. weinreich , eds),of _ cism courses and lectures _ , pages 229290 .
springer - verlag , wien , 1995 .
hirschberg , a.  , van  de laar , r.  w.  a. , marrou - maurires , j.  p. , wijnands , a.  p.  j. , dane , h.  j. , kruijswijk , s.  g. and houtsma , a.  j.  m.  a quasi - stationary model of air flow in the reed channel of single - reed woodwind instruments  .
_ acustica _ , 70 , 1990 , 146154 .
dalmont , j .-
p . ,  gilbert , j.  and ollivier , s.  .
nonlinear characteristics of single - reed instruments : quasistatic volume flow and reed opening measurements , j. acoust .
_ , _ 114 , 2003 , 22532262 .


studies of cp - violation in @xmath0 decays must involve final state interaction ( fsi ) effects .
unfortunately , a reliable estimate of such effects is very hard to achieve . in the analyses of @xmath9
decays ( @xmath10 - pseudoscalar meson ) only some intermediate states , believed to provide nonnegligible contributions , are usually taken into account .
many authors restrict their studies to elastic rescatterring @xmath11 only . in regge language
this is described in terms of a pomeron exchange .
although in @xmath12 quasi - elastic rescatterring at @xmath13 contributions from other nonleading regge exchanges are much smaller , they are not completely negligible and have been included in various analyses .    the main problem , however , is posed by the sequence @xmath14 involving _ inelastic _ rescatterring processes @xmath15 .
arguments have been given that it is these inelastic processes that actually constitute the main source of soft fsi phases @xcite .
it has been also pointed out @xcite that nonzero inelasticity strongly affects the extraction of fsi phases in models based on quasi - elastic rescatterring .
thus , inelastic events affect model predictions even if rescatterring is of quasi - elastic type only .
on the other hand , fsi phases are often attributed directly to short - distance ( sd ) quark - line diagrams in the hope that this will take into account all inelastic production phenomena .
this belief persists despite justified scepticism about the dominance of short - distance qcd in fsi of @xmath0 decays ( see eg.@xcite ) .
in fact , it is known that the resulting prescription strongly violates such tenets of strong interactions as isospin symmetry @xcite ( see also @xcite ) .
the origin of the problem pointed out in ref .
@xcite is the lack of any correlation between the spectator quark and the products of @xmath16 quark decay . by its very nature
such correlation can not be provided by sd dynamics .
what must be involved here is a long - distance ( ld ) mechanism which ensures that quarks `` know '' about each other .
thus we are led to hadrons , hadron - level dynamics , and inelastic rescatterring effects .
in this paper we perform a simplified analysis of corrections which should be introduced by inelastic rescatterring into the sd - based description of some nonleptonic decays of @xmath0 mesons .
of course , any such analysis must be half - qualitative in nature , because - in the presence of many decay channels - it is well beyond our ability to take them accurately into account .
since it is believed that @xmath0 decays into @xmath17 and @xmath3 may provide some handle on the determination of angle @xmath18 of the unitarity triangle , we shall concentrate on these decays : it is important to know the effects of fsi here . as shown in ref .
@xcite , inclusion of coupled - channel quasi - elastic effects ( @xmath19 ) generates an effective long - distance penguin amplitude comparable in size to the short - distance one .
one may expect that inelastic channels will also contribute to this effect .
we start with an sd - based model of nonleptonic @xmath0 decays . on the basis of standard tree - dominated mechanism for these decays , enriched with the related and well - established models of semileptonic @xmath0 decays , we qualitatively estimate the types and the number of states produced in the first stage of the nonleptonic decay . as in other existing models , we take these states as composed of two ( flying apart ) resonances @xmath4 ( section 2 ) . these resonances are assumed to rescatter into @xmath1 through regge exchange .    in order to provide the basis for an estimate of this rescatterring
, we recall how in a regge picture the unitarity relation involving @xmath20 and other @xmath21 processes looks like .
this enables us to make a rough estimate as to what part of all inelastic @xmath22 processes is due to the @xmath20 transitions and , consequently , how much the situation deviates from the case of ( quasi-)elastic rescatterring . using a rough estimate for the contribution @xmath23 from an average single inelastic intermediate channel @xmath24
, we estimate the number of inelastic channels involved ( section 3 ) .    in section 4
we analyse the behaviour of the @xmath25 and @xmath26 amplitudes as a function of the number of intermediate states considered .
we show that although effects due to each single intermediate state are small , the combined effect of all intermediate states is large . in particular , amplitudes of decays into @xmath3
become significantly larger than those estimated through short - distance penguin amplitudes , to the point of being comparable to the @xmath7 amplitudes .
our conclusions are given in section 5 .
in short - distance approaches to nonleptonic weak @xmath0 decays , the relevant amplitudes are usually expressed as sums of amplitudes corresponding to different types of quark diagrams ( @xmath27 - tree , @xmath28 - colour - suppressed , @xmath29 - @xmath30-exchange , @xmath10 - penguin , @xmath31 - annihilation , @xmath32 - penguin annihilation ) . in this paper
we concentrate on @xmath33 decays of @xmath34 mesons which are initiated by a @xmath35 transition . in these
decays one expects that the dominant contribution comes from the tree diagram @xmath27 , with the main corrections to it provided by the colour - suppressed and penguin diagrams @xmath28 and @xmath10 @xcite .
accordingly , neglecting contributions from other diagrams , the sd amplitudes @xmath36 for @xmath37 decays into a pair of octet pseudoscalar mesons @xmath38 with total isospin @xmath39 are expressed in the su(3 ) symmetry case as @xmath40 inclusion of @xmath41 and @xmath42 into our considerations is mandatory if we want to maintain su(3 ) symmetry @xcite .    from the above formulas one may find sd amplitudes @xmath43 from state @xmath44 into states @xmath45 of given isospin @xmath39 belonging to definite representations @xmath46 of su(3 ) : @xmath47 with @xmath48\nonumber \\
{ \bf w}_0^t&=&[\langle ( \pi \pi ) _ 0 | w |
b^0 \rangle , \langle ( k \bar{k})_0|
w | b^0 \rangle , \langle ( \eta _ 8 \eta _ 8)_0 | w | b^0 \rangle]\end{aligned}\ ] ] and the matrices @xmath49 given by : @xmath50 @xmath51\ ] ] and @xmath52\ ] ] where rows correspond ( from top to bottom ) to @xmath53 for @xmath54 ; @xmath55 for @xmath56 ; and @xmath57 for @xmath58 .
one expects that sd amplitudes @xmath28 and @xmath10 constitute a 10 - 20% correction @xcite . indeed , for the factor @xmath59 in the relation @xmath60 , the short - distance qcd corrections give the value of @xmath61 ( @xmath62 , @xmath63 are wilson coefficients ) , while estimates in ref .
@xcite yield @xmath64 in the range of @xmath65 .
dominance of @xmath27 amplitudes is expected to hold for other decays initiated by @xmath66 as well . for larger invariant masses of @xmath67 and @xmath68 systems
these quark - level states should hadronize mainly as @xmath69 states @xcite .
it is the rescatterring from these states to the final pp state that is of interest to us . since we want to estimate rescatterring at hadron level , we need to know what hadronic states are produced in the first stage of the decay .
there are two groups of hadrons produced : one comes from the decay @xmath70 , the other one originates from the recombination of the @xmath71 quark with spectator @xmath72 ( fig .
1 ) . the values of invariant mass squares @xmath73 and @xmath74 ( see fig
. 1 for definitions ) are not large .
the well - known probability distribution of @xmath16-quark decay is @xmath75 which falls for increasing @xmath76 , the average @xmath77 being @xmath78 .
the values of @xmath74 are smaller .
estimates obtained in various models for semileptonic and nonleptonic decays yield @xmath79 around @xmath80 or so , with the distribution of @xmath81 extending to around 2.0 or 2.5 gev @xcite .
when such a low - mass quark - antiquark state hadronizes , a resonance is produced . in the isgw2 model @xcite for semileptonic decays , exclusive partial widths for the production of lowest - in - mass resonances have been predicted .
the resonances considered were the ground - state mesons ( pseudoscalars @xmath10 : @xmath82 and vector mesons @xmath83 : @xmath84 ) , the p - wave mesons ( tensor mesons @xmath27 : @xmath85 , axial mesons @xmath31 : @xmath86 and @xmath0 : @xmath87 , and scalar mesons @xmath88 : @xmath89 ) , and the @xmath90 and @xmath91 states .
the highest ( nonstrange ) resonance mass explicitly considered in @xcite is below 1.5 gev .
the total partial width for the production of all these resonances is 5 - 6 times larger than the partial width for the production of a pseudoscalar meson .
thus , the average partial width into a resonance is smaller by a factor of 0.7 or so than that for the production of a pseudoscalar meson .
the resonances explicitly considered in @xcite do not saturate the inclusive decay rate which is still about two times larger @xcite .
thus , several other resonances of masses below 2.0 or 2.5 gev should be added to the list given above .
assuming that the average contribution from each one is similar to that just estimated , one expects that the number of types of `` average '' resonances produced should be of the order of 15 .
similar or even larger number of resonance types is expected from the hadronization of the @xmath67 created from the @xmath30-boson .
thus , in nonleptonic b decays , apart from the @xmath1 state , many other resonance pairs must be produced : @xmath4 = @xmath92 , ... , @xmath32 , @xmath93 , @xmath94 , @xmath95 , @xmath96, ... etc .
exact counting of the number of all these two - resonance states is not important for our purposes .
however , it is fairly easy to give an estimate : limiting oneself to resonances of mass smaller than 2 gev , this number will definitely be greater than 10 , probably of the order of a few tens .    clearly , there is no hope that one can reliably calculate the contribution from rescatterring into @xmath1 from each of these intermediate states .
however , one may try to estimate their overall contribution in an average way . for the first stage of the decay process
, we will assume that the only important amplitude is the tree amplitude @xmath27 and that this amplitude is approximately the same for all intermediate states considered ( for low values of @xmath77 that we shall be concerned with later , this is is indeed the case in eq .
( [ eq : bdecay ] ) ) .
although this may seem a very rough assumption for @xmath74 , we shall further see that our general results should be fairly independent of it as long as there is a rather large number of two - resonance states with production amplitudes scatterred around the average .
the next question is how to describe transitions @xmath97 ( or vice versa ) in an average way .
this is what we shall discuss in section 3 .
elastic scatterring and multiparticle production processes are related to each other through unitarity of the s - matrix .
since @xmath0 has spin @xmath98 , we shall work with @xmath99 sector only . in the following
we shall suppress the subscript @xmath100 and its value . in refs .
@xcite it was shown that in the su(3 ) symmetry case , with the effects of coupled channels included , one should work with states @xmath101 belonging to definite representations @xmath102 of su(3 ) .
the unitarity relation for @xmath1 scattering in the @xmath103 partial wave and in su(3 ) representation @xmath46 is : latexmath:[\[\label{eq : unitarity }     @xmath46 , @xmath105 labels states different from @xmath101 .
matrix elements occuring in eq .
( [ eq : unitarity ] ) may be expressed in terms of argand amplitudes @xmath106 as follows @xmath107 apart from the pomeron , there are other regge trajectories , whose exchange in the @xmath108-channel contributes to quasi - elastic scatterring @xmath109 . when the leading non - pomeron exchange - degenerate regge trajectories @xmath110 , @xmath111 , @xmath112 , and @xmath113 and their su(3)-symmetric partners are taken into account , one obtains the @xmath103 partial wave amplitudes @xcite : @xmath114 with @xmath115 where for the leading non - pomeron regge trajectory we use @xmath116 ie .
we have put @xmath117 and , consequently , in eqs ( [ eq : as ] ) and further on both @xmath118 and @xmath108 are in @xmath119 .
the amplitudes @xmath120 are independent of isospin @xmath39 as shown in @xcite , where their sizes at @xmath13 have also been estimated .
the su(3)-symmetric regge residue @xmath121 is fixed from experiment as @xcite @xmath122 while for the pomeron one has @xcite @xmath123 using the above values in eqs ( [ eq : as ] ) one finds @xmath124 while the leading non - pomeron regge contributions _ alone _ are : @xmath125 contribution from elastic scatterring ( the pomeron in regge language ) is independent of @xmath46 : @xmath126 ( cf .
@xmath127 in ref .
@xcite ) . omitting the reggeon - pomeron interference term , the value of the contribution from the leading non - pomeron regge exchange to the unitarity relation of eq .
( [ eq : unitarity ] ) , after averaging over representations @xmath46 , is equal to @xmath128 neglecting reggeon contribution to @xmath12 one obtains from eq .
( [ eq : unitarity ] ) : @xmath129 one may conjecture that contributions to the above sum from reggeon - exchange - induced processes @xmath130 ( where @xmath131 denotes low - lying resonance ) will be of size similar or smaller than @xmath132 . thus , if all inelastic channels @xmath133 were two - resonance states @xmath4 , one would obtain from eq .
( [ eq : numericalunitarity ] ) the number of @xmath134 as the number of states contributing to the unitarity relation .
this should be compared with the estimate of a few tens obtained in the previous section for the number of two - resonance states produced in weak decays of @xmath0 meson .
of course , the estimate of eq .
( [ eq : totalnumber ] ) is probably too low : contributions from transitions @xmath135 for heavier resonances @xmath131 are likely to be smaller and the total number of states may be larger .
for example , with average value of @xmath106 going down by a factor of 0.6 - 0.7 from 0.08 to 0.05 ( section 2 ) , the above estimate of the number of average quasi - two - body states increases from 25 to over 60 . on the other hand , two - resonance states
need not saturate eq.([eq : numericalunitarity ] ) .
thus , there are two important questions which should be answered :    1 .
how many of states @xmath133 in the unitarity relation of eqs ( [ eq : unitarity],[eq : numericalunitarity ] ) are indeed of the form @xmath4 , so that they can contribute to rescatterring in @xmath0 decays , and 2 .
can the effect of all of these states be described in some average way .    at present , elastic and quasi - elastic contributions from long - distance fsi in @xmath0 decays are usually evaluated using the old language of regge theory .
this language was used in the past also for the description of resonance and multiparticle production processes that are of interest to us . in particular , the question of the build - up of elastic scatterring ( pomeron exchange ) as a shadow of inelastic multiparticle production processes ( and thus the very content of the unitarity relation of eqs ( [ eq : unitarity],[eq : numericalunitarity ] ) ) was discussed extensively .
therefore , we must recall the essential elements of an approach which dealt with that problem .
the approach , predominantly occupied with the issue of the unitarization of dual models ( for reviews see @xcite ) , was based on general properties of dual string models and on phenomenological analysis of resonance and multiparticle production data .    in this approach ,
multiparticle production processes occurring in hadron - hadron collisions at high energies are pictured as proceeding through the production of resonances or clusters in a multiperipheral model ( fig .
2 ) with leading non - pomeron reggeons exchanged in between the clusters .
it is the _ sum _ over _ all types _ and _ numbers _ of resonances produced in this way that saturates the s - matrix unitarity relation eq .
( [ eq : unitarity ] ) . for @xmath118 above the inelastic threshold but still small enough , one may limit oneself to the production of just one pair of resonances @xmath4 . with increasing energy @xmath118 , the number of resonances produced in a single collision increases on average .
although cross - section for a production of a particular number of resonances goes down at sufficiently high energy , the sum of cross sections over all possible numbers of resonances remains approximately constant .
this constancy of the total inelastic cross section is ensured by a sufficiently fast increase in the number of all possible quark - line diagrams , ie . in the number of all possible states @xmath133 ( and ways in which they are produced ) @xcite .    in our case , at @xmath136 the model @xcite predicts that states composed of just two resonances are produced in the fraction of @xmath137 cases approximately .
further 35% comes from production of three resonances , etc .
although these numbers are obtained in @xcite for resonances of any mass , the main contribution comes from production of resonances of invariant mass squared smaller than @xmath138 .
( contribution from the production of objects of mass @xmath139 is suppressed as @xmath140 for larger @xmath139 @xcite ) .
the average mass of a resonance produced may be estimated in various ways to be around @xmath141 ( @xmath142 in ref .
@xcite ) , in good agreement with mass @xmath143 or so @xcite , expected for average @xmath81 in the sd - based models of weak decays .
contributions from rescatterring of states with larger values of @xmath77 will be suppressed because production of such states in @xmath1 collisions is not likely .
thus , in a rescatterring process @xmath144 one may expect that the dominant contribution will indeed come from the rescatterring of states composed of two low - mass resonances @xmath24 . translating the above expectation of a 50% share of @xmath4 states in the unitarity relation of eqs ( [ eq : unitarity],[eq : numericalunitarity ] ) into a number @xmath145 of contributing channels which may connect to the state originally produced by the sd dynamics ( ie .
@xmath146 ) , we conclude that this number should be around @xmath147 for average @xmath148 or @xmath149 for average @xmath150 . the latter estimate is probably more realistic since the average size of a contribution from a single @xmath4 channel should diminish with growing resonance masses .
if one accepts that final state interactions can not modify the probability of the original sd weak decay , it follows that vector @xmath151 representing the fsi - corrected amplitudes is related to vector @xmath152 of the original sd amplitudes through @xcite @xmath153 indeed , in the basis of @xmath154-matrix eigenstates @xmath155 the above equation reduces to @xmath156 , ie .
the condition of unchanged probability ( @xmath157 ) admits watson phases only .
the s - matrix may be written in terms of the matrix @xmath158 of amplitudes @xmath106 : @xmath159 we assume that we may treat the fsi - induced corrections to the sd decay amplitudes in a perturbative fashion .
this is in agreement with the ideas of the dominance of sd dynamics .
if this assumption is incorrect , obtaining even half - quantitative predictions will be almost impossible ( cf .
although this assumption may be questioned , it has an important advantage : one may study what happens when the number of contributing two - resonance intermediate states is increased to its expected share ( @xmath160 ) . in agreement with the assumption of a perturbative treatment of rescatterring ( ie . small contribution from @xmath158 ) , we expand the square root in @xmath161 and keep only the first term .
this leads to @xmath162 this is in fact the k - matrix prescription for the estimate of rescatterring effects @xcite . in this prescription ,
final state interactions modify probabilities of the sd amplitudes as can be seen explicitly from eq .
( [ eq : kmatrix ] ) written in the basis of s - matrix eigenstates .
our calculations will be based on eq.([eq : kmatrix ] ) .    in previous sections we considered a simplified picture of contribution from resonance pairs in which all amplitudes @xmath163 were equal in absolute magnitude while their phases were arbitrary . indeed , when one row of the unitarity condition ( ie .
( [ eq : unitarity ] ) ) is discussed , no knowledge of phases is needed . for the purpose of studies of cp - violation the question of phases
is important , however .
therefore , we have to make a very rough estimate of the fsi phases appearing in the strong rescatterring amplitudes @xmath164 in @xmath165 thanks to @xmath166 invariance of strong interactions , the @xmath164 amplitudes are symmetric ( as is the @xmath167-matrix ) . in order to estimate them we have to recall what are the predictions of dual string models for the production of two resonances in high energy @xmath1 collisions . in the appendix of ref .
@xcite it is shown that the dual string model predicts that the amplitude for the @xmath135 production through an uncrossed diagram of fig .
3a will pick up a rotating regge phase resulting from the expression @xmath168 similarly , for the crossed diagrams of fig .
3b one has to remove the @xmath169 sign in the above expression , ie .
the amplitude is real .
thus , the phase - generating factor differs from the familiar one in @xmath12 scatterring ( ie .
@xmath170 ) only by a different scaling factor ( @xmath171 ) in the denominator .
such a dependence of regge amplitudes on the masses of produced resonances has been confirmed in analyses of experimental data @xcite .    at this point
we have to take into account the fact that regge amplitudes describe scatterring of two colliding resonances @xmath172 and @xmath173 in a state of definite momenta into a similar @xmath1 state .
in particular , the produced @xmath1 state is a superposition of partial waves , while we are interested in the @xmath99 sector of the @xmath88-matrix only .
restriction to the @xmath98 sector is achieved by integrating the rescatterring amplitudes @xmath163 with @xmath174 ( with @xmath175 being the angular momentum of the @xmath1 pair ) over the allowed range of scatterring angle @xmath176 , or over the corresponding range of momentum transfer @xmath177 .
angular momentum conservation will then admit states @xmath24 with total angular momentum @xmath98 only .
this will have no effect on the assumed form of the sd decay amplitude since in the sd mechanism the decaying @xmath16 quark does not `` know '' about the spectator quark , and , consequently , about the value of @xmath100 for the whole system .
one calculates that for @xmath178 the minimal value of @xmath108 is @xmath179 , while for @xmath180 one may assume @xmath181 . projecting the @xmath97 regge amplitude onto the @xmath98 sector , ie . integrating regge expressions over @xmath182
, we obtain , up to a common overall normalization factor , the following phase factors :    1 .   for the uncrossed diagram ( fig .
3a ) @xmath183 2 .   for the crossed diagram ( fig .
3b ) @xmath184 where @xmath185 .
the overall normalization of contributions from amplitudes @xmath97 was fixed in section 3 .
for the regge description to be valid , the value of @xmath186 must be large . with @xmath187
, the product @xmath171 should not be greater than @xmath188 , possibly @xmath189 , corresponding to the minimum value of @xmath190 equal to @xmath191 .
for @xmath77 equal to @xmath74 , this corresponds to the maximum value of resonance mass being around 1.6 - 1.8 gev .
these are still reasonable numbers when compared with our previous estimates of @xmath192 . in fig .
4 we show dependence of the phase of @xmath193 on @xmath194 for @xmath195 . if one approximates @xmath196 by their average values of around @xmath197 , one finds that the average value of the product @xmath194 is close to @xmath198 . from fig .
4 we see that for @xmath199 the phase of @xmath200 is zero .
although the phase of @xmath200 at smaller and larger values of @xmath194 deviates from zero quite significantly , these deviations are on average of the order of @xmath201 and , consequently , it is still meaningful to talk about coherent superposition ( ie . with approximately similar phases ) of regge contributions from various intermediate states .    because the quark - line structure of @xmath12 and @xmath202 amplitudes is the same , the eigenvalues with which the above phase factors enter into expressions for @xmath203 in definite su(3 ) representations @xmath46 are the same as those in amplitudes @xmath109 .
the latter were calculated in ref .
@xcite .    on the basis of ref .
@xcite we have therefore @xmath204 the fsi - corrected amplitudes for @xmath0-meson decays into states in definite representations @xmath46 of su(3 ) and with isospin @xmath39 are @xmath205 where @xmath206 are sd weak decay amplitudes into @xmath24 , all assumed approximately equal to @xmath43 ( eq .
( [ eq : pptor ] ) , section 2 ) . as in section 3
, we assume that _ all _ amplitudes @xmath207 are equal to some average value @xmath208 .
although the phase of @xmath200 changes over the range of corresponding @xmath194 , it is very instructive first to approximate it everywhere by a constant , namely its average evaluated at , say , @xmath209 , where @xmath200 is real .
the approximate reality of average @xmath200 is a consequence of the particular value of @xmath118 ( being here equal to @xmath210 ) and not an @xmath118-independent feature of the approach . at @xmath211
we get @xmath212 the sum over all @xmath145 two - resonance states will then yield the contribution from inelastic rescatterring : @xmath213 note that for fixed @xmath214 the average amplitude @xmath215 is inversely proportional to @xmath216 . thus , contribution from inelastic events in eq .
( [ eq : constphaseapprox ] ) is proportional to @xmath216 : the smaller the value of @xmath208 , the larger the summed contribution from all two - resonance states , provided their contribution in the unitarity relation ( eq .
( [ eq : unitarity ] ) ) is fixed by the same value of @xmath214 .
obviously , this is a general feature of any perturbative treatment of rescatterring contribution from several channels , the reason being the _ linear _ nature of perturbatively treated fsi as compared to the _ quadratic _ nature of the unitarity relation .
thus , estimates of fsi effects given here are most likely estimates from below : although amplitudes @xmath217 corresponding to rescatterring from states composed of resonances of larger masses are expected to be smaller , their combined rescatterring effect should be relatively larger .
the amplitudes for @xmath0 decays into @xmath218 , @xmath219 are calculated from the inverse of eq .
( [ eq : pptor ] ) : @xmath220 in table 1 we give the values of @xmath221 calculated from the above equation with the approximation of average phase ( eq .
( [ eq : constphaseapprox ] ) ) for @xmath222 and @xmath223 ( @xmath224 respectively ) .
predictions of the model without inelastic rescatterring , ie . with quasi - elastic fsi ( @xmath1 intermediate states )
only , are also given for comparison .
amplitudes are given in units of input tree amplitude @xmath27 .    from table 1
one can see that quasi - elastic fsi do not affect the decays @xmath225 strongly .
the amplitude moduli are somewhat smaller than those without the fsi .
this is due mainly to the @xmath226~$ ] factor originating from pomeron exchange .
phase changes are small . for the @xmath227 decays , quasi - elastic fsi
affect the amplitude significantly : the amplitude driven in sd - dynamics by the penguin diagram ( here vanishing ) receives contribution from the coupled - channel - generated long - distance penguin ( hereafter denoted @xmath228 ) @xcite .
using @xmath229 to estimate the effective ld penguin , the size of @xmath228 is ( as in ref .
@xcite ) of the order of 5% of the tree amplitude @xmath27 , thus permitting significant interference effects with the short - distance penguin amplitude when the latter is taken from standard sd estimates .
the main results of this work are given in the last two columns of table 1 .
* one can see that the @xmath225 amplitudes , when compared with their estimates taking into account quasi - elastic fsi only , increase in absolute magnitude by a factor between 1 and 2 .
this is due to the additional contribution coming from rescatterring chain @xmath230 .
* an important change can be seen in phase sizes : they are now one order of magnitude larger than in the quasi - elastic case .
the origin of this effect is as follows .
the @xmath30 amplitude is composed of two parts : one ( approximately real ) is the sd amplitude @xmath231 weakly suppressed by elastic ( quasi - elastic ) rescatterring , while the other contains the contribution from the inelastic @xmath232 rescatterring .
the latter part , being proportional to @xmath233 , is mainly imaginary ( for approximately real @xmath106 ) . with a large contribution from inelastic rescatterring
, the resulting phase must therefore be large .
* finally , we see that the @xmath234 amplitude becomes much larger than in the case of quasi - elastic rescatterring only , and is dominantly imaginary .
since the @xmath234 amplitude is fed from no - hidden - strangeness states @xmath235 , one might be tempted to compare this with ref .
namely , it was shown there for the case of a simple two - channel s - matrix that the phase of amplitude in channel 1 is large if the particle originally decays to channel 2 . in our case , although a similar result holds , it is not general - it depends on the value of @xmath118 at which amplitudes @xmath163 are evaluated .
note that inelastic rescatterring renders the @xmath227 amplitude larger than the sd penguin estimate : in fact , it is comparable to the @xmath236 amplitude .    in our approach the full ( fsi - corrected ) amplitudes for @xmath237
decays vanish since : + 1)we have assumed that only the tree amplitude @xmath27 is nonzero , and + 2 ) the rescatterring from the isospin @xmath238 state of @xmath239 , which in principle might feed the final @xmath240 channel , is zero .
state of @xmath4 is antisymmetric ( ie . of the form @xmath241 ) , while the rescatterring contribution due to the uncrossed diagram of fig .
3a is zero when evaluated in between the antisymmetric state @xmath242 and the symmetric state @xmath243 ( for the definitions of states , see eq .
( 1 ) of ref .
@xcite ) . rescatterring through diagram of fig .
3b can not change the type of quarks and , consequently , can not induce a transition from the no - hidden - strangeness state @xmath244 into @xmath245 . ]
+ the dependence of the ld rescatterring - induced effective penguin on the isospin channel - with vanishing ( large ) effects in i=1 ( 0 ) channel - should be compared with the sd mechanism which assigns the same size and phase to penguin amplitude @xmath10 in both the @xmath246 and @xmath240 decay channels . this is a general feature of long - distance dynamics : the size and phase of a quark - level diagram depend on what isospin ( su(3 ) ) amplitude it contributes to @xcite .
although the above expectation of fsi effect larger than that naively expected is fairly general , we have to analyse in some detail our assumption of replacing the @xmath171-dependent phase of @xmath200 with an average ( and vanishing ) phase .
such an assumption would be well justified if a large fraction of @xmath24 states led to phases close to the average . since for a given value of the product @xmath194 the phase is fixed , we need to know the density of two - resonance states as a function of @xmath247 .
regge models @xcite predict that the dependence of @xmath248 on @xmath249 is proportional to @xmath250 for larger values of @xmath251 .
this fall of distribution for larger mass values should be clearly visible already at the beginning of the region where using regge description becomes sensible , ie . at @xmath252 or so .
the distribution of @xmath74 in sd models vanishes even faster ( it is fairly negligible above @xmath253 ) .
since we want to estimate corrections to sd models , direct use of mass distributions generated in sd - models might seem to be the simplest and most natural choice .
furthermore , as rescatterring of states with larger values of @xmath77 is suppressed by the small size of the contribution from transitions @xmath254 , one might use the sd distribution of @xmath74 as that of @xmath77 as well .
the problem is , however , that a significant part of the sd distribution of @xmath74 ( or @xmath77 ) corresponds to @xmath255 . on the other hand ,
regge amplitudes assume the form of @xmath256 only for @xmath257 . for @xmath258 , one should replace @xmath259 with @xmath260 .
thus , if regge phases are to be reasonably evaluated we must use an appropriately modified distribution of @xmath74 .
we model this situation in the simplest possible way : by assuming that the distribution of @xmath261 vanishes below @xmath262 and above @xmath263 , while in between these values it is given by @xmath264 this yields the average value of `` effective '' @xmath139 equal to @xmath265 , which is also the value obtained from the accmm distribution of @xmath81 ( see eg . ref .
@xcite ) if contributions from @xmath266 are replaced with contributions at @xmath267 . using @xmath268 as in eq . ( [ eq : mmdistr ] )
, it is straightforward to evaluate the distribution @xmath269 of @xmath171 .
this distribution ( fig .
5 ) is peaked at @xmath270 , its median being at @xmath271 , and the average @xmath272 at @xmath273 .
the tail above @xmath189 contributes a few percent only .
thus , the assumption of @xmath274 used in our previous discussion appears quite reasonable .    with two - resonance states
@xmath133 spread in @xmath171 according to distribution @xmath269 ( larger values of @xmath275 correspond to states with appropriately larger values of @xmath194 ) , it is simple to analyse the predictions of the model numerically .
we are interested in the question how the fsi effects change when heavier and heavier intermediate states @xmath133 are included .
for decays into states in definite representations of su(3 ) , the amplitudes of interest to us are therefore : @xmath276 where @xmath277 is the number of intermediate inelastic states considered ( @xmath278 ) and @xmath279 are given in eq .
( [ eq : deltas ] ) .    in fig .
6 we present predictions of the @xmath277-dependent version of eq .
( [ eq : inverse ] ) for the absolute values and phases of the amplitudes of @xmath0 decays into @xmath280 and @xmath246 states .
these predictions are given as a function of the number @xmath277 of intermediate states considered .
we show the case with @xmath281 ( for @xmath282 one obtains very similar plots with features discussed below being even more pronounced ) .
the point most to the right in each plot ( ie . at @xmath283 ) corresponds to a large value of @xmath171 , where regge approximation breaks down .
consequently , this point should be discarded .
one can see that for @xmath284 the size of the amplitude does not depend very strongly on @xmath277 and is close to the value of the input fsi - free amplitude @xmath285 , which is @xmath286 .
furthermore , the fsi - induced phase is still relatively small ( of the order of @xmath287 ) .
this can not be said of the @xmath227 process . here , the absolute value of the amplitude grows fast with the increasing number of intermediate states .
in addition , already at @xmath288 or so , the phase becomes close to @xmath289 , not far from the previous result of @xmath290 ( table 1 ) obtained for constant phases .
one can also see that when the number of intermediate states approaches its maximum allowed value , the size of the @xmath227 amplitude becomes significant when compared to the tree @xmath284 amplitude . as the number of intermediate states increases , the long - distance - induced penguin amplitude starts to dominate over the sd one ( which was estimated at 0.04 to 0.20 of the tree amplitude , @xcite ) .
furthermore , assuming applicability of eq .
( [ eq : qdbdecays ] ) , one estimates from fig . 6 that at large @xmath277 the effective long - distance penguin amplitude @xmath228 is around @xmath291 .
this should be compared with the input or effective @xmath27 amplitudes which are around @xmath292 .
thus , the long - distance effective penguin amplitudes become really large .
large rescatterring effects obtained above stem from an _ approximately coherent superposition of contributions from several intermediate channels_. one might argue that in the calculation of this paper sizes of individual rescatterring amplitudes @xmath106 ( ie .
contributions from individual intermediate channels ) are overestimated .
however , the trend of the results expected for the case of smaller amplitudes @xmath106 may be seen from table 1 and eqs .
( [ eq : constphaseapprox],[eq : constphaseapprox1 ] ) .
the column of @xmath282 corresponds to a smaller value of average amplitude @xmath293 .
the connection existing between the values of @xmath145 to @xmath294 stems from the assumption that at energy @xmath13 a two - resonance state is produced in approximately half of all inelastic @xmath1 collisions .
as long as this fraction ( @xmath214 ) is kept constant and the perturbative treatment of fsi is valid , smaller values of amplitudes @xmath106 - after summing over all intermediate channels - result in rescatterring corrections to weak decays larger than naively expected .
one may also try to estimate roughly the rescatterring contribution in a direct way ( ie . without using the value of @xmath214 ) by just adding rescatterring contributions with amplitudes @xmath217 assumed to be of the order of @xmath295 or so ( in agreement with experiment ) .
then , for the number of average intermediate states taken as equal to the number of final states in standard models of sd decays , ie . of the order of @xmath296 or more , one is bound to obtain a large rescatterring effect .
please note that our whole approach starts with the generally accepted features of the sd decay ( and resonance production ) mechanism and assumes that subsequent fsi may be treated perturbatively .
we have shown that even in this case the corrections tend to become large .
of course , if they are too large , the whole perturbative scheme of their estimation ( as well as the sd mechanism for the description of @xmath0 decays ) ceases to be viable .
we have discussed multichannel inelastic rescatterring effects in @xmath0 decays into @xmath17 and @xmath297 .
generally accepted features of short - distance decay mechanism have been assumed as part of our input .
these assumptions included estimates of the types and number of resonances produced in two - resonance states initiated by the @xmath298 transition .
rescatterring of these two - resonance states into the final state consisting of one pair @xmath1 of pseudoscalar mesons was evaluated under the assumption that such fsi may be treated perturbatively .
the basis for this evaluation was provided by existing knowledge about how the inelastic multiparticle ( resonance ) production in @xmath1 collisions is correlated with elastic @xmath1 ( pomeron exchange ) scatterring .
this knowledge permitted us to estimate that at @xmath13 @xmath1 scatter into a two - resonance @xmath4 state in approximately 50% of cases .
thus , the total size of rescatterring from @xmath4 to @xmath1 was fixed .    using the regge model for the description of the @xmath97 processes
, we have shown that the rescatterring contributions from the individual intermediate channels add approximately coherently . as a result , the combined effect of rescatterring through many two - resonance intermediate states was shown to be quite large .
this was demonstrated under the assumption that fsi may be treated perturbatively . if that assumption is overoptimistic , reliable estimate of ( presumably even larger ) fsi effects will almost certainly be much more difficult .    in our calculations
the amplitude for the decays @xmath234 was induced by ld rescatterring from no - hidden - strangeness isospin 0 states produced via short - distance tree diagram .
this fsi - induced amplitude was shown to be larger than its short - distance penguin counterpart .
the phase of the ld amplitude was estimated to be around @xmath289 .
the author would like to thank prof .
p. singer and the department of physics at the israel institute of technology , haifa , where some ideas of this paper have originated , for their hospitality .
j. f. donoghue et al . , phys .
77 , 2178 ( 1996 ) .
m. suzuki and l. wolfenstein , phys . rev .
d60:074019 ( 1999 ) .
m. suzuki , _ final state interactions in heavy hadron decay _ , talk given at 3rd international conference on b physics and cp violation , taipei , taiwan , 3 - 7 dec 1999 , hep - ph/0001170 .
a. n. kamal and c. w. luo , phys .
d57 , 4275 ( 1998 ) .
a. n. kamal , phys .
d60:094018 ( 1999 ) .
grard , j. weyers , eur .
j. c7 , 1 ( 1999 ) .
l. wolfenstein , phys .
d52 , 537 ( 1995 ) .
p. enczykowski , phys .
b460 , 390 ( 1999 ) .
m. gronau , o. f. hernandez , d. london , j. l. rosner , phys .
d52 , 6356 ( 1995 ) .
p. enczykowski , acta phys .
b28 , 1605 ( 1997 ) .
m. gronau , phys .
b300 , 163 ( 1993 ) .
a. v. dobrovolskaya , a. b. kaidalov , k. a. ter - martirosyan , and v. r. zoller , phys .
b229 , 293 ( 1989 ) .
l. wolfenstein , phys .
d43 , 151 ( 1991 ) .
g. altarelli et al . , nucl .
b208 , 365 ( 1982 ) .
v. barger , c. s. kim , and r. j. n. phillips , phys .
b251 , 629 ( 1990 ) .
d. scora and n. isgur , phys .
d52 , 2783 ( 1995 ) .
i. bigi , r. d. dikeman , n. uraltsev , eur .
j. c4 , 453 ( 1998 ) .
m. battaglia , `` study of @xmath299 decays with an inclusive generator '' hep - ex/0002040 .
hong - mo chan , sheung tsun tsou `` dual unitarization : a new approach to hadron reactions '' , rl-76 - 080 , sep.1976 , published in bielefeld sum .
1976 , 83 ; g. f. chew , c. rosenzweig , phys . rept . 41 , 263 ( 1978 ) .
huan lee , phys .
30 , 719 ( 1973 ) ; g. veneziano , phys . lett .
b43 , 413 ( 1973 ) .
hong - mo chan , j. e. paton , sheung tsun tsou , sing wai ng , nucl .
b92 , 13 ( 1975 ) .
hong - mo chan , j. e. paton , sheung tsun tsou , nucl .
b86 , 479 ( 1975 ) .
n. sakai , nucl .
b99 , 167 ( 1975 ) .
p. hoyer , r. g. roberts , and d. p. roy , nucl .
b56 , 173 ( 1973 ) a. c. irving , r. p. worden , phys .
34 , 117 ( 1977 )                fig.6 + dependence of fsi - induced effects on the number of intermediate channels included
. +   + table 1 .
+ effects of inelastic rescatterring on @xmath304 amplitudes in average phase approximation .
amplitudes are in units of input tree amplitude @xmath27 .