Normalized to: Price, M.
[1]
oai:arXiv.org:2004.07855 [pdf] - 2079363
Spherical Bayesian mass-mapping with uncertainties: full sky
observations on the celestial sphere
Submitted: 2020-04-16
To date weak gravitational lensing surveys have typically been restricted to
small fields of view, such that the $\textit{flat-sky approximation}$ has been
sufficiently satisfied. However, with Stage IV surveys ($\textit{e.g. LSST}$
and $\textit{Euclid}$) imminent, extending mass-mapping techniques to the
sphere is a fundamental necessity. As such, we extend the sparse hierarchical
Bayesian mass-mapping formalism presented in previous work to the spherical
sky. For the first time, this allows us to construct $\textit{maximum a
posteriori}$ spherical weak lensing dark-matter mass-maps, with principled
Bayesian uncertainties, without imposing or assuming Gaussianty. We solve the
spherical mass-mapping inverse problem in the analysis setting adopting a
sparsity promoting Laplace-type wavelet prior, though this theoretical
framework supports all log-concave posteriors. Our spherical mass-mapping
formalism facilitates principled statistical interpretation of reconstructions.
We apply our framework to convergence reconstruction on high resolution N-body
simulations with pseudo-Euclid masking, polluted with a variety of realistic
noise levels, and show a dramatic increase in reconstruction fidelity compared
to standard approaches. Furthermore we perform the largest joint reconstruction
to date of all publically available shear observational datasets (combining
DESY1, KiDS450 and CFHTLens) and find that our formalism recovers a convergence
map with significantly enhanced small-scale detail. Within our Bayesian
framework we validate, in a statistically rigourous manner, the community's
intuition regarding the need to smooth spherical Kaiser-Squies estimates to
provide physically meaningful convergence maps. Such approaches cannot reveal
the small-scale physical structures that we recover within our framework.
[2]
oai:arXiv.org:1812.04017 [pdf] - 2030436
Sparse Bayesian mass-mapping with uncertainties: local credible
intervals
Submitted: 2018-12-10
Until recently mass-mapping techniques for weak gravitational lensing
convergence reconstruction have lacked a principled statistical framework upon
which to quantify reconstruction uncertainties, without making strong
assumptions of Gaussianity. In previous work we presented a sparse hierarchical
Bayesian formalism for convergence reconstruction that addresses this
shortcoming. Here, we draw on the concept of local credible intervals (cf.
Bayesian error bars) as an extension of the uncertainty quantification
techniques previously detailed. These uncertainty quantification techniques are
benchmarked against those recovered via Px-MALA - a state of the art proximal
Markov Chain Monte Carlo (MCMC) algorithm. We find that typically our recovered
uncertainties are everywhere conservative, of similar magnitude and highly
correlated (Pearson correlation coefficient $\geq 0.85$) with those recovered
via Px-MALA. Moreover, we demonstrate an increase in computational efficiency
of $\mathcal{O}(10^6)$ when using our sparse Bayesian approach over MCMC
techniques. This computational saving is critical for the application of
Bayesian uncertainty quantification to large-scale stage IV surveys such as
LSST and Euclid.
[3]
oai:arXiv.org:1812.04014 [pdf] - 1795991
Sparse Bayesian mass-mapping with uncertainties: hypothesis testing of
structure
Submitted: 2018-12-10
A crucial aspect of mass-mapping, via weak lensing, is quantification of the
uncertainty introduced during the reconstruction process. Properly accounting
for these errors has been largely ignored to date. We present results from a
new method that reconstructs maximum a posteriori (MAP) convergence maps by
formulating an unconstrained Bayesian inference problem with Laplace-type
$\ell_1$-norm sparsity-promoting priors, which we solve via convex
optimization. Approaching mass-mapping in this manner allows us to exploit
recent developments in probability concentration theory to infer theoretically
conservative uncertainties for our MAP reconstructions, without relying on
assumptions of Gaussianity. For the first time these methods allow us to
perform hypothesis testing of structure, from which it is possible to
distinguish between physical objects and artifacts of the reconstruction. Here
we present this new formalism, demonstrate the method on illustrative examples,
before applying the developed formalism to two observational datasets of the
Abel-520 cluster. In our Bayesian framework it is found that neither Abel-520
dataset can conclusively determine the physicality of individual local massive
substructure at significant confidence. However, in both cases the recovered
MAP estimators are consistent with both sets of data.
[4]
oai:arXiv.org:1812.04018 [pdf] - 1966620
Sparse Bayesian mass-mapping with uncertainties: peak statistics and
feature locations
Submitted: 2018-12-10
Weak lensing convergence maps - upon which higher order statistics can be
calculated - can be recovered from observations of the shear field by solving
the lensing inverse problem. For typical surveys this inverse problem is
ill-posed (often seriously) leading to substantial uncertainty on the recovered
convergence maps. In this paper we propose novel methods for quantifying the
Bayesian uncertainty in the location of recovered features and the uncertainty
in the cumulative peak statistic - the peak count as a function of signal to
noise ratio (SNR). We adopt the sparse hierarchical Bayesian mass-mapping
framework developed in previous work, which provides robust reconstructions and
principled statistical interpretation of reconstructed convergence maps without
the need to assume or impose Gaussianity. We demonstrate our uncertainty
quantification techniques on both Bolshoi N-body (cluster scale) and Buzzard
V-1.6 (large scale structure) N-body simulations. For the first time, this
methodology allows one to recover approximate Bayesian upper and lower limits
on the cumulative peak statistic at well defined confidence levels.
[5]
oai:arXiv.org:1612.05060 [pdf] - 1533451
Survival of the impactor during hypervelocity collisions II: An analogue
for high porosity targets
Submitted: 2016-12-15
We investigated how a target's porosity affects the outcome of a collision,
with respect to the impactor's fate. Laboratory impact experiments using
peridot projectiles were performed at a speed range between 0.3 and 3.0 km/s,
onto high porosity water-ice (40%) and fine-grained calcium carbonate (70%)
targets. We report that the amount of implanted material in the target body
increases with increasing target's porosity, while the size frequency
distribution of the projectile's ejecta fragments becomes steeper. A
supplementary Raman study showed no sign of change of the Raman spectra of the
recovered olivine projectile fragments indicate minimal physical change.
[6]
oai:arXiv.org:1512.03262 [pdf] - 1347789
Survival of the impactor during hypervelocity collisions I: An analogue
for low porosity targets
Submitted: 2015-12-10
Recent observations of asteroidal surfaces indicate the presence of materials
that do not match the bulk lithology of the body. A possible explanation for
the presence of these exogenous materials is that they are products of
inter-asteroid impacts in the Main Belt, and thus interest has increased in
understanding the fate of the projectile during hypervelocity impacts. In order
to gain insight into the fate of impactor we have carried out a laboratory
programme, covering the velocity range of 0.38 - 3.50 km/s, devoted to
measuring the survivability, fragmentation and final state of the impactor.
Forsterite olivine and synthetic basalt projectiles were fired onto low
porosity (<10%) pure water-ice targets using the University of Kent's Light Gas
Gun (LGG). We developed a novel method to identify impactor fragments which
were found in ejecta and implanted into the target. We applied astronomical
photometry techniques, using the SOURCE EXTRACTOR software, to automatically
measure the dimensions of thousands of fragments. This procedure enabled us to
estimate the implanted mass on the target body, which was found to be a few
percent of the initial mass of the impactor. We calculated an order of
magnitude difference in the energy density of catastrophic disruption, Q*,
between peridot and basalt projectiles. However, we found very similar
behaviour of the size frequency distributions for the hypervelocity shots (>1
km/s). After each shot, we examined the largest peridot fragments with Raman
spectroscopy and no melt or alteration in the final state of the projectile was
observed.
[7]
oai:arXiv.org:0808.3968 [pdf] - 15766
XAX: a multi-ton, multi-target detection system for dark matter, double
beta decay and pp solar neutrinos
Arisaka, K.;
Wang, H.;
Smith, P. F.;
Cline, D.;
Teymourian, A.;
Brown, E.;
Ooi, W.;
Aharoni, D.;
Lam, C. W.;
Lung, K.;
Davies, S.;
Price, M.
Submitted: 2008-08-28, last modified: 2009-01-06
A multi-target detection system XAX, comprising concentric 10 ton targets of
136Xe and 129/131Xe, together with a geometrically similar or larger target of
liquid Ar, is described. Each is configured as a two-phase
scintillation/ionization TPC detector, enhanced by a full 4pi array of
ultra-low radioactivity Quartz Photon Intensifying Detectors (QUPIDs) replacing
the conventional photomultipliers for detection of scintillation light. It is
shown that background levels in XAX can be reduced to the level required for
dark matter particle (WIMP) mass measurement at a 10^-10 pb WIMP-nucleon cross
section, with single-event sensitivity below 10^-11 pb. The use of multiple
target elements allows for confirmation of the A^2 dependence of a coherent
cross section, and the different Xe isotopes provide information on the
spin-dependence of the dark matter interaction. The event rates observed by Xe
and Ar would modulate annually with opposite phases from each other for WIMP
mass >~100 GeV/c^2. The large target mass of 136Xe and high degree of
background reduction allow neutrinoless double beta decay to be observed with
lifetimes of 10^27-10^28 years, corresponding to the Majorana neutrino mass
range 0.01-0.1 eV, the most likely range from observed neutrino mass
differences. The use of a 136Xe-depleted 129/131Xe target will also allow
measurement of the pp solar neutrino spectrum to a precision of 1-2%.
[8]
oai:arXiv.org:astro-ph/9403059 [pdf] - 1234204
A Search for Extragalactic Methanol Masers
Submitted: 1994-03-27
A sensitive search for 6.7--GHz methanol maser emission has been made towards
10 galaxies that have yielded detectable microwave molecular--line transitions.
These include several which show OH megamaser or superluminous \water\/ maser
emission. Within the Galaxy, \methanol\/ and OH masers often occur in the same
star formation regions and, in most cases, the \methanol\/ masers have a
greater peak flux density than their OH counterparts. Thus we might expect
\methanol\/ masers to be associated with extragalactic OH maser sources. We
failed to detect any emission or absorption above our 60--mJy detection limit.
We conclude that if the physical conditions exist to produce \methanol\/
megamaser emission, they are incompatible with the conditions which produce OH
megamaser emission.