Source author record

Prasenjit Saha

Prasenjit Saha appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

38works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

38 published item(s)

preprint2023arXiv

What are the parities of photon-ring images near a black hole?

Light that grazes a black-hole event horizon can loop around one or more times before escaping again, resulting for distance observers in an infinite sequence of ever fainter and more delayed images near the black hole shadow. In the case of the M87 and Sgr A$^*$ back holes, the first of these so-called photon-ring images have now been observed. A question then arises: are such images minima, maxima, or saddle-points in the sense of Fermat's principle in gravitational lensing? or more briefly, the title question above. In the theory of lensing by weak gravitational fields, image parities are readily found by considering the time-delay surface (also called the Fermat potential or the arrival-time surface). In this work, we extend the notion of the time delay surface to strong gravitational fields and compute the surface for a Schwarzschild black hole. The time-delay surface is the difference of two wavefronts, one travelling forward from the source and one travelling backwards from the observer. Image parities are read off from the topography of the surface, exactly as in the weak-field regime, but the surface itself is more complicated. Of the images, furthest from the black hole and similar to the weak-field limit, are a minimum and a saddle point. The strong field repeats the pattern, corresponding to light taking one or more loops around the back hole. In between, there are steeply-rising walls in the time-delay surface, which can be interpreted as maxima and saddle points that are infinitely delayed and not observable -- these correspond to light rays taking a U-turn around the black hole.

preprint2022arXiv

Observed versus Simulated Halo c-Mvir Relations

The concentration - virial mass relation is a well-defined trend that reflects the formation of structure in an expanding Universe. Numerical simulations reveal a marked correlation that depends on the collapse time of dark matter halos and their subsequent assembly history. However, observational constraints are mostly limited to the massive end via X-ray emission of the hot diffuse gas in clusters. An alternative approach, based on gravitational lensing over galaxy scales, revealed an intriguingly high concentration at Milky Way-sized halos. This letter focuses on the robustness of these results by adopting a bootstrapping approach that combines stellar and lensing mass profiles. We also apply the identical methodology to simulated halos from EAGLE to assess any systematic. We bypass several shortcomings of ensemble type lens reconstruction and conclude that the mismatch between observed and simulated concentration-to-virial-mass relations are robust, and need to be explained either invoking a lensing-related sample selection bias, or a careful investigation of the evolution of concentration with assembly history. For reference, at a halo mass of $10^{12} M_\odot$, the concentration of observed lenses is $c_{12}\sim 40\pm 5$, whereas simulations give $c_{12}\sim 15\pm1$.

preprint2022arXiv

Simulations of astrometric planet detection in Alpha Centauri by intensity interferometry

Recent dynamical studies indicate that the possibility of an Earth-like planet around $α\;$Cen A or B should be taken seriously. Such a planet, if it exists, would perturb the orbital astrometry by $<10 \ μ\rm as$, which is $10^{-6}$ of the separation between the two stars. We assess the feasibility of detecting such perturbations using ground-based intensity interferometry. We simulate a dedicated setup consisting of four 40-cm telescopes equipped with photon counters and correlators with time resolution $0.1\,\rm ns$, and a sort of matched filter implemented through an aperture mask. The astrometric error from one night of observing $α\;$Cen AB is $\approx0.5\,\rm mas$. The error decreases if longer observing times and multiple spectral channels are used, as $(\hbox{channels}\times\hbox{nights})^{-1/2}$.

preprint2021arXiv

The Hubble constant from eight time-delay galaxy lenses

We present a determination of the Hubble constant from the joint, free-form analysis of 8 strongly, quadruply lensing systems. In the concordance cosmology, we find $H_0 = 71.8^{+3.9}_{-3.3}\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$ with a precision of $4.97\%$. This is in agreement with the latest measurements from Supernovae Type Ia and Planck observations of the cosmic microwave background. Our precision is lower compared to these and other recent time-delay cosmography determinations, because our modelling strategies reflect the systematic uncertainties of lensing degeneracies. We furthermore are able to find reasonable lensed image reconstructions by constraining to either value of $H_0$ from local and early Universe measurements. This leads us to conclude that current lensing constraints on $H_0$ are not strong enough to break the "Hubble tension" problem of cosmology.

preprint2020arXiv

Lessons from a blind study of simulated lenses: image reconstructions do not always reproduce true convergence

In the coming years, strong gravitational lens discoveries are expected to increase in frequency by two orders of magnitude. Lens-modelling techniques are being developed to prepare for the coming massive influx of new lens data, and blind tests of lens reconstruction with simulated data are needed for validation. In this paper we present a systematic blind study of a sample of 15 simulated strong gravitational lenses from the EAGLE suite of hydrodynamic simulations. We model these lenses with a free-form technique and evaluate reconstructed mass distributions using criteria based on shape, orientation, and lensed image reconstruction. Especially useful is a lensing analogue of the Roche potential in binary star systems, which we call the $\textit{lensing Roche potential}$. This we introduce in order to factor out the well-known problem of steepness or mass-sheet degeneracy. Einstein radii are on average well recovered with a relative error of ${\sim}5\%$ for quads and ${\sim}25\%$ for doubles; the position angle of ellipticity is on average also reproduced well up to $\pm10^{\circ}$, but the reconstructed mass maps tend to be too round and too shallow. It is also easy to reproduce the lensed images, but optimising on this criterion does not guarantee better reconstruction of the mass distribution.

preprint2020arXiv

The theory of intensity interferometry revisited

With the current revival of interest in astronomical intensity interferometry, it is interesting to revisit the associated theory, which was developed in the 1950s and 1960s. This paper argues that intensity interferometry can be understood as an extension of Fraunhofer diffraction to incoherent light. Interference patterns are still produced, but they are speckle-like and transient, changing on a time scale of $1/Δν$ (where $Δν$ is the frequency bandwidth) known as the coherence time. Bright fringes average less than one photon per coherence time, hence fringes change before they can be observed. But very occasionally, two or even more photons may be detected from an interference pattern within a coherence time. These rare coincident photons provide information about the underlying transient interference pattern, and hence about the source brightness distribution. Thinking in terms of transient sub-photon interference patterns makes it easy to see why intensity interferometry will have large optical-path tolerance, and be immune to atmospheric seeing. The unusual signal-to-noise properties also become evident. We illustrate the unobservable but conceptually useful transient interference patterns, and their observable correlation signal, with three simulated examples: (i) an elongated source like Achernar, (ii) a three-star system like Algol, and (iii) a crescent source that roughly mimics an exoplanet transit or perhaps the M87 black hole environment. Of these, (i) and (ii) are good targets for currently-planned setups, while (iii) is interesting to think about for the longer term.

preprint2020arXiv

Towards a Polarisation Prediction for LISA via Intensity Interferometry

Compact Galactic binary systems with orbital periods of a few hours are expected to be detected in gravitational waves (GW) by LISA or a similar mission. At present, these so-called verification binaries provide predictions for GW frequency and amplitude. A full polarisation prediction would provide a new method to calibrate LISA and other GW observatories, but requires resolving the orientation of the binary on the sky, which is not currently possible. We suggest a method to determine the elusive binary orientation and hence predict the GW polarisation, using km-scale optical intensity interferometry. The most promising candidate is CD-30$^{\circ}$ 11223, consisting of a hot helium subdwarf with $m_B = 12$ and a much fainter white dwarf companion, in a nearly edge-on orbit with period 70.5 min. We estimate that the brighter star is tidally stretched by 6%. Resolving the tidal stretching would provide the binary orientation. The resolution needed is far beyond any current instrument, but not beyond current technology. We consider scenarios where an array of telescopes with km-scale baselines and/or the Very Large Telescope (VLT) and Extremely Large Telescope (ELT) are equipped with recently-developed kilo-pixel sub-ns single-photon counters and used for intensity interferometry. We estimate that a team-up of the VLT and ELT could measure the orientation to $\pm 1^{\circ}$ at 2$σ$ confidence in 24 hours of observation.

preprint2019arXiv

Astronomy and the new SI

In 2019 the International System of units (SI) conceptually re-invented itself. This was necessary because quantum-electronic devices had become so precise that the old SI could no longer calibrate them. The new system defines values of fundamental constants (including $c,h,k,e$ but not $G$) and allows units to be realized from the defined constants through any applicable equation of physics. In this new and more abstract SI, units can take on new guises --- for example, the kilogram is at present best implemented as a derived electrical unit. Relevant to astronomy, however, is that several formerly non-SI units, such as electron-volts, light-seconds, and what we may call "gravity seconds" $GM/c^3$, can now be interpreted not as themselves units, but as shorthand for volts and seconds being used with particular equations of physics. Moreover, the classical astronomical units have exact and rather convenient equivalents in the new SI: zero AB magnitude amounts to $\simeq5\times10^{10}$ photons $\rm m^{-2}\,s^{-1}$ per logarithmic frequency or wavelength interval, $\rm 1\,au\simeq 500$ light-seconds, $\rm 1\,pc\simeq 10^8$ light-seconds, while a solar mass $\simeq5$ gravity-seconds. As a result, the unit conversions ubiquitous in astrophysics can now be eliminated, without introducing other problems, as the old-style SI would have done. We review a variety of astrophysical processes illustrating the simplifications possible with the new-style SI, with special attention to gravitational dynamics, where care is needed to avoid propagating the uncertainty in $G$. Well-known systems (GPS satellites, GW170817, and the M87 black hole) are used as examples wherever possible.

preprint2016arXiv

Microlensing as a possible probe of event-horizon structure in quasars

In quasars which are lensed by galaxies, the point-like images sometimes show sharp and uncorrelated brightness variations (microlensing). These brightness changes are associated with the innermost region of the quasar passing through a complicated pattern of caustics produced by the stars in the lensing galaxy. In this paper, we study whether the universal properties of optical caustics could enable extraction of shape information about the central engine of quasars. We present a toy model with a crescent-shaped source crossing a fold caustic. The silhouette of a black hole over an accretion disk tends to produce roughly crescent sources. When a crescent-shaped source crosses a fold caustic, the resulting light curve is noticeably different from the case of a circular luminosity profile or Gaussian source. With good enough monitoring data, the crescent parameters, apart from one degeneracy, can be recovered.

preprint2016arXiv

Quantifying substructures in {\it Hubble Frontier Field} clusters: comparison with $ΛCDM$ simulations

The Hubble Frontier Fields (HFF) are six clusters of galaxies, all showing indications of recent mergers, which have recently been observed for lensed images. As such they are the natural laboratories to study the merging history of galaxy clusters. In this work, we explore the 2D power spectrum of the mass distribution $P_{\rm M}(k)$ as a measure of substructure. We compare $P_{\rm M}(k)$ of these clusters (obtained using strong gravitational lensing) to that of $Λ$CDM simulated clusters of similar mass. To compute lensing $P_{\rm M}(k)$, we produced free-form lensing mass reconstructions of HFF clusters, without any light traces mass (LTM) assumption. The inferred power at small scales tends to be larger if (i)~the cluster is at lower redshift, and/or (ii)~there are deeper observations and hence more lensed images. In contrast, lens reconstructions assuming LTM show higher power at small scales even with fewer lensed images; it appears the small scale power in the LTM reconstructions is dominated by light information, rather than the lensing data. The average lensing derived $P_{\rm M}(k)$ shows lower power at small scales as compared to that of simulated clusters at redshift zero, both dark-matter only and hydrodynamical. The possible reasons are: (i)~the available strong lensing data are limited in their effective spatial resolution on the mass distribution, (ii)~HFF clusters have yet to build the small scale power they would have at $z\sim 0$, or (iii)~simulations are somehow overestimating the small scale power.

preprint2016arXiv

Strong Gravitational Lensing and the Stellar IMF of Early-type Galaxies

Systematic variations of the IMF in early-type galaxies, and their connection with possible drivers such as velocity dispersion or metallicity, have been much debated in recent years. Strong lensing over galaxy scales combined with photometric and spectroscopic data provides a powerful method to constrain the stellar mass-to-light ratio and hence the functional form of the IMF. We combine photometric and spectroscopic constraints from the latest set of population synthesis models of Charlot & Bruzual, including a varying IMF, with a non-parametric analysis of the lens masses of 18 ETGs from the SLACS survey, with velocity dispersions in the range 200-300 km/s. We find that very bottom-heavy IMFs are excluded. However, the upper limit to the bimodal IMF slope ($μ\lesssim 2.2$, accounting for a dark matter fraction of 20-30%, where $μ=1.3$ corresponds to a Kroupa-like IMF) is compatible at the $1σ$ level with constraints imposed by gravity-sensitive line strengths. A two-segment power law parameterisation of the IMF (Salpeter-like for high masses) is more constrained ($Γ\lesssim 1.5$, where $Γ$ is the power index at low masses) but requires a dark matter contribution of $\gtrsim 25\%$ to reconcile the results with a Salpeter IMF. For a standard Milky Way-like IMF to be applicable, a significant dark matter contribution is required within $1R_e$. Our results reveal a large range of allowed IMF slopes, which, when interpreted as intrinsic scatter in the IMF properties of ETGs, could explain the recent results of Smith et al., who find Milky Way-like IMF normalisations in a few massive lensing ETGs.

preprint2016arXiv

Testing light-traces-mass in Hubble Frontier Fields Cluster MACS-J0416.1-2403

We reconstruct the projected mass distribution of a massive merging Hubble Frontier Fields cluster MACSJ0416 using the genetic algorithm based free-form technique called Grale. The reconstructions are constrained by 149 lensed images identified by Jauzac et al. using HFF data. No information about cluster galaxies or light is used, which makes our reconstruction unique in this regard. Using visual inspection of the maps, as well as galaxy-mass correlation functions we conclude that overall light does follow mass. Furthermore, the fact that brighter galaxies are more strongly clustered with mass is an important confirmation of the standard biasing scenario in galaxy clusters. On the smallest scales, approximately less than a few arcseconds, the resolution afforded by 149 images is still not sufficient to confirm or rule out galaxy-mass offsets of the kind observed in ACO 3827. We also compare the mass maps of MACSJ0416 obtained by three different groups: Grale, and two parametric Lenstool reconstructions from the CATS and Sharon/Johnson teams. Overall, the three agree well; one interesting discrepancy between Grale and Lenstool galaxy-mass correlation functions occurs on scales of tens of kpc and may suggest that cluster galaxies are more biased tracers of mass than parametric methods generally assume.

preprint2015arXiv

Geometrical vs wave optics under gravitational waves

We present some new derivations of the effect of a plane gravitational wave on a light ray. A simple interpretation of the results is that a gravitational wave causes a phase modulation of electromagnetic waves. We arrive at this picture from two contrasting directions, namely null geodesics and Maxwell's equations, or, geometric and wave optics. Under geometric optics, we express the geodesic equations in Hamiltonian form and solve perturbatively for the effect of gravitational waves. We find that the well-known time-delay formula for light generalizes trivially to massive particles. We also recover, by way of a Hamilton-Jacobi equation, the phase modulation obtained under wave optics. Turning then to wave optics - rather than solving Maxwell's equations directly for the fields, as in most previous approaches - we derive a perturbed wave equation (perturbed by the gravitational wave) for the electromagnetic four-potential. From this wave equation it follows that the four-potential and the electric and magnetic fields all experience the same phase modulation. Applying such a phase modulation to a superposition of plane waves corresponding to a Gaussian wave packet leads to time delays.

preprint2015arXiv

Gravitational lens modelling in a citizen science context

We develop a method to enable collaborative modelling of gravitational lenses and lens candidates, that could be used by non-professional lens enthusiasts. It uses an existing free-form modelling program (glass), but enables the input to this code to be provided in a novel way, via a user-generated diagram that is essentially a sketch of an arrival-time surface. We report on an implementation of this method, SpaghettiLens, which has been tested in a modelling challenge using 29 simulated lenses drawn from a larger set created for the Space Warps citizen science strong lens search. We find that volunteers from this online community asserted the image parities and time ordering consistently in some lenses, but made errors in other lenses depending on the image morphology. While errors in image parity and time ordering lead to large errors in the mass distribution, the enclosed mass was found to be more robust: the model-derived Einstein radii found by the volunteers were consistent with those produced by one of the professional team, suggesting that given the appropriate tools, gravitational lens modelling is a data analysis activity that can be crowd-sourced to good effect. Ideas for improvement are discussed, these include (a) overcoming the tendency of the models to be shallower than the correct answer in test cases, leading to systematic overestimation of the Einstein radius by 10 per cent at present, and (b) detailed modelling of arcs.

preprint2015arXiv

Space Warps II. New Gravitational Lens Candidates from the CFHTLS Discovered through Citizen Science

We report the discovery of 29 promising (and 59 total) new lens candidates from the CFHT Legacy Survey (CFHTLS) based on about 11 million classifications performed by citizen scientists as part of the first Space Warps lens search. The goal of the blind lens search was to identify lens candidates missed by robots (the RingFinder on galaxy scales and ArcFinder on group/cluster scales) which had been previously used to mine the CFHTLS for lenses. We compare some properties of the samples detected by these algorithms to the Space Warps sample and find them to be broadly similar. The image separation distribution calculated from the Space Warps sample shows that previous constraints on the average density profile of lens galaxies are robust. SpaceWarps recovers about 65% of known lenses, while the new candidates show a richer variety compared to those found by the two robots. This detection rate could be increased to 80% by only using classifications performed by expert volunteers (albeit at the cost of a lower purity), indicating that the training and performance calibration of the citizen scientists is very important for the success of Space Warps. In this work we present the SIMCT pipeline, used for generating in situ a sample of realistic simulated lensed images. This training sample, along with the false positives identified during the search, has a legacy value for testing future lens finding algorithms. We make the pipeline and the training set publicly available.

preprint2015arXiv

Space Warps: I. Crowd-sourcing the Discovery of Gravitational Lenses

We describe Space Warps, a novel gravitational lens discovery service that yields samples of high purity and completeness through crowd-sourced visual inspection. Carefully produced colour composite images are displayed to volunteers via a web- based classification interface, which records their estimates of the positions of candidate lensed features. Images of simulated lenses, as well as real images which lack lenses, are inserted into the image stream at random intervals; this training set is used to give the volunteers instantaneous feedback on their performance, as well as to calibrate a model of the system that provides dynamical updates to the probability that a classified image contains a lens. Low probability systems are retired from the site periodically, concentrating the sample towards a set of lens candidates. Having divided 160 square degrees of Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) imaging into some 430,000 overlapping 82 by 82 arcsecond tiles and displaying them on the site, we were joined by around 37,000 volunteers who contributed 11 million image classifications over the course of 8 months. This Stage 1 search reduced the sample to 3381 images containing candidates; these were then refined in Stage 2 to yield a sample that we expect to be over 90% complete and 30% pure, based on our analysis of the volunteers performance on training images. We comment on the scalability of the SpaceWarps system to the wide field survey era, based on our projection that searches of 10$^5$ images could be performed by a crowd of 10$^5$ volunteers in 6 days.

preprint2015arXiv

The behaviour of dark matter associated with 4 bright cluster galaxies in the 10kpc core of Abell 3827

Galaxy cluster Abell 3827 hosts the stellar remnants of four almost equally bright elliptical galaxies within a core of radius 10kpc. Such corrugation of the stellar distribution is very rare, and suggests recent formation by several simultaneous mergers. We map the distribution of associated dark matter, using new Hubble Space Telescope imaging and VLT/MUSE integral field spectroscopy of a gravitationally lensed system threaded through the cluster core. We find that each of the central galaxies retains a dark matter halo, but that (at least) one of these is spatially offset from its stars. The best-constrained offset is 1.62+/-0.48kpc, where the 68% confidence limit includes both statistical error and systematic biases in mass modelling. Such offsets are not seen in field galaxies, but are predicted during the long infall to a cluster, if dark matter self-interactions generate an extra drag force. With such a small physical separation, it is difficult to definitively rule out astrophysical effects operating exclusively in dense cluster core environments - but if interpreted solely as evidence for self-interacting dark matter, this offset implies a cross-section sigma/m=(1.7+/-0.7)x10^{-4}cm^2/g x (t/10^9yrs)^{-2}, where t is the infall duration.

preprint2014arXiv

Clocks around Sgr A*

The S stars near the Galactic centre and any pulsars that may be on similar orbits, can be modelled in a unified way as clocks orbiting a black hole, and hence are potential probes of relativistic effects, including black hole spin. The high eccentricities of many S stars mean that relativistic effects peak strongly around pericentre; for example, orbit precession is not a smooth effect but almost a kick at pericentre. We argue that concentration around pericentre will be an advantage when analysing redshift or pulse-arrival data to measure relativistic effects, because cumulative precession will be drowned out by Newtonian perturbations from other mass in the Galactic-centre region. Wavelet decomposition may be a way to disentangle relativistic effects from Newton perturbations. Assuming a plausible model for Newtonian perturbations on S2, relativity appears to be strongest in a two-year interval around pericentre, in wavelet modes of timescale approximately 6 months.

preprint2014arXiv

Feasibility of observing Hanbury Brown and Twiss phase

The interferometers of Hanbury Brown and collaborators in the 1950s and 60s, and their modern descendants now being developed (intensity interferometers) measure the spatial power spectrum of the source from intensity correlations at two points. The quantum optical theory of the Hanbury Brown and Twiss (HBT) effect shows that more is possible, in particular the phase information can be recovered by correlating intensities at three points (bispectrum). In this paper we argue that such 3 point measurements are possible for bright stars such as Sirius and Betelgeuse using off the shelf single photon counters with collecting areas of the order of 100m2. It seems possible to map individual features on the stellar surface. Simple diameter measurements would be possible with amateur class telescopes.

preprint2014arXiv

Gravitational Lens Recovery with GLASS: Measuring the mass profile and shape of a lens

We use a new non-parametric gravitational modelling tool -- \Glass{} -- to determine what quality of data (strong lensing, stellar kinematics, and/or stellar masses) are required to measure the circularly averaged mass profile of a lens and its shape. \Glass{} uses an under-constrained adaptive grid of mass pixels to model the lens, searching through thousands of models to marginalise over model uncertainties. Our key findings are as follows: (i) for pure lens data, multiple sources with wide redshift separation give the strongest constraints as this breaks the well-known mass-sheet or steepness degeneracy; (ii) a single quad with time delays also performs well, giving a good recovery of both the mass profile and its shape; (iii) stellar masses -- for lenses where the stars dominate the central potential -- can also break the steepness degeneracy, giving a recovery for doubles almost as good as having a quad with time delay data, or multiple source redshifts; (iv) stellar kinematics provide a robust measure of the mass at the half light radius of the stars $r_{1/2}$ that can also break the steepness degeneracy if the Einstein radius $r_E \neq r_{1/2}$; and (v) if $r_E \sim r_{1/2}$, then stellar kinematic data can be used to probe the stellar velocity anisotropy $β$ -- an interesting quantity in its own right. Where information on the mass distribution from lensing and/or other probes becomes redundant, this opens up the possibility of using strong lensing to constrain cosmological models.

preprint2014arXiv

Lensing time delays as a substructure constraint: a case study with the cluster SDSS~J1004+4112

Gravitational lensing time delays are well known to depend on cosmological parameters, but they also depend on the details of the mass distribution of the lens. It is usual to model the mass distribution and use time-delay observations to infer cosmological parameters, but it is naturally also possible to take the cosmological parameters as given and use time delays as constraints on the mass distribution. This paper develops a method to isolate what exactly those constraints are, using a principal-components analysis of ensembles of free-form mass models. We find that time delays provide tighter constraints on the distribution of matter in the very high dense regions of the lensing clusters. We apply it to the cluster lens SDSS J1004+4112, whose rich lensing data includes two time delays. We find, assuming a concordance cosmology, that the time delays constrain the central region of the cluster to be rounder and less lopsided than would be allowed by lensed images alone. This detailed information about the distribution of the matter is very useful for studying the dense regions of the galaxy clusters which are very difficult to study with direct measurements. A further time-delay measurement, which is expected, will make this system even more interesting.

preprint2014arXiv

Mass-Galaxy offsets in Abell 3827, 2218 and 1689: intrinsic properties or line-of-sight substructures?

We have made mass maps of three strong-lensing clusters, Abell 3827, Abell 2218 and Abell 1689, in order to test for mass-light offsets. The technique used is GRALE, which enables lens reconstruction with minimal assumptions, and specifically with no information about the cluster light being given. In the first two of these clusters, we find local mass peaks in the central regions that are displaced from the nearby galaxies by a few to several kpc. These offsets {\em could\/} be due to line of sight structure unrelated to the clusters, but that is very unlikely, given the typical levels of chance line-of-sight coincidences in $ΛCDM$ simulations --- for Abell 3827 and Abell 2218 the offsets appear to be intrinsic. In the case of Abell 1689, we see no significant offsets in the central region, but we do detect a possible line of sight structure: it appears only when sources at $z\ga 3$ are used for reconstructing the mass. We discuss possible origins of the mass-galaxy offsets in Abell 3827 and Abell 2218: these include pure gravitational effects like dynamical friction, but also non-standard mechanisms like self-interacting dark-matter.

preprint2014arXiv

Spacecraft Clocks and Relativity: Prospects for Future Satellite Missions

The successful miniaturization of extremely accurate atomic clocks invites prospects for satellite missions to perform precise timing experiments. This will allow effects predicted by general relativity to be detected in Earth's gravitational field. In this paper we introduce a convenient formalism for studying these effects, and compute the fractional timing differences generated by them for the orbit of a satellite capable of accurate time transfer to a terrestrial receiving station on Earth, as proposed by planned missions. We find that (1) Schwarzschild perturbations would be measurable through their effects both on the orbit and on the signal propagation, (2) frame-dragging of the orbit would be readily measurable, and (3) in optimistic scenarios, the spin-squared metric effects may be measurable for the first time ever. Our estimates suggest that a clock with a fractional timing inaccuracy of $10^{-16}$ on a highly eccentric Earth orbit will measure all these effects, while for a low Earth circular orbit like that of the Atomic Clock Ensemble in Space Mission, detection will be more challenging.

preprint2014arXiv

Testing General Relativity and Alternative Theories of Gravity with Space-based Atomic Clocks and Atom Interferometers

The successful miniaturisation of extremely accurate atomic clocks and atom interferometers invites prospects for satellite missions to perform precision experiments. We discuss the effects predicted by general relativity and alternative theories of gravity that can be detected by a clock, which orbits the Earth. Our experiment relies on the precise tracking of the spacecraft using its observed tick-rate. The spacecraft's reconstructed four-dimensional trajectory will reveal the nature of gravitational perturbations in Earth's gravitational field, potentially differentiating between different theories of gravity. This mission can measure multiple relativistic effects all during the course of a single experiment, and constrain the Parametrized Post-Newtonian Parameters around the Earth. A satellite carrying a clock of fractional timing inaccuracy of $Δf/f \sim 10^{-16}$ in an elliptic orbit around the Earth would constrain the PPN parameters $|β-1|, |γ-1| \lesssim 10^{-6}$. We also briefly review potential constraints by atom interferometers on scalar tensor theories and in particular on Chameleon and dilaton models.

preprint2013arXiv

Intensity interferometry with more than two detectors?

The original intensity interferometers were instruments built in the 1950s and 60s by Hanbury Brown and collaborators, achieving milli-arcsec resolutions in visible light without optical-quality mirrors. They exploited a then-novel physical effect, now known as HBT correlation after the experiments of Hanbury Brown and Twiss, and nowadays considered fundamental in quantum optics. Now a new generation of inten- sity interferometers is being designed, raising the possibility of measuring intensity correlations with three or more detectors. Quantum optics predicts some interesting features in higher-order HBT. One is that HBT correlation increases combinatorially with the number of detectors. Signal to noise considerations suggest, that many-detector HBT correlations would be mea- surable for bright masers, but very difficult for thermal sources. But the more modest three-detector HBT correlation seems measurable for bright stars, and would provide image information (namely the bispectrum) not present in standard HBT.

preprint2012arXiv

Diagnostics of Baryonic Cooling in Lensing Galaxies

Theoretical studies of structure formation find an inverse proportionality between the concentration of dark matter haloes and virial mass. This trend has been recently confirmed for virial masses Mvir > ~6e12 Msun by the observation of the X-ray emission from the hot halo gas. We present an alternative approach to this problem, exploring the concentration of dark matter haloes over galaxy scales on a sample of 18 early-type systems. Our c-Mvir relation is consistent with the X-ray analysis, extending towards lower virial masses, covering the range from ~4e11 Msun up to 5e12 Msun. A combination of the lensing analysis along with photometric data allows us to constrain the baryon fraction within a few effective radii, which is compared with prescriptions for adiabatic contraction (AC) of the dark matter haloes. We find that the standard methods for AC are strongly disfavored, requiring additional mechanisms -- such as mass loss during the contraction process -- to play a role during the phases following the collapse of the haloes.

preprint2011arXiv

Light/Mass Offsets in the Lensing Cluster Abell 3827: Evidence for Collisional Dark Matter?

If dark matter has a non-zero self-interaction cross-section, then dark matter halos of individual galaxies in cluster cores should experience a drag force from the ambient dark matter of the cluster, which will not affect the stellar components of galaxies, and thus will lead to a separation between the stellar and dark matter. If the cross-section is only a few decades below its current astrophysically determined upper limit, then kpc-scale separations should result. However, such separations will be observable only under very favorable conditions. Abell 3827 is a nearby late stage cluster merger with four massive central ellipticals within 20 kpc of each other. The ten strong lensing images tightly surrounding the ellipticals provide an excellent set of constraints for a free-form lens reconstruction. Our free-form mass maps show a massive dark extended clump, about 6 kpc from one of the ellipticals. The robustness of this result has been tested with many reconstructions, and confirmed with experiments using synthetic lens mass distributions. Interpreted in terms of dark matter collisionality, our result yields sigma/m ~ 4.5 10^{-7} (t/10^{10} yr)^{-2} cm^2/g, where t is the merger's age.

preprint2011arXiv

Resolving the baryon-fraction profile in lensing galaxies

We show the radial dependence of stellar baryon fraction curves derived for 21 lensing galaxies from the CfA-Arizona Space Telescope LEns Survey by means of stellar population synthesis and pixel-based mass reconstruction. The sample covers a stellar mass range of Ms~2x10^9-3x10^11 Msol (solar masses) which corresponds to a total mass range of ML~7x10^9-3x10^12 Msol on scales from 0.25 to 5 Re (effective radii). By examining the Ms-to-ML dependence on radial distance to the center of each galaxy we find pairs of lenses on small to intermediate mass scales which approach at large radii the same values for their enclosed total mass but exhibit very different stellar masses and stellar baryon fractions. This behaviour subsides for the most massive lenses. All baryon fraction profiles show that the dark matter halo overtakes the stellar content between 1.5 and 2.5 Re. We find evidence for a stellar baryon fraction steadily declining over the full mass range. We shed light on the Fundamental Plane puzzle by showing that the slope of the ML(<R)-to-Ms(<R) relation approaches the mass-to-light relation of recent Fundamental Plane studies at large radii. Less massive dark matter halos turn out to be influenced by the distribution of stellar matter on resolved scales below 10 kpc. The ongoing study of resolved baryon fraction profiles will make it possible to evaluate the validity of star formation models as well as adiabatic contraction prescriptions commonly used in simulations. [abridged]

preprint2010arXiv

An analytic solution for weak-field Schwarzschild geodesics

It is well known that the classical gravitational two body problem can be transformed into a spherical harmonic oscillator by regularization. We find that a modification of the regularization transformation has a similar result to leading order in general relativity. In the resulting harmonic oscillator, the leading-order relativistic perturbation is formally a negative centrifugal force. The net centrifugal force changes sign at three Schwarzschild radii, which interestingly mimics the innermost stable circular orbit (ISCO) of the full Schwarzschild problem. Transforming the harmonic-oscillator solution back to spatial coordinates yields, for both timelike and null weak-field Schwarzschild geodesics, a solution for $t,r,ϕ$ in terms of elementary functions of a variable that can be interpreted as a generalized eccentric anomaly. The textbook expressions for relativistic precession and light deflection are easily recovered. We suggest how this solution could be combined with additional perturbations into numerical methods suitable for applications such as relativistic accretion or dynamics of the Galactic-centre stars.

preprint2010arXiv

Constraining the low-mass end of the Initial Mass Function with Gravitational Lensing

The low-mass end of the stellar Initial Mass Function (IMF) is constrained by focusing on the baryon-dominated central regions of strong lensing galaxies. We study in this letter the Einstein Cross (Q2237+0305), a z=0.04 barred galaxy whose bulge acts as lens on a background quasar. The positions of the four quasar images constrain the surface mass density on the lens plane, whereas the surface brightness (H-band NICMOS/HST imaging) along with deep spectroscopy of the lens (VLT/FORS1) allow us to constrain the stellar mass content, for a range of IMFs. We find that a classical single power law (Salpeter IMF) predicts more stellar mass than the observed lensing estimates. This result is confirmed at the 99% confidence level, and is robust to systematic effects due to the choice of population synthesis models, the presence of dust, or the complex disk/bulge population mix. Our non-parametric methodology is more robust than kinematic estimates, as we do not need to make any assumptions about the dynamical state of the galaxy or its decomposition into bulge and disk. Over a range of low-mass power law slopes (with Salpeter being Gamma=+1.35) we find that at a 90% confidence level, slopes with Gamma>0 are ruled out.

preprint2010arXiv

Relativistic redshift effects and the Galactic-center stars

The high pericenter velocities (up to a few percent of light) of the S stars around the Galactic-center black hole suggest that general relativistic effects may be detectable through the time variation of the redshift during pericenter passage. Previous work has computed post-Newtonian perturbations to the stellar orbits. We study the additional redshift effects due to perturbations of the light path (what one may call "post-Minkowskian'' effects), a calculation that can be elegantly formulated as a boundary-value problem. The post-Newtonian and post-Minkowskian redshift effects are comparable: both are O(beta^3) and amount to a few km/s at pericenter for the star S2. On the other hand, the post-Minkowskian redshift contribution of spin is O(beta^5) and much smaller than the O(beta^4) post-Newtonian effect, which would be approximately 0.1km/s for S2.

preprint2010arXiv

Testing GR with Galactic-centre Stars

The Galactic Centre S-stars orbiting the central supermassive black hole reach velocities of a few percent of the speed of light. The GR-induced perturbations to the redshift enter the dynamics via two distinct channels. The post-Newtonian regime perturbs the orbit from the Keplerian (Zucker et al., 2006, Kannan & Saha 2009), and the photons from the Minkowski (Angelil & Saha 2010). The inclusion of gravitational time dilation at order v^2 marks the first departure of the redshift from the line-of-sight velocities. The leading-order Schwarzschild terms curve space, and enter at order v^3. The classical Keplerian phenomenology dominates the total redshift. Spectral measurements of sufficient resolution will allow for the detection of these post-Newtonian effects. We estimate the spectral resolution required to detect each of these effects by fitting the redshift curve via the five keplerian elements plus black hole mass to mock data. We play with an exaggerated S2 orbit - one with a semi-major axis a fraction of that of the real S2. This amplifies the relativistic effects, and allows clear visual distinctions between the relativistic terms. We argue that spectral data of S2 with a dispersion of about 10km/s would allow for a clear detection of gravitational redshift, and about 1 km/s would suffice for leading-order space curvature detection.

preprint2010arXiv

Towards relativistic orbit fitting of Galactic center stars and pulsars

The S stars orbiting the Galactic center black hole reach speeds of up to a few percent the speed of light during pericenter passage. This makes, for example, S2 at pericenter much more relativistic than known binary pulsars, and opens up new possibilities for testing general relativity. This paper develops a technique for fitting nearly-Keplerian orbits with perturbations from Schwarzschild curvature, frame dragging, and spin-induced torque, to redshift measurements distributed along the orbit but concentrated around pericenter. Both orbital and light-path effects are taken into account. It turns out that absolute calibration of rest-frame frequency is not required. Hence, if pulsars on orbits similar to the S stars are discovered, the technique described here can be applied without change, allowing the much greater accuracies of pulsar timing to be taken advantage of. For example, pulse timing of 3 microsec over one hour amounts to an effective redshift precision of 30 cm/s, enough to measure frame dragging and the quadrupole moment from an S2-like orbit, provided problems like the Newtonian "foreground" due to other masses can be overcome. On the other hand, if stars with orbital periods of order a month are discovered, the same could be accomplished with stellar spectroscopy from the E-ELT at the level of 1 km/s.

preprint2010arXiv

Understanding micro-image configurations in quasar microlensing

The micro-arcsecond scale structure of the seemingly point-like images in lensed quasars, though unobservable, is nevertheless much studied theoretically, because it affects the observable (or macro) brightness, and through that provides clues to substructure in both source and lens. A curious feature is that, while an observable macro-image is made up of a very large number of micro-images, the macro flux is dominated by a few micro-images. Micro minima play a key role, and the well-known broad distribution of macro magnification can be decomposed into narrower distributions with 0,1,2,3,... micro minima. This paper shows how the dominant micro-images exist alongside the others, using the ideas of Fermat's principle and arrival-time surfaces, alongside simulations.

preprint2009arXiv

Interpreting the Kustaanheimo-Stiefel transform in gravitational dynamics

The Kustaanheimo-Stiefel transform turns a gravitational two-body problem into a harmonic oscillator, by going to four dimensions. In addition to the mathematical-physics interest, the KS transform has proved very useful in N-body simulations, where it helps handle close encounters. Yet the formalism remains somewhat arcane, with the role of the extra dimension being especially mysterious. This paper shows how the basic transformation can be interpreted as a rotation in three dimensions. For example, if we slew a telescope from zenith to a chosen star in one rotation, we can think of the rotation axis and angle as the KS transform of the star. The non-uniqueness of the rotation axis encodes the extra dimension. This geometrical interpretation becomes evident on writing KS transforms in quaternion form, which also helps derive concise expressions for regularized equations of motion.

preprint2009arXiv

Weak Microlensing

A nearby star having a near-transit of a galaxy will cause a time-dependent weak lensing of the galaxy. Because the effect is small, we refer to this as weak microlensing. This could provide a useful method to weigh low-mass stars and brown dwarfs. We examine the feasibility of measuring masses in this way and we find that a star causes measurable weak microlensing in a galaxy even at 10 Einstein radii away. Of order one magnitude I < 25 galaxy comes close enough to one or other of the ~100 nearest stars per year.

preprint2008arXiv

Frame-dragging and the kinematics of Galactic-Center stars

We calculate the effects of frame dragging on the Galactic-Center stars. Assuming the stars are only slightly relativistic, we derive an approximation to the Kerr metric, which turns out to be a weak field Schwarzschild metric plus a frame dragging term. By numerically integrating the resulting geodesic equations, we compute the effect on keplerian elements and the kinematics. We find that the kinematic effect at pericenter passage is proportional to (a(1-e^2))^{-2}. For known Galactic-center stars it is of order 10 m/s. If observed this would provide a measurement of the spin of the black hole.

preprint1997arXiv

Non-parametric reconstruction of the galaxy-lens in PG1115+080

We describe a new, non-parametric, method for reconstructing lensing mass distributions in multiple-image systems, and apply it to PG1115, for which time delays have recently been measured. It turns out that the image positions and the ratio of time delays between different pairs of images constrain the mass distribution in a linear fashion. Since observational errors on image positions and time delay ratios are constantly improving, we use these data as a rigid constraint in our modelling. In addition, we require the projected mass distributions to be inversion-symmetric and to have inward-pointing density gradients. With these realistic yet non-restrictive conditions it is very easy to produce mass distributions that fit the data precisely. We then present models, for $H_0=42$, 63 and 84 \kmsmpc, that in each case minimize mass-to-light variations while strictly obeying the lensing constraints. (Only a very rough light distribution is available at present.) All three values of $H_0$ are consistent with the lensing data, but require quite different morphologies for the lensing galaxy. If $H_0$ is low, the main lensing galaxy could be an early type. If $H_0$ is high, the galaxy is reconstructed as a late elliptical or nearly edge-on disc; a binary or merging system is also plausible.