title: A generalized Kasner transition for bouncing Bianchi I models in modified gravity theories

abstract: We derive transition rules for Kasner exponents in bouncing Bianchi I models with generic perfect fluid matter fields for a broad class of modified gravity theories where cosmological singularities are resolved and replaced by a non-singular bounce. This is a generalization of results obtained previously in limiting curvature mimetic gravity and loop quantum cosmology. A geometric interpretation is provided for the transition rule as a linear map in the Kasner plane. We show that the general evolution of anisotropies in a Bianchi I universe -- including during the bounce phase -- is equivalent to the motion of a point particle on a sphere, where the sphere is the one-point compactification of the Kasner plane. In addition, we study the evolution of anisotropies in a large family of bouncing Bianchi I space-times. We also present a novel explicit solution to the Einstein equations for a Bianchi I universe with ekpyrotic matter with a constant equation of state $\omega=3$.


title: A Deep Learning Framework for Detection of Targets in Thermal Images to Improve Firefighting

abstract: Intelligent detection and processing capabilities can be instrumental to improving the safety, efficiency, and successful completion of rescue missions conducted by firefighters in emergency first response settings. The objective of this research is to create an automated system that is capable of real-time, intelligent object detection and recognition and facilitates the improved situational awareness of firefighters during an emergency response. We have explored state of the art machine/deep learning techniques to achieve this objective. The goal for this work is to enhance the situational awareness of firefighters by effectively exploiting the information gathered from infrared cameras carried by firefighters. To accomplish this, we use a trained deep Convolutional Neural Network (CNN) system to classify and identify objects of interest from thermal imagery in real time. In the midst of those critical circumstances created by structure fire, this system is able to accurately inform the decision making process of firefighters with real-time up-to-date scene information by extracting, processing, and analyzing crucial information. With the new information produced by the framework, firefighters are able to make more informed inferences about the circumstances for their safe navigation through such hazardous and potentially catastrophic environments.


title: Leavers and remainers: Galaxies split by group-exit

abstract: The disruption of substructure in galaxy clusters likely plays an important role in shaping the cluster population as a significant fraction of cluster galaxies today have spent time in a previous host system, and thus may have been pre-processed. Once inside the cluster, group galaxies face the combined environmental effects from group and cluster - so called 'post-processing'. We investigate these concepts, by tracking the evolution of satellites and their hosts after entering the cluster and find that tidal forces during their first pericentric passage are very efficient at breaking up groups, preferentially removing satellites at larger distances from their hosts. 92.2% of satellites whose host has passed pericentre will leave their host by $z=0$, typically no later than half a Gyr after pericentric passage. We find satellites leave with high velocities, and quickly separate to large distances from their hosts, making their identification within the cluster population challenging. Those few satellites ($\sim$7.8%) that remain bound to their hosts after a pericentric passage are typically found close to their host centres. This implies that substructure seen in clusters today is very likely on first infall into the cluster, and yet to pass pericentre. This is even more likely if the substructure is extended, with satellites beyond R$_{200}$ of their host. We find the cluster dominates the tidal mass loss and destruction of satellites, and is responsible for rapidly halting the accretion of new satellites onto hosts once they reach 0.6-0.8 R$_{200}$ radii from the cluster.


title: Resilience of the Rank of Random Matrices

abstract: Let $M$ be an $n \times m$ matrix of independent Rademacher ($\pm 1$) random variables. It is well known that if $n \leq m$, then $M$ is of full rank with high probability. We show that this property is resilient to adversarial changes to $M$. More precisely, if $m \geq n + n^{1-\varepsilon/6}$, then even after changing the sign of $(1-\varepsilon)m/2$ entries, $M$ is still of full rank with high probability. Note that this is asymptotically best possible as one can easily make any two rows proportional with at most $m/2$ changes. Moreover, this theorem gives an asymptotic solution to a slightly weakened version of a conjecture made by Van Vu.


title: Receding Horizon Curiosity

abstract: Sample-efficient exploration is crucial not only for discovering rewarding experiences but also for adapting to environment changes in a task-agnostic fashion. A principled treatment of the problem of optimal input synthesis for system identification is provided within the framework of sequential Bayesian experimental design. In this paper, we present an effective trajectory-optimization-based approximate solution of this otherwise intractable problem that models optimal exploration in an unknown Markov decision process (MDP). By interleaving episodic exploration with Bayesian nonlinear system identification, our algorithm takes advantage of the inductive bias to explore in a directed manner, without assuming prior knowledge of the MDP. Empirical evaluations indicate a clear advantage of the proposed algorithm in terms of the rate of convergence and the final model fidelity when compared to intrinsic-motivation-based algorithms employing exploration bonuses such as prediction error and information gain. Moreover, our method maintains a computational advantage over a recent model-based active exploration (MAX) algorithm, by focusing on the information gain along trajectories instead of seeking a global exploration policy. A reference implementation of our algorithm and the conducted experiments is publicly available.


title: On smooth mesoscopic linear statistics of the eigenvalues of random permutation matrices

abstract: We study the limiting behavior of smooth linear statistics of the spectrum of random permutation matrices in the mesoscopic regime, when the permutation follows one of the Ewens measures on the symmetric group. If we apply a smooth enough test function $f$ to all the determinations of the eigenangles of the permutations, we get a convergence in distribution when the order of the permutation tends to infinity. Two distinct kinds of limit appear: if $f(0)\neq 0$, we have a central limit theorem with a logarithmic variance, and if $f(0) = 0$, the convergence holds without normalization and the limit involves a scale-invariant Poisson point process.


title: Mobile App Privacy in Software Engineering Research: A Systematic Mapping Study

abstract: Mobile applications (apps) have become deeply personal, constantly demanding access to privacy-sensitive information in exchange for more personalized user experiences. Such privacy-invading practices have generated major multidimensional and unconventional privacy concerns among app users. To address these concerns, the research on mobile app privacy has experienced rapid growth over the past decade. In general, this line of research is aimed at systematically exposing the privacy practices of apps and proposing solutions to protect the privacy of mobile app users. In this survey paper, we conduct a systematic mapping study of 54 Software Engineering (SE) primary studies on mobile app privacy. Our objectives are to a) explore trends in SE app privacy research, b) categorize existing evidence, and c) identify potential directions for future research. Our results show that existing literature can be divided into four main categories: privacy policy, requirements, user perspective, and leak detection. Furthermore, our survey reveals an imbalance between these categories; majority of existing research focuses on proposing tools for detecting privacy leaks, with less studies targeting privacy requirements and policy and even less on user perspective. Finally, our survey exposes several gaps in existing research and suggests areas for improvement.


title: Invariable generation of finite classical groups

abstract: A subset of a group invariably generates the group if it generates even when we replace the elements by any of their conjugates. In a 2016 paper, Pemantle, Peres and Rivin show that the probability that four randomly selected elements invariably generate $S_n$ is bounded away from zero by an absolute constant for all $n$. Subsequently, Eberhard, Ford and Green have shown that the probability that three randomly selected elements invariably generate $S_n$ tends to zero as $n \rightarrow \infty$. In this paper, we prove an analogous result for the finite classical groups. More precisely, let $G_r(q)$ be a finite classical group of rank $r$ over $\mathbb{F}_q$. We show that for $q$ large enough, the probability that four randomly selected elements invariably generate $G_r(q)$ is bounded away from zero by an absolute constant for all $r$, and for three elements the probability tends to zero as $q \rightarrow \infty$ and $r \rightarrow \infty$. We use the fact that most elements in $G_r(q)$ are separable and the well-known correspondence between classes of maximal tori containing separable elements in classical groups and conjugacy classes in their Weyl groups.


title: SmoothFool: An Efficient Framework for Computing Smooth Adversarial Perturbations

abstract: Deep neural networks are susceptible to adversarial manipulations in the input domain. The extent of vulnerability has been explored intensively in cases of $\ell_p$-bounded and $\ell_p$-minimal adversarial perturbations. However, the vulnerability of DNNs to adversarial perturbations with specific statistical properties or frequency-domain characteristics has not been sufficiently explored. In this paper, we study the smoothness of perturbations and propose SmoothFool, a general and computationally efficient framework for computing smooth adversarial perturbations. Through extensive experiments, we validate the efficacy of the proposed method for both the white-box and black-box attack scenarios. In particular, we demonstrate that: (i) there exist extremely smooth adversarial perturbations for well-established and widely used network architectures, (ii) smoothness significantly enhances the robustness of perturbations against state-of-the-art defense mechanisms, (iii) smoothness improves the transferability of adversarial perturbations across both data points and network architectures, and (iv) class categories exhibit a variable range of susceptibility to smooth perturbations. Our results suggest that smooth APs can play a significant role in exploring the vulnerability extent of DNNs to adversarial examples.


title: On relations between transportation cost spaces and $\ell_1$

abstract: The present paper deals with some structural properties of transportation cost spaces, also known as Arens-Eells spaces, Lipschitz-free spaces and Wasserstein spaces. The main results of this work are: (1) A necessary and sufficient condition on an infinite metric space $M$, under which the transportation cost space on $M$ contains an isometric copy of $\ell_1$. The obtained condition is applied to answer the open questions asked by C\'uth and Johanis (2017) concerning several specific metric spaces. (2) The description of the transportation cost space of a weighted finite graph $G$ as the quotient $\ell_1(E(G))/Z(G)$, where $E(G)$ is the edge set and $Z(G)$ is the cycle space of $G$. This is a generalization of the previously known result to the case of any finite metric space.


title: In the Presence of a Wrecking Ball: Orbital Stability in the HR 5183 System

abstract: Discoveries of exoplanets using the radial velocity method are progressively reaching out to increasingly longer orbital periods as the duration of surveys continues to climb. The improving sensitivity to potential Jupiter analogs is revealing a diversity of orbital architectures that are substantially different from that found in our solar system. An excellent example of this is the recent discovery of HR 5183b; a giant planet on a highly eccentric ($e = 0.84$) $\sim$75~year orbit. The presence of such giant planet orbits are intrinsically interesting from the perspective of the dynamical history of planetary systems, and also for examining the implications of on-going dynamical stability and habitability of these systems. In this work, we examine the latter, providing results of dynamical simulations that explore the stable regions that the eccentric orbit of the HR 5183 giant planet allows to exist within the Habitable Zone of the host star. Our results show that, despite the incredible perturbing influence of the giant planet, there remain a narrow range of locations within the Habitable Zone where terrestrial planets may reside in long-term stable orbits. We discuss the effects of the giant planet on the potential habitability of a stable terrestrial planet, including the modulation of terrestrial planet eccentricities and the periodically spectacular view of the giant planet from the terrestrial planet location.


title: Controlling Costs: Feature Selection on a Budget

abstract: The traditional framework for feature selection treats all features as costing the same amount. However, in reality, a scientist often has considerable discretion regarding which variables to measure, and the decision involves a tradeoff between model accuracy and cost (where cost can refer to money, time, difficulty, or intrusiveness). In particular, unnecessarily including an expensive feature in a model is worse than unnecessarily including a cheap feature. We propose a procedure, which we call cheap knockoffs, for performing feature selection in a cost-conscious manner. The key idea behind our method is to force higher cost features to compete with more knockoffs than cheaper features. We derive an upper bound on the weighted false discovery proportion associated with this procedure, which corresponds to the fraction of the feature cost that is wasted on unimportant features. We prove that this bound holds simultaneously with high probability over a path of selected variable sets of increasing size. A user may thus select a set of features based, for example, on the overall budget, while knowing that no more than a particular fraction of feature cost is wasted. We investigate, through simulation and a biomedical application, the practical importance of incorporating cost considerations into the feature selection process.


title: The role of mainstreamness and interdisciplinarity for the relevance of scientific papers

abstract: There is demand from science funders, industry, and the public that science should become more risk-taking, more out-of-the-box, and more interdisciplinary. Is it possible to tell how interdisciplinary and out-of-the-box scientific papers are, or which papers are mainstream? Here we use the bibliographic coupling network, derived from all physics papers that were published in the Physical Review journals in the past century, to try to identify them as mainstream, out-of-the-box, or interdisciplinary. We show that the network clusters into scientific fields. The position of individual papers with respect to these clusters allows us to estimate their degree of mainstreamness or interdisciplinary. We show that over the past decades the fraction of mainstream papers increases, the fraction of out-of-the-box decreases, and the fraction of interdisciplinary papers remains constant. Studying the rewards of papers, we find that in terms of absolute citations, both, mainstream and interdisciplinary papers are rewarded. In the long run, mainstream papers perform less than interdisciplinary ones in terms of citation rates. We conclude that to avoid a trend towards mainstreamness a new incentive scheme is necessary.


title: Symmetry enhancement in 4d Spin(n) gauge theories and compactification from 6d

abstract: We consider a known sequence of dualities involving $4d$ ${\cal N}=1$ theories with $Spin(n)$ gauge groups and use it to construct a new sequence of models exhibiting IR symmetry enhancement. Then, motivated by the observed pattern of IR symmetries we conjecture six-dimensional theories the compactification of which on a Riemann surface yields the $4d$ sequence of models along with their symmetry enhancements, and put them to several consistency checks.


title: An introduction to a general records theory both for dependent data and high dimensions: revisited 2

abstract: The probabilistic investigation on record values and record times of a sequence of random variables defined on the same probability space has received much attention from 1952 to now. A great deal of such theory focused on \textit{iid} or independent real-valued random variables. There exists a few results for real-valued dependent random variables. Some papers deal also with multivariate random variables. But a large theory regarding vectors and dependent data has yet to be done. In preparation of that, the probability laws of records are investigated here, without any assumption on the dependence structure. The results are extended sequences with values in partially ordered spaces whose order is compatible with measurability. The general characterizations are checked in known cases mostly for \textit{iid} sequences. The frame is ready for undertaking a vast study of records theory in high dimensions and for types of dependence.


title: ConvoSource: Radio-Astronomical Source-Finding with Convolutional Neural Networks

abstract: Finding and classifying astronomical sources is key in the scientific exploitation of radio surveys. Source-finding usually involves identifying the parts of an image belonging to an astronomical source, against some estimated background. This can be problematic in the radio regime, owing to the presence of correlated noise, which can interfere with the signal from the source. In the current work, we present ConvoSource, a novel method based on a deep learning technique, to identify the positions of radio sources, and compare the results to a Gaussian-fitting method. Since the deep learning approach allows the generation of more training images, it should perform well in the source-finding task. We test the source-finding methods on artificial data created for the data challenge of the Square Kilometer Array (SKA). We investigate sources that are divided into three classes: star forming galaxies (SFGs) and two classes of active galactic nuclei (AGN). The~artificial data are given at two different frequencies (560~MHz and 1400~MHz), three total integration times (8 h, 100 h, 1000 h), and three signal-to-noise ratios (SNRs) of 1, 2, and 5. At~lower SNRs, ConvoSource tends to outperform a Gaussian-fitting approach in the recovery of SFGs and all sources, although at the lowest SNR of one, the better performance is likely due to chance matches. The~Gaussian-fitting method performs better in the recovery of the AGN-type sources at lower SNRs. At~a higher SNR, ConvoSource performs better on average in the recovery of AGN sources, whereas the Gaussian-fitting method performs better in the recovery of SFGs and all sources. ConvoSource usually performs better at shorter total integration times and detects more true positives and misses fewer sources compared to the Gaussian-fitting method; however, it detects more false positives.


title: Distilling Importance Sampling for Likelihood Free Inference

abstract: Likelihood-free inference involves inferring parameter values given observed data and a simulator model. The simulator is computer code which takes parameters, performs stochastic calculations, and outputs simulated data. In this work, we view the simulator as a function whose inputs are (1) the parameters and (2) a vector of pseudo-random draws. We attempt to infer all these inputs conditional on the observations. This is challenging as the resulting posterior can be high dimensional and involve strong dependence. We approximate the posterior using normalizing flows, a flexible parametric family of densities. Training data is generated by likelihood-free importance sampling with a large bandwidth value epsilon, which makes the target similar to the prior. The training data is "distilled" by using it to train an updated normalizing flow. The process is iterated, using the updated flow as the importance sampling proposal, and slowly reducing epsilon so the target becomes closer to the posterior. Unlike most other likelihood-free methods, we avoid the need to reduce data to low dimensional summary statistics, and hence can achieve more accurate results. We illustrate our method in two challenging examples, on queuing and epidemiology.


title: Multi-Vehicle Interaction Scenarios Generation with Interpretable Traffic Primitives and Gaussian Process Regression

abstract: Generating multi-vehicle interaction scenarios can benefit motion planning and decision making of autonomous vehicles when on-road data is insufficient. This paper presents an efficient approach to generate varied multi-vehicle interaction scenarios that can both adapt to different road geometries and inherit the key interaction patterns in real-world driving. Towards this end, the available multi-vehicle interaction scenarios are temporally segmented into several interpretable fundamental building blocks, called traffic primitives, via the Bayesian nonparametric learning. Then, the changepoints of traffic primitives are transformed into the desired road to generate collision-free interaction trajectories through a sampling-based path planning algorithm. The Gaussian process regression is finally introduced to control the variance and smoothness of the generated multi-vehicle interaction trajectories. Experiments with simulation results of three typical multi-vehicle trajectories at different road conditions are carried out. The experimental results demonstrate that our proposed method can generate a bunch of human-like multi-vehicle interaction trajectories that can fit different road conditions remaining the key interaction patterns of agents in the provided scenarios, which is import to the development of autonomous vehicles.


title: Prose for a Painting

abstract: Painting captions are often dry and simplistic which motivates us to describe a painting creatively in the style of Shakespearean prose. This is a difficult problem, since there does not exist a large supervised dataset from paintings to Shakespearean prose. Our solution is to use an intermediate English poem description of the painting and then apply language style transfer which results in Shakespearean prose describing the painting. We rate our results by human evaluation on a Likert scale, and evaluate the quality of language style transfer using BLEU score as a function of prose length. We demonstrate the applicability and limitations of our approach by generating Shakespearean prose for famous paintings. We make our models and code publicly available.


title: On vertex-edge and independent vertex-edge domination

abstract: Given a graph $G = (V,E)$, a vertex $u \in V$ ve-dominates all edges incident to any vertex of $N_G[u]$. A set $S \subseteq V$ is a ve-dominating set if for all edges $e\in E$, there exists a vertex $u \in S$ such that $u$ ve-dominates $e$. Lewis [Ph.D. thesis, 2007] proposed a linear time algorithm for ve-domination problem for trees. In this paper, first we have constructed an example where the proposed algorithm fails. Then we have proposed a linear time algorithm for ve-domination problem in block graphs, which is a superclass of trees. We have also proved that finding minimum ve-dominating set is NP-complete for undirected path graphs. Finally, we have characterized the trees with equal ve-domination and independent ve-domination number.


title: Reconstructing Functions and Estimating Parameters with Artificial Neural Networks: A Test with the Hubble Parameter and SNe Ia

abstract: In this work, we propose a new nonparametric approach for reconstructing a function from observational data using an Artificial Neural Network (ANN), which has no assumptions about the data and is a completely data-driven approach. We test the ANN method by reconstructing functions of the Hubble parameter measurements $H(z)$ and the distance-redshift relation $D_L(z)$ of Type Ia supernovae. We find that both $H(z)$ and $D_L(z)$ can be reconstructed with high accuracy. Furthermore, we estimate cosmological parameters using the reconstructed functions of $H(z)$ and $D_L(z)$ and find the results are consistent with those obtained using the observational data directly. Therefore, we propose that the function reconstructed by ANN can represent the actual distribution of observational data and can be used for parameter estimation in further cosmological research. In addition, we present a new strategy for training and evaluating the neural network, and a code for reconstructing functions using ANN has been developed and is available for download.


title: Decoding Quantum Field Theory with Machine Learning

abstract: We demonstrate how one can use machine learning techniques to bypass the technical difficulties of designing an experiment and translating its outcomes into concrete claims about fundamental features of quantum fields. In practice, all measurements of quantum fields are carried out through local probes. Despite measuring only a small portion of the field, such local measurements have the capacity to reveal many of the field's global features. This is because, when in equilibrium with their environments, quantum fields store global information locally, albeit in a scrambled way. We show that neural networks can be trained to unscramble this information from data generated from a very simple one-size-fits-all local measurement protocol. To illustrate this general claim we will consider three non-trivial features of the field as case studies: a) how, as long as the field is in a stationary state, a particle detector can learn about the field's boundary conditions even before signals have time to propagate from the boundary to the detector, b) how detectors can determine the temperature of the quantum field even without thermalizing with it, and c) how detectors can distinguish between Fock states and coherent states even when the first and second moments of all their quadrature operators match. Each of these examples uses the exact same simple fixed local measurement protocol and machine-learning ansatz successfully. This supports the claim that the framework proposed here can be applied to nearly any kind of local measurement on a quantum field to reveal nearly any of the field's global properties in a one-size-fits-all manner.


title: Bregman Proximal Framework for Deep Linear Neural Networks

abstract: A typical assumption for the analysis of first order optimization methods is the Lipschitz continuity of the gradient of the objective function. However, for many practical applications this assumption is violated, including loss functions in deep learning. To overcome this issue, certain extensions based on generalized proximity measures known as Bregman distances were introduced. This initiated the development of the Bregman proximal gradient (BPG) algorithm and an inertial variant (momentum based) CoCaIn BPG, which however rely on problem dependent Bregman distances. In this paper, we develop Bregman distances for using BPG methods to train Deep Linear Neural Networks. The main implications of our results are strong convergence guarantees for these algorithms. We also propose several strategies for their efficient implementation, for example, closed form updates and a closed form expression for the inertial parameter of CoCaIn BPG. Moreover, the BPG method requires neither diminishing step sizes nor line search, unlike its corresponding Euclidean version. We numerically illustrate the competitiveness of the proposed methods compared to existing state of the art schemes.


title: Dark Energy from Topology

abstract: The concordance model of cosmology suffers from the major theoretical problems surrounding the observed value and recent emergence of a cosmological constant. In this paper we present a novel approach, which explains more naturally its value than that based on quantum vacuum energy, in the form of topological invariants characteristic classes, included as Lagrange multipliers in the action. The approach draws from topological as well as dynamical system consideration, generating as a byproduct an effective cosmological constant. General Relativity is recovered by canceling the torsion in a region containing the observable Universe, which boundary constraints the invariants, thus yielding the effective cosmological constant's form. As that form's denominator contains the total volume of the average black hole, calculated from a geometrical mean on the estimated black hole mass distribution and directly associated to the ratio of the total volume boundary of the space-time manifold and the dominant term in its Euler characteristic. The constant's small estimated value compared to the Planck scale is therefore natural and our evaluation fits remarkably well with the observed value.


title: MAMS-A*: Multi-Agent Multi-Scale A*

abstract: We present a multi-scale forward search algorithm for distributed agents to solve single-query shortest path planning problems. Each agent first builds a representation of its own search space of the common environment as a multi-resolution graph, it communicates with the other agents the result of its local search, and it uses received information from other agents to refine its own graph and update the local inconsistency conditions. As a result, all agents attain a common subgraph that includes a provably optimal path in the most informative graph available among all agents, if one exists, without necessarily communicating the entire graph. We prove the completeness and optimality of the proposed algorithm, and present numerical results supporting the advantages of the proposed approach.


title: Linking emotions to behaviors through deep transfer learning

abstract: Human behavior refers to the way humans act and interact. Understanding human behavior is a cornerstone of observational practice, especially in psychotherapy. An important cue of behavior analysis is the dynamical changes of emotions during the conversation. Domain experts integrate emotional information in a highly nonlinear manner, thus, it is challenging to explicitly quantify the relationship between emotions and behaviors. In this work, we employ deep transfer learning to analyze their inferential capacity and contextual importance. We first train a network to quantify emotions from acoustic signals and then use information from the emotion recognition network as features for behavior recognition. We treat this emotion-related information as behavioral primitives and further train higher level layers towards behavior quantification. Through our analysis, we find that emotion-related information is an important cue for behavior recognition. Further, we investigate the importance of emotional-context in the expression of behavior by constraining (or not) the neural networks' contextual view of the data. This demonstrates that the sequence of emotions is critical in behavior expression. To achieve these frameworks we employ hybrid architectures of convolutional networks and recurrent networks to extract emotion-related behavior primitives and facilitate automatic behavior recognition from speech.


title: From Convex Ideal Polyhedra to Fundamental Domains in H^3

abstract: Our goal is to better understand the relationship between the polyhedron and the group associated with a fundamental domain in H^3. In this paper, we will study torsion-free groups and determine a formula for how many edge classes a given abstract polyhedron must have. We will use that result to classify all fundamental domains on the cube with torsion-free groups, including a discussion of the explicit groups associated to those domains. We will then turn to more general fundamental domains and prove a series of results about how properties of the group place restrictions on the edge classes in the quotient manifold. These results give insight into how the polyhedron and the group associated to a fundamental domain interact, as well as offering concrete tools to find fundamental domains.


title: Variance reduction for Markov chains with application to MCMC

abstract: In this paper we propose a novel variance reduction approach for additive functionals of Markov chains based on minimization of an estimate for the asymptotic variance of these functionals over suitable classes of control variates. A distinctive feature of the proposed approach is its ability to significantly reduce the overall finite sample variance. This feature is theoretically demonstrated by means of a deep non asymptotic analysis of a variance reduced functional as well as by a thorough simulation study. In particular we apply our method to various MCMC Bayesian estimation problems where it favourably compares to the existing variance reduction approaches.

categories: math.ST cs.LG math.PR stat.CO stat.ML stat.TH


title: Stochastic triangular mesh mapping: A terrain mapping technique for autonomous mobile robots

abstract: For mobile robots to operate autonomously in general environments, perception is required in the form of a dense metric map. For this purpose, we present the stochastic triangular mesh (STM) mapping technique: a 2.5-D representation of the surface of the environment using a continuous mesh of triangular surface elements, where each surface element models the mean plane and roughness of the underlying surface. In contrast to existing mapping techniques, a STM map models the structure of the environment by ensuring a continuous model, while also being able to be incrementally updated with linear computational cost in the number of measurements. We reduce the effect of uncertainty in the robot pose (position and orientation) by using landmark-relative submaps. The uncertainty in the measurements and robot pose are accounted for by the use of Bayesian inference techniques during the map update. We demonstrate that a STM map can be used with sensors that generate point measurements, such as light detection and ranging (LiDAR) sensors and stereo cameras. We show that a STM map is a more accurate model than the only comparable online surface mapping technique$\unicode{x2014}$a standard elevation map$\unicode{x2014}$and we also provide qualitative results on practical datasets.


title: Span-core Decomposition for Temporal Networks: Algorithms and Applications

abstract: When analyzing temporal networks, a fundamental task is the identification of dense structures (i.e., groups of vertices that exhibit a large number of links), together with their temporal span (i.e., the period of time for which the high density holds). In this paper we tackle this task by introducing a notion of temporal core decomposition where each core is associated with two quantities, its coreness, which quantifies how densely it is connected, and its span, which is a temporal interval: we call such cores \emph{span-cores}. For a temporal network defined on a discrete temporal domain $T$, the total number of time intervals included in $T$ is quadratic in $|T|$, so that the total number of span-cores is potentially quadratic in $|T|$ as well. Our first main contribution is an algorithm that, by exploiting containment properties among span-cores, computes all the span-cores efficiently. Then, we focus on the problem of finding only the \emph{maximal span-cores}, i.e., span-cores that are not dominated by any other span-core by both their coreness property and their span. We devise a very efficient algorithm that exploits theoretical findings on the maximality condition to directly extract the maximal ones without computing all span-cores. Finally, as a third contribution, we introduce the problem of \emph{temporal community search}, where a set of query vertices is given as input, and the goal is to find a set of densely-connected subgraphs containing the query vertices and covering the whole underlying temporal domain $T$. We derive a connection between this problem and the problem of finding (maximal) span-cores. Based on this connection, we show how temporal community search can be solved in polynomial-time via dynamic programming, and how the maximal span-cores can be profitably exploited to significantly speed-up the basic algorithm.


title: Topological charge pumping of bound bosonic pairs

abstract: Experiments with bosonic atoms in optical superlattices allow for the interesting possibility to study the adiabatic quantized pumping of bosonic atoms in the presence of interactions. We investigate this exotic phenomenon for bound bosonic pairs in the paradigmatic Su-Schrieffer-Heeger model where the ground state exhibits topological phase transitions due to dimerized hoppings. At unit filling we show that there exist crossovers and phase transitions to bond-order phases of paired bosons known as pair-bond-order phase as a function of attractive interactions. The pair bond order phase is found to exhibit effective topological properties such as the presence of polarized paired edge states. This is further analyzed by studying the emergence and breakdown of the Thouless charge pumping of this bound bosonic pairs by a parametric extension to the famous Rice-Mele model. Finally we discuss how the pumping of paired bosons or different regimes of breakdown of charge pumping can be probed by state-of-the art experiments with repulsively bound bosons.


title: On the fibration method for rational points

abstract: We study weak approximation on rationally connected varieties under an assumption of strong approximation for a "simple" variety or under Schinzel's hypothesis. We also get some unconditional results.


title: Meta-Transfer Learning through Hard Tasks

abstract: Meta-learning has been proposed as a framework to address the challenging few-shot learning setting. The key idea is to leverage a large number of similar few-shot tasks in order to learn how to adapt a base-learner to a new task for which only a few labeled samples are available. As deep neural networks (DNNs) tend to overfit using a few samples only, typical meta-learning models use shallow neural networks, thus limiting its effectiveness. In order to achieve top performance, some recent works tried to use the DNNs pre-trained on large-scale datasets but mostly in straight-forward manners, e.g., (1) taking their weights as a warm start of meta-training, and (2) freezing their convolutional layers as the feature extractor of base-learners. In this paper, we propose a novel approach called meta-transfer learning (MTL) which learns to transfer the weights of a deep NN for few-shot learning tasks. Specifically, meta refers to training multiple tasks, and transfer is achieved by learning scaling and shifting functions of DNN weights for each task. In addition, we introduce the hard task (HT) meta-batch scheme as an effective learning curriculum that further boosts the learning efficiency of MTL. We conduct few-shot learning experiments and report top performance for five-class few-shot recognition tasks on three challenging benchmarks: miniImageNet, tieredImageNet and Fewshot-CIFAR100 (FC100). Extensive comparisons to related works validate that our MTL approach trained with the proposed HT meta-batch scheme achieves top performance. An ablation study also shows that both components contribute to fast convergence and high accuracy.


title: Darondeau-Pragacz formulas in complex cobordism

abstract: In this paper, we generalize the push-forward (Gysin) formulas for flag bundles in the ordinary cohomology theory, which are due to Darondeau-Pragacz, to the complex cobordism theory. Then we introduce the {\it universal quadratic Schur functions}, which are a generalization of the (ordinary) quadratic Schur functions introduced by Darondeau-Pragacz, and establish some Gysin formulas for the universal quadratic Schur functions as an application of our Gysin formulas.


title: Multi-Head Attention for Multi-Modal Joint Vehicle Motion Forecasting

abstract: This paper presents a novel vehicle motion forecasting method based on multi-head attention. It produces joint forecasts for all vehicles on a road scene as sequences of multi-modal probability density functions of their positions. Its architecture uses multi-head attention to account for complete interactions between all vehicles, and long short-term memory layers for encoding and forecasting. It relies solely on vehicle position tracks, does not need maneuver definitions, and does not represent the scene with a spatial grid. This allows it to be more versatile than similar model while combining any forecasting capabilities, namely joint forecast with interactions, uncertainty estimation, and multi-modality. The resulting prediction likelihood outperforms state-of-the-art models on the same dataset.


title: Improved Bounds for Two Query Adaptive Bitprobe Schemes Storing Five Elements

abstract: In this paper, we study two-bitprobe adaptive schemes storing five elements. For these class of schemes, the best known lower bound is m^{1/2} due to Alon and Feige [SODA 2009]. Recently, it was proved by Kesh [FSTTCS 2018] that two-bitprobe adaptive schemes storing three elements will take at least m^{2/3} space, which also puts a lower bound on schemes storing five elements. In this work, we have improved the lower bound to m^{3/4}. We also present a scheme for the same that takes O(m^{5/6}) space. This improves upon the O(m^{18/19})-scheme due to Garg [Ph.D. Thesis] and the O(m^{10/11})-scheme due to Baig et al. [WALCOM 2019].


title: Canadian Investigations of the Interstellar Medium

abstract: The interstellar medium mediates galactic evolution as the reservoir of material for future star formation and the repository of energy and matter output by stellar processes. Canadians have played leading roles in ISM science for decades. The Canadian Galactic Plane Survey identified a wealth of small-scale structure in H I emission as well as self-absorption and in the structure of polarized emission. These observations demonstrated that no phase of the ISM, including the transition from atomic gas to star formation, can be understood in isolation. Canadians have also played leading roles in the characterization of dust with Planck and balloon-borne telescopes. Canadians have also used pulsar scintillometry to measure structure in the ISM at the smallest scales, below 1 AU. The 2020s offer many opportunities for ISM science in Canada. A major but cost-effective upgrade to the Synthesis Telescope with broadband (400-1800 MHz) single-pixel feeds would enable broadband polarimetry as well as wide-area, arcminute surveys of radio recombination lines. The next generation of balloon-borne telescopes will investigate magnetic fields and dust properties. Large single dishes, particularly the Green Bank Telescope, remain essential for our understanding of the diffuse structure which characterizes the ISM. Very long baseline interferometry capability enables parallax measurements of pulsars and masers and for further progress in scintillometry. ISM astronomers will continue to participate in cosmological experiments including CHIME and CHORD. Protecting quiet radio frequency interference environments will be ever more important as broadband observations are ever more central to ISM science. Computational capability is essential both for numerical work and for handling the observational data.


title: Schauder estimates for degenerate stable Kolmogorov equations

abstract: We provide here global Schauder-type estimates for a chain of integro-partial differential equations (IPDE) driven by a degenerate stable Ornstein-Uhlenbeck operator possibly perturbed by a deterministic drift, when the coefficients lie in some suitable anisotropic H{\"o}lder spaces. Our approach mainly relies on a perturbative method based on forward parametrix expansions and, due to the low regularizing properties on the degenerate variables and to some integrability constraints linked to the stability index, it also exploits duality results between appropriate Besov Spaces. In particular, our method also applies in some super-critical cases. Thanks to these estimates, we show in addition the well-posedness of the considered IPDE in a suitable functional space.


title: Analytical Models of Phonon-Point Defect Scattering

abstract: Point defects exist widely in engineering materials and are known to scatter vibrational modes to reduce thermal conductivity. The Klemens description of point defect scattering is the most prolific analytical model for this effect. This work reviews the essential physics of the model and compares its predictions to first principles results for isotope and alloy scattering, demonstrating the model as a useful materials design metric. A treatment of the scattering parameter ($\Gamma$) for a multiatomic lattice is recommended and compared to other treatments presented in literature, which have been at times misused to yield incomplete conclusions about the system's scattering mechanisms. Additionally, we demonstrate a reduced sensitivity of the model to the full phonon dispersion and discuss its origin. Finally, a simplified treatment of scattering in alloy systems with vacancies and interstitial defects is demonstrated to suitably describe the potent scattering strength of these off-stoichiometric defects.


title: Executing Instructions in Situated Collaborative Interactions

abstract: We study a collaborative scenario where a user not only instructs a system to complete tasks, but also acts alongside it. This allows the user to adapt to the system abilities by changing their language or deciding to simply accomplish some tasks themselves, and requires the system to effectively recover from errors as the user strategically assigns it new goals. We build a game environment to study this scenario, and learn to map user instructions to system actions. We introduce a learning approach focused on recovery from cascading errors between instructions, and modeling methods to explicitly reason about instructions with multiple goals. We evaluate with a new evaluation protocol using recorded interactions and online games with human users, and observe how users adapt to the system abilities.


title: Bayesian open games

abstract: This paper generalises the treatment of compositional game theory as introduced by Ghani et al. in 2018, where games are modelled as morphisms of a symmetric monoidal category. From an economic modelling perspective, the notion of a game in the work by Ghani et al. is not expressive enough for many applications. This includes stochastic environments, stochastic choices by players, as well as incomplete information regarding the game being played. The current paper addresses these three issues all at once.


title: Calculating Opacity in Hot, Dense Matter using Second-Order Electron-Photon and Two-Photon Transitions to Approximate Line Broadening

abstract: Calculations of the opacity of hot, dense matter require models for plasma line broadening. However, the most general theories are too complex to calculate directly and some approximation is inevitably required. The most widely-used approaches focus on the line centre, where a Lorentzian shape is obtained. Here, we demonstrate that in the opposite limit, far from the line centre, the opacity can be expressed in terms of second-order transitions, such as electron-photon and two-photon processes. We suggest that this insight could form the basis for a new approach to improve calculations of opacity in hot, dense matter. Preliminary calculations suggest that this approach could yield increased opacity away from absorption lines.


title: Simulations of particle tracking in the oligociliated mouse node and implications for left-right symmetry breaking mechanics

abstract: The concept of internal anatomical asymmetry is familiar; usually in humans the heart is on the left and the liver is on the right, however how does the developing embryo know to produce this consistent laterality? Symmetry breaking initiates with left-right asymmetric cilia-driven fluid mechanics in a small fluid-filled structure called the ventral node in mice. However the question of what converts this flow into left-right asymmetric development remains unanswered. A leading hypotheses is that flow transports morphogen containing vesicles within the node, the absorption of which results in asymmetrical gene expression. To investigate how vesicle transport might result in the situs patterns observed in wildtype and mutant experiments, we extend the open source Stokes flow package, NEAREST, to consider the hydrodynamic and Brownian motion of particles in a mouse model with flow driven by one, two, and 112 beating cilia. Three models for morphogen-containing particle released are simulated to assess their compatibility with observed results in oligociliated and wildtype mouse embryos: uniformly random release, localised cilium stress induced release, and localised release from motile cilia themselves. Only the uniformly random release model appears consistent with the data, with neither localised-release model resulting in significant transport in the oligociliated embryo.


title: Link Prediction Under Imperfect Detection: Collaborative Filtering for Ecological Networks

abstract: Matrix completion based collaborative filtering is considered scalable and effective for online service link prediction (e.g., movie recommendation) but does not meet the challenges of link prediction in ecological networks. A unique challenge of ecological networks is that the observed data are subject to systematic imperfect detection, due to the difficulty of accurate field sampling. In this work, we propose a new framework customized for ecological bipartite network link prediction. Our approach starts with incorporating the Poisson N-mixture model, a widely used framework in statistical ecology for modeling imperfect detection of a single species in field sampling. Despite its extensive use for single species analysis, this model has never been considered for link prediction between different species, perhaps because of the complex nature of both link prediction and N-mixture model inference. By judiciously combining the Poisson N-mixture model with a probabilistic nonnegative matrix factorization (NMF) model in latent space, we propose an intuitive statistical model for the problem of interest. We also offer a scalable and convergence-guaranteed optimization algorithm to handle the associated maximum likelihood identification problem. Experimental results on synthetic data and two real-world ecological networks data are employed to validate our proposed approach.


title: Penalized regression via the restricted bridge estimator

abstract: This article is concerned with the Bridge Regression, which is a special family in penalized regression with penalty function $\sum_{j=1}^{p}|\beta_j|^q$ with $q>0$, in a linear model with linear restrictions. The proposed restricted bridge (RBRIDGE) estimator simultaneously estimates parameters and selects important variables when a prior information about parameters are available in either low dimensional or high dimensional case. Using local quadratic approximation, the penalty term can be approximated around a local initial values vector and the RBRIDGE estimator enjoys a closed-form expression which can be solved when $q>0$. Special cases of our proposal are the restricted LASSO ($q=1$), restricted RIDGE ($q=2$), and restricted Elastic Net ($1< q < 2$) estimators. We provide some theoretical properties of the RBRIDGE estimator under for the low dimensional case, whereas the computational aspects are given for both low and high dimensional cases. An extensive Monte Carlo simulation study is conducted based on different prior pieces of information and the performance of the RBRIDGE estiamtor is compared with some competitive penalty estimators as well as the ORACLE. We also consider four real data examples analysis for comparison sake. The numerical results show that the suggested RBRIDGE estimator outperforms outstandingly when the prior is true or near exact


title: Intimal Growth in Cylindrical Arteries: Impact of Anisotropic Growth on Glagov Remodeling

abstract: In this paper we investigate the effect of anisotropic growth on Glagov remodeling in different cases: pure radial, pure circumferential, pure axial and general anisotropic growth. We use the theory of morphoelasticity on an axisymmetric arterial domain. For each case we explore their specific effect on the Glagov curves and stress and provide the changes in collagen fibers angles in the intima, media and adventitia. In addition, we compare the strain energy produced by growth in radial, circumferential and axial direction and deduce that anisotropic growth generally leads to lower strain energy than isotropic growth. Therefore, we explore an anisotropic growth regime and use the resulting model to simulate vessel remodeling. We compare the Glagov curves, stress, energies and fiber angles in the anisotropic case with those of the isotropic case. Our results show that the anisotropic growth produces a remodeling curve more consistent with Glagov's experimental data with gentler outward remodeling and more realistic stress profiles.


title: An Adversarial Model for Attack Vector Vulnerability Analysis on Power and Gas Delivery Operations

abstract: Power systems often rely on natural gas pipeline networks to supply fuel for gas-fired generation. Market inefficiencies and a lack of formal coordination between the wholesale power and gas delivery infrastructures may magnify the broader impact of a cyber-attack on a natural gas pipeline. In this study we present a model that can be used to quantify the impact of cyber-attacks on electricity and gas delivery operations. We model activation of cyber-attack vectors that attempt to gain access to pipeline gas compressor controls using a continuous-time Markov chain over a state space based on the gas operator Industrial Control System firewall zone partition. Our approach evaluates the operating states and decision-making in the networks using physically realistic and operationally representative models. We summarize these models, the sequence of analyses used to quantify the impacts of a cyber-incident, and propose a Monte Carlo simulation approach to quantify the resulting effect on the reliability of the bulk power system by the increase in operational cost. The methodology is applied to a case study of interacting power, gas, and cyber test networks.


title: Revisiting the Late-Time Growth of Single-mode Rayleigh-Taylor Instability and the Role of Vorticity

abstract: Growth of the single-fluid single-mode Rayleigh-Taylor instability (RTI) is revisited in 2D and 3D using fully compressible high-resolution simulations. We conduct a systematic analysis of the effects of perturbation Reynolds number ($Re_p$) and Atwood number ($A$) on RTI's late-time growth. Contrary to the common belief that single-mode RTI reaches a terminal bubble velocity, we show that the bubble re-accelerates when $Re_p$ is sufficiently large, consistent with [Ramaparabhu et al. 2006, Wei and Livescu 2012]. However, unlike in [Ramaparabhu et al. 2006], we find that for a sufficiently high $Re_p$, the bubble's late-time acceleration is persistent and does not vanish. Analysis of vorticity dynamics shows a clear correlation between vortices inside the bubble and re-acceleration. Due to symmetry around the bubble and spike (vertical) axes, the self-propagation velocity of vortices points in the vertical direction. If viscosity is sufficiently small, the vortices persist long enough to enter the bubble tip and accelerate the bubble [Wei and Livescu 2012]. A similar effect has also been observed in ablative RTI [Betti and Sanz 2006]. As the spike growth increases relative to that of the bubble at higher $A$, vorticity production shifts downward, away from the centerline and toward the spike tip. We modify the Betti-Sanz model for bubble velocity by introducing a vorticity efficiency factor $\eta=0.45$ to accurately account for re-acceleration caused by vorticity in the bubble tip. It had been previously suggested that vorticity generation and the associated bubble re-acceleration are suppressed at high $A$. However, we present evidence that if the large $Re_p$ limit is taken first, bubble re-acceleration is still possible. Our results also show that re-acceleration is much easier to occur in 3D than 2D, requiring smaller $Re_p$ thresholds.


title: Invariant Vector and Majorana Star Representations of Qutrit States

abstract: The Qutrit state density matrix is of order 3 and depends on 8 parameters in general case. Visualization of this 8-dimensional state space is practically impossible using 8-dimensional vectors commonly used. Recently a 3-dimensional vector representation of the Qutrit state space (also called Invariant Vector Representation) has been proposed [1]. In this work we present the time-evolution of the IVR vectors for the Qutrit Cascade or {\Xi}-model to emphasize the advantages of the IVR and also relate IVR vectors to well-known stars of the Majorana State Representation (MSR).


title: Canadian Astronomy on Maunakea: On Respecting Indigenous Rights

abstract: (Abridged) Canadian astronomy has, for decades, benefited from access to observatories and participating in international consortia on one of the best astronomical sites in the world: Maunakea. However, Maunakea is part of the unceded territory of the Native Hawaiian peoples and has always been of special significance to Hawaiian culture. The use of the summit and its science reserve has created tensions in the past decade, particularly with the development of the Thirty Meter Telescope. A meaningful and respectful response from the International astronomy community is still lacking. It is expected that the LRP 2020 will continue to support Canadian astronomy on Maunakea so a better official statement on the position and involvement of CASCA should be prepared. In this paper we present recommendations, based on the United Nation Declaration for the Rights of Indigenous Peoples, for the Canadian astronomical community to better support Indigenous rights on Maunakea and Hawaii while providing clear guidelines for the astronomical community to participate in activities conducted on Indigenous land. This framework is designed to motivate conversations with Indigenous communities regarding our place on Indigenous lands and our roles, and responsibilities toward the communities we are working with. Furthermore, we propose this framework as a basis for engaging with communities around the world regarding consent for astronomical facilities.


title: Hadamard renormalization for a charged scalar field

abstract: The Hadamard representation of the Green's function of a quantum field on a curved space-time is a powerful tool for computations of renormalized expectation values. We study the Hadamard form of the Feynman Green's function for a massive charged complex scalar field in an arbitrary number of space-time dimensions. Explicit expressions for the coefficients in the Hadamard parametrix are given for two, three and four space-time dimensions. We then develop the formalism for the Hadamard renormalization of the expectation values of the scalar field condensate, current and stress-energy tensor. These results will have applications in the computation of renormalized expectation values for a charged quantum scalar field on a charged black hole space-time, and hence in addressing issues such as the quantum stability of the inner horizon.


title: REFUGE Challenge: A Unified Framework for Evaluating Automated Methods for Glaucoma Assessment from Fundus Photographs

abstract: Glaucoma is one of the leading causes of irreversible but preventable blindness in working age populations. Color fundus photography (CFP) is the most cost-effective imaging modality to screen for retinal disorders. However, its application to glaucoma has been limited to the computation of a few related biomarkers such as the vertical cup-to-disc ratio. Deep learning approaches, although widely applied for medical image analysis, have not been extensively used for glaucoma assessment due to the limited size of the available data sets. Furthermore, the lack of a standardize benchmark strategy makes difficult to compare existing methods in a uniform way. In order to overcome these issues we set up the Retinal Fundus Glaucoma Challenge, REFUGE (\url{https://refuge.grand-challenge.org}), held in conjunction with MICCAI 2018. The challenge consisted of two primary tasks, namely optic disc/cup segmentation and glaucoma classification. As part of REFUGE, we have publicly released a data set of 1200 fundus images with ground truth segmentations and clinical glaucoma labels, currently the largest existing one. We have also built an evaluation framework to ease and ensure fairness in the comparison of different models, encouraging the development of novel techniques in the field. 12 teams qualified and participated in the online challenge. This paper summarizes their methods and analyzes their corresponding results. In particular, we observed that two of the top-ranked teams outperformed two human experts in the glaucoma classification task. Furthermore, the segmentation results were in general consistent with the ground truth annotations, with complementary outcomes that can be further exploited by ensembling the results.


title: The Poisson-Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients

abstract: We present the theory and implementation of a Poisson-Boltzmann implicit solvation model for electrolyte solutions. This model can be combined with arbitrary electronic structure methods that provide an accurate charge density of the solute. A hierarchy of approximations for this model includes a linear approximation for weak electrostatic potentials, finite size of the mobile electrolyte ions and a Stern-layer correction. Recasting the Poisson-Boltzmann equations into Euler-Lagrange equations then significantly simplifies the derivation of the free energy of solvation for these approximate models. The parameters of the model are then either fit directly to experimental observables, e.g. the finite ion size, or optimized for agreement with experimental results. Experimental data for this optimization is available in the form of Sechenov coefficients that describe the linear dependence of the salting-out effect of solutes with respect to the electrolyte concentration. In the final part we rationalize the qualitative disagreement of the finite ion size modification to the Poisson-Boltzmann model with experimental observations by taking into account the electrolyte concentration dependence of the Stern layer. A route towards a revised model that captures the experimental observations while including the finite ion size effects is then outlined. This implementation paves the way for the study of electrochemical and electrocatalytic processes of molecules and cluster models with accurate electronic structure methods.


title: On the feasibility of parsimonious variable selection for Hotelling's $T^2$-test

abstract: Hotelling's $T^2$-test for the mean of a multivariate normal distribution is one of the triumphs of classical multivariate analysis. It is uniformly most powerful among invariant tests, and admissible, proper Bayes, and locally and asymptotically minimax among all tests. Nonetheless, investigators often prefer non-invariant tests, especially those obtained by selecting only a small subset of variables from which the $T^2$-statistic is to be calculated, because such reduced statistics are more easily interpretable for their specific application. Thus it is relevant to ask the extent to which power is lost when variable selection is limited to very small subsets of variables, e.g. of size one (yielding univariate Student-$t^2$ tests) or size two (yielding bivariate $T^2$-tests). This study presents some evidence, admittedly fragmentary and incomplete, suggesting that in some cases no power may be lost over a wide range of alternatives.

categories: math.ST econ.EM q-bio.QM q-fin.ST stat.TH


title: Acoustokinetics: Crafting force landscapes from sound waves

abstract: Factoring the pressure field of a harmonic sound wave into its amplitude and phase profiles provides the foundation for an analytical framework for studying acoustic forces that not only provides novel insights into the forces exerted by specified sound waves, but also addresses the inverse problem of designing sound waves to implement desired force landscapes. We illustrate the benefits of this acoustokinetic framework through case studies of purely nonconservative force fields, standing waves, pseudo-standing waves, and tractor beams.


title: Large deviations for the KPZ equation from the KP equation

abstract: Recently, Quastel and Remenik \cite{QRKP} [arXiv:1908.10353] found a remarkable relation between some solutions of the finite time Kardar-Parisi-Zhang (KPZ) equation and the Kadomtsev-Petviashvili (KP) equation. Using this relation we obtain the large deviations at large time and at short time for the KPZ equation with droplet initial conditions, and at short time with half-Brownian initial conditions. It is consistent with previous results and allows to obtain sub-leading corrections, as well as results at intermediate time. In addition, we find that the appropriate generating function associated to the full Brownian initial condition also satisfies the KP equation. Finally, generating functions for some linear statistics of the Airy point process are also found to satisfy the KP property, and consequences are discussed.

categories: cond-mat.dis-nn cond-mat.stat-mech math-ph math.MP nlin.SI


title: Quantum Error Correction

abstract: Quantum error correction is a set of methods to protect quantum information--that is, quantum states--from unwanted environmental interactions (decoherence) and other forms of noise. The information is stored in a quantum error-correcting code, which is a subspace in a larger Hilbert space. This code is designed so that the most common errors move the state into an error space orthogonal to the original code space while preserving the information in the state. It is possible to determine whether an error has occurred by a suitable measurement and to apply a unitary correction that returns the state to the code space, without measuring (and hence disturbing) the protected state itself. In general, codewords of a quantum code are entangled states. No code that stores information can protect against all possible errors; instead, codes are designed to correct a specific error set, which should be chosen to match the most likely types of noise. An error set is represented by a set of operators that can multiply the codeword state. Quantum error correction is used to protect information in quantum communication (where quantum states pass through noisy channels) and quantum computation (where quantum states are transformed through a sequence of imperfect computational steps in the presence of environmental decoherence to solve a computational problem). In quantum computation, error correction is just one component of fault-tolerant design.


title: Progress towards practical qubit computation using approximate Gottesman-Kitaev-Preskill codes

abstract: Encoding a qubit in the continuous degrees of freedom of an oscillator is a promising path to error-corrected quantum computation. One advantageous way to achieve this is through Gottesman-Kitaev-Preskill (GKP) grid states, whose symmetries allow for the correction of any small continuous error on the oscillator. Unfortunately, ideal grid states have infinite energy, so it is important to find finite-energy approximations that are realistic, practical, and useful for applications. In the first half of this work we investigate the impact of imperfect GKP states on computational circuits independently of the physical architecture. To this end, we analyze the behaviour of the physical and logical content of normalizable GKP states through several figures of merit, employing a recently-developed modular subsystem decomposition. By tracking the errors that enter into the computational circuit due to imperfections in the GKP states, we are able to gauge the utility of these states for NISQ (Noisy Intermediate-Scale Quantum) devices. In the second half, we focus on a state preparation approach in the photonic domain wherein photon-number-resolving measurements on some modes of Gaussian states produce non-Gaussian states in others. We produce detailed numerical results for the preparation of GKP states alongside estimating the resource requirements in practical settings and probing the quality of the resulting states with the tools we develop. Our numerical experiments indicate that we can generate any state in the GKP Bloch sphere with nearly equal resources, which has implications for magic state preparation overheads.


title: Stone pseudovarieties

abstract: Profinite algebras are the residually finite compact algebras; their underlying topological spaces are Stone spaces. We extend the theory of profinite algebras to a more general setting of Stone topological algebras. We introduce Stone pseudovarieties, that is, classes of Stone topological algebras of a fixed topological signature that are closed under taking Stone quotients, closed subalgebras and finite direct products. Looking at Stone spaces as the dual spaces of Boolean algebras, we find a simple characterization of when the dual space admits a natural structure of topological algebra. This provides a new approach to duality theory which culminates in the proof that a Stone quotient of a Stone topological algebra that is residually in a given Stone pseudovariety is also residually in it, thereby extending the corresponding result of M. Gehrke for the Stone pseudovariety of all finite algebras over discrete signatures. The residual closure of a Stone pseudovariety is thus a Stone pseudovariety, and these are precisely the Stone analogues of varieties. A Birkhoff type theorem for Stone varieties is also established.


title: Estimands and Inference in Cluster-Randomized Vaccine Trials

abstract: Cluster-randomized trials are often conducted to assess vaccine effects. Defining estimands of interest before conducting a trial is integral to the alignment between a study's objectives and the data to be collected and analyzed. This paper considers estimands and estimators for overall, indirect, and total vaccine effects in trials where clusters of individuals are randomized to vaccine or control. The scenario is considered where individuals self-select whether to participate in the trial and the outcome of interest is measured on all individuals in each cluster. Unlike the overall, indirect, and total effects, the direct effect of vaccination is shown in general not to be estimable without further assumptions, such as no unmeasured confounding. An illustrative example motivated by a cluster-randomized typhoid vaccine trial is provided.


title: Representation Learning with Statistical Independence to Mitigate Bias

abstract: Presence of bias (in datasets or tasks) is inarguably one of the most critical challenges in machine learning applications that has alluded to pivotal debates in recent years. Such challenges range from spurious associations between variables in medical studies to the bias of race in gender or face recognition systems. Controlling for all types of biases in the dataset curation stage is cumbersome and sometimes impossible. The alternative is to use the available data and build models incorporating fair representation learning. In this paper, we propose such a model based on adversarial training with two competing objectives to learn features that have (1) maximum discriminative power with respect to the task and (2) minimal statistical mean dependence with the protected (bias) variable(s). Our approach does so by incorporating a new adversarial loss function that encourages a vanished correlation between the bias and the learned features. We apply our method to synthetic data, medical images (containing task bias), and a dataset for gender classification (containing dataset bias). Our results show that the learned features by our method not only result in superior prediction performance but also are unbiased. The code is available at https://github.com/QingyuZhao/BR-Net/.


title: Baseline predictions of elliptic flow and fluctuations at the RHIC Beam Energy Scan using response coefficients

abstract: Currently the RHIC Beam Energy Scan is exploring a new region of the Quantum Chromodynamic phase diagram at large baryon densities that approaches nuclear astrophysics regimes. This provides an opportunity to study relativistic hydrodynamics in a regime where the net conserved charges of baryon number, strangeness, and electric charge play a role, which will significantly change the theoretical approach to simulating the baryon-dense Quark-Gluon Plasma. Here we detail many of the important changes needed to adapt both initial conditions and the medium to baryon-rich matter. Then, we make baseline predictions for the elliptical flow and fluctuations based on extrapolating the physics at LHC and top RHIC energies to support future analyses of where and how the new baryon-dense physics causes these extrapolations to break down. First we compare eccentricities across beam energies, exploring their underlying assumptions; we find the the extrapolated initial state is predicted to be nearly identical to that at AuAu $\sqrt{s_{NN}}=200$ GeV. Then the final flow harmonic predictions are based on linear+cubic response. We discuss preliminary STAR results in order to determine the implications that they have for linear+cubic response coefficients at the lowest beam energy of AuAu $\sqrt{s_{NN}}=7$ GeV.


title: Unfolding the Structure of a Document using Deep Learning

abstract: Understanding and extracting of information from large documents, such as business opportunities, academic articles, medical documents and technical reports, poses challenges not present in short documents. Such large documents may be multi-themed, complex, noisy and cover diverse topics. We describe a framework that can analyze large documents and help people and computer systems locate desired information in them. We aim to automatically identify and classify different sections of documents and understand their purpose within the document. A key contribution of our research is modeling and extracting the logical and semantic structure of electronic documents using deep learning techniques. We evaluate the effectiveness and robustness of our framework through extensive experiments on two collections: more than one million scholarly articles from arXiv and a collection of requests for proposal documents from government sources.


title: Performance Impact of Memory Channels on Sparse and Irregular Algorithms

abstract: Graph processing is typically considered to be a memory-bound rather than compute-bound problem. One common line of thought is that more available memory bandwidth corresponds to better graph processing performance. However, in this work we demonstrate that the key factor in the utilization of the memory system for graph algorithms is not necessarily the raw bandwidth or even the latency of memory requests. Instead, we show that performance is proportional to the number of memory channels available to handle small data transfers with limited spatial locality. Using several widely used graph frameworks, including Gunrock (on the GPU) and GAPBS \& Ligra (for CPUs), we evaluate key graph analytics kernels using two unique memory hierarchies, DDR-based and HBM/MCDRAM. Our results show that the differences in the peak bandwidths of several Pascal-generation GPU memory subsystems aren't reflected in the performance of various analytics. Furthermore, our experiments on CPU and Xeon Phi systems demonstrate that the number of memory channels utilized can be a decisive factor in performance across several different applications. For CPU systems with smaller thread counts, the memory channels can be underutilized while systems with high thread counts can oversaturate the memory subsystem, which leads to limited performance. Finally, we model the potential performance improvements of adding more memory channels with narrower access widths than are found in current platforms, and we analyze performance trade-offs for the two most prominent types of memory accesses found in graph algorithms, streaming and random accesses.


title: Fingerprints of Majorana bound states in quantum-rings

abstract: In this work, we investigate a quantum-ring coupled to a topological superconductor, hosting Majorana bound states (MBSs). We study the MBSs effects over the spectrum and persistent current along the quantum-ring. To obtain physical quantities, we employ the Green's function formalism. We found that the presence of MBSs leads to dramatic changes in the quantum-ring spectrum by inducing particle-hole symmetry. On the other hand, we have obtained a substantial reduction of the persistent current in the strong coupling limit between ring and MBSs depending on the ring parity. These behaviors provide a measurable signature of MBSs. Our findings could be used as additional information on MBSs presence in these quantum systems.


title: Large time behavior of solutions to Schr\"odinger equation with complex-valued potential

abstract: We study the large-time behavior of the solutions to the Schr\"odinger equation associated with a quickly decaying potential in dimension three. We establish the resolvent expansions at threshold zero and near positive resonances. The large-time expansions of solutions are obtained under different conditions, including the existence of positive resonances and zero resonance or/and zero eigenvalue.


title: Lippmann-Schwinger equation and the connection between the scattering operator and the scattering amplitude in the relativistic case

abstract: In this paper, we consider two types of the scattering problems (relativistic case), namely, the stationary scattering problem, where the distance $r$ tends to infinity, and the dynamical scattering problem, where the time $t$ tends to infinity. Using our results on Lippmann-Schwinger equation in the relativistic case, we found the connection between the stationary scattering problem (the scattering amplitude) and the dynamical scattering problem (the scattering operator). This result is the quantum mechanical analog of the ergodic formulas in the classical mechanics.


title: Schwinger-based QCD formulation's derivation of elastic pp scattering

abstract: Using previously described functional techniques for some non-perturbative, gauge invariant, renormalized QCD processes, a simplified version of the amplitudes - in which forms akin to Pomerons naturally appear - provides fits to ISR and LHC-TOTEM pp elastic scattering data. Those amplitudes rely on a specific function $\phi(b)$ which describes the fluctuations of the transverse position of quarks inside hadrons.


title: On the sensitivity of the optimal partition for parametric second-order conic optimization

abstract: In this paper, using an optimal partition approach, we study the parametric analysis of a second-order conic optimization problem, where the objective function is perturbed along a fixed direction. We characterize the notions of so-called invariancy set and nonlinearity interval, which serve as stability regions of the optimal partition. We then propose, under the strict complementarity condition, an iterative procedure to compute a nonlinearity interval of the optimal partition. Furthermore, under primal and dual nondegeneracy conditions, we show that a boundary point of a nonlinearity interval can be numerically identified from a nonlinear reformulation of the parametric second-order conic optimization problem. Our theoretical results are supported by numerical experiments.


title: An Exact and Scalable Problem Decomposition for Security-Constrained Optimal Power Flow

abstract: In this paper, we present decomposition techniques for solving large-scale instances of the security-constrained optimal power flow (SCOPF) problem with primary response. Specifically, under each contingency state, we require that the nodal demands are met and that the synchronized units generating below their limits follow a linear model for primary response. The resulting formulation is a mixed-integer linear program since the primary response model introduces disjunctions to the SCOPF problem. Unfortunately, exact methods relying on traditional Benders' decomposition do not scale well. As an alternative, we propose a decomposition scheme based on the column-and-constraint-generation algorithm where we iteratively add disjunctions and cuts. We provide procedures for preprocessing dedicated cuts and for numerically determining the post-contingency responses based on the master problem solutions. We also discuss heuristics to generate high-quality primal solutions and upper bounds for the method. Finally, we demonstrate the efficiency of the proposed method on large-scale systems.


title: New Problems and Solutions in IoT Security and Privacy

abstract: In a previous article for S&P magazine, we made a case for the new intellectual challenges in the Internet of Things security research. In this article, we revisit our earlier observations and discuss a few results from the computer security community that tackle new issues. Using this sampling of recent work, we identify a few broad general themes for future work.


title: Singular Perturbation-based Large-Signal Order Reduction of Microgrids for Stability and Accuracy Synthesis with Control

abstract: With the increasing penetration of distributed energy resources (DERs), it is of vital importance to study the dynamic stability of microgrids (MGs) with external control inputs in the electromagnetic transient (EMT) time scale. This requires detailed models of the underlying control structure of MGs and results in a high-order nonlinear MG control system. Higher-level controller design and stability analysis of such high-order systems are usually intractable and computation-costly. To overcome these challenges, this paper proposes a large-signal order reduction (LSOR) method for MGs with considerations of external control inputs and the detailed dynamics of underlying control levels based on singular perturbation theory (SPT). Specially, we innovatively proposed and strictly proved a general stability and accuracy assessment theorem that allows us to analyze the dynamic stability of a full-order nonlinear system by only leveraging its corresponding reduced-order model (ROM) and boundary layer model (BLM). Moreover, this theorem also theoretically provides a set of conditions under which the developed ROM is accurate. Finally, by embedding such a theorem into the SPT, we propose a novel LSOR approach with guaranteed accuracy and stability analysis equivalence. The proposed LSOR method is generic and can be applied to arbitrary dynamic systems. Multiple case studies are conducted on MG systems to show the effectiveness of the proposed approach.


title: Investigating the Temperature Distribution of Diatomic Carbon in Comets using the Swan Bands

abstract: We present high spectral-resolution observations of comets 122P/de Vico and 153P/Ikeya-Zhang obtained with the Tull Coud\'{e} spectrograph on the 2.7m Harlan J. Smith telescope of McDonald Observatory. We used these data to study the distribution of the lines of the $\mathrm{d} ^3\Pi_g - \mathrm{a} ^3\Pi_u$ C$_2$ (Swan) bands. We show that the data are best represented with two rotational temperatures, with the lowest energy lines being at a relatively cool temperature and the higher energy lines being at a higher temperature. We discuss the implications of this two temperature distribution and suggest future work.


title: Regularity for minimizers of a class of non-autonomous functionals with sub-quadratic growth

abstract: We consider a class of integral functionals with convex integrand with respect to the gradient variable, assuming that the function that measures the oscillation of the integrand with respect to the x variable belongs to a suitable Sobolev space W^{1;q}. We prove a result of higer differentiability for the minimizers. We also infer a result of Lipschitz regularity of minimizers if q > n, and a result of higher integrability for the gradient if q = n. The novelty here is that we deal with integrands satisfying subquadratic growth conditions with respect to gradient variable.


title: The Semicontinuity of Attractors for Closed Relations on Compact Hausdorff Spaces

abstract: We show that attractors are semicontinuous for closed relations on compact Hausdorff spaces. Semicontinuity is what guarantees that small changes to a system do not result in massive growth of certain features, notably attractors. That is, there is a certain preservation of structure. When it comes to flows, semiflows, and maps, it is well established that attractors are semicontinuous. In [2], relations were established as a way to generalize maps, and a formal definition of attractors was established. Relations (in the dynamical systems sense) represent discrete time systems, which may lack uniqueness (or existence) in forward time.


title: Time optimal observability for Grushin Schr\"odinger equation

abstract: We consider two dimensional Grushin Schr\"odinger equation posed on a finite cylinder $\Omega=(-1,1)_x\times \T_y$ with Dirichlet boundary condition. We obtain the sharp observability by any horizontal strip, with the optimal time $T_*>0$ depending on the size of the strip. Consequently, we prove the exact controllability for Grushin Schr\"odinger equation. By exploiting the concentration of eigenfunctions of harmonic oscillator at $x=0$, we also show that the observability fails for any $T\leq T_*$.


title: Balanced Viscosity solutions to a rate-independent coupled elasto-plastic damage system

abstract: A rate-independent model coupling small strain associative elasto-plasticity and damage is studied via a 'vanishing-viscosity' analysis with respect to all the variables describing the system. This extends the analysis performed for the same system in [Crismale-Lazzaroni 2016], where a vanishing-viscosity regularization involving only the damage variable was set forth. In the present work, an additional approximation featuring vanishing plastic hardening is introduced in order to deal with the vanishing viscosity in the plastic variable. Different regimes are considered, leading to different notions of Balanced Viscosity solutions for the perfectly plastic damage system, and for its version with hardening.


title: A dynamically load-balanced parallel $ p $-adaptive implicit high-order flux reconstruction method for under-resolved turbulence simulation

abstract: We present a dynamically load-balanced parallel $ p $-adaptive implicit high-order flux reconstruction method for under-resolved turbulence simulation. The high-order explicit first stage, singly diagonal implicit Runge-Kutta (ESDIRK) method is employed to circumvent the restriction on the time step size. The pseudo transient continuation is coupled with the matrix-free restarted generalized minimal residual (GMRES) method to solve the nonlinear equations at each stage, except the first one, of ESDIRK. We use the spectral decay smoothness indicator as the refinement/coarsening indicator for $ p $-adaptation. A dynamic load balancing technique is developed with the aid of the open-source library ParMETIS. The trivial cost, compared to implicit time stepping, of mesh repartitioning and data redistribution enables us to conduct $ p $-adaptation and load balancing every time step. An isentropic vortex propagation case is employed to study the impact of element weights used in mesh repartitioning on parallel efficiency. We apply the $ p$-adaptive solver for implicit large eddy simulation (ILES) of the transitional flows over a cylinder when Reynolds number (Re) is $ 3900 $ and the SD7003 wing when Re is $ 60000 $. Numerical experiments demonstrate that a significant reduction in the run time (up to $70\%$) and total number of solution points (up to $76\%$) can be achieved with $ p $-adaptation.


title: Breakdown of the equivalence between gravitational mass and energy due to quantum effects

abstract: We review our recent theoretical results about inequivalence between passive and active gravitational masses and energy in semiclassical variant of general relativity, where gravitational field is not quantized but matter is quantized. To this end, we consider the simplest quantum body with internal degrees of freedom - a hydrogen atom. We concentrate our attention on the following physical effects, related to electron mass. The first one is inequivalence between passive gravitational mass and energy at microscopic level. Indeed, quantum measurement of gravitational mass can give result, which is different from the expected, $m \neq m_e + \frac{E_1}{c^2}$, where electron is initially in its ground state; $m_e$ is the bare electron mass. The second effect is that the expectation values of both passive and active gravitational masses of stationary quantum states are equivalent to the expectation value of energy. The most spectacular effects are inequivalence of passive and active gravitational masses and energy at macroscopic level for ensemble of coherent superpositions of stationary quantum states. We show that, for such superpositions, the expectation values of passive and active gravitational masses are not related to the expectation value of energy by the famous Einstein's equation, $m \neq \frac{E}{c^2}$. In this review, we also improve several drawbacks of the original pioneering works.


title: NADS-Net: A Nimble Architecture for Driver and Seat Belt Detection via Convolutional Neural Networks

abstract: A new convolutional neural network (CNN) architecture for 2D driver/passenger pose estimation and seat belt detection is proposed in this paper. The new architecture is more nimble and thus more suitable for in-vehicle monitoring tasks compared to other generic pose estimation algorithms. The new architecture, named NADS-Net, utilizes the feature pyramid network (FPN) backbone with multiple detection heads to achieve the optimal performance for driver/passenger state detection tasks. The new architecture is validated on a new data set containing video clips of 100 drivers in 50 driving sessions that are collected for this study. The detection performance is analyzed under different demographic, appearance, and illumination conditions. The results presented in this paper may provide meaningful insights for the autonomous driving research community and automotive industry for future algorithm development and data collection.


title: A preliminary study on a multi-resolution-level inverse planning algorithm for Gamma Knife radiosurgery

abstract: Manual forward planning for GK radiosurgery is complicated and time-consuming, particularly for cases with large or irregularly shaped targets. Inverse planning eases GK planning by solving an optimization problem. However, due to the vast search space, most inverse planning algorithms have to decouple the planning process to isocenter preselection and sector duration optimization. This sequential scheme does not necessarily lead to optimal isocenter locations and hence optimal plans. In this study, we attempt to optimize the isocenter positions, beam shapes and durations simultaneously by proposing a multi-resolution-level (MRL) strategy to handle the large-scale GK optimization problem. In our approach, several rounds of optimizations were performed with a progressively increased spatial resolution for isocenter candidate selection. The isocenters selected from last round and their neighbors on a finer resolution were used as new isocenter candidates for next round of optimization. After plan optimization, shot sequencing was performed to group the optimized sectors to deliverable shots supported by GK treatment units. We have tested our algorithm on 6 GK cases previously treated in our institution (2 meningioma cases, 3 cases with single metastasis and 1 case with 6 metastases). Compared with manual planning, achieving same coverage and similar selectivity, our algorithm improved the gradient index averagely from 3.1 to 2.9 and reduced the maximum dose of brainstem from 8.0Gy to 5.6Gy. The average beam-on time was also reduced by from 103.8 mins to 87.4 mins. Our method was also compared with the inverse planning algorithm provided in Leksell GammaPlan planning system, and outperformed it with better plan quality for all the 6 cases.This preliminary study has demonstrated the effectiveness and feasibility of our MRL inverse planning approach for GK radiosurgery.


title: An epitaxial graphene platform for zero-energy edge state nanoelectronics

abstract: Graphene's original promise to succeed silicon faltered due to pervasive edge disorder in lithographically patterned deposited graphene and the lack of a new electronics paradigm. Here we demonstrate that the annealed edges in conventionally patterned graphene epitaxially grown on a silicon carbide substrate (epigraphene) are stabilized by the substrate and support a protected edge state. The edge state has a mean free path that is greater than 50 microns, 5000 times greater that the bulk states and involves a theoretically unexpected Majorana-like zero-energy non-degenerate quasiparticle that does not produce a Hall voltage. In seamless integrated structures, the edge state forms a zero-energy one-dimensional ballistic network with essentially dissipationless nodes at ribbon-ribbon junctions. Seamless device structures offer a variety of switching possibilities including quantum coherent devices at low temperatures. This makes epigraphene a technologically viable graphene nanoelectronics platform that has the potential to succeed silicon nanoelectronics.


title: AutoML using Metadata Language Embeddings

abstract: As a human choosing a supervised learning algorithm, it is natural to begin by reading a text description of the dataset and documentation for the algorithms you might use. We demonstrate that the same idea improves the performance of automated machine learning methods. We use language embeddings from modern NLP to improve state-of-the-art AutoML systems by augmenting their recommendations with vector embeddings of datasets and of algorithms. We use these embeddings in a neural architecture to learn the distance between best-performing pipelines. The resulting (meta-)AutoML framework improves on the performance of existing AutoML frameworks. Our zero-shot AutoML system using dataset metadata embeddings provides good solutions instantaneously, running in under one second of computation. Performance is competitive with AutoML systems OBOE, AutoSklearn, AlphaD3M, and TPOT when each framework is allocated a minute of computation. We make our data, models, and code publicly available.


title: Structural, Rheological and Dynamic Aspects of Hydrogen-Bonding Molecular Liquids: Aqueous Solutions of Hydrotropic tert-Butyl Alcohol

abstract: Hypothesis: The structural details, viscosity trends and dynamic phenomena in t-butanol/water solutions are closely related on the molecular scales across the entire composition range. Utilizing the experimental small- and wide-angle x-ray scattering (SWAXS) method, molecular dynamics (MD) simulations and the complemented-system approach method developed in our group it is possible to comprehensively describe the structure-viscosity-dynamics relationship in such structurally versatile hydrogen-bonded molecular liquids, as well as in similar, self-assembling systems with pronounced molecular and supramolecular structures at the intra-, inter-, and supra-molecular scales. Experiments: The SWAXS and x-ray diffraction experiments and MD simulations were performed for aqueous t-butanol solutions at 25 {\deg}C. Literature viscosity and self-diffusion data were also used. Findings: The interpretive power of the proposed scheme was demonstrated by the extensive and diverse results obtained for aqueous t-butanol solutions across the whole concentration range. Four composition ranges with qualitatively different structures and viscosity trends were revealed. The experimental and calculated zero-shear viscosities and molecular self-diffusion coefficients were successfully related to the corresponding structural details. The hydrogen bonds, that were, along with hydrophobic effects, recognized as the most important driving force for the formation of t-butanol aggregates, show intriguing lifetime trends and thermodynamic properties of their formation.


title: Scalable Memdiodes Exhibiting Rectification and Hysteresis for Neuromorphic Computing

abstract: Metal-Nb$_{2}$O$_{5-x}$-metal memdiodes exhibiting rectification, hysteresis, and capacitance are demonstrated for applications in neuromorphic circuitry. These devices do not require any post-fabrication treatments such as filament creation by electroforming that would impede circuit scalability. Instead these devices operate due to Poole-Frenkel defect controlled transport where the high defect density is inherent to the Nb$_{2}$O$_{5-x}$ deposition rather than post-fabrication treatments. Temperature dependent measurements reveal that the dominant trap energy is 0.22 eV suggesting it results from the oxygen deficiencies in the amorphous Nb$_{2}$O$_{5-x}$. Rectification occurs due to a transition from thermionic emission to tunneling current and is present even in thick devices (> 100 nm) due to charge trapping which controls the tunneling distance. The turn-on voltage is linearly proportional to the Schottky barrier height and, in contrast to traditional metal-insulator-metal diodes, is logarithmically proportional to the device thickness. Hysteresis in the I-V curve occurs due to the current limited filling of traps.