in this article we consider self - similar solutions to smoluchowski s mean - field model for coagulation @xcite that is given by the equation @xmath12\ , , \label{s1e1}\end{split}\ ] ] where @xmath13 denotes the number density of clusters of size @xmath14 at time @xmath15 .
the kernel @xmath16 describes the rate of coalescence of clusters of size @xmath14 and @xmath17 and subsumes all the microscopic aspects of the coagulation process .
it is well - known , cf .
for example the review article @xcite and references therein , that if the kernel grows at most linearly , then the initial value problem for ( [ s1e1 ] ) is well - posed for initial data with finite mass and the mass is conserved for all times .
for homogeneous kernels it is furthermore expected that solutions converge to a self - similar form for large times .
in fact , for the kernels @xmath18 and @xmath19 , the constant and additive kernel respectively , this issue has by now been completely solved .
it has been established @xcite that besides explicitly known exponentially decaying self - similar solutions also solutions with algebraic decay exist . in @xcite
their domains of attraction could also be completely characterized .
however , not much is known about self - similar solutions for other kernels than these solvable ones .
existence of fast - decaying self - similar solutions has been established for a large class of kernels @xcite , and local properties of such solutions have been investigated in @xcite , but it is not known whether they are unique in the class of solutions with finite mass .
it is also not clear , not even on the formal level , whether self - similar solutions with algebraic decay exist .    in this article
we will consider mass - conserving self - similar solutions to coagulation equations with the following product kernel of homogeneity @xmath20 : @xmath21 existence of a solution for this kernel has been established in @xcite .
the purpose of this paper is to give an asymptotic description of such a solution in the regime @xmath22 .
let us describe briefly what is known for self - similar solutions for kernels as in ( [ s1e2 ] ) and compare it to results for the so - called sum kernels . to fix ideas , we restrict ourselves to kernels of the form @xmath23 even though most of the results that we mention also apply to more general kernels with the same growth behavior as the ones in ( [ s1e2b ] ) .
we will need in particular to distinguish between the case @xmath24 , the sum kernel , and the case @xmath25 , the product kernel .
it has long been predicted @xcite that self - similar solutions for kernels as in ( [ s1e2b ] ) exhibit a singular power - law behavior of the form @xmath26 with @xmath27 in the case @xmath24 , and @xmath28 for the case @xmath25 .
this has been rigorously proved for the case @xmath24 in @xcite and for the case @xmath25 in @xcite . as has been pointed out in @xcite , the next order behavior for small clusters in the case @xmath24
can then easily be established and is as predicted by the physicists . in the case
@xmath25 , however , the next order behavior has not been known @xcite , and only numerical simulations suggested that it is oscillatory @xcite .    our goal in this paper is to describe formally how to construct mass - conserving self - similar solutions for kernels as in ( [ s1e2 ] ) in the limit @xmath22 and by this also describe their asymptotic behavior on the whole positive real line . while we believe that such self - similar solutions are unique
, our approach only gives the construction of one such solution .
in forthcoming work @xcite we will also show how to make this construction rigorous .
we will see that in the limit @xmath22 the oscillatory character of the solutions becomes very explicit and can be interpreted in terms of simple odes coupled with explicitly solvable equations .    in the limit @xmath29
the product kernel becomes close to the sum kernel @xmath30 .
as described above the power law for the sum kernel is different from the one for the product kernel .
this is due to the fact that , thinking in terms of coagulating particles , the physical behaviour of these particles is different in these two cases .
indeed , the power law obtained for the sum case shows that particles with small @xmath8 coalesce mostly with particles that are much bigger than themselves . on the contrary ,
in the case of the product kernel small particles interact mostly with the ones having a comparable size .
the point @xmath31 can be thought as a bifurcation point where both kind of behaviours take place simultaneously .
we are now going to derive the equations that are solved by self - similar solutions .
it is known that the coagulation equation ( [ s1e1 ] ) can be written in the following conservative form that makes the conservation of the number of monomers @xmath32 transparent .
@xmath33    using the self - similar variables @xmath34 as well as the form of the kernel ( [ s1e2 ] ) , equation ( [ s2e1 ] ) becomes @xmath35 integrating this equation , assuming decay of the solutions as @xmath36 and imposing absence of particle fluxes , we obtain @xmath37 the balance of the terms on both sides of ( [ s2e2 ] ) suggests the following power law behaviour for @xmath38 near the origin , @xmath39 where @xmath40 is the classical beta function @xcite .
for further reference we notice that @xmath41    in order to remove the power law behaviour we introduce the new function @xmath42 then @xmath9 solves @xmath43 in the rest of the paper we will study solutions of ( [ s2e5 ] ) in the limit @xmath22 that satisfy @xmath44 as @xmath45 .
notice that ( [ s2e5 ] ) has the explicit constant solution @xmath46 that corresponds to the well - known solution with infinite mass @xmath47 in the formulation ( [ s2e2 ] ) .
a main idea in our approach is to reformulate the problem of finding solutions of ( [ s2e5 ] ) as a shooting problem for a volterra integro - differential equation .
indeed , we can rewrite ( [ s2e5 ] ) as @xmath48 differentiating ( [ s3e5 ] ) we obtain the equations @xmath49 thus the solution of ( [ s2e5 ] ) solves ( [ s3e7 ] ) with @xmath9 as in ( [ s3e4 ] ) , that has the advantage of being a volterra integro - differential equation .
however , the problem ( [ s3e4 ] ) , ( [ s3e7 ] ) admits solutions that in general do not decay as @xmath50 and therefore they do not provide solutions of ( [ s2e5 ] ) . therefore , among all the solutions of ( [ s3e4 ] ) , ( [ s3e7 ] ) we must select those satisfying @xmath51      we now sketch the asymptotics of the solution that we are going to construct .
it passes , roughly speaking , through three stages that are described in detail in sections [ s.zero ] , [ combregion ] and [ s.asymptotics ] respectively .
section [ s.zero ] discusses the behavior of the solution near @xmath52 .
as pointed out before , near the origin the solution is oscillatory and behaves as @xmath53 where @xmath54 and @xmath55 as @xmath22 .
our analysis consists of constructing a solution by starting with @xmath9 as in ( [ oscillations ] ) and using @xmath56 as a shooting parameter .
however , since equation ( [ s2e5 ] ) is invariant under the rescaling @xmath57 for @xmath58 , we can restrict the range of @xmath59 to @xmath60 ( see also the comment in section [ ss.variables ] ) . in this regime near the origin
we can approximate ( [ s3e4 ] ) and ( [ s3e5 ] ) by a nonlinear ode ( see ( [ s5e6 ] ) ) , that is a perturbation of a simple ode system . for the perturbed system we can use an adiabatic approximation to compute the increase of an associated energy @xmath61 along trajectories . as we discuss in section [ ss.validity ] this first ode approximation is valid as long as @xmath62 . when @xmath63 we enter a new regime that is described in section [ combregion ] . in this regime
@xmath9 develops peaks of height @xmath10 and width of order one that are separated by wide regions in which @xmath9 is small .
the regime where @xmath9 is small is again described by an ode ( cf .
section [ ss.ode ] ) , while the peaks are described by an integro - differential equation ( see section [ ss.integro ] ) .
the analysis of these respective regimes is done in sections [ ss.intanalysis ] and [ ss.odeanalysis ] , while their coupling is described in section [ asymptpeaks ] . in section [ s.shooting ]
we will then show by a continuity argument that there exists a shooting parameter such that the corresponding solution @xmath9 converges to zero as @xmath45 . in the third regime , that is discussed in section [ s.asymptotics ] , @xmath9 then decays exponentially fast to zero .
it has been observed in @xcite that the self - similar solutions associated with ( [ s1e1 ] ) exhibit oscillations for small values of @xmath65 this oscillatory behaviour can be seen as follows . if we make the ansatz @xmath66 with some @xmath67 , plug this into ( [ s2e5 ] ) and keep only the leading order terms , we obtain that @xmath68 must satisfy @xmath69 after some computations , we find that this is equivalent to @xmath70    one can show that @xmath71 has two complex conjugate roots @xmath72 with positive real part and nonzero imaginary part .
indeed , in the limit @xmath6 this can be seen easily as the dominating terms in ( [ s3e3 ] ) are , assuming @xmath73 , @xmath74 and thus the roots of @xmath71 satisfy @xmath75    consequently we expect the following asymptotics for the solutions of ( [ s2e5]):@xmath76 for suitable real constants @xmath77 and @xmath78    notice that the asymptotics ( [ s2e6 ] ) indicate that two consecutive local maxima of @xmath9 , called @xmath79 and @xmath80 , satisfy @xmath81 for small @xmath82 in comparison , to double the amplitude of ( [ s2e6 ] ) we must multiply @xmath8 by a number of order @xmath83 . thus ,
the behaviour ( [ s2e6 ] ) is essentially oscillatory , with a growth in the amplitude that takes place on a much larger scale .      for our forthcoming analysis
it will be more convenient to use the following variables : @xmath84 then ( [ s3e4 ] ) , ( [ s3e7 ] ) becomes @xmath85 \ , , \label{s3e10}\\ \frac{du}{dx } &    = -\left (   1{-}\lambda\right )   u   + \left ( 1{-}\lambda\right )   h\ , , \qquad \frac{dv}{dx }     = \lambda v   -\lambda h\ , , \label{s3e12}\ ] ] where @xmath86   \left (   x\right )   = h_{\lambda}\int_{-\infty}^{x}dye^{\left (   1-\lambda\right )   \left (   y - x\right )   } h\left (   y\right ) \int_{\log\left (   e^{x}-e^{y}\right )   } ^{x}dze^{-\lambda\left (   z - x\right ) } h\left (   z\right)\ , .
\label{s3e13}\ ] ]    the system ( [ s3e10])-([s3e13 ] ) has two explicit solutions , namely @xmath87 @xmath88 and @xmath89 .
our goal is to construct a solution of ( [ s3e10])-([s3e13 ] ) for small values of @xmath90 with the property @xmath91 we will see that a solution that satisfies ( [ s4eb ] ) decays exponentially fast to zero .
this in particular implies that the corresponding self - similar solution @xmath38 has finite mass .
equations ( [ s3e8 ] ) , ( [ s3e9 ] ) and ( [ s2e6 ] ) yield the following asymptotics of the solutions as @xmath92 : @xmath93 where @xmath94 is as in ( [ d1e2 ] ) and @xmath95    the complex number @xmath59 determines the solution of ( [ s3e10])-([s3e13 ] ) uniquely .
however , we notice that ( [ s3e10])-([s3e13 ] ) is invariant under translations @xmath96 and thus all the complex numbers in the spiral @xmath97  yield the same solution up to translations .
thus we can identify the solutions of ( [ s3e10])-([s3e13 ] ) satisfying ( [ s4ea ] ) , up to translations , with the set of positive real numbers contained between two consecutive intersections of @xmath98 with the real axis .
in particular , given @xmath99 we have , since @xmath100 , that the next consecutive point in @xmath101 is @xmath102 therefore , there is a one - to - one correspondence between the points in the interval @xmath103 and the solutions of ( [ s3e10])-([s3e13 ] ) satisfying ( [ s4ea ] ) .
the main result of this paper is that we show , using asymptotic arguments , that for small @xmath90 a value of @xmath104 exists such that ( [ s4eb ] ) holds .
the key idea in computing the asymptotics of the solutions of ( [ s3e10])-([s3e13 ] ) is to obtain suitable approximations of the operator @xmath105   \left (   x\right )   $ ] in ( [ s3e13 ] ) for small @xmath90 .
the formula for @xmath106 ( cf .
( [ s2e3 ] ) ) , combined with ( [ s2e3b ] ) and ( [ s3e13 ] ) , suggests the approximation @xmath86
\left (   x\right )   = \lambda\int_{-\infty}^{x}e^{\left ( y - x\right )   } h\left (   y\right )   dy\int_{\log\left (   e^{x}-e^{y}\right )   } ^{x}h\left (   z\right )   dz\ , . \label{e1}\ ] ] we will discuss the consistency of this assumption in section [ s.consistency ] .    the operator @xmath107   \left (   x\right )   $ ] can be further approximated if @xmath108 changes significantly faster than @xmath109 , as we saw is the case for @xmath110 .
then it is natural to approximate @xmath111 by @xmath112 and we obtain    @xmath86   \left (   x\right )   \sim\lambda h\left (   x\right ) \int_{-\infty}^{x}e^{\left (   y - x\right )   } dy\int_{\log\left (   e^{x}-e^{y}\right )   } ^{x}h\left (   z\right )   dz\ \ ] ] and changing the order of the integrals @xmath86   \left (   x\right )   \sim\lambda h\left (   x\right ) \int_{-\infty}^{x}h\left (   z\right )   e^{z - x}dz\ , .
\label{s5e1}\ ] ]    from ( [ s3e12 ] ) we have @xmath113 which together with ( [ s5e1 ] ) leads to @xmath86     \sim\lambda h   u   \ , .
\label{s5e2}\ ] ] using ( [ s5e2 ] ) as well as ( [ s2e3 ] ) and ( [ s2e3b ] ) we obtain the following approximation of ( [ s3e10 ] ) up to order @xmath90 : @xmath114 and plugging this approximation into ( [ s3e12 ] ) , ( [ s3e13 ] ) we obtain@xmath115 as ( [ d1e2 ] ) and ( [ s4e0 ] ) imply that the changes of @xmath116 compared to those of @xmath117 are of order @xmath118 , we make the change of variables @xmath119 that transform ( [ s5e3 ] ) , ( [ s5e4 ] ) up to order @xmath118 into @xmath120      let us denote by @xmath121 any solution of ( [ s5e6 ] ) .
we first notice that the asymptotics of the solutions of ( [ s5e6 ] ) as @xmath122 agree with those obtained in ( [ s4e0 ] ) .
indeed , in the limit @xmath123 we can compute the asymptotics of @xmath124 as @xmath125 using ( [ s4e0 ] ) , ( [ d1e2 ] ) and ( [ s5e5 ] ) : @xmath126 linearizing ( [ s5e6 ] ) near @xmath127 we obtain the same asymptotics for the functions @xmath128 as @xmath129 this gives the desired matching between the solutions of the linearization of ( [ s3e10])-([s3e13 ] ) around @xmath130 and the solutions of the approximated problem ( [ s5e6 ] ) .
the approximation ( [ s6e2 ] ) is valid as long as @xmath131 is small .
however , a detailed analysis of the nonlinear problem ( [ s5e6 ] ) is needed if @xmath132 becomes of order one . in order to describe the solutions of ( [ s5e6 ] ) in this regime
, we notice that this equation is a perturbation of the equation @xmath133 this equation can be explicitly integrated , since the following quantity is conserved along trajectories : @xmath134    the solutions of ( [ s6e3 ] ) are periodic , but this behavior is not compatible with the exponential growth obtained in ( [ s6e2 ] ) .
this growth is due to the increase in @xmath61 produced by the terms of order @xmath135 in ( [ s5e6 ] ) .
we now compute this change of energy for values of @xmath124 of order one . using ( [ s5e6 ] ) and ( [ s6e4 ] ) we obtain @xmath136   \label{s6e5}\,.\ ] ]    since the solutions of ( [ s5e6 ] )
are close to those of ( [ s6e3 ] ) during finite time intervals , we can adiabatically compute the change in @xmath137 more precisely , the solutions of ( [ s6e3 ] ) satisfying ( [ s6e4 ] ) are periodic with a period @xmath138 given by @xmath139 where the functions @xmath140 are defined as the roots of the equation @xmath141 for any given value of @xmath142 and @xmath143 integrating ( [ s6e5 ] ) and using the adiabatic approximation , we can then approximate the change of @xmath61 during an interval of length @xmath144 as @xmath145   d\xi=\phi\left (   e\right ) \label{s6e7}\ ] ] where @xmath146 is a solution of ( [ s6e3 ] ) satisfying ( [ s6e4 ] ) .
it turns out that the function @xmath147 is positive for any @xmath148 this follows , using the second equation in ( [ s6e3 ] ) , via @xmath149   d\omega -\int_{-\sqrt{2e}}^{\sqrt{2e}}\left [   \left (   1-u_{+}\left (   e,\omega\right ) \right )   + \omega^{2}\right ]   d\omega\,,\ ] ] whence @xmath150    we can compute the asymptotics of ( [ s6e8 ] ) as @xmath151 using ( [ s6e6 ] ) .
notice that the leading contribution to the integral is given by @xmath152 since @xmath153 then @xmath154      we have obtained that , as long as the approximation ( [ s5e6 ] ) is valid , we can approximate the functions @xmath155 by @xmath156 while @xmath61 is computed via the iterative recursion ( [ s6e7 ] ) .
due to ( [ s6e8a ] ) these values of @xmath61 increase an amount of order @xmath157 in each period of size @xmath158 we now proceed to discuss the range of validity of the different approximations that have been used to derive ( [ s6e7 ] ) and ( [ s5e6 ] ) .    to examine the range of values of @xmath61 for which the adiabatic approximation is valid , we need to show that the value of @xmath159 remains close to @xmath160 for @xmath161 due to ( [ s6e6 ] ) and ( [ s6e8 ] ) the adiabatic approximation is valid as long as @xmath162 let us remark that ( [ s6e9 ] ) does not follow immediately from the inequality @xmath163 because the right - hand side of ( [ s6e5 ] ) does not have a sign and therefore it is not clear that @xmath164 for some @xmath165 a priori it is not possible to rule out the possibility of big changes of @xmath166 along each cycle balanced in such a way that the overall change of @xmath61 along the cycle is just @xmath167 however , as shown in appendix [ a.energy ] , it turns out that ( [ s6e10 ] ) holds true .    using ( [ s6e7 ] ) , ( [ s6e8a ] ) and ( [ s6e10 ] ) it follows that the adiabatic approximation is valid as long as @xmath168 and this is satisfied as long as @xmath169    we now study the range of validity of the approximation ( [ s5e2 ] ) that is the main ingredient in deriving the approximate problem ( [ s5e6 ] ) .
the main assumption used in the derivation of ( [ s5e2 ] ) is that the characteristic length scale for which @xmath170 has changes comparable to itself is much larger than one . due to ( [ s5e3a ] )
it follows that as long as @xmath171 remains close to @xmath172 the characteristic length scale for @xmath173 is the same as for @xmath174 notice that @xmath175 is small if @xmath176 since we have @xmath177 , the condition @xmath178 holds if ( [ s6e11 ] ) is true , and so the characteristic length scales for @xmath173 and @xmath179 are the same under the assumption ( [ s6e11 ] ) . for @xmath180 of order one
, the evolution of ( [ s6e3 ] ) takes place in the length scale @xmath181 or equivalently for changes in @xmath182 of order @xmath183 therefore , the condition in ( [ s6e11 ] ) that allows us to obtain ( [ s5e2 ] ) is immediately satisfied and we can restrict our analysis to the case @xmath184 the first equation in ( [ s6e3 ] ) combined with ( [ s5e5 ] ) implies that @xmath185    the left - hand side measures the relative variations in @xmath186 with respect to changes in @xmath187 the characteristic length scale that describes the changes in this quantity is large as long as @xmath188 which is again true if ( [ s6e11 ] ) holds .
therefore , if the condition ( [ s6e11 ] ) is satisfied , we can apply simultaneously all the approximations that lead to ( [ s5e6 ] ) , the adiabatic approximation yielding ( [ s6e7 ] ) as well as the condition @xmath189 however , the three assumptions break down simultaneously if @xmath61  becomes of order @xmath190 .
since ( [ s6e7 ] ) and ( [ s6e8a ] ) imply that @xmath61 increases to arbitrarily large values , the failure of ( [ s6e11 ] ) happens for every solution of ( [ s3e10])-([s3e13 ] ) satisfying ( [ s4e0 ] ) .
in order to describe the solutions of ( [ s3e10])-([s3e13 ] ) when @xmath191 we introduce a new group of variables .
it turns out to be convenient to split the regime into the one where @xmath186 ( or @xmath192 ) is of order one or smaller , and the one where @xmath186 ( or @xmath192 ) is large .
notice that ( [ s6e6 ] ) implies that for a given value of @xmath193 the minimum of @xmath186 scales as @xmath194 , and the maximum scales as @xmath10 .
we remark that the forthcoming analysis will show that as long as @xmath171 remains of order one and is not small the values of @xmath192 and @xmath186 have the same order of magnitude .
this assumption will be made implicitly in all the remaining computations and will be justified by the self - consistency of the derived asymptotics .
we begin by studying the case in which @xmath186 ( and @xmath192 ) are large .
since we are interested in the case @xmath195 , equation ( [ s6e6 ] ) suggests the rescaling @xmath196 then ( [ s6e12 ] ) yields a characteristic length scale for @xmath182 of order one .
this suggests introducing the new set of variables @xmath197 plugging ( [ s7e1 ] ) into ( [ s3e10])-([s3e13 ] ) we obtain @xmath198\ , ,   \label{s7e4}\\ \frac{d\mathcal{u}}{dx }   &   = -\left (   1{-}\lambda\right )   \mathcal{u }    + \left (   1{-}\lambda\right )   \mathcal{h}\ , , \ , \qquad \frac{d\mathcal{v}}{dx }     = \lambda\mathcal{v } -\mathcal{h}\ , ,   \label{s7e6}\ ] ] where the operator @xmath199   $ ] is as in ( [ s3e13 ] ) . using the approximation ( [ e1 ] ) that will be seen to be still valid in this region we obtain @xmath200   \left (   x\right )   = \left ( 1+o\left (   \lambda\right )   \right )   \int_{-\infty}^{x}e^{\left (   y - x\right ) } \mathcal{h}\left (   y\right )   dy\int_{\log\left (   e^{x}-e^{y}\right )   } ^{x}\mathcal{h}\left (   z\right )   dz\,.\ ] ]    taking the limit @xmath201 in ( [ s7e4])-([s7e6 ] )
we obtain @xmath202      if @xmath191 but @xmath186 ( and @xmath192 ) remains bounded we can approximate ( [ s3e10])-([s3e13 ] ) as follows . keeping just the leading terms in ( [ s3e10 ] ) , ( [ s3e13 ] ) we arrive at the approximation @xmath203 . plugging this approximation into ( [ s3e12 ] ) we obtain @xmath204
we are going to compute the solutions of ( [ s7e7 ] ) , ( [ s7e9 ] ) explicitly using laplace transforms . in order to bring ( [ s7e7 ] ) , ( [ s7e9 ] ) into the form of a convolution equation we introduce the set of variables @xmath205 that transforms ( [ s7e7 ] ) , ( [ s7e9 ] ) into @xmath206    notice that this system of equations can be obtained from ( [ s3e4 ] ) , ( [ s3e7 ] ) by means of the rescaling @xmath207 and taking the limit @xmath208 we can compute explicitly a family of solutions of ( [ s8e1 ] ) , ( [ s8e2 ] ) that will be used to describe the solutions of ( [ s3e10])-([s3e13 ] ) in the limit @xmath209    [ solint ]    for any given @xmath210 and @xmath211 there exists a solution to ( [ s8e1 ] ) , ( [ s8e2 ] ) that satifies @xmath212 it is given by @xmath213 \frac{e^{\zeta}}{\zeta}d\zeta\label{f7e5}\ ] ] where @xmath214 is a contour in the complex plane contained in the domain @xmath215 satisfying @xmath216 as @xmath217 and @xmath218 as @xmath219 with @xmath220    if @xmath221 we have @xmath222 and @xmath223 for any @xmath224    the following asymptotics hold : @xmath225@xmath226@xmath227    if @xmath228@xmath229    the solutions with different values of @xmath230 are essentially the same up to rescaling
. notice also that ( [ f7e7 ] ) implies that @xmath231 becomes negative for @xmath8 sufficiently large if @xmath232    some comments about the choice of the range of values of @xmath233 are in order .
we need the restriction @xmath234 to avoid the singular points @xmath235 crossing the imaginary axis .
it would also be possible to include the value @xmath236 in the second formula of ( [ g1e1 ] ) and ( [ g1e2 ] ) since the right - hand side vanishes .
however this would result in asymptotic formulas like @xmath237 that do not have a precise meaning .
finally the constraint @xmath238 in the second formula of ( [ g1e3 ] ) is strictly needed , since the asymptotics of @xmath239 for @xmath240 is proportional to @xmath241    it is interesting to note that the solutions described in theorem [ solint ] for @xmath242 are just the self - similar solutions with algebraic decay that have been obtained for the coagulation equation with constant kernel in @xcite .
this correspondence can be seen as follows .
the solutions we obtained in theorem [ solint ] satisfy @xmath243 with @xmath244 using this formula to eliminate @xmath245 from ( [ s8e1 ] ) we obtain @xmath246    differentiating ( [ a1 ] ) we obtain the equation that is satisfied by the self - similar solutions of @xmath247   -\frac{\left (   1-a\right ) } { \tau}f\ ] ] having the form @xmath248 with @xmath249   $ ] as in ( [ s1e1])@xmath250 equation ( [ a2 ] ) can be transformed into ( [ s1e1 ] ) using the change of variables @xmath251 and @xmath252 and correspondingly the solutions obtained in theorem [ solint ] are transformed into those obtained in @xcite if @xmath253    integrating the first equation in ( [ s8e2 ] ) we arrive at @xmath254 using the second equation in ( [ s8e2 ] ) we obtain @xmath255 assuming that all the integrals involved are convergent .
therefore , using ( [ s8e2a ] ) and ( [ s8e4 ] ) we find @xmath256    eliminating @xmath257 and @xmath258 from ( [ s8e1 ] ) , the second equation in ( [ s8e2 ] ) , as well as ( [ s8e2a ] ) , ( [ s8e3 ] ) , ( [ s8e4 ] ) , implies @xmath259    integrating by parts in the last equation we obtain @xmath260    equation ( [ s8e6 ] ) is a convolution equation . in order to solve it we introduce the laplace transform of @xmath261 @xmath262
then ( [ s8e6 ] ) becomes @xmath263 .
this equation can be integrated explicitly .
its only solution that does not have singularities and takes real values along the line @xmath264 is @xmath265 \ , .
\label{s8e8}\ ] ] inverting the laplace transform we obtain@xmath266   \frac{e^{\zeta x}}{\zeta}d\zeta\ , .
\label{s8e9}\ ] ]    the change of variables @xmath267 transforms this formula into ( [ f7e5 ] ) .
the choice of the contour of integration ensures the exponential convergence of the integral as well as the absence of singularities of the integrand along the contour of integration . taking the limits @xmath268 and @xmath50 we obtain ( [ f7e7 ] ) .
the negativity of @xmath269 for @xmath221 can be proved by differentiating ( [ s8e9 ] ) and deforming the contour of integration to a double half - line following @xmath270 in a positive and negative direction with @xmath271 and @xmath272 respectively .
then @xmath273    the asymptotics in ( [ g1e1 ] ) can be obtained rewriting ( [ s8e9 ] ) as @xmath274 and@xmath275 differentiating these formulas we obtain @xmath276 and using the change of variables @xmath267 in both integrals it follows that @xmath277    taking the limit @xmath278 in the first formula and @xmath36 in the second and using the fact that @xmath279 we obtain ( [ g1e1 ] ) . using then ( [ s8e2 ] ) we obtain ( [ g1e2 ] ) , ( [ g1e3 ] ) .
formula ( [ g1e4 ] ) follows by integrating by residues .
we can reformulate the results in theorem [ solint ] in terms of the original functions @xmath280 in the following we choose @xmath281 as a suitable normalization .
[ thintas]for any @xmath234 there exists a solution of ( [ s7e7 ] ) , ( [ s7e9 ] ) satisfying @xmath282 .
it is given by @xmath283@xmath284@xmath285    if @xmath221
we have @xmath286 for any @xmath287 .
we have the asymptotics @xmath288@xmath289@xmath290    if @xmath228@xmath291      equation ( [ s7e10 ] ) that approximates the behaviour of the solutions of ( [ s3e10])-([s3e13 ] ) if @xmath191 and @xmath186 is of order one , can be solved explicitly since @xmath292 is preserved along trajectories .
the solutions of ( [ s7e10 ] ) will be used to match the solutions of ( [ s7e7 ] ) , ( [ s7e9 ] ) obtained in the previous section for @xmath293 for this range of values @xmath171 is almost constant if @xmath294 moreover , since we are away from the ode regime described in section [ ode ] we can assume that @xmath295 in ( [ s6e6 ] ) is of order one and therefore @xmath296 is of order one , whence also @xmath297 is of order one .
suppose that we consider a matching region ( to be made precise later ) where @xmath296 is of order one , @xmath171 is almost constant and satisfies @xmath298 .
then the first equation in ( [ s7e10 ] ) implies that @xmath186 decreases exponentially .
more preciesely , let us assume that @xmath299 for @xmath300 then , the asymptotics of @xmath301 are @xmath302 for some @xmath303 .
the asymptotics ( [ s9e2 ] ) will be shown to match with suitable solutions among the ones described in theorem [ solint ] .
moreover , notice that ( [ s7e10 ] ) implies that these asymptotics are valid as long as @xmath186 remains larger than some small number @xmath304 that we can assume to be of order one , although small .
the exponential decay of @xmath186 in ( [ s9e2 ] ) implies that the transition between @xmath299 and @xmath305 takes place on a length scale @xmath306 of order one .
we now describe the solution when @xmath186 is less than @xmath307 . for this purpose
let us define by @xmath308 and @xmath309 the two roots of the equation @xmath310 .
then , in the limit @xmath29 , the curves given by ( [ s9e1 ] ) in the plane @xmath311 are approximately the two horizontal lines @xmath312 and @xmath313 connected by a curve that is approximately vertical . along this connecting curve @xmath171
changes an amount of order one and @xmath314 continues to be of order one since @xmath315 is preserved .
when @xmath186 is less than @xmath304 we can approximate the second equation in ( [ s7e10 ] ) as @xmath316 then @xmath317 increases exponentially .
the first equation in ( [ s7e10 ] ) indicates that for the solutions of ( [ s7e10 ] ) under consideration , @xmath186 remains small while @xmath171 changes from @xmath318 to@xmath319 then we can use the approximation @xmath320 where @xmath321 is the range of @xmath182 that it takes for @xmath186 to decrease from @xmath299 to @xmath322 we recall that @xmath321 is of order one .    notice that the structure of the curve ( [ s9e1 ] ) implies that eventually @xmath186 becomes again of order @xmath304 with @xmath171 close to @xmath323 it then follows from ( [ s7e10 ] ) that for small @xmath324 @xmath186 reaches again the value @xmath172 at @xmath325 after an additional length @xmath326 of order one .
notice that , using ( [ s9e3 ] ) , we have @xmath327 and @xmath186 has the behaviour , for @xmath328 , @xmath329 we need to estimate the relation between @xmath308 and @xmath330 from ( [ s9e1 ] ) we deduce that to leading order @xmath331    on the other hand ( [ s9e4 ] ) allows us to obtain an approximation for the transition length @xmath332 indeed , we have @xmath333 whence , since @xmath334 are of order one and @xmath304 can be made arbitrarily small @xmath335 as @xmath209    we have described the transition of the solution @xmath336 of ( [ s3e10])-([s3e13 ] ) for small @xmath90 and @xmath337 the following result will play a crucial role in describing the solutions , because it will show that the amplitude of the oscillations increases with @xmath187    [ l1]suppose that @xmath338 @xmath339 satisfy ( [ s9e6 ] )
. then @xmath340 .
this follows from the inequality @xmath341 that is just a consequence of @xmath342 .
we now describe the behaviour of the solutions of ( [ s3e10])-([s3e13 ] ) for @xmath201 and @xmath343 with @xmath61 as in ( [ s6e6 ] ) .
it turns out that these solutions can be described for this range of values by a sequence of intervals where the solutions can be approximated alternately by solutions of ( [ s7e7 ] ) , ( [ s7e9 ] ) or by solutions of ( [ s7e10 ] ) with both types of regions connected with suitable matching regions .
more precisely , suppose that @xmath344 is a point where @xmath345 .
notice that at such a point we can expect to be able to use the approximation ( [ s7e10 ] ) .
this will be seen matching the ode - integrodifferential equation regime with the ode regime described in section [ odelolkavolterra ] .
we will assume for the moment that this approximation is valid .
then @xmath192 is also of order one and we can approximate it to the leading order as @xmath203 ( cf .
section [ ss.ode ] ) .
we can also use the approximation ( [ s9e5 ] ) with @xmath346 and approximate @xmath171 by a constant .
then we obtain for @xmath347 of order one that @xmath348    we can now match the asymptotics ( [ t2e1 ] ) with those obtained for the solutions of ( [ s7e7 ] ) , ( [ s7e9 ] ) .
we use ( [ s7e1 ] ) to obtain the matching condition @xmath349    this behaviour must be matched with the one for the solutions of ( [ s7e7 ] ) , ( [ s7e9 ] ) obtained in theorem [ thintas ] . more precisely we will match these behaviours with the ones of @xmath350 for a suitable choice of @xmath351 and @xmath352 the asymptotics of @xmath353 combined with the one in ( [ t1e3 ] ) yield @xmath354 since in the matching region we expect to have @xmath355 on the other hand ( [ t1e5 ] ) , ( [ t1e6 ] ) give the choice @xmath356 with these choices of @xmath357 and @xmath358 we obtain , using ( [ t1e3 ] ) , ( [ t1e5 ] ) , ( [ t1e6 ] ) and ( [ t2e2 ] ) , a matching to the first order between ( [ t2e1 ] ) and @xmath359 in the common region of validity where @xmath360 and @xmath361 .
we can then use the approximation @xmath359 to describe the solutions of ( [ s3e10])-([s3e13 ] ) for small @xmath90 and @xmath362 of order one .
this approximation breaks down for @xmath363 sufficiently large .
indeed , ( [ t1e5 ] ) and ( [ t1e6 ] ) imply that @xmath364 and @xmath365 converge exponentially to zero as @xmath366 however , the approximation ( [ s7e7 ] ) , ( [ s7e9 ] ) is only valid if @xmath367 in the region where @xmath368 becomes of order @xmath90 we must use again the approximation ( [ s7e10 ] ) . in particular this region can be described using the analysis in section [ ss.odeanalysis ] .
the asymptotics of the solutions for @xmath186 of order one is as in ( [ s9e2 ] ) for suitable choices of @xmath369 more precisely , since this matching is made for @xmath370 we can use the asymptotics ( [ t1e3 ] ) , ( [ t1e5 ] ) , ( [ t1e6 ] ) for @xmath371 we then obtain , using also the rescaling ( [ s7e1 ] ) , the matching conditions @xmath372 we now match these formulas with ( [ s9e2 ] ) ( combined with the approximation @xmath203 ) . to this end
we must choose @xmath373    we remark that all this analysis will be meaningful only if @xmath374 therefore @xmath375    the region where @xmath376 can then be described using the analysis in section [ ss.odeanalysis ] .
the conclusion of this analysis is that @xmath377 moves close to the point @xmath378 for some suitable @xmath379 in a characteristic length given by ( [ s9e7 ] ) . more precisely if we define ( cf .
( [ s9e6 ] ) ) @xmath380 @xmath381 we can approximate @xmath301 by means of ( [ t2e1 ] ) with @xmath344 replaced by @xmath382    combining ( [ t2e3 ] ) , ( [ t2e4 ] ) and ( [ t2e6 ] ) we obtain to the leading order@xmath383      we first remark that equation ( [ s5e6 ] ) agrees with the approximating equations ( [ s7e7 ] ) , ( [ s7e9 ] ) and ( [ s7e10 ] ) if @xmath384 in their respective regimes of validity .    indeed , we can rewrite ( [ s5e6 ] ) using ( [ s5e5 ] ) as @xmath385 and these equations yield the same behaviour as ( [ s7e10 ] ) if @xmath186 is bounded .    on the other hand , in order to approximate the solutions of ( [ s7e7 ] ) , ( [ s7e9 ] ) for @xmath384 we use the fact that for small @xmath386 the solutions of ( [ s7e7 ] ) , ( [ s7e9 ] ) in theorem [ thintas ] have the characteristic length scale @xmath387 then , for @xmath233 small , which corresponds to the matching region indicated above , we have a characteristic length very large compared with one .
then , as in ( [ s5e2 ] ) we can approximate @xmath388 \sim
\lambda \mathcal{h } \mathcal{u}$ ] .
since @xmath389 we have @xmath390 and we obtain from ( [ s7e7 ] ) the approximation @xmath391 and plugging this formula into ( [ s7e9 ] ) and using the rescaling ( [ s7e1 ] ) we obtain @xmath392 that agrees with ( [ t2e7 ] ) in the range of values @xmath393    the previous computations show that the equations used in both transition regimes agree in the intermediate matching regime . actually it is possible to check that the main characteristics of the computed solution also agree , as could be expected .
notice that the length scale @xmath394 can be computed using the function @xmath138 in ( [ s6e5a ] ) for @xmath384 : @xmath395    in the computation of this integral , we have split the region of integration into the two subsets @xmath396 and @xmath397 where @xmath398 is a small but fixed number .
it turns out that the contribution of the second integral is much smaller than the one due to the first integrand , exactly in the same manner that the contributions in @xmath399 computed in subsection [ iteration ] are due mostly to the region with small values of @xmath400 using the rescaling ( [ s5e5 ] ) we then obtain the approximation @xmath401    on the other hand , we can compute the same length using ( [ t2e6a ] ) . taking into account that in the matching region @xmath402 are small
, we obtain from ( [ t2e5 ] ) that @xmath403 using ( [ t2e6a ] ) we then obtain to leading order that @xmath404 next , equation ( [ s6e6 ] ) and the definition of @xmath357 by means of @xmath405 imply @xmath406 and plugging ( [ t3e4a ] ) into ( [ t3e4 ] ) we obtain the sought - for matching with ( [ t3e2 ] ) .
finally we obtain the matching for the formula of the change of energy in each iteration .
we have obtained using the ode approximation that the change of energy in each cycle is given by ( [ s6e7 ] ) where in the matching region @xmath407 we can use the approximation ( [ s6e8a ] ) .
then @xmath408 on the other hand , using again @xmath409 and plugging it in ( [ t3e3 ] ) we obtain again ( [ t3e5 ] ) and this yields the desired matching .
the structure of the solutions of ( [ s3e10])-([s3e13 ] ) described in section [ asymptpeaks ] by means of alternate approximate solutions of ( [ s7e7])-([s7e9 ] ) and ( [ s7e10 ] ) yields several sequences of numbers @xmath410 that measure the amplitude of the oscillations , the positions of the points where @xmath411 and the position of the peaks where @xmath186 and @xmath192 are of order @xmath190 respectively .
it is relevant to notice that the sequence @xmath412 is increasing due to lemma [ l1 ] .
the values of @xmath357 approach zero in the matching region ( see the approximation ( [ t3e4a ] ) for @xmath407 ) , but they become of order one for @xmath413 with increases in each iteration of order one . eventually , this sequence of points becomes larger than or equal to one .
a detailed study of the solutions of ( [ s3e10])-([s3e13 ] ) shows that there are several possibilities . if @xmath357 becomes strictly larger than one , theorem [ thintas ] shows that during the peak regime @xmath171 becomes strictly negative .
in particular @xmath171 reaches the value @xmath414 at some finite @xmath415 if @xmath171 remains positive for larger values of @xmath182 , something that could happen if @xmath416 it would be possible to approximate ( [ s3e10])-([s3e13 ] ) by means of ( [ s7e10 ] ) . therefore @xmath171 would increase , and the value of @xmath357 would be increased to a new value @xmath417 at the next iteration .
the only possibility of not having such a behaviour would be with @xmath418 as @xmath419 such approximation to zero can not be described by the approximate equation ( [ s7e10 ] ) , because in such a regime the integral terms in ( [ s3e10])-([s3e13 ] ) must be taken into account in full detail .
we now recall that for any @xmath420 there exists a unique solution of ( [ s3e10])-([s3e13 ] ) satisfying ( [ s4e0 ] ) .
the discussion in section [ ss.variables ] shows that all the solutions of ( [ s3e10])-([s3e13 ] ) with such a behaviour can be obtained , up to rescaling , by choosing @xmath59 in the interval @xmath421 a continuity argument taking into account the behaviour of the trajectories shows that there exists a value of @xmath422 in this interval for which the corresponding solution of ( [ s3e10])-([s3e13 ] ) satisfies @xmath423 this also implies @xmath424 and @xmath425 using ( [ s3e10 ] ) , ( [ s3e12 ] ) .
more precisely , the inversion of ( [ s3e10 ] ) allows us to write @xmath426 in terms of @xmath427 then @xmath428 and ( [ s3e12 ] ) yields @xmath429
we finally describe the asymptotics of the solution of ( [ s3e10])-([s3e13 ] ) satisfying ( [ s4e0 ] ) with @xmath422 described in the previous section .
we have seen that @xmath431 to leading order the behaviour of this trajectory is described by the solution of ( [ s7e7 ] ) , ( [ s7e9 ] ) given in theorem [ thintas ] with @xmath432 the asymptotics @xmath433 is valid as long as @xmath434 , where , by assumption , to the leading order @xmath435 it is convenient to reformulate ( [ t4e1 ] ) in the original variables ( cf .
( [ s3e8 ] ) , ( [ s3e9 ] ) ) @xmath436 that is valid for @xmath437    we can match ( [ t4e2 ] ) with the following exponential asymptotics for the solutions of ( [ s2e5 ] ) as @xmath50 : @xmath438    the asymptotics ( [ t4e2 ] ) , ( [ t4e3 ] ) match for small @xmath90 if @xmath439 this gives the desired exponential decay for the obtained solution .
we now show that the solution we obtained is self - consistent with the approximations made for the integral operator @xmath107   \left ( x\right )   $ ] defined in ( [ s3e13 ] ) .
we have made two main approximations .
the first one is ( [ e1 ] ) , that approximates @xmath107   $ ] by an operator with a simple dependence on @xmath82 the second approximation is ( [ s5e2 ] ) and it allows us to replace the integral operator by a much simpler local operator .    concerning the validity of ( [ e1 ] ) we remark that its precise meaning is that the error made in the approximation is smaller than the right - hand side . in order to check its validity
we distinguish between the regions where @xmath192 is of order @xmath190 and the regions where @xmath192 is smaller than that quantity .
notice that for any region we can expect the errors made in the approximation to be of order @xmath440 and we can expect these two terms to be very small compared with the right - hand side of ( [ e1 ] ) .
indeed , in the case of @xmath441 , the contribution due to the region @xmath442 is small compared with the term in ( [ e1 ] ) if @xmath304 is small . on the other hand , if @xmath443 we can use the smallness of the exponential factor @xmath444 to obtain estimates for the corresponding terms in @xmath441 as @xmath445 and since @xmath446
it then follows that this contribution is exponentially small . concerning @xmath447
we can argue similarly for @xmath448 smaller than @xmath304 . if @xmath449 , since @xmath450 it follows that the size of the region of integration can be bounded by @xmath451 and therefore , it gives also a very small contribution .    concerning the approximations of @xmath107   $ ] by means of local terms that have been made in the derivation of the ode approximations ( [ s5e6 ] ) or ( [ s7e10 ] )
two main ingredients are required . in the derivation of ( [ s5e6 ] )
we have used just the fact that @xmath452 has significant changes over distances much longer than one . since
the order of magnitude of @xmath192 is roughly the same for the range of values described by means of ( [ s5e6 ] ) no special care is required to estimate the values of @xmath452 with @xmath453 on the other hand , in the case of the approximation ( [ s7e10 ] ) we use the fact that @xmath107   $ ] is a corrective term that is completely ignored in ( [ s7e10 ] ) .
notice that @xmath452 takes values for @xmath454 much larger than the ones of @xmath426 in the region where the approximation ( [ s7e10 ] ) is used .
however , due to the exponential factor @xmath444 as well as the exponential decay of @xmath452 in the region described by the integro - differential equation ( cf .
subsection [ iteration ] ) the corresponding contribution of such values of @xmath452 is negligible .
during the part of the dynamics dominated by a perturbation of a conservative ode it has been assumed that the energy @xmath61 does not change significantly during each cycle of length @xmath158 this is the justification for the adiabatic approximation and we will check now that , indeed , the adiabatic approximation yields such smallness for the variation of the energy during the trajectory and thus the assumption is self - consistent .
let us denote by @xmath455 a value where @xmath456 our goal is to estimate this quantity for @xmath457where @xmath458 is the next value of @xmath14 where @xmath459 and @xmath460 again .
due to the invariance of ( [ s5e6 ] ) we can assume without loss of generality that @xmath461 due to ( [ s6e5 ] ) and assuming the adiabatic approximation we obtain the following approximation for the change of the energy along a trajectory : @xmath462   ds=\sigma\left (   \xi\right)\ , .   \label{a1e1}\ ] ] since @xmath463\ ] ] and given the form of the curves @xmath464 satisfying ( [ s6e4 ] ) we obtain the existence of numbers @xmath465 satisfying @xmath466 as well as @xmath467 moreover @xmath468 , @xmath469 and @xmath470    the value of the total change of the energy during each cycle , namely @xmath471 has been computed in ( [ s6e7 ] ) , ( [ s6e8 ] ) and approximated asymptotically as @xmath151 ( cf . ( [ s6e8a ] ) ) .
notice that @xmath472 , assuming that @xmath473    due to ( [ a1e2 ] ) , ( [ a1e3 ] ) we need to compute @xmath474 in order to estimate the range of variation of @xmath475 .
we then have @xmath476 where we use that @xmath477 in order to estimate these quantities we define an auxiliary quantity @xmath478 via @xmath479 that measures the maximum value of @xmath480 for a given value of @xmath61 .
it will be more convenient to compute @xmath481 using @xmath480 as independent variable .
to this end we define @xmath482 by means of the roots of the equation @xmath483    taking into account ( [ s6e3 ] ) we can use @xmath480 as variable of integration instead of @xmath484 in order to compute @xmath485 for @xmath486 then , with @xmath487 we obtain @xmath488 d\omega_{0}\ \ , \ \
k=1,2,3\ , , \label{a1e4}\ ] ] where the integration is made along the curve defined by ( [ u2e3 ] ) .
notice that we use @xmath489    in order to approximate @xmath490 for small @xmath478 we use taylor in ( [ u2e3 ] ) to obtain the approximation @xmath491 that combined with ( [ a1e4 ] ) yields @xmath492   d\omega_{0}\ \mbox { as } \omega_{0}\rightarrow0 \,.\ ] ] on the other hand , we can approximate @xmath493 by @xmath494 for @xmath495 whence @xmath496 as @xmath497 for @xmath495 and @xmath498          since the region where @xmath179 is of order one gives a small contribution to the integrals , we obtain the approximations @xmath505 combining this with ( [ s6e8a ] ) we obtain that @xmath506 for some @xmath507 independent of @xmath502 it then follows that the change of energy during each cycle can be estimated by the final change . therefore
, this justifies the adiabatic approximation .
* acknowledgment : * bn and jjlv gratefully acknowledge the warm atmosphere at the isaac newton institute for mathematical sciences where part of this work was done during the program on pde in kinetic theories .
this work was also supported by the epsrc science and innovation award to the oxford centre for nonlinear pde ( ep / e035027/1 ) and through the dges grant mtm2007 - 61755 .
r.  l. drake , _ a general mathematical survey of the coagulation equation _ , in : topics in current aerosol research ( part 2 ) .
international reviews in aerosol physics and chemistry , oxford .
pergamon press ( 1972 ) , 203 - 376 .
p. laurenot and s. mischler , _ on coalescence equations and related models _ ,
modeling and computational methods for kinetic equations , eds .
p. degond , l. pareschi , g. russo , series modeling and simulation in science , engineering and technology , birkhuser , ( 2004 ) , 321 - 356 .