Source author record

Stella Seitz

Stella Seitz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

34 published item(s)

preprint2016arXiv

Anomaly detection for machine learning redshifts applied to SDSS galaxies

We present an analysis of anomaly detection for machine learning redshift estimation. Anomaly detection allows the removal of poor training examples, which can adversely influence redshift estimates. Anomalous training examples may be photometric galaxies with incorrect spectroscopic redshifts, or galaxies with one or more poorly measured photometric quantity. We select 2.5 million 'clean' SDSS DR12 galaxies with reliable spectroscopic redshifts, and 6730 'anomalous' galaxies with spectroscopic redshift measurements which are flagged as unreliable. We contaminate the clean base galaxy sample with galaxies with unreliable redshifts and attempt to recover the contaminating galaxies using the Elliptical Envelope technique. We then train four machine learning architectures for redshift analysis on both the contaminated sample and on the preprocessed 'anomaly-removed' sample and measure redshift statistics on a clean validation sample generated without any preprocessing. We find an improvement on all measured statistics of up to 80% when training on the anomaly removed sample as compared with training on the contaminated sample for each of the machine learning routines explored. We further describe a method to estimate the contamination fraction of a base data sample.

preprint2016arXiv

Stacking for machine learning redshifts applied to SDSS galaxies

We present an analysis of a general machine learning technique called 'stacking' for the estimation of photometric redshifts. Stacking techniques can feed the photometric redshift estimate, as output by a base algorithm, back into the same algorithm as an additional input feature in a subsequent learning round. We shown how all tested base algorithms benefit from at least one additional stacking round (or layer). To demonstrate the benefit of stacking, we apply the method to both unsupervised machine learning techniques based on self-organising maps (SOMs), and supervised machine learning methods based on decision trees. We explore a range of stacking architectures, such as the number of layers and the number of base learners per layer. Finally we explore the effectiveness of stacking even when using a successful algorithm such as AdaBoost. We observe a significant improvement of between 1.9% and 21% on all computed metrics when stacking is applied to weak learners (such as SOMs and decision trees). When applied to strong learning algorithms (such as AdaBoost) the ratio of improvement shrinks, but still remains positive and is between 0.4% and 2.5% for the explored metrics and comes at almost no additional computational cost.

preprint2016arXiv

The evolution of metallicity and metallicity gradients from z=2.7-0.6 with KMOS3D

We present measurements of the [NII]/Ha ratio as a probe of gas-phase oxygen abundance for a sample of 419 star-forming galaxies at z=0.6-2.7 from the KMOS3D near-IR multi-IFU survey. The mass-metallicity relation (MZR) is determined consistently with the same sample selection, metallicity tracer, and methodology over the wide redshift range probed by the survey. We find good agreement with long-slit surveys in the literature, except for the low-mass slope of the relation at z~2.3, where this sample is less biased than previous samples based on optical spectroscopic redshifts. In this regime we measure a steeper slope than some literature results. Excluding the AGN contribution from the MZR reduces sensitivity at the high mass end, but produces otherwise consistent results. There is no significant dependence of the [NII]/Ha ratio on SFR or environment at fixed redshift and stellar mass. The IFU data allow spatially resolved measurements of [NII]/Ha, from which we can infer abundance gradients for 180 galaxies, thus tripling the current sample in the literature. The observed gradients are on average flat, with only 15 gradients statistically offset from zero at >3sigma. We have modelled the effect of beam-smearing, assuming a smooth intrinsic radial gradient and known seeing, inclination and effective radius for each galaxy. Our seeing-limited observations can recover up to 70% of the intrinsic gradient for the largest, face-on disks, but only 30% for the smaller, more inclined galaxies. We do not find significant trends between observed or corrected gradients and any stellar population, dynamical or structural galaxy parameters, mostly in agreement with existing studies with much smaller sample sizes. In cosmological simulations, strong feedback is generally required to produce flat gradients at high redshift.

preprint2016arXiv

Tuning target selection algorithms to improve galaxy redshift estimates

We showcase machine learning (ML) inspired target selection algorithms to determine which of all potential targets should be selected first for spectroscopic follow up. Efficient target selection can improve the ML redshift uncertainties as calculated on an independent sample, while requiring less targets to be observed. We compare the ML targeting algorithms with the Sloan Digital Sky Survey (SDSS) target order, and with a random targeting algorithm. The ML inspired algorithms are constructed iteratively by estimating which of the remaining target galaxies will be most difficult for the machine learning methods to accurately estimate redshifts using the previously observed data. This is performed by predicting the expected redshift error and redshift offset (or bias) of all of the remaining target galaxies. We find that the predicted values of bias and error are accurate to better than 10-30% of the true values, even with only limited training sample sizes. We construct a hypothetical follow-up survey and find that some of the ML targeting algorithms are able to obtain the same redshift predictive power with 2-3 times less observing time, as compared to that of the SDSS, or random, target selection algorithms. The reduction in the required follow up resources could allow for a change to the follow-up strategy, for example by obtaining deeper spectroscopy, which could improve ML redshift estimates for deeper test data.

preprint2015arXiv

Accurate photometric redshift probability density estimation - method comparison and application

We introduce an ordinal classification algorithm for photometric redshift estimation, which significantly improves the reconstruction of photometric redshift probability density functions (PDFs) for individual galaxies and galaxy samples. As a use case we apply our method to CFHTLS galaxies. The ordinal classification algorithm treats distinct redshift bins as ordered values, which improves the quality of photometric redshift PDFs, compared with non-ordinal classification architectures. We also propose a new single value point estimate of the galaxy redshift, that can be used to estimate the full redshift PDF of a galaxy sample. This method is competitive in terms of accuracy with contemporary algorithms, which stack the full redshift PDFs of all galaxies in the sample, but requires orders of magnitudes less storage space. The methods described in this paper greatly improve the log-likelihood of individual object redshift PDFs, when compared with a popular Neural Network code (ANNz). In our use case, this improvement reaches 50\% for high redshift objects ($z \geq 0.75$). We show that using these more accurate photometric redshift PDFs will lead to a reduction in the systematic biases by up to a factor of four, when compared with less accurate PDFs obtained from commonly used methods. The cosmological analyses we examine and find improvement upon are the following: gravitational lensing cluster mass estimates, modelling of angular correlation functions, and modelling of cosmic shear correlation functions.

preprint2015arXiv

Characterization and correction of charge-induced pixel shifts in DECam

Interaction of charges in CCDs with the already accumulated charge distribution causes both a flux dependence of the point-spread function (an increase of observed size with flux, also known as the brighter/fatter effect) and pixel-to-pixel correlations of the Poissonian noise in flat fields. We describe these effects in the Dark Energy Camera (DECam) with charge dependent shifts of effective pixel borders, i.e. the Antilogus et al. (2014) model, which we fit to measurements of flat-field Poissonian noise correlations. The latter fall off approximately as a power-law r^-2.5 with pixel separation r, are isotropic except for an asymmetry in the direct neighbors along rows and columns, are stable in time, and are weakly dependent on wavelength. They show variations from chip to chip at the 20% level that correlate with the silicon resistivity. The charge shifts predicted by the model cause biased shape measurements, primarily due to their effect on bright stars, at levels exceeding weak lensing science requirements. We measure the flux dependence of star images and show that the effect can be mitigated by applying the reverse charge shifts at the pixel level during image processing. Differences in stellar size, however, remain significant due to residuals at larger distance from the centroid.

preprint2015arXiv

CLASH-X: A Comparison of Lensing and X-ray Techniques for Measuring the Mass Profiles of Galaxy Clusters

We present profiles of temperature (Tx), gas mass, and hydrostatic mass estimated from new and archival X-ray observations of CLASH clusters. We compare measurements derived from XMM and Chandra observations with one another and compare both to gravitational lensing mass profiles derived with CLASH HST and ground-based lensing data. Radial profiles of Chandra and XMM electron density and enclosed gas mass are nearly identical, indicating that differences in hydrostatic masses inferred from X-ray observations arise from differences in Tx measurements. Encouragingly, cluster Txs are consistent with one another at ~100-200 kpc radii but XMM Tx systematically decline relative to Chandra Tx at larger radii. The angular dependence of the discrepancy suggests additional investigation on systematics such as the XMM point spread function correction, vignetting and off-axis responses. We present the CLASH-X mass-profile comparisons in the form of cosmology-independent and redshift-independent circular-velocity profiles. Ratios of Chandra HSE mass profiles to CLASH lensing profiles show no obvious radial dependence in the 0.3-0.8 Mpc range. However, the mean mass biases inferred from the WL and SaWLens data are different. e.g., the weighted-mean value at 0.5 Mpc is <b> = 0.12 for the WL comparison and <b> = -0.11 for the SaWLens comparison. The ratios of XMM HSE mass profiles to CLASH lensing profiles show a pronounced radial dependence in the 0.3-1.0 Mpc range, with a weighted mean mass bias of value rising to <b>~0.3 at ~1 Mpc for the WL comparison and <b> of 0.25 for SaWLens comparison. The enclosed gas mass profiles from both Chandra and XMM rise to a value 1/8 times the total-mass profiles inferred from lensing at 0.5 Mpc and remain constant outside of that radius, suggesting that [8xMgas] profiles may be an excellent proxy for total-mass profiles at >0.5 Mpc in massive galaxy clusters.

preprint2015arXiv

Data augmentation for machine learning redshifts applied to SDSS galaxies

We present analyses of data augmentation for machine learning redshift estimation. Data augmentation makes a training sample more closely resemble a test sample, if the two base samples differ, in order to improve measured statistics of the test sample. We perform two sets of analyses by selecting 800k (1.7M) SDSS DR8 (DR10) galaxies with spectroscopic redshifts. We construct a base training set by imposing an artificial r band apparent magnitude cut to select only bright galaxies and then augment this base training set by using simulations and by applying the K-correct package to artificially place training set galaxies at a higher redshift. We obtain redshift estimates for the remaining faint galaxy sample, which are not used during training. We find that data augmentation reduces the error on the recovered redshifts by 40% in both sets of analyses, when compared to the difference in error between the ideal case and the non augmented case. The outlier fraction is also reduced by at least 10% and up to 80% using data augmentation. We finally quantify how the recovered redshifts degrade as one probes to deeper magnitudes past the artificial magnitude limit of the bright training sample. We find that at all apparent magnitudes explored, the use of data augmentation with tree based methods provide a estimate of the galaxy redshift with a negligible bias, although the error on the recovered values increases as we probe to deeper magnitudes. These results have applications for surveys which have a spectroscopic training set which forms a biased sample of all photometric galaxies, for example if the spectroscopic detection magnitude limit is shallower than the photometric limit.

preprint2015arXiv

Environment-based selection effects of Planck clusters

We investigate whether the large scale structure environment of galaxy clusters imprints a selection bias on Sunyaev Zel'dovich (SZ) catalogs. Such a selection effect might be caused by line of sight (LoS) structures that add to the SZ signal or contain point sources that disturb the signal extraction in the SZ survey. We use the Planck PSZ1 union catalog (Planck Collab- oration et al. 2013a) in the SDSS region as our sample of SZ selected clusters. We calculate the angular two-point correlation function (2pcf) for physically correlated, foreground and background structure in the RedMaPPer SDSS DR8 catalog with respect to each cluster. We compare our results with an optically selected comparison cluster sample and with theoretical predictions. In contrast to the hypothesis of no environment-based selection, we find a mean 2pcf for background structures of -0.049 on scales of $\lesssim 40'$, significantly non-zero at $\sim 4 σ$, which means that Planck clusters are more likely to be detected in regions of low background density. We hypothesize this effect arises either from background estimation in the SZ survey or from radio sources in the background. We estimate the defect in SZ signal caused by this effect to be negligibly small, of the order of $\sim 10^{-4}$ of the signal of a typical Planck detection. Analogously, there are no implications on X-ray mass measurements. However, the environ- mental dependence has important consequences for weak lensing follow up of Planck galaxy clusters: we predict that projection effects account for half of the mass contained within a 15' radius of Planck galaxy clusters. We did not detect a background underdensity of CMASS LRGs, which also leaves a spatially varying redshift dependence of the Planck SZ selection function as a possible cause for our findings.

preprint2015arXiv

Feature importance for machine learning redshifts applied to SDSS galaxies

We present an analysis of importance feature selection applied to photometric redshift estimation using the machine learning architecture Decision Trees with the ensemble learning routine Adaboost (hereafter RDF). We select a list of 85 easily measured (or derived) photometric quantities (or `features') and spectroscopic redshifts for almost two million galaxies from the Sloan Digital Sky Survey Data Release 10. After identifying which features have the most predictive power, we use standard artificial Neural Networks (aNN) to show that the addition of these features, in combination with the standard magnitudes and colours, improves the machine learning redshift estimate by 18% and decreases the catastrophic outlier rate by 32%. We further compare the redshift estimate using RDF with those from two different aNNs, and with photometric redshifts available from the SDSS. We find that the RDF requires orders of magnitude less computation time than the aNNs to obtain a machine learning redshift while reducing both the catastrophic outlier rate by up to 43%, and the redshift error by up to 25%. When compared to the SDSS photometric redshifts, the RDF machine learning redshifts both decreases the standard deviation of residuals scaled by 1/(1+z) by 36% from 0.066 to 0.041, and decreases the fraction of catastrophic outliers by 57% from 2.32% to 0.99%.

preprint2015arXiv

The Detection and Statistics of Giant Arcs Behind CLASH Clusters

We developed an algorithm to find and characterize gravitationally lensed galaxies (arcs) to perform a comparison of the observed and simulated arc abundance. Observations are from the Cluster Lensing And Supernova survey with Hubble (CLASH). Simulated CLASH images are created using the MOKA package and also clusters selected from the high resolution, hydrodynamical simulations, MUSIC, over the same mass and redshift range as the CLASH sample. The algorithm' s arc elongation accuracy, completeness and false positive rate are determined and used to compute an estimate of the true arc abundance. We derive a lensing efficiency of $4 \pm 1$ arcs (with length $\ge 6"$ and length-to-width ratio $\ge 7$) per cluster for the X-ray selected CLASH sample, $4 \pm 1$ arcs per cluster for the MOKA simulated sample and $3 \pm 1$ arcs per cluster for the MUSIC simulated sample. The observed and simulated arc statistics are in full agreement. We measure the photometric redshifts of all detected arcs and find a median redshift $z_s = 1.9$ with 33% of the detected arcs having $z_s > 3$. We find that the arc abundance does not depend strongly on the source redshift distribution but is sensitive to the mass distribution of the dark matter halos (e.g. the $c-M$ relation). Our results show that consistency between the observed and simulated distributions of lensed arc sizes and axial ratios can be achieved by using cluster-lensing simulations that are carefully matched to the selection criteria used in the observations.

preprint2014arXiv

A Consistent Study of Metallicity Evolution at 0.8 < z < 2.6

We present the correlations between stellar mass, star formation rate (SFR) and [NII]/Ha flux ratio as indicator of gas-phase metallicity for a sample of 222 galaxies at 0.8 < z < 2.6 and log(M*/Msun)=9.0-11.5 from the LUCI, SINS/zC-SINF and KMOS3D surveys. This sample provides a unique analysis of the mass-metallicity relation (MZR) over an extended redshift range using consistent data analysis techniques and strong-line metallicity indicator. We find a constant slope at the low-mass end of the relation and can fully describe its redshift evolution through the evolution of the characteristic turnover mass where the relation begins to flatten at the asymptotic metallicity. At fixed mass and redshift, our data do not show a correlation between the [NII]/Ha ratio and SFR, which disagrees with the 0.2-0.3dex offset in [NII]/Ha predicted by the "fundamental relation" between stellar mass, SFR and metallicity discussed in recent literature. However, the overall evolution towards lower [NII]/Ha at earlier times does broadly agree with these predictions.

preprint2014arXiv

CLASH-VLT: Constraints on the Dark Matter Equation of State from Accurate Measurements of Galaxy Cluster Mass Profiles

A pressureless scenario for the Dark Matter (DM) fluid is a widely adopted hypothesis, despite the absence of a direct observational evidence. According to general relativity, the total mass-energy content of a system shapes the gravitational potential well, but different test particles perceive this potential in different ways depending on their properties. Cluster galaxy velocities, being $\ll$c, depend solely on the gravitational potential, whereas photon trajectories reflect the contributions from the gravitational potential plus a relativistic-pressure term that depends on the cluster mass. We exploit this phenomenon to constrain the Equation of State (EoS) parameter of the fluid, primarily DM, contained in galaxy clusters. We use the complementary information provided by the kinematic and lensing mass profiles of the galaxy cluster MACS 1206.2-0847 at $z=0.44$, as obtained in an extensive imaging and spectroscopic campaign within the CLASH survey. The unprecedented high quality of our data-set and the properties of this cluster are well suited to determine the EoS parameter of the cluster fluid. Since baryons contribute at most $15\%$ to the total mass in clusters and their pressure is negligible, the EoS parameter we derive describes the behavior of the DM fluid. We obtain the most stringent constraint on the DM EoS parameter to date, $w=(p_r+2\,p_t)/(3\,c^2ρ)=0.00\pm0.15\mathrm{(stat)}\pm0.08\mathrm{(syst)}$, averaged over the radial range $0.5\,\mathrm{Mpc}\leq$$r$$\leq$$r_{200}$, where $p_r$ and $p_t$ are the radial and tangential pressure, and $ρ$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program.

preprint2014arXiv

CLASH: Extreme Emission Line Galaxies and Their Implication on Selection of High-Redshift Galaxies

We utilize the CLASH (Cluster Lensing And Supernova survey with Hubble) observations of 25 clusters to search for extreme emission-line galaxies (EELGs). The selections are carried out in two central bands: F105W (Y105) and F125W (J125), as the flux of the central bands could be enhanced by the presence of [O III] 4959, 5007 at redshift of about 0.93-1.14 and 1.57-1.79, respectively. The multi-band observations help to constrain the equivalent widths of emission lines. Thanks to cluster lensing, we are able to identify 52 candidates down to an intrinsic limiting magnitude of 28.5 and to a rest-frame [O III] 4959,5007 equivalent width of about 3737 angstrom. Our samples include a number of EELGs at lower luminosities that are missed in other surveys, and the extremely high equivalent width can be only found in such faint galaxies. These EELGs can mimic the dropout feature similar to that of high redshift galaxies and contaminate the color-color selection of high redshift galaxies when the S/N ratio is limited or the band coverage is incomplete. We predict that the fraction of EELGs in the future high redshift galaxy selections cannot be neglected.

preprint2014arXiv

CLASH: Weak-Lensing Shear-and-Magnification Analysis of 20 Galaxy Clusters

We present a joint shear-and-magnification weak-lensing analysis of a sample of 16 X-ray-regular and 4 high-magnification galaxy clusters at 0.19<z<0.69 selected from the Cluster Lensing And Supernova survey with Hubble (CLASH). Our analysis uses wide-field multi-color imaging, taken primarily with Suprime-Cam on the Subaru Telescope. From a stacked shear-only analysis of the X-ray-selected subsample, we detect the ensemble-averaged lensing signal with a total signal-to-noise ratio of ~25 in the radial range of 200 to 3500kpc/h. The stacked tangential-shear signal is well described by a family of standard density profiles predicted for dark-matter-dominated halos in gravitational equilibrium, namely the Navarro-Frenk-White (NFW), truncated variants of NFW, and Einasto models. For the NFW model, we measure a mean concentration of $c_{200c}=4.01^{+0.35}_{-0.32}$ at $M_{200c}=1.34^{+0.10}_{-0.09} 10^{15}M_{\odot}$. We show this is in excellent agreement with Lambda cold-dark-matter (LCDM) predictions when the CLASH X-ray selection function and projection effects are taken into account. The best-fit Einasto shape parameter is $α_E=0.191^{+0.071}_{-0.068}$, which is consistent with the NFW-equivalent Einasto parameter of $\sim 0.18$. We reconstruct projected mass density profiles of all CLASH clusters from a joint likelihood analysis of shear-and-magnification data, and measure cluster masses at several characteristic radii. We also derive an ensemble-averaged total projected mass profile of the X-ray-selected subsample by stacking their individual mass profiles. The stacked total mass profile, constrained by the shear+magnification data, is shown to be consistent with our shear-based halo-model predictions including the effects of surrounding large-scale structure as a two-halo term, establishing further consistency in the context of the LCDM model.

preprint2014arXiv

The 64 Mpixel wide field imager for the Wendelstein 2m Telescope: Design and Calibration

The Wendelstein Observatory of Ludwig Maximilians University of Munich has recently been upgraded with a modern 2m robotic telescope. One Nasmyth port of the telescope has been equipped with a wide-field corrector which preserves the excellent image quality (< 0.8" median seeing) of the site (Hopp et al. 2008) over a field of view of 0.7 degrees diameter. The available field is imaged by an optical imager (WWFI, the Wendelstein Wide Field Imager) built around a customized 2 $\times$ 2 mosaic of 4k $\times$ 4k 15 μm e2v CCDs from Spectral Instruments. This paper provides an overview of the design and the WWFI's performance. We summarize the system mechanics (including a structural analysis), the electronics (and its electromagnetic interference (EMI) protection) and the control software. We discuss in detail detector system parameters, i.e. gain and readout noise, quantum efficiency as well as charge transfer efficiency (CTE) and persistent charges. First on sky tests yield overall good predictability of system throughput based on lab measurements.

preprint2014arXiv

The M31 Near-Infrared Period-Luminosity Relation and its non-linearity for $δ$ Cep Variables with $0.5 \leq \log(P) \leq 1.7$

We present the largest M31 near-infrared (F110W (close to J band), F160W (H band)) Cepheid sample so far. The sample consists of 371 Cepheids with photometry obtained from the HST PHAT program. The sample of 319 fundamental mode Cepheids, 16 first overtone Cepheids and 36 type II Cepheids, was identified using the median absolute deviation (MAD) outlier rejection method we develop here. This method does not rely on priors and allows us to obtain this clean Cepheid sample without rejecting a large fraction of Cepheids. The obtained Period-Luminosity relations (PLRs) have a very small dispersion, i.e. 0.155 mag in F160W, despite using random phased observations. This remarkably small dispersion allows us to determine that the PLRs are significantly better described by a broken slope at ten days than a linear slope. The use of our sample as an anchor to determine the Hubble constant gives a $3.2\%$ larger Hubble constant compared to the Riess et al. (2012) sample.

preprint2014arXiv

Three Gravitationally Lensed Supernovae Behind CLASH Galaxy Clusters

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CLO12Car (z = 1.28), SN CLN12Did (z = 0.85), and SN CLA11Tib (z = 1.14), are located behind three different clusters, MACSJ1720.2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and Abell 383 (z = 0.187), respectively. Each SN was detected in Hubble Space Telescope (HST) optical and infrared images. Based on photometric classification, we find that SNe CLO12Car and CLN12Did are likely to be Type Ia supernovae (SNe Ia), while the classification of SN CLA11Tib is inconclusive. Using multi-color light-curve fits to determine a standardized SN Ia luminosity distance, we infer that SN CLO12Car was approximately 1.0 +/- 0.2 mag brighter than field SNe Ia at a similar redshift and ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is approximately 0.2 +/- 0.2 mag brighter than field SNe Ia. We derive independent estimates of the predicted magnification from CLASH strong+weak lensing maps of the clusters: 0.83 +/- 0.16 mag for SN CLO12Car, 0.28 +/- 0.08 mag for SN CLN12Did, and 0.43 +/- 0.11 mag for SN CLA11Tib. The two SNe Ia provide a new test of the cluster lens model predictions: we find that the magnifications based on the SN Ia brightness and those predicted by the lens maps are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe Ia (from the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct determination of absolute magnifications.

preprint2013arXiv

CLASH: Complete Lensing Analysis of the Largest Cosmic Lens MACS J0717.5+3745 and Surrounding Structures

The galaxy cluster MACS J0717.5+3745 (z=0.55) is the largest known cosmic lens, with complex internal structures seen in deep X-ray, Sunyaev-Zel'dovich effect and dynamical observations. We perform a combined weak and strong lensing analysis with wide-field BVRi'z' Subaru/Suprime-Cam observations and 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). We find consistent weak distortion and magnification measurements of background galaxies, and combine these signals to construct an optimally estimated radial mass profile of the cluster and its surrounding large-scale structure out to 5 Mpc/h. We find consistency between strong-lensing and weak-lensing in the region where these independent data overlap, <500 kpc/h. The two-dimensional weak-lensing map reveals a clear filamentary structure traced by distinct mass halos. We model the lensing shear field with 9 halos, including the main cluster, corresponding to mass peaks detected above 2.5σ_κ. The total mass of the cluster as determined by the different methods is M_{vir}=(2.8\pm0.4) \times 10^15 M_sun. Although this is the most massive cluster known at z>0.5, in terms of extreme value statistics we conclude that the mass of MACS J0717.5+3745 by itself is not in serious tension with LambdaCDM, representing only a ~2σ departure above the maximum simulated halo mass at this redshift.

preprint2013arXiv

Galaxy halo truncation and Giant Arc Surface Brightness Reconstruction in the Cluster MACSJ1206.2-0847

In this work we analyze the mass distribution of MACSJ1206.2-0847, especially focusing on the halo properties of its cluster members. The cluster appears relaxed in its X-ray emission, but has significant amounts of intracluster light which is not centrally concentrated, suggesting that galaxy-scale interactions are still ongoing despite the overall relaxed state. The cluster lenses 12 background galaxies into multiple images and one galaxy at $z=1.033$ into a giant arc and its counterimage. The multiple image positions and the surface brightness distribution (SFB) of the arc which is bent around several cluster members are sensitive to the cluster galaxy halo properties. We model the cluster mass distribution with a NFW profile and the galaxy halos with two parameters for the mass normalization and extent of a reference halo assuming scalings with their observed NIR--light. We match the multiple image positions at an r.m.s. level of $0.85\arcsec$ and can reconstruct the SFB distribution of the arc in several filters to a remarkable accuracy based on this cluster model. The length scale where the enclosed galaxy halo mass is best constrained is about 5 effective radii -- a scale in between those accessible to dynamical and field strong lensing mass estimates on one hand and galaxy--galaxy weak lensing results on the other hand. The velocity dispersion and halo size of a galaxy with $m_{\rm 160W,AB}=19.2$ or $M_{\rm B,Vega}=-20.7$ are $σ=150 \rm kms^{-1}$ and $r\approx 26\pm 6 \rm kpc$, indicating that the halos of the cluster galaxies are tidally stripped. We also reconstruct the unlensed source (which is smaller by a factor of $\sim5.8$ in area), demonstrating the increase of morphological information due to lensing and conclude that this galaxy has likely star--forming spiral arms with a red (older) central component.

preprint2013arXiv

Photometric Redshifts and Systematic Variations in the SEDs of Luminous Red Galaxies from the SDSS DR7

We describe the construction of a template set of spectral energy distributions (SEDs) for the estimation of photometric redshifts of luminous red galaxies (LRGs) with a Bayesian template fitting method. By examining the color properties of several publicly available SED sets within a redshift range of 0<z<0.5 and comparing them to SDSS DR7 data, we show that only some of the investigated SEDs approximately match the colors of the LRG data throughout the redshift range, however not at the quantitative level required for precise photometric redshifts. We generate new SEDs by superposing model SEDs of composite stellar populations with a burst model, allowing both components to be reddened by dust, in order to match the data in five different redshift bins. We select a set of SEDs which represents the LRG data in color space within five redshift bins, thus defining our new SED template set for photometric redshift estimates. The results we get with the new template set and our Bayesian template fitting photometric redshift code (PhotoZ) are nearly unbiased, with a scatter of σ(Δz)=0.027 (including outliers), and a fraction of catastrophic outliers (|z_phot-z_spec|/(1+z_spec)>0.15) of 0.12%. We show that templates that optimally describe the brightest galaxies (-24.5<M_R<-22.7) indeed vary from z=0.1 to z=0.5, consistent with aging of the stellar population. Furthermore, we find that templates that optimally describe galaxies at z<0.1 strongly differ as a function of the absolute magnitude of the galaxies, indicating an increase in star formation activity for less luminous galaxies. Our findings based on the photometry of the SDSS LRGs and our SED template fitting are supported by comparison to the average SDSS LRG spectra in different luminosity and redshift bins.

preprint2013arXiv

Properties of M31. II: A Cepheid disk sample derived from the first year of PS1 PAndromeda data

We present a sample of Cepheid variable stars towards M31 based on the first year of regular M31 observations of the PS1 survey in the r_P1 and i_P1 filters. We describe the selection procedure for Cepheid variable stars from the overall variable source sample and develop an automatic classification scheme using Fourier decomposition and the location of the instability strip. We find 1440 fundamental mode (classical δ) Cep stars, 126 Cepheids in the first overtone mode, and 147 belonging to the Population II types. 296 Cepheids could not be assigned to one of these classes and 354 Cepheids were found in other surveys. These 2009 Cepheids constitute the largest Cepheid sample in M31 known so far and the full catalog is presented in this paper. We briefly describe the properties of our sample in its spatial distribution throughout the M31 galaxy, in its age properties, and we derive an apparent period-luminosity relation (PLR) in our two bands. The Population I Cepheids nicely follow the dust pattern of the M31 disk, whereas the 147 Type II Cepheids are distributed throughout the halo of M31. We outline the time evolution of the star formation in the major ring found previously and find an age gradient. A comparison of our PLR to previous results indicates a curvature term in the PLR.

preprint2013arXiv

The Contribution of Halos with Different Mass Ratios to the Overall Growth of Cluster-Sized Halos

We provide a new observational test for a key prediction of the ΛCDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $σ$ agreement with ΛCDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated.

preprint2012arXiv

A highly magnified candidate for a young galaxy seen when the Universe was 500 Myrs old

The early Universe at redshift z\sim6-11 marks the reionization of the intergalactic medium, following the formation of the first generation of stars. However, those young galaxies at a cosmic age of \lesssim 500 million years (Myr, at z \gtrsim 10) remain largely unexplored as they are at or beyond the sensitivity limits of current large telescopes. Gravitational lensing by galaxy clusters enables the detection of high-redshift galaxies that are fainter than what otherwise could be found in the deepest images of the sky. We report the discovery of an object found in the multi-band observations of the cluster MACS1149+22 that has a high probability of being a gravitationally magnified object from the early universe. The object is firmly detected (12 sigma) in the two reddest bands of HST/WFC3, and not detected below 1.2 μm, matching the characteristics of z\sim9 objects. We derive a robust photometric redshift of z = 9.6 \pm 0.2, corresponding to a cosmic age of 490 \pm 15Myr (i.e., 3.6% of the age of the Universe). The large number of bands used to derive the redshift estimate make it one of the most accurate estimates ever obtained for such a distant object. The significant magnification by cluster lensing (a factor of \sim15) allows us to analyze the object's ultra-violet and optical luminosity in its rest-frame, thus enabling us to constrain on its stellar mass, star-formation rate and age. If the galaxy is indeed at such a large redshift, then its age is less than 200 Myr (at the 95% confidence level), implying a formation redshift of zf \lesssim 14. The object is the first z>9 candidate that is bright enough for detailed spectroscopic studies with JWST, demonstrating the unique potential of galaxy cluster fields for finding highly magnified, intrinsically faint galaxies at the highest redshifts.

preprint2012arXiv

Bright radio emission from an ultraluminous stellar-mass microquasar in M31

A subset of ultraluminous X-ray sources (those with luminosities < 10^40 erg/s) are thought to be powered by the accretion of gas onto black holes with masses of ~5-20 M_solar, probably via an accretion disc. The X-ray and radio emission are coupled in such Galactic sources, with the radio emission originating in a relativistic jet thought to be launched from the innermost regions near the black hole, with the most powerful emission occurring when the rate of infalling matter approaches a theoretical maximum (the Eddington limit). Only four such maximal sources are known in the Milky Way, and the absorption of soft X-rays in the interstellar medium precludes determining the causal sequence of events that leads to the ejection of the jet. Here we report radio and X-ray observations of a bright new X-ray source whose peak luminosity can exceed 10^39 erg/s in the nearby galaxy, M31. The radio luminosity is extremely high and shows variability on a timescale of tens of minutes, arguing that the source is highly compact and powered by accretion close to the Eddington limit onto a stellar mass black hole. Continued radio and X-ray monitoring of such sources should reveal the causal relationship between the accretion flow and the powerful jet emission.

preprint2012arXiv

CLASH: Mass Distribution in and around MACS J1206.2-0847 from a Full Cluster Lensing Analysis

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dlnΣ/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

preprint2012arXiv

CLASH: Precise New Constraints on the Mass Profile of Abell 2261

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with ΛCDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent ΛCDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

preprint2012arXiv

CLASH: Three Strongly Lensed Images of a Candidate z ~ 11 Galaxy

We present a candidate for the most distant galaxy known to date with a photometric redshift z = 10.7 +0.6 / -0.4 (95% confidence limits; with z < 9.5 galaxies of known types ruled out at 7.2-sigma). This J-dropout Lyman Break Galaxy, named MACS0647-JD, was discovered as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). We observe three magnified images of this galaxy due to strong gravitational lensing by the galaxy cluster MACSJ0647.7+7015 at z = 0.591. The images are magnified by factors of ~8, 7, and 2, with the brighter two observed at ~26th magnitude AB (~0.15 uJy) in the WFC3/IR F160W filter (~1.4 - 1.7 um) where they are detected at >~ 12-sigma. All three images are also confidently detected at >~ 6-sigma in F140W (~1.2 - 1.6 um), dropping out of detection from 15 lower wavelength HST filters (~0.2 - 1.4 um), and lacking bright detections in Spitzer/IRAC 3.6um and 4.5um imaging (~3.2 - 5.0 um). We rule out a broad range of possible lower redshift interlopers, including some previously published as high redshift candidates. Our high redshift conclusion is more conservative than if we had neglected a Bayesian photometric redshift prior. Given CLASH observations of 17 high mass clusters to date, our discoveries of MACS0647-JD at z ~ 10.8 and MACS1149-JD1 at z ~ 9.6 are consistent with a lensed luminosity function extrapolated from lower redshifts. This would suggest that low luminosity galaxies could have reionized the universe. However given the significant uncertainties based on only two galaxies, we cannot yet rule out the sharp drop off in number counts at z >~ 10 suggested by field searches.

preprint2012arXiv

Cluster Lensing And Supernova survey with Hubble (CLASH): An Overview

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated.

preprint2012arXiv

Golden gravitational lensing systems from the Sloan Lens ACS Survey. II. SDSS J1430+4105: A precise inner total mass profile from lensing alone

We study the SLACS strong lensing system SDSSJ1430+4105 at z=0.285. The lensed source (z=0.575) of this system has a complex morphology with several subcomponents. Its subcomponents span a radial range from 4 to 10kpc in the lens plane. Therefore we can constrain the slope of the total projected mass profile around the Einstein radius (R_E) from lensing alone. We measure a density profile that is slightly but not significantly shallower than isothermal at R_E. We decompose the mass of the lensing galaxy into a de Vaucouleurs (deV) component to trace the stars and an additional dark component. The spread of multiple image components over a large radial range also allows us to determine the amplitude of the deV and dark matter components separately. We get a mass to light ratio of M_deV/L_B~5.5\pm1.5M/L_sun,B and a dark matter fraction within R_E of ~20 to 40%. Modelling the star formation history assuming composite stellar populations at solar metallicity to the galaxy's photometry yields a mass to light ratio of M_star,salp/L_B~4.0_{-1.3}^{+0.6}M/L_sun,B and M_star,chab/L_B~2.3_{-0.8}^{+0.3}M/L_sun,B for Salpeter and Chabrier IMFs, respectively. Hence, the mass to light ratio derived from lensing is more Salpeter-like, in agreement with results for massive Coma galaxies and other nearby massive early type galaxies. We examine the consequences of the galaxy group in which the lensing galaxy is embedded, showing that it has little influence on the mass to light ratio obtained for the deV component of the lensing galaxy. Finally, we decompose the projected, azimuthally averaged 2D density distribution of the deV and dark matter component of the lensing signal into spherically averaged 3D density profiles. We can show that the 3D dark and luminous matter density within R_E~0.6R_eff of this SLACS galaxy is similar to the values of Coma galaxies with the same velocity dispersions.

preprint2011arXiv

Measuring Lensing Magnification of Quasars by Large Scale Structure using the Variability-Luminosity Relation

We introduce a technique to measure gravitational lensing magnification using the variability of type I quasars. Quasars' variability amplitudes and luminosities are tightly correlated, on average. Magnification due to gravitational lensing increases the quasars' apparent luminosity, while leaving the variability amplitude unchanged. Therefore, the mean magnification of an ensemble of quasars can be measured through the mean shift in the variability-luminosity relation. As a proof of principle, we use this technique to measure the magnification of quasars spectroscopically identified in the Sloan Digital Sky Survey, due to gravitational lensing by galaxy clusters in the SDSS MaxBCG catalog. The Palomar-QUEST Variability Survey, reduced using the DeepSky pipeline, provides variability data for the sources. We measure the average quasar magnification as a function of scaled distance (r/R200) from the nearest cluster; our measurements are consistent with expectations assuming NFW cluster profiles, particularly after accounting for the known uncertainty in the clusters' centers. Variability-based lensing measurements are a valuable complement to shape-based techniques because their systematic errors are very different, and also because the variability measurements are amenable to photometric errors of a few percent and to depths seen in current wide-field surveys. Given the data volume expected from current and upcoming surveys, this new technique has the potential to be competitive with weak lensing shear measurements of large scale structure.

preprint2011arXiv

The BOSS Emission-Line Lens Survey (BELLS). I. A large spectroscopically selected sample of Lens Galaxies at redshift ~ 0.5

We present a catalog of 25 definite and 11 probable strong galaxy-galaxy gravitational lens systems with redshifts 0.4 \lesssim z \lesssim 0.7, discovered spectroscopically by the presence of higher redshift emission-lines within the Baryon Oscillation Spectroscopic Survey (BOSS) of luminous galaxies, and confirmed with high-resolution Hubble Space Telescope (HST) images of 44 candidates. Our survey extends the methodology of the Sloan Lens ACS Survey (SLACS: Bolton et al. 2006; 2008) to higher redshift. We describe the details of the BOSS spectroscopic candidate detections, our HST Adanced Camera for Surveys (ACS) image processing and analysis methods, and our strong gravitational lens modeling procedure. We report BOSS spectroscopic parameters and ACS photometric parameters for all candidates, and mass-distribution parameters for the best-fit singular isothermal ellipsoid models of definite lenses. Our sample to date was selected using only the first six months of BOSS survey-quality spectroscopic data. The full 5-year BOSS database should produce a sample of several hundred strong galaxy-galaxy lenses and in combination with SLACS lenses at lower redshift, strongly constrain the redshift evolution of the structure of elliptical, bulge-dominated galaxies as a function of luminosity, stellar mass, and rest-frame color, thereby providing a powerful test for competing theories of galaxy formation and evolution.

preprint2010arXiv

A pan-chromatic view of the galaxy cluster XMMU J1230.3+1339 at z=0.975 - Observing the assembly of a massive system

We present a comprehensive galaxy cluster study of XMMU J1230.3+1339 based on a joint analysis of X-ray data, optical imaging and spectroscopy observations, weak lensing results, and radio properties for achieving a detailed multi-component view of this newly discovered system at z=0.975. We find an optically very rich and massive system with M200$\simeq$(4.2$\pm$0.8)$\times$10^14 M$\sun$, Tx$\simeq$5.3(+0.7--0.6)keV, and Lx$\simeq$(6.5$\pm$0.7)$\times$10^44 erg/s, for which various widely used mass proxies are measured and compared. We have identified multiple cluster-related components including a central fly-through group close to core passage with associated marginally extended 1.4GHz radio emission possibly originating from the turbulent wake region of the merging event. On the cluster outskirts we see evidence for an on-axis infalling group with a second Brightest Cluster Galaxy (BCG) and indications for an additional off-axis group accretion event. We trace two galaxy filaments beyond the nominal cluster radius and provide a tentative reconstruction of the 3D-accretion geometry of the system. In terms of total mass, ICM structure, optical richness, and the presence of two dominant BCG-type galaxies, the newly confirmed cluster XMMU J1230.3+1339 is likely the progenitor of a system very similar to the local Coma cluster, differing by 7.6 Gyr of structure evolution.

preprint1994arXiv

Light Propagation in Arbitrary Spacetimes and the Gravitational Lens Approximation

We reformulate the transport equation which determines the size, shape and orientation of infinitesimal light beams in arbitrary spacetimes. The behaviour of such light beams near vertices and conjugate points is investigated, with special attention to the singular behaviour of the optical scalars. We then specialize the general transport equation to the case of an approximate metric of an inhomogeneous universe, which is a Friedmann metric `on average' with superposed isolated weak matter inhomogeneities. In a series of well-defined approximations, the equations of gravitational lens theory are derived. Finally, we derive a relative optical focusing equation which describes the focusing of light beams relative to the case that the beam is unaffected by matter inhomogeneities in the universe, from which it follows immediately that no beam can be focused less than one which is unaffected by matter clumps, before it propagates through its first conjugate point.