Showing 1–50 of 85 results for author: Lazar, A

Decorated square paths at q=-1

Authors: Sylvie Corteel, Alexander Lazar, Anna Vanden Wyngaerd

Abstract: The valley Delta square conjecture states that the symmetric function $\frac{[n-k]_q}{[n]_q}Δ_{e_{n-k}}ω(p_n)$ can be expressed as the enumerator of a certain class of decorated square paths with respect to the bistatistic (dinv,area). Inspired by recent positivity results of Corteel, Josuat-Vergès, and Vanden Wyngaerd, we study the evaluation of this enumerator at $q=-1$. By considering a cyclic… ▽ More The valley Delta square conjecture states that the symmetric function $\frac{[n-k]_q}{[n]_q}Δ_{e_{n-k}}ω(p_n)$ can be expressed as the enumerator of a certain class of decorated square paths with respect to the bistatistic (dinv,area). Inspired by recent positivity results of Corteel, Josuat-Vergès, and Vanden Wyngaerd, we study the evaluation of this enumerator at $q=-1$. By considering a cyclic group action on the decorated square paths which we call cutting and pasting, we show that $\left.\left\langle \frac{[n-k]_q}{[n]_q}Δ_{e_{n-k}}ω(p_n), h_1^n\right\rangle\right|_{q=-1}$ is $0$ whenever $n-k$ is even, and is a positive polynomial related to the Euler numbers when $n-k$ is odd. We also show that the combinatorics of this enumerator is closely connected to that of the Dyck path enumerator for $\langleΔ_{e_{n-k-1}}'e_n,h_1^n\rangle$ considered by Corteel-Josuat Vergès-Vanden Wyngaerd. △ Less

Submitted 20 August, 2024; originally announced August 2024.

Elevated UV luminosity density at Cosmic Dawn explained by non-evolving, weakly-mass dependent star formation efficiency

Authors: Robert Feldmann, Michael Boylan-Kolchin, James S. Bullock, Onur Çatmabacak, Claude-André Faucher-Giguère, Christopher C. Hayward, Dušan Kereš, Alexandres Lazar, Lichen Liang, Jorge Moreno, Pascal A. Oesch, Eliot Quataert, Xuejian Shen, Guochao Sun

Abstract: Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR… ▽ More Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox-HR re-simulates the cosmic volume (L = 22.1 cMpc) of the original FIREbox run with eight times higher mass resolution (m_b ~ 7800 M_sun), but with identical physics, down to z ~ 6. FIREbox-HR predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at z ~ 6 - 14, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox-HR, the SFE - halo mass relation for intermediate mass halos (M_halo ~ 10^9 - 10^11 M_sun) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE - halo mass relation lead to a larger contribution from lower mass halos at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE - halo mass relation inferred from FIREbox-HR allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at z > 12 will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 26 pages, 14 figures, 5 tables, submitted to MNRAS, comments welcome

Voronoi topology in two dimensions: theory, algorithms, and applications

Authors: Emanuel A. Lazar, Jiayin Lu, Chris H. Rycroft, Deborah Schwartz

Abstract: This paper introduces a new approach toward characterizing local structural features of two-dimensional particle systems. The approach can accurately identify and characterize defects in high-temperature crystals, distinguish a wide range of nominally disordered systems, and robustly describe complex structures such as grain boundaries. This paper also introduces two-dimensional functionality into… ▽ More This paper introduces a new approach toward characterizing local structural features of two-dimensional particle systems. The approach can accurately identify and characterize defects in high-temperature crystals, distinguish a wide range of nominally disordered systems, and robustly describe complex structures such as grain boundaries. This paper also introduces two-dimensional functionality into the open-source software program VoroTop which automates this analysis. This software package is built on a recently-introduced multithreaded version of VORO++, enabling the analysis of systems with billions of particles on high-performance computer architectures. △ Less

Submitted 1 June, 2024; originally announced June 2024.

Comments: 22 pages, 9 figures

A Language Model for Particle Tracking

Authors: Andris Huang, Yash Melkani, Paolo Calafiura, Alina Lazar, Daniel Thomas Murnane, Minh-Tuan Pham, Xiangyang Ju

Abstract: Particle tracking is crucial for almost all physics analysis programs at the Large Hadron Collider. Deep learning models are pervasively used in particle tracking related tasks. However, the current practice is to design and train one deep learning model for one task with supervised learning techniques. The trained models work well for tasks they are trained on but show no or little generalization… ▽ More Particle tracking is crucial for almost all physics analysis programs at the Large Hadron Collider. Deep learning models are pervasively used in particle tracking related tasks. However, the current practice is to design and train one deep learning model for one task with supervised learning techniques. The trained models work well for tasks they are trained on but show no or little generalization capabilities. We propose to unify these models with a language model. In this paper, we present a tokenized detector representation that allows us to train a BERT model for particle tracking. The trained BERT model, namely TrackingBERT, offers latent detector module embedding that can be used for other tasks. This work represents the first step towards developing a foundational model for particle detector understanding. △ Less

Submitted 14 February, 2024; originally announced February 2024.

Comments: 7 pages, 3 figures, A Proceeding of the Connecting the Dots Workshop (CTD 2023)

Report number: PROC-CTD2023-33

Shuffle theorems and sandpiles

Authors: Michele D'Adderio, Mark Dukes, Alessandro Iraci, Alexander Lazar, Yvan Le Borgne, Anna Vanden Wyngaerd

Abstract: We provide an explicit description of the recurrent configurations of the sandpile model on a family of graphs $\widehat{G}_{μ,ν}$, which we call clique-independent graphs, indexed by two compositions $μ$ and $ν$. Moreover, we define a delay statistic on these configurations, and we show that, together with the usual level statistic, it can be used to provide a new combinatorial interpretation of… ▽ More We provide an explicit description of the recurrent configurations of the sandpile model on a family of graphs $\widehat{G}_{μ,ν}$, which we call clique-independent graphs, indexed by two compositions $μ$ and $ν$. Moreover, we define a delay statistic on these configurations, and we show that, together with the usual level statistic, it can be used to provide a new combinatorial interpretation of the celebrated shuffle theorem of Carlsson and Mellit. More precisely, we will see how to interpret the polynomials $\langle \nabla e_n, e_μh_ν\rangle$ in terms of these configurations. △ Less

Submitted 15 January, 2024; v1 submitted 12 January, 2024; originally announced January 2024.

Comments: 12 pages, 2 figures. Comments are welcome!

MSC Class: 05E05; 82B20

An analytic surface density profile for $Λ$CDM halos and gravitational lensing studies

Authors: Alexandres Lazar, James S. Bullock, Anna Nierenberg, Leonidas Moustakas, Michael Boylan-Kolchin

Abstract: We introduce an analytic surface density profile for dark matter halos that accurately reproduces the structure of simulated halos of mass $M_{\rm vir} = 10^{7-11}\ M_\odot$, making it useful for modeling line-of-sight perturbers in strong gravitational lensing models. The two-parameter function has an analytic deflection potential and is more accurate than the projected Navarro, Frenk & White (NF… ▽ More We introduce an analytic surface density profile for dark matter halos that accurately reproduces the structure of simulated halos of mass $M_{\rm vir} = 10^{7-11}\ M_\odot$, making it useful for modeling line-of-sight perturbers in strong gravitational lensing models. The two-parameter function has an analytic deflection potential and is more accurate than the projected Navarro, Frenk & White (NFW) profile commonly adopted at this mass scale for perturbers, especially at the small radii of most relevant for lensing perturbations. Using a characteristic radius, $R_{-1}$, where the log slope of surface density is equal to $-1$, and an associated surface density, $Σ_{-1}$, we can represent the expected lensing signal from line-of-sight halos statistically, for an ensemble of halo orientations, using a distribution of {\em projected concentration} parameters, $\mathcal{C}_{\rm vir} := r_{\rm vir}/ R_{-1}$. Though an individual halo can have a projected concentration that varies with orientation with respect to the observer, the range of projected concentrations correlates with the usual three-dimensional halo concentration in a way that enables ease of use. △ Less

Submitted 3 January, 2024; v1 submitted 21 April, 2023; originally announced April 2023.

Comments: 15 pages, 13 figures, Accepted to MNRAS

Identification of high-reliability regions of machine learning predictions in materials science using transparent conducting oxides and perovskites as examples

Authors: Evan M. Askanazi, Emanuel A. Lazar, Ilya Grinberg

Abstract: Progress in the application of machine learning (ML) methods to materials design is hindered by the lack of understanding of the reliability of ML predictions, in particular for the application of ML to small data sets often found in materials science. Using ML prediction for transparent conductor oxide formation energy and band gap, dilute solute diffusion, and perovskite formation energy, band g… ▽ More Progress in the application of machine learning (ML) methods to materials design is hindered by the lack of understanding of the reliability of ML predictions, in particular for the application of ML to small data sets often found in materials science. Using ML prediction for transparent conductor oxide formation energy and band gap, dilute solute diffusion, and perovskite formation energy, band gap and lattice parameter as examples, we demonstrate that 1) analysis of ML results by construction of a convex hull in feature space that encloses accurately predicted systems can be used to identify regions in feature space for which ML predictions are highly reliable 2) analysis of the systems enclosed by the convex hull can be used to extract physical understanding and 3) materials that satisfy all well-known chemical and physical principles that make a material physically reasonable are likely to be similar and show strong relationships between the properties of interest and the standard features used in ML. We also show that similar to the composition-structure-property relationships, inclusion in the ML training data set of materials from classes with different chemical properties will not be beneficial and will slightly decrease the accuracy of ML prediction and that reliable results likely will be obtained by ML model for narrow classes of similar materials even in the case where the ML model will show large errors on the dataset consisting of several classes of materials. Our work suggests that analysis of the error distributions of ML predictions will be beneficial for the further development of the application of ML methods in material science. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Understanding Context to Capture when Reconstructing Meaningful Spaces for Remote Instruction and Connecting in XR

Authors: Hanuma Teja Maddali, Amanda Lazar

Abstract: Recent technological advances are enabling HCI researchers to explore interaction possibilities for remote XR collaboration using high-fidelity reconstructions of physical activity spaces. However, creating these reconstructions often lacks user involvement with an overt focus on capturing sensory context that does not necessarily augment an informal social experience. This work seeks to understan… ▽ More Recent technological advances are enabling HCI researchers to explore interaction possibilities for remote XR collaboration using high-fidelity reconstructions of physical activity spaces. However, creating these reconstructions often lacks user involvement with an overt focus on capturing sensory context that does not necessarily augment an informal social experience. This work seeks to understand social context that can be important for reconstruction to enable XR applications for informal instructional scenarios. Our study involved the evaluation of an XR remote guidance prototype by 8 intergenerational groups of closely related gardeners using reconstructions of personally meaningful spaces in their gardens. Our findings contextualize physical objects and areas with various motivations related to gardening and detail perceptions of XR that might affect the use of reconstructions for remote interaction. We discuss implications for user involvement to create reconstructions that better translate real-world experience, encourage reflection, incorporate privacy considerations, and preserve shared experiences with XR as a medium for informal intergenerational activities. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 26 pages, 5 figures, 4 tables

ACM Class: H.5.2; I.2.1; I.3.5; I.3.7; I.3.8; K.8.0; K.3.1; K.4.2

Sociality and Skill Sharing in the Garden

Authors: Hanuma Teja Maddali, Amanda Lazar

Abstract: Gardening is an activity that involves a number of dimensions of increasing interest to HCI and CSCW researchers, including recreation, sustainability, and engagement with nature. This paper considers the garden setting in order to understand the role that collaborative and social computing technologies might play for practitioners engaging in outdoor skilled activities. We conducted participant o… ▽ More Gardening is an activity that involves a number of dimensions of increasing interest to HCI and CSCW researchers, including recreation, sustainability, and engagement with nature. This paper considers the garden setting in order to understand the role that collaborative and social computing technologies might play for practitioners engaging in outdoor skilled activities. We conducted participant observations with nine experienced gardeners aged 22-71 years. Through this process, we find that gardeners continuously configure their environments to accommodate their preferences for sociality. They share embodied skills and help others attune to sensory information in person, but also influence learning through the features in their garden that are observed by others. This paper provides an understanding of sociality in the garden, highlights skill sharing as a key domain for design in this space, and contributes design considerations for collaborative technologies in outdoor settings. △ Less

Submitted 19 November, 2022; originally announced November 2022.

Comments: 13 pages

ACM Class: H.5.0; K.3.1; K.4.2

Journal ref: Proceedings of the 2020 CHI Conference on Human Factors in Computing Systems (CHI '20). Association for Computing Machinery, New York, NY, USA, 1-13

Investigating the Potential of Artificial Intelligence Powered Interfaces to Support Different Types of Memory for People with Dementia

Authors: Hanuma Teja Maddali, Emma Dixon, Alisha Pradhan, Amanda Lazar

Abstract: There has been a growing interest in HCI to understand the specific technological needs of people with dementia and supporting them in self-managing daily activities. One of the most difficult challenges to address is supporting the fluctuating accessibility needs of people with dementia, which vary with the specific type of dementia and the progression of the condition. Researchers have identifie… ▽ More There has been a growing interest in HCI to understand the specific technological needs of people with dementia and supporting them in self-managing daily activities. One of the most difficult challenges to address is supporting the fluctuating accessibility needs of people with dementia, which vary with the specific type of dementia and the progression of the condition. Researchers have identified auto-personalized interfaces, and more recently, Artificial Intelligence or AI-driven personalization as a potential solution to making commercial technology accessible in a scalable manner for users with fluctuating ability. However, there is a lack of understanding on the perceptions of people with dementia around AI as an aid to their everyday technology use and its role in their overall self-management systems, which include other non-AI technology, and human assistance. In this paper, we present future directions for the design of AI-based systems to personalize an interface for dementia-related changes in different types of memory, along with expectations for AI interactions with the user with dementia. △ Less

Submitted 19 November, 2022; originally announced November 2022.

Comments: 7 pages

ACM Class: H.5.2; I.2.1; K.8.0

Journal ref: Extended Abstracts of the 2022 CHI Conference on Human Factors in Computing Systems (CHI EA '22). Association for Computing Machinery, New York, NY, USA, Article 226, 1-7

Pair correlation function based on Voronoi topology

Authors: Vasco M. Worlitzer, Gil Ariel, Emanuel A. Lazar

Abstract: The pair correlation function (PCF) has proven an effective tool for analyzing many physical systems due to its simplicity and its applicability to simulated and experimental data. However, as an averaged quantity, the PCF can fail to capture subtle structural differences in particle arrangements, even when those differences can have a major impact on system properties. Here, we use Voronoi topolo… ▽ More The pair correlation function (PCF) has proven an effective tool for analyzing many physical systems due to its simplicity and its applicability to simulated and experimental data. However, as an averaged quantity, the PCF can fail to capture subtle structural differences in particle arrangements, even when those differences can have a major impact on system properties. Here, we use Voronoi topology to introduce a discrete version of the PCF that highlights local inter-particle topological configurations. The advantages of the Voronoi PCF are demonstrated in several examples including crystalline, hyperuniform, and active systems showing clustering and giant number fluctuations. △ Less

Submitted 13 December, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 8 pages, 9 figures

Journal ref: Phys. Rev. E 108, 064115 (2023)

An extension to VORO++ for multithreaded computation of Voronoi cells

Authors: Jiayin Lu, Emanuel A. Lazar, Chris H. Rycroft

Abstract: VORO++ is a software library written in C++ for computing the Voronoi tessellation, a technique in computational geometry that is widely used for analyzing systems of particles. VORO++ was released in 2009 and is based on computing the Voronoi cell for each particle individually. Here, we take advantage of modern computer hardware, and extend the original serial version to allow for multithreaded… ▽ More VORO++ is a software library written in C++ for computing the Voronoi tessellation, a technique in computational geometry that is widely used for analyzing systems of particles. VORO++ was released in 2009 and is based on computing the Voronoi cell for each particle individually. Here, we take advantage of modern computer hardware, and extend the original serial version to allow for multithreaded computation of Voronoi cells via the OpenMP application programming interface. We test the performance of the code, and demonstrate that we can achieve parallel efficiencies greater than 95% in many cases. The multithreaded extension follows standard OpenMP programming paradigms, allowing it to be incorporated into other programs. We provide an example of this using the VoroTop software library, performing a multithreaded Voronoi cell topology analysis of up to 102.4 million particles. △ Less

Submitted 8 July, 2023; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: Fix typo and section numbers

What Makes You Hold on to That Old Car? Joint Insights from Machine Learning and Multinomial Logit on Vehicle-level Transaction Decisions

Authors: Ling Jin, Alina Lazar, Caitlin Brown, Bingrong Sun, Venu Garikapati, Srinath Ravulaparthy, Qianmiao Chen, Alexander Sim, Kesheng Wu, Tin Ho, Thomas Wenzel, C. Anna Spurlock

Abstract: What makes you hold on that old car? While the vast majority of the household vehicles are still powered by conventional internal combustion engines, the progress of adopting emerging vehicle technologies will critically depend on how soon the existing vehicles are transacted out of the household fleet. Leveraging a nationally representative longitudinal data set, the Panel Study of Income Dynamic… ▽ More What makes you hold on that old car? While the vast majority of the household vehicles are still powered by conventional internal combustion engines, the progress of adopting emerging vehicle technologies will critically depend on how soon the existing vehicles are transacted out of the household fleet. Leveraging a nationally representative longitudinal data set, the Panel Study of Income Dynamics, this study examines how household decisions to dispose of or replace a given vehicle are: (1) influenced by the vehicle's attributes, (2) mediated by households' concurrent socio-demographic and economic attributes, and (3) triggered by key life cycle events. Coupled with a newly developed machine learning interpretation tool, TreeExplainer, we demonstrate an innovative use of machine learning models to augment traditional logit modeling to both generate behavioral insights and improve model performance. We find the two gradient-boosting-based methods, CatBoost and LightGBM, are the best performing machine learning models for this problem. The multinomial logistic model can achieve similar performance levels after its model specification is informed by TreeExplainer. Both machine learning and multinomial logit models suggest that while older vehicles are more likely to be disposed of or replaced than newer ones, such probability decreases as the vehicles serve the family longer. We find that married families, families with higher education levels, homeowners, and older families tend to keep their vehicles longer. Life events such as childbirth, residential relocation, and change of household composition and income are found to increase vehicle disposal and/or replacement. We provide additional insights on the timing of vehicle replacement or disposal, in particular, the presence of children and childbirth events are more strongly associated with vehicle replacement among younger parents. △ Less

Submitted 13 May, 2022; originally announced May 2022.

MyMove: Facilitating Older Adults to Collect In-Situ Activity Labels on a Smartwatch with Speech

Authors: Young-Ho Kim, Diana Chou, Bongshin Lee, Margaret Danilovich, Amanda Lazar, David E. Conroy, Hernisa Kacorri, Eun Kyoung Choe

Abstract: Current activity tracking technologies are largely trained on younger adults' data, which can lead to solutions that are not well-suited for older adults. To build activity trackers for older adults, it is crucial to collect training data with them. To this end, we examine the feasibility and challenges with older adults in collecting activity labels by leveraging speech. Specifically, we built My… ▽ More Current activity tracking technologies are largely trained on younger adults' data, which can lead to solutions that are not well-suited for older adults. To build activity trackers for older adults, it is crucial to collect training data with them. To this end, we examine the feasibility and challenges with older adults in collecting activity labels by leveraging speech. Specifically, we built MyMove, a speech-based smartwatch app to facilitate the in-situ labeling with a low capture burden. We conducted a 7-day deployment study, where 13 older adults collected their activity labels and smartwatch sensor data, while wearing a thigh-worn activity monitor. Participants were highly engaged, capturing 1,224 verbal reports in total. We extracted 1,885 activities with corresponding effort level and timespan, and examined the usefulness of these reports as activity labels. We discuss the implications of our approach and the collected dataset in supporting older adults through personalized activity tracking technologies. △ Less

Submitted 31 March, 2022; originally announced April 2022.

Comments: To appear at ACM CHI 2022. 21 pages, 3 figures, 7 tables. For the NSF funded project, visit https://mymove-collective.github.io

ACM Class: H.5.2; H.5.1; I.2.1

Reconstruction of Large Radius Tracks with the Exa.TrkX pipeline

Authors: Chun-Yi Wang, Xiangyang Ju, Shih-Chieh Hsu, Daniel Murnane, Paolo Calafiura, Steven Farrell, Maria Spiropulu, Jean-Roch Vlimant, Adam Aurisano, V Hewes, Giuseppe Cerati, Lindsey Gray, Thomas Klijnsma, Jim Kowalkowski, Markus Atkinson, Mark Neubauer, Gage DeZoort, Savannah Thais, Alexandra Ballow, Alina Lazar, Sylvain Caillou, Charline Rougier, Jan Stark, Alexis Vallier, Jad Sardain

Abstract: Particle tracking is a challenging pattern recognition task at the Large Hadron Collider (LHC) and the High Luminosity-LHC. Conventional algorithms, such as those based on the Kalman Filter, achieve excellent performance in reconstructing the prompt tracks from the collision points. However, they require dedicated configuration and additional computing time to efficiently reconstruct the large rad… ▽ More Particle tracking is a challenging pattern recognition task at the Large Hadron Collider (LHC) and the High Luminosity-LHC. Conventional algorithms, such as those based on the Kalman Filter, achieve excellent performance in reconstructing the prompt tracks from the collision points. However, they require dedicated configuration and additional computing time to efficiently reconstruct the large radius tracks created away from the collision points. We developed an end-to-end machine learning-based track finding algorithm for the HL-LHC, the Exa.TrkX pipeline. The pipeline is designed so as to be agnostic about global track positions. In this work, we study the performance of the Exa.TrkX pipeline for finding large radius tracks. Trained with all tracks in the event, the pipeline simultaneously reconstructs prompt tracks and large radius tracks with high efficiencies. This new capability offered by the Exa.TrkX pipeline may enable us to search for new physics in real time. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 5 pages, 3 figures. Proceedings of 20th International Workshop on Advanced Computing and Analysis Techniques in Physics Research

Accelerating the Inference of the Exa.TrkX Pipeline

Authors: Alina Lazar, Xiangyang Ju, Daniel Murnane, Paolo Calafiura, Steven Farrell, Yaoyuan Xu, Maria Spiropulu, Jean-Roch Vlimant, Giuseppe Cerati, Lindsey Gray, Thomas Klijnsma, Jim Kowalkowski, Markus Atkinson, Mark Neubauer, Gage DeZoort, Savannah Thais, Shih-Chieh Hsu, Adam Aurisano, V Hewes, Alexandra Ballow, Nirajan Acharya, Chun-yi Wang, Emma Liu, Alberto Lucas

Abstract: Recently, graph neural networks (GNNs) have been successfully used for a variety of particle reconstruction problems in high energy physics, including particle tracking. The Exa.TrkX pipeline based on GNNs demonstrated promising performance in reconstructing particle tracks in dense environments. It includes five discrete steps: data encoding, graph building, edge filtering, GNN, and track labelin… ▽ More Recently, graph neural networks (GNNs) have been successfully used for a variety of particle reconstruction problems in high energy physics, including particle tracking. The Exa.TrkX pipeline based on GNNs demonstrated promising performance in reconstructing particle tracks in dense environments. It includes five discrete steps: data encoding, graph building, edge filtering, GNN, and track labeling. All steps were written in Python and run on both GPUs and CPUs. In this work, we accelerate the Python implementation of the pipeline through customized and commercial GPU-enabled software libraries, and develop a C++ implementation for inferencing the pipeline. The implementation features an improved, CUDA-enabled fixed-radius nearest neighbor search for graph building and a weakly connected component graph algorithm for track labeling. GNNs and other trained deep learning models are converted to ONNX and inferenced via the ONNX Runtime C++ API. The complete C++ implementation of the pipeline allows integration with existing tracking software. We report the memory usage and average event latency tracking performance of our implementation applied to the TrackML benchmark dataset. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: Proceedings submission to ACAT2021 Conference, 7 pages

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

Authors: Jorge Moreno, Shany Danieli, James S. Bullock, Robert Feldmann, Philip F. Hopkins, Onur Catmabacak, Alexander Gurvich, Alexandres Lazar, Courtney Klein, Cameron B. Hummels, Zachary Hafen, Francisco J. Mercado, Sijie Yu, Fangzhou Jiang, Coral Wheeler, Andrew Wetzel, Daniel Angles-Alcazar, Michael Boylan-Kolchin, Eliot Quataert, Claude-Andre Faucher-Giguere, Dusan Keres

Abstract: The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark ma… ▽ More The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark matter and modified gravity. Apprehension grew after several cosmological simulations using the conventional model failed to form adequate numerical analogues with comparable internal characteristics (stellar masses, sizes, velocity dispersions and morphologies). Here we show that the standard paradigm naturally produces galaxies lacking dark matter with internal characteristics in agreement with observations. Using a state-of-the-art cosmological simulation and a meticulous galaxy-identification technique, we find that extreme close encounters with massive neighbours can be responsible for this. We predict that approximately 30 percent of massive central galaxies (with at least 1e11 solar masses in stars) harbour at least one dark-matter-deficient satellite (with 1e8 - 1e9 solar masses in stars). This distinctive class of galaxies provides an additional layer in our understanding of the role of interactions in shaping galactic properties. Future observations surveying galaxies in the aforementioned regime will provide a crucial test of this scenario. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: 55 pages, 4 figures, 13 supplementary figures. 29 pages, 4 figures. Accepted for publication in Nature Astronomy. To appear on 14-February-2022. Published version: https://www.nature.com/articles/s41550-021-01598-4

Voronoi cell analysis: The shapes of particle systems

Authors: Emanuel A. Lazar, Jiayin Lu, Chris H. Rycroft

Abstract: Many physical systems can be studied as collections of particles embedded in space, evolving through deterministic evolution equations. Natural questions arise concerning how to characterize these arrangements - are they ordered or disordered? If they are ordered, how are they ordered and what kinds of defects do they possess? Originally introduced to study problems in pure mathematics, Voronoi te… ▽ More Many physical systems can be studied as collections of particles embedded in space, evolving through deterministic evolution equations. Natural questions arise concerning how to characterize these arrangements - are they ordered or disordered? If they are ordered, how are they ordered and what kinds of defects do they possess? Originally introduced to study problems in pure mathematics, Voronoi tessellations have become a powerful and versatile tool for analyzing countless problems in pure and applied physics. In this paper we explain the basics of Voronoi tessellations and the shapes they produce, and describe how they can be used to study many physical systems. △ Less

Submitted 1 June, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: 14 pages, 12 figures

Journal ref: Am. J. Phys. 90:469 (2022)

Irreducible Quotient Maps From Locally Compact Separable Metric Spaces

Authors: Aldo J. Lazar, Douglas W. B. Somerset

Abstract: Let X be a Hausdorff quotient of a standard space (that is of a locally compact separable metric space). It is shown that the following are equivalent: (i) X is the image of an irreducible quotient map from a standard space; (ii) X has a sequentially dense subset satisfying two technical conditions involving double sequences; (iii) whenever q : Y\to X is a quotient map from a standard space Y , th… ▽ More Let X be a Hausdorff quotient of a standard space (that is of a locally compact separable metric space). It is shown that the following are equivalent: (i) X is the image of an irreducible quotient map from a standard space; (ii) X has a sequentially dense subset satisfying two technical conditions involving double sequences; (iii) whenever q : Y\to X is a quotient map from a standard space Y , the restriction q_{\st}|V is an irreducible quotient map from V onto X (where q_{\st} : Y_{\st}\to X is the pure quotient derived from q, and V is the closure of the set of singleton fibres of Y_{\st}). The proof uses extensions of the theorems of Whyburn and Zarikian from compact to locally compact standard spaces. The results are new even for quotients of locally compact subsets of the real line. △ Less

Submitted 18 January, 2022; originally announced January 2022.

Comments: 26 pages

MSC Class: 54B15

Scale Reliant Inference

Authors: Michelle Pistner Nixon, Kyle C. McGovern, Jeffrey Letourneau, Lawrence A. David, Nicole A. Lazar, Sayan Mukherjee, Justin D. Silverman

Abstract: Scientific fields such as genomics, ecology, and political science often collect multivariate count data. In these fields, the data are often sufficiently noisy such that inferences regarding the total size of the measured systems have substantial uncertainty. This uncertainty can hinder downstream analyses, such as differential analysis in case-control studies. There have historically been two ap… ▽ More Scientific fields such as genomics, ecology, and political science often collect multivariate count data. In these fields, the data are often sufficiently noisy such that inferences regarding the total size of the measured systems have substantial uncertainty. This uncertainty can hinder downstream analyses, such as differential analysis in case-control studies. There have historically been two approaches to this problem: one considers the data as compositional and the other as counts that can be normalized. In this article, we use the framework of partially identified models to rigorously study the types of scientific questions (estimands) that can be answered (estimated) using these data. We prove that satisfying Frequentist inferential criteria is impossible for many estimation problems. In contrast, we find that the criteria for Bayesian inference can be satisfied, yet it requires a particular type of model called a Bayesian partially identified model. We introduce Scale Simulation Random Variables as a flexible and computationally efficient form of Bayesian partially identified models for analyzing these data. We use simulations and data analysis to validate our theory. △ Less

Submitted 5 April, 2024; v1 submitted 10 January, 2022; originally announced January 2022.

Ferrers Graphs, D-Permutations, and Surjective Staircases

Authors: Alexander Lazar

Abstract: We introduce a new family of hyperplane arrangements inspired by the homogenized Linial arrangement (which was recently introduced by Hetyei), and show that the intersection lattices of these arrangements are isomorphic to the bond lattices of Ferrers graphs. Using recent work of Lazar and Wachs we are able to give combinatorial interpretations of the characteristic polynomials of these arrangemen… ▽ More We introduce a new family of hyperplane arrangements inspired by the homogenized Linial arrangement (which was recently introduced by Hetyei), and show that the intersection lattices of these arrangements are isomorphic to the bond lattices of Ferrers graphs. Using recent work of Lazar and Wachs we are able to give combinatorial interpretations of the characteristic polynomials of these arrangements in terms of permutation enumeration. For certain infinite families of these hyperplane arrangements, we are able to give generating function formulas for their characteristic polynomials. To do so, we develop a generalization of Dumont's surjective staircases, and introduce a polynomial which enumerates these generalized surjective staircases according to several statistics. We prove a recurrence for these polynomials and show that in certain special cases this recurrence can be solved explicitly to yield a generating function. We also prove refined versions of several of these results using the theory of complex hyperplane arrangements. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 31 pages, 12 figures. Expanded version of Chapter 5 of the author's Ph.D. thesis at the University of Miami. Some results were also previously announced in the extended abstract arXiv:1811.06882

MSC Class: 52C35 (Primary); 05A05; 05A15; 05B35; 06A07

Probing the initial mass function of the first stars with transients

Authors: Alexanders Lazar, Volker Bromm

Abstract: The emergence of the first, so-called Population III (Pop III), stars shaped early cosmic history in ways that crucially depends on their initial mass function (IMF). However, because of the absence of direct observational constraints, the detailed IMF remains elusive. Nevertheless, numerical simulations agree in broad terms that the first stars were typically massive and should often end their li… ▽ More The emergence of the first, so-called Population III (Pop III), stars shaped early cosmic history in ways that crucially depends on their initial mass function (IMF). However, because of the absence of direct observational constraints, the detailed IMF remains elusive. Nevertheless, numerical simulations agree in broad terms that the first stars were typically massive and should often end their lives in violent, explosive deaths. These fates include extremely luminous pair-instability supernovae (PISNe) and bright gamma-ray bursts (GRBs), the latter arising from the collapse of rapidly rotating progenitor stars into black holes. These high-redshift transients are expected to be within the detection limits of upcoming space telescope missions, allowing to place effective constraints on the shape of the primordial IMF that is not easily accessible with other probes. This paper presents a framework to probe the Pop III IMF, utilizing the cosmological source densities of high-redshift PISNe and GRBs. Considering these transients separately could provide useful constraints on the Pop III IMF, but tighter bounds are obtainable by combining PISN and GRB counts. This combined diagnostic is more robust as it is independent of the underlying Pop III star formation rate density, an unknown prior. Future surveys promise to capture most high-redshift GRBs across the entire sky, but high-redshift PISN searches with future telescopes, e.g. Roman Space Telescope, will likely be substantially incomplete. Nevertheless, we demonstrate that even such lower bounds on the PISN count will be able to provide key constraints on the primordial IMF, in particular, if it is top-heavy or not. △ Less

Submitted 17 January, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: 11 pages, 5 figures, accepted to MNRAS

Pure quotients and Morita's theorem for $k_ω-spaces

Authors: Aldo J. Lazar, Douglas W. B. Somerset

Abstract: A $k_ω$-space $X$ is a Hausdorff quotient of a locally compact, $σ$-compact Hausdorff space. A theorem of Morita's describes the structure of $X$ when the quotient map is closed, but in 2010 a question of Arkhangel'skii's highlighted the lack of a corresponding theorem for non-closed quotient maps (even from subsets of $\mathbb{R}^n$). Arkhangel'skii's specific question had in fact been answered b… ▽ More A $k_ω$-space $X$ is a Hausdorff quotient of a locally compact, $σ$-compact Hausdorff space. A theorem of Morita's describes the structure of $X$ when the quotient map is closed, but in 2010 a question of Arkhangel'skii's highlighted the lack of a corresponding theorem for non-closed quotient maps (even from subsets of $\mathbb{R}^n$). Arkhangel'skii's specific question had in fact been answered by Siwiec in 1976, but a general structure theorem for $k_ω$-spaces is still lacking. We introduce pure quotient maps, extend Morita's theorem to these, and use Fell's topology to show that every quotient map can be 'purified' (and thus every $k_ω$-space is the image of a pure quotient map). This clarifies the structure of arbitrary $k_ω$-spaces and gives a fuller answer to Arkhangel'skii's question. △ Less

Submitted 28 July, 2021; originally announced August 2021.

Comments: To appear in the Canadian Mathematical Bulletin. arXiv admin note: substantial text overlap with arXiv:2010.03741

MSC Class: 54B15 (primary); 54D45 secondary

Opportunities and challenges in partitioning the graph measure space of real-world networks

Authors: Máté Józsa, Alpár S. Lázár, Zsolt I. Lázár

Abstract: Based on a large dataset containing thousands of real-world networks ranging from genetic, protein interaction, and metabolic networks to brain, language, ecology, and social networks we search for defining structural measures of the different complex network domains (CND). We calculate 208 measures for all networks and using a comprehensive and scrupulous workflow of statistical and machine learn… ▽ More Based on a large dataset containing thousands of real-world networks ranging from genetic, protein interaction, and metabolic networks to brain, language, ecology, and social networks we search for defining structural measures of the different complex network domains (CND). We calculate 208 measures for all networks and using a comprehensive and scrupulous workflow of statistical and machine learning methods we investigated the limitations and possibilities of identifying the key graph measures of CNDs. Our approach managed to identify well distinguishable groups of network domains and confer their relevant features. These features turn out to be CND specific and not unique even at the level of individual CNDs. The presented methodology may be applied to other similar scenarios involving highly unbalanced and skewed datasets. △ Less

Submitted 20 June, 2021; originally announced June 2021.

Comments: 11 pages, 6 figures

High-resolution population estimation using household survey data and building footprints

Authors: Gianluca Boo, Edith Darin, Douglas R Leasure, Claire A Dooley, Heather R Chamberlain, Attila N Lázár, Kevin Tschirhart, Cyrus Sinai, Nicole A Hoff, Trevon Fuller, Kamy Musene, Arly Batumbo, Anne W Rimoin, Andrew J Tatem

Abstract: The national census is an essential data source to support decision-making in many areas of public interest. However, this data may become outdated during the intercensal period, which can stretch up to several decades. We developed a Bayesian hierarchical model leveraging recent household surveys with probabilistic sampling designs and building footprints to produce up-to-date population estimate… ▽ More The national census is an essential data source to support decision-making in many areas of public interest. However, this data may become outdated during the intercensal period, which can stretch up to several decades. We developed a Bayesian hierarchical model leveraging recent household surveys with probabilistic sampling designs and building footprints to produce up-to-date population estimates. We estimated population totals and age and sex breakdowns with associated uncertainty measures within grid cells of approximately 100m in five provinces of the Democratic Republic of the Congo, a country where the last census was completed in 1984. The model exhibited a very good fit, with an R^2 value of 0.79 for out-of-sample predictions of population totals at the microcensus-cluster level and 1.00 for age and sex proportions at the province level. The results confirm the benefits of combining household surveys and building footprints for high-resolution population estimation in countries with outdated censuses. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Comments: 27 pages, 6 figures

MSC Class: 62 ACM Class: I.6.5

On the $q$-Enumeration of Barely Set-Valued Tableaux and Plane Partitions

Authors: Sam Hopkins, Alexander Lazar, Svante Linusson

Abstract: Barely set-valued tableaux are a variant of Young tableaux in which one box contains two numbers as its entry. It has recently been discovered that there are product formulas enumerating certain classes of barely set-valued tableaux. We give some $q$-analogs of these product formulas by introducing a version of major index for these tableaux. We also give product formulas and $q$-analogs for barel… ▽ More Barely set-valued tableaux are a variant of Young tableaux in which one box contains two numbers as its entry. It has recently been discovered that there are product formulas enumerating certain classes of barely set-valued tableaux. We give some $q$-analogs of these product formulas by introducing a version of major index for these tableaux. We also give product formulas and $q$-analogs for barely set-valued plane partitions. Many of the results are stated in the generality of $P$-partitions that then specialize to particularly nice formulas for rectangles and minuscule posets. The proofs use several probability distributions on the set of order ideals of a poset, depending on the real parameter $q>0$, which we think could be of independent interest. △ Less

Submitted 5 June, 2023; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: 34 pages, 6 tables, 3 figures; v2: Rewrote proof outline in Introduction, rewrote proof of Corollary 2.10, several other minor revisions at the recommendation of referees to improve exposition. To appear in European Journal of Combinatorics

MSC Class: 05A15 (Primary) 05A19; 06A07 (Secondary)

Journal ref: European Journal of Combinatorics, 113, 2023

Incentivizing Efficient Equilibria in Traffic Networks with Mixed Autonomy

Authors: Erdem Bıyık, Daniel A. Lazar, Ramtin Pedarsani, Dorsa Sadigh

Abstract: Traffic congestion has large economic and social costs. The introduction of autonomous vehicles can potentially reduce this congestion by increasing road capacity via vehicle platooning and by creating an avenue for influencing people's choice of routes. We consider a network of parallel roads with two modes of transportation: (i) human drivers, who will choose the quickest route available to them… ▽ More Traffic congestion has large economic and social costs. The introduction of autonomous vehicles can potentially reduce this congestion by increasing road capacity via vehicle platooning and by creating an avenue for influencing people's choice of routes. We consider a network of parallel roads with two modes of transportation: (i) human drivers, who will choose the quickest route available to them, and (ii) a ride hailing service, which provides an array of autonomous vehicle route options, each with different prices, to users. We formalize a model of vehicle flow in mixed autonomy and a model of how autonomous service users make choices between routes with different prices and latencies. Developing an algorithm to learn the preferences of the users, we formulate a planning optimization that chooses prices to maximize a social objective. We demonstrate the benefit of the proposed scheme by comparing the results to theoretical benchmarks which we show can be efficiently calculated. △ Less

Submitted 5 May, 2021; originally announced June 2021.

Comments: 12 pages, 7 figures, 2 tables. To appear at IEEE Transactions on Control of Network Systems (TCNS). arXiv admin note: substantial text overlap with arXiv:1904.02209

Emergent Prosociality in Multi-Agent Games Through Gifting

Authors: Woodrow Z. Wang, Mark Beliaev, Erdem Bıyık, Daniel A. Lazar, Ramtin Pedarsani, Dorsa Sadigh

Abstract: Coordination is often critical to forming prosocial behaviors -- behaviors that increase the overall sum of rewards received by all agents in a multi-agent game. However, state of the art reinforcement learning algorithms often suffer from converging to socially less desirable equilibria when multiple equilibria exist. Previous works address this challenge with explicit reward shaping, which requi… ▽ More Coordination is often critical to forming prosocial behaviors -- behaviors that increase the overall sum of rewards received by all agents in a multi-agent game. However, state of the art reinforcement learning algorithms often suffer from converging to socially less desirable equilibria when multiple equilibria exist. Previous works address this challenge with explicit reward shaping, which requires the strong assumption that agents can be forced to be prosocial. We propose using a less restrictive peer-rewarding mechanism, gifting, that guides the agents toward more socially desirable equilibria while allowing agents to remain selfish and decentralized. Gifting allows each agent to give some of their reward to other agents. We employ a theoretical framework that captures the benefit of gifting in converging to the prosocial equilibrium by characterizing the equilibria's basins of attraction in a dynamical system. With gifting, we demonstrate increased convergence of high risk, general-sum coordination games to the prosocial equilibrium both via numerical analysis and experiments. △ Less

Submitted 13 May, 2021; originally announced May 2021.

Comments: 9 pages, 6 figures, IJCAI 2021

Taking Stock of the Present and Future of Smart Technologies for Older Adults and Caregivers

Authors: Christina N. Harrington, Ben Jelen, Amanda Lazar, Aqueasha Martin-Hammond, Alisha Pradhan, Blaine Reeder, Katie Siek

Abstract: Technology has the opportunity to assist older adults as they age in place, coordinate caregiving resources, and meet unmet needs through access to resources. Currently, older adults use consumer technologies to support everyday life, however these technologies are not always accessible or as useful as they can be. Indeed, industry has attempted to create smart home technologies with older adults… ▽ More Technology has the opportunity to assist older adults as they age in place, coordinate caregiving resources, and meet unmet needs through access to resources. Currently, older adults use consumer technologies to support everyday life, however these technologies are not always accessible or as useful as they can be. Indeed, industry has attempted to create smart home technologies with older adults as a target user group, however these solutions are often more focused on the technical aspects and are short lived. In this paper, we advocate for older adults being involved in the design process - from initial ideation to product development to deployment. We encourage federally funded researchers and industry to create compensated, diverse older adult advisory boards to address stereotypes about aging while ensuring their needs are considered. We envision artificial intelligence systems that augment resources instead of replacing them - especially in under-resourced communities. Older adults rely on their caregiver networks and community organizations for social, emotional, and physical support; thus, AI should be used to coordinate resources better and lower the burden of connecting with these resources. Although sociotechnical smart systems can help identify needs of older adults, the lack of affordable research infrastructure and translation of findings into consumer technology perpetuates inequities in designing for diverse older adults. In addition, there is a disconnect between the creation of smart sensing systems and creating understandable, actionable data for older adults and caregivers to utilize. We ultimately advocate for a well-coordinated research effort across the United States that connects older adults, caregivers, community organizations, and researchers together to catalyze innovative and practical research for all stakeholders. △ Less

Submitted 31 March, 2021; originally announced April 2021.

Comments: A Computing Community Consortium (CCC) white paper, 6 pages

Report number: ccc2021whitepaper_3

The Role of Differentiation in Tolling of Traffic Networks with Mixed Autonomy

Authors: Daniel A. Lazar, Ramtin Pedarsani

Abstract: With autonomous vehicles now sharing roads with human drivers, the era of mixed autonomy brings new challenges in dealing with congestion. One cause of congestion is when vehicle users choose their routes selfishly to minimize their personal travel delay rather than a global travel delay, and prior works address this phenomenon using tolling to influence routing choices, but do not address the set… ▽ More With autonomous vehicles now sharing roads with human drivers, the era of mixed autonomy brings new challenges in dealing with congestion. One cause of congestion is when vehicle users choose their routes selfishly to minimize their personal travel delay rather than a global travel delay, and prior works address this phenomenon using tolling to influence routing choices, but do not address the setting of mixed autonomy. Tolls may be differentiated, meaning different users of a road experience different tolls, or they may be anonymous; the latter is desirable to allay concerns of fairness and privacy, as well as logistical challenges. In this work we examine the role of differentiation in traffic networks with mixed autonomy. Specifically, we first establish differentiated tolls which completely eliminate inefficiency due to selfish routing. We then show the fundamental limitations of anonymous tolls in our setting, and we provide anonymous tolls with mild performance guarantees. We show that in parallel networks, an infinitesimal differentiation in tolls is enough to guarantee optimality, and finally we establish a lower bound on the inefficiency of variable marginal cost tolling in the mixed autonomy setting. △ Less

Submitted 3 August, 2021; v1 submitted 24 March, 2021; originally announced March 2021.

Performance of a Geometric Deep Learning Pipeline for HL-LHC Particle Tracking

Authors: Xiangyang Ju, Daniel Murnane, Paolo Calafiura, Nicholas Choma, Sean Conlon, Steve Farrell, Yaoyuan Xu, Maria Spiropulu, Jean-Roch Vlimant, Adam Aurisano, V Hewes, Giuseppe Cerati, Lindsey Gray, Thomas Klijnsma, Jim Kowalkowski, Markus Atkinson, Mark Neubauer, Gage DeZoort, Savannah Thais, Aditi Chauhan, Alex Schuy, Shih-Chieh Hsu, Alex Ballow, and Alina Lazar

Abstract: The Exa.TrkX project has applied geometric learning concepts such as metric learning and graph neural networks to HEP particle tracking. Exa.TrkX's tracking pipeline groups detector measurements to form track candidates and filters them. The pipeline, originally developed using the TrackML dataset (a simulation of an LHC-inspired tracking detector), has been demonstrated on other detectors, includ… ▽ More The Exa.TrkX project has applied geometric learning concepts such as metric learning and graph neural networks to HEP particle tracking. Exa.TrkX's tracking pipeline groups detector measurements to form track candidates and filters them. The pipeline, originally developed using the TrackML dataset (a simulation of an LHC-inspired tracking detector), has been demonstrated on other detectors, including DUNE Liquid Argon TPC and CMS High-Granularity Calorimeter. This paper documents new developments needed to study the physics and computing performance of the Exa.TrkX pipeline on the full TrackML dataset, a first step towards validating the pipeline using ATLAS and CMS data. The pipeline achieves tracking efficiency and purity similar to production tracking algorithms. Crucially for future HEP applications, the pipeline benefits significantly from GPU acceleration, and its computational requirements scale close to linearly with the number of particles in the event. △ Less

Submitted 21 September, 2021; v1 submitted 11 March, 2021; originally announced March 2021.

Investigating Underlying Drivers of Variability in Residential Energy Usage Patterns with Daily Load Shape Clustering of Smart Meter Data

Authors: Ling Jin, C. Anna Spurlock, Sam Borgeson, Alina Lazar, Daniel Fredman, Annika Todd, Alexander Sim, Kesheng Wu

Abstract: Residential customers have traditionally not been treated as individual entities due to the high volatility in residential consumption patterns as well as a historic focus on aggregated loads from the utility and system feeder perspective. Large-scale deployment of smart meters has motivated increasing studies to explore disaggregated daily load patterns, which can reveal important heterogeneity a… ▽ More Residential customers have traditionally not been treated as individual entities due to the high volatility in residential consumption patterns as well as a historic focus on aggregated loads from the utility and system feeder perspective. Large-scale deployment of smart meters has motivated increasing studies to explore disaggregated daily load patterns, which can reveal important heterogeneity across different time scales, weather conditions, as well as within and across individual households. This paper aims to shed light on the mechanisms by which electricity consumption patterns exhibit variability and the different constraints that may affect demand-response (DR) flexibility. We systematically evaluate the relationship between daily time-of-use patterns and their variability to external and internal influencing factors, including time scales of interest, meteorological conditions, and household characteristics by application of an improved version of the adaptive K-means clustering method to profile "household-days" of a summer peaking utility. We find that for this summer-peaking utility, outdoor temperature is the most important external driver of the load shape variability relative to seasonality and day-of-week. The top three consumption patterns represent approximately 50% of usage on the highest temperature days. The variability in summer load shapes across customers can be explained by the responsiveness of the households to outside temperature. Our results suggest that depending on the influencing factors, not all the consumption variability can be readily translated to consumption flexibility. Such information needs to be further explored in segmenting customers for better program targeting and tailoring to meet the needs of the rapidly evolving electricity grid. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: 11 pages, 11 figures

Incentivizing Routing Choices for Safe and Efficient Transportation in the Face of the COVID-19 Pandemic

Authors: Mark Beliaev, Erdem Bıyık, Daniel A. Lazar, Woodrow Z. Wang, Dorsa Sadigh, Ramtin Pedarsani

Abstract: The COVID-19 pandemic has severely affected many aspects of people's daily lives. While many countries are in a re-opening stage, some effects of the pandemic on people's behaviors are expected to last much longer, including how they choose between different transport options. Experts predict considerably delayed recovery of the public transport options, as people try to avoid crowded places. In t… ▽ More The COVID-19 pandemic has severely affected many aspects of people's daily lives. While many countries are in a re-opening stage, some effects of the pandemic on people's behaviors are expected to last much longer, including how they choose between different transport options. Experts predict considerably delayed recovery of the public transport options, as people try to avoid crowded places. In turn, significant increases in traffic congestion are expected, since people are likely to prefer using their own vehicles or taxis as opposed to riskier and more crowded options such as the railway. In this paper, we propose to use financial incentives to set the tradeoff between risk of infection and congestion to achieve safe and efficient transportation networks. To this end, we formulate a network optimization problem to optimize taxi fares. For our framework to be useful in various cities and times of the day without much designer effort, we also propose a data-driven approach to learn human preferences about transport options, which is then used in our taxi fare optimization. Our user studies and simulation experiments show our framework is able to minimize congestion and risk of infection. △ Less

Submitted 17 February, 2021; v1 submitted 28 December, 2020; originally announced December 2020.

Comments: ICCPS 2021. 11 pages, 4 figures

Out of sight, out of mind? The impact of correlated clustering in substructure lensing

Authors: Alexandres Lazar, James S. Bullock, Michael Boylan-Kolchin, Robert Feldmann, Onur Çatmabacak, Leonidas Moustakas

Abstract: A promising route for revealing the existence of dark matter structures on mass scales smaller than the faintest galaxies is through their effect on strong gravitational lenses. We examine the role of local, lens-proximate clustering in boosting the lensing probability relative to contributions from substructure and unclustered line-of-sight (LOS) halos. Using two cosmological simulations that can… ▽ More A promising route for revealing the existence of dark matter structures on mass scales smaller than the faintest galaxies is through their effect on strong gravitational lenses. We examine the role of local, lens-proximate clustering in boosting the lensing probability relative to contributions from substructure and unclustered line-of-sight (LOS) halos. Using two cosmological simulations that can resolve halo masses of $M_{\rm halo} \simeq 10^{9}\ M_{\odot}$ (in a simulation box of length $L_{\rm box}{\sim}100\,{\rm Mpc}$) and $10^{7}\ M_{\odot}$ ($L_{\rm box}\sim20\,{\rm Mpc}$), we demonstrate that clustering in the vicinity of the lens host produces a clear enhancement relative to an assumption of unclustered halos that persists to $> 20\,R_{\rm vir}$. This enhancement exceeds estimates that use a two-halo term to account for clustering, particularly within $2-5\,R_{\rm vir}$. We provide an analytic expression for this excess, clustered contribution. We find that local clustering boosts the expected count of $10^9 \ M_\odot$ perturbing halos by ${\sim}35\%$ compared to substructure alone, a result that will significantly enhance expected signals for low-redshift ($z_l \simeq 0.2$) lenses, where substructure contributes substantially compared to LOS halos. We also find that the orientation of the lens with respect to the line of sight (e.g., whether the line of sight passes through the major axis of the lens) can also have a significant effect on the lensing signal, boosting counts by an additional $\sim 50\%$ compared to a random orientations. This could be important if discovered lenses are biased to be oriented along their principal axis. △ Less

Submitted 14 February, 2021; v1 submitted 7 December, 2020; originally announced December 2020.

Comments: Accepted to MNRAS

Voronoi chains, blocks, and clusters in perturbed square lattices

Authors: Emanuel A. Lazar, Amir Shoan

Abstract: Perturbed lattices provide simple models for studying many physical systems. In this paper we study the distribution of Voronoi chains, blocks, and clusters with prescribed combinatorial features in the perturbed square lattice, generalizing earlier work. In particular, we obtain analytic results for the presence of hexagonally-ordered regions within a square-ordered phase. Connections to high-tem… ▽ More Perturbed lattices provide simple models for studying many physical systems. In this paper we study the distribution of Voronoi chains, blocks, and clusters with prescribed combinatorial features in the perturbed square lattice, generalizing earlier work. In particular, we obtain analytic results for the presence of hexagonally-ordered regions within a square-ordered phase. Connections to high-temperature crystals and to non-equilibrium phase transitions are considered. In an appendix, we briefly consider the site-percolation threshold for this system. △ Less

Submitted 24 October, 2020; originally announced October 2020.

Comments: 15 pages, 7 figures

Journal ref: J. Stat. Mech. (2020) 103204

The structure of k_ω spaces and a question of Arhangel'skii's

Authors: A. J. Lazar, D. W. B. Somerset

Abstract: In 2010 a question of Arhangel'skii's highlighted a gap in the knowledge of k_ω-spaces. His specific question had in fact been answered by Siwiec in 1976, but the highlighted gap still remains. We introduce the simple idea of pure quotient maps, extend Morita's theorem to these, and use Fell's topology to show that every quotient map onto a k_ω-space can be 'purified'; and thus fill the gap, eluci… ▽ More In 2010 a question of Arhangel'skii's highlighted a gap in the knowledge of k_ω-spaces. His specific question had in fact been answered by Siwiec in 1976, but the highlighted gap still remains. We introduce the simple idea of pure quotient maps, extend Morita's theorem to these, and use Fell's topology to show that every quotient map onto a k_ω-space can be 'purified'; and thus fill the gap, elucidate the structure of k_ω-spaces, and obtain a fuller answer to Arhangel'skii's question. △ Less

Submitted 13 October, 2021; v1 submitted 7 October, 2020; originally announced October 2020.

Comments: The paper was submitted under another title. The readers should instead refer to arXiv:2108.05202

MSC Class: 54B15

Optimal Tolling for Multitype Mixed Autonomous Traffic Networks

Authors: Daniel A. Lazar, Ramtin Pedarsani

Abstract: When selfish users share a road network and minimize their individual travel costs, the equilibrium they reach can be worse than the socially optimal routing. Tolls are often used to mitigate this effect in traditional congestion games, where all vehicle contribute identically to congestion. However, with the proliferation of autonomous vehicles and driver-assistance technology, vehicles become he… ▽ More When selfish users share a road network and minimize their individual travel costs, the equilibrium they reach can be worse than the socially optimal routing. Tolls are often used to mitigate this effect in traditional congestion games, where all vehicle contribute identically to congestion. However, with the proliferation of autonomous vehicles and driver-assistance technology, vehicles become heterogeneous in how they contribute to road latency. This magnifies the potential inefficiencies due to selfish routing and invalidates traditional tolling methods. To address this, we consider a network of parallel roads where the latency on each road is an affine function of the quantity of flow of each vehicle type. We provide tolls (which differentiate between vehicle types) which are guaranteed to minimize social cost at equilibrium. The tolls are a function of a calculated optimal routing; to enable this tolling, we prove that some element in the set of optimal routings has a lack of cycles in a graph representing the way vehicles types share roads. We then show that unless a planner can differentiate between vehicle types in the tolls given, the resulting equilibrium can be unboundedly worse than the optimal routing, and that marginal cost tolling fails in our setting. △ Less

Submitted 31 August, 2020; originally announced September 2020.

The time-scales probed by star formation rate indicators for realistic, bursty star formation histories from the FIRE simulations

Authors: José A. Flores Velázquez, Alexander B. Gurvich, Claude-André Faucher-Giguère, James S. Bullock, Tjitske K. Starkenburg, Jorge Moreno, Alexandres Lazar, Francisco J. Mercado, Jonathan Stern, Martin Sparre, Christopher C. Hayward, Andrew Wetzel, Kareem El-Badry

Abstract: Understanding the rate at which stars form is central to studies of galaxy formation. Observationally, the star formation rates (SFRs) of galaxies are measured using the luminosity in different frequency bands, often under the assumption of a time-steady SFR in the recent past. We use star formation histories (SFHs) extracted from cosmological simulations of star-forming galaxies from the FIRE pro… ▽ More Understanding the rate at which stars form is central to studies of galaxy formation. Observationally, the star formation rates (SFRs) of galaxies are measured using the luminosity in different frequency bands, often under the assumption of a time-steady SFR in the recent past. We use star formation histories (SFHs) extracted from cosmological simulations of star-forming galaxies from the FIRE project to analyze the time-scales to which the H$α$ and far-ultraviolet (FUV) continuum SFR indicators are sensitive. In these simulations, the SFRs are highly time variable for all galaxies at high redshift, and continue to be bursty to z=0 in dwarf galaxies. When FIRE SFHs are partitioned into their bursty and time-steady phases, the best-fitting FUV time-scale fluctuates from its ~10 Myr value when the SFR is time-steady to >~100 Myr immediately following particularly extreme bursts of star formation during the bursty phase. On the other hand, the best-fitting averaging time-scale for H$α$ is generally insensitive to the SFR variability in the FIRE simulations and remains ~5 Myr at all times. These time-scales are shorter than the 100 Myr and 10 Myr time-scales sometimes assumed in the literature for FUV and H$α$, respectively, because while the FUV emission persists for stellar populations older than 100 Myr, the time-dependent luminosities are strongly dominated by younger stars. Our results confirm that the ratio of SFRs inferred using H$α$ vs. FUV can be used to probe the burstiness of star formation in galaxies. △ Less

Submitted 15 December, 2020; v1 submitted 19 August, 2020; originally announced August 2020.

Comments: 14 pages, 10 figures, accepted to MNRAS

Distribution of Topological Types in Grain-Growth Microstructures

Authors: Emanuel A. Lazar, Jeremy K. Mason, Robert D. MacPherson, David J. Srolovitz

Abstract: An open question in studying normal grain growth concerns the asymptotic state to which microstructures converge. In particular, the distribution of grain topologies is unknown. We introduce a thermodynamic-like theory to explain these distributions in two- and three-dimensional systems. In particular, a bending-like energy $E_i$ is associated to each grain topology $t_i$, and the probability of o… ▽ More An open question in studying normal grain growth concerns the asymptotic state to which microstructures converge. In particular, the distribution of grain topologies is unknown. We introduce a thermodynamic-like theory to explain these distributions in two- and three-dimensional systems. In particular, a bending-like energy $E_i$ is associated to each grain topology $t_i$, and the probability of observing that particular topology is proportional to $\frac{1}{s(t_i)}e^{-βE_i}$, where $s(t_i)$ is the order of an associated symmetry group and $β$ is a thermodynamic-like constant. We explain the physical origins of this approach, and provide numerical evidence in support. △ Less

Submitted 4 July, 2020; originally announced July 2020.

Comments: 6 pages, 5 figures

Journal ref: Phys. Rev. Lett. 125, 015501 (2020)

Hypothesis Testing for Shapes using Vectorized Persistence Diagrams

Authors: Chul Moon, Nicole A. Lazar

Abstract: Topological data analysis involves the statistical characterization of the shape of data. Persistent homology is a primary tool of topological data analysis, which can be used to analyze topological features and perform statistical inference. In this paper, we present a two-stage hypothesis test for vectorized persistence diagrams. The first stage filters vector elements in the vectorized persiste… ▽ More Topological data analysis involves the statistical characterization of the shape of data. Persistent homology is a primary tool of topological data analysis, which can be used to analyze topological features and perform statistical inference. In this paper, we present a two-stage hypothesis test for vectorized persistence diagrams. The first stage filters vector elements in the vectorized persistence diagrams to enhance the power of the test. The second stage consists of multiple hypothesis tests, with false positives controlled by false discovery rates. We demonstrate the flexibility of our method by applying it to a variety of simulated and real-world data types. Our results show that the proposed hypothesis test enables accurate and informative inferences on the shape of data compared to the existing hypothesis testing methods for persistent homology. △ Less

Submitted 28 February, 2023; v1 submitted 9 June, 2020; originally announced June 2020.

A dark matter profile to model diverse feedback-induced core sizes of $Λ$CDM haloes

Authors: Alexandres Lazar, James S. Bullock, Michael Boylan-Kolchin, T. K. Chan, Philip F. Hopkins, Andrew S. Graus, Andrew Wetzel, Kareem El-Badry, Coral Wheeler, Maria C. Straight, Dušan Kereš, Claude-André Faucher-Giguère, Alex Fitts, Shea Garrison-Kimmel

Abstract: We analyze the cold dark matter density profiles of 54 galaxy halos simulated with FIRE-2 galaxy formation physics, each resolved within $0.5\%$ of the halo virial radius. These halos contain galaxies with masses that range from ultra-faint dwarfs ($M_\star \simeq 10^{4.5} M_{\odot}$) to the largest spirals ($M_\star \simeq 10^{11} M_{\odot}$) and have density profiles that are both cored and cusp… ▽ More We analyze the cold dark matter density profiles of 54 galaxy halos simulated with FIRE-2 galaxy formation physics, each resolved within $0.5\%$ of the halo virial radius. These halos contain galaxies with masses that range from ultra-faint dwarfs ($M_\star \simeq 10^{4.5} M_{\odot}$) to the largest spirals ($M_\star \simeq 10^{11} M_{\odot}$) and have density profiles that are both cored and cuspy. We characterize our results using a new analytic density profile that extends the standard Einasto form to allow for a pronounced constant-density core in the resolved innermost radius. With one additional core-radius parameter, $r_{c}$, this "core-Einasto" profile is able to characterize the shape and normalization of our feedback-impacted dark matter halos. In order to enable comparisons with observations, we provide fitting functions for $r_{c}$ and other profile parameters as a function of both $M_\star$ and $M_{\star}/M_{\rm halo}$. In agreement with similar studies done in the literature, we find that dark matter core formation is most efficient at the characteristic stellar-mass to halo-mass ratio $M_\star/M_{\rm halo} \simeq 5 \times 10^{-3}$, or $M_{\star} \sim 10^9 \, M_{\odot}$, with cores that are roughly the size of the galaxy half-light radius, $r_{c} \simeq 1-5$ kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with $M_{\star}/M_{\rm halo} < 5\times 10^{-4}$ or $M_\star \lesssim 10^6 \, M_{\odot}$. For Milky Way-size galaxies, baryonic contraction often makes halos significantly more concentrated and dense at the stellar half-light radius than dark matter only runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of $\simeq 0.5-2$ kpc in size. Recent evidence for a ${\sim} 2$ kpc core in the Milky Way's dark matter halo is consistent with this expectation. △ Less

Submitted 8 July, 2020; v1 submitted 22 April, 2020; originally announced April 2020.

Comments: 27 pages; 19 figures; Accepted by MNRAS

Partition and Cohen-Macaulay Extenders

Authors: Joseph Doolittle, Bennet Goeckner, Alexander Lazar

Abstract: If a pure simplicial complex is partitionable, then its $h$-vector has a combinatorial interpretation in terms of any partitioning of the complex. Given a non-partitionable complex $Δ$, we construct a complex $Γ\supseteq Δ$ of the same dimension such that both $Γ$ and the relative complex $(Γ,Δ)$ are partitionable. This allows us to rewrite the $h$-vector of any pure simplicial complex as the diff… ▽ More If a pure simplicial complex is partitionable, then its $h$-vector has a combinatorial interpretation in terms of any partitioning of the complex. Given a non-partitionable complex $Δ$, we construct a complex $Γ\supseteq Δ$ of the same dimension such that both $Γ$ and the relative complex $(Γ,Δ)$ are partitionable. This allows us to rewrite the $h$-vector of any pure simplicial complex as the difference of two $h$-vectors of partitionable complexes, giving an analogous interpretation of the $h$-vector of a non-partitionable complex. By contrast, for a given complex $Δ$ it is not always possible to find a complex $Γ$ such that both $Γ$ and $(Γ,Δ)$ are Cohen-Macaulay. We characterize when this is possible, and we show that the construction of such a $Γ$ in this case is remarkably straightforward. We end with a note on a similar notion for shellability and a connection to Simon's conjecture on extendable shellability for uniform matroids. △ Less

Submitted 28 October, 2021; v1 submitted 28 November, 2019; originally announced November 2019.

Comments: 14 pages, 4 figures. (V3): Corrected the remark on nonpure partitionability and expanded it into a new Section 5; clarified the proof of Proposition 6.2 and fixed an off-by-one error; other minor corrections (all based on referee suggestions). To appear in the European Journal of Combinatorics

MSC Class: 05E45; 13F55

The Homogenized Linial Arrangement and Genocchi Numbers

Authors: Alexander Lazar, Michelle L. Wachs

Abstract: We study the intersection lattice of a hyperplane arrangement recently introduced by Hetyei who showed that the number of regions of the arrangement is a median Genocchi number. Using a different method, we refine Hetyei's result by providing a combinatorial interpretation of the coefficients of the characteristic polynomial of the intersection lattice of this arrangement. We also show that the Mö… ▽ More We study the intersection lattice of a hyperplane arrangement recently introduced by Hetyei who showed that the number of regions of the arrangement is a median Genocchi number. Using a different method, we refine Hetyei's result by providing a combinatorial interpretation of the coefficients of the characteristic polynomial of the intersection lattice of this arrangement. We also show that the Möbius invariant of the intersection lattice is a (nonmedian) Genocchi number. The Genocchi numbers count a class of permutations known as Dummont permutations and the median Genocchi numbers count the derangements in this class. We show that the signless coefficients of the characteristic polynomial count Dumont-like permutations with a given number of cycles. This enables us to derive formulas for the generating function of the characteristic polynomial, which reduce to known formulas for the generating functions of the Genocchi numbers and the median Genocchi numbers. As a byproduct of our work, we obtain new models for the Genocchi and median Genocchi numbers. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: 33 pages, 10 figures. Proofs of Type A results announced in the extended abstract arXiv:1811.06882, plus some additional results

MSC Class: 52C35 (Primary); 05A05; 05A15; 05B35; 06A07; 11B68 (Secondary)

Sparse Identification of Contrast Gain Control in the Fruit Fly Photoreceptor and Amacrine Cell Layer

Authors: Aurel A. Lazar, Nikul H. Ukani, Yiyin Zhou

Abstract: The fruit fly's natural visual environment is often characterized by light intensities ranging across several orders of magnitude and by rapidly varying contrast across space and time. Fruit fly photoreceptors robustly transduce and, in conjunction with amacrine cells, process visual scenes and provide the resulting signal to downstream targets. Here we model the first step of visual processing in… ▽ More The fruit fly's natural visual environment is often characterized by light intensities ranging across several orders of magnitude and by rapidly varying contrast across space and time. Fruit fly photoreceptors robustly transduce and, in conjunction with amacrine cells, process visual scenes and provide the resulting signal to downstream targets. Here we model the first step of visual processing in the photoreceptor-amacrine cell layer. We propose a novel divisive normalization processor (DNP) for modeling the computation taking place in the photoreceptor-amacrine cell layer. The DNP explicitly models the photoreceptor feedforward and temporal feedback processing paths and the spatio-temporal feedback path of the amacrine cells. We then formally characterize the contrast gain control of the DNP and provide sparse identification algorithms that can efficiently identify each the feedforward and feedback DNP components. The algorithms presented here are the first demonstration of tractable and robust identification of the components of a divisive normalization processor. The sparse identification algorithms can be readily employed in experimental settings, and their effectiveness is demonstrated with several examples. △ Less

Submitted 3 October, 2019; originally announced October 2019.

Optimal Tolling for Heterogeneous Traffic Networks with Mixed Autonomy

Authors: Daniel A. Lazar, Samuel Coogan, Ramtin Pedarsani

Abstract: When people pick routes to minimize their travel time, the total experienced delay, or social cost, may be significantly greater than if people followed routes assigned to them by a social planner. This effect is accentuated when human drivers share roads with autonomous vehicles. When routed optimally, autonomous vehicles can make traffic networks more efficient, but when acting selfishly, the in… ▽ More When people pick routes to minimize their travel time, the total experienced delay, or social cost, may be significantly greater than if people followed routes assigned to them by a social planner. This effect is accentuated when human drivers share roads with autonomous vehicles. When routed optimally, autonomous vehicles can make traffic networks more efficient, but when acting selfishly, the introduction of autonomous vehicles can actually worsen congestion. We seek to mitigate this effect by influencing routing choices via tolling. We consider a network of parallel roads with affine latency functions that are heterogeneous, meaning that the increase in capacity due to to the presence of autonomous vehicles may vary from road to road. We show that if human drivers and autonomous users have the same tolls, the social cost may be arbitrarily worse than when optimally routed. We then prove qualities of the optimal routing and use them to design tolls that are guaranteed to minimize social cost at equilibrium. To the best of our knowledge, this is the first tolling scheme that yields a unique socially optimal equilibrium for parallel heterogeneous network with affine latency functions. △ Less

Submitted 25 September, 2019; originally announced September 2019.

Comments: Accepted to Conference on Decision and Control, 2019

Method of Moments Histograms

Authors: James S. Weber, Nicole A. Lazar

Abstract: Uniform bin width histograms are widely used so this data graphic should represent data as correctly as possible. Method of moments based on familiar mean, variance and Fisher-Pearson skewness cure this problem. Uniform bin width histograms are widely used so this data graphic should represent data as correctly as possible. Method of moments based on familiar mean, variance and Fisher-Pearson skewness cure this problem. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Comments: Updates: Weber, J. S. (2016) - What Can We Learn from Correct Calculation of Histograms? - In JSM Proceedings, Stat. Graphics Sect.. Alexandria, VA: Amer. Stat. Assoc. 1893-1913. Weber, J. S. (2016) - Calculating Method of Moments Uniform Bin Width Histograms. - arXiv:1606.04891

MSC Class: 62-07; 62-09

Learning How to Dynamically Route Autonomous Vehicles on Shared Roads

Authors: Daniel A. Lazar, Erdem Bıyık, Dorsa Sadigh, Ramtin Pedarsani

Abstract: Road congestion induces significant costs across the world, and road network disturbances, such as traffic accidents, can cause highly congested traffic patterns. If a planner had control over the routing of all vehicles in the network, they could easily reverse this effect. In a more realistic scenario, we consider a planner that controls autonomous cars, which are a fraction of all present cars.… ▽ More Road congestion induces significant costs across the world, and road network disturbances, such as traffic accidents, can cause highly congested traffic patterns. If a planner had control over the routing of all vehicles in the network, they could easily reverse this effect. In a more realistic scenario, we consider a planner that controls autonomous cars, which are a fraction of all present cars. We study a dynamic routing game, in which the route choices of autonomous cars can be controlled and the human drivers react selfishly and dynamically. As the problem is prohibitively large, we use deep reinforcement learning to learn a policy for controlling the autonomous vehicles. This policy indirectly influences human drivers to route themselves in such a way that minimizes congestion on the network. To gauge the effectiveness of our learned policies, we establish theoretical results characterizing equilibria and empirically compare the learned policy results with best possible equilibria. We prove properties of equilibria on parallel roads and provide a polynomial-time optimization for computing the most efficient equilibrium. Moreover, we show that in the absence of these policies, high demand and network perturbations would result in large congestion, whereas using the policy greatly decreases the travel times by minimizing the congestion. To the best of our knowledge, this is the first work that employs deep reinforcement learning to reduce congestion by indirectly influencing humans' routing decisions in mixed-autonomy traffic. △ Less

Submitted 3 June, 2021; v1 submitted 9 September, 2019; originally announced September 2019.

Comments: Accepted to Transportation Research Part C

Accurate mass estimates from the proper motions of dispersion-supported galaxies

Authors: Alexandres Lazar, James S. Bullock

Abstract: Starting with the spherical Jeans equation, we show that there exists a radius where the mass enclosed depends only on the projected tangential velocity dispersion, assuming that the anisotropy profile slowly varies. This is well-approximated at the radius where the log-slope of the stellar tracer profile is $-2$: $r_{-2}$. The associated mass is… ▽ More Starting with the spherical Jeans equation, we show that there exists a radius where the mass enclosed depends only on the projected tangential velocity dispersion, assuming that the anisotropy profile slowly varies. This is well-approximated at the radius where the log-slope of the stellar tracer profile is $-2$: $r_{-2}$. The associated mass is $M(r_{-2}) = 2 G^{-1} \langle σ_{\mathcal{T}}^{2}\rangle^{*} r_{-2}$ and the circular velocity is $V^{2}({r_{-2}}) = 2\langle σ_{\mathcal{T}}^{2}\rangle^{*}$. For a Plummer profile $r_{-2} \simeq 4R_e/5$. Importantly, $r_{-2}$ is smaller than the characteristic radius for line-of-sight velocities derived by Wolf et al. 2010. Together, the two estimators can constrain the mass profiles of dispersion-supported galaxies. We illustrate its applicability using published proper motion measurements of dwarf galaxies Draco and Sculptor, and find that they are consistent with inhabiting cuspy NFW subhalos of the kind predicted in CDM but we cannot rule out a core. We test our combined mass estimators against previously-published, non-spherical cosmological dwarf galaxy simulations done in both CDM and SIDM. For CDM, the estimates for the dynamic rotation curves are found to be accurate to $10\%$ while SIDM are accurate to $15\%$. Unfortunately, this level of accuracy is not good enough to measure slopes at the level required to distinguish between cusps and cores of the type predicted in viable SIDM models without stronger priors. However, we find that this provides good enough accuracy to distinguish between the normalization differences predicted at small radii ($r \simeq r_{-2} < r_{\rm core}$) for interesting SIDM models. As the number of galaxies with internal proper motions increases, mass estimators of this kind will enable valuable constraints on SIDM and CDM models. △ Less

Submitted 9 March, 2020; v1 submitted 20 July, 2019; originally announced July 2019.

Comments: Accepted to MNRAS

Dark and luminous satellites of LMC-mass galaxies in the FIRE simulations

Authors: Ethan D. Jahn, Laura V. Sales, Andrew Wetzel, Michael Boylan-Kolchin, T. K. Chan, Kareem El-Badry, Alexandres Lazar, James S. Bullock

Abstract: Within $Λ$CDM, dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhalos, several of which should host faint dwarf companions. Recent Gaia proper motions confirm new members of the LMC-system in addition to the previously known SMC, including two classical dwarf galaxies ($M_\ast$ > $10^5$ M$_{\odot}$; Carina and Fornax) as well as several ultra-faint… ▽ More Within $Λ$CDM, dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhalos, several of which should host faint dwarf companions. Recent Gaia proper motions confirm new members of the LMC-system in addition to the previously known SMC, including two classical dwarf galaxies ($M_\ast$ > $10^5$ M$_{\odot}$; Carina and Fornax) as well as several ultra-faint dwarfs (Car2, Car3, Hor1, and Hyd1). We use the Feedback In Realistic Environments (FIRE) simulations to study the dark and luminous (down to ultrafaint masses, $M_\ast$ ~ $6\times10^3$ M$_{\odot}$) substructure population of isolated LMC-mass hosts ($M_\text{200}$ = $1-3\times10^{11}$ M$_{\odot}$) and place the Gaia + DES results in a cosmological context. By comparing number counts of subhalos in simulations with and without baryons, we find that, within 0.2 $r_\text{200}$, LMC-mass hosts deplete ~30% of their substructure, significantly lower than the ~70% of substructure depleted by Milky Way (MW) mass hosts. For our highest resolution runs ($m_\text{bary}$ = 880 M$_{\odot}$), ~5-10 subhalos form galaxies with $M_\ast$ > $10^4$ M$_{\odot}$, in agreement with the 7 observationally inferred pre-infall LMC companions. However, we find steeper simulated luminosity functions than observed, hinting at observation incompleteness at the faint end. The predicted DM content for classical satellites in FIRE agrees with observed estimates for Carina and Fornax, supporting the case for an LMC association. We predict that tidal stripping within the LMC potential lowers the inner dark matter density of ultra faint companions of the LMC. Thus, in addition to their orbital consistency, the low densities of dwarfs Car2, Hyd1, and Hyd2 reinforce their likelihood of Magellanic association. △ Less

Submitted 31 August, 2019; v1 submitted 5 July, 2019; originally announced July 2019.

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

The Green Choice: Learning and Influencing Human Decisions on Shared Roads

Authors: Erdem Bıyık, Daniel A. Lazar, Dorsa Sadigh, Ramtin Pedarsani

Abstract: Autonomous vehicles have the potential to increase the capacity of roads via platooning, even when human drivers and autonomous vehicles share roads. However, when users of a road network choose their routes selfishly, the resulting traffic configuration may be very inefficient. Because of this, we consider how to influence human decisions so as to decrease congestion on these roads. We consider a… ▽ More Autonomous vehicles have the potential to increase the capacity of roads via platooning, even when human drivers and autonomous vehicles share roads. However, when users of a road network choose their routes selfishly, the resulting traffic configuration may be very inefficient. Because of this, we consider how to influence human decisions so as to decrease congestion on these roads. We consider a network of parallel roads with two modes of transportation: (i) human drivers who will choose the quickest route available to them, and (ii) ride hailing service which provides an array of autonomous vehicle ride options, each with different prices, to users. In this work, we seek to design these prices so that when autonomous service users choose from these options and human drivers selfishly choose their resulting routes, road usage is maximized and transit delay is minimized. To do so, we formalize a model of how autonomous service users make choices between routes with different price/delay values. Developing a preference-based algorithm to learn the preferences of the users, and using a vehicle flow model related to the Fundamental Diagram of Traffic, we formulate a planning optimization to maximize a social objective and demonstrate the benefit of the proposed routing and learning scheme. △ Less

Submitted 9 April, 2019; v1 submitted 3 April, 2019; originally announced April 2019.

Comments: Submitted to CDC 2019