Showing 1–50 of 115 results for author: Taylor, J E

The infall region as a complementary probe to cluster abundance

Authors: Charlie T. Mpetha, James E. Taylor, Yuba Amoura, Roan Haggar

Abstract: Galaxy cluster abundance measurements provide a classic test of cosmology. They are most sensitive to the evolved amplitude of fluctuations, usually expressed as $S_8 = σ_8\sqrt{Ω_m/0.3}$. Thus, abundance constraints exhibit a strong degeneracy between $σ_8$ and $Ω_{\rm m}$, as do other similar low-redshift tests such as cosmic shear. The mass distribution in the infall region around galaxy cluste… ▽ More Galaxy cluster abundance measurements provide a classic test of cosmology. They are most sensitive to the evolved amplitude of fluctuations, usually expressed as $S_8 = σ_8\sqrt{Ω_m/0.3}$. Thus, abundance constraints exhibit a strong degeneracy between $σ_8$ and $Ω_{\rm m}$, as do other similar low-redshift tests such as cosmic shear. The mass distribution in the infall region around galaxy clusters, where material is being accreted from the surrounding field, also exhibits a cosmological dependence, but in this case it is nearly orthogonal to the $S_8$ direction in the $Ω_m$--$σ_8$ plane, making it highly complementary to halo abundance or cosmic shear studies. We explore how weak lensing measurements of the infall region might be used to complement abundance studies, considering three different tests. The splashback radius is a prominent feature of the infall region; we show that detection of this feature in lensing data from the Euclid survey could independently constrain $Ω_{\rm m}$ and $σ_8$ to $\pm 0.05$. Another feature, the depletion radius where the bias reaches a minimum, also shows cosmological dependence, though it is challenging to observe in practice. The strongest constraints come from direct measurements of the shear profile in the infall region at $2$--$4\,r_{200{\rm c}}$. Combining the latter with abundance constraints such as those reported from SRG$/$eROSITA should reduce the area of the error contours by an estimated factor of $1.2$ using a sample of clusters observed by the UNIONS survey, or a factor of $3$ using clusters observed by the Euclid Wide survey over a broader range of redshift. △ Less

Submitted 24 July, 2024; v1 submitted 1 July, 2024; originally announced July 2024.

Comments: 14 pages, 10 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 532, Issue 2, August 2024, Pages 2521-2533

Constraining cosmological parameters using the splashback radius of galaxy clusters

Authors: Roan Haggar, Yuba Amoura, Charlie T. Mpetha, James E. Taylor, Kris Walker, Chris Power

Abstract: Cosmological parameters such as $Ω_{\rm{M}}$ and $σ_{8}$ can be measured indirectly using various methods, including galaxy cluster abundance and cosmic shear. These measurements constrain the composite parameter $S_{8}$, leading to degeneracy between $Ω_{\rm{M}}$ and $σ_{8}$. However, some structural properties of galaxy clusters also correlate with cosmological parameters, due to their dependenc… ▽ More Cosmological parameters such as $Ω_{\rm{M}}$ and $σ_{8}$ can be measured indirectly using various methods, including galaxy cluster abundance and cosmic shear. These measurements constrain the composite parameter $S_{8}$, leading to degeneracy between $Ω_{\rm{M}}$ and $σ_{8}$. However, some structural properties of galaxy clusters also correlate with cosmological parameters, due to their dependence on a cluster's accretion history. In this work, we focus on the splashback radius, an observable cluster feature that represents a boundary between a cluster and the surrounding Universe. Using a suite of cosmological simulations with a range of values for $Ω_{\rm{M}}$ and $σ_{8}$, we show that the position of the splashback radius around cluster-mass halos is greater in cosmologies with smaller values of $Ω_{\rm{M}}$ or larger values of $σ_{8}$. This variation breaks the degeneracy between $Ω_{\rm{M}}$ and $σ_{8}$ that comes from measurements of the $S_{8}$ parameter. We also show that this variation is, in principle, measurable in observations. As the splashback radius can be determined from the same weak lensing analysis already used to estimate $S_{8}$, this new approach can tighten low-redshift constraints on cosmological parameters, either using existing data, or using upcoming data such as that from Euclid and LSST. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 12 pages, 7 figures, 1 table, accepted for publication in ApJ

Reconsidering the dynamical states of galaxy clusters using PCA and UMAP

Authors: Roan Haggar, Federico De Luca, Marco De Petris, Elizaveta Sazonova, James E. Taylor, Alexander Knebe, Meghan E. Gray, Frazer R. Pearce, Ana Contreras-Santos, Weiguang Cui, Ulrike Kuchner, Robert A. Mostoghiu Paun, Chris Power

Abstract: Numerous metrics exist to quantify the dynamical state of galaxy clusters, both observationally and within simulations. Many of these correlate strongly with one another, but it is not clear whether all of these measures probe the same intrinsic properties. In this work, we use two different statistical approaches -- principal component analysis (PCA) and uniform manifold approximation and project… ▽ More Numerous metrics exist to quantify the dynamical state of galaxy clusters, both observationally and within simulations. Many of these correlate strongly with one another, but it is not clear whether all of these measures probe the same intrinsic properties. In this work, we use two different statistical approaches -- principal component analysis (PCA) and uniform manifold approximation and projection (UMAP) -- to investigate which dynamical properties of a cluster are in fact the best descriptors of its dynamical state. We use measurements taken directly from The Three Hundred suite of galaxy cluster simulations, as well as morphological properties calculated using mock X-ray and SZ maps of the same simulated clusters. We find that four descriptions of dynamical state naturally arise, and although correlations exist between these, a given cluster can be "dynamically relaxed" according to all, none, or some of these four descriptions. These results demonstrate that it is highly important for future observational and theoretical studies to consider in which sense clusters are dynamically relaxed. Cluster dynamical states are complex and multi-dimensional, and so it is not meaningful to classify them simply as "relaxed" and "unrelaxed" based on a single linear scale. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 18 pages, 7 figures, 3 tables, accepted for publication in MNRAS

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, S. Alvi, A. Amara , et al. (1115 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 24 September, 2024; v1 submitted 22 May, 2024; originally announced May 2024.

Comments: Accepted for publication in the A&A special issue`Euclid on Sky'

Halo Growth and Merger Rates as a Cosmological Test

Authors: Yuba Amoura, Nicole E. Drakos, Anael Berrouet, James E. Taylor

Abstract: Dark matter haloes grow at a rate that depends on the value of the cosmological parameters $σ_8$ and $Ω_{\rm m}$ through the initial power spectrum and the linear growth factor. While halo abundance is routinely used to constrain these parameters, through cluster abundance studies, the halo growth rate is not. In recent work, we proposed constraining the cosmological parameters using observational… ▽ More Dark matter haloes grow at a rate that depends on the value of the cosmological parameters $σ_8$ and $Ω_{\rm m}$ through the initial power spectrum and the linear growth factor. While halo abundance is routinely used to constrain these parameters, through cluster abundance studies, the halo growth rate is not. In recent work, we proposed constraining the cosmological parameters using observational estimates of the overall dynamical "age" of clusters, expressed, for instance, by their half-mass assembly redshift $z_{50}$. Here we explore the prospects for using the instantaneous growth rate, as estimated from the halo merger rate, from the average growth rate over the last dynamical time, or from the fraction of systems with recent episodes of major growth. We show that the merger rate is mainly sensitive to the amplitude of fluctuations $σ_8$, while the rates of recent growth provide constraints in the $Ω_{\rm m}$-$σ_8$ plane that are almost orthogonal to those provided by abundance studies. Data collected for forthcoming cluster abundance studies, or studies of the galaxy merger rate in current and future galaxy surveys, may thus provide additional constraints on the cosmological parameters complementary to those already derived from halo abundance. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 16 pages, 12 figures, 1 table and two figures in appendix

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 2, January 2024, Pages 3459-3473

Testing the Surface Brightness Fluctuation Method on Dwarf Galaxies in the COSMOS Field

Authors: Lauren M. Foster, James E. Taylor, John P. Blakeslee

Abstract: Dwarf galaxies are important tracers of small-scale cosmological structure, yet much of our knowledge about these systems comes from the limited sample of dwarf galaxies within the Local Group. To make a comprehensive inventory of dwarf populations in the local Universe, we require effective methods for deriving distance estimates for large numbers of faint, low surface brightness objects. Here we… ▽ More Dwarf galaxies are important tracers of small-scale cosmological structure, yet much of our knowledge about these systems comes from the limited sample of dwarf galaxies within the Local Group. To make a comprehensive inventory of dwarf populations in the local Universe, we require effective methods for deriving distance estimates for large numbers of faint, low surface brightness objects. Here we test the surface brightness fluctuation (SBF) method, traditionally applied to brighter early-type galaxies, on a sample of 20 nearby dwarf galaxies detected in the COSMOS field. These objects are partially resolved in HST ACS images, and have confirmed redshift distances in the range 17-130 Mpc. We discuss the many model choices required in applying the SBF method, and explore how these affect the final distance estimates. Amongst other variations on the method, when applying the SBF method, we alter the standard equation to include a term accounting for the power spectrum of the background, greatly improving our results. For the most robust modelling choices, we find a roughly Gaussian SBF signal that correlates linearly with distance out to distances of 50-100 Mpc, but with only a fraction of the power expected. At larger distances, there is excess power relative to that predicted, probably from undetected point sources. Overall, obtaining accurate SBF distances to faint, irregular galaxies remains challenging, but may yet prove possible with the inclusion of more information about galaxy properties and point source populations, and the use of more advanced techniques. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 18 pages, 26 figures, accepted by MNRAS

Unbiased Estimation of Structured Prediction Error

Authors: Kevin Fry, Jonathan E. Taylor

Abstract: Many modern datasets, such as those in ecology and geology, are composed of samples with spatial structure and dependence. With such data violating the usual independent and identically distributed (IID) assumption in machine learning and classical statistics, it is unclear a priori how one should measure the performance and generalization of models. Several authors have empirically investigated c… ▽ More Many modern datasets, such as those in ecology and geology, are composed of samples with spatial structure and dependence. With such data violating the usual independent and identically distributed (IID) assumption in machine learning and classical statistics, it is unclear a priori how one should measure the performance and generalization of models. Several authors have empirically investigated cross-validation (CV) methods in this setting, reaching mixed conclusions. We provide a class of unbiased estimation methods for general quadratic errors, correlated Gaussian response, and arbitrary prediction function $g$, for a noise-elevated version of the error. Our approach generalizes the coupled bootstrap (CB) from the normal means problem to general normal data, allowing correlation both within and between the training and test sets. CB relies on creating bootstrap samples that are intelligently decoupled, in the sense of being statistically independent. Specifically, the key to CB lies in generating two independent "views" of our data and using them as stand-ins for the usual independent training and test samples. Beginning with Mallows' $C_p$, we generalize the estimator to develop our generalized $C_p$ estimators (GC). We show at under only a moment condition on $g$, this noise-elevated error estimate converges smoothly to the noiseless error estimate. We show that when Stein's unbiased risk estimator (SURE) applies, GC converges to SURE as in the normal means problem. Further, we use these same tools to analyze CV and provide some theoretical analysis to help understand when CV will provide good estimates of error. Simulations align with our theoretical results, demonstrating the effectiveness of GC and illustrating the behavior of CV methods. Lastly, we apply our estimator to a model selection task on geothermal data in Nevada. △ Less

Submitted 16 October, 2023; originally announced October 2023.

Comments: 28 pages, 13 figures

Do assumptions about the central density of subhaloes affect dark matter annihilation and lensing calculations?

Authors: Nicole E. Drakos, James E. Taylor, Andrew J. Benson

Abstract: A growing body of evidence suggests that the central density of cuspy dark matter subhaloes is conserved in minor mergers. However, empirical models of subhalo evolution, calibrated from simulations, often assume a drop in the central density. Since empirical models of subhaloes are used in galaxy-galaxy lensing studies and dark matter annihilation calculations, we explore the consequences of assu… ▽ More A growing body of evidence suggests that the central density of cuspy dark matter subhaloes is conserved in minor mergers. However, empirical models of subhalo evolution, calibrated from simulations, often assume a drop in the central density. Since empirical models of subhaloes are used in galaxy-galaxy lensing studies and dark matter annihilation calculations, we explore the consequences of assuming different subhalo models. We find that dark matter annihilation calculations are very sensitive to the assumed subhalo mass profile, and different models can give more than a magnitude difference in the J-factor and boost factor in individual haloes. On the other hand, the shear and convergence profiles used in galaxy-galaxy lensing are sensitive to the initial profile assumed (e.g., NFW versus Einato) but are otherwise well-approximated by a simple model in which the original profile is sharply truncated. We conclude that since the innermost parts of haloes are difficult to resolve in simulations, it is important to have a theoretical understanding of how subhaloes evolve to make accurate predictions of the dark matter annihilation signal. △ Less

Submitted 22 September, 2024; v1 submitted 1 August, 2023; originally announced August 2023.

Comments: 13 figures. To be submitted to JCAP

Exposing flaws of generative model evaluation metrics and their unfair treatment of diffusion models

Authors: George Stein, Jesse C. Cresswell, Rasa Hosseinzadeh, Yi Sui, Brendan Leigh Ross, Valentin Villecroze, Zhaoyan Liu, Anthony L. Caterini, J. Eric T. Taylor, Gabriel Loaiza-Ganem

Abstract: We systematically study a wide variety of generative models spanning semantically-diverse image datasets to understand and improve the feature extractors and metrics used to evaluate them. Using best practices in psychophysics, we measure human perception of image realism for generated samples by conducting the largest experiment evaluating generative models to date, and find that no existing metr… ▽ More We systematically study a wide variety of generative models spanning semantically-diverse image datasets to understand and improve the feature extractors and metrics used to evaluate them. Using best practices in psychophysics, we measure human perception of image realism for generated samples by conducting the largest experiment evaluating generative models to date, and find that no existing metric strongly correlates with human evaluations. Comparing to 17 modern metrics for evaluating the overall performance, fidelity, diversity, rarity, and memorization of generative models, we find that the state-of-the-art perceptual realism of diffusion models as judged by humans is not reflected in commonly reported metrics such as FID. This discrepancy is not explained by diversity in generated samples, though one cause is over-reliance on Inception-V3. We address these flaws through a study of alternative self-supervised feature extractors, find that the semantic information encoded by individual networks strongly depends on their training procedure, and show that DINOv2-ViT-L/14 allows for much richer evaluation of generative models. Next, we investigate data memorization, and find that generative models do memorize training examples on simple, smaller datasets like CIFAR10, but not necessarily on more complex datasets like ImageNet. However, our experiments show that current metrics do not properly detect memorization: none in the literature is able to separate memorization from other phenomena such as underfitting or mode shrinkage. To facilitate further development of generative models and their evaluation we release all generated image datasets, human evaluation data, and a modular library to compute 17 common metrics for 9 different encoders at https://github.com/layer6ai-labs/dgm-eval. △ Less

Submitted 30 October, 2023; v1 submitted 7 June, 2023; originally announced June 2023.

Comments: NeurIPS 2023. 53 pages, 29 figures, 12 tables. Code at https://github.com/layer6ai-labs/dgm-eval, reviews at https://openreview.net/forum?id=08zf7kTOoh

Journal ref: Thirty-seventh Conference on Neural Information Processing Systems (2023)

A universal model for the evolution of tidally stripped systems

Authors: Nicole E. Drakos, James E. Taylor, Andrew J. Benson

Abstract: Accurate models of the structural evolution of dark matter subhaloes, as they orbit within larger systems, are fundamental to understanding the detailed distribution of dark matter at the present day. Numerical simulations of subhalo evolution support the idea that the mass loss associated with tidal stripping is most naturally understood in energy space, with the particles that are the least boun… ▽ More Accurate models of the structural evolution of dark matter subhaloes, as they orbit within larger systems, are fundamental to understanding the detailed distribution of dark matter at the present day. Numerical simulations of subhalo evolution support the idea that the mass loss associated with tidal stripping is most naturally understood in energy space, with the particles that are the least bound being removed first. Starting from this premise, we recently proposed a zero-parameter "energy-truncation model" for subhalo evolution. We tested this model with simulations of tidal stripping of satellites with initial NFW profiles, and showed that the energy-truncation model accurately predicts both the mass loss and density profiles. In this work, we apply the model to a variety of Hernquist, Einasto and King profiles. We show that it matches the simulation results quite closely in all cases, indicating that it may serve as a universal model to describe tidally stripped collisionless systems. A key prediction of the energy-truncation model is that the central density of dark matter subhaloes is conserved as they lose mass; this has important implications for dark matter annihilation calculations, and for other observational tests of dark matter. △ Less

Submitted 29 July, 2022; originally announced July 2022.

Comments: 19 pages, 14 figures. Accepted for publication in MNRAS

A hierarchical clustering method for quantifying satellite abundance

Authors: Chengyu Xi, James E. Taylor

Abstract: We present a new method for quantifying the abundance of satellites around field galaxies and in groups. The method is designed to work with samples, such as local photometric redshift catalogues, that do not have full spectroscopic coverage, but for which some redshift or distance information is available. It consists of identifying the galaxies most likely to be centrals, and using the clusterin… ▽ More We present a new method for quantifying the abundance of satellites around field galaxies and in groups. The method is designed to work with samples, such as local photometric redshift catalogues, that do not have full spectroscopic coverage, but for which some redshift or distance information is available. It consists of identifying the galaxies most likely to be centrals, and using the clustering signal around them as a template to iteratively decompose the full population into satellite and central populations. In that sense it is similar to performing crowded-field photometry, after having first used isolated stars to determine the point spread function of the image. The method does not identify individual satellites or centrals conclusively, but assigns a probability to each galaxy of being one or the other. Averaged over a large sample, it provides a statistical estimate of satellite abundance, even in crowded fields with large redshift uncertainties. We test the method using data from the COSMOS field, which includes a large set of local objects with accurate photometric redshifts. We measure satellite abundance as a function of central stellar or halo mass, as well as the satellite luminosity function, and find results consistent with previous studies, but extending over a broader range of central masses. We also consider a number of possible systematic uncertainties in the method, and show that they are generally smaller than our random errors. Having presented the method in this paper, we will use it to study the properties of the satellite populations in a forthcoming one. △ Less

Submitted 20 December, 2021; originally announced December 2021.

Comments: 17 pages, 17 figures

Journal ref: 2021MNRAS.503.4976X

Cluster Assembly Times as a Cosmological Test

Authors: Yuba Amoura, Nicole E. Drakos, Anael Berrouet, James E. Taylor

Abstract: The abundance of galaxy clusters in the low-redshift universe provides an important cosmological test, constraining a product of the initial amplitude of fluctuations and the amount by which they have grown since early times. The degeneracy of the test with respect to these two factors remains a limitation of abundance studies. Clusters will have different mean assembly times, however, depending o… ▽ More The abundance of galaxy clusters in the low-redshift universe provides an important cosmological test, constraining a product of the initial amplitude of fluctuations and the amount by which they have grown since early times. The degeneracy of the test with respect to these two factors remains a limitation of abundance studies. Clusters will have different mean assembly times, however, depending on the relative importance of initial fluctuation amplitude and subsequent growth. Thus, structural probes of cluster age such as concentration, shape or substructure may provide a new cosmological test that breaks the main degeneracy in number counts. We review analytic predictions for how mean assembly time should depend on cosmological parameters, and test these predictions using cosmological simulations. Given the overall sensitivity expected, we estimate the cosmological parameter constraints that could be derived from the cluster catalogues of forthcoming surveys such as Euclid, the Nancy Grace Roman Space Telescope, eROSITA, or CMB-S4. We show that by considering the structural properties of their cluster samples, such surveys could easily achieve errors of $Δσ_8$ = 0.01 or better. △ Less

Submitted 18 September, 2021; originally announced September 2021.

Comments: 18 pages, 17 figures

The volume-of-tube method for Gaussian random fields with inhomogeneous variance

Authors: Satoshi Kuriki, Akimichi Takemura, Jonathan E. Taylor

Abstract: The tube method or the volume-of-tube method approximates the tail probability of the maximum of a smooth Gaussian random field with zero mean and unit variance. This method evaluates the volume of a spherical tube about the index set, and then transforms it to the tail probability. In this study, we generalize the tube method to a case in which the variance is not constant. We provide the volume… ▽ More The tube method or the volume-of-tube method approximates the tail probability of the maximum of a smooth Gaussian random field with zero mean and unit variance. This method evaluates the volume of a spherical tube about the index set, and then transforms it to the tail probability. In this study, we generalize the tube method to a case in which the variance is not constant. We provide the volume formula for a spherical tube with a non-constant radius in terms of curvature tensors, and the tail probability formula of the maximum of a Gaussian random field with inhomogeneous variance, as well as its Laplace approximation. In particular, the critical radius of the tube is generalized for evaluation of the asymptotic approximation error. As an example, we discuss the approximation of the largest eigenvalue distribution of the Wishart matrix with a non-identity matrix parameter. The Bonferroni method is the tube method when the index set is a finite set. We provide the formula for the asymptotic approximation error for the Bonferroni method when the variance is not constant. △ Less

Submitted 9 September, 2021; v1 submitted 4 August, 2021; originally announced August 2021.

Comments: 30 pages, 3 figures

MSC Class: 62H10 (Primary); 60G60 (Secondary)

How do we generalize?

Authors: Jessica Elizabeth Taylor, Aurelio Cortese, Helen C. Barron, Xiaochuan Pan, Masamichi Sakagami, Dagmar Zeithamova

Abstract: Humans and animals are able to generalize or transfer information from previous experience so that they can behave appropriately in novel situations. What mechanisms--computations, representations, and neural systems--give rise to this remarkable ability? The members of this Generative Adversarial Collaboration (GAC) come from a range of academic backgrounds but are all interested in uncovering th… ▽ More Humans and animals are able to generalize or transfer information from previous experience so that they can behave appropriately in novel situations. What mechanisms--computations, representations, and neural systems--give rise to this remarkable ability? The members of this Generative Adversarial Collaboration (GAC) come from a range of academic backgrounds but are all interested in uncovering the mechanisms of generalization. We started out this GAC with the aim of arbitrating between two alternative conceptual accounts: (1) generalization stems from integration of multiple experiences into summary representations that reflect generalized knowledge, and (2) generalization is computed on-the-fly using separately stored individual memories. Across the course of this collaboration, we found that--despite using different terminology and techniques, and although some of our specific papers may provide evidence one way or the other--we in fact largely agree that both of these broad accounts (as well as several others) are likely valid. We believe that future research and theoretical synthesis across multiple lines of research is necessary to help determine the degree to which different candidate generalization mechanisms may operate simultaneously, operate on different scales, or be employed under distinct conditions. Here, as the first step, we introduce some of these candidate mechanisms and we discuss the issues currently hindering better synthesis of generalization research. Finally, we introduce some of our own research questions that have arisen over the course of this GAC, that we believe would benefit from future collaborative efforts. △ Less

Submitted 27 August, 2021; v1 submitted 2 April, 2021; originally announced April 2021.

Evolution of subhalo orbits in a smoothly-growing host halo potential

Authors: Go Ogiya, James E. Taylor, Michael J. Hudson

Abstract: The orbital parameters of dark matter (DM) subhaloes play an essential role in determining their mass-loss rates and overall spatial distribution within a host halo. Haloes in cosmological simulations grow by a combination of relatively smooth accretion and more violent mergers, and both processes will modify subhalo orbits. To isolate the impact of the smooth growth of the host halo from other re… ▽ More The orbital parameters of dark matter (DM) subhaloes play an essential role in determining their mass-loss rates and overall spatial distribution within a host halo. Haloes in cosmological simulations grow by a combination of relatively smooth accretion and more violent mergers, and both processes will modify subhalo orbits. To isolate the impact of the smooth growth of the host halo from other relevant mechanisms, we study subhalo orbital evolution using numerical calculations in which subhaloes are modelled as massless particles orbiting in a time-varying spherical potential. We find that the radial action of the subhalo orbit decreases over the first few orbits, indicating that the response to the growth of the host halo is not adiabatic during this phase. The subhalo orbits can shrink by a factor of $\sim$1.5 in this phase. Subsequently, the radial action is well conserved and orbital contraction slows down. We propose a model accurately describing the orbital evolution. Given these results, we consider the spatial distribution of the population of subhaloes identified in high-resolution cosmological simulations. We find that it is consistent with this population having been accreted at z < 3, indicating that any subhaloes accreted earlier are unresolved in the simulations. We also discuss tidal stripping as a formation scenario for NGC1052-DF2, an ultra diffuse galaxy significantly lacking DM, and find that its expected DM mass could be consistent with observational constraints if its progenitor was accreted early enough, z > 1.5, although it should still be a relatively rare object. △ Less

Submitted 4 February, 2021; originally announced February 2021.

Comments: 15 pages, 11 figures, 1 table, accepted for publication in MNRAS

Void Engineering in Epitaxially Regrown GaAs-Based Photonic Crystal Surface Emitting Lasers by Grating Profile Design

Authors: Adam F. McKenzie, Ben C. King, Katherine J. Rae, Stephen Thoms, Neil D. Gerrard, Jonathan Orchard, Kenishi Nishi, Keizo Takemasa, Mitsuru Sugawara, Richard J. E. Taylor, David T. D. Childs, Donald A. McLaren, Richard A. Hogg

Abstract: We report the engineering of air-voids embedded in GaAs-based photonic crystal surface emitting lasers realised by metalorganic vapour-phase epitaxy regrowth. Two distinct void geometries are obtained by modifying the photonic crystal grating profile within the reactor prior to regrowth. The mechanism of void formation is inferred from scanning transmission electron microscopy analysis, with the e… ▽ More We report the engineering of air-voids embedded in GaAs-based photonic crystal surface emitting lasers realised by metalorganic vapour-phase epitaxy regrowth. Two distinct void geometries are obtained by modifying the photonic crystal grating profile within the reactor prior to regrowth. The mechanism of void formation is inferred from scanning transmission electron microscopy analysis, with the evolution of the growth front illustrated though the use of an AlAs/GaAs superlattice structure. Competition between rapid lateral growth of the (100) surface and slow diffusion across higher index planes is exploited in order to increase void volume, leading to an order of magnitude reduction in threshold current and an increase in output power through an increase in the associated grating coupling strength. △ Less

Submitted 9 November, 2020; originally announced November 2020.

Comments: 7 pages, 4 figures

Coherent Power Scaling in Photonic Crystal Surface Emitting Laser Arrays

Authors: Ben C. King, Katherine J. Rae, Adam F. McKenzie, Aleksander Boldin, Daehyun Kim, Neil D. Gerrard, Guangrui Li, Kenishi Nishi, Mitsuru Sugawara, Richard J. E. Taylor, David T. D. Childs, Richard A. Hogg

Abstract: A key benefit of photonic crystal surface emitting lasers (PCSELs) is the abillity to increase output power through scaling the emission area while mainting high quality single mode emission, allowing them to close the brightness gap which exists between semiconductor lasers and gas and fibre lasers. However, there are practical limits to the size, and hence power, of an individual PCSEL device an… ▽ More A key benefit of photonic crystal surface emitting lasers (PCSELs) is the abillity to increase output power through scaling the emission area while mainting high quality single mode emission, allowing them to close the brightness gap which exists between semiconductor lasers and gas and fibre lasers. However, there are practical limits to the size, and hence power, of an individual PCSEL device and there are trade-offs between single-mode stability and parasitic in-plane losses with increasing device size. In this paper we discuss 2D coherent arrays as an approach to area and coherent power scaling of PCSELs. We demonstrate in two and three element PCSEL arrays an increase in the differential efficiency of the system due to a reduction in in-plane loss. △ Less

Submitted 9 November, 2020; originally announced November 2020.

Comments: Paper submitted. 9 pages, 4 figures

Identifying and interpreting tuning dimensions in deep networks

Authors: Nolan S. Dey, J. Eric Taylor, Bryan P. Tripp, Alexander Wong, Graham W. Taylor

Abstract: In neuroscience, a tuning dimension is a stimulus attribute that accounts for much of the activation variance of a group of neurons. These are commonly used to decipher the responses of such groups. While researchers have attempted to manually identify an analogue to these tuning dimensions in deep neural networks, we are unaware of an automatic way to discover them. This work contributes an unsup… ▽ More In neuroscience, a tuning dimension is a stimulus attribute that accounts for much of the activation variance of a group of neurons. These are commonly used to decipher the responses of such groups. While researchers have attempted to manually identify an analogue to these tuning dimensions in deep neural networks, we are unaware of an automatic way to discover them. This work contributes an unsupervised framework for identifying and interpreting "tuning dimensions" in deep networks. Our method correctly identifies the tuning dimensions of a synthetic Gabor filter bank and tuning dimensions of the first two layers of InceptionV1 trained on ImageNet. △ Less

Submitted 7 December, 2020; v1 submitted 5 November, 2020; originally announced November 2020.

Comments: 15 pages, 12 figures, Camera-ready for Shared Visual Representations in Human & Machine Intelligence NeurIPS Workshop 2020

ACM Class: I.2.10

1.5 μm Epitaxially Regrown Photonic Crystal Surface Emitting Laser Diode

Authors: Zijun Bian, Katherine J. Rae, Adam F. McKenzie, Ben C. King, Nasser Babazadeh, Guangrui Li, Jonathan R. Orchard, Neil D. Gerrard, Stephen Thoms, Donald A. McLaren, Richard J. E. Taylor, David Childs, Richard A. Hogg

Abstract: We present an InP-based epitaxially regrown photonic crystal surface emitting laser diode, lasing in quasi- CW conditions at 1523nm. We present an InP-based epitaxially regrown photonic crystal surface emitting laser diode, lasing in quasi- CW conditions at 1523nm. △ Less

Submitted 19 October, 2020; originally announced October 2020.

Comments: 4 pages, 5 figures, journal submission for review

Mass loss in tidally stripped systems; the energy-based truncation method

Authors: Nicole E. Drakos, James E. Taylor, Andrew J. Benson

Abstract: The ability to accurately predict the evolution of tidally stripped haloes is important for understanding galaxy formation and testing the properties of dark matter. Most studies of substructure evolution make predictions based on empirical models of tidal mass loss that are calibrated using numerical simulations. This approach can be accurate in the cases considered, but lacks generality and does… ▽ More The ability to accurately predict the evolution of tidally stripped haloes is important for understanding galaxy formation and testing the properties of dark matter. Most studies of substructure evolution make predictions based on empirical models of tidal mass loss that are calibrated using numerical simulations. This approach can be accurate in the cases considered, but lacks generality and does not provide a physical understanding of the processes involved. Recently, we demonstrated that truncating NFW distribution functions sharply in energy results in density profiles that resemble those of tidally stripped systems, offering a path to constructing physically motivated models of tidal mass loss. In this work, we review calculations of mass loss based on energy truncation alone, and then consider what secondary effects may modulate mass loss beyond this. We find that a combination of dependence on additional orbital parameters and variations in individual particle energies over an orbit results in a less abrupt truncation in energy space as a subhalo loses mass. Combining the energy truncation approach with a simple prediction for the mass-loss rate, we construct a full model of mass loss that can accurately predict the evolution of a subhalo in terms of a single parameter $η_{\rm eff}$. This parameter can be fully determined from the initial orbital and halo properties, and does not require calibration with numerical simulations. △ Less

Submitted 20 March, 2020; originally announced March 2020.

Comments: 20 pages, 14 figures, accepted to MNRAS

The Assembly of the Virgo Cluster, Traced by its Galaxy Halos

Authors: James E. Taylor, Jihye Shin, Nathalie N. -Q. Ouellette, Stéphane Courteau

Abstract: Kinematic studies have produced accurate measurements of the total dark matter mass and mean dark matter density within the optical extent of galaxies, for large samples of objects. Here we consider theoretical predictions for the latter quantity, $\barρ_{dm}$, measured within the isophotal radius $R_{23.5}$, for isolated halos with universal density profiles. Through a combination of empirical sc… ▽ More Kinematic studies have produced accurate measurements of the total dark matter mass and mean dark matter density within the optical extent of galaxies, for large samples of objects. Here we consider theoretical predictions for the latter quantity, $\barρ_{dm}$, measured within the isophotal radius $R_{23.5}$, for isolated halos with universal density profiles. Through a combination of empirical scaling relations, we show that $\barρ_{dm}$ is expected to depend weakly on halo mass and strongly on redshift. When galaxy halos fall into larger groups or clusters they become tidally stripped, reducing their total dark matter mass, but this process is expected to preserve central density until an object is close to disruption. We confirm this with collisonless simulations of cluster formation, finding that subhalos have values of $\barρ_{dm}$ close to the analytic predictions. This suggests that $\barρ_{dm}$ may be a useful indicator of infall redshift onto the cluster. We test this hypothesis with data from the SHIVir survey, which covers a reasonable fraction of the Virgo cluster. We find that galaxies with high $\barρ_{dm}$ do indeed trace the densest regions of the cluster, with a few notable exceptions. Samples selected by environment have higher densities at a significance of 3.5-4$σ$, while samples selected by density are more clustered at 3-3.5$σ$ significance. We conclude that halo density can be a powerful tracer of the assembly history of clusters and their member galaxies. △ Less

Submitted 18 June, 2019; originally announced June 2019.

Comments: 17 pages, 18 figures; accepted for publication in MNRAS

Major mergers between dark matter haloes -- II. Profile and concentration changes

Authors: Nicole E. Drakos, James E. Taylor, Anael Berrouet, Aaron S. G. Robotham, Chris Power

Abstract: Several lines of evidence suggest that as dark matter haloes grow their scale radius increases, and that the density in their central region drops. Major mergers seem an obvious mechanism to explain both these phenomena, and the resulting patterns in the concentration--mass--redshift relation. To test this possibility, we have simulated equal-mass mergers between haloes with a variety of cosmologi… ▽ More Several lines of evidence suggest that as dark matter haloes grow their scale radius increases, and that the density in their central region drops. Major mergers seem an obvious mechanism to explain both these phenomena, and the resulting patterns in the concentration--mass--redshift relation. To test this possibility, we have simulated equal-mass mergers between haloes with a variety of cosmological density profiles, placed on various different orbits. The remnants typically have higher densities than the initial conditions, but differ only slightly from self-similar scaling predictions. They are reasonably well fit by Einasto profiles, but have parameters distinct from those of the initial conditions. The net internal energy available to the merger remnant, relative to the internal energy of the initial conditions, $κ$, has the greatest influence on the properties of the final mass distribution. As expected, energetic encounters produce more extended remnants while mergers of strongly bound systems produce compact remnants. Surprisingly, however, the scale radius of the density profile shows the opposite trend, {\it increasing} in the remnants of low-energy encounters relative to energetic ones. Also even in the most energetic encounters, the density within the scale radius decreases only slightly (by 10--20\%), while for very low-energy systems it increases significantly after the merger. We conclude that while major mergers can produce remnants that are more diffuse at large radii, they are relatively ineffective at changing the central densities of haloes, and seem unlikely to explain the mean trends in the concentration--mass--redshift relation. △ Less

Submitted 27 May, 2019; v1 submitted 30 November, 2018; originally announced November 2018.

Comments: 19 pages, 16 figures; accepted for publication, MNRAS

Major mergers between dark matter haloes -- I. Predictions for size, shape, and spin

Authors: Nicole E. Drakos, James E. Taylor, Anael Berrouet, Aaron S. G. Robotham, Chris Power

Abstract: The structural properties of individual dark matter haloes, including shape, spin, concentration, and substructure, are linked to the halo's growth history, but the exact connection between the two is unclear. One open question, in particular, is the effect of major mergers on halo structure. We have performed a large set of simulations of binary equal-mass mergers between isolated haloes with var… ▽ More The structural properties of individual dark matter haloes, including shape, spin, concentration, and substructure, are linked to the halo's growth history, but the exact connection between the two is unclear. One open question, in particular, is the effect of major mergers on halo structure. We have performed a large set of simulations of binary equal-mass mergers between isolated haloes with various density profiles, to map out the relationship between the initial conditions and merger parameters and the structure of the final remnant. In this paper we describe our initial set-up and analysis methods, and report on the results for the size, shape, and spin of the merger remnant. The outcomes of mergers are most easily understood in terms of a scaled dimensionless energy parameter $κ$ and an angular momentum (or spin) parameter $λ$. We find that the axis ratio $c/a$ scales roughly linearly with energy $κ$ while the axis ratio $c/b$ scales linearly with spin $λ$. Qualitatively, mergers on radial orbits produce prolate remnants, while mergers on tangential orbits produce oblate remnants. The spin of the remnant can be predicted from angular momentum conservation, while the overall size changes as $\sim κ^{-5}$, as expected from self-similar scaling at constant mean density. We discuss potential cosmological applications for these simple patterns. △ Less

Submitted 27 May, 2019; v1 submitted 30 November, 2018; originally announced November 2018.

Comments: 18 pages, 20 figures; accepted for publication, MNRAS

Quantifying the abundance of faint, low-redshift satellite galaxies in the COSMOS survey

Authors: ChengYu Xi, James E. Taylor, Richard J. Massey, Jason Rhodes, Anton Koekemoer, Mara Salvato

Abstract: Faint dwarf satellite galaxies are important as tracers of small-scale structure, but remain poorly characterized outside the Local Group, due to the difficulty of identifying them consistently at larger distances. We review a recently proposed method for estimating the average satellite population around a given sample of nearby bright galaxies, using a combination of size and magnitude cuts (to… ▽ More Faint dwarf satellite galaxies are important as tracers of small-scale structure, but remain poorly characterized outside the Local Group, due to the difficulty of identifying them consistently at larger distances. We review a recently proposed method for estimating the average satellite population around a given sample of nearby bright galaxies, using a combination of size and magnitude cuts (to select low-redshift dwarf galaxies preferentially) and clustering measurements (to estimate the fraction of true satellites in the cut sample). We test this method using the high-precision photometric redshift catalog of the COSMOS survey, exploring the effect of specific cuts on the clustering signal. The most effective of the size-magnitude cuts considered recover the clustering signal around low-redshift primaries (z < 0.15) with about two-thirds of the signal and 80\%\ of the signal-to-noise ratio obtainable using the full COSMOS photometric redshifts. These cuts are also fairly efficient, with more than one third of the selected objects being clustered satellites. We conclude that structural selection represents a useful tool in characterizing dwarf populations to fainter magnitudes and/or over larger areas than are feasible with spectroscopic surveys.} In reviewing the low-redshift content of the COSMOS field, we also note the existence of several dozen objects that appear resolved or partially resolved in the HST imaging, and are confirmed to be local (at distances of $\sim$250 Mpc or less) by their photometric or spectroscopic redshifts. This underlines the potential for future space-based surveys to reveal local populations of intrinsically faint galaxies through imaging alone. △ Less

Submitted 18 May, 2018; originally announced May 2018.

Comments: 22 pages, 21 figures; accepted for publication in MNRAS

Searching for dark matter annihilation from individual halos: uncertainties, scatter and signal-to-noise ratios

Authors: Chiamaka Okoli, James E. Taylor, Niayesh Afshordi

Abstract: Individual extragalactic dark matter halos, such as those associated with nearby galaxies and galaxy clusters, are promising targets for searches for gamma-rays from dark matter annihilation. We review the predictions for the annihilation flux from individual halos, focusing on the effect of current uncertainties in the concentration-mass relation and the contribution from halo substructure, and a… ▽ More Individual extragalactic dark matter halos, such as those associated with nearby galaxies and galaxy clusters, are promising targets for searches for gamma-rays from dark matter annihilation. We review the predictions for the annihilation flux from individual halos, focusing on the effect of current uncertainties in the concentration-mass relation and the contribution from halo substructure, and also estimating the intrinsic halo-to-halo scatter expected. After careful consideration of recent simulation results, we conclude that the concentrations of the smallest halos, while well determined at high redshift, are still uncertain by a factor of 4-6 when extrapolated to low redshift. This in turn produces up to two orders of magnitude uncertainty in the predicted annihilation flux for any halo mass above this scale. Substructure evolution, the small-scale cutoff to the power spectrum, cosmology, and baryonic effects all introduce smaller, though cumulative, uncertainties. We then consider intrinsic variations from halo to halo. These arise from variations in concentration and substructure, leading to a scatter of $\sim$ 2.5 in the predicted annihilation luminosity. Finally, we consider the problem of detecting gamma-rays from annihilation, given the expected contributions from other sources. We estimate the signal-to-noise ratio for gamma-ray detection as a function of halo mass, assuming that cosmic rays from star formation are the main noise source in the detection. This calculation suggests that group-scale halos, individually or in stacks, may be a particularly interesting target for the next generation of annihilation searches. △ Less

Submitted 15 August, 2018; v1 submitted 14 November, 2017; originally announced November 2017.

Comments: 30 pages, 9 figures, updated to published version

Journal ref: JCAP08(2018)019

The Spectroscopy and H-band Imaging of Virgo cluster galaxies (SHIVir) Survey: Scaling Relations and the Stellar-to-Total Mass Relation

Authors: Nathalie N. -Q. Ouellette, Stéphane Courteau, Jon A. Holtzman, Aaron A. Dutton, Michele Cappellari, Julianne J. Dalcanton, Michael McDonald, Joel C. Roediger, James E. Taylor, R. Brent Tully, Patrick Côté, Laura Ferrarese, Eric W. Peng

Abstract: We present here parameter distributions and fundamental scaling relations for 190 galaxies as part of the Spectroscopy and H-bang Imaging of Virgo cluster galaxies (SHIVir) survey. We find the distribution of galaxy velocities to be bimodal about $V_{\rm circ} \sim 125$ km ${\rm s^{-1}}$, hinting at the existence of dynamically unstable modes in the inner regions of galaxies. An analysis of the Tu… ▽ More We present here parameter distributions and fundamental scaling relations for 190 galaxies as part of the Spectroscopy and H-bang Imaging of Virgo cluster galaxies (SHIVir) survey. We find the distribution of galaxy velocities to be bimodal about $V_{\rm circ} \sim 125$ km ${\rm s^{-1}}$, hinting at the existence of dynamically unstable modes in the inner regions of galaxies. An analysis of the Tully-Fisher relation (TFR) of late-type galaxies (LTGs) and fundamental plane (FP) of early-type galaxies (ETGs) is also presented, yielding a compendium of galaxy scaling relations. The slope and zero-point of the Virgo TFR match those of field galaxies, while scatter differences likely reflect distinct evolutionary histories. The velocities minimizing scatter for the TFR and FP are measured at large apertures where the baryonic fraction becomes subdominant. While TFR residuals remain independent of any galaxy parameters, FP residuals (i.e. the FP "tilt") correlate strongly with the dynamical-to-stellar mass ratio, yielding stringent galaxy formation constraints. Furthermore, we construct a stellar-to-total mass relation (STMR) for ETGs and LTGs and find linear but distinct trends over the range $M_{*} = 10^{8-11} M_{\odot}$. Stellar-to-halo mass relations (SHMRs), which probe the extended dark matter halo, can be scaled down to masses estimated within the optical radius, showing a tight match with the Virgo STMR at low masses; however, possibly inadequate halo abundance matching prescriptions and broad radial scalings complicate this comparison at all masses. While ETGs appear to be more compact than LTGs of the same stellar mass in projected space, their mass-size relations in physical space are identical. The trends reported here call for validation through well-resolved numerical simulations. △ Less

Submitted 30 May, 2017; originally announced May 2017.

Comments: Accepted to ApJ. 28 pages, 16 figures

The phase-space structure of tidally stripped halos

Authors: Nicole E. Drakos, James E. Taylor, Andrew J. Benson

Abstract: We propose a new method for generating equilibrium models of spherical systems of collisionless particles that are finite in extent, but whose central regions resemble dark matter halos from cosmological simulations. This method involves iteratively removing unbound particles from a Navarro-Frenk-White profile truncated sharply at some radius. The resulting models are extremely stable, and thus pr… ▽ More We propose a new method for generating equilibrium models of spherical systems of collisionless particles that are finite in extent, but whose central regions resemble dark matter halos from cosmological simulations. This method involves iteratively removing unbound particles from a Navarro-Frenk-White profile truncated sharply at some radius. The resulting models are extremely stable, and thus provide a good starting point for N-body simulations of isolated halos. We provide a code to generate such models for NFW and a variety of other common density profiles. We then develop an analytic approximation to this truncated distribution function. Our method proceeds by analogy with the King model, truncating and shifting the original distribution function of an infinitely extended Navarro-Frenk-White profile in energy space. We show that the density profiles of our models closely resemble the tidally truncated density profiles seen previously in studies of satellite evolution. Pursuing this analogy further with a series of simulations of tidal mass loss, we find that our models provide a good approximation to the full distribution function of tidally stripped systems, thus allowing theoretically motivated phase-space calculations for such systems. △ Less

Submitted 22 March, 2017; originally announced March 2017.

The Next Generation Virgo Cluster Survey (NGVS). XXIV. The Red Sequence to $\sim$10$^6$ L$_{\odot}$ and Comparisons with Galaxy Formation Models

Authors: Joel C. Roediger, Laura Ferrarese, Patrick Côté, Lauren A. MacArthur, Rúben Sánchez-Janssen, John P. Blakeslee, Eric W. Peng, Chengze Liu, Roberto Munoz, Jean-Charles Cuillandre, Stephen Gwyn, Simona Mei, Samuel Boissier, Alessandro Boselli, Michele Cantiello, Stéphane Courteau, Pierre-Alain Duc, Ariane Lançon, J. Christopher Mihos, Thomas H. Puzia, James E. Taylor, Patrick R. Durrell, Elisa Toloba, Puragra Guhathakurta, Hongxin Zhang

Abstract: We use deep optical photometry from the Next Generation Virgo Cluster Survey [NGVS] to investigate the color-magnitude diagram for the galaxies inhabiting the core of this cluster. The sensitivity of the NGVS imaging allows us to continuously probe galaxy colors over a factor of $\sim 2 \times 10^5$ in luminosity, from brightest cluster galaxies to scales overlapping classical satellites of the Mi… ▽ More We use deep optical photometry from the Next Generation Virgo Cluster Survey [NGVS] to investigate the color-magnitude diagram for the galaxies inhabiting the core of this cluster. The sensitivity of the NGVS imaging allows us to continuously probe galaxy colors over a factor of $\sim 2 \times 10^5$ in luminosity, from brightest cluster galaxies to scales overlapping classical satellites of the Milky Way [$M_{g^{\prime}}$ $\sim$ $-$9; $M_{*}$ $\sim 10^6$ M$_{\odot}$], within a single environment. Remarkably, we find the first evidence that the RS flattens in all colors at the faint-magnitude end [starting between $-$14 $\le$ $M_{g^{\prime}}$ $\le$ $-$13, around $M_{*}$ $\sim 4 \times 10^7$ M$_{\odot}$], with the slope decreasing to $\sim$60% or less of its value at brighter magnitudes. This could indicate that the stellar populations of faint dwarfs in Virgo's core share similar characteristics [e.g. constant mean age] over $\sim$3 mags in luminosity, suggesting that these galaxies were quenched coevally, likely via pre-processing in smaller hosts. We also compare our results to galaxy formation models, finding that the RS in model clusters have slopes at intermediate magnitudes that are too shallow, and in the case of semi-analytic models, do not reproduce the flattening seen at both extremes [bright/faint] of the Virgo RS. Deficiencies in the chemical evolution of model galaxies likely contribute to the model-data discrepancies at all masses, while overly efficient quenching may also be a factor at dwarf scales. Deep UV and near-IR photometry are required to unambiguously diagnose the cause of the faint-end flattening. △ Less

Submitted 28 October, 2016; originally announced October 2016.

Comments: 19 pages, 11 figures, accepted to ApJ

Urban Energy Flux: Human Mobility as a Predictor for Spatial Changes

Authors: Neda Mohammadi, John E Taylor

Abstract: As a key energy challenge, we urgently require a better understanding of how growing urban populations interact with municipal energy systems and the resulting impact on energy demand across city neighborhoods, which are dense hubs of both consumer population and CO2 emissions. Currently, the physical characteristics of urban infrastructure are the main determinants in predictive modeling of the d… ▽ More As a key energy challenge, we urgently require a better understanding of how growing urban populations interact with municipal energy systems and the resulting impact on energy demand across city neighborhoods, which are dense hubs of both consumer population and CO2 emissions. Currently, the physical characteristics of urban infrastructure are the main determinants in predictive modeling of the demand side of energy in our rapidly growing urban areas; overlooking influence related to fluctuating human activities. Here, we show how applying intra-urban human mobility as an indicator for interactions of the population with local energy systems can be translated into spatial imprints to predict the spatial distribution of energy use in urban settings. Our findings establish human mobility as an important element in explaining the spatial structure underlying urban energy flux and demonstrate the utility of a human mobility driven approach for predicting future urban energy demand with implications for CO2 emission strategies. △ Less

Submitted 5 September, 2016; originally announced September 2016.

The Next Generation Virgo Cluster Survey (NGVS). XIII. The Luminosity and Mass Function of Galaxies in the Core of the Virgo Cluster and the Contribution from Disrupted Satellites

Authors: Laura Ferrarese, Patrick Cote, Ruben Sanchez-Janssen, Joel Roediger, Alan W. McConnachie, Patrick R. Durrell, Lauren A. MacArthur, John P. Blakeslee, Pierre-Alain Duc, S. Boissier, Alessandro Boselli, Stephane Courteau, Jean-Charles Cuillandre, Eric Emsellem, S. D. J. Gwyn, Puragra Guhathakurta, Andres Jordan, Ariane Lancon, Chengze Liu, Simona Mei, J. Christopher Mihos, Thomas H. Puzia, James E. Taylor, Hongxin Zhang

Abstract: We present measurements of the galaxy luminosity and stellar mass function in a 3.71 deg$^2$ (0.3 Mpc$^2$) area in the core of the Virgo cluster, based on $ugriz$ data from the Next Generation Virgo Cluster Survey (NGVS). The galaxy sample consists of 352 objects brighter than $M_g=-9.13$ mag, the 50% completeness limit of the survey. Using a Bayesian analysis, we find a best-fit faint end slope o… ▽ More We present measurements of the galaxy luminosity and stellar mass function in a 3.71 deg$^2$ (0.3 Mpc$^2$) area in the core of the Virgo cluster, based on $ugriz$ data from the Next Generation Virgo Cluster Survey (NGVS). The galaxy sample consists of 352 objects brighter than $M_g=-9.13$ mag, the 50% completeness limit of the survey. Using a Bayesian analysis, we find a best-fit faint end slope of $α=-1.33 \pm 0.02$ for the g-band luminosity function; consistent results are found for the stellar mass function as well as the luminosity function in the other four NGVS bandpasses. We discuss the implications for the faint-end slope of adding 92 ultra compact dwarfs galaxies (UCDs) -- previously compiled by the NGVS in this region -- to the galaxy sample, assuming that UCDs are the stripped remnants of nucleated dwarf galaxies. Under this assumption, the slope of the luminosity function (down to the UCD faint magnitude limit, $M_g = -9.6$ mag) increases dramatically, up to $α= -1.60 \pm 0.06$ when correcting for the expected number of disrupted non-nucleated galaxies. We also calculate the total number of UCDs and globular clusters that may have been deposited in the core of Virgo due to the disruption of satellites, both nucleated and non-nucleated. We estimate that ~150 objects with $M_g\lesssim-9.6$ mag and that are currently classified as globular clusters, might, in fact, be the nuclei of disrupted galaxies. We further estimate that as many as 40% of the (mostly blue) globular clusters in the core of Virgo might once have belonged to such satellites; these same disrupted satellites might have contributed ~40% of the total luminosity in galaxies observed in the core region today. Finally, we use an updated Local Group galaxy catalog to provide a new measurement of the luminosity function of Local Group satellites, $α=-1.21\pm0.05$. △ Less

Submitted 21 April, 2016; originally announced April 2016.

Comments: ApJ, accepted

The Intrinsic Geometry of Some Random Manifolds

Authors: Sunder Ram Krishnan, Jonathan E. Taylor, Robert J. Adler

Abstract: We study the a.s. convergence of a sequence of random embeddings of a fixed manifold into Euclidean spaces of increasing dimensions. We show that the limit is deterministic. As a consequence, we show that many intrinsic functionals of the embedded manifolds also converge to deterministic limits. Particularly interesting examples of these functionals are given by the Lipschitz-Killing curvatures, f… ▽ More We study the a.s. convergence of a sequence of random embeddings of a fixed manifold into Euclidean spaces of increasing dimensions. We show that the limit is deterministic. As a consequence, we show that many intrinsic functionals of the embedded manifolds also converge to deterministic limits. Particularly interesting examples of these functionals are given by the Lipschitz-Killing curvatures, for which we also prove unbiasedness, using the Gaussian kinematic formula. △ Less

Submitted 17 December, 2015; originally announced December 2015.

Comments: 11 pages

MSC Class: 60G15; 57N35; 60D05; 60G60; 70G45

LoCuSS: Testing hydrostatic equilibrium in galaxy clusters

Authors: G. P. Smith, P. Mazzotta, N. Okabe, F. Ziparo, S. L. Mulroy, A. Babul, A. Finoguenov, I. G. McCarthy, M. Lieu, Y. Bahe, H. Bourdin, A. E. Evrard, T. Futamase, C. P. Haines, M. Jauzac, D. P. Marrone, R. Martino, P. E. May, J. E. Taylor, K. Umetsu

Abstract: We test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at $0.15<z<0.3$ from the Local Cluster Substructure Survey (LoCuSS). Our weak-lensing measurements of $M_{500}$ control systematic biases to sub-4 per cent, and our hydrostatic measurements of the same achieve excellent agreement between XMM-Newton and Chandra. The mean ratio of X-ray to… ▽ More We test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at $0.15<z<0.3$ from the Local Cluster Substructure Survey (LoCuSS). Our weak-lensing measurements of $M_{500}$ control systematic biases to sub-4 per cent, and our hydrostatic measurements of the same achieve excellent agreement between XMM-Newton and Chandra. The mean ratio of X-ray to lensing mass for these 50 clusters is $β_{\rm X}=0.95\pm0.05$, and for the 44 clusters also detected by Planck, the mean ratio of Planck mass estimate to LoCuSS lensing mass is $β_{\rm P}=0.95\pm0.04$. Based on a careful like-for-like analysis, we find that LoCuSS, the Canadian Cluster Comparison Project (CCCP), and Weighing the Giants (WtG) agree on $β_{\rm P}\simeq0.9-0.95$ at $0.15<z<0.3$. This small level of hydrostatic bias disagrees at $\sim5σ$ with the level required to reconcile Planck cosmology results from the cosmic microwave background and galaxy cluster counts. △ Less

Submitted 5 November, 2015; originally announced November 2015.

Comments: 5 pages, 1 figure, accepted by MNRAS

Selective inference in regression models with groups of variables

Authors: Joshua R. Loftus, Jonathan E. Taylor

Abstract: We provide a general mathematical framework for selective inference with supervised model selection procedures characterized by quadratic forms in the outcome variable. Forward stepwise with groups of variables is an important special case as it allows models with categorical variables or factors. Models can be chosen by AIC, BIC, or a fixed number of steps. We provide an exact significance test f… ▽ More We provide a general mathematical framework for selective inference with supervised model selection procedures characterized by quadratic forms in the outcome variable. Forward stepwise with groups of variables is an important special case as it allows models with categorical variables or factors. Models can be chosen by AIC, BIC, or a fixed number of steps. We provide an exact significance test for each group of variables in the selected model based on an appropriately truncated $χ$ or $F$ distribution for the cases of known and unknown $σ^2$ respectively. An efficient software implementation is available as a package in the R statistical programming language. △ Less

Submitted 4 November, 2015; originally announced November 2015.

Comments: 16 pages, 4 figures

MSC Class: 62J05

Selective inference with a randomized response

Authors: Xiaoying Tian, Jonathan E. Taylor

Abstract: Inspired by sample splitting and the reusable holdout introduced in the field of differential privacy, we consider selective inference with a randomized response. We discuss two major advantages of using a randomized response for model selection. First, the selectively valid tests are more powerful after randomized selection. Second, it allows consistent estimation and weak convergence of selectiv… ▽ More Inspired by sample splitting and the reusable holdout introduced in the field of differential privacy, we consider selective inference with a randomized response. We discuss two major advantages of using a randomized response for model selection. First, the selectively valid tests are more powerful after randomized selection. Second, it allows consistent estimation and weak convergence of selective inference procedures. Under independent sampling, we prove a selective (or privatized) central limit theorem that transfers procedures valid under asymptotic normality without selection to their corresponding selective counterparts. This allows selective inference in nonparametric settings. Finally, we propose a framework of inference after combining multiple randomized selection procedures. We focus on the classical asymptotic setting, leaving the interesting high-dimensional asymptotic questions for future work. △ Less

Submitted 30 November, 2016; v1 submitted 24 July, 2015; originally announced July 2015.

The Next Generation Virgo Cluster Survey. IX. Estimating the Efficiency of Galaxy Formation on the Lowest-Mass Scales

Authors: Jonathan Grossauer, James E. Taylor, Laura Ferrarese, Lauren A. MacArthur, Patrick Cote, Joel Roediger, Stephane Courteau, Jean-Charles Cuillandre, Pierre-Alain Duc, Patrick R. Durrell, S. D. J. Gwyn, Andres Jordan, Simona Mei, Eric W. Peng

Abstract: The Next Generation Virgo Cluster Survey has recently determined the luminosity function of galaxies in the core of the Virgo cluster down to unprecedented magnitude and surface brightness limits. Comparing simulations of cluster formation to the derived central stellar mass function, we attempt to estimate the stellar-to-halo-mass ratio (SHMR) for dwarf galaxies, as it would have been before they… ▽ More The Next Generation Virgo Cluster Survey has recently determined the luminosity function of galaxies in the core of the Virgo cluster down to unprecedented magnitude and surface brightness limits. Comparing simulations of cluster formation to the derived central stellar mass function, we attempt to estimate the stellar-to-halo-mass ratio (SHMR) for dwarf galaxies, as it would have been before they fell into the cluster. This approach ignores several details and complications, e.g., the contribution of ongoing star formation to the present-day stellar mass of cluster members, and the effects of adiabatic contraction and/or violent feedback on the subhalo and cluster potentials. The final results are startlingly simple, however; we find that the trends in the SHMR determined previously for bright galaxies appear to extend down in a scale-invariant way to the faintest objects detected in the survey. These results extend measurements of the formation efficiency of field galaxies by two decades in halo mass, or five decades in stellar mass, down to some of the least massive dwarf galaxies known, with stellar masses of $\sim 10^5 M_\odot$. △ Less

Submitted 9 July, 2015; originally announced July 2015.

Comments: 18 pages, 12 figures; published in ApJ July 1st 2015

Journal ref: ApJ 807,88 (2015)

Selective inference with unknown variance via the square-root LASSO

Authors: Xiaoying Tian, Joshua R. Loftus, Jonathan E. Taylor

Abstract: There has been much recent work on inference after model selection when the noise level is known, however, $σ$ is rarely known in practice and its estimation is difficult in high-dimensional settings. In this work we propose using the square-root LASSO (also known as the scaled LASSO) to perform selective inference for the coefficients and the noise level simultaneously. The square-root LASSO has… ▽ More There has been much recent work on inference after model selection when the noise level is known, however, $σ$ is rarely known in practice and its estimation is difficult in high-dimensional settings. In this work we propose using the square-root LASSO (also known as the scaled LASSO) to perform selective inference for the coefficients and the noise level simultaneously. The square-root LASSO has the property that choosing a reasonable tuning parameter is scale-free, namely it does not depend on the noise level in the data. We provide valid p-values and confidence intervals for the coefficients after selection, and estimates for model specific variance. Our estimates perform better than other estimates of $σ^2$ in simulation. △ Less

Submitted 9 February, 2017; v1 submitted 29 April, 2015; originally announced April 2015.

Communication-efficient sparse regression: a one-shot approach

Authors: Jason D. Lee, Yuekai Sun, Qiang Liu, Jonathan E. Taylor

Abstract: We devise a one-shot approach to distributed sparse regression in the high-dimensional setting. The key idea is to average "debiased" or "desparsified" lasso estimators. We show the approach converges at the same rate as the lasso as long as the dataset is not split across too many machines. We also extend the approach to generalized linear models. We devise a one-shot approach to distributed sparse regression in the high-dimensional setting. The key idea is to average "debiased" or "desparsified" lasso estimators. We show the approach converges at the same rate as the lasso as long as the dataset is not split across too many machines. We also extend the approach to generalized linear models. △ Less

Submitted 11 August, 2015; v1 submitted 14 March, 2015; originally announced March 2015.

Comments: 29 pages, 3 figures

Convergence of the Reach for a Sequence of Gaussian-Embedded Manifolds

Authors: Robert J. Adler, Sunder Ram Krishnan, Jonathan E. Taylor, Shmuel Weinberger

Abstract: Motivated by questions of manifold learning, we study a sequence of random manifolds, generated by embedding a fixed, compact manifold $M$ into Euclidean spheres of increasing dimension via a sequence of Gaussian mappings. One of the fundamental smoothness parameters of manifold learning theorems is the reach, or critical radius, of $M$. Roughly speaking, the reach is a measure of a manifold's dep… ▽ More Motivated by questions of manifold learning, we study a sequence of random manifolds, generated by embedding a fixed, compact manifold $M$ into Euclidean spheres of increasing dimension via a sequence of Gaussian mappings. One of the fundamental smoothness parameters of manifold learning theorems is the reach, or critical radius, of $M$. Roughly speaking, the reach is a measure of a manifold's departure from convexity, which incorporates both local curvature and global topology. This paper develops limit theory for the reach of a family of random, Gaussian-embedded, manifolds, establishing both almost sure convergence for the global reach, and a fluctuation theory for both it and its local version. The global reach converges to a constant well known both in the reproducing kernel Hilbert space theory of Gaussian processes, as well as in their extremal theory. △ Less

Submitted 7 August, 2016; v1 submitted 5 March, 2015; originally announced March 2015.

Comments: 38 pages. Removed motivational material from the previous version

MSC Class: 60G15; 57N35; 60D05; 60G60

The Mass-Concentration Relation and the Stellar-to-Halo Mass Ratio in the CFHT Stripe 82 Survey

Authors: Huanyuan Shan, Jean-Paul Kneib, Ran Li, Johan Comparat, Thomas Erben, Martin Makler, Bruno Moraes, Ludovic Van Waerbeke, James E. Taylor, Aldée Charbonnier, Maria E. S. Pereira

Abstract: We present a new measurement of the mass-concentration relation and the stellar-to-halo mass ratio over the halo mass range $5\times 10^{12}$ to $2\times 10^{14}M_{\odot}$. To achieve this, we use weak lensing measurements from the CFHT Stripe 82 Survey (CS82), combined with the central galaxies from the redMaPPer cluster catalogue and the LOWZ/CMASS galaxy sample of the Sloan Digital Sky Survey-I… ▽ More We present a new measurement of the mass-concentration relation and the stellar-to-halo mass ratio over the halo mass range $5\times 10^{12}$ to $2\times 10^{14}M_{\odot}$. To achieve this, we use weak lensing measurements from the CFHT Stripe 82 Survey (CS82), combined with the central galaxies from the redMaPPer cluster catalogue and the LOWZ/CMASS galaxy sample of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The stacked lensing signals around these samples are modelled as a sum of contributions from the central galaxy, its dark matter halo, and the neighboring halos, as well as a term for possible centering errors. We measure the mass-concentration relation: $c_{200c}(M)=A(\frac{M_{200c}}{M_0})^{B}$ with $A=5.24\pm1.24, B=-0.13\pm0.10$ for $0.2<z<0.4$ and $A=6.61\pm0.75, B=-0.15\pm0.05$ for $0.4<z<0.6$. These amplitudes and slopes are completely consistent with predictions from recent simulations. We also measure the stellar-to-halo mass ratio for our samples, and find results consistent with previous measurements from lensing and other techniques. △ Less

Submitted 23 September, 2017; v1 submitted 1 February, 2015; originally announced February 2015.

Comments: 10 pages, 3 figures, 3 tables

Journal ref: ApJ, 840, 104 (2017)

Brightest X-ray clusters of galaxies in the CFHTLS wide fields: Catalog and optical mass estimator

Authors: M. Mirkazemi, A. Finoguenov, M. J. Pereira, M. Tanaka, M. Lerchster, F. Brimioulle, E. Egami, K. Kettula, G. Erfanianfar, H. J. McCracken, Y. Mellier, J. P. Kneib, E. Rykoff, S. Seitz, T. Erben, J. E. Taylor

Abstract: The CFHTLS presents a unique data set for weak lensing studies, having high quality imaging and deep multi-band photometry. We have initiated an XMM-CFHTLS project to provide X-ray observations of the brightest X-ray selected clusters within the wide CFHTLS area. Performance of these observations and the high quality of CFHTLS data, allows us to revisit the identification of X-ray sources, introdu… ▽ More The CFHTLS presents a unique data set for weak lensing studies, having high quality imaging and deep multi-band photometry. We have initiated an XMM-CFHTLS project to provide X-ray observations of the brightest X-ray selected clusters within the wide CFHTLS area. Performance of these observations and the high quality of CFHTLS data, allows us to revisit the identification of X-ray sources, introducing automated reproducible algorithms, based on the multi-color red sequence finder. We have also introduced a new optical mass proxy. We provide the calibration of the red sequence observed in the CFHT filters and compare the results with the traditional single color red sequence and photoz. We test the identification algorithm on the subset of highly significant XMM clusters and identify 100% of the sample. We find that the integrated z-band luminosity of the red sequence galaxies correlates well with the X-ray luminosity with a surprisingly small scatter of 0.20 dex. We further use the multi-color red sequence to reduce spurious detections in the full XMM and RASS data sets, resulting in catalogs of 196 and 32 clusters, respectively. We made spectroscopic follow-up observations of some of these systems with HECTOSPEC and in combination with BOSS DR9 data. We also describe the modifications needed to the source detection algorithm in order to keep high purity of extended sources in the shallow X-ray data. We also present the scaling relation between X-ray luminosity and velocity dispersion. △ Less

Submitted 12 November, 2014; originally announced November 2014.

Comments: 29 pages, 20 figures, accepted for publication in ApJ

The Next Generation Virgo Cluster Survey. VIII. The Spatial Distribution of Globular Clusters in the Virgo Cluster

Authors: Patrick R. Durrell, Patrick Côté, Eric W. Peng, John P. Blakeslee, Laura Ferrarese, J. Christopher Mihos, Thomas H. Puzia, Ariane Lançon, Chengze Liu, Hongxin Zhang, Jean-Charles Cuillandre, Alan McConnachie, Andrés Jordan, Katharine Accetta, Samual Boissier, Alessandro Boselli, Stéphane Courteau, Pierre-Alain Duc, Eric Emsellem, Stephen Gwyn, Simona Mei, James E. Taylor

Abstract: We report on a large-scale study of the distribution of globular clusters (GCs) throughout the Virgo cluster, based on photometry from the Next Generation Virgo Cluster Survey, a large imaging survey covering Virgo's primary subclusters to their virial radii. Using the g', (g'-i') color-magnitude diagram of unresolved and marginally-resolved sources, we constructed 2-D maps of the GC distribution.… ▽ More We report on a large-scale study of the distribution of globular clusters (GCs) throughout the Virgo cluster, based on photometry from the Next Generation Virgo Cluster Survey, a large imaging survey covering Virgo's primary subclusters to their virial radii. Using the g', (g'-i') color-magnitude diagram of unresolved and marginally-resolved sources, we constructed 2-D maps of the GC distribution. We present the clearest evidence to date showing the difference in concentration between red and blue GCs over the extent of the cluster, where the red (metal-rich) GCs are largely located around the massive early-type galaxies, whilst the blue (metal-poor) GCs have a more extended spatial distribution, with significant populations present beyond 83' (215 kpc) along the major axes of M49 and M87. The GC distribution around M87 and M49 shows remarkable agreement with the shape, ellipticity and boxiness of the diffuse light surrounding both galaxies. We find evidence for spatial enhancements of GCs surrounding M87 that may be indicative of recent interactions or an ongoing merger history. We compare the GC map to the locations of Virgo galaxies and the intracluster X-ray gas, and find good agreement between these baryonic structures. The Virgo cluster contains a total population of 67300$\pm$14400 GCs, of which 35% are located in M87 and M49 alone. We compute a cluster-wide specific frequency S_N,CL=$2.8\pm0.7$, including Virgo's diffuse light. The GC-to-baryonic mass fraction is e_b=$5.7\pm1.1\times10^{-4} $and the GC-to-total cluster mass formation efficiency is e_t=$2.9\pm0.5\times10^{-5}$, values slightly lower than, but consistent with, those derived for individual galactic halos. Our results show that the production of the complex structures in the unrelaxed Virgo cluster core (including the diffuse intracluster light) is an ongoing process.(abridged) △ Less

Submitted 16 September, 2014; v1 submitted 12 August, 2014; originally announced August 2014.

Comments: 23 pages, 17 figures. Accepted for publication in the Astrophysical Journal. Figure 1 has reduced resolution. Revised version with updated references, corrected typos -- no changes to results

Topological consistency via kernel estimation

Authors: Omer Bobrowski, Sayan Mukherjee, Jonathan E. Taylor

Abstract: We introduce a consistent estimator for the homology (an algebraic structure representing connected components and cycles) of level sets of both density and regression functions. Our method is based on kernel estimation. We apply this procedure to two problems: (1) inferring the homology structure of manifolds from noisy observations, (2) inferring the persistent homology (a multi-scale extension… ▽ More We introduce a consistent estimator for the homology (an algebraic structure representing connected components and cycles) of level sets of both density and regression functions. Our method is based on kernel estimation. We apply this procedure to two problems: (1) inferring the homology structure of manifolds from noisy observations, (2) inferring the persistent homology (a multi-scale extension of homology) of either density or regression functions. We prove consistency for both of these problems. In addition to the theoretical results, we demonstrate these methods on simulated data for binary regression and clustering applications. △ Less

Submitted 29 September, 2016; v1 submitted 20 July, 2014; originally announced July 2014.

Comments: Published at http://dx.doi.org/10.3150/15-BEJ744 in the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm)

Report number: IMS-BEJ-BEJ744

Journal ref: Bernoulli 2017, Vol. 23, No. 1, 288-328

A significance test for forward stepwise model selection

Authors: Joshua R. Loftus, Jonathan E. Taylor

Abstract: We apply the methods developed by Lockhart et al. (2013) and Taylor et al. (2013) on significance tests for penalized regression to forward stepwise model selection. A general framework for selection procedures described by quadratic inequalities includes a variant of forward stepwise with grouped variables, allowing us to handle categorical variables and factor models. We provide an algorithm to… ▽ More We apply the methods developed by Lockhart et al. (2013) and Taylor et al. (2013) on significance tests for penalized regression to forward stepwise model selection. A general framework for selection procedures described by quadratic inequalities includes a variant of forward stepwise with grouped variables, allowing us to handle categorical variables and factor models. We provide an algorithm to compute a new statistic with an exact null distribution conditional on the outcome of the model selection procedure. This new statistic, which we denote $Tχ$, has a truncated $χ$ distribution under the global null. We apply this test in forward stepwise iteratively on the residual after each step. The resulting method has the computational strengths of stepwise selection and addresses the problem of invalid test statistics due to model selection. We illustrate the flexibility of this method by applying it to several specialized applications of forward stepwise including a hierarchical interactions model and a recently described additive model that adaptively chooses between linear and nonlinear effects for each variable. △ Less

Submitted 15 May, 2014; originally announced May 2014.

Valid post-correction inference for censored regression problems

Authors: Yuekai Sun, Jonathan E. Taylor

Abstract: Two-step estimators often called upon to fit censored regression models in many areas of science and engineering. Since censoring incurs a bias in the naive least-squares fit, a two-step estimator first estimates the bias and then fits a corrected linear model. We develop a framework for performing valid /post-correction inference/ with two-step estimators. By exploiting recent results on post-sel… ▽ More Two-step estimators often called upon to fit censored regression models in many areas of science and engineering. Since censoring incurs a bias in the naive least-squares fit, a two-step estimator first estimates the bias and then fits a corrected linear model. We develop a framework for performing valid /post-correction inference/ with two-step estimators. By exploiting recent results on post-selection inference, we obtain valid confidence intervals and significance tests for the fitted coefficients. △ Less

Submitted 13 March, 2014; originally announced March 2014.

Comments: 20 pages, 6 figues

Exact Post Model Selection Inference for Marginal Screening

Authors: Jason D Lee, Jonathan E Taylor

Abstract: We develop a framework for post model selection inference, via marginal screening, in linear regression. At the core of this framework is a result that characterizes the exact distribution of linear functions of the response $y$, conditional on the model being selected (``condition on selection" framework). This allows us to construct valid confidence intervals and hypothesis tests for regression… ▽ More We develop a framework for post model selection inference, via marginal screening, in linear regression. At the core of this framework is a result that characterizes the exact distribution of linear functions of the response $y$, conditional on the model being selected (``condition on selection" framework). This allows us to construct valid confidence intervals and hypothesis tests for regression coefficients that account for the selection procedure. In contrast to recent work in high-dimensional statistics, our results are exact (non-asymptotic) and require no eigenvalue-like assumptions on the design matrix $X$. Furthermore, the computational cost of marginal regression, constructing confidence intervals and hypothesis testing is negligible compared to the cost of linear regression, thus making our methods particularly suitable for extremely large datasets. Although we focus on marginal screening to illustrate the applicability of the condition on selection framework, this framework is much more broadly applicable. We show how to apply the proposed framework to several other selection procedures including orthogonal matching pursuit, non-negative least squares, and marginal screening+Lasso. △ Less

Submitted 27 February, 2014; v1 submitted 23 February, 2014; originally announced February 2014.

Quantifying Human Mobility Perturbation and Resilience in Natural Disasters

Authors: Qi Wang, John E. Taylor

Abstract: Human mobility is influenced by environmental change and natural disasters. Researchers have used trip distance distribution, radius of gyration of movements, and individuals' visited locations to understand and capture human mobility patterns and trajectories. However, our knowledge of human movements during natural disasters is limited owing to both a lack of empirical data and the low precision… ▽ More Human mobility is influenced by environmental change and natural disasters. Researchers have used trip distance distribution, radius of gyration of movements, and individuals' visited locations to understand and capture human mobility patterns and trajectories. However, our knowledge of human movements during natural disasters is limited owing to both a lack of empirical data and the low precision of available data. Here, we studied human mobility using high-resolution movement data from individuals in New York City during and for several days after Hurricane Sandy in 2012. We found the human movements followed truncated power-law distributions during and after Hurricane Sandy, although the β value was noticeably larger during the first 24 hours after the storm struck. Also, we examined two parameters: the center of mass and the radius of gyration of each individual's movements. We found that their values during perturbation states and steady states are highly correlated, suggesting human mobility data obtained in steady states can possibly predict the perturbation state. Our results demonstrate that human movement trajectories experienced significant perturbations during hurricanes, but also exhibited high resilience. We expect the study will stimulate future research on the perturbation and inherent resilience of human mobility under the influence of natural disasters. For example, mobility patterns in coastal urban areas could be examined as tropical cyclones approach, gain or dissipate in strength, and as the path of the storm changes. Understanding nuances of human mobility under the influence of disasters will enable more effective evacuation, emergency response planning and development of strategies and policies to reduce fatality, injury, and economic loss. △ Less

Submitted 9 February, 2014; originally announced February 2014.

Comments: 11 pages, 3 figures

Exact post-selection inference, with application to the lasso

Authors: Jason D. Lee, Dennis L. Sun, Yuekai Sun, Jonathan E. Taylor

Abstract: We develop a general approach to valid inference after model selection. At the core of our framework is a result that characterizes the distribution of a post-selection estimator conditioned on the selection event. We specialize the approach to model selection by the lasso to form valid confidence intervals for the selected coefficients and test whether all relevant variables have been included in… ▽ More We develop a general approach to valid inference after model selection. At the core of our framework is a result that characterizes the distribution of a post-selection estimator conditioned on the selection event. We specialize the approach to model selection by the lasso to form valid confidence intervals for the selected coefficients and test whether all relevant variables have been included in the model. △ Less

Submitted 3 May, 2016; v1 submitted 25 November, 2013; originally announced November 2013.

Comments: Published at http://dx.doi.org/10.1214/15-AOS1371 in the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOS-AOS1371

Journal ref: Annals of Statistics 2016, Vol. 44, No. 3, 907-927

First Measurement of the Cross-Correlation of CMB Lensing and Galaxy Lensing

Authors: Nick Hand, Alexie Leauthaud, Sudeep Das, Blake D. Sherwin, Graeme E. Addison, J. Richard Bond, Erminia Calabrese, Aldée Charbonnier, Mark J. Devlin, Joanna Dunkley, Thomas Erben, Amir Hajian, Mark Halpern, Joachim Harnois-Déraps, Catherine Heymans, Hendrik Hildebrandt, Adam D. Hincks, Jean-Paul Kneib, Arthur Kosowsky, Martin Makler, Lance Miller, Kavilan Moodley, Bruno Moraes, Michael D. Niemack, Lyman A. Page , et al. (11 additional authors not shown)

Abstract: We measure the cross-correlation of cosmic microwave background lensing convergence maps derived from Atacama Cosmology Telescope data with galaxy lensing convergence maps as measured by the Canada-France-Hawaii Telescope Stripe 82 Survey. The CMB-galaxy lensing cross power spectrum is measured for the first time with a significance of 4.2σ, which corresponds to a 12% constraint on the amplitude o… ▽ More We measure the cross-correlation of cosmic microwave background lensing convergence maps derived from Atacama Cosmology Telescope data with galaxy lensing convergence maps as measured by the Canada-France-Hawaii Telescope Stripe 82 Survey. The CMB-galaxy lensing cross power spectrum is measured for the first time with a significance of 4.2σ, which corresponds to a 12% constraint on the amplitude of density fluctuations at redshifts ~ 0.9. With upcoming improved lensing data, this novel type of measurement will become a powerful cosmological probe, providing a precise measurement of the mass distribution at intermediate redshifts and serving as a calibrator for systematic biases in weak lensing measurements. △ Less

Submitted 11 December, 2014; v1 submitted 24 November, 2013; originally announced November 2013.

Comments: 11 pages, 7 figures. Submitted to PRD. Updates to lensing maps and analysis pipeline, with revised results, increased significance

Weak lensing mass map and peak statistics in Canada-France-Hawaii Telescope Stripe 82 survey

Authors: HuanYuan Shan, Jean-Paul Kneib, Johan Comparat, Eric Jullo, Aldee Charbonnier, Thomas Erben, Martin Makler, Bruno Moraes, Ludovic Van Waerbeke, Frederic Courbin, Georges Meylan, Charling Tao, James E. Taylor

Abstract: We present a weak lensing mass map covering ~124 square degrees of the Canada-France-Hawaii Telescope Stripe 82 Survey (CS82). We study the statistics of rare peaks in the map, including peak abundance, the peak-peak correlation functions and the tangential-shear profiles around peaks. We find that the abundance of peaks detected in CS82 is consistent with predictions from a Lambda-CDM cosmologica… ▽ More We present a weak lensing mass map covering ~124 square degrees of the Canada-France-Hawaii Telescope Stripe 82 Survey (CS82). We study the statistics of rare peaks in the map, including peak abundance, the peak-peak correlation functions and the tangential-shear profiles around peaks. We find that the abundance of peaks detected in CS82 is consistent with predictions from a Lambda-CDM cosmological model, once noise effects are properly included. The correlation functions of peaks with different signal-to-noise ratio (SNR) are well described by power laws, and there is a clear cross-correlation between the Sloan Digital Sky Survey III/Constant Mass galaxies and high SNR peaks. The tangential-shear profiles around peaks increase with peak SNR. We fit analytical models to the tangential-shear profiles, including a projected singular isothermal sphere (SIS) model and a projected Navarro, Frenk & White (NFW) model, plus a two-halo term. For the high SNR peaks, the SIS model is rejected at ~3-sigma. The NFW model plus a two-halo term gives more acceptable fits to the data. Some peaks match the positions of optically detected clusters, while others are relatively dark. Comparing dark and matched peaks, we find a difference in lensing signal of a factor of 2, suggesting that about half of the dark peaks are false detections. △ Less

Submitted 18 July, 2014; v1 submitted 6 November, 2013; originally announced November 2013.

Comments: 10 pages, 12 figures, Published in MNRAS

Journal ref: MNRAS 442, 2534 (2014)

Faint Dwarfs in Nearby Groups

Authors: Ryan Speller, James E. Taylor

Abstract: The number and distribution of dwarf satellite galaxies remain a critical test of cold dark matter-dominated structure formation on small scales. Until recently, observational information about galaxy formation on these scales has been limited mainly to the Local Group. We have searched for faint analogues of Local Group dwarfs around nearby bright galaxies, using a spatial clustering analysis of… ▽ More The number and distribution of dwarf satellite galaxies remain a critical test of cold dark matter-dominated structure formation on small scales. Until recently, observational information about galaxy formation on these scales has been limited mainly to the Local Group. We have searched for faint analogues of Local Group dwarfs around nearby bright galaxies, using a spatial clustering analysis of the photometric catalog of the Sloan Digital Sky Survey (SDSS) Data Release 8. Several other recent searches of SDSS have detected clustered satellite populations down to $Δm_r \equiv ({m}_{r,\, {\rm sat}} -\, {m}_{r,\, {\rm main}}) \sim 6$-$8$, using photometric redshifts to reduce background contamination. SDSS photometric redshifts are relatively imprecise, however, for faint and nearby galaxies. Instead we use angular size to select potential nearby dwarfs, and consider only the nearest isolated bright galaxies as primaries. As a result, we are able to detect an excess clustering signal from companions down to $Δm_r = 12$, four magnitudes fainter than most recent studies. We detect an over-density of objects at separations $< 400$ kpc, corresponding to about $4.6 \pm 0.5$ satellites per central galaxy, consistent with the satellite abundance expected from the Local Group given our selection function. Although the sample of satellites detected is incomplete by construction, since it excludes the least and most compact dwarfs, this detection provides a lower bound on the average satellite luminosity function, down to luminosities corresponding to the faintest "classical" dwarfs of the Local Group. △ Less

Submitted 29 October, 2013; originally announced October 2013.

Comments: 12 pages, 7 figures, submitted to ApJ