Enhancing Video Streaming in Vehicular Networks via Resource Slicing. Vehicle to everything (V2X) communication is a key enabler that connects vehicles to neighboring vehicles, infrastructure and pedestrians. In the past few years, multimedia services have seen an enormous growth and it is expected to increase as more devices will utilize infotainment services in the future i.e. vehicular devices. Therefore, it is important to focus on user centric measures i.e. quality of experience (QoE) such as video quality (resolution) and fluctuations therein. In this paper, a novel joint video quality selection and resource allocation technique is proposed for increasing the QoE of vehicular devices. The proposed approach exploits the queuing dynamics and channel states of vehicular devices, to maximize the QoE while ensuring seamless video playback at the end users with high probability. The network wide QoE maximization problem is decoupled into two subparts. First, a network slicing based clustering algorithm is applied to partition the vehicles into multiple logical networks. Secondly, vehicle scheduling and quality selection is formulated as a stochastic optimization problem which is solved using the Lyapunov drift plus penalty method. Numerical results show that the proposed algorithm ensures high video quality experience compared to the baseline. Simulation results also show that the proposed technique achieves low latency and high reliability communication.


Artist and style exposure bias in collaborative filtering based music recommendations. Algorithms have an increasing influence on the music that we consume and understanding their behavior is fundamental to make sure they give a fair exposure to all artists across different styles. In this on going work we contribute to this research direction analyzing the impact of collaborative filtering recommendations from the perspective of artist and music style exposure given by the system. We first analyze the distribution of the recommendations considering the exposure of different styles or genres and compare it to the users' listening behavior. This comparison suggests that the system is reinforcing the popularity of the items. Then, we simulate the effect of the system in the long term with a feedback loop. From this simulation we can see how the system gives less opportunity to the majority of artists, concentrating the users on fewer items. The results of our analysis demonstrate the need for a better evaluation methodology for current music recommendation algorithms, not only limited to user focused relevance metrics.


p adic families of d th Shintani liftings. In this note we give a detailed construction of a Lambda adic mathfrak d th Shintani lifting. We derive a Lambda adic version of Kohnen's formula relating Fourier coefficients of half integral weight modular forms and special values of twisted L series. As a by product, we derive a mild generalization of such classical formulae, and also point out a relation between Fourier coefficients of Lambda adic d th Shintani liftings and Stark Heegner points.


Regular Model Checking Revisited (Technical Report). In this contribution we revisit regular model checking, a powerful framework that has been successfully applied for the verification of infinite state systems, especially parameterized systems (concurrent systems with an arbitrary number of processes). We provide a reformulation of regular model checking with length preserving transducers in terms of existential second order theory over automatic structures. We argue that this is a natural formulation that enables us tap into powerful synthesis techniques that have been extensively studied in the software verification community. More precisely, in this formulation the first order part represents the verification conditions for the desired correctness property (for which we have complete solvers), whereas the existentially quantified second order variables represent the relations to be synthesized. We show that many interesting correctness properties can be formulated in this way, examples being safety, liveness, bisimilarity, and games. More importantly, we show that this new formulation allows new interesting benchmarks (and old regular model checking benchmarks that were previously believed to be difficult), especially in the domain of parameterized system verification, to be solved.


New Steiner systems from old ones by paramodifications. Techniques of producing new combinatorial structures from old ones are commonly called trades. The switching principle applies for a broad class of designs it is a local transformation that modifies two columns of the incidence matrix. In this paper, we present a construction, which is a generalization of the switching transform for the class of Steiner 2 designs. We call this construction paramodification of Steiner 2 designs, since it modifies the parallelism of a subsystem. We study in more detail the paramodifications of affine planes, Steiner triple systems, and abstract unitals. Computational results show that paramodification can construct many new unitals.


Modeling water relative permeability in uniform grain packs. Accurate estimation of water relative permeability has been of great interest in various research areas because of its broad applications in soil physics and hydrology as well as oil and gas production and recovery. Critical path analysis (CPA), a promising technique from statistical physics, is well known to be applicable to heterogeneous media with broad conductance or pore size distribution (PSD). By heterogeneity, we mean variations in the geometrical properties of pore space. In this study, we demonstrate that CPA is also applicable to packings of spheres of the same size, known as homogeneous porous media. More specifically, we apply CPA to model water relative permeability (krw) in mono sized sphere packs whose PSDs are fairly narrow. We estimate the krw from (1) the PSD and (2) the PSD and saturation dependent electrical conductivity ( sigma r) for both drainage and imbibition processes. We show that the PSD of mono sized sphere packs approximately follows the log normal probability density function. Comparison with numerical simulations indicate that both the imbibition and drainage krw are estimated from the PSD and sigma r data more accurately than those from the PSD. We show that CPA can estimate krw in mono sized sphere packs precisely.


Bounds for the Petersson norms of the pullbacks of Saito Kurokawa lifts. Using the amplification technique, we prove that mass' of the pullback of the Saito Kurokawa lift of a Hecke eigen form g in S 2k is bounded by k 1 frac 1 210 epsilon . This improves the previously known bound k for this quantity.


The Wellposedness of Path dependent Multidimensional Forward backward SDE. We study in this paper the wellposedness of path dependent multidimensional forward backward stochastic differential equations (FBSDE). By path dependent we mean that the coefficients of the forward backward SDE at time t can depend on the whole path of the forward process up to time t. These kinds of forward backward SDE appear when solving path dependent stochastic control problem by means of variational calculus. At the heart of our analysis is the construction of a decoupling random field on the path space. We first prove the existence and the uniqueness of decoupling field on small time interval. Then by introducing the characteristic BSDE, we show that a global decoupling field can be constructed by patching local solutions together as long as the solution of the characteristic BSDE remains bounded. Finally, we show that the solution of a path dependent forward backward SDE is stable.


Computational Experimental Investigation of a Fission Thermal Probe in TREAT. The development of nuclear fuels requires unsteady transient testing for design process and qualification under postulated accident conditions. Breach, rupture, fracture, melting, and other fuel failure modes may occur during the use of nuclear fuels or when they are exposed to extreme overpower conditions. Therefore, it is crucial for the fuel to retain reasonable structural integrity and coolable geometry. The experimental facility of Transient Reactor Test (TREAT) was designed and constructed in November 1958 to conduct transient testing of fuels and structural materials. The magnitude of nuclear heating is one of the most important key parameters in test design, analysis, and data interpretation in TREAT. Some steady state tests are able to measure heating directly via enthalpy rise of coolant in thermally isolated loops, but an enormous number of important tests must rely upon nuclear modeling to correlate core operating parameters to specimen nuclear heating. Uncertainties of these models and their inputs can prevent experimenters from achieving desired conditions and hamper advanced models from describing the phenomena to the level of confidence needed. The situation can be even more difficult for unsteady transient tests where nuclear heating is intentionally varied over short time scales. This research develops a novel nuclear heating sensor technology which directly measures the parameters of interest using spatially resolved real time thermometry of fissionable instrument materials, demonstrate the instruments' use in TREAT, and perform data comparisons to nuclear models to facilitate an advanced understanding of transient power coupling. Here we present our results of the designed probe and its testing outcome.


Rapid measurement of tension in multi wire arrays using free damped oscillations. We explore the possibility of using free (as opposed to driven) oscillations to verify wire tension in large, open frame multiwire planes. Using finite element simulation we predict the signal when two wires are set in motion by applying and holding a DC voltage between them (electrostatic excitation). The resulting current signal is marginal. However, the program relates voltage, oscillation amplitude and signal, and it is easy to get enough amplitude (of order microns) by gently tapping the wire frame. We perform an end to end experimental test. Four 1.5m wires, at various tensions around 5N, are mounted in an aluminum channel. They are connected to a low noise, balanced current to voltage preamp designed for good common mode rejection. The 150V bias, from 9V batteries in series, is built into the preamp. Data are taken with a Tektronix 1202B scope, stored on a USB flash drive, and transferred to a laptop computer for analysis. Ten measurements, each with 2500 samples spanning 1s, are taken for each of the three wire pairs. The wires are excited by tapping the frame with a small mallet. Because the wire tensions are different, the waveforms exhibit beats. The data are analyzed by a program which decomposes the waveform step by step to extract the base and beat frequencies and thereby, two resonant frequencies per wire pair. Finally, those frequencies are disambiguated to obtain the frequency for each wire, which is then converted to a tension. The four tensions are then found mechanically by measuring the g mm required to deflect each wire at its midpoint. The two values of tension agree well. We conclude that the method verifies tension to approximately 2 (1 sigma) in a little over 1s for 1.5m wires, and would be faster for shorter wires.


Unified Gas kinetic Wave Particle Methods II Multiscale Simulation on Unstructured Mesh. In this paper, we present a unified gas kinetic wave particle (UGKWP) method on unstructured mesh for multiscale simulation of continuum and rarefied flow. Inheriting from the multicale transport in the unified gas kinetic scheme (UGKS), the integral solution of kinetic model equation is employed in the construction of UGKWP method to model the flow physics in the cell size and time step scales. A novel wave particle adaptive formulation is introduced in the UGKWP method to describe the flow dynamics in each control volume. The local gas evolution is constructed through the dynamical interaction of the deterministic hydrodynamic wave and the stochastic kinetic particle. Within the resolution of cell size and time step, the decomposition, interaction, and evolution of the hydrodynamic wave and the kinetic particle depend on the ratio of the time step to the local particle collision time. In the rarefied flow regime, the flow physics is mainly recovered by the discrete particles and the UGKWP method performs as a stochastic particle method. In the continuum flow regime, the flow behavior is solely followed by macroscopic variable evolution and the UGKWP method becomes a gas kinetic hydrodynamic flow solver for the viscous and heat conducting Navier Stokes solutions. In different flow regimes, many numerical test cases are computed to validate the UGKWP method on unstructured mesh. The UGKWP method can get the same UGKS solutions in all Knudsen regimes without the requirement of the time step and mesh size being less than than the particle collision time and mean free path. With an automatic wave particle decomposition, the UGKWP method becomes very efficient. For example, at Mach number 30 and Knudsen number 0.1, in comparison with UGKS several order of magnitude reductions in computational cost and memory requirement have been achieved by UGKWP.


Inertial Three Operator Splitting Method and Applications. We introduce an inertial variant of the forward Douglas Rachford splitting and analyze its convergence. We specify an instance of the proposed method to the three composite convex minimization template. We provide practical guidance on the selection of the inertial parameter based on the adaptive starting idea. Finally, we illustrate the practical performance of our method in various machine learning applications.


A neural network closure for the Euler Poisson system based on kinetic simulations. This work deals with the modeling of plasmas, which are charged particle fluids. Thanks to machine leaning, we construct a closure for the one dimensional Euler Poisson system valid for a wide range of collision regimes. This closure, based on a fully convolutional neural network called V net, takes as input the whole spatial density, mean velocity and temperature and predicts as output the whole heat flux. It is learned from data coming from kinetic simulations of the Vlasov Poisson equations. Data generation and preprocessings are designed to ensure an almost uniform accuracy over the chosen range of Knudsen numbers (which parametrize collision regimes). Finally, several numerical tests are carried out to assess validity and flexibility of the whole pipeline.


A Deep Learning Pipeline for Patient Diagnosis Prediction Using Electronic Health Records. Augmentation of disease diagnosis and decision making in healthcare with machine learning algorithms is gaining much impetus in recent years. In particular, in the current epidemiological situation caused by COVID 19 pandemic, swift and accurate prediction of disease diagnosis with machine learning algorithms could facilitate identification and care of vulnerable clusters of population, such as those having multi morbidity conditions. In order to build a useful disease diagnosis prediction system, advancement in both data representation and development of machine learning architectures are imperative. First, with respect to data collection and representation, we face severe problems due to multitude of formats and lack of coherency prevalent in Electronic Health Records (EHRs). This causes hindrance in extraction of valuable information contained in EHRs. Currently, no universal global data standard has been established. As a useful solution, we develop and publish a Python package to transform public health dataset into an easy to access universal format. This data transformation to an international health data format facilitates researchers to easily combine EHR datasets with clinical datasets of diverse formats. Second, machine learning algorithms that predict multiple disease diagnosis categories simultaneously remain underdeveloped. We propose two novel model architectures in this regard. First, DeepObserver, which uses structured numerical data to predict the diagnosis categories and second, ClinicalBERT Multi, that incorporates rich information available in clinical notes via natural language processing methods and also provides interpretable visualizations to medical practitioners. We show that both models can predict multiple diagnoses simultaneously with high accuracy.


NASS Optimizing Secure Inference via Neural Architecture Search. Due to increasing privacy concerns, neural network (NN) based secure inference (SI) schemes that simultaneously hide the client inputs and server models attract major research interests. While existing works focused on developing secure protocols for NN based SI, in this work, we take a different approach. We propose NASS, an integrated framework to search for tailored NN architectures designed specifically for SI. In particular, we propose to model cryptographic protocols as design elements with associated reward functions. The characterized models are then adopted in a joint optimization with predicted hyperparameters in identifying the best NN architectures that balance prediction accuracy and execution efficiency. In the experiment, it is demonstrated that we can achieve the best of both worlds by using NASS, where the prediction accuracy can be improved from 81.6 to 84.6 , while the inference runtime is reduced by 2x and communication bandwidth by 1.9x on the CIFAR 10 dataset.


On the Importance of Delexicalization for Fact Verification. In this work we aim to understand and estimate the importance that a neural network assigns to various aspects of the data while learning and making predictions. Here we focus on the recognizing textual entailment (RTE) task and its application to fact verification. In this context, the contributions of this work are as follows. We investigate the attention weights a state of the art RTE method assigns to input tokens in the RTE component of fact verification systems, and confirm that most of the weight is assigned to POS tags of nouns (e.g., NN, NNP etc.) or their phrases. To verify that these lexicalized models transfer poorly, we implement a domain transfer experiment where a RTE component is trained on the FEVER data, and tested on the Fake News Challenge (FNC) dataset. As expected, even though this method achieves high accuracy when evaluated in the same domain, the performance in the target domain is poor, marginally above chance.To mitigate this dependence on lexicalized information, we experiment with several strategies for masking out names by replacing them with their semantic category, coupled with a unique identifier to mark that the same or new entities are referenced between claim and evidence. The results show that, while the performance on the FEVER dataset remains at par with that of the model trained on lexicalized data, it improves significantly when tested in the FNC dataset. Thus our experiments demonstrate that our strategy is successful in mitigating the dependency on lexical information.


Doubling nodal solutions to the Yamabe equation in mathbb R n with maximal rank. We construct a new family of entire solutions to the Yamabe equation Delta u frac n(n 2) 4 u frac 4 n 2 u mbox in mathcal D 1,2 ( mathbb R n). If n 3 , our solutions have maximal rank, being the first example in odd dimension. Our construction has analogies with the doubling of the equatorial spheres in the construction of minimal surfaces in S 3(1) .


Multi patch epidemic models with general exposed and infectious periods. We study multi patch epidemic models where individuals may migrate from one patch to another in either of the susceptible, exposed latent, infectious and recovered states. We assume that infections occur both locally with a rate that depends on the patch as well as "from distance" from all the other patches. The exposed and infectious periods have general distributions, and are not affected by the possible migrations of the individuals. The migration processes in either of the three states are assumed to be Markovian, and independent of the exposed and infectious periods. We establish a functional law of large number (FLLN) and a function central limit theorem (FCLT) for the susceptible, exposed latent, infectious and recovered processes. In the FLLN, the limit is determined by a set of Volterra integral equations. In the special case of deterministic exposed and infectious periods, the limit becomes a system of ODEs with delays. In the FCLT, the limit is given by a set of stochastic Volterra integral equations driven by a sum of independent Brownian motions and continuous Gaussian processes with an explicit covariance structure.


Contrastive Learning for Sequential Recommendation. Sequential recommendation methods play a crucial role in modern recommender systems because of their ability to capture a user's dynamic interest from her his historical interactions. Despite their success, we argue that these approaches usually rely on the sequential prediction task to optimize the huge amounts of parameters. They usually suffer from the data sparsity problem, which makes it difficult for them to learn high quality user representations. To tackle that, inspired by recent advances of contrastive learning techniques in the computer version, we propose a novel multi task model called textbf C ontrastive textbf L earning for textbf S equential textbf Rec ommendation ( textbf CL4SRec ). CL4SRec not only takes advantage of the traditional next item prediction task but also utilizes the contrastive learning framework to derive self supervision signals from the original user behavior sequences. Therefore, it can extract more meaningful user patterns and further encode the user representation effectively. In addition, we propose three data augmentation approaches to construct self supervision signals. Extensive experiments on four public datasets demonstrate that CL4SRec achieves state of the art performance over existing baselines by inferring better user representations.


Biclustering and Boolean Matrix Factorization in Data Streams. We study the clustering of bipartite graphs and Boolean matrix factorization in data streams. We consider a streaming setting in which the vertices from the left side of the graph arrive one by one together with all of their incident edges. We provide an algorithm that, after one pass over the stream, recovers the set of clusters on the right side of the graph using sublinear space to the best of our knowledge, this is the first algorithm with this property. We also show that after a second pass over the stream, the left clusters of the bipartite graph can be recovered and we show how to extend our algorithm to solve the Boolean matrix factorization problem (by exploiting the correspondence of Boolean matrices and bipartite graphs). We evaluate an implementation of the algorithm on synthetic data and on real world data. On real world datasets the algorithm is orders of magnitudes faster than a static baseline algorithm while providing quality results within a factor 2 of the baseline algorithm. Our algorithm scales linearly in the number of edges in the graph. Finally, we analyze the algorithm theoretically and provide sufficient conditions under which the algorithm recovers a set of planted clusters under a standard random graph model.


SoK Blockchain Technology and Its Potential Use Cases. Bitcoin's success has led to significant interest in its underlying components, particularly Blockchain technology. Over 10 years after Bitcoin's initial release, the community still suffers from a lack of clarity regarding what properties defines Blockchain technology, its relationship to similar technologies, and which of its proposed use cases are tenable and which are little more than hype. In this paper we answer four common questions regarding Blockchain technology (1) what exactly is Blockchain technology, (2) what capabilities does it provide, and (3) what are good applications for Blockchain technology, and (4) how does it relate to other approache distributed technologies (e.g., distributed databases). We accomplish this goal by using grounded theory (a structured approach to gathering and analyzing qualitative data) to thoroughly analyze a large corpus of literature on Blockchain technology. This method enables us to answer the above questions while limiting researcher bias, separating thought leadership from peddled hype and identifying open research questions related to Blockchain technology. The audience for this paper is broad as it aims to help researchers in a variety of areas come to a better understanding of Blockchain technology and identify whether it may be of use in their own research.


On the weak Lefschetz property for almost complete intersections generated by uniform powers of general linear forms. In 2012, Migliore, the first author, and Nagel conjectured that, for all n geq 4 , the artinian ideal I (L 0 d, ldots,L 2n 1 d) subset R k x 0, ldots,x 2n generated by the d th powers of 2n 2 general linear forms fails to have the weak Lefschetz property if and only if d 1 . This paper is entirely devoted to prove partially this conjecture. More precisely, we prove that R I fails to have the weak Lefschetz property, provided 4 leq n leq 8, d geq 4 or d 2r, 1 leq r leq 8, 4 leq n leq 2r(r 2) 1 .


CDPM Convolutional Deformable Part Models for Semantically Aligned Person Re identification. Part level representations are essential for robust person re identification. However, common errors that arise during pedestrian detection frequently result in severe misalignment problems for body parts, which degrade the quality of part representations. Accordingly, to deal with this problem, we propose a novel model named Convolutional Deformable Part Models (CDPM). CDPM works by decoupling the complex part alignment procedure into two easier steps first, a vertical alignment step detects each body part in the vertical direction, with the help of a multi task learning model second, a horizontal refinement step based on attention suppresses the background information around each detected body part. Since these two steps are performed orthogonally and sequentially, the difficulty of part alignment is significantly reduced. In the testing stage, CDPM is able to accurately align flexible body parts without any need for outside information. Extensive experimental results demonstrate the effectiveness of the proposed CDPM for part alignment. Most impressively, CDPM achieves state of the art performance on three large scale datasets Market 1501, DukeMTMC ReID,and CUHK03.


Counting critical subgraphs in k critical graphs. Gallai asked in 1984 if any k critical graph on n vertices contains at least n distinct (k 1) critical subgraphs. The answer is trivial for k leq 3 . Improving a result of Stiebitz, Abbott and Zhou proved in 1995 that for all k geq 4 , such graph contains Omega(n 1 (k 1) ) distinct (k 1) critical subgraphs. Since then no progress had been made until very recently, Hare resolved the case k 4 by showing that any 4 critical graph on n vertices contains at least (8n 29) 3 odd cycles. In this paper, we mainly focus on 4 critical graphs and develop some novel tools for counting cycles of specified parity. Our main result shows that any 4 critical graph on n vertices contains Omega(n 2) odd cycles, which is tight up to a constant factor by infinite many graphs. As a crucial step, we prove the same bound for 3 connected non bipartite graphs, which may be of independent interest. Using the tools, we also give a very short proof for the case k 4 . Moreover, we improve the longstanding lower bound of Abbott and Zhou to Omega(n 1 (k 2) ) for the general case k geq 5 . We will also discuss some related problems on k critical graphs in the final section.


A Least Squares Finite Element Reduced Basis Method. We present a reduced basis (RB) method for parametrized linear elliptic partial differential equations (PDEs) in a least squares finite element framework. A rigorous and reliable error estimate is developed, and is shown to bound the error with respect to the exact solution of the PDE, in contrast to estimates that measure error with respect to a finite dimensional (high fidelity) approximation. It is shown that the first order formulation of the least squares finite element is a key ingredient. The method is demonstrated using numerical examples.


Error and Tamper Tolerant State Estimation for Discrete Event Systems under Cost Constraints. This paper deals with the state estimation problem in discrete event systems modeled with nondeterministic finite automata, partially observed via a sensor measuring unit whose measurements (reported observations) may be vitiated by a malicious attacker. The attacks considered in this paper include arbitrary deletions, insertions, or substitutions of observed symbols by taking into account a bounded number of attacks or, more generally, a total cost constraint (assuming that each deletion, insertion, or substitution bears a positive cost to the attacker). An efficient approach is proposed to describe possible sequences of observations that match the one received by the measuring unit, as well as their corresponding state estimates and associated total costs. We develop an algorithm to obtain the least cost matching sequences by reconstructing only a finite number of possible sequences, which we subsequently use to efficiently perform state estimation. We also develop a technique for verifying tamper tolerant diagnosability under attacks that involve a bounded number of deletions, insertions, and substitutions (or, more generally, under attacks of bounded total cost) by using a novel structure obtained by attaching attacks and costs to the original plant. The overall construction and verification procedure have complexity that is of O( X 2 C 2),where X is the number of states of the given finite automaton and C is the maximum total cost that is allowed for all the deletions, insertions, and substitutions. We determine the minimum value of C such that the attacker can coordinate its tampering action to keep the observer indefinitely confused while utilizing a finite number of attacks. Several examples are presented to demonstrate the proposed methods.


DA Transformer Distance aware Transformer. Transformer has achieved great success in the NLP field by composing various advanced models like BERT and GPT. However, Transformer and its existing variants may not be optimal in capturing token distances because the position or distance embeddings used by these methods usually cannot keep the precise information of real distances, which may not be beneficial for modeling the orders and relations of contexts. In this paper, we propose DA Transformer, which is a distance aware Transformer that can exploit the real distance. We propose to incorporate the real distances between tokens to re scale the raw self attention weights, which are computed by the relevance between attention query and key. Concretely, in different self attention heads the relative distance between each pair of tokens is weighted by different learnable parameters, which control the different preferences on long or short term information of these heads. Since the raw weighted real distances may not be optimal for adjusting self attention weights, we propose a learnable sigmoid function to map them into re scaled coefficients that have proper ranges. We first clip the raw self attention weights via the ReLU function to keep non negativity and introduce sparsity, and then multiply them with the re scaled coefficients to encode real distance information into self attention. Extensive experiments on five benchmark datasets show that DA Transformer can effectively improve the performance of many tasks and outperform the vanilla Transformer and its several variants.


On the geometric structure of currents tangent to smooth distributions. It is well known that a k dimensional smooth surface in a Euclidean space cannot be tangent to a non involutive distribution of k dimensional planes. In this paper we discuss the extension of this statement to weaker notions of surfaces, namely integral and normal currents. We find out that integral currents behave to this regard exactly as smooth surfaces, while the behaviour of normal currents is rather multifaceted. This issue is strictly related to a geometric property of the boundary of currents, which is also discussed in details.


Guidelines for Experimental Algorithmics in Network Analysis. The field of network science is a highly interdisciplinary area for the empirical analysis of network data, it draws algorithmic methodologies from several research fields. Hence, research procedures and descriptions of the technical results often differ, sometimes widely. In this paper we focus on methodologies for the experimental part of algorithm engineering for network analysis an important ingredient for a research area with empirical focus. More precisely, we unify and adapt existing recommendations from different fields and propose universal guidelines including statistical analyses for the systematic evaluation of network analysis algorithms. This way, the behavior of newly proposed algorithms can be properly assessed and comparisons to existing solutions become meaningful. Moreover, as the main technical contribution, we provide SimexPal, a highly automated tool to perform and analyze experiments following our guidelines. To illustrate the merits of SimexPal and our guidelines, we apply them in a case study we design, perform, visualize and evaluate experiments of a recent algorithm for approximating betweenness centrality, an important problem in network analysis. In summary, both our guidelines and SimexPal shall modernize and complement previous efforts in experimental algorithmics they are not only useful for network analysis, but also in related contexts.


Albert, an intermediate smart contract language for the Tezos blockchain. Tezos is a smart contract blockchain. Tezos smart contracts are written in a low level stack based language called Michelson. In this article we present Albert, an intermediate language for Tezos smart contracts which abstracts Michelson stacks as linearly typed records. We also describe its compiler to Michelson, written in Coq, that targets Mi Cho Coq, a formal specification of Michelson implemented in Coq.


Estimates for commutators of fractional differential operators via harmonic extension. This master thesis is based on the paper "Sharp commutator estimates via harmonic extensions" by Lenzmann and Schikorra, in which they proposed a method to prove estimates for commutators involving Riesz transforms, fractional Laplacians and Riesz potentials, see arXiv 1609.08547. These proofs only involve harmonic extensions to the upper half space and integration by parts next to some elementary transfromations, since the deeper theory is concentrated in a variety of trace characterization results which can be used as a blackbox. In the first half of this thesis, after collecting some elementary results for the s harmonic extension by Caffarelli and Silvestre, we use this method to prove a variety of commutator estimates, closely following Lenzmann and Schikorra except for shortening some proofs. In the second half, we prove generalized versions of the blackbox estimates listed by Lenzmann and Schikorra and discuss the different building blocks which make up these blackbox estimates, including Triebel Lizorkin and Besov Lipschitz space characterizations as well as square function estimates.


Reflexive sheaves, Hermitian Yang Mills connections, and tangent cones. In this paper we give a complete algebro geometric characterization of analytic tangent cones of admissible Hermitian Yang Mills connections over any reflexive sheaves.


The spherical image of singular varieties of bounded mean curvature. In this paper we deal with singular varieties of bounded mean curvature in the viscosity sense. They contain all varifolds of bounded generalized mean curvature. In the first part we investigate the second order properties of these varieties, obtaining results that are new also in the varifold's setting. In particular we prove that the generalized normal bundle of these varieties satisfies a natural Lusin (N) condition, which allows to extend the classical Coarea formula for the Gauss map of smooth varieties, and to introduce for all integral varifolds of bounded mean curvature a natural definition of second fundamental form, whose trace equals the generalized varifold mean curvature. In the second part, we use this machinery to extend a sharp geometric inequality of Almgren to all compact varieties of bounded mean curvature in the viscosity sense and we characterize the equality case. As a consequence we formulate sufficient conditions to conclude that the area blow up set is empty for sequences of varifolds whose first variation is controlled.


Optimal Learning Dynamics of Multi Agents in Restless Multiarmed Bandit Game. Social learning is learning through the observation of or interaction with other individuals it is critical in the understanding of the collective behaviors of humans in social physics. We study the learning process of agents in a restless multiarmed bandit (rMAB). The binary payoff of each arm changes randomly and agents maximize their payoffs by exploiting an arm with payoff 1, searching the arm at random (individual learning), or copying an arm exploited by other agents (social learning). The system has Pareto and Nash equilibria in the mixed strategy space of social and individual learning. We study several models in which agents maximize their expected payoffs in the strategy space, and demonstrate analytically and numerically that the system converges to the equilibria. We also conducted an experiment and investigated whether human participants adopt the optimal strategy. In this experiment, three participants play the game. If the reward of each group is proportional to the sum of the payoffs, the median of the social learning rate almost coincides with that of the Pareto equilibrium.


An Automorphic Classification of Real Cubic Curves. The action of ring automorphisms of the polynomial ring in two variables over the real numbers on real plane curves is considered. The orbits containing degree three polynomials are computed, with one representative per orbit being selected.


A Syllable Structured, Contextually Based Conditionally Generation of Chinese Lyrics. This paper presents a novel, syllable structured Chinese lyrics generation model given a piece of original melody. Most previously reported lyrics generation models fail to include the relationship between lyrics and melody. In this work, we propose to interpret lyrics melody alignments as syllable structural information and use a multi channel sequence to sequence model with considering both phrasal structures and semantics. Two different RNN encoders are applied, one of which is for encoding syllable structures while the other for semantic encoding with contextual sentences or input keywords. Moreover, a large Chinese lyrics corpus for model training is leveraged. With automatic and human evaluations, results demonstrate the effectiveness of our proposed lyrics generation model. To the best of our knowledge, there is few previous reports on lyrics generation considering both music and linguistic perspectives.


Efficient Evaluation of Exact Exchange for Periodic Systems via Concentric Atomic Density Fitting. The evaluation of exact (Hartree Fock, HF) exchange operator is a crucial ingredient for the accurate description of electronic structure in periodic systems through ab initio and hybrid density functional approaches. An efficient formulation of periodic HF exchange in LCAO representation presented here is based on the concentric atomic density fitting (CADF) approximation, a domain free local density fitting approach in which the product of two atomic orbitals (AOs) is approximated using a linear combination of fitting basis functions centered at the same nuclei as the AOs in that product. Significant reduction in the computational cost of exact exchange is demonstrated relative to the conventional approach due to avoiding the need to evaluate four center two electron integrals, with sub millihartree atom errors in absolute Hartree Fock energies and good cancellation of fitting errors in relative energies. Novel aspects of the evaluation of the Coulomb contribution to the Fock operator, such as the use of real two center multipole expansions and spheropole compensated unit cell densities are also described.


X ToM Explaining with Theory of Mind for Gaining Justified Human Trust. We present a new explainable AI (XAI) framework aimed at increasing justified human trust and reliance in the AI machine through explanations. We pose explanation as an iterative communication process, i.e. dialog, between the machine and human user. More concretely, the machine generates sequence of explanations in a dialog which takes into account three important aspects at each dialog turn (a) human's intention (or curiosity) (b) human's understanding of the machine and (c) machine's understanding of the human user. To do this, we use Theory of Mind (ToM) which helps us in explicitly modeling human's intention, machine's mind as inferred by the human as well as human's mind as inferred by the machine. In other words, these explicit mental representations in ToM are incorporated to learn an optimal explanation policy that takes into account human's perception and beliefs. Furthermore, we also show that ToM facilitates in quantitatively measuring justified human trust in the machine by comparing all the three mental representations. We applied our framework to three visual recognition tasks, namely, image classification, action recognition, and human body pose estimation. We argue that our ToM based explanations are practical and more natural for both expert and non expert users to understand the internal workings of complex machine learning models. To the best of our knowledge, this is the first work to derive explanations using ToM. Extensive human study experiments verify our hypotheses, showing that the proposed explanations significantly outperform the state of the art XAI methods in terms of all the standard quantitative and qualitative XAI evaluation metrics including human trust, reliance, and explanation satisfaction.


Attention Beam An Image Captioning Approach. The aim of image captioning is to generate textual description of a given image. Though seemingly an easy task for humans, it is challenging for machines as it requires the ability to comprehend the image (computer vision) and consequently generate a human like description for the image (natural language understanding). In recent times, encoder decoder based architectures have achieved state of the art results for image captioning. Here, we present a heuristic of beam search on top of the encoder decoder based architecture that gives better quality captions on three benchmark datasets Flickr8k, Flickr30k and MS COCO.


When Virtual Therapy and Art Meet A Case Study of Creative Drawing Game in Virtual Environments. There have been a resurge lately on virtual therapy and other virtual and tele medicine services due to the new normal of practicing 'shelter at home'. In this paper, we propose a creative drawing game for virtual therapy and investigate user's comfort and movement freedom in a pilot study. In a mixed design study, healthy participants (N 16, 8 females) completed one of the easy or hard trajectories of the virtual therapy game in standing and seated arrangements using a virtual reality headset. The results from participants' movement accuracy, task completion time, and usability questionnaires indicate that participants had significant performance differences on two levels of the game based on its difficulty (between subjects factor), but no difference in seated and standing configurations (within subjects factor). Also, the hard mode was more favorable among participants. This work offers implications on virtual reality and 3D interactive systems, with specific contributions to virtual therapy, and serious games for healthcare applications.


Generator Versus Segmentor Pseudo healthy Synthesis. This paper investigates the problem of pseudo healthy synthesis that is defined as synthesizing a subject specific pathology free image from a pathological one. Recent approaches based on Generative Adversarial Network (GAN) have been developed for this task. However, these methods will inevitably fall into the trade off between preserving the subject specific identity and generating healthy like appearances. To overcome this challenge, we propose a novel adversarial training regime, Generator versus Segmentor (GVS), to alleviate this trade off by a divide and conquer strategy. We further consider the deteriorating generalization performance of the segmentor throughout the training and develop a pixel wise weighted loss by muting the well transformed pixels to promote it. Moreover, we propose a new metric to measure how healthy the synthetic images look. The qualitative and quantitative experiments on the public dataset BraTS demonstrate that the proposed method outperforms the existing methods. Besides, we also certify the effectiveness of our method on datasets LiTS. Our implementation and pre trained networks are publicly available at https github.com Au3C2 Generator Versus Segmentor.


Knowledge Integration Networks for Action Recognition. In this work, we propose Knowledge Integration Networks (referred as KINet) for video action recognition. KINet is capable of aggregating meaningful context features which are of great importance to identifying an action, such as human information and scene context. We design a three branch architecture consisting of a main branch for action recognition, and two auxiliary branches for human parsing and scene recognition which allow the model to encode the knowledge of human and scene for action recognition. We explore two pre trained models as teacher networks to distill the knowledge of human and scene for training the auxiliary tasks of KINet. Furthermore, we propose a two level knowledge encoding mechanism which contains a Cross Branch Integration (CBI) module for encoding the auxiliary knowledge into medium level convolutional features, and an Action Knowledge Graph (AKG) for effectively fusing high level context information. This results in an end to end trainable framework where the three tasks can be trained collaboratively, allowing the model to compute strong context knowledge efficiently. The proposed KINet achieves the state of the art performance on a large scale action recognition benchmark Kinetics 400, with a top 1 accuracy of 77.8 . We further demonstrate that our KINet has strong capability by transferring the Kinetics trained model to UCF 101, where it obtains 97.8 top 1 accuracy.


Federated Learning for Edge Networks Resource Optimization and Incentive Mechanism. Recent years have witnessed a rapid proliferation of smart Internet of Things (IoT) devices. IoT devices with intelligence require the use of effective machine learning paradigms. Federated learning can be a promising solution for enabling IoT based smart applications. In this paper, we present the primary design aspects for enabling federated learning at network edge. We model the incentive based interaction between a global server and participating devices for federated learning via a Stackelberg game to motivate the participation of the devices in the federated learning process. We present several open research challenges with their possible solutions. Finally, we provide an outlook on future research.


A classification of isogeny torsion graphs of mathbb Q isogeny classes of elliptic curves. Let mathcal E be a mathbb Q isogeny class of elliptic curves defined over mathbb Q . The isogeny graph associated to mathcal E is a graph which has a vertex for each elliptic curve in the mathbb Q isogeny class mathcal E , and an edge for each cyclic mathbb Q isogeny of prime degree between elliptic curves in the isogeny class, with the degree recorded as a label of the edge. In this paper, we define an isogeny torsion graph to be an isogeny graph where, in addition, we label each vertex with the abstract group structure of the torsion subgroup over mathbb Q of the corresponding elliptic curve. Then, the main result of the article is a classification of all the possible isogeny torsion graphs that occur for mathbb Q isogeny classes of elliptic curves defined over the rationals.


Behavior Sequence Transformer for E commerce Recommendation in Alibaba. Deep learning based methods have been widely used in industrial recommendation systems (RSs). Previous works adopt an Embedding MLP paradigm raw features are embedded into low dimensional vectors, which are then fed on to MLP for final recommendations. However, most of these works just concatenate different features, ignoring the sequential nature of users' behaviors. In this paper, we propose to use the powerful Transformer model to capture the sequential signals underlying users' behavior sequences for recommendation in Alibaba. Experimental results demonstrate the superiority of the proposed model, which is then deployed online at Taobao and obtain significant improvements in online Click Through Rate (CTR) comparing to two baselines.


Obstructions to representations up to homotopy and ideals. This paper considers the Pontryagin characters of graded vector bundles of finite rank, in the cohomology vector spaces of a Lie algebroid over the same base. These Pontryagin characters vanish if the graded vector bundle carries a representation up to homotopy of the Lie algebroid. As a consequence, this gives a strong obstruction to the existence of a representation up to homotopy on a graded vector bundle of finite rank. In particular, if a graded vector bundle E 0 oplus F 1 to M carries a 2 term representation up to homotopy of a Lie algebroid A to M , then all the (classical) A Pontryagin classes of E and F must coincide. This paper generalises as well Bott's vanishing theorem to the setting of Lie algebroid representations (up to homotopy) on arbitrary vector bundles. As an application, the main theorems induce new obstructions to the existence of infinitesimal ideal systems in a given Lie algebroid.


Bifurcation analysis of a free boundary model of plaque formation associated with the cholesterol ratio. The low density lipoprotein (LDL) high density lipoprotein (HDL) cholesterol ratio has been shown a high correlation with the cardiovascular risk assessment. Is it possible to quantify the correlation mathematically? In this paper, we develop a bifurcation analysis for a mathematical model of the plaque formation with a free boundary in the early stage of atherosclerosis. This bifurcation analysis, to the ratio of LDL HDL, is based on the explicit formulation of radially symmetric steady state solutions. By performing the perturbation analysis to these solutions, we establish the existence of bifurcation branches and derive a theoretical condition that a bifurcation occurs for different modes. Then we also analyze the stability of radially symmetric steady state solutions and conduct numerical simulations to verify all the theoretical results.


Hyperfine structure in the H 2 and HD molecular ions at m alpha 6 order. A complete effective Hamiltonian for relativistic corrections at orders m alpha 6 and m alpha 6(m M) in a one electron molecular system is derived from the NRQED Lagrangian. It includes spin independent corrections to the energy levels and spin spin scalar interactions contributing to the hyperfine splitting, both of which had been studied previously. In addition, corrections to electron spin orbit and spin spin tensor interactions are newly obtained. This allows improving the hyperfine structure theory in the hydrogen molecular ions. Improved values of the spin orbit hyperfine coefficient are calculated for a few transitions of current experimental interest.


Controlled edge dependent stacking of WS2 WS2 Homo and WS2 WSe2 Hetero structures A Computational Study. Transition Metal Dichalcogenides (TMDs) are one of the most studied two dimensional materials in the last 5 10 years due to their extremely interesting layer dependent properties. Despite the presence of vast research work on TMDs, the complex relationship between the electrochemical and physical properties make them the subject of further research. Our main objective is to provide a better insight into the electronic structure of TMDs. This will help us better understand the stability of the bilayer post growth homo hetero products based on the various edge termination, and different stacking of the two layers. In this regard, two Tungsten (W) based non periodic chalcogenide flakes (sulfides and selenides) were considered. An in depth analysis of their different edge termination and stacking arrangement was performed via Density Functional Theory method using VASP software. Our finding indicates the preference of chalcogenide (c ) terminated structures over the metal (m ) terminated structures for both homo and hetero layers, and thus strongly suggests the nonexistence of the m terminated TMDs bilayer products.


CACHE SNIPER Accurate timing control of cache evictions. Microarchitectural side channel attacks have been very prominent in security research over the last few years. Caches have been an outstanding covert channel, as they provide high resolution and generic cross core leakage even with simple user mode code execution privileges. To prevent these generic cross core attacks, all major cryptographic libraries now provide countermeasures to hinder key extraction via cross core cache attacks, for instance avoiding secret dependent access patterns and prefetching data. In this paper, we show that implementations protected by 'good enough' countermeasures aimed at preventing simple cache attacks are still vulnerable. We present a novel attack that uses a special timing technique to determine when an encryption has started and then evict the data precisely at the desired instant. This new attack does not require special privileges nor explicit synchronization between the attacker and the victim. One key improvement of our attack is a method to evict data from the cache with a single memory access and in absence of shared memory by leveraging the transient capabilities of TSX and relying on the recently reverse engineered L3 replacement policy. We demonstrate the efficiency by performing an asynchronous last level cache attack to extract an RSA key from the latest wolfSSL library, which has been especially adapted to avoid leaky access patterns, and by extracting an AES key from the S Box implementation included in OpenSSL bypassing the per round prefetch intended as a protection against cache attacks.


LSTM vs. GRU vs. Bidirectional RNN for script generation. Scripts are an important part of any TV series. They narrate movements, actions and expressions of characters. In this paper, a case study is presented on how different sequence to sequence deep learning models perform in the task of generating new conversations between characters as well as new scenarios on the basis of a script (previous conversations). A comprehensive comparison between these models, namely, LSTM, GRU and Bidirectional RNN is presented. All the models are designed to learn the sequence of recurring characters from the input sequence. Each input sequence will contain, say "n" characters, and the corresponding targets will contain the same number of characters, except, they will be shifted one character to the right. In this manner, input and output sequences are generated and used to train the models. A closer analysis of explored models performance and efficiency is delineated with the help of graph plots and generated texts by taking some input string. These graphs describe both, intraneural performance and interneural model performance for each model.


Stable functors and cohomology theory in Grothendieck categories. In this paper, we introduce and study a relative complete cohomology theory in Grothendieck categories. Some properties of this cohomology including vanishing and balancedness are given. As an application of our approach, we prove a new balanced result for Tate cohomology in Grothendieck categories.


Learning Bayesian networks from demographic and health survey data. Child mortality from preventable diseases such as pneumonia and diarrhoea in low and middle income countries remains a serious global challenge. We combine knowledge with available Demographic and Health Survey (DHS) data from India, to construct Causal Bayesian Networks (CBNs) and investigate the factors associated with childhood diarrhoea. We make use of freeware tools to learn the graphical structure of the DHS data with score based, constraint based, and hybrid structure learning algorithms. We investigate the effect of missing values, sample size, and knowledge based constraints on each of the structure learning algorithms and assess their accuracy with multiple scoring functions. Weaknesses in the survey methodology and data available, as well as the variability in the CBNs generated by the different algorithms, mean that it is not possible to learn a definitive CBN from data. However, knowledge based constraints are found to be useful in reducing the variation in the graphs produced by the different algorithms, and produce graphs which are more reflective of the likely influential relationships in the data. Furthermore, valuable insights are gained into the performance and characteristics of the structure learning algorithms. Two score based algorithms in particular, TABU and FGES, demonstrate many desirable qualities a) with sufficient data, they produce a graph which is similar to the reference graph, b) they are relatively insensitive to missing values, and c) behave well with knowledge based constraints. The results provide a basis for further investigation of the DHS data and for a deeper understanding of the behaviour of the structure learning algorithms when applied to real world settings.


Non equilibrium statistical physics, transitory epigenetic landscapes, and cell fate decision dynamics. Statistical physics provides a useful perspective for the analysis of many complex systems it allows us to relate microscopic fluctuations to macroscopic observations. Developmental biology, but also cell biology more generally, are examples where apparently robust behaviour emerges from highly complex and stochastic sub cellular processes. Here we attempt to make connections between different theoretical perspectives to gain qualitative insights into the types of cell fate decision making processes that are at the heart of stem cell and developmental biology. We discuss both dynamical systems as well as statistical mechanics perspectives on the classical Waddington or epigenetic landscape. We find that non equilibrium approaches are required to overcome some of the shortcomings of classical equilibrium statistical thermodynamics or statistical mechanics in order to shed light on biological processes, which, almost by definition, are typically far from equilibrium.


A Scalable Molecular Force Field Parameterization Method Based on Density Functional Theory and Quantum Level Machine Learning. Fast and accurate molecular force field (FF) parameterization is still an unsolved problem. Accurate FFs are not generally available for all molecules, like novel druglike molecules. While methods based on quantum mechanics (QM) exist to parameterize them with better accuracy, they are computationally expensive and slow, which limits applicability to a small number of molecules. Here, we present an automated FF parameterization method which can utilize either DFT calculations or approximate QM energies produced by different neural network potentials (NNPs), to obtain improved parameters for molecules. We demonstrate that for the case of torchani ANI 1x NNP, we can parameterize small molecules in a fraction of time compared with an equivalent parameterization using DFT QM calculations while producing more accurate parameters than FF (GAFF2). We expect our method to be of critical importance in computational structure based drug discovery. The current version is available at PlayMolecule (www.playmolecule.org) and implemented in HTMD, allowing to parameterize molecules with different QM and NNP options.


Forensic Scanner Identification Using Machine Learning. Due to the increasing availability and functionality of image editing tools, many forensic techniques such as digital image authentication, source identification and tamper detection are important for forensic image analysis. In this paper, we describe a machine learning based system to address the forensic analysis of scanner devices. The proposed system uses deep learning to automatically learn the intrinsic features from various scanned images. Our experimental results show that high accuracy can be achieved for source scanner identification. The proposed system can also generate a reliability map that indicates the manipulated regions in an scanned image.


Long wavelength infrared photovoltaic heterodyne receivers using patch antenna quantum cascade detectors. Quantum cascade detectors (QCD) are unipolar infrared devices where the transport of the photo excited carriers takes place through confined electronic states, without an applied bias. In this photovoltaic mode, the detector's noise is not dominated by a dark shot noise process, therefore, performances are less degraded at high temperature with respect to photoconductive detectors. This work describes a 9 um QCD embedded into a patch antenna metamaterial which operates with state of the art performances. The metamaterial gathers photons on a collection area, Acoll, much bigger than the geometrical area of the detector, improving the signal to noise ratio up to room temperature. The background limited detectivity at 83 K is 5.5 x 10 10 cm Hz 1 2 W 1, while at room temperature, the responsivity is 50 mA W at 0 V bias. Patch antenna QCD is an ideal receiver for a heterodyne detection set up, where a signal at a frequency 1.4 GHz and T 295 K is reported as first demonstration of uncooled 9um photovoltaic receivers with GHz electrical bandwidth. These findings guide the research towards uncooled IR quantum limited detection.


Twenty Years of Network Science A Bibliographic and Co Authorship Network Analysis. Two decades ago three pioneering papers turned the attention to complex networks and initiated a new era of research, establishing an interdisciplinary field called network science. Namely, these highly cited seminal papers were written by Watts Strogatz, Barab 'asi Albert, and Girvan Newman on small world networks, on scale free networks and on the community structure of complex networks, respectively. In the past 20 years due to the multidisciplinary nature of the field a diverse but not divided network science community has emerged. In this paper, we investigate how this community has evolved over time with respect to speed, diversity and interdisciplinary nature as seen through the growing co authorship network of network scientists (here the notion refers to a scholar with at least one paper citing at least one of the three aforementioned milestone papers). After providing a bibliographic analysis of 31,763 network science papers, we construct the co authorship network of 56,646 network scientists and we analyze its topology and dynamics. We shed light on the collaboration patterns of the last 20 years of network science by investigating numerous structural properties of the co authorship network and by using enhanced data visualization techniques. We also identify the most central authors, the largest communities, investigate the spatiotemporal changes, and compare the properties of the network to scientometric indicators.


KNN DBSCAN a DBSCAN in high dimensions. Clustering is a fundamental task in machine learning. One of the most successful and broadly used algorithms is DBSCAN, a density based clustering algorithm. DBSCAN requires epsilon nearest neighbor graphs of the input dataset, which are computed with range search algorithms and spatial data structures like KD trees. Despite many efforts to design scalable implementations for DBSCAN, existing work is limited to low dimensional datasets, as constructing epsilon nearest neighbor graphs is expensive in high dimensions. In this paper, we modify DBSCAN to enable use of kappa nearest neighbor graphs of the input dataset. The kappa nearest neighbor graphs are constructed using approximate algorithms based on randomized projections. Although these algorithms can become inaccurate or expensive in high dimensions, they possess a much lower memory overhead than constructing epsilon nearest neighbor graphs ( mathcal O (nk) vs. mathcal O (n 2) ). We delineate the conditions under which k NN DBSCAN produces the same clustering as DBSCAN. We also present an efficient parallel implementation of the overall algorithm using OpenMP for shared memory and MPI for distributed memory parallelism. We present results on up to 16 billion points in 20 dimensions, and perform weak and strong scaling studies using synthetic data. Our code is efficient in both low and high dimensions. We can cluster one billion points in 3D in less than one second on 28K cores on the Frontera system at the Texas Advanced Computing Center (TACC). In our largest run, we cluster 65 billion points in 20 dimensions in less than 40 seconds using 114,688 x86 cores on TACC's Frontera system. Also, we compare with a state of the art parallel DBSCAN code on 20d 4M point dataset, our code is up to 37 times faster.


CAFE Coarse to Fine Neural Symbolic Reasoning for Explainable Recommendation. Recent research explores incorporating knowledge graphs (KG) into e commerce recommender systems, not only to achieve better recommendation performance, but more importantly to generate explanations of why particular decisions are made. This can be achieved by explicit KG reasoning, where a model starts from a user node, sequentially determines the next step, and walks towards an item node of potential interest to the user. However, this is challenging due to the huge search space, unknown destination, and sparse signals over the KG, so informative and effective guidance is needed to achieve a satisfactory recommendation quality. To this end, we propose a CoArse to FinE neural symbolic reasoning approach (CAFE). It first generates user profiles as coarse sketches of user behaviors, which subsequently guide a path finding process to derive reasoning paths for recommendations as fine grained predictions. User profiles can capture prominent user behaviors from the history, and provide valuable signals about which kinds of path patterns are more likely to lead to potential items of interest for the user. To better exploit the user profiles, an improved path finding algorithm called Profile guided Path Reasoning (PPR) is also developed, which leverages an inventory of neural symbolic reasoning modules to effectively and efficiently find a batch of paths over a large scale KG. We extensively experiment on four real world benchmarks and observe substantial gains in the recommendation performance compared with state of the art methods.


Direct Vlasov solvers. In these proceedings we will describe the theory and practical steps required to build Vlasov solvers such as those commonly used to compute coherent instabilities in synchrotrons. Thanks to a Hamiltonian formalism, we will derive a compact and general form of the linearized Vlasov equation, written using Poisson brackets. This in turn will be the basis of a procedure to build Vlasov solvers, applied to the specific example of transverse instabilities arising from beam coupling impedance.


On the correctness of Egalitarian Paxos. This paper identifies a problem in both the TLA specification and the implementation of the Egalitarian Paxos protocol. It is related to how replicas switch from one ballot to another when computing the dependencies of a command. The problem may lead replicas to diverge and break the linearizability of the replicated service.


A Data Driven and Integrated Evaluation of Area wide Impacts of Double Parking Using Macroscopic and Microscopic Models. Double parking that often negatively affects traffic operations and safety is not a new phenomenon on urban streets. This study proposes a novel data driven integrated framework for estimating the actual frequency of double parking so that both microscopic and macroscopic models can be utilized to quantify area wide impacts in the presence of double parking. The findings of the study can provide transportation agencies with useful insights on identifying locations that will experience the greatest benefits by removing problematic double parking. As a result, various parking enforcement and management strategies can be planned more effectively.


Putting Together the Pieces A Concept for Holistic Industrial Intrusion Detection. Besides the advantages derived from the ever present communication properties, it increases the attack surface of a network as well. As industrial protocols and systems were not designed with security in mind, spectacular attacks on industrial systems occurred over the last years. Most industrial communication protocols do not provide means to ensure authentication or encryption. This means attackers with access to a network can read and write information. Originally not meant to be connected to public networks, the use cases of Industry 4.0 require interconnectivity, often through insecure public networks. This lead to an increasing interest in information security products for industrial applications. In this work, the concept for holistic intrusion detection methods in an industrial context is presented. It is based on different works considering several aspects of industrial environments and their capabilities to identify intrusions as an anomaly in network or process data. These capabilities are based on preceding experiments on real and synthetic data. In order to justify the concept, an overview of potential and actual attack vectors and attacks on industrial systems is provided. It is shown that different aspects of industrial facilities, e.g. office IT, shop floor OT, firewalled connections to customers and partners are analysed as well as the different layers of the automation pyramid require different methods to detect attacks. Additionally, the singular steps of an attack on industrial applications are characterised. Finally, a resulting concept for integration of these methods is proposed, providing the means to detect the different stages of an attack by different means.


Relevant Region Exploration On General Cost maps For Sampling Based Motion Planning. Asymptotically optimal sampling based planners require an intelligent exploration strategy to accelerate convergence. After an initial solution is found, a necessary condition for improvement is to generate new samples in the so called "Informed Set". However, Informed Sampling can be ineffective in focusing search if the chosen heuristic fails to provide a good estimate of the solution cost. This work proposes an algorithm to sample the "Relevant Region" instead, which is a subset of the Informed Set. The Relevant Region utilizes cost to come information from the planner's tree structure, reduces dependence on the heuristic, and further focuses the search. Benchmarking tests in uniform and general cost space settings demonstrate the efficacy of Relevant Region sampling.


Two Layer Neural Networks for Partial Differential Equations Optimization and Generalization Theory. The problem of solving partial differential equations (PDEs) can be formulated into a least squares minimization problem, where neural networks are used to parametrize PDE solutions. A global minimizer corresponds to a neural network that solves the given PDE. In this paper, we show that the gradient descent method can identify a global minimizer of the least squares optimization for solving second order linear PDEs with two layer neural networks under the assumption of over parametrization. We also analyze the generalization error of the least squares optimization for second order linear PDEs and two layer neural networks, when the right hand side function of the PDE is in a Barron type space and the least squares optimization is regularized with a Barron type norm, without the over parametrization assumption.


Stated skein modules of marked 3 manifolds surfaces, a survey. We give a survey of some old and new results about the stated skein modules algebras of 3 manifolds surfaces. For generic quantum parameter, we discuss the splitting homomorphism for the 3 manifold case, general structures of the stated skein algebras of marked surfaces (or bordered punctured surfaces) and their embeddings into quantum tori. For roots of 1 quantum parameter, we discuss the Frobenius homomorphism (for both marked 3 manifolds and marked surfaces), describe the center of the skein algebra of marked surfaces, the dimension of the skein algebra over the center, and the representation theory of the skein algebra. In particular, we show that the skein algebra of non closed marked surface at any root of 1 is a maximal order. We give a full description of the Azumaya locus of the skein algebra of the puncture torus and give partial results for closed surfaces.


On certain self orthogonal AG codes with applications to Quantum error correcting codes. In this paper a construction of quantum codes from self orthogonal algebraic geometry codes is provided. Our method is based on the CSS construction as well as on some peculiar properties of the underlying algebraic curves, named Swiss curves. Several classes of well known algebraic curves with many rational points turn out to be Swiss curves. Examples are given by Castle curves, GK curves, generalized GK curves and the Abdon Bezerra Quoos maximal curves. Applications of our method to these curves are provided. Our construction extends a previous one due to Hernando, McGuire, Monserrat, and Moyano Fernandez.


AoI Delay Tradeoff in Mobile Edge Caching with Freshness Aware Content Refreshing. Mobile edge caching can effectively reduce service delay but may introduce information staleness, calling for timely content refreshing. However, content refreshing consumes additional transmission resources and may degrade the delay performance of mobile systems. In this work, we propose a freshness aware refreshing scheme to balance the service delay and content freshness measured by Age of Information (AoI). Specifically, the cached content items will be refreshed to the up to date version upon user requests if the AoI exceeds a certain threshold (named as refreshing window). The average AoI and service delay are derived in closed forms approximately, which reveals an AoI delay tradeoff relationship with respect to the refreshing window. In addition, the refreshing window is optimized to minimize the average delay while meeting the AoI requirements, and the results indicate to set a smaller refreshing window for the popular content items. Extensive simulations are conducted on the OMNeT platform to validate the analytical results. The results indicate that the proposed scheme can restrain frequent refreshing as the request arrival rate increases, whereby the average delay can be reduced by around 80 while maintaining the AoI below one second in heavily loaded scenarios.


Lightlike manifolds and Cartan geometries. Lightlike Cartan geometries are introduced as Cartan geometries modelled on the future lightlike cone in Lorentz Minkowski spacetime. Then, we provide an approach to the study of lightlike manifolds from this point of view. It is stated that every lightlike Cartan geometry on a manifold N provides a lightlike metric h with radical distribution globally spanned by a vector field Z . For lightlike hypersurfaces of a Lorentz manifold, we give the condition that characterizes when the pull back of the Levi Civita connection form of the ambient manifold is a lightlike Cartan connection on such hypersurface. In the particular case that a lightlike hypersurface is properly totally umbilical, this construction essentially returns the original lightlike metric. From the intrinsic point of view, starting from a given lightlike manifold (N,h) , we show a method to construct a family of ambient Lorentzian manifolds that realize (N,h) as a hypersurface. This method is inspired on the Feffermann Graham ambient metric construction in conformal geometry and provides a lightlike Cartan geometry on the original manifold when (N,h) is generic.


An Expert System Approach for determine the stage of UiTM Perlis Palapes Cadet Performance and Ranking Selection. The palapes cadets are one of the uniform organizations in UiTM Perlis for extra curricular activities. The palapes cadets arrange their organization in a hierarchy according to grade. Senior uniform officer (SUO) is the highest rank, followed by a junior uniform officer (JUO), sergeant, corporal, lance corporal, and lastly, cadet officer, which is the lowest rank. The palapes organization has several methods to measure performance toward promotion to a higher rank, whether individual performance or in a group. Cadets are selected for promotion based on demonstrated leadership abilities, acquired skills, physical fitness, and comprehension of information as measured through standardized testing. However, this method is too complicated when manually assessed by a trainer or coach. Therefore, this study will propose an expert system, which is one of the artificial intelligence techniques that can recognize the readiness and progression of a palapes cadet.


Universal Algebra in UniMath. We present an ongoing effort to implement Universal Algebra in the UniMath system. Our aim is to develop a general framework for formalizing and studying Universal Algebra in a proof assistant. By constituting a formal system for isolating the invariants of the theory we are interested in that is, general algebraic structures modulo isomorphism Univalent Mathematics seems to provide a suitable environment to carry on our endeavour.


Monochromatic Equilateral Triangles in the Unit Distance Graph. Let chi Delta ( mathbb R n ) denote the minimum number of colors needed to color mathbb R n so that there will not be a monochromatic equilateral triangle with side length 1 . Using the slice rank method, we reprove a result of Frankl and Rodl, and show that chi Delta left( mathbb R n right) grows exponentially with n . This technique substantially improves upon the best known quantitative lower bounds for chi Delta left( mathbb R n right) , and we obtain chi Delta left( mathbb R n right) (1.01446 o(1)) n .


Evaluation of Surrogate Models for Multi fin Flapping Propulsion Systems. The aim of this study is to develop surrogate models for quick, accurate prediction of thrust forces generated through flapping fin propulsion for given operating conditions and fin geometries. Different network architectures and configurations are explored to model the training data separately for the lead fin and rear fin of a tandem fin setup. We progressively improve the data representation of the input parameter space for model predictions. The models are tested on three unseen fin geometries and the predictions validated with computational fluid dynamics (CFD) data. Finally, the orders of magnitude gains in computational performance of these surrogate models, compared to experimental and CFD runs, vs their tradeoff with accuracy is discussed within the context of this tandem fin configuration.


How far do Activated Random Walkers spread from a single source?. Unlike many particle systems, Activated Random Walk has nontrivial behavior even in one spatial dimension. We prove inner and outer bounds on the spread of n activated random walkers from a single source in Z. The inner bound involves a comparison with the stationary distribution of activated random walkers on a finite interval, while the outer bound involves a comparison with the stabilization of an infinite Bernoulli configuration of activated random walkers on Z.


Concentration versus oscillation effects in brittle damage. This work is concerned with an asymptotic analysis, in the sense of Gamma convergence, of a sequence of variational models of brittle damage in the context of linearized elasticity. The study is performed as the damaged zone concentrates into a set of zero volume and, at the same time and to the same order varepsilon , the stiffness of the damaged material becomes small. Three main features make the analysis highly nontrivial at varepsilon fixed, minimizing sequences of each brittle damage model oscillate and develop microstructures as varepsilon to 0 , concentration and saturation of damage are favoured and the competition of these phenomena translates into a degeneration of the growth of the elastic energy, which passes from being quadratic (at varepsilon fixed) to being of linear growth type (in the limit). Consequently, homogenization effects interact with singularity formation in a nontrivial way, which requires new methods of analysis. In particular, the interaction of homogenization with singularity formation in the framework of linearized elasticity appears to not have been considered in the literature so far. We explicitly identify the Gamma limit in two and three dimensions for isotropic Hooke tensors. The expression of the limit effective energy turns out to be of Hencky plasticity type. We further consider the regime where the divergence remains square integrable in the limit, which leads to a Tresca type model.


Asymptotic stability of controlled differential equations. Part I Young integrals. We provide a unified analytic approach to study stationary states of controlled differential equations driven by rough paths, using the framework of random dynamical systems and random attractors. Part I deals with driving paths of finite p variations with 1 leq p 2 so that the integrals are interpreted in the Young sense. Our method helps to generalize recent results cite GAKLBSch2010 , cite ducGANSch18 , cite duchongcong18 on the existence of the global pullback attractors for the generated random dynamical systems. We also prove sufficient conditions for the attractor to be a singleton, thus the pathwise convergence is in both pullback and forward senses.


An ensemble Kalman filter approach based on level set parameterization for acoustic source identification using multiple frequency information. The spatial dependent unknown acoustic source is reconstructed according noisy multiple frequency data on a remote closed surface. Assume that the unknown function is supported on a bounded domain. To determine the support, we present a statistical inversion algorithm, which combines the ensemble Kalman filter approach with level set technique. Several numerical examples show that the proposed method give good numerical reconstruction.


Recurrence relation for plethysm. We describe a recurrence formula for the plethysm h 3 h n . The proof is based on the original formula by Thrall.


How Safe are European Safe Bonds? An Analysis from the Perspective of Modern Portfolio Credit Risk Models. Several proposals for the reform of the euro area advocate the creation of a market in synthetic securities backed by portfolios of sovereign bonds. Most debated are the so called European Safe Bonds or ESBies proposed by Brunnermeier, Langfield, Pagano,Reis, Van Nieuwerburgh and Vayanos (2017). The potential benefits of ESBies and other bond backed securities hinge on the assertion that these products are really safe. In this paper we provide a comprehensive quantitative study of the risks associated with ESBies and related products, using an affine credit risk model with regime switching as vehicle for our analysis. We discuss a recent proposal of Standard and Poors for the rating of ESBies, we analyse the impact of model parameters and attachment points on the size and the volatility of the credit spread of ESBies and we consider several approaches to assess the market risk of ESBies. Moreover, we compare ESBies to synthetic securities created by pooling the senior tranche of national bonds as suggested by Leandro and Zettelmeyer(2019). The paper concludes with a brief discussion of the policy implications from our analysis.


A Common Gene Expression Signature Analysis Method for Multiple Types of Cancer. Mining gene expression profiles has proven valuable for identifying signatures serving as surrogates of cancer phenotypes. However, the similarities of such signatures across different cancer types have not been strong enough to conclude that they represent a universal biological mechanism shared among multiple cancer types. Here we describe a network based approach that explores gene to gene connections in multiple cancer datasets while maximizing the overall association of the subnetwork with clinical outcomes. With the dataset of The Cancer Genome Atlas (TCGA), we studied the characteristics of common gene expression of three types of cancers Rectum adenocarcinoma (READ), Breast invasive carcinoma (BRCA) and Colon adenocarcinoma (COAD). By analyzing several pairs of highly correlated genes after filtering and clustering work, we found that the co expressed genes across multiple types of cancers point to particular biological mechanisms related to cancer cell progression , suggesting that they represent important attributes of cancer in need of being elucidated for potential applications in diagnostic, prognostic and therapeutic products applicable to multiple cancer types.


A Morpho Syntactically Informed LSTM CRF Model for Named Entity Recognition. We propose a morphologically informed model for named entity recognition, which is based on LSTM CRF architecture and combines word embeddings, Bi LSTM character embeddings, part of speech (POS) tags, and morphological information. While previous work has focused on learning from raw word input, using word and character embeddings only, we show that for morphologically rich languages, such as Bulgarian, access to POS information contributes more to the performance gains than the detailed morphological information. Thus, we show that named entity recognition needs only coarse grained POS tags, but at the same time it can benefit from simultaneously using some POS information of different granularity. Our evaluation results over a standard dataset show sizable improvements over the state of the art for Bulgarian NER.


Hunt's Formula for SU q(N) and U q(N) . We provide a Hunt type formula for the infinitesimal generators of L 'evy process on the quantum groups SU q(N) and U q(N) . In particular, we obtain a decomposition of such generators into a gaussian part and a jump type' part determined by a linear functional that resembles the functional induced by the L 'evy measure. The jump part on SU q(N) decomposes further into parts that live on the quantum subgroups SU q(n) , n le N . Like in the classical Hunt formula for locally compact Lie groups, the ingredients become unique once a certain projection is chosen. There are analogous result for U q(N) .


Transformer Encoder Detector Module Using Context to Improve Robustness to Adversarial Attacks on Object Detection. Deep neural network approaches have demonstrated high performance in object recognition (CNN) and detection (Faster RCNN) tasks, but experiments have shown that such architectures are vulnerable to adversarial attacks (FFF, UAP) low amplitude perturbations, barely perceptible by the human eye, can lead to a drastic reduction in labeling performance. This article proposes a new context module, called textit Transformer Encoder Detector Module , that can be applied to an object detector to (i) improve the labeling of object instances and (ii) improve the detector's robustness to adversarial attacks. The proposed model achieves higher mAP, F1 scores and AUC average score of up to 13 compared to the baseline Faster RCNN detector, and an mAP score 8 points higher on images subjected to FFF or UAP attacks due to the inclusion of both contextual and visual features extracted from scene and encoded into the model. The result demonstrates that a simple ad hoc context module can improve the reliability of object detectors significantly.


A hierarchical approach to the a posteriori error estimation of isogeometric Kirchhoff plates and Kirchhoff Love shells. This work focuses on the development of a posteriori error estimates for fourth order, elliptic, partial differential equations. In particular, we propose a novel algorithm to steer an adaptive simulation in the context of Kirchhoff plates and Kirchhoff Love shells by exploiting the local refinement capabilities of hierarchical B splines. The method is based on the solution of an auxiliary residual like variational problem, formulated by means of a space of localized spline functions. This space is characterized by C 1 continuous B splines with compact support on each active element of the hierarchical mesh. We demonstrate the applicability of the proposed estimator to Kirchhoff plates and Kirchhoff Love shells by studying several benchmark problems which exhibit both smooth and singular solutions. In all cases, we obtain optimal asymptotic rates of convergence for the error measured in the energy norm and an excellent approximation of the true error.


Asymptotic expansions of solutions of the Yamabe equation and the sigma k Yamabe equation near isolated singular points. We study asymptotic behaviors of positive solutions to the Yamabe equation and the sigma k Yamabe equation near isolated singular points and establish expansions up to arbitrary orders. Such results generalize an earlier pioneering work by Caffarelli, Gidas, and Spruck, and a work by Korevaar, Mazzeo, Pacard, and Schoen, on the Yamabe equation, and a work by Han, Li, and Teixeira on the sigma k Yamabe equation. The study is based on a combination of classification of global singular solutions and an analysis of linearized operators at these global singular solutions. Such linearized equations are uniformly elliptic near singular points for 1 leq k leq n 2 and become degenerate for n 2 k leq n . In a significant portion of the paper, we establish a degree 1 expansion for the sigma k Yamabe equation for n 2 k n , generalizing a similar result for k 1 by Korevaar, Mazzeo, Pacard, and Schoen and for 2 leq k leq n 2 by Han, Li, and Teixeira.


ConfigTron Tackling network diversity with heterogeneous configurations. The web serving protocol stack is constantly changing and evolving to tackle technological shifts in networking infrastructure and website complexity. As a result of this evolution, the web serving stack includes a plethora of protocols and configuration parameters that enable the web serving stack to address a variety of realistic network conditions. Yet, today, most content providers have adopted a "one size fits all" approach to configuring the networking stack of their user facing web servers (or at best employ moderate tuning), despite the significant diversity in end user networks and devices. In this paper, we revisit this problem and ask a more fundamental question Are there benefits to tuning the network stack? If so, what system design choices and algorithmic ensembles are required to enable modern content provider to dynamically and flexibly tune their protocol stacks. We demonstrate through substantial empirical evidence that this "one size fits all" approach results in sub optimal performance and argue for a novel framework that extends existing CDN architectures to provide programmatic control over the configuration options of the CDN serving stack. We designed ConfigTron a data driven framework that leverages data from all connections to identify their network characteristics and learn the optimal configuration parameters to improve end user performance. ConfigTron uses contextual multi arm bandit based learning algorithm to find optimal configurations in minimal time, enabling a content providers to systematically explore heterogeneous configurations while improving end user page load time by as much as 19 (upto 750ms) on median.


Feature Engineering for Mid Price Prediction with Deep Learning. Mid price movement prediction based on limit order book (LOB) data is a challenging task due to the complexity and dynamics of the LOB. So far, there have been very limited attempts for extracting relevant features based on LOB data. In this paper, we address this problem by designing a new set of handcrafted features and performing an extensive experimental evaluation on both liquid and illiquid stocks. More specifically, we implement a new set of econometrical features that capture statistical properties of the underlying securities for the task of mid price prediction. Moreover, we develop a new experimental protocol for online learning that treats the task as a multi objective optimization problem and predicts i) the direction of the next price movement and ii) the number of order book events that occur until the change takes place. In order to predict the mid price movement, the features are fed into nine different deep learning models based on multi layer perceptrons (MLP), convolutional neural networks (CNN) and long short term memory (LSTM) neural networks. The performance of the proposed method is then evaluated on liquid and illiquid stocks, which are based on TotalView ITCH US and Nordic stocks, respectively. For some stocks, results suggest that the correct choice of a feature set and a model can lead to the successful prediction of how long it takes to have a stock price movement.


On quadratic approximation for hyperquadratic continued fractions. We study quadratic approximations for two families of hyperquadratic continued fractions in the field of Laurent series over a finite field. As the first application, we give the answer to a question of the second author concerning Diophantine exponents for algebraic Laurent series. As the second application, we determine the degrees of these families in particular case.


GPU based parallelism for ASP solving. Answer Set Programming (ASP) has become, the paradigm of choice in the field of logic programming and non monotonic reasoning. Thanks to the availability of efficient solvers, ASP has been successfully employed in a large number of application domains. The term GPU computing indicates a recent programming paradigm aimed at enabling the use of modern parallel Graphical Processing Units (GPUs) for general purpose computing. In this paper we describe an approach to ASP solving that exploits GPU parallelism. The design of a GPU based solver poses various challenges due to the peculiarities of GPUs' software and hardware architectures and to the intrinsic nature of the satisfiability problem.


An O( log 3 2 n) Parallel Time Population Protocol for Majority with O( log n) States. In population protocols, the underlying distributed network consists of n nodes (or agents), denoted by V , and a scheduler that continuously selects uniformly random pairs of nodes to interact. When two nodes interact, their states are updated by applying a state transition function that depends only on the states of the two nodes prior to the interaction. The efficiency of a population protocol is measured in terms of both time (which is the number of interactions until the nodes collectively have a valid output) and the number of possible states of nodes used by the protocol. By convention, we consider the parallel time cost, which is the time divided by n . In this paper we consider the majority problem, where each node receives as input a color that is either black or white, and the goal is to have all of the nodes output the color that is the majority of the input colors. We design a population protocol that solves the majority problem in O( log 3 2 n) parallel time, both with high probability and in expectation, while using O( log n) states. Our protocol improves on a recent protocol of Berenbrink et al. that runs in O( log 5 3 n) parallel time, both with high probability and in expectation, using O( log n) states.


Dynamical convexity and closed orbits on symmetric spheres. The main theme of this paper is the dynamics of Reeb flows with symmetries on the standard contact sphere. We introduce the notion of strong dynamical convexity for contact forms invariant under a group action, supporting the standard contact structure, and prove that in dimension 2n 1 any such contact form satisfying a condition slightly weaker than strong dynamical convexity has at least n 1 simple closed Reeb orbits. For contact forms with antipodal symmetry, we prove that strong dynamical convexity is a consequence of ordinary convexity. In dimension five or greater, we construct examples of antipodally symmetric dynamically convex contact forms which are not strongly dynamically convex, and thus not contactomorphic to convex ones via a contactomorphism commuting with the antipodal map. Finally, we relax this condition on the contactomorphism furnishing a condition that has non empty C 1 interior.


Symmetric contact systems of segments, pseudotriangulations and inductive constructions for corresponding surface graphs. We characterise the quotient surface graphs arising from symmetric contact systems of line segments in the plane and also from symmetric pointed pseudotriangulations in the case where the group of symmetries is generated by a translation or a rotation of finite order. These results generalise well known results of Thomassen, in the case of line segments, and of Streinu and Haas et al., in the case of pseudotriangulations. Our main tool is a new inductive characterisation of the appropriate classes of surface graphs. We also discuss some consequences of our results in the area of geometric rigidity theory.


Abstraction based branch and bound approach to Q learning for hybrid optimal control. In this paper, we design a theoretical framework allowing to apply model predictive control on hybrid systems. For this, we develop a theory of approximate dynamic programming by leveraging the concept of alternating simulation. We show how to combine these notions in a branch and bound algorithm that can further refine the Q functions using Lagrangian duality. We illustrate the approach on a numerical example.