SciELO - Scientific Electronic Library Online

 
vol.53 número2LX Leo: A High Mass-Ratio Totally Eclipsing W-type W UMa System índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de astronomía y astrofísica

versión impresa ISSN 0185-1101

Rev. mex. astron. astrofis vol.53 no.2 Ciudad de México oct. 2017  Epub 21-Oct-2019

 

Artículos

Elements of a Gravitational Lens by Assuming an Elliptical Galaxy Model

U. E. Molina1 

P. Viloria2 

I. Steffanell3 

1 Departamento de Física, Universidad del Atlántico, Barranquilla, Colombia (ubaldomolina@mail.unialtlantico.edu.co).

2 Departamento de Ciencias Naturales y Exactas, Universidad de la Costa, CUC, Barranquilla, Colombia (pviloria@cuc.edu.co).

3 Facultad de Ingeniería Universidad Libre, Barranquilla, Colombia (isteffanell@unilibrebaq.edu.co).


ABSTRACT

In this work we study elements for gravitational lenses per galaxies such as the equation of lens, deviation angle, deflection potential and time delay, modeling the distribution of the volumetric mass of the lens in elliptical form. The function of volumetric distribution of mass in the deflecting galaxy ρ has a nucleus with radius a in its core, form-free parameter b (b > a), and volumetric density in its nucleus ρ 0 . Through the distribution of volumetric mass density ρ, we initially find surface mass density Σ (projected on plane of the lens), followed by elements of a gravitational lens which are completely general and expressed in terms of the geometric parameters a and b. These are related by the adimensional factor n= bα >1. Results are applied to a galaxy-specific lens system to conduct an analysis based on the temporal delay between two images and to observe the conditions with which parameter n must comply.

Key Words: Gravitational lens; the lens equation; deflection angle; deflecting potential; time delay

RESUMEN

En éste trabajo se estudian algunos elementos para lentes gravitacionales por galaxias como son ecuación de la lente, ángulo de desviación, potencial de desviación y retardo temporal, modelando la distribución volumétrica de masa de la lente en forma elíptica. La función de distribución volumétrica de masa en la galaxia deflectora ρ, contiene en su centro un núcleo de radio a, una densidad volumétrica en su núcleo ρ 0 and un parámetro libre de forma b, (b > a). Mediante la distribución de densidad volumétrica de masa se encuentran inicialmente la densidad superficial de masa Σ, (proyectada en el plano de la lente), para después hallar los elementos de una lente gravitacional los cuales son totalmente generales y quedan en términos de los parámetros geométricos α and b, que se han relacionado mediante un factor adimensional n= bα >1. Los resultados encontrados se aplican a un sistema de lente por galaxia particular para hacer un análisis basado en el retardo temporal entre dos imágenes y ver las condiciones que debe cumplir el parámetro n.

1. INTRODUCTION

In the study of gravitational lenses (GL) volumetric distribution of mass density of the deflector (ρ), can be projected onto a plane perpendicular to the line of vision between the observer and the light source. This plane is called the lens plane, according to Narayan, R. and Bartelmann, M. (1997). The lens is considered thin and the volumetric distribution of lens mass is replaced by a plane on which the surface mass density is (Σ), constituting the so-called approximation of the thin lens, in accordance with Schneider, P. et al. (1992).

In general, GL require a few basic elements such as: surface mass density (Σ), lens equation, deviation angle (α), deflection potential (Ψ) and time delay (∆t), which comprise basic tools for their application in the study of some lens systems in astrophysics. These analytical expressions are specified when applied to a model of mass volume distribution in a particular galaxy.

Some observational data shows that astronomical objects that act as lenses are modeled in diverse ways, according to Cohen, A. S. and Hewitt, J. N. (2000). In systems of gravitational lenses per galaxy, we observationally mea- sure parameters such as: dispersion speed of the material particles that make up the lens (σ p ), angular position of the images (θ), red shifts of lens (z L ), source (z S ) and time delay (∆t); with these it becomes possible to study and assume models in said systems.

The distances from the source to the lens and from the lens to the observer are ∼ 1pc. In this way, systems that constitute the lenses, such as deflecting galaxy, source and observer, are too far away, which causes light to travel in free space most of the time, this light only locally becoming deviated when it passes through the lens. We must therefore model the universe through which the ray of light should pass; to do this we need to set cosmological parameters, Used, for example, by the authors Adler et al. (1975), Ciufolini et al. (1995), Foster et al. (1994), Kenion (1995): vacuum density (Ω v ); matter density (Ω m ); softness parameter ( α~) and Hubble constant (H 0 ).

By applying the properties of GL to a specific system, using observational values and setting cosmological parameters, it is possible to study different galactic models.

In our case we will assume a model of a galaxy which is elliptical in its volumetric distribution of mass (ρ), with which we determine the analytical expressions of GL, used, e.g., by Brainerd et al. (1996) and Golse et al. (2002). Here we follow Hjorth et al. (1997) and Molina et al. (2006), who propose a distribution of volumetric mass density useful for models of elliptical galaxies acting as gravitational lenses.

2. LENS ELEMENTS AND ELLIPTIC MODEL IN THE GALAXY DEFLECTOR

2.1 Lens elements

In GL literature, the approximation of a flat lens is characterized by a surface mass density given by the projection operator, according to Miranda, C., Molina, U and Viloria, P. (2014),

Σ(R) =ρ(R , z)dz (1)

where R is a radius vector in the lens plane as shown in Figure 1, and ρ is the volumetric distribution of the lens mass. The radius vector R, called impact parameter, can be written R = ξ0x. The quantity ξ0 is known as scale parameter or scale factor, and is defined according to the lens model being used. In our case it is defined below in equation (10). According to Schneider et al., 1992 (p. 231), matter within the disc of radius around the mass contributes to the deflection of the ray of light, and matter outside of this disc (x’ < x) does not contribute importantly to deflection, in such a way that the skew angle follows the expression,

α(x) =2xx'κ(x')dx' (2)

Fig. 1 Illustration of a gravitational lens system. The angular separations of the source and the image from the optic axis as seen by the observer are β and θ, respectively. The angular diameter distances between the observer and the source, the observer and the lens, and the lens and the source are DS , DL and DLS , respectively, by Narayan et al. (1997)  

when x’ < x. The quantity κ defined as κ(x) = Σ(x)cr is the so-called convergence which indicates the existence of a minimal or critical surface density for the GL phenomenon to occur, in accordance with Narayan et al. (1997) and Schneider et al. (1992). We define surface density of critical mass as: Σcr = c2 Ds4. π.G.DLDLS, where D S , D L and D LS , are the angular diameter distances between observer-source, observer-lens and lens-source, respectively, as shown in Figure 1.

The deflection potential of the lens is according to Narayan et al. (1997),

ψ =1πκ(x') ln |x - x'| dx'2 (3)

The time delay ∆t between two light beams detected by an observer is given by Molina et al. (2006) and Narayan et al. (1997),

Δt=1+zLc ξ02DsDLDLS12α22-α12-ψ2-ψ1(4)

where z L is the redshift of the deflecting galaxy, D L is the angular diameter distance of observer-lens, D S is the angular diameter distance of observer- source, D LS is the angular diameter distance of lens-source.

The expression (4), contains the geometric delay described in equation (2), and the gravitational potential given by equation (3).

Furthermore, the relationship between source position (β), positions of images ( θ =ξDL) and deflection angle can be written in accordance with equation (2.15a) in Schneider et al., 1992 (p.31) as,

β=ξ0DLx-DLSDSαx (5)

which is called the lens equation. Equations (2), (3), (4) and (5) make up the group of basic elements for a GL study.

2.2 Elliptical galaxy model deflector

For our study, we modeled the distribution of lens mass as an elliptical galaxy, following Hjorth and Kneib (1997), who proposed a distribution of volumetric density which is useful for elliptical galaxy models acting as gravi- tational lenses. With these distributions we find the analytical expressions of GL. This volumetric distribution of mass is,

ρr=ρ01+r2α21+r2b2 (6)

The model of deflecting galaxy contains a central nucleus with radius a, a free-form parameter acting as scale b (b > a), and a volumetric density of fixed mass in the nucleus ρ0. Introducing the volumetric density given in equation (6), in the projection operator defined in equation (1) and making the change of variable, z 2 = r 2R 2 (see Fig. 1), we obtain the surface mass density of the lens. The distance z is typically much smaller than the distances between observer and lens and between lens and source, so after developing the integrals the surface mass density of the lens takes the following form,

n, R=0n2αn2-11α2+R2-1n2α2+R2 (7)

where n=bα>1 is the adimensional parameter and Σ0 = πρ0 a is the surface density contained in the nucleus.

To find the analytical expressions of the deviation angle and the deflection potential, first we change R = ξ0x, in equation (7), in such a way that the surface density of the flat lens is,

n, x=0n2αn2-1ξ01A2+x2-1n2A2+x2 (8)

where A=αξ0 a new parameter. The convergence factor k=ΣΣcr (defined above), now has the form,

kn,x=121A2+x2-1n2A2+x2 (9)

where we choose a scale factor such as,

ξ0=8πGDLDLSΣ0n2αc2DSn2-1 (10)

The scale factor is set by the dispersion speed of the components of the deflecting galaxy σ p (this is explained in the following paragraphs). We observe that the scale factor is expressed in terms of the central radius a, and the adimensional parameter n=bα>1. We see that when the adimensional parameter approaches one, that is, n ≈ 1, the scale factor becomes infinite, which prompts one to think that it would be preferable choose a scale that depends only on the radius of the nucleus a.

In the time delay expressed in equation (4), we see the arbitrary scale factor ξ 0 , which we choose for the distribution of elliptical mass as shown in equation (10). The central volumetric density ρ 0 is not a quantity that is observationally measurable in a gravitational lens system, but it is possible to establish it by knowing the rate of dispersion σ p of the matter components of the deflecting galaxy. To do this, through the expression proposed by Hjorth and Kneib (1997),

σρ2=2GΣ(R)RMrr2ρ(r)r2-R2dr (11)

it is possible to obtain the dispersion rate.

In this expression (11) we see that if the dispersion velocity of the matter components of the deflecting galaxy is known, the central volumetric density can be set, and with this, the scale factor (10) is also set.

2.3 Deviation angle

Expressed in terms of impact parameter R = ξ 0 x and using the convergence factor in equation (9), the deviation angle (2) takes the form,

αn,R=n-1αR+1+α2R2- 1+n2α2R2 (12)

which is given in terms of central radius a and adimensional parameter n =bα>1. We observe that if it is the case that this parameter is approximate to one, that is, n ≈ 1 and selecting a set scale, the deviation angle becomes α(R) ≈ 0.

2.4 Deflection Potential

We find the deflection potential of the lens that we are considering by re- placing equation (9) in equation (3) and writing in terms of impact parameter R = ξ 0 x, we obtain,

ψn,R=αξ0n-1+ln2αξ0-nln2nαξ0+1+R2α2 -n2+R2α2-lnα1+1+R2α2ξ0+nlnαn+n2R2α2ξ0 (13)

given that n > 1. We note that if this adimensional parameter approximates one, that is, n ≈ 1, and if we select a set scale, the deflection potential becomes zero, ψ(R) ≈ 0.

2.5 Time delay

Differentiating the square of the deviation angle for two images, α22-α12 and subsequently differentiating the deflection potential ψ 2ψ 1, and then replacing these two differentiations in equation (4), we get the time delay as,

Δt=1+zLξ02DScDLDLSh(n,α)2-αg(n,α)ξ0 (14)

Equation (13) allows us to determine the time delay between two images, where a the radius of the nucleus of lens mass distribution is and n is the adimensional parameter.

Furthermore, we have defined two new functions: h(n, a) and g(n, a), de- fined below, depend on how the nucleus radius a is set and and on how much these two new functions vary from the adimensional parameter n.

Function h(n, a) is defined in the form,

hn,α=n2-n+1α2R22-α2R12

+(n-1)α2R22+α4R24-α2R22+n2α4R24

+(n-1)α2R12+n2α4R14-α2R12+α4R14

-1+α2R22+n2α2R22+α4R24

+1+α2R12+n2α2R12+α4R14 (15)

Similarly, g(n, a) is defined as,

gn,α= 1+R22α2-1+R12α2

+n2+R12α2-n2+R22α2

+ln1+1+R12α21+1+R22α2

+nlnn+n2+R22α21+n2+R12α2

h(n, a) and g(n, a), are normalized so that when they are introduced into the time delay equation (14), their units are expressed in seconds.

2.6 Lens Equation

By substituting the deflection angle (12) in equation (5), we get,

β=RDL-DLSDS n-1αR-DLSDS[1+α2R2-1+n2α2R2](17)

This formula represents the equation of the lens in the elliptical lens model we propose, which is a function of the impact parameter and adimensional parameter n.

3. APPLICATION OF THE PROPOSED GRAVITATIONAL LENS MODEL TO GALACTIC LENS B0218 + 357

The expressions articulated in the previous section are completely general and can be applied to any gravitational lens system. In our case, by way of example, we choose to apply them to the B0218 + 357 lens system, to analyze the effectiveness of the results we obtain.

Some researchers, such as, notably, Wucknitz et al. (2004) who, in their work Models for the Lens and Source of B0218 + 357 determine the Hubble constant H0, and discuss different models for the B0218 + 357 galactic lens. In addition to this, in order to obtain an estimate of the Hubble constant,

A.D. Biggs et al. (1999) model the B0218 + 357 system using the lens model described by Kormann et al. (1994) as a Singular Isothermal Ellipsoid (SIE) mass. More information on the morphology of the received images can be found in the work of C. Spingola et al. (2015). Observational data for the B0218 + 357 lens system deposited in the CASTLES Survey, according to Cohen, A. S., and Hewitt, J. N., (2000) is summarized in table 1 below: the dispersion speed σp, the angular positions of two images θ 1 and θ 2 , the redshifts for lens z L and the source z S ; and the difference in time delay between the two images ∆t, as follows:

Table 1 Data observed in B0218 +357 lens system 

σ p θ 1 θ 2 z L z S t
150km/s 40mas 290mas 0.68 0.94 10.5±0.2days

For our case, we choose cosmological parameters within the range of values most accepted in literature; we thus rely on the work of several authors, and among these point especially to: Kessler et al. (2009), Boughn and Crittenden (2001), Bartelmann et al. (1997), Weinberg (1972), Grogin et al. (1996)) and others. The parameters we choose for our elliptical lens model are: Hubble Constant, H0 = 76 km/spc, Vacuum density Ω v = 0.7, matter density Ω d = 0.3, and softness parameter α~ = 0.5. This softness parameter of matter in the universe is smoothly distributed (that is, it is not bound up in galaxies), note Dyer and Roeder (1973), see also P. Schneider et al. (1992) (p.138) and Xi Yang et al. (2013).

According to Dyer, C. C. (1973) and Wucknitz et al. (2004), by using, for the B0218 + 357 galactic lens system, the previous cosmological parameters and the following values: cross section DLDLSDs, observer-lens distance D L , impact parameters of the two images R 1 and R 2 , we obtain: DLDLSDs=2.25×107pc, DL=1.364×109pc, R1=264.58pc and R2=1918.21pc.

With these quantities, we can find the elements of the proposed lens model.

3.1. Scale factor, surface mass density of the lens, deviation angle and deflection potential.

Because there are two images in the B0218 + 357 lens system, the previous calculations on the two impact parameters R 1 = 264·58pc and R 2 = 1918·21pc allow us to set the radius of the nucleus in the range of values of these two parameters (R 1 and R 2 ). In this work, we give the radius of the nucleus ap- proximately the value of the minor impact parameter, expressing the radius of the nucleus approximately as a = 264pc, which corresponds to the impact parameter with least deviation. To facilitate analysis we define the adimensional quantity λ as λ = R α> 1; by setting the nucleus radius and varying the impact parameter according to this nucleus radius, we can also set the surface density in the nucleus as the approximate value of Σ0 = 58.62kg/m 2 .

Furthermore, the scale factor in equation (10), is in terms of the adimen- sional parameter n, that is,

ξ0 = (200·16pc)n2n2-1 (18)

being pc = 3.086 × 1016 m. The surface mass density of the lens, expressed in equation (8) is then found to be,

Σn, λ=58.62n2n2-111+ λ2-1n2+λ2kg/m2 (19)

where for this particular lens, the impact parameter satisfies the condition 1 ≤ λ ≤ 7.25 and n > 1. When the value of n is fixed, expression (19) allows us to estimate the surface mass density of the lens as a function of λ, in the given interval.

After making the replacements required by our proposed system, the de- viation angle represented in equation (12), takes the form,

αn,λ=n-1λ+ 1+1λ2-1+n2λ2 (20)

in which we know that 1 ≤ λ ≤ 7.25 and n > 1.

At the same time, the deflection potential expressed in equation (13), takes the new form,

ψn, λ=αξ0[n-1+ln2αξ0-nln2nαξ0+1+λ2-n2+λ2-lnα1+1+λ2ξ0+nlnαn+n2+λ2ξ0] (21)

which depends upon the scale factor represented in equation (18) and the following conditions: 1 ≤ λ ≤ 7.25 and n > 1.

3.2. Time Delay Model

By using the values obtained from the B0218 + 357 lens system, we reduce the time delay stated in equation (14) is reduced to,

Δt=5.68daysn2n2-12hn2-1.32n2-1gnn2 (22)

Quantities h and g, and equations (15) and (16), depend only upon the adimensional parameter n. This allows us to establish values for the time delay. We take advantage of the fact that the time delay between the two images is measured observationally, as shown in table 1 (10.5 days), and a series of different values of n are explored in equation (22), until observed value is reached. The time delay of equation (22), for some values of n, is as follows:

When we compare the observational time delay ∆t = 10.5 ± 0.2days, as shown in table 1, with equation (22), we see that for the B0218 + 357 lens system, the adimensional parameter has a range of certainty of 2 ≤ n ≤ 2.2. We can appreciate the most approximate value to n = 2.1 in Table 2.

Table 2 Time delay values, for some values of n, in the B0218 +357 lens system, according to equation (22). 

n 2 2.1 2.2
t(days) 9.67 10.3 10.9

Establishing the value of parameter n allowed us to estimate the geometric parameters of the elliptical lens, a and b, and afterwards find the numerical values for the other basic elements such as the surface density of the lens, the deviation angle, the scale factor and the deflection potential, as indicated in the following section.

3.3. Estimation of deviation angle, deflection potential and lens equation

Given that for the B0218 + 357 lens system, the approximate radius of the nucleus is, the adimensional parameter is a = 264pc, the cross section is DLDLSDL and the observer-lens distance is D L , we can estimate the deviation angle, deflection potential and the lens equation in this model:

  1. Using as our premise equation (19), taking into account that the value of the adimensional parameter is n = 2.1, the surface density diminishes as the impact parameter increases, and because this is in the interval 1 ≤ λ ≤ 7.25, the surface mass density of the lens for the proposed system is estimated in the range of 0.32kg/m 2 ≤ Σ ≤ 32kg/m 2 .

  2. In accordance with the equation (20), and taking into account that the parameter λ is in the range 1 ≤ λ ≤ 7.25, we deduce that deviation angle oscillates between 120mas and 220mas.

  3. In accordance with equation (18), and with the values of adimensional parameter n, represented in section 3.1, we can establish an approximate value of the scale factor, that is, ξ 0 = 260pc. Furthermore, as the values of parameter λ oscillate between 1 ≤ λ ≤ 7.25, then the values of the deflection potential, according to equation (21), oscillate between 150mas and 930mas.

4. CONCLUSIONS

This work is based on a volumetric mass distribution which describes a fast relaxation scenario, similar to a model of an SIS isothermal sphere, according to Kenion, I. R. (1995), where the mass of the lens is considered spherically symmetric. The model of elliptical density contains a central nucleus with radius a, mass density ρ 0 in the central nucleus and free-form parameter a, a mass density in the central core ρ 0 and also a free shape parameter b (b > a). From this volumetric mass distribution we find new analytical expressions of the lens elements. These new elements are equations for: lens surface density, deviation angle, deflection potential and time delay. These expressions depend upon the impact parameter of the images and on geometric lens elements a and b, related by the adimensional parameter n = bα > 1.

The analytical expressions of the surface mass density of the lens, deviation angle, deflection potential and time delay used in this work to describe our proposed model can be used to analyze any other system of galaxy lens. Our equations are quite general and for their application to study a specific lens system only require the observational measurements indicated in Table 1.

To implement the results of section (2), one basically needs a system of lenses per galaxy whose mass density distribution fits the elliptical model we describe in this paper, the analytical expressions found here being a good point of departure for further research. In this work the proposed gravitational lens model is applied specifically to the B0218 + 357 lens system, using observational values shown in Table 1. To do this, we have used the cosmological parameters most widely accepted in literature, used, e.g., by authors Bartelmann, M., at al. (1997), Boughn et al. (2001), Foster et al. (1994), Kessler et al. (2009), Schneider et al. (1992),such as the Hubble constant, vacuum density, matter density and softness parameter contained in the angular dia- metric distances of observer-lens, observer-source and lens-source. This allows us to find the approximate radius of the nucleus of the deflecting galaxy, the angular diametrical distances, the scale factor synthesized in equation (18), and afterwards, adjust the adimensional parameter n.

The adimensional parameter n is adjusted through the time delay represented in equation (22), and compared to the observed time delay ∆t obs = 10.5days, shown in Table 1. To do this, it was necessary to write the impact parameter R, in terms of the nucleus radius ; in the form of λ = Rα  1.

Because the impact parameter is expressed in terms of the radius of the nucleus, we were able, once we have set cosmological parameters, determine the theoretical time delay given in equation (22), depending only on the adimensional parameter n. Thus, when we compare the theoretical time delay given in equation (22), with the observationally measured time delay in Table 1, we find that the range of values of the adimensional parameter n is 2 ≤ n ≤ 2.2.

Finally, as evidenced in section (3.3), by adopting the approximate value n = 2.1 for the adimensional parameter and the interval for the impact parameter 1 ≤ λ ≤ 7.25, we were able to estimate numerical values for the surface density of the lens, deviation angle, scale factor and deflection potential in our proposed system.

REFERENCES

Adler, R., Schiffer, M., and Bazin, J., 1975, Introduction to General Relativity, McGraw-Hill co, New YorK. [ Links ]

Bartelmann, M., et al., 1997, On Arcs and Omega. arXiv:astro-ph/9709229v1. [ Links ]

Biggs A. D., Browne I. W. A., Helbig P., Koopmans L. V. E., Wilkinson, P. N., and Perley R. A., 1999, MNRAS, 304, 349 [ Links ]

Boughn, S., and Crittenden, R., 2001, arxiv.org/pdf/0908.4274v1.pdf. [ Links ]

Brainerd, T. G., Blandford, R. D., and Smail, I., 1996, Weak Gravitational Lensing by Galaxies. Apj, 466. [ Links ]

Ciufolini, I., and Archibald, J., 1995, Gravitation and Inertia, Ed. Princeton Series in Physics, New Jersey. [ Links ]

Cohen, A. S., and Hewitt, J. N., 2000, Further Investigation of the Time Delay, Magnification Ratios, and Variability in the Gravitational Lens 0218+357. Apj, 545. [ Links ]

Dyer, C. C., and Roeder, R. C., 1973, The Distance-Redshift Relation for Universes with some Intergalactic Medium. Apj, 180, L31. [ Links ]

Foster, J., and Nightingale, J. D., 1994, A Short Course in General Relativity, Ed. Springer-Verlag, NY. [ Links ]

Golse, G., and Kneib, J., 2002, Pseudo Elliptical NFW Lensing Mass Model: application to the NFW mass distribution, Astronomy Astrophysics 390, 821-827. [ Links ]

Grogin, N. A., and Narayan, R., 1996, A New Model of the Gravitational Lens 0957+561 and a Limit on the Hubble Constant, Apj, 464. [ Links ]

Hjorth, J., andKneib, J. , 1997, Elliptical galaxies as gravitational lenses, P. Apj. [ Links ]

Kenion, I. R., 1995, General Relativity, Oxford Univrsity Press Inc., NY. [ Links ]

Kessler, R. et al., 2009, arxiv.org/pdf/0908.4274v1.pdf. [ Links ]

Kormann, R., Schneider, P., and Bartelmann, M. , 1994, Astron. Strophys. 284-299. [ Links ]

Meneghetti, M., Bartelmann, M. , and Moscardini, L., 2002, astro-ph/0201501. [ Links ]

Miranda, C., Molina, U and Viloria, P., 2014, Retardo temporal en las lentes por galaxias en el contexto de Reissner - Nordstrom., Revista Mexicana de Física 60 (2014). [ Links ]

Molina, U., Tejeiro, J., andViloria, P., 2006, Modelo Elíptico a partir del tiempo de retardo observacional., Rev. Colombiana de Física. Vol. 38, No.3. [ Links ]

Narayan, R., and Bartelmann, M. , 1997, Lectures on Gravitacional Lensing, Max Planck-Institut für Astrophisik. [ Links ]

Schneider, P. , Ehlers, J., and Falco. E., 1992, Gravitacional lenses, Ed. Springer- Verlag, NY. [ Links ]

Schneider, P. , Kochanek, C., and Wambsganss, J., 2006, Gravitational Lensing: Strong, weak and micro. Springer-Verlag, Berlin Heidelberg. [ Links ]

Spingola, C., Dallacasa, D., Orienti, M., Giroletti, M., McKean, J. P., Cheung, C.C., Hovatta, T., Ciprini, S., Ammando, F. D., Falco, E., Larsson, S., W. Max- Moerbeck, W., Ojha, R., Readhead, A. C., Richards, J. L., and Scargle, J., 2015, MNRAS 000, 19. [ Links ]

Weinberg, S., 1972, Gravitation and Comology, Wiley, New York. [ Links ]

Wucknitz, O., Biggs, A. D., and Browne, W. A., 2004, Models for the lens and source of B0218+357: a LENSCLEAN approach to determine H0. Mon. Not. R. Astron. Soc. 349, 1430. [ Links ]

Xi Yang., Hao-Ran Yu and Tong-Jie Zhang., 2013, Constraining smoothness parameter and the DD relation of Dyer-Roeder equation with supernovae, Astrophys. arXiv:1305.6989v2. [ Links ]

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License