Advanced Topics

Advanced Topics

Gravitational Lensing

Einstein’s theory of General Relativity predicts that massive objects can bend the path of light. This effect is known as gravitational lensing. A gravitational lens is matter that bends light from a distant source as it travels towards an observer. This can be a star, black hole, galaxy, or even a cluster of galaxies.

Gravitational lensing can be categorized into three main types:

  • Strong lensing is when there are easily visible distortions, such as the formation of Einstein rings, arcs, or multiple images of the same source.
  • Weak lensing does not produce arcs or multiple images, but only flattens the shape of the source.
  • Microlensing is when no distortion in shape can be seen, but the amount of light received from a background object changes.

The impact parameter ξ\xi is the shortest separation between the center of the lens and the path of a particular light ray. Here, L is the massive object (lens) at a distance DLD_L from the observer, which bends the light emitted from the source S at a distance DSD_S. The observed images are S1_1 and S2_2. Angles θ\theta, ϕ\phi, and β\beta are very small. The deflection is given by

ϕ=4GMc21ξ(5.5.1)\tag{5.5.1} \phi = \frac{4GM}{c^2} \frac{1}{\xi}

where MM is the mass of the lens.

Gravitational lensing by a massive object

From the diagram, we have that ξ=θDL\xi = \theta D_L, and

βDS=θDSϕ(DSDL) \beta D_S = \theta D_S - \phi (D_S - D_L)

Since ϕ=4GMc21θDL\phi = \frac{4GM}{c^2} \frac{1}{\theta D_L},

βDS=θDS4GMc21θDL(DSDL)(5.5.2)\tag{5.5.2} \beta D_S = \theta D_S - \frac{4GM}{c^2} \frac{1}{\theta D_L} (D_S - D_L)

When the source is perfectly aligned with the lens (β=0\beta = 0), a ring-like image is formed, called the Einstein ring, having angular radius (called the Einstein radius)

θE=4GMc2DSDLDSDL(5.5.3)\tag{5.5.3} \theta_E = \sqrt{\frac{4GM}{c^2} \frac{D_S - D_L}{D_S D_L}}
Einstein Ring
An example of an Einstein ring formed by gravitational lensing.

Using θE\theta_E, equation (5.5.2) can be rewritten as

β=θ[1θE2θ2](5.5.4)\tag{5.5.4} \beta = \theta \left[ 1 - \frac{\theta_E^2}{\theta^2} \right]

which gives two solutions for θ\theta, corresponding to the two images formed

θ=β2[1±1+4θE2β2](5.5.5)\tag{5.5.5} \theta = \frac{\beta}{2} \left[1 \pm \sqrt{1 + \frac{4 \theta_E^2}{\beta^2}} \right]

A property of gravitational lensing is that the surface brightness of the images formed remains the same as that of the original object. Thus the apparent brightness of the image is directly proportional to its area. From this, we can find out the amplification μ\mu of light due to gravitational lensing

μ=θδθβδβ \mu = \frac{\theta \delta \theta}{\beta \delta \beta}

Evaluating this, we get

μ=1(θEθ)41(5.5.6)\tag{5.5.6} \mu = \left| 1 - \left( \frac{\theta_E}{\theta} \right)^4 \right|^{-1}

Cosmic Distance Ladder

There are several methods to find distances to celestial objects. Closer objects (moons, planets, sun) require a different technique to find their distance than distant objects like quasars. The cosmic distance ladder is a succession of methods by which astronomers determine the distances to celestial objects. Each method is built on the previous one, allowing for a more accurate measurement of distance.

Trigonometric Parallax

This is based on the apparent yearly back and forth movement of stars in the sky, caused by the orbital motion of Earth. This can measure distances up to ~30pc.

Statistical Parallax

The sun moves through the LSR with a speed of 20 km/s. Combining the apparent proper motion with observed radial motion allows the estimation of distance to a star cluster. This has been used to find the distance to Hyades.

Main Sequence Fitting

If the distance to a cluster is known, its HR diagram can be plotted. Another cluster whose distance is to be determined can be plotted on the same HR diagram and compared to get the distance modulus.

Standard Candles

Standard candles are classes of objects that have a calibrated absolute magnitude MVM_V.

Type Ia supernova is a type of supernova that occurs in binary systems in which one of the stars is a white dwarf. It occurs when the white dwarf exceeds the Chandrasekhar limit. The peak absolute magnitude of every Type Ia supernova is the same, MV=19.3M_V = -19.3.

Cepheid variables display a unique property known as the period-luminosity relationship. All cepheids exhibit a linear relationship between their luminosity and period. Another class of stars, RR Lyrae, also exhibit this property. The absolute magnitude of RR Lyrae stars is MV=0.6M_V = 0.6.

Hubble’s Law

Edwin Hubble discovered that galaxies move away from us with velocities proportional to their separation, i.e. that the universe is expanding. In terms of the redshift z=Δλ/λz = \Delta \lambda / \lambda,

v=cz=H0r(5.5.7)\tag{5.5.7} v = cz = H_0 r

where H0=70±7H_0 = 70 \pm 7 km/s/Mpc is the Hubble constant, and rr is the distance to the galaxy.

Apparent Superluminal Motion

Superluminal motion is the apparently faster-than-light motion seen in some radio galaxies. Bursts of energy moving out along the relativistic jets emitted from these objects can have a proper motion that appears greater than the speed of light.

It occurs usually when the source itself is moving at relativistic speeds, in a direction close to the line of sight.

Superluminal motion is most often observed in two opposing jets emanating from the core of a star or black hole. In this case, one jet is moving away from and one towards the Earth. If Doppler shifts are observed in both sources, the velocity and the distance can be determined independently of other observations.

Consider at time t1t_1, the source is at AA and at time t2=t1+δtt_2 = t_1 + \delta t it moves to BB. Then the photons emitted at AA at time t1t_1 are detected by the observer at

t1=t1+DL+vδtcosθc(5.5.8)\tag{5.5.8} t_1' = t_1 + \frac{D_L + v \delta t \cos \theta}{c}

where DLD_L is the distance from the source to the observer.

Superluminal motion is the apparent faster-than-light motion observed in some radio galaxies.

Photons emitted at BB at time t2t_2 are detected at

t2=t2+(DL2+v2δt2cos2θ)1/2ct2+DLc(5.5.9)\tag{5.5.9} t_2' = t_2 + \frac{(D_L^2 + v^2 \delta t^2 \cos^2 \theta)^{1/2}}{c} \approx t_2 + \frac{D_L}{c}

where vδtDLv \delta t \ll D_L. Therefore the time difference between the two photons is

Δt=t2t1=δt(1βcosθ)(5.5.10)\tag{5.5.10} \Delta t = t_2' - t_1' = \delta t (1 - \beta \cos \theta)

where β=v/c\beta = v/c. The angle θ\theta is the angle between the line of sight and the direction of motion of the source.

During the time Δt\Delta t, the radio source will move by an angle Δϕ\Delta \phi. Hence the proper motion μ\mu is

μ=ΔϕΔt=vsinθd(1βcosθ)(5.5.11)\tag{5.5.11} \mu = \frac{\Delta \phi}{\Delta t} = \frac{v \sin \theta}{d(1 - \beta \cos \theta)}

This gives the transverse velocity

βt=μdc=βsinθ1βcosθ(5.5.12)\tag{5.5.12} \beta_t = \frac{\mu d}{c} = \frac{\beta \sin \theta}{1 - \beta \cos \theta}

The maximum value of βt\beta_t for a given β\beta occurs when θ=cos1β\theta = \cos^{-1} \beta.

βtmax=β1β2=βγ(5.5.13)\tag{5.5.13} \beta_t^\text{max} = \frac{\beta}{\sqrt{1 - \beta^2}} = \beta \gamma

If γ1\gamma \gg 1 (i.e. when velocity of jet is close to the velocity of light) then βtmax>1 \beta_t^\text{max} > 1 , means that the apparent transverse velocity along CB is larger than the velocity of light in vacuum, i.e. the motion is apparently superluminal.

Recombination and Decoupling

Recombination refers to the process in which free electrons combine with protons to form neutral hydrogen atoms. The epoch of recombination is the time at which the baryonic component of the universe goes from being ionized to neutral.

Decoupling refers to the process by which photons cease to interact significantly with matter. The epoch of photon decoupling is the time at which the rate at which the photons scatter from electrons becomes less than the Hubble parameter. This is when the universe becomes transparent to light. The epoch of last scattering is the time at which a typical CMB photon underwent its last scattering from an electron.

In the early universe, atoms and ions were in statistical equilibrium with each other

H+γp+eH + \gamma \leftrightharpoons p + e

while the photons are still coupled to the baryonic component. In statistical equilibrium at temperature TT, the number density nXn_X of particles with mass mXm_X is given by the Boltzmann distribution

nX=gX(mXkT2π2)3/2exp(mXc2kT)(5.5.14)\tag{5.5.14} n_X = g_X \left( \frac{m_X kT}{2 \pi \hbar^2} \right)^{3/2} \, \exp(-\frac{m_X c^2}{kT})

when kTmXc2kT \ll m_X c^2. Here gXg_X is the statistical weight of the particle XX. For electrons, protons, and neutrons, g=2g = 2. For Hydrogen, gH=gpge=4g_H = g_p g_e = 4.

Since mHmpm_H \approx m_p and defining Q=(mp+memH)c213.6eVQ = (m_p + m_e - m_H)c^2 \approx 13.6 \, \mathrm{eV} as the ionization energy, we have

nHnenp=(mekT2π2)3/2exp(QkT)(5.5.15)\tag{5.5.15} \frac{n_H}{n_e n_p} = \left( \frac{m_e kT}{2 \pi \hbar^2} \right)^{-3/2} \, \exp(-\frac{Q}{kT})

This is called the Saha equation.

We define the fraction ionization XX as

X=npnp+nH=nenbary=nenbary(5.5.16)\tag{5.5.16} X = \frac{n_p}{n_p + n_H} = \frac{n_e}{n_\text{bary}} = \frac{n_e}{n_\text{bary}}

From this we get that nH=1XXnpn_H = \frac{1-X}{X} n_p. For overall charge neutrality, we have ne=npn_e = n_p. Therefore

1XX=np(mekT2π2)3/2exp(QkT)(5.5.17)\tag{5.5.17} \frac{1-X}{X} = n_p \left( \frac{m_e kT}{2 \pi \hbar^2} \right)^{-3/2} \, \exp(-\frac{Q}{kT})

The ratio of baryonic number density and photon number density is

η=nbarynγ=npXnγ\eta = \frac{n_\text{bary}}{n_\gamma} = \frac{n_p}{X n_\gamma}

For a blackbody spectrum, nγ=0.243(kTc)3n_\gamma = 0.243 \left( \frac{kT}{\hbar c} \right)^3.

For X=0.5X = 0.5 and η=5.5×1010\eta = 5.5 \times 10^{-10}, we get Trec=3740KT_\text{rec} = 3740 \, \mathrm{K} and zrec=1370z_\text{rec} = 1370, which is t=240,000yrst = 240,000 \, \mathrm{yrs} after big bang.

The rate of photon scattering is

Γ=cλ=neσec=X(1+z)3nbary,0σec(5.5.18)\tag{5.5.18} \Gamma = \frac{c}{\lambda} = n_e\, \sigma_e\, c= X \, (1+z)^3 \, n_{\text{bary},0} \, \sigma_e \, c

where σe\sigma_e is the Thomson cross-section. Since H2H02=Ωm,0(1+z)3\frac{H^2}{H_0^2} = \Omega_{\text{m},0} (1+z)^3, we get

1+zdec43Xzdec3/2(5.5.19)\tag{5.5.19} 1+z_\text{dec} \approx \frac{43}{X z_\text{dec}^{3/2}}

The redshift of photon decoupling was found to be zdec=1120z_\text{dec} = 1120. However, the exact value is somewhat smaller, about zdec=1090z_\text{dec} = 1090, for which Tdec=2970KT_\text{dec} = 2970 \, \mathrm{K} and t=371,000yrst = 371,000 \, \mathrm{yrs}. The last scattering happened too around this time, zlszdec1100z_\text{ls} \approx z_\text{dec} \approx 1100.