Rigid local systems and the multiplicative eigenvalue problem. We give a construction which produces irreducible complex rigid local systems on Bbb P Bbb C 1 p 1, dots,p s via quantum Schubert calculus and strange duality. These local systems are unitary and arise from a study of vertices in the polytopes controlling the multiplicative eigenvalue problem for the special unitary groups operatorname SU (n) (i.e., determination of the possible eigenvalues of a product of unitary matrices given the eigenvalues of the matrices). Roughly speaking, we show that the strange duals of the simplest vertices of these polytopes give all possible unitary irreducible rigid local systems. As a consequence we obtain that the ranks of unitary irreducible rigid local systems, including those with finite global monodromy, on Bbb P 1 S are bounded above if we fix the cardinality of the set S p 1, dots,p s and require that the local monodromies have orders which divide n , for a fixed n . Answering a question of N. Katz, we show that there are no irreducible rigid local systems of rank greater than one, with finite global monodromy, all of whose local monodromies have orders dividing n , when n is a prime number. We also show that all unitary irreducible rigid local systems on Bbb P 1 Bbb C S with finite local monodromies arise as solutions to the Knizhnik Zamalodchikov equations on conformal blocks for the special linear group. Along the way, generalising previous works of the author and J. Kiers, we give an inductive mechanism for determining all vertices in the multiplicative eigenvalue problem for operatorname SU (n) .


Hardy Littlewood Sobolev inequalities for a class of non symmetric and non doubling hypoelliptic semigroups. In his seminal 1934 paper on Brownian motion and the theory of gases Kolmogorov introduced a second order evolution equation which displays some challenging features. In the opening of his 1967 hypoellipticity paper H "ormander discussed a general class of degenerate Ornstein Uhlenbeck operators that includes Kolmogorov's as a special case. In this note we combine semigroup theory with a nonlocal calculus for these hypoelliptic operators to establish new inequalities of Hardy Littlewood Sobolev type in the situation when the drift matrix has nonnegative trace. Our work has been influenced by ideas of E. Stein and Varopoulos in the framework of symmetric semigroups. One of our objectives is to show that such ideas can be pushed to successfully handle the present degenerate non symmetric setting.


Representation Formula for Viscosity Solutions to Parabolic PDEs with Sublinear Operators. We provide a representation formula for viscosity solutions to a class of nonlinear second order parabolic PDE problem involving sublinear operators. This is done through a dynamic programming principle derived from 8 . The formula can be seen as a nonlinear extension of the Feynman Kac formula and is based on the backward stochastic differential equations theory.


A natural extension of Mittag Leffler function associated with a triple infinite series. We establish a new natural extension of Mittag Leffler function with three variables which is so called "trivariate Mittag Leffler function". The trivariate Mittag Leffler function can be expressed via complex integral representation by putting to use of the eminent Hankel's integral. We also investigate Laplace integral relation and convolution result for a univariate version of this function. Moreover, we present fractional derivative of trivariate Mittag Leffler function in Caputo type and we also discuss Riemann Liouville type fractional integral and derivative of this function. The link of trivariate Mittag Leffler function with fractional differential equation systems involving different fractional orders is necessary on certain applications in physics. Thus, we provide an exact analytic solutions of homogeneous and inhomogeneous multi term fractional differential equations by means of a newly defined trivariate Mittag Leffler functions.


Cubic vertex transitive graphs of girth six. In this paper, a complete classification of finite simple cubic vertex transitive graphs of girth 6 is obtained. It is proved that every such graph, with the exception of the Desargues graph on 20 vertices, is either a skeleton of a hexagonal tiling of the torus, the skeleton of the truncation of an arc transitive triangulation of a closed hyperbolic surface, or the truncation of a 6 regular graph with respect to an arc transitive dihedral scheme. Cubic vertex transitive graphs of girth larger than 6 are also discussed.


Generalized Planning Non Deterministic Abstractions and Trajectory Constraints. We study the characterization and computation of general policies for families of problems that share a structure characterized by a common reduction into a single abstract problem. Policies mu that solve the abstract problem P have been shown to solve all problems Q that reduce to P provided that mu terminates in Q. In this work, we shed light on why this termination condition is needed and how it can be removed. The key observation is that the abstract problem P captures the common structure among the concrete problems Q that is local (Markovian) but misses common structure that is global. We show how such global structure can be captured by means of trajectory constraints that in many cases can be expressed as LTL formulas, thus reducing generalized planning to LTL synthesis. Moreover, for a broad class of problems that involve integer variables that can be increased or decreased, trajectory constraints can be compiled away, reducing generalized planning to fully observable non deterministic planning.


Accept Synthetic Objects as Real End to End Training of Attentive Deep Visuomotor Policies for Manipulation in Clutter. Recent research demonstrated that it is feasible to end to end train multi task deep visuomotor policies for robotic manipulation using variations of learning from demonstration (LfD) and reinforcement learning (RL). In this paper, we extend the capabilities of end to end LfD architectures to object manipulation in clutter. We start by introducing a data augmentation procedure called Accept Synthetic Objects as Real (ASOR). Using ASOR we develop two network architectures implicit attention ASOR IA and explicit attention ASOR EA. Both architectures use the same training data (demonstrations in uncluttered environments) as previous approaches. Experimental results show that ASOR IA and ASOR EA succeed ina significant fraction of trials in cluttered environments where previous approaches never succeed. In addition, we find that both ASOR IA and ASOR EA outperform previous approaches even in uncluttered environments, with ASOR EA performing better even in clutter compared to the previous best baseline in an uncluttered environment.


Sense Deliberate Act Cognitive Agents for Sense Compute Control Applications in the Internet of Things Services. In this paper, we advocate Agent Oriented Software Engi neering (AOSE) through employing Belief Desire Intention (BDI) intel ligent agents for developing Sense Compute Control (SCC) applications in the Internet of Things and Services (IoTS). We argue that not only the agent paradigm, in general, but also cognitive BDI agents with sense deliberate act cycle, in particular, fit very well to the nature of SCC applications in the IoTS. However, considering the highly constrained heterogeneous devices that are prevalent in the IoTS, existing BDI agent frameworks, even those especially created for Wireless Sensor Networks (WSNs), do not work. We elaborate on the challenges and propose pos sible approaches to address them.


Direct observation of structural heterogeneity and tautomerization of single hypericin molecules. Tautomerization is a fast chemical reaction where structures of the reactants differ only in the position of a proton and a double bond. Tautomerization often occurs in natural substances and is a fundamental process in organic and biochemistry. However, studying the optical properties of tautomeric species is challenging due to ensemble averaging. Many molecules, such as porphines, porphycenes or phenanthroperylene quinones, exhibit a reorientation of the transition dipole moment (TDM) during tautomerization, which can be directly observed in a single molecule experiment. A prominent phenanthroperylene quinone is hypericin showing antiviral, antidepressive, and photodynamical properties. Here, we study single hypericin molecules by using confocal microscopy combined with higher order laser modes. Observing abrupt flipping of the image pattern allows to draw conclusions about the coexistence of different tautomers and their conversion path. Time dependent density functional theory calculations show that hypericin is cycling between the four most stable tautomers. This approach allows to unambiguously assign a TDM orientation to a specific tautomer and enables to determine the chemical structure in situ. Additionally, tautomerization can not only be observed by the image pattern orientation, but also as intermittency in the fluorescence emission of a single molecule. Time correlated single photon counting enables to determine the excited state lifetimes of the hypericin tautomers. Our approach is not only limited to hypericin, but can be applied to other molecules showing a TDM reorientation during tautomerization, helping to get a deeper understanding of this important process.


Chaotic wave dynamics in weakly magnetised spherical Couette flows. Direct numerical simulations of a liquid metal filling the gap between two concentric spheres are presented. The flow is governed by the interplay between the rotation of the inner sphere (measured by the Reynolds number Re) and a weak externally applied axial magnetic field (measured by the Hartmann number Ha). By varying the latter a rich variety of flow features, both in terms of spatial symmetry and temporal dependence, is obtained. Flows with two or three independent frequencies describing their time evolution are found as a result of Hopf bifurcations. They are stable on a sufficiently large interval of Hartmann numbers where regions of multistability of two, three and even four types of these different flows are detected. The temporal character of the solutions is analysed by means of an accurate frequency analysis and Poincar 'e sections. An unstable branch of flows undergoing a period doubling cascade and frequency locking of three frequency solutions is described as well.


Convexification for a 3D inverse scattering problem with the moving point source. For the first time, we develop in this paper the globally convergent convexification numerical method for a Coefficient Inverse Problem for the 3D Helmholtz equation for the case when the backscattering data are generated by a point source running along an interval of a straight line and the wavenumber is fixed. Thus, by varying the wavenumber, one can reconstruct the dielectric constant depending not only on spatial variables but the wavenumber (i.e. frequency) as well. Our approach relies on a new derivation of a boundary value problem for a system of coupled quasilinear elliptic partial differential equations. This is done via an application of a special truncated Fourier like method. First, we prove the Lipschitz stability estimate for this problem via a Carleman estimate. Next, using the Carleman Weight Function generated by that estimate, we construct a globally strictly convex cost functional and prove the global convergence to the exact solution of the gradient projection method. Finally, our theoretical finding is verified via several numerical tests with computationally simulated data. These tests demonstrate that we can accurately recover all three important components of targets of interest locations, shapes and dielectric constants. In particular, large target background contrasts in dielectric constants (up to 10 1) can be accurately calculated.


A Comparative Study of Asynchronous Many Tasking Runtimes Cilk, Charm , ParalleX and AM . We evaluate and compare four contemporary and emerging runtimes for high performance computing(HPC) applications Cilk, Charm , ParalleX and AM . We compare along three bases programming model, execution model and the implementation on an underlying machine model. The comparison study includes a survey of each runtime system's programming models, their corresponding execution models, their stated features, and performance and productivity goals. We first qualitatively compare these runtimes with programmability in mind. The differences in expressivity and programmability arising from their syntax and semantics are clearly enunciated through examples common to all runtimes. Then, the execution model of each runtime, which acts as a bridge between the programming model and the underlying machine model, is compared and contrasted to that of the others. We also evaluate four mature implementations of these runtimes, namely Intel Cilk , Charm 6.5.1, AM and HPX, that embody the principles dictated by these models. With the emergence of the next generation of supercomputers, it is imperative for parallel programming models to evolve and address the integral challenges introduced by the increasing scale. Rather than picking a winner out of these four models under consideration, we end with a discussion on lessons learned, and how such a study is instructive in the evolution of parallel programming frameworks to address the said challenges.


Generalized Zero Shot Learning Via Over Complete Distribution. A well trained and generalized deep neural network (DNN) should be robust to both seen and unseen classes. However, the performance of most of the existing supervised DNN algorithms degrade for classes which are unseen in the training set. To learn a discriminative classifier which yields good performance in Zero Shot Learning (ZSL) settings, we propose to generate an Over Complete Distribution (OCD) using Conditional Variational Autoencoder (CVAE) of both seen and unseen classes. In order to enforce the separability between classes and reduce the class scatter, we propose the use of Online Batch Triplet Loss (OBTL) and Center Loss (CL) on the generated OCD. The effectiveness of the framework is evaluated using both Zero Shot Learning and Generalized Zero Shot Learning protocols on three publicly available benchmark databases, SUN, CUB and AWA2. The results show that generating over complete distributions and enforcing the classifier to learn a transform function from overlapping to non overlapping distributions can improve the performance on both seen and unseen classes.


Results from the Alfv ' e n Eigenmode Active Diagnostic during the 2019 2020 JET deuterium campaign. This paper presents results of extensive analysis of mode excitation observed during the operation of the Alfv ' e n Eigenmode Active Diagnostic (AEAD) in the JET tokamak during the 2019 2020 deuterium campaign. Six of eight toroidally spaced antennas, each with independent power and phasing, were successful in actively exciting stable MHD modes in 479 plasmas. In total, 4768 magnetic resonances were detected with up to fourteen fast magnetic probes. In this work, we present the calculations of resonant frequencies f 0 , damping rates gamma 0 , and toroidal mode numbers n , spanning the parameter range f 0 approx 30 250 kHz, gamma approx 0 13 kHz, and vert n vert leq 30 . In general, good agreement is seen between the resonant and the calculated toroidal Alfv ' e n Eigenmode frequencies, and between the toroidal mode numbers applied by the AEAD and estimated of the excited resonances. We note several trends in the database the probability of resonance detection decreases with plasma current and external heating power the normalized damping rate increases with edge safety factor but decreases with external heating. These results provide key information to prepare future experimental campaigns and to better understand the physics of excitation and damping of Alfv ' e n Eigenmodes in the presence of alpha particles during the upcoming DT campaign, thereby extrapolating with confidence to future tokamaks.


On the integration of transitive Lie algebroids. We revisit the problem of integrating Lie algebroids A Rightarrow M to Lie groupoids G rightrightarrows M , for the special case that the Lie algebroid A is transitive. We obtain a geometric explanation of the Crainic Fernandes obstructions for this situation, and an explicit construction of the integration whenever these obstructions vanish. We also indicate an extension of this approach to regular Lie algebroids.


Supervised Learning Algorithms for Controlling Underactuated Dynamical Systems. Control of underactuated dynamical systems has been studied for decades in robotics, and is now emerging in other fields such as neuroscience. Most of the advances have been in model based control theory, which has limitations when the system under study is very complex and it is not possible to construct a model. This calls for data driven control methods like machine learning, which has spread to many fields in the recent years including control theory. However, the success of such algorithms has been dependent on availability of large datasets. Moreover, due to their black box nature, it is challenging to analyze how such algorithms work, which may be crucial in applications where failure is very costly. In this paper, we develop two related novel supervised learning algorithms. The algorithms are powerful enough to control a wide variety of complex underactuated dynamical systems, and yet have a simple and intelligent structure that allows them to work with a sparse data set even in the presence of noise. Our algorithms output a bang bang (binary) control input by taking in feedback of the state of the dynamical system. The algorithms learn this control input by maximizing a reward function in both short and long time horizons. We demonstrate the versatility of our algorithms by applying them to a diverse range of applications including switching between bistable states, changing the phase of an oscillator, desynchronizing a population of synchronized coupled oscillators, and stabilizing an unstable fixed point. For most of these applications we are able to reason why our algorithms work by using traditional dynamical systems and control theory. We also compare our learning algorithms with some traditional control algorithms, and reason why our algorithms work better.


Robust Output Tracking for a Room Temperature Model with Distributed Control and Observation. We consider robust output regulation of a partial differential equation model describing temperature evolution in a room. More precisely, we examine a two dimensional room model with the velocity field and temperature evolution governed by the incompressible steady state Navier Stokes and advection diffusion equations, respectively, which coupled together form a simplification of the Boussinesq equations. We assume that the control and observation operators of our system are distributed, whereas the disturbance acts on a part of the boundary of the system. We solve the robust output regulation problem using a finite dimensional low order controller, which is constructed using model reduction on a finite element approximation of the model. Through numerical simulations, we compare performance of the reduced order controller to that of the controller without model reduction as well as to performance of a low gain robust controller.


Analysis and a Solution of Momentarily Missed Detection for Anchor based Object Detectors. The employment of convolutional neural networks has led to significant performance improvement on the task of object detection. However, when applying existing detectors to continuous frames in a video, we often encounter momentary miss detection of objects, that is, objects are undetected exceptionally at a few frames, although they are correctly detected at all other frames. In this paper, we analyze the mechanism of how such miss detection occurs. For the most popular class of detectors that are based on anchor boxes, we show the followings i) besides apparent causes such as motion blur, occlusions, background clutters, etc., the majority of remaining miss detection can be explained by an improper behavior of the detectors at boundaries of the anchor boxes and ii) this can be rectified by improving the way of choosing positive samples from candidate anchor boxes when training the detectors.


On the Estimation of Mutual Information. In this paper we focus on the estimation of mutual information from finite samples ( mathcal X times mathcal Y ) . The main concern with estimations of mutual information is their robustness under the class of transformations for which it remains invariant i.e. type I (coordinate transformations), III (marginalizations) and special cases of type IV (embeddings, products). Estimators which fail to meet these standards are not textit robust in their general applicability. Since most machine learning tasks employ transformations which belong to the classes referenced in part I, the mutual information can tell us which transformations are most optimal cite Carrara Ernst . There are several classes of estimation methods in the literature, such as non parametric estimators like the one developed by Kraskov et. al cite KSG , and its improved versions cite LNC . These estimators are extremely useful, since they rely only on the geometry of the underlying sample, and circumvent estimating the probability distribution itself. We explore the robustness of this family of estimators in the context of our design criteria.


Excited state diabatization on the cheap using DFT Photoinduced electron and hole transfer. Excited state electron and hole transfer underpin fundamental steps in processes such as exciton dissociation at photovoltaic heterojunctions, photoinduced charge transfer at electrodes, and electron transfer in photosynthetic reaction centers. Diabatic states corresponding to charge or excitation localized species, such as locally excited and charge transfer states, provide a physically intuitive framework to simulate and understand these processes. However, obtaining accurate diabatic states and their couplings from adiabatic electronic states generally leads to inaccurate results when combined with low tier electronic structure methods, such as time dependent density functional theory (TDDFT), and exorbitant computational cost when combined with high level wavefunction based methods. Here we introduce a DFT based diabatization scheme, Delta ALMO(MSDFT2), which directly constructs the diabatic states using absolutely localized molecular orbitals (ALMOs). We demonstrate that our method, which combines ALMO calculations with the Delta SCF technique to construct electronically excited diabatic states and obtains their couplings with charge transfer states using our MSDFT2 scheme, gives accurate results for excited state electron and hole transfer in both charged and uncharged systems that underlie DNA repair, charge separation in donor acceptor dyads, chromophore to solvent electron transfer, and singlet fission. This framework for the accurate and efficient construction of excited state diabats and evaluation of their couplings directly from DFT thus offers a route to simulate and elucidate photoinduced electron and hole transfer in large disordered systems, such as those encountered in the condensed phase.


Neural Re Simulation for Generating Bounces in Single Images. We introduce a method to generate videos of dynamic virtual objects plausibly interacting via collisions with a still image's environment. Given a starting trajectory, physically simulated with the estimated geometry of a single, static input image, we learn to 'correct' this trajectory to a visually plausible one via a neural network. The neural network can then be seen as learning to 'correct' traditional simulation output, generated with incomplete and imprecise world information, to obtain context specific, visually plausible re simulated output, a process we call neural re simulation. We train our system on a set of 50k synthetic scenes where a virtual moving object (ball) has been physically simulated. We demonstrate our approach on both our synthetic dataset and a collection of real life images depicting everyday scenes, obtaining consistent improvement over baseline alternatives throughout.


Polytopic Discontinuous Galerkin methods for the numerical modelling of flow in porous media with networks of intersecting fractures. We present a numerical approximation of Darcy's flow through a porous medium that incorporates networks of fractures with non empty intersection. Our scheme employs PolyDG methods, i.e. discontinuous Galerkin methods on general polygonal and polyhedral (polytopic, for short) grids, featuring elements with edges faces that may be in arbitrary number (potentially unlimited) and whose measure may be arbitrarily small. Our approach is then very well suited to tame the geometrical complexity featured by most of applications in the computational geoscience field. From the modelling point of view, we adopt a reduction strategy that treats fractures as manifolds of codimension one and we employ the primal version of Darcy's law to describe the flow in both the bulk and the fracture network. In addition, some physically consistent conditions couple the two problems, allowing for jump of pressure at their interface, and they as well prescribe the behaviour of the fluid along the intersections, imposing pressure continuity and flux conservation. Both the bulk and fracture discretizations are obtained employing the Symmetric Interior Penalty DG method extended to the polytopic setting. The key instrument to obtain a polyDG approximation of the problem in the fracture network is the generalization of the concepts of jump and average at the intersection, so that the contribution from all the fractures is taken into account. We prove the well posedness of the discrete formulation and perform an error analysis obtaining a priori hp error estimates. All our theoretical results are validated performing preliminary numerical tests with known analytical solution.


Manifest Contracts with Intersection Types. We present a manifest contract system PCFv Delta H with intersection types. A manifest contract system is a typed functional calculus in which software contracts are integrated into a refinement type system and consistency of contracts is checked by combination of compile and run time type checking. Intersection types naturally arise when a contract is expressed by a conjunction of smaller contracts. Run time contract checking for conjunctive higher order contracts in an untyped language has been studied but our typed setting poses an additional challenge due to the fact that an expression of an intersection type tau 1 wedge tau 2 may have to perform different run time checking whether it is used as tau 1 or tau 2 . We build PCFv Delta H on top of the Delta calculus, a Church style intersection type system by Liquori and Stolze. In the Delta calculus, a canonical expression of an intersection type is a strong pair, whose elements are the same expressions except for type annotations. To address the challenge above, we relax strong pairs so that expressions in a pair are the same except for type annotations and casts, which are a construct for run time checking. We give a formal definition of PCFv Delta H and show its basic properties as a manifest contract system preservation, progress, and value inversion. Furthermore, we show that run time checking does not affect essential computation.


A Bayesian multiscale CNN framework to predict local stress fields in structures with microscale features. Multiscale computational modelling is challenging due to the high computational cost of direct numerical simulation by finite elements. To address this issue, concurrent multiscale methods use the solution of cheaper macroscale surrogates as boundary conditions to microscale sliding windows. The microscale problems remain a numerically challenging operation both in terms of implementation and cost. In this work we propose to replace the local microscale solution by an Encoder Decoder Convolutional Neural Network that will generate fine scale stress corrections to coarse predictions around unresolved microscale features, without prior parametrisation of local microscale problems. We deploy a Bayesian approach providing credible intervals to evaluate the uncertainty of the predictions, which is then used to investigate the merits of a selective learning framework. We will demonstrate the capability of the approach to predict equivalent stress fields in porous structures using linearised and finite strain elasticity theories.


On the equivalence of the Hermitian eigenvalue problem and hypergraph edge elimination. It is customary to identify sparse matrices with the corresponding adjacency or incidence graph. For the solution of linear systems of equations using Gaussian elimination, the representation by its adjacency graph allows a symbolic computation that can be used to predict memory footprints and enables the determination of near optimal elimination orderings based on heuristics. The Hermitian eigenvalue problem on the other hand seems to evade such treatment at first glance due to its inherent iterative nature. In this paper we prove this assertion wrong by showing the equivalence of the Hermitian eigenvalue problem with a symbolic edge elimination procedure. A symbolic calculation based on the incidence graph of the matrix can be used in analogy to the symbolic phase of Gaussian elimination to develop heuristics which reduce memory footprint and computations. Yet, we also show that the question of an optimal elimination strategy remains NP hard, in analogy to the linear systems case.


Generalized Framework for Nonlinear Acceleration. Nonlinear acceleration algorithms improve the performance of iterative methods, such as gradient descent, using the information contained in past iterates. However, their efficiency is still not entirely understood even in the quadratic case. In this paper, we clarify the convergence analysis by giving general properties that share several classes of nonlinear acceleration Anderson acceleration (and variants), quasi Newton methods (such as Broyden Type I or Type II, SR1, DFP, and BFGS) and Krylov methods (Conjugate Gradient, MINRES, GMRES). In particular, we propose a generic family of algorithms that contains all the previous methods and prove its optimal rate of convergence when minimizing quadratic functions. We also propose multi secants updates for the quasi Newton methods listed above. We provide a Matlab code implementing the algorithm.


Multi modal gated recurrent units for image description. Using a natural language sentence to describe the content of an image is a challenging but very important task. It is challenging because a description must not only capture objects contained in the image and the relationships among them, but also be relevant and grammatically correct. In this paper a multi modal embedding model based on gated recurrent units (GRU) which can generate variable length description for a given image. In the training step, we apply the convolutional neural network (CNN) to extract the image feature. Then the feature is imported into the multi modal GRU as well as the corresponding sentence representations. The multi modal GRU learns the inter modal relations between image and sentence. And in the testing step, when an image is imported to our multi modal GRU model, a sentence which describes the image content is generated. The experimental results demonstrate that our multi modal GRU model obtains the state of the art performance on Flickr8K, Flickr30K and MS COCO datasets.


Robustness Issues of Timing and Synchronization for Free Electron Lasers. Free electron lasers (FEL) require strict time relations for the electron bunch and RF field interaction, which must be precise and deterministic in time. This is guaranteed by the timing and synchronization systems that should be robust under the situations like a power cycle in the master oscillator, timing master, reference frequency distribution devices or low level radio frequency (LLRF) devices. After the power cycles, the time relations should be kept or be capable to recover quickly to improve the availability of the FEL machine. This article focuses on the robustness of the timing and synchronization systems, such as the time uncertainty of the RF pulse related with the trigger, the phase uncertainty of frequency dividers after power cycles and the race condition between the trigger and clock. The possible solutions to achieve a robust design of the timing and synchronization systems are also discussed.


Flows Over Periodic Hills of Parameterized Geometries A Dataset for Data Driven Turbulence Modeling From Direct Simulations. Computational fluid dynamics models based on Reynolds averaged Navier Stokes equations with turbulence closures still play important roles in engineering design and analysis. However, the development of turbulence models has been stagnant for decades. With recent advances in machine learning, data driven turbulence models have become attractive alternatives worth further explorations. However, a major obstacle in the development of data driven turbulence models is the lack of training data. In this work, we survey currently available public turbulent flow databases and conclude that they are inadequate for developing and validating data driven models. Rather, we need more benchmark data from systematically and continuously varied flow conditions (e.g., Reynolds number and geometry) with maximum coverage in the parameter space for this purpose. To this end, we perform direct numerical simulations of flows over periodic hills with varying slopes, resulting in a family of flows over periodic hills which ranges from incipient to mild and massive separations. We further demonstrate the use of such a dataset by training a machine learning model that predicts Reynolds stress anisotropy based on a set of mean flow features. We expect the generated dataset, along with its design methodology and the example application presented herein, will facilitate development and comparison of future data driven turbulence models.


Microtubule based actin transport and localization in a spherical cell. The interaction between actin filaments and microtubules is crucial for many eukaryotic cellular processes, such as, among others, cell polarization, cell motility and cellular wound healing. The importance of this interaction has long been recognised, yet very little is understood about both the underlying mechanisms and the consequences for the spatial (re)organization of the cellular cytoskeleton. At the same time, understanding the causes and the consequences of the interaction between different biomolecular components are key questions for emph in vitro research involving reconstituted biomolecular systems, especially in the light of current interest in creating minimal synthetic cells. In this light, recent emph in vitro experiments have shown that the actin microtubule interaction mediated by the cytolinker TipAct, which binds to actin lattice and microtubule tip, causes the directed transport of actin filaments. We develop an analytical theory of dynamically unstable microtubules, nucleated from the center of a spherical cell, in interaction with actin filaments. We show that, depending on the balance between the diffusion of unbound actin filaments and propensity to bind microtubules, actin is either concentrated in the center of the cell, where the density of microtubules is highest, or becomes localized to the cell cortex.


Importance of suppression and mitigation measures in managing COVID 19 outbreaks. I employ a simple mathematical model of an epidemic process to evaluate how four basic quantities the reproduction number (R), the numbers of sensitive (S) and infectious individuals(I), and total community size (N) affect strategies to control COVID 19. Numerical simulations show that strict suppression measures at the beginning of an epidemic can create low infectious numbers, which thereafter can be managed by mitigation measures over longer periods to flatten the epidemic curve. The stronger the suppression measure, the faster it achieves the low numbers of infections that are conducive to subsequent management. We discuss the predictions of this analysis and how it fits into longer term sequences of measures, including using the herd immunity concept to leverage acquired immunity.


A stochastic kinetic scheme for multi scale plasma transport with uncertainty quantification. In this paper, a physics oriented stochastic kinetic scheme will be developed that includes random inputs from both flow and electromagnetic fields via a hybridization of stochastic Galerkin and collocation methods. Based on the BGK type relaxation model of the multi component Boltzmann equation, a scale dependent kinetic central upwind flux function is designed in both physical and particle velocity space, and the governing equations in the discrete temporal spatial random domain are constructed. By solving Maxwell's equations with the wave propagation method, the evolutions of ions, electrons and electromagnetic field are coupled throughout the simulation. We prove that the scheme is formally asymptotic preserving in the Vlasov, magnetohydrodynamical, and neutral Euler regimes with the inclusion of random variables. Therefore, it can be used for the study of multi scale and multi physics plasma system under the effects of uncertainties, and provide scale adaptive physical solutions under different ratios among numerical cell size, particle mean free path and gyroradius (or time step, local particle collision time and plasma period). Numerical experiments including one dimensional Landau Damping, the two stream instability and the Brio Wu shock tube problem with one to three dimensional velocity settings, and each under stochastic initial conditions with one dimensional uncertainty, will be presented to validate the scheme.


CNN based Lidar Point Cloud De Noising in Adverse Weather. Lidar sensors are frequently used in environment perception for autonomous vehicles and mobile robotics to complement camera, radar, and ultrasonic sensors. Adverse weather conditions are significantly impacting the performance of lidar based scene understanding by causing undesired measurement points that in turn effect missing detections and false positives. In heavy rain or dense fog, water drops could be misinterpreted as objects in front of the vehicle which brings a mobile robot to a full stop. In this paper, we present the first CNN based approach to understand and filter out such adverse weather effects in point cloud data. Using a large data set obtained in controlled weather environments, we demonstrate a significant performance improvement of our method over state of the art involving geometric filtering. Data is available at https github.com rheinzler PointCloudDeNoising.


The vimentin cytoskeleton When polymer physics meets cell biology. The proper functions of tissues depend on the ability of cells to withstand stress and maintain shape. Central to this process is the cytoskeleton, comprised of three polymeric networks F actin, microtubules, and intermediate filaments. Intermediate filament proteins are among the most abundant cytoskeletal proteins in cells yet they remain one of the least understood. Their structure and function deviate from those of their cytoskeletal partners, F actin and microtubules. Intermediate filament networks show a unique combination of extensibility, flexibility and toughness that confers mechanical resilience to the cell. Vimentin is an intermediate filament protein expressed in mesenchymal cells. This review highlights exciting new results on the physical biology of vimentin intermediate filaments and their role in allowing whole cells and tissues to cope with stress.


Multi frequency coherent emission from superstructure thermal emitters. Long range spatial coherence can be induced in thermal emitters by embedding a periodic grating into a material supporting propagating polaritons or dielectric modes. However, the emission angle and frequency cannot be defined simultaneously and uniquely, resulting in emission at unusable angles or frequencies. Here, we explore superstructure gratings (SSGs) to control the spatial and spectral properties of thermal emitters. SSGs have long range periodicity, but a unit cell that provides tailorable Bragg components to interact with light. These Bragg components allow simultaneous launching of polaritons with different frequencies wavevectors in a single grating, manifesting as additional spatial and spectral bands upon the emission profile. As the unit cell period approaches the spatial coherence length, the coherence properties of the superstructure will be lost. Whilst the 1D k space representation of the grating provides insights into the emission, the etch depth of the grating can result in strong polariton polariton interactions. An emergent effect of these interactions is the creation of polaritonic band gaps, and defect states that can have a well defined frequency and emission angle. In all, our results show experimentally how even in simple 1D gratings there is significant design flexibility for engineering the profile of thermal emitters, bound by finite coherence length.


High efficiency resonant metasurface germanium photodetector with ultra thin intrinsic layer. Photodetectors at telecommunications band with high efficiency and high speed are becoming increasingly important as the booming of big data, 5G, internet of things, cloud computing, artificial intelligence and relevent applications. Silicon based Germanium photodetectors exhibit great potential in reducing the cost and power dissipation, due to its compatibility of monolithic integration with signal processing electronics. We report the first demonstration of normal incident resonant metasurface germanium photodetector, to address the trade off between quantum efficiency, bandwidth and wavelength coverage for free space detectors. With an ultra thin intrinsic layer thickness of 350 nm, a high external quantum efficiency of more than 60 and clearly open eyes at the speed of 20 Gbps are achieved, for a 30 mu m diameter device. The photodetector employs multiple trapped mode resonances to enhance the localized electromagnetic field, which not only enhances the external quantum efficiency by more than 300 at 1550 nm, but also extends to the whole C band. Further simulation and measurement with small incident spot size show the feasibility to achieve 50 GHz 3 dB bandwidth by simply reducing the mesa size, with minor sacrifice of the enhanced absorption. This work paves the way for the future development of low cost, high efficiency normal incident germanium photodetectors operating at data rates of 50 Gbps or higher.


The existence of a nontrivial weak solution to a double critical problem involving fractional Laplacian in R n with a Hardy term. In this paper, we consider the existence of nontrivial weak solutions to a double critical problem involving fractional Laplacian with a Hardy term begin equation label eq0.1 ( Delta) s u gamma frac u x 2s frac u 2 s ( beta) 2 u x beta big I mu F alpha ( cdot,u) big (x)f alpha (x,u), u in dot H s( R n ) end equation where s in(0,1) , 0 leq alpha, beta 2s n , mu in (0,n) , gamma gamma H , I mu (x) x mu , F alpha (x,u) frac u(x) 2 mu ( alpha) x delta mu ( alpha) , f alpha (x,u) frac u(x) 2 mu ( alpha) 2 u(x) x delta mu ( alpha) , 2 mu ( alpha) (1 frac mu 2n ) cdot 2 s ( alpha) , delta mu ( alpha) (1 frac mu 2n ) alpha , 2 s ( alpha) frac 2(n alpha) n 2s and gamma H 4 s frac Gamma 2( frac n 2s 4 ) Gamma 2( frac n 2s 4 ) . We show that problem ( ref eq0.1 ) admits at least a weak solution under some conditions. To prove the main result, we develop some useful tools based on a weighted Morrey space. To be precise, we discover the embeddings begin equation label eq0.2 dot H s( R n ) hookrightarrow L 2 s ( alpha) ( R n , y alpha ) hookrightarrow L p, frac n 2s 2 p pr ( R n , y pr ) end equation where s in (0,1) , 0 alpha 2s n , p in 1,2 s ( alpha)) , r frac alpha 2 s ( alpha) We also establish an improved Sobolev inequality. By using mountain pass lemma along with an improved Sobolev inequality, we obtain a nontrivial weak solution to problem ( ref eq0.1 ) in a direct way. It is worth while to point out that the improved Sobolev inequality could be applied to simplify the proof of the main results in cite NGSS and cite RFPP .


Sylvester type quaternion matrix equations with arbitrary equations and arbitrary unknowns. In this paper, we prove a conjecture which was presented in a recent paper Linear Algebra Appl. 2016 496 549 593 . We derive some practical necessary and sufficient conditions for the existence of a solution to a system of coupled two sided Sylvester type quaternion matrix equations with arbitrary equations and arbitrary unknowns A i X i B i C i X i 1 D i E i , i overline 1,k . As an application, we give some practical necessary and sufficient conditions for the existence of an eta Hermitian solution to the system of quaternion matrix equations A i X i A eta i C i X i 1 C eta i E i in terms of ranks, i overline 1,k .


Applications of entropy of product systems higher rank graphs. We consider C algebras of finite higher rank graphs along with their rotational action. We show how the entropy theory of product systems with finite frames applies to identify the phase transitions of the dynamics. We compute the positive inverse temperatures where symmetry breaks, and in particular we identify the subharmonic parts of the gauge invariant equilibrium states. Our analysis applies to positively weighted rotational actions through a recalibration of the entropies.


A computer simulation protocol to assess the accuracy of a Radio Stereometric Analysis (RSA) image processor according to the ISO 5725. Radio Stereometric Analysis and x ray fluoroscopy are radiological techniques that require dedicated software to process data. The accurate calibration of these software is therefore critical. The aim of this work is to produce a protocol for evaluating the softwares' accuracy according to the ISO 5725. A series of computer simulations of the radiological setup and images were employed. The noise level of the images was also changed to evaluate the accuracy with different image qualities. The protocol was tested on a custom software developed by the authors. Radiological scene reconstruction accuracy was of (0.092 0.14) mm for tube position, and (0.38 0.31) mm (2.09 1.39) deg for detectors oriented in a direction other than the source detector direction. In the source detector direction the accuracy was of (2.68 3.08) mm for tube position, and of (0.16 0.27) mm (0.075 1.16) deg for the detectors. These disparate results are widely discussed in the literature. Model positioning and orientation was also highly accurate (0.22 0.46) mm (0.26 0.22) deg. Accuracy was not affected by the noise level. The protocol was able to assess the accuracy of the RSA system. It was also useful to detect and fix hidden bugs. It was also useful to detect and resolve hidden bugs in the software, and in optimizing the algorithms.


Unsupervised Intrusion Detection System for Unmanned Aerial Vehicle with Less Labeling Effort. Along with the importance of safety, an IDS has become a significant task in the real world. Prior studies proposed various intrusion detection models for the UAV. Past rule based approaches provided a concrete baseline IDS model, and the machine learning based method achieved a precise intrusion detection performance on the UAV with supervised learning models. However, previous methods have room for improvement to be implemented in the real world. Prior methods required a large labeling effort on the dataset, and the model could not identify attacks that were not trained before. To jump over these hurdles, we propose an IDS with unsupervised learning. As unsupervised learning does not require labeling, our model let the practitioner not to label every type of attack from the flight data. Moreover, the model can identify an abnormal status of the UAV regardless of the type of attack. We trained an autoencoder with the benign flight data only and checked the model provides a different reconstruction loss at the benign flight and the flight under attack. We discovered that the model produces much higher reconstruction loss with the flight under attack than the benign flight thus, this reconstruction loss can be utilized to recognize an intrusion to the UAV. With consideration of the computation overhead and the detection performance in the wild, we expect our model can be a concrete and practical baseline IDS on the UAV.


Delayed Feedback Modeling for the Entire Space Conversion Rate Prediction. Estimating post click conversion rate (CVR) accurately is crucial in E commerce. However, CVR prediction usually suffers from three major challenges in practice i) data sparsity compared with impressions, conversion samples are often extremely scarce ii) sample selection bias conventional CVR models are trained with clicked impressions while making inference on the entire space of all impressions iii) delayed feedback many conversions can only be observed after a relatively long and random delay since clicks happened, resulting in many false negative labels during training. Previous studies mainly focus on one or two issues while ignoring the others. In this paper, we propose a novel neural network framework ESDF to tackle the above three challenges simultaneously. Unlike existing methods, ESDF models the CVR prediction from a perspective of entire space, and combines the advantage of user sequential behavior pattern and the time delay factor. Specifically, ESDF utilizes sequential behavior of user actions on the entire space with all impressions to alleviate the sample selection bias problem. By sharing the embedding parameters between CTR and CVR networks, data sparsity problem is greatly relieved. Different from conventional delayed feedback methods, ESDF does not make any special assumption about the delay distribution. We discretize the delay time by day slot and model the probability based on survival analysis with deep neural network, which is more practical and suitable for industrial situations. Extensive experiments are conducted to evaluate the effectiveness of our method. To the best of our knowledge, ESDF is the first attempt to unitedly solve the above three challenges in CVR prediction area.


Temporal LiDAR Frame Prediction for Autonomous Driving. Anticipating the future in a dynamic scene is critical for many fields such as autonomous driving and robotics. In this paper we propose a class of novel neural network architectures to predict future LiDAR frames given previous ones. Since the ground truth in this application is simply the next frame in the sequence, we can train our models in a self supervised fashion. Our proposed architectures are based on FlowNet3D and Dynamic Graph CNN. We use Chamfer Distance (CD) and Earth Mover's Distance (EMD) as loss functions and evaluation metrics. We train and evaluate our models using the newly released nuScenes dataset, and characterize their performance and complexity with several baselines. Compared to directly using FlowNet3D, our proposed architectures achieve CD and EMD nearly an order of magnitude lower. In addition, we show that our predictions generate reasonable scene flow approximations without using any labelled supervision.


Pixel Level Cycle Association A New Perspective for Domain Adaptive Semantic Segmentation. Domain adaptive semantic segmentation aims to train a model performing satisfactory pixel level predictions on the target with only out of domain (source) annotations. The conventional solution to this task is to minimize the discrepancy between source and target to enable effective knowledge transfer. Previous domain discrepancy minimization methods are mainly based on the adversarial training. They tend to consider the domain discrepancy globally, which ignore the pixel wise relationships and are less discriminative. In this paper, we propose to build the pixel level cycle association between source and target pixel pairs and contrastively strengthen their connections to diminish the domain gap and make the features more discriminative. To the best of our knowledge, this is a new perspective for tackling such a challenging task. Experiment results on two representative domain adaptation benchmarks, i.e. GTAV rightarrow Cityscapes and SYNTHIA rightarrow Cityscapes, verify the effectiveness of our proposed method and demonstrate that our method performs favorably against previous state of the arts. Our method can be trained end to end in one stage and introduces no additional parameters, which is expected to serve as a general framework and help ease future research in domain adaptive semantic segmentation. Code is available at https github.com kgl prml Pixel Level Cycle Association.


Real Time Multi Diver Tracking and Re identification for Underwater Human Robot Collaboration. Autonomous underwater robots working with teams of human divers may need to distinguish between different divers, e.g. to recognize a lead diver or to follow a specific team member. This paper describes a technique that enables autonomous underwater robots to track divers in real time as well as to reidentify them. The approach is an extension of Simple Online Realtime Tracking (SORT) with an appearance metric (deep SORT). Initial diver detection is performed with a custom CNN designed for realtime diver detection, and appearance features are subsequently extracted for each detected diver. Next, realtime tracking by detection is performed with an extension of the deep SORT algorithm. We evaluate this technique on a series of videos of divers performing human robot collaborative tasks and show that our methods result in more divers being accurately identified during tracking. We also discuss the practical considerations of applying multi person tracking to on board autonomous robot operations, and we consider how failure cases can be addressed during on board tracking.


Sensitivity Analysis of Submodular Function Maximization. We study the recently introduced idea of worst case sensitivity for monotone submodular maximization with cardinality constraint k , which captures the degree to which the output argument changes on deletion of an element in the input. We find that for large classes of algorithms that non trivial sensitivity of o(k) is not possible, even with bounded curvature, and that these results also hold in the distributed framework. However, we also show that in the regime k Omega(n) that we can obtain O(1) sensitivity for sufficiently low curvature.


Machine Learning based English Sentiment Analysis. Sentiment analysis or opinion mining aims to determine attitudes, judgments and opinions of customers for a product or a service. This is a great system to help manufacturers or servicers know the satisfaction level of customers about their products or services. From that, they can have appropriate adjustments. We use a popular machine learning method, being Support Vector Machine, combine with the library in Waikato Environment for Knowledge Analysis (WEKA) to build Java web program which analyzes the sentiment of English comments belongs one in four types of woman products. That are dresses, handbags, shoes and rings. We have developed and test our system with a training set having 300 comments and a test set having 400 comments. The experimental results of the system about precision, recall and F measures for positive comments are 89.3 , 95.0 and 92,.1 for negative comments are 97.1 , 78.5 and 86.8 and for neutral comments are 76.7 , 86.2 and 81.2 .


On the stochastic nonlinear Schr "odinger equations at critical regularities. We consider the Cauchy problem for the defocusing stochastic nonlinear Schr "odinger equations (SNLS) with an additive noise in the mass critical and energy critical settings. By adapting the probabilistic perturbation argument employed in the context of the random data Cauchy theory by the first author with B 'enyi and Pocovnicu (2015) to the current stochastic PDE setting, we present a concise argument to establish global well posedness of the mass critical and energy critical SNLS.


Probabilistic Predictions of People Perusing Evaluating Metrics of Language Model Performance for Psycholinguistic Modeling. By positing a relationship between naturalistic reading times and information theoretic surprisal, surprisal theory (Hale, 2001 Levy, 2008) provides a natural interface between language models and psycholinguistic models. This paper re evaluates a claim due to Goodkind and Bicknell (2018) that a language model's ability to model reading times is a linear function of its perplexity. By extending Goodkind and Bicknell's analysis to modern neural architectures, we show that the proposed relation does not always hold for Long Short Term Memory networks, Transformers, and pre trained models. We introduce an alternate measure of language modeling performance called predictability norm correlation based on Cloze probabilities measured from human subjects. Our new metric yields a more robust relationship between language model quality and psycholinguistic modeling performance that allows for comparison between models with different training configurations.


math.NA cs.NA physics.comp-ph physics.plasm-ph

A Multilevel Block Preconditioner for the HDG Trace System Applied to Incompressible Resistive MHD. We present a scalable block preconditioning strategy for the trace system coming from the high order hybridized discontinuous Galerkin (HDG) discretization of incompressible resistive magnetohydrodynamics (MHD). We construct the block preconditioner with a least squares commutator (BFBT) approximation for the inverse of the Schur complement that segregates out the pressure unknowns of the trace system. The remaining velocity, magnetic field, and Lagrange multiplier unknowns form a coupled nodal unknown block (the upper block), for which a system algebraic multigrid (AMG) is used for the approximate inverse. The complexity of the MHD equations together with the algebraic nature of the statically condensed HDG trace system makes the choice of smoother in the system AMG part critical for the convergence and performance of the block preconditioner. Our numerical experiments show GMRES preconditioned by ILU(0) of overlap zero as a smoother inside system AMG performs best in terms of robustness, time per nonlinear iteration and memory requirements. With several transient test cases in 2D and 3D including the island coalescence problem at high Lundquist number we demonstrate the robustness and parallel scalability of the block preconditioner. Additionally for the upper block a preliminary study of an alternate nodal block system solver based on a multilevel approximate nested dissection is presented. On a 2D island coalescence problem the multilevel approximate nested dissection preconditioner shows better scalability with respect to mesh refinement than the system AMG, but is relatively less robust with respect to Lundquist number scaling.


Homotopical decompositions of simplicial and Vietoris Rips complexes. Motivated by applications in Topological Data Analysis, we consider decompositions of a simplicial complex induced by a cover of its vertices. We study how the homotopy type of such decompositions approximates the homotopy of the simplicial complex itself. The difference between the simplicial complex and such an approximation is quantitatively measured by means of the so called obstruction complexes. Our general machinery is then specialized to clique complexes, Vietoris Rips complexes and Vietoris Rips complexes of metric gluings. For the latter we give metric conditions which allow to recover the first and zero th homology of the gluing from the respective homologies of the components.


State space construction of Hybrid Petri nets with multiple stochastic firings. Hybrid Petri nets have been extended to include general transitions that fire after a randomly distributed amount of time. With a single general one shot transition the state space and evolution over time can be represented either as a Parametric Location Tree or as a Stochastic Time Diagram. Recent work has shown that both representations can be combined and then allow multiple stochastic firings. This work presents an algorithm for building the Parametric Location Tree with multiple general transition firings and shows how its transient probability distribution can be computed using multi dimensional integration. We discuss the (dis )advantages of an interval arithmetic and a geometric approach to compute the areas of integration. Furthermore, we provide details on how to perform a Monte Carlo integration either directly on these intervals or convex polytopes, or after transformation to standard simplices. A case study on a battery backup system shows the feasibility of the approach and discusses the performance of the different integration approaches.


Part per billion measurement of the 4 2S 1 2 rightarrow 3 2D 5 2 electric quadrupole transition isotope shifts between 42,44,48 Ca and 40 Ca . We report a precise measurement of the isotope shifts in the 4 2 S 1 2 rightarrow 3 2 D 5 2 electric quadrupole transition at 729 nm in the 40 42,44,48 Ca . The measurement has been made via high resolution laser spectroscopy of co trapped ions, finding measured shifts of 2,771,872,467.6(7.6), 5,340,887,394.6(7.8), and 9,990,381,870.0(6.3) Hz between 42,44,48 Ca and 40 Ca , respectively. By exciting the two isotopes simultaneously using frequency sidebands derived from a single laser systematic uncertainties resulting from laser frequency drifts are eliminated. This permits far greater precision than similar previously published measurements in other alkaline earth systems. The resulting measurement precision provides a benchmark for tests of theoretical isotope shift calculations, and also offers a step towards probing New Physics via isotope shift spectroscopy.


Continuity of extensions of Lipschitz maps. We establish the sharp rate of continuity of extensions of mathbb R m valued 1 Lipschitz maps from a subset A of mathbb R n to a 1 Lipschitz maps on mathbb R n . We consider several cases when there exists a 1 Lipschitz extension with preserved uniform distance to a given 1 Lipschitz map. We prove that if m 1 then a given map is 1 Lipschitz and affine if and only if such distance preserving extension exists for any 1 Lipschitz map defined on any subset of mathbb R n . This shows a striking difference from the case m 1 , where any 1 Lipschitz function has such property. Another example where we prove it is possible to find an extension with the same Lipschitz constant and the same uniform distance to another Lipschitz map v is when the difference between the two maps takes values in a fixed one dimensional subspace of mathbb R m and the set A is geodesically convex with respect to a Riemannian pseudo metric associated with v .


High Harmonic Generation and Correlated Electron Emission from Relativistic Plasma Mirrors at 1 kHz Repetition Rate. We report evidence for the first generation of XUV spectra from relativistic surface high harmonic generation (SHHG) on plasma mirrors at a kilohertz repetition rate, emitted simultaneously with energetic electrons. SHHG spectra and electron angular distributions are measured as a function of the experimentally controlled plasma density gradient scale length L mathrm g for three increasingly short and intense driving pulses 24 fs and a 0 1.1 , 8 fs and a 0 1.6 , and finally 4 fs and a 0 approx2.1 . For all driver pulses, we observe correlated relativistic SHHG and electron emission in the range L mathrm g in lambda 25, lambda 6 , with an optimum gradient scale length of L mathrm g approx lambda 20 .


Even A cycles have the edge Erd H o s P 'osa property. I prove that even A cycles have the edge Erd H o s P 'osa property.


Twisted conjugacy and commensurability invariance. A group G is said to have property R infty if for every automorphism varphi in rm Aut (G) , the cardinality of the set of varphi twisted conjugacy classes is infinite. Many classes of groups are known to have such property. However, very few examples are known for which R infty is it geometric , i.e., if G has property R infty then any group quasi isometric to G also has property R infty . In this paper, we give examples of groups and conditions under which R infty is preserved under commensurability. The main tool is to employ the Bieri Neumann Strebel invariant.


Treewidth Pliability and PTAS for Max CSPs. We identify a sufficient condition, treewidth pliability, that gives a polynomial time approximation scheme (PTAS) for a large class of Max 2 CSPs parametrised by the class of allowed constraint graphs (with arbitrary constraints on an unbounded alphabet). Our result applies more generally to the maximum homomorphism problem between two rational valued structures. The condition unifies the two main approaches for designing PTASes. One is Baker's layering technique, which applies to sparse graphs such as planar or excluded minor graphs. The other is based on Szemer ' e di's regularity lemma and applies to dense graphs. We extend the applicability of both techniques to new classes of Max CSPs. Treewidth pliability turns out to be a robust notion that can be defined in several equivalent ways, including characterisations via size, treedepth, or the Hadwiger number. We show connections to the notions of fractional treewidth fragility from structural graph theory, hyperfiniteness from the area of property testing, and regularity partitions from the theory of dense graph limits. These may be of independent interest. In particular we show that a monotone class of graphs is hyperfinite if and only if it is fractionally treewidth fragile and has bounded degree.


SAWNet A Spatially Aware Deep Neural Network for 3D Point Cloud Processing. Deep neural networks have established themselves as the state of the art methodology in almost all computer vision tasks to date. But their application to processing data lying on non Euclidean domains is still a very active area of research. One such area is the analysis of point cloud data which poses a challenge due to its lack of order. Many recent techniques have been proposed, spearheaded by the PointNet architecture. These techniques use either global or local information from the point clouds to extract a latent representation for the points, which is then used for the task at hand (classification segmentation). In our work, we introduce a neural network layer that combines both global and local information to produce better embeddings of these points. We enhance our architecture with residual connections, to pass information between the layers, which also makes the network easier to train. We achieve state of the art results on the ModelNet40 dataset with our architecture, and our results are also highly competitive with the state of the art on the ShapeNet part segmentation dataset and the indoor scene segmentation dataset. We plan to open source our pre trained models on github to encourage the research community to test our networks on their data, or simply use them for benchmarking purposes.


Euclidean Bottleneck Bounded Degree Spanning Tree Ratios. Inspired by the seminal works of Khuller et al. (STOC 1994) and Chan (SoCG 2003) we study the bottleneck version of the Euclidean bounded degree spanning tree problem. A bottleneck spanning tree is a spanning tree whose largest edge length is minimum, and a bottleneck degree K spanning tree is a degree K spanning tree whose largest edge length is minimum. Let beta K be the supremum ratio of the largest edge length of the bottleneck degree K spanning tree to the largest edge length of the bottleneck spanning tree, over all finite point sets in the Euclidean plane. It is known that beta 5 1 , and it is easy to verify that beta 2 ge 2 , beta 3 ge sqrt 2 , and beta 4 1.175 . It is implied by the Hamiltonicity of the cube of the bottleneck spanning tree that beta 2 le 3 . The degree 3 spanning tree algorithm of Ravi et al. (STOC 1993) implies that beta 3 le 2 . Andersen and Ras (Networks, 68(4) 302 314, 2016) showed that beta 4 le sqrt 3 . We present the following improved bounds beta 2 ge sqrt 7 , beta 3 le sqrt 3 , and beta 4 le sqrt 2 . As a result, we obtain better approximation algorithms for Euclidean bottleneck degree 3 and degree 4 spanning trees. As parts of our proofs of these bounds we present some structural properties of the Euclidean minimum spanning tree which are of independent interest.


Self training for Few shot Transfer Across Extreme Task Differences. Most few shot learning techniques are pre trained on a large, labeled "base dataset". In problem domains where such large labeled datasets are not available for pre training (e.g., X ray, satellite images), one must resort to pre training in a different "source" problem domain (e.g., ImageNet), which can be very different from the desired target task. Traditional few shot and transfer learning techniques fail in the presence of such extreme differences between the source and target tasks. In this paper, we present a simple and effective solution to tackle this extreme domain gap self training a source domain representation on unlabeled data from the target domain. We show that this improves one shot performance on the target domain by 2.9 points on average on the challenging BSCD FSL benchmark consisting of datasets from multiple domains. Our code is available at https github.com cpphoo STARTUP.


Prescribing the curvature to Killing Graphs. In this work we prove the existence and uniqueness of Killing graphs with prescribed mean curvature considering functions which are not necessarily constant along the flow lines of the Killing vector field.


Three New Fundamental Problems in Physics of the Aftershocks. Recently, the physics of aftershocks has been enriched by three new problems. We will conditionally call them dynamic, inverse, and morphological problems. Dynamic problem is to search for the cumulative effect of a round the world seismic echo appearing after the main shock of earthquake. According to the theory, a converging surface seismic wave excited by the main shock returns to the epicenter about 3 hours after the main shock and stimulates the excitation of a strong aftershock. The second problem is an adequate description of the averaged evolution of the aftershock flow. We introduced a new concept about the coefficient of deactivation of the earthquake source, and proposed an equation that describes the evolution of aftershocks. Based on the equation of evolution, we posed and solved the inverse problem of the physics of the earthquake source. The third task is to model the spatio temporal distribution of aftershocks. Its solution refines our understanding of the structure of the earthquake source.


"The Perfect One" Understanding Communication Practices and Challenges with Animated GIFs. Animated GIFs are increasingly popular in text based communication. Finding the perfect GIF can make conversations funny, interesting, and engaging, but GIFs also introduce potentials for miscommunication. Through 24 in depth qualitative interviews, this empirical, exploratory study examines the nuances of communication practices with animated GIFs to better understand why and how GIFs can send unintentional messages. We find participants leverage contexts like source material and interpersonal relationship to find the perfect GIFs for different communication scenarios, while these contexts are also the primary reason for miscommunication and some technical usability issues in GIFs. This paper concludes with a discussion of the important role that different types of context play in the use and interpretations of GIFs, and argues that nonverbal communication tools should account for complex contexts and common ground that communication media rely on.


Clamp cell with in situ pressure monitoring for low temperature neutron scattering measurements. A clamp pressure cell for neutron scattering experiments at low temperatures and in external magnetic fields under pressure up to 2 GPa has been fabricated and tested. The cell provides optical access to the sample space that allows instantaneous pressure determination during sample loading, cooling and measuring using ruby and or samarium doped strontium tetraborate fluorescence monitoring. A new calibration curve of the pressure induced shift of the 7D 0 5F 0 (0 0) line in the fluorescent spectrum of SrB 4 O 7 Sm 2 for moderate pressures, P leqslant 2 GPa, is given.


On the Potential of Smarter Multi layer Maps. The most common way for robots to handle environmental information is by using maps. At present, each kind of data is hosted on a separate map, which complicates planning because a robot attempting to perform a task needs to access and process information from many different maps. Also, most often correlation among the information contained in maps obtained from different sources is not evaluated or exploited. In this paper, we argue that in robotics a shift from single source maps to a multi layer mapping formalism has the potential to revolutionize the way robots interact with knowledge about their environment. This observation stems from the raise in metric semantic mapping research, but expands to include in its formulation also layers containing other information sources, e.g., people flow, room semantic, or environment topology. Such multi layer maps, here named hypermaps, not only can ease processing spatial data information but they can bring added benefits arising from the interaction between maps. We imagine that a new research direction grounded in such multi layer mapping formalism for robots can use artificial intelligence to process the information it stores to present to the robot task specific information simplifying planning and bringing us one step closer to high level reasoning in robots.


A Distributed Black Box' Audit Trail Design Specification for Connected and Automated Vehicle Data and Software Assurance. Automotive software is increasingly complex and critical to safe vehicle operation, and related embedded systems must remain up to date to ensure long term system performance. Update mechanisms and data modification tools introduce opportunities for malicious actors to compromise these cyber physical systems, and for trusted actors to mistakenly install incompatible software versions. A distributed and stratified "black box" audit trail for automotive software and data provenance is proposed to assure users, service providers, and original equipment manufacturers (OEMs) of vehicular software integrity and reliability. The proposed black box architecture is both layered and diffuse, employing distributed hash tables (DHT), a parity system and a public blockchain to provide high resilience, assurance, scalability, and efficiency for automotive and other high assurance systems.


Voxel2Mesh 3D Mesh Model Generation from Volumetric Data. CNN based volumetric methods that label individual voxels now dominate the field of biomedical segmentation. However, 3D surface representations are often required for proper analysis. They can be obtained by post processing the labeled volumes which typically introduces artifacts and prevents end to end training. In this paper, we therefore introduce a novel architecture that goes directly from 3D image volumes to 3D surfaces without post processing and with better accuracy than current methods. We evaluate it on Electron Microscopy and MRI brain images as well as CT liver scans. We will show that it outperforms state of the art segmentation methods.


The Pohozaev Schoen Identity on Asymptotically Euclidean Manifolds Conservation Identities and their applications. The aim of this paper is to present a version of the generalized Pohozaev Schoen identity in the context of asymptotically euclidean manifolds. Since these kind of geometric identities have proven to be a very powerful tool when analysing different geometric problems for compact manifolds, we will present a variety of applications within this new context. Among these applications, we will show some rigidity results for asymptotically euclidean Ricci solitons and Codazzi solitons. Also, we will present an almost Schur type inequality valid in this non compact setting which does not need restrictions on the Ricci curvature. Finally, we will show how some rigidity results related with static potentials also follow from these type of conservation principles.


Type Centric Kotlin Compiler Fuzzing Preserving Test Program Correctness by Preserving Types. Kotlin is a relatively new programming language from JetBrains its development started in 2010 with release 1.0 done in early 2016. The Kotlin compiler, while slowly and steadily becoming more and more mature, still crashes from time to time on the more tricky input programs, not least because of the complexity of its features and their interactions. This makes it a great target for fuzzing, even the basic forms of which can find a significant number of Kotlin compiler crashes. There is a problem with fuzzing, however, closely related to the cause of the crashes generating a random, non trivial and semantically valid Kotlin program is hard. In this paper, we talk about type centric compiler fuzzing in the form of type centric enumeration, an approach inspired by skeletal program enumeration and based on a combination of generative and mutation based fuzzing, which solves this problem by focusing on program types. After creating the skeleton program, we fill the typed holes with fragments of suitable type, created via generation and enhanced by semantic aware mutation. We implemented this approach in our Kotlin compiler fuzzing framework called Backend Bug Finder (BBF) and did an extensive evaluation, not only testing the real world feasibility of our approach, but also comparing it to other compiler fuzzing techniques. The results show our approach to be significantly better compared to other fuzzing approaches at generating semantically valid Kotlin programs, while creating more interesting crash inducing inputs at the same time. We managed to find more than 50 previously unknown compiler crashes, of which 18 were considered important after their triage by the compiler team.


Dont Even Look Once Synthesizing Features for Zero Shot Detection. Zero shot detection, namely, localizing both seen and unseen objects, increasingly gains importance for large scale applications, with large number of object classes, since, collecting sufficient annotated data with ground truth bounding boxes is simply not scalable. While vanilla deep neural networks deliver high performance for objects available during training, unseen object detection degrades significantly. At a fundamental level, while vanilla detectors are capable of proposing bounding boxes, which include unseen objects, they are often incapable of assigning high confidence to unseen objects, due to the inherent precision recall tradeoffs that requires rejecting background objects. We propose a novel detection algorithm Dont Even Look Once (DELO), that synthesizes visual features for unseen objects and augments existing training algorithms to incorporate unseen object detection. Our proposed scheme is evaluated on Pascal VOC and MSCOCO, and we demonstrate significant improvements in test accuracy over vanilla and other state of art zero shot detectors


Resilient Optimal Estimation Using Measurement Prior. This paper considers the problem of optimal estimation for linear system with the measurement vector subject to arbitrary corruption by an adversarial agent. This problem is relevant to cyber physical systems where, due to the tight coupling of physics, communication and computation, a malicious agent is able to exploit multiple inherent vulnerabilities in order to inject stealthy signals into the measurement process. These malicious signals are calculated to serve the attack objectives of causing false situation awareness and or triggering a sequence of cascading effects leading to an ultimate system failure. We assume that the attacker can only compromise a portion, but not all, of the measurement channels simultaneously. However, once a channel is compromised, the attacker is free to modify the corresponding measurement arbitrarily. Consequently, the problem is formulated as a compressive sensing problem with additional prior information model. The prior information considered is a set inclusion constraint on the measurement vector. It is shown that if the prior set satisfies certain conditions, the resulting recovery error bound is much stronger. The approach is applied to the problem of resilient sate estimation of a power system. For this application, Gaussian Process is used to build a prior generative probabilistic regression model from historical data. The resulting Gaussian Process Regression model recursively maps energy market information to emph iid Gaussian distributions on the relevant system measurements. An optimization based resilient state estimator is then developed using a re weighted ell 1 minimization scheme. The developed algorithm is evaluated through a numerical simulation example of the IEEE 14 bus system mapped to the New York Independent System Operator (NYISO) grid data.


Deep correction of breathing related artifacts in real time MR thermometry. Real time MR imaging has been clinically adapted for monitoring thermal therapies since it can provide on the fly temperature maps simultaneously with anatomical information. However, proton resonance frequency based thermometry of moving targets remains challenging since temperature artifacts are induced by the respiratory as well as physiological motion. If left uncorrected, these artifacts lead to severe errors in temperature estimates and impair therapy guidance. In this study, we evaluated deep learning for on line correction of motion related errors in abdominal MR thermometry. For this, a convolutional neural network (CNN) was designed to learn the apparent temperature perturbation from images acquired during a preparative learning stage prior to hyperthermia. The input of the designed CNN is the most recent magnitude image and no surrogate of motion is needed. During the subsequent hyperthermia procedure, the recent magnitude image is used as an input for the CNN model in order to generate an on line correction for the current temperature map. The method's artifact suppression performance was evaluated on 12 free breathing volunteers and was found robust and artifact free in all examined cases. Furthermore, thermometric precision and accuracy was assessed for in vivo ablation using high intensity focused ultrasound. All calculations involved at the different stages of the proposed workflow were designed to be compatible with the clinical time constraints of a therapeutic procedure.


L infty BMO bounds for pseudo multipliers associated to the harmonic oscillator. In this note we investigate some conditions of H "ormander Mihlin type in order to assure the L infty textnormal BMO boundedness for pseudo multipliers of the harmonic oscillator. The textnormal H 1 L 1 continuity for Hermite multipliers also is investigated. The final version of this paper will appear in Rev. Colombiana Mat.


Critical elliptic equations on non compact Finsler manifolds. In the present paper, we deal with a quasilinear elliptic equation involving a critical Sobolev exponent on non compact Randers spaces. Under very general assumptions on the perturbation, we prove the existence of a non trivial solution. The approach is based on the direct methods of the calculus of variations. One of the key steps is to prove that the energy functional associated with the problem is weakly lower semicontinuous on small balls of the Sobolev space, which is provided by a general inequality. In the end, we prove Hardy type inequalities on Finsler manifolds as an application of this inequality.


Zero one laws for existential first order sentences of bounded quantifier depth. For any fixed positive integer k , let alpha k denote the smallest alpha in (0,1) such that the random graph sequence left G left(n, n alpha right) right does not satisfy the zero one law for the set mathcal E k of all existential first order sentences that are of quantifier depth at most k . This paper finds upper and lower bounds on alpha k , showing that as k rightarrow infty , we have alpha k left(k 2 t(k) right) 1 for some function t(k) Theta(k 2 ) . We also establish the precise value of alpha k when k 4 .


Confluence by Critical Pair Analysis Revisited (Extended Version). We present two methods for proving confluence of left linear term rewrite systems. One is hot decreasingness, combining the parallel development closedness theorems with rule labelling based on a terminating subsystem. The other is critical pair closing system, allowing to boil down the confluence problem to confluence of a special subsystem whose duplicating rules are relatively terminating.


Dynamic Role Based Access Control for Decentralized Applications. Access control management is an integral part of maintaining the security of an application. Although there has been significant work in the field of cloud access control mechanisms, however, with the advent of Distributed Ledger Technology (DLT), on chain access control management frameworks hardly exist. Existing access control management mechanisms are tightly coupled with the business logic, resulting in governance issues, non coherent with existing Identity Management Solutions, low security, and compromised usability. We propose a novel framework to implement dynamic role based access control for decentralized applications (dApps). The framework allows for managing access control on a dApp, which is completely decoupled from the business application and integrates seamlessly with any dApps. The smart contract architecture allows for the independent management of business logic and execution of access control policies. It also facilitates secure, low cost, and a high degree of flexibility of access control management. The proposed framework promotes decentralized governance of access control policies and efficient smart contract upgrades. We also provide quantitative and qualitative metrics for the efficacy and efficiency of the framework. Any Turing complete smart contract programming language is an excellent fit to implement the framework. We expect this framework to benefit enterprise and non enterprise dApps and provide greater access control flexibility and effective integration with traditional and state of the art identity management solutions.


Non local model for surface tension in fluid fluid simulations. We propose a non local model for surface tension obtained in the form of an integral of a molecular force like function with support 3.5 varepsilon added to the Navier Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference P varepsilon, in P varepsilon, out satisfies the Young Laplace law for droplet radius greater than 3.5 varepsilon and otherwise deviates from the Young Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier Stokes model subject to the Young Laplace boundary condition at the fluid fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid size dependent.


Martingale interpretation of weakly cancelling differential operators. We provide martingale analogs of weakly cancelling differential operators and prove a Sobolev type embedding theorem for these operators in the martingale setting.


A Review of Winograd Schema Challenge Datasets and Approaches. The Winograd Schema Challenge is both a commonsense reasoning and natural language understanding challenge, introduced as an alternative to the Turing test. A Winograd schema is a pair of sentences differing in one or two words with a highly ambiguous pronoun, resolved differently in the two sentences, that appears to require commonsense knowledge to be resolved correctly. The examples were designed to be easily solvable by humans but difficult for machines, in principle requiring a deep understanding of the content of the text and the situation it describes. This paper reviews existing Winograd Schema Challenge benchmark datasets and approaches that have been published since its introduction.


Collimated gamma ray beams from structured laser irradiated targets how to increase the efficiency without increasing the laser intensity. Using three dimensional kinetic simulations, we examine the emission of collimated gamma ray beams from structured laser irradiated targets with a pre filled cylindrical channel. The channel guides the incident laser pulse, enabling generation of a slowly evolving azimuthal plasma magnetic field that serves two key functions to enhance laser driven electron acceleration and to induce emission of gamma rays by the energetic electrons. Our main finding is that the conversion efficiency of the laser energy into a beam of gamma rays ( 5 circ opening angle) can be significantly increased without increasing the laser intensity by utilizing channels with an optimal density. The conversion efficiency into multi MeV photons increases roughly linearly with the incident laser power P , as we increase P from 1 PW to 4 PW while keeping the laser peak intensity fixed at 5 times 10 22 W cm 2 . This scaling is achieved by using an optimal range of plasma densities in the channel between 10 and 20 n cr , where n cr is the classical cutoff density for electromagnetic waves. The corresponding number of photons scales as P 2 . One application that directly benefits from such a strong scaling is the pair production via two photon collisions, with the number of generated pairs increasing as P 4 at fixed laser intensity.


The Hilbert scheme of a pair of linear spaces is a Mori dream space. Let H a,b n denote the component of the Hilbert scheme whose general point parameterizes an a plane union a b plane meeting transversely in mathbf P n . We describe the effective and nef cones of H a,b n and determine when the component is Fano. Moreover, we show that H a,b n is a Mori dream space for all values of a,b,n .


Unsupervised Transfer Learning with Self Supervised Remedy. Generalising deep networks to novel domains without manual labels is challenging to deep learning. This problem is intrinsically difficult due to unpredictable changing nature of imagery data distributions in novel domains. Pre learned knowledge does not transfer well without making strong assumptions about the learned and the novel domains. Different methods have been studied to address the underlying problem based on different assumptions, e.g. from domain adaptation to zero shot and few shot learning. In this work, we address this problem by transfer clustering that aims to learn a discriminative latent space of the unlabelled target data in a novel domain by knowledge transfer from labelled related domains. Specifically, we want to leverage relative (pairwise) imagery information, which is freely available and intrinsic to a target domain, to model the target domain image distribution characteristics as well as the prior knowledge learned from related labelled domains to enable more discriminative clustering of unlabelled target data. Our method mitigates nontransferrable prior knowledge by self supervision, benefiting from both transfer and self supervised learning. Extensive experiments on four datasets for image clustering tasks reveal the superiority of our model over the state of the art transfer clustering techniques. We further demonstrate its competitive transferability on four zero shot learning benchmarks.


Hybrid Approaches to Detect Comments Violating Macro Norms on Reddit. In this dataset paper, we present a three stage process to collect Reddit comments that are removed comments by moderators of several subreddits, for violating subreddit rules and guidelines. Other than the fact that these comments were flagged by moderators for violating community norms, we do not have any other information regarding the nature of the violations. Through this procedure, we collect over 2M comments removed by moderators of 100 different Reddit communities, and publicly release the data. Working with this dataset of removed comments, we identify 8 macro norms norms that are widely enforced on most parts of Reddit. We extract these macro norms by employing a hybrid approach classification, topic modeling, and open coding on comments identified to be norm violations within at least 85 out of the 100 study subreddits. Finally, we label over 40K Reddit comments removed by moderators according to the specific type of macro norm being violated, and make this dataset publicly available. By breaking down a collection of removed comments into more granular types of macro norm violation, our dataset can be used to train more nuanced machine learning classifiers for online moderation.


Injury risk increases minimally over a large range of changes in activity level in children. Background Limited research exists on the association between changes in physical activity levels and injury in children. Objective To assess how well different variations of the acute chronic workload ratio (ACWR), a measure of change in activity, predict injury in children. Methods We conducted a prospective cohort study using data from 1670 Danish schoolchildren measured over 5.5 years (2008 to 2014). Coupled 4 week, uncoupled 4 week, and uncoupled 5 week ACWRs were calculated using activity frequency in the past week as the acute load (numerator), and average weekly activity frequency in the past 4 or 5 weeks as the chronic load (denominator). We modelled the relationship between different ACWR variations and injury using generalized linear and generalized additive models, with and without accounting for repeated measures. Results The prognostic relationship between the ACWR and injury risk was best represented using a generalized additive mixed model for the uncoupled 5 week ACWR. It predicted an injury risk of 3 for ACWRs between 0.8 (activity level decreased by 20 ) and 1.5 (activity level increased by 50 ). When activity decreased by more than 20 (ACWR 0.8), injury risk was lower (minimum of 1.5 at ACWR 0). When activity increased by more than 50 (ACWR 1.5), injury risk was higher (maximum of 6 at ACWR 5). Girls were at significantly higher risk of injury than boys. Conclusion Increases in physical activity in children are associated with much lower injury risks compared to previous results in adults.


Existence of an equilibrium for pure exchange economy with fuzzy preferences. This paper focuses on a new model to reach the existence of equilibrium in a pure exchange economy with fuzzy preferences (PXE FP). The proposed model integrates exchange, consumption and the agent's fuzzy preference in the consumption set. We set up a new fuzzy binary relation on the consumption set to evaluate the fuzzy preferences. Also, we prove that there exists a continuous fuzzy order preserving function in the consumption set under certain conditions. The existence of a fuzzy competitive equilibrium for the PXE FP is confirmed through a new result on the existence of fuzzy Nash equilibrium for fuzzy non cooperative games. The payoffs of all strategy profiles for any agent are fuzzy numbers in fuzzy non cooperative games. Finally, we show that the fuzzy competitive equilibrium could be characterized as a solution to an associated quasi variational inequality, giving rise to an equilibrium solution.


Primigraph spaces. In this paper we introduce primigraph spaces, which are topological spaces together with a sheaf of C algebras that can be covered by some Prim A's, that is, by the primitive spectra of some C algebras endowed with Jacobson topology and together with the sheaf of bounded continuous functions on them. This notion is analoguous, in some sense, to the ones of schemes and perfectoid spaces. Our first main result here is that every topological space gives rise to a primigraph space this will imply, thanks to Theorem 4.1, that primigraph spaces also constitute a new kind of topological invariant.