in a couple of experiments , the measurement of the frequencies of collective modes in trapped fermi gases revealed a lot of interesting information on the equation of state , the validity of superfluid hydrodynamics , and the superfluid - normal phase transition @xcite .
more recently , the transition from the hydrodynamic to the collisionless regime in the normal phase was also studied @xcite . however , since the frequencies of the collective modes depend on the trap frequencies , a precise knowledge of the latter is required for a meaningful interpretation of the collective - mode data .    a possibility to determine the trap frequency with high precision is the measurement of the frequency of the sloshing mode @xcite , which is an oscillation of the center of mass of the trapped atom cloud . in a harmonic trap , this oscillation ( also called kohn mode ) is undamped and its frequency coincides exactly with the trap frequency in the corresponding direction , independently of the number of atoms @xmath0 , or their interaction , their temperature @xmath1 etc .
this is a consequence of the kohn theorem @xcite and follows from the fact that the center - of - mass oscillation decouples completely from the internal dynamics of the gas if the interaction is translationally invariant and the external potential is harmonic @xcite .
however , in practice the trap potential is never exactly harmonic . in optical dipole traps
@xcite , the potential is typically gaussian ( corresponding to the intensity profile of the laser beam ) . because of this anharmonicity , the frequency of the sloshing mode
is shifted , and the shift depends on the system parameters such as @xmath0 , @xmath1 , the scattering length @xmath2 characterizing the interaction strength , etc .
furthermore , the sloshing mode is no longer undamped .
the anharmonicity of the trap complicates considerably the analysis of collective - mode experiments . in ref .
@xcite the measured frequencies were corrected for the anharmonicity effects by giving them in units of the measured frequency of the sloshing mode .
however , it is clear that the damping rates of the modes can not be corrected in this way . in ref .
@xcite , the damping rate of the sloshing mode was used to estimate the increase of the damping rate of other modes due to the anharmonicity .
both corrections are ad - hoc prescriptions without rigorous justification .
it is therefore strongly desirable to get a better understanding of the anharmonicity effects on the sloshing mode .
the aim of the present paper is to describe theoretically the frequency shift and the damping rate of the sloshing mode in an anharmonic trap .
we will compare our results with the experimental data available from the innsbruck group @xcite and with the numerical results by wu and zhang @xcite .
we will also discuss in detail the damping mechanism of the sloshing mode . in an anharmonic potential ,
the center - of - mass motion is no longer decoupled from the internal degrees of freedom of the cloud .
we will see that the damping of the transverse sloshing mode is a consequence of its coupling to other damped collective modes , in particular to the radial dipole mode and the bending mode .
the framework of our study is the boltzmann equation , including mean field @xcite and in - medium cross - section @xcite .
especially the mean field is expected to be important in the present context , because it can have a sizable effect on the density profile , i.e. , on how far the cloud extends into the anharmonic region of the trap potential .
the boltzmann equation is solved approximately with the help of the phase - space moments method .
this method , when extended beyond the lowest order , has proven to be in very good quantitative agreement with the results of a full numerical simulation @xcite .
it has also been quite successful for the description of the experimental results for the frequency and damping rate of the radial quadrupole mode @xcite . in the present case of the sloshing mode
, we include phase - space moments of first and third order .
the paper is organized as follows .
the general formalism is briefly presented in sec .
[ sec : formalism ] .
then we specify our model for the experimental trap potential in sec .
[ sec : potential ] . in sec .
[ sec : firstorder ] , we give a formula for the frequency shift of the sloshing mode within the first - order moments method .
then we extend the ansatz to third order in sec .
[ sec : thirdorder ] .
the physical contents of the extended ansatz and the numerical results are discussed . in sec .
[ sec : cubic ] the third - order ansatz is used to describe also the radial dipole , radial octupole , and bending modes . finally , in sec .
[ sec : cl ] , we will conclude .    throughout the paper
, we use units with @xmath3 .
in this section , we give a short summary of the formalisms of refs .
more details can be found there .      we consider a balanced two - component ( @xmath4 ) fermi gas of atoms with mass @xmath5 and interspecies attractive interaction ( scattering length @xmath6 ) , trapped in a potential @xmath7 .
the framework we use to describe the collective dynamics of the system in the normal - fluid phase is the boltzmann equation .
we assume that the two components move in phase , so that only one distribution function @xmath8 is needed .
it is normalized to @xmath9 , and expectation values of one - body operators are given by @xmath10 as in ref .
@xcite , we include in - medium effects in both the transport and the collision parts of the boltzmann equation : a mean field like potential @xmath11 and the in - medium modified cross - section @xmath12 , respectively .    within the thomas - fermi or local - density approximation ( lda )
, the equilibrium distribution function reads @xmath13 } + 1},\ ] ] where @xmath14 is the mean field in equilibrium , @xmath15 is the chemical potential , and @xmath16 is the inverse temperature . as in ref .
@xcite , we obtain @xmath14 from the single - particle self - energy in ladder approximation , evaluated at the fermi level . in the weak coupling limit
, this reduces to the hartree term @xmath17 , where @xmath18 is the density per spin state , and it remains finite for all interaction strengths up to the unitary limit , @xmath19 . in equilibrium ,
the main effect of the mean field is to enhance the density in the center of the trap , as shown in fig . 3 of ref .
@xcite .    for the study of collective oscillations , it is sufficient to consider small deviations from equilibrium and to linearize the boltzmann equation with respect to @xmath20
. if we write the variation of the distribution function in the form @xcite @xmath21 with @xmath22 , the linearized boltzmann equation can be written as @xmath23\ , , \label{eq : boltzlin}\end{gathered}\ ] ] where @xmath24 denotes the poisson bracket , @xmath25 is the perturbation of the trap potential that is used to excite the collective mode , @xmath26 is the variation of the mean field due to the variation of the density , and @xmath27 $ ] is the linearized collision term . since we want to calculate the so - called response function , we take the perturbation to be a pulse , @xmath28 as in ref .
@xcite , we approximate the variation of the mean field by @xmath29 the linearized collision integral reads @xmath30 = \int \frac{d^3p_1}{(2\pi)^3 } \int d \omega\ ,     \frac{d\sigma}{d\omega } \frac{|{{\bm{\mathrm{p}}}}-{{\bm{\mathrm{p}}}}_1|}{m } { f_{\mathit{eq}}}f_{{\mathit{eq}}1 }     { \bar{f}_{\mathit{eq}}}^\prime \bar{f}_{{\mathit{eq}}1}^\prime\\      \times ( \phi+\phi_1-\phi^\prime-\phi_1^\prime)\,,\end{gathered}\ ] ] where an obvious notation for the different @xmath31 and @xmath32 at the momenta before ( @xmath33 , @xmath34 ) and after the collision ( @xmath35 , @xmath36 ) has been used . note that , especially near the critical temperature , the in - medium cross - section @xmath12 can differ strongly from the free one @xcite .      as in ref .
@xcite , we are looking for a semi - analytical solution of the boltzmann equation ( [ eq : boltzlin ] ) by using the method of phase - space moments . in this section
, we will generalize the formalism of that paper to the case of an arbitrary trap potential @xmath37 and with mean field @xmath11 .
the basic idea is to approximate the function @xmath32 by a polynomial in the components of @xmath38 and @xmath33 with time - dependent coefficients @xmath39 , @xmath40 where the @xmath41 are suitable basis functions , e.g. , monomials in the components of @xmath38 and @xmath33 .
multiplying the linearized boltzmann equation ( [ eq : boltzlin ] ) by @xmath42 and integrating over phase space , one obtains , after a fourier transform with respect to @xmath43 , a set of @xmath44 coupled linear algebraic equations for the @xmath44 coefficients @xmath45 . in matrix form , they read @xmath46 where @xmath47\ , , \label{eq : acollij}\end{gathered}\ ] ] and @xmath48 the variation of the mean field in eq .
( [ eq : adeltauij ] ) is defined as @xmath49 equations ( [ eq : atransij ] ) and ( [ eq : ai ] ) correspond to eqs .
( 14 ) and ( 16 ) in ref .
@xcite which have been simplified by integration by parts . comparing eqs .
( [ eq : ac = a])-([eq : ai ] ) with the analogous ones of ref .
@xcite , eqs .
( 11)-(16 ) , one sees that the mean field gives rise to a new term , @xmath50 .
however , one should keep in mind that implicitly all terms depend on the mean field since it modifies the equilibrium distribution @xmath51 .
one can see that the matrices @xmath52 and @xmath53 are symmetric , @xmath54 is antisymmetric , whereas @xmath50 has no defined symmetry . in practice , the calculation of the matrices is straight - forward , but tedious , and we made use of the mathematica software to express the numerous matrix elements in terms of a smaller number of integrals over equilibrium quantities .      without an external perturbation @xmath55 , i.e. , for @xmath56 , eq . ( [ eq : ac = a ] ) has a solution with non - vanishing coefficients @xmath39 only if @xmath57 . the frequencies @xmath58 for which this happens
are obviously given by the eigenvalues of the matrix @xmath59 .
if they are well separated , it is possible to interpret them as the frequencies of the eigenmodes of the system . in general , they are complex , and their imaginary part describes the damping rate of the corresponding mode @xcite .
however , as discussed in ref .
@xcite , when the moments method is extended to higher order , there can be many eigenvalues belonging to a single collective mode . in this case , the scattering of the eigenvalues , which goes over into a continuous spectrum in the limit of an infinite number of moments @xcite , corresponds to a new contribution to the damping ( landau damping ) in addition to the imaginary parts coming from the collision term . in order to obtain the mode frequency and damping rate in this case
, it is useful to look at the response function which contains the contributions of all eigenvalues .
we denote by @xmath60 the fourier transform of the expectation value @xmath61 of some operator @xmath62 after the perturbation .
the so - called response function is equal to @xmath63 in the special case that the excitation operator is @xmath64 . in all excitations considered in this paper
, @xmath65 will be zero .
the strength function is proportional to the imaginary part of the response function .    in order to calculate the response function ,
we need the coefficients @xmath45 of the ansatz ( [ generalansatz ] ) . by diagonalizing the matrix @xmath66 with @xmath67
, we can write them as @xmath68 then , the response function can be easily obtained as @xmath69 , where @xmath70    the frequency and damping of the mode can be extracted from the response @xmath60 , e.g. , by fitting the peak in the strength function @xmath71 corresponding to the mode under consideration with a lorentzian .    using the solution for the coefficients @xmath45 at the peak
, we can also obtain the velocity field of the corresponding collective mode : @xmath72
in this section we recall the shape of a typical optical dipole trap .
we concentrate on the configuration used by the innsbruck group , consisting of a focused - beam trap @xcite with additional magnetic confinement in the axial direction . following ref .
@xcite , we parametrize the experimental trap potential as @xmath73\\    + \frac{1}{2 } m { \omega_{z,\mathit{mag}}}^2 z^2 \ , , \label{eq : vtrap}\end{gathered}\ ] ] where @xmath74 , @xmath75 is the trap depth , @xmath76 is the minimal waist of the laser beam , @xmath77 is the rayleigh length , @xmath78 the laser wavelength and @xmath79 defines the magnetic trapping in the @xmath80 direction .    in order to compare our results with the data from refs .
@xcite , we focus on the setup of that experiment .
the trap was of the type ( [ eq : vtrap ] ) , with strong anharmonicity in the @xmath81-@xmath82 plane due to the gaussian shape of the laser beams .
the confinement along the @xmath80 direction was practically harmonic . in our calculations
, we will expand eq .
( [ eq : vtrap ] ) up to second order in @xmath80 and up to sixth order in @xmath83 .
the result can be written as @xmath84 where @xmath85 and @xmath86 are the radial and axial trap frequencies , respectively .
we use eq .
( [ eq : vsix ] ) instead of eq .
( [ eq : vtrap ] ) because it simplifies the calculations in the sense that it leads to a density that tends to zero for @xmath87 .
the necessity to go beyond fourth order is also clear from eq .
( [ eq : vsix ] ) : otherwise the potential would be unbound from below . for illustration , fig .
[ fig : vtrap ]     for @xmath88 .
the fermi energy for the parameters of ref .
@xcite is indicated by the thin dashes .
the potentials are in units of @xmath89 and @xmath81 is in units @xmath90 . ]
shows the real trap potential ( [ eq : vtrap ] ) , its harmonic approximation , and the one we will use , eq .
( [ eq : vsix ] ) .
in the present paper , we concentrate on axially symmetric traps . nevertheless , we write @xmath91 and @xmath92 instead of @xmath93 for the trap frequencies in @xmath81 and @xmath82 direction if the formulas can be generalized to the triaxial case .    in our numerical calculations , we use as an example the parameters of the innsbruck experiment @xcite .
as mentioned in ref .
@xcite , the sloshing mode was studied with the same parameters as the compression mode , i.e. , @xmath94 hz , @xmath95 hz , and @xmath96 @xmath97k .
the trap was loaded with @xmath98 atoms of @xmath99li in the unitary limit , @xmath100 , and the temperature was varied between @xmath101 and @xmath102 . since we can not describe the superfluid phase , we limit ourselves to temperatures above @xmath103 @xcite .
moreover , we approximate the unitary limit numerically by setting @xmath104 .    as pointed out in ref .
@xcite , the anharmonicity effects depend mainly upon the ratio @xmath105 . defining the fermi energy as usual by @xmath106 ,
one obtains @xmath107 , i.e. , at low temperature the anharmonicity effects are relatively weak . for illustration , @xmath108
is indicated in fig .
[ fig : vtrap ] by the thin dashed line .
we see that , for the present choice of parameters , the atoms start to feel the anharmonicity of the potential when their energy exceeds @xmath108 , while the sixth - order approximation to the potential stays very precise up to about five times the fermi energy .
in order to describe the sloshing motion , say , along the @xmath81 direction , the ansatz ( [ generalansatz ] ) has to contain at least two basis functions : @xmath109 the second term describes the collective velocity of the cloud in @xmath81 direction , while the first one corresponds approximately to a displacement of the center of mass .
the latter statement becomes exact in the special case of a purely harmonic trap without mean field ( @xmath110 ) . in this case , the ansatz ( [ ansatzfirst ] ) is closed with respect to the operator @xmath111 on the left - hand side of the boltzmann equation ( [ eq : boltzlin ] ) , and since the collision term does not contribute ( @xmath112=i[p_x]=0 $ ] because of momentum conservation in a collision ) , the ansatz ( [ ansatzfirst ] ) becomes exact .
however , if the trap is not purely harmonic , this is no longer true , since the gradient of the trap potential generates new terms . and
even if the trap is harmonic , it is not true for interacting atoms , which feel a mean field @xmath113 in addition to the trap potential @xmath37 . in the next section
, we will therefore extend the ansatz ( [ ansatzfirst ] ) , but let us first look what happens at first order .
we now apply the formalism described in sec .
[ sec : formalism ] to @xmath114 .
the calculation of the matrix @xmath115 is very easy in this case .
after some algebra , we obtain a quadratic equation for the frequency of the sloshing mode by imposing @xmath116 .
its solution reads @xmath117 the parameters @xmath118 and @xmath119 characterize , respectively , the strength of the mean field and of the anharmonicity effects . equation ( [ eq : wk1st ] ) has two important features : ( a ) if @xmath37 is harmonic ( @xmath120 ) , then @xmath121 , independently of the interaction , in accordance with the kohn theorem @xcite .
this point was already discussed in ref .
@xcite and shows the consistency of our approach , in particular of eq .
( [ deltau ] ) .
( b ) if @xmath37 is not purely harmonic ( @xmath122 ) , then @xmath123 depends on both the anharmonicity of the trap and on the interaction .
another property of eq .
( [ eq : wk1st ] ) is that the sloshing mode is undamped ( i.e. , @xmath123 is real ) even in the case of an anharmonic potential .
this is , however , only a consequence of the first - order approximation .    a similar calculation was performed in the appendix b of ref .
@xcite . in that reference ,
only the leading anharmonicity correction was kept . if we combine eqs .
( [ eq : vsix ] ) , ( [ eq : wk1st ] ) , and ( [ eq : cse ] ) and keep only the leading order in @xmath124 , we find @xmath125 this result differs from eq .
( b2 ) of ref .
@xcite in two respects .
first , there is an additional factor of three in front of the @xmath126 term ( misprint in ref .
second , the anharmonicity correction is enhanced by a factor @xmath127 .
this factor is missing in ref .
@xcite because there the mean field was not considered , although the expectation values @xmath126 and @xmath128 were calculated with the density profiles of an interacting gas .
as in ref .
@xcite , we will extend the ansatz ( [ ansatzfirst ] ) by including higher - order terms . at the next higher order
, the ansatz contains 18 terms @xmath129 where : @xmath130 the first two terms of @xmath131 coincide with @xmath114 , and the subsequent ones are all possible terms of third order in the components of @xmath38 and @xmath33 satisfying the symmetries of the mode : odd under @xmath132 and even under @xmath133 and @xmath134 .    note that some of the third - order terms , e.g. , @xmath135 , generate a velocity field that is quadratic in the coordinates , whereas the velocity field corresponding to the first - order ansatz is spatially constant .
other terms , e.g. , @xmath136 or @xmath137 , describe momentum - sphere distortions .
they give rise to a non - vanishing collision integral so that they generate a damping .
the momentum - independent third - order terms such as @xmath138 describe deformations of the cloud in coordinate space that will be discussed in more detail in the next subsection .    as explained in sec .
[ sec : response ] , the frequency and damping rate should now be extracted from the response function .
the suitable excitation operator is @xmath139 because it gives all atoms a constant kick in @xmath81 direction , @xmath140 , at the moment of the excitation @xmath141 .
the response function is defined as the expectation value @xmath142 .
we checked numerically that , in a purely harmonic potential @xmath143 , the strength function has only a single sharp peak ( no damping ) at @xmath144 within the precision of our calculation ( better than @xmath145 ) , independently of the temperature and of the interaction strength .    in fig .
[ fig : resp_sloshing ]     for the sloshing mode ( excitation @xmath146 ) in the anharmonic trap for various temperatures .
the system parameters are those of ref .
[ fig : resp_sloshing ] ]    we show the strength functions for the sloshing mode in the anharmonic trap for various temperatures . in all cases
the strength is concentrated in a single peak and the width of the peak corresponds to the damping rate .
the position of the peak is always below @xmath144 , and with increasing temperature , this shift gets stronger .
this is easily understood : with increasing temperature , more and more atoms reach the peripheral region where the anharmonic potential is flatter than the harmonic one ( cf .
[ fig : vtrap ] ) .    since in all the cases shown in fig .
[ fig : resp_sloshing ] the strength is concentrated in a single peak , this allows us to extract the frequency and damping rate by fitting it with a single lorentzian .
( this is not clear a priori : for instance , if we do the calculation without mean field , the strength function for low temperature has two separate peaks , and it is not evident how one has to define the average frequency and total width . ) the results for the frequency and damping rate as functions of temperature are shown in fig .
[ fig : sloshing](a )    . solid line : third - order result ; dashed line : first - order result ; dash - dotted line : first - order result obtained without mean field .
the experimental data points are taken from ref .
( b ) temperature dependence of the damping rate in units of @xmath91 .
solid line : full third - order result ; dashed line : third - order result obtained with the free instead of the in - medium cross section .
the filled and empty data points are measured damping rates of the sloshing modes in @xmath81 and @xmath82 direction , respectively , from ref .
[ fig : sloshing],title="fig : " ] . solid line : third - order result ; dashed line : first - order result ; dash - dotted line : first - order result obtained without mean field .
the experimental data points are taken from ref .
( b ) temperature dependence of the damping rate in units of @xmath91 .
solid line : full third - order result ; dashed line : third - order result obtained with the free instead of the in - medium cross section .
the filled and empty data points are measured damping rates of the sloshing modes in @xmath81 and @xmath82 direction , respectively , from ref .
[ fig : sloshing],title="fig : " ]    and ( b ) as the solid lines . for comparison , the dashed line in fig .
[ fig : sloshing](a ) shows the first - order result for the frequency .
the dash - dotted line corresponds to the first - order result obtained without mean field .
as expected , the mean field reduces the shift of the frequency , but its effect is moderate and most pronounced at lower temperature .    the experimental results from fig . 2 of ref .
@xcite are also shown in fig .
[ fig : sloshing](a ) . in that paper ,
the measured empirical temperatures @xmath147 were already converted into real temperatures @xmath148 as explained in @xcite , so that the results can immediately be compared with our calculation .
the trend is correctly reproduced , and the third - order result is clearly in better agreement with the data than the first - order one , but the theoretical frequency shift is still stronger than the experimental shifts reported in ref .
however , one should keep in mind that while @xmath149 was precisely determined , the ratio @xmath150 depends also on @xmath89 , which was not very well known in that experiment .
actually , @xmath89 was deduced from @xmath149 under the assumption that the anharmonicity effects become negligible at @xmath151 @xcite . if we suppose that the true trap frequency @xmath89 was just @xmath152% higher , the data points for @xmath150 are slightly shifted downwards and the agreement between theory and data becomes much better .    the damping rate , shown in fig .
[ fig : sloshing](b ) , increases strongly with the temperature : at @xmath153 the sloshing mode survives for several hundreds of oscillation periods , whereas at @xmath154 its amplitude decreases by a factor of @xmath155 after @xmath156 oscillations .
the in - medium modification of the cross - section has only a small effect on the damping of the sloshing mode . the experimental data from ref .
@xcite are in quite good agreement with our theoretical result , especially at the high temperature .
the strongly scattered damping rates at low temperature are probably an artefact due to a beat caused by the small residual ellipticity of the trap potential in the experiment .
we also compared our results with those shown in fig .
5(a ) of the recent paper @xcite by wu and zhang , where they solve the boltzmann equation numerically within the relaxation - time approximation with a local relaxation time @xmath157 taken from ref . @xcite . as we discussed in our previous work @xcite , the inclusion of higher - order moments accounts in an approximate way for the spatial dependence of the relaxation time . for the sake of comparison
, we use the same trap parameters as wu and zhang , namely , @xmath158 hz , @xmath159 hz , and @xmath160 @xmath97k , and neglect , as it is done in ref .
@xcite , the mean field and in - medium effects on the cross section .
because of the larger trap depth , the frequency shifts are somewhat weaker than those shown in fig . [
fig : sloshing](a ) .
since in the calculation by wu and zhang the frequency depends on the oscillation amplitude , we have to compare with the results obtained for the smallest amplitude [ red squares in fig .
5(a ) of ref . @xcite ] . if we extrapolate these results , shown only up to @xmath153 in ref .
@xcite , to higher temperatures , they are in quite good agreement with our results above @xmath161 . at lower temperatures , however , our frequency shift is slightly stronger than that by wu and zhang , probably because they use boltzmann instead of fermi distributions so that their density is more concentrated in the trap center at low temperatures .      in a harmonic trap ,
a characteristic feature of the kohn mode @xcite is that , as long as the interaction is translationally invariant , the cloud moves as a whole without any change in size or shape @xcite , as illustrated in figs .
[ fig : sloshskew](a )     direction ( @xmath162 being the density integrated over @xmath82 and @xmath80 ) for ( a ) the kohn mode , in which the cloud oscillates as a whole , and ( b ) the radial compressional dipole mode , where the cloud shape oscillates , while the center of mass of the cloud stays at rest .
all density profiles have the same normalization ( number of atoms ) and the same width . [
fig : sloshskew],title="fig:",width=151 ]   direction ( @xmath162 being the density integrated over @xmath82 and @xmath80 ) for ( a ) the kohn mode , in which the cloud oscillates as a whole , and ( b ) the radial compressional dipole mode , where the cloud shape oscillates , while the center of mass of the cloud stays at rest .
all density profiles have the same normalization ( number of atoms ) and the same width .
[ fig : sloshskew],title="fig:",width=151 ]    and [ fig : sloshbend](a ) .
, title="fig:",width=151 ] , title="fig:",width=151 ]    in an anharmonic trap , this is no longer the case and the cloud will be slightly distorted during the oscillation . in other words , because of the anharmonicity , the pure center - of - mass motion gets coupled to other collective modes .
these are damped as usual by collisions , which results in a damping of the sloshing mode .
what kind of distortions of the cloud shape can one expect ?
let us first note that , to first order in the perturbation , the cloud width can not be changed ( neither in the transverse nor in the axial direction ) , because the excitation operator is odd with respect to @xmath163 .    however , the density profile can become skewed in the direction of the oscillation ( i.e. , in the @xmath81 direction in our case ) . for illustration , fig . [
fig : sloshskew](b ) shows how the density profile changes during a radial compressional dipole oscillation . to derive the corresponding operator ,
let us start from the usual definition of the skewness , which is proportional to @xmath164 .
keeping only terms linear in the variations @xmath165 and @xmath166 , and using @xmath167 , we find that for small oscillations , this is equal to @xmath168 . if the trap is axially symmetric , collective modes of different multipolarity in the radial direction do not mix .
it is therefore useful to decompose the operator into dipolar and octupolar parts : @xmath169 with @xmath170 only @xmath171 can get a non - vanishing expectation value during the sloshing oscillation .
another possible distortion involves the axial ( @xmath80 ) direction .
let us consider again the sloshing mode in @xmath81 direction .
then one could imagine that the cloud near @xmath172 , owing to its larger radial size , is more sensitive to the anharmonicity of the trap and therefore oscillates slightly more slowly than the parts of the cloud at large @xmath173 .
hence , after a few oscillations , the oscillations near @xmath172 and at large @xmath173 become out of phase , which results in a bending of the cloud in the @xmath81-@xmath80 plane , as illustrated in fig .
[ fig : sloshbend](b ) . a suitable measure for such
a deformation is the expectation value of the operator @xmath174 like @xmath171 and @xmath175 , this operator is defined in such a way that it is insensitive to a translation of the cloud as a whole in @xmath81 direction .    within the first - order ansatz ( [ ansatzfirst ] ) ,
it is easy to see that in a harmonic trap without mean field and to linear order in the perturbation , the density oscillates as @xmath176 , with @xmath177 , i.e. , the cloud shape remains unchanged .
however , in the presence of a mean field @xmath11 , the cloud starts to deform during the oscillation .
for illustration , the measures discussed above for the cloud skewness , @xmath178 , and for the bending @xmath179 , obtained with the first - order ansatz in the case of a harmonic trap ( @xmath180 ) , are displayed in fig .
[ fig : skewbend](a ) by the dash - dotted and dotted lines .     and bending @xmath181 of the cloud during the sloshing motion excited by @xmath182 at @xmath153 . in order to have comparable orders of magnitude , we multiplied the bending strengths by a factor @xmath183 .
( a ) in a harmonic trap , the unphysical skewness and bending of the first - order approximation ( dash - dotted and dotted lines ) are strongly suppressed when the third order moments are included ( solid and dashed lines ) . since in a harmonic trap
there is no damping , we replaced @xmath58 by @xmath184 to generate this plot .
( b ) skewness and bending during the sloshing motion in the anharmonic trap within the third - order approximation .
[ fig : skewbend ] ]    the existence of these distortions is in contradiction to the kohn theorem and an unphysical consequence of the crude ansatz . at third order ,
owing to the additional degrees of freedom of the ansatz ( [ ansatzthird ] ) , the unphysical change of the cloud shape in the harmonic potential is strongly suppressed , as shown by the solid and dashed lines in fig .
[ fig : skewbend](a ) .    while unphysical in a harmonic trap , a change of the cloud shape during the sloshing motion is expected in an anharmonic trap . from the preceding discussion it is clear that at least the third - order ansatz is required to get a meaningful description of this effect . in fig .
[ fig : skewbend](b ) , third - order results for the expectation values of @xmath171 and @xmath185 during the sloshing motion are displayed .
we note that , compared to the third - order results for a harmonic trap , both responses are now significantly increased , so that the peaks correspond to a real physical effect . moreover , the observed shapes are typical of what one can obtain in a schematic model by coupling two damped modes ( dipole and bending ) to an undamped one ( sloshing )
. finally one can also see that the peaks are centered around the sloshing frequency , which is lower than in the harmonic case .
in the previous section we have discussed coupling of the sloshing mode to other collective modes due to the trap anharmonicity . in particular , the ansatz ( [ ansatzthird ] ) contains the necessary terms to describe the transverse compressional dipole and bending modes corresponding to the operators ( [ eq : qdip ] ) and ( [ eq : qbend ] ) . in this section
, we will study these modes in more detail .
however , one should keep in mind that , since the excitation operators are already third - order ones , the ansatz ( [ ansatzthird ] ) can only give the leading - order result for these modes .
let us start with the radial dipole oscillation excited by the operator ( [ eq : qdip ] ) .
in an ideal fermi gas in an harmonic trap , this operator would excite two modes with frequencies @xmath186 and @xmath187 .    in fig .
[ fig : resp_dip ] ,     for the radial dipole mode ( excitation @xmath188 ) in the anharmonic trap for various temperatures .
the parameters are the same as in fig .
[ fig : sloshing ] .
[ fig : resp_dip ] ]    the corresponding response function is shown for various temperatures .
we see that in the collisionless regime , i.e. , at high temperature , both modes exist but their frequencies are significantly lowered by the anharmonicity of the trap .
even at @xmath154 both modes are still very strongly damped ( damping rate @xmath189 ) . at lower temperature ,
the damping gets so strong that both modes disappear completely : the hydrodynamic regime is not reached at any temperature . by artificially increasing the collision cross section to reach the hydrodynamic regime at @xmath190 , we found that there would be a mode at @xmath191 ( @xmath192 in the case of a harmonic trap ) , but in order to reach this regime in an experiment one would have to use a trap with a much lower radial frequency ( i.e. , a more spherical or even pancake - shaped trap ) . at the frequency of the sloshing mode ,
the response exhibits a characteristic dip due to the coupling between the dipole mode and the almost undamped sloshing mode ( the sloshing is damped mainly because of this coupling ) .    in order to get a better understanding of the character of the two modes in the collisionless regime , we display their velocity fields in fig . [
fig : v_dip ] .
-@xmath82 plane , averaged over @xmath80 , corresponding to the dipole modes at @xmath193 ( a ) and @xmath194 ( b ) at @xmath195 ( collisionless regime ) .
the background color indicates the density per spin state integrated over @xmath80 in units of @xmath196 , with @xmath197 .
[ fig : v_dip],title="fig : " ] -@xmath82 plane , averaged over @xmath80 , corresponding to the dipole modes at @xmath193 ( a ) and @xmath194 ( b ) at @xmath195 ( collisionless regime ) .
the background color indicates the density per spin state integrated over @xmath80 in units of @xmath196 , with @xmath197 .
[ fig : v_dip],title="fig : " ]    although both modes are excited by the same operator , their velocity fields are quite different . the velocity field of the low - lying mode ( @xmath198 at @xmath195 ) , cf . fig .
[ fig : v_dip](a ) , confirms that this mode comprises both center - of - mass and skewness oscillations .
the velocity field of the high - lying mode ( @xmath199 at @xmath195 ) , shown in fig .
[ fig : v_dip](b ) , resembles that of the dipole compression modes in atomic nuclei , shown e.g. , in fig .
. @xcite .
although in an axially symmetric trap the octupole mode does not couple to the sloshing mode , let us for the sake of completeness briefly discuss this mode .
like the dipole operator ( [ eq : qdip ] ) , the octupole operator ( [ eq : qoct ] ) would excite two modes at @xmath200 and @xmath187 if the gas was collisionless and in a harmonic trap .
the octupole response function for the realistic case is displayed in fig .
[ fig : octupole](a ) .     for the radial octupole mode ( excitation @xmath201 ) in the anharmonic trap for various temperatures .
the parameters are the same as in fig .
[ fig : sloshing ] .
( b ) velocity field corresponding to the high - lying mode at @xmath202 ( see caption of fig . [
fig : v_dip ] for details ) .
[ fig : octupole],title="fig : " ]   for the radial octupole mode ( excitation @xmath201 ) in the anharmonic trap for various temperatures .
the parameters are the same as in fig .
[ fig : sloshing ] .
( b ) velocity field corresponding to the high - lying mode at @xmath202 ( see caption of fig . [
fig : v_dip ] for details ) .
[ fig : octupole],title="fig : " ]    as in the dipole response , both modes are considerably shifted downwards because of the trap anharmonicity .
the damping is even stronger than in the case of the dipole modes .
both octupole modes have a similar velocity field , as an example fig .
[ fig : octupole](b ) shows that of the higher - lying mode . in a nuclear physics context
it was found that their difference lies mainly in the non - diagonal pressure tensor ( quadrupole moments in momentum space ) @xcite . at @xmath190 , there is only a very broad peak around the hydrodynamic frequency @xmath203 predicted in ref .
@xcite for a harmonic trap .
however the strong damping shows that one is still far from the hydrodynamic regime . by increasing again artificially the cross - section ,
the hydrodynamic frequency is found to be lowered to @xmath204 due to anharmonicity .    finally , in fig .
[ fig : bending ]     for the bending mode ( excitation @xmath205 ) in the anharmonic trap for various temperatures .
the parameters are the same as in fig .
[ fig : sloshing ] .
( b ) velocity field in the @xmath81-@xmath80 plane , averaged over @xmath82 , corresponding to the bending mode at @xmath206 .
the background color indicates the density integrated over @xmath82 in units of @xmath207 with @xmath208 .
[ fig : bending],title="fig : " ]   for the bending mode ( excitation @xmath205 ) in the anharmonic trap for various temperatures .
the parameters are the same as in fig .
[ fig : sloshing ] .
( b ) velocity field in the @xmath81-@xmath80 plane , averaged over @xmath82 , corresponding to the bending mode at @xmath206 .
the background color indicates the density integrated over @xmath82 in units of @xmath207 with @xmath208 .
[ fig : bending],title="fig : " ]    we display the strength function and the velocity field of the bending mode . in a non - interacting gas in a harmonic trap , there would be three modes at @xmath209 and @xmath210 . in the realistic case
, these peaks can not be resolved and one sees a damped oscillation with a frequency close to @xmath89 .
in fact , the behavior of its frequency and damping rate is qualitatively similar to that of the sloshing mode , but the temperature dependence of the frequency is stronger . above @xmath211 , the damping is weaker than that of the sloshing mode , but this is probably an artefact of our approximation : we expect that the bending mode will receive an additional damping ( analogous to that of the sloshing mode ) from its coupling to a radial dipole mode modulated in @xmath80 [ operator @xmath212 , which is not included in the third - order ansatz . except at the lowest temperature ,
the mixing between sloshing and bending mode is weaker than that between sloshing and dipole mode , and the bending response does not show the characteristic dip that was found in the dipole response .
the velocity field , shown in fig . [ fig : bending](b ) confirms that this mode corresponds to a bending of the elongated cloud in the @xmath81-@xmath80 plane .
in this paper , we studied the frequency shift and damping of the sloshing mode of a normal - fluid fermi gas in a realistic anharmonic trap potential .
we used the moments method in order to find approximate solutions of the boltzmann equation with mean field and in - medium cross - section . already at first order ,
the moments method predicts a downwards shift of the frequency of the sloshing mode .
however , the first order ansatz is insufficient to describe the coupling between the center - of - mass and the internal degrees of freedom of the gas , which is responsible for the damping of the sloshing mode in an anharmonic trap .
we therefore extended the ansatz to include also third - order phase - space moments .
the third - order ansatz contains not only the sloshing mode , but also more complicated radial dipole and bending modes .
the important mechanism for the damping of the sloshing mode is its coupling to these modes ( especially to the dipole mode ) due to the trap anharmonicity .
in addition , the third - order ansatz allowed us to discuss the radial octupole mode .
the comparison with the available experimental data for the sloshing mode shows that the essential features of the frequency shift can be reproduced .
quantitatively , the calculated ratios @xmath150 are somewhat lower than those given in ref . @xcite .
since both the mean field and the third - order moments give only small corrections to the frequency shift predicted by the first - order moments method , further improvements of the theory will probably give even smaller corrections . an effect which has not been considered in our linear - response study
is the dependence of the sloshing frequency on the amplitude of the mode .
however , this effect goes into the wrong direction , since the restoring force gets weaker with increasing amplitude ( see fig .
[ fig : vtrap ] ) and therefore the frequency of the sloshing mode will be further reduced ( cf . also @xcite ) . as discussed in sec .
[ sec : extendedansatz ] , a possible explanation of the discrepancy is that the experimental trap frequency @xmath89 may have been slightly underestimated in the analysis of ref .
actually , our results could be used to determine the true value of @xmath89 .
the agreement between the calculated and measured damping rates of the sloshing mode is satisfactory and makes us confident that our approach contains the essential physics to describe the coupling between the sloshing and internal degrees of freedom .
this may also help to better understand other effects where the dissipation of center - of - mass kinetic energy plays a role , such as the heating of the cloud due to laser - beam - pointing noise discussed in ref .
@xcite .
in addition to the sloshing mode , we studied the radial dipole and octupole modes and the bending mode in the @xmath81-@xmath80 plane .
however , one should keep in mind that for these modes the third - order ansatz contains only the leading order , and higher - order corrections can be important .
we found that , with the present trap parameters , the dipole and octupole modes do not behave hydrodynamically at any temperature . the frequencies of the dipole , octupole , and bending modes are more strongly affected by the anharmonicity than the frequency of the sloshing mode .
for the dipole and octupole modes , this can be intuitively understood , since these higher - order modes are more sensitive to what happens at larger @xmath213 , where the atoms feel the anharmonicity of the trap .
the same argument should also apply to the radial quadrupole and scissors modes not studied in the present paper .
one should therefore be careful when comparing theoretical results obtained for a harmonic trap with experimental ones ( as done in @xcite ) , since the anharmonicity effects can not be completely absorbed in a renormalized radial trap frequency .    as a final remark we note that a similar third - order ansatz as the one used here [ eq .
( [ ansatzthird ] ) with @xmath214 can be used to describe one of the higher - order axial modes studied recently by the innsbruck group @xcite ( the @xmath215 mode in the notation of @xcite , corresponding to a motion as shown in fig .
[ fig : sloshskew](b ) but in @xmath80 instead of @xmath81 direction ) .
having more nodes than the usually considered collective modes , these kinds of modes are very interesting since they are closer to sound waves in uniform systems .
we thank l. sidorenkov and r. grimm for useful explanations concerning the experimental data of @xcite and for the permission to use the unpublished data of @xcite .
is supported by the _ fundao para a cincia e a tecnologia _ ( fct , portugal ) and the _
european social fund _ ( esf ) via the post - doctoral grant sfrh / bpd/64405/2009 . s.c . also acknowledges partial support by qren / feder , the compete programme , and fct under the project no .
ptdc / fis/113292/2009 .    * j. kinast , s.l .
hemmer , m.e .
gehm , a. turlapov , and j.e .
thomas , phys .
lett . * 92 * , 150402 ( 2004 ) .
m. bartenstein , a. altmeyer , s. riedl , s. jochim , c. chin , j. hecker denschlag , and r. grimm , phys .
92 , 203201 ( 2004 ) ; a. altmeyer , s. riedl , c. kohstall , m.j .
wright , j. hecker denschlag , and r. grimm , e - print cond - mat/0611285 ( 2006 ) .
j. kinast , a. turlapov , and j.e .
thomas , phys .
a * 70 * 051401(r ) ( 2004 ) ; j. kinast , a. turlapov , and j.e .
thomas , phys .
a * 71 * , 029901(e ) ( 2005 ) .
a. altmeyer , s. riedl , c. kohstall , m.j .
wright , r. geursen , m. bartenstein , c. chin , j. hecker denschlag , and r. grimm , phys .
lett . * 98 * , 040401 ( 2007 ) .
a. altmeyer , s. riedl , m.j .
wright , c. kohstall , j. hecker denschlag , and r. grimm , phys .
a * 76 * , 033610 ( 2007 ) .
wright , s. riedl , a. altmeyer , c. kohstall , e.r .
snchez guajardo , j. hecker denschlag , and r. grimm phys .
99 , 150403 ( 2007 ) .
s. riedl , e. r. snchez guajardo , c. kohstall , a. altmeyer , m. j. wright , j. hecker denschlag , and r. grimm , g. m. bruun and h. smith , phys .
a * 78 * , 053609 ( 2008 ) . w. kohn , phys
. rev . * 123 * , 1242 ( 1961 ) .
l. brey , n.f .
johnson , and b.i .
halperin , phys .
b * 40 * , 10647 ( 1989 ) .
dobson , phys .
. lett . * 73 * , 2244 ( 1994 ) .
r. grimm , m. weidemller , and y.b .
ovchinnikov , adv . at .
. phys . * 42 * , 95 ( 2000 ) .
l. wu and y. zhang , phys .
a * 85 * , 045601 ( 2012 ) .
l. sidorenkov and r. grimm , unpublished ; the present damping data have been obtained by reanalyzing the raw data of ref .
l. wu and y. zhang , europhys .
lett . * 97 * , 16003 ( 2012 ) .
s. chiacchiera , t. lepers , d. davesne , and m. urban , phys . rev .
a * 79 * , 033613 ( 2009 ) .
bruun and h. smith , phys .
a * 75 * , 043612 ( 2007 ) .
t. lepers , d. davesne , s. chiacchiera , and m. urban , phys . rev . a * 82 * , 023609 ( 2010 ) .
s. chiacchiera , t. lepers , d. davesne , and m. urban , phys . rev .
a * 84 * , 043634 ( 2011 ) .
e.m . lifshitz and l.p .
pitaevskii , _ physical kinetics _ , landau - lifshitz course of theoretical physics , vol
. 10 ( pergamon , oxford , 1980 ) .
da providncia , nucl .
a * 489 * , 111 ( 1988 ) .
a. altmeyer , phd thesis ( innsbruck , 2007 ) .
j. kinast , a. turlapov , j.e .
thomas , q. chen , j. stajic , and k. levin , science * 307 * , 1296 ( 2005 ) .
q. chen , j. stajic , and k. levin , phys .
rev . lett . * 95 * , 260405 ( 2005 ) .
p. massignan , g.m .
bruun , and h. smith , phys .
a * 71 * , 033607 ( 2005 ) .
serr , t.s .
dumitrescu , t. suzuki , and c.h .
dasso , nucl .
phys . a * 404 * , 359 ( 1983 ) .
h. kohl , p. schuck , and s. stringari , nucl .
a * 459 * , 265 ( 1986 ) .
a. griffin , w .- c .
wu , and s. stringari , phys .
* 78 * , 1838 ( 1997 ) .
savard , k.m .
ohara , and j.e .
thomas , phys .
rev . a * 56 * , r1095 ( 1997 ) .
snchez guajardo , phd thesis ( innsbruck , 2011 ) .


the thermodynamic analogue of the casimir effect originating from the confinement of vacuum fluctuations @xcite is the critical casimir effect due to the presence of long - ranged thermal fluctuations in a fluid close to its critical point at @xmath0 .
the corresponding critical casimir forces have been predicted theoretically by fisher and de gennes in 1978 @xcite .
experimentally , critical casimir forces have been studied only during the last decade when first measurements were performed indirectly via monitoring the thickness of wetting films upon approaching a critical end point @xcite .
later on , the critical casimir force has been measured directly by using colloidal particles suspended in a binary liquid mixture @xcite .
generically , the surfaces , which confine a binary liquid mixture , preferentially attract one of its two components leading to either positive [ @xmath1 or negative [ @xmath2 values of the scalar order parameter @xmath3 which describes the difference between the local concentration of one of the two components and its critical value .
this generic preference of the surfaces confining the liquid can be described by effective surface fields . upon approaching @xmath4 , the critical adsorption profiles , which describe the concentration enhancement near the surface , become long ranged due to the concomitant divergence of the bulk correlation length @xcite . in semi - infinite systems , the transition from the phase in which only the region near the single surface is ordered to the one
in which also the bulk is ordered is known as the extraordinary or normal transition @xcite . for two surfaces opposing each other , depending on the mutual combinations of boundary conditions ( bcs ) critical casimir forces are either attractive or repulsive .
their range is set by the bulk correlation length @xmath5 .
whereas in fluids @xmath5 is typically of molecular size , it attains values of the order of micrometers upon approaching the critical point @xcite . thus
, critical fluctuations may induce effective interaction potentials with a strength of several @xmath6 at the nanometer and micrometer scale .
moreover , critical casimir forces are universal in character : due to the divergence of the correlation length molecular details of the confined binary liquid mixture become irrelevant and only a few gross features of the system determine the main characteristics of the critical casimir forces @xcite .     from a lower , parallel , crenellated wall which exhibits a periodic pattern with period @xmath7 consisting of rectangular crenels of width @xmath8 and depth @xmath9 and merlons of width @xmath10 and height @xmath9 .
the system is spatially invariant in the @xmath11-direction . for later reference
, the box at the bottom summarizes the definitions of the various scaling variables which the scaling function of the critical casimir force depends on ; @xmath12 is the reduced temperature @xmath13 and @xmath14 is the nonuniversal amplitude of the bulk correlation length @xmath15 in the disordered phase , which defines the universal critical exponent @xmath16 . ,
title="fig:",width=302 ]   from a lower , parallel , crenellated wall which exhibits a periodic pattern with period @xmath7 consisting of rectangular crenels of width @xmath8 and depth @xmath9 and merlons of width @xmath10 and height @xmath9 .
the system is spatially invariant in the @xmath11-direction . for later reference
, the box at the bottom summarizes the definitions of the various scaling variables which the scaling function of the critical casimir force depends on ; @xmath12 is the reduced temperature @xmath13 and @xmath14 is the nonuniversal amplitude of the bulk correlation length @xmath15 in the disordered phase , which defines the universal critical exponent @xmath16 . , title="fig:",width=321 ]    [ fig : sketch ]    in view of nano- and micro - electromechanical devices , nowadays various experimental techniques are available to endow solid surfaces with precisely defined geometrical or chemical structures ( see , e.g. , refs .
critical casimir forces for _ chemically _ structured confinements have been studied theoretically @xcite as well as in experiments with colloidal particles @xcite .
it has been demonstrated that such patterns induce _ lateral _ critical casimir forces which can be used to trap particles reversibly along the lateral direction in a designated way @xcite .
moreover , a suitable combination of attractive and repulsive critical casimir forces may even lead to stable levitation @xcite .    here
, we study critical casimir forces for _ geometrically _ structured confinements .
the quantum - electrodynamic casimir effect in the presence of geometrically structured surfaces has been studied theoretically and experimentally for various surface topographies @xcite .
two opposing surface gratings are subject to lateral quantum electrodynamic casimir forces @xcite . in this context ,
experimental studies focus on spherical particles near crenellated surfaces , i.e. , forming grooves with rectangular cross - sections @xcite .
typically , the radius of the spherical particles is much larger than the period of the pattern of the crenellated surface such that , effectively , in the region of closest approach the system mimics the geometrical setup of a planar wall near a parallel , crenellated surface .
critical adsorption and critical casimir forces for geometrically structured confinements have been studied for structures shaped as wedges and ridges with triangular cross - section @xcite , as well as for truncated wedges @xcite .
it was found that , for large distances between a sawtooth - shaped wall and a planar wall , the critical casimir force effectively corresponds to the one between two planar walls ; on the other hand , for short distances between the two surfaces , the tips of the ridges dominate the order parameter profile and the characteristic power law behavior of the critical casimir force differs from that for planar surfaces @xcite .
the present study extends these previous investigations @xcite into various directions .
we consider a crenellated substrate close to a planar substrate at minimal surface - to - surface distance @xmath17 , as shown in fig .
[ fig : sketch ] .
the details of the geometry as well as the finite - size scaling of the critical casimir phenomena are described in section  [ sec : theory ] .
we calculate universal scaling functions for the critical casimir forces and for the order parameter profiles within mean - field theory and for identical chemical bcs at both walls ( @xmath18 configuration ) . in section  [ sec : profiles ] we first study the universal features of the order parameter profiles close to @xmath4 for the geometry under consideration .
second , in section  [ sec : force ] we study the critical casimir forces acting on such geometrically structured substrates .
finally , in section  [ sec : summary ] we summarize our main findings .
according to the theory of finite - size scaling , the singular contribution to the critical casimir force is described by a universal scaling function , which is independent of the molecular details of the binary liquid mixture and depends only on the bulk universality class of the associated critical point @xcite . here , we focus on the ising universality class characterized by a scalar order parameter @xmath3 , which encompasses the experimentally relevant binary liquid mixtures and simple fluids . upon approaching the critical point of the fluid ,
the bulk correlation length diverges as @xmath19 , where @xmath20 in spatial dimension @xmath21 and @xmath22 in spatial dimension @xmath23 @xcite ; @xmath24 are nonuniversal amplitudes characterized by the universal ratio @xmath25 in @xmath21 and @xmath26 in @xmath23 .
the sign of the reduced temperature @xmath27 is chosen such that @xmath28 corresponds to the mixed ( disordered ) phase of the fluid , whereas @xmath29 corresponds to the ordered phase , corresponding to spontaneous phase separation .
for an upper critical point the homogeneous phase is found at high temperatures , and one has @xmath13
. however , many experimentally relevant binary liquid mixtures exhibit a lower critical point ; in this case @xmath30 .    in general , the sign and the amplitude of the critical casimir force depend on the types of effective chemical bcs at the walls and on the geometry of the confining surfaces . here
we focus on the case of equal , symmetry - breaking @xmath31 bcs , which corresponds to the generic case of preferential adsorption of one of the two species of a binary liquid mixture .
this leads to an attractive critical casimir force @xcite .
inspired by experiments encompassing binary liquid mixtures of water and lutidine ( with a lower critical point at @xmath32 ) @xcite , here a binary liquid mixture with a lower critical point is considered at fixed pressure and at its critical composition .
first , we briefly review the film geometry . in this case
the liquid is confined between two parallel , macroscopically extended walls at a distance @xmath33 . according to renormalization group theory
the critical casimir force @xmath34 per area of one wall , which is acting on the parallel walls ( @xmath35 ) , scales as @xcite @xmath36 the scaling function @xmath37 depends only on a single scaling variable given by the film thickness @xmath33 in units of @xmath38 , raised to the power @xmath39 . for equal chemical bcs , as discussed here , @xmath40 is negative , so that the critical casimir force is attractive . for @xmath41
the scaling function of the critical casimir force reduces to a universal constant value , the so - called critical casimir amplitude @xcite : @xmath42 accordingly , at @xmath4 the critical casimir force decays algebraicly @xmath43 .
away from criticality , the critical casimir force decays exponentially as a function of @xmath44 .
the scaling function @xmath37 has been calculated exactly in @xmath45 @xcite , for @xmath46 using a perturbative field - theoretical method @xcite or a local - functional method @xcite , and in @xmath21 numerically via monte carlo simulations @xcite .      in the following we consider a crenellated wall located at a minimal surface - to - surface distance @xmath17 from a planar wall as sketched in fig .
[ fig : sketch ] .
the width and depth of the crenels are given by @xmath8 and @xmath9 , respectively , and the structure is periodic along the lateral @xmath47-direction with period @xmath7 , so that the width of the merlons , i.e. , the surface - to - surface separation between two neighboring crenels , is given by @xmath10 .
accordingly , the corresponding universal contribution to the critical casimir force @xmath48 per area of the planar wall scales as @xmath49 where the geometrical parameters form the following scaling variables : @xmath50    the critical casimir force between a crenellated and a planar wall attains the value of the corresponding force between two planar walls in various limits .
for very shallow crenels , i.e. , in the limit @xmath51 or @xmath52 , the following relation holds : @xmath53 as will be discussed in more detail below , we find that for narrow crenels with @xmath54 , the order parameter profile attains the form of the planar wall geometry , _ independent _ of the value of @xmath55 : @xmath56 similarly , for @xmath57 , one has @xmath58    for very broad crenels compared with the film thickness , i.e. , @xmath59 , the limiting behavior of the scaling function of the critical casimir force is given by the average of the scaling functions for two planar walls at distances @xmath17 and @xmath60 , respectively .
this corresponds to the so - called proximity force approximation ( pfa ) , which we discuss in the following . within the pfa
the surfaces are subdivided into infinitely small elements parallel to each other , and the resulting force is obtained by pairwise adding the individual contributions to the force .
the derjaguin approximation is a special case of the pfa for smoothly curved surfaces .
detailed comparisons of the derjaguin approximation with experimental and theoretical results revealed a good agreement for a large range of parameters @xcite . for the present geometry the resulting critical casimir force @xmath61 per area [ eq .
] acts along the @xmath62-direction and is the sum of two contributions : ( i ) the force between the fraction @xmath63 of the upper flat wall and the top surfaces of the merlons separated by the distance @xmath17 , and ( ii ) the force between the fraction @xmath64 of the upper flat wall and the bottom surfaces of the crenels separated by the distance @xmath60 : @xmath65 accordingly , the scaling function of the critical casimir force within the pfa is given by @xmath66 the scaling function @xmath67 is independent of @xmath68 because within pfa the effective interactions between the steps and the upper wall are ignored .
that is , as long the values of @xmath64 and @xmath55 are fixed , steps may be arbitrarily added to or removed from the structured substrate ( e.g. , via a transformation @xmath69 and @xmath70 , where @xmath71 ) without changing the force within pfa , independent of the value of @xmath68 .
analogous to the case of chemically striped surfaces , discussed in detail in refs .
@xcite , every isolated geometrical step gives rise to a contribution to the scaling function of the critical casimir force per area which is proportional to @xmath72 .
the asymptotic behavior for @xmath59 of the universal scaling function for the critical casimir force between a planar wall and a crenellated surface is therefore given by @xmath73 where we define @xmath74 as the universal contribution of a pair of geometrical steps as sketched in fig .
[ fig : sketch_step](a ) . for a fixed value @xmath75 , in eq .
@xmath76 actually does not depend on @xmath64 because in the limit @xmath59 the distance @xmath17 between the steps and the upper wall is much smaller than the step - step distances @xmath8 and @xmath10 , respectively , so that step - step interactions are negligible in this limit . in the following
we therefore define @xmath77 and consider this as the generic case . for very deep crenels
@xmath78 attains a @xmath79-dependent value which corresponds to the contribution to the critical casimir force between a pair of top corners of the right - angled edges of the merlons opposite to a planar wall ( see the sketch in fig .
[ fig : sketch_step](b ) ) .
( eq .  ) for @xmath80 at fixed values of @xmath81 and @xmath82 .
( a ) the surfaces of the lower , geometrically structured substrate which are parallel to the film ( indicated by a green line ) give rise to the pfa contribution @xmath67 to the universal scaling function @xmath83 of the critical casimir force . in the limit @xmath59 , the steps consisting of the upper and lower edges as well as the sidewalls ( indicated by the red lines and dots ) generate the step contribution @xmath84 to the critical casimir force ( see eq .  ) . for fixed @xmath75 , the step - step interaction vanishes and @xmath76 becomes independent of @xmath64 .
( b ) for @xmath85 the step contribution @xmath86 to the critical casimir force effectively corresponds to the contribution of two right - angled corners opposite to a planar substrate and attains a limiting value independent of @xmath55 .
, title="fig:",width=302 ]   ( eq .  ) for @xmath80 at fixed values of @xmath81 and @xmath82 . ( a ) the surfaces of the lower , geometrically structured substrate which are parallel to the film ( indicated by a green line ) give rise to the pfa contribution @xmath67 to the universal scaling function @xmath83 of the critical casimir force . in the limit @xmath59 ,
the steps consisting of the upper and lower edges as well as the sidewalls ( indicated by the red lines and dots ) generate the step contribution @xmath84 to the critical casimir force ( see eq .  ) . for fixed @xmath75 , the step - step interaction vanishes and @xmath76 becomes independent of @xmath64 .
( b ) for @xmath85
the step contribution @xmath86 to the critical casimir force effectively corresponds to the contribution of two right - angled corners opposite to a planar substrate and attains a limiting value independent of @xmath55 . , title="fig:",width=321 ]    [ fig : sketch_step ]      the standard landau - ginzburg - wilson fixed - point effective hamiltonian describing critical phenomena of the ising bulk universality class is given by @xcite @xmath87=\int_v\,{{\ensuremath{\textrm{d}}}}^d{r}\,\left\ {          \frac{1}{2}(\nabla\phi)^2          + \frac{{\hat{\tau}}}{2}\phi^2         + \frac{u}{4!}\phi^4         -h\phi               \right\},\ ] ] where @xmath88 is proportional to the order parameter describing the fluid , which completely fills the accessible volume @xmath89 in @xmath90-dimensional space .
the statistical weight of a configuration @xmath88 is proportional to @xmath91 .
the parameter @xmath92 in eq .   is proportional to @xmath12 , and @xmath93 is a coupling constant . the last term in eq .
vanishes for the case considered here ( i.e. , @xmath94 ) , which corresponds to the situation that the concentrations of the species forming the binary liquid mixture are fixed to their critical values . in a finite - size system
the bulk hamiltonian @xmath95 $ ] is supplemented by appropriate surface and curvature ( edge ) contributions @xcite .
this surface contribution , which adds to eq .  , is given by @xcite @xmath96=\int_{\partial v}\,{{\ensuremath{\textrm{d}}}}^{(d-1)}{r}\,\left\ {              \frac{c}{2}\phi^2              -h_1\phi               \right\},\ ] ] where @xmath97 is the so - called surface enhancement and @xmath98 is a surface field ; @xmath99 is the surface of the volume @xmath89 . in the strong adsorption limit @xcite ,
as discussed in the present study , these contributions generate boundary conditions for the order parameter such that @xmath100 corresponding to @xmath31 bcs .
thus , the use of eq .
together with additional surface contributions can be replaced by applying the appropriate bc to @xmath3 and by using eq .   throughout the bulk .
the mean - field order parameter profile minimizes the hamiltonian , i.e. , @xmath101/\updelta\phi|_{\phi=\langle\phi\rangle}=0 $ ] . in the _ bulk _ ,
the mean - field order parameter is spatially constant and attains the values @xmath102 for @xmath103 and @xmath104 for @xmath105 , where , besides @xmath14 , @xmath106 is the only additional independent nonuniversal amplitude appearing in the description of bulk critical phenomena ; @xmath107 is a standard critical exponent . within mean - field theory ( mft ) the following relations
hold : @xmath108 and @xmath109 .    for the film geometry
, the mft scaling function for the critical casimir force can be determined analytically @xcite .
one finds [ see eq .  ] for the case of the same strong adsorption at both surfaces the critical casimir amplitude @xmath110 ^ 4/u\simeq-283.61\times u^{-1}$ ] , where @xmath111 is the complete elliptic integral of the first kind @xcite .
renormalization group arguments tell that mft provides the correct universal properties of critical phenomena for spatial dimensions above the upper critical dimension , i.e. , @xmath112 , up to logarithmic corrections in @xmath113 @xcite .
moreover , mft provides the lowest - order contribution to universal quantities within an expansion in terms of @xmath114 .
( see eqs . and ) at @xmath115 ( i.e. , @xmath116 ) for various values of @xmath82 and @xmath117 with @xmath118 fixed ( see fig .  [
fig : sketch ] ) : ( a ) @xmath119 ( deep crenels and widely spaced merlons ) , ( b ) @xmath120 , @xmath121 ( shallow crenels and widely spaced merlons ) , and ( c ) and @xmath122 ( shallow crenels and closely spaced merlons ) .
the local values of @xmath123 are indicated by the color code given by the side bars .
large values of @xmath124 indicate strong adsorption of the fluid close to the surfaces . from the shape of the contour lines for @xmath125 near the structured wall ( lower white lines ) the dependence of the critical adsorption profile on the scaling variables @xmath55 and @xmath68 is clearly visible . whereas for widely spaced merlons with @xmath121 the adsorption profile follows the shape of the crenelated surface , for closely spaced merlons with @xmath126 the crenels are filled with the adsorbed fluid and for increasing normal distances @xmath127 the order parameter profile rapidly adopts an effective planar wall behavior corresponding to rather straight contour lines . , title="fig:",width=321 ]   ( see eqs . and ) at @xmath115 ( i.e. , @xmath116 ) for various values of @xmath82 and @xmath117 with @xmath118 fixed ( see fig .  [
fig : sketch ] ) : ( a ) @xmath119 ( deep crenels and widely spaced merlons ) , ( b ) @xmath120 , @xmath121 ( shallow crenels and widely spaced merlons ) , and ( c ) and @xmath122 ( shallow crenels and closely spaced merlons ) .
the local values of @xmath123 are indicated by the color code given by the side bars .
large values of @xmath124 indicate strong adsorption of the fluid close to the surfaces . from the shape of the contour lines for @xmath125 near the structured wall ( lower white lines )
the dependence of the critical adsorption profile on the scaling variables @xmath55 and @xmath68 is clearly visible .
whereas for widely spaced merlons with @xmath121 the adsorption profile follows the shape of the crenelated surface , for closely spaced merlons with @xmath126 the crenels are filled with the adsorbed fluid and for increasing normal distances @xmath127 the order parameter profile rapidly adopts an effective planar wall behavior corresponding to rather straight contour lines .
, title="fig:",width=321 ]   ( see eqs . and ) at @xmath115 ( i.e. , @xmath116 ) for various values of @xmath82 and @xmath117 with @xmath118 fixed ( see fig .  [
fig : sketch ] ) : ( a ) @xmath119 ( deep crenels and widely spaced merlons ) , ( b ) @xmath120 , @xmath121 ( shallow crenels and widely spaced merlons ) , and ( c ) and @xmath122 ( shallow crenels and closely spaced merlons ) .
the local values of @xmath123 are indicated by the color code given by the side bars .
large values of @xmath124 indicate strong adsorption of the fluid close to the surfaces . from the shape of the contour lines for @xmath125 near the structured wall ( lower white lines ) the dependence of the critical adsorption profile on the scaling variables @xmath55 and @xmath68 is clearly visible . whereas for widely spaced merlons with @xmath121 the adsorption profile follows the shape of the crenelated surface , for closely spaced merlons with @xmath126 the crenels are filled with the adsorbed fluid and for increasing normal distances @xmath127 the order parameter profile rapidly adopts an effective planar wall behavior corresponding to rather straight contour lines . ,
title="fig:",width=321 ]   ( see eqs . and ) at @xmath115 ( i.e. , @xmath116 ) for various values of @xmath82 and @xmath117 with @xmath118 fixed ( see fig .  [
fig : sketch ] ) : ( a ) @xmath119 ( deep crenels and widely spaced merlons ) , ( b ) @xmath120 , @xmath121 ( shallow crenels and widely spaced merlons ) , and ( c ) and @xmath122 ( shallow crenels and closely spaced merlons ) .
the local values of @xmath123 are indicated by the color code given by the side bars .
large values of @xmath124 indicate strong adsorption of the fluid close to the surfaces . from the shape of the contour lines for @xmath125 near the structured wall ( lower white lines ) the dependence of the critical adsorption profile on the scaling variables @xmath55 and @xmath68 is clearly visible .
whereas for widely spaced merlons with @xmath121 the adsorption profile follows the shape of the crenelated surface , for closely spaced merlons with @xmath126 the crenels are filled with the adsorbed fluid and for increasing normal distances @xmath127 the order parameter profile rapidly adopts an effective planar wall behavior corresponding to rather straight contour lines .
, title="fig:",width=321 ]   ( see eqs . and ) at @xmath115 ( i.e. , @xmath116 ) for various values of @xmath82 and @xmath117 with @xmath118 fixed ( see fig .  [
fig : sketch ] ) : ( a ) @xmath119 ( deep crenels and widely spaced merlons ) , ( b ) @xmath120 , @xmath121 ( shallow crenels and widely spaced merlons ) , and ( c ) and @xmath122 ( shallow crenels and closely spaced merlons ) .
the local values of @xmath123 are indicated by the color code given by the side bars .
large values of @xmath124 indicate strong adsorption of the fluid close to the surfaces . from the shape of the contour lines for @xmath125 near the structured wall ( lower white lines ) the dependence of the critical adsorption profile on the scaling variables @xmath55 and @xmath68 is clearly visible . whereas for widely spaced merlons with @xmath121 the adsorption profile follows the shape of the crenelated surface , for closely spaced merlons with @xmath126 the crenels are filled with the adsorbed fluid and for increasing normal distances @xmath127 the order parameter profile rapidly adopts an effective planar wall behavior corresponding to rather straight contour lines . ,
title="fig:",width=321 ]   ( see eqs . and ) at @xmath115 ( i.e. , @xmath116 ) for various values of @xmath82 and @xmath117 with @xmath118 fixed ( see fig .  [ fig : sketch ] ) : ( a ) @xmath119 ( deep crenels and widely spaced merlons ) , ( b ) @xmath120 , @xmath121 ( shallow crenels and widely spaced merlons ) , and ( c ) and @xmath122 ( shallow crenels and closely spaced merlons ) . the local values of @xmath123 are indicated by the color code given by the side bars .
large values of @xmath124 indicate strong adsorption of the fluid close to the surfaces . from the shape of the contour lines for @xmath125 near the structured wall ( lower white lines )
the dependence of the critical adsorption profile on the scaling variables @xmath55 and @xmath68 is clearly visible . whereas for widely spaced merlons with @xmath121
the adsorption profile follows the shape of the crenelated surface , for closely spaced merlons with @xmath126 the crenels are filled with the adsorbed fluid and for increasing normal distances @xmath127 the order parameter profile rapidly adopts an effective planar wall behavior corresponding to rather straight contour lines . ,
title="fig:",width=321 ]    [ fig : op ]    the order parameter @xmath3 exhibits the following scaling properties ( see subsec . 2.5 in ref .
@xcite ) : @xmath128 or equivalently @xmath129 with universal scaling functions @xmath130 and @xmath131 .
the bulk order parameter varies as @xmath132 .
this implies @xmath133 within mft one has @xmath134 so that with the definition @xmath135 eq .
renders the following mft approximation for the scaling function @xmath136 : @xmath137    in the following we present these mft results which we have obtained by minimizing numerically @xmath95 $ ] using a finite element method in order to obtain the ( spatially inhomogeneous ) profile @xmath138 for the geometries under consideration . here , we focus on the case of strong adsorption and the same chemical bcs at the two surfaces . for distances from the surface of a substrate which are small compared to @xmath38 , or for @xmath139 ,
the order parameter varies algebraicly .
in order to obtain a bc for the numerical calculations we use a short distance expansion ( see , e.g. , ref .
@xcite and reference therein ) with @xmath140 at the surfaces of the two walls shown in fig .
[ fig : sketch ] .    in fig .
[ fig : op ] the order parameter profile of a fluid confined between a planar and a crenellated wall at the bulk critical point @xmath115 is shown for @xmath141 ; @xmath138 depends on @xmath47 and @xmath62 and is invariant along the @xmath11-direction . from fig .
[ fig : op ] we can infer that for @xmath142 the order parameter profile follows the shape of the crenellated surface [ figs .
[ fig : op ] ( a ) and ( b ) ] , whereas for @xmath126 the space between the merlons corresponds to high values of the order parameter , i.e. , the crenels are `` filled '' with the adsorbed fluid [ fig .  [
fig : op ] ( c ) ] until an almost straight contour line has formed separating the fluid with high order parameter , which fills the crenels , from the fluid with lower order parameter in the middle of the slit .
hence for narrow crenels , for increasing @xmath127 the order parameter profiles rapidly approach the ones of a corresponding film of thickness @xmath17 .
the critical casimir forces are calculated directly from the numerically obtained mean - field order parameter profiles using the stress tensor @xcite . as in sec .  [
sec : profiles ] , here throughout we focus on the case @xmath143 .
we estimate the numerical error of the present method to be less than 1% .     of the critical casimir force as a function of @xmath144 for a fixed ratio @xmath143 and various values of @xmath82 and @xmath117 .
in addition to the numerically obtained full mft scaling functions ( full curves ) , the corresponding scaling functions @xmath67 obtained within the pfa [ eq .  ] ( dashed curves ) and the scaling function for two planar walls ( pw ) at distance @xmath17 ( dashed - dotted curves ) are shown . for @xmath119 , in ( a ) the mft scaling function
does neither agree with the one obtained within pfa nor with the one for the pw case , but lies roughly in between these two curves .
however , for shallow crenels and widely spaced merlons with @xmath120 and @xmath121 , in ( b ) the pfa is closer to the full mft scaling function . for shallow crenels and
closely spaced merlons with @xmath122 , in ( c ) the full mft scaling function is similar to the one for two planar walls ( pw ) .
figures [ fig : op ] and [ fig : force ] allow a direct comparison between the order parameter distribution and the critical casimir force at @xmath4 , i.e. , for @xmath145 . here
the scaling function @xmath83 is normalized with @xmath146 so that the ratio does not depend on the coupling constant @xmath147 , which within mft is undetermined .
, title="fig:",width=302 ]   of the critical casimir force as a function of @xmath144 for a fixed ratio @xmath143 and various values of @xmath82 and @xmath117 .
in addition to the numerically obtained full mft scaling functions ( full curves ) , the corresponding scaling functions @xmath67 obtained within the pfa [ eq .  ] ( dashed curves ) and the scaling function for two planar walls ( pw ) at distance @xmath17 ( dashed - dotted curves ) are shown . for @xmath119 , in ( a ) the mft scaling function
does neither agree with the one obtained within pfa nor with the one for the pw case , but lies roughly in between these two curves .
however , for shallow crenels and widely spaced merlons with @xmath120 and @xmath121 , in ( b ) the pfa is closer to the full mft scaling function . for shallow crenels and
closely spaced merlons with @xmath122 , in ( c ) the full mft scaling function is similar to the one for two planar walls ( pw ) .
figures [ fig : op ] and [ fig : force ] allow a direct comparison between the order parameter distribution and the critical casimir force at @xmath4 , i.e. , for @xmath145 . here
the scaling function @xmath83 is normalized with @xmath146 so that the ratio does not depend on the coupling constant @xmath147 , which within mft is undetermined .
, title="fig:",width=321 ]    [ fig : force ]    in fig .
[ fig : force ] the scaling function @xmath148 of the critical casimir force is shown as a function of @xmath144 for various values of @xmath82 and @xmath117 .
in addition to the numerically obtained full mft scaling functions ( solid curves ) , the corresponding scaling functions @xmath149 obtained within the pfa [ eq .  ] ( dashed curves ) , and the scaling function @xmath150 for two planar walls
( pw ) at distance @xmath17 ( dashed - dotted curves ) are shown . for @xmath119 , in fig .
[ fig : force](a ) the mft scaling function does neither agree with the one obtained within the pfa nor with the one for the pw case , but lies roughly in between these two curves . the corresponding mft order parameter profile at @xmath4 is shown in fig .
[ fig : op](a ) . for @xmath120 and @xmath121 , in fig .
[ fig : force](b ) the scaling function obtained within the pfa is closer to the full mft scaling function which corresponds to the critical mft order parameter profile shown in fig .
[ fig : op](b ) . for these shallow crenels and widely spaced merlons the relative contribution from the the right - angled steps of the merlons to the total critical casimir force is smaller than for the case shown in fig .  [
fig : force](a ) . for shallow crenels and closely spaced merlons with @xmath122 , in fig .
[ fig : force](c ) the mft scaling function is similar to the one for two planar walls ( pw ) .
this corresponds to the case in which the crenels are filled by a fluid with a high value of the order parameter ( see fig .  [
fig : op](a ) for @xmath115 ) .    as can be inferred from fig .
[ fig : force ] , in general the pfa deviates from the full mft results because critical phenomena do not allow for linear superposition . in order to study the deviations of the results for the force from the corresponding ones following from the assumption of pairwise additivity
, we study the ratio @xmath151 of the scaling function @xmath83 obtained within full mft and the one ( @xmath152 ) obtained within pfa . from eqs .
( [ eq : limit - omega-0 ] ) and ( [ eq : derjaguin - scaling ] ) we find for @xmath116 , i.e. , @xmath115 : @xmath153     of the scaling function @xmath83 of the critical casimir force obtained within full mft and the one ( @xmath152 ) obtained within pfa , for @xmath116 ( i.e. , @xmath115 ) and @xmath154 ( eq .  ) .
the ratio varies between @xmath155 and @xmath156 corresponding to deviations of the actual force from the force calculated within pfa between 0   in ( a ) @xmath151 is shown as a function of @xmath117 for several values of @xmath82 .
for widely spaced merlons ( @xmath157 ) , according to eq .
the ratio @xmath151 approaches @xmath155 .
the corresponding decay @xmath158 is shown as dashed lines .
for very shallow crenels ( @xmath52 ) , the ratio @xmath151 approaches @xmath155 , so that pfa is valid . in ( b ) @xmath151
is shown as a function of @xmath82 for various values of @xmath68 .
we find that , for deep crenels ( @xmath85 ) , @xmath151 attains a plateau . for large values of @xmath68
this corresponds to the limiting behavior given in eq .  ; ultimately , for @xmath59 pfa is valid . on the other hand for closely spaced merlons ( @xmath159 )
the critical casimir force approaches its value for an effective film geometry with the two parallel walls at distance @xmath17 . in this limit
the ratio @xmath151 is given by eq .
with @xmath154 and @xmath23 ( dashed line , `` pw '' ) . , title="fig:",width=302 ]   of the scaling function @xmath83 of the critical casimir force obtained within full mft and the one ( @xmath152 ) obtained within pfa , for @xmath116 ( i.e. , @xmath115 ) and @xmath154 ( eq .  ) . the ratio varies between @xmath155 and @xmath156 corresponding to deviations of the actual force from the force calculated within pfa between 0   in ( a ) @xmath151 is shown as a function of @xmath117 for several values of @xmath82 . for widely spaced merlons ( @xmath157 ) , according to eq .
the ratio @xmath151 approaches @xmath155 .
the corresponding decay @xmath158 is shown as dashed lines .
for very shallow crenels ( @xmath52 ) , the ratio @xmath151 approaches @xmath155 , so that pfa is valid . in ( b ) @xmath151
is shown as a function of @xmath82 for various values of @xmath68 .
we find that , for deep crenels ( @xmath85 ) , @xmath151 attains a plateau . for large values of @xmath68
this corresponds to the limiting behavior given in eq .  ; ultimately , for @xmath59 pfa is valid . on the other hand for closely spaced merlons ( @xmath159 )
the critical casimir force approaches its value for an effective film geometry with the two parallel walls at distance @xmath17 . in this limit
the ratio @xmath151 is given by eq .
with @xmath154 and @xmath23 ( dashed line , `` pw '' ) . , title="fig:",width=321 ]    [ fig : compare ]    figure  [ fig : compare ] shows the ratio @xmath151 as function of @xmath55 and @xmath68 for fixed values @xmath116 and @xmath154 ( eq .  ) .
the ratio @xmath151 varies between @xmath155 and @xmath156 which corresponds to deviations of the actual force from the one calculated within pfa between 0% and 100% , respectively . in fig .
[ fig : compare](a ) @xmath151 is shown as a function of @xmath68 for various values of @xmath55 .
this graph tells that the limiting behavior for @xmath59 given in eq .   already holds for @xmath160 .
the amplitude function @xmath161 [ eq .  ] has been determined via a least square fit to the numerical data . as expected on physical grounds , for very shallow crenels
the critical casimir force can be approximated by the corresponding pfa expression and hence @xmath162 for @xmath52 . in fig .
[ fig : compare](b ) we show @xmath151 as a function of @xmath55 for various values of @xmath68 . for deep crenels with @xmath163 this ratio reaches a plateau , i.e. , for @xmath85 the strength of the critical casimir force is not affected by the depth of the crenels so that for @xmath59 pfa becomes valid . for @xmath159
the critical casimir force approaches its value for an effective film geometry with the two parallel walls ( dashed line , `` pw '' ) at distance @xmath17 ( see eq .   for @xmath154 and @xmath23 ) .
this means that for closely spaced merlons the critical casimir force reduces to that between two parallel flat surfaces at separation @xmath17 .
[ see eq .
] of the step contribution to the critical casimir force for @xmath116 as a function of the reduced depth @xmath82 of the crenels . for @xmath164 a plateau value @xmath165
is reached and the step contribution to the critical casimir force becomes independent of the depth of the crenels . in ( b ) @xmath166 is shown as a function of @xmath79 as obtained for @xmath167 .
according to our analysis , for @xmath167 the limiting behavior for @xmath168 is almost reached ( see ( a ) ) .
( according to eq .
@xmath76 does not depend on @xmath64 as long as the value of the latter is not @xmath169 or @xmath155 . ) , title="fig:",width=302 ]   [ see eq .
] of the step contribution to the critical casimir force for @xmath116 as a function of the reduced depth @xmath82 of the crenels . for @xmath164 a plateau value @xmath165
is reached and the step contribution to the critical casimir force becomes independent of the depth of the crenels . in ( b ) @xmath166 is shown as a function of @xmath79 as obtained for @xmath167 . according to our analysis , for @xmath167 the limiting behavior for @xmath168
is almost reached ( see ( a ) ) .
( according to eq .
@xmath76 does not depend on @xmath64 as long as the value of the latter is not @xmath169 or @xmath155 . ) , title="fig:",width=321 ]    [ fig : a ]    as can be inferred from eq .   and fig .
[ fig : compare](a ) for large values of @xmath68 the critical casimir force acting between a flat and a crenellated wall can be described reliably by the sum of the pfa - contribution and a contribution @xmath170 .
the latter contribution stems from effects due to the presence of geometrical steps ( see fig .  [
fig : sketch_step ] ) . in fig .
[ fig : a](a ) we show @xmath171 as a function of @xmath55 as obtained numerically within full mft via least square fits to the data for @xmath116 . as expected , for @xmath52
, i.e. , @xmath172 , @xmath171 vanishes for the case of chemically homogeneous boundary conditions as considered here .
( for chemically inhomogeneous boundary conditions in lateral direction a line contribution arises due to the presence of a _ chemical _ step even for a planar substrate @xcite . ) for @xmath85 , i.e. , deep crenels , @xmath171 attains negative values and depends on @xmath79 only .
@xmath86 corresponds to the contribution to the critical casimir force of two infinitely extended right - angled corners opposite to a planar substrate ( see fig .  [
fig : sketch_step](b ) ) .
we find that this plateau value is reached for @xmath164 . in this limit the lower parts of the crenels including the right - angled wedges at the bottom do not contribute to the critical casimir force .
figure  [ fig : a](b ) shows the amplitude @xmath171 of the step contribution as a function of @xmath79 for deep crenels with @xmath167 , normalized by the absolute value @xmath146 of the critical casimir amplitude . for the chosen value @xmath167
we find that @xmath173 holds [ see eq .  ] .
we find that the functional shapes of this amplitude of the step contribution and of @xmath174 are similar but the minima of the two scaling functions are displaced such that the minimum of @xmath76 as a function of @xmath79 is located closer to the critical point at @xmath116 .
for comparison we note that the functional shape of the amplitude of the step contribution is very close to the numerically obtained mft scaling function of the attractive critical casimir force acting between a thin cylinder and a planar substrate , where the symmetry axis of the cylinder is parallel to the substrate @xcite . indeed , the lower white contour line in fig .
[ fig : op](a ) indicates that one may approximate each cross - section of the two right - angled edges of the merlons by a small inscribed quarter circle , keeping in mind that in the limit @xmath175 the critical adsorption profile near a thin cylinder is independent of the radius @xmath136 of the cylinder @xcite .
within mean field theory ( mft ) we have calculated the critical casimir force between a planar and a crenellated surface , i.e. , a periodic structure of geometric steps forming crenels and merlons ( fig .
[ fig : sketch ] ) . to this end
we have numerically calculated the critical order parameter profiles of a fluid confined by such a structure .
we have focused on the experimentally relevant case of binary liquid mixtures or simple liquids , which belong to the ising bulk universality class , and on @xmath31 boundary conditions ( bcs ) which correspond to the generic case of strong critical adsorption . using the stress tensor
we have calculated the universal scaling function of the singular contribution to the critical casimir force , acting on the confining walls along the normal direction upon approaching the critical point .
we have obtained the following main results :    1 .   according to finite size scaling ,
the universal scaling function of the critical casimir force depends on only a few scaling variables which describe the geometry under consideration [ eqs .  , , and
[ fig : sketch ] ] . in the limits , that the width @xmath8 or the depth @xmath9 of the crenels are small compared to the film width or the period of the geometric structure , the force reduces to the one acting between two parallel flat substrates at a distance @xmath17 [ eqs .    ] . on the other hand , for fixed reduced crenel depth @xmath55 and for reduced crenel width @xmath176 ,
the expression for the critical casimir force approaches the one obtained via the proximity force approximation ( pfa ) [ eq .
] as @xmath177 , which can be interpreted as a step contribution [ see eq .   and fig .
[ fig : sketch_step ] ] .
2 .   in sec .
[ sec : profiles ] we have shown generic examples of order parameter profiles ( i.e. , the universal contribution to the profiles of the deviations of the local concentration or density from its critical value in the bulk ) as obtained within mft at the critical point .
whereas for crenel sizes comparable with the film thickness the order parameter profile follows the geometrical structure [ figs .  [
fig : op](a ) and ( b ) ] , for shallow crenels and closely spaced merlons the order parameter profile resembles the one in between two planar walls [ fig .
[ fig : op](c ) ] . in this latter case
strong critical adsorption effectively suppresses the influence of the merlons .
these properties of the order parameter profiles agree with our results for the universal scaling function @xmath83 [ eq .  ] of the critical casimir force presented in section  [ sec : force ] : for deep crenels and widely spaced merlons the full mft scaling function of the critical casimir force as a function of the scaling variable @xmath79 does neither agree with the planar wall limit ( pw ) nor with the pfa [ fig .
[ fig : force](a ) ] . on the other hand ,
for shallow crenels and widely spaced merlons the pfa describes the actual full mft behavior better [ fig .
[ fig : force](b ) ] , whereas for shallow crenels and closely spaced merlons the pw limit is approached [ fig .
[ fig : force](c ) ] . thus , depending on the values of the scaling variables @xmath82 and @xmath117 , the critical casimir force interpolates between its limiting behaviors described by the pfa , which always underestimates the strength of the force , and the pw limit , which overestimates the critical casimir force
this can also be seen in fig .
[ fig : compare](b ) which shows the ratio of the universal scaling function of the critical casimir force and its corresponding pfa value as a function of @xmath55 .
4 .   upon increasing @xmath68 the scaling function @xmath67 as obtained within pfa
is approached as @xmath177 [ fig .
[ fig : compare]a ) ] , which can be interpreted as a step contribution with @xmath77 de facto independent of @xmath64 [ eq .  ] .
thus , for widely spaced merlons the critical casimir force can be described as the linear superposition of the planar wall contributions at distances @xmath17 and @xmath60 , respectively , plus additional step contributions stemming from the geometrical steps . for large values of @xmath55
this step contribution approaches a specific negative value , independent of the crenel depth @xmath9 [ fig .
[ fig : a](a ) ] .    to conclude , our numerical calculations within mft extend previous investigations of the critical casimir force caused by the confinement due to structured substrates .
in particular , we have analyzed in detail the crossover of the universal scaling function of the critical casimir force from an effective planar wall limit for small roughness to the pfa limit .
this may not only be useful for the analysis of critical casimir forces in such experimentally relevant , designed geometries , but could also help in understanding unavoidable roughness corrections in the case of planar geometries .
moreover , our results show that the first - order correction to the pfa can be expressed in terms of a contribution stemming from the individual geometrical steps .
this is analogous to similar situations involving chemically structured surfaces @xcite , which are discussed in ref .
@xcite . for comparison
we note that such step contributions arising from crenellated surfaces are also of importance for non - critical fluids consisting of rod - like particles close to the isotropic to nematic phase transition of the bulk fluid @xcite .
the authors thank nikolas brunner for very helpful contributions .    77 ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]
 + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty @noop * * , ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop _ _  ( , ) @noop _ _  ( ,  ) @noop * * ,   ( ) link:\doibase 10.1103/physrevlett.83.1187 [ * * , ( ) ] @noop * * ,   ( ) link:\doibase 10.1103/physrevlett.94.135702 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.97.075301 [ * * ,   ( ) ] link:\doibase 10.1016/j.physa.2007.07.072 [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1038/nature06443 [ * * ,   ( ) ] @noop * * ,   ( ) link:\doibase 10.1209/0295 - 5075/88/26001 [ * * ,   ( ) ] in  @noop _ _ ,  vol .  ,  ( ,  )  p.   in  @noop _ _ ,  vol .  ,
( ,  ) p.   @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) link:\doibase    10.1209/0295 - 5075/88/40004 [ * * ,   ( ) ] link:\doibase    10.1063/1.3464770 [ * * ,   ( ) ] link:\doibase 10.1039/c0sm00635a [ * * , ( ) ] http://stacks.iop.org/1742-5468/2010/i=11/a=p11003 [ ( ) ] @noop * * ,   ( ) @noop ( ) link:\doibase 10.1039/c3sm52858h [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.101.208301 [ * * ,   ( ) ] http://www.informaworld.com/10.1080/00268976.2011.553639 [ * * ,   ( ) ] @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ( ) link:\doibase 10.1103/physrevlett.101.160403 [ * * ,   ( ) ] http://dx.doi.org/10.1038/nphoton.2011.39 [ * * , ( ) ] link:\doibase 10.1103/physreva.86.062502 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.88.101801 [ * * , ( ) ] @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ( ) @noop * * ( ) link:\doibase 10.1103/physrevb.81.115417 [ * * ,   ( ) ] @noop * * ,   ( ) ; _ ibid . _ * 107 * , 019901(e ) ( 2011 ) link:\doibase 10.1103/physrevlett.105.250402 [ * * ,   ( ) ] link:\doibase    10.1103/physreva.87.052514 [ * * ,   ( ) ] http://dx.doi.org/10.1038/ncomms3515 [ * * ,   ( ) ] @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) link:\doibase 10.1063/1.2977999 [ * * ,   ( ) ] , preprint arxiv:1406.0962 ( 2014 ) in  @noop _ _ ,  vol .  ,  ( ,  )  p.   in  @noop _ _ ,  ( ,  )  p.
link:\doibase 10.1016/s0370 - 1573(02)00219 - 3 [ * * ,   ( ) ] @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) link:\doibase 10.1103/physrevb.49.8842 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.101.125702 [ * * ,   ( ) ] link:\doibase 10.1209/0295 - 5075/80/60009 [ * * ,   ( ) ] link:\doibase 10.1103/physreve.79.041142 [ * * ,   ( ) ] ; _ ibid . _
* 80 * , 039902(e ) ( 2009 ) link:\doibase 10.1103/physrevb.82.104425 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.82.174434 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.82.174433 [ * * ,   ( ) ] @noop * * ,   ( ) @noop * * ,   ( ) link:\doibase 10.1103/physrevb.47.5841 [ * * ,   ( ) ] ; _ ibid . _ 6740 ( 1993 ) @noop _ _ ,  ed .
( ,  ) link:\doibase 10.1063/1.3259188 [ * * ,   ( ) ] link:\doibase 10.1103/physreve.59.5081 [ * * ,   ( ) ] @noop * * ,   ( )