title: Branch Prediction Is Not a Solved Problem: Measurements, Opportunities, and Future Directions

abstract: Modern branch predictors predict the vast majority of conditional branch instructions with near-perfect accuracy, allowing superscalar, out-of-order processors to maximize speculative efficiency and thus performance. However, this impressive overall effectiveness belies a substantial missed opportunity in single-threaded instructions per cycle (IPC). For example, we show that correcting the mispredictions made by the state-of-the-art TAGE-SC-L branch predictor on SPECint 2017 would improve IPC by margins similar to an advance in process technology node. In this work, we measure and characterize these mispredictions. We find that they categorically arise from either (1) a small number of systematically hard-to-predict (H2P) branches; or (2) rare branches with low dynamic execution counts. Using data from SPECint 2017 and additional large code footprint applications, we quantify the occurrence and IPC impact of these two categories. We then demonstrate that increasing the resources afforded to existing branch predictors does not alone address the root causes of most mispredictions. This leads us to reexamine basic assumptions in branch prediction and to propose new research directions that, for example, deploy machine learning to improve pattern matching for H2Ps, and use on-chip phase learning to track long-term statistics for rare branches.


title: Effectiveness of Data Augmentation in Cellular-based Localization Using Deep Learning

abstract: Recently, deep learning-based positioning systems have gained attention due to their higher performance relative to traditional methods. However, obtaining the expected performance of deep learning-based systems requires large amounts of data to train model. Obtaining this data is usually a tedious process which hinders the utilization of such deep learning approaches. In this paper, we introduce a number of techniques for addressing the data collection problem for deep learning-based cellular localization systems. The basic idea is to generate synthetic data that reflects the typical pattern of the wireless data as observed from a small collected dataset. Evaluation of the proposed data augmentation techniques using different Android phones in a cellular localization case study shows that we can enhance the performance of the localization systems in both indoor and outdoor scenarios by 157% and 50.5%, respectively. This highlights the promise of the proposed techniques for enabling deep learning-based localization systems.


title: MediaPipe: A Framework for Building Perception Pipelines

abstract: Building applications that perceive the world around them is challenging. A developer needs to (a) select and develop corresponding machine learning algorithms and models, (b) build a series of prototypes and demos, (c) balance resource consumption against the quality of the solutions, and finally (d) identify and mitigate problematic cases. The MediaPipe framework addresses all of these challenges. A developer can use MediaPipe to build prototypes by combining existing perception components, to advance them to polished cross-platform applications and measure system performance and resource consumption on target platforms. We show that these features enable a developer to focus on the algorithm or model development and use MediaPipe as an environment for iteratively improving their application with results reproducible across different devices and platforms. MediaPipe will be open-sourced at https://github.com/google/mediapipe.


title: Write-Optimized and Consistent RDMA-based NVM Systems

abstract: In order to deliver high performance in cloud computing, we generally exploit and leverage RDMA (Remote Direct Memory Access) in networking and NVM (Non-Volatile Memory) in end systems. Due to no involvement of CPU, one-sided RDMA becomes efficient to access the remote memory, and NVM technologies have the strengths of non-volatility, byte-addressability and DRAM-like latency. In order to achieve end-to-end high performance, many efforts aim to synergize one-sided RDMA and NVM. Due to the need to guarantee Remote Data Atomicity (RDA), we have to consume extra network round-trips, remote CPU participation and double NVM writes. In order to address these problems, we propose a zero-copy log-structured memory design for Efficient Remote Data Atomicity, called Erda. In Erda, clients directly transfer data to the destination address at servers via one-sided RDMA writes without redundant copy and remote CPU consumption. To detect the incompleteness of fetched data, we verify a checksum without client-server coordination. We further ensure metadata consistency by leveraging an 8-byte atomic update in the hash table, which also contains the address information for the stale data. When a failure occurs, the server properly restores to a consistent version. Experimental results show that compared with Redo Logging (a CPU involvement scheme) and Read After Write (a network dominant scheme), Erda reduces NVM writes approximately by 50%, as well as significantly improves throughput and decreases latency.


title: Optical follow-up study of 32 high-redshift galaxy cluster candidates from Planck with the William Herschel Telescope

abstract: The Planck satellite has detected cluster candidates via the Sunyaev Zel'dovich (SZ) effect, but the optical follow-up required to confirm these candidates is still incomplete, especially at high redshifts and for SZ detections at low significance. In this work we present our analysis of optical observations obtained for 32 Planck cluster candidates using ACAM on the 4.2-m William Herschel Telescope. These cluster candidates were preselected using SDSS, WISE, and Pan-STARRS images to likely represent distant clusters at redshifts $z \gtrsim 0.7$. We obtain photometric redshift and richness estimates for all of the cluster candidates from a red-sequence analysis of $r$-, $i$-, and $z$-band imaging data. In addition, long-slit observations allow us to measure the redshifts of a subset of the clusters spectroscopically. The optical richness is often lower than expected from the inferred SZ mass when compared to scaling relations previously calibrated at low redshifts. This likely indicates the impact of Eddington bias and projection effects or noise-induced detections, especially at low SZ-significance. Thus, optical follow-up not only provides redshift measurements, but also an important independent verification method. We find that 18 (7) of the candidates at redshifts $z > 0.5$ ($z > 0.8$) are at least half as rich as expected from scaling relations, thereby clearly confirming these candidates as massive clusters. While the complex selection function of our sample due to our preselection hampers its use for cosmological studies, we do provide a validation of massive high-redshift clusters particularly suitable for further astrophysical investigations.


title: Identities in Brandt semigroups, revisited

abstract: We present a new proof for the main claim made in the author's paper "On the identity bases of Brandt semigroups" (Ural. Gos. Univ. Mat. Zap. 14, no.1 (1985), 38--42); this claim provides an identity basis for an arbitrary Brandt semigroup over a group of finite exponent. We also show how to fill a gap in the original proof of the claim in loc. cit.


title: MaGPoS -- A novel decentralized consensus mechanism combining magnetism and proof of stake

abstract: We describe MaGPoS, a novel consensus mechanism which is well suited to decentralized blockchain based protocols. MaGPoS is based on a combination of the well known physics of nano-scale magnetism, and previous implementations of proof of stake. This system has been studied by hundreds of thousands of scientists worldwide for over a hundred years, giving it an extreme level of reliability that is needed for a consensus mechanism. We start by explaining the physics, and study the properties that make it particularly beneficial for use in a consensus mechanism. We then show how to apply the physical model to a decentralized network of nodes, each with their own copy of a blockchain. After this, we describe some example calculations that a node in the decentralized network would make, and provide pseudo code for implementation. Finally, we discuss the how the model achieves all of the important properties that one expects of a consensus mechanism.


title: AI-enabled Blockchain: An Outlier-aware Consensus Protocol for Blockchain-based IoT Networks

abstract: A new framework for a secure and robust consensus in blockchain-based IoT networks is proposed using machine learning. Hyperledger fabric, which is a blockchain platform developed as part of the Hyperledger project, though looks very apt for IoT applications, has comparatively low tolerance for malicious activities in an untrustworthy environment. To that end, we propose AI-enabled blockchain (AIBC) with a 2-step consensus protocol that uses an outlier detection algorithm for consensus in an IoT network implemented on hyperledger fabric platform. The outlier-aware consensus protocol exploits a supervised machine learning algorithm which detects anomaly activities via a learned detector in the first step. Then, the data goes through the inherent Practical Byzantine Fault Tolerance (PBFT) consensus protocol in the hyperledger fabric for ledger update. We measure and report the performance of our framework with respect to the various delay components. Results reveal that our implemented AIBC network (2-step consensus protocol) improves hyperledger fabric performance in terms of fault tolerance by marginally compromising the delay performance.


title: Strategy Representation by Decision Trees with Linear Classifiers

abstract: Graph games and Markov decision processes (MDPs) are standard models in reactive synthesis and verification of probabilistic systems with nondeterminism. The class of $\omega$-regular winning conditions; e.g., safety, reachability, liveness, parity conditions; provides a robust and expressive specification formalism for properties that arise in analysis of reactive systems. The resolutions of nondeterminism in games and MDPs are represented as strategies, and we consider succinct representation of such strategies. The decision-tree data structure from machine learning retains the flavor of decisions of strategies and allows entropy-based minimization to obtain succinct trees. However, in contrast to traditional machine-learning problems where small errors are allowed, for winning strategies in graph games and MDPs no error is allowed, and the decision tree must represent the entire strategy. In this work we propose decision trees with linear classifiers for representation of strategies in graph games and MDPs. We have implemented strategy representation using this data structure and we present experimental results for problems on graph games and MDPs, which show that this new data structure presents a much more efficient strategy representation as compared to standard decision trees.


title: Equivariant higher twisted K-theory of SU(n) for exponential functor twists

abstract: We prove that each exponential functor on the category of finite-dimensional complex inner product spaces and isomorphisms gives rise to an equivariant higher (ie. non-classical) twist of $K$-theory over $G=SU(n)$. This twist is represented by a Fell bundle $\mathcal{E} \to \mathcal{G}$, which reduces to the basic gerbe for the top exterior power functor. The groupoid $\mathcal{G}$ comes equipped with a $G$-action and an augmentation map $\mathcal{G} \to G$, that is an equivariant equivalence. The $C^*$-algebra $C^*(\mathcal{E})$ associated to $\mathcal{E}$ is stably isomorphic to the section algebra of a locally trivial bundle with stabilised strongly self-absorbing fibres. Using a version of the Mayer-Vietoris spectral sequence we compute the equivariant higher twisted $K$-groups $K^G_*(C^*(\mathcal{E}))$ for arbitrary exponential functor twists over $SU(2)$, and also over $SU(3)$ after rationalisation.


title: Standalone and RTK GNSS on 30,000 km of North American Highways

abstract: There is a growing need for vehicle positioning information to support Advanced Driver Assistance Systems (ADAS), Connectivity (V2X), and Automated Driving (AD) features. These range from a need for road determination (<5 meters), lane determination (<1.5 meters), and determining where the vehicle is within the lane (<0.3 meters). This work examines the performance of Global Navigation Satellite Systems (GNSS) on 30,000 km of North American highways to better understand the automotive positioning needs it meets today and what might be possible in the near future with wide area GNSS correction services and multi-frequency receivers. This includes data from a representative automotive production GNSS used primarily for turn-by-turn navigation as well as an Inertial Navigation System which couples two survey grade GNSS receivers with a tactical grade Inertial Measurement Unit (IMU) to act as ground truth. The latter utilized networked Real-Time Kinematic (RTK) GNSS corrections delivered over a cellular modem in real-time. We assess on-road GNSS accuracy, availability, and continuity. Availability and continuity are broken down in terms of satellite visibility, satellite geometry, position type (RTK fixed, RTK float, or standard positioning), and RTK correction latency over the network. Results show that current automotive solutions are best suited to meet road determination requirements at 98% availability but are less suitable for lane determination at 57%. Multi-frequency receivers with RTK corrections were found more capable with road determination at 99.5%, lane determination at 98%, and highway-level lane departure protection at 91%.


title: Engineering stable quantum currents at bulk boundaries

abstract: We study transport properties of discrete quantum dynamical systems on the lattice, in particular Coined Quantum Walks and the Chalker--Coddington model. We prove existence of a non trivial charge transport and that the absolutely continuous spectrum covers the whole unit circle under mild assumptions. For Quantum Walks we exhibit explicit constructions of coins which imply existence of stable directed quantum currents along classical curves. The results are of topological nature and independent of the details of the model.


title: Observing the Effect of Polarization Mode Dispersion on Nonlinear Interference Generation in Wide-Band Optical Links

abstract: With the extension of the spectral exploitation of optical fibers beyond the C-band, accurate modeling and simulation of nonlinear interference (NLI) generation is of the utmost performance. Models and numerical simulation tools rely on the widely used Manakov equation (ME): however, this approach when considering also the effect of polarization mode dispersion (PMD) is formally valid only over a narrow optical bandwidth. In order to analyze the range of validity of the ME and its applicability to future wide-band systems, we present numerical simulations, showing the interplay between NLI generation and PMD over long dispersion-uncompensated optical links, using coherent polarization division multiplexing (PDM) quadrature amplitude modulation (QAM) formats. Using a Monte-Carlo analysis of different PMD realizations based on the coupled nonlinear Schr\"{o}dinger equations, we show that PMD has a negligible effect on NLI generation, independently from the total system bandwidth. Based on this, we give strong numerical evidence that the ME can be safely used to estimate NLI generation well beyond its bandwidth of validity that is limited to the PMD coherence bandwidth.


title: Some remarks on Petty projection of log-concave functions

abstract: In this note we study the Petty projection of a log-concave function, which has been recently introduced in [9]. Moreover, we present some new inequalities involving this new notion, partly complementing and correcting some results from [9].


title: A planarity estimate for pinched solutions of mean curvature flow

abstract: We show that the blow-ups of compact solutions to the mean curvature flow in $\mathbb{R}^N$ initially satisfying the pinching condition $|A|^2 < c |H|^2$ for a suitable constant $c = c(n)$ must be codimension one.


title: Quantum Cosmological Backreactions II: Purely Homogeneous Quantum Cosmology

abstract: This is the second paper in a series of four in which we use space adiabatic methods in order to incorporate backreactions among the homogeneous and between the homogeneous and inhomogeneous degrees of freedom in quantum cosmological perturbation theory. The purpose of the present paper is twofold. On the one hand, it illustrates the formalism of space adiabatic perturbation theory (SAPT) for two simple quantum mechanical toy models. On the other, it proves the main point, namely that backreactions lead to additional correction terms in effective Hamiltonians that one would otherwise neglect in a crude Born-Oppenheimer approximation. The first model that we consider is a harmonic oscillator coupled to an anharmonic oscillator. We chose it because it displays many similarities with the more interesting second model describing the coupling between an inflaton and gravity restricted to the purely homogeneous and isotropic sector. These results have potential phenomenological consequences in particular for quantum cosmological theories describing big bounces such as Loop Quantum Cosmology (LQC).


title: Modifying AMG coarse spaces with weak approximation property to exhibit approximation in energy norm

abstract: Algebraic multigrid (AMG) coarse spaces are commonly constructed so that they exhibit the so-called weak approximation (WAP) property which is necessary and sufficient condition for uniform two-grid convergence. This paper studies a modification of such coarse spaces so that the modified ones provide approximation in energy norm. Our modification is based on the projection in energy norm onto an orthogonal complement of original coarse space. This generally leads to dense modified coarse space matrices which is hence computationally infeasible. To remedy this, based on the fact that the projection involves inverse of a well-conditioned matrix, we use polynomials to approximate the projection and, therefore, obtain a practical, sparse modified coarse matrix and prove that the modified coarse space maintains computationally feasible approximation in energy norm. We present some numerical results for both, PDE discretization matrices as well as graph Laplacian ones, which are in accordance with our theoretical results.


title: The Fornax3D project: Tracing the assembly history of the cluster from the kinematic and line-strength maps

abstract: The 31 brightest galaxies (m_B < 15 mag) inside the virial radius of the Fornax cluster were observed from the centres to the outskirts with the Multi Unit Spectroscopic Explorer on the Very Large Telescope. These observations provide detailed high-resolution maps of the line-of-sight kinematics and line strengths of the stars and ionised gas reaching 2-3 Re for 21 early-type galaxies and 1-2 Re for 10 late-type galaxies. The majority of the galaxies are regular rotators, with eight hosting a kinematically distinct core. Only two galaxies are slow rotators. The mean age, total metallicity, and [Mg/Fe] abundance ratio in the bright central region inside 0.5 Re and in the galaxy outskirts are presented. Extended emission-line gas is detected in 13 galaxies, most of them are late-type objects with wide-spread star formation. The measured structural properties are analysed in relation to the galaxies' position in the projected phase space of the cluster. This shows that the Fornax cluster appears to consist of three main groups of galaxies inside the virial radius: the old core; a clump of galaxies, which is aligned with the local large-scale structure and was accreted soon after the formation of the core; and a group of galaxies that fell in more recently.


title: Critical metrics for Log-determinant functionals in conformal geometry

abstract: We consider critical points of a class of functionals on compact four-dimensional manifolds arising from Regularized Determinants for conformally covariant operators, whose explicit form was derived in [10], extending Polyakov's formula. These correspond to solutions of elliptic equations of Liouville type that are quasilinear, of mixed orders and of critical type. After studying existence, asymptotic behaviour and uniqueness of fundamental solutions, we prove a quantization property under blow-up, and then derive existence results via critical point theory.


title: Reward Prediction Error as an Exploration Objective in Deep RL

abstract: A major challenge in reinforcement learning is exploration, when local dithering methods such as epsilon-greedy sampling are insufficient to solve a given task. Many recent methods have proposed to intrinsically motivate an agent to seek novel states, driving the agent to discover improved reward. However, while state-novelty exploration methods are suitable for tasks where novel observations correlate well with improved reward, they may not explore more efficiently than epsilon-greedy approaches in environments where the two are not well-correlated. In this paper, we distinguish between exploration tasks in which seeking novel states aids in finding new reward, and those where it does not, such as goal-conditioned tasks and escaping local reward maxima. We propose a new exploration objective, maximizing the reward prediction error (RPE) of a value function trained to predict extrinsic reward. We then propose a deep reinforcement learning method, QXplore, which exploits the temporal difference error of a Q-function to solve hard exploration tasks in high-dimensional MDPs. We demonstrate the exploration behavior of QXplore on several OpenAI Gym MuJoCo tasks and Atari games and observe that QXplore is comparable to or better than a baseline state-novelty method in all cases, outperforming the baseline on tasks where state novelty is not well-correlated with improved reward.


title: Control What You Can: Intrinsically Motivated Task-Planning Agent

abstract: We present a novel intrinsically motivated agent that learns how to control the environment in the fastest possible manner by optimizing learning progress. It learns what can be controlled, how to allocate time and attention, and the relations between objects using surprise based motivation. The effectiveness of our method is demonstrated in a synthetic as well as a robotic manipulation environment yielding considerably improved performance and smaller sample complexity. In a nutshell, our work combines several task-level planning agent structures (backtracking search on task graph, probabilistic road-maps, allocation of search efforts) with intrinsic motivation to achieve learning from scratch.


title: Low-energy neutrino scattering in experiment and astrophysics

abstract: We review the relevance of neutrino-nucleus interactions at energy transfers below 100 MeV for accelerator-based experiments, experiments at lower energies and for astrophysical neutrinos. The impact of low-energy scattering processes in the energy reconstruction analysis of oscillation experiments is investigated. We discuss the modeling of coherent scattering processes and compare its strength to that of inelastic interactions. The presented results are obtained within a continuum random phase approximation approach.


title: Exploiting Weak Head Movements for Camera-based Respiration Detection

abstract: In recent years, considerable progress has been made in the non-contact based detection of the respiration rate from video sequences. Common techniques either directly assess the movement of the chest due to breathing or are based on analyzing subtle color changes that occur as a result of hemodynamic properties of the skin tissue by means of remote photoplethysmography (rPPG). However, extracting hemodynamic parameters from rPPG is often difficult especially if the skin is not visible to the camera. In contrast, extracting respiratory signals from chest movements turned out to be a robust method. However, the detectability of chest regions cannot be guaranteed in any application scenario, for instance if the camera setting is optimized to provide close-up images of the head. In such a case an alternative method for respiration detection is required. It is reasonable to assume that the mechanical coupling between chest and head induces minor movements of the head which, like in rPPG, can be detected from subtle color changes as well. Although the strength of these movements is expected to be much smaller in scale, sensing these intensity variations could provide a reasonably suited respiration signal for subsequent respiratory rate analysis. In order to investigate this coupling we conducted an experimental study with 12 subjects and applied motion- and rPPG-based methods to estimate the respiratory frequency from both head regions and chest. Our results show that it is possible to derive signals correlated to chest movement from facial regions. The method is a feasible alternative to rPPG-based respiratory rate estimations when rPPG-signals cannot be derived reliably and chest movement detection cannot be applied as well.


title: Generic quantized zero-bias conductance peaks in superconductor-semiconductor hybrid structures

abstract: We show theoretically that quantized zero-bias conductance peaks should be ubiquitous in superconductor-semiconductor hybrids by employing a zero-dimensional random matrix model with continuous tuning parameters. We demonstrate that a normal metal-superconductor (NS) junction conductance spectra can be generically obtained in this model replicating all features seen in recent experimental results. The theoretical quantized conductance peaks, which explicitly do not arise from spatially isolated Majorana zero modes, are easily found by preparing a contour plot of conductance over several independent tuning parameters, mimicking the effect of Zeeman splitting and voltages on gates near the junction. This suggests that, even stable apparently quantized conductance peaks need not correspond to isolated Majorana modes; rather, the $\textit{a priori}$ expectation should be that such quantized peaks generically occupy a significant fraction of the high-dimensional tuning parameter space that characterizes the NS tunneling experiments.


title: Quantum Cosmological Backreactions III: Deparametrised Quantum Cosmological Perturbation Theory

abstract: This is the third paper in a series of four in which we use space adiabatic methods in order to incorporate backreactions among the homogeneous and between the homogeneous and inhomogeneous degrees of freedom in quantum cosmological perturbation theory. In this paper we consider a particular kind of cosmological perturbation theory which starts from a gauge fixed version of General Relativity. The gauge fixing is performed using a material reference system called Gaussian dust. The resulting system has no constraints any more but possesses a physical Hamiltonian that drives the dynamics of both geometry and matter. As observable matter content we restrict to a scalar field (inflaton). We then explore the sector of that theory which is purely homogeneous and isotropic with respect to the geometry degrees of freedom but contains inhomogeneous perturbations up to second order of the scalar field. The purpose of this paper is to explore the quantum field theoretical challenges of the space adiabatic framework in a cosmological model of inflation which is technically still relatively simple. We compute the quantum backreaction effects from every energy band of the inhomogeneous matter modes on the evolution of the homogeneous geometry up to second order in the adiabatic parameter. These contributions turn out to be significant due to the infinite number of degrees of freedom and are very sensitive to the choice of Fock representation chosen for the inhomogeneous matter modes.


title: Multi-shot Echo Planar Imaging for accelerated Cartesian MR Fingerprinting: an alternative to conventional spiral MR Fingerprinting

abstract: Purpose: To develop an accelerated Cartesian MRF implementation using a multi-shot EPI sequence for rapid simultaneous quantification of T1 and T2 parameters. Methods: The proposed Cartesian MRF method involved the acquisition of highly subsampled MR images using a 16-shot EPI readout. A linearly varying flip angle train was used for rapid, simultaneous T1 and T2 quantification. The accuracy of parametric map estimations were improved by using an iterative projection algorithm. The results were compared to a conventional spiral MRF implementation. The acquisition time per slice was 8s and this method was validated on a phantom and a healthy volunteer brain in vivo. Results: Joint T1 and T2 estimations using the 16-shot EPI readout are in good agreement with the spiral implementation using the same acquisition parameters (deviation less than 3% for T1 and less than 4% for T2) for the healthy volunteer brain. The T1 and T2 values also agree with the conventional values previously reported in the literature. The visual quality of the multi-parametric maps generated by the multi-shot EPI-MRF and spiral-MRF implementations were comparable. Conclusion: The multi-shot EPI-MRF method generated accurate quantitative multi-parametric maps similar to conventional Spiral - MRF. This multi-shot approach achieved provides an alternative for performing MRF using an accelerated Cartesian readout, thereby increasing the potential usability of MRF.


title: Anticanonical tropical cubic del Pezzos contain exactly 27 lines

abstract: The classical statement of Cayley-Salmon that there are 27 lines on every smooth cubic surface in P^3 fails to hold under tropicalization: a tropical cubic surface in TP^3 often contains infinitely many tropical lines. Under mild genericity assumptions, we show that when embedded using the Eckardt triangles in the anticanonical system, tropical cubic del Pezzo surfaces contain exactly 27 tropical lines. In the non-generic case, which we identify explicitly, we find up to 27 extra lines, no multiple of which lifts to a curve on the cubic surface. We realize the moduli space of stable anticanonical tropical cubics as a four-dimensional fan in R^40 with an action of the Weyl group W(E_6). In the absence of Eckardt points, we show the combinatorial types of these tropical surfaces are determined by the boundary arrangement of 27 metric trees corresponding to the tropicalization of the classical 27 lines on the smooth algebraic cubic surfaces. Tropical convexity and the combinatorics of the root system E_6 play a central role in our analysis.


title: Parameter sensitivity analysis of dynamic ice sheet models

abstract: The velocity field and the height at the surface of a dynamic ice sheet are observed. The ice sheets are modeled by the full Stokes equations and shallow shelf/shelfy stream approximations. Time dependence is introduced by a kinematic free surface equation which updates the surface elevation using the velocity solution. The sensitivity of the observed quantities at the ice surface to parameters in the models, for example the basal topography and friction coefficients, is analyzed by first deriving the time dependent adjoint equations. Using the adjoint solutions, the effect of a perturbation in a parameter is obtained showing the importance of including the time dependence, in particular when the height is observed. The adjoint equations are solved analytically and numerically and the sensitivity of the desired parameters is determined in several examples in two dimensions. A closed form of the analytical solutions to the adjoint equations is given for a two dimensional grounding line migration problem in steady state under the shallow shelf approximation.


title: Robust Clustering Using Tau-Scales

abstract: K means is a popular non-parametric clustering procedure introduced by Steinhaus (1956) and further developed by MacQueen (1967). It is known, however, that K means does not perform well in the presence of outliers. Cuesta-Albertos et al (1997) introduced a robust alternative, trimmed K means, which can be tuned to be robust or efficient, but cannot achieve these two properties simultaneously in an adaptive way. To overcome this limitation we propose a new robust clustering procedure called K Tau Centers, which is based on the concept of Tau scale introduced by Yohai and Zamar (1988). We show that K Tau Centers performs well in extensive simulation studies and real data examples. We also show that the centers found by the proposed method are consistent estimators of the "true" centers defined as the minimizers of the the objective function at the population level.


title: Wannier basis method for KAM effect in quantum mechanics

abstract: The effect of Kolmogorov-Arnold-Moser (KAM) theorem in quantum systems is manifested in dividing eigenstates into regular and irregular states. We propose an effective method based on Wannier basis in phase space to illustrate this division of eigenstates. The quantum kicked-rotor model is used to illustrate this method, which allows us to define the area and effective dimension of each eigenstate to distinguish quantitatively regular and irregular eigenstates. This Wannier basis method also allows us to define the length of a Planck cell in the spectrum that measures how many Planck cells the system will traverse if it starts at the given Planck cell. Moreover, with this Wannier approach, we are able to clarify the distinction between KAM effect and Anderson localization.


title: Automated Definition of Skeletal Disease Burden in Metastatic Prostate Carcinoma: a 3D analysis of SPECT/CT images

abstract: To meet the current need for skeletal tumor-load estimation in prostate cancer (mCRPC), we developed a novel approach, based on adaptive bone segmentation. In this study, we compared the program output with existing estimates and with the radiological outcome. Seventy-six whole-body 99mTc-DPD-SPECT/CT from mCRPC patients were analyzed. The software identified the whole skeletal volume (SVol) and classified it voxels metastases (MVol) or normal bone (BVol). SVol was compared with the estimation of a commercial software. MVol was compared with manual assessment and with PSA-level. Counts/voxel were extracted from MVol and BVol. After six cycles of 223RaCl2-therapy every patient was re-evaluated as progressing (PD), stabilized (SD) or responsive (PR). SVol correlated with the one of the commercial software (R=0,99, p<0,001). MVol correlated with manually-counted lesions (R=0,61, p<0,001) and PSA (R=0,46, p<0.01). PD had a lower counts/voxel in MVol than PR/SD (715 \pm 190 Vs. 975 \pm 215 and 1058 \pm 255, p<0,05 and p<0,01) and in BVol (PD 275 \pm 60, PR 515 \pm 188 and SD 528 \pm 162 counts/voxel, p<0,001). Segmentation-based tumor load correlated with radiological/laboratory indices. Uptake was linked with the clinical outcome, suggesting that metastases in PD patients have a lower affinity for bone-seeking radionuclides and might benefit less from bone-targeted radioisotope therapies.


title: Rotation sensing in two coupled whispering-gallery-mode resonators with loss and gain

abstract: We theoretically propose a scheme to realize rotation sensing based on two coupled whispering-gallery-mode resonators with loss and gain. We consider that the active resonator with gain is rotated while the passive one with loss is stationary. The rotation will induce Sagnac effect and we show that the eigenfrequencies of the supermodes are sensitive to the Sagnac-Fizeau shift. Therefore, we can measure the average photon number in the steady state or the fluctuation spectrum of the output fields to detect the angular velocity of the rotation. We hope that our investigation will be useful in the design of quantum gyroscope based on spinning resonators.


title: A Grasping-centered Analysis for Cloth Manipulation

abstract: Compliant and soft hands have gained a lot of attention in the past decade because of their ability to adapt to the shape of the objects, increasing their effectiveness for grasping. However, when it comes to grasping highly flexible objects such as textiles, we face the dual problem: it is the object that will adapt to the shape of the hand or gripper. In this context, the classic grasp analysis or grasping taxonomies are not suitable for describing textile objects grasps. This work proposes a novel definition of textile object grasps that abstracts from the robotic embodiment or hand shape and recovers concepts from the early neuroscience literature on hand prehension skills. This framework enables us to identify what grasps have been used in literature until now to perform robotic cloth manipulation, and allows for a precise definition of all the tasks that have been tackled in terms of manipulation primitives based on regrasps. In addition, we also review what grippers have been used. Our analysis shows how the vast majority of cloth manipulations have relied only on one type of grasp, and at the same time we identify several tasks that need more variety of grasp types to be executed successfully. Our framework is generic, provides a classification of cloth manipulation primitives and can inspire gripper design and benchmark construction for cloth manipulation.


title: Thermodynamics of weakly coherent collisional models

abstract: We introduce the idea of weakly coherent collisional models, where the elements of an environment interacting with a system of interest are prepared in states that are approximately thermal, but have an amount of coherence proportional to a short system-environment interaction time in a scenario akin to well-known collisional models. We show that, in the continuous-time limit, the model allows for a clear formulation of the first and second laws of thermodynamics, which are modified to include a non-trivial contribution related to quantum coherence. Remarkably, we derive a bound showing that the degree of such coherence in the state of the elements of the environment represents a resource, which can be consumed to convert heat into an ordered (unitary-like) energy term in the system, even though no work is performed in the global dynamics. Our results therefore represent an instance where thermodynamics can be extended beyond thermal systems, opening the way for combining classical and quantum resources.


title: A Novel DDoS Attack Detection Method Using Optimized Generalized Multiple Kernel Learning

abstract: Distributed Denial of Service (DDoS) attack has become one of the most destructive network attacks which can pose a mortal threat to Internet security. Existing detection methods can not effectively detect early attacks. In this paper, we propose a detection method of DDoS attacks based on generalized multiple kernel learning (GMKL) combining with the constructed parameter R. The super-fusion feature value (SFV) and comprehensive degree of feature (CDF) are defined to describe the characteristic of attack flow and normal flow. A method for calculating R based on SFV and CDF is proposed to select the combination of kernel function and regularization paradigm. A DDoS attack detection classifier is generated by using the trained GMKL model with R parameter. The experimental results show that kernel function and regularization parameter selection method based on R parameter reduce the randomness of parameter selection and the error of model detection, and the proposed method can effectively detect DDoS attacks in complex environments with higher detection rate and lower error rate.


title: Absolute dimensions of the low-mass eclipsing binary system NSVS 10653195

abstract: Low-mass stars in eclipsing binary systems show radii larger and effective temperatures lower than theoretical stellar models predict for isolated stars with the same masses. Eclipsing binaries with low-mass components are hard to find due to their low luminosity. As a consequence, the analysis of the known low-mass eclipsing systems is key to understand this behavior. We developed a physical model of the LMDEB system NSVS 10653195 to accurately measure the masses and radii of the components. We obtained several high-resolution spectra in order to fit a spectroscopic orbit. Standardized absolute photometry was obtained to measure reliable color indices and to measure the mean Teff of the system in out-of-eclipse phases. We observed and analyzed optical VRI and infrared JK band differential light-curves which were fitted using PHOEBE. A Markov-Chain Monte Carlo (MCMC) simulation near the solution found provides robust uncertainties for the fitted parameters. NSVS 10653195 is a detached eclipsing binary composed of two similar stars with masses of M1=0.6402+/-0.0052 Msun and M2=0.6511+/-0.0052 Msun and radii of R1=0.687^{+0.017}_{-0.024} Rsun and R2=0.672^{+0.018}_{-0.022} Rsun. Spectral types were estimated to be K6V and K7V. These stars rotate in a circular orbit with an orbital inclination of i=86.22+/-0.61 degrees and a period of P=0.5607222(2) d. The distance to the system is estimated to be d=135.2^{+7.6}_{-7.9} pc, in excellent agreement with the value from Gaia. If solar metallicity were assumed, the age of the system would be older than log(age)~8 based on the Mbol-log Teff diagram. NSVS 10653195 is composed of two oversized and active K stars. While their radii is above model predictions their Teff are in better agreement with models.


title: Killings of social leaders in the Colombian post-conflict: Data analysis for investigative journalism

abstract: After the peace agreement of 2016 with FARC, the killings of social leaders have emerged as an important post-conflict challenge for Colombia. We present a data analysis based on official records obtained from the Colombian General Attorney's Office spanning the time period from 2012 to 2017. The results of the analysis show a drastic increase in the officially recorded number of killings of democratically elected leaders of community organizations, in particular those belonging to Juntas de Acci\'on Comunal [Community Action Boards]. These are important entities that have been part of the Colombian democratic apparatus since 1958, and enable communities to advocate for their needs. We also describe how the data analysis guided a journalistic investigation that was motivated by the Colombian government's denial of the systematic nature of social leaders killings.


title: Variational Fair Clustering

abstract: We propose a general variational framework of fair clustering, which integrates an original Kullback-Leibler (KL) fairness term with a large class of clustering objectives, including prototype or graph based. Fundamentally different from the existing combinatorial and spectral solutions, our variational multi-term approach enables to control the trade-off levels between the fairness and clustering objectives. We derive a general tight upper bound based on a concave-convex decomposition of our fairness term, its Lipschitz-gradient property and the Pinsker's inequality. Our tight upper bound can be jointly optimized with various clustering objectives, while yielding a scalable solution, with convergence guarantee. Interestingly, at each iteration, it performs an independent update for each assignment variable. Therefore, it can be easily distributed for large-scale datasets. This scalability is important as it enables to explore different trade-off levels between the fairness and clustering objectives. Unlike spectral relaxation, our formulation does not require computing its eigenvalue decomposition. We report comprehensive evaluations and comparisons with state-of-the-art methods over various fair-clustering benchmarks, which show that our variational formulation can yield highly competitive solutions in terms of fairness and clustering objectives.


title: Clock synchronization over networks -- Identifiability of the sawtooth model

abstract: In this paper, we analyze the two-node joint clock synchronization and ranging problem. We focus on the case of nodes that employ time-to-digital converters to determine the range between them precisely. This specific design choice leads to a sawtooth model for the captured signal, which has not been studied before from an estimation theoretic standpoint. In the study of this model, we recover the basic conclusion of a well-known article by Freris, Graham, and Kumar in clock synchronization. More importantly, we discover a surprising identifiability result on the sawtooth signal model: noise improves the theoretical condition of the estimation of the phase and offset parameters. To complete our study, we provide performance references for joint clock synchronization and ranging using the sawtooth signal model by presenting an exhaustive simulation study on basic estimation strategies under different realistic conditions. With our contributions in this paper, we enable further research in the estimation of sawtooth signal models and pave the path towards their industrial use for clock synchronization and ranging.


title: Parameter sensitivity analysis of dynamic ice sheet models-Numerical computations

abstract: The friction coefficient and the base topography of a stationary and a dynamic ice sheet are perturbed in two models for the ice: the full Stokes equations and the shallow shelf approximation. The sensitivity to the perturbations of the velocity and the height at the surface is quantified by solving the adjoint equations of the stress and the height equations providing weights for the perturbed data. The adjoint equations are solved numerically and the sensitivity is computed in several examples in two dimensions. Comparisons are made with analytical solutions to simplified problems.


title: Silicate features in the circumstellar envelopes of the Class~I binary driving source of HH250

abstract: We investigate the silicate feature of the two Class I components of HH250-IRS, a resolved binary system with a separation of $0''53$ driving a Herbig-Haro flow. Each component has its own circumstellar envelope, and the system is surrounded by a circumbinary disk. We have carried out low resolution spectroscopy in the 8-13$\mu$m range using VISIR at ESO's Very Large Telescope. The silicate features of both sources are clearly different. The NW component has a broad, smooth absorption profile lacking structure. The SE component shows the silicate feature in emission, with structure longwards of 9.5$\mu$m indicating the presence of crystalline dust in the dominant form of forsterite. The apparent lack of an absorption feature caused by foreground dust is probably due to the filling of the band with emission by amorphous silicates in the envelope of the object. Despite their virtually certain coevality, the differences in the components of the HH250-IRS binary are most likely due to markedly different circumstellar environments. The NW component displays an unevolved envelope, whereas dust growth and crystallization has taken place in the SE component. The weak or absent signatures of enstatite in the latter are fairly unusual among envelopes with crystalline dust, and we tentatively relate it to a possible wide gap or an inner truncation of the disk already hinted in previous observations by a drop in the $L'$-band flux, which might indicate that the SE component could actually be a very close binary. We speculate that the clear differences between the silicate feature spectra of both components of HH250-IRS may be due either to disk evolution sped up by multiplicity, or by accretion variability leading to episodes of crystal formation.


title: Optical Wireless cabin communication system

abstract: Offering a high data rate to users inside the cabin is attractive for the aircraft building companies. This paper presents an optical wireless system that uses visible light for transmitting data. The reading light for each seat is used as a transmitter in this work to offer high data rates. Red, yellow, green, and blue (RYGB) laser diode (LD) are utilised in each reading light to obtain a high modulation bandwidth. Two types of reading lights based on angle diversity transmitter (ADT) units are utilised in this paper for illumination and communication. Additionally, two kinds of optical receivers are used: four branches angle diversity receiver (ADR) and 25 pixels imaging receiver (ImR). The delay spread and the signal to interference and noise ratio (SINR) are evaluated. The proposed system can offer high data rates up to 22.8 Gbps for each user by using simple on-off-keying (OOK) modulation.


title: Co-existence of Micro, Pico and Atto Cells in Optical Wireless Communication

abstract: Interference between cells or users can have a significant impact on the quality of optical wireless communication (OWC) links. This paper studies the co-existence of infrared based Micro cells with Visible light communication based Pico and Atto Cells for downlink communication. The signal to noise ratio (SNR) of each cell and the signal to interference and noise ratio (SINR) between cells are evaluated when an angle diversity receiver and an imaging receiver are used. The results show that Atto cell systems provide higher SNR and data rates compared to the Pico cell systems which provide higher SNR and data rates compared to Micro cell systems. The Micro cell systems however provide mobility as they provide a larger coverage area, followed by the Pico cell systems.


title: Progress in global helioseismology: a new light on the solar modelling problem and its implications for solar-like stars

abstract: Since the first observations of solar oscillations, helioseismology has been one of the most successful fields of astrophysics. Data of high quality were obtained through the implementation of networks of ground-based observatories such as the GONG project or the BiSON network, coupled with space-based telescopes such as SOHO and SDO missions. Besides the improvement of observational data, solar seismologists developed sophisticated techniques to infer the internal structure of the Sun. These methods, then already extensively used in the field of Geophysics, are called inversion techniques. They allowed to determine the position of the solar convective envelope, its helium abundance and the internal radial profiles of thermodynamic quantities. Back in 1990s these comparisons showed a very high agreement between solar models and the Sun. However, the downward revision of the CNO surface abundances in the Sun induced a drastic reduction of this agreement leading to the solar modelling problem. More than ten years later, in the era of the space-based photometry missions which have established asteroseismology of solar-like stars as a standard approach to obtain their masses, radii and ages, the solar modelling problem still awaits a solution. We will present the results of new helioseismic inversions, discuss the current uncertainties of solar models as well as some possible solutions to the solar problem. We will show how helioseismology can help us grasp what is amiss in our models. We will also show that, far from being an argument about details of solar models, the solar problem has significant implications for seismology of solar-like stars, on the main sequence and beyond, impacting asteroseismology as a whole as well as the fields requiring precise and accurate knowledge of stellar masses, radii and ages, such as Galactic archaeology and exoplanetology.


title: Weak-field limit of Kaluza-Klein model with non-linear perfect fluid

abstract: The main purpose of our paper is to construct a viable Kaluza-Klein model satisfying the observable constraints. To this end, we investigate the six-dimensional model with spherical compactification of the internal space. Background matter is considered in the form of a perfect fluid with non-linear equations of state both in the external/our and internal spaces and the model is set to include an additional bare cosmological constant $\Lambda_6$. In the weak-field approximation, the background is perturbed by pressureless gravitating mass that is a static point-like particle. The non-linearity of the equations of state of a perfect fluid makes it possible to solve simultaneously a number of problems. The demand that the parameterized post-Newtonian parameter $\gamma$ be equal to 1 in this configuration, first, ensures compatibility with gravitational tests in the Solar system (deflection of light and time delay of radar echoes) at the same level of accuracy as General Relativity. Second, it translates into the absence of internal space variations so that the gravitational potential coincides exactly with the Newtonian one, securing the absence of the fifth force. Third, the gravitating mass remains pressurless in the external space as in the standard approach to non-relativistic astrophysical objects and meanwhile, acquires effective tension in the internal space.


title: Bayesian Learning from Sequential Data using Gaussian Processes with Signature Covariances

abstract: We develop a Bayesian approach to learning from sequential data by using Gaussian processes (GPs) with so-called signature kernels as covariance functions. This allows to make sequences of different length comparable and to rely on strong theoretical results from stochastic analysis. Signatures capture sequential structure with tensors that can scale unfavourably in sequence length and state space dimension. To deal with this, we introduce a sparse variational approach with inducing tensors. We then combine the resulting GP with LSTMs and GRUs to build larger models that leverage the strengths of each of these approaches and benchmark the resulting GPs on multivariate time series (TS) classification datasets. Code available at https://github.com/tgcsaba/GPSig.


title: Free action and cyclic sieving on skew semi-standard Young tableaux

abstract: In this note, we provide a short proof of Theorem 3.3 in the paper titled \emph{Crystals, semistandard tableaux and cyclic sieving phenomenon}, by Y.-T.~Oh and E.~Park, which concerns a cyclic sieving phenomenon on semi-standard Young tableaux. We also extend the result to skew shapes.


title: Cosmic Conundra Explained by Thermal History and Primordial Black Holes

abstract: A universal mechanism may be responsible for several unresolved cosmic conundra. The sudden drop in the pressure of relativistic matter at $W^{\pm}/Z^{0}$ decoupling, the quark--hadron transition and $e^{+}e^{-}$ annihilation enhances the probability of primordial black hole (PBH) formation in the early Universe. Assuming the amplitude of the primordial curvature fluctuations is approximately scale-invariant, this implies a multi-modal PBH mass spectrum with peaks at $10^{-6}$, 1, 30, and $10^{6}\,M_{\odot}$. This suggests a unified PBH scenario which naturally explains the dark matter and recent microlensing observations, the LIGO/Virgo black hole mergers, the correlations in the cosmic infrared and X-ray backgrounds, and the origin of the supermassive black holes in galactic nuclei at high redshift. A distinctive prediction of our model is that LIGO/Virgo should observe black hole mergers in the mass gaps between 2 and $5\,M_{\odot}$ (where no stellar remnants are expected) and above $65\,M_{\odot}$ (where pair-instability supernovae occur) and low-mass-ratios in between. Therefore the recent detection of events GW190425, GW190814 and GW190521 with these features is striking confirmation of our prediction and may indicate a primordial origin for the black holes. In this case, the exponential sensitivity of the PBH abundance to the equation of state would offer a unique probe of the QCD phase transition. The detection of PBHs would also offer a novel way to probe the existence of new particles or phase transitions with energy between $1\,{\rm MeV}$ and $10^{10}\,$GeV.


title: Exponential asymptotics for the eigenvalues in the broken PT-symmetric region

abstract: Stemming from the seminal work of Bender & Boettcher in 1998 (Phys. Rev. Lett. vol. 80 pp. 5243-5246), there has been great interest in the study of PT-symmetric models of quantum mechanics, where the primary focus is with the study of non-Hermitian Hamiltonians that nevertheless produce countably infinite sets of real-valued eigenvalues. One of the fundamental models of such a system is governed by the Hamiltonian $H = \hat{p}^2 + x^2(ix)^{\varepsilon}$. In their work, Bender & Boettcher proposed a WKB methodology for the prediction of the discrete eigenvalues in the so-called unbroken region of $\varepsilon > 0$. However, the authors noted that this methodology fails to predict those 'broken' eigenvalues for $\varepsilon < 0$. Here, we shall explain why the traditional WKB methodology fails, and we shall demonstrate how eigenvalues for all relevant values of $\varepsilon$ can be predicted using techniques in exponential asymptotics. These predictions provide excellent agreement to exact numerical results over nearly the entire range of values. Moreover, such techniques can be extended to a much wider range PT-symmetric problems.


title: NGTS-7Ab: An ultra-short period brown dwarf transiting a tidally-locked and active M dwarf

abstract: We present the discovery of NGTS-7Ab, a high mass brown dwarf transiting an M dwarf with a period of 16.2 hours, discovered as part of the Next Generation Transit Survey (NGTS). This is the shortest period transiting brown dwarf around a main or pre-main sequence star to date. The M star host (NGTS-7A) has an age of roughly 55 Myr and is in a state of spin-orbit synchronisation, which we attribute to tidal interaction with the brown dwarf acting to spin up the star. The host star is magnetically active and shows multiple flares across the NGTS and follow up lightcurves, which we use to probe the flare-starspot phase relation. The host star also has an M star companion at a separation of 1.13 arcseconds with very similar proper motion and systemic velocity, suggesting the NGTS-7 system is a hierarchical triple. The combination of tidal synchronisation and magnetic braking is expected to drive ongoing decay of the brown dwarf orbit, with a remaining lifetime of only 5-10 Myr.


title: Non-Relativistic Four Dimensional p-Brane Supersymmetric Theories and Lie Algebra Expansion

abstract: We apply the Lie algebra expansion method to the $\mathcal{N}=1$ super-Poincar\'e algerba in four dimensions. We define a set of p-brane projectors that induce a decomposition of the super-Poincar\'e algebra preparatory for the expansion. We show that starting from the $\mathcal{N}=1$ supergravity action in four dimensions it is possible to obtain two non-relativistic supersymmetric theories, one describing strings, the other membranes.


title: New constraints on Lorentz Invariance violation from Crab Nebula spectrum beyond $100$ TeV

abstract: Recently two collaborations, Tibet and HAWC, presented new measurements of gamma-ray spectrum from Crab Nebula \cite{Amenomori:2019rjd},\cite{Abeysekara:2019edl} which continues beyond $100$ TeV. We use these data to establish two-sided constraints on parameters of Lorentz Invariance violation in quantum electrodynamics. The limits on Lorentz violating mass scale for quartic dispersion relation are $4.1\times 10^{14}\ \mbox{GeV}$ (photon splitting) and $1.9\times 10^{13}\ \mbox{GeV}$ (photon decay) for superluminal case, and $1.4\times 10^{12}$ GeV (suppression of shower formation) for subluminal case.


title: Deep Fuzzy Systems

abstract: An investigation of deep fuzzy systems is presented in this paper. A deep fuzzy system is represented by recursive fuzzy systems from an input terminal to output terminal. Recursive fuzzy systems are sequences of fuzzy grade memberships obtained using fuzzy transmition functions and recursive calls to fuzzy systems. A recursive fuzzy system which calls a fuzzy system n times includes fuzzy chains to evaluate the final grade membership of this recursive system. A connection matrix which includes recursive calls are used to represent recursive fuzzy systems.


title: Measurement of the total neutron scattering cross section ratios of noble gases of natural isotopic composition using a pulsed neutron beam

abstract: Precision measurements of slow neutron cross sections with atoms have several scientific applications. In particular the n-$^{4}$He s-wave scattering length is important to know both for helping to constrain the nuclear three-body interaction and for the proper interpretation of several ongoing slow neutron experiments searching for other types of neutron-atom interactions. We present new measurements of the ratios of the neutron differential scattering cross sections for natural isotopic-abundance mixtures of the noble gases He, Ar, Kr, and Xe to natural isotopic abundance Ne. These measurements were performed using a recently developed neutron scattering apparatus for gas samples located on a pulsed slow neutron beamline which was designed to search for possible exotic neutron-atom interactions and employs both neutron time of flight information and a position-sensitive neutron detector for scattering event reconstruction. We found agreement with the literature values of scattering cross sections inferred from Ar/Ne, Kr/Ne and Xe/Ne differential cross section ratios over the $q$ range of $1 - 7$ nm$^{-1}$. However for the case of He/Ne we find that the cross section inferred differs by 11.3% (7.6 $\sigma$) from previously-reported values inferred from neutron phase shift measurements, but is in reasonable agreement with values from other measurements. The very large discrepancy in the He/Ne ratio calls for a new precision measurement of the n-$^{4}$He scattering length using neutron interferometry.


title: Phonon scattering induced carrier resistivity in twisted double bilayer graphene

abstract: In this work we carry out a theoretical study of the phonon-induced resistivity in twisted double bilayer graphene (TDBG), in which two Bernal-stacked bilayer graphene devices are rotated relative to each other by a small angle $\theta$. We show that at small twist angles ($\theta\sim 1^\circ$) the effective mass of the TDBG system is greatly enhanced, leading to a drastically increased phonon-induced resistivity in the high-temperature limit where phonon scattering leads to a linearly increasing resistivity with increasing temperature. We also discuss possible implications of our theory on superconductivity in such a system, and provide an order of magnitude estimation of the superconducting transition temperature.


title: Scale-invariant estimates and vorticity alignment for Navier-Stokes in the half-space with no-slip boundary conditions

abstract: This paper is concerned with geometric regularity criteria for the Navier-Stokes equations in $\mathbb{R}^3_{+}\times (0,T)$ with no-slip boundary condition, with the assumption that the solution satisfies the `ODE blow-up rate' Type I condition. More precisely, we prove that if the vorticity direction is uniformly continuous on subsets of $$\bigcup_{t\in(T-1,T)} \big(B(0,R)\cap\mathbb{R}^3_{+}\big)\times {\{t\}},\,\,\,\,\,\, R=O(\sqrt{T-t})$$ where the vorticity has large magnitude, then $(0,T)$ is a regular point. This result is inspired by and improves the regularity criteria given by Giga, Hsu and Maekawa (2014). We also obtain new local versions for suitable weak solutions near the flat boundary. Our method hinges on new scaled Morrey estimates, blow-up and compactness arguments and `persistence of singularites' on the flat boundary. The scaled Morrey estimates seem to be of independent interest.


title: Unsupervised State Representation Learning in Atari

abstract: State representation learning, or the ability to capture latent generative factors of an environment, is crucial for building intelligent agents that can perform a wide variety of tasks. Learning such representations without supervision from rewards is a challenging open problem. We introduce a method that learns state representations by maximizing mutual information across spatially and temporally distinct features of a neural encoder of the observations. We also introduce a new benchmark based on Atari 2600 games where we evaluate representations based on how well they capture the ground truth state variables. We believe this new framework for evaluating representation learning models will be crucial for future representation learning research. Finally, we compare our technique with other state-of-the-art generative and contrastive representation learning methods. The code associated with this work is available at https://github.com/mila-iqia/atari-representation-learning


title: Local Bures-Wasserstein Transport: A Practical and Fast Mapping Approximation

abstract: Optimal transport (OT)-based methods have a wide range of applications and have attracted a tremendous amount of attention in recent years. However, most of the computational approaches of OT do not learn the underlying transport map. Although some algorithms have been proposed to learn this map, they rely on kernel-based methods, which makes them prohibitively slow when the number of samples increases. Here, we propose a way to learn an approximate transport map and a parametric approximation of the Wasserstein barycenter. We build an approximated transport mapping by leveraging the closed-form of Gaussian (Bures-Wasserstein) transport; we compute local transport plans between matched pairs of the Gaussian components of each density. The learned map generalizes to out-of-sample examples. We provide experimental results on simulated and real data, comparing our proposed method with other mapping estimation algorithms. Preliminary experiments suggest that our proposed method is not only faster, with a factor 80 overall running time, but it also requires fewer components than state-of-the-art methods to recover the support of the barycenter. From a practical standpoint, it is straightforward to implement and can be used with a conventional machine learning pipeline.


title: Enhanced Ehlers Transformation and the Majumdar-Papapetrou-NUT Spacetime

abstract: The transformation which adds (or removes) NUT charge when it is applied to electrovacuum, axisymmetric and stationary space-times is studied. After analysing the Ehlers and the Reina-Treves transformations we propose a new one, more precise in the presence of the Maxwell electromagnetic field. The enhanced Ehlers transformation proposed turns out to act as a gravitomagnetic duality, analogously to the electromagnetic duality, but for gravity: it rotates the mass charge into the gravomagnetic (or NUT) charge. As an example the Kerr-Newman-NUT black hole is obtained with the help of this enhanced transformation. Moreover a new analytical exact solution is built adding the NUT charge to a double charged black hole, at equilibrium. It describes the non-extremal generalisation of the Majumdar-Papapetrou-NUT solution. From the near-horizon analysis, its microscopic entropy, according to the Kerr/CFT correspondence, is found and the second law of black hole thermodynamics is discussed.


title: Mathematical modeling of a Cosserat method in finite-strain holonomic plasticity

abstract: This article deals with the mathematical derivation and the validation over benchmark examples of a numerical method for the solution of a finite-strain holonomic (rate-independent) Cosserat plasticity problem for materials, possibly with microstructure. Two improvements are made in contrast to earlier approaches: First, the micro-rotations are parameterized with the help of an Euler-Rodrigues formula related to quaternions. Secondly, as main result, a novel two-pass preconditioning scheme for searching the energy-minimizing solutions based on the limited memory Broyden-Fletcher-Goldstein-Shanno quasi-Newton method is proposed that consists of a predictor step and a corrector-iteration. After outlining the necessary adaptations to the model, numerical simulations compare the performance and efficiency of the new and the old algorithm. The proposed numerical model can be effectively employed for studying the mechanical response of complicated materials featuring large size effects.


title: Evaluating Protein Transfer Learning with TAPE

abstract: Protein modeling is an increasingly popular area of machine learning research. Semi-supervised learning has emerged as an important paradigm in protein modeling due to the high cost of acquiring supervised protein labels, but the current literature is fragmented when it comes to datasets and standardized evaluation techniques. To facilitate progress in this field, we introduce the Tasks Assessing Protein Embeddings (TAPE), a set of five biologically relevant semi-supervised learning tasks spread across different domains of protein biology. We curate tasks into specific training, validation, and test splits to ensure that each task tests biologically relevant generalization that transfers to real-life scenarios. We benchmark a range of approaches to semi-supervised protein representation learning, which span recent work as well as canonical sequence learning techniques. We find that self-supervised pretraining is helpful for almost all models on all tasks, more than doubling performance in some cases. Despite this increase, in several cases features learned by self-supervised pretraining still lag behind features extracted by state-of-the-art non-neural techniques. This gap in performance suggests a huge opportunity for innovative architecture design and improved modeling paradigms that better capture the signal in biological sequences. TAPE will help the machine learning community focus effort on scientifically relevant problems. Toward this end, all data and code used to run these experiments are available at https://github.com/songlab-cal/tape.


title: Holistic evaluation of XML queries with structural preferences on an annotated strong dataguide

abstract: With the emergence of XML as de facto format for storing and exchanging information over the Internet, the search for ever more innovative and effective techniques for their querying is a major and current concern of the XML database community. Several studies carried out to help solve this problem are mostly oriented towards the evaluation of so-called exact queries which, unfortunately, are likely (especially in the case of semi-structured documents) to yield abundant results (in the case of vague queries) or empty results (in the case of very precise queries). From the observation that users who make requests are not necessarily interested in all possible solutions, but rather in those that are closest to their needs, an important field of research has been opened on the evaluation of preferences queries. In this paper, we propose an approach for the evaluation of such queries, in case the preferences concern the structure of the document. The solution investigated revolves around the proposal of an evaluation plan in three phases: rewriting-evaluation-merge. The rewriting phase makes it possible to obtain, from a partitioning -transformation operation of the initial query, a hierarchical set of preferences path queries which are holistically evaluated in the second phase by an instrumented version of the algorithm TwigStack. The merge phase is the synthesis of the best results.


title: Representability of Chow groups of codimension three cycles

abstract: In this note we are going to prove that if we have a fibration of smooth projective varieties $X\to S$ over a surface $S$ such that $X$ is of dimension four and that the geometric generic fiber has finite dimensional motive and the first \'etale cohomology of the geometric generic fiber with respect to $\mathbb {Q}_l$ coefficients is zero and the second \'etale cohomology is spanned by divisors, then $A^3(X)$ (codimension three algebraically trivial cycles modulo rational equivalence) is dominated by finitely many copies of $A_0(S)$. Meaning that there exists finitely many correspondences $\Gamma_i$ on $S\times X$, such that $\sum_i \Gamma_i$ is surjective from $\oplus A^2(S)$ to $A^3(X)$.


title: On a question of Colliot-Th\'{e}l\`{e}ne on Chow group mod $n$

abstract: In this note we define the notion of Tate-Shafarevich group and Selmer group of the Chow group of an abelian variety defined over a number field. In this context we give positive answer to the question of Colliot-Th\'{e}l\`{e}ne that the Chow group of zero cycles modulo $n$ is finite over any number field.


title: On edge-ordered Ramsey numbers

abstract: An edge-ordered graph is a graph with a linear ordering of its edges. Two edge-ordered graphs are equivalent if their is an isomorphism between them preserving the ordering of the edges. The edge-ordered Ramsey number $r_{edge}(H; q)$ of an edge-ordered graph $H$ is the smallest $N$ such that there exists an edge-ordered graph $G$ on $N$ vertices such that, for every $q$-coloring of the edges of $G$, there is a monochromatic subgraph of $G$ equivalent to $H$. Recently, Balko and Vizer announced that $r_{edge}(H;q)$ exists. However, their proof uses the Graham-Rothschild theorem and consequently gives an enormous upper bound on these numbers. We give a new proof giving a much better bound. We prove that for every edge-ordered graph $H$ on $n$ vertices, we have $r_{edge}(H;q) \leq 2^{c^qn^{2q-2}\log^q n}$, where $c$ is an absolute constant. We also explore the edge-ordered Ramsey number of sparser graphs and prove a polynomial bound for edge-ordered graphs of bounded degeneracy. We also prove a strengthening for edge-labeled graphs, graphs where every edge is given a label and the labels do not necessary have an ordering.


title: The Prospects of Observing Tidal Disruption Events with the LSST

abstract: The upcoming Large Synoptic Survey Telescope (LSST) will observe 18 000 deg$^2$ of the Southern sky and is expected to discover thousands of transients every night due to its large coverage of the sky and its observing strategy. In this work we address the prospects of the LSST in discovering Tidal Disruption Events (TDEs) and in probing the supermassive black hole (SMBH) mass distribution in the Universe. We used the LSST simulation framework and defined TDE catalogs on 20 fields of 20.25 deg$^2$ size. TDE properties were defined by randomly chosen impact factors and SMBH masses drawn from six different mass distributions. Observations of TDEs in 10 years of LSST operations were simulated by querying the simulated observing strategy database minion_1016. Based on the results of our simulations we estimate that the LSST should discover between 35 000 and 80 000 TDEs in 10 years of operations, depending on the assumed SMBH mass distribution. We also find that probing the SMBH mass distribution with TDE observations will not be straightforward due to the fact that TDEs caused by low mass black holes ($10^ 5 M_\odot$) are expected to be less luminous and shorter than TDEs by heavier SMBHs ($> 10^6 M_\odot$), and therefore will mostly be missed by the irregular LSST cadence minion_1016.


title: PyRobot: An Open-source Robotics Framework for Research and Benchmarking

abstract: This paper introduces PyRobot, an open-source robotics framework for research and benchmarking. PyRobot is a light-weight, high-level interface on top of ROS that provides a consistent set of hardware independent mid-level APIs to control different robots. PyRobot abstracts away details about low-level controllers and inter-process communication, and allows non-robotics researchers (ML, CV researchers) to focus on building high-level AI applications. PyRobot aims to provide a research ecosystem with convenient access to robotics datasets, algorithm implementations and models that can be used to quickly create a state-of-the-art baseline. We believe PyRobot, when paired up with low-cost robot platforms such as LoCoBot, will reduce the entry barrier into robotics, and democratize robotics. PyRobot is open-source, and can be accessed via https://pyrobot.org.


title: XLNet: Generalized Autoregressive Pretraining for Language Understanding

abstract: With the capability of modeling bidirectional contexts, denoising autoencoding based pretraining like BERT achieves better performance than pretraining approaches based on autoregressive language modeling. However, relying on corrupting the input with masks, BERT neglects dependency between the masked positions and suffers from a pretrain-finetune discrepancy. In light of these pros and cons, we propose XLNet, a generalized autoregressive pretraining method that (1) enables learning bidirectional contexts by maximizing the expected likelihood over all permutations of the factorization order and (2) overcomes the limitations of BERT thanks to its autoregressive formulation. Furthermore, XLNet integrates ideas from Transformer-XL, the state-of-the-art autoregressive model, into pretraining. Empirically, under comparable experiment settings, XLNet outperforms BERT on 20 tasks, often by a large margin, including question answering, natural language inference, sentiment analysis, and document ranking.


title: Anionic nanoparticle-lipid membrane interactions: the protonation of anionic ligands at the membrane surface reduces membrane disruption

abstract: Monolayer-protected gold nanoparticles (Au NPs) are promising biomedical tools with applications to diagnosis and therapy, thanks to their biocompatibility and versatility. Here we show how the NP surface functionalization can drive the mechanism of interaction with lipid membranes. In particular, we show that the spontaneous protonation of anionic carboxylic groups on the NP surface can make the NP-membrane interaction faster and less disruptive.


title: Reduced I/O Latency with Futures

abstract: Task parallelism research has traditionally focused on optimizing computation-intensive applications. Due to the proliferation of commodity parallel processors, there has been recent interest in supporting interactive applications. Such interactive applications frequently rely on I/O operations that may incur significant latency. In order to increase performance, when a particular thread of control is blocked on an I/O operation, ideally we would like to hide this latency by using the processing resources to do other ready work instead of blocking or spin waiting on this I/O. There has been limited prior work on hiding this latency and only one result that provides a theoretical bound for interactive applications that use I/Os. In this work, we propose a method for hiding the latency of I/O operations by using the futures abstraction. We provide a theoretical analysis of our algorithm that shows our algorithm provides better execution time guarantees than prior work. We also implemented the algorithm in a practically efficient prototype library that runs on top of the Cilk-F runtime, a runtime system that supports futures within the context of the Cilk Plus language, and performed experiments that demonstrate the efficiency of our implementation.


title: Neural Point-Based Graphics

abstract: We present a new point-based approach for modeling the appearance of real scenes. The approach uses a raw point cloud as the geometric representation of a scene, and augments each point with a learnable neural descriptor that encodes local geometry and appearance. A deep rendering network is learned in parallel with the descriptors, so that new views of the scene can be obtained by passing the rasterizations of a point cloud from new viewpoints through this network. The input rasterizations use the learned descriptors as point pseudo-colors. We show that the proposed approach can be used for modeling complex scenes and obtaining their photorealistic views, while avoiding explicit surface estimation and meshing. In particular, compelling results are obtained for scene scanned using hand-held commodity RGB-D sensors as well as standard RGB cameras even in the presence of objects that are challenging for standard mesh-based modeling.


title: Provable Gradient Variance Guarantees for Black-Box Variational Inference

abstract: Recent variational inference methods use stochastic gradient estimators whose variance is not well understood. Theoretical guarantees for these estimators are important to understand when these methods will or will not work. This paper gives bounds for the common "reparameterization" estimators when the target is smooth and the variational family is a location-scale distribution. These bounds are unimprovable and thus provide the best possible guarantees under the stated assumptions.


title: Comment on "Configuration interaction calculations and excitation rates of X-ray and EUV transitions in sulfurlike manganese" by El-Maaref et al. [J. Elect. Spectrosc. Related Phen. 215 (2017) 22]

abstract: In one of the papers, El-Maaref et al. [J. Elect. Spectrosc. Related Phen. 215 (2017) 22] have reported results for energy levels, oscillator strengths (f-values), radiative rates (A-values), collision strengths ($\Omega$), and excitation rates for transitions in S-like Mn~X. However, except for energy levels their results are restricted to only a few transitions and hence have limited application. Furthermore, most of their results have scope for improvement, but for $\Omega$ are not correct as these have been found to be overestimated by an order of magnitude. In this comment we discuss the discrepancies and deficiencies of their results and recommend that a fresh calculation should be performed for this ion.


title: A simple argument that small hydrogen may exist

abstract: This paper discusses a possible existence of small hydrogen, which may have been created during the Big Bang or large explosions in Universe.


title: Predicting Patent Citations to measure Economic Impact of Scholarly Research

abstract: A crucial goal of funding research and development has always been to advance economic development. On this basis, a consider-able body of research undertaken with the purpose of determining what exactly constitutes economic impact and how to accurately measure that impact has been published. Numerous indicators have been used to measure economic impact, although no single indicator has been widely adapted. Based on patent data collected from Altmetric we predict patent citations through various social media features using several classification models. Patents citing a research paper implies the potential it has for direct application inits field. These predictions can be utilized by researchers in deter-mining the practical applications for their work when applying for patents.


title: Detection of 4765 MHz OH Emission in a Pre-Planetary Nebula -- CRL 618

abstract: Jets and outflows are ubiquitous phenomena in astrophysics, found in our Galaxy in diverse environments, from the formation of stars to late-type stellar objects. We present observations conducted with the 305m Arecibo Telescope of the pre-planetary nebula CRL 618 (Westbrook Nebula) - a well studied late-type star that has developed bipolar jets. The observations resulted in the first detection of 4765 MHz OH in a late-type stellar object. The line was narrow (FWHM ~ 0.6 km/s) and ~40 km/s blueshifted with respect to the systemic velocity, which suggests association with the expanding jets/bullets in CRL 618. We also report non-detection at Arecibo of any other OH transition between 1 and 9 GHz. The non-detections were obtained during the observations in 2008, when the 4765 MHz OH line was first discovered, and also in 2015 when the 4765 MHz OH line was not detected. Our data indicate that the 4765 MHz OH line was a variable maser. Modeling of the 4765 MHz OH detection and non-detection of the other transitions is consistent with the physical conditions expected in CRL 618. The 4765 MHz OH maser could originate from dissociation of H2O by shocks after sublimation of icy objects in this dying carbon-rich stellar system, although other alternatives such as OH in an oxygen-rich circumstellar region associated with a binary companion are also possible.


title: Cortical Divisive Normalization from Wilson-Cowan Neural Dynamics

abstract: Divisive Normalization and the Wilson-Cowan equations are influential models of neural interaction and saturation [Carandini and Heeger Nat.Rev.Neurosci. 2012; Wilson and Cowan Kybernetik 1973]. However, they have not been analytically related yet. In this work we show that Divisive Normalization can be obtained from the Wilson-Cowan model. Specifically, assuming that Divisive Normalization is the steady state solution of the Wilson-Cowan differential equation, we find that the kernel that controls neural interactions in Divisive Normalization depends on the Wilson-Cowan kernel but also has a signal-dependent contribution. A standard stability analysis of a Wilson-Cowan model with the parameters obtained from our relation shows that the Divisive Normalization solution is a stable node. This stability demonstrates the consistency of our steady state assumption. The proposed theory provides a physiological foundation (a relation to a dynamical network with fixed wiring among neurons) for the functional suggestions that have been done on the need of signal-dependent Divisive Normalization [e.g. in Coen-Cagli et al., PLoS Comp.Biol. 2012]. Moreover, this theory explains the modifications that had to be introduced ad-hoc in Gaussian kernels of Divisive Normalization in [Martinez et al. Front. Neurosci. 2019] to reproduce contrast responses. The proposed relation implies that the Wilson-Cowan dynamics also reproduces visual masking and subjective image distortion metrics, which up to now had been mainly explained via Divisive Normalization. Finally, this relation allows to apply to Divisive Normalization the methods which up to now had been developed for dynamical systems such as Wilson-Cowan networks.


title: The 8-parameter Fisher-Bingham distribution on the sphere

abstract: The Fisher-Bingham distribution ($\mathrm{FB}_8$) is an eight-parameter family of probability density functions (PDF) on $S^2$ that, under certain conditions, reduce to spherical analogues of bivariate normal PDFs. Due to difficulties in computing its overall normalization constant, applications have been mainly restricted to subclasses of $\mathrm{FB}_8$, such as the Kent ($\mathrm{FB}_5$) or von Mises-Fisher (vMF) distributions. However, these subclasses often do not adequately describe directional data that are not symmetric along great circles. The normalizing constant of $\mathrm{FB}_8$ can be numerically integrated, and recently Kume and Sei showed that it can be computed using an adjusted holonomic gradient method. Both approaches, however, can be computationally expensive. In this paper, I show that the normalization of $\mathrm{FB}_8$ can be expressed as an infinite sum consisting of hypergeometric functions, similar to that of the $\mathrm{FB}_5$. This allows the normalization to be computed under summation with adequate stopping conditions. I then fit the $\mathrm{FB}_8$ to a synthetic dataset using a maximum-likelihood approach and show its improvements over a fit with the more restrictive $\mathrm{FB}_5$ distribution.


title: Proposal for a new quantum theory of gravity II: Spectral equation of motion for the atom of space-time-matter

abstract: In the first paper of this series, we have introduced the concept of an atom of space-time-matter [STM], which is described by the spectral action of non-commutative geometry, corresponding to a classical theory of gravity. In the present work, we use the Connes time parameter along with the spectral action, to incorporate gravity into trace dynamics. We then derive the spectral equation of motion for the STM atom, which turns out to be the Dirac equation on a non-commutative space.


title: Mutually Orthogonal Latin Squares based on Cellular Automata

abstract: We investigate sets of Mutually Orthogonal Latin Squares (MOLS) generated by Cellular Automata (CA) over finite fields. After introducing how a CA defined by a bipermutive local rule of diameter $d$ over an alphabet of $q$ elements generates a Latin square of order $q^{d-1}$, we study the conditions under which two CA generate a pair of orthogonal Latin squares. In particular, we prove that the Latin squares induced by two Linear Bipermutive CA (LBCA) over the finite field $\mathbb{F}_q$ are orthogonal if and only if the polynomials associated to their local rules are relatively prime. Next, we enumerate all such pairs of orthogonal Latin squares by counting the pairs of coprime monic polynomials with nonzero constant term and degree $n$ over $\mathbb{F}_q$. Finally, we present a construction of MOLS generated by LBCA with irreducible polynomials and prove the maximality of the resulting sets, as well as a lower bound which is asymptotically close to their actual number.


title: Electron optics in phosphorene pn junctions: Negative reflection and anti super-Klein tunneling

abstract: Ballistic electrons in phosphorene $pn$ junctions show optical-like phenomena. Phosphorene is modeled by a tight-binding Hamiltonian that describes its electronic structure at low energies, where the electrons behave in the zigzag direction as massive Dirac fermions and in the orthogonal armchair direction as Schr\"odinger electrons. Applying the continuum approximation, we derive the electron optics laws in phosphorene $pn$ junctions, which show very particular and unusual properties. Due to the anisotropy of the electronic structure, these laws depend strongly on the orientation of the junction with respect to the sublattice. Negative and anomalous reflection are observed for tilted junctions, while the typical specular reflection is found only, if the junction is parallel to the zigzag or armchair edges. Moreover, omni-directional total reflection, called anti-super Klein tunneling, is observed if the junction is parallel to the armchair edge. Applying the nonequilibrium Green's function method on the tight-binding model, we calculate numerically the current flow. The good agreement of both approaches confirms the atypical transport properties, which can be used in nano-devices to collimate and filter the electron flow, or to switch its direction.


title: Capturing strong correlations in spin, electron and local moment systems

abstract: We address the question of identifying degrees of freedom for quantum systems. Typically, quasi-particle descriptions of correlated matter are based upon the canonical algebras of bosons or fermions. Here we highlight that a special class of non-canonical algebras also offer useful quantum degrees of freedom, allowing for the development of quasi-particle descriptions which go beyond the weakly correlated paradigm. We give a broad overview of such algebras for spin, electron and local moment systems, and outline important test problems upon which to develop the framework.


title: Entanglement Entropy in Generalised Quantum Lifshitz Models

abstract: We compute universal finite corrections to entanglement entropy for generalised quantum Lifshitz models in arbitrary odd spacetime dimensions. These are generalised free field theories with Lifshitz scaling symmetry, where the dynamical critical exponent $z$ equals the number of spatial dimensions $d$, and which generalise the 2+1-dimensional quantum Lifshitz model to higher dimensions. We analyse two cases: one where the spatial manifold is a $d$-dimensional sphere and the entanglement entropy is evaluated for a hemisphere, and another where a $d$-dimensional flat torus is divided into two cylinders. In both examples the finite universal terms in the entanglement entropy are scale invariant and depend on the compactification radius of the scalar field.


title: When to Trust Your Model: Model-Based Policy Optimization

abstract: Designing effective model-based reinforcement learning algorithms is difficult because the ease of data generation must be weighed against the bias of model-generated data. In this paper, we study the role of model usage in policy optimization both theoretically and empirically. We first formulate and analyze a model-based reinforcement learning algorithm with a guarantee of monotonic improvement at each step. In practice, this analysis is overly pessimistic and suggests that real off-policy data is always preferable to model-generated on-policy data, but we show that an empirical estimate of model generalization can be incorporated into such analysis to justify model usage. Motivated by this analysis, we then demonstrate that a simple procedure of using short model-generated rollouts branched from real data has the benefits of more complicated model-based algorithms without the usual pitfalls. In particular, this approach surpasses the sample efficiency of prior model-based methods, matches the asymptotic performance of the best model-free algorithms, and scales to horizons that cause other model-based methods to fail entirely.


title: Resilient Leader-Follower Consensus with Time-Varying Leaders in Discrete-Time Systems

abstract: The problem of consensus in the presence of adversarially behaving agents has been studied extensively in the literature. The proposed algorithms typically guarantee that the consensus value lies within the convex hull of initial normal agents' states. In leader-follower consensus problems however, the objective for normally behaving agents is to track a time-varying reference state that may take on values outside of this convex hull. In this paper we present a method for agents with discrete-time dynamics to resiliently track a set of leaders' common time-varying reference state despite a bounded subset of the leaders and followers behaving adversarially. The efficacy of our results are demonstrated through simulations.


title: Training on test data: Removing near duplicates in Fashion-MNIST

abstract: MNIST and Fashion MNIST are extremely popular for testing in the machine learning space. Fashion MNIST improves on MNIST by introducing a harder problem, increasing the diversity of testing sets, and more accurately representing a modern computer vision task. In order to increase the data quality of FashionMNIST, this paper investigates near duplicate images between training and testing sets. Near-duplicates between testing and training sets artificially increase the testing accuracy of machine learning models. This paper identifies near-duplicate images in Fashion MNIST and proposes a dataset with near-duplicates removed.


title: Steinhaus Filtration and Stable Paths in the Mapper

abstract: Two central concepts from topological data analysis are persistence and the Mapper construction. Persistence employs a sequence of objects built on data called a filtration. A Mapper produces insightful summaries of data, and has found widespread applications in diverse areas. We define a new filtration called the cover filtration built from a single cover based on a generalized Steinhaus distance, which is a generalization of Jaccard distance. We prove a stability result: the cover filtrations of two covers are $\alpha/m$ interleaved, where $\alpha$ is a bound on bottleneck distance between covers and $m$ is the size of smallest set in either cover. We also show our construction is equivalent to the Cech filtration under certain settings, and the Vietoris-Rips filtration completely determines the cover filtration in all cases. We then develop a theory for stable paths within this filtration. Unlike standard results on stability in topological persistence, our definition of path stability aligns exactly with the above result on stability of cover filtration. We demonstrate how our framework can be employed in a variety of applications where a metric is not obvious but a cover is readily available. First we present a new model for recommendation systems using cover filtration. For an explicit example, stable paths identified on a movies data set represent sequences of movies constituting gentle transitions from one genre to another. As a second application in explainable machine learning, we apply the Mapper for model induction, providing explanations in the form of paths between subpopulations. Stable paths in the Mapper from a supervised machine learning model trained on the FashionMNIST data set provide improved explanations of relationships between subpopulations of images.


title: Safe and Near-Optimal Policy Learning for Model Predictive Control using Primal-Dual Neural Networks

abstract: In this paper, we propose a novel framework for approximating the explicit MPC law for linear parameter-varying systems using supervised learning. In contrast to most existing approaches, we not only learn the control policy, but also a "certificate policy", that allows us to estimate the sub-optimality of the learned control policy online, during execution-time. We learn both these policies from data using supervised learning techniques, and also provide a randomized method that allows us to guarantee the quality of each learned policy, measured in terms of feasibility and optimality. This in turn allows us to bound the probability of the learned control policy of being infeasible or suboptimal, where the check is performed by the certificate policy. Since our algorithm does not require the solution of an optimization problem during run-time, it can be deployed even on resource-constrained systems. We illustrate the efficacy of the proposed framework on a vehicle dynamics control problem where we demonstrate a speedup of up to two orders of magnitude compared to online optimization with minimal performance degradation.


title: Characterization of a silicon photo-multiplier array with summing board as a photo-multiplier tube replacement in organic scintillator assemblies

abstract: We report on the energy, timing, and pulse-shape discrimination performance of cylindrical 5 cm diameter x 5 cm thick and 7 cm diameter x 7 cm thick {\it trans}-stilbene crystals read out with the passively summed output of three different commercial silicon photo-multiplier arrays. Our results indicate that using the summed output of an 8x8 array of SiPMs provides performance competitive with photo-multiplier tubes for many neutron imaging and correlated particle measurements: for the 5x5 cm crystal read out with SensL's ArrayJ-60035\_64P-PCB, which had the best overall properties, we measure the energy resolution as 13.6$\pm$1.8\% at 341 keVee, the timing resolution in the 100--400 keVee range as 277$\pm$34 ps, and the pulse-shape discrimination figure-of-merit as 2.21$\pm$0.03 in the 230--260 keVee energy range. These results enable many scintillator-based instruments to enjoy the size, robustness, and power benefits of silicon photo-multiplier arrays as replacement for the photo-multiplier tubes that are predominantly used today.