Skeleton DML Deep Metric Learning for Skeleton Based One Shot Action Recognition. One shot action recognition allows the recognition of human performed actions with only a single training example. This can influence human robot interaction positively by enabling the robot to react to previously unseen behaviour. We formulate the one shot action recognition problem as a deep metric learning problem and propose a novel image based skeleton representation that performs well in a metric learning setting. Therefore, we train a model that projects the image representations into an embedding space. In embedding space the similar actions have a low euclidean distance while dissimilar actions have a higher distance. The one shot action recognition problem becomes a nearest neighbor search in a set of activity reference samples. We evaluate the performance of our proposed representation against a variety of other skeleton based image representations. In addition, we present an ablation study that shows the influence of different embedding vector sizes, losses and augmentation. Our approach lifts the state of the art by 3.3 for the one shot action recognition protocol on the NTU RGB D 120 dataset under a comparable training setup. With additional augmentation our result improved over 7.7 .


Unification and combination of a class of traversal strategies made with pattern matching and fixed points. Motivated by an ongoing project on computer aided derivation of asymptotic models governed by partial differential equations, we introduce a class of term transformations that consists of traversal strategies and insertion of contexts. We define unification and combination operations on this class which amount to merging transformations in order to obtain more complex ones. We show that the unification and combination operations enjoy nice algebraic properties like associativity, congruence and the existence of neutral elements. The main part of this paper is devoted to proving that the unification and combination operations are correct.


An optimization problem in thermal insulation with Robin boundary conditions. We study thermal insulating of a bounded body Omega subset mathbb R n . Under a prescribed heat source f geq 0 , we consider a model of heat transfer between Omega and the environment determined by convection this corresponds, before insulation, to Robin boundary conditions. The body is then surrounded by a layer of insulating material of thickness of size varepsilon 0 , and whose conductivity is also proportional to varepsilon . This corresponds to the case of a small amount of insulating material, with excellent insulating properties. We then compute the Gamma limit of the energy functional F varepsilon and prove that this is a functional F whose minimizers still satisfy an elliptic PDEs system with a non uniform Robin boundary condition depending on the distribution of insulating layer around Omega . In a second step we study the maximization of heat content (which measures the goodness of the insulation) among all the possible distributions of insulating material with fixed mass, and prove an optimal upper bound in terms of geometric properties. Eventually we prove a conjecture which states that the ball surrounded by a uniform distribution of insulating material maximizes the heat content.


Stationary fully nonlinear mean field games. In this paper we examine fully nonlinear mean field games associated with a minimization problem. The variational setting is driven by a functional depending on its argument through its Hessian matrix. We work under fairly natural conditions and establish improved (sharp) regularity for the solutions in Sobolev spaces. Then, we prove the existence of minimizers for the variational problem and the existence of solutions to the mean field games system. We also investigate a unidimensional example and unveil new information on the explicit solutions. Our findings can be generalized to a larger class of operators, yielding information on a broader range of examples.


Coherent control of reactive scattering at low temperatures Signatures of quantum interference in the differential cross sections for F H2 and F HD. Fundamental entanglement related challenges have prevented quantum interference based control (i.e. coherent control) of collisional cross sections from being implemented in the laboratory. Here, differential cross sections for reactive scattering at low temperatures are shown to provide a unique opportunity to display such interference based control by forming coherent superpositions of degenerate rotational states of reactant molecules jmi with different m. In particular, we identify and quantify a unique signature of coherent control in reactive scattering with applications to F H2 ! H HF and HF D F HD ! HD F at 11 K. Control is shown to be extensive.


Rainbow matchings for 3 uniform hypergraphs. K " u hn, Osthus, and Treglown and, independently, Khan proved that if H is a 3 uniform hypergraph with n vertices such that n in 3 mathbb Z and large, and delta 1(H) n 1 choose 2 2n 3 choose 2 , then H contains a perfect matching. In this paper, we show that for n in 3 mathbb Z sufficiently large, if F 1, ldots, F n 3 are 3 uniform hypergrapghs with a common vertex set and delta 1(F i) n 1 choose 2 2n 3 choose 2 for i in n 3 , then F 1, dots, F n 3 admits a rainbow matching, i.e., a matching consisting of one edge from each F i . This is done by converting the rainbow matching problem to a perfect matching problem in a special class of uniform hypergraphs.


Teaching Cultural Astronomy to Undergraduates with an Interdisciplinary Frame. Cultural Astronomy is interdisciplinary connecting the arts, humanities, social physical sciences. Data collection methods and theories are used from many disciplines and meld with methods and theories within cultural astronomy. The burden on the student is that to do cultural astronomy research it is necessary to be widely read within and across disciplines. I developed a series of courses that divided cultural astronomy content into broad regions such as Africa, North America, and the Pacific. The courses were structured to accommodate students from all parts of the university, but had to have enough mathematics and science to serve as a general science requirement. The majority of the grade for the course lay with the final project, which included a presentation and a written document. Thus, the course was designed to give the students a foundation for doing this final project that had to be original research. The students rarely opted to collect their own data, instead they reanalyzed existing materials. The students learned critical thinking, formulating hypotheses, and how to test their hypotheses, as well as how to understand cultural astronomy data.


Biologically Plausible Learning of Text Representation with Spiking Neural Networks. This study proposes a novel biologically plausible mechanism for generating low dimensional spike based text representation. First, we demonstrate how to transform documents into series of spikes spike trains which are subsequently used as input in the training process of a spiking neural network (SNN). The network is composed of biologically plausible elements, and trained according to the unsupervised Hebbian learning rule, Spike Timing Dependent Plasticity (STDP). After training, the SNN can be used to generate low dimensional spike based text representation suitable for text document classification. Empirical results demonstrate that the generated text representation may be effectively used in text classification leading to an accuracy of 80.19 on the bydate version of the 20 newsgroups data set, which is a leading result amongst approaches that rely on low dimensional text representations.


Geometric realizations of regular abstract polyhedra with automorphism group H 3 . A textit geometric realization of an abstract polyhedron mathcal P is a mapping rho mathcal P to mathbb E 3 that sends an i face to an open set of dimension i . This work adapts a method based on Wythoff construction to generate a full rank realization of a regular abstract polyhedron from its automorphism group Gamma . The method entails finding a real orthogonal representation of Gamma of degree 3 and applying its image to suitably chosen open sets in space. To demonstrate the use of the method, we apply it to the abstract polyhedra whose automorphism groups are isomorphic to the non crystallographic Coxeter group H 3 .


How to Distribute Computation in Networks. In network function computation is as a means to reduce the required communication flow in terms of number of bits transmitted per source symbol. However, the rate region for the function computation problem in general topologies is an open problem, and has only been considered under certain restrictive assumptions (e.g. tree networks, linear functions, etc.). In this paper, we propose a new perspective for distributing computation, and formulate a flow based delay cost minimization problem that jointly captures the costs of communications and computation. We introduce the notion of entropic surjectivity as a measure to determine how sparse the function is and to understand the limits of computation. Exploiting Little's law for stationary systems, we provide a connection between this new notion and the computation processing factor that reflects the proportion of flow that requires communications. This connection gives us an understanding of how much a node (in isolation) should compute to communicate the desired function within the network without putting any assumptions on the topology. Our analysis characterizes the functions only via their entropic surjectivity, and provides insight into how to distribute computation. We numerically test our technique for search, MapReduce, and classification tasks, and infer for each task how sensitive the processing factor to the entropic surjectivity is.


Encounter rates between bacteria and small sinking particles. Bacteria in aquatic environments often interact with particulate matter. A key example is bacterial degradation of marine snow responsible for carbon export from the upper ocean in the biological pump. The ecological interaction between bacteria and sinking particles is regulated by their encounter rate, which is therefore important to predict accurately in models of bacteria particle interactions. Models available to date cover the diffusive encounter regime, valid for sinking particles larger than the typical run length of a bacterium. The majority of sinking particles, however, are small, and the encounter process is then ballistic rather than diffusive. In the ballistic regime, the shear generated by the particle's motion can be important in reorienting bacteria and thus determining the encounter rate, yet the effect of shear is not captured in current encounter rate models. Here, we combine analytical and numerical calculations to quantify the encounter rate between sinking particles and non motile or motile microorganisms in the ballistic regime, explicitly accounting for the hydrodynamic shear created by the particle and its coupling with microorganism shape. We complement results with selected experiments on non motile diatoms. We find that the shape shear coupling has a considerable effect on the encounter rate and encounter location through the mechanisms of hydrodynamic focusing and screening, whereby elongated microorganisms preferentially orient normally to the particle surface downstream of the particle (focusing) and tangentially to the particle surface upstream of the particle (screening). We study these mechanisms as a function of the key dimensionless parameters the ratio of particle sinking speed to microorganism swimming speed, the ratio of particle radius to microorganism length, and the microorganism's aspect ratio.


Impact of Network and Host Characteristics on the Keystroke Pattern in Remote Desktop Sessions. Authentication based on keystroke dynamics is a convenient biometric approach, easy in use, transparent, and cheap as it does not require a dedicated sensor. Keystroke authentication, as part of multi factor authentication, can be used in remote display access to guarantee the security of use of remote connectivity systems during the access control phase or throughout the session. This paper investigates how network conditions and additional host interaction may impact the behavioural pattern of keystrokes when used in a remote desktop application scenario. We focus on the timing of adjacent keys and investigate this impact by calculating the variations of the Euclidean distance between a reference profile and resulting profiles following such impairments. The experimental results indicate that variations of congestion latency, whether produced by adjacent traffic sources or by additional remote desktop interactions, have a substantive impact on the Euclidian distance, which in turn may affect the effectiveness of the biometric authentication algorithm. Results also indicate that data flows within remote desktop protocol are not prioritized and therefore additional traffic will have a significant impact on the keystroke timings, which renders continuous authentication less effective for remote access and more appropriate for one time login.


Liouville type theorems on manifolds with nonnegative curvature and strictly convex boundary. We prove some Liouville type theorems on smooth compact Riemannian manifolds with nonnegative sectional curvature and strictly convex boundary. This gives a nonlinear generalization in low dimension of the recent sharp lower bound of the first Steklov eigenvalue by Xia Xiong and verifies partially a conjecture by the third author. As a consequence, we derive several sharp Sobolev trace inequalities on these manifolds.


Integrating Deep Neural Networks with Full waveform Inversion Reparametrization, Regularization, and Uncertainty Quantification. Full waveform inversion (FWI) is an accurate imaging approach for modeling velocity structure by minimizing the misfit between recorded and predicted seismic waveforms. However, the strong non linearity of FWI resulting from fitting oscillatory waveforms can trap the optimization in local minima. We propose a neural network based full waveform inversion method (NNFWI) that integrates deep neural networks with FWI by representing the velocity model with a generative neural network. Neural networks can naturally introduce spatial correlations as regularization to the generated velocity model, which suppresses noise in the gradients and mitigates local minima. The velocity model generated by neural networks is input to the same partial differential equation (PDE) solvers used in conventional FWI. The gradients of both the neural networks and PDEs are calculated using automatic differentiation, which back propagates gradients through the acoustic PDEs and neural network layers to update the weights of the generative neural network. Experiments on 1D velocity models, the Marmousi model, and the 2004 BP model demonstrate that NNFWI can mitigate local minima, especially for imaging high contrast features like salt bodies, and significantly improves the inversion in the presence of noise. Adding dropout layers to the neural network model also allows analyzing the uncertainty of the inversion results through Monte Carlo dropout. NNFWI opens a new pathway to combine deep learning and FWI for exploiting both the characteristics of deep neural networks and the high accuracy of PDE solvers. Because NNFWI does not require extra training data and optimization loops, it provides an attractive and straightforward alternative to conventional FWI.


On the Demystification of Knowledge Distillation A Residual Network Perspective. Knowledge distillation (KD) is generally considered as a technique for performing model compression and learned label smoothing. However, in this paper, we study and investigate the KD approach from a new perspective we study its efficacy in training a deeper network without any residual connections. We find that in most of the cases, non residual student networks perform equally or better than their residual versions trained on raw data without KD (baseline network). Surprisingly, in some cases, they surpass the accuracy of baseline networks even with the inferior teachers. After a certain depth of non residual student network, the accuracy drop, coming from the removal of residual connections, is substantial, and training with KD boosts the accuracy of the student up to a great extent however, it does not fully recover the accuracy drop. Furthermore, we observe that the conventional teacher student view of KD is incomplete and does not adequately explain our findings. We propose a novel interpretation of KD with the Trainee Mentor hypothesis, which provides a holistic view of KD. We also present two viewpoints, loss landscape, and feature reuse, to explain the interplay between residual connections and KD. We substantiate our claims through extensive experiments on residual networks.


Energy Expenditure Estimation Through Daily Activity Recognition Using a Smart phone. This paper presents a 3 step system that estimates the real time energy expenditure of an individual in a non intrusive way. First, using the user's smart phone's sensors, we build a Decision Tree model to recognize his physical activity ( textit running , textit standing , ...). Then, we use the detected physical activity, the time and the user's speed to infer his daily activity ( textit watching TV , textit going to the bathroom , ...) through the use of a reinforcement learning environment, the Partially Observable Markov Decision Process framework. Once the daily activities are recognized, we translate this information into energy expenditure using the compendium of physical activities. By successfully detecting 8 physical activities at 90 , we reached an overall accuracy of 80 in recognizing 17 different daily activities. This result leads us to estimate the energy expenditure of the user with a mean error of 26 of the expected estimation.


Twisted Dolbeault cohomology of nilpotent Lie algebras. It is well known that cohomology of any non trivial 1 dimensional local system on a nilmanifold vanishes (this result is due to L. Alaniya). A complex nilmanifold is a quotient of a nilpotent Lie group equipped with a left invariant complex structure by an action of a discrete, co compact subgroup. We prove a Dolbeault version of Alaniya's theorem, showing that the Dolbeault cohomology of a nilpotent Lie algebra with coefficients in any non trivial 1 dimensional local system vanishes. Note that the Dolbeault cohomology of the corresponding local system on the manifold is not necessarily zero. This implies that the twisted version of Console Fino theorem is false (Console Fino proved that the Dolbeault cohomology of a complex nilmanifold is equal to the Dolbeault cohomology of its Lie algebra). As an application, we give a new proof of a theorem due to H. Sawai, who obtained an explicit description of LCK nilmanifolds. An LCK structure on a manifold M is a K "ahler structure on its cover tilde M such that the deck transform map acts on tilde M by homotheties. We show that any complex nilmanifold admitting an LCK structure is Vaisman, and is obtained as a compact quotient of the product of a Heisenberg group and the real line.


Watchlist Risk Assessment using Multiparametric Cost and Relative Entropy. This paper addresses the facial biometric enabled watchlist technology in which risk detectors are mandatory mechanisms for early detection of threats, as well as for avoiding offense to innocent travelers. We propose a multiparametric cost assessment and relative entropy measures as risk detectors. We experimentally demonstrate the effects of mis identification and impersonation under various watchlist screening scenarios and constraints. The key contributions of this paper are the novel techniques for design and analysis of the biometric enabled watchlist and the supporting infrastructure, as well as measuring the impersonation impact on e border performance.


Interpreting and Improving Adversarial Robustness of Deep Neural Networks with Neuron Sensitivity. Deep neural networks (DNNs) are vulnerable to adversarial examples where inputs with imperceptible perturbations mislead DNNs to incorrect results. Despite the potential risk they bring, adversarial examples are also valuable for providing insights into the weakness and blind spots of DNNs. Thus, the interpretability of a DNN in the adversarial setting aims to explain the rationale behind its decision making process and makes deeper understanding which results in better practical applications. To address this issue, we try to explain adversarial robustness for deep models from a new perspective of neuron sensitivity which is measured by neuron behavior variation intensity against benign and adversarial examples. In this paper, we first draw the close connection between adversarial robustness and neuron sensitivities, as sensitive neurons make the most non trivial contributions to model predictions in the adversarial setting. Based on that, we further propose to improve adversarial robustness by constraining the similarities of sensitive neurons between benign and adversarial examples which stabilizes the behaviors of sensitive neurons towards adversarial noises. Moreover, we demonstrate that state of the art adversarial training methods improve model robustness by reducing neuron sensitivities which in turn confirms the strong connections between adversarial robustness and neuron sensitivity as well as the effectiveness of using sensitive neurons to build robust models. Extensive experiments on various datasets demonstrate that our algorithm effectively achieves excellent results.


Positive solutions for a Minkowski curvature equation with indefinite weight and super exponential nonlinearity. We investigate the existence of positive solutions for a class of Minkowski curvature equations with indefinite weight and nonlinear term having superlinear growth at zero and super exponential growth at infinity. As an example, for the equation begin equation Biggl ( dfrac u' sqrt 1 (u') 2 Biggr ) ' a(t) bigl ( e u p 1 bigr ) 0, end equation where p 1 and a(t) is a sign changing function satisfying the mean value condition int 0 T a(t) , mathrm d t 0 , we prove the existence of a positive solution for both periodic and Neumann boundary conditions. The proof relies on a topological degree technique.


Generating Rationales in Visual Question Answering. Despite recent advances in Visual QuestionAnswering (VQA), it remains a challenge todetermine how much success can be attributedto sound reasoning and comprehension ability.We seek to investigate this question by propos ing a new task ofrationale generation. Es sentially, we task a VQA model with generat ing rationales for the answers it predicts. Weuse data from the Visual Commonsense Rea soning (VCR) task, as it contains ground truthrationales along with visual questions and an swers. We first investigate commonsense un derstanding in one of the leading VCR mod els, ViLBERT, by generating rationales frompretrained weights using a state of the art lan guage model, GPT 2. Next, we seek to jointlytrain ViLBERT with GPT 2 in an end to endfashion with the dual task of predicting the an swer in VQA and generating rationales. Weshow that this kind of training injects com monsense understanding in the VQA modelthrough quantitative and qualitative evaluationmetrics


An efficient method for approximating resonance curves of weakly damped nonlinear mechanical systems. A method is presented for tracing the locus of a specific peak in the frequency response under variation of a parameter. It is applicable to periodic, steady state vibrations of harmonically forced nonlinear mechanical systems. It operates in the frequency domain and its central idea is to assume a constant phase lag between forcing and response. The method is validated for a two degree of freedom oscillator with cubic spring and a bladed disk with shroud contact. The method provides superior computational efficiency, but is limited to weakly damped systems. Finally, the capability to reveal isolated solution branches is highlighted.


Explicit solutions for a class of nonlinear backward stochastic differential equations and their nodal sets. In this paper, we investigate a class of nonlinear backward stochastic differential equations (BSDEs) arising from financial economics, and give specific information about the nodal sets of the related solutions. As applications, we are able to obtain the explicit solutions to an interesting class of nonlinear BSDEs including the k ignorance BSDE arising from the modeling of ambiguity of asset pricing.


Semantic Change Pattern Analysis. Change detection is an important problem in vision field, especially for aerial images. However, most works focus on traditional change detection, i.e., where changes happen, without considering the change type information, i.e., what changes happen. Although a few works have tried to apply semantic information to traditional change detection, they either only give the label of emerging objects without taking the change type into consideration, or set some kinds of change subjectively without specifying semantic information. To make use of semantic information and analyze change types comprehensively, we propose a new task called semantic change pattern analysis for aerial images. Given a pair of co registered aerial images, the task requires a result including both where and what changes happen. We then describe the metric adopted for the task, which is clean and interpretable. We further provide the first well annotated aerial image dataset for this task. Extensive baseline experiments are conducted as reference for following works. The aim of this work is to explore high level information based on change detection and facilitate the development of this field with the publicly available dataset.


Looking Beyond Sentence Level Natural Language Inference for Downstream Tasks. In recent years, the Natural Language Inference (NLI) task has garnered significant attention, with new datasets and models achieving near human level performance on it. However, the full promise of NLI particularly that it learns knowledge that should be generalizable to other downstream NLP tasks has not been realized. In this paper, we study this unfulfilled promise from the lens of two downstream tasks question answering (QA), and text summarization. We conjecture that a key difference between the NLI datasets and these downstream tasks concerns the length of the premise and that creating new long premise NLI datasets out of existing QA datasets is a promising avenue for training a truly generalizable NLI model. We validate our conjecture by showing competitive results on the task of QA and obtaining the best reported results on the task of Checking Factual Correctness of Summaries.


Fourier holographic endoscopy for label free imaging through a narrow and curved passage. Ultrathin lensless fibre endoscopes offer minimally invasive investigation, but they mostly operate as a rigid type due to the need for prior calibration of a fibre probe. Furthermore, most implementations work in fluorescence mode rather than label free imaging mode, making them unsuitable for medicine and industry. Herein, we report a fully flexible ultrathin fibre endoscope taking 3D holographic images of unstained tissues with 0.87 mu m spatial resolution. Using a bare fibre bundle as thin as 200 mu m diameter, we design a lensless Fourier holographic imaging configuration to selectively detect weak reflections from biological tissues, a critical step for stain free reflectance imaging. A unique algorithm is developed for calibration free holographic image reconstruction, allowing us to image through a narrow and curved passage regardless of fibre bending. We demonstrate endoscopic reflectance imaging of unstained rat intestine tissues that are completely invisible to conventional endoscopes. The proposed endoscope will expedite more accurate and earlier diagnosis than before with minimal complications.


Non Local Pearson diffusions. In this paper we focus on strong solutions of some heat like problems with a non local derivative in time induced by a Bernstein function and an elliptic operator given by the generator or the Fokker Planck operator of a Pearson diffusion. Such kind of non local equations naturally arise in the treatment of particle motion in heterogeneous media. In particular, we use spectral decomposition results for the usual Pearson diffusion to exploit explicit solutions of the aforementioned equations. Moreover, we provide stochastic representation of such solutions in terms of time changed Pearson diffusions. Finally, we exploit some further properties of these processes, such as limit distributions and long short range dependence.


D2RL Deep Dense Architectures in Reinforcement Learning. While improvements in deep learning architectures have played a crucial role in improving the state of supervised and unsupervised learning in computer vision and natural language processing, neural network architecture choices for reinforcement learning remain relatively under explored. We take inspiration from successful architectural choices in computer vision and generative modelling, and investigate the use of deeper networks and dense connections for reinforcement learning on a variety of simulated robotic learning benchmark environments. Our findings reveal that current methods benefit significantly from dense connections and deeper networks, across a suite of manipulation and locomotion tasks, for both proprioceptive and image based observations. We hope that our results can serve as a strong baseline and further motivate future research into neural network architectures for reinforcement learning. The project website with code is at this link https sites.google.com view d2rl home.


Local extension estimates for the hyperbolic hyperboloid in three dimensions. We establish Fourier extension estimates for compact subsets of the hyperbolic hyperboloid in three dimensions via polynomial partitioning.


Psybrackets, Pseudoknots and Singular Knots. We introduce algebraic structures known as psybrackets and use them to define invariants of pseudoknots and singular knots and links. Psybrackets are Niebrzydowski tribrackets with additional structure inspired by the Reidemeister moves for pseudoknots and singular knots. Examples and computations are provided.


Stochastic Matching with Few Queries (1 varepsilon) Approximation. Suppose that we are given an arbitrary graph G (V, E) and know that each edge in E is going to be realized independently with some probability p . The goal in the stochastic matching problem is to pick a sparse subgraph Q of G such that the realized edges in Q , in expectation, include a matching that is approximately as large as the maximum matching among the realized edges of G . The maximum degree of Q can depend on p , but not on the size of G . This problem has been subject to extensive studies over the years and the approximation factor has been improved from 0.5 to 0.5001 to 0.6568 and eventually to 2 3 . In this work, we analyze a natural sampling based algorithm and show that it can obtain all the way up to (1 epsilon) approximation, for any constant epsilon 0 . A key and of possible independent interest component of our analysis is an algorithm that constructs a matching on a stochastic graph, which among some other important properties, guarantees that each vertex is matched independently from the vertices that are sufficiently far. This allows us to bypass a previously known barrier towards achieving (1 epsilon) approximation based on existence of dense Ruzsa Szemer 'edi graphs.


Federated Learning for Emoji Prediction in a Mobile Keyboard. We show that a word level recurrent neural network can predict emoji from text typed on a mobile keyboard. We demonstrate the usefulness of transfer learning for predicting emoji by pretraining the model using a language modeling task. We also propose mechanisms to trigger emoji and tune the diversity of candidates. The model is trained using a distributed on device learning framework called federated learning. The federated model is shown to achieve better performance than a server trained model. This work demonstrates the feasibility of using federated learning to train production quality models for natural language understanding tasks while keeping users' data on their devices.


Discriminative Sleep Patterns of Alzheimer's Disease via Tensor Factorization. Sleep change is commonly reported in Alzheimer's disease (AD) patients and their brain wave studies show decrease in dreaming and non dreaming stages. Although sleep disturbance is generally considered as a consequence of AD, it might also be a risk factor of AD as new biological evidence shows. Leveraging National Sleep Research Resource (NSRR), we built a unique cohort of 83 cases and 331 controls with clinical variables and EEG signals. Supervised tensor factorization method was applied for this temporal dataset to extract discriminative sleep patterns. Among the 30 patterns extracted, we identified 5 significant patterns (4 patterns for AD likely and 1 pattern for normal ones) and their visual patterns provide interesting linkage to sleep with repeated wakefulness, insomnia, epileptic seizure, and etc. This study is preliminary but findings are interesting, which is a first step to provide quantifiable evidences to measure sleep as a risk factor of AD.


Improved regularity for the porous medium equation along zero level sets. In the present work we establish sharp regularity estimates for the solutions of the porous medium equation, along their zero level sets. We work under a proximity regime on the exponent governing the nonlinearity of the problem. Then, we prove that solutions are locally of class mathcal C 1 , frac 1 2 along free boundary points x 0 in partial left lbrace u 0 right rbrace , both in time and space. Our argument consists of importing information from the heat equation, through approximation and localization methods.


On proper holomorphic maps between bounded symmetric domains. We study proper holomorphic maps between bounded symmetric domains D and Omega . In particular, when D and Omega are of the same rank ge 2 such that all irreducible factors of D are of rank ge 2 , we prove that any proper holomorphic map from D to Omega is a totally geodesic holomorphic isometric embedding with respect to certain canonical K "ahler metrics of D and Omega . We also obtain some results regarding holomorphic maps F D to Omega which map minimal disks of D properly into rank 1 characteristic symmetric subspaces of Omega . On the other hand, we obtain new rigidity results regarding semi product proper holomorphic maps between D and Omega under a certain rank condition on D and Omega .


Peterzil Steinhorn subgroups and mu stabilizers in ACF. We consider G , a linear group defined over k , an algebraically closed field. By considering k as an embedded residue field of an algebraically closed valued field K , we can associate to it a compact G space S mu G(k) , consisting of mu types on G . We showed that for each p mu in S mu G(k) , text Stab mu(p) text Stab (p mu) is a solvable infinite algebraic group when p mu is centered at infinity and residually algebraic. Moreover we give a description of the dimension text Stab (p mu) in terms of dimension of p .


Compatible Paths on Labelled Point Sets. Let P and Q be finite point sets of the same cardinality in mathbb R 2 , each labelled from 1 to n . Two noncrossing geometric graphs G P and G Q spanning P and Q , respectively, are called compatible if for every face f in G P , there exists a corresponding face in G Q with the same clockwise ordering of the vertices on its boundary as in f . In particular, G P and G Q must be straight line embeddings of the same connected n vertex graph. Deciding whether two labelled point sets admit compatible geometric paths is known to be NP complete. We give polynomial time algorithms to find compatible paths or report that none exist in three scenarios O(n) time for points in convex position O(n 2) time for two simple polygons, where the paths are restricted to remain inside the closed polygons and O(n 2 log n) time for points in general position if the paths are restricted to be monotone.


Finite Time Performance of Distributed Two Time Scale Stochastic Approximation. Two time scale stochastic approximation is a popular iterative method for finding the solution of a system of two equations. Such methods have found broad applications in many areas, especially in machine learning and reinforcement learning. In this paper, we propose a distributed variant of this method over a network of agents, where the agents use two graphs representing their communication at different speeds due to the nature of their two time scale updates. Our main contribution is to provide a finite time analysis for the performance of the proposed method. In particular, we establish an upper bound for the convergence rates of the mean square errors at the agents to zero as a function of the step sizes and the network topology.


Incompressible limit for a two species tumour model with coupling through Brinkman's law in one dimension. We present a two species model with applications in tumour modelling. The main novelty is the coupling of both species through the so called Brinkman law which is typically used in the context of visco elastic media, where the velocity field is linked to the total population pressure via an elliptic equation. The same model for only one species has been studied by Perthame and Vauchelet in the past. The first part of this paper is dedicated to establishing existence of solutions to the problem, while the second part deals with the incompressible limit as the stiffness of the pressure law tends to infinity. Here we present a novel approach in one spatial dimension that differs from the kinetic reformulation used in the aforementioned study and, instead, relies on uniform BV estimates.


Moments of Dirichlet L functions with prime conductors over function fields. We compute the second moment in the family of quadratic Dirichlet L functions with prime conductors over mathbb F q x when the degree of the discriminant goes to infinity, obtaining one of the lower order terms. We also obtain an asymptotic formula with the leading order term for the mean value of the derivatives of L functions associated to quadratic twists of a fixed elliptic curve over mathbb F q(t) by monic irreducible polynomials, which allows us to show that there exists a monic irreducible polynomial such that the analytic rank of the corresponding twisted elliptic curve is equal to 1 .


The Ka Band High Power Klystron Amplifier Design Program of INFN. In the framework of the Compact Light XLS project, a short ultra high gradient linearizer working on the third harmonic of the main LINAC frequency is requested. Increasing gradients and reducing dimensions are requirements for XLS and all next generation linear accelerators. Actually, ultra compact normal conducting accelerating structures, operating in the Ka band regime ranging from 100 to 150 MV m are required to achieve ultra high gradients for research, industrial and medical applications. To fulfill these strong requirements, the R D of a proper Ka band klystron with RF power output and a high efficiency is mandatory. This contribution reports the design of a possible klystron amplifier tube that provides and output signal at 36 GHz to feed the phase space linearizer, while receiving at the input, a low level signal oscillating at 12 GHz which is the LINAC frequency. The proposed structure reaches an efficiency of 42 and a 20 MW RF power output. This contribution discusses also the high power DC gun, the beam focusing channel and the RF beam dynamics.


Convergence to a common fixed point of a finite family of nonexpansive mappings on the unit sphere of a Hilbert space. In this paper, we consider the CQ projection method and the shrinking projection method given by a finite family of nonexpansive mappings of the unit sphere of a real Hilbert space, and prove strong convergence theorems to common fixed points of the mappings.


Possibility to construct a clock on non 0 0 transition. In this letter, we discuss the possibility of non 0 0 atomic lines as the reference transitions to design a new atomic clock. The proposed clock operates in high C field regime offering an interesting alternative to conventional clock based on 0 0 transition as the former can also eliminate the influence of the first order Zeeman effect. Using the Breit Rabi formula, we theoretically calculate the magnetically insensitive frequency and the corresponding C field value for different alkali atom species. We also provide an estimate of the uncertainty of frequency shift to change of C field. We show that a high C field clock based on Cs 3, 1 4,0, 3,0 4, 1, and 3, 1 4, 1, can achieve a bias of 10 14 with a C field of 10 6 uncertainty. The free first order Zeeman effect C field for 3, 1 4,0 or 3,0 4, 1 transition is about one half that for 3, 1 4, 1 transition.


CalibreNet Calibration Networks for Multilingual Sequence Labeling. Lack of training data in low resource languages presents huge challenges to sequence labeling tasks such as named entity recognition (NER) and machine reading comprehension (MRC). One major obstacle is the errors on the boundary of predicted answers. To tackle this problem, we propose CalibreNet, which predicts answers in two steps. In the first step, any existing sequence labeling method can be adopted as a base model to generate an initial answer. In the second step, CalibreNet refines the boundary of the initial answer. To tackle the challenge of lack of training data in low resource languages, we dedicatedly develop a novel unsupervised phrase boundary recovery pre training task to enhance the multilingual boundary detection capability of CalibreNet. Experiments on two cross lingual benchmark datasets show that the proposed approach achieves SOTA results on zero shot cross lingual NER and MRC tasks.


Learning When to Trust a Dynamics Model for Planning in Reduced State Spaces. When the dynamics of a system are difficult to model and or time consuming to evaluate, such as in deformable object manipulation tasks, motion planning algorithms struggle to find feasible plans efficiently. Such problems are often reduced to state spaces where the dynamics are straightforward to model and evaluate. However, such reductions usually discard information about the system for the benefit of computational efficiency, leading to cases where the true and reduced dynamics disagree on the result of an action. This paper presents a formulation for planning in reduced state spaces that uses a classifier to bias the planner away from state action pairs that are not reliably feasible under the true dynamics. We present a method to generate and label data to train such a classifier, as well as an application of our framework to rope manipulation, where we use a Virtual Elastic Band (VEB) approximation to the true dynamics. Our experiments with rope manipulation demonstrate that the classifier significantly improves the success rate of our RRT based planner in several difficult scenarios which are designed to cause the VEB to produce incorrect predictions in key parts of the environment.


Independence over arbitrary sets in NSOP 1 theories. We study Kim independence over arbitrary sets. Assuming that forking satisfies existence, we establish Kim's lemma for Kim dividing over arbitrary sets in an NSOP 1 theory. We deduce symmetry of Kim independence and the independence theorem for Lascar strong types.


On Chip Concentration and Patterning of Biological Cells Using Interplay of Electrical and Thermal Fields. We demonstrate a method of concentrating and patterning of biological cells on a chip, exploiting the confluence of electric and thermal fields, without necessitating the use of any external heating or illuminating source. The technique simply employs two parallel plate electrodes and an insulating layer over the bottom electrode, with a drilled insulating layer for inducing localized variations in the thermal field. A strong induced electric field, in the process, penetrates through the narrow hole and generates highly non uniform heating, which in turn, results in gradients in electrical properties and induces mobile charges to impose directional fluid flow. The toroidal vortices, induced by secondary electrokinetic forces originating out of temperature dependent electrical property variations, transport the suspended cells towards a hot spot site of the chip, for rapid concentrating and patterning into different shaped clusters based on pre designed conditions, without exceeding safe temperature limits that do not result in damage of thermally labile biological samples. We characterize the efficacy of the cell trapping process for two different biological entities, namely, Escherichia coli bacteria and yeast cell. These results may be of profound importance towards developing novel biomedical microdevices for drug discovery, antibiotic resistance assessment and medical diagnostics.


A Benchmarking Study to Evaluate Apache Spark on Large Scale Supercomputers. As dataset sizes increase, data analysis tasks in high performance computing (HPC) are increasingly dependent on sophisticated dataflows and out of core methods for efficient system utilization. In addition, as HPC systems grow, memory access and data sharing are becoming performance bottlenecks. Cloud computing employs a data processing paradigm typically built on a loosely connected group of low cost computing nodes without relying upon shared storage and or memory. Apache Spark is a popular engine for large scale data analysis in the cloud, which we have successfully deployed via job submission scripts on production clusters. In this paper, we describe common parallel analysis dataflows for both Message Passing Interface (MPI) and cloud based applications. We developed an effective benchmark to measure the performance characteristics of these tasks using both types of systems, specifically comparing MPI C based analyses with Spark. The benchmark is a data processing pipeline representative of a typical analytics framework implemented using map reduce. In the case of Spark, we also consider whether language plays a role by writing tests using both Python and Scala, a language built on the Java Virtual Machine (JVM). We include performance results from two large systems at Argonne National Laboratory including Theta, a Cray XC40 supercomputer on which our experiments run with 65,536 cores (1024 nodes with 64 cores each). The results of our experiments are discussed in the context of their applicability to future HPC architectures. Beyond understanding performance, our work demonstrates that technologies such as Spark, while typically aimed at multi tenant cloud based environments, show promise for data analysis needs in a traditional clustering supercomputing environment.


A note on watchman's walks in de Bruijn graphs. The watchman's walk problem in a digraph calls for finding a minimum length closed dominating walk, where direction of arcs is respected. The watchman's walk of a de Bruijn graph of order k is described by a de Bruijn sequence of order k 1 . This idea is extended to certain subdigraphs of de Bruijn graphs.


A conditional proof of the ISP for quasinilpotent operators. The invariant subspace problem (ISP) is a well known unsolved problem in funtional analysis. While many partial results are known, the general case for complex, infinite dimensional separable Hilbert spaces is still open. It has been shown that the problem can be reduced to the case of operators which are norm limits of nilpotents. One of the most important subcases is the one of quasinilpotent operators, for which the problem has been extensively studied for many years. In this paper, we will introduce a new conjecture (supported by a heuristic argument), and we will prove conditionally that every quasinilpotent operator has a nontrivial invariant subspace. We will conclude by posing an open problem which would have deep implications regarding the ISP.


An Adaptive Parareal Algorithm. In this paper, we consider the problem of accelerating the numerical simulation of time dependent problems by time domain decomposition. The available algorithms enabling such decompositions present severe efficiency limitations and are an obstacle for the solution of large scale and high dimensional problems. Our main contribution is the improvement of the parallel efficiency of the parareal in time method. The parareal method is based on combining predictions made by a numerically inexpensive solver (with coarse physics and or coarse resolution) with corrections coming from an expensive solver (with high fidelity physics and high resolution). At convergence, the algorithm provides a solution that has the fine solver's high fidelity physics and high resolution. In the classical version, the fine solver has a fixed high accuracy which is the major obstacle to achieve a competitive parallel efficiency. In this paper, we develop an adaptive variant that overcomes this obstacle by dynamically increasing the accuracy of the fine solver across the parareal iterations. We theoretically show that the parallel efficiency becomes very competitive in the ideal case where the cost of the coarse solver is small, thus proving that the only remaining factors impeding full scalability become the cost of the coarse solver and communication time. The developed theory has also the merit of setting a general framework to understand the success of several extensions of parareal based on iteratively improving the quality of the fine solver and re using information from previous parareal steps. We illustrate the actual performance of the method in stiff ODEs, which are a challenging family of problems since the only mechanism for adaptivity is time and efficiency is affected by the cost of the coarse solver.


Conjunctive Queries with Theta Joins Under Updates. Modern application domains such as Composite Event Recognition (CER) and real time Analytics require the ability to dynamically refresh query results under high update rates. Traditional approaches to this problem are based either on the materialization of subresults (to avoid their recomputation) or on the recomputation of subresults (to avoid the space overhead of materialization). Both techniques have recently been shown suboptimal instead of materializing results and subresults, one can maintain a data structure that supports efficient maintenance under updates and can quickly enumerate the full query output, as well as the changes produced under single updates. Unfortunately, these data structures have been developed only for aggregate join queries composed of equi joins, limiting their applicability in domains such as CER where temporal joins are commonplace. In this paper, we present a new approach for dynamically evaluating queries with multi way theta joins under updates that is effective in avoiding both materialization and recomputation of results, while supporting a wide range of applications. To do this we generalize Dynamic Yannakakis, an algorithm for dynamically processing acyclic equi join queries. In tandem, and of independent interest, we generalize the notions of acyclicity and free connexity to arbitrary theta joins and show how to compute corresponding join trees. We instantiate our framework to the case where theta joins are only composed of equalities and inequalities and experimentally compare our algorithm to state of the art CER systems as well as incremental view maintenance engines. Our approach performs consistently better than the competitor systems with up to two orders of magnitude improvements in both time and memory consumption.


Structural Color 3D Printing By Shrinking Photonic Crystals. The rings, spots and stripes found on some butterflies, Pachyrhynchus weevils, and many chameleons are notable examples of natural organisms employing photonic crystals to produce colorful patterns. Despite advances in nanotechnology, we still lack the ability to print arbitrary colors and shapes in all three dimensions at this microscopic length scale. Commercial nanoscale 3D printers based on two photon polymerization are incapable of patterning photonic crystal structures with the requisite 300 nm lattice constant to achieve photonic stopbands bandgaps in the visible spectrum and generate colors. Here, we introduce a means to produce 3D printed photonic crystals with a 5x reduction in lattice constants (periodicity as small as 280 nm), achieving sub 100 nm features with a full range of colors. The reliability of this process enables us to engineer the bandstructures of woodpile photonic crystals that match experiments, showing that observed colors can be attributed to either slow light modes or stopbands. With these lattice structures as 3D color volumetric elements (voxels), we printed 3D microscopic scale objects, including the first multi color microscopic model of the Eiffel Tower measuring only 39 microns tall with a color pixel size of 1.45 microns. The technology to print 3D structures in color at the microscopic scale promises the direct patterning and integration of spectrally selective devices, such as photonic crystal based color filters, onto free form optical elements and curved surfaces.


Improved Bounds for (b,k) hashing. For fixed integers b geq k , a problem of relevant interest in computer science and combinatorics is that of determining the asymptotic growth, with n , of the largest set for which a (b, k) hash family of n functions exists. Equivalently, determining the asymptotic growth of a largest subset of 1,2, ldots,b n such that, for any k distinct elements in the set, there is a coordinate where they all differ. An important asymptotic upper bound for general b, k , was derived by Fredman and Koml 'os in the '80s and improved for certain b neq k by K "orner and Marton and by Arikan. Only very recently better bounds were derived for the general b,k case by Guruswami and Riazanov while stronger results for small values of b k were obtained by Arikan, by Dalai, Guruswami and Radhakrishnan and by Costa and Dalai. In this paper, we both show how some of the latter results extend to b neq k and further strengthen the bounds for some specific small values of b and k . The method we use, which depends on the reduction of an optimization problem to a finite number of cases, shows that further results might be obtained by refined arguments at the expense of higher complexity which could be reduced by using more sophisticated and optimized algorithmic approaches.


Solving Partial Differential Equations on Closed Surfaces with Planar Cartesian Grids. We present a general purpose method for solving partial differential equations on a closed surface, based on a technique for discretizing the surface introduced by Wenjun Ying and Wei Cheng Wang J. Comput. Phys. 252 (2013), pp. 606 624 which uses projections on coordinate planes. Assuming it is given as a level set, the surface is represented by a set of points at which it intersects the intervals between grid points in a three dimensional grid. They are designated as primary or secondary. Discrete functions on the surface have independent values at primary points, with values at secondary points determined by an equilibration process. Each primary point and its neighbors have projections to regular grid points in a coordinate plane where the equilibration is done and finite differences are computed. The solution of a p.d.e. can be reduced to standard methods on Cartesian grids in the coordinate planes, with the equilibration allowing seamless transition from one system to another. We observe second order accuracy in examples with a variety of equations, including surface diffusion determined by the Laplace Beltrami operator and the shallow water equations on a sphere.


Attentive Weakly Supervised land cover mapping for object based satellite image time series data with spatial interpretation. Nowadays, modern Earth Observation systems continuously collect massive amounts of satellite information. The unprecedented possibility to acquire high resolution Satellite Image Time Series (SITS) data (series of images with high revisit time period on the same geographical area) is opening new opportunities to monitor the different aspects of the Earth Surface but, at the same time, it is raising up new challenges in term of suitable methods to analyze and exploit such huge amount of rich and complex image data. One of the main task associated to SITS data analysis is related to land cover mapping where satellite data are exploited via learning methods to recover the Earth Surface status aka the corresponding land cover classes. Due to operational constraints, the collected label information, on which machine learning strategies are trained, is often limited in volume and obtained at coarse granularity carrying out inexact and weak knowledge that can affect the whole process. To cope with such issues, in the context of object based SITS land cover mapping, we propose a new deep learning framework, named TASSEL (aTtentive weAkly Supervised Satellite image time sEries cLassifier), that is able to intelligently exploit the weak supervision provided by the coarse granularity labels. Furthermore, our framework also produces an additional side information that supports the model interpretability with the aim to make the black box gray. Such side information allows to associate spatial interpretation to the model decision via visual inspection.


Perturbation theories for symmetry protected bound states in the continuum on two dimensional periodic structures. On dielectric periodic structures with a reflection symmetry in a periodic direction, there can be antisymmetric standing waves (ASWs) that are symmetry protected bound states in the continuum (BICs). The BICs have found many applications, mainly because they give rise to resonant modes of extremely large quality factors ( Q factors). The ASWs are robust to symmetric perturbations of the structure, but they become resonant modes if the perturbation is non symmetric. The Q factor of a resonant mode on a perturbed structure is typically O(1 delta 2) where delta is the amplitude of the perturbation, but special perturbations can produce resonant modes with larger Q factors. For two dimensional (2D) periodic structures with a 1D periodicity, we derive conditions on the perturbation profile such that the Q factors are O(1 delta 4) or O(1 delta 6) . For the unperturbed structure, an ASW is surrounded by resonant modes with a nonzero Bloch wave vector. For 2D periodic structures, the Q factors of nearby resonant modes are typically O(1 beta 2) , where beta is the Bloch wavenumber. We show that the Q factors can be O(1 beta 6) if the ASW satisfies a simple condition.


A closed linkage mechanism having the shape of a discrete M "obius strip. A closed linkage mechanism in three dimensional space is an object comprising rigid bodies connected with hinges in a circular form like a rosary. Such linkages include Bricard6R and Bennett4R. To design such a closed linkage, it is necessary to solve a high degree algebraic equation, which is generally difficult. In this lecture, the author proposes a new family of closed linkage mechanisms with an arbitrary number of hinges as an extension of a certain Bricard6R. They have singular properties, such as one dimensional degree of freedom (1 DOF), and certain energies taking a constant value regardless of the state. These linkage mechanisms can be regarded as discrete M "obius strips and may be of interest in the context of pure mathematics as well. However, many of the properties described here have been confirmed only numerically, with no rigorous mathematical proof, and should be interpreted with caution.


Predicting risk of dyslexia with an online gamified test. Dyslexia is a specific learning disorder related to school failure. Detection is both crucial and challenging, especially in languages with transparent orthographies, such as Spanish. To make detecting dyslexia easier, we designed an online gamified test and a predictive machine learning model. In a study with more than 3,600 participants, our model correctly detected over 80 of the participants with dyslexia. To check the robustness of the method we tested our method using a new data set with over 1,300 participants with age customized tests in a different environment a tablet instead of a desktop computer reaching a recall of over 72 for the class with dyslexia for children 9 years old or older. Our work shows that dyslexia can be screened using a machine learning approach. An online screening tool based on our methods has already been used by more than 200,000 people.


Hierarchical Deep Learning of Multiscale Differential Equation Time Steppers. Nonlinear differential equations rarely admit closed form solutions, thus requiring numerical time stepping algorithms to approximate solutions. Further, many systems characterized by multiscale physics exhibit dynamics over a vast range of timescales, making numerical integration computationally expensive due to numerical stiffness. In this work, we develop a hierarchy of deep neural network time steppers to approximate the flow map of the dynamical system over a disparate range of time scales. The resulting model is purely data driven and leverages features of the multiscale dynamics, enabling numerical integration and forecasting that is both accurate and highly efficient. Moreover, similar ideas can be used to couple neural network based models with classical numerical time steppers. Our multiscale hierarchical time stepping scheme provides important advantages over current time stepping algorithms, including (i) circumventing numerical stiffness due to disparate time scales, (ii) improved accuracy in comparison with leading neural network architectures, (iii) efficiency in long time simulation forecasting due to explicit training of slow time scale dynamics, and (iv) a flexible framework that is parallelizable and may be integrated with standard numerical time stepping algorithms. The method is demonstrated on a wide range of nonlinear dynamical systems, including the Van der Pol oscillator, the Lorenz system, the Kuramoto Sivashinsky equation, and fluid flow pass a cylinder audio and video signals are also explored. On the sequence generation examples, we benchmark our algorithm against state of the art methods, such as LSTM, reservoir computing, and clockwork RNN. Despite the structural simplicity of our method, it outperforms competing methods on numerical integration.


Construct b symplectic toric manifolds from toric manifolds. In cite btoric , Guillemin et al. proved a Delzant type theorem which classifies b symplectic toric manifolds. More generally, in cite torus they proved a similar convexity result for general Hamiltonian torus action on b symplectic manifolds. In this paper, we provide a new way to construct b symplectic toric manifolds from usual toric manifolds. Conversely, through this way, we can also decompose a b symplectic toric manifolds to usual toric manifolds. Finally, we will try to prove that this kind of decomposition is useful, although the symplectic structure for our decomposition or construction is not canonical.


Some Inequalities for Continuous Algebra Multiplications on a Banach Space. In this short note, we first consider some inequalities for comparison of some algebraic properties of two continuous algebra multiplications on an arbitrary Banach space and then, as an application, we consider some very basic observations on the space of all continuous algebra multiplications on a Banach space.


Efficient Analysis of Test beam Data with the Corryvreckan Framework. Stringent requirements are posed on the the next generations of vertex and tracking detectors for high energy physics experiments to reach the foreseen physics goals. A large variety of silicon pixel sensors targeting the specific needs of each use case are developed and tested both in laboratory and test beam measurement campaigns. Corryvreckan is a flexible, fast and lightweight test beam data reconstruction and analysis framework based on a modular concept of the reconstruction chain. It is designed to fulfil the requirements for offline event building in complex data taking environments combining detectors with different readout schemes. Its modular architecture separates the framework core from the implementation of reconstruction, analysis and detector specific algorithms. In this paper, a brief overview of the software framework and the reconstruction and analysis chain is provided. This is complemented by an example analysis of a data set using the offline event building capabilities of the framework and an improved event building scheme allowing for a more efficient usage of test beam data exploiting the pivot pixel information of the Mimosa26 sensors.


Phase transitions for spatially extended pinning. We consider a directed polymer of length N interacting with a linear interface. The monomers carry i.i.d. random charges ( omega i) i 1 N taking values in mathbb R with mean zero and variance one. Each monomer i contributes an energy ( beta omega i h) varphi(S i) to the interaction Hamiltonian, where S i in mathbb Z is the height of monomer i with respect to the interface, varphi mathbb Z to 0, infty) is the interaction potential, beta in 0, infty) is the inverse temperature, and h in mathbb R is the charge bias parameter. The configurations of the polymer are weighted according to the Gibbs measure associated with the interaction Hamiltonian, where the reference measure is given by a Markov chain on mathbb Z . We study both the quenched and the annealed free energy per monomer in the limit as N to infty . We show that each exhibits a phase transition along a critical curve in the ( beta, h) plane, separating a localized phase (where the polymer stays close to the interface) from a delocalized phase (where the polymer wanders away from the interface). We derive variational formulas for the critical curves and we obtain upper and lower bounds on the quenched critical curve in terms of the annealed critical curve. In addition, for the special case where the reference measure is given by a Bessel random walk, we derive the scaling limit of the annealed free energy as beta, h downarrow 0 in three different regimes for the tail exponent of varphi .


Solving a Special Case of the Intensional vs Extensional Conjecture in Probabilistic Databases. We consider the problem of exact probabilistic inference for Union of Conjunctive Queries (UCQs) on tuple independent databases. For this problem, two approaches currently coexist. In the extensional method, query evaluation is performed by exploiting the structure of the query, and relies heavily on the use of the inclusion exclusion principle. In the intensional method, one first builds a representation of the lineage of the query in a tractable formalism of knowledge compilation. The chosen formalism should then ensure that the probability can be efficiently computed using simple disjointness and independence assumptions, without the need of performing inclusion exclusion. The extensional approach has long been thought to be strictly more powerful than the intensional approach, the reason being that for some queries, the use of inclusion exclusion seemed unavoidable. In this paper we introduce a new technique to construct lineage representations as deterministic decomposable circuits in polynomial time. We prove that this technique applies to a class of UCQs that had been conjectured to separate the complexity of the two approaches. In essence, we show that relying on the inclusion exclusion formula can be avoided by using negation. This result brings back hope to prove that the intensional approach can handle all tractable UCQs.


An AGI with Time Inconsistent Preferences. This paper reveals a trap for artificial general intelligence (AGI) theorists who use economists' standard method of discounting. This trap is implicitly and falsely assuming that a rational AGI would have time consistent preferences. An agent with time inconsistent preferences knows that its future self will disagree with its current self concerning intertemporal decision making. Such an agent cannot automatically trust its future self to carry out plans that its current self considers optimal.


Single use register automata for data words. Our starting point are register automata for data words, in the style of Kaminski and Francez. We study the effects of the single use restriction, which says that a register is emptied immediately after being used. We show that under the single use restriction, the theory of automata for data words becomes much more robust. The main results are (a) five different machine models are equivalent as language acceptors, including one way and two way single use register automata (b) one can recover some of the algebraic theory of languages over finite alphabets, including a version of the Krohn Rhodes Theorem (c) there is also a robust theory of transducers, with four equivalent models, including two way single use transducers and a variant of streaming string transducers for data words. These results are in contrast with automata for data words without the single use restriction, where essentially all models are pairwise non equivalent.


A curvature flow approach to Lp Christoffel Minkowski problem for p 1. We study the motion of smooth, closed, strictly convex hypersurfaces in Rn 1 expanding in the direction of their normal vector field with speed depending on the k th elementary symmetric polynomial of the principal radii of curvature. As an application, we gives a unifed fow approach to Lp Christoffel Minkowski problem for p 1.


Plug and Play Anomaly Detection with Expectation Maximization Filtering. Anomaly detection in crowds enables early rescue response. A plug and play smart camera for crowd surveillance has numerous constraints different from typical anomaly detection the training data cannot be used iteratively there are no training labels and training and classification needs to be performed simultaneously. We tackle all these constraints with our approach in this paper. We propose a Core Anomaly Detection (CAD) neural network which learns the motion behavior of objects in the scene with an unsupervised method. On average over standard datasets, CAD with a single epoch of training shows a percentage increase in Area Under the Curve (AUC) of 4.66 and 4.9 compared to the best results with convolutional autoencoders and convolutional LSTM based methods, respectively. With a single epoch of training, our method improves the AUC by 8.03 compared to the convolutional LSTM based approach. We also propose an Expectation Maximization filter which chooses samples for training the core anomaly detection network. The overall framework improves the AUC compared to future frame prediction based approach by 24.87 when crowd anomaly detection is performed on a video stream. We believe our work is the first step towards using deep learning methods with autonomous plug and play smart cameras for crowd anomaly detection.


Joint CLT for top eigenvalues of sample covariance matrices of separable high dimensional long memory processes. For N,n in mathbb N , consider the sample covariance matrix S N(T) frac 1 N XX from a data set X C N 1 2 ZT n 1 2 , where Z (Z i,j ) is a N times n matrix having i.i.d. entries with mean zero and variance one, and C N, T n are deterministic positive semi definite Hermitian matrices of dimension N and n , respectively. We assume that (C N) N is bounded in spectral norm, and T n is a Toeplitz matrix with its largest eigenvalues diverging to infinity. The matrix X can be viewed as a data set of an N dimensional long memory stationary process having separable dependence structure. As N,n to infty and Nn 1 to r in (0, infty) , we establish the asymptotics and the joint CLT for ( lambda 1(S N(T)), cdots, lambda m(S N(T)) top where lambda j(S N(T)) denotes the j th largest eigenvalue of S N(T) , and m is a fixed integer. For the CLT, we first study the case where the entries of Z are Gaussian, and then we generalize the result to some more generic cases. This result substantially extends our previous result in Merlev ede et al. 2019, where we studied lambda 1(S N(T)) in the case where m 1 and X ZT n 1 2 with Z having Gaussian entries. In order to establish this CLT, we are led to study the first order asymptotics of the largest eigenvalues and the associated eigenvectors of some deterministic Toeplitz matrices. We are specially interested in the autocovariance matrices of long memory stationary processes. We prove multiple spectral gap properties for the largest eigenvalues and a delocalization property for their associated eigenvectors.


Algebraic commutators with respect to subnormal subgroups in division rings. Let D be a division ring and K a subfield of D which is not necessarily contained in the center F of D . In this paper, we study the structure of D under the condition of left algebraicity of certain subsets of D over K . Among results, it is proved that if D contains a noncentral normal subgroup which is left algebraic over K of bounded degree d , then D F le d 2 . In case K F , the obtained results show that if either all additive commutators or all multiplicative commutators with respect to a noncentral subnormal subgroup of D are algebraic of bounded degree d over F , then D F le d 2 .


Probabilistic Software Modeling A Data driven Paradigm for Software Analysis. Software systems are complex, and behavioral comprehension with the increasing amount of AI components challenges traditional testing and maintenance strategies.The lack of tools and methodologies for behavioral software comprehension leaves developers to testing and debugging that work in the boundaries of known scenarios. We present Probabilistic Software Modeling (PSM), a data driven modeling paradigm for predictive and generative methods in software engineering. PSM analyzes a program and synthesizes a network of probabilistic models that can simulate and quantify the original program's behavior. The approach extracts the type, executable, and property structure of a program and copies its topology. Each model is then optimized towards the observed runtime leading to a network that reflects the system's structure and behavior. The resulting network allows for the full spectrum of statistical inferential analysis with which rich predictive and generative applications can be built. Applications range from the visualization of states, inferential queries, test case generation, and anomaly detection up to the stochastic execution of the modeled system. In this work, we present the modeling methodologies, an empirical study of the runtime behavior of software systems, and a comprehensive study on PSM modeled systems. Results indicate that PSM is a solid foundation for structural and behavioral software comprehension applications.


Powered Multifocal THz Radiation by Mixing of Two Skew (ChG) Cosh Gaussian Laser Beams in Collisional Plasma. Considering a realistic situation, where electron neutral collisions persist in plasma, analytical calculations are carried out for the THz radiation generation by beating of two Skew Cosh Gaussian (SG) lasers of parameters n s. The competency of these lasers over Gaussian lasers is discussed in detail with respect to the effects of collision and beam width on the THz field amplitude and efficiency of the mechanism. A critical transverse distance of the peak of the THz field is defined that shows a dependence on the skewness parameter of skew ChG lasers. Although electron neutral collisions and larger beam width lead to the drastic reduction in the THz field when the skew ChG lasers are used in the plasma, the efficiency of the mechanism remains much larger than the case of Gaussian lasers. Moreover, the skewChG lasers produce stronger and multifocal THz radiation.


Perturbation and Stability of Continuous Operator Frames in Hilbert C ast Modules. Frame Theory has a great revolution in recent years. This Theory have been extended from Hilbert spaces to Hilbert C ast modules. In this paper we consider the stability of continuous operator frame and continuous K operator frames in Hilbert C ast Modules under perturbation and we establish some properties


Global solvability criteria for some classes of nonlinear second order ordinary differential equations. The Riccati equation method is used to establish some global solvability criteria for some classes of second order nonlinear ordinary differential equations. Two oscillation theorems are proved. The results are applied to the Emden Fowler equation and to the Van der Pol type equation.


Sparse Noncommutative Polynomial Optimization. This article focuses on optimization of polynomials in noncommuting variables, while taking into account sparsity in the input data. A converging hierarchy of semidefinite relaxations for eigenvalue and trace optimization is provided. This hierarchy is a noncommutative analogue of results due to Lasserre SIAM J. Optim. 17(3) (2006), pp. 822 843 and Waki et al. SIAM J. Optim. 17(1) (2006), pp. 218 242 . The Gelfand Naimark Segal (GNS) construction is applied to extract optimizers if flatness and irreducibility conditions are satisfied. Among the main techniques used are amalgamation results from operator algebra. The theoretical results are utilized to compute lower bounds on minimal eigenvalue and trace of noncommutative polynomials from the literature.


TransForm Formally Specifying Transistency Models and Synthesizing Enhanced Litmus Tests. Memory consistency models (MCMs) specify the legal ordering and visibility of shared memory accesses in a parallel program. Traditionally, instruction set architecture (ISA) MCMs assume that relevant program visible memory ordering behaviors only result from shared memory interactions that take place between user level program instructions. This assumption fails to account for virtual memory (VM) implementations that may result in additional shared memory interactions between user level program instructions and both 1) system level operations (e.g., address remappings and translation lookaside buffer invalidations initiated by system calls) and 2) hardware level operations (e.g., hardware page table walks and dirty bit updates) during a user level program's execution. These additional shared memory interactions can impact the observable memory ordering behaviors of user level programs. Thus, memory transistency models (MTMs) have been coined as a superset of MCMs to additionally articulate VM aware consistency rules. However, no prior work has enabled formal MTM specifications, nor methods to support their automated analysis. To fill the above gap, this paper presents the TransForm framework. First, TransForm features an axiomatic vocabulary for formally specifying MTMs. Second, TransForm includes a synthesis engine to support the automated generation of litmus tests enhanced with MTM features (i.e., enhanced litmus tests, or ELTs) when supplied with a TransForm MTM specification. As a case study, we formally define an estimated MTM for Intel x86 processors, called x86t elt, that is based on observations made by an ELT based evaluation of an Intel x86 MTM implementation from prior work and available public documentation. Given x86t elt and a synthesis bound as input, TransForm's synthesis engine successfully produces a set of ELTs including relevant ELTs from prior work.


A 1 connected components of ruled surfaces. A conjecture of Morel asserts that the sheaf of mathbb A 1 connected components of a space is mathbb A 1 invariant. Using purely algebro geometric methods, we determine the sheaf of mathbb A 1 connected components of a smooth projective surface, which is birationally ruled over a curve of genus 0 . As a consequence, we show that Morel's conjecture holds for all smooth projective surfaces over an algebraically closed field of characteristic 0 .


Automatic Video Object Segmentation via Motion Appearance Stream Fusion and Instance aware Segmentation. This paper presents a method for automatic video object segmentation based on the fusion of motion stream, appearance stream, and instance aware segmentation. The proposed scheme consists of a two stream fusion network and an instance segmentation network. The two stream fusion network again consists of motion and appearance stream networks, which extract long term temporal and spatial information, respectively. Unlike the existing two stream fusion methods, the proposed fusion network blends the two streams at the original resolution for obtaining accurate segmentation boundary. We develop a recurrent bidirectional multiscale structure with skip connection for the stream fusion network to extract long term temporal information. Also, the multiscale structure enables to obtain the original resolution features at the end of the network. As a result of two stream fusion, we have a pixel level probabilistic segmentation map, which has higher values at the pixels belonging to the foreground object. By combining the probability of foreground map and objectness score of instance segmentation mask, we finally obtain foreground segmentation results for video sequences without any user intervention, i.e., we achieve successful automatic video segmentation. The proposed structure shows a state of the art performance for automatic video object segmentation task, and also achieves near semi supervised performance.


The Fraction of Broken Waves in Natural Surf Zones. This paper presents a novel quantification of the fraction of broken waves (Qb) in natural surf zones using data from seven wave dominated Australian beaches. Qb is a critical, but rarely quantified, parameter for parametric surf zone energy dissipation models which are commonly used as coastal management tools. Here, Qb is quantified using a combination of remote sensing and in situ data. These data and machine learning techniques enable quantification of Qb for a substantial dataset ( 350,000 waves). The results show that Qb is a highly variable parameter with a high degree of inter and intra beach variability. Such variance could be explained (at least partially) by correlations between Qb and environmental parameters. Tidal variations drive changes in Qb of up to 70 for a given local water depth (h) on steep beaches, and increased infragravity energy levels decreased terminal values of Qb by about 20 . The links between Qb and environmental forcing lead to the development of a correspondence between Qb and the Australian beach morphodynamic model. Qb is larger for a given normalized depth (h Hs, where Hs is offshore wave height) for dissipative beaches than for intermediate beaches. Finally, when comparing data to existing models, three commonly used theoretical formulations for Qb are observed to be poor predictors with errors of the order of 40 . Existing theoretical Qb models are shown to improve (revised errors of the order of 10 ) if the Rayleigh probability distribution that describes the wave height is in these models is replaced by the Weibull distribution.


Spontaneous and Stimulated Electron Photon Interactions in Nanoscale Plasmonic Near Fields. We demonstrate spatially resolved measurements of spontaneous and stimulated electron photon interactions in nanoscale optical near fields using electron energy loss spectroscopy (EELS), cathodoluminescence spectroscopy (CL), and photon induced near field electron microscopy (PINEM). Specifically, we study resonant surface plasmon modes that are tightly confined to the tip apexes of an Au nanostar, enabling a direct correlation of EELS, CL, and PINEM on the same physical structure at the nanometer length scale. Complemented by numerical electromagnetic boundary element method calculations, we discuss the spontaneous and stimulated electron photon interaction strength and spatial dependence of our EELS, CL and PINEM distributions. We demonstrate that in the limit of an isolated tip mode, spatial variations in the electron near field coupling are fully determined by the modal electric field profile, irrespective of the spontaneous (in EELS and CL) or stimulated nature (in PINEM) of the process. Yet we show that coupling to the tip modes crucially depends on the incident electron energy with a maximum at a few keV, depending on the proximity of the interaction to the tip apex. Our results provide elementary insights into spontaneous and stimulated electron light matter interactions at the nanoscale that have key implications for research on (quantum) coherent optical phenomena in electron microscopy.


Some partition and analytical identities arising from Alladi, Andrews, Gordon bijections. In a work of 1995, Alladi, Andrews, and Gordon provided a generalization of the two Capparelli identities involving certain classes of integer partitions. Inspired by that contribution, in particular as regards the general setting and the tools the authors employed, we obtain new partition identities by identifying further sets of partitions that can be explicitly put into a one to one correspondence by the method described in the 1995 paper. As a further result, although of a different nature, we obtain an analytical identity of Rogers Ramanujan type, involving generating functions, for a class of partition identities already found in that paper and that generalize the first Capparelli identity and include it as a particular case. To achieve this, we apply the same strategy as Kanade and Russel did in a recent paper. This method relies on the use of jagged partitions that can be seen as a more general kind of integer partitions.


Average tempered stable subordinators with applications. In this paper the running average of a subordinator with a tempered stable distribution is considered. We investigate a family of previously unexplored infinite activity subordinators induced by the probability distribution of the running average process and determine their jump intensity measures. Special cases including gamma processes and inverse Gaussian processes are discussed. Then we derive easily implementable formulas for the distribution functions, cumulants, and moments, as well as provide explicit estimates for their asymptotic behaviors. Numerical experiments are conducted for illustrating the applicability and efficiency of the proposed formulas. Two important extensions of the running average process and its equi distributed subordinator are examined with concrete applications to structural degradation modeling and financial derivatives pricing, where their advantages relative to several existing models are highlighted together with the mention of Euler discretization and compound Poisson approximation techniques.


SMART Paths for Latency Reduction in ReRAM Processing In Memory Architecture for CNN Inference. This research work proposes a design of an analog ReRAM based PIM (processing in memory) architecture for fast and efficient CNN (convolutional neural network) inference. For the overall architecture, we use the basic hardware hierarchy such as node, tile, core, and subarray. On the top of that, we design intra layer pipelining, inter layer pipelining, and batch pipelining to exploit parallelism in the architecture and increase overall throughput for the inference of an input image stream. We also optimize the performance of the NoC (network on chip) routers by decreasing hop counts using SMART (single cycle multi hop asynchronous repeated traversal) flow control. Finally, we experiment with weight replications for different CNN layers in VGG (A E) for large scale data set ImageNet. In our simulation, we achieve 40.4027 TOPS (tera operations per second) for the best case performance, which corresponds to over 1029 FPS (frames per second). We also achieve 3.5914 TOPS W (tera operaions per second per watt) for the best case energy efficiency. In addition, the architecture with aggressive pipelining and weight replications can achieve 14X speedup compared to the baseline architecture with basic pipelining, and SMART flow control achieves 1.08X speedup in this architecture compared to the baseline. Last but not least, we also evaluate the performance of SMART flow control using synthetic traffic.


Hyperbolic metrics on open subsests of Ptolemaic spaces with sharp parameter bounds. It is shown that a construction of Z. Zhang and Y. Xiao on open subsets of Ptolemaic spaces yields, when the subset has boundary containing at least two points, metrics that are Gromov hyperbolic with parameter log 2 and strongly hyperbolic with parameter 1 with no further conditions on the open set. A class of examples is constructed on Hadamard manifolds showing these estimates of the parameters are sharp.


Vertex distinction with subgraph centrality a proof of Estrada's conjecture and some generalizations. Centrality measures are used in network science to identify the most important vertices for transmission of information and dynamics on a graph. One of these measures, introduced by Estrada and collaborators, is the beta subgraph centrality, which is based on the exponential of the matrix beta A , where A is the adjacency matrix of the graph and beta is a real parameter ("inverse temperature"). We prove that for algebraic beta , two vertices with equal beta subgraph centrality are necessarily cospectral. We further show that two such vertices must have the same degree and eigenvector centralities. Our results settle a conjecture of Estrada and a generalization of it due to Kloster, Kr 'al and Sullivan. We also discuss possible extensions of our results.


On the Coexistence of Stability and Incentive Compatibility in Fractional Matchings. Stable matchings have been studied extensively in social choice literature. The focus has been mostly on integral matchings, in which the nodes on the two sides are wholly matched. A fractional matching, which is a convex combination of integral matchings, is a natural extension of integral matchings. The topic of stability of fractional matchings has started receiving attention only very recently. Further, incentive compatibility in the context of fractional matchings has received very little attention. With this as the backdrop, our paper studies the important topic of incentive compatibility of mechanisms to find stable fractional matchings. We work with preferences expressed in the form of cardinal utilities. Our first result is an impossibility result that there are matching instances for which no mechanism that produces a stable fractional matching can be incentive compatible or even approximately incentive compatible. This provides the motivation to seek special classes of matching instances for which there exist incentive compatible mechanisms that produce stable fractional matchings. Our study leads to a class of matching instances that admit unique stable fractional matchings. We first show that a unique stable fractional matching for a matching instance exists if and only if the given matching instance satisfies the conditional mutual first preference (CMFP) property. To this end, we provide a polynomial time algorithm that makes ingenious use of envy graphs to find a non integral stable matching whenever the preferences are strict and the given instance is not a CMFP matching instance. For this class of CMFP matching instances, we prove that every mechanism that produces the unique stable fractional matching is (a) incentive compatible and further (b) resistant to coalitional manipulations.


Atoms of the matching measure. We prove that the matching measure of an infinite vertex transitive connected graph has no atoms. Generalizing the results of Salez, we show that for an ergodic non amenable unimodular random rooted graph with uniformly bounded degrees, the matching measure has only finitely many atoms. Ku and Chen proved the analogue of the Gallai Edmonds structure theorem for non zero roots of the matching polynomial for finite graphs. We extend their results for infinite graphs. We also show that the corresponding Gallai Edmonds decomposition is compatible with the zero temperature monomer dimer model.


On the number of alternating paths in random graphs. In the noisy channel model from coding theory, we wish to detect errors introduced during transmission by optimizing various parameters of the code. Bennett, Dudek, and LaForge framed a variation of this problem in the language of alternating paths in edge colored complete bipartite graphs in 2016. Here, we extend this problem to the random graph mathbb G (n,p) . We seek the alternating connectivity, kappa r, ell (G) , which is the maximum t such that there is an r edge coloring of G such that any pair of vertices is connected by t internally disjoint and alternating (i.e. no consecutive edges of the same color) paths of length ell . We have three main results about how this parameter behaves in mathbb G (n,p) that basically cover all ranges of p one for paths of length two, one for the dense case, and one for the sparse case. For paths of length two, we found that kappa r, ell (G) is essentially the codegree of a pair of vertices. For the dense case when p is constant, we were able to achieve the natural upper bounds of minimum degree (minus some intersection) or the total number of disjoint paths between a pair of vertices. For the sparse case, we were able to find colorings that achieved the natural obstructions of minimum degree or (in a slightly less precise result) the total number of paths of a certain length in a graph. We broke up this sparse case into ranges of p corresponding to when mathbb G (n,p) has diameter k or k 1 . We close with some remarks about a similar parameter and a generalization to pseudorandom graphs.


Cloth in the Wind A Case Study of Physical Measurement through Simulation. For many of the physical phenomena around us, we have developed sophisticated models explaining their behavior. Nevertheless, measuring physical properties from visual observations is challenging due to the high number of causally underlying physical parameters including material properties and external forces. In this paper, we propose to measure latent physical properties for cloth in the wind without ever having seen a real example before. Our solution is an iterative refinement procedure with simulation at its core. The algorithm gradually updates the physical model parameters by running a simulation of the observed phenomenon and comparing the current simulation to a real world observation. The correspondence is measured using an embedding function that maps physically similar examples to nearby points. We consider a case study of cloth in the wind, with curling flags as our leading example a seemingly simple phenomena but physically highly involved. Based on the physics of cloth and its visual manifestation, we propose an instantiation of the embedding function. For this mapping, modeled as a deep network, we introduce a spectral layer that decomposes a video volume into its temporal spectral power and corresponding frequencies. Our experiments demonstrate that the proposed method compares favorably to prior work on the task of measuring cloth material properties and external wind force from a real world video.


On the number of generators of an algebra over a commutative ring. A theorem of O. Forster says that if R is a noetherian ring of Krull dimension d , then any projective R module of rank n can be generated by d n elements. S. Chase and R. Swan subsequently showed that this bound is sharp there exist examples that cannot be generated by fewer than d n elements. We view projective R modules as R forms of the non unital R algebra where the product of any two elements is 0 . The first two authors generalized Forster's theorem to forms of other algebras (not necessarily commutative, associative or unital) A. Shukla and the third author then showed that this generalized Forster bound is optimal for 'etale algebras. In this paper, we prove new upper and lower bound on the number of generators of an R form of a k algebra, where k is an infinite field and R has finite transcendence degree d over k . In particular, we show that, contrary to expectations, for most types of algebras, the generalized Forster bound is far from optimal. Our results are particularly detailed in the case of Azumaya algebras. Our proofs are based on reinterpreting the problem as a question about approximating the classifying stack BG , where G is the automorphism group of the algebra in question, by algebraic spaces of a certain form.


An infinity operad of normalized cacti. We show that normalized cacti form an infty operad in the form of a dendroidal space satisfying a weak Segal condition. To do this, we introduce a new topological operad of bracketed trees and an enrichment of the dendroidal category Omega .


The Cry Wolf Effect in Evacuation a Game Theoretic Approach. In today's terrorism prone and security focused world, evacuation emergencies, drills, and false alarms are becoming more and more common. Compliance to an evacuation order made by an authority in case of emergency can play a key role in the outcome of an emergency. In case an evacuee experiences repeated emergency scenarios which may be a false alarm (e.g., an evacuation drill, a false bomb threat, etc.) or an actual threat, the Aesop's cry wolf effect (repeated false alarms decrease order compliance) can severely affect his her likelihood to evacuate. To analyse this key unsolved issue of evacuation research, a game theoretic approach is proposed. Game theory is used to explore mutual best responses of an evacuee and an authority. In the proposed model the authority obtains a signal of whether there is a threat or not and decides whether to order an evacuation or not. The evacuee, after receiving an evacuation order, subsequently decides whether to stay or leave based on posterior beliefs that have been updated in response to the authority's action. Best responses are derived and Sequential equilibrium and Perfect Bayesian Equilibrium are used as solution concepts (refining equilibria with the intuitive criterion). Model results highlight the benefits of announced evacuation drills and suggest that improving the accuracy of threat detection can prevent large inefficiencies associated with the cry wolf effect.