BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Liliya L. R. Williams

Liliya L. R. Williams appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-phcond-mat.stat-mechastro-ph.IM

Catalog footprint

What is connected

36works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

36 published item(s)

preprint2025arXiv

Transient star B/R ratio and star formation in $z\gtrsim 1$ lensed galaxies

The extreme magnification from galaxy clusters and microlenses therein allows the detection of individual, luminous stars in lensed galaxies as transient events, and hence provides a valuable window into the high mass stellar population in $z>1$ galaxies. As these bright stars can only be formed at specific ages, the relative abundance of transient events at blue (B) and red (R) optical wavelengths ($B/R$ ratio) can provide insights into the recent star formation history of galaxies that are not well constrained by their spectral energy distributions (SEDs). Here, we forward model the transient detection rates in an idealized mock scenario to find that the $B/R$ ratio of strongly lensed $z>1$ galaxies decreases quickly with increasing age. This ratio has moderate sensitivity to metallicity and comparatively low sensitivity to dust attenuation, with no significant dependency on the stellar initial mass function. Fitting model stellar populations to either the SED or $B/R$ ratio alone of ``Warhol'' arc ($z = 0.94$), we find that neither a simple single starburst nor a more complex star formation can simultaneously reproduce both constraints. We then demonstrate that a best-fit model constrained by both the B/R ratio and SED requires a star-formation rate that has varied quite dramatically over the past $\sim$50 Myr, for which the total stellar mass formed over this time is a factor of 10 (with $2-3σ$ significance) different from the best-fit models to the SED alone. Our work shows that the transient $B/R$ ratio can be used as an additional powerful constraint on the recent star formation history of higher-redshift galaxies in future works that are strongly lensed by galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Hubble constant from eight time-delay galaxy lenses

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Extended lens reconstructions with Grale: exploiting time domain, substructural and weak-lensing information

The information about the mass density of galaxy clusters provided by the gravitational lens effect has inspired many inversion techniques. In this article, updates to the previously introduced method in Grale are described, and explored in a number of examples. The first looks into a different way of incorporating time delay information, not requiring the unknown source position. It is found that this avoids a possible bias that leads to "over-focusing" the images, i.e. providing source position estimates that lie in a considerably smaller region than the true positions. The second is inspired by previous reconstructions of the cluster of galaxies MACS J1149.6+2223, where a multiply-imaged background galaxy contained a supernova, SN Refsdal, of which four additional images were produced by the presence of a smaller cluster galaxy. The inversion for the cluster as a whole, was not able to recover sufficient detail interior to this quad. We show how constraints on such different scales, from the entire cluster to a single member galaxy, can now be used, allowing such small scale substructures to be resolved. Finally, the addition of weak lensing information to this method is investigated. While this clearly helps recover the environment around the strong lensing region, the mass sheet degeneracy may make a full strong and weak inversion difficult, depending on the quality of the ellipticity information at hand. We encounter ring-like structure at the boundary of the two regimes, argued to be the result of combining strong and weak lensing constraints, possibly affected by degeneracies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Free-form Grale lens inversion of galaxy clusters with up to 1000 multiple images

In the near future, ultra deep observations of galaxy clusters with HST or JWST will uncover $300-1000$ lensed multiple images, increasing the current count per cluster by up to an order of magnitude. This will further refine our view of clusters, leading to a more accurate and precise mapping of the total and dark matter distribution in clusters, and enabling a better understanding of background galaxy population and their luminosity functions. However, to effectively use that many images as input to lens inversion will require a re-evaluation of, and possibly upgrades to the existing methods. In this paper we scrutinize the performance of the free-form lens inversion method Grale in the regime of $150-1000$ input images, using synthetic massive galaxy clusters. Our results show that with an increasing number of input images, Grale produces improved reconstructed mass distributions, with the fraction of the lens plane recovered at better than $10\%$ accuracy increasing from $40-50\%$ for $\sim\!\!150$ images to $65\%$ for $\sim\!1000$ images. The reconstructed time delays imply a more precise measurement of $H_0$, with $\lesssim 1\%$ bias. While the fidelity of the reconstruction improves with the increasing number of multiple images used as model constraints, $\sim 150$ to $\sim 1000$, the lens plane rms deteriorates from $\sim 0.11''$ to $\sim 0.28''$. Since lens plane rms is not necessarily the best indicator of the quality of the mass reconstructions, looking for an alternative indicator is warranted.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Power law pseudo phase-space density profiles of dark matter halos: fluke of physics?

It has been known for nearly 20 years that the pseudo phase-space density profile of equilibrium simulated dark matter halos, $ρ(r)/σ^3(r)$, is well described by a power law over 3 decades in radius, even though both the density $ρ(r)$, and the velocity dispersion $σ(r)$ deviate significantly from power laws. The origin of this scale-free behavior is not understood. It could be an inherent property of self-gravitating collisionless systems, or it could be a mere coincidence. To address the question we work with equilibrium halos, and more specifically, the second derivative of the Jeans equation, which, under the assumptions of (i) Einasto density profile, (ii) linear velocity anisotropy - density slope relation, and (iii) $ρ/σ^3\propto r^{-α}$, can be transformed from a differential equation to a cubic algebraic equation. Relations (i)-(iii) are all observed in numerical simulations, and are well parametrized by a total of 4 or 6 model parameters. We do not consider dynamical evolution of halos; instead, taking advantage of the fact that the algebraic Jeans equation for equilibrium halos puts relations (i)-(iii) on the same footing, we study the (approximate) solutions of this equation in the 4 and 6 dimensional spaces. We argue that the distribution of best solutions in these parameter spaces is inconsistent with $ρ/σ^3\propto r^{-α}$ being an fundamental property of gravitational evolution, and conclude that the scale-free nature of this quantity is likely to be a fluke.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The BUFFALO HST Survey

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Two-component mass models of the lensing galaxy in the quadruply imaged supernova iPTF16geu

The first resolved, multiply imaged supernova Type Ia, iPTF16geu, was observed 4 years ago, five decades after such systems were first envisioned. Because of the unique properties of the source, these systems hold a lot of promise for the study of galaxy structure and cosmological parameters. However, this very first example presented modelers with a few puzzles. It was expected that to explain image fluxes a contribution from microlensing by stars would be required, but to accommodate the magnitude of microlensing, the density slope of the elliptical power law lens model had to be quite shallow, $ρ_{2D} \propto r^{-0.7}$. Furthermore, the center of mass had to be displaced from that of observed light by ~0.1 kpc, and the position angle of light distribution was misaligned with that of mass by ~40 degrees. In this paper we present mass models that resolve the first two problems, and suggest a resolution of the third. Motivated by observations of local ellipticals, and some recent analysis of galaxy-scale lenses, our mass models consist of two offset (baryonic) mass components. The resulting mass distributions have a single centroid, but are lopsided, and have isodensity contours that are not purely elliptical and not self-similar with radius. For many of our models the microlensing requirements are modest, and the ring formed by the extended supernova host galaxy resembles the observed one.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Image Flux Ratios of Gravitationally Lensed HS 0810+2554 with High Resolution Infrared Imaging

We report near simultaneous imaging using LMIRCam on the LBTI of the quadruply imaged lensed quasar HS 0810+2554 at wavelengths of 2.16, 3.7 and $4.78~μ$m with a Full Width Half Max (FWHM) spatial resolution of $0^{\prime\prime}\!\!.13$, $0^{\prime\prime}\!\!.12$ and $0^{\prime\prime}\!\!.15$ respectively, comparable to HST optical imaging. In the $\rm{z} = 1.5$ rest frame of the quasar, the observed wavelengths correspond to 0.86, 1.48, and $1.91~μ$m respectively. The two brightest images in the quad, A and B, are clearly resolved from each other with a separation of $0.187^{\prime\prime}$. The flux ratio of these two images (A/B) trends from 1.79 to 1.23 from 2.16 to $4.78~μ$m. The trend in flux ratio is consistent with the $2.16~μ$m flux originating from a small sized accretion disk in the quasar that experiences only microlensing. The excess flux above the contribution from the accretion disk at the two longer wavelengths originates from a larger sized region that experiences no microlensing. A simple model employing multiplicative factors for image B due to stellar microlensing $(m)$ and sub-structure millilensing $(M)$ is presented. The result is tightly constrained to the product $m\times M=1.79$. Given the observational errors, the 60\% probability contour for this product stretches from $m= 2.6$, $M = 0.69$ to $m= 1.79$, $M = 1.0$, where the later is consistent with microlensing only.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Model-independent and model-based local lensing properties of B0128+437 from resolved quasar images

The galaxy-scale gravitational lens B0128+437 generates a quadrupole-image configuration of a background quasar that shows milli-arcsecond-scale subcomponents in the multiple images observed with VLBI. As this multiple-image configuration including the subcomponents has eluded a parametric lens-model characterisation so far, we determine local lens properties at the positions of the multiple images with our model-independent approach. Using PixeLens, we also succeed in setting up a global free-form mass density reconstruction including all subcomponents as constraints. We compare the model-independent local lens properties with those obtained by PixeLens and those obtained by the parametric modelling algorithm Lensmodel. A comparison of all three approaches and a model-free analysis based on the relative polar angles of the multiple images corroborate the hypothesis that elliptically symmetric models are too simplistic to characterise the asymmetric mass density distribution of this lenticular or late-type galaxy. In addition, the model-independent approach efficiently determines local lens properties on the scale of the quasar subcomponents, which are computationally intensive to obtain by free-form model-based approaches. As only 40% of the small-scale subcomponent local lens properties overlap within the 1-$σ$ confidence bounds, mass density gradients on milli-arcsecond scales cannot be excluded. Hence, aiming at a global reconstruction of the deflecting mass density distribution, increasingly detailed observations require flexible free-form models that allow for density fluctuations on milli-arcsecond scale to replace parametric ones, especially for asymmetric lenses or lenses with localised inhomogeneities like B0128.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Lens Models Under the Microscope: Comparison of Hubble Frontier Field Cluster Magnification Maps

Using the power of gravitational lensing magnification by massive galaxy clusters, the Hubble Frontier Fields provide deep views of six patches of the high redshift Universe. The combination of deep Hubble imaging and exceptional lensing strength has revealed the greatest numbers of multiply-imaged galaxies available to constrain models of cluster mass distributions. However, even with $\mathcal{O}(100)$ images per cluster, the uncertainties associated with the reconstructions are not negligible. The goal of this paper is to show the diversity of model magnification predictions. We examine 7 and 9 mass models of Abell 2744 and MACS J0416, respectively, submitted to the Mikulski Archive for Space Telescopes for public distribution in September 2015. The dispersion between model predictions increases from 30% at common low magnifications ($μ\sim2$) to 70% at rare high magnifications ($μ\sim40$). MACS J0416 exhibits smaller dispersions than Abell 2744 for $2<μ<10$. We show that magnification maps based on different lens inversion techniques typically differ from each other by more than their quoted statistical errors. This suggests that some models underestimate the true uncertainties, which are primarily due to various lensing degeneracies. Though the exact mass sheet degeneracy is broken, its generalized counterpart is not broken at least in Abell 2744. Other, local degeneracies are also present in both clusters. Our comparison of models is complementary to the comparison of reconstructions of known synthetic mass distributions. By focusing on observed clusters, we can identify those that are best constrained, and therefore provide the clearest view of the distant Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Testing DARKexp against energy and density distributions of Millennium-II halos

We test the DARKexp model for relaxed, self-gravitating, collisionless systems against equilibrium dark matter halos from the Millennium-II simulation. While limited tests of DARKexp against simulations and observations have been carried out elsewhere, this is the first time the testing is done with a large sample of simulated halos spanning a factor of ~ 50 in mass, and using independent fits to density and energy distributions. We show that DARKexp, a one shape parameter family, provides very good fits to the shapes of density profiles, ρ(r), and differential energy distributions, N(E), of individual simulated halos. The best fit shape parameter $ϕ_{0}$ obtained from the two types of fits are correlated, though with scatter. Our most important conclusions come from ρ(r) and N(E) that have been averaged over many halos. These show that the bulk of the deviations between DARKexp and individual Millennium-II halos come from halo-to-halo fluctuations, likely driven by substructure, and other density perturbations. The average ρ(r) and N(E) are quite smooth and follow DARKexp very closely. The only deviation that remains after averaging is small, and located at most bound energies for N(E) and smallest radii for ρ(r). Since the deviation is confined to 3-4 smoothing lengths, and is larger for low mass halos, it is likely due to numerical resolution effects.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ubiquity of density slope oscillations in the central regions of galaxy and cluster-sized systems

One usually thinks of a radial density profile as having a monotonically changing logarithmic slope, such as in NFW or Einasto profiles. However, in two different classes of commonly used systems, this is often not the case. These classes exhibit non-monotonic changes in their density profile slopes which we call oscillations for short. We analyze these two unrelated classes separately. Class 1 consists of systems that have density oscillations and that are defined through their distribution function $f(E)$, or differential energy distribution $N(E)$, such as isothermal spheres, King profiles, or DARKexp, a theoretically derived model for relaxed collisionless systems. Systems defined through $f(E)$ or $N(E)$ generally have density slope oscillations. Class 1 system oscillations can be found at small, intermediate, or large radii but we focus on a limited set of Class 1 systems that have oscillations in the central regions, usually at $\log(r/r_{-2})\lesssim -2$, where $r_{-2}$ is the largest radius where $d\log(ρ)/d\log(r)=-2$. We show that the shape of their $N(E)$ can roughly predict the amplitude of oscillations. Class 2 systems which are a product of dynamical evolution, consist of observed and simulated galaxies and clusters, and pure dark matter halos. Oscillations in the density profile slope seem pervasive in the central regions of Class 2 systems. We argue that in these systems, slope oscillations are an indication that a system is not fully relaxed. We show that these oscillations can be reproduced by small modifications to $N(E)$ of DARKexp. These affect a small fraction of systems' mass and are confined to $\log(r/r_{-2})\lesssim 0$. The size of these modifications serves as a potential diagnostic for quantifying how far a system is from being relaxed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Illuminating a Dark Lens : A Type Ia Supernova Magnified by the Frontier Fields Galaxy Cluster Abell 2744

SN HFF14Tom is a Type Ia Supernova (SN) discovered at z = 1.3457 +- 0.0001 behind the galaxy cluster Abell 2744 (z = 0.308). In a cosmology-independent analysis, we find that HFF14Tom is 0.77 +- 0.15 magnitudes brighter than unlensed Type Ia SNe at similar redshift, implying a lensing magnification of mu_obs = 2.03 +- 0.29. This observed magnification provides a rare opportunity for a direct empirical test of galaxy cluster lens models. Here we test 17 lens models, 13 of which were generated before the SN magnification was known, qualifying as pure "blind tests". The models are collectively fairly accurate: 8 of the models deliver median magnifications that are consistent with the measured mu to within 1-sigma. However, there is a subtle systematic bias: the significant disagreements all involve models overpredicting the magnification. We evaluate possible causes for this mild bias, and find no single physical or methodological explanation to account for it. We do find that model accuracy can be improved to some extent with stringent quality cuts on multiply-imaged systems, such as requiring that a large fraction have spectroscopic redshifts. In addition to testing model accuracies as we have done here, Type Ia SN magnifications could also be used as inputs for future lens models of Abell 2744 and other clusters, providing valuable constraints in regions where traditional strong- and weak-lensing information is unavailable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Model-free analysis of quadruply imaged gravitationally lensed systems and substructured galaxies

Multiple image gravitational lens systems, and especially quads are invaluable in determining the amount and distribution of mass in galaxies. This is usually done by mass modeling using parametric or free-form methods. An alternative way of extracting information about lens mass distribution is to use lensing degeneracies and invariants. Where applicable, they allow one to make conclusions about whole classes of lenses without model fitting. Here, we use approximate, but observationally useful invariants formed by the three relative polar angles of quad images around the lens center to show that many smooth elliptical+shear lenses can reproduce the same set of quad image angles within observational error. This result allows us to show in a model-free way what the general class of smooth elliptical+shear lenses looks like in the three dimensional (3D) space of image relative angles, and that this distribution does not match that of the observed quads. We conclude that, even though smooth elliptical+shear lenses can reproduce individual quads, they cannot reproduce the quad population. What is likely needed is substructure, with clump masses larger than those responsible for flux ratio anomalies in quads, or luminous or dark nearby perturber galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-universality of dark-matter halos: cusps, cores, and the central potential

Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps with mean values of $d\logρ/d\log r \approx -1$, and some dispersion, which is generally parametrized by the varying index $α$ in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for collisionless self-gravitating systems. We demonstrate that DARKexp - which has only one shape parameter, $ϕ_0$ - fits the dispersion in profile shapes of massive simulated halos as well as observed clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf spheroidal galaxies. We provide approximate analytical relations between DARKexp $ϕ_0$, Einasto $α$, or the central logarithmic slope in the Dehnen-Tremaine analytical $γ$-models. The range in halo parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistical mechanics of collisionless orbits. IV. Distribution of angular momentum

It has been shown in previous work that DARKexp, which is a theoretically derived, maximum entropy, one shape parameter model for isotropic collisionless systems, provides very good fits to simulated and observed dark-matter halos. Specifically, it fits the energy distribution, N(E), and the density profiles, including the central cusp. Here, we extend DARKexp N(E) to include the distribution in angular momentum, L^2, for spherically symmetric systems. First, we argue, based on theoretical, semi-analytical, and simulation results, that while dark-matter halos are relaxed in energy, they are not nearly as relaxed in angular momentum, which precludes using maximum entropy to uniquely derive N(E,L^2). Instead, we require that when integrating N(E,L^2) over squared angular momenta one retrieves the DARKexp N(E). Starting with a general expression for N(E,L^2) we show how the distribution of particles in L^2 is related to the shape of the velocity distribution function, VDF, and velocity anisotropy profile, β(r). We then demonstrate that astrophysically realistic halos, as judged by the VDF shape and β(r), must have linear or convex distributions in L^2, for each separate energy bin. The distribution in energy of the most bound particles must be nearly flat, and become more tilted in favor of radial orbits for less bound particles. These results are consistent with numerical simulations and represent an important step towards deriving the full distribution function for spherically symmetric dark-matter halos.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Collisionless dynamics in Globular Clusters

Since globular clusters (GCs) are old, low-N systems their dynamics is widely believed to be fully dominated by collisional two-body processes, and their surface brightness profiles are fit by King models. However, for many GCs, especially those with HST-resolved central regions, and `extra-tidal' features, King models provide poor fits. We suggest that this is partly because collisionless dynamics is also important and contribute to shaping the cluster properties. We show using time-scale and length-scale arguments that except for the very centers of clusters, collisionless dynamics should be more important than collisional. We then fit 38 GCs analyzed by Noyola and Gebhardt (2006) with (collisional) King and (collisionless) DARKexp models over the full available radial range, and find that the latter provide a better fit to 29 GCs; for six of these the fit is at least ~5x better in term of rms. DARKexp models are theoretically derived maximum entropy equilibrium states of self-gravitating collisionless systems and have already been shown to fit the results of dark matter N-body simulations. (We do not attempt fits with ad hoc fitting functions.)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Entropy Production in Collisionless Systems. II. Arbitrary Phase-Space Occupation Numbers

We present an analysis of two thermodynamic techniques for determining equilibria of self-gravitating systems. One is the Lynden-Bell entropy maximization analysis that introduced violent relaxation. Since we do not use the Stirling approximation which is invalid at small occupation numbers, our systems have finite mass, unlike Lynden-Bell's isothermal spheres. (Instead of Stirling, we utilize a very accurate smooth approximation for $\ln{x!}$.) The second analysis extends entropy production extremization to self-gravitating systems, also without the use of the Stirling approximation. In addition to the Lynden-Bell (LB) statistical family characterized by the exclusion principle in phase-space, and designed to treat collisionless systems, we also apply the two approaches to the Maxwell-Boltzmann (MB) families, which have no exclusion principle and hence represent collisional systems. We implicitly assume that all of the phase-space is equally accessible. We derive entropy production expressions for both families, and give the extremum conditions for entropy production. Surprisingly, our analysis indicates that extremizing entropy production rate results in systems that have maximum entropy, in both LB and MB statistics. In other words, both thermodynamic approaches lead to the same equilibrium structures.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Entropy Production in Collisionless Systems. I. Large Phase-Space Occupation Numbers

Certain thermal non-equilibrium situations, outside of the astrophysical realm, suggest that entropy production extrema, instead of entropy extrema, are related to stationary states. In an effort to better understand the evolution of collisionless self-gravitating systems, we investigate the role of entropy production and develop expressions for the entropy production rate in two particular statistical families that describe self-gravitating systems. From these entropy production descriptions, we derive the requirements for extremizing the entropy production rate in terms of specific forms for the relaxation function in the Boltzmann equation. We discuss some implications of these relaxation functions and point to future work that will further test this novel thermodynamic viewpoint of collisionless relaxation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Surface Brightness and Intrinsic Luminosity of Ellipticals

{Abridged} We show that the surface brightness (SB) profiles of elliptical galaxies can be parametrized using a linear superposition of 2-3 components, described by functions developed in Dhar & Williams as the 2D projections of a 3D Einasto profile. For a sample of 23 ellipticals with -24 < Mv < -15, our multi-component models span a range of up to 10^6 in SB and 10^5 in radius, have a median rms of 0.032 mag arcsec^-2, and are statistically justified at >3σ. Our models indicate that i) the central component is more concentrated than the outer component; and ii) the central component of 'core' galaxies is much more luminous, extended and concentrated than that of 'cuspy' galaxies, with their near exponential central profiles indicating disk-like systems whose existence must be verified spectroscopically. While such central excess components are not necessarily contrary to the notion of a mass deficit in 'core' galaxies, we show that the existence, amount, radial extent and sign of mass deficits disagree substantially in the literature, both for a given galaxy and on an average over a sample. We discuss possible implications and suggest that SMBH binaries are unlikely to be the sole mechanism for producing the large 'cores'. We also deduce conditions under which the 3D light density can be described with a multi-component Einasto model for both cuspy and core galaxies; indicating an universality in the functional form of the 3D density of light in galaxies and dark matter in LCDM N-body haloes. Finally, we show that our result - the outer component of the SB profiles of massive galaxies has 5 < n < 8 - could imply i) a common feature of collisionless systems; and ii) that galaxies with such n for their outer component are dark matter dominated.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Surface mass density of the Einasto family of dark matter haloes: Are they Sersic-like?

Recent advances in N-body simulations of dark matter halos have shown that three-parameter models, in particular the Einasto profile characterized by d ln ρ(r)/d ln r / r with a shape parameter α < 0.3, are able to produce better fits to the 3D spatial density profiles than two-parameter models like the Navarro, Frenk and White (NFW), and Moore et al. profiles. In this paper, we present for the first time an analytically motivated form for the 2D surface mass density of the Einasto family of dark matter haloes, in terms of the 3D spatial density parameters for a wide range of the shape parameter 0.1 < α < 1. Our model describes a projected (2D) Einasto profile remarkably well between 0 and (3 - 5) r_{200}, with errors less than 0.3 per cent for α < 0.3 and less than 2 per cent for α as large as 1. This model (in 2D) can thus be used to fit strong and weak lensing observations of galaxies and clusters whose total spatial (3D) density distributions are believed to be Einasto-like. Further, given the dependence of our model on the 3D parameters, one can reliably estimate structural parameters of the spatial (3D) density from 2D observations. We also consider a Sersic-like parametrization for the above family of projected Einasto profiles and observe that fits with a Sersic profile are sensitive to whether one fits the projected density in linear scale or logarithmic scale and yield widely varying results. Structural parameters of Einasto-like systems, inferred from fits with a Sersic profile, should be used with caution.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Fundamental Surface of Quad Lenses

In a quadruply imaged lens system the angular distribution of images around the lens center is completely described by three relative angles. We show empirically that in the 3D space of these angles, spanning 180 x 180 x 90 degrees, quads from simple two-fold symmetric lenses of arbitrary radial density profile and arbitrary radially dependent ellipticity or external shear define a nearly invariant 2D surface. We give a fitting formula for the surface using SIS+elliptical lensing potential. Various circularly symmetric mass distributions with shear up to 0.4 deviate from it by typically, rms~0.1 deg, while elliptical mass distributions with ellipticity of up 0.4 deviate from it by rms~1.5 deg. The existence of a near invariant surface gives a new insight into the lensing theory and provides a framework for studying quads. It also allows one to gain information about the lens mass distribution from the image positions alone, without any recourse to mass modeling. As an illustration, we show that about 3/4 of observed galaxy-lens quads do not belong to this surface within observational error, and so require additional external shear or substructure to be modeled adequately.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Statistical mechanics of collisionless orbits. I. Origin of central cusps in dark-matter halos

We present an equilibrium statistical mechanical theory of collisionless self-gravitational systems with isotropic velocity distributions. Compared to existing standard theories, we introduce two changes: (1) the number of possible microstates is computed in energy (orbit) space rather than phase space and (2) low occupation numbers are treated more appropriately than using Stirling's approximation. Combined, the two modifications predict that the relaxed parts of collisionless self-gravitating systems, such as dark-matter halos, have a differential energy distribution N(E) ~ [exp(phi_0 - E) - 1], dubbed "DARKexp". Such systems have central power-law density cusps rho(r) ~ r^-1, which suggests a statistical mechanical origin of cusps in simulated dark-matter halos.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Statistical mechanics of collisionless orbits. II. Structure of halos

In this paper, we present the density, ρ, velocity dispersion, σ, and ρ/σ^3 profiles of isotropic systems which have the energy distribution, N(E)\propto[\exp(ϕ_0-E)-1], derived in Paper I. This distribution, dubbed DARKexp, is the most probable final state of a collisionless self-gravitating system, which is relaxed in terms of particle energies, but not necessarily in terms of angular momentum. We compare the DARKexp predictions with the results obtained using the extended secondary infall model (ESIM). The ESIM numerical scheme is optimally suited for the purpose because (1) it relaxes only through energy redistribution, leaving shell/particle angular momenta unaltered, and (2) being a shell code with radially increasing shell thickness it has very good mass resolution in the inner halo, where the various theoretical treatments give different predictions. The ESIM halo properties, and especially their energy distributions, are very well fit by DARKexp, implying that the techniques of statistical mechanics can be used to explain the structure of relaxed self-gravitating systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Statistical mechanics of collisionless orbits. III. Comparison with N-body simulations

We compare the DARKexp differential energy distribution, N(E) \propto \exp(ϕ_0-E)-1, obtained from statistical mechanical considerations, to the results of N-body simulations of dark matter halos. We first demonstrate that if DARKexp halos had anisotropic velocity distributions similar to those of N-body simulated halos, their density and energy distributions could not be distinguished from those of isotropic DARKexp halos. We next carry out the comparison in two ways, using (1) the actual energy distribution extracted from simulations, and (2) N-body fitting formula for the density distribution as well as N(E) computed from the density using the isotropic Eddington formula. Both the methods independently agree that DARKexp N(E) with ϕ_0\approx 4-5 is an excellent match to N-body N(E). Our results suggest (but do not prove) that statistical mechanical principles of maximum entropy can be used to explain the equilibrated final product of N-body simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Understanding micro-image configurations in quasar microlensing

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dark Matter Angular Momentum Profile from the Jeans Equation

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 $λ$ increases as a function of radius.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Radial Orbit Instability in Collisionless N-Body Simulations

Using a suite of self-gravitating, collisionless N-body models, we systematically explore a parameter space relevant to the onset and behavior of the radial orbit instability (ROI), whose strength is measured by the systemic axis ratios of the models. We show that a combination of two initial conditions, namely the velocity anisotropy and the virial ratio, determines whether a system will undergo ROI and exactly how triaxial the system will become. A third initial condition, the radial shape of the density profile, plays a smaller, but noticeable role. Regarding the dynamical development of the ROI, the instability a) begins after systems collapse to their most compact configuration and b) evolves fastest when a majority of the particles have radially anisotropic orbits while there is a lack of centrally-concentrated isotropic orbits. We argue that this is further evidence that self-reinforcing torques are the key to the onset of the ROI. Our findings support the idea that a separate orbit instability plays a role in halting the ROI.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Density Profiles of Collisionless Equilibria. II. Anisotropic Spherical Systems

It has long been realized that dark matter halos formed in cosmological N-body simulations are characterized by density profiles rho(r) that, when suitably scaled, have similar shapes. Additionally, combining the density and velocity dispersion profiles sigma(r), each of which have decidedly nonpower-law shapes, leads to quantity rho/(sigma^3) that is a power-law in radius over 3 orders of magnitude in radius. Halos' velocity anisotropy profiles beta(r) vary from isotropic near the centers of halos to quite radially anisotropic near the virial radius. Finally, there appears to be a nearly linear correlation between beta and the logarithmic density slope gamma for a wide variety of halos. This work is part of a continuing investigation of the above interrelationships and their origins using analytical and semi-analytical techniques. Our finding suggest that the nearly linear beta--gamma relationship is not just another expression of scale-free rho/(sigma^3) behavior. We also note that simultaneously reproducing density and anisotropy profiles like those found in simulations requires beta(r) and gamma(r) to have similar shapes, leading to nearly linear beta--gamma correlations. This work suggests that the beta--gamma and power-law rho/(sigma^3) relations have distinct physical origins.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Angular Distribution of Gamma-ray Bursts and Weak Lensing

We investigate whether Gamma-Ray Bursts (GRBs) from the Current BATSE Catalog have been affected by weak lensing by the nearby large scale structure. The redshift distribution of GRBs is believed to be broad, extending to z~5, so most events can be assumed to be at large redshifts, and hence subject to weak lensing, which would betray itself as projected (anti-)correlations between GRB events and galaxies or clusters that trace the intervening mass. Given the observed distribution of GRBs in fluence, and statistical positional error, e, we predict that most subsets drawn from BATSE Catalog will be anti-correlated with the foreground structure due to weak lensing, i.e. will show negative magnification bias. We find that GRBs are indeed anti-correlated with the APM galaxies (z~0.2-0.3) in the sense that galaxy density in circles of radii 1-1.5 deg (15-20 Mpc at z~0.3) centered on e<1 GRBs is about 10% lower than expected from a random distribution; the significance of GRB-APM anti-correlations reaches 99.7%. Cross-correlation between GRBs and distant rich Abell-Corwin-Olowin clusters is also negative. Standard cosmological models with Omega_matter ~ 0.3, Omega_Lambda ~ 0.7, and matter distribution on large scales following observed APM galaxy distribution with the biasing parameter of around 1 are not able to reproduce our GRB-APM anti-correlations. We propose a speculative model that does account for these anti-correlations as well as positive correlations found previously, between QSOs and APM galaxies. We briefly discuss if the proposed scheme is in conflict with observations of cosmic microwave background, galaxy surveys, cosmic velocity flows, and weak shear lensing.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Image Separation vs. Redshift of Lensed QSOs: Implications for Galaxy Mass Profiles

Recently, Park and Gott reported an interesting observation: image separation of lensed QSOs declines with QSO redshift more precipitously than expected in any realistic world model, if the lenses are taken to be either singular isothermal spheres or point masses. In this Letter I propose that the observed trend arises naturally if the lensing galaxies have logarithmic surface mass density profiles that gradually change with radius. If the observed lack of central (odd) images is also taken into account, the data favor a universal dark matter density profile over an isothermal sphere with a core. Since the trend of image separation vs. source redshift is mostly a reflection of galaxy properties, it cannot be straightforwardly used as a test of cosmological models. Furthermore, the current upper limits on the cosmological constant may have to be revised.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

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

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

0 citations0 reviewsNo code linkedOpen access