in a remarkable series of experiments in the nineties , sophisticated cooling techniques were used to expose the quantum nature of ultracold dilute gases . for bosons , the realization of bose - einstein condensation ( bec ) in 1995 @xcite revealed the importance of the underlying aggregating statistics for the existence of a broken symmetry phase in weakly interacting systems . for fermions ,
the statistics has just the opposite effect . due to the pauli exclusion principle ,
low energy collisions of fermionic particles are strongly suppressed , worsening the thermalization between the atoms and , therefore , making the cooling techniques more involved .
nevertheless , a few years later another milestone experiment was realized in which a quantum degenerate fermi gas was observed @xcite . since then , research with both types of degenerate quantum gases has experienced a strong growth , and different trends have been pursued , each showing interesting and relevant connections to other fields of physics such as atomic , molecular and condensed matter physics .
rotating quantum gases constitute one of these fast growing areas . in the case of bosons in a harmonic trap ,
intensive investigations have been carried out , both theoretically @xcite and experimentally @xcite . if the rotation frequency @xmath1 is of the order of ( but smaller than ) the radial trapping frequency @xmath2 , i.e. @xmath3 , the rotation is called slow . in this regime , intrinsic properties of rotating systems like the temperature dependence of the moment of inertia @xcite , rotation induced excitations @xcite , and the formation of vortices @xcite have been considered .
if one rotates faster , vortex arrays are created @xcite where small - amplitude oscillations , so called tkachenko modes , have been observed @xcite and theoretically investigated @xcite .
in addition , fast rotating ( @xmath4 ) becs have been linked to fascinating many - body phenomena such as the integer quantum hall effect @xcite or its fractional version @xcite , and are suspected to display quantum phase transitions to non superfluid strongly correlated states @xcite .    in the fast rotation regime ,
experiments in harmonic traps are limited by the loss of confinement due to centrifugal forces when the rotation frequency approaches the axial harmonic trapping frequency . to overcome this problem ,
a. l. fetter suggested adding a quartic term to the trap potential @xcite .
the successful experimental implementation of such a quartic term was carried out in the group of j. dalibard in paris @xcite . in their
set - up a laser is used to stir the quartically trapped condensates , making possible the exploration of many properties of the system such as , for instance , the vortex distribution , the optical density @xcite , and the shape oscillations @xcite .
the experimental realization triggered further theoretical research focusing on bose - einstein condensates ( bec ) .
for instance , their collective oscillations @xcite , thermodynamical @xcite and dynamical @xcite properties have been considered in the presence of the quartic term .
thus , the question arises , what new features come about through the interplay between rotation and fermionic statistics . in a harmonic trap ,
some properties have already been investigated .
for instance , the moment of inertia and its relation to quadrupole oscillations , in a framework valid for bosons and fermions @xcite , and the collective excitations @xcite have been studied . in the superfluid phase ,
formation of vortices @xcite has been observed and a rotation - induced phase separation has been proposed @xcite . in addition , the influence of the landau levels on the density profiles of two - dimensional polarized fermions in the context of the integer quantum hall effect @xcite and the experimental feasibility of fractional quantum hall states @xcite have also been theoretically considered .    in the case of an anharmonic trap , apart from studying the breathing mode in the bardeen - cooper - schrieffer ( bcs)-bec crossover regime @xcite , fermions have not been considered and most of the interesting aspects of the physics of rotating systems remain restricted to pure harmonic confinement . at low temperatures , polarized atomic fermi gases are essentially non - interacting due to pauli blocking and the fact that p - wave interactions are negligible in such systems .
however , due to the fermi pressure , particles tend to avoid the center of the trap , so that , effectively , a statistical interaction arises even in the case of an ideal gas , which bears some analogies with repulsively interacting bosons @xcite .
in this work we are concerned with rotating ideal fermi gases of spin polarized particles rotating with frequency @xmath1 .
we consider the same trap configuration as in the paris experiment , which will be described in detail in section [ sec:2 ] .
in section [ sec:3 ] we outline a general approach to formally justify the semiclassical approximation for calculating the density of states in anharmonic traps .
furthermore , we apply this technique to the paris trap and see that a treatment within the semiclassical approximation is justified . in section [ sec:4 ] we consider ground - state properties of such a system in a similar semiclassical approach .
we start by evaluating the fermi energy for arbitrary rotation , then we obtain the particle density , the momentum distribution , and free expansion time - of - flight absorption pictures . in section [ sec:5 ] temperature effects
are considered and we derive expressions for the grand - canonical thermodynamic properties in terms of incomplete fermi functions .
in addition we analyze analytic low- and high - temperature expansions . for the latter ,
a virial theorem approach is also presented .
section [ sec:6 ] discusses the main results of this work emphasizing the effects of rotation and anharmonicity in comparison to a non - rotating harmonic trap and presents an outlook for the field .
in the following we consider a trap potential which is anharmonic and consists of harmonic axially symmetric magnetic trapping supplemented by an additional gaussian trapping term in the @xmath5-plane .
the gaussian laser beam is approximated by its fourth order taylor series @xmath6 where @xmath7 and @xmath8 are the laser intensity and width , respectively .
the radial distance to the origin is denoted by @xmath9 .
the first term in is just a shift in the energy scale and will be neglected in the following .
the second order term will be absorbed into the harmonic trapping term .
the last term is the one responsible for the quartic confinement and is characterized by the anharmonicity constant @xmath10 .
it is important to notice that the harmonic trapping in the @xmath11-direction is not affected by the gaussian laser potential and is characterized by its harmonic oscillator frequency @xmath12 .
an issue in the field of rotating quantum gases is that the rotation frequency , with which the condensate effectively rotates , might be different from the one of the stirring laser , due to changes in the shape dependent rotation properties of the gas such as the moment of inertia @xcite .
we neglect these effects and assume that a stirring laser can adiabatically add angular momentum to the particles , so that the gas rotates with frequency @xmath1 , yielding an additional centrifugal harmonic term to the trap potential in the co - rotating frame @xcite . under these considerations ,
the rotating anharmonic potential takes the form @xcite @xmath13 where @xmath14 and @xmath15 are the axial harmonic energy and length scales , respectively , @xmath16 denotes the anisotropy of the trap , @xmath17 represents the anharmonicity of the trap , and @xmath18 is the mass of the atomic species under consideration .
the rotation is taken into account with the parameter @xmath19 . adopting the values given for the paris experiment @xcite , the resulting rotating trap
is then described by the radial harmonic trapping frequency @xmath20 , the axial harmonic trapping frequency @xmath21 , the force constant @xmath22 of the quartic trapping term , the anisotropy @xmath23 , and a tunable rotation frequency @xmath1 induced by the stirring laser .    in the paris experiment ,
a bec was obtained with the bosonic isotope @xmath24 , for which the dimensionless anharmonicity takes the value @xmath25 . since fermionic species like @xmath26 and @xmath27
have been sympathetic cooled by using the bosonic atoms @xmath24 @xcite and @xmath28 @xcite , respectively , we consider them likely candidates for experiments with ultracold degenerate fermi gases in the paris trap . for definiteness
, we evaluate all results of this paper for the parameters of the trap of the paris experiment and the @xmath26 species , for which we have @xmath29 ( see table [ tab:1 ] ) .
let us discuss the trap potential and some of its peculiarities concerning rotation .
for the moment , we restrict ourselves to the @xmath5-plane . as depicted in fig .
[ fig:1 ] , for slow rotations , i.e. @xmath30 , the minimum of the potential is located at the origin .
rotating faster one eventually reaches the point where the harmonic confinement vanishes ( @xmath31 ) and the trap is purely quartic .
hereafter this will be called critical rotation .
if the rotation is made even faster , one achieves the supercritical regime , where the harmonic confinement becomes centrifugal and the trap minimum is moved to @xmath32 and , with respect to the origin , has a depth given by @xmath33 which naturally characterizes the anharmonic trap .
the lowest available potential energy is , thus , @xmath34 , in the supercritical case .
, for the parameters given at table [ tab:1 ] at both sub- and supercritical rotations . for the latter , a minimum at @xmath35 appears due to the centrifugal nature of the harmonic term.[fig:1],width=288 ]    llllll @xmath36 ( hz ) & @xmath37 ( hz ) & @xmath38 & @xmath39 ( jm@xmath40 ) & @xmath41 & @xmath42 + @xmath43 & @xmath44 & @xmath45 & @xmath46 & @xmath47 & @xmath48 +
in this section we briefly discuss the meaning and the implementation of the semiclassical approximation for a system of particles possessing one - body hamiltonian @xmath49 with the eigenvalue equation @xmath50 where @xmath51 and @xmath52 denote the @xmath53-th single - particle eigenstate and eigenvalue , respectively .    in view of studying the various thermodynamic properties of the system ,
one is usually interested in evaluating sums of the form @xmath54 where @xmath55 is some function of the single - particle energy levels .
the semiclassical approximation consists in considering the energy levels so close to one another that the sum over them can be taken as an integral over a phase space density : @xmath56 to this end we define the density of states @xmath57 where @xmath58 is the dirac delta distribution . in the following
we investigate for the parameters of the table [ tab:1 ] whether it is justified to approximate the density of states semiclassically according to @xmath59,\label{eq : g-1b}\ ] ] where @xmath60 denotes the single - particle classical hamiltonian . in systems where the exact energy spectrum is known ,
e.g. a harmonic trap , the validity of this approximation can be determined by comparing and . however , in the case that the energy spectrum is unknown , as for the present anharmonic potential , we must take a different approach to justify the semiclassical approximation . to this end , we determine not only the leading semiclassical approximation @xmath61 , but also its subleading correction @xmath62 . after setting the relevant energy scale of the system , the validity of the semiclassical approximation can be determined by the _ relative _ magnitude of the correction .
we first treat the general case and , then , specialize to the anharmonic trap of the paris experiment .      to obtain the semiclassical expansion of the density of states of the anharmonic potential , we follow the general procedure worked out in ref .
@xcite , which derives the density of states in form of a gradient expansion @xmath63 the first term corresponds to the semiclassical approximation @xmath61 in , which yields after integration over the momenta @xmath64\dfrac{\left[e - v({\bf x})\right]^{\frac{d}{2}-1}}{\gamma(d/2)}.\label{eq : approx}\ ] ] the second term is the subleading correction to the density of states @xmath65^{\frac{d}{2}}\dfrac{1}{\gamma(d/2 - 2)}\nonumber\\ & & \times\int d^{d}x\,\theta\left[e - v({\bf x})\right]\nabla^{2}v({\bf x})\left[e - v({\bf x})\right]^{\frac{d}{2}-3}.\nonumber\\ \label{eq : g1}\end{aligned}\ ] ] we note that , for the physically relevant dimension @xmath66 , the spatial integral diverges due to the behavior of the term @xmath67^{\frac{d}{2}-3}$ ] at the classical turning points .
this problem can be solved through the method of analytic continuation as demonstrated later by applying to the anharmonic potential .
it remains to determine the energy scales for which the semiclassical approach is supposed to be valid .
at high enough temperatures , the particles are classical .
when the temperature is gradually lowered , or other parameters are changed , classical statistical mechanics is no longer exact . at first
, a semiclassical approach may work , but for certain conditions it may completely fail .
the temperature , particle number , and trap potential parameter scales in which the semiclassical approximation is valid depend on the underlying statistics . for bosons in the gas phase ,
the thermal distribution is characterized by the thermal energy @xmath68 . if this energy is sufficiently larger than the the energy of the single - particle ground - state , the leading semiclassical approximation is valid . at the bec critical temperature @xmath69
, however , the ground - state occupancy becomes macroscopic .
this effect must be treated outside of semiclassics because it corresponds to a broken symmetry state associated to the divergence of the bose distribution @xmath70 for @xmath71 .
we will see explicitly in the case of the anharmonic oscillator that the subleading correction is not negligible for energies near the ground - state energy .
the thermal energy at the critical temperature is thus the relevant energy scale for determining the validity of this approach .    for fermions
the considerations are slightly different . due to the pauli exclusion principle , as @xmath72 the particles in the system do not all drop to the ground - state but rather fill every energy level up to the fermi energy @xmath73 , which is defined as the energy of the highest occupied state .
as the density of states @xmath74 typically increases with a power of @xmath75 , the majority of particles in the system will have energies much larger than the ground - state value and their average energy is of order @xmath73 , even at @xmath76 , provided @xmath42 is large .
therefore , in principle , a spin - polarized fermi gas can be well described at all temperatures in the semiclassical approximation provided the particle number is large enough .
the relevant energy determining the validity of this assertion is simply the fermi energy @xmath73 .
numerical calculations indicate that this is indeed the case @xcite .      to consider the anharmonic potential
, we must first determine explicitly the leading and subleading terms of the semiclassical density of states .
we can , then , apply the leading term to calculate the fermi energy , and with this result return to verify the a priori assumption of the validity of the approximation .
the form of the semiclassical density of states in the anharmonic potential can be found from . performing the spatial integral
we obtain the result @xmath77 where it is meant in this and other equations , that the function takes on nonzero values for @xmath78 only for super - critical rotations .
the density of states for the parameters of the anharmonic trap described in table [ tab:1 ] is shown at several rotation frequencies @xmath1 in fig .
[ fig:2 ] . for subcritical rotations ,
the density of states vanishes for @xmath78 . for supercritical rotations ,
the density of states is nonzero also for @xmath79 , due to the dip in the potential .     for the anisotropy and the anharmonicity of the trap .
the solid lines are from bottom to top at rotations of @xmath80 to @xmath81 in steps of @xmath82 .
the dashed line is the density of states at the critical rotation @xmath83.[fig:2 ] ]    there are two important limiting cases which shed some light on the behavior of the system .
the first one is the harmonic limit @xmath84 , for which we have from @xmath85 which corresponds to the known result @xcite .
note that the harmonic limit is valid only for subcritical rotations ; for supercritical rotations the density of states becomes negative corresponding to an instability due to the fact that particles escape the trap .
the second physically important case is the pure quartic limit in radial confinement .
it is realized by either @xmath86 , i.e. a critical rotation , or @xmath87 , so that the density of states takes the form @xmath88 examining and , we see that the rotation parameter @xmath89 appears in the density of states only in ratios to the anharmonicity parameter @xmath90 . for supercritical rotations , we will find that there is a symmetry between the range of rotations , @xmath91 ( @xmath92 ) to @xmath86 ( critical rotation ) , and the range of anharmonicities from small @xmath90 to @xmath93 .
the supercritical properties of the system can thus be seen as an interpolation between the harmonic limit and the pure quartic limit .
we now calculate the subleading semiclassical approximation @xmath62 as given in for the particular case of the anharmonic potential . a direct integration in for @xmath66 yields a divergent result . however , we can rewrite the spatial integral in terms of the beta function @xmath94 for which in ref .
@xcite provides the identity @xmath95 thus , the formally divergent integrals can be treated as finite based on the analytic continuation of the gamma function through the relationship @xmath96 .
to evaluate the integral we need the laplacian of the anharmonic potential , which reads : @xmath97 with this and the full potential given in , we obtain from the first correction to the semiclassical density of states : @xmath98 where , again , the lower line is valid only for supercritical rotations .    to determine the effects of the trap parameters on the relevance of the correction term , we examine the ratio @xmath99 .
we first see that it is in general reduced in both the sub- and supercritical regimes as @xmath100 increases .
likewise , the importance of the correction term decreases as the anharmonicity @xmath90 decreases .
as the correction is largest for low energies , the most interesting case for supercritical rotations is when @xmath78 , where we have    @xmath101^{2 } , \;\;\;\;-\delta e
< e<0.\label{eq : correction - ratio}\ ] ]    from we see that for a fixed @xmath102 , the relevance of the correction increases linearly with the rotation parameter @xmath89 .
furthermore , we see that the correction diverges as @xmath103 .
however , this divergence is only of the form @xmath104 , and so any interesting integral of the form over the density of states will still have a finite contribution from this correction , provided @xmath105 is not itself divergent . for bose systems , @xmath105 contains a factor of @xmath106^{-1}$ ] .
the phase transition to bose - einstein condensation happens as @xmath107 , where @xmath108 is the lowest energy eigenvalue .
thus , in the condensate phase , an integral of the form has an improper lower limit for @xmath105 blows up as @xmath109 , since it is the lowest available energy .
this problem can be overcome by dividing the gas explicitly into a condensate plus thermal particles .
we can now check the validity of the semiclassical approximation in the case of the anharmonic trap .
both the leading and the subleading term are plotted in fig .
[ fig:3 ] for the parameters of table [ tab:1 ] , where we see that even for supercritical rotations the subleading correction becomes negligible on energy scales of more than a few @xmath110 . in the following section we determine the fermi energy @xmath73 of the system , from which we will see that , for the chosen parameters , the difference between @xmath73 and the minimum potential is orders of magnitude larger than @xmath110 . with this
, we justify the semiclassical approximation to the density of states .    .
[ fig:3],width=302 ]
though not fully accessible from the experimental point of view , the ground - state and its properties are crucial for understanding the system .
the fermi energy , e.g. , carries information about aspects like dimensionality , trap geometry , and particle number and , thus , providing information necessary to handle important physical issues .    in this section , we consider some properties of the ground - state like the dependence of the fermi energy on the particle number , the ground - state energy , the spatial and momentum distribution and , lastly , the ballistic expansion of the system .
the fermi distribution is given by @xmath111 where @xmath112 is the chemical potential and @xmath113 with @xmath114 the boltzmann constant and @xmath115 the temperature of the system . at @xmath76 ,
i.e. , @xmath116 , the fermi distribution becomes a step function at the fermi energy , which is defined as the chemical potential at @xmath76 : @xmath117    the ground - state particle number is therefore given by @xmath118 evaluating this integral with the semiclassical density of states from yields @xmath119 \label{eq : n-1a}\ ] ] for the first case , in which one has @xmath120 .
for the second case , in which @xmath121 , the particle number is given by @xmath122 the inverse of ( [ eq : n-1a]-[eq : n-1b ] ) determines the fermi energy as a function of the trap parameters and the particle number . in the harmonic @xcite and pure quartic limits we have @xmath123 respectively .
we note that the fermi temperature @xmath124 is of the same order of magnitude as the critical temperature for non - interacting bose - einstein condensates in the quartic regime @xcite @xmath125^{\frac{2}{5 } } , \label{eq : eflimiting_bec}\ ] ] with @xmath126 .
it should be pointed out that negative fermi energies are allowed by ( [ eq : n-1a]-[eq : n-1b ] ) for supercritical rotations , in which case the gas is confined entirely to the dip in the potential , creating a hole in the center of the trap analogous to the one found in two - dimensional becs @xcite . in this interesting regime ( @xmath121 )
we see that the fermi energy also takes on a simple analytic form , @xmath127 where the superscript do stands for donut , in reference to the shape that the gas assumes when all the particles are in this energy interval , as discussed below .
notice that in the ` donut ' regime , for a fixed particle number the fermi energy is set a distance above the trap minimum that is independent of the rotation frequency , but is directly related to the anharmonicity and to the number of particles .
we now have three cases where the fermi energy takes a simple analytic form corresponding to the three interesting rotation regimes : the harmonic limit for subcritical rotations , the pure quartic limit for critical rotations , and the ` donut ' regime for supercritical rotations .
although we can not obtain from ( [ eq : n-1a]-[eq : n-1b ] ) a general explicit expression for the fermi energy , we can determine its behavior from these special cases . in fig .
[ fig:4 ] we see that for low @xmath73 , the general supercritical behavior approaches the harmonic oscillator behavior . as rotation increases to the critical value , the curves approach the pure quartic limit .
as the rotation becomes supercritical , the curves are offset to more negative energies and the ` donut ' regime eventually sets up .    in fig .
[ fig:4 ] we also see that for any fermi energy possible , the quantum correction term to the semiclassical approximation is negligible , thus justifying our a priori use of the semiclassical approximation .    in the following
we use the fermi energy at zero rotation as a convenient energy scale . for our parameters
it is@xmath128 the fermi temperature @xmath129 at zero rotation provides a corresponding temperature scale , @xmath130      the total internal energy of the system in the ground - state is given by @xmath131 evaluating this integral with the density of states from yields @xmath132\\ w(e_{f } ) \end{cases } \label{eq : u-1}\ ] ] with the remark that the upper line is valid for @xmath120 and the lower one for @xmath133 and the abbreviation @xmath134 in the harmonic and purely quartic limits , respectively , the internal energy it is related to the particle number and fermi energy by @xmath135 there is also a simple relation in the ` donut ' regime , @xmath136    we see in fig .
[ fig:5 ] that the internal energy for subcritical rotations is bounded by the @xmath84 harmonic oscillator limit and the @xmath93 pure quartic limit .
this is due to the symmetry between the anharmonicity and the rotation , which causes the latter to enhance the former .     for fixed @xmath42 and the parameters of table [ tab:1 ] .
the dashed horizontal line is the harmonic limit and the dot - dashed horizontal line is the quartic limit .
the vertical lines indicate the critical rotation ( @xmath137 ) and the start of the ` donut ' regime , where @xmath138.[fig:5],width=302 ]      to obtain the particle density and the momentum distribution in the semiclassical approximation , we start from the wigner function @xcite in the semiclassical approximation @xmath139\ } } , \label{wig_func}\ ] ] which delivers the particle density through @xmath140 and the momentum distribution through @xmath141    at @xmath76 the wigner function becomes a heaviside theta function and is given by @xmath142 .
\label{wig_func_zerot}\ ] ]    in this limit , both the particle density and the momentum distribution assume a well - defined shape .
the particle density is given for a general velocity independent potential @xmath143 by @xmath144^{\frac{3}{2}}}{3\pi^{2}a_{z}^{3}\epsilon_{z}^{\frac{3}{2}}}\theta\left[e_{f}-v({\bf x})\right].\label{eq : config - dist-1}\ ] ] from , we see that a given rotation @xmath145 , at which the hole in the trap center appears , the system looks like a ` donut ' , so to speak .
the value of @xmath145 can , therefore , be determined from by setting @xmath146 and yields @xmath147 in an effectively two - dimensional ( uniform over an axial length @xmath148 ) , interacting bec trapped by the radial components of the anharmonic potential , rotation can lead as well to a hole in the center of the trap @xcite , where , in the thomas - fermi approximation , a similar analytic expression can be written down : @xmath149 where @xmath150 denotes the s - wave scattering length . in the case of non - interacting bosons ,
@xmath151 coincides with @xmath152 , but this is not true for fermions , where an effective interaction arises from the fermi statistics and @xmath145 does not equal @xmath2 even in the absence of interactions .
a discussion of this pauli `` pseudopotential '' in the case of harmonically trapped fermi gases is provided by ref .
the value of @xmath145 could also be obtained from the length scales @xmath153 set up by the fermi energy , describing the extent of the gas in space , @xmath154\right|_{x_{i}=r_{i}^{\rm tf}}=0 , \label{eq : fermir-1}\ ] ] where @xmath155 range over @xmath156 .
these length scales describe the maximum extent of the coordinate space distributions at @xmath76 and are called the thomas - fermi radii @xmath157 . for the anharmonic potential , they are given by @xmath158 } , \label{eq : r - tf}\end{aligned}\ ] ] where the positive root @xmath159 gives the outer radius of the distribution and the negative root has a physical meaning only in the ` donut ' regime , where it gives the inner radius .
we note that for supercritical rotations ( @xmath138 ) , the density distribution has the shape of a shell , of which we then designate the outermost limit in the @xmath160-plane as the thomas - fermi radius @xmath159 and the largest extention in the @xmath11-direction as @xmath157 . in the harmonic limit ,
the thomas - fermi radii are @xmath161 in particular , we are interested in the optical density which can be directly measured by absorption images . performing an integration over the axial @xmath11-direction , we obtain @xmath162^{2}}{4\pi a_{z}^{2}\epsilon_{z}^{2}}.\label{eq : config - dist-2}\ ] ] the rotation dependence of the @xmath76 column density is shown in fig .
[ fig:6 ] for the anharmonic potential . just above the critical rotation
the distribution becomes nearly homogeneous , and at high supercritical rotations the gas eventually enters the ` donut ' regime , where the gas is entirely confined to the dip in the trap .
in contrast , without the quartic trapping the gas spreads out to become an infinite pancake as the rotation approaches the critical rotation , as we can see from the divergence of the harmonic oscillator thomas - fermi radius at the critical rotation .
spatial density distribution for the parameters of the paris experiment at increasing rotations .
the anharmonic term prevents the formation of a ` pancake ' at near - critical rotations and causes the gas to take on a ` donut ' shape for large supercritical rotations.[fig:6],width=298 ]      the fermi energy also sets a momentum scale through the fermi momentum @xmath163}$ ] .
it represents the maximum extent of the momentum distribution at @xmath76 . for the anharmonic trap we have @xmath164/{\epsilon_{z}}}$ ] .
a simple dependence on the particle number emerges again in the harmonic limit @xmath165    at @xmath76 the momentum space distribution for the anharmonic potential reads @xmath166+\sqrt{q } ; & q\geq0\\ ( q+1)\pi ; & q\leq0\end{cases},\nonumber\\ \label{eq : mom - distr}\ ] ] with the dimensionless variable @xmath167 .
\label{eq : xdef}\ ] ] we see that , in contrast to bec , where the momentum distribution of the condensed atoms inherits the anisotropy of the trap , the fermionic momentum distribution of a non - interacting gas remains isotropic despite the trap anharmonicity .
this feature is drastically affected by presence of the anisotropic dipole - dipole interaction , as then the fock term of the two - body mean - field energy induces a characteristic deformation of the momentum distribution @xcite .
note that we have only considered the co - rotating frame ; as we will see , the rotation leads to an anisotropy in the momentum distribution in the lab frame .
measurements in cold gases are usually carried out by turning off the trap and shining a resonant laser on the sample after some time of free expansion .
these are called time - of - flight measurements ( tof ) .
so far we have treated the problem in the corotating frame . in considering tof measurements ,
it is necessary to consider the laboratory frame , where the gas expands ballistically .
the transformation of the momenta from the corotating frame to the lab frame yields @xmath168 where @xmath169 is the tangential unit vector . since the trap is radially symmetric , the coordinate transformation of the position does not affect the distribution .
we have then @xmath170 thus although the momentum distribution does not inherit the anisotropy of the trap , the trap rotation leads nonetheless to an anisotropic momentum distribution .    to study the expansion in detail we apply the ballistic substitution @xmath171 which yields a relation between the lab frame distribution and the corotating distribution : @xmath172,\ ] ] where @xmath173 and @xmath174 denote respectively the radius and tangential vector at the position @xmath175 , and can be expressed as @xmath176 with @xmath177 the unit vector in the axial direction .
of particular experimental interest is the density distribution at a time @xmath178 after the trap is turned off , @xmath179.\label{eq : ballistic}\end{aligned}\ ] ] for the case of harmonic trapping and zero rotation , the kinetic term can be combined with the harmonic term and the integral over momentum space in takes nearly the same form as the integration over momentum space in .
unfortunately , in the case of the anharmonic potential with rotation , the integration can not be cast into a simple form and does not provide an insightful analytic expression .    a measure of the spatial density is given by the aspect ratio @xmath180 with the respective widths @xmath181^{\frac{1}{2}}.\label{eq : wi}\ ] ] for harmonically confined fermi gases with zero rotation , it is known that anisotropic samples become asymptotically isotropic for large expansion time @xcite , since their widths scale in time as @xmath182 and the aspect ratio as @xmath183    in the case of a rotating gas with additional quartic trap , it is interesting to see how the rotation affects the asymptotic value of the aspect ratio , which reflects the lab frame momentum distribution . in fig .
[ fig:7 ] the aspect ratio has been evaluated numerically for the anharmonic potential . comparing to the harmonic trap at zero rotation , we see that the anisotropy of the momentum distribution increases with the rotation frequency for the anharmonic trap .
it is worth mentioning that the value of @xmath184 for large asymptotically large @xmath178 depends on the rotation , so that the aspect ratio could be used to measure the actual rotation frequency in an experiment .     from top to bottom.[fig:7],width=302 ]
after having studied the effects of rotation on the ground state of an anharmonically trapped polarized fermi gas , we turn our attention to the thermodynamic properties of such a system .
the standard thermodynamical recipe is specialized to the case of the anharmonic potential and the consequences for different rotation frequencies are explored in both low- and high - temperature regimes .      in this section
we describe the thermodynamical properties of the system within the framework of the grand - canonical ensemble , where the most fundamental quantity is the grand - canonical partition function . for particles obeying fermi - dirac statistics
it is given by @xmath185 , \label{eq : zg}\ ] ] where @xmath52 , as before , is the eigenvalue corresponding to the eigenstate @xmath186 .
the grand - canonical free energy is obtained from the grand - canonical partition function through @xmath187 inserting in and making use of the semiclassical approximation , i.e. , writing the discrete sum as an integral over the density of states , we obtain @xmath188de.\label{eq : fg-3}\end{aligned}\ ] ] performing a partial integration of , one sees that the integral for the grand - canonical free energy contains the fermi distribution : @xmath189 where @xmath190 denotes the number of single - particle states with energy less than @xmath75 @xmath191 in view of using the semiclassical approximation for the density of states in , we define the incomplete fermi integral @xmath192 which is called incomplete due to the fact that its integration range is smaller than the one of the original fermi integral , defined by @xmath193 for convenience we introduce the further shorthand @xmath194 with and , we write the grand - canonical free energy of a system contained in the anharmonic potential in terms of the following fermi integrals : @xmath195 } { 2\beta^{\frac{7}{2 } } \epsilon_{z}^{\frac{5}{2 } } \kappa^{\frac{1}{2 } } }   + \dfrac{\lambda^{2}\eta\zeta_{3}^{+,{\rm inc } } } { 2\beta^{3 } \epsilon_{z}^{2}\kappa } \nonumber \\ &   &   - \dfrac{\lambda^{6 } \eta^{3}\zeta_{2}^{+,{\rm inc } } } { 24\beta^{2 } \epsilon_{z } \kappa^{2 } }   + \dfrac{\lambda^{10 } \eta^{5}\zeta_{1}^{+,{\rm inc } } } { 320\beta\kappa^{3 } } , \label{eq : fg-5}\end{aligned}\ ] ] where @xmath196 is the sign of @xmath89 , defined as @xmath197    to obtain the other thermodynamical quantities , we express the grand - canonical free energy in terms of a legendre transformation @xmath198 with the differential @xmath199 therefore , the particle number @xmath42 and entropy @xmath200 can thus be determined from the grand - canonical free energy according to @xmath201 for the particle number we obtain @xmath202 } { 2\beta^{\frac{5}{2 } } \epsilon_{z}^{\frac{5}{2 } } \kappa^{\frac{1}{2}}}-\dfrac{\lambda^{2}\eta\zeta_{2}^{+,{\rm inc}}}{2\beta^{2}\epsilon_{z}^{2}\kappa}\nonumber\\ & &   + \dfrac{\lambda^{6}\eta^{3}\zeta_{1}^{+,{\rm inc}}}{24\beta\epsilon_{z}\kappa^{2 } } , \label{eq : n}\end{aligned}\ ] ] and for the entropy we get @xmath203    the internal energy can then be obtained from yielding @xmath204 } { 4\beta^{\frac{7}{2 } } \epsilon_{z}^{\frac{5}{2}}\kappa^{\frac{1}{2}}}-\dfrac{\lambda^{2}\eta\zeta_{3}^{+,{\rm inc}}}{\beta^{3}\epsilon_{z}^{2}\kappa } \nonumber\\ & &   + \dfrac{\lambda^{6}\eta^{3}\zeta_{2}^{+,{\rm inc}}}{6\beta^{2}\epsilon_{z}}-n\delta e. \label{eq : u}\end{aligned}\ ] ]    notice that , in the pure quartic limit @xmath205 , a simple relation between the free and the internal energies arises @xmath206    the heat capacity of the system at constant particle number and trap parameters is given by @xmath207 to evaluate this derivative , we must , as usual , obtain the temperature derivative of the chemical potential from the with we then obtain    @xmath208}{8\beta^{\frac{5}{2 } } \epsilon_{z}^{\frac{5}{2 } } \kappa^{\frac{1}{2 } } } - \dfrac{3\lambda^{2}\eta\zeta_{3}^{+,{\rm inc } } } { \beta^{2}\epsilon_{z}^{2 } \kappa } + \dfrac{\lambda^{6 } \eta^{3}\zeta_{2}^{+,{\rm inc } } } { 12\beta \epsilon_{z}\kappa^{2 } } \nonumber\\ &   & -\dfrac{\left\{\dfrac{5\sqrt{\pi}\left [ 2\theta(-\eta ) \zeta_{\frac{5}{2}}^{+ } + \sigma(\eta ) \zeta_{\frac{5}{2}}^{+,{\rm inc } } \right]}{\beta^{\frac{5}{2}}\epsilon_{z}^{\frac{5}{2}}\kappa^{\frac{1}{2 } } } - \dfrac{4\lambda^{2}\eta\zeta_{2}^{+,{\rm inc } } } { \beta^{2 } \epsilon_{z}^{2}\kappa } + \dfrac{\lambda^{6 } \eta^{3}\zeta_{1}^{+,{\rm inc } } } { 6\beta \epsilon_{z } \kappa^{2 } } \right\}^{2 } } { \dfrac { 8\sqrt{\pi}\left[2\theta(-\eta ) \zeta_{\frac{3}{2}}^{+ } + \sigma(\eta)\zeta_{\frac{3}{2}}^{+,{\rm inc } } \right ] } { \beta^{\frac{5}{2 } } \epsilon_{z}^{\frac{5}{2 } } \kappa^{\frac{1}{2 } } }   - \dfrac{8\lambda^{2 } \eta\zeta_{1}^{+,{\rm inc}}}{\beta^2\epsilon_{z}^{2}\kappa}}\nonumber\\ \label{eq : cv - sub}\end{aligned}\ ] ]    since this analytic expression is not particularly illuminating , we analyze the behavior of for a range of temperatures from fig .
[ fig:8 ] . at low temperatures
the heat capacity is nearly linear in temperature , and at high temperatures the heat capacity approaches @xmath209 for a fixed rotation . at high temperatures
a dip appears with a minimum occurring at a rotation that increases with the temperature .
this interesting case and the general high - temperature behavior of the heat capacity will be discussed in detail below .    , the fermi energy of the system at zero rotation ( @xmath210 . the temperature scale is the fermi temperature , @xmath211 .
the angular frequency scale is the critical frequency @xmath2.[fig:8],width=302 ]      in the low - temperature regime , we can understand temperature effects by using the sommerfeld expansion as sketched in ref .
@xcite . to this end
we consider the incomplete fermi function with @xmath212 @xmath213 where @xmath214 is the fermi distribution . a partial integration yields @xmath215,\label{eq : fermi - som-1a}\end{aligned}\ ] ] where the derivative of the fermi distribution is given by @xmath216\left[e^{-\beta(e - x ) } + 1 \right ] } .\label{eq : df}\ ] ]
consider the case that @xmath217 .
we see that the derivative becomes very small . if @xmath218 , then the lower limit of integration in can be extended to @xmath219 and @xmath220 . if @xmath221 then the integrand vanishes and @xmath222 .
we now expand @xmath223 in a taylor series about @xmath224 , @xmath225\frac{(e - x ) ^{n}}{\gamma(n+1)}\nonumber\\ \label{eq : fermi - taylor}\end{aligned}\ ] ] and integrate with to obtain @xmath226 it is worth remarking that the sommerfeld expansion ( [ eq : fermi - inc - low ] ) is valide for @xmath227 and vanishes for @xmath228 .
furthermore , the @xmath229 are constants given in terms of the riemann zeta function @xmath230 as @xmath231\zeta(2n).\ ] ] for the case of the complete fermi integral , eq .   reduces to the usual sommerfeld expansion given by @xmath232 for @xmath227 .
although the series in and terminate for @xmath233 , the expansion is not exact since the limits of integration were extended in .
it is useful to obtain an explicit equation for the temperature dependence of the chemical potential using the sommerfeld expansion and the particle number equation . to this end , we substitute @xmath234 and keep terms linear in @xmath58 as well as quadratic in temperature to obtain up to third order contributions @xmath235 ( k_{b}t)^{2 } , \label{eq : mulow}\ ] ] where @xmath236 is the density of states evaluated at the fermi energy .
we see that for low temperatures the chemical potential decreases with the square of the temperature .
this result is compared to the numerical inversion of the particle number equation with both the exact evaluation of the fermi integrals and with the evaluation using terms quartic in temperature from the sommerfeld expansion in fig .
[ fig:9 ] .
.[fig:9],width=302 ]      in the high - temperature limit we expect quantum statistics for both bosons and fermions to approach classical statistics .
here we derive the high - temperature expansions of the thermodynamic properties of a fermi gas in the anharmonic trap of the paris experiment .
we also consider in more detail the unique rotation dependence of the heat capacity in the high - temperature limit .
a high - temperature expansion of the incomplete fermi integrals is obtained as follows .
we rewrite as @xmath237 where @xmath238 . as seen in fig .
[ fig:9 ] , if the temperature is high enough , @xmath239 and , thus , @xmath240
so that we can use the binomial expansion for the factor @xmath241 . then we obtain @xmath242 with the definition @xmath243 where @xmath244 is the gamma function .    for @xmath245 in the high - temperature limit ,
it is sufficient to keep only the first term in the sum , yielding @xmath246 for positive integer @xmath53 , @xmath247 can be expanded exactly : @xmath248    the high - temperature expansion is compared to the exact evaluation of the fermi integrals in fig
.  [ fig:9 ] in the context of a numerical inversion of the particle number equation to solve for the chemical potential .
the approximation is in general very accurate for temperatures larger than half the fermi temperature .    from
we can determine the high - temperature behavior of the general expressions for the thermodynamic properties of the system .
for the particle number , we get @xmath249^{\frac{1}{2 } }   \left[2 \theta(-\eta ) + \sigma(\eta)q\left(\frac{5}{2},\beta{\delta e}\right)\right]%\right\ } + \cdots .
\label{eq : hight - n}\ ] ] it should be remarked that ( [ eq : hight - n ] ) is to be understood as giving the temperature dependence of the chemical potential in the high - temperature regime , so that @xmath112 is changed to maintain a constant value of @xmath42 on the left - hand side of the equation .
note that the high - temperature chemical potential can be determined explicitly by taking the logarithm of .
now we are in a position to write down the high - temperature versions of the internal energy @xmath250}+\cdots \label{eq : hight - u}\ ] ] and of the heat capacity @xmath251\dfrac{\sigma(\eta)(\beta\delta e)^{\frac{1}{2}}\pi^{-1/2 } e^{-\beta{\delta e}}}{\left[2\theta(-\eta)+\sigma(\eta)q\left(\frac{5}{2},\beta{\delta e}\right)\right]}\nonumber\\ & & -\dfrac{\beta{\delta e}e^{-2\beta{\delta e}}}{\pi\left[2\theta(-\eta)+\sigma(\eta)q\left(\frac{5}{2},\beta{\delta e}\right)\right]^{2 } } + \cdots\;. \label{eq : hight - cv}\end{aligned}\ ] ]    the high - temperature behavior of the heat capacity is shown in fig .  [ fig:10 ] .
for subcritical rotations we have @xmath252 , for critical rotations @xmath253 , and for supercritical rotations @xmath254 . for a fixed rotation ,
the heat capacity approaches in all cases @xmath255 as @xmath256 . for a large constant temperature
, the heat capacity also approaches @xmath255 in the high rotation limit , @xmath257 .
the most interesting feature of the high - temperature heat capacity is the dip found in the supercritical region
. differentiation of with respect to @xmath1 shows that the minimum heat capacity occurs at @xmath258 .
interestingly , the value of the minimum heat capacity , @xmath259 , is independent of temperature and rotation .
the relation between the temperature of the system and the rotation @xmath260 corresponding to the minimum heat capacity can be written more clearly as    @xmath261    the rotation @xmath260 at which the minimum heat capacity is located increases monotonically with temperature and has the asymptotic form @xmath262 .
thus , although the minimum is present at any temperature , it occurs at increasingly high rotations so that the heat capacity flattens to @xmath255 as @xmath256 for any finite rotation .
furthermore , the similarity of to and , which give the critical rotation for the transition to the zero temperature donut regime , suggests that the minimum of the heat capacity is related to the transition of the gas into a configuration where the particles are mostly contained in the minimum of the supercritical potential . recalling that @xmath263 is the potential height at the trap mininum
, one should expect that the `` donut regime '' is realized when the thermal energy of the particles is of the order of @xmath263 .
thus , setting @xmath264 should provide an estimate of @xmath265 .
indeed , this gives @xmath266 which is not bad at all .     and
@xmath267 with . from eq .
( [ eq : cv - min ] ) , one obtains that the minimum of the heat capacity at @xmath267 is located at @xmath268.[fig:10],width=302 ]      at high enough temperatures , the quantum mechanical effects may be neglected .
thus , the high - temperature properties of systems in an anharmonic trap can also be explained in the context of classical statistical mechanics .
to do so , one has to average the desired quantities over an ensemble . for a function @xmath269 of the phase space points such an average @xmath270
is given explicitly through @xmath271 where @xmath272 is again the classical hamiltonian .    in this context , the classical virial theorem states that for @xmath42 particles in thermal equilibrium at temperature @xmath115 in a potential @xmath273 , the classical ensemble averaged potential energy @xmath274 and translational kinetic energy @xmath275 of a particle in the system are related according to @xcite : @xmath276 where @xmath277 are the components of the position vector of a particle and @xmath278 .
for a cylinder symmetrical potential of the form @xmath279 , one has @xmath280 the anharmonic potential , as given in , can be decomposed according to @xmath281 with the terms @xmath282,\quad%\label{eq : pot-3}\\ v_{z}(z ) = \dfrac{\epsilon_{z}z^{2}}{2a_{z}^{2}}.\label{eq : pot-3b}\end{aligned}\ ] ] using , we obtain for the @xmath11-degrees of freedom the familiar harmonic potential result @xmath283 for the radial degrees of freedom we find correspondingly @xmath284 the total internal energy @xmath285 of a system of @xmath42 non - interacting particles in the potential @xmath273 is then given by @xmath286 before we consider the general result , let us make some remarks about eq .  .
by using the radial symmetry of the problem , one can express the average square radius through @xmath287 with the momenta @xmath288}.\ ] ]    with the help of ref .
@xcite , the limiting cases can be analytically explored . in the case of a pure quartic trap ( @xmath86 ) , the average radius is given by @xmath289 showing , as expected , that @xmath90 sets up the important length scale .
for the case of a pure harmonic trap , i. e. @xmath290 , one has @xmath291    from , we conclude that the average position of the particles diverges as @xmath89 approaches @xmath292 .
nevertheless , the energy per particle , as given from , remains @xmath293 .
in general we obtain @xmath294\sqrt{\pi}}\right\}.\end{aligned}\ ] ] substituting this into we recover .
we , thus , see that the classical statistical treatment of the system is equivalent to the high - temperature expansion which we have developed by keeping only the leading term in .
we have considered non - interacting spin - polarized fermi particles in an anharmonic trap in the presence of a centrifugal potential due to rotation of the system . applying a semiclassical treatment to the density of states ,
we have worked out the temperature dependence of several thermodynamical properties . at absolute zero
we have studied the density of particles and the momentum distribution for an anharmonically trapped system .
the analysis of the particle density shows that a well defined geometry characterizes the ground state and has revealed a tendency for the particles to move away from the center .
an analytical expression was obtained for the value of the rotation frequency @xmath0 at which a hole is created at the center of the trap , giving rise to a regime where the gas is ` donut'-shaped .
the frequency @xmath0 was shown to depart from @xmath2 , which characterizes the stability limit for harmonic confinements , due to the interplay between the fermionic statistical interaction and the anharmonicity , in close analogy to repulsively interacting bosons . at finite temperature ,
expressions for the chemical potential , the internal energy , and the heat capacity were given .
these were followed by discussions of their dependence on both anharmonicity and rotation frequency . at high - temperature ,
a virial theorem approach was presented and compared to previous results as well as shown to be particularly useful in the intuitive limits of pure harmonic and quartic radial confinements .
the most natural extension of this work would be the inclusion of interactions . in the low - temperature regime ,
short - range interactions are strongly suppressed for fully polarized fermions . to overcome this problem
, one could include other zeeman states by using an optical trap @xcite , yielding a spinor fermi gas .
issues as ( anti-)ferromagnetism and magnetization dependence of the different quantities @xcite could be addressed in the presence of the additional quartic confinement . as for long - range interactions , no further zeeman states are necessarily required .
due to its anisotropy , the dipole - dipole interaction is particularly interesting .
atoms with non - zero nuclear spin are possible constituents for systems with large magnetic dipole moment @xmath295 .
degenerate stable , for @xmath296cr ( @xmath297 ) @xcite , and meta - stable states , for @xmath298yb ( @xmath299 ) @xcite , are believed to be achievable . to date , only studies in two - dimensional systems consisting of a few atoms have already been performed @xcite , so that a three - dimensional many - body analog could be thought of .    a further , and even simpler extension of the present work could be achieved by investigating the two - dimensional limit ( @xmath300 ) of the present system . by an appropriate tuning of the quartic constant , a quantum hall regime could be realized making possible the observation of underlying landau - levels through a plateau structure in the particle density when @xmath301 @xcite . for higher rotation frequencies
, one would expect this structure to disappear and a hole to be created at a frequency value analogous to @xmath0 .
we acknowledge financial support from the german academic exchange service ( daad ) and from the german research foundation ( dfg ) within the collaborative research center sfb / tr12 symmetries and universalities in mesoscopic systems .


a _ unit disk graph _
@xmath7 is a graph whose @xmath0 vertices can be mapped to points in the plane and whose @xmath8 edges are defined by pairs of points within euclidean distance at most @xmath9 from one another .
alternatively , one can regard the vertices of @xmath7 as mapped to coplanar disks of unit diameter , so that two vertices are adjacent whenever the corresponding disks intersect .    a _ dominating set _
@xmath10 is a subset of the vertices of a graph such that every vertex not in @xmath10 is adjacent to some vertex in @xmath10
independent set _ is a subset of mutually non - adjacent vertices .
an _ independent dominating set _ is a dominating set which is also an independent set . note that any maximal independent set is an independent dominating set .    dominating sets in unit disk graphs
are widely studied due to their application in wireless ad - hoc networks  @xcite .
since it is * np*-hard to compute a minimum dominating set of a unit disk graph  @xcite , several approximation algorithms have been proposed  @xcite .
such algorithms are of two main types . _ graph - based algorithms _ receive as input the adjacency representation of the graph and assume no knowledge of the point coordinates , whereas _
geometric algorithms _ work in the real ram model of computation and receive solely the vertex coordinates as input .
if the coordinates of the @xmath0 disk centers are known , the @xmath8 edges of the corresponding graph @xmath7 can be figured out easily .
it can be done in @xmath2 time under the real ram model with floor function and constant - time hashing , and in @xmath11 time without those operations  @xcite .
thus , for the price of a conversion step , graph - based algorithms can be used when the input is a unit disk realization of @xmath7 . however ,
unless * * p**=*np * , no efficient algorithm exists to decide whether a given graph admits a unit disk realization  @xcite , let alone exhibit one . as a consequence , geometric algorithms can not be efficiently transformed into graph - based algorithms . in this paper , we introduce approximation algorithms of both types , benefiting from the same approximation factor analysis .
the proposed graph - based algorithm runs in @xmath2 time , and the geometric algorithm runs in @xmath4 time regardless of @xmath8 .
[ [ previous - algorithms ] ] previous algorithms + + + + + + + + + + + + + + + + + + +    a graph - based @xmath1-approximation algorithm that runs in @xmath2 time was presented in  @xcite .
the algorithm computes a maximal independent set , which turns out to be a @xmath1-approximation because unit disk graphs contain no @xmath12 as induced subgraphs , as shown in that same paper .
consists of a vertex with @xmath13 pendant neighbors . ]    polynomial - time approximation schemes ( ptas ) were first presented as geometric algorithms  @xcite and later as graph - based algorithms  @xcite .
also , a graph - based ptas for the more general disk graphs was proposed in  @xcite .
unfortunately , the complexities of the existing ptass are high - degree polynomials . for example , the ptas presented in  @xcite takes @xmath14 time to obtain a @xmath1-approximation ( the analysis from  @xcite ) .
although its analysis is not tight , the running time is too high even for moderately large graphs .
the reason is that these ptass invoke a subroutine that verifies by brute force whether a graph admits a dominating set with @xmath15 vertices .
the verification takes @xmath16 time , and it is unlikely that this can be improved ( unless * * fpt**=*w[1 ] * , as proved in  @xcite ) . such a subroutine is applied to several subgraphs , and the value of @xmath15 grows as the approximation error decreases .
a similar strategy was used in  @xcite to obtain a ptas for the minimum independent dominating set .
the lack of fast algorithms with approximation factor less than @xmath1 was recently noticed in  @xcite , where geometric algorithms with approximation factors of @xmath17 and @xmath18 and running times respectively @xmath19 and @xmath20 were presented .
while a significant step towards approximating large instances , those algorithms require the geometric representation of the graph , and their running times are still polynomials of rather high degrees .
linear and near - linear - time approximation algorithms constitute an active topic of research , even for problems that can be solved exactly in polynomial time , such as maximum flow and maximum matching  @xcite .
it is useful to contrast the minimum dominating set problem with the maximum independent set problem . while a maximal independent set is a @xmath1-approximation to both problems , it is easy to obtain a geometric @xmath18-approximation to the maximum independent set problem in @xmath4 time  @xcite . in the graph - based version ,
a related strategy takes roughly @xmath21 time , though .
no similar results are known for the minimum dominating set problem .
the existing ptass for the minimum dominating set problem in unit disk graphs are based on some packing constraints that apply to unit disk graphs .
one of these constraints is the _ bounded growth property _ :
the size of an independent set formed by vertices within distance @xmath22 of a given vertex , in a unit disk graph , is at most @xmath23 .
note , however , that the bounded growth property is not tight . for example
, for @xmath24 , it gives an upper bound of @xmath25 vertices where the actual maximum size is @xmath1 .
since the bounded growth property is strongly connected to the problem of packing circles in a circle , obtaining exact values for all @xmath22 seems unlikely  @xcite .
[ [ our - contribution ] ] our contribution + + + + + + + + + + + + + + + +    our main result consists of the two approximation algorithms given in section  [ s : algorithm ] : a graph - based algorithm , which runs in linear @xmath2 time , and its geometric counterpart , which runs in @xmath26 time in the real ram model , regardless of the number of edges .
the approximation factor of both algorithms is @xmath3 .
the strategy in both cases is to construct a @xmath1-approximate solution using the algorithm from  @xcite , and then perform local improvements to that initial dominating set .
our main lemma ( lemma  [ l : irreducible ] ) uses forbidden subgraphs to show that a solution that admits no local improvement is a @xmath3-approximation . since the dominating sets produced by our algorithms are independent sets , the same approximation factor holds for the independent dominating set problem .
proving that a certain graph is _ not _ a unit disk graph ( and is therefore a forbidden induced subgraph ) is no easy feat .
we make use of an assortment of results from discrete geometry in order to prove properties of unit disk graphs that are interesting _ per se_. for example , we use universal covers and disk packings to show that the neighborhood of a clique in a unit disk graph contains at most @xmath27 independent vertices .
these properties , along with a tighter version of the bounded growth property , are collected in section  [ s : forbidden ] , and allow us to show that certain graphs are not unit disk graphs .
consequently , the analyses of our algorithms employ a broader set of forbidden subgraphs which include , but are not limited to , the @xmath12 .
additionally , in section  [ s : partial ] , we show that a possible , somewhat natural refinement to our graph - based algorithm leads to a tighter @xmath5-approximation , albeit for the price of an extra @xmath28 multiplying factor in the time complexity of the algorithm .
in this section , we introduce some lemmas about the structure of unit disk graphs . these lemmas will be applied to prove our approximation factors in sections  [ s : algorithm ] and [ s : partial ] .
we start by stating three previous results from the area of discrete geometry .
the first lemma comes from the study of universal covers ( for a recent survey see  @xcite ) .
[ l : universal_cover ] if a set of points @xmath29 has diameter @xmath9 , then @xmath29 can be enclosed by a circle of radius @xmath30 .
packing congruent disks in a circle is a well - studied problem .
exact bounds on the radius of the smallest circle packing @xmath15 unitary disks are known for some small values of @xmath15 , namely @xmath31 and @xmath32  @xcite .
the bound for @xmath33 will be useful to us .
[ l : pack13 ] the radius of the smallest circle enclosing @xmath34 points with mutual distances at least @xmath9 is @xmath35 .    the _ density _ of a packing is the ratio between the covered area and the total area .
the following general upper bound is useful when no exact bound is known .
[ l : density ] every packing of two or more congruent disks in a convex region has density at most @xmath36 .    given a graph @xmath37 and a vertex @xmath38 , let @xmath39 denote the _ open neighborhood _ of @xmath40 and let @xmath41 = n(v ) \cup \{v\}$ ] denote the _ closed neighborhood _ of  @xmath40 . more generally , the _
open @xmath22-neighborhood _ of a vertex @xmath40 is the set of vertices @xmath42 such that the distance between @xmath40 and @xmath42 in @xmath7 is exactly @xmath22 , while the _ closed @xmath22-neighborhood _ of a vertex @xmath40 is the set of vertices @xmath42 such that the distance between @xmath40 and @xmath42 in @xmath7 is at most @xmath22 . for a set @xmath43 ,
we let @xmath44 and @xmath45 = n[v ] \cap s$ ] . finally , given
a subgraph @xmath46 of @xmath7 , the closed neighborhood of @xmath46 , denoted @xmath47 $ ] , is the set of vertices that belong to the closed neighborhood of some vertex of @xmath46 .
the following two lemmas concern neighborhoods in unit disk graphs .
[ l : clique ] the closed neighborhood of a clique in a unit disk graph contains at most @xmath27 independent vertices .    by lemma
[ l : universal_cover ] , the points which define a clique in a unit disk graph can be enclosed by a circle of radius @xmath30 .
therefore , the points corresponding to the closed neighborhood of such a clique are enclosed by a circle of radius @xmath48 . by lemma  [ l : pack13 ] ,
a circle enclosing @xmath34 points with mutual distances at least @xmath9 has radius at least @xmath35 . since @xmath49 , the lemma follows .
[ l:2neighborhood ] given an integer @xmath50 , the closed @xmath22-neighborhood of a vertex in a unit disk graph contains at most @xmath51 independent vertices .
all @xmath0 disks of diameter @xmath9 corresponding to the closed @xmath22-neighborhood of a vertex @xmath40 must be enclosed by a circle @xmath52 of radius @xmath53 centered on @xmath40 .
each disk of diameter @xmath9 has area @xmath54 and @xmath52 has area @xmath55 . using lemma  [ l : density ] , we have @xmath56 , and the lemma follows .
we say that a graph @xmath7 is _ @xmath57-pendant _ if there is a vertex @xmath40 in @xmath7 with @xmath15  vertices of degree  @xmath9 in the open neighborhood of @xmath40 and @xmath58 vertices of degree  @xmath9 in the open @xmath59-neighborhood of @xmath40 .
we refer to @xmath40 as a _ generator _ of the @xmath57-pendant graph .
the following lemma bounds the value of the parameter @xmath58 for @xmath60-pendant unit disk graphs .
[ l:4l - pendant ] if @xmath7 is a @xmath60-pendant unit disk graph , then @xmath61 .
let @xmath40 be a generator of @xmath7 .
since @xmath12 is a forbidden induced subgraph  @xcite and @xmath40 has @xmath17 neighbors of degree @xmath9 , we have that the remaining neighbors of @xmath40 together with @xmath40 itself form a clique . by lemma  [ l : clique ] , it follows that @xmath62 .
next , we consider the case of @xmath63-pendant unit disk graphs .
[ l:3l - pendant ] if @xmath7 is a @xmath63-pendant unit disk graph , then @xmath64 .
let @xmath40 be a generator of @xmath7 . since two vertices are adjacent if and only if their euclidean distance is at most @xmath9 , the closed neighborhood of @xmath40 lies inside a circle of radius @xmath9 centered at @xmath40 .
let @xmath65 be a neighbor of @xmath40 with degree @xmath9 .
we divide the aforementioned circle into six congruent sectors @xmath66 in such a way that @xmath65 sits on the boundary of two adjacent sectors @xmath67 .
since the diameter of each sector is @xmath9 and @xmath65 has degree @xmath9 , we have that @xmath68 only contains vertex @xmath65 . by the same argument ,
the remaining two neighbors of @xmath40 that have degree @xmath9 are contained in two of the remaining four sectors , which we call @xmath69 ( note that the sectors do not necessarily appear in the order @xmath66 ) .
therefore , only sectors @xmath70 may contain the neighbors of @xmath40 that have degree greater than @xmath9 .
notice that the neighbors of the vertices in @xmath71 should be located within euclidean distance at most @xmath9 from @xmath71 .
consequently , the disks of unit diameter corresponding to such vertices should be completely contained in the region defined as the minkowski sum of @xmath72 and a disk of radius @xmath73 , i.e.  the region within distance at most @xmath73 from @xmath72 .
we compute upper bounds to the area of such a region by considering two cases , depending on whether @xmath74 and @xmath75 are opposite sectors ( bounded by the same pair of straight lines ) .
first , we consider the case when @xmath74 and @xmath75 are _ not _ opposite sectors . in this scenario , @xmath72 is contained in a semicircle of radius @xmath9 centered at the generator @xmath40 , as represented in figure  [ f : sectors](a ) .
we define a region @xmath76 as the locus of points within distance at most @xmath73 from the aforementioned semicircle of radius @xmath9 .
the area @xmath77 of @xmath76 is therefore @xmath78 .
by lemma  [ l : density ] , the number of disks of unit diameter contained in @xmath76 is at most @xmath79 since @xmath80 of these disks are the degree-1 neighbors of @xmath40 , we have @xmath81 .
second , we consider the case when @xmath74 and @xmath75 are opposite sectors . in this scenario ,
the region @xmath72 is not convex . in order to obtain a convex region @xmath82
, we define @xmath82 as the locus of points within distance at most @xmath73 from the convex hull of @xmath72 ( see figure  [ f : sectors](b ) ) .
the area @xmath83 of @xmath82 is therefore @xmath84 .
by lemma  [ l : density ] , the number of disks of unit diameter contained in @xmath82 is at most @xmath85 since @xmath80 of these disks are the degree-1 neighbors of @xmath40 , we have @xmath86 .    finally , the following lemma holds for the general case .
[ l : kl - pendant ] if @xmath7 is a @xmath57-pendant unit disk graph , then @xmath87 .    immediately from lemma  [ l:2neighborhood ] with @xmath88 .
in this section , we present two 44/9-approximation algorithms . the key property to analyze the approximation factor
is presented in lemma  [ l : irreducible ] , while the running time analyses are presented in sections  [ s : graph - alg ] and  [ s : geo - alg ] .
hereafter , let @xmath37 be a unit disk graph , and let @xmath89 be an independent dominating set of @xmath7 .
if @xmath90 and @xmath91 , we say that @xmath40 _ dominates _ @xmath65 and , conversely , that @xmath65 _ is dominated _ by @xmath40 .
as already mentioned , unit disk graphs are free of induced @xmath12 .
therefore , at most @xmath1 vertices of @xmath10 may belong to the closed neighborhood of any given vertex @xmath38 .
corona _ is a set @xmath92 consisting of exactly @xmath1 neighbors of some vertex @xmath93 .
such a vertex @xmath94 , which is not necessarily unique , is called a _ core _ of the corona @xmath95 , whereas the @xmath1 vertices of the corona are referred to as the corona s _
petals_. notice that the subgraph induced by a corona @xmath95 and a corresponding core @xmath94 is a @xmath96 .    given a dominating set @xmath10 ,
a corona @xmath92 is said to be _
reducible _ if there is a core @xmath94 of @xmath95 such that @xmath97 is a dominating set .
we refer to the operation that converts @xmath10 into the smaller dominating set @xmath97 as a _ reduction _ of @xmath95 with respect to @xmath94 .
if there is no core allowing for a reduction of @xmath95 , than @xmath95 is dubbed _
irreducible_. if @xmath95 is an irreducible corona , then , for every core @xmath94 of @xmath95 , there must exist a vertex @xmath98 , such that :    a.   @xmath42 is not adjacent to @xmath94 ; b.   @xmath42 is only dominated , in @xmath10 , by vertices that belong to @xmath95 .
we call @xmath42 a _ witness _ of @xmath94 , conveying the idea that the corona having @xmath94 as a core can not be reduced with respect to @xmath94 due to the existence of @xmath42 .
[ l : irreducible ] let @xmath37 be a unit disk graph , @xmath10 an independent dominating set in @xmath7 , and @xmath99 a minimum dominating set of @xmath7 .
if all coronas in @xmath10 are irreducible , then @xmath100 .
we use a charging argument to bound the ratio between the cardinalities of @xmath10 and @xmath99 . consider that each vertex @xmath101 splits a _ unit charge _ evenly among the vertices in the closed neighborhood @xmath102 $ ] .
the function @xmath103}$ ] below corresponds to the total charge assigned to each vertex @xmath104 , accumulating the ( fractional ) charges that @xmath105 receives from the vertices in @xmath106 $ ] :    @xmath107 } \frac{1}{\left|n_{d^*}[u]\right|}.\end{aligned}\ ] ]    note that , since @xmath10 and @xmath99 are dominating sets , neither @xmath102 $ ] nor @xmath108 $ ] are ever empty , and @xmath109|$ ] .
such function @xmath110 allows us to write the cardinality of @xmath10 as @xmath111    since @xmath112 is precisely the average value of @xmath113 over the elements of @xmath99 , we obtain the desired bound @xmath114 by showing that the existence of vertices @xmath115 in @xmath99 with @xmath116 is counterbalanced by a sufficiently large number of vertices @xmath117 in @xmath99 with @xmath118 .
before we continue , note that @xmath116 means that @xmath119 , because the summation in ( [ eq : fv ] ) has at most @xmath1 terms , all of which are of the form @xmath120 for some integer @xmath121 .
thus , let @xmath115 be a vertex in @xmath99 with @xmath119 . clearly , @xmath122 , otherwise @xmath123| = 1 $ ] , because @xmath10 is an independent set .
moreover , @xmath115 must have exactly @xmath1 neighbors in @xmath10 , since a larger number of neighbors in @xmath10 would imply the existence of an induced @xmath12 in @xmath7 , which is not possible , and a smaller number would imply @xmath123| \leq 4 $ ] , a contradiction .
therefore , vertex @xmath115 is a core .
now let @xmath92 be the corona of which @xmath115 is a core . because there are no reducible coronas in @xmath10
, the core @xmath115 must have a witness @xmath42 .
note that , for all petals @xmath124 , the only vertex in @xmath102 $ ] is the core  @xmath115 .
otherwise , the contribution of some @xmath124 in the summation yielding @xmath125  given by ( [ eq : fv ] )  would be at most @xmath126 , and @xmath125 would be at most @xmath127 , a contradiction .
in particular , the above implies that the witness @xmath42 , which is adjacent to at least one vertex in @xmath95 , can not belong to @xmath99 .
but @xmath99 is a dominating set , so there must exist a vertex @xmath128 that is adjacent to @xmath42 , and @xmath129 because a witness @xmath42 is not adjacent to the corresponding core by definition .
we call @xmath117 a _ reliever _ of  @xmath115
. figure  [ f : witness ] illustrates this situation .
we now show that @xmath130| \leq 4 $ ] . for sake of contradiction ,
assume @xmath130| > 4 $ ] . because @xmath7 is free of induced @xmath12
, such number must be exactly @xmath1 , so that @xmath117 is the core of a corona @xmath131 .
such a corona must be disjoint from corona @xmath95 , otherwise there would be a vertex in adjacent to more than one vertex in @xmath99 , namely @xmath115 and @xmath117 , contradicting the fact that the only neighbor in @xmath99 of any petal of @xmath95 is the core  @xmath115 . since , by definition , the witness @xmath42 is only dominated in @xmath10 by vertices of @xmath95 , we have @xmath132 .
hence , @xmath133 is an independent set of @xmath7 , constituting , along with the core @xmath117 , an induced @xmath12 in @xmath7 .
this is a contradiction , because @xmath7 is a unit disk graph .
thus , @xmath130| \leq 4 $ ] . since @xmath134|$ ] , we have @xmath118 .
we have just shown that the existence of a vertex @xmath115 in @xmath99 with implies the existence of a vertex @xmath128 such that @xmath118 .
were this correspondence one - to - one , we would be able to state that the average of @xmath113 over the elements of @xmath99 was no greater than @xmath135 .
unfortunately , this correspondence is not necessarily one - to - one , as exemplified by the graph in figure  [ f : badgraph ] , for which a disk model is given in figure  [ f : badgraph_model ] with coordinates presented in table  [ t : coordinates ] .
still , the lemmas in section  [ s : forbidden ] allow us to bound the ratio between the number of vertices @xmath115 with @xmath119 and the number of vertices @xmath117 for which the values of @xmath110 are significantly lower .
let @xmath136 be a reliever .
in order to obtain the claimed bound , we consider two cases depending on the size of @xmath137 $ ] :    \(i ) @xmath138| \leq 3.}$ ] by lemma  [ l:2neighborhood ] , the closed @xmath17-neighborhood of @xmath117 contains at most @xmath139 independent vertices .
since each corona contains @xmath1 independent vertices , at most @xmath140|)/5 \rfloor = 14 $ ] cores may share a common reliever .
to derive an upper bound , let @xmath141 denote such cores .
if @xmath130| \leq 3 $ ] , then the average value of @xmath113 among @xmath142 is at most @xmath143    \(ii ) @xmath138| = 4.}$ ] by lemma  [ l:4l - pendant ] , if @xmath130| = 4 $ ] , then at most @xmath144 cores @xmath145 may have @xmath117 as their common reliever , for otherwise we obtain a @xmath146-pendant graph , which can not be a unit disk graph
. thus , the average value of @xmath113 among @xmath147 is at most @xmath148    the worst case is @xmath130| = 4 $ ] , and therefore @xmath114 , concluding the proof .
@xmath149 & @xmath150 & + @xmath151 & @xmath152 & + @xmath153 & @xmath154 & + @xmath155 & @xmath156 & +   +  @xmath157 & @xmath158 & @xmath159 +  @xmath160 & @xmath161 & @xmath162 +  @xmath163 & @xmath164 & @xmath165 +  @xmath166 & @xmath167 & @xmath168 +      by lemma  [ l : irreducible ] , an independent dominating set with no reducible coronas is a @xmath3-approximation to the minimum dominating set .
we now describe how to obtain such a set in linear time given the adjacency list representation of the graph .
we can easily compute a maximal independent set @xmath10 , which is a @xmath1-approximation to the minimum dominating set  @xcite , in @xmath2 time . an independent dominating set with no reducible coronas
can then be obtained by iteratively performing reductions .
however , naively performing such reductions leads to a running time of @xmath169 , since ( i ) there are @xmath170 candidate cores for a reducible corona , ( ii ) detecting whether a vertex @xmath40 is in fact the core of a reducible corona by inspecting the @xmath18-neighborhood of @xmath40 takes @xmath171 time , and ( iii ) we may need to reduce a total of @xmath170 coronas . fortunately , the following algorithm modifies the set @xmath10 and returns an independent dominating set with no reducible coronas in @xmath2 time .    1 .   for each vertex @xmath172 , compute @xmath173 .
2 .   for each vertex @xmath172 , if @xmath174 , add @xmath173 to the list of coronas @xmath175 ( unless it is already there ) .
3 .   let @xmath176 . for each corona
@xmath177 , if there is a vertex @xmath94 such that @xmath97 is a dominating set , then add @xmath94 to the set @xmath178 .
4 .   choose a maximal subset @xmath179 of @xmath178 such that the pairwise distance of the vertices in @xmath179 is at least @xmath1 .
5 .   for each vertex @xmath180 , perform a reduction @xmath181 .
repeat all the steps above until @xmath182 .
the algorithm is correct since all changes made to @xmath10 along its execution preserve the property that @xmath10 is an independent dominating set .
notice that , in step ( 4 ) , we only reduce coronas that are sufficiently far from each other , in order to guarantee that we do not reduce a corona that may have ceased to be reducible due to a previous reduction . moreover , the algorithm always terminates because the size of @xmath10 decreases at every iteration , except for the last one .
next , we show that the running time is @xmath2 .
step ( 1 ) can be easily implemented to run in @xmath2 time . to execute step ( 2 ) in @xmath2 time
, we must determine in constant time whether a corona is already in the list @xmath175 .
this can be achieved by indexing each corona @xmath95 by an arbitrary vertex @xmath183 ( say , the one with the lowest index ) , and by storing with @xmath40 a list of coronas that are in @xmath184 and whose index is @xmath40 .
note that , because of the packing constraints inherent to unit disk graphs , the number of coronas that contain a given vertex is constant .
step ( 3 ) can be implemented as follows ( for each corona @xmath177 ) :    1 .   let @xmath185 be the union of the open neighborhoods of the @xmath1 petals of @xmath95 .
2 .   let @xmath186 be the set of vertices dominated only be vertices of @xmath95 , i.e. , @xmath187 contains every vertex @xmath188 such that @xmath189 .
3 .   let @xmath190 be the intersection of the closed neighborhoods @xmath41 $ ] of all @xmath191 .
4 .   if @xmath192 , then add an arbitrary vertex of @xmath190 to the set @xmath178 .
the steps above take @xmath2 total time when executed for all coronas @xmath177 , because the number of coronas that contain or are adjacent to a given vertex is also constant by packing constraints .
it is easy to perform steps ( 4 ) and ( 5 ) in linear time .
it remains to show that the whole process is only repeated for a constant number of iterations .
let @xmath193 denote the set of reducible coronas at each iteration of the algorithm with @xmath194 .
note that the reductions performed in step ( 5 ) never create a new reducible corona .
therefore @xmath195 .
let @xmath95 denote a corona that was reduced in the last iteration @xmath15 .
if @xmath95 was not reduced during a previous iteration @xmath196 , then another corona within distance @xmath1 from @xmath95 was reduced at that very iteration @xmath197 . since , again by packing constraints , the maximum number of coronas within constant distance from @xmath95 is itself a constant , we have @xmath198 .
the following theorem summarizes the result from section  [ s : graph - alg ] .
[ thm : graph - alg ] given the adjacency list representation of a unit disk graph with @xmath0 vertices and @xmath8 edges , it is possible to find a @xmath3-approximation to the minimum dominating set problem in @xmath199 time .      in this section , we describe how to obtain an independent dominating set with no reducible corona in @xmath4 time given the geometric representation of the graph .
the input is therefore a set @xmath29 of @xmath0 points . without loss of generality , we assume that the corresponding unit disk graph is connected ( otherwise , we can compute the connected components in @xmath4 time using a delaunay triangulation  @xcite ) .
we use terms related to vertices of the graph and to the corresponding points interchangeably .
for example , we say a set of points is independent if all pairwise distances are greater than @xmath9 .
we want the points of @xmath29 to be structured in a suitable fashion .
thus , as a preliminary step , we sort the points by @xmath200-coordinates and by @xmath201-coordinates separately ( such orderings will also be useful later on ) , and we partition the points of @xmath29 according to an infinite grid with unitary square cells by performing two sweeps on the sorted points . without loss of generality , we assume that no point lies on the boundary of a grid cell . given @xmath202 , let @xmath203 denote the grid cell that contains  @xmath204 .
we refer to the set of at most @xmath144 non - empty grid cells surrounding a cell @xmath205
as the _ open vicinity _ of @xmath205 , denoted @xmath206 , and to the union of @xmath205 and its open vicinity as the _ closed vicinity _ of @xmath205 , denoted @xmath207 $ ] .
note that a point @xmath204 can only be adjacent to points in the closed vicinity of @xmath203 , that is , @xmath208
\subset n[\sigma(p)]$ ] .
each point @xmath202 stores a pointer to its containing cell @xmath203 . also , each cell stores the list of points it contains and pointers to the cells in its open vicinity .
since the graph is connected , the diameter of the point set is at most @xmath209 , and thus this whole step can be done in @xmath4 time .
we are now able to show how to compute a maximal independent set @xmath10 efficiently .
we begin by making a copy @xmath210 of @xmath29 , and by letting @xmath211 .
then we repeat the two following steps while set @xmath210 is non - empty .
( i )  choose an arbitrary point @xmath212 and add it to set @xmath10 .
( ii )  for each point @xmath213 in the closed vicinity of @xmath203 , remove @xmath213 from @xmath210 if @xmath214 . when @xmath210 becomes empty
, @xmath10 is an independent dominating set .
this process takes @xmath170 time due to the two following facts .
first , a cell belongs to the closed vicinity of a constant number of cells .
second , the maximum number of points inside a cell with pairwise distances greater than @xmath9 is also a constant .
we now have that @xmath10 is a maximal independent set , and therefore a @xmath1-approximation to the minimum dominating set .
next , we show how to modify @xmath10 in order to produce an independent dominating set with no reducible corona , therefore a @xmath3-approximation to the minimum dominating set .
the algorithm mirrors the one in section  [ s : graph - alg ] , but each step takes no more than @xmath4 time using the geometric representation of the graph .
since @xmath10 is an independent set and a grid cell @xmath205 has side @xmath9 , a simple packing argument shows that @xmath215 .
we store the set @xmath216 in the corresponding cell @xmath205 . in order to compute @xmath217
, it suffices to inspect the at most @xmath218 points in @xmath216 for @xmath219 $ ] .
we can then build a list of coronas in @xmath170 time ( steps ( 1 ) and ( 2 ) of section  [ s : graph - alg ] ) .    to perform step ( 3 ) , we need to find out whether there is a core @xmath94 such that @xmath97 is a dominating set , for each corona @xmath220 .
first , we make @xmath185 the union of @xmath221 for @xmath222 .
then , we let @xmath186 be the subset of @xmath185 containing only the points @xmath204 with @xmath223 .
these first two steps are similar to steps ( 3a ) and ( 3b ) in section  [ s : graph - alg ] .
the remaining sub - steps of step ( 3 ) are significantly different , though .
we proceed by making @xmath224 .
we need to determine whether there is a point @xmath225 that is adjacent to all points in @xmath190 . for each @xmath225 ,
let @xmath226 denote the disk of radius @xmath9 centered at @xmath204 .
let @xmath76 denote the convex region defined by the intersection of @xmath226 for all @xmath225 .
a point @xmath204 is adjacent to all points in @xmath190 if and only if @xmath227 .
we can compute the region @xmath76 in @xmath228 time using divide - and - conquer in a manner analogous to half - plane intersection  @xcite .
we can then test whether each point @xmath225 belongs to the region @xmath76 in logarithmic time using binary search ( remember the points were previously sorted ) .
if there is at least one point @xmath229 , then we add @xmath204 to the set @xmath178 .
therefore , the whole step ( 3 ) takes @xmath4 time .    in step ( 4 ) of the geometric algorithm , we choose an alternative set @xmath230 which can be computed in @xmath170 time as follows . for each @xmath231 , we add @xmath204 to @xmath179 and then remove from @xmath178 all points that are contained in the cells within euclidean distance at most @xmath17 of @xmath203 .
since by packing constraints there are @xmath232 points in the intersection of @xmath10 and the closed vicinity of a cell , we can easily perform step ( 5 ) in @xmath170 time . finally , the number of repetitions triggered by step ( 6 ) is constant by an argument identical to the one given for the graph - based algorithm .
the following theorem summarizes the result from section  [ s : geo - alg ] .
[ thm : geo - alg ] given a set of @xmath0 points representing a unit disk graph , it is possible to find a @xmath3-approximation to the minimum dominating set problem in @xmath4 time in the real ram model of computation .
in the previous section , a 44/9-approximation was obtained by reducing coronas of a maximal independent set @xmath10 of graph @xmath7 , that is , by subsequently replacing @xmath1 petals with @xmath9 core in @xmath10 as long as that operation preserved dominance . a natural step to tighten the approximation factor
is to allow for _ weak reductions _ ,
whereby the @xmath1 petals of a corona @xmath95 are removed from the independent dominating set @xmath10 , yet not only is a core @xmath94 of @xmath95 inserted into @xmath10 but also some mutually non - adjacent witnesses of @xmath94 , as long as their number is no greater than @xmath18 and they dominate all witnesses of @xmath94 . by doing so , the weak reduction of @xmath95 ( with respect to @xmath94 ) preserves dominance and still shaves off at least one unit from the size of @xmath10 .
if such operation is possible on a corona @xmath95 , then @xmath95 is said to be _
weakly reducible_. a core @xmath94 which has @xmath17 ( or more ) mutually non - adjacent witnesses is said to be an _ overwhelmed _ core .
is an overwhelmed core of a corona @xmath95 , then it might be the case that a weak reduction on @xmath95 with respect to @xmath94 is still possible .
if the subgraph @xmath233 $ ] induced by the set @xmath234 of witnesses of @xmath94 admits an independent dominating set @xmath235 of size no greater than @xmath18 , then @xmath236 is still an independent dominating set of @xmath7 , and its cardinality is strictly less than @xmath237 .
however , one can not decide in ( close to ) linear time whether such a dominating set exists , and that will not be required by our algorithm . ]    we consider the graph - based algorithm presented in section  [ s : graph - alg ] with some modifications to cope with weak reductions .
the whole modified algorithm becomes :    1 .   for each vertex @xmath172
, compute @xmath173 .
2 .   for each vertex @xmath172 , if @xmath174 , add @xmath238 to the list of coronas @xmath175 ( unless it is already there ) , and add @xmath40 to the list @xmath239 containing the cores of @xmath95 .
3 .   let @xmath178 be an initially empty mapping of cores onto sets of witnesses . for each corona
@xmath177 , if there is a vertex @xmath240 and an independent set @xmath241 with at most @xmath18 witnesses of @xmath94 such that @xmath242 is a dominating set , then add @xmath243 to @xmath178 .
4 .   choose a maximal subset @xmath179 of the cores in @xmath178 such that the pairwise distance of the vertices in @xmath179 is at least @xmath1 .
5 .   for each vertex @xmath180 , perform a ( weak ) reduction @xmath244 . 6 .
repeat all the steps above until @xmath182 .
the new step ( 3 ) can be implemented as follows ( for each corona @xmath177 ) :    1 .
let @xmath185 be the union of the open neighborhoods of the @xmath1 petals of @xmath95 .
2 .   let @xmath186 be the set of all vertices @xmath188 with @xmath189 .
3 .   for each core @xmath240 , greedily obtain a maximal independent set @xmath241 of @xmath245 $ ] . if @xmath246 , add @xmath243 to @xmath178 and break .
because each core is evaluated separately as to whether it allows for a weak reduction ( whereas in the algorithm of section  [ s : graph - alg ] a constant number of set operations per corona was executed ) , step ( 3c ) dominates the complexity of the whole sequence of steps ( 1 ) to ( 5 ) , with a total @xmath247 time for running on all coronas @xmath177 . another important difference , as far as time complexity goes , is that in this modified algorithm the number of iterations of the main loop  steps ( 1 ) to ( 5 )  is no longer @xmath232 , due to the fact that new ( weakly ) reducible coronas can be created , as illustrated in figure  [ f : newcoronas ] . however , the number of iterations is certainly @xmath170 , because the size of @xmath10 decreases by at least @xmath9 in each iteration .
hence the overall time complexity of the modified algorithm is @xmath6 .
next , we establish the approximation factor of the modified algorithm .
note that step ( 3c ) asserts that the coronas that are not ( weakly ) reduced by the algorithm present only overwhelmed cores .
[ l : irreducible_partial ] let @xmath37 be a unit disk graph , @xmath10 an independent dominating set in @xmath7 , and @xmath99 a minimum dominating set of @xmath7 . if all coronas in @xmath10 have only overwhelmed cores , then @xmath248 .
we follow a strategy similar to that in the proof of lemma  [ l : irreducible ] ( also using the concept of reliever defined therein ) : by employing the same function @xmath103}$ ] defined in ( [ eq : fv ] ) , we prove there is an appropriate balance among vertices with high ( @xmath249 ) and low ( @xmath250 ) images under @xmath110 , thus yielding an average value for @xmath113 that is no greater than @xmath5  the claimed approximation factor .
let @xmath115 be a vertex in @xmath99 with @xmath251 . because @xmath252 , and all ( at most @xmath1 ) terms in the summation that yields @xmath110 are either @xmath9 or no greater than @xmath253 , we have that @xmath251 implies @xmath119 .
so @xmath115 is a core . by hypothesis ,
@xmath115 is overwhelmed .
hence , @xmath115 possesses at least @xmath17 mutually non - adjacent witnesses @xmath254 , each one implying the existence of a reliever @xmath128 with @xmath118 .
surely , such relievers need not be all distinct , and different cores may share a common reliever .
still , geometric properties of unit disk graphs allow us to derive upper bounds for the core - to - reliever ratio .
we can thus obtain an upper bound to the average value of @xmath113 in @xmath99 , and consequently to the approximation factor of the algorithm .
suppose there are @xmath255 cores @xmath256 such that @xmath257 . for @xmath258 ,
we let @xmath259 be the corona having @xmath260 as a core , @xmath261 be a set of ( at least 4 ) mutually non - adjacent witnesses of @xmath260 , and @xmath262 be the set of relievers of @xmath260 .
now , we construct a bipartite multigraph @xmath263 as follows .
the parts of @xmath264 are @xmath265 and @xmath266 .
the multiset @xmath267 contains , between each core @xmath268 and reliever @xmath269 , a number @xmath270 of parallel edges that is equal to the number of petals of @xmath259 adjacent to witnesses @xmath42 ( of @xmath260 ) such that @xmath42 is a neighbor of @xmath271 .
thus , the total number of edges incident to a core @xmath268 is    @xmath272    analogously , the total number of edges incident to a reliever @xmath269 is    @xmath273    we now obtain an upper bound @xmath274 for the average value of @xmath113 over @xmath275 .
observe that @xmath276 contains all vertices @xmath115 of @xmath99 such that @xmath251 , hence the average value @xmath277 of @xmath113 over the whole set @xmath278 can not be any greater .
of course the average value we are interested in depends on the core - to - reliever ratio @xmath279 in @xmath264 : the more cores ( respectively , the fewer relievers ) in @xmath276 ( in @xmath82 ) , the greater the average . therefore , in order to obtain the desired upper bound @xmath280 , we must consider the case in which the elements of @xmath276 ( respectively , of @xmath82 ) have degrees in @xmath264 that are as low ( as high ) as possible .    first ,
notice that if @xmath115 is a core of corona @xmath95 and @xmath119 , then it is not possible that more than two non - adjacent witnesses of @xmath115 sharing a common reliever @xmath117 are adjacent to the same petal @xmath281 .
otherwise , let @xmath282 and @xmath283 be such witnesses .
since @xmath204 and @xmath117 are non - adjacent ( due to the image of @xmath115 under @xmath110 being @xmath1 ) , the subgraph of @xmath7 induced by @xmath284 is a @xmath285 .
this is a contradiction , because the graph @xmath285 is not a unit disk graph  @xcite .
thus , since @xmath115 has at least four witnesses , the number of petals of @xmath95 adjacent to witnesses of @xmath115  and therefore the degree of each core in @xmath276  is at least @xmath286 .
the above lower bound can be improved for cores @xmath287 having a with exactly four neighbors in @xmath10 .
if a reliever @xmath117 has four neighbors in ( the independent set ) @xmath10 , then the remaining neighbors of @xmath117 form a clique .
since there is a witness @xmath42 of @xmath115 among these remaining neighbors of @xmath117 , then the other ( at least ) three mutually non - adjacent witnesses of @xmath115 must have relievers distinct from @xmath117 .
moreover , by an argument analogous to the one used in the previous paragraph , those witnesses must be adjacent to at least @xmath288 petals of the corona having @xmath115 as a core .
this means there are at least two edges connecting @xmath115 to its relievers in @xmath289 , plus one edge connecting @xmath115 to @xmath117 .
hence , the degree of @xmath115 in @xmath264 is at least @xmath18 .
the maximum degrees of vertices @xmath290 depend on the number @xmath291 of neighbors of @xmath117 in @xmath10 , and now we employ some of the geometric lemmas of section  [ s : forbidden ] to infer suitable upper bounds . recall that any set of petals is an independent set .    *
if @xmath292 ( implying @xmath118 ) , then , by lemma  [ l:4l - pendant ] , there are at most @xmath144 petals in the @xmath59-neighborhood of @xmath117 . hence , @xmath293 . * if @xmath294 ( implying @xmath295 ) , then , by lemma  [ l:3l - pendant ] , there are at most @xmath296 petals in the @xmath59-neighborhood of @xmath117 .
hence , @xmath297 . * if @xmath298 ( implying @xmath299 ) , then , by lemma  [ l : kl - pendant ] , there are at most @xmath300 petals in the @xmath59-neighborhood of @xmath117
. hence , @xmath301 . *
if @xmath302 ( implying @xmath303 ) , then , by lemma  [ l : kl - pendant ] , there are at most @xmath304 petals in the @xmath59-neighborhood of @xmath117 .
hence , @xmath305 .    for @xmath306 ,
let @xmath307 denote the number of relievers in @xmath7 containing exactly @xmath15 neighbors in @xmath10 , so that @xmath308 .
let also @xmath309 be the set of cores having at least one reliever with exactly @xmath17 neighbors in @xmath10,and @xmath310 be the set of cores whose relievers have at most three neighbors in @xmath10 . finally , we let @xmath311 and @xmath312 , so that @xmath313 .    since @xmath264 is bipartite , the number of edges incident to @xmath276 and the number of edges incident to @xmath82 are the same .
consequently ,    @xmath314 and @xmath315    the same reasoning holds for the subgraph of @xmath264 induced by the cores in @xmath316 and their relievers in @xmath82 , so @xmath317 by dividing both sides of ( [ eq:3 ] ) by @xmath59 and adding it to ( [ eq:2 ] ) , we obtain @xmath318 as for the average @xmath280 of @xmath113 over the elements of @xmath264 , we can write    @xmath319    now , substituting @xmath320 in the expression above by the upper bound obtained from ( [ eq:4 ] ) , we have    @xmath321    thus , @xmath280 is bounded by a multivariate rational function whose maximum is @xmath322 , achieved when @xmath323 .
the following theorem summarizes the result from section  [ s : partial ] .
[ thm : graph - alg - partial ] given the adjacency list representation of a unit disk graph with @xmath0 vertices and @xmath8 edges , it is possible to find a @xmath5-approximation to the minimum dominating set problem in @xmath6 time .
we introduced novel efficient algorithms for approximating the minimum dominating set and minimum independent dominating set in unit disk graphs .
on one hand , a linear - time algorithm was devised attaining a sub-5 approximation factor , namely @xmath324 nevertheless , the best lower bound we know for the proposed algorithm is @xmath325 , which corresponds to the unit disk graph given in figure  [ f : badgraph ] . closing this gap would likely require the development of new tools to prove that certain graphs are not unit disk graphs , for which computer generated proofs may be useful .         on the other hand , an enhanced approximation factor of @xmath327
was obtained by allowing for more local replacements , yet a lower bound of  @xmath328 , corresponding to the unit disk graph given in figure  [ f : lowerbound425 ] ( with coordinates in table  [ t : coordinates2 ] ) , is the best we are aware of . notwithstanding the @xmath28 factor increase on its time complexity , such a modified algorithm is still much faster than , say , the state - of - the - art 4-approximation algorithm from  @xcite , which requires a geometric model as input and runs in @xmath329 time . moreover , since the number of ( weak ) reductions that are performed remains linear , it may be possible to conceive either a refined analysis or a smarter implementation .
p. christiano , j.a .
kelner , a. madry , d.a .
spielman , s .- h .
teng , electrical flows , laplacian systems , and faster approximation of maximum flow in undirected graphs , in : proc .
43rd annual acm symp . on theory of computing ( stoc ) , 2011 , pp
. 273282 .
t. erlebach , m. mihalk , a ( @xmath330)-approximation for the minimum - weight dominating set problem in unit disk graphs , in : proc .
7th workshop on approximation and online algorithms , vol .
5893 of lecture notes in computer science , 2010 , pp . 135146 .
m. gibson , i. pirwani , algorithms for dominating set in disk graphs : breaking the @xmath331 barrier , in : proc .
18th annual european symposium on algorithms ( esa ) , vol .
6346 of lecture notes in computer science , 2010 , pp . 243254 .
hunt iii , m.v .
marathe , v. radhakrishnan , s. ravi , d.j .
rosenkrantz , r.e .
stearns , nc - approximation schemes for np- and pspace - hard problems for geometric graphs , journal of algorithms 26 ( 1998 ) 238274 .
d. marx , parameterized complexity of independence and domination on geometric graphs , in : proc .
2nd international workshop on parameterized and exact computation ( iwpec ) , vol .
4169 of lecture notes in computer science , 2006 , pp . 154165 .
f. zou , y. wang , x .- h .
xu , x. li , h. du , p. wan , w. wu , new approximations for minimum - weighted dominating sets and minimum - weighted connected dominating sets on unit disk graphs , theoretical computer science 412  ( 3 ) ( 2011 ) 198208 .