Search

We explore the universality of neural encodings in convolutional neural networks trained on image classification tasks. We develop a procedure to directly compare the learned weights rather than their representations. It is based on a factorization of spatial and channel dimensions and measures the similarity of aligned weight covariances. We show that, for a range of layers of VGG-type networks, the learned eigenvectors appear to be universal across different natural image datasets. Our results suggest the existence of a universal neural encoding for natural images. They explain, at a more fundamental level, the success of transfer learning. Our work shows that, instead of aiming at maximizing the performance of neural networks, one can alternatively attempt to maximize the universality of the learned encoding, in order to build a principled foundation model., Comment: Appeared at the ICLR 2024 Workshop on Representational Alignment (Re-Align), 13 pages, 5 figures

How does the human visual cortex encode sensory information? To address this question, we explore the covariance structure of neural representations. We perform a cross-decomposition analysis of fMRI responses to natural images in multiple individuals from the Natural Scenes Dataset and find that neural representations systematically exhibit a power-law covariance spectrum over four orders of magnitude in ranks. This scale-free structure is found in multiple regions along the visual hierarchy, pointing to the existence of a generic encoding strategy in visual cortex. We also show that, up to a rotation, a large ensemble of principal axes of these population codes are shared across subjects, showing the existence of a universal high-dimensional representation. This suggests a high level of convergence in how the human brain learns to represent natural scenes despite individual differences in neuroanatomy and experience. We further demonstrate that a spectral approach is critical for characterizing population codes in their full extent, and in doing so, we reveal a vast space of uncharted dimensions that have been out of reach for conventional variance-weighted methods. A global view of neural representations thus requires embracing their high-dimensional nature and understanding them statistically rather than through visual or semantic interpretation of individual dimensions., Comment: 22 pages, 4 main figures, 7 supplementary figures

As weak lensing surveys go deeper, there is an increasing need for reliable characterization of non-Gaussian structures at small angular scales. Here we present the first cosmological constraints with weak lensing scattering transform, a statistical estimator that combines efficiency, robustness, and interpretability. With the Hyper Suprime-Cam survey (HSC) year 1 data, we obtain $\Omega_\text{m}=0.29_{-0.03}^{+0.04}$, $S_8\equiv \sigma_8(\Omega_\text{m}/0.3)^{0.5}=0.83\pm0.02$, and intrinsic alignment strength $A_\text{IA}=1.0\pm0.4$ through simulation-based forward modeling. Our constraints are consistent with those derived from Planck. The error bar of $\Omega_\text{m}$ is 2 times tighter than that obtained from the power spectrum when the same scale range is used. This constraining power is on par with that of convolutional neural networks, suggesting that further investment in spatial information extraction may not yield substantial benefits. We also point out an internal tension of $S_8$ estimates linked to a redshift bin around z ~ 1 in the HSC data. We found that discarding that bin leads to a consistent decrease of $S_8$ from 0.83 to 0.79, for all statistical estimators. We argue that photometric redshift estimation is now the main limitation in the estimation of $S_8$ using HSC. This limitation is likely to affect other ground-based weak lensing surveys reaching redshifts greater than one. Alternative redshift estimation techniques, like clustering redshifts, may help alleviate this limitation., Comment: 8 pages, 6 figures

Physicists routinely need probabilistic models for a number of tasks such as parameter inference or the generation of new realizations of a field. Establishing such models for highly non-Gaussian fields is a challenge, especially when the number of samples is limited. In this paper, we introduce scattering spectra models for stationary fields and we show that they provide accurate and robust statistical descriptions of a wide range of fields encountered in physics. These models are based on covariances of scattering coefficients, i.e. wavelet decomposition of a field coupled with a point-wise modulus. After introducing useful dimension reductions taking advantage of the regularity of a field under rotation and scaling, we validate these models on various multi-scale physical fields and demonstrate that they reproduce standard statistics, including spatial moments up to 4th order. These scattering spectra provide us with a low-dimensional structured representation that captures key properties encountered in a wide range of physical fields. These generic models can be used for data exploration, classification, parameter inference, symmetry detection, and component separation., Comment: 11 pages, 6 figures, plus appendices, updated to published version

A central question in deep learning is to understand the functions learned by deep networks. What is their approximation class? Do the learned weights and representations depend on initialization? Previous empirical work has evidenced that kernels defined by network activations are similar across initializations. For shallow networks, this has been theoretically studied with random feature models, but an extension to deep networks has remained elusive. Here, we provide a deep extension of such random feature models, which we call the rainbow model. We prove that rainbow networks define deterministic (hierarchical) kernels in the infinite-width limit. The resulting functions thus belong to a data-dependent RKHS which does not depend on the weight randomness. We also verify numerically our modeling assumptions on deep CNNs trained on image classification tasks, and show that the trained networks approximately satisfy the rainbow hypothesis. In particular, rainbow networks sampled from the corresponding random feature model achieve similar performance as the trained networks. Our results highlight the central role played by the covariances of network weights at each layer, which are observed to be low-rank as a result of feature learning., Comment: 59 pages, 10 figures. To appear at JMLR

Extracting information from stochastic fields or textures is a ubiquitous task in science, from exploratory data analysis to classification and parameter estimation. From physics to biology, it tends to be done either through a power spectrum analysis, which is often too limited, or the use of convolutional neural networks (CNNs), which require large training sets and lack interpretability. In this paper, we advocate for the use of the scattering transform (Mallat 2012), a powerful statistic which borrows mathematical ideas from CNNs but does not require any training, and is interpretable. We show that it provides a relatively compact set of summary statistics with visual interpretation and which carries most of the relevant information in a wide range of scientific applications. We present a non-technical introduction to this estimator and we argue that it can benefit data analysis, comparison to models and parameter inference in many fields of science. Interestingly, understanding the core operations of the scattering transform allows one to decipher many key aspects of the inner workings of CNNs., Comment: 18 pages, 16 figures

As weak lensing surveys become deeper, they reveal more non-Gaussian aspects of the convergence field which can only be extracted using statistics beyond the power spectrum. In Cheng et al. (2020) we showed that the scattering transform, a novel statistic borrowing mathematical concepts from convolutional neural networks, is a powerful tool for cosmological parameter estimation in the non-Gaussian regime. Here, we extend that analysis to explore its sensitivity to dark energy and neutrino mass parameters with weak lensing surveys. We first use image synthesis to show visually that, compared to the power spectrum and bispectrum, the scattering transform provides a better statistical vocabulary to characterize the perceptual properties of lensing mass maps. We then show that it is also better suited for parameter inference: (i) it provides higher sensitivity in the noiseless regime, and (ii) at the noise level of Rubin-like surveys, though the constraints are not significantly tighter than those of the bispectrum, the scattering coefficients have a more Gaussian sampling distribution, which is an important property for likelihood parametrization and accurate cosmological inference. We argue that the scattering coefficients are preferred statistics considering both constraining power and likelihood properties., Comment: 9 pages, 6 figures. Accepted to MNRAS

Scattering of seismic waves can reveal subsurface structures but usually in a piecemeal way focused on specific target areas. We used a manifold learning algorithm called "the Sequencer" to simultaneously analyze thousands of seismograms of waves diffracting along the core-mantle boundary and obtain a panoptic view of scattering across the Pacific region. In nearly half of the diffracting waveforms, we detected seismic waves scattered by three-dimensional structures near the core-mantle boundary. The prevalence of these scattered arrivals shows that the region hosts pervasive lateral heterogeneity. Our analysis revealed loud signals due to a plume root beneath Hawaii and a previously unrecognized ultralow-velocity zone beneath the Marquesas Islands. These observations illustrate how approaches flexible enough to detect robust patterns with little to no user supervision can reveal distinctive insights into the deep Earth., Comment: 13 pages, 4 figures. Supplement available at: http://www.geol.umd.edu/facilities/seismology/wp-content/uploads/2013/02/Kim_et_al_2020_SOM.pdf

As cosmic structures form, matter density fluctuations collapse gravitationally and baryonic matter is shock-heated and thermalized. We therefore expect a connection between the mean gravitational potential energy density of collapsed halos, $\Omega_{W}^{\rm halo}$, and the mean thermal energy density of baryons, $\Omega_{\rm th}$. These quantities can be obtained using two fundamentally different estimates: we compute $\Omega_{W}^{\rm halo}$ using the theoretical framework of the halo model which is driven by dark matter statistics, and measure $\Omega_{\rm th}$ using the Sunyaev-Zeldovich (SZ) effect which probes the mean thermal pressure of baryons. First, we derive that, at the present time, about 90% of $\Omega_{W}^{\rm halo}$ originates from massive halos with $M>10^{13}\,M_\odot$. Then, using our measurements of the SZ background, we find that $\Omega_{\rm th}$ accounts for about 80% of the kinetic energy of the baryons available for pressure in halos at $z\lesssim 0.5$. This constrains the amount of non-thermal pressure, e.g., due to bulk and turbulent gas motion sourced by mass accretion, to be about $\Omega_{\rm non-th}\simeq 0.4\times 10^{-8}$ at $z=0$., Comment: 11 pages + references, 4 figures, 2 tables. (v2) Expanded discussion on the modelling uncertainty. (v3) Fixed a minor typo in Eq.(22). Accepted for publication in the Astrophysical Journal. Julia codes to reproduce the figures and tables are available in https://github.com/komatsu5147/OmegaGrav.jl

The cosmic thermal history, quantified by the evolution of the mean thermal energy density in the universe, is driven by the growth of structures as baryons get shock heated in collapsing dark matter halos. This process can be probed by redshift-dependent amplitudes of the thermal Sunyaev-Zeldovich (SZ) effect background. To do so, we cross-correlate eight sky intensity maps in the $\it{Planck}$ and Infrared Astronomical Satellite missions with two million spectroscopic redshift references in the Sloan Digital Sky Surveys. This delivers snapshot spectra for the far-infrared to microwave background light as a function of redshift up to $z\sim3$. We decompose them into the SZ and thermal dust components. Our SZ measurements directly constrain $\langle bP_{\rm e} \rangle$, the halo bias-weighted mean electron pressure, up to $z\sim 1$. This is the highest redshift achieved to date, with uncorrelated redshift bins thanks to the spectroscopic references. We detect a threefold increase in the density-weighted mean electron temperature $\bar{T}_{\rm{e}}$ from $7\times 10^5~{\rm K}$ at $z=1$ to $2\times 10^6~{\rm K}$ today. Over $z=1$-$0$, we witness the build-up of nearly $70\%$ of the present-day mean thermal energy density $\rho_{\rm{th}}$, with the corresponding density parameter $\Omega_{\rm th}$ reaching $1.5 \times10^{-8}$. We find the mass bias parameter of $\it{Planck}$'s universal pressure profile of $B=1.27$ (or $1-b=1/B=0.79$), consistent with the magnitude of non-thermal pressure in gas motion and turbulence from mass assembly. We estimate the redshift-integrated mean Compton parameter $y\sim1.2\times10^{-6}$, which will be tested by future spectral distortion experiments. More than half of which originates from the large-scale structure at $z<1$, which we detect directly., Comment: 22 pages, 11 figures. Accepted for publication in ApJ

Scientists aim to extract simplicity from observations of the complex world. An important component of this process is the exploration of data in search of trends. In practice, however, this tends to be more of an art than a science. Among all trends existing in the natural world, one-dimensional trends, often called sequences, are of particular interest as they provide insights into simple phenomena. However, some are challenging to detect as they may be expressed in complex manners. We present the Sequencer, an algorithm designed to generically identify the main trend in a dataset. It does so by constructing graphs describing the similarities between pairs of observations, computed with a set of metrics and scales. Using the fact that continuous trends lead to more elongated graphs, the algorithm can identify which aspects of the data are relevant in establishing a global sequence. Such an approach can be used beyond the proposed algorithm and can optimize the parameters of any dimensionality reduction technique. We demonstrate the power of the Sequencer using real-world data from astronomy, geology as well as images from the natural world. We show that, in a number of cases, it outperforms the popular t-SNE and UMAP dimensionality reduction techniques. This approach to exploratory data analysis, which does not rely on training nor tuning of any parameter, has the potential to enable discoveries in a wide range of scientific domains. The source code is available on github and we provide an online interface at \url{http://sequencer.org}., Comment: Code available at: https://github.com/dalya/Sequencer. Comments are welcome!

Parameter estimation with non-Gaussian stochastic fields is a common challenge in astrophysics and cosmology. In this paper, we advocate performing this task using the scattering transform, a statistical tool sharing ideas with convolutional neural networks (CNNs) but requiring no training nor tuning. It generates a compact set of coefficients, which can be used as robust summary statistics for non-Gaussian information. It is especially suited for fields presenting localized structures and hierarchical clustering, such as the cosmological density field. To demonstrate its power, we apply this estimator to a cosmological parameter inference problem in the context of weak lensing. On simulated convergence maps with realistic noise, the scattering transform outperforms classic estimators and is on a par with state-of-the-art CNN. It retains the advantages of traditional statistical descriptors, has provable stability properties, allows to check for systematics, and importantly, the scattering coefficients are interpretable. It is a powerful and attractive estimator for observational cosmology and the study of physical fields in general., Comment: 13 pages, 7 figures; accepted to MNRAS

Double white dwarf (double-WD) binaries may merge within a Hubble time and produce high-mass WDs. Compared to other high-mass WDs, the double-WD merger products have higher velocity dispersion because they are older. With the power of Gaia data, we show strong evidence for double-WD merger products among high-mass WDs by analyzing the transverse-velocity distribution of more than a thousand high-mass WDs (0.8--1.3 $M_\odot$). We estimate that the fraction of double-WD merger products in our sample is about 20 %. We also obtain a precise double-WD merger rate and its mass dependence. Our merger rate estimates are close to binary population synthesis results and support the idea that double-WD mergers may contribute to a significant fraction of type Ia supernovae., Comment: 9 pages, 5 figures; published in ApJ

Recently, the power of Gaia data has revealed an enhancement of high-mass white dwarfs (WDs) on the Hertzsprung--Russell diagram, called the Q branch. This branch is located at the high-mass end of the recently identified crystallization branch. Investigating its properties, we find that the number density and velocity distribution on the Q branch cannot be explained by the cooling delay of crystallization alone, suggesting the existence of an extra cooling delay. To quantify this delay, we statistically compare two age indicators -- the dynamical age inferred from transverse velocity, and the photometric isochrone age -- for more than one thousand high-mass WDs (1.08--1.23 $M_\odot$) selected from Gaia Data Release 2. We show that about 6 % of the high-mass WDs must experience an 8 Gyr extra cooling delay on the Q branch, in addition to the crystallization and merger delays. This cooling anomaly is a challenge for WD cooling models. We point out that $^{22}$Ne settling in C/O-core WDs could account for this extra cooling delay., Comment: 18 pages, 11 figures; published in ApJ; discussions of double-WD merger rate have been moved to arxiv:1910.09558

The scaling relations between supermassive black holes and their host galaxy properties are of fundamental importance in the context black hole-host galaxy co-evolution throughout cosmic time. In this work, we use a novel algorithm that identifies smooth trends in complex datasets and apply it to a sample of 2,000 type I active galactic nuclei (AGN) spectra. We detect a sequence in emission line shapes and strengths which reveals a correlation between the narrow L([OIII])/L(H$\beta$) line ratio and the width of the broad H$\alpha$. This scaling relation ties the kinematics of the gas clouds in the broad line region to the ionisation state of the narrow line region, connecting the properties of gas clouds kiloparsecs away from the black hole to material gravitationally bound to it on sub-parsec scales. This relation can be used to estimate black hole masses from narrow emission lines only. It therefore enables black hole mass estimation for obscured type II AGN and allows us to explore the connection between black holes and host galaxy properties for thousands of objects, well beyond the local Universe. Using this technique, we present the $M_{\mathrm{BH}} - \sigma$ and $M_{\mathrm{BH}} - M_{*}$ scaling relations for a sample of about 10,000 type II AGN from SDSS. These relations are remarkably consistent with those observed for type I AGN, suggesting that this new method may perform as reliably as the classical estimate used in non-obscured type I AGN. These findings open a new window for studies of black hole-host galaxy co-evolution throughout cosmic time., Comment: to be submitted to MNRAS, comments are welcome!

Cosmic photons can be efficiently collected by broadband intensity mapping but information on their emission redshift and frequency is largely lost. We introduce a technique to statistically recover these otherwise collapsed dimensions by exploiting information in spatial fluctuations and apply it to the Galaxy Evolution Explorer (GALEX) All Sky and Medium Imaging Surveys. By spatially cross-correlating photons in the GALEX far-UV (1500A) and near-UV (2300A) bands with a million spectroscopic objects in the Sloan Digital Sky Survey as a function of redshift, we robustly detect the redshift-dependent intensity of the UV background (UVB) modulated by its clustering bias up to $z\sim2$. These measurements clearly reveal the imprints of UVB spectral features redshifting through the filters. Using a simple parameterization, we simultaneously fit a UVB emissivity and clustering bias factor to these observations and constrain the main spectral features of the UVB spectrum: (i) the Lyman break, (ii) the non-ionizing UV continuum, which agrees with the Haardt & Madau model but does not rely on any assumption regarding the nature of the sources, and (iii) the Ly$\alpha$ emission, whose luminosity density is consistent with estimates of the combined galaxy and AGN contributions at $z\sim1$. Because the technique probes the total background including low surface brightness emission, we place constraints on the amount of UV light originating from the diffuse intergalactic medium (IGM). Finally, the clustering bias of UV photons is found to be chromatic and evolving. Our frequency- and redshift-dependent UVB measurement delivers a summary statistic of the universe's net radiation output from stars, black holes, and the IGM combined., Comment: Accepted for publication in the Astrophysical Journal

Luminous red galaxies (LRGs) are the most massive galaxies at $z\sim 0.5$ and, by selection, have negligible star formation. These objects have halo masses between those of $L_{*}$ galaxies, whose circumgalactic media (CGM) are observed to have large masses of cold gas, and clusters of galaxies, which primarily contain hot gas. Here we report detections of strong and extended metal (CIII 977) and HI lines in the CGM of two LRGs. The CIII lines have equivalent widths of $\sim 1.8$ \r{A} and $\sim 1.2$ \r{A} , and velocity spreads of $\sim 796$ km s$^{-1}$ and $\sim 1245$ km s$^{-1}$, exceeding all such measurements from local $\sim L_{*}$ galaxies (maximal CIII equivalent widths $\sim 1$ \r{A}). The data demonstrate that a subset of halos hosting very massive, quenched galaxies contain significant complexes of cold gas. Possible scenarios to explain our observations include that the LRGs' CGM originate from past activity (e.g., star formation or active galactic nuclei driven outflows) or from the CGM of galaxies in overlapping subhalos. We favor the latter scenario, in which the properties of the CGM are more tightly linked to the underlying dark matter halo than properties of the targeted galaxies (e.g., star formation)., Comment: 9 pages, 5 figures. Accepted for publication in ApJ

Extragalactic astronomy relies on the accurate estimation of source photometry corrected for Milky Way dust extinction. This has motivated the creation of a number of "Galactic" dust maps. We investigate whether these maps are contaminated by extragalactic signals using the clustering-redshift technique, i.e., by measuring a set of angular cross-correlations with spectroscopic objects as a function of redshift. Our tomographic analysis reveals imprints of extragalactic large-scale structure patterns in nine out of 10 Galactic dust maps, including all infrared-based maps as well as "stellar" reddening maps. When such maps are used for extinction corrections, this extragalactic contamination introduces redshift- and scale-dependent biases in photometric estimates at the millimagnitude level. It can affect both object-based analyses, such as the estimation of the Hubble diagram with supernovae, as well as spatial statistics. The bias can be appreciable when measuring angular correlation functions with low amplitudes, such as lensing-induced correlations or angular correlations for sources distributed over a broad redshift range. As expected, we do not detect any extragalactic contamination for the dust map inferred from 21cm HI observations. Such a map provides an alternative to widely used infrared-based maps but relies on the assumption of a constant dust-to-gas ratio. We note that, using the WISE 12 micron map sensitive to polycyclic aromatic hydrocarbons (PAH), an indirect dust tracer, we detect the diffuse extragalactic PAH background up to $z\sim2$. Finally, we provide a procedure to minimize the level of biased magnitude corrections in maps with extragalactic imprints., Comment: Accepted for publication in the Astrophysical Journal

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress., Comment: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop., Comment: 50 pages, 24 figures, more details at http://spherex.caltech.edu

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project., Comment: 47 pages, 18 figures; accepted for publication in the Astronomical Journal

We investigate the limitations of statistical absorption measurements with the SDSS optical spectroscopic surveys. We show that changes in the data reduction strategy throughout different data releases have led to a better accuracy at long wavelengths, in particular for sky line subtraction, but a degradation at short wavelengths with the emergence of systematic spectral features with an amplitude of about one percent. We show that these features originate from inaccuracy in the fitting of modeled F-star spectra used for flux calibration. The best-fit models for those stars are found to systematically over-estimate the strength of metal lines and under-estimate that of Lithium. We also identify the existence of artifacts due to masking and interpolation procedures at the wavelengths of the hydrogen Balmer series leading to the existence of artificial Balmer $\alpha$ absorption in all SDSS optical spectra. All these effects occur in the rest-frame of the standard stars and therefore present Galactic longitude variations due to the rotation of the Galaxy. We demonstrate that the detection of certain weak absorption lines reported in the literature are solely due to calibration effects. Finally, we discuss new strategies to mitigate these issues., Comment: 11 pages, 7 figures, MNRAS accepted

We present a morphological study of the 17 lensed Lyman-$\alpha$ emitter (LAE) galaxies of the Baryon Oscillation Spectroscopic Survey Emission-Line Lens Survey (BELLS) for the GALaxy-Ly$\alpha$ EmitteR sYstems (BELLS GALLERY) sample. This analysis combines the magnification effect of strong galaxy-galaxy lensing with the high resolution of the Hubble Space Telescope ($HST$) to achieve a physical resolution of $\sim$80 pc for this $2

We use the latest version of the metal line absorption catalogue of Zhu & M\'enard (2013) to study the clustering of MgII absorbers around massive galaxies (~10^11.5 M_sun), quasars and radio-loud AGN with redshifts between 0.4 and 0.75. Clustering is evaluated in two dimensions, by binning absorbers both in projected radius and in velocity separation. Excess MgII is detected around massive galaxies out to R_p=20 Mpc. At projected radii less than 800 kpc, the excess extends out to velocity separations of 10,000 km/s. The extent of the high velocity tail within this radius is independent of the mean stellar age of the galaxy and whether or not it harbours an active galactic nucleus. We interpret our results using the publicly available Illustris and Millennium simulations. Models where the MgII absorbers trace the dark matter particle or subhalo distributions do not fit the data. They overpredict the clustering on small scales and do not reproduce the excess high velocity separation MgII absorbers seen within the virial radius of the halo. The Illustris simulations which include thermal, but not mechanical feedback from AGN, also do not provide an adequate fit to the properties of the cool halo gas within the virial radius. We propose that the large velocity separation MgII absorbers trace gas that has been pushed out of the dark matter halos, possibly by multiple episodes of AGN-driven mechanical feedback acting over long timescales., Comment: 10 pages, 11 figures, accepted in MNRAS

Using a sample of nearly half million galaxies, intersected by over 7 million lines of sight from the Sloan Digital Sky Survey Data Release 12, we trace H$\alpha$ + [N{\small II}] emission from a galactocentric projected radius, $r_p$, of 5 kpc to more than 100 kpc. The emission flux surface brightness is $\propto r_p^{-1.9 \pm 0.4}$. We obtain consistent results using only the H$\alpha$ or [N{\small II}] flux. We measure a stronger signal for the bluer half of the target sample than for the redder half on small scales, $r_p <$ 20 kpc. We obtain a $3\sigma$ detection of H$\alpha$ + [N{\small II}] emission in the 50 to 100 kpc $r_p$ bin. The mean emission flux within this bin is $(1.10 \pm 0.35) \times 10^{-20}$ erg cm$^{-2}$ s$^{-1}$ \AA$^{-1}$, which corresponds to $1.87 \times 10^{-20}$ erg cm$^{-2}$ s$^{-1}$ arcsec$^{-2}$ or 0.0033 Rayleigh. This detection is 34 times fainter than a previous strict limit obtained using deep narrow-band imaging. The faintness of the signal demonstrates why it has been so difficult to trace recombination radiation out to large radii around galaxies. This signal, combined with published estimates of n$_{\rm H}$, lead us to estimate the temperature of the gas to be 12,000 K, consistent with independent empirical estimates based on metal ion absorption lines and expectations from numerical simulations., Comment: 12 pages, 13 figures

We present \textsl{Hubble Space Telescope} (\textsl{HST}) F606W-band imaging observations of 21 galaxy-Ly$\alpha$ emitter lens candidates in the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$\alpha$ EmitteR sYstems (BELLS GALLERY) survey. 17 systems are confirmed to be definite lenses with unambiguous evidence of multiple imaging. The lenses are primarily massive early-type galaxies (ETGs) at redshifts of approximately $0.55$, while the lensed sources are Ly$\alpha$ emitters (LAEs) at redshifts from 2 to 3. Although the \textsl{HST} imaging data are well fit by smooth lens models consisting of singular isothermal ellipsoids in an external shear field, a thorough exploration of dark substructures in the lens galaxies is required. The Einstein radii of the BELLS GALLERY lenses are on average $60\%$ larger than those of the BELLS lenses because of the much higher source redshifts which will allow a detailed investigation of the radius evolution of the mass profile in ETGs. With the aid of the average $\sim 13 \times$ lensing magnification, the LAEs are resolved to comprise individual star-forming knots of a wide range of properties with characteristic sizes from less than 100 pc to several kpc, rest-frame far UV apparent AB magnitudes from 29.6 to 24.2, and typical projected separations of 500 pc to 2 kpc., Comment: 15 pages, 4 figures, minor edits to match the ApJ published version

We compare predictions of a number of empirical models and numerical simulations of galaxy formation to the conditional stellar mass functions (CSMF)of galaxies in groups of different masses obtained recently by Lan et al. to test how well different models accommodate the data. The observational data clearly prefer a model in which star formation in low-mass halos changes behavior at a characteristic redshift $z_c\sim 2$. There is also tentative evidence that this characteristic redshift depends on environment, becoming $z_c\sim 4$ in regions that eventually evolve into rich clusters of galaxies. The constrained model is used to understand how galaxies form and evolve in dark matter halos, and to make predictions for other statistical properties of the galaxy population, such as the stellar mass functions of galaxies at high $z$, the star formation and stellar mass assembly histories in dark matter halos. A comparison of our model predictions with those of other empirical models shows that different models can make vastly different predictions, even though all of them are tuned to match the observed stellar mass functions of galaxies., Comment: 17 pages, 10 figures, accepted for publication in MNRAS

We introduce the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$\alpha$ EmitteR sYstems (BELLS GALLERY) Survey, which is a Hubble Space Telescope program to image a sample of galaxy-scale strong gravitational lens candidate systems with high-redshift Ly$\alpha$ emitters (LAEs) as the background sources. The goal of the BELLS GALLERY Survey is to illuminate dark substructures in galaxy-scale halos by exploiting the small-scale clumpiness of rest-frame far-UV emission in lensed LAEs, and to thereby constrain the slope and normalization of the substructure-mass function. In this paper, we describe in detail the spectroscopic strong-lens selection technique, which is based on methods adopted in the previous Sloan Lens ACS (SLACS) Survey, BELLS, and SLACS for the Masses Survey. We present the BELLS GALLERY sample of the 21 highest-quality galaxy--LAE candidates selected from $\approx 1.4 \times 10^6$ galaxy spectra in the BOSS of the Sloan Digital Sky Survey III. These systems consist of massive galaxies at redshifts of approximately 0.5 strongly lensing LAEs at redshifts from 2--3. The compact nature of LAEs makes them an ideal probe of dark substructures, with a substructure-mass sensitivity that is unprecedented in other optical strong-lens samples. The magnification effect from lensing will also reveal the structure of LAEs below 100 pc scales, providing a detailed look at the sites of the most concentrated unobscured star formation in the universe. The source code used for candidate selection is available for download as a part of this release., Comment: 14 pages, 5 figures, accepted for publication in the ApJ (ApJ, 824, 86). Minor edits to match the ApJ published version

We apply clustering-based redshift inference to all extended sources from the Sloan Digital Sky Survey photometric catalogue, down to magnitude r = 22. We map the relationships between colours and redshift, without assumption of the sources' spectral energy distributions (SED). We identify and locate star-forming, quiescent galaxies, and AGN, as well as colour changes due to spectral features, such as the 4000 \AA{} break, redshifting through specific filters. Our mapping is globally in good agreement with colour-redshift tracks computed with SED templates, but reveals informative differences, such as the need for a lower fraction of M-type stars in certain templates. We compare our clustering-redshift estimates to photometric redshifts and find these two independent estimators to be in good agreement at each limiting magnitude considered. Finally, we present the global clustering-redshift distribution of all Sloan extended sources, showing objects up to z ~ 0.8. While the overall shape agrees with that inferred from photometric redshifts, the clustering redshift technique results in a smoother distribution, with no indication of structure in redshift space suggested by the photometric redshift estimates (likely artifacts imprinted by their spectroscopic training set). We also infer a higher fraction of high redshift objects. The mapping between the four observed colours and redshift can be used to estimate the redshift probability distribution function of individual galaxies. This work is an initial step towards producing a general mapping between redshift and all available observables in the photometric space, including brightness, size, concentration, and ellipticity., Comment: 12 pages, 9 figures, accepted to MNRAS

We characterize the luminosity functions of galaxies residing in $z\sim0$ groups and clusters over the broadest ranges of luminosity and mass reachable by the Sloan Digital Sky Survey. Our measurements cover four orders of magnitude in luminosity, down to about $M_r=-12$ mag or $L=10^7\,L_\odot$, and three orders of magnitude in halo mass, from $10^{12}$ to $10^{15} \, {\rm M}_\odot$. We find a characteristic scale, $M_r\sim-18$ mag or $L\sim10^9\, L_\odot$, below which the slope of the luminosity function becomes systematically steeper. This trend is present for all halo masses and originates mostly from red satellites. This ubiquitous faint-end upturn suggests that it is formation, rather than halo-specific environmental effect, that plays a major role in regulating the stellar masses of faint satellites. We show that the satellite luminosity functions can be described in a simple manner by a double Schechter function with amplitudes scaling with halo mass over the entire range of observables. Combining these conditional luminosity functions with the dark matter halo mass function, we accurately recover the entire field luminosity function over 10 visual magnitudes and reveal that satellite galaxies dominate the field luminosity function at magnitudes fainter than $-17$. We find that the luminosity functions of blue and red satellite galaxies show distinct shapes and we present estimates of the stellar mass fraction as a function of halo mass and galaxy type. Finally, using a simple model, we demonstrate that the abundances and the faint-end slopes of blue and red satellite galaxies can be interpreted in terms of their formation history, with two distinct modes separated by some characteristic time., Comment: 24 pages, 17 figures. Figure 7 shows the main result. Matched the published version. The data can be found at http://www.pha.jhu.edu/~tlan/research/CLFs/

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.92.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS., Comment: 39 pages, 13 figures, accepted AJ, revised for consistency with accepted version

The Two-Micron All-Sky Survey (2MASS) has mapped out the low-redshift Universe down to $K_S\sim14$ mag. As its near-infrared photometry primarily probes the featureless Rayleigh-Jeans tail of galaxy spectral energy distributions, colour-based redshift estimation is rather uninformative. Until now, redshift estimates for this dataset have relied on optical follow-up suffering from selection biases. Here we use the newly-developed technique of clustering-based redshift estimation to infer the redshift distribution of the 2MASS sources regardless of their optical properties. We characterise redshift distributions of objects from the Extended Source Catalogue as a function of near-infrared colours and brightness and report some observed trends. We also apply the clustering redshift technique to dropout populations, sources with non-detections in one or more near-infrared bands, and present their redshift distributions. Combining all extended sources, we confirm with clustering redshifts that the distribution of this sample extends up to $z\sim0.35$. % We perform a similar analysis with the Point Source Catalogue and show that it can be separated into stellar and extragalactic contributions with galaxies reaching $z\sim0.7$. We estimate that the Point Source Catalogue contains 1.6 million extragalactic objects: as many as in the Extended Source Catalogue but probing a cosmic volume ten times larger., Comment: 10 pages, 7 figures, supplementary data available online. Accepted to MNRAS

The Prime Focus Spectrograph (PFS) is an optical/near-infrared multifiber spectrograph with 2394 science fibers distributed across a 1.3-deg diameter field of view at the Subaru 8.2-m telescope. The wide wavelength coverage from 0.38 {\mu}m to 1.26 {\mu}m, with a resolving power of 3000, simultaneously strengthens its ability to target three main survey programs: cosmology, galactic archaeology and galaxy/AGN evolution. A medium resolution mode with a resolving power of 5000 for 0.71 {\mu}m to 0.89 {\mu}m will also be available by simply exchanging dispersers. We highlight some of the technological aspects of the design. To transform the telescope focal ratio, a broad-band coated microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of the cable system, optimizing overall throughput; a fiber with low focal ratio degradation is selected for the fiber-positioner and fiber-slit components, minimizing the effects of fiber movements and fiber bending. Fiber positioning will be performed by a positioner consisting of two stages of piezo-electric rotary motors. The positions of these motors are measured by taking an image of artificially back-illuminated fibers with the metrology camera located in the Cassegrain container; the fibers are placed in the proper location by iteratively measuring and then adjusting the positions of the motors. Target light reaches one of the four identical fast-Schmidt spectrograph modules, each with three arms. The PFS project has passed several project-wide design reviews and is now in the construction phase., Comment: 19 pages, 11 figures

We map out calcium II & sodium I absorption (Fraunhofer H, K & D lines) induced by both the interstellar medium and the circumgalactic medium of the Milky Way. Our measurements cover more than $9000$ deg$^2$ and make use of about $300,000$ extragalactic spectra from the Sloan Digital Sky Survey. We present absorption maps for these two species and then compare their distributions to those of neutral hydrogen and dust. We show that the abundance of Na I with respect to neutral hydrogen stays roughly constant in different environments, while that of Ca II decreases with hydrogen column density. Studying how these tracers vary as a function of velocity, we show that, on average, the N(Na I)/N(Ca II) ratio decreases at higher velocity with respect to the local standard of rest, similar to the local Routly-Spitzer effect but seen on Galactic scale. We show that it is likely caused by higher gas/dust density at lower velocity. Finally, we show that Galactic Ca II and Na I absorption needs to be taken into account for precision photometry and, more importantly, for photometric redshift estimation with star forming galaxies. Our maps of Ca II and Na I absorption are publicly available., Comment: Accepted to MNRAS

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra., Comment: DR12 data are available at http://www.sdss3.org/dr12. 30 pages. 11 figures. Accepted to ApJS

Using spectroscopically selected galaxies from the Sloan Digital Sky Survey we present a detection of reddening due to dust in the circumgalactic medium of galaxies. We detect the mean change in the colors of "standard crayons" correlated with the presence of foreground galaxies at z ~0.05 as a function of angular separation. Following Peek & Graves (2010), we create standard crayons using passively evolving galaxies corrected for Milky Way reddening and color-redshift trends, leading to a sample with as little as 2% scatter in color. We devise methods to ameliorate possible systematic effects related to the estimation of colors, and we find an excess reddening induced by foreground galaxies at a level ranging from 10 to 0.5 millimagnitudes on scales ranging from 30 kpc to 1 Mpc. We attribute this effect to a large-scale distribution of dust around galaxies similar to the findings of Menard et al. 2010. We find that circumgalactic reddening is a weak function of stellar mass over the range $6 \times 10^9 M_\odot$ -- $6 \times 10^{10} M_\odot$ and note that this behavior appears to be consistent with recent results on the distribution of metals in the gas phase., Comment: Submitted to ApJ

The Prime Focus Spectrograph (PFS) is an optical/near-infrared multi-fiber spectrograph with 2394 science fibers, which are distributed in 1.3 degree diameter field of view at Subaru 8.2-meter telescope. The simultaneous wide wavelength coverage from 0.38 um to 1.26 um, with the resolving power of 3000, strengthens its ability to target three main survey programs: cosmology, Galactic archaeology, and galaxy/AGN evolution. A medium resolution mode with resolving power of 5000 for 0.71 um to 0.89 um also will be available by simply exchanging dispersers. PFS takes the role for the spectroscopic part of the Subaru Measurement of Images and Redshifts project, while Hyper Suprime-Cam works on the imaging part. To transform the telescope plus WFC focal ratio, a 3-mm thick broad-band coated glass-molded microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of cable system, while one with a better FRD performance is selected for the fiber-positioner and fiber-slit components, given the more frequent fiber movements and tightly curved structure. Each Fiber positioner consists of two stages of piezo-electric rotary motors. Its engineering model has been produced and tested. Fiber positioning will be performed iteratively by taking an image of artificially back-illuminated fibers with the Metrology camera located in the Cassegrain container. The camera is carefully designed so that fiber position measurements are unaffected by small amounts of high special-frequency inaccuracies in WFC lens surface shapes. Target light carried through the fiber system reaches one of four identical fast-Schmidt spectrograph modules, each with three arms. Prototype VPH gratings have been optically tested. CCD production is complete, with standard fully-depleted CCDs for red arms and more-challenging thinner fully-depleted CCDs with blue-optimized coating for blue arms., Comment: 14 pages, 12 figures, submitted to "Ground-based and Airborne Instrumentation for Astronomy V, Suzanne K. Ramsay, Ian S. McLean, Hideki Takami, Editors, Proc. SPIE 9147 (2014)"

We investigate the potential and accuracy of clustering-based redshift estimation using the method proposed by M\'enard et al. (2013). This technique enables the inference of redshift distributions from measurements of the spatial clustering of arbitrary sources, using a set of reference objects for which redshifts are known. We apply it to a sample of spectroscopic galaxies from the Sloan Digital Sky Survey and show that, after carefully controlling the sampling efficiency over the sky, we can estimate redshift distributions with high accuracy. Probing the full colour space of the SDSS galaxies, we show that we can recover the corresponding mean redshifts with an accuracy ranging from $\delta$z=0.001 to 0.01. We indicate that this mapping can be used to infer the redshift probability distribution of a single galaxy. We show how the lack of information on the galaxy bias limits the accuracy of the inference and show comparisons between clustering redshifts and photometric redshifts for this dataset. This analysis demonstrates, using real data, that clustering-based redshift inference provides a powerful data-driven technique to explore the redshift distribution of arbitrary datasets, without any prior knowledge on the spectral energy distribution of the sources., Comment: 13 pages. Submitted to MNRAS. Comments welcome

Cross-correlating the Planck High Frequency Instrument (HFI) maps against quasars from the Sloan Digital Sky Survey (SDSS) DR7, we estimate the intensity distribution of the Cosmic Infrared Background (CIB) over the redshift range 0 < z < 5.We detect redshift-dependent spatial cross-correlations between the two datasets using the 857, 545 and 353 GHz channels and we obtain upper limits at 217 GHz consistent with expectations. At all frequencies with detectable signal we infer a redshift distribution peaking around z ~ 1.2 and find the recovered spectrum to be consistent with emission arising from star forming galaxies. By assuming simple modified blackbody and Kennicutt relations, we estimate dust and star formation rate density as a function of redshift, finding results consistent with earlier multiwavelength measurements over a large portion of cosmic history. However, we note that, lacking mid-infrared coverage, we are not able to make an accurate determination of the mean temperature for the dust responsible for the CIB. Our results demonstrate that clustering-based redshift inference is a valuable tool for measuring the entire evolution history of the cosmic star formation rate from a single and homogeneous dataset., Comment: 14 pages, 10 figures, published in MNRAS. Some revisions and an additional figure added since version 1

We use star, galaxy and quasar spectra taken by the Sloan Digital Sky Survey to map out the distribution of diffuse interstellar bands (DIBs) induced by the Milky Way. We show that, after carefully removing the intrinsic spectral energy distribution of each source, it is possible to measure statistical flux fluctuations at the 1e-3 level, detect more than 20 DIBs and measure their strength as a function of position on the sky. We create a map of DIB absorption covering about 5000 square degrees and measure correlations with various tracers of the interstellar medium: atomic and molecular hydrogen, dust and polycyclic aromatic hydrocarbons (PAHs). After recovering known correlations, we show that each DIB has a different dependence on atomic and molecular hydrogen: while they are all positively correlated with N(HI), they exhibit a range of behaviours with N(H2) showing positive, negative or no correlation. We show that a simple parametrization involving only N(HI) and N(H2) applied to all the DIBs is sufficient to reproduce a large collection of observational results reported in the literature: it allows us to naturally describe the relations between DIB strength and dust reddening (including the so-called skin effect), the related scatter, DIB pair-wise correlations & families, the affinity for $\sigma/\zeta$-type environments and other correlations related to molecules. Our approach allows us to characterize DIB dependencies in a simple manner and provides us with a metric to characterize the similarity between different DIBs., Comment: Match the published version. The DIB map can be viewed interactively at http://www.pha.jhu.edu/~tlan/DIB_SDSS/

We explore the distribution of cool (~$10^4$K) gas around galaxies and its dependence on galaxy properties. By cross-correlating about 50,000 MgII absorbers with millions of sources from the SDSS (optical), WISE (IR), and GALEX (UV) surveys we effectively extract about 2,000 galaxy-absorber pairs at z~0.5 and probe relations between absorption strength and galaxy type, impact parameter and azimuthal angle. We find that cool gas traced by MgII absorbers exists around both star-forming and passive galaxies with a similar incidence rate on scales greater than 100 kpc but each galaxy type exhibits a different behavior on smaller scales: MgII equivalent width does not correlate with the presence of passive galaxies whereas stronger MgII absorbers tend to be found in the vicinity of star-forming galaxies. This effect is preferentially seen along the minor axis of these galaxies, suggesting that some of the gas is associated with outflowing material. In contrast, the distribution of cool gas around passive galaxies is consistent with being isotropic on the same scales. We quantify the average excess MgII equivalent width $<\delta W_{0}^{\rm MgII}>$ as a function of galaxy properties and find $<\delta W_0^{\rm MgII}>\propto SFR^{1.2}, sSFR^{0.5}$ and $M_\ast^{0.4}$ for star-forming galaxies. This work demonstrates that the dichotomy between star-forming and passive galaxies is reflected in the CGM traced by low-ionized gas. We also measure the covering fraction of MgII absorption and find it to be about 2-10 times higher for star-forming galaxies than passive ones within 50 kpc. We estimate the amount of neutral gas in the halo of $<\log M_\ast/{\rm M_\odot}>$~10.8 galaxies to be a few x$10^9 \rm M_\odot$ for both types of galaxies. Finally, we find that correlations between absorbers and sources detected in the UV and IR lead to physical trends consistent with those measured in the optical., Comment: 14 pages, 11 figures, 1 table, updated to match the published version

We estimate the abundance of dust in damped Lyman-{\alpha} absorbers (DLAs) by statistically measuring the excess reddening they induce on their background quasars. We detect systematic reddening behind DLA consistent with the SMC type reddening curve, but it is inconsistent with the Milky Way type reddening. We find that the derived dust-to-gas ratio is, on average, inversely proportional to the column density of neutral hydrogen, implying that the amount of dust is constant, irrespective of the column density of hydrogen. It means that the average metallicity is inversely proportional to the column density of hydrogen, unless the average dust-to-metal ratio varies with the hydrogen column density. This indicates that the prime origin of metals seen in DLAs is not by in situ star formation, with which Z ~ N_{HI}^{0.4} is expected from the empirical star formation law, contrary to our observation. We interpret the metals observed in absorbers being deposited dominantly from nearby galaxies by galactic winds ubiquitous in intergalactic space. When extrapolating the relation between dust-to-gas ratio and HI column density to lower column density, we find a value which is consistent with what is observed for Mg II absorbers., Comment: 8 pages, 5 figures. Version accepted for publication in the Astrophysical Journal

We have obtained spectra of 163 quasars at $z_\mathrm{em}>4.4$ with the Gemini Multi Object Spectrometers on the Gemini North and South telescopes, the largest publicly available sample of high-quality, low-resolution spectra at these redshifts. From this homogeneous data set, we generated stacked quasar spectra in three redshift intervals at $z\sim 5$. We have modelled the flux below the rest-frame Lyman limit ($\lambda_\mathrm{r}<912$\AA) to assess the mean free path $\lambda_\mathrm{mfp}^{912}$ of the intergalactic medium to HI-ionizing radiation. At mean redshifts $z_\mathrm{q}=4.56$, 4.86 and 5.16, we measure $\lambda_\mathrm{mfp}^{912}=(22.2\pm 2.3, 15.1\pm 1.8, 10.3\pm 1.6)h_{70}^{-1}$ proper Mpc with uncertainties dominated by sample variance. Combining our results with $\lambda_\mathrm{mfp}^{912}$ measurements from lower redshifts, the data are well modelled by a simple power-law $\lambda_\mathrm{mfp}^{912}=A[(1+z)/5]^\eta$ with $A=(37\pm 2)h_{70}^{-1}$ Mpc and $\eta = -5.4\pm 0.4$ between $z=2.3$ and $z=5.5$. This rapid evolution requires a physical mechanism -- beyond cosmological expansion -- which reduces the cosmic effective Lyman limit opacity. We speculate that the majority of HI Lyman limit opacity manifests in gas outside galactic dark matter haloes, tracing large-scale structures (e.g. filaments) whose average density (and consequently neutral fraction) decreases with cosmic time. Our measurements of the strongly redshift-dependent mean free path shortly after the completion of HI reionization serve as a valuable boundary condition for numerical models thereof. Having measured $\lambda_\mathrm{mfp}^{912}\approx 10$ Mpc at $z=5.2$, we confirm that the intergalactic medium is highly ionized by that epoch and that the redshift evolution of the mean free path does not show a break that would indicate a recent end to HI reionization., Comment: 16 pages, 11 figures, plus 20 pages of online material; submitted to MNRAS

We present a measurement of the correlation function between luminous red galaxies and cool gas traced by Mg II \lambda \lambda 2796, 2803 absorption, on scales ranging from about 30 kpc to 20 Mpc. The measurement is based on cross-correlating the positions of about one million red galaxies at z~0.5 and the flux decrements induced in the spectra of about 10^5 background quasars from the Sloan Digital Sky Survey. We find that: (i) This galaxy-gas correlation reveals a change of slope on scales of about 1 Mpc, consistent with the expected transition from a dark matter halo dominated environment to a regime where clustering is dominated by halo-halo correlations. Assuming that, on average, the distribution of Mg II gas follows that of dark matter up to a gas-to-mass ratio, we find the standard halo model to provide an accurate description of the gas distribution over three orders of magnitude in scale. Within this framework we estimate the average host halo mass of luminous red galaxies to be about 10^{13.5} M_solar, in agreement with other methods. We also find the Mg II gas-to-mass ratio around LRGs to be consistent with the cosmic value estimated on Mpc scales. Combining our galaxy-gas correlation and the galaxy-mass correlation function from galaxy-galaxy lensing analyses we can directly measure the Mg II gas-to-mass ratio as a function of scale and reach the same conclusion. (ii) From line-width estimates, we show that the velocity dispersion of the gas clouds also shows the expected 1- and 2-halo behaviors. On large scales the gas distribution follows the Hubble flow, whereas on small scales we observe the velocity dispersion of the Mg II gas clouds to be lower than that of collisionless dark matter particles within their host halo. This is in line with the fact that cool clouds are subject to the pressure of the virialized hot gas., Comment: 18 pages, 11 figures, 1 table, submitted to MNRAS

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees., Comment: 15 figures; 1 table. Accepted to ApJS. DR10 is available at http://www.sdss3.org/dr10 v3 fixed 3 diacritic markings in the arXiv HTML listing of the author names

We present a measurement of the mean density profile of Ca II gas around galaxies out to ~ 200 kpc, traced by Fraunhofer's H & K absorption lines. The measurement is based on cross-correlating the positions of about one million foreground galaxies at z ~ 0.1 and the flux decrements induced in the spectra of about 10^5 background quasars from the Sloan Digital Sky Survey. This technique allows us to trace the total amount of Ca II absorption induced by the circumgalactic medium, including absorbers too weak to be detected in individual spectra. We can statistically measure Ca II rest equivalent widths down to several mA, corresponding to column densities of about 5x10^10 cm^{-2}. We find that the Ca II column density distribution follows N ~ rp^{-1.4} and the mean Ca II mass in the halo within 200 kpc is ~ 5x10^3 Msolar, averaged over the foreground galaxy sample with median mass ~ 10^10.3 Msolar. This is about an order-of-magnitude larger than the Ca II mass in the interstellar medium of the Milky Way, suggesting more than 90% of Ca II in the Universe is in the circum- and inter-galactic environments. Our measurements indicate that the amount of Ca II in halos is larger for galaxies with higher stellar mass and higher star formation rate. For edge-on galaxies we find Ca II to be more concentrated along the minor axis, i.e. in the polar direction. This suggests that bipolar outflows induced by star formation must have played a significant role in producing Ca II in galaxy halos., Comment: 10 pages, 9 figures, 2 tables. submitted to ApJ. Introduction includes a review of Fraunhofer's line labeling. Comments are welcome

We present a data-driven method to infer the redshift distribution of an arbitrary dataset based on spatial cross-correlation with a reference population and we apply it to various datasets across the electromagnetic spectrum to show its potential and limitations. Our approach advocates the use of clustering measurements on all available scales, in contrast to previous works focusing only on linear scales. We also show how its accuracy can be enhanced by optimally sampling a dataset within its photometric space rather than applying the estimator globally. We show that the ultimate goal of this technique is to characterize the mapping between the space of photometric observables and redshift space as this characterization then allows us to infer the clustering-redshift p.d.f. of a single galaxy. We apply this technique to estimate the redshift distributions of luminous red galaxies and emission line galaxies from the SDSS, infrared sources from WISE and radio sources from FIRST. We show that consistent redshift distributions are found using both quasars and absorber systems as reference populations. This technique brings valuable information on the third dimension of astronomical datasets. It is widely applicable to a large range of extra-galactic surveys., Comment: 10 pages, 5 figures. Improved description of the formalism

Determining accurate redshift distributions for very large samples of objects has become increasingly important in cosmology. We investigate the impact of extending cross-correlation based redshift distribution recovery methods to include small scale clustering information. The major concern in such work is the ability to disentangle the amplitude of the underlying redshift distribution from the influence of evolving galaxy bias. Using multiple simulations covering a variety of galaxy bias evolution scenarios, we demonstrate reliable redshift recoveries using linear clustering assumptions well into the non-linear regime for redshift distributions of narrow redshift width. Including information from intermediate physical scales balances the increased information available from clustering and the residual bias incurred from relaxing of linear constraints. We discuss how breaking a broad sample into tomographic bins can improve estimates of the redshift distribution, and present a simple bias removal technique using clustering information from the spectroscopic sample alone., Comment: 13 pages, 9 figures. Submitted to MNRAS