arxiv:astro-ph/ v1 29 Oct 2003

Transcription

1 Astronomy & Astrophysics manuscript no. ms418 February, 8 (DOI: will be inserte by han later) Approimate angular iameter istance in a locally inhomogeneous universe with nonero cosmological constant M.Demianski 1,, R. e Ritis 3,4, A.A.Marino 4, E. Pieipalumbo 3,5 arxiv:astro-ph/3183v1 9 Oct 3 1 Institute for Theoretical Physics, University of Warsaw, Warsaw, Polan Department of Astronomy, Williams College, Williamstown, MA 167 USA 3 Dipartimento i Sciene Fisiche, Università i Napoli Feerico II, Compl. Univ. Monte S. Angelo, 816 Napoli, Italia 4 Istituto Naionale i Fisica Nucleare, Se. Napoli, Via Cinthia, Compl. Univ. Monte S. Angelo, 816 Napoli, Italia 5 Osservatorio Astronomico i Capoimonte, Via Moiariello, Napoli, Italia Receive / Accepte Abstract. We iscuss the general an approimate angular iameter istance in the Frieman-Robertson-Walker cosmological moels with nonero cosmological constant. We moify the equation for the angular iameter istance by taking into account the fact that locally the istribution of matter is non homogeneous. We present eact solutions of this equation in a few special cases. We propose an approimate analytic solution of this equation which is simple enough an sufficiently accurate to be useful in practical applications. Key wors. Cosmology, angular iameter istance, gravitational lensing 1. Introuction Recent observations of the type Ia supernovae an CMB anisotropy strongly inicate that the total matter-energy ensity of the universe is now ominate by some kin of vacuum energy also calle ark energy or the cosmological constant Λ (Perlmutter 1997, Perlmutter et al.1999, Riess et al. 1998, Riess ). The origin an nature of this vacuum energy remains unknown. There are several review articles proviing thorough iscussion of the history, interpretations an problems connecte with the vacuum energy an observational constrains (Zel ovich 1967, Weinberg 1989, Carroll et al.199). The type Ia supernovae have been alreay observe at reshifts > 1. It is well known from galay surveys that galaies an clusters of galaies up to a scale of 1Gpc are istribute non homogeneously forming filaments, walls an unerense vois. This inicates that on similar scales also the ark matter is istribute non homogeneously. In this paper we analye the influence of local non homogeneities on the angular iameter istance in a universe with non ero cosmological constant. The angular iameter istance in a locally non homogeneous universe was iscusse by Zelovich (Dashveski & Zelovich 1965) see Sen offprint requests to: E.Pieipalumbo, ester@na.infn.it also (Dashveski & Slysh 1966) an (Weinberg 1989, Kayser et al.1997). Later Dyer& Roeer (197) use the so calle empty beam approimation to erive an equation for the angular iameter istance; for more etaile references see also (Schneier, Ehlers & Falco 199), (Kantowski 1998), an (Tomita & Futamase 1999). We follow the Dyer & Roeer metho to erive the equation for angular iameter istance in a locally non homogeneous universe with a cosmological constant. In the general case the equation for the angular iameter istance oes not have analytical solutions, it can be solve only numerically (Kantowski, Kao & Thomas ). We have foun an analytic approimate solution of this equation, which is simple an accurate enough to be useful in practical applications. This allowe us to fin an approimate epenence of the angular iameter istance on the basic cosmological parameters. The paper is organie as follows: we begin with the general form of Sachs equations escribing light propagation in an arbitrary spacetime an using the empty beam approimation we erive the equation for angular iameter istance. Then, we iscuss the angular iameter istance as observe from a location not at the origin (=). Finally, after iscussing some properties of the analytical solutions of the equation for the angular iameter istance in a locally non homogeneous universe, we propose an analytic approimation solution of this

2 M.Demianski et al.: Approimate angular iameter istance.. equation vali in a wie reshift interval (, 1). We have been motivate by the great avances in observing further an further objects an the construction of very powerful telescopes that provie the possibility to observe gravitational lensing by clusters of galaies, an supernovae at large istances. It is therefore necessary to evelop a more accurate formalism to escribe the istance reshift relation an gravitational lensing by high reshift objects. To achieve this goal we use cosmological moels with realistic values of the basic cosmological parameters such as the Hubble constant an the average matter-energy ensity. In concluing remarks we summarie our results an iscuss some perspectives.. Propagation of light an local non homogeneities.1. General consierations Let us consier a beam of light emanating from a source S in an arbitrary spacetime escribe by the metric tensor g αβ. The light rays propagate along a null surface Σ, which is etermine by the eikonal equation g αβ Σ, α Σ, β =. (1) A light ray is ientifie with a null geoesic on Σ with the tangent vector k α = Σ, α. The light rays in the beam can be escribe by α = α (v, y a ), where v is an affine parameter, an y a (a=1,, 3) are three parameters specifying ifferent rays. The vector fiel tangent to the light ray congruence, k α = α v = Σ, α, etermines two optical scalars, the epansion θ an the shear σ, which are efine by θ 1 kα ;α, σ k α;β m α m β, () where m α = 1 (ξ α iη α ) is a comple vector spanning the spacelike -space (the screen space) orthogonal to k α (k α m α = ). Since k α = Σ, α the vorticity connecte with the light beam is ero in all our consierations; therefore in this case σ an θ fully characterie the congruence. These two optical scalars satisfy the Sachs (Sachs & Kristian 1966) propagation equations θ + θ + σ = 1 R αβk α k β, (3) σ + θσ = 1 C αβγδ m α k β m γ k δ, (4) where ot enotes the erivative with respect to v, R αβ is the Ricci tensor, an C αβγδ is the Weyl tensor. Eqs. (3) an (4) follow from the Ricci ientity. The optical scalars θ an σ escribe the relative rate of change of an infinitesimal area A of the cross section of the beam of light rays an its istortion. The epansion θ is relate to the relative change of an infinitesimal surface area A of the beam s cross section by θ = 1 lna v. (5) Let us use these equations to stuy the propagation of light in the Frieman-Robertson-Walker (FRW) spacetime. The FRW spacetime is conformally flat, so that C αβγδ =. From Eq. (4) it follows that, in the FRW spacetime if the shear of the null ray congruence is initially ero, then it always vanishes. Therefore, assuming that the light beam emanating from the source S has vanishing shear, we can isregar the shear parameter altogether (empty beam approimation). Using Eq. (5) we can rewrite equation (3) in the form 1 A + R αβk α k β A =. (6) An observer moving with the 4-velocity vector u α will associate with the light-ray a circular frequency ω = cu α k α. Different observers assign ifferent frequencies to the same light ray. The shift of frequency, as measure by an observer comoving with the source an an arbitrary observer, is relate to their relative velocity or the reshift by 1 + = ω = c u α k α, (7) ω o ω o where ω is the frequency measure by a comoving observer an ω by an arbitrary observer. Differentiating this equation with respect to the affine parameter v, we obtain v = c ω o k α k β u α;β. (8) Since the angular iameter istance D is proportional to A we can rewrite Eq.(6) using D instea of A an, at the same time, replace the affine parameter v by the reshift ; we obtain ( ) ( D ) v + D v + 4πG c 4 T αβk α k β D =, (9) where we have use the Einstein equations to replace the Ricci tensor by the energy-momentum tensor. A solution of Eq. (9) is relate to the angular iameter istance if it satisfies the following initial conitions: D() = =, D() = = c. H.. The angular iameter istance between two objects at ifferent reshifts (1) In orer to use solutions of the Eq. (9) to escribe gravitational lenses we have to introuce the angular iameter istance between the source an the lens D( l, s ), where l an s enote corresponingly the reshifts of the lens an the source. Let D( 1, ) enote the angular iameter istance between a fictitious observer at 1 an a source at ; of course, D(, ) = D(). Suppose that we know the general solution of Eq. (9) for D(), which satisfies the initial conitions (1); let us efine D( 1, ) by D( 1, ) = c (1 + 1 )D( 1 )D() H 1 D ( )g( ), (11)

3 M.Demianski et al.: Approimate angular iameter istance.. 3 subject to the conitions D( 1, ) =1 =, D( 1, ) =1 = sign( 1 ) c H g( 1 ). (1) We see that D( 1, ) satisfies Eq. (9), if the function g() is a solution of the equation lng() = so that g() = g ep v ( v ), (13) v ( v ), (14) where g is an arbitrary constant of integration. It is easy to relate the function g() to the Hubble function H(). In fact, using t an = (1 + )H(), (15) H() = (1 + ), (16) v H it follows that Eq. (13) can be rewritten as [ ln (1 + ) H() H lng() = so, ], (17) g() = (1 + ) H() H, (18) where we have impose the initial conition g() = 1. Let us remin that one of the Frieman cosmological equations can be rewritten as H() = H Ωm (1 + ) 3 + Ω k (1 + ) + Ω Λ, (19) where H is the present value of the Hubble constant, Ω m = 8πG m 3H is the matter ensity parameter, Ω k = kc R H is the curvature parameter an Ω Λ = Λc is the 3H cosmological constant parameter. The omega parameters satisfy a simple algebraic relation Ω m + Ω k + Ω Λ = 1. () The initial conitions (13) epresse in terms of H() become D( 1, ) =1 =, (1) D( c 1, ) =1 = sign( 1 ) (1 + 1 )H( 1 ). ().3. The angular iameter istance To apply the angular iameter istance to the realistic universe it is necessary to take into account local inhomogeneities in the istribution of matter. Unfortunately, so far an acceptable averaging proceure for smoothing out local inhomogeneities has not been evelope (Krasinski 1997). Following Dyer & Roeer (Dyer & Roeer 197), we introuce a phenomenological parameter α = 1 clumps < > calle the clumpiness parameter, which is relate to the amount of matter in clumps relative to the amount of matter istribute uniformly (Dashveski & Zelovich 1965, Tomita et al.1999). Therefore, Eq. (9), in the FRW case, can be rewritten in the form ( ) ( D ) v + D v + 3 αω m(1 + ) 5 D =. (3) Let us note however that this equation oes not fully escribe the influence of non homogeneities in the istribution of matter on the propagation of light. It takes into account only the fact that in locally non homogeneous universe a light beam encounters less matter than in the FRW moel with the same average matter ensity. Equation (3) oes not inclue the effects of shear, which in general is non ero in a locally non homogeneous universe (see, for eample SEF Chapter 4 an 11). The effects of shear have been estimate with the help of numerical simulations (see, for eample Schneier & Weiss 1988, Liner 1998, Jarosynski 1991, Watanabe & Sasaki 199) but satisfactory analytic treatment is still lacking. It is customary to measure cosmological istances in units of c/h, therefore we introuce the imensionless angular iameter istance r = DH /c. Using Eq. (16) an Eq. (19), we finally obtain (1 + ) [ Ω m (1 + ) 3 + Ω k (1 + ) ] r + Ω Λ + (4) ( ) 7 r Ω m(1 + ) 3 + 3Ω k (1 + ) + Ω Λ + 3 αω m(1 + ) r =, with initial conitions r() = =, (5) r() = = 1. (6) Eq. (4) can be cast into a ifferent form by using u = 1/(1 + ) as an inepenent variable instea of ; we get u (Ω m + Ω k u + Ω Λ u 3 ) r u u(3 Ω m + Ω k u) r u + 3 αω mr =. (7) On the other han, using the cosmic time t as an inepenent variable, Eq. (3) assumes the form r t H(t)r t + 4πGα m(t)r =. (8)

4 4 M.Demianski et al.: Approimate angular iameter istance.. This equation was for the first time introuce by Dashevski & Zelovich (Dashveski & Zelovich 1965) (see also Dashveski & Slysh 1966), an more recently Kayser et al. (Kayser et al.1997) an (Kantowski 1998) have use it to erive an equation similar to Eq. (5). In Eq. (8) the clumpiness parameter α is usually consiere as a constant. However, in (Dashveski & Zelovich 1965 an Kayser et al.1997), α is allowe to vary with time, but only the case α = constant is really consiere. For a iscussion of the general case when α epens on see, for eample, the paper by Liner (Liner 1988). To give an eample of our proceure of constructing the angular iameter istance between a fictitious observer at 1 an a source at, let us consier the simple case of α = constant, Ω Λ =, an Ω k =. In this case it is easy to integrate the equation (), we obtain r() = (1 + )β (1 + ) β β(1 + ) 5/4, (9) where β = α, is the general solution of Eq. (5) with the impose initial conitions. This solution can be foun, for eample, in SEF (Schneier, Ehlers & Falco 199), (see Eq. (4.56) there). The function g() can be easily compute in the flat Λ = Frieman-Robertson-Walker cosmological moel; in fact, in this case we have H() = H (1 + ) + 1, (3) so that, accoring to the Eq. (18), we have g() = (1 + ) (31) Thus using (11), we get the familiar solution [ ] D( l, s ) = 1 (1 + ) β 5 4 (1 + 1) β 1 4, (3) β (1 + 1 ) β+ 1 4 (1 + 1 ) β+5 4 foun in SEF. 3. Eact solutions In the general form the Eq. (4) is very complicate. From the mathematical point of view this equation is of Fuchsian type with four regular singular points an a regular singular point at infinity (Ince 1964). General properties of this equation have been etensively stuie by Kantowski (Kantowski 1998, Kantowski, Kao & Thomas, Kantowski & Thomas 1). Following Kantowski let us cast the Eq. (4) into a ifferent form replacing 1 + = an introucing r() = h() 1+, we obtain ( Ωm 3 + Ω k ) h + Ω Λ + ( 3 (α 1)Ω m Ω k ( ) 3 Ω m + Ω k h + ) h =. (33) When Ω k by rescaling, = Ω k Ω m this equation can be turne into the Heun equation ( Ω ΛΩ m Ω 3 ) h k + (3 + 1) h + [ 3 (α 1) 1]h =. (34) The angular iameter istance can be epresse in terms of basic solutions of the Heun equation as r() = h(ωm Ω k (1 + )), (35) 1 + provie that h( Ωm Ω k ) = an h ( Ωm Ω k ) = Ω k Ω m. The generic solution for the angular iameter istance epresse in terms of the Heun functions is very complicate, in an eplicit form it is given by Kantowski (Kantowski 1998). When Ω k =, iviing Eq. (33) by Ω m, we obtain ( 3 + Ω Λ ) h Ω m + 3 h + 3 (α 1)h =. (36) By changing variables to η() = 1 + Ω Λ an introucing new inepenent variable P() = 3/4 h() this Ω m 3 equation can be transforme into the associate Legenre equation (1 η ) P η ηp η 1 5 4α ( (1 η )P =. (37) ) In this case (Ω k = ) the angular iameter istance is given by r() = Γ(1 δ 6 )Γ(1 + δ 6 ) (1 + ) 1/4 1 + Ω Λ (38) δ Ω m [ ( ) ( ) P δ P δ Ω Λ 1 + Ω m (1 + ) 3 ( ) ( 1 + Ω Λ P δ 6 Ω 1 m 6 P δ Ω Λ Ω m Ω Λ Ω m (1 + ) 3 )] where δ = 5 4α, an P µ ν () enotes the associate Legenre function of the first kin. In section 5 we will propose a simple approimate analytic solution of the general equation for the angular iameter istance which is reproucing reasonably well the eact numerical solution in the range of reshifts (, 1) far eceeing the range of reshifts of observe supernovae of type Ia. 4. Angular iameter istances in the gravitational lensing theory. The angular iameter istances, or their combinations, appear in the main equations an in the most important observational quantities of the gravitational lensing (GL) theory. This fact, together with the relation between the angular iameter istance (D A ) an the luminosity istance (D L ) D L = (1 + ) D A, (39),

5 M.Demianski et al.: Approimate angular iameter istance.. 5 makes the stuy of the equation for the angular iameter istance still more important, also from the point of view of better interpretation of observational ata. In this section we escribe very briefly the role of angular iameter istance in the gravitational lensing theory, mentioning only some effects in which it plays an important role. Let us begin with the general epression for the time elay between ifferent light rays reaching the observer { } D D s c t = (1 + ) (θ β) ψ(ξ) + constant, (4) D s where D, D s an D s are corresponingly the angular iameter istances to the eflector, to the source an between the eflector an the source. The term (θ β) in Eq. (4) represents the geometrical time elay an the other term is connecte with the non homogeneous istribution of matter (Schneier, Ehlers & Falco 199). The first term in Eq. (4) is proportional to an important combination of angular iameter istances namely, to D D s D s. This combination of angular iameter istances was introuce for the first time by Refsal (Refsal 1966), an was use to obtain an estimate of H from the time elay measurements of multiply image quasars. It is possible to rewrite it in the following way (1 + ) D D s D s = c H [χ( ) χ( s )] 1, (41) where the function χ is given by χ(, Ω m, Ω Λ, Ω k, α) r g() = H r ()(1 + ) H(). (4) So, χ is connecte with the general solution of Eq. (5), an it irectly appears in the epression for time elay. Let us now consier the cosmological lens equation. The metric escribing FRW cosmological moels is conformally flat. In this case the simplest way to erive the lens equation is to use the Fermat principle, so we have t =, (43) θ or β = θ R g ch (1 + )[χ( ) χ( s )] ψ θ, (44) where R g is the Schwarschil raius of the eflector. By enotingξ = D θ, η = D s β, an α(ξ) = R g ψ, we can D θ transform Eq.(44) into η = D s D ξ D s α(ξ), (45) which is formally ientical with the lens equation for << 1. As is apparent, in the lens equation (45) another imensionless combination of angular iameter istances, D s /D appears, besies the angular iameter istance itself, D s. In other wors, in the equations escribing gravitational lensing we fin quantities, which can be written in terms of the angular iameter istance D() an the χ() function. In the general case D() an χ() can be evaluate only numerically. In the net section we propose a simple analytical approimation for both functions. We hope that these approimate forms can be useful, for eample, in big numerical coes use to erive basic cosmological parameters from observational ata an to stuy for instance, statistical lensing, weak lensing, microlensing of QSO, etc. 5. The approimate epression for the angular iameter istance In the generic case the equation (4) oes not have analytical solutions (Kantowski 1998, Kantowski, Kao & Thomas ). From the mathematical point of view this equation is of the Fuchsian type (Ince 1964) with four regular singular points an a regular singular point at infinity. The solutions near each of the singular points, incluing the point at infinity, are given by the Riemann P-symbol an in general the solutions can be epane in a series of hypergeometric functions (Tricomi 1961) or epresse in terms of the Heun functions (Kantowski 1998). In practice, in the general case, the equation (4) with appropriate initial conitions is solve numerically. In Fig.1 we show a numerical solution of the D-R equation for α =.8, Ω Λ =.7, Ω k = an Ω m =.3. Using r Fig.1. The imensionless angular iameter istance as a function of the reshift in a flat cosmological moel with Ω m =.3, Ω Λ =.7 an α =.8. Mathematica we iscovere that there is a simple function r() which quite accurately reprouces the eact numerical solutions of the equation (4) for up to 1, it has the form r() =, (46) 1 + ( )

6 6 M.Demianski et al.: Approimate angular iameter istance.. where 1, an 3 are constants which epen on the parameters that specify the consiere cosmological moel. First of all please note that the function (46) automatically satisfies the impose initial conitions, so r() = an r () = 1. To epress the constants 1, an 3 in terms of Ω Λ, Ω m, Ω k an α we inserte the propose form of r() into the equation (4) an require that it be satisfie to the highest orer. In this way we obtaine that: 1 = ( Ω k + 3 Ω m ) 4 14 Ω k + (8 + 5 α) Ω m 6 Ω k + 7 Ω m + 4 Ω Λ, + 7 Ω k + (13 + α) Ω m (47) = 1 + Ω k + 3 Ω m, (48) 4 3 = 1 ( ( Ω k + 3Ω m ) 4 (5 Ω k + 9 Ω m ) + (49) 16 ) 8 (14 Ω k + (8 + 5α) Ω m ). 6 Ω k + 7 Ω m + 4 Ω Λ In Fig. we show the approimate solution of the equation (4) for α =.8, Ω Λ =.7, Ω k =, an Ω m =.3 an for comparison we also plot the eact solution. the eact numerical solution for χ() an for comparison we plot the approimate solution with the best fit value of γ = 6. χ Fig. 3. The eact imensionless function χ (soli line) for a moel with α =.8, Ω Λ =.7, Ω k =, an Ω m =.3. With ashe line we plot the analytic approimation with the best fit value of the parameter γ = 6..4 r Fig.. The eact imensionless angular iameter istance r is plotte (soli line) an compare with the analytic approimation (ashe line) in a moel with α =.8, Ω Λ =.7, Ω k = an Ω m =.3. In the consiere range of reshifts the relative error is less than 1%. Using Mathematica we have also foun a useful approimate analytical form for the χ function 1 χ() = ( 1 + ), (5) γ where γ is a constant. Unfortunately we have not been able to analytically relate γ to α an other cosmological parameters an the appropriate value of γ shoul be obtaine by the stanar fitting proceure. In Fig. 3 we show 6. Conclusions In this paper we iscuss the angular iameter istance in the Frieman-Robertson-Walker cosmological moels an consier the case when the cosmological constant an the curvature of space coul be ifferent from ero. The effects of local non homogeneous istribution of matter are escribe by a phenomenological parameter α, consistently with the so calle empty beam approimation. Unfortunately, at the moment there are no generally accepte moels that escribe the istribution of baryonic an ark matter an therefore the influence of inhomogeneities of matter istribution can be inclue only at this approimate level. In the generic case the equation (4) is of a Fuchsian type, with four regular singular points an one regular singular point at infinity. The general solution of this type of orinary ifferential equation is given in terms of the Heun functions (Kantowski 1998, Demianski et al. ). However the eact solution is so complicate that it is useless in practical applications (Kantowski 1998, Kantowski, Kao & Thomas ). Therefore we have propose an approimate analytic solution, simple enough to be use in many applications an at the same time it is sufficiently accurate, at least in the interesting range of reshifts ( 1). In Fig. we compare the eact numerical solution of the equation for the angular iameter istance with the approimate one. The approimate

7 M.Demianski et al.: Approimate angular iameter istance.. 7 solution reprouces the eact curve quite well an the relative error oes not ecee 1%. Following SEF, we have foun the function χ which appears in the epression for time elay as well as in the lens equation, an which naturally appears in the epression for angular iameter istance between two arbitrary objects at reshifts 1 an ( see Eq. (11)). We have also propose an approimate analytical form of the function χ which epens only on one parameter γ but unfortunately γ has to be fie by a stanar fitting proceure (see Fig. 3). Our approimations have been alreay applie in complicate coes use to stuy the statistical lensing (Perrotta et al.1). Finally we woul like to stress that from our analysis it follows that variations in the angular iameter istance cause by the presence of cosmological constant are quite similar to the effects of a non homogeneous istribution of matter escribe here by the clumpiness parameter α. In Fig. 4 we plot the angular iameter istance for two moels, one with a homogeneous istribution of matter (α = 1) an Λ an another with a non homogeneous istribution of matter (α =.7) an Λ =. We see that an inhomogeneous istribution of matter can mimic the effect of a non ero cosmological constant. This is an important observation in view of the recent conclusions base on observations of high reshift type Ia supernovae that the cosmological constant is ifferent form ero (Perlmutter 1997, Riess )..4.3 r Fig. 4. The angular iameter istance for two moels, one with a homogeneous istribution of matter (α = 1) an Λ (soli line) an another with a non homogeneous istribution of matter (α =.7) an Λ = (ashe line). 7. Acknowlegments It is a pleasure to thank V.F. Carone, G. Covone, C. Rubano, P. Scuellaro, an M. Sereno for the iscussions we ha on the manuscript. E.P. is partially supporte by the P.R.I.N. SINTESI. This work was supporte in part by EC network HPRN-CT--14. M.D. wants to acknowlege partial support of the Polish State Committee for Scientific Research through grant Nr. -P3D an to thank the Theoretical Astrophysics Center in Copenhagen for warm hospitality. 8. Appeni As was alreay note the equation for the angular iameter istance in a locally non homogeneous Frieman- Robertson-Walker universe with non ero cosmological constant can be reuce to the Heun equation (see Eq. 8) (Heun, 1889). The Heun ifferential equation generalies the Gauss hypergeometric equation, it has one more finite regular singular point (see Ince, 1964). This equation often appears in physical problems in particular in stuying of iffusion, wave propagation, heat an mass transfer, an magnetohyroynamics (see Ronveau, 1995). In the general form the Heun equation can be written as H y + ( a1 y y 1 + a y y + a 3 y y 3 ) H y + (a 4 y q) H =, (51) (y y 1 )(y y )(y y 3 ) where a 1, a, a 3, a 4, q, y 1, y an y 3 are constants, or in a canonical form as (Bateman & Erélyi 1955) ( 1)( a) H + { (α + β + 1) [α + β a(γ + δ) δ] + aγ } H + (αβ q)h =, (5) where α, β, γ, δ, a, an q are constants. When γ is not an integer, the general solution of the Heun equation can be written as: H() = C 1 F(a, q, α, β, γ, δ, ) + (53) C 1 γ F(a, q 1, α γ + 1, β γ + 1, γ, δ, ), where C 1 an C are constants an q 1 = q + (α γ + 1)(β γ + 1) αβ + δ(γ 1). Following our success in fining an approimate solution of the equation for angular iameter istance we woul like to present an approimate solution of the Heun equation. It has the form: H() = B ( ) 1 (54) 1 ( ( ) ( ) ), ( ) where the parameter represents the point where the initial conitions are specifie an 1,, 3, B are constants. For the Heun equation we aopt the following initial conitions 1 : H[ ] =, (55) 1 To get the general initial conitions i.e. H( ) = A it is enough to a to the (55) a constant H + A.

8 8 M.Demianski et al.: Approimate angular iameter istance.. H = B. (56) Please note that the form (55) automatically satisfies the impose initial conitions. To epress the constants 1, an 3 in terms of α, β, γ, δ, a, an q, it is necessary to insert the propose form of H() into the canonical Heun equation (5) an require that it be satisfie to the highest orer. The constants 1, an 3 obtaine in this way are quite complicate. It will take about 3 pages to give them here. They are shown in an eplicit form at ester h Fig. 5. Comparison of the eact (soli line) with the approimate solution (ashe line) of the Heun equation. ε Fig. 6. Relative error of the approimate solution of the Heun equation as a function of. In Fig. 5 we compare an eact numerical solution of the Heun equation with the approimate solution for the same values of α, β, γ, δ, a an q. References Bateman, H., an Erélyi, A., Higher Transceental Functions, Vol. 3, 1955, McGraw-Hill Book Comp., New York Carroll, S.M., Press, W.H., Turner, E.L., 199, Ann. Rev. Astro. Astrophys. 3, 499 Dashveski V.M., Zelovich Ya.B., 1965, Soviet Astr., 8, 854 Dashveski V.M., Slysh V.I.,1966, Soviet Astr., 9, 671 Dyer C.C., Roeer R.C., 197, ApJ, 174, L115 Demianski M., e Ritis R.,Marino A.A., Pieipalumbo E., astro-ph/4376 Dyer C.C.,Oattes L.M., 1988, ApJ, 36,5 Dashevski V.M., Zelovich Y.B., 1965,Soviet. Astr., 8, 854. Etherington J.M.H., 1933, Phil.Mag., 15, 761 Heun, K., 1889, Math.Ann., 33, 161 Ince E.L, Orinary ifferential equation, 1964, Dover, New York Kayser R., Helbig P., Schramm T., 1997, A&A, 318, 68 Kayser R., Refsal S., 1988, A&A, 197, 63K Kantowski R., 1998, ApJ, 57,483 Kantowski R., Kao J.K., Thomas, R.C.,, ApJ, 545, 549 Kantowski R., an Thomas, R.C., 1, ApJ, 561, 491 Krasinski A., Inhomogeneous Cosmological Moels, 1997, Cambrige University Press, Cambrige Jarosynski M., 1991, MNRAS, 49, 43 Liner E.V.,1988, A&A, 6, 19 Liner E.V.,1998, ApJ, 497, 8 Nakamura T., 1997, Phys. Rev. D, 4, 5 Perlmutter S. et al., 1997, ApJ, 483, 565 Perlmutter S. et al., 1999, ApJ, 517, 565 Perrotta F., Baccigalupi C., Bartelmann M., De Zotti G., Granato G. L.,, MNRAS, 39, 445; astro-ph/1143 Riess et al., 1998, AJ, 116, 19 Riess, A.G.,, PASP, 11, 184 Refsal S., 1966, MNRAS, 13,15R Refsal S., 197, ApJ, 159, 357R Seit S.,Schneier P., Ehlers J., 1994, Class.Quant.Grav Ronveau, A., Heun s Differential Equation, 1995, Ofor University Press, Ofor Sachs R.K., Kristian, J., 1966, ApJ, 143, 379 Schneier P., Weiss A., 1988, ApJ,37, 56 Schneier P., Ehlers J., Falco E.E., 199, Gravitational lenses, Springer- Verlag, Berlin Tomita K., Futamase T., Gravitational lensing an the highreshift universe, 1999, Prog. Theor. Phys., 133 Tomita K., Premai P., Nakamura T., 1999 Prog. Theor. Phys., 133, 85 Tricomi F.G., Differential equations, 1961, Blackie, Lonon Watanabe K., Sasaki M., 199, Publ. Astron. Soc. Japan, 4, L33 Weinberg S., 1976, ApJ, 8, L1 Weinberg S., 1989, Rev. Mo. Phys., 61, 1 Zelovich, Ya., B., 1967, Pis ma Zh. Eksp. Teor. Fi. 6, 883; (JETP Lett. 6,316 (1967)).