We show in detail that the recently derived expression for evaluating the integrated Sachs-Wolfe (ISW) temperature shift in the cosmic microwave background (CMB) caused by individual embedded (compensated) lenses is equivalent to the conventional approach for flat background cosmologies. The conventional approach requires evaluating an integral of the time derivative of the lensing potential, whereas the new Fermat potential approach is simpler and only requires taking a derivative of the potential part of the time delay.

Recently Graphics Processing Units (GPUs) have been used to speed up very CPU-intensive gravitational microlensing simulations. In this work, we use the Xeon Phi coprocessor to accelerate such simulations and compare its performance on a microlensing code with that of NVIDIA's GPUs. For the selected set of parameters evaluated in our experiment, we find that the speedup by Intel's Knights Corner coprocessor is comparable to that by NVIDIA's Fermi family of GPUs with compute capability 2.0, but less significant than GPUs with higher compute capabilities such as the Kepler. However, the very recently released second generation Xeon Phi, Knights Landing, is about 5.8 times faster than the Knights Corner, and about 2.9 times faster than the Kepler GPU used in our simulations. We conclude that the Xeon Phi is a very promising alternative to GPUs for modern high performance microlensing simulations.

Active galactic nuclei (AGNs) and quasars are important astrophysical objects to understand. Recently, microlensing observations have constrained the size of the quasar X-ray emission region to be of the order of 10 gravitational radii of the central supermassive black hole. For distances within a few gravitational radii, light paths are strongly bent by the strong gravity field of the central black hole. If the central black hole has nonzero angular momentum (spin), a photon's polarization plane will be rotated by the gravitational Faraday effect. The observed X-ray flux and polarization will then be influenced significantly by the strong gravity field near the source. Consequently, linear gravitational lensing theory is inadequate for such extreme circumstances. We present simple algorithms computing strong lensing effects of Kerr black holes, including effects on polarization. Our algorithms are realized in a program "KERTAP" in two versions: MATLAB and Python. The key ingredients of KERTAP are: a graphic user interface, a {\it backward} ray-tracing algorithm, a polarization propagator dealing with gravitational Faraday rotation, and algorithms computing observables such as flux magnification and polarization angles. Our algorithms can be easily realized in other programming languages such as FORTRAN, C, and C++. The MATLAB version of KERTAP is parallelized using the MATLAB Parallel Computing Toolbox and the Distributed Computing Server. The Python code was sped up using Cython and supports full implementation of MPI using 'mpi4py' package. As an example, we investigate the inclination angle dependence of the observed polarization and the strong lensing magnification of AGN X-ray emission. We conclude that it is possible to perform complex numerical-relativity-related computations using interpreted languages such as MATLAB and Python.

The photometry profile of the integrated Sachs-Wolfe (ISW) effect recently obtained by the Planck consortium by stacking patches of Cosmic Microwave Background (CMB) sky maps around a large number of cosmic voids, contains a cold ring at about half the void's effective radius surrounded by a hot ring near the void's boundary. The source of the temperature structure is assumed to be the ISW effect but the exact cause of the ringed structure is not currently well understood, particularly the outer hot ring. Numerical simulations have suggested that hot/cold ring structures can be produced by motions associated with nonlinear growths of cosmic structures whose gravitational potentials produce the ISW effect. We have recently developed the embedded lens theory and the Fermat potential formalism which can be used to model the ISW effect caused by intervening individual lens inhomogeneities evolving arbitrarily. This theory only requires knowledge of the void's projected mass profile as a function of the passing CMB photons' impact radius and the rate of change of that mass distribution at passage. We present two simple embedded void lens models with evolving mass densities and investigate the ISW effect caused by these lenses. Both models posses expanding mass shells which produce hot rings around central cold regions, consistent with the recent observations. By adding a small over-density at the void's center we can produce the slight positive temperature excess hinted at in Planck's photometric results. We conclude that the embedded lens theory and the Fermat potential formalism is well suited for modeling the ISW effect.

We describe a method using the integrated Sachs-Wolfe (ISW) effect caused by individual inhomogeneities to determine the cosmological parameters, $H_0$, $\Omega_{\rm m}$, and $\Omega_\Lambda$, etc. This ISW-redshift test requires detailed knowledge of the internal kinematics of a set of individual density perturbations, e.g., galaxy clusters and/or cosmic voids, in particular their density and velocity profiles, and their mass accretion rates. It assumes the density perturbations are isolated and imbedded (equivalently compensated) and makes use of the newly found relation between the ISW temperature perturbation of the CMB and the Fermat potential of the lens. Given measurements of the amplitudes of the temperature variations in the CMB caused by such clusters or voids at various redshifts and estimates of their angular sizes or masses, one can constrain the cosmological parameters. More realistically, the converse is more likely, i.e., if the background cosmology is sufficiently constrained, measurement of ISW profiles of clusters and voids (e.g., hot and cold spots and rings) can constrain dynamical properties of the dark matter, including accretion, associated with such lenses and thus constrain the evolution of these objects with redshift.

We derive the Fermat potential for a spherically symmetric lens embedded in an FLRW cosmology and use it to investigate the late-time integrated Sachs-Wolfe (ISW) effect, i.e., secondary temperature fluctuations in the cosmic microwave background (CMB) caused by individual large scale clusters and voids. We present a simple analytical expression for the temperature fluctuation in the CMB across such a lens as a derivative of the lens' Fermat potential. This formalism is applicable to both linear and nonlinear density evolution scenarios, to arbitrarily large density contrasts, and to all open and closed background cosmologies. It is much simpler to use and makes the same predictions as conventional approaches. In this approach the total temperature fluctuation can be split into a time-delay part and an evolutionary part. Both parts must be included for cosmic structures that evolve and both can be equally important. We present very simple ISW models for cosmic voids and galaxy clusters to illustrate the ease of use of our formalism. We use the Fermat potentials of simple cosmic void models to compare predicted ISW effects with those recently extracted from WMAP and \emph{Planck} data by stacking large cosmic voids using the aperture photometry method. If voids in the local universe with large density contrasts are no longer evolving we find that the time delay contribution alone predicts values consistent with the measurements. However, we find that for voids still evolving linearly, the evolutionary contribution cancels a significant part of the time delay contribution and results in predicted signals that are much smaller than recently observed.

We present a simple gravitational lens model to illustrate the ease of using the embedded lensing theory when studying cosmic voids. It confirms the previously used repulsive lensing models for deep voids. We start by estimating magnitude fluctuations and weak lensing shears of background sources lensed by large voids. We find that sources behind large ($\sim$$90\,\rm Mpc$) and deep voids (density contrast about $-0.9$) can be magnified or demagnified with magnitude fluctuations of up to $\sim$$0.05\,\rm mag$ and that the weak-lensing shear can be up to the $\sim$$10^{-2}$ level in the outer regions of large voids. Smaller or shallower voids produce proportionally smaller effects. We investigate the "wiggling" of the primary cosmic microwave background (CMB) temperature anisotropies caused by intervening cosmic voids. The void-wiggling of primary CMB temperature gradients is of the opposite sign to that caused by galaxy clusters. Only extremely large and deep voids can produce wiggling amplitudes similar to galaxy clusters, $\sim$$\rm 15\,\mu K$ by a large void of radius $\sim$$4^\circ$ and central density contrast $-0.9$ at redshift 0.5 assuming a CMB background gradient of $\sim$$\rm10\,\mu K\, arcmin^{-1}$. The dipole signal is spread over the entire void area, and not concentrated at the lens' center as it is for clusters. Finally we use our model to simulate CMB sky maps lensed by large cosmic voids. Our embedded theory can easily be applied to more complicated void models and used to study gravitational lensing of the CMB, to probe dark-matter profiles, to reduce the lensing-induced systematics in supernova Hubble diagrams, as well as study the integrated Sachs-Wolfe effect.

We present a simplified version of the lowest-order embedded point mass gravitational lens theory and then make the extension of this theory to any embedded transparent lens. Embedding a lens effectively reduces the gravitational potential's range, i.e., partially shields the lensing potential because the lens mass is made a contributor to the mean mass density of the universe and not simply superimposed upon it. We give the time-delay function for the embedded point mass lens from which we can derive the simplified lens equation by applying Fermat's least-time principle. Even though rigorous derivations are only made for the point mass in a flat background, the generalization of the lens equation to lowest-order for any distributed lens in any homogeneous background is obvious. We find from this simplified theory that embedding can introduce corrections above the few percent level in weak lensing shears caused by large clusters but only at large impacts. The potential part of the time delay is also affected in strong lensing at the few percent level. Additionally we again confirm that the presence of a cosmological constant alters the gravitational deflection of passing photons.

Recent quasar microlensing observations have constrained the sizes of X-ray emission regions to be within about 10 gravitational radii of the central supermassive black hole. Therefore, the X-ray emission from lensed quasars is first strongly lensed by the black hole before it is lensed by the foreground galaxy and star fields. We present a scheme that combines the initial strong lensing of a Kerr black hole with standard linearized microlensing by intervening stars. We find that X-ray microlensed light curves incorporating Kerr strong gravity can differ significantly from standard curves. The amplitude of the fluctuations in the light curves can increase or decrease by ~0.65-0.75 mag by including Kerr strong gravity. Larger inclination angles give larger amplitude fluctuations in the microlensing light curves. Consequently, current X-ray microlensing observations might have under or overestimated the sizes of the X-ray emission regions. We estimate this bias using a simple metric based on the amplitude of magnitude fluctuations. The half light radius of the X-ray emission region can be underestimated up to ~50% or overestimated up to ~20%. Underestimates are found in most situations we have investigated. The only exception is for a disk with large spin, radially flat emission profile, and observed nearly face on, where an overestimate is found. Thus, more accurate microlensing size constraints should be obtainable by including Kerr lensing. The caustic crossing time can differ by months after including Kerr strong gravity. A simultaneous monitoring of gravitational lensed quasars in both X-ray and optical bands with densely sampled X-ray light curves might reveal this feature. We conclude that it should be possible to constrain important parameters such as inclination angles and black hole spins from combined Kerr and microlensing effects.

We give analytic expressions for image properties of objects seen around point mass lenses embedded in a flat $\Lambda$CDM universe. An embedded lens in an otherwise homogeneous universe offers a more realistic representation of the lens's gravity field and its associated deflection properties than does the conventional linear superposition theory. Embedding reduces the range of the gravitational force acting on passing light beams thus altering all quantities such as deflection angles, amplifications, shears and Einstein ring sizes. Embedding also exhibits the explicit effect of the cosmological constant on these same lensing quantities. In this paper we present these new results and demonstrate how they can be used. The effects of embedding on image properties, although small i.e., usually less than a fraction of a percent, have a more pronounced effect on image distortions in weak lensing where the effects can be larger than 10%. Embedding also introduces a negative surface mass density for both weak and strong lensing, a quantity altogether absent in conventional Schwarzschild lensing. In strong lensing we find only one additional quantity, the potential part of the time delay, which differs from conventional lensing by as much as 4%, in agreement with our previous numerical estimates.

We give the lens equation for light deflections caused by point mass condensations in an otherwise spatially homogeneous and flat universe. We assume the signal from a distant source is deflected by a single condensation before it reaches the observer. We call this deflector an embedded lens because the deflecting mass is part of the mean density. The embedded lens equation differs from the conventional lens equation because the deflector mass is not simply an addition to the cosmic mean. We prescribe an iteration scheme to solve this new lens equation and use it to compare our results with standard linear lensing theory. We also compute analytic expressions for the lowest order corrections to image amplifications and distortions caused by incorporating the lensing mass into the mean.

We compute time delays for gravitational lensing in a flat LambdaCDM Swiss cheese universe. We assume a primary and secondary pair of light rays are deflected by a single point mass condensation described by a Kottler metric (Schwarzschild with Lambda) embedded in an otherwise homogeneous cosmology. We find that the cosmological constant's effect on the difference in arrival times is non-linear and at most around 0.002% for a large cluster lens; however, we find differences from time delays predicted by conventional linear lensing theory that can reach ~4% for these large lenses. The differences in predicted delay times are due to the failure of conventional lensing to incorporate the lensing mass into the mean mass density of the universe.

We compute the deflection angle to order (m/r_0)^2 and m/r_0*Lambda r_0^2 for a light ray traveling in a flat LambdaCDM cosmology which encounters a completely condensed mass region. We use a Swiss cheese model for the inhomogeneities and find that the most significant correction to the Einstein angle occurs not because of the non-linear terms but instead occurs because the condensed mass is embedded in a background cosmology.The Swiss cheese model predicts a decrease in the deflection angle of ~2% for weakly lensed galaxies behind the rich cluster A1689, and that the reduction can be as large as ~5% for similar rich clusters at z\approx 1. Weak lensing deflection angles caused by galaxies can likewise be reduced by as much as ~4%. We show that the lowest order correction in which Lambda appears is proportional to m/r_0*\sqrt{Lambda r_0^2} and could cause as much as a ~0.02% increase in the deflection angle for light that passes through a rich cluster. The lowest order non-linear correction in the mass is proportional to m/r_0*\sqrt{m/r_0} and can increase the deflection angle by ~0.005% for weak lensing by galaxies.

We show that it is possible to equate the intensity reduction of a light wave caused by weak absorption with a geometrical reduction in intensity caused by a "transverse" conformal transformation of the spacetime metric in which the wave travels. We are consequently able to modify Gordon's optical metric to account for electromagnetic properties of ponderable material whose properties include both refraction and absorption. Unlike refraction alone however, including absorption requires a modification of the optical metric that depends on the eikonal of the wave itself. We derive the distance-redshift relation from the modified optical metric for Friedman-Lema\^itre-Robertson-Walker spacetimes whose cosmic fluid has associated refraction and absorption coefficients. We then fit the current supernovae data and provide an alternate explanation (other than dark energy) of the apparent acceleration of the universe.

We show that Gordon's optical metric on a curved spacetime can be generalized to include absorption by allowing the metric to become complex. We demonstrate its use in the realm of geometrical optics by giving three simple examples. We use one of these examples to compute corrected distance-redshift relations for Friedman-Lema\^itre-Robertson-Walker models in which the cosmic fluid has an associated complex index of refraction that represents grey extinction. We then fit this corrected Hubble curve to the type Ia supernovae data and provide a possible explanation (other than dark energy) of the deviation of these observations from dark matter predictions.

We review Gordon's optical metric and the transport equations for the amplitude and polarization of a geometrical optics wave traveling in a gravity field. We apply the theory to the FLRW cosmologies by associating a refraction index with the cosmic fluid. We then derive an expression for the accumulated effect of a refraction index on the distance redshift relations and fit the Hubble curve of current supernova observations with a non-accelerating cosmological model. We also show that some observational effects caused by inhomogeneities, e.g. the Sachs-Wolfe effect, can be interpreted as being caused by an effective index of refraction, and hence this theory could extend to other speed of light communications such as gravitational radiation and neutrino fluxes.

We derive the radiative transfer equation for arbitrary stationary relativistic flows in stationary spacetimes, i.e. for steady-state transfer problems. We show how the standard characteristics method of solution developed by Mihalas and used throughout the radiative transfer community can be adapted to multi-dimensional applications with isotropic sources. Because the characteristics always coincide with geodesics and can always be specified by constants, direct integration of the characteristics derived from the transfer equation as commonly done in 1-D applications is not required. The characteristics are known for a specified metric from the geodesics. We give details in both flat and static spherically symmetric spacetimes. This work has direct application in 3-dimensional simulations of supernovae, gamma-ray bursts, and active galactic nuclei, as well as in modeling neutron star atmospheres.

A description of the light curves of microlensed Type Ia supernovae (SNe Ia) as extended and expanding sources is presented. We give examples of what microlensing by stellar-mass Chang-Refsdal lenses can do to a small percentage of supernova light curves. We find that in addition to overall brightening, significant changes in light--curve shapes can also occur. Peaks can be distorted, plateaus can appear, and even secondary peaks can be formed. The effects of both the relative motion of the lens and the supernova and the expansion of the supernova are given and compared. The effects of relative motion are more pronounced when a distant supernova ($z_s \sim 1$) impacts well within the Einstein ring of a nearby microlens ($z_d \sim 0.05$) and are less important for more distant deflectors. We also find that the increase in shear that comes with increased deflector distance tends to reduce the time variability of microlensing. We briefly discuss the probability of observing these effects.

A brief description of the deformed spectra of microlensed SNe Ia is presented. We show that microlensing amplification can have significant effects on line profiles. The resonance-scattering code SYNOW is used to compute the intensity profile in the rest frame of the supernova. The observed (microlensed) spectral lines are predicted assuming a simple stellar-size deflector, and are compared to unlensed cases to show the effects microlensing by solar-size deflectors can have on spectral lines. We limit our work to spherically symmetric deflectors.

A detailed description of the apparent light curves of microlensed SNe Ia as extended and expanding sources is presented. We show that microlensing amplification can have significant effects on supernova observation. A model light curve is used to compare lensed and unlensed cases and we find that significant changes in shape can occur because of microlensing. We briefly discuss the probability of observing such effects as well. Throughout the paper we consider the redshift of the source to be 1.0 and we limit our work to spherically symmetric deflectors.

The differential equation governing distance-redshift for partially filled-beam optics in pressure-freeFriedmann-Lema\^\i tre-Robertson-Walker Cosmology (FLRW) is shown to be the Lame$^{\prime}$ equation. The distance-redshift, $D(z)$, discussed is appropriate for observations in inhomogeneous cosmologies for which lensing by masses external to the observing beam are negligible and for which lensing by transparent matter within the beam can be approximated by a homogeneous mass density expanding with the FLRW background. Some solutions of the derived Lame$^{\prime}$ equation are given in terms of Weierstrass elliptic integrals. A new simplified and useful expression for filled-beam $D(z)$ in standard flat FLRW is also given.

If quantum fields exist in extra compact dimensions, they will give rise to a quantum vacuum or Casimir energy. That vacuum energy will manifest itself as a cosmological constant. The fact that supernova and cosmic microwave background data indicate that the cosmological constant is of the same order as the critical mass density to close the universe supplies a lower bound on the size of the extra dimensions. Recent laboratory constraints on deviations from Newton's law place an upper limit. The allowed region is so small as to suggest that either extra compact dimensions do not exist, or their properties are about to be tightly constrained by experimental data.

Distance--redshift relations are given in terms of associated Legendre functions for partially filled beam observations inspatially flat Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmologies. These models are dynamically pressure-free, flat FLRW on large scales but, due to mass inhomogeneities, differ in their optical properties. The partially filled beam area-redshift equation is a Lame$^{\prime}$ equation for arbitrary FLRW and is shown to simplify to the associated Legendre equation for the spatially flat, i.e. $\Omega_0=1$ case. We fit these new analytic Hubble curves to recent supernovae (SNe) data in an attempt to determine both the mass parameter $\Omega_m$ and the beam filling parameter $\nu$. We find that current data are inadequate to limit $\nu$. However, we are able to estimate what limits are possible when the number of observed SNe is increased by factor of 10 or 100, sample sizes achievable in the near future with the proposed SuperNova Acceleration Probe satellite.

We present compact, analytic expressions for the age-redshift relation $\tau(z)$ for standard Friedmann-Lema\^ \itre-Robertson-Walker (FLRW) cosmology. The new expressions are given in terms of incomplete Legendre elliptic integrals and evaluate much faster than by direct numerical integration.

We give distance--redshift relations in terms of elliptic integrals for three different mass distributions of the Friedmann-Lema\^\i tre-Robertson-Walker (FLRW) cosmology. These models are dynamically pressure free FLRW on large scales but, due to mass inhomogeneities, differ in their optical properties. They are the filled-beam model (standard FLRW), the empty-beam model (no mass density exists in the observing beams) and the 2/3 filled-beam model. For special $\OM$-- $\OL$ values the elliptic integrals reduce to more familiar functions. These new expressions for distance-redshift significantly reduce computer evaluation times.

Analytic expressions for distance-redshift relations which have been corrected for the effects of inhomogeneities in the Friedmann-Lema\^itre-Robertson-Walker (FLRW) mass density are given in terms of Heun functions and are used to illustrate the significance of inhomogeneities on a determination of the mass parameter $\Omega_m$ and the cosmological constant $\Lambda$. The values of these parameters inferred from a given set of observations depend on the fractional amount of matter in inhomogeneities and can significantly differ from those obtained by using the standard magnitude-redshift ($m$-$z$) result for pure dust FLRW models. As an example a determination of $\Omega_m$ made by applying the homogeneous distance-redshift relation to SN 1997ap at $z=0.83$ could be as much as 50% lower than its true value.

Modifications in Friedmann-Lemaitre-Robertson-Walker (FLRW) Hubble diagrams caused by mass density inhomogeneities are used to illustrate possible effects on a determination of the mass parameter $\Omega_m$ and the cosmological constant $\Lambda$. The values of these parameters inferred from a given set of observations depend on the fractional amount of matter in inhomogeneities and can differ from those obtained by using standard FLRW Hubble diagrams by as much as a factor of two.

Analytic expressions for distance-redshift relations which have been corrected for the effects of some inhomogeneities in the Friedmann-Robertson-Walker (FRW) mass density are used to illustrate the significance of inhomogeneities on a determination of $q_0$ made by using Type Ia supernovae. The value of $q_0$ inferred from a given set of observations depends on the fractional amount of matter in inhomogeneities and is up to 50 percent larger than that obtained by using the standard Mattig $m$-$z$ result for pure dust FRW models.