title: A simple generalization of Prandtl-Tomlinson model to study nanoscale rolling friction

abstract: Prandtl-Tomlinson (PT) model has been very successful in explaining nanoscale friction in a variety of situations. However, the simplistic PT model, on account of having a point mass being dragged across a sinusoidal force field, cannot be used for studying rolling friction at nanoscales. In this manuscript, we generalize the PT model as a collection of point particles arranged in a circle of radius $R$. The resulting ``rigid body'' is driven in a composite force field by a moving spring (of stiffness $k$) connected to the center of mass of the rigid body in presence of damping. The force field is a product of the familiar sinusoidal function used in the PT model with a parametrically controlled ($\lambda$) exponentially varying function that is dependent on the vertical coordinates of the particles. Our generalized model degenerates to the standard PT model if $R \ll 1$ and $\lambda \to 0$. With $R \sim 1$ and $\lambda \to 0$, the model undergoes a transition from sticky dynamics to smooth dynamics as $k$ is increased to a critical value. The analytical expression agrees well with the simulation results. Similar analytical expressions have been derived for $ \lambda \neq 0$ as well. In this scenario, the sticky dynamics is experienced in both $x$ and $y$ directions, and our numerical results agree with the analytical solution for $x$ direction. The dynamics, investigated numerically for the general case of $R \sim 1$ and $\lambda \neq 0$, reveals several interesting aspects of nanoscale tribology including the regimes where energy dissipation due to friction is minimum. Further the results from our proposed model are in qualitative agreement with those from MD simulations as well. We believe that the simplicity of our model along with its similarity to the PT model may make it a popular tool for analyzing complicated nanotribological regimes.


title: Convolutional neural networks model improvements using demographics and image processing filters on chest x-rays

abstract: Purpose: The purpose of this study was to observe change in accuracies of convolutional neural networks (CNN) models (ratio of correct classifications to total predictions) on thoracic radiological images by creating different binary classification models based on age, gender, and image pre-processing filters on 14 pathologies. Methodology: This is a quantitative research exploring variation in CNN model accuracies. Radiological thoracic images were divided by age and gender and pre-processed by various image processing filters. Findings: We found partial support for enhancement to model accuracies by segregating modeling images by age and gender and applying image processing filters even though image processing filters are sometimes thought of as information filters. Research limitations: This study may be biased because it is based on radiological images by another research that tagged the images using an automated process that was not checked by a human. Practical implications: Researchers may want to focus on creating models segregated by demographics and pre-process the modeling images using image processing filters. Practitioners developing assistive technologies for thoracic diagnoses may benefit from incorporating demographics and employing multiple models simultaneously with varying statistical likelihood. Originality/value: This study uses demographics in model creation and utilizes image processing filters to improve model performance. Keywords: Convolutional Neural Network (CNN), Chest X-Ray, ChestX-ray14, Lung, Atelectasis, Cardiomegaly, Consolidation, Edema, Effusion, Emphysema, Infiltration, Mass, Nodule, Pleural Thickening, Pneumonia, Pneumathorax


title: Fuzzy approach on modelling cyber attacks patterns on data transfer in industrial control systems

abstract: Cybersecurity of industrial control system is a very complex and challenging research topic, due to the integration of these systems in national critical infrastructures. The control systems are now interconnected in industrial networks and frequently to the Internet. In this context they are becoming targets of various cyber attacks conducted by malicious people such as hackers, script kiddies, industrial spies and even foreign armies and intelligence agencies. In this paper the authors propose a way to model the most frequent attacker profiles and to estimate the success rate of an attack conducted in given conditions. The authors use a fuzzy approach for generating attacker profiles based on attacker attributes such as knowledge, technical resources and motivation. The attack success rate is obtained by using another fuzzy inference system that analyzes the attacker profile and system intrinsic characteristics.


title: Propagation and stability of relativistic jets

abstract: A simple look at the steady high-energy Universe reveals a clear correlation with outflows generated around compact objects (winds and jets). In the case of relativistic jets, they are thought to be produced as a consequence of the extraction of rotational energy from a Kerr black hole (Blandford-Znajek), or from the disc (Blandford-Payne). A fraction of the large energy budget provided by accretion and/or black hole rotational energy is invested into jet formation. After formation, the acceleration and collimation of these outflows allow them to propagate to large distances away from the compact object. The synchrotron cooling times demand that re-acceleration of particles takes place along the jets to explain high-energy and very-high-energy emission from kiloparsec scales. At these scales, jets in radio galaxies are divided in two main morphological/luminosity types, namely, Fanaroff-Riley type I and II (FRI, FRII), the latter being more luminous, collimated and edge-brightened than the former, which show clear hints of decollimation and deceleration. In this contribution, I summarise a set of mechanisms that may contribute to dissipate magnetic and kinetic energy: Magnetohydrodynamic instabilities or jet-obstacle interactions trigger shocks, shearing and mixing, which are plausible scenarios for particle acceleration. I also derive an expression for the expected distance in which the entrainment by stellar winds starts to be relevant, which is applicable to FRI jets. Finally, I discuss the differences in the evolutionary scenarios and the main dissipative mechanisms that take place in extragalactic and microquasar jets.


title: A simple explanation for the "shuffling phenomenon'' for lozenge tilings of dented hexagons

abstract: In a recent paper, Lai and Rohatgi proved a "shuffling theorem" for lozenge tilings of a hexagon with "dents" (i.e., missing triangles). Here, we shall point out that this follows immediately from the enumeration of Gelfand--Tsetlin patterns with given bottom row. This observation is also contained in a recent preprint of Byun.


title: Introduction to Computational Thinking: a new high school curriculum using CodeWorld

abstract: The Louisiana Department of Education partnered with the Gordon A. Cain Center at LSU to pilot a Computing High School Graduation Pathway. The first course in the pathway, Introduction to Computational Thinking (ICT), is designed to teach programming and reinforce mathematical practice skills of nine-grade students, with an emphasis on promoting higher order thinking. In 2017-18, about 200 students and five teachers participated in the pilot, in 2018-2019 the participation increased to 400 students, and in the current 2019-2020 year about 800 students in 11 schools are involved. Professional development starts with a five-week intensive summer institute, which is complemented with follow-up Saturday sessions and coaching support during the academic year. After describing the course content and briefly the teacher training, we discuss the data we have collected in the last two years. The overall student reception of the course has been positive, but the course was categorized by most students as hard. However, the Computing Attitude Survey analysis indicates that the difficulty of the course did not demotivate the students. The pre-post test content assessments show that students learned not only the language, but also general principles of programming, logic and modeling, as well as use of variables, expressions and functions. Lessons learned during the pilot phase motivated changes, such as emphasizing during PD the need to provide timely feedback to students, provide detailed rubrics for the projects and reorganize the lessons to increase the initial engagement with the material. After two years of running pilots, the course is becoming student-centered, where most of the code and image samples provided in the lessons are based on code created by previous students.


title: Visualizing structural balance in signed networks

abstract: Network visualization has established as a key complement to network analysis since the large variety of existing network layouts are able to graphically highlight different properties of networks. However, signed networks, i.e., networks whose edges are labeled as friendly (positive) or antagonistic (negative), are target of few of such layouts and none, to our knowledge, is able to show structural balance, i.e., the tendency of cycles towards including an even number of negative edges, which is a well-known theory for studying friction and polarization. In this work we present Structural-balance-viz: a novel visualization method showing whether a connected signed network is balanced or not and, in the latter case, how close the network is to be balanced. Structural-balance-viz exploits spectral computations of the signed Laplacian matrix to place network's nodes in a Cartesian coordinate system resembling a balance (a scale). Moreover, it uses edge coloring and bundling to distinguish positive and negative interactions. The proposed visualization method has characteristics desirable in a variety of network analysis tasks: Structural-balance-viz is able to provide indications of balance/polarization of the whole network and of each node, to identify two factions of nodes on the basis of their polarization, and to show their cumulative characteristics. Moreover, the layout is reproducible and easy to compare. Structural-balance-viz is validated over synthetic-generated networks and applied to a real-world dataset about political debates confirming that it is able to provide meaningful interpretations.


title: Neural language modeling of free word order argument structure

abstract: Neural language models trained with a predictive or masked objective have proven successful at capturing short and long distance syntactic dependencies. Here, we focus on verb argument structure in German, which has the interesting property that verb arguments may appear in a relatively free order in subordinate clauses. Therefore, checking that the verb argument structure is correct cannot be done in a strictly sequential fashion, but rather requires to keep track of the arguments' cases irrespective of their orders. We introduce a new probing methodology based on minimal variation sets and show that both Transformers and LSTM achieve a score substantially better than chance on this test. As humans, they also show graded judgments preferring canonical word orders and plausible case assignments. However, we also found unexpected discrepancies in the strength of these effects, the LSTMs having difficulties rejecting ungrammatical sentences containing frequent argument structure types (double nominatives), and the Transformers tending to overgeneralize, accepting some infrequent word orders or implausible sentences that humans barely accept.


title: Quantum gravity on a torus

abstract: Causal Dynamical Triangulations (CDT) is a non-perturbative lattice approach to quantum gravity where one assumes space-time foliation into spatial hyper-surfaces of fixed topology. Most of the CDT results were obtained for the spatial topology of the 3-sphere. It was shown that CDT has rich phase structure, including the semiclassical phase consistent with Einstein's general relativity. Some of the phase transitions were found to be second (or higher) order which makes a possibility of taking continuum limit viable. Here we present new results of changing the spatial topology to that of the 3-torus. We argue that the topology change does not change the phase structure nor the order of the phase transitions. Therefore CDT results seem to be universal independent of the topology chosen.


title: Heuristic derivation of Casimir effect in minimal length theories

abstract: We propose a heuristic derivation of Casimir effect in the context of minimal length theories based on a Generalized Uncertainty Principle (GUP). By considering a GUP with only a quadratic term in the momentum, we compute corrections to the standard formula of Casimir energy for the parallel-plate geometry, the sphere and the cylindrical shell. For the first configuration, we show that our result is consistent with the one obtained via more rigorous calculations in Quantum Field Theory. Experimental developments are finally discussed.


title: $H_0$ tension and the String Swampland

abstract: [Abridged] We realize the Agrawal-Obied-Vafa (AOV) swampland proposal of fading dark matter by the model of Salam-Sezgin and its string realization of Cvetic-Gibbons-Pope. The model describes a compactification of 6-dimensional supergravity with a monopole background on a 2-sphere. In 4 dimensions, there are 2 scalar fields, $X$ and $Y $, and the effective potential in the Einstein frame is an exponential with respect to $Y$ times a quadratic polynomial in the field $e^{-X}$. When making the volume of the 2-sphere large, namely for large values of $Y$, there appears a tower of states, which according to the infinite distance swampland conjecture becomes exponentially massless. If the standard model fields are confined on Neveu-Schwarz 5-branes the 6-dimensional gauge couplings are independent of the string dilaton in the string frame, and upon compactification to 4 dimensions the 4-dimensional gauge couplings depend on $X$ (rather than the dilaton $Y$) which is fixed at the minimum of the potential. This avoids direct couplings of the dilaton to matter suppressing extra forces competing with gravity. We show that this set up has the salient features of the AOV models, and ergo can potentially ameliorate the tension between local distance ladder and cosmic microwave background estimates of the Hubble constant $H_0$. Although the AOV proposal does not fully resolve the tension in $H_0$ measurements, it provides a dynamical dark energy model of cosmology that satisfies the de Sitter swampland conjecture. We comment on a viable solution to overcome the tension between low- and high-redshift observations within the AOV background and discuss the implications for the swampland program.


title: The Polarization of X-rays from Warped Black Hole Accretion Disks

abstract: It is commonly assumed that in black hole accretion disks the angular momenta of the disk and the black hole are aligned. However, for a significant fraction of stellar mass black holes and supermassive black holes, the momenta may not be aligned. In such systems, the interplay of disk viscosity and general relativistic frame dragging can cause the disk to warp or break into two (or more) distinct planes; this is called the Bardeen-Petterson effect. We have developed a general relativistic ray-tracing code to find the energy spectra and polarization of warped accretion disks, accounting for the emission from the disk and for photons reflecting one or multiple times off the warped accretion disk segments. We find that polarization angle can be used to give a lower limit on the misalignment angle when a previous measurement of the jet, which is thought be aligned with the black hole angular momentum, can be spatially resolved.


title: A Machine Learning Approach to Adaptive Robust Utility Maximization and Hedging

abstract: We investigate the adaptive robust control framework for portfolio optimization and loss-based hedging under drift and volatility uncertainty. Adaptive robust problems offer many advantages but require handling a double optimization problem (infimum over market measures, supremum over the control) at each instance. Moreover, the underlying Bellman equations are intrinsically multi-dimensional. We propose a novel machine learning approach that solves for the local saddle-point at a chosen set of inputs and then uses a nonparametric (Gaussian process) regression to obtain a functional representation of the value function. Our algorithm resembles control randomization and regression Monte Carlo techniques but also brings multiple innovations, including adaptive experimental design, separate surrogates for optimal control and the local worst-case measure, and computational speed-ups for the sup-inf optimization. Thanks to the new scheme we are able to consider settings that have been previously computationally intractable and provide several new financial insights about learning and optimal trading under unknown market parameters. In particular, we demonstrate the financial advantages of adaptive robust framework compared to adaptive and static robust alternatives.


title: Scalable Graph Algorithms

abstract: Processing large complex networks recently attracted considerable interest. Complex graphs are useful in a wide range of applications from technological networks to biological systems like the human brain. Sometimes these networks are composed of billions of entities that give rise to emerging properties and structures. Analyzing these structures aids us in gaining new insights about our surroundings. As huge networks become abundant, there is a need for scalable algorithms to perform analysis. A prominent example is the PageRank algorithm, which is one of the measures used by web search engines such as Google to rank web pages displayed to the user. In order to find these patterns, massive amounts of data have to be acquired and processed. Designing and evaluating scalable graph algorithms to handle these data sets is a crucial task on the road to understanding the underlying systems. This habilitation thesis is a summary a broad spectrum of scalable graph algorithms that I developed over the last six years with many coauthors. In general, this research is based on four pillars: multilevel algorithms, practical kernelization, parallelization and memetic algorithms that are highly interconnected. Experiments conducted indicate that our algorithms find better solutions and/or are much more scalable than the previous state-of-the-art.


title: From 2-d Polyakov Action to the 4-d Pseudo-Conformal Field Theory

abstract: The characteristic property of the 2-dimensional Polyakov action is its independence on the metric tensor, without being topological. A renormalizable 4-dimensional action is found satisfying this fundamental property. The fundamental quantity of this pseudo-conformal field theory (PCFT) is the lorentzian Cauchy-Riemann (LCR) structure. This action describes all current phenomenology: 1) The Poincar\'e group is determined. 2) Stable solitonic LCR-tetrads are found, which belong to representations of the Poincar\'e group and they are determined by the irreducible and reducible algebraic quadratic surfaces of CP3. 3) The static (irreducible) LCR-structure implies the Kerr-Newman manifold with g=2 gyromagnetic ratio and it is identified with the electron. The stationary (reducible) LCR-structure is identified with the neutrino. The antiparticles have conjugate LCR-structures. The Hawking-Penrose singularity theorems are bypassed in the electron LCR-manifold. 4) The LCR-tetrad defines Einstein's metric and the U(2) electroweak connection. 5) An effective leptonic standard model action is derived using the Bogoliubov-Scharf recursive procedure. 6) The three generations of flavors are implied by the limited number (for curved spacetime) of permitted algebraic surfaces of CP3. 7) For every LCR-structure there exists a solitonic distributional gauge field configuration, identified with the corresponding quark, which explains the lepton-quark correspondence. It is explicitly computed for the static LCR-structure. 8) The derivation of a proper geometric SU(3) Cartan connection opens up the possibility to achieve Einstein's goal to derive all interactions from the pure geometric LCR-structure.


title: Harmonic measure of the outer boundary of colander sets

abstract: We present two companion results: Phragm\'en-Lindel\"of type tight bounds on the minimal possible growth of subharmonic functions with recurrent zero set, and tight bounds on the maximal possible decay of the harmonic measure of the outer boundary of colander sets.


title: Hierarchic Controllability for a nonlinear parabolic equation in one dimension

abstract: This paper deals with the hierarchical control of the parabolic equation.We use Stackelberg{Nash strategies. As usual, we consider one leader and two followers. To each leader we associate a Nash equilibrium corresponding to a bi-objective optimal control problem, then, we look for a leader that solves null controllability e with trajectories problem. We consider linear and nonlinear equations in dimension 1.


title: High-temperature expansion of the grand thermodynamic potential for scalar particles in crossed electromagnetic fields

abstract: The problem of a scalar particle in a constant crossed electromagnetic field ($\mathbf{E}\perp\mathbf{H}$ and $|\mathbf{E}|=|\mathbf{H}|$) is examined. The high-temperature expansion of the grand thermodynamic potential and vacuum energy with account for non-perturbative corrections are derived. The contribution from particles and antiparticles is considered separately. It is shown that the non-perturbative corrections depend on boundary conditions but do not depend on fields.


title: Linear Stability in the Inner Heliosphere: Helios Re\"evaluated

abstract: Wave-particle instabilities driven by departures from local thermodynamic equilibrium have been conjectured to play a role in governing solar wind dynamics. We calculate the statistical variation of linear stability over a large subset of Helios I and II fast solar wind observations using a numerical evaluation of the Nyquist stability criterion, accounting for multiple sources of free energy associated with protons and helium including temperature anisotropies and relative drifts. We find that 88\% of the surveyed intervals are linearly unstable. The median growth rate of the unstable modes is within an order of magnitude of the turbulent transfer rate, fast enough to potentially impact the turbulent scale-to-scale energy transfer. This rate does not significantly change with radial distance, though the nature of the unstable modes, and which ion components are responsible for driving the instabilities, does vary. The effect of ion-ion collisions on stability is found to be significant; collisionally young wind is much more unstable than collsionally old wind, with very different kinds of instabilities present in the two kinds of wind.


title: Smooth Fictitious Play in $N\times 2$ Potential Games

abstract: The paper shows that smooth fictitious play converges to a neighborhood of a pure-strategy Nash equilibrium with probability 1 in almost all $N\times 2$ ($N$-player, two-action) potential games. The neighborhood of convergence may be made arbitrarily small by taking the smoothing parameter to zero. Simple proof techniques are furnished by considering regular potential games.


title: Emergence of extended Newtonian gravity from thermodynamics

abstract: Discovery of a novel thermodynamic aspect of nonrelativistic gravity is reported. Here, initially, an unspecified scalar field potential is considered and treated not as an externally applied field but as a thermodynamic variable on an equal footing with the fluid variables. It is shown that the second law of thermodynamics imposes a stringent constraint on the field, and, quite remarkably, the allowable field turns out to be only of gravity. The resulting field equation for the gravitational potential derived from the analysis of the entropy production rate contains a dissipative term due to irreversibility. It is found that the system relaxes to the conventional theory of Newtonian gravity up to a certain spatial scale, whereas on the larger scale there emerges non-Newtonian gravity described by a nonlinear field equation containing a single coefficient. A comment is made on an estimation of the coefficient that has its origin in the thermodynamic property of the system.

categories: cond-mat.stat-mech astro-ph.GA gr-qc physics.class-ph


title: Idle Time Optimization for Target Assignment and Path Finding in Sortation Centers

abstract: In this paper, we study the one-shot and lifelong versions of the Target Assignment and Path Finding problem in automated sortation centers, where each agent needs to constantly assign itself a sorting station, move to its assigned station without colliding with obstacles or other agents, wait in the queue of that station to obtain a parcel for delivery, and then deliver the parcel to a sorting bin. The throughput of such centers is largely determined by the total idle time of all stations since their queues can frequently become empty. To address this problem, we first formalize and study the one-shot version that assigns stations to a set of agents and finds collision-free paths for the agents to their assigned stations. We present efficient algorithms for this task based on a novel min-cost max-flow formulation that minimizes the total idle time of all stations in a fixed time window. We then demonstrate how our algorithms for solving the one-shot problem can be applied to solving the lifelong problem as well. Experimentally, we believe to be the first researchers to consider real-world automated sortation centers using an industrial simulator with realistic data and a kinodynamic model of real robots. On this simulator, we showcase the benefits of our algorithms by demonstrating their efficiency and effectiveness for up to 350 agents.


title: Averaging Essential and Fundamental Matrices in Collinear Camera Settings

abstract: Global methods to Structure from Motion have gained popularity in recent years. A significant drawback of global methods is their sensitivity to collinear camera settings. In this paper, we introduce an analysis and algorithms for averaging bifocal tensors (essential or fundamental matrices) when either subsets or all of the camera centers are collinear. We provide a complete spectral characterization of bifocal tensors in collinear scenarios and further propose two averaging algorithms. The first algorithm uses rank constrained minimization to recover camera matrices in fully collinear settings. The second algorithm enriches the set of possibly mixed collinear and non-collinear cameras with additional, "virtual cameras," which are placed in general position, enabling the application of existing averaging methods to the enriched set of bifocal tensors. Our algorithms are shown to achieve state of the art results on various benchmarks that include autonomous car datasets and unordered image collections in both calibrated and unclibrated settings.


title: Do Metal-Rich Stars Make Metal-Rich Planets? New Insights on Giant Planet Formation from Host Star Abundances

abstract: The relationship between the compositions of giant planets and their host stars is of fundamental interest in understanding planet formation. The solar system giant planets are enhanced above solar composition in metals, both in their visible atmospheres and bulk compositions. A key question is whether the metal enrichment of giant exoplanets is correlated with that of their host stars. Thorngren et al. (2016) showed that in cool (Teq < 1000 K) giant exoplanets, the total heavy-element mass increases with total Mp and the heavy element enrichment relative to the parent star decreases with total Mp. In their work, the host star metallicity was derived from literature [Fe/H] measurements. Here we conduct a more detailed and uniform study to determine whether different host star metals (C, O, Mg, Si, Fe, and Ni) correlate with the bulk metallicity of their planets, using correlation tests and Bayesian linear fits. We present new host star abundances of 19 cool giant planet systems, and combine these with existing host star data for a total of 22 cool giant planet systems (24 planets). Surprisingly, we find no clear correlation between stellar metallicity and planetary residual metallicity (the relative amount of metal versus that expected from the planet mass alone), which is in conflict with common predictions from formation models. We also find a potential correlation between residual planet metals and stellar volatile-to-refractory element ratios. These results provide intriguing new relationships between giant planet and host star compositions for future modeling studies of planet formation.


title: Rational Polyhedral Outer-Approximations of the Second-Order Cone

abstract: It is well-known that the second-order cone can be outer-approximated to an arbitrary accuracy $\epsilon$ by a polyhedral cone of compact size defined by irrational data. In this paper, we propose two rational polyhedral outer-approximations of compact size retaining the same guaranteed accuracy $\epsilon$. The first outer-approximation has the same size as the optimal but irrational outer-approximation from the literature. In this case,we provide a practical approach to obtain such an approximation defined by the smallest integer coefficients possible, which requires solving a few, small-size integer quadratic programs. The second outer-approximation has a size larger than the optimal irrational outer-approximation by a linear additive factor in the dimension of the second-order cone. However, in this case, the construction is explicit, and it is possible to derive an upper bound on the largest coefficient, which is sublinear in $\epsilon$ and logarithmic in the dimension. We also propose a third outer-approximation, which yields the best possible approximation accuracy given an upper bound on the size of its coefficients. Finally, we discuss two theoretical applications in which having a rational polyhedral outer-approximation is crucial, and run some experiments which explore the benefits of the formulations proposed in this paper from a computational perspective.


title: Stationary polyhedral varifolds minimize area

abstract: We prove that every stationary polyhedral varifold minimizes area in the following senses: (1) its area cannot be decreased by a one-to-one Lipschitz ambient deformation that coincides with the identity outside of a compact set, and (2) it is the varifold associated to a mass-minimizing flat chain with coefficients in a certain metric abelian group. NOTE: After this paper was posted, I learned that (1) and (2) were already proved by Choe and Morgan, respectively. Thus this paper is an exposition of their results.


title: Measuring out-of-time-ordered correlation functions with a single impurity qubit in a bosonic Josephson junction

abstract: We calculate the out-of-time-ordered correlation function (OTOC) of a single impurity qubit coupled to fully a connected many-particle system such as a bosonic Josephson junction or spins with long-range interactions. In these systems the qubit OTOC can be used to detect both ground state and excited state quantum phase transitions (QPTs), making it a robust order parameter that is considerably more sensitive than the standard one-body correlation function. Finite size scaling exponents for an $N$ body system can also be accurately extracted from the long-time OTOC dynamics, however, for short times there is a discrepancy due to the fact that the qubit has not had enough time to couple to the larger system. Our results show that the OTOC of even the smallest probe is enough to diagnose a QPT in fully connected models but, like a continuous measurement, can still cause a backaction effect which leads to weakly chaotic dynamics and gradual information scrambling.


title: Asymptotic Mean Value Laplacian in Metric Measure Spaces

abstract: We use the mean value property in an asymptotic way to provide a notion of a pointwise Laplacian, called AMV Laplacian, that we study in several contexts including the Heisenberg group and weighted Lebesgue measures. We focus especially on a class of metric measure spaces including intersecting submanifolds of $\mathbb{R}^n$, a context in which our notion brings new insights; the Kirchhoff law appears as a special case. In the general case, we also prove a maximum and comparison principle, as well as a Green-type identity for a related operator.


title: Transferable Force-Torque Dynamics Model for Peg-in-hole Task

abstract: We present a learning-based force-torque dynamics to achieve model-based control for contact-rich peg-in-hole task using force-only inputs. Learning the force-torque dynamics is challenging because of the ambiguity of the low-dimensional 6-d force signal and the requirement of excessive training data. To tackle these problems, we propose a multi-pose force-torque state representation, based on which a dynamics model is learned with the data generated in a sample-efficient offline fashion. In addition, by training the dynamics model with peg-and-holes of various shapes, scales, and elasticities, the model could quickly transfer to new peg-and-holes after a small number of trials. Extensive experiments show that our dynamics model could adapt to unseen peg-and-holes with 70% fewer samples required compared to learning from scratch. Along with the learned dynamics, model predictive control and model-based reinforcement learning policies achieve over 80% insertion success rate. Our video is available at https://youtu.be/ZAqldpVZgm4.


title: The Ooguri-Vafa space as a moduli space of framed wild harmonic bundles

abstract: The Ooguri-Vafa space is a 4-dimensional incomplete hyperk\"{a}hler manifold, defined on the total space of a singular torus fibration with one singular nodal fiber. It has been proposed that the Ooguri-Vafa hyperk\"ahler metric should be part of the local model of the hyperk\"{a}hler metric of the Hitchin moduli spaces, near the most generic kind of singular locus of the Hitchin fibration. In order to relate the Ooguri-Vafa space with the Hitchin moduli spaces, we show that the Ooguri-Vafa space can be interpreted as a set of rank 2, framed wild harmonic bundles over $\mathbb{C}P^1$, with one irregular singularity. Along the way we show that a certain twistor family of holomorphic Darboux coordinates, which describes the hyperk\"{a}hler geometry of the Ooguri-Vafa space, has an interpretation in terms of Stokes data associated to our framed wild harmonic bundles.


title: Image segmentation of liver stage malaria infection with spatial uncertainty sampling

abstract: Global eradication of malaria depends on the development of drugs effective against the silent, yet obligate liver stage of the disease. The gold standard in drug development remains microscopic imaging of liver stage parasites in in vitro cell culture models. Image analysis presents a major bottleneck in this pipeline since the parasite has significant variability in size, shape, and density in these models. As with other highly variable datasets, traditional segmentation models have poor generalizability as they rely on hand-crafted features; thus, manual annotation of liver stage malaria images remains standard. To address this need, we develop a convolutional neural network architecture that utilizes spatial dropout sampling for parasite segmentation and epistemic uncertainty estimation in images of liver stage malaria. Our pipeline produces high-precision segmentations nearly identical to expert annotations, generalizes well on a diverse dataset of liver stage malaria parasites, and promotes independence between learned feature maps to model the uncertainty of generated predictions.


title: Survival analysis for AdVerse events with VarYing follow-up times (SAVVY): Rationale and statistical concept of a meta-analytic study

abstract: The assessment of safety is an important aspect of the evaluation of new therapies in clinical trials, with analyses of adverse events being an essential part of this. Standard methods for the analysis of adverse events such as the incidence proportion, i.e. the number of patients with a specific adverse event out of all patients in the treatment groups, do not account for both varying follow-up times and competing risks. Alternative approaches such as the Aalen-Johansen estimator of the cumulative incidence function have been suggested. Theoretical arguments and numerical evaluations support the application of these more advanced methodology, but as yet there is to our knowledge only insufficient empirical evidence whether these methods would lead to different conclusions in safety evaluations. The Survival analysis for AdVerse events with VarYing follow-up times (SAVVY) project strives to close this gap in evidence by conducting a meta-analytical study to assess the impact of the methodology on the conclusion of the safety assessment empirically. Here we present the rationale and statistical concept of the empirical study conducted as part of the SAVVY project. The statistical methods are presented in unified notation and examples of their implementation in R and SAS are provided.


title: Towards Efficient Integration of Blockchain for IoT Security: The Case Study of IoT Remote Access

abstract: The booming Internet of Things (IoT) market has drawn tremendous interest from cyber attackers. The centralized cloud-based IoT service architecture has serious limitations in terms of security, availability, and scalability, and is subject to single points of failure (SPOF). Recently, accommodating IoT services on blockchains has become a trend for better security, privacy, and reliability. However, blockchain's shortcomings of high cost, low throughput, and long latency make it unsuitable for IoT applications. In this paper, we take a retrospection of existing blockchain-based IoT solutions and propose a framework for efficient blockchain and IoT integration. Following the framework, we design a novel blockchain-assisted decentralized IoT remote accessing system, RS-IoT, which has the advantage of defending IoT devices against zero-day attacks without relying on any trusted third-party. By introducing incentives and penalties enforced by smart contracts, our work enables "an economic approach" to thwarting the majority of attackers who aim to achieve monetary gains. Our work presents an example of how blockchain can be used to ensure the fairness of service trading in a decentralized environment and punish misbehaviors objectively. We show the security of RS-IoT via detailed security analyses. Finally, we demonstrate its scalability, efficiency, and usability through a proof-of-concept implementation on the Ethereum testnet blockchain.


title: Data-driven molecular modeling with the generalized Langevin equation

abstract: The complexity of molecular dynamics simulations necessitates dimension reduction and coarse-graining techniques to enable tractable computation. The generalized Langevin equation (GLE) describes coarse-grained dynamics in reduced dimensions. In spite of playing a crucial role in non-equilibrium dynamics, the memory kernel of the GLE is often ignored because it is difficult to characterize and expensive to solve. To address these issues, we construct a data-driven rational approximation to the GLE. Building upon previous work leveraging the GLE to simulate simple systems, we extend these results to more complex molecules, whose many degrees of freedom and complicated dynamics require approximation methods. We demonstrate the effectiveness of our approximation by testing it against exact methods and comparing observables such as autocorrelation and transition rates.


title: Superconducting proximity effect in flat band systems

abstract: We study theoretically proximity-induced superconductivity and its inverse effect in dice lattice flat band model by considering Josephson junction with an s-wave pairing in the superconducting leads. Using self-consistent tight-binding Bogoliubov-de Gennes method, we show that there is a critical value for chemical potential of the superconductors depending on paring interaction strength over which for undoped normal region the proximity effect is enhanced. Whereas if the superconductor chemical potential is less than the critical one the proximity effect decreases regardless of normal region doping and in the meanwhile, the pairing amplitude of superconducting region increases significantly. Furthermore, we unveil that the supercurrent passing through the junction is large (vanishingly small) when the superconductor chemical potential is smaller (larger) than the critical value which increases as a function of normal region chemical potential.


title: Fluid models with phase transition for kinetic equations in swarming

abstract: We concentrate on kinetic models for swarming with individuals interacting through self-propelling and friction forces, alignment and noise. We assume that the velocity of each individual relaxes to the mean velocity. In our present case, the equilibria depend on the density and the orientation of the mean velocity, whereas the mean speed is not anymore a free parameter and a phase transition occurs in the homogeneous kinetic equation. We analyze the profile of equilibria for general potentials identifying a family of potentials leading to phase transitions. Finally, we derive the fluid equations when the interaction frequency becomes very large.


title: WLS-ENO Remap: Superconvergent and Non-Oscillatory Weighted Least Squares Data Transfer on Surfaces

abstract: Data remap between non-matching meshes is a critical step in multiphysics coupling using a partitioned approach. The data fields being transferred often have jumps in function values or derivatives. It is important but very challenging to avoid spurious oscillations (a.k.a. the Gibbs Phenomenon) near discontinuities and at the same time to achieve high-order accuracy away from discontinuities. In this work, we introduce a new approach, called WLS-ENOR, or Weighted-Least-Squares-based Essentially Non-Oscillatory Remap, to address this challenge. Based on the WLS-ENO reconstruction technique proposed by Liu and Jiao (J. Comput. Phys. vol 314, pp 749--773, 2016), WLS-ENOR differs from WLS-ENO and other WENO schemes in that it resolves not only the O(1) oscillations due to C 0 discontinuities, but also the accumulated effect of O(h) oscillations due to C 1 discontinuities. To this end, WLS-ENOR introduces a robust detector of discontinuities and a new weighting scheme for WLS-ENO near discontinuities. We also optimize the weights at smooth regions to achieve superconvergence. As a result, WLS-ENOR is more than fifth-order accurate and highly conservative in smooth regions, while being non-oscillatory and minimally diffusive near discontinuities. We also compare WLS-ENOR with some commonly used methods based on L 2 projection, moving least squares, and radial basis functions.


title: Optimal forest rotation under carbon pricing and forest damage risk

abstract: Forests will have two notable economic roles in the future: providing renewable raw material and storing carbon to mitigate climate change. The pricing of forest carbon leads to longer rotation times and consequently larger carbon stocks, but also exposes landowners to a greater risk of forest damage. This paper investigates optimal forest rotation under carbon pricing and forest damage risk. I provide the optimality conditions for this problem and illustrate the setting with numerical calculations representing boreal forests under a range of carbon prices and damage probabilities. The relation between damage probability and carbon price towards the optimal rotation length is nearly linear, with carbon pricing having far greater impact. As such, increasing forest carbon stocks by lengthening rotations is an economically attractive method for climate change mitigation, despite the forest damage risk. Carbon pricing also increases land expectation value and reduces the economic risks of the landowner. The production possibility frontier under optimal rotation suggests that significantly larger forests carbon stocks are achievable, but imply lower harvests. However, forests' societally optimal role between these two activities is not yet clear-cut; but rests on the future development of relative prices between timber, carbon and other commodities dependent on land-use.


title: Challenges for automated spike sorting: beware of pharmacological manipulations

abstract: The advent of large-scale and high-density extracellular recording devices allows simultaneous recording from thousands of neurons. However, the complexity and size of the data makes it mandatory to develop robust algorithms for fully automated spike sorting. Here it is shown that limitations imposed by biological constraints such as changes in spike waveforms induced under different drug regimes should be carefully taken into consideration in future developments.


title: Biometrics Recognition Using Deep Learning: A Survey

abstract: Deep learning-based models have been very successful in achieving state-of-the-art results in many of the computer vision, speech recognition, and natural language processing tasks in the last few years. These models seem a natural fit for handling the ever-increasing scale of biometric recognition problems, from cellphone authentication to airport security systems. Deep learning-based models have increasingly been leveraged to improve the accuracy of different biometric recognition systems in recent years. In this work, we provide a comprehensive survey of more than 120 promising works on biometric recognition (including face, fingerprint, iris, palmprint, ear, voice, signature, and gait recognition), which deploy deep learning models, and show their strengths and potentials in different applications. For each biometric, we first introduce the available datasets that are widely used in the literature and their characteristics. We will then talk about several promising deep learning works developed for that biometric, and show their performance on popular public benchmarks. We will also discuss some of the main challenges while using these models for biometric recognition, and possible future directions to which research in this area is headed.


title: On Multi-Cascade Influence Maximization: Model, Hardness and Algorithmic Framework

abstract: This paper studies the multi-cascade influence maximization problem, which explores strategies for launching one information cascade in a social network with multiple existing cascades. With natural extensions to the classic models, we first propose the independent multi-cascade model where the diffusion process is governed by the so-called activation function. We show that the proposed model is sufficiently flexible as it generalizes most of the existing cascade-based models. We then study the multi-cascade influence maximization problem under the designed model and provide approximation hardness under common complexity assumptions, namely Exponential Time Hypothesis and $NP \subseteq DTIME(n^{\poly \log n})$. Given the hardness results, we build a framework for designing heuristic seed selection algorithms with a testable data-dependent approximation ratio. The designed algorithm leverages upper and lower bounds, which reveal the key combinatorial structure behind the multi-cascade influence maximization problem. The performance of the framework is theoretically analyzed and practically evaluated through extensive simulations. The superiority of the proposed solution is supported by encouraging experimental results, in terms of effectiveness and efficiency.


title: Extended Nestohedra and their Face Numbers

abstract: Nestohedra are a family of convex polytopes that includes permutohedra, associahedra, and graph associahedra. In this paper, we study an extension of such polytopes, called extended nestohedra. We show that these objects are indeed the boundaries of simple polytopes, answering a question of Lam and Pylyavskyy. We also study the duals of (extended) nestohedra, giving a complete characterization of isomorphisms (as simplicial complexes) between the duals of extended nestohedra and a partial characterization of isomorphisms between the duals of nestohedra and extended nestohedra. In addition, we give formulas for their $f$-, $h$-, and $\gamma$-vectors. This includes showing that the $f$-vectors of the extended nestohedron corresponding to a forest $F$ and the nestohedron corresponding to the line graph of $F$ are the same, as well as showing that all flag extended nestohedra have nonnegative $\gamma$-vectors, thus proving Gal's conjecture for a large class of flag simple polytopes. We also relate the $f$- and $h$-vectors of the nestohedra and extended nestohedra, as well as give explicit formulas for the $h$- and $\gamma$-vectors in terms of descent statistics for a certain class of flag extended nestohedra. Finally, we define a partial ordering on partial permutations that is a join semilattice quotient of the weak Bruhar order on the symmetric group, and such that any linear extension of the partial order provides a shelling of the dual of the stellohedron.


title: Fast Neutrino Flavor Instability in the Neutron-star Convection Layer of Three-dimensional Supernova Models

abstract: Neutrinos from a supernova (SN) might undergo fast flavor conversions near the collapsed stellar core. We perform a detailed study of this intriguing possibility, analyzing time-dependent state-of-the-art 3D SN models of 9 and 20 Msun. Both models were computed with multi-D three-flavor neutrino transport based on a two-moment solver, and both exhibit the presence of the lepton-number emission self-sustained asymmetry (LESA). The transport solution does not provide the angular distributions of the neutrino fluxes, which are crucial to track the fast flavor instability. To overcome this limitation, we use a recently proposed approach based on the angular moments of the energy-integrated electron lepton-number distribution. With this method we find the possibility of fast neutrino flavor instability at radii <~20 km, which is well interior to the neutrinosphere. Our results confirm recent observations in a 2D SN model and in 2D/3D models with fixed matter background, which were computed with Boltzmann neutrino transport. However, the flavor unstable locations are not isolated points as discussed previously, but thin skins surrounding volumes where electron antineutrinos are more abundant than electron neutrinos. These volumes grow with time and appear first in the convective layer of the proto-neutron star (PNS), where a decreasing electron fraction (Ye) and high temperatures favor the occurrence of regions with negative neutrino chemical potential. Since Ye remains higher in the LESA dipole direction, where convective lepton-number transport out from the nonconvective PNS core slows down the deleptonization, flavor unstable conditions become more widespread in the opposite hemisphere. This interesting phenomenon deserves further investigation, since its impact on SN modeling and possible consequences for SN dynamics and neutrino observations are presently unclear. (abridged)


title: On the Graph Laplacian and the Rankability of Data

abstract: Recently, Anderson et al. (2019) proposed the concept of rankability, which refers to a dataset's inherent ability to produce a meaningful ranking of its items. In the same paper, they proposed a rankability measure that is based on a integer program for computing the minimum number of edge changes made to a directed graph in order to obtain a complete dominance graph, i.e., an acyclic tournament graph. In this article, we prove a spectral-degree characterization of complete dominance graphs and apply this characterization to produce a new measure of rankability that is cost-effective and more widely applicable. We support the details of our algorithm with several results regarding the conditioning of the Laplacian spectrum of complete dominance graphs and the Hausdorff distance between their Laplacian spectrum and that of an arbitrary directed graph with weights between zero and one. Finally, we analyze the rankability of datasets from the world of chess and college football.


title: Turbulence, String Theory and Ising Model

abstract: We advance the vortex cell approach to turbulence \cite{TSVS} by elaborating the Clebsch field dynamics on the surface of vortex cells. We argue that resulting statistical system can be described as 3D Ising model interacting with Compactified Bosonic String on the Riemann surface describing phase boundary. The cells come out as $\sigma=1$ phase of the Ising model in strong magnetic field. At large circulation around large loop the dominant Ising configuration correspond to minimal volume cells covering the area inside the loop, which naturally leads to the Area law and exponential PDF tails.


title: Branching epistemic uncertainty and thickness of tails

abstract: This is an epistemological approach to errors in both inference and risk management, leading to necessary structural properties for the probability distribution. Many mechanisms have been used to show the emergence of fat tails. Here we follow an alternative route, the epistemological one, using counterfactual analysis, and show how nested uncertainty, that is, errors on the error in estimation of parameters, lead to fattailedness of the distribution. The results have relevant implications for forecasting, dealing with model risk and generally all statistical analyses. The more interesting results are as follows: + The forecasting paradox: The future is fatter tailed than the past. Further, out of sample results should be fatter tailed than in-sample ones. + Errors on errors can be explosive or implosive with different consequences. Infinite recursions can be easily dealt with, pending on the structure of the errors. We also present a method to perform counterfactual analysis without the explosion of branching counterfactuals.


title: R-L Relation in Realistic FRADO Model

abstract: In Failed Radiatively Accelerated Dusty Outflow (FRADO) model which provides the source of material above the accretion disk (AD) as an option to explain the formation mechanism of Broad Line Region (BLR) in AGNs, the BLR inner radius ($\rm{BLR}_{in}$ hereafter) is set by the condition that the dust evaporates immediately upon departure from the AD surface. On the other hand, the location of BLR clouds obtained observationaly via reverberation mapping shows some scaling with the source luminosity, so-called RL relation. We assume $\rm{BLR}_{in}$ to be the location of BLR clouds, then using a realistic expression for the radiation pressure of an AD, and having included the proper values of dust opacity, and shielding effect as well, we report our numerical results on calculation of $\rm{BLR}_{in}$ based on FRADO model. We investigate how it scales with monochromatic luminosity at 5100 angstrom for a grid of blackhole masses and Eddington ratios to compare along with the FRADO analytically predicted RL directly to observational data.


title: Generalized diffusion of quantum Brownian motion

abstract: This article discusses the numerical result predicted by the quantum Langevin equation of the generalized diffusion function of a Brownian particle immersed in an Ohmic quantum bath of harmonic oscillators. The time dependence of the standard deviation of the reduced Wiener function of the system, obtained by integrating the whole function in the momentum space, is determined as well. The complexity of the equations leads to resort to a much simpler calculation based in the position correlation function. They are done for the three possible regimes of the system, namely, periodic, aperiodic and overdamped. It is found in the periodic case that the generalized diffusion is a discontinuous function exhibiting negative values during short time periods of time. This counterintuitive result, found theoretically in other systems and waiting for its experimental confirmation, can be perfectly explained in the framework of the quantum Langevin equation. Its oscillatory behavior is primordially due to the response to the external field while its quantum origin contributes only in its magnitude. The results are compared with those of the continuum limit which exhibits similar behavior.


title: Morphing and Sampling Network for Dense Point Cloud Completion

abstract: 3D point cloud completion, the task of inferring the complete geometric shape from a partial point cloud, has been attracting attention in the community. For acquiring high-fidelity dense point clouds and avoiding uneven distribution, blurred details, or structural loss of existing methods' results, we propose a novel approach to complete the partial point cloud in two stages. Specifically, in the first stage, the approach predicts a complete but coarse-grained point cloud with a collection of parametric surface elements. Then, in the second stage, it merges the coarse-grained prediction with the input point cloud by a novel sampling algorithm. Our method utilizes a joint loss function to guide the distribution of the points. Extensive experiments verify the effectiveness of our method and demonstrate that it outperforms the existing methods in both the Earth Mover's Distance (EMD) and the Chamfer Distance (CD).


title: Many-body wavefunctions for quantum impurities out of equilibrium. I. The nonequilibrium Kondo model

abstract: We present here the details of a method [A. B. Culver and N. Andrei, Phys. Rev. B 103, L201103 (2021)] for calculating the time-dependent many-body wavefunction that follows a local quench. We apply the method to the voltage-driven nonequilibrium Kondo model to find the exact time-evolving wavefunction following a quench where the dot is suddenly attached to the leads at $t=0$. The method, which does not use Bethe ansatz, also works in other quantum impurity models and may be of wider applicability. We show that the long-time limit (with the system size taken to infinity first) of the time-evolving wavefunction of the Kondo model is a current-carrying nonequilibrium steady state that satisfies the Lippmann-Schwinger equation. We show that the electric current in the time-evolving wavefunction is given by a series expression that can be expanded either in weak coupling or in strong coupling, converging to all orders in the steady-state limit in either case. The series agrees to leading order with known results in the well-studied regime of weak antiferromagnetic coupling and also reveals a universal regime of strong ferromagnetic coupling with Kondo temperature $T_K^{(F)} = D e^{-\frac{3\pi^2}{8} \rho |J|}$ ($J<0$, $\rho|J|\to\infty$). In this regime, the differential conductance $dI/dV$ reaches the unitarity limit $2e^2/h$ asymptotically at large voltage or temperature.


title: Counterexamples to Siegel's Conjecture

abstract: We prove that the intersection of a Hirsch polytope and a cube may be a non-Hirsch polytope.


title: Interpreting Deep Learning Features for Myoelectric Control: A Comparison with Handcrafted Features

abstract: The research in myoelectric control systems primarily focuses on extracting discriminative representations from the electromyographic (EMG) signal by designing handcrafted features. Recently, deep learning techniques have been applied to the challenging task of EMG-based gesture recognition. The adoption of these techniques slowly shifts the focus from feature engineering to feature learning. However, the black-box nature of deep learning makes it hard to understand the type of information learned by the network and how it relates to handcrafted features. Additionally, due to the high variability in EMG recordings between participants, deep features tend to generalize poorly across subjects using standard training methods. Consequently, this work introduces a new multi-domain learning algorithm, named ADANN, which significantly enhances (p=0.00004) inter-subject classification accuracy by an average of 19.40% compared to standard training. Using ADANN-generated features, the main contribution of this work is to provide the first topological data analysis of EMG-based gesture recognition for the characterisation of the information encoded within a deep network, using handcrafted features as landmarks. This analysis reveals that handcrafted features and the learned features (in the earlier layers) both try to discriminate between all gestures, but do not encode the same information to do so. Furthermore, using convolutional network visualization techniques reveal that learned features tend to ignore the most activated channel during gesture contraction, which is in stark contrast with the prevalence of handcrafted features designed to capture amplitude information. Overall, this work paves the way for hybrid feature sets by providing a clear guideline of complementary information encoded within learned and handcrafted features.


title: Who Asked Us? How the Theory of Computing Answers Questions about Analysis

abstract: Algorithmic fractal dimensions -- constructs of computability theory -- have recently been used to answer open questions in classical geometric measure theory, questions of mathematical analysis whose statements do not involve computability theory or logic. We survey these developments and the prospects for future such results.


title: Permanental sequences that are related to a Markov chain example of Kolmogorov

abstract: Permanental sequences with non-symmetric kernels that are generalization of the potentials of a Markov chain with state space $\{0,1/2, \ldots, 1/n,\ldots\}$ that was introduced by Kolmogorov, are studied. Depending on a parameter in the kernels we obtain an exact rate of divergence of the sequence at $0$, an exact local modulus of continuity of the sequence at $0$, or a precise bounded discontinuity for the sequence at $0$.


title: Training Distributed Deep Recurrent Neural Networks with Mixed Precision on GPU Clusters

abstract: In this paper, we evaluate training of deep recurrent neural networks with half-precision floats. We implement a distributed, data-parallel, synchronous training algorithm by integrating TensorFlow and CUDA-aware MPI to enable execution across multiple GPU nodes and making use of high-speed interconnects. We introduce a learning rate schedule facilitating neural network convergence at up to $O(100)$ workers. Strong scaling tests performed on clusters of NVIDIA Pascal P100 GPUs show linear runtime and logarithmic communication time scaling for both single and mixed precision training modes. Performance is evaluated on a scientific dataset taken from the Joint European Torus (JET) tokamak, containing multi-modal time series of sensory measurements leading up to deleterious events called plasma disruptions, and the benchmark Large Movie Review Dataset~\cite{imdb}. Half-precision significantly reduces memory and network bandwidth, allowing training of state-of-the-art models with over 70 million trainable parameters while achieving a comparable test set performance as single precision.


title: Nuclear Fission Dynamics: Past, Present, Needs, and Future

abstract: Recent developments in theoretical modeling and in computational power have allowed us to make significant progress on a goal not achieved yet in nuclear theory: a fully microscopic theory of nuclear fission. The complete microscopic description remains a computationally demanding task, but the information that can be provided by current calculations can be extremely useful to guide and constrain phenomenological approaches. First, a truly microscopic framework that can describe the real-time dynamics of the fissioning system can justify or rule out assumptions and approximations incompatible with an accurate quantum treatment or with our understanding of the inter nucleon interactions. Second, the microscopic approach can be used to obtain trends such as: the excitation energy sharing mechanism between fission fragments (FFs) with increasing excitation energy of the fissioning system, the angular momentum content of the FFs, or even to compute observables that cannot be otherwise calculated in phenomenological approaches or even measured, as in the case of astronomical environments. Merely the characterization of the trends would be of great importance for various application. We present here arguments that a truly microscopic approach to fission does not support the assumption of adiabaticity of the large amplitude collective motion in fission, particularly starting from the outer saddle down to the scission configuration.


title: Towards end-to-end verifiable online voting: adding verifiability to established voting systems

abstract: Online voting for independent elections is generally supported by trusted election providers. Typically these providers do not offer any way in which a voter can verify their vote, so the providers are trusted with ballot privacy and ensuring correctness. Despite the desire to offer online voting for political elections, this lack of transparency and verifiability is often seen as a significant barrier to the large-scale adoption of online elections. Adding verifiability to an online election increases transparency and integrity, allowing voters to verify that their vote has been recorded correctly and included in the tally. However, replacing existing online systems with those that provide verifiable voting requires new algorithms and code to be deployed, and this presents a significant business risk to commercial election providers. In this paper we present the first step in an incremental approach which minimises the business risk but demonstrates the advantages of verifiability, by developing an implementation of key elements of a Selene-based verifiability layer and adding it to an operational online voting system. Selene is a verifiable voting protocol that uses trackers to enable voters to confirm that their votes have been captured correctly while protecting voter anonymity. This results in a system where even the election authority running the system cannot change the result in an undetectable way, and gives stronger guarantees on the integrity of the election than were previously present. We explore the challenges presented by adding a verifiability layer to an operational system. We describe the results of two initial trials, which obtained that survey respondents found this form of verifiability easy to use and that they broadly appreciated it. We conclude by outlining the further steps in the road-map towards the deployment of a fully trustworthy online voting system.


title: A Programmatic and Semantic Approach to Explaining and DebuggingNeural Network Based Object Detectors

abstract: Even as deep neural networks have become very effective for tasks in vision and perception, it remains difficult to explain and debug their behavior. In this paper, we present a programmatic and semantic approach to explaining, understanding, and debugging the correct and incorrect behaviors of a neural network-based perception system. Our approach is semantic in that it employs a high-level representation of the distribution of environment scenarios that the detector is intended to work on. It is programmatic in that scenario representation is a program in a domain-specific probabilistic programming language which can be used to generate synthetic data to test a given perception module. Our framework assesses the performance of a perception module to identify correct and incorrect detections, extracts rules from those results that semantically characterizes the correct and incorrect scenarios, and then specializes the probabilistic program with those rules in order to more precisely characterize the scenarios in which the perception module operates correctly or not. We demonstrate our results using the SCENIC probabilistic programming language and a neural network-based object detector. Our experiments show that it is possible to automatically generate compact rules that significantly increase the correct detection rate (or conversely the incorrect detection rate) of the network and can thus help with understanding and debugging its behavior.


title: XGBOD: Improving Supervised Outlier Detection with Unsupervised Representation Learning

abstract: A new semi-supervised ensemble algorithm called XGBOD (Extreme Gradient Boosting Outlier Detection) is proposed, described and demonstrated for the enhanced detection of outliers from normal observations in various practical datasets. The proposed framework combines the strengths of both supervised and unsupervised machine learning methods by creating a hybrid approach that exploits each of their individual performance capabilities in outlier detection. XGBOD uses multiple unsupervised outlier mining algorithms to extract useful representations from the underlying data that augment the predictive capabilities of an embedded supervised classifier on an improved feature space. The novel approach is shown to provide superior performance in comparison to competing individual detectors, the full ensemble and two existing representation learning based algorithms across seven outlier datasets.


title: A sub-Neptune sized planet transiting the M2.5-dwarf G 9-40: Validation with the Habitable-zone Planet Finder

abstract: We validate the discovery of a 2 Earth radii sub-Neptune-size planet around the nearby high proper motion M2.5-dwarf G 9-40 (EPIC 212048748), using high-precision near-infrared (NIR) radial velocity (RV) observations with the Habitable-zone Planet Finder (HPF), precision diffuser-assisted ground-based photometry with a custom narrow-band photometric filter, and adaptive optics imaging. At a distance of $d=27.9\mathrm{pc}$, G 9-40b is the second closest transiting planet discovered by K2 to date. The planet's large transit depth ($\sim$3500ppm), combined with the proximity and brightness of the host star at NIR wavelengths (J=10, K=9.2) makes G 9-40b one of the most favorable sub-Neptune-sized planet orbiting an M-dwarf for transmission spectroscopy with JWST, ARIEL, and the upcoming Extremely Large Telescopes. The star is relatively inactive with a rotation period of $\sim$29 days determined from the K2 photometry. To estimate spectroscopic stellar parameters, we describe our implementation of an empirical spectral matching algorithm using the high-resolution NIR HPF spectra. Using this algorithm, we obtain an effective temperature of $T_{\mathrm{eff}}=3404\pm73$K, and metallicity of $\mathrm{[Fe/H]}=-0.08\pm0.13$. Our RVs, when coupled with the orbital parameters derived from the transit photometry, exclude planet masses above $11.7 M_\oplus$ with 99.7% confidence assuming a circular orbit. From its radius, we predict a mass of $M=5.0^{+3.8}_{-1.9} M_\oplus$ and an RV semi-amplitude of $K=4.1^{+3.1}_{-1.6}\mathrm{m\:s^{-1}}$, making its mass measurable with current RV facilities. We urge further RV follow-up observations to precisely measure its mass, to enable precise transmission spectroscopic measurements in the future.


title: Lattice QCD estimate of the quark-gluon plasma photon emission rate

abstract: We present a computation of the photon emission rate of the quark-gluon plasma from two-flavor lattice QCD at a temperature of 254 MeV, which follows up on the work presented in [1]. We perform a continuum extrapolation of the vector-current correlator, and consider a linear combination of the Lorentz indices corresponding to a UV-finite spectral function. To extract the spectral function from the lattice correlators, an ill-posed inverse problem, we model the spectral function with a Pad\'e ansatz. We further constrain our analysis by simultaneously fitting data with different momenta. We present results for a multi-momentum fit including the three smallest momenta available from our lattice analysis.


title: Energy transfer from large to small scales in turbulence by multi-scale nonlinear strain and vorticity interactions

abstract: An intrinsic feature of turbulent flows is an enhanced rate of mixing and kinetic energy dissipation due to the rapid generation of small-scale motions from large-scale excitation. The transfer of kinetic energy from large to small scales is commonly attributed to the stretching of vorticity by the strain-rate, but strain self-amplification also plays a role. Previous treatments of this connection are phenomenological or inexact, or cannot distinguish the contribution of vorticity stretching from that of strain self-amplification. In this paper, an exact relationship is derived which quantitatively establishes how intuitive multi-scale mechanisms such as vorticity stretching and strain self-amplification together actuate the inter-scale transfer of energy in turbulence. Numerical evidence validates this result and uses it to demonstrate that the contribution of strain self-amplification to energy transfer is higher than that of vorticity stretching, but not overwhelmingly so.


title: Immersion in Rn by Complex irreducible Spinors

abstract: The first time that the connection between isometric immersion of surfaces and solutions of the Dirac equation appeared in the literature was in the seminal paper of Thomas Friedrich in 1998. In consequence of that, several authors contributed to this topic hereafter, by obtaining the spinorial representation of Spin manifolds with arbitrary dimension and also by presenting a generalization of the Weierstrass representation map, for example. All these results assume that the manifolds and bundles involved carry a Spin structure, however this hypothesis is somehow restrictive, as for instance, if we consider complex manifolds, it is more natural to consider SpinC structures. There exist an alternative to adapt this result to SpinC manifolds, where the idea is to use the left regular representation of a Clifford algebra in itself to build the spinor bundles, but unlike the original works of Friedrich and Morel, this representation is not irreducible. Thus, this paper aims to present the spinorial representation of SpinC manifolds into Euclidian space with arbitrary dimensions using spinors that came from an irreducible representation of a complex Clifford algebra.


title: Efficient Estimation of Mixture Cure Frailty Model for Clustered Current Status Data

abstract: Current status data abounds in the field of epidemiology and public health, where the only observable data for a subject is the random inspection time and the event status at inspection. Motivated by such a current status data from a periodontal study where data are inherently clustered, we propose a unified methodology to analyze such complex data. We allow the time-to-event to follow the semiparametric GOR model with a cure fraction, and develop a unified estimation scheme powered by the EM algorithm. The within-subject correlation is accounted for by a random (frailty) effect, and the non-parametric component of the GOR model is approximated via penalized splines, with a set of knot points that increases with the sample size. Proposed methodology is accompanied by a rigorous asymptotic theory, and the related semiparametric efficiency. The finite sample performance of our model parameters are assessed via simulation studies. Furthermore, the proposed methodology is illustrated via application to the oral health data, accompanied by diagnostic checks to identify influential observations. An easy to use R package CRFCSD is also available for implementation.


title: Image Based Identification of Ghanaian Timbers Using the XyloTron: Opportunities, Risks and Challenges

abstract: Computer vision systems for wood identification have the potential to empower both producer and consumer countries to combat illegal logging if they can be deployed effectively in the field. In this work, carried out as part of an active international partnership with the support of UNIDO, we constructed and curated a field-relevant image data set to train a classifier for wood identification of $15$ commercial Ghanaian woods using the XyloTron system. We tested model performance in the laboratory, and then collected real-world field performance data across multiple sites using multiple XyloTron devices. We present efficacies of the trained model in the laboratory and in the field, discuss practical implications and challenges of deploying machine learning wood identification models, and conclude that field testing is a necessary step - and should be considered the gold-standard - for validating computer vision wood identification systems.


title: Nonstandard approach to Hausdorff outer measure

abstract: We use nonstandard techniques, in the sense of Abraham Robinson, to give the exact Hausdorff outer measure.


title: Demonstration of Minisuperspace Quantum Cosmology Using Quantum Computational Algorithms on IBM Quantum Computer

abstract: Quantum computers promise to efficiently solve important problems that are intractable on a conventional computer. Quantum computational algorithms have the potential to be an exciting new way of studying quantum cosmology. In quantum cosmology, we learn about the dynamics of the universe without constructing a complete theory of quantum gravity. Since the universal wavefunction exists in an infinite-dimensional superspace over all possible 3D metrics and modes of matter configurations, we take minisuperspaces for our work by constraining the degrees of freedom to particular 3D metrics and uniform scalar field configurations. Here, we consider a wide variety of cosmological models. We begin by analyzing an anisotropic universe with cosmological constant and classical radiation. We then study the results for higher derivatives, Kaluza-Klein theories and string dilaton in quantum cosmology. We use IBM's Quantum Information Science Kit (QISKit) python library and the Variational Quantum Eigensolver (VQE) algorithm for studying these systems. The VQE algorithm is a hybrid algorithm that uses the variational approach and interleaves quantum and classical computations in order to find the minimum eigenvalue of the Hamiltonian for a given system.


title: HIJING can describe the anisotropy-scaled charge-dependent correlations at the Relativistic Heavy Ion Collider

abstract: The experimentally measured charge-depdendent correlations in heavy ion collisions have been suggested as a signature of the chiral magenetic effect (CME). Early model studies could not reproduce the measurement. For example, the Hijing model yielded far smaller magnitude for the charge-dependent correlation than observed in data. This led to the conclusion that the CME had to be invoked to explain the observed correlations in heavy ion collisions. In this paper we show that this conclusion of the CME interpretation is premature. We show that the reason that Hijing predicts a far smaller correlation than data is because the elliptic anisotrpy ($v_{2}$) parameter in Hijing is too small. When properly scaled, the Hijing model can reproduce in entirety the measured correlations. We also employ the AMPT model, which has a large enough $v_2$, to demonstrate that the measured data can be easily accommodated by models without invoking the CME.


title: Clustering in $^{18}$O -- absolute determination of branching ratios via high-resolution particle spectroscopy

abstract: The determination of absolute branching ratios for high-energy states in light nuclei is an important and useful tool for probing the underlying nuclear structure of individual resonances: for example, in establishing the tendency of an excited state towards $\alpha$-cluster structure. Difficulty arises in measuring these branching ratios due to similarities in available decay channels, such as ($\mathbf{^{18}}$O,$\mathbf{n}$) and ($\mathbf{^{18}}$O,$\mathbf{2n}$), as well as differences in geometric efficiencies due to population of bound excited levels in daughter nuclei. Methods are presented using Monte Carlo techniques to overcome these issues.


title: Every Planar Set has a Conformally Removable Subset with the Same Hausdorff Dimension

abstract: In this paper we show that given any compact set $E \subset \hat{\mathbb{C}}$, we can always find a conformally removable subset with the same Hausdorff dimension as $E$.


abstract: In this paper, we compute sub-Riemannian limits of Gaussian curvature for a Euclidean $C^2$-smooth surface in the affine group and the group of rigid motions of the Minkowski plane away from characteristic points and signed geodesic curvature for Euclidean $C^2$-smooth curves on surfaces. We get Gauss-Bonnet theorems in the affine group and the group of rigid motions of the Minkowski plane.


title: MANELA: A Multi-Agent Algorithm for Learning Network Embeddings

abstract: Playing an essential role in data mining, machine learning has a long history of being applied to networks on multifarious tasks and has played an essential role in data mining. However, the discrete and sparse natures of networks often render it difficult to apply machine learning directly to networks. To circumvent this difficulty, one major school of thought to approach networks using machine learning is via network embeddings. On the one hand, this network embeddings have achieved huge success on aggregated network data in recent years. On the other hand, learning network embeddings on distributively stored networks still remained understudied: To the best of our knowledge, all existing algorithms for learning network embeddings have hitherto been exclusively centralized and thus cannot be applied to these networks. To accommodate distributively stored networks, in this paper, we proposed a multi-agent model. Under this model, we developed the multi-agent network embedding learning algorithm (MANELA) for learning network embeddings. We demonstrate MANELA's advantages over other existing centralized network embedding learning algorithms both theoretically and experimentally. Finally, we further our understanding in MANELA via visualization and exploration of its relationship to DeepWalk.


title: Borrowing From the Future: An Attempt to Address Double Sampling

abstract: For model-free reinforcement learning, one of the main difficulty of stochastic Bellman residual minimization is the double sampling problem, i.e., while only one single sample for the next state is available in the model-free setting, two independent samples for the next state are required in order to perform unbiased stochastic gradient descent. We propose new algorithms for addressing this problem based on the idea of borrowing extra randomness from the future. When the transition kernel varies slowly with respect to the state, it is shown that the training trajectory of new algorithms is close to the one of unbiased stochastic gradient descent. Numerical results for policy evaluation in both tabular and neural network settings are provided to confirm the theoretical findings.


title: Deconfinement and chiral phase transitions in quark matter with a strong electric field

abstract: The deconfinement and chiral phase transitions are studied in the context of the electrized quark matter at finite temperature in the two-flavor Polyakov-Nambu--Jona-Lasinio model. Using the mean field approximation and an electric field independet regularization we show that the effect of temperature and/or electric fields is to partially restore the chiral symmetry. The deconfinement phase transition is slightly affected by the magnitude of the electric field. To this end we show how the effective quark masses and the expectation value of the Polyakov Loop are affected by the electric fields at finite temperatures. As a very interesting result, the pseudocritical temperatures for chiral symmetry restoration and deconfinement decrease as we increase the magnitude of the electric fields, however, both start to increase after some critical value of the electric field.


title: Efficient adjustment sets for population average treatment effect estimation in non-parametric causal graphical models

abstract: The method of covariate adjustment is often used for estimation of population average treatment effects in observational studies. Graphical rules for determining all valid covariate adjustment sets from an assumed causal graphical model are well known. Restricting attention to causal linear models, a recent article derived two novel graphical criteria: one to compare the asymptotic variance of linear regression treatment effect estimators that control for certain distinct adjustment sets and another to identify the optimal adjustment set that yields the least squares treatment effect estimator with the smallest asymptotic variance among consistent adjusted least squares estimators. In this paper we show that the same graphical criteria can be used in non-parametric causal graphical models when treatment effects are estimated by contrasts involving non-parametrically adjusted estimators of the interventional means. We also provide a graphical criterion for determining the optimal adjustment set among the minimal adjustment sets, which is valid for both linear and non-parametric estimators. We provide a new graphical criterion for comparing time dependent adjustment sets, that is, sets comprised by covariates that adjust for future treatments and that are themselves affected by earlier treatments. We show by example that uniformly optimal time dependent adjustment sets do not always exist. In addition, for point interventions, we provide a sound and complete graphical criterion for determining when a non-parametric optimally adjusted estimator of an interventional mean, or of a contrast of interventional means, is as efficient as an efficient estimator of the same parameter that exploits the information in the conditional independencies encoded in the non-parametric causal graphical model.


title: Estimating the population mean for a vertical profile of energy dissipation rate

abstract: Energy dissipation rates are an important characteristic of turbulence; however, their magnitude in observational profiles can be incorrectly determined owing to their irregular appearance during vertical evolution. By analysing the data obtained from oceanic turbulence measurements, we demonstrate that the vertical sequences of energy dissipation rates exhibit a scaling property. Utilising this property, we propose a method to estimate the population mean for a profile. For scaling in the observed profiles, we demonstrate that our data exhibit a statistical property consistent with that exhibited by the universal multifractal model. Meanwhile, the population mean and its uncertainty can be estimated by inverting the probability distribution obtained by Monte Carlo simulations of a cascade model; to this end, observational constraints from several moments are imposed over each vertical sequence. This approach enables us to determine, to some extent, whether a profile shows an occasionally large mean or whether the population mean itself is large. Thus, it will contribute to the refinement of the regional estimation of the ocean energy budget, where only a small amount of turbulence observation data is available.


title: Exploiting Motion Information from Unlabeled Videos for Static Image Action Recognition

abstract: Static image action recognition, which aims to recognize action based on a single image, usually relies on expensive human labeling effort such as adequate labeled action images and large-scale labeled image dataset. In contrast, abundant unlabeled videos can be economically obtained. Therefore, several works have explored using unlabeled videos to facilitate image action recognition, which can be categorized into the following two groups: (a) enhance visual representations of action images with a designed proxy task on unlabeled videos, which falls into the scope of self-supervised learning; (b) generate auxiliary representations for action images with the generator learned from unlabeled videos. In this paper, we integrate the above two strategies in a unified framework, which consists of Visual Representation Enhancement (VRE) module and Motion Representation Augmentation (MRA) module. Specifically, the VRE module includes a proxy task which imposes pseudo motion label constraint and temporal coherence constraint on unlabeled videos, while the MRA module could predict the motion information of a static action image by exploiting unlabeled videos. We demonstrate the superiority of our framework based on four benchmark human action datasets with limited labeled data.


title: A P\'{o}sa-type condition of potentially $_3C_\ell$-graphic sequences

abstract: A non-increasing sequence $\pi=(d_1,\ldots,d_n)$ of nonnegative integers is said to be graphic if it is realizable by a simple graph $G$ on $n$ vertices. A graphic sequence $\pi=(d_1,\ldots,d_n)$ is said to be potentially $_3C_\ell$-graphic if there is a realization of $\pi$ containing cycles of every length $r$, $3\le r\le \ell$. It is well-known that if the non-increasing degree sequence $(d_1,\ldots,d_\ell)$ of a graph $G$ on $\ell$ vertices satisfies the P\'{o}sa condition that $d_{\ell+1-i}\ge i+1$ for every $i$ with $1\le i<\frac{\ell}{2}$, then $G$ is either pancyclic or bipartite. In this paper, we obtain a P\'{o}sa-type condition of potentially $_3C_\ell$-graphic sequences, that is, we prove that if $\ell\ge 5$ is an integer, $n\ge \ell$ and $\pi=(d_1,\ldots,d_n)$ is a graphic sequence with $d_{\ell+1-i}\ge i+1$ for every $i$ with $1\le i<\frac{\ell}{2}$, then $\pi$ is potentially $_3C_\ell$-graphic. This result improves a Dirac-type condition of potentially $_3C_\ell$-graphic sequences due to Yin et al. [Appl. Math. Comput., 353 (2019) 88--94], and asymptotically answers a problem due to Li et al. [Adv. Math., 33 (2004) 273--283]. As an application, this result also completely implies the value $\sigma(C_\ell,n)$ for $\ell\ge 5$ and $n\ge \ell$, improving the result of Lai [J. Combin. Math. Combin. Comput., 49 (2004) 57--64].


title: Coarse-graining in micromagnetic simulations of dynamic hysteresis loops

abstract: Micromagnetic simulations based on the stochastic Landau-Lifshitz-Gilbert equation are used to calculate dynamic magnetic hysteresis loops relevant to magnetic hyperthermia. With the goal to effectively simulate room-temperature loops for large iron-oxide-based systems at relatively slow sweep rates on the order of 1 Oe/ns or less, a previously derived renormalization group approach for coarse-graining (Grinstein and Koch, Phys. Rev. Lett. 20, 207201, 2003) is modified and applied to calculating loops for a magnetite nanorod. The nanorod modelled is the building block for larger nanoparticles that were employed in preclinical studies (Dennis et al., Nanotechnology 20, 395103, 2009). The scaling algorithm is shown to produce nearly identical loops over several decades in the model grain size. Sweep-rate scaling involving the Gilbert damping parameter is also demonstrated to allow orders of magnitude speed-up of the loop calculations.


title: Learning to Relate from Captions and Bounding Boxes

abstract: In this work, we propose a novel approach that predicts the relationships between various entities in an image in a weakly supervised manner by relying on image captions and object bounding box annotations as the sole source of supervision. Our proposed approach uses a top-down attention mechanism to align entities in captions to objects in the image, and then leverage the syntactic structure of the captions to align the relations. We use these alignments to train a relation classification network, thereby obtaining both grounded captions and dense relationships. We demonstrate the effectiveness of our model on the Visual Genome dataset by achieving a recall@50 of 15% and recall@100 of 25% on the relationships present in the image. We also show that the model successfully predicts relations that are not present in the corresponding captions.


title: Embodiment in Socially Interactive Robots

abstract: Physical embodiment is a required component for robots that are structurally coupled with their real-world environments. However, most socially interactive robots do not need to physically interact with their environments in order to perform their tasks. When and why should embodied robots be used instead of simpler and cheaper virtual agents? This paper reviews the existing work that explores the role of physical embodiment in socially interactive robots. This class consists of robots that are not only capable of engaging in social interaction with humans, but are using primarily their social capabilities to perform their desired functions. Socially interactive robots provide entertainment, information, and/or assistance; this last category is typically encompassed by socially assistive robotics. In all cases, such robots can achieve their primary functions without performing functional physical work. To comprehensively evaluate the existing body of work on embodiment, we first review work from established related fields including psychology, philosophy, and sociology. We then systematically review 65 studies evaluating aspects of embodiment published from 2003 to 2017 in major peer-reviewed robotics publication venues. We examine relevant aspects of the selected studies, focusing on the embodiments compared, tasks evaluated, social roles of robots, and measurements. We introduce three taxonomies for the types of robot embodiment, robot social roles, and human-robot tasks. These taxonomies are used to deconstruct the design and interaction spaces of socially interactive robots and facilitate analysis and discussion of the reviewed studies. We use this newly-defined methodology to critically discuss existing works, revealing topics within embodiment research for social interaction, assistive robotics, and service robotics.


title: On Morawetz estimates with time-dependent weights for the Klein-Gordon equation

abstract: We obtain some new Morawetz estimates for the Klein-Gordon flow of the form \begin{equation*} \big\||\nabla|^{\sigma} e^{it \sqrt{1-\Delta}}f \big\|_{L^2_{x,t}(|(x,t)|^{-\alpha})} \lesssim \|f\|_{H^s} \end{equation*} where $\sigma,s\geq0$ and $\alpha>0$. The conventional approaches to Morawetz estimates with $|x|^{-\alpha}$ are no longer available in the case of time-dependent weights $|(x,t)|^{-\alpha}$. Here we instead apply the Littlewood-Paley theory with Muckenhoupt $A_2$ weights to frequency localized estimates thereof that are obtained by making use of the bilinear interpolation between their bilinear form estimates which need to carefully analyze some relevant oscillatory integrals according to the different scaling of $\sqrt{1-\Delta}$ for low and high frequencies.


title: ACE -- An Anomaly Contribution Explainer for Cyber-Security Applications

abstract: In this paper, we introduce Anomaly Contribution Explainer or ACE, a tool to explain security anomaly detection models in terms of the model features through a regression framework, and its variant, ACE-KL, which highlights the important anomaly contributors. ACE and ACE-KL provide insights in diagnosing which attributes significantly contribute to an anomaly by building a specialized linear model to locally approximate the anomaly score that a black-box model generates. We conducted experiments with these anomaly detection models to detect security anomalies on both synthetic data and real data. In particular, we evaluate performance on three public data sets: CERT insider threat, netflow logs, and Android malware. The experimental results are encouraging: our methods consistently identify the correct contributing feature in the synthetic data where ground truth is available; similarly, for real data sets, our methods point a security analyst in the direction of the underlying causes of an anomaly, including in one case leading to the discovery of previously overlooked network scanning activity. We have made our source code publicly available.


title: Topic-aware chatbot using Recurrent Neural Networks and Nonnegative Matrix Factorization

abstract: We propose a novel model for a topic-aware chatbot by combining the traditional Recurrent Neural Network (RNN) encoder-decoder model with a topic attention layer based on Nonnegative Matrix Factorization (NMF). After learning topic vectors from an auxiliary text corpus via NMF, the decoder is trained so that it is more likely to sample response words from the most correlated topic vectors. One of the main advantages in our architecture is that the user can easily switch the NMF-learned topic vectors so that the chatbot obtains desired topic-awareness. We demonstrate our model by training on a single conversational data set which is then augmented with topic matrices learned from different auxiliary data sets. We show that our topic-aware chatbot not only outperforms the non-topic counterpart, but also that each topic-aware model qualitatively and contextually gives the most relevant answer depending on the topic of question.


title: Equivariant cohomology for differentiable stacks

abstract: We construct and analyse models of equivariant cohomology for differentiable stacks with Lie group actions extending classical results for smooth manifolds due to Borel, Cartan and Getzler. We also derive various spectral sequences for the equivariant cohomology of a differentiable stack generalising among others Bott's spectral sequence which converges to the cohomology of the classifying space of a Lie group.