Degeneracy of semi-analytic solution for axially symmetric lens equations

We derive a simple semi-analytical approximation for lens equations with an arbitrary radially symmetric mass density ρ( r ), when   r /ξ₀≪ 1  and ξ₀ is the scalelength of the density profile. At the strong lensing regime, which is mostly constrained by the inner part of the mass density profile, we assume ρ∝ r ^α.
A dark matter (DM) haloes (GNFW model) are parametrized through a shape parameter α, a concentration parameter c ₁ and the total mass M . We apply our semi-analytical model to show how the solutions of the axially symmetric lens equations are degenerated in respect to the parameters α and c ₁.
In the case of an asymmetric dual image lens system, similar effective degeneracy is produced when the geometry of the lens is relaxed. Because it is impossible to determine the exact location of the source image, a family of solutions is acquired when the mass of the lens object and location of the observed images are fixed.
Our results indicate that the amount of degeneration is only weakly affected by the asymmetry in the lensing geometry set-up, e.g. the observational effective degeneracy is very close to the true physical degeneracy of the Einstein ring solutions. Basically with high-enough values for the concentration parameter, the degeneracy spawns the whole range for the shape parameter  α=[−2.0, −1.0]  .     

Non-parametric strong lens inversion of Cl 0024+1654: illustrating the monopole degeneracy

The cluster lens Cl 0024+1654 is undoubtedly one of the most beautiful examples of strong gravitational lensing, providing five large images of a single source with well-resolved substructure. Using the information contained in the positions and the shapes of the images, combined with the null space information, a non-parametric technique is used to infer the strong lensing mass map of the central region of this cluster. This yields a strong lensing mass of  1.60 × 10¹⁴ M_⊙  within a 0.5  arcmin radius around the cluster centre. This mass distribution is then used as a case study of the monopole degeneracy, which may be one of the most important degeneracies in gravitational lensing studies and which is extremely hard to break. We illustrate the monopole degeneracy by adding circularly symmetric density distributions with zero total mass to the original mass map of Cl 0024+1654. These redistribute mass in certain areas of the mass map without affecting the observed images in any way. We show that the monopole degeneracy and the mass-sheet degeneracy together lie at the heart of the discrepancies between different gravitational lens reconstructions that can be found in the literature for a given object, and that many images/sources, with an overall high image density in the lens plane, are required to construct an accurate, high-resolution mass map based on strong lensing data.     

Gravitational lens magnification by Abell 1689: distortion of the background galaxy luminosity function

Gravitational lensing magnifies the observed flux of galaxies behind the lens. We use this effect to constrain the total mass in the cluster Abell 1689 by comparing the lensed luminosities of background galaxies with the luminosity function of an undistorted field. Under the assumption that these galaxies are a random sample of luminosity space, this method is not limited by clustering noise. We use photometric redshift information to estimate galaxy distance and intrinsic luminosity. Knowing the redshift distribution of the background population allows us to lift the mass/background degeneracy common to lensing analysis. In this paper we use nine filters observed over 12 h with the Calar Alto 3.5-m telescope to determine the redshifts of 1000 galaxies in the field of Abell 1689. Using a complete sample of 146 background galaxies we measure the cluster mass profile. We find that the total projected mass interior to 0.25  h ⁻¹ Mpc is M _2D(<0.25  h ⁻¹ Mpc)=(0.48±0.16)×10¹⁵  h ⁻¹ M_⊙, where our error budget includes uncertainties from the photometric redshift determination, the uncertainty in the offset calibration and finite sampling. This result is in good agreement with that found by number-count and shear-based methods and provides a new and independent method to determine cluster masses.     

Degeneracies and scaling relations in general power-law models for gravitational lenses

The time-delay in gravitational lenses can be used to derive the Hubble constant in a relatively simple way. The results of this method are less dependent on astrophysical assumptions than in many other methods. For systems with accurately measured positions and time-delays, the most important uncertainty is related to the mass model used. Simple parametric models like isothermal ellipsoidal mass distributions seem to provide consistent results with a reasonably small scatter when applied to several lens systems. We discuss a family of models with a separable radial power law and an arbitrary angular dependence for the potential   ψ = r ^β F ( θ )  . Isothermal potentials are a special case of these models with   β =1  . An additional external shear is used to take into account perturbations from other galaxies. Using a simple linear formalism for quadruple lenses, we can derive H ₀ as a function of the observables and the shear. If the latter is fixed, the result depends on the assumed power-law exponent according to   H ₀∝(2- β )/ β   . The effect of external shear is quantified by introducing a 'critical shear' γ _c as a measure for the amount of shear that changes the result significantly. The analysis shows that in the general case H ₀ and γ _c do not depend on the position of the lens galaxy. Spherical lens models with images close to the Einstein radius with fitted external shear differ by a factor of   β /2  from shearless models, leading to   H ₀∝2- β   in this case. We discuss these results and compare them with numerical models for a number of real lens systems.     

The growth of linear perturbations in generic defect models for structure formation

We study the growth of linear perturbations induced by a generic causal scaling source as a function of the cosmological parameters h ,     and     . We show that for wavenumbers k ≳0.01  h  Mpc⁻¹ the spectrum of density and velocity perturbations scales in a similar way to that found in inflationary models with primordial perturbations. We show that this result is independent of the more-or-less incoherent nature of the source, the small-scale power spectrum of the source and of deviations from scaling that naturally occur at late times if     .     

Evidence for accretion disc reprocessing in QSO 0957+561

We use archival g - and r -band photometry of the gravitational lens system QSO 0957+561A, B to estimate the intrinsic variability of the quasar during 1996 February–June. The light curves span 234 d with temporal resolutions of about 2.5 d. Both light curves display a single large-amplitude event, of ∼0.1 mag (max-to-min) in about 100 d, followed by small-amplitude variations of ∼0.02 mag on time-scales of tens of days. We find the r -band variations lag those at g by 3.4_−1.4^+1.5 d for the large-amplitude event. This lag is greater than zero at no less than 98 per cent confidence. The delayed coupling of the rest-frame UV intrinsic variations strongly suggests the existence of a stratified reprocessing region extending ∼light-days from the putative central black hole source. The observed lag is consistent with that expected from a reverberation within an irradiated accretion disc structure. However, any definitive statement requires further detailed theoretical modelling and high-quality, signal-to-noise ratio of about 100, optical/IR simultaneous monitoring with about 3-d resolution for approximately 6 months.     

Principal-component analysis of the cosmic microwave background anisotropies: revealing the tensor degeneracy

A principal-component analysis of cosmic microwave background (CMB) anisotropy measurements is used to investigate degeneracies among cosmological parameters. The results show that a degeneracy with tensor modes – the 'tensor degeneracy'– dominates uncertainties in estimates of the baryon and cold dark matter densities,   ω _b=Ω_b  h ²  ,   ω _c=Ω_c  h ²  , ¹ from an analysis of CMB anisotropies alone. The principal-component analysis agrees well with a maximum-likelihood analysis of the observations, identifying the main degeneracy directions and providing an impression of the effective dimensionality of the parameter space.     

Combining the baryon budget with cosmic microwave background radiation measurements

Measurements of the cosmic microwave background radiation (CMBR) provide a powerful tool with which to measure the primary cosmological parameters. However, there is a large degree of parameter degeneracy in simultaneous measurements of the matter density, Ω_m, and the Hubble parameter, H ₀. In the present paper we use the currently available CMBR data together with measurements of the cosmological baryon-to-photon ratio, η , from big bang nucleosynthesis, and the relative mass fraction of baryons in clusters to break the parameter degeneracy in measuring Ω_m and H ₀. We find that present data are inconsistent with the standard Ω=1, matter-dominated model. Our analysis favours a medium-density universe with a rather low Hubble parameter. This is compatible with new measurements of Type Ia supernovae, and the joint estimate of the two parameters is     and     . We stress that the upper bound on the Hubble parameter is likely to be much more uncertain than indicated here, because of the limited number of free parameters in our analysis.     

Limits on isocurvature perturbations from non-Gaussianity in WMAP temperature anisotropy

We study the effect of primordial isocurvature perturbations on non-Gaussian properties of cosmic microwave background (CMB) temperature anisotropies. We consider generic forms of the non-linearity of isocurvature perturbations which can be applied to a wide range of theoretical models. We derive analytical expressions for the bispectrum and the Minkowski Functionals for CMB temperature fluctuations to describe the non-Gaussianity from isocurvature perturbations. We find that the isocurvature non-Gaussianity in the quadratic isocurvature model, where the isocurvature perturbation S is written as a quadratic function of the Gaussian variable  σ,  S =σ²−〈σ²〉  , can give the same signal-to-noise ratio as   f _NL= 30  even if we impose the current observational limit on the fraction of isocurvature perturbations contained in the primordial power spectrum α. We give constraints on isocurvature non-Gaussianity from Minkowski Functionals using the Wilkinson Microwave Anisotropy Probe ( WMAP ) 5-year data. We do not find a significant signal of isocurvature non-Gaussianity. For the quadratic isocurvature model, we obtain a stringent upper limit on the isocurvature fraction  α < 0.070  (95 per cent CL) for a scale-invariant spectrum which is comparable to the limit obtained from the power spectrum.     

Perturbative signature of substructures in strong gravitational lenses

In the perturbative approach, substructures in the lens can be reduced to their effect on the two perturbative fields f ₁ and  d f ₀/dθ  . A simple generic model of an elliptical lens with a substructure situated near the critical radius is investigated in detail. Analytical expressions are derived for each perturbative field, and basic properties are analysed. The power spectrum of the fields is well approximated by a power law, resulting in significant tails at high frequencies. Another feature of the perturbation by a substructure is that the ratio of the power spectrum at order n of the two fields   R _n   is almost one. The ratio   R _n ≃ 1  is specific to substructures, for instance a higher order distortion  ( n > 2)  but with autosimilar isophotes will result in   R _n ∝ 1/ n ²  . Finally, the problem of reconstructing the perturbative field is investigated. Local field models are implemented and fitted to maximize image similarity in the source plane. The non-linear optimization is greatly facilitated, since in the perturbative approach the circular source solution is always known. Examples of image distortions in the subcritical regime due to substructures are presented, and analysed for different source shapes. Provided enough images and signal are available, the substructure field can be identified confidently. These results suggest that the perturbative method is an efficient tool to estimate the contribution of substructures to the mass distribution of lenses.     

Pixel lensing as a way to detect extrasolar planets in M31

We study the possibility to detect extrasolar planets in M31 through pixel-lensing observations. Using a Monte Carlo approach, we select the physical parameters of the binary lens system, a star hosting a planet, and we calculate the pixel-lensing light curve taking into account the finite source effects. Indeed, their inclusion is crucial since the sources in M31 microlensing events are mainly giant stars. Light curves with detectable planetary features are selected by looking for significant deviations from the corresponding Paczyński shapes. We find that the time-scale of planetary deviations in light curves increase (up to 3–4 d) as the source size increases. This means that only few exposures per day, depending also on the required accuracy, may be sufficient to reveal in the light curve a planetary companion. Although the mean planet mass for the selected events is about     , even small mass planets  ( M _P < 20 M_⊕)  can cause significant deviations, at least in the observations with large telescopes. However, even in the former case, the probability to find detectable planetary features in pixel-lensing light curves is at most a few per cent of the detectable events, and therefore many events have to be collected in order to detect an extrasolar planet in M31. Our analysis also supports the claim that the anomaly found in the candidate event PA-99-N2 towards M31 can be explained by a companion object orbiting the lens star.     

The new gravitational lens system B1030 + 074

We report the discovery of a new double-image gravitational lens system, B1030 + 074, which was found during the Jodrell Bank–VLA Astrometric Survey (JVAS). We have collected extensive radio data on the system using the VLA, MERLIN, the EVN and the VLBA, and optical observations using WFPC2 on the HST . The lensed images are separated by 1.56 arcsec and their flux density ratio at centimetric wavelengths is approximately 14:1, although the ratio is slightly frequency-dependent and the images appear to be time-variable. The HST pictures show both the lensed images and the lensing galaxy close to the weaker image. The lensing galaxy has substructure which could be a spiral arm or an interacting galaxy. We have modelled B1030 + 074 using a singular isothermal ellipsoid which yields a time delay of 156/ h ₅₀ d. This lens is likely to be suitable for the measurement of the Hubble constant.     

The lens and source of the optical Einstein ring gravitational lens ER 0047–2808

We perform a detailed analysis of the optical gravitational lens ER 0047–2808 imaged with the Wide Field Planetary Camera 2 on the Hubble Space Telescope . Using software specifically designed for the analysis of resolved gravitational lens systems, we focus on how the image alone can constrain the mass distribution in the lens galaxy. We find that the data are of sufficient quality to strongly constrain the lens model with no a priori assumptions about the source. Using a variety of mass models, we find statistically acceptable results for elliptical isothermal-like models with an Einstein radius of 1.17 arcsec. An elliptical power-law model  (Σ∝ R ^−β)  for the surface mass density favours a slope slightly steeper than isothermal with  β= 1.08 ± 0.03  . Other models including a constant mass-to-light ratio (M/L), pure Navarro, Frenk & White halo and (surprisingly) an isothermal sphere with external shear are ruled out by the data. We find the galaxy light profile can only be fit with a Sérsic plus point-source model. The resulting total  M/L _B   contained within the images is  4.7  h ₆₅± 0.3  . In addition, we find the luminous matter is aligned with the total mass distribution within a few degrees. This is the first time a resolved optical gravitational lens image has been quantitatively reproduced using a non-parametric source.
The source, reconstructed by the software, is revealed to have two bright regions, with an unresolved component inside the caustic and a resolved component straddling a fold caustic. The angular size of the entire source is ∼0.1 arcsec and its (unlensed) Lyα flux is  3 × 10⁻¹⁷ erg s⁻¹ cm⁻²  .     

Debris disc stirring by secular perturbations from giant planets

Detectable debris discs are thought to require dynamical excitation ('stirring'), so that planetesimal collisions release large quantities of dust. We investigate the effects of the secular perturbations of a planet, which may lie at a significant distance from the planetesimal disc, to see if these perturbations can stir the disc, and if so over what time-scale. The secular perturbations cause orbits at different semimajor axes to precess at different rates, and after some time   t _cross  initially non-intersecting orbits begin to cross. We show that   t _cross∝ a ^9/2_disc/( m _pl e _pl a ³_pl)  , where   m _pl, e _pl  and   a _pl  are the mass, eccentricity and semimajor axis of the planet, and   a _disc  is the semimajor axis of the disc. This time-scale can be faster than that for the growth of planetesimals to Pluto's size within the outer disc. We also calculate the magnitude of the relative velocities induced among planetesimals and infer that a planet's perturbations can typically cause destructive collisions out to 100 s of au. Recently formed planets can thus have a significant impact on planet formation in the outer disc which may be curtailed by the formation of giant planets much closer to the star. The presence of an observed debris disc does not require the presence of Pluto-sized objects within it, since it can also have been stirred by a planet not in the disc. For the star ε Eridani, we find that the known radial velocity planet can excite the planetesimal belt at 60 au sufficiently to cause destructive collisions of bodies up to 100 km in size, on a time-scale of 40 Myr.     

W Crv: the shortest-period Algol with non-degenerate components?

Radial velocity data for both components of W Crv are presented. In spite of providing full radial-velocity information, the new data are not sufficient to establish the configuration of this important system because of large seasonal light curve perturbations, which prevent a combined light curve/radial-velocity solution. It is noted that the primary minimum is free of photometric perturbations, and this property may help to explain the elusive source of these perturbations. Photometrically, the system appears to be a contact binary with poor or absent energy exchange, but such an explanation – in view of the presence of light curve perturbations – is no more plausible than any one of the semi-detached configurations, with either the more-massive or less-massive components filling the associated Roche lobes. Lengthening of the orbital period and the size of the less-massive component above its main-sequence value, both suggest that the system is the shortest-period (0.388 d) known Algol with non-degenerate components.     

Gravitational lensing of anisotropic sources

In strong gravitational lensing, the multiple images we see correspond to light rays that leave the source in slightly different directions. If the source emission is anisotropic, the images may differ from conventional lensing predictions (which assume isotropy). To identify scales on which source anisotropy may be important, we study the angle δ between the light rays emerging from the source, for different lensing configurations. If the lens has a power-law profile   M ∝ R ^γ  , the angle δ initially increases with lens redshift and then either diverges (for a steep profile  γ < 1  ), remains constant (for an isothermal profile  γ= 1  ), or vanishes (for a shallow profile  γ > 1  ) as   z _l→ z _s  . The scaling with lens mass is roughly  δ∝ M ^1/(2−γ)  . The results for an Navarro–Frenk–White (NFW) profile are qualitatively similar to those for a shallow power law, with δ peaking at about half the redshift of the source (not half the distance). In practice, beaming could modify the statistics of beamed sources lensed by massive clusters: for an opening angle  θ_jet  , there is a probability as high as   P ∼ 0.02–0.07(θ_jet/0.5°)⁻¹  that one of the lensed images may be missed (for  2 ≲ z _s≲ 6  ). Differential absorption within active galactic nuclei (AGNs) could modify the flux ratios of AGNs lensed by clusters; a sample of AGNs lensed by clusters could provide further constraints on the sizes of absorbing regions. Source anisotropy is not likely to be a significant effect in galaxy-scale strong lensing.     

Distribution of caustic-crossing intervals for Galactic binary-lens microlensing events

Detection of caustic crossings of binary-lens gravitational microlensing events is important because by detecting them one can obtain useful information about both the lens and the source star. In this paper, we compute the distribution of the intervals between two successive caustic crossings, f ( t _cc), for Galactic bulge binary-lens events to investigate the observational strategy for the optimal detection and resolution of caustic crossings. From this computation, we find that the distribution is highly skewed towards short t _cc and peaks at t _cc∼1.5 d. For the maximal detection of caustic crossings, therefore, prompt initiation of follow-up observations for intensive monitoring of events will be important. We estimate that, under the strategy of the current follow-up observations with a second caustic-crossing preparation time of ∼2 d, the fraction of events with resolvable caustic crossing is ∼80 per cent. We find that if the follow-up observations can be initiated within 1 d after the first caustic crossing by adopting more aggressive observational strategies, the detection rate can be improved to ∼90 per cent.     

Structures in a class of magnetized scale-free discs

We construct analytically stationary global configurations for both aligned and logarithmic spiral coplanar magnetohydrodynamics (MHD) perturbations in an axisymmetric background MHD disc with a power-law surface mass density  Σ₀∝ r ^−α  , a coplanar azimuthal magnetic field   B ₀∝ r ^−γ  , a consistent self-gravity and a power-law rotation curve   v ₀∝ r ^−β  , where v ₀ is the linear azimuthal gas rotation speed. The barotropic equation of state  Π∝Σ ⁿ   is adopted for both MHD background equilibrium and coplanar MHD perturbations where Π is the vertically integrated pressure and n is the barotropic index. For a scale-free background MHD equilibrium, a relation exists among  α, β, γ  and n such that only one parameter (e.g. β) is independent. For a linear axisymmetric stability analysis, we provide global criteria in various parameter regimes. For non-axisymmetric aligned and logarithmic spiral cases, two branches of perturbation modes (i.e. fast and slow MHD density waves) can be derived once β is specified. To complement the magnetized singular isothermal disc analysis of Lou, we extend the analysis to a wider range of  −1/4 < β < 1/2  . As an illustrative example, we discuss specifically the  β= 1/4  case when the background magnetic field is force-free. Angular momentum conservation for coplanar MHD perturbations and other relevant aspects of our approach are discussed.     

Determination of cosmology and evolution from K-band magnitude–redshift and number count observations

We determine cosmological and evolutionary parameters from the 3CR K -band Hubble diagram and K -band number counts, assuming that the galaxies in question undergo pure luminosity evolution. Separately the two data sets are highly degenerate with respect to choice of cosmological and evolutionary parameters, but in combination the degeneracy is resolved. Of models that either are flat or have  Ω_Λ=0  , the preferred ones are close to the canonical case  Ω_{cold  matter}=1  ,  Ω_Λ=0  , with luminosity evolution amounting to 1 mag brighter at   z =1  .     

How selection and weighting of astrometric observations influence the impact probability. The case of asteroid (99942) Apophis

The aim of this paper is to show that in the case of a low probability of asteroid collision with the Earth, the appropriate selection and weighting of the data are crucial for the impact investigation and for analysing the impact possibilities using extensive numerical simulations. By means of the Monte Carlo special method, a large number of 'clone' orbits have been generated. A full range of orbital elements in the six-dimensional parameter space, that is, in the entire confidence region allowed by the observational material, has been examined. On the basis of 1000 astrometric observations of (99942) Apophis, the best solutions for the geocentric encounter distance of  6.065 ± 0.081 R_⊕  (without perturbations by asteroids) or  6.064 ± 0.095 R_⊕  (including perturbations by the four largest asteroids) were derived for the close encounter with the Earth on 2029 April 13. The present uncertainties allow for special configurations ('keyholes') during this encounter that may lead to very close encounters in future approaches of Apophis. Two groups of keyholes are connected with the close encounter with the Earth in 2036 (within the minimal distance of  5.7736−5.7763 R_⊕  on 2029 April 13) and 2037 (within the minimal distance of  6.3359–6.3488 R_⊕  ). The nominal orbits for our most accurate models run almost exactly in the middle of these two impact keyhole groups. A very small keyhole for the impact in 2076 has been found between these groups at the minimal distance of 5.97347   R_⊕  . This keyhole is close to the nominal orbit. The present observations are not sufficiently accurate to eliminate definitely the possibility of impact with the Earth in 2036 and for many years after.