title: Optimal Bayesian clustering using non-negative matrix factorization

abstract: Bayesian model-based clustering is a widely applied procedure for discovering groups of related observations in a dataset. These approaches use Bayesian mixture models, estimated with MCMC, which provide posterior samples of the model parameters and clustering partition. While inference on model parameters is well established, inference on the clustering partition is less developed. A new method is developed for estimating the optimal partition from the pairwise posterior similarity matrix generated by a Bayesian cluster model. This approach uses non-negative matrix factorization (NMF) to provide a low-rank approximation to the similarity matrix. The factorization permits hard or soft partitions and is shown to perform better than several popular alternatives under a variety of penalty functions.


title: FFT Convolutions are Faster than Winograd on Modern CPUs, Here is Why

abstract: Winograd-based convolution has quickly gained traction as a preferred approach to implement convolutional neural networks (ConvNet) on various hardware platforms because it requires fewer floating point operations than FFT-based or direct convolutions. This paper compares three highly optimized implementations (regular FFT--, Gauss--FFT--, and Winograd--based convolutions) on modern multi-- and many--core CPUs. Although all three implementations employed the same optimizations for modern CPUs, our experimental results with two popular ConvNets (VGG and AlexNet) show that the FFT--based implementations generally outperform the Winograd--based approach, contrary to the popular belief. To understand the results, we use a Roofline performance model to analyze the three implementations in detail, by looking at each of their computation phases and by considering not only the number of floating point operations, but also the memory bandwidth and the cache sizes. The performance analysis explains why, and under what conditions, the FFT--based implementations outperform the Winograd--based one, on modern CPUs.


title: Modelling cell-cell collision and adhesion with the Filament Based Lamellipodium Model

abstract: We extend the live-cell motility Filament Based Lamellipodium Model to incorporate the forces exerted on the lamellipodium of the cells due to cell-cell collision and cadherin induced cell-cell adhesion. We take into account the nature of these forces via physical and biological constraints and modelling assumptions. We investigate the effect these new components have in the migration and morphology of the cells through particular experiments. We exhibit moreover the similarities between our simulated cells and HeLa cancer cells.


title: On Folding and Twisting (and whatknot): towards a characterization of workspaces in syntax

abstract: Syntactic theory has traditionally adopted a constructivist approach, in which a set of atomic elements are manipulated by combinatory operations to yield derived, complex elements. Syntactic structure is thus seen as the result or discrete recursive combinatorics over lexical items which get assembled into phrases, which are themselves combined to form sentences. This view is common to European and American structuralism (e.g., Benveniste, 1971; Hockett, 1958) and different incarnations of generative grammar, transformational and non-transformational (Chomsky, 1956, 1995; and Kaplan & Bresnan, 1982; Gazdar, 1982). Since at least Uriagereka (2002), there has been some attention paid to the fact that syntactic operations must apply somewhere, particularly when copying and movement operations are considered. Contemporary syntactic theory has thus somewhat acknowledged the importance of formalizing aspects of the spaces in which elements are manipulated, but it is still a vastly underexplored area. In this paper we explore the consequences of conceptualizing syntax as a set of topological operations applying over spaces rather than over discrete elements. We argue that there are empirical advantages in such a view for the treatment of long-distance dependencies and cross-derivational dependencies: constraints on possible configurations emerge from the dynamics of the system.


title: Excitation of positronium from 2S to 20P state

abstract: We report the observation of positronium excitation from the 2S to the 20P state. The Rydberg positronium atoms fly a distance of 40 mm before being field ionized and detected in a micro-channel plate. The time of flight can thus be measured and the velocity distribution of the atoms excited in the 2S state is reconstructed. This is used as an input to the model of the line-shape in order to properly take into account the second order Doppler shift which is the main systematic uncertainty in the 1S-2S measurements of positronium.


title: Energy-Efficient Mechanism for Smart Communication in Cellular Networks

abstract: Internet-of-things (IoT) is gaining popularity in recent times. Cellular network can be an appropriate solution for the IoT connectivity. With the exponential increase in the number of connected IoT devices, energy consumption is a bottleneck for wireless connectivity solutions. In this work, we proposed a scheme to decrease the power consumption of a base station of a cell. The maximum power consumption of a base station is associated with its data processing, routing, and transmission. Our scheme reduce the size of the transmission data which ultimately reduce the power consumption associated with it. Multiple data reduction techniques have been used for this purpose. System level simulations show the efficiency of our proposed system.


title: Inferring short-term volatility indicators from Bitcoin blockchain

abstract: In this paper, we study the possibility of inferring early warning indicators (EWIs) for periods of extreme bitcoin price volatility using features obtained from Bitcoin daily transaction graphs. We infer the low-dimensional representations of transaction graphs in the time period from 2012 to 2017 using Bitcoin blockchain, and demonstrate how these representations can be used to predict extreme price volatility events. Our EWI, which is obtained with a non-negative decomposition, contains more predictive information than those obtained with singular value decomposition or scalar value of the total Bitcoin transaction volume.

categories: q-fin.ST cs.CE cs.SI physics.soc-ph stat.ML


title: In-Edge AI: Intelligentizing Mobile Edge Computing, Caching and Communication by Federated Learning

abstract: Recently, along with the rapid development of mobile communication technology, edge computing theory and techniques have been attracting more and more attentions from global researchers and engineers, which can significantly bridge the capacity of cloud and requirement of devices by the network edges, and thus can accelerate the content deliveries and improve the quality of mobile services. In order to bring more intelligence to the edge systems, compared to traditional optimization methodology, and driven by the current deep learning techniques, we propose to integrate the Deep Reinforcement Learning techniques and Federated Learning framework with the mobile edge systems, for optimizing the mobile edge computing, caching and communication. And thus, we design the "In-Edge AI" framework in order to intelligently utilize the collaboration among devices and edge nodes to exchange the learning parameters for a better training and inference of the models, and thus to carry out dynamic system-level optimization and application-level enhancement while reducing the unnecessary system communication load. "In-Edge AI" is evaluated and proved to have near-optimal performance but relatively low overhead of learning, while the system is cognitive and adaptive to the mobile communication systems. Finally, we discuss several related challenges and opportunities for unveiling a promising upcoming future of "In-Edge AI".


title: Shouji: A Fast and Efficient Pre-Alignment Filter for Sequence Alignment

abstract: Motivation: The ability to generate massive amounts of sequencing data continues to overwhelm the processing capability of existing algorithms and compute infrastructures. In this work, we explore the use of hardware/software co-design and hardware acceleration to significantly reduce the execution time of short sequence alignment, a crucial step in analyzing sequenced genomes. We introduce Shouji, a highly-parallel and accurate pre-alignment filter that remarkably reduces the need for computationally-costly dynamic programming algorithms. The first key idea of our proposed pre-alignment filter is to provide high filtering accuracy by correctly detecting all common subsequences shared between two given sequences. The second key idea is to design a hardware accelerator that adopts modern FPGA (Field-Programmable Gate Array) architectures to further boost the performance of our algorithm. Results: Shouji significantly improves the accuracy of pre-alignment filtering by up to two orders of magnitude compared to the state-of-the-art pre-alignment filters, GateKeeper and SHD. Our FPGA-based accelerator is up to three orders of magnitude faster than the equivalent CPU implementation of Shouji. Using a single FPGA chip, we benchmark the benefits of integrating Shouji with five state-of-the-art sequence aligners, designed for different computing platforms. The addition of Shouji as a pre-alignment step reduces the execution time of the five state-of-the-art sequence aligners by up to 18.8x. Shouji can be adapted for any bioinformatics pipeline that performs sequence alignment for verification. Unlike most existing methods that aim to accelerate sequence alignment, Shouji does not sacrifice any of the aligner capabilities, as it does not modify or replace the alignment step. Availability: https://github.com/CMU-SAFARI/Shouji


title: Post-disaster 4G/5G Network Rehabilitation using Drones: Solving Battery and Backhaul Issues

abstract: Drone-based communications is a novel and attractive area of research in cellular networks. It provides several degrees of freedom in time (available on demand), space (mobile) and it can be used for multiple purposes (self-healing, offloading, coverage extension or disaster recovery). This is why the wide deployment of drone-based communications has the potential to be integrated in the 5G standard. In this paper, we utilize a grid of drones to provide cellular coverage to disaster-struck regions where the terrestrial infrastructure is totally damaged due to earthquake, flood, etc. We propose solutions for the most challenging issues facing drone networks which are limited battery energy and limited backhauling. Our proposed solution based mainly on using three types of drones; tethered backhaul drone (provides high capacity backhauling), untethered powering drone (provides on the fly battery charging) and untethered communication drone (provides cellular connectivity). Hence, an optimization problem is formulated to minimize the energy consumption of drones in addition to determining the placement of these drones and guaranteeing a minimum rate for the users. The simulation results show that we can provide unlimited cellular service to the disaster-affected region under certain conditions with a guaranteed minimum rate for each user.


title: Revenue Maximization in an Optical Router Node Using Multiple Wavelengths

abstract: In this paper, an optical router node with multiple wavelengths is considered. We introduce revenue for successful transmission and study the ensuing revenue maximization problem. We present an efficient and accurate heuristic procedure for solving the NP-hard revenue maximization problem and investigate the advantage offered by having multiple wavelengths.


title: Machine Learning for Forecasting Mid Price Movement using Limit Order Book Data

abstract: Forecasting the movements of stock prices is one the most challenging problems in financial markets analysis. In this paper, we use Machine Learning (ML) algorithms for the prediction of future price movements using limit order book data. Two different sets of features are combined and evaluated: handcrafted features based on the raw order book data and features extracted by ML algorithms, resulting in feature vectors with highly variant dimensionalities. Three classifiers are evaluated using combinations of these sets of features on two different evaluation setups and three prediction scenarios. Even though the large scale and high frequency nature of the limit order book poses several challenges, the scope of the conducted experiments and the significance of the experimental results indicate that Machine Learning highly befits this task carving the path towards future research in this field.


title: A Traffic-Aware Medium Access Control Mechanism for Energy-Efficient Wireless Network-on-Chip Architectures

abstract: Wireless interconnection has emerged as an energy efficient solution to the challenges of multi-hop communication over the wireline paths in conventional Networks-on-Chips (NoCs). However, to ensure the full benefits of this novel interconnect technology, design of simple, fair and efficient Medium Access Control (MAC) mechanism to grant access to the on-chip wireless communication channel is needed. Moreover, to adapt to the varying traffic demands from the applications running on a multicore environment, MAC mechanisms should dynamically adjust the transmission slots of the wireless interfaces (WIs). Such dynamic adjustment in transmission slots will result in improving the utilization of the wireless medium in a Wireless NoC (WiNoC). In this paper we present the design of two dynamic MAC mechanisms that adjust the transmission slots of the WIs based on predicted traffic demands and allow partial packet transfer. Through system level simulations, we demonstrate that the traffic aware MAC mechanisms are more energy efficient as well as capable of sustaining higher data bandwidth in WiNoCs.


title: Decentralized dynamic task allocation for UAVs with limited communication range

abstract: We present the Limited-range Online Routing Problem (LORP), which involves a team of Unmanned Aerial Vehicles (UAVs) with limited communication range that must autonomously coordinate to service task requests. We first show a general approach to cast this dynamic problem as a sequence of decentralized task allocation problems. Then we present two solutions both based on modeling the allocation task as a Markov Random Field to subsequently assess decisions by means of the decentralized Max-Sum algorithm. Our first solution assumes independence between requests, whereas our second solution also considers the UAVs' workloads. A thorough empirical evaluation shows that our workload-based solution consistently outperforms current state-of-the-art methods in a wide range of scenarios, lowering the average service time up to 16%. In the best-case scenario there is no gap between our decentralized solution and centralized techniques. In the worst-case scenario we manage to reduce by 25% the gap between current decentralized and centralized techniques. Thus, our solution becomes the method of choice for our problem.


title: Enabling Ultra-Low Delay Teleorchestras using Software Defined Networking

abstract: Ultra-low delay sensitive applications can afford delay only at the level of msec. An example of this application class are the Networked Music Performance (NMP) systems that describe a live music performance by geographically separate musicians over the Internet. The present work proposes a novel architecture for NMP systems, where the key-innovation is the close collaboration between the network and the application. Using SDN principles, the applications are enabled to adapt their internal audio signal processing, in order to cope with network delay increase. Thus, affordable end-to-end delay is provided to NMP users, even under considerable network congestion.


title: Mean Field Game Systems with Common Noise and Markovian Latent Processes

abstract: In many stochastic games stemming from financial models, the environment evolves with latent factors and there may be common noise across agents' states. Two classic examples are: (i) multi-agent trading on electronic exchanges, and (ii) systemic risk induced through inter-bank lending/borrowing. Moreover, agents' actions often affect the environment, and some agent's may be small while others large. Hence sub-population of agents may act as minor agents, while another class may act as major agents. To capture the essence of such problems, here, we introduce a general class of non-cooperative heterogeneous stochastic games with one major agent and a large population of minor agents where agents interact with an observed common process impacted by the mean field. A latent Markov chain and a latent Wiener process (common noise) modulate the common process, and agents cannot observe them. We use filtering techniques coupled with a convex analysis approach to (i) solve the mean field game limit of the problem, (ii) demonstrate that the best response strategies generate an $\epsilon$-Nash equilibrium for finite populations, and (iii) obtain explicit characterisations of the best response strategies.


title: Constructions and uses of incomplete pairwise balanced designs

abstract: We give explicit constructions for incomplete pairwise balanced designs IPBD$((v;w),K)$, or, equivalently, edge-decompositions of a difference of two cliques $K_v \setminus K_w$ into cliques whose sizes belong to the set $K$. Our constructions produce such designs whenever $v$ and $w$ satisfy the usual divisibility conditions, have ratio $v/w$ bounded away from the smallest value in $K$ minus one, say $v/w > k-1+\epsilon$, for $k =\min K$ and $\epsilon>0$, and are sufficiently large (depending on $K$ and $\epsilon$). As a consequence, some new results are obtained on many related designs, including class-uniformly resolvable designs, incomplete mutually orthogonal latin squares, and group divisible designs. We also include several other applications that illustrate the power of using IPBDs as `templates'.


title: Biological plausibility and stochasticity in scalable VO2 active memristor neurons

abstract: Neuromorphic networks of artificial neurons and synapses can solve computational hard problems with energy efficiencies unattainable for von Neumann architectures. For image processing, silicon neuromorphic processors outperform graphic processing units (GPUs) in energy efficiency by a large margin, but they deliver much lower chip-scale throughput. The performance-efficiency dilemma for silicon processors may not be overcome by Moore's law scaling of complementary metal-oxide-semiconductor (CMOS) field-effect transistors. Scalable and biomimetic active memristor neurons and passive memristor synapses form a self-sufficient basis for a transistorless neural network. However, previous demonstrations of memristor neurons only showed simple integrate-and-fire (I&F) behaviors and did not reveal the rich dynamics and computational complexity of biological neurons. Here we show that neurons built with nanoscale vanadium dioxide active memristors possess all three classes of excitability and most of the known biological neuronal dynamics, and are intrinsically stochastic. With the favorable size and power scaling, there is a path toward an all-memristor neuromorphic cortical computer.


title: Fusion hierarchies, $T$-systems and $Y$-systems for the $A_2^{(1)}$ models

abstract: The family of $A^{(1)}_2$ models on the square lattice includes a dilute loop model, a $15$-vertex model and, at roots of unity, a family of RSOS models. The fused transfer matrices of the general loop and vertex models are shown to satisfy $s\ell(3)$-type fusion hierarchies. We use these to derive explicit $T$- and $Y$-systems of functional equations. At roots of unity, we further derive closure identities for the functional relations and show that the universal $Y$-system closes finitely. The $A^{(1)}_2$ RSOS models are shown to satisfy the same functional and closure identities but with finite truncation.

categories: math-ph cond-mat.stat-mech hep-th math.MP


title: Current Status of Very-Large-Basis Hamiltonian Diagonalizations for Nuclear Physics

abstract: Today there are a plethora of many-body techniques for calculating nuclear wave functions and matrix elements. I review the status of that reliable workhorse, the interacting shell model, a.k.a. configuration-interaction methods, a.k.a. Hamiltonian diagonalization, and survey its advantages and disadvantages. With modern supercomputers one can tackle dimensions up to about 20 billion! I discuss how we got there and where we might go in the near future.


title: Suspended Load Path Tracking Control Using a Tilt-rotor UAV Based on Zonotopic State Estimation

abstract: This work addresses the problem of path tracking control of a suspended load using a tilt-rotor UAV. The main challenge in controlling this kind of system arises from the dynamic behavior imposed by the load, which is usually coupled to the UAV by means of a rope, adding unactuated degrees of freedom to the whole system. Furthermore, to perform the load transportation it is often needed the knowledge of the load position to accomplish the task. Since available sensors are commonly embedded in the mobile platform, information on the load position may not be directly available. To solve this problem in this work, initially, the kinematics of the multi-body mechanical system are formulated from the load's perspective, from which a detailed dynamic model is derived using the Euler-Lagrange approach, yielding a highly coupled, nonlinear state-space representation of the system, affine in the inputs, with the load's position and orientation directly represented by state variables. A zonotopic state estimator is proposed to solve the problem of estimating the load position and orientation, which is formulated based on sensors located at the aircraft, with different sampling times, and unknown-but-bounded measurement noise. To solve the path tracking problem, a discrete-time mixed $\mathcal{H}_2/\mathcal{H}_\infty$ controller with pole-placement constraints is designed with guaranteed time-response properties and robust to unmodeled dynamics, parametric uncertainties, and external disturbances. Results from numerical experiments, performed in a platform based on the Gazebo simulator and on a Computer Aided Design (CAD) model of the system, are presented to corroborate the performance of the zonotopic state estimator along with the designed controller.


title: On the supersymmetry invariance of flat supergravity with boundary

abstract: The supersymmetry invariance of flat supergravity (i.e., supergravity in the absence of any internal scale in the Lagrangian) in four dimensions on a manifold with non-trivial boundary is explored. Using a geometric approach we find that the supersymmetry invariance of the Lagrangian requires to add appropriate boundary terms. This is achieved by considering additional gauge fields to the boundary without modifying the bulk Lagrangian. We also construct an enlarged supergravity model from which, in the vanishing cosmological constant limit, flat supergravity with a non-trivial boundary emerges properly.


title: Ultralow Thermal Conductivity in a Two-Dimensional Material due to Surface Enhanced Resonant Bonding

abstract: Crystalline materials with ultralow thermal conductivity are highly desirable for thermoelectric applications. Many known crystalline materials with low thermal conductivity, including PbTe and Bi2Te3, possess a special kind of chemical bond called "resonant bond". Resonant bonds consist of superposition of degenerate bonding configurations that leads to structural instability, anomalous long-range interatomic interaction and soft optical phonons. These factors contribute to large lattice anharmonicity and strong phonon-phonon scattering, which result in low thermal conductivity. In this work, we use first-principles simulation to investigate the effect of resonant bonding in two dimensions (2D), where resonant bonds are in proximity to the surface. We find that the long-range interatomic interaction due to resonant bonding becomes more prominent in 2D due to reduced screening of the atomic-displacement-induced charge density distortion. To demonstrate this effect, we analyze the phonon properties of quasi-2D Bi2PbTe4 with an ultralow thermal conductivity of 0.74 W/mK at 300K. By comparing the interatomic force constants of quasi-2D Bi2PbTe4 and its bulk counterpart, and the properties of resonant bonds near the surface and in the bulk, we conclude that resonant bonds are significantly enhanced in reduced dimensions and are more effective in reducing the lattice thermal conductivity. Our results will provide new clues to searching for thermal insulators in low-dimensional materials.


title: Astrophysical Neutrinos with IceCube

abstract: The 1 km$^3$ IcCube neutrino observatory was built to find high-energy neutrinos that are associated with the sources of ultra-high energy cosmic rays. Its 5,160 optical sensors detect Cherenkov light from the charged particles produced when neutrinos interact in the ice. In this talk, I will describe the techniques that IceCube has used to search for astrophysical neutrinos. An emphasis will be given to diffuse neutrinos (not associated with a specific source), including analyses of contained events and energetic through-going neutrinos. At the end, I will discuss multi-messenger astronomy, and present an intriguing correlation between a high-energy IceCube neutrino and a blazar in outburst.


title: Task-Driven Estimation and Control via Information Bottlenecks

abstract: Our goal is to develop a principled and general algorithmic framework for task-driven estimation and control for robotic systems. State-of-the-art approaches for controlling robotic systems typically rely heavily on accurately estimating the full state of the robot (e.g., a running robot might estimate joint angles and velocities, torso state, and position relative to a goal). However, full state representations are often excessively rich for the specific task at hand and can lead to significant computational inefficiency and brittleness to errors in state estimation. In contrast, we present an approach that eschews such rich representations and seeks to create task-driven representations. The key technical insight is to leverage the theory of information bottlenecks}to formalize the notion of a "task-driven representation" in terms of information theoretic quantities that measure the minimality of a representation. We propose novel iterative algorithms for automatically synthesizing (offline) a task-driven representation (given in terms of a set of task-relevant variables (TRVs)) and a performant control policy that is a function of the TRVs. We present online algorithms for estimating the TRVs in order to apply the control policy. We demonstrate that our approach results in significant robustness to unmodeled measurement uncertainty both theoretically and via thorough simulation experiments including a spring-loaded inverted pendulum running to a goal location.


title: Can high-density human collective motion be forecasted by spatiotemporal fluctuations?

abstract: Concerts, protests, and sporting events are occurring with increasing frequency and magnitude. The extreme physical conditions common to these events are known to cause injuries and loss-of-life due to the emergence of collective motion such as crowd crush, turbulence, and density waves. Mathematical models of human crowds aimed at enhancing crowd safety by understanding these phenomena are developed with input from a variety of disciplines. However, model validation is challenged by a lack of high-quality empirical data and ethical constraints surrounding human crowd research. Consequently, generalized model-based approach for real-time monitoring/risk-assessment of crowd collective motion remains an open problem. Here, we take a model-free approach to crowd analysis and show that emergent collective motion can be forecasted directly from video data. We use mode analysis methods from material science and concepts from non-equilibrium physics to study footage of a human crowd at an Oasis rock concert. We analyze the attendees positional fluctuations during a period of crowd turbulence to predict the spatial patterns of an emergent human density wave. In addition to predicting spatial patterns of collective motion, we also identify and measure temporal patterns that precede the density wave and forecast its appearance by 1~s. Looking ahead, widening this forecasting window beyond 1~s will enable new computer vision technologies for real-time risk-assessment of emergent human collective motion.


title: Tiling Billards on Triangle Tilings, and Interval Exchange Transformations

abstract: We consider the dynamics of light rays in triangle tilings where triangles are transparent and adjacent triangles have equal but opposite indices of refraction. We find that the behavior of a trajectory on a triangle tiling is described by an orientation-reversing three-interval exchange transformation on the circle, and that the behavior of all the trajectories on a given triangle tiling is described by a polygon exchange transformation. We show that, for a particular choice of triangle tiling, certain trajectories approach the Rauzy fractal, under rescaling.


title: Lepto-philic 2-HDM + singlet scalar portal induced fermionic dark matter

abstract: We explore the possibility that the discrepancy in the observed anomalous magnetic moment of the muon $\Delta a _\mu$ and the predicted relic abundance of Dark Matter by Planck data, can be explained in a lepto-philic 2-HDM augmented by a real SM singlet scalar of mass $\sim$ 10-80 GeV. We constrain the model from the observed Higgs Decay width at LHC, LEP searches for low mass exotic scalars and anomalous magnetic moment of an electron $\Delta a_e$. This constrained light singlet scalar serves as a portal for the fermionic Dark Matter, which contributes to the required relic density of the universe. A large region of model parameter space is found to be consistent with the present observations from the Direct and Indirect DM detection experiments.


title: Learning Quickly to Plan Quickly Using Modular Meta-Learning

abstract: Multi-object manipulation problems in continuous state and action spaces can be solved by planners that search over sampled values for the continuous parameters of operators. The efficiency of these planners depends critically on the effectiveness of the samplers used, but effective sampling in turn depends on details of the robot, environment, and task. Our strategy is to learn functions called "specializers" that generate values for continuous operator parameters, given a state description and values for the discrete parameters. Rather than trying to learn a single specializer for each operator from large amounts of data on a single task, we take a modular meta-learning approach. We train on multiple tasks and learn a variety of specializers that, on a new task, can be quickly adapted using relatively little data -- thus, our system "learns quickly to plan quickly" using these specializers. We validate our approach experimentally in simulated 3D pick-and-place tasks with continuous state and action spaces. Visit http://tinyurl.com/chitnis-icra-19 for a supplementary video.


title: Statistical analysis of astro-geodetic data through principal component analysis, linear modelling and bootstrap based inference

abstract: The paper demonstrates the application of statistical based methodology for the analysis of the vertical deviation angle. The studied data set contains astro-geodetic observations. The Principal Component Analysis and the Multiple Linear Regression models are embedded within a bootstrap procedure, in order to overcome the difficulties related to data correlation, while taking advantage of all the information provided. The methodology is applied on real data. The obtained results indicate that the pressure, the temperature and the humidity are variables that may influence the measure of the vertical deviation.


title: Teaching Social Behavior through Human Reinforcement for Ad hoc Teamwork -The STAR Framework

abstract: As AI technology continues to develop, more and more agents will become capable of long term autonomy alongside people. Thus, a recent line of research has studied the problem of teaching autonomous agents the concept of ethics and human social norms. Most existing work considers the case of an individual agent attempting to learn a predefined set of rules. In reality, however, social norms are not always pre-defined and are very difficult to represent algorithmically. Moreover, the basic idea behind the social norms concept is ensuring that one's actions do not negatively influence others' utilities, which is inherently a multiagent concept. Thus, here we investigate a way to teach agents, as a team, how to act according to human social norms. In this research, we introduce the STAR framework used to teach an ad hoc team of agents to act in accordance with human social norms. Using a hybrid team (agents and people), when taking an action considered to be socially unacceptable, the agents receive negative feedback from the human teammate(s) who has(have) an awareness of the team's norms. We view STAR as an important step towards teaching agents to act more consistently with respect to human morality.


title: Triple correlation and long gaps in the spectrum of flat tori

abstract: We evaluate the triple correlation of eigenvalues of the Laplacian on generic flat tori in an averaged sense. As two consequence we show that (a) the limit inferior (resp. limit superior) of the triple correlation is almost surely at most (resp. at least) Poissonian, and that (b) almost all flat tori contain infinitely many gaps in their spectrum that are at least 2.006 times longer than the average gap. The significance of the constant 2.006 lies in the fact that there exist sequences with Poissonian pair correlation and with consecutive spacings bounded uniformly from above by 2, as we also prove in this paper. Thus our result goes beyond what can be deduced solely from the Poissonian behavior of the pair correlation.


title: Uncertainty Aware AI ML: Why and How

abstract: This paper argues the need for research to realize uncertainty-aware artificial intelligence and machine learning (AI\&ML) systems for decision support by describing a number of motivating scenarios. Furthermore, the paper defines uncertainty-awareness and lays out the challenges along with surveying some promising research directions. A theoretical demonstration illustrates how two emerging uncertainty-aware ML and AI technologies could be integrated and be of value for a route planning operation.


title: Differential Dynamic Programming on Lie Groups: Derivation, Convergence Analysis and Numerical Results

abstract: We develop a discrete-time optimal control framework for systems evolving on Lie groups. Our work generalizes the original Differential Dynamic Programming method, by employing a coordinate-free, Lie-theoretic approach for its derivation. A key element lies, specifically, in the use of quadratic expansion schemes for cost functions and dynamics defined on manifolds. The obtained algorithm iteratively optimizes local approximations of the control problem, until reaching a (sub)optimal solution. On the theoretical side, we also study the conditions under which convergence is attained. Details about the behavior and implementation of our method are provided through a simulated example on T SO(3).


title: Absence of singular continuous spectrum for perturbed discrete Schr\"odinger operators

abstract: We show that the spectral measure of discrete Schr\"odinger operators $ (Hu)(n)= u({n+1})+u({n-1})+V(n)u(n)$ does not have singular continuous component if the potential $V(n)=O(n^{-1})$.


title: Meromorphic Painlev\'e III transcendents and the Joukowski correspondence

abstract: We study a twistor correspondence based on the Joukowski map reduced from one for stationary-axisymmetric self-dual Yang-Mills and adapt it to the Painlev\'e III equation. A natural condition on the geometry (axissimplicity) leads to solutions that are meromorphic at the fixed singularity at the origin. We show that it also implies a quantisation condition for the parameter in the equation. From the point of view of generalized monodromy data, the condition is equivalent to triviality of the Stokes matrices and half-integral exponents of formal monodromy. We obtain canonically defined representations in terms of a Birkhoff factorization whose entries are related to the data at the origin and the Painlev\'e constants.


title: Emission of intense resonant radiation by dispersion-managed Kerr cavity solitons

abstract: We report on an experimental and numerical study of temporal Kerr cavity soliton dynamics in dispersion-managed fiber ring resonators. We find that dispersion management can significantly magnify the Kelly-like resonant radiation sidebands emitted by the solitons. Because of the underlying phase-matching conditions, the sideband amplitudes tend to increase with increasing pump-cavity detuning, ultimately limiting the range of detunings over which the solitons can exist. Our experimental findings show excellent agreement with numerical simulations.


title: Closeness of Solutions for Singularly Perturbed Systems via Averaging

abstract: This paper studies the behavior of singularly perturbed nonlinear differential equations with boundary-layer solutions that do not necessarily converge to an equilibrium. Using the average of the fast variable and assuming the boundary layer solutions converge to a bounded set, results on the closeness of solutions of the singularly perturbed system to the solutions of the reduced average and boundary layer systems over a finite time interval are presented. The closeness of solutions error is shown to be of order O(\sqrt(\epsilon)), where \epsilon is the perturbation parameter.


title: Probabilistic Logic Programming with Beta-Distributed Random Variables

abstract: We enable aProbLog---a probabilistic logical programming approach---to reason in presence of uncertain probabilities represented as Beta-distributed random variables. We achieve the same performance of state-of-the-art algorithms for highly specified and engineered domains, while simultaneously we maintain the flexibility offered by aProbLog in handling complex relational domains. Our motivation is that faithfully capturing the distribution of probabilities is necessary to compute an expected utility for effective decision making under uncertainty: unfortunately, these probability distributions can be highly uncertain due to sparse data. To understand and accurately manipulate such probability distributions we need a well-defined theoretical framework that is provided by the Beta distribution, which specifies a distribution of probabilities representing all the possible values of a probability when the exact value is unknown.


title: Predicting the Argumenthood of English Prepositional Phrases

abstract: Distinguishing between arguments and adjuncts of a verb is a longstanding, nontrivial problem. In natural language processing, argumenthood information is important in tasks such as semantic role labeling (SRL) and prepositional phrase (PP) attachment disambiguation. In theoretical linguistics, many diagnostic tests for argumenthood exist but they often yield conflicting and potentially gradient results. This is especially the case for syntactically oblique items such as PPs. We propose two PP argumenthood prediction tasks branching from these two motivations: (1) binary argument-adjunct classification of PPs in VerbNet, and (2) gradient argumenthood prediction using human judgments as gold standard, and report results from prediction models that use pretrained word embeddings and other linguistically informed features. Our best results on each task are (1) $acc.=0.955$, $F_1=0.954$ (ELMo+BiLSTM) and (2) Pearson's $r=0.624$ (word2vec+MLP). Furthermore, we demonstrate the utility of argumenthood prediction in improving sentence representations via performance gains on SRL when a sentence encoder is pretrained with our tasks.


title: A non-iterative method for optimizing linear functions in convex linearly constrained spaces

abstract: This document introduces a strategy to solve linear optimization problems. The strategy is based on the bounding condition each constraint produces on each one of the problem's dimension. The solution of a linear optimization problem is located at the intersection of the constraints defining the extreme vertex. By identifying the constraints that limit the growth of the objective function value, we formulate linear equations system leading to the optimization problem's solution.The most complex operation of the algorithm is the inversion of a matrix sized by the number of dimensions of the problem. Therefore, the algorithm's complexity is comparable to the corresponding to the classical Simplex method and the more recently developed Linear Programming algorithms. However, the algorithm offers the advantage of being non-iterative.


title: The asymptotics of purely atomic one-dimensional Levy approximations

abstract: For arbitrary Borel probability measures on the real line, necessary and sufficient conditions are presented that characterize best purely atomic approximations relative to the classical Levy probability metric, given any number of atoms, and allowing for additional constraints regarding locations or weights of atoms. The precise asymptotics (as the number of atoms goes to infinity) of the approximation error is identified for the important special cases of best uniform (i.e., all atoms having equal weight) and best (unconstrained) approximations, respectively. When compared to similar results known for other probability metrics, the results for Levy approximations are more complete and require fewer assumptions.


title: Error Analysis of the Stochastic Linear Feedback Particle Filter

abstract: This paper is concerned with the convergence and long-term stability analysis of the feedback particle filter (FPF) algorithm. The FPF is an interacting system of $N$ particles where the interaction is designed such that the empirical distribution of the particles approximates the posterior distribution. It is known that in the mean-field limit ($N=\infty$), the distribution of the particles is equal to the posterior distribution. However little is known about the convergence to the mean-field limit. In this paper, we consider the FPF algorithm for the linear Gaussian setting. In this setting, the algorithm is similar to the ensemble Kalman-Bucy filter algorithm. Although these algorithms have been numerically evaluated and widely used in applications, their convergence and long-term stability analysis remains an active area of research. In this paper, we show that, (i) the mean-field limit is well-defined with a unique strong solution; (ii) the mean-field process is stable with respect to the initial condition; (iii) we provide conditions such that the finite-$N$ system is long term stable and we obtain some mean-squared error estimates that are uniform in time.


title: Solving Large Extensive-Form Games with Strategy Constraints

abstract: Extensive-form games are a common model for multiagent interactions with imperfect information. In two-player zero-sum games, the typical solution concept is a Nash equilibrium over the unconstrained strategy set for each player. In many situations, however, we would like to constrain the set of possible strategies. For example, constraints are a natural way to model limited resources, risk mitigation, safety, consistency with past observations of behavior, or other secondary objectives for an agent. In small games, optimal strategies under linear constraints can be found by solving a linear program; however, state-of-the-art algorithms for solving large games cannot handle general constraints. In this work we introduce a generalized form of Counterfactual Regret Minimization that provably finds optimal strategies under any feasible set of convex constraints. We demonstrate the effectiveness of our algorithm for finding strategies that mitigate risk in security games, and for opponent modeling in poker games when given only partial observations of private information.


title: Towards the determination of heavy-quark transport coefficients in quark-gluon plasma

abstract: Several transport models have been employed in recent years to analyze heavy-flavor meson spectra in high-energy heavy-ion collisions. Heavy-quark transport coefficients extracted from these models with their default parameters vary, however, by up to a factor of 5 at high momenta. To investigate the origin of this large theoretical uncertainty, a systematic comparison of heavy-quark transport coefficients is carried out between various transport models. Within a common scheme devised for the nuclear modification factor of charm quarks in a brick medium of a quark-gluon plasma, the systematic uncertainty of the extracted drag coefficient among these models is shown to be reduced to a factor of 2, which can be viewed as the smallest intrinsic systematical error band achievable at present time. This indicates the importance of a realistic hydrodynamic evolution constrained by bulk hadron spectra and of heavy-quark hadronization for understanding the final heavy-flavor hadron spectra and extracting heavy-quark drag coefficient. The transverse transport coefficient is less constrained due to the influence of the underlying mechanism for heavy-quark medium interaction. Additional constraints on transport models such as energy loss fluctuation and transverse-momentum broadening can further reduce theoretical uncertainties in the extracted transport coefficients.


title: Large-Scale Video Classification with Feature Space Augmentation coupled with Learned Label Relations and Ensembling

abstract: This paper presents the Axon AI's solution to the 2nd YouTube-8M Video Understanding Challenge, achieving the final global average precision (GAP) of 88.733% on the private test set (ranked 3rd among 394 teams, not considering the model size constraint), and 87.287% using a model that meets size requirement. Two sets of 7 individual models belonging to 3 different families were trained separately. Then, the inference results on a training data were aggregated from these multiple models and fed to train a compact model that meets the model size requirement. In order to further improve performance we explored and employed data over/sub-sampling in feature space, an additional regularization term during training exploiting label relationship, and learned weights for ensembling different individual models.


title: 3D Move to See: Multi-perspective visual servoing for improving object views with semantic segmentation

abstract: In this paper, we present a new approach to visual servoing for robotics, referred to as 3D Move to See (3DMTS), based on the principle of finding the next best view using a 3D camera array and a robotic manipulator to obtain multiple samples of the scene from different perspectives. The method uses semantic vision and an objective function applied to each perspective to sample a gradient representing the direction of the next best view. The method is demonstrated within simulation and on a real robotic platform containing a custom 3D camera array for the challenging scenario of robotic harvesting in a highly occluded and unstructured environment. It was shown on a real robotic platform that by moving the end effector using the gradient of an objective function leads to a locally optimal view of the object of interest, even amongst occlusions. The overall performance of the 3DMTS method obtained a mean increase in target size by 29.3% compared to a baseline method using a single RGB-D camera, which obtained 9.17%. The results demonstrate qualitatively and quantitatively that the 3DMTS method performed better in most scenarios, and yielded three times the target size compared to the baseline method. The increased target size in the final view will improve the detection of key features of the object of interest for further manipulation, such as grasping and harvesting.


title: Modalities, Cohesion, and Information Flow

abstract: It is informally understood that the purpose of modal type constructors in programming calculi is to control the flow of information between types. In order to lend rigorous support to this idea, we study the category of classified sets, a variant of a denotational semantics for information flow proposed by Abadi et al. We use classified sets to prove multiple noninterference theorems for modalities of a monadic and comonadic flavour. The common machinery behind our theorems stems from the the fact that classified sets are a (weak) model of Lawvere's theory of axiomatic cohesion. In the process, we show how cohesion can be used for reasoning about multi-modal settings. This leads to the conclusion that cohesion is a particularly useful setting for the study of both information flow, but also modalities in type theory and programming languages at large.


title: MOA-2016-BLG-319Lb: Microlensing Planet Subject to Rare Minor-Image Perturbation Degeneracy in Determining Planet Parameter

abstract: We present the analysis of the planetary microlensing event MOA-2016-BLG-319. The event light curve is characterized by a brief ($\sim 3$ days) anomaly near the peak produced by minor-image perturbations. From modeling, we find two distinct solutions that describe the observed light curve almost equally well. From the investigation of the lens-system configurations, we find that the confusion in the lensing solution is caused by the degeneracy between the two solutions resulting from the source passages on different sides of the planetary caustic. These degeneracies can be severe for major-image perturbations but it is known that they are considerably less severe for minor-image perturbations. From the comparison of the lens-system configuration with those of two previously discovered planetary events, for which similar degeneracies were reported, we find that the degeneracies are caused by the special source trajectories that passed the star-planet axes at approximately right angles. By conducting a Bayesian analysis, it is estimated that the lens is a planetary system in which a giant planet with a mass $M_{\rm p}=0.62^{+1.16}_{-0.33}~M_{\rm J}$ ($0.65^{+1.21}_{-0.35}~M_{\rm J}$) is orbiting a low-mass M-dwarf host with a mass $M_{\rm h}=0.15^{+0.28}_{-0.08}~M_\odot$. Here the planet masses in and out of the parentheses represent the masses for the individual degenerate solutions. The projected host-planet separations are $a_\perp\sim 0.95$ au and $\sim 1.05$ au for the two solutions. The identified degeneracy indicates the need to check similar degeneracies in future analyses of planetary lensing events with minor-image perturbations.


title: Measuring K$^0_{\rm S}$K$^{\rm{\pm}}$ interactions using pp collisions at $\sqrt{s}=7$ TeV

abstract: We present the first measurements of femtoscopic correlations between the K$^0_{\rm S}$ and K$^{\rm \pm}$ particles in pp collisions at $\sqrt{s}=7$ TeV measured by the ALICE experiment. The observed femtoscopic correlations are consistent with final-state interactions proceeding solely via the $a_0(980)$ resonance. The extracted kaon source radius and correlation strength parameters for K$^0_{\rm S}$K$^{\rm -}$ are found to be equal within the experimental uncertainties to those for K$^0_{\rm S}$K$^{\rm +}$. Results of the present study are compared with those from identical-kaon femtoscopic studies also performed with pp collisions at $\sqrt{s}=7$ TeV by ALICE and with a K$^0_{\rm S}$K$^{\rm \pm}$ measurement in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV. Combined with the Pb-Pb results, our pp analysis is found to be compatible with the interpretation of the $a_0(980)$ having a tetraquark structure instead of that of a diquark.


title: Thermomagnetic stability in isotropic type-II superconductors under multicomponent magnetic fields

abstract: The goal of this research is the study of the thermomagnetic consequences in isotropic type-II superconductors, subjected to multi-component magnetic fields $\mathbf{H}_a = H_{ay}\hat{y}+H_{az}\hat{z}$, because the instability field $\mathbf{H}_{fi}$ is closely related with a flux jump occurrence. At the critical-state model framework, once the Lorentz $\mathbf{F}_L$ and pinning forces $\mathbf{F}_P$ are at equilibrium, the current density reaches a critical value $j_c$ and a stationary magnetic induction distribution $\mathbf{B}$ is established. The equilibrium of forces is analytically solved considering that the pinning force is mainly affected by temperature increments; the energy dissipation is incorporated throughout the heat equation at the adiabatic regime. The theory is able to obtain the instability field according to the thermal bath and applied field values; moreover, it provides of instability field branches comprising both partial and full penetrates states. With this information is possible to construct a field-temperature map. The results are compared with already published experimental data, finding a qualitatively agreement between them. This theoretical study works with a first order perturbation, then the perturbation presents a periodical behavior along the thickness direction; considering this environment it is constructed the magnetic induction distributions which resemble flexible cantilever structures.


title: Integrated Analysis of Performance and Resource of Large-Scale Quantum Computing

abstract: To see the feasibility of a large-scale quantum computing, it is required to accurately analyze the performance and the quantum resource. However, most of the analysis reported so far have focused on the statistical examination, i.e., simply calculating the performance and resource based on individual data, and even worse usually only a few components have been considered. In this work, to achieve more exact analysis, we propose an integrated analysis method for a practical quantum computing model with three components (\textit{algorithm}, \textit{error correction} and \textit{device}) under a realistic quantum computer system architecture. To implement the above method, we develop a quantum computing framework composed of three functional layers: compile, system and building block. This framework can support, for the first time, the mapping of quantum algorithm from physical qubit level to system architecture level with a given fault-tolerant scheme. Therefore, the proposed method can measure the effect of dynamic situation when the quantum computer practically runs. By using our method, we found that Shor algorithm to factorize 512-bit integer requires $8.78\times 10^ 5$ hours. We also show how the proposed method can be used for analyzing optimal concatenation level and code distance of fault-tolerant quantum computing.


title: Charge Order in the Holstein Model on a Honeycomb Lattice

abstract: The effect of electron-electron interactions on Dirac fermions, and the possibility of an intervening spin liquid phase between the semi-metal and antiferromagnetic (AF) regimes, has been a focus of intense quantum simulation effort over the last five years. We use determinant quantum Monte Carlo (DQMC) to study the Holstein model on a Honeycomb lattice and explore the role of electron-phonon interactions on Dirac fermions. We show that they give rise to charge density wave (CDW) order, and present evidence that this occurs only above a finite critical interaction strength. We evaluate the temperature for the transition into the CDW which, unlike the AF transition, can occur at finite values owing to the discrete nature of the broken symmetry.


title: Charge-Density-Wave Transitions of Dirac Fermions Coupled to Phonons

abstract: The spontaneous generation of charge-density-wave order in a Dirac fermion system via the natural mechanism of electron-phonon coupling is studied in the framework of the Holstein model on the honeycomb lattice. Using two independent and unbiased quantum Monte Carlo methods, the phase diagram as a function of temperature and coupling strength is determined. It features a quantum critical point as well as a line of thermal critical points. Finite-size scaling appears consistent with fermionic Gross-Neveu-Ising universality for the quantum phase transition, and bosonic Ising universality for the thermal phase transition. The critical temperature has a maximum at intermediate couplings. Our findings motivate experimental efforts to identify or engineer Dirac systems with sufficiently strong and tunable electron-phonon coupling.

categories: cond-mat.str-el cond-mat.supr-con physics.atom-ph quant-ph


title: Automatic Rule Learning for Autonomous Driving Using Semantic Memory

abstract: This paper presents a novel approach for automatic rule learning applicable to an autonomous driving system using real driving data.


title: Subgroup Identification Using the personalized Package

abstract: A plethora of disparate statistical methods have been proposed for subgroup identification to help tailor treatment decisions for patients. However a majority of them do not have corresponding R packages and the few that do pertain to particular statistical methods or provide little means of evaluating whether meaningful subgroups have been found. Recently, the work of Chen, Tian, Cai, and Yu (2017) unified many of these subgroup identification methods into one general, consistent framework. The goal of the personalized package is to provide a corresponding unified software framework for subgroup identification analyses that provides not only estimation of subgroups, but evaluation of treatment effects within estimated subgroups. The personalized package allows for a variety of subgroup identification methods for many types of outcomes commonly encountered in medical settings. The package is built to incorporate the entire subgroup identification analysis pipeline including propensity score diagnostics, subgroup estimation, analysis of the treatment effects within subgroups, and evaluation of identified subgroups. In this framework, different methods can be accessed with little change in the analysis code. Similarly, new methods can easily be incorporated into the package. Besides familiar statistical models, the package also allows flexible machine learning tools to be leveraged in subgroup identification. Further estimation improvements can be obtained via efficiency augmentation.


title: A study of X-ray emission of galaxies hosting molecular outflows (MOX sample)

abstract: We have carried out an extensive X-ray spectral analysis of a sample of galaxies exhibiting molecular outflows (MOX sample), to characterize the X-ray properties and investigate the effect of AGN on the dynamical properties of the molecular outflows. We find that the X-ray bolometric correction $(L_{2-10\rm keV}/L_{\rm AGN})$ of these sources ranges from $\sim10^{-4.5}$ to $10^{-0.5}$, with $\sim 70\%$ of the sources below $10^{-2}$, implying a weak X-ray emission relative to the AGN bolometric luminosity ($L_{\rm AGN}$). However, the upper limit on the $2-10\rm keV$ luminosity ($L_{\rm 2-10 keV, \,12\mu m}$) obtained from $12\mu$m flux, following the correlation derived by Asmus et al., are $\sim 0.5-3$ orders of magnitude larger than the $L_{2-10\rm keV}$ values estimated using X-ray spectroscopy, implying a possibility that the MOX sources host normal AGN (not X-ray weak), and their X-ray spectra are extremely obscured. We find that both $L_{2-10\rm keV}$, and $L_{\rm AGN}$ correlates strongly with the molecular outflow velocity as well as the mass outflow rates ($\dot{M}_{\rm out}$), implying that the central AGN plays an important role in driving these massive outflows. However, we also find statistically significant positive correlations between the starburst emission and MO mass outflow rate, $L_{\rm Starburst}$ vs $\dot{M}_{\rm out}$, and $L_{0.6-2\rm keV}$ vs $\dot{M}_{\rm out}$, which implies that starbursts can generate and drive the molecular outflows. The correlations of MO velocity and $\dot{M}_{\rm out}$ with AGN luminosities are found to be stronger compared to those with the starburst luminosities. We conclude that both starbursts and AGN play crucial role in driving the large scale MO.


title: A Unified Framework for the Teleoperation of Surgical Robots in Constrained Workspaces

abstract: In adult laparoscopy, robot-aided surgery is a reality in thousands of operating rooms worldwide, owing to the increased dexterity provided by the robotic tools. Many robots and robot control techniques have been developed to aid in more challenging scenarios, such as pediatric surgery and microsurgery. However, the prevalence of case-specific solutions, particularly those focused on non-redundant robots, reduces the reproducibility of the initial results in more challenging scenarios. In this paper, we propose a general framework for the control of surgical robotics in constrained workspaces under teleoperation, regardless of the robot geometry. Our technique is divided into a slave-side constrained optimization algorithm, which provides virtual fixtures, and with Cartesian impedance on the master side to provide force feedback. Experiments with two robotic systems, one redundant and one non-redundant, show that smooth teleoperation can be achieved in adult laparoscopy and infant surgery.


title: Two-level systems and growth-induced metastability in hydrogenated amorphous silicon

abstract: Specific heat measurements from 2 to 300 K of hydrogenated amorphous silicon prepared by hot-wire chemical vapor deposition show a large excess specific heat at low temperature, significantly larger than the Debye specific heat calculated from the sound velocity. The as-prepared films have a Schottky anomaly that is associated with metastable hydrogen in the amorphous network, as well as large linear and excess cubic term commonly associated with tunneling two-level systems in amorphous solids. Annealing at 200 {\deg}C, a temperature that enables hydrogen mobility but not evaporation, irreversibly reduces the heat capacity, eliminating the Schottky anomaly and leaving a reduced linear heat capacity. A non-monotonic dependence on growth temperature and H content is observed in all effects, except for sound velocity, which suggests that the tunneling two-level systems and the Schottky anomaly are associated with atomic hydrogen and require low density regions to form, while sound velocity is associated with the silicon network and improves with increasing growth temperature.


title: On the existence and multiplicity of solutions to fractional Lane-Emden elliptic systems involving measures

abstract: We study positive solutions to the fractional Lane-Emden system \begin{equation*} \tag{S}\label{S} \left\{ \begin{aligned} (-\Delta)^s u &= v^p+\mu \quad &&\text{in } \Omega \\ (-\Delta)^s v &= u^q+\nu \quad &&\text{in } \Omega\\ u = v &= 0 \quad &&\text{in } \Omega^c={\mathbb R}^N \setminus \Omega, \end{aligned} \right. \end{equation*} where $\Omega$ is a $C^2$ bounded domains in ${\mathbb R}^N$, $s\in(0,1)$, $N>2s$, $p>0$, $q>0$ and $\mu,\, \nu$ are positive measures in $\Omega$. We prove the existence of the minimal positive solution of the above system under a smallness condition on the total mass of $\mu$ and $\nu$. Furthermore, if $p,q \in (1,\frac{N+s}{N-s})$ and $0 \leq \mu,\, \nu\in L^r(\Omega)$ for some $r>\frac{N}{2s}$ then we show the existence of at least two positive solutions of the above system. We also discuss the regularity of the solutions.


title: Simple Local Computation Algorithms for the General Lovasz Local Lemma

abstract: We consider the task of designing Local Computation Algorithms (LCA) for applications of the Lov\'{a}sz Local Lemma (LLL). LCA is a class of sublinear algorithms proposed by Rubinfeld et al.~\cite{Ronitt} that have received a lot of attention in recent years. The LLL is an existential, sufficient condition for a collection of sets to have non-empty intersection (in applications, often, each set comprises all objects having a certain property). The ground-breaking algorithm of Moser and Tardos~\cite{MT} made the LLL fully constructive, following earlier results by Beck~\cite{beck_lll} and Alon~\cite{alon_lll} giving algorithms under significantly stronger LLL-like conditions. LCAs under those stronger conditions were given in~\cite{Ronitt}, where it was asked if the Moser-Tardos algorithm can be used to design LCAs under the standard LLL condition. The main contribution of this paper is to answer this question affirmatively. In fact, our techniques yield LCAs for settings beyond the standard LLL condition.


title: Hydrodynamic Limit for the SSEP with a Slow Membrane

abstract: In this paper we consider a symmetric simple exclusion process (SSEP) on the $d$-dimensional discrete torus $\mathbb{T}^d_N$ with a spatial non-homogeneity given by a slow membrane. The slow membrane is defined here as the boundary of a smooth simple connected region $\Lambda$ on the continuous $d$-dimensional torus $\mathbb{T}^d$. In this setting, bonds crossing the membrane have jump rate $\alpha/N^\beta$ and all other bonds have jump rate one, where $\alpha>0$, $\beta\in[0,\infty]$, and $N\in \mathbb{N}$ is the scaling parameter. In the diffusive scaling we prove that the hydrodynamic limit presents a dynamical phase transition, that is, it depends on the regime of $\beta$. For $\beta\in[0,1)$, the hydrodynamic equation is given by the usual heat equation on the continuous torus, meaning that the slow membrane has no effect in the limit. For $\beta\in(1,\infty]$, the hydrodynamic equation is the heat equation with Neumann boundary conditions, meaning that the slow membrane $\partial \Lambda$ divides $\mathbb{T}^d$ into two isolated regions $\Lambda$ and $\Lambda^\complement$. And for the critical value $\beta=1$, the hydrodynamic equation is the heat equation with certain Robin boundary conditions related to the Fick's Law.


title: Cloud-Based Approximate Constrained Shortest Distance Queries Over Encrypted Graphs With Privacy Protection

abstract: Constrained shortest distance (CSD) querying is one of the fundamental graph query primitives, which finds the shortest distance from an origin to a destination in a graph with a constraint that the total cost does not exceed a given threshold. CSD querying has a wide range of applications, such as routing in telecommunications and transportation. With an increasing prevalence of cloud computing paradigm, graph owners desire to outsource their graphs to cloud servers. In order to protect sensitive information, these graphs are usually encrypted before being outsourced to the cloud. This, however, imposes a great challenge to CSD querying over encrypted graphs. Since performing constraint filtering is an intractable task, existing work mainly focuses on unconstrained shortest distance queries. CSD querying over encrypted graphs remains an open research problem. In this paper, we propose Connor, a novel graph encryption scheme that enables approximate CSD querying. Connor is built based on an efficient, tree-based ciphertext comparison protocol, and makes use of symmetric-key primitives and the somewhat homomorphic encryption, making it computationally efficient. Using Connor, a graph owner can first encrypt privacy-sensitive graphs and then outsource them to the cloud server, achieving the necessary privacy without losing the ability of querying. Extensive experiments with real-world datasets demonstrate the effectiveness and efficiency of the proposed graph encryption scheme.


title: The quantization of the standard triadic Cantor distribution

abstract: The quantization scheme in probability theory deals with finding a best approximation of a given probability distribution by a probability distribution that is supported on finitely many points. For a given $k\geq 2$, let $\{S_j : 1\leq j\leq k\}$ be a set of $k$ contractive similarity mappings such that $S_j(x)=\frac 1 {2k-1} x +\frac{2 (j-1)} {2k-1}$ for all $x\in \mathbb R$, and let $P= \frac 1 k \sum_{j=1}^kP\circ S_j^{-1}$. Then, $P$ is a unique Borel probability measure on $\mathbb R$ such that $P$ has support the Cantor set generated by the similarity mappings $S_j$ for $1\leq j\leq k$. In this paper, for the probability measure $P$, when $k=3$, we investigate the optimal sets of $n$-means and the $n$th quantization errors for all $n\geq 2$. We further show that the quantization coefficient does not exist though the quantization dimension exists.


title: Secure Phrase Search for Intelligent Processing of Encrypted Data in Cloud-Based IoT

abstract: Phrase search allows retrieval of documents containing an exact phrase, which plays an important role in many machine learning applications for cloud-based IoT, such as intelligent medical data analytics. In order to protect sensitive information from being leaked by service providers, documents (e.g., clinic records) are usually encrypted by data owners before being outsourced to the cloud. This, however, makes the search operation an extremely challenging task. Existing searchable encryption schemes for multi-keyword search operations fail to perform phrase search, as they are unable to determine the location relationship of multiple keywords in a queried phrase over encrypted data on the cloud server side. In this paper, we propose P3, an efficient privacy-preserving phrase search scheme for intelligent encrypted data processing in cloud-based IoT. Our scheme exploits the homomorphic encryption and bilinear map to determine the location relationship of multiple queried keywords over encrypted data. It also utilizes a probabilistic trapdoor generation algorithm to protect users search patterns. Thorough security analysis demonstrates the security guarantees achieved by P3. We implement a prototype and conduct extensive experiments on real-world datasets. The evaluation results show that compared with existing multikeyword search schemes, P3 can greatly improve the search accuracy with moderate overheads.


title: Diffusion of charm and beauty in the Glasma

abstract: Relativistic nuclear collisions offer a unique way to study strong interactions at very high energy. The collision process can be described within the gluon saturation framework as the interaction of two colored glasses, and because of this interaction strong longitudinal gluon fields, namely the Glasma, are produced immediately after the collision. Besides, heavy quarks are also produced in the very early stage and because of their large mass and small concentration, their motion does not affect the evolution of the Glasma, thus behaving as ideal probes of the Glasma itself. We study the evolution of the heavy quarks in the Glasma allegedly produced in high energy p-Pb collisions by solving consistently the equations of motion of the quarks in the evolving Glasma fields. We find that this motion can be understood in terms of diffusion in momentum space, similarly to the random motion of a heavy probe in a hot thermalized medium. We show how the diffusion of heavy probes affects the nuclear modification factor of D and B mesons in p-Pb collisions.


title: Decentralized Optimal Merging Control for Connected and Automated Vehicles

abstract: This paper addresses the optimal control of Connected and Automated Vehicles (CAVs) arriving from two roads at a merging point where the objective is to jointly minimize the travel time and energy consumption of each CAV. The solution guarantees that a speed-dependent safety constraint is always satisfied, both at the merging point and everywhere within a control zone which precedes it. We first analyze the case of no active constraints and prove that under certain conditions the safety constraint remains inactive, thus significantly simplifying the determination of an explicit decentralized solution. When these conditions do not apply, an explicit solution is still obtained that includes intervals over which the safety constraint is active. Our analysis allows us to study the tradeoff between the two objective function components (travel time and energy within the control zone). Simulation examples are included to compare the performance of the optimal controller to a baseline with human-driven vehicles with results showing improvements in both metrics.


title: Adaptive O-CNN: A Patch-based Deep Representation of 3D Shapes

abstract: We present an Adaptive Octree-based Convolutional Neural Network (Adaptive O-CNN) for efficient 3D shape encoding and decoding. Different from volumetric-based or octree-based CNN methods that represent a 3D shape with voxels in the same resolution, our method represents a 3D shape adaptively with octants at different levels and models the 3D shape within each octant with a planar patch. Based on this adaptive patch-based representation, we propose an Adaptive O-CNN encoder and decoder for encoding and decoding 3D shapes. The Adaptive O-CNN encoder takes the planar patch normal and displacement as input and performs 3D convolutions only at the octants at each level, while the Adaptive O-CNN decoder infers the shape occupancy and subdivision status of octants at each level and estimates the best plane normal and displacement for each leaf octant. As a general framework for 3D shape analysis and generation, the Adaptive O-CNN not only reduces the memory and computational cost, but also offers better shape generation capability than the existing 3D-CNN approaches. We validate Adaptive O-CNN in terms of efficiency and effectiveness on different shape analysis and generation tasks, including shape classification, 3D autoencoding, shape prediction from a single image, and shape completion for noisy and incomplete point clouds.


title: On the Markus conjecture in convex case

abstract: In this paper, we show that any convex affine domain with a nonempty limit sets on the boundary under the action of the identity component of the automorphism group cannot cover a compact affine manifold with a parallel volume, which is a positive answer to the Markus conjecture for convex case. Consequently, we show that the Markus conjecture is true for convex affine manifolds of dimension $\leq 5$.


title: Some interior regularity estimates for solutions of complex Monge-Amp\`ere equations on a ball

abstract: In this paper, we consider the Dirichlet problem of a complex Monge-Amp\`ere equation on a ball in $\mathbb C^n$. With $\mathcal C^{1,\alpha}$ (resp. $\mathcal C^{0,\alpha}$) data, we prove an interior $\mathcal C^{1,\alpha}$ (resp. $\mathcal C^{0,\alpha}$) estimate for the solution. These estimates are generalized versions of the Bedford-Taylor interior $\mathcal C^{1,1}$ estimate.


title: The distribution of Weierstrass points on a tropical curve

abstract: We show that on a metric graph of genus $g$, a divisor of degree $n$ generically has $g(n-g+1)$ Weierstrass points. For a sequence of generic divisors on a metric graph whose degrees grow to infinity, we show that the associated Weierstrass points become distributed according to the Zhang canonical measure. In other words, the limiting distribution is determined by effective resistances on the metric graph. This distribution result has an analogue for complex algebraic curves, due to Neeman, and for curves over non-Archimedean fields, due to Amini.


title: Adversarial 3D Human Pose Estimation via Multimodal Depth Supervision

abstract: In this paper, a novel deep-learning based framework is proposed to infer 3D human poses from a single image. Specifically, a two-phase approach is developed. We firstly utilize a generator with two branches for the extraction of explicit and implicit depth information respectively. During the training process, an adversarial scheme is also employed to further improve the performance. The implicit and explicit depth information with the estimated 2D joints generated by a widely used estimator, in the second step, are together fed into a deep 3D pose regressor for the final pose generation. Our method achieves MPJPE of 58.68mm on the ECCV2018 3D Human Pose Estimation Challenge.


title: The baryon loading effect on relativistic astrophysical jet transport in the interstellar medium

abstract: The composition of the astrophysical relativistic jets remains uncertain. By kinetic particle-in-cell simulations, we show that the baryon component in the jet, or the so-called baryon loading effect (BLE), heavily affects relativistic jets transport dynamics in the interstellar medium. On the one hand, with the BLE, relativistic jets can transport in a much longer distance, because jet electrons draw a significant amount of energy from jet baryons via the Buneman-induced electrostatic waves and the Weibel-mediated collisionless shock; on the other hand, the jet electron phase space distribution may transform from a bottom-wide-single-peak structure to a center-wide-multiple-peak one by increasing the BLE, which largely influences the observed jet morphology. Implications for related astrophysical studies are also discussed.


title: A semi-model structure for Grothendieck weak 3-groupoids

abstract: In this paper we apply some tools developed in our previous work on Grothendieck $\infty$-groupoids to the finite-dimensional case of weak 3-groupoids. We obtain a semi-model structure on the category of Grothendieck 3-groupoids of suitable type, thanks to the construction of an endofunctor $\mathbb{P}$ that has enough structure to behave like a path object. This makes use of a recognition principle we prove here that characterizes globular theories whose models can be viewed as Grothendieck $n$-groupoids (for $0\leq n \leq \infty$). Finally, we prove that the obstruction in arbitrary dimension (possibly infinite) only resides in the construction of (slightly less than) a path object on a suitable category of Grothendieck (weak) $n$-categories with weak inverses. This also gives a sufficient condition for endowing an $n$-groupoid \`a la Batanin with the structure of a Grothendieck $n$-groupoid.


title: Sensitivity of future lepton colliders to the search for charged lepton flavor violation

abstract: The observation of lepton flavor violation indicates new physics beyond the Standard Model. Lepton colliders are ideal facilities to probe charged lepton flavor violation (CLFV) signals induced by new physics at high energy. In this work we perform a comprehensive study of the sensitivity of future lepton colliders to charged lepton flavor violation. We consider the most general renormalizable Lagrangian coupling two leptons to new bosonic particles, involving both $\Delta L=0$ and $\Delta L=2$ interactions. The CLFV processes are introduced by the exchange of off-shell new particles at tree-level. We find that CEPC, ILC, FCC-ee, and CLIC each provide a complementary probe of CLFV couplings to low-energy precision experiments for $\tau$ lepton(s) in final states, while low-energy precision experiments are more sensitive in the absence of $\tau$ leptons.


title: Transport of relativistic jet in the magnetized intergalactic medium

abstract: The kinetic effects of magnetic fields on the transport of relativistic jet in the intergalactic medium remain uncertain, especially for their perpendicular component. By particle-in-cell simulations, we find that when only jet electrons are fully magnetized, they are directly deflected by the magnetic field, but jet protons are mainly dragged by collective charge-separation electric field. However, when both electrons and protons are fully magnetized, the contrary is the case. Their balance tremendously distorts the jet density and electromagnetic fields, leading to enormous energy exchange between different species and fields. As a result, the electron spectrum energy distribution (SED) gets reshaped and the power law slope of the SED decreases as the magnetic field strength increases. In other words, we may infer that magnetic fields around the relativistic jet play a crucial role in shaping the observed SED.


title: Experimental realization of an intrinsic magnetic topological insulator

abstract: Intrinsic magnetic topological insulator (TI) is a stoichiometric magnetic compound possessing both inherent magnetic order and topological electronic states. Such a material can provide a shortcut to various novel topological quantum effects but remains elusive experimentally so far. Here, we report the experimental realization of high-quality thin films of an intrinsic magnetic TI---MnBi$_2$Te$_4$---by alternate growth of a Bi$_2$Te$_3$ quintuple-layer and a MnTe bilayer with molecular beam epitaxy. The material shows the archetypical Dirac surface states in angle-resolved photoemission spectroscopy and is demonstrated to be an antiferromagnetic topological insulator with ferromagnetic surfaces by magnetic and transport measurements as well as first-principles calculations. The unique magnetic and topological electronic structures and their interplays enable the material to embody rich quantum phases such as quantum anomalous Hall insulators and axion insulators in a well-controlled way.


title: S3BD: Secure Semantic Search over Encrypted Big Data in the Cloud

abstract: Cloud storage is a widely utilized service for both personal and enterprise demands. However, despite its advantages, many potential users with enormous amounts of sensitive data (big data) refrain from fully utilizing the cloud storage service due to valid concerns about data privacy. An established solution to the cloud data privacy problem is to perform encryption on the client-end. This approach, however, restricts data processing capabilities (eg, searching over the data). Accordingly, the research problem we investigate is how to enable real-time searching over the encrypted big data in the cloud. In particular, semantic search is of interest to clients dealing with big data. To address this problem, in this research, we develop a system (termed S3BD) for searching big data using cloud services without exposing any data to cloud providers. To keep real-time response on big data, S3BD proactively prunes the search space to a subset of the whole dataset. For that purpose, we propose a method to cluster the encrypted data. An abstract of each cluster is maintained on the client-end to navigate the search operation to appropriate clusters at the search time. Results of experiments, carried out on real-world big datasets, demonstrate that the search operation can be achieved in real-time and is significantly more efficient than other counterparts. In addition, a fully functional prototype of S3BD is made publicly available.


title: A Prospect Theoretic Approach for Trust Management in IoT Networks under Manipulation Attacks

abstract: As Internet of Things (IoT) and Cyber-Physical systems become more ubiquitous and an integral part of our daily lives, it is important that we are able to trust the data aggregate from such systems. However, the interpretation of trustworthiness is contextual and varies according to the risk tolerance attitude of the concerned application and varying levels of uncertainty associated with the evidence upon which trust models act. Hence, the data integrity scoring mechanisms should have provisions to adapt to varying risk attitudes and uncertainties. In this paper, we propose a Bayesian inference model and a prospect theoretic framework for data integrity scoring that quantify the trustworthiness of data collected from IoT devices in the presence of an adversaries who manipulate the data. We consider an imperfect anomaly monitoring mechanism that monitors the data being sent from each device and classifies the outcome as not compromised, compromised, and cannot be inferred. These outcomes are conceptualized as a multinomial hypothesis of a Bayesian inference model with three parameters which are then used for calculating a utility value on how reliable the aggregate data is. We use a prospect theory inspired approach to quantify this data integrity score and evaluate the trustworthiness of the aggregate data from the IoT framework. Furthermore, we also model the system using the traditionally used expected utility theory and compare the results with that obtained using prospect theory. As decisions are based on how the data is fused, we propose two measuring models: one optimistic and another conservative. The proposed framework is validated using extensive simulation experiments. We show how data integrity scores vary under a variety of system factors like attack intensity and inaccurate detection.


title: A Time Domain Acoustic Model for the Production of Rodent Ultrasonic Vocalizations

abstract: Mammalian ultrasonic vocalization (USV) has been a subject of interest for decades. This interest has mainly been driven by the intelligence of dolphins and other odontocetes. However the semantic content of odontocete USV and its mechanism of production remain poorly understood. Serendipitously however many rodent species have convergently evolved the ability to produce USVs in a similar manner. In this paper we use rodent USV as a model process to help us gain insight into the production mechanism for mammalian USV as a whole. We derive a model that describes the production of rodent USVs by considering the interaction of an unstable jet, emerging from the vocal folds, with the passive resonance modes of the upper vocal tract. Thus our model is also a solution to a special case of the jet susceptibility problem. The derived model takes the form of a set of coupled nonlinear time domain ODEs, whose solutions are controlled by biologically relevant parameters such as subglottal pressure and vocal fold radius. In our analysis of the model we find the existence of a subglottal blowing pressure threshold ($p \approx 710$ Pa), above which steady acoustic oscillations occur. Furthermore we also reproduce the $22$ kHz rat alarm call at realistic blowing pressures ($p \approx 1500$ Pa)


title: Infrared spectra of Rg1,2 - C6H6 complexes, Rg = He, Ne, Ar

abstract: Infrared spectra of Rg1,2 - C6H6 complexes (Rg = He, Ne, Ar) are observed in the region of the nu12 fundamental of C6H6 using a pulsed supersonic jet expansion and a tunable optical parametric oscillator laser source. The mixed trimer He - Ne - C6H6 is also detected. Four bands are analyzed for each complex, namely nu12 itself (~3048 cm-1) and three linked combination bands (~3079, 3100, and 3102 cm-1). The results are consistent with previous ultraviolet and microwave results, with Ne2 - C6H6 and He - Ne - C6H6 being analyzed spectroscopically here for the first time.


title: A Geometric Observer for Scene Reconstruction Using Plenoptic Cameras

abstract: This paper proposes an observer for generating depth maps of a scene from a sequence of measurements acquired by a two-plane light-field (plenoptic) camera. The observer is based on a gradient-descent methodology. The use of motion allows for estimation of depth maps where the scene contains insufficient texture for static estimation methods to work. A rigourous analysis of stability of the observer error is provided, and the observer is tested in simulation, demonstrating convergence behaviour.


title: Double Grothendieck Polynomials for Symplectic and Odd Orthogonal Grassmannians

abstract: We study the double Grothendieck polynomials of Kirillov--Naruse for the symplectic and odd orthogonal Grassmannians. These functions are explicitly written as sums of Pfaffian and are identified with the stable limits of the fundamental classes of Schubert varieties in the torus equivariant connective K-theory of these isotropic Grassmannians. We also provide a combinatorial description of the ring formally spanned by double Grothendieck polynomials.


title: Semi De Morgan logic properly displayed

abstract: In the present paper, we endow a family of axiomatic extensions of semi De Morgan logic with proper multi-type display calculi which are sound, complete, conservative, and enjoy cut elimination and subformula property. Our proposal builds on an algebraic analysis of semi De Morgan algebras and its subvarieties and applies the guidelines of the multi-type methodology in the design of display calculi.


title: Residual dimension of nilpotent groups

abstract: The functions $F_{G}(n)$ measures the asymptotic behavior of residual finiteness for a finitely generated group $G$. In previous work \cite{Pengitore_1}, the author claimed a characterization for $F_{N}(n)$ when $N$ is a finitely generated nilpotent group. However, a counterexample to the above claim was communicated to the author, and subsequently, the statement of the asymptotic characterization of $F_{N}(n)$ is incorrect. In this article, we introduce new tools to provide lower asymptotic bounds for $F_{N}(n)$ when $N$ is a finitely generated nilpotent group. Moreover, we introduce a class of finitely generated nilpotent groups for which the upper bound of \cite{Pengitore_1} can be improved. Finally, we construct of a class of finitely generated nilpotent groups $N$ for which the asymptotic behavior of $F_N(n)$ can be fully characterized.


title: A Deep Near-Infrared [Fe II]+[Si I] Emission Line Image of the Supernova Remnant Cassiopeia A

abstract: We present a long-exposure (~10 hr) image of the supernova (SN) remnant Cassiopeia A (Cas A) obtained with the UKIRT 3.8-m telescope using a narrow band filter centered at 1.644 um emission. The passband contains [Fe II] 1.644 um and [Si I] 1.645 um lines, and our `deep [Fe II]+[Si I] image' provides an unprecedented panoramic view of Cas A, showing both shocked and unshocked SN ejecta together with shocked circumstellar medium at subarcsec (~0.7 arcsec or 0.012 pc) resolution. The diffuse emission from the unshocked SN ejecta has a form of clumps, filaments, and arcs, and their spatial distribution correlates well with that of the Spitzer [Si II] infrared emission, suggesting that the emission is likely due to [Si I] line not [Fe II] line as in shocked material. The structure of the optically-invisible western area of Cas A is clearly seen for the first time. The area is filled with many Quasi-Stationary Flocculi (QSFs) and fragments of the disrupted ejecta shell. We suggest that the anomalous radio properties in this area could be due to the increased number of such dense clumps. We identified 309 knots in the deep [Fe II]+[Si I] image and classified them into QSFs and fast-moving knots (FMKs). The total H+He mass of QSFs is ~0.23 Msun, implying that the mass fraction of dense clumps in the progenitor's red-supergiant wind is 4--13%. The spatial distribution of QSFs suggests that there had been a highly asymmetric mass loss $10^4$--$10^5$ yr before the SN explosion. The mass of the [Fe II] line-emitting, shocked dense Fe ejecta is ~3x$10^{-5}$ Msun. The comparison with the ionic S-line dominated Hubble Space Telescope WFC3/IR image suggests that the outermost FMKs in the southeastern area are Fe-rich.


title: Computational modelling of cardiac ischaemia using a variable-order fractional Laplacian

abstract: Heart failure is one of the most common causes of death in the western world. Many heart problems are linked to disturbances in cardiac electrical activity, such as wave re-entry caused by ischaemia. In terms of mathematical modelling, the monodomain equation is widely used to model electrical activity in the heart. Recently, Bueno-Orovio et al. [J. R. Soc. Interface 11: 20140352, 2014] pioneered the use of a fractional Laplacian operator in the monodomain equation to account for the complex heterogeneous structures in heart tissue. In this work we consider how to extend this approach to apply to hearts with regions of damaged tissue. This requires the use of a fractional Laplacian operator whose fractional order varies spatially. We develop efficient numerical methods capable of solving this challenging problem on domains ranging from simple one-dimensional intervals with uniform meshes, through to full three-dimensional geometries on unstructured meshes. Results are presented for several test problems in one dimension, demonstrating the effects of different fractional orders in regions of healthy and damaged tissue. Then we showcase some new results for a three-dimensional fractional monodomain equation with a Beeler-Reuter ionic current model on a rabbit heart mesh. These simulation results are found to exhibit wave re-entry behaviour, brought about only by varying the value of the fractional order in a region representing damaged tissue.


title: Bugs in Infrastructure as Code

abstract: Infrastructure as code (IaC) scripts are used to automate the maintenance and configuration of software development and deployment infrastructure. IaC scripts can be complex in nature, containing hundreds of lines of code, leading to defects that can be difficult to debug, and lead to wide-scale system discrepancies such as service outages at scale. Use of IaC scripts is getting increasingly popular, yet the nature of defects that occur in these scripts have not been systematically categorized. A systematic categorization of defects can inform practitioners about process improvement opportunities to mitigate defects in IaC scripts. The goal of this paper is to help software practitioners improve their development process of infrastructure as code (IaC) scripts by categorizing the defect categories in IaC scripts based upon a qualitative analysis of commit messages and issue report descriptions. We mine open source version control systems collected from four organizations namely, Mirantis, Mozilla, Openstack, and Wikimedia Commons to conduct our research study. We use 1021, 3074, 7808, and 972 commits that map to 165, 580, 1383, and 296 IaC scripts, respectively, collected from Mirantis, Mozilla, Openstack, and Wikimedia Commons. With 89 raters we apply the defect type attribute of the orthogonal defect classification (ODC) methodology to categorize the defects. We also review prior literature that have used ODC to categorize defects, and compare the defect category distribution of IaC scripts with 26 non-IaC software systems. Respectively, for Mirantis, Mozilla, Openstack, and Wikimedia Commons, we observe (i) 49.3%, 36.5%, 57.6%, and 62.7% of the IaC defects to contain syntax and configuration-related defects; (ii) syntax and configuration-related defects are more prevalent amongst IaC scripts compared to that of previously-studied non-IaC software.


title: Asteroseismic detection of latitudinal differential rotation in 13 Sun-like stars

abstract: The differentially rotating outer layers of stars are thought to play a role in driving their magnetic activity, but the underlying mechanisms that generate and sustain differential rotation are poorly understood. We report the measurement of latitudinal differential rotation in the convection zones of 40 Sun-like stars using asteroseismology. For the most significant detections, the stars' equators rotate approximately twice as fast as their mid-latitudes. The latitudinal shear inferred from asteroseismology is much larger than predictions from numerical simulations.


title: Matrix product solution of the stationary states of two-species open zero range processes

abstract: Using the matrix product ansatz, we obtain solutions of the steady-state distribution of the two-species open one-dimensional zero range process. Our solution is based on a conventionally employed constraint on the hop rates, which eventually allows us to simplify the constituent matrices of the ansatz. It is shown that the matrix at each site is given by the tensor product of two sets of matrices and the steady-state distribution assumes an inhomogeneous factorized form. Our method can be generalized to the cases of more than two species of particles.


title: Printing-while-moving: a new paradigm for large-scale robotic 3D Printing

abstract: Building and Construction have recently become an exciting application ground for robotics. In particular, rapid progress in materials formulation and in robotics technology has made robotic 3D Printing of concrete a promising technique for in-situ construction. Yet, scalability remains an important hurdle to widespread adoption: the printing systems (gantry- based or arm-based) are often much larger than the structure to be printed, hence cumbersome. Recently, a mobile printing system - a manipulator mounted on a mobile base - was proposed to alleviate this issue: such a system, by moving its base, can potentially print a structure larger than itself. However, the proposed system could only print while being stationary, imposing thereby a limit on the size of structures that can be printed in a single take. Here, we develop a system that implements the printing-while-moving paradigm, which enables printing single-piece structures of arbitrary sizes with a single robot. This development requires solving motion planning, localization, and motion control problems that are specific to mobile 3D Printing. We report our framework to address those problems, and demonstrate, for the first time, a printing-while-moving experiment, wherein a 210 cm x 45 cm x 10 cm concrete structure is printed by a robot arm that has a reach of 87 cm.


title: LIDAR-Camera Fusion for Road Detection Using Fully Convolutional Neural Networks

abstract: In this work, a deep learning approach has been developed to carry out road detection by fusing LIDAR point clouds and camera images. An unstructured and sparse point cloud is first projected onto the camera image plane and then upsampled to obtain a set of dense 2D images encoding spatial information. Several fully convolutional neural networks (FCNs) are then trained to carry out road detection, either by using data from a single sensor, or by using three fusion strategies: early, late, and the newly proposed cross fusion. Whereas in the former two fusion approaches, the integration of multimodal information is carried out at a predefined depth level, the cross fusion FCN is designed to directly learn from data where to integrate information; this is accomplished by using trainable cross connections between the LIDAR and the camera processing branches. To further highlight the benefits of using a multimodal system for road detection, a data set consisting of visually challenging scenes was extracted from driving sequences of the KITTI raw data set. It was then demonstrated that, as expected, a purely camera-based FCN severely underperforms on this data set. A multimodal system, on the other hand, is still able to provide high accuracy. Finally, the proposed cross fusion FCN was evaluated on the KITTI road benchmark where it achieved excellent performance, with a MaxF score of 96.03%, ranking it among the top-performing approaches.