# First-Order Light Deflection by Einstein-Strauss Vacuole Method

Authors: Debasish Saha; Amarjit Tamang; Ramil Izmailov; Carlo Cattani; Kamal K. Nandi
Journal: Journal of Gravity
 (2013)

---

## Abstract

We resolve here an outstanding problem plaguing conformal gravity in its
role in making consistent astrophysical predictions. Though its static
spherically symmetric solution incorporates all the successes of
Schwarzschild gravity, the fit to observed galactic rotation curves requiresγ>0, while the observed increase in the Schwarzschild light
deflection by galaxies appears to demand γ<0. Here we show that,
contrary to common knowledge, there is an increase in the
Schwarzschild deflection angle in the vicinity of galaxies due purely to the
effect of γ>0, when the idea of the Einstein-Strauss vacuole model
is employed. With the inconsistency now out of the way, conformal gravity should be
regarded as a good theory explaining light deflection by galaxies.


---

## Body

## 1. Introduction

The metric exterior to a static spherically symmetric distribution in Weyl conformal gravity has been obtained by Mannheim and Kazanas [1]. Recently, the solution has been used to fit rotation curves of many galaxy samples [2] as well as to predict the maximal size of galaxies [3]. The metric, which we call Mannheim-Kazanas-de Sitter (MKdS) metric, reads (G=c=1)
(1)dτ2=-B(r)dt2+1B(r)dr2+r2(dθ2+sin2θdφ2),(2)B(r)=1-2Mr+γr-kr2,
where M is the central mass and k and γ are arbitrary constants that could be appropriately fixed by using the fit to rotation curves. For distances neither too small nor too large, the above mentioned metric is a good approximation. Now, there could be three possible ways to calculate light deflection in the above spacetime. First, the conventional calculations for light deflection show that the constant k does not appear in the relevant equations, leading finally to the two way deflection as [4]
(3)2ϵ=4Mr0-γr0,
where r0 is the distance of closest approach. The difficulty is that the fit to observed rotation curve requires γ>0, and for consistency all other astrophysical observations should respect this sign. Now, the observed light deflection by a galaxy is always more than the Schwarzschild value 4M/r0, and hence to avoid the negative contribution in (3), one must demand γ<0 [4]. Thus there appears an inconsistency from the usual method. The second option is to use the Rindler-Ishak method [5], which is based on the realization that conventional methods do not apply to asymptotically nonflat spacetimes as the limit r→∞ makes no sense in it. Their original method of invariant angle is most appropriate in such situations, but it has an as yet unnoticed difficulty on the galactic scales, as explained below. Ishak et al. [6] thereafter improved the calculations using the Einstein-Strauss vacuole model and this provides us with the third and best option in our opinion. The purpose of this paper is to use the vacuole model to show that there is an increase in the Schwarzschild deflection angle in the vicinity of galaxies due purely to the effect of γ>0.


## 2. Rindler-Ishak Original Method

The bending angle in general is defined byϵ=ψ-φ. Rindler and Ishak considered the case φ=0 so that the deflection angle is ϵ=ψ given by
(4)tanψ=rB(r)|A|,
where A(r,φ)=dr/dφ. With u=1/r, the photon trajectory from (2) is given by
(5)d2udφ2=-u+3Mu2-γ2.
As evident, k has disappeared from the above equation. Solving perturbatively to first order in M, we have
(6)u=1r=sinφR-γ2+M4R2[6+3R2γ2-3Rγ(π-2φ)cosφnnnnnnnnnnnnnnnnnnnn+2cos2φ-6RγsinφR2γ2],
where the parameter R is related to the closest approach distance r0 at φ=π/2. Note that the usual Schwarzschild orbit equation u=(1/R)sinφ+(M/2R2)(3+cos2φ) is obtained at γ=0. From (6), we can find r at φ=0 as
(7)r=4R2X,
where
(8)X≡8M-3πRMγ-2R2γ+3MR2γ2.
Also, at φ=0, we find that
(9)|A|=16R3X2.
It can be seen again that, at γ=0, we recover the Schwarzschild values r=R2/2M and |A|=R3/4M2. Using the value of r from (7) and |A| from (9), we get from (4) the required one way deflection angle for small ψ in the first order in γ as
(10)ψ≃rB(r)|A|=2MRY+f(M,R,Y)γ,Y≡1-4M2R2-kR44M2.
For the solar system, M⊙=1.48×105 cm, R⊙=6.96×1010 cm, and taking k=9.54×10-54 cm-2 [2], we get kR4/4M2=1.15×10-22, which is far less than unity allowing a series expansion of Y. The expansion up to correct second-order contribution also follows. One also obtains a positive contribution to deflection due to the Mannheim value of γ (= 3.06×10-30 cm), namely, +13M2γ/R≃1.25×10-29, a very insignificant effect. So there is no difficulty in the local scale; the difficulty arises when one goes to the galactic scale. Typically, say for UGC0128, the disc plus bulge plus hydrogen gas mass M=5.16×1015 cm, and the scale length R0=2.07×1022 cm [2], we have kR04/4M2=7.40, which makes ψ imaginary and the expansion terms from Y unreliable! Thus the method needed improvement and that is exactly what was done by Ishak et al. [6] in their subsequent implementation of the Einstein-Strauss vacuole. We shall follow this below.


## 3. Einstein-Strauss Vacuole: Improved Rindler-Ishak Method

The suggestion of a cut-off boundary was first made by Edery and Paranjape [4]. We shall assume that the spacetime (2) (or MKdS vacuole) is exactly embedded in a FLRW universe, and the matching between the two spacetimes occurs at a radius rb, the boundary of the vacuole. We shall further assume that light deflection occurs only inside the MKdS vacuole and that once the light transitions are out of the vacuole and into FLRW universe, light ceases to bend. The junction conditions are prescribed by Einstein and Strauss [7] as
(11)rbinMKdS=a(t)rbinFLRW,MMKdS=4π3rbinMKdS3×ρmatterinFLRW.
Note that conformal gravity was used to fit galactic rotation curves without assuming any hypothetical dark matter, and the global positive constants γ and k can be assigned suitable numerical values to obtain the required fit [2]. Therefore, by MMKdS, we mean only the luminous mass contributed by the bulge, disc, and hydrogen gas inside the vacuole. In general, the vacuole radius rb would also change due to cosmic expansion, but we shall consider rb at that particular instant t0 of cosmic epoch when the light ray happens to pass the point of the closest approach r0 to the lens. Current observations suggest that the universe is spatially flat so that we can take ρmatterinFLRW=ρcritical=3H02/8π=1.1×10-29(H0/75km/sec/Mpc)2 gm/cm3. The bending is not sensitive to the exact value of ρ. Normalizing the scale factor to a(t0)=1 and dropping suffixes, the above prescription translates to
(12)rb=(3M4πρ)1/3,whereM is the luminous lens mass expressed in units of sun’s mass M⊙=1.989×1033 gm. Clearly, for any given lens, rb is much larger than the closest approach distance r0≃R, often expressed in terms of impact parameter.Assume a small entry angleφb radian at the vacuole boundary radius rb such that sinφb≃φb and cosφb≃1. Then, (6) gives
(13)1rb=-γ2+φbR+M[2R2-3πγR+3γ24].
Similarly,
(14)Ab(rb,φb)=rb2R-2Mφbrb2R2.
The small one way bending angle ψ becomes
(15)ψ=rbB(rb)|Ab|-φb,
which expands in first order in γ as
(16)ψ=-φb+(2M+RφbR-2Mφb)Z+g(M,R,φb,Z)γ≡ψ0+gγ,
where
(17)Z≡1-4M2R2-2MφbR-kR4(2M+Rφb)2,(18)ψ0≡-φb+(2M+RφbR-2Mφb)Z.
Given any M in grams of a galaxy or cluster lens, one can easily calculate its rb using (12). Then using the value of scale length R of the lens, and putting the value of rb in the equation (rewritten from (13))
(19)φb=R[1rb+γ2-M(2R2-3πγ4R)],
it can be verified that φb is indeed very small justifying its small approximation at the vacuole boundary entry. Furthermore, the term involving k in (17) is also much smaller than unity thus allowing for a legitimate expansion of Z. The difficulty involving the k dependent term mentioned in Section 2 now disappears, as the numerical calculations with typical galactic values show. We also note that for any galaxy lens the contribution due to gγ is negative but extremely small and is of no interest. Expanding the remaining term ψ0 in (18), we get
(20)ψ0≃2MR+M(φb2R+kRφb2-kR)-kR22φb.
With φbrb≃R, the last term can be recognized as the term -ΛRrb/6 obtained in [6] contributing negatively. The middle term coupled to mass M and γ (via φb) is larger than the second-order Schwarzschild term (15πM2/8R2) with the leading order term Mφb2/R contributing positively to the Schwarzschild term 2M/R. This contribution is independent of k and purely due to the other constant γ. However, note that this additional term is always positive, no matter what the sign of γ is but the value is slightly more when γ is taken with a positive sign, especially when the impact parameter R≫R0 inside the vacuole, which indicates that a positive sign γ is preferable. The other contributions proportional to negative k in (20) are much less and negligible for typical galactic values. The fact that the vacuole model provides a positive contribution to bending proportional to γ is the main result of this paper.


## 4. Conclusion

The present paper has been inspired by the fact that, in recent times, MKdS solution considered by Weyl conformal gravity has seen remarkable successes in explaining the flat rotation curves and maximal sizes of galaxies without requiring the hypothesis of dark matter. However, with regard to light bending by the galaxies, when light has to traverse attractive dark matter, one would expect that light should bend more than the Schwarzschild value. This increase has actually been confirmed by experimental observations on gravitational lensing by galactic clusters [4]. The difficulty with the standard method is that it is not suitable for asymptotically nonflat spacetimes such as MKdS. Also, the original Rindler-Ishak method leads to imaginary values when applied to galactic scales. Thus, the improved method developed by them (using the concept of Einstein-Strauss vacuole) shows that the contribution due to γ term is always positive thereby making the total first-order contribution in excess of the Schwarzschild value, as desired. The result of this paper places conformal gravity on a firmer footing with regard to astrophysical observations.Finally, we should mention that the above result is supported also by a similar result in Cattani et al. [8] about the first-order bending being in excess of the Schwarzschild value. Their method is based on the original Rindler-Ishak method as developed in [5] appropriate to asymptotically nonflat situations. However, Sultana and Kazanas [9] used the standard method, which seems rather inappropriate in such situations. Moreover, their actual calculations are disfigured by certain errors that lead to a decrease from the Schwarzschild value, which is unphysical as the galactic halo is attractive.



---

*Source: 686950-2013-08-01.xml*