title: Guessing Individual Sequences: Generating Randomized Guesses Using Finite-State Machines

abstract: Motivated by earlier results on universal randomized guessing, we consider an individual-sequence approach to the guessing problem: in this setting, the goal is to guess a secret, individual (deterministic) vector $x^n=(x_1,\ldots,x_n)$, by using a finite-state machine that sequentially generates randomized guesses from a stream of purely random bits. We define the finite-state guessing exponent as the asymptotic normalized logarithm of the minimum achievable moment of the number of randomized guesses, generated by any finite-state machine, until $x^n$ is guessed successfully. We show that the finite-state guessing exponent of any sequence is intimately related to its finite-state compressibility (due to Lempel and Ziv), and it is asymptotically achieved by the decoder of (a certain modified version of) the 1978 Lempel-Ziv data compression algorithm (a.k.a. the LZ78 algorithm), fed by purely random bits. The results are also extended to the case where the guessing machine has access to a side information sequence, $y^n=(y_1,\ldots,y_n)$, which is also an individual sequence.


title: Perfect anomalous splitter by acoustic meta-grating

abstract: As an inversely designed artificial device, metasurface usually means densely arranged meta-atoms with complex substructures. In acoustics, those meta-atoms are usually constructed by multi-folded channels or multi-connected cavities of deep sub-wavelength feature, which limits their implementation in pragmatic applications. We propose here a comprehensive concept of a perfect anomalous splitter based on an acoustic meta-grating. The beam splitter is designed by etching only two or four straight-walled grooves per period on a planar hard surface. Different from the recently reported reflectors or splitters, our device can perfectly split an incident wave into different desired directions with arbitrary power flow partition. In addition, because ultrathin substructures with thin walls and narrow channels are avoided in our design procedure, the proposed beam splitter can be used for waves with much shorter wavelength compared to the previous suggested systems. The design is established by rigorous formulae developed under the framework of the grating theory and a genetic optimization algorithm. Numerical simulation and experimental evidence are provided to discuss the involved physical mechanism and to give the proof-of-concept for the proposed perfect anomalous acoustic splitter.


title: End-to-End Emotional Speech Synthesis Using Style Tokens and Semi-Supervised Training

abstract: This paper proposes an end-to-end emotional speech synthesis (ESS) method which adopts global style tokens (GSTs) for semi-supervised training. This model is built based on the GST-Tacotron framework. The style tokens are defined to present emotion categories. A cross entropy loss function between token weights and emotion labels is designed to obtain the interpretability of style tokens utilizing the small portion of training data with emotion labels. Emotion recognition experiments confirm that this method can achieve one-to-one correspondence between style tokens and emotion categories effectively. Objective and subjective evaluation results show that our model outperforms the conventional Tacotron model for ESS when only 5\% of training data has emotion labels. Its subjective performance is close to the Tacotron model trained using all emotion labels.


title: Lawn: an Unbound Low Latency Timer Data Structure for Large Scale, High Throughput Systems

abstract: As demand for Real-Time applications rises among the general public, the importance of enabling large-scale, unbound algorithms to solve conventional problems with low to no latency is critical for product viability. Timer algorithms are prevalent in the core mechanisms behind operating systems, network protocol implementation, stream processing, and several database capabilities. This paper presents a field-tested algorithm for low latency, unbound range timer structure, based upon the well excepted Timing Wheel algorithm. Using a set of queues hashed by TTL, the algorithm allows for a simpler implementation, minimal overhead no overflow and no performance degradation in comparison to the current state of the algorithms under typical use cases.


title: Assessing Post Deletion in Sina Weibo: Multi-modal Classification of Hot Topics

abstract: Widespread Chinese social media applications such as Weibo are widely known for monitoring and deleting posts to conform to Chinese government requirements. In this paper, we focus on analyzing a dataset of censored and uncensored posts in Weibo. Despite previous work that only considers text content of posts, we take a multi-modal approach that takes into account both text and image content. We categorize this dataset into 14 categories that have the potential to be censored on Weibo, and seek to quantify censorship by topic. Specifically, we investigate how different factors interact to affect censorship. We also investigate how consistently and how quickly different topics are censored. To this end, we have assembled an image dataset with 18,966 images, as well as a text dataset with 994 posts from 14 categories. We then utilized deep learning, CNN localization, and NLP techniques to analyze the target dataset and extract categories, for further analysis to better understand censorship mechanisms in Weibo. We found that sentiment is the only indicator of censorship that is consistent across the variety of topics we identified. Our finding matches with recently leaked logs from Sina Weibo. We also discovered that most categories like those related to anti-government actions (e.g. protest) or categories related to politicians (e.g. Xi Jinping) are often censored, whereas some categories such as crisis-related categories (e.g. rainstorm) are less frequently censored. We also found that censored posts across all categories are deleted in three hours on average.


title: Highly efficient acoustic refractive metasurfaces by harnessing near field coupling

abstract: Typical acoustic refractive metasurfaces governed by generalized Snell law require several types of subwavelength subunits to provide an extra phase gradient along the surface. This design strategy, however, has several kinds of drawback. For instance, the inevitable viscous loss brought out by the complex subwavelength subunits, and the negligence of the coupling between adjacent subunits which leads to low-efficiency in wavefront manipulation, especially for large angles. To overcome these limitations, we propose a new type of refractive metasurface composed of only one straight channel and several surface-etched grooves per period. By harnessing the nonlocal coupling between the channel/grooves, and the evanescent modes inside them, superiorly feasible acoustic transmission manipulations can be achieved. Nearly perfect acoustic bending with transmission efficiency up to 95% is demonstrated with theory and experiments for an extremely large angle of 81 degree. The reported results introduce a novel concept of acoustic metasurfaces and offer a real leap towards the development of high-efficient acoustic devices for wavefront manipulation.


title: Cross-correlation between the 21-cm signal and [OIII] emitters during early cosmic reionization

abstract: We study statistics of the 21-cm signal from the epoch of reionization. We propose to use [OIII] line emitting galaxies to cross-correlate with the 21-cm signal from $z = 7 - 10$. To this aim, we employ simulations of reionization obtained post-processing the high-resolution cosmological hydrodynamic simulation Massive Black-II with the three-dimensional (3D) radiative transfer code CRASH to follow the propagation of ionizing photons from a variety of sources. We show that, during the early phases of reionization, the 21-cm signal is positively correlated with the spatial distribution of the [OIII] emitters on large scales $(k < 1h~\rm cMpc^{-1})$. This positive correlation is generated by the temperature-galaxy correlation and it is a few times larger than when we assume that the heating is saturated. As the reionized regions expand, the correlation changes its sign to negative from $z = 10$ to 8. The signals at this epoch can be detected by combining the Square Kilometre Array (SKA) and a wide-field [OIII] emitter survey. We also calculate the cross-power spectrum with a 3D [OIII] intensity field, aiming at exploiting future intensity mapping observations. We conclude that high-redshift [OIII] line emitters can be used to probe the reionization process when the intergalactic medium is largely neutral.


title: Cross-correlated Contrast Source Inversion

abstract: In this paper, we improved the performance of the contrast source inversion (CSI) method by incorporating a so-called cross-correlated cost functional, which interrelates the state error and the data error in the measurement domain. The proposed method is referred to as the cross-correlated CSI. It enables better robustness and higher inversion accuracy than both the classical CSI and multiplicative regularized CSI (MR-CSI). In addition, we show how the gradient of the modified cost functional can be calculated without significantly increasing the computational burden. The advantages of the proposed algorithms are demonstrated using a 2-D benchmark problem excited by a transverse magnetic wave as well as a transverse electric wave, respectively, in comparison to classical CSI and MR-CSI.


title: The Syntax of the Accounting Language: A First Step

abstract: We review and interpret two basic propositions published by Ellerman (2014). The propositions address the algebraic structure of T accounts and double entry bookkeeping (DEB). The paper builds on this previous contribution with the view of reconciling the two, apparently dichotomous, perspectives of accounting measurement: the one that focuses preferably on the stock of wealth and to the one that focuses preferably on the flow of income. The paper claims that T-accounts and DEB have an underlying algebraic structure suitable for approaching measurement from either or both perspectives. Accountants preferences for stocks or flows can be framed in ways which are mutually consistent. The paper is a first step in addressing this consistency issue. It avoids the difficult mathematics of abstract algebra by applying the concept of syntax to accounting numbers such that the accounting procedure qualifies as a formal language with which accountants convey meaning.


title: $\Omega$-symmetric measures and related singular integrals

abstract: Let $\mathbb{S} \subset \mathbb{C}$ be the circle in the plane, and let $\Omega: \mathbb{S} \to \mathbb{S}$ be an odd bi-Lipschitz map with constant $1+\delta_\Omega$, where $\delta_\Omega>0$ is small. Assume also that $\Omega$ is twice continuously differentiable. Motivated by a question raised by Mattila and Preiss in [MP95], we prove the following: if a Radon measure $\mu$ has positive lower density and finte upper density almost everywhere, and the limit $$ \lim_{\epsilon \downarrow 0} \int_{\mathbb{C} \setminus B(x,\epsilon)} \frac{\Omega\left((x-y)/|x-y|\right)}{|x-y|} \, d\mu(y) $$ exists $\mu$-almost everywhere, then $\mu$ is $1$-rectifiable. To achieve this, we prove first that if an Ahlfors-David 1-regular measure $\mu$ is symmetric with respect to $\Omega$, that is, if $$ \int_{B(x,r)} |x-y|\Omega\left(\frac{x-y}{|x-y|}\right) \, d\mu(y) = 0 \mbox{ for all } x \in \mbox{spt}(\mu) \mbox{ and } r>0, $$ then $\mu$ is flat, or, in other words, there exists a constant $c>0$ and a line $L$ so that $\mu= c \mathcal{H}^{1}|_{L}$.


title: Conclusive nonlinear phase sensitivity limit for a Mach-Zehnder interferometer with single-mode non-vacuum inputs

abstract: Many works have stated that nonlinear interactions can improve phase sensitivity beyond the Heisenberg limit scaling of $1/N$ with $N$ being the mean photon number. This raises some open questions---among them the conclusive sensitivity limits with respect to single-mode inputs. Namely, when one of two inputs is vacuum, is there a shot-noise-style sensitivity bound on a nonlinear Mach-Zehnder interferometer? Within the reach of second-order nonlinear phase shifts, we make an attempt to provide an answer to this question. Based upon phase-averaging approach, this puzzle is partially resolved with careful calculations of the quantum Fisher information regarding three kinds of common inputs: Gaussian states, squeezed number states, and Schr\"odinger cat states. The results suggest that shot-noise-style sensitivity limit is no longer available, and the ideal candidate is squeezed vacuum.


title: Level inversion in kaonic nuclei and the high-density nuclear equation of state

abstract: It is very difficult for any nuclear model to pin down the saturation property and high-density equation of state (EOS) simultaneously because of high nonlinearity of the nuclear many-body problem. In this work, we propose, for the first time, to use the special property of light kaonic nuclei to characterize the relation between saturation property and high-density EOS. With a series of relativistic mean-field models, this special property is found to be the level inversion between orbitals $2S_{1/2}$ and $1D_{5/2}$ in light kaonic nuclei. This level inversion can serve as a theoretical laboratory to group the incompressibility at saturation density and the EOS at supra-normal densities simultaneously.


title: Piecewise polynomial approximation of probability density functions with application to uncertainty quantification for stochastic PDEs

abstract: The probability density function (PDF) associated with a given set of samples is approximated by a piecewise-linear polynomial constructed with respect to a binning of the sample space. The kernel functions are a compactly supported basis for the space of such polynomials, i.e. finite element hat functions, that are centered at the bin nodes rather than at the samples, as is the case for the standard kernel density estimation approach. This feature naturally provides an approximation that is scalable with respect to the sample size. On the other hand, unlike other strategies that use a finite element approach, the proposed approximation does not require the solution of a linear system. In addition, a simple rule that relates the bin size to the sample size eliminates the need for bandwidth selection procedures. The proposed density estimator has unitary integral, does not require a constraint to enforce positivity, and is consistent. The proposed approach is validated through numerical examples in which samples are drawn from known PDFs. The approach is also used to determine approximations of (unknown) PDFs associated with outputs of interest that depend on the solution of a stochastic partial differential equation.


title: Small projective spaces and Stillman uniformity for sheaves

abstract: We prove an analogue of Ananyan--Hochster's small subalgebra theorem in the context of sheaves on projective space, and deduce from this a version of Stillman's Conjecture for cohomology tables of sheaves. The main tools in the proof are Draisma's GL-noetherianity theorem and the BGG correspondence.


title: Discrete Symmetries and Nonlocal Reductions

abstract: We show that nonlocal reductions of systems of integrable nonlinear partial differential equations are the special discrete symmetry transformations.


title: Gaussian concentration bound and Ensemble equivalence in generic quantum many-body systems including long-range interaction

abstract: This work explores fundamental statistical and thermodynamic properties of short-and long-range-interacting systems. The purpose of this study is twofold. Firstly, we rigorously prove that the probability distribution of arbitrary few-body observables is restricted by a Gaussian concentration bound (or Chernoff--Hoeffding inequality) above some threshold temperature. This bound is then derived for arbitrary Gibbs states of systems that include long-range interactions Secondly, we establish a quantitative relationship between the concentration bound of the Gibbs state and the equivalence of canonical and micro-canonical ensembles. We then evaluate the difference in the averages of thermodynamic properties between the canonical and the micro-canonical ensembles. Under the assumption of the Gaussian concentration bound on the canonical ensemble, the difference between the ensemble descriptions is upper-bounded by $\left[n^{-1} \log (n^{3/2}\Delta^{-1})\right]^{1/2}$ with $n$ being the system size and $\Delta$ being the width of the energy shell of the micro-canonical ensemble This limit gives a non-trivial upper bound \textit{exponentially small energy width} with respect to the system size. By combining these two results, we prove the ensemble equivalence as well as the weak eigenstate thermalization in arbitrary long-range-interacting systems above a threshold temperature.

categories: cond-mat.stat-mech math-ph math.MP quant-ph


title: Making Smartphone Application Permissions Meaningful for the Average User

abstract: Smartphones hold important private information, yet users routinely expose this information to questionable applications written by developers they know nothing about. Users may be tempted to think of smartphones as old-style dumb phones, not as powerful network-connected computers, and this opens a gap between the permissions-based security paradigm (offered by platforms like Android) and what users expect. This makes it easy to fool users into installing applications that steal their information. Not surprisingly, Android is now a more favored target for hackers than Windows. We propose an approach for closing this gap, based on the observation that the current permissions system--rooted in good ol' UNIX-style thinking--is both too coarse and too fine grained, because it uses the wrong axes for defining the permissions space. We argue for replacing the paradigm in which "an app accesses device resources" (which is foreign to most non-geeks) with a paradigm in which "an app accesses user-tangible services." By using a simple piece of middleware, we can wrap this view of application control around today's permission system, and, by doing so, no conceptual refactoring of applications is required.


title: Well-posedness and regularity for a fractional tumor growth model

abstract: In this paper, we study a system of three evolutionary operator equations involving fractional powers of selfadjoint, monotone, unbounded, linear operators having compact resolvents. This system constitutes a generalization of a phase field system of Cahn-Hilliard type modelling tumor growth that has been proposed in Hawkins-Daarud et al. (Int. J. Numer. Math. Biomed. Eng. 28 (2012), 3-24) and investigated in recent papers co-authored by the present authors and E. Rocca. The model consists of a Cahn-Hilliard equation for the tumor cell fraction, coupled to a reaction-diffusion equation for a function S representing the nutrient-rich extracellular water volume fraction. Effects due to fluid motion are neglected. The generalization investigated in this paper is motivated by the possibility that the diffusional regimes governing the evolution of the different constituents of the model may be of different (e.g., fractional) type. Under rather general assumptions, well-posedness and regularity results are shown. In particular, by writing the equation governing the evolution of the chemical potential in the form of a general variational inequality, also singular or nonsmooth contributions of logarithmic or of double obstacle type to the energy density can be admitted.


title: A Linear Method for Shape Reconstruction based on the Generalized Multiple Measurement Vectors Model

abstract: In this paper, a novel linear method for shape reconstruction is proposed based on the generalized multiple measurement vectors (GMMV) model. Finite difference frequency domain (FDFD) is applied to discretized Maxwell's equations, and the contrast sources are solved iteratively by exploiting the joint sparsity as a regularized constraint. Cross validation (CV) technique is used to terminate the iterations, such that the required estimation of the noise level is circumvented. The validity is demonstrated with an excitation of transverse magnetic (TM) experimental data, and it is observed that, in the aspect of focusing performance, the GMMV-based linear method outperforms the extensively used linear sampling method (LSM).


title: Auxiliary Interference Speaker Loss for Target-Speaker Speech Recognition

abstract: In this paper, we propose a novel auxiliary loss function for target-speaker automatic speech recognition (ASR). Our method automatically extracts and transcribes target speaker's utterances from a monaural mixture of multiple speakers speech given a short sample of the target speaker. The proposed auxiliary loss function attempts to additionally maximize interference speaker ASR accuracy during training. This will regularize the network to achieve a better representation for speaker separation, thus achieving better accuracy on the target-speaker ASR. We evaluated our proposed method using two-speaker-mixed speech in various signal-to-interference-ratio conditions. We first built a strong target-speaker ASR baseline based on the state-of-the-art lattice-free maximum mutual information. This baseline achieved a word error rate (WER) of 18.06% on the test set while a normal ASR trained with clean data produced a completely corrupted result (WER of 84.71%). Then, our proposed loss further reduced the WER by 6.6% relative to this strong baseline, achieving a WER of 16.87%. In addition to the accuracy improvement, we also showed that the auxiliary output branch for the proposed loss can even be used for a secondary ASR for interference speakers' speech.


title: Security Rating Metrics for Distributed Wireless Systems

abstract: The paper examines quantitative assessment of wireless distribution system security, as well as an assessment of risks from attacks and security violations. Furthermore, it describes typical security breach and formal attack models and five methods for assessing security. The proposed normalized method for assessing the degree of security assurance operates with at least three characteristics, which allows comparatively analyze heterogeneous information systems. The improved calculating formulas have been proposed for two security assessment methods, and the elements of functional-cost analysis have been applied to calculate the degree of security. To check the results of the analysis, the coefficient of concordance was calculated, which gives opportunity to determine the quality of expert assessment. The simultaneous use of several models to describe attacks and the effectiveness of countering them allows us to create a comprehensive approach to countering modern security threats to information networks at the commercial enterprises and critical infrastructure facilities.


title: Men-in-the-Middle Attack Simulation on Low Energy Wireless Devices using Software Define Radio

abstract: The article presents a method of organizing men-in-the-middle attack and penetration test on Bluetooth Low Energy devices and ZigBee packets using software define radio with sniffing and spoofing packets, capture and analysis techniques on wireless waves with the focus on Bluetooth. The paper contains the analysis of the latest scientific work in this area, provides a comparative analysis of SDRs and the rationale for the choice of hardware, gives the sequence of actions for collecting wireless data packets and data collection from ZigBee and BLE devices, and analyzes ways to improve captured wireless packet analysis techniques. For the study collected experimental setup, the results of which are analyzed in real time. The collected wireless data packets are compared with those sent. The result of the experiment shows the weaknesses of local wireless networks.


title: N-Body Simulations of Gas-free Disc Galaxies with SMBH Seed in Binary Systems

abstract: We have shown the outcome of N-body simulations of the interactions of two disc galaxies without gas with the same mass. Both disc galaxies have halos of dark matter, central bulges and initial supermassive black hole (SMBH) seeds at their centers. The purpose of this work is to study the mass and dynamical evolution of the initial SMBH seed during a Hubble cosmological time. It is a complementation of our previous paper with different initial orbit conditions and by introducing the SMBH seed in the initial galaxy. The disc of the secondary galaxy has coplanar or polar orientation in relation to the disc of the primary galaxy and their initial orbit are eccentric and prograde. The primary and secondary galaxies have mass and size of Milky Way with an initial SMBH seed. We have found that the merger of the primary and secondary discs can result in a final normal disc or a final warped disc. After the fusion of discs, the final one is thicker and larger than the initial disc. The tidal effects are very important, modifying the evolution of the SMBH in the primary and secondary galaxy differently. The mass of the SMBH of the primary galaxy have increased by a factor ranging from 52 to 64 times the initial seed mass, depending on the experiment. However, the mass of the SMBH of the secondary galaxy have increased by a factor ranging from 6 to 33 times the initial SMBH seed mass, depending also on the experiment. Most of the accreted particles have come from the bulge and from the halo, depleting their particles. This could explain why the observations show that the SMBH with masses of approximately $10^6 M_\odot$ are found in many bulgeless galaxies. Only a small number of the accreted particles has come from the disc. In some cases of final merging stage of the two galaxies, the final SMBH of the secondary galaxy was {ejected out of the galaxy}.


title: New Explicit Lorentzian Einstein-Weyl Structures in 3-Dimensions

abstract: On a $3$D manifold, a Weyl geometry consists of pairs $(g, A) =$ (metric, $1$-form) modulo gauge $\widehat{g} = {\rm e}^{2\varphi} g$, $\widehat{A} = A + {\rm d}\varphi$. In 1943, Cartan showed that every solution to the Einstein-Weyl equations $R_{(\mu\nu)} - \frac{1}{3} R g_{\mu\nu} = 0$ comes from an appropriate $3$D leaf space quotient of a $7$D connection bundle associated with a 3$^{\rm rd}$ order ODE $y''' = H(x,y,y',y'')$ modulo point transformations, provided $2$ among $3$ primary point invariants vanish $$ \text{W\"unschmann}(H) \equiv 0\equiv \text{Cartan}(H). $$We find that point equivalence of a single PDE $z_y = F(x,y,z,z_x)$ with para-CR integrability $DF := F_x + z_x F_z \equiv 0$ leads to a completely similar $7$D Cartan bundle and connection. Then magically, the (complicated) equation $\text{W\"unschmann}(H) \equiv 0$ becomes $$0\equiv\text{Monge}(F):=9F_{pp}^2F_{ppppp}-45F_{pp}F_{ppp}F_{pppp}+40F_{ppp}^3,\qquad p:=z_x, $$ whose solutions are just conics in the $\{p, F\}$-plane. As an ansatz, we take $$F(x,y,z,p):= \frac{\alpha(y)(z-xp)^2+\beta(y)(z-xp)p+\gamma(y)(z-xp) +\delta(y)p^2+\varepsilon(y)p+\zeta(y)}{\lambda(y)(z-xp)+\mu(y) p+\nu(y)}, $$ with $9$ arbitrary functions $\alpha, \dots, \nu$ of $y$. This $F$ satisfies $DF \equiv 0 \equiv \text{Monge}(F)$, and we show that the condition $\text{Cartan}(H) \equiv 0 $ passes to a certain $\boldsymbol{K}(F) \equiv 0$ which holds for any choice of $\alpha(y), \dots, \nu(y)$. Descending to the leaf space quotient, we gain $\infty$-dimensional functionally parametrized and explicit families of Einstein-Weyl structures $\big[ (g, A) \big]$ in $3$D. These structures are nontrivial in the sense that ${\rm d}A \not\equiv 0$ and $\text{Cotton}([g]) \not \equiv 0$.


title: Automatic Hierarchical Classification of Kelps using Deep Residual Features

abstract: Across the globe, remote image data is rapidly being collected for the assessment of benthic communities from shallow to extremely deep waters on continental slopes to the abyssal seas. Exploiting this data is presently limited by the time it takes for experts to identify organisms found in these images. With this limitation in mind, a large effort has been made globally to introduce automation and machine learning algorithms to accelerate both classification and assessment of marine benthic biota. One major issue lies with organisms that move with swell and currents, like kelps. This paper presents an automatic hierarchical classification method (local binary classification as opposed to the conventional flat classification) to classify kelps in images collected by autonomous underwater vehicles. The proposed kelp classification approach exploits learned feature representations extracted from deep residual networks. We show that these generic features outperform the traditional off-the-shelf CNN features and the conventional hand-crafted features. Experiments also demonstrate that the hierarchical classification method outperforms the traditional parallel multi-class classifications by a significant margin (90.0% vs 57.6% and 77.2% vs 59.0%) on Benthoz15 and Rottnest datasets respectively. Furthermore, we compare different hierarchical classification approaches and experimentally show that the sibling hierarchical training approach outperforms the inclusive hierarchical approach by a significant margin. We also report an application of our proposed method to study the change in kelp cover over time for annually repeated AUV surveys.


title: Automatic Co-Registration of Aerial Imagery and Untextured Model Data Utilizing Average Shading Gradients

abstract: The comparison of current image data with existing 3D model data of a scene provides an efficient method to keep models up to date. In order to transfer information between 2D and 3D data, a preliminary co-registration is necessary. In this paper, we present a concept to automatically co-register aerial imagery and untextured 3D model data. To refine a given initial camera pose, our algorithm computes dense correspondence fields using SIFT flow between gradient representations of the model and camera image, from which 2D-3D correspondences are obtained. These correspondences are then used in an iterative optimization scheme to refine the initial camera pose by minimizing the reprojection error. Since it is assumed that the model does not contain texture information, our algorithm is built up on an existing method based on Average Shading Gradients (ASG) to generate gradient images based on raw geometry information only. We apply our algorithm for the co-registering of aerial photographs to an untextured, noisy mesh model. We have investigated different magnitudes of input error and show that the proposed approach can reduce the final reprojection error to a minimum of 1.27 plus-minus 0.54 pixels, which is less than 10 % of its initial value. Furthermore, our evaluation shows that our approach outperforms the accuracy of a standard Iterative Closest Point (ICP) implementation.


title: Branched splines

abstract: Spline functions have long been used in numerical solution of differential equations. Recently it revives as isogeometric analysis, which offers integration of finite element analysis and NURBS based CAD into a single unified process. Usually many NURBS pieces are needed to build geometrically continuous CAD models. In this paper, we introduce some splines defined on branched covering of sphere, torus or general domains of $R^2$, which are called branched splines in this paper. A single piece of such splines is enough to build some complex CAD models. Multiresolution analysis on surfaces of high genus built from such splines can be carried out naturally. CAD and FEA are integrated directly on such models. A theoretical framework is presented in this paper, together with some simple examples.


title: Christian Horrebow's Sunspot Observations -- I. Life and Published Writings

abstract: Between 1761 and 1776, Christian Horrebow made regular observations of sunspots from Rundetaarn in Copenhagen. Based on these observations he writes in 1775 that "it appears that after the course of a certain number of years, the appearance of the Sun repeats itself with respect to the number and size of the spots". Thus, Horrebow hypothesized the idea of a cyclic Sun several decades before Heinrich Schwabe discovered the solar cycle and estimated its period. This proves the ability of Horrebow as a sunspot observer. In this article, we present a general overview of the work of Christian Horrebow, including a brief biography and a complete bibliography. We also present a translation from Danish to English of his writings on sunspots in the Dansk Historisk Almanak. These writings include tables of daily sunspot measurements of which we discuss the completeness.


title: FatPaths: Routing in Supercomputers and Data Centers when Shortest Paths Fall Short

abstract: We introduce FatPaths: a simple, generic, and robust routing architecture that enables state-of-the-art low-diameter topologies such as Slim Fly to achieve unprecedented performance. FatPaths targets Ethernet stacks in both HPC supercomputers as well as cloud data centers and clusters. FatPaths exposes and exploits the rich ("fat") diversity of both minimal and non-minimal paths for high-performance multi-pathing. Moreover, FatPaths uses a redesigned "purified" transport layer that removes virtually all TCP performance issues (e.g., the slow start), and incorporates flowlet switching, a technique used to prevent packet reordering in TCP networks, to enable very simple and effective load balancing. Our design enables recent low-diameter topologies to outperform powerful Clos designs, achieving 15% higher net throughput at 2x lower latency for comparable cost. FatPaths will significantly accelerate Ethernet clusters that form more than 50% of the Top500 list and it may become a standard routing scheme for modern topologies.


title: Joint Multi-frame Detection and Segmentation for Multi-cell Tracking

abstract: Tracking living cells in video sequence is difficult, because of cell morphology and high similarities between cells. Tracking-by-detection methods are widely used in multi-cell tracking. We perform multi-cell tracking based on the cell centroid detection, and the performance of the detector has high impact on tracking performance. In this paper, UNet is utilized to extract inter-frame and intra-frame spatio-temporal information of cells. Detection performance of cells in mitotic phase is improved by multi-frame input. Good detection results facilitate multi-cell tracking. A mitosis detection algorithm is proposed to detect cell mitosis and the cell lineage is built up. Another UNet is utilized to acquire primary segmentation. Jointly using detection and primary segmentation, cells can be fine segmented in highly dense cell population. Experiments are conducted to evaluate the effectiveness of our method, and results show its state-of-the-art performance.


title: Spatial Transformer for 3D Point Clouds

abstract: Deep neural networks are widely used for understanding 3D point clouds. At each point convolution layer, features are computed from local neighborhoods of 3D points and combined for subsequent processing in order to extract semantic information. Existing methods adopt the same individual point neighborhoods throughout the network layers, defined by the same metric on the fixed input point coordinates. This common practice is easy to implement but not necessarily optimal. Ideally, local neighborhoods should be different at different layers, as more latent information is extracted at deeper layers. We propose a novel end-to-end approach to learn different non-rigid transformations of the input point cloud so that optimal local neighborhoods can be adopted at each layer. We propose both linear (affine) and non-linear (projective and deformable) spatial transformers for 3D point clouds. With spatial transformers on the ShapeNet part segmentation dataset, the network achieves higher accuracy for all categories, with 8\% gain on earphones and rockets in particular. Our method also outperforms the state-of-the-art on other point cloud tasks such as classification, detection, and semantic segmentation. Visualizations show that spatial transformers can learn features more efficiently by dynamically altering local neighborhoods according to the geometry and semantics of 3D shapes in spite of their within-category variations. Our code is publicly available at https://github.com/samaonline/spatial-transformer-for-3d-point-clouds.


title: European Option Pricing of electricity under exponential functional of L\'evy processes with Price-Cap principle

abstract: We propose a new model for electricity pricing based on the price cap principle. The particularity of the model is that the asset price is an exponential functional of a jump L\'evy process. This model can capture both mean reversion and jumps which are observed in electricity market. It is shown that the value of an European option of this asset is the unique viscosity solution of a partial integro-differential equation (PIDE). A numerical approximation of this solution by the finite differences method is provided. The consistency, stability and convergence results of the scheme are given. Numerical simulations are performed under a smooth initial condition.


title: Dynamics of reverse annealing for the fully-connected $p$-spin model

abstract: Reverse annealing is a relatively new variant of quantum annealing, in which one starts from a classical state and increases and then decreases the amplitude of the transverse field, in the hope of finding a better classical state than the initial state for a given optimization problem. We numerically study the unitary quantum dynamics of reverse annealing for the mean-field-type $p$-spin model and show that the results are consistent with the predictions of equilibrium statistical mechanics. In particular, we corroborate the equilibrium analysis prediction that reverse annealing provides an exponential speedup over conventional quantum annealing in terms of solving the $p$-spin model. This lends support to the expectation that equilibrium analyses are effective at revealing essential aspects of the dynamics of quantum annealing. We also compare the results of quantum dynamics with the corresponding classical dynamics, to reveal their similarities and differences. We distinguish between two reverse annealing protocols we call adiabatic and iterated reverse annealing. We further show that iterated reverse annealing, as has been realized in the D-Wave device, is ineffective in the case of the $p$-spin model, but note that a recently-introduced protocol ("$h$-gain"), which implements adiabatic reverse annealing, may lead to improved performance.


title: Multi-scale Mining of Kinematic Distributions with Wavelets

abstract: Typical LHC analyses search for local features in kinematic distributions. Assumptions about anomalous patterns limit them to a relatively narrow subset of possible signals. Wavelets extract information from an entire distribution and decompose it at all scales, simultaneously searching for features over a wide range of scales. We propose a systematic wavelet analysis and show how bumps, bump-dip combinations, and oscillatory patterns are extracted. Our kinematic wavelet analysis kit KWAK provides a publicly available framework to analyze and visualize general distributions.


title: On the performance of residual block design alternatives in convolutional neural networks for end-to-end audio classification

abstract: Residual learning is a recently proposed learning framework to facilitate the training of very deep neural networks. Residual blocks or units are made of a set of stacked layers, where the inputs are added back to their outputs with the aim of creating identity mappings. In practice, such identity mappings are accomplished by means of the so-called skip or residual connections. However, multiple implementation alternatives arise with respect to where such skip connections are applied within the set of stacked layers that make up a residual block. While ResNet architectures for image classification using convolutional neural networks (CNNs) have been widely discussed in the literature, few works have adopted ResNet architectures so far for 1D audio classification tasks. Thus, the suitability of different residual block designs for raw audio classification is partly unknown. The purpose of this paper is to analyze and discuss the performance of several residual block implementations within a state-of-the-art CNN-based architecture for end-to-end audio classification using raw audio waveforms. For comparison purposes, we analyze as well the performance of the residual blocks under a similar 2D architecture using a conventional time-frequency audio represen-tation as input. The results show that the achieved accuracy is considerably dependent, not only on the specific residual block implementation, but also on the selected input normalization.


title: On a class of reaction-diffusion equations with aggregation

abstract: In this paper, global-in-time existence and blow up results are shown for a reaction-diffusion equation appearing in the theory of aggregation phenomena (including chemotaxis). Properties of the corresponding steady-state problem are also presented. Moreover, the stability around constant equilibria and the non-existence of non-constant solutions are studied in certain cases.


title: Privacy-Preserving Blockchain-Based Federated Learning for IoT Devices

abstract: Home appliance manufacturers strive to obtain feedback from users to improve their products and services to build a smart home system. To help manufacturers develop a smart home system, we design a federated learning (FL) system leveraging the reputation mechanism to assist home appliance manufacturers to train a machine learning model based on customers' data. Then, manufacturers can predict customers' requirements and consumption behaviors in the future. The working flow of the system includes two stages: in the first stage, customers train the initial model provided by the manufacturer using both the mobile phone and the mobile edge computing (MEC) server. Customers collect data from various home appliances using phones, and then they download and train the initial model with their local data. After deriving local models, customers sign on their models and send them to the blockchain. In case customers or manufacturers are malicious, we use the blockchain to replace the centralized aggregator in the traditional FL system. Since records on the blockchain are untampered, malicious customers or manufacturers' activities are traceable. In the second stage, manufacturers select customers or organizations as miners for calculating the averaged model using received models from customers. By the end of the crowdsourcing task, one of the miners, who is selected as the temporary leader, uploads the model to the blockchain. To protect customers' privacy and improve the test accuracy, we enforce differential privacy on the extracted features and propose a new normalization technique. We experimentally demonstrate that our normalization technique outperforms batch normalization when features are under differential privacy protection. In addition, to attract more customers to participate in the crowdsourcing FL task, we design an incentive mechanism to award participants.


title: Speed limit of quantum dynamics near the event horizon of black holes

abstract: Quantum mechanics imposes a fundamental bound on the minimum time required for the quantum systems to evolve between two states of interest. This bound introduces a limit on the speed of the dynamical evolution of the systems, known as the quantum speed limit. We show that black holes can drastically affect the speed limit of a two-level fermionic quantum system subjected to an open quantum dynamics. As we demonstrate, the quantum speed limit can enhance at the vicinity of a black hole's event horizon in the Schwarzschild spacetime.


title: Christian Horrebow's Sunspot Observations -- II. Construction of a Record of Sunspot Positions

abstract: The number of spots on the surface of the Sun is one of the best tracers of solar variability we have. The sunspot number is not only known to change in phase with the 11-year solar cycles, but also to show variability on longer time scales. It is, however, not only the sunspot number that changes in connection with solar variability. The location of the spots on the solar surface is also known to change in phase with the 11-year solar cycle. This has traditionally been visualised in the so-called butterfly diagram, but this is only well constrained from the beginning of the 19th century. This is unfortunate, as knowledge about the butterfly diagram could aid our understanding of the variability and the Sun-Earth connection. As part of a larger review of the work done on sunspots by the Danish astronomer Christian Horrebow, we here present a reanalysis of Christian Horrebow's notebooks covering the years 1761 and 1764 - 1777. These notebooks have been analysed in at least three earlier studies by Thiele (Astron. Nachr. 50, 257, 1859), d'Arrest (published in Wolf, Astron. Mitt. Eidgenoss. Sternwarte Zur. 4, 77, 1873) and Hoyt and Schatten (Solar Phys. 160, 387, 1995). In this article, we construct a complete record of sunspot positions covering the years 1761 and 1764 - 1777. The resulting butterfly diagram shows the characteristic structure known from observations in the 19th and 20th century. We do see some indications of equatorial sunspots in the observations we have from Cycle 1. However, in Cycle 2, which has much better coverage, we do not see such indications.


title: Visual Anomaly Detection in Event Sequence Data

abstract: Anomaly detection is a common analytical task that aims to identify rare cases that differ from the typical cases that make up the majority of a dataset. When applied to the analysis of event sequence data, the task of anomaly detection can be complex because the sequential and temporal nature of such data results in diverse definitions and flexible forms of anomalies. This, in turn, increases the difficulty in interpreting detected anomalies. In this paper, we propose an unsupervised anomaly detection algorithm based on Variational AutoEncoders (VAE) to estimate underlying normal progressions for each given sequence represented as occurrence probabilities of events along the sequence progression. Events in violation of their occurrence probability are identified as abnormal. We also introduce a visualization system, EventThread3, to support interactive exploration and interpretations of anomalies within the context of normal sequence progressions in the dataset through comprehensive one-to-many sequence comparison. Finally, we quantitatively evaluate the performance of our anomaly detection algorithm and demonstrate the effectiveness of our system through a case study.


title: Task-Driven Common Representation Learning via Bridge Neural Network

abstract: This paper introduces a novel deep learning based method, named bridge neural network (BNN) to dig the potential relationship between two given data sources task by task. The proposed approach employs two convolutional neural networks that project the two data sources into a feature space to learn the desired common representation required by the specific task. The training objective with artificial negative samples is introduced with the ability of mini-batch training and it's asymptotically equivalent to maximizing the total correlation of the two data sources, which is verified by the theoretical analysis. The experiments on the tasks, including pair matching, canonical correlation analysis, transfer learning, and reconstruction demonstrate the state-of-the-art performance of BNN, which may provide new insights into the aspect of common representation learning.


title: Structural and electronic properties of the incommensurate host-guest Bi-III phase

abstract: At high pressure, bismuth acquires a complex incommensurate host-guest structure, only recently discovered. Characterizing the structure and properties of this incommensurate phase from first principles is challenging owing to its non-periodic nature. In this study we use large scale DFT calculations to model commensurate approximants to the Bi-III phase, and in particular to describe the atomic modulations with respect to their ideal positions, shown here to strongly affect the electronic structure of the lattice and its stability. The equation of state and range of stability of Bi-III are reproduced in excellent agreement with experiment using a fully relativistic model. We demonstrate the importance of employing large unit-cells for the accurate description of the geometric and electronic configuration of Bi-III. In contrast, accurate description of the equation of state of bismuth is found to be primarily sensitive to the choice of pseudopotential and exchange-correlation function, while almost completely insensitive to the commensurate approximation.


title: Nonfactorizable QCD Effects in Higgs Boson Production via Vector Boson Fusion

abstract: We discuss nonfactorizable QCD corrections to Higgs boson production in vector boson fusion at the Large Hadron Collider. We point out that these corrections can be computed in the eikonal approximation retaining all the terms that are not suppressed by the ratio of the transverse momenta of the tagging jets to the total center-of-mass energy. Our analysis shows that in certain kinematic distributions the nonfactorizable corrections can be as large as a percent making them quite comparable to their factorizable counter-parts.


title: A Linear Model for Microwave Imaging of Highly Conductive Scatterers

abstract: In this paper, a linear model based on multiple measurement vectors model is proposed to formulate the inverse scattering problem of highly conductive objects at one single frequency. Considering the induced currents which are mostly distributed on the boundaries of the scatterers, joint sparse structure is enforced by a sum-of-norm regularization. Since no \textit{a priori} information is required and no approximation of the scattering model has been made, the proposed method is versatile. Imaging results with transverse magnetic and transverse electric polarized synthetic data and Fresnel data demonstrate its higher resolving ability than both linear sampling method and its improved version with higher, but acceptable, computational complexity.


title: Frequency response of time-delay interferometry for space-based gravitational wave antennas

abstract: Space-based gravitational wave detectors cannot keep rigid structures and precise arm length equality, so the precise equality of detector arms which is required in a ground-based interferometer to cancel the overwhelming laser noise is impossible. The time-delay interferometry method is applied to unequal arm lengths to cancel the laser frequency noise. We give analytical formulas of the averaged response functions for tensor, vector, breathing and longitudinal polarizations in different TDI combinations, and obtain their asymptotic behaviors. At low frequencies, $f\ll f_*$, the averaged response functions of all TDI combinations increase as $f^2$ for all six polarizations. The one exception is the averaged response functions of $\zeta$ for all six polarizations increase as $f^4$ in the equilateral-triangle case. At high frequencies, $f\gg f_*$, the averaged response functions of all TDI combinations for the tensor and breathing modes fall off as $1/f^2$, the averaged response functions of all TDI combinations for the vector mode fall off as $\ln(f)/f^2$ , and the averaged response functions of all TDI combinations for the longitudinal mode fall as $1/f$. We also give LISA and TianQin sensitivity curves in different TDI combinations for tensor, vector, breathing and longitudinal polarizations.


title: Quantum ammeter

abstract: We present the theoretical model of the "quantum ammeter", a device that is able to measure the full counting statistics of an electron current at quantum time scales. It consists of an Ohmic contact, perfectly coupled to chiral quantum Hall channels, and of a quantum dot attached to one of the outgoing channels. At energies small compared to its charging energy, the Ohmic contact fractionalizes each incoming electron and redistributes it between outgoing channels. By monitoring the resonant tunneling current through the quantum dot, one gets an access to the moment generator of the current in one of the incoming channels at time scales comparable to its correlation time.


title: Double parton distributions for a positronium-like bound state using light-front wave functions

abstract: We investigate the double parton distributions (DPDs) for a positronium-like bound state using light-front QED. We incorporate the higher Fock three particle component of the state, that includes a photon. We obtain the overlap representation of the DPDs in terms of the three-particle light-front wave functions (LFWFs). Our calculation explores the correlations between the momentum fractions of the particles probed and the transverse distance between them, without any assumption of factorization between them. We also investigate the behavior of the DPDs near the kinematical boundary when the sum of the momentum fractions is close to one.


title: Witnessing incompatibility of quantum channels

abstract: We introduce the notion of incompatibility witness for quantum channels, defined as an affine functional that is non-negative on all pairs of compatible channels and strictly negative on some incompatible pair. This notion extends the recent definition of incompatibility witnesses for quantum measurements. We utilize the general framework of channels acting on arbitrary finite dimensional von Neumann algebras, thus allowing us to investigate incompatibility witnesses on measurement-measurement, measurement-channel and channel-channel pairs. We prove that any incompatibility witness can be implemented as a state discrimination task in which some intermediate classical information is obtained before completing the task. This implies that any incompatible pair of channels gives an advantage over compatible pairs in some such state discrimination task.


title: Suppressed Auger scattering and tunable light emission of Landau-quantized massless Kane electrons

abstract: The Landau level laser has been proposed a long time ago as a unique source of monochromatic radiation, widely tunable in the THz and infrared spectral ranges using an externally applied magnetic field. In spite of decades of efforts, this appealing concept never resulted in the design of a reliable device. This is due to efficient Auger scattering of Landau-quantized electrons, which is an intrinsic non-radiative recombination channel that eventually gains over cyclotron emission in all materials studied so far: in conventional semiconductors with parabolic bands, but also in graphene with massless electrons. The Auger processes are favored in these systems by Landau levels (or their subsets) equally spaced in energy. Here we show that this scheme does not apply to massless Kane electrons in gapless HgCdTe alloy, in which undesirable Auger scattering is strongly suppressed and the sizeable cyclotron emission observed, for the first time in the case of massless particles. The gapless HgCdTe thus appears as a material of choice for future technology of Landau level lasers.


title: Improved H\"older regularity for strongly elliptic PDEs

abstract: We establish surprising improved Schauder regularity properties for solutions to the Leray-Lions divergence type equation in the plane. The results are achieved by studying the nonlinear Beltrami equation and making use of special new relations between these two equations. In particular, we show that solutions to an autonomous Beltrami equation enjoy a quantitative improved degree of H\"older regularity, higher than what is given by the classical exponent $1/K$.


title: Leveraging Text Repetitions and Denoising Autoencoders in OCR Post-correction

abstract: A common approach for improving OCR quality is a post-processing step based on models correcting misdetected characters and tokens. These models are typically trained on aligned pairs of OCR read text and their manually corrected counterparts. In this paper we show that the requirement of manually corrected training data can be alleviated by estimating the OCR errors from repeating text spans found in large OCR read text corpora and generating synthetic training examples following this error distribution. We use the generated data for training a character-level neural seq2seq model and evaluate the performance of the suggested model on a manually corrected corpus of Finnish newspapers mostly from the 19th century. The results show that a clear improvement over the underlying OCR system as well as previously suggested models utilizing uniformly generated noise can be achieved.


title: Prediction Poisoning: Towards Defenses Against DNN Model Stealing Attacks

abstract: High-performance Deep Neural Networks (DNNs) are increasingly deployed in many real-world applications e.g., cloud prediction APIs. Recent advances in model functionality stealing attacks via black-box access (i.e., inputs in, predictions out) threaten the business model of such applications, which require a lot of time, money, and effort to develop. Existing defenses take a passive role against stealing attacks, such as by truncating predicted information. We find such passive defenses ineffective against DNN stealing attacks. In this paper, we propose the first defense which actively perturbs predictions targeted at poisoning the training objective of the attacker. We find our defense effective across a wide range of challenging datasets and DNN model stealing attacks, and additionally outperforms existing defenses. Our defense is the first that can withstand highly accurate model stealing attacks for tens of thousands of queries, amplifying the attacker's error rate up to a factor of 85$\times$ with minimal impact on the utility for benign users.


title: Probabilistic Two-Ray Model for Air-to-Air Channel in Built-Up Areas

abstract: In this paper, we present a probabilistic two-ray (PTR) path loss model for air-to-air (AA) propagation channel in built-up areas. Based on the statistical model of city deployment, the PTR path loss model can be applied to suburban, urban, dense urban, and high-rise urban. The path loss is optimally fitted as the Weibull distribution and its fluctuation is fitted as the Normal distribution in ray-tracing simulations. The good agreements between our model and ray tracing indicate the proposed model can provide a useful tool for accurate and quick prediction for aerial platforms. As an extended research of PTR model, we extract the shadowing factor by numerous simulations and propose the altitude-dependent shadowing model. The result shows that the proposed shadowing model has very good consistent with the measurement-based model, which indicates that our research performs well in the extensibility and generality.


title: Creating A Neural Pedagogical Agent by Jointly Learning to Review and Assess

abstract: Machine learning plays an increasing role in intelligent tutoring systems as both the amount of data available and specialization among students grow. Nowadays, these systems are frequently deployed on mobile applications. Users on such mobile education platforms are dynamic, frequently being added, accessing the application with varying levels of focus, and changing while using the service. The education material itself, on the other hand, is often static and is an exhaustible resource whose use in tasks such as problem recommendation must be optimized. The ability to update user models with respect to educational material in real-time is thus essential; however, existing approaches require time-consuming re-training of user features whenever new data is added. In this paper, we introduce a neural pedagogical agent for real-time user modeling in the task of predicting user response correctness, a central task for mobile education applications. Our model, inspired by work in natural language processing on sequence modeling and machine translation, updates user features in real-time via bidirectional recurrent neural networks with an attention mechanism over embedded question-response pairs. We experiment on the mobile education application SantaTOEIC, which has 559k users, 66M response data points as well as a set of 10k study problems each expert-annotated with topic tags and gathered since 2016. Our model outperforms existing approaches over several metrics in predicting user response correctness, notably out-performing other methods on new users without large question-response histories. Additionally, our attention mechanism and annotated tag set allow us to create an interpretable education platform, with a smart review system that addresses the aforementioned issue of varied user attention and problem exhaustion.


title: A note on flatness of some fiber type contractions

abstract: We discuss the flatness property of some fiber type contractions of complex smooth projective varieties of arbitrary dimensions. We relate the flatness of some morphisms having one-dimensional fibers with their conic bundles structures, also in the general case in which some mild singularities of the varieties are admitted.


title: Coarse graining hadronic scattering

abstract: We show that it makes sense to coarse grain hadronic interactions such as $\pi\pi$ and $\pi N$ reactions following previous work on NN scattering. Moreover, if the interaction is taken to be given by chiral dynamics at long distances above a given value $r > r_c$ larger than the elementary radii of the interaction hadrons the unknown short distance region $r< r_c$ is characterized by a {\it finite} number of fitting parameters. This number of independent parameters needed for a presumably complete description of scattering data for a CM energy below $\sqrt{s}$ has been found to be given by $N_{\rm Par} = N_S \times N_I \times (p r_c )^2 /2 $ with $N_S$ and $N_I$ the number of spin and isospin channels, and $p$ the CM momentum respectively. Therefore, for an experiment (or sets of experiments) with a total number of data $N_{\rm Dat}$ the number of degrees of freedom involved in a $\chi^2$-fit is given by $\nu = N_{\rm Dat}-N_{\rm Par}$ and confidence levels can be obtained accordingly by standard means. Namely a $1 \sigma$ confidence level corresponds to $\chi_{\rm min}^2/\nu \in (1- \sqrt{2/\nu},1+\sqrt{2/\nu})$. We discuss the approach for $\pi\pi$ and $\pi N$ with an eye put on a data selection program and the eventual validation of chiral symmetry.


title: Singularities of rational inner functions in higher dimensions

abstract: We study the boundary behavior of rational inner functions (RIFs) in dimensions three and higher from both analytic and geometric viewpoints. On the analytic side, we use the critical integrability of the derivative of a rational inner function of several variables to quantify the behavior of a RIF near its singularities, and on the geometric side we show that the unimodular level sets of a RIF convey information about its set of singularities. We then specialize to three-variable degree $(m,n,1)$ RIFs and conduct a detailed study of their derivative integrability, zero set and unimodular level set behavior, and non-tangential boundary values. Our results, coupled with constructions of non-trivial RIF examples, demonstrate that much of the nice behavior seen in the two-variable case is lost in higher dimensions.


title: Possible glass-like random singlet magnetic state in 1T-TaS2

abstract: Two-dimensional layered transition-metal-dichalcogenide compound 1T-TaS2 shows the rare coexistence of charge density wave (CDW) and electron correlation driven Mott transition. In addition, atomic-cluster spins on the triangular lattice of the CDW state of 1T-TaS2 give rise to the possibility of the exotic spin-singlet state in which quantum fluctuations of spins are strong enough to prevent any long range magnetic ordering down to absolute zero ( 0 K). We present here the evidences of a glass-like random singlet magnetic state in 1T-TaS2 at low temperatures through a study of temperature and time dependence of magnetization. Comparing the experimental results with a representative canonical spin-glass system Au(1.8%Mn), we show that this glass-like state is distinctly different from the well established canonical spin-glass state.


title: First search for 2$\varepsilon$ and $\varepsilon\beta^+$ processes in $^{168}$Yb

abstract: The double-electron capture and the electron capture with positron emission in $^{168}$Yb have been investigated for the first time at the STELLA facility of the Gran Sasso underground laboratory (Italy) measuring 371 g of highly purified ytterbium oxide placed on the end-cap of a 465 cm$^3$ ultra-low-background high purity Germanium detector (HPGe). No gamma associated to double beta processes in $^{168}$Yb have been observed after 2074 h of data taking. This has allowed setting the half-life limits on the level of $\lim T_{1/2}\sim$ $10^{14}-10^{18}$ yr at 90\% C.L. Particularly, a lower half-life limit on a possible resonant neutrinoless double-electron capture in $^{168}$Yb to the $(2)^-$ 1403.7 keV excited state of $^{168}$Er is set as $T_{1/2}\geq1.9\times 10^{18}$ yr at 90\% C.L. Half-life limits $T_{1/2}^{2\nu(0\nu)}\geq 4.5(4.3)\times10^{16}$ yr were set on the $2\nu(0\nu)2\beta^-$ decay of $^{176}$Yb to the $2^+$ 84.3 keV first excited level of $^{176}$Hf.


title: Smart metastructure method for increasing TC of Bi(Pb)SrCaCuO high-temperature superconductors

abstract: Improving the critical transition temperature (TC) of Bi(Pb)SrCaCuO (B(P)SCCO) high-temperature superconductors is important, however, considerable challenges exist. In this study, on the basis of the metamaterial structure and the idea that the injecting energy will promote the formation of Cooper pairs, a smart meta-superconductor B(P)SCCO consisting of B(P)SCCO microparticles and Y2O3:Eu3++Ag or Y2O3:Eu3+ luminophor was designed. In the applied electric field, the Y2O3:Eu3++Ag or Y2O3:Eu3+ luminophor generates an electroluminescence (EL), thereby promoting the TC via EL energy injection. A series of Y2O3:Eu3++Ag topological luminophor-doped B(P)SCCO samples was prepared. Results showed that Y2O3:Eu3++Ag was dispersed around B(P)SCCO particles, forming a metastructure. Accordingly, the onset transition temperature (T_(C,on)) and zero resistance transition temperature (T_(C,0)) of B(P)SCCO increased. Meanwhile, the B(P)SCCO sample doped with 0.2 wt% Y2O3 or Y2O3:Sm3+ nonluminous inhomogeneous phase was also prepared to further prove the influence of EL on the T_C rather than the rare earth effect. Results indicated that the TC of the Y2O3 or Y2O3:Sm3+ doping sample decreased. However, the TC of the 0.2 wt% Y2O3:Eu3++Ag or Y2O3:Eu3+ luminophor-doped sample improved. This outcome further demonstrated that the smart metastructure method can improve the TC of B(P)SCCO.


title: On the correctness of Egalitarian Paxos

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


title: Towards Empathic Deep Q-Learning

abstract: As reinforcement learning (RL) scales to solve increasingly complex tasks, interest continues to grow in the fields of AI safety and machine ethics. As a contribution to these fields, this paper introduces an extension to Deep Q-Networks (DQNs), called Empathic DQN, that is loosely inspired both by empathy and the golden rule ("Do unto others as you would have them do unto you"). Empathic DQN aims to help mitigate negative side effects to other agents resulting from myopic goal-directed behavior. We assume a setting where a learning agent coexists with other independent agents (who receive unknown rewards), where some types of reward (e.g. negative rewards from physical harm) may generalize across agents. Empathic DQN combines the typical (self-centered) value with the estimated value of other agents, by imagining (by its own standards) the value of it being in the other's situation (by considering constructed states where both agents are swapped). Proof-of-concept results in two gridworld environments highlight the approach's potential to decrease collateral harms. While extending Empathic DQN to complex environments is non-trivial, we believe that this first step highlights the potential of bridge-work between machine ethics and RL to contribute useful priors for norm-abiding RL agents.


title: Thermal conductance of interleaving fins

abstract: Interleaving fins can significantly increase the heat transfer by increasing the effective area per unit base area. The fins are separated uniformly by a gap, which is filled with a flow medium to control the heat flux. The heat flux between the plates depends strongly on the thermal conductivity of the fin material and the medium between them as well as the dimensions. In earlier studies empirical a fitting method is used to determine the total effectiveness of the fins. However, it required complete characterization of the fins for each new set of operating conditions. In this paper, a simplified analytical model, but still preserving the main physical traits of the problem is developed. This model reveals the dimensionless parameter group containing both material properties and the fin geometry that govern the heat transfer. Rigorously testing of the model using a numerical finite element model shows an accuracy within 2% over a large parameter space, varying both dimensions and material properties. Lastly, this model is put to test with previously measured experimental data and a good agreement is obtained.


title: Control variate selection for Monte Carlo integration

abstract: Monte Carlo integration with variance reduction by means of control variates can be implemented by the ordinary least squares estimator for the intercept in a multiple linear regression model with the integrand as response and the control variates as covariates. Even without special knowledge on the integrand, significant efficiency gains can be obtained if the control variate space is sufficiently large. Incorporating a large number of control variates in the ordinary least squares procedure may however result in (i) a certain instability of the ordinary least squares estimator and (ii) a possibly prohibitive computation time. Regularizing the ordinary least squares estimator by preselecting appropriate control variates via the Lasso turns out to increase the accuracy without additional computational cost. The findings in the numerical experiment are confirmed by concentration inequalities for the integration error.


title: Clustering piecewise stationary processes

abstract: The problem of time-series clustering is considered in the case where each data-point is a sample generated by a piecewise stationary ergodic process. Stationary processes are perhaps the most general class of processes considered in non-parametric statistics and allow for arbitrary long-range dependence between variables. Piecewise stationary processes studied here for the first time in the context of clustering, relax the last remaining assumption in this model: stationarity. A natural formulation is proposed for this problem and a notion of consistency is introduced which requires the samples to be placed in the same cluster if and only if the piecewise stationary distributions that generate them have the same set of stationary distributions. Simple, computationally efficient algorithms are proposed and are shown to be consistent without any additional assumptions beyond piecewise stationarity.


title: Challenges for Security Assessment of Enterprises in the IoT Era

abstract: For years, attack graphs have been an important tool for security assessment of enterprise networks, but IoT devices, a new player in the IT world, might threat the reliability of this tool. In this paper, we review the challenges that must be addressed when using attack graphs to model and analyze enterprise networks that include IoT devices. In addition, we propose novel ideas and countermeasures aimed at addressing these challenges.


title: A crossinggram for random fields on lattices

abstract: The modeling of risk situations that occur in a space-time framework can be done using max-stable random fields on lattices. Although the summary coefficients for the spatial and temporal behaviour do not characterize the finite-dimensional distributions of the random field, they have the advantage of being immediate to interpret and easier to estimate. The coefficients that we propose give us information about the tendency of a random field for local oscillations of its values in relation to real valued high levels. It is not the magnitude of the oscillations that is being evaluated, but rather the greater or lesser number of oscillations, that is, the tendency of the trajectories to oscillate. We can observe surface trajectories more smooth over a region according to higher crossinggram value. It takes value in $[0,1]$ and increases with the concordance of the variables of the random field.


title: Interpretable Question Answering on Knowledge Bases and Text

abstract: Interpretability of machine learning (ML) models becomes more relevant with their increasing adoption. In this work, we address the interpretability of ML based question answering (QA) models on a combination of knowledge bases (KB) and text documents. We adapt post hoc explanation methods such as LIME and input perturbation (IP) and compare them with the self-explanatory attention mechanism of the model. For this purpose, we propose an automatic evaluation paradigm for explanation methods in the context of QA. We also conduct a study with human annotators to evaluate whether explanations help them identify better QA models. Our results suggest that IP provides better explanations than LIME or attention, according to both automatic and human evaluation. We obtain the same ranking of methods in both experiments, which supports the validity of our automatic evaluation paradigm.


title: FA-Harris: A Fast and Asynchronous Corner Detector for Event Cameras

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


title: Global rigidity of 2-dimensional linearly constrained frameworks

abstract: A linearly constrained framework in $\mathbb{R}^d$ is a point configuration together with a system of constraints which fixes the distances between some pairs of points and additionally restricts some of the points to lie in given affine subspaces. It is globally rigid if the configuration is uniquely defined by the constraint system, and is rigid if it is uniquely defined within some small open neighbourhood. Streinu and Theran characterised generic rigidity of linearly constrained frameworks in $\mathbb{R}^2$ in 2010. We obtain an analagous characterisation for generic global rigidity in $\mathbb{R}^2$. More precisely we show that a generic linearly constrained framework in $\mathbb{R}^2$ is globally rigid if and only if it is redundantly rigid and `balanced'. For generic frameworks which are not balanced, we determine the precise number of solutions to the constraint system whenever the underlying rigidity matroid of the given framework is connected. We also obtain a stress matrix sufficient condition and a Hendrickson type necessary condition for a generic linearly constrained framework to be globally rigid in $\mathbb{R}^d$.


title: A low-energy compact Shanghai-Wuhan electron beam ion trap for extraction of highly charged ions

abstract: A low-energy, compact and superconducting electron beam ion trap (the Shanghai-Wuhan EBIT or SW-EBIT) for extraction of highly charged ions is presented. The magnetic field in the central drift tube of the SW-EBIT is approximately 0.21 T produced by a pair of high-temperature superconducting coils. The electron-beam energy of the SW-EBIT is in the range of 30-4000 eV, and the maximum electron-beam current is up to 9 mA. Acting as a source of highly charged ions, the ion-beam optics for extraction is integrated, including an ion extractor and an einzel lens. A Wien filter is then used to measure the charge-state distribution of the extracted ions. In this work, the tungsten ions below the charge state of 15 have been produced, extracted, and analyzed. The charge-state distributions and spectra in the range of 530-580 nm of tungsten ions have been measured simultaneously with the electron-beam energy of 279 eV and 300 eV, which preliminarily indicates that the 549.9 nm line comes from $W^{14+}$.


title: A wrinkle in time: A case study in DNS poisoning

abstract: The Domain Name System (DNS) provides a translation between readable domain names and IP addresses. The DNS is a key infrastructure component of the Internet and a prime target for a variety of attacks. One of the most significant threat to the DNS's wellbeing is a DNS poisoning attack, in which the DNS responses are maliciously replaced, or poisoned, by an attacker. To identify this kind of attack, we start by an analysis of different kinds of response times. We present an analysis of typical and atypical response times, while differentiating between the different levels of DNS servers' response times, from root servers down to internal caching servers. We successfully identify empirical DNS poisoning attacks based on a novel method for DNS response timing analysis. We then present a system we developed to validate our technique that does not require any changes to the DNS protocol or any existing network equipment. Our validation system tested data from different architectures including LAN and cloud environments and real data from an Internet Service Provider (ISP). Our method and system differ from most other DNS poisoning detection methods and achieved high detection rates exceeding 99%. These findings suggest that when used in conjunction with other methods, they can considerably enhance the accuracy of these methods.


title: Entanglement Certification $-$ From Theory to Experiment

abstract: Entanglement is an important resource that allows quantum technologies to go beyond the classically possible. There are many ways quantum systems can be entangled, ranging from the archetypal two-qubit case to more exotic scenarios of entanglement in high dimensions or between many parties. Consequently, a plethora of entanglement quantifiers and classifiers exist, corresponding to different operational paradigms and mathematical techniques. However, for most quantum systems, exactly quantifying the amount of entanglement is extremely demanding, if at all possible. This is further exacerbated by the difficulty of experimentally controlling and measuring complex quantum states. Consequently, there are various approaches for experimentally detecting and certifying entanglement when exact quantification is not an option, with a particular focus on practically implementable methods and resource efficiency. The applicability and performance of these methods strongly depends on the assumptions one is willing to make regarding the involved quantum states and measurements, in short, on the available prior information about the quantum system. In this review we discuss the most commonly used paradigmatic quantifiers of entanglement. For these, we survey state-of-the-art detection and certification methods, including their respective underlying assumptions, from both a theoretical and experimental point of view.


title: Landau-Khalatnikov-Fradkin transformation and the mystery of even $\zeta$-values in Euclidean massless correlators

abstract: The Landau-Khalatnikov-Fradkin (LKF) transformation is a powerful and elegant transformation allowing to study the gauge dependence of the propagator of charged particles interacting with gauge fields. With the help of this transformation, we derive a non-perturbative identity between massless propagators in two different gauges. From this identity, we find that the corresponding perturbative series can be exactly expressed in terms of a hatted transcendental basis that eliminates all even Euler $\zeta$-functions. This explains the mystery of even $\zeta$-values observed in multi-loop calculations of Euclidean massless correlators for almost three decades now. Our construction further allows us to derive an exact formula relating hatted and standard $\zeta$-functions to all orders of perturbation theory.


title: Linearized 3-D Electromagnetic Contrast Source Inversion and Its Applications to Half-space Configurations

abstract: One of the main computational drawbacks in the application of 3-D iterative inversion techniques is the requirement of solving the field quantities for the updated contrast in every iteration. In this paper, the 3-D electromagnetic inverse scattering problem is put into a discretized finite-difference frequency-domain scheme and linearized into a cascade of two linear functionals. To deal with the nonuniqueness effectively, the joint structure of the contrast sources is exploited using a sum-of-$\ell_1$-norm optimization scheme. A cross-validation technique is used to check whether the optimization process is accurate enough. The total fields are, then, calculated and used to reconstruct the contrast by minimizing a cost functional defined as the sum of the data error and state error. In this procedure, the total fields in the inversion domain are computed only once, while the quality and accuracy of the obtained reconstructions are maintained. The novel method is applied to ground-penetrating radar imaging and through-the-wall imaging, in which the validity and efficiency of the method is demonstrated.


title: Study of phases in a holographic QCD model

abstract: Witten-Sakai-Sugimoto model is used to study two flavour Yang-Mills theory with large number of colours at finite temperature and in presence of chemical potential for baryon number and isospin. Sources for $U(1)_B$ and $U(1)_3$ gauge fields on the flavour 8-branes are D4-branes wrapped on $S^4$ part of the background. Here, gauge symmetry on the flavour branes has been decomposed as $U(2) \equiv U(1)_B \times SU(2)$ and $U(1)_3$ is within $SU(2)$ and generated by the diagonal generator. We show various brane configurations, along with the phases in the boundary theory they correspond to, and explore the possibility of phase transition between various pairs of phases.


title: Dual representations for systemic risk measures based on acceptance sets

abstract: We establish dual representations for systemic risk measures based on acceptance sets in a general setting. We deal with systemic risk measures of both "first allocate, then aggregate" and "first aggregate, then allocate" type. In both cases, we provide a detailed analysis of the corresponding systemic acceptance sets and their support functions. The same approach delivers a simple and self-contained proof of the dual representation of utility-based risk measures for univariate positions.


title: A double mean field equation related to a curvature prescription problem

abstract: We study a double mean field-type PDE related to a prescribed curvature problem on compacts surfaces with boundary. We provide a general blow-up analysis, then a Moser-Trudinger inequality, which gives energy-minimizing solutions for some range of parameters. Finally, we provide existence of min-max solutions for a wider range of parameters, which is dense in the plane if $\S$ is not simply connected.


title: Solving Statistical Mechanics on Sparse Graphs with Feedback Set Variational Autoregressive Networks

abstract: We propose a method for solving statistical mechanics problems defined on sparse graphs. It extracts a small Feedback Vertex Set (FVS) from the sparse graph, converting the sparse system to a much smaller system with many-body and dense interactions with an effective energy on every configuration of the FVS, then learns a variational distribution parameterized using neural networks to approximate the original Boltzmann distribution. The method is able to estimate free energy, compute observables, and generate unbiased samples via direct sampling without auto-correlation. Extensive experiments show that our approach is more accurate than existing approaches for sparse spin glasses. On random graphs and real-world networks, our approach significantly outperforms the standard methods for sparse systems such as the belief-propagation algorithm; on structured sparse systems such as two-dimensional lattices our approach is significantly faster and more accurate than recently proposed variational autoregressive networks using convolution neural networks.

categories: cond-mat.stat-mech cond-mat.dis-nn cs.LG stat.ML


title: Cyclic Flats of Binary Matroids

abstract: In this paper, first steps are taken towards characterising lattices of cyclic flats $\mathcal{Z}(M)$ that belong to matroids $M$ that can be represented over a prescribed finite field $\mathbb{F}_q$. Two natural maps from $\mathcal{Z}(M)$ to the lattice of cyclic flats of a minor of $M$ are given. Binary matroids are characterised via their lattice of cyclic flats. It is shown that the lattice of cyclic flats of a simple binary matroid without isthmuses is atomic.


title: Multiplicity and concentration results for a magnetic Schr\"{o}dinger equation with exponential critical growth in $\mathbb{R}^{2}$

abstract: In this paper we study the following nonlinear Schr\"{o}dinger equation with magnetic field \[ \Big(\frac{\varepsilon}{i}\nabla-A(x)\Big)^{2}u+V(x)u=f(| u|^{2})u,\quad x\in\mathbb{R}^{2}, \] where $\varepsilon>0$ is a parameter, $V:\mathbb{R}^{2}\rightarrow \mathbb{R}$ and $A: \mathbb{R}^{2}\rightarrow \mathbb{R}^{2}$ are continuous potentials and $f:\mathbb{R}\rightarrow \mathbb{R}$ has exponential critical growth. Under a local assumption on the potential $V$, by variational methods, penalization technique, and Ljusternick-Schnirelmann theory, we prove multiplicity and concentration of solutions for $\varepsilon$ small.


title: Design of patchy rhombi: from close-packed tilings to open lattices

abstract: In the realm of functional materials, the production of two-dimensional structures with tuneable porosity is of paramount relevance for many practical applications: surfaces with regular arrays of pores can be used for selective adsorption or immobilization of guest units that are complementary in shape and/or size to the pores, thus achieving, for instance, selective filtering or well-defined responses to external stimuli. The principles that govern the formation of such structures are valid at both the molecular and the colloidal scale. Here we provide simple design directions to combine the anisotropic shape of the building units -- either molecules or colloids -- and selective directional bonding. Using extensive computer simulations we show that regular rhombic platelets decorated with attractive and repulsive interaction sites lead to specific tilings, going smoothly from close-packed arrangements to open lattices. The rationale behind the rich tiling scenario observed can be described in terms of steric incompatibilities, unsatisfied bonding geometries and interplays between local and long-range order.


title: Some modular abelian surfaces

abstract: We use the main theorem of Boxer-Calegari-Gee-Pilloni (arXiv:1812.09269) to give explicit examples of modular abelian surfaces $A$ over $\mathbf{Q}$ without extra endomorhpisms such that $A$ has good reduction outside the primes 2, 3, 5, and 7.


title: Clausius inequality versus quantum coherence

abstract: In this study, we model a harmonic oscillator that enters an interferometer partially coupled to a thermal bath of oscillatory fields by employing a Brownian-type Lindblad master equation. More specifically, we investigate the dynamics and the variations of the thermodynamic quantities of the system at different temperatures. We recognize that although the system can remain coherent during its interaction with the thermal bath in the low-temperature limit, the system's entropy production violates the Clausius inequality. Furthermore, we argue that the system's coherence is the source of this violation, rather than the entanglement degree of system-environment, as reported in previous studies.


title: Investigating `dark' energy in the solar corona using forward modeling of MHD waves

abstract: It is now well established that the Alfv\'enic waves are ubiquitous in the solar corona. However, the Alfv\'enic wave energy estimated from the Doppler velocity measurements in the corona was found to be four orders of magnitude less than that estimated from non-thermal line widths. McIntosh & De Pontieu (2012) suggested that this discrepancy in energy might be due to the line-of-sight (LOS) superposition of the several oscillating structures, which can lead to an underestimation of the Alfv\'enic wave amplitudes and energies. McIntosh & De Pontieu (2012) termed this coronal `dark' or `hidden' energy. However, their simulations required the use of an additional, unknown source of Alfv\'enic wave energy to provide agreement with measurements of the coronal non-thermal line widths. In this study, we investigate the requirement of this unknown source of additional `dark' energy in the solar corona using gravitationally stratified 3D magnetohydrodynamic (MHD) simulations of propagating waves. We excite the transverse MHD waves and generate synthetic observations for the Fe XIII emission line. We establish that the LOS superposition greatly reduces the Doppler velocity amplitudes and increases the non-thermal line widths. Importantly, our model generates the observed wedge-shaped correlation between Doppler velocities and non-thermal line widths. We find that the observed wave energy is only 0.2-1\% of the true wave energy which explains 2-3 orders of magnitude of the energy discrepancy. We conclusively establish that the true wave energies are hidden in the non-thermal line widths. Hence, our results rule out the requirement for an additional `dark' energy in the solar corona.


title: Direct Quantification of Water Surface Charge by Phase-Sensitive Second Harmonic Spectroscopy

abstract: We develop and verify a phase-sensitive second harmonic generation spectroscopic scheme that allows for direct determination of the absolute surface charge density and surface potential of a water interface without need of prior interfacial information. The method relies on selective probing of surface-field-induced reorientation order of water molecules in the electrical double layer and is, hence, independent of the interfacial molecular bonding structure. Application of this technique to a mixed surfactant monolayer on water suggests the manifest effect of the chain-chain interactions among the monolayer on adsorption of soluble ionic surfactants. We also deduce the third-order nonlinear susceptibility of bulk water and prove its applicability to analysis of charges of various water interfaces. In addition, we show that the Debye-H\"uckle theory should be avoided in the spectroscopic analysis for its potential significant error, as evidenced experimentally and theoretically.


title: Heuristic Approach Towards Countermeasure Selection using Attack Graphs

abstract: Selecting the optimal set of countermeasures is a challenging task that involves various considerations and tradeoffs such as prioritizing the risks to mitigate and costs. The vast majority of studies for selecting a countermeasure deployment are based on a limited risk assessment procedure that utilizes the common vulnerability scoring system (CVSS). Such a risk assessment procedure does not necessarily consider the prerequisites and exploitability of a specific asset, cannot distinguish insider from outsider threat actor, and does not express the consequences of exploiting a vulnerability as well as the attacker's lateral movements. Other studies applied a more extensive risk assessment procedure that relies on manual work and repeated assessment. These solutions however, do not consider the network topology and do not specify the optimal position for deploying the countermeasures, and therefore are less practical. In this paper we suggest a heuristic search approach for selecting the optimal countermeasure deployment under a given budget limitation. The proposed method expresses the risk of the system using an extended attack graph modeling, which considers the prerequisites and consequences of exploiting a vulnerability, examines the attacker's potential lateral movements, and express the physical network topology as well as vulnerabilities in network protocols. In addition, unlike previous studies which utilizes attack graph for countermeasure planning, the proposed methods does not require re-generating the attack graph at each stage of the procedure, which is computationally heavy, and therefore it provides a more accurate and practical countermeasure deployment planning process.


title: A High-Performance Implementation of a Robust Preconditioner for Heterogeneous Problems

abstract: We present an efficient implementation of the highly robust and scalable GenEO preconditioner in the high-performance PDE framework DUNE. The GenEO coarse space is constructed by combining low energy solutions of a local generalised eigenproblem using a partition of unity. In this paper we demonstrate both weak and strong scaling for the GenEO solver on over 15,000 cores by solving an industrially motivated problem with over 200 million degrees of freedom. Further, we show that for highly complex parameter distributions arising in certain real-world applications, established methods become intractable while GenEO remains fully effective. The purpose of this paper is two-fold: to demonstrate the robustness and high parallel efficiency of the solver and to document the technical details that are crucial to the efficiency of the code.


title: Sliding Mode Control Techniques and Artificial Potential Field for Dynamic Collision Avoidance in Rendezvous Maneuvers

abstract: The paper considers autonomous rendezvous maneuver and proximity operations of two spacecraft in presence of obstacles. A strategy that combines guidance and control algorithms is analyzed. The proposed closed-loop system is able to guarantee a safe path in a real environment, as well as robustness with respect to external disturbances and dynamic obstacles. The guidance strategy exploits a suitably designed Artificial Potential Field (APF), while the controller relies on Sliding Mode Control (SMC), for both position and attitude tracking of the spacecraft. As for the position control, two different first order SMC methods are considered, namely the component-wise and the simplex-based control techniques. The proposed integrated guidance and control strategy is validated by extensive simulations performed with a six degree-of-freedom (DOF) orbital simulator and appears suitable for real-time control with minimal on-board computational effort. Fuel consumption and control effort are evaluated, including different update frequencies of the closed-loop software.


title: A unifying computations of Whittle's Index for Markovian bandits

abstract: The multi-armed restless bandit framework allows to model a wide variety of decision-making problems in areas as diverse as industrial engineering, computer communication, operations research, financial engineering, communication networks etc. In a seminal work, Whittle developed a methodology to derive well-performing (Whittle's) index policies that are obtained by solving a relaxed version of the original problem. However, the computation of Whittle's index itself is a difficult problem and hence researchers focused on calculating Whittle's index numerically or with a problem dependent approach. In our main contribution we derive an analytical expression for Whittle's index for any Markovian bandit with both finite and infinite transition rates. We derive sufficient conditions for the optimal solution of the relaxed problem to be of threshold type, and obtain conditions for the bandit to be indexable, a property assuring the existence of Whittle's index. Our solution approach provides a unifying expression for Whittle's index, which we highlight by retrieving known indices from literature as particular cases. The applicability of finite rates is illustrated with the machine repairmen problem, and that of infinite rates by an example of communication networks where transmission rates react instantaneously to packet losses.