title: Predicting with High Correlation Features

abstract: It has been shown that instead of learning actual object features, deep networks tend to exploit non-robust (spurious) discriminative features that are shared between training and test sets. Therefore, while they achieve state of the art performance on such test sets, they achieve poor generalization on out of distribution (OOD) samples where the IID (independent, identical distribution) assumption breaks and the distribution of non-robust features shifts. In this paper, we consider distribution shift as a shift in the distribution of input features during test time that exhibit low correlation with targets in the training set. Under this definition, we evaluate existing robust feature learning methods and regularization methods and compare them against a baseline designed to specifically capture high correlation features in training set. As a controlled test-bed, we design a colored MNIST (C-MNIST) dataset and find that existing methods trained on this set fail to generalize well on an OOD version this dataset, showing that they overfit the low correlation color features. This is avoided by the baseline method trained on the same C-MNIST data, which is designed to learn high correlation features, and is able to generalize on the test sets of vanilla MNIST, MNIST-M and SVHN datasets. Our code is available at \url{https://github.com/salesforce/corr_based_prediction}.


title: Kloosterman sums over finite Frobenius rings

abstract: We study Kloosterman sums in a generalized ring-theoretic context, that of finite commutative Frobenius rings. We prove a number of identities for twisted Kloosterman sums, loosely clustered around moment computations.


title: Consistency Analysis of the Simplified Refined Instrumental Variable Method for Continuous-time Systems

abstract: In this paper, we analyse the consistency of the Simplified Refined Instrumental Variable method for Continuous-time systems (SRIVC). It is well known that the intersample behaviour of the input signal influences the quality and accuracy of the results when estimating and simulating continuous-time models. Here, we present a comprehensive analysis on the consistency of the SRIVC estimator while taking into account the intersample behaviour of the input signal. The main result of the paper shows that, under some mild conditions, the SRIVC estimator is generically consistent. We also describe some conditions when consistency is not achieved, which is important from a practical standpoint. The theoretical results are supported by simulation examples.


title: LRP2020: Signposts of planet formation in protoplanetary disks

abstract: Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.


title: Leaky Forcing: A New Variation of Zero Forcing

abstract: Zero forcing is a one-player game played on a graph. The player chooses some set of vertices to color, then iteratively applies a color change rule: If all but one of a colored vertex's neighbors are colored, color (i.e. "force") the remaining uncolored neighbor. Generally, the goal is to find the minimum number of vertices to initially color such that all vertices eventually become colored. Recently, equivalent formations of zero forcing have been developed in different settings including sensor allocation to solve linear systems (K.-Lin 2018), controllability in follower-leader dynamics (Monshizadeh-Zhang-Camlibel 2014), and edge covering in specific hypergraphs (Brimkov-Fast-Hicks 2016). While many variations of zero forcing are motivated by an associated minimum rank problem, these new formulations give new inspiration for new meaningful zero forcing variants. In our case, we study a new variation based on the linear algebraic interpretation mentioned above. In particular, what if there is a juncture in a network that has a leak, and, hence, is unreliable to facilitate solving a linear system on the network? In the context of zero forcing this corresponds to the following variation we call $\ell$-forcing: Given $\ell$, find a set of vertices such that for any set of $\ell$ vertices that are unable to force, all vertices will still be colored. We compute the $\ell$-forcing number for selected families of graphs including grid graphs. Perhaps surprisingly, we find examples where additional edges make the graph more "resilient" to these leaks. Further, we also implement known computational methods for our new leaky forcing variation.


title: The Next Generation Virgo Cluster Survey. XVII. A Search for Planetary Nebulae in Virgo Cluster Globular Clusters

abstract: The occurrence of planetary nebulae (PNe) in globular clusters (GCs) provides an excellent chance to study low-mass stellar evolution in a special (low-metallicity, high stellar density) environment. We report a systematic spectroscopic survey for the [O{\sc iii}] 5007 emission line of PNe in 1469 Virgo GCs and 121 Virgo ultra-compact dwarfs (UCDs), mainly hosted in the giant elliptical galaxies M87, M49, M86, and M84. We detected zero PNe in our UCD sample and discovered one PN ($M_{5007} = -4.1$ mag) associated with an M87 GC. We used the [O{\sc iii}] detection limit for each GC to estimate the luminosity-specific frequency of PNe, $\alpha$, and measured $\alpha$ in the Virgo cluster GCs to be $\alpha \sim 3.9_{-0.7}^{+5.2}\times 10^{-8}\mathrm{PN}/L_\odot$. $\alpha$ in Virgo GCs is among the lowest values reported in any environment, due in part to the large sample size, and is 5--6 times lower than that for the Galactic GCs. We suggest that $\alpha$ decreases towards brighter and more massive clusters, sharing a similar trend as the binary fraction, and the discrepancy between the Virgo and Galactic GCs can be explained by the observational bias in extragalactic surveys toward brighter GCs. This low but non-zero efficiency in forming PNe may highlight the important role played by binary interactions in forming PNe in GCs. We argue that a future survey of less massive Virgo GCs will be able to determine whether PN production in Virgo GCs is governed by internal process (mass, density, binary fraction), or is largely regulated by external environment.


title: How To Catch A Lion In The Desert -- On The Solution Of The Coverage Directed Generation (CDG) Problem

abstract: The testing and verification of a complex hardware or software system, such as modern integrated circuits (ICs) found in everything from smartphones to servers, can be a difficult process. One of the most difficult and time-consuming tasks a verification team faces is reaching coverage closure, or hitting all events in the coverage space. Coverage-directed-generation (CDG), or the automatic generation of tests that can hit hard-to-hit coverage events, and thus provide coverage closure, holds the potential to save verification teams significant simulation resources and time. In this paper, we propose a new approach to the CDG problem by formulating the CDG problem as a noisy derivative free optimization (DFO) problem. However, this formulation is complicated by the fact that derivatives of the objective function are unavailable, and the objective function evaluations are corrupted by noise. We solve this noisy optimization problem by utilizing techniques from direct optimization coupled with a robust noise estimator, and by leveraging techniques from inverse problems to estimate the gradient of the noisy objective function. We demonstrate the efficiency and reliability of this new approach through numerical experiments with an abstract model of part of IBM's NorthStar processor, a superscalar in-order processor designed for servers.


title: A novel compound synapse using probabilistic spin-orbit-torque switching for MTJ based deep neural networks

abstract: Analog electronic non-volatile memories mimicking synaptic operations are being explored for the implementation of neuromorphic computing systems. Compound synapses consisting of ensembles of stochastic binary elements are alternatives to analog memory synapses to achieve multilevel memory operation. Among existing binary memory technologies, magnetic tunneling junction (MTJ) based Magnetic Random Access Memory (MRAM) technology has matured to the point of commercialization. More importantly for this work, stochasticity is natural to the MTJ switching physics e.g devices referred as p-bits which mimic binary stochastic neurons. In this article, we experimentally demonstrate a novel compound synapse that uses stochastic spin-orbit torque (SOT) switching of an ensemble of nano-magnets that are located on one shared spin Hall effect (SHE) material channel, i.e. tantalum. By using a properly chosen pulse scheme, we are able to demonstrate linear potentiation and depression in the synapse, as required for many neuromorphic architectures. In addition to this experimental effort, we also performed circuit simulations on an SOT-MRAM based 784*200*10 deep belief network (DBN) consisting of p-bit based neurons and compound synapses. MNIST pattern recognition was used to evaluate the system performance, and our findings indicate that a significant reduction in recognition error rates can be achieved when using our incremental pulse scheme rather than a non-linear potentiation and depression as obtained when employing identical pulses.


title: Unconditionally energy stable DG schemes for the Swift-Hohenberg equation

abstract: The Swift-Hohenberg equation as a central nonlinear model in modern physics has a gradient flow structure. Here we introduce fully discrete discontinuous Galerkin (DG) schemes for a class of fourth order gradient flow problems, including the nonlinear Swift-Hohenberg equation, to produce free-energy-decaying discrete solutions, irrespective of the time step and the mesh size. We exploit and extend the mixed DG method introduced in [H. Liu and P. Yin, J. Sci. Comput., 77: 467--501, 2018] for the spatial discretization, and the "Invariant Energy Quadratization" method for the time discretization. The resulting IEQ-DG algorithms are linear, thus they can be efficiently solved without resorting to any iteration method. We actually prove that these schemes are unconditionally energy stable. We present several numerical examples that support our theoretical results and illustrate the efficiency, accuracy and energy stability of our new algorithm. The numerical results on two dimensional pattern formation problems indicate that the method is able to deliver comparable patterns of high accuracy.


title: Finite time blowup of 2D Boussinesq and 3D Euler equations with $C^{1,\alpha}$ velocity and boundary

abstract: Inspired by the numerical evidence of a potential 3D Euler singularity by Luo-Hou [30,31] and the recent breakthrough by Elgindi [11] on the singularity formation of the 3D Euler equation without swirl with $C^{1,\alpha}$ initial velocity, we prove the finite time singularity for the 2D Boussinesq and the 3D axisymmetric Euler equations with boundary and $C^{1,\alpha}$ initial data for the velocity (and density in the case of Boussinesq equations). Our finite time blowup solution for the 3D Euler equations and the singular solution considered in [30,31] share many essential features, including the symmetry properties of the solution, the flow structure, and the sign of the solution in each quadrant, except that we use $C^{1,\alpha}$ initial velocity. We use a dynamic rescaling formulation and follow the general framework of analysis developed in [11]. We also use some strategy proposed in our recent joint work with Huang in [7] and adopt several methods of analysis in [11] to establish the linear and nonlinear stability of an approximate self-similar profile. The nonlinear stability enables us to prove that the solution of the 3D Euler or the 2D Boussinesq equations with $C^{1,\alpha}$ initial data will develop a finite time singularity. Moreover, the velocity field has finite energy before the singularity time. In the previous version of this paper, we proved the blowup results for the 3D axisymmetric Euler equations with initial data $(u_0^{\theta})^2, u_0^r, u_0^z \in C^{1,\alpha}$. Though the velocity $u^r, u^z$ in the axisymmetric setting is $C^{1,\alpha}$, our interpretation that the velocity is $C^{1,\alpha}$ is not correct since the velocity in 3D also depends on $u^{\theta}$, which is not $C^{1,\alpha}$. This oversight can be fixed easily with a minor change in the construction of the approximate steady state and a minor modification to localize the approximate steady state.


title: Randomized Ablation Feature Importance

abstract: Given a model $f$ that predicts a target $y$ from a vector of input features $\pmb{x} = x_1, x_2, \ldots, x_M$, we seek to measure the importance of each feature with respect to the model's ability to make a good prediction. To this end, we consider how (on average) some measure of goodness or badness of prediction (which we term "loss" $\ell$), changes when we hide or ablate each feature from the model. To ablate a feature, we replace its value with another possible value randomly. By averaging over many points and many possible replacements, we measure the importance of a feature on the model's ability to make good predictions. Furthermore, we present statistical measures of uncertainty that quantify how confident we are that the feature importance we measure from our finite dataset and finite number of ablations is close to the theoretical true importance value.


title: Higher Orbit Integrals, Cyclic Cocyles, and K-theory of Reduced Group C*-algebra

abstract: Let G be a connected real reductive group. Orbit integrals define traces on the group algebra of G. We introduce a construction of higher orbit integrals in the direction of higher cyclic cocycles on the Harish-Chandra Schwartz algebra of G. We analyze these higher orbit integrals via Fourier transform by expressing them as integrals on the tempered dual of G. We obtain explicit formulas for the pairing between the higher orbit integrals and the K-theory of the reduced group C*-algebra, and discuss their applications to representation theory and K-theory.


title: Oxygen deficiency and migration mediated electric polarization in Fe,Co-substituted SrTiO$_{3-\delta}$

abstract: We use density functional theory (DFT) calculations to show that oxygen vacancies ($v_\mathrm{O}$) and mobility induce noncentrosymmetric polar structures in SrTi$_{1-x-y}$Fe$_{x}$Co$_{y}$O$_{3-\delta}$ ($x=y=0.125$) with $\delta = \{0.125, 0.25\}$, enhance the saturation magnetization and give rise to large changes in the electric polarization $\vert\Delta P\vert$. We present an intuitive set of rules for SrTiFeCoO$_{3-\delta}$ (STFC), which are based on the interplay between (Co/Fe)-$v_\mathrm{O}$ defects, magnetic cations coordination and topological vacancy disorder. STFC structures convey layered crystals with sheets of linear organized O$_{4,5,6}$-coordinated Fe-Co pairs, sandwiched with layers of O$_{5}$-coordinated Ti. Co,Fe-$v_\mathrm{O}$ defects are the source of the crystal distortions, cations off-centering and bending of the oxygen octahedra, which added to the charge redistribution mediated by $v_\mathrm{O}$, the cations electronegativity and valence states trigger an effective electric polarization. Oxygen migrations for $\delta=0.125$ provides us with $\vert\Delta \mathbf{P}\vert$ $>\sim10 \mu$C/cm$^2$ due to a quantum-of-polarization differences between $\delta=0.125$ structures. Increasing the deficiency to $\delta=0.25$ yields $\vert\Delta \mathbf{P}\vert$ whose O-migration resolved polarization for $\delta=0.25$ is $>\sim3 \mu$C/cm$^2$ in the worst case scenario. Magnetism is dominated by the Fe,Co spin states for $\delta=0.125$ while there is a raid of Ti magnetic moments ($\sim1\mu_{B}$) for $\delta=0.25$. Magnetic and electric order parameters change for variations of $\delta$ or oxygen migrations for a given deficiency. Our results capture characteristics observed in the end-members of the series SrTi(Co,Fe)O$_{3}$, and suggest the existence of a broader set of rules for oxygen deficient multiferroic oxides.


title: Advantage-Weighted Regression: Simple and Scalable Off-Policy Reinforcement Learning

abstract: In this paper, we aim to develop a simple and scalable reinforcement learning algorithm that uses standard supervised learning methods as subroutines. Our goal is an algorithm that utilizes only simple and convergent maximum likelihood loss functions, while also being able to leverage off-policy data. Our proposed approach, which we refer to as advantage-weighted regression (AWR), consists of two standard supervised learning steps: one to regress onto target values for a value function, and another to regress onto weighted target actions for the policy. The method is simple and general, can accommodate continuous and discrete actions, and can be implemented in just a few lines of code on top of standard supervised learning methods. We provide a theoretical motivation for AWR and analyze its properties when incorporating off-policy data from experience replay. We evaluate AWR on a suite of standard OpenAI Gym benchmark tasks, and show that it achieves competitive performance compared to a number of well-established state-of-the-art RL algorithms. AWR is also able to acquire more effective policies than most off-policy algorithms when learning from purely static datasets with no additional environmental interactions. Furthermore, we demonstrate our algorithm on challenging continuous control tasks with highly complex simulated characters.


title: NGEMM: Optimizing GEMM for Deep Learning via Compiler-based Techniques

abstract: Quantization has emerged to be an effective way to significantly boost the performance of deep neural networks (DNNs) by utilizing low-bit computations. Despite having lower numerical precision, quantized DNNs are able to reduce both memory bandwidth and computation cycles with little losses of accuracy. Integer GEMM (General Matrix Multiplication) is critical to running quantized DNN models efficiently, as GEMM operations often dominate the computations in these models. Various approaches have been developed by leveraging techniques such as vectorization and memory layout to improve the performance of integer GEMM. However, these existing approaches are not fast enough in certain scenarios. We developed NGEMM, a compiler-based GEMM implementation for accelerating lower-precision training and inference. NGEMM has better use of the vector units by avoiding unnecessary vector computation that is introduced during tree reduction. We compared NGEMM's performance with the state-of-art BLAS libraries such as MKL. Our experimental results showed that NGEMM outperformed MKL non-pack and pack version by an average of 1.86x and 1.16x, respectively. We have applied NGEMM to a number of production services in Microsoft.


title: Thermal Conductivity and Magnetic Phase Diagram of CuB2O4

abstract: We have measured temperature and magnetic field dependences of the thermal conductivity along the c-axis, kc, and that along the [110] direction, k110, of CuB2O4 single crystals in zero field and magnetic fields along the c-axis and along the [110] direction. It has been found that the thermal conductivity is nearly isotropic and very large in zero field and that the thermal conductivity due to phonons is dominant in CuB2O4. The temperature and field dependences of kc and k110 have markedly changed at phase boundaries in the magnetic phase diagram, which has been understood to be due to the change of the mean free path of phonons caused by the change of the phonon-spin scattering rate at the phase boundaries. It has been concluded that thermal conductivity measurements are very effective for detecting magnetic phase boundaries.


title: Photometry of the Uranian Satellites with Keck and the Search for Mab

abstract: We present photometric properties of six small (radii < 100 km) satellites of Uranus based on 32 H-(1.49-1.78 um) band images taken on August 29, 2015 from the Keck II Telescope on Maunakea, Hawaii with the near-infrared camera NIRC2 coupled to the adaptive optics system. The sub-observer latitude of our observations was 31 deg, i.e., we view much of the satellites' north poles, in contrast to the 1986 Voyager measurements. We derive reflectivities based on mean-stacking measurements of these six minor moons of Uranus. We find that the small satellites are significantly brighter than in previous observations, which we attribute to albedo variations between hemispheres. We also search for Mab, a small satellite with an unknown surface composition, orbiting between Puck and Miranda. Despite the significantly improved signal-to-noise ratio we achieved, we could not detect Mab. We suggest that Mab is more similar to Miranda, an icy body, than to the inner rocky moons. Assuming Mab is spherical with a radius of 6 km, as derived from Hubble Space Telescope (HST) observations if its reflectivity is ~0.46, we derive a 3{\sigma} upper limit to its reflectivity of 0.14 [I/F] at 1.6 um.


title: Rigid connections on $\mathbb{P}^1$ via the Bruhat-Tits building

abstract: We apply the theory of fundamental strata of Bremer and Sage to find cohomologically rigid $G$-connections on the projective line, generalising the work of Frenkel and Gross. In this theory, one studies the leading term of a formal connection with respect to the Moy-Prasad filtration associated to a point in the Bruhat-Tits building. If the leading term is regular semisimple with centraliser a (not necessarily split) maximal torus $S$, then we have an $S$-toral connection. In this language, the irregular singularity of the Frenkel-Gross connection gives rise to the homogenous toral connection of minimal slope associated to the Coxeter torus $\mathcal{C}$. In the present paper, we consider connections on $\mathbb{G}_m$ which have an irregular homogeneous $\mathcal{C}$-toral singularity at zero of slope $i/h$, where $h$ is the Coxeter number and $i$ is a positive integer coprime to $h$, and a regular singularity at infinity with unipotent monodromy. Our main result is the characterisation of all such connections which are rigid.


title: Atomic Imaging of Mechanically Induced Topological Transition of Ferroelectric Vortices

abstract: Ferroelectric vortices formed through complex lattice-charge interactions have great potential in applications for future nanoelectronics such as memories. For practical applications, it is crucial to manipulate these topological states under external stimuli. Here, we apply mechanical loads to locally manipulate the vortices in a PbTiO3-SrTiO3 superlattice via atomically resolved in situ scanning transmission electron microscopy. The vortices undergo a transition to the a-domain with in-plane polarization under external compressive stress and spontaneously recover after removal of the stress. We reveal the detailed transition process at the atomic scale and reproduce this numerically using phase-field simulations. These findings provide new pathways to control the exotic topological ferroelectric structures for future nanoelectronics and also valuable insights into understanding of lattice-charge interactions at nanoscale.


title: Bearing-Only Consensus and Formation Control under Directed Topologies

abstract: We address the problems of bearing-only consensus and formation control, where each agent can only measure the relative bearings of its neighbors and relative distances are not available. We provide stability results for the Filippov solutions of two gradient-descent laws from non-smooth Lyapunov functions in the context of differential inclusion. For the consensus and formation control problems with undirected sensing topologies, we prove finite-time and asymptotic convergence of the proposed non-smooth gradient flows. For the directed consensus problem, we prove asymptotic convergence using a different non-smooth Lyapunov function given that the sensing graph has a globally reachable node. Finally, For the directed formation control problem we prove asymptotic convergence for directed cycles and directed acyclic graphs and also introduce a new notion of bearing persistence which guarantees convergence to the desired bearings.


title: Weyl nodes and magnetostructural instability in antiperovskite Mn$_3$ZnC

abstract: The ferromagnetic phase of the cubic antiperovskite Mn$_3$ZnC is suggested from first-principles calculation to be a nodal line Weyl semimetal. Features in the electronic structure that are the hallmark of a nodal line Weyl state, a large density of linear band crossings near the Fermi level, can also be interpreted as signatures of a structural and/or magnetic instability. Indeed, it is known that Mn$_3$ZnC undergoes transitions upon cooling from a paramagnetic to a cubic ferromagnetic state under ambient conditions and then further into a non-collinear ferrimagnetic tetragonal phase at a temperature between 250$\,$K and 200$\,$K. The existence of Weyl nodes and their destruction via structural and magnetic ordering is likely to be relevant to a range of magnetostructurally coupled materials.


title: Predicting Alzheimer's Disease by Hierarchical Graph Convolution from Positron Emission Tomography Imaging

abstract: Imaging-based early diagnosis of Alzheimer Disease (AD) has become an effective approach, especially by using nuclear medicine imaging techniques such as Positron Emission Topography (PET). In various literature it has been found that PET images can be better modeled as signals (e.g. uptake of florbetapir) defined on a network (non-Euclidean) structure which is governed by its underlying graph patterns of pathological progression and metabolic connectivity. In order to effectively apply deep learning framework for PET image analysis to overcome its limitation on Euclidean grid, we develop a solution for 3D PET image representation and analysis under a generalized, graph-based CNN architecture (PETNet), which analyzes PET signals defined on a group-wise inferred graph structure. Computations in PETNet are defined in non-Euclidean, graph (network) domain, as it performs feature extraction by convolution operations on spectral-filtered signals on the graph and pooling operations based on hierarchical graph clustering. Effectiveness of the PETNet is evaluated on the Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset, which shows improved performance over both deep learning and other machine learning-based methods.


title: Matching complexes of polygonal line tilings

abstract: The matching complex of a simple graph $G$ is a simplicial complex consisting of the matchings on $G$. Jeli\'c Milutinovi\'c et al. studied the matching complexes of the polygonal line tilings, and they gave a lower bound for the connectivity of the matching complexes of polygonal line tilings. In this paper, we determine the homotopy types of the matching complexes of polygonal line tilings recursively, and determine their connectivities.


title: A model of debt with bankruptcy risk and currency devaluation

abstract: The paper studies a system of Hamilton-Jacobi equations, arising from a stochastic optimal debt management problem in an infinite time horizon with exponential discount, modeled as a noncooperative interaction between a borrower and a pool of risk-neutral lenders. In this model, the borrower is a sovereign state that can decide how much to devaluate its currency and which fraction of its income should be used to repay the debt. Moreover, the borrower has the possibility of going bankrupt at a random time and must declare bankruptcy if the debt reaches a threshold x*. When bankruptcy occurs, the lenders only recover a fraction of their capital. To offset the possible loss of part of their investment, the lenders buy bonds at a discounted price which is not given a priori. This leads to a nonstandard optimal control problem. We establish an existence result of solutions to this system and in turn recover optimal feedback payment strategy $u*(x)$ and currency devaluation $v*(x)$. In addition, the behavior of $(u*,v*)$ near $0$ and x* is studied.


title: Beyond Textual Issues: Understanding the Usage and Impact of GitHub Reactions

abstract: Recently, GitHub introduced a new social feature, named reactions, which are "pictorial characters" similar to emoji symbols widely used nowadays in text-based communications. Particularly, GitHub users can use a pre-defined set of such symbols to react to issues and pull requests. However, little is known about the real usage and impact of GitHub reactions. In this paper, we analyze the reactions provided by developers to more than 2.5 million issues and 9.7 million issue comments, in order to answer an extensive list of nine research questions about the usage and adoption of reactions. We show that reactions are being increasingly used by open source developers. Moreover, we also found that issues with reactions usually take more time to be handled and have longer discussions.


title: The Non-IID Data Quagmire of Decentralized Machine Learning

abstract: Many large-scale machine learning (ML) applications need to perform decentralized learning over datasets generated at different devices and locations. Such datasets pose a significant challenge to decentralized learning because their different contexts result in significant data distribution skew across devices/locations. In this paper, we take a step toward better understanding this challenge by presenting a detailed experimental study of decentralized DNN training on a common type of data skew: skewed distribution of data labels across devices/locations. Our study shows that: (i) skewed data labels are a fundamental and pervasive problem for decentralized learning, causing significant accuracy loss across many ML applications, DNN models, training datasets, and decentralized learning algorithms; (ii) the problem is particularly challenging for DNN models with batch normalization; and (iii) the degree of data skew is a key determinant of the difficulty of the problem. Based on these findings, we present SkewScout, a system-level approach that adapts the communication frequency of decentralized learning algorithms to the (skew-induced) accuracy loss between data partitions. We also show that group normalization can recover much of the accuracy loss of batch normalization.


title: The energy-spectrum of bicompatible sequences

abstract: Background: Genotype-phenotype maps provide a meaningful filtration of sequence space and RNA secondary structures are particular such phenotypes. Compatible sequences i.e.~sequences that satisfy the base pairing constraints of a given RNA structure play an important role in the context of neutral networks and inverse folding. Sequences satisfying the constraints of two structures simultaneously are called bicompatible and phenotypic change, induced by erroneously replicating populations of RNA sequences, is closely connected to bicompatibility. Furthermore, bicompatible sequences are relevant for riboswitch sequences, beacons of evolution, realizing two distinct phenotypes. Results: We present a full loop energy model Boltzmann sampler of bicompatible sequences for pairs of structures. The novel dynamic programming algorithm is based on a topological framework encapsulating the relations between loops. We utilize our sequence sampler to study the energy spectra and density of bicompatible sequences, the rankings of the structures and key properties for evolutionary transitions. Conclusion: Our analysis of riboswitch sequences shows that key properties of bicompatible sequences depend on the particular pair of structures. While there always exist bicompatible sequences for random structure pairs, they are less suited to facilitate transitions. We show that native riboswitch sequences exhibit a distinct signature with regards to the ranking of their two phenotypes relative to the minimum free energy, suggesting a new criterion for identifying native sequences and sequences subjected to evolutionary pressure.


title: Continuous integration logic

abstract: We combine continuous and integral logics and found a logical framework for metric measure spaces equipped with a family of continuous relations and operations. We prove the ultraproduct theorem and deduce compactness and other usual results. We also give applications of the compactness theorem in metric measure theory.


title: Writing habits and telltale neighbors: analyzing clinical concept usage patterns with sublanguage embeddings

abstract: Natural language processing techniques are being applied to increasingly diverse types of electronic health records, and can benefit from in-depth understanding of the distinguishing characteristics of medical document types. We present a method for characterizing the usage patterns of clinical concepts among different document types, in order to capture semantic differences beyond the lexical level. By training concept embeddings on clinical documents of different types and measuring the differences in their nearest neighborhood structures, we are able to measure divergences in concept usage while correcting for noise in embedding learning. Experiments on the MIMIC-III corpus demonstrate that our approach captures clinically-relevant differences in concept usage and provides an intuitive way to explore semantic characteristics of clinical document collections.


title: Parallel Algorithm for Approximating Nash Equilibrium in Multiplayer Stochastic Games with Application to Naval Strategic Planning

abstract: Many real-world domains contain multiple agents behaving strategically with probabilistic transitions and uncertain (potentially infinite) duration. Such settings can be modeled as stochastic games. While algorithms have been developed for solving (i.e., computing a game-theoretic solution concept such as Nash equilibrium) two-player zero-sum stochastic games, research on algorithms for non-zero-sum and multiplayer stochastic games is limited. We present a new algorithm for these settings, which constitutes the first parallel algorithm for multiplayer stochastic games. We present experimental results on a 4-player stochastic game motivated by a naval strategic planning scenario, showing that our algorithm is able to quickly compute strategies constituting Nash equilibrium up to a very small degree of approximation error.


title: Improved Word Sense Disambiguation Using Pre-Trained Contextualized Word Representations

abstract: Contextualized word representations are able to give different representations for the same word in different contexts, and they have been shown to be effective in downstream natural language processing tasks, such as question answering, named entity recognition, and sentiment analysis. However, evaluation on word sense disambiguation (WSD) in prior work shows that using contextualized word representations does not outperform the state-of-the-art approach that makes use of non-contextualized word embeddings. In this paper, we explore different strategies of integrating pre-trained contextualized word representations and our best strategy achieves accuracies exceeding the best prior published accuracies by significant margins on multiple benchmark WSD datasets. We make the source code available at https://github.com/nusnlp/contextemb-wsd.


title: How noise affects the Hessian spectrum in overparameterized neural networks

abstract: Stochastic gradient descent (SGD) forms the core optimization method for deep neural networks. While some theoretical progress has been made, it still remains unclear why SGD leads the learning dynamics in overparameterized networks to solutions that generalize well. Here we show that for overparameterized networks with a degenerate valley in their loss landscape, SGD on average decreases the trace of the Hessian of the loss. We also generalize this result to other noise structures and show that isotropic noise in the non-degenerate subspace of the Hessian decreases its determinant. In addition to explaining SGDs role in sculpting the Hessian spectrum, this opens the door to new optimization approaches that may confer better generalization performance. We test our results with experiments on toy models and deep neural networks.


title: Unique Dirac and Triple point fermiology in simple transition metals and their binary alloys

abstract: Noble metal surfaces (Au, Ag and Cu etc.) have been extensively studied for the Shockley type surface states (SSs). Very recently, some of these Shockley SSs have been understood from the topological consideration, with the knowledge of global properties of electronic structure. In this letter, we show the existence of Dirac like excitations in the elemental noble metal Ru, Re and Os based on symmetry analysis and first principle calculations. The unique SSs driven Fermi arcs have been investigated in details for these metals. Our calculated SSs and Fermi arcs are consistent with the previous transport and photo-emission results. We attribute these Dirac excitation mediated Fermi arc topology to be the possible reasons behind several existing transport anomalies, such as large non-saturating magneto resistance, anomalous Nernst electromotive force and its giant oscillations, magnetic breakdown etc. We further show that the Dirac like excitations in these elemental metal can further be tuned to three component Fermionic excitations, using symmetry allowed alloy mechanism.


title: UltraShare: FPGA-based Dynamic Accelerator Sharing and Allocation

abstract: Despite all the available commercial and open-source frameworks to ease deploying FPGAs in accelerating applications, the current schemes fail to support sharing multiple accelerators among various applications. There are three main features that an accelerator sharing scheme requires to support: exploiting dynamic parallelism of multiple accelerators for a single application, sharing accelerators among multiple applications, and providing a non-blocking congestion-free environment for applications to invoke the accelerators. In this paper, we developed a scalable fully functional hardware controller, called UltraShare, with a supporting software stack that provides a dynamic accelerator sharing scheme through an accelerators grouping mechanism. UltraShare allows software applications to fully utilize FPGA accelerators in a non-blocking congestion-free environment. Our experimental results for a simple scenario of a combination of three streaming accelerators invocation show an improvement of up to 8x in throughput of the accelerators by removing accelerators idle times.


title: Hidden Supersymmetries of Deformed Supersymmetric Mechanics

abstract: We consider quantum models corresponding to superymmetrizations of the two-dimensional harmonic oscillator based on worldline $d=1$ realizations of the supergroup SU$(\,{\cal N}/2\,|1)$, where the number of supersymmetries ${\cal N}$ is arbitrary even number. Constructed models possess the hidden supersymmetry SU$(\,{\cal N}/2\,|2)$. Degeneracies of energy levels are spanned by representations of the hidden supersymmetry group.


title: Saliency is a Possible Red Herring When Diagnosing Poor Generalization

abstract: Poor generalization is one symptom of models that learn to predict target variables using spuriously-correlated image features present only in the training distribution instead of the true image features that denote a class. It is often thought that this can be diagnosed visually using attribution (aka saliency) maps. We study if this assumption is correct. In some prediction tasks, such as for medical images, one may have some images with masks drawn by a human expert, indicating a region of the image containing relevant information to make the prediction. We study multiple methods that take advantage of such auxiliary labels, by training networks to ignore distracting features which may be found outside of the region of interest. This mask information is only used during training and has an impact on generalization accuracy depending on the severity of the shift between the training and test distributions. Surprisingly, while these methods improve generalization performance in the presence of a covariate shift, there is no strong correspondence between the correction of attribution towards the features a human expert has labelled as important and generalization performance. These results suggest that the root cause of poor generalization may not always be spatially defined, and raise questions about the utility of masks as "attribution priors" as well as saliency maps for explainable predictions.


title: Monitoring the Superorbital Period Variation and the Spin Period Evolution of SMC X-1

abstract: The X-ray pulsar SMC X-1 shows a superorbital modulation with an unstable cycle length in the X-ray band. We present its timing behaviors, including the spin, orbital, and superorbital modulations, beyond the end of the Rossi X-ray Timing Explorer mission. The time-frequency maps derived by the wavelet Z-transform and the Hilbert-Huang transform suggest that a new superorbital period excursion event occurred in ~MJD 57100 (2015 March). This indicates the excursion is recurrent and probably (quasi)periodic. The hardness ratio obtained with the Monitor of All-sky X-ray Image (MAXI) suggests increased absorption during the transition from the high to the low state in the superorbital cycle. Compared to the regular epochs, the superorbital profile during the excursion epochs has a shallower and narrower valley, likely caused by a flatter warp. By tracking the spin period evolution with the MAXI Gas Slit Camera in 2-20 keV, we derive an averaged spin-up rate of $\dot{\nu}=2.515(3)\times10^{-11}$ s$^{-2}$ during the period between MJD 55141 (2009 November) and 58526 (2019 February). We obtain no positive correlation between the spin frequency residual and the superorbital frequency, but a torque change accompanying the superorbital period excursion is possible. We suggest that the accretion torque on the neutron star could be changed by various mechanisms, including the change of mass accretion rate and the warp angle. We update the value of the orbital decay as $\dot{P}_{\rm{orb}}/P_{\rm{orb}}=-3.380(6)\times10^{-6}$ yr$^{-1}$. Finally, we reconfirm the detection of the superorbital modulation in the optical band and its coherence in phase with the X-ray modulation.


title: Blending Diverse Physical Priors with Neural Networks

abstract: Machine learning in context of physical systems merits a re-examination of the learning strategy. In addition to data, one can leverage a vast library of physical prior models (e.g. kinematics, fluid flow, etc) to perform more robust inference. The nascent sub-field of \emph{physics-based learning} (PBL) studies the blending of neural networks with physical priors. While previous PBL algorithms have been applied successfully to specific tasks, it is hard to generalize existing PBL methods to a wide range of physics-based problems. Such generalization would require an architecture that can adapt to variations in the correctness of the physics, or in the quality of training data. No such architecture exists. In this paper, we aim to generalize PBL, by making a first attempt to bring neural architecture search (NAS) to the realm of PBL. We introduce a new method known as physics-based neural architecture search (PhysicsNAS) that is a top-performer across a diverse range of quality in the physical model and the dataset.


title: Theta series and number fields: theorems and experiments

abstract: We construct certain $\theta$-series associated to number fields and prove that for number fields of degree less than equal to 4, these $\theta$-series are number field invariants. We also investigate whether or not the collection of $\theta$-series associated to number fields of the same degree and discriminant are linearly independent. This is known to be true if the degree of the number field is less than or equal to 3. We do not prove in this paper that they are linearly independent in general but we do give computational and heuristic evidence that we would expect them to be.


title: Analyzing Sentence Fusion in Abstractive Summarization

abstract: While recent work in abstractive summarization has resulted in higher scores in automatic metrics, there is little understanding on how these systems combine information taken from multiple document sentences. In this paper, we analyze the outputs of five state-of-the-art abstractive summarizers, focusing on summary sentences that are formed by sentence fusion. We ask assessors to judge the grammaticality, faithfulness, and method of fusion for summary sentences. Our analysis reveals that system sentences are mostly grammatical, but often fail to remain faithful to the original article.


title: TriMap: Large-scale Dimensionality Reduction Using Triplets

abstract: We introduce "TriMap"; a dimensionality reduction technique based on triplet constraints, which preserves the global structure of the data better than the other commonly used methods such as t-SNE, LargeVis, and UMAP. To quantify the global accuracy of the embedding, we introduce a score that roughly reflects the relative placement of the clusters rather than the individual points. We empirically show the excellent performance of TriMap on a large variety of datasets in terms of the quality of the embedding as well as the runtime. On our performance benchmarks, TriMap easily scales to millions of points without depleting the memory and clearly outperforms t-SNE, LargeVis, and UMAP in terms of runtime.


title: Optimizing Information Freshness using Low-Power Status Updates via Sleep-Wake Scheduling

abstract: In this paper, we consider the problem of optimizing the freshness of status updates that are sent from a large number of low-power source nodes to a common access point. The source nodes utilize carrier sensing to reduce collisions and adopt an asychronized sleep-wake strategy to achieve an extended battery lifetime (e.g., 10-25 years). We use age of information (AoI) to measure the freshness of status updates, and design the sleep-wake parameters for minimizing the weighted-sum peak AoI of the sources, subject to per-source battery lifetime constraints. When the sensing time is zero, this sleep-wake design problem can be solved by resorting to nested convex optimization; however, for positive sensing times, the problem is non-convex. We devise a low-complexity solution to solve this problem and prove that, for practical sensing times, the solution is within a small gap from the optimum AoI performance. Our numerical and NS-3 simulation results show that our solution can indeed elongate the batteries lifetime of information sources, while providing a competitive AoI performance.


title: Classification of toric log del Pezzo surfaces with few singular points

abstract: We give a classification of toric log del Pezzo surfaces with two or three singular points.


title: A basis for a quotient of symmetric polynomials

abstract: Consider the ring $\mathcal{S}$ of symmetric polynomials in $k$ variables over an arbitrary base ring $\mathbf{k}$. Fix $k$ scalars $a_{1},a_{2},\ldots,a_{k}\in\mathbf{k}$. Let $I$ be the ideal of $\mathcal{S}$ generated by $h_{n-k+1}-a_{1},h_{n-k+2}-a_{2},\ldots,h_{n}-a_{k}$, where $h_{i}$ is the $i$-th complete homogeneous symmetric polynomial. The quotient ring $\mathcal{S}/I$ generalizes both the usual and the quantum cohomology of the Grassmannian. We show that $\mathcal{S}/I$ has a $\mathbf{k}$-module basis consisting of (residue classes of) Schur polynomials fitting into an $\left( n-k\right) \times k$-rectangle; and that its multiplicative structure constants satisfy the same $S_3$-symmetry as those of the Grassmannian cohomology. We prove a Pieri rule and a "rim hook algorithm", and conjecture a positivity property generalizing that of Gromov-Witten invariants. We construct two further bases of $\mathcal{S}/I$ as well. We also study the quotient of the whole polynomial ring (not just the symmetric polynomials) by the ideal generated by the same $k$ polynomials as $I$.


title: Finite-time stabilization control of quantum systems

abstract: The finite-time control problem of quantum systems is investigated in this paper. We first define finite-time stability and present a finite-time Lyapunov stability criterion for finite-dimensional quantum systems in coherence vector representation. Then, for two-level quantum systems, we design a continuous non-smooth control law with a state-dependent fractional power and prove the uniqueness of solutions of the system dynamics with the controller via the concept of transversality. By combining the finite-time Lyapunov stability criterion with the homogeneity theory, the finite-time convergence of the system to an eigenstate of its internal Hamiltonian is proved. Numerical results on a spin-1/2 system demonstrate the effectiveness of the proposed finite-time stabilization control scheme.


title: Orders generated by character values

abstract: Let $K:=\mathbb{Q}(G)$ be the number field generated by the complex character values of a finite group $G$. Let $\mathbb{Z}_K$ be the ring of integers of $K$. In this paper we investigate the suborder $\mathbb{Z}[G]$ of $\mathbb{Z}_K$ generated by the character values of $G$. We prove that every prime divisor of the order of the finite abelian group $\mathbb{Z}_K/\mathbb{Z}[G]$ divides $|G|$. Moreover, if $G$ is nilpotent, we show that the exponent of $\mathbb{Z}_K/\mathbb{Z}[G]$ is a proper divisor of $|G|$ unless $G=1$. We conjecture that this holds for arbitrary finite groups $G$.


title: Enormous enhancement of p-orbital magnetism and band gap in the lightly doped carbyne

abstract: This paper presents a path to tailor adapted magnetic and optical properties in carbyne. Although p-orbital magnetism is generally much weaker than d-orbital magnetism, we demonstrate that the charge fluctuation of the free radical electrons triggered by a time-varying electric dipole moment leads to enormous p-orbital magnetism. By introducing 25% arsenic and 12.5% fluorine into the monoatomic carbon chain, the magnetic moment of the arsenic atom reaches 2.9 Bohr Magneton, which is ~1.3 times stronger than magnetic moment of bulk Fe. This magnetically optimized carbyne composite carries an exchange-correlation energy of 22meV (~270K). On the other hand, we convert the carbyne (in beta-form) from metallic to a semiconducting state by using anionic dopants. After doping 12.5% nitrogen and 12.5% oxygen into the beta-carbyne, the semiconducting gap of this composite is optimized at 1.6eV, which is 1.4 times larger than the band gap of bulk silicon.


title: Reinforcement Learning for Multi-Objective Optimization of Online Decisions in High-Dimensional Systems

abstract: This paper describes a purely data-driven solution to a class of sequential decision-making problems with a large number of concurrent online decisions, with applications to computing systems and operations research. We assume that while the micro-level behaviour of the system can be broadly captured by analytical expressions or simulation, the macro-level or emergent behaviour is complicated by non-linearity, constraints, and stochasticity. If we represent the set of concurrent decisions to be computed as a vector, each element of the vector is assumed to be a continuous variable, and the number of such elements is arbitrarily large and variable from one problem instance to another. We first formulate the decision-making problem as a canonical reinforcement learning (RL) problem, which can be solved using purely data-driven techniques. We modify a standard approach known as advantage actor critic (A2C) to ensure its suitability to the problem at hand, and compare its performance to that of baseline approaches on the specific instance of a multi-product inventory management task. The key modifications include a parallelised formulation of the decision-making task, and a training procedure that explicitly recognises the quantitative relationship between different decisions. We also present experimental results probing the learned policies, and their robustness to variations in the data.


title: Generalizations of forest fires with ignition at origin

abstract: We study generalizations of the Forest Fire model introduced in [van den Berg, J., and J\'arai, A. A. "On the asymptotic density in a one-dimensional self-organized critical forest-fire model". Comm. Math. Phys. 253 (2005)] and [Volkov, Stanislav. "Forest fires on $\mathbb{Z}_+$ with ignition only at 0". ALEA 6 (2009)] by allowing the rates at which the tree grow to depend on their location, introducing long-range burning, as well as continuous-space generalization of the model. We establish that in all the models in consideration the time required to reach site at distance $x$ from the origin is of order at most $(\log x)^{(\log 2)^{-1}+\delta}$ for any $\delta>0$.


title: Joint Estimation of the Non-parametric Transitivity and Preferential Attachment Functions in Scientific Co-authorship Networks

abstract: We propose a statistical method to estimate simultaneously the non-parametric transitivity and preferential attachment functions in a growing network, in contrast to conventional methods that either estimate each function in isolation or assume some functional form for them. Our model is shown to be a good fit to two real-world co-authorship networks and be able to bring to light intriguing details of the preferential attachment and transitivity phenomena that would be unavailable under traditional methods. We also introduce a method to quantify the amount of contributions of those phenomena in the growth process of a network based on the probabilistic dynamic process induced by the model formula. Applying this method, we found that transitivity dominated PA in both co-authorship networks. This suggests the importance of indirect relations in scientific creative processes. The proposed methods are implemented in the R package FoFaF.

categories: physics.soc-ph cs.SI physics.data-an stat.AP stat.CO


title: Automatic Throughput and Critical Path Analysis of x86 and ARM Assembly Kernels

abstract: Useful models of loop kernel runtimes on out-of-order architectures require an analysis of the in-core performance behavior of instructions and their dependencies. While an instruction throughput prediction sets a lower bound to the kernel runtime, the critical path defines an upper bound. Such predictions are an essential part of analytic (i.e., white-box) performance models like the Roofline and Execution-Cache-Memory (ECM) models. They enable a better understanding of the performance-relevant interactions between hardware architecture and loop code. The Open Source Architecture Code Analyzer (OSACA) is a static analysis tool for predicting the execution time of sequential loops. It previously supported only x86 (Intel and AMD) architectures and simple, optimistic full-throughput execution. We have heavily extended OSACA to support ARM instructions and critical path prediction including the detection of loop-carried dependencies, which turns it into a versatile cross-architecture modeling tool. We show runtime predictions for code on Intel Cascade Lake, AMD Zen, and Marvell ThunderX2 micro-architectures based on machine models from available documentation and semi-automatic benchmarking. The predictions are compared with actual measurements.


title: Noisy Guesses

abstract: We consider the problem of guessing a random, finite-alphabet, secret $n$-vector, where the guesses are transmitted via a noisy channel. We provide a single-letter formula for the best achievable exponential growth rate of the $\rho$--th moment of the number of guesses, as a function of $n$. This formula exhibits a fairly clear insight concerning the penalty due to the noise. We describe two different randomized schemes that achieve the optimal guessing exponent. One of them is fully universal in the sense of being independent of source (that governs the vector to be guessed), the channel (that corrupts the guesses), and the moment power $\rho$. Interestingly, it turns out that, in general, the optimal guessing exponent function exhibits a phase transition when it is examined either as a function of the channel parameters, or as a function of $\rho$: as long as the channel is not too distant (in a certain sense to be defined precisely) from the identity channel (i.e., the clean channel), or equivalently, as long $\rho$ is larger than a certain critical value, $\rho_{\mbox{\tiny c}}$, there is no penalty at all in the guessing exponent, compared to the case of noiseless guessing.


title: Revisiting Fine-tuning for Few-shot Learning

abstract: Few-shot learning is the process of learning novel classes using only a few examples and it remains a challenging task in machine learning. Many sophisticated few-shot learning algorithms have been proposed based on the notion that networks can easily overfit to novel examples if they are simply fine-tuned using only a few examples. In this study, we show that in the commonly used low-resolution mini-ImageNet dataset, the fine-tuning method achieves higher accuracy than common few-shot learning algorithms in the 1-shot task and nearly the same accuracy as that of the state-of-the-art algorithm in the 5-shot task. We then evaluate our method with more practical tasks, namely the high-resolution single-domain and cross-domain tasks. With both tasks, we show that our method achieves higher accuracy than common few-shot learning algorithms. We further analyze the experimental results and show that: 1) the retraining process can be stabilized by employing a low learning rate, 2) using adaptive gradient optimizers during fine-tuning can increase test accuracy, and 3) test accuracy can be improved by updating the entire network when a large domain-shift exists between base and novel classes.


title: On Randomized Approximation of Scattered Data

abstract: Let $A$ be a set and $V$ a real Hilbert space. Let $H$ be a real Hilbert space of functions $f:A\to V$ and assume $H$ is continuously embedded in the Banach space of bounded functions. For $i=1,\cdots,n$, let $(x_i,y_i)\in A\times V$ comprise our dataset. Let $0<q<1$ and $f^*\in H$ be the unique global minimizer of the functional \begin{equation*} u(f) = \frac{q}{2}\Vert f\Vert_{H}^{2} + \frac{1-q}{2n}\sum_{i=1}^{n}\Vert f(x_i)-y_i\Vert_{V}^{2}. \end{equation*} For $x\in A$ and $v\in V$ let $\Phi(x,v)\in H$ be the unique element such that $(\Phi(x,v),f)_{H}=(f(x),v)_{V}$ for all $f\in H$. In this paper we show that for each $k\in\mathbb{N}$, $k\geq 2$ one has a random function $F_{k}\in H$ with the structure \begin{equation*} F_{k} = \sum_{h=1}^{N_k} \Lambda_{k, h} \Phi(x_{I_h}, \mathcal{E}_{h}) \end{equation*} (where $0\leq N_k\leq k-1$ are Binomially distributed with success probability $1-q$, $\Lambda_{k, h}\in\mathbb{R}$ are random coefficients, $1\leq I_{h}\leq n$ are independent and uniformly distributed and $\mathcal{E}_{h}\in V$ are random vectors) such that asymptotically for large $k$ we have \begin{equation*} E\left[ \Vert F_{k}-f^*\Vert_{H}^{2} \right] = O(\frac{1}{k}). \end{equation*} Thus we achieve the Monte Carlo type error estimate with no metric or measurability structure on $A$, possibly infinite dimensional $V$ and the ingredients of approximating functions are just the Riesz representatives $\Phi(x,v)\in H$. We obtain this result by considering the stochastic gradient descent sequence in the Hilbert space $H$ to minimize the functional $u$.


title: Influence of chirp, jitter and relaxation oscillations on probabilistic properties of laser pulse interference

abstract: Interference of laser pulses is an essential ingredient of quantum randomness; therefore, probabilistic properties of laser pulses gains new relevance. Here, we consider in detail the combined influence of the three effects - chirp, jitter and relaxation oscillations - on the probability density function of the interference of pulses from a gain-switched semiconductor laser. We develop a rigorous model based on laser rate equations and demonstrate that only consideration of all these three effects together allows describing the interference statistics properly. We supplement our theoretical calculations with corresponding measurements at various pump currents. Experimental results demonstrate perfect agreement with predictions of the model and are well reproduced by Monte-Carlo simulations.


title: Quantum noise extraction from the interference of laser pulses in optical quantum random number generator

abstract: We propose a method for quantum noise extraction from the interference of laser pulses with random phase. Our technique is based on the calculation of a parameter, which we called the quantum reduction factor, and which allows determining the contributions of quantum and classical noises in the assumption that classical fluctuations exhibit Gaussian distribution. To the best of our knowledge, the concept of the quantum reduction factor is introduced for the first time. We use such an approach to implement the post-processing-free optical quantum random number generator with the random bit generation rate of 2 Gbps.


title: The Nash Equilibrium with Inertia in Population Games

abstract: In the traditional game-theoretic set up, where agents select actions and experience corresponding utilities, an equilibrium is a configuration where no agent can improve their utility by unilaterally switching to a different action. In this work, we introduce the novel notion of inertial Nash equilibrium to account for the fact that, in many practical situations, action changes do not come for free. Specifically, we consider a population game and introduce the coefficients $c_{ij}$ describing the cost an agent incurs by switching from action $i$ to action $j$. We define an inertial Nash equilibrium as a distribution over the action space where no agent benefits in moving to a different action, while taking into account the cost of this change. First, we show that the set of inertial Nash equilibria contains all the Nash equilibria, but is in general not convex. Second, we argue that classical algorithms for computing Nash equilibria cannot be used in the presence of switching costs. We then propose a natural better-response dynamics and prove its convergence to an inertial Nash equilibrium. We apply our results to predict the drivers' distribution of an on-demand ride-hailing platform.


title: Optimal two-qubit states for quantum teleportation vis-\`a-vis state properties

abstract: Quantum teleportation with a two-qubit state can be suitably characterized in terms of maximal fidelity and fidelity deviation, where the former is the maximal value of the average fidelity achievable within the standard protocol and local unitary operations and the latter is the standard deviation of fidelity over all input states. In this paper, we consider the problem of characterizing two-qubit states that are optimal for quantum teleportation for a given value of some state property. The optimal states are defined as those states that, for a given value of the state property under consideration, achieve the largest maximal fidelity and also exhibit zero fidelity deviation. We provide a complete characterization of optimal states for a given linear entropy, maximum mean value of the Bell-CHSH observable, and concurrence, respectively. We find that for a given linear entropy or Bell-CHSH violation, the largest maximal fidelity states are optimal, but for a given concurrence, the optimal states form a strict subset of the largest maximal fidelity states.


title: The Role of NewSpace in Furthering Canadian Astronomy

abstract: [Highly abridged, from executive summary] As much as NewSpace presents opportunities, there are significant challenges that must be overcome, requiring engagement with policy makers to influence domestic and international space governance. Failure to do so could result in a range of long-lasting negative outcomes for science and space stewardship. How will the Canadian astronomical community engage with NewSpace? What are the implications for NewSpace on the astro-environment, including Earth orbits, lunar and cis-lunar orbits, and surfaces of celestial bodies? This white paper analyzes the rapid changes in space use and what those changes could mean for Canadian astronomers. Our recommendations are as follows: Greater cooperation between the astronomical and the Space Situational Awareness communities is needed. Build closer ties between the astronomical community and Global Affairs Canada (GAC). Establish a committee for evaluating the astro-environmental impacts of human space use, including on and around the Moon and other bodies. CASCA and the Tri-Council should coordinate to identify programs that would enable Canadian astronomers to participate in pay-for-use services at appropriate funding levels. CASCA should continue to foster a relationship with CSA, but also build close ties to the private space industry. Canadian-led deep space missions are within Canada's capabilities, and should be pursued.


title: An improved analysis and unified perspective on deterministic and randomized low rank matrix approximations

abstract: We introduce a Generalized LU-Factorization (\textbf{GLU}) for low-rank matrix approximation. We relate this to past approaches and extensively analyze its approximation properties. The established deterministic guarantees are combined with sketching ensembles satisfying Johnson-Lindenstrauss properties to present complete bounds. Particularly good performance is shown for the sub-sampled randomized Hadamard transform (SRHT) ensemble. Moreover, the factorization is shown to unify and generalize many past algorithms. It also helps to explain the effect of sketching on the growth factor during Gaussian Elimination.


title: Atoms in separated resonators can jointly absorb a single photon

abstract: The coherent nonlinear process where a single photon simultaneously excites two or more two-level systems (qubits) in a single-mode resonator has recently been theoretically predicted. Here we explore the case where the two qubits are placed in different resonators in an array of two or three weakly coupled resonators. Investigating different setups and excitation schemes, we show that this process can still occur with a probability approaching one under specific conditions. The obtained results provide interesting insights into subtle causality issues underlying the simultaneous excitation processes of qubits placed in different resonators.


title: Mistakes in Games

abstract: We define a new concept of "mistake" strategies and actions for strategic-form and extensive-form games, analyze the relationship to prior main game-theoretic solution concepts, study algorithms for computation, and explore practicality. This concept has potential applications to cybersecurity, for example detecting whether a human player is illegally using real-time assistance in games like poker.


title: Proper time evolution of magnetic susceptibility in a magnetized quark-gluon plasma

abstract: In ultrarelativistic heavy-ion collisions, enormous magnetic fields are generated because of fast-moving charged particles. In the presence of these magnetic fields, the spin of particles is aligned either in the parallel or in the antiparallel direction with respect to the direction of the magnetic field. A finite magnetization is thus produced. It is known that a finite magnetic susceptibility, $\chi_{m}$, changes the evolution of the energy density of the quark-gluon plasma (QGP), which is believed to be created in these collisions. Depending on whether the system under consideration is a paramagnetic ($\chi_{m}>0$) or diamagnetic ($\chi_{m}<0$) fluid, it slows down or speeds up the decay of the energy density, and affects other thermodynamic quantities. In general, one expects that the magnetic susceptibility depends on the magnetic field and temperature. Bearing in mind that these parameters evolve with the evolution of the fluid, a nonuniform magnetic susceptibility in this system is thus expected. In this work, we first determine $\chi_{m}$ by using a certain analogy to the standard anisotropic kinetic theory, where the one-particle distribution function is replaced by the corresponding anisotropic distribution function. We then determine the proper time dependence of the magnetic susceptibility in the framework of the ideal transverse magnetohydrodynamics. We also study the effect of dissipation on the evolution of $\chi_{m}$.


title: Small and Moderate Aperture Telescopes for Research and Education

abstract: In this white paper (WP), we highlight several examples of small and moderate aperture telescopes that are being used for education and/or research. We further discuss potential costs for establishing new, small observatories, as well as joining existing international consortia. The WP includes a brief overview of select observing sites, with a discussion on how small telescopes at exceptional observing locations can be competitive, under certain circumstances, with larger and more expensive facilities located at poorer sites. In addition to research, these facilities enable many different types of educational experiences for wide range of people, from high school students to undergraduates to graduate students to postdocs. Canada should remain committed to partnering with large, international observatories such as CFHT, Gemini, and TMT, but it should also negotiate international agreements and commit funding to expand the use of small and moderate research observatories at domestic and international sites through coordination with the NRC, the Tri-Council, and the Canadian Foundation for Innovation. Both capital and operational costs (with site rental costs allowed) need to be included in support possibilities. CASCA should establish and maintain a small to moderate telescope expression of interest database that would help to facilitate Canadian institutions in organizing consortia, particularly for smaller institutions. The astronomical community should work with the NRC to make existing facilities more accessible to the astronomical community for research. This could involve, for example, automating the Plaskett and/or providing travel funds for supporting classical observing modes. Finally, a small to moderate aperture facility in the Arctic would be a world-class observatory and should be advanced over the next decade.


title: Regularity of the solution of the scalar Signorini problem in polygonal domains

abstract: The Signorini problem for the Laplace operator is considered in a general polygonal domain. It is proved that the coincidence set consists of a finite number of boundary parts plus isolated points. The regularity of the solution is described. In particular, we show that the leading singularity is in general $r_i^{\pi/(2\alpha_i)}$ at transition points of Signorini to Dirichlet or Neumann conditions but $r_i^{\pi/\alpha_i}$ at kinks of the Signorini boundary, with $\alpha_i$ being the internal angle of the domain at these critical points.


title: Confidence intervals for median absolute deviations

abstract: The median absolute deviation (MAD) is a robust measure of scale that is simple to implement and easy to interpret. Motivated by this, we introduce interval estimators of the MAD to make reliable inferences for dispersion for a single population and ratios and differences of MADs for comparing two populations. Our simulation results show that the coverage probabilities of the intervals are very close to the nominal coverage for a variety of distributions. We have used partial influence functions to investigate the robustness properties of the difference and ratios of independent MADs.


title: Exploring Self-Assembling Behaviors in a Swarm of Bio-micro-robots using Surrogate-Assisted MAP-Elites

abstract: Swarms of molecular robots are a promising approach to create specific shapes at the microscopic scale through self-assembly. However, controlling their behavior is a challenging problem as it involves complex non-linear dynamics and high experimental variability. Hand-crafting a molecular controller will often be time-consuming and give sub-optimal results. Optimization methods, like the bioNEAT algorithm, were previously employed to partially overcome these difficulties, but they still had to cope with deceptive high-dimensional search spaces and computationally expensive simulations. Here, we describe a novel approach to solve this problem by using MAP-Elites, an algorithm that searches for both high-performing and diverse solutions. We then apply it to a molecular robotic framework we recently introduced that allows sensing, signaling and self-assembly at the micro-scale and show that MAP-Elites outperforms previous approaches. Additionally, we propose a surrogate model of micro-robots physics and chemical reaction dynamics to reduce the computational costs of simulation. We show that the resulting methodology is capable of optimizing controllers with similar accuracy as when using only a full-fledged realistic model, with half the computational budget.


title: On the equality between the infimum obtained by solving various Plateau's problem

abstract: In this paper we will compare the Plateau's problem with \v{C}ech and singular homological boundary conditions, we also compare these with the size minimizing problem for integral currents with a given boundary. Finally we get the agreement on the infimum values for these Plateau's problem.


title: Partition Function of the Bose-Einstein Condensed Dark Matter and the Modified Gross-Pitaevskii Equation

abstract: Intending to describe the dark matter of dwarf galaxies, we concentrate on one model of the slowly rotating and gravitating Bose-Einstein condensate. For a deeper understanding of its properties, we calculate the partition function and compare the characteristics derived from it with the results based on the solution of Gross-Pitaevskii equation. In our approach, which uses the Green's functions of spatial evolution operators, we formulate in a unified way the boundary conditions, important for applying the Thomas-Fermi approximation. Taking this into account, we revise some of the results obtained earlier. We also derive the spatial particle distribution, similar to the model with rotation, by using the deformation of commutation relations for a macroscopic wave function and modifying the Gross-Pitaevskii equation. It is shown that such an approach leads to the entropy inhomogeneity and makes the distribution dependent on temperature.


title: Scalar model of flocking dynamics on complex social networks

abstract: We investigate the effects of long-range social interactions in flocking dynamics by studying the dynamics of a scalar model of collective motion embedded in a complex network representing a pattern of social interactions, as observed in several social species. In this scalar model we find a phenomenology analogous to that observed in the classic Vicsek model: In networks with low heterogeneity, a phase transition separates an ordered from a disordered phase. At high levels of heterogeneity, instead, the transition is suppressed and the system is always ordered. This observation is backed up analytically by the solution of a modified scalar model within an heterogeneous mean-field approximation. Our work extends the understanding of the effects of social interactions in flocking dynamics and opens the path to the analytical study of more complex topologies of social ties.


title: Probing the Electroweak Sphaleron with Gravitational Waves

abstract: We present the relation between the sphaleron energy and the gravitational wave signals from a first order electroweak phase transition. The crucial ingredient is the scaling law between the sphaleron energy at the temperature of the phase transition and that at zero temperature. We estimate the baryon number preservation criterion, and observe that for a sufficiently strong phase transition, it is possible to probe the electroweak sphaleron using measurements of future space-based gravitational wave detectors.


title: Rich dynamics and anticontrol of extinction in a prey-predator system

abstract: This paper reveals some new and rich dynamics of a two-dimensional prey-predator system and to anticontrol the extinction of one of the species. For a particular value of the bifurcation parameter, one of the system variable dynamics is going to extinct, while another remains chaotic. To prevent the extinction, a simple anticontrol algorithm is applied so that the system orbits can escape from the vanishing trap. As the bifurcation parameter increases, the system presents quasiperiodic, stable, chaotic and also hyperchaotic orbits. Some of the chaotic attractors are Kaplan-Yorke type, in the sense that the sum of its Lyapunov exponents is positive. Also, atypically for undriven discrete systems, it is numerically found that, for some small parameter ranges, the system seemingly presents strange nonchaotic attractors. It is shown both analytically and by numerical simulations that the original system and the anticontrolled system undergo several Neimark-Sacker bifurcations. Beside the classical numerical tools for analyzing chaotic systems, such as phase portraits, time series and power spectral density, the 0-1 test is used to differentiate regular attractors from chaotic attractors.


title: Magnetostriction in magnetic gels and elastomers as a function of the internal structure and particle distribution

abstract: Magnetic gels and elastomers are promising candidates to construct reversibly excitable soft actuators, triggered from outside by magnetic fields. These magnetic fields induce or alter the magnetic interactions between discrete rigid particles embedded in a soft elastic polymeric matrix, leading to overall deformations. It is a major challenge in theory to correctly predict from the discrete particle configuration the type of deformation resulting for a finite-sized system. Considering an elastic sphere, we here present such an approach. The method is in principle exact, at least within the framework of linear elasticity theory and for large enough interparticle distances. Different particle arrangements are considered. We find, for instance, that regular simple cubic configurations show elongation of the sphere along the magnetization if oriented along a face or space diagonal of the cubic unit cell. Contrariwise, with the magnetization along the edge of the cubic unit cell, they contract. The opposite is true in this geometry for body- and face-centered configurations. Remarkably, for the latter configurations but the magnetization along a face or space diagonal of the unit cell, contraction was observed to revert to expansion with decreasing Poisson ratio of the elastic material. Randomized configurations were considered as well. They show a tendency of elongating the sphere along the magnetization, which is more pronounced for compressible systems. Our results can be tested against actual experiments for spherical samples. Moreover, our approach shall support the search of optimal particle distributions for a maximized effect of actuation.


title: Sampling Unknown Decision Functions to Build Classifier Copies

abstract: Copies have been proposed as a viable alternative to endow machine learning models with properties and features that adapt them to changing needs. A fundamental step of the copying process is generating an unlabelled set of points to explore the decision behavior of the targeted classifier throughout the input space. In this article we propose two sampling strategies to produce such sets. We validate them in six well-known problems and compare them with two standard methods. We evaluate our proposals in terms of both their accuracy performance and their computational cost.


title: Polynomial functions on rings of dual numbers over residue class rings of the integers

abstract: The ring of dual numbers over a ring $R$ is $R[\alpha] = R[x]/(x^2)$, where $\alpha$ denotes $x+(x^2)$. For any finite commutative ring $R$, we characterize null polynomials and permutation polynomials on $R[\alpha]$ in terms of the functions induced by their coordinate polynomials ($f_1,f_2\in R[x]$, where $f=f_1+\alpha f_2$) and their formal derivatives on $R$. We derive explicit formulas for the number of polynomial functions and the number of polynomial permutations on $\mathbb{Z}_{p^n}[\alpha]$ for $n\le p$ ($p$ prime).


title: A note on the cycle of curves in a product of pairs

abstract: We prove that the cycle-valued logarithmic Gromov--Witten theory of a product of simple normal crossings pairs $X\times Y$ decomposes into a product of pieces coming from $X$ and $Y$, provided that the decomposition is considered over a blowup of the moduli space of curves. The result is established by applying the toroidal semistable reduction theorem to the stabilization morphism on the stack of maps to the Artin fans. Since products of smooth pairs are naturally simple normal crossings pairs, the result provides a direct avenue of access to questions in logarithmic Gromov--Witten theory via relative Gromov--Witten theory.


title: About the Bloch-Connelly-Henderson Theorem on the simplexwise linear homeomorphisms of a convex 2-disk

abstract: This paper gives an improved version of the original proof of Bloch-Connelly-Henderson's theorem about the space of $SL$ homeomorphisms of a convex 2-disk. A major improvement is related to the main lemma of the original paper.


title: Optimal Dyck Reachability for Data-Dependence and Alias Analysis

abstract: A fundamental algorithmic problem at the heart of static analysis is Dyck reachability. The input is a graph where the edges are labeled with different types of opening and closing parentheses, and the reachability information is computed via paths whose parentheses are properly matched. We present new results for Dyck reachability problems with applications to alias analysis and data-dependence analysis. Our main contributions, that include improved upper bounds as well as lower bounds that establish optimality guarantees, are as follows. First, we consider Dyck reachability on bidirected graphs, which is the standard way of performing field-sensitive points-to analysis. Given a bidirected graph with $n$ nodes and $m$ edges, we present: (i)~an algorithm with worst-case running time $O(m + n \cdot \alpha(n))$, where $\alpha(n)$ is the inverse Ackermann function, improving the previously known $O(n^2)$ time bound; (ii)~a matching lower bound that shows that our algorithm is optimal wrt to worst-case complexity; and (iii)~an optimal average-case upper bound of $O(m)$ time, improving the previously known $O(m \cdot \log n)$ bound. Second, we consider the problem of context-sensitive data-dependence analysis, where the task is to obtain analysis summaries of library code in the presence of callbacks. Our algorithm preprocesses libraries in almost linear time, after which the contribution of the library in the complexity of the client analysis is only linear, and only wrt the number of call sites. Third, we prove that combinatorial algorithms for Dyck reachability on general graphs with truly sub-cubic bounds cannot be obtained without obtaining sub-cubic combinatorial algorithms for Boolean Matrix Multiplication, which is a long-standing open problem. We also show that the same hardness holds for graphs of constant treewidth.


title: Geometric effects on the electronic structure and the bound states in annular corrugated wires

abstract: In the spirit of the thin-layer quantization scheme, we give the effective Hamiltonian describing the noninteracting electrons confined to an annular corrugated surface, and find that the geometrically induced potential is considerably influenced by corrugations. By using numerical calculation, we investigate the eigenenergies and the corresponding eigenstates, and find that the transition energies can be sufficiently improved by adding corrugations. Particularly, the transition energy between the adjacent eigenstates corresponds to energy levels difference based on the wavefunction of annular wire, and the number of the energy levels is equal to the number of corrugations. And the larger magnitude of corrugations is capable of increasing the number of bound states. In addition, the distribution of ground state probability density is reconstructed by the corrugations, and the energy shift is generated.


title: Extension of Lipschitz-type operators on Banach function spaces

abstract: We study extension theorems for Lipschitz-type operators acting on metric spaces and with values on spaces of integrable functions. Pointwise domination is not a natural feature of such spaces, and so almost everywhere inequalities and other measure-theoretic notions are introduced.%more appropriate. Thus, we adapt the classical definition of Lipschitz map to the context of spaces of integrable functions by introducing such elements. We analyze Lipschitz type inequalities in two fundamental cases. The first concerns a.e. pointwise inequalities, while the second considers dominations involving integrals. These Lipschitz type inequalities provide the suitable frame to work with operators that take values on Banach function spaces. In the last part of the paper we use some interpolation procedures to extend our study to interpolated Banach function spaces.


title: Impartial SWIPT-Assisted User Cooperation Schemes

abstract: In this paper, we propose an impartial simultaneous wireless information and power transfer (SWIPT)-assisted cooperation mechanism for a non-orthogonal multiple access (NOMA) downlink scenario. Specifically, both a cell-center user and a cell-edge user apply the power-splitting technique and utilize the harvested energy to forward the other user's information on the premise of successful decoding of their own information. Both analytical and numerical results show that the proposed impartial user cooperation mechanism (IUCM) outperforms the traditional partial cooperation mechanism in terms of outage probability, diversity order and diversity-multiplexing trade-off (DMT). For comparison, we further incorporate the IUCM into an orthogonal frequency-division multiple access (OFDMA) framework, which is shown to preserve the same diversity order, while has a worse but more flexible DMT performance in comparison with the IUCM in the NOMA framework. Although the IUCM in OFDMA has a worse outage performance, it is proved that it has the same optimal system outage probability with the IUCM in NOMA when the relaying channel between the two users is error-free.


title: Long-distance dissipation-assisted transport of entangled states via a chiral waveguide

abstract: Quantum networks provide a prominent platform for realizing quantum information processing and quantum communication, with entanglement being a key resource in such applications. Here, we describe the dissipative transport protocol for entangled states, where entanglement stored in the first node of quantum network can be transported with high fidelity to the second node via a 1D chiral waveguide. In particular, we exploit the directional asymmetry in chirally-coupled single-mode ring resonators to transport entangled states. For the fully chiral waveguide, Bell states, multipartite $W$-states and and Dicke states can be transported with fidelity as high as $0.954$, despite the fact that the communication channel is noisy. Our proposal can be utilized for long-distance distribution of multipartite entangled states between the quantum nodes of the open quantum network.


title: MTab: Matching Tabular Data to Knowledge Graph using Probability Models

abstract: This paper presents the design of our system, namely MTab, for Semantic Web Challenge on Tabular Data to Knowledge Graph Matching (SemTab 2019). MTab combines the voting algorithm and the probability models to solve critical problems of the matching tasks. Results on SemTab 2019 show that MTab obtains promising performance for the three matching tasks.


title: Formation of self-propelling clusters starting from randomly dispersed Brownian particles

abstract: We present a simple chemical strategy for the formation of a self-propelling cluster via the process of capture and assembly of passive colloids on the surface of a chemically active colloid. The two species of colloids that are isotropic and Brownian otherwise interact to form propelling cluster. With the help of coarse-grained numerical simulations, we show that a chemically active colloid can induce diffusiophoretic motility to nearby chemically inert colloids towards itself. This propulsion and then self-assembly can then lead to the formation of active cluster. We observe the formation of propelling dimers, trimers, tetramers, etc. depending on the chemical activity and size of the colloids.


title: Decoy-state round-robin differential-phase-shift quantum key distribution with source errors

abstract: As a promising quantum key distribution (QKD), most of the existing round-robin differential-phase-shift quantum key distribution (RRDPS-QKD) protocols have adopted the decoy-state method and have assumed the source states are exactly controlled. However, the precise manipulation of source states is impossible for any practical experiment, and the RRDPS-QKD with source errors has an unignorable impact on the performance of the protocol. In the paper, we study the three-intensity decoy-state RRDPS-QKD protocol with source errors, formulate the secure generation key rate of the proposed protocol, do the numerical simulations to testify the deductions. The results show that our evaluation can estimate the influence of source errors.


title: Population inversion and entanglement in single and double glassy Jaynes-Cummings models

abstract: We find that a suppression of the collapse and revival of population inversion occurs in response to insertion of Gaussian quenched disorder in atom-cavity interaction strength in the Jaynes-Cummings model. The character of suppression can be significantly different in the presence of non-Gaussian disorder, which we uncover by studying the cases when the disorder is uniform, discrete, and Cauchy-Lorentz. Interestingly, the quenched averaged atom-photon entanglement keeps displaying nontrivial oscillations even after the population inversion has been suppressed. Subsequently, we show that disorder in atom-cavity interactions helps to avoid sudden death of atom-atom entanglement in the double Jaynes-Cummings model. We identify the minimal disorder strengths required to eliminate the possibility of sudden death. We also investigate the response of entanglement sudden death in the disordered double Jaynes-Cummings model in the presence of atom-atom coupling.


title: Reduction of multiple scattering of high-energy positively charged particles during channeling in single crystals

abstract: We present the experimental observation of the reduction of multiple scattering of high-energy positively charged particles during channeling in single crystals. According to our measurements the rms angle of multiple scattering in the plane orthogonal to the plane of the channeling is less than half that for non-channeled particles moving in the same crystal. In the experiment we use focusing bent single crystals. Such crystals have a variable thickness in the direction of beam propagation. This allows us to measure rms angles of scattering as a function of thickness for channeled and non-channeled particles. The behaviour with thickness of non-channeled particles is in agreement with expectations whereas the behaviour of channeled particles has unexpected features. We give a semi-quantitative explanation of the observed effect.


title: Local geometry of Jordan classes in semisimple algebraic groups

abstract: We prove that the closure of every Jordan class J in a semisimple simply connected complex group G at a point x with Jordan decomposition x = rv is smoothly equivalent to the union of closures of those Jordan classes in the centraliser of r that are contained in J and contain x in their closure. For x unipotent we also show that the closure of J around x is smoothly equivalent to the closure of a Jordan class in Lie(G) around exp^{-1}x. For G simple we apply these results in order to determine a (non-exhaustive) list of smooth sheets in G, the complete list of regular Jordan classes whose closure is normal and Cohen-Macaulay, and to prove that all sheets and Lusztig's strata in SL(n,C) are smooth.


title: A coarse-grained deep neural network model for liquid water

abstract: We introduce a coarse-grained deep neural network model (CG-DNN) for liquid water that utilizes 50 rotational and translational invariant coordinates, and is trained exclusively against energies of ~30,000 bulk water configurations. Our CG-DNN potential accurately predicts both the energies and molecular forces of water; within 0.9 meV/molecule and 54 meV/angstrom of a reference (coarse-grained bond-order potential) model. The CG-DNN water model also provides good prediction of several structural, thermodynamic, and temperature dependent properties of liquid water, with values close to that obtained from the reference model. More importantly, CG-DNN captures the well-known density anomaly of liquid water observed in experiments. Our work lays the groundwork for a scheme where existing empirical water models can be utilized to develop fully flexible neural network framework that can subsequently be trained against sparse data from high-fidelity albeit expensive beyond-DFT calculations.


title: Toward canonical convex functions in Alexandrov spaces

abstract: We construct for every finite-dimensional Alexandrov space $A$ and every point $p \in A$ a $2$-convex function $f_p$ in a small neighborhood around $p$, which approximates $\operatorname{dist}_p^2$ up to second order. Moreover, the function $f_p$ can be lifted to Gromov-Hausdorff close Alexandrov spaces of the same dimension.


title: Multilingual End-to-End Speech Translation

abstract: In this paper, we propose a simple yet effective framework for multilingual end-to-end speech translation (ST), in which speech utterances in source languages are directly translated to the desired target languages with a universal sequence-to-sequence architecture. While multilingual models have shown to be useful for automatic speech recognition (ASR) and machine translation (MT), this is the first time they are applied to the end-to-end ST problem. We show the effectiveness of multilingual end-to-end ST in two scenarios: one-to-many and many-to-many translations with publicly available data. We experimentally confirm that multilingual end-to-end ST models significantly outperform bilingual ones in both scenarios. The generalization of multilingual training is also evaluated in a transfer learning scenario to a very low-resource language pair. All of our codes and the database are publicly available to encourage further research in this emergent multilingual ST topic.


title: Approach to equilibrium via Tsallis distributions in a realistic ionic--crystal model and in the FPU model

abstract: In Statistical Mechanics, Tsallis distributions were apparently conceived in connection with systems presenting long--range interactions. In fact, they were observed in numerical computations for models of such a type, as occurring in the approach to equilibrium, i.e., to a Maxwell--Boltzmann distribution. Here we exhibit two apparently new results. The first one is that Tsallis distributions occur also in an ionic--crystal model with long--range Coulomb forces, which is so realistic as to reproduce in an impressively good way the experimental infrared spectra. Thus such distributions may be expected to be actual physical features of crystals. The second result is that Tsallis distributions occur in the standard short--range FPU model too, so that the presence of long--range interactions is not a necessary condition for Tsallis distributions to occur. In fact, this is in agreement with a previous result of the first author in connection with the statistics of return times for the classical FPU model. We thus confirm the thesis advanced by Tsallis himself, that the relevant property for a dynamical system to present Tsallis distributions is that its dynamics should be not fully chaotic, a property which is known to actually pertain to long--range systems.


title: A review of problem- and team-based methods for teaching statistics in Higher Education

abstract: The teaching of statistics in higher education in the UK is still largely lecture-based. This is despite recommendations such as those given by the American Statistical Association's GAISE report that more emphasis should be placed on active learning strategies where students take more responsibility for their own learning. One possible model is that of collaborative learning, where students learn in groups through carefully crafted `problems', which has long been suggested as a strategy for teaching statistics. In this article, we review two specific approaches that fall under the collaborative learning model: problem- and team-based learning. We consider the evidence for changing to this model of teaching in statistics, as well as give practical suggestions on how this could be implemented in typical statistics classes in Higher Education.