A (Slightly) Improved Approximation Algorithm for Metric TSP. For some epsilon 10 36 we give a 3 2 epsilon approximation algorithm for metric TSP.


On the Endomorphism Semigroups of Extra special p groups and Automorphism Orbits. For an odd prime p and a positive integer n , it is well known that there are two types of extra special p groups of order p 2n 1 , first one is the Heisenberg group which has exponent p and the second one is of exponent p 2 . In this article, a new way of representing the extra special p group of exponent p 2 is given. These representations facilitate an explicit way of finding formulae for any endomorphism and any automorphism of an extra special p group G for both the types. Based on these formulae, the endomorphism semigroup End(G) and the automorphism group Aut(G) are described. The endomorphism semigroup image of any element in G is found and the orbits under the action of the automorphism group Aut(G) are determined. As a consequence it is deduced that, under the notion of degeneration of elements in G , the endomorphism semigroup End(G) induces a partial order on the automorphism orbits when G is the Heisenberg group and does not induce when G is the extra special p group of exponent p 2 . Finally we prove that the cardinality of isotropic subspaces of any fixed dimension in a non degenerate symplectic space is a polynomial in p with non negative integer coefficients. Using this fact we compute the cardinality of End(G) .


Almost Sure Central Limit Theorems in Stochastic Geometry. We prove an almost sure central limit theorem on the Poisson space, which is perfectly tailored for stabilizing functionals emerging in stochastic geometry. As a consequence, we provide almost sure central limit theorems for (i) the total edge length of the k nearest neighbors random graph, (ii) the clique count in random geometric graphs, (iii) the volume of the set approximation via the Poisson Voronoi tessellation.


Lattice PUF A Strong Physical Unclonable Function Provably Secure against Machine Learning Attacks. We propose a strong physical unclonable function (PUF) provably secure against machine learning (ML) attacks with both classical and quantum computers. Its security is derived from cryptographic hardness of learning decryption functions of public key cryptosystems. Our design compactly realizes the decryption function of the learning with errors (LWE) cryptosystem. Due to the fundamental connection of LWE to lattice problems, we call the construction the lattice PUF. Lattice PUF is constructed using a physically obfuscated key (POK), an LWE decryption function block, and a linear feedback shift register (LFSR) as a pseudo random number generator. The POK provides the secret key of the LWE decryption function its stability is ensured by a fuzzy extractor (FE). To reduce the challenge size, we exploit distributional relaxations of space efficient LWEs. That allows only a small challenge seed to be transmitted with the full length challenge generated by the LFSR, resulting in a 100X reduction of communication cost. To prevent an active challenge manipulation attack, a self incrementing counter is embedded into the challenge seed. We prototyped the lattice PUF with 2 136 challenge response pairs (CRPs) on a Spartan 6 FPGA, which required 45 slices for the PUF logic proper and 233 slices for the FE. Simulation based evaluation shows the mean (std) of uniformity to be 49.98 (1.58 ), of uniqueness to be 50.00 (1.58 ), and of reliability to be 1.26 (2.88 ). The LWE concrete hardness estimator guarantees that a successful ML attack of the lattice PUF will require the infeasible 2 128 CPU operations. Several classes of empirical ML attacks, including support vector machine, logistic regression, and deep neural networks, are used in all attacks, the prediction error remains above 49.76 after 1 million training CRPs.


Bordism und projective space bundles. We prove that total spaces of CP 2 bundles generate the oriented cobordism ring Omega . Combined with the surgery lemma this yields a somewhat different proof of Gromov's and Lawson's theorem that all simply connected non spin manifolds of dimension geq 5 carry a metric of positive scalar curvature.


Uncertainty Weighted Causal Graphs. Causality has traditionally been a scientific way to generate knowledge by relating causes to effects. From an imaginery point of view, causal graphs are a helpful tool for representing and infering new causal information. In previous works, we have generated automatically causal graphs associated to a given concept by analyzing sets of documents and extracting and representing the found causal information in that visual way. The retrieved information shows that causality is frequently imperfect rather than exact, feature gathered by the graph. In this work we will attempt to go a step further modelling the uncertainty in the graph through probabilistic improving the management of the imprecision in the quoted graph.


Can 100 Machines Agree?. Agreement protocols have been typically deployed at small scale, e.g., using three to five machines. This is because these protocols seem to suffer from a sharp performance decay. More specifically, as the size of a deployment i.e., degree of replication increases, the protocol performance greatly decreases. There is not much experimental evidence for this decay in practice, however, notably for larger system sizes, e.g., beyond a handful of machines. In this paper we execute agreement protocols on up to 100 machines and observe on their performance decay. We consider well known agreement protocols part of mature systems, such as Apache ZooKeeper, etcd, and BFT Smart, as well as a chain and a novel ring based agreement protocol which we implement ourselves. We provide empirical evidence that current agreement protocols execute gracefully on 100 machines. We observe that throughput decay is initially sharp (consistent with previous observations) but intriguingly as each system grows beyond a few tens of replicas the decay dampens. For chain and ring based replication, this decay is slower than for the other systems. The positive takeaway from our evaluation is that mature agreement protocol implementations can sustain out of the box 300 to 500 requests per second when executing on 100 replicas on a wide area public cloud platform. Chain and ring based replication can reach between 4K and 11K (up to 20x improvements) depending on the fault assumptions.


From code to market Network of developers and correlated returns of cryptocurrencies. "Code is law" is the funding principle of cryptocurrencies. The security, transferability, availability and other properties of a crypto asset are determined by the code through which it is created. If code is open source, as it happens for most cryptocurrencies, this principle would prevent manipulations and grant transparency to users and traders. However, this approach considers cryptocurrencies as isolated entities thus neglecting possible connections between them. Here, we show that 4 of developers contribute to the code of more than one cryptocurrency and that the market reflects these cross asset dependencies. In particular, we reveal that the first coding event linking two cryptocurrencies through a common developer leads to the synchronisation of their returns in the following months. Our results identify a clear link between the collaborative development of cryptocurrencies and their market behaviour. More broadly, our work reveals a so far overlooked systemic dimension for the transparency of code based ecosystems and we anticipate it will be of interest to researchers, investors and regulators.


An Emergent Autonomous Flow for Mean Field Spin Glasses. We study the dynamics of symmetric and asymmetric spin glass models of size N . The analysis is in terms of the double empirical process this contains both the spins, and the field felt by each spin, at a particular time (without any knowledge of the correlation history). It is demonstrated that in the large N limit, the dynamics of the double empirical process becomes deterministic and autonomous over finite time intervals. This does not contradict the well known fact that SK spin glass dynamics is non Markovian (in the large N limit) because the empirical process has a topology that does not discern correlations in individual spins at different times. In the large N limit, the evolution of the density of the double empirical process approaches a nonlocal autonomous PDE operator Phi t . Because the emergent dynamics is autonomous, in future work one will be able to apply PDE techniques to analyze bifurcations in Phi t . Preliminary numerical results for the SK Glauber dynamics suggest that the glassy dynamical phase transition' occurs when a stable fixed point of the flow operator Phi t destabilizes.


CLEVR Parser A Graph Parser Library for Geometric Learning on Language Grounded Image Scenes. The CLEVR dataset has been used extensively in language grounded visual reasoning in Machine Learning (ML) and Natural Language Processing (NLP) domains. We present a graph parser library for CLEVR, that provides functionalities for object centric attributes and relationships extraction, and construction of structural graph representations for dual modalities. Structural order invariant representations enable geometric learning and can aid in downstream tasks like language grounding to vision, robotics, compositionality, interpretability, and computational grammar construction. We provide three extensible main components parser, embedder, and visualizer that can be tailored to suit specific learning setups. We also provide out of the box functionality for seamless integration with popular deep graph neural network (GNN) libraries. Additionally, we discuss downstream usage and applications of the library, and how it accelerates research for the NLP research community.


Bounds on the minimum distance of algebraic geometry codes defined over some families of surfaces. We prove lower bounds for the minimum distance of algebraic geometry codes over surfaces whose canonical divisor is either nef or anti strictly nef and over surfaces without irreducible curves of small genus. We sharpen these lower bounds for surfaces whose arithmetic Picard number equals one, surfaces without curves with small self intersection and fibered surfaces. Finally we specify our bounds to the case of surfaces of degree d geq 3 embedded in mathbb P) 3 .


On A Class of Degenerate And Singular Monge Amp ere Equations. In this paper we shall prove the existence, uniqueness and global H ddot o lder continuity for the Dirichlet problem of a class of Monge Amp ere type equations which may be degenerate and singular on the boundary of convex domains. We will establish a relation of the H ddot o lder exponent for the solutions with the convexity for the domains.


The general spectral radius and majorization theorem of t cone graphs with given degree sequences. The general spectral radius of a graph G , denoted by Theta(G, alpha) , is the maximal eigenvalue of M alpha (G) A(G) alpha D(G) ( alpha geq 0) , where A(G) and D(G) are the adjacency matrix and the diagonal matrix of vertex degrees of G , respectively. A graph G is called Theta alpha maximal in a class of connected simple graphs mathcal G if Theta(G, alpha) is maximal among all graphs of mathcal G . A t cone c cyclic graph is the join of a complete graph K t and a c cyclic connected simple graph. Let pi big(d 1,d 2, ldots,d n big) and pi' big(d' 1,d' 2, ldots,d' n big) be two non increasing degree sequences of t cone c cyclic graphs with n vertices. We say pi is strictly majorized by pi' , denoted by pi lhd pi' , if pi neq pi' , sum i 1 n d i sum i 1 n d i' , and sum i 1 k d i leq sum i 1 k d i' for k 1,2, ldots,n 1 . Denote by Gamma( pi,t c) the class of t cone c cyclic graphs with pi as its degree sequence. In this paper, we determine some properties of Theta alpha maximal graphs of Gamma( pi,t c) and characterize the unique Theta alpha maximal graph of Gamma( pi,t 0) big(resp. Gamma( pi,t 1) and Gamma( pi,t 2) big). Moreover, we prove that if pi lhd pi' , G and G' are the Theta alpha maximal graphs of Gamma( pi,t c) and Gamma( pi',t c) respectively, then Theta(G, alpha) Theta(G', alpha) for c in big 0,1 big , and we also consider the similar result for c 2 .


Improving VQA and its Explanations by Comparing Competing Explanations. Most recent state of the art Visual Question Answering (VQA) systems are opaque black boxes that are only trained to fit the answer distribution given the question and visual content. As a result, these systems frequently take shortcuts, focusing on simple visual concepts or question priors. This phenomenon becomes more problematic as the questions become complex that requires more reasoning and commonsense knowledge. To address this issue, we present a novel framework that uses explanations for competing answers to help VQA systems select the correct answer. By training on human textual explanations, our framework builds better representations for the questions and visual content, and then reweights confidences in the answer candidates using either generated or retrieved explanations from the training set. We evaluate our framework on the VQA X dataset, which has more difficult questions with human explanations, achieving new state of the art results on both VQA and its explanations.


A Pollution Control Problem for the Aluminum Production in Eastern Siberia Differential Game Approach. In this paper, we apply a dynamic game theoretic model and analyze the problem of pollution control in Eastern Siberia region of Russia. When carrying out the analysis we use real numerical values of parameters. It is shown that cooperation between the major pollutants can be beneficial not only for the nature but also for the respective companies.


Multimodal Routing Improving Local and Global Interpretability of Multimodal Language Analysis. The human language can be expressed through multiple sources of information known as modalities, including tones of voice, facial gestures, and spoken language. Recent multimodal learning with strong performances on human centric tasks such as sentiment analysis and emotion recognition are often black box, with very limited interpretability. In this paper we propose Multimodal Routing, which dynamically adjusts weights between input modalities and output representations differently for each input sample. Multimodal routing can identify relative importance of both individual modalities and cross modality features. Moreover, the weight assignment by routing allows us to interpret modality prediction relationships not only globally (i.e. general trends over the whole dataset), but also locally for each single input sample, meanwhile keeping competitive performance compared to state of the art methods.


ADOL Markovian approximation of rough lognormal model. In this paper we apply Markovian approximation of the fractional Brownian motion (BM), known as the Dobric Ojeda (DO) process, to the fractional stochastic volatility model where the instantaneous variance is modelled by a lognormal process with drift and fractional diffusion. Since the DO process is a semi martingale, it can be represented as an Ito diffusion. It turns out that in this framework the process for the spot price S t is a geometric BM with stochastic instantaneous volatility sigma t , the process for sigma t is also a geometric BM with stochastic speed of mean reversion and time dependent colatility of volatility, and the supplementary process calV t is the Ornstein Uhlenbeck process with time dependent coefficients, and is also a function of the Hurst exponent. We also introduce an adjusted DO process which provides a uniformly good approximation of the fractional BM for all Hurst exponents H in 0,1 but requires a complex measure. Finally, the characteristic function (CF) of log S t in our model can be found in closed form by using asymptotic expansion. Therefore, pricing options and variance swaps (by using a forward CF) can be done via FFT, which is much easier than in rough volatility models.


UniquID A Quest to Reconcile Identity Access Management and the Internet of Things. The Internet of Things (IoT) has caused a revolutionary paradigm shift in computer networking. After decades of human centered routines, where devices were merely tools that enabled human beings to authenticate themselves and perform activities, we are now dealing with a device centered paradigm the devices themselves are actors, not just tools for people. Conventional identity access management (IAM) frameworks were not designed to handle the challenges of IoT. Trying to use traditional IAM systems to reconcile heterogeneous devices and complex federations of online services (e.g., IoT sensors and cloud computing solutions) adds a cumbersome architectural layer that can become hard to maintain and act as a single point of failure. In this paper, we propose UniquID, a blockchain based solution that overcomes the need for centralized IAM architectures while providing scalability and robustness. We also present the experimental results of a proof of concept UniquID enrolment network, and we discuss two different use cases that show the considerable value of a blockchain based IAM.


A Neural Symbolic Framework for Mental Simulation. We present a neural symbolic framework for observing the environment and continuously learning visual semantics and intuitive physics to reproduce them in an interactive simulation. The framework consists of five parts, a neural symbolic hybrid network based on capsules for inverse graphics, an episodic memory to store observations, an interaction network for intuitive physics, a meta learning agent that continuously improves the framework and a querying language that acts as the framework's interface for simulation. By means of lifelong meta learning, the capsule network is expanded and trained continuously, in order to better adapt to its environment with each iteration. This enables it to learn new semantics using a few shot approach and with minimal input from an oracle over its lifetime. From what it learned through observation, the part for intuitive physics infers all the required physical properties of the objects in a scene, enabling predictions. Finally, a custom query language ties all parts together, which allows to perform various mental simulation tasks, such as navigation, sorting and simulation of a game environment, with which we illustrate the potential of our novel approach.


Attacks on Image Encryption Schemes for Privacy Preserving Deep Neural Networks. Privacy preserving machine learning is an active area of research usually relying on techniques such as homomorphic encryption or secure multiparty computation. Recent novel encryption techniques for performing machine learning using deep neural nets on images have recently been proposed by Tanaka and Sirichotedumrong, Kinoshita, and Kiya. We present new chosen plaintext and ciphertext only attacks against both of these proposed image encryption schemes and demonstrate the attacks' effectiveness on several examples.


Cosmos QA Machine Reading Comprehension with Contextual Commonsense Reasoning. Understanding narratives requires reading between the lines, which in turn, requires interpreting the likely causes and effects of events, even when they are not mentioned explicitly. In this paper, we introduce Cosmos QA, a large scale dataset of 35,600 problems that require commonsense based reading comprehension, formulated as multiple choice questions. In stark contrast to most existing reading comprehension datasets where the questions focus on factual and literal understanding of the context paragraph, our dataset focuses on reading between the lines over a diverse collection of people's everyday narratives, asking such questions as "what might be the possible reason of ...?", or "what would have happened if ..." that require reasoning beyond the exact text spans in the context. To establish baseline performances on Cosmos QA, we experiment with several state of the art neural architectures for reading comprehension, and also propose a new architecture that improves over the competitive baselines. Experimental results demonstrate a significant gap between machine (68.4 ) and human performance (94 ), pointing to avenues for future research on commonsense machine comprehension. Dataset, code and leaderboard is publicly available at https wilburone.github.io cosmos.


Generation and characterization of few pulse attosecond pulse trains at 100kHz repetition rate. Many experiments in attosecond science will benefit from attosecond pulses at high repetition rates with sufficient photon flux for pump probe experiments. We use 7fs, 800nm pulses from a non collinear optical parametric chirped pulse amplification (NOPCPA) laser system to generate few pulse attosecond pulse trains (APTs) in the extreme ultraviolet (XUV) at a repetition rate of 100kHz. The pulse trains have been fully characterized by recording FROG CRAB (Frequency Resolved Optical Gating for Complete Reconstruction of Attosecond Bursts) traces with a velocity map imaging spectrometer.


Out the Window A Crowd Sourced Dataset for Activity Classification in Security Video. The Out the Window (OTW) dataset is a crowdsourced activity dataset containing 5,668 instances of 17 activities from the NIST Activities in Extended Video (ActEV) challenge. These videos are crowdsourced from workers on the Amazon Mechanical Turk using a novel scenario acting strategy, which collects multiple instances of natural activities per scenario. Turkers are instructed to lean their mobile device against an upper story window overlooking an outdoor space, walk outside to perform a scenario involving people, vehicles and objects, and finally upload the video to us for annotation. Performance evaluation for activity classification on VIRAT Ground 2.0 shows that the OTW dataset provides an 8.3 improvement in mean classification accuracy, and a 12.5 improvement on the most challenging activities involving people with vehicles.


A Computer Assisted Study of Red Coral Population Dynamics. We consider a 13 dimensional age structured discrete red coral population model varying with respect to a fitness parameter. Our numerical results give a bifurcation diagram of both equilibria and stable invariant curves of orbits. We observe that not only for low levels of fitness, but also for high levels of fitness, populations are extremely vulnerable, in that they spend long time periods near extinction. We then use computer assisted proofs techniques to rigorously validate the set of regular and bifurcation fixed points that have been found numerically.


Enhancement of the Coloring Invariant for Folded Molecular Chains. Folded linear molecular chains are ubiquitous in biology. Folding is mediated by intra chain interactions that "glue" two or more regions of a chain. The resulting fold topology is widely believed to be a determinant of biomolecular properties and function. Recently, knot theory has been extended to describe the topology of folded linear chains such as proteins and nucleic acids. To classify and distinguish chain topologies, algebraic structure of quandles has been adapted and applied. However, the approach is limited as apparently distinct topologies may end up having the same number of colorings. Here, we enhance the resolving power of the quandle coloring approach by introducing Boltzmann weights. We demonstrate that the enhanced coloring invariants can distinguish fold topologies with an improved resolution.


Lineability and modes of convergence. In this paper we look for the existence of large linear and algebraic structures of sequences of measurable functions with different modes of convergence. Concretely, the algebraic size of the family of sequences that are convergent in measure but not a.e. pointwise, uniformly but not pointwise convergent, and uniformly convergent but not in L 1 norm, are analyzed. These findings extend and complement a number of earlier results by several authors.


Matching numbers and dimension of edge ideals. Let G be a finite simple graph on the vertex set V(G) x 1 , ldots, x n and match (G) , min match (G) and ind match (G) the matching number, minimum matching number and induced matching number of G , respectively. Let K V(G) K x 1 , ldots, x n denote the polynomial ring over a field K and I(G) subset K V(G) the edge ideal of G . The relationship between these graph theoretic invariants and ring theoretic invariants of the quotient ring K V(G) I(G) has been studied. In the present paper, we study the relationship between match (G) , min match (G) , ind match (G) and dim K V(G) I(G) .


Towards Generalization and Data Efficient Learning of Deep Robotic Grasping. Deep reinforcement learning (DRL) has been proven to be a powerful paradigm for learning complex control policy autonomously. Numerous recent applications of DRL in robotic grasping have successfully trained DRL robotic agents end to end, mapping visual inputs into control instructions directly, but the amount of training data required may hinder these applications in practice. In this paper, we propose a DRL based robotic visual grasping framework, in which visual perception and control policy are trained separately rather than end to end. The visual perception produces physical descriptions of grasped objects and the policy takes use of them to decide optimal actions based on DRL. Benefiting from the explicit representation of objects, the policy is expected to be endowed with more generalization power over new objects and environments. In addition, the policy can be trained in simulation and transferred in real robotic system without any further training. We evaluate our framework in a real world robotic system on a number of robotic grasping tasks, such as semantic grasping, clustered object grasping, moving object grasping. The results show impressive robustness and generalization of our system.


SMART Simultaneous Multi Agent Recurrent Trajectory Prediction. We propose advances that address two key challenges in future trajectory prediction (i) multimodality in both training data and predictions and (ii) constant time inference regardless of number of agents. Existing trajectory predictions are fundamentally limited by lack of diversity in training data, which is difficult to acquire with sufficient coverage of possible modes. Our first contribution is an automatic method to simulate diverse trajectories in the top view. It uses pre existing datasets and maps as initialization, mines existing trajectories to represent realistic driving behaviors and uses a multi agent vehicle dynamics simulator to generate diverse new trajectories that cover various modes and are consistent with scene layout constraints. Our second contribution is a novel method that generates diverse predictions while accounting for scene semantics and multi agent interactions, with constant time inference independent of the number of agents. We propose a convLSTM with novel state pooling operations and losses to predict scene consistent states of multiple agents in a single forward pass, along with a CVAE for diversity. We validate our proposed multi agent trajectory prediction approach by training and testing on the proposed simulated dataset and existing real datasets of traffic scenes. In both cases, our approach outperforms SOTA methods by a large margin, highlighting the benefits of both our diverse dataset simulation and constant time diverse trajectory prediction methods.


Extension of a Diophantine triple with the property D(4) . In this paper, we give an upper bound on the number of extensions of a triple to a quadruple for the Diophantine m tuples with the property D(4) . We also confirm the conjecture of the uniqueness of such an extension in some special cases.


Pioneer dataset and automatic recognition of Urdu handwritten characters using a deep autoencoder and convolutional neural network. Automatic recognition of Urdu handwritten digits and characters, is a challenging task. It has applications in postal address reading, bank's cheque processing, and digitization and preservation of handwritten manuscripts from old ages. While there exists a significant work for automatic recognition of handwritten English characters and other major languages of the world, the work done for Urdu lan guage is extremely insufficient. This paper has two goals. Firstly, we introduce a pioneer dataset for handwritten digits and characters of Urdu, containing samples from more than 900 individuals. Secondly, we report results for automatic recog nition of handwritten digits and characters as achieved by using deep auto encoder network and convolutional neural network. More specifically, we use a two layer and a three layer deep autoencoder network and convolutional neural network and evaluate the two frameworks in terms of recognition accuracy. The proposed framework of deep autoencoder can successfully recognize digits and characters with an accuracy of 97 for digits only, 81 for characters only and 82 for both digits and characters simultaneously. In comparison, the framework of convolutional neural network has accuracy of 96.7 for digits only, 86.5 for characters only and 82.7 for both digits and characters simultaneously. These frameworks can serve as baselines for future research on Urdu handwritten text.


The 1 5 th Rule with Rollbacks On Self Adjustment of the Population Size in the (1 ( lambda, lambda)) GA. Self adjustment of parameters can significantly improve the performance of evolutionary algorithms. A notable example is the (1 ( lambda, lambda)) genetic algorithm, where the adaptation of the population size helps to achieve the linear runtime on the OneMax problem. However, on problems which interfere with the assumptions behind the self adjustment procedure, its usage can lead to performance degradation compared to static parameter choices. In particular, the one fifth rule, which guides the adaptation in the example above, is able to raise the population size too fast on problems which are too far away from the perfect fitness distance correlation. We propose a modification of the one fifth rule in order to have less negative impact on the performance in scenarios when the original rule reduces the performance. Our modification, while still having a good performance on OneMax, both theoretically and in practice, also shows better results on linear functions with random weights and on random satisfiable MAX SAT instances.


Near Optimal Primal Dual Algorithms for Quantity Based Network Revenue Management. We study the canonical quantity based network revenue management (NRM) problem where the decision maker must irrevocably accept or reject each arriving customer request with the goal of maximizing the total revenue given limited resources. The exact solution to the problem by dynamic programming is computationally intractable due to the well known curse of dimensionality. Existing works in the literature make use of the solution to the deterministic linear program (DLP) to design asymptotically optimal algorithms. Those algorithms rely on repeatedly solving DLPs to achieve near optimal regret bounds. It is, however, time consuming to repeatedly compute the DLP solutions in real time, especially in large scale problems that may involve hundreds of millions of demand units. In this paper, we propose innovative algorithms for the NRM problem that are easy to implement and do not require solving any DLPs. Our algorithm achieves a regret bound of O( log k) , where k is the system size. To the best of our knowledge, this is the first NRM algorithm that (i) has an o( sqrt k ) asymptotic regret bound, and (ii) does not require solving any DLPs.


A Framework for Data Driven Computational Dynamics Based on Nonlinear Optimization. In this article, we present an extension of the formulation recently developed by the authors (A Framework for Data Driven Computational Mechanics Based on Nonlinear Optimization, arXiv 1910.12736 math.NA ) to the structural dynamics setting. Inspired by a structure preserving family of variational integrators, our new formulation relies on a discrete balance equation that establishes the dynamic equilibrium. From this point of departure, we first derive an "exact" discrete continuous nonlinear optimization problem that works directly with data sets. We then develop this formulation further into an "approximate" nonlinear optimization problem that relies on a general constitutive model. This underlying model can be identified from a data set in an offline phase. To showcase the advantages of our framework, we specialize our methodology to the case of a geometrically exact beam formulation that makes use of all elements of our approach. We investigate three numerical examples of increasing difficulty that demonstrate the excellent computational behavior of the proposed framework and motivate future research in this direction.


PC DARTS Partial Channel Connections for Memory Efficient Architecture Search. Differentiable architecture search (DARTS) provided a fast solution in finding effective network architectures, but suffered from large memory and computing overheads in jointly training a super network and searching for an optimal architecture. In this paper, we present a novel approach, namely, Partially Connected DARTS, by sampling a small part of super network to reduce the redundancy in exploring the network space, thereby performing a more efficient search without comprising the performance. In particular, we perform operation search in a subset of channels while bypassing the held out part in a shortcut. This strategy may suffer from an undesired inconsistency on selecting the edges of super net caused by sampling different channels. We alleviate it using edge normalization, which adds a new set of edge level parameters to reduce uncertainty in search. Thanks to the reduced memory cost, PC DARTS can be trained with a larger batch size and, consequently, enjoys both faster speed and higher training stability. Experimental results demonstrate the effectiveness of the proposed method. Specifically, we achieve an error rate of 2.57 on CIFAR10 with merely 0.1 GPU days for architecture search, and a state of the art top 1 error rate of 24.2 on ImageNet (under the mobile setting) using 3.8 GPU days for search. Our code has been made available at https github.com yuhuixu1993 PC DARTS.


Dynamic Resource Allocation in Co Located and Cell Free Massive MIMO. In this paper, we study joint power control and scheduling in uplink massive multiple input multiple output (MIMO) systems with randomly arriving data traffic. We consider both co located and Cell Free (CF) Massive MIMO, where the difference lies in whether the antennas are co located at the base station or spread over a wide network area. The data is generated at each user according to an individual stochastic process. Using Lyapunov optimization techniques, we develop a dynamic scheduling algorithm (DSA), which decides at each time slot the amount of data to admit to the transmission queues and the transmission rates over the wireless channel. The proposed algorithm optimizes the long term user throughput under various fairness policies while keeping the transmission queues stable. Simulation results show that the state of the art power control schemes developed for Massive MIMO with infinite backlogs can fail to stabilize the system even when the data arrival rates are within the network capacity region. Our proposed DSA shows advantage in providing finite delay with performance optimization whenever the network can be stabilized.


Commutators and images of noncommutative polynomials. Let A be an algebra and let f be a nonconstant noncommutative polynomial. In the first part of the paper, we consider the relationship between A,A , the linear span of commutators in A , and span f(A) , the linear span of the image of f in A . In particular, we show that A,A A implies span f(A) A . In the second part, we establish some Waring type results for images of polynomials. For example, we show that if C is a commutative unital algebra over a field F of characteristic 0 , A is the matrix algebra M n(C) , and the polynomial f is neither an identity nor a central polynomial of M n(F) , then every commutator in A can be written as a difference of two elements, each of which is a sum of 7788 elements from f(A) (if C F is an algebraically closed field, then 4 elements suffice). Similar results are obtained for some other algebras, in particular for the algebra B(H) of all bounded linear operators on a Hilbert space H .


Some conditions from a finite regular semigroup to a finite inverse semigroup. In 1978, Munn proved that a bounded complex representation of an inverse semigroup is semiunitary and completely reducible. We consider the converse question in the finite case. We provide some sufficient conditions from a finite regular semigroup to a finite inverse semigroup.


Augmenting Machine Learning with Information Retrieval to Recommend Real Cloned Code Methods for Code Completion. Software developers frequently reuse source code from repositories as it saves development time and effort. Code clones accumulated in these repositories hence represent often repeated functionalities and are candidates for reuse in an exploratory or rapid development. In previous work, we introduced DeepClone, a deep neural network model trained by fine tuning GPT 2 model over the BigCloneBench dataset to predict code clone methods. The probabilistic nature of DeepClone output generation can lead to syntax and logic errors that requires manual editing of the output for final reuse. In this paper, we propose a novel approach of applying an information retrieval (IR) technique on top of DeepClone output to recommend real clone methods closely matching the predicted output. We have quantitatively evaluated our strategy, showing that the proposed approach significantly improves the quality of recommendation.


Reactive miscible displacement of light oil in porous media. We develop a theory for the problem of high pressure air injection into deep reservoirs containing light oil. Under these conditions, the injected fluid (oxygen inert components) is completely miscible with the oil in the reservoir. Moreover, exothermic reactions between dissolved oxygen and oil are possible. We use Koval's model to account for the miscibility of the components, such that the fractional flow functions resemble the ones from Buckley Leverett flow. This allows to decompose the solution of this problem into a series of waves. We then proceed to obtain full analytical solutions in each wave. Of particular interest is the case where the combustion wave presents a singularity in its internal wave profile. Evaluation of the variables of the problem at the singular point determines the macroscopic parameters of the wave, i.e., combustion temperature, wave speed and downstream oil fraction. The waves structure was observed previously for reactive immiscible displacement and we describe it here for the first time for reactive miscible displacement of oil. We validate the developed theory using numerical simulations.


On the existence of Hopf bifurcations in the sequential and distributive double phosphorylation cycle. Protein phosphorylation cycles are important mechanisms of the post translational modification of a protein and as such an integral part of intracellular signaling and control. We consider the sequential phosphorylation and dephosphorylation of a protein at two binding sites. While it is known that proteins where phosphorylation is processive and dephosphorylation is distributive admit oscillations (for some value of the rate constants and total concentrations) it is not known whether or not this is the case if both phosphorylation and dephosphorylation are distributive. We study four simplified mass action models of sequential and distributive phosphorylation and show that for each of those there do not exist rate constants and total concentrations where a Hopf bifurcation occurs. To arrive at this result we use convex parameters to parameterize the steady state and Hurwitz matrices.


Scikit network Graph Analysis in Python. Scikit network is a Python package inspired by scikit learn for the analysis of large graphs. Graphs are represented by their adjacency matrix in the sparse CSR format of SciPy. The package provides state of the art algorithms for ranking, clustering, classifying, embedding and visualizing the nodes of a graph. High performance is achieved through a mix of fast matrix vector products (using SciPy), compiled code (using Cython) and parallel processing. The package is distributed under the BSD license, with dependencies limited to NumPy and SciPy. It is compatible with Python 3.6 and newer. Source code, documentation and installation instructions are available online.


Partial Quantifier Elimination With Learning. We consider a modification of the Quantifier Elimination (QE) problem called Partial QE (PQE). In PQE, only a small part of the formula is taken out of the scope of quantifiers. The appeal of PQE is that many verification problems, e.g. equivalence checking and model checking, reduce to PQE and the latter is much easier than complete QE. Earlier, we introduced a PQE algorithm based on the machinery of D sequents. A D sequent is a record stating that a clause is redundant in a quantified CNF formula in a specified subspace. To make this algorithm efficient, it is important to reuse learned D sequents. However, reusing D sequents is not as easy as conflict clauses in SAT solvers because redundancy is a structural rather than a semantic property. Earlier, we modified the definition of D sequents to enable their safe reusing. In this paper, we present a PQE algorithm based on new D sequents. It is different from its predecessor in two aspects. First, the new algorithm can learn and reuse D sequents. Second, it proves clauses redundant one by one and thus backtracks as soon as the current target clause is proved redundant in the current subspace. This makes the new PQE algorithm similar to a SAT solver that backtracks as soon as just one clause is falsified. We show experimentally that the new PQE algorithm outperforms its predecessor.


An Integration of UTAUT and Task Technology Fit Frameworks for Assessing the Acceptance of Clinical Decision Support Systems in the Context of a Developing Country. This paper is to create a basis of theoretical contribution for a new PhD thesis in the area of Clinical Decision Support Systems (CDSS) acceptance. Over the past three years, we conducted qualitative research into three distinctive phases to develop an extended Task Technology Fit (TTF) Framework. These phases are for initiating requirement generation of the framework, discovering the factors of the framework through perspectives, and evaluating the new proposed framework. The new condition is related to developing country in which various sectors such as healthcare is mostly under attention. We conduct a new inspective for assisting decisions to support technology and its usefulness in this sector to integrate with other frameworks for assisting the value, use, and how can be better accepted in the context of healthcare professionals.


Arbitrary positive powers of semicirculant and r circulant matrices. We provide a novel recursive method, which does not require any assumption, to compute the entries of the kth power of a semicirculant matrix. As an application, a method for computing the entries of the kth power of r circulant matrices is also presented.


The irreducibility of some Wronskian Hermite polynomials. We study the irreducibility of Wronskian Hermite polynomials labelled by partitions. It is known that these polynomials factor as a power of x times a remainder polynomial. We show that the remainder polynomial is irreducible for the partitions (n, m) with m 2, and (n, n) when n 1 is a square. Our main tools are two theorems that we prove for all partitions. The first result gives a sharp upper bound for the slope of the edges of the Newton polygon for the remainder polynomial. The second result is a Schur type congruence for Wronskian Hermite polynomials. We also explain how irreducibility determines the number of real zeros of Wronskian Hermite polynomials, and prove Veselov's conjecture for partitions of the form (n, k, k 1, ..., 1).


Event Based Angular Velocity Regression with Spiking Networks. Spiking Neural Networks (SNNs) are bio inspired networks that process information conveyed as temporal spikes rather than numeric values. A spiking neuron of an SNN only produces a spike whenever a significant number of spikes occur within a short period of time. Due to their spike based computational model, SNNs can process output from event based, asynchronous sensors without any pre processing at extremely lower power unlike standard artificial neural networks. This is possible due to specialized neuromorphic hardware that implements the highly parallelizable concept of SNNs in silicon. Yet, SNNs have not enjoyed the same rise of popularity as artificial neural networks. This not only stems from the fact that their input format is rather unconventional but also due to the challenges in training spiking networks. Despite their temporal nature and recent algorithmic advances, they have been mostly evaluated on classification problems. We propose, for the first time, a temporal regression problem of numerical values given events from an event camera. We specifically investigate the prediction of the 3 DOF angular velocity of a rotating event camera with an SNN. The difficulty of this problem arises from the prediction of angular velocities continuously in time directly from irregular, asynchronous event based input. Directly utilising the output of event cameras without any pre processing ensures that we inherit all the benefits that they provide over conventional cameras. That is high temporal resolution, high dynamic range and no motion blur. To assess the performance of SNNs on this task, we introduce a synthetic event camera dataset generated from real world panoramic images and show that we can successfully train an SNN to perform angular velocity regression.


COVID 19 On the quarantine duration after short visits to high risk regions. A simple Monte Carlo method will be put forward herein, to enable the extraction of an estimate for the quarantine duration, applicable to visitors to high risk regions. Results will be obtained on the basis of an analysis of the upper tail of the cumulative distribution function of the time span between the departure of the travellers from the place where the infection occurs and the time instant when COVID 19 infections may currently be detected. As expected, the quarantine duration is a decreasing function of the fraction of the infected travellers, which one is prepared to identify as acceptable risk'. The analysis suggests that a maximal 5 risk (of new infections originating from subjects who become infective after their quarantine is lifted) may be associated with a minimal quarantine duration of about eight days, 1 with about twelve, and 0.1 with about sixteen. Unless the distribution of the duration of short (typically, up to three weeks) travels departs significantly from the shape assumed in this study, the results of the present analysis do not provide support for the plans to shorten the quarantine duration of about ten days to two weeks, which currently applies to travellers entering most European countries from regions with a high risk of infection.


Hadamard powers of rank two, doubly nonnegative matrices. We study ranks of the r textrm th Hadamard powers of doubly nonnegative matrices and show that the matrix A circ r is positive definite for every n times n doubly nonnegative matrix A and for every r n 2 if and only if no column of A is a scalar multiple of any other column of A. A particular emphasis is given to the study of rank, positivity and monotonicity of Hadamard powers of rank two, positive semidefinite matrices that have all entries positive.


Complementary problems with polynomial data. Given polynomial maps f, g colon mathbb R n to mathbb R n, we consider the em polynomial complementary problem of finding a vector x in mathbb R n such that begin equation f(x) ge 0, quad g(x) ge 0, quad textrm and quad langle f(x), g(x) rangle 0. end equation In this paper, we present various properties on the solution set of the problem, including genericity, nonemptiness, compactness, uniqueness as well as error bounds with exponents explicitly determined. These strengthen and generalize some previously known results, and hence broaden the boundary knowledge of nonlinear complementarity problems as well.


The effect of the stochasticity of photoionization on 3D streamer simulations. Positive streamer discharges require a source of free electrons ahead of them for their growth. In air, these electrons are typically provided by photoionization. Here we investigate how stochastic fluctuations due to the discreteness of ionizing photons affect positive streamers in air. We simulate positive streamers between two planar electrodes with a 3D plasma fluid model, using both a stochastic and a continuum method for photoionization. With stochastic photoionization, fluctuations are visible in the streamer's direction, maximal electric field, velocity, and electron density. The streamers do not branch, and we find good agreement between the averaged stochastic results and the results with continuum photoionization. The streamers stay roughly axisymmetric, and we show that results obtained with an axisymmetric model indeed agree well with the 3D results. However, we find that positive streamers are sensitive to the amount of photoionization. When the amount of photoionization is doubled, there is even better agreement between the stochastic and continuum results, but with half the amount of photoionization, stochastic fluctuations become more important and streamer branching starts to occur.


Robust Bang Off Bang Low Thrust Guidance Using Model Predictive Static Programming. Model Predictive Static Programming (MPSP) was always used under the assumption of continuous control, which impedes it for applications with bang off bang control directly. In this paper, MPSP is employed for the first time as a guidance scheme for low thrust transfers with bang off bang control where the fuel optimal trajectory is used as the nominal solution. In our method, dynamical equations in Cartesian coordinates are augmented by the mass costate equation, while the unconstrained velocity costate vector is used as control variable, and is expressed as a combination of Fourier basis functions with corresponding weights. A two loop MPSP algorithm is designed where the weights and the initial mass costate are updated in the inner loop and continuation is conducted on the outer loop in case of large perturbations. The sensitivity matrix (SM) is recursively calculated using analytical derivatives and SM at switching points is compensated based on calculus of variations. An sample interplanetary CubeSat mission to an asteroid is used as study case to illustrate the effectiveness of the method developed.


The shape of a free boundary driven by a line of fast diffusion. We complete the description, initiated in 6 , of a free boundary travelling at constant speed in a half plane, where the propagation is controlled by a line having a large diffusion on its own. The main result of this work is that the free boundary is asymptotic to a line at infinity, whose angle to the horizontal is dicatated by the velocity of the wave imposed by the line. This helps understanding some rather counter intuitive numerical simulations of 8 .


Dual Energy Subtraction Efficiency Development of an Objective Quality Metric for Tissue Subtraction Radiography. We propose a new objective numerical figure of merit to aid in the evaluation and comparison of tissue selective images generated from dual energy radiography systems. A metric is developed through identification of the requirements of a successful objective metric and analyzed in a variety of scenarios. The proposed Dual Energy Subtraction Efficiency (DSE) metric adequately describes a multitude of image properties for all simulated image scenarios, indicating its ability to accurately and objectively describe image quality. The DSE and its measurement method described can become a tool to characterize dual energy radiographic image quality objectively and quantitatively, allowing for improved system comparison, development, and optimization.


Metrics and methods for a systematic comparison of fairness aware machine learning algorithms. Understanding and removing bias from the decisions made by machine learning models is essential to avoid discrimination against unprivileged groups. Despite recent progress in algorithmic fairness, there is still no clear answer as to which bias mitigation approaches are most effective. Evaluation strategies are typically use case specific, rely on data with unclear bias, and employ a fixed policy to convert model outputs to decision outcomes. To address these problems, we performed a systematic comparison of a number of popular fairness algorithms applicable to supervised classification. Our study is the most comprehensive of its kind. It utilizes three real and four synthetic datasets, and two different ways of converting model outputs to decisions. It considers fairness, predictive performance, calibration quality, and speed of 28 different modelling pipelines, corresponding to both fairness unaware and fairness aware algorithms. We found that fairness unaware algorithms typically fail to produce adequately fair models and that the simplest algorithms are not necessarily the fairest ones. We also found that fairness aware algorithms can induce fairness without material drops in predictive power. Finally, we found that dataset idiosyncracies (e.g., degree of intrinsic unfairness, nature of correlations) do affect the performance of fairness aware approaches. Our results allow the practitioner to narrow down the approach(es) they would like to adopt without having to know in advance their fairness requirements.


Fine grained Classification via Categorical Memory Networks. Motivated by the desire to exploit patterns shared across classes, we present a simple yet effective class specific memory module for fine grained feature learning. The memory module stores the prototypical feature representation for each category as a moving average. We hypothesize that the combination of similarities with respect to each category is itself a useful discriminative cue. To detect these similarities, we use attention as a querying mechanism. The attention scores with respect to each class prototype are used as weights to combine prototypes via weighted sum, producing a uniquely tailored response feature representation for a given input. The original and response features are combined to produce an augmented feature for classification. We integrate our class specific memory module into a standard convolutional neural network, yielding a Categorical Memory Network. Our memory module significantly improves accuracy over baseline CNNs, achieving competitive accuracy with state of the art methods on four benchmarks, including CUB 200 2011, Stanford Cars, FGVC Aircraft, and NABirds.


Bi Parametric Operator Preconditioning. We extend the operator preconditioning framework R. Hiptmair, Comput. Math. with Appl. 52 (2006), pp. 699 706 to Petrov Galerkin methods while accounting for parameter dependent perturbations of both variational forms and their preconditioners, as occurs when performing numerical approximations. By considering different perturbation parameters for the original form and its preconditioner, our bi parametric abstract setting leads to robust and controlled schemes. For Hilbert spaces, we derive exhaustive linear and super linear convergence estimates for iterative solvers, such as h independent convergence bounds, when preconditioning with low accuracy or, equivalently, with highly compressed approximations.


Optimal liquidation for a risk averse investor in a one sided limit order book driven by a Levy process. In a one sided limit order book, satisfying some realistic assumptions, where the unaffected price process follows a Levy process, we consider a market agent that wants to liquidate a large position of shares. We assume that the agent has constant absolute risk aversion and aims at maximising the expected utility of the cash position at the end of time. The agent is then faced with the problem of balancing the market risk and the cost of a rapid execution. In particular we are interested in how the agent should go about optimally submitting orders. Since liquidation normally takes place within a short period of time, modelling the risk as a Levy process should provide a realistic model with good statistical fit to observed market data, and thus the model should provide a realistic reflection of the agent's market risk. We reduce the optimisation problem to a deterministic two dimensional singular problem, to which we are able to derive an explicit solution in terms of the model data. In particular we find an expression for the optimal intervention boundary, which completely characterise the optimal liquidation strategy.


Decomposing planar graphs into graphs with degree restrictions. Given a graph G , a decomposition of G is a partition of its edges. A graph is (d, h) decomposable if its edge set can be partitioned into a d degenerate graph and a graph with maximum degree at most h . For d le 4 , we are interested in the minimum integer h d such that every planar graph is (d,h d) decomposable. It was known that h 3 le 4 , h 2 le 8 , and h 1 infty . This paper proves that h 4 1, h 3 2 , and 4 le h 2 le 6 .


Deciding Differential Privacy for Programs with Finite Inputs and Outputs. Differential privacy is a de facto standard for statistical computations over databases that contain private data. The strength of differential privacy lies in a rigorous mathematical definition that guarantees individual privacy and yet allows for accurate statistical results. Thanks to its mathematical definition, differential privacy is also a natural target for formal analysis. A broad line of work uses logical methods for proving privacy. However, these methods are not complete, and only partially automated. A recent and complementary line of work uses statistical methods for finding privacy violations. However, the methods only provide statistical guarantees (but no proofs). We propose the first decision procedure for checking the differential privacy of a non trivial class of probabilistic computations. Our procedure takes as input a program P parametrized by a privacy budget epsilon , and either proves differential privacy for all possible values of epsilon or generates a counterexample. In addition, our procedure applies both to epsilon differential privacy and ( epsilon, delta) differential privacy. Technically, the decision procedure is based on a novel and judicious encoding of the semantics of programs in our class into a decidable fragment of the first order theory of the reals with exponentiation. We implement our procedure and use it for (dis)proving privacy bounds for many well known examples, including randomized response, histogram, report noisy max and sparse vector.


Big Green at WNUT 2020 Shared Task 1 Relation Extraction as Contextualized Sequence Classification. Relation and event extraction is an important task in natural language processing. We introduce a system which uses contextualized knowledge graph completion to classify relations and events between known entities in a noisy text environment. We report results which show that our system is able to effectively extract relations and events from a dataset of wet lab protocols.


An order adaptive compact approximation Taylor method for systems of conservation laws. We present a new family of high order shock capturing finite difference numerical methods for systems of conservation laws. These methods, called Adaptive Compact Approximation Taylor (ACAT) schemes, use centered (2p 1) point stencils, where p may take values in 1, 2, dots, P according to a new family of smoothness indicators in the stencils. The methods are based on a combination of a robust first order scheme and the Compact Approximate Taylor (CAT) methods of order 2p order, p 1,2, dots, P so that they are first order accurate near discontinuities and have order 2p in smooth regions, where (2p 1) is the size of the biggest stencil in which large gradients are not detected. CAT methods, introduced in cite CP2019 , are an extension to nonlinear problems of the Lax Wendroff methods in which the Cauchy Kovalesky (CK) procedure is circumvented following the strategy introduced in cite ZBM2017 that allows one to compute time derivatives in a recursive way using high order centered differentiation formulas combined with Taylor expansions in time. The expression of ACAT methods for 1D and 2D systems of balance laws are given and the performance is tested in a number of test cases for several linear and nonlinear systems of conservation laws, including Euler equations for gas dynamics.


Image Deconvolution with Deep Image and Kernel Priors. Image deconvolution is the process of recovering convolutional degraded images, which is always a hard inverse problem because of its mathematically ill posed property. On the success of the recently proposed deep image prior (DIP), we build an image deconvolution model with deep image and kernel priors (DIKP). DIP is a learning free representation which uses neural net structures to express image prior information, and it showed great success in many energy based models, e.g. denoising, super resolution, inpainting. Instead, our DIKP model uses such priors in image deconvolution to model not only images but also kernels, combining the ideas of traditional learning free deconvolution methods with neural nets. In this paper, we show that DIKP improve the performance of learning free image deconvolution, and we experimentally demonstrate this on the standard benchmark of six standard test images in terms of PSNR and visual effects.


Relative g noncommuting graph of finite groups. Let G be a finite group. For a fixed element g in G and a given subgroup H of G , the relative g noncommuting graph of G is a simple undirected graph whose vertex set is G and two vertices x and y are adjacent if x in H or y in H and x,y neq g, g 1 . We denote this graph by Gamma H, G g . In this paper, we obtain computing formulae for degree of any vertex in Gamma H, G g and characterize whether Gamma H, G g is a tree, star graph, lollipop or a complete graph together with some properties of Gamma H, G g involving isomorphism of graphs. We also present certain relations between the number of edges in Gamma H, G g and certain generalized commuting probabilities of G which give some computing formulae for the number of edges in Gamma H, G g . Finally, we conclude this paper by deriving some bounds for the number of edges in Gamma H, G g .


Model Driven Engineering for Data Protection and Privacy Application and Experience with GDPR. In Europe and indeed worldwide, the General Data Protection Regulation (GDPR) provides protection to individuals regarding their personal data in the face of new technological developments. GDPR is widely viewed as the benchmark for data protection and privacy regulations that harmonizes data privacy laws across Europe. Although the GDPR is highly beneficial to individuals, it presents significant challenges for organizations monitoring or storing personal information. Since there is currently no automated solution with broad industrial applicability, organizations have no choice but to carry out expensive manual audits to ensure GDPR compliance. In this paper, we present a complete GDPR UML model as a first step towards designing automated methods for checking GDPR compliance. Given that the practical application of the GDPR is influenced by national laws of the EU Member States, we suggest a two tiered description of the GDPR, generic and specialized. In this paper, we provide (1) the GDPR conceptual model we developed with complete traceability from its classes to the GDPR, (2) a glossary to help understand the model, (3) the plain English description of 35 compliance rules derived from GDPR along with their encoding in OCL, and (4) the set of 20 variations points derived from GDPR to specialize the generic model. We further present the challenges we faced in our modeling endeavor, the lessons we learned from it, and future directions for research.


Lower bound on the size of a quasirandom forcing set of permutations. A set S of permutations is forcing if for any sequence Pi i i in mathbb N of permutations where the density d( pi, Pi i) converges to frac 1 pi ! for every permutation pi in S , it holds that Pi i i in mathbb N is quasirandom. Graham asked whether there exists an integer k such that the set of all permutations of order k is forcing this has been shown to be true for any k ge 4 . In particular, the set of all twenty four permutations of order 4 is forcing. We provide the first non trivial lower bound on the size of a forcing set of permutations every forcing set of permutations (with arbitrary orders) contains at least four permutations.


Lowest degree robust finite element scheme for a fourth order elliptic singular perturbation problem on rectangular grids. In this paper, a piecewise quadratic nonconforming finite element method on rectangular grids for a fourth order elliptic singular perturbation problem is presented. This proposed method is robustly convergent with respect to the perturbation parameter. Numerical results are presented to verify the theoretical findings. The new method uses piecewise quadratic polynomials, and is of the lowest degree possible. Optimal order approximation property of the finite element space is proved by means of a locally averaged interpolation operator newly constructed. This interpolator, however, is not a projection. Indeed, we establish a general theory and show that no locally defined interpolation associated with the locally supported basis functions can be projective for the finite element space in use. Particularly, the general theory gives an answer to a long standing open problem presented in Demko, J. Approx. Theory, bf 43 (2) 151 156, 1985 .


Explainable AI for Interpretable Credit Scoring. With the ever growing achievements in Artificial Intelligence (AI) and the recent boosted enthusiasm in Financial Technology (FinTech), applications such as credit scoring have gained substantial academic interest. Credit scoring helps financial experts make better decisions regarding whether or not to accept a loan application, such that loans with a high probability of default are not accepted. Apart from the noisy and highly imbalanced data challenges faced by such credit scoring models, recent regulations such as the right to explanation' introduced by the General Data Protection Regulation (GDPR) and the Equal Credit Opportunity Act (ECOA) have added the need for model interpretability to ensure that algorithmic decisions are understandable and coherent. An interesting concept that has been recently introduced is eXplainable AI (XAI), which focuses on making black box models more interpretable. In this work, we present a credit scoring model that is both accurate and interpretable. For classification, state of the art performance on the Home Equity Line of Credit (HELOC) and Lending Club (LC) Datasets is achieved using the Extreme Gradient Boosting (XGBoost) model. The model is then further enhanced with a 360 degree explanation framework, which provides different explanations (i.e. global, local feature based and local instance based) that are required by different people in different situations. Evaluation through the use of functionallygrounded, application grounded and human grounded analysis show that the explanations provided are simple, consistent as well as satisfy the six predetermined hypotheses testing for correctness, effectiveness, easy understanding, detail sufficiency and trustworthiness.


Proxy Anchor Loss for Deep Metric Learning. Existing metric learning losses can be categorized into two classes pair based and proxy based losses. The former class can leverage fine grained semantic relations between data points, but slows convergence in general due to its high training complexity. In contrast, the latter class enables fast and reliable convergence, but cannot consider the rich data to data relations. This paper presents a new proxy based loss that takes advantages of both pair and proxy based methods and overcomes their limitations. Thanks to the use of proxies, our loss boosts the speed of convergence and is robust against noisy labels and outliers. At the same time, it allows embedding vectors of data to interact with each other in its gradients to exploit data to data relations. Our method is evaluated on four public benchmarks, where a standard network trained with our loss achieves state of the art performance and most quickly converges.


Conditions for Open Ended Evolution in Immigration Games. The Immigration Game (invented by Don Woods in 1971) extends the solitaire Game of Life (invented by John Conway in 1970) to enable two player competition. The Immigration Game can be used in a model of evolution by natural selection, where fitness is measured with competitions. The rules for the Game of Life belong to the family of semitotalistic rules, a family with 262,144 members. Woods' method for converting the Game of Life into a two player game generalizes to 8,192 members of the family of semitotalistic rules. In this paper, we call the original Immigration Game the Life Immigration Game and we call the 8,192 generalizations Immigration Games (including the Life Immigration Game). The question we examine here is, what are the conditions for one of the 8,192 Immigration Games to be suitable for modeling open ended evolution? Our focus here is specifically on conditions for the rules, as opposed to conditions for other aspects of the model of evolution. In previous work, it was conjectured that Turing completeness of the rules for the Game of Life may have been necessary for the success of evolution using the Life Immigration Game. Here we present evidence that Turing completeness is a sufficient condition on the rules of Immigration Games, but not a necessary condition. The evidence suggests that a necessary and sufficient condition on the rules of Immigration Games, for open ended evolution, is that the rules should allow growth.


A BK inequality for random matchings. Let G (S,T,E) be a bipartite graph. For a matching M of G , let V(M) be the set of vertices covered by M , and let B(M) be the symmetric difference of V(M) and S . We prove that if M is a uniform random matching of G , then B(M) satisfies the BK inequality for increasing events.


Sequences of well distributed vertices on graphs and spectral bounds on optimal transport. Given a graph G (V,E) , suppose we are interested in selecting a sequence of vertices (x j) j 1 n such that left x 1, dots, x k right is well distributed' uniformly in k . We describe a greedy algorithm motivated by potential theory and corresponding developments in the continuous setting. The algorithm performs nicely on graphs and may be of use for sampling problems. We can interpret the algorithm as trying to greedily minimize a negative Sobolev norm we explain why this is related to Wasserstein distance by establishing a purely spectral bound on the Wasserstein distance on graphs that mirrors R. Peyre's estimate in the continuous setting. We illustrate this with many examples and discuss several open problems.


Properties and applications of some special integer number sequences. In this paper, we provide properties and applications of some special integer sequences. We generalize and give some properties of Pisano period. Moreover, we provide a new application in Cryptography and applications of some quaternion elements.


Phylogenetic Study of 2019 nCoV by Using Alignment Free Method (Evolutionary Bifurcation of Novel Coronavirus Mutants). The phylogenetic tree of SARS CoV 2 (nCov 19) viruses is reconstructed according to the similarity of genome sequences. The tree topology of Betacoronavirus is remarkably consistent with biologist's systematics. Because the tree construction contains enough information about virus mutants, it is suitable to study the evolutionary relationship between novel coronavirus mutants transmitted among humans. The emergences of 14 kinds of main mutants are studied and these strains can be classified as eight bifurcations of the phylogenetic tree. It is found that there exist three types of virus mutations, namely, the mutation among sub branches of the same branch, the off root mutation and the root oriented mutation between large branches of the tree. From the point of the relation between viral mutation and host selection we found that individuals with low immunity provide a special environment for the positive natural selection of virus evolution. It gives a mechanism to explain why large mutations between two distant branches generally occur in the nCov 19 phylogenetic tree. The finding is helpful to formulate strategies to control the spread of COVID 19.


Assessing Demographic Bias in Named Entity Recognition. Named Entity Recognition (NER) is often the first step towards automated Knowledge Base (KB) generation from raw text. In this work, we assess the bias in various Named Entity Recognition (NER) systems for English across different demographic groups with synthetically generated corpora. Our analysis reveals that models perform better at identifying names from specific demographic groups across two datasets. We also identify that debiased embeddings do not help in resolving this issue. Finally, we observe that character based contextualized word representation models such as ELMo results in the least bias across demographics. Our work can shed light on potential biases in automated KB generation due to systematic exclusion of named entities belonging to certain demographics.


On Uniquely Registrable Networks. Consider a network with N nodes in d dimensional Euclidean space, and M subsets of these nodes P 1, cdots,P M . Assume that the nodes in a given P i are observed in a local coordinate system. The registration problem is to compute the coordinates of the N nodes in a global coordinate system, given the information about P 1, cdots,P M and the corresponding local coordinates. The network is said to be uniquely registrable if the global coordinates can be computed uniquely (modulo Euclidean transforms). We formulate a necessary and sufficient condition for a network to be uniquely registrable in terms of rigidity of the body graph of the network. A particularly simple characterization of unique registrability is obtained for planar networks. Further, we show that k vertex connectivity of the body graph is equivalent to quasi k connectivity of the bipartite correspondence graph of the network. Along with results from rigidity theory, this helps us resolve a recent conjecture due to Sanyal et al. (IEEE TSP, 2017) that quasi 3 connectivity of the correspondence graph is both necessary and sufficient for unique registrability in two dimensions. We present counterexamples demonstrating that while quasi (d 1) connectivity is necessary for unique registrability in any dimension, it fails to be sufficient in three and higher dimensions.


Signals induced on electrodes by moving charges, a general theorem for Maxwell's equations based on Lorentz reciprocity. We discuss a signal theorem for charged particle detectors where the finite propagation time of the electromagnetic waves produced by a moving charge cannot be neglected. While the original Ramo Shockley theorem and related extensions are all based on electrostatic or quasi electrostatic approximations, the theorem presented in this report is based on the full extent of Maxwell's equations and does account for all electrodynamic effects. It is therefore applicable to all devices that detect fields and radiation from charged particles.


Traveling wave solutions for a class of discrete diffusive SIR epidemic model. This paper is concerned with the conditions of existence and nonexistence of traveling wave solutions (TWS) for a class of discrete diffusive epidemic models. We find that the existence of TWS is determined by the so called basic reproduction number and the critical wave speed When the basic reproduction number R0 greater than 1, there exists a critical wave speed c 0, such that for each c c the system admits a nontrivial TWS and for c c there exists no nontrivial TWS for the system. In addition, the boundary asymptotic behaviour of TWS is obtained by constructing a suitable Lyapunov functional and employing Lebesgue dominated convergence theorem. Finally, we apply our results to two discrete diffusive epidemic models to verify the existence and nonexistence of TWS.


A new approach for the study of limit cycles. We prove that star like limit cycles of any planar polynomial system can also be seen either as solutions defined on a given interval of a new associated planar non autonomous polynomial system or as heteroclinic solutions of a 3 dimensional polynomial system. We illustrate these points of view with several examples. One of the key ideas in our approach is to decompose the periodic solutions as the sum of two suitable functions. As a first application we use these new approaches to prove that all star like reversible limit cycles are algebraic. As a second application we introduce a function whose zeroes control the periodic orbits that persist as limit cycles when we perturb a star like reversible center. As far as we know this is the first time that this question is solved in full generality. Somehow, this function plays a similar role that an Abelian integral for studying perturbations of Hamiltonian systems.


The impact of students' epistemological framing on a task requiring representational consistency. The ability to flexibly transform between different representations (e.g., from mathematical to graphical representations) of the same concept is a hallmark of expertise. Prior research suggests that many introductory physics students show lack of representational consistency, e.g., they may construct two representations of the same concept in the same situation that are inconsistent with one another. In this case study, we asked students to construct two representations for the electric field for a situation involving Gauss's law with spherical symmetry (charge conducting sphere surrounded by charged conducting spherical shell). Prior research also suggests that this type of problem results in many students constructing representations that are not consistent with one another. Here we present findings from individual interviews with three students about this problem which suggest that students' lack of representational consistency may be partly attributed to the type of knowledge that the graphical and mathematical representations trigger. In the epistemic games framework terminology, the two representations students are asked to construct (mathematical vs. graphical) in the problem may contribute to their lack of representational consistency.


Matrix Valued Mean Field Type Games Risk Sensitive, Adversarial, and Risk Neutral Linear Quadratic Case. In this paper we study a class of matrix valued linear quadratic mean field type games for both the risk neutral, risk sensitive and robust cases. Non cooperation, full cooperation and adversarial between teams are treated. We provide a semi explicit solution for both problems by means of a direct method. The state dynamics is described by a matrix valued linear jump diffusion regime switching system of conditional mean field type where the conditioning is with respect to common noise which is a regime switching process. The optimal strategies are in a state and conditional mean field feedback form. Semi explicit solutions of equilibrium costs and strategies are also provided for the full cooperative, adversarial teams, risk sensitive full cooperative and risk sensitive adversarial team cases. It is shown that full cooperation increases the well posedness domain under risk sensitive decision makers by means of population risk sharing. Finally, relationship between risk sensitivity and robustness are established in the mean field type context.


End to End Speech Translation with Knowledge Distillation. End to end speech translation (ST), which directly translates from source language speech into target language text, has attracted intensive attentions in recent years. Compared to conventional pipeline systems, end to end ST models have advantages of lower latency, smaller model size and less error propagation. However, the combination of speech recognition and text translation in one model is more difficult than each of these two tasks. In this paper, we propose a knowledge distillation approach to improve ST model by transferring the knowledge from text translation model. Specifically, we first train a text translation model, regarded as a teacher model, and then ST model is trained to learn output probabilities from teacher model through knowledge distillation. Experiments on English French Augmented LibriSpeech and English Chinese TED corpus show that end to end ST is possible to implement on both similar and dissimilar language pairs. In addition, with the instruction of teacher model, end to end ST model can gain significant improvements by over 3.5 BLEU points.


Comparative analysis of word embeddings in assessing semantic similarity of complex sentences. Semantic textual similarity is one of the open research challenges in the field of Natural Language Processing. Extensive research has been carried out in this field and near perfect results are achieved by recent transformer based models in existing benchmark datasets like the STS dataset and the SICK dataset. In this paper, we study the sentences in these datasets and analyze the sensitivity of various word embeddings with respect to the complexity of the sentences. We build a complex sentences dataset comprising of 50 sentence pairs with associated semantic similarity values provided by 15 human annotators. Readability analysis is performed to highlight the increase in complexity of the sentences in the existing benchmark datasets and those in the proposed dataset. Further, we perform a comparative analysis of the performance of various word embeddings and language models on the existing benchmark datasets and the proposed dataset. The results show the increase in complexity of the sentences has a significant impact on the performance of the embedding models resulting in a 10 20 decrease in Pearson's and Spearman's correlation.


Towards Robust Explanations for Deep Neural Networks. Explanation methods shed light on the decision process of black box classifiers such as deep neural networks. But their usefulness can be compromised because they are susceptible to manipulations. With this work, we aim to enhance the resilience of explanations. We develop a unified theoretical framework for deriving bounds on the maximal manipulability of a model. Based on these theoretical insights, we present three different techniques to boost robustness against manipulation training with weight decay, smoothing activation functions, and minimizing the Hessian of the network. Our experimental results confirm the effectiveness of these approaches.


Representations and deformations of Hom Lie Yamaguti superalgebras. Let (L, alpha) be a Hom Lie Yamaguti superalgebra. We first introduce the representation and cohomology theory of Hom Lie Yamaguti superalgebras. Furthermore, we introduce the notions of generalized derivations and representations of (L, alpha) and present some properties. Finally, we investigate the deformations of (L, alpha) by choosing some suitable cohomology.


Littlewood Paley Stein inequalities on textup RCD (K, infty) spaces. The L p boundedness on vertical Littlewood Paley square functions for heat flows on textup RCD (K, infty) spaces with K in mathbb R is proved. With regards to the proof, for 1 p leq 2 , Stein's analytical method is applied, while for 2 p infty , the probabilistic approach in the sense of Ba n uelos Bogdan Luks introduced recently is employed.


Lorentzian polynomials from polytope projections. Lorentzian polynomials, recently introduced by Br "and 'en and Huh, generalize the notion of log concavity of sequences to homogeneous polynomials whose supports are integer points of generalized permutahedra. Br "and 'en and Huh show that normalizations of polynomials equaling integer point transforms of generalized permutahedra are Lorentzian moreover, normalizations of certain projections of integer point transforms of generalized permutahedra with zero one vertices are also Lorentzian. Taking this polytopal perspective further, we show that normalizations of certain projections of integer point transforms of flow polytopes (which, before projection, are not Lorentzian), are also Lorentzian.


TD implies CCR. Assuming mathrm ZF , we prove that Turing determinacy ( mathrm TD ) implies countable choice axiom for sets of reals ( mathrm CCR ).


Sample driven optimal stopping From the secretary problem to the i.i.d. prophet inequality. We take a unifying approach to single selection optimal stopping problems with random arrival order and independent sampling of items. In the problem we consider, a decision maker (DM) initially gets to sample each of N items independently with probability p , and can observe the relative rankings of these sampled items. Then, the DM faces the remaining items in an online fashion, observing the relative rankings of all revealed items. While scanning the sequence the DM makes irrevocable stop continue decisions and her reward for stopping the sequence facing the item with rank i is Y i . The goal of the DM is to maximize her reward. We start by studying the case in which the values Y i are known to the DM, and then move to the case in which these values are adversarial. For the former case, we write the natural linear program that captures the performance of an algorithm, and take its continuous limit. We prove a structural result about this continuous limit, which allows us to reduce the problem to a relatively simple real optimization problem. We establish that the optimal algorithm is given by a sequence of thresholds t 1 le t 2 le cdots such that the DM should stop if seeing an item with current ranking i after time t i . Additionally we are able to recover several classic results in the area such as those for secretary problem and the minimum ranking problem. For the adversarial case, we obtain a similar linear program with an additional stochastic dominance constraint. Using the same machinery we are able to pin down the optimal competitive ratios for all values of p . Notably, we prove that as p approaches 1, our guarantee converges linearly to 0.745, matching that of the i.i.d. prophet inequality. Also interesting is the case p 1 2 , where our bound evaluates to 0.671 , which improves upon the state of the art.


Interdicting Restructuring Networks with Applications in Illicit Trafficking. We consider a new class of max flow network interdiction problems, where the defender is able to introduce new arcs to the network after the attacker has made their interdiction decisions. We prove properties of when this restructuring will not increase the value of the minimum cut, which has important practical interpretations for problems of disrupting drug trafficking networks. In particular, it demonstrates that disrupting lower levels of these networks will not impact their operations when replacing the disrupted participants is easy. For the bilevel mixed integer linear programming formulation of this problem, we devise a column and constraint generation (C CG) algorithm to solve it. Our approach uses partial information on the feasibility of restructuring plans and is shown to be orders of magnitude faster than previous C CG methods. We demonstrate that applying decisions from standard max flow network interdiction problems can result in significantly higher flows than interdictions that account for the restructuring.


A degree sequence version of the K "uhn Osthus tiling theorem. A fundamental result of K "uhn and Osthus The minimum degree threshold for perfect graph packings, Combinatorica, 2009 determines up to an additive constant the minimum degree threshold that forces a graph to contain a perfect H tiling. We prove a degree sequence version of this result which allows for a significant number of vertices to have lower degree.


MAD VAE Manifold Awareness Defense Variational Autoencoder. Although deep generative models such as Defense GAN and Defense VAE have made significant progress in terms of adversarial defenses of image classification neural networks, several methods have been found to circumvent these defenses. Based on Defense VAE, in our research we introduce several methods to improve the robustness of defense models. The methods introduced in this paper are straight forward yet show promise over the vanilla Defense VAE. With extensive experiments on MNIST data set, we have demonstrated the effectiveness of our algorithms against different attacks. Our experiments also include attacks on the latent space of the defensive model. We also discuss the applicability of existing adversarial latent space attacks as they may have a significant flaw.


Metrics Towards Measuring Cyber Agility. In cyberspace, evolutionary strategies are commonly used by both attackers and defenders. For example, an attacker's strategy often changes over the course of time, as new vulnerabilities are discovered and or mitigated. Similarly, a defender's strategy changes over time. These changes may or may not be in direct response to a change in the opponent's strategy. In any case, it is important to have a set of quantitative metrics to characterize and understand the effectiveness of attackers' and defenders' evolutionary strategies, which reflect their em cyber agility . Despite its clear importance, few systematic metrics have been developed to quantify the cyber agility of attackers and defenders. In this paper, we propose the first metric framework for measuring cyber agility in terms of the effectiveness of the dynamic evolution of cyber attacks and defenses. The proposed framework is generic and applicable to transform any relevant, quantitative, and or conventional static security metrics (e.g., false positives and false negatives) into dynamic metrics to capture dynamics of system behaviors. In order to validate the usefulness of the proposed framework, we conduct case studies on measuring the evolution of cyber attacks and defenses using two real world datasets. We discuss the limitations of the current work and identify future research directions.


Optimal execution with rough path signatures. We present a method for obtaining approximate solutions to the problem of optimal execution, based on a signature method. The framework is general, only requiring that the price process is a geometric rough path and the price impact function is a continuous function of the trading speed. Following an approximation of the optimisation problem, we are able to calculate an optimal solution for the trading speed in the space of linear functions on a truncation of the signature of the price process. We provide strong numerical evidence illustrating the accuracy and flexibility of the approach. Our numerical investigation both examines cases where exact solutions are known, demonstrating that the method accurately approximates these solutions, and models where exact solutions are not known. In the latter case, we obtain favourable comparisons with standard execution strategies.


Arithmetic topology in Ihara theory II Milnor invariants, dilogarithmic Heisenberg coverings and triple power residue symbols. We introduce mod l Milnor invariants of a Galois element associated to Ihara's Galois representation on the pro l fundamental group of a punctured projective line ( l being a prime number), as arithmetic analogues of Milnor invariants of a pure braid. We then show that triple quadratic (resp. cubic) residue symbols of primes in the rational (resp. Eisenstein) number field are expressed by mod 2 (resp. mod 3 ) triple Milnor invariants of Frobenius elements. For this, we introduce dilogarithmic mod l Heisenberg ramified covering cal D (l) of mathbb P 1 , which may be regarded as a higher analog of the dilogarithmic function, for the gerbe associated to the mod l Heisenberg group, and we study the monodromy transformations of certain functions on cal D (l) along the pro l longitudes of Frobenius elements for l 2,3 .


Drifts and Shifts Characterizing the Evolution of Users Interests on Reddit. Selective exposure is the main driver for the economy of attention when consuming online content. We select information adhering to our system of beliefs and ignore dissenting information. However, even personal interest is likely to play a role in determining our attention patterns. To understand in more detail the dynamics of interest driven choices, we address the evolution of users' interest on Reddit by means of an analysis on more than 7 million of users. We observe that the lifetime distribution of users on subreddits (online thematic communities) is 'short' with respect to the observation period. Furthermore, users tend to be active on a very limited number of subreddits with respect to the wide offer of the platform. These aspects indicate the presence of a migrating behavior of users among subreddits. To characterize this phenomenon we propose a metric based on a geometrical encoding of the 'interest space' of the user. The movement of users across subreddits is characterized by a bursty trend, made of sudden variations. The most frequent of them take place between recreational subreddits and those more related to news and politics. We describe this kind of activity with two behaviors called interest drift and interest shift. Our results suggest that selective exposure and personal interest coexist in driving content consumption.


FA Harris A Fast and Asynchronous Corner Detector for Event Cameras. Recently, the emerging bio inspired event cameras have demonstrated potentials for a wide range of robotic applications in dynamic environments. In this paper, we propose a novel fast and asynchronous event based corner detection method which is called FA Harris. FA Harris consists of several components, including an event filter, a Global Surface of Active Events (G SAE) maintaining unit, a corner candidate selecting unit, and a corner candidate refining unit. The proposed G SAE maintenance algorithm and corner candidate selection algorithm greatly enhance the real time performance for corner detection, while the corner candidate refinement algorithm maintains the accuracy of performance by using an improved event based Harris detector. Additionally, FA Harris does not require artificially synthesized event frames and can operate on asynchronous events directly. We implement the proposed method in C and evaluate it on public Event Camera Datasets. The results show that our method achieves approximately 8x speed up when compared with previously reported event based Harris detector, and with no compromise on the accuracy of performance.


2D Anisotropic KPZ at stationarity scaling, tightness and non triviality. In this work we focus on the two dimensional anisotropic KPZ (aKPZ) equation, which is formally given by begin equation partial t h frac nu 2 Delta h lambda(( partial 1 h) 2 ( partial 2 h) 2) nu frac 1 2 xi, end equation where xi denotes a noise which is white in both space and time, and lambda and nu are positive constants. Due to the wild oscillations of the noise and the quadratic nonlinearity, the previous equation is classically ill posed. It is not possible to linearise it via the Cole Hopf transformation and the pathwise techniques for singular SPDEs (the theory of Regularity Structures by M. Hairer or the paracontrolled distributions approach of M. Gubinelli, P. Imkeller, N. Perkowski) are not applicable. In the present work, we consider a regularised version of aKPZ which preserves its invariant measure. We show that in order to have subsequential limits once the regularisation is removed, it is necessary to suitably renormalise lambda and nu . Moreover, we prove that, in the regime suggested by the (non rigorous) renormalisation group computations of D.E. Wolf, "Kinetic roughening of vicinal surfaces'', Phys. Rev. Lett., 1991 , i.e. nu constant and the coupling constant lambda converging to 0 as the inverse of the square root logarithm, any limit differs from the solution to the linear equation obtained by simply dropping the nonlinearity in aKPZ.


A mathematical model for cell polarization in zebrafish primordial germ cells. Blebs are cell protrusions generated by local membrane cortex detachments followed by expansion of the plasma membrane. Blebs are formed by some migrating cells, for example primordial germ cells of the zebrafish. While blebs occur randomly at each part of the membrane in unpolarized cells, a polarization process guarantees the occurrence of blebs at a preferential site and thereby facilitates migration towards a specified direction. Little is known about the factors involved in development and maintenance of a polarized state, yet recent studies revealed the influence of an intracellular flow and the stabilizing role of the membrane cortex linker molecule Ezrin. Based on this information, we develop and analyse a coupled bulk surface model describing a potential cellular mechanism by which a bleb could be induced at a controlled site. The model rests upon intracellular Darcy flow and a diffusion advection reaction system, describing the temporal evolution from an unpolarized to a stable polarized Ezrin distribution. We prove the well posedness of the mathematical model and show that simulations qualitatively correspond to experimental observations, suggesting that indeed the interaction of an intracellular flow with membrane proteins can be the cause of the cell polarization.