Showing 1–50 of 53 results for author: Stein, G

Anticipatory Task and Motion Planning

Authors: Roshan Dhakal, Duc M. Nguyen, Tom Silver, Xuesu Xiao, Gregory J. Stein

Abstract: We consider a sequential task and motion planning (tamp) setting in which a robot is assigned continuous-space rearrangement-style tasks one-at-a-time in an environment that persists between each. Lacking advance knowledge of future tasks, existing (myopic) planning strategies unwittingly introduce side effects that impede completion of subsequent tasks: e.g., by blocking future access or manipula… ▽ More We consider a sequential task and motion planning (tamp) setting in which a robot is assigned continuous-space rearrangement-style tasks one-at-a-time in an environment that persists between each. Lacking advance knowledge of future tasks, existing (myopic) planning strategies unwittingly introduce side effects that impede completion of subsequent tasks: e.g., by blocking future access or manipulation. We present anticipatory task and motion planning, in which estimates of expected future cost from a learned model inform selection of plans generated by a model-based tamp planner so as to avoid such side effects, choosing configurations of the environment that both complete the task and minimize overall cost. Simulated multi-task deployments in navigation-among-movable-obstacles and cabinet-loading domains yield improvements of 32.7% and 16.7% average per-task cost respectively. When given time in advance to prepare the environment, our learning-augmented planning approach yields improvements of 83.1% and 22.3%. Both showcase the value of our approach. Finally, we also demonstrate anticipatory tamp on a real-world Fetch mobile manipulator. △ Less

Submitted 18 July, 2024; originally announced July 2024.

Benchmarking Robust Self-Supervised Learning Across Diverse Downstream Tasks

Authors: Antoni Kowalczuk, Jan Dubiński, Atiyeh Ashari Ghomi, Yi Sui, George Stein, Jiapeng Wu, Jesse C. Cresswell, Franziska Boenisch, Adam Dziedzic

Abstract: Large-scale vision models have become integral in many applications due to their unprecedented performance and versatility across downstream tasks. However, the robustness of these foundation models has primarily been explored for a single task, namely image classification. The vulnerability of other common vision tasks, such as semantic segmentation and depth estimation, remains largely unknown.… ▽ More Large-scale vision models have become integral in many applications due to their unprecedented performance and versatility across downstream tasks. However, the robustness of these foundation models has primarily been explored for a single task, namely image classification. The vulnerability of other common vision tasks, such as semantic segmentation and depth estimation, remains largely unknown. We present a comprehensive empirical evaluation of the adversarial robustness of self-supervised vision encoders across multiple downstream tasks. Our attacks operate in the encoder embedding space and at the downstream task output level. In both cases, current state-of-the-art adversarial fine-tuning techniques tested only for classification significantly degrade clean and robust performance on other tasks. Since the purpose of a foundation model is to cater to multiple applications at once, our findings reveal the need to enhance encoder robustness more broadly. Our code is available at ${github.com/layer6ai-labs/ssl-robustness}$. △ Less

Submitted 18 July, 2024; v1 submitted 17 July, 2024; originally announced July 2024.

Comments: Accepted at the ICML 2024 Workshop on Foundation Models in the Wild

TabPFGen -- Tabular Data Generation with TabPFN

Authors: Junwei Ma, Apoorv Dankar, George Stein, Guangwei Yu, Anthony Caterini

Abstract: Advances in deep generative modelling have not translated well to tabular data. We argue that this is caused by a mismatch in structure between popular generative models and discriminative models of tabular data. We thus devise a technique to turn TabPFN -- a highly performant transformer initially designed for in-context discriminative tabular tasks -- into an energy-based generative model, which… ▽ More Advances in deep generative modelling have not translated well to tabular data. We argue that this is caused by a mismatch in structure between popular generative models and discriminative models of tabular data. We thus devise a technique to turn TabPFN -- a highly performant transformer initially designed for in-context discriminative tabular tasks -- into an energy-based generative model, which we dub TabPFGen. This novel framework leverages the pre-trained TabPFN as part of the energy function and does not require any additional training or hyperparameter tuning, thus inheriting TabPFN's in-context learning capability. We can sample from TabPFGen analogously to other energy-based models. We demonstrate strong results on standard generative modelling tasks, including data augmentation, class-balancing, and imputation, unlocking a new frontier of tabular data generation. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Analysis of the BraTS 2023 Intracranial Meningioma Segmentation Challenge

Authors: Dominic LaBella, Ujjwal Baid, Omaditya Khanna, Shan McBurney-Lin, Ryan McLean, Pierre Nedelec, Arif Rashid, Nourel Hoda Tahon, Talissa Altes, Radhika Bhalerao, Yaseen Dhemesh, Devon Godfrey, Fathi Hilal, Scott Floyd, Anastasia Janas, Anahita Fathi Kazerooni, John Kirkpatrick, Collin Kent, Florian Kofler, Kevin Leu, Nazanin Maleki, Bjoern Menze, Maxence Pajot, Zachary J. Reitman, Jeffrey D. Rudie , et al. (96 additional authors not shown)

Abstract: We describe the design and results from the BraTS 2023 Intracranial Meningioma Segmentation Challenge. The BraTS Meningioma Challenge differed from prior BraTS Glioma challenges in that it focused on meningiomas, which are typically benign extra-axial tumors with diverse radiologic and anatomical presentation and a propensity for multiplicity. Nine participating teams each developed deep-learning… ▽ More We describe the design and results from the BraTS 2023 Intracranial Meningioma Segmentation Challenge. The BraTS Meningioma Challenge differed from prior BraTS Glioma challenges in that it focused on meningiomas, which are typically benign extra-axial tumors with diverse radiologic and anatomical presentation and a propensity for multiplicity. Nine participating teams each developed deep-learning automated segmentation models using image data from the largest multi-institutional systematically expert annotated multilabel multi-sequence meningioma MRI dataset to date, which included 1000 training set cases, 141 validation set cases, and 283 hidden test set cases. Each case included T2, T2/FLAIR, T1, and T1Gd brain MRI sequences with associated tumor compartment labels delineating enhancing tumor, non-enhancing tumor, and surrounding non-enhancing T2/FLAIR hyperintensity. Participant automated segmentation models were evaluated and ranked based on a scoring system evaluating lesion-wise metrics including dice similarity coefficient (DSC) and 95% Hausdorff Distance. The top ranked team had a lesion-wise median dice similarity coefficient (DSC) of 0.976, 0.976, and 0.964 for enhancing tumor, tumor core, and whole tumor, respectively and a corresponding average DSC of 0.899, 0.904, and 0.871, respectively. These results serve as state-of-the-art benchmarks for future pre-operative meningioma automated segmentation algorithms. Additionally, we found that 1286 of 1424 cases (90.3%) had at least 1 compartment voxel abutting the edge of the skull-stripped image edge, which requires further investigation into optimal pre-processing face anonymization steps. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 16 pages, 11 tables, 10 figures, MICCAI

When Medical Imaging Met Self-Attention: A Love Story That Didn't Quite Work Out

Authors: Tristan Piater, Niklas Penzel, Gideon Stein, Joachim Denzler

Abstract: A substantial body of research has focused on developing systems that assist medical professionals during labor-intensive early screening processes, many based on convolutional deep-learning architectures. Recently, multiple studies explored the application of so-called self-attention mechanisms in the vision domain. These studies often report empirical improvements over fully convolutional approa… ▽ More A substantial body of research has focused on developing systems that assist medical professionals during labor-intensive early screening processes, many based on convolutional deep-learning architectures. Recently, multiple studies explored the application of so-called self-attention mechanisms in the vision domain. These studies often report empirical improvements over fully convolutional approaches on various datasets and tasks. To evaluate this trend for medical imaging, we extend two widely adopted convolutional architectures with different self-attention variants on two different medical datasets. With this, we aim to specifically evaluate the possible advantages of additional self-attention. We compare our models with similarly sized convolutional and attention-based baselines and evaluate performance gains statistically. Additionally, we investigate how including such layers changes the features learned by these models during the training. Following a hyperparameter search, and contrary to our expectations, we observe no significant improvement in balanced accuracy over fully convolutional models. We also find that important features, such as dermoscopic structures in skin lesion images, are still not learned by employing self-attention. Finally, analyzing local explanations, we confirm biased feature usage. We conclude that merely incorporating attention is insufficient to surpass the performance of existing fully convolutional methods. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 10 pages, 2 figures, 5 tables, presented at VISAPP 2024

Journal ref: Proceedings of the 19th International Joint Conference on Computer Vision, Imaging and Computer Graphics Theory and Applications - Volume 2: VISAPP (2024), ISBN 978-989-758-679-8, ISSN 2184-4321, SciTePress, pages 149-158

Reducing Bias in Pre-trained Models by Tuning while Penalizing Change

Authors: Niklas Penzel, Gideon Stein, Joachim Denzler

Abstract: Deep models trained on large amounts of data often incorporate implicit biases present during training time. If later such a bias is discovered during inference or deployment, it is often necessary to acquire new data and retrain the model. This behavior is especially problematic in critical areas such as autonomous driving or medical decision-making. In these scenarios, new data is often expensiv… ▽ More Deep models trained on large amounts of data often incorporate implicit biases present during training time. If later such a bias is discovered during inference or deployment, it is often necessary to acquire new data and retrain the model. This behavior is especially problematic in critical areas such as autonomous driving or medical decision-making. In these scenarios, new data is often expensive and hard to come by. In this work, we present a method based on change penalization that takes a pre-trained model and adapts the weights to mitigate a previously detected bias. We achieve this by tuning a zero-initialized copy of a frozen pre-trained network. Our method needs very few, in extreme cases only a single, examples that contradict the bias to increase performance. Additionally, we propose an early stopping criterion to modify baselines and reduce overfitting. We evaluate our approach on a well-known bias in skin lesion classification and three other datasets from the domain shift literature. We find that our approach works especially well with very few images. Simple fine-tuning combined with our early stopping also leads to performance benefits for a larger number of tuning samples. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 12 pages, 12 figures, presented at VISAPP 2024

Journal ref: Proceedings of the 19th International Joint Conference on Computer Vision (2024), Imaging and Computer Graphics Theory and Applications - Volume 2: VISAPP, ISBN 978-989-758-679-8, ISSN 2184-4321, SciTePress, pages 90-101

Team Coordination on Graphs: Problem, Analysis, and Algorithms

Authors: Manshi Limbu, Yanlin Zhou, Gregory Stein, Xuan Wang, Daigo Shishika, Xuesu Xiao

Abstract: Team Coordination on Graphs with Risky Edges (TCGRE) is a recently emerged problem, in which a robot team collectively reduces graph traversal cost through support from one robot to another when the latter traverses a risky edge. Resembling the traditional Multi-Agent Path Finding (MAPF) problem, both classical and learning-based methods have been proposed to solve TCGRE, however, they lacked eith… ▽ More Team Coordination on Graphs with Risky Edges (TCGRE) is a recently emerged problem, in which a robot team collectively reduces graph traversal cost through support from one robot to another when the latter traverses a risky edge. Resembling the traditional Multi-Agent Path Finding (MAPF) problem, both classical and learning-based methods have been proposed to solve TCGRE, however, they lacked either computation efficiency or optimality assurance. In this paper, we reformulate TCGRE as a constrained optimization and perform rigorous mathematical analysis. Our theoretical analysis shows the NP-hardness of TCGRE by reduction from the Maximum 3D Matching problem and that efficient decomposition is a key to tackle this combinatorial optimization problem. Further more, we design three classes of algorithms to solve TCGRE, i.e., Joint State Graph (JSG) based, coordination based, and receding-horizon sub-team based solutions. Each of these proposed algorithms enjoy different provable optimality and efficiency characteristics that are demonstrated in our extensive experiments. △ Less

Submitted 23 March, 2024; originally announced March 2024.

Comments: 8 pages, 4 figures

Active Information Gathering for Long-Horizon Navigation Under Uncertainty by Learning the Value of Information

Authors: Raihan Islam Arnob, Gregory J. Stein

Abstract: We address the task of long-horizon navigation in partially mapped environments for which active gathering of information about faraway unseen space is essential for good behavior. We present a novel planning strategy that, at training time, affords tractable computation of the value of information associated with revealing potentially informative regions of unseen space, data used to train a grap… ▽ More We address the task of long-horizon navigation in partially mapped environments for which active gathering of information about faraway unseen space is essential for good behavior. We present a novel planning strategy that, at training time, affords tractable computation of the value of information associated with revealing potentially informative regions of unseen space, data used to train a graph neural network to predict the goodness of temporally-extended exploratory actions. Our learning-augmented model-based planning approach predicts the expected value of information of revealing unseen space and is capable of using these predictions to actively seek information and so improve long-horizon navigation. Across two simulated office-like environments, our planner outperforms competitive learned and non-learned baseline navigation strategies, achieving improvements of up to 63.76% and 36.68%, demonstrating its capacity to actively seek performance-critical information. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: Submitted at IROS'24. arXiv admin note: text overlap with arXiv:2307.14501

Embracing the black box: Heading towards foundation models for causal discovery from time series data

Authors: Gideon Stein, Maha Shadaydeh, Joachim Denzler

Abstract: Causal discovery from time series data encompasses many existing solutions, including those based on deep learning techniques. However, these methods typically do not endorse one of the most prevalent paradigms in deep learning: End-to-end learning. To address this gap, we explore what we call Causal Pretraining. A methodology that aims to learn a direct mapping from multivariate time series to th… ▽ More Causal discovery from time series data encompasses many existing solutions, including those based on deep learning techniques. However, these methods typically do not endorse one of the most prevalent paradigms in deep learning: End-to-end learning. To address this gap, we explore what we call Causal Pretraining. A methodology that aims to learn a direct mapping from multivariate time series to the underlying causal graphs in a supervised manner. Our empirical findings suggest that causal discovery in a supervised manner is possible, assuming that the training and test time series samples share most of their dynamics. More importantly, we found evidence that the performance of Causal Pretraining can increase with data and model size, even if the additional data do not share the same dynamics. Further, we provide examples where causal discovery for real-world data with causally pretrained neural networks is possible within limits. We argue that this hints at the possibility of a foundation model for causal discovery. △ Less

Submitted 14 February, 2024; originally announced February 2024.

Comments: AAAI Workshop (AI4TS) 2024

MSC Class: 68T07

Self-supervised Representation Learning From Random Data Projectors

Authors: Yi Sui, Tongzi Wu, Jesse C. Cresswell, Ga Wu, George Stein, Xiao Shi Huang, Xiaochen Zhang, Maksims Volkovs

Abstract: Self-supervised representation learning~(SSRL) has advanced considerably by exploiting the transformation invariance assumption under artificially designed data augmentations. While augmentation-based SSRL algorithms push the boundaries of performance in computer vision and natural language processing, they are often not directly applicable to other data modalities, and can conflict with applicati… ▽ More Self-supervised representation learning~(SSRL) has advanced considerably by exploiting the transformation invariance assumption under artificially designed data augmentations. While augmentation-based SSRL algorithms push the boundaries of performance in computer vision and natural language processing, they are often not directly applicable to other data modalities, and can conflict with application-specific data augmentation constraints. This paper presents an SSRL approach that can be applied to any data modality and network architecture because it does not rely on augmentations or masking. Specifically, we show that high-quality data representations can be learned by reconstructing random data projections. We evaluate the proposed approach on a wide range of representation learning tasks that span diverse modalities and real-world applications. We show that it outperforms multiple state-of-the-art SSRL baselines. Due to its wide applicability and strong empirical results, we argue that learning from randomness is a fruitful research direction worthy of attention and further study. △ Less

Submitted 20 March, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Published as a conference paper of ICLR 2024. https://openreview.net/pdf?id=EpYnZpDpsQ

Improving Reliable Navigation under Uncertainty via Predictions Informed by Non-Local Information

Authors: Raihan Islam Arnob, Gregory J. Stein

Abstract: We improve reliable, long-horizon, goal-directed navigation in partially-mapped environments by using non-locally available information to predict the goodness of temporally-extended actions that enter unseen space. Making predictions about where to navigate in general requires non-local information: any observations the robot has seen so far may provide information about the goodness of a particu… ▽ More We improve reliable, long-horizon, goal-directed navigation in partially-mapped environments by using non-locally available information to predict the goodness of temporally-extended actions that enter unseen space. Making predictions about where to navigate in general requires non-local information: any observations the robot has seen so far may provide information about the goodness of a particular direction of travel. Building on recent work in learning-augmented model-based planning under uncertainty, we present an approach that can both rely on non-local information to make predictions (via a graph neural network) and is reliable by design: it will always reach its goal, even when learning does not provide accurate predictions. We conduct experiments in three simulated environments in which non-local information is needed to perform well. In our large scale university building environment, generated from real-world floorplans to the scale, we demonstrate a 9.3\% reduction in cost-to-go compared to a non-learned baseline and a 14.9\% reduction compared to a learning-informed planner that can only use local information to inform its predictions. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: IROS 2023

Guided Sampling-Based Motion Planning with Dynamics in Unknown Environments

Authors: Abhish Khanal, Hoang-Dung Bui, Gregory J. Stein, Erion Plaku

Abstract: Despite recent progress improving the efficiency and quality of motion planning, planning collision-free and dynamically-feasible trajectories in partially-mapped environments remains challenging, since constantly replanning as unseen obstacles are revealed during navigation both incurs significant computational expense and can introduce problematic oscillatory behavior. To improve the quality of… ▽ More Despite recent progress improving the efficiency and quality of motion planning, planning collision-free and dynamically-feasible trajectories in partially-mapped environments remains challenging, since constantly replanning as unseen obstacles are revealed during navigation both incurs significant computational expense and can introduce problematic oscillatory behavior. To improve the quality of motion planning in partial maps, this paper develops a framework that augments sampling-based motion planning to leverage a high-level discrete layer and prior solutions to guide motion-tree expansion during replanning, affording both (i) faster planning and (ii) improved solution coherence. Our framework shows significant improvements in runtime and solution distance when compared with other sampling-based motion planners. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: 8 Pages, 8 Figures, IEEE International Conference on Automation Science and Engineering (CASE) 2023

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

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

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

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

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

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

Anticipatory Planning: Improving Long-Lived Planning by Estimating Expected Cost of Future Tasks

Authors: Roshan Dhakal, Md Ridwan Hossain Talukder, Gregory J. Stein

Abstract: We consider a service robot in a household environment given a sequence of high-level tasks one at a time. Most existing task planners, lacking knowledge of what they may be asked to do next, solve each task in isolation and so may unwittingly introduce side effects that make subsequent tasks more costly. In order to reduce the overall cost of completing all tasks, we consider that the robot must… ▽ More We consider a service robot in a household environment given a sequence of high-level tasks one at a time. Most existing task planners, lacking knowledge of what they may be asked to do next, solve each task in isolation and so may unwittingly introduce side effects that make subsequent tasks more costly. In order to reduce the overall cost of completing all tasks, we consider that the robot must anticipate the impact its actions could have on future tasks. Thus, we propose anticipatory planning: an approach in which estimates of the expected future cost, from a graph neural network, augment model-based task planning. Our approach guides the robot towards behaviors that encourage preparation and organization, reducing overall costs in long-lived planning scenarios. We evaluate our method on blockworld environments and show that our approach reduces the overall planning costs by 5% as compared to planning without anticipatory planning. Additionally, if given an opportunity to prepare the environment in advance (a special case of anticipatory planning), our planner improves overall cost by 11%. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Exploring the Non-Gaussianity of the Cosmic Infrared Background and Its Weak Gravitational Lensing

Authors: Jaemyoung Lee, J. Richard Bond, Pavel Motloch, Alexander van Engelen, George Stein

Abstract: Gravitational lensing deflects the paths of photons, altering the statistics of cosmic backgrounds and distorting their information content. We take the Cosmic Infrared Background (CIB), which provides plentiful information about galaxy formation and evolution, as an example to probe the effect of lensing on non-Gaussian statistics. Using the Websky simulations, we first quantify the non-Gaussiani… ▽ More Gravitational lensing deflects the paths of photons, altering the statistics of cosmic backgrounds and distorting their information content. We take the Cosmic Infrared Background (CIB), which provides plentiful information about galaxy formation and evolution, as an example to probe the effect of lensing on non-Gaussian statistics. Using the Websky simulations, we first quantify the non-Gaussianity of the CIB, revealing additional detail on top of its well-measured power spectrum. To achieve this, we use needlet-like multipole-band-filters to calculate the variance and higher-point correlations. Using our simulations, we show the 2-point, 3-point and 4-point spectra, and compare our calculated power spectra and bispectra to Planck values. We then lens the CIB, shell-by-shell with corresponding convergence maps, to capture the broad redshift extent of both the CIB and its lensing convergence. The lensing of the CIB changes the 3-point and 4-point functions by a few tens of percent at large scales, unlike with the power spectrum, which changes by less than two percent. We expand our analyses to encompass the full intensity probability distribution functions (PDFs) involving all n-point correlations as a function of scale. In particular, we use the relative entropy between lensed and unlensed PDFs to create a spectrum of templates that can allow estimation of lensing. The underlying CIB model is missing the important role of star-bursting, which we test by adding a stochastic log-normal term to the intensity distributions. The novel aspects of our filtering and lensing pipeline should prove useful for any radiant background, including line intensity maps. △ Less

Submitted 23 February, 2024; v1 submitted 14 April, 2023; originally announced April 2023.

Comments: 16 pages, 16 figures, accepted by MNRAS

Data-Efficient Policy Selection for Navigation in Partial Maps via Subgoal-Based Abstraction

Authors: Abhishek Paudel, Gregory J. Stein

Abstract: We present a novel approach for fast and reliable policy selection for navigation in partial maps. Leveraging the recent learning-augmented model-based Learning over Subgoals Planning (LSP) abstraction to plan, our robot reuses data collected during navigation to evaluate how well other alternative policies could have performed via a procedure we call offline alt-policy replay. Costs from offline… ▽ More We present a novel approach for fast and reliable policy selection for navigation in partial maps. Leveraging the recent learning-augmented model-based Learning over Subgoals Planning (LSP) abstraction to plan, our robot reuses data collected during navigation to evaluate how well other alternative policies could have performed via a procedure we call offline alt-policy replay. Costs from offline alt-policy replay constrain policy selection among the LSP-based policies during deployment, allowing for improvements in convergence speed, cumulative regret and average navigation cost. With only limited prior knowledge about the nature of unseen environments, we achieve at least 67% and as much as 96% improvements on cumulative regret over the baseline bandit approach in our experiments in simulated maze and office-like environments. △ Less

Submitted 1 August, 2023; v1 submitted 3 April, 2023; originally announced April 2023.

Comments: 8 pages, 5 figures. Accepted at IROS 2023

Learning Augmented, Multi-Robot Long-Horizon Navigation in Partially Mapped Environments

Authors: Abhish Khanal, Gregory J. Stein

Abstract: We present a novel approach for efficient and reliable goal-directed long-horizon navigation for a multi-robot team in a structured, unknown environment by predicting statistics of unknown space. Building on recent work in learning-augmented model based planning under uncertainty, we introduce a high-level state and action abstraction that lets us approximate the challenging Dec-POMDP into a tract… ▽ More We present a novel approach for efficient and reliable goal-directed long-horizon navigation for a multi-robot team in a structured, unknown environment by predicting statistics of unknown space. Building on recent work in learning-augmented model based planning under uncertainty, we introduce a high-level state and action abstraction that lets us approximate the challenging Dec-POMDP into a tractable stochastic MDP. Our Multi-Robot Learning over Subgoals Planner (MR-LSP) guides agents towards coordinated exploration of regions more likely to reach the unseen goal. We demonstrate improvement in cost against other multi-robot strategies; in simulated office-like environments, we show that our approach saves 13.29% (2 robot) and 4.6% (3 robot) average cost versus standard non-learned optimistic planning and a learning-informed baseline. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: 7 pages, 7 figures, ICRA2023

Comparison of Model-Free and Model-Based Learning-Informed Planning for PointGoal Navigation

Authors: Yimeng Li, Arnab Debnath, Gregory J. Stein, Jana Kosecka

Abstract: In recent years several learning approaches to point goal navigation in previously unseen environments have been proposed. They vary in the representations of the environments, problem decomposition, and experimental evaluation. In this work, we compare the state-of-the-art Deep Reinforcement Learning based approaches with Partially Observable Markov Decision Process (POMDP) formulation of the poi… ▽ More In recent years several learning approaches to point goal navigation in previously unseen environments have been proposed. They vary in the representations of the environments, problem decomposition, and experimental evaluation. In this work, we compare the state-of-the-art Deep Reinforcement Learning based approaches with Partially Observable Markov Decision Process (POMDP) formulation of the point goal navigation problem. We adapt the (POMDP) sub-goal framework proposed by [1] and modify the component that estimates frontier properties by using partial semantic maps of indoor scenes built from images' semantic segmentation. In addition to the well-known completeness of the model-based approach, we demonstrate that it is robust and efficient in that it leverages informative, learned properties of the frontiers compared to an optimistic frontier-based planner. We also demonstrate its data efficiency compared to the end-to-end deep reinforcement learning approaches. We compare our results against an optimistic planner, ANS and DD-PPO on Matterport3D dataset using the Habitat Simulator. We show comparable, though slightly worse performance than the SOTA DD-PPO approach, yet with far fewer data. △ Less

Submitted 17 December, 2022; originally announced December 2022.

Comments: arXiv admin note: text overlap with arXiv:2211.07898

Learning-Augmented Model-Based Planning for Visual Exploration

Authors: Yimeng Li, Arnab Debnath, Gregory Stein, Jana Kosecka

Abstract: We consider the problem of time-limited robotic exploration in previously unseen environments where exploration is limited by a predefined amount of time. We propose a novel exploration approach using learning-augmented model-based planning. We generate a set of subgoals associated with frontiers on the current map and derive a Bellman Equation for exploration with these subgoals. Visual sensing a… ▽ More We consider the problem of time-limited robotic exploration in previously unseen environments where exploration is limited by a predefined amount of time. We propose a novel exploration approach using learning-augmented model-based planning. We generate a set of subgoals associated with frontiers on the current map and derive a Bellman Equation for exploration with these subgoals. Visual sensing and advances in semantic mapping of indoor scenes are exploited for training a deep convolutional neural network to estimate properties associated with each frontier: the expected unobserved area beyond the frontier and the expected timesteps (discretized actions) required to explore it. The proposed model-based planner is guaranteed to explore the whole scene if time permits. We thoroughly evaluate our approach on a large-scale pseudo-realistic indoor dataset (Matterport3D) with the Habitat simulator. We compare our approach with classical and more recent RL-based exploration methods. Our approach surpasses the greedy strategies by 2.1% and the RL-based exploration methods by 8.4% in terms of coverage. △ Less

Submitted 9 August, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted to IROS 2023

Rapid Development of a Mobile Robot Simulation Environment

Authors: Gordon Stein, Chan-Jin Chung

Abstract: Robotics simulation provides many advantages during the development of an intelligent ground vehicle (IGV) such as testing the software components in varying scenarios without requiring a complete physical robot. This paper discusses a 3D simulation environment created using rapid application development and the Unity game engine to enable testing during a mobile robotics competition. Our experien… ▽ More Robotics simulation provides many advantages during the development of an intelligent ground vehicle (IGV) such as testing the software components in varying scenarios without requiring a complete physical robot. This paper discusses a 3D simulation environment created using rapid application development and the Unity game engine to enable testing during a mobile robotics competition. Our experience shows that the simulation environment contributed greatly to the development of software for the competition. The simulator also contributed to the hardware development of the robot. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: Presented at AUVSI XPONENTIAL 2016, May 2-5, 2016, Ernest N. Morial Convention Center, New Orleans

A Probabilistic Autoencoder for Type Ia Supernovae Spectral Time Series

Authors: George Stein, Uros Seljak, Vanessa Bohm, G. Aldering, P. Antilogus, C. Aragon, S. Bailey, C. Baltay, S. Bongard, K. Boone, C. Buton, Y. Copin, S. Dixon, D. Fouchez, E. Gangler, R. Gupta, B. Hayden, W. Hillebrandt, M. Karmen, A. G. Kim, M. Kowalski, D. Kusters, P. F. Leget, F. Mondon, J. Nordin , et al. (15 additional authors not shown)

Abstract: We construct a physically-parameterized probabilistic autoencoder (PAE) to learn the intrinsic diversity of type Ia supernovae (SNe Ia) from a sparse set of spectral time series. The PAE is a two-stage generative model, composed of an Auto-Encoder (AE) which is interpreted probabilistically after training using a Normalizing Flow (NF). We demonstrate that the PAE learns a low-dimensional latent sp… ▽ More We construct a physically-parameterized probabilistic autoencoder (PAE) to learn the intrinsic diversity of type Ia supernovae (SNe Ia) from a sparse set of spectral time series. The PAE is a two-stage generative model, composed of an Auto-Encoder (AE) which is interpreted probabilistically after training using a Normalizing Flow (NF). We demonstrate that the PAE learns a low-dimensional latent space that captures the nonlinear range of features that exists within the population, and can accurately model the spectral evolution of SNe Ia across the full range of wavelength and observation times directly from the data. By introducing a correlation penalty term and multi-stage training setup alongside our physically-parameterized network we show that intrinsic and extrinsic modes of variability can be separated during training, removing the need for the additional models to perform magnitude standardization. We then use our PAE in a number of downstream tasks on SNe Ia for increasingly precise cosmological analyses, including automatic detection of SN outliers, the generation of samples consistent with the data distribution, and solving the inverse problem in the presence of noisy and incomplete data to constrain cosmological distance measurements. We find that the optimal number of intrinsic model parameters appears to be three, in line with previous studies, and show that we can standardize our test sample of SNe Ia with an RMS of $0.091 \pm 0.010$ mag, which corresponds to $0.074 \pm 0.010$ mag if peculiar velocity contributions are removed. Trained models and codes are released at \href{https://github.com/georgestein/suPAErnova}{github.com/georgestein/suPAErnova} △ Less

Submitted 15 July, 2022; originally announced July 2022.

Comments: 23 pages, 8 Figures, 1 Table. Accepted to ApJ

Self-supervised similarity search for large scientific datasets

Authors: George Stein, Peter Harrington, Jacqueline Blaum, Tomislav Medan, Zarija Lukic

Abstract: We present the use of self-supervised learning to explore and exploit large unlabeled datasets. Focusing on 42 million galaxy images from the latest data release of the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys, we first train a self-supervised model to distill low-dimensional representations that are robust to symmetries, uncertainties, and noise in each image. We then us… ▽ More We present the use of self-supervised learning to explore and exploit large unlabeled datasets. Focusing on 42 million galaxy images from the latest data release of the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys, we first train a self-supervised model to distill low-dimensional representations that are robust to symmetries, uncertainties, and noise in each image. We then use the representations to construct and publicly release an interactive semantic similarity search tool. We demonstrate how our tool can be used to rapidly discover rare objects given only a single example, increase the speed of crowd-sourcing campaigns, and construct and improve training sets for supervised applications. While we focus on images from sky surveys, the technique is straightforward to apply to any scientific dataset of any dimensionality. The similarity search web app can be found at https://github.com/georgestein/galaxy_search △ Less

Submitted 30 November, 2021; v1 submitted 25 October, 2021; originally announced October 2021.

Comments: 5 pages, 2 figures. The similarity search web app can be found at https://github.com/georgestein/galaxy_search. Accepted to the Fourth Workshop on Machine Learning and the Physical Sciences (NeurIPS 2021). ArXiv admin note: text overlap with arXiv:2110.00023

Mining for Strong Gravitational Lenses with Self-supervised Learning

Authors: George Stein, Jacqueline Blaum, Peter Harrington, Tomislav Medan, Zarija Lukic

Abstract: We employ self-supervised representation learning to distill information from 76 million galaxy images from the Dark Energy Spectroscopic Instrument Legacy Imaging Surveys' Data Release 9. Targeting the identification of new strong gravitational lens candidates, we first create a rapid similarity search tool to discover new strong lenses given only a single labelled example. We then show how train… ▽ More We employ self-supervised representation learning to distill information from 76 million galaxy images from the Dark Energy Spectroscopic Instrument Legacy Imaging Surveys' Data Release 9. Targeting the identification of new strong gravitational lens candidates, we first create a rapid similarity search tool to discover new strong lenses given only a single labelled example. We then show how training a simple linear classifier on the self-supervised representations, requiring only a few minutes on a CPU, can automatically classify strong lenses with great efficiency. We present 1192 new strong lens candidates that we identified through a brief visual identification campaign, and release an interactive web-based similarity search tool and the top network predictions to facilitate crowd-sourcing rapid discovery of additional strong gravitational lenses and other rare objects: https://github.com/georgestein/ssl-legacysurvey. △ Less

Submitted 21 June, 2022; v1 submitted 30 September, 2021; originally announced October 2021.

Comments: 24 Pages, 15 figures, published in ApJ, data at github.com/georgestein/ssl-legacysurvey

Journal ref: The Astrophysical Journal, Volume 932, Number 2, 2022

Superclustering with the Atacama Cosmology Telescope and Dark Energy Survey: I. Evidence for thermal energy anisotropy using oriented stacking

Authors: M. Lokken, R. Hložek, A. van Engelen, M. Madhavacheril, E. Baxter, J. DeRose, C. Doux, S. Pandey, E. S. Rykoff, G. Stein, C. To, T. M. C. Abbott, S. Adhikari, M. Aguena, S. Allam, F. Andrade-Oliveira, J. Annis, N. Battaglia, G. M. Bernstein, E. Bertin, J. R. Bond, D. Brooks, E. Calabrese, A. Carnero Rosell, M. Carrasco Kind , et al. (82 additional authors not shown)

Abstract: The cosmic web contains filamentary structure on a wide range of scales. On the largest scales, superclustering aligns multiple galaxy clusters along inter-cluster bridges, visible through their thermal Sunyaev-Zel'dovich signal in the Cosmic Microwave Background. We demonstrate a new, flexible method to analyze the hot gas signal from multi-scale extended structures. We use a Compton-$y$ map from… ▽ More The cosmic web contains filamentary structure on a wide range of scales. On the largest scales, superclustering aligns multiple galaxy clusters along inter-cluster bridges, visible through their thermal Sunyaev-Zel'dovich signal in the Cosmic Microwave Background. We demonstrate a new, flexible method to analyze the hot gas signal from multi-scale extended structures. We use a Compton-$y$ map from the Atacama Cosmology Telescope (ACT) stacked on redMaPPer cluster positions from the optical Dark Energy Survey (DES). Cutout images from the $y$ map are oriented with large-scale structure information from DES galaxy data such that the superclustering signal is aligned before being overlaid. We find evidence for an extended quadrupole moment of the stacked $y$ signal at the 3.5$σ$ level, demonstrating that the large-scale thermal energy surrounding galaxy clusters is anisotropically distributed. We compare our ACT$\times$DES results with the Buzzard simulations, finding broad agreement. Using simulations, we highlight the promise of this novel technique for constraining the evolution of anisotropic, non-Gaussian structure using future combinations of microwave and optical surveys. △ Less

Submitted 18 July, 2022; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: 37 pages, 23 figures, 4 tables. Added explanatory figure, table, covariance matrix equations, discussion of CIB impact. Matches the version published in ApJ

Learning and Planning for Temporally Extended Tasks in Unknown Environments

Authors: Christopher Bradley, Adam Pacheck, Gregory J. Stein, Sebastian Castro, Hadas Kress-Gazit, Nicholas Roy

Abstract: We propose a novel planning technique for satisfying tasks specified in temporal logic in partially revealed environments. We define high-level actions derived from the environment and the given task itself, and estimate how each action contributes to progress towards completing the task. As the map is revealed, we estimate the cost and probability of success of each action from images and an enco… ▽ More We propose a novel planning technique for satisfying tasks specified in temporal logic in partially revealed environments. We define high-level actions derived from the environment and the given task itself, and estimate how each action contributes to progress towards completing the task. As the map is revealed, we estimate the cost and probability of success of each action from images and an encoding of that action using a trained neural network. These estimates guide search for the minimum-expected-cost plan within our model. Our learned model is structured to generalize across environments and task specifications without requiring retraining. We demonstrate an improvement in total cost in both simulated and real-world experiments compared to a heuristic-driven baseline. △ Less

Submitted 28 April, 2021; v1 submitted 21 April, 2021; originally announced April 2021.

Comments: 7 Pages, 7 Figures, Accepted to ICRA 2021

The LHC Olympics 2020: A Community Challenge for Anomaly Detection in High Energy Physics

Authors: Gregor Kasieczka, Benjamin Nachman, David Shih, Oz Amram, Anders Andreassen, Kees Benkendorfer, Blaz Bortolato, Gustaaf Brooijmans, Florencia Canelli, Jack H. Collins, Biwei Dai, Felipe F. De Freitas, Barry M. Dillon, Ioan-Mihail Dinu, Zhongtian Dong, Julien Donini, Javier Duarte, D. A. Faroughy, Julia Gonski, Philip Harris, Alan Kahn, Jernej F. Kamenik, Charanjit K. Khosa, Patrick Komiske, Luc Le Pottier , et al. (22 additional authors not shown)

Abstract: A new paradigm for data-driven, model-agnostic new physics searches at colliders is emerging, and aims to leverage recent breakthroughs in anomaly detection and machine learning. In order to develop and benchmark new anomaly detection methods within this framework, it is essential to have standard datasets. To this end, we have created the LHC Olympics 2020, a community challenge accompanied by a… ▽ More A new paradigm for data-driven, model-agnostic new physics searches at colliders is emerging, and aims to leverage recent breakthroughs in anomaly detection and machine learning. In order to develop and benchmark new anomaly detection methods within this framework, it is essential to have standard datasets. To this end, we have created the LHC Olympics 2020, a community challenge accompanied by a set of simulated collider events. Participants in these Olympics have developed their methods using an R&D dataset and then tested them on black boxes: datasets with an unknown anomaly (or not). This paper will review the LHC Olympics 2020 challenge, including an overview of the competition, a description of methods deployed in the competition, lessons learned from the experience, and implications for data analyses with future datasets as well as future colliders. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 108 pages, 53 figures, 3 tables

Estimating Galactic Distances From Images Using Self-supervised Representation Learning

Authors: Md Abul Hayat, Peter Harrington, George Stein, Zarija Lukić, Mustafa Mustafa

Abstract: We use a contrastive self-supervised learning framework to estimate distances to galaxies from their photometric images. We incorporate data augmentations from computer vision as well as an application-specific augmentation accounting for galactic dust. We find that the resulting visual representations of galaxy images are semantically useful and allow for fast similarity searches, and can be succ… ▽ More We use a contrastive self-supervised learning framework to estimate distances to galaxies from their photometric images. We incorporate data augmentations from computer vision as well as an application-specific augmentation accounting for galactic dust. We find that the resulting visual representations of galaxy images are semantically useful and allow for fast similarity searches, and can be successfully fine-tuned for the task of redshift estimation. We show that (1) pretraining on a large corpus of unlabeled data followed by fine-tuning on some labels can attain the accuracy of a fully-supervised model which requires 2-4x more labeled data, and (2) that by fine-tuning our self-supervised representations using all available data labels in the Main Galaxy Sample of the Sloan Digital Sky Survey (SDSS), we outperform the state-of-the-art supervised learning method. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Statistical exploration of halo anisotropic clustering and intrinsic alignments with the mass-Peak Patch algorithm

Authors: Bruno Regaldo-Saint Blancard, Sandrine Codis, J. Richard Bond, George Stein

Abstract: The anisotropy or triaxiality of massive dark matter haloes largely defines the structure of the cosmic web, in particular the filaments that join the haloes together. Here we investigate such oriented correlations in mass-Peak Patch halo catalogues by using the initial strain tensor of spherical proto-halo regions to orient the haloes. To go beyond the spherically averaged two-point correlation f… ▽ More The anisotropy or triaxiality of massive dark matter haloes largely defines the structure of the cosmic web, in particular the filaments that join the haloes together. Here we investigate such oriented correlations in mass-Peak Patch halo catalogues by using the initial strain tensor of spherical proto-halo regions to orient the haloes. To go beyond the spherically averaged two-point correlation function of haloes we use oriented stacks to compute oriented two-point correlations: we explicitly break isotropy by imposing a local frame set by the strain tensor of the reference halo before stacking neighbouring haloes. Beyond the exclusion zone of the reference halo, clustering is found to be strongly enhanced along the major direction of the strain tensor as expected. This anisotropic clustering of haloes along filaments is further quantified by using a spherical harmonics decomposition. Furthermore, we compute the evolution of cluster-scale halo principal directions relative to those of their neighbours and show that there are strong correlations extending up to very large scales. In order to provide calculations more suitable to observational confrontations, we also utilize 2D projected versions of some equivalent correlation functions. Finally, we show that the multipole structure of the mass-peak patch halo's anisotropic clustering can be qualitatively captured in an analytic treatment based on peak theory. Though highly informative, analytic evaluation involves extensive use of Monte Carlo methods, which is also what the simulated catalogue uses, taking into account as they do the adaptive nature of the mass-peak patch mass hierarchy and all non-local complexities associated with the exclusion of smaller haloes overlapping with larger ones: there is no substitute for the mass-Peak Patch simulation-based determination of oriented and anisotropic correlations. △ Less

Submitted 30 March, 2021; v1 submitted 5 January, 2021; originally announced January 2021.

Comments: 22 pages, 15 figures, accepted by MNRAS

Self-Supervised Representation Learning for Astronomical Images

Authors: Md Abul Hayat, George Stein, Peter Harrington, Zarija Lukić, Mustafa Mustafa

Abstract: Sky surveys are the largest data generators in astronomy, making automated tools for extracting meaningful scientific information an absolute necessity. We show that, without the need for labels, self-supervised learning recovers representations of sky survey images that are semantically useful for a variety of scientific tasks. These representations can be directly used as features, or fine-tuned… ▽ More Sky surveys are the largest data generators in astronomy, making automated tools for extracting meaningful scientific information an absolute necessity. We show that, without the need for labels, self-supervised learning recovers representations of sky survey images that are semantically useful for a variety of scientific tasks. These representations can be directly used as features, or fine-tuned, to outperform supervised methods trained only on labeled data. We apply a contrastive learning framework on multi-band galaxy photometry from the Sloan Digital Sky Survey (SDSS) to learn image representations. We then use them for galaxy morphology classification, and fine-tune them for photometric redshift estimation, using labels from the Galaxy Zoo 2 dataset and SDSS spectroscopy. In both downstream tasks, using the same learned representations, we outperform the supervised state-of-the-art results, and we show that our approach can achieve the accuracy of supervised models while using 2-4 times fewer labels for training. △ Less

Submitted 8 April, 2021; v1 submitted 23 December, 2020; originally announced December 2020.

Comments: The codes, trained models, and data can be found at https://portal.nersc.gov/project/dasrepo/self-supervised-learning-sdss

Journal ref: The Astrophysical Journal Letters, Volume 911 (2021), Number 2, Letter 33

Unsupervised in-distribution anomaly detection of new physics through conditional density estimation

Authors: George Stein, Uros Seljak, Biwei Dai

Abstract: Anomaly detection is a key application of machine learning, but is generally focused on the detection of outlying samples in the low probability density regions of data. Here we instead present and motivate a method for unsupervised in-distribution anomaly detection using a conditional density estimator, designed to find unique, yet completely unknown, sets of samples residing in high probability… ▽ More Anomaly detection is a key application of machine learning, but is generally focused on the detection of outlying samples in the low probability density regions of data. Here we instead present and motivate a method for unsupervised in-distribution anomaly detection using a conditional density estimator, designed to find unique, yet completely unknown, sets of samples residing in high probability density regions. We apply this method towards the detection of new physics in simulated Large Hadron Collider (LHC) particle collisions as part of the 2020 LHC Olympics blind challenge, and show how we detected a new particle appearing in only 0.08% of 1 million collision events. The results we present are our original blind submission to the 2020 LHC Olympics, where it achieved the state-of-the-art performance. △ Less

Submitted 21 December, 2020; originally announced December 2020.

Comments: Accepted to NeurIPS Machine Learning and the Physical Sciences workshop. See arXiv:2007.00674 for further methods

Stabilizing Transformer-Based Action Sequence Generation For Q-Learning

Authors: Gideon Stein, Andrey Filchenkov, Arip Asadulaev

Abstract: Since the publication of the original Transformer architecture (Vaswani et al. 2017), Transformers revolutionized the field of Natural Language Processing. This, mainly due to their ability to understand timely dependencies better than competing RNN-based architectures. Surprisingly, this architecture change does not affect the field of Reinforcement Learning (RL), even though RNNs are quite popul… ▽ More Since the publication of the original Transformer architecture (Vaswani et al. 2017), Transformers revolutionized the field of Natural Language Processing. This, mainly due to their ability to understand timely dependencies better than competing RNN-based architectures. Surprisingly, this architecture change does not affect the field of Reinforcement Learning (RL), even though RNNs are quite popular in RL, and time dependencies are very common in RL. Recently, Parisotto et al. 2019) conducted the first promising research of Transformers in RL. To support the findings of this work, this paper seeks to provide an additional example of a Transformer-based RL method. Specifically, the goal is a simple Transformer-based Deep Q-Learning method that is stable over several environments. Due to the unstable nature of Transformers and RL, an extensive method search was conducted to arrive at a final method that leverages developments around Transformers as well as Q-learning. The proposed method can match the performance of classic Q-learning on control environments while showing potential on some selected Atari benchmarks. Furthermore, it was critically evaluated to give additional insights into the relation between Transformers and RL. △ Less

Submitted 18 December, 2020; v1 submitted 23 October, 2020; originally announced October 2020.

Comments: Transformers, Reinforcement Learning, 8 pages, AAAI format

Nonlinear 3D Cosmic Web Simulation with Heavy-Tailed Generative Adversarial Networks

Authors: Richard M. Feder, Philippe Berger, George Stein

Abstract: Fast and accurate simulations of the non-linear evolution of the cosmic density field are a major component of many cosmological analyses, but the computational time and storage required to run them can be exceedingly large. For this reason, we use generative adversarial networks (GANs) to learn a compressed representation of the 3D matter density field that is fast and easy to sample, and for the… ▽ More Fast and accurate simulations of the non-linear evolution of the cosmic density field are a major component of many cosmological analyses, but the computational time and storage required to run them can be exceedingly large. For this reason, we use generative adversarial networks (GANs) to learn a compressed representation of the 3D matter density field that is fast and easy to sample, and for the first time show that GANs are capable of generating samples at the level of accuracy of other conventional methods. Using sub-volumes from a suite of GADGET-2 N-body simulations, we demonstrate that a deep-convolutional GAN can generate samples that capture both large- and small-scale features of the matter density field, as validated through a variety of n-point statistics. The use of a data scaling that preserves high-density features and a heavy-tailed latent space prior allow us to obtain state of the art results for fast 3D cosmic web generation. In particular, the mean power spectra from generated samples agree to within 5% up to k=3 and within 10% for k<5 when compared with N-body simulations, and similar accuracy is obtained for a variety of bispectra. By modeling the latent space with a heavy-tailed prior rather than a standard Gaussian, we better capture sample variance in the high-density voxel PDF and reduce errors in power spectrum and bispectrum covariance on all scales. Furthermore, we show that a conditional GAN can smoothly interpolate between samples conditioned on redshift. Deep generative models, such as the ones described in this work, provide great promise as fast, low-memory, high-fidelity forward models of large-scale structure. △ Less

Submitted 13 November, 2020; v1 submitted 6 May, 2020; originally announced May 2020.

Comments: 19 pages, 17 figures. v3: Reflects changes in version published in PRD

Journal ref: Phys. Rev. D 102, 103504 (2020)

Enabling Topological Planning with Monocular Vision

Authors: Gregory J. Stein, Christopher Bradley, Victoria Preston, Nicholas Roy

Abstract: Topological strategies for navigation meaningfully reduce the space of possible actions available to a robot, allowing use of heuristic priors or learning to enable computationally efficient, intelligent planning. The challenges in estimating structure with monocular SLAM in low texture or highly cluttered environments have precluded its use for topological planning in the past. We propose a robus… ▽ More Topological strategies for navigation meaningfully reduce the space of possible actions available to a robot, allowing use of heuristic priors or learning to enable computationally efficient, intelligent planning. The challenges in estimating structure with monocular SLAM in low texture or highly cluttered environments have precluded its use for topological planning in the past. We propose a robust sparse map representation that can be built with monocular vision and overcomes these shortcomings. Using a learned sensor, we estimate high-level structure of an environment from streaming images by detecting sparse vertices (e.g., boundaries of walls) and reasoning about the structure between them. We also estimate the known free space in our map, a necessary feature for planning through previously unknown environments. We show that our mapping technique can be used on real data and is sufficient for planning and exploration in simulated multi-agent search and learned subgoal planning applications. △ Less

Submitted 31 March, 2020; originally announced March 2020.

Comments: 7 pages (6 for content + 1 for references), 5 figures. Accepted to the 2020 IEEE International Conference on Robotics and Automation

The Websky Extragalactic CMB Simulations

Authors: George Stein, Marcelo A. Alvarez, J. Richard Bond, Alexander van Engelen, Nicholas Battaglia

Abstract: We present a new pipeline for the efficient generation of synthetic observations of the extragalactic microwave sky, tailored to large ground-based CMB experiments such as the Simons Observatory, Advanced ACTPol, SPT-3G, and CMB-S4. Such simulated observations are a key technical challenge in cosmology because of the dynamic range and accuracy required. The first part of the pipeline generates a r… ▽ More We present a new pipeline for the efficient generation of synthetic observations of the extragalactic microwave sky, tailored to large ground-based CMB experiments such as the Simons Observatory, Advanced ACTPol, SPT-3G, and CMB-S4. Such simulated observations are a key technical challenge in cosmology because of the dynamic range and accuracy required. The first part of the pipeline generates a random cosmological realization in the form of a dark matter halo catalog and matter displacement field, as seen from a given position. The halo catalog and displacement field are modeled with ellipsoidal collapse dynamics and Lagrangian perturbation theory, respectively. In the second part, the cosmological realization is converted into a set of intensity maps over the range 10 - 10^3 GHz using models based on existing observations and hydrodynamical simulations. These maps include infrared emission from dusty star forming galaxies (CIB), Comptonization of CMB photons by hot gas in groups and clusters through the thermal Sunyaev-Zel'dovich effect (tSZ), Doppler boosting by Thomson scattering of the CMB by bulk flows through the kinetic Sunyaev-Zel'dovich effect (kSZ), and weak gravitational lensing of primary CMB anisotropies by the large-scale distribution of matter in the universe. After describing the pipeline and its implementation, we present the Websky maps, created from a realization of the cosmic web on our past light cone in the redshift interval 0<z<4.6 over the full-sky and a volume of ~(600 Gpc/h)^3 resolved with ~10^12 resolution elements. The Websky maps and halo catalog are publicly available at mocks.cita.utoronto.ca/websky. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: 38 pages, 11 figures, 2 tables. Submitted to JCAP

A narrow-band sodium-resonant fiber-coupled single photon source

Authors: Guilherme Stein, Vladislav Bushmakin, Yijun Wang, Andreas W. Schell, Ilja Gerhardt

Abstract: Quantum technology requires the creation and control over single photons as an important resource. We present a single photon source based on a single molecule which is attached to the end-facet of an optical fiber. To realize a narrow linewidth, the system is cooled down to liquid-helium temperatures. The molecule is optically excited and its fluorescence is collected through the fiber. We have r… ▽ More Quantum technology requires the creation and control over single photons as an important resource. We present a single photon source based on a single molecule which is attached to the end-facet of an optical fiber. To realize a narrow linewidth, the system is cooled down to liquid-helium temperatures. The molecule is optically excited and its fluorescence is collected through the fiber. We have recorded an excitation spectrum, a saturation curve and analyzed the contributions of Raman background fluorescence. This presents to date the crucial limit for the introduced device. The single photon nature is proven by an anti-bunched auto-correlation recording, which also shows coherent Rabi oscillations. △ Less

Submitted 18 September, 2019; originally announced September 2019.

Comments: 8 pages, 7 figures

Journal ref: Phys. Rev. Applied 13, 054042 (2020)

The Simons Observatory: Astro2020 Decadal Project Whitepaper

Authors: The Simons Observatory Collaboration, Maximilian H. Abitbol, Shunsuke Adachi, Peter Ade, James Aguirre, Zeeshan Ahmed, Simone Aiola, Aamir Ali, David Alonso, Marcelo A. Alvarez, Kam Arnold, Peter Ashton, Zachary Atkins, Jason Austermann, Humna Awan, Carlo Baccigalupi, Taylor Baildon, Anton Baleato Lizancos, Darcy Barron, Nick Battaglia, Richard Battye, Eric Baxter, Andrew Bazarko, James A. Beall, Rachel Bean , et al. (258 additional authors not shown)

Abstract: The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021… ▽ More The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4. △ Less

Submitted 16 July, 2019; originally announced July 2019.

Comments: Astro2020 Decadal Project Whitepaper. arXiv admin note: text overlap with arXiv:1808.07445

Journal ref: Bull. Am. Astron. Soc. 51 (2019) 147

Conditioning of Reinforcement Learning Agents and its Policy Regularization Application

Authors: Arip Asadulaev, Igor Kuznetsov, Gideon Stein, Andrey Filchenkov

Abstract: The outcome of Jacobian singular values regularization was studied for supervised learning problems. It also was shown that Jacobian conditioning regularization can help to avoid the ``mode-collapse'' problem in Generative Adversarial Networks. In this paper, we try to answer the following question: Can information about policy conditioning help to shape a more stable and general policy of reinfor… ▽ More The outcome of Jacobian singular values regularization was studied for supervised learning problems. It also was shown that Jacobian conditioning regularization can help to avoid the ``mode-collapse'' problem in Generative Adversarial Networks. In this paper, we try to answer the following question: Can information about policy conditioning help to shape a more stable and general policy of reinforcement learning agents? To answer this question, we conduct a study of Jacobian conditioning behavior during policy optimization. To the best of our knowledge, this is the first work that research condition number in reinforcement learning agents. We propose a conditioning regularization algorithm and test its performance on the range of continuous control tasks. Finally, we compare algorithms on the CoinRun environment with separated train end test levels to analyze how conditioning regularization contributes to agents' generalization. △ Less

Submitted 13 July, 2020; v1 submitted 12 June, 2019; originally announced June 2019.

Deconfusing intensity maps with neural networks

Authors: Daniel N. Pfeffer, Patrick C. Breysse, George Stein

Abstract: Line intensity maps (LIMs) are in principle sensitive to a large amount of information about faint, distant galaxies which are invisible to conventional surveys. However, actually extracting that information from a confused, foreground-contaminated map can be challenging. In this work we present the first application of convolutional neural network (CNN) to directly determine the underlying lumino… ▽ More Line intensity maps (LIMs) are in principle sensitive to a large amount of information about faint, distant galaxies which are invisible to conventional surveys. However, actually extracting that information from a confused, foreground-contaminated map can be challenging. In this work we present the first application of convolutional neural network (CNN) to directly determine the underlying luminosity function of a LIM, including a treatment of extragalactic foregrounds and instrumental noise. We apply the CNN to simulations of mock Carbon Monoxide (CO) line intensity maps similar to those which will be produced by the currently-active COMAP experiment. We evaluate the trained CNN on a number of noise scenarios in order to determine how robust the network predictions are for application to realistic data. We find that, in the ideal case where the mock data capture all of the features of the real data, the CNN performs comparably to or better than conventional analyses. However, the network's accuracy degrades considerably when tested on signals and systematics outside of those it was trained on. For both intensity mapping and cosmology as a whole, this motivates a broad-based study of whether simulated data can ever be generated with sufficient detail to realize the enormous potential of machine learning methods. △ Less

Submitted 24 May, 2019; originally announced May 2019.

Comments: 15 pages, 12 figures, 4 tables

Graphene Induced Large Shift of Surface Plasmon Resonances of Gold Films: Effective Medium Theory for Atomically Thin Materials

Authors: Md Kamrul Alam, Chao Niu, Yanan Wang, Wei Wang, Yang Li, Chong Dai, Tian Tong, Xiaonan Shan, Earl Charlson, Steven Pei, Xiang-Tian Kong, Yandi Hu, Alexey Belyanin, Gila Stein, Zhaoping Liu, Jonathan Hu, Zhiming Wang, Jiming Bao

Abstract: Despite successful modeling of graphene as a 0.34-nm thick optical film synthesized by exfoliation or chemical vapor deposition (CVD), graphene induced shift of surface plasmon resonance (SPR) of gold films has remained controversial. Here we report the resolution of this controversy by developing a clean CVD graphene transfer method and extending Maxwell-Garnet effective medium theory (EMT) to 2D… ▽ More Despite successful modeling of graphene as a 0.34-nm thick optical film synthesized by exfoliation or chemical vapor deposition (CVD), graphene induced shift of surface plasmon resonance (SPR) of gold films has remained controversial. Here we report the resolution of this controversy by developing a clean CVD graphene transfer method and extending Maxwell-Garnet effective medium theory (EMT) to 2D materials. A SPR shift of 0.24 is obtained and it agrees well with 2D EMT in which wrinkled graphene is treated as a 3-nm graphene/air layered composite, in agreement with the average roughness measured by atomic force microscope. Because the anisotropic built-in boundary condition of 2D EMT is compatible with graphene's optical anisotropy, graphene can be modelled as a film thicker than 0.34-nm without changing its optical property; however, its actual roughness, i.e., effective thickness will significantly alter its response to strong out-of-plane fields, leading to a larger SPR shift. △ Less

Submitted 17 April, 2019; originally announced April 2019.

Comments: 18 pages, 5 figures

Journal ref: Phys. Rev. Research 2, 013008 (2020)

Measurement of the Splashback Feature around SZ-selected Galaxy Clusters with DES, SPT and ACT

Authors: T. Shin, S. Adhikari, E. J. Baxter, C. Chang, B. Jain, N. Battaglia, L. Bleem, S. Bocquet, J. DeRose, D. Gruen, M. Hilton, A. Kravtsov, T. McClintock, E. Rozo, E. S. Rykoff, T. N. Varga, R. H. Wechsler, H. Wu, S. Aiola, S. Allam, K. Bechtol, B. A. Benson, E. Bertin, J. R. Bond, M. Brodwin , et al. (85 additional authors not shown)

Abstract: We present a detection of the splashback feature around galaxy clusters selected using their Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, $r_{\rm sp}$, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that $r_{\rm sp}$ inferred from the ob… ▽ More We present a detection of the splashback feature around galaxy clusters selected using their Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, $r_{\rm sp}$, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that $r_{\rm sp}$ inferred from the observed radial distribution of galaxies is affected by selection effects related to the optical cluster-finding algorithms. We test this possibility by measuring the splashback feature in clusters selected via the SZ effect in data from the South Pole Telescope SZ survey and the Atacama Cosmology Telescope Polarimeter survey. The measurement is accomplished by correlating these clusters with galaxies detected in the Dark Energy Survey Year 3 data. The SZ observable used to select clusters in this analysis is expected to have a tighter correlation with halo mass and to be more immune to projection effects and aperture-induced biases than optically selected clusters. We find that the measured $r_{\rm sp}$ for SZ-selected clusters is consistent with the expectations from simulations, although the small number of SZ-selected clusters makes a precise comparison difficult. In agreement with previous work, when using optically selected redMaPPer clusters, $r_{\rm sp}$ is $\sim$ $2σ$ smaller than in the simulations. These results motivate detailed investigations of selection biases in optically selected cluster catalogs and exploration of the splashback feature around larger samples of SZ-selected clusters. Additionally, we investigate trends in the galaxy profile and splashback feature as a function of galaxy color, finding that blue galaxies have profiles close to a power law with no discernible splashback feature, which is consistent with them being on their first infall into the cluster. △ Less

Submitted 24 May, 2019; v1 submitted 14 November, 2018; originally announced November 2018.

Comments: 17 pages, 13 figures, published in MNRAS

The mass-Peak Patch algorithm for fast generation of deep all-sky dark matter halo catalogues and its N-Body validation

Authors: George Stein, Marcelo A. Alvarez, J. Richard Bond

Abstract: We present a detailed description and validation of our massively-parallel update to the mass-Peak Patch method, a fully predictive initial-space algorithm to quickly generate dark matter halo catalogues in very large cosmological volumes. We perform an extensive systematic comparison to a suite of N-body simulations covering a broad range of redshifts and simulation resolutions, and find that, wi… ▽ More We present a detailed description and validation of our massively-parallel update to the mass-Peak Patch method, a fully predictive initial-space algorithm to quickly generate dark matter halo catalogues in very large cosmological volumes. We perform an extensive systematic comparison to a suite of N-body simulations covering a broad range of redshifts and simulation resolutions, and find that, without any parameter fitting, our method is able to generally reproduce N-body results while typically using over 3 orders of magnitude less CPU time, and a fraction of the memory cost. Instead of calculating the full non-linear gravitational collapse determined by an N-body simulation, the mass-Peak Patch method finds an overcomplete set of just-collapsed structures around a hierarchy of density-peak points by coarse-grained (homogeneous) ellipsoidal dynamics. A complete set of mass-peaks, or halos, is then determined by exclusion of overlapping patches, and second order Lagrangian displacements are used to move the halos to their final positions and to give their flow velocities. Our results show that the mass-Peak Patch method is well-suited for creating large ensembles of halo catalogues to mock cosmological surveys, and to aid in complex statistical interpretations of cosmological models. △ Less

Submitted 17 October, 2018; originally announced October 2018.

Comments: 17 pages, 12 figures. Submitted to MNRAS. Comments welcome!

Cross-correlating Carbon Monoxide Line-intensity Maps with Spectroscopic and Photometric Galaxy Surveys

Authors: Dongwoo T. Chung, Marco P. Viero, Sarah E. Church, Risa H. Wechsler, Marcelo A. Alvarez, J. Richard Bond, Patrick C. Breysse, Kieran A. Cleary, Hans K. Eriksen, Marie K. Foss, Joshua O. Gundersen, Stuart E. Harper, Håvard T. Ihle, Laura C. Keating, Norman Murray, Hamsa Padmanabhan, George F. Stein, Ingunn K. Wehus

Abstract: Line-intensity mapping (LIM or IM) is an emerging field of observational work, with strong potential to fit into a larger effort to probe large-scale structure and small-scale astrophysical phenomena using multiple complementary tracers. Taking full advantage of such complementarity means, in part, undertaking line-intensity surveys with galaxy surveys in mind. We consider the potential for detect… ▽ More Line-intensity mapping (LIM or IM) is an emerging field of observational work, with strong potential to fit into a larger effort to probe large-scale structure and small-scale astrophysical phenomena using multiple complementary tracers. Taking full advantage of such complementarity means, in part, undertaking line-intensity surveys with galaxy surveys in mind. We consider the potential for detection of a cross-correlation signal between COMAP and blind surveys based on photometric redshifts (as in COSMOS) or based on spectroscopic data (as with the HETDEX survey of Lyman-$α$ emitters). We find that obtaining $σ_z/(1+z)\lesssim0.003$ accuracy in redshifts and $\gtrsim10^{-4}$ sources per Mpc$^3$ with spectroscopic redshift determination should enable a CO-galaxy cross spectrum detection significance at least twice that of the CO auto spectrum. Either a future targeted spectroscopic survey or a blind survey like HETDEX may be able to meet both of these requirements. △ Less

Submitted 17 January, 2019; v1 submitted 12 September, 2018; originally announced September 2018.

Comments: 19 pages + appendix (31 pages total), 16 figures, 6 tables; accepted for publication in ApJ

Journal ref: ApJ, 872, 186 (2019)

Joint power spectrum and voxel intensity distribution forecast on the CO luminosity function with COMAP

Authors: Håvard Tveit Ihle, Dongwoo Chung, George Stein, Marcelo Alvarez, J. Richard Bond, Patrick C. Breysse, Kieran A. Cleary, Hans Kristian Eriksen, Marie Kristine Foss, Joshua Ott Gundersen, Stuart Harper, Norman Murray, Hamsa Padmanabhan, Marco P. Viero, Ingunn Katerine Wehus

Abstract: We develop a framework for joint constraints on the CO luminosity function based on power spectra (PS) and voxel intensity distributions (VID), and apply this to simulations of COMAP, a CO intensity mapping experiment. This Bayesian framework is based on a Markov chain Monte Carlo (MCMC) sampler coupled to a Gaussian likelihood with a joint PS + VID covariance matrix computed from a large number o… ▽ More We develop a framework for joint constraints on the CO luminosity function based on power spectra (PS) and voxel intensity distributions (VID), and apply this to simulations of COMAP, a CO intensity mapping experiment. This Bayesian framework is based on a Markov chain Monte Carlo (MCMC) sampler coupled to a Gaussian likelihood with a joint PS + VID covariance matrix computed from a large number of fiducial simulations, and re-calibrated with a small number of simulations per MCMC step. The simulations are based on dark matter halos from fast peak patch simulations combined with the $L_\text{CO}(M_\text{halo})$ model of Li et al. (2016). We find that the relative power to constrain the CO luminosity function depends on the luminosity range of interest. In particular, the VID is more sensitive at both small and large luminosities, while the PS is more sensitive at intermediate luminosities. The joint analysis is superior to using either observable separately. When averaging over CO luminosities ranging between $L_\text{CO} = 10^4-10^7L_\odot$, and over 10 cosmological realizations of COMAP Phase 2, the uncertainties (in dex) are larger by 58 % and 30 % for the PS and VID, respectively, when compared to the joint analysis (PS + VID). This method is generally applicable to any other random field, with a complicated likelihood, as long a fast simulation procedure is available. △ Less

Submitted 20 March, 2019; v1 submitted 22 August, 2018; originally announced August 2018.

Comments: 13 pages, 5 figures. As accepted to ApJ

Journal ref: ApJ, 871, 1 (2019)

The Simons Observatory: Science goals and forecasts

Authors: The Simons Observatory Collaboration, Peter Ade, James Aguirre, Zeeshan Ahmed, Simone Aiola, Aamir Ali, David Alonso, Marcelo A. Alvarez, Kam Arnold, Peter Ashton, Jason Austermann, Humna Awan, Carlo Baccigalupi, Taylor Baildon, Darcy Barron, Nick Battaglia, Richard Battye, Eric Baxter, Andrew Bazarko, James A. Beall, Rachel Bean, Dominic Beck, Shawn Beckman, Benjamin Beringue, Federico Bianchini , et al. (225 additional authors not shown)

Abstract: The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225… ▽ More The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $μ$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $σ(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $μ$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources. △ Less

Submitted 1 March, 2019; v1 submitted 22 August, 2018; originally announced August 2018.

Comments: This paper presents an overview of the Simons Observatory science goals, details about the instrument will be presented in a companion paper. The author contribution to this paper is available at https://simonsobservatory.org/publications.php (Abstract abridged) -- matching version published in JCAP

Journal ref: JCAP 1902 (2019) 056

CCAT-prime: Science with an Ultra-widefield Submillimeter Observatory at Cerro Chajnantor

Authors: G. J. Stacey, M. Aravena, K. Basu, N. Battaglia, B. Beringue, F. Bertoldi, J. R. Bond, P. Breysse, R. Bustos, S. Chapman, D. T. Chung, N. Cothard, J. Erler, M. Fich, S. Foreman, P. Gallardo, R. Giovanelli, U. U. Graf, M. P. Haynes, R. Herrera-Camus, T. L. Herter, R. Hložek, D. Johnstone, L. Keating, B. Magnelli , et al. (15 additional authors not shown)

Abstract: We present the detailed science case, and brief descriptions of the telescope design, site, and first light instrument plans for a new ultra-wide field submillimeter observatory, CCAT-prime, that we are constructing at a 5600 m elevation site on Cerro Chajnantor in northern Chile. Our science goals are to study star and galaxy formation from the epoch of reionization to the present, investigate th… ▽ More We present the detailed science case, and brief descriptions of the telescope design, site, and first light instrument plans for a new ultra-wide field submillimeter observatory, CCAT-prime, that we are constructing at a 5600 m elevation site on Cerro Chajnantor in northern Chile. Our science goals are to study star and galaxy formation from the epoch of reionization to the present, investigate the growth of structure in the Universe, improve the precision of B-mode CMB measurements, and investigate the interstellar medium and star formation in the Galaxy and nearby galaxies through spectroscopic, polarimetric, and broadband surveys at wavelengths from 200 um to 2 mm. These goals are realized with our two first light instruments, a large field-of-view (FoV) bolometer-based imager called Prime-Cam (that has both camera and an imaging spectrometer modules), and a multi-beam submillimeter heterodyne spectrometer, CHAI. CCAT-prime will have very high surface accuracy and very low system emissivity, so that combined with its wide FoV at the unsurpassed CCAT site our telescope/instrumentation combination is ideally suited to pursue this science. The CCAT-prime telescope is being designed and built by Vertex Antennentechnik GmbH. We expect to achieve first light in the spring of 2021. △ Less

Submitted 11 July, 2018; originally announced July 2018.

Comments: Presented at SPIE Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy IX, June 14th, 2018

Comparing approximate methods for mock catalogues and covariance matrices III: Bispectrum

Authors: Manuel Colavincenzo, Emiliano Sefusatti, Pierluigi Monaco, Linda Blot, Martin Crocce, Martha Lippich, Ariel G. Sánchez, Marcelo A. Alvarez, Aniket Agrawal, Santiago Avila, Andrés Balaguera-Antolínez, Richard Bond, Sandrine Codis, Claudio Dalla Vecchia, Antonio Dorta, Pablo Fosalba, Albert Izard, Francisco-Shu Kitaura, Marcos Pellejero-Ibanez, George Stein, Mohammadjavad Vakili, Gustavo Yepes

Abstract: We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in… ▽ More We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method. △ Less

Submitted 8 October, 2018; v1 submitted 25 June, 2018; originally announced June 2018.

Comments: Additional results with respect to v1, new section and new figures added. 25 pages, 1 table 18 figures

Comparing approximate methods for mock catalogues and covariance matrices II: Power spectrum multipoles

Authors: Linda Blot, Martin Crocce, Emiliano Sefusatti, Martha Lippich, Ariel G. Sánchez, Manuel Colavincenzo, Pierluigi Monaco, Marcelo A. Alvarez, Aniket Agrawal, Santiago Avila, Andrés Balaguera-Antolínez, Richard Bond, Sandrine Codis, Claudio Dalla Vecchia, Antonio Dorta, Pablo Fosalba, Albert Izard, Francisco-Shu Kitaura, Marcos Pellejero-Ibanez, George Stein, Mohammadjavad Vakili, Gustavo Yepes

Abstract: We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo densi… ▽ More We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance. △ Less

Submitted 18 February, 2019; v1 submitted 25 June, 2018; originally announced June 2018.

Comments: 20 pages, 16 figures, replaced to match accepted MNRAS version. Results on parameter errors changed

Comparing approximate methods for mock catalogues and covariance matrices I: correlation function

Authors: Martha Lippich, Ariel G. Sánchez, Manuel Colavincenzo, Emiliano Sefusatti, Pierluigi Monaco, Linda Blot, Martin Crocce, Marcelo A. Alvarez, Aniket Agrawal, Santiago Avila, Andrés Balaguera-Antolínez, Richard Bond, Sandrine Codis, Claudio Dalla Vecchia, Antonio Dorta, Pablo Fosalba, Albert Izard, Francisco-Shu Kitaura, Marcos Pellejero-Ibanez, George Stein, Mohammadjavad Vakili, Gustavo Yepes

Abstract: This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that fol… ▽ More This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others. △ Less

Submitted 13 May, 2019; v1 submitted 25 June, 2018; originally announced June 2018.

Comments: 23 pages, 11 figures. Replaced to match accepted MNRAS version. Included Kullback-Leibler divergence

A volumetric deep Convolutional Neural Network for simulation of mock dark matter halo catalogues

Authors: Philippe Berger, George Stein

Abstract: For modern large-scale structure survey techniques it has become standard practice to test data analysis pipelines on large suites of mock simulations, a task which is currently prohibitively expensive for full N-body simulations. Instead of calculating this costly gravitational evolution, we have trained a three-dimensional deep Convolutional Neural Network (CNN) to identify dark matter protohalo… ▽ More For modern large-scale structure survey techniques it has become standard practice to test data analysis pipelines on large suites of mock simulations, a task which is currently prohibitively expensive for full N-body simulations. Instead of calculating this costly gravitational evolution, we have trained a three-dimensional deep Convolutional Neural Network (CNN) to identify dark matter protohalos directly from the cosmological initial conditions. Training on halo catalogues from the Peak Patch semi-analytic code, we test various CNN architectures and find they generically achieve a Dice coefficient of ~92% in only 24 hours of training. We present a simple and fast geometric halo finding algorithm to extract halos from this powerful pixel-wise binary classifier and find that the predicted catalogues match the mass function and power spectra of the ground truth simulations to within ~10%. We investigate the effect of long-range tidal forces on an object-by-object basis and find that the network's predictions are consistent with the non-linear ellipsoidal collapse equations used explicitly by the Peak Patch algorithm. △ Less

Submitted 19 November, 2018; v1 submitted 11 May, 2018; originally announced May 2018.

Comments: 12 pages, 8 figures, 1 table. Accepted to MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 482, Issue 3, p.2861-2871, 2019

Weak-Lensing Mass Calibration of ACTPol Sunyaev-Zel'dovich Clusters with the Hyper Suprime-Cam Survey

Authors: Hironao Miyatake, Nicholas Battaglia, Matt Hilton, Elinor Medezinski, Atsushi J. Nishizawa, Surhud More, Simone Aiola, Neta Bahcall, J. Richard Bond, Erminia Calabrese, Steve K. Choi, Mark J. Devlin, Joanna Dunkley, Rolando Dunner, Brittany Fuzia, Patricio Gallardo, Megan Gralla, Matthew Hasselfield, Mark Halpern, Chiaki Hikage, J. Colin Hill, Adam D. Hincks, Renée Hložek, Kevin Huffenberger, John P. Hughes , et al. (35 additional authors not shown)

Abstract: We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking brig… ▽ More We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements. △ Less

Submitted 16 April, 2018; originally announced April 2018.

Comments: 19 pages, 11 figures, 2 tables, comments are welcome