the theory of inflation provides an explanation for many observed properties of our visible universe , in particular its homogeneity and isotropy ( see e.g.  @xcite and  @xcite for reviews ) .
it also gives an explanation for the origin of small density perturbations which eventually gave rise to the structures in our universe .
these aspects are well understood . however , from the particle physics point of view , a suitable candidate for the inflaton field has yet to be found  @xcite . among the many candidates , hybrid inflation  @xcite is a particularly well motivated one , especially in supersymmetric models  @xcite@xcite .
there has been a lot of recent work investigating the possibility of realising inflation within string theory ( see e.g.  @xcite@xcite for reviews ) .
the most natural framework for brane inflation is the brane - antibrane system where an attractive force leads to a potential for the interbrane distance  @xcite@xcite .
inflation ends with an open string tachyonic instability similar to hybrid inflation .
this scenario has been extensively studied and leads to interesting phenomena such as @xmath4-string formation at the end of inflation .
one particular issue in these models comes from the need for moduli stabilisation during inflation .
if this is not achieved , the potential for the moduli is of runaway type destroying the existence of slow roll inflation .
this problem can be circumvented using the kklt scenario whereby the kahler moduli are stabilised once the complex moduli have been stabilised and non - perturbative gaugino condensation has occurred on d7 branes  @xcite . the coupling to inflation has been studied in kklmmt where a fine - tuning of the inflaton superpotential has been advocated  @xcite .
one can also use a shift symmetry to alleviate the @xmath2 problem  @xcite@xcite .
another way of realising hybrid inflation in string models can be obtained using the d3/d7 system , see e.g.  @xcite@xcite . in this case
the end of inflation happens at a lower scale than the planck scale , when the charged open string fields between the branes condense .
the d3/d7 system can be modelled using an @xmath5 supersymmetric @xmath4-term inflation model  @xcite .
the interbrane distance plays the role of the inflaton and the charged open strings between the branes are the waterfall fields .
when the supergravity corrections are neglected , i.e.  in global supersymmetry , and assuming that the kahler moduli have been stabilised , the inflaton direction is flat .
slow roll inflation is driven by the one loop logarithmic corrections as the one loop quadratic divergences are inflaton - independent .    in this paper
, we are investigating @xmath0-term hybrid inflation with the inflationary superpotential @xmath6 where @xmath7 is the @xmath8 gauge coupling , in the presence of additional moduli fields .
it would be very interesting to derive the origin of @xmath9 from fundamental string theory  @xcite .
we treat the supergravity case going beyond the global supersymmetry analysis .
we show that the treatment of fayet - iliopoulos @xmath4-terms ( i.e.  when @xmath10 ) needs extra care in supergravity . to include such a term would require additional fields and the study of their stabilisation . for simplicity we therefore do not include inflationary @xmath4-terms in our models .
inflation is then driven by the presence of the non - vanishing @xmath9 .
we show that supergravity corrections , and supergravity - induced interactions with a moduli sector induce a tree level slope and mass for the inflaton  @xcite .
moreover there are additional contributions from the quadratic divergence part of the one - loop corrections , since they are inflaton - dependent in supergravity , as opposed to inflaton - independent in global supersymmetry . in general
these effects spoil the flatness of the inflationary potential . we show that most of the tree - level problems can be removed by including shift symmetries in the inflaton sector , and taking the moduli sector to be no - scale .
it then appears that a viable slow - roll inflation model can be found with limited fine - tuning of the moduli sector .
however we also find that the coupling between the inflaton and the moduli has the effect of inducing a small variation of the moduli during inflation . despite the small size of this variation
, it gives a significant contribution to the inflaton slope and violates slow - roll inflation .
the paper is arranged as follows . in section  [ sec : v ] , we construct the full potential for our combined model of inflation and moduli stabilisation .
we calculate tree and one loop level parts of the potential , and describe symmetries which help maintain the flatness of the inflaton potential . in section  [ sec : mod ] we determine the role of moduli fields during inflation , and show their effect is significant .
we find the constraints corresponding to the cobe normalisation and the wmap3 spectral index constraint in section  [ sec : con ] .
our arguments and analysis apply to a general moduli sector , although as a specific example , we also give explicit results for the kklt model , which has just one moduli field @xmath11 . in section
[ sec : kklt ] we study particular choices of racetrack superpotentials for the kklt scenario , and find that the running of the moduli field during inflation typically leads to a large value of @xmath2 .
a small value of @xmath2 could only be achieved if the moduli superpotential is heavily fine - tuned .
finally in section [ sec : nogo ] we consider models where a minkowski vacuum after inflation is obtained without the need for a lifting term . a model with this property
is also interesting as it may give a small gravitino mass , as is required to have an low energy sparticle spectrum .
we derive a no - go theorem which states that , for no - scale models with no lifting term , no such stable supersymmetry breaking minkowski vacua can be found .
we conclude in section  [ sec : conc ] .
as already discussed , many string theory inflation scenarios can be modelled at low energy with a supergravity description of hybrid inflation .
the system uses three fields , the inflaton @xmath12 measuring the inter - brane distance and two charged fields @xmath13 ( which in the d3/d7 system would represent the open strings between the two types of branes ) .
the fields interact according to the superpotential  ( [ inflasec ] ) with @xmath10 , which has been studied extensively in the literature . in  @xcite it was shown how the fayet - iliopoulos @xmath4-term arose from fluxes on the d7 brane , allowing @xmath4-term inflation in this system . here
we will focus almost exclusively on hybrid inflation which is driven only by @xmath0-terms ( so @xmath14 ) , although exactly how such terms are embedded in string theory remains to be found . including inflationary @xmath4-terms in a supergravity theory
creates extra complications , as we will discuss in the following subsection .      using  ( [ inflasec ] ) and taking to @xmath15
we obtain @xmath16 this potential has an @xmath17 origin which was derived in  @xcite . during inflation , when the waterfall fields @xmath18 , the potential is @xmath19 which gives a positive contribution to the inflation energy .    extending this model to a supergravity theory , and working in units with @xmath20
, we replace the above @xmath0-term by @xmath21 where @xmath22 , and @xmath23 is the kahler potential for @xmath12 and @xmath13 .
the simplest possibility would be @xmath24 .
it was originally thought that extending the @xmath4-term in the above model to sugra would be straightforward .
however it was later discovered that in order to have a fayet - iliopoulos term @xmath25 appearing in a supergravity theory , the superpotential must have charge @xmath26  @xcite .
thus the extension of  ( [ susyinf ] ) to sugra described in the earlier paper  @xcite requires reconsideration . to preserve the properties of @xmath12 and its role as the inflaton
, it should be uncharged  @xcite .
this is only possible when @xmath10 and the charges of @xmath13 are shifted to @xmath27 , where @xmath28 .
now all the contributions to @xmath29 from any other sectors of the theory would also have to be charged .
this would be for instance the case with the mssm superpotential .
since they are , _ a priori _ , unrelated to inflation , this seems unnatural .
hence we have an incompatibility between a non - vanishing @xmath9 and a fayet - iliopoulos term @xmath25 . in the usual inflationary scenario ,
inflation is driven by the fayet - iliopoulos term .
one possibility is to introduce a fayet - iliopoulos term with a gauge invariant superpotential .
this can be achieved in supergravity with a @xmath4-term by making the @xmath30 pseudo - anomalous .
we extend the model to include one modulus @xmath31 and matter fields @xmath32 , all of which are charged under the pseudo - anomalous @xmath30 gauge group .
the relevant part of the kahler potential for these fields is @xmath33 where the gauge field , @xmath34 has been included ( the usual kahler potential appearing in the scalar potential expression is obtained by putting @xmath35 formally ) .    the potential arising from the @xmath4-term is then @xmath36 where we have taken the gauge coupling function to be @xmath37 in order to cancel the @xmath30 anomaly .
we see that a fayet - iliopoulos term can give rise to inflation when @xmath10 , assuming that @xmath31 and @xmath32 can be stabilised by a suitable superpotential .    during inflation ,
the waterfall fields @xmath13 vanish .
this implies that the @xmath30 @xmath4-term and the @xmath9 dependent part of the @xmath0-term potential give rise to the inflation energy @xmath38 when @xmath10 and the extra fields are stable , this resembles the usual inflation energy in hybrid inflation .
if on the other hand we just use the @xmath0-term to give inflation , there is no need to introduce the additional fields @xmath32 and @xmath31 .
the fayet - iliopoulos @xmath4-term must then be zero , due to the gauge invariance of @xmath29 , so @xmath39 .
we will take @xmath40 and assume this is the case for the rest of the paper .
the string theory origin for @xmath9 has yet to be found  @xcite .
we are interested in the hybrid inflation scenario based on the superpotential above  ( [ inflasec ] ) .
we are also taking into account the effects of a separate sector , so the total superpotential is @xmath41 the @xmath42 are moduli and matter fields , which are usually assumed to be stable during inflation . with one modulus in @xmath0-theory
no stabilisation was obtained  @xcite .
here we consider the most general case in supergravity with more than one moduli field and generic superpotential @xmath43 .
the moduli dependent part of the superpotential @xmath43 is then responsible for stabilising the moduli fields .
the total kahler potential has the form @xmath44 with @xmath45 we have imposed two shift symmetries in the @xmath46 sector to alleviate the @xmath2 problem of supergravity inflation .
this is an extension to the usual @xmath12 shift symmetry  @xcite@xcite . in string theory
, this corresponds to the translational invariance of the brane system . for the full superpotential  ( [ wfull ] ) ,
the @xmath0-term part of the potential can be expanded as @xmath47 \
, , \end{aligned}\ ] ] where @xmath48 denotes the real part of @xmath49 , and we have defined @xmath50 there is also a contribution to @xmath51 from the @xmath4-term @xmath52}\left(|\phi^+|^2 - |\phi^-|^2\right)^2 \ , .\ ] ]    in this paper we will be particularly interested in kahler potentials with the no - scale property , such as @xmath53 in the kklt scenario . in this case
the above expressions reduce to @xmath54    at the end of inflation , the charged fields @xmath13 condense and the inflaton @xmath12 is zero .
the @xmath30 @xmath4-term part of the potential is then zero , and @xmath0-term part reduces to @xmath55 .
this part of @xmath56 , which comes from the moduli dependent part of the superpotential @xmath43 , is then responsible for stabilising the moduli fields . for @xmath57
its minimum is an @xmath58 vacuum where supersymmetry is preserved , i.e.  @xmath59 for all the fields .
since we want a minkowski vacuum at the end of inflation , it is necessary to add an additional lifting term , @xmath60 to the potential . in the kklt scenario a non - supersymmetric lifting term
is used @xmath61 where @xmath62 depending on the origin of the term , i.e.  @xmath63 for anti d3-branes and @xmath64 for fluxes on @xmath65 branes .
an alternative possibility is that the minimum is raised by an uplifting @xmath4-term  @xcite .
this allows one to lift the @xmath58 vacuum while avoiding the need for a non - supersymmetric potential . in this case
the set of fields @xmath42 comprises one modulus @xmath11 and matter fields @xmath66 , all of which are charged under a second , pseudo - anomalous @xmath67 gauge group .
the kahler potential in the moduli sector can be written as @xmath68 unlike the kklt case it is not no - scale , although it is straightforward to write down a no - scale version .
the potential arising from the @xmath4-term for the additional @xmath67 is then @xmath69 if @xmath70 and all @xmath71 are positive , this potential will be nonzero and provide a suitable lifting term . as explained in  @xcite , an @xmath58 vacuum in supergravity can not normally be lifted by @xmath4-terms .
one can circumvent this argument by using non - analytic superpotentials with positive definite @xmath4-terms , as described in  @xcite .
the total moduli stabilisation potential is then @xmath72 where @xmath73 is given by , for example , ( [ nslift ] ) or ( [ dlift ] ) .
after inflation , the potential possesses a minkowski vacuum thanks to the presence of the lifting term .    for the rest of the paper we will only concentrate on the case where there are no additional matter fields @xmath66 .
this implies that the gauge invariance of the superpotential can only be maintained if @xmath11 is gauge invariant , and we must use the non - supersymmetric lifting term  ( [ nslift ] ) .
this is the minimal setting to study hybrid inflation coupled to moduli .
the structure of this potential is rather rich due to the possible interplay between the moduli and other fields . during hybrid inflation ,
the fields @xmath13 vanish , and the inflaton field is real , so @xmath74 .
the full ( tree - level ) potential then reads @xmath75 where we have identified @xmath76 and @xmath77 in addition to the terms above , there will be contributions from loop corrections , which we will discuss in the next subsection .
note that the inflaton is much smaller than the planck mass , i.e.@xmath78 is such that @xmath79 . during inflation ,
the minkowski vacuum is lifted by @xmath80 .
we see that in general we will get large , tree - level contributions to the slope and mass of the inflaton  @xcite ( from @xmath81 and @xmath82 respectively ) . in other words
the model has @xmath1 and @xmath2 problems .
this new @xmath2 problem comes from the inclusion of the moduli sector , and is distinct from the usual @xmath2 problem which comes from embedding inflation in supergravity . in the case of no - scale models
@xmath83 , and so this @xmath2 problem can be removed by a symmetry of the kahler potential @xmath84 . to obtain a small inflaton slope
, we will need to fine tune @xmath81 .
the form of the potential discussed so far was based on the hybrid inflationary epoch , in which the fields @xmath85 vanish .
inflation ends when the mass matrix acquires a tachyonic direction .
this usually occurs when one of the two charged fields @xmath13 ( or a combination thereof ) has a negative mass , i.e.  when @xmath86 . without moduli fields
, this instability occurs at @xmath87 .
however , the inclusion of the moduli sector modifies this result . for simplicity
we take @xmath83 in what follows ( as is the case for a no - scale @xmath84 ) .
working out the mass matrix for the @xmath13 sector , we find @xmath88 where we have assumed @xmath29 to be real .
if @xmath89 , we recover the standard result @xmath87
. however , due to the presence of the moduli fields , the value of @xmath90 is altered .
note that if we had used a more general @xmath23 , which did not have the shift symmetries , then @xmath90 would have received far more corrections .
taking a canonical kahler potential , and neglecting the contributions of order @xmath81 in the regime @xmath91 , the lowest mass is given by @xmath92\ , .\ ] ] we see that if @xmath93 , inflation never ends as the mass is always positive .
we naturally expect that @xmath29 will be much larger than the inflation scale @xmath9 in order to stabilise moduli during inflation .
we see that both shift symmetries must be included in @xmath94 for inflation to be viable .    in the above discussion
we have considered how the moduli sector will alter inflation .
of course we must not forget that the inflationary sector can also interfere with the moduli stabilisation , since the inflation potential @xmath80 acts as a perturbation of the moduli potential @xmath95 . in the kklt scenario
the minimum of @xmath96 is separated from @xmath97 ( which also has @xmath98 ) by a potential barrier . denoting the height of the barrier by @xmath99
, we see that the value of the above potential  ( [ inflapotential ] ) must be less than @xmath99 during inflation , or the modulus field @xmath11 will roll off to infinity , and inflation will not finish in the correct vacuum .
the slope terms in the above potential must be small ( in order to have slow - roll inflation ) , so we need @xmath100 .
this guarantees that the minimum of the moduli potential is only shifted during inflation and not completely destroyed .
typically @xmath101 , where @xmath102 is the energy of the ads minimum of @xmath51 ( if @xmath73 were absent ) .
this provides a rough estimate of the barrier height . in this case
we need @xmath103 , which is reasonable if the inflation scale is far below the moduli stabilisation scale .
the one - loop corrections to the effective potential for a general theory , with cut - off @xmath104 , are  @xcite @xmath105 the supertrace is @xmath106 .
expressions for the boson and fermion mass matrices are given in the appendix .
as long as the masses are below the cut - off scale , the log corrections are smaller than the quadratic divergences .
this is guaranteed provided the potential @xmath43 is lower than the cut - off , i.e.  lower than the planck scale in practice .
this is also a phenomenological requirement in order to obtain a hierarchy between the planck mass and the gravitino mass .
we therefore concentrate on the quadratic divergences in what follows .
we see there are corrections to all parts of the inflationary potential  ( [ inflapotential ] ) .
in particular we find @xmath107\right ) e^k \ , .\ ] ] in general , the corrections to the other parts of the potential are complicated , and not particularly illuminating . for simplicity
we will explicitly consider the case for one modulus field only , i.e.@xmath108 .
we will use @xmath109 . in the minkowski background @xmath110 during inflation
we have instead @xmath111 we see that the loop corrections are @xmath12 dependent .
note that we have included not only the contributions from @xmath13 , but also those from the moduli sector .    for the kklt scenario
, we find the tree - level potential to be ( after taking @xmath11 to be real ) @xmath112 where @xmath113 .
adding on the loop corrections  ( [ strmink ] ) gives the full potential for the minkowski background @xmath114 - \frac{w w'}{2t^2}\left[1+\frac{\lambda^2}{8 \pi^2}\right ] -\frac{w^2}{t^3 } \frac{\lambda^2}{32 \pi^2 } \\ \nonumber & & { } + \frac{e}{(2t)^n}\left[1+\frac{\lambda^2 n(n+1)}{48 \pi^2}\right ] \ , .
\label{kkltvs}\end{aligned}\ ] ] adding the loop corrections will change the value of @xmath11 at the minimum .
furthermore the minimum will no longer be minkowski .
the parameter @xmath115 needs to be re - tuned to fix this .    during inflation , i.e.  @xmath12 non - vanishing ,
we get @xmath116 to which we need to add the loop corrections  ( [ strinf ] ) .    in this paper , we take the cut off at the planck scale . hence ,
as long as @xmath117 and @xmath29 do not vanish during inflation , the radiative corrections induce a slope for the inflationary potential , which is dependent on @xmath29 , and which can in principle be a danger to the flatness of the inflationary potential .
this effect is only present in supergravity , as in global supersymmetry the quadratic divergences are @xmath12 independent . for one modulus field @xmath11 ,
the full potential during inflation now reads @xmath118 for @xmath11 real , the full moduli stabilisation potential is given above  ( [ kkltvs ] ) .
we denote the value of @xmath11 at the minimum of @xmath119 by @xmath120 .
similarly @xmath121 - w \frac{\lambda^2}{4 \pi^2 } \label{kkltv1}\ ] ] @xmath122 as it is a no - scale model , the only contribution to @xmath123 comes from loop corrections .
in this section we will discuss the form of the effective inflationary potential , taking into account the moduli stabilisation mechanism .
it is usually assumed that moduli fields are stabilised before inflation begins .
the energy scale of the corresponding potential is much higher than that of inflation , and so the effects of the inflaton will not significantly alter the values of the moduli fields .
it therefore seems reasonable to assume the moduli are fixed during inflation .
as we will show , this is not the case .
let us now analyse the moduli stabilisation mechanism when the inflation sector is present .
we denote the values of the moduli fields at the minimum of @xmath124 by @xmath125 , and define @xmath126 to be the deviation of the moduli fields from this value . expanding the inflationary potential around this minimum
we find @xmath127 where @xmath128 are independent combinations of @xmath129 and @xmath130 which satisfy @xmath131 @xmath132 .
minimising this with respect to @xmath128 , we find @xmath133 so we see that the moduli fields are now stabilised along a valley depending on @xmath12 . if @xmath134 ( which is perfectly reasonable ) , then @xmath135 , and the variation of the moduli fields is tiny , as is usually assumed .
the correction to the inflationary potential is tiny too .
however the slope of @xmath51 is also tiny , and so the contribution of the moduli fields can be comparable to the other effects we are considering . hence it is not reasonable to ignore their variation during inflation after all .
now if we substitute the above expression  ( [ yeq ] ) back into the potential  ( [ vy ] ) , we obtain the effective inflationary potential @xmath136 with slope @xmath137 and effective mass @xmath138 \ , .\ ] ] we remind the reader that loop corrections are included in these expressions .
hence , we see for a wide range of models ( including the kklt scenario ) that @xmath139 , and so the inflaton is tachyonic .
this tachyonic instability is not necessarily a problem if @xmath12 is small , or the mass term is subdominant .
since we need the field to roll towards zero during inflation we also require that @xmath140 is positive .    it is tempting to use a moduli sector whose superpotential satisfies @xmath141 at the minimum , since @xmath142 in this case , and so we automatically avoid any large contributions to the slow - roll parameter @xmath1 .
this type of model has other appealing features , such as a small gravitino mass  @xcite .
if we consider just one moduli field @xmath11 , and take a no - scale @xmath84 , then we see from ( [ kkltv2 ] ) that loop corrections give a small negative inflaton mass .
the effects of the moduli variation will make it even more negative .
even if this mass is small , it will still dominate the behaviour of @xmath12 , and cause it to roll away from zero .
hence we see that a moduli potential with @xmath143 at its minimum ( i.e.  unbroken supersymmetry ) can not be successfully combined with @xmath0-term hybrid inflation .
we have established that during inflation the effective potential has the leading order form @xmath144 with @xmath140 positive and @xmath145 negative . the potential has the form of an inverted parabola .
since we want inflation to end by the tachyonic instability of the @xmath146-field , @xmath12 has to roll towards zero .
this means that during inflation @xmath147 .
thus , an additional constraint on the model parameters is that @xmath148 , as given in  ( [ phiend ] ) . to get an idea about the allowed parameter ranges for the theory , making only modest assumptions about the moduli sector , we provide some general constraints on potentials of the form  ( [ workingmodel ] ) , leaving @xmath149 , @xmath140 and @xmath145 as free parameters , but with @xmath145 negative . note that @xmath149 , @xmath140 and @xmath145 depend on the underlying theory and are _ not _ independent . as shown in the previous section , they depend on the details of the moduli sector as well as the parameters @xmath7 and @xmath9 .
an inflationary model has to satisfy certain constraints in order to agree with observations .
firstly , in order to realise a period of slow - roll inflation , the two slow - roll parameters @xmath150 and @xmath151 have to be small ( we have approximated @xmath152 during inflation ) .    also , in order to get the right amplitude for the density perturbations from inflation , the cobe normalisation must hold  @xcite @xmath153 the current wmap+2dfgrs constraint on the spectral index is  @xcite @xmath154 ( a weaker bound has been found in  @xcite ) . in the case
that @xmath155 , the above two constraints are fulfilled when @xmath156 and @xmath157 the constant @xmath140 can be eliminated from the last expression , using  ( [ norm ] ) .
since @xmath145 is negative and @xmath158 , this restricts the parameters to be in the region ( in planck units ) @xmath159 implying an upper bound for the inflationary scale of order @xmath160 gev .
the allowed region in the @xmath161 plane is shown in figure  1 .
[ ] [ ] [ 2]@xmath149 [ ] [ ] [ 2]@xmath145    the parameters @xmath80 , @xmath140 and @xmath145 depend on the details of the moduli superpotential as well as on @xmath7 and @xmath9 .
in order to generate slow - roll inflation , we require that @xmath162 and @xmath163 note that if we had not included the contribution from the variation of the moduli fields , we would instead have @xmath164 .
in general , the above contributions to @xmath2 from the moduli sector are likely to be large .
we see this type of model may have another @xmath2 problem , which is not removed by the symmetries we used in section  [ sec : v ] to avoid the other two @xmath2 problems .    using  ( [ v0con ] ) , which follows from the cobe and spectral index constraints ( [ cobe ] ) and ( [ nsconstr ] )
, we find that @xmath165 this is consistent with @xmath166 , as we have assumed ( and is required , for example , for @xmath96 to retain its minimum during inflation ) .
finally , for the maximal allowed value of @xmath167 gev , the relative contribution of gravitational waves @xmath168 is 13 percent .
however , it can be easily shown that the field variation @xmath169 during inflation , related to the e - fold number by @xmath170 is larger than @xmath171 ( with @xmath172 ) , as long as @xmath173 gev . since field variations below the planck scale are preferred theoretically , its very likely that the predicted amount of gravitational waves is much lower .      here
we study the stability of the modulus during inflation .
in particular , the previous approximation of assuming that the moduli potential is quadratic when perturbed by the inflaton is only valid if the inflaton does not remove the potential barrier between the inflation minimum and infinity .
we have to ensure that the term linear in @xmath12 in the potential is not too big , otherwise this term will completely destabilise the modulus @xmath11 . the condition for this to happen is that for at least 60 e - folds of inflation the term linear in @xmath12 in  ( [ vinf ] ) is smaller than the barrier height , @xmath99 .
typically @xmath174 which implies @xmath175 using  ( [ nfold ] ) to eliminate @xmath169 from this equation , we get @xmath176 where we have used @xmath177 ( ignoring loop corrections ) , and have dropped numbers of order one . using the cobe normalisation
, we get @xmath178 where @xmath179 .
this is the condition guaranteeing that the modulus is not destabilised during inflation .
it states that the gravitino mass is never small . in terms of the soft breaking terms and
the sparticle spectrum obtained with a large gravitino mass , this implies that the scalar sparticles have very large masses .
this is the type of scalar spectrum advocated in split susy models  @xcite .
of course if one is to observe scalars at the lhc , the scalars can not be very massive and the gravitino mass must be reduced .
this requires a superpotential which gives a small gravitino mass , but which still produces a stabilisation potential with a high barrier .
the model in subsection  [ ssec : eg2 ] has these properties .
it is even possible to have @xmath180 and still stabilise the moduli in a minkowski vacuum , as described in  @xcite . in this case
there is no lifting term @xmath73 , and the vacuum is supersymmetric .
we could then obtain a gravitino mass from additional supersymmetry breaking effects . as we will show in section  [ sec : nogo ] ,
if the vacuum is be minkowski , the susy - breaking must involve a lifting term .
otherwise the vacuum will be unstable .
as has already been noted , the parameters @xmath149 , @xmath140 and @xmath145 are not independent , but depend on details of the moduli sector as well as on @xmath7 and @xmath9 . as a simple example we will discuss the kklt scenario  @xcite , which has just one modulus field @xmath11 and uses a non - supersymmetric lifting term  ( [ nslift ] ) . during inflation the full potential is then given by  ( [ kkltvs ] ) and ( [ vinf])([kkltv2 ] ) .
the stabilisation potential s minimum is at @xmath181 , so @xmath182 .
this allows us to eliminate @xmath115 , and one other parameter .
following  @xcite , we will use a racetrack superpotential @xmath183 where @xmath184 , @xmath185 , @xmath186 and @xmath187 are real constants .
the above superpotential comes from gaugino condensation .
we take @xmath188 and @xmath189 to guarantee the validity of the supergravity approximation and the weak gauge coupling limit .
let us first focus on the kklt model obtained with @xmath190 and @xmath191 . to
leading order this gives @xmath192 and @xmath193\ , , \qquad \tilde v_2 = -\frac{\lambda^2}{4 \pi^2 } \ , , \nonumber \\ & & ( e^k\tilde v_1)'\approx -\frac{a a^2}{4t_0 ^ 2 } e^{-at_0 } \left[1+\frac{\lambda^2}{8\pi^2}\right ]   \ , , \qquad \lambda \approx \frac{a^2 a^4}{3t_0}e^{-2at_0 } \left[1+\frac{\lambda^2}{8\pi^2}\right ] \ , , \label{eg1}\end{aligned}\ ] ] where @xmath194 .    during inflation
we also require the slow roll parameter @xmath195 ^ 2 < 1 $ ] .
we see above that @xmath196 , and so we need @xmath197 .
we observe that @xmath123 is small , suggesting @xmath2 is small too .
however there are also corrections to @xmath2 from the modulus field variation during inflation ( see section  [ sec : mod ] ) , which can not be ignored .
these imply @xmath198')^2}{e^k \lambda}\right ) \ , .
\label{neweta}\ ] ] we see that for the above choice of @xmath0-term model , the slow - roll parameter is found to be negative , independent of the parameters , and large : @xmath199 .
this value is not compatible with the slow - roll condition and , hence , not compatible with observational constraints .
in general there could also be corrections to @xmath2 from variations of the imaginary part of @xmath11 .
however for the examples in this section we find that the potential is always minimised at @xmath200 , even during inflation .
another constraint on the model comes from requiring that the minimum of @xmath201 does not disappear during inflation .
this would be the case if @xmath202 . for the above example
we find that the maximum of the potential is at @xmath203 , with @xmath204/[n+(5n-6)\lambda^2/(24\pi^2)]$ ] .
the barrier height is then @xmath205 \ , .\ ] ] we see that @xmath206 , as we assumed in the previous section .
we need @xmath207 , implying @xmath208 .
however this contradicts @xmath209 , which was required to have small @xmath1 .
hence the moduli superpotential in this subsection fails on both slow - roll parameters .
as a second case , consider the race - track models with @xmath210 , and take the limit @xmath211 . in this case
we find @xmath212 and @xmath213 \ , , \qquad \tilde v_2 \approx -\frac{\lambda^2}{4 \pi^2 } \
, , \nonumber \\ & & ( e^k \tilde v_1 ) ' \approx -a(a - b ) e^{-a t_0 } \frac{a}{4t_0 ^ 2 } \left[1+\frac{\lambda^2}{8\pi^2}\right ]   \ , , \nonumber \\ & & \lambda \approx a^2 ( a - b)^2 e^{-2at_0 } \frac{a^2}{3t_0 } \left[1+\frac{\lambda^2}{8\pi^2}\right ] \ , .
\label{eg2}\end{aligned}\ ] ]    small @xmath1 requires @xmath214 which , as in the previous subsection , implies @xmath215 .
this suggests that it will again be impossible to have small @xmath1 and a stable modulus @xmath11 .    using an asymptotic analysis of @xmath201
, we find that this time the maximum of the potential occurs at @xmath216 , and the height of the barrier is @xmath217 \ , .\ ] ] in contrast to the previous subsection we find @xmath99 is not of order @xmath218 . in this case
it is possible to have @xmath219 and @xmath220 .
unfortunately when we calculate @xmath2 , we run into the same problem as before .
if the modulus field @xmath11 were fixed at @xmath120 , @xmath2 would be small .
however this does not occur , and  ( [ neweta ] ) applies instead .
we find again @xmath199 , ruling out this case too
.      we will now use approximate analytic arguments to determine if any choices of @xmath29 could give viable slow - roll inflation . using the expressions for @xmath201  ( [ kkltvs ] ) , and @xmath221  ( [ kkltv1 ] )
, we obtain @xmath222 v_1 + \frac{(n-4)}{24 t^3 } v_1 ^ 2 + \mathcal{o}\left(\frac{\lambda^2}{8 \pi^2}\right ) = 0 \label{approx}\end{aligned}\ ] ] at @xmath181 .
now the requirement that @xmath223 implies that @xmath224 , and since @xmath166 we can take @xmath225 to find approximate solutions to  ( [ approx ] ) . noting that @xmath226 , the effect of the loop corrections will also be small , allowing us to neglect them . for small @xmath81  ( [ approx ] ) is only satisfied if either @xmath29 or @xmath227 is very small .
for the former we obtain @xmath228 using this , the expression for @xmath229 simplifies to @xmath230 , to leading order in @xmath81 . from  ( [ kkltv1 ] )
we obtain @xmath231 to leading order in @xmath81 this is @xmath232 . substituting all this into the expression for @xmath2  ( [ neweta ] ) , which includes the effect of varying @xmath11
, we find @xmath233')^2}{e^k \lambda } \approx -6 \ , .\ ] ] so we retrieve an @xmath2 problem along this branch of solutions .
the models discussed in the previous two subsections fall into this category .
the second approximate solution to  ( [ approx ] ) has @xmath234 this again allows us to simplify @xmath235 , and we find @xmath236 . expanding in @xmath81 , ( [ dv1app ] )
implies @xmath237 .
combining all this , ( [ neweta ] ) implies @xmath238 noting that @xmath239 in order to get successful inflation , we need both @xmath240 and @xmath241 at the minimum @xmath181 . if this fine - tuning is not satisfied then slow roll inflation does not occur .
this fine - tuning of the moduli superpotential is in addition to the fine - tuning required for the stabilisation potential @xmath96 to have a minkowski minimum .
furthermore , @xmath29 must also be such that the maximum of @xmath96 is greater than the inflationary scale ( so that minimum of the moduli potential does not disappear ) .
this all adds up to a huge amount of fine - tuning , and the two - term racetrack potential used above  ( [ racetwo ] ) does not seem to have enough parameters to satisfy all the conditions simultaneously .
this suggests that a heavily fine - tuned , three - term @xmath29 would be needed to give viable inflation .
it is hard to see how such a fine - tuned model could arise naturally .
in this section we consider minkowski vacua with a non - vanishing gravitino mass .
this is the traditional framework used in particle physics . at the end of inflation ,
the vacuum is a non - supersymmetric configuration with a vanishing cosmological constant .
one would also like to impose @xmath242 in order to have a hierarchy between the planck scale and the sparticle masses .
supersymmetry breaking could be achieved with @xmath0-terms , although as we will show this leads to an unstable vacuum .
it is therefore necessary to include a lifting term in the theory , such as those discussed at the end of subsection  [ ssec : vv ] .
consider a theory with only one modulus field , no lifting term , and the following no - scale kahler potential @xmath243 and @xmath108 .
it is convenient to transform the field according to @xmath244 in terms of the new variable @xmath49 , one obtains @xmath245 and @xmath246 the kahler transformation  ( [ kahlertrafo ] ) preserves @xmath247 , which is the only relevant combination in supergravity . using the new variable @xmath49 one obtains @xmath248 where primes denote here @xmath249 . imposing that the potential has an extremum with a vanishing cosmological constant leads to the following relations @xmath250 and @xmath251 at the extremum when @xmath57 and @xmath252 , i.e.  when supersymmetry is broken at the extremum .
the extremum equations are also satisfied for @xmath253 .
the mass at the extremum involves @xmath254 when @xmath57 and @xmath252 , the extremum equations give @xmath255 at the extremum @xmath256 . as @xmath257 is generically non - zero , this implies that the extremum is a saddle point where the direction along @xmath258 ( for @xmath257 chosen real positive ) is tachyonic .
this is a generic property of no - scale models and implies that the modulus can not be stabilised with a vanishing potential and broken supersymmetry .    on the contrary when supersymmetry is preserved , we find @xmath259 as long as @xmath227 does not vanish , guaranteeing the stability of the model  @xcite .
this argument can be generalised to an arbitrary number of fields in the no - scale case @xmath260 where the @xmath66 s are , for example , matter fields . indeed putting @xmath261 and collectively @xmath262 ,
the kahler superpotential becomes @xmath263 and the superpotential @xmath264 is turned into @xmath265 with these redefinitions , the scalar potential reads @xmath266 the vanishing of the cosmological constant at the extremum implies that @xmath267 the minimum equation reads @xmath268 when supersymmetry is broken . in the supersymmetric case , the minimum equations are @xmath269 .
the stability of the extremum depends on the mass matrix which reads @xmath270 \ , .\end{aligned}\ ] ] using the extremum equation we find that @xmath271 in the supersymmetry breaking case .
this mass matrix has nice algebraic properties .
when @xmath272 has zero eigenvalues , the mass matrix @xmath273 vanishes in these directions . in the directions orthogonal to the zero eigenstates of @xmath272 , all the eigenvalues of @xmath273 are positive but one which vanishes along the eigenvector @xmath274 where @xmath275 .
hence the mass matrix @xmath273 vanishes along @xmath276 and the zero eigenstates of @xmath272 .    in the supersymmetric case
, we find @xmath277 \ , .\ ] ] defining a vector space basis with @xmath278 where @xmath279 , and orthogonal vectors @xmath280 , we can use the sum rule @xmath281 to find that @xmath282 \ , .\end{aligned}\ ] ] which is a positive definite hermitian matrix as @xmath283 for any @xmath284 provided the matrix @xmath272 does not have zero eigenstates corroborating the results in  @xcite .    hence we have found that the mass matrix @xmath285 has zero eigenstates when supersymmetry is broken .
now for a generic supersymmetry breaking model , the holomorphic mass matrix is such that @xmath286 . expanding along the direction @xmath287 where @xmath288 is the minimum leads to a potential @xmath289 notice that this leads to a tachyonic direction as before , implying that the supersymmetry breaking extremum with zero cosmological constant is not stable .
hence we see that a supersymmetry breaking minkowski vacuum from a no - scale theory which does not use a lifting term can not be stable .
since such a setup does not even stabilise the moduli fields , it is certainly not suitable to use in a hybrid inflation model .
in this paper we considered @xmath0-term hybrid inflation in supergravity , when moduli fields are present . in the literature
it is usually assumed that moduli fields are stabilised before the period of slow - roll inflation begins and that they therefore have no impact on the dynamics of the inflaton field .
we have shown that this is not necessarily the case .
firstly , the tree - level coupling of moduli and inflaton fields in supergravity induces a slope and , in general , a mass term in the potential for the inflaton field . on top of this
, loop corrections from both the inflation and moduli sectors of theory contribute to the slope and curvature of the potential .
a further problem comes from the evolution of the moduli fields during inflation .
since they are stabilised by a steep potential , they are roughly constant during inflation . however , despite the small size of their variation , we have shown that they still give a significant contribution to the effective inflaton potential .
this generally leads to a large and tachyonic mass for the inflaton .
we find that this class of model has up to three different @xmath2 problems .
firstly there is the usual one coming from embedding hybrid inflation in supergravity ( which is avoided with the help of shift symmetries ) .
secondly there are tree - level contributions to the inflaton mass from the coupling of the inflaton to the moduli sector ( although these are absent for no - scale models ) .
finally there is the negative contribution to the effective inflaton mass from the rolling of the moduli fields during inflation .
it is not obvious that this third @xmath2 problem can be removed with natural additional symmetries .
the model also has an @xmath1 problem due to the tree - level contributions to the inflaton slope from the moduli sector .
we found that in the popular model for moduli stabilisation , the kklt scenario , the induced mass term for the inflaton field is generally too big to give a period of slow - roll inflation .
it seems that only with a heavily fine - tuned superpotential for the moduli sector may viable slow - roll inflation take place .    to conclude
, we have shown that the coupling of moduli fields to hybrid inflation can have an important impact on the effective inflationary potential .
neglecting this coupling is not consistent and , as we have seen , generally leads to incorrect conclusions about the inflationary dynamics .
these considerations are also relevant to other inflationary scenarios .
in each case it should be checked that the moduli fields , even if thought to be stabilised , do not spoil the flatness of the inflationary potential , and therefore the dynamics of the inflaton field .    we would like to thank k. dasgupta and r. kallosh for comments on a previous version of the manuscript and n. chatillon and u. ellwanger for discussions .
acd thanks cea saclay for hospitality while some of this work was in progress .
acd and cvdb were supported in part by pparc .
scd was supported by the swiss science foundation and the netherlands organisation for scientific research ( nwo ) .
phb acknowledges support from rtn european programme mrn - ct-2004 - 503369 .
we have latexmath:[$m^2_{\rm ( boson ) } = 2 \mathcal{k}^{i\bar \jmath } \partial_i \partial_{\bar \jmath } v + 6f_r^{-1 } \sum_{i=\pm }    assuming the only gauge group in the theory is @xmath30 .
if there were additional gauge groups , the additional gauge sector would also give contributions to the mass matrices .
for the fermions @xmath291 where @xmath292 with @xmath293 , and @xmath294 where @xmath295 , and @xmath296 otherwise .
define @xmath297 \ , .\ ] ] the diagonal entries of the boson mass matrix for the minkowski background are @xmath298 and non - zero elements of fermion mass matrix are @xmath299 hence we obtain @xmath300 for the inflationary background , the diagonal entries of the boson mass matrix are @xmath301e^k \ , , \ \nonumber \\ m^2_{\pm \bar \pm } & = & [ 2(g^2x^4 + g^2[2x^4 + 2x^2 + 1]\phi^2 -\sqrt{2}g \phi \re(w)[1 + 2x^2 ] + |w|^2 ) \nonumber \\ & & \hspace{1 in } { }   + v_s - 2\sqrt{2}gx^2 \phi \re ( v_1)+2 g^2x^4\phi^2 v_2]e^k \ , , \
m^2_{i \bar j } & = & 2g^2x^4 e^k ( k_{i\bar j}+k_i k_{\bar j } ) + ( e^k v_s + v_\mathrm{lift})_{,i\bar j } \nonumber \\ & & \hspace{1 in } { } - 2\sqrt{2}gx^2 \phi ( e^k \re v_1)_{,i\bar j } + 2 g^2 x^4 \phi^2 ( e^k v_2)_{,i\bar j } \ , , \end{aligned}\ ] ] while fermion mass matrix has @xmath302 - w ) \ , , \nonumber \\
m_{\phi j } & = & -e^{k/2}\sqrt{2 } g x^2 k_j \ , , \qquad m_{\lambda\lambda } = e^{k/2 } \frac{\partial_i f k^{i\bar j}}{2f_r } ( d_{\bar j } \bar w-\sqrt{2}g x^2 \phi k_{\bar j } ) \ , , \nonumber \\
m_{ij } & = & m^{(0)}_{ij } -\sqrt{2 } g x^2\phi m^{(1)}_{ij } \ , , \qquad m_{3/2 } = e^{k/2 } ( w-\sqrt{2}g x^2 \phi ) \ , .\end{aligned}\ ] ] this implies @xmath303g^2x^4e^k \nonumber \\ & & { } - 4\sqrt{2}g x^2\phi \re \biggl[(w+3 v_1 ) e^k + k^{i\bar j } ( e^k v_1)_{,i\bar j } \nonumber \\ & & \hspace{0.3 in } { } - m^{(0)}_{ij}k^{i\bar a } \bar m^{(1)}_{\bar a\bar b}k^{j\bar   b } -\frac{e^k}{4 f_r^2}(k^{i\bar j } k_{\bar j } \partial_i f ) ( k^{b\bar a}d_b w
\partial_{\bar a } \bar f ) \biggr ] \nonumber \\ & &   { } + 4g^2x^4\phi^2 \biggl[(1+v_2 ) e^k + k^{i\bar j}(e^k v_2)_{,i\bar j } \nonumber \\ & & \hspace{1 in } { } - |m^{(1)}_{ik}k^{k\bar j}|^2 -\frac{e^k}{4f_r^2}|k^{i\bar j } k_{\bar j } \partial_i f|^2\biggr ] \ , .\end{aligned}\ ] ]
99 a.  d.  linde , particle physics and inflationary cosmology , contemp .
concepts phys .   * 5 * ( 2005 ) 1 [ hep - th/0503203 ] .


the gas and dust that is ejected by dying stars and from which new stars will form , is constantly being processed chemically as well as physically ( _ e.g. _ , * ? ? ? * ) .
a particularly rich and diverse route is followed by carbon . in the cool surroundings of carbon stars  evolved stars that have dredged up large amounts of carbon to the stellar surface
more than 60 individual molecular species have been identified , including benzene , polyynes and cyanopolyynes @xcite .
in addition , the spectra reveal diverse dusty materials such as hydrogenated amorphous carbon ( hac or a - c : h ) and silicon carbide ( sic ) . as the mass - loss process strips the outer layers of the star ,
a hot ( @xmath0 k ) central white dwarf becomes exposed , whose strong uv irradiation further processes the ejecta and makes them visible as a planetary nebula ( pn ) . a key component of the carbon dust inventory in pne are polycyclic aromatic hydrocarbons ( pahs ) and pah - like species , a class of large and hardy carbonaceous species whose formation mechanisms are unclear @xcite .
although not a single pah member has been identified in space , their characteristic spectral features at infrared wavelengths are observed throughout the universe , from which it is inferred that they reside ubiquitously in space and could lock up as much as 15% of the cosmic carbon ( _ e.g. _ , * ? ? ?
two other large aromatic species have been recently identified in space : c@xmath1 and c@xmath2 .
these are the best - known members of the family of fullerenes , a class of molecules made of hexagonal and pentagonal aromatic carbon rings , fused in the shape of a hollow sphere or ellipsoid .
the most abundant of these molecules , the buckminsterfullerene c@xmath3 , has the structure of an old - fashioned black and white soccer ball .
the near - spherical carbon configuration of these two molecules , which have a closed surface , a closed - shell electron distribution and an almost unstrained network , results in a very high stability against dissociation and prolonged exposure to high temperature .
fullerenes , in particular c@xmath3 , have peculiar and appealing photochemical , electrochemical and physical properties , which can be exploited in various fields , from nanotechnology to medicine , to space science .
c@xmath3 and c@xmath4 were discovered during laser ablation experiments on graphite targets , aiming to study long carbon chains in interstellar clouds @xcite .
because of their remarkable stability , fullerenes appeared particularly suited to survive the harsh conditions of the interstellar medium .
their unique properties indicate that they may play an important role in organic astrochemistry and astrobiology .
cosmic fullerenes remained elusive until the recent discovery of c@xmath3 and c@xmath4 in the planetary nebula tc 1 by @xcite .
fullerenes have since been confirmed in many more evolved star environments @xcite , as well as in young stellar objects @xcite , reflection nebulae ( rne ) @xcite and photodissociation regions @xcite .
fullerenes can be efficiently synthesized in the laboratory by the vaporization of carbon rods in an electric arc @xcite and from hydrocarbon combustion under optimised conditions @xcite .
however , the formation routes of fullerenes in space are still unknown .
we review the currently known formation mechanisms of fullerenes in the laboratory (   [ earth_sec ] ) , showing why these methods would not work in space (   [ not_working_sec ] ) .
we propose an alternative pathway , consistent with astrophysical conditions , based on the photo - processing of a family of carbonaceous species which we name `` arophatics '' (   [ our_mechanism_sec ] ) .
we discuss the observations of fullerenes in pne (   [ observations_sec ] ) and finally we summarize our conclusions (   [ conclusions_sec ] ) .
the original discovery of c@xmath3 and c@xmath4 dates back to 1985 with the experiment of @xcite . in this experiment ,
carbon species were vaporized from a graphite target into a he flow with tunable pressure / density .
the vaporization was done by a pulsed laser , the resulting carbon clusters were expanded in a supersonic molecular beam , photoionized and their masses measured by a time - of - flight mass spectrometer .
vaporization at low he pressure ( less than 10 torr ) leads to a broad distribution of clusters , with 38 to 120 carbon atoms ( always even numbers ) , with c@xmath3 and c@xmath4 present but not the dominant species . at a pressure of 760 torr
the c@xmath3 and c@xmath4 peaks clearly dominated .
when the expansion of the he+carbon cluster mixture was delayed , the resulting mass distribution was completely dominated by c@xmath3 - 50% of total large cluster abundance , and c@xmath5 , with 5 % of total large cluster abundance .
these distributions have been interpreted as due to the increasing number of he - cluster and cluster - cluster collisions , resulting in the `` survival of the fittest '' , i.e. the more stable clusters .
the famous icosahedral soccer ball structure was proposed for the first time for c@xmath6 , and later confirmed by spectroscopic studies @xcite .
however , the elementary reaction mechanisms occurring during the `` thermalization '' process remained unknown .
@xcite experimentally studied the gas - phase formation of carbonaceous compounds .
their first experiment was essentially a replication of the kroto experiment , i.e. laser vaporization of a graphite target in a quenching atmosphere of he or he / h@xmath7 with pressures between 3.3 and 26.7 mbar .
the vibrational temperature of the laser - induced plasma was estimated to be between 4000 and 6000 k for the laser power densities used . in a second set of experiments ,
the laser - induced pyrolysis of c@xmath7h@xmath7 and c@xmath8h@xmath8 using a pulsed laser with high power densities was studied .
the resulting condensation temperatures were above 3500 k , comparable with the vaporization experiment . in a third set of experiments , laser - induced pyrolysis of the same hydrocarbon precursors
was carried out , but this time at a much lower temperature ( max 1700 k ) .
the analysis of the condensation products showed a striking effect of the temperature . during the high temperature experiments , only fullerene - like soot and fullerenes were produced . during the low temperature experiments , only soot and pahs
were formed ( 100 % pahs at @xmath9 1000 k ) .
the results clearly tell that the temperature determines which kind of condensates will be formed .
moreover , the two pathways seem to be mutually exclusive , at least under these experimental conditions : fullerenes and pahs can not be formed together .
it is important to note that high pressure has been used during the experiments to concentrate the condensation within a small volume , i.e. to mantain a high density of the condensing species .
the laboratory experiments described in the previous sections do not tell us how fullerene is formed from graphite and hydrocarbons , but only under which conditions this happens .
the elementary reaction mechanisms involved in the formation of fullerenes can be investigated performing quantum chemical molecular dynamics simulations @xcite . in one set of simulations
a hot carbon vapor was reproduced by putting 40 randomly oriented c@xmath7 molecules in a tiny 20  cubic box @xcite .
the system was let to evolve , and extra c@xmath7 molecules were added at fixed times in order to simulate an open environment . during the first stage ,
the system was kept at a constant temperature of 2000 k and giant carbon cages self - assembled . during the second stage of the simulation
the temperature was raised to 3000 k , and the shrinking of these hot giant fullerenes down to c@xmath10 was observed @xcite .
the formation of the giant cages starts with nucleation of polycyclic structures ( hexagons and pentagons ) from entangled polyyne chains .
this is followed by growth of the structure and finally cage closure similar to the self - capping of open - ended single - walled nanotubes .
because of the rapid gain of energy due to cage closure , the giant fullerenic cages are produced in a vibrationally highly excited state .
the excess energy has to be dissipated , either by unimolecular dissociation or collision with other carbon clusters or carrier gas atoms .
the simulations show that the newly produced giant fullerenes inevitably shrink to smaller sizes .
all the road maps to fullerene formation proposed by previous models were associated with intermediate structures that are in thermodynamic equilibrium , while a hot carbon vapor is indeed a system far from thermodynamic equilibrium . in the model developed by irle , morokuma and collaborators , for the first time
the dynamic self - assembly of fullerene molecules occurs as an irreversible process emerging naturally under the nonequilibrium conditions typical of a hot carbon vapor .
the experimental results of @xcite have presented evidence of the spontaneous tendency of electron - irradiated graphite to include pentagons in its hexagonal network , hence form curved structures .
this has been further confirmed by @xcite , who have shown that fullerenes can form directly from graphene in a similar fashion . in the experiment of chuvilin and co - workers ,
a sheet of graphene was exposed to a beam of energetic 80 kev electrons .
the energy transferred to the graphene results in the loss of carbon atoms from the edges of the sheet .
if the size of the graphene flake is not too big ( less than a few hundreds of carbon atoms ) , the loss of carbon atoms at the edge will lead to the formation of pentagons , which triggers the curving of graphene into a bowl - shaped structure .
further carbon removal from the edges using the electron beam will reduce the size of the curved structure until it is sufficiently small to zip up its open edges and isomerize into a closed fullerene structure .
@xcite claim having identified the fundamental processes leading to the formation of fullerenes in a recent experimental work , based on laser vaporization of pure c@xmath3 in the presence of carbon vapor . according to this study
, fullerenes would self - assemble bottom - up through a closed network growth ( cng ) mechanism based on the ingestion of atomic carbon and c@xmath7 clusters .
it should be noted that the work from @xcite provides experimental evidence of the growth of larger fullerenes from pre - existing c@xmath3 only .
because the experiments show that the cng of larger fullerenes does not result in the production of c@xmath6 , it is deduced that c@xmath3 formation must be a result of cng from smaller fullerenes
. however , the initial formation mechanism of such smaller fullerenes is still under debate .
although all the chemical ingredients and the required temperatures for graphite vaporization and hydrocarbon combustion can be found in astrophysical environments , the densities are far too low to proceed to fullerene formation on reasonable timescales . in space , the most favourable conditions for fullerene formation via vaporization / combustion are found in a post - shock gas . in such environments the fullerene building blocks , i.e. c@xmath11 groups and polycyclic species coming from the fragmentation of pahs @xcite , ought to exist in the required vibrationally excited states .
however , the post - shock carbon densities are low and represent a serious obstacle to fullerene formation .
we have derived the scaling rule ( eq .  [ scaling_eq ] ) relating the time @xmath12 , required for fullerene condensation , with the initial density @xmath13 of the carbon gas .
the details of the calculation are reported in the appendix .
this scaling rule has then been used to estimate the time necessary for the coagulation of fullerenes from a carbon gas with interstellar densities ( @xmath14 10@xmath15  carbon  @xmath16 ) .
the scaling rule is : @xmath17 where @xmath18 and @xmath19 are the condensation timescales corresponding to densities @xmath20 and @xmath21 respectively . knowing the density and timescale in simulations , @xmath20 and @xmath18 , and the density of the interstellar gas , @xmath21
, we can then calculate the fullerene condensation time in space , @xmath19 .
we adopt @xmath20  =  10@xmath22  carbon  @xmath16 and @xmath18  =  50 ps @xcite .
the hydrogen density of the post - shock gas is @xmath9  1  @xmath16 , but carbon is 10@xmath15 less abundant , so @xmath21  =  10@xmath15  carbon  @xmath16 . from eq .
[ scaling_eq ] we obtain @xmath19  =  1.6  @xmath23  10@xmath24 yr , which is longer than the age of the universe .
we also estimate the carbon density @xmath21 required to form fullerenes in a shocked gas before the temperature drops below 1000 k. at this temperature the condensation of fullerenes can not occur . for a typical interstellar shock with a velocity of 125  km  s@xmath25 the cooling time ( after which the temperature of the post - shock gas drops below 1000 k ) is about 2@xmath2310@xmath26 yr . under these conditions ,
the required carbon atom density for the condensation of fullerenes is @xmath21  =  3@xmath2310@xmath27  carbon  @xmath16 , and in space such densities do not exist outside of stellar atmospheres .
@xcite propose a scheme for the direct formation of c@xmath6 from pahs in the reflection nebula ngc 7023 based on graphene dissociation experiments from @xcite . in ngc 7023 fullerenes and pahs appear to coexist in the same location close to the central star , while only pahs are detected further away from the star .
according to @xcite , a similar phenomenon to the one described by @xcite occurs in ngc 7023 , but induced by photo processing instead of electron bombardment .
the absorption of uv photons from the central star first completely dehydrogenates the pahs @xcite , producing graphene flakes .
further photon absorption induces the loss of carbons from the edges of the graphene sheet .
this , according to @xcite , gives rise to pentagonal defects which will induce the curvature of the sheet and the subsequent closure of the cage to form a c@xmath3 molecule .
there are some issues related to the applicability of the results from such experiments to interstellar conditions .
first of all , it is not possible in space to tune the number of carbons removed from the graphene flake in order to get the right size required for cage closure .
the number of ejected carbons depends on the rather slow photodissociation process , which is governed by the dissociation parameter @xmath28 .
this parameter describes the dissociation rate of a highly excited pah molecule using an arrhenius law @xcite , and to a first approximation can be seen as the energy that needs to be concentrated over a bond in order to break it and release a fragment .
the dissociation probability ( and hence the dissociation rate ) decreases for increasing @xmath28 because more energy is required in the bond that has to be broken .
the value of @xmath28 has been derived from experimental data only for very small pahs ( up to 24 c - atoms ) with a very open carbon skeleton .
the extrapolation of such results to larger ( @xmath9 50 c - atoms ) and compact astrophysically relevant pahs dissociating under interstellar conditions is very uncertain , leading to values of @xmath28 ranging from 3.65 to 5.6 ev @xcite .
this corresponds to values for the dissociation probability going from 0.5 to 3@xmath2310@xmath29 respectively . to form fullerenes during the lifetime of ngc 7023 ,
the mechanism proposed by @xcite needs a specific tuning of the parameters involved .
bacause the size of pahs can not be controlled as during experiments , the number of carbons in the precursor molecule is severely constrained to about 70 atoms , very close to the 60 carbons required to form c@xmath3 .
more importantly , the most favourable conditions for dissociation are chosen , i.e. , the lowest value for the dissociation parameter , @xmath28 = 3.65 ev . however
, this choice does not appear physically motivated .
indeed , if another legitimate value is adopted for the dissociation parameter , e.g. , @xmath28 = 4.6 ev ( the value falling in the middle of the range determined for @xmath28  see * ? ? ?
* ) , the decrease of the dissociation rate is such that the fragmentation process is far too slow and the conclusion suggested in @xcite no longer holds : pahs with 70 carbon atoms require longer than the lifetime of ngc 7023 to fragment , therefore they can not be the precursors of the c@xmath3 observed at the present epoch .
thus , none of these mechanisms are feasible in the astrophysical settings where fullerenes have been detected , and hence , alternative routes are needed .
trace amounts of fullerenes have been identified among the dissociation products of uv - irradiated hacs @xcite .
together with fullerene emission , the three pne tc 1 , smp lmc 56 and smp smc 16 show broad emission plateaus between 6  9 @xmath30 m and between 10  13 @xmath30 m , compatible with emission from h - rich amorphous carbon nanoparticles @xcite .
this suggests that a viable formation route in space could exist that starts from hac processing ( _ e.g. _ , * ? ? ?
* ; * ? ? ?
in this study we focus on the formation of fullerenes in the circumstellar environments of proto - planetary nebulae ( ppne ) and pne .
hacs are an abundant carbonaceous dust component in both these types of objects @xcite . in addition ,
one of the main characteristics of pne is their strong radiation field , which will process the hacs .      in their experiments ,
@xcite and @xcite showed that hac / a - c : h decomposition leads to a low density material ( @xmath31 g  @xmath16 ) .
in particular , @xcite and @xcite showed that the decomposition end - product is an aerogel - like material consisting of weakly - connected aromatic `` proto - graphitic '' or pah - like clusters in a `` friable network '' .
the dehydrogenation of small carbonaceous particles , containing less than a few hundred c atoms , will likely lead to structures that resemble the particles seen by @xcite in their experiments and/or the locally aromatic polycyclic hydrocarbons ( laphs ) proposed by @xcite .
the particles in the @xcite experiments are found to be small alkanes , unsaturated carbon - chain radicals and small dehydrogenated pah - like species .
after the release of these lighter molecules @xcite observed a strong mass peak near 500amu ( @xmath32 c atom , proto - graphitic , aromatic clusters ) having ir spectra similar to that of hac / a - c(:h ) with absorption and emission features at 3.3 , 3.4 and @xmath33 m .
the results of these experiments and recent modelling @xcite therefore imply an intimate mix of aromatic clusters and aliphatic groups in small carbonaceous species .
we note that the laph structures studied by @xcite have compositions in the range c@xmath34h@xmath35 to c@xmath36h@xmath37 and contain aromatic domains with only a few rings per domain , which are linked by bridging aliphatic ring systems , and should have ir spectral properties similar to hac / a - c(:h ) materials ( _ e.g. _ , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the conditions required for hac dehydrogenation are common in space , we can thus expect the formation of laph - like species where the aromatic components are bridged by olefinic carbon - containing rings and chains .
the key component in these structures is the presence of the aliphatic bridging groups linking the aromatic clusters , in what we here call `` arophatic '' structures .
such structures , analogous to an aliphatic periphery on pah molecules ( _ e.g. _ , * ? ? ?
* ) , could provide sites for molecular hydrogen formation in pdrs @xcite .
size - dependent @xmath38 band considerations indicate that the intrinsic aromatic clusters are indeed most likely to consist of ` isolated ' two- and three - ring systems in such small hydrocarbon particles @xcite .
the non - aromatic carbon atoms are in the form of `` more flexible '' short olefinic and/or aliphatic bridging structures such as ch = ch and ch@xmath7ch@xmath7. this generic type of structure was also proposed by @xcite , who suggested that aliphatic cyclo - hydrocarbons are the precursors to aromatic cluster formation in hac materials .
this kind of configuration , with aliphatic chains attached to aromatic sites , is also found in aromatic - aliphatic linked hydrocarbon ( aalh ) structures , which appear to have an important role in the nucleation of nanoparticles in soot formation in flames @xcite .
it should be noted that aalhs are small and relatively simple molecules , such as 1,2-di(naphthalen- 2-yl)ethane ( dne ) , while our proposed `` arophatic '' clusters are bigger and more complex structures .
the aliphatic bridging structures must be quasi - orthogonal to the aromatic species that they bridge and will be unable to transform into aromatic structures without compromising ( _ i.e. _ , destroying ) the entire structure . these active sites could , likely , interchange between @xmath39 olefinic and radical @xmath40 states with dangling bonds , in a kind of isomerisation `` flipping '' process , which could play a key role in cage closure and fullerene formation in the solid state .
the work of @xcite indicates that small clusters within a - c(:h ) materials are chain - like , while the larger clusters are three - dimensional , cage - like , @xmath39 structures and as @xcite show these cage - like structures have a natural tendency to close - up and form fullerene - like and fullerene structures .
the `` arophatic '' clusters are structures evolving from larger , amorphous hydrocarbon particles and , from a topological point of view , they will inevitably end up with a tube - like , cup - like or cage - like appearance ( _ e.g. _ , see the two right hand species in fig .  1 of * ? ? ?
such structures are ` curled - up ' or ` folded - over ' , defected graphite ( dg ) or graphene sheets @xcite , and we can therefore represent them by their ` flattened - out ' equivalents as per fig .
12 in @xcite but with aromatic domains reduced in size to only a few rings per aromatic cluster ( fig .
[ fig_arophatics_0 ] ) .
[ b ]        the aromatisation of the aliphatic bridges to form pentagons can provide a viable route to fullerene formation in the solid phase via these `` arophatic '' species .
the key process is the dehydrogenation of the aliphatic and olefinic `` bridges '' which naturally leads to aromatic pentagon formation , as shown in fig .
[ fig_schemes ] for small a - c : h particle fragments .
it appears that there is a link between small carbonaceous particle dehydrogenation in intense uv radiation fields and the formation of fullerenes , which could explain the observation of fullerenes in circumstellar regions associated with small hydrocarbon grains ( _ e.g. _ , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the mechanism that we propose for fullerene formation in space is the following .
the starting point is large , nm - sized , h - rich a - c : h particles containing more than a hundred carbon atoms .
uv photolysis by photons with energy above 10 ev dehydrogenates these grains @xcite and leads to further changes in the molecular structure  a process that we call * stir * ( * s*tructural * t*ransformations in * i*ntense * r*adiation fields ) . stir processing requires constant heating ( by photon absorption ) to keep the particle in a vibrationally excited state that promotes structural transformations leading to the closing of cage - like structures within the a - c : h nanoparticle .
amorphous carbon nanoparticles can be maintained at thermal - equilibrium temperatures of the order of 100 - 150  k by photon absorption ( for h - rich a - c : h particles with band gaps larger than 1 ev ) even at thousands of astronomical units from the central stars of pne , due to their low emissivity at the long wavelengths where they radiate thermally @xcite .
however , at these same distances , almost completely dehydrogenated nanoparticles have temperatures of only tens of kelvin , as expected .    in order to form fullerenes ,
the arophatic clusters resulting from stir processing must be large , containing more than a hundred carbon atoms , as exemplified in figs .
[ fig_arophatics_2 ] and [ fig_movie2_shot ] .
a key step in the route to fullerenes is the formation of 5-membered rings .
this can occur via dehydrogenation during the stir processing of the parent a - c : h or of the emerging arophatic cluster , in a similar fashion as shown in fig .
[ fig_schemes ] for small a - c : h particles .
dehydrogenation is a viable way to introduce 5-membered rings , as indicated by @xcite .
this work shows that structures with a high degree of curvature result from dehydrogenation reactions occurring in aromatic hydrocarbon compounds .
these reactions involve the loss of hydrogen with rearrangement of the structure and ring closure which leads to five - membered rings , and hence curvature .
uv irradiation of hac nanoparticles not only causes dehydrogenation ( * ? ? ? * and references therein ) , but also forces the structure to curl up because of the introduction of pentagonal rings .
in fact , the uv - induced dehydrogenation _ is _ the cause of the formation of 5-membered rings .
thus , in a single step , dehydrogenation fulfills two essential conditions for the formation of fullerene molecules : 1 ) removing hydrogen atoms and 2 ) triggering pentagonal ring formation in the evolving structure .
h@xmath41 ) , precursor of cosmic fullerenes .
carbon is represented in dark grey , hydrogen in light gray . the color version of this figure ( carbon in red , hydrogen in gray ) is available in the electronic edition of the _
astrophysical journal_. ]    h@xmath41 , precursor of cosmic fullerenes and schematically represented in fig .
[ fig_arophatics_2 ] .
carbon is in dark grey , hydrogen in light gray .
the movie , with carbon in red and hydrogen in gray , is available in the electronic edition of the _
astrophysical journal_. ]    the fullerene precursor is a big and complex arophatic cluster , with the aromatic domains being non - co - planar .
the 3-d structure favors the formation of 5-membered aromatic rings , and then curvature , even in such a big structure , via photo - induced dehydrogenation .
the particle emerging from the dehydrogenation process is a giant fullerenic cage similar to the ones formed in simulations @xcite .
the initial , stir - driven ( catastrophic ) photolysis and the following cage closure induced by the formation of pentagonal rings will result in a highly excited structure .
the cage will release the excess energy via ejection of c@xmath42 molecules , shrinking down to the islands of stability represented by c@xmath3 and c@xmath4 @xcite .
the down - sizing occurs through unimolecular decomposition of the vibrationally excited cage , resulting in the ejection of carbon molecules . to model the unimolecular dissociation we adopt an arrhenius law @xcite , which describes the evolution of the dissociation rate as a function of the temperature of the dissociating species : @xmath43 the quantity @xmath44 is the rate for carbon ejection measured in carbon atom lost per second ( c / s ) , @xmath45 the vibrational temperature of the fullerenic cage and @xmath46 the boltzmann constant .
the unimolecular decomposition is governed by two parameters : the pre - exponential factor @xmath47 and the dissociation energy @xmath28 .
to estimate the value of these two parameters , we use two sets of simulations which follow the shrinking of giant fullerenic cages at two different temperatures , 2000 and 3000 k , and provide the corresponding dissociation rates .
the average dissociation rate @xmath44 is 1.5@xmath2310@xmath27 c / s and 3@xmath2310@xmath27 c / s for @xmath45  =  2000 k and @xmath45
=  3000 k respectively @xcite .
we obtain the following values for @xmath47 and @xmath28 : @xmath48    the time @xmath49 required for a big fullerenic cage to shrink down to c@xmath1 is given by : @xmath50 where @xmath51 is the number of carbons that have to be ejected from the cage and @xmath44 is the ejection rate from eq .  [
eq : arrhenius_generic ] .
we consider a fullerenic cage containing 144 carbon atoms , i.e. with size comparable to the size of the particles considered in simulations @xcite , therefore , @xmath51 = 84 .
the shrinking timescales for various temperatures of the initial cage are reported in table  [ tab : shrink_time ] .
.[tab : shrink_time ] timescale , @xmath49 , for shrinking of the fullerenic cage c@xmath52 to c@xmath1 .
[ cols="^ , < , < " , ]     this process requires relatively short time scales : to shrink a large ( 144 carbon atoms ) fullerene cage to c@xmath1 requires only about 0.8 s at 200 k , and increases to 27 years at 100 k.    to the best of our knowledge , our extrapolation is the first attempt to characterize the shrinking of a giant fullerenic cage under interstellar conditions ( e.g. , very low temperatures with respect to the usual laboratory conditions ) .
the value of 0.36 ev that we have calculated for @xmath28 should be considered as a lower limit , having been determined under the most favourable conditions for fullerene formation on earth , as indicated by simulations ( nevertheless the optimized conditions for industrial production of fullerene are quite similar ) .
the issue that we are facing here is another example of the well known problem related to the extrapolation to extraterrestrial conditions of the results from statistical unimolecular dissociation theories @xcite .
the proper evaluation of @xmath28 and of the uncertainty associated with it requires the development of a specific model which will be the subject of a following paper .
our present scope is to provide a first indication of the timescales associated with the shrinking process in space . in the case of pah dissociation under extraterrestrial conditions ,
the detailed evaluation of the dissociation parameter @xmath28 resulted in a factor of 10 uncertainty on the calculated fragmentation timescales @xcite .
because of the similarities between large pahs and fullerenes @xcite , we can consider a @xmath53 uncertainty for the fullerene shrinking timescale as a very conservative assumption .
we obtain then a maximum shrinking time of 2670 years at @xmath45=100 k ( the least favourable condition ) .
this is fully consistent with the age of young pne where fullerenes have been observed .
c@xmath3 and c@xmath4 are the only , and therefore smallest , c@xmath54 fullerenes with @xmath55 70 satisfying the isolated pentagon rule ( ipr ) .
this implies ( almost ) spherical , thermodynamically favored , closed shell structures particularly suited to survive harsh conditions .
the shrinking stops at c@xmath1 because of the very high arrhenius dissociation energy @xmath28 for c@xmath11 elimination to c@xmath56 , @xmath57  10  ev @xcite .
the estimated temperatures of a few hundreds of kelvin for fullerenes detected in planetary nebulae @xcite are far from the approximately 1000 k required to overcome this decomposition barrier @xcite .
our model proposes the formation of a single c@xmath3 ( c@xmath5 ) molecule from a single large arophatic cluster .
this implies that the parent a - c : h particles must be nm - sized , and therefore contain a hundred or more carbon atoms , because of the dissociative loss of carbon ( and hydrogen ) atoms that will occur during structural transformation .
note that the large arophatic clusters can be further fragmented by the same radiation field that caused their formation from the parent a - c : h grain .
most of the formed arophatics will fragment into structures other than fullerenes , resulting in a fullerene formation process that is rather inefficient , consistent with the low fullerene abundance observed in space so far .
the peculiarities of the pne where fullerenes have been detected are their intense uv radiation fields , sufficient to sustain photo - chemical processing of hacs , and the ablation of carbonaceous material from the dense circumstellar matter formed during their agb phase .
thus , in these objects abundant , nm - sized , h - rich , a - c : h dust is being liberated and undergoing catastrophic uv - photolysis , which leads to stir processing . in tc 1 , c@xmath3 and c@xmath4
have been detected far for the central star ( 8000 au  * ? ? ?
* ) , but even at this distance the uv irradiation is sufficient to keep the particles in a vibrationally excited state ( @xmath58 100 - 130  k ) and to trigger and sustain the stir processing of the a - c : h dust and the shrinking of the resulting giant fullerenic cages down to c@xmath3 and c@xmath4 .
it appears that the young pne where fullerenes have been observed are the optimised factories for their formation because their stellar emission peaks in the required 10 ev photon energy regime . in their pre - pne stages
these objects do not have sufficient uv photons to dehydrogenate and heat the dust , and in their later stages the radiation field is so much harder that dissociative losses will destroy the h - rich carbonaceous particles before they can re - structure to form fullerenes .
we studied the formation of fullerenes c@xmath3 and c@xmath4 in planetary nebulae , for which a satisfactory formation path is still missing .    in space ,
the low density of the gas precludes the formation of fullerene materials following known vaporization or combustion synthesis routes even on astronomical timescales .
the direct formation of c@xmath3 through dissociation - induced curvature of dehydrogenated pahs requires a very specific tuning of the dissociation parameters , and thus appears unlikely .
the scheme for fullerene formation in space that we propose is based on a top - down approach .
the top of the tree is represented by large a - c : h / hac particles , i.e. materials whose presence in space has been firmly established .
uv photolysis promotes structural transformations in the parent a - c : h particle with the formation of what we call `` arophatic '' clusters : 3-d , hollow structures characterized by the presence of aromatic clusters linked by aliphatic bridging groups .
uv - induced dehydrogenation of the emerging arophatic cluster ( or the parent a - c : h ) introduces pentagonal rings , and hence curvature , in the evolving structure .    the result is a large , vibrationally - excited fullerenic cage that shrinks down to the stable molecules c@xmath3 and c@xmath4 through emission of c@xmath7 groups . the very high energy barrier for c@xmath7 ejection from c@xmath3 prevents further dissociation .
the young pne where fullerenes have been observed appear as the optimised factories for the formation of these molecules . in their pre - pne stages
the formation of fullerene - like structure is not favoured because these objects lack sufficient uv photons to dehydrogenate and heat the dust . in their later stages
the radiation field is so much harder that it will destroy h - rich carbonaceous particles before they can re - structure to form fullerenes .
we would like to thank g. lavaux for the realization of the movie and fig .
[ fig_arophatics_2 ] .
e.r.m , j.c . and e.p .
acknowledge the support from the national science and engineering research council of canada ( nserc ) .
wishes to acknowledge the support from a marie curie intra - european fellowship within the 7th european community framework program under project number 272820 .
this research has made use of nasa s astrophysics data system .
52 natexlab#1#1    , j. , cami , j. , peeters , e. , et  al .
2012 , , 757 , 41    , o. , & tielens , a.  g.  g.  m. 2012 , proceedings of the national academy of science , 109 , 401    , j. , bernard - salas , j. , peeters , e. , & malek , s.  e. 2010 , science , 329 , 1180    , j. , bernard - salas , j. , peeters , e. , & malek , s.  e. 2011 , in iau symposium , vol .
280 , iau symposium , 216227    , j. , heras , a.  m. , tielens , a.  g.  g.  m. , et  al .
2001 , , 546 , l123    , i. , le teuff , y.  h. , williams , p.  m. , & tielens , a.  g.  g.  m. 2000 , , 357 , 572    , s .- h .
, & violi , a. 2011 , proc .
inst . , 33 , 693    chuvilin , a. , kaiser , u. , bichoutskaia , e. , besley , n.  a. , & khlobystov , a.  n. 2010 , nature chemistry , 2 , 450    , g.  c. , de marco , o. , whitney , b.  a. , et  al .
2011 , , 142 , 54    dunk , p.  w. , kaiser , n.  k. , hendrickson , c.  l. , et  al .
2012 , nat .
commun . , 3 , 855    , a. , van loon , j.  t. , woodward , c.  e. , et  al .
2012 , , 421 , l92    , a.  c. , & robertson , j. 2004 , phil .
lond . a , 362 , 2477    , k.  a.  k. , mutschke , h. , & jger , c. 2011 , , 528 , a56    , d.  a. , kameswara rao , n. , & lambert , d.  l. 2011 , , 729 , 126    , d.  a. , manchado , a. , garca - lario , p. , et  al . 2010 , , 724 , l39    , d.  a. , iglesias - groth , s. , acosta - pulido , j.  a. , et  al .
2011 , , 737 , l30    , c. , cami , j. , bouwman , j. , peeters , e. , & min , m. 2011 , , 536 , a54    , m. , gaessler , w. , hayano , y. , et  al .
2003 , , 589 , 419    , v.  i. , & duley , w.  w. 2000 , , 543 , l85    , v.  i. , tereszchuk , k. , duley , w.  w. , & bernath , p. 2001
, , 558 , l129    , j.  b. , mckinnon , j.  t. , makarovsky , y. , lafleur , a.  l. , & johnson , m.  e. 1991 , , 352 , 139    irle , s. , zheng , g. , wang , z. , & morokuma , k. 2006 , the journal of physical chemistry b , 110 , 14531 , pmid : 16869552    , c. , huisken , f. , mutschke , h. , jansa , i.  l. , & henning , t. 2009 , , 696 , 706    , a.  p. 2012
, , 540 , a1 +    . 2012 , , 540 , a2 +    . 2012 , , 542 , a98 +    , w. , fostiropoulos , k. , & huffman , d.  r. 1990 , chem .
lett . , 170 , 167    , w. , lamb , l.  d. , fostiropoulos , k. , & huffman , d.  r. 1990 , , 347 , 354    , h.  w. , heath , j.  r. , obrien , s.  c. , curl , r.  f. , & smalley , r.  e. 1985 , , 318 , 162    , s. 2009 , international journal of astrobiology , 8 , 161    , s.  d. , goodwin , d.  g. , anderson , m.  s. , & anderson , j.  r. 1995 , , 51 , 9973    , e.  r. , jones , a.  p. , & tielens , a.  g.  g.  m. 2010 , , 510 , a36
+    . 2010 , , 510 , a37 +    .
2004 , masters of the flame : industrial production of fullerenes becomes a reality ( http://www.nano-c.com/technologies.as )    , j.  r. , cernicharo , j. , goicoechea , j.  r. , gulin , m. , & asensio ramos , a. 2007 , , 661 , 250    , e. , tielens , a.  g.  g.  m. , allamandola , l.  j. , & wolfire , m.  g. 2012 , , 747 , 44    , s. , stranger , r. , & duley , w.  w. 2003 , , 594 , 869    , e. , & hornekr , l. 2008 , , 679 , 531    , k.  r.  g. , smith , k.  t. , & sarre , p.  j. 2012 , , 421 , 3277    , r.  h. , simpson , j.  p. , odell , c.  r. , et  al .
2011 , , 410 , 1320    , a. , & duley , w.  w. 1996 , , 472 , l123 +    , a. , duley , w.  w. , & pinho , g.  p. 1997
, , 489 , l193 +    , k. , werner , m.  w. , ingalls , j.  g. , et  al .
2010 , , 722 , l54    , f.  w. 1984 , journal of applied physics , 55 , 764    taylor , r. , hare , j.  p. , abdul - sada , a.  k. , & kroto , h.  w. 1990 , j. chem .
, 1423    , a.  g.  g.  m. 2005 , the physics and chemistry of the interstellar medium ( cambridge univ .
press , cambridge )    , s. , andersen , j.  u. , gottrup , c. , hvelplund , p. , & pedersen , u.  v. 2001 , phys .
, 87 , 073401    , d. 1992 , , 359 , 707    , a. 2004 , combustion and flame , 139 , 279    , y. , & kwok , s. 2011 , , 730 , 126    zheng , g. , irle , s. , & morokuma , k. 2005 , the journal of chemical physics , 122 , 014708    , g. , wang , z. , irle , s. , & morokuma , k. 2007 , j. nanosci .
nanotechnol . , 7 , 1662
we report here the details of the calculation of the scaling rule ( eq .  [ scaling_eq ] ) , which relates the time @xmath12 required for fullerene condensation with the initial density @xmath13 of the carbon gas .
consider an ensemble of @xmath59 particles that can occupy @xmath60 cells .
the probability of having one particle in one cell is @xmath61 .
the probability @xmath62 of having two particles in the same cell at the same time is described by the poisson statistics : @xmath63 we have that @xmath64 , where @xmath65 is the total volume and @xmath66 is the volume of a single cell .
we can then rewrite the single occupancy probability as @xmath67 .
the probability @xmath68 of having interaction between two particles at time @xmath69 can be written as : @xmath70 where @xmath71 , @xmath72 being the time and @xmath73 the characteristic time of a single interaction ( depending on the nature of the interaction itself ) .
the term @xmath74 represents the probability of interaction between two particles at a generic time @xmath72 .
the probability of having interaction between two particles up to the time @xmath69 is then given by the following equation : @xmath75 if @xmath74 is small , we can rewrite eq .
[ tilde2_eq ] as : @xmath76 the probability @xmath74 can be written in the following way :
@xmath77 the term @xmath78 includes the characteristics intrinsic to the interaction ( which we do not need to specify ) while @xmath79 is the probability from eq .
[ poisson_eq ] .
the approximation resulting in the right - hand term of the equation is valid for @xmath80 . combining eqs .
[ tilde3_eq ] and [ lambda_eq ] , with @xmath81 , we obtain @xmath82    from eq .
[ tilde4_eq ] we can derive the timescale for condensation : @xmath83 because both @xmath73 and @xmath84 are independent from the initial particle density @xmath13 , we can write the following scaling rule @xmath85 where @xmath18 and @xmath19 are the condensation timescales corresponding to densities @xmath20 and @xmath21 respectively .