the presented work takes place in a more general framework : the modelling of unsteady mixed flows in any kind of closed pipe taking into account the cavitation problem and air entrapment .
we are interested in flows occurring in closed pipes with non uniform sections , where some parts of the flow can be free surface ( it means that only a part of the pipe is filled ) and other parts are pressurized ( it means that the pipe is full ) .
the transition phenomenon between the two types of flows occurs in many situations such as storm sewers , waste or supply pipes in hydroelectric installations .
it can be induced by sudden changes in the boundary conditions as failure pumping . during this process
, the pressure can reach severe values and may cause damages .
the simulation of such a phenomenon is thus a major challenge and a great amount of works was devoted to it these last years ( see @xcite,@xcite,@xcite,@xcite , for instance ) .
the classical shallow water equations are commonly used to describe free surface flows in open channels .
they are also used in the study of mixed flows using the preissman slot artefact ( see for example @xcite ) . however
, this technic does not take into account the depressurisation phenomenon which occurs during a water hammer except in recent works @xcite .
on the other hand the allievi equations , commonly used to describe pressurized flows , are written in a non - conservative form which is not well adapted to a natural coupling with the shallow water equations .    a model for the unsteady mixed water flows in closed pipes and a finite volume discretisation have been previously studied by two of the authors @xcite and a kinetic formulation has been proposed in @xcite .
we propose here the * pfs*-model which tends to extend naturally the work in @xcite in the case of a closed pipe with non uniform section .
for the sake of simplicity , we do not deal with the deformation of the domain induced by the change of pressure .
we will consider only an infinitely rigid pipe .
the paper is organized as follows .
section 2 is devoted to the derivation of the free surface model from the @xmath0d _ incompressible _ euler equations which are written in a suitable local reference frame ( following @xcite ) in order to take into account the local effects produced by the changes of section and the slope variation .
the construction of the free surface model is done by a formal asymptotic analysis .
seeking for an approximation at first order gives the model called * fs*-model . in section  [ sectionpfmodel ] , we adapt the derivation of the * fs*-model to derive the pressurized model , called * p*-model , from the @xmath0d _ compressible _ euler equations .
writing the source terms of these two models , * p * and * fs*-model , into a unified form and using the same couple of conservative unknowns as in @xcite , we propose in section [ sectioncoupling ] a model for mixed flows , that we call * pfs*-model .
we state some mathematical properties of this model .
section [ pfsroe ] is devoted to the extension of the vfroe scheme described in @xcite that was used for the case of uniform pipes . in section [ sectionwb ]
, we show how to construct a convection matrix in order to get an exactly well - balanced scheme .
several numerical tests are presented in section [ numericaltests ] .
the classical shallow water equations are used to describe physical situations like rivers , coastal domains , oceans and sedimentation problems .
these equations are obtained from the incompressible euler system ( see e.g. @xcite ) or from the incompressible navier - stokes system ( see for instance @xcite ) by several techniques ( e.g. by direct integration or asymptotic analysis ) .
we adapt here the derivation in @xcite to get a new unidirectional shallow water model .
we start from the @xmath0d incompressible euler equations where we neglect the acceleration following the @xmath1-axis supposing the existence of a privileged main flow axis .
we write then the euler equations in the local serret - frenet reference frame in order to take into account the local effects produced by the changes of section and the slope variation .
then we derive a shallow water model by a formal asymptotic analysis ( done in subsection [ sectionfsmodelasymptotic ] ) .
let us consider the cartesian reference frame @xmath2 .
in the corresponding coordinate system @xmath3 , the @xmath0d incompressible euler system writes : @xmath4 where @xmath5 denotes the velocity with components @xmath6 , @xmath7 is the isotropic pressure tensor , @xmath8 the density of the fluid at atmospheric pressure @xmath9 and @xmath10 is the exterior strength of gravity .
we close classically system  ( [ e3dicartesian ] ) using a kinematic law for the evolution of the free surface : _ any free surface particle is advected by the fluid velocity _ @xmath11 and on the wet boundary , we assume the no - leak condition @xmath12 where @xmath13 is the outward unit normal vector to the wet part of the pipe ( see f. [ supplementaire ] ) .
we set the pressure @xmath14 to @xmath15 at the free surface .
we define the domain @xmath16 of the flow at time @xmath17 as the union of sections @xmath18 ( assumed to be simply connected compact sets ) orthogonal to some plane curve @xmath19 lying in @xmath20 to follow the privileged main flow axis .
we choose the parametrization @xmath21 in the cartesian reference frame @xmath2 where @xmath22 follows the vertical direction ; @xmath23 is then the elevation of the point @xmath24 over the plane @xmath25 ( see f. [ oxoz ] ) .
we define a local reference frame as follows : we introduce the curvilinear variable defined by : @xmath26 where @xmath27 is an arbitrary abscissa .
we set @xmath28 and we denote by @xmath29 the @xmath30-coordinate of any fluid particle @xmath31 in the serret - frenet reference frame @xmath32 at point @xmath24 with @xmath33 the tangent vector @xmath34 , the normal and @xmath30 the binormal vector ( see f. [ oxoz ] and f. [ oyoz ] for the notations ) .
@xmath30 is normal to @xmath19 in the vertical plane @xmath20 .
then , at each point @xmath35 , @xmath36 is defined by the set : @xmath37,\ , y \in [ \alpha(x , z),\beta(x , z)]\right\}\ ] ] where @xmath38 denotes the radius , @xmath39 the physical water height at section @xmath36 .
we denote @xmath40 ( respectively @xmath41 y - coordinate of the left ( respectively right ) boundary point of the domain at altitude @xmath29 , @xmath42 ( see f. [ oyoz ] ) .
we denote also @xmath43 by @xmath44 which is the @xmath29-coordinate of the water level .
[ oxoz ]    in the sequel , we will use a curvilinear map which will be an admissible transformation under the geometrical hypothesis on the domain :    @xmath45 : :    let @xmath46 be the algebraic curvature radius of    the plane curve @xmath47 .
we assume that :    @xmath48      following the work in @xcite , we write system  ( [ e3dicartesian ] ) in the serret - frenet reference frame @xmath32 at point @xmath49 by the transformation @xmath50 using the divergence chain rule lemma that we recall here :    [ lemma ]    let @xmath51 be a @xmath52 diffeomorphism and + @xmath53 the jacobian matrix of the transformation with determinant @xmath54 .
then , for any vector field @xmath55 , one has : @xmath56 and , for any scalar function @xmath57 , one has : @xmath58 where @xmath59 stands for the transpose of the matrix @xmath60 .
let @xmath61 be the components of the velocity vector in the @xmath62 coordinates defined as @xmath63 where @xmath64 is the matrix @xmath65 where we denote by @xmath66 the angle @xmath67 in the @xmath68 plane .
using lemma [ lemma ] , the incompressible euler system in the variables @xmath62 reads : @xmath69 where @xmath70 is the determinant of the transformation and @xmath71 the interested reader can find the details of the calculus in @xcite .
we have denoted by @xmath72 the derivative with respect to the space variable @xmath73 of any function @xmath74 .
on the wet boundary , the no - leak condition reads : @xmath75    [ remarkadmissibilitycondition ]    notice that @xmath76 is the algebraic curvature of the axis at point @xmath77 and the function @xmath78 depends only on the variables @xmath79 .
moreover , under the hypothesis @xmath45 , we have @xmath80 in @xmath81 .
consequently , @xmath82 defines a diffeomorphism and thus the performed transformation is admissible .      in this section ,
we perform a formal asymptotic analysis on system  ( [ e3dicurvilinear ] ) . according to the work in @xcite
, the shallow water equations can be obtained from the incompressible navier - stokes equations with particular boundary conditions . here
, we perform this analysis directly on the incompressible euler system in order to get @xmath83 for some small parameter @xmath84 .
let us introduce the usual small parameter @xmath85 where @xmath86 ( the height ) and @xmath87 ( the length ) are two characteristics dimensions along the @xmath30 and @xmath33 axis respectively .
moreover , we assume that the characteristic dimension along the @xmath88 axis is the same as for the @xmath22 axis .
we introduce the other characteristics dimensions @xmath89 for time , pressure and velocity respectively and the dimensionless quantities as follows : @xmath90 @xmath91 in the sequel , we set @xmath92 and @xmath93 ( i.e. we only consider laminar flows ) .    under these hypotheses , we have @xmath94 . thus , the rescaled system  ( [ e3dicurvilinear ] ) reads : @xmath95 where @xmath96 @xmath97 @xmath98 is the froude number along the @xmath33 axis and the @xmath30 or @xmath34 axis where @xmath31 is any generic variable equal to @xmath87 or @xmath86 .
formally , when @xmath99 vanishes , system ( [ e3dicurvilinearrescaled ] ) reduces to : @xmath100    let us introduce the conservative variables @xmath101 and @xmath102 representing respectively the wet area and the discharge defined as : @xmath103 where @xmath104 is the mean value of the velocity : @xmath105 we integrate the preceding system ( [ e3dicurvilinearrescaledepsilonvanishes1]-[e3dicurvilinearrescaledepsilonvanishes2]-[e3dicurvilinearrescaledepsilonvanishes3 ] ) along the cross - section with the approximation @xmath106 and @xmath107 . then , returning to the physical variables , the free surface model , that we call * fs*-model , reads : @xmath108 where @xmath109 and @xmath110 are respectively the classical term of hydrostatic pressure and the pressure source term defined by : @xmath111 which are obtained from the integration of the pressure term @xmath112 in equation ( [ e3dicurvilinearrescaledepsilonvanishes1 ] ) with @xmath113 ( obtained from equation ( [ e3dicurvilinearrescaledepsilonvanishes3 ] ) ) .    in these formulas @xmath114
is the width of the cross - section at position @xmath73 and at height @xmath29 . the additional term @xmath115 is defined by @xmath116 .
it is the @xmath29-coordinate of the center of mass : @xmath117 in system ( [ fsmodel ] ) , we may add a friction term @xmath118 to take into account the dissipation of energy .
we have chosen this term @xmath119 as the one given by the manning - strickler law ( see e.g. @xcite ) : @xmath120 the term @xmath121 is defined by : @xmath122 , @xmath123 is the strickler coefficient of roughness depending on the material , @xmath124 is the hydraulic radius and @xmath125 is the perimeter of the wet surface area ( length of the part of the channel s section in contact with the water ) .
in this section , we present a new set of unidirectional shallow water like equations to describe pressurized flows in closed non uniform water pipes .
this model is constructed to be coupled in natural way with the * fs*-model  ( [ fsmodel ] ) .
starting from the @xmath0d compressible euler equations in cartesian coordinates , @xmath126 @xmath127 where @xmath5 and @xmath128 denotes the velocity with components @xmath6 and the density respectively .
@xmath129 is the scalar pressure and @xmath10 the exterior strength of gravity .
we define the pressurized domain of the flow as the continuous extension of @xmath81 ( see subsection  [ subsectionframework ] ) defined by some plane curve @xmath19 with parametrization @xmath21 in the cartesian reference frame @xmath2 ; we recall that @xmath23 is then the elevation of the point @xmath35 over the plane @xmath25 ( see f. [ oxoz ] ) .
the curve may be , for instance , the axis spanned by the center of mass of each orthogonal section @xmath130 to the main mean flow axis , particularly in the case of a piecewise cone - shaped pipe .
notice that we consider only the case of infinitely rigid pipes , thus the sections @xmath131 are only @xmath132-dependent .
we then write equations ( [ euler3d_mass_conservation]-[euler3d_momentum_conservation ] ) in the @xmath62 coordinates introduced in subsection  [ subsectionframework ] . as we want a unidirectional model
, we suppose that the mean flow follows the @xmath73-axis . to this end
, we neglect the second and third equation for the conservation of the momentum .
+ by a straightforward computation , the mass and the first momentum conservation equation in the @xmath62 coordinates reads :    @xmath133    applying the same asymptotic analysis developed in subsection [ sectionfsmodelasymptotic ] , equations ( [ euler3d_mass_conservation]-[euler3d_momentum_conservation ] ) read :    @xmath134    we choose the linearized pressure law : @xmath135 ( see e.g. @xcite ) in which @xmath8 represents the density of the fluid at atmospheric pressure @xmath9 , @xmath136 is some function set to zero and @xmath137 the water compressibility coefficient ( equal to @xmath138 in practice ) .
the sonic speed is then given by @xmath139 and thus @xmath140 .
+ for @xmath141 , @xmath142 is the outward unit vector at the point @xmath143 in the @xmath144-plane and @xmath145 stands for the vector * @xmath146 * ( as displayed on f. [ oyoz ] ) .
following the section - averaging method performed in subsection [ sectionfsmodelasymptotic ] , we integrate system  ( [ e3dcrescaledcurvilinear ] ) over the cross - section @xmath144 . noting the averaged values over @xmath144 by the overlined letters ( except @xmath147 ) , and using the approximations @xmath148 the shallow water like equations read : @xmath149 where @xmath150 is the velocity in the @xmath151-plane .
we denote by @xmath152 the area of the cross - section of the pipe at position @xmath73 .
the integral terms appearing in ( [ stv_mass ] ) and ( [ stv_momentum ] ) vanish , as the pipe is infinitely rigid , i.e. @xmath153 ( see @xcite for the dilatable case ) .
it follows the non - penetration condition ( see f. [ restriction ] ) : @xmath154    omitting the overlined letters ( except @xmath147 ) , we introduce the conservative variables @xmath155 and dividing equations ( [ stv_mass])-([stv_momentum ] ) by @xmath8 we get : @xmath156 as introduced previously for the * fs*-model in section ( [ sectionfsmodelasymptotic ] ) , we may introduce the friction term @xmath157 given by the manning - strickler law ( see e.g. @xcite ) : @xmath158 where @xmath159 is defined by : @xmath160 , @xmath123 is the strickler coefficient of roughness depending on the material and @xmath161 is the hydraulic radius where @xmath125 is the perimeter of the wet surface area ( length of the part of the channel s section in contact with the water , equal to @xmath162 in the case of circular pipe ) .
[ [ section ] ]    this choice of variables is motivated by the fact that this system is formally close to the * fs*-model ( [ fsmodel ] ) where the terms @xmath163 , @xmath164 , @xmath115 are respectively the counterparts of @xmath165 , @xmath166 , @xmath167 in system  ( [ pmodel ] ) .
let us remark that the term @xmath147 is continuous through the change of state ( pressurized to free surface or free surface to pressurized state ) when the same curve plane is chosen ( in practice , the main axis of the pipe ) .
then , we are motivated to connect `` continuously '' system ( [ fsmodel ] ) and ( [ pmodel ] ) through transition points ( through the change of state ) by defining a continuous pressure law .
it leads to a `` natural '' coupling between the pressurized and free surface model as we will see in section  [ sectioncoupling ] .
the formulations of the * fs*-model  ( [ fsmodel ] ) and * p*-model  ( [ pmodel ] ) are very close to each other .
the main difference comes from the pressure law . in order to build a coupling between the two models
, we have to define a pressure that ensures its continuity through transition points in the same spirit of @xcite .
as pointed out in the previous section , we will use the common couple of unknowns @xmath168 and the same plane curve @xmath19 ( see remark [ remarkrestrictionplanecurve ] ) to get a continuous model for mixed flows .
[ remarkrestrictionplanecurve ] the plane curve with parametrization @xmath21 is chosen as the main pipe axis in the axisymmetric case .
actually this choice is the more convenient for pressurized flows while the bottom line is adapted to free surface flows .
thus we must assume small variations of the section ( @xmath169 small ) or equivalently small angle @xmath170 as displayed on f. [ restriction ] .
we introduce a state indicator @xmath171 ( see f. [ depression ] ) such that : @xmath172    next , we define the _ physical wet area _ @xmath173 by : @xmath174 and a modified pressure law ( see f. [ depression ] ) which ensures its continuity through the change of state by : @xmath175    *   * indeed , when a change of state occurs we have : @xmath176 which ensures the continuity of the pressure . *
the flux gradient @xmath177 is discontinuous through the change of state since @xmath178    finally , from the * p*-model ( [ pmodel ] ) , the * fs*-model ( [ fsmodel ] ) , the definition of @xmath171 ( [ definitione ] ) , the definition of @xmath173 ( [ definitionphysicalwetarea ] ) and the pressure law ( [ pfsmodelpressureinflux ] ) , the * pfs*-model for unsteady mixed flows can be simply expressed into a single formulation as : @xmath179 where @xmath180 , @xmath181 , and @xmath182 denotes respectively the friction , the pressure source and the geometry source term defined as follows : @xmath183 and @xmath184 stands for @xmath185 .
@xmath186 represents the @xmath29-coordinate of the water level : @xmath187    [ rquerecoveringbothmodels ]    setting @xmath188 in system ( [ pfs ] ) , we obtain obviously the free surface model ( [ fsmodel ] ) . for all pressurized states , when @xmath189 , the pressure law ( [ pfsmodelpressureinflux ] ) reads , for instance , in the case of circular pipe : @xmath190 which is not exactly the pressure law of the * p*-model ( [ pmodel ] ) .
indeed , the derivation of the * p*-model is done with the linearized pressure law ( [ linearizedpressurelaw ] ) ( see section [ sectionpfmodel ] ) with @xmath191 .
thus , the property of the continuity of models ( [ pmodel])-([fsmodel ] ) through a change of state is obtained if and only if @xmath136 is chosen as @xmath192 which is the hydrostatic pressure corresponding to a full section .
[ depression ]   ( top ) , pressurized state with overpressure @xmath193 ( bottom left ) , pressurized state with depression @xmath194 ( bottom right).,title="fig : " ]    the * pfs*-model ( [ pfs ] ) satisfies the following properties :    [ thmpfsmodel ]    1 .    2 .
the right eigenvalues of system ( [ pfs ] ) are given by : @xmath195 with @xmath196 where @xmath197 is the width of the free surface ( see f. [ oyoz ] ) .
+ then , system ( [ pfs ] ) is strictly hyperbolic on the set : @xmath198 3 .   for smooth solutions , the mean velocity @xmath199 satisfies
@xmath200 the quantity @xmath201 is called the total head .
the still water steady state reads : @xmath202 5 .
it admits a mathematical entropy @xmath203 which satisfies the entropy relation for smooth solutions @xmath204    notice that the total head and @xmath205 are defined continuously through the transition points .
the term @xmath206 is also called `` corrective term '' since it allows to write the equations ( [ thmpfsequationforu ] ) and ( [ thmpfsentropy ] ) with ( [ thmpfsmathematicalentropy ] ) .
* proof of theorem  [ thmpfsmodel ] : * the results  ( [ thmpfsequationforu ] ) and ( [ thmpfsentropy ] ) are obtained in a classical way . indeed ,
equation ( [ thmpfsequationforu ] ) is obtained by subtracting the result of the multiplication of the mass equation by @xmath207 to the momentum equation . then multiplying the mass equation by @xmath208 and adding the result of the multiplication of equation ( [ thmpfsequationforu ] ) by @xmath209 , we get : @xmath210 we see that the term @xmath211 is identically @xmath15 since we have @xmath212 when the flow is free surface whereas @xmath213 when the flow is pressurized .
moreover , from the last inequality , when @xmath212 , we have the classical entropy inequality ( see @xcite ) with @xmath205 : @xmath214 while in the pressurized case , it is : @xmath215 finally , the entropy for the * pfs*-model reads : @xmath216 let us remark that the term @xmath217 makes @xmath205 continuous through transition points and it permits also to write the entropy flux under the classical form @xmath218 .
in this section , we adapt the vfroe scheme described in @xcite .
the new terms appearing in the * pfs*-model related to the curvature and the section variation are upwinded in the same spirit of @xcite .
the numerical scheme is adapted to discontinuities of the flux gradient occurring in the treatment of the transitions between free surface and pressurized states .
the spatial domain is a pipe of length @xmath87 .
the main axis of the pipe is divided in cells @xmath220,\
1\leq i\leq n$ ] .
@xmath221 denotes the timestep and we set @xmath222 .
+ the discrete unknowns are @xmath223 . for the sake of simplicity ,
the boundary conditions are not treated ( the interested reader can find this treatment in details in @xcite ) .
we propose to extend the finite volume discretisation @xcite to the * pfs*-model using the upwinding of the new source terms : the curvature and section variation of the pipe . in what follows , we do not write the @xmath171 dependency .    first , following leroux
_ et al . _
@xcite we use piecewise constant functions to approximate @xmath224 @xmath225 as well as the term @xmath226 and the cross section area @xmath152 . adding the equations @xmath227 , @xmath228 and @xmath229 , the * pfs*-model can be written under a non - conservative form with the variable @xmath230 : @xmath231 where @xmath232 and @xmath233 where we have written the pressure source term due to the geometry @xmath234 as @xmath235 .
for instance , for a circular cross - section pipe we have : @xmath236 let @xmath237 be an approximation of the mean value of @xmath238 on the mesh @xmath239 at time @xmath240 . since the values of @xmath241 are known , integrating the above equations over @xmath242x_{i-1/2},x_{i+\frac{1}{2}}[\times [ t_n , t_{n+1 } [ $ ] , we can write a finite volume scheme as follows :    @xmath243    with @xmath244    @xmath245 is the exact or an approximate solution to the riemann problem at interface @xmath246 associated to the left and right states @xmath247 and @xmath248 .
let us also remark that the term @xmath249 does not appear explicitly in this formulation since @xmath184 , @xmath250 and @xmath251 are null on @xmath242x_{i-1/2},x_{i+\frac{1}{2}}[$ ] but contributes to the computation of the numerical flux .
the computation of the interface quantities @xmath252 will depend on two types of interfaces located at the point @xmath253 : the first one is a non transition point , when the flow on both sides of the interface is of the same type .
the second one is a transition point , when the flow changes of type through the interface .
we recall the approach used in @xcite and adapt it here to the new terms . according to the type of interface
, we have to solve two different linearised riemann problems .      expanding the term @xmath254 in the non - conservative equations ( [ pfsncv ] )
, the quasilinear formulation of the * pfs*-model ( [ pfs ] ) reads : @xmath255 with @xmath256 the convection matrix defined by @xmath257    where @xmath258 and @xmath259 denotes the speed of the water .
@xmath260 is then @xmath261 for the pressurized flow or @xmath262 for the free surface flow .
let us remark that , since @xmath263 , the pressure source term @xmath264 does not appear in equation ( [ fnc ] ) .    to compute the interface quantities denoted by @xmath265 for the left hand side and @xmath266 for the right hand side ( see figure [ qamp ] below ) , we solve the following linearised riemann problem : @xmath267 with @xmath268 and @xmath269 where @xmath270 is some approximate state of the left @xmath271 and the right @xmath272 state .
[ remarkchoiceoftilded ]    we will see in section [ sectionwb ] that the classical approximation @xmath273 of the roe matrix @xmath274 defined by @xmath275 is not a suitable choice to preserve the still water steady state .
however , we propose in section [ sectionwb ] a new approximation of @xmath276 which maintains it perfectly .
we have then @xmath277
.    the eigenvalues of the matrix @xmath276 are @xmath278 , @xmath279 , @xmath280 , @xmath281 , @xmath282 and the associated right eigenvectors : @xmath283 @xmath284    ) .
the lines @xmath285 corresponds to the characteristic lines @xmath286 , for @xmath287 . ]
we denote @xmath14 the transition matrix associated to the right eigenvectors of @xmath288 and @xmath289 its inverse . setting @xmath290 ,
the solution of the riemann problem are constant states connected by shocks propagating along the characteristic lines @xmath286 . the jump associated to the eigenvectors
@xmath291 is then equal to @xmath292 . in particular , the discharge is continuous through the line @xmath293 since the fifth component of vectors @xmath294 , @xmath295 and @xmath296 are null .
it follows that the equation associated to the wet area @xmath297 remains conservative .
thus , for instance in the subcritical case ( when @xmath298 ) , we have : @xmath299 where @xmath300 is the upwinded source term @xmath301 .
the friction term can also be upwinded in the same way .
writing the friction term under a conservative form @xmath302 ( for some arbitrary @xmath303 ) allows us to write the `` static '' slope @xmath224 as a `` dynamic '' one as follows : @xmath304 that we denote again @xmath224 .
thus , the upwinding of the dynamic slope @xmath305 reads : @xmath306 which is equal to : @xmath307 since @xmath297 and @xmath209 are constant on each cells .
the terminology `` dynamic '' and `` static '' slope is used since one is @xmath308-dependent while the other is only @xmath132-dependent .      in the case of a transition point , we assume that the propagation of the interface ( pressurized - free surface or free surface - pressurized ) has a constant speed @xmath309 during a time step .
the half line @xmath310 is then the discontinuity line of @xmath311 .
+ let us now consider @xmath312 and @xmath313 the ( unknown ) states respectively on the left and on the right hand side of the line @xmath310 with @xmath314 . both states @xmath315 and @xmath316 ( resp .
@xmath317 and @xmath318 ) correspond to the same type of flow .
thus it makes sense to define the averaged matrices in each zone as follows :    * for @xmath319 , we set @xmath320 for some approximation @xmath321 which connects the state @xmath271 and @xmath322 ( see remark [ remarkchoiceoftilded ] ) . * for @xmath323 , we set @xmath324 for some approximation @xmath321 which connects the state @xmath325 and @xmath272 ( see remark [ remarkchoiceoftilded ] ) .
then we formally solve two riemann problems and use the rankine - hugoniot jump conditions through the line @xmath310 which writes : @xmath326 with @xmath327 . according to the left ( @xmath316 , @xmath328 ) and right unknowns ( @xmath318 , @xmath329 ) at the interface @xmath330 and the sign of the speed @xmath309
, we have to deal with four cases :    * pressure state propagating downstream , * pressure state propagating upstream , * free surface state propagating downstream , * free surface state propagating upstream .
we can next consider two couples of `` twin cases '' : pressure state is propagating downstream ( or upstream ) as shown in the figure [ plus1 ] and free surface state propagating downstream ( or upstream ) as shown in the figure [ plus2 ] .
moreover , for all existing transition case , the upwinded altitude term @xmath331 in @xcite are replaced by @xmath332 .    [ [ pressure - state - propagating - downstream ] ] pressure state propagating downstream : + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    it is the case when on the left hand side of the line @xmath333 , we have a pressurized flow and on the right hand side we have a free surface flow : the speed @xmath309 of the transition point being positive . following song @xcite ( see also @xcite ) , an equivalent stationary hydraulic jump must occur from a supercritical to a subcritical condition and thus the characteristics speed satisfies the inequalities : @xmath334 where @xmath261 is the sound speed for the pressure flow , @xmath335 , @xmath336 , and @xmath337 are defined by the same formula obtained in the case of a non transition point but according to @xmath338 and @xmath339 .        therefore , only the characteristic lines drawn with solid lines are taken into account
. indeed they are related to incoming waves with respect to the corresponding space - time area @xmath340 .
conversely , the dotted line @xmath341 , for instance , related to the free surface zone but drawn in the area of pressurized flow is a `` ghost wave '' and is not considered .
thus @xmath342 and @xmath315 , @xmath316 are connected through the jumps across the characteristics @xmath343 and @xmath344 .
eliminating @xmath309 in the rankine - hugoniot jump relations ( [ rh1])-([rh2 ] ) , we get @xmath316 as the solution to the nonlinear system :    @xmath345    finally , we obtain : @xmath346    [ [ free - surface - state - propagating - downstream ] ] free surface state propagating downstream : + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    on the left hand side of the line @xmath333 we have a free surface flow while on the right hand side , we have a pressurized flow ( the speed @xmath309 of the transition point being positive ) .
following song @xcite again , the characteristic speed satisfies the inequalities : @xmath347        there are two incoming characteristic lines with respect to the free surface area @xmath340 ( actually three with @xmath343 ) and they can connect the given left state @xmath315 with any arbitrary free surface state @xmath328 . thus only one characteristic line ( @xmath348 ) gives any information ( it is the equation ( [ slchp1 ] ) above ) as an incoming characteristic line with respect to the pressurized zone @xmath349 . from the jump relations through the characteristic @xmath350 , and after the elimination of @xmath309 in the rankine - hugoniot jump relations ( [ rh1]),([rh2 ] ) we get another equation , namely equation ( [ slchp2 ] ) above .
it remains to close the system of four unknowns @xmath351 .
firstly , we use a jump relation across the transition point ( with speed @xmath309 ) for the total head @xmath352 , from equation ( [ thmpfsequationforu ] ) , which writes : @xmath353 finally , we use the relation : @xmath354 we have then to solve the nonlinear system : @xmath355    the states @xmath328 et @xmath329 are then obtained by the following identities : @xmath356 \dsp ap =   am + \frac{g\,\at_l\,\psi_l^r}{\ut_l^2-\ct_l^2}\\[5 mm ] \dsp qm = qp = qmp = q_l + \frac{g\,\at_l\,\psi_l^r}{2\,\ct_l } \;+ \\
\hspace{4.8cm}\dsp + \;\frac{\ut_l^2-\ct_l^2}{2\,\ct_l}\,(a^- - a_l ) - \frac{\ut_l - \ct_l}{2\,\ct_l}\,(q^- - q_l ) \end{array}\ ] ] finally , the update state @xmath357 are obtained by the same relation as in the case of a non transition point .
using equations ( [ fvscheme ] ) we update the values of @xmath358 with a standard stability condition of courant - friedrich - levy controlling the time step size @xmath221 .
to update the state @xmath171 in the cell @xmath239 ( see f. [ updatea ] ) , we use a discrete version of the state indicator @xmath171 defined by ( [ definitione ] ) equal to @xmath360 for a pressurized flow and @xmath15 otherwise . following @xcite , after the computation of the wet area @xmath361 we predict the state of the flow in the cell @xmath239 by the following criterion :    * if @xmath362 then : + if @xmath363 then @xmath364 , else @xmath365 , * if @xmath366 : + if @xmath367 then @xmath365 , else @xmath368 .
indeed , if @xmath367 it is clear that the state of the flow in the cell @xmath239 becomes pressurized , on the other hand if @xmath369 in a mesh previously pressurized , we do not know _ a priori _ if the new state is free surface ( @xmath370 and the value of the wetted area is less than @xmath371 ) or pressurized ( in depression , with @xmath372 and the value of the wetted area is equal to @xmath373 : see remark [ remarkdepressionoverpressurestates ] and f. [ rhocriterion ] ) .
+ so far as we do not take into account complex phenomena such that entrapment of air pockets or cavitation and keeping in mind that the cfl condition ensures that a transition point crosses at most one mesh at each time step , we postulate that :    1 .
if the state of the flow in the cell @xmath239 is free surface at time @xmath240 , its state at time @xmath374 is only determined by the value of @xmath361 and it can not become in depression .
2 .   if the state of the flow in the cell @xmath239 is pressurized at time @xmath240 and if @xmath363 , it becomes free surface if and only if at least one adjacent cell was free surface at time @xmath240 .
this is exactly the discrete version of the continuous @xmath375 criterion in remark [ remarkdepressionoverpressurestates ] and displayed on f. [ rhocriterion ] .    [
cols="^,^ " , ]     [ remarkdepressionoverpressurestates ]    as we do not take into account complex phenomena such that entrapment of air pockets , each connected component of the pressurized area is simply connected ( see f. [ rhocriterion ] ) .
moreover , for each depression area @xmath256 , its closure @xmath376 is a strict subset of the pressurized set .
it follows that when @xmath377 on each pressurized area , we observe a depression as displayed on f. [ rhocriterion ] .
moreover , we may also use a visual depression indicator given by the function @xmath375 : the case @xmath212 corresponds to a free surface state while @xmath378 to an overpressure state and @xmath379 to a depression state . on f.
[ rhocriterion ] , we draw the behavior of the interface speed @xmath309 in the @xmath380-plane and the graph of the function @xmath375 at fixed time @xmath381 .
this section is devoted to the construction of an exactly well - balanced scheme in the sense that it maintains perfectly the still water steady states .
this scheme , noted ewbs , is obtained by a suitable definition of the convection matrix .
the numerical approximation of the * pfs*-model ( [ pfsmodelpressureinflux ] ) reads : @xmath382 where @xmath383 stands for @xmath384 and @xmath385 .
for instance , in the subcritical case , the interface quantities reads : @xmath386 where the upwinded source term reads : @xmath387 and @xmath388 stands for @xmath389 with @xmath390 given by @xmath391 and the approximation of @xmath392 to be specified .    starting from a discrete state @xmath393 at time @xmath240 such that :    ( h ) : :    let @xmath394 such that : @xmath395 ,    @xmath396 and @xmath397 satisfy the discrete    still water steady state equation ( according to equation    ( [ thmpfsequationforu ] ) ) :    @xmath398    we will say that :    [ definitionwbs ]    1 .    2 .   the numerical scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) for some approximations of the terms @xmath399 is @xmath400 well - balanced ( also denoted by @xmath400-wb ) if : @xmath401 with @xmath402 , @xmath403 the well - balanced order of the numerical scheme ( [ numericalschemea ] ) and ( [ numericalschemeq ] ) respectively .
3 .   the numerical scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) for some approximations of the terms @xmath399 is exactly well - balanced ( also denoted by ewb ) if : @xmath404    we will denote by @xmath400-wbs the @xmath400 well - balanced scheme and ewbs the exactly well - balanced scheme .    ( sf ) : :    in the rest of this paper , we assume @xmath405    constant , the radius @xmath406 and @xmath224 are given    affine functions , the angle @xmath407 is constant which    implies that the jumps across the interface @xmath246 :    @xmath408 ,    @xmath409 are constant and    @xmath410 is null .    then under this simplified framework , we show that :    [ thmwb ]    1 .    2 .
* the numerical scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) with the classical choice @xmath411 and non constant section @xmath152 is not well - balanced in the sense of definition [ definitionwbs ] ( we have @xmath412 ) .
* for constant section and @xmath413 , the numerical scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) with ( [ demisomme ] ) is ewb .
3 .   under a suitable choice of @xmath414 ,
the numerical scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) is ewb .    the following section deals with the numerical scheme with @xmath392 defined as ( [ demisomme ] ) where we show the first point of theorem [ thmwb ] . the second point is studied in section [ subsectionewbs ] where a convenient definition of @xmath392 leads to an ewbs
the simpler choice of definition for the convection matrix @xmath273 is the one obtained by the approximation of the mean value of the roe matrix @xmath415 .
this approximation is given by @xmath416 that we call `` the classical approximation '' .
thus , defining @xmath392 as follows : @xmath417 provides the classical approximation .
but , it is not suitable to preserve the still water steady state : we will see that the numerical scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) with ( [ demisomme ] ) defines a non well - balanced scheme in the sense of definition ( [ definitionwbs ] ) since @xmath412 .    to this end , let us assume * ( h ) * and * ( sf ) * at time @xmath240 , equations ( [ interfacequantitiessubcriticalcase ] ) read for every @xmath418 : @xmath419    denoting @xmath420 by @xmath421 , we have : @xmath422 where @xmath423 denoting then @xmath424 and observing that @xmath425 and @xmath426 we deduce : @xmath427 it follows that : @xmath428    then , we denote @xmath429 by @xmath430 with    @xmath431 and @xmath432 where @xmath433 and @xmath434 as the term @xmath435 is at least of order @xmath436 since @xmath437 we have : @xmath438 it follows that : @xmath439    for constant section @xmath152 and @xmath413 , it is easy to see that the numerical scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) with ( [ demisomme ] ) is ewb .
although the scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) with ( [ demisomme ] ) have an order @xmath412 for non constant section , the still water steady state for the pressurized case is very well maintained for great value of the sonic speed @xmath261 .
but it is not the case for the free surface numerical scheme ( see f. [ statiosl ] with @xmath440 , @xmath441 , uniform pipe with diameter @xmath360 ) .
we plot in f. [ comparaisonmeanmatrixstationnarymatrixc30 ] and f. [ comparaisonmeanmatrixstationnarymatrixc100 ] a still water steady state for two values of @xmath261 for a given @xmath442 ( with @xmath443 ): the absolute error obtained is @xmath444 for @xmath445 and @xmath446 for @xmath447 .
indeed , let us consider , for the sake of simplicity , the case @xmath448 .
at the continuous level , the still water steady state equation reads : @xmath449 with the hypothesis * ( h ) * , in particular using equation ( [ discretesteadystateequation ] ) , we write : @xmath450 given @xmath451 , for great value of @xmath261 , we can approximate @xmath452 by : @xmath453 then , replacing the right hand - side of @xmath452 in ( [ diffainterface ] ) , we have @xmath454 as we can see @xmath455 is always equal to @xmath360 but the constant @xmath456 plays the role of a smoothing term which helps the scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) with ( [ demisomme ] ) to stabilize rapidly towards the equilibrium .
physically , @xmath261 is approximatively @xmath457 ( for a pressurized flow without air ) , thus @xmath458 . since @xmath459
, this feature is not observed for the free surface numerical scheme .
this section is devoted to the construction of an ewbs .
we have seen in subsection [ subsectionsteadystateswithclassicalapproximation ] that the classical approximation of with @xmath414 ( [ demisomme ] ) is not appropriate to preserve the still water steady state .
thus , we have to find a suitable definition for @xmath414 to obtain an exactly well - balanced scheme . for this purpose ,
let us assume * ( sf ) * and start with :    at the discrete level , the still water steady state is perfectly maintained ( see f. [ discretrepresentationdomain ] ) : it exists @xmath394 such that for every @xmath418 , if @xmath460 and @xmath461    a1 : : :    @xmath462 , a2 : : :    @xmath463 , a3 : : :    @xmath464 ,    then , for all @xmath465 the conditions a1 , a2 and a3 holds .    . ]
the condition a2 is satisfied if and only if @xmath466    the condition a3 is satisfied if and only if @xmath467    the condition a1 is satisfied for pressurized flows if and only if @xmath468    the condition a1 is satisfied for free surface flows if and only if @xmath469 for circular cross - section pipe , @xmath452 is computed by : @xmath470 with @xmath471 is the angle displayed on f. [ computationomega ] .    ]
thus , the discrete still water steady state is perfectly maintained if and only if , @xmath472 is the solution of the non - linear system : @xmath473 where we have replaced the expression of @xmath452 in ( [ conditiona2]-[conditiona3 ] ) by ( [ conditiona1pf ] ) for pressurized by ( [ conditiona1fs ] ) for free surface flows : @xmath474 with @xmath475 .
for circular pipe , @xmath476 is given by ( [ forcircularpipe ] ) .
finally , the numerical scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) with @xmath477 as the solution of the non linear system ( [ snl ] ) defines an exactly well - balanced scheme .    for uniform pipe and pressurized flow ,
the previous system simply writes : @xmath478 the solution is easily obtained by : @xmath479 let us also remark that using the relation @xmath480 , we have : @xmath481 it follows that @xmath392 can be expressed as the mean value of ( [ atildebyai ] ) and ( [ atildebyai+1 ] ) as follows : @xmath482 for small @xmath451 , we have @xmath483 it follows that the scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) with ( [ demisomme ] ) is the zero order approximation of the solution given by the ewbs .
the same analysis shows that the free surface numerical scheme with ( [ demisomme ] ) is also the zero order approximation of the solution given by the ewbs .
[ [ section-1 ] ]    on f. [ comparaisonmeanmatrixstationnarymatrixc30 ] , f. [ comparaisonmeanmatrixstationnarymatrixc100 ] , f. [ statiosl ] , we display the still water steady pressurized and free surface state computed by the ewbs and the scheme with the approximation ( [ demisomme ] ) .
the steady state with the ewbs is numerically well preserved while as pointed out before ( see subsection [ subsectionsteadystateswithclassicalapproximation ] ) the classical approach is not convenient .
we also display an unsteady simulation on f. [ testnonstationnaire ] where the results of the two methods are very well reproduced .
the previous sections deals with the well - balanced property of the numerical scheme ( [ numericalschemea]-[numericalschemeq])-([definitionoftildebthetasq ] ) for free surface and pressurized flows . to use the well - balanced scheme developed in subsections ( [ subsectionsteadystateswithclassicalapproximation ] ) and ( [ subsectionewbs ] )
, we start from the discrete representation of a mixed still water steady state ( as displayed on f. [ discretrepresentationdomain ] ) .
assume that there exists @xmath484 such that all cells @xmath239 on the left hand side of the interface @xmath485 are free surface while the other are pressurized : @xmath486    the interface @xmath485 is such that the speed of propagation of the interface is null : @xmath487 ( see f. [ wzero ] ) .
therefore , @xmath488 and @xmath489 .
for example on f. [ wzero ] , we apply the free surface numerical scheme on the left hand side of the interface and the pressurized one on the right cells .
as the ewbs preserves the pressurized and free surface still water steady state , it also preserves the mixed still water steady state .
the numerical validation for pipes with constant section and slope has been previously studied by two of the authors in @xcite and thus are not presented in this paper . since experimental data for mixed flows in any pipe are not available , we focus on the behavior of our method for several circular cross - section contracting and expanding pipe .
notice that , the equivalent pipe method is not relevant for the mixed flows as pointed out by @xcite for instance .
the mixed flow case is numerically performed on a water hammer test . starting from an horizontal free surface still water steady state ,
the water hammer occurs immediately after the increase of the upstream piezometric head while the downstream discharge is set to @xmath15 .
the prescribed hydrograph produces a travelling wave which produces a pressurized state propagating from upstream to downstream end .
physically an trapped air pocket may appear : it is not taken into account in the * pfs*-model .
actually , the trapped air pockets vanish or move ; some parts of these pockets undergo condensation / vaporisation and others parts move and lead to a two phase flow .
consequently the sound speed decreases .
as our model does not take into account these phenomena , the value of @xmath261 is assumed to be constant .
moreover we should have to deal with the entrapment of air bubbles which have a non negligible effect ( see @xcite for instance ) .
the numerical experiments are performed in the case of a @xmath490 @xmath143 long closed circular pipe at altitude @xmath491 with @xmath15 slope which corresponds to the elevation and slope of the main pipe axis ( we have @xmath492 ) .
the manning roughness coefficient is @xmath493 .
the simulation starts from a steady state as a free surface flow with a discharge @xmath494 .
the upstream boundary condition is a prescribed hydrograph ( see f. [ boundarycondition ] ) while the downstream discharge is kept constant to @xmath495 ( as displayed on f. [ boundarycondition ] ) .
we compare then the results obtained for uniform , contracting an expanding pipes . for each test , the parameters are the same except the downstream diameter : the upstream diameter is kept constant to @xmath496 .
the contracting pipe is chosen for @xmath497 and the expanding one for @xmath498 ( where @xmath256 denotes the downstream diameter ) .
let us recall that the zero water level corresponds to the main pipe axis .
the piezometric head is defined as @xmath499 : @xmath500    results are then represented on f. [ x50 ] .
the sudden raise of the upstream piezometric level produces a pressurized state with a travelling wave . a water hammer
is then observed since the downstream discharge is null . a careful analysis of the flow which is performed by the variable @xmath171 or equivalently by @xmath375 ( see remark [ remarkdepressionoverpressurestates ] and f. [ rhocriterion ] ) shows that after this transition point , the flow is pressurized but in depression which starts approximatively at time @xmath501 for the contracting pipe , @xmath502 for the uniform pipe and @xmath503 for the expanding one .
we display the piezometric line for different times around the depression time ( see f. [ depress ] ) and also the graph of the function @xmath375 which confirm that the observed times correspond exactly to a depression state for each pipes ( see f. [ depressx ] ) .
we also observe a little smoothing effect and absorption due to the first order discretisation type .
we have derived a free surface and a pressurized model which have been coupled using a common set of variables and a suitable pressure law .
we have thus obtained a mathematical model for unsteady mixed flows in non uniform water pipes , that we have called * pfs*-model .
this model takes into account the local perturbation of the section and of the slope .
the * pfs * model is numerically solved by a vfroe scheme using the interfacial upwind to include the source terms into the numerical fluxes .
we have shown that the classical approximation of the convection matrix ( the roe matrix approximation ) is not suitable to preserve the still water steady state ( except for the pressurized case where the value of @xmath261 helps the scheme to maintain this state ) . moreover ,
we have proposed a manner to obtain an exactly well - balanced scheme .    as mentioned in @xcite this numerical method with the classical approximation of the convection matrix , for constant section ,
reproduces correctly laboratory tests for uniform pipes and can deal with multiple points of transition between the two types of flows . as pointed out before , due to the lack of experimental data for non uniform pipes , we have only shown the behavior of the piezometric line which seems reasonable ( at less no major difference was observed ) .    as a well - known feature on approximate godunov scheme ,
the upwinding of the source terms may introduce stationary waves with a vanishing denominator when critical flows occurs i.e. @xmath504 .
moreover , in its actual form , the presented numerical scheme is not able to deal with drying and flooding area .
nevertheless , it may be possible to introduce a cut - off function to avoid division by zero for each problems : critical stationary waves , drying area and flooding area .
but it is not the better choice that we can propose , since , truncation of the wet area induces a loss of water mass leading to the non - conservativity of the mass .
nevertheless , at the present time , we are interested in a mathematical kinetic formulation of the * pfs * model and the construction of a numerical kinetic scheme that avoids all these drawbacks @xcite .
the next step is to take into account the air entrainment which may have non negligible effects on the behaviour of the piezometric head .
a first approach has been derived in the case of perfect fluid and perfect gas modelised as a bilayer model based on the * pfs*-model @xcite .
nguyen  trieu dong .
sur une mthode numrique de calcul des coulements non permanents soit  surface libre , soit en charge , soit partiellement  surface libre et partiellement en charge . , 2:149158 , 1990 .
f.  kerger , p.  archambeau , s.  erpicum , b.  j. dewals , and m.  pirotton .
improved one - dimensional numerical simulation of transient mixed flow in water pipe . in _ proceedings of 4th int .
conf . on advanced computational methods in engineering liege , belgium . _ , 2008 .
f.  kerger , s.  detrembleur , p.  archambeau , s.  erpicum , b.  j. dewals , and m.  pirotton .
an experimental analysis of effects induced by moving bodies in shallow water . in _ proceedings of 2nd int .
junior researcher and engineer workshop on hydraulic structures , pisa , italy .
_ , 2008 .
p.l . roe . some contributions
to the modelling of discontinuous flow . in b.
e. engquist , s.  osher , and r.  c.  j. somerville , editors , _ large - scale computations in fluid mechanics .
proceedings of the fifteenth ams - siam summer seminar on applied mathematics held at scripps institution of oceanography , la jolla , calif . , june 27-july 8 , 1983 _ , volume  22 of _ lectures in applied mathematics _ , pages 163193 .
american mathematical society , 1985 .


in this paper we demonstrate that the profile of angular distribution of electrons emitted in the process of a single - photon ionization from a strongly polarizable target ( in particular , a metal cluster anion ) exposed to a laser field can be noticeably modified due to the action of the field @xmath2 of alternating dipole moment @xmath3 induced at the residue by the laser field .
the field @xmath2 acts on the escaping electron and , being dependent on the polarization vector of the laser field , brings additional dependence of the cross section on the emission angle . the degree of this dependence
is determined by the magnitude of @xmath3 which , in turn , is proportional to the residue s dynamic polarizability .
therefore , for a target with large polarizability the effect can be very pronounced even for laser fields of low intensities .
the effect of the laser - field  induced dipole moment in atomic photoionization process was discussed in @xcite ( where one finds also the references to the earlier works by the authors ) and , more recently , in @xcite .
both papers were devoted to the study of _ multiphoton _ detachment of electrons from atomic negative ions . in @xcite
the process was analyzed by means of the perturbation theory .
numerical analysis was carried out for halogen and alkali metal anions , and it demonstrated , that the induced dipole moment influences the magnitude of total cross section of the many - photon detachment . the paper @xcite dealt with the process of a _ strong - field _ many - photon detachment of an atomic anion .
the non - perturbative adiabatic theory , developed in @xcite , was modified to account for the field of the induced dipole moment .
the numerical calculations , performed for rb@xmath4 and cs@xmath4 , showed the increase of the detachment rate and essential modification of the angular distribution .    in our paper
we investigate the influence of the induced dipole moment on the angular and spectral distributions of photoelectrons in the process of _ single - photon _ ionization ( or detachment ) occurring in the _ weak - field _ regime .
numerical results are presented for sodium cluster anions na@xmath5 with @xmath6 .
the corresponding residue clusters are characterized by large values of the dynamic dipole polarizability @xmath7 , which in the vicinity of plasmon resonances , i.e. at @xmath8 ev , can be as large as @xmath9 a.u . , which by far exceed the maximum values of @xmath7 of individual atoms ( e.g. , @xmath10 a.u . for cs at the photon energy 0.47 ev @xcite ) .
we demonstrate that dramatic changes to the profile of angular distribution occur in vicinities of the plasmon resonances and for comparatively low intensities @xmath11 w/@xmath1 of the laser field .
these are exactly the ranges used in recent experimental studies @xcite of the photodetachment from various metal cluster anions carried out by means of angle resolved photoelectron spectroscopy .
therefore , the experimental test of the predicted effects seems to be feasible .    the atomic system of units , @xmath12 is used throughout the paper .
consider the process of photoionization of a strongly polarizable target ( a cluster , a fullerene , a complex molecule , etc ) by spatially homogenous linearly polarized laser field @xmath13 we assume , that the field wavelength @xmath14 greatly exceeds the characteristic size @xmath15 of the target , @xmath16 .
thus , the process can be treated within the dipole - photon limit .
the laser field ionizes the target and polarizes the residue .
therefore , in addition to other fields ( the static potential of the residue , the long - range polarizational potential ) the escaping photoelectron feels the field due to the oscillating dipole moment @xmath3 induced at the residue : @xmath17    to simplify the consideration we assume , that the residue is spherically symmetric , so that its tensor of polarizability contains only a scalar part which is notated as @xmath7 .
then , the moment @xmath3 reads @xmath18 where @xmath19 .
the field of the induced dipole can be written as follows : @xmath20 with @xmath21 .
this field is dependent on the direction of @xmath22 , therefore , the angular distribution of photoelectrons acquires additional dependence on the polarization of the laser field .
the degree of this dependence is determined by the magnitude of @xmath3 , which is proportional to the residue dynamic dipole polarizability .
hence , for a target with large polarizability the modification of the angular distribution due to the induced dipole moment can become noticeable even for a laser field of a comparatively low intensity .
let us clearly state the conditions which we impose on the parameters of the laser field .
these include the following inequalities : @xmath23 the first condition implies that the laser field strength is much smaller than the strength @xmath24 of the internal field in the target .
this condition allows one to ignore the laser field when constructing the wavefunction of the bound electron with the binding energy @xmath25 .
numerical results presented in section [ numericalresults ] refer to the laser field strength @xmath26 v / cm .
the second condition states that the energy @xmath27 of the escaping electron is large compared to the ponderomotive shift due to the wiggling in the laser field .
owing to this condition we do not `` dress '' the photoelectron wavefunction with the laser field .
the last inequality in ( [ conditions ] ) indicates that the photon energy exceeds the ionization threshold of the bound electron .
therefore , the process of a _ single - photon _ ionization is allowed for which the energy conservation law reads @xmath28    combining the conditions ( i ) and ( iii ) one demonstrates that the following inequality is met : @xmath29 where @xmath30 is a so - called keldysh parameter , and @xmath31 .
the limit @xmath32 defines a weak field regime ( see , e.g. @xcite ) .
thus , in this paper we investigate the influence of the induced dipole moment @xmath3 on the angular and spectral distributions of photoelectrons in the process of _ single - photon _ ionization ( or detachment ) occurring in the weak field regime .
earlier the role of @xmath3 was studied in the process of many - photon detachment ( i.e. , when @xmath33 ) from atomic negative ions .
it was done in the weak - field @xcite and in the strong - field @xcite limits .
let us note , that the weak field condition ( ii ) does not immediately imply the applicability of the perturbative approach with respect to @xmath2 .
indeed , from ( [ problem.2])-([problem2.3a ] ) follows that @xmath34 exceeds @xmath35 by a factor @xmath36 . for a strongly polarizable target
this factor can be large enough in a wide range of radial distances @xmath37 , and this might lead to a non - perturbative treatment of the action of @xmath2 .    to carry out a quantitative analysis
we first postulate that under certain conditions the field ( [ problem2.3a ] ) can be considered as a time - dependent perturbation which modifies the wavefunction of the escaping electron .
the criterion of applicability of this approach will be formulated in the course of calculations .
the matrix element @xmath38 , which describes a dipole - photon transition of the electron from the initial bound state @xmath39 to the final state @xmath40 with the asymptotic momentum @xmath41 and energy @xmath42 , can be written in the following form : @xmath43 the right - hand side of this equation represents the power series in @xmath2 .
the structure of the terms of the series is illustrated by diagrams in figure [ diagram.fig ] .
( indicated with the dashed vertical lines ) .
the double line represents the initial bound state .
the solid line stands for the escaping electron whose wavefunction is `` dressed '' with all static fields . ]
the zeroth order term @xmath44 describes the direct process of photoionization standing for the matrix element of the operator @xmath45 : @xmath46 where @xmath47 we assume that the wavefunction @xmath48 ( whose asymptotic form is the plane wave + the incoming wave ) corresponds to the state `` dressed '' with all _ static _ fields .
the delta - function on the right - hand side of ( [ m0.1 ] ) expresses the energy conservation law ( [ m0.3 ] ) .
the terms @xmath49 with @xmath50 are the corrections to the transition amplitude due to the @xmath51-times interaction of the photoelectron with @xmath2 . since this field explicitly depends on @xmath52 ( see ( [ problem2.3a ] ) ) then in each act of the interaction the electron energy is changed ( increased or decreased ) by @xmath53 . as a result ,
not all terms on the right - hand side of ( [ problem.3 ] ) represent the correction to the direct amplitude @xmath44 , which implies the validity of ( [ m0.3 ] ) .
for example , evaluating @xmath54 ( the second diagram in figure [ diagram.fig ] ) one finds , that it contains the terms proportional to @xmath55 and to @xmath56 , which conflict with the conservation law ( [ m0.3 ] ) .
therefore , the term @xmath54 must be ignored when constructing the correction to @xmath44 .
the lowest - order correction originates from the amplitude @xmath57 ( the last diagram in figure [ diagram.fig ] ) , which contains the terms proportional to @xmath58 .
applying the standard technique of the time - dependent perturbation theory , one derives the following expression for this contribution : @xmath59 where @xmath60 here @xmath61 is the green s function of the photoelectron and @xmath62 is the unit vector of the laser field polarization .
hence , within the second second order perturbation theory in @xmath2 , the photoionization amplitude is written as : @xmath63    using the amplitude ( [ ampl.1 ] ) one can write the following expression for the spectra - angular distribution : @xmath64 here @xmath65 stands for the differential cross section calculated without the influence of the induced dipole moment . within the dipole - photon approximation
the differential cross section @xmath66 can be presented in the well - known general form ( see , e.g. , @xcite ) : @xmath67 where @xmath68 is the total cross section of photoionization , @xmath69 stands for the angle between @xmath41 and @xmath70 , @xmath71 is the legendre polynomial of the second order .
the quantity @xmath72 stands for the angular asymmetry parameter @xcite .
the factor @xmath73 defines the profile of the angular distribution .
equation ( [ cs.3 ] ) indicates , that the profile is modified by the factor @xmath74 , when the field ( [ problem2.3a ] ) is taken into account .
the correction term @xmath75 is defined as follows : @xmath76    formula ( [ cs.3 ] ) was obtained by means of the perturbation theory .
thus , it is implied that the absolute value of the correction term satisfies the condition @xmath77 .
if otherwise , the field @xmath2 must be accounted for in all orders of the perturbation series ( [ problem.3 ] ) .
it follows from ( [ ampl.3 ] ) , that @xmath75 is proportional to the squared induced dipole moment , @xmath78 .
therefore , for a target of a small polarizability , one has to apply very intensive laser fields in order to achieve a noticeable correction to the photoelectron angular distribution , i.e. @xmath79 .
on the contrary , for a strongly polarizable target ( e.g. , a metallic cluster or / and its anion ) a strong effect can be expected for laser fields of a moderate strength . in section [ numericalresults ]
we demonstrate , that dramatic change of the angular distribution profile can occur for @xmath80 v / cm .
it is important to note , that the correction term ( [ cs.5 ] ) depends on the ratio of the amplitudes @xmath81 and @xmath82 .
thus , one can expect that it is less sensitive to the approximation chosen to describe the photoelectron wavefunction than the cross section @xmath83 .
therefore , to estimate the influence of the target polarization one can calculate @xmath75 within the framework of a comparatively simple approximation , which is described below .
in this section we derive an approximate expression for the correction term @xmath75 using a simple , although realistic and physically clear , model .    to start with
let us assume , that the ionized target is a spherically - symmetric anion . then , firstly , the ground state wavefunction @xmath84 does not depend on the angles of the position vector : @xmath85 .
secondly , the wavefunction @xmath86 of the photoelectron does not contain the contribution of a long - range coulomb field of the residue .
these assumptions , being not of a principal nature , allow one to simplify the intermediate algebra .
making use of the partial wave expansion @xmath87 where @xmath88 are the spherical harmonics , @xmath89 are the scattering phaseshift , and @xmath90 are the radial wavefunctions , one writes the amplitude @xmath82 in the following form : @xmath91    another simplification is based on the approximation that the field @xmath2 acts on the escaping electron only outside the target , i.e. at the distances @xmath92 .
additionally , one can assume , that the wavefunction @xmath93 is concentrated in the interior of the target , i.e. at @xmath94 .
therefore , when evaluating the right - hand side of ( [ ampl.3 ] ) one takes into account that @xmath95 .
this allows one to relate @xmath81 to the amplitude @xmath82 . to do this
let us approximate the green s function @xmath96 as follows : @xmath97 here @xmath98 is the largest / smallest of @xmath99 , @xmath100 stands for the irregular solutions of the radial schrdinger equation for the energy @xmath42 .
the quantity @xmath101 is the exact radial green s function corresponding to the multipolarity @xmath102 .
the first equality in ( [ longdistance.2 ] ) states the partial - wave expansion of the green s function .
the approximate equality is based on the assumption @xmath103 .    using ( [ longdistance.2 ] ) in ( [ ampl.3 ] ) one notices that the radial integral over @xmath104 reduces to that from ( [ m0_born.1 ] ) .
therefore , the amplitude @xmath81 acquires the form : @xmath105 with @xmath106 here the integration is carried out over the spatial region @xmath107 .
the correction term ( [ cs.5 ] ) is related to @xmath108 as @xmath109 .
let us comment on the evaluation of the right - hand side of ( [ longdistance.8 ] ) . introducing the expansion ( [ m0.1a ] ) and representing the green s functions @xmath110 and @xmath111 in terms of the partial - wave series similar to ( [ longdistance.2 ] ) , one carries out the angular integration by means of the standard methods ( see , e.g. , @xcite ) and arrives at : @xmath112 { { \rm im}}\ , b_{ll^{\prime}}. \label{longdistance1.5}\end{aligned}\ ] ] the kronecker symbols @xmath113 , which reduce the allowed values of the orbital momenta to @xmath114 and @xmath115 , reflect the dipole selection rules applicable to each of the vertices in the last diagram in figure [ diagram.fig ] .
the ratio @xmath116 reduces to @xmath117 with the help of the recurrence relation for the legendre polynomials .
the quantities @xmath118 stand for the following radial integrals : @xmath119 for any realistic target the exact evaluation of these integrals can be done by numerical means only , implying the computation of the regular @xmath90 and the irregular @xmath120 functions as well as of the green s functions @xmath121 and @xmath122 .
to carry out the approximate integration one can substitute these quantities with the leading terms of their asymptotic representations . for @xmath90 and @xmath120
this results in @xmath123 and @xmath124 .
the leading terms in the expansions of the radial green s functions are : @xmath125 here @xmath126 , @xmath127 , @xmath128 and @xmath129 are the elements of the scattering matrix . for the real argument
@xmath130 this quantity reads as @xmath131 .    by using the asymptotic formulae one carries out analytical evaluation of the integrals from ( [ longdistance1.5a ] ) . the final result ( the leading term ) for @xmath132 is as follows : @xmath133 _ { l=1,3 \atop l^{\prime}=0,2 } , \label{functionb3.3}\end{aligned}\ ] ] where @xmath134 , @xmath135 and @xmath136 .
it follows from the approximations formulated above , that formula ( [ functionb3.3 ] ) is valid if the photoelectron momentum @xmath137 and the ground state energy @xmath138 satisfy the conditions : @xmath139    letting @xmath140 in ( [ functionb3.3 ] ) one finds the plane - wave born limit of @xmath132 .
the presence of oscillatory terms on the right - hand side of ( [ functionb3.3 ] ) is physically clear . indeed , in the limit @xmath141 the wavelength of the photoelectron is less than the size of the target , which is supposed to have a well - pronounced edge ( e.g. , the edge of an ionic core in a cluster ) .
thus , the electron experiences a diffraction at the edge when leaving the core .
this diffraction leads to the oscillatory character of the cross section .
this phenomenon was discussed previously in connection with various processes involving metallic clusters and fullerenes : the radiative electron capture by metallic clusters @xcite , the electron scattering @xcite , the polarizational bremsstrahlung @xcite .
carrying out obvious summations in ( [ longdistance1.5 ] ) , one writes correction factor as follows : @xmath142 where @xmath143    formula ( [ correctionfactor2.1 ] ) stresses the proportionality of @xmath75 to the laser field intensity ( which is proportional to @xmath144 ) and exhibits explicitly its dependence on the emission angle @xmath69 , which enters , as well as in ( [ cs.4a ] ) , via the second order legendre polynomial .    as a function of @xmath69 the modulus of @xmath75 attains its maximum either at @xmath145 ( and @xmath146 ) or at @xmath147 , depending on the sign of the product @xmath148 and on the relative magnitudes of @xmath149 and @xmath150 . for further reference
let us define the quantity @xmath151 which allows one to estimate the laser field strength needed to achieve a noticeable correction to the angular distribution ( see the discussion in section [ numericalresults ] ) .    substituting ( [ correctionfactor2.1 ] ) into ( [ cs.3 ] ) , taking into account ( [ cs.4a ] ) and expressing @xmath152 via @xmath71 and @xmath153 , one derives the following formula for the differential cross section within the lowest order of perturbation theory in @xmath2 : @xmath154 the coefficients @xmath155 are related to the asymmetry parameter @xmath156 and to @xmath157 via @xmath158 , \quad a_4 = { 18 \over 35}\,e_0 ^ 2\,\beta\,\kappa_2 .
\label{modification2.9}\end{aligned}\ ] ]    the right - hand side of ( [ modification2.3 ] ) suggests that not only the profile of the angular distribution can be changed , but the magnitude of the cross section , integrated over the emission angles , is scaled by the factor @xmath159 which can be greater or lower than one depending on the signs and magnitudes of @xmath156 and @xmath157 .
if the action of @xmath2 is ignored then @xmath160 , @xmath161 and @xmath162 , and equation ( [ modification2.3 ] ) reduces to ( [ cs.4a ] ) .
numerical analysis of the influence of the induced dipole moment on the profile of the angular distribution was carried out for the process of photodetachment of sodium cluster anions na@xmath163 .
the results presented below in this section refer to @xmath164 and @xmath165 .
the corresponding neutral clusters are spherically - symmetric which makes applicable the theory described in section [ framework ] .
the parameters of sodium clusters and their anions , used in our numerical calculations , are summarized in table [ nan.table ] .
let us comment on these data .
the ionic core radius was calculated using the standard relation @xmath166 with the wigner - seitz radius @xmath167 set to 4 a.u .
the values of electron affinities @xmath168 are taken from a recent work by kostko @xcite which contains , in particular , an excellent collection of reference data on various metal clusters .
the other data presented in the table refer to the polarizability of the core and the parameters of the plasmon peak .
we considered the dynamic dipole polarizability @xmath7 of the clusters within the framework of the resonance plasmon approximation ( see , e.g. , @xcite ) : @xmath169 here @xmath170 is the static polarizability , @xmath171 is the surface plasmon energy , and @xmath172 is the linewidth of the plasmon resonance .    due to the spill - out effect the static polarizability of a metallic cluster of a radius @xmath15 exceeds the classical value @xmath173 , which characterizes a metallic ball of the same radius . in the table
we present the ratio @xmath174 calculated within the rpa @xcite and the lda @xcite ) as well as the experimentally measured values @xcite .
the @xmath171 values were obtained with the help of the relation proposed in @xcite to account for the size - dependence of the plasmon energy : @xmath175 , where @xmath176 ev is the mie plasmon energy in a sodium metallic sphere , and @xmath177  takes into account the spill - out effect .
@r|rrrrrrrr & & & & & & + & ( a.u . ) & ( ev ) & ( a.u . ) & @xcite & @xcite & @xcite & ( ev ) & ( ev ) + na@xmath178 & 8.00 & 0.915 & 512 & 1.48 & 1.44 & 1.72 & 2.64 & 0.2 + na@xmath179 & 10.86 & 1.337 & 1280 & 1.41 & 1.37 & 1.67 & 2.81 & 0.4 + na@xmath180 & 13.68 & 1.509 & 2560 & 1.10 & 1.32 & 1.60 & 2.90 & 0.7 + na@xmath181 & 15.48 & 1.719 & 3712 & 1.24 & 1.22 & - & 2.95 & 0.4 + na@xmath182 & 18.06 & 1.847 & 5888 & 1.22 & 1.20 & - & 2.99 & 0.2 +    [ nan.table ]    the size dependence of the width @xmath172 of the plasmon resonance for neutral sodium clusters with @xmath183 up to several hundreds was analyzed in @xcite within the rpa scheme and in @xcite within a semiclassical approach .
these two methods , being consistent for @xmath184 , deviate noticeably for lower @xmath183 values .
the rpa calculations reproduce the @xmath185 scaling law which , is inadequate when compared to the experimental data for @xmath186 ( see the discussion in @xcite ) . the semiclassical calculations @xcite seem to be consistent with the experimental data in the whole region @xmath187 .
the values of @xmath172 , presented in the table , were deduced from figure 1 in @xcite .
it is seen from the table , that for all @xmath183 the plasmon energy exceeds the electron affinity .
this fact justifies the choice of cluster anions in a view of the present study .
indeed , the inequality @xmath188 leads to the additional enhancement of @xmath75 owing to the pronounced increase of the modulus of the dynamic polarizability in the region @xmath189 .    to conclude the discussion on the parameters of the chosen targets
let us note , that both criteria from ( [ functionb3.3a ] ) can be easily met for na@xmath163 .
the first inequality fails only in the narrow near - threshold region @xmath190 but is well fulfilled for @xmath189 .
the validity of the second condition one proves by introducing @xmath191 .     and
for na@xmath192 as indicated .
+ in each graph the upper panel represents the dependence of @xmath193 on @xmath53 ; the middle panel - dependences @xmath194 and @xmath195 , eqs .
( [ correctionfactor2.3])-([correctionfactor2.1a ] ) ; the lower panel - dependence of @xmath196 , eq .
( [ correctionfactor2.1b ] ) .
parameters @xmath171 and @xmath172 are as indicated in table [ nan.table ] .
the data refer to the @xmath197 ratio equal to 1.45 for na@xmath198 and to 1.30 for na@xmath180 .
, title="fig : " ] +   and for na@xmath192 as indicated .
+ in each graph the upper panel represents the dependence of @xmath193 on @xmath53 ; the middle panel - dependences @xmath194 and @xmath195 , eqs .
( [ correctionfactor2.3])-([correctionfactor2.1a ] ) ; the lower panel - dependence of @xmath196 , eq .
( [ correctionfactor2.1b ] ) .
parameters @xmath171 and @xmath172 are as indicated in table [ nan.table ] .
the data refer to the @xmath197 ratio equal to 1.45 for na@xmath198 and to 1.30 for na@xmath180 . ,
title="fig : " ]     but for na@xmath199 and na@xmath200 . for both clusters the ratio @xmath197 is set to 1.20 .
, title="fig : " ] +   but for na@xmath199 and na@xmath200 . for both clusters the ratio @xmath197 is set to 1.20 .
, title="fig : " ]    the results of calculation of the quantities relevant to the correction term @xmath75 are presented in figures [ na08_40.fig ] and [ na58_92.fig ] . for each anion
the calculations were performed for the photon energies above the ionization threshold @xmath168 and in the vicinity of the plasmon peak .
the latter manifests itself as a resonance in the dependence of @xmath193 on @xmath53 , as it is seen from the upper panel in each graph .
the dynamic polarizabilities were calculated via ( [ numerical.1 ] ) with the parameters @xmath171 and @xmath172 as indicated in table [ nan.table ] . given the differences between the calculated and experimentally measured values of the static polarizabilities @xmath170 , as well as the absence of the experimental data for na@xmath181 and na@xmath182 , for the ratio @xmath197 we used the values close to the theoretical data ( see the figures captions ) .
the dependences of the coefficients @xmath149 and @xmath150 on @xmath53 are presented in the middle panels .
they were computed from ( [ correctionfactor2.3 ] ) and ( [ functionb3.3 ] ) . to calculate the right - hand side of the latter formula beyond the plane - wave born approximation one must know the scattering phaseshifts @xmath201 and @xmath202 . to determine the phaseshifts we solved numerically the phase equation ( see , e.g. , @xcite ) for an electron moving in a local potential .
the model potential , which was used to describe the interaction between the photoelectron and the neutral residue , accounted for the following three terms : ( a ) the short - range potential due to the coulomb interaction with the core s electrons and the ionic jellium core , ( b ) the exchange - correlation term , which was treated within the gunnarsson - lundqvist approximation @xcite , ( c ) the polarization potential was considered in the form @xmath203 for @xmath37 and @xmath204 if otherwise .
inclusion of the polarization term is important for an adequate description the low - energy electron scattering from metallic clusters ( e.g. , @xcite ) . a more rigorous treatment accounts for the many - body correlations , which may be considered , for example , within the dyson equation scheme @xcite .
however , for the purpose of our paper it is sufficient to treat the polarization potential in the simple form described above .    comparing the middle panels in the graphs from figures [ na08_40.fig ] and [ na58_92.fig ] one
notices that the behaviour of @xmath194 ( and of @xmath195 as well ) exhibits some general trends along with the specific features determined by the dynamic and geometrical properties of an individual target .
common to all clusters is the existence of extrema for @xmath194 at @xmath189 , which is due to the factor @xmath193 , see ( [ correctionfactor2.3 ] ) .
the extremum can be clearly defined ( see the curves @xmath149 for na@xmath198 and @xmath150 for na@xmath182 which have distinct maxima ) or less pronounced ( both curves for na@xmath181 ) or barely seen ( the case of na@xmath180 ) .
the variation in the shape of the extremum , as well as the behaviour of the curves in the regions below and above the plasmon resonance , can be understood by analyzing the @xmath53-dependence of the right - hand sides in ( [ correctionfactor2.3 ] ) and ( [ functionb3.3 ] ) .
the non - oscillatory part of the latter depends on the photon energy mainly through the factor @xmath205 .
this leads to @xmath206 .
therefore , a symmetric resonance profile of @xmath193 is distorted due to the factor @xmath207 . for a powerful and narrow resonance ( as in the case of na@xmath182 )
the distortion is comparatively small , whereas for a wide resonance the disproportion between the low- and the high-@xmath53 shoulders may completely change the behaviour of the curve in the vicinity of the resonance , as it happens for na@xmath180 . in either case the factor @xmath207 results in a faster decrease of @xmath208 with @xmath53 beyond the resonance . in the region
@xmath209 the curves @xmath194 , being enhanced by the factor , exhibit a non - monotonous oscillatory character ( most clearly seen in na@xmath180 , na@xmath181 and na@xmath182 graphs ) .
this is due to the second and the third terms on the right - hand side of ( [ functionb3.3 ] ) , which contain the oscillating factors @xmath210 and @xmath211 , dependent on the cluster size , and the factors @xmath212 , @xmath213 , dependent on the phaseshifts .
the solid curves in the middle panels of each graph represent the dependences @xmath195 defined in ( [ correctionfactor2.1a ] ) . for a given photon energy the quantity @xmath214 stands for the maximum value of @xmath215 within the interval @xmath216 $ ] .
it was already mentioned , that the correction term must satisfy the condition @xmath77 .
thus , taking into account that @xmath217 , one can introduce a critical value @xmath218 , which defines the upper boundary for the laser field strength consistent with the perturbative approach adopted in this paper : @xmath219 the corresponding intensity of the laser field ( in w/@xmath1 ) can be calculated from the relation @xmath220 , where @xmath221 is measured the atomic units .    in the limit
@xmath222 a strong inequality @xmath223 is valid , so that the photoionization process is not affected by the induced dipole moment .
the influence of the latter can become pronounced for @xmath224 resulting in @xmath79 .
@l|crrccccc & @xmath225 & & & @xmath226 & @xmath227 & @xmath228 & @xmath229 & @xmath230 + & ( ev ) & ( @xmath231 ) & ( @xmath231 ) & ( @xmath231 ) & ( @xmath232 a.u . ) & ( @xmath233w/@xmath1 ) & ( @xmath232 a.u . ) & ( @xmath232 ev ) + na@xmath178 & 2.64 & -0.470 & -0.031 & 0.501 & 1.41 & 0.66 & 17.5 & 1.44 + na@xmath179 & 2.73 & 0.021 & 0.022 & 0.043 & 4.82 & 7.67 & 30.8 & 15.7 + na@xmath180 & 2.72 & 0.022 & -0.100 & 0.078 & 3.58 & 4.23 & 37.0 & 8.64 + na@xmath181 & 2.89 & 0.151 & -0.160 & 0.231 & 2.08 & 1.43 & 44.9 & 2.61 + na@xmath182 & 2.98 & 0.128 & -1.589 & 1.461 & 0.83 & 0.23 & 50.1 & 0.39 +    [ nan.table2 ]    dependences @xmath218 are plotted in the lower panels in figures [ na08_40.fig ] and [ na58_92.fig ] . for all clusters
this function attains its minimum in the point @xmath234 in the vicinity of the plasmon resonance .
the values of @xmath225 , as well as of other quantities calculated in this point , are presented in table [ nan.table2 ] .
it is instructive to compare the critical field @xmath235 with the strength @xmath229 of the internal anionic field .
the latter can be related to the ionization potential @xmath168 through @xmath236 .
it is seen that @xmath229 exceeds @xmath235 by the order of magnitude .
the strong inequality @xmath237 , being in accordance with the first condition from ( [ conditions ] ) , justifies the approximation adopted in this paper to neglect the modification of the ground state orbitals due to the laser field .
the last column in the table contains the values of ponderomotive energy of the photoelectron in the critical laser field . comparing the ponderomotive shifts and the kinetic energies @xmath238 ev
, one notices that the latter is larger by several orders of magnitude .
therefore , the second condition in ( [ conditions ] ) is also well fulfilled .
hence , one can expect , that for the intensities @xmath239 w/@xmath1 the modification of the spectral - angular distribution due to the action of the field of the induced dipole moment can be quite pronounced . to carry out the corresponding quantitative estimates one can use the formulae obtained within the framework of perturbation theory .     of the angular distribution ( [ profile.1 ] ) calculated for different sodium cluster anions and for the three values of asymmetry parameter @xmath156 ( as indicated ) , and for several values of the ratio @xmath240 ( see the legend in the graph @xmath241 for na@xmath242 ) .
the data refer to the photon energies @xmath225 indicated in table [ nan.table2 ] .
thick black curves in each graph represent the field free distribution @xmath243 , eq .
( [ profile.0 ] ) .
the length of line segment connecting the origin and a given point on the curves is equal to the value of @xmath244 or @xmath243 ( in absolute units ) in the corresponding direction .
the emission angle @xmath69 is measured with respect to @xmath22 , which is shown in the graphs with @xmath245 . in each graph
the spacing between ticks is equal to 1 for both directions .
, title="fig : " ] +   of the angular distribution ( [ profile.1 ] ) calculated for different sodium cluster anions and for the three values of asymmetry parameter @xmath156 ( as indicated ) , and for several values of the ratio @xmath240 ( see the legend in the graph @xmath241 for na@xmath242 ) .
the data refer to the photon energies @xmath225 indicated in table [ nan.table2 ] .
thick black curves in each graph represent the field free distribution @xmath243 , eq .
( [ profile.0 ] ) .
the length of line segment connecting the origin and a given point on the curves is equal to the value of @xmath244 or @xmath243 ( in absolute units ) in the corresponding direction .
the emission angle @xmath69 is measured with respect to @xmath22 , which is shown in the graphs with @xmath245 . in each graph
the spacing between ticks is equal to 1 for both directions .
, title="fig : " ] +   of the angular distribution ( [ profile.1 ] ) calculated for different sodium cluster anions and for the three values of asymmetry parameter @xmath156 ( as indicated ) , and for several values of the ratio @xmath240 ( see the legend in the graph @xmath241 for na@xmath242 ) .
the data refer to the photon energies @xmath225 indicated in table [ nan.table2 ] .
thick black curves in each graph represent the field free distribution @xmath243 , eq .
( [ profile.0 ] ) .
the length of line segment connecting the origin and a given point on the curves is equal to the value of @xmath244 or @xmath243 ( in absolute units ) in the corresponding direction .
the emission angle @xmath69 is measured with respect to @xmath22 , which is shown in the graphs with @xmath245 . in each graph
the spacing between ticks is equal to 1 for both directions .
, title="fig : " ] +   of the angular distribution ( [ profile.1 ] ) calculated for different sodium cluster anions and for the three values of asymmetry parameter @xmath156 ( as indicated ) , and for several values of the ratio @xmath240 ( see the legend in the graph @xmath241 for na@xmath242 ) .
the data refer to the photon energies @xmath225 indicated in table [ nan.table2 ] .
thick black curves in each graph represent the field free distribution @xmath243 , eq .
( [ profile.0 ] ) .
the length of line segment connecting the origin and a given point on the curves is equal to the value of @xmath244 or @xmath243 ( in absolute units ) in the corresponding direction .
the emission angle @xmath69 is measured with respect to @xmath22 , which is shown in the graphs with @xmath245 . in each graph
the spacing between ticks is equal to 1 for both directions .
, title="fig : " ] +   of the angular distribution ( [ profile.1 ] ) calculated for different sodium cluster anions and for the three values of asymmetry parameter @xmath156 ( as indicated ) , and for several values of the ratio @xmath240 ( see the legend in the graph @xmath241 for na@xmath242 ) .
the data refer to the photon energies @xmath225 indicated in table [ nan.table2 ] .
thick black curves in each graph represent the field free distribution @xmath243 , eq .
( [ profile.0 ] ) .
the length of line segment connecting the origin and a given point on the curves is equal to the value of @xmath244 or @xmath243 ( in absolute units ) in the corresponding direction .
the emission angle @xmath69 is measured with respect to @xmath22 , which is shown in the graphs with @xmath245 . in each graph
the spacing between ticks is equal to 1 for both directions .
, title="fig : " ]    the calculated profiles of the angular distributions are presented in figure [ profiles.fig ] . in the present paper
we have not made an attempt to compute the angular asymmetry parameter @xmath246 for sodium anion clusters ( neither have we found its values in the literature ) .
instead , for each target we carried out the calculations for the three values , @xmath247 , aiming to demonstrate that the modification of the profile takes place for _ any _ @xmath156 from the allowable interval @xmath248 $ ] .    in each graph from figure [ profiles.fig ]
the thick solid curve represents the unperturbed profile @xmath249 the emission angle @xmath69 is measured with respect to the vector @xmath22 , which is indicated in the graphs with @xmath245 .
other curves correspond to the modified profiles @xmath250 ( see ( [ cs.3 ] ) , ( [ correctionfactor2.1 ] ) and ( [ correctionfactor2.1b ] ) ) calculated for several values of the ratio @xmath251 . for each anion
the data refer to the photon energy @xmath225 .
figure [ profiles.fig ] demonstrates , that even for comparatively low intensities ( corresponding to @xmath252 ) the changes in the profiles are already visible . for @xmath253
they become much more pronounced .
it is seen that variation of the profile depends , both qualitatively and quantitatively , on the value of @xmath156 and on the choice of the target anion which enters via the coefficients @xmath254 and @xmath255 .    in the case
@xmath245 the field - free profile @xmath243 is proportional to @xmath256 yielding @xmath257 as the minimum and @xmath258 as the maximum values within the interval @xmath259 $ ] .
the modified profile can be written as @xmath260 with @xmath261 and @xmath262 .
one easily shows , that for @xmath263 the function @xmath264 acquires additional extremum .
however , this inequality can not be met for the considered targets ( see the values of @xmath254 , @xmath255 and @xmath265 in table [ nan.table2 ] ) and for @xmath253 .
therefore , the pattern of the modified angular distribution is similar to that obtained without the influence of the induced dipole .
namely , the strongest emission occurs in the directions perpendicular to the field @xmath22 and there is no emission for @xmath145 and @xmath266 .
on the whole , the distortion increases with @xmath267 , attaining the largest values at @xmath147 .
the trend of the largest distortion depends on the sign of @xmath254 : if @xmath268 ( as in the cases of all anions but na@xmath178 ) then @xmath269 whereas negative @xmath254 leads to a shrinking of the modified profiles ( as it happens for na@xmath242 ) .    for @xmath270
the unperturbed angular distribution is isotropic .
the modified profile reads as @xmath271 , with @xmath272 and @xmath273 defined as above .
the way the profile evolves with the increase of @xmath267 to a great extent depends on the sign of @xmath255 . for @xmath274
the emission along @xmath275 dominates over the emission in the perpendicular directions ( see the case of na@xmath276 ) , whereas for negative @xmath255 the perpendicular direction becomes preferable .
the graphs from the central column in figure [ profiles.fig ] demonstrate , that for @xmath277 the distortion can become dramatic both in form and in magnitude
.    in the limit @xmath241 the modified profile behaves as @xmath278 .
the additional extrema appear if the following two conditions , @xmath279 and @xmath280 , are met .
these conditions can be satisfied in the limit @xmath281 for all considered anions but na@xmath276 .
most explicitly it is seen in na@xmath192 and na@xmath200 graphs , where the curves with @xmath282 represent the four - petal angular distributions with the maximum emission at @xmath283 and @xmath284 .
the same pattern characterizes na@xmath285 , although in this case it is hardly visible due to a strong decrease in the magnitude of @xmath286 for @xmath282 .
( see ( [ spectra.1 ] ) ) calculated for different sodium cluster anions and for the three values of asymmetry parameter @xmath156 ( as indicated in the top left graph ) .
the data refer to the laser field intensity @xmath287 w/@xmath1 . ]    finally , let us briefly discuss the influence of the induced dipole moment on the cross section integrated over the emission angles . recalling ( [ modification2.3 ] ) and ( [ modification2.9 ] ) , one finds the following formula for the ratio of the cross sections with , @xmath288 , and without , @xmath289 , account for the correction : @xmath290 figure [ spectra.fig ] presents the ratios @xmath291 versus @xmath53 . for each anion
the ratio was calculated using the dependences @xmath149 and @xmath149 , presented in figures [ na08_40.fig ] and [ na58_92.fig ] , and for three values of the asymmetry parameter ( as indicated ) .
the intensity of the laser field was set to @xmath287 w/@xmath1 , which corresponds to the field strength @xmath292 a.u . comparing the latter value to the critical field @xmath218 ( see figures [ na08_40.fig ] and [ na58_92.fig ] ) one notices , that for all targets and in the whole range of @xmath53 the parameter @xmath293 is less than one .
thus , formula ( [ spectra.1 ] ) is applicable for the quantitative estimate of the correction due to the induced moment .
as one can see from figure [ spectra.fig ] , the largest deviation ( up to @xmath294 ) of @xmath288 from @xmath289 occurs for na@xmath285 and na@xmath200 in the vicinities of the plasmon resonances , where @xmath267 reaches its maximum values : @xmath295 for na@xmath285 and @xmath296 for na@xmath200 . in the case of na@xmath200
a strong inequality @xmath297 is valid ( see figure [ na58_92.fig ] or / and table [ nan.table2 ] ) .
therefore , the magnitude of the correction term from the right - hand side of ( [ spectra.1 ] ) becomes sensitive to the choice of the asymmetry parameter . as a result
, the ratio @xmath298 strongly varies with @xmath156 . on the contrary , for na@xmath285
the relation @xmath299 is valid , leading to a weak dependence of the ratio on @xmath156 .    for the anions na@xmath192 and na@xmath199 the largest values of @xmath267 equal to 0.05 and 0.15 , correspondingly .
hence , the intensity @xmath300 is small compared to @xmath301 . as a result , the deviation of @xmath288 from the field - free value @xmath289 is less pronounced .
in summary , we have demonstrated that the dipole moment , induced at the residue by the incoming laser field , can strongly influence the photoionization process .
the additional time - dependent long - range field due to the induced moment , acting on the escaping electron , modifies the spectral - angular and spectral distributions of photoelectrons .
the effect itself , as well as its quantitative treatment , is quite sensitive to the choice of the target and to the parameters of the laser field .
we have demonstrated , that for a strongly polarizable target ( e.g. , metallic cluster anions ) strong modifications of the characteristics of a single - photon ionization process in the vicinity of the plasmon resonance ( typically , @xmath302 ev ) can be achieved by applying laser fields of comparatively low intensities @xmath0 w/@xmath1 .    from a theoretical viewpoint ,
the advantage of this regime is that the field of induced dipole moment can be treated perturbatively , whereas the influence of the laser field on the photoelectron ( a ponderomotive effect ) and on the ground state of the target can be ignored .
we have shown , that even in this weak - field regime one can expect dramatic changes in the profile of angular distribution of photoelectrons as well as in the @xmath53 dependence of the photoionization cross section . to a great extent , these changes depend on the parameters of the target , like its radius , the static and dynamic polarizabilities , the energy of the plasmon resonance peak and the affinity ( or , the ionization potential in the case of a neutral target ) .
additional variety can be obtained by changing the intensity and the frequency of the laser field .
the analysis , which we have carried out , is based on simple models adopted for the description of the photodetachment process ( the plasmon resonance approximation ) and of the interaction of the escaping electron with the static field of the residue ( via the lda + long - range polarization potential ) .
we have also restricted ourselves to the case of spherically - symmetric targets .
a more accurate quantitative treatment must include more sophisticated approaches , which take into account the many - body correlation effects intrinsic to both the photodetachment and the electron escape . however , these will leave untouched the physics behind the basic phenomenon discussed in our paper .
it is possible to study the predicted phenomena by means of modern experimental techniques .
the latter include the production and manipulation of the beams of cluster anions of various types and sizes , the energy and angle resolved photoelectron spectroscopy , the availability of the laser fields of various intensity , frequency and pulse duration @xcite .
we believe , that with all these components at hand , soon it will become possible to test the theory against experimental data .
we are grateful to roman polozkov for the help in numerical calculation of the ground states of sodium clusters , and to steven manson and bernd von issendorff for the helpful discussions .
this work was supported by the european commission within the network of excellence project excell ( project number 515703 ) .
99    p.a .
golovinski , i.yu .
kiyan , 1990 , _ sov .
- uspekhi _ * 33 * , 453 .
chernov , i.yu .
kiyan , h.  helm , b.a .
zon , _ phys .
rev . _ a * 71 * , 033410 ( 2005 ) .
gribakin , m.yu .
kuchiev , _ phys .
rev . _ a * 55 * , 3760 ( 1997 ) .
o.  kostko , c.  bartels , j.  schwbel , c.  hock , b. v.  issendorff , _ j. phys .
series _ * 88 * , 012034 ( 2007 ) .
o.  kostko , 2007 , _ photoelectron spectroscopy of mass - selected sodium , coinage metal and divalent metal cluster anions .
_ inaugural - dissertation ( universitt freiburg ) .
gordon w.  f.  drake ( ed . ) , 2006 , _ springer handbook of atomic , molecular and optical physics _
( germany : springer ) .