We analyze the Miyamoto--Nagai substitution, which was introduced over forty years ago to build models of thick disks and flattened elliptical galaxies. Through it, any spherical potential can be converted to an axisymmetric potential via the replacement of spherical polar $r^2$ with $R^2 + ( a + \!\sqrt{z^2+b^2} )^2$, where ($R,z$) are cylindrical coordinates and $a$ and $b$ are constants. We show that if the spherical potential has everywhere positive density, and satisfies some straightforward constraints, then the transformed model also corresponds to positive density everywhere. This is in sharp contradistinction to substitutions like $r^2 \rightarrow R^2 + z^2/q^2$, which leads to simple potentials but can give negative densities. We use the Miyamoto--Nagai substitution to generate a number of new flattened models with analytic potential--density pairs. These include (i) a flattened model with an asymptotically flat rotation curve, which (unlike Binney's logarithmic model) is always non-negative for a wide-range of axis ratios, (ii) flattened generalizations of the hypervirial models which include Satoh's disk as a limiting case and (iii) a flattened analogue of the Navarro--Frenk--White halo which has the cosmologically interesting density fall-off of (distance)$^{-3}$. Finally, we discuss properties of the prolate and triaxial generalizations of the Miyamoto-Nagai substitution.

Isolated, steady-state galaxies correspond to equilibrium solutions of the Poisson--Vlasov system. We show that (i) all galaxies with a distribution function depending on energy alone must be spherically symmetric and (ii) all axisymmetric galaxies with a distribution function depending on energy and the angular momentum component parallel to the symmetry axis must also be reflection-symmetric about the plane $z=0$. The former result is Lichtenstein's Theorem, derived here by a method exploiting symmetries of solutions of elliptic partial differential equations, while the latter result is new. These results are subsumed into the Symmetry Theorem, which specifies how the symmetries of the distribution function in configuration or velocity space can control the planes of reflection symmetries of the ensuing stellar system.

In this paper, we study the cosmic constraint to $w$CDM model via $118$ strong gravitational lensing systems which are complied from SLACS, BELLS, LSD and SL2S surveys, where the ratio between two angular diameter distances $D^{obs} = D_A(z_l,z_s)/D_A(0,z_s)$ is taken as a cosmic observable. To obtain this ratio, we adopt two strong lensing models: one is the singular isothermal sphere model (SIS), the other one is the power-law density profile (PLP) model. Via the Markov Chain Mote Carlo method, the posterior distribution of the cosmological model parameters space is obtained. The results show that the cosmological model parameters are not sensitive to the parameterized forms of the power-law index $\gamma$. Furthermore, the PLP model gives a relative tighter constraint to the cosmological parameters than that of the SIS model. The predicted value of $\Omega_m=0.31^{+0.44}_{-0.24}$ by SIS model is compatible with that obtained by {\it Planck}2015: $\Omega_{m}=0.313\pm0.013$. However, the value of $\Omega_m=0.15^{+0.13}_{-0.11}$ based on the PLP model is smaller and has $1.25\sigma$ tension with that obtained by {\it Planck}2015 result. This discrepancy maybe come from the systematic errors.

We show that, provided the principal axes of the second velocity moment tensor of a stellar population are generally unequal and are oriented perpendicular to a set of orthogonal surfaces at each point, then those surfaces must be confocal quadric surfaces and the potential must be separable or Stackel. This is true under the mild assumption that the even part of the distribution function is invariant under time reversal $v_i \rightarrow -v_i$ of each velocity component. In particular, if the second velocity moment tensor is everywhere exactly aligned in spherical polar coordinates, then the potential must be of Stackel form (excepting degenerate cases where two or more of the semiaxes of ellipsoid are everywhere the same). The theorem also has consequences for alignment in cylindrical polar coordinates, which is used in the popular Jeans Anisotropic Models (JAM). We analyse data on the radial velocities and proper motions of a sample of $\sim 7400$ stars in the stellar halo of the Milky Way. We provide the distributions of the tilt angles or misalignments from the spherical polar coordinate systems. We show that in this sample the misalignment is always small (usually within $3^\circ$) for radii between 7 and 12 kpc. The velocity anisotropy is very radially biased ($\beta \approx 0.7$), and almost invariant across the volume in our study. Finally, we construct a triaxial stellar halo in a triaxial NFW dark matter halo using a made-to-measure method. Despite the triaxiality of the potential, the velocity ellipsoid of the stellar halo is nearly spherically aligned within $\sim6^\circ$ for large regions of space, particularly outside the scale radius of the stellar halo. We conclude that the second velocity moment ellipsoid can be close to spherically aligned for a much wider class of potentials than the strong constraints that arise from exact alignment might suggest.

We prove that, if the time-independent distribution function $F(v;x)$ of a steady-state stellar system is symmetric under velocity inversion such that $F(-v_1,v_2,v_3;x)=F(v_1,v_2,v_3;x)$ and the same for $v_2$ and $v_3$, where $(v_1,v_2,v_3)$ is the velocity component projected onto an orthogonal frame, then the potential within which the system is in equilibrium must be separable (i.e. the Staeckel potential). Furthermore, we find that the Jeans equations imply that, if all mixed second moments of the velocity vanish, that is, $\langle v_iv_j\rangle=0$ for any $i\ne j$, in some Staeckel coordinate system and the only non-vanishing fourth moments in the same coordinate are those in the form of $\langle v_i^4\rangle$ or $\langle v_i^2v_j^2\rangle$, then the potential must be separable in the same coordinates. Finally we also show that all second and fourth velocity moments of tracers with an odd power to the radial component $v_r$ being zero is a sufficient condition to guarantee the potential to be of the form $\Phi=f(r)+r^{-2}g(\theta,\phi)$.

A family of spherical halo models with flat circular velocity curves is presented. This includes models in which the rotation curve has a finite central value but declines outwards (like the Jaffe model). It includes models in which the rotation curve is rising in the inner parts, but flattens asymptotically (like the Binney model). The family encompasses models with both finite and singular (cuspy) density profiles. The self-consistent distribution function depending on binding energy $E$ and angular momentum $L$ is derived and the kinematical properties of the models discussed. These really describe the properties of the total matter (both luminous and dark). For comparison with observations, it is better to consider tracer populations of stars. These can be used to represent elliptical galaxies or the spheroidal components of spiral galaxies. Accordingly, we study the properties of tracers with power-law or Einasto profiles moving in the doubloon potential. Under the assumption of spherical alignment, we provide a simple way to solve the Jeans equations for the velocity dispersions. This choice of alignment is supported by observations on the stellar halo of the Milky Way. Power-law tracers have prolate spheroidal velocity ellipsoids everywhere. However, this is not the case for Einasto tracers, for which the velocity ellipsoids change from prolate to oblate spheroidal near the pole. Asymptotic forms of the velocity distributions close to the escape speed are also derived, with an eye to application to the high velocity stars in the Milky Way. Power-law tracers have power-law or Maxwellian velocity distributions tails, whereas Einasto tracers have super-exponential cut-offs.

We study microlensing light curves by a triple lens, in particular, by a primary star plus two planets. A four-fold degeneracy is confirmed in the light curves, similar to the close and wide degeneracy found in a double lens. Furthermore, we derive a set of equations for triple-lens in the external shear approximation. By using these external shear equations, we identify two kinds of continuous degeneracies which may confused double and triple lenses, i.e. the continuous external shear degeneracy among triple-lens systems and the double-triple lens degeneracy. These degeneracies are particularly important in high magnification events, and thus some caution needs to be applied when one infers the fraction of stars hosting multiple planets from microlensing. We study the dependence of the degeneracies on the lensing parameters (e.g., source trajectory) and give recipes about how the degeneracies should be explored with real data.

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

An axially symmetric potential psi(R,z)=psi(r,theta) is completely separable if the ratio s:k is constant. Here r*s=d^2(r^2*psi)/dr/d(theta) and k=d^2(psi)/dR/dz. If beta=s/k, then the potential admits an integral of the form of I=(L^2+beta*v_z^2)/2+xi where xi is some function of positions determined by the potential psi. More generally, an axially symmetric potential respects the third axisymmetric integral of motion -- in addition to the classical integrals of the Hamiltonian and the axial component of the angular momentum -- if there exist three real constants a,b,c (not all simultaneously zero, a^2+b^2+c^2>0) such that a*s+b*h+c*k=0 where r*h=d^2(r*psi)/d(sigma)/d(tau) and (sigma,tau) is the parabolic coordinate in the meridional plane such that sigma^2=r+z and tau^2=r-z.

We present a unified parameterization of the fitting functions suitable for density profiles of dark matter haloes or elliptical galaxies. A notable feature is that the classical Einasto profile appears naturally as the continuous limiting case of the cored subfamily amongst the double power-law profiles of Zhao (1996). Based on this, we also argue that there is basically no qualitative difference between halo models well-fitted by the Einasto profile and the standard NFW model. This may even be the case quantitatively unless the resolutions of simulations and the precisions of fittings are sufficiently high to make meaningful distinction possible.

We use distant Blue Horizontal Branch stars with Galactocentric distances 16 < r/kpc < 48 as kinematic tracers of the Milky Way dark halo. We model the tracer density as an oblate, power-law embedded within a spherical power-law potential. Using a distribution function method, we estimate the overall power-law potential and the velocity anisotropy of the halo tracers. We measure the slope of the potential to be gamma ~ 0.4 and the overall mass within 50 kpc is ~ 4 x 10^11 M_sol. The tracer velocity anisotropy is radially biased with beta ~ 0.5, which is in good agreement with local solar neighbourhood studies. Our results provide an accurate outer circular velocity profile for the Milky Way and suggest a relatively high concentration dark matter halo (c_vir ~ 20).

Assuming the separable augmented density, it is always possible to construct a distribution function of a spherical population with any given density and anisotropy. We consider under what conditions the distribution constructed as such is in fact non-negative everywhere in the accessible phase-space. We first generalize known necessary conditions on the augmented density using fractional calculus. The condition on the radius part R(r^2) (whose logarithmic derivative is the anisotropy parameter) is equivalent to the complete monotonicity of R(1/w)/w. The condition on the potential part on the other hand is given by its derivative up to any order not greater than (3/2-beta) being non-negative where beta is the central anisotropy parameter. We also derive a specialized inversion formula for the distribution from the separable augmented density, which leads to sufficient conditions on separable augmented densities for the non-negativity of the distribution. The last generalizes the similar condition derived earlier for the generalized Cuddeford system to arbitrary separable systems.

Forthcoming astrometric missions such as the Gaia satellite will bring to the fore the problem of estimating the enclosed mass from a set of positions, radial velocities and proper motions of tracer stars. Here, we show how to construct the tracer mass estimator when the proper motion data are available in addition to the usual line-of-sight velocity data. Notably, the mass estimators do not require any assumption on the anisotropy, as it is naturally incorporated through the different components of proper motions. In addition, the separate treatment of the proper motions and the line-of-sight velocities is desirable because they are observationally independent and thus the propagation of the combined uncertainties is rather straightforward. The extension to projected data is also sketched, together with a possible application of measuring the masses of Galactic globular clusters.

A static binary star solution for OGLE-2008-BLG-513 has been found. The mass ratio and separation for this solution are q~0.28, s~0.86. This model has a chi^2 comparable to the chi^2 of the model with a planetary mass ratio and orbital motion that was presented in this paper. Because the model with a static, binary star lens is simpler than the model originally presented in this paper, Occam's Razor says it is likely to be the correct solution. Consequently, this paper has been withdrawn from ApJ.

If the augmented density of a spherical anisotropic system is assumed to be multiplicatively separable to functions of the potential and the radius, the radial function, which can be completely specified by the behavior of the anisotropy parameter alone, also fixes the anisotropic ratios of every higher-order velocity moment. It is inferred from this that the non-negativity of the distribution function necessarily limits the allowed behaviors of the radial function. This restriction is translated into the constraints on the behavior of the anisotropy parameter. We find that not all radial variations of the anisotropy parameter satisfy these constraints and thus that there exist anisotropy profiles that cannot be consistent with any separable augmented density.

We report the discovery of a planet with a high planet-to-star mass ratio in the microlensing event MOA-2009-BLG-387, which exhibited pronounced deviations over a 12-day interval, one of the longest for any planetary event. The host is an M dwarf, with a mass in the range 0.07 M_sun < M_host < 0.49M_sun at 90% confidence. The planet-star mass ratio q = 0.0132 +- 0.003 has been measured extremely well, so at the best-estimated host mass, the planet mass is m_p = 2.6 Jupiter masses for the median host mass, M = 0.19 M_sun. The host mass is determined from two "higher order" microlensing parameters. One of these, the angular Einstein radius \theta_E = 0.31 +- 0.03 mas, is very well measured, but the other (the microlens parallax \pi_E, which is due to the Earth's orbital motion) is highly degenate with the orbital motion of the planet. We statistically resolve the degeneracy between Earth and planet orbital effects by imposing priors from a Galactic model that specifies the positions and velocities of lenses and sources and a Kepler model of orbits. The 90% confidence intervals for the distance, semi-major axis, and period of the planet are 3.5 kpc < D_L < 7.9 kpc, 1.1 AU < a < 2.7AU, and 3.8 yr < P < 7.6 yr, respectively.

We consider the problem of estimating the virial mass of a dark halo from the positions and velocities of a tracer population. Although a number of general tools are available, more progress can be made if we are able to specify the functional form of the halo potential (although not its normalization). Here, we consider the particular case of the cosmologically motivated Navarro-Frenk-White (NFW) halo and develop two simple estimators. We demonstrate their effectiveness against numerical simulations and use them to provide new mass estimates of Carina, Fornax, Sculptor, and Sextans dSphs.

This paper presents a set of new conditions on the augmented density of a spherical anisotropic system that is necessary for the underlying two-integral phase-space distribution function to be non-negative. In particular, it is shown that the partial derivatives of the Abel transformations of the augmented density must be non-negative. Applied for the separable augmented densities, this recovers the result of van Hese et al. (2011).

We show how to estimate the enclosed mass from the observed motions of an ensemble of test particles. Traditionally, this problem has been attacked through virial or projected mass estimators. Here, we examine and extend these systematically, and show how to construct an optimal estimator for any given assumption as to the potential. The estimators do not explicitly depend on any properties of the density of the test objects, which is desirable as in practice such information is dominated by selection effects. As particular examples, we also develop estimators tailored for the problem of estimating the mass of the Hernquist or NFW dark matter haloes from the projected positions and velocities of stars.

We present a family of robust tracer mass estimators to compute the enclosed mass of galaxy haloes from samples of discrete positional and kinematical data of tracers, such as halo stars, globular clusters and dwarf satellites. The data may be projected positions, distances, line of sight velocities or proper motions. Forms of the estimator tailored for the Milky Way galaxy and for M31 are given. Monte Carlo simulations are used to quantify the uncertainty as a function of sample size. For the Milky Way, the satellite sample consists of 26 galaxies with line-of-sight velocities. We find that the mass of the Milky Way within 300 kpc is ~ 0.9 x 10^12 solar masses assuming velocity isotropy. However, the mass estimate is sensitive to the anisotropy and could plausibly lie between 0.7 - 3.4 x 10^12 solar masses. Incorporating the proper motions of 6 Milky Way satellites into the dataset, we find ~ 1.4 x 10^12 solar masses. The range here if plausible anisotropies are used is still broader, from 1.2 - 2.7 x 10^12 solar masses. For M31, there are 23 satellite galaxies with measured line-of-sight velocities, but only M33 and IC 10 have proper motions. We use the line of sight velocities and distances of the satellite galaxies to estimate the mass of M31 within 300 kpc as ~ 1.4 x 10^12 solar masses assuming isotropy. There is only a modest dependence on anisotropy, with the mass varying between 1.3 -1.6 x 10^12 solar masses. Given the uncertainties, we conclude that the satellite data by themselves yield no reliable insights into which of the two galaxies is actually the more massive.

It is shown that, if the tracer population is supported by a spherical dark halo with a core or a cusp diverging more slowly than that of a singular isothermal sphere, the logarithmic cusp slope 'g' of the tracers must be given exactly by g=2b where b is their velocity anisotropy parameter at the center unless the same tracers are dynamically cold at the center. If the halo cusp diverges faster than that of the singular isothermal sphere, the velocity dispersion of the tracers must diverge at the center too. In particular, if the logarithmic halo cusp slope is larger than two, the diverging velocity dispersion also traces the behavior of the potential. The implication of our theorem on projected quantities is also discussed. We argue that our theorem should be understood as a warning against interpreting results based on simplifying assumptions such as isotropy and spherical symmetry.

Cosmological simulations of dark matter structures have shown that the equilibrated dark matter structures have a fairly small angular momentum. It appears from these N-body simulations that the radial profile of the angular momentum has an almost universal behavior, even if the different dark matter structures have experienced very different formation and merger histories. We suggest a perturbed Jeans equation, which includes a rotational term. This is done under a reasonable assumed form of the change in the distribution function. By conjecturing that the (new) subdominant rotation term must be proportional to the (old) dominant mass term, we find a clear connection, which is in rather good agreement with the results of recent high resolution simulations. We also present a new connection between the radial profiles of the angular momentum and the velocity anisotropy, which is also in fair agreement with numerical findings. Finally we show how the spin parameter $\lambda$ increases as a function of radius.

A sample of roughly 1,800 halo subdwarf stars with radial velocities and proper motions is assembled, using the repeated multi-band Sloan Digital Sky Survey photometric measurements in Stripe 82. Our sample of halo subdwarfs is extracted via a reduced proper motion diagram and distances are obtained using photometric parallaxes, thus giving full phase space information. The tilt of the velocity ellipsoid with respect to the spherical polar coordinate system is computed and found to be consistent with zero for two of the three tilt angles, and very small for the third. We prove that if the inner halo is in a steady-state and the triaxial velocity ellipsoid is everywhere aligned in spherical polar coordinates, then the potential must be spherically symmetric. The detectable, but very mild, misalignment with spherical polars is consistent with the perturbative effects of the Galactic disk on a spherical dark halo. Banana orbits are generated at the 1:1 resonance (in horizontal and vertical frequency) by the disk. They populate Galactic potentials at the typical radii of our subdwarf sample, along with the much more dominant short-axis tubes. However, on geometric grounds alone, the tilt cannot vanish for the banana orbits and this leads to a slight, but detectable, misalignment. We argue that the tilt of the stellar halo velocity ellipsoid therefore provides a hitherto largely neglected but important line of argument that the Milky Way's dark halo, which dominates the potential, must be nearly spherical.

We consider the problem of determining the structure of the dark halo of nearby dwarf spheroidal galaxies (dSphs) from the spherical Jeans equations. Whether the dark halos are cusped or cored at the centre is an important strategic problem in modern astronomy. The observational data comprise the line-of-sight velocity dispersion of a luminous tracer population. We show that when such data are analysed to find the dark matter density with the spherical Poisson and Jeans equations, then the generic solution is a dark halo density that is cusped like an isothermal. Although milder cusps (like the Navarro-Frenk-White 1/r cusp and even cores are possible, they are not generic. Such solutions exist only if the anisotropy parameter beta and the logarithmic slope of the stellar density gamma satisfy the constraint gamma = 2 x beta at the centre or if the radial velocity dispersion falls to zero at the centre. This surprisingly strong statement is really a consequence of the assumption of spherical symmetry, and the consequent coordinate singularity at the origin. So, for example, a dSph with an exponential light profile can exist in Navarro-Frenk- White halo and have a flat velocity dispersion, but anisotropy in general drives the dark halo solution to an isothermal cusp. The identified cusp or core is therefore a consequence of the assumptions (particularly of spherical symmetry and isotropy), and not the data.

A number of recent surveys for gravitational lenses have found examples of double Einstein rings. Here, we investigate analytically the occurrence of multiple Einstein rings. We prove, under very general assumptions, that at most one Einstein ring can arise from a mass distribution in a single plane lensing a single background source. Two or more Einstein rings can therefore only occur in multi-plane lensing. Surprisingly, we show that it is possible for a single source to produce more than one Einstein ring. If two point masses (or two isothermal spheres) in different planes are aligned with observer and source on the optical axis, we show that there are up to three Einstein rings. We also discuss the image morphologies for these two models if axisymmetry is broken, and give the first instances of magnification invariants in the case of two lens planes.

It is shown that gravitational lensing by point masses arranged in an infinitely extended regular lattice can be studied analytically using the Weierstrass functions. In particular, we draw the critical curves and the caustic networks for the lenses arranged in regular-polygonal -- square, equilateral triangle, regular hexagon -- grids. From this, the mean number of positive parity images as a function of the average optical depth is derived and compared to the case of the infinitely extended field of randomly distributed lenses. We find that the high degree of the symmetry in the lattice arrangement leads to a significant bias towards canceling of the shear caused by the neighboring lenses on a given lens position and lensing behaviour that is qualitatively distinct from the random star field. We also discuss some possible connections to more realistic lensing scenarios.

This paper provides a complete theoretical treatment of the point-mass lens perturbed by constant external shear, often called the Chang-Refsdal lens. We show that simple invariants exist for the products of the (complex) positions of the four images, as well as moment sums of their signed magnifications. The image topographies and equations of the caustics and critical curves are also studied. We derive the fully analytic expressions for precaustics, which are the loci of non-critical points that map to the caustics under the lens mapping. They constitute boundaries of the region in the image domain that maps onto the interior of the caustics. The areas under the critical curves, caustics and precaustics are all evaluated, which enables us to calculate the mean magnification of the source within the caustics. Additionally, the exact analytic expression for the magnification distribution for the source in the triangular caustics is derived, as well as a useful approximate expression. Finally, we find that the Chang-Refsdal lens with the convergence greater than unity can exhibit third-order critical behaviour, if the reduced shear is exactly equal to \sqrt{3}/2, and that the number of images for N-point masses with non-zero constant shear cannot be greater than 5N-1.

For a wide class of self-gravitating systems, we show that if the density is cusped like 1/r^{gamma} near the center, then the limiting value of the anisotropy parameter beta = 1 - <v_T^2>/(2<v_r^2>) at the center may not be greater than (gamma/2). Here, <v_r^2> and <v_T^2> are the radial and tangential velocity second moments. This follows from the non-negativity of the phase space density. We compare this theorem to other proposed relations between the cusp slope and the central anisotropy to clarify their applicabilities and underlying assumptions. The extension of this theorem to tracer populations in an externally imposed potential is also derived. In particular, for stars moving in the vicinity of a central black hole, this reduces to gamma >= beta+(1/2), indicating that an isotropic system in Keplerian potential should be cusped at least as steep as 1/r^{0.5}. Similar limits have been noticed before for specific forms of the distribution function, but here we establish this as a general result.

There is good evidence from N-body simulations that the velocity distribution in the outer parts of halos is radially anisotropic, with the kinetic energy in the radial direction roughly equal to the sum of that in the two tangential directions. We provide a simple algorithm to generate such cosmologically important distribution functions. Introducing r_E(E), the radius of the largest orbit of a particle with energy E, we show how to write down almost trivially a distribution function of the form f(E,L)=g(r_E)/L for any spherical model -- including the NFW profile. We in addition give the generic form of the distribution function for any model with a local density power-law index and anisotropy parameter, and provide limiting forms appropriate for the central parts and envelopes of dark matter halos. From those, we argue that, regardless of the anisotropy, the density fall-off at large radii must evolve to 1/r^4 or steeper ultimately.

A family of cusped potential-density pairs is presented for modelling galaxies and dark haloes. The density profile is cusped like 1/r^(2-p) at small radii. The distribution function is simple and takes the form L^(p-2) E^[(3p+1)/2] where E is the binding energy and L is the angular momentum. The models all possess the remarkable property that the virial theorem holds locally, from which they earn their name as the hypervirial family. Famously, this property was first discovered by Eddington to hold for the Plummer model in 1916. In fact, the seductive properties of the Plummer model extend to the whole hypervirial family, including the members possessing the cosmologically important cusps of 1/r or 1/r^1.5 or 1/r^1.33. The intrinsic and projected properties of the family of models are discussed in some detail.

This paper provides two families of flexible and simple galaxy models. Representatives of many of the families possess the important cosmological cusps, with the density behaving like 1/r or 1/r^1.33 or 1/r^1.5 at small radii. The density falls off between 1/r^3 and 1/r^5 at large radii. We provide analytic and anisotropic distribution functions for all the models. Unlike many existing methods, our algorithm can yield tangentially anisotropic velocity dispersions in the outer parts, and so is useful for modeling populations of satellite galaxies and substructure in host galaxy halos. As an application, we demonstrate the degeneracy between mass and anisotropy for the satellite galaxy population of the Milky Way. This can introduce a factor of ~3 uncertainty in the mass of the Milky Way as inferred from the kinematics of the satellite population.

We introduce a simple family of models for representing the dark matter in galaxies. The potential and phase space distribution function are all elementary, while the density is cusped. The models are all hypervirial, that is, the virial theorem holds locally, as well as globally. As an application to dark matter studies, we compute some of the properties of gamma-ray sources caused by neutralino self-annihilation in dark matter clumps.

The POINT-AGAPE collaboration is carrying out a search for gravitational microlensing toward M31 to reveal galactic dark matter in the form of MACHOs (Massive Astrophysical Compact Halo Objects) in the halos of the Milky Way and M31. A high-threshold analysis of 3 years of data yields 6 bright, short--duration microlensing events, which are confronted to a simulation of the observations and the analysis. The observed signal is much larger than expected from self lensing alone and we conclude, at the 95% confidence level, that at least 20% of the halo mass in the direction of M31 must be in the form of MACHOs if their average mass lies in the range 0.5-1 M$_\odot$. This lower bound drops to 8% for MACHOs with masses $\sim 0.01$ M$_\odot$. In addition, we discuss a likely binary microlensing candidate with caustic crossing. Its location, some 32' away from the centre of M31, supports our conclusion that we are detecting a MACHO signal in the direction of M31.

We report the results of our statistical study of the POINT-AGAPE CNe catalogue. The first Sloan $r'$ and $i'$ analysis of the MMRD relationship in M31 is performed using the resulting POINT-AGAPE CN catalogue. Good fits are produced to the MMRD in the two filters. However, we are unable to verify the $t_{15}$ relationship for either Sloan filter. The subsequent analysis of our automated CN detection pipeline has provided us with the most thorough knowledge of the completeness of a CN survey to-date. In addition, the large field of view of the survey has permitted us to probe the outburst rate well into the galactic disk, unlike previous CCD imaging surveys. Using this analysis we are able to probe the CN distribution of M31 and evaluate the global nova rate. Using models of the galactic surface brightness of M31, we show that the observed CN distribution consists of a separate bulge and disk population. We also show that the M31 bulge CN eruption rate per unit $r'$ flux is more than five times greater than that of the disk. We deduce a global M31 CN rate of $65^{+16}_{-15}$ year$^{-1}$, a value much higher than found by previous surveys. Using the global rate, we derive a M31 bulge rate of $38^{+15}_{-12}$ year$^{-1}$ and a disk rate of $27^{+19}_{-15}$ year$^{-1}$. Given our understanding of the completeness and an analysis of other sources of error, we conclude that the true global nova rate of M31 is at least 50% higher than was previously thought and this has consequent implications for the presumed CN rate in the Milky Way. We deduce a Galactic bulge rate of $14^{+6}_{-5}$ year$^{-1}$, a disk rate of $20^{+14}_{-11}$ year$^{-1}$ and a global Galactic rate of $34^{+15}_{-12}$ year$^{-1}$.

Perhaps as many as 30 parallax microlensing events are known, thanks to the efforts of the MACHO, OGLE, EROS and MOA experiments monitoring the bulge. Using Galactic models, we construct mock catalogues of microlensing light curves towards the bulge, allowing for the uneven sampling and observational error bars of the OGLE-II experiment. The fraction of parallax events with delta chi^2 > 50 in the OGLE-II database is around ~1%, though higher fractions are reported by some other surveys. This is in accord with expectations from standard Galactic models. The fraction of parallax events depends strongly on the Einstein crossing time (t_E), being less than 5% at t_E = 50 days but rising to 50% at t_E > 1 yr. We find that the existence of parallax signatures is essentially controlled by the acceleration of the observer normalised to the projected Einstein radius on the observer plane divided by t_E^2. The properties of the parallax events - time-scales, projected velocities, source and lens locations - in our mock catalogues are analysed. Typically, ~38% of parallax events are caused by a disk star microlensing a bulge source, while \~33% are caused by a disk star microlensing a disk source (of these disk sources, one sixth are at a distance of 5 kpc or less). There is a significant shift in mean time-scale from 32 d for all events to ~130d for our parallax events. There are corresponding shifts for other parameters, such as the lens-source velocity projected onto the observer plane (~1110 km/s for all events versus ~80 km/s for parallax events) and the lens distance (6.7 kpc versus 3.7 kpc). We also assess the performance of parallax mass estimators and investigate whether our mock catalogue can reproduce events with features similar to a number of conjectured `black hole' lens candidates.

The POINT-AGAPE collaboration surveyed M31 with the primary goal of optical detection of microlensing events, yet its data catalogue is also a prime source of lightcurves of variable and transient objects, including classical novae (CNe). A reliable means of identification, combined with a thorough survey of the variable objects in M31, provides an excellent opportunity to locate and study an entire galactic population of CNe. This paper presents a set of 440 neural networks, working in 44 committees, designed specifically to identify fast CNe. The networks are developed using training sets consisting of simulated novae and POINT-AGAPE lightcurves, in a novel variation on K-fold cross-validation. They use the binned, normalised power spectra of the lightcurves as input units. The networks successfully identify 9 of the 13 previously identified M31 CNe within their optimal working range (and 11 out of 13 if the network error bars are taken into account). They provide a catalogue of 19 new candidate fast CNe, of which 4 are strongly favoured.

When the gravitational lensing potential can be approximated by that of circularly symmetric system affected by weak perturbations, it is found that the shape of the resulting (tangential) caustics is entirely specified by the local azimuthal behaviour of the affecting perturbations. This provides a common mathematical groundwork for understanding problems such as the close-wide (d <-> 1/d) separation degeneracy of binary lens microlensing lightcurves and the shear-ellipticity degeneracy of quadruple image lens modelling.

An automated search is carried out for microlensing events using a catalogue of 44554 variable superpixel lightcurves derived from our three-year monitoring program of M31. Each step of our candidate selection is objective and reproducible by a computer. Our search is unrestricted, in the sense that it has no explicit timescale cut. So, it must overcome the awkward problem of distinguishing long-timescale microlensing events from long-period stellar variables. The basis of the selection algorithm is the fitting of the superpixel lightcurves to two different theoretical models, using variable star and blended microlensing templates. Only if microlensing is preferred is an event retained as a possible candidate. Further cuts are made with regard to (i) sampling, (ii) goodness of fit of the peak to a Paczynski curve, (iii) consistency of the microlensing hypothesis with the absence of a resolved source, (iv) achromaticity, (v) position in the colour-magnitude diagram and (vi) signal-to-noise ratio. Our results are reported in terms of first-level candidates, which are the most trustworthy, and second-level candidates, which are possible microlensing but have lower signal-to-noise and are more questionable. The pipeline leaves just 3 first-level candidates, all of which have very short full-width half-maximum timescale (<5 days) and 3 second-level candidates, which have timescales of 31, 36 and 51 days respectively. We also show 16 third-level lightcurves, as an illustration of the events that just fail the threshold for designation as microlensing candidates. They are almost certainly mainly variable stars. Two of the 3 first-level candidates correspond to known events (PA 00-S3 and PA 00-S4) already reported by the POINT-AGAPE project. The remaining first-level candidate is new.

With their excellent photometric precision and dramatic increase in monitoring frequency, future microlensing survey experiments are expected to be sensitive to very short time-scale, isolated events caused by free-floating and wide-separation planets with mass as low as a few lunar masses. We estimate the probability of measuring the Einstein radius \theta_E for bound and free-floating planets. We carry out detailed simulations of the planetary events expected in next-generation surveys and estimate the resulting uncertainty in \theta_E for these events. We show that, for main-sequence sources and Jupiter-mass planets, the caustic structure of wide-separation planets with projected separations of < 20 AU substantially increases the probability of measuring the dimensionless source size and thus determining \theta_E compared to the case of unbound planets. In this limit where the source is much smaller than the caustic, the effective cross-section to measure \theta_E to 10% is ~25% larger than the full width of the caustic. Measurement of the lens parallax is possible for low-mass planetary events by combined observations from the ground and a satellite located in an L2 orbit; this would complete the mass measurements for such wide-separation planets. Finally, short-duration events caused by bound planets can be routinely distinguished from those caused by free-floating planets for planet-star separations < 20 AU from either the deviations due to the planetary caustic or (more often) the low-amplitude bump from the magnification due to the parent star.

The POINT-AGAPE survey is an optical search for gravitational microlensing events towards the Andromeda Galaxy (M31). As well as microlensing, the survey is sensitive to many different classes of variable stars and transients. Here we describe the automated detection and selection pipeline used to identify M31 classical novae (CNe) and we present the resulting catalogue of 20 CN candidates observed over three seasons. CNe are observed both in the bulge region as well as over a wide area of the M31 disk. Nine of the CNe are caught during the final rise phase and all are well sampled in at least two colours. The excellent light-curve coverage has allowed us to detect and classify CNe over a wide range of speed class, from very fast to very slow. Among the light-curves is a moderately fast CN exhibiting entry into a deep transition minimum, followed by its final decline. We have also observed in detail a very slow CN which faded by only 0.01 mag day$^{-1}$ over a 150 day period. We detect other interesting variable objects, including one of the longest period and most luminous Mira variables. The CN catalogue constitutes a uniquely well-sampled and objectively-selected data set with which to study the statistical properties of classical novae in M31, such as the global nova rate, the reliability of novae as standard-candle distance indicators and the dependence of the nova population on stellar environment. The findings of this statistical study will be reported in a follow-up paper.

The POINT-AGAPE collaboration has been monitoring M31 for three seasons with the Wide Field Camera on the Isaac Newton Telescope. In each season, data are taken for one hour per night for roughly sixty nights during the six months that M31 is visible. The two fields of view straddle the central bulge, northwards and southwards. We have calculated the locations, periods and amplitudes of 35414 variable stars in M31 as a by-product of our microlensing search. The variables are classified according to their period and amplitude of variation. They are classified into population I and II Cepheids, Miras and semi-regular long-period variables. The population I Cepheids are associated with the spiral arms, while the central concentration of the Miras and long-period variables varies noticeably, the stars with brighter (and shorter) variations being much more centrally concentrated. A crucial role in the microlensing experiment is played by the asymmetry signal. It was initially assumed that the variable stars would not be a serious problem as their distributions would be symmetric. We demonstrate that this assumption is not correct. We find that differential extinction associated with the dust lanes causes the variable star distributions to be asymmetric. The size and direction of the asymmetry of the variable stars is measured as a function of period and amplitude of variation. The implications of this discovery for the successful completion of the microlensing experiments towards M31 are discussed. (Abridged)

The lightcurve of PA-99-N2, one of the recently announced microlensing candidates towards M31, shows small deviations from the standard Paczynski form. We explore a number of possible explanations, including correlations with the seeing, the parallax effect and a binary lens. We find that the observations are consistent with an unresolved RGB or AGB star in M31 being microlensed by a binary lens. We find that the best fit binary lens mass ratio is about one hundredth, which is one of most extreme values found for a binary lens so far. If both the source and lens lie in the M31 disk, then the standard M31 model predicts the probable mass range of the system to be 0.02-3.6 solar masses (95 % confidence limit). In this scenario, the mass of the secondary component is therefore likely to be below the hydrogen-burning limit. On the other hand, if a compact halo object in M31 is lensing a disk or spheroid source, then the total lens mass is likely to lie between 0.09-32 solar masses, which is consistent with the primary being a stellar remnant and the secondary a low mass star or brown dwarf. The optical depth (or alternatively the differential rate) along the line of sight toward the event indicates that a halo lens is more likely than a stellar lens provided that dark compact objects comprise no less than 15 per cent (or 5 per cent) of haloes.

We present a new 0.9" resolution 3.29 micron narrowband image of the reflection nebula NGC 7023. We find that the 3.29 micron IEF in NGC 7023 is brightest in narrow filaments NW of the illuminating star. These filaments have been seen in images of K', molecular hydrogen emission lines, the 6.2 and 11.3 micron IEFs, and HCO+. We also detect 3.29 micron emission faintly but distinctly between the filaments and the star. The 3.29 micron image is in contrast to narrowband images at 2.09, 2.14, and 2.18 micron, which show an extended emission peak midway between the filaments and the star, and much fainter emission near the filaments. The [2.18]-[3.29] color shows a wide variation, ranging from 3.4-3.6 mag at the 2 micron continuum peak to 5.5 mag in the filaments. We observe [2.18]-[3.29] to increase smoothly with increasing distance from the star, up until the filament, suggesting that the main difference between the spatial distributions of the 2 micron continuum and the the 3.29 micron emission is related to the incident stellar flux. Our result suggests that the 3.29 micron IEF carriers are likely to be distinct from, but related to, the 2 micron continuum emitters. Our finding also imply that, in NGC 7023, the 2 micron continuum emitters are mainly associated with HI, while the 3.29 micron IEF carriers are primarily found in warm molecular hydrogen, but that both can survive in HI or molecular hydrogen. (abridged)

The POINT-AGAPE collaboration is currently searching for massive compact halo objects (MACHOs) towards the Andromeda galaxy (M31). The survey aims to exploit the high inclination of the M31 disk, which causes an asymmetry in the spatial distribution of M31 MACHOs. Here, we investigate the effects of halo velocity anisotropy and flattening on the asymmetry signal using simple halo models. For a spherically symmetric and isotropic halo, we find that the underlying pixel-lensing rate in far-disk M31 MACHOs is more than 5 times the rate of near-disk events. We find that the asymmetry is increased further by about 30% if the MACHOs occupy radial orbits rather than tangential orbits, but is substantially reduced if the MACHOs lie in a flattened halo. However, even for haloes with a minor-to-major axis ratio q = 0.3, the numbers of M31 MACHOs in the far-side outnumber those in the near-side by a factor of ~2. We show that, if positional information is exploited in addition to number counts, then the number of candidate events required to confirm asymmetry for a range of flattened and anisotropic halo models is achievable, even with significant contamination by variable stars and foreground microlensing events. For pixel-lensing surveys which probe a representative portion of the M31 disk, a sample of around 50 candidates is likely to be sufficient to detect asymmetry within spherical haloes, even if half the sample is contaminated, or to detect asymmetry in haloes as flat as q = 0.3 provided less than a third of the sample comprises contaminants. We also argue that, provided its mass-to-light ratio is less than 100, the recently observed stellar stream around M31 is not problematic for the detection of asymmetry. (Abstract slightly abridged.)

We obtain high-precision limb-darkening measurements in five bands (V, V_E, I_E, I, and H) for the K3 III (Teff=4200 K, [Fe/H]=+0.3, log(g)=2.3) source of the Galactic bulge microlensing event EROS BLG-2000-5. These measurements are inconsistent with the predictions of atmospheric models at >10 sigma. While the disagreement is present in all bands, it is most apparent in I, I_E and V_E, in part because the data are better and in part because the intrinsic disagreement is stronger. We find that when limb-darkening profiles are normalized to have unit total flux, the I-band models for a broad range of temperatures all cross each other at a common point. The solar profile also passes through this point. However, the profile as measured by microlensing does not. We conjecture that the models have incorporated some aspect of solar physics that is not shared by giant atmospheres.

We analyze PLANET photometric observations of the caustic-crossing binary-lens microlensing event, EROS BLG-2000-5, and find that modeling the observed light curve requires incorporation of the microlens parallax and the binary orbital motion. The projected Einstein radius (~r_E = 3.61 +/- 0.11 AU) is derived from the measurement of the microlens parallax, and we are also able to infer the angular Einstein radius (theta_E = 1.38 +/- 0.12 mas) from the finite source effect on the light curve, combined with an estimate of the angular size of the source given by the source position in a color-magnitude diagram. The lens mass, M = 0.612 +/- 0.057 M_sun, is found by combining these two quantities. This is the first time that parallax effects are detected for a caustic-crossing event and also the first time that the lens mass degeneracy has been completely broken through photometric monitoring alone. The combination of the projected Einstein radius and the angular Einstein radius also allows us to conclude that the lens lies in the near side of the disk, within 2.6 kpc of the Sun, while the radial velocity measurement indicates that the source is a Galactic bulge giant.

We analyze PLANET and MACHO observations of MACHO 99-BLG-47, the first nearly-normal microlensing event for which high signal-to-noise-ratio data reveal a well-covered, short-duration anomaly. This anomaly occurs near the the peak of the event. Short-duration anomalies near the peak of otherwise normal events are expected to arise both from extreme-separation (either very close or very wide), roughly equal-mass binary lenses, and from planetary systems. We show that the lens of MACHO 99-BLG-47 is in fact an extreme-separation binary, not a planetary system, thus demonstrating for the first time that these two important classes of events can be distinguished in practice. However, we find that the wide-binary and close-binary lens solutions fit the data equally well, and cannot be distinguished even at Delta-chi^2=1. This degeneracy is qualitatively much more severe than the one identified for MACHO 98-SMC-1 because the present degeneracy spans two rather than one dimension in the magnification field and does not require significantly different blending fractions. In the appendix, we explore this result, and show that it is related to the symmetry in the lens equation.

In this paper, we investigate the effect of a wide binary companion of the lens on the astrometric behavior of Galactic gravitational microlensing events and compare it to the effect on the photometric behavior. We find that the wide binary companion of the lens can affect the centroid motion of images substantially even if the corresponding light curve appears to be the one of a standard single point-mass lens event. The relatively significant effect of the wide binary lens on the astrometric lensing behavior, on one side, calls for careful consideration of the lens binarity in analyzing the future astrometric lensing data. On the other side, larger astrometric effect of the companion makes astrometric lensing an efficient method to detect binary lenses over a broad range of separations.

We show that one can measure the effects of microlens parallax for binary microlensing events with three well-measured peaks -- two caustic crossings plus a cusp approach, and hence derive the projected Einstein radius. Since the angular Einstein radius is measurable from finite source effects for almost any well-observed caustic crossing, triple-peak events can yield the determination of the lens mass. We note that, to a certain extent, rotation of the binary can mimic the effects of parallax, but it should often be possible to disentangle parallax from rotation by making use of the late-time light curve.

Blended light is an important source of degeneracy in the characterization of microlensing events, particularly in binary-lens and high magnification events. We show how the techniques of image subtraction can be applied to form an image of the blend with the source removed. In many cases, it should be possible to construct images with very high signal-to-noise ratio. Analysis of these images can help distinguish between competing models that have different blend fractions, and in some cases should allow direct detection of the lens.

We analyze five years of PLANET photometry of microlensing events toward the Galactic bulge to search for the short-duration deviations from single lens light curves that are indicative of the presence of planetary companions to the primary microlenses. Using strict event selection criteria, we construct a well defined sample of 43 intensively monitored events. We search for planetary perturbations in these events over a densely sampled region of parameter space spanning two decades in mass ratio and projected separation, but find no viable planetary candidates. By combining the detection efficiencies of the events, we find that, at 95% confidence, less than 25% of our primary lenses have companions with mass ratio q=0.01 and separations in the lensing zone, 0.6-1.6 Einstein ring radii. Using a model of the mass, velocity and spatial distribution of bulge lenses, we infer that the majority of our lenses are likely M dwarfs in the Galactic bulge. We conclude that less than 33% of M-dwarfs in the Galactic bulge have Jupiter-mass companions between 1.5 and 4 AU, and less than 45% have 3 Jupiter-mass companions between 1 and 7 AU, the first significant limits on planetary companions to M-dwarfs. We consider the effects of the finite size of the source stars and changing our detection criterion, but find that these do not alter our conclusions substantially.

We present near-infrared J and H surface photometry of 24 of the nearby Seyfert 1.8, 1.9 and 2 galaxies from the CfA Seyfert sample. The excellent angular resolution of HST probes spatial scales as small as tens of parsecs in most of these AGN and is sensitive to the presence of nuclear bars and other potential signatures of the AGN fueling process that channels host galaxy gas and dust to the nuclear region. We have used elliptical isophote techniques to search for nuclear bars in all of these galaxies and have employed a two-dimensional fitting technique to model the nuclear point source and surface brightness distribution of a bright subsample of these galaxies in an attempt to alleviate the impact of the nuclear point source on our sensitivity to nuclear bars. We find stellar nuclear bar candidates in four of these galaxies: Mrk 471, Mrk 270, Mrk 573, and NGC 5929, nearly 20% of the total sample. The percentage rises to ~30% when systems with disturbed morphologies or high inclinations are excluded. The nuclear bars in Mrk 573 and Mrk 270 exhibit some evidence for dust lanes along their leading edges, analogous to the structures seen in host galaxy bars, while the dust lanes in Mrk 471 and NGC 5929 exhibit a more complex morphology. The fact that most of these AGN do not appear to contain stellar nuclear bars suggests that they are not the fueling mechanism for most low-luminosity AGN.

In this paper, we show that the pattern of microlensing light curve anomalies induced by multiple planets are well described by the superposition of those of the single-planet systems where the individual planet-primary binary pairs act as independent lens systems. Since the outer deviation regions around the planetary caustics of the individual planets occur in general at different locations, we find that the pattern of anomalies in these regions are hardly affected by the existence of other planet(s). This implies that even if an event is caused by a multiple planetary system, a simple single-planet lensing model is good enough for the description of most anomalies caused by the source passage of the outer deviation regions. Detection of the anomalies resulting from the source trajectory passing both the outer deviation regions caused by more than two planets will provide a new channel of detecting multiple planets.

We search for signatures of planets in 43 intensively monitored microlensing events that were observed between 1995 and 1999. Planets would be expected to cause a short duration (~1 day) deviation on the smooth, symmetric light curve produced by a single-lens. We find no such anomalies and infer that less than 1/3 of the ~0.3 M_sun stars that typically comprise the lens population have Jupiter-mass companions with semi-major axes in the range of 1.5 AU <a < 4 AU. Since orbital periods of planets at these radii are 3-15 years, the outer portion of this region is currently difficult to probe with any other technique.

We present VLT FORS1 spectroscopy that temporally resolves the second caustic crossing of the Bulge K giant source of microlensing event EROS 2000-BLG-5, the first time this has been accomplished for several phases of a caustic transit. The ~1 angstrom H-alpha equivalent width of the source star increases slightly as the center of the star egresses the caustic and then plummets by 30% during the final limb crossing. These changes are not seen in contemporaneous spectra of control stars in the FORS1 slit, but are qualitatively consistent with expectations from stellar atmosphere models as the caustic differentially magnifies different portions of the stellar face of the target. Observations such as these in a variety of stellar lines are equivalent to atmospheric tomography and are expected to provide a direct test of stellar models.

We present PLANET observations of OGLE-1999-BUL-23, a binary-lens microlensing event towards the Galactic bulge. PLANET observations in the I and V bands cover the event from just before the first caustic crossing until the end of the event. In particular, a densely-sampled second caustic crossing enables us to derive the linear limb-darkening coefficients of the source star; c_V = 0.786^{+0.080}_{-0.078} and c_I = 0.632^{+0.047}_{-0.037}. Combined analysis of the light curve and the color-magnitude diagram suggests that the source star is a G/K subgiant in the Galactic bulge (T_eff = 4800 K). The resulting linear limb-darkening coefficient of the source is consistent with theoretical predictions, although it is likely that non-linearity of the stellar surface brightness profile complicates the interpretation, especially for the I band. The global light curve fit to the data indicates that the event is due to a binary lens of a mass ratio q = 0.39 and a projected separation d = 2.42. The lens/source relative proper motion is (22.8 +_ 1.5) km/s/kpc, typical of bulge/bulge or bulge/disk events.

Microlensing is the only technique likely, within the next 5 years, to constrain the frequency of Jupiter-analogs. The PLANET collaboration has monitored nearly 100 microlensing events of which more than 20 have sensitivity to the perturbations that would be caused by a Jovian-mass companion to the primary lens. No clear signatures of such planets have been detected. These null results indicate that Jupiter mass planets with separations of 1.5-3 AU occur in less than 1/3 of systems. A similar limit applies to planets of 3 Jupiter masses between 1-4 AU.