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.

In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated datasets presented in paper I. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20\%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ($ |A| < 2\%$) and precision ($< 6\%$ per system), even in the presence of poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix, and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges.

LDN 1157, is one of the several clouds situated in the cloud complex, LDN 1147/1158, represents a coma-shaped morphology with a well-collimated bipolar outflow emanating from a Class 0 protostar, LDN 1157-mm. The main goals of this work are (a) to map the inter-cloud magnetic field (ICMF) geometry of the region surrounding LDN 1157 to investigate its relationship with the cloud morphology, with the outflow direction and with the core magnetic field (CMF) geometry inferred from the mm- and sub-mm polarization results from the literature, and (b) to investigate the kinematic structure of the cloud. We carried out R-band polarization observations of the stars projected on the cloud to map the pc-scale magnetic field geometry and made spectroscopic observations of the entire cloud in 12CO, C18O and N2H+ (J=1-0) lines to investigate its kinematic structure. We obtained a distance of 340$\pm$3 pc to the LDN 1147/1158, complex based on the Gaia DR2 parallaxes and proper motion values of the three YSOs associated with the complex. A single filament of $\sim1.2$ pc in length and $\sim0.09$ pc in width is found to run all along the coma-shaped cloud. Based on the relationships between the ICMF, CMF, filament orientations, outflow direction, and the presence of an hour-glass morphology of the magnetic field, it is likely that the magnetic field had played an important role in the star formation process in LDN 1157. Combining the proper motions of the YSOs and the radial velocity of LDN 1147/1158 and another complex LDN 1172/1174 which is situated $\sim2$\dgr~east of it, we found that both the complexes are moving collectively toward the Galactic plane. The filamentary morphology of the east-west segment of LDN 1157 may have formed as a result of mass lost by ablation due to the interaction of the moving cloud with the ambient interstellar medium.

One of the important parameters in cosmology is the parameter characterizing the equation of state (EoS) of the sources driving the cosmic expansion. Epochs that are dominated by radiation, matter, or scalar fields, whether they are probed either directly or indirectly, can be characterised by a unique value of this parameter. However, the EoS parameter during reheating---a phase succeeding inflation which is supposed to rapidly defrost our universe---remains to be understood satisfactorily. In order to circumvent the complexity of defining an instantaneous EoS parameter during reheating, an effective parameter $w_\mathrm{eff}$, which is an average of the EoS parameter over the duration of reheating, is usually considered. The value of $w_\mathrm{eff}$ is often chosen arbitrarily to lie in the range $-1/3 \leq w_\mathrm{eff} \leq 1$. In this work, we consider the time evolution of the EoS parameter during reheating and relate it to inflationary potentials $V(\phi)$ that behave as $\phi^p$ around the minimum, a proposal which can be applied to a wide class of inflationary models. We find that, given the index $p$, the effective EoS parameter $w_{\rm eff}$ is determined uniquely. We discuss the corresponding effects on the reheating temperature and its implications.

The region surrounding the well-known reflection nebula, NGC 7023, illuminated by a Herbig Be star, HD 200775, located in the dark cloud L1174 is studied in this work. Based on the distances and proper motion values from Gaia DR2 of 20 previously known young stellar object candidates, we obtained a distance of 335$\pm$11 pc to the cloud complex L1172/1174. Using polarization measurements of the stars projected on the cloud complex, we show additional evidence for the cloud to be at $\sim$ 335 pc distance. Using this distance and proper motion values of the YSO candidates, we searched for additional comoving sources in the vicinity of HD 200775 and found 20 new sources which show low infrared excess emission and are of age $\sim$ 1 Myr. Among these, 10 YSO candidates and 4 newly identified comoving sources are found to show X-ray emission. Three of the four new sources for which we have obtained optical spectra show H$\alpha$ in emission. About 80% of the total sources are found within $\sim$ 1 pc distance from HD 200775. Spatial correlation of some of the YSO candidates with the Herschel dust column density peaks suggests that star formation is still active in the region and may have been triggered by HD 200775.

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.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

Supernova remnants (SNRs) have a variety of overall morphology as well as rich structures over a wide range of scales. Quantitative study of these structures can potentially reveal fluctuations of density and magnetic field originating from the interaction with ambient medium and turbulence in the expanding ejecta. We have used $1.5$GHz (L band) and $5$GHz (C band) VLA data to estimate the angular power spectrum $C_{\ell}$ of the synchrotron emission fluctuations of the Kepler SNR. This is done using the novel, visibility based, Tapered Gridded Estimator of $C_{\ell}$. We have found that, for $\ell = (1.9 - 6.9) \times 10^{4}$, the power spectrum is a broken power law with a break at $\ell = 3.3 \times 10^{4}$, and power law index of $-2.84\pm 0.07$ and $-4.39\pm 0.04$ before and after the break respectively. The slope $-2.84$ is consistent with 2D Kolmogorov turbulence and earlier measurements for the Tycho SNR. We interpret the break to be related to the shell thickness of the SNR ($0.35 $ pc) which approximately matches $\ell = 3.3 \times 10^{4}$ (i.e., $0.48$ pc). However, for $\ell > 6.9 \times 10^{4}$, the estimated $C_{\ell}$ of L band is likely to have dominant contribution from the foregrounds while for C band the power law slope $-3.07\pm 0.02$ is roughly consistent with $3$D Kolmogorov turbulence like that observed at large $\ell$ for Cas A and Crab SNRs.

Distinguishing conformally coupled frames from the tree-level perturbative observables (scalar spectral index $n_{\rm s}$ and tensor-to-scalar ratio $r$) is challenging in cosmology as they are nearly identical. However, since the background evolution in these two frames differs significantly, we can look for potential signatures in the reheating constraints to discriminate these frames. In this work, we study the reheating phase in these frames and find that the difference in the inflationary energy scales in these frames contributes to a significant difference in the reheating e-folding number and hence, different reheating temperature. This difference will eventually lead to a contrasting thermal history in the two frames, which may have a potential observational signature in future observations. This study will open up an avenue for distinguishing various conformally connected otherwise indistinguishable frames and may finally lead us to the correct theory of gravity for our Universe.

We have studied the preheating phase for a class of plateau inflationary model considering the four-legs interaction term $(1/2)g^2\phi^2\chi^2$ between the inflaton $(\phi)$ and reheating field $(\chi)$. We specifically focus on the effects of a parameter $\phi_*$ that controls inflationary dynamics and the shape of the inflaton potential. For $\phi_* < M_p$, the departure of the inflaton potential from the usual power-law behavior $\phi^n$ significantly modifies the microscopic behavior of the preheating dynamics. We analyze and compare the efficiency of production, thermalization and the final equation of state of the system for different models under consideration with $n=2,4,6$ for two different values of $\phi_*$. Most importantly as we increase $n$, or decrease $\phi_*$, the preheating occurs very efficiently with the final equation of state to be that of the radiation, $w=1/3$. Specially for $n=2$, the final equation of state turned out to be $w\simeq 0.2$. However, a complete decay of inflaton could not be achieved with the four-legs interaction for any model under consideration. Therefore, in order to complete the reheating process, we perform the perturbative analysis for the second stage of the reheating phase. Taking the end product of the preheating phase as an initial condition we have solved the homogeneous Boltzmann equations for both the fields supplemented by the constraints coming from the subsequent entropy conservation. In so doing, we can calculate the reheating temperature which is otherwise ill-defined right after the end of preheating. The temperature can be uniquely fixed for a given inflaton decay constant and the CMB temperature. We also compare our results with the conventional reheating constraint analysis and discuss the limit of inflaton decay constant from the field theory perspective.

In microlensing of a Galactic star by a brown dwarf or other compact object, the amplified image really consists of two unresolved images with slightly different light-travel times. The difference (of order a microsecond) is GM/c^3 times a dimensionless factor depending on the total magnification. Since magnification is well-measured in microlensing events, a single time-delay measurement would provide the mass of the lens, without degeneracies. The challenge is to find an observable that varies on sub-microsecond time scales. This paper notes that the narrow-band intensity of the unresolved image pair will show photon bunching (the Hanbury Brown and Twiss effect), and argues that the lensed intensity will have an auto-correlation peak at the lensing time delay. The ultrafast photon-counting technology needed for this type of measurement exists, but the photon numbers required to give sufficient signal-to-noise appear infeasible at present. Preliminary estimates suggest time-delayed photon bunching may be measurable for lensed early-type main-sequence stars at 10 kpc, with the help of 30 m-class telescopes.

We explore the effect of relic gravitational waves on the primordial phase change from the quark-gluon plasma into the low density hadron gas that occurred approximately $10^{-5}$s after the beginning. We show that the gravitational wave, through doing work on the fluid, modulates the local volumes, causing a pressure modulation which either suppresses or boosts the bubble cavitation rate. The boosted rate is significant, implying that the phase transition could have occurred earlier than if this interaction is not considered.

With the growing consensus on simple power law inflation models not being favored by the PLANCK observation, dynamics for the non-standard form of the inflaton potential gain significant interest in the recent past. In this paper, we analyze in great detail classes of phenomenologically motivated inflationary models with non-polynomial potential which are the generalization of the potential introduced in \cite{mhiggs}. After the end of inflation, inflaton field will coherently oscillate around its minimum. Depending upon the initial amplitude of the oscillation and coupling parameters standard parametric resonance phenomena will occur. Therefore, we will study how the inflationary model parameters play an important role in understanding the resonant structure of our model under study. Subsequently, the universe will go through the perturbative reheating phase. However, without any specific model consideration, we further study the constraints on our models based on model independent reheating constraint analysis.

With the growing consensus on simple power law inflation models not being favored by the PLANCK observations, dynamics for the non-standard inflation gain significant interest in the recent past. In this paper, we analyze in detail a class of supergravity inspired phenomenological inflationary models with non-polynomial potential based on Maity (Nucl.Phys. B919 (2017) 560), and compare the model predictions with the currently most favored Starobinsky and its generalized $\alpha$-attractor models in the ($n_s$ , $r$) plane constrained by PLANCK. Importantly for a wide range of parameter space, our model provides successful inflation in the sub-Planckian regime. We also have performed model independent analysis of reheating in terms of the effective equation of state parameter. In particular, we consider two stages of reheating dynamics with generalized inflaton equation of state in the initial and relativistic equation of state in the later phase. Finally, we show how our generalized reheating analysis constrains the inflation models under consideration

With the Planck 2015 result, most of the well known canonical large field inflation models turned out to be strongly disfavored. Axion inflation is one of such models which is becoming marginalized with the increasing precession of CMB data. In this paper, we have shown that with a simple Galileon type modification to the marginally favored axion model calling G-axion, we can turn them into one of the most favored models with its detectable prediction of $r$ and $n_s$ within its PLANCK $1\sigma $ range for a wide range of parameters. Interestingly it is this modification which plays the important role in turning the inflationary predictions to be independent of the explicit value of axion decay constant $f$. However, dynamics after the inflation turned out to have a non-trivial dependence on $f$. For each G-axion model there exists a critical value of $f_c$ such that for $f>f_c$ we have the oscillating phase after inflation and for $f<f_c$ we have non-oscillatory phase. Therefore, we obtained a range of sub-Planckian value of model parameters which give rise to consistent inflation. However for sub-Planckian axion decay constant the inflaton field configuration appeared to be singular after the end of inflation. To reheat the universe we, therefore, employ the instant preheating mechanism at the instant of first zero crossing of the inflaton. To our surprise, the instant preheating mechanism turned out to be inefficient as opposed to usual non-oscillatory quintessence model. For another class of G-axion model with super-Planckian axion decay constant, we performed in detail the reheating constraints analysis considering the latest PLANCK result.

The 21 cm hyperfine transition of the atomic hydrogen (HI) in ground state is a powerful probe of the neutral gas content of the universe. This radio frequency transition has been used routinely for decades to observe, both in emission and absorption, HI in the Galactic interstellar medium as well as in extragalactic sources. In general, however, it is not trivial to derive the physically relevant parameters like temperature, density or column density from these observations. Here, we have considered the issue of column density estimation from the HI 21 cm emission spectrum for sightlines with a non-negligible optical depth and a mix of gas at different temperatures. The complicated radiative transfer and a lack of knowledge about the relative position of gas clouds along the sightline often make it impossible to uniquely separate the components, and hinders reliable estimation of column densities in such cases. Based on the observed correlation between the 21 cm brightness temperature and optical depth, we propose a method to get an unbiased estimate of the HI column density using only the 21 cm emission spectrum. This formalism is further used for a large sample to study the spin temperature of the neutral interstellar medium.

We report modelling follow-up of recently-discovered gravitational-lens candidates in the Canada France Hawaii Telescope Legacy Survey. Lens modelling was done by a small group of specially-interested volunteers from the SpaceWarps citizen-science community who originally found the candidate lenses. Models are categorised according to seven diagnostics indicating (a) the image morphology and how clear or indistinct it is, (b) whether the mass map and synthetic lensed image appear to be plausible, and (c) how the lens-model mass compares with the stellar mass and the abundance-matched halo mass. The lensing masses range from ~10^11 Msun to >10^13 Msun. Preliminary estimates of the stellar masses show a smaller spread in stellar mass (except for two lenses): a factor of a few below or above ~10^11 Msun. Therefore, we expect the stellar-to-total mass fraction to decline sharply as lensing mass increases. The most massive system with a convincing model is J1434+522 (SW05). The two low-mass outliers are J0206-095 (SW19) and J2217+015 (SW42); if these two are indeed lenses, they probe an interesting regime of very low star-formation efficiency. Some improvements to the modelling software (SpaghettiLens), and discussion of strategies regarding scaling to future surveys with more and frequent discoveries, are included.

Pulsars around the Massive Black Hole (MBH) in the Galactic Center (GC) are expected to be revealed by the incoming facilities (e.g., the Square Kilometre Array). Under a full relativistic framework with the pulsar approximated as a test particle, we investigate the constraints on the spinning of the MBH by monitoring the timing of surrounding pulsars. For GC pulsars orbiting closely around the MBH (e.g., $\lesssim1000$AU), we find that full relativistic treatment in modeling accurately their timing signals can be necessary, as the relativistic signals are orders of magnitude larger than the time of arrival measurement accuracies. Although usually there are near-degeneracies among MBH spin parameters, the constraints on the spinning of the MBH are still very tight. By continuously monitoring a normal pulsar in orbits with a period of $\sim2.6$yr and an eccentricity of $0.3-0.9$ under timing precision of $1-5$ms, within $\sim 8$yr the spin magnitude and the orientations of the GC MBH can be constrained with $2\sigma$ error of $10^{-3}-10^{-2}$ and $10^{-1}-10^\circ$, respectively. Even for pulsars in orbits similar to the detected star S2/S0-2 or S0-102, we find that the spinning of the MBH can still be constrained within $4-8$yr, with the most significant constraints provided near the pericenter passage. If the proper motion of the pulsars with astrometric accuracy of $10\mu$as can also be collected along with the timing measurement, then the position, velocity, mass and the distance to the Solar System of the MBH can be constrained about $\sim10\mu$as, $\sim1\mu$as$/$yr, $\sim 10 M_\odot$ and $\sim1$pc, respectively.

We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The new spectroscopic data enable better subtraction of foreground light, and better identification of multiple background images. The inferred distribution of dark matter is consistent with being centered on the galaxies, as expected by LCDM. Each galaxy's dark matter also appears to be symmetric. Whilst we do not find an offset between mass and light (suggestive of self-interacting dark matter) as previously reported, the numerical simulations that have been performed to calibrate Abell 3827 indicate that offsets and asymmetry are still worth looking for in collisions with particular geometries. Meanwhile, ALMA proves exceptionally useful for strong lens image identifications.

Satellite tracking involves sending electromagnetic signals to Earth. Both the orbit of the spacecraft and the electromagnetic signals themselves are affected by the curvature of spacetime. The arrival time of the pulses is compared to the ticks of local clocks to reconstruct the orbital path of the satellite to high accuracy, and to implicitly measure general relativistic effects. In particular, Schwarzschild space curvature (static) and frame-dragging (stationary) due to the planet's spin affect the satellite's orbit. The dominant relativistic effect on the path of the signal photons is Shapiro delay due to static space curvature. We compute these effects for some current and proposed space missions, using a Hamiltonian formulation in four dimensions. For highly eccentric orbits, such as in the Juno mission and in the Cassini Grand Finale, the relativistic effects have a kick-like nature, which could be advantageous for detecting them if their signatures are properly modeled as functions of time. Frame-dragging appears, in principle, measurable by Juno and Cassini, though not by Galileo 5 and 6. Practical measurement would require disentangling frame-dragging from the Newtonian 'foreground' such as the gravitational quadrupole which has an impact on both the spacecraft's orbit and the signal propagation. The foreground problem remains to be solved.

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.

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 ($\mu \lesssim 2.2$, accounting for a dark matter fraction of 20-30%, where $\mu=1.3$ corresponds to a Kroupa-like IMF) is compatible at the $1\sigma$ 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 ($\Gamma \lesssim 1.5$, where $\Gamma$ 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.

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.

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 $\Lambda$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.

We measure the projected density profile, shape and alignment of the stellar and dark matter mass distribution in 11 strong-lens galaxies. We find that the projected dark matter density profile - under the assumption of a Chabrier stellar initial mass function - shows significant variation from galaxy to galaxy. Those with an outermost image beyond $\sim 10$ kpc are very well fit by a projected NFW profile; those with images within 10 kpc appear to be more concentrated than NFW, as expected if their dark haloes contract due to baryonic cooling. We find that over several half-light radii, the dark matter haloes of these lenses are rounder than their stellar mass distributions. While the haloes are never more elliptical than $e_{dm} = 0.2$, their stars can extend to $e_* > 0.2$. Galaxies with high dark matter ellipticity and weak external shear show strong alignment between light and dark; those with strong shear ($\gamma \gtrsim 0.1$) can be highly misaligned. This is reassuring since isolated misaligned galaxies are expected to be unstable. Our results provide a new constraint on galaxy formation models. For a given cosmology, these must explain the origin of both very round dark matter haloes and misaligned strong-lens systems.

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.

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.

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.

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.

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.

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 $\Delta f/f \sim 10^{-16}$ in an elliptic orbit around the Earth would constrain the PPN parameters $|\beta -1|, |\gamma-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.

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.

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 $\beta$ -- 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.

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.

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.

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 $\Lambda 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.

In this paper we focus on a type of inverse problem in which the data is expressed as an unknown function of the sought and unknown model function (or its discretised representation as a model parameter vector). In particular, we deal with situations in which training data is not available. Then we cannot model the unknown functional relationship between data and the unknown model function (or parameter vector) with a Gaussian Process of appropriate dimensionality. A Bayesian method based on state space modelling is advanced instead. Within this framework, the likelihood is expressed in terms of the probability density function ($pdf$) of the state space variable and the sought model parameter vector is embedded within the domain of this $pdf$. As the measurable vector lives only inside an identified sub-volume of the system state space, the $pdf$ of the state space variable is projected onto the space of the measurables, and it is in terms of the projected state space density that the likelihood is written; the final form of the likelihood is achieved after convolution with the distribution of measurement errors. Application motivated vague priors are invoked and the posterior probability density of the model parameter vectors, given the data is computed. Inference is performed by taking posterior samples with adaptive MCMC. The method is illustrated on synthetic as well as real galactic data.

We develop a novel statistical strong lensing approach to probe the cosmological parameters by exploiting multiple redshift image systems behind galaxies or galaxy clusters. The method relies on free-form mass inversion of strong lenses and does not need any additional information other than gravitational lensing. Since in free-form lensing the solution space is a high-dimensional convex polytope, we consider Bayesian model comparison analysis to infer the cosmological parameters. The volume of the solution space is taken as a tracer of the probability of the underlying cosmological assumption. In contrast to parametric mass inversions, our method accounts for the mass-sheet degeneracy, which implies a degeneracy between the steepness of the profile and the cosmological parameters. Parametric models typically break this degeneracy, introducing hidden priors to the analysis that contaminate the inference of the parameters. We test our method with synthetic lenses, showing that it is able to infer the assumed cosmological parameters. Applied to the CLASH clusters, the method might be competitive with other probes.

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.

We present 3D smoothed particle hydrodynamics simulations of the collapse of clumps formed through gravitational instability in the outer part of a protoplanetary disc. The initial conditions are taken directly from a global disc simulation, and a realistic equation of state is used to follow the clumps as they contract over several orders of magnitude in density, approaching the molecular hydrogen dissociation stage. The effects of clump rotation, asymmetries, and radiative cooling are studied. Rotation provides support against fast collapse, but non-axisymmetric modes develop and efficiently transport angular momentum outward, forming a circumplanetary disc. This transport helps the clump reach the dynamical collapse phase, resulting from molecular hydrogen dissociation, on a thousand-year timescale, which is smaller than timescales predicted by some previous spherical 1D collapse models. Extrapolation to the threshold of the runaway hydrogen dissociation indicates that the collapse timescales can be shorter than inward migration timescales, suggesting that clumps could survive tidal disruption and deliver a proto-gas giant to distances of even a few AU from the central star.

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.

We present accurate time delays for the quadruply imaged quasar HE 0435-1223. The delays were measured from 575 independent photometric points obtained in the R-band between January 2004 and March 2010. With seven years of data, we clearly show that quasar image A is affected by strong microlensing variations and that the time delays are best expressed relative to quasar image B. We measured Delta_t(BC) = 7.8+/-0.8 days, Delta_t(BD) = -6.5+/-0.7 days and Delta_t_CD = -14.3+/-0.8 days. We spacially deconvolved HST NICMOS2 F160W images to derive accurate astrometry of the quasar images and to infer the light profile of the lensing galaxy. We combined these images with a stellar population fitting of a deep VLT spectrum of the lensing galaxy to estimate the baryonic fraction, $f_b$, in the Einstein radius. We measured f_b = 0.65+0.13-0.10 if the lensing galaxy has a Salpeter IMF and f_b = 0.45+0.04-0.07 if it has a Kroupa IMF. The spectrum also allowed us to estimate the velocity dispersion of the lensing galaxy, sigma_ap = 222+/-34 km/s. We used f_b and sigma_ap to constrain an analytical model of the lensing galaxy composed of an Hernquist plus generalized NFW profile. We solve the Jeans equations numerically for the model and explored the parameter space under the additional requirement that the model must predict the correct astrometry for the quasar images. Given the current error bars on f_b and sigma_ap, we did not constrain H0 yet with high accuracy, i.e., we found a broad range of models with chi^2 < 1. However, narrowing this range is possible, provided a better velocity dispersion measurement becomes available. In addition, increasing the depth of the current HST imaging data of HE 0435-1223 will allow us to combine our constraints with lens reconstruction techniques that make use of the full Einstein ring that is visible in this object.

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]

The S-Stars in the Galactic-center region are found to be on near-perfect Keplerian orbits around presumably a supermassive black hole, with periods of 15-50 yr. Since these stars reach a few percent of light speed at pericenter, various relativistic effects are expected, and have been discussed in the literature. We argue that an elegant test of the Einstein equivalence principle should be possible with existing instruments, through spectroscopic monitoring of an S-star concentrated during the months around pericenter, supplemented with an already-adequate astrometric determination of the inclination. In essence, the spectrum of an S-star can be considered a heterogeneous ensemble of clocks in a freely-falling frame, which near pericenter is moving at relativistic speeds.

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.

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.

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.

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.

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.

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,\phi$ 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.

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.

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.

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.

Stars near the Galactic center reach a few percent of light speed during pericenter passage, which makes post-Newtonian effects potentially detectable. We formulate the orbit equations in Hamiltonian form such that the $O(v^2/c^2)$ and $O(v^3/c^3)$ post-Newtonian effects of the Kerr metric appear as a simple generalization of the Kepler problem. A related perturbative Hamiltonian applies to photon paths. We then derive a symplectic integrator with adaptive time-steps, for fast and accurate numerical calculation of post-Newtonian effects. Using this integrator, we explore relativistic effects. Taking the star S2 as an example, we find that general relativity would contribute tenths of mas in astrometry and tens of $\rm km s^{-1}$ in kinematics. (For eventual comparison with observations, redshift and time-delay contributions from the gravitational field on light paths will need to be calculated, but we do attempt these in the present paper.) The contribution from stars, gas, and dark matter in the Galactic center region is still poorly constrained observationally, but current models suggest that the resulting Newtonian perturbation on the orbits could plausibly be of the same order as the relativistic effects for stars with semi-major axes $\gtrsim 0.01$ pc (or 250 mas). Nevertheless, the known and distinctive {\it time dependence} of the relativistic perturbations may make it possible to disentangle and extract both effects from observations.

Gravitationally lensed quasars can be used to map the mass distribution in lensing galaxies and to estimate the Hubble constant H0 by measuring the time delays between the quasar images. Here we report the measurement of two independent time delays in the quadruply imaged quasar WFI J2033-4723 (z = 1.66). Our data consist of R-band images obtained with the Swiss 1.2 m EULER telescope located at La Silla and with the 1.3 m SMARTS telescope located at Cerro Tololo. The light curves have 218 independent epochs spanning 3 full years of monitoring between March 2004 and May 2007, with a mean temporal sampling of one observation every 4th day. We measure the time delays using three different techniques, and we obtain Dt(B-A) = 35.5 +- 1.4 days (3.8%) and Dt(B-C) = 62.6 +4.1/-2.3 days (+6.5%/-3.7%), where A is a composite of the close, merging image pair. After correcting for the time delays, we find R-band flux ratios of F_A/F_B = 2.88 +- 0.04, F_A/F_C = 3.38 +- 0.06, and F_A1/F_A2 = 1.37 +- 0.05 with no evidence for microlensing variability over a time scale of three years. However, these flux ratios do not agree with those measured in the quasar emission lines, suggesting that longer term microlensing is present. Our estimate of H0 agrees with the concordance value: non-parametric modeling of the lensing galaxy predicts H0 = 67 +13/-10 km s-1 Mpc-1, while the Single Isothermal Sphere model yields H0 = 63 +7/-3 km s-1 Mpc-1 (68% confidence level). More complex lens models using a composite de Vaucouleurs plus NFW galaxy mass profile show twisting of the mass isocontours in the lensing galaxy, as do the non-parametric models. As all models also require a significant external shear, this suggests that the lens is a member of the group of galaxies seen in field of view of WFI J2033-4723.

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.

ACO 1703 is a cluster recently found to have a variety of strongly lensed objects: there is a quintuply-imaged system at z=0.888 and several other lensed objects from z=2.2 to 3.0 (the cluster itself is at z=0.28). It is not difficult to model the lens, as previous work has already done. However, lens models are generically non-unique. We generate ensembles of models to explore the non-uniqueness. When the full range of source redshifts is included, all models are close to \rho \propto r^{-1} out to 200 kpc. But if the quint is omitted, both shallower and steeper models (e.g., \rho \propto r^{-2}) are possible. The reason is that the redshift contrast between the quint and the other sources gives a good measurement of the enclosed mass at two different radii, thus providing a good estimate of the mass profile in between. This result supports universal profiles and explains why single-model approaches can give conflicting results. The mass map itself is elongated in the NW-SE direction, like the galaxy distribution. An overdensity in both mass and light is also apparent to the SE, which suggests meso-structure.

We present a spatially resolved comparison of the stellar-mass and total-mass surface distributions of nine early-type galaxies. The galaxies are a subset of the Sloan Lens ACS survey (or SLACS; Bolton et al. 2006). The total-mass distributions are obtained by exploring pixelated mass models that reproduce the lensed images. The stellar-mass distributions are derived from population synthesis models fit to the photometry of the lensing galaxies. Uncertainties - mainly model degeneracies - are also computed. Stars can account for all the mass in the inner regions. A Salpeter IMF actually gives too much stellar mass in the inner regions and hence appears ruled out. Dark matter becomes significant by the half-light radius and becomes increasingly dominant at larger radii. The stellar and dark components are closely aligned, but the actual ellipticities are not correlated. Finally, we attempt to intuitively summarize the results by rendering the density, stellar-vs-dark ratio, and uncertainties as false-colour maps.

We present the distribution of luminous and dark matter in a set of strong lensing (early-type) galaxies. By combining two independent techniques - stellar population synthesis and gravitational lensing - we can compare the baryonic and dark matter content in these galaxies within the regions that can be probed using the images of the lensed background source. Two samples were studied, extracted from the CASTLES and SLACS surveys. The former probes a wider range of redshifts and allows us to explore the mass distribution out to ~5Re. The high resolution optical images of the latter (using HST/ACS) are used to show a pixellated map of the ratio between total and baryonic matter. We find dark matter to be absent in the cores of these galaxies, with an increasing contribution at projected radii R>Re. The slopes are roughly compatible with an isothermal slope (better interpreted as an adiabatically contracted NFW profile), but a large scatter in the slope exists among galaxies. There is a trend suggesting most massive galaxies have a higher content of dark matter in the regions probed by this analysis.

We present non-parametric radial mass profiles for ten QSO strong lensing galaxies. Five of the galaxies have profiles close to $\rho(r)\propto r^{-2}$, while the rest are closer to r^{-1}, consistent with an NFW profile. The former are all relatively isolated early-types and dominated by their stellar light. The latter --though the modeling code did not know this-- are either in clusters, or have very high mass-to-light, suggesting dark-matter dominant lenses (one is a actually pair of merging galaxies). The same models give $H_0^{-1} = 15.2_{-1.7}^{+2.5}\Gyr$ ($H_0 = 64_{-9}^{+8} \legacy$), consistent with a previous determination. When tested on simulated lenses taken from a cosmological hydrodynamical simulation, our modeling pipeline recovers both H_0 and $\rho(r)$ within estimated uncertainties. Our result is contrary to some recent claims that lensing time delays imply either a low H_0 or galaxy profiles much steeper than r^{-2}. We diagnose these claims as resulting from an invalid modeling approximation: that small deviations from a power-law profile have a small effect on lensing time-delays. In fact, as we show using using both perturbation theory and numerical computation from a galaxy-formation simulation, a first-order perturbation of an isothermal lens can produce a zeroth-order change in the time delays.

We map substructure in three strong lensing systems having particularly good image data: the galaxy lens MG J0414+053 and the clusters SDSS J1004+411 and ACO 1689. Our method is to first reconstruct the lens as a pixelated mass map and then substract off the symmetric part (in the galaxy case) or a projected NFW (for the cluster lenses). In all three systems we find extended irregular structures, or meso-structures, having of order 10% of the total mass. In J0414+053, the meso-structure suggests a tidal tail connecting the main lens with a nearby galaxy; however this interpretation is tentative. In the clusters the identification of meso-structure is more secure, especially in ACO 1689 where two independent sets of lensed images imply very similar meso-structure. In all three cases the meso-structures are correlated with galaxies but much more extended and massive than the stellar components of single galaxies. Such extended structures cannot plausibly persist in such high-density regions without being mixed; the crossing times are too short. The meso-structures therefore appear to be merging or otherwise dynamically evolving systems.

In the galaxy clusters SSDS J1004+411 and ACO 1689, a large number of multiply-imaged background objects have recently been observed. We use these data to map the projected mass distribution in the inner regions of these clusters. The source redshifts span a large range, which eliminates the degeneracies that plague nearly all lensing work. As a result the mass maps are very well-constrained, despite very weak prior assumptions. ACO 1689 lenses so many objects that we can afford to map it twice using disjoint sets of images, thus verifying our internal error estimates. We then deproject the mass maps (pretending for this purpose that they are spherical) and obtain inner profiles consistent with rho propto r^{-1} and indistinguishable from recent cold dark matter simulations. Assuming that baryons make only a small difference to the profile outside sim 10 kpc, these results support the prediction of a universal inner profile for dark matter structures, independent of any parametrization of that profile.

In gravitational lensing, steeper mass profiles generically produce longer time delays but smaller magnifications, without necessarily changing the image positions or magnification ratios between different images. This is well known. We find in this paper, however, that even if steepness is fixed, time delays can still have significant model dependence, which we attribute to shape modeling degeneracies. This conclusion follows from numerical experiments with models of 35 galaxy lenses. We suggest that varying and twisting ellipticities, features that are explored by pixelated lens models but not so far by parametric models, have an important effect on time delays.

We present a simultaneous analysis of 10 galaxy lenses having time-delay measurements. For each lens we derive a detailed free-form mass map, with uncertainties, and with the additional requirement of a shared value of the Hubble parameter across all the lenses. We test the prior involved in the lens reconstruction against a galaxy-formation simulation. Assuming a concordance cosmology, we obtain 1/H_0 = 13.5 (+2.5/-1.3) Gyr

(Abridged) We present our VLT/FORS1 deep spectroscopic observations of the gravitationally lensed quasar SDSS J0924+0219, as well as archival HST/NICMOS and ACS images of the same object. The two-epoch spectra, obtained in the Multi Object Spectroscopy (MOS) mode, allow for very accurate flux calibration, spatial deconvolution of the data, and provide the redshift of the lensing galaxy z=0.394 +/- 0.001. These spectra, taken 15 days apart, show only slight continuum variations, while the broad emission lines display obvious changes in the red wing of the Mg II line, in the Fe II bands, and in the central part of the C III] line. Even though variations in the line profiles are present, we do not see any significant differences between the continuum and emission line flux ratios of images A and B of the quasar. Spatial deconvolution of the HST images reveals a double Einstein ring. One ring corresponds to the lensed quasar host galaxy at z=1.524 and a second bluer one, is the image either of a star-forming region in the host galaxy, or of another unrelated lower redshift object. We find that a broad range of lens models gives a satisfactory fit to the data. However, they predict very different time delays, making SDSS J0924+0219 an object of particular interest for photometric monitoring. In addition, the lens models reconstructed using exclusively the constraints from the Einstein rings, or using exclusively the astrometry of the quasar images, are not compatible. This suggests that substructures play an important role in SDSS J0924+0219.

Aims: To predict time delays for a sample of gravitationally lensed quasars and to evaluate the accuracy that can be realistically achieved on the value of H_0. Methods: We consider 14 lensed quasars that are candidates for time-delay monitoring and model them in detail using pixelized lens models. For each system, we provide a mass map, arrival-time surface and the distribution of predicted time-delays in a concordance cosmology, assuming H_0^{-1}=14 Gyr (H_0=70 in local units). Based on the predicted time-delays and on the observational circumstances, we rate each lens as `excellent' or `good' or `unpromising' for time-delay monitoring. Finally, we analyze simulated time delays for the 11 lens rated excellent or good, and show that H_0 can be recovered to a precision of 5%. Results: In combination with COSMOGRAIL paper I on the temporal sampling of lensed quasar light curves, the present work will help design monitoring campaigns of lensed quasars.

We consider the problem of tidal disruption of stars in the centre of a galaxy containing a supermassive binary black hole with unequal masses. We assume that over the separation distance between the black holes the gravitational potential is dominated by the more massive black hole. Also, we assume that the number density of stars is concentric with the primary black hole and has a power law cusp.We show that the bulk of stars with a small angular-momentum component normal to the black hole binary orbit can reach a small value of total angular momentum through secular evolution in the gravitational field of the binary, and hence they can be tidally disrupted. This effect is analogous to the so-called Kozai effect (Kozai, 1962, Lidov, 1961,1962) well known in celestial mechanics. We develop an analytical theory of secular evolution of the stellar orbits and calculate the rate of tidal disruption. We confront our analytical theory with a simple numerical model and find very good agreement. Our results show that for primary black-hole mass \sim 10^{6}-10^{7}M_{\odot}, the black- hole mass ratio q > 10^{-2}, cusp size \sim 1pc, the tidal disruption rate can be as large as \sim 10^{-2}-1M_{\odot}/yr. This is at least 10^{2}-10^{4} times larger than estimated for the case of a single supermassive black hole. The duration of the phase of enhanced tidal disruption is determined by the dynamical friction time scale, and it is rather short: \sim 10^{5}yr. The dependence of the tidal disruption rate on the mass ratio, as well as on the size of the cusp, is also discussed.

We implement Lie transform perturbation theory to second order for the planar spin-orbit problem. The perturbation parameter is the asphericity of the body, with the orbital eccentricity entering as an additional parameter. We study first and second order resonances for different values of these parameters. For nearly spherical bodies like Mercury and the Moon first order perturbation theory is adequate, whereas for highly aspherical bodies like Hyperion the spin is mostly chaotic and perturbation theory is of limited use. However, in between, we identify a parameter range where second order perturbation theory is useful and where as yet unidentified objects may be in second order resonances.

For 18 well-observed gravitationally lensed QSOs, we compare new non-parametric mass profiles for the lensing galaxies with stellar-population models derived from published HST photometry. The large volume of parameter space searched -- with respect to the possible star formation histories -- allows us to infer robust estimates and uncertainties for the stellar masses. The most interesting results are: (1) the transition from little or no dark matter in the inner regions (<r_e) to dark matter dominating on the ~5r_e scale (~20 kpc) is clearly seen in massive ellipticals; (2) Such a trend is not seen in lower-mass galaxies, so that the stellar content dominates the mass budget out to ~5r_e; (3) the radial gradient in the dark-matter fraction for these intermediate redshift galaxies agrees with published data on nearby galaxies. This result can help reconcile the discrepancies found in recent estimates of dark matter in elliptical galaxies using different techniques (e.g. Planetary nebulae versus X-ray). The observed trend suggests the stellar component in massive galaxies extends further out in terms of the dark matter scale radius.

We use numerical simulations to test a broad range of plausible observational strategies designed to measure the time delay between the images of gravitationally lensed quasars. Artificial quasar light curves are created along with Monte-Carlo simulations in order to determine the best temporal sampling to adopt when monitoring the photometric variations of systems with time delays between 5 and 120 days, i.e., always shorter than the visibility window across the year. Few and realistic assumptions are necessary on the quasar photometric variations (peak-to-peak amplitude and time-scale of the variations) and on the accuracy of the individual photometric points. The output of the simulations is the (statistical) relative error made on the time delay measurement, as a function of 1- the object visibility over the year, 2- the temporal sampling of the light curves and 3- the time delay. Also investigated is the effect of long term microlensing variations which must be below the 5 % level (either intrinsically or by subtraction) if the goal is to measure time delays with an accuracy of 1-2 %. However, while microlensing increases the random error on the time delay, it does not significantly increase the systematic error, which is always a factor 5 to 10 smaller than the random error. Finally, it is shown that, when the time delay is comparable to the visibility window of the object, a logarithmic sampling can significantly improve the time delay determination. All results are presented in the form of compact plots to be used to optimize the observational strategy of future monitoring programs.

A recently discovered quadruply-imaged QSO, SDSS J1004+4112 (Inada et al. 2003; Oguri et al. 2004) in the core of a $z=0.68$ galaxy cluster has an unprecedented image separation of ~13''. This lens gives us a unique opportunity to study the detailed mass distribution in the central regions of this cluster. We present free-form reconstructions of the lens using recently developed methods. The projected mass within 100 kpc is well-constrained as 5+/-1 x 10^{13} M_solar, consistent with previous simpler models. Unlike previous work, however, we are able to detect structures in the lens associated with cluster galaxies. We estimate the mass associated with these galaxies, and show that they contribute not more than about 10% of the total cluster mass within 100 kpc. Typical galaxy masses, combined with typical luminosities yield a rough estimate of their mass-to-light ratio, which is ~<10, implying that these galaxies consist mostly of stars, and possess little dark matter.

SDSS J1004+4112 is an unprecedented object. It looks much like several quadruple quasars lensed by individual galaxies, only it is ~10 times larger, and the lens is a cluster dominated by dark matter. We present free-form reconstructions of the lens using recently-developed methods. The projected cluster mass profile is consistent with being shallow, r^{-0.3...-0.5}, and can be fit with either an NFW or a flat-cored 3 dimensional mass distribution. However, we cannot rule out projected profiles as steep as r^{-1.3}. The projected mass within 100 kpc is well-constrained as 5+/-1 10^{13}M_sun, consistent with previous simpler models. Unlike previous work, however, we are able to detect structures in the lens associated with cluster galaxies. We estimate the mass associated with these galaxies, and show that they contribute not more than about 10% of the total cluster mass within 100 kpc. Typical galaxy masses, combined with typical luminosities yield a rough estimate of their mass-to-light ratio, which is in the single digits. Finally, we discuss implications for time-delay measurements in this system, and possibilities for a partial Einstein ring.

We introduce and implement two novel ideas for modeling lensed quasars. The first idea is to require different lenses to agree about H_0. This means that some models for one lens can be ruled out by data on a different lens. We explain using two worked examples. One example models 1115+080, 1608+656 (time-delay quads) and 1933+503 (a prospective time-delay system) all together, yielding time-delay predictions for the third lens and a 90%-confidence estimate of 1/H_0=14.6_{-1.7}^{+9.4} Gyr (H_0=67_{-26}^{+9} km/s/Mpc) assuming Omega_M=0.3, Omega_Lambda=0.7. The other example models the time-delay doubles 1520+530, 1600+434, 1830-211, and 2149-275, which gives 1/H_0=14.5_{-1.5}^{+3.3} Gyr (H_0=67_{-13}^{+8} km/s/Mpc). Our second idea is to write the whole modeling software as a highly interactive Java applet, which can be used both for coarse-grained results inside a browser and for fine-grained results on a workstation. Several obstacles come up in trying to implement a numerically-intensive method thus, but we overcome them.

We can think of a lensed quasar as taking the Hubble time, shrinking it by \~10^{-11}, and then presenting the result to us as a time delay; the shrinking factor is of the order of fractional sky-area that the lens occupies. This cute fact is a straightforward consequence of lensing theory, and enables a simple rescaling of time delays. Observed time delays have a 40-fold range, but after rescaling the range reduces to 5-fold. The latter range depends on details of the lens and lensing configuration--for example, quads have systematically shorter rescaled time delays than doubles--and is as expected from a simple model. The hypothesis that observed time-delay lenses all come from a generalized-isothermal family can be ruled out. But there is no indication of drastically different populations either.

B1422+231 is a quadruply-imaged QSO with an exceptionally large lensing contribution from group galaxies other than the main lensing galaxy. We detect diffuse X-rays from the galaxy group in archival Chandra observations; the inferred temperature is consistent with the published velocity dispersion. We then explore the range of possible mass maps that would be consistent with the observed image positions, radio fluxes, and ellipticities. Under plausible but not very restrictive assumptions about the lensing galaxy, predicted time delays involving the faint fourth image are fairly well constrained around 7/h days.

We show that some characteristics of multiply-imaged QSO systems are very model-independent and can be deduced accurately by simply scrutinizing the relative positions of images and galaxy-lens center. These include the time-ordering of the images, the orientation of the lens potential, and the rough morphology of any ring. Other features can differ considerably between specific models; H_0 is an example. Surprisingly, properties inherited from a circularly symmetric lens system are model-dependent, whereas features that arise from the breaking of circular symmetry are model-independent. We first develop these results from some abstract geometrical ideas, then illustrate them for some well-known systems (the quads Q2237+030, H1413+117, HST14113+5211, PG1115+080, MG0414+0534, B1608+656, B1422+231, and RXJ0911+0551, and the ten-image system B1933+507), and finally remark on two systems (B1359+154 and PMN J0134-0931) where the lens properties are more complex. We also introduce a Java applet which produces simple lens systems, and helps further illustrate the concepts.

Massive structures, such as galaxies, act as strong gravitational lenses on background sources. When the background source is a quasar, several lensed images are seen, as magnified or de-magnified versions of the same object. The detailed study of the image configuration and the measurement of ``time-delays'' between the images yield estimates of the Hubble parameter H_0. We describe in a simple way the phenomenon of strong lensing and review recent progress made in the field, including microlensing by stars in the main lensing galaxy.

We explain how an approximation to the rings formed by the host galaxies in lensed QSOs can be inferred from the QSO data alone. A simple ring image can be made from any lens model by a simple piece of computer graphics: just plot a contour map of the arrival-time surface with closely-spaced contours. We go on to explain that rings should be (a) sensitive to time-delay ratios between different pairs of images, but (b) very insensitive to H_0. We illustrate this for the well-known quads 1115+080 and 1608+656.

We estimate the mass of the central black hole in our Galaxy from stellar kinematical data published by Ghez et al. (1998) and Genzel et al. (2000). For this we develop a method, related to Merritt (1993), for non-parametrically reconstructing the mass profile and the stellar distribution function in the central region of the Galaxy from discrete kinematic data, including velocity errors. Models are searched using the Metropolis algorithm. We assume that the mass distribution is spherical and the stellar velocity distribution is isotropic, and devise a test of this assumption. From proper motions we obtain an enclosed mass of $2.0\pm{0.7}\times10^6{\rm M}_{\odot}$ within the inner $0.0044\rm pc$, from radial velocities we obtain a mass of $2.2^{+1.6}_{-1.0}\times10^6{\rm M}_{\odot}$ within 0.046pc and from three-dimensional velocities we obtain $1.8^{+0.4}_{-0.3}\times10^6{\rm M}_{\odot}$ within 0.046pc.

This paper shows that the mass-sheet degeneracy and other degeneracies in lensing have simple geometrical interpretations: they are mostly rescalings of the arrival-time surface. Different degeneracies appear in Local Group lensing and in cosmological lensing, because in the former the absolute magnification is measured but the image structure is not resolved, whereas in the latter the reverse usually applies. The most dangerous of these is a combination we may call the `mass-disk degeneracy' in multiply-imaging galaxy lenses, which may lead to large systematic uncertainties in estimates of cosmological parameters from these systems.

Evidence for a bar at the center of the Milky Way triggered a renewed enthusiasm for dynamical modelling of the Galactic bar-bulge. Our goal is to compare the kinematics of a sample of tracers, planetary nebulae, widely distributed over the bulge with the corresponding kinematics for a range of models of the inner Galaxy. Three of these models are N-body barred systems arising from the instabilities of a stellar disk (Sellwood, Fux and Kalnajs), and one is a Schwarzschild system constructed to represent the 3D distribution of the COBE/DIRBE near-IR light and then evolved as an N-body system for a few dynamical times (Zhao). For the comparison of our data with the models, we use a new technique developed by Saha (1998). The procedure finds the parameters of each model, i.e. the solar galactocentric distance R_o in model units, the orientation angle phi, the velocity scale (in km/s per model unit), and the solar tangential velocity which best fit the data.

Observed time delays between images of a lensed QSO lead to the determination of the Hubble constant by Refsdal's method, provided the mass distribution in the lensing galaxy is reasonably well known. Since the two or four QSO images usually observed are woefully inadequate by themselves to provide a unique reconstruction of the galaxy mass, most previous reconstructions have been limited to simple parameterized models, which may lead to large systematic errors in the derived H_0 by failing to consider enough possibilities for the mass distribution of the lens. We use non-parametric modeling of galaxy lenses to better explore physically plausible but not overly constrained galaxy mass maps, all of which reproduce the lensing observables exactly, and derive the corresponding distribution of H_0's. Blind tests - where one of us simulated galaxy lenses, lensing observables, and a value for H_0, and the other applied our modeling technique to estimate H_0 indicate that our procedure is reliable. For four simulated lensed QSOs the distribution of inferred H_0 have an uncertainty of \simeq 10% at 90% confidence. Application to published observations of the two best constrained time-delay lenses, PG1115+080 and B1608+656, yields H_0=61 +/- 11 km/s/Mpc at 68% confidence and 61 +/- 18 km/s/Mpc at 90% confidence.

The largest source of uncertainty in the determination of H_0 from a multiply-imaged QSO system is the unknown mass distribution in the lensing galaxy. Parametric models severely restrict the shape of the galaxy thereby underestimating errors and biasing the best estimate of H_0. We present a method that explores the whole of the model space allowed by the image observables and a few general properties of galaxies. We describe blind test of the method and then apply it to PG1115+080 and B1608+656 which yield H_0 between 45 and 80 km/s/Mpc at 90% confidence level, with the best estimate of 60 km/s/Mpc.

The lensing information provided by multiple images, arclets, and statistical distortions can all be formulated as linear constraints on the arrival-time surface, and hence on the mass distribution. This reduces cluster lens reconstruction from combined strong and weak lensing to a standard type of inversion problem. Moreover, the mass sheet degeneracy is broken if there are sources at different redshifts.

We present a version of the Fourier Bessel method first introduced by Fisher et al (1994) and Zaroubi et al (1994) with two extensions: (a) we amend the formalism to allow a generic galaxy weight which can be constant rather than the more conventional overweighting of galaxies at high distances, and (b) we correct for the masked zones by extrapolation of Fourier Bessel modes rather than by cloning from the galaxy distribution in neighbouring regions. We test the procedure extensively on N-body simulations and find that it gives generally unbiased results but that the reconstructed velocities tend to be overpredicted in high-density regions. Applying the formalism to the PSZz redshift catalog, we find that beta = 0.7 +/- 0.5 from a comparison of the reconstructed Local Group velocity to the CMB dipole. From an anisotropy test of the velocity field, we find that beta = 1 CDM models models normalized to the current cluster abundance can be excluded with 90% confidence. The density and velocity fields reconstructed agree with the fields found by Branchini et al (1998) in most points. We find a back-infall into the Great Attractor region (Hydra-Centaurus region) but tests suggest that this may be an artifact. We identify all the major clusters in our density field and confirm the existence of some previously identified possible ones.

We determine the values of parameters of an N-body model for the Galaxy developed by Fux via comparison with an unbiased, homogeneous sample of OH/IR stars. Via Monte-Carlo simulation, we find the plausibilities of the best-fitting models, as well as their errors. The parameters that are constrained best by these projected data are the total mass of the model and the viewing angle of the central Bar, although the distribution of the latter has multiple maxima. The best model has a viewing angle of 44 degrees, semi-major axis of 2.5 kpc, a bar mass of 1.7E10 solar masses and a tangential velocity of the local standard of rest of 171 km/s . We argue that the lower values that are commonly found from stellar data for the viewing angle (around 25 degrees) arise when too few coordinates are available, when the longitude range is too narrow or when low latitudes are excluded from the fit. The new constraints on the viewing angle of the galactic Bar from stellar line-of-sight velocities decrease further the ability of the Bar's distribution to account for the observed micro-lensing optical depth toward Baade's window : our model reproduces only half the observed value. The signal of triaxiality diminishes quickly with increasing latitude, fading within approximately one scaleheight. This suggests that Baade's window is not a very appropriate region to sample Bar properties.

We apply a new non-parametric technique to reconstruct, with uncertainties, the projected mass distribution of the inner region of Abell 2218, using combined strong and weak lensing constraints from multiple-image systems and arclets with known redshifts. The reconstructed mass map broadly resembles previous, less detailed, parametric models, but when examined in detail shows several sub-structures, not necessarily associated with light but strongly required by the lensing data. In particular, the highest mass peak is offset by about 30 kpc/h50 from the main light peak, and projected mass-to-light in the directions of different cluster galaxies varies by at least a factor of 10. On comparing with mass estimates from models of the X-ray emitting gas, we find that the X-ray models under-predict the enclosed mass profile by, at least a factor of 2.5; the discrepancy gets worse if we assume that mass traces light to the extent allowed by the lensing constraints.

Peebles (1989) showed that in the gravitational instability picture galaxy orbits can be traced back in time from a knowledge of their current positions, via a variational principle. We modify this variational principle so that galaxy redshifts can be input instead of distances, thereby recovering the distances. As a test problem, we apply the new method to a Local Group model. We infer M = 4 to 8 x 10^{12} M_sun depending on cosmology, implying that the dynamics of the outlying Local Group dwarves are consistent with the timing argument. Some algorithmic issues need to be addressed before the method can be applied to recover nonlinear evolution from large redshift surveys, but there are no more difficulties in principle.

We describe a new non-parametric technique for reconstructing the mass distribution in galaxy clusters with strong lensing, i.e., from multiple images of background galaxies. The observed positions and redshifts of the images are considered as rigid constraints and through the lens (ray-trace) equation they provide us with linear constraint equations. These constraints confine the mass distribution to some allowed region, which is then found by linear programming. Within this allowed region we study in detail the mass distribution with minimum mass-to-light variation; also some others, such as the smoothest mass distribution. The method is applied to the extensively studied cluster Abell 370, which hosts a giant luminous arc and several other multiply imaged background galaxies. Our mass maps are constrained by the observed positions and redshifts (spectroscopic or model-inferred by previous authors) of the giant arc and multiple image systems. The reconstructed maps obtained for \a370 reveal a detailed mass distribution, with substructure quite different from the light distribution. The method predicts the bimodal nature of the cluster and that the projected mass distribution is indeed elongated along the axis defined by the two dominant cD galaxies. But the peaks in the mass distribution appear to be offset from the centres of the cDs. We also present an estimate for the total mass of the central region of the cluster. This is in good agreement with previous mass determinations. The total mass of the central region is M=(2.0-2.7) 10^14 Msun/h50, depending on the solution chosen.

This paper describes a method for quantitatively comparing an N-body model with a sample of discrete kinematic data. The comparison has two stages: (i) finding the optimum scaling and orientation of the model relative to the data; and (ii) calculating a goodness of fit, and hence assessing the plausibility of the model in vew of the data. The method derives from considering the data and model both as samples from some underlying binned distribution function, and applying probability theory arguments. As an example, I consider a published N-body model for the Galactic bulge and disc, and fictitious l,b,v measurements, and recover (with error estimates) the spatial and velocity scales of the model and the orientation of the bar. The fictitious data are actually derived from the model by assuming the mass scale and the solar position, but their size and extent mimics a recent survey of OH/IR stars. The results indicate that mass of the bulge and our viewing angle of the bar are usefully estimable from current surveys.

We propose a new technique to reconstruct non-parametrically the projected mass distribution of galaxy clusters from their gravitational lens effect on background galaxies. The beauty of our technique, is that it combines information from multiple imaging (strong lensing) and shear (weak lensing) as linear constraints on the projected mass distribution. Moreover, our technique overcomes all the drawbacks of non-linear methods. The method is applied to the first cluster-lens A370 and to the cluster-lens A2218 which is exceptionally rich in multiple images and arc(let)s. The reconstructed maps for both cases revealed unprecedented levels of detail.

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.

We describe how long-term solar system orbit integration could be implemented on a parallel computer. The interesting feature of our algorithm is that each processor is assigned not to a planet or a pair of planets but to a time-interval. Thus, the 1st week, 2nd week,..., 1000th week of an orbit are computed concurrently. The problem of matching the input to the (n+1)-st processor with the output of the n-th processor can be solved efficiently by an iterative procedure. Our work is related to the so-called waveform relaxation methods in the computational mathematics literature, but is specialized to the Hamiltonian and nearly integrable nature of solar system orbits. Simulations on serial machines suggest that, for the reasonable accuracy requirement of 1" per century, our preliminary parallel algorithm running on a 1000-processor machine would be about 50 times faster than the fastest available serial algorithm, and we have suggestions for further improvements in speed.

We reanalyze published kinematic and photometric data for the cool star population in the central 10 pc (240 arcsec) of the Galaxy, while (a) isolating the photometric data appropriate to this population, and (b) properly allowing for projection effects. Under the assumptions that the system is spherical and isotropic, we find that M/L_K varies from <=1 outside a radius of 0.8 pc to >2 at 0.35 pc. We suggest that such a varying M/L_K could be due to an increasing concentration of stellar remnants towards the Galactic center. Our derived mass-radius curve confirms the existence of ~3x10^6 M_sol within 0.35 pc of the Galactic Center, and ~1.5x10^6 M_sol with 0.2 pc. However, the latter estimate is subject to the uncertain distribution of cool stars in the region. We also consider the dynamics of the hot star population close to the Galactic center and show that the velocity dispersion of He I stars and the surface brightness distribution of the hot stars are consistent with the mass distribution inferred from cool stars.

We argue that individual images in multiple-image QSOs could easily have substructure at the level of $0.1''$ (i.e., unresolvable even with HST); microlensing within such substructure would cause centroid shifts, observable even from the ground as pseudo proper motions. We present a model of the four-image system 2237+0305 in which Image B shows such ``improper motions" of order of $0.01''$ over a few years.

We describe an algorithm for long-term planetary orbit integrations, including the dominant post-Newtonian effects, that employs individual timesteps for each planet. The algorithm is symplectic and exhibits short-term errors that are $O(\epsilon\Omega^2\tau^2)$ where $\tau$ is the timestep, $\Omega$ is a typical orbital frequency, and $\epsilon\ll1$ is a typical planetary mass in solar units. By a special starting procedure long-term errors over an integration interval $T$ can be reduced to $O(\epsilon^2\Omega^3\tau^2T)$. A sample 0.8 Myr integration of the nine planets illustrates that Pluto can have a timestep more than 100 times Mercury's, without dominating the positional error. Our algorithm is applicable to other $N$-body systems.