Showing 1–50 of 146 results for author: Brunner, R

A Triangulation of the Flow Polytope of the Zigzag Graph

Authors: Rachel Brunner, Christopher R. H. Hanusa

Abstract: We show that the dual graph of the triangulation of the flow polytope of the zigzag graph adorned with the length-reverse-length framing is a subgraph of a grid graph. Through Mészáros, Morales, and Striker's bijection between simplices of the triangulation, integer flows of a different, supplemental flow polytope, we provide a simple numerical characterization of the adjacency between the triangu… ▽ More We show that the dual graph of the triangulation of the flow polytope of the zigzag graph adorned with the length-reverse-length framing is a subgraph of a grid graph. Through Mészáros, Morales, and Striker's bijection between simplices of the triangulation, integer flows of a different, supplemental flow polytope, we provide a simple numerical characterization of the adjacency between the triangulation's simplices in terms of their corresponding integer flows. The proofs result from the development of Postnikov and Stanley's sequences of noncrossing bipartite trees as combinatorial objects we call groves. We propose two new statistics derived from this construction that we conjecture recover the $h^*$-polynomial of the flow polytope of the zigzag graph. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: 20 pages, 12 figures, 1 table

MSC Class: Primary: 05A19; 05C21; 52B12; 52B05; 52C07; Secondary: 05C05; 05C20; 05C30; 52A38; 52B11

Forensic Data Analytics for Anomaly Detection in Evolving Networks

Authors: Li Yang, Abdallah Moubayed, Abdallah Shami, Amine Boukhtouta, Parisa Heidari, Stere Preda, Richard Brunner, Daniel Migault, Adel Larabi

Abstract: In the prevailing convergence of traditional infrastructure-based deployment (i.e., Telco and industry operational networks) towards evolving deployments enabled by 5G and virtualization, there is a keen interest in elaborating effective security controls to protect these deployments in-depth. By considering key enabling technologies like 5G and virtualization, evolving networks are democratized,… ▽ More In the prevailing convergence of traditional infrastructure-based deployment (i.e., Telco and industry operational networks) towards evolving deployments enabled by 5G and virtualization, there is a keen interest in elaborating effective security controls to protect these deployments in-depth. By considering key enabling technologies like 5G and virtualization, evolving networks are democratized, facilitating the establishment of point presences integrating different business models ranging from media, dynamic web content, gaming, and a plethora of IoT use cases. Despite the increasing services provided by evolving networks, many cybercrimes and attacks have been launched in evolving networks to perform malicious activities. Due to the limitations of traditional security artifacts (e.g., firewalls and intrusion detection systems), the research on digital forensic data analytics has attracted more attention. Digital forensic analytics enables people to derive detailed information and comprehensive conclusions from different perspectives of cybercrimes to assist in convicting criminals and preventing future crimes. This chapter presents a digital analytics framework for network anomaly detection, including multi-perspective feature engineering, unsupervised anomaly detection, and comprehensive result correction procedures. Experiments on real-world evolving network data show the effectiveness of the proposed forensic data analytics solution. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: Electronic version of an article published as [Book Series: World Scientific Series in Digital Forensics and Cybersecurity, Volume 2, Innovations in Digital Forensics, 2023, Pages 99-137] [DOI:10.1142/9789811273209_0004] \c{opyright} copyright World Scientific Publishing Company [https://doi.org/10.1142/9789811273209_0004]

MSC Class: 68T01 ACM Class: I.2.6; C.2.0

Multi-Perspective Content Delivery Networks Security Framework Using Optimized Unsupervised Anomaly Detection

Authors: Li Yang, Abdallah Moubayed, Abdallah Shami, Parisa Heidari, Amine Boukhtouta, Adel Larabi, Richard Brunner, Stere Preda, Daniel Migault

Abstract: Content delivery networks (CDNs) provide efficient content distribution over the Internet. CDNs improve the connectivity and efficiency of global communications, but their caching mechanisms may be breached by cyber-attackers. Among the security mechanisms, effective anomaly detection forms an important part of CDN security enhancement. In this work, we propose a multi-perspective unsupervised lea… ▽ More Content delivery networks (CDNs) provide efficient content distribution over the Internet. CDNs improve the connectivity and efficiency of global communications, but their caching mechanisms may be breached by cyber-attackers. Among the security mechanisms, effective anomaly detection forms an important part of CDN security enhancement. In this work, we propose a multi-perspective unsupervised learning framework for anomaly detection in CDNs. In the proposed framework, a multi-perspective feature engineering approach, an optimized unsupervised anomaly detection model that utilizes an isolation forest and a Gaussian mixture model, and a multi-perspective validation method, are developed to detect abnormal behaviors in CDNs mainly from the client Internet Protocol (IP) and node perspectives, therefore to identify the denial of service (DoS) and cache pollution attack (CPA) patterns. Experimental results are presented based on the analytics of eight days of real-world CDN log data provided by a major CDN operator. Through experiments, the abnormal contents, compromised nodes, malicious IPs, as well as their corresponding attack types, are identified effectively by the proposed framework and validated by multiple cybersecurity experts. This shows the effectiveness of the proposed method when applied to real-world CDN data. △ Less

Submitted 23 July, 2021; originally announced July 2021.

Comments: Accepted and to Appear in IEEE Transactions on Network and Service Management

MSC Class: 68T01 ACM Class: I.2.6; C.2.0

A Spitzer survey of Deep Drilling Fields to be targeted by the Vera C. Rubin Observatory Legacy Survey of Space and Time

Authors: M. Lacy, J. A. Surace, D. Farrah, K. Nyland, J. Afonso, W. N. Brandt, D. L. Clements, C. D. P. Lagos, C. Maraston, J. Pforr, A. Sajina, M. Sako, M. Vaccari, G. Wilson, D. R. Ballantyne, W. A. Barkhouse, R. Brunner, R. Cane, T. E. Clarke, M. Cooper, A. Cooray, G. Covone, C. D'Andrea, A. E. Evrard, H. C. Ferguson , et al. (38 additional authors not shown)

Abstract: The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will observe several Deep Drilling Fields (DDFs) to a greater depth and with a more rapid cadence than the main survey. In this paper, we describe the ``DeepDrill'' survey, which used the Spitzer Space Telescope Infrared Array Camera (IRAC) to observe three of the four currently defined DDFs in two bands, centered on 3.6 $μ$m and… ▽ More The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will observe several Deep Drilling Fields (DDFs) to a greater depth and with a more rapid cadence than the main survey. In this paper, we describe the ``DeepDrill'' survey, which used the Spitzer Space Telescope Infrared Array Camera (IRAC) to observe three of the four currently defined DDFs in two bands, centered on 3.6 $μ$m and 4.5 $μ$m. These observations expand the area which was covered by an earlier set of observations in these three fields by the Spitzer Extragalactic Representative Volume Survey (SERVS). The combined DeepDrill and SERVS data cover the footprints of the LSST DDFs in the Extended Chandra Deep Field-South field (ECDFS), the ELAIS-S1 field (ES1), and the XMM Large-Scale Structure Survey field (XMM-LSS). The observations reach an approximate $5σ$ point-source depth of 2 $μ$Jy (corresponding to an AB magnitude of 23.1; sufficient to detect a 10$^{11} M_{\odot}$ galaxy out to $z\approx 5$) in each of the two bands over a total area of $\approx 29\,$deg$^2$. The dual-band catalogues contain a total of 2.35 million sources. In this paper we describe the observations and data products from the survey, and an overview of the properties of galaxies in the survey. We compare the source counts to predictions from the SHARK semi-analytic model of galaxy formation. We also identify a population of sources with extremely red ([3.6]$-$[4.5] $>1.2$) colours which we show mostly consists of highly-obscured active galactic nuclei. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Comments: 20 pages, 9 figures; MNRAS in press

Cost-optimal V2X Service Placement in Distributed Cloud/Edge Environment

Authors: Abdallah Moubayed, Abdallah Shami, Parisa Heidari, Adel Larabi, Richard Brunner

Abstract: Deploying V2X services has become a challenging task. This is mainly due to the fact that such services have strict latency requirements. To meet these requirements, one potential solution is adopting mobile edge computing (MEC). However, this presents new challenges including how to find a cost efficient placement that meets other requirements such as latency. In this work, the problem of cost-op… ▽ More Deploying V2X services has become a challenging task. This is mainly due to the fact that such services have strict latency requirements. To meet these requirements, one potential solution is adopting mobile edge computing (MEC). However, this presents new challenges including how to find a cost efficient placement that meets other requirements such as latency. In this work, the problem of cost-optimal V2X service placement (CO-VSP) in a distributed cloud/edge environment is formulated. Additionally, a cost-focused delay-aware V2X service placement (DA-VSP) heuristic algorithm is proposed. Simulation results show that both CO-VSP model and DA-VSP algorithm guarantee the QoS requirements of all such services and illustrates the trade-off between latency and deployment cost. △ Less

Submitted 14 October, 2020; originally announced October 2020.

Comments: 6 pages, 4 figures, 1 table, 1 algorithm pseudocode, Accepted & presented in IEEE 16th International Conference on Wireless and Mobile Computing, Networking and Communications (WiMob 2020)

The 1st Agriculture-Vision Challenge: Methods and Results

Authors: Mang Tik Chiu, Xingqian Xu, Kai Wang, Jennifer Hobbs, Naira Hovakimyan, Thomas S. Huang, Honghui Shi, Yunchao Wei, Zilong Huang, Alexander Schwing, Robert Brunner, Ivan Dozier, Wyatt Dozier, Karen Ghandilyan, David Wilson, Hyunseong Park, Junhee Kim, Sungho Kim, Qinghui Liu, Michael C. Kampffmeyer, Robert Jenssen, Arnt B. Salberg, Alexandre Barbosa, Rodrigo Trevisan, Bingchen Zhao , et al. (17 additional authors not shown)

Abstract: The first Agriculture-Vision Challenge aims to encourage research in developing novel and effective algorithms for agricultural pattern recognition from aerial images, especially for the semantic segmentation task associated with our challenge dataset. Around 57 participating teams from various countries compete to achieve state-of-the-art in aerial agriculture semantic segmentation. The Agricultu… ▽ More The first Agriculture-Vision Challenge aims to encourage research in developing novel and effective algorithms for agricultural pattern recognition from aerial images, especially for the semantic segmentation task associated with our challenge dataset. Around 57 participating teams from various countries compete to achieve state-of-the-art in aerial agriculture semantic segmentation. The Agriculture-Vision Challenge Dataset was employed, which comprises of 21,061 aerial and multi-spectral farmland images. This paper provides a summary of notable methods and results in the challenge. Our submission server and leaderboard will continue to open for researchers that are interested in this challenge dataset and task; the link can be found here. △ Less

Submitted 23 April, 2020; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: CVPR 2020 Workshop

Machine Learning for Performance-Aware Virtual Network Function Placement

Authors: Dimitrios Michael Manias, Manar Jammal, Hassan Hawilo, Abdallah Shami, Parisa Heidari, Adel Larabi, Richard Brunner

Abstract: With the growing demand for data connectivity, network service providers are faced with the task of reducing their capital and operational expenses while simultaneously improving network performance and addressing the increased connectivity demand. Although Network Function Virtualization (NFV) has been identified as a solution, several challenges must be addressed to ensure its feasibility. In th… ▽ More With the growing demand for data connectivity, network service providers are faced with the task of reducing their capital and operational expenses while simultaneously improving network performance and addressing the increased connectivity demand. Although Network Function Virtualization (NFV) has been identified as a solution, several challenges must be addressed to ensure its feasibility. In this paper, we address the Virtual Network Function (VNF) placement problem by developing a machine learning decision tree model that learns from the effective placement of the various VNF instances forming a Service Function Chain (SFC). The model takes several performance-related features from the network as an input and selects the placement of the various VNF instances on network servers with the objective of minimizing the delay between dependent VNF instances. The benefits of using machine learning are realized by moving away from a complex mathematical modelling of the system and towards a data-based understanding of the system. Using the Evolved Packet Core (EPC) as a use case, we evaluate our model on different data center networks and compare it to the BACON algorithm in terms of the delay between interconnected components and the total delay across the SFC. Furthermore, a time complexity analysis is performed to show the effectiveness of the model in NFV applications. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Comments: 6 pages, 6 figures, 1 table, 9 equations, 18 references, Conference

Edge-enabled V2X Service Placement for Intelligent Transportation Systems

Authors: Abdallah Moubayed, Abdallah Shami, Parisa Heidari, Adel Larabi, Richard Brunner

Abstract: Vehicle-to-everything (V2X) communication and services have been garnering significant interest from different stakeholders as part of future intelligent transportation systems (ITSs). This is due to the many benefits they offer. However, many of these services have stringent performance requirements, particularly in terms of the delay/latency. Multi-access/mobile edge computing (MEC) has been pro… ▽ More Vehicle-to-everything (V2X) communication and services have been garnering significant interest from different stakeholders as part of future intelligent transportation systems (ITSs). This is due to the many benefits they offer. However, many of these services have stringent performance requirements, particularly in terms of the delay/latency. Multi-access/mobile edge computing (MEC) has been proposed as a potential solution for such services by bringing them closer to vehicles. Yet, this introduces a new set of challenges such as where to place these V2X services, especially given the limit computation resources available at edge nodes. To that end, this work formulates the problem of optimal V2X service placement (OVSP) in a hybrid core/edge environment as a binary integer linear programming problem. To the best of our knowledge, no previous work considered the V2X service placement problem while taking into consideration the computational resource availability at the nodes. Moreover, a low-complexity greedy-based heuristic algorithm named "Greedy V2X Service Placement Algorithm" (G-VSPA) was developed to solve this problem. Simulation results show that the OVSP model successfully guarantees and maintains the QoS requirements of all the different V2X services. Additionally, it is observed that the proposed G-VSPA algorithm achieves close to optimal performance while having lower complexity. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Comments: 13 pages, 16 figures (including 5 bio pictures), accepted and to be published in IEEE Transactions on Mobile Computing

Agriculture-Vision: A Large Aerial Image Database for Agricultural Pattern Analysis

Authors: Mang Tik Chiu, Xingqian Xu, Yunchao Wei, Zilong Huang, Alexander Schwing, Robert Brunner, Hrant Khachatrian, Hovnatan Karapetyan, Ivan Dozier, Greg Rose, David Wilson, Adrian Tudor, Naira Hovakimyan, Thomas S. Huang, Honghui Shi

Abstract: The success of deep learning in visual recognition tasks has driven advancements in multiple fields of research. Particularly, increasing attention has been drawn towards its application in agriculture. Nevertheless, while visual pattern recognition on farmlands carries enormous economic values, little progress has been made to merge computer vision and crop sciences due to the lack of suitable ag… ▽ More The success of deep learning in visual recognition tasks has driven advancements in multiple fields of research. Particularly, increasing attention has been drawn towards its application in agriculture. Nevertheless, while visual pattern recognition on farmlands carries enormous economic values, little progress has been made to merge computer vision and crop sciences due to the lack of suitable agricultural image datasets. Meanwhile, problems in agriculture also pose new challenges in computer vision. For example, semantic segmentation of aerial farmland images requires inference over extremely large-size images with extreme annotation sparsity. These challenges are not present in most of the common object datasets, and we show that they are more challenging than many other aerial image datasets. To encourage research in computer vision for agriculture, we present Agriculture-Vision: a large-scale aerial farmland image dataset for semantic segmentation of agricultural patterns. We collected 94,986 high-quality aerial images from 3,432 farmlands across the US, where each image consists of RGB and Near-infrared (NIR) channels with resolution as high as 10 cm per pixel. We annotate nine types of field anomaly patterns that are most important to farmers. As a pilot study of aerial agricultural semantic segmentation, we perform comprehensive experiments using popular semantic segmentation models; we also propose an effective model designed for aerial agricultural pattern recognition. Our experiments demonstrate several challenges Agriculture-Vision poses to both the computer vision and agriculture communities. Future versions of this dataset will include even more aerial images, anomaly patterns and image channels. More information at https://www.agriculture-vision.com. △ Less

Submitted 19 March, 2020; v1 submitted 5 January, 2020; originally announced January 2020.

Comments: CVPR 2020

Unsupervised Star Galaxy Classification with Cascade Variational Auto-Encoder

Authors: Hao Sun, Jiadong Guo, Edward J. Kim, Robert J. Brunner

Abstract: The increasing amount of data in astronomy provides great challenges for machine learning research. Previously, supervised learning methods achieved satisfactory recognition accuracy for the star-galaxy classification task, based on manually labeled data set. In this work, we propose a novel unsupervised approach for the star-galaxy recognition task, namely Cascade Variational Auto-Encoder (CasVAE… ▽ More The increasing amount of data in astronomy provides great challenges for machine learning research. Previously, supervised learning methods achieved satisfactory recognition accuracy for the star-galaxy classification task, based on manually labeled data set. In this work, we propose a novel unsupervised approach for the star-galaxy recognition task, namely Cascade Variational Auto-Encoder (CasVAE). Our empirical results show our method outperforms the baseline model in both accuracy and stability. △ Less

Submitted 30 October, 2019; originally announced October 2019.

Extended Isolation Forest

Authors: Sahand Hariri, Matias Carrasco Kind, Robert J. Brunner

Abstract: We present an extension to the model-free anomaly detection algorithm, Isolation Forest. This extension, named Extended Isolation Forest (EIF), resolves issues with assignment of anomaly score to given data points. We motivate the problem using heat maps for anomaly scores. These maps suffer from artifacts generated by the criteria for branching operation of the binary tree. We explain this proble… ▽ More We present an extension to the model-free anomaly detection algorithm, Isolation Forest. This extension, named Extended Isolation Forest (EIF), resolves issues with assignment of anomaly score to given data points. We motivate the problem using heat maps for anomaly scores. These maps suffer from artifacts generated by the criteria for branching operation of the binary tree. We explain this problem in detail and demonstrate the mechanism by which it occurs visually. We then propose two different approaches for improving the situation. First we propose transforming the data randomly before creation of each tree, which results in averaging out the bias. Second, which is the preferred way, is to allow the slicing of the data to use hyperplanes with random slopes. This approach results in remedying the artifact seen in the anomaly score heat maps. We show that the robustness of the algorithm is much improved using this method by looking at the variance of scores of data points distributed along constant level sets. We report AUROC and AUPRC for our synthetic datasets, along with real-world benchmark datasets. We find no appreciable difference in the rate of convergence nor in computation time between the standard Isolation Forest and EIF. △ Less

Submitted 8 July, 2020; v1 submitted 5 November, 2018; originally announced November 2018.

Comments: 12 pages; 21 figures, Published. Open source code in https://github.com/sahandha/eif

Structure and Content of the Visible Darknet

Authors: Georgia Avarikioti, Roman Brunner, Aggelos Kiayias, Roger Wattenhofer, Dionysis Zindros

Abstract: In this paper, we analyze the topology and the content found on the "darknet", the set of websites accessible via Tor. We created a darknet spider and crawled the darknet starting from a bootstrap list by recursively following links. We explored the whole connected component of more than 34,000 hidden services, of which we found 10,000 to be online. Contrary to folklore belief, the visible part of… ▽ More In this paper, we analyze the topology and the content found on the "darknet", the set of websites accessible via Tor. We created a darknet spider and crawled the darknet starting from a bootstrap list by recursively following links. We explored the whole connected component of more than 34,000 hidden services, of which we found 10,000 to be online. Contrary to folklore belief, the visible part of the darknet is surprisingly well-connected through hub websites such as wikis and forums. We performed a comprehensive categorization of the content using supervised machine learning. We observe that about half of the visible dark web content is related to apparently licit activities based on our classifier. A significant amount of content pertains to software repositories, blogs, and activism-related websites. Among unlawful hidden services, most pertain to fraudulent websites, services selling counterfeit goods, and drug markets. △ Less

Submitted 7 November, 2018; v1 submitted 4 November, 2018; originally announced November 2018.

Star-galaxy classification in the Dark Energy Survey Y1 dataset

Authors: I. Sevilla-Noarbe, B. Hoyle, M. J. Marchã, M. T. Soumagnac, K. Bechtol, A. Drlica-Wagner, F. Abdalla, J. Aleksić, C. Avestruz, E. Balbinot, M. Banerji, E. Bertin, C. Bonnett, R. Brunner, M. Carrasco-Kind, A. Choi, T. Giannantonio, E. Kim, O. Lahav, B. Moraes, B. Nord, A. J. Ross, E. S. Rykoff, B. Santiago, E. Sheldon , et al. (53 additional authors not shown)

Abstract: We perform a comparison of different approaches to star-galaxy classification using the broad-band photometric data from Year 1 of the Dark Energy Survey. This is done by performing a wide range of tests with and without external `truth' information, which can be ported to other similar datasets. We make a broad evaluation of the performance of the classifiers in two science cases with DES data th… ▽ More We perform a comparison of different approaches to star-galaxy classification using the broad-band photometric data from Year 1 of the Dark Energy Survey. This is done by performing a wide range of tests with and without external `truth' information, which can be ported to other similar datasets. We make a broad evaluation of the performance of the classifiers in two science cases with DES data that are most affected by this systematic effect: large-scale structure and Milky Way studies. In general, even though the default morphological classifiers used for DES Y1 cosmology studies are sufficient to maintain a low level of systematic contamination from stellar mis-classification, contamination can be reduced to the O(1%) level by using multi-epoch and infrared information from external datasets. For Milky Way studies the stellar sample can be augmented by ~20% for a given flux limit. Reference catalogs used in this work will be made available upon publication. △ Less

Submitted 30 October, 2018; v1 submitted 7 May, 2018; originally announced May 2018.

Comments: Reference catalogs used in this work will be made available upon publication

Report number: Fermilab-PUB-18-112-AE-PPD

Journal ref: Monthly Notices of the Royal Astronomical Society, 2018, Volume 481, Issue 4, p.5451-5469

An NFV and Microservice Based Architecture for On-the-fly Component Provisioning in Content Delivery Networks

Authors: Narjes Tahghigh Jahromi, Roch H. Glitho, Adel Larabi, Richard Brunner

Abstract: Content Delivery Networks (CDNs) deliver content (e.g. Web pages, videos) to geographically distributed end-users over the Internet. Some contents do sometimes attract the attention of a large group of end-users. This often leads to flash crowds which can cause major issues such as outage in the CDN. Microservice architectural style aims at decomposing monolithic systems into smaller components wh… ▽ More Content Delivery Networks (CDNs) deliver content (e.g. Web pages, videos) to geographically distributed end-users over the Internet. Some contents do sometimes attract the attention of a large group of end-users. This often leads to flash crowds which can cause major issues such as outage in the CDN. Microservice architectural style aims at decomposing monolithic systems into smaller components which can be independently deployed, upgraded and disposed. Network Function Virtualization (NFV) is an emerging technology that aims to reduce costs and bring agility by decoupling network functions from the underlying hardware. This paper leverages the NFV and microservice architectural style to propose an architecture for on-the-fly CDN component provisioning to tackle issues such as flash crowds. In the proposed architecture, CDN components are designed as sets of microservices which interact via RESTFul Web services and are provisioned as Virtual Network Functions (VNFs), which are deployed and orchestrated on-the-fly. We have built a prototype in which a CDN surrogate server, designed as a set of microservices, is deployed on-the-fly. The prototype is deployed on SAVI, a Canadian distributed test bed for future Internet applications. The performance is also evaluated. △ Less

Submitted 13 October, 2017; originally announced October 2017.

Comments: 7 pages, 6 figures, 2 tables

Workload Analysis of Blue Waters

Authors: Matthew D. Jones, Joseph P. White, Martins Innus, Robert L. DeLeon, Nikolay Simakov, Jeffrey T. Palmer, Steven M. Gallo, Thomas R. Furlani, Michael Showerman, Robert Brunner, Andry Kot, Gregory Bauer, Brett Bode, Jeremy Enos, William Kramer

Abstract: Blue Waters is a Petascale-level supercomputer whose mission is to enable the national scientific and research community to solve "grand challenge" problems that are orders of magnitude more complex than can be carried out on other high performance computing systems. Given the important and unique role that Blue Waters plays in the U.S. research portfolio, it is important to have a detailed unders… ▽ More Blue Waters is a Petascale-level supercomputer whose mission is to enable the national scientific and research community to solve "grand challenge" problems that are orders of magnitude more complex than can be carried out on other high performance computing systems. Given the important and unique role that Blue Waters plays in the U.S. research portfolio, it is important to have a detailed understanding of its workload in order to guide performance optimization both at the software and system configuration level as well as inform architectural balance tradeoffs. Furthermore, understanding the computing requirements of the Blue Water's workload (memory access, IO, communication, etc.), which is comprised of some of the most computationally demanding scientific problems, will help drive changes in future computing architectures, especially at the leading edge. With this objective in mind, the project team carried out a detailed workload analysis of Blue Waters. △ Less

Submitted 2 March, 2017; originally announced March 2017.

Comments: 107 pages, >100 figures (figure sizes reduced to save space, contact authors for version with full resolution)

MSC Class: 68M14; 68M20; 68U20 ACM Class: I.6.3; J.2; J.3; J.4; J.5; K.6.4

Vizic: A Jupyter-based Interactive Visualization Tool for Astronomical Catalogs

Authors: W. Yu, M. Carrasco Kind, R. J. Brunner

Abstract: The ever-growing datasets in observational astronomy have challenged scientists in many aspects, including an efficient and interactive data exploration and visualization. Many tools have been developed to confront this challenge. However, they usually focus on displaying the actual images or focus on visualizing patterns within catalogs in a predefined way. In this paper we introduce Vizic, a Pyt… ▽ More The ever-growing datasets in observational astronomy have challenged scientists in many aspects, including an efficient and interactive data exploration and visualization. Many tools have been developed to confront this challenge. However, they usually focus on displaying the actual images or focus on visualizing patterns within catalogs in a predefined way. In this paper we introduce Vizic, a Python visualization library that builds the connection between images and catalogs through an interactive map of the sky region. Vizic visualizes catalog data over a custom background canvas using the shape, size and orientation of each object in the catalog. The displayed objects in the map are highly interactive and customizable comparing to those in the images. These objects can be filtered by or colored by their properties, such as redshift and magnitude. They also can be sub-selected using a lasso-like tool for further analysis using standard Python functions from inside a Jupyter notebook. Furthermore, Vizic allows custom overlays to be appended dynamically on top of the sky map. We have initially implemented several overlays, namely, Voronoi, Delaunay, Minimum Spanning Tree and HEALPix grid layers, which are helpful for visualizing large-scale structure. All these overlays can be generated, added or removed interactively with one line of code. The catalog data is stored in a non-relational database, and the interfaces were developed in JavaScript and Python to work within Jupyter Notebook, which allows to create custom widgets, user generated scripts to analyze and plot the data selected/displayed in the interactive map. This unique design makes Vizic a very powerful and flexible interactive analysis tool. Vizic can be adopted in variety of exercises, for example, data inspection, clustering analysis, galaxy alignment studies, outlier identification or simply large-scale visualizations. △ Less

Submitted 6 May, 2017; v1 submitted 5 January, 2017; originally announced January 2017.

Comments: 14 pages, 13 figures, revised for Astronomy and Computing

Journal ref: Astron. Comput. 20 (2017) 128-139

Star-galaxy Classification Using Deep Convolutional Neural Networks

Authors: Edward J. Kim, Robert J. Brunner

Abstract: Most existing star-galaxy classifiers use the reduced summary information from catalogs, requiring careful feature extraction and selection. The latest advances in machine learning that use deep convolutional neural networks allow a machine to automatically learn the features directly from data, minimizing the need for input from human experts. We present a star-galaxy classification framework tha… ▽ More Most existing star-galaxy classifiers use the reduced summary information from catalogs, requiring careful feature extraction and selection. The latest advances in machine learning that use deep convolutional neural networks allow a machine to automatically learn the features directly from data, minimizing the need for input from human experts. We present a star-galaxy classification framework that uses deep convolutional neural networks (ConvNets) directly on the reduced, calibrated pixel values. Using data from the Sloan Digital Sky Survey (SDSS) and the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we demonstrate that ConvNets are able to produce accurate and well-calibrated probabilistic classifications that are competitive with conventional machine learning techniques. Future advances in deep learning may bring more success with current and forthcoming photometric surveys, such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope (LSST), because deep neural networks require very little, manual feature engineering. △ Less

Submitted 13 October, 2016; v1 submitted 15 August, 2016; originally announced August 2016.

Comments: 13 page, 13 figures. Accepted for publication in the MNRAS. Code available at https://github.com/EdwardJKim/dl4astro

Teaching Data Science

Authors: Robert J. Brunner, Edward J. Kim

Abstract: We describe an introductory data science course, entitled Introduction to Data Science, offered at the University of Illinois at Urbana-Champaign. The course introduced general programming concepts by using the Python programming language with an emphasis on data preparation, processing, and presentation. The course had no prerequisites, and students were not expected to have any programming exper… ▽ More We describe an introductory data science course, entitled Introduction to Data Science, offered at the University of Illinois at Urbana-Champaign. The course introduced general programming concepts by using the Python programming language with an emphasis on data preparation, processing, and presentation. The course had no prerequisites, and students were not expected to have any programming experience. This introductory course was designed to cover a wide range of topics, from the nature of data, to storage, to visualization, to probability and statistical analysis, to cloud and high performance computing, without becoming overly focused on any one subject. We conclude this article with a discussion of lessons learned and our plans to develop new data science courses. △ Less

Submitted 25 April, 2016; originally announced April 2016.

Comments: 10 pages, 4 figures, International Conference on Computational Science (ICCS 2016)

Galaxy clustering with photometric surveys using PDF redshift information

Authors: J. Asorey, M. Carrasco Kind, I. Sevilla-Noarbe, R. J. Brunner, J. Thaler

Abstract: Photometric surveys produce large-area maps of the galaxy distribution, but with less accurate redshift information than is obtained from spectroscopic methods. Modern photometric redshift (photo-z) algorithms use galaxy magnitudes, or colors, that are obtained through multi-band imaging to produce a probability density function (PDF) for each galaxy in the map. We used simulated data to study the… ▽ More Photometric surveys produce large-area maps of the galaxy distribution, but with less accurate redshift information than is obtained from spectroscopic methods. Modern photometric redshift (photo-z) algorithms use galaxy magnitudes, or colors, that are obtained through multi-band imaging to produce a probability density function (PDF) for each galaxy in the map. We used simulated data to study the effect of using different photo-z estimators to assign galaxies to redshift bins in order to compare their effects on angular clustering and galaxy bias measurements. We found that if we use the entire PDF, rather than a single-point (mean or mode) estimate, the deviations are less biased, especially when using narrow redshift bins. When the redshift bin widths are $Δz=0.1$, the use of the entire PDF reduces the typical measurement bias from 5%, when using single point estimates, to 3%. △ Less

Submitted 13 June, 2016; v1 submitted 3 January, 2016; originally announced January 2016.

Comments: Matches the MNRAS published version. 19 pages, 19 Figures

The Dark Energy Survey: more than dark energy - an overview

Authors: Dark Energy Survey Collaboration, T. Abbott, F. B. Abdalla, J. Aleksic, S. Allam, A. Amara, D. Bacon, E. Balbinot, M. Banerji, K. Bechtol, A. Benoit-Levy, G. M. Bernstein, E. Bertin, J. Blazek, C. Bonnett, S. Bridle, D. Brooks, R. J. Brunner, E. Buckley-Geer, D. L. Burke, G. B. Caminha, D. Capozzi, J. Carlsen, A. Carnero-Rosell, M. Carollo , et al. (116 additional authors not shown)

Abstract: This overview article describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is using a wide-field camera (DECam) on the 4m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion the survey is expected to have generated a catalogue of… ▽ More This overview article describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is using a wide-field camera (DECam) on the 4m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion the survey is expected to have generated a catalogue of 300 million galaxies with photometric redshifts and 100 million stars. In addition, a time-domain survey search over 27 sq deg is expected to yield a sample of thousands of Type Ia supernovae and other transients. The main goals of DES are to characterise dark energy and dark matter, and to test alternative models of gravity; these goals will be pursued by studying large scale structure, cluster counts, weak gravitational lensing and Type Ia supernovae. However, DES also provides a rich data set which allows us to study many other aspects of astrophysics. In this paper we focus on additional science with DES, emphasizing areas where the survey makes a difference with respect to other current surveys. The paper illustrates, using early data (from `Science Verification', and from the first, second and third seasons of observations), what DES can tell us about the solar system, the Milky Way, galaxy evolution, quasars, and other topics. In addition, we show that if the cosmological model is assumed to be Lambda+ Cold Dark Matter (LCDM) then important astrophysics can be deduced from the primary DES probes. Highlights from DES early data include the discovery of 34 Trans Neptunian Objects, 17 dwarf satellites of the Milky Way, one published z > 6 quasar (and more confirmed) and two published superluminous supernovae (and more confirmed). △ Less

Submitted 19 August, 2016; v1 submitted 3 January, 2016; originally announced January 2016.

Comments: 32 pages, 15 figures; a revised Figure 1 and minor changes, to match the published MNRAS version

Report number: DES 2015-0085; FERMILAB-PUB-16-003-AE

Journal ref: Mon.Not.Roy.Astron.Soc. 460 (2016) no. 2, 1270-1299

Observation and Confirmation of Six Strong Lensing Systems in The Dark Energy Survey Science Verification Data

Authors: B. Nord, E. Buckley-Geer, H. Lin, H. T. Diehl, J. Helsby, N. Kuropatkin, A. Amara, T. Collett, S. Allam, G. Caminha, C. De Bom, S. Desai, H. Dúmet-Montoya, M. Elidaiana da S. Pereira, D. A. Finley, B. Flaugher, C. Furlanetto, H. Gaitsch, M. Gill, K. W. Merritt, A. More, D. Tucker, E. S. Rykoff, E. Rozo, F. B. Abdalla , et al. (67 additional authors not shown)

Abstract: We report the observation and confirmation of the first group- and cluster-scale strong gravitational lensing systems found in Dark Energy Survey (DES) data. Through visual inspection of data from the Science Verification (SV) season, we identified 53 candidate systems. We then obtained spectroscopic follow-up of 21 candidates using the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South t… ▽ More We report the observation and confirmation of the first group- and cluster-scale strong gravitational lensing systems found in Dark Energy Survey (DES) data. Through visual inspection of data from the Science Verification (SV) season, we identified 53 candidate systems. We then obtained spectroscopic follow-up of 21 candidates using the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South telescope and the Inamori-Magellan Areal Camera and Spectrograph (IMACS) at the Magellan/Baade telescope. With this follow-up, we confirmed six candidates as gravitational lenses: Three of the systems are newly discovered, and the remaining three were previously known. Of the 21 observed candidates, the remaining 15 were either not detected in spectroscopic observations, were observed and did not exhibit continuum emission (or spectral features), or were ruled out as lensing systems. The confirmed sample consists of one group-scale and five galaxy cluster-scale lenses. The lensed sources range in redshift z ~ 0.80-3.2, and in i-band surface brightness i_{SB} ~ 23-25 mag/sq.-arcsec. (2" aperture). For each of the six systems, we estimate the Einstein radius and the enclosed mass, which have ranges ~ 5.0 - 8.6" and ~ 7.5 x 10^{12} - 6.4 x 10^{13} solar masses, respectively. △ Less

Submitted 9 December, 2015; originally announced December 2015.

Comments: 17 pages, 7 figures, 4 tables; submitted to ApJ

Report number: Fermilab-PUB-15-509-AE

Creating updated, scientifically-calibrated mosaic images for the RC3 catalogue

Authors: Jung Lin Lee, Robert J. Brunner

Abstract: The Third Reference Catalogue of Bright Galaxies (RC3) is a reasonably complete listing of 23,011 nearby, large, bright galaxies. By using the final imaging data release from the Sloan Digital Sky Survey, we generate scientifically-calibrated FITS mosaics by using the montage program for all SDSS imaging bands for all RC3 galaxies that lie within the survey footprint. We further combine the SDSS g… ▽ More The Third Reference Catalogue of Bright Galaxies (RC3) is a reasonably complete listing of 23,011 nearby, large, bright galaxies. By using the final imaging data release from the Sloan Digital Sky Survey, we generate scientifically-calibrated FITS mosaics by using the montage program for all SDSS imaging bands for all RC3 galaxies that lie within the survey footprint. We further combine the SDSS g, r, and i band FITS mosaics for these galaxies to create color-composite images by using the STIFF program. We generalized this software framework to make FITS mosaics and color-composite images for an arbitrary catalog and imaging data set. Due to positional inaccuracies inherent in the RC3 catalog, we employ a recursive algorithm in our mosaicking pipeline that first determines the correct location for each galaxy, and subsequently applies the mosaicking procedure. As an additional test of this new software pipeline and to obtain mosaic images of a larger sample of RC3 galaxies, we also applied this pipeline to photographic data taken by the Second Palomar Observatory Sky Survey with $B_J$, $R_F$, and $I_N$ plates. We publicly release all generated data, accessible via a web search form, and the software pipeline to enable others to make galaxy mosaics by using other catalogs or surveys. △ Less

Submitted 3 December, 2015; originally announced December 2015.

Comments: 11 pages, 13 figures

Machine Learning and Cosmological Simulations II: Hydrodynamical Simulations

Authors: Harshil M. Kamdar, Matthew J. Turk, Robert J. Brunner

Abstract: We extend a machine learning (ML) framework presented previously to model galaxy formation and evolution in a hierarchical universe using N-body + hydrodynamical simulations. In this work, we show that ML is a promising technique to study galaxy formation in the backdrop of a hydrodynamical simulation. We use the Illustris Simulation to train and test various sophisticated machine learning algorit… ▽ More We extend a machine learning (ML) framework presented previously to model galaxy formation and evolution in a hierarchical universe using N-body + hydrodynamical simulations. In this work, we show that ML is a promising technique to study galaxy formation in the backdrop of a hydrodynamical simulation. We use the Illustris Simulation to train and test various sophisticated machine learning algorithms. By using only essential dark matter halo physical properties and no merger history, our model predicts the gas mass, stellar mass, black hole mass, star formation rate, $g-r$ color, and stellar metallicity fairly robustly. Our results provide a unique and powerful phenomenological framework to explore the galaxy-halo connection that is built upon a solid hydrodynamical simulation. The promising reproduction of the listed galaxy properties demonstrably place ML as a promising and a significantly more computationally efficient tool to study small-scale structure formation. We find that ML mimics a full-blown hydrodynamical simulation surprisingly well in a computation time of mere minutes. The population of galaxies simulated by ML, while not numerically identical to Illustris, is statistically and physically robust and follows the same fundamental observational constraints. Machine learning offers an intriguing and promising technique to create quick mock galaxy catalogs in the future. △ Less

Submitted 12 January, 2016; v1 submitted 26 October, 2015; originally announced October 2015.

Comments: 20 pages, 27 figures, 6 tables. Accepted to MNRAS

Machine Learning and Cosmological Simulations I: Semi-Analytical Models

Authors: Harshil M. Kamdar, Matthew J. Turk, Robert J. Brunner

Abstract: We present a new exploratory framework to model galaxy formation and evolution in a hierarchical universe by using machine learning (ML). Our motivations are two-fold: (1) presenting a new, promising technique to study galaxy formation, and (2) quantitatively analyzing the extent of the influence of dark matter halo properties on galaxies in the backdrop of semi-analytical models (SAMs). We use th… ▽ More We present a new exploratory framework to model galaxy formation and evolution in a hierarchical universe by using machine learning (ML). Our motivations are two-fold: (1) presenting a new, promising technique to study galaxy formation, and (2) quantitatively analyzing the extent of the influence of dark matter halo properties on galaxies in the backdrop of semi-analytical models (SAMs). We use the influential Millennium Simulation and the corresponding Munich SAM to train and test various sophisticated machine learning algorithms (k-Nearest Neighbors, decision trees, random forests and extremely randomized trees). By using only essential dark matter halo physical properties for haloes of $M>10^{12} M_{\odot}$ and a partial merger tree, our model predicts the hot gas mass, cold gas mass, bulge mass, total stellar mass, black hole mass and cooling radius at z = 0 for each central galaxy in a dark matter halo for the Millennium run. Our results provide a unique and powerful phenomenological framework to explore the galaxy-halo connection that is built upon SAMs and demonstrably place ML as a promising and a computationally efficient tool to study small-scale structure formation. △ Less

Submitted 21 October, 2015; originally announced October 2015.

Comments: Accepted for publication in MNRAS. 19 pages, 20 figures, 4 tables

Journal ref: MNRAS Vol. 455 642-658 (2016)

Galaxy clustering, photometric redshifts and diagnosis of systematics in the DES Science Verification data

Authors: M. Crocce, J. Carretero, A. H. Bauer, A. J. Ross, I. Sevilla-Noarbe, T. Giannantonio, F. Sobreira, J. Sanchez, E. Gaztanaga, M. Carrasco Kind, C. Sanchez, C. Bonnett, A. Benoit-Levy, R. J. Brunner, A. Carnero Rosell, R. Cawthon, P. Fosalba, W. Hartley, E. J. Kim, B. Leistedt, R. Miquel, H. V. Peiris, W. J. Percival, R. Rosenfeld, E. S. Rykoff , et al. (62 additional authors not shown)

Abstract: We study the clustering of galaxies detected at $i<22.5$ in the Science Verification observations of the Dark Energy Survey (DES). Two-point correlation functions are measured using $2.3\times 10^6$ galaxies over a contiguous 116 deg$^2$ region in five bins of photometric redshift width $Δz = 0.2$ in the range $0.2 < z < 1.2.$ The impact of photometric redshift errors are assessed by comparing res… ▽ More We study the clustering of galaxies detected at $i<22.5$ in the Science Verification observations of the Dark Energy Survey (DES). Two-point correlation functions are measured using $2.3\times 10^6$ galaxies over a contiguous 116 deg$^2$ region in five bins of photometric redshift width $Δz = 0.2$ in the range $0.2 < z < 1.2.$ The impact of photometric redshift errors are assessed by comparing results using a template-based photo-$z$ algorithm (BPZ) to a machine-learning algorithm (TPZ). A companion paper (Leistedt et al 2015) presents maps of several observational variables (e.g. seeing, sky brightness) which could modulate the galaxy density. Here we characterize and mitigate systematic errors on the measured clustering which arise from these observational variables, in addition to others such as Galactic dust and stellar contamination. After correcting for systematic effects we measure galaxy bias over a broad range of linear scales relative to mass clustering predicted from the Planck $Λ$CDM model, finding agreement with CFHTLS measurements with $χ^2$ of 4.0 (8.7) with 5 degrees of freedom for the TPZ (BPZ) redshifts. We test a "linear bias" model, in which the galaxy clustering is a fixed multiple of the predicted non-linear dark-matter clustering. The precision of the data allow us to determine that the linear bias model describes the observed galaxy clustering to $2.5\%$ accuracy down to scales at least $4$ to $10$ times smaller than those on which linear theory is expected to be sufficient. △ Less

Submitted 15 December, 2015; v1 submitted 19 July, 2015; originally announced July 2015.

Comments: 23 pages, 18 figures, matches the version published in MNRAS. MNRAS 455, 4301-4324 (2015)

Robust micro-magnet design for fast electrical manipulations of single spins in quantum dots

Authors: Jun Yoneda, Tomohiro Otsuka, Tatsuki Takakura, Michel Pioro-Ladrière, Roland Brunner, Hong Lu, Takashi Nakajima, Toshiaki Obata, Akito Noiri, Christopher J. Palmstrøm, Arthur C. Gossard, Seigo Tarucha

Abstract: Tailoring spin coupling to electric fields is central to spintronics and spin-based quantum information processing. We present an optimal micromagnet design that produces appropriate stray magnetic fields to mediate fast electrical spin manipulations in nanodevices. We quantify the practical requirements for spatial field inhomogeneity and tolerance for misalignment with spins, and propose a desig… ▽ More Tailoring spin coupling to electric fields is central to spintronics and spin-based quantum information processing. We present an optimal micromagnet design that produces appropriate stray magnetic fields to mediate fast electrical spin manipulations in nanodevices. We quantify the practical requirements for spatial field inhomogeneity and tolerance for misalignment with spins, and propose a design scheme to improve the spin-rotation frequency (to exceed 50MHz in GaAs nanostructures). We then validate our design by experiments in separate devices. Our results will open a route to rapidly control solid-state electron spins with limited lifetimes and to study coherent spin dynamics in solids. △ Less

Submitted 7 July, 2015; originally announced July 2015.

Journal ref: Appl. Phys. Express 8, 084401 (2015)

A Hybrid Ensemble Learning Approach to Star-Galaxy Classification

Authors: Edward J. Kim, Robert J. Brunner, Matias Carrasco Kind

Abstract: There exist a variety of star-galaxy classification techniques, each with their own strengths and weaknesses. In this paper, we present a novel meta-classification framework that combines and fully exploits different techniques to produce a more robust star-galaxy classification. To demonstrate this hybrid, ensemble approach, we combine a purely morphological classifier, a supervised machine learn… ▽ More There exist a variety of star-galaxy classification techniques, each with their own strengths and weaknesses. In this paper, we present a novel meta-classification framework that combines and fully exploits different techniques to produce a more robust star-galaxy classification. To demonstrate this hybrid, ensemble approach, we combine a purely morphological classifier, a supervised machine learning method based on random forest, an unsupervised machine learning method based on self-organizing maps, and a hierarchical Bayesian template fitting method. Using data from the CFHTLenS survey, we consider different scenarios: when a high-quality training set is available with spectroscopic labels from DEEP2, SDSS, VIPERS, and VVDS, and when the demographics of sources in a low-quality training set do not match the demographics of objects in the test data set. We demonstrate that our Bayesian combination technique improves the overall performance over any individual classification method in these scenarios. Thus, strategies that combine the predictions of different classifiers may prove to be optimal in currently ongoing and forthcoming photometric surveys, such as the Dark Energy Survey and the Large Synoptic Survey Telescope. △ Less

Submitted 14 July, 2015; v1 submitted 8 May, 2015; originally announced May 2015.

Comments: 15 pages, 18 figures. Accepted for publication in MNRAS. Code available at https://github.com/EdwardJKim/astroclass

Journal ref: MNRAS 2015 453 (1): 507-521

On the Clustering of Compact Galaxy Pairs in Dark Matter Haloes

Authors: Yiran Wang, R. J. Brunner

Abstract: We analyze the clustering of photometrically selected galaxy pairs by using the halo-occupation distribution (HOD) model. We measure the angular two-point auto-correlation function, $ω(θ)$, for galaxies and galaxy pairs in three volume-limited samples and develop an HOD to model their clustering. Our results are successfully fit by these HOD models, and we see the separation of "1-halo" and "2-hal… ▽ More We analyze the clustering of photometrically selected galaxy pairs by using the halo-occupation distribution (HOD) model. We measure the angular two-point auto-correlation function, $ω(θ)$, for galaxies and galaxy pairs in three volume-limited samples and develop an HOD to model their clustering. Our results are successfully fit by these HOD models, and we see the separation of "1-halo" and "2-halo" clustering terms for both single galaxies and galaxy pairs. Our clustering measurements and HOD model fits for the single galaxy samples are consistent with previous results. We find that the galaxy pairs generally have larger clustering amplitudes than single galaxies, and the quantities computed during the HOD fitting, e.g., effective halo mass, $M_{eff}$, and linear bias, $b_{g}$, are also larger for galaxy pairs. We find that the central fractions for galaxy pairs are significantly higher than single galaxies, which confirms that galaxy pairs are formed at the center of more massive dark matter haloes. We also model the clustering dependence of the galaxy pair correlation function on redshift, galaxy type, and luminosity. We find early-early pairs (bright galaxy pairs) cluster more strongly than late-late pairs (dim galaxy pairs), and that the clustering does not depend on the luminosity contrast between the two galaxies in the compact group. △ Less

Submitted 30 July, 2014; originally announced July 2014.

Comments: Accepted by MNRAS, 17 pages, 12 figures

Photometric redshift analysis in the Dark Energy Survey Science Verification data

Authors: C. Sánchez, M. Carrasco Kind, H. Lin, R. Miquel, F. B. Abdalla, A. Amara, M. Banerji, C. Bonnett, R. Brunner, D. Capozzi, A. Carnero, F. J. Castander, L. A. N. da Costa, C. Cunha, A. Fausti, D. Gerdes, N. Greisel, J. Gschwend, W. Hartley, S. Jouvel, O. Lahav, M. Lima, M. A. G. Maia, P. Martí, R. L. C. Ogando , et al. (52 additional authors not shown)

Abstract: We present results from a study of the photometric redshift performance of the Dark Energy Survey (DES), using the early data from a Science Verification (SV) period of observations in late 2012 and early 2013 that provided science-quality images for almost 200 sq.~deg.~at the nominal depth of the survey. We assess the photometric redshift performance using about 15000 galaxies with spectroscopic… ▽ More We present results from a study of the photometric redshift performance of the Dark Energy Survey (DES), using the early data from a Science Verification (SV) period of observations in late 2012 and early 2013 that provided science-quality images for almost 200 sq.~deg.~at the nominal depth of the survey. We assess the photometric redshift performance using about 15000 galaxies with spectroscopic redshifts available from other surveys. These galaxies are used, in different configurations, as a calibration sample, and photo-$z$'s are obtained and studied using most of the existing photo-$z$ codes. A weighting method in a multi-dimensional color-magnitude space is applied to the spectroscopic sample in order to evaluate the photo-$z$ performance with sets that mimic the full DES photometric sample, which is on average significantly deeper than the calibration sample due to the limited depth of spectroscopic surveys. Empirical photo-$z$ methods using, for instance, Artificial Neural Networks or Random Forests, yield the best performance in the tests, achieving core photo-$z$ resolutions $σ_{68} \sim 0.08$. Moreover, the results from most of the codes, including template fitting methods, comfortably meet the DES requirements on photo-$z$ performance, therefore, providing an excellent precedent for future DES data sets. △ Less

Submitted 14 October, 2014; v1 submitted 12 June, 2014; originally announced June 2014.

Comments: Published in MNRAS. This version accounts for minor comments in the journal review

Journal ref: MNRAS 445, 1482-1506 (2014)

Sparse Representation of Photometric Redshift PDFs: Preparing for Petascale Astronomy

Authors: M. Carrasco Kind, R. J. Brunner

Abstract: One of the consequences of entering the era of precision cosmology is the widespread adoption of photometric redshift probability density functions (PDFs). Both current and future photometric surveys are expected to obtain images of billions of distinct galaxies. As a result, storing and analyzing all of these PDFs will be non-trivial and even more severe if a survey plans to compute and store mul… ▽ More One of the consequences of entering the era of precision cosmology is the widespread adoption of photometric redshift probability density functions (PDFs). Both current and future photometric surveys are expected to obtain images of billions of distinct galaxies. As a result, storing and analyzing all of these PDFs will be non-trivial and even more severe if a survey plans to compute and store multiple different PDFs. In this paper we propose the use of a sparse basis representation to fully represent individual photo-$z$ PDFs. By using an Orthogonal Matching Pursuit algorithm and a combination of Gaussian and Voigt basis functions, we demonstrate how our approach is superior to a multi-Gaussian fitting, as we require approximately half of the parameters for the same fitting accuracy with the additional advantage that an entire PDF can be stored by using a 4-byte integer per basis function, and we can achieve better accuracy by increasing the number of bases. By using data from the CFHTLenS, we demonstrate that only ten to twenty points per galaxy are sufficient to reconstruct both the individual PDFs and the ensemble redshift distribution, $N(z)$, to an accuracy of 99.9% when compared to the one built using the original PDFs computed with a resolution of $δz = 0.01$, reducing the required storage of two hundred original values by a factor of ten to twenty. Finally, we demonstrate how this basis representation can be directly extended to a cosmological analysis, thereby increasing computational performance without losing resolution nor accuracy. △ Less

Submitted 25 April, 2014; originally announced April 2014.

Comments: 12 pages, 10 figures. Accepted for publication in MNRAS. The code can be found at http://lcdm.astro.illinois.edu/code/pdfz.html

Exhausting the Information: Novel Bayesian Combination of Photometric Redshift PDFs

Authors: M. Carrasco Kind, R. J. Brunner

Abstract: The estimation and utilization of photometric redshift probability density functions (photo-$z$ PDFs) has become increasingly important over the last few years and currently there exist a wide variety of algorithms to compute photo-$z$'s, each with their own strengths and weaknesses. In this paper, we present a novel and efficient Bayesian framework that combines the results from different photo-… ▽ More The estimation and utilization of photometric redshift probability density functions (photo-$z$ PDFs) has become increasingly important over the last few years and currently there exist a wide variety of algorithms to compute photo-$z$'s, each with their own strengths and weaknesses. In this paper, we present a novel and efficient Bayesian framework that combines the results from different photo-$z$ techniques into a more powerful and robust estimate by maximizing the information from the photometric data. To demonstrate this we use a supervised machine learning technique based on random forest, an unsupervised method based on self-organizing maps, and a standard template fitting method but can be easily extend to other existing techniques. We use data from the DEEP2 and the SDSS surveys to explore different methods for combining the predictions from these techniques. By using different performance metrics, we demonstrate that we can improve the accuracy of our final photo-$z$ estimate over the best input technique, that the fraction of outliers is reduced, and that the identification of outliers is significantly improved when we apply a Naïve Bayes Classifier to this combined information. Our more robust and accurate photo-$z$ PDFs will allow even more precise cosmological constraints to be made by using current and future photometric surveys. These improvements are crucial as we move to analyze photometric data that push to or even past the limits of the available training data, which will be the case with the Large Synoptic Survey Telescope. △ Less

Submitted 4 June, 2014; v1 submitted 28 February, 2014; originally announced March 2014.

Comments: 21 pages, 19 figures, minor corrections, accepted for publication to MNRAS

SOMz: photometric redshift PDFs with self organizing maps and random atlas

Authors: M. Carrasco Kind, R. J. Brunner

Abstract: In this paper we explore the applicability of the unsupervised machine learning technique of Self Organizing Maps (SOM) to estimate galaxy photometric redshift probability density functions (PDFs). This technique takes a spectroscopic training set, and maps the photometric attributes, but not the redshifts, to a two dimensional surface by using a process of competitive learning where neurons compe… ▽ More In this paper we explore the applicability of the unsupervised machine learning technique of Self Organizing Maps (SOM) to estimate galaxy photometric redshift probability density functions (PDFs). This technique takes a spectroscopic training set, and maps the photometric attributes, but not the redshifts, to a two dimensional surface by using a process of competitive learning where neurons compete to more closely resemble the training data multidimensional space. The key feature of a SOM is that it retains the topology of the input set, revealing correlations between the attributes that are not easily identified. We test three different 2D topological mapping: rectangular, hexagonal, and spherical, by using data from the DEEP2 survey. We also explore different implementations and boundary conditions on the map and also introduce the idea of a random atlas where a large number of different maps are created and their individual predictions are aggregated to produce a more robust photometric redshift PDF. We also introduced a new metric, the $I$-score, which efficiently incorporates different metrics, making it easier to compare different results (from different parameters or different photometric redshift codes). We find that by using a spherical topology mapping we obtain a better representation of the underlying multidimensional topology, which provides more accurate results that are comparable to other, state-of-the-art machine learning algorithms. Our results illustrate that unsupervised approaches have great potential for many astronomical problems, and in particular for the computation of photometric redshifts. △ Less

Submitted 18 December, 2013; originally announced December 2013.

Comments: 14 pages, 8 figures. Accepted for publication in MNRAS. The code can be found at http://lcdm.astro.illinois.edu/research/SOMZ.html

Spectroscopic Needs for Imaging Dark Energy Experiments: Photometric Redshift Training and Calibration

Authors: J. Newman, A. Abate, F. Abdalla, S. Allam, S. Allen, R. Ansari, S. Bailey, W. Barkhouse, T. Beers, M. Blanton, M. Brodwin, J. Brownstein, R. Brunner, M. Carrasco-Kind, J. Cervantes-Cota, E. Chisari, M. Colless, J. Comparat, J. Coupon, E. Cheu, C. Cunha, A. de la Macorra, I. Dell'Antonio, B. Frye, E. Gawiser , et al. (36 additional authors not shown)

Abstract: Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training mak… ▽ More Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments. △ Less

Submitted 20 September, 2013; originally announced September 2013.

Comments: White paper for the "Dark Energy and CMB" working group for the American Physical Society's Division of Particles and Fields long-term planning exercise ("Snowmass")

Narrow absorption line variability in repeat quasar observations from the Sloan Digital Sky Survey

Authors: Troy L. Hacker, Robert J. Brunner, Britt F. Lundgren, Donald G. York

Abstract: We present the results from a time domain study of absorption lines detected in quasar spectra with repeat observations from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Beginning with over 4500 unique time separation baselines of various absorption line species identified in the SDSS DR7 quasar spectra, we create a catalogue of 2522 quasar absorption line systems with two to eight repe… ▽ More We present the results from a time domain study of absorption lines detected in quasar spectra with repeat observations from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Beginning with over 4500 unique time separation baselines of various absorption line species identified in the SDSS DR7 quasar spectra, we create a catalogue of 2522 quasar absorption line systems with two to eight repeat observations, representing the largest collection of unbiased and homogeneous multi-epoch absorption systems ever published. To investigate these systems for time variability of narrow absorption lines, we refine this sample based on the reliability of the system detection, the proximity of pixels with bright sky contamination to individual absorption lines, and the quality of the continuum fit. Variability measurements of this sub-sample based on the absorption line equivalent widths yield a total of 33 systems with indications of significantly variable absorption strengths on time-scales ranging from one day to several years in the rest frame of the absorption system. Of these, at least 10 are from a class known as intervening absorption systems caused by foreground galaxies along the line of sight to the background quasar. This is the first evidence of possible absorption line variability detected in intervening systems, and their short time-scale variations suggest that small-scale structures (~10-100 au) are likely to exist in their host foreground galaxies. △ Less

Submitted 30 July, 2013; originally announced July 2013.

Comments: 77 pages, 52 figures, accepted for publication in MNRAS

Dark energy with gravitational lens time delays

Authors: T. Treu, P. J. Marshall, F. -Y. Cyr-Racine, C. D. Fassnacht, C. R. Keeton, E. V. Linder, L. A. Moustakas, M. Bradac, E. Buckley-Geer, T. Collett, F. Courbin, G. Dobler, D. A. Finley, J. Hjorth, C. S. Kochanek, E. Komatsu, L. V. E. Koopmans, G. Meylan, P. Natarajan, M. Oguri, S. H. Suyu, M. Tewes, K. C. Wong, A. I. Zabludoff, D. Zaritsky , et al. (13 additional authors not shown)

Abstract: Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays… ▽ More Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes. △ Less

Submitted 5 June, 2013; originally announced June 2013.

Comments: White paper submitted to SNOWMASS2013

TPZ : Photometric redshift PDFs and ancillary information by using prediction trees and random forests

Authors: M. Carrasco Kind, R. J. Brunner

Abstract: With the growth of large photometric surveys, accurately estimating photometric redshifts, preferably as a probability density function (PDF), and fully understanding the implicit systematic uncertainties in this process has become increasingly important. In this paper, we present a new, publicly available, parallel, machine learning algorithm that generates photometric redshift PDFs by using pred… ▽ More With the growth of large photometric surveys, accurately estimating photometric redshifts, preferably as a probability density function (PDF), and fully understanding the implicit systematic uncertainties in this process has become increasingly important. In this paper, we present a new, publicly available, parallel, machine learning algorithm that generates photometric redshift PDFs by using prediction trees and random forest techniques, which we have named TPZ. This new algorithm incorporates measurement errors into the calculation while also dealing efficiently with missing values in the data. In addition, our implementation of this algorithm provides supplementary information regarding the data being analyzed, including unbiased estimates of the accuracy of the technique without resorting to a validation data set, identification of poor photometric redshift areas within the parameter space occupied by the spectroscopic training data, a quantification of the relative importance of the variables used to construct the PDF, and a robust identification of outliers. This extra information can be used to optimally target new spectroscopic observations and to improve the overall efficacy of the redshift estimation. We have tested TPZ on galaxy samples drawn from the SDSS main galaxy sample and from the DEEP2 survey, obtaining excellent results in each case. We also have tested our implementation by participating in the PHAT1 project, which is a blind photometric redshift contest, finding that TPZ performs comparable to if not better than other empirical photometric redshift algorithms. Finally, we discuss the various parameters that control the operation of TPZ, the specific limitations of this approach and an application of photometric redshift PDFs. △ Less

Submitted 28 March, 2013; originally announced March 2013.

Comments: 21 pages, 15 figures, Accepted for publication in MNRAS. TPZ code at http://lcdm.astro.illinois.edu/research/TPZ.html

The SDSS Galaxy Angular Two-Point Correlation Function

Authors: Y. Wang, R. J. Brunner, J. C. Dolence

Abstract: We present the galaxy two-point angular correlation function for galaxies selected from the seventh data release of the Sloan Digital Sky Survey. The galaxy sample was selected with $r$-band apparent magnitudes between 17 and 21; and we measure the correlation function for the full sample as well as for the four magnitude ranges: 17-18, 18-19, 19-20, and 20-21. We update the flag criteria to selec… ▽ More We present the galaxy two-point angular correlation function for galaxies selected from the seventh data release of the Sloan Digital Sky Survey. The galaxy sample was selected with $r$-band apparent magnitudes between 17 and 21; and we measure the correlation function for the full sample as well as for the four magnitude ranges: 17-18, 18-19, 19-20, and 20-21. We update the flag criteria to select a clean galaxy catalog and detail specific tests that we perform to characterize systematic effects, including the effects of seeing, Galactic extinction, and the overall survey uniformity. Notably, we find that optimally we can use observed regions with seeing $< 1\farcs5$, and $r$-band extinction < 0.13 magnitudes, smaller than previously published results. Furthermore, we confirm that the uniformity of the SDSS photometry is minimally affected by the stripe geometry. We find that, overall, the two-point angular correlation function can be described by a power law, $ω(θ) = A_ωθ^{(1-γ)}$ with $γ\simeq 1.72$, over the range $0\fdg005$--$10\degr$. We also find similar relationships for the four magnitude subsamples, but the amplitude within the same angular interval for the four subsamples is found to decrease with fainter magnitudes, in agreement with previous results. We find that the systematic signals are well below the galaxy angular correlation function for angles less than approximately $5\degr$, which limits the modeling of galaxy angular correlations on larger scales. Finally, we present our custom, highly parallelized two-point correlation code that we used in this analysis. △ Less

Submitted 12 March, 2013; v1 submitted 11 March, 2013; originally announced March 2013.

Comments: 22 pages, 17 figures, accepted by MNRAS

Bring out your codes! Bring out your codes! (Increasing Software Visibility and Re-use)

Authors: Alice Allen, Bruce Berriman, Robert Brunner, Dan Burger, Kimberly DuPrie, Robert J. Hanisch, Robert Mann, Jessica Mink, Christer Sandin, Keith Shortridge, Peter Teuben

Abstract: Progress is being made in code discoverability and preservation, but as discussed at ADASS XXI, many codes still remain hidden from public view. With the Astrophysics Source Code Library (ASCL) now indexed by the SAO/NASA Astrophysics Data System (ADS), the introduction of a new journal, Astronomy & Computing, focused on astrophysics software, and the increasing success of education efforts such a… ▽ More Progress is being made in code discoverability and preservation, but as discussed at ADASS XXI, many codes still remain hidden from public view. With the Astrophysics Source Code Library (ASCL) now indexed by the SAO/NASA Astrophysics Data System (ADS), the introduction of a new journal, Astronomy & Computing, focused on astrophysics software, and the increasing success of education efforts such as Software Carpentry and SciCoder, the community has the opportunity to set a higher standard for its science by encouraging the release of software for examination and possible reuse. We assembled representatives of the community to present issues inhibiting code release and sought suggestions for tackling these factors. The session began with brief statements by panelists; the floor was then opened for discussion and ideas. Comments covered a diverse range of related topics and points of view, with apparent support for the propositions that algorithms should be readily available, code used to produce published scientific results should be made available, and there should be discovery mechanisms to allow these to be found easily. With increased use of resources such as GitHub (for code availability), ASCL (for code discovery), and a stated strong preference from the new journal Astronomy & Computing for code release, we expect to see additional progress over the next few years. △ Less

Submitted 9 December, 2012; originally announced December 2012.

Comments: Birds of a Feather session at ADASS XXII (Champaign, IL; November, 2012) for proceedings; 4 pages. Organized by the Astrophysics Source Code Library (ASCL), which is available at ascl.net Unedited notes taken at the session are available here: http://asterisk.apod.com/wp/?p=192

The SDSS DR7 Galaxy Angular Power Spectrum: Volume-Limits and Galaxy Morphology

Authors: Brett Hayes, Robert Brunner

Abstract: We use a quadratic estimator with KL-compression to calculate the angular power spectrum of a volume-limited Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample out to l = 200. We also determine the angular power spectrum of selected subsamples with photometric redshifts z < 0.3 and 0.3 < z < 0.4 to examine the possible evolution of the angular power spectrum, as well as early-type… ▽ More We use a quadratic estimator with KL-compression to calculate the angular power spectrum of a volume-limited Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample out to l = 200. We also determine the angular power spectrum of selected subsamples with photometric redshifts z < 0.3 and 0.3 < z < 0.4 to examine the possible evolution of the angular power spectrum, as well as early-type and late-type galaxy subsamples to examine the relative linear bias. In addition, we calculate the angular power spectrum of the SDSS DR7 main galaxy sample in a ~ 53.7 square degree area out to l = 1600 to determine the SDSS DR7 angular power spectrum to high multipoles. We perform a χ^2 fit to compare the resulting angular power spectra to theoretical nonlinear angular power spectra to extract cosmological parameters and the linear bias. We find the best-fit cosmological parameters of Ω_m = 0.267 +- 0.038 and Ω_b = 0.045 +- 0.012. We find an overall linear bias of b = 1.075 +- 0.056, an early-type bias of b_e = 1.727 +- 0.065, and a late-type bias of b_l = 1.256 +- 0.051. Finally, we present evidence of a selective misclassification of late-type galaxies as stars by the SDSS photometric data reduction pipeline in areas of high stellar density (e.g., at low Galactic latitudes). △ Less

Submitted 6 November, 2012; originally announced November 2012.

Comments: 11 pages, 12 figures, accepted for publication in MNRAS

The SDSS DR7 Galaxy Angular Power Spectrum

Authors: Brett Hayes, Robert Brunner, Ashley Ross

Abstract: We calculate the angular power spectrum of galaxies selected from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) by using a quadratic estimation method with KL-compression. The primary data sample includes over 18 million galaxies covering more than 5,700 square degrees after masking areas with bright objects, reddening greater than 0.2 magnitudes, and seeing of more than 1.5 arcseconds.… ▽ More We calculate the angular power spectrum of galaxies selected from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) by using a quadratic estimation method with KL-compression. The primary data sample includes over 18 million galaxies covering more than 5,700 square degrees after masking areas with bright objects, reddening greater than 0.2 magnitudes, and seeing of more than 1.5 arcseconds. We test for systematic effects by calculating the angular power spectrum by SDSS stripe and find that these measurements are minimally affected by seeing and reddening. We calculate the angular power spectrum for l \leq 200 multipoles by using 40 bandpowers for the full sample, and l \leq 1000 multipoles using 50 bandpowers for individual stripes. We also calculate the angular power spectrum for this sample separated into 3 magnitude bins with mean redshifts of z = 0.171, z = 0.217, and z = 0.261 to examine the evolution of the angular power spectrum. We determine the theoretical linear angular power spectrum by projecting the 3D power spectrum to two dimensions for a basic comparison to our observational results. By minimizing the χ^2 fit between these data and the theoretical linear angular power spectrum we measure a loosely-constrained fit of Ω_m = 0.31^{+0.18}_{-0.11} with a linear bias of b = 0.94 \pm 0.04. △ Less

Submitted 24 December, 2011; originally announced December 2011.

Two-qubit Gate of Combined Single Spin Rotation and Inter-dot Spin Exchange in a Double Quantum Dot

Authors: R. Brunner, Y. -S. Shin, T. Obata, M. Pioro-Ladrière, T. Kubo, K. Yoshida, T. Taniyama, Y. Tokura, S. Tarucha

Abstract: A crucial requirement for quantum information processing is the realization of multiple-qubit quantum gates. Here, we demonstrate an electron spin based all-electrical two-qubit gate consisting of single spin rotations and inter-dot spin exchange in a double quantum dot. A partially entangled output state is obtained by the application of the two-qubit gate to an initial, uncorrelated state. We fi… ▽ More A crucial requirement for quantum information processing is the realization of multiple-qubit quantum gates. Here, we demonstrate an electron spin based all-electrical two-qubit gate consisting of single spin rotations and inter-dot spin exchange in a double quantum dot. A partially entangled output state is obtained by the application of the two-qubit gate to an initial, uncorrelated state. We find that the degree of entanglement is controllable by the exchange operation time. The approach represents a key step towards the realization of universal multiple qubit gates. △ Less

Submitted 15 September, 2011; originally announced September 2011.

Comments: accepted for publication, in press Phys.Rev.Lett. 2011

Evolution of the Clustering of Photometrically Selected SDSS Galaxies

Authors: Ashley J. Ross, Will J. Percival, Robert J. Brunner

Abstract: We measure the angular auto-correlation functions (w) of SDSS galaxies selected to have photometric redshifts 0.1 < z < 0.4 and absolute r-band magnitudes Mr < -21.2. We split these galaxies into five overlapping redshift shells of width 0.1 and measure w in each subsample in order to investigate the evolution of SDSS galaxies. We find that the bias increases substantially with redshift - much m… ▽ More We measure the angular auto-correlation functions (w) of SDSS galaxies selected to have photometric redshifts 0.1 < z < 0.4 and absolute r-band magnitudes Mr < -21.2. We split these galaxies into five overlapping redshift shells of width 0.1 and measure w in each subsample in order to investigate the evolution of SDSS galaxies. We find that the bias increases substantially with redshift - much more so than one would expect for a passively evolving sample. We use halo-model analysis to determine the best-fit halo-occupation-distribution (HOD) for each subsample, and the best-fit models allow us to interpret the change in bias physically. In order to properly interpret our best-fit HODs, we convert each halo mass to its z = 0 passively evolved bias (bo), enabling a direct comparison of the best-fit HODs at different redshifts. We find that the minimum halo bo required to host a galaxy decreases as the redshift decreases, suggesting that galaxies with Mr < -21.2 are forming in halos at the low-mass end of the HODs over our redshift range. We use the best-fit HODs to determine the change in occupation number divided by the change in mass of halos with constant bo and we find a sharp peak at bo ~ 0.9 - corresponding to an average halo mass of ~ 10^12Msol/h. We thus present the following scenario: the bias of galaxies with Mr < -21.2 decreases as the Universe evolves because these galaxies form in halos of mass ~ 10^12Msol/h (independent of redshift), and the bias of these halos naturally decreases as the Universe evolves. △ Less

Submitted 26 April, 2010; v1 submitted 8 February, 2010; originally announced February 2010.

Comments: 17 pages, 14 figures, matches version accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.407:420,2010

LSST Science Book, Version 2.0

Authors: LSST Science Collaboration, Paul A. Abell, Julius Allison, Scott F. Anderson, John R. Andrew, J. Roger P. Angel, Lee Armus, David Arnett, S. J. Asztalos, Tim S. Axelrod, Stephen Bailey, D. R. Ballantyne, Justin R. Bankert, Wayne A. Barkhouse, Jeffrey D. Barr, L. Felipe Barrientos, Aaron J. Barth, James G. Bartlett, Andrew C. Becker, Jacek Becla, Timothy C. Beers, Joseph P. Bernstein, Rahul Biswas, Michael R. Blanton, Joshua S. Bloom , et al. (223 additional authors not shown)

Abstract: A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south… ▽ More A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy. △ Less

Submitted 1 December, 2009; originally announced December 2009.

Comments: 596 pages. Also available at full resolution at http://www.lsst.org/lsst/scibook

Halo-model Analysis of the Clustering of Photometrically Selected Galaxies from SDSS

Authors: Ashley J Ross, Robert J. Brunner

Abstract: We measure the angular 2-point correlation functions of galaxies in a volume limited, photometrically selected galaxy sample from the fifth data release of the Sloan Digital Sky Survey. We split the sample both by luminosity and galaxy type and use a halo-model analysis to find halo-occupation distributions that can simultaneously model the clustering of all, early-, and late-type galaxies in a… ▽ More We measure the angular 2-point correlation functions of galaxies in a volume limited, photometrically selected galaxy sample from the fifth data release of the Sloan Digital Sky Survey. We split the sample both by luminosity and galaxy type and use a halo-model analysis to find halo-occupation distributions that can simultaneously model the clustering of all, early-, and late-type galaxies in a given sample. Our results for the full galaxy sample are generally consistent with previous results using the SDSS spectroscopic sample, taking the differences between the median redshifts of the photometric and spectroscopic samples into account. We find that our early- and late- type measurements cannot be fit by a model that allows early- and late-type galaxies to be well-mixed within halos. Instead, we introduce a new model that segregates early- and late-type galaxies into separate halos to the maximum allowed extent. We determine that, in all cases, it provides a good fit to our data and thus provides a new statistical description of the manner in which early- and late-type galaxies occupy halos. △ Less

Submitted 26 June, 2009; originally announced June 2009.

Comments: Accepted to MNRAS 11 pages, 6 figures

Journal ref: Mon.Not.Roy.Astron.Soc.399:878,2009

Data Mining and Machine Learning in Astronomy

Authors: Nicholas M. Ball, Robert J. Brunner

Abstract: We review the current state of data mining and machine learning in astronomy. 'Data Mining' can have a somewhat mixed connotation from the point of view of a researcher in this field. If used correctly, it can be a powerful approach, holding the potential to fully exploit the exponentially increasing amount of available data, promising great scientific advance. However, if misused, it can be littl… ▽ More We review the current state of data mining and machine learning in astronomy. 'Data Mining' can have a somewhat mixed connotation from the point of view of a researcher in this field. If used correctly, it can be a powerful approach, holding the potential to fully exploit the exponentially increasing amount of available data, promising great scientific advance. However, if misused, it can be little more than the black-box application of complex computing algorithms that may give little physical insight, and provide questionable results. Here, we give an overview of the entire data mining process, from data collection through to the interpretation of results. We cover common machine learning algorithms, such as artificial neural networks and support vector machines, applications from a broad range of astronomy, emphasizing those where data mining techniques directly resulted in improved science, and important current and future directions, including probability density functions, parallel algorithms, petascale computing, and the time domain. We conclude that, so long as one carefully selects an appropriate algorithm, and is guided by the astronomical problem at hand, data mining can be very much the powerful tool, and not the questionable black box. △ Less

Submitted 10 August, 2010; v1 submitted 11 June, 2009; originally announced June 2009.

Comments: Published in IJMPD. 61 pages, uses ws-ijmpd.cls. Several extra figures, some minor additions to the text

Journal ref: Int.J.Mod.Phys.D19:1049-1106,2010

Clustering of Low-Redshift (z <= 2.2) Quasars from the Sloan Digital Sky Survey

Authors: Nicholas P. Ross, Yue Shen, Michael A. Strauss, Daniel E. Vanden Berk, Andrew J. Connolly, Gordon T. Richards, Donald P. Schneider, David H. Weinberg, Patrick B. Hall, Neta A. Bahcall, Robert J. Brunner

Abstract: We present measurements of the quasar two-point correlation function, ξ_{Q}, over the redshift range z=0.3-2.2 based upon data from the SDSS. Using a homogeneous sample of 30,239 quasars with spectroscopic redshifts from the DR5 Quasar Catalogue, our study represents the largest sample used for this type of investigation to date. With this redshift range and an areal coverage of approx 4,000 deg… ▽ More We present measurements of the quasar two-point correlation function, ξ_{Q}, over the redshift range z=0.3-2.2 based upon data from the SDSS. Using a homogeneous sample of 30,239 quasars with spectroscopic redshifts from the DR5 Quasar Catalogue, our study represents the largest sample used for this type of investigation to date. With this redshift range and an areal coverage of approx 4,000 deg^2, we sample over 25 h^-3 Gpc^3 (comoving) assuming the current LCDM cosmology. Over this redshift range, we find that the redshift-space correlation function, xi(s), is adequately fit by a single power-law, with s_{0}=5.95+/-0.45 h^-1 Mpc and γ_{s}=1.16+0.11-0.16 when fit over s=1-25 h^-1 Mpc. Using the projected correlation function we calculate the real-space correlation length, r_{0}=5.45+0.35-0.45 h^-1 Mpc and γ=1.90+0.04-0.03, over scales of rp=1-130 h^-1 Mpc. Dividing the sample into redshift slices, we find very little, if any, evidence for the evolution of quasar clustering, with the redshift-space correlation length staying roughly constant at s_{0} ~ 6-7 h^-1 Mpc at z<2.2 (and only increasing at redshifts greater than this). Comparing our clustering measurements to those reported for X-ray selected AGN at z=0.5-1, we find reasonable agreement in some cases but significantly lower correlation lengths in others. We find that the linear bias evolves from b~1.4 at z=0.5 to b~3 at z=2.2, with b(z=1.27)=2.06+/-0.03 for the full sample. We compare our data to analytical models and infer that quasars inhabit dark matter haloes of constant mass M ~2 x 10^12 h^-1 M_Sol from redshifts z~2.5 (the peak of quasar activity) to z~0. [ABRIDGED] △ Less

Submitted 18 March, 2009; originally announced March 2009.

Comments: 28 pages, 26 figures, ApJ accepted. Online materials (including source code, catalogues and high-resolution figures) can be found at http://www.astro.psu.edu/users/npr/DR5/

A Cross-Correlation Analysis of Mg II Absorption Line Systems and Luminous Red Galaxies from the SDSS DR5

Authors: Britt F. Lundgren, Robert J. Brunner, Donald G. York, Ashley J. Ross, Jean M. Quashnock, Adam D. Myers, Donald P. Schneider, Yusra AlSayyad, Neta Bahcall

Abstract: We analyze the cross-correlation of 2,705 unambiguously intervening Mg II (2796,2803A) quasar absorption line systems with 1,495,604 luminous red galaxies (LRGs) from the Fifth Data Release of the Sloan Digital Sky Survey within the redshift range 0.36<=z<=0.8. We confirm with high precision a previously reported weak anti-correlation of equivalent width and dark matter halo mass, measuring the… ▽ More We analyze the cross-correlation of 2,705 unambiguously intervening Mg II (2796,2803A) quasar absorption line systems with 1,495,604 luminous red galaxies (LRGs) from the Fifth Data Release of the Sloan Digital Sky Survey within the redshift range 0.36<=z<=0.8. We confirm with high precision a previously reported weak anti-correlation of equivalent width and dark matter halo mass, measuring the average masses to be log M_h(M_[solar]h^-1)=11.29 [+0.36,-0.62] and log M_h(M_[solar]h^-1)=12.70 [+0.53,-1.16] for systems with W[2796A]>=1.4A and 0.8A<=W[2796A]<1.4A, respectively. Additionally, we investigate the significance of a number of potential sources of bias inherent in absorber-LRG cross-correlation measurements, including absorber velocity distributions and the weak lensing of background quasars, which we determine is capable of producing a 20-30% bias in angular cross-correlation measurements on scales less than 2'. We measure the Mg II - LRG cross-correlation for 719 absorption systems with v<60,000 km s^-1 in the quasar rest frame and find that these associated absorbers typically reside in dark matter haloes that are ~10-100 times more massive than those hosting unambiguously intervening Mg II absorbers. Furthermore, we find evidence for evolution of the redshift number density, dN/dz, with 2-sigma significance for the strongest (W>2.0A) absorbers in the DR5 sample. This width-dependent dN/dz evolution does not significantly affect the recovered equivalent width-halo mass anti-correlation and adds to existing evidence that the strongest Mg II absorption systems are correlated with an evolving population of field galaxies at z<0.8, while the non-evolving dN/dz of the weakest absorbers more closely resembles that of the LRG population. △ Less

Submitted 26 May, 2009; v1 submitted 23 February, 2009; originally announced February 2009.

Comments: 21 pages, 19 figures; Published in Astrophysical Journal

Journal ref: Astrophys.J.698:819-839,2009

The 2dF-SDSS LRG and QSO Survey: The spectroscopic QSO catalogue

Authors: Scott M. Croom, Gordon T. Richards, Tom Shanks, Brian J. Boyle, Robert G. Sharp, Joss Bland-Hawthorn, Terry Bridges, Robert J. Brunner, Russell Cannon, Daniel Carson, Kuenley Chiu, Matthew Colless, Warrick Couch, Roberto De Propris, Michael J. Drinkwater, Alastair Edge, Stephen Fine, Jon Loveday, Lance Miller, Adam D. Myers, Robert C. Nichol, Phil Outram, Kevin Pimbblet, Isaac Roseboom, Nicholas Ross , et al. (5 additional authors not shown)

Abstract: We present the final spectroscopic QSO catalogue from the 2dF-SDSS LRG and QSO (2SLAQ) Survey. This is a deep, 18<g<21.85 (extinction corrected), sample aimed at probing in detail the faint end of the broad line AGN luminosity distribution at z<2.6. The candidate QSOs were selected from SDSS photometry and observed spectroscopically with the 2dF spectrograph on the Anglo-Australian Telescope. Th… ▽ More We present the final spectroscopic QSO catalogue from the 2dF-SDSS LRG and QSO (2SLAQ) Survey. This is a deep, 18<g<21.85 (extinction corrected), sample aimed at probing in detail the faint end of the broad line AGN luminosity distribution at z<2.6. The candidate QSOs were selected from SDSS photometry and observed spectroscopically with the 2dF spectrograph on the Anglo-Australian Telescope. This sample covers an area of 191.9 deg^2 and contains new spectra of 16326 objects, of which 8764 are QSOs, and 7623 are newly discovered (the remainder were previously identified by the 2QZ and SDSS surveys). The full QSO sample (including objects previously observed in the SDSS and 2QZ surveys) contains 12702 QSOs. The new 2SLAQ spectroscopic data set also contains 2343 Galactic stars, including 362 white dwarfs, and 2924 narrow emission line galaxies with a median redshift of z=0.22. We present detailed completeness estimates for the survey, based on modelling of QSO colours, including host galaxy contributions. This calculation shows that at g~21.85 QSO colours are significantly affected by the presence of a host galaxy up to redshift z~1 in the SDSS ugriz bands. In particular we see a significant reddening of the objects in g-i towards fainter g-band magnitudes. This reddening is consistent with the QSO host galaxies being dominated by a stellar population of age at least 2-3 Gyr. The full catalogue, including completeness estimates, is available on-line at http://www.2slaq.info/ △ Less

Submitted 27 October, 2008; originally announced October 2008.

Comments: 28 pages, 23 figures. Accepted for publication in MNRAS. Survey data, including catalogue and spectra available from http://www.2slaq.info/

Journal ref: Mon.Not.Roy.Astron.Soc.392:19-44,2009

Quasar Clustering from SDSS DR5: Dependences on Physical Properties

Authors: Yue Shen, Michael A. Strauss, Nicholas P. Ross, Patrick B. Hall, Yen-Ting Lin, Gordon T. Richards, Donald P. Schneider, David H. Weinberg, Andrew J. Connolly, Xiaohui Fan, Joseph F. Hennawi, Francesco Shankar, Daniel E. Vanden Berk, Neta A. Bahcall, Robert J. Brunner

Abstract: Using a homogenous sample of 38,208 quasars with a sky coverage of $4000 {\rm deg^2}$ drawn from the SDSS Data Release Five quasar catalog, we study the dependence of quasar clustering on luminosity, virial black hole mass, quasar color, and radio loudness. At $z<2.5$, quasar clustering depends weakly on luminosity and virial black hole mass, with typical uncertainty levels $\sim 10%$ for the me… ▽ More Using a homogenous sample of 38,208 quasars with a sky coverage of $4000 {\rm deg^2}$ drawn from the SDSS Data Release Five quasar catalog, we study the dependence of quasar clustering on luminosity, virial black hole mass, quasar color, and radio loudness. At $z<2.5$, quasar clustering depends weakly on luminosity and virial black hole mass, with typical uncertainty levels $\sim 10%$ for the measured correlation lengths. These weak dependences are consistent with models in which substantial scatter between quasar luminosity, virial black hole mass and the host dark matter halo mass has diluted any clustering difference, where halo mass is assumed to be the relevant quantity that best correlates with clustering strength. However, the most luminous and most massive quasars are more strongly clustered (at the $\sim 2σ$ level) than the remainder of the sample, which we attribute to the rapid increase of the bias factor at the high-mass end of host halos. We do not observe a strong dependence of clustering strength on quasar colors within our sample. On the other hand, radio-loud quasars are more strongly clustered than are radio-quiet quasars matched in redshift and optical luminosity (or virial black hole mass), consistent with local observations of radio galaxies and radio-loud type 2 AGN. Thus radio-loud quasars reside in more massive and denser environments in the biased halo clustering picture. Using the Sheth et al.(2001) formula for the linear halo bias, the estimated host halo mass for radio-loud quasars is $\sim 10^{13} h^{-1}M_\odot$, compared to $\sim 2\times 10^{12} h^{-1}M_\odot$ for radio-quiet quasar hosts at $z\sim 1.5$. △ Less

Submitted 13 December, 2008; v1 submitted 22 October, 2008; originally announced October 2008.

Comments: Updated version; accepted for publication in ApJ

Eight-Dimensional Mid-Infrared/Optical Bayesian Quasar Selection

Authors: Gordon T. Richards, Rajesh P. Deo, Mark Lacy, Adam D. Myers, Robert C. Nichol, Nadia L. Zakamska, Robert J. Brunner, W. N. Brandt, Alexander G. Gray, John K. Parejko, Andrew Ptak, Donald P. Schneider, Lisa J. Storrie-Lombardi, Alexander S. Szalay

Abstract: We explore the multidimensional, multiwavelength selection of quasars from mid-IR (MIR) plus optical data, specifically from Spitzer-IRAC and the Sloan Digital Sky Survey (SDSS). We apply modern statistical techniques to combined Spitzer MIR and SDSS optical data, allowing up to 8-D color selection of quasars. Using a Bayesian selection method, we catalog 5546 quasar candidates to an 8.0 um dept… ▽ More We explore the multidimensional, multiwavelength selection of quasars from mid-IR (MIR) plus optical data, specifically from Spitzer-IRAC and the Sloan Digital Sky Survey (SDSS). We apply modern statistical techniques to combined Spitzer MIR and SDSS optical data, allowing up to 8-D color selection of quasars. Using a Bayesian selection method, we catalog 5546 quasar candidates to an 8.0 um depth of 56 uJy over an area of ~24 sq. deg; ~70% of these candidates are not identified by applying the same Bayesian algorithm to 4-color SDSS optical data alone. Our selection recovers 97.7% of known type 1 quasars in this area and greatly improves the effectiveness of identifying 3.5<z<5 quasars. Even using only the two shortest wavelength IRAC bandpasses, it is possible to use our Bayesian techniques to select quasars with 97% completeness and as little as 10% contamination. This sample has a photometric redshift accuracy of 93.6% (Delta Z +/-0.3), remaining roughly constant when the two reddest MIR bands are excluded. While our methods are designed to find type 1 (unobscured) quasars, as many as 1200 of the objects are type 2 (obscured) quasar candidates. Coupling deep optical imaging data with deep mid-IR data could enable selection of quasars in significant numbers past the peak of the quasar luminosity function (QLF) to at least z~4. Such a sample would constrain the shape of the QLF and enable quasar clustering studies over the largest range of redshift and luminosity to date, yielding significant gains in our understanding of quasars and the evolution of galaxies. △ Less

Submitted 25 February, 2009; v1 submitted 20 October, 2008; originally announced October 2008.

Comments: 49 pages, 14 figures, 7 tables. AJ, accepted

Journal ref: Astron.J.137:3884,2009