6TiSCH with Bluetooth 5 and Concurrent Transmissions. Targeting dependable communications for industrial Internet of Things applications, IETF 6TiSCH provides mechanisms for efficient scheduling, routing, and forwarding of IPv6 traffic across low power mesh networks. Yet, despite an overwhelming body of literature covering both centralized and distributed scheduling schemes for 6TiSCH, an effective control solution for large scale multi hop mesh networks remains an open challenge. This paper addresses this with a novel approach that eliminates much of the routing and link layer overhead incurred by centralized schedulers, and provides a robust mechanism for data dissemination synchronization within 6TiSCH. Specifically, we leverage the physical layer (PHY) switching capabilities of modern low power wireless platforms to build on recent work demonstrating the viability of Concurrent Transmission (CT) based flooding protocols across the Bluetooth 5 (BT 5) PHYs. By switching the PHY and MAC layer at runtime, we inject a BT 5 based CT flood within a standard IEEE 802.15.4 TSCH slotframe, thereby providing a reliable, low latency scheme for 6TiSCH control messaging. We present an analytical model and experimental evaluation showing how our solution not only exploits the BT 5 high data rate PHY layers for rapid data dissemination, but can also provide reliable 6TiSCH association and synchronization even under external radio interference. We further discuss how the proposed technique can be used to address other open challenges within the standard.


Exploration of Fine Grained Parallelism for Load Balancing Eager K truss on GPU and CPU. In this work we present a performance exploration on Eager K truss, a linear algebraic formulation of the K truss graph algorithm. We address performance issues related to load imbalance of parallel tasks in symmetric, triangular graphs by presenting a fine grained parallel approach to executing the support computation. This approach also increases available parallelism, making it amenable to GPU execution. We demonstrate our fine grained parallel approach using implementations in Kokkos and evaluate them on an Intel Skylake CPU and an Nvidia Tesla V100 GPU. Overall, we observe between a 1.261. 48x improvement on the CPU and a 9.97 16.92x improvement on the GPU due to our fine grained parallel formulation.


Projective normality of torus quotients of flag varieties. Let G SL n( mathbb C) and T be a maximal torus in G . We show that the quotient T backslash backslash G P alpha 1 cap P alpha 2 is projectively normal with respect to the descent of a suitable line bundle, where P alpha i is the maximal parabolic subgroup in G associated to the simple root alpha i , i 1,2 . We give a degree bound of the generators of the homogeneous coordinate ring of T backslash backslash (G 3,6 ) ss T( mathcal L 2 varpi 3 ) . If G Spin 7 , we give a degree bound of the generators of the homogeneous coordinate ring of T backslash backslash (G P alpha 2 ) ss T( mathcal L 2 varpi 2 ) whereas we prove that the quotient T backslash backslash (G P alpha 3 ) ss T( mathcal L 4 varpi 3 ) is projectively normal with respect to the descent of the line bundles mathcal L 4 varpi 3 .


Time smoothing for parabolic variational problems in metric measure spaces. In 2013, Masson and Siljander determined a method to prove that the p minimal upper gradient g f varepsilon for the time mollification f varepsilon , varepsilon 0 , of a parabolic Newton Sobolev function f in L p mathrm loc (0, tau N 1,p mathrm loc ( Omega)) , with tau 0 and Omega open domain in a doubling metric measure space ( mathbb X ,d, mu) supporting a weak (1,p) Poincar 'e inequality, p in(1, infty) , is such that g f f varepsilon to0 as varepsilon to0 in L p mathrm loc ( Omega tau) , Omega tau being the parabolic cylinder Omega tau Omega times(0, tau) . Their approach involved the use of Cheeger's differential structure, and therefore exhibited some limitations here, we shall see that the definition and the formal properties of the parabolic Sobolev spaces themselves allow to find a more direct method to show such convergence, which relies on p weak upper gradients only and which is valid regardless of structural assumptions on the ambient space, also in the limiting case when p 1 .


Concentration of Distortion The Value of Extra Voters in Randomized Social Choice. We study higher statistical moments of Distortion for randomized social choice in a metric implicit utilitarian model. The Distortion of a social choice mechanism is the expected approximation factor with respect to the optimal utilitarian social cost (OPT). The k th moment of Distortion is the expected approximation factor with respect to the k th power of OPT. We consider mechanisms that elicit alternatives by randomly sampling voters for their favorite alternative. We design two families of mechanisms that provide constant (with respect to the number of voters and alternatives) k th moment of Distortion using just k samples if all voters can then participate in a vote among the proposed alternatives, or 2k 1 samples if only the sampled voters can participate. We also show that these numbers of samples are tight. Such mechanisms deviate from a constant approximation to OPT with probability that drops exponentially in the number of samples, independent of the total number of voters and alternatives. We conclude with simulations on real world Participatory Budgeting data to qualitatively complement our theoretical insights.


Nonexistence of variational minimizers related to a quasilinear singular problem in metric measure spaces. In this article we consider a variational problem related to a quasilinear singular problem and obtain a nonexistence result in a metric measure space with a doubling measure and a Poincar 'e inequality. Our method is purely variational and to the best of our knowledge, this is the first work concerning singular problems in a general metric setting.


Hybridized Summation By Parts Finite Difference Methods. We present a hybridization technique for summation by parts finite difference methods with weak enforcement of interface and boundary conditions for second order, linear elliptic partial differential equations. The method is based on techniques from the hybridized discontinuous Galerkin literature where local and global problems are defined for the volume and trace grid points, respectively. By using a Schur complement technique the volume points can be eliminated, which drastically reduces the system size. We derive both the local and global problems, and show that the linear systems that must be solved are symmetric positive definite. The theoretical stability results are confirmed with numerical experiments as is the accuracy of the method.


The degree of the central curve in semidefinite, linear, and quadratic programming. The Zariski closure of the central path which interior point algorithms track in convex optimization problems such as linear, quadratic, and semidefinite programs is an algebraic curve. The degree of this curve has been studied in relation to the complexity of these interior point algorithms, and for linear programs it was computed by De Loera, Sturmfels, and Vinzant in 2012. We show that the degree of the central curve for generic semidefinite programs is equal to the maximum likelihood degree of linear concentration models. New results from the intersection theory of the space of complete quadrics imply that this is a polynomial in the size of semidefinite matrices with degree equal to the number of constraints. Besides its degree we explore the arithmetic genus of the same curve. We also compute the degree of the central curve for generic linear programs with different techniques which extend to bounding the same degree for generic quadratic programs.


Improving Tractability of Real Time Control Schemes via Simplified mathcal S Lemma. Various control schemes rely on a solution of a convex optimization problem involving a particular robust quadratic constraint, which can be reformulated as a linear matrix inequality using the well known mathcal S lemma. However, the computational effort required to solve the resulting semidefinite program may be prohibitively large for real time applications requiring a repeated solution of such a problem. We use some recent advances in robust optimization that allow us to reformulate such a robust constraint as a set of linear and second order cone constraints, which are computationally better suited to real time applications. A numerical example demonstrates a huge speedup that can be obtained using the proposed reformulation.


On Non uniqueness of H " o lder continuous globally dissipative Euler flows. We show that for any al frac 17 there exist al H "older continuous weak solutions of the three dimensional incompressible Euler equation, which satisfy the local energy inequality and strictly dissipate the total kinetic energy. The proof relies on the convex integration scheme and the main building blocks of the solution are various Mikado flows with disjoint supports in space and time.


Achievements and Challenges in Explaining Deep Learning based Computer Aided Diagnosis Systems. Remarkable success of modern image based AI methods and the resulting interest in their applications in critical decision making processes has led to a surge in efforts to make such intelligent systems transparent and explainable. The need for explainable AI does not stem only from ethical and moral grounds but also from stricter legislation around the world mandating clear and justifiable explanations of any decision taken or assisted by AI. Especially in the medical context where Computer Aided Diagnosis can have a direct influence on the treatment and well being of patients, transparency is of utmost importance for safe transition from lab research to real world clinical practice. This paper provides a comprehensive overview of current state of the art in explaining and interpreting Deep Learning based algorithms in applications of medical research and diagnosis of diseases. We discuss early achievements in development of explainable AI for validation of known disease criteria, exploration of new potential biomarkers, as well as methods for the subsequent correction of AI models. Various explanation methods like visual, textual, post hoc, ante hoc, local and global have been thoroughly and critically analyzed. Subsequently, we also highlight some of the remaining challenges that stand in the way of practical applications of AI as a clinical decision support tool and provide recommendations for the direction of future research.


Influence of Electromagnetic Fields on Nuclear Processes. Although the energies associated with nuclear reactions are due primarily to interactions involving nuclear forces, the rates and probabilities associated with those reactions are effectively governed by electromagnetic forces. Charges in the local environment can modulate the Coulomb barrier, and thereby change the rates of nuclear processes. Experiments are presented in which low temperature electrons are attached to high density rotating neutrals to form negative ions. The steady state quiescent rotating plasma generated here lends itself to first prove the principle that low temperature systems can yield MeV fusion particles. It allows the use of high density of neutrals interacting with the wall to yield gain greater than unity. It also demonstrates that instabilities can be avoided with high neutral densities. Collective dynamic interactions within this steady state quiescent plasma result in an arrangement of negative charges that lowers the effective Coulomb barrier to nuclear reactions at a solid wall of reactants. MeV alpha particles are synchronously observed with externally imposed pulses as evidence of fusion being enabled by Coulomb fields. Impacts on fusion, the source of energy in the universe, will be discussed.


SNR adaptive OCT angiography enabled by statistical characterization of intensity and decorrelation with multi variate time series model. In OCT angiography (OCTA), decorrelation computation has been widely used as a local motion index to identify dynamic flow from static tissues, but its dependence on SNR severely degrades the vascular visibility, particularly in low SNR regions. To mathematically characterize the decorrelation SNR dependence of OCT signals, we developed a multi variate time series (MVTS) model. Based on the model, we derived a universal asymptotic linear relation of decorrelation to inverse SNR (iSNR), with the variance in static and noise regions determined by the average kernel size. Accordingly, with the population distribution of static and noise voxels being explicitly calculated in the iSNR and decorrelation (ID) space, a linear classifier is developed by removing static and noise voxels at all SNR, to generate a SNR adaptive OCTA, termed as ID OCTA. Then, flow phantom and human skin experiments were performed to validate the proposed ID OCTA. Both qualitative and quantitative assessments demonstrated that ID OCTA offers a superior visibility of blood vessels, particularly in the deep layer. Finally, implications of this work on both system design and hemodynamic quantification are further discussed.


On the Strength of Connectivity of Inhomogeneous Random K out Graphs. Random graphs are an important tool for modelling and analyzing the underlying properties of complex real world networks. In this paper, we study a class of random graphs known as the inhomogeneous random K out graphs which were recently introduced to analyze heterogeneous sensor networks secured by the pairwise scheme. In this model, first, each of the n nodes is classified as type 1 (respectively, type 2) with probability 0 mu 1 (respectively, 1 mu) independently from each other. Next, each type 1 (respectively, type 2) node draws 1 arc towards a node (respectively, K n arcs towards K n distinct nodes) selected uniformly at random, and then the orientation of the arcs is ignored. From the literature on homogeneous K out graphs wherein all nodes select K n neighbors (i.e., mu 0 ), it is known that when K n geq2 , the graph is K n connected asymptotically almost surely (a.a.s.) as n gets large. In the inhomogeneous case (i.e., mu 0 ), it was recently established that achieving even 1 connectivity a.a.s. requires K n omega(1) . Here, we provide a comprehensive set of results to complement these existing results. First, we establish a sharp zero one law for k connectivity, showing that for the network to be k connected a.a.s., we need to set K n frac 1 1 mu ( log n (k 2) log log n omega(1)) for all k 2, 3, ldots . Despite such large scaling of K n being required for k connectivity, we show that the trivial condition of K n geq 2 for all n is sufficient to ensure that inhomogeneous K out graph has a connected component of size n O(1) whp.


Beto, Bentz, Becas The Surprising Cross Lingual Effectiveness of BERT. Pretrained contextual representation models (Peters et al., 2018 Devlin et al., 2018) have pushed forward the state of the art on many NLP tasks. A new release of BERT (Devlin, 2018) includes a model simultaneously pretrained on 104 languages with impressive performance for zero shot cross lingual transfer on a natural language inference task. This paper explores the broader cross lingual potential of mBERT (multilingual) as a zero shot language transfer model on 5 NLP tasks covering a total of 39 languages from various language families NLI, document classification, NER, POS tagging, and dependency parsing. We compare mBERT with the best published methods for zero shot cross lingual transfer and find mBERT competitive on each task. Additionally, we investigate the most effective strategy for utilizing mBERT in this manner, determine to what extent mBERT generalizes away from language specific features, and measure factors that influence cross lingual transfer.


Explicit and recursive estimates of the Lambert W function. Solutions to a wide variety of transcendental equations can be expressed in terms of the Lambert mathrm W function. The mathrm W function, occurring frequently in applications, is a non elementary, but now standard mathematical function implemented in all major technical computing systems. In this work, we discuss some approximations of the two real branches, mathrm W 0 and mathrm W 1 . On the one hand, we present some analytic lower and upper bounds on mathrm W 0 for large arguments that improve on some earlier results in the literature. On the other hand, we analyze two logarithmic recursions, one with linear, and the other with quadratic rate of convergence. We propose suitable starting values for the recursion with quadratic rate that ensure convergence on the whole domain of definition of both real branches. We also provide a priori, simple, explicit and uniform estimates on its convergence speed that enable guaranteed, high precision approximations of mathrm W 0 and mathrm W 1 at any point. Finally, as an application of the mathrm W 0 function, we settle a conjecture about the growth rate of the positive non trivial solutions to the equation x y y x .


Microlocally accurate solution of the inverse source problem of thermoacoustic tomography. We consider the inverse source problem of thermo and photoacoustic tomography, with data registered on an open surface partially surrounding the source of acoustic waves. Under the assumption of constant speed of sound we develop an explicit non iterative reconstruction procedure that recovers the Radon transform of the sought source, up to an infinitely smooth additive error term. The source then can be found by inverting the Radon transform. Our analysis is microlocal in nature and does not provide a norm estimate on the error in the so obtained image. However, numerical simulations show that this error is quite small in practical terms. We also present an asymptotically fast implementation of this procedure for the case when the data are given on a circular arc in 2D.


Asymptotics for Hecke eigenvalues of automorphic forms on compact arithmetic quotients. In this paper, we describe the asymptotic distribution of Hecke eigenvalues in the Laplace eigenvalue aspect for certain families of Hecke Maass forms on compact arithmetic quotients. Instead of relying on the trace formula, which was the primary tool in preceding studies on the subject, we use Fourier integral operator methods. This allows us to treat not only spherical, but also non spherical Hecke Maass forms with corresponding remainder estimates. Our asymptotic formulas are available for arbitrary simple and connected algebraic groups over number fields with cocompact arithmetic subgroups.


SuperFront From Low resolution to High resolution Frontal Face Synthesis. Advances in face rotation, along with other face based generative tasks, are more frequent as we advance further in topics of deep learning. Even as impressive milestones are achieved in synthesizing faces, the importance of preserving identity is needed in practice and should not be overlooked. Also, the difficulty should not be more for data with obscured faces, heavier poses, and lower quality. Existing methods tend to focus on samples with variation in pose, but with the assumption data is high in quality. We propose a generative adversarial network (GAN) based model to generate high quality, identity preserving frontal faces from one or multiple low resolution (LR) faces with extreme poses. Specifically, we propose SuperFront GAN (SF GAN) to synthesize a high resolution (HR), frontal face from one to many LR faces with various poses and with the identity preserved. We integrate a super resolution (SR) side view module into SF GAN to preserve identity information and fine details of the side views in HR space, which helps model reconstruct high frequency information of faces (i.e., periocular, nose, and mouth regions). Moreover, SF GAN accepts multiple LR faces as input, and improves each added sample. We squeeze additional gain in performance with an orthogonal constraint in the generator to penalize redundant latent representations and, hence, diversify the learned features space. Quantitative and qualitative results demonstrate the superiority of SF GAN over others.


Normal Reflection Subgroups. We study normal reflection subgroups of complex reflection groups. Our point of view leads to a refinement of a theorem of Orlik and Solomon to the effect that the generating function for fixed space dimension over a reflection group is a product of linear factors involving generalized exponents. Our refinement gives a uniform proof and generalization of a recent theorem of the second author.


Towards an efficient and exact algorithm for dynamic dedicated path protection. We present a novel algorithm for dynamic routing with dedicated path protection which, as the presented simulation results suggest, can be efficient and exact. We present the algorithm in the setting of optical networks, but it should be applicable to other networks, where services have to be protected, and where the network resources are finite and discrete, e.g., wireless radio or networks capable of advance resource reservation. To the best of our knowledge, we are the first to propose an algorithm for this long standing fundamental problem, which can be efficient and exact, as suggested by simulation results. The algorithm can be efficient because it can solve large problems, and it can be exact because its results are optimal, as demonstrated and corroborated by simulations. We offer a worst case analysis to argue that the search space is polynomially upper bounded. Network operations, management, and control require efficient and exact algorithms, especially now, when greater emphasis is placed on network performance, reliability, softwarization, agility, and return on investment. The proposed algorithm uses our generic Dijkstra algorithm on a search graph generated "on the fly" based on the input graph. We corroborated the optimality of the results of the proposed algorithm with brute force enumeration for networks up to 15 nodes large. We present the extensive simulation results of dedicated path protection with signal modulation constraints for elastic optical networks of 25, 50, and 100 nodes, and with 160, 320, and 640 spectrum units. We also compare the bandwidth blocking probability with the commonly used edge exclusion algorithm. We had 48,600 simulation runs with about 41 million searches.


Hardy Littlewood inequality and L p L q Fourier multipliers on compact hypergroups. This paper deals with the inequalities devoted to the comparison between the norm of a function on a compact hypergroup and the norm of its Fourier coefficients. We prove the classical Paley inequality in the setting of compact hypergroups which further gives the Hardy Littlewood and Hausdorff Young Paley (Pitt) inequalities in the noncommutative context. We establish H "ormander's L p L q Fourier multiplier theorem on compact hypergroups for 1 p leq 2 leq q infty as an application of Hausdorff Young Paley inequality. We examine our results for the hypergroups constructed from the conjugacy classes of compact Lie groups and for a class of countable compact hypergroups.


Using circular dichroism to control energy transfer in multi photon ionization. Chirality causes symmetry breaks in a large variety of natural phenomena ranging from particle physics to biochemistry. We investigate one of the simplest conceivable chiral systems, a laser excited, oriented, effective one electron Li target. Prepared in a polarized p state with m 1 in an optical trap, the atoms are exposed to co and counter rotating circularly polarized femtosecond laser pulses. For a field frequency near the excitation energy of the oriented initial state, a strong circular dichroism is observed and the photoelectron energies are significantly affected by the helicity dependent Autler Townes splitting. Besides its fundamental relevance, this system is suited to create spin polarized electron pulses with a reversible switch on a femtosecond timescale at an energy resolution of a few meV.


Geometric Variations of Local Systems and Elliptic Surfaces. Geometric variations of local systems are families of variations of Hodge structure they typically correspond to fibrations of K " a hler manifolds for which each fibre itself is fibred by codimension one K " a hler manifolds. In this article, we introduce the formalism of geometric variations of local systems and then specialize the theory to study families of elliptic surfaces. We interpret a construction of twisted elliptic surface families used by Besser Livn ' e in terms of the middle convolution functor, and use explicit methods to calculate the variations of Hodge structure underlying the universal families of M N polarized K3 surfaces. Finally, we explain the connection between geometric variations of local systems and geometric isomonodromic deformations, which were originally considered by the first author in 1999.


Efficient Alternating Least Squares Algorithms for Low Multilinear Rank Approximation of Tensors. The low multilinear rank approximation, also known as the truncated Tucker decomposition, has been extensively utilized in many applications that involve higher order tensors. Popular methods for low multilinear rank approximation usually rely directly on matrix SVD, therefore often suffer from the notorious intermediate data explosion issue and are not easy to parallelize, especially when the input tensor is large. In this paper, we propose a new class of truncated HOSVD algorithms based on alternating least squares (ALS) for efficiently computing the low multilinear rank approximation of tensors. The proposed ALS based approaches are able to eliminate the redundant computations of the singular vectors of intermediate matrices and are therefore free of data explosion. Also, the new methods are more flexible with adjustable convergence tolerance and are intrinsically parallelizable on high performance computers. Theoretical analysis reveals that the ALS iteration in the proposed algorithms is q linear convergent with a relatively wide convergence region. Numerical experiments with large scale tensors from both synthetic and real world applications demonstrate that ALS based methods can substantially reduce the total cost of the original ones and are highly scalable for parallel computing.


A Novel Privacy Preserved Recommender System Framework based on Federated Learning. Recommender System (RS) is currently an effective way to solve information overload. To meet users' next click behavior, RS needs to collect users' personal information and behavior to achieve a comprehensive and profound user preference perception. However, these centrally collected data are privacy sensitive, and any leakage may cause severe problems to both users and service providers. This paper proposed a novel privacy preserved recommender system framework (PPRSF), through the application of federated learning paradigm, to enable the recommendation algorithm to be trained and carry out inference without centrally collecting users' private data. The PPRSF not only able to reduces the privacy leakage risk, satisfies legal and regulatory requirements but also allows various recommendation algorithms to be applied.


Towards combined optical coherence tomography and multi spectral imaging with MHz a scan rates for endoscopy. We demonstrate a preliminary setup of a combined MHz OCT and RGB narrowband reflection microscope and investigate the performance of the new RGB branch and different display modes of colored OCT data sets.


SDRSAC Semidefinite Based Randomized Approach for Robust Point Cloud Registration without Correspondences. This paper presents a novel randomized algorithm for robust point cloud registration without correspondences. Most existing registration approaches require a set of putative correspondences obtained by extracting invariant descriptors. However, such descriptors could become unreliable in noisy and contaminated settings. In these settings, methods that directly handle input point sets are preferable. Without correspondences, however, conventional randomized techniques require a very large number of samples in order to reach satisfactory solutions. In this paper, we propose a novel approach to address this problem. In particular, our work enables the use of randomized methods for point cloud registration without the need of putative correspondences. By considering point cloud alignment as a special instance of graph matching and employing an efficient semi definite relaxation, we propose a novel sampling mechanism, in which the size of the sampled subsets can be larger than minimal. Our tight relaxation scheme enables fast rejection of the outliers in the sampled sets, resulting in high quality hypotheses. We conduct extensive experiments to demonstrate that our approach outperforms other state of the art methods. Importantly, our proposed method serves as a generic framework which can be extended to problems with known correspondences.


Geometry Inspired Top k Adversarial Perturbations. The brittleness of deep image classifiers to small adversarial input perturbations has been extensively studied in the last several years. However, the main objective of existing perturbations is primarily limited to change the correctly predicted Top 1 class by an incorrect one, which does not intend to change the Top k prediction. In many digital real world scenarios Top k prediction is more relevant. In this work, we propose a fast and accurate method of computing Top k adversarial examples as a simple multi objective optimization. We demonstrate its efficacy and performance by comparing it to other adversarial example crafting techniques. Moreover, based on this method, we propose Top k Universal Adversarial Perturbations, image agnostic tiny perturbations that cause the true class to be absent among the Top k prediction for the majority of natural images. We experimentally show that our approach outperforms baseline methods and even improves existing techniques of finding Universal Adversarial Perturbations.


The Typical Structure of Sets with Small Sumset. In this paper we determine the number and typical structure of sets of integers with bounded doubling. In particular, improving recent results of Green and Morris, and of Mazur, we show that the following holds for every fixed lambda 2 and every k geqslant ( log n) 4 if omega to infty as n to infty (arbitrarily slowly), then almost all sets A subset n with A k and A A leqslant lambda k are contained in an arithmetic progression of length lambda k 2 omega .


Inferring and Improving Street Maps with Data Driven Automation. Street maps are a crucial data source that help to inform a wide range of decisions, from navigating a city to disaster relief and urban planning. However, in many parts of the world, street maps are incomplete or lag behind new construction. Editing maps today involves a tedious process of manually tracing and annotating roads, buildings, and other map features. Over the past decade, many automatic map inference systems have been proposed to automatically extract street map data from satellite imagery, aerial imagery, and GPS trajectory datasets. However, automatic map inference has failed to gain traction in practice due to two key limitations high error rates (low precision), which manifest in noisy inference outputs, and a lack of end to end system design to leverage inferred data to update existing street maps. At MIT and QCRI, we have developed a number of algorithms and approaches to address these challenges, which we combined into a new system we call Mapster. Mapster is a human in the loop street map editing system that incorporates three components to robustly accelerate the mapping process over traditional tools and workflows high precision automatic map inference, data refinement, and machine assisted map editing. Through an evaluation on a large scale dataset including satellite imagery, GPS trajectories, and ground truth map data in forty cities, we show that Mapster makes automation practical for map editing, and enables the curation of map datasets that are more complete and up to date at less cost.


Specializing Coherence, Consistency, and Push Pull for GPU Graph Analytics. This work provides the first study to explore the interaction of update propagation with and without fine grained synchronization (push vs. pull), emerging coherence protocols (GPU vs. DeNovo coherence), and software centric consistency models (DRF0, DRF1, and DRFrlx) for graph workloads on emerging integrated GPU CPU systems with native unified shared memory. We study 6 graph applications with 6 graph inputs for a total of 36 workloads running on 12 system (hardware software) configurations reflecting the above design space of update propagation, coherence, and memory consistency. We make three key contributions. First, we show that there is no single best system configuration for all workloads, motivating systems with flexible coherence and consistency support. Second, we develop a model to accurately predict the best system configuration this model can be used by software designers to decide on push vs. pull and the consistency model and by flexible hardware to invoke the appropriate coherence and consistency configuration for the given workload. Third, we show that the design dimensions explored here are inter dependent, reinforcing the need for software hardware co design in the above design dimensions. For example, software designers deciding on push vs. pull must consider the consistency model supported by hardware in some cases, push maybe better if hardware supports DRFrlx while pull may be better if hardware does not support DRFrlx.


CAIL2019 SCM A Dataset of Similar Case Matching in Legal Domain. In this paper, we introduce CAIL2019 SCM, Chinese AI and Law 2019 Similar Case Matching dataset. CAIL2019 SCM contains 8,964 triplets of cases published by the Supreme People's Court of China. CAIL2019 SCM focuses on detecting similar cases, and the participants are required to check which two cases are more similar in the triplets. There are 711 teams who participated in this year's competition, and the best team has reached a score of 71.88. We have also implemented several baselines to help researchers better understand this task. The dataset and more details can be found from https github.com china ai law challenge CAIL2019 tree master scm.


Spatial Social Network (SSN) Hot Spot Detection Scan Methods for Non Planar Networks. Moving window and hot spot detection analyses are statistical methods used to analyze point patterns within a given area. Such methods have been used to successfully detect clusters of point events such as car thefts or incidences of cancer. Yet, these methods do not account for the connections between individual events, such as social ties within a neighborhood. This paper presents two GIS methods, EdgeScan and NDScan, for capturing areas with high and low levels of local social connections. Both methods are moving window processes that count the number of edges and network density, respectively, in a given focal area (window area). The focal window attaches resultant EdgeScan and NDScan statistics to nodes at the center of the focal window area. We implement these methods on a case study of 1960s connections between members of the Mafia in New York City. We use various definitions of a focal neighborhood including Euclidean, Manhattan and K Nearest Neighbor (KNN) definitions. We find that KNN tends to overstate the values of local networks, and that there is more variation in outcome values for nodes on the periphery of the study area. We find that, location wise, EdgeScan and NDScan hot spots differ from traditional spatial hot spots in the study area. These methods can be extended to future studies that detect local triads and motifs, which can capture the local network structure in more detail.


An Information theoretic Approach to Unsupervised Feature Selection for High Dimensional Data. In this paper, we propose an information theoretic approach to design the functional representations to extract the hidden common structure shared by a set of random variables. The main idea is to measure the common information between the random variables by Watanabe's total correlation, and then find the hidden attributes of these random variables such that the common information is reduced the most given these attributes. We show that these attributes can be characterized by an exponential family specified by the eigen decomposition of some pairwise joint distribution matrix. Then, we adopt the log likelihood functions for estimating these attributes as the desired functional representations of the random variables, and show that such representations are informative to describe the common structure. Moreover, we design both the multivariate alternating conditional expectation (MACE) algorithm to compute the proposed functional representations for discrete data, and a novel neural network training approach for continuous or high dimensional data. Furthermore, we show that our approach has deep connections to existing techniques, such as Hirschfeld Gebelein R ' e nyi (HGR) maximal correlation, linear principal component analysis (PCA), and consistent functional map, which establishes insightful connections between information theory and machine learning. Finally, the performances of our algorithms are validated by numerical simulations.


Forced convection in a fluid saturated anisotropic porous channel with isoflux boundaries. Fully developed forced convective flow in a channel filled with porous material bounded by two impermeable walls subject to constant heat flux is considered. We consider Brinkman Forchhimer equation to govern the flow inside the porous medium which accounts for the presence of inertial term. We assume that the porous medium is anisotropic in nature and permeability is varying along all the directions, correspondingly, permeability will be a matrix which appear as a positive semi definite matrix in the momentum equation. The anisotropic nature of the porous medium is characterized by anisotropic ratio and angle. We have obtained the velocity, temperature and Nusselt number numerically, due to the presence of the non linear quadratic term in the momentum equation. The effect of anisotropic permeability ratio and anisotropic angle on hydrodynamic and heat transfer is shown. We see that Nusselt number is altered significantly with the anisotropic ratio and angle. We present a fresh analysis about the inclusion of the permeability matrix in the Brinkman Forchhimer extended Darcy momentum equation. To the best of the authors knowledge, this is the first attempt to analyze general anisotropy corresponding to Brinkman Forchhimer extended Darcy model.


Using an Educational IoT Lab Kit and Gamification for Energy Awareness in European Schools. The use of maker community tools and IoT technologies inside classrooms is spreading in an increasing number of education and science fields. GAIA is a European research project focused on achieving behavior change for sustainability and energy awareness in schools. In this work, we report on how a large IoT deployment in a number of educational buildings and real world data from this infrastructure, are utilized to support a "maker" lab kit activity inside the classroom, together with a serious game. We also provide some insights to the integration of these activities in the school curriculum, along with a discussion on our feedback so far from a series of workshop activities in a number of schools. Our initial results show strong acceptance by the school community.


Reconstructing Perceived Images from Brain Activity by Visually guided Cognitive Representation and Adversarial Learning. Reconstructing visual stimulus (image) only from human brain activity measured with functional Magnetic Resonance Imaging (fMRI) is a significant and meaningful task in Human AI collaboration. However, the inconsistent distribution and representation between fMRI signals and visual images cause the heterogeneity gap. Moreover, the fMRI data is often extremely high dimensional and contains a lot of visually irrelevant information. Existing methods generally suffer from these issues so that a satisfactory reconstruction is still challenging. In this paper, we show that it is possible to overcome these challenges by learning visually guided cognitive latent representations from the fMRI signals, and inversely decoding them to the image stimuli. The resulting framework is called Dual Variational Autoencoder Generative Adversarial Network (D VAE GAN), which combines the advantages of adversarial representation learning with knowledge distillation. In addition, we introduce a novel three stage learning approach which enables the (cognitive) encoder to gradually distill useful knowledge from the paired (visual) encoder during the learning process. Extensive experimental results on both artificial and natural images have demonstrated that our method could achieve surprisingly good results and outperform all other alternatives.


Congruence properties of coefficients of the eighth order mock theta function V 0(q) . We study the divisibility properties of the partition function associated with the eighth order mock theta function V 0(q) , introduced by Gordon and McIntosh. We obtain congruences modulo powers of 2 for certain coefficients of the partition function, akin to Ramanujan's partition congruences. Further, we also present several infinite families of congruences molulo 13, 25 and 27.


A Multiplicity Result for Orthogonal Geodesic Chords in Finsler disks. In this paper, we study the existence and multiplicity problems for orthogonal Finsler geodesic chords in a manifold with boundary which is homeomorphic to a N dimensional disk. Under a suitable assumption, which is weaker than convexity, we prove that, if the Finsler metric is reversible, then there are at least N orthogonal Finsler geodesic chords that are geometrically distinct. If the reversibility assumption does not hold, then there are at least two orthogonal Finsler geodesic chords with different values of the energy functional.


Suprathermal Proton Spectra at Interplanetary Shocks in Hybrid Simulations. Interplanetary shocks are one of the proposed sources of suprathermal ion populations (i.e., ions with energies of a few times the solar wind energy). Here, we present results from a series of three dimensional hybrid simulations of collisionless shocks in the solar wind. We focus on the influence of the shock normal angle, theta Bn , and the shock speed, V s , on producing protons with energies a few to hundreds of times the thermal energy of the upstream plasma. The combined effects of theta Bn and V s result in shocks with Alfv 'en Mach numbers in the range 3.0 to 6.0 and fast magnetosonic Mach numbers in the range 2.5 to 5.0, representing moderate to strong interplanetary shocks. We find that theta Bn largely organizes the shape of proton energy spectra while shock speed controls acceleration efficiency. All shocks accelerate protons at the shock front but the spectral evolution depends on theta Bn . Shocks with theta Bn geq 60 circ produce isolated bursts of suprathermal protons at the shock front while shocks with theta Bn leq 45 circ create suprathermal beams upstream of the shock. Downstream proton energy spectra have exponential or smoothed broken power law forms when theta Bn geq 45 circ , and a single power law form when theta Bn leq 30 circ . Protons downstream of the strongest shocks have energies at least 100 times the upstream thermal energy, with theta Bn leq 30 circ shocks producing the highest energy protons and theta Bn geq 60 circ shocks producing the largest number of protons with energies at least a few times the thermal energy.


Multi scale sparse domination. We prove a bilinear form sparse domination theorem that applies to many multi scale operators beyond Calder 'on Zygmund theory, and also establish necessary conditions. Among the applications, we cover large classes of Fourier multipliers, maximal functions, square functions and variation norm operators.


Non destructive three dimensional measurement of hand vein based on self supervised network. At present, supervised stereo methods based on deep neural network have achieved impressive results. However, in some scenarios, accurate three dimensional labels are inaccessible for supervised training. In this paper, a self supervised network is proposed for binocular disparity matching (SDMNet), which computes dense disparity maps from stereo image pairs without disparity labels In the self supervised training, we match the stereo images densely to approximate the disparity maps and use them to warp the left and right images to estimate the right and left images we build the loss function between estimated images and original images for self supervised training, which adopts perceptual loss to help improve the quality of disparity maps in both detail and structure. Then, we use SDMNet to obtain disparities of hand vein. SDMNet has achieved excellent results on KITTI 2012, KITTI 2015, simulated vein dataset and real vein dataset, outperforming many state of the art supervised matching methods.


Multi Step Inference for Reasoning Over Paragraphs. Complex reasoning over text requires understanding and chaining together free form predicates and logical connectives. Prior work has largely tried to do this either symbolically or with black box transformers. We present a middle ground between these two extremes a compositional model reminiscent of neural module networks that can perform chained logical reasoning. This model first finds relevant sentences in the context and then chains them together using neural modules. Our model gives significant performance improvements (up to 29 relative error reduction when comfibined with a reranker) on ROPES, a recently introduced complex reasoning dataset.


Free stream preserving finite difference schemes for ideal magnetohydrodynamics on curvilinear meshes. In this paper, a high order free stream preserving finite difference weighted essentially non oscillatory (WENO) scheme is developed for the ideal magnetohydrodynamic (MHD) equations on curvilinear meshes. Under the constrained transport framework, magnetic potential evolved by a Hamilton Jacobi (H J) equation is introduced to control the divergence error. In this work, we use the alternative formulation of WENO scheme 10 for the nonlinear hyperbolic conservation law, and design a novel method to solve the magnetic potential. Theoretical derivation and numerical results show that the scheme can preserve free stream solutions of MHD equations, and reduce error more effectively than the standard finite difference WENO schemes for such problems.


SChuBERT Scholarly Document Chunks with BERT encoding boost Citation Count Prediction. Predicting the number of citations of scholarly documents is an upcoming task in scholarly document processing. Besides the intrinsic merit of this information, it also has a wider use as an imperfect proxy for quality which has the advantage of being cheaply available for large volumes of scholarly documents. Previous work has dealt with number of citations prediction with relatively small training data sets, or larger datasets but with short, incomplete input text. In this work we leverage the open access ACL Anthology collection in combination with the Semantic Scholar bibliometric database to create a large corpus of scholarly documents with associated citation information and we propose a new citation prediction model called SChuBERT. In our experiments we compare SChuBERT with several state of the art citation prediction models and show that it outperforms previous methods by a large margin. We also show the merit of using more training data and longer input for number of citations prediction.


Range separated double hybrid density functional theory with coupled cluster and random phase approximations. We construct range separated double hybrid schemes which combine coupled cluster or random phase approximations with a density functional based on a two parameter Coulomb attenuating method like decomposition of the electron electron interaction. We find that the addition of a fraction of short range electron electron interaction in the wave function part of the calculation is globally beneficial for the range separated double hybrid scheme involving a variant of the random phase approximation with exchange terms. Even though the latter scheme is globally as accurate as the corresponding scheme employing only second order M o ller Plesset perturbation theory for atomization energies, reaction barrier heights, and weak intermolecular interactions of small molecules, it is more accurate for the more complicated case of the benzene dimer in the stacked configuration. The present range separated double hybrid scheme employing a random phase approximation thus represents a new member in the family of double hybrids with minimal empiricism which could be useful for general chemical applications.


The Mobility Network of Scientists Analyzing Temporal Correlations in Scientific Careers. The mobility of scientists between different universities and countries is important to foster knowledge exchange. At the same time, the potential mobility is restricted by geographic and institutional constraints, which leads to temporal correlations in the career trajectories of scientists. To quantify this effect, we extract 3.5 million career trajectories of scientists from two large scale bibliographic data sets and analyze them applying a novel method of higher order networks. We study the effect of temporal correlations at three different levels of aggregation, universities, cities and countries. We find strong evidence for such correlations for the top 100 universities, i.e. scientists move likely between specific institutions. These correlations also exist at the level of countries, but cannot be found for cities. Our results allow to draw conclusions about the institutional path dependence of scientific careers and the efficiency of mobility programs.


Essential self adjointness of perturbed quadharmonic operators on Riemannian manifolds with an application to the separation problem. We consider perturbed quadharmonic operators, Delta 4 V , acting on sections of a Hermitian vector bundle over a complete Riemannian manifold, with the potential V satisfying a bound from below by a non positive function depending on the distance from a point. Under a bounded geometry assumption on the Hermitian vector bundle and the underlying Riemannian manifold, we give a sufficient condition for the essential self adjointness of such operators. We then apply this to prove the separation property in L 2 when the perturbed operator acts on functions.


Mapping the Interplanetary Filesystem. The Interplanetary Filesystem (IPFS) is a distributed data storage service frequently used by blockchain applications and for sharing content in a censorship resistant manner. Data is distributed within an open set of peers using a Kademlia based distributed hash table (DHT). In this paper, we study the structure of the resulting overlay network, as it significantly influences the robustness and performance of IPFS. We monitor and systematically crawl IPFS' DHT towards mapping the IPFS overlay network. Our measurements found an average of 44474 nodes at every given time. At least 52.19 of these reside behind a NAT and are not reachable from the outside, suggesting that a large share of the network is operated by private individuals on an as needed basis. Based on our measurements and our analysis of the IPFS code, we conclude that the topology of the IPFS network is, in its current state, closer to an unstructured overlay network than it is to a classical DHT. While such a structure has benefits for robustness and the resistance against Sybil attacks, it leaves room for improvement in terms of performance and query privacy.


Generative Partial Multi View Clustering. Nowadays, with the rapid development of data collection sources and feature extraction methods, multi view data are getting easy to obtain and have received increasing research attention in recent years, among which, multi view clustering (MVC) forms a mainstream research direction and is widely used in data analysis. However, existing MVC methods mainly assume that each sample appears in all the views, without considering the incomplete view case due to data corruption, sensor failure, equipment malfunction, etc. In this study, we design and build a generative partial multi view clustering model, named as GP MVC, to address the incomplete multi view problem by explicitly generating the data of missing views. The main idea of GP MVC lies at two fold. First, multi view encoder networks are trained to learn common low dimensional representations, followed by a clustering layer to capture the consistent cluster structure across multiple views. Second, view specific generative adversarial networks are developed to generate the missing data of one view conditioning on the shared representation given by other views. These two steps could be promoted mutually, where learning common representations facilitates data imputation and the generated data could further explores the view consistency. Moreover, an weighted adaptive fusion scheme is implemented to exploit the complementary information among different views. Experimental results on four benchmark datasets are provided to show the effectiveness of the proposed GP MVC over the state of the art methods.


Trend to equilibrium for granular media equations under non convex potential and application to log Coulomb gases. We investigate the long time asymptotics of solutions of granular media equations with non convex potentials, by establishing a HWI inequality. Our result extends the work of Carrillo et al.. Moreover, we consider unidimensional log Coulomb gases. We prove existence and uniqueness of an equilibrium for the entropy and show exponential convergence of symmetric solutions towards this equilibrium under non convex polynomial potential, with respect to the Wasserstein metric. Finally, we study the case of the confining and the non confining quartic potentials, and answer conjectures made by Li et al. and Biane et al..


Impact of aperture, depth, and acoustic clutter on performance of Coherent Multi Transducer Ultrasound imaging. Transducers with larger aperture size are desirable in ultrasound imaging to improve resolution and image quality. A coherent multi transducer ultrasound imaging system (CoMTUS) enables an extended effective aperture through coherent combination of multiple transducers. In this study, the discontinuous extended aperture created by CoMTUS and its performance for deep imaging and through layered media are investigated by both simulations and experiments. Typical image quality metrics resolution, contrast and contrast to noise ratio are evaluated and compared with a standard single probe imaging system. Results suggest that, the image performance of CoMTUS depends on the relative spatial location of the arrays. The resulting effective aperture significantly improves resolution, while the separation between the arrays may degrade contrast. For a limited gap in the effective aperture (less than few centimetres), CoMTUS provides benefits to image quality compared to standard single probe imaging system. Overall, CoMTUS shows higher sensitivity and reduced loss of resolution with imaging depth. In general, CoMTUS imaging performance was unaffected when imaging through a layered medium with different speed of sound values and resolution improved up to 80 at large imaging depths.


Avoidable paths in graphs. We prove a recent conjecture of Beisegel et al. that for every positive integer k, every graph containing an induced P k also contains an avoidable P k. Avoidability generalises the notion of simpliciality best known in the context of chordal graphs. The conjecture was only established for k in 1,2 (Ohtsuki et al. 1976, and Beisegel et al. 2019, respectively). Our result also implies a result of Chv 'atal et al. 2002, which assumed cycle restrictions. We provide a constructive and elementary proof, relying on a single trick regarding the induction hypothesis. In the line of previous works, we discuss conditions for multiple avoidable paths to exist.


Deep ahead of threat virtual patching. Many applications have security vulnerabilities that can be exploited. It is practically impossible to find all of them due to the NP complete nature of the testing problem. Security solutions provide defenses against these attacks through continuous application testing, fast patching of vulnerabilities, automatic deployment of patches, and virtual patching detection techniques deployed in network and endpoint security tools. These techniques are limited by the need to find vulnerabilities before the black hats. We propose an innovative technique to virtually patch vulnerabilities before they are found. We leverage testing techniques for supervised learning data generation, and show how artificial intelligence techniques can use this data to create predictive deep neural network models that read an application's input and predict in real time whether it is a potential malicious input. We set up an ahead of threat experiment in which we generated data on old versions of an application, and then evaluated the predictive model accuracy on vulnerabilities found years later. Our experiments show ahead of threat detection on LibXML2 and LibTIFF vulnerabilities with 91.3 and 93.7 accuracy, respectively. We expect to continue work on this field of research and provide ahead of threat virtual patching for more libraries. Success in this research can change the current state of endless racing after application vulnerabilities and put the defenders one step ahead of the attackers


A Model Counter's Guide to Probabilistic Systems. In this paper, we systematize the modeling of probabilistic systems for the purpose of analyzing them with model counting techniques. Starting from unbiased coin flips, we show how to model biased coins, correlated coins, and distributions over finite sets. From there, we continue with modeling sequential systems, such as Markov chains, and revisit the relationship between weighted and unweighted model counting. Thereby, this work provides a conceptual framework for deriving SAT encodings for probabilistic inference.


On the Numerical Analysis and Visualisation of Implicit Ordinary Differential Equations. We discuss how the geometric theory of differential equations can be used for the numerical integration and visualisation of implicit ordinary differential equations, in particular around singularities of the equation. The Vessiot theory automatically transforms an implicit differential equation into a vector field distribution on a manifold and thus reduces its analysis to standard problems in dynamical systems theory like the integration of a vector field and the determination of invariant manifolds. For the visualisation of low dimensional situations we adapt the streamlines algorithm of Jobard and Lefer to 2.5 and 3 dimensions. A concrete implementation in Matlab is discussed and some concrete examples are presented.


Incremental Refinements and Multiple Descriptions with Feedback. It is well known that independent (separate) encoding of K correlated sources may incur some rate loss compared to joint encoding, even if the decoding is done jointly. This loss is particularly evident in the multiple descriptions problem, where the sources are repetitions of the same source, but each description must be individually good. We observe that under mild conditions about the source and distortion measure, the rate ratio Rindependent(K) Rjoint goes to one in the limit of small rate high distortion. Moreover, we consider the excess rate with respect to the rate distortion function, Rindependent(K, M) R(D), in M rounds of K independent encodings with a final distortion level D. We provide two examples a Gaussian source with mean squared error and an exponential source with one sided error for which the excess rate vanishes in the limit as the number of rounds M goes to infinity, for any fixed D and K. This result has an interesting interpretation for a multi round variant of the multiple descriptions problem, where after each round the encoder gets a (block) feedback regarding which of the descriptions arrived In the limit as the number of rounds M goes to infinity (i.e., many incremental rounds), the total rate of received descriptions approaches the rate distortion function. We provide theoretical and experimental evidence showing that this phenomenon is in fact more general than in the two examples above.


Comparison of Monge Amp ere capacities. Let (X, omega) be a compact K "ahler manifold. We prove that all Monge Amp ere capacities are comparable. Using this we give an alternative direct proof of the integration by parts formula for non pluripolar products recently proved by M. Xia.


A kernel independent sum of Gaussians method by de la Vall 'ee Poussin sums. Approximation of interacting kernels by sum of Gaussians (SOG) is frequently required in many applications of scientific and engineering computing in order to construct efficient algorithms for kernel summation or convolution problems. In this paper, we propose a kernel independent SOG method by introducing the de la Vall 'ee Poussin sum and Chebyshev polynomials. The SOG works for general interacting kernels and the lower bound of Gaussian bandwidths is tunable and thus the Gaussians can be easily summed by fast Gaussian algorithms. The number of Gaussians can be further reduced via the model reduction based on the balanced truncation based on the square root method. Numerical results on the accuracy and model reduction efficiency show attractive performance of the proposed method.


Regularity of weak supersolutions to elliptic and parabolic equations lower semicontinuity and pointwise behavior. We demonstrate a measure theoretical approach to the local regularity of weak supersolutions to elliptic and parabolic equations in divergence form. In the first part, we show that weak supersolutions become lower semicontinuous after redefinition on a set of measure zero. The proof relies on a general principle, i.e. the De Giorgi type lemma, which offers a unified approach for a wide class of elliptic and parabolic equations, including an anisotropic elliptic equation, the parabolic p Laplace equation, and the porous medium equation. In the second part, we shall show that for parabolic problems the lower semicontinuous representative at an instant can be recovered pointwise from the "essliminf" of past times. We also show that it can be recovered by the limit of certain integral average of past times. The proof hinges on the expansion of positivity for weak supersolutions. Our results are structural properties of partial differential equations, independent of any kind of comparison principle.


Revisited bilinear Schr " o dinger estimates with applications to generalized Boussinesq equations. In this paper, our goal is to improve the local well posedness theory for certain generalized Boussinesq equations by revisiting bilinear estimates related to the Schr " o dinger equation. Moreover, we propose a novel, automated procedure to handle the summation argument for these bounds.


On the complexity of classes of uncountable structures trees on aleph 1 . We analyse the complexity of the class of (special) Aronszajn, Suslin and Kurepa trees in the projective hierarchy of the higher Baire space omega 1 omega 1 . First, we will show that none of these classes have the Baire property (unless they are empty). Moreover, under (V L) , (a) the class of Aronszajn and Suslin trees is Pi 1 1 complete, (b) the class of special Aronszajn trees is Sigma 1 1 complete, and (c) the class of Kurepa trees is Pi 1 2 complete. We achieve these results by finding nicely definable reductions that map subsets X of omega 1 to trees T X so that T X is in a given tree class mathcal T if and only if X is stationary non stationary (depending on the class mathcal T ). Finally, we present models of CH where these classes have lower projective complexity.


Globally simple Heffter arrays H(n k) when k equiv 0,3 mod 4 . Square Heffter arrays are n times n arrays such that each row and each column contains k filled cells, each row and column sum is divisible by 2nk 1 and either x or x appears in the array for each integer 1 leq x leq nk . Archdeacon noted that a Heffter array, satisfying two additional conditions, yields a face 2 colourable embedding of the complete graph K 2nk 1 on an orientable surface, where for each colour, the faces give a k cycle system. Moreover, a cyclic permutation on the vertices acts as an automorphism of the embedding. These necessary conditions pertain to cyclic orderings of the entries in each row and each column of the Heffter array and are (1) for each row and each column the sequential partial sums determined by the cyclic ordering must be distinct modulo 2nk 1 (2) the composition of the cyclic orderings of the rows and columns is equivalent to a single cycle permutation on the entries in the array. We construct Heffter arrays that satisfy condition (1) whenever (a) k equiv 0 mod 4 or (b) n equiv 1 mod 4 and k equiv 3 mod 4 or (c) n equiv 0 mod 4 , k equiv3 mod 4 and n gg k . As corollaries to the above we obtain pairs of orthogonal k cycle decompositions of K 2nk 1 .


New Techniques and Fine Grained Hardness for Dynamic Near Additive Spanners. Maintaining and updating shortest paths information in a graph is a fundamental problem with many applications. As computations on dense graphs can be prohibitively expensive, and it is preferable to perform the computations on a sparse skeleton of the given graph that roughly preserves the shortest paths information. Spanners and emulators serve this purpose. This paper develops fast dynamic algorithms for sparse spanner and emulator maintenance and provides evidence from fine grained complexity that these algorithms are tight. Under the popular OMv conjecture, we show that there can be no decremental or incremental algorithm that maintains an n 1 o(1) edge (purely additive) n delta emulator for any delta 1 2 with arbitrary polynomial preprocessing time and total update time m 1 o(1) . Also, under the Combinatorial k Clique hypothesis, any fully dynamic combinatorial algorithm that maintains an n 1 o(1) edge (1 epsilon,n o(1) ) spanner or emulator must either have preprocessing time mn 1 o(1) or amortized update time m 1 o(1) . Both of our conditional lower bounds are tight. As the above fully dynamic lower bound only applies to combinatorial algorithms, we also develop an algebraic spanner algorithm that improves over the m 1 o(1) update time for dense graphs. For any constant epsilon in (0,1 , there is a fully dynamic algorithm with worst case update time O(n 1.529 ) that whp maintains an n 1 o(1) edge (1 epsilon,n o(1) ) spanner. Our new algebraic techniques and spanner algorithms allow us to also obtain (1) a new fully dynamic algorithm for All Pairs Shortest Paths (APSP) with update and path query time O(n 1.9 ) (2) a fully dynamic (1 epsilon) approximate APSP algorithm with update time O(n 1.529 ) (3) a fully dynamic algorithm for near 2 approximate Steiner tree maintenance.


Maximum of the integer valued Gaussian free field. We investigate the order of the maximum of the integer valued Gaussian free field in two dimensions, and show that it grows logarithmically with the size of the box. Our treatment follows closely that of a recent paper by Kharash and Peled on the Fr " o hlich Spencer proof of the Berezinskii Kosterlitz Thouless transition.


Automatic Detection of Satire in Bangla Documents A CNN Approach Based on Hybrid Feature Extraction Model. Widespread of satirical news in online communities is an ongoing trend. The nature of satires is so inherently ambiguous that sometimes it's too hard even for humans to understand whether it's actually satire or not. So, research interest has grown in this field. The purpose of this research is to detect Bangla satirical news spread in online news portals as well as social media. In this paper, we propose a hybrid technique for extracting features from text documents combining Word2Vec and TF IDF. Using our proposed feature extraction technique, with standard CNN architecture we could detect whether a Bangla text document is satire or not with an accuracy of more than 96 .


On simple modules of cyclotomic quiver Hecke algebras of type A. We give a classification of the graded simple modules of cyclotomic quiver Hecke algebras of type A using the diagram calculus of the diagrammatic Cherednik algebra. We also obtain a non trivial lower bound for the dimension of the simple modules which is independent of the field characteristic.


On Lieb Thirring inequalities for one dimensional non self adjoint Jacobi and Schr " o dinger operators. We study to what extent Lieb Thirring inequalities are extendable from self adjoint to general (possibly non self adjoint) Jacobi and Schr " o dinger operators. Namely, we prove the conjecture of Hansmann and Katriel from Complex Anal. Oper. Theory 5, No. 1 (2011), 197 218 and answer another open question raised therein. The results are obtained by means of asymptotic analysis of eigenvalues of discrete Schr " o dinger operators with rectangular barrier potential and complex coupling. Applying the ideas in the continuous setting, we also solve a similar open problem for one dimensional Schr " o dinger operators with complex valued potentials published by Demuth, Hansmann, and Katriel in Integral Equations Operator Theory 75, No. 1 (2013), 1 5 .


Miura ori tube metamaterial with tunable dynamic property. Structures and or materials with engineered functionality, capable of achieving targeted mechanical responses reacting to changes in external excitation, have various potential engineering applications, e.g. aerospace, oceanographic engineering, soft robot, and several others. Yet tunable mechanical performance is normally realized through carefully designing the architecture of structures, which is usually porous, leading to the complexity of the fabrication of the structures even using the recently emerged 3D printing technique. In this study we show that origami technique can provide an alternative solution to achieving the aim by carefully stacking the classical Miura sheets into the Miura ori tube metamaterial and tuning the geometric parameters of the origami metamaterial. By combining numerical and experimental methods, we have demonstrated that an extremely broad range of natural frequency and dynamic response of the metamaterial can be achieved. The proposed structure can be easily fabricated from a single thin sheet made of one material and simultaneously owns better mechanical properties than the Miura sheet.


Non invariant deformations of left invariant complex structures on compact Lie groups. We consider deformations of left invariant complex structures on simply connected semisimple compact Lie groups which are a priori non invariant. Computing their cohomologies, we show that they are not actually biholomorphic to left invariant (right invariant) complex structures.


Modeling Engagement in Long Term, In Home Socially Assistive Robot Interventions for Children with Autism Spectrum Disorders. Socially assistive robotics (SAR) has great potential to provide accessible, affordable, and personalized therapeutic interventions for children with autism spectrum disorders (ASD). However, human robot interaction (HRI) methods are still limited in their ability to autonomously recognize and respond to behavioral cues, especially in atypical users and everyday settings. This work applies supervised machine learning algorithms to model user engagement in the context of long term, in home SAR interventions for children with ASD. Specifically, we present two types of engagement models for each user (i) generalized models trained on data from different users and (ii) individualized models trained on an early subset of the user's data. The models achieved approximately 90 accuracy (AUROC) for post hoc binary classification of engagement, despite the high variance in data observed across users, sessions, and engagement states. Moreover, temporal patterns in model predictions could be used to reliably initiate re engagement actions at appropriate times. These results validate the feasibility and challenges of recognition and response to user disengagement in long term, real world HRI settings. The contributions of this work also inform the design of engaging and personalized HRI, especially for the ASD community.


Fine Grained Dynamic Head for Object Detection. The Feature Pyramid Network (FPN) presents a remarkable approach to alleviate the scale variance in object representation by performing instance level assignments. Nevertheless, this strategy ignores the distinct characteristics of different sub regions in an instance. To this end, we propose a fine grained dynamic head to conditionally select a pixel level combination of FPN features from different scales for each instance, which further releases the ability of multi scale feature representation. Moreover, we design a spatial gate with the new activation function to reduce computational complexity dramatically through spatially sparse convolutions. Extensive experiments demonstrate the effectiveness and efficiency of the proposed method on several state of the art detection benchmarks. Code is available at https github.com StevenGrove DynamicHead.


Occupation densities of Ensembles of Branching Random Walks. We study the limiting occupation density process for a large number of critical and driftless branching random walks. We show that the rescaled occupation densities of lfloor sN rfloor branching random walks, viewed as a function valued, increasing process g s N s ge 0 , converges weakly to a pure jump process in the Skorohod space mathbb D( 0, infty), mathcal C 0 ( mathbb R)) , as N to infty . Moreover, the jumps of the limiting process consist of i.i.d. copies of an Integrated super Brownian Excursion (ISE) density, rescaled and weighted by the jump sizes in a real valued stable 1 2 subordinator.


A thermodynamically consistent pseudo potential lattice Boltzmann model for multi component, multiphase, partially miscible mixtures. Current multi component, multiphase pseudo potential lattice Boltzmann models have thermodynamic inconsistencies that prevent them to correctly predict the thermodynamic phase behavior of partially miscible multi component mixtures, such as hydrocarbon mixtures. This paper identifies these inconsistencies and attempts to design a thermodynamically consistent multi component, multiphase pseudo potential lattice Boltzmann model that allows mass transfer across the phase interfaces and is capable to predict the phase behavior of typically partially miscible hydrocarbon mixtures. The designed model defines the total interaction force for the entire phase and split the force into individual components. Through a properly derived force split factor associated with the volatility of each component, the model can achieve precise thermodynamic consistency in multi component hydrocarbon mixtures, which is described by the iso fugacity rule.


TEASER Fast and Certifiable Point Cloud Registration. We propose the first fast and certifiable algorithm for the registration of two sets of 3D points in the presence of large amounts of outlier correspondences. We first reformulate the registration problem using a Truncated Least Squares (TLS) cost that is insensitive to a large fraction of spurious correspondences. Then, we provide a general graph theoretic framework to decouple scale, rotation, and translation estimation, which allows solving in cascade for the three transformations. Despite the fact that each subproblem is still non convex and combinatorial in nature, we show that (i) TLS scale and (component wise) translation estimation can be solved in polynomial time via adaptive voting, (ii) TLS rotation estimation can be relaxed to a semidefinite program (SDP) and the relaxation is tight, even in the presence of extreme outlier rates, and (iii) the graph theoretic framework allows drastic pruning of outliers by finding the maximum clique. We name the resulting algorithm TEASER (Truncated least squares Estimation And SEmidefinite Relaxation). While solving large SDP relaxations is typically slow, we develop a second fast and certifiable algorithm, named TEASER , that uses graduated non convexity to solve the rotation subproblem and leverages Douglas Rachford Splitting to efficiently certify global optimality. For both algorithms, we provide theoretical bounds on the estimation errors, which are the first of their kind for robust registration problems. Moreover, we test their performance on standard, object detection, and the 3DMatch benchmarks, and show that (i) both algorithms dominate the state of the art and are robust to more than 99 outliers, (ii) TEASER can run in milliseconds, and (iii) TEASER is so robust it can also solve problems without correspondences, where it largely outperforms ICP and it is more accurate than Go ICP while being orders of magnitude faster.


Summation arithmetic functions with bounded terms, having a limit normal distribution law. The paper considers the properties of pseudo stationarity in a broad sense and pseudo strong mixing for sequences of random variables corresponding to arithmetic functions. Assertions on this topic have been proven. The implementation of these properties for known arithmetic functions has been verified. The article proves a statement about sufficient conditions under which a summation arithmetic function with bounded terms has a limit normal distribution law. The fulfillment of the specified sufficient conditions for known arithmetic functions is considered.


The poset of copies for automorphism groups of countable relational structures. Let mathrm G be a subgroup of the symmetric group mathfrak S(U) of all permutations of a countable set U . Let overline mathrm G be the topological closure of mathrm G in the function topology on U U . We initiate the study of the poset overline mathrm G U f U mid f in overline mathrm G of images of the functions in overline mathrm G , being ordered under inclusion. This set overline mathrm G U of subsets of the set U will be called the emph poset of copies for the group mathrm G . A denomination being justified by the fact that for every subgroup mathrm G of the symmetric group mathfrak S(U) there exists a homogeneous relational structure R on U such that overline G is the set of embeddings of the homogeneous structure R into itself and overline mathrm G U is the set of copies of R in R and that the set of bijections overline G cap mathfrak S(U) of U to U forms the group of automorphisms of mathrm R .


Locating the source of interacting signal in complex networks. We investigate the problem of locating the source of a self interacting signal spreading in a complex networks. We use a well known rumour model as an example of the process with self interaction. According to this model based on the SIR epidemic dynamics, the infected nodes may interact and discourage each other from gossiping with probability alpha . We compare three algorithms of source localization Limited Pinto Thiran Vettarli (LPTV), Gradient Maximum Likelihood (GMLA) and one based on Pearson correlation between time and distance. The results of numerical simulations show that additional interactions between infected nodes decrease the quality of LPTV and Pearson. GMLA is the most resistant to harmful effects of the self interactions, which is especially visible for medium and high level of stochasticity of the process, when spreading rate is below 0.5. The reason for this may be the fact that GMLA uses only the nearest observers, which are much less likely affected by the interactions between infected nodes, because these contacts become important as the epidemics develops and the number of infected agents increases.


Team Optimal Control of Coupled Subsystems with Mean Field Sharing. We investigate team optimal control of stochastic subsystems that are weakly coupled in dynamics (through the mean field of the system) and are arbitrary coupled in the cost. The controller of each subsystem observes its local state and the mean field of the state of all subsystems. The system has a non classical information structure. Exploiting the symmetry of the problem, we identify an information state and use that to obtain a dynamic programming decomposition. This dynamic program determines a globally optimal strategy for all controllers. Our solution approach works for arbitrary number of controllers and generalizes to the setup when the mean field is observed with noise. The size of the information state is time invariant thus, the results generalize to the infinite horizon control setups as well. In addition, when the mean field is observed without noise, the size of the corresponding information state increases polynomially (rather than exponentially) with the number of controllers which allows us to solve problems with moderate number of controllers. We illustrate our approach by an example motivated by smart grids that consists of 100 coupled subsystems.


A refinement of Gorenstein flat dimension via the flat cotorsion theory. We introduce a refinement of the Gorenstein flat dimension for complexes over an associative ring the Gorenstein flat cotorsion dimension and prove that it, unlike the Gorenstein flat dimension, behaves as one expects of a homological dimension without extra assumptions on the ring. Crucially, we show that it coincides with the Gorenstein flat dimension for complexes where the latter is finite, and for complexes over right coherent rings the setting where the Gorenstein flat dimension is known to behave as expected.


Query guided End to End Person Search. Person search has recently gained attention as the novel task of finding a person, provided as a cropped sample, from a gallery of non cropped images, whereby several other people are also visible. We believe that i. person detection and re identification should be pursued in a joint optimization framework and that ii. the person search should leverage the query image extensively (e.g. emphasizing unique query patterns). However, so far, no prior art realizes this. We introduce a novel query guided end to end person search network (QEEPS) to address both aspects. We leverage a most recent joint detector and re identification work, OIM 37 . We extend this with i. a query guided Siamese squeeze and excitation network (QSSE Net) that uses global context from both the query and gallery images, ii. a query guided region proposal network (QRPN) to produce query relevant proposals, and iii. a query guided similarity subnetwork (QSimNet), to learn a query guided reidentification score. QEEPS is the first end to end query guided detection and re id network. On both the most recent CUHK SYSU 37 and PRW 46 datasets, we outperform the previous state of the art by a large margin.


Scalable parallel algorithm for solving non stationary systems of linear inequalities. In this paper, a scalable iterative projection type algorithm for solving non stationary systems of linear inequalities is considered. A non stationary system is understood as a large scale system of inequalities in which coefficients and constant terms can change during the calculation process. The proposed parallel algorithm uses the concept of pseudo projection which generalizes the notion of orthogonal projection. The parallel pseudo projection algorithm is implemented using the parallel BSF skeleton. An analytical estimation of the algorithm scalability boundary is obtained on the base of the BSF cost metric. The large scale computational experiments were performed on a cluster computing system. The obtained results confirm the efficiency of the proposed approach.


The Proceedings of First Work in Progress Session of The CSI International Symposium on Real Time and Embedded Systems and Technologies. The present volume contains the proceedings of RTEST WiP 2018, chaired by Marco Caccamo, University of Illinois at Urbana Champaign. This event has been organized by the School of Electrical and Computer Engineering at the University of Tehran, in conjunction with the Department of Computer Engineering at Sharif University of Technology, Tehran, Iran. The topics of interest in RTEST WiP span over all theoretical and application oriented aspects, reporting design, analysis, implementation, evaluation, and empirical results, of real time and embedded systems, internet of things, and cyber physical systems. The program committee of RTEST 2018 consists of 54 top researchers in the mentioned fields from top universities, industries, and research centers around the world. RTEST 2018 has received a total of 41 submissions, out of which we have accepted 14 regular papers and 4 work in progress papers. Each submission has been reviewed by 3 to 5 independent referees, for its quality, originality, contribution, clarity of presentation, and relevance to the symposium topics.


A Hybrid Achromatic Metalens. Metamaterials and metasurfaces are widely used to manipulate electromagnetic waves over a broad range of wavelengths. Several recent efforts have focused on metalenses, ultra thin optical elements that focus light using subwavelength structures. Compared to their refractive counterparts, metalenses offer reduced size and weight, improved manufacturability, and new functionality such as polarization control. However, metalenses that correct chromatic aberration also suffer from markedly reduced focusing efficiency. Here we introduce a Hybrid Achromatic Metalens (HAML), that overcomes this trade off and offers improved focusing efficiency over a broad wavelength range from 1000 1800 nm. Fabricated HAMLs demonstrated diffraction limited performance for numerical apertures (NA) of 0.27, 0.11, and 0.06 with average focusing efficiencies 60 and maximum efficiencies 80 . HAMLs can be designed by combining recursive ray tracing and simulated phase libraries rather than computational intensive global search algorithms. Moreover, HAMLs can be fabricated in low refractive index materials using multi photon lithography for customization or using molding for mass production.


Time restricted sensitivity and entropy. In this paper, we consider measure theoretical restricted sensitivity and topological restricted sensitivities by restricting the first sensitive time. For a given topological dynamical system, we define measure theoretical restricted asymptotic rate with respect to sensitivity, and obtain that it equal to the reciprocal of the Brin Katok local entropy for almost every point. For topological version we have similar definitions and conclusions.


Optical theorems and physical bounds on absorption in lossy media. Two different versions of an optical theorem for a scattering body embedded inside a lossy background medium are derived in this paper. The corresponding fundamental upper bounds on absorption are then obtained in closed form by elementary optimization techniques. The first version is formulated in terms of polarization currents (or equivalent currents) inside the scatterer and generalizes previous results given for a lossless medium. The corresponding bound is referred to here as a variational bound and is valid for an arbitrary geometry with a given material property. The second version is formulated in terms of the T matrix parameters of an arbitrary linear scatterer circumscribed by a spherical volume and gives a new fundamental upper bound on the total absorption of an inclusion with an arbitrary material property (including general bianisotropic materials). The two bounds are fundamentally different as they are based on different assumptions regarding the structure and the material property. Numerical examples including homogeneous and layered (core shell) spheres are given to demonstrate that the two bounds provide complimentary information in a given scattering problem.


Novel therapeutic targets in chronic myeloid leukaemia through a discrete time discrete Markov chain model of BCR ABL1 interactions. Chronic Myeloid Leukaemia (CML) is a blood derived proliferative disorder, which is highly associated to a translocation of chromosomes 9 and 22 or the creation of Philadelphia chromosome Ph( ) cases, inducing the synthesis of a chimeric fusion protein, namely BCR ABL1 (Breakpoint Cluster Region Abelson 1 chimeric protein), which is known for driving the pathophysiology of the disease, however variants of CML are also recognized as CML Ph( ), these nonetheless account for a small percentage of the overall CML patients posing thus the question whether BCR ABL1 fusion protein is required for the whole of the pathophysiology of CML. Hereof, through a stochastic description, a discrete time discrete Markov chain depicts the various protein protein interactions of BCR ABL1 to better understand signalling pathways and time dependent evolution of these pathways, as well as to provide prospective therapeutic protein targets to improve both the specificity of the treatment and the life expectancy of the patients.


Long range Rydberg molecule Rb 2 Two electron textit R matrix calculations at intermediate internuclear distances. The adiabatic potential energy curves of Rb 2 in the long range Rydberg electronic states are calculated using the two electron textit R matrix method M. Tarana, R. v C ur ' i k, Phys. Rev. A textbf 93 , 012515 (2016) for the intermediate internuclear separations between 35 a.u. and 200 a.u. The results are compared with the zero range models to find a region of the internuclear distances where the Fermi's pseudopotential approach provides accurate energies. A finite range potential model of the atomic perturber is used to calculate the wave functions of the Rydberg electron and their features specific for the studied range of internuclear distances are identified.


On the Capacity of Channels with Deletions and States. We consider the class of channels formed from the concatenation of a deletion channel and a finite state channel. For this class of channels, we show that the operational capacity is equal to the stationary capacity, which can be approached by a sequence of Markov processes with increasing Markovian orders. As a by product, we show that the polar coding scheme constructed by Tal, Pfister, Fazeli and Vardy achieves the capacity of the deletion channel.


Learning Robust 3D Face Reconstruction and Discriminative Identity Representation. 3D face reconstruction from a single 2D image is a very important topic in computer vision. However, the current reconstruction methods are usually non sensitive to face identities and over sensitive to facial poses, which may result in similar 3D geometries for faces of different identities, or obtain different shapes for the same identity with different poses. When such methods are applied practically, their 3D estimates are either changeable for different photos of the same subject or over regularized and generic to distinguish face identities. In this paper, we propose a robust solution to solve this problem by carefully designing a novel Siamese Convolutional Neural Network (SCNN). Specifically, regarding the 3D Morphable face Model (3DMM) parameters of the same individual as the same class, we employ the contrastive loss to enlarge the inter class distance and meanwhile reduce the intra class distance for the output 3DMM parameters. We also propose an identity loss to preserve the identity information for the same individual in the feature space. Training with these two losses, our SCNN could learn representations that are more discriminative for face identity and generalizable for pose variants. Experiments on the challenging database 300W LP and AFLW2000 3D have shown the effectiveness of our method by comparing with state of the arts.


Definitions and Semantic Simulations Based on Object Oriented Analysis and Modeling. We have proposed going beyond traditional ontologies to use rich semantics implemented in programming languages for modeling. In this paper, we discuss the application of executable semantic models to two examples, first a structured definition of a waterfall and second the cardiopulmonary system. We examine the components of these models and the way those components interact. Ultimately, such models should provide the basis for direct representation.


Inequity exemplified by underrepresented minority students receiving lower grades and having higher rates of attrition across STEM disciplines. Underrepresented minority (URM) students are subjected to historically rooted inequities when pursuing an education, especially in STEM disciplines with little diversity. In order to make STEM education equitable and inclusive, evidence for how students from different racial ethnic demographics are faring is necessary. We use 10 years of institutional data at a large public university to investigate trends in the majors that Asian, URM, and White students declare, drop after declaring, and earn degrees in as well as the GPA of the students who drop or earn a degree. We find that higher percentages of the URM students drop most majors compared to other students and these trends are particularly pronounced in physics and economics. Moreover, we find alarming GPA trends in that the URM students consistently earn lower grades than their Asian and White peers. Furthermore, in some STEM disciplines, the URM students who earn a degree are earning the same grades as the Asian and White students who dropped the major. This troubling trend may signify lack of sufficient support, mentoring, and guidance to ensure excellence of the URM students who are already severely disadvantaged. These quantitative findings call for making learning environments equitable and inclusive so that many URM students who come to college with severe disadvantages are appropriately supported and can excel similar to other students.


Persian Signature Verification using Fully Convolutional Networks. Fully convolutional networks (FCNs) have been recently used for feature extraction and classification in image and speech recognition, where their inputs have been raw signal or other complicated features. Persian signature verification is done using conventional convolutional neural networks (CNNs). In this paper, we propose to use FCN for learning a robust feature extraction from the raw signature images. FCN can be considered as a variant of CNN where its fully connected layers are replaced with a global pooling layer. In the proposed manner, FCN inputs are raw signature images and convolution filter size is fixed. Recognition accuracy on UTSig database, shows that FCN with a global average pooling outperforms CNN.


Computing Shortest Paths and Diameter in the Hybrid Network Model. The mathsf HYBRID model, introduced in Augustine et al., SODA '20 , provides a theoretical foundation for networks that allow multiple communication modes. The model follows the principles of synchronous message passing, whereas nodes are allowed to use textit two fundamentally different communication modes. First, a local mode where nodes may exchange arbitrary information per round over edges of a local communication graph G (akin to the mathsf LOCAL model). Second, a global mode where every node may exchange O( log n) messages of size O( log n) bits per round with arbitrary nodes in the network. The mathsf HYBRID model intends to reflect the conditions of many real hybrid networks, where high bandwidth but inherently local communication is combined with highly flexible global communication with restricted bandwidth. We continue to explore the power and limitations of the mathsf HYBRID model by investigating the complexity of computing shortest paths and diameter of the local communication graph G . We improve on the known upper bound for the exact all pairs shortest paths problem (APSP) from Augustine et al., SODA '20 and provide algorithms to approximate solutions for the k source shortest paths problem ( k SSSP). We demonstrate that our results for APSP and k SSP are almost tight (up to poly logarithmic factors). Furthermore, we give an improved algorithm for the exact single source shortest paths problem for graphs with large diameter. For the problem of approximating the diameter of the local communication network we give the first non trivial upper bound. This upper bound is complemented by a lower bound for the exact diameter problem.


Identifying Vulnerabilities of Industrial Control Systems using Evolutionary Multiobjective Optimisation. In this paper we propose a novel methodology to assist in identifying vulnerabilities in a real world complex heterogeneous industrial control systems (ICS) using two evolutionary multiobjective optimisation (EMO) algorithms, NSGA II and SPEA2. Our approach is evaluated on a well known benchmark chemical plant simulator, the Tennessee Eastman (TE) process model. We identified vulnerabilities in individual components of the TE model and then made use of these to generate combinatorial attacks to damage the safety of the system, and to cause economic loss. Results were compared against random attacks, and the performance of the EMO algorithms were evaluated using hypervolume, spread and inverted generational distance (IGD) metrics. A defence against these attacks in the form of a novel intrusion detection system was developed, using a number of machine learning algorithms. Designed approach was further tested against the developed detection methods. Results demonstrate that EMO algorithms are a promising tool in the identification of the most vulnerable components of ICS, and weaknesses of any existing detection systems in place to protect the system. The proposed approach can be used by control and security engineers to design security aware control, and test the effectiveness of security mechanisms, both during design, and later during system operation.