Gravitational microlensing: A parallel,large-data implementation

Gravitational lensing allows us to probe the structure of matter on a broad range of astronomical scales, and as light from a distant source traverses an intervening galaxy, compact matter such as planets, stars, and black holes act as individual lenses. The magnification from such microlensing results in rapid brightness fluctuations which reveal not only the properties of the lensing masses, but also the surface brightness distribution in the source. However, while the combination of deflections due to individual stars is linear, the resulting magnifications are highly non-linear, leading to significant computational challenges which currently limit the range of problems which can be tackled. This paper presents a new and novel implementation of a numerical approach to gravitational microlensing, increasing the scale of the problems that can be tackled by more than two orders of magnitude, opening up a new regime of astrophysically interesting problems.     

The Δθ–zs relation for gravitational lenses as a cosmological test

Recently, Park &38; Gott claimed that there is a statistically significant, strong, negative correlation between the image separation Δθ and source redshift z _s for gravitational lenses. This is somewhat puzzling if one believes in a flat ( k  = 0) universe, since in this case the typical image separation is expected to be independent of the source redshift, while one expects a negative correlation in a k  = −1 universe and a positive one in a k  = +1 universe. Park &38; Gott explored several effects that could cause the observed correlation, but no combination of these can explain the observations with a realistic scenario. Here, I explore this test further in three ways. First, I show that in an inhomogeneous universe a negative correlation is expected regardless of the value of k . Secondly, I test whether the Δθ– z _s relation can be used as a test to determine λ₀ and Ω₀, rather than just the sign of k . Thirdly, I compare the results of the test from the Park &38; Gott sample with those using other samples of gravitational lenses, which can illuminate (unknown) selection effects and probe the usefulness of the Δθ– z _s relation as a cosmological test.     

基于天文开源软件的卫星干涉测量数据处理系统

VLBI (Very Long Baseline Interferometry)技术观测卫星需要对干涉测量数据进行相关和后处理,通过相关、时延校准、条纹搜索,最终得到卫星的基线几何时延.基于天文开源软件建立起一套卫星干涉测量数据处理系统.该系统可工作在实时和事后两种状态,实现相关、中性大气、电离层、钟模型以及仪器硬件的时延校准、条纹搜索、生成基线时延和时延率序列.使用该系统处理北斗GEO (Geosynchronous Earth Orbit)卫星的干涉测量试验数据,得到了精度在1–2 ns量级的卫星基线时延序列.     

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⁻²  .     

Rotation in gravitational lenses

Gravitational lensing deflects light. A single lens deflector can only shear images, but cannot induce rotations. Multiple lens planes can induce rotations. Such rotations can be observed in quadruply imaged sources, and can be used to distinguish between two proposed solutions of the flux anomaly problem: substructures in lensing galaxies versus large-scale structure. We predict the expected amount of rotation due to large-scale structure in strong lensing systems, and show how this effect can be measured using ∼mas very long baseline interferometry astrometry of quadruple lenses with extended source structures. The magnitude of rotation is around 1°. The biggest theoretical uncertainty is the power spectrum of dark matter on very small scales. This procedure can potentially be turned around to measure the dark matter power spectrum on very small scales. We list the predicted rms rotation angles for several quadruple lenses with known lens and source redshifts.     

Reflected iron line from a source above a Kerr black hole accretion disc

In this paper we present a fully relativistic approach to modelling both the continuum emission and the reflected fluorescent iron line from a primary X-ray source near a Kerr black hole. The X-ray source is located above an accretion disc orbiting around the black hole. The source is assumed to be a static point source located on an arbitrary position above the disc, on or off the axis of rotation. We carry out Monte Carlo simulations in order to estimate the iron line spectrum as well as its equivalent width. Because of the gravitational lensing effect, an enhancement of the iron line is expected when the primary source is located close to the central black hole. We find that for a source located on the axis of rotation the enhancement is relatively modest. An observer at inclination 30° would measure an equivalent width of ∼300 eV in the extreme case of a maximally rotating black hole and a source located at height 1.5 gravitational radii from the centre. This corresponds to an equivalent width enhancement factor of about 2 compared with the classical value where no lensing effect comes into play. However, when allowing the source to be located off the axis of rotation, much stronger enhancement can be obtained. In the extreme case of a maximally rotating black hole and a source located just above the approaching side of the disc, an observer at inclination 30° could measure an equivalent width as high as ∼1.5 keV (i.e., ∼10 times the classical value). We also find that observers located at high inclination angles observe a stronger line than observers at low inclination angles.     

Baryonic acoustic oscillations in 21-cm emission: a probe of dark energy out to high redshifts

Low-frequency observatories are currently being constructed with the goal of detecting redshifted 21-cm emission from the epoch of reionization. These observatories will also be able to detect intensity fluctuations in the cumulative 21-cm emission after reionization, from hydrogen in unresolved damped Lyα absorbers (such as gas-rich galaxies) down to a redshift z ∼ 3.5. The inferred power spectrum of 21-cm fluctuations at all redshifts will show acoustic oscillations, whose comoving scale can be used as a standard ruler to infer the evolution of the equation of state for the dark energy. We find that the first generation of low-frequency experiments (such as MWA or LOFAR) will be able to constrain the acoustic scale to within a few per cent in a redshift window just prior to the end of the reionization era, provided that foregrounds can be removed over frequency bandpasses of ≳8 MHz. This sensitivity to the acoustic scale is comparable to the best current measurements from galaxy redshift surveys, but at much higher redshifts. Future extensions of the first-generation experiments (involving an order of magnitude increase in the antennae number of the MWA) could reach sensitivities below 1 per cent in several redshift windows and could be used to study the dark energy in the unexplored redshift regime of 3.5 ≲ z ≲ 12. Moreover, new experiments with antennae designed to operate at higher frequencies would allow precision measurements (≲1 per cent) of the acoustic peak to be made at more moderate redshifts (1.5 ≲ z ≲ 3.5), where they would be competitive with ambitious spectroscopic galaxy surveys covering more than 1000 deg². Together with other data sets, observations of 21-cm fluctuations will allow full coverage of the acoustic scale from the present time out to z ∼ 12.     

Probing subparsec structure in the Lyman α forest with gravitational microlensing

We present the results of microlens ray-tracing simulations showing the effect of absorbing material between a source quasar and a lensing galaxy in a gravitational lens system. We find that, in addition to brightness fluctuations due to microlensing, the strength of the absorption line relative to the continuum varies with time, with the properties of the variations depending on the structure of the absorbing material. We conclude that such variations will be measurable via ultraviolet spectroscopy of image A of the gravitationally lensed quasar Q2237+0305 if the Lyman α clouds between the quasar and the lensing galaxy possess structure on scales smaller than ∼0.1 pc. The time-scale for the variations is on the order of years to decades, although very short-term variability can occur. While the Lyman α lines may not be accessible at all wavelengths, this approach is applicable to any absorption system, including metal lines.     

Infrared observations of gravitational lensing in Abell 2219 with CIRSI

We present the first detection of a gravitational depletion signal at near-infrared wavelengths, based on deep panoramic images of the cluster Abell 2219 ( z =0.22) taken with the Cambridge Infrared Survey Instrument (CIRSI) at the prime focus of the 4.2-m William Herschel Telescope. Infrared studies of gravitational depletion offer a number of advantages over similar techniques applied at optical wavelengths, and can provide reliable total masses for intermediate-redshift clusters. Using the maximum-likelihood technique developed by Schneider, King & Erben, we detect the gravitational depletion at the 3 confidence level. By modelling the mass distribution as a singular isothermal sphere and ignoring the uncertainty in the unlensed number counts, we find an Einstein radius of (66 per cent confidence limit). This corresponds to a projected velocity dispersion of _v 800 km s¹, in agreement with constraints from strongly lensed features. For a Navarro, Frenk & White mass model, the radial dependence observed indicates a best-fitting halo scalelength of 125 h ¹ kpc. We investigate the uncertainties arising from the observed fluctuations in the unlensed number counts, and show that clustering is the dominant source of error. We extend the maximum-likelihood method to include the effect of incompleteness, and discuss the prospects of further systematic studies of lensing in the near-infrared band.     

Simulating a stochastic background of gravitational waves from neutron star formation at cosmological distances

We develop a temporal simulation of the potentially detectable gravitational wave background from neutron star formation at cosmological distances. By using a recent model for the evolving star formation rate, we investigate the statistical distribution of gravitational wave amplitudes due to supernovae that result in neutron star formation in the Einstein–de Sitter cosmology. We find that the gravitational wave amplitude distribution in our frame is highly skewed, with skewness related to the distribution of sources, and that the potentially detectable gravitational wave strain is dominated by sources at a redshift of     . Time traces of the simulation, using selected waveforms, are presented graphically and are also made available as web-based audio files. The method developed can readily be extended to different cosmologies, as well as to incorporate other waveforms and source types. This type of simulation will be useful in testing and optimizing detection strategies for gravitational wave backgrounds due to various types of individually undetectable astrophysical sources.     

On a systematic bias in surface brightness fluctuations based distances due to gravitational microlensing

The effect of gravitational microlensing on the determination of extragalactic distances using the surface brightness fluctuations (SBF) technique is considered and a method to calculate SBF amplitudes in the presence of microlensing is presented. With a simple approximation for the magnification power spectrum at low optical depth, the correction to the SBF-based luminosity distance is calculated. The results suggest the effect can be safely neglected at present but may become important for SBF-based Hubble diagrams at luminosity distances of about 1 Gpc and beyond.     

The statistics of weak lensing at small angular scales: probability distribution function

Weak gravitational lensing surveys have the potential to probe mass density fluctuation in the Universe directly. Recent studies have shown that it is possible to model the statistics of the convergence field at small angular scales by modelling the statistics of the underlying density field in the highly non-linear regime. We propose a new method to model the complete probability distribution function of the convergence field as a function of smoothing angle and source redshift. The model relies on a hierarchical ansatz for the behaviour of higher order correlations of the density field. We compare our results with ray-tracing simulations and find very good agreement over a range of smoothing angles. Whereas the density probability distribution function is not sensitive to the cosmological model, the probability distribution function for the convergence can be used to constrain both the power spectrum and cosmological parameters.     

Extended source effects in astrometric gravitational microlensing

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 a circularly symmetric limb-darkening profile, the centroid motion can be expressed as a one-dimensional integral, which can be evaluated numerically. We find 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 for determining stellar parameters, such as stellar radius, to excellent accuracy.     

General polytropic spheres as gravitational lenses

We propose hydrostatic polytropic spheres governed by the Lane-Emden equation (LEE) of index $n$ as a novel set of physical models for axially averaged gravitational lenses anywhere in the Universe, alternative to the familiar singular isothermal sphere (SIS) and the Navarro–Frenk–White (NFW) profile, as such general polytropic spheres are conceptually simple, versatile in representing a series of equations of state, and able to address both the inner core and cusp features. As LEE is nonlinear, there exist several distinct classes of LEE solutions to serve as physical lens models. With a few scaling parameters, the complete problem can be readily reconstructed with full physical dimensions. A given mass density profile satisfying LEE produces lensing effects that are solely determined by a dimensionless parameter $q$ which contains geometric and kinematic information about the source-lens-observer system. The lens mapping and tangential shear or distortion profile are derived, first analytically for special cases and then asymptotically at the outskirts or near the edge of the lens. Numerical procedures for calculating full lensing profiles of a general lens are developed. Our results include the analytical “singular polytropic sphere” (SPS) profile which generalizes the SIS model and may outperform the latter in modeling dark matter halos among others. We further point out that dynamic models of general polytropic spheres in self-similar evolution can serve as several broad classes of gravitational lenses and produce time-dependent lensing effects slow or fast depending on the pertinent time scales. Astrophysical sources that can be lensed include electromagnetic wave sources in the entire frequency band, gravitational wave sources in the entire frequency band, gravitons even possibly with finite masses, neutrino sources of three different types, neutron sources, and ultra high energy cosmic rays (UHECRs) of electrically charged particles which can also interact with magnetic fields. We discuss and elabrate applications to dark matter halos, hypermassive black holes and supermassive black holes in the entire Universe including the early Universe, magnetized supermassive stars, static and dynamically evolving spherical and cylindrical lenses in contexts of astrophysics and cosmology.     

Planetesimals to protoplanets – II. Effect of debris on terrestrial planet formation

In this paper, we extend our numerical method for simulating terrestrial planet formation to include dynamical friction from the unresolved debris component. In the previous work, we implemented a rubble pile planetesimal collision model into direct N -body simulations of terrestrial planet formation. The new collision model treated both accretion and erosion of planetesimals but did not include dynamical friction from debris particles smaller than the resolution limit for the simulation. By extending our numerical model to include dynamical friction from the unresolved debris, we can simulate the dynamical effect of debris produced during collisions and can also investigate the effect of initial debris mass on terrestrial planet formation. We find that significant initial debris mass, 10 per cent or more of the total disc mass, changes the mode of planetesimal growth. Specifically, planetesimals in this situation do not go through a runaway growth phase. Instead, they grow concurrently, similar to oligarchic growth. The dynamical friction from the unresolved debris damps the eccentricities of the planetesimals, reducing the mean impact speeds and causing all collisions to result in merging with no mass loss. As a result, there is no debris production. The mass in debris slowly decreases with time. In addition to including the dynamical friction from the unresolved debris, we have implemented particle tracking as a proxy for monitoring compositional mixing. Although there is much less mixing due to collisions and gravitational scattering when dynamical friction of the background debris is included, there is significant inward migration of the largest protoplanets in the most extreme initial conditions (for which the initial mass in unresolved debris is at least equal to the mass in resolved planetesimals).     

OJ287双黑洞轨道计算:3.5阶后牛顿近似

大质量双黑洞OJ287是一个强引力辐射源.为了探测其引力波信号,需要知道波形,而这主要是由轨道运动所决定.为此,从广义相对论3.5阶后牛顿近似的运动方程出发对OJ287的轨道进行仔细研究,取大黑洞位置固定作为近似,给出了后牛顿近似下3.5阶的次黑洞轨道解,比他人2.5阶的工作高了一阶.次黑洞撞击吸积盘面到光学爆发存在时间延迟,这对于确定轨道参数有很大影响.利用径向距离与爆发时间关系的线性模型,对最近7次爆发时刻的观测值拟合,给出了更精确的OJ287双黑洞的轨道参数及其运动轨道.分析了计算结果,研究运动特征,并且发现了两个新性质:次黑洞进动在初期增加,在晚期接近并和时,进动达到最大值,然后减小并越过0而趋于负值.尚不能确定晚期的这个行为是否由3.5阶近似不够准确所造成.运动方程中耗散性的辐射项,后牛顿2.5阶和3.5阶的系数具有相反的符号.这意味着3.5阶项反而是从外界吸收能量.但2.5阶与3.5阶之和仍然是向外辐射引力波的,体系能量变化率为负.这个工作的计算结果可以用来更精确地计算OJ287的引力辐射.     

Theoretical estimates of intrinsic galaxy alignment

It has recently been argued that the observed ellipticities of galaxies may be determined at least in part by the primordial tidal gravitational field in which the galaxy formed. Long-range correlations in the tidal field could thus lead to an ellipticity–ellipticity correlation for widely separated galaxies. We present a new model relating ellipticity to angular momentum, which can be calculated in linear theory. We use this model to calculate the angular power spectrum of intrinsic galaxy shape correlations. We show that, for low-redshift galaxy surveys, our model predicts that intrinsic correlations will dominate correlations induced by weak lensing, in good agreement with previous theoretical work and observations. We find that our model produces ' E -mode' correlations enhanced by a factor of 3.5 over B -modes on small scales, making it harder to disentangle intrinsic correlations from those induced by weak gravitational lensing.     

The statistics of wide-separation lensed quasars

The absence of any wide-separation gravitational lenses in the Large Bright Quasar Survey is used to place limits on the population of cluster-sized haloes in the universe, and hence constrain a number of cosmological parameters. The results agree with previous investigations in strongly ruling out the standard cold dark matter model but they are consistent with low-density universes in which the primordial fluctuation spectrum matches both cluster abundances and cosmic microwave background measurements. These conclusions are essentially independent of the cosmological constant, which is in stark contrast to the statistics of galaxy lenses. The constraints presented here are nullified if clusters have core radii of ≳10 kpc, but are free of a number of potential systematic errors, owing to the homogeneity of the data.     

Two new large-separation gravitational lenses from SDSS

We present discovery images, together with follow-up imaging and spectroscopy, of two large-separation gravitational lenses found by our survey for wide arcs [the CAmbridge Sloan Survey Of Wide ARcs in the skY (CASSOWARY)]. The survey exploits the multicolour photometry of the Sloan Digital Sky Survey to find multiple blue components around red galaxies. CASSOWARY 2 (or 'the Cheshire Cat') is composed of two massive early-type galaxies at   z = 0.426  and 0.432, respectively, lensing two background sources, the first a star-forming galaxy at   z = 0.97  and the second a high -redshift galaxy  ( z > 1.4)  . There are at least three images of the former source and probably four or more of the latter, arranged in two giant arcs. The mass enclosed within the larger arc of radius ∼11 arcsec is  ∼33 × 10¹² M_⊙  . CASSOWARY 3 comprises an arc of three bright images of a   z = 0.725  source, lensed by a foreground elliptical at   z = 0.274  . The radius of the arc is ∼4 arcsec and the enclosed mass is  ∼2.5 × 10¹² M_⊙  . Together with earlier discoveries like the Cosmic Horseshoe and the 8 o'clock Arc, these new systems, with separations intermediate between the arcsecond-separation lenses of typical strong galaxy lensing and arcminute-separation cluster lenses, probe the very high end of the galaxy mass function.     

A genetic algorithm for the non-parametric inversion of strong lensing systems

We present a non-parametric technique to infer the projected mass distribution of a gravitational lens system with multiple strong-lensed images. The technique involves a dynamic grid in the lens plane on which the mass distribution of the lens is approximated by a sum of basis functions, one per grid cell. We used the projected mass densities of Plummer spheres as basis functions. A genetic algorithm then determines the mass distribution of the lens by forcing images of a single source, projected back on to the source plane, to coincide as well as possible. Averaging several tens of solutions removes the random fluctuations that are introduced by the reproduction process of genomes in the genetic algorithm and highlights those features common to all solutions. Given the positions of the images and the redshifts of the sources and the lens, we show that the mass of a gravitational lens can be retrieved with an accuracy of a few percent and that, if the sources sufficiently cover the caustics, the mass distribution of the gravitational lens can also be reliably retrieved. A major advantage of the algorithm is that it makes full use of the information contained in the radial images, unlike methods that minimize the residuals of the lens equation, and is thus able to accurately reconstruct also the inner parts of the lens.