Kinematic deprojection and mass inversion of spherical systems of known velocity anisotropy

Traditionally, the mass / velocity anisotropy degeneracy (MAD) inherent in the spherical, stationary, non-streaming Jeans equation has been handled by assuming a mass profile and fitting models to the observed kinematical data. Here, the opposite approach is considered: the equation of anisotropic kinematic projection is inverted for known arbitrary anisotropy to yield the space radial velocity dispersion profile in terms of an integral involving the radial profiles of anisotropy and isotropic dynamical pressure. Then, through the Jeans equation, the mass profile is derived in terms of double integrals of observable quantities. Single integral formulas for both deprojection and mass inversion are provided for several simple anisotropy models (isotropic, radial, circular, general constant, Osipkov-Merritt, Mamon-Lokas and Diemand-Moore-Stadel). Tests of the mass inversion on NFW models with these anisotropy models yield accurate results in the case of perfect observational data, and typically better than 70% (in 4 cases out of 5) accurate mass profiles for the sampling errors expected from current observational data on clusters of galaxies. For the NFW model with mildly increasing radial anisotropy, the mass is found to be insensitive to the adopted anisotropy profile at 7 scale radii and to the adopted anisotropy radius at 3 scale radii. This anisotropic mass inversion method is a useful complementary tool to analyze the mass and anisotropy profiles of spherical systems. It provides the practical means to lift the MAD in quasi-spherical systems such as globular clusters, round dwarf spheroidal and elliptical galaxies, as well as groups and clusters of galaxies, when the anisotropy of the tracer is expected to be linearly related to the slope of its density.Comment: Accepted in MNRAS. 19 pages. Minor changes from previous version: Table 1 of nomenclature, some math simplifications, paragraph in Discussion on alternative deprojection method by deconvolution. 19 pages. 6 figure

Universal inversion formulas for recovering a function from spherical means

The problem of reconstruction a function from spherical means is at the heart of several modern imaging modalities and other applications. In this paper we derive universal back-projection type reconstruction formulas for recovering a function in arbitrary dimension from averages over spheres centered on the boundary an arbitrarily shaped smooth convex domain. Provided that the unknown function is supported inside that domain, the derived formulas recover the unknown function up to an explicitly computed smoothing integral operator. For elliptical domains the integral operator is shown to vanish and hence we establish exact inversion formulas for recovering a function from spherical means centered on the boundary of elliptical domains in arbitrary dimension.Comment: [20 pages, 2 figures] Compared to the previous versions I corrected some typo

Full field inversion in photoacoustic tomography with variable sound speed

Recently, a novel measurement setup has been introduced to photoacoustic tomography, that collects data in the form of projections of the full 3D acoustic pressure distribution at a certain time instant. Existing imaging algorithms for this kind of data assume a constant speed of sound. This assumption is not always met in practice and thus leads to erroneous reconstructions. In this paper, we present a two-step reconstruction method for full field detection photoacoustic tomography that takes variable speed of sound into account. In the first step, by applying the inverse Radon transform, the pressure distribution at the measurement time is reconstructed point-wise from the projection data. In the second step, one solves a final time wave inversion problem where the initial pressure distribution is recovered from the known pressure distribution at the measurement time. For the latter problem, we derive an iterative solution approach, compute the required adjoint operator, and show its uniqueness and stability

Force dipoles and stable local defects on fluid vesicles

An exact description is provided of an almost spherical fluid vesicle with a fixed area and a fixed enclosed volume locally deformed by external normal forces bringing two nearby points on the surface together symmetrically. The conformal invariance of the two-dimensional bending energy is used to identify the distribution of energy as well as the stress established in the vesicle. While these states are local minima of the energy, this energy is degenerate; there is a zero mode in the energy fluctuation spectrum, associated with area and volume preserving conformal transformations, which breaks the symmetry between the two points. The volume constraint fixes the distance $S$, measured along the surface, between the two points; if it is relaxed, a second zero mode appears, reflecting the independence of the energy on $S$; in the absence of this constraint a pathway opens for the membrane to slip out of the defect. Logarithmic curvature singularities in the surface geometry at the points of contact signal the presence of external forces. The magnitude of these forces varies inversely with $S$ and so diverges as the points merge; the corresponding torques vanish in these defects. The geometry behaves near each of the singularities as a biharmonic monopole, in the region between them as a surface of constant mean curvature, and in distant regions as a biharmonic quadrupole. Comparison of the distribution of stress with the quadratic approximation in the height functions points to shortcomings of the latter representation. Radial tension is accompanied by lateral compression, both near the singularities and far away, with a crossover from tension to compression occurring in the region between them.Comment: 26 pages, 10 figure

ADAM: a general method for using various data types in asteroid reconstruction

We introduce ADAM, the All-Data Asteroid Modelling algorithm. ADAM is simple and universal since it handles all disk-resolved data types (adaptive optics or other images, interferometry, and range-Doppler radar data) in a uniform manner via the 2D Fourier transform, enabling fast convergence in model optimization. The resolved data can be combined with disk-integrated data (photometry). In the reconstruction process, the difference between each data type is only a few code lines defining the particular generalized projection from 3D onto a 2D image plane. Occultation timings can be included as sparse silhouettes, and thermal infrared data are efficiently handled with an approximate algorithm that is sufficient in practice due to the dominance of the high-contrast (boundary) pixels over the low-contrast (interior) ones. This is of particular importance to the raw ALMA data that can be directly handled by ADAM without having to construct the standard image. We study the reliability of the inversion by using the independent shape supports of function series and control-point surfaces. When other data are lacking, one can carry out fast nonconvex lightcurve-only inversion, but any shape models resulting from it should only be taken as illustrative global-scale ones.Comment: 11 pages, submitted to A&