Extended source effects can be seen in gravitational lensing events when sources cross critical lines. Those events probe the stellar intensity profile and could be used to measure limb darkening coefficients to test stellar model predictions. A data base of accurately measured stellar profiles will be needed to correctly subtract the stellar flux in planetary transient events. The amount of data that is being and will be produced in current and future microlensing surveys, from the space and the ground, requires algorithms that can quickly compute light curves for different source-lens configurations. Based on the convolution method we describe a general formalism to compute those curves for single lenses. We develop approximations in terms of quadratures of elliptic integrals that we integrate by solving the associated first order differential equations. We construct analytic solutions for a limb darkening and, for the first time, for a parabolic profile that are accurate at the $\sim 1-3\%$ and $0.5\%$ level, respectively. These solutions can be computed orders of magnitude faster than other integration routines. They can be implemented in pipelines processing large data sets to extract stellar parameters in real time.

Time delays in gravitational lenses can be used to determine the Hubble constant and the lens potential. In future surveys, many gravitational lenses can be discovered, and their time delays and image positions can in principle be measured. Using an elliptical power-law potential, we show that combinations of image positions and time delays for quadruple lenses yield simple analytical expressions that are connected with observable quantities. These relations can be used to obtain the approximate axis ratio q, the Einstein radius and the slope. We apply this method to RX J1131-1231, and show that our analytical results match the full numerical determinations approximately. Our approach can quickly determine rough values of lens parameters, which can then be used as initial guesses for further refinement through numerical modelling and may be useful for automated lens search in large surveys.

We study three-dimensional microlensing where two lenses are located at different distances along the line of sight. We formulate the lens equation in complex notations and recover several previous results. There are in total either 4 or 6 images, with an equal number of images with positive and negative parities. We find that the sum of signed magnifications for six image configurations is unity. Furthermore, we show that the light curves can be qualitatively different from those for binary lensing in a single plane. In particular, the magnifications between a `U'-shaped caustic crossing can be close to unity, rather than having a minimum magnification of 3 as in the single plane binary lensing. There is only a small probability three-dimensional microlensing events will be seen in microlensing toward the Galactic centre. It is more likely they will be first seen in cosmological microlensing.

Most galaxies host central supermassive black holes. As two galaxies merge, the black holes also merge. The final single black hole may suffer a kick due to asymmetric gravitational radiation and may not be at the centre of the galaxy; off-centre black holes may also be produced by other means such as sustained acceleration due to asymmetric jet power. We model the main galaxy as a singular isothermal sphere and the black hole as an off-centre point lens, and study the critical curves and caustics using complex notation. We identify the critical parameters that govern the transitions in the topology of critical curves, caustics and pseudo-caustics, and find the number of images can be two, three, four and five. We show examples of image configurations, including cases where three highly de-magnified images are found close to the centre. The perturbation on the image magnification due to the black hole scales linearly with its mass in the off-centre case, and quadratically when the black hole is at the centre. Such images are difficult to observe unless high-contrast and high-resolution imaging facilities (e.g., the Square Kilometer Array in the radio) become available.

The observables in a strong gravitational lens are usually just the image positions and sometimes the flux ratios. We develop a new and simple algorithm which allows a set of models to be fitted exactly to the observations. Taking our cue from the strong body of evidence that early-type galaxies are close to isothermal, we assume that the lens is scale-free with a flat rotation curve. However, our algorithm allows full flexibility as regards the angular structure of the lensing potential. Importantly, all the free parameters enter linearly into the model and so the lens and flux ratio equations can always be solved by straightforward matrix inversion. We use this new algorithm to examine some of the claims made for substructure on the basis of "anomalous flux ratios". We demonstrate by explicit construction that some of the lenses for which substructure has been claimed can be well-fit by smooth lens models. This is especially the case when the systematic errors in the flux ratios (caused by microlensing or differential extinction or scatter-broadening) are taken into account. However, there is certainly one system (B1422+231) for which the existing smooth models are definitely inadequate and for which substructure may be implicated. Within a few tens of kpc of galaxy centres, dynamical friction is efficient at dissolving any substructure. So, little substructure lies projected within the Einstein radius of most strong lenses. The numbers of strong lenses for which substructure is currently being claimed may be so large that this contradicts rather than supports cold dark matter theories.

We study gravitational lensing by the family of scale-free galaxies with flat rotation curves. The models are defined by a shape function, which prescribes the radius of the isophote as a function of position angle from the major axis. The critical curves are analytic, while the caustic network is reducible to a simple quadrature. The cusps are always located at the turning points of the shape function. We show that the models with exactly elliptic isophotes never admit butterfly or swallowtail cusps and so have at most 4 (or 5) images. Higher order imaging is brought about by deviations of the isophotes from pure ellipses -- such as pointedness caused by embedded disks or boxiness caused by recent merging. The criteria for the onset of sextuple and octuple imaging can be calculated analytically in terms of the ellipticity and the fourth-order Fourier coefficients (a_4 and b_4) used by observers to parametrise the isophote shapes. The 6 or 8 images are arranged roughly in a circle, which appears as an incomplete Einstein ring if inadequately resolved. Using data on the shapes of elliptical galaxies and merger remnants, we estimate that 1% of all multiply imaged quasars may be sextuplet systems or higher. Forthcoming satellites like the Global Astrometric Interferometer for Astrophysics (GAIA) will provide datasets of roughly 4000 multiply imaged systems, and so about 40 will show sextuple imaging or higher.

Recent observations indicate that many if not all galaxies host massive central black holes. In this paper we explore the influence of black holes on the lensing properties. We model the lens as an isothermal ellipsoid with a finite core radius plus a central black hole. We show that the presence of the black hole substantially changes the critical curves and caustics. If the black hole mass is above a critical value, then it will completely suppress the central images for all source positions. Realistic central black holes likely have masses below this critical value. Even in such sub-critical cases, the black hole can suppress the central image when the source is inside a zone of influence, which depends on the core radius and black hole mass. In the sub-critical cases, an additional image may be created by the black hole in some regions, which for some radio lenses may be detectable with high-resolution and large dynamic-range VLBI maps. The presence of central black holes should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.

We study gravitational lens time delays for a general family of lensing potentials, which includes the popular singular isothermal elliptical potential and singular isothermal elliptical density distribution but allows general angular structure. Using a novel approach, we show that the time delay can be cast in a very simple form, depending only on the observed image positions. Including an external shear changes the time delay proportional to the shear strength, and varying the radial profile of the potential changes the time delay approximately linearly. These analytic results can be used to obtain simple estimates of the time delay and the Hubble constant in observed gravitational lenses. The naive estimates for four of five time delay lenses show surprising agreement with each other and with local measurements of H_0; the complicated Q 0957+561 system is the only outlier. The agreement suggests that it is reasonable to use simple isothermal lens models to infer H_0, although it is still important to check this conclusion by examining detailed models and by measuring more lensing time delays.

We study the problem of gravitational lensing by an isothermal elliptical density galaxy in the presence of a tidal perturbation. When the perturbation is fairly strong and oriented near the galaxy's minor axis, the lens can produce image configurations with six or even eight highly magnified images lying approximately on a circle. We classify the caustic structures in the model and identify the range of models that can produce such lenses. Sextuple and octuple lenses are likely to be rare because they require special lens configurations, but a full calculation of the likelihood will have to include both the existence of lenses with multiple lens galaxies and the strong magnification bias that affects sextuple and octuple lenses. At optical wavelengths these lenses would probably appear as partial or complete Einstein rings, but at radio wavelengths the individual images could probably be resolved.

We present a new method of studying quadruple lenses in elliptical power-law potentials parameterized by $\psi(x,y) \propto (x^2+y^2/q^2)^{\beta/2}/\beta (0 \leq \beta < 2)$. For this potential, the moments of the four image positions weighted by signed magnifications (magnification times parity) have very simple properties. In particular, we find that the zeroth moment -- the sum of four signed magnifications satisfies $\simeq 2/(2-\beta)$; the relation is exact for $\beta=0$ (point-lens) and $\beta=1$ (isothermal potential), independent of the axial ratio. Similar relations can be derived when a shear is present along the major or minor axes. These relations, however, do not hold well for the closely-related elliptical density distributions. For a singular isothermal elliptical density distribution without shear, the sum of signed magnifications for quadruple lenses is $\approx 2.8$, again nearly independent of the ellipticity. For the same distribution with shear, the total signed magnification is around 2-3 for most cases, but can be significantly different for some combinations of the axial ratio and shear where more than four images can appear.

Extended source size effects have been detected in photometric monitoring of gravitational microlensing events. We study similar effects in the centroid motion of an extended source lensed by a point mass. We show that the centroid motion of a source with uniform surface brightness can be obtained analytically. For a source with circularly symmetric limb-darkening profile, the centroid motion can be expressed as a one-dimensional integral, which can be evaluated numerically. We found that when the impact parameter is comparable to the source radius, the centroid motion is significantly modified by the finite source size. In particular, when the impact parameter is smaller than the source radius, the trajectories become clover-leaf like. Such astrometric motions can be detected using space interferometers such as the Space Interferometry Mission. Such measurements offer exciting possibilities of determining stellar parameters such as stellar radius to excellent accuracy.

We study the critical curves, caustics, and multiple imaging due to one of the simplest many-body gravitational lens configurations: equal-mass point masses on the vertices of a regular polygon. Some examples of the critical curves and caustics are also displayed. We pose the problem of determining the maximum number of lensed images due to regular-polygon and general point mass configurations. Our numerical simulations suggest a maximum that is linear, rather than quadratic, in the number of point masses.

We study a general elliptical potential of the form $\psi(x^2+y^2/q^2)~ (0<q\le 1)$ plus an additional shear (with an arbitrary direction) as models for the observed quadruple lenses. It is shown that a minimum additional shear is needed even just to reproduce the observed positions alone. We also obtain the dependence of the axial ratio, $q$, on the orientation of the major axis of potential. A general relation also exists between the shear, the position angle and axial ratio of the lensing galaxy. The relation shows a generic degeneracy in modelling quadruple lenses. In particular, it shows that only the ratio of the ellipticity, $\epsilon\equiv (1-q^2)/(1+q^2)$, to the magnitude of shear, $\gamma$ can be determined. All these results are valid regardless of the radial profile of the potential. Our formalism applies when the galaxy position is observed, which is the case for seven of the eight known quadruple lenses. Application to these seven cases reveals two quadruple lenses CLASS 1608+656 and HST 12531-2914, requiring highly significant shear with magnitude $\approx 0.2$. For HST 12531-2914, there must be a misalignment between the major axis of light and the major axis of potential (mass). We conclude that detailed modelling of quadruple lenses can yield valuable quantitative information about the shape of lensing galaxies and their dark matter halos.

We discuss the optical/UV/Xray spectra of AGN within the frame of the corona model of Witt, Czerny and Zycki (1996). In this model both the disk and the corona accrete and release energy through the viscosity. The relative strength the the disk and the corona is therefore determined by the model. Translated into the spectra, this model allows to predict the spectral index alpha_{ox}, measuring the relative strenth of the big blue bump with respect to the hard X-ray power law. Comparison of the predicted and observed distributions of alpha_{ox} allow us to conclude that the luminosity to the Eddington luminosity ratio in quasars cover typically the range 0.01 - 0.1 and it is broader in Seyferts, 0.001 - 0.3, or even higher. We identify Narrow Line Seyfert galaxies as objects having high luminosity to the Eddington luminosity ratio although, in principle, the second branch of high big blue bump objects for this ratio being below 0.001 is also predicted.

We show here that quadruple lenses can be useful laboratories to probe whether the potential of the lensing galaxy is purely elliptical or whether an additional distortion is present in the deflector plane. For this test we only have to know the relative image positions of the quadruple lens system and the (relative) center of light position of the lensing galaxy. Furthermore we introduce new methods which easily allow us to determine the location (rotation angle relative to the image positions) of the major axis of the lensing galaxy. In due course we can determine the parity of the four images as well. We apply these methods to the 8 currently known quadruple lenses and find that in the case of MG 0414+0534, CLASS 1608+656 and HST 12531-2914 it is impossible to accommodate the relative image positions and the galaxy position with any elliptical potential whereas the other five cases can be described very well with a simple elliptical potential. This method will have important impacts for $\chi^2$-fits and the reconstruction of galaxy models for quadruple lenses.

We study an accretion disc with a hot continuous corona. We assume that the corona itself accretes and therefore it is powered directly by the release of the gravitational energy and cooled by radiative interaction with the disc. We consider the vertical structure of such a corona and show that the radial infall is accompanied by strong vertical outflow. Such a model has two important consequences: (i) at a given radius the corona forms only for accretion rate larger than a limiting value and the fraction of energy dissipated in the corona decreases with increasing accretion rate, and (ii) the disc spectra are significantly softer in the optical/UV band in comparison with the predictions of standard accretion discs due to the mass loss and the decrease of internal dissipation in the disc. Both trends correspond well to the mean spectra of radio quiet AGN and observed luminosity states in galactic black hole candidates.