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

Adam Amara

Adam Amara 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.IMastro-ph.GAastro-ph.EPastro-phphysics.comp-phApplicationsastro-ph.HEastro-ph.SRComputationMathematical SoftwareProgramming Languages

Catalog footprint

What is connected

50works
12topics
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

50 published item(s)

preprint2022arXiv

Probing the Cosmic Dark Ages with the Lunar Crater Radio Telescope

The Cosmic Dark Ages represent the period in the early evolution of the Universe, starting immediately after the decoupling of CMB photons from matter, and ending with the formation of the first stars and galaxies. The HI signal from the neutral hydrogen atoms is the only mechanism for us to understand this crucial phase in the cosmological history of the Universe and answer fundamental questions about the validity of the standard cosmological model, dark matter physics, and inflation. Due to cosmological redshift, this signal is now only observable in the 3-30~MHz frequency band, which is blocked from reaching the surface of the Earth by the ionosphere. In this paper, we present the design of the Lunar Crater Radio Telescope that intends to carry out unprecedented measurements of this signal by deploying a kilometer-sized parabolic reflector mesh inside a lunar crater on the far side of the Moon and suspending a receiver at its focus.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

SOAR/Goodman Spectroscopic Assessment of Candidate Counterparts of the LIGO-Virgo Event GW190814

On 2019 August 14 at 21:10:39 UTC, the LIGO/Virgo Collaboration (LVC) detected a possible neutron star-black hole merger (NSBH), the first ever identified. An extensive search for an optical counterpart of this event, designated GW190814, was undertaken using the Dark Energy Camera (DECam) on the 4m Victor M. Blanco Telescope at the Cerro Tololo Inter-American Observatory. Target of Opportunity interrupts were issued on 8 separate nights to observe 11 candidates using the 4.1m Southern Astrophysical Research (SOAR) telescope's Goodman High Throughput Spectrograph in order to assess whether any of these transients was likely to be an optical counterpart of the possible NSBH merger. Here, we describe the process of observing with SOAR, the analysis of our spectra, our spectroscopic typing methodology, and our resultant conclusion that none of the candidates corresponded to the gravitational wave merger event but were all instead other transients. Finally, we describe the lessons learned from this effort. Application of these lessons will be critical for a successful community spectroscopic follow-up program for LVC observing run 4 (O4) and beyond.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Likelihood-free Forward Modeling for Cluster Weak Lensing and Cosmology

Likelihood-free inference provides a rigorous approach to preform Bayesian analysis using forward simulations only. The main advantage of likelihood-free methods is its ability to account for complex physical processes and observational effects in forward simulations. Here we explore the potential of likelihood-free forward modeling for Bayesian cosmological inference using the redshift evolution of the cluster abundance combined with weak-lensing mass calibration. We use two complementary likelihood-free methods, namely Approximate Bayesian Computation (ABC) and Density-Estimation Likelihood-Free Inference (DELFI), to develop an analysis procedure for inference of the cosmological parameters $(Ω_\mathrm{m},σ_8)$ and the mass scale of the survey sample. Adopting an eROSITA-like selection function and a 10-percent scatter in the observable-mass relation in a flat $Λ$CDM cosmology with $Ω_\mathrm{m}=0.286$ and $σ_8=0.82$, we create a synthetic catalog of observable-selected NFW clusters in a survey area of 50 deg$^2$. The stacked tangential shear profile and the number counts in redshift bins are used as summary statistics for both methods. By performing a series of forward simulations, we obtain convergent solutions for the posterior distribution from both methods. We find that ABC recovers broader posteriors than DELFI, especially for the $Ω_\mathrm{m}$ parameter. For a weak-lensing survey with a source density of $n_\mathrm{g}=20$ arcmin$^{-2}$, we obtain posterior constraints on $S_8=σ_8(Ω_\mathrm{m}/0.3)^{0.3}$ of $0.836 \pm 0.032$ and $0.810 \pm 0.019$ from ABC and DELFI, respectively. The analysis framework developed in this study will be particularly powerful for cosmological inference with ongoing cluster cosmology programs, such as the XMM-XXL survey and the eROSITA all-sky survey, in combination with wide-field weak-lensing surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cross-correlating 21 cm and galaxy surveys: implications for cosmology and astrophysics

We forecast astrophysical and cosmological parameter constraints from synergies between 21 cm intensity mapping and wide field optical galaxy surveys (both spectroscopic and photometric) over $z \sim 0-3$. We focus on the following survey combinations in this work: (i) a CHIME-like and DESI-like survey in the northern hemisphere, (ii) an LSST-like and SKA I MID-like survey and (ii) a MeerKAT-like and DES-like survey in the southern hemisphere. We work with the $Λ$CDM cosmological model having parameters $\{h, Ω_m, n_s, Ω_b, σ_8\}$, parameters $v_{c,0}$ and $β$ representing the cutoff and slope of the HI-halo mass relation in the previously developed HI halo model framework, and a parameter $Q$ that represents the scale dependence of the optical galaxy bias. Using a Fisher forecasting framework, we explore (i) the effects of the HI and galaxy astrophysical uncertainties on the cosmological parameter constraints, assuming priors from the present knowledge of the astrophysics, (ii) the improvements on astrophysical constraints over their current priors in the three configurations considered, (ii) the tightening of the constraints on the parameters relative to the corresponding HI auto-correlation surveys alone.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spectro-Imaging Forward Model of Red and Blue Galaxies

For the next generation of spectroscopic galaxy surveys, it is important to forecast their performances and to accurately interpret their large data sets. For this purpose, it is necessary to consistently simulate different populations of galaxies, in particular Emission Line Galaxies (ELGs), less used in the past for cosmological purposes. In this work, we further the forward modeling approach presented in Fagioli et al. 2018, by extending the spectra simulator Uspec to model galaxies of different kinds with improved parameters from Tortorelli et al. 2020. Furthermore, we improve the modeling of the selection function by using the image simulator Ufig. We apply this to the Sloan Digital Sky Survey (SDSS), and simulate $\sim157,000$ multi-band images. We pre-process and analyse them to apply cuts for target selection, and finally simulate SDSS/BOSS DR14 galaxy spectra. We compute photometric, astrometric and spectroscopic properties for red and blue, real and simulated galaxies, finding very good agreement. We compare the statistical properties of the samples by decomposing them with Principal Component Analysis (PCA). We find very good agreement for red galaxies and a good, but less pronounced one, for blue galaxies, as expected given the known difficulty of simulating those. Finally, we derive stellar population properties, mass-to-light ratios, ages and metallicities, for all samples, finding again very good agreement. This shows how this method can be used not only to forecast cosmology surveys, but it is also able to provide insights into studies of galaxy formation and evolution.

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

The PAU Survey: Background light estimation with deep learning techniques

In any imaging survey, measuring accurately the astronomical background light is crucial to obtain good photometry. This paper introduces BKGnet, a deep neural network to predict the background and its associated error. BKGnet has been developed for data from the Physics of the Accelerating Universe Survey (PAUS), an imaging survey using a 40 narrow-band filter camera (PAUCam). Images obtained with PAUCam are affected by scattered light: an optical effect consisting of light multiply that deposits energy in specific detector regions contaminating the science measurements. Fortunately, scattered light is not a random effect, but it can be predicted and corrected for. We have found that BKGnet background predictions are very robust to distorting effects, while still being statistically accurate. On average, the use of BKGnet improves the photometric flux measurements by 7% and up to 20% at the bright end. BKGnet also removes a systematic trend in the background error estimation with magnitude in the i-band that is present with the current PAU data management method. With BKGnet, we reduce the photometric redshift outlier rate

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

The PAU Survey: A Forward Modeling Approach for Narrow-band Imaging

Weak gravitational lensing is a powerful probe of the dark sector, once measurement systematic errors can be controlled. In Refregier & Amara (2014), a calibration method based on forward modeling, called MCCL, was proposed. This relies on fast image simulations (e.g., UFig; Berge et al. 2013) that capture the key features of galaxy populations and measurement effects. The MCCL approach has been used in Herbel et al. (2017) to determine the redshift distribution of cosmological galaxy samples and, in the process, the authors derived a model for the galaxy population mainly based on broad-band photometry. Here, we test this model by forward modeling the 40 narrow-band photometry given by the novel PAU Survey (PAUS). For this purpose, we apply the same forced photometric pipeline on data and simulations using Source Extractor (Bertin & Arnouts 1996). The image simulation scheme performance is assessed at the image and at the catalogues level. We find good agreement for the distribution of pixel values, the magnitudes, in the magnitude-size relation and the interband correlations. A principal component analysis is then performed, in order to derive a global comparison of the narrow-band photometry between the data and the simulations. We use a `mixing' matrix to quantify the agreement between the observed and simulated sets of Principal Components (PCs). We find good agreement, especially for the first three most significant PCs. We also compare the coefficients of the PCs decomposition. While there are slight differences for some coefficients, we find that the distributions are in good agreement. Together, our results show that the galaxy population model derived from broad-band photometry is in good overall agreement with the PAUS data. This offers good prospect for incorporating spectral information to the galaxy model by adjusting it to the PAUS narrow-band data using forward modeling.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A new method to measure galaxy bias by combining the density and weak lensing fields

We present a new method to measure the redshift-dependent galaxy bias by combining information from the galaxy density field and the weak lensing field. This method is based on Amara et al. (2012), where they use the galaxy density field to construct a bias-weighted convergence field kg. The main difference between Amara et al. (2012) and our new implementation is that here we present another way to measure galaxy bias using tomography instead of bias parameterizations. The correlation between kg and the true lensing field k allows us to measure galaxy bias using different zero-lag correlations, such as <kgk>/<kk> or <kgkg>/<kgk>. Our method measures the linear bias factor on linear scales under the assumption of no stochasticity between galaxies and matter. We use the MICE simulation to measure the linear galaxy bias for a flux-limited sample (i < 22.5) in tomographic redshift bins using this method. This paper is the first that studies the accuracy and systematic uncertainties associated with the implementation of the method, and the regime where it is consistent with the linear galaxy bias defined by projected 2-point correlation functions (2PCF). We find that our method is consistent with linear bias at the percent level for scales larger than 30 arcmin, while nonlinearities appear at smaller scales. This measurement is a good complement to other measurements of bias, since it does not depend strongly on sigma8 as the 2PCF measurements. We apply this method to the Dark Energy Survey Science Verification data in a follow-up paper.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An Integrated System at the Bleien Observatory for Mapping the Galaxy

We describe the design and performance of the hardware system at the Bleien Observatory. The system is designed to deliver a map of the Galaxy for studying the foreground contamination of low-redshift (z=0.13--0.43) H$_{\rm I}$ intensity mapping experiments as well as other astronomical Galactic studies. This hardware system is composed of a 7m parabolic dish, a dual-polarization corrugated horn feed, a pseudo correlation receiver, a Fast Fourier Transform spectrometer, and an integrated control system that controls and monitors the progress of the data collection. The main innovative designs in the hardware are (1) the pseudo correlation receiver and the cold reference source within (2) the high dynamic range, high frequency resolution spectrometer and (3) the phase-switch implementation of the system. This is the first time these technologies are used together for a L-band radio telescope to achieve an electronically stable system, which is an essential first step for wide-field cosmological measurements. This work demonstrates the prospects and challenges for future H$_{\rm I}$ intensity mapping experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmology from large scale galaxy clustering and galaxy-galaxy lensing with Dark Energy Survey Science Verification data

We present cosmological constraints from the Dark Energy Survey (DES) using a combined analysis of angular clustering of red galaxies and their cross-correlation with weak gravitational lensing of background galaxies. We use a 139 square degree contiguous patch of DES data from the Science Verification (SV) period of observations. Using large scale measurements, we constrain the matter density of the Universe as Omega_m = 0.31 +/- 0.09 and the clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/- 0.13 after marginalizing over seven nuisance parameters and three additional cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16} = 0.74 +/- 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 = 0.78 +/- 0.09 using two bins over the range 0.2 <z< 0.5. We study the robustness of the results under changes in the data vectors, modelling and systematics treatment, including photometric redshift and shear calibration uncertainties, and find consistency in the derived cosmological parameters. We show that our results are consistent with previous cosmological analyses from DES and other data sets and conclude with a joint analysis of DES angular clustering and galaxy-galaxy lensing with Planck CMB data, Baryon Accoustic Oscillations and Supernova type Ia measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Integrated approach to cosmology: Combining CMB, large-scale structure and weak lensing

Recent observational progress has led to the establishment of the standard $Λ$CDM model for cosmology. This development is based on different cosmological probes that are usually combined through their likelihoods at the latest stage in the analysis. We implement here an integrated scheme for cosmological probes, which are combined in a common framework starting at the map level. This treatment is necessary as the probes are generally derived from overlapping maps and are thus not independent. It also allows for a thorough test of the cosmological model and of systematics through the consistency of different physical tracers. As a first application, we combine current measurements of the Cosmic Microwave Background (CMB) from the Planck satellite, and galaxy clustering and weak lensing from SDSS. We consider the spherical harmonic power spectra of these probes including all six auto- and cross-correlations along with the associated full Gaussian covariance matrix. This provides an integrated treatment of different analyses usually performed separately including CMB anisotropies, cosmic shear, galaxy clustering, galaxy-galaxy lensing and the Integrated Sachs-Wolfe (ISW) effect with galaxy and shear tracers. We derive constraints on $Λ$CDM parameters that are compatible with existing constraints and highlight tensions between data sets, which become apparent in this integrated treatment. We discuss how this approach provides a complete and powerful integrated framework for probe combination and how it can be extended to include other tracers in the context of current and future wide field cosmological surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantifying Concordance in Cosmology

Quantifying the concordance between different cosmological experiments is important for testing the validity of theoretical models and systematics in the observations. In earlier work, we thus proposed the Surprise, a concordance measure derived from the relative entropy between posterior distributions. We revisit the properties of the Surprise and describe how it provides a general, versatile, and robust measure for the agreement between datasets. We also compare it to other measures of concordance that have been proposed for cosmology. As an application, we extend our earlier analysis and use the Surprise to quantify the agreement between WMAP 9, Planck 13 and Planck 15 constraints on the $Λ$CDM model. Using a principle component analysis in parameter space, we find that the large Surprise between WMAP 9 and Planck 13 (S = 17.6 bits, implying a deviation from consistency at 99.8% confidence) is due to a shift along a direction that is dominated by the amplitude of the power spectrum. The Planck 15 constraints deviate from the Planck 13 results (S = 56.3 bits), primarily due to a shift in the same direction. The Surprise between WMAP and Planck consequently disappears when moving to Planck 15 (S = -5.1 bits). This means that, unlike Planck 13, Planck 15 is not in tension with WMAP 9. These results illustrate the advantages of the relative entropy and the Surprise for quantifying the disagreement between cosmological experiments and more generally as an information metric for cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Simulating the Large-Scale Structure of HI Intensity Maps

Intensity mapping of neutral hydrogen (HI) is a promising observational probe of cosmology and large-scale structure. We present wide field simulations of HI intensity maps based on N-body simulations of a $2.6\, {\rm Gpc / h}$ box with $2048^3$ particles (particle mass $1.6 \times 10^{11}\, {\rm M_\odot / h}$). Using a conditional mass function to populate the simulated dark matter density field with halos below the mass resolution of the simulation ($10^{8}\, {\rm M_\odot / h} < M_{\rm halo} < 10^{13}\, {\rm M_\odot / h}$), we assign HI to those halos according to a phenomenological halo to HI mass relation. The simulations span a redshift range of 0.35 < z < 0.9 in redshift bins of width $Δz \approx 0.05$ and cover a quarter of the sky at an angular resolution of about 7'. We use the simulated intensity maps to study the impact of non-linear effects and redshift space distortions on the angular clustering of HI. Focusing on the autocorrelations of the maps, we apply and compare several estimators for the angular power spectrum and its covariance. We verify that these estimators agree with analytic predictions on large scales and study the validity of approximations based on Gaussian random fields, particularly in the context of the covariance. We discuss how our results and the simulated maps can be useful for planning and interpreting future HI intensity mapping surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The mass-sheet degeneracy and time-delay cosmography: Analysis of the strong lens RXJ1131-1231

We present extended modeling of the strong lens system RXJ1131-1231 with archival data in two HST bands in combination with existing line-of-sight contribution and velocity dispersion estimates. Our focus is on source size and its influence on time-delay cosmography. We therefore examine the impact of mass-sheet degeneracy and especially the degeneracy pointed out by Schneider & Sluse (2013) using the source reconstruction scale. We also extend on previous work by further exploring the effects of priors on the kinematics of the lens and the external convergence in the environment of the lensing system. Our results coming from RXJ1131-1231 are given in a simple analytic form so that they can be easily combined with constraints coming from other cosmological probes. We find that the choice of priors on lens model parameters and source size are subdominant for the statistical errors for $H_0$ measurements of this systems. The choice of prior for the source is sub-dominant at present (2% uncertainty on $H_0$) but may be relevant for future studies. More importantly, we find that the priors on the kinematic anisotropy of the lens galaxy have a significant impact on our cosmological inference. When incorporating all the above modeling uncertainties, we find $H_0 = 86.6^{+6.8}_{-6.9}$ km s$^{-1}$Mpc$^{-1}$, when using kinematic priors similar to other studies. When we use a different kinematic prior motivated by Barnabè et al. (2012) but covering the same anisotropic range, we find $H_0 = 74.5^{+8.0}_{-7.8}$ km s$^{-1}$Mpc$^{-1}$. This means that the choice of kinematic modeling and priors have a significant impact on cosmographic inferences. The way forward is either to get better velocity dispersion measures which would down weight the impact of the priors or to construct physically motivated priors for the velocity dispersion model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Approximate Bayesian Computation for Forward Modeling in Cosmology

Bayesian inference is often used in cosmology and astrophysics to derive constraints on model parameters from observations. This approach relies on the ability to compute the likelihood of the data given a choice of model parameters. In many practical situations, the likelihood function may however be unavailable or intractable due to non-gaussian errors, non-linear measurements processes, or complex data formats such as catalogs and maps. In these cases, the simulation of mock data sets can often be made through forward modeling. We discuss how Approximate Bayesian Computation (ABC) can be used in these cases to derive an approximation to the posterior constraints using simulated data sets. This technique relies on the sampling of the parameter set, a distance metric to quantify the difference between the observation and the simulations and summary statistics to compress the information in the data. We first review the principles of ABC and discuss its implementation using a Population Monte-Carlo (PMC) algorithm and the Mahalanobis distance metric. We test the performance of the implementation using a Gaussian toy model. We then apply the ABC technique to the practical case of the calibration of image simulations for wide field cosmological surveys. We find that the ABC analysis is able to provide reliable parameter constraints for this problem and is therefore a promising technique for other applications in cosmology and astrophysics. Our implementation of the ABC PMC method is made available via a public code release.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Beam calibration of radio telescopes with drones

We present a multi-frequency far-field beam map for the 5m dish telescope at the Bleien Observatory measured using a commercially available drone. We describe the hexacopter drone used in this experiment, the design of the flight pattern, and the data analysis scheme. This is the first application of this calibration method to a single dish radio telescope in the far-field. The high signal-to-noise data allows us to characterise the beam pattern with high accuracy out to at least the 4th side-lobe. The resulting 2D beam pattern is compared with that derived from a more traditional calibration approach using an astronomical calibration source. We discuss the advantages of this method compared to other beam calibration methods. Our results show that this drone-based technique is very promising for ongoing and future radio experiments, where the knowledge of the beam pattern is key to obtaining high-accuracy cosmological and astronomical measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Calibrated Ultra Fast Image Simulations for the Dark Energy Survey

Weak lensing by large-scale structure is a powerful technique to probe the dark components of the universe. To understand the measurement process of weak lensing and the associated systematic effects, image simulations are becoming increasingly important. For this purpose we present a first implementation of the $\textit{Monte Carlo Control Loops}$ ($\textit{MCCL}$; Refregier & Amara 2014), a coherent framework for studying systematic effects in weak lensing. It allows us to model and calibrate the shear measurement process using image simulations from the Ultra Fast Image Generator (UFig; Berge et al. 2013). We apply this framework to a subset of the data taken during the Science Verification period (SV) of the Dark Energy Survey (DES). We calibrate the UFig simulations to be statistically consistent with DES images. We then perform tolerance analyses by perturbing the simulation parameters and study their impact on the shear measurement at the one-point level. This allows us to determine the relative importance of different input parameters to the simulations. For spatially constant systematic errors and six simulation parameters, the calibration of the simulation reaches the weak lensing precision needed for the DES SV survey area. Furthermore, we find a sensitivity of the shear measurement to the intrinsic ellipticity distribution, and an interplay between the magnitude-size and the pixel value diagnostics in constraining the noise model. This work is the first application of the $\textit{MCCL}$ framework to data and shows how it can be used to methodically study the impact of systematics on the cosmic shear measurement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Confirmation and characterization of the protoplanet HD100546 b - Direct evidence for gas giant planet formation at 50 au

We present the first multi-wavelength, high-contrast imaging study confirming the protoplanet embedded in the disk around the Herbig Ae/Be star HD100546. The object is detected at $L'$ ($\sim 3.8\,μm$) and $M'$ ($\sim 4.8\,μm$), but not at $K_s$ ($\sim 2.1\,μm$), and the emission consists of a point source component surrounded by spatially resolved emission. For the point source component we derive apparent magnitudes of $L'=13.92\pm0.10$ mag, $M'=13.33\pm0.16$ mag, and $K_s>15.43\pm0.11$ mag (3$σ$ limit), and a separation and position angle of $(0.457\pm0.014)"$ and $(8.4\pm1.4)^\circ$, and $(0.472\pm0.014)"$ and $(9.2\pm1.4)^\circ$ in $L'$ and $M'$, respectively. We demonstrate that the object is co-moving with HD100546 and can reject any (sub-)stellar fore-/background object. Fitting a single temperature blackbody to the observed fluxes of the point source component yields an effective temperature of $T_{eff}=932^{+193}_{-202}$ K and a radius for the emitting area of $R=6.9^{+2.7}_{-2.9}$ R$_{\rm Jupiter}$. The best-fit luminosity is $L=(2.3^{+0.6}_{-0.4})\cdot 10^{-4}\,L_{\rm Sun}$. We quantitatively compare our findings with predictions from evolutionary and atmospheric models for young, gas giant planets, discuss the possible existence of a warm, circumplanetary disk, and note that the de-projected physical separation from the host star of $(53\pm2)$ au poses a challenge standard planet formation theories. Considering the suspected existence of an additional planet orbiting at $\sim$13--14 au, HD100546 appears to be an unprecedented laboratory to study the formation of multiple gas giant planets empirically.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gravitational lens modeling with basis sets

We present a strong lensing modeling technique based on versatile basis sets for the lens and source planes. Our method uses high performance Monte Carlo algorithms, allows for an adaptive build up of complexity and bridges the gap between parametric and pixel based reconstruction methods. We apply our method to a HST image of the strong lens system RXJ1131-1231 and show that our method finds a reliable solution and is able to detect substructure in the lens and source planes simultaneously. Using mock data we show that our method is sensitive to sub-clumps with masses four orders of magnitude smaller than the main lens, which corresponds to about $10^8 M_{\odot}$, without prior knowledge on the position and mass of the sub-clump. The modelling approach is flexible and maximises automation to facilitate the analysis of the large number of strong lensing systems expected in upcoming wide field surveys. The resulting search for dark sub-clumps in these systems, without mass-to-light priors, offers promise for probing physics beyond the standard model in the dark matter sector.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Information Gains from Cosmic Microwave Background Experiments

To shed light on the fundamental problems posed by Dark Energy and Dark Matter, a large number of experiments have been performed and combined to constrain cosmological models. We propose a novel way of quantifying the information gained by updates on the parameter constraints from a series of experiments which can either complement earlier measurements or replace them. For this purpose, we use the Kullback-Leibler divergence or relative entropy from information theory to measure differences in the posterior distributions in model parameter space from a pair of experiments. We apply this formalism to a historical series of Cosmic Microwave Background experiments ranging from Boomerang to WMAP, SPT, and Planck. Considering different combinations of these experiments, we thus estimate the information gain in units of bits and distinguish contributions from the reduction of statistical errors and the `surprise' corresponding to a significant shift of the parameters' central values. For this experiment series, we find individual relative entropy gains ranging from about 1 to 30 bits. In some cases, e.g. when comparing WMAP and Planck results, we find that the gains are dominated by the surprise rather than by improvements in statistical precision. We discuss how this technique provides a useful tool for both quantifying the constraining power of data from cosmological probes and detecting the tensions between experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

3-dimensional spherical analyses of cosmological spectroscopic surveys

Spectroscopic redshift surveys offer great prospects for constraining the dark sector in cosmology. Future surveys will however be both deep and wide and will thus require an analysis in 3-dimensional spherical geometry. We review and compare several methods which have been proposed in the literature for this purpose, focusing in particular on implementations of the spherical harmonic tomography (SHT) power spectrum $C^{i j}_{l}$ and the spherical Fourier Bessel (SFB) power spectrum $C_{l} (k, k')$. Using a Fisher analysis, we compare the forecasted constraints on cosmological parameters using these statistics. These constraints typically rely on approximations such as the Limber approximation and make specific choices in the numerical implementation of each statistic. Using a series of toy models, we explore the applicability of these approximations and study the sensitivity of the SHT and SFB statistics to the details of their implementation. In particular, we show that overlapping redshift bins may improve cosmological constraints using the SHT statistic when the number of bins is small, and that the SFB constraints are quite robust to changes in the assumed distance-redshift relation. We also find that the SHT can be tailored to be more sensitive to modes at redshifts close to the survey boundary, while the SFB appears better suited to capture information beyond the smooth shape of the power spectrum. In this context, we discuss the pros and cons of the different techniques and their impact on the design and analysis of future wide field spectroscopic surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A simple model linking galaxy and dark matter evolution

We construct a simple phenomenological model for the evolving galaxy population by incorporating pre-defined baryonic prescriptions into a dark matter hierarchical merger tree. Specifically the model is based on the simple gas-regulator model introduced by Lilly et al. 2013 coupled with the empirical quenching rules of Peng et al. 2010/12. The simplest model already does quite well in reproducing, without re-adjusting the input parameters, many observables including the Main Sequence sSFR-mass relation, the faint end slope of the galaxy mass function and the shape of the star-forming and passive mass functions. Compared with observations and/or the recent phenomenological model of Behroozi et al. 2013 based on epoch-dependent abundance-matching, our model also qualitatively reproduces the evolution of the Main Sequence sSFR(z) and SFRD(z) star formation rate density relations, the $M_s - M_h$ stellar-to-halo mass relation and also the $SFR - M_h$ relation. Quantitatively the evolution of sSFR(z) and SFRD(z) is not steep enough, the $M_s - M_h$ relation is not quite peaked enough and, surprisingly, the ratio of quenched to star-forming galaxies around M* is not quite high enough. We show that these deficiencies can simultaneously be solved by ad hoc allowing galaxies to re-ingest some of the gas previously expelled in winds, provided that this is done in a mass-dependent and epoch-dependent way. These allow the model galaxies to reduce an inherent tendency to saturate their star-formation efficiency. This emphasizes how efficient galaxies around M* are in converting baryons into stars and highlights the fact that quenching occurs just at the point when galaxies are rapidly approaching the maximum possible efficiency of converting baryons into stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Way Forward for Cosmic Shear: Monte-Carlo Control Loops

Weak lensing by large scale structure or 'cosmic shear' is a potentially powerful cosmological probe to shed new light on Dark Matter, Dark Energy and Modified Gravity. It is based on the weak distortions induced by large-scale structures on the observed shapes of distant galaxies through gravitational lensing. While the potentials of this purely gravitational effect are great, results from this technique have been hampered because the measurement of this weak effect is difficult and limited by systematics effects. In particular, a demanding step is the measurement of the weak lensing shear from wide field CCD images of galaxies. We describe the origin of the problem and propose a way forward for cosmic shear. Our proposed approach is based on Monte-Carlo Control Loops and draws upon methods widely used in particle physics and engineering. We describe the control loop scheme and show how it provides a calibration method based on fast image simulations tuned to reproduce the statistical properties of a specific cosmic shear data set. Through a series of iterative loops and diagnostic tests, the Monte Carlo image simulations are made robust to perturbations on modeling input parameters and thus to systematic effects. We discuss how this approach can make the problem tractable and unleash to full potential of cosmic shear for cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Baryonic effects on weak-lensing two-point statistics and its cosmological implications

We develop an extension of \textit{the Halo Model} that describes analytically the corrections to the matter power spectrum due to the physics of baryons. We extend these corrections to the weak-lensing shear angular power spectrum. Within each halo, our baryonic model accounts for: 1) a central galaxy, the major stellar component whose properties are derived from abundance matching techniques; 2) a hot plasma in hydrostatic equilibrium and 3) an adiabatically-contracted dark matter component. This analytic approach allows us to compare our model to the dark-matter-only case. Our basic assumptions are tested against the hydrodynamical simulations of Martizzi et. al. (2014), with which a remarkable agreement is found. Our baryonic model has only one free parameter, $M_{\rm crit}$, the critical halo mass that marks the transition between feedback-dominated halos, mostly devoid of gas, and gas rich halos, in which AGN feedback effects become weaker. We explore the entire cosmological parameter space, using the angular power spectrum in three redshift bins as the observable, assuming a Euclid-like survey. We derive the corresponding constraints on the cosmological parameters, as well as the possible bias introduced by neglecting the effects of baryonic physics. We find that, up to $\ell_{max}$=4000, baryonic physics plays very little role in the cosmological parameters estimation. However, if one goes up to $\ell_{max}$=8000, the marginalized errors on the cosmological parameters can be significantly reduced, but neglecting baryonic physics can lead to bias in the recovered cosmological parameters up to 10$σ$. These biases are removed if one takes into account the main baryonic parameter, $M_{\rm crit}$, which can also be determined up to 1-2\%, along with the other cosmological parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Discovery of a Companion Candidate in the HD169142 Transition Disk and the Possibility of Multiple Planet Formation

We present L' and J-band high-contrast observations of HD169142, obtained with the VLT/NACO AGPM vector vortex coronagraph and the Gemini Planet Imager, respectively. A source located at 0".156+/-0".032 north of the host star (PA=7.4+/-11.3 degrees) appears in the final reduced L' image. At the distance of the star (~145 pc), this angular separation corresponds to a physical separation of 22.7+/-4.7 AU, locating the source within the recently resolved inner cavity of the transition disk. The source has a brightness of L'=12.2+/-0.5 mag, whereas it is not detected in the J band (J>13.8 mag). If its L' brightness arose solely from the photosphere of a companion and given the J-L' color constraints, it would correspond to a 28-32 MJupiter object at the age of the star, according to the COND models. Ongoing accretion activity of the star suggests, however, that gas is left in the inner disk cavity from which the companion could also be accreting. In this case the object could be lower in mass and its luminosity enhanced by the accretion process and by a circumplanetary disk. A lower mass object is more consistent with the observed cavity width. Finally, the observations enable us to place an upper limit on the L'-band flux of a second companion candidate orbiting in the disk annular gap at ~50 AU, as suggested by millimeter observations. If the second companion is also confirmed, HD169142 might be forming a planetary system, with at least two companions opening gaps and possibly interacting with each other.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-contrast Imaging with Spitzer: Deep Observations of Vega, Fomalhaut, and epsilon Eridani

Stars with debris disks are intriguing targets for direct imaging exoplanet searches, both due to previous detections of wide planets in debris disk systems, as well as commonly existing morphological features in the disks themselves that may be indicative of a planetary influence. Here we present observations of three of the most nearby young stars, that are also known to host massive debris disks: Vega, Fomalhaut, and eps Eri. The Spitzer Space Telescope is used at a range of orientation angles for each star, in order to supply a deep contrast through angular differential imaging combined with high-contrast algorithms. The observations provide the opportunity to probe substantially colder bound planets (120--330 K) than is possible with any other technique or instrument. For Vega, some apparently very red candidate point sources detected in the 4.5 micron image remain to be tested for common proper motion. The images are sensitive to ~2 Mjup companions at 150 AU in this system. The observations presented here represent the first search for planets around Vega using Spitzer. The upper 4.5 micron flux limit on Fomalhaut b could be further constrained relative to previous data. In the case of eps Eri, planets below both the effective temperature and the mass of Jupiter could be probed from 80 AU and outwards, although no such planets were found. The data sensitively probe the regions around the edges of the debris rings in the systems where planets can be expected to reside. These observations validate previous results showing that more than an order of magnitude improvement in performance in the contrast-limited regime can be acquired with respect to conventional methods by applying sophisticated high-contrast techniques to space-based telescopes, thanks to the high degree of PSF stability provided in this environment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

HOPE: A Python Just-In-Time compiler for astrophysical computations

The Python programming language is becoming increasingly popular for scientific applications due to its simplicity, versatility, and the broad range of its libraries. A drawback of this dynamic language, however, is its low runtime performance which limits its applicability for large simulations and for the analysis of large data sets, as is common in astrophysics and cosmology. While various frameworks have been developed to address this limitation, most focus on covering the complete language set, and either force the user to alter the code or are not able to reach the full speed of an optimised native compiled language. In order to combine the ease of Python and the speed of C++, we developed HOPE, a specialised Python just-in-time (JIT) compiler designed for numerical astrophysical applications. HOPE focuses on a subset of the language and is able to translate Python code into C++ while performing numerical optimisation on mathematical expressions at runtime. To enable the JIT compilation, the user only needs to add a decorator to the function definition. We assess the performance of HOPE by performing a series of benchmarks and compare its execution speed with that of plain Python, C++ and the other existing frameworks. We find that HOPE improves the performance compared to plain Python by a factor of 2 to 120, achieves speeds comparable to that of C++, and often exceeds the speed of the existing solutions. We discuss the differences between HOPE and the other frameworks, as well as future extensions of its capabilities. The fully documented HOPE package is available at http://hope.phys.ethz.ch and is published under the GPLv3 license on PyPI and GitHub.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Model Breaking Measure for Cosmological Surveys

Recent observations have led to the establishment of the concordance LCDM model for cosmology. A number of experiments are being planned to shed light on dark energy, dark matter, inflation and gravity, which are the key components of the model. To optimize and compare the reach of these surveys, several figures of merit have been proposed. They are based on either the forecasted precision on the LCDM model and its expansion, or on the expected ability to distinguish two models. We propose here another figure of merit that quantifies the capacity of future surveys to rule out the LCDM model. It is based on a measure of the difference in volume of observable space that the future surveys will constrain with and without imposing the model. This model breaking figure of merit is easy to compute and can lead to different survey optimizations than other metrics. We illustrate its impact using a simple combination of supernovae and BAO mock observations and compare the respective merit of these probes to challenge LCDM. We discuss how this approach would impact the design of future cosmological experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

PynPoint Code for Exoplanet Imaging

We announce the public release of PynPoint, a Python package that we have developed for analysing exoplanet data taken with the angular differential imaging observing technique. In particular, PynPoint is designed to model the point spread function of the central star and to subtract its flux contribution to reveal nearby faint companion planets. The current version of the package does this correction by using a principal component analysis method to build a basis set for modelling the point spread function of the observations. We demonstrate the performance of the package by reanalysing publicly available data on the exoplanet beta Pictoris b, which consists of close to 24,000 individual image frames. We show that PynPoint is able to analyse this typical data in roughly 1.5 minutes on a Mac Pro, when the number of images is reduced by co-adding in sets of 5. The main computational work parallelises well as a result of a reliance on SciPy and NumPy functions. For this calculation the peak memory load is 6Gb, which can be run comfortably on most workstations. A simpler calculation, by co-adding over 50, takes 3 seconds with a peak memory usage of 600 Mb. This can be performed easily on a laptop. In developing the package we have modularised the code so that we will be able to extend functionality in future releases, through the inclusion of more modules, without it affecting the users application programming interface. We distribute the PynPoint package through the central PyPi sever, and the documentation is available online (http://pynpoint.ethz.ch).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A young protoplanet candidate embedded in the circumstellar disk of HD100546

We present high-contrast observations of the circumstellar environment of the Herbig Ae/Be star HD100546. The final 3.8 micron image reveals an emission source at a projected separation of 0.48"+-0.04" (corresponding to ~47+-4 AU at a position angle of 8.9+-0.9 degree. The emission appears slightly extended with a point source component with an apparent magnitude of 13.2+-0.4 mag. The position of the source coincides with a local deficit in polarization fraction in near-infrared polarimetric imaging data, which probes the surface of the well-studied circumstellar disk of HD100546. This suggests a possible physical link between the emission source and the disk. Assuming a disk inclination of ~47 degree the de-projected separation of the object is ~68 AU. Assessing the likelihood of various scenarios we favor an interpretation of the available high-contrast data with a planet in the process of forming. Follow-up observations in the coming years can easily distinguish between the different possible scenarios empirically. If confirmed, HD100546 "b" would be a unique laboratory to study the formation process of a new planetary system, with one giant planet currently forming in the disk and a second planet possibly orbiting in the disk gap at smaller separations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

CosmoHammer: Cosmological parameter estimation with the MCMC Hammer

We study the benefits and limits of parallelised Markov chain Monte Carlo (MCMC) sampling in cosmology. MCMC methods are widely used for the estimation of cosmological parameters from a given set of observations and are typically based on the Metropolis-Hastings algorithm. Some of the required calculations can however be computationally intensive, meaning that a single long chain can take several hours or days to calculate. In practice, this can be limiting, since the MCMC process needs to be performed many times to test the impact of possible systematics and to understand the robustness of the measurements being made. To achieve greater speed through parallelisation, MCMC algorithms need to have short auto-correlation times and minimal overheads caused by tuning and burn-in. The resulting scalability is hence influenced by two factors, the MCMC overheads and the parallelisation costs. In order to efficiently distribute the MCMC sampling over thousands of cores on modern cloud computing infrastructure, we developed a Python framework called CosmoHammer which embeds emcee, an implementation by Foreman-Mackey et al. (2012) of the affine invariant ensemble sampler by Goodman and Weare (2010). We test the performance of CosmoHammer for cosmological parameter estimation from cosmic microwave background data. While Metropolis-Hastings is dominated by overheads, CosmoHammer is able to accelerate the sampling process from a wall time of 30 hours on a dual core notebook to 16 minutes by scaling out to 2048 cores. Such short wall times for complex data sets opens possibilities for extensive model testing and control of systematics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An Ultra Fast Image Generator (UFig) for wide-field astronomy

Simulated wide-field images are becoming an important part of observational astronomy, either to prepare for new surveys or to test measurement methods. In order to efficiently explore vast parameter spaces, the computational speed of simulation codes is a central requirement to their implementation. We introduce the Ultra Fast Image Generator (UFig) which aims to bring wide-field imaging simulations to the current limits of computational capabilities. We achieve this goal through: (1) models of galaxies, stars and observational conditions, which, while simple, capture the key features necessary for realistic simulations, and (2) state-of-the-art computational and implementation optimizations. We present the performances of UFig and show that it is faster than existing public simulation codes by several orders of magnitude. It allows us to produce images more quickly than SExtractor needs to analyze them. For instance, it can simulate a typical 0.25 deg^2 Subaru SuprimeCam image (10k x 8k pixels) with a 5-sigma limiting magnitude of R=26 in 30 seconds on a laptop, yielding an average simulation time for a galaxy of 30 microseconds. This code is complementary to end-to-end simulation codes and can be used as a fast, central component of observational methods relying on simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraining primordial non-Gaussianity with future galaxy surveys

We study the constraining power on primordial non-Gaussianity of future surveys of the large-scale structure of the Universe for both near-term surveys (such as the Dark Energy Survey - DES) as well as longer term projects such as Euclid and WFIRST. Specifically we perform a Fisher matrix analysis forecast for such surveys, using DES-like and Euclid-like configurations as examples, and take account of any expected photometric and spectroscopic data. We focus on two-point statistics and we consider three observables: the 3D galaxy power spectrum in redshift space, the angular galaxy power spectrum, and the projected weak-lensing shear power spectrum. We study the effects of adding a few extra parameters to the basic LCDM set. We include the two standard parameters to model the current value for the dark energy equation of state and its time derivative, w_0, w_a, and we account for the possibility of primordial non-Gaussianity of the local, equilateral and orthogonal types, of parameter fNL and, optionally, of spectral index n_fNL. We present forecasted constraints on these parameters using the different observational probes. We show that accounting for models that include primordial non-Gaussianity does not degrade the constraint on the standard LCDM set nor on the dark-energy equation of state. By combining the weak lensing data and the information on projected galaxy clustering, consistently including all two-point functions and their covariance, we find forecasted marginalised errors sigma (fNL) ~ 3, sigma (n_fNL) ~ 0.12 from a Euclid-like survey for the local shape of primordial non-Gaussianity, while the orthogonal and equilateral constraints are weakened for the galaxy clustering case, due to the weaker scale-dependence of the bias. In the lensing case, the constraints remain instead similar in all configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Measurement and Calibration of Noise Bias in Weak Lensing Galaxy Shape Estimation

Weak gravitational lensing has the potential to constrain cosmological parameters to high precision. However, as shown by the Shear TEsting Programmes (STEP) and GRavitational lEnsing Accuracy Testing (GREAT) Challenges, measuring galaxy shears is a nontrivial task: various methods introduce different systematic biases which have to be accounted for. We investigate how pixel noise on the image affects the bias on shear estimates from a Maximum-Likelihood forward model-fitting approach using a sum of co-elliptical Sérsic profiles, in complement to the theoretical approach of an an associated paper. We evaluate the bias using a simple but realistic galaxy model and find that the effects of noise alone can cause biases of order 1-10% on measured shears, which is significant for current and future lensing surveys. We evaluate a simulation-based calibration method to create a bias model as a function of galaxy properties and observing conditions. This model is then used to correct the simulated measurements. We demonstrate that this method can effectively reduce noise bias so that shear measurement reaches the level of accuracy required for estimating cosmic shear in upcoming lensing surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Noise bias in weak lensing shape measurements

Weak lensing experiments are a powerful probe of cosmology through their measurement of the mass distribution of the universe. A challenge for this technique is to control systematic errors that occur when measuring the shapes of distant galaxies. In this paper we investigate noise bias, a systematic error that arises from second order noise terms in the shape measurement process. We first derive analytical expressions for the bias of general Maximum Likelihood Estimators (MLEs) in the presence of additive noise. We then find analytical expressions for a simplified toy model in which galaxies are modeled and fitted with a Gaussian with its size as a single free parameter. Even for this very simple case we find a significant effect. We also extend our analysis to a more realistic 6-parameter elliptical Gaussian model. We find that the noise bias is generically of the order of the inverse-squared signal-to-noise ratio (SNR) of the galaxies and is thus of the order of a percent for galaxies of SNR of 10, i.e. comparable to the weak lensing shear signal. This is nearly two orders of magnitude greater than the systematics requirements for future all-sky weak lensing surveys. We discuss possible ways to circumvent this effect, including a calibration method using simulations discussed in an associated paper.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

PynPoint: An Image Processing Package for Finding Exoplanets

We present the scientific performance results of PynPoint, our Python-based software package that uses principle component analysis to detect and estimate the flux of exoplanets in two dimensional imaging data. Recent advances in adaptive optics and imaging technology at visible and infrared wavelengths have opened the door to direct detections of planetary companions to nearby stars, but image processing techniques have yet to be optimized. We show that the performance of our approach gives a marked improvement over what is presently possible using existing methods such as LOCI. To test our approach, we use real angular differential imaging (ADI) data taken with the adaptive optics assisted high resolution near-infrared camera NACO at the VLT. These data were taken during the commissioning of the apodising phase plate (APP) coronagraph. By inserting simulated planets into these data, we test the performance of our method as a function of planet brightness for different positions on the image. We find that in all cases PynPoint has a detection threshold that is superior to that given by our LOCI analysis when assessed in a common statistical framework. We obtain our best improvements for smaller inner working angles (IWA). For an IWA of ~ 0.29" we find that we achieve a detection sensitivity that is a factor of 5 better than LOCI. We also investigate our ability to correctly measure the flux of planets. Again, we find improvements over LOCI, with PynPoint giving more stable results. Finally, we apply our package to a non-APP dataset of the exoplanet beta Pictoris b and reveal the planet with high signal-to-noise. This confirms that PynPoint can potentially be applied with high fidelity to a wide range of high-contrast imaging datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Space-quality data from balloon-borne telescopes: the High Altitude Lensing Observatory (HALO)

We present a method for attaining sub-arcsecond pointing stability during sub- orbital balloon flights, as designed for in the High Altitude Lensing Observatory (HALO) concept. The pointing method presented here has the potential to perform near-space quality optical astronomical imaging at 1-2% of the cost of space-based missions. We also discuss an architecture that can achieve sufficient thermomechanical stability to match the pointing stability. This concept is motivated by advances in the development and testing of Ultra Long Duration Balloon (ULDB) flights which promise to allow observation campaigns lasting more than three months. The design incorporates a multi-stage pointing architecture comprising: a gondola coarse azimuth control system, a multi-axis nested gimbal frame structure with arcsecond stability, a telescope de-rotator to eliminate field rotation, and a fine guidance stage consisting of both a telescope mounted angular rate sensor and guide CCDs in the focal plane to drive a fast-steering mirror. We discuss the results of pointing tests together with a preliminary thermo-mechanical analysis required for sub-arcsecond pointing at high altitude. Possible future applications in the areas of wide-field surveys and exoplanet searches are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational Lensing Accuracy Testing 2010 (GREAT10) Challenge Handbook

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

iCosmo: an Interactive Cosmology Package

Aims: The interactive software package iCosmo, designed to perform cosmological calculations is described. Methods: iCosmo is a software package to perform interactive cosmological calculations for the low redshift universe. Computing distance measures, the matter power spectrum, and the growth factor is supported for any values of the cosmological parameters. It also computes derived observed quantities for several cosmological probes such as cosmic shear, baryon acoustic oscillations and type Ia supernovae. The associated errors for these observables can be derived for customised surveys, or for pre-set values corresponding to current or planned instruments. The code also allows for the calculation of cosmological forecasts with Fisher matrices which can be manipulated to combine different surveys and cosmological probes. The code is written in the IDL language and thus benefits from the convenient interactive features and scientific library available in this language. iCosmo can also be used as an engine to perform cosmological calculations in batch mode, and forms a convenient adaptive platform for the development of further cosmological modules. With its extensive documentation, it may also serve as a useful resource for teaching and for newcomers in the field of cosmology. Results: The iCosmo package is described with various examples and command sequences. The code is freely available with documentation at http://www.icosmo.org, along with an interactive web interface and is part of the Initiative for Cosmology, a common archive for cosmological resources.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Image processing challenges in weak gravitational lensing

The field of weak gravitational lensing, which measures the basic properties of the Universe by studying the way that light from distant galaxies is perturbed as it travels towards us, is a very active field in astronomy. This short article presents a broad overview of the field, including some of the important questions that cosmologists are trying to address, such as understanding the nature of dark energy and dark matter. To do this, there is an increasing feeling within the weak lensing community that other disciplines, such as computer science, machine learning, signal processing and image processing, have the expertise that would bring enormous advantage if channelled into lensing studies. To illustrate this point, the article below outlines some of the key steps in a weak lensing analysis chain. The challenges are distinct at each step, but each could benefit from ideas developed in the signal processing domain. This article also gives a brief overview of current and planned lensing experiments that will soon bring about an influx of data sets that are substantially larger than those analysed to date. It is, therefore, inevitable that current techniques are likely to be insufficient, thus leading to an exciting era where new methods will become crucial for the continued success of the field.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On Point Spread Function modelling: towards optimal interpolation

Point Spread Function (PSF) modeling is a central part of any astronomy data analysis relying on measuring the shapes of objects. It is especially crucial for weak gravitational lensing, in order to beat down systematics and allow one to reach the full potential of weak lensing in measuring dark energy. A PSF modeling pipeline is made of two main steps: the first one is to assess its shape on stars, and the second is to interpolate it at any desired position (usually galaxies). We focus on the second part, and compare different interpolation schemes, including polynomial interpolation, radial basis functions, Delaunay triangulation and Kriging. For that purpose, we develop simulations of PSF fields, in which stars are built from a set of basis functions defined from a Principal Components Analysis of a real ground-based image. We find that Kriging gives the most reliable interpolation, significantly better than the traditionally used polynomial interpolation. We also note that although a Kriging interpolation on individual images is enough to control systematics at the level necessary for current weak lensing surveys, more elaborate techniques will have to be developed to reach future ambitious surveys' requirements.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Probe combination in large galaxy surveys : Application of Fisher information and Shannon entropy to weak lensing

This paper is aimed at developing a better understanding of the structure of the information that in contained in galaxy surveys, so as to find optimal ways to combine observables from such surveys. We first show how Jaynes' Maximal Entropy Principle allows us, in the general case, to express the Fisher information content of data sets in terms of the curvature of the Shannon entropy surface with respect to the relevant observables. This allows us to understand the Fisher information content of a data set, once a physical model is specified, independently of the specific way that the data will be processed, and without any assumptions of Gaussianity. This includes as a special case the standard Fisher matrix prescriptions for Gaussian variables widely used in the cosmological community, for instance for power spectra extraction. As an application of this approach, we evaluate the prospects of a joint analysis of weak lensing tracers up to second order in the shapes distortions, in the case that the noise in each probe can be effectively treated as model independent. These include the magnification, the two ellipticity and the four flexion fields. At the two point level, we show that the only effect of treating these observables in combination is a simple scale dependent decrease of the noise contaminating the accessible spectrum of the lensing E-mode. We provide simple bounds to its extraction by a combination of such probes, as well as its quantitative evaluation when the correlations between the noise variables for any two such probes can be neglected.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Figures of Merit for Testing Standard Models: Application to Dark Energy Experiments in Cosmology

Given a standard model to test, an experiment can be designed to: (i) measure the standard model parameters; (ii) extend the standard model; or (iii) look for evidence of deviations from the standard model. To measure (or extend) the standard model, the Fisher matrix is widely used in cosmology to predict expected parameter errors for future surveys under Gaussian assumptions. In this article, we present a frame- work that can be used to design experiments such that it maximises the chance of finding a deviation from the standard model. Using a simple illustrative example, discussed in the appendix, we show that the optimal experimental configuration can depend dramatically on the optimisation approach chosen. We also show some simple cosmology calculations, where we study Baryonic Acoustic Oscillation and Supernove surveys. In doing so, we also show how external data, such as the positions of the CMB peaks measured by WMAP, and theory priors can be included in the analysis. In the cosmological cases that we have studied (DETF Stage III), we find that the three optimisation approaches yield similar results, which is reassuring and indicates that the choice of optimal experiment is fairly robust at this level. However, this may not be the case as we move to more ambitious future surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Optimal capture of non-Gaussianity in weak lensing surveys: power spectrum, bispectrum and halo counts

We compare the efficiency of weak lensing-selected galaxy clusters counts and of the weak lensing bispectrum at capturing non-Gaussian features in the dark matter distribution. We use the halo model to compute the weak lensing power spectrum, the bispectrum and the expected number of detected clusters, and derive constraints on cosmological parameters for a large, low systematic weak lensing survey, by focusing on the $Ω_m$-$σ_8$ plane and on the dark energy equation of state. We separate the power spectrum into the resolved and the unresolved parts of the data, the resolved part being defined as detected clusters, and the unresolved part as the rest of the field. We consider four kinds of clusters counts, taking into account different amount of information : signal-to-noise ratio peak counts; counts as a function of clusters' mass; counts as a function of clusters' redshift; and counts as a function of clusters' mass and redshift. We show that when combined with the power spectrum, those four kinds of counts provide similar constraints, thus allowing one to perform the most direct counts, signal-to-noise peaks counts, and get percent level constraints on cosmological parameters. We show that the weak lensing bispectrum gives constraints comparable to those given by the power spectrum and captures non-Gaussian features as well as clusters counts, its combination with the power spectrum giving errors on cosmological parameters that are similar to, if not marginally smaller than, those obtained when combining the power spectrum with cluster counts. We finally note that in order to reach its potential, the weak lensing bispectrum must be computed using all triangle configurations, as equilateral triangles alone do not provide useful information.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Weak lensing forecasts for dark energy, neutrinos and initial conditions

Weak gravitational lensing provides a sensitive probe of cosmology by measuring the mass distribution and the geometry of the low redshift universe. We show how an all-sky weak lensing tomographic survey can jointly constrain different sets of cosmological parameters describing dark energy, massive neutrinos (hot dark matter), and the primordial power spectrum. In order to put all sectors on an equal footing, we introduce a new parameter $β$, the second order running spectral index. Using the Fisher matrix formalism with and without CMB priors, we examine how the constraints vary as the parameter set is enlarged. We find that weak lensing with CMB priors provides robust constraints on dark energy parameters and can simultaneously provide strong constraints on all three sectors. We find that the dark energy sector is largely insensitive to the inclusion of the other cosmological sectors. Implications for the planning of future surveys are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Photo-z Performance for Precision Cosmology

Current and future weak lensing surveys will rely on photometrically estimated redshifts of very large numbers of galaxies. In this paper, we address several different aspects of the demanding photo-z performance that will be required for future experiments, such as the proposed ESA Euclid mission. It is first shown that the proposed all-sky near-infrared photometry from Euclid, in combination with anticipated ground-based photometry (e.g. PanStarrs-2 or DES) should yield the required precision in individual photo-z of sigma(z) < 0.05(1+z) at I_AB < 24.5. Simple a priori rejection schemes based on the photometry alone can be tuned to recognise objects with wildly discrepant photo-z and to reduce the outlier fraction to < 0.25% with only modest loss of otherwise usable objects. Turning to the more challenging problem of determining the mean redshift <z> of a set of galaxies to a precision of 0.002(1+z) we argue that, for many different reasons, this is best accomplished by relying on the photo-z themselves rather than on the direct measurement of <z> from spectroscopic redshifts of a representative subset of the galaxies. A simple adaptive scheme based on the statistical properties of the photo-z likelihood functions is shown to meet this stringent systematic requirement. We also examine the effect of an imprecise correction for Galactic extinction and the effects of contamination by fainter over-lapping objects in photo-z determination. The overall conclusion of this work is that the acquisition of photometrically estimated redshifts with the precision required for Euclid, or other similar experiments, will be challenging but possible. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Combined analysis of weak lensing and X-ray blind surveys

We present a joint weak lensing and X-ray analysis of 4 deg$^2$ from the CFHTLS and XMM-LSS surveys. Our weak lensing analysis is the first analysis of a real survey using shapelets, a new generation weak lensing analysis method. We create projected mass maps of the images, and extract 6 weak-lensing-detected clusters of galaxies. We show that their counts can be used to constrain the power spectrum normalisation $σ_8 =0.92_{-0.30}^{+0.26}$ for $Ω_m=0.24$. We show that despite the large scatter generally observed in the M-T relation derived from lensing masses, tight constraints on both its slope and normalisation $M_*$ can be obtained with a moderate number of sources provided that the covered mass range is large enough. Adding clusters from Bardeau et al. (2007) to our sample, we measure $M_* = 2.71_{-0.61}^{+0.79} 10^{14} h^{-1} M_\odot$. Although they are dominated by shot noise and sample variance, our measurements are consistent with currently favoured values, and set the stage for future surveys. We thus investigate the dependence of those estimates on survey size, depth, and integration time, for joint weak lensing and X-ray surveys. We show that deep surveys should be dedicated to the study of the physics of clusters and groups of galaxies. For a given exposure time, wide surveys provide a larger number of detected clusters and are therefore preferred for the measurement of cosmological parameters such as $σ_8$ and $M_*$. We show that a wide survey of a few hundred square degrees is needed to improve upon current measurements of these parameters. More ambitious surveys covering 7000 deg$^2$ will provide the 1% accuracy in the estimation of the power spectrum and the M-T relation normalisations.

0 citations0 reviewsNo code linkedOpen access