Performance of internal Covariance Estimators for Cosmic Shear Correlation Functions

Data re-sampling methods such as the delete-one jackknife are a common tool for estimating the covariance of large scale structure probes. In this paper we investigate the concepts of internal covariance estimation in the context of cosmic shear two-point statistics. We demonstrate how to use log-normal simulations of the convergence field and the corresponding shear field to carry out realistic tests of internal covariance estimators and find that most estimators such as jackknife or sub-sample covariance can reach a satisfactory compromise between bias and variance of the estimated covariance. In a forecast for the complete, 5-year DES survey we show that internally estimated covariance matrices can provide a large fraction of the true uncertainties on cosmological parameters in a 2D cosmic shear analysis. The volume inside contours of constant likelihood in the $\Omega_m$-$\sigma_8$ plane as measured with internally estimated covariance matrices is on average $\gtrsim 85\%$ of the volume derived from the true covariance matrix. The uncertainty on the parameter combination $\Sigma_8 \sim \sigma_8 \Omega_m^{0.5}$ derived from internally estimated covariances is $\sim 90\%$ of the true uncertainty.Comment: submitted to mnra

Weak Lensing Reconstruction and Power Spectrum Estimation: Minimum Variance Methods

Large-scale structure distorts the images of background galaxies, which allows one to measure directly the projected distribution of dark matter in the universe and determine its power spectrum. Here we address the question of how to extract this information from the observations. We derive minimum variance estimators for projected density reconstruction and its power spectrum and apply them to simulated data sets, showing that they give a good agreement with the theoretical minimum variance expectations. The same estimator can also be applied to the cluster reconstruction, where it remains a useful reconstruction technique, although it is no longer optimal for every application. The method can be generalized to include nonlinear cluster reconstruction and photometric information on redshifts of background galaxies in the analysis. We also address the question of how to obtain directly the 3-d power spectrum from the weak lensing data. We derive a minimum variance quadratic estimator, which maximizes the likelihood function for the 3-d power spectrum and can be computed either from the measurements directly or from the 2-d power spectrum. The estimator correctly propagates the errors and provides a full correlation matrix of the estimates. It can be generalized to the case where redshift distribution depends on the galaxy photometric properties, which allows one to measure both the 3-d power spectrum and its time evolution.Comment: revised version, 36 pages, AAS LateX, submitted to Ap

A Comparison of Simple Mass Estimators for Galaxy Clusters

High-resolution N-body simulations are used to investigate systematic trends in the mass profiles and total masses of clusters as derived from 3 simple estimators: (1) the weak gravitational lensing shear field under the assumption of an isothermal cluster potential, (2) the dynamical mass obtained from the measured velocity dispersion under the assumption of an isothermal cluster potential, and (3) the classical virial estimator. The clusters consist of order 2.5e+05 particles of mass m_p \simeq 10^{10} \Msun, have triaxial mass distributions, and significant substructure exists within their virial radii. Not surprisingly, the level of agreement between the mass profiles obtained from the various estimators and the actual mass profiles is found to be scale-dependent. The virial estimator yields a good measurement of the total cluster mass, though it is systematically underestimated by of order 10%. This result suggests that, at least in the limit of ideal data, the virial estimator is quite robust to deviations from pure spherical symmetry and the presence of substructure. The dynamical mass estimate based upon a measurement of the cluster velocity dispersion and an assumption of an isothermal potential yields a poor measurement of the total mass. The weak lensing estimate yields a very good measurement of the total mass, provided the mean shear used to determine the equivalent cluster velocity dispersion is computed from an average of the lensing signal over the entire cluster (i.e. the mean shear is computed interior to the virial radius). [abridged]Comment: Accepted for publication in The Astrophysical Journal. Complete paper, including 3 large colour figures can also be obtained from http://bu-ast.bu.edu/~brainerd/preprints

Determination of the Earth's pole tide Love number k₂ from observations of polar motion using an adaptive Kalman filter approach

The geophysical interpretation of observed time series of Earth rotation parameters (ERP) is commonly based on numerical models that describe and balance variations of angular momentum in various subsystems of the Earth. Naturally, models are dependent on geometrical, rheological and physical parameters. Many of these are weakly determined from other models or observations. In our study we present an adaptive Kalman filter approach for the improvement of parameters of the dynamic Earth system model DyMEG which acts as a simulator of ERP. In particular we focus on the improvement of the pole tide Love number k(2). In the frame of a sensitivity analysis k(2) has been identified as one of the most crucial parameters of DyMEG since it directly influences the modeled Chandler oscillation. At the same time k(2) is one of the most uncertain parameters in the model. Our simulations with DyMEG cover a period of 60 years after which a steady state of k(2) is reached. The estimate for k(2), accounting for the anelastic response of the Earth's mantle and the ocean, is 0.3531 + 0.0030i. We demonstrate that the application of the improved parameter k(2) in DyMEG leads to significantly better results for polar motion than the original value taken from the Conventions of the International Earth Rotation and Reference Systems Service (IERS)

Automatic detection of arcs and arclets formed by gravitational lensing

We present an algorithm developed particularly to detect gravitationally lensed arcs in clusters of galaxies. This algorithm is suited for automated surveys as well as individual arc detections. New methods are used for image smoothing and source detection. The smoothing is performed by so-called anisotropic diffusion, which maintains the shape of the arcs and does not disperse them. The algorithm is much more efficient in detecting arcs than other source finding algorithms and the detection by eye.Comment: A&A in press, 12 pages, 16 figure

Rank 3 permutation characters and maximal subgroups

In this paper we classify all maximal subgroups M of a nearly simple primitive rank 3 group G of type L=Omega_{2m+1}(3), m > 3; acting on an L-orbit E of non-singular points of the natural module for L such that 1_P^G <=1_M^G where P is a stabilizer of a point in E. This result has an application to the study of minimal genera of algebraic curves which admit group actions.Comment: 41 pages, to appear in Forum Mathematicu

The GALATEA Test-Facility for High Purity Germanium Detectors

GALATEA is a test facility designed to investigate bulk and surface effects in high purity germanium detectors. A vacuum tank houses an infrared screened volume with a cooled detector inside. A system of three stages allows an almost complete scan of the detector. The main feature of GALATEA is that there is no material between source and detector. This allows the usage of alpha and beta sources as well as of a laser beam to study surface effects. A 19-fold segmented true-coaxial germanium detector was used for commissioning

Probing the Universe with Weak Lensing

Gravitational lenses can provide crucial information on the geometry of the Universe, on the cosmological scenario of formation of its structures as well as on the history of its components with look-back time. In this review, I focus on the most recent results obtained during the last five years from the analysis of the weak lensing regime. The interest of weak lensing as a probe of dark matter and the for study of the coupling between light and mass on scales of clusters of galaxies, large scale structures and galaxies is discussed first. Then I present the impact of weak lensing for the study of distant galaxies and of the population of lensed sources as function of redshift. Finally, I discuss the potential interest of weak lensing to constrain the cosmological parameters, either from pure geometrical effects observed in peculiar lenses, or from the coupling of weak lensing with the CMB.Comment: To appear Annual Review of Astronomy and Astrophysiscs Vol. 37. Latex and psfig.sty. Version without figure, 54 pages, 73Kb. Complete version including 13 figures (60 pages) available on ftp.iap.fr anonymous account in /pub/from_users/mellier/AnnualReview ; file ARAAmellier.ps.gz 1.6 M

Simplifying Random Satisfiability Problem by Removing Frustrating Interactions

How can we remove some interactions in a constraint satisfaction problem (CSP) such that it still remains satisfiable? In this paper we study a modified survey propagation algorithm that enables us to address this question for a prototypical CSP, i.e. random K-satisfiability problem. The average number of removed interactions is controlled by a tuning parameter in the algorithm. If the original problem is satisfiable then we are able to construct satisfiable subproblems ranging from the original one to a minimal one with minimum possible number of interactions. The minimal satisfiable subproblems will provide directly the solutions of the original problem.Comment: 21 pages, 16 figure