title: Control Flow Graph Modifications for Improved RF-Based Processor Tracking Performance

abstract: Many dedicated embedded processors do not have memory or computational resources to coexist with traditional (host-based) security solutions. As a result, there is interest in using out-of-band analog side-channel measurements and their analyses to accurately monitor and analyze expected program execution. In this paper, we describe an approach to this problem using externally observable multi-band radio frequency (RF) measurements to make inferences about a program's execution. Because it is very difficult to identify individual instructions solely from their RF emissions, we compare RF measurements with the constrained execution logic of the program so that multiple RF measurements over time can effectively track program execution dynamically. In our approach, a program's execution is modeled by control flow graphs (CFG) and transitions between nodes of such graphs. We demonstrate that tracking performance can be improved through applications program modifications such as changing basic block transition properties and/or adding new basic blocks that are highly observable. In addition to demonstrating these principled approaches on some simple programs, we present initial results on the complexity and structure of real-world applications programs, namely gzip and md5sum, in this modeling framework.


title: Non-Orthogonal Multiple Access: Common Myths and Critical Questions

abstract: Non-orthogonal multiple access (NOMA) has received tremendous attention for the design of radio access techniques for fifth generation (5G) wireless networks and beyond. The basic concept behind NOMA is to serve more than one user in the same resource block, e.g., a time slot, subcarrier, spreading code, or space. With this, NOMA promotes massive connectivity, lowers latency, improves user fairness and spectral efficiency, and increases reliability compared to orthogonal multiple access (OMA) techniques. While NOMA has gained significant attention from the communications community, it has also been subject to several widespread misunderstandings, such as $``\textit{NOMA is based on allocating higher power to users with worse channel conditions. As such, cell-edge users receive more power in NOMA and due to this biased power allocation toward cell-edge users inter-cell interference is more severe in NOMA compared to OMA. NOMA also compromises security for spectral efficiency.}''$ The above statements are actually false, and this paper aims at identifying such common myths about NOMA and clarifying why they are not true. We also pose critical questions that are important for the effective adoption of NOMA in 5G and beyond and identify promising research directions for NOMA, which will require intense investigation in the future.


title: Deterministic limit of temporal difference reinforcement learning for stochastic games

abstract: Reinforcement learning in multiagent systems has been studied in the fields of economic game theory, artificial intelligence and statistical physics by developing an analytical understanding of the learning dynamics (often in relation to the replicator dynamics of evolutionary game theory). However, the majority of these analytical studies focuses on repeated normal form games, which only have a single environmental state. Environmental dynamics, i.e., changes in the state of an environment affecting the agents' payoffs has received less attention, lacking a universal method to obtain deterministic equations from established multistate reinforcement learning algorithms. In this work we present a novel methodological extension, separating the interaction from the adaptation time scale, to derive the deterministic limit of a general class of reinforcement learning algorithms, called temporal difference learning. This form of learning is equipped to function in more realistic multistate environments by using the estimated value of future environmental states to adapt the agent's behavior. We demonstrate the potential of our method with the three well established learning algorithms Q learning, SARSA learning and Actor-Critic learning. Illustrations of their dynamics on two multiagent, multistate environments reveal a wide range of different dynamical regimes, such as convergence to fixed points, limit cycles, and even deterministic chaos.


title: Critical parameters for reaction-diffusion equations involving space-time fractional derivatives

abstract: We will look at reaction-diffusion type equations of the following type, $$\partial^\beta_tV(t,x)=-(-\Delta)^{\alpha/2} V(t,x)+I^{1-\beta}_t[V(t,x)^{1+\eta}].$$ We first study the equation on the whole space by making sense of it via an integral equation. Roughly speaking, we will show that when $0<\eta\leq\eta_c$, there is no global solution other than the trivial one while for $\eta>\eta_c$, non-trivial global solutions do exist. We also study the equation on a bounded domain with Dirichlet boundary condition and show that the presence of the time derivative induces a significant change in the behaviour of the solution.


title: A Template-based $\gamma$-ray Reconstruction Method for Air Shower Arrays

abstract: We introduce a new Monte Carlo template-based reconstruction method for air shower arrays, with a focus on shower core and energy reconstruction of $\gamma$-ray induced air showers. The algorithm fits an observed lateral amplitude distribution of an extensive air shower against an expected probability distribution using a likelihood approach. A full Monte Carlo air shower simulation in combination with the detector simulation is used to generate the expected probability distributions. The goodness of fit can be used to discriminate between $\gamma$-ray and hadron induced air showers. As an example, we apply this method to the High Altitude Water Cherenkov $\gamma$-ray Observatory and its recently installed high-energy upgrade. The performance of this method and the applicability to air shower arrays with mixed detector types makes it a promising reconstruction approach for current and future instruments.


title: Out-of-Time-Order Correlation Functions for Unitary Minimal Models

abstract: We analytically study the Out-of-Time-Order Correlation functions (OTOC) for two spatially separated primary operators in two-dimensional unitary minimal models. Besides giving general arguments using the conformal symmetry, we also use the Coulomb gas formalism to explicitly calculate the OTOC across the full time regime. In contrast to large-$N$ chaotic systems, these models do not display a separation of time scales, due to the lack of a large parameter. We find that the physics at early times ($0<t-x\ll\beta$) and late times ($t-x\gg\beta$) are controlled by different OPE channels, which are related to each other via the braiding matrix. The normalized OTOC obeys a power-law decay with a fractional power at early times and approaches a generically nonzero value in an exponential way at late times. The late time value is related to the modular $S$-matrix and is in agreement with earlier calculations. All of the results above are readily generalized to rational conformal field theories.


title: Spin signatures in the electrical response of graphene nanogaps

abstract: We analyse the electrical response of narrow graphene nanogaps in search for transport signatures stemming from spin-polarized edge states. We find that the electrical transport across graphene nanogaps having perfectly defined zigzag edges does not carry any spin-related signature. We also analyse the magnetic and electrical properties of nanogaps whose electrodes have wedges that possibly occur in the currently fabricated nanogaps. These wedges can host spin polarized wedge low-energy states due to the bipartite nature of the graphene lattice. We find that these spin-polarized low-energy modes give rise to low-voltage signatures in the differential conductance and to distinctive features in the stability diagrams. These are caused by fully spin-polarized currents.


title: Sensitivity Function Trade-offs for Networks with a String Topology

abstract: We present two sensitivity function trade-offs that apply to a class of networks with a string topology. In particular we show that a lower bound on the H-infinity norm and a Bode sensitivity relation hold for an entire family of sensitivity functions associated with growing the network. The trade-offs we identify are a direct consequence of growing the network, and can be used to explain why poorly regulated low frequency behaviours emerge in long vehicle platoons even when using dynamic feedback.


title: Overall constitutive description of symmetric layered media based on scattering of oblique SH waves

abstract: This papers investigates the scattering of oblique shear horizontal (SH) waves off finite periodic media made of elastic and viscoelastic layers. It further considers whether a Willis-type constitutive matrix (in temporal and spatial Fourier domain) may reproduce the scattering matrix (SM) of such a system. In answering this question the procedure to determine the relevant overall constitutive parameters for such a medium is presented. To do this, first the general form of the dispersion relation and impedances for oblique SH propagation in such coupled Willis-type media are developed. The band structure and scattering of layered media are calculated using the transfer matrix (TM) method. The dispersion relation may be derived based on the eigensolutions of an infinite periodic domain. The wave impedances associated with the exterior surfaces of a finite thickness slab are extracted from the scattering of such a system. Based on reciprocity and available symmetries of the structure and each constituent layer, the general form of the dispersion and impedances may be simplified. The overall quantities may be extracted by equating the scattering data from TM with those expected from a Willis-type medium. It becomes evident that a Willis-type coupled constitutive tensor with components that are assumed independent of wave vector is unable to reproduce all oblique scattering data. Therefore, non-unique wave vector dependent formulations are introduced, whose SM matches that of the layered media exactly. It is further shown that the dependence of the overall constitutive tensors of such systems on the wave vector is not removable even at very small frequencies and incidence angles and that analytical considerations significantly limit the potential forms of the spatially dispersive constitutive tensors.


title: Modelling the data and not the images in FMRI

abstract: The standard approach to the analysis of functional magnetic resonance imaging (FMRI) data applies various preprocessing steps to the original FMRI. These preprocessings lead to a general underestimation of residual variance in the downstream analysis. This negatively impacts the type I error of statistical tests and increases the risk for reporting false positive results. A genuine approach to the statistical analysis of FMRI data of brain scans is derived from first principles that is deeply rooted in statistical test theory. The method combines all preprocessing steps of the standard approach into one single modelling step, enabling valid statistical tests to be constructed. On population level, BOLD effects are modelled by random effects meta regression models. This acknowledges that subjects are random entities, and it acknowledges that the accuracy of the BOLD signal is estimated with various certainty in an FMRI. The concept of a reference scalar field is introduced that enables individual effect sizes to be related to each other with respect to a common unit. In particular, multicentre studies will gain interpretability and power by its use.


title: Local moduli of scalar-flat K\"ahler ALE surfaces

abstract: In this article, we give a survey of our construction of a local moduli space of scalar-flat K\"ahler ALE metrics in complex dimension $2$. We also prove an explicit formula for the dimension of this moduli space on a scalar-flat K\"ahler ALE surface which deforms to the minimal resolution of $\mathbb{C}^2/\Gamma$, where $\Gamma$ is a finite subgroup of ${\rm{U}}(2)$ without complex reflections, in terms of the embedding dimension of the singularity.


title: Unsupervised cross-lingual matching of product classifications

abstract: Unsupervised cross-lingual embeddings mapping has provided a unique tool for completely unsupervised translation even for languages with different scripts. In this work we use this method for the task of unsupervised cross-lingual matching of product classifications. Our work also investigates limitations of unsupervised vector alignment and we also suggest two other techniques for aligning product classifications based on their descriptions: using hierarchical information and translations.


title: Approximating microswimmer dynamics by active Brownian motion: Energetics and efficiency

abstract: We consider the dynamics of a microswimmer and show that they can be approximated by active Brownian motion. The swimmer is modeled by coupled overdamped Langevin equations with periodic driving. We compare the energy dissipation of the real swimmer to that of the active Brownian motion model, finding that the latter can massively underestimate the complete dissipation. This discrepancy is related to the inability to infer the full dissipation from partial observation of the complete system. We introduce an efficiency that measures how much of the dissipated energy is spent on forward propulsion.


title: Counterexample to Equivalent Nodal Analysis for Voltage Stability Assessment

abstract: Existing literature claims that the L-index for voltage instability detection is inaccurate and proposes an improved index quantifying voltage stability through system equivalencing. The proposed stability condition is claimed to be exact in determining voltage instability.We show the condition is incorrect through simple arguments accompanied by demonstration on a two-bus system counterexample.


title: A global weak solution to the Lorentzian harmonic map flow

abstract: We investigate a parabolic-elliptic system which is related to a harmonic map from a compact Riemann surface with a smooth boundary into a Lorentzian manifold with a warped product metric. We prove that there exists a unique global weak solution for this system which is regular except for at most finitely many singular points.


title: Pose Estimation for Non-Cooperative Spacecraft Rendezvous Using Convolutional Neural Networks

abstract: On-board estimation of the pose of an uncooperative target spacecraft is an essential task for future on-orbit servicing and close-proximity formation flying missions. However, two issues hinder reliable on-board monocular vision based pose estimation: robustness to illumination conditions due to a lack of reliable visual features and scarcity of image datasets required for training and benchmarking. To address these two issues, this work details the design and validation of a monocular vision based pose determination architecture for spaceborne applications. The primary contribution to the state-of-the-art of this work is the introduction of a novel pose determination method based on Convolutional Neural Networks (CNN) to provide an initial guess of the pose in real-time on-board. The method involves discretizing the pose space and training the CNN with images corresponding to the resulting pose labels. Since reliable training of the CNN requires massive image datasets and computational resources, the parameters of the CNN must be determined prior to the mission with synthetic imagery. Moreover, reliable training of the CNN requires datasets that appropriately account for noise, color, and illumination characteristics expected in orbit. Therefore, the secondary contribution of this work is the introduction of an image synthesis pipeline, which is tailored to generate high fidelity images of any spacecraft 3D model. The proposed technique is scalable to spacecraft of different structural and physical properties as well as robust to the dynamic illumination conditions of space. Through metrics measuring classification and pose accuracy, it is shown that the presented architecture has desirable robustness and scalable properties.


title: Exact combinatorial approach to finite coagulating systems through recursive equations

abstract: This work outlines an exact combinatorial approach to finite coagulating systems through recursive equations and use of generating function method. In the classic approach the mean-field Smoluchowski coagulation is used. However, the assumptions of the mean-field theory are rarely met in real systems which limits the accuracy of the solution. In our approach, cluster sizes and time are discrete, and the binary aggregation alone governs the time evolution of the systems. By considering the growth histories of all possible clusters and applying monodisperse initial conditions, the exact expression for the probability of finding a coagulating system with an arbitrary kernel in a given cluster configuration is derived. Then, the average number of such clusters and the standard deviation of these solutions can be calculated. In this work, recursive equations for all possible growth histories of clusters are introduced. The correctness of our expressions was proved based on the comparison with numerical results obtained for systems with constant, multiplicative and additive kernels. For the first time the exact solutions for the multiplicative and additive kernels were obtained with this framework. In addition, our results were compared with the results arising from the solutions to the mean-field Smoluchowski equation. Our theoretical predictions outperform the classic approach.


title: Graph magnitude homology via algebraic Morse theory

abstract: We compute magnitude homology of various graphs using algebraic Morse theory. Specifically, we (1) give an alternative proof that trees are diagonal, (2) identify a new class of diagonal graphs, (3) prove that the icosahedral graph is diagonal, and (4) compute the magnitude homology of cycles. These results answer several questions of Hepworth and Willerton [HW17].


title: Energy identity for a class of approximate Dirac-harmonic maps from surfaces with boundary

abstract: For a sequence of coupled fields $\{(\phi_n,\psi_n)\}$ from a compact Riemann surface $M$ with smooth boundary to a general compact Riemannian manifold with uniformly bounded energy and satisfying the Dirac-harmonic system up to some uniformly controlled error terms, we show that the energy identity holds during a blow-up process near the boundary. As an application to the heat flow of Dirac-harmonic maps from surfaces with boundary, when such a flow blows up at infinite time, we obtain an energy identity.


title: TESS Discovery of an ultra-short-period planet around the nearby M dwarf LHS 3844

abstract: Data from the newly-commissioned \textit{Transiting Exoplanet Survey Satellite} (TESS) has revealed a "hot Earth" around LHS 3844, an M dwarf located 15 pc away. The planet has a radius of $1.32\pm 0.02$ $R_\oplus$ and orbits the star every 11 hours. Although the existence of an atmosphere around such a strongly irradiated planet is questionable, the star is bright enough ($I=11.9$, $K=9.1$) for this possibility to be investigated with transit and occultation spectroscopy. The star's brightness and the planet's short period will also facilitate the measurement of the planet's mass through Doppler spectroscopy.


title: A threshold for cutoff in two-community random graphs

abstract: In this paper, we are interested in the impact of communities on the mixing behavior of the non-backtracking random walk. We consider sequences of sparse random graphs of size $N$ generated according to a variant of the classical configuration model which incorporates a two-community structure. The strength of the bottleneck is measured by a parameter $\alpha$ which roughly corresponds to the fraction of edges that go from one community to the other. We show that if $\alpha\gg \frac{1}{\log N}$, then the non-backtracking random walk exhibits cutoff at the same time as in the one-community case, but with a larger cutoff window, and that the distance profile inside this window converges to the Gaussian tail function. On the other hand, if $\alpha \ll \frac{1}{\log N}$ or $\alpha \asymp \frac{1}{\log N}$, then the mixing time is of order $1/\alpha$ and there is no cutoff.


title: Prime-Residue-Class of Uniform Charges on the Integers

abstract: There is a probability charge on the power set of the integers that gives probability $1/p$ to every residue class modulo a prime $p$. There exists such a charge that gives probability $w$ to the set of prime numbers iff $w \in [0,1/2]$. Similarly, there is such a charge that gives probability $x$ to a residue class modulo $c$, where $c$ is composite, iff $x \in [0,1/y]$, where $y$ is the largest prime factor of $c$.


title: Fuzzy Based Wellness Analyzer

abstract: Social health and emotional wellness is a matter of concern in today's urban world. Being the part of a metropolis has an effect on mental health through the influence of increased stressors and factors such as overcrowded and polluted environment, high levels of violence, and reduced social support. It is important to realize that only healthy citizens can constitute together a smart city. In this paper, we present a fuzzy-based approach for analyzing the well being of a person. We track the general day to day activities of a person and analyze its performance. To do so, we divide the factors affecting the wellness of a person into three components which are the physical, productive and social. Using these parameters, we output a coefficient for the overall well being of a person.


title: The qualitative behavior at the free boundary for approximate harmonic maps from surfaces

abstract: Let $\{u_n\}$ be a sequence of maps from a compact Riemann surface $M$ with smooth boundary to a general compact Riemannian manifold $N$ with free boundary on a smooth submanifold $K\subset N$ satisfying \[ \sup_n \ \left(\|\nabla u_n\|_{L^2(M)}+\|\tau(u_n)\|_{L^2(M)}\right)\leq \Lambda, \] where $\tau(u_n)$ is the tension field of the map $u_n$. We show that the energy identity and the no neck property hold during a blow-up process. The assumptions are such that this result also applies to the harmonic map heat flow with free boundary, to prove the energy identity at finite singular time as well as at infinity time. Also, the no neck property holds at infinity time.


title: Time to Stop Telling Biophysics Students that Light is Primarily a Wave

abstract: Standard pedagogy introduces optics as though it were a consequence of Maxwell's equations, and only grudgingly admits, usually in a rushed aside, that light has a particulate character that can somehow be reconciled with the wave picture. Recent revolutionary advances in optical imaging, however, make this approach more and more unhelpful: How are we to describe two-photon imaging, FRET, localization microscopy, and a host of related techniques to students who think of light primarily as a wave? I was surprised to find that everything I wanted my biophysics students to know about light, including image formation, x-ray diffraction, and even Bessel beams, could be expressed as well (or better) from the quantum viewpoint pioneered by Richard Feynman. Even my undergraduate students grasp this viewpoint as well as (or better than) the traditional one, and by mid-semester they are already well positioned to integrate the latest advances into their understanding. Moreover, I have found that this approach clarifies my own understanding of new techniques.


title: Algorithmic aspects of broadcast independence

abstract: An independent broadcast on a connected graph $G$ is a function $f:V(G)\to \mathbb{N}_0$ such that, for every vertex $x$ of $G$, the value $f(x)$ is at most the eccentricity of $x$ in $G$, and $f(x)>0$ implies that $f(y)=0$ for every vertex $y$ of $G$ within distance at most $f(x)$ from $x$. The broadcast independence number $\alpha_b(G)$ of $G$ is the largest weight $\sum\limits_{x\in V(G)}f(x)$ of an independent broadcast $f$ on $G$. We describe an efficient algorithm that determines the broadcast independence number of a given tree. Furthermore, we show NP-hardness of the broadcast independence number for planar graphs of maximum degree four, and hardness of approximation for general graphs. Our results solve problems posed by Dunbar, Erwin, Haynes, Hedetniemi, and Hedetniemi (2006), Hedetniemi (2006), and Ahmane, Bouchemakh, Sopena (2018).


title: The Measure Aspect of Quantum Uncertainty, of Entanglement, and the Associated Entropies

abstract: Indeterminacy associated with probing of a quantum state is commonly expressed through spectral distances (metric) featured in the outcomes of repeated experiments. Here we express it as an effective amount (measure) of distinct outcomes instead. The resulting $\mu$-uncertainties are described by the effective number theory [1] whose central result, the existence of a minimal amount, leads to a well-defined notion of intrinsic irremovable uncertainty. We derive $\mu$-uncertainty formulas for arbitrary set of commuting operators, including the cases with continuous spectra. The associated entropy-like characteristics, the $\mu$-entropies, convey how many degrees of freedom are effectively involved in a given measurement process. In order to construct quantum $\mu$-entropies, we are led to quantum effective numbers designed to count independent, mutually orthogonal states effectively comprising a density matrix. This concept is basis-independent and leads to a measure-based characterization of entanglement.

categories: quant-ph cond-mat.dis-nn hep-lat hep-th math-ph math.MP


title: Wild pseudohyperbolic attractor in a four-dimensional Lorenz system

abstract: We present an example of a new strange attractor which, as we show, belongs to a class of wild pseudohyperbolic spiral attractors. We find this attractor in a four-dimensional system of differential equations which can be represented as an extension of the Lorenz system.


title: Exploiting atomic layer deposition for fabricating sub-10 nm X-ray lenses

abstract: Moving towards significantly smaller nanostructures, direct structuring techniques such as electron beam lithography approach fundamental limitations in feature size and aspect ratios. Application of nanostructures like diffractive X-ray lenses requires feature sizes of below 10 nm to enter a new regime in high resolution X-raymicroscopy. As such dimensions are difficult to obtain using conventional electron beam lithography, we pursue a line-doubling approach. We demonstrate that thismethod yields structure sizes as small as 6.4 nm. X-ray lenses fabricated in this way are tested for their efficiency and microscopic resolution. In addition, the line-doubling technique is successfully extended to a six-fold scheme, where each line in a template structure written by electron beam lithography evolves into six metal lines.


title: How Much SETI Has Been Done? Finding Needles in the n-Dimensional Cosmic Haystack

abstract: Many articulations of the Fermi Paradox have as a premise, implicitly or explicitly, that humanity has searched for signs of extraterrestrial radio transmissions and concluded that there are few or no obvious ones to be found. Tarter et al. (2010) and others have argued strongly to the contrary: bright and obvious radio beacons might be quite common in the sky, but we would not know it yet because our search completeness to date is so low, akin to having searched a drinking glass's worth of seawater for evidence of fish in all of Earth's oceans. Here, we develop the metaphor of the multidimensional "Cosmic Haystack" through which SETI hunts for alien "needles" into a quantitative, eight-dimensional model and perform an analytic integral to compute the fraction of this haystack that several large radio SETI programs have collectively examined. Although this model haystack has many qualitative differences from the Tarter et al. (2010) haystack, we conclude that the fraction of it searched to date is also very small: similar to the ratio of the volume of a large hot tub or small swimming pool to that of the Earth's oceans. With this article we provide a Python script to calculate haystack volumes for future searches and for similar haystacks with different boundaries. We hope this formalism will aid in the development of a common parameter space for the computation of upper limits and completeness fractions of search programs for radio and other technosignatures.


title: Coarse-graining via the fluctuation-dissipation theorem and large-deviation theory

abstract: The fluctuation-dissipation theorem is a central result in statistical mechanics and is usually formulated for systems described by diffusion processes. In this paper, we propose a generalization for a wider class of stochastic processes, namely the class of Markov processes that satisfy detailed balance and a large-deviation principle. The generalized fluctuation-dissipation theorem characterizes the deterministic limit of such a Markov process as a generalized gradient flow, a mathematical tool to model a purely irreversible dynamics via a dissipation potential and an entropy function: these are expressed in terms of the large-deviation dynamic rate function of the Markov process and its stationary distribution. We exploit the generalized fluctuation-dissipation theorem to develop a new method of coarse-graining and test it in the context of the passage from the diffusion in a double-well potential to the jump process that describes the simple reaction $A \rightleftarrows B$ (Kramers' escape problem).


title: Distributionally Robust Chance Constrained Optimal Power Flow Assuming Unimodal Distributions with Misspecified Modes

abstract: Chance constrained optimal power flow (CC-OPF) formulations have been proposed to minimize operational costs while controlling the risk arising from uncertainties like renewable generation and load consumption. To solve CC-OPF, we often need access to the (true) joint probability distribution of all uncertainties, which is rarely known in practice. A solution based on a biased estimate of the distribution can result in poor reliability. To overcome this challenge, recent work has explored distributionally robust chance constraints, in which the chance constraints are satisfied over a family of distributions called the ambiguity set. Commonly, ambiguity sets are only based on moment information (e.g., mean and covariance) of the random variables; however, specifying additional characteristics of the random variables reduces conservatism and cost. Here, we consider ambiguity sets that additionally incorporate unimodality information. In practice, it is difficult to estimate the mode location from the data and so we allow it to be potentially misspecified. We formulate the problem and derive a separation-based algorithm to efficiently solve it. Finally, we evaluate the performance of the proposed approach on a modified IEEE-30 bus network with wind uncertainty and compare with other distributionally robust approaches. We find that a misspecified mode significantly affects the reliability of the solution and the proposed model demonstrates a good trade-off between cost and reliability.


title: Shape of Dark Matter Haloes in the Illustris Simulation: Effects of Baryons

abstract: We study the effect of baryonic processes on the shapes of dark matter (DM) haloes from Illustris, a suite of hydrodynamical (Illustris) and DM-only (Illustris-Dark) cosmological simulations performed with the moving-mesh code {\sc arepo}. DM halo shapes are determined using an iterative method based on the inertia tensor for a wide range of $z=0$ masses ($M_{200} = 1 \times 10^{11} - 3 \times 10^{14} M_\odot$). Convergence tests shows that the local DM shape profiles are converged only for $r > 9\epsilon$, $\epsilon$ being the Plummer-equivalent softening length, larger than expected. Haloes from non-radiative simulations (i.e. neglecting radiative processes, star formation, and feedback) exhibit no alteration in shapes from their DM-only counterparts: thus moving-mesh hydrodynamics alone is insufficient to cause differences in DM shapes. With the full galaxy-physics implementation, condensation of baryons results in significantly rounder and more oblate haloes, with the median minor-to-major axis ratio $\left< s\equiv c/a \right> \approx 0.7$, almost invariant throughout the halo and across halo masses. This somewhat improves the agreement between simulation predictions and observational estimates of the Milky Way halo shape. Consistently, the velocity anisotropy of DM is also reduced in Illustris, across halo masses and radii. Within the inner halo ($r=0.15 R_{200}$), both $s$ and $q$ (intermediate-to-major axis ratio) exhibit non-monotonicity with galaxy mass, peaking at $m_* \approx 10^{10.5-11} M_\odot$, which we find is due to the strong dependence of inner halo shape with galaxy formation efficiency. Baryons in Illustris affect the correlation of halo shape with halo properties, leading to a positive correlation of sphericity of MW-mass haloes with halo formation time and concentration, the latter being mildly more pronounced than in Illustris-Dark.


title: Audio Based Disambiguation Of Music Genre Tags

abstract: In this paper, we propose to infer music genre embeddings from audio datasets carrying semantic information about genres. We show that such embeddings can be used for disambiguating genre tags (identification of different labels for the same genre, tag translation from a tag system to another, inference of hierarchical taxonomies on these genre tags). These embeddings are built by training a deep convolutional neural network genre classifier with large audio datasets annotated with a flat tag system. We show empirically that they makes it possible to retrieve the original taxonomy of a tag system, spot duplicates tags and translate tags from a tag system to another.


title: MTLE: A Multitask Learning Encoder of Visual Feature Representations for Video and Movie Description

abstract: Learning visual feature representations for video analysis is a daunting task that requires a large amount of training samples and a proper generalization framework. Many of the current state of the art methods for video captioning and movie description rely on simple encoding mechanisms through recurrent neural networks to encode temporal visual information extracted from video data. In this paper, we introduce a novel multitask encoder-decoder framework for automatic semantic description and captioning of video sequences. In contrast to current approaches, our method relies on distinct decoders that train a visual encoder in a multitask fashion. Our system does not depend solely on multiple labels and allows for a lack of training data working even with datasets where only one single annotation is viable per video. Our method shows improved performance over current state of the art methods in several metrics on multi-caption and single-caption datasets. To the best of our knowledge, our method is the first method to use a multitask approach for encoding video features. Our method demonstrates its robustness on the Large Scale Movie Description Challenge (LSMDC) 2017 where our method won the movie description task and its results were ranked among other competitors as the most helpful for the visually impaired.


title: DPPy: Sampling DPPs with Python

abstract: Determinantal point processes (DPPs) are specific probability distributions over clouds of points that are used as models and computational tools across physics, probability, statistics, and more recently machine learning. Sampling from DPPs is a challenge and therefore we present DPPy, a Python toolbox that gathers known exact and approximate sampling algorithms for both finite and continuous DPPs. The project is hosted on GitHub and equipped with an extensive documentation.


title: Quasilocal Smarr relation for an asymptotically flat spacetime

abstract: We investigate the thermodynamics of Einstein-Maxwell(-Dilaton) theory for an asymptotically flat spacetime in a quasilocal frame. We firstly define a quasilocal thermodynamic potential via the Euclidean on-shell action and formulate a quasilocal Smarr relation from Eulerian theorem. Then we calculate quasilocal energy and surface pressure by employing Brown-York quasilocal method along with Mann-Marolf counterterm and find entropy from the quasilocal thermodynamic potential. These quasilocal variables are consistent with Tolman temperature and the entropy in a quasilocal frame turns out to be same as the Bekenstein-Hawking entropy. As a result, we found that a surface pressure term and its conjugate variable, a quasilocal area, do not participate in a quasilocal thermodynamic potential, but should present in a quasilocal Smarr relation and the quasilocal first law of black hole thermodynamics. For dyonic black hole solutions having dynamic dilaton field, non-trivial dilaton contribution should take part in the quasilocal first law but not in the quasilocal Smarr relation.


title: Bayesian functional optimisation with shape prior

abstract: Real world experiments are expensive, and thus it is important to reach a target in minimum number of experiments. Experimental processes often involve control variables that changes over time. Such problems can be formulated as a functional optimisation problem. We develop a novel Bayesian optimisation framework for such functional optimisation of expensive black-box processes. We represent the control function using Bernstein polynomial basis and optimise in the coefficient space. We derive the theory and practice required to dynamically adjust the order of the polynomial degree, and show how prior information about shape can be integrated. We demonstrate the effectiveness of our approach for short polymer fibre design and optimising learning rate schedules for deep networks.


title: Search for High-Mass Protostellar Objects in Cold IRAS Sources

abstract: We present the results of CS J=2-1 mapping observations towards 39 massive star-forming regions selected from the previous CO line survey of cold IRAS sources with high-velocity CO flows along the Galactic plane (Yang et al. 2002). All sources are detected in CS J=2-1 showing the existence of CS clumps around the IRAS sources. However, one-third of the sources are not deeply embedded in the dense clumps by comparison of the central powering IRAS sources and the morphologies of CS clumps. Physical parameters of the dense molecular clumps are presented. We have identified 12 high-mass protostellar object (HMPO) candidates by checking the association between the dense cores and the IRAS sources, the detection of water maser, and the radio properties towards the IRAS sources. We find that the HMPO sources are characterized by low FIR luminosity to virial mass ratios since they are in very early evolutionary stages when the massive protostars have not reached their full luminosities, which are typical for zero-age main sequence stars. Large turbulent motion in the HMPO sources may be largely due to the large kinetic energy ejected by the central protostars formed in the dense clumps. However, alternative means or undetected outflows may also be responsible for the turbulence in the clumps.


title: A Novel Warehouse Multi-Robot Automation System with Semi-Complete and Computationally Efficient Path Planning and Adaptive Genetic Task Allocation Algorithms

abstract: We consider the problem of warehouse multi-robot automation system in discrete-time and discrete-space configuration with focus on the task allocation and conflict-free path planning. We present a system design where a centralized server handles the task allocation and each robot performs local path planning distributively. A genetic-based task allocation algorithm is firstly presented, with modification to enable heuristic learning. A semi-complete potential field based local path planning algorithm is then proposed, named the recursive excitation/relaxation artificial potential field (RERAPF). A mathematical proof is also presented to show the semi-completeness of the RERAPF algorithm. The main contribution of this paper is the modification of conventional artificial potential field (APF) to be semi-complete while computationally efficient, resolving the traditional issue of incompleteness. Simulation results are also presented for performance evaluation of the proposed path planning algorithm and the overall system.


title: Effects of the number of layers on the vibrational, electronic and optical properties of alpha lead oxide

abstract: We have investigated the effects of number of layers on the structural, vibrational, electronic and optical properties of $\alpha$-PbO using first principles calculations. Our theoretical calculations have shown that four Raman active modes of $\alpha$-PbO tend to red-shift from bulk to monolayer due to decreasing of force constants and increasing of bond lengths. It has been shown that while bulk and multilayer $\alpha$-PbO have an indirect band gap, monolayer form has a direct band gap value of 2.59 eV. Although lead atoms have 5d states, spin-orbit coupling does not significantly affect the band structure of $\alpha$-PbO. The computed cleavage energy value (0.67 J/m$^2$) confirms that monolayer PbO can be easily obtained from its bulk counterpart by exfoliation methods. In addition to the band structure, we also calculated the optical properties and absorbed photon flux $J_{abs}$ of $\alpha$-PbO structures to investigate the possibility of solar absorption. Our calculations reveal that while monolayer and bilayer PbO have relatively larger band gaps and lower absorption coefficients, their $J_{abs}$ values are not ideal for solar absorption devices. In contrast, the multilayer and bulk phases of the $\alpha$-PbOs show good overlap with the solar spectrum and yield high electrical current values. Our calculations have indicated that ultrathin films of $\alpha$-PbO (such as 3nm thickness) could be excellent candidates for solar cells. We believe that our work can be utilized to improve electronic and optical devices based on lead oxide structures.


title: Stability of the Cosine-Sine Functional Equation on amenable groups

abstract: In this paper we establish the stability of the functional equation \begin{equation*}f(xy)=f(x)g(y)+g(x)f(y)+h(x)h(y),\;x,y\in G,\end{equation*} where $G$ is an amenable group.


title: Markov selection for the stochastic compressible Navier--Stokes system

abstract: We analyze the Markov property of solutions to the compressible Navier--Stokes system perturbed by a general multiplicative stochastic forcing. We show the existence of an almost sure Markov selection to the associated martingale problem. Our proof is based on the abstract framework introduced in [F. Flandoli, M. Romito: Markov selections for the 3D stochastic Navier--Stokes equations. Probab. Theory Relat. Fields 140, 407--458. (2008)]. A major difficulty arises from the fact, different from the incompressible case, that the velocity field is not continuous in time. In addition, it cannot be recovered from the variables whose time evolution is described by the Navier--Stokes system, namely, the density and the momentum. We overcome this issue by introducing an auxiliary variable into the Markov selection procedure.


title: Black hole memory effect

abstract: We compute the memory effect produced at the black hole horizon by a transient gravitational shockwave. As shown by Hawking, Perry, and Strominger (HPS) such a gravitational wave produces a deformation of the black hole geometry which from future null infinity is seen as a Bondi-Metzner-Sachs (BMS) supertranslation. This results in a diffeomorphic but physically distinct geometry which differs from the original black hole by their charges at infinity. Here we give the complementary description of this physical process in the near-horizon region as seen by an observer hovering just outside the event horizon. From this perspective, in addition to a supertranslation the shockwave also induces a horizon superrotation. We compute the associated superrotation charge and show that its form agrees with the one obtained by HPS at infinity. In addition, there is a supertranslation contribution to the horizon charge, which measures the entropy change in the process. We then turn to electrically and magnetically charged black holes and generalize the near-horizon asymptotic symmetry analysis to Einstein-Maxwell theory. This reveals an additional infinite-dimensional current algebra that acts non-trivially on the horizon superrotations. Finally, we generalize the black hole memory effect to Reissner-Nordstr\"{o}m black holes.


title: LFRic: Meeting the challenges of scalability and performance portability in Weather and Climate models

abstract: This paper describes LFRic: the new weather and climate modelling system being developed by the UK Met Office to replace the existing Unified Model in preparation for exascale computing in the 2020s. LFRic uses the GungHo dynamical core and runs on a semi-structured cubed-sphere mesh. The design of the supporting infrastructure follows object orientated principles to facilitate modularity and the use of external libraries where possible. In particular, a `separation of concerns' between the science code and parallel code is imposed to promote performance portability. An application called PSyclone, developed at the STFC Hartree centre, can generate the parallel code enabling deployment of a single source science code onto different machine architectures. This paper provides an overview of the scientific requirement, the design of the software infrastructure, and examples of PSyclone usage. Preliminary performance results show strong scaling and an indication that hybrid MPI/OpenMP performs better than pure MPI.


title: Long-term analysis of exponential integrators for highly oscillatory conservative systems

abstract: In this paper, we investigate the long-time near-conservations of energy and kinetic energy by the widely used exponential integrators to highly oscillatory conservative systems. The modulated Fourier expansions of two kinds of exponential integrators have been constructed and the long-time numerical conservations of energy and kinetic energy are obtained by deriving two almost-invariants of the expansions. Practical examples of the methods are given and the theoretical results are confirmed and demonstrated by a numerical experiment.


title: Towards Dialogue-based Navigation with Multivariate Adaptation driven by Intention and Politeness for Social Robots

abstract: Service robots need to show appropriate social behaviour in order to be deployed in social environments such as healthcare, education, retail, etc. Some of the main capabilities that robots should have are navigation and conversational skills. If the person is impatient, the person might want a robot to navigate faster and vice versa. Linguistic features that indicate politeness can provide social cues about a person's patient and impatient behaviour. The novelty presented in this paper is to dynamically incorporate politeness in robotic dialogue systems for navigation. Understanding the politeness in users' speech can be used to modulate the robot behaviour and responses. Therefore, we developed a dialogue system to navigate in an indoor environment, which produces different robot behaviours and responses based on users' intention and degree of politeness. We deploy and test our system with the Pepper robot that adapts to the changes in user's politeness.


title: Holographic quantization of linearized higher-spin gravity in the de Sitter causal patch

abstract: We study the dS/CFT duality between minimal type-A higher-spin gravity and the free Sp(2N) vector model. We consider the bulk spacetime as "elliptic" de Sitter space dS_4/Z_2, in which antipodal points have been identified. We apply a technique from arXiv:1509.05890, which extracts the quantum-mechanical commutators (or Poisson brackets) of the linearized bulk theory in an *observable patch* of dS_4/Z_2 directly from the boundary 2-point function. Thus, we construct the Lorentzian commutators of the linearized bulk theory from the Euclidean CFT. In the present paper, we execute this technique for the entire higher-spin multiplet, using a higher-spin-covariant language, which provides a promising framework for the future inclusion of bulk interactions. Aside from its importance for dS/CFT, our construction of a Hamiltonian structure for a bulk causal region should be of interest within higher-spin theory itself. The price we pay is a partial symmetry breaking, from the full dS group (and its higher-spin extension) to the symmetry group of an observable patch. While the boundary field theory plays a role in our arguments, the results can be fully expressed within a boundary *particle mechanics*. Bulk fields arise from this boundary mechanics via a version of second quantization.


title: A novel technique to solve the fuzzy system of equations

abstract: The aim of this research is to apply a novel technique based on the embedding method to solve the n*n fuzzy system of linear equations (FSLEs). By using this method, the strong fuzzy number solutions of FSLEs can be obtained by transforming the n*n FSLE to the crisp system. In this paper, Ezzati's method to solve the FSLEs is modified and improved. Several theorems are proved to show the number of operations for presented method are less than the methods of Friedman and Ezzati. In order to show the advantages of scheme, two applicable algorithms are presented and several examples are solved by applying them. Also, some graphs of obtained results are demonstrated which show the solutions are in the fuzzy form.


title: Grid invariants in universally tight lens spaces

abstract: We define combinatorial invariants of Legendrian and transverse links in universally tight lens spaces using grid diagrams, generalizing [OST08] and prove that they are equivalent to the invariants defined in [BVVV13] and [LOSS09]. We use these combinatorial invariants to characterize index one grid diagrams for knots in lens spaces which admit surgeries to the 3-sphere and discuss a potential application to the Berge conjecture.


title: Self-calibration method for II and GI types of intrinsic alignments of galaxies

abstract: We introduce a self-calibration method that can be applied to the intrinsic ellipticity--intrinsic ellipticity (II) and gravitational shear -- intrinsic ellipticity (GI) types of intrinsic alignment of galaxies. The method combines previous self-calibration techniques with modifications to one of them in order to use auto-spectra in addition to cross-spectra between redshift bins. This allows one to use the self-calibration while preserving all the constraining power of cosmic shear from surveys. We show that the new method provides more flexibility in using various redshift bin widths. We perform cosmological parameter constraint forecast when this method is applied to the Large Synoptic Survey Telescope (LSST). Compared to the original self-calibration, we find that the new method provides further significant reduction in any residual shift in the cosmological parameters (e.g. factors of $2-4$ for the dark energy equation of state) which is promising for accurate cosmology.


title: A closer look at the transition between fully convective and partly radiative low mass stars

abstract: Recently, Jao et al. (2018) discovered a gap in the mid-M dwarfs main sequence revealed by the analysis of Gaia data Release 2. They suggested the feature is linked to the onset of full convection in M dwarfs. Following the announcement of this discovery, MacDonald & Gizis (2018) proposed an explanation based on standard stellar evolution models. In this paper we re-examine the explanation suggested by MacDonald & Gizis (2018). We confirm that nuclear burning and mixing process of $^3$He provide the best explanation for the observed feature. We also find that a change in the energy transport from convection to radiation does not induce structural changes that could be visible. Regarding the very details of the process, however, we disagree with MacDonald & Gizis (2018) and propose a different explanation.


title: On the reaction rate constants that enable multistationarity in the two-site phosphorylation cycle

abstract: Parametrized polynomial ordinary differential equation systems are broadly used for modeling, specially in the study of biochemical reaction networks under the assumption of mass-action kinetics. Understanding the qualitative behavior of the solutions with respect to the parameter values gives rise to complex problems within real algebraic geometry, concerning the study of the signs of multivariate polynomials over the positive orthant. In this work we provide further insight into the number of positive steady states of a benchmark model, namely the two-site phosphorylation cycle. In particular, we provide new conditions on the reaction rate constants for the existence of one or three positive steady states, partially filling a gap left in previous works.


title: Music Mood Detection Based On Audio And Lyrics With Deep Neural Net

abstract: We consider the task of multimodal music mood prediction based on the audio signal and the lyrics of a track. We reproduce the implementation of traditional feature engineering based approaches and propose a new model based on deep learning. We compare the performance of both approaches on a database containing 18,000 tracks with associated valence and arousal values and show that our approach outperforms classical models on the arousal detection task, and that both approaches perform equally on the valence prediction task. We also compare the a posteriori fusion with fusion of modalities optimized simultaneously with each unimodal model, and observe a significant improvement of valence prediction. We release part of our database for comparison purposes.


title: Dolbeault cohomologies of blowing up complex manifolds II: bundle-valued case

abstract: We use a sheaf-theoretic approach to obtain a blow-up formula for Dolbeault cohomology groups with values in the holomorphic vector bundle over a compact complex manifold. As applications, we present several positive (or negative) examples associated to the vanishing theorems of Girbau, Kawamata-Viehweg and Green-Lazarsfeld in a uniform manner and study the blow-up invariance of some classical holomorphic invariants.


title: Dependence of alpha-particle-driven Alfv\'en eigenmode linear stability on device magnetic field strength and consequences for next-generation tokamaks

abstract: Recently-proposed tokamak concepts use magnetic fields up to 12 T, far higher than in conventional devices, to reduce size and cost. Theoretical and computational study of trends in plasma behavior with increasing field strength is critical to such proposed devices. This paper considers trends in Alfv\'en eigenmode (AE) stability. Energetic particles, including alphas from D-T fusion, can destabilize AEs, possibly causing loss of alpha heat and damage to the device. AEs are sensitive to device magnetic field via the field dependence of resonances, alpha particle beta, and alpha orbit width. We describe the origin and effect of these dependences analytically and by using recently-developed numerical techniques (Rodrigues et al. 2015 Nucl. Fusion 55 083003). The work suggests high-field machines where fusion-born alphas are sub-Alfv\'enic or nearly sub-Alfv\'enic may partially cut off AE resonances, reducing growth rates of AEs and the energy of alphas interacting with them. High-field burning plasma regimes have non-negligible alpha particle beta and AE drive, but faster slowing down time, provided by high electron density, and higher field strength reduces this drive relative to low-field machines with similar power densities. The toroidal mode number of the most unstable modes will tend to be higher in high magnetic field devices. The work suggests that high magnetic field devices have unique, and potentially advantageous, AE instability properties at both low and high densities.


title: Feedback Control of a Cassie Bipedal Robot: Walking, Standing, and Riding a Segway

abstract: The Cassie bipedal robot designed by Agility Robotics is providing academics a common platform for sharing and comparing algorithms for locomotion, perception, and navigation. This paper focuses on feedback control for standing and walking using the methods of virtual constraints and gait libraries. The designed controller was implemented six weeks after the robot arrived at the University of Michigan and allowed it to stand in place as well as walk over sidewalks, grass, snow, sand, and burning brush. The controller for standing also enables the robot to ride a Segway. A model of the Cassie robot has been placed on GitHub and the controller will also be made open source if the paper is accepted.


title: Measurement of photon-jet transverse momentum correlations in 5.02 TeV Pb+Pb and $pp$ collisions with ATLAS

abstract: Jets created in association with a photon can be used as a calibrated probe to study energy loss in the medium created in nuclear collisions. Measurements of the transverse momentum balance between isolated photons and inclusive jets are presented using integrated luminosities of 0.49 nb$^{-1}$ of Pb+Pb collision data at $\sqrt{s_\mathrm{NN}}=5.02$ TeV and 25 pb$^{-1}$ of $pp$ collision data at $\sqrt{s}=5.02$ TeV recorded with the ATLAS detector at the LHC. Photons with transverse momentum $63.1 < p_\mathrm{T}^{\gamma} < 200$ GeV and $\left|\eta^{\gamma}\right| < 2.37$ are paired inclusively with all jets in the event that have $p_\mathrm{T}^\mathrm{jet} > 31.6$ GeV and pseudorapidity $\left|\eta^\mathrm{jet}\right| < 2.8$. The transverse momentum balance given by the jet-to-photon $p_\mathrm{T}$ ratio, $x_\mathrm{J\gamma}$, is measured for pairs with azimuthal opening angle $\Delta\phi > 7\pi/8$. Distributions of the per-photon jet yield as a function of $x_\mathrm{J\gamma}$, $(1/N_\gamma)(\mathrm{d}N/\mathrm{d}x_\mathrm{J\gamma})$, are corrected for detector effects via a two-dimensional unfolding procedure and reported at the particle level. In $pp$ collisions, the distributions are well described by Monte Carlo event generators. In Pb+Pb collisions, the $x_\mathrm{J\gamma}$ distribution is modified from that observed in $pp$ collisions with increasing centrality, consistent with the picture of parton energy loss in the hot nuclear medium. The data are compared with a suite of energy-loss models and calculations.


title: Triboelectrification of KCl and ZnS particles in approximated exoplanet environments

abstract: When mobilized, granular materials become charged as grains undergo collisions and frictional interactions. On Earth, this process, known as triboelectrification, has been recognized in volcanic plumes and sandstorms. Yet, frictional charging almost certainly exists on other worlds, both in our own Solar System (such as Mars, the Moon, and Venus) and exosolar planets. Indeed, observations suggest that numerous planets in the galaxy are enshrouded by optically-thick clouds or hazes. Triboelectric charging within these clouds may contribute to global electric circuits of these worlds, providing mechanisms to generate lightning, drive chemical processes in the atmospheres, and, perhaps, influence habitability. In this work, we explore the frictional electrification of potassium chloride and zinc sulfide, two substances proposed to make up the clouds of giant exo-planets with >50x solar metallicities, including the widely-studied super-Earth GJ 1214b, super-earth HD 97658b, Neptune-sized GJ 436b, and hot- Jupiter WASP-31b. We find that both materials become readily electrified when mobilized, attaining charge densities similar to those found on volcanic ash particles. Thus, if these worlds do indeed host collections of mineral particles in their atmospheres, these clouds are likely electrified and may be capable of producing lightning or corona discharge.


title: Modeling Online Discourse with Coupled Distributed Topics

abstract: In this paper, we propose a deep, globally normalized topic model that incorporates structural relationships connecting documents in socially generated corpora, such as online forums. Our model (1) captures discursive interactions along observed reply links in addition to traditional topic information, and (2) incorporates latent distributed representations arranged in a deep architecture, which enables a GPU-based mean-field inference procedure that scales efficiently to large data. We apply our model to a new social media dataset consisting of 13M comments mined from the popular internet forum Reddit, a domain that poses significant challenges to models that do not account for relationships connecting user comments. We evaluate against existing methods across multiple metrics including perplexity and metadata prediction, and qualitatively analyze the learned interaction patterns.


title: Rapid Electron Backscatter Diffraction Mapping: Painting by Numbers

abstract: Microstructure characterisation has been greatly enhanced through the use of electron backscatter diffraction (EBSD), where rich maps are generated through analysis of the crystal phase and orientation in the scanning electron microscope (SEM). Conventional EBSD analysis involves raster scanning of the electron beam and the serial analysis of each diffraction pattern in turn. For grain shape, crystallographic texture, and microstructure analysis this can be inefficient. In this work, we present Rapid EBSD, a data fusion approach combining forescatter electron (FSE) imaging with static sparse sampling of EBSD patterns. We segment the FSE image into regions of similar colour (i.e. phase and crystal orientation) and then collect representative EBSD data for each segmented region. This enables microstructural assessment to be performed at the spatial resolution of the (fast) FSE imaging whilst including orientation and phase information from EBSD analysis of representative points. We demonstrate the Rapid EBSD technique on samples of a cobalt based superalloy and a strained dual phase titanium alloy, comparing the results with conventional analysis. Rapid EBSD is advantageous for assessing grain size distributions in time-limited experiments.


title: A question of Norton-Sullivan in the analytic case

abstract: In 1996, A. Norton and D. Sullivan asked the following question: If $f:\mathbb{T}^2\rightarrow\mathbb{T}^2$ is a diffeomorphism, $h:\mathbb{T}^2\rightarrow\mathbb{T}^2$ is a continuous map homotopic to the identity, and $h f=T_{\rho} h$ where $\rho\in\mathbb{R}^2$ is a totally irrational vector and $T_{\rho}:\mathbb{T}^2\rightarrow\mathbb{T}^2,\, z\mapsto z+\rho$ is a translation, are there natural geometric conditions (e.g. smoothness) on $f$ that force $h$ to be a homeomorphism? In [ J. Wang and Z. Zhang, GAFA 2018 ], the first author and Z. Zhang gave a negative answer to the above question in the $C^{\infty}$ category: In general, not even the infinite smoothness condition can force $h$ to be a homeomorphism. In this article, we give a negative answer in the $C^{\omega}$ category: We construct a real-analytic conservative and minimal totally irrational pseudo-rotation of $\mathbb{T}^2$ that is semi-conjugate to a translation but not conjugate to a translation, which simultaneously answers a question raised in [ J. Wang and Z. Zhang, GAFA 2018 ].


title: Cosmic ray Nickel nuclei spectrum by the NUCLEON experiment

abstract: The NUCLEON experiment is designed to measure chemical composition of cosmic rays with charges from Z=1 to 30 in an energy region from 5*10^11 to 10^15 eV. In this article the data analysis algorithm and spectra of Ni and Fe nuclei, measured in the NUCLEON experiment, are presented.


title: Isocurvature Perturbations of Dark Energy and Dark Matter from the Swampland Conjecture

abstract: We point out that the recently proposed Swampland conjecture on the potential gradient can lead to isocurvature perturbations of dark energy, if the quintessence field acquires large quantum fluctuations during high-scale inflation preferred by the conjecture. Also, if the quintessence field is coupled to a dark sector that contains dark matter, isocurvature perturbation of dark matter is similarly induced. Both isocurvature perturbations can be suppressed if the quintessence potential allows a tracker solution in the early Universe. We find that a vector field of mass $\lesssim {\cal O}(1)$ meV is an excellent dark matter candidate in this context, not only because the right abundance is known to be produced by quantum fluctuations during high-scale inflation without running afoul of isocurvature bounds, but also because its coupling to the quintessence does not spoil the flatness of the potential.


title: A Physical Perspective on Control Points and Polar Forms: B\'ezier Curves, Angular Momentum and Harmonic Oscillators

abstract: Bernstein polynomials and B\'ezier curves play an important role in computer-aided geometric design and numerical analysis, and their study relates to mathematical fields such as abstract algebra, algebraic geometry and probability theory. We describe a theoretical framework that incorporates the different aspects of the Bernstein-B\'ezier theory, based on concepts from theoretical physics. We relate B\'ezier curves to the theory of angular momentum in both classical and quantum mechanics, and describe physical analogues of various properties of B\'ezier curves -- such as their connection with polar forms -- in the context of quantum spin systems. This previously unexplored relationship between geometric design and theoretical physics is established using the mathematical theory of Hamiltonian mechanics and geometric quantization. An alternative description of spin systems in terms of harmonic oscillators serves as a physical analogue of P\'olya's urn models for B\'ezier curves. We relate harmonic oscillators to Poisson curves and the analytical blossom as well. We present an overview of the relevant mathematical and physical concepts, and discuss opportunities for further research.


title: Time-varying Projected Dynamical Systems with Applications to Feedback Optimization of Power Systems

abstract: This paper is concerned with the study of continuous-time, non-smooth dynamical systems which arise in the context of time-varying non-convex optimization problems, as for example the feedback-based optimization of power systems. We generalize the notion of projected dynamical systems to time-varying, possibly non-regular, domains and derive conditions for the existence of so-called Krasovskii solutions. The key insight is that for trajectories to exist, informally, the time-varying domain can only contract at a bounded rate whereas it may expand discontinuously. This condition is met, in particular, by feasible sets delimited via piecewise differentiable functions under appropriate constraint qualifications. To illustrate the necessity and usefulness of such a general framework, we consider a simple yet insightful power system example, and we discuss the implications of the proposed conditions for the design of feedback optimization schemes.


title: Critical role of phenyl substitution and catalytic substrate on the surface-assisted polymerization of dibromobianthracene derivatives

abstract: Understanding the nature and hierarchy of on surface reactions is a major chal- lenge for designing coordination and covalent nanostructures by means of multistep synthetic routes. In particular, intermediates and final products are hard to predict since reaction paths and their activation windows depend on the choice of both the molecular precursor design and the substrate. Here we report a systematic study of the effect of the catalytic metal surface to reveal how a single precursor can give rise to very distinct polymers that range from coordination and covalent non planar polymer chains of distinct chirality, to atomically precise graphene nanoribbons and nanoporous graphene. Our precursor consists on adding two phenyl substituents to 10,10'-dibromo 9,9'-bianthracene, a well-studied precursor in the on-surface synthesis of graphene nanoribbons. The critical role of the monomer design on the reaction paths is inferred from the fact that the phenyl substitution leads to very distinct products in each one of the studied metallic substrates.


title: Higgs bundles and geometric structures on manifolds

abstract: Geometric structures on manifolds became popular when Thurston used them in his work on the geometrization conjecture. They were studied by many people and they play an important role in higher Teichm\"uller theory. Geometric structures on a manifold are closely related with representations of the fundamental group and with flat bundles. Higgs bundles can be very useful in describing flat bundles explicitly, via solutions of Hitchin's equations. Baraglia has shown in his Ph.D. Thesis that Higgs bundles can also be used to construct geometric structures in some interesting cases. In this paper, we will explain the main ideas behind this theory and we will survey some recent results in this direction, which are joint work with Qiongling Li.


title: Interpretable Textual Neuron Representations for NLP

abstract: Input optimization methods, such as Google Deep Dream, create interpretable representations of neurons for computer vision DNNs. We propose and evaluate ways of transferring this technology to NLP. Our results suggest that gradient ascent with a gumbel softmax layer produces n-gram representations that outperform naive corpus search in terms of target neuron activation. The representations highlight differences in syntax awareness between the language and visual models of the Imaginet architecture.


title: Online control of the false discovery rate in biomedical research

abstract: Modern biomedical research frequently involves testing multiple related hypotheses, while maintaining control over a suitable error rate. In many applications the false discovery rate (FDR), which is the expected proportion of false positives among the rejected hypotheses, has become the standard error criterion. Procedures that control the FDR, such as the well-known Benjamini-Hochberg procedure, assume that all p-values are available to be tested at a single time point. However, this ignores the sequential nature of many biomedical experiments, where a sequence of hypotheses is tested without having access to future p-values or even the number of hypotheses. Recently, the first procedures that control the FDR in this online manner have been proposed by Javanmard and Montanari (Ann. Stat. 2018), and built upon by Ramdas et al. (NIPS 2017, ICML 2018). In this paper, we compare and contrast these proposed procedures, with a particular focus on the setting where the p-values are dependent. We also propose a simple modification of the procedures for when there is an upper bound on the number of hypotheses to be tested. Using comprehensive simulation scenarios and case studies, we provide recommendations for which procedures to use in practice for online FDR control.


title: Sectional monodromy groups of projective curves

abstract: Fix a degree $d$ projective curve $X \subset \mathbb{P}^r$ over an algebraically closed field $K$. Let $U \subset (\mathbb{P}^r)^*$ be a dense open subvariety such that every hyperplane $H \in U$ intersects $X$ in $d$ smooth points. Varying $H \in U$ produces the monodromy action $\varphi: \pi_1^{\text{\'et}}(U) \to S_d$. Let $G_X := \mathrm{im}(\varphi)$. The permutation group $G_X$ is called the sectional monodromy group of $X$. In characteristic zero $G_X$ is always the full symmetric group, but sectional monodromy groups in characteristic $p$ can be smaller. For a large class of space curves ($r \geqslant 3$) we classify all possibilities for the sectional monodromy group $G$ as well as the curves with $G_X=G$. We apply similar methods to study a particular family of rational curves in $\mathbb{P}^2$, which enables us to answer an old question about Galois groups of generic trinomials.


title: Generative Adversarial Network in Medical Imaging: A Review

abstract: Generative adversarial networks have gained a lot of attention in the computer vision community due to their capability of data generation without explicitly modelling the probability density function. The adversarial loss brought by the discriminator provides a clever way of incorporating unlabeled samples into training and imposing higher order consistency. This has proven to be useful in many cases, such as domain adaptation, data augmentation, and image-to-image translation. These properties have attracted researchers in the medical imaging community, and we have seen rapid adoption in many traditional and novel applications, such as image reconstruction, segmentation, detection, classification, and cross-modality synthesis. Based on our observations, this trend will continue and we therefore conducted a review of recent advances in medical imaging using the adversarial training scheme with the hope of benefiting researchers interested in this technique.


title: Time Synchronization in modular collaborative robots

abstract: A new generation of robot systems which are modular, flexible and safe for human-robot interaction are needed. Existing cobots seem to meet only the later and require a modular approach to improve their reconfigurability and interoperability. We propose a new sub-class of cobots named M-cobots which tackle these problems. In particular, we discuss the relevance of synchronization for these systems, analyze it and demonstrate how with a properly configured M-cobot, we are able to obtain a) distributed sub-microsecond clock synchronization accuracy among modules, b) timestamping accuracy of ROS 2.0 messages under 100 microseconds and c) millisecond-level end-to-end communication latencies, even when disturbed with networking overloads of up to 90% of the network capacity.


title: Evolving SDN for Low-Power IoT Networks

abstract: Software Defined Networking (SDN) offers a flexible and scalable architecture that abstracts decision making away from individual devices and provides a programmable network platform. However, implementing a centralized SDN architecture within the constraints of a low-power wireless network faces considerable challenges. Not only is controller traffic subject to jitter due to unreliable links and network contention, but the overhead generated by SDN can severely affect the performance of other traffic. This paper addresses the challenge of bringing high-overhead SDN architecture to IEEE 802.15.4 networks. We explore how traditional SDN needs to evolve in order to overcome the constraints of low-power wireless networks, and discuss protocol and architectural optimizations necessary to reduce SDN control overhead - the main barrier to successful implementation. We argue that interoperability with the existing protocol stack is necessary to provide a platform for controller discovery and coexistence with legacy networks. We consequently introduce {\mu}SDN, a lightweight SDN framework for Contiki, with both IPv6 and underlying routing protocol interoperability, as well as optimizing a number of elements within the SDN architecture to reduce control overhead to practical levels. We evaluate {\mu}SDN in terms of latency, energy, and packet delivery. Through this evaluation we show how the cost of SDN control overhead (both bootstrapping and management) can be reduced to a point where comparable performance and scalability is achieved against an IEEE 802.15.4-2012 RPL-based network. Additionally, we demonstrate {\mu}SDN through simulation: providing a use-case where the SDN configurability can be used to provide Quality of Service (QoS) for critical network flows experiencing interference, and we achieve considerable reductions in delay and jitter in comparison to a scenario without SDN.


title: Analyzing Uncertainty Matrices Associated with Multiple S-parameter Measurements

abstract: This paper presents a detailed analysis of uncertainty matrices (i.e. measurement covariance matrices) associated with multiple complex-valued microwave scattering S-parameter measurements. The analysis is based on forming combinations of (2 X 2) sub-matrices from selected elements in the uncertainty matrix and using these sub-matrices to define uncertainty ellipses in two-dimensional measurement planes. These uncertainty ellipses can be used to assess the quality, accuracy and validity of the measurements. A simple example is given using measurements made on a waveguide mismatched line section at frequencies from 75 GHz to 110 GHz. The analysis is also useful when multiple S-parameter measurements are used as input quantities for measurement models used to determine other measurement quantities (i.e. as output quantities).


title: Periodic Functions, Lattices and Their Projections

abstract: Functions whose symmetries form a crystallographic group in particular have a lattice of periods, and the set of their level curves forms a periodic pattern. We show how after projecting these functions, one obtains new functions with a lattice of periods that is not the projection of the initial lattice. We also characterise all the crystallographic groups in three dimensions that are symmetry groups of patterns whose projections have periods in a given two-dimensional lattice. The particular example of patterns that after projection have a hexagonal lattice of periods is discussed in detail.


title: The Warm-starting Sequential Selection Problem and its Multi-round Extension

abstract: In the Sequential Selection Problem (SSP), immediate and irrevocable decisions need to be made as candidates randomly arrive for a job interview. Standard SSP variants, such as the well-known secretary problem, begin with an empty selection set (cold-start) and perform the selection process once over a single candidate set (single-round). In this paper we address these two limitations. First, we introduce the novel Warm-starting SSP (WSSP) setting which considers at hand a reference set, a set of previously selected items of a given quality, and tries to update optimally that set by (re-)assigning each job at most once. We adopt a cutoff-based approach to optimize a rank-based objective function over the final assignment of the jobs. In our technical contribution, we provide analytical results regarding the proposed WSSP setting, we introduce the algorithm Cutoff-based Cost Minimization (CCM) (and the low failures-CCM, which is more robust to high rate of resignations) that adapts to changes in the quality of the reference set thanks to the translation method we propose. Finally, we implement and test CCM in a multi-round setting that is particularly interesting for real-world application scenarios.


title: Pricing American Options by Exercise Rate Optimization

abstract: We present a novel method for the numerical pricing of American options based on Monte Carlo simulation and the optimization of exercise strategies. Previous solutions to this problem either explicitly or implicitly determine so-called optimal exercise regions, which consist of points in time and space at which a given option is exercised. In contrast, our method determines the exercise rates of randomized exercise strategies. We show that the supremum of the corresponding stochastic optimization problem provides the correct option price. By integrating analytically over the random exercise decision, we obtain an objective function that is differentiable with respect to perturbations of the exercise rate even for finitely many sample paths. The global optimum of this function can be approached gradually when starting from a constant exercise rate. Numerical experiments on vanilla put options in the multivariate Black-Scholes model and a preliminary theoretical analysis underline the efficiency of our method, both with respect to the number of time-discretization steps and the required number of degrees of freedom in the parametrization of the exercise rates. Finally, we demonstrate the flexibility of our method through numerical experiments on max call options in the classical Black-Scholes model, and vanilla put options in both the Heston model and the non-Markovian rough Bergomi model.


title: An Inter-Laboratory Comparison of NVNA Measurements

abstract: A comparison of nonlinear vector network analyser (NVNA) measurements has been carried out involving four organisations (National Physical Laboratory, UK, University of Surrey, UK, Chalmers University of Technology, Sweden and Keysight Technologies, Denmark). Three nonlinear devices consisting of two amplifiers and a nonlinear verification device (NLVD) were measured by each of the organisations. Results are presented which show generally good agreement between the measurements and give some indication of the typical amount of variability to be expected in measurements of this type.


title: Occurrence of Mixed Phase in Bi0.5Sr0.5Mn0.9Cr0.1O3 bulk sample: Electron Paramagnetic Resonance and Magnetization Studies

abstract: We study the effects of 10% Cr substitution in Mn sites of Bi0.5Sr0.5MnO3 on the antiferromagnetic (AFM) (TN ~ 110 K) transition using structural, magnetic and electron paramagnetic resonance (EPR) techniques. Field cooled (FC) and zero field cooled (ZFC) magnetization measurements done from 400 K down to 4 K show that the compound is in the paramagnetic (PM) phase till 50 K where it undergoes a transition to a short range ferromagnetic phase (FM). Electron paramagnetic resonance measurements performed in the temperature range 300 K till 80 K conform with the magnetization measurements as symmetric signals are observed owing to the paramagnetic phase. Below 80 K, signals become asymmetric. Electron paramagnetic resonance intensity peaks at ~ 110 K, the decreasing intensity below this temperature confirming the presence of antiferromagnetism. We conclude that below 50 K the magnetization and EPR results are consistent with a cluster glass phase of BSMCO, where ferromagnetic clusters coexist with an antiferromagnetic background.


title: Collision of rarefaction waves in Bose-Einstein condensates

abstract: We consider the problem of expansion of Bose-Einstein condensate released form a box. On the contrary to the standard situation of release from a harmonic trap, in this case the dynamics is complicated by a process of collision of two rarefaction waves propagating to the center of the initially uniform distribution. Complete analytical solution of this problem is obtained by Riemann method in hydrodynamic dispersionless approximation and the results are compared with the exact numerical solution of the Gross-Pitaevskii equation.


title: Evolution of circumbinary protoplanetary disks with photoevaporative winds driven by External Far Ultraviolet Radiation

abstract: Lifetimes of protoplanetary disks (PPDs) are believed to be severely constrained by material depleting mechanisms, including photoevaporative winds due to the host star radiation or external radiation sources. Most previous studies focused on exploring the role of the winds in the exposed PPDs with a single star; however, exploring the evolution of the circumbinary disks with the photoevaporative winds driven by the host star radiation and external radiation sources deserves further investigation. In this study, we investigate the evolution of the circumbinary PPDs with the photoevaporative winds induced by external far ultraviolet (FUV) radiation field. We show that this mass-loss process can significantly constrain properties of a circumbinary PPD, including its lifetime, mass and radius. The lifetime of a circumbinary PPD, for instance, is found by a factor of about two longer than a similar circumstellar disk and this enhancement strongly depends on the viscosity parameter. But our model shows that viscosity dependence of the disk lifetime in the circumbinary case is more pronounced compared to the circumstellar case. We also show that dispersal of a circumbinary PPD occurs over a longer time as the disk temperature distribution becomes steeper. Our results also imply that dead zone in a photoevaporative circumbinary PPD extends over a larger radial range in comparison to a circumstellar disk counterpart. We also show that our calculations are in agreement with the observed circumbinary PPDs orbiting equal-mass binaries.


title: Emulation of CPU-demanding reactive transport models: comparison of Gaussian processes, polynomial chaos expansion and deep neural networks

abstract: This paper presents a detailed comparison between 3 methods for emulating CPU-intensive reactive transport models (RTMs): Gaussian processes (GPs), polynomial chaos expansion (PCE) and deep neural networks (DNNs). Besides direct emulation of the simulated uranium concentration time series, replacing the original RTM by its emulator is also investigated for global sensitivity analysis (GSA), uncertainty propagation (UP) and probabilistic calibration using Markov chain Monte Carlo (MCMC) sampling. The selected DNN is found to be superior to both GPs and PCE in reproducing the input - output behavior of the considered 8-dimensional and 13-dimensional CPU-intensive RTMs. Furthermore, the two used PCE variants: standard PCE and sparse PCE (sPCE) appear to always provide the least accuracy while not differing much in performance. As a consequence of its better emulation capabilities, the used DNN outperforms the two other methods for UP. In addition, DNNs and GPs offer equally good approximations to the true first-order and total-order Sobol sensitivity indices while PCE does somewhat less well. Most surprisingly, despite its superior emulation skills the DNN approach leads to the worst solution of the considered synthetic inverse problem which involves 1224 measurement data with low noise. This apparently contradicting behavior is at least partially due to the small but complicated deterministic noise that affects the DNN-based predictions. Indeed, this complex error structure can drive the emulated solutions far away from the true posterior distribution. Overall, our findings indicate that when the available training set is relatively small (75 - to 500 input - output examples) and fixed beforehand, DNNs can well emulate RTMs but are not suited to emulation-based inversion. In contrast, GPs perform fairly well across all considered tasks: direct emulation, GSA, UP, and inversion.


title: Clustering students' open-ended questionnaire answers

abstract: Open responses form a rich but underused source of information in educational data mining and intelligent tutoring systems. One of the major obstacles is the difficulty of clustering short texts automatically. In this paper, we investigate the problem of clustering free-formed questionnaire answers. We present comparative experiments on clustering ten sets of open responses from course feedback queries in English and Finnish. We also evaluate how well the main topics could be extracted from clusterings with the HITS algorithm. The main result is that, for English data, affinity propagation performed well despite frequent outliers and considerable overlapping between real clusters. However, for Finnish data, the performance was poorer and none of the methods clearly outperformed the others. Similarly, topic extraction was very successful for the English data but only satisfactory for the Finnish data. The most interesting discovery was that stemming could actually deteriorate the clustering quality significantly.


title: A Game-Theoretic Analysis of Shard-Based Permissionless Blockchains

abstract: Low transaction throughput and poor scalability are significant issues in public blockchain consensus protocols such as Bitcoins. Recent research efforts in this direction have proposed shard-based consensus protocols where the key idea is to split the transactions among multiple committees (or shards), which then process these shards or set of transactions in parallel. Such a parallel processing of disjoint sets of transactions or shards by multiple committees significantly improves the overall scalability and transaction throughout of the system. However, one significant research gap is a lack of understanding of the strategic behavior of rational processors within committees in such shard-based consensus protocols. Such an understanding is critical for designing appropriate incentives that will foster cooperation within committees and prevent free-riding. In this paper, we address this research gap by analyzing the behavior of processors using a game-theoretic model, where each processor aims at maximizing its reward at a minimum cost of participating in the protocol. We first analyze the Nash equilibria in an N-player static game model of the sharding protocol. We show that depending on the reward sharing approach employed, processors can potentially increase their payoff by unilaterally behaving in a defective fashion, thus resulting in a social dilemma. In order to overcome this social dilemma, we propose a novel incentive-compatible reward sharing mechanism to promote cooperation among processors. Our numerical results show that achieving a majority of cooperating processors (required to ensure a healthy state of the blockchain network) is easier to achieve with the proposed incentive-compatible reward sharing mechanism than with other reward sharing mechanisms.


title: Ordinal Synchronization: Using ordinal patterns to capture interdependencies between time series

abstract: We introduce Ordinal Synchronization ($OS$) as a new measure to quantify synchronization between dynamical systems. $OS$ is calculated from the extraction of the ordinal patterns related to two time series, their transformation into $D$-dimensional ordinal vectors and the adequate quantification of their alignment. $OS$ provides a fast and robust-to noise tool to assess synchronization without any implicit assumption about the distribution of data sets nor their dynamical properties, capturing in-phase and anti-phase synchronization. Furthermore, varying the length of the ordinal vectors required to compute $OS$ it is possible to detect synchronization at different time scales. We test the performance of $OS$ with data sets coming from unidirectionally coupled electronic Lorenz oscillators and brain imaging datasets obtained from magnetoencephalographic recordings, comparing the performance of $OS$ with other classical metrics that quantify synchronization between dynamical systems.


title: On generalized J\o{}rgensen inequality in ${\rm SL}(2, \mathbb C)$

abstract: Wang, Jiang and Cao have obtained a generalized version of the J\o{}rgensen inequality in Proc. Indian Acad. Sci. Math. Sci., 123(2):245--251, 2013, for two generator subgroups of ${\rm SL}(2, \mathbb C)$ where one of the generators is loxodromic. We prove that their inequality is strict.


title: Combinatorial and Structural Results for gamma-Psi-dimensions

abstract: This article deals with the generalization performance of margin multi-category classifiers, when minimal learnability hypotheses are made. In that context, the derivation of a guaranteed risk is based on the handling of capacity measures belonging to three main families: Rademacher/Gaussian complexities, metric entropies and scale-sensitive combinatorial dimensions. The scale-sensitive combinatorial dimensions dedicated to the classifiers of this kind are the gamma-Psi-dimensions. We introduce the combinatorial and structural results needed to involve them in the derivation of guaranteed risks and establish the corresponding upper bounds on the metric entropies and the Rademacher complexity. Two major conclusions can be drawn: 1. the gamma-Psi-dimensions always bring an improvement compared to the use of the fat-shattering dimension of the class of margin functions; 2. thanks to their capacity to take into account basic features of the classifier, they represent a promising alternative to performing the transition from the multi-class case to the binary one with covering numbers.


title: Accelerating Flash Calculation through Deep Learning Methods

abstract: In the past two decades, researchers have made remarkable progress in accelerating flash calculation, which is very useful in a variety of engineering processes. In this paper, general phase splitting problem statements and flash calculation procedures using the Successive Substitution Method are reviewed, while the main shortages are pointed out. Two acceleration methods, Newton's method and the Sparse Grids Method are presented afterwards as a comparison with the deep learning model proposed in this paper. A detailed introduction from artificial neural networks to deep learning methods is provided here with the authors' own remarks. Factors in the deep learning model are investigated to show their effect on the final result. A selected model based on that has been used in a flash calculation predictor with comparison with other methods mentioned above. It is shown that results from the optimized deep learning model meet the experimental data well with the shortest CPU time. More comparison with experimental data has been conducted to show the robustness of our model.


title: An Information Matrix Approach for State Secrecy

abstract: This paper studies the problem of remote state estimation in the presence of a passive eavesdropper. A sensor measures a linear plant's state and transmits it to an authorized user over a packet-dropping channel, which is susceptible to eavesdropping. Our goal is to design a coding scheme such that the eavesdropper cannot infer the plant's current state, while the user successfully decodes the sent messages. We employ a novel class of codes, termed State-Secrecy Codes, which are fast and efficient for dynamical systems. They apply linear time-varying transformations to the current and past states received by the user. In this way, they force the eavesdropper's information matrix to decrease with asymptotically the same rate as in the open-loop prediction case, i.e. when the eavesdropper misses all messages. As a result, the eavesdropper's minimum mean square error (mmse) for the unstable states grows unbounded, while the respective error for the stable states converges to the open-loop prediction one. These secrecy guarantees are achieved under minimal conditions, which require that, at least once, the user receives the corresponding packet while the eavesdropper fails to intercept it. Meanwhile, the user's estimation performance remains optimal. The theoretical results are illustrated in simulations.


title: Symmetric Shannon capacity is the independence number minus 1

abstract: A symmetric variant of Shannon capacity is defined and computed.


title: Efficient and Privacy-Preserving Ride Sharing Organization for Transferable and Non-Transferable Services

abstract: Ride-sharing allows multiple persons to share their trips together in one vehicle instead of using multiple vehicles. This can reduce the number of vehicles in the street, which consequently can reduce air pollution, traffic congestion and transportation cost. However, a ride-sharing organization requires passengers to report sensitive location information about their trips to a trip organizing server (TOS) which creates a serious privacy issue. In addition, existing ride-sharing schemes are non-flexible, i.e., they require a driver and a rider to have exactly the same trip to share a ride. Moreover, they are non-scalable, i.e., inefficient if applied to large geographic areas. In this paper, we propose two efficient privacy-preserving ride-sharing organization schemes for Non-transferable Ride-sharing Services (NRS) and Transferable Ride-sharing Services (TRS). In the NRS scheme, a rider can share a ride from its source to destination with only one driver whereas, in TRS scheme, a rider can transfer between multiple drivers while en route until he reaches his destination. In both schemes, the ride-sharing area is divided into a number of small geographic areas, called cells, and each cell has a unique identifier. Each driver/rider should encrypt his trip's data and send an encrypted ride-sharing offer/request to the TOS. In NRS scheme, Bloom filters are used to compactly represent the trip information before encryption. Then, the TOS can measure the similarity between the encrypted trips data to organize shared rides without revealing either the users' identities or the location information. In TRS scheme, drivers report their encrypted routes, an then the TOS builds an encrypted directed graph that is passed to a modified version of Dijkstra's shortest path algorithm to search for an optimal path of rides that can achieve a set of preferences defined by the riders.


title: Analog Coding Frame-work

abstract: Analog coding is a low-complexity method to combat erasures, based on linear redundancy in the signal space domain. Previous work examined "band-limited discrete Fourier transform (DFT)" codes for Gaussian channels with erasures or impulses. We extend this concept to source coding with "erasure side-information" at the encoder and show that the performance of band-limited DFT can be significantly improved using irregular spectrum, and more generally, using equiangular tight frames (ETF). Frames are overcomplete bases and are widely used in mathematics, computer science, engineering, and statistics since they provide a stable and robust decomposition. Design of frames with favorable properties of random subframes is motivated in variety of applications, including code-devision multiple access (CDMA), compressed sensing and analog coding. We present a novel relation between deterministic frames and random matrix theory. We show empirically that the MANOVA ensemble offers a universal description of the spectra of randomly selected subframes with constant aspect ratios, taken from deterministic near-ETFs. Moreover, we derive an analytic framework and bring a formal validation for some of the empirical results, specifically that the asymptotic form for the moments of high orders of subsets of ETF agree with that of MANOVA. Finally, when exploring over-complete bases, the Welch bound is a lower bound on the root mean square cross correlation between vectors. We extend the Welch bound to an erasure setting, in which a reduced frame, composed of a random subset of Bernoulli selected vectors, is of interest. The lower bound involves moment of the reduced frame, and it is tight for ETFs and asymptotically coincides with the MANOVA moments. This result offers a novel perspective on the superiority of ETFs over other frames.