universality is the amazing phenomenon when the critical behavior arising in different physical systems is characterized by a common set of universal quantities .
all systems belonging to a certain universality class share the same critical exponents , universal amplitude ratios and universal scaling functions .
typically , universality classes are determined by spatial dimension , symmetry and the range of interaction potentials . at times ,
two systems belonging to one universality class appear to be entirely disparate at first sight . in this paper
we discuss such a remarkable case , _ viz_. that of driven surface diffusion under detailed balance conditions and crystalline solids near their shear instability .
we show that these respectively two and three dimensional systems are described by a common field - theoretic hamiltonian .
we calculate the scaling behavior of this universality class using field - theoretic renormalization group ( rg ) methods .
surface growth problems constitute an important class of generic nonequilibrium phenomena @xcite .
particles are deposited on a surface and can diffuse around on it . in ballistic deposition processes such as molecular beam epitaxy ( mbe ) , desorption and bulk defect formation can often be neglected . after an initial transient , a steady state
is established which is characterized by time - independent macroscopic properties , provided that a suitable reference frame is chosen .
if one is primarily interested in universal , large - scale , long - time characteristics , the dynamical evolution can often be formulated in terms of a langevin equation which then can be analyzed using the methods of renormalized field theory .
several years ago , schmittmann , pruessner and one of us ( h.k.j . )
@xcite , denoted in the following as spj , extended a model by marsili , maritan , toigo , and banavar ( mmtb ) @xcite , to describe ballistic deposition processes with oblique particle incidence . in broad terms , their approach can be described as follows : it is assumed that there is no shadowing , i.e. , that there are no overhangs so that the growing surface can be described in terms of a height variable @xmath0 ( so - called monge representation ) . furthermore , it is assumed that particle desorption and bulk defect formation can be neglected , so that the ( deterministic ) surface relaxation processes are particle - number conserving .
then , the time rate of change of the surface - hight is formulated in an idealized continuum description  @xcite in terms of the divergence of a current .
since the absolute hight of the surface is irrelevant , these currents are constructed from the slope of the height @xmath1 and its derivatives .
next , shot noise modeling the random deposition of the particles is incorporated into the growth equation by adding a non - conserving stochastic term which dominates the conserving diffusional noise . since the oblique particle beam selects a preferred ( longitudinal ) direction in the substrate plane , the resulting langevin equation is necessarily anisotropic , and the primary relaxation mechanism is driven surface diffusion .
the interplay of interatomic interactions and kinetic effects , such as ehrlich - schwoebel barriers , generates an anisotropic effective surface tension which can become very small or even negative . due to the anisotropy , these mechanisms lead to a phase space with four different regimes , and the potentially scale - invariant behaviors in these regimes were studied by spj .
the most interesting case turned out to have a transversal instability which leads to the generation of longitudinal ripples via a continuous phase transition in contrast to a longitudinal instability which leads to the generation of transversal ripples by a ( fluctuation induced ) first order transition .
such ripple structures are commonly found in sandy deserts where they are generated under the influence of a steady wind as the driving force , see fig .
[ fig : ripples ] .
it is worth noting that mmtb and spj formulated their model under the assumption of an invariance with respect to tilts of the surface .
however , in a real setup , the orientation of the basic surface relative to the incident particle beam is generally fixed , and hence it is useful to generalize their model to go beyond the limiting tilt - invariant case . here
, in the present work , we do _ not _ assume tilt invariance
. we will see as we move along that the resulting growth equation is in detailed balance when an appropriate co - moving coordinate system is used and a certain model parameter vanish . note that we use the term detailed balance  not in a microscopic sense @xcite .
even if the microscopic laws lack detailed balance , coarse graining and the following neglect of irrelevant terms can lead to detailed balance in the asymptotic semi - macroscopic equation of motion for a driven system  @xcite . under the assumption that this is the case , we derive a quasi - static hamiltonian @xmath2 which describes the stationary distribution , and whose basic ingredients are derivatives of the surface height field @xmath3 .    the elastic free energy of crystalline solids with hexagonal or tetragonal ( including cubic ) symmetry shows an instability against shear when the elastic constant @xmath4 goes to zero .
the crystallographic unit cell undergoes an elastic deformation with a certain linear combination of components of the strain tensor as an order parameter . in the case
that no third order invariants of the order parameter are possible ( or vanish by chance in the cubic case ) , the transition to the structurally distorted crystal is continuous , and the order parameter can be reduced to the spatial derivatives of the displacement field @xmath5 along the crystallographic symmetry axis @xcite . augmenting the stretching energy with stabilizing bending terms and then reducing it to its parts that are relevant in the rg sense , we set up a generalized elastic energy that is of the same form as the aforementioned quasi - static hamiltonian @xmath2 with the displacement field @xmath6 taking on the role of the height field @xmath7 .
note , however , there is one important distinction between the two systems .
whereas in the driven diffusion problem , both the longitudinal component @xmath8 and the transversal components @xmath9 lie in the plane of the substrate , i.e. , in the plane orthogonal to the height @xmath1 , the coordinate @xmath8 in the elastic problem lies in the direction of the crystal axis and is therefore parallel to the displacement @xmath5 .
hence , the physical dimension of @xmath9 is @xmath10 in driven surface diffusion and @xmath11 in the elastic transition .
the outline of the present paper is as follows : in sec .
[ sec : model ] we set up field theoretic models for driven surface diffusion with detailed balance and elastic phase transitions , respectively , through augmenting and modifying the existing models for these systems .
we unify the two field theoretic models by formulating a hamiltonian that commonly describes both systems . in sec .
[ sec : rpt ] we present our renormalized perturbation theory .
we sketch our diagrammatic calculation , discuss the resulting renormalization group and its flow . in sec .
[ sec : scalingandexponents ] , we derive the scaling forms for several experimentally measurable quantities . for driven surface diffusion ,
we calculate critical exponents , and for elastic transitions we calculate logarithmic corrections to the leading mean - field behavior . in sec  [ sec : concludingremarks ] we give concluding remarks .
appendix  [ app : calcdiagrams1loop ] contains some details on our calculation of feynman diagrams .
appendix  [ app : addren ] discusses some intricacies related to the fact that the calculation of the specific heat for the elastic transitions requires additive renormalizations .
here , we set up our model for driven surface diffusion following the work of mmtb  @xcite and spj .
as a starting point , let us review the langevin equation presented by spj . in the to - interpretation , this langevin equation describes the hight @xmath0 of a @xmath12-dimensional growing hyper - surface as @xmath13 the incident particle beam has components @xmath14 and @xmath15 normal and parallel to the substrate plane , respectively .
as discussed in detail by mmtb and spj , the oblique particle incidence induces a fundamental anisotropy into the system which breaks the full rotational symmetry within the @xmath12-dimensional space of the substrate . to account for this anisotropy , it is useful to choose a coordinate system in which a particular axis , say the @xmath8-axis , is aligned with @xmath15 and then to write the spatial coordinates as @xmath16 .
the two uniform growth terms , the first two terms on the right hand side of eq .
( [ genlang ] ) , can be removed by a galilei transformation @xmath17 . in the following
, we will work in this co - moving frame .
@xmath18 is the surface current , and we model its longitudinal and transversal components following spj by @xmath19 note that under the assumptions made by mmtb and spj , @xmath20 . as we will discuss in more detail below
, we will work under different assumptions here .
the langevin noise @xmath21 accounts for the randomness in the particle deposition .
because this shot noise adds particles to the surface , it is commonly modeled as non - conserving noise with correlations of the form @xmath22 in principle , conserved diffusional noise could also be incorporated into the model . however , it turns out to be irrelevant compared to the non - conserved noise .
the character of the surface produced by the process depends on the values of the critical parameters @xmath23 and @xmath24 .
our model leads to qualitatively the same phase diagram in the space spanned by @xmath23 and @xmath24 as the spj model , see fig .  1 of ref
.  @xcite .
this phase diagram contains , in particular , a continuous transition from isotropic , edwards - wilkinson  @xcite type behavior to a state characterized by longitudinal ripples for @xmath25 and @xmath26 . in the following ,
we will focus on this transition .
the other transitions in the phase diagram which includes a fluctuation induced discontinuous transition to a surface with transversal ripples for @xmath27 and @xmath28 and a multi critical point for @xmath27 and @xmath29 are set aside for future work .
now , we review the symmetries and the physical contents of the langevin equation  ( [ genlang ] ) with its ingredients ( [ parcurr ] ) to ( [ shot - noise ] ) in some more depths . the work by mmtp and spj assumed that the physics of the system is invariant under an infinitesimal tilt of the surface by an angle @xmath30 , @xmath31 .
the assumption of this invariance enforces the relation @xmath20 mentioned above  @xcite .
however , this tilt - invariance is non - physical since the orientation of the substrate plane is fixed by the normal component @xmath14 of the incident current . in the following , we will refrain from imposing this condition . stability demands positive @xmath32 and @xmath33 as well as @xmath34 . at the transition for @xmath25 and @xmath26 that we are interested in , @xmath32 and @xmath35 are irrelevant .
the positive parameters @xmath33 and @xmath23 can be set equal to @xmath10 and a dimensionless positive parameter @xmath36 , respectively , via a simple rescaling of the coordinates @xmath16 and @xmath37 .
note that after this transformation , the coordinates , the time and the `` temperature '' scale as @xmath38 , @xmath39 , @xmath40 and @xmath41 , respectively , where @xmath42 is an inverse external length scale appropriate to discuss the infrared limit that will be used in the renormalization program to follow .
it thus follows from the langevin equations , eqs .
( [ genlang ] ) and ( [ shot - noise ] ) , that the shot noise @xmath21 and the height @xmath1 scale as @xmath43 and @xmath44 .
furthermore , the coupling constants scale as @xmath45 and @xmath46 which shows that @xmath47 is irrelevant and that the upper critical dimension is @xmath48 .
returning to the tilt - invariance and the fact that we do not assume it in the present work , we note that there is a possible contribution to @xmath49 of the form @xmath50 that was omitted by mmtb and spj on grounds of tilt - invariance , and we have to consider it here .
power counting reveals that @xmath51 , i.e. , this term is relevant .
even if we omitted it at the onset , it would be generated by the rg .
here , we include it in our model at this stage .
we are interested here in driven surface diffusion with detailed balance so that the statistical properties of the steady state can be derived from a free energyor hamiltonian @xmath52 .
for this to hold true , we have to be able to express the components of the surface current through functional derivatives of @xmath52 . because the absolute height of the surface should not matter , i.e.
, @xmath53 with a constant @xmath54 should be a symmetry , the functional variation of @xmath52 with respect to @xmath1 should be of the form @xmath55 which defines the surface current @xmath18 up to divergence - free additive part .
it is easy to see that this is possible if we equalize the coupling constants @xmath56 , and in the following we will do so .
then the hamiltonian is given by@xmath57 for further details and background on detailed balance in the context of driven surface diffusion , we refer to ref .
@xcite .
there is another symmetry present that we have not mentioned yet , viz .  the so - called up - down symmetry .
this symmetry , which physically means that mounds and valleys have the same form , translates for mbe with oblique incidence to a form - invariance under the simultaneous inversion @xmath58 and @xmath59 .
this type of up - down symmetry was assumed by mmtb and spj , and we also assume it here .
however , even though it is frequently assumed in mbe theory , up - down symmetry has to be taken with a grain of salt because the particle beam sets a preferred direction in space .
therefore , the up - down symmetry is strictly speaking more appropriate for interface than for surface models .
we keep this symmetry here firstly because we are primarily interested in the connection between and the mutual properties of driven surface diffusion and elastic phase transitions and secondly because the up - down symmetry simplifies the model . in a sense , this simplified model can be viewed as a leading order approximation for mbe with oblique incidence that can be improved later by including terms that break the up - down symmetry .
in fact , we plan to study the effects of breaking the up - down symmetry in the future .
note that this symmetry has been discussed in the literature  @xcite as a prerequisite for having a free energyin this kind of system .
however , we can add a further coupling proportional to @xmath60 to the hamiltonian ( [ hamsurf ] ) which explicitly destroys up - down symmetry . only in the one - dimensional case
this coupling is a total gradient and therefore vanishes after integration over the full one - dimensional space .
the langevin equation ( [ genlang ] ) now takes the form @xmath61 in the comoving - frame introduced above . to generate correlation and response functions of this stochastic process
, it lends itself to reformulate the langevin equation as a dynamical response functional @xmath62 so that averages can be derived from path integrals with the statistical weight @xmath63 . via introducing a response field @xmath64 , the functional @xmath62 ( in to - discretization ) that we extract with standard methods  @xcite from eqs .
( [ speclang ] ) and ( [ shot - noise ] ) reads @xmath65=\int d^{d}x\int dt\ \lambda \big\{\tilde{h}\big(\lambda^{-1}\partial_{t}h+\frac{\delta\mathcal{h}_{\text{surf}}[h]}{\delta h}\big)-\tilde{h}^{2}\big\}\ , , \label{def_j}\ ] ] and shows that our driven surface diffusion model obeys detailed balance @xcite .
the fluctuation - dissipation theorem ( fdt ) holds , and all stationary expectation values can be calculated from path - integrals with the statistical weight @xmath66 .    it is a well established fact that in all dynamical field theories that are invariant under a field shift of the type @xmath53 with a constant @xmath54 , the dynamical exponent is given by @xmath67 , where @xmath68 is the anomalous scaling exponent of the order parameter field @xmath1 .
as it should , the same applies here .
this can be verified in detail by noting that in the full dynamical field theory based on the response functional  ( [ def_j ] ) one does not need to introduce renormalization factors for the @xmath69 and @xmath70 terms when using minimal renormalization , i.e. , dimensional regularization in conjunction with minimal subtraction .
hence , we can calculate all of the critical exponents of driven surface diffusion with detailed balance from @xmath52 including the dynamical exponent .
an elastic phase transitions is a structural phase transition in which the crystallographic unit cell undergoes an elastic deformation when a normal mode becomes soft in response to the vanishing of an elastic constant or combination of elastic constants . in this section ,
we essentially follow the early work of cowley @xcite , and folk , iro , and schwabl ( fis )  @xcite . for background information on elastic phase transitions ,
we refer to these original papers as well as reviews by bruce and cowley  @xcite , rao and rao  @xcite , and schwabl and tuber @xcite .    in general ,
the order parameter for such a transition is a certain component of the lagrange strain tensor or a linear combination of its components .
depending on the symmetry of a crystal , its lagrange elastic energy density can be unpleasantly complicated even at harmonic order due to the presence of a host of different elastic constants . here , we restrict ourselves to crystals with tetragonal or higher symmetry , i.e. , crystals whose harmonic lagrange energy density can be obtained from the tetragonal case by imposing relations between certain elastic constants . as we will see below ,
crystals with lower symmetry such as orthorhombic or monoclinic symmetry have simple mean field behavior at their respective elastic transition and are hence less interesting from the standpoint of critical phenomena than the higher - symmetry crystals we focus on here .
when generalized to @xmath12 spatial dimensions , the lagrange energy density of tetragonal solids to second order in the strains reads @xmath71 where @xmath72 , and @xmath73 for a physical solid .
the lagrange strain tensor ( see , e.g. , refs .
@xcite ) is defined via the derivatives of the displacement field @xmath74 with respect to the reference space coordinates of the undeformed solid , @xmath75 with @xmath76 .
if we specialize the elastic constants @xmath77 to @xmath78 , @xmath79 , @xmath80 or @xmath81 , eq .  ( [ fren ] ) reduces to the harmonic energy density of a cubic or hexagonal crystal , respectively .
if all these conditions hold , eq .  ( [ fren ] ) describes an isotropic solid .
moreover , the hexagonal case with the additional relation @xmath82 reduces to the isotropic case after a suitable rescaling of the coordinates .
because of crystal anisotropy , sound velocities depend on the direction of propagation , and modes become soft only in certain selected directions , so - called soft sectors , which may be one- or two - dimensional in 3@xmath12 crystals .
the type of the soft modes depends on which elastic constant or combination of elastic constants vanishes . in the following , we focus on the shear instability that occurs when @xmath4 goes to zero .
the accompanying elastic phase transition is continuous if there are no third order couplings of the corresponding order parameter .
we assume that this is the case .
the upper critical dimension of the phase transition is @xmath83 , where @xmath84 is the dimension of the soft mode sector in fourier space , as was shown by fis .
this implies that , from the standpoint of critical phenomena , @xmath85 is the most interesting case since otherwise the upper critical dimension is lower then the physical dimension @xmath86 , and the transition is simply be described by mean - field - theory . in the remainder , we focus on @xmath85 . in this case , the @xmath87-dimensional soft sector is orthogonal to the polarization of the corresponding critical sound waves , i.e. , to the direction of the symmetry - axis in fourier space .
following fis , we radically reduce the displacement field to the critical displacement alone@xmath88 this step is justified if there are no third order couplings in the critical strain .
such couplings can be ruled out on symmetry grounds for the crystals we focus on except for the cubic case . for cubic crystals ,
we assume that third order couplings are absent for some given reason , e.g. , chance .
otherwise , the elastic transition becomes discontinuous .
we will comment on this issue and possible couplings to the non - critical displacements @xmath89 further below .    to figure out which parts of the elastic energy are relevant in the sense of the rg , it is useful to rank the different terms according to the order at which they contribute in the long - wavelength limit . to this end , let us introduce the smallness parameter @xmath90 , and set the derivatives to@xmath91 with the reduction ( [ reduz ] ) , the strain tensor becomes@xmath92 and we obtain the free energy density ( [ fren])@xmath93 where we have set @xmath94 and @xmath95 .
now we see that it is appropriate to rescale the vanishing elastic constant @xmath4 by@xmath96 and reduce @xmath97 to the leading terms proportional to @xmath98:@xmath99 of course , the anisotropic last term in @xmath100 with hypercubic symmetry vanishes in the hexagonal case .    at this stage ,
a comment regarding higher than harmonic terms in the initial lagrange energy density  ( [ fren ] ) is warranted .
as long as we ignore the uncritical displacement @xmath101 , the addition to the right hand side of eq .
( [ fren ] ) of third and fourth order terms in the strains leads to exactly the same form of @xmath102 as given in eq .
( [ redelaen ] ) .
the only effect of the higher order strain terms is to modify the constants decorating the individual terms in eq .
( [ redelaen ] ) , and hence their omission is justified .
however , if we included the noncritical displacements which scale as @xmath103 , then already the harmonic strain terms in the lagrange elastic energy density would lead to couplings of the form @xmath104 , @xmath105 and @xmath106 which contribute for @xmath107 at the same order in @xmath90 as the terms retained in eq .
( [ redelaen ] ) .
thus far , we only included stretching in our discussion of elastic energy densities . however , since we are interested in the limit @xmath108 where there is an elastic instability , bending terms will be important as they stabilize the crystal .
the bending terms that are relevant here are @xmath109 with @xmath110 in the hexagonal case . adding stretching and bending contributions , we obtain the model elastic energy or hamiltonian @xmath111 which is , up to the anisotropic cubic terms which are absent in the hexagonal case , of the same form as the hamiltonian ( [ hamsurf ] ) for driven surface diffusion . this is the hamiltonian that we will study , after suitable rescalings and renamings , in our renormalized perturbation theory .
form this point on , we will disregard the hyper - cubic terms , i.e. , the last term in the stretching energy density  ( [ redelaen ] ) and the term proportional to @xmath112 in the bending energy density  ( [ redbenden ] ) . as mentioned above , crystals with cubic symmetry are expected to have anomalous scaling behavior and are hence certainly interesting , however , the presence of the hyper - cubic terms makes the diagrammatic perturbation calculation significantly more difficult . though we have done this calculation , we feel that a useful discussion of the cubic case goes beyond the scope of the present paper
. we will cover elastic transitions in cubic crystals in a future paper  @xcite .
before we move on , let us make some comments on the relation of our hamiltonian @xmath113 to other , prominent model elastic energies that have been discussed in the literature .
first , for @xmath114 ( and @xmath115 ) , @xmath113 reduces to one of the two models studied by fis ( their model ii )
. we will refer to this model as the fis model .
secondly , in the limiting case of a fully isotropic crystal with @xmath116 , @xmath117 , @xmath118 , we have@xmath119^{2}\ , , \label{isotr}\ ] ] with@xmath120=\big\{\frac{c_{44}}{c_{11}}+\bigl(\partial_{\vert}u\bigr)+\frac{1}{2}\bigl(\nabla_{\bot}u\bigr)^{2}\big\}\ , , \label{invar}\ ] ] up to a total derivative .
@xmath121 $ ] is invariant under a type of a galilean transformation , @xmath122 with infinitesimal transversal @xmath123 .
when augmented with bending the bending term proportional to @xmath124 , this isotropic elastic energy density leads to a hamiltonian @xmath125 that is equivalent to that considered by grinstein und pelcovits ( gp ) in their seminal work on the anomalous elasticity of smectic liquid crystals  @xcite .
we will refer to this model as the gp model .
jointly , we will refer to the fis and gp models as the special models .      to further emphasize that both driven surface diffusion with detailed balance and elastic transition are described by the same field theoretic hamiltonian , we rescale the coordinates and @xmath5 in eq .
( [ hel ] ) . after renaming the order parameter field and the coupling constants in the rescaled version of eq .
( [ hel ] ) as well as renaming the order parameter field @xmath1 in eq .
( [ hamsurf ] ) , both our model systems are described by the hamiltonian @xmath126 where @xmath127 is the common order parameter field .
next , let us discuss the symmetry transformations of this hamiltonian and the combinations of its parameters that are invariant under the respective transformation .
knowing these combination will come in handy later on when we analyze the rg flow .
this hamiltonian is invariant under the longitudinal scale transformation @xmath128 the following combinations of coupling constants are invariant under this transformation : @xmath129 a further symmetry of the hamiltonian @xmath2 is the @xmath8-dependent field shift @xmath130 the galilean transformation with transversal @xmath123 @xmath131 is a symmetry if only and only if the conditions @xmath132 hold . of course , if @xmath133 , the hamiltonian is invariant under the inversion @xmath134 .
in this section , we outline our renormalized perturbation theory . for background on the methods that we use ,
we refer to the textbooks  @xcite .
the feyman diagrams of our perturbation theory consist of the gaussian propagator @xmath135 and two vertices , see fig .
[ fig : statelem ] .
the @xmath136-term in the denominator of the propagator will be omitted in the following because an expansion of any diagram in powers of this term will inevitably lead to unnecessary contributions that are convergent in the ultraviolet ( uv ) .
there is a 3-leg vertex @xmath137 \label{3-vert}\ ] ] and a 4-leg vertex @xmath138\ , .
\label{4-vertex}\ ] ]    we use dimensional regularization about @xmath139 dimensions in conjunction with minimal subtraction , i.e. , we use minimal renormalization .
our renormalization scheme reads @xmath140 where @xmath141 we choose the additive renormalization @xmath142 so that @xmath143 .
the longitudinal tilt - transformation ( [ longtilt ] ) leads to the exact relation@xmath144 note that formally @xmath145 in minimal renormalization . in general ,
i.e. , beyond minimal renormalization , @xmath146 is a positive or negative constant . in a momentum shell rg , for example , it depends on the momentum cut - off and on the coupling constants .
we calculate the 1-loop contributions to the superficially uv - divergent vertex functions @xmath147 , @xmath148 and @xmath149 .
some details of this calculation can be found in appendix  [ app : calcdiagrams1loop ] .
for the counter - terms that render their respective vertex function finite at 1-loop order , we obtain @xmath150 form these counter - terms , we extract the 1-loop renomalization factors @xmath151 note that if the system is invariant under the galilean transformation  ( [ galtraf ] ) which implies @xmath152 , the conditions ( [ galsym ] ) lead to the following two _ exact _ relations between the renormalization factors : @xmath153      now , we set up our rg equations and discuss the rg flow . to simplify the involved equations , we use in the following the shorthands @xmath154 for the renormalized versions of the invariant combination of coupling constants given in eq .
( [ invcoup ] ) . in terms of @xmath155 and @xmath156 ,
the rg @xmath157-functions are defined as @xmath158 where @xmath159 and @xmath160 denote the @xmath161-residues of the respective renormalization factors . to 1-loop order
, we obtain @xmath162 inspecting the relations ( [ galbez ] ) , in particular for @xmath163 , we figure that it is useful to define @xmath164 this quantity has a perturbation expansion in @xmath155 and @xmath156 which is to lowest order given by @xmath165 from the rg @xmath157-functions , we get the wilson and gell - mann  low functions @xmath166 note that for @xmath163 , which corresponds to the gp model where the ward identities ( [ galbez ] ) do hold , we have @xmath167 to arbitrary order in perturbation theory .
our 1-loop results are consistent with these relations and hence satisfy an important check . having these rg functions , we can write down the rg differential operator . applying @xmath168 to cumulants ( green functions ) amounts to to rg equation ( rge ) @xmath169 with the rg differential operator @xmath170    now , we are in position to discuss the fixed points of the rg . after setting @xmath171 ,
the flow of @xmath155 and @xmath156 is given by the ordinary differential equations @xmath172 with @xmath173 and @xmath174 .
these fixed points , which are approached asymptotically for @xmath107 , follow from the zeros of the gell - mann  low functions , @xmath175 with these functions given to @xmath10-loop order by @xmath176      in addition to the unstable trivial fixed point @xmath178 , there are three more fixed points for @xmath177 .
( 1 ) the galilean - invariant fixed point @xmath179 that corresponds to the gp model , @xmath163
. we will refer to this fixed point as the gp fixed point .
( 2 ) the inversion - invariant fixed point @xmath180 that corresponds to the fis model .
we will refer to this fixed point as the fis fixed point .
( 3 ) the new fixed point @xmath181 we will see below that this is the only stable fixed point .    next , we discuss the stability of the fixed points which is determined by the stability matrix whose components are given by the partial derivatives with respect to @xmath155 and @xmath156 , respectively , @xmath182 its eigenvalues ( the wegner correction exponents ) and the corresponding right eigenvectors which indicate the flow near the respective fixed points are as follows : ( 1 ) for the gp fixed point , we have @xmath183 with @xmath184 and @xmath185 with @xmath186 .
thus , this fixed point is stable only in one direction .
( 2 ) for the fis fixed point , we have @xmath187 with @xmath188 and @xmath189 with @xmath190 . hence
, this fixed point is also stable in only one direction .
( 3 ) for the new fixed point , we have @xmath191 with @xmath192 and @xmath193 with @xmath194 .
therefore , the new fixed point is stable in any direction .
the picture of the rg flow that emerges from the three invariant lines @xmath195 , @xmath163 , and @xmath196 , as well as from the fixed points and their respective stability is sketched in fig .
( [ fig : flussbild ] ) and can be described as follows : the domain of attraction of the stable fixed point is given by the positive quadrant @xmath197 , @xmath198 , that is bounded by the lines @xmath163 ( gp model ) and @xmath195 ( fis model ) .
the unstable fixed points that correspond to the gp and fis models , respectively , are located on these lines . for negative @xmath156
, the rg flow runs off to infinity .
this suggests that we are dealing with a fluctuation - induced first order transition in the case @xmath199 .
of course , the trivial fixed point @xmath178 is stable for @xmath200 .
thus , for dimensions @xmath201 , we find mean - field behavior .
right at the critical dimension , the mean - field power laws have logarithmic corrections .
this case is of particular importance here because for the elastic phase transitions the upper critical dimension and the physical dimension are the same . in the next section
, we will discuss these logarithmic corrections in detail . as a prelude , we study here the asymptotic flow of @xmath155 and @xmath156 to zero for @xmath202 .    for @xmath202 , the flow equations ( [ flow ] ) in conjunction with gell - mann  low functions ( [ beta1l ] ) lead to the differential equations@xmath203 solving these equations for @xmath204
, we obtain@xmath205 the amplitudes are given in general by @xmath206 the special models lead to@xmath207
in this section , we first study the scaling behavior of the cumulants of our theory .
then , we extract from this general result the critical behavior of several experimentally relevant quantities for driven surface diffusion and elastic transitions , respectively .
now , we determine scaling properties of the cumulants @xmath208 . besides the order - parameter - field @xmath209 and its spatial derivatives , we consider the energy - field @xmath210 .
this field is a composite field in the language of renormalization theory with its own scaling behavior .
this scaling behavior can be studied via so - called insertions that are generated by variations @xmath211 of the cumulants with respect to the local temperature@xmath212 .
in other words , the objects of our study here are the expectation values @xmath213^{n}[e]^{m}\rangle ^{(cum)}=\big[-\frac{\delta}{\delta\tau}\big]^{m}g_{n,0 } \ , , \label{kum}\ ] ] which we have written here in a somewhat compressed yet evident form .    to extract the scaling properties of the green functions , we first note that simple dimensional analysis gives @xmath214 for the leading , mean - field behavior .
next , we solve the rge for the green functions to get the anomalous contributions to scaling . for @xmath215 ,
this rge reads @xmath216 as usual , this partial differential equation can be solved by the method of characteristics .
the characteristics read @xmath217 , @xmath218 and @xmath219 , where the @xmath220-factors are the solutions of the ordinary differential equations@xmath221 collecting the solutions of these equations , the flowing couplings @xmath222 that solve eq .
( [ flow ] ) and the mean - field scaling as given in eq .
( [ kumscal ] ) , we obtain the scalings forms @xmath223 this result is correct as long as @xmath215 .    in the case
@xmath224 , @xmath225 or @xmath11 one has to account for additive renormalizations which become essential for the specific heat for @xmath226 .
the additive renormalizations and its consequences are the theme of appendix  [ app : addren ] .
we apply the scalings ( [ kumscal ] ) to the result ( [ g - fin ] ) and obtain@xmath227^{m-1}\ , , \label{enkorr}\ ] ] with @xmath228 given by eq .
( [ i - l ] ) .      in experiments on surface growth ,
one often measures the roughness of a surface and the associated roughness exponents .
these exponents are defined through the height - height correlation function@xmath229 for surface diffusion , the physical dimension is @xmath230 .
hence , we are interested in the scaling behavior of the height - height correlation and its roughening properties for @xmath177 . the asymptotic solutions of the rg - equations ( [ ampldiff ] ) are@xmath231 where @xmath232 , @xmath233 , and @xmath234 are the values of the rg - functions at the fixed point @xmath235 , @xmath236 . for the stable fixed point as well as for the gp fixed point , we have @xmath237 .
hence we get from the definitions ( [ rgfu ] ) the exact relation@xmath238 the scaling behavior of the height - height correlation function  ( [ hh - corr ] ) can be deduced from the general result ( [ genscal ] ) : @xmath239 where we have defined the anisotropy and the correlation length exponents@xmath240 which are related ( with exception of the fis - model ) by@xmath241 due to the presence of strong anisotropy , we have to dealing with two different roughness exponents @xmath242 and @xmath243 defined by    @xmath244    the two exponents    @xmath245    are read off immediately . in the dimensional - expansion about three surface - dimensions ,
we get@xmath246 to first order in @xmath161 .      as mentioned above , the physical and the upper critical dimensions coincide for elastic transitions .
therefore we have @xmath226 , and we are interested in logarithmic corrections .
the typical measurable quantities are the susceptibility @xmath247 of the critical modes @xmath248 , and the critical part of the specific heat @xmath249 that , as the second line eq .
( [ defspez ] ) indicates , can be calculated from the energy correlation function . here
, we need the asymptotic solutions of the rg - equations ( [ ampldiff ] ) for the amplitudes using the assymtotic flow ( [ asympt ] ) .
we get@xmath250 with the constants @xmath251 and @xmath252 given in eqs .
( [ ampl ] ) and ( [ specampl ] ) . from the general scaling properties of the greens - functions ( [ genscal ] ) and bysetting@xmath253 ( we have set the scale so that @xmath254 to simplify the conclusions )
, we obtain @xmath255 asymptotically .
the correction exponent @xmath256 reads @xmath257 whereas the special models lead to@xmath258 next , using the scaling result for the energy correlations ( [ enkorr ] ) we get @xmath259 with the correction exponent @xmath260 .
we recognize that the term originating from the additive renormalization overwhelms the first one asymptotically .
the correction exponent reads @xmath261 the special models lead to@xmath262    note that our results for the correction exponents for the susceptibility and the specific heat of the fis model are in full agreement with those obtained by fis originally .
hence , our calculation passes an another consistency check .
in summary , we have studied driven surface diffusion at the phase transition from isotropic , edwards - wilkinson type behavior to a state characterized by longitudinal ripples and continuous elastic phase transitions .
we have shown that the two transitions belong under appropriate conditions to the same universality class .
we derived the field theoretic hamiltonian for this universality class which has an upper critical dimension of @xmath263 .
the dimension @xmath12 plays different roles in the two systems . in driven surface diffusion
, @xmath12 is the dimension of the substrate whereas for the elastic phase transition , @xmath12 is the full dimension of space . hence , physically ,
@xmath230 in the former problem and @xmath73 in the latter .
thus , the critical behavior in the surface diffusion problem is described by power laws with anomalous scaling exponents whereas that of the elastic phase transitions under consideration is described by power laws with mean - field exponents decorated by logarithmic corrections .
we calculated these anomalous exponents and logarithmic corrections , respectively , for several experimentally relevant quantities .    as far as elastic phase transitions were concerned , we restricted ourselves to crystals with tetragonal or higher symmetry , i.e. , crystals whose harmonic lagrange energy density can be obtained from the tetragonal case by imposing relations between certain elastic constants because crystals with lower symmetry such as orthorhombic or monoclinic symmetry have simple mean field behavior .
in addition to that , we have left aside crystals with cubic symmetry in the transversal subspace . these do have anomalous scaling behavior , however , the diagrammatic perturbation calculation gets significantly more involved when cubic terms are present and is beyond the scope or the present paper .
we will present our results for the cubic case in a future publication  @xcite .
our model hamiltonian differs from the spj model in that we do not assume tilt invariance and that we do assume detailed balance , at least on a semi - macroscopic scale .
our result for the anisotropy exponent @xmath264 is smaller and those for the roughness exponents are larger than the corresponding 1-loop results produced by the spj model .
hence , in comparison , our model leads to a less anisotropic but slightly rougher surface .
its restriction to detailed balance makes our model a special case of a more general model .
our results , moreover , show that our model corresponds to a rg fixed point of this more general model . whether or not this fixed point is the stable fixed point of the more general model
, we can not judge at this point , and we will address this question in future work .
this question is also important from an experimental standpoint .
if our fixed point turns out being stable , this would mean that the asymptotic behavior in a typical experiment on driven surface diffusion with oblique incidence would show detailed balance as characterized by our fixed point .
our model hamiltonian differs from the fis model in that it includes a relevant third order coupling that does not destroy the continuous character of the transition .
the presence of this term modifies the exponents of the logarithmic corrections to the susceptibility and the specific heat . for the former ,
the difference is quantitative , for the latter it is even qualitative .
the fis model predicts that the specific heat diverges as the transition is approached , whereas we find that it vanishes .
as the fis model , our model entirely neglects the noncritical transversal components of the elastic displacement field . if the order parameter field corresponded directly to a component of the strain tensor , one could simply integrate out the noncritical strain components .
this would produce additional negative contributions to the positive ( on grounds of stability ) fourth order coupling constant of our model .
if the additional negative contributions exceed the positive contribution in magnitude , the transition becomes first order .
however , here we have to integrate out the transversal displacement components and not entire components of the strain tensor .
the transversal components of the displacement are contained together with its critical components in different components of the strain tensor , and it is not clear how to integrate them out .
we will return to the question of the influence of the noncritical displacements in future work .    from a conceptual standpoint , it is worthwhile to reflect on the different roles played by reference and target space symmetries in elastic phase transitions and related problems  @xcite .
invariance under rigid rotations in target space , unless broken by external fields , is a symmetry that applies likewise to any crystal , or for that matter to any elastic system that can be described by lagrange elasticity theory , including smectic liquid crystals , liquid crystal elastomers etc .
reference space symmetry , on the other hand , is determined by the crystallographic point group of a given elastic system .
it has been an established approach for setting up field theoretic models for elastic systems for about three decades to employ gradient expansion and then to retain only those terms that are deemed relevant based on dimensional analysis .
examples where this approach has been used include , to name a few , the original work on nematic liquid crystals by gp , isotropic  @xcite and anisotropic  @xcite tethered membranes , sliding columnar phases  @xcite , nematic elastomers  @xcite and their membranes  @xcite and so on , and we also used this approach here . for systems with comparatively high reference space symmetry ,
vestiges of the target space rotational invariance remain in the final field theoretic hamiltonian .
the nonlinear part of the strain @xmath265 in the gp model , for example , ensures that the target space symmetry is realized at least linearly .
the situation is similar in the other examples just mentioned . for the elastic transition studied in the present paper , with exception of the limiting isotropic case that coincides with the gp model ,
the reference space symmetry is comparatively low , and there are no remains of the target space rotational invariance in the final hamiltonian .
it transpires that the extend to which the latter symmetry survives when setting up field theories depends on the former and that , in a sense , reference space symmetry is superior to target space symmetry .
it is desirable to get a deeper understanding of the interplay of the two symmetries in field theories , and we see here an interesting opportunity for future work .
this work was supported in part ( os ) by the national science foundation ( dmr1104707 and mrsec dmr-1120901 ) .
os thanks t.c .
lubensky for spurring his interest in elastic phase transitions and useful discussions .
here , we present a few cornerstones of our calculation of feynman diagrams . the 1-loop diagrams of our theory are depicted in figs .
[ fig:2-leg ] to [ fig:4-leg ] . a typical integral
encountered in calculating these diagrams reads @xmath266 where @xmath267 . in all of our 1-loop integrals
, we can set the external longitudinal momentum to zero . via setting @xmath268 , the integration over the internal longitudinal momentum
can then be performed using the identity @xmath269^{2}(2c)^{2n-1}}\ , .
\label{hilfe}\ ] ] we obtain @xmath270 the remaining integrations are simplified by switching integration variables to @xmath271 and @xmath272 and defining @xmath273 : @xmath274 using dimensional regularization of the integral with @xmath139 , we finally obtain @xmath275 where @xmath276 .
all other feyman integrals at 1-loop order can be calculated by similar means .
we get:@xmath277@xmath278@xmath279 with @xmath280.@xmath281@xmath282    the @xmath11-leg diagrams are shown in fig .
( [ fig:2-leg ] ) . with the results
above we obtain their singular parts that constitute the counter - terms@xmath283@xmath284 and@xmath285    our calculation of the @xmath86-leg diagrams ( fig .
( [ fig:3-leg ] ) ) yields @xmath286 and@xmath287 after symmetrization of the external leg moments . from the @xmath288-leg diagrams , fig .
( [ fig:4-leg ] ) we obtain@xmath289@xmath290 and@xmath291 again after symmetrization .
as pointed out in sec .  [ sec : scalingandexponents ] , the green functions @xmath292 and @xmath293 can not be renormalized entirely by multiplicative renormalization alone . rather , they require extra additive renormalization which leads to the inhomogeneous rges ( @xmath294 ) @xmath295^{m-1}\ , , \label{rge02}\ ] ] where @xmath296 is a function of the invariant coupling constants @xmath155 and @xmath156 only .
to lowest order , the additive renormalization arises from the uv - divergence of the diagram ( 4a ) in fig .
( [ fig:4-leg ] ) with removed external legs and coupling constants .
this divergence is cancelled by the additive counter - term@xmath297 which leads after applying the operator @xmath298 to @xmath299 at lowest order .    to proceed towards its solution ,
let s discus the form of the rge ( [ rge02 ] ) in some more detail .
this rge is a partial differential equation that describes the dependence of @xmath300 upon the variables @xmath301 . as far as the rge is concerned , the spatial coordinate plays solely the role of a constant parameter .
the general form of the rge is a quasi - linear partial differential equation ( pde ) , see _
e.g. _ @xcite@xmath302 to find its solution , we consider @xmath303 as the solution of the implicit equation @xmath304 and write the pde as@xmath305 this homogeneous pde can be solved by employing the characteristics @xmath306 , @xmath307 determined as usual by the ordinary differential equations@xmath308 with @xmath309 , and initial conditions @xmath310 , @xmath311 . since @xmath312 remaines constant along the characteristic flow , the solution is@xmath313 and @xmath307 is of the same functional form as that of @xmath314 as a function of @xmath315 , @xmath316 in our case , the characteristics @xmath317 , @xmath218 , @xmath318 , @xmath319 , @xmath320 are defined by the differential equations ( [ floweq ] ) and ( [ ampldiff ] ) , and their asymptotic solutions in the case @xmath226 , @xmath107 are given in equations ( [ asympt ] ) and ( [ ampllog ] ) .
it remains to solve the differential equation for @xmath307@xmath321 this equation has the solution @xmath322 where @xmath228 is the integral@xmath323 since @xmath324 and @xmath325 , we finally get @xmath326^{m-1 } \label{g - fin}\ ] ] from eq .
( [ u - l ] ) .
janssen , z.  phys .
b * 23 * , 377 ( 1976 ) ; r.  bausch , h.k .
janssen , and h.  wagner , z.  phys .
b * 24 * , 113 ( 1976 ) ; h.k .
janssen , in _ dynamical critical phenomena and related topics _
( _ lecture notes in physics _ , vol .
104 ) , edited by c.p .
enz , ( springer , heidelberg , 1979 ) .


positronium ( ps ) , the bound state of electron and positron , was extensively investigated since its discovery @xcite contributing to the development of bound state qed ( see e.g. @xcite for a review on the current status of this field ) .
furthermore , this system provided stringent limits on possible deviation from the standard model that could indicate new physics @xcite ( see also @xcite for a comprehensive review of former experiments ) . moreover , in the field of materials science ps found various applications due to its unique properties ( see e.g. @xcite for modern reviews on this subject ) .
one recent example is the characterization via positron annihilation lifetime spectroscopy ( pals ) of low - k dielectrics that are potential candidates for the next generation of integrated circuits @xcite .
the motivation of the work presented in this paper is to understand if mesoporous silica with an interconnected pore network could be used for producing a high fraction of positronium at low temperatures .
this would open the door for a new generation of experiments in fundamental research .
cold positronium could be used to improve the precision of spectroscopic studies of ps @xcite or to perform the first spectroscopy of the ps@xmath2 molecule @xcite .
furthermore , it could provide an alternative to the methods that are used for anti - hydrogen formation @xcite-@xcite . as it was suggested sometime ago @xcite-@xcite , anti - hydrogen
could be formed using charge - exchange of ps with anti - protons .
this was demonstrated for the charge conjugate reaction @xcite and in the atrap collaboration resonant charge exchange collisions of positrons with rydberg cs atoms were used to form rydberg ps that via charge exchange with the anti - protons produced anti - hydrogen in rydberg states @xcite .
recently , two experiments @xcite were proposed to perform an anti - gravity test using this process . in both experiments ,
one of the main issues is that a high fraction of ps at low temperatures should be available
. another interesting application would be the possibility to perform an experiment in order to confirm the interpretation of the recent dama / libra annual modulation signal @xcite as generated by mirror type dark matter @xcite .
a step further would be to achieve bose - einstein condensation of ps @xcite .
this would allow the exploration for the first time of the effects of the collective properties of a matter - antimatter system .    in mesoporous silica , ps
is produced by injecting positrons into the film and the distribution of the implantation depth follows a makhovian profile @xcite . in the following we solely consider the long lived triplet spin state ( called ortho - positronium with 142 ns lifetime , o - ps ) because the singlet spin state ( called para - positronium , p - ps ) has a very short lifetime of 125 ps and can be considered as annihilating in the target .
the o - ps ( for simplicity we will refer to it as ps ) that diffuses into the pores loses its kinetic energy via scattering .
if the pores are interconnected , the ps has a probability to tunnel from one pore to another .
a fraction of the ps reaches the film surface and exits into vacuum .
a classical model of thermalization process was developed by nagashima et al .
their calculations reproduce very well the behavior for sio@xmath2 aerogel with pore sizes of about 100 nm . however , a classical approach is not expected to give reliable predictions for ps confined in few nm pores because quantum mechanical effects become relevant . as a matter of fact , in this regime
the de broglie thermal wavelength of ps is comparable with the size of the pores .
recently , mariazzi et al .
@xcite considered phonon scattering to reproduce the thermalization process of ps in a box ( closed porosity ) showing that the minimal energy is not that of the lowest accessible level because the momentum phonons can exchange is fixed . in the same paper , they pointed out that this is not the case in rectangular channels because in one direction the side of the potential well tends to infinity ( z - axis ) . therefore , the magnitudes of the @xmath3 and @xmath4 momentum are quantized but the @xmath5 tends to be a continuum and the minimal energy that ps can reach is given by the ground state in the x - y components .
thus measuring the ps emission energy provides a method to distinguish between the two different pores architectures ( see section [ sec : discussion ] ) .    in previous studies of ps emission in vacuum using time of flight ( tof )
, many interesting effects , such as the emission from the surface of different materials @xcite-@xcite were investigated .
recently , this technique was applied to study mesoporous and hybrid silica films @xcite-@xcite , to evaluate the continuity barrier @xcite and the effect of thermalization for pore surfaces decorated with different groups @xcite .
however , the influence of the temperature on ps emission in vacuum from mesoporous thin films was never studied in detail .
to our knowledge , only a work of mills et al .
@xcite with fumed silica revealed a dependence with the temperature of the ps emission in vacuum .
the fraction of ps in the low energy tail was estimated but no quantitative estimate on the value of the minimal emission energy was obtained .
very recently cassidy et al .
@xcite measured the emission energy of ps from mesoporous films using doppler spectroscopy .
this is a different technique from that used in the present study .
results consistent with the one presented in this paper ( though not sample temperature dependencies ) have been obtained .
in this paper , we study the ps yield and emission energy in vacuum as a function of the target temperature for two different kinds of mesoporous thin films with the same tetraethoxysilane ( teos ) mineral source for the silica network skeleton precursor : ctacl - teos and f127-teos .
the density of the c sample is approximately 1.2 g/@xmath6 and of the f sample is 1.5 g/@xmath6 .
both samples were spin - coated on glass similar to the ones we measured in @xcite .
the c samples are prepared via a sol - gel process using cetyl trimethyl ammonium chloride ( ctacl ) cationic surfactants as the organic pore generator ( porogen ) agent @xcite .
a pure aqueous method is used .
the ctacl / teos molar ratio for the films prepared is 0.22 . after deposition ,
the ctacl - teos / glass samples are treated at 130 @xmath7c and stored in air .
the f samples use non - ionic pluronic f-127 triblock copolymer ( eo@xmath8po@xmath9eo@xmath8 ) as surfactant and were prepared in the same way as described in @xcite .
both samples were calcinated for 15 minutes at 450 @xmath7c in air immediately before the e@xmath10 measurements .
the recorded x - ray diffraction patterns indicate no symmetry in the pore organization .
the ethz slow positron beam used for these measurements is described in greater detail in @xcite .
the positrons flux is 25000 e@xmath10/s .
the slow positron beam is stopped in the sio@xmath2 target .
the positrons can either form positronium , i.e. o - ps or p - ps , or annihilate into @xmath11 s .
the detection with an mcp ( micro - channel plate ) of the secondary electrons ( se ) emitted when the positrons hit the target serves for tagging the positronium formation .
the se leave the target accelerated to 1 - 11 kev by the same voltage applied to the target relative to the grounded transport tube that is used to implant the positrons in the positronium converter .
the se are then transported by a magnetic field in the backward direction , as shown in figure [ traj ] .
the electrons move along the magnetic field line in spirals and are deflected to the mcp region by the @xmath12 filter .
the tagging efficiency varies from 70% to 30% in the energy range 1 - 10 kev .
the samples are mounted on a cryocooler head to allow the possibility of varying their temperature in the range of 50 - 400 k.      the start time t@xmath13 for the detectors is triggered by the mcp detecting the se emitted when the positrons hit the target . in both pals and tof detectors
the stop is given by one ( or more ) annihilation photons depositing some energy in the calorimeter ( ecal ) .
both ecal are composed of bgo crystals with hexagonal shape , 61 mm external diameter and 200 mm length .
the time resolution of the system mcp - ecal was measured to be around 5 ns ( fwhm ) .
the typical energy resolution of the crystals is about 25 - 30@xmath14 ( fwhm ) .
the ecal for the tof is placed behind a lead slit at a distance z from the target that can be varied ( as sketched in figure [ traj ] ) .
the detector is screened by four half cylinders of lead surrounding the beam pipe .
the thickness of the shielding is 70 mm , the width of the slit is set to 5 mm and its position with respect to the target can be adjusted . for this measurement in order to maximize the signal to background ratio the center of the slit was placed at a distance of 18 mm from the target .
the main contribution to the background is given by photons coming from direct positrons and p - ps annihilations in the target ( so called prompt peak ) .
some of those photons can be detected in the crystals after compton scattering in the lead shield .
the determination of the slit position with respect to the target was done by scanning the slit position in 0.1 mm steps and recording the maximum of the 511 kev annihilation peak .
the scans were performed 30 min after the temperature was set on the sample .
this technique was very important in order to correct the slit position as a function of the temperature .
in fact , at 50 k the contraction of the cryocooler head was measured to be @xmath15 mm in agreement with the prediction of finite elements calculations performed with the comsol package @xcite .
the pals detector was designed to have a large acceptance to provide a uniform efficiency for detecting the ps emitted in vacuum @xcite .
in order to design the detectors and interpret the data , simulations served as a powerful tool . in the monte carlo
( mc ) simulations of our setup , the 3d eb - fields were calculated with the comsol multi - physics program and the positron / electron trajectories in the beam were simulated with geant4 .
the simulation of the photon detection in the apparatus was based on the same package @xcite .
new classes were written in order to simulate the o - ps production , propagation in the beam pipe , reflection on pipe walls that were assumed to be knudsen - like ( ps is reflected isotropically ) , pick - off effect or decay ( ps is not a standard particle in geant4 ) .
the events for the o - ps@xmath16 process were generated taking into account the decay matrix element @xcite .
the geometries of the beam transport pipe , photon detector , positron tagging system and its material were coded into simulations .
the results were cross - checked with our experimental measurements for both photon detection @xcite and particle transport in the eb - fields @xcite .
the measurements were taken in a clean vacuum of @xmath17 mbar . to avoid water contamination of the film during the cooling down ,
the target was kept at room temperature for an hour using a heater before lowering its temperature .
the cooling cycles were repeated and the data confirmed the reproducibility of the results . the signals from the bgo s photomultipliers are split to record both energy and timing with a charge - to - digital converter ( qdc caen v792 ) and a time - to - digital converter ( tdc caen v775 ) .
a cut on the energy deposited in the bgo between @xmath18 kev was applied to optimize the signal to background ratio suppressing compton scattering events in the collimator from direct and p - ps annihilations in the target .      in this section
we present the results we obtained from the pals measurements .
the spectra are analyzed using the lt9 @xcite program .
the lifetimes of the decay components and its fractions are resolved by fitting the pals spectra .
the program finds 3 exponentials convoluted with a 5 ns fwhm resolution function of the spectrometer .
the shortest exponential ( less than 4 ns ) is originated by direct positron and p - ps annihilations .
it is disregarded because we are interested only in contribution of o - ps .
we define as ( @xmath19 , @xmath20 ) , and ( @xmath21 , @xmath22 ) the intensities and the lifetimes of the two longer exponentials . to determine the yield @xmath23 of ps emitted in vacuum and the lifetime @xmath24 in the pores of the film we used a model of ps escaping in vacuum @xcite . according to this model the lifetime @xmath24 in the pores of the film
is defined as : @xmath25 the yield @xmath23 of ps emitted in vacuum was calculated according to : @xmath26 where @xmath27 is the escape rate of free ps into vacuum .
we extracted the films thicknesses for the two samples by fitting the total ps yield as a function of the positron implantation energy to the makhovian profile .
these thicknesses are reported in table [ tbl1 ] in which we also present the values of @xmath24 and @xmath28 calculated with the expressions above .
we use the results of the fits of the pals spectra at 6 kev for the c sample and 10 kev for the f sample at 50 k and 300 k. these implantation energies will serve as a reference for the rest of the paper . we choose those values in order to maximize the amount of thermalized ps . at these energies the majority of the positrons are still implanted within the films ( no significant drop of the total ps yield @xmath29 is observed , see fig .
[ palsdata ] ) and the emission energy is in the constant region ( see next section ) . with @xmath24 one can calculate the effective pore size @xmath30 in the films applying the gidley et al .
@xcite extension of the tau - eldrup model @xcite-@xcite ( we will call it hereafter rte model ) . as one can see , at lower temperature the lifetime in the film increases as predicted by the rte model .
this can be understood in the following way : the overlap of the ps wave function with the volume contained within a distance @xmath31 nm @xcite from the walls for which the annihilation rate is assumed to increase is less for ps confined in the pores occupying the ground state than for ps in excited states .
since at lower temperatures the population of the ground state is higher the pick - off rate decreases . however , the measured lifetimes are lower than what is expected by the calculations using the rte model . for the c sample the discrepancy is less than 10% but for the f sample the measured value is 30% lower suggesting that the pick - off rate increases .
deviations from the rte model were already observed in previous measurements using the 2 to 3@xmath32 ratio technique @xcite and in a measurement using pals @xcite .
different factors could be responsible for that as pointed out in those papers and deserves further studies .
since the parameter @xmath33 of the rte model was calibrated at room temperature we use the lifetimes at 300 k to determine the effective pore size that ps experiences .
we report the results at the end of the paper in table [ tbl3 ] for both rectangular channels and cubic boxes .
interestingly , we found that the yield of ps emitted in vacuum can be considered as independent of the temperature of the target ( see fig .
[ palsdata ] ) .    ) and total yield of ps ( @xmath34 ) for f ( upper plot ) and c samples ( lower plot ) at room temperature and 50 k.__,scaledwidth=50.0% ]    this may seem to contradict the fact that the lifetime in the films increases with the temperature .
however , our measurements indicate that the escape rate in vacuum is smaller at lower temperature ( see table [ tbl1 ] ) , explaining the observation that the yield is constant .
a possible explanation of this effect can be found considering that ps tunnels from one pore to the other . in this case , the tunneling probability will decrease with temperature explaining our measurements .
.[tbl1 ] film thickness ( z ) , ps lifetime in the pores @xmath35 and escape rate @xmath28 at 50 k and 300 k. [ cols="^,^,^,^,^,^",options="header " , ]      figure [ tof_fd ] shows the acquired tof spectra at different implantation energies for the f sample . for every run the time spectra of each crystal
are calibrated finding the position of the peak arising from the annihilations in the target that defines t@xmath13 .
first we extract the mean energy of ps emission in vacuum applying the analysis method proposed in @xcite as follows .
the background due to the annihilations in the target is subtracted by fitting the measured time spectra with the resolution function of our detector determined using targets ( aluminum and kapton ) in which the ps formation is negligible .
the tof spectra are corrected for the ps decays and the time spent in front of the detector with the factor @xmath36 .
the maximum of the peak distribution defines the mean ps emission energy in the direction perpendicular to the film surface .
let us note that , neglecting the reflection of ps in the beam pipe , with the time - of - flight method one measures only the mean of the energy component perpendicular to the collimator that we define as the z - axis ( we define it as @xmath37 ) .
the triangles in fig .
[ opsevsv ] , represent it as a function of the implantation energy .
as a function of the implantation voltage at a target temperature of 300 k. the triangles represents the energy extracted with the maximum analysis method ( ma ) , the squares are after the subtraction of the non - thermalized part ( maximum analysis method after correction of the spectra , ma - cs).__,title="fig:",scaledwidth=50.0% ]   as a function of the implantation voltage at a target temperature of 300 k. the triangles represents the energy extracted with the maximum analysis method ( ma ) , the squares are after the subtraction of the non - thermalized part ( maximum analysis method after correction of the spectra , ma - cs).__,title="fig:",scaledwidth=50.0% ]    one can see that starting from 3 kev for the c sample ( this was confirmed by other measurements on similar samples @xcite ) and 4 kev for the f sample the value of the energy emitted in vacuum tends to be constant . clearly , the emission energy of ps calculated in this way is not the minimal energy due to the confinement in the pores because one has different contributions given by the convolution of the emission energy with the implantation profile . in order to isolate the thermalized part from the tof spectra , one has to subtract the contributions of non - thermalized ps components ( see fig . [ components ] ) . to estimate these non - thermalized contributions , we suppose that the shape of their distribution , nt(t ) , is represented by the tof spectrum obtained at a low implantation energy . to select the implantation energy ( from the ones we had measured ) and the scaling factor of nt(t ) for the subtraction , we relied on the mc .
we used the values for which the best fit between the mc and the spectra obtained after the subtraction of nt(t ) was achieved ( for more details see @xcite ) .
we found that the best fits were obtained for 2 kev with l@xmath38=200 nm in the case of the c sample , and 3 kev and l@xmath39=350 nm for the f sample .
for a given implantation energy e@xmath40 , nt(t ) is scaled down by the fraction of positrons implanted at depths smaller than l@xmath38 and l@xmath39 .
this was determined by using a makhovian profile .
the results of this analysis are shown as the squares in fig .
[ opsevsv ] .
a parabolic fit is used to determine the position of the maximum .
the statistical error of the fit is typically @xmath419 ns for the f and @xmath416 ns for the c sample . the uncertainty on the determination of the slit position of @xmath42 mm results in a systematic error of the order of @xmath43 mev in the determination of the ps mean emission energy @xmath37 ( for implantation energies above 3 kev ) .
the subtraction procedure described above introduces a systematic error that we estimated analyzing the data using different values of l@xmath44 nm and l@xmath45 nm .
the estimated error is @xmath46 mev for the c and @xmath47 mev for the f sample .
thus the combined statistical and systematic error is at a level of @xmath412.9 mev for the f and @xmath413.0 mev for the c sample .    the results for the c and the f samples at room temperature and at 50 k are shown in figs .
[ opsevsv_rt]-[opsevsv_lt ] .
for implantation energies higher than 4 kev for the c and 5 kev for the f sample the values @xmath37 of the mean emission energy are constant . in the c sample
, @xmath37 reaches its constant value at lower implantation voltages because the pore size is smaller than in the f sample ( see fig .
[ opsevsv_rt ] ) .
those values are higher than the thermal energy that ps will have if it would thermalize at the temperature of the film .
as expected in the presence of confinement in the pores , the mean emission energy is higher for the c sample with pore sizes smaller than the f sample .     as a function of the positron implantation energy for c and f samples at 300 k.__,scaledwidth=50.0% ]     for implantation energies
higher than 2 kev at 50 k and 300 k for the c sample .
lower plot : ps mean energy as a function of the implantation energies higher than 4 kev at 50 k and 300 k for the f sample.__,scaledwidth=50.0% ]    the tof technique measures @xmath37 , the mean energy of the ps atoms in the z - direction . to find the mean emission energy of ps in vacuum , @xmath37 should be multiplied by a factor @xmath48 that takes into account the angular distribution .
assuming that the ps is emitted mono - energetically and isotropically from the surface , one can calculate that @xmath48 is equal 2 . in this estimation , the reflection of ps in the beam pipe and the detector acceptance , i.e. the fact that a fraction of events decaying before or after the collimator aperture are detected , are not taken into account .
therefore , to determine @xmath48 considering these effects we used the mc simulation we described in section [ sec : mc ] .
as shown in the upper plot of fig .
[ mc ] , a satisfactory agreement between the data and the mc ( adding the spectra at 2 kev that takes into account the non - thermalized ps ) is achieved .
we attribute the difference between 40 and 100 ns to the approximation used in the subtraction method where the contribution of the non - thermalized ps is underestimated since only a spectrum of a defined energy is used for this correction .
the fact that the ps is emitted with an angular spread is clearly supported by the data .
as one can see in the lower plot of fig .
[ mc ] for ps emitted with no angular spread the data are not reproduced .
the physical interpretation is that in the films we studied the pores have no organization thus they are expected to be randomly aligned .
the value of @xmath48 estimated with the mc is 1.7 .
this is consistent with the expectation of the analytical result and the values reported in previous experiments @xcite .
a detailed scan shows that the mean ps energy ( @xmath49 ) decreases with the sample temperature down to a minimum level ( see fig .
[ opsenergyvst ] ) . for the c sample
this value is basically constant ( 73@xmath415 mev ) in the range of temperature in which we performed our measurement .
this can be understood by the fact that in this sample the confinement energy is much higher than the thermal energy at room temperature ( @xmath50 mev ) thus almost all ps is in the ground state . for the f sample
there is a weak dependence on the temperature . due to the bigger pore size compared to the c sample ,
the energy of the ground state is only twice the thermal energy at room temperature .
therefore , the probability to find the ps occupying an excited state is higher .
as expected , this probability decreases with the temperature thus the minimal energy reaches its constant value of 48@xmath415 mev .
the time that ps spends in the films before being emitted in vacuum was not considered in our determination of the emission energy .
the measurements presented in fig .
[ opsenergyvst ] are in a regime in which a classical approach is not expected to give reliable results .
some theoretical work to develop a full quantum mechanical picture of the emission process is required to address this problem ( as pointed out in @xcite as well ) .      to understand the behavior of the value of the minimal energy as a function of the film temperature we present a simple model of ps in thermodynamic equilibrium at a temperature t in rectangular boxes .
the expectation value @xmath51 of the hamiltonian operator for ps confined in 1 dimensional infinite well in contact with a reservoir at a temperature @xmath52 is given by : @xmath53 where z is the partition function defined as @xmath54 where @xmath30 is the dimension of the well , @xmath55 is the ps mass , @xmath56 is the principal quantum number and @xmath57 and @xmath58 are the planck and the boltzmann constants . to calculate the mean value @xmath59 of the energy for the 3d case we can use @xmath60 where to calculate @xmath61 and @xmath62 one can substitute the pore side length @xmath30 in eq .
[ eq : z ] with @xmath63 and @xmath64 . for ps confined in rectangular pores ,
we thus obtain : @xmath65 for the case of a cubic box one can set @xmath66 .    to compare the prediction of the pore size that one can extract from the tof measurements with the pals results we fit the data using eq .
[ eq : h ] . to construct the function used for the fit we kept only the first 50 terms of the sum .
this is very conservative , the probability for ps to occupy a state higher than n@xmath6710 for the kind of target and the temperatures we used in this study is already negligible .
the solid lines in fig .
[ opsenergyvst ] represent the fit to the data where the pore side length ( @xmath30,@xmath63,@xmath64 ) are left as free parameters .
the fits were repeated assuming cubic box pores and the results are shown as the dashed lines in fig .
[ opsenergyvst ] .
we used migrad from the minuit package @xcite as a minimization procedure .
as one can see , the fit to the data suggests that the pores of both samples are better modeled as rectangular pores than cubic boxes . as proposed in @xcite , to compare the values obtained from the fit with the ones extracted from the pals measurements , one has to add twice the parameter @xmath33 ( see previous section ) .
the pore side lengths obtained in this way are reported in table [ tbl3 ] .
in particular , the fit supports the idea that the pores are better modeled by rectangular channels .
the pore side length in one direction ( @xmath64 ) obtained from the fit is much longer than the substrate thickness . since
, the same values for the side length @xmath30 and @xmath63 are obtained fitting the data assuming rectangular channels instead of rectangular boxes we do not report this value in table [ tbl3 ] as it has no physical meaning .     with the pore side lengths @xmath30,@xmath63,@xmath64 left as free parameters .
the dashed lines were obtained fitting with eq .
[ eq : h ] with a single side length free @xmath68 ( cubic box pores).__,scaledwidth=50.0% ]    ccc + & tof & pals + @xmath69 [ nm ] & 3.3@xmath410.1 & 4.2@xmath410.5 + @xmath70 [ nm ] & ( 2.7@xmath410.8 , 2.7@xmath410.8 ) & ( 3.3@xmath410.5 , 3.3@xmath410.5 ) + @xmath71 [ mev ] & 73@xmath415 +   + & tof & pals + @xmath69 [ nm ] & 4.1@xmath410.1 & 6.4@xmath410.5 + @xmath70 [ nm ] & ( 2.9@xmath410.8 , 3.6@xmath411.8 ) & ( 4.3@xmath410.5 , 5.3@xmath410.5 ) + @xmath71 [ mev ] & 48@xmath415 +    according to the quantum mechanical model for ps thermalization of mariazzi et al .
@xcite the index @xmath56 in the sum of eq .
[ eq : z ] could differ from 1 for cubic pores if the level separation of two close ps energy levels is higher than the maximum momentum that a single phonon can exchange .
we performed fits with different @xmath72 but the results did not improve , supporting the idea that the pores are not very well modeled by cubic boxes . to summarize , for the c sample with the tof data the best fit was obtained for pores modeled as square channels of 2.7@xmath410.8 nm side length while the lifetime method gives 3.3@xmath410.5 nm .
the best fit to the data for the f sample was obtained for rectangular channels with a cross section of 2.9@xmath410.8 nm by 3.6@xmath411.8 nm that has to be compared with 4.3@xmath410.5 by 5.3@xmath410.5 of the pals measurement .
this pore size was extracted by applying the rte model to reproduce the measured lifetime of 74 ns assuming the ratio @xmath73 to be the same as in the tof results .
both values obtained with the tof measurement are systematically lower than the pals results . nevertheless , considering the approximation used in our model , the assumptions that the pores can be treated as rectangular channels and the uncertainty in the determination of the pore size we conclude that our results ( summarized in table [ tbl3 ] ) are in reasonable agreement with expectations .
in this paper , we show that the yield of ps emitted in vacuum measured with the pals technique is independent of the temperature of the mesoporous thin films .
the lifetime in the films increases with a decrease in temperature but the measured values are lower than expected by the rte model .
this suggests that another source of pick - off which depends on the temperature should be invoked .
further studies are needed to investigate the origin of this effect .
the escape rate in vacuum decreases explaining the observation that the yield of ps emitted in vacuum is the same at 50 and 300 k.    furthermore , we show that due to quantum mechanical confinement in the pores the ps emission energy into vacuum has a minimal value .
the minimal energy is higher for the sample with smaller pores and this constant value is reached at a lower positron implantation energy .
our results are in fair agreement with a model of ps confined in rectangular channels in thermal equilibrium with the sample .
the measured minimal energy for the c sample was found to be 73@xmath415 mev while for the f sample it is 48@xmath415 mev .
these experimental results provide a solid ground to understand how to produce ps at lower temperature and could serve to develop a more realistic model to interpret the data .
this work was supported by the swiss national foundation , the eth zurich and the cnpq ( brazil ) .
p.c . would like to thank c. lenz cesar , p. lotti and b. barbiellini for the very useful discussions .
m. deutsch , phys .
82 , 455 ( 1951 ) .
s. karshenboim , int . j. mod
. phys . a 19 , 3879 ( 2004 ) . a. badertscher , p. crivelli , w. fetscher , u. gendotti , s. gninenko , v. postoev , a. rubbia , v. samoylenko , d. sillou , phys . rev .
d 75 , 032004 ( 2007 ) [ arxiv : hep - ex/0609059 ] .
a. rubbia , int .
phys . a 19 , 3961 ( 2004 ) [ arxiv : hep - ph/0402151 ] . s.w.h .
eijt , b. barbiellini , a.j .
houtepen , d. vanmaekelbergh , p.e .
mijnarends , a. bansil , phys .
( c ) 4 , 3883 ( 2007 ) , arxiv:0705.1113v1 .
d. w. gidley , h. g. peng , r. s. vallery annual review of materials research 36 , 49 ( 2006 ) .
d. w. gidley , w. e. frieze , t. l. dull , a. f. yee , e. t. ryan , and h. m. ho , phys . rev .
b 60 , r5157 ( 1999 ) .
d. w. gidley , w. e. frieze , t. l. dull , j. sun , a. f. yee , c. v. nguyen , and d. y. yoon , appl .
76 , 1282 ( 2000 ) .
s.  chu , a.  p.  mills , j. l. hall , phys .
52 , 1689 ( 1984 ) .
m.  s.  fee , a.  p.  mills , s.  chu , e.  d.  shaw , k.  danzmann , r.  j.  chichester and d.  m.  zuckerman , phys .
70 , 1397 ( 1993 ) .
d. b. cassidy , h. w. k. tom and a. p. mills , aip conf .
1037 , 66 ( 2008 ) .
m. amoretti , et al .
[ athena collaboration ] , nature ( london ) 419 , 456 ( 2002 ) .
g. gabrielse , n. bowden , p. oxley , a. speck , c. h. storry , j. n. tan , m. wessels , d. grzonka , w. oelert , g. schepers , t. sefzick , j. walz , h. pittner , t. w. hnsch , and e. a. hessels , phys .
89 , 213401 ( 2002 ) .
g. b. andresen , [ alpha collaboration ] , phys .
b 685 , 141 ( 2010 ) .
j. w. humberston , m.  charlton , f.  m.  jacobsen and b.  i.  deutch , j. phys .
b 20 , l25 ( 1987 ) .
m. charlton , phys .
a 143 , 143 ( 1990 ) . b. i. deutch , m. charlton , m. h. holzscheiter , p. hvelplund , l. v. jorgensen , h. knudsen , g. laricchia , j. p. merrison and m. r. poulsen , hyperfine interaction 76 , 153 ( 1993 ) .
j. p. merrison , h. bluhme , j. chevallier , b. i. deutch , p. hvelplund , l. v. jorgensen , h. knudsen , and m. r. poulsen , m. charlton , phys .
78 , 2728 ( 1996 ) . c. h. storry , a. speck , d. le sage , n. guise , g. gabrielse , d. grzonka , w. oelert , g. schepers , t. sefzick , h. pittner , m. herrmann , j. walz , t. w. hnsch , d. comeau , e. a. hessels , phys .
93 , 263401 ( 2004 ) .
a. kellerbauer et al .
[ aegis collaboration ] , nucl .
instr . and meth .
b 266 , 351 ( 2008 ) .
p. perez , l. liszkay , m. f. barthe , j .-
boilot , m. carty , c. corbel , p. crivelli , a. curtoni , o. delferrierre , p. desgardin , m. etienne , u. gendotti , m. hassan , a. mohri , t. ohdaira , j .-
m . rey , a. rubbia , h. saitoh , y. sauce , r. suzuki , a. walcarius , and y. yamazaki , aip conf .
1037 , 35 ( 2008 ) .
r. bernabei et al .
[ dama collaboration ] , eur .
j. c 56 , 333 ( 2008 ) .
[ arxiv:0804.2741 [ astro - ph ] ] .
r. foot , phys .
d 78 , 043529 ( 2008 ) [ arxiv:0804.4518 [ hep - ph ] ] .
a.  badertscher , a.  belov , p.  crivelli , s.  n.  gninenko , j.  p.  peigneux , a.  rubbia and d.  sillou , int .
j.  mod .
phys .  a 19 , 3833 ( 2004 ) . p. m. platzmann and a. p. mills jr . , phys
b 49 , 454 ( 1994 ) .
s. valkealahti and r. nieminen , appl .
phys . a 32 , 95 ( 1983 ) .
m. tuomisaari , r. h. howell , and t. mcmullen , phys .
b 40 , 2060 ( 1989 ) .
y. nagashima et al . , phys . rev .
b 58 , 12676 ( 1998 ) .
a. p. mills , jr .
, e. d. shaw , r. j. chichester , and d. m. zuckerman , phys . rev .
b 40 , 2045 ( 1989 ) .
r. s. yu , t. ohdaira , r. suzuki , k. ito , k. hirata , k. sato , y. kobayashi , and j. xu , appl .
83 , 4966 ( 2003 ) .
k. ito , r. s. yu , k. sato , k. hirata , y. kobayashi , t. kurihara , m. egami , h.arao , a. nakashima , and m. komatsu , j. appl .
98 , 094307 ( 2005 ) . c. he , m. muramatsu , n. oshima , t. ohdaira , a. kinomura , and r. suzuki , phys .
lett . a 355 , 73 ( 2006 ) .
m. p. petkov , m. h. weber , k. g. lynn , and k. p. rodbell , appl .
79 , 3884 ( 2001 ) .
r. suzuki , t. ohdaira , y. kobayashi , k. ito , r. s. yu , y. shioya , h. ichikawa , h. hosomi , k. ishikiriyama , h. shirataki , s. matsuno , and j. xu , mater .
forum 445 - 446 , 224 ( 2004 ) .
j. xu , a. p. mills , jr . , and c. case , appl .
87 , 54105 ( 2005 ) .
h. k. m. tanaka , t. kurihara , and a. p. mills , jr .
b 72 , 193408 ( 2005 ) . c. he ,
t. ohdaira , m. muramatsu , n. oshima , a. kinomura , r. suzuki , t. oka and y. kobayashi phys .
b 75 , 195404 ( 2007 ) .
d. b. cassidy , p. crivelli , t. h. hisakado , l. liszkay , v. e. meligne , p. perez , h. w. k. tom , and a. p. mills jr .
a 81 , 012715 ( 2010 ) .
l. liszkay , c. corbel , p. perez , p. desgardin , m .- f .
barthe , t. ohdaira , r. suzuki , p. crivelli , u. gendotti , a. rubbia , m. etienne , and a. walcarius , appl .
lett . , 92 , 063114 ( 2008 ) .
y. cohen , k. landskron , n. tetretault , s. fournier - bidoz , b. hatton , g. a. ozin , adv .
15 , 593 ( 2005 ) .
w. yantasee , y. lin , x. li , g. e. fryxell , t. s. zemanian and v. v. viswanathan , analyst , 128 , 899 ( 2003 ) .
n. alberola , t. anthonioz , a. badertscher , c. bas , a.s .
belov , p. crivelli , s.n .
gninenko , n.a .
golubev , m.m .
kirsanov , a. rubbia and d. sillou , nucl .
instr . method a 560 , 224 ( 2006 ) .
p. crivelli , u. gendotti , l. liszkay , p. perez , a. rubbia , mat .
forum 607 , 251 ( 2009 ) .
agostinelli et al [ geant4 collaboration ] , nucl .
instr . method a 506 , 250 - 303 ( 2003 ) .
a. ore and j. l. powell , phys .
75 , 1696 ( 1949 ) .
a. badertscher , p. crivelli , m. felcini , s.n .
gninenko , n.a .
goloubev , p. nedelec , j.p .
peigneux , v.postoev , a. rubbia and d. sillou phys .
p. crivelli ,  a new search for invisible decays of ortho - positronium " , phd thesis , no.16617 , eth zurich , switzerland ( 2006 ) .
l. liszkay , c. corbel , l. raboin , j .-
boilot , p. perez , a. brunet - bruneau , p. crivelli , u. gendotti , a. rubbia , t. ohdaira and r. suzuki , appl .
95 , 124103 ( 2009 ) .
j. kansy , nucl .
instrum . meth . a 374 , 235 ( 1996 ) .
l. dull , w. e. frieze , d. w. gidley , j. sun , and a. f. yee , j. phys .
b 105 , 4657 ( 2001 ) .
s. j. tau , chem .
56 , 5499 , ( 1972 ) .
m. eldrup , d. lightbody amd j. n. sherwood , chem .
phys . lett .
63 , 51 , ( 1981 ) .
t. goworek , k. ciesielski , b. jasinska , and j. wawryszczuk , chem .
272 , 91 ( 1997 ) .
k. ito , h. nakanishi and y. ujihira j. phys .
b 103 , 4555 - 4558 ( 1999 ) .
s. mariazzi , l. toniutti , n. patel , r.s .
brusa , applied surface science 255 , 191 - 193 ( 2008 ) . c. g. fischer , m. h. weber , c. l. wang , s. p. mcneil , and k. g.lynn , phys . rev .
b 71 , 180102r ( 2005 ) .
s. thraenert , e. m. hassan , d. enke , d. fuerst and r. krause - rehberg , phys .
status solidi c 4 , 3819 ( 2007 ) .
u. gendotti , ph.d thesis unpublished , eth zurich , switzerland ( 2010 ) .
f. james , m. roos , minuit , function minimization and error analysis , cern d506 ( long writeup ) .