Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: it has long been recognized that when an atom is situated in a material medium , the local electromagnetic field acting on it differs from the macroscopic one due to the gaps between the atom and the surrounding medium atoms , which are ignored on a coarse - grained macroscopic scale @xcite . accounting for the difference between the two fields hence requires a correction  the local - field correction . classical calculations of local - field effects , which typically have their origin in near dipole - dipole interactions , can be found in textbooks ( see , e.g. , ref . @xcite ) . in quantum theory , investigations of local - field effects are often related to the problem of the spontaneous decay of an excited guest atom embedded in a ( dielectric ) host . local - field effects in spontaneous decay have been studied , e.g. , for crystals @xcite and disordered dielectrics @xcite on the basis on microscopic models for coupled atomic dipoles .    in macroscopic descriptions , local - field effects are frequently taken into account by regarding the guest atom as being enclosed in a virtual @xcite or real ( spherical ) cavity @xcite surrounded by the medium , with the cavity size being small compared to the relevant transition wavelength . the cavity in the former model is virtual in the sense that it does not perturb the macroscopic field . it is filled by the atoms comprising the medium , which produce no net effect at the central position in important special cases such as cubic or random structures . in the latter model , the field is modified by the presence of the cavity , which is an empty region containing only the guest atom . microscopic models often tend to agree with the virtual - cavity results @xcite , while many recent experiments on spontaneous emission in dielectrics support the real - cavity model @xcite . it has been presumed that while the virtual - cavity model applies to interstitial atoms , the real - cavity model is specific to substitutional atoms , and that the case of substitutional atoms occurs prevalently for impurity atoms in disordered dielectrics  @xcite .    since in the approaches to the local - field effects , the host medium has typically been assumed to be a bulk medium that extends homogeneously to infinity , the question of the effect of the size and shape of the host medium on the local field has arisen . in a macroscopic approach , the spontaneous - decay rate of an excited atom in some free - space region can be given in terms of the imaginary part of the green tensor of the macroscopic maxwell equations , which characterizes the ( macroscopic ) environment of the atom . this relation in principle allows for including local - field corrections for atoms embedded in arbitrary material configurations by assuming a real ( spherical ) cavity surrounding the atom and calculating the corresponding green tensor . using the real - cavity model , toma @xcite studied the local - field correction to the spontaneous - decay rate of an excited atom , which is located at the center of a dispersing and absorbing dielectric sphere . reformulating the result by representing it in a form , which does not explicitly refer to the highly symmetric system considered , he made the conjecture that it may also remain valid beyond the specific example and hence also apply to other locations of the atom and other shapes of the host body .    based on a numerical computation of the respective green tensors , rahmani and bryant @xcite considered the case of an atom at an arbitrary location within a dielectric sphere or a dielectric cube . comparing their results for the dielectric sphere including the local - field correction with earlier results disregarding the local - field correction @xcite , they suggested a rate formula , which in the case of weakly absorbing material corresponds to toma s formula . however , since their approach relies heavily on numerical calculations , it can not produce explicit expressions for the quantities that are related to the local - field correction . the exact analytical evaluations of the green tensors of realistic systems which have finite sizes and include a cavity can be very cumbersome . in this paper , we present an attempt to overcome this difficulty by writing the green tensor as a born series in terms of the susceptibility , where in many situations one can restrict oneself to several leading - order terms . in particular , we show that to linear order the spontaneous - decay rate of a guest atom in a dielectric host body can be separated into a term representing the local - field correction to the decay rate in free space and a term related to the scattering green tensor of the body without the atom  a result which exactly corresponds to toma s conjecture mentioned above . furthermore , we show that for atoms which are situated at the interior of a macroscopic body toma s conjecture remains valid beyond the linear order . we illustrate the theory by discussing in detail the case of an atom embedded in a spherical dielectric body .    the paper is organized as follows . the basic equations for the spontaneous - decay rate and the born expansion of the green tensor determining the rate are given in sec . [ sec_be ] . they are used in sec . [ subsec_rc ] to study the problem of the local - field corrected decay rate within the frame of the real - cavity model , and a proof of toma s formula is given . the examples of an atom embedded in a bulk dielectric medium or in a dielectric sphere are examined in sec . [ sec_ex ] , followed by a summary ( sec .  [ sum ] ) . consider an excited two - level electric - dipole emitter , henceforth referred to as an atom , which is positioned at @xmath0 and surrounded by dispersing and absorbing dielectric bodies . the spontaneous - decay rate can be given in the form of @xcite @xmath1 where @xmath2 and @xmath3 are the ( real ) dipole matrix element and ( shifted ) frequency of the relevant atomic transition , respectively , and @xmath4 @xmath5 @xmath6 . the green tensor of the bodies , @xmath7 , satisfies the equation @xmath8 ( @xmath9 , unit tensor ) together with the boundary condition @xmath10 where @xmath11 is the frequency- and space - dependent complex permittivity which satisfies the kramers  kronig relations . note that satisfaction of the boundary condition  ( [ e4b ] ) is ensured by assuming @xmath12 @xmath13 @xmath14 . equation ( [ e3 ] ) always applies when the atom is placed in some free - space region . separating the green tensor into bulk and scattering parts @xmath15 and @xmath16 , respectively , @xmath17 and taking into account that in the case where the bulk part refers to free space , the relation @xmath18 holds ( see , e.g. , ref . @xcite ) , we may rewrite eq . ( [ e3 ] ) as @xmath19 where @xmath20 is the spontaneous - decay rate in free space .    if the atom is embedded in a body , application of eq . ( [ e3 ] ) requires special care in two respects . firstly , the coincidence limit of the bulk part of the green tensor diverges when the permittivity of the body is complex , as is the case in general . only if material absorption can be neglected so that the permittivity can be regarded as being real , @xmath21 @xmath22 @xmath23 , this limit exists , @xmath24 ( see , e.g. , ref . @xcite ) , leading to @xmath25 secondly , the green tensor of the macroscopic maxwell equations does not account for the fact that the local field felt by the atom is different from the macroscopic one in general . that is , even if absorption is neglected , the rate formula ( [ eq.10 ] ) is not complete because it lacks the local - field corrections .      to calculate the ( scattering part of the ) green tensor for an arbitrary arrangement of dielectric bodies , it may be helpful to use an appropriate born expansion . decomposing the permittivity as @xmath26 and assuming that the solution @xmath27 to the equation @xmath28 is known [ where @xmath29 is defined as in eq . ( [ e5 ] ) with @xmath30 instead of @xmath31 , the green tensor can be written in the form of a born series , @xmath32\nonumber\\ & \quad\times   \overline{\bm{g}}(\mathbf{r},\mathbf{s}_1,\omega)\!\cdot\ !   \overline{\bm{g}}(\mathbf{s}_1,\mathbf{s}_2,\omega)\!\cdots\ !   \overline{\bm{g}}(\mathbf{s}_k,\mathbf{r}',\omega),\end{aligned}\ ] ] as can be verified using the relationships @xmath33,\\ \label{e11 } \hat{h}\,\delta_k\bm{g}(\mathbf{r},\mathbf{r}',\omega ) = & \;\frac{\omega^2}{c^2}\chi(\mathbf{r},\omega )   \bigl[\delta_{k-1}{\bm{g}}(\mathbf{r},\mathbf{r}',\omega )   \nonumber\\ & \ ; -\delta_k\bm{g}(\mathbf{r},\mathbf{r}',\omega)\bigr ]   \qquad\mbox{for } k>1.\end{aligned}\ ] ]    the expansion ( [ e8 ] ) of the green tensor is valid for arbitrarily spatially varying @xmath34 and @xmath35 . obviously , it may be very useful when @xmath36 can be regarded as being a ( small ) perturbation to @xmath30 such that one makes only a small error by disregarding the higher - order terms . in particular , this is the case if the bodies are weakly polarizable , as we shall assume in the following .      for weakly polarizable bodies , it is natural to regard the susceptibilities of the bodies as a small perturbation to the free - space permittivity @xmath34 @xmath5 @xmath37 @xmath38 @xmath39 , i.e. , @xmath40 @xmath41 @xmath42 . this means that we focus on frequencies that are sufficiently far from a resonance frequency of the dielectric material . note that the small ( positive ) imaginary part of @xmath30 ensures that @xmath43 fulfills the boundary condition according to eq . ( [ e4b ] ) so that the spatial integrals in eq . ( [ e9 ] ) converge . we have ( see , e.g. , ref . @xcite ) @xmath44 where @xmath45 @xmath46 @xmath5 @xmath47 @xmath48 @xmath49 , @xmath50 @xmath5 @xmath51 , @xmath52 @xmath5 @xmath53 , @xmath54 @xmath5 @xmath55 .    separating the green tensor into bulk and scattering parts in accordance with eq . ( [ e14.2 ] ) , assuming the atom to be located in a free - space region such that @xmath56 [ cf . ( [ eq.6 ] ) ] , and applying eq . ( [ e8 ] ) , we may represent the scattering part of the green tensor in the rate formula  ( [ e14.4 ] ) as @xmath57 for @xmath58 and @xmath59 belonging to a free - space region , substitution of eq . ( [ e13 ] ) into eq . ( [ e9 ] ) yields the first- and second - order terms in the born expansion , @xmath60 and @xmath61 , respectively , as follows : @xmath62   e^{i\sqrt{1+i\eta}\,(q+q')}\end{aligned}\ ] ] [ @xmath46 @xmath5 @xmath47 @xmath48 @xmath63 , @xmath54 @xmath5 @xmath55 , @xmath64 @xmath5 @xmath65 , @xmath66 @xmath5 @xmath67 ; @xmath68 @xmath5 @xmath63 @xmath48 @xmath49 , @xmath69 @xmath5 @xmath70 , @xmath71 @xmath5 @xmath72 , @xmath73 @xmath5 @xmath74 , @xmath75   e^{i\sqrt{1+i\eta}\,(q+q')}\nonumber\\ & + \frac{k^7}{64\pi^3}\!\int\!\mathrm{d}^3s_1\ ,   \chi(\mathbf{s}_1,\omega)\ ! \int\!\mathrm{d}^3s_2\,\chi(\mathbf{s}_2,\omega )   e^{i\sqrt{1+i\eta}\,(q_1+q_{12}+q_2 ) }   \nonumber\\ & \quad\times\bigl[a_1a_{12}a_2\bm{i }   -a_1a_{12}b_2\tilde{\mathbf{u}}_2\tilde{\mathbf{u}}_2   -a_1b_{12}a_2\tilde{\mathbf{u}}_{12}\tilde{\mathbf{u}}_{12 }   \nonumber\\ & \qquad -b_1a_{12}a_2\tilde{\mathbf{u}}_1\tilde{\bm{u}}_1   + a_1b_{12}b_2(\tilde{\mathbf{u}}_{12 }   \!\cdot\!\tilde{\mathbf{u}}_2 )   \tilde{\mathbf{u}}_{12}\tilde{\mathbf{u}}_2\nonumber\\ & \qquad + b_1a_{12}b_2(\tilde{\mathbf{u}}_1   \!\cdot\!\tilde{\mathbf{u}}_2 )   \tilde{\mathbf{u}}_1\tilde{\mathbf{u}}_2   + b_1b_{12}a_2(\tilde{\mathbf{u}}_1   \!\cdot\!\tilde{\mathbf{u}}_{12 } )   \tilde{\mathbf{u}}_1\tilde{\mathbf{u}}_{12 }   \nonumber\\ & \qquad -b_1b_{12}b_2(\tilde{\mathbf{u}}_1   \!\cdot\!\tilde{\mathbf{u}}_{12 } )   ( \tilde{\mathbf{u}}_{12 }   \!\cdot\!\tilde{\mathbf{u}}_2 )   \tilde{\mathbf{u}}_1\tilde{\mathbf{u}}_2\bigr]\end{aligned}\ ] ] [ @xmath76 @xmath5 @xmath47 @xmath48 @xmath77 , @xmath78 @xmath5 @xmath77 @xmath48 @xmath79 , @xmath80 @xmath5 @xmath79 @xmath48 @xmath49 ; @xmath81 @xmath5 @xmath82 , @xmath83 @xmath5 @xmath84 , @xmath85 @xmath5 @xmath86 for @xmath87 @xmath88 @xmath89 . consider an excited two - level atom embedded in an arbitrary dispersing and absorbing dielectric body characterized by @xmath11 . in order to find the spontaneous - decay rate including local - field corrections , we employ the real - cavity model , that is to say , we assume that the atom is located at the center of an empty - space region of the form of a spherical cavity of radius @xmath90 , as sketched in fig . [ schematic](a ) , where we have denoted the cavity volume by @xmath91 , the volume of the host body without the cavity by @xmath92 ( the overall volume of the host body being @xmath93 @xmath5 @xmath94 ) , and all the remaining space by @xmath95 . hence , the permittivity of the system changes to @xmath96     into ( b ) @xmath97 and ( c ) @xmath98 ; the dot indicates the position of the guest atom . ] the cavity radius @xmath90 is a model parameter representing an average distance from the atom to the nearest neighboring atoms constituting the host body ; it has to be determined from other ( preferably microscopic ) calculations or experiments . note that the real - cavity model is applicable provided that the unperturbed host body is homogeneous ( and isotropic ) in the region where the guest atom is implanted , @xmath99      restricting our attention to the first - order term in the born expansion ( [ e14.6 ] ) , the ( scattering ) green tensor corresponding to @xmath100 can be calculated from eq . ( [ e14.1 ] ) , where after some manipulations one obtains @xmath101\nonumber\\ & \qquad = \bm{g}_\mathrm{c}^{(1 ) }   ( \mathbf{r}_\mathrm{a},\mathbf{r}_\mathrm{a},\omega_\mathrm{a } )   + \bm{g}_\mathrm{b}^{(1 ) }   ( \mathbf{r}_\mathrm{a},\mathbf{r}_\mathrm{a},\omega_\mathrm{a})\end{aligned}\ ] ] [ @xmath46 @xmath5 @xmath102 @xmath48 @xmath63 , @xmath54 @xmath5 @xmath103 , where , again to linear order in the susceptibility , @xmath104\end{gathered}\ ] ] is the scattering green tensor of a spherical cavity embedded in a bulk medium of susceptibility @xmath105 @xmath106 @xmath107 [ cf . [ schematic](b ) ] , and @xmath108   \nonumber\\ & -\frac{k_\mathrm{a}^4\chi(\omega_\mathrm{a})}{16\pi^2 }   \int_{b\cup v }   \mathrm{d}^3s\ ,   e^{(2i-\eta)q}\bigl[a^2\bm{i}+(b^2 - 2ab )   \tilde{\mathbf{u}}\tilde{\mathbf{u}}\bigr]\end{aligned}\ ] ] is the scattering part of the green tensor of the host body without the cavity [ cf . [ schematic](c ) ] . the result can be verified by applying eq . ( [ e14.2 ] ) [ together with eqs . ( [ e8 ] ) and ( [ e14.1 ] ) ] to the two systems mentioned and using eq . ( [ e15.2 ] ) . equations ( [ e15.3])([e15.5 ] ) show that @xmath16 in the rate formula ( [ e14.4 ] ) can be written as the sum of two terms , where the first term , @xmath97 , only depends on the cavity radius and the local permittivity of the host body at the position of the atom , whereas the second term , @xmath98 , is determined by the properties of the host body in the absence of the atom . hence , the term in the decay rate which is proportional to @xmath109 can be regarded as being the local - field correction to the uncorrected term proportional to @xmath110 . the fact that the local - field correction additively enters the rate formula is due to the linear expansion in @xmath36 . inspection of the second - order term in the born expansion , eq . ( [ e16 ] ) , indicates that in general , terms depending on both the cavity and the ( unperturbed ) host body will appear which may lead to a breakdown of the additivity . however , as shown in sec . [ subsec_con ] , there are situations where a generalization beyond the linear approximation is possible .    equations ( [ e15.4 ] ) and ( [ e15.5 ] ) can be further evaluated by introducing a spherical coordinate system whose origin coincides with the position of the atom , @xmath111 where @xmath112 and @xmath113 , respectively , refer to the inner and outer boundary areas of the integration volumes sketched in figs . [ schematic](b ) and [ schematic](c ) . performing the radial integral in eq . ( [ e15.4 ] ) , we find , on recalling eq . ( [ e14 ] ) , @xmath114 ( @xmath115 @xmath5 @xmath116 ) with @xmath117   = -\frac{k_\mathrm{a}\chi(\omega_\mathrm{a})}{16\pi^2 }   \int_0^{2\pi}\mathrm{d}\phi\int_0^\pi\mathrm{d}\theta\sin\theta   \biggl\{e^{(2i-\eta)q } \nonumber\\ & \times\biggl[\biggl(\frac{1}{3q^3}-\frac{2i}{3q^2 }   -\frac{5}{3q}+\frac{i}{2}\biggr)\bm{i }   + \biggl(\frac{1}{q^3}-\frac{2i}{q^2}+\frac{3}{q }   -\frac{i}{2}\biggr )   \tilde{\mathbf{s}}\tilde{\mathbf{s}}\biggr ]   \nonumber\\ & \hspace{9ex } + 4i\,\mathrm{ei}[(2i\!-\!\eta)q]\biggl(\frac{1}{3}\bm{i }   -\tilde{\mathbf{s}}\tilde{\mathbf{s}}\biggr)\biggr\}\end{aligned}\ ] ] [ @xmath118 @xmath5 @xmath119 ; @xmath120 @xmath5 @xmath121 ; @xmath122 , exponential integral ] . using the fact that @xmath115 is independent of @xmath123 and @xmath124 as well as the relation @xmath125 from eq . ( [ e19 ] ) we obtain @xmath126 in particular , when the radius of the cavity is much smaller than the atomic transition wavelength , @xmath127 then eq .  ( [ e19.1 ] ) reduces to @xmath128\bm{i}.\ ] ] substitution of this result together with eq . ( [ e19b1 ] ) into eq . ( [ e15.3 ] ) reveals that to linear order in @xmath36 , @xmath129\bm{i }   \\ + \bm{g}_\mathrm{b}^{(1 ) }   ( \mathbf{r}_\mathrm{a},\mathbf{r}_\mathrm{a},\omega_\mathrm{a } )   + o(k_\mathrm{a}r_\mathrm{c}).\end{gathered}\ ] ]    for a homogeneous host body which is star - shaped w.r.t . the position of the atom ( i.e. , every point on the outer boundary area can be connected to the atomic position by a straight line that lies entirely within the body ) , @xmath98 as given by eq . ( [ e15.5 ] ) can be evaluated in a similar manner . using again spherical coordinates and recalling eq . ( [ e14 ] ) , we may evaluate the radial integrals in eq . ( [ e15.5 ] ) to obtain @xmath130   -\bm{f}_0(0)\nonumber\\ & \quad -\bigl[\lim_{\eta\to 0_+}\lim_{q\to\infty }   \bm{f}_\eta(q)-\bm{f}_0(0)\bigr ]   = \lim_{\eta\to 0_+}\bm{f}_\eta[q_\mathrm{o}(\phi,\theta)]\end{aligned}\ ] ] [ @xmath131 @xmath5 @xmath132 @xmath5 @xmath133 , where henceforth @xmath113 denotes the outer boundaries of the host body .    substituting eq . ( [ e18 ] ) into eq . ( [ e14.4 ] ) , we find that to linear order in @xmath36 the local - field corrected spontaneous - decay rate of an atom within a body can be given as follows : @xmath134 where @xmath135 is the free - space decay rate as given in eq . ( [ e14.5 ] ) , @xmath136 and @xmath137 equations ( [ e24])([e24.2 ] ) are valid for an atom embedded in a weakly polarizable , dispersing and absorbing body of arbitrary size and shape . note that only @xmath138 takes into account the effect of local - field correction , while @xmath139 , which is simply determined by the scattering part of the green tensor of the host body without the cavity , reflects the uncorrected influence of the size and shape of the body on the decay rate . in particular , the first two terms in eq . ( [ e24.1 ] ) obviously result from irreversible energy transfer from the atom to the surrounding matter . if these two terms are dominant over the last one , effectively no radiation is emitted . we conclude this section by making two remarks concerning the necessary conditions under which eq . ( [ e24 ] ) provides a good approximation to the spontaneous - decay rate . ( i ) the permittivity often appears in the green tensor as a common factor @xmath140 [ cf . ( [ e46 ] ) in sec .  [ subsec_oc ] ] . to be on the conservative side , one should then take @xmath141 [ with @xmath142 @xmath5 @xmath143 rather than @xmath36 as the small parameter . ( ii ) equation ( [ e24.1 ] ) already contains terms of the order of @xmath144 . thus @xmath144 , not @xmath145 , should be much smaller than unity for the linear approximation to yield a good estimation of the decay rate . the cavity radius @xmath90 represents the average distance between the atom and the constituents of the body . for , say , @xmath146 @xmath147 @xmath148 , it is sufficient to require that , so that .      toma @xcite has found that in the special case of an atom situated at the center of a homogenous , dielectric sphere , the real - cavity model ( fig .  [ sphere ] in the case where ) leads to        the relation @xmath149\biggr\}\bm{i}\nonumber\\ & \quad+\biggl(\frac{3\varepsilon}{2\varepsilon+1}\biggr)^2   \bm{g}_\mathrm{b}^{(1 ) }   ( \mathbf{r}_\mathrm{a},\mathbf{r}_\mathrm{a},\omega_\mathrm{a } )   + o(k_\mathrm{a}r_\mathrm{c}),\end{aligned}\ ] ] @xmath150 @xmath5 @xmath151 , where @xmath152 is the scattering green tensor of the homogeneous dielectric sphere , whose geometry dependence is  for @xmath58 @xmath5 @xmath49 @xmath5 @xmath102 at the sphere center  entirely given by its dependence on the sphere radius @xmath153 . as already mentioned in sec . [ sec_intro ] , he made the conjecture that this relation might be more generally valid for ( i ) bodies of arbitrary sizes and shapes and ( ii ) arbitrary positions of the atom , provided that the atom is not on the body surface . let us study the validity of this conjecture in more detail . it can easily be seen that to linear order in @xmath36 , eq . ( [ e21 ] ) obviously reduces to eq . ( [ e18 ] ) , so the results of sec .  [ linear ] show that in this order the conjecture is true even under the more general conditions of inhomogeneous host bodies , provided that the requirement ( [ e15.2 ] ) for the applicability of the real - cavity model is valid . moreover , we will demonstrate in the following that eq . ( [ e21 ] ) remains valid beyond the linear order in @xmath36 , provided that the respective host body can be regarded as being homogeneous in the vicinity of the atom , i.e. , @xmath154 with @xmath155 being some small positive number . we begin with the case of a homogeneous body , in which case satisfaction of the condition ( [ e15.2b ] ) simply ensures that the entire cavity lies inside the body . we first recall that @xmath156 determines the electric field reaching the point @xmath0 , in which it has originated , after being scattered at the surfaces of inhomogeneity [ cf . [ schematic](a ) ] . at the cavity surface the electric field can be separated into two parts , namely , one part that is ( multiply ) reflected at this surface and eventually returns to the point @xmath0 , and one part that is eventually transmitted to the exterior of the cavity . the contribution due to the reflected part is obviously given by @xmath157 [ cf . [ schematic](b ) ] , which reads @xcite @xmath158\biggr\}\bm{i}\nonumber\\ & \quad+o(k_\mathrm{a}r_\mathrm{c}).\end{aligned}\ ] ] in an infinitely extended body [ cf .  fig . [ schematic](b ) ] , the transmitted part at a point @xmath58 outside the cavity is determined by the green tensor @xcite @xmath159 [ @xmath46 @xmath5 @xmath47 @xmath48 @xmath102 , @xmath54 @xmath5 @xmath160 , @xmath161 @xmath5 @xmath162 , and @xmath64 and @xmath66 according to eq . ( [ e14 ] ) ] . the coefficient @xmath163 ( which corresponds to @xmath164 in ref . @xcite ) is given by @xmath165 '   -[z_0j_1(z_0)]'h_1^{(1)}(z_0)}{j_1(z_0)[z_1h_1^{(1)}(z_1 ) ] '   -\varepsilon[z_0j_1(z_0)]'h_1^{(1)}(z_1)}\,,\ ] ] where @xmath166 @xmath5 @xmath116 , @xmath167 @xmath5 @xmath168 ( the primes denote derivatives w.r.t . @xmath166 and @xmath167 ) , and @xmath169 and @xmath170 are the spherical bessel and hankel functions , respectively . inspection of eq . ( [ e21e ] ) shows that @xmath171 substitution of this result into eq . ( [ e21d ] ) leads to @xmath172 where @xmath173 is the green tensor of the infinite body without the cavity [ according to eqs . ( [ e13 ] ) and ( [ e14 ] ) with @xmath174 . in other words , the electric field transmitted through the cavity surface to a point outside the cavity is equal to the field that would be transmitted to the same point in the absence of the cavity , multiplied by a global factor . by means of the general symmetry property ( see , e.g. , ref . @xcite ) @xmath175 we also have @xmath176 i.e. , the electric field transmitted through the cavity surface from an arbitrary point outside the cavity is also equal to the corresponding result in the absence of the cavity , multiplied by the same factor . for a finite body , the electric field will be ( multiply ) reflected from the body s outer surface , eventually giving rise to a field at @xmath0 . without the cavity , processes of this kind are taken into account by replacing the infinite - body green tensor @xmath177 with its finite - body counterpart @xmath178 . combining this observation with eqs . ( [ e21 g ] ) and ( [ e21i ] ) , we conclude , on recalling the linearity of maxwell s equations , that the electric field which is transmitted though the cavity surface , scattered at the outer body surface , and finally retransmitted into the cavity is given by @xmath179 so combining eqs . ( [ e21c ] ) and ( [ e21j ] ) , we arrive at eq . ( [ e21 ] ) .    in this derivation , we have disregarded processes involving scattering of the field at the cavity surface from the outside . processes of this kind can indeed be neglected , because their contributions are of orders higher than @xmath180 . to see this , note that for a cavity in bulk material [ fig . [ schematic](b ) ] the field at a point @xmath58 outside the cavity originating in a point @xmath59 outside the cavity via scattering at the cavity surface , is characterized by @xcite @xmath181,\end{gathered}\ ] ] where @xmath182 '   -[z_0j_m(z_0)]'j_m(z_1)}{j_m(z_0)[z_1h_m^{(1)}(z_1 ) ] '   -[z_0j_m(z_0)]'h_m^{(1)}(z_1)}\,,\\ & b_m^n=-\frac{j_m(z_0)[z_1j_m(z_1 ) ] '   -\varepsilon[z_0j_m(z_0)]'j_m(z_1)}{j_m(z_0)[z_1h_m^{(1)}(z_1 ) ] '   -\varepsilon[z_0j_m(z_0)]'h_m^{(1)}(z_1)},\end{aligned}\ ] ] and the tensors @xmath183 and @xmath184 do not depend on @xmath90 . using the relations @xcite @xmath185,\\ \label{e21n } & h_m^{(1)}(z ) = -\frac{i(2m-1)!!}{z^{m+1 } }   \bigl[1+o(z^2)\bigr],\end{aligned}\ ] ] one can easily show that @xmath186,\quad b_m^n = o\bigl[(k_\mathrm{a}r_\mathrm{c})^{2m+1}\bigr],\ ] ] so eq . ( [ e21k ] ) leads to @xmath187.\ ] ] processes involving reflections of the field at the cavity surface from the outside hence start to contribute in third order of @xmath146 and therefore do not need to be included in eq . ( [ e21 ] ) . so far we have demonstrated the validity of the relation ( [ e21 ] ) for homogeneous dielectric bodies of arbitrary shapes provided that the atom is situated at the interior of the body , such that the condition ( [ e15.2b ] ) is satisfied . note that this condition practically coincides with the condition ( [ e15.2 ] ) for the applicability of the real - cavity model . the arguments given above can be extended to inhomogeneous bodies . if the condition ( [ e15.2b ] ) is satisfied , one can divide such a body into a more or less small homogeneous part containing the cavity plus an inhomogeneous rest . equations  ( [ e21 g ] ) , ( [ e21i ] ) , and ( [ e21p ] ) then still describe the propagation of the electric field inside the homogeneous part of the body , and the effect of the inhomogeneous part can be taken into account by the scattering at the ( fictitious ) surface dividing the two parts . consequently , we are again left with eq . ( [ e21 ] ) .    substituting eq . ( [ e21 ] ) into eq . ( [ e14.4 ] ) and recalling eq . ( [ e14.5 ] ) , we can again cast the local - field corrected spontaneous - decay rate in the form of of eq . ( [ e24 ] ) , where now @xmath188\biggr\},\end{aligned}\ ] ] and @xmath189 \!\cdot\mathbf{d}_\mathrm{a}\ ] ] [ @xmath150 @xmath5 @xmath190 . recall that @xmath98 is the ( unperturbed ) scattering part of the green tensor of the host body without the cavity . needless to say that to linear order in  @xmath36 , eqs . ( [ e22.2 ] ) and ( [ e22.2b ] ) reduce to eqs . ( [ e24.1 ] ) and ( [ e24.2 ] ) , respectively . in particular in the case of weakly absorbing matter , it may be sufficient to retain in eqs . ( [ e22.2 ] ) and ( [ e22.2b ] ) only terms to linear order in @xmath191 . from eqs . ( [ e24 ] ) , ( [ e22.2 ] ) , and ( [ e22.2b ] ) it then follows that when @xmath192 then @xmath193 can be given in the form @xmath194 where @xmath195 and @xmath196 is the uncorrected decay rate as given by eq . ( [ eq.10 ] ) with @xmath156 @xmath106 @xmath197 . rahmani and bryant @xcite concluded from a numerical computation of the green tensor of a dielectric sphere that contains a ( small ) empty sphere at an arbitrary position inside the material that the local - field corrected spontaneous decay rate has the form @xmath193 @xmath5 @xmath198 @xmath38 @xmath199 , where the shift term @xmath199 is due to absorption . comparing with eq . ( [ eq.61 ] ) [ together with eq . ( [ eq.62 ] ) ] , we see that this is indeed the case when the effect of absorption is sufficiently weak and in particular , the inequality ( [ eq.60 ] ) holds . however , from eqs . ( [ e24 ] ) , ( [ e22.2 ] ) , and ( [ e22.2b ] ) it is clearly seen that in general , @xmath193 can not be given in the form assumed in ref . note that already the analytical solution for the special case considered in ref . @xcite implies that this form can not be true in general . in the case of bulk material , we have @xmath98 @xmath5 @xmath14 , so eq .  ( [ e24 ] ) simplifies to @xmath200 where @xmath201 is given by eq . ( [ e22.2 ] ) @xcite , which to linear order in @xmath36 reduces to eq . ( [ e24.1 ] ) . in particular , when material absorption can be neglected , @xmath150 @xmath22 @xmath202 , we simply have  @xcite @xmath203 [ which follows directly from eq . ( [ eq.61 ] ) with @xmath196 @xmath5 @xmath204 , cf . ( [ eq.10 ] ) ] . note that in this case the virtual - cavity model leads to @xmath205 ( see , e.g. , ref . it is not difficult to see that to linear order in @xmath36 both eq .  ( [ e0 ] ) and eq .  ( [ e0a ] ) lead to @xmath206\gamma_0,\ ] ] i.e. , the first - order theory for the spontaneous decay rate does not distinguish between the virtual- and the real - cavity model @xcite , provided that absorption can be disregarded . let us consider an atom in a homogeneous dielectric sphere as sketched in fig . [ sphere ] and first calculate the spontaneous - decay rate to linear order in @xmath36 . in particular , we need to calculate @xmath139 as given by eq . ( [ e24.2 ] ) with @xmath98 from eq . ( [ e19b2 ] ) together with eq . ( [ e19 ] ) . noting that @xmath207\ ] ] ( cf . [ sphere ] ) is independent of @xmath124 , performing the @xmath124-integral , and using the fact that @xmath208 @xmath5 @xmath209 for a radially ( @xmath210 ) oriented dipole and @xmath208 @xmath5 @xmath211 for a tangentially ( @xmath212 ) oriented dipole , we find @xmath213,\end{aligned}\ ] ] where @xmath214   e^{(2i-\eta)q_\mathrm{o}}\nonumber\\ & + 4i\biggl(\frac{1}{3}-z\biggr )   \mathrm{ei}[(2i-\eta)q_\mathrm{o}],\\[1ex ] \label{e28 } f_\eta^\parallel(q_\mathrm{o},z )   = & \ ; f_\eta^\perp[q_\mathrm{o},(1-z)/2],\end{aligned}\ ] ] and @xmath131 according to eq . ( [ e25 ] ) , with @xmath215 being replaced by  @xmath216 .    in order to exactly calculate @xmath139 [ eq .  ( [ e22.2b ] ) ] , we make use of the exact green tensor for a dielectric sphere as given in ref . @xcite , leading to @xmath217 ^ 2 \biggr\},\end{gathered}\ ] ] @xmath218 }   { ( nk_\mathrm{a}l_\mathrm{a})\mathrm{d}l_\mathrm{a}}\biggr)^2   \biggr]\biggr\},\end{gathered}\ ] ] where @xmath219 '   -[z_1h_m^{(1)}(z_1)]'h_m^{(1)}(z_0 ) }   { \varepsilon j_m^{(1)}(z_1)[z_0h_m^{(1)}(z_0 ) ] '   -[z_1j_m^{(1)}(z_1)]'h_m^{(1)}(z_0)}\,,\\ \label{e37 } & c^m_m=-\frac{h_m^{(1)}(z_1)[z_0h_m^{(1)}(z_0 ) ] '   -[z_1h_m^{(1)}(z_1)]'h_m^{(1)}(z_0 ) }   { j_m^{(1)}(z_1)[z_0h_m^{(1)}(z_0 ) ] '   -[z_1j_m^{(1)}(z_1)]'h_m^{(1)}(z_0)}\,,\end{aligned}\ ] ] with @xmath166 @xmath106 @xmath220 , @xmath167 @xmath106 @xmath221 .      to compare the exact decay rate [ eq . ( [ e24 ] ) together with eqs . ( [ e14.5 ] ) , and ( [ e22.2 ] ) as well as eqs . ( [ e34 ] ) and/or ( [ e35 ] ) ] with the decay rate obtained in the linear born approximation [ eq . ( [ e24 ] ) together with eqs . ( [ e14.5 ] ) , ( [ e24.1 ] ) , and ( [ e26 ] ) ] , let us consider , for simplicity , the case where the atom is positioned at the center of the sphere . setting @xmath222 @xmath5 @xmath14 and hence , in eqs . ( [ e26])([e28 ] ) , and carrying out the @xmath216 integral in eq . ( [ e26 ] ) , we derive @xmath223 with @xmath224e^{(2i-\eta)k_\mathrm{a}r}\biggr\},\end{gathered}\ ] ] giving @xmath225e^{2ik_\mathrm{a}r}\biggr\}\end{gathered}\ ] ] for any finite radius @xmath153 , and approaching zero in the limit @xmath153 @xmath226 @xmath227 ( which has to be taken before performing the limit @xmath228 @xmath226 @xmath14 ) . when material absorption is small enough such that @xmath229 @xmath41 @xmath230 , and the radius of the sphere is large @xmath231 @xmath232 @xmath42 , then eq .  ( [ e24 ] ) together with eqs . ( [ e24.1 ] ) and ( [ e44.1 ] ) reduces to [ @xmath105 @xmath22 @xmath233 @xmath234\gamma_0   -\frac{\chi(\omega_\mathrm{a})}{2}\cos(2k_\mathrm{a}r)\gamma_0.\ ] ] note that eq . ( [ e2 ] ) differs from eq . ( [ e1 ] ) in the second term , which reflects the finite size of the sphere .    when the atom is on center , only the terms in eqs . ( [ e34 ] ) and ( [ e35 ] ) contribute to the exact @xmath235 , cf . ( [ e21 m ] ) and ( [ e21n ] ) , and hence we find @xmath236.\ ] ] combining eqs . ( [ e24 ] ) , ( [ e22.2 ] ) , and ( [ e45 ] ) , we arrive at @xmath237 which can be shown to agree with eq . ( [ e44.1 ] ) to linear order in @xmath36 . in a more sophisticated linearization of eq . ( [ eq.80 ] ) with respect to @xmath36 , it may be advantageous to leave the @xmath36-dependence in the exponentials appearing in the coefficient @xmath238 unchanged . in particular , for sufficiently small material absorption , @xmath229 @xmath41 @xmath230 , and for a sufficiently large sphere , @xmath231 @xmath232 @xmath42 , this kind of linear approximation leads to [ @xmath105 @xmath22 @xmath233 @xmath239\gamma_0\\   -\frac{\chi(\omega_\mathrm{a})}{2}\ ,   \cos\bigl[2\mathrm{re}(n)k_\mathrm{a}r\bigr ]   e^{-2\mathrm{im}(n)k_\mathrm{a}r}\gamma_0,\end{gathered}\ ] ] recall eqs . ( [ e2.1 ] ) and ( [ e47 ] ) . @xmath5 @xmath240 and ( b ) @xmath150 @xmath5 @xmath241 as a function of the sphere radius @xmath153 . the radius of the real ( empty - space ) cavity is @xmath90 @xmath5 @xmath242 . for comparison , the decay rate in a bulk medium of the same permittivity ( dotted lines ) is also plotted . ] comparison of eq . ( [ e46 ] ) with eq . ( [ e2 ] ) indicates that with increasing size of the sphere , the validity of the linear born approximation to the decay rate becomes less satisfying in two respects . ( i ) the oscillating term is not damped , and ( ii ) there is an accumulated error in the phase . the former effect is insignificant as long as @xmath243 @xmath41 @xmath42 . this restriction still allows for a large range of sphere sizes ; for example , for @xmath145 @xmath22 @xmath244 and @xmath245 @xmath5 @xmath246 @xmath22 @xmath247 , the condition follows . figure  [ cnt ] shows the dependence on the sphere radius of the local - field corrected spontaneous - decay rates according to the exact equation ( [ eq.80 ] ) and the linear born expansion [ eq . ( [ e24 ] ) together with eqs . ( [ e24.1 ] ) and ( [ e44.1 ] ) ] for two values of the permittivity . for the parameters used , we have @xmath144 @xmath5 @xmath248 , so that the first two terms in the curly brackets in eq . ( [ eq.80 ] ) , which arise from absorption , are negligibly small compared to the last one . as expected from eqs . ( [ e2 ] ) and ( [ e46 ] ) , the decay rate oscillates with the sphere radius around the bulk value . for small spheres , the linear born expansion is seen to be in good agreement with the exact result . it is further seen that , as the sphere radius increases , an increasing phase shift develops between the exact rate and the rate in the linear born approximation . a comparison of figs . [ cnt](a ) and ( b ) reveals that this phase shift is larger for larger values @xmath249 . @xmath5 @xmath250 as a function of @xmath145 ( @xmath251 @xmath5 @xmath252 ) , for real - cavity radii @xmath90 @xmath5 @xmath242 and @xmath253 . ] the dependence of the decay rate on the imaginary part of the susceptibility is illustrated in fig . [ ci ] for two values of the real - cavity radius . it can be seen that the disagreement between the curves representing the linear born expansion and the exact result increases with increasing material absorption or decreasing cavity radius , i.e. , with an increase in the combined factor @xmath144 . in particular , if in the case where @xmath90 @xmath5 @xmath242 , @xmath145 changes from @xmath254 to @xmath255 , then @xmath256 changes from @xmath257 to @xmath37 , and the agreement worsens from being very good to being moderately good . for the larger cavity radius @xmath90 @xmath5 @xmath258 , which physically means a more dilute medium , @xmath256 varies from @xmath259 to @xmath260 for the same variation of @xmath145 . throughout this range , the combined factor remains much smaller than one and a good agreement is observed . basing on a numerical computation of @xmath261 as given by eq . ( [ e26 ] ) in the linear born approximation , we have also studied the case when the atom is localized at an arbitrary position inside the sphere .     @xmath5 @xmath240 and radius and ( b ) @xmath153 @xmath5 @xmath262 . the real - cavity radius is @xmath90 @xmath5 @xmath242 . the solid ( dashed ) curves refer to radially ( tangentially ) oriented transition dipole moments . ] figure [ offcnt ] compares the position dependences of the local - field corrected decay rate , as given by eq . ( [ e24 ] ) together with eqs . ( [ e24.1 ] ) and ( [ e26 ] ) , and the uncorrected rate according to eq . ( [ eq.10 ] ) ( which is valid if absorption is neglected ) , with the second term being given by eq . ( [ e26 ] ) , for two sphere radii ( for the uncorrected rate beyond the linear born approximation , see also ref . @xcite ) . from the figure it is seen that when @xmath222 @xmath13 @xmath14 , i.e. , when the atom is not localized at the center of the sphere , radial and tangential dipole orientations have to be distinguished , especially when the sphere radius exceeds the wavelength of the radiation emitted by the atom , such that interference effects begin to play a role . note that the whispering gallery resonances which may give rise to strong enhancement near the rim of the sphere are not manifested here . the existence of these resonances requires larger values of the real part of the permittivity or larger sphere radii . unfortunately , in such cases the linear born expansion is too rough an approximation to the decay rate . from the figure it is further seen that for the parameters used , the local - field correction increases the decay rate ; the amount of increase is determined essentially by the last term in eq . ( [ e24.1 ] ) .    in fig . [ dp ] the dependence of the local - field corrected decay rate on the sphere radius is illustrated . it can be seen that the further the atom is displaced from the sphere center , the smaller the amplitudes of oscillation of the decay rate become . @xmath5 @xmath263 and ( b ) @xmath222 @xmath5 @xmath264 off center of a dielectric sphere of permittivity as a function of the sphere radius in linear born approximation . the cavity radius is @xmath146 @xmath5 @xmath148 . the solid ( dashed ) curves refer to radially ( tangentially ) oriented transition dipole moments . for comparison , the curves for an atom at the center are also shown . ] expressing the spontaneous - decay rate of an excited atom in the presence of dielectric bodies in terms of the scattering part of the associated green tensor of the macroscopic maxwell equations and expanding the green tensor in a born series enables one to systematically approach arbitrary geometries . in particular , in the case of weakly polarizable bodies it may be possible to neglect higher - order terms in the born expansion and approximate the scattering part of the green tensor by the term linear in the susceptibility , and the decay rate accordingly .    using the real - cavity model of the local - field correction , we have applied the theory to the problem of the local - field correction to the spontaneous - decay rate of an excited atom embedded in a dispersing and absorbing dielectric body of arbitrary size and shape . in this way , we have derived a rate formula , which , within the linear born approximation , applies to atoms in arbitrary dielectric bodies , and have given explicit conditions of its validity . to illustrate the results , we have considered the case of an atom at an arbitrary position inside spherical body in more detail . it has surprisingly turned out that the rate formula found in linear born approximation agrees , to linear order in the susceptibility , with a rate formula suggested by toma @xcite from his study of the real - cavity model in the special and analytically solvable case of a spherically symmetric system . we have then shown that this quite general formula indeed applies to atoms in dielectric bodies of arbitrary sizes and shapes , provided that the atoms are not in the very vicinity of the surfaces of the bodies . so it can be shown that the scattering part of the green tensor that enters the basic formula for the decay rate can always be decomposed into a term that only depends on the properties of the local environment of the guest atom and the local - field corrected scattering part of the green tensor of the host body without the guest atom , where the correction simply appears in form of a factor . in particular , the former term can be substantially determined by the absorptive properties of the host body , thereby giving rise to a shift of the decay rate . finally , we note that in the same spirit as in the treatment of the spontaneous decay , the born expansion can also be employed in studying other phenomena of the atom  field interaction in realistic systems whose green tensors are not known or analytically too involved . typical examples may be the casimir - polder interaction of a ground - state atom with an inhomogeneous body @xcite , the resonance fluorescence of an atom near such a body , or the resonant energy transfer between atoms embedded in realistic bodies . we would like to thank m.  s.  toma for discussions . thanks the alexander von humboldt stiftung and the national program for basic research of vietnam for support . this work was supported by the deutsche forschungsgemeinschaft .
Abstract:

Answer: the influence of the size and shape of a dispersing and absorbing dielectric body on the local - field corrected spontaneous - decay of an excited atom embedded in the body is studied on the basis of the real - cavity model . by means of a born expansion of the green tensor of the system it is shown that to linear order in the susceptibility of the body the decay rate exactly follows toma s formula found for the special case of an atom at the center of a homogeneous dielectric sphere [ phys . rev . a * 63 * , 053811 ( 2001 ) ] . it is further shown that for an atom situated at the interior of an arbitrary dielectric body this formula remains valid beyond the linear order . the case of an atom embedded in a weakly polarizable sphere is discussed in detail .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: nowadays , very high redshift galaxies and quasars are being found . however , quasars still offer the opportunity of a ( relatively ) easy construction of high - z samples , due to their high luminosity , quasi - stellar shape ( allowing accurate photometry ) , and to their emission lines ( allowing spectroscopical identification ) . therefore , many large optical quasar surveys are in progress or in project . even if high energy emission is quite a common quasar property , x - ray samples still need the detection of emission lines in the optical domain for the purposes of redshift determination . to improve the efficiency of all quasar studies , one has to increase both the size and the upper redshift limit of the sample , and to achieve a maximum completeness over a redshift range as large as possible . these surveys are essentially based on the photometric preselection of quasar candidates , in order to optimize the efficiency of telescope time . preselection makes the completeness questionable both at low and at high redshift ( see discussion in e.g. miller & mitchell , 1998 ) . even multiple color - color diagrams @xcite still reveal a bias in the range @xmath1 , exactly where a reversing of evolution is thought to occur , making completeness specially important . at low redshift , it has been shown that the usual morphologic prerequisite for _ quasi - stellar _ candidates results in a noticeable deficiency @xcite . this trend may increase the apparent evolution of quasars . the present paper describes a joint method both for distinguishing between quasars and stars / galaxies by their photometry _ and _ for obtaining an estimate of the photometric redshift of the quasar candidates as well . in  [ state ] we present the current state of the art in the optical search for quasars .  [ method ] gives a brief description of our method , which has been first validated with simulations , as shown in  [ simul ] , before being applied to real samples (  [ eis ] ) . the last section of this paper (  [ discuss ] ) contains a summary of our conclusions and a brief discussion . one of the most common problems in optical surveys is the selection criteria one has to apply in order to distinguish between stars and quasars , since both categories of objects have ( in most cases ) a point  like morphology . quasars are spatially resolved only at low redshifts and only when high quality imaging data are available . selection of qso candidates among stellar objects from photometric data is becoming a standard method . the uvx technique ( uv excess ) separates stars from quasars with a redshift @xmath2 . this technique is still currently being used , e.g. in the _ chile - uk quasar survey _ , and the _ aat 2df qso survey_.    to avoid both the redshift limitation @xmath3 and the bias towards blue objects , other equivalent techniques are also being applied ( for example brx ) , with quite a high efficiency in other redshift ranges . many multicolor surveys use an increased number of filters , especially towards the red , with special care paid to high - z quasars : examples are found in boyle et al . ( 1991 ) , kennefick et al . ( 1997 ) , osmer et al . ( 1998 ) and the references therein , jarvis & mac alpine ( 1998 ) and the _ dposs _ ( digital poss ii ) . in all these cases , the efficiency never exceeds 35% ( @xmath4 3 candidates selected per one real quasar ) , depending on magnitude and redshift . for all the selection methods described above , the working space is the n @xmath0 2d - manifold of n color - color planes . the locus of quasars at different redshifts in each of these planes is compared to the locus expected for main sequence stars . for quasars , the allowed multi - color space is defined through simulated template spectra , set to different redshifts . basically all the selection procedures use the positions of candidates in the multi - color space either to compare with the expected locus for quasars or to compute the distance to the main sequence ( newberg & yanny , 1997 ; sloan digital sky survey - fan , 1999 ; krisciunas et al . , 1998 ) . quasar candidates are among the objects at the largest distances with respect to this main sequence `` snake '' . however the selection of candidates with quasar spectra is basically a 1d problem : the fitting between the observed spectral energy distribution ( sed ) and the equivalent one obtained from templates ( quasars / stars ) , all of them computed within the same photometric system . two quasar surveys make use of such a fit to select candidates and to estimate their photometric redshifts . the calar alto deep imaging survey ( _ cadis _ ) uses specialized filters and a multicolor classification algorithm to distinguish between stars , galaxies and quasars , and the identification of high - z quasars in particular is based on the photometric redshift method @xcite . the _ large zenith telescope _ , a 6-meter liquid - mirror telescope , will provide photometric redshifts for more than 1 million of galaxies , thanks to a series of 40 medium - band filters specially designed , and will be able to select quasar candidates quite efficiently . for the purposes of this paper , even a reduced set of broad band photometric data on a given object shall be considered as a very low resolution spectrum ( _ r_@xmath5 1 ! ) , which could be fitted by synthetic templates according to the standard methods commonly used on photometric redshift estimates . assumptions on the morphology of the objects could be avoided in principle , in such a way that both point - like and spatially resolved sources could be examined through the same pipeline . furthermore , it is interesting to have a prior idea on the redshift range of each candidate , in order to perform the spectroscopic follow up in the suitable wavelength range ( visible or ir ) . the software described in the present paper has been developed for the needs of preparing the quasars / agn sample to be assembled in the future virmos survey . it can however be applied to all multicolor catalogues . the forthcoming virmos spectroscopic galaxy survey will consist in 2 surveys : a deep one and a shallow one . the virmos - shallow spectroscopic sample ( statistically 1/4 to 1/3 complete to @xmath6 , no other preselection ) will ensure random completeness , without any color or morphology preselection . the virmos - deep survey ( complete for all objects with @xmath7 ) will provide the opportunity of assembling a quasar sample avoiding the drawbacks of any preselection , therefore basically free of the usual biases ( colors , redshift etc ) . so far , very few such samples are available : three examples come from the cfrs ( i=22.5 , schade et al . , 1995 ) , the faint galaxy redshift survey ( b=24 , glazebrook et al . , 1995 ) and more recently cohen et al . ( 1999 ) to r=24 in the hdf , involving , however , very few quasars . we are examining various possibilities of assembling a third quasar sample in the limits of the virmos - shallow survey , trying to get as much as possible of the spectra of the several thousands of quasars present in this field by selecting photometric candidates . comparing the three virmos samples should provide a unique opportunity of examining the biases introduced by a photometric and/or morphological preselection on the quasar samples . our method used here is closely based on _ hyperz _ , a public photometric redshift code presently under development at the observatoire midi - pyrnes . details on _ hyperz _ will be given in a forthcoming paper ( bolzonella et al . , 2000 ; see also miralles & pell , 1998 ; pell et al . , 1999 ) . it is basically a standard sed fitting procedure : the observed photometric sed of a given object , obtained through @xmath8 filters , is compared to the sed computed for a set of template spectra . the aim is to find the best fit between the observed and the model photometry through a standard @xmath9 minimization procedure :    @xmath10    where @xmath11 and @xmath12 are respectively the observed and the template fluxes in the @xmath13 band , and @xmath14 is the error on the observed flux in this band . fluxes are normalized to an arbitrary reference filter . this is a quite general method . the photometric redshift is the redshift @xmath15 corresponding to the minimum @xmath9 value . as expected , the accuracy on the photometric redshift determination depends strongly on the set of filters used ( number and wavelength coverage ) ; this is illustrated in the next paragraph .    in order to apply this method to quasars , the template set is made by series of simulated quasar spectra . these spectra have been constructed by varying the slope of the power - law spectra ( spectral index ) in the optical ( between 0.0 and 1.0 , in 3 steps ) while keeping the ultraviolet index , @xmath16 , constant at 1.76 ( a value compatible with wang et al . , 1998 ) . emission lines are included ( ly@xmath17 , ly@xmath18 , ciii , civ , mgii , siiv , h@xmath17 , h@xmath18 and h@xmath19 ) with gaussian profiles and typical intensities ( e.g. peterson , 1997 ) , as well as the small blue bump , centered at 3000  . the redshift value is set between 0 and 7 . totally , there are 213 quasar spectra , with 71 redshift steps ( d@xmath20 ) , with 3 simulated spectra on each step , obtained by varying the optical spectral index . the ly@xmath17 forest has been modelized according to madau ( 1995 ) , while no reddening has been included . we also included the stellar library of pickles ( 1998 ) , which contains 108 spectra of stellar types ranging from o to m , including dwarfs and giant stars . additionally , we included 22 white dwarf and carbon stars spectra , which are known to contaminate the quasar candidates catalogues .    when applying the method to a given sample , the discrimination between stars and quasars is performed in two steps . a first distinction between stars and quasars is made by identifying as star ( quasars ) the objects that show a better fit to one of the stellar ( quasar ) spectra . a further selection on the remaining sample can be made , based on the value of the reduced @xmath9 , with stars identified as objects excluded as quasars at the 95% or 99% confidence levels . in the present paper , the terms `` @xmath9 selection '' or `` no @xmath9 selection '' will apply for such a distinction . when an object is identified as quasar candidate , a photometric redshift estimation is automatically provided . the number of degrees of freedom , @xmath21 , is equal to the number of filters minus 2 : we lose one degree by normalizing the flux to a reference filter and one degree because of the `` best fit '' procedure . in this way , the separation between stars and quasars is reduced to a one  dimensional fit between the data ( the observed sed ) and the spectral templates .    as expected , and as confirmed by the simulations , the quality of the quasar / star discrimination depends on the filter combination , as it happens for the photometric redshift estimate in the case of galaxies . not only the number of filters is important , but also the wavelength range covered . another important point is the width of the filters which must be suited to the width of quasar spectral features . a realistic determination of the photometric errors is the key to a successful issue for this method . the first test of this method has been done on mixed ( quasars + stars ) simulated catalogues . in this paragraph we intend to show the importance of the set of filters used for optimizing the results . first , let us define an arbitrary and useful quantity which is aimed to quantify the efficiency of this procedure . let @xmath22 be the number of quasar candidates resulting from the identification technique , @xmath23 the number of actual quasars found as candidates ( real quasars effectively selected ) , and @xmath24 the number of expected quasars . in principle , @xmath24 is a quantity which can be estimated from previous surveys . we define the `` _ efficiency _ '' of the method as follows : + _ efficiency_=_completeness _ @xmath0 _ confirmation rate _ , + where _ completeness_@xmath25/@xmath26 and _ confirmation rate _ @xmath25/@xmath27 . obviously , both the _ completeness _ and the _ confirmation rate _ take values in the interval [ 0,1 ] and so does the _ efficiency_.    the number of stars in our simulated catalogues was given by the galaxy model of robin et al . , 1995 ( r95 ) . the star counts correspond to the six deep multicolor survey ( hereafter dms ) fields , since we primarily tested our results with the results of this survey . the fields with their surface ( in square arc minutes ) and the estimated star numbers are given in table [ gal_lat ] . quasar numbers were based on previous surveys ( hartwick & schade , 1990 , hereafter hs90 ) . catalogues contain @xmath28100 quasars and @xmath283000 stars , with a magnitude limit of @xmath29=22.2 and photometric errors scaling with magnitude , with typical values of the order of 0.1 magnitudes . quasars are uniformly distributed in redshift , in order to obtain reasonable statistical results at all redshifts . .the six dms fields , galactic coordinates , surfaces and estimated star numbers . [ cols="<,^,^,^",options="header " , ]     a closer look at fig . [ chi2sel ] and table [ filter_comb ] reveals that the filter combination is , most of the time , more important than the number of filters . in general , filters like u or the infrared filters h and/or k are very useful : they add much photometric information constraining the number of candidates . neighboring filters like b , v , r , when used all together , decrease the _ efficiency _ because , as mentioned above , they increase the number of the degrees of freedom without adding any essential information . a comparison of the two figures reveals several trends . the _ completeness _ is almost not influenced by a @xmath9 selection in the case of 3-filter combinations ( such as ubv , bri , bvr , bvi , or ubr ) . it is , however , strongly influenced by a @xmath9 selection when the filter u is present , in sets of 4 or more filters . note that this quantity decreases significantly , when a @xmath9 selection is imposed , for filter combinations with index higher than 26 , which are all sets containing the u - filter . furthermore , the _ confirmation rate _ is generally increasing with increasing number of filters . the decision whether to make a @xmath9 selection or not depends on the available time for spectroscopy and the scientific objectives of the project .    for this purpose , we made an approximative estimation of the time needed for making the spectroscopy in order to cover an area of 1 square degree , @xmath30 , for the same 40 filter combinations , given a detector area , @xmath31 , a magnitude limit , @xmath29 , a quasar candidate density , @xmath32 ( depending on this limit as well as on the filter combination ) , a number of slits , @xmath33 and an exposure time , @xmath34 , with and without any @xmath9 selection . we calculate @xmath30 as : @xmath35    fig . [ teltime ] illustrates @xmath30 ( in hours ) for @xmath36 , s=0.06 square degrees , g=10 and @xmath373 hours , without any preselection ( solid line ) and with a @xmath9 selection ( dashed line ) . note that the y - axis is logarithmic . as expected , there is an anticorrelation between @xmath30 and the _ confirmation rate _ , as it appears from the comparison of fig . [ chi2sel ] and [ teltime ] . the highest the _ confirmation rate _ , the less time is required . the problem of the quasar photometric redshift ( hereafter @xmath38 ) estimation is at least as complex as the one for galaxies . quasar spectral features ( mainly emission lines ) are rather narrow and the overall shape of the ultra - violet and optical continuum is a single power low , while the filters commonly used for photometry are wide - band filters , not very suitable for the quasar characteristics . all the above are in the origin of the dispersed values of the estimated quasar @xmath38 , when compared to the spectroscopic ones , especially for quasars with redshifts lower than @xmath39 . the filter combination is very important for the accurate determination of the @xmath38 . figure [ zphotsimps ] shows the photometric redshift versus the model redshift , for a catalogue of 1540 simulated quasars , uniformly distributed over the redshift range [ 0,6 ] , and for 4 different filter combinations . from top to bottom and left to right the following filter combinations are illustrated : bvri ( a ) , bvik ( b ) , ubvri ( c ) , ubvrik ( d ) . 0.2 cm    the use of the filter set bvri is clearly a bad choice : these four filters cover a relatively short wavelength range and their effective wavelengths lie close to each other . the @xmath38 in this case is poorly determined for @xmath40 , while there is also a large dispersion at @xmath41 . if we replace the r - filter by the k - filter ( figure [ zphotsimps]b ) , the dispersion becomes more important at redshifts higher than 4 , since this is the redshift at which ly@xmath42 should enter the r - band , but it decreases the dispersion at redshifts lower than 3 , due to the very large wavelength range covered . adding the filter u to the bvri combination decreases the dispersion at both high and low redshifts ( figure [ zphotsimps]c ) , but one can still feel the absence of an infrared filter , whose importance is seen in figure [ zphotsimps]d . the periodic accumulations of points at different redshifts are due to the passage of ly@xmath42 from one filter to the other . as a general rule , the determination of the photometric redshifts should be rather accurate in the redshift range [ 3,5 ] even without the u or an infrared filter , but their presence is crucial for all other redshift ranges . the quality of the @xmath38 determination is likely to improve if real quasar spectra are used . however , this implies that a spectral library is required , containing some hundreds of quasar spectra , with redshifts spanning from 0 to @xmath43 and with a spectral coverage spanning form the ultra - violet to the infrared .      the _ completeness _ , the _ confirmation rate _ and , therefore , the _ efficiency _ of the method depend on the redshift . the demonstration of this dependence is easy to show for the _ completeness _ , but not for the _ confirmation rate_. the _ completeness _ within redshift bins has the same definition as for the entire sample : the ratio between the number of real quasars with spectroscopic redshifts ( hereafter @xmath44 ) within the examined bin , identified as candidates even if their @xmath38 is significantly different from their @xmath44 , over the estimated total number of quasars . figure [ cvsz ] illustrates the _ completeness _ versus the redshift , without @xmath9 selection , for 4 different filter sets : the solid , dotted , dashed and dashed - dotted lines correspond to bvik , ubvrih , uvk and bvri , respectively . the _ completeness _ over the whole redshift range for this 4 filter combinations is in all cases more than 85% , as shown in figure [ chi2sel]a . figure [ cvsz ] reveals that , independently of the filter combination , the _ completeness _ is always lower in the redshift range [ 2,3 ] due , as already mentioned , to the stellar  like colors of quasars . this fact could certainly increases the number of quasar candidates within this range but it could also be a possible source of incompleteness , since some quasars will inevitably be classified as stars . furthermore , different filter combinations offer a different _ completeness _ at different redshifts . the filter set bvik assures a _ completeness _ of over 70% in the critical redshift range [ 2,3 ] , which however remains low up to @xmath45 . uvk offers the possibility of a _ completeness _ of @xmath2880% for @xmath46 , while ubvrih results give a _ completeness _ over 90% for high-@xmath15 quasar samples ( @xmath47 ) . things become more complicated if one tries to make an equivalent evaluation of the _ confirmation rate _ within redshift bins . this time , the _ confirmation rate _ is the ratio between the number of real quasars with @xmath38 in the examined redshift bin , over the number of quasar candidates ( quasars + possibly stars ) with @xmath38 that belong to this particular redshift interval . obviously , this definition depends on the accuracy of the determination of the @xmath38 as well as on the width of the redshift bins . therefore , the dependence of the confirmation rate versus the redshift can be evaluated qualitatively but not quantitatively , and so it goes for the _ efficiency _ as well . if the final aim of a survey is assembling a quasar sample with redshifts within a restrained interval , all candidates with a @xmath38 that lies outside this interval will , in principle , be rejected . in the more realistic case of an all redshift survey , the most interesting strategy is to simply optimize the spectroscopy towards the optical or the infrared wavelengths . the spectroscopy of quasars with redshifts up to @xmath283 is much more interesting in the optical band , while higher redshift quasars should be spectroscopically observed in the infrared , since most of their emission lines are shifted to these wavelengths . the simulations have shown that , although the determination of quasar @xmath38 is sometimes dispersed at @xmath48 and @xmath49 ( figure [ zphotsimps ] ) , this effect hardly translates into high-@xmath15 quasars ( @xmath50 ) being wrongly identified as low-@xmath15 ones ( @xmath51 ) or vice versa . the method presented above has been tested on a series of samples found in the literature , giving rather satisfactory results . depending on the sample and the available information , the _ confirmation rate _ and/or the _ completeness _ have been calculated . the deep multicolor survey is an imaging survey covering 6 fields on a total area of 0.83 square degrees on the northern sky at high galactic latitude , carried out at the 4 m mayall telescope of the kitt peak national observatory . 6 filters have been used , covering the wavelength range from 3000 to 10000  : u , b , v , r , i75 and i86 , and a photometric catalogue of @xmath28 21000 stellar  like objects has been assembled , with the aim of conducting a multicolor quasar search ( osmer et al . ( 1998 ) and the references therein ) . 137 stars , 49 compact narrow emission line galaxies ( hereafter cnelgs ) , and 54 quasars with redshifts spanning from 0.3 to 4.3 have already been spectroscopically identified . these objects were the 240 quasar candidates selected by hall et al . according to their position on color - color plots , a high fraction of them ( 194 objects ) selected by the uvx technique ( hall et al . , 1996 ) . this selection favored quasars with redshifts less than 2.2 . the redshift distribution of these objects can be seen in fig . [ histo ] , where the spectroscopic sample is plotted by a dotted line .    for our purposes , we cut the catalogue at b=22.3 , a limit at which the sample is almost 100% complete and only marginally contaminated ( @xmath28 8% ) by galaxies , and we only used objects for which photometry was available in all 6 filters . hence , 3 of the spectroscopically confirmed quasars were excluded , and we ended up with a catalogue of 3720 stellar - like objects . for the magnitude limit and the area covered by the survey , approximately 106@xmath52 quasars are expected to be found among the objects of the catalogue ( hs90 ) , resulting a _ completeness _ of at least 80% for the usual multicolor selection techniques ( hall , private communication ) and @xmath5388% for our method . lcccccccccccc & & & & & + survey&surface&mag.&filters&&&&&candidates & & + & covered&limit&&exp&cand&exp&95%&99%&no @xmath9&conf +   + dms&0.83sd&b=22.3&u , b , v , r,&@xmath54&3332&@xmath55&295&348&388&51 + & & & i75,i86&(r95)&&(hs90)&(40)&(45)&(45 ) & + hdf@xmath56&1.00sd&b=21.0&b , v , r&&&30&26&30&30&30 + parkes&&h=19.6&b , v , r , i,&&&129&99&107&125&129 + & & & j , h , k & & & & & & & +   +    table [ real_samples ] summarizes the characteristics of the different catalogues examined in this section : surface covered , magnitude limits , filters used , number of quasars , and number of candidates for different confidence levels . the expected number of stars on the dms field has been computed from r95 . in the case of the dms field , the number of spectroscopically confirmed quasars is also given between brackets .    among the 3720 objects in this catalogue , we identified 3332 as stars because their photometry shows a better match to the stellar templates . with our method , 6 of the 51 spectroscopically confirmed quasars ( dms0059 - 0056 , dms1358 - 0054 , dms1714 + 5012 , dms1714 + 5003 , dms1714 - 4959 and dms1358 - 0055 ) are misclassified as stars . only 7 of the 97 b @xmath57 candidates , spectroscopically identified as stars , but all 34 b @xmath57 cnelgs belong to our candidate list , if no @xmath9 selection is made . at a 99% confidence level , 25 cnelgs are classified as quasar candidates .    to summarize , we find 388 quasar candidates for 106 actual quasars , with an expected _ completeness _ of 45/51=0.88 and an expected _ confirmation rate _ of ( 106 @xmath0 0.88)/388=0.24 , if no @xmath9 selection is made . the equivalent values with 99% confidence level selection are 0.88 and 0.27 , respectively , for 348 quasar candidates ( see also table [ efficiency ] ) .    fig . [ histo ] displays with a solid line the estimated redshift distribution of our 388 quasar candidates . the accumulation of identifications around @xmath58 is due to the cnelgs , with @xmath44 in the range [ 0.1,0.8 ] typically , identified by our procedure as quasar candidates with @xmath59 . the broad distribution at @xmath60 $ ] is probably due to the stellar  like colors of quasars at these redshifts , which makes the distinction between stars and quasars particularly difficult by photometric means , increasing the number of quasar candidates , but it also corresponds to the expected maximum of the quasar space distribution . for comparison , we plot the redshift distribution of the 51 spectroscopically confirmed quasars ( dotted line ) . the above results are summarized in the first 6 columns of table [ efficiency ] , where the values of the _ completeness _ ( c ) , the _ confirmation rate _ ( cr ) and the _ efficiency _ ( e ) are shown in columns 3 and 4 . the values predicted by the simulations are also given for comparison in columns 5 and 6 , using the same filter combinations as for the real data , and keeping the filter responses as close as possible to those actually used . column 2 gives the same quantities for the dms and the usual multicolor selection techniques . note that the values given for the _ completeness _ and the _ efficiency _ are the lower estimations . the two methods give comparable values for all 3 quantities ( c , cr and e ) , with a slightly better efficiency for the multicolor technique when a @xmath9 criterion is imposed . l|ccccc|cccccccc & & & & & & & + & & 99%&no @xmath9&99%&no @xmath9&99%&no @xmath9&99%&no @xmath9&99%&no @xmath9&99%&no @xmath9 + c&0.80&0.88&0.88&0.79&0.86&1.00&1.00&0.74&0.76&0.83&0.97&0.52&0.81 + cr&0.23&0.27&0.24&0.20&0.21 + e&0.18&0.24&0.21&0.16&0.18 +    in order to estimate the contamination due to compact galaxies , we have added to the template library the spectra of blue starburst galaxies , obtained from the public model starburst99 code by leitherer et al . ( 1999 ) , as well as a delta - burst model taken from the new bruzual & charlot evolutionary code ( gissel98 , bruzual & charlot 1993 ) , and a set of 4 empirical seds compiled by coleman , wu and weedman ( 1980 ) ( hereafter cww ) to represent the local population of galaxies ( e , sbc , scd and i m ) . spectra from the starburst99 code were selected with solar metallicity and salpeter imf , taking different ages for an instantaneous burst ranging from 1 to 20 myr . in the case of the gissel98 delta - burst , the imf is that of miller & scalo ( 1979 ) , with solar metallicity and 51 different ages ranging from 0.001 to 20 gyr . cww spectra were extended to wavelengths @xmath61  and @xmath62  using the equivalent gissel98 spectra .    when applying the method to the dms sample of 3523 stellar objects ( we remind here that we only use objects with photometry in all 6 filters ) , using the template library extended with galaxy spectra , we find that @xmath63 of them are better fitted by stellar or quasar templates rather than by galaxy spectra . this fraction is only 56% for the subsample of 97 objects spectroscopically identified as stars . among the subsample of the 51 objects spectroscopically identified as quasars , there are 38% that show a better fit to the galaxy spectra , and 80% among them correspond in particular to genuine young starbursts according to the best - fit spectra , with ages ranging between 1 and 20 myrs . the result is different when the same extended template library is applied to the 34 b @xmath57 dms cnelgs . in this case , 85% of the sample ( 29 objects ) is well fitted by galaxy spectra rather than quasar or stellar templates . the fraction of best - fit galaxy templates corresponding to young starbursts is @xmath64 . generally speaking , a typical fraction between @xmath65 and @xmath66 of the quasar candidates is also compatible with young starbursts , and cnelgs are hardly selected as quasars , when supplementary galaxy templates are used . to summarize , adding star - bursting and early - type galaxy templates does not give a better solution to the problem of selecting quasar candidates : if star - burst candidates are included in the list , the number of candidates is multiplied by a factor of 4 to 5 . on the other side , objects like cnelgs will not be mixed with the quasar candidates but we take the risk to miss @xmath67 of the quasars , identified as star - bursting galaxies . however , for the purposes of a survey , all objects will be tested through the same , general pipeline . some objects ( like quasars and compact blue galaxies ) will inevitably end up as candidates in more than one list , depending on the different selection criteria . figure [ zmzpps ] left presents the @xmath38 versus @xmath44 for the 51 dms quasars . for 23 of them , @xmath68 . we note an important dispersion for redshifts lower than typically 2.5 : quasars with a @xmath69 $ ] are better fitted by templates of quasars with @xmath70 $ ] . this is also the tendency arising from simulations and shown on figure [ zphotsimps ] ( upper right ) , for the filter combination ubvri . figure [ zmzpps ] right presents the @xmath38 versus @xmath44 for a new sample of 46 spectroscopically identified quasars ( hall et al . , 2000 ) . 17 of them have @xmath68 . this time we note that quasars with a @xmath71 $ ] are better fitted by quasars with @xmath72 $ ] , and this trend is also seen on the simulations ( figure [ zphotsimps ] upper right ) for the same set of filters . the _ completeness _ has also been tested on two recently published quasar catalogues . liu et al . ( 1999 ) published a list of 30 spectroscopically confirmed quasars with their photometry in b , v and r. these objects have @xmath74 and @xmath75 , and were selected on a one square degree field centered on the hdf north . they were firstly identified as quasar candidates by their ( b - v ) ( v - r ) colors , along with 31 other candidates that , after the spectroscopy , turned out to be stars .    with our method we identified all 30 quasars as quasar candidates . however , the limited number of filters did not allow us to reduce the number of stars selected as quasar candidates , 30 out of 31 were in our own list . we also applied our method to the parkes subsample published by francis et al . only point - like objects , with the photometry in all 7 filters ( b , v , r , i , j , h , k@xmath76 ) were kept , thus 129 out of the 157 quasars of the initial sample remain in our catalogue . at a confidence level of 99% , 107 among these objects were classified as quasar candidates ( _ completeness_=0.83 ) . if no @xmath77 selection is applied , the _ completeness _ is equal to 0.97 . columns 7 to 14 in table [ efficiency ] summarize these results . simulation results are also given for comparison . note that in both cases , the _ completeness _ for the real samples is found to be even higher than expected according to the simulations . _ confirmation rate _ and _ efficiency _ can not be computed for these two samples , since the entire point - like object catalogues ( including stars ) are not available . an attempt to apply our method to the eis  wide ( patch b ) data was made . the results are approximatively the same as those obtained in multicolor procedures . this was somewhat expected : the multicolor approach proposed here , when applied to a 3 filter data set , is equivalent to a color  selection technique . note that 3 filters have been used to carry out this survey : b , v and i. on a ( b - v ) ( v - i ) diagram , quasar evolutionary tracks up to a redshift of @xmath78 and the main sequence stars mix together . the quasar multicolor selection technique presented in this paper has some advantages as compared to the `` traditional '' color - color methods . first , it has the advantage of selecting quasar candidates even at redshifts where quasar and stellar colors are very much alike ( 2.5  3.0 ) . however , the _ efficiency _ varies with the redshift . independently of the filter set , some quasars with redshifts lying in the interval [ 2.5,3 ] will be skipped . furthermore , the number of quasar candidates with estimated @xmath38 within this range will be far higher in comparison to other redshift bins . the set of filters used plays an important role in these results . adding a filter is not necessarily equivalent to adding information and does not always improve the results , but imposes more stringent constraints by increasing the number of degrees of freedom . for this reason , a @xmath9 selection does not improve the results in all cases . this technique gives also an estimate of the ( photometric ) redshift of the candidates , and it is thus able to improve the spectroscopic follow up of a given sample . in some cases , which depend ( again ) on the filter combination and the redshift of the objects , this determination can be accurate ( @xmath79 ) . in most cases , quasars at redshift @xmath40 will not be misidentified at @xmath80 and vice versa , therefore the spectroscopy can be reliably targeted towards the optical or the infrared bands .    for the purposes of a survey aiming in assembling a quasar catalogue with a minimum cost , adding star - bursting and early - type galaxy templates does not appear as a priority . on the contrary , this addition will cause either the increase of the candidates by a factor of almost 5 ( if star - burst candidates are considered as possible quasars ) or the decrease of the _ completeness _ @xmath81 ( because an important fraction of real quasars will be better fitted by star - burst galaxies templates ) . for practical purposes ( e.g the virmos survey ) , when all the sources examined through the same pipeline , a fraction of objects identified as quasars will be also found within the `` starburst '' sample , but this does not affect the present selection .    the same method can be extended towards different directions . first , it can be adapted to particular objects like very red quasars , by adding the suited templates . however , this has to be done thoroughly , to prevent from a possible and undesirable increase of the quasar candidates . it can also be used on spatially resolved objects , like seyfert galaxies , because by construction , the software is independent of the morphology of the objects . this can be done by adding a stellar component to the agn simulated spectra already used . last but not least , it can be extended to very high redshifts ( higher than 6 ) , after an appropriate modeling of the intergalactic absorption . as already mentioned , this software has been developed for the needs of preparing the quasars / agn sample assembled in the future virmos survey . however , the results and conclusions of this paper apply to all multicolor surveys in general . a cross - identification between optical and x - ray sources can be extremely useful , when x - ray data are available , and can increase the _ efficiency _ , as quasar candidates with a compact x - ray counterpart will have an even higher chance of being real quasars . we would like to thank m. bolzonella , j. p. picat and m. elvis for useful discussions . we would also like to thank the referee , dr . p. hall for his very useful comments that led to the improvement of our work . bolzonella , m. , miralles , j. m. & pell , r. , 2000 , submitted to a&a , pre - print astro - ph/0003380 boyle b.j . ; jones , l.r . , shanks , t. , 1991 , mnras , 251 , 482 bruzual , g. , charlot , s. 1993 , apj 405 , 538 cuqs : http://sa1.star.uclan.ac.uk/  prn / cuqs.html cohen , j. , hogg , d. , blandford , r. , cowie , l. , hu , e. , songaila , a. , shopbell , p. & richberg , k. , 1999 , accepted for publication in apj , astro - ph/9912048 coleman , d.g . , wu , c.c . , weedman , d.w . 1980 , apjs 43 , 393 fan , x. , 1999 , aj , 117 , 2528 glazebrook , k. , ellis , r. , colless , m. , broadhurst , t. , allington - smith , j. & tanvir , n. , 1995 , mnras , 273 , 157 hall , p. b. , osmer , p. s. , green , r. f. , porter , a. c. , warren , s. j. 1996 , apjss , 104 , 185 hall , p. b. , osmer , p. s. , green , r. f. , porter , a. c. , warren , s. j. 1996 , apj , 462 , 614 hartwick f. & schade d. , 1990 , ara&a , 28 , 437 jarvis , r.m . & macalpine , g.m . , 1998 , aj , 116 , 2624 kennefick j.d . , osmer , p.s . , hall p.b . , green , r.f . , 1997 , aj , 114 , 2269 krisciunas k. , margon , b. , szkody , p. , 1998 , pasp , 110 , 1342 leitherer , c. , schaerer , d. , goldader , j. d. , delgado , r. m. , robert , c. , kune , d. f. , de mello , d.f . , devost , d. , heckman , t. m. , 1999 , apjs 123 , 3 madau p. , 1995 , apj , 441 , 18 miller , g.e . , scalo , j.m . 1979 , apjs 41 , 513 miller , l. & mitchell , p. s. , 1988 , asp conference series , volume 2 , proceedings of a workshop on optical surveys for quasars , ( san francisco : asp ) , ed . by patrick osmer and m.m . phillips , p. 114 miralles j.m . & pell , r. , 1998 , astro - ph/9801062 newberg , h.j . & yanny , b. , 1997 , apjss , 113 , 89 osmer , p. , kennefick j.d . , hall , p.b . , green , r.f . , 1998 , apjs , 119 , 189 , and http://www.astronomy.ohio-state.edu/  posmer / dms/ pell , r. , kneib , j .- bolzonella , m. , miralles , j .- , 1999 , proceedings of the `` photometric redshifts and high redshift galaxies '' , pasadena , pasp conf . 191 , p 247 . [ astro - ph : 9907054 ] peri , f. , iovino , a. & hickson , p. , 1997 , `` quasar detection using multi - narrow - band photometry '' , in _ science with liquid mirrors telescopes _ , proceedings of the marseille meeting on lmts , april 1997 , in press . peterson , b. m. , 1997 , in _ an introduction to active galactic nuclei _ , cambridge university press pickles a. j. , 1998 , pasp , 110 , 863 prandoni i. , wichmann r. , da costa l. , benoist , c. , mendez , r. , nonino , m. , olsen , l. f. , wicenec , a. , zaggia , s. , bertin , e. , deul , e. , erben , t. , guarnieri , m. d. , hook , i. , hook , r. , scodeggio , m. , slijkhuis , r. , 1999 , a&a , 345 , 448 robin , a. , haywood , m. , gazelle , f. , bienaym , o. , crz , m. , oblak , e. , guglielmo , f. 1995 , `` model of stellar population synthesis of the galaxy '' , http://www.obs-besancon.fr/www/modele/modele_ang.html http://www.sdss.org/news/releases/19981208.qso.html wang , t. g. , lu , y. j. & zhou , y. y. , 1998 , apj , 493 , 1 warren s.j . , hewett , p.c . , osmer , p.s . 1994 , apj , 421 , 412 wolf , c. , meisenheimer , k. , roeser , h .- j . , 1999 , a&a , 343 , 399
Abstract:

Answer: we present a quasar candidate identification technique based on multicolor photometry . the traditional multi - dimensional method ( 2 @xmath0 n dimensions , where n is the number of the color - color diagrams ) is reduced to a one - dimensional technique , which consists in a standard fitting procedure , where the observed spectral energy distributions are compared to quasar simulated spectra and stellar templates . this new multicolor approach is firstly applied to simulated catalogues and its _ efficiency _ is examined in various redshift ranges , as a function of the filter combination and the available observing time for spectroscopy . we conclude that this method is better suited than the usual multicolor selection techniques to quasar identification , especially for high - redshift quasars . the application of the method to real quasar samples found in the literature results in an _ efficiency _ comparable to the one obtained from the use of color - color diagrams . the major advantage of the new method is the estimation of the _ photometric _ redshift of quasar candidates , enabling , in almost all cases , spectroscopy to be targeted to best suited wavelength ranges .