Showing 1–50 of 97 results for author: Paranjape, A

Dynamics of the response of dark matter halo to galaxy evolution in IllustrisTNG

Authors: Premvijay Velmani, Aseem Paranjape

Abstract: We present the dynamical evolution of the dark matter's relaxation response to galaxies and their connection to the astrophysical properties as simulated in the IllustrisTNG suite of cosmological hydrodynamical simulations. Our results show that the radially-dependent linear relaxation relation model from our previous work is applicable at least from redshift $z=5$. We focus on the offset paramete… ▽ More We present the dynamical evolution of the dark matter's relaxation response to galaxies and their connection to the astrophysical properties as simulated in the IllustrisTNG suite of cosmological hydrodynamical simulations. Our results show that the radially-dependent linear relaxation relation model from our previous work is applicable at least from redshift $z=5$. We focus on the offset parameter $q_0$, which characterizes the relaxation of dark matter shells without changing the enclosed mass. We perform multiple time-series analyses to determine the possible causal connections between the relaxation mechanism and astrophysical processes such as star formation and associated feedback processes, as well as feedback due to active galactic nuclei. We show that star formation activity significantly influences the halo relaxation response throughout its evolutionary history, with essentially immediate effects in the inner haloes and delayed effects of 2 to 3 Gyr in the outer regions. Metal content shows a weaker connection to relaxation than star formation rates, but the accumulated wind from feedback processes exhibits a stronger correlation. These findings enhance our understanding of halo relaxation mechanisms. Our estimates of the time-scales relevant for dark matter relaxation can potentially improve the description of halo profiles in existing baryonification schemes and semi-analytical galaxy formation models. Our results also show how the relaxation response of dark haloes can probe the evolutionary history of the galaxies they host. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 22 pages, 11 figures; prepared for submission to JCAP

Eigen-decomposition of Covariance matrices: An application to the BAO Linear Point

Authors: Jaemyoung Jason Lee, Farnik Nikakhtar, Aseem Paranjape, Ravi K. Sheth

Abstract: The Baryon Acoustic Oscillation (BAO) feature in the two-point correlation function (TPCF) of discrete tracers such as galaxies is an accurate standard ruler. The covariance matrix of the TPCF plays an important role in determining how the precision of this ruler depends on the number density and clustering strength of the tracers, as well as the survey volume. An eigen-decomposition of this matri… ▽ More The Baryon Acoustic Oscillation (BAO) feature in the two-point correlation function (TPCF) of discrete tracers such as galaxies is an accurate standard ruler. The covariance matrix of the TPCF plays an important role in determining how the precision of this ruler depends on the number density and clustering strength of the tracers, as well as the survey volume. An eigen-decomposition of this matrix provides an objective way to separate the contributions of cosmic variance from those of shot-noise to the statistical uncertainties. For the signal-to-noise levels that are expected in ongoing and next-generation surveys, the cosmic variance eigen-modes dominate. These modes are smooth functions of scale, meaning that: they are insensitive to the modest changes in binning that are allowed if one wishes to resolve the BAO feature in the TPCF; they provide a good description of the correlated residuals which result from fitting smooth functional forms to the measured TPCF; they motivate a simple but accurate approximation for the uncertainty on the Linear Point (LP) estimate of the BAO distance scale. This approximation allows one to quantify the precision of the BAO distance scale estimate without having to generate a large ensemble of mock catalogs and explains why: the uncertainty on the LP does not depend on the functional form fitted to the TPCF or the binning used; the LP is more constraining than the peak or dip scales in the TPCF; the evolved TPCF is less constraining than the initial one, so that reconstruction schemes can yield significant gains in precision. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 10 pages, 8 figures, submitted to Physical Review D

Optimization of cosmic filament finders and unbiased recovery of filament phase space profiles using mock filaments

Authors: Saee Dhawalikar, Aseem Paranjape

Abstract: Cosmic filaments, the most prominent features of the cosmic web, possibly hold untapped potential for cosmological inference. While it is natural to expect the structure of filaments to show universality similar to that seen in dark matter halos, the lack of agreement between different filament finders on what constitutes a filament has hampered progress on this topic. We initiate a programme to s… ▽ More Cosmic filaments, the most prominent features of the cosmic web, possibly hold untapped potential for cosmological inference. While it is natural to expect the structure of filaments to show universality similar to that seen in dark matter halos, the lack of agreement between different filament finders on what constitutes a filament has hampered progress on this topic. We initiate a programme to systematically investigate and uncover possible universal features in the phase space structure of cosmic filaments, by generating particle realizations of mock filaments with $\textit{a priori}$ known properties. Using these, we identify an important source of bias in the extraction of radial density profiles, which occurs when the local curvature $κ$ of the spine exceeds a threshold determined by the filament thickness. This bias exists even for perfectly determined spines, thus affecting $\textit{all}$ filament finders. We show that this bias can be nearly eliminated by simply discarding the regions with the highest $κ$, with little loss of precision. An additional source of bias is the noise generated by the filament finder when identifying the spine, which depends on both the finder algorithm as well as intrinsic properties of the individual filament. We find that, to mitigate this bias, it is essential not only to smooth the estimated spine, but to $\textit{optimize}$ this smoothing separately for each filament. We propose a novel optimization based on minimizing the estimated filament thickness, along with Fourier space smoothing. We implement these techniques using two tools, $\texttt{FilGen}$ which generates mock filaments and $\texttt{FilAPT}$ which analyses and processes them. We expect these tools to be useful in calibrating the performance of filament finders, thereby enabling searches for filament universality. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: 35 pages, 13 figures; prepared for submission to JCAP

Voice-Based Smart Assistant System for Vehicles using RASA

Authors: Aditya Paranjape, Yash Patwardhan, Vedant Deshpande, Aniket Darp, Jayashree Jagdale

Abstract: Conversational AIs, or chatbots, mimic human speech when conversing. Smart assistants facilitate the automation of several tasks that needed human intervention earlier. Because of their accuracy, absence of dependence on human resources, and accessibility around the clock, chatbots can be employed in vehicles too. Due to people's propensity to divert their attention away from the task of driving w… ▽ More Conversational AIs, or chatbots, mimic human speech when conversing. Smart assistants facilitate the automation of several tasks that needed human intervention earlier. Because of their accuracy, absence of dependence on human resources, and accessibility around the clock, chatbots can be employed in vehicles too. Due to people's propensity to divert their attention away from the task of driving while engaging in other activities like calling, playing music, navigation, and getting updates on the weather forecast and latest news, road safety has declined and accidents have increased as a result. It would be advantageous to automate these tasks using voice commands rather than carrying them out manually. This paper focuses on the development of a voice-based smart assistance application for vehicles based on the RASA framework. The smart assistant provides functionalities like navigation, communication via calls, getting weather forecasts and the latest news updates, and music that are completely voice-based in nature. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: 6 pages, 4 figures, accepted at IEEE International Conference on Computational Intelligence, Networks and Security ICCINS-2023

A self-similar model of galaxy formation and dark halo relaxation

Authors: Premvijay Velmani, Aseem Paranjape

Abstract: We develop a spherical self-similar model for the formation of a galaxy through gas collapsing in an isolated self-gravitating dark matter halo. We improve upon the existing literature on self-similar collapse in two ways. First, we include the effects of radiative cooling and the formation of a pseudo-disk at the center of collapse, in a parametrised manner. More importantly, we solve for the evo… ▽ More We develop a spherical self-similar model for the formation of a galaxy through gas collapsing in an isolated self-gravitating dark matter halo. We improve upon the existing literature on self-similar collapse in two ways. First, we include the effects of radiative cooling and the formation of a pseudo-disk at the center of collapse, in a parametrised manner. More importantly, we solve for the evolution of gas and dark matter simultaneously and self-consistently using a novel iterative approach. As a result, our model produces shell trajectories of both gas and dark matter that qualitatively agree with the results of full hydrodynamical simulations. We discuss the impact of various ingredients such as the accretion rate, gas equation of state, disk radius and cooling rate amplitude on the evolution of the gas shells. The self-consistent evolution of gas and dark matter allows us to study the response of the dark matter trajectories to the presence of collapsing gas, an effect that has gained increasing importance recently in the context of precision estimates of small-scale statistics like the matter power spectrum. Our default configuration produces a relaxation relation in qualitative agreement with that seen in cosmological hydrodynamical simulations, and further allows us to easily study the impact of the model ingredients mentioned above. As an initial application, we vary one ingredient at a time and find that the accretion rate and gas equation of state have the largest impact on the relaxation relation, while the cooling amplitude plays only a minor role. Our model thus provides a convenient framework to rapidly explore the coupled nonlinear impact of multiple astrophysical processes on the mass and velocity profiles of dark matter in galactic halos. (Abridged) △ Less

Submitted 3 April, 2024; v1 submitted 23 November, 2023; originally announced November 2023.

Comments: 23 pages, 9 figures; v2: added discussion, accepted in JCAP

A GPR-Based Emulator for Semi-numerical Reionization Code SCRIPT: Parameter Inference from 21 cm Data

Authors: T. Roy Choudhury, A. Paranjape, B. Maity

Abstract: Semi-numerical models of reionization typically involve a large number of unknown parameters whose values are constrained by comparing with observations. Increasingly often, exploring this parameter space using semi-numerical simulations can become computationally intensive, thus necessitating the use of emulators. In this work, we present a likelihood emulator based on Gaussian Process Regression… ▽ More Semi-numerical models of reionization typically involve a large number of unknown parameters whose values are constrained by comparing with observations. Increasingly often, exploring this parameter space using semi-numerical simulations can become computationally intensive, thus necessitating the use of emulators. In this work, we present a likelihood emulator based on Gaussian Process Regression (GPR) for our semi-numerical reionization code, SCRIPT, and use it for parameter inference using mock 21 cm power spectrum data and Bayesian MCMC analysis. A unique aspect of our methodology is the utilization of coarse resolution simulations to identify high-probability regions within the parameter space, employing only a moderate amount of computational time. Samples drawn from these high-probability regions are used to construct the training set for the emulator. The subsequent MCMC using this GPR-trained emulator is found to provide parameter posteriors that agree reasonably well with those obtained using conventional MCMC. The computing time for the analysis, which includes both generation of training sets and training the emulator, is reduced by approximately an order of magnitude. This methodology is particularly advantageous in scenarios where one wants to use different parametrizations of reionization models and/or needs to start with broad prior distributions on the parameters, offering an efficient and effective means of parameter inference. △ Less

Submitted 2 February, 2024; v1 submitted 8 November, 2023; originally announced November 2023.

Comments: Accepted for publication in JCAP

Covariance matrices for the Lyman-$α$ forest using the lognormal approximation

Authors: Bhaskar Arya, Aseem Paranjape, Tirthankar Roy Choudhury

Abstract: We investigate the nature of correlations in the small-scale flux statistics of the Lyman-$α$ (Ly$α$) forest across redshift bins. Understanding these correlations is important for unbiased cosmological and astrophysical parameter inference using the Ly$α$ forest. We focus on the 1-dimensional flux power spectrum (FPS) and mean flux ($\bar F$) simulated using the semi-numerical lognormal model we… ▽ More We investigate the nature of correlations in the small-scale flux statistics of the Lyman-$α$ (Ly$α$) forest across redshift bins. Understanding these correlations is important for unbiased cosmological and astrophysical parameter inference using the Ly$α$ forest. We focus on the 1-dimensional flux power spectrum (FPS) and mean flux ($\bar F$) simulated using the semi-numerical lognormal model we developed in earlier work. The lognormal model can capture the effects of long wavelength modes with relative ease as compared to full smoothed particle hydrodynamical (SPH) simulations that are limited by box volume. For a single redshift bin of size $Δz\simeq 0.1$, we show that the lognormal model predicts positive cross-correlations between $k$-bins in the FPS, and a negative correlation for $\bar F\times$ FPS, in qualitative agreement with SPH simulations and theoretical expectations. For measurements across two neighbouring redshift bins of width $Δz$ each (obtained by 'splitting' skewers of length $2Δz$ in half), the lognormal model predicts an anti-correlation for FPS $\times$ FPS and a positive correlation for $\bar F\times$ FPS, caused by long wavelength modes. This is in contrast to SPH simulations which predict a negligible magnitude for cross-redshift correlations derived from such `split' skewers, and we discuss possible reasons for this difference. Finally, we perform a preliminary test of the impact of neglecting long wavelength modes on parameter inference, finding that whereas the correlation structure of neighbouring redshift bins has relatively little impact, the absence of long wavelength modes in the model can lead to $\gtrsim2-σ$ biases in the inference of astrophysical parameters. Our results motivate a more careful treatment of long wavelength modes in analyses that rely on the small scale Ly$α$ forest for parameter inference. △ Less

Submitted 12 March, 2024; v1 submitted 25 October, 2023; originally announced October 2023.

Comments: 18 pages, 6 figures; v2: minor changes to match version accepted in JCAP

A modified lognormal approximation of the Lyman-$α$ forest: comparison with full hydrodynamic simulations at $2\leq z\leq 2.7$

Authors: Bhaskar Arya, Tirthankar Roy Choudhury, Aseem Paranjape, Prakash Gaikwad

Abstract: Observations of the Lyman-$α$ forest in distant quasar spectra with upcoming surveys are expected to provide significantly larger and higher-quality datasets. To interpret these datasets, it is imperative to develop efficient simulations. One such approach is based on the assumption that baryonic densities in the intergalactic medium (IGM) follow a lognormal distribution. We extend our earlier wor… ▽ More Observations of the Lyman-$α$ forest in distant quasar spectra with upcoming surveys are expected to provide significantly larger and higher-quality datasets. To interpret these datasets, it is imperative to develop efficient simulations. One such approach is based on the assumption that baryonic densities in the intergalactic medium (IGM) follow a lognormal distribution. We extend our earlier work to assess the robustness of the lognormal model of the Lyman-$α$ forest in recovering the parameters characterizing IGM state, namely, the mean-density IGM temperature ($T_0$), the slope of the temperature-density relation ($γ$), and the hydrogen photoionization rate ($Γ_{12}$), by comparing with high-resolution Sherwood SPH simulations across the redshift range $2 \leq z \leq 2.7$. These parameters are estimated through a Markov Chain Monte Carlo technique, using the mean and power spectrum of the transmitted flux. We find that the usual lognormal distribution of IGM densities cannot recover the parameters of the SPH simulations. This limitation arises from the fact that the SPH baryonic density distribution cannot be described by a simple lognormal form. To address this, we extend the model by scaling the linear density contrast by a parameter $ν$. While the resulting baryonic density is still lognormal, the additional parameter gives us extra freedom in setting the variance of density fluctuations. With this extension, values of $T_0$ and $γ$ implied in the SPH simulations are recovered at $\sim 1-σ$ ($\lesssim$ 10%) of the median (best-fit) values for most redshifts bins. However, this extended lognormal model cannot recover $Γ_{12}$ reliably, with the best-fit value discrepant by $\gtrsim 3-σ$ for $z > 2.2$. Despite this limitation in the recovery of $Γ_{12}$, we argue that the model remains useful for constraining cosmological parameters. △ Less

Submitted 12 March, 2024; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 20 pages, 8 figures; v2: minor changes to match version accepted in JCAP

Lost in the Middle: How Language Models Use Long Contexts

Authors: Nelson F. Liu, Kevin Lin, John Hewitt, Ashwin Paranjape, Michele Bevilacqua, Fabio Petroni, Percy Liang

Abstract: While recent language models have the ability to take long contexts as input, relatively little is known about how well they use longer context. We analyze the performance of language models on two tasks that require identifying relevant information in their input contexts: multi-document question answering and key-value retrieval. We find that performance can degrade significantly when changing t… ▽ More While recent language models have the ability to take long contexts as input, relatively little is known about how well they use longer context. We analyze the performance of language models on two tasks that require identifying relevant information in their input contexts: multi-document question answering and key-value retrieval. We find that performance can degrade significantly when changing the position of relevant information, indicating that current language models do not robustly make use of information in long input contexts. In particular, we observe that performance is often highest when relevant information occurs at the beginning or end of the input context, and significantly degrades when models must access relevant information in the middle of long contexts, even for explicitly long-context models. Our analysis provides a better understanding of how language models use their input context and provides new evaluation protocols for future long-context language models. △ Less

Submitted 20 November, 2023; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: 18 pages, 16 figures. Accepted for publication in Transactions of the Association for Computational Linguistics (TACL), 2023

A fast method of reionization parameter space exploration using GPR trained SCRIPT

Authors: Barun Maity, Aseem Paranjape, Tirthankar Roy Choudhury

Abstract: Efficient exploration of parameter spaces is crucial to extract physical information about the Epoch of Reionization from various observational probes. To this end, we propose a fast technique based on Gaussian Process Regression (GPR) training applied to a semi-numerical photon-conserving reionization model, SCRIPT. Our approach takes advantage of the numerical convergence properties of SCRIPT an… ▽ More Efficient exploration of parameter spaces is crucial to extract physical information about the Epoch of Reionization from various observational probes. To this end, we propose a fast technique based on Gaussian Process Regression (GPR) training applied to a semi-numerical photon-conserving reionization model, SCRIPT. Our approach takes advantage of the numerical convergence properties of SCRIPT and constructs a training set based on low-cost, coarse-resolution simulations. A likelihood emulator is then trained using this set to produce results in approximately two orders of magnitude less computational time than a full MCMC run, while still generating reasonable 68% and 95% confidence contours. Furthermore, we conduct a forecasting study using simulated data to demonstrate the applicability of this technique. This method is particularly useful when full MCMC analysis is not feasible due to expensive likelihood computations. △ Less

Submitted 28 September, 2023; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS

Impact of tidal environment on galaxy clustering in GAMA

Authors: Shadab Alam, Aseem Paranjape, John A. Peacock

Abstract: We constrain models of the galaxy distribution in the cosmic web using data from the Galaxy and Mass Assembly (GAMA) survey. We model the redshift-space behaviour of the 2-point correlation function (2pcf) and the recently proposed Voronoi volume function (VVF) -- which includes information beyond 2-point statistics. We extend the standard halo model using extra satellite degrees of freedom and tw… ▽ More We constrain models of the galaxy distribution in the cosmic web using data from the Galaxy and Mass Assembly (GAMA) survey. We model the redshift-space behaviour of the 2-point correlation function (2pcf) and the recently proposed Voronoi volume function (VVF) -- which includes information beyond 2-point statistics. We extend the standard halo model using extra satellite degrees of freedom and two assembly bias parameters, $α_{\rm cen}$ and $α_{\rm sat}$, which respectively correlate the occupation numbers of central and satellite galaxies with their host halo's tidal environment. We measure $α_{\rm sat}=1.44^{+0.25}_{-0.43}$ and $α_{\rm cen}=-0.79^{+0.29}_{-0.11}$ using a combination of 2pcf and VVF measurements, representing a detection of assembly bias at the 3.3$σ$ (2.4$σ$) significance level for satellite (central) galaxies. This result remains robust to possible anisotropies in the halo-centric distribution of satellites as well as technicalities of estimating the data covariance. We show that the growth rate ($fσ_8$) deduced using models with assembly bias is about 7\% (i.e. $1.5σ$) lower than if assembly bias is ignored. When projected onto the $Ω_m$-$σ_8$ plane, the model constraints without assembly bias overlap with Planck expectations, while allowing assembly bias introduces significant tension with Planck, preferring either a lower $Ω_m$ or a lower $σ_8$. Finally, we find that the all-galaxy weak lensing signal is unaffected by assembly bias, but the central and satellite sub-populations individually show significantly different signals in the presence of assembly bias. Our results illustrate the importance of accurately modelling galaxy formation for cosmological inference from future surveys. △ Less

Submitted 5 December, 2023; v1 submitted 2 May, 2023; originally announced May 2023.

Comments: 18 pages, 14 figures, 3 tables, minor changes in text and one additional table matches the accepted version

Model-agnostic cosmological constraints from the baryon acoustic oscillation feature in redshift space

Authors: Aseem Paranjape, Ravi K. Sheth

Abstract: We develop a framework for self-consistently extracting cosmological information from the clustering of tracers in redshift space, $\textit{without}$ relying on model-dependent templates to describe the baryon acoustic oscillation (BAO) feature. Our approach uses the recently proposed Laguerre reconstruction technique for the BAO feature and its linear point $r_{\rm LP}$, and substantially extends… ▽ More We develop a framework for self-consistently extracting cosmological information from the clustering of tracers in redshift space, $\textit{without}$ relying on model-dependent templates to describe the baryon acoustic oscillation (BAO) feature. Our approach uses the recently proposed Laguerre reconstruction technique for the BAO feature and its linear point $r_{\rm LP}$, and substantially extends it to simultaneously model the multipoles $\ell=0,2,4$ of the anisotropic galaxy 2-point correlation function (2pcf). The approach is `model-agnostic': it assumes that the non-linear growth of structure smears the BAO feature by an approximately Gaussian kernel with a smearing scale $σ_{\rm v}$, but does not assume any fiducial cosmology for describing the shape of the feature itself. Using mock observations for two realistic survey configurations assuming $Λ$ cold dark matter ($Λ$CDM), combined with Bayesian parameter inference, we show that the linear point $r_{\rm LP}$ and smearing scale $σ_{\rm v}$ can be accurately recovered by our method in both existing and upcoming surveys. The precision of the recovery of $r_{\rm LP}$ is always better than $1\%$, while $σ_{\rm v}$ can be recovered with $\lesssim10\%$ uncertainty provided the linear galaxy bias $b$ is separately constrained, e.g., using weak lensing observations. Our method is also sensitive to the linear growth rate $f$, albeit with larger uncertainties and systematic errors, especially for upcoming surveys such as DESI. We discuss how our model can be modified to improve the recovery of $f$, such that the resulting constraints on $\{f,σ_{\rm v},r_{\rm LP}\}$ can potentially be used as a test of cosmological models including and beyond $Λ$CDM. △ Less

Submitted 8 September, 2023; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 18 pages, 6 figures; v2: minor changes, accepted in MNRAS

A test of invariance of dark matter halo surface density using multiwavelength mock galaxy catalogues

Authors: Gopika K., Shantanu Desai, Aseem Paranjape

Abstract: A large number of observations have shown that the dark matter halo surface density, given by the product of halo core radius and core density is nearly constant for a diverse suite of galaxies. Although this invariance of the halo surface density is violated at galaxy cluster and group scales, it is still an open question on whether the aforementioned constancy on galactic scales can be explained… ▽ More A large number of observations have shown that the dark matter halo surface density, given by the product of halo core radius and core density is nearly constant for a diverse suite of galaxies. Although this invariance of the halo surface density is violated at galaxy cluster and group scales, it is still an open question on whether the aforementioned constancy on galactic scales can be explained within $Λ$CDM. For this purpose, we probe the variation of halo surface density as a function of mass using multi-wavelength mock galaxy catalogs from $Λ$CDM simulations, where the adiabatic contraction of dark matter halos in the presence of baryons has been taken into account. We find that these baryonified $Λ$CDM halos were best fitted with a generalized-NFW profile, and the halo surface density from these halos has a degeneracy with respect to both the halo mass and the virial concentration. We find that the correlation with mass when averaged over concentration is consistent with a constant halo surface density. However, a power-law dependence as a function of halo mass also cannot be ruled out. △ Less

Submitted 26 May, 2023; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: 10 pages, 11 figures. Accepted for publications in MNRAS

Loss of Distributed Coverage Using Lazy Agents Operating Under Discrete, Local, Event-Triggered Communication

Authors: Edward Vickery, Aditya A. Paranjape

Abstract: Continuous surveillance of a spatial region using distributed robots and sensors is a well-studied application in the area of multi-agent systems. This paper investigates a practically-relevant scenario where robotic sensors are introduced asynchronously and inter-robot communication is discrete, event-driven, local and asynchronous. Furthermore, we work with lazy robots; i.e., the robots seek to… ▽ More Continuous surveillance of a spatial region using distributed robots and sensors is a well-studied application in the area of multi-agent systems. This paper investigates a practically-relevant scenario where robotic sensors are introduced asynchronously and inter-robot communication is discrete, event-driven, local and asynchronous. Furthermore, we work with lazy robots; i.e., the robots seek to minimize their area of responsibility by equipartitioning the domain to be covered. We adapt a well-known algorithm which is practicable and known to generally work well for coverage problems. For a specially chosen geometry of the spatial domain, we show that there exists a non-trivial sequence of inter-robot communication events which leads to an instantaneous loss of coverage when the number of robots exceeds a certain threshold. The same sequence of events preserves coverage and, further, leads to an equipartition of the domain when the number of robots is smaller than the threshold. This result demonstrates that coverage guarantees for a given algorithm might be sensitive to the number of robots and, therefore, may not scale in obvious ways. It also suggests that when such algorithms are to be verified and validated prior to field deployment, the number of robots or sensors used in test scenarios should match that deployed on the field. △ Less

Submitted 17 February, 2023; originally announced February 2023.

Comments: 14 pages, 4 figures, extended abstract in AAMAS 2023

MSC Class: 93A16; 68T40

Evaluating Human-Language Model Interaction

Authors: Mina Lee, Megha Srivastava, Amelia Hardy, John Thickstun, Esin Durmus, Ashwin Paranjape, Ines Gerard-Ursin, Xiang Lisa Li, Faisal Ladhak, Frieda Rong, Rose E. Wang, Minae Kwon, Joon Sung Park, Hancheng Cao, Tony Lee, Rishi Bommasani, Michael Bernstein, Percy Liang

Abstract: Many real-world applications of language models (LMs), such as writing assistance and code autocomplete, involve human-LM interaction. However, most benchmarks are non-interactive in that a model produces output without human involvement. To evaluate human-LM interaction, we develop a new framework, Human-AI Language-based Interaction Evaluation (HALIE), that defines the components of interactive… ▽ More Many real-world applications of language models (LMs), such as writing assistance and code autocomplete, involve human-LM interaction. However, most benchmarks are non-interactive in that a model produces output without human involvement. To evaluate human-LM interaction, we develop a new framework, Human-AI Language-based Interaction Evaluation (HALIE), that defines the components of interactive systems and dimensions to consider when designing evaluation metrics. Compared to standard, non-interactive evaluation, HALIE captures (i) the interactive process, not only the final output; (ii) the first-person subjective experience, not just a third-party assessment; and (iii) notions of preference beyond quality (e.g., enjoyment and ownership). We then design five tasks to cover different forms of interaction: social dialogue, question answering, crossword puzzles, summarization, and metaphor generation. With four state-of-the-art LMs (three variants of OpenAI's GPT-3 and AI21 Labs' Jurassic-1), we find that better non-interactive performance does not always translate to better human-LM interaction. In particular, we highlight three cases where the results from non-interactive and interactive metrics diverge and underscore the importance of human-LM interaction for LM evaluation. △ Less

Submitted 5 January, 2024; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: Authored by the Center for Research on Foundation Models (CRFM) at the Stanford Institute for Human-Centered Artificial Intelligence (HAI)

Exploring the origin of stronger survival of polarized vortex beams through scattering media

Authors: Atharva Paranjape, Shyamal Guchhait, Athira B S, Nirmalya Ghosh

Abstract: Laguerre-Gaussian (LG) beams carrying orbital angular momentum (OAM) have shown promise in deep tissue imaging, medical diagnostics, and optical communication due to their robust propagation properties through scattering media. However, an exact model that provides a mechanism for the enhanced scattering properties of LG beams over Gaussian beams has not been established till date. Here, we examin… ▽ More Laguerre-Gaussian (LG) beams carrying orbital angular momentum (OAM) have shown promise in deep tissue imaging, medical diagnostics, and optical communication due to their robust propagation properties through scattering media. However, an exact model that provides a mechanism for the enhanced scattering properties of LG beams over Gaussian beams has not been established till date. Here, we examine this issue by studying the propagation of polarized vortex beams transmitted through tissue-like turbid scattering media. We demonstrate that the intensity profile has a much more profound effect on depolarization than the phase profile for LG beams. Our results indicate that the observed stronger propagation for the higher order LG beams is due to a higher anisotropy factor g as seen by the incident beam. This insight is expected to contribute towards building a complete picture of light transport in the presence of scattering, as well as guide optimization of the intensity and polarization structure of light for use in biomedical applications. △ Less

Submitted 21 December, 2022; v1 submitted 14 December, 2022; originally announced December 2022.

Bayesian evidence comparison for distance scale estimates

Authors: Aseem Paranjape, Ravi K. Sheth

Abstract: Constraints on cosmological parameters are often distilled from sky surveys by fitting templates to summary statistics of the data that are motivated by a fiducial cosmological model. However, recent work has shown how to estimate the distance scale using templates that are more generic: the basis functions used are not explicitly tied to any one cosmological model. We describe a Bayesian framewor… ▽ More Constraints on cosmological parameters are often distilled from sky surveys by fitting templates to summary statistics of the data that are motivated by a fiducial cosmological model. However, recent work has shown how to estimate the distance scale using templates that are more generic: the basis functions used are not explicitly tied to any one cosmological model. We describe a Bayesian framework for (i) determining how many basis functions to use and (ii) comparing one basis set with another. Our formulation provides intuition into how (a) one's degree of belief in different basis sets, (b) the fact that the choice of priors depends on basis set, and (c) the data set itself, together determine the derived constraints. We illustrate our framework using measurements in simulated datasets before applying it to real data. △ Less

Submitted 14 October, 2022; v1 submitted 1 September, 2022; originally announced September 2022.

Comments: 10 pages, 7 figures; v2: minor clarifications, fixed some typos, accepted in MNRAS

When can I Speak? Predicting initiation points for spoken dialogue agents

Authors: Siyan Li, Ashwin Paranjape, Christopher D. Manning

Abstract: Current spoken dialogue systems initiate their turns after a long period of silence (700-1000ms), which leads to little real-time feedback, sluggish responses, and an overall stilted conversational flow. Humans typically respond within 200ms and successfully predicting initiation points in advance would allow spoken dialogue agents to do the same. In this work, we predict the lead-time to initiati… ▽ More Current spoken dialogue systems initiate their turns after a long period of silence (700-1000ms), which leads to little real-time feedback, sluggish responses, and an overall stilted conversational flow. Humans typically respond within 200ms and successfully predicting initiation points in advance would allow spoken dialogue agents to do the same. In this work, we predict the lead-time to initiation using prosodic features from a pre-trained speech representation model (wav2vec 1.0) operating on user audio and word features from a pre-trained language model (GPT-2) operating on incremental transcriptions. To evaluate errors, we propose two metrics w.r.t. predicted and true lead times. We train and evaluate the models on the Switchboard Corpus and find that our method outperforms features from prior work on both metrics and vastly outperforms the common approach of waiting for 700ms of silence. △ Less

Submitted 7 August, 2022; originally announced August 2022.

Comments: SIGDIAL 2022

Neural Generation Meets Real People: Building a Social, Informative Open-Domain Dialogue Agent

Authors: Ethan A. Chi, Ashwin Paranjape, Abigail See, Caleb Chiam, Trenton Chang, Kathleen Kenealy, Swee Kiat Lim, Amelia Hardy, Chetanya Rastogi, Haojun Li, Alexander Iyabor, Yutong He, Hari Sowrirajan, Peng Qi, Kaushik Ram Sadagopan, Nguyet Minh Phu, Dilara Soylu, Jillian Tang, Avanika Narayan, Giovanni Campagna, Christopher D. Manning

Abstract: We present Chirpy Cardinal, an open-domain social chatbot. Aiming to be both informative and conversational, our bot chats with users in an authentic, emotionally intelligent way. By integrating controlled neural generation with scaffolded, hand-written dialogue, we let both the user and bot take turns driving the conversation, producing an engaging and socially fluent experience. Deployed in the… ▽ More We present Chirpy Cardinal, an open-domain social chatbot. Aiming to be both informative and conversational, our bot chats with users in an authentic, emotionally intelligent way. By integrating controlled neural generation with scaffolded, hand-written dialogue, we let both the user and bot take turns driving the conversation, producing an engaging and socially fluent experience. Deployed in the fourth iteration of the Alexa Prize Socialbot Grand Challenge, Chirpy Cardinal handled thousands of conversations per day, placing second out of nine bots with an average user rating of 3.58/5. △ Less

Submitted 16 January, 2023; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: SIGDIAL '22

Lognormal semi-numerical simulations of the Lyman-$α$ forest: comparison with full hydrodynamic simulations

Authors: Bhaskar Arya, Tirthankar Roy Choudhury, Aseem Paranjape, Prakash Gaikwad

Abstract: Observations of the Lyman-$α$ (Ly$α$) forest in spectra of distant quasars enable us to probe the matter power spectrum at relatively small scales. With several upcoming surveys, it is expected that there will be a many-fold increase in the quantity and quality of data, and hence it is important to develop efficient simulations to forward model these data sets. One such semi-numerical method is ba… ▽ More Observations of the Lyman-$α$ (Ly$α$) forest in spectra of distant quasars enable us to probe the matter power spectrum at relatively small scales. With several upcoming surveys, it is expected that there will be a many-fold increase in the quantity and quality of data, and hence it is important to develop efficient simulations to forward model these data sets. One such semi-numerical method is based on the assumption that the baryonic densities in the intergalactic medium (IGM) follow a lognormal distribution. In this work, we test the robustness of the lognormal model of the Ly$α$ forest in recovering a set of IGM parameters by comparing with high-resolution Sherwood SPH simulations. We study the recovery of the parameters $T_0$ (temperature of the mean-density IGM), $γ$ (slope of the temperature-density relation) and $Γ_{12}$ (hydrogen photoionization rate) at $z \sim 2.5$ using a Markov Chain Monte Carlo (MCMC) technique for parameter estimation. Using three flux statistics, the probability distribution, the mean flux and the power spectrum, values of all three parameters, $T_0$, $γ$ and $Γ_{12}$ implied in the SPH simulations are recovered within $1 - σ$ ($\sim$ 9, 4 and 1% respectively) of the median (best-fit) values. We verify the validity of our results at different baryon smoothing filter, SNR, box size & resolution, and data seed and confirm that the lognormal model can be used as an efficient tool for modelling the Ly$α$ transmitted flux at $z \sim 2.5$. △ Less

Submitted 1 February, 2023; v1 submitted 16 June, 2022; originally announced June 2022.

Comments: 16 pages, 11 figures; v2: revised version accepted in MNRAS

The quasi-adiabatic relaxation of haloes in the IllustrisTNG and EAGLE cosmological simulations

Authors: Premvijay Velmani, Aseem Paranjape

Abstract: The dark matter content of a gravitationally bound halo is known to be affected by the galaxy and gas it hosts. We characterise this response for haloes spanning over four orders of magnitude in mass in the hydrodynamical simulation suites IllustrisTNG and EAGLE. We present simple fitting functions in the spherically averaged quasi-adiabatic relaxation framework that accurately capture the dark ma… ▽ More The dark matter content of a gravitationally bound halo is known to be affected by the galaxy and gas it hosts. We characterise this response for haloes spanning over four orders of magnitude in mass in the hydrodynamical simulation suites IllustrisTNG and EAGLE. We present simple fitting functions in the spherically averaged quasi-adiabatic relaxation framework that accurately capture the dark matter response over the full range of halo mass and halo-centric distance we explore. We show that commonly employed schemes, which consider the relative change in radius $r_f/r_i-1$ of a spherical dark matter shell to be a function of only the relative change in its mass $M_i/M_f-1$, do not accurately describe the measured response of most haloes in IllustrisTNG and EAGLE. Rather, $r_f/r_i$ additionally explicitly depends upon halo-centric distance $r_f/R_{\rm vir}$ for haloes with virial radius $R_{\rm vir}$, being very similar between IllustrisTNG and EAGLE and across halo mass. We also account for a previously unmodelled effect, likely driven by feedback-related outflows, in which shells having $r_f/r_i\simeq1$ (i.e., no relaxation) have $M_i/M_f$ significantly different from unity. Our results are immediately applicable to a number of semi-analytical tools for modelling galactic and large-scale structure. We also study the dependence of this response on several halo and galaxy properties beyond total mass, finding that it is primarily related to halo concentration and star formation rate. We discuss possible extensions of these results to build a deeper physical understanding of the small-scale connection between dark matter and baryons. △ Less

Submitted 25 January, 2023; v1 submitted 15 June, 2022; originally announced June 2022.

Comments: 21 pages, 17 figures; v2: added discussion, accepted in MNRAS

A simulated annealing approach to parameter inference with expensive likelihoods

Authors: Aseem Paranjape

Abstract: We present a new approach to parameter inference targeted on generic situations where the evaluation of the likelihood $\mathcal{L}$ (i.e., the probability to observe the data given a fixed model configuration) is numerically expensive. Inspired by ideas underlying simulated annealing, the method first evaluates $χ^2=-2\ln\mathcal{L}$ on a sparse sequence of Latin hypercubes of increasing density… ▽ More We present a new approach to parameter inference targeted on generic situations where the evaluation of the likelihood $\mathcal{L}$ (i.e., the probability to observe the data given a fixed model configuration) is numerically expensive. Inspired by ideas underlying simulated annealing, the method first evaluates $χ^2=-2\ln\mathcal{L}$ on a sparse sequence of Latin hypercubes of increasing density in parameter (eigen)space. The semi-stochastic choice of sampling points accounts for anisotropic gradients of $χ^2$ and rapidly zooms in on the minimum of $χ^2$. The sampled $χ^2$ values are then used to train an interpolator which is further used in a standard Markov Chain Monte Carlo (MCMC) algorithm to inexpensively explore the parameter space with high density, similarly to emulator-based approaches now popular in cosmological studies. Comparisons with example linear and non-linear problems show gains in the number of likelihood evaluations of factors of 10 to 100 or more, as compared to standard MCMC algorithms. As a specific implementation, we publicly release the code PICASA: Parameter Inference using Cobaya with Anisotropic Simulated Annealing, which combines the minimizer (of a user-defined $χ^2$) with Gaussian Process Regression for training the interpolator and a subsequent MCMC implementation using the COBAYA framework. Being agnostic to the nature of the observable data and the theoretical model, our implementation is potentially useful for a number of emerging problems in cosmology, astrophysics and beyond. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 14 pages, 6 figures; submitted to MNRAS; comments welcome! PICASA is publicly available at https://bitbucket.org/aparanjape/picasa/

You Only Need One Model for Open-domain Question Answering

Authors: Haejun Lee, Akhil Kedia, Jongwon Lee, Ashwin Paranjape, Christopher D. Manning, Kyoung-Gu Woo

Abstract: Recent approaches to Open-domain Question Answering refer to an external knowledge base using a retriever model, optionally rerank passages with a separate reranker model and generate an answer using another reader model. Despite performing related tasks, the models have separate parameters and are weakly-coupled during training. We propose casting the retriever and the reranker as internal passag… ▽ More Recent approaches to Open-domain Question Answering refer to an external knowledge base using a retriever model, optionally rerank passages with a separate reranker model and generate an answer using another reader model. Despite performing related tasks, the models have separate parameters and are weakly-coupled during training. We propose casting the retriever and the reranker as internal passage-wise attention mechanisms applied sequentially within the transformer architecture and feeding computed representations to the reader, with the hidden representations progressively refined at each stage. This allows us to use a single question answering model trained end-to-end, which is a more efficient use of model capacity and also leads to better gradient flow. We present a pre-training method to effectively train this architecture and evaluate our model on the Natural Questions and TriviaQA open datasets. For a fixed parameter budget, our model outperforms the previous state-of-the-art model by 1.0 and 0.7 exact match scores. △ Less

Submitted 28 October, 2022; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: EMNLP 2022 (main)

A simulation driven optimization algorithm for scheduling sorting center operations

Authors: Supratim Ghosh, Aritra Pal, Prashant Kumar, Ankush Ojha, Aditya Paranjape, Souvik Barat, Harshad Khadilkar

Abstract: Parcel sorting operations in logistics enterprises aim to achieve a high throughput of parcels through sorting centers. These sorting centers are composed of large circular conveyor belts on which incoming parcels are placed, with multiple arms known as chutes for sorting the parcels by destination, followed by packing into roller cages and loading onto outbound trucks. Modern sorting systems need… ▽ More Parcel sorting operations in logistics enterprises aim to achieve a high throughput of parcels through sorting centers. These sorting centers are composed of large circular conveyor belts on which incoming parcels are placed, with multiple arms known as chutes for sorting the parcels by destination, followed by packing into roller cages and loading onto outbound trucks. Modern sorting systems need to complement their hardware innovations with sophisticated algorithms and software to map destinations and workforce to specific chutes. While state of the art systems operate with fixed mappings, we propose an optimization approach that runs before every shift, and uses real-time forecast of destination demand and labor availability in order to maximize throughput. We use simulation to improve the performance and robustness of the optimization solution to stochasticity in the environment, through closed-loop tuning of the optimization parameters. △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: 12 pages, Winter Simulation Conference 2021

The phenomenology of the external field effect in cold dark matter models

Authors: Aseem Paranjape, Ravi K. Sheth

Abstract: In general relativity (GR), the internal dynamics of a self-gravitating system under free-fall in an external gravitational field should not depend on the external field strength. Recent work has claimed a statistical detection of an `external field effect' (EFE) using galaxy rotation curve data. We show that large uncertainties in rotation curve analyses and inaccuracies in published simulation-b… ▽ More In general relativity (GR), the internal dynamics of a self-gravitating system under free-fall in an external gravitational field should not depend on the external field strength. Recent work has claimed a statistical detection of an `external field effect' (EFE) using galaxy rotation curve data. We show that large uncertainties in rotation curve analyses and inaccuracies in published simulation-based external field estimates compromise the significance of the claimed EFE detection. We further show analytically that a qualitatively similar statistical signal is, in fact, expected in a $Λ$-cold dark matter ($Λ$CDM) universe without any violation of the strong equivalence principle. Rather, such a signal arises simply because of the inherent correlations between galaxy clustering strength and intrinsic galaxy properties. We explicitly demonstrate the effect in a baryonified mock catalog of a $Λ$CDM universe. Although the detection of an EFE-like signal is not, by itself, evidence for physics beyond GR, our work shows that the $\textit{sign}$ of the EFE-like correlation between the external field strength and the shape of the radial acceleration relation can be used to probe new physics: e.g., in MOND, the predicted sign is opposite to that in our $Λ$CDM mocks. △ Less

Submitted 20 September, 2022; v1 submitted 30 November, 2021; originally announced December 2021.

Comments: 11 pages, 6 figures; v2: minor changes and added discussion, accepted in MNRAS

Hindsight: Posterior-guided training of retrievers for improved open-ended generation

Authors: Ashwin Paranjape, Omar Khattab, Christopher Potts, Matei Zaharia, Christopher D. Manning

Abstract: Many text generation systems benefit from using a retriever to retrieve passages from a textual knowledge corpus (e.g., Wikipedia) which are then provided as additional context to the generator. For open-ended generation tasks (like generating informative utterances in conversations) many varied passages may be equally relevant and we find that existing methods that jointly train the retriever and… ▽ More Many text generation systems benefit from using a retriever to retrieve passages from a textual knowledge corpus (e.g., Wikipedia) which are then provided as additional context to the generator. For open-ended generation tasks (like generating informative utterances in conversations) many varied passages may be equally relevant and we find that existing methods that jointly train the retriever and generator underperform: the retriever may not find relevant passages even amongst the top-10 and hence the generator may not learn a preference to ground its generated output in them. We propose using an additional guide retriever that is allowed to use the target output and "in hindsight" retrieve relevant passages during training. We model the guide retriever after the posterior distribution Q of passages given the input and the target output and train it jointly with the standard retriever and the generator by maximizing the evidence lower bound (ELBo) in expectation over Q. For informative conversations from the Wizard of Wikipedia dataset, with posterior-guided training, the retriever finds passages with higher relevance in the top-10 (23% relative improvement), the generator's responses are more grounded in the retrieved passage (19% relative improvement) and the end-to-end system produces better overall output (6.4% relative improvement). △ Less

Submitted 20 October, 2021; v1 submitted 14 October, 2021; originally announced October 2021.

The Life and Science of Thanu Padmanabhan

Authors: Jasjeet Singh Bagla, Krishnakanta Bhattacharya, Sumanta Chakraborty, Sunu Engineer, Valerio Faraoni, Sanved Kolekar, Dawood Kothawala, Kinjalk Lochan, Sujoy Modak, V. Parameswaran Nair, Aseem Paranjape, Krishnamohan Parattu, Sarada G. Rajeev, Bibhas Ranjan Majhi, Tirthankar Roy Choudhury, Mohammad Sami, Sudipta Sarkar, Sandipan Sengupta, T. R. Seshadri, S. Shankaranarayanan, Suprit Singh, Tejinder P. Singh, L. Sriramkumar, Urjit Yajnik

Abstract: Thanu Padmanabhan was a renowned Indian theoretical physicist known for his research in general relativity, cosmology, and quantum gravity. In an extraordinary career spanning forty-two years, he published more than three hundred research articles, wrote ten highly successful technical and popular books, and mentored nearly thirty graduate students and post-doctoral fellows. He is best known for h… ▽ More Thanu Padmanabhan was a renowned Indian theoretical physicist known for his research in general relativity, cosmology, and quantum gravity. In an extraordinary career spanning forty-two years, he published more than three hundred research articles, wrote ten highly successful technical and popular books, and mentored nearly thirty graduate students and post-doctoral fellows. He is best known for his deep work investigating gravitation as an emergent thermodynamic phenomenon. He was an outstanding teacher, and an indefatigable populariser of science, who travelled very widely to motivate and inspire young students. Paddy, as he was affectionately known, was also a close friend to his students and collaborators, treating them as part of his extended academic family. On September 17, 2021 Paddy passed away very unexpectedly, at the age of sixty-four and at the height of his research career, while serving as a Distinguished Professor at the Inter-University Centre for Astronomy and Astrophysics, Pune. His untimely demise has come as a shock to his family and friends and colleagues. In this article, several of them have come together to pay their tributes and share their fond memories of Paddy. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: In Memoriam of Prof. T. Padmanabhan (1957-2021)

The distribution of HI velocity profiles in a $Λ$CDM universe

Authors: Aseem Paranjape, R. Srianand, Tirthankar Roy Choudhury, Ravi K. Sheth

Abstract: We model the distribution of the observed profiles of 21 cm line emission from neutral hydrogen (HI) in central galaxies selected from a statistically representative mock catalog of the local Universe in the Lambda-cold dark matter framework. The distribution of these HI velocity profiles (specifically, their widths $W_{50}$) has been observationally constrained, but has not been systematically st… ▽ More We model the distribution of the observed profiles of 21 cm line emission from neutral hydrogen (HI) in central galaxies selected from a statistically representative mock catalog of the local Universe in the Lambda-cold dark matter framework. The distribution of these HI velocity profiles (specifically, their widths $W_{50}$) has been observationally constrained, but has not been systematically studied theoretically. Our model profiles derive from rotation curves of realistically baryonified haloes in an N-body simulation, including the quasi-adiabatic relaxation of the dark matter profile of each halo in response to its baryons. We study the predicted $W_{50}$ distribution using a realistic pipeline applied to noisy profiles extracted from our luminosity-complete mock catalog with an ALFALFA-like survey geometry and redshift selection. Our default mock is in good agreement with observed ALFALFA results for $W_{50}\gtrsim700$ km/s, being incomplete at lower widths due to the intrinsic threshold of $M_r\leq-19$. Variations around the default model show that the velocity width function at $W_{50}\gtrsim300$ km/s is most sensitive to a possible correlation between galaxy inclination and host concentration, followed by the physics of quasi-adiabatic relaxation. We also study the excess kurtosis of noiseless velocity profiles, obtaining a distribution which tightly correlates with $W_{50}$, with a shape and scatter that depend on the properties of the turbulent HI disk. Our results open the door towards using the shapes of HI velocity profiles as a novel statistical probe of the baryon-dark matter connection. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 22 pages, 14 figures, submitted to MNRAS

Transient and Asymptotic Properties of Robust Adaptive Controllers in the Presence of Non-Coercive Lyapunov Functions

Authors: Aditya A. Paranjape, Vivek Natarajan, Supratim Ghosh

Abstract: Adaptive control architectures often make use of Lyapunov functions to design adaptive laws. We are specifically interested in adaptive control methods, such as the well-known L1 adaptive architecture, which employ a parameter observer for this purpose. In such architectures, the observation error plays a critical role in determining analytical bounds on the tracking error as well as robustness. I… ▽ More Adaptive control architectures often make use of Lyapunov functions to design adaptive laws. We are specifically interested in adaptive control methods, such as the well-known L1 adaptive architecture, which employ a parameter observer for this purpose. In such architectures, the observation error plays a critical role in determining analytical bounds on the tracking error as well as robustness. In this paper, we show how the non-existence of coercive Lyapunov operators can impact the analytical bounds, and with it the performance and the robustness of such adaptive systems. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 13 pages,1 figure

MSC Class: 93C20; 93C25; 93C40; 93D05

Human-like informative conversations: Better acknowledgements using conditional mutual information

Authors: Ashwin Paranjape, Christopher D. Manning

Abstract: This work aims to build a dialogue agent that can weave new factual content into conversations as naturally as humans. We draw insights from linguistic principles of conversational analysis and annotate human-human conversations from the Switchboard Dialog Act Corpus to examine humans strategies for acknowledgement, transition, detail selection and presentation. When current chatbots (explicitly p… ▽ More This work aims to build a dialogue agent that can weave new factual content into conversations as naturally as humans. We draw insights from linguistic principles of conversational analysis and annotate human-human conversations from the Switchboard Dialog Act Corpus to examine humans strategies for acknowledgement, transition, detail selection and presentation. When current chatbots (explicitly provided with new factual content) introduce facts into a conversation, their generated responses do not acknowledge the prior turns. This is because models trained with two contexts - new factual content and conversational history - generate responses that are non-specific w.r.t. one of the contexts, typically the conversational history. We show that specificity w.r.t. conversational history is better captured by Pointwise Conditional Mutual Information ($\text{pcmi}_h$) than by the established use of Pointwise Mutual Information ($\text{pmi}$). Our proposed method, Fused-PCMI, trades off $\text{pmi}$ for $\text{pcmi}_h$ and is preferred by humans for overall quality over the Max-PMI baseline 60% of the time. Human evaluators also judge responses with higher $\text{pcmi}_h$ better at acknowledgement 74% of the time. The results demonstrate that systems mimicking human conversational traits (in this case acknowledgement) improve overall quality and more broadly illustrate the utility of linguistic principles in improving dialogue agents. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: NAACL 2021

The radial acceleration relation in a $Λ$CDM universe

Authors: Aseem Paranjape, Ravi K. Sheth

Abstract: We study the radial acceleration relation (RAR) between the total ($a_{\rm tot}$) and baryonic ($a_{\rm bary}$) centripetal acceleration profiles of central galaxies in the cold dark matter (CDM) paradigm. We analytically show that the RAR is intimately connected with the physics of the quasi-adiabatic relaxation of dark matter in the presence of baryons in deep potential wells. This cleanly demon… ▽ More We study the radial acceleration relation (RAR) between the total ($a_{\rm tot}$) and baryonic ($a_{\rm bary}$) centripetal acceleration profiles of central galaxies in the cold dark matter (CDM) paradigm. We analytically show that the RAR is intimately connected with the physics of the quasi-adiabatic relaxation of dark matter in the presence of baryons in deep potential wells. This cleanly demonstrates how the mean RAR and its scatter emerge in the low-acceleration regime ($10^{-12}\,{\rm m\,s}^{-2}\lesssim a_{\rm bary}\lesssim10^{-10}\,{\rm m\,s}^{-2}$) from an interplay between baryonic feedback processes and the distribution of CDM in dark halos. Our framework allows us to go further and study both higher and lower accelerations in detail, using analytical approximations and a realistic mock catalog of $\sim342,000$ low-redshift central galaxies with $M_r\leq-19$. We show that, while the RAR in the baryon-dominated, high-acceleration regime ($a_{\rm bary}\gtrsim10^{-10}\,{\rm m\,s}^{-2}$) is very sensitive to details of the relaxation physics, a simple `baryonification' prescription matching the relaxation results of hydrodynamical CDM simulations is remarkably successful in reproducing the observed RAR without any tuning. And in the (currently unobserved) ultra-low-acceleration regime ($a_{\rm bary}\lesssim 10^{-12}\,{\rm m\,s}^{-2}$), the RAR is sensitive to the abundance of diffuse gas in the halo outskirts, with our default model predicting a distinctive break from a simple power-law-like relation for HI-deficient, diffuse gas-rich centrals. Our mocks also show that the RAR provides more robust, testable predictions of the $Λ$CDM paradigm at galactic scales, with implications for alternative gravity theories, than the baryonic Tully-Fisher relation. △ Less

Submitted 23 July, 2021; v1 submitted 25 February, 2021; originally announced February 2021.

Comments: 20 pages, 16 figures; v2 - added discussion, clarifications and references, conclusions unchanged, accepted in MNRAS

Multi-wavelength mock galaxy catalogs of the low-redshift Universe

Authors: Aseem Paranjape, Tirthankar Roy Choudhury, Ravi K. Sheth

Abstract: We present a new suite of mock galaxy catalogs mimicking the low-redshift Universe, based on an updated halo occupation distribution (HOD) model and a scaling relation between optical properties and the neutral hydrogen (HI) content of galaxies. Our algorithm is constrained by observations of the luminosity function and luminosity- and colour-dependent clustering of SDSS galaxies, as well as the H… ▽ More We present a new suite of mock galaxy catalogs mimicking the low-redshift Universe, based on an updated halo occupation distribution (HOD) model and a scaling relation between optical properties and the neutral hydrogen (HI) content of galaxies. Our algorithm is constrained by observations of the luminosity function and luminosity- and colour-dependent clustering of SDSS galaxies, as well as the HI mass function and HI-dependent clustering of massive HI-selected galaxies in the ALFALFA survey. Mock central and satellite galaxies with realistic values of $r$-band luminosity, $g-r$ and $u-r$ colour, stellar mass and HI mass are populated in an $N$-body simulation, inheriting a number of properties of the density and tidal environment of their host halos. The host halo of each central galaxy is also `baryonified' with realistic spatial distributions of stars as well as hot and cold gas, along with the corresponding rotation curve. Our default HOD assumes that galaxy properties are a function of group halo mass alone, and can optionally include effects such as galactic conformity and colour-dependent galaxy assembly bias. The mocks predict the relation between the stellar mass and HI mass of massive HI galaxies, as well as the 2-point cross-correlation function of spatially co-located optical and HI-selected samples. They enable novel null tests for galaxy assembly bias, provide predictions for the HI velocity width function, and clarify the origin and universality of the radial acceleration relation in the $Λ$CDM framework. △ Less

Submitted 8 March, 2021; v1 submitted 7 January, 2021; originally announced January 2021.

Comments: 17 pages, 13 figures; v2: minor changes and clarifications, accepted in MNRAS

Mass accretion rates and multi-scale halo environment in cold and warm dark matter cosmologies

Authors: Payaswinee Dhoke, Aseem Paranjape

Abstract: We study the evolving environment dependence of mass accretion by dark haloes in simulations of cold and warm dark matter (CDM and WDM) cosmologies. The latter allows us to probe the nature of halo growth at scales below the WDM half-mode mass, which form an extreme regime of nonlinear collisionless dynamics and offer an excellent test-bed for ideas relating to hierarchical growth. As environmenta… ▽ More We study the evolving environment dependence of mass accretion by dark haloes in simulations of cold and warm dark matter (CDM and WDM) cosmologies. The latter allows us to probe the nature of halo growth at scales below the WDM half-mode mass, which form an extreme regime of nonlinear collisionless dynamics and offer an excellent test-bed for ideas relating to hierarchical growth. As environmental proxies, we use the local halo-centric matter density $δ$ and tidal anisotropy $α$, as well as large-scale halo bias $b_1$. Our analysis, while reproducing known trends for environment-dependent accretion in CDM, as well as the comparison between accretion in CDM and WDM, reveals several interesting new features. As expected from excursion set models, WDM haloes have higher specific accretion rates, dominated by the accretion of diffuse mass, as compared to CDM haloes. For low-mass WDM haloes, we find that the environment-dependence of both diffuse mass accretion as well as accretion by mergers is almost fully explained by $α$. For the other cases, $δ$ plays at least a comparable role. We detect, for the first time, a significant and evolving assembly bias due to diffuse mass accretion for low-mass CDM and WDM haloes (after excluding splashback objects), with a $z=0$ strength higher than with almost all known secondary variables and largely explained by $α$. Our results place constraints on semi-analytical merger tree algorithms, which in turn could affect the predictions of galaxy evolution models based on them. △ Less

Submitted 25 February, 2022; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: 16 pages, 11 figures, submitted to MNRAS

Journal ref: MNRAS 508 (2021) 852-867

Mock halo catalogs: assigning unresolved halo properties using correlations with local halo environment

Authors: Sujatha Ramakrishnan, Aseem Paranjape, Ravi K. Sheth

Abstract: Large-scale sky surveys require companion large volume simulated mock catalogs. To ensure precision cosmology studies are unbiased, the correlations in these mocks between galaxy properties and their large-scale environments must be realistic. Since galaxies are embedded in dark matter halos, an important first step is to include such correlations -- sometimes called assembly bias -- for dark matt… ▽ More Large-scale sky surveys require companion large volume simulated mock catalogs. To ensure precision cosmology studies are unbiased, the correlations in these mocks between galaxy properties and their large-scale environments must be realistic. Since galaxies are embedded in dark matter halos, an important first step is to include such correlations -- sometimes called assembly bias -- for dark matter halos. However, galaxy properties correlate with smaller scale physics in halos which large simulations struggle to resolve. We describe an algorithm which addresses and largely mitigates this problem. Our algorithm exploits the fact that halo assembly bias is unchanged as long as correlations between halo property $c$ and the intermediate-scale tidal environment $α$ are preserved. Therefore, knowledge of $α$ is sufficient to assign small-scale, otherwise unresolved properties to a halo in a way which preserves its large-scale assembly bias accurately. We demonstrate this explicitly for halo internal properties like formation history (concentration $c_{\rm 200b}$), shape $c/a$, dynamics $c_{v}/a_{v}$, velocity anisotropy $β$ and angular momentum (spin $λ$). Our algorithm increases a simulation's reach in halo mass and number density by an order of magnitude, with improvements in the bias signal as large as 45% for 30-particle halos, thus significantly reducing the cost of mocks for future weak lensing and redshift space distortion studies. △ Less

Submitted 22 February, 2021; v1 submitted 18 December, 2020; originally announced December 2020.

Comments: 13 pages, 10 figures, minor revisions, accepted in MNRAS

Sub-Optimality of a Dyadic Adaptive Control Architecture

Authors: Aditya A. Paranjape, Soon-Jo Chung

Abstract: The dyadic adaptive control architecture evolved as a solution to the problem of designing control laws for nonlinear systems with unmatched nonlinearities, disturbances and uncertainties. A salient feature of this framework is its ability to work with infinite as well as finite dimensional systems, and with a wide range of control and adaptive laws. In this paper, we consider the case where a con… ▽ More The dyadic adaptive control architecture evolved as a solution to the problem of designing control laws for nonlinear systems with unmatched nonlinearities, disturbances and uncertainties. A salient feature of this framework is its ability to work with infinite as well as finite dimensional systems, and with a wide range of control and adaptive laws. In this paper, we consider the case where a control law based on the linear quadratic regulator theory is employed for designing the control law. We benchmark the closed-loop system against standard linear quadratic control laws as well as those based on the state-dependent Riccati equation. We pose the problem of designing a part of the control law as a Nehari problem. We obtain analytical expressions for the bounds on the sub-optimality of the control law. △ Less

Submitted 18 October, 2020; originally announced October 2020.

Comments: 13 pages, 1 figure

MSC Class: 93C40; 93C15; 93C20; 93C25

Thermodynamics of Near BPS Black Holes in AdS$_4$ and AdS$_7$

Authors: Finn Larsen, and Shruti Paranjape

Abstract: We develop the thermodynamics of black holes in AdS$_4$ and AdS$_7$ near their BPS limit. In each setting we study the two distinct deformations orthogonal to the BPS surface as well as their nontrivial interplay with each other and with BPS properties. Our results illuminate recent microscopic calculations of the BPS entropy. We show that these microscopic computations can be leveraged to also de… ▽ More We develop the thermodynamics of black holes in AdS$_4$ and AdS$_7$ near their BPS limit. In each setting we study the two distinct deformations orthogonal to the BPS surface as well as their nontrivial interplay with each other and with BPS properties. Our results illuminate recent microscopic calculations of the BPS entropy. We show that these microscopic computations can be leveraged to also describe the near BPS regime, by generalizing the boundary conditions imposed on states. △ Less

Submitted 8 October, 2020; originally announced October 2020.

Comments: 45 pages, 1 figure

A sea of tides: star formation and the central-satellite dichotomy in a continuum of tidal environments

Authors: Jolanta Zjupa, Aseem Paranjape, Oliver Hahn, Rüdiger Pakmor

Abstract: The environment-dependent bimodality of the distribution of stellar mass ($M_\ast$) and specific star formation rate (sSFR) of galaxies, and its explanation in terms of the central-satellite dichotomy, form a cornerstone of our current understanding of galaxy evolution in the hierarchical structure formation paradigm. We revisit this framework in the IllustrisTNG simulation in the context of the m… ▽ More The environment-dependent bimodality of the distribution of stellar mass ($M_\ast$) and specific star formation rate (sSFR) of galaxies, and its explanation in terms of the central-satellite dichotomy, form a cornerstone of our current understanding of galaxy evolution in the hierarchical structure formation paradigm. We revisit this framework in the IllustrisTNG simulation in the context of the most extreme local tidal anisotropy $α_{\rm peak}$ experienced by each galaxy over cosmic time, which is an excellent proxy for environmental influence. We show that, while sharing a common monotonic $M_\ast$-$v_{\rm peak}$ relation, central, satellite and `splashback' galaxies define a hierarchy of increasing $α_{\rm peak}$. We also find that the sSFR of objects in small haloes unaffected by feedback from an active nucleus typically decreases with increasing $α_{\rm peak}$. Our results support an alternate viewpoint in which a galaxy can be identified by the value of $α_{\rm peak}$; i.e., rather than being placed on the central-satellite dichotomy, a galaxy is better classified by its location in a continuum of tidal environments. This conceptual shift can potentially yield a more robust understanding of galaxy evolution and the galaxy-dark matter connection, e.g., in accurately modelling subtle effects such as sSFR-induced secondary clustering. △ Less

Submitted 7 September, 2020; originally announced September 2020.

Comments: 5 pages, 3 figures, submitted to MNRAS, comments welcome!

Neural Generation Meets Real People: Towards Emotionally Engaging Mixed-Initiative Conversations

Authors: Ashwin Paranjape, Abigail See, Kathleen Kenealy, Haojun Li, Amelia Hardy, Peng Qi, Kaushik Ram Sadagopan, Nguyet Minh Phu, Dilara Soylu, Christopher D. Manning

Abstract: We present Chirpy Cardinal, an open-domain dialogue agent, as a research platform for the 2019 Alexa Prize competition. Building an open-domain socialbot that talks to real people is challenging - such a system must meet multiple user expectations such as broad world knowledge, conversational style, and emotional connection. Our socialbot engages users on their terms - prioritizing their interests… ▽ More We present Chirpy Cardinal, an open-domain dialogue agent, as a research platform for the 2019 Alexa Prize competition. Building an open-domain socialbot that talks to real people is challenging - such a system must meet multiple user expectations such as broad world knowledge, conversational style, and emotional connection. Our socialbot engages users on their terms - prioritizing their interests, feelings and autonomy. As a result, our socialbot provides a responsive, personalized user experience, capable of talking knowledgeably about a wide variety of topics, as well as chatting empathetically about ordinary life. Neural generation plays a key role in achieving these goals, providing the backbone for our conversational and emotional tone. At the end of the competition, Chirpy Cardinal progressed to the finals with an average rating of 3.6/5.0, a median conversation duration of 2 minutes 16 seconds, and a 90th percentile duration of over 12 minutes. △ Less

Submitted 5 September, 2020; v1 submitted 27 August, 2020; originally announced August 2020.

Comments: Published in 3rd Proceedings of Alexa Prize (Alexa Prize 2019)

Separate Universe calibration of the dependence of halo bias on cosmic web anisotropy

Authors: Sujatha Ramakrishnan, Aseem Paranjape

Abstract: We use the Separate Universe technique to calibrate the dependence of linear and quadratic halo bias $b_1$ and $b_2$ on the local cosmic web environment of dark matter haloes. We do this by measuring the response of halo abundances at fixed mass and cosmic web tidal anisotropy $α$ to an infinite wavelength initial perturbation. We augment our measurements with an analytical framework developed in… ▽ More We use the Separate Universe technique to calibrate the dependence of linear and quadratic halo bias $b_1$ and $b_2$ on the local cosmic web environment of dark matter haloes. We do this by measuring the response of halo abundances at fixed mass and cosmic web tidal anisotropy $α$ to an infinite wavelength initial perturbation. We augment our measurements with an analytical framework developed in earlier work which exploits the near-Lognormal shape of the distribution of $α$ and results in very high precision calibrations. We present convenient fitting functions for the dependence of $b_1$ and $b_2$ on $α$ over a wide range of halo mass for redshifts $0\leq z\leq1$. Our calibration of $b_2(α)$ is the first demonstration to date of the dependence of non-linear bias on the local web environment. Motivated by previous results which showed that $α$ is the primary indicator of halo assembly bias for a number of halo properties beyond halo mass, we then extend our analytical framework to accommodate the dependence of $b_1$ and $b_2$ on any such secondary property which has, or can be monotonically transformed to have, a Gaussian distribution. We demonstrate this technique for the specific case of halo concentration, finding good agreement with previous results. Our calibrations will be useful for a variety of halo model analyses focusing on galaxy assembly bias, as well as analytical forecasts of the potential for using $α$ as a segregating variable in multi-tracer analyses. △ Less

Submitted 26 September, 2020; v1 submitted 7 July, 2020; originally announced July 2020.

Comments: 16 pages, 11 figures; v2 - added discussion and references, results unchanged, accepted in MNRAS, calibration codes available at https://github.com/rsujatha/CAB

Analytical halo models of cosmic tidal fields

Authors: Aseem Paranjape

Abstract: The non-linear cosmic web environment of dark matter haloes plays a major role in shaping their growth and evolution, and potentially also affects the galaxies that reside in them. We develop an analytical (halo model) formalism to describe the tidal field of anisotropic halo-centric density distributions, as characterised by the halo-centric tidal tensor $\langle T_{ij} \rangle(<R)$ spherically a… ▽ More The non-linear cosmic web environment of dark matter haloes plays a major role in shaping their growth and evolution, and potentially also affects the galaxies that reside in them. We develop an analytical (halo model) formalism to describe the tidal field of anisotropic halo-centric density distributions, as characterised by the halo-centric tidal tensor $\langle T_{ij} \rangle(<R)$ spherically averaged on scale $R\sim4R_{\rm vir}$ for haloes of virial radius $R_{\rm vir}$. We focus on axisymmetric anisotropies, which allows us to explore simple and intuitive toy models of (sub)halo configurations that exemplify some of the most interesting anisotropies in the cosmic web. We build our models around the spherical Navarro-Frenk-White (NFW) profile after describing it as a Gaussian mixture, which leads to almost fully analytical expressions for the `tidal anisotropy' scalar $α(<4R_{\rm vir})$ extracted from the tidal tensor. Our axisymmetric examples include (i) a spherical halo at the axis of a cylindrical filament, (ii) an off-centred satellite in a spherical host halo and (iii) an axisymmetric halo. Using these, we demonstrate several interesting results. For example, the tidal tensor at the axis of a pure cylindrical filament gives $α^{\rm (fil)}(<R)=1/2$ exactly, for any $R$. Also, $α(<4R_{\rm vir,sat})$ for a satellite of radius $R_{\rm vir,sat}$ as a function of its host-centric distance is a sensitive probe of dynamical mass loss of the satellite in its host environment. Finally, we discuss a number of potentially interesting extensions and applications of our formalism that can deepen our understanding of the multi-scale phenomenology of the cosmic web. △ Less

Submitted 4 February, 2021; v1 submitted 24 June, 2020; originally announced June 2020.

Comments: 18 pages, 8 figures; v2: minor changes, references added. Accepted in MNRAS. Code to implement the models along with example usage can be downloaded from https://bitbucket.org/aparanjape/gaussmixnfw

Three- and two-point spatial correlations of IGM at $z\sim 2$: Cloud based analysis using simulations

Authors: Soumak Maitra, Raghunathan Srianand, Prakash Gaikwad, Tirthankar Roy Choudhury, Aseem Paranjape, Patrick Petitjean

Abstract: Ly$α$ forest absorption spectra decomposed into multiple Voigt profile components (clouds) allow us to study clustering of intergalactic medium (IGM) as a function of HI column density ($N_{\rm HI}$). Here, we explore the transverse three-point correlation ($ζ$) of these Ly$α$ clouds using mock triplet spectra obtained from hydrodynamical simulations at $z \sim 2$ on scales of 1-5 $h^{-1}$cMpc. We… ▽ More Ly$α$ forest absorption spectra decomposed into multiple Voigt profile components (clouds) allow us to study clustering of intergalactic medium (IGM) as a function of HI column density ($N_{\rm HI}$). Here, we explore the transverse three-point correlation ($ζ$) of these Ly$α$ clouds using mock triplet spectra obtained from hydrodynamical simulations at $z \sim 2$ on scales of 1-5 $h^{-1}$cMpc. We find $ζ$ to depend strongly on $N_{\rm HI}$ and scale and weakly on angle ($θ$) of the triplet configuration. We show that the "hierarchical ansatz" is applicable for scales $\ge~ 3h^{-1}$cMpc, and obtain a median reduced three-point correlation (Q) in the range 0.2-0.7. We show, $ζ$ is influenced strongly by the thermal and ionization state of the gas. As found in the case of galaxies, the influence of physical parameters on Q is weaker compared to that of $ζ$. We show difference in $ζ$ and Q between different simulations are minimized if we use appropriate $N_{\rm HI}$ cut-offs corresponding to a given baryon over-density ($Δ$) using the measured $N_{\rm HI}~vs~Δ$ relationship obtained from individual simulations. Additionally, we see the effect of pressure broadening on $ζ$ in a model with artificially boosted heating rates. However, for models with realistic thermal and ionization histories the effect of pressure broadening on $ζ$ is weak and sub-dominant compared to other local effects. We find strong redshift evolution shown by $ζ$, mainly originating from the redshift evolution of thermal and ionization state of the IGM. We discuss the observational requirements for the detection of three-point correlation, specifically, in small intervals of configuration parameters and redshift. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: 26 pages, 23 figures. Submitted to MNRAS

Studying the Lyman-$α$ optical depth fluctuations at $z \sim 5.5$ using fast semi-numerical methods

Authors: T. Roy Choudhury, Aseem Paranjape, Sarah E. I. Bosman

Abstract: We present a computationally efficient and fast semi-numerical technique for simulating the Lyman-$α$ (Ly$α$) absorption optical depth in presence of neutral hydrogen "islands" left over from reionization at redshifts $5 \lesssim z \lesssim 6$. The main inputs to the analysis are (i) a semi-numerical photon-conserving model of ionized regions during reionization (named SCRIPT) along with a prescri… ▽ More We present a computationally efficient and fast semi-numerical technique for simulating the Lyman-$α$ (Ly$α$) absorption optical depth in presence of neutral hydrogen "islands" left over from reionization at redshifts $5 \lesssim z \lesssim 6$. The main inputs to the analysis are (i) a semi-numerical photon-conserving model of ionized regions during reionization (named SCRIPT) along with a prescription for simulating the shadowing by neutral islands and (ii) the fluctuating Gunn-Peterson approximation to model the Ly$α$ absorption. Our model is then used for simulating the large-scale fluctuations in the effective optical depth as observed along sight lines towards high$-z$ quasars. Our model is fully described by five parameters. By setting two of them to default values and varying the other three, we obtain the constraints on reionization history at $5 \lesssim z \lesssim 6$ as allowed by the data. We confirm that reionization is \emph{not} complete before $z \sim 5.6$ at $\gtrsim 2σ$ confidence, with the exact confidence limits depending on how the non-detections of the flux in the data are treated. We also confirm that the completion of reionization can be as late as $z \sim 5.2$. With further improvements in the model and with more sight lines at $z \sim 6$, we can take advantage of the computational efficiency of our analysis to obtain more stringent constraints on the ionization fraction at the tail-end of reionization. △ Less

Submitted 17 December, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: Accepted in MNRAS. Matches the accepted version

Voronoi volume function: A new probe of cosmology and galaxy evolution

Authors: Aseem Paranjape, Shadab Alam

Abstract: We study the Voronoi volume function (VVF) -- the distribution of cell volumes (or inverse local number density) in the Voronoi tessellation of any set of cosmological tracers (galaxies/haloes). We show that the shape of the VVF of biased tracers responds sensitively to physical properties such as halo mass, large-scale environment, substructure and redshift-space effects, making this a hitherto u… ▽ More We study the Voronoi volume function (VVF) -- the distribution of cell volumes (or inverse local number density) in the Voronoi tessellation of any set of cosmological tracers (galaxies/haloes). We show that the shape of the VVF of biased tracers responds sensitively to physical properties such as halo mass, large-scale environment, substructure and redshift-space effects, making this a hitherto unexplored probe of both primordial cosmology and galaxy evolution. Using convenient summary statistics -- the width, median and a low percentile of the VVF as functions of average tracer number density -- we explore these effects for tracer populations in a suite of N-body simulations of a range of dark matter models. Our summary statistics sensitively probe primordial features such as small-scale oscillations in the initial matter power spectrum (as arise in models involving collisional effects in the dark sector), while being largely insensitive to a truncation of initial power (as in warm dark matter models). For vanilla cold dark matter (CDM) cosmologies, the summary statistics display strong evolution and redshift-space effects, and are also sensitive to cosmological parameter values for realistic tracer samples. Comparing the VVF of galaxies in the GAMA survey with that of abundance matched CDM (sub)haloes tentatively reveals environmental effects in GAMA beyond halo mass (modulo unmodelled satellite properties). Our exploratory analysis thus paves the way for using the VVF as a new probe of galaxy evolution physics as well as the nature of dark matter and dark energy. △ Less

Submitted 23 May, 2020; v1 submitted 23 January, 2020; originally announced January 2020.

Comments: 20 pages, 15 figures; v2 - minor changes to match version accepted in MNRAS

Cosmic web anisotropy is the primary indicator of halo assembly bias

Authors: Sujatha Ramakrishnan, Aseem Paranjape, Oliver Hahn, Ravi K. Sheth

Abstract: The internal properties of dark matter haloes correlate with the large-scale halo clustering strength at fixed halo mass $-$ an effect known as assembly bias $-$ and are also strongly affected by the local, non-linear cosmic web. Characterising a halo's local web environment by its tidal anisotropy $α$ at scales $\sim4$ x the halo radius, we demonstrate that these multi-scale correlations represen… ▽ More The internal properties of dark matter haloes correlate with the large-scale halo clustering strength at fixed halo mass $-$ an effect known as assembly bias $-$ and are also strongly affected by the local, non-linear cosmic web. Characterising a halo's local web environment by its tidal anisotropy $α$ at scales $\sim4$ x the halo radius, we demonstrate that these multi-scale correlations represent two distinct statistical links: one between the internal property and $α$, and the other between $α$ and large-scale ( $>30h^{-1}$Mpc) halo bias $b_1$. We focus on scalar internal properties of haloes related to formation time (concentration $c_{\rm vir}$), shape (mass ellipsoid asphericity $c/a$), velocity dispersion structure (velocity ellipsoid asphericity $c_v/a_v$ and velocity anisotropy $β$) and angular momentum (dimensionless spin $λ$) in the mass range $8\times10^{11}< M_{\rm vir}/(h^{-1}M_\odot)<5\times10^{14}$. Using conditional correlation coefficients and other detailed tests, we show that the joint distribution of $α$, $b_1$ and any of the internal properties $c\in\{β,c_v/a_v,c/a,c_{\rm vir},λ\}$ is consistent with $p(α,b_1,c)\simeq p(α)p(b_1|α)p(c|α)$, at all but the largest masses. $\textit{Thus, the assembly bias trends $c-b_1$ reflect the two fundamental correlations $c-α$ and $b_1-α$.}$ Our results are unaffected by the exclusion of haloes with recent major merger events or splashback objects, although the latter are distinguished by the fact that $α$ does not explain their assembly bias trends. The overarching importance of $α$ provides a new perspective on the nature of assembly bias of distinct haloes, with potential ramifications for incorporating realistic assembly bias effects into mock catalogs of future large-scale structure surveys and for detecting galaxy assembly bias. △ Less

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

Comments: 21 pages, 14 figures; v2: Added discussion, results strengthened. Conclusions unchanged. Accepted by MNRAS

Global analysis of luminosity- and colour-dependent galaxy clustering in the Sloan Digital Sky Survey

Authors: Niladri Paul, Isha Pahwa, Aseem Paranjape

Abstract: We present a Halo Occupation Distribution (HOD) analysis of the luminosity- and colour-dependent galaxy clustering in the Sloan Digital Sky Survey. A novelty of our technique is that it uses a combination of clustering measurements in luminosity bins to perform a global likelihood analysis, simultaneously constraining the HOD parameters for a range of luminosity thresholds. We present simple, smoo… ▽ More We present a Halo Occupation Distribution (HOD) analysis of the luminosity- and colour-dependent galaxy clustering in the Sloan Digital Sky Survey. A novelty of our technique is that it uses a combination of clustering measurements in luminosity bins to perform a global likelihood analysis, simultaneously constraining the HOD parameters for a range of luminosity thresholds. We present simple, smooth fitting functions which accurately describe the resulting luminosity dependence of the best-fit HOD parameters. To minimise systematic halo modelling effects, we use theoretical halo 2-point correlation functions directly measured and tabulated from a suite of $N$-body simulations spanning a large enough dynamic range in halo mass and spatial separation. Thus, our modelling correctly accounts for non-linear and scale-dependent halo bias as well as any departure of halo profiles from universality, and we additionally account for halo exclusion using the hard sphere approximation. Using colour-dependent clustering information, we constrain the satellite galaxy red fraction in a model-independent manner which does not rely on any group-finding algorithm. We find that the resulting luminosity dependence of the satellite red fraction is significantly shallower than corresponding measurements from galaxy group catalogues, and we provide a simple fitting function to describe this dependence. Our fitting functions are readily usable in generating low-redshift mock galaxy catalogues, and we discuss some potentially interesting applications as well as possible extensions of our technique. △ Less

Submitted 27 June, 2019; v1 submitted 5 February, 2019; originally announced February 2019.

Comments: 16 pages, 16 figures. (Accepted version in MNRAS with minor changes)

Photon number conservation and the large-scale 21 cm power spectrum in semi-numerical models of reionization

Authors: T. Roy Choudhury, Aseem Paranjape

Abstract: Semi-numerical models of the reionization of neutral hydrogen (HI) based on the excursion set (ES) approach are known to violate photon number conservation at the few per cent level. In this work, we highlight a more severe, previously unrecognized shortcoming of ES models: the large-scale 21 cm power spectrum (equivalently, HI bias b_HI) is a relatively strong function of the spatial resolution u… ▽ More Semi-numerical models of the reionization of neutral hydrogen (HI) based on the excursion set (ES) approach are known to violate photon number conservation at the few per cent level. In this work, we highlight a more severe, previously unrecognized shortcoming of ES models: the large-scale 21 cm power spectrum (equivalently, HI bias b_HI) is a relatively strong function of the spatial resolution used to generate ES ionization maps. We trace this problem to the fact that photon non-conservation in these models arises from a resolution-dependent mixture of spatially resolved, photon non-conserving bubbles, and partially ionized grid cells which are perfectly photon-conserving by construction. We argue that this inevitably leads to a resolution-dependence of b_HI, with the correct, converged value only emerging at very coarse resolution. Quantitatively, we find that b_HI can be non-converged by as much as ~20-25% in conservative ES implementations with grid sizes Delta x=5-10 cMpc/h, even when photon non-conservation is as small as ~3-4%. Thus, although numerically efficient, ES ionization maps coarse enough to produce a converged HI bias would wash out all topological features of the ionization field at scales k>~0.05 h/cMpc. We therefore present a new, explicitly photon conserving (PC) semi-numerical algorithm which distributes photons isotropically around sources while also accounting for anisotropic overlaps between nearby bubbles. Our PC algorithm predicts a resolution-independent value of b_HI consistent with the result of low-resolution ES maps, thus serving as a useful compromise between standard ES implementations and more expensive radiative transfer simulations. △ Less

Submitted 14 September, 2018; v1 submitted 2 July, 2018; originally announced July 2018.

Comments: Accepted for publication in MNRAS

The dependence of galaxy clustering on tidal environment in the Sloan Digital Sky Survey

Authors: Aseem Paranjape, Oliver Hahn, Ravi K. Sheth

Abstract: The influence of the Cosmic Web on galaxy formation and evolution is of great observational and theoretical interest. We investigate whether the Cosmic Web leaves an imprint in the spatial clustering of galaxies in the Sloan Digital Sky Survey (SDSS), using the group catalog of Yang et al. and tidal field estimates at $\sim2h^{-1}$Mpc scales from the Mass-Tides-Velocity data set of Wang et al. We… ▽ More The influence of the Cosmic Web on galaxy formation and evolution is of great observational and theoretical interest. We investigate whether the Cosmic Web leaves an imprint in the spatial clustering of galaxies in the Sloan Digital Sky Survey (SDSS), using the group catalog of Yang et al. and tidal field estimates at $\sim2h^{-1}$Mpc scales from the Mass-Tides-Velocity data set of Wang et al. We use the $\textit{tidal anisotropy}$ $α$ (Paranjape et al.) to characterise the tidal environment of groups, and measure the redshift-space 2-point correlation function (2pcf) of group positions and the luminosity- and colour-dependent clustering of group galaxies using samples segregated by $α$. We find that all the 2pcf measurements depend strongly on $α$, with factors of $\sim20$ between the large-scale 2pcf of objects in the most and least isotropic environments. To test whether these strong trends imply `beyond halo mass' effects for galaxy evolution, we compare our results with corresponding 2pcf measurements in mock catalogs constructed using a halo occupation distribution that only uses halo mass as an input. We find that this prescription qualitatively reproduces $\textit{all}$ observed trends, and also quantitatively matches many of the observed results. Although there are some statistically significant differences between our `halo mass only' mocks and the data -- in the most and least isotropic environments -- which deserve further investigation, our results suggest that if the tidal environment induces additional effects on galaxy properties other than those inherited from their host halos, then these must be weak. △ Less

Submitted 16 January, 2018; v1 submitted 14 January, 2018; originally announced January 2018.

Comments: 11 pages, 6 figures, submitted to MNRAS

Halo models of HI selected galaxies

Authors: Niladri Paul, Tirthankar Roy Choudhury, Aseem Paranjape

Abstract: Modelling the distribution of neutral hydrogen (HI) in dark matter halos is important for studying galaxy evolution in the cosmological context. We use a novel approach to infer the HI-dark matter connection at the massive end ($m_{\rm HI} > 10^{9.8} M_{\odot}$) from radio HI emission surveys, using optical properties of low-redshift galaxies as an intermediary. In particular, we use a previously… ▽ More Modelling the distribution of neutral hydrogen (HI) in dark matter halos is important for studying galaxy evolution in the cosmological context. We use a novel approach to infer the HI-dark matter connection at the massive end ($m_{\rm HI} > 10^{9.8} M_{\odot}$) from radio HI emission surveys, using optical properties of low-redshift galaxies as an intermediary. In particular, we use a previously calibrated optical HOD describing the luminosity- and colour-dependent clustering of SDSS galaxies and describe the HI content using a statistical scaling relation between the optical properties and HI mass. This allows us to compute the abundance and clustering properties of HI-selected galaxies and compare with data from the ALFALFA survey. We apply an MCMC-based statistical analysis to constrain the free parameters related to the scaling relation. The resulting best-fit scaling relation identifies massive HI galaxies primarily with optically faint blue centrals, consistent with expectations from galaxy formation models. We compare the HI-stellar mass relation predicted by our model with independent observations from matched HI-optical galaxy samples, finding reasonable agreement. As a further application, we make some preliminary forecasts for future observations of HI and optical galaxies in the expected overlap volume of SKA and Euclid/LSST. △ Less

Submitted 10 June, 2018; v1 submitted 12 December, 2017; originally announced December 2017.

Comments: 12 pages, 11 figures, v2: substantial revision, matches accepted version

Halo assembly bias and the tidal anisotropy of the local halo environment

Authors: Aseem Paranjape, Oliver Hahn, Ravi K. Sheth

Abstract: We study the role of the local tidal environment in determining the assembly bias of dark matter haloes. Previous results suggest that the anisotropy of a halo's environment (i.e, whether it lies in a filament or in a more isotropic region) can play a significant role in determining the eventual mass and age of the halo. We statistically isolate this effect using correlations between the large-sca… ▽ More We study the role of the local tidal environment in determining the assembly bias of dark matter haloes. Previous results suggest that the anisotropy of a halo's environment (i.e, whether it lies in a filament or in a more isotropic region) can play a significant role in determining the eventual mass and age of the halo. We statistically isolate this effect using correlations between the large-scale and small-scale environments of simulated haloes at $z=0$ with masses between $10^{11.6}\lesssim (m/h^{-1}M_{\odot})\lesssim10^{14.9}$. We probe the large-scale environment using a novel halo-by-halo estimator of linear bias. For the small-scale environment, we identify a variable $α_R$ that captures the $\textit{tidal anisotropy}$ in a region of radius $R=4R_{\textrm{200b}}$ around the halo and correlates strongly with halo bias at fixed mass. Segregating haloes by $α_R$ reveals two distinct populations. Haloes in highly isotropic local environments ($α_R\lesssim0.2$) behave as expected from the simplest, spherically averaged analytical models of structure formation, showing a $\textit{negative}$ correlation between their concentration and large-scale bias at $\textit{all}$ masses. In contrast, haloes in anisotropic, filament-like environments ($α_R\gtrsim0.5$) tend to show a $\textit{positive}$ correlation between bias and concentration at any mass. Our multi-scale analysis cleanly demonstrates how the overall assembly bias trend across halo mass emerges as an average over these different halo populations, and provides valuable insights towards building analytical models that correctly incorporate assembly bias. We also discuss potential implications for the nature and detectability of galaxy assembly bias. △ Less

Submitted 25 February, 2018; v1 submitted 29 June, 2017; originally announced June 2017.

Comments: 19 pages, 15 figures; v2: revised in response to referee comments, added references and discussion, conclusions unchanged. Accepted in MNRAS

Motifs in Temporal Networks

Authors: Ashwin Paranjape, Austin R. Benson, Jure Leskovec

Abstract: Networks are a fundamental tool for modeling complex systems in a variety of domains including social and communication networks as well as biology and neuroscience. Small subgraph patterns in networks, called network motifs, are crucial to understanding the structure and function of these systems. However, the role of network motifs in temporal networks, which contain many timestamped links betwe… ▽ More Networks are a fundamental tool for modeling complex systems in a variety of domains including social and communication networks as well as biology and neuroscience. Small subgraph patterns in networks, called network motifs, are crucial to understanding the structure and function of these systems. However, the role of network motifs in temporal networks, which contain many timestamped links between the nodes, is not yet well understood. Here we develop a notion of a temporal network motif as an elementary unit of temporal networks and provide a general methodology for counting such motifs. We define temporal network motifs as induced subgraphs on sequences of temporal edges, design fast algorithms for counting temporal motifs, and prove their runtime complexity. Our fast algorithms achieve up to 56.5x speedup compared to a baseline method. Furthermore, we use our algorithms to count temporal motifs in a variety of networks. Results show that networks from different domains have significantly different motif counts, whereas networks from the same domain tend to have similar motif counts. We also find that different motifs occur at different time scales, which provides further insights into structure and function of temporal networks. △ Less

Submitted 29 December, 2016; originally announced December 2016.

Journal ref: Proceedings of the Tenth ACM International Conference on Web Search and Data Mining, 2017