Normalized to: Balan, S.

**15** article(s) in total. **91** co-authors, from **1** to **6** common article(s). Median position in authors list is **2,0**.

[1] oai:arXiv.org:1811.02578 [pdf] - 1949312

Upper Bound of Neutrino Masses from Combined Cosmological Observations
and Particle Physics Experiments

Comments: 5 pages, 2 figures

Submitted: **2018-11-06**, last modified: **2019-08-27**

We investigate the impact of prior models on the upper bound of the sum of
neutrino masses, $\sum m_{\nu}$. We use data from Large Scale Structure of
galaxies, Cosmic Microwave Background, Type Ia SuperNovae, and Big Bang
Nucleosynthesis. We probe physically motivated neutrino mass models (respecting
oscillation experiment constraints) and compare them to constraints using
standard cosmological approximations. The former give a consistent upper bound
of $\sum m_{\nu} \lesssim 0.26$ eV ($95\%$ CI) and yields a strong competitive
upper bound for the lightest neutrino mass species, $m_0^{\nu} < 0.086$ eV
($95\%$ CI). By contrast one of the approximations, which is somewhat
inconsistent with oscillation experiments, yields an upper bound of $\sum
m_{\nu} \lesssim 0.15$ eV ($95\%$ CI), which differs substantially from the
former upper bound. We, therefore, argue that cosmological neutrino mass and
hierarchy determination should be pursued using physically motivated models
since approximations might lead to incorrect and nonphysical upper bounds.

[2] oai:arXiv.org:1809.07204 [pdf] - 1949306

Cosmological Measurements from Angular Power Spectra Analysis of BOSS
DR12 Tomography

Comments: 28 pages, 28 figures, 3 appendices

Submitted: **2018-09-19**, last modified: **2019-08-27**

We constrain cosmological parameters by analysing the angular power spectra
of the Baryon Oscillation Spectroscopic Survey DR12 galaxies, a spectroscopic
follow-up of around 1.3 million SDSS galaxies over 9,376 deg$^2$ with an
effective volume of $\sim 6.5$ (Gpc $h^{-1}$)$^3$ in the redshift range $0.15
\leq z < 0.80$. We split this sample into 13 tomographic bins ($\Delta z =
0.05$); angular power spectra were calculated using a Pseudo-$C_{\ell}$
estimator, and covariance matrices were estimated using log-normal simulated
maps. Cosmological constraints obtained from these data were combined with
constraints from Planck CMB experiment as well as the JLA supernovae
compilation. Considering a $w$CDM cosmological model measured on scales up to
$k_{max} = 0.07h$ Mpc$^{-1}$, we constrain a constant dark energy
equation-of-state with a $\sim 4\%$ error at the 1-$\sigma$ level: $w_0 =
-0.993^{+0.046}_{-0.043}$, together with $\Omega_m = 0.330\pm 0.012$, $\Omega_b
= 0.0505 \pm 0.002$, $S_8 \equiv \sigma_8 \sqrt{\Omega_m/0.3} = 0.863 \pm
0.016$, and $h = 0.661 \pm 0.012$. For the same combination of datasets, but
now considering a $\Lambda$CDM model with massive neutrinos and the same scale
cut, we find: $\Omega_m = 0.328 \pm 0.009$, $\Omega_b =
0.05017^{+0.0009}_{-0.0008}$, $S_8 = 0.862 \pm 0.017$, and $h =
0.663^{+0.006}_{-0.007}$ and a 95\% credible interval (CI) upper limit of $\sum
m_{\nu} < 0.14$ eV for a normal hierarchy. These results are competitive if not
better than standard analyses with the same dataset, and demonstrate this
should be a method of choice for future surveys, opening the door for their
full exploitation in cross-correlations probes.

[3] oai:arXiv.org:1805.06799 [pdf] - 1777870

Radio Galaxy Shape Measurement with Hamiltonian Monte Carlo in the
Visibility Domain

Comments: 15 pages, 10 figures, 3 tables

Submitted: **2018-05-17**, last modified: **2018-11-01**

Radio weak lensing, while a highly promising complementary probe to optical
weak lensing, will require incredible precision in the measurement of galaxy
shape parameters. In this paper, we extend the Bayesian Inference for Radio
Observations model fitting approach to measure galaxy shapes directly from
visibility data of radio continuum surveys, instead of from image data. We
apply a Hamiltonian Monte Carlo (HMC) technique for sampling the posterior,
which is more efficient than the standard Monte Carlo Markov Chain method when
dealing with a large dimensional parameter space. Adopting the exponential
profile for galaxy model fitting allows us to analytically calculate the
likelihood gradient required by HMC, allowing a faster and more accurate
sampling. The method is tested on SKA1-MID simulated observations at 1.4 GHz of
a field containing up to 1000 star-forming galaxies. It is also applied to a
simulated observation of the weak lensing precursor survey SuperCLASS. In both
cases we obtain reliable measurements of the galaxies' ellipticity and size for
all sources with SNR $\ge 10$, and we also find relationships between the
convergence properties of the HMC technique and some source parameters. Direct
shape measurement in the visibility domain achieves high accuracy at the
expected source number densities of the current and next SKA precursor
continuum surveys. The proposed method can be easily extended for the fitting
of other galaxy and scientific parameters, as well as simultaneously
marginalising over systematic and instrumental effects.

[4] oai:arXiv.org:1612.00307 [pdf] - 1525851

A Joint Analysis for Cosmology and Photometric Redshift Calculation
Using Cross Correlations

Comments: 12 Pages, 6 figures

Submitted: **2016-12-01**

We present a method of calibrating the properties of photometric redshift
bins as part of a larger Markov Chain Monte Carlo (MCMC) analysis for the
inference of cosmological parameters. The redshift bins are characterised by
their mean and variance, which are varied as free parameters and marginalised
over when obtaining the cosmological parameters. We demonstrate that the
likelihood function for cross-correlations in an angular power spectrum
framework tightly constrains the properties of bins such that they may be well
determined, reducing their influence on cosmological parameters and avoiding
the bias from poorly estimated redshift distributions. We demonstrate that even
with only three photometric and three spectroscopic bins, we can recover
accurate estimates of the mean redshift of a bin to within $\Delta\mu \approx
3-4 \times10^{-3}$ and the width of the bin to $\Delta\sigma \approx
1\times10^{-3}$ for galaxies near $z = 1$. This indicates that we may be able
to bring down the photometric redshift errors to a level which is in line with
the requirements for the next generation of cosmological experiments.

[5] oai:arXiv.org:1602.01839 [pdf] - 1426241

Large scale distribution of total mass versus luminous matter from
Baryon Acoustic Oscillations: First search in the SDSS-III BOSS Data Release
10

Comments: 6 pages, 2 figures, accepted for publication in PRL

Submitted: **2016-02-04**

Baryon Acoustic Oscillations (BAOs) in the early Universe are predicted to
leave an as yet undetected signature on the relative clustering of total mass
versus luminous matter. A detection of this effect would provide an important
confirmation of the standard cosmological paradigm and constrain alternatives
to dark matter as well as non-standard fluctuations such as Compensated
Isocurvature Perturbations (CIPs). We conduct the first observational search
for this effect, by comparing the number-weighted and luminosity-weighted
correlation functions, using the SDSS-III BOSS Data Release 10 CMASS sample.
When including CIPs in our model, we formally obtain evidence at $3.2\sigma$ of
the relative clustering signature and a limit that matches the existing upper
limits on the amplitude of CIPs. However, various tests suggest that these
results are not yet robust, perhaps due to systematic biases in the data. The
method developed in this Letter, used with more accurate future data such as
that from DESI, is likely to confirm or disprove our preliminary evidence.

[6] oai:arXiv.org:1210.3578 [pdf] - 684311

Hyper-efficient model-independent Bayesian method for the analysis of
pulsar timing data

Comments: 9 pages, 4 figures

Submitted: **2012-10-12**, last modified: **2013-06-25**

A new model independent method is presented for the analysis of pulsar timing
data and the estimation of the spectral properties of an isotropic
gravitational wave background (GWB). We show that by rephrasing the likelihood
we are able to eliminate the most costly aspects of computation normally
associated with this type of data analysis. When applied to the International
Pulsar Timing Array Mock Data Challenge data sets this results in speedups of
approximately 2 to 3 orders of magnitude compared to established methods. We
present three applications of the new likelihood. In the low signal to noise
regime we sample directly from the power spectrum coefficients of the GWB
signal realization. In the high signal to noise regime, where the data can
support a large number of coefficients, we sample from the joint probability
density of the power spectrum coefficients for the individual pulsars and the
GWB signal realization. Critically in both these cases we need make no
assumptions about the form of the power spectrum of the GWB, or the individual
pulsars. Finally we present a method for characterizing the spatial correlation
between pulsars on the sky, making no assumptions about the form of that
correlation, and therefore providing the only truly general Bayesian method of
confirming a GWB detection from pulsar timing data.

[7] oai:arXiv.org:1202.5254 [pdf] - 630580

Image Analysis for Cosmology: Results from the GREAT10 Galaxy Challenge

Comments: Accepted to MNRAS

Submitted: **2012-02-23**, last modified: **2013-02-25**

In this paper we present results from the weak lensing shape measurement
GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. This
marks an order of magnitude step change in the level of scrutiny employed in
weak lensing shape measurement analysis. We provide descriptions of each method
tested and include 10 evaluation metrics over 24 simulation branches. GREAT10
was the first shape measurement challenge to include variable fields; both the
shear field and the Point Spread Function (PSF) vary across the images in a
realistic manner. The variable fields enable a variety of metrics that are
inaccessible to constant shear simulations including a direct measure of the
impact of shape measurement inaccuracies, and the impact of PSF size and
ellipticity, on the shear power spectrum. To assess the impact of shape
measurement bias for cosmic shear we present a general pseudo-Cl formalism,
that propagates spatially varying systematics in cosmic shear through to power
spectrum estimates. We also show how one-point estimators of bias can be
extracted from variable shear simulations. The GREAT10 Galaxy Challenge
received 95 submissions and saw a factor of 3 improvement in the accuracy
achieved by shape measurement methods. The best methods achieve sub-percent
average biases. We find a strong dependence in accuracy as a function of
signal-to-noise, and indications of a weak dependence on galaxy type and size.
Some requirements for the most ambitious cosmic shear experiments are met above
a signal-to-noise ratio of 20. These results have the caveat that the simulated
PSF was a ground-based PSF. Our results are a snapshot of the accuracy of
current shape measurement methods and are a benchmark upon which improvement
can continue. This provides a foundation for a better understanding of the
strengths and limitations of shape measurement methods.

[8] oai:arXiv.org:1201.0668 [pdf] - 605254

A Uniformly Derived Catalogue of Exoplanets from Radial Velocities

Comments: 16 pages, 6 figures, 6 tables

Submitted: **2012-01-03**

A new catalogue of extrasolar planets is presented by re-analysing a
selection of published radial velocity data sets using EXOFIT (Balan & Lahav
2009). All objects are treated on an equal footing within a Bayesian framework,
to give orbital parameters for 94 exoplanetary systems. Model selection
(between one- and two-planet solutions) is then performed, using both a visual
flagging method and a standard chi-square analysis, with agreement between the
two methods for 99% of the systems considered. The catalogue is to be made
available online, and this 'proof of concept' study may be maintained and
extended in the future to incorporate all systems with publicly available
radial velocity data, as well as transit and microlensing data.

[9] oai:arXiv.org:1009.0779 [pdf] - 445530

Gravitational Lensing Accuracy Testing 2010 (GREAT10) Challenge Handbook

Comments: Published in at http://dx.doi.org/10.1214/11-AOAS484 the Annals of
Applied Statistics (http://www.imstat.org/aoas/) by the Institute of
Mathematical Statistics (http://www.imstat.org)

Submitted: **2010-09-03**, last modified: **2011-11-30**

GRavitational lEnsing Accuracy Testing 2010 (GREAT10) is a public image
analysis challenge aimed at the development of algorithms to analyze
astronomical images. Specifically, the challenge is to measure varying image
distortions in the presence of a variable convolution kernel, pixelization and
noise. This is the second in a series of challenges set to the astronomy,
computer science and statistics communities, providing a structured environment
in which methods can be improved and tested in preparation for planned
astronomical surveys. GREAT10 extends upon previous work by introducing
variable fields into the challenge. The "Galaxy Challenge" involves the precise
measurement of galaxy shape distortions, quantified locally by two parameters
called shear, in the presence of a known convolution kernel. Crucially, the
convolution kernel and the simulated gravitational lensing shape distortion
both now vary as a function of position within the images, as is the case for
real data. In addition, we introduce the "Star Challenge" that concerns the
reconstruction of a variable convolution kernel, similar to that in a typical
astronomical observation. This document details the GREAT10 Challenge for
potential participants. Continually updated information is also available from
http://www.greatchallenges.info.

[10] oai:arXiv.org:1105.1150 [pdf] - 402430

Bayesian evidence for two companions orbiting HIP 5158

Comments: 5 pages, 4 figures, 3 tables. Added a couple of figures showing the
residuals after one and two companion fits. Accepted for publication in MNRAS
Letters

Submitted: **2011-05-05**, last modified: **2011-07-14**

We present results of a Bayesian analysis of radial velocity (RV) data for
the star HIP 5158, confirming the presence of two companions and also
constraining their orbital parameters. Assuming Keplerian orbits, the
two-companion model is found to be e^{48} times more probable than the
one-planet model, although the orbital parameters of the second companion are
only weakly constrained. The derived orbital periods are 345.6 +/- 2.0 d and
9017.8 +/- 3180.7 d respectively, and the corresponding eccentricities are 0.54
+/- 0.04 and 0.14 +/- 0.10. The limits on planetary mass (m \sin i) and
semimajor axis are (1.44 +/- 0.14 M_{J}, 0.89 +/- 0.01 AU) and (15.04 +/- 10.55
M_{J}, 7.70 +/- 1.88 AU) respectively. Owing to large uncertainty on the mass
of the second companion, we are unable to determine whether it is a planet or a
brown dwarf. The remaining `noise' (stellar jitter) unaccounted for by the
model is 2.28 +/- 0.31 m/s. We also analysed a three-companion model, but found
it to be e^{8} times less probable than the two-companion model.

[11] oai:arXiv.org:1012.5129 [pdf] - 402410

Detecting extrasolar planets from stellar radial velocities using
Bayesian evidence

Comments: 12 pages, 6 figures, 12 tables. Added results obtained by applying
the method to 47 Ursae Majoris and HD 10180. Accepted for publication in
MNRAS

Submitted: **2010-12-22**, last modified: **2011-05-04**

Stellar radial velocity (RV) measurements have proven to be a very successful
method for detecting extrasolar planets. Analysing RV data to determine the
parameters of the extrasolar planets is a significant statistical challenge
owing to the presence of multiple planets and various degeneracies between
orbital parameters. Determining the number of planets favoured by the observed
data is an even more difficult task. Bayesian model selection provides a
mathematically rigorous solution to this problem by calculating marginal
posterior probabilities of models with different number of planets, but the use
of this method in extrasolar planetary searches has been hampered by the
computational cost of the evaluating Bayesian evidence. Nonetheless, Bayesian
model selection has the potential to improve the interpretation of existing
observational data and possibly detect yet undiscovered planets. We present a
new and efficient Bayesian method for determining the number of extrasolar
planets, as well as for inferring their orbital parameters, without having to
calculate directly the Bayesian evidence for models containing a large number
of planets. Instead, we work iteratively and at each iteration obtain a
conservative lower limit on the odds ratio for the inclusion of an additional
planet into the model. We apply this method to simulated data-sets containing
one and two planets and successfully recover the correct number of planets and
reliable constraints on the orbital parameters. We also apply our method to RV
measurements of HD 37124, 47 Ursae Majoris and HD 10180. For HD 37124, we
confirm that the current data strongly favour a three-planet system. We find
strong evidence for the presence of a fourth planet in 47 Ursae Majoris, but
its orbital period is suspiciously close to one year, casting doubt on its
validity. For HD 10180 we find strong evidence for a six-planet system.

[12] oai:arXiv.org:0908.0945 [pdf] - 1003324

Results of the GREAT08 Challenge: An image analysis competition for
cosmological lensing

Comments: 10 figures, 11 pages of main text, 9 pages of appendices. Submitted
to MNRAS. Challenge website at http://www.great08challenge.info

Submitted: **2009-08-07**

We present the results of the GREAT08 Challenge, a blind analysis challenge
to infer weak gravitational lensing shear distortions from images. The primary
goal was to stimulate new ideas by presenting the problem to researchers
outside the shear measurement community. Six GREAT08 Team methods were
presented at the launch of the Challenge and five additional groups submitted
results during the 6 month competition. Participants analyzed 30 million
simulated galaxies with a range in signal to noise ratio, point-spread function
ellipticity, galaxy size, and galaxy type. The large quantity of simulations
allowed shear measurement methods to be assessed at a level of accuracy
suitable for currently planned future cosmic shear observations for the first
time. Different methods perform well in different parts of simulation parameter
space and come close to the target level of accuracy in several of these. A
number of fresh ideas have emerged as a result of the Challenge including a
re-examination of the process of combining information from different galaxies,
which reduces the dependence on realistic galaxy modelling. The image
simulations will become increasingly sophisticated in future GREAT challenges,
meanwhile the GREAT08 simulations remain as a benchmark for additional
developments in shear measurement algorithms.

[13] oai:arXiv.org:0907.3613 [pdf] - 26516

EXOFIT: Bayesian Estimation of Orbital Parameters of Extrasolar Planets

Comments: To appear in the Proceedings of Molecules in the Atmospheres of
Extrasolar Planets. 7 pages, 2 tables and 4 Figures

Submitted: **2009-07-21**

We introduce EXOFIT, a Bayesian tool for estimating orbital parameters of
extrasolar planets from radial velocity measurements. EXOFIT can search for
either one or two planets at present. EXOFIT employs Markov Chain Monte Carlo
method implemented in an object oriented manner. As an example we re-analyze
the orbital solution of HD155358 and the results are compared with that of the
published orbital parameters. In order to check the agreement of the EXOFIT
orbital parameters with the published ones we examined radial velocity data of
30 stars taken randomly from www.exoplanet.eu. We show that while orbital
periods agree in both methods, EXOFIT prefers lower eccentricity solutions for
planets with higher (e >=0.5) orbital eccentricities.

[14] oai:arXiv.org:0802.1214 [pdf] - 9903

Handbook for the GREAT08 Challenge: An image analysis competition for
cosmological lensing

Comments: Published in at http://dx.doi.org/10.1214/08-AOAS222 the Annals of
Applied Statistics (http://www.imstat.org/aoas/) by the Institute of
Mathematical Statistics (http://www.imstat.org)

Submitted: **2008-02-11**, last modified: **2009-06-15**

The GRavitational lEnsing Accuracy Testing 2008 (GREAT08) Challenge focuses
on a problem that is of crucial importance for future observations in
cosmology. The shapes of distant galaxies can be used to determine the
properties of dark energy and the nature of gravity, because light from those
galaxies is bent by gravity from the intervening dark matter. The observed
galaxy images appear distorted, although only slightly, and their shapes must
be precisely disentangled from the effects of pixelisation, convolution and
noise. The worldwide gravitational lensing community has made significant
progress in techniques to measure these distortions via the Shear TEsting
Program (STEP). Via STEP, we have run challenges within our own community, and
come to recognise that this particular image analysis problem is ideally
matched to experts in statistical inference, inverse problems and computational
learning. Thus, in order to continue the progress seen in recent years, we are
seeking an infusion of new ideas from these communities. This document details
the GREAT08 Challenge for potential participants. Please visit
http://www.great08challenge.info for the latest information.

[15] oai:arXiv.org:0805.3532 [pdf] - 12878

ExoFit: Orbital Parameters of Extra-solar Planets from Radial Velocities

Comments: Accepted to MNRAS, 11 pages, 10 figures. Minor changes in Section 2,
Appendix removed . ExoFit can be downloaded from
http://www.star.ucl.ac.uk/~lahav/exofit.html

Submitted: **2008-05-23**, last modified: **2009-03-16**

Retrieval of orbital parameters of extrasolar planets poses considerable
statistical challenges.Due to sparse sampling, measurement errors, parameters
degeneracy and modelling limitations, there are no unique values of basic
parameters, such as period and eccentricity. Here, we estimate the orbital
parameters from radial velocity data in a Bayesian framework by utilizing
Markov Chain Monte Carlo (MCMC) simulations with the Metropolis-Hastings
algorithm. We follow a methodology recently proposed by Gregory and Ford. Our
implementation of MCMC is based on the object-oriented approach outlined by
Graves. We make our resulting code, ExoFit, publicly available with this paper.
It can search for either one or two planets as illustrated on mock data. As an
example we re-analysed the orbital solution of companions to HD 187085 and HD
159868 from the published radial velocity data. We confirm the degeneracy
reported for orbital parameters of the companion to HD 187085, and show that a
low-eccentricity orbit is more probable for this planet. For HD 159868, we
obtained slightly different orbital solution and a relatively high 'noise'
factor indicating the presence of an unaccounted signal in the radial velocity
data. ExoFit is designed in such a way that it can be extended for a variety of
probability models, including different Bayesian priors.