title: Power spectrum of turbulent convection in the solar photosphere

abstract: The solar photosphere provides us with a laboratory for understanding turbulence in a layer where the fundamental processes of transport vary rapidly and a strongly superadiabatic region lies very closely to a subadiabatic layer. Our tools for probing the turbulence are high-resolution spectropolarimetric observations such as have recently been obtained with the two sunrise missions, and numerical simulations. Our aim is to study photospheric turbulence with the help of Fourier power spectra that we compute from observations and simulations. We also attempt to explain some properties of the photospheric overshooting flow with the help of its governing equations and simulations. We find that quiet-Sun observations and smeared simulations exhibit a power-law behavior in the subgranular range of their Doppler velocity power spectra with an index of$~\approx -2$. The unsmeared simulations exhibit a power-law index of$~\approx -2.25$. The smearing considerably reduces the extent of the power-law-like portion of the spectra. Therefore, the limited spatial resolution in some observations might eventually result in larger uncertainties in the estimation of the power-law indices. The simulated vertical velocity power spectra as a function of height show a rapid change in the power-law index from the solar surface to $300$~km above it. A scale-dependent transport of the vertical momentum occurs. At smaller scales, the vertical momentum is more efficiently transported sideways than at larger scales. This results in less vertical velocity power transported upward at small scales than at larger scales and produces a progressively steeper vertical velocity power law below $180$ km. Above this height, the gravity work progressively gains importance at all scales, making the atmosphere progressively more hydrostatic and resulting in a gradually less steep power law.


title: Backscattering in Nonlinear Microring Resonators Via A Gaussian Treatment of Coupled Cavity Modes

abstract: Systems of coupled cavity modes have the potential to provide bright quantum optical states of light in a highly versatile manner. Microring resonators for instance are highly scalable candidates for photon sources thanks to CMOS fabrication techniques, their small footprint and the relative ease of coupling many such microrings together, however, surface roughness of the wave-guides, and defects in the coupler geometry routinely induce splitting of the cavity modes due to backscattering and backcoupling. The parasitic back-propagating mode in the microring leads to hybridisation of the modes, altering the linear and nonlinear properties of this system of coupled cavity modes, and ultimately constraining the fidelity of quantum light sources that can be produced. In this paper, we derive a comprehensive general model for Gaussian nonlinear processes in systems of coupled cavity modes, based on an effective field Hamiltonian and a dispersive input-output model. The resulting dynamics of the equations of motion are evaluated in a Gaussian process formalism via the symplectic transformations on the optical modes. We then use this framework to numerically model and explore the problem of backscattering in microring resonators in physically relevant parameter regimes, involving the splitting of various resonances, we calculate the consequent impurity and heralding efficiency of various heralded photon schemes, we explore a perturbative explanation of the observations and assess the correspondence between spontaneous and stimulated processes in these systems.


title: On licenses for [Open] Hardware

abstract: This document explains the basic concepts related to software and hardware licenses, and it summarizes the most popular licenses that are currently used for hardware projects. Two case studies of hardware projects at different levels of abstraction are also presented, together with a discussion of license applicability, commercial issues, code protection, and related concerns. This paper intends to help the reader understand how to release open hardware with the most appropriate license, and to answer questions that are of current interest. We have been mainly motivated by the growing influence of the open RISC-V ISA, but trying to address a wider hardware point of view.


title: Studying Visualization Guidelines According to Grounded Theory

abstract: Visualization guidelines, if defined properly, are invaluable to both practical applications and the theoretical foundation of visualization. In this paper, we present a collection of research activities for studying visualization guidelines according to Grounded Theory (GT). We used the discourses at VisGuides, which is an online discussion forum for visualization guidelines, as the main data source for enabling data-driven research processes as advocated by the grounded theory methodology. We devised a categorization scheme focusing on observing how visualization guidelines were featured in different threads and posts at VisGuides, and coded all 248 posts between September 27, 2017 (when VisGuides was first launched) and March 13, 2019. To complement manual categorization and coding, we used text analysis and visualization to help reveal patterns that may have been missed by the manual effort and summary statistics. To facilitate theoretical sampling and negative case analysis, we made an in-depth analysis of the 148 posts (with both questions and replies) related to a student assignment of a visualization course. Inspired by two discussion threads at VisGuides, we conducted two controlled empirical studies to collect further data to validate specific visualization guidelines. Through these activities guided by grounded theory, we have obtained some new findings about visualization guidelines.


title: Evaluation of a Vision-to-Audition Substitution System that Provides 2D WHERE Information and Fast User Learning

abstract: Vision to audition substitution devices are designed to convey visual information through auditory input. The acceptance of such systems depends heavily on their ease of use, training time, reliability and on the amount of coverage of online auditory perception of current auditory scenes. Existing devices typically require extensive training time or complex and computationally demanding technology. The purpose of this work is to investigate the learning curve for a vision to audition substitution system that provides simple location features. Forty-two blindfolded users participated in experiments involving location and navigation tasks. Participants had no prior experience with the system. For the location task, participants had to locate 3 objects on a table after a short familiarisation period (10 minutes). Then once they understood the manipulation of the device, they proceeded to the navigation task: participants had to walk through a large corridor without colliding with obstacles randomly placed on the floor. Participants were asked to repeat the task 5 times. In the end of the experiment, each participant had to fill out a questionnaire to provide feedback. They were able to perform localisation and navigation effectively after a short training time with an average of 10 minutes. Their navigation skills greatly improved across the trials.


title: Addressing Variance Shrinkage in Variational Autoencoders using Quantile Regression

abstract: Estimation of uncertainty in deep learning models is of vital importance, especially in medical imaging, where reliance on inference without taking into account uncertainty could lead to misdiagnosis. Recently, the probabilistic Variational AutoEncoder (VAE) has become a popular model for anomaly detection in applications such as lesion detection in medical images. The VAE is a generative graphical model that is used to learn the data distribution from samples and then generate new samples from this distribution. By training on normal samples, the VAE can be used to detect inputs that deviate from this learned distribution. The VAE models the output as a conditionally independent Gaussian characterized by means and variances for each output dimension. VAEs can therefore use reconstruction probability instead of reconstruction error for anomaly detection. Unfortunately, joint optimization of both mean and variance in the VAE leads to the well-known problem of shrinkage or underestimation of variance. We describe an alternative approach that avoids this variance shrinkage problem by using quantile regression. Using estimated quantiles to compute mean and variance under the Gaussian assumption, we compute reconstruction probability as a principled approach to outlier or anomaly detection. Results on simulated and Fashion MNIST data demonstrate the effectiveness of our approach. We also show how our approach can be used for principled heterogeneous thresholding for lesion detection in brain images.


title: Berkovich curves and Schottky uniformization

abstract: This text is an exposition of non-Archimedean curves and Schottky uniformization from the point of view of Berkovich geometry. It consists of two parts, the first one of an introductory nature, and the second one more advanced. The first part is meant to be an introduction to the theory of Berkovich spaces focused on the case of the affine line. We define the Berkovich affine line and present its main properties, with many details: classification of points, path-connectedness, metric structure, variation of rational functions, etc. Contrary to many other introductory texts, we do not assume that the base field is algebraically closed. The second part is devoted to the theory of Mumford curves and Schottky uniformization. We start by briefly reviewing the theory of Berkovich curves, then introduce Mumford curves in a purely analytic way (without using formal geometry). We define Schottky groups acting on the Berkovich projective line, highlighting how geometry and group theory come together to prove that the quotient by the action of a Schottky group is an analytic Mumford curve. Finally, we present an analytic proof of Schottky uniformization, showing that any analytic Mumford curves can be described as a quotient of this kind. The guiding principle of our exposition is to stress notions and fully prove results in the theory of non-Archimedean curves that, to our knowledge, are not fully treated in other texts.


title: Attenuation of electromagnetic radiation in Nuclear Track Detectors

abstract: A systematic study of the attenuation of electromagnetic radiation in Nuclear Track Detectors (NTDs) is carried out. The attenuation of gamma-ray, X-ray, UV, visible, and infrared radiation in NTDs are investigated using NaI(Tl) detector, Gas Electron Multiplier (GEM) detector, UV-Vis spectrophotometer, and FTIR spectrophotometer respectively. The values of some important parameters (e.g., optical depth, attenuation coefficient, etc.) of three commercially available NTDs (PET, Makrofol r and CR-39 r ), at the relevant region of the electromagnetic spectrum, is determined. The details of the experimental techniques and the results are also presented in this paper.


title: On Threshold Solutions of equivariant Chern-Simons-Schr\"odinger Equation

abstract: We consider the self-dual Chern-Simons-Schr\"odinger model in two spatial dimensions. This problem is $L^2$-critical. Under equivariant setting, global wellposedness and scattering were proved in [Liu-Smith, 2016] for solution with initial charge below certain threshold given by the ground state. In this work, we show that the only non-scattering solutions with threshold charge are exactly the ground state up to scaling, phase rotation and the pseudoconformal transformation. We also obtain partial result for non-self-dual system.


title: Unsupervised Neural Machine Translation for Low-Resource Domains via Meta-Learning

abstract: Unsupervised machine translation, which utilizes unpaired monolingual corpora as training data, has achieved comparable performance against supervised machine translation. However, it still suffers from data-scarce domains. To address this issue, this paper presents a novel meta-learning algorithm for unsupervised neural machine translation (UNMT) that trains the model to adapt to another domain by utilizing only a small amount of training data. We assume that domain-general knowledge is a significant factor in handling data-scarce domains. Hence, we extend the meta-learning algorithm, which utilizes knowledge learned from high-resource domains, to boost the performance of low-resource UNMT. Our model surpasses a transfer learning-based approach by up to 2-4 BLEU scores. Extensive experimental results show that our proposed algorithm is pertinent for fast adaptation and consistently outperforms other baseline models.


title: Inverse Dynamics Control of Compliant Hybrid Zero Dynamic Walking

abstract: We present a trajectory planning and control architecture for bipedal locomotion at a variety of speeds on a highly underactuated and compliant bipedal robot. A library of compliant walking trajectories are planned offline, and stored as compact arrays of polynomial coefficients for tracking online. The control implementation uses a floating-base inverse dynamics controller which generates dynamically consistent feedforward torques to realize walking using information obtained from the trajectory optimization. The effectiveness of the controller is demonstrated in simulation and on hardware for walking both indoors on flat terrain and over unplanned disturbances outdoors. Additionally, both the controller and optimization source code are made available on GitHub.


title: Vacuum semiclassical gravity does not leave space for safe singularities

abstract: General relativity predicts its own demise at singularities, but also appears to conveniently shield itself from the catastrophic consequences of such singularities, making them safe. For instance, if strong cosmic censorship were ultimately satisfied, spacetime singularities, although present, would not pose any practical problems to predictability. Here we argue that under semiclassical effects the situation should be rather different: the potential singularities which could appear in the theory will generically affect predictability and so one will be forced to analyse whether there is a way to regularise them. For these possible regularisations, the presence and behaviour of matter during gravitational collapse and stabilisation into new structures will play a key role. First we show that the static semiclassical counterparts to the Schwarzschild and Reissner-Nordstr\"om geometries have singularities which are no longer hidden behind horizons. Then we argue that in dynamical scenarios of formation and evaporation of black holes, we are left with only three possible outcomes which could avoid singularities and eventual predictability issues. We briefly analyse the viability of each one of them within semiclassical gravity, and discuss the expected characteristic timescales of their evolution.


title: Testing the Quasar Hubble Diagram with LISA Standard Sirens

abstract: Quasars have recently been used as an absolute distance indicator, extending the Hubble diagram to high redshift to reveal a deviation from the expansion history predicted for the standard, $\Lambda$CDM cosmology. Here we show that the Laser Interferometer Space Antenna (LISA) will efficiently test this claim with standard sirens at high redshift, defined by the coincident gravitational wave (GW) and electromagnetic (EM) observations of the merger of massive black hole binaries (MBHBs). Assuming a fiducial $\Lambda$CDM cosmology for generating mock standard siren datasets, the evidence for the $\Lambda$CDM model with respect to an alternative model inferred from quasar data is investigated. By simulating many realizations of possible future LISA observations, we find that for $50\%$ of these realizations (median result) 4 MBHB standard siren measurements will suffice to strongly differentiate between the two models, while 14 standard sirens will yield a similar result in $95\%$ of the realizations. In addition, we investigate the measurement precision of cosmological parameters as a function of the number of observed LISA MBHB standard sirens, finding that 15 events will on average achieve a relative precision of 5\% for $H_0$, reducing to 3\% and 2\% with 25 and 40 events, respectively. Our investigation clearly highlights the potential of LISA as a cosmological probe able to accurately map the expansion of the universe at $z\gtrsim 2$, and as a tool to cross-check and cross-validate cosmological EM measurements with complementary GW observations.


title: Driving enhanced quantum sensing in partially accessible many-body systems

abstract: The Ground-state criticality of many-body systems is a resource for quantum-enhanced sensing, namely the Heisenberg precision limit, provided that one has access to the whole system. We show that for partial accessibility, the sensing capabilities of a block of spins in the ground state reduces to the sub-Heisenberg limit. To compensate for this, we drive the hamiltonian periodically and use a local steady-state for quantum sensing. Remarkably, the steady-state sensing shows a significant enhancement in precision compared to the ground state and even achieves super-Heisenberg scaling for low frequencies. The origin of this precision enhancement is related to the closing of the Floquet quasienergy gap. It is in close correspondence with the vanishing of the energy gap at criticality for ground state sensing with global accessibility. The proposal is general to all the integrable models and can be implemented on existing quantum devices.


title: Long time universality of black-hole lasers

abstract: For flowing quantum gases, it has been found that at long times an initial black-hole laser (BHL) configuration exhibits only two possible states: the ground state or a periodic self-oscillating state of continuous emission of solitons. So far, all the works on this subject are based on a highly idealized model, quite difficult to implement experimentally. Here we study the instability spectrum and the time evolution of a recently proposed realistic model of a BHL, thus providing a useful theoretical tool for the clear identification of black-hole lasing in future experiments. We further confirm the existence of a well-defined phase diagram at long times, which bespeaks universality in the long-time behavior of a BHL. Additionally, we develop a complementary model in which the same potential profile is applied to a subsonic homogeneous flowing condensate that, despite not forming a BHL, evolves towards the same phase diagram as the associated BHL model. This result reveals an even stronger form of robustness in the long-time behavior with respect to the transient, which goes beyond what has been described in the previous literature.


title: Diagrammatic study of optical excitations in correlated systems

abstract: The optical conductivity contains relevant information on the properties of correlated electron systems. In infinite dimensions, where dynamical mean field theory becomes exact, vertex corrections can be neglected and the conductivity computed from particle-hole bubbles. An interesting question concerns the nature and effect of the most relevant vertex corrections in finite-dimensional systems. A recent numerical study showed that the dominant vertex correction near an ordering instability with wave vector {\pi} comes from a vertical ladder, analogous to the Maki-Thompson diagram. Since the RPA version of this ladder diagram, dubbed {\pi}-ton, can be easily evaluated, this suggests a simple procedure for incorporating antiferromagnetic or charge density wave fluctuations into dynamical mean field estimates of the optical conductivity and related susceptibilities. We implement this procedure for the half-filled Hubbard model, considering the {\pi}-ton and a double-ladder extension of the {\pi}-ton, and reveal the spectral signatures of these vertex corrections.


title: On regularization of the Heun functions

abstract: In the paper we consider the Heun functions, which are solutions of the equation introduced by Karl Heun in 1889. The Heun functions generalize many known special functions and appear in many fields of modern physics. Evaluation of the functions was described in [1]. It is based on local power series solutions near the origin, derived by the Frobenius method, and analytic continuation to the whole complex plane with branch cuts. However, exceptional cases can occur at integer values of an exponent-related parameter $\gamma$ of the equation, when one of the two local solutions should include a logarithmic term. This also means singular behavior of the Heun functions as $\gamma$ approaches the integer values. Here we suggest a method of regularization and redefine the Heun functions in some vicinities of the integer values of $\gamma$, where the new functions depend smoothly on the parameter.


title: Model-free conventions in multi-agent reinforcement learning with heterogeneous preferences

abstract: Game theoretic views of convention generally rest on notions of common knowledge and hyper-rational models of individual behavior. However, decades of work in behavioral economics have questioned the validity of both foundations. Meanwhile, computational neuroscience has contributed a modernized 'dual process' account of decision-making where model-free (MF) reinforcement learning trades off with model-based (MB) reinforcement learning. The former captures habitual and procedural learning while the latter captures choices taken via explicit planning and deduction. Some conventions (e.g. international treaties) are likely supported by cognition that resonates with the game theoretic and MB accounts. However, convention formation may also occur via MF mechanisms like habit learning; though this possibility has been understudied. Here, we demonstrate that complex, large-scale conventions can emerge from MF learning mechanisms. This suggests that some conventions may be supported by habit-like cognition rather than explicit reasoning. We apply MF multi-agent reinforcement learning to a temporo-spatially extended game with incomplete information. In this game, large parts of the state space are reachable only by collective action. However, heterogeneity of tastes makes such coordinated action difficult: multiple equilibria are desirable for all players, but subgroups prefer a particular equilibrium over all others. This creates a coordination problem that can be solved by establishing a convention. We investigate start-up and free rider subproblems as well as the effects of group size, intensity of intrinsic preference, and salience on the emergence dynamics of coordination conventions. Results of our simulations show agents establish and switch between conventions, even working against their own preferred outcome when doing so is necessary for effective coordination.


title: Large-Scale Maintenance and Unit Commitment: A Decentralized Subgradient Approach

abstract: Unit Commitment (UC) is a fundamental problem in power system operations. When coupled with generation maintenance, the joint optimization problem poses significant computational challenges due to coupling constraints linking maintenance and UC decisions. Obviously, these challenges grow with the size of the network. With the introduction of sensors for monitoring generator health and condition-based maintenance(CBM), these challenges have been magnified. ADMM-based decentralized methods have shown promise in solving large-scale UC problems, especially in vertically integrated power systems. However, in their current form, these methods fail to deliver similar computational performance and scalability when considering the joint UC and CBM problem. This paper provides a novel decentralized optimization framework for solving large-scale, joint UC and CBM problems. Our approach relies on the novel use of the subgradient method to temporally decouple various subproblems of the ADMM-based formulation of the joint problem along the maintenance horizon. By effectively utilizing multithreading, our decentralized subgradient approach delivers superior computational performance and eliminates the need to move sensor data thereby alleviating privacy and security concerns. Using experiments on large scale test cases, we show that our framework can provide a speedup of upto 50x as compared to various state of the art benchmarks without compromising on solution quality.


title: Social-VRNN: One-Shot Multi-modal Trajectory Prediction for Interacting Pedestrians

abstract: Prediction of human motions is key for safe navigation of autonomous robots among humans. In cluttered environments, several motion hypotheses may exist for a pedestrian, due to its interactions with the environment and other pedestrians. Previous works for estimating multiple motion hypotheses require a large number of samples which limits their applicability in real-time motion planning. In this paper, we present a variational learning approach for interaction-aware and multi-modal trajectory prediction based on deep generative neural networks. Our approach can achieve faster convergence and requires significantly fewer samples comparing to state-of-the-art methods. Experimental results on real and simulation data show that our model can effectively learn to infer different trajectories. We compare our method with three baseline approaches and present performance results demonstrating that our generative model can achieve higher accuracy for trajectory prediction by producing diverse trajectories.


title: Asynchronous Optimization over Graphs: Linear Convergence under Error Bound Conditions

abstract: We consider convex and nonconvex constrained optimization with a partially separable objective function: agents minimize the sum of local objective functions, each of which is known only by the associated agent and depends on the variables of that agent and those of a few others. This partitioned setting arises in several applications of practical interest. We propose what is, to the best of our knowledge, the first distributed, asynchronous algorithm with rate guarantees for this class of problems. When the objective function is nonconvex, the algorithm provably converges to a stationary solution at a sublinear rate whereas linear rate is achieved when the objective satisfies under the renowned Luo-Tseng error bound condition (which is less stringent than strong convexity). Numerical results on matrix completion and LASSO problems show the effectiveness of our method.


title: Coregular submanifolds and Poisson submersions

abstract: We analyze \emph{submersions with Poisson fibres}. These are submersions whose total space carries a Poisson structure, on which the ambient Poisson structure pulls back, as a Dirac structure, to Poisson structures on each individual fibre. Our ``Poisson-Dirac viewpoint'' is prompted by natural examples of Poisson submersions with Poisson fibers -- in toric geometry and Poisson-Lie groups -- whose analysis was not possible using the existing tools in the Poisson literature. The first part of the paper studies the Poisson-Dirac perspective of inducing Poisson structures on submanifolds. This is a rich landscape, in which subtle behaviours abound -- as illustrated by a surprising ``jumping phenomenon'' concerning the complex relation between the induced and the ambient symplectic foliations, which we discovered here. These pathologies, however, are absent from the well-behaved and abundant class of \emph{coregular} submanifolds, with which we are mostly concerned here. The second part of the paper studies Poisson submersions with Poisson fibres -- the natural Poisson generalization of flat symplectic bundles. These Poisson submersions have coregular Poisson-Dirac fibres, and behave functorially with respect to such submanifolds. We discuss the subtle collective behavior of the Poisson fibers of such Poisson fibrations, and explain their relation to pencils of Poisson structures. The third and final part applies the theory developed to Poisson submersions with Poisson fibres which arise in Lie theory. We also show that such submersions are a convenient setting for the associated bundle construction, and we illustrate this by producing new Poisson structures with a finite number of symplectic leaves. Some of the points in the paper being fairly new, we illustrate the many fine issues that appear with an abundance of (counter-)examples.


title: Discrete Empirical Interpolation and unfitted mesh FEMs: application in PDE-constrained optimization

abstract: In this work, we investigate the performance CutFEM as a high fidelity solver as well as we construct a competent and economical reduced order solver for PDE-constrained optimization problems in parametrized domains that live in a fixed background geometry and mesh. Its effectiveness and reliability will be assessed through its application for the numerical solution of quadratic optimization problems with elliptic equations as constraints, examining an archetypal case. The reduction strategy will be via Proper Orthogonal Decomposition of suitable FE snapshots, using an aggregated state and adjoint test space, while the efficiency of the offline-online decoupling will be ensured by means of Discrete Empirical Interpolation of the optimality system matrix and right-hand side, enabling thus a rapid resolution of the reduced order model for each new spatial configuration.


title: Transitional millisecond pulsars

abstract: Millisecond pulsars in tight binaries have recently opened new challenges in our understanding of physical processes governing the evolution of binaries and the interaction between astrophysical plasma and electromagnetic fields. Transitional systems that showed changes from rotation-powered to accretion powered states and vice versa have bridged the populations of radio and accreting millisecond pulsars, eventually demonstrating the tight evolutionary link envisaged by the recycling scenario. A decade of discoveries and theoretical efforts have just grasped the complex phenomenology of transitional millisecond pulsars from the radio to the gamma-ray band. This review summarises the main properties of the three transitional millisecond pulsars discovered so far, as well as of candidates and related systems, discussing the various models proposed to cope with the multifaceted behaviour.


title: The Third Data Release of the KODIAQ Survey

abstract: We present and make publicly available the third data release (DR3) of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ) survey. KODIAQ DR3 consists of a fully-reduced sample of 727 quasars at 0.1 < z < 6.4 observed with ESI at moderate resolution (4000 < R < 10000). DR3 contains 872 spectra available in flux calibrated form, representing a sum total exposure time of approximately 2.8 megaseconds. These co-added spectra arise from a total of 2753 individual exposures of quasars taken from the Keck Observatory Archive (KOA) in raw form and uniformly processed using a data reduction package made available through the XIDL distribution. DR3 is publicly available to the community, housed as a higher level science product at the KOA and in the igmspec database.


title: The Nondipole Part of the Main Geomagnetic Field and the Large Scale Topographical Heterogeneities of the Core-Mantle Boundary

abstract: Based on satellite data, both large and small-scale flows were identified on the core surface .At the same time, the observed rapid decrease in the dipole component of the main geomagnetic field does not find a satisfactory explanation. In the previous papers we carried out the comparative studies of the structure of the lowest mantle with the motion path of the small size current loops approximating the MGF's small-scale anomalies. The hypothesis was stated that heterogeneities in the lowest mantle structure and the topographic irregularities of the core-mantle boundary associated with ancient subduction zones is one of the primarily responsible for the small-scale vortices formation and the local MGF variation. In this paper, we carry out the analogous comparison for the large scale sources approximating the nondipole part of the geomagnetic field.


title: Enabling Fast Differentially Private SGD via Just-in-Time Compilation and Vectorization

abstract: A common pain point in differentially private machine learning is the significant runtime overhead incurred when executing Differentially Private Stochastic Gradient Descent (DPSGD), which may be as large as two orders of magnitude. We thoroughly demonstrate that by exploiting powerful language primitives, including vectorization, just-in-time compilation, and static graph optimization, one can dramatically reduce these overheads, in many cases nearly matching the best non-private running times. These gains are realized in two frameworks: JAX and TensorFlow. JAX provides rich support for these primitives as core features of the language through the XLA compiler. We also rebuild core parts of TensorFlow Privacy, integrating features from TensorFlow 2 as well as XLA compilation, granting significant memory and runtime improvements over the current release version. These approaches allow us to achieve up to 50x speedups in comparison to the best alternatives. Our code is available at https://github.com/TheSalon/fast-dpsgd.


title: Massive particles in the Einstein-Lovelock-anti-de Sitter black hole spacetime

abstract: An interpretation to the physics of stable geodesics of massive particles in black hole backgrounds has been recently proposed in the context of the AdS/CFT correspondence. It was argued that the existence of stable orbits indicates that the dual state does not thermalize on a thermal time scale and the bulk excitations can be interpreted as metastable states in the dual field theory. Here we study motion of massive particles in the background of the D-dimensional asymptotically anti-de Sitter (AdS) black holes in the Einstein-Lovelock theory. We show that, unlike the asymptotically flat case, for any kind of higher curvature Lovelock corrections there appear a stable circular orbit at a distance from the black hole. We find the general analytical expressions for the frequencies of distant circular orbits and radial oscillation frequencies. We show that the corresponding correction has the same power as in the Schwarzschild-AdS case, implying a universal scaling with the temperature for any Lovelock theory.


title: Long-time behavior of scalar conservation laws with critical dissipation

abstract: The critical Burgers equation $\partial_t u + u \partial_x u + \Lambda u = 0$ is a toy model for the competition between transport and diffusion with regard to shock formation in fluids. It is well known that smooth initial data does not generate shocks in finite time. Less is known about the long-time behavior for `shock-like' initial data: $u_0 \to \pm a$ as $x \to \mp \infty$. We describe this long-time behavior in the general setting of multidimensional critical scalar conservation laws $\partial_t u + \text{div}f(u) + \Lambda u = 0$ when the initial data has limits at infinity. The asymptotics are given by certain self-similar solutions, whose stability we demonstrate with the optimal diffusive rates.


title: Exploiting Context for Robustness to Label Noise in Active Learning

abstract: Several works in computer vision have demonstrated the effectiveness of active learning for adapting the recognition model when new unlabeled data becomes available. Most of these works consider that labels obtained from the annotator are correct. However, in a practical scenario, as the quality of the labels depends on the annotator, some of the labels might be wrong, which results in degraded recognition performance. In this paper, we address the problems of i) how a system can identify which of the queried labels are wrong and ii) how a multi-class active learning system can be adapted to minimize the negative impact of label noise. Towards solving the problems, we propose a noisy label filtering based learning approach where the inter-relationship (context) that is quite common in natural data is utilized to detect the wrong labels. We construct a graphical representation of the unlabeled data to encode these relationships and obtain new beliefs on the graph when noisy labels are available. Comparing the new beliefs with the prior relational information, we generate a dissimilarity score to detect the incorrect labels and update the recognition model with correct labels which result in better recognition performance. This is demonstrated in three different applications: scene classification, activity classification, and document classification.


title: Multimodal semantic forecasting based on conditional generation of future features

abstract: This paper considers semantic forecasting in road-driving scenes. Most existing approaches address this problem as deterministic regression of future features or future predictions given observed frames. However, such approaches ignore the fact that future can not always be guessed with certainty. For example, when a car is about to turn around a corner, the road which is currently occluded by buildings may turn out to be either free to drive, or occupied by people, other vehicles or roadworks. When a deterministic model confronts such situation, its best guess is to forecast the most likely outcome. However, this is not acceptable since it defeats the purpose of forecasting to improve security. It also throws away valuable training data, since a deterministic model is unable to learn any deviation from the norm. We address this problem by providing more freedom to the model through allowing it to forecast different futures. We propose to formulate multimodal forecasting as sampling of a multimodal generative model conditioned on the observed frames. Experiments on the Cityscapes dataset reveal that our multimodal model outperforms its deterministic counterpart in short-term forecasting while performing slightly worse in the mid-term case.


title: Differentiation of subjects of the Russian Federation according to the main parameters of socio-economic development

abstract: This article presents the results of a cluster analysis of the regions of the Russian Federation in terms of the main parameters of socio-economic development according to the data presented in the official data sources of the Federal State Statistics Service (Rosstat). Studied and analyzed the domestic and foreign (Eurostat) methodology for assessing the socio-economic development of territories. The aim of the study is to determine the main parameters of territorial differentiation and to identify key indicators that affect the socio-economic development of Russian regions. The authors have carried out a classification of the constituent entities of the Russian Federation not in terms of territorial location and geographical features, but in terms of the specifics and key parameters of the socio-economic situation.


title: Littlewood and Duffin--Schaeffer-type problems in diophantine approximation

abstract: Gallagher's theorem describes the multiplicative diophantine approximation rate of a typical vector. We establish a fully-inhomogeneous version of Gallagher's theorem, a diophantine fibre refinement, and a sharp and unexpected threshold for Liouville fibres. Along the way, we prove an inhomogeneous version of the Duffin--Schaeffer conjecture for a class of non-monotonic approximation functions.


title: Catalytic transformations with finite-size environments: applications to cooling and thermometry

abstract: The laws of thermodynamics are usually formulated under the assumption of infinitely large environments. While this idealization facilitates theoretical treatments, real physical systems are always finite and their interaction range is limited. These constraints have consequences for important tasks such as cooling, not directly captured by the second law of thermodynamics. Here, we study catalytic transformations that cannot be achieved when a system exclusively interacts with a finite environment. Our core result consists of constructive conditions for these transformations, which include the corresponding global unitary operation and the explicit states of all the systems involved. From this result we present various findings regarding the use of catalysts for cooling. First, we show that catalytic cooling is always possible if the dimension of the catalyst is sufficiently large. In particular, the cooling of a qubit using a hot qubit can be maximized with a catalyst as small as a three-level system. We also identify catalytic enhancements for tasks whose implementation is possible without a catalyst. For example, we find that in a multiqubit setup catalytic cooling based on a three-body interaction outperforms standard (non-catalytic) cooling using higher order interactions. Another advantage is illustrated in a thermometry scenario, where a qubit is employed to probe the temperature of the environment. In this case, we show that a catalyst allows to surpass the optimal temperature estimation attained only with the probe.


title: Oscillations for order statistics of some discrete processes

abstract: Suppose $k$ balls are dropped into $n$ boxes independently with uniform probability, where $n, k$ are large with ratio approximately equal to some positive real $\lambda$. The maximum box count has a counterintuitive behavior: first of all, with high probability it takes at most two values $m_n$ or $m_n+1$, where $m_n$ is roughly $\frac{\ln n}{\ln \ln n}$. Moreover, it oscillates between these two values with an unusual periodicity. In order to prove this statement and various generalizations, it is first shown that for $X_1,...,X_n$ independent and identically distributed discrete random variables with common distribution $F$, under mild conditions, the limiting distribution of their maximum oscillates in three possible families, depending on the tail of the distribution. The result stated at the beginning follows from the equivalence of ensemble for the order statistics in various allocations problems, obtained via conditioning limit theory. Results about the number of ties for the maximum, as well as applications, are also provided.


title: UoB at SemEval-2020 Task 1: Automatic Identification of Novel Word Senses

abstract: Much as the social landscape in which languages are spoken shifts, language too evolves to suit the needs of its users. Lexical semantic change analysis is a burgeoning field of semantic analysis which aims to trace changes in the meanings of words over time. This paper presents an approach to lexical semantic change detection based on Bayesian word sense induction suitable for novel word sense identification. This approach is used for a submission to SemEval-2020 Task 1, which shows the approach to be capable of the SemEval task. The same approach is also applied to a corpus gleaned from 15 years of Twitter data, the results of which are then used to identify words which may be instances of slang.


title: Handling Bidirectional Control Flow: Technical Report

abstract: Pressed by the difficulty of writing asynchronous, event-driven code, mainstream languages have recently been building in support for a variety of advanced control-flow features. Meanwhile, experimental language designs have suggested effect handlers as a unifying solution to programmer-defined control effects, subsuming exceptions, generators, and async--await. Despite these trends, complex control flow---in particular, control flow exhibiting a bidirectional pattern---remains challenging to manage. We introduce bidirectional algebraic effects, a new programming abstraction that supports bidirectional control transfer in a more natural way. Handlers of bidirectional effects can raise further effects to transfer control back to the site where the initiating effect was raised, and can use themselves to handle their own effects. We present applications of this expressive power, which falls out naturally as we push toward the unification of effectful programming with object-oriented programming. We pin down the mechanism and the unification formally using a core language that generalizes to effect operations and effect handlers. The usual propagation semantics of control effects such as exceptions conflicts with modular reasoning in the presence of effect polymorphism---it breaks parametricity. Bidirectionality exacerbates the problem. Hence, we set out to show the core language, which builds on the existing tunneling semantics for algebraic effects, is not only type-safe (no effects go unhandled), but also abstraction-safe (no effects are accidentally handled). We devise a step-indexed logical-relations model, and construct its parametricity and soundness proofs. These core results are fully mechanized in Coq. Preliminary experiments show that as a first-class language feature, bidirectional handlers can be implemented efficiently.


title: The Duopoly Analysis of Graphics Card Market

abstract: By analyzing the duopoly market of computer graphics cards, we categorized the effects of enterprise's technological progress into two types, namely, cost reduction and product diversification. Our model proved that technological progress is the most effective means for enterprises in this industry to increase profits. Due to the technology-intensive nature of this industry, monopolistic enterprises face more intense competition compared with traditional manufacturing. Therefore, they have more motivation for technological innovation. Enterprises aiming at maximizing profits have incentives to reduce costs and achieve a higher degree of product differentiation through technological innovation.


title: Bayesian Quantum Multiphase Estimation Algorithm

abstract: Quantum phase estimation (QPE) is the key subroutine of several quantum computing algorithms as well as a central ingredient in quantum computational chemistry and quantum simulation. While QPE strategies have focused on the estimation of a single phase, applications to the simultaneous estimation of several phases may bring substantial advantages; for instance, in the presence of spatial or temporal constraints. In this work, we study a Bayesian algorithm for the parallel (simultaneous) estimation of multiple arbitrary phases. The protocol gives access to correlations in the Bayesian multi-phase distribution resulting in covariance matrix elements scaling inversely proportional to the square of the total number of quantum resources. The parallel estimation allows to surpass the sensitivity of sequential single-phase estimation strategies for optimal linear combinations of phases. Furthermore, the algorithm proves a certain noise resilience and can be implemented using single photons and standard optical elements in currently accessible experiments.


title: RADIATE: A Radar Dataset for Automotive Perception in Bad Weather

abstract: Datasets for autonomous cars are essential for the development and benchmarking of perception systems. However, most existing datasets are captured with camera and LiDAR sensors in good weather conditions. In this paper, we present the RAdar Dataset In Adverse weaThEr (RADIATE), aiming to facilitate research on object detection, tracking and scene understanding using radar sensing for safe autonomous driving. RADIATE includes 3 hours of annotated radar images with more than 200K labelled road actors in total, on average about 4.6 instances per radar image. It covers 8 different categories of actors in a variety of weather conditions (e.g., sun, night, rain, fog and snow) and driving scenarios (e.g., parked, urban, motorway and suburban), representing different levels of challenge. To the best of our knowledge, this is the first public radar dataset which provides high-resolution radar images on public roads with a large amount of road actors labelled. The data collected in adverse weather, e.g., fog and snowfall, is unique. Some baseline results of radar based object detection and recognition are given to show that the use of radar data is promising for automotive applications in bad weather, where vision and LiDAR can fail. RADIATE also has stereo images, 32-channel LiDAR and GPS data, directed at other applications such as sensor fusion, localisation and mapping. The public dataset can be accessed at http://pro.hw.ac.uk/radiate/.


title: A Spatial-Temporal Graph Based Hybrid Infectious Disease Model with Application to COVID-19

abstract: As the COVID-19 pandemic evolves, reliable prediction plays an important role for policy making. The classical infectious disease model SEIR (susceptible-exposed-infectious-recovered) is a compact yet simplistic temporal model. The data-driven machine learning models such as RNN (recurrent neural networks) can suffer in case of limited time series data such as COVID-19. In this paper, we combine SEIR and RNN on a graph structure to develop a hybrid spatio-temporal model to achieve both accuracy and efficiency in training and forecasting. We introduce two features on the graph structure: node feature (local temporal infection trend) and edge feature (geographic neighbor effect). For node feature, we derive a discrete recursion (called I-equation) from SEIR so that gradient descend method applies readily to its optimization. For edge feature, we design an RNN model to capture the neighboring effect and regularize the landscape of loss function so that local minima are effective and robust for prediction. The resulting hybrid model (called IeRNN) improves the prediction accuracy on state-level COVID-19 new case data from the US, out-performing standard temporal models (RNN, SEIR, and ARIMA) in 1-day and 7-day ahead forecasting. Our model accommodates various degrees of reopening and provides potential outcomes for policymakers.


title: Incorporating Count-Based Features into Pre-Trained Models for Improved Stance Detection

abstract: The explosive growth and popularity of Social Media has revolutionised the way we communicate and collaborate. Unfortunately, this same ease of accessing and sharing information has led to an explosion of misinformation and propaganda. Given that stance detection can significantly aid in veracity prediction, this work focuses on boosting automated stance detection, a task on which pre-trained models have been extremely successful on, as on several other tasks. This work shows that the task of stance detection can benefit from feature based information, especially on certain under performing classes, however, integrating such features into pre-trained models using ensembling is challenging. We propose a novel architecture for integrating features with pre-trained models that address these challenges and test our method on the RumourEval 2019 dataset. This method achieves state-of-the-art results with an F1-score of 63.94 on the test set.


title: Graphite: GRAPH-Induced feaTure Extraction for Point Cloud Registration

abstract: 3D Point clouds are a rich source of information that enjoy growing popularity in the vision community. However, due to the sparsity of their representation, learning models based on large point clouds is still a challenge. In this work, we introduce Graphite, a GRAPH-Induced feaTure Extraction pipeline, a simple yet powerful feature transform and keypoint detector. Graphite enables intensive down-sampling of point clouds with keypoint detection accompanied by a descriptor. We construct a generic graph-based learning scheme to describe point cloud regions and extract salient points. To this end, we take advantage of 6D pose information and metric learning to learn robust descriptions and keypoints across different scans. We Reformulate the 3D keypoint pipeline with graph neural networks which allow efficient processing of the point set while boosting its descriptive power which ultimately results in more accurate 3D registrations. We demonstrate our lightweight descriptor on common 3D descriptor matching and point cloud registration benchmarks and achieve comparable results with the state of the art. Describing 100 patches of a point cloud and detecting their keypoints takes only ~0.018 seconds with our proposed network.


title: Poisoned classifiers are not only backdoored, they are fundamentally broken

abstract: Under a commonly-studied backdoor poisoning attack against classification models, an attacker adds a small trigger to a subset of the training data, such that the presence of this trigger at test time causes the classifier to always predict some target class. It is often implicitly assumed that the poisoned classifier is vulnerable exclusively to the adversary who possesses the trigger. In this paper, we show empirically that this view of backdoored classifiers is incorrect. We describe a new threat model for poisoned classifier, where one without knowledge of the original trigger, would want to control the poisoned classifier. Under this threat model, we propose a test-time, human-in-the-loop attack method to generate multiple effective alternative triggers without access to the initial backdoor and the training data. We construct these alternative triggers by first generating adversarial examples for a smoothed version of the classifier, created with a procedure called Denoised Smoothing, and then extracting colors or cropped portions of smoothed adversarial images with human interaction. We demonstrate the effectiveness of our attack through extensive experiments on high-resolution datasets: ImageNet and TrojAI. We also compare our approach to previous work on modeling trigger distributions and find that our method are more scalable and efficient in generating effective triggers. Last, we include a user study which demonstrates that our method allows users to easily determine the existence of such backdoors in existing poisoned classifiers. Thus, we argue that there is no such thing as a secret backdoor in poisoned classifiers: poisoning a classifier invites attacks not just by the party that possesses the trigger, but from anyone with access to the classifier.


title: A Ramsey Characterisation of Eventually Periodic Words

abstract: A factorisation $x = u_1 u_2 \cdots$ of an infinite word $x$ on alphabet $X$ is called `monochromatic', for a given colouring of the finite words $X^*$ on alphabet $X$, if each $u_i$ is the same colour. Wojcik and Zamboni proved that the word $x$ is periodic if and only if for every finite colouring of $X^*$ there is a monochromatic factorisation of $x$. On the other hand, it follows from Ramsey's theorem that, for \textit{any} word $x$, for every finite colouring of $X^*$ there is a suffix of $x$ having a monochromatic factorisation. A factorisation $x = u_1 u_2 \cdots$ is called `super-monochromatic' if each word $u_{k_1} u_{k_2} \cdots u_{k_n}$, where $k_1 < \cdots < k_n$, is the same colour. Our aim in this paper is to show that a word $x$ is eventually periodic if and only if for every finite colouring of $X^*$ there is a suffix of $x$ having a super-monochromatic factorisation. Our main tool is a Ramsey result about alternating sums that may be of independent interest.


title: Fully-heavy tetraquark spectra and production at hadron colliders

abstract: Motivated by the observation of exotic structure around 6900MeV in the $J/\psi$-pair mass spectrum using proton-proton collision data by the LHCb collaboration, we study the spectra of fully-heavy tetraquarks within Bethe-Salpeter equation and Regge trajectory relation. The $X(6900)$ may be explained as a radially excited state with quark content $cc\bar{c}\bar{c}$ and spin-parity $0^{++}(3S)$ or $2^{++}(3S)$ or an orbitally excited $2P$ state. New $cc\bar{c}\bar{c}$ structures around 6.0GeV, 6.5GeV, and 7.1GeV are predicted together. Other $bb\bar{b}\bar{b}$ and $bc\bar{b}\bar{c}$ structures which may be experimentally prominent are discussed. On the other hand, the fully-heavy S-wave tetraquark productions at hadron colliders are investigated and their cross-sections are obtained.


title: Large fluctuations and transport properties of the L\'evy-Lorentz gas

abstract: The L\'evy-Lorentz gas describes the motion of a particle on the real line in the presence of a random array of scattering points, whose distances between neighboring points are heavy-tailed i.i.d. random variables with finite mean. The motion is a continuous-time, constant-speed interpolation of the simple symmetric random walk on the marked points. In this paper we study the large fluctuations of the continuous-time process and the resulting transport properties of the model, both annealed and quenched, confirming and extending previous work by physicists that pertain to the annealed framework. Specifically, focusing on the particle displacement, and under the assumption that the tail distribution of the interdistances between scatterers is regularly varying at infinity, we prove a precise large deviation principle for the annealed fluctuations and present the asymptotics of annealed moments, demonstrating annealed superdiffusion. Then, we provide an upper large deviation estimate for the quenched fluctuations and the asymptotics of quenched moments, showing that the asymptotic diffusive regime conditional on a typical arrangement of the scatterers is normal diffusion, and not superdiffusion. Although the L\'evy-Lorentz gas seems to be accepted as a model for anomalous diffusion, our findings suggest that superdiffusion is a transient behavior which develops into normal diffusion on long timescales, and raise a new question about how the transition from the quenched normal diffusion to the annealed superdiffusion occurs.

categories: math.PR cond-mat.stat-mech math-ph math.MP


title: Gait Recognition using Multi-Scale Partial Representation Transformation with Capsules

abstract: Gait recognition, referring to the identification of individuals based on the manner in which they walk, can be very challenging due to the variations in the viewpoint of the camera and the appearance of individuals. Current methods for gait recognition have been dominated by deep learning models, notably those based on partial feature representations. In this context, we propose a novel deep network, learning to transfer multi-scale partial gait representations using capsules to obtain more discriminative gait features. Our network first obtains multi-scale partial representations using a state-of-the-art deep partial feature extractor. It then recurrently learns the correlations and co-occurrences of the patterns among the partial features in forward and backward directions using Bi-directional Gated Recurrent Units (BGRU). Finally, a capsule network is adopted to learn deeper part-whole relationships and assigns more weights to the more relevant features while ignoring the spurious dimensions. That way, we obtain final features that are more robust to both viewing and appearance changes. The performance of our method has been extensively tested on two gait recognition datasets, CASIA-B and OU-MVLP, using four challenging test protocols. The results of our method have been compared to the state-of-the-art gait recognition solutions, showing the superiority of our model, notably when facing challenging viewing and carrying conditions.


title: Logistic wavelets and logistic function: An application to model the spread of SARS-CoV-2 virus infections

abstract: In the present paper, we model the cumulative number of persons reported to be infected by the SARS-CoV-2 virus, in a country or a region, by a sum of logistic functions. For a given logistic function, using Eulerian numbers, we find the zeros of its successive derivatives and their relationship with the saturation level of this function. In a given time series, having potentially the logistic trend, we use its second differences to determine points corresponding to these zeros. To estimate the parameters of the approximating logistic function, we define and use logistic wavelets. Then we apply the theory to the cases of SARS-CoV-2 infections in the United States and the United Kingdom.


title: Blockchain Based Decentralized Replay Attack Detection for Large Scale Power Systems

abstract: Large scale power systems are comprised of regional utilities with assets that stream sensor readings in real time. In order to detect cyberattacks, the globally acquired, real time sensor data needs to be analyzed in a centralized fashion. However, owing to operational constraints, such a centralized sharing mechanism turns out to be a major obstacle. In this paper, we propose a blockchain based decentralized framework for detecting coordinated replay attacks with full privacy of sensor data. We develop a Bayesian inference mechanism employing locally reported attack probabilities that is tailor made for a blockchain framework. We compare our framework to a traditional decentralized algorithm based on the broadcast gossip framework both theoretically as well as empirically. With the help of experiments on a private Ethereum blockchain, we show that our approach achieves good detection quality and significantly outperforms gossip driven approaches in terms of accuracy, timeliness and scalability.


title: Wireless Control for Smart Manufacturing: Recent Approaches and Open Challenges

abstract: Smart manufacturing aims to overcome the limitations of today's rigid assembly lines by making the material flow and manufacturing process more flexible, versatile, and scalable. The main economic drivers are higher resource and cost efficiency as the manufacturers can more quickly adapt to changing market needs and also increase the lifespan of their production sites. The ability to close feedback loops fast and reliably over long distances among mobile robots, remote sensors, and human operators is a key enabler for smart manufacturing. Thus, this article provides a perspective on control and coordination over wireless networks. Based on an analysis of real-world use cases, we identify the main technical challenges that need to be solved to close the large gap between the current state of the art in industry and the vision of smart manufacturing. We discuss to what extent existing control-over-wireless solutions in the literature address those challenges, including our own approach toward a tight integration of control and wireless communication. In addition to a theoretical analysis of closed-loop stability, practical experiments on a cyber-physical testbed demonstrate that our approach supports relevant smart manufacturing scenarios. The article concludes with a discussion of open challenges and future research directions.


title: Energy-based error bound of physics-informed neural network solutions in elasticity

abstract: An energy-based a posteriori error bound is proposed for the physics-informed neural network solutions of elasticity problems. An admissible displacement-stress solution pair is obtained from a mixed form of physics-informed neural networks, and the proposed error bound is formulated as the constitutive relation error defined by the solution pair. Such an error estimator provides an upper bound of the global error of neural network discretization. The bounding property, as well as the asymptotic behavior of the physics-informed neural network solutions, are studied in a demonstrating example.


title: Training Stronger Baselines for Learning to Optimize

abstract: Learning to optimize (L2O) has gained increasing attention since classical optimizers require laborious problem-specific design and hyperparameter tuning. However, there is a gap between the practical demand and the achievable performance of existing L2O models. Specifically, those learned optimizers are applicable to only a limited class of problems, and often exhibit instability. With many efforts devoted to designing more sophisticated L2O models, we argue for another orthogonal, under-explored theme: the training techniques for those L2O models. We show that even the simplest L2O model could have been trained much better. We first present a progressive training scheme to gradually increase the optimizer unroll length, to mitigate a well-known L2O dilemma of truncation bias (shorter unrolling) versus gradient explosion (longer unrolling). We further leverage off-policy imitation learning to guide the L2O learning, by taking reference to the behavior of analytical optimizers. Our improved training techniques are plugged into a variety of state-of-the-art L2O models, and immediately boost their performance, without making any change to their model structures. Especially, by our proposed techniques, an earliest and simplest L2O model can be trained to outperform the latest complicated L2O models on a number of tasks. Our results demonstrate a greater potential of L2O yet to be unleashed, and urge to rethink the recent progress. Our codes are publicly available at: https://github.com/VITA-Group/L2O-Training-Techniques.


title: A Framework for Multiphase Galactic Wind Launching using TIGRESS

abstract: Galactic outflows have density, temperature, and velocity variations at least as large as that of the multiphase, turbulent interstellar medium (ISM) from which they originate. We have conducted a suite of parsec-resolution numerical simulations using the TIGRESS framework, in which outflows emerge as a consequence of interaction between supernovae (SNe) and the star-forming ISM. The outflowing gas is characterized by two distinct thermal phases, cool (T<10^4 K) and hot (T>10^6 K), with most mass carried by the cool phase and most energy and newly-injected metals carried by the hot phase. Both components have a broad distribution of outflow velocity, and especially for cool gas this implies a varying fraction of escaping material depending on the halo potential. Informed by the TIGRESS results, we develop straightforward analytic formulae for the joint probability density functions (PDFs) of mass, momentum, energy, and metal loading as distributions in outflow velocity and sound speed. The model PDFs have only two parameters, SFR surface density \Sigma_SFR and the metallicity of the ISM, and fully capture the behavior of the original TIGRESS simulation PDFs over \Sigma_SFR~(10^{-4},1)M_sun/kpc^2/yr. Employing PDFs from resolved simulations will enable galaxy formation subgrid model implementations with wind velocity and temperature (as well as total loading factors) that are based on theoretical predictions rather than empirical tuning. This is a critical step to incorporate advances from TIGRESS and other high-resolution simulations in future cosmological hydrodynamics and semi-analytic galaxy formation models. We release a python package to prototype our model and to ease its implementation.


title: Colourings of $(m, n)$-coloured mixed graphs

abstract: A mixed graph is, informally, an object obtained from a simple undirected graph by choosing an orientation for a subset of its edges. A mixed graph is $(m, n)$-coloured if each edge is assigned one of $m \geq 0$ colours, and each arc is assigned one of $n \geq 0$ colours. Oriented graphs are $(0, 1)$-coloured mixed graphs, and 2-edge-coloured graphs are $(2, 0)$-coloured mixed graphs. We show that results of Sopena for vertex colourings of oriented graphs, and of Kostochka, Sopena and Zhu for vertex colourings oriented graphs and 2-edge-coloured graphs, are special cases of results about vertex colourings of $(m, n)$-coloured mixed graphs. Both of these can be regarded as a version of Brooks' Theorem.


title: View-Invariant Gait Recognition with Attentive Recurrent Learning of Partial Representations

abstract: Gait recognition refers to the identification of individuals based on features acquired from their body movement during walking. Despite the recent advances in gait recognition with deep learning, variations in data acquisition and appearance, namely camera angles, subject pose, occlusions, and clothing, are challenging factors that need to be considered for achieving accurate gait recognition systems. In this paper, we propose a network that first learns to extract gait convolutional energy maps (GCEM) from frame-level convolutional features. It then adopts a bidirectional recurrent neural network to learn from split bins of the GCEM, thus exploiting the relations between learned partial spatiotemporal representations. We then use an attention mechanism to selectively focus on important recurrently learned partial representations as identity information in different scenarios may lie in different GCEM bins. Our proposed model has been extensively tested on two large-scale CASIA-B and OU-MVLP gait datasets using four different test protocols and has been compared to a number of state-of-the-art and baseline solutions. Additionally, a comprehensive experiment has been performed to study the robustness of our model in the presence of six different synthesized occlusions. The experimental results show the superiority of our proposed method, outperforming the state-of-the-art, especially in scenarios where different clothing and carrying conditions are encountered. The results also revealed that our model is more robust against different occlusions as compared to the state-of-the-art methods.


title: On the irreducibility of the extensions of Burau and Gassner representations

abstract: We study the $n$th degree representations $\hat{\rho_G}$ of $Cb_{n}$ and $\hat{\rho_B}$ of $C_{n}$, defined by Valerij G. Bardakov, where $Cb_{n}$ is the group of basis conjugating automorphisms and $C_n$ is the group of conjugating automorphisms. We prove that $\hat{\rho_G}$ is reducible and its $(n-1)$th degree composition factor $\hat{\phi_G}$ is irreducible if and only if $t_i\neq 1$ for all $1\leq i \leq n$. Also we prove that $\hat{\rho_B}$ is reducible and its $(n-1)$th degree composition factor $\hat{\phi_B}$ is irreducible if and only if $t\neq 1$. Moreover, for $n=3$, we prove that $\hat{\phi_G}(t_1,t_2,t_3) \otimes \hat{\phi_G}(m_1,m_2,m_3)$ is irreducible if and only if $(t_1,t_2,t_3)$ and $(m_1,m_2,m_3)$ are distinct vectors, and the representation $\hat{\phi_B}(t) \otimes \hat{\phi_B}(m)$ is irreducible if and only if $t \neq m$.


title: Multi-Agent Reinforcement Learning in NOMA-aided UAV Networks for Cellular Offloading

abstract: A novel framework is proposed for cellular offloading with the aid of multiple unmanned aerial vehicles (UAVs), while the non-orthogonal multiple access (NOMA) technique is employed at each UAV to further improve the spectrum efficiency of the wireless network. The optimization problem of joint three-dimensional (3D) trajectory design and power allocation is formulated for maximizing the throughput. Since ground mobile users are considered as roaming continuously, the UAVs need to be re-deployed timely based on the movement of users. In an effort to solve this pertinent dynamic problem, a K-means based clustering algorithm is first adopted for periodically partitioning users. Afterward, a mutual deep Q-network (MDQN) algorithm is proposed to jointly determine the optimal 3D trajectory and power allocation of UAVs. In contrast to the conventional DQN algorithm, the MDQN algorithm enables the experience of multi-agent to be input into a shared neural network to shorten the training time with the assistance of state abstraction. Numerical results demonstrate that: 1) the proposed MDQN algorithm is capable of converging under minor constraints and has a faster convergence rate than the conventional DQN algorithm in the multi-agent case; 2) The achievable sum rate of the NOMA enhanced UAV network is 23% superior to the case of orthogonal multiple access (OMA); 3) By designing the optimal 3D trajectory of UAVs with the aid of the MDON algorithm, the sum rate of the network enjoys 142% and 56% gains than that of invoking the circular trajectory and the 2D trajectory, respectively.


title: First-principles computation of boron-nitride-based ultrathin UV-C light emitting diodes

abstract: Short wavelength ultraviolet (UV-C) light deactivates DNA of any germs, including multiresistive bacteria and viruses like COVID-19. Two-dimensional (2D) material-based UV-C light emitting diodes can potentially be integrated into arbitrary surfaces to allow for shadow-free surface disinfection. In this work, we perform a series of first-principles calculations to identify the core components of ultrathin LEDs based on hexagonal boron nitride (hBN). The electrons and holes are predicted to be confined in multiple quantum wells (MQWs) by combining hBN layers with different stacking orders. Various p- and n-doping candidates for hBN are assessed, and the relative p- and n-type metal contacts with low Schottky barrier heights are identified. The findings are summarized in a concrete UV-C LED structure proposal.


title: On Near-Linear-Time Algorithms for Dense Subset Sum

abstract: In the Subset Sum problem we are given a set of $n$ positive integers $X$ and a target $t$ and are asked whether some subset of $X$ sums to $t$. Natural parameters for this problem that have been studied in the literature are $n$ and $t$ as well as the maximum input number $\rm{mx}_X$ and the sum of all input numbers $\Sigma_X$. In this paper we study the dense case of Subset Sum, where all these parameters are polynomial in $n$. In this regime, standard pseudo-polynomial algorithms solve Subset Sum in polynomial time $n^{O(1)}$. Our main question is: When can dense Subset Sum be solved in near-linear time $\tilde{O}(n)$? We provide an essentially complete dichotomy by designing improved algorithms and proving conditional lower bounds, thereby determining essentially all settings of the parameters $n,t,\rm{mx}_X,\Sigma_X$ for which dense Subset Sum is in time $\tilde{O}(n)$. For notational convenience we assume without loss of generality that $t \ge \rm{mx}_X$ (as larger numbers can be ignored) and $t \le \Sigma_X/2$ (using symmetry). Then our dichotomy reads as follows: - By reviving and improving an additive-combinatorics-based approach by Galil and Margalit [SICOMP'91], we show that Subset Sum is in near-linear time $\tilde{O}(n)$ if $t \gg \rm{mx}_X \Sigma_X/n^2$. - We prove a matching conditional lower bound: If Subset Sum is in near-linear time for any setting with $t \ll \rm{mx}_X \Sigma_X/n^2$, then the Strong Exponential Time Hypothesis and the Strong k-Sum Hypothesis fail. We also generalize our algorithm from sets to multi-sets, albeit with non-matching upper and lower bounds.


title: Goodwillie tower of the norm functor

abstract: In this note we compare the fracture squares from genuine equivariant stable homotopy theory and the fracture squares which appear in the Goodwillie tower for the norm functor.


title: Real-time Quadrotor Navigation Through Planning in Depth Space in Unstructured Environments

abstract: This paper addresses the problem of real-time vision-based autonomous obstacle avoidance in unstructured environments for quadrotor UAVs. We assume that our UAV is equipped with a forward facing stereo camera as the only sensor to perceive the world around it. Moreover, all the computations are performed onboard. Feasible trajectory generation in this kind of problems requires rapid collision checks along with efficient planning algorithms. We propose a trajectory generation approach in the depth image space, which refers to the environment information as depicted by the depth images. In order to predict the collision in a look ahead robot trajectory, we create depth images from the sequence of robot poses along the path. We compare these images with the depth images of the actual world sensed through the forward facing stereo camera. We aim at generating fuel optimal trajectories inside the depth image space. In case of a predicted collision, a switching strategy is used to aggressively deviate the quadrotor away from the obstacle. For this purpose we use two closed loop motion primitives based on Linear Quadratic Regulator (LQR) objective functions. The proposed approach is validated through simulation and hardware experiments.


title: Decentralized and Secure Generation Maintenance with Differential Privacy

abstract: Decentralized methods are gaining popularity for data-driven models in power systems as they offer significant computational scalability while guaranteeing full data ownership by utility stakeholders. However, decentralized methods still require sharing information about network flow estimates over public facing communication channels, which raises privacy concerns. In this paper we propose a differential privacy driven approach geared towards decentralized formulations of mixed integer operations and maintenance optimization problems that protects network flow estimates. We prove strong privacy guarantees by leveraging the linear relationship between the phase angles and the flow. To address the challenges associated with the mixed integer and dynamic nature of the problem, we introduce an exponential moving average based consensus mechanism to enhance convergence, coupled with a control chart based convergence criteria to improve stability. Our experimental results obtained on the IEEE 118 bus case demonstrate that our privacy preserving approach yields solution qualities on par with benchmark methods without differential privacy. To demonstrate the computational robustness of our method, we conduct experiments using a wide range of noise levels and operational scenarios.


title: Two-photon amplitude interferometry for precision astrometry

abstract: Improved quantum sensing of photons from astronomical objects could provide high resolution observations in the optical benefiting numerous fields, including general relativity, dark matter studies, and cosmology. It has been recently proposed that stations in optical interferometers would not require a phase-stable optical link if instead sources of quantum-mechanically entangled pairs could be provided to them, potentially enabling hitherto prohibitively long baselines. A new refinement of this idea is developed, in which two photons from different sources are interfered at two separate and decoupled stations, requiring only a slow classical information link between them. We rigorously calculate the observables and contrast this new interferometric technique with the Hanbury Brown & Twiss intensity interferometry. We argue this technique could allow robust high-precision measurements of the relative astrometry of the two sources. A basic calculation suggests that angular precision on the order of $10$~microarcsecond in the relative opening angle could be achieved in a single night's observation of two bright stars.


title: The Convergence of AI code and Cortical Functioning -- a Commentary

abstract: Neural nets, one of the oldest architectures for AI programming, are loosely based on biological neurons and their properties. Recent work on language applications has made the AI code closer to biological reality in several ways. This commentary examines this convergence and, in light of what is known of neocortical structure, addresses the question of whether ``general AI'' looks attainable with these tools.


title: Variational Capsule Encoder

abstract: We propose a novel capsule network based variational encoder architecture, called Bayesian capsules (B-Caps), to modulate the mean and standard deviation of the sampling distribution in the latent space. We hypothesized that this approach can learn a better representation of features in the latent space than traditional approaches. Our hypothesis was tested by using the learned latent variables for image reconstruction task, where for MNIST and Fashion-MNIST datasets, different classes were separated successfully in the latent space using our proposed model. Our experimental results have shown improved reconstruction and classification performances for both datasets adding credence to our hypothesis. We also showed that by increasing the latent space dimension, the proposed B-Caps was able to learn a better representation when compared to the traditional variational auto-encoders (VAE). Hence our results indicate the strength of capsule networks in representation learning which has never been examined under the VAE settings before.


title: Unsupervised Foveal Vision Neural Networks with Top-Down Attention

abstract: Deep learning architectures are an extremely powerful tool for recognizing and classifying images. However, they require supervised learning and normally work on vectors the size of image pixels and produce the best results when trained on millions of object images. To help mitigate these issues, we propose the fusion of bottom-up saliency and top-down attention employing only unsupervised learning techniques, which helps the object recognition module to focus on relevant data and learn important features that can later be fine-tuned for a specific task. In addition, by utilizing only relevant portions of the data, the training speed can be greatly improved. We test the performance of the proposed Gamma saliency technique on the Toronto and CAT2000 databases, and the foveated vision in the Street View House Numbers (SVHN) database. The results in foveated vision show that Gamma saliency is comparable to the best and computationally faster. The results in SVHN show that our unsupervised cognitive architecture is comparable to fully supervised methods and that the Gamma saliency also improves CNN performance if desired. We also develop a topdown attention mechanism based on the Gamma saliency applied to the top layer of CNNs to improve scene understanding in multi-object images or images with strong background clutter. When we compare the results with human observers in an image dataset of animals occluded in natural scenes, we show that topdown attention is capable of disambiguating object from background and improves system performance beyond the level of human observers.


title: Quantum cellular automata and quantum field theory in two spatial dimensions

abstract: Quantum walks on lattices can give rise to one-particle relativistic wave equations in the long-wavelength limit. In going to multiple particles, quantum cellular automata (QCA) are natural generalizations of quantum walks. In one spatial dimension, the quantum walk can be "promoted" to a QCA that, in the long-wavelength limit, gives rise to the Dirac quantum field theory (QFT) for noninteracting fermions. This QCA/QFT correspondence has both theoretical and practical applications, but there are obstacles to similar constructions in two or more spatial dimensions. Here we show that a method of construction employing distinguishable particles confined to the completely antisymmetric subspace yields a QCA in two spatial dimensions that gives rise to the 2D Dirac QFT. Generalizing to 3D will entail some additional complications, but no conceptual barriers. We examine how this construction evades the "no go" results in earlier work.


title: Movement-induced Priors for Deep Stereo

abstract: We propose a method for fusing stereo disparity estimation with movement-induced prior information. Instead of independent inference frame-by-frame, we formulate the problem as a non-parametric learning task in terms of a temporal Gaussian process prior with a movement-driven kernel for inter-frame reasoning. We present a hierarchy of three Gaussian process kernels depending on the availability of motion information, where our main focus is on a new gyroscope-driven kernel for handheld devices with low-quality MEMS sensors, thus also relaxing the requirement of having full 6D camera poses available. We show how our method can be combined with two state-of-the-art deep stereo methods. The method either work in a plug-and-play fashion with pre-trained deep stereo networks, or further improved by jointly training the kernels together with encoder-decoder architectures, leading to consistent improvement.


title: Robust Learning under Strong Noise via SQs

abstract: This work provides several new insights on the robustness of Kearns' statistical query framework against challenging label-noise models. First, we build on a recent result by \cite{DBLP:journals/corr/abs-2006-04787} that showed noise tolerance of distribution-independently evolvable concept classes under Massart noise. Specifically, we extend their characterization to more general noise models, including the Tsybakov model which considerably generalizes the Massart condition by allowing the flipping probability to be arbitrarily close to $\frac{1}{2}$ for a subset of the domain. As a corollary, we employ an evolutionary algorithm by \cite{DBLP:conf/colt/KanadeVV10} to obtain the first polynomial time algorithm with arbitrarily small excess error for learning linear threshold functions over any spherically symmetric distribution in the presence of spherically symmetric Tsybakov noise. Moreover, we posit access to a stronger oracle, in which for every labeled example we additionally obtain its flipping probability. In this model, we show that every SQ learnable class admits an efficient learning algorithm with OPT + $\epsilon$ misclassification error for a broad class of noise models. This setting substantially generalizes the widely-studied problem of classification under RCN with known noise rate, and corresponds to a non-convex optimization problem even when the noise function -- i.e. the flipping probabilities of all points -- is known in advance.


title: Conformal prediction for time series

abstract: We develop a general framework for constructing distribution-free prediction intervals for time series. Theoretically, we establish explicit bounds on conditional and marginal coverage gaps of estimated prediction intervals, which asymptotically converge to zero under additional assumptions. We obtain similar bounds on the size of set differences between oracle and estimated prediction intervals. Methodologically, we introduce a computationally efficient algorithm called \texttt{EnbPI} that wraps around ensemble predictors, which is closely related to conformal prediction (CP) but does not require data exchangeability. \texttt{EnbPI} avoids data-splitting and is computationally efficient by avoiding retraining and thus scalable to sequentially producing prediction intervals. We perform extensive simulation and real-data analyses to demonstrate its effectiveness compared with existing methods. We also discuss the extension of \texttt{EnbPI} on various other applications.


title: Bridging the gap between Markowitz planning and deep reinforcement learning

abstract: While researchers in the asset management industry have mostly focused on techniques based on financial and risk planning techniques like Markowitz efficient frontier, minimum variance, maximum diversification or equal risk parity, in parallel, another community in machine learning has started working on reinforcement learning and more particularly deep reinforcement learning to solve other decision making problems for challenging task like autonomous driving, robot learning, and on a more conceptual side games solving like Go. This paper aims to bridge the gap between these two approaches by showing Deep Reinforcement Learning (DRL) techniques can shed new lights on portfolio allocation thanks to a more general optimization setting that casts portfolio allocation as an optimal control problem that is not just a one-step optimization, but rather a continuous control optimization with a delayed reward. The advantages are numerous: (i) DRL maps directly market conditions to actions by design and hence should adapt to changing environment, (ii) DRL does not rely on any traditional financial risk assumptions like that risk is represented by variance, (iii) DRL can incorporate additional data and be a multi inputs method as opposed to more traditional optimization methods. We present on an experiment some encouraging results using convolution networks.


title: Technical Brief: Finite Element Modeling of Tight Elastic Knots

abstract: We present a methodology to simulate the mechanics of knots in elastic rods using geometrically nonlinear, full three-dimensional (3D) finite element analysis. We focus on the mechanical behavior of knots in tight configurations, for which the full 3D deformation must be taken into account. To set up the topology of our knotted structures, we apply a sequence of prescribed displacement steps to the centerline of an initially straight rod that is meshed with 3D solid elements. Self-contact is enforced with a normal penalty force combined with Coulomb friction. As test cases, we investigate both overhand and figure-of-eight knots. Our simulations are validated with precision model experiments, combining rod fabrication and X-ray tomography. Even if the focus is given to the methods, our results reveal that 3D deformation of tight elastic knots is central to their mechanical response. These findings contrast to a previous analysis of loose knots, for which 1D centerline-based rod theories sufficed for a predictive understanding. Our method serves as a robust framework to access complex mechanical behavior of tightly knotted structures that are not readily available through experiments nor existing reduced-order theories.


title: Functional strong laws of large numbers for Euler characteristic processes of extreme sample clouds

abstract: To recover the topology of a manifold in the presence of heavy tailed or exponentially decaying noise, one must understand the behavior of geometric complexes whose points lie in the tail of these noise distributions. This study advances this line of inquiry, and demonstrates functional strong laws of large numbers for the Euler characteristic process of random geometric complexes formed by random points outside of an expanding ball in $\mathbb{R}^d$. When the points are drawn from a heavy tailed distribution with a regularly varying tail, the Euler characteristic process grows at a regularly varying rate, and the scaled process converges uniformly and almost surely to a smooth function. When the points are drawn from a distribution with an exponentially decaying tail, the Euler characteristic process grows logarithmically, and the scaled process converges to another smooth function in the same sense. All of the limit theorems take place when the points inside the expanding ball are densely distributed, so that the simplex counts outside of the ball of all dimensions contribute to the Euler characteristic process.


title: Quantifier completions, choice principles and applications

abstract: We contribute to the knowledge of the quantifier completions and their applications by using the language of doctrines. This algebraic presentation allows us to properly analyse the behaviour of the existential and universal quantifiers. We wish to convey the following points: the first is that these completions preserve the lattice structure and the distributive lattice structure of the fibres under opportune hypotheses which turn out to be preserved as well; the second regards the applications, in particular to the dialectica construction; the third is that these free constructions carry on some relevant choice principles.


title: BBB-Voting: 1-out-of-k Blockchain-Based Boardroom Voting

abstract: Voting is a means to agree on a collective decision based on available choices (e.g., candidates), where participants agree to abide by their outcome. To improve some features of e-voting, decentralized blockchain-based solutions can be employed, where the blockchain represents a public bulletin board that in contrast to a centralized bulletin board provides extremely high availability, censorship resistance, and correct code execution. A blockchain ensures that all entities in the voting system have the same view of the actions made by others due to its immutability and append-only features. The existing remote blockchain-based boardroom voting solution called Open Voting Network (OVN) provides the privacy of votes, universal & End-to-End verifiability, and perfect ballot secrecy; however, it supports only two choices and lacks robustness enabling recovery from stalling participants. We present BBB-Voting, an equivalent blockchain-based approach for decentralized voting such as OVN, but in contrast to it, BBB-Voting supports 1-out-of-$k$ choices and provides robustness that enables recovery from stalling participants. We make a cost-optimized implementation using an Ethereum-based environment respecting Ethereum Enterprise Alliance standards, which we compare with OVN and show that our work decreases the costs for voters by 13.5% in normalized gas consumption. Finally, we show how BBB-Voting can be extended to support the number of participants limited only by the expenses paid by the authority and the computing power to obtain the tally.


title: Disinformation in the Online Information Ecosystem: Detection, Mitigation and Challenges

abstract: With the rapid increase in access to internet and the subsequent growth in the population of online social media users, the quality of information posted, disseminated and consumed via these platforms is an issue of growing concern. A large fraction of the common public turn to social media platforms and in general the internet for news and even information regarding highly concerning issues such as COVID-19 symptoms. Given that the online information ecosystem is extremely noisy, fraught with misinformation and disinformation, and often contaminated by malicious agents spreading propaganda, identifying genuine and good quality information from disinformation is a challenging task for humans. In this regard, there is a significant amount of ongoing research in the directions of disinformation detection and mitigation. In this survey, we discuss the online disinformation problem, focusing on the recent 'infodemic' in the wake of the coronavirus pandemic. We then proceed to discuss the inherent challenges in disinformation research, and then elaborate on the computational and interdisciplinary approaches towards mitigation of disinformation, after a short overview of the various directions explored in detection efforts.


title: Free energy upper bound for mean-field vector spin glasses

abstract: We consider vector spin glasses whose energy function is a Gaussian random field with covariance given in terms of the matrix of scalar products. For essentially any model in this class, we give an upper bound for the limit free energy, which is expected to be sharp. The bound is expressed in terms of an infinite-dimensional Hamilton-Jacobi equation.


title: Shortest Paths Among Obstacles in the Plane Revisited

abstract: Given a set of pairwise disjoint polygonal obstacles in the plane, finding an obstacle-avoiding Euclidean shortest path between two points is a classical problem in computational geometry and has been studied extensively. The previous best algorithm was given by Hershberger and Suri [FOCS 1993, SIAM J. Comput. 1999] and the algorithm runs in $O(n\log n)$ time and $O(n\log n)$ space, where $n$ is the total number of vertices of all obstacles. The algorithm is time-optimal because $\Omega(n\log n)$ is a lower bound. It has been an open problem for over two decades whether the space can be reduced to $O(n)$. In this paper, we settle it by solving the problem in $O(n\log n)$ time and $O(n)$ space, which is optimal in both time and space; we achieve this by modifying the algorithm of Hershberger and Suri. Like their original algorithm, our new algorithm can build a shortest path map for a source point $s$ in $O(n\log n)$ time and $O(n)$ space, such that given any query point $t$, the length of a shortest path from $s$ to $t$ can be computed in $O(\log n)$ time and a shortest path can be produced in additional time linear in the number of edges of the path.


title: Topic Recommendation for Software Repositories using Multi-label Classification Algorithms

abstract: Many platforms exploit collaborative tagging to provide their users with faster and more accurate results while searching or navigating. Tags can communicate different concepts such as the main features, technologies, functionality, and the goal of a software repository. Recently, GitHub has enabled users to annotate repositories with topic tags. It has also provided a set of featured topics, and their possible aliases carefully curated with the help of the community. This creates the opportunity to use this initial seed of topics to automatically annotate all remaining repositories, by training models that recommend high-quality topic tags to developers. In this work, we study the application of multi-label classification techniques to predict software repositories' topics. First, we map the large space of user-defined topics to those featured by GitHub. The core idea is to derive more information from projects' available documentation. Our data contains about $152$K GitHub repositories and $228$ featured topics. Then, we apply supervised models on repositories' textual information such as descriptions, README files, wiki pages, and file names. We assess the performance of our approach both quantitatively and qualitatively. Our proposed model achieves Recall@5 and LRAP scores of $0.890$ and $0.805$, respectively. Moreover, based on users' assessment, our approach is highly capable of recommending a correct and complete set of topics. Finally, we use our models to develop an online tool named \texttt{Repository Catalogue}, that automatically predicts topics for GitHub repositories and is publicly available.