title: Geometric Sparsification of Closeness Relations: Eigenvalue Clustering for Computing Matrix Functions

abstract: We show how to efficiently solve a clustering problem that arises in a method to evaluate functions of matrices. The problem requires finding the connected components of a graph whose vertices are eigenvalues of a real or complex matrix and whose edges are pairs of eigenvalues that are at most \delta away from each other. Davies and Higham proposed solving this problem by enumerating the edges of the graph, which requires at least $\Omega(n^{2})$ work. We show that the problem can be solved by computing the Delaunay triangulation of the eigenvalues, removing from it long edges, and computing the connected components of the remaining edges in the triangulation. This leads to an $O(n\log n)$ algorithm. We have implemented both algorithms using CGAL, a mature and sophisticated computational-geometry software library, and we demonstrate that the new algorithm is much faster in practice than the naive algorithm. We also present a tight analysis of the naive algorithm, showing that it performs $\Theta(n^{2})$ work, and correct a misrepresentation in the original statement of the problem. To the best of our knowledge, this is the first application of computational geometry to solve a real-world problem in numerical linear algebra.


title: robROSE: A robust approach for dealing with imbalanced data in fraud detection

abstract: A major challenge when trying to detect fraud is that the fraudulent activities form a minority class which make up a very small proportion of the data set. In most data sets, fraud occurs in typically less than 0.5% of the cases. Detecting fraud in such a highly imbalanced data set typically leads to predictions that favor the majority group, causing fraud to remain undetected. We discuss some popular oversampling techniques that solve the problem of imbalanced data by creating synthetic samples that mimic the minority class. A frequent problem when analyzing real data is the presence of anomalies or outliers. When such atypical observations are present in the data, most oversampling techniques are prone to create synthetic samples that distort the detection algorithm and spoil the resulting analysis. A useful tool for anomaly detection is robust statistics, which aims to find the outliers by first fitting the majority of the data and then flagging data observations that deviate from it. In this paper, we present a robust version of ROSE, called robROSE, which combines several promising approaches to cope simultaneously with the problem of imbalanced data and the presence of outliers. The proposed method achieves to enhance the presence of the fraud cases while ignoring anomalies. The good performance of our new sampling technique is illustrated on simulated and real data sets and it is shown that robROSE can provide better insight in the structure of the data. The source code of the robROSE algorithm is made freely available.


title: Analysis of Greenhouse Gases

abstract: Climate change is a result of a complex system of interactions of greenhouse gases (GHG), the ocean, land, ice, and clouds. Large climate change models use several computers and solve several equations to predict the future climate. The equations may include simple polynomials to partial differential equations. Because of the uptake mechanism of the land and ocean, greenhouse gas emissions can take a while to affect the climate. The IPCC has published reports on how greenhouse gas emissions may affect the average temperature of the troposphere and the predictions show that by the end of the century, we can expect a temperature increase from 0.8 C to 5 C. In this article, I use Linear Regression (LM), Quadratic Regression and Gaussian Process Regression (GPR) on monthly GHG data going back several years and try to predict the temperature anomalies based on extrapolation. The results are quite similar to the IPCC reports.


title: Adversarial Examples and the Deeper Riddle of Induction: The Need for a Theory of Artifacts in Deep Learning

abstract: Deep learning is currently the most widespread and successful technology in artificial intelligence. It promises to push the frontier of scientific discovery beyond current limits. However, skeptics have worried that deep neural networks are black boxes, and have called into question whether these advances can really be deemed scientific progress if humans cannot understand them. Relatedly, these systems also possess bewildering new vulnerabilities: most notably a susceptibility to "adversarial examples". In this paper, I argue that adversarial examples will become a flashpoint of debate in philosophy and diverse sciences. Specifically, new findings concerning adversarial examples have challenged the consensus view that the networks' verdicts on these cases are caused by overfitting idiosyncratic noise in the training set, and may instead be the result of detecting predictively useful "intrinsic features of the data geometry" that humans cannot perceive (Ilyas et al., 2019). These results should cause us to re-examine responses to one of the deepest puzzles at the intersection of philosophy and science: Nelson Goodman's "new riddle" of induction. Specifically, they raise the possibility that progress in a number of sciences will depend upon the detection and manipulation of useful features that humans find inscrutable. Before we can evaluate this possibility, however, we must decide which (if any) of these inscrutable features are real but available only to "alien" perception and cognition, and which are distinctive artifacts of deep learning-for artifacts like lens flares or Gibbs phenomena can be similarly useful for prediction, but are usually seen as obstacles to scientific theorizing. Thus, machine learning researchers urgently need to develop a theory of artifacts for deep neural networks, and I conclude by sketching some initial directions for this area of research.


title: Partially volume expanding diffeomorphisms

abstract: We call a partially hyperbolic diffeomorphism \emph{partially volume expanding} if the Jacobian restricted to any hyperplane that contains the unstable bundle $E^u$ is larger than $1$. This is a $C^1$ open property. We show that any $C^{1+}$ partially volume expanding diffeomorphisms admits finitely many physical measures, the union of whose basins has full volume.


title: Counterfactual Policy Evaluation for Decision-Making in Autonomous Driving

abstract: Learning-based approaches, such as reinforcement and imitation learning are gaining popularity in decision-making for autonomous driving. However, learned policies often fail to generalize and cannot handle novel situations well. Asking and answering questions in the form of "Would a policy perform well if the other agents had behaved differently?" can shed light on whether a policy has seen similar situations during training and generalizes well. In this work, a counterfactual policy evaluation is introduced that makes use of counterfactual worlds - worlds in which the behaviors of others are non-actual. If a policy can handle all counterfactual worlds well, it either has seen similar situations during training or it generalizes well and is deemed to be fit enough to be executed in the actual world. Additionally, by performing the counterfactual policy evaluation, causal relations and the influence of changing vehicle's behaviors on the surrounding vehicles becomes evident. To validate the proposed method, we learn a policy using reinforcement learning for a lane merging scenario. In the application-phase, the policy is only executed after the counterfactual policy evaluation has been performed and if the policy is found to be safe enough. We show that the proposed approach significantly decreases the collision-rate whilst maintaining a high success-rate.


title: A simple Stochastic SIR model for COVID 19 Infection Dynamics for Karnataka: Learning from Europe

abstract: In this short note we model the region-wise trends of the evolution to COVID-19 infections using a stochastic SIR model. The SIR dynamics are expressed using \textit{It\^o-stochastic differential equations}. We first derive the parameters of the model from the available daily data from European regions based on a 24-day history of infections, recoveries and deaths. The derived parameters have been aggregated to project future trends for the Indian subcontinent, which is currently at an early stage in the infection cycle. The projections are meant to serve as a guideline for strategizing the socio-political counter measures to mitigate COVID-19.


title: A Modified Global Climate Simulation Model

abstract: In this paper, we incorporate seasonal variations of insolation into the global climate model C-GOLDSTEIN. We use a new approach for modelling insolation from the space perspective presented in the authors' earlier work and build it into the existing climate model. Realistic monthly temperature distributions have been obtained after running C-GOLDSTEIN with the new insolation component. Also, the average accuracy of modelling the insolation within the model has been increased by 2%. In addition, new types of experiments can now be performed with C-GOLDSTEIN, such as the investigation of consequences of random variations of insolation on temperature etc.


title: Rat big, cat eaten! Ideas for a useful deep-agent protolanguage

abstract: Deep-agent communities developing their own language-like communication protocol are a hot (or at least warm) topic in AI. Such agents could be very useful in machine-machine and human-machine interaction scenarios long before they have evolved a protocol as complex as human language. Here, I propose a small set of priorities we should focus on, if we want to get as fast as possible to a stage where deep agents speak a useful protolanguage.


title: Time series and machine learning to forecast the water quality from satellite data

abstract: Managing the quality of water for present and future generations of coastal regions should be a central concern of both citizens and public officials. Remote sensing can contribute to the management and monitoring of coastal water and pollutants. Algal blooms are a coastal pollutant that is a cause of concern. Many satellite data, such as MODIS, have been used to generate water-quality products to detect the blooms such as chlorophyll a (Chl-a), a photosynthesis index called fluorescence line height (FLH), and sea surface temperature (SST). It is important to characterize the spatial and temporal variations of these water quality products by using the mathematical models of these products. However, for monitoring, pollution control boards will need nowcasts and forecasts of any pollution. Therefore, we aim to predict the future values of the MODIS Chl-a, FLH, and SST of the water. This will not be limited to one type of water but, rather, will cover different types of water varying in depth and turbidity. This is very significant because the temporal trend of Chl-a, FLH, and SST is dependent on the geospatial and water properties. For this purpose, we will decompose the time series of each pixel into several components: trend, intra-annual variations, seasonal cycle, and stochastic stationary. We explore three such time series machine learning models that can characterize the non-stationary time series data and predict future values, including the Seasonal ARIMA (Auto Regressive Integrated Moving Average) (SARIMA), regression, and neural network. The results indicate that all these methods are effective at modelling Chl-a, FLH, and SST time series and predicting the values reasonably well. However, regression and neural network are found to be the best at predicting Chl-a in all types of water (turbid and shallow). Meanwhile, the SARIMA model provides the best prediction of FLH and SST.


title: Exploring the Sharma-Mittal HDE models with different diagnostic tools

abstract: In this paper, we have examined the Sharma-Mittal holographic dark energy model (SMHDE) in the framework of an isotropic and spatially homogeneous flat Friedmann-Robertson-Walker(FRW) Universe by considering different values of parameter $\delta$ and $R$, where the infrared cut-off is taken care by the Hubble horizon. We examined the SMHDE model through the analysis of Statefinder hierarchy and the growth rate of perturbation. The evolutionary trajectories of the statefinder hierarchy $S_3^1$, $S_3^2$ $S_4^1$, $S_4^2$ versus redshift z, show satisfactory behavior throughout the Universe evaluation. One promising tool for investigating the dark energy models is the composite null diagnostic(CND) $\{ S_3^1 - \epsilon\}$, where the evolutionary trajectories of the $S_3^1 - \epsilon$ pair present different property and the departure from $\Lambda$CDM could be well evaluated. Additionally, we investigated the dynamical analysis of the model by $\omega_{D}-\omega^{'}_{D}$ pair analysis.


title: Probabilistic magnetometry with two-spin system in diamond

abstract: Solid-state magnetometers like the Nitrogen-Vacancy center in diamond have been of paramount importance for the development of quantum sensing with nanoscale spatial resolution. The basic protocol is a Ramsey sequence, that imprints an external static magnetic field into phase of the quantum sensor, which is subsequently readout. In this work we show that the hyperfine coupling between the Nitrogen-Vacancy and a nearby Carbon-13 can be used to set a post-selection protocol that leads to an enhancement of the sensitivity under realistic experimental conditions. We found that for an isotopically purified sample the detection of weak magnetic fields in the $\mu$T range can be achieved with a sensitivity of few nTHz$^{-1/2}$ at cryogenic temperature ($4$ K), and $0.1$ $\mu$THz$^{-1/2}$ at room temperature.


title: Twisting on pre-Lie algebras and quasi-pre-Lie bialgebras

abstract: We study (quasi-)twilled pre-Lie algebras and the associated $L_\infty$-algebras and differential graded Lie algebras. Then we show that certain twisting transformations on (quasi-)twilled pre-Lie algbras can be characterized by the solutions of Maurer-Cartan equations of the associated differential graded Lie algebras ($L_\infty$-algebras). Furthermore, we show that $\mathcal{O}$-operators and twisted $\mathcal{O}$-operators are solutions of the Maurer-Cartan equations. As applications, we study (quasi-)pre-Lie bialgebras using the associated differential graded Lie algebras ($L_\infty$-algebras) and the twisting theory of (quasi-)twilled pre-Lie algebras. In particular, we give a construction of quasi-pre-Lie bialgebras using symplectic Lie algebras, which is parallel to that a Cartan $3$-form on a semi-simple Lie algebra gives a quasi-Lie bialgebra.


title: Improving data-driven global weather prediction using deep convolutional neural networks on a cubed sphere

abstract: We present a significantly-improved data-driven global weather forecasting framework using a deep convolutional neural network (CNN) to forecast several basic atmospheric variables on a global grid. New developments in this framework include an offline volume-conservative mapping to a cubed-sphere grid, improvements to the CNN architecture, and the minimization of the loss function over multiple steps in a prediction sequence. The cubed-sphere remapping minimizes the distortion on the cube faces on which convolution operations are performed and provides natural boundary conditions for padding in the CNN. Our improved model produces weather forecasts that are indefinitely stable and produce realistic weather patterns at lead times of several weeks and longer. For short- to medium-range forecasting, our model significantly outperforms persistence, climatology, and a coarse-resolution dynamical numerical weather prediction (NWP) model. Unsurprisingly, our forecasts are worse than those from a high-resolution state-of-the-art operational NWP system. Our data-driven model is able to learn to forecast complex surface temperature patterns from few input atmospheric state variables. On annual time scales, our model produces a realistic seasonal cycle driven solely by the prescribed variation in top-of-atmosphere solar forcing. Although it is currently less accurate than operational weather forecasting models, our data-driven CNN executes much faster than those models, suggesting that machine learning could prove to be a valuable tool for large-ensemble forecasting.


title: Zero-Assignment Constraint for Graph Matching with Outliers

abstract: Graph matching (GM), as a longstanding problem in computer vision and pattern recognition, still suffers from numerous cluttered outliers in practical applications. To address this issue, we present the zero-assignment constraint (ZAC) for approaching the graph matching problem in the presence of outliers. The underlying idea is to suppress the matchings of outliers by assigning zero-valued vectors to the potential outliers in the obtained optimal correspondence matrix. We provide elaborate theoretical analysis to the problem, i.e., GM with ZAC, and figure out that the GM problem with and without outliers are intrinsically different, which enables us to put forward a sufficient condition to construct valid and reasonable objective function. Consequently, we design an efficient outlier-robust algorithm to significantly reduce the incorrect or redundant matchings caused by numerous outliers. Extensive experiments demonstrate that our method can achieve the state-of-the-art performance in terms of accuracy and efficiency, especially in the presence of numerous outliers.


title: Peano continua with self regenerating fractals

abstract: We deal with the question of Masayoshi Hata: is every Peano continuum a topological fractal? A compact space $X$ is a topological fractal if there exists $\mathcal{F}$ a finite family of self-maps on $X$ such that $X=\bigcup_{f\in\mathcal{F}}f(X)$ and for every open cover $\mathcal{U}$ of $X$ there is $n\in\mathbb{N}$ such that for all maps $f_1,\dots,f_n\in\mathcal{F}$ the set $f_1\circ\dots\circ f_n(X)$ is contained in some set $U\in\mathcal{U}$. In the paper we present some idea how to extend a topological fractal and we show that a Peano continuum is a topological fractal if it contains so-called self regenerating fractal with nonempty interior. A Hausdorff topological space $A$ is a self regenerating fractal if for every non-empty open subset $U$, $A$ is a topological fractal for some family of maps constant on $A\setminus U$. The notion of self regenerating fractal much better reflects the intuitive perception of self-similarity. We present some classical fractals which are self regenerating.


title: Beaming as an explanation of the repetition/width relation in FRBs

abstract: It is currently not known if repeating fast radio bursts (FRBs) are fundamentally different from those that have not been seen to repeat. One striking difference between repeaters and apparent non-repeaters in the CHIME sample is that the once-off events are typically shorter in duration than sources that have been detected two or more times. We offer a simple explanation for this discrepancy based on a selection effect due to beamed emission, in which highly-beamed FRBs are less easily observed to repeat, but are abundant enough to detect often as once-off events. The explanation predicts that there is a continuous distribution of burst duration---not a static bimodal one---with a correlation between repetition rate and width. Pulse width and opening angle may be related by relativistic effects in shocks, where short-duration bursts have small solid angles due to a large common Lorentz factor. Alternatively, the relationship could be a geometric effect where narrow beams sweep past the observer more quickly, as with pulsars. Our model has implications for the FRB emission mechanism and energy scale, volumetric event rates, and the application of FRBs to cosmology.


title: Triad State Space Construction for Chaotic Signal Classification with Deep Learning

abstract: Inspired by the well-known permutation entropy (PE), an effective image encoding scheme for chaotic time series, Triad State Space Construction (TSSC), is proposed. The TSSC image can recognize higher-order temporal patterns and identify new forbidden regions in time series motifs beyond the Bandt-Pompe probabilities. The Convolutional Neural Network (ConvNet) is widely used in image classification. The ConvNet classifier based on TSSC images (TSSC-ConvNet) are highly accurate and very robust in the chaotic signal classification.


title: Properties of heavy ion beams produced by a PW sub-picosecond laser

abstract: The paper reports the results of numerical studies on the generation of heavy ions from a sub-micrometre gold target irradiated by a high-energy (60 J) high-intensity (~ 2 x 10^20 W/cm2) sub-ps laser pulse. The properties of the heavy ion beam produced from a 0.5-um gold target with and without hydrogen contaminants are investigated using a multi-dimensional (2D3V) particle-in-cell code. It was found that in the case of the pure (without the contaminants) target, the ion beam contains over 20 ion species with a charge state > 30 and a total number of high-energy (> 20 MeV) ions reaching values well above 10^12. The ion energy spectrum is broad, with maximum energies > 1 GeV and mean energies > 100 MeV. At a distance ~ 1 mm from the target, the intensity of the ion beam is ~ 10^15 W/cm2, the ion fluence reaches ~ 10^16 ions/cm2, and the ion pulse duration is ~ 100 ps. The contaminants significantly reduce the ion beam parameters and both the mean ion energy and the intensity and fluence of the ion beam generated from the contaminated gold target are almost an order of magnitude lower than those produced from the pure target. Heavy ion beams with the parameters obtained in the case of the pure gold target are barely achievable in conventional RF-driven accelerators, so they can open the door to new areas of research, in particular in high energy-density physics.


title: Spin textures in chiral magnetic monolayers with suppressed nearest-neighbor exchange

abstract: High tunability of two dimensional magnetic materials (by strain, gating, heterostructuring or otherwise) provides unique conditions for studying versatile magnetic properties and controlling emergent magnetic phases. Expanding the scope of achievable magnetic phenomena in such materials is important for both fundamental and technological advances. Here we perform atomistic spin-dynamics simulations to explore the (chiral) magnetic phases of atomic monolayers in the limit of suppressed first-neighbors exchange interaction. We report the rich phase diagram of exotic magnetic configurations, obtained for both square and honeycomb lattice symmetries, comprising coexistence of ferromagnetic and antiferromagnetic spin-cycloids, as well as multiple types of magnetic skyrmions. We perform a minimum-energy path analysis for the skyrmion collapse to evaluate the stability of such topological objects, and reveal that magnetic monolayers could be good candidates to host the antiferromagnetic skyrmions that are experimentally evasive to date.


title: Stability Analysis of Droop-Controlled Inverter-Based Power Grids via Timescale Separation

abstract: We consider the problem of stability analysis for distribution grids with droop-controlled inverters and dynamic distribution power lines. The inverters are modeled as voltage sources with controllable frequency and amplitude. This problem is very challenging for large networks as numerical simulations and detailed eigenvalue analysis are impactical. Motivated by the above limitations, we present in this paper a systematic and computationally efficient framework for stability analysis of inverter-based distribution grids. To design our framework, we use tools from singular perturbation and Lyapunov theories. Interestingly, we show that stability of the fast dynamics of the power grid depends only on the voltage droop gains of the inverters while, stability of the slow dynamics, depends on both voltage and frequency droop gains. Finally, by leveraging these timescale separation properties, we derive sufficient conditions on the frequency and voltage droop gains of the inverters that warrant stability of the full system. We illustrate our theoretical results through a numerical example on the IEEE 13-bus distribution grid.


title: Full counting statistics and fluctuation-dissipation relation for periodically driven two-state systems

abstract: We derive the fluctuation theorem for a stochastic and periodically driven system coupled to two reservoirs with the aid of a master equation. We write down the cumulant generating functions for both the current and entropy production in closed compact forms so as to treat the adiabatic and nonadiabatic contributions systematically. We derive the fluctuation theorem by taking into account the time reversal symmetry and the property that the instantaneous currents flowing into the left and the right reservoir are not equal. It is found that the fluctuation-dissipation relation derived from the fluctuation theorem involves an expansion with respect to the time derivative of the affinity.


title: Cryptography using generalized Fibonacci matrices with Affine-Hill cipher

abstract: In this article, we have proposed a public key cryptography using Affine-Hill cipher with a generalized Fibonacci matrix(called multinacci matrix). Also proposed a key establishment(exchange of key matrix $K=Q_{\lambda}^{k}$ of order $\lambda\times\lambda$ for encryption-decryption) scheme with the help of multinacci sequences under prime modulo. In this scheme, instead of exchanging key matrix, we need to exchange the only pair of numbers $(\lambda, k)$, which reduces the time complexity as well as space complexity and comes with a large key-space.


title: Rotation-activity relations and flares of M dwarfs with K2 long- and short-cadence data

abstract: Using light curves obtained by the K2 mission, we study the relation between stellar rotation and magnetic activity with special focus on stellar flares. Our sample comprises 56 bright and nearby M dwarfs observed by K2 during campaigns C0-C18 in long- and short-cadence mode. We derive rotation periods for 46 M dwarfs and measure photometric activity indicators such as amplitude of the rotational signal, standard deviation of the light curves, and the basic flare properties (flare rate, flare energy, flare duration, and flare amplitude). We found 1662 short-cadence flares, 363 of which have a long-cadence counterpart with flare energies of up to $5.6\cdotp10^{34}$erg. The flare amplitude, duration, and frequency derived from the short-cadence light curves differ significantly from those derived from the long-cadence data. The analysis of the short-cadence light curves results in a flare rate that is 4.6 times higher than the long-cadence data. We confirm the abrupt change in activity level in the rotation-activity relation at a critical period of ~10d when photometric activity diagnostics are used. This change is most drastic in the flare duration and frequency for short-cadence data. Our flare studies revealed that the highest flare rates are not found among the fastest rotators and that stars with the highest flare rates do not show the most energetic flares. We found that the superflare frequency ($E\geq5\cdotp10^{34}$erg) for the fast-rotating M stars is twice higher than for solar like stars in the same period range. By fitting the cumulative FFD, we derived a power-law index of $\alpha=1.84 \pm 0.14$, consistent with previous M dwarf studies and the value found for the Sun.


title: Toward source region tomography with inter-source interferometry: Shear wave velocity from 2018 West Bohemia swarm earthquakes

abstract: The concept of seismic interferometry embraces the construction of waves traveling between receivers or sources with cross-correlation techniques. In the present study cross-correlations of coda waves are used to measure travel times of shear waves between earthquake locations for five event clusters of the 2018 West Bohemia earthquake swarm. With the help of a high quality earthquake catalog, I was able to determine the shear wave velocity in the region of the five clusters separately. The shear wave velocities range between 3.5 km/s and 4.2 km/s. The resolution of this novel method is given by the extent of the clusters and better than for a comparable classical tomography. It is suggested to use the method in a tomographic inversion and map the shear wave velocity in the source region with unprecedented resolution. Furthermore, the influence of focal mechanisms and the attenuation properties on the polarity and location of the maxima in the cross-correlation functions is discussed. The intra-cluster ratio of P wave to S wave velocity is approximately fixed at 1.68.


title: Advances in Bayesian Probabilistic Modeling for Industrial Applications

abstract: Industrial applications frequently pose a notorious challenge for state-of-the-art methods in the contexts of optimization, designing experiments and modeling unknown physical response. This problem is aggravated by limited availability of clean data, uncertainty in available physics-based models and additional logistic and computational expense associated with experiments. In such a scenario, Bayesian methods have played an impactful role in alleviating the aforementioned obstacles by quantifying uncertainty of different types under limited resources. These methods, usually deployed as a framework, allows decision makers to make informed choices under uncertainty while being able to incorporate information on the the fly, usually in the form of data, from multiple sources while being consistent with the physical intuition about the problem. This is a major advantage that Bayesian methods bring to fruition especially in the industrial context. This paper is a compendium of the Bayesian modeling methodology that is being consistently developed at GE Research. The methodology, called GE's Bayesian Hybrid Modeling (GEBHM), is a probabilistic modeling method, based on the Kennedy and O'Hagan framework, that has been continuously scaled-up and industrialized over several years. In this work, we explain the various advancements in GEBHM's methods and demonstrate their impact on several challenging industrial problems.


title: A general framework for causal classification

abstract: In many applications, there is a need to predict the effect of an intervention on different individuals from data. For example, which customers are persuadable by a product promotion? which patients should be treated with a certain type of treatment? These are typical causal questions involving the effect or the change in outcomes made by an intervention. The questions cannot be answered with traditional classification methods as they only use associations to predict outcomes. For personalised marketing, these questions are often answered with uplift modelling. The objective of uplift modelling is to estimate causal effect, but its literature does not discuss when the uplift represents causal effect. Causal heterogeneity modelling can solve the problem, but its assumption of unconfoundedness is untestable in data. So practitioners need guidelines in their applications when using the methods. In this paper, we use causal classification for a set of personalised decision making problems, and differentiate it from classification. We discuss the conditions when causal classification can be resolved by uplift (and causal heterogeneity) modelling methods. We also propose a general framework for causal classification, by using off-the-shelf supervised methods for flexible implementations. Experiments have shown two instantiations of the framework work for causal classification and for uplift (causal heterogeneity) modelling, and are competitive with the other uplift (causal heterogeneity) modelling methods.


title: AliExpress Learning-To-Rank: Maximizing Online Model Performance without Going Online

abstract: Learning-to-rank (LTR) has become a key technology in E-commerce applications. Most existing LTR approaches follow a supervised learning paradigm from offline labeled data collected from the online system. However, it has been noticed that previous LTR models can have a good validation performance over offline validation data but have a poor online performance, and vice versa, which implies a possible large inconsistency between the offline and online evaluation. We investigate and confirm in this paper that such inconsistency exists and can have a significant impact on AliExpress Search. Reasons for the inconsistency include the ignorance of item context during the learning, and the offline data set is insufficient for learning the context. Therefore, this paper proposes an evaluator-generator framework for LTR with item context. The framework consists of an evaluator that generalizes to evaluate recommendations involving the context, and a generator that maximizes the evaluator score by reinforcement learning, and a discriminator that ensures the generalization of the evaluator. Extensive experiments in simulation environments and AliExpress Search online system show that, firstly, the classic data-based metrics on the offline dataset can show significant inconsistency with online performance, and can even be misleading. Secondly, the proposed evaluator score is significantly more consistent with the online performance than common ranking metrics. Finally, as the consequence, our method achieves a significant improvement (\textgreater$2\%$) in terms of Conversion Rate (CR) over the industrial-level fine-tuned model in online A/B tests.


title: Towards Backward-Compatible Representation Learning

abstract: We propose a way to learn visual features that are compatible with previously computed ones even when they have different dimensions and are learned via different neural network architectures and loss functions. Compatible means that, if such features are used to compare images, then "new" features can be compared directly to "old" features, so they can be used interchangeably. This enables visual search systems to bypass computing new features for all previously seen images when updating the embedding models, a process known as backfilling. Backward compatibility is critical to quickly deploy new embedding models that leverage ever-growing large-scale training datasets and improvements in deep learning architectures and training methods. We propose a framework to train embedding models, called backward-compatible training (BCT), as a first step towards backward compatible representation learning. In experiments on learning embeddings for face recognition, models trained with BCT successfully achieve backward compatibility without sacrificing accuracy, thus enabling backfill-free model updates of visual embeddings.


title: Averaging Principle on Infinite Intervals for Stochastic Ordinary Differential Equations

abstract: In contrast to existing works on stochastic averaging on finite intervals, we establish an averaging principle on the whole real axis, i.e. the so-called second Bogolyubov theorem, for semilinear stochastic ordinary differential equations in Hilbert space with Poisson stable (in particular, periodic, quasi-periodic, almost periodic, almost automorphic etc) coefficients. Under some appropriate conditions we prove that there exists a unique recurrent solution to the original equation, which possesses the same recurrence property as the coefficients, in a small neighborhood of the stationary solution to the averaged equation, and this recurrent solution converges to the stationary solution of averaged equation uniformly on the whole real axis when the time scale approaches zero.


title: From geometry to arithmetic of compact hyperbolic Coxeter polytopes

abstract: We establish some geometric constraints on compact Coxeter polytopes in hyperbolic spaces and show that these constraints can be a very useful tool for the classification problem of reflective anisotropic Lorentzian lattices and cocompact arithmetic hyperbolic reflection groups.


title: Quantum Semantic Learning by Reverse Annealing an Adiabatic Quantum Computer

abstract: Boltzmann Machines constitute a class of neural networks with applications to image reconstruction, pattern classification and unsupervised learning in general. Their most common variants, called Restricted Boltzmann Machines (RBMs) exhibit a good trade-off between computability on existing silicon-based hardware and generality of possible applications. Still, the diffusion of RBMs is quite limited, since their training process proves to be hard. The advent of commercial Adiabatic Quantum Computers (AQCs) raised the expectation that the implementations of RBMs on such quantum devices could increase the training speed with respect to conventional hardware. To date, however, the implementation of RBM networks on AQCs has been limited by the low qubit connectivity when each qubit acts as a node of the neural network. Here we demonstrate the feasibility of a complete RBM on AQCs, thanks to an embedding that associates its nodes to virtual qubits, thus outperforming previous implementations based on incomplete graphs. Moreover, to accelerate the learning, we implement a semantic quantum search which, contrary to previous proposals, takes the input data as initial boundary conditions to start each learning step of the RBM, thanks to a reverse annealing schedule. Such an approach, unlike the more conventional forward annealing schedule, allows sampling configurations in a meaningful neighborhood of the training data, mimicking the behavior of the classical Gibbs sampling algorithm. We show that the learning based on reverse annealing quickly raises the sampling probability of a meaningful subset of the set of the configurations. Even without a proper optimization of the annealing schedule, the RBM semantically trained by reverse annealing achieves better scores on reconstruction tasks.


title: Effective transmission conditions for reaction-diffusion processes in domains separated by thin channels

abstract: We consider a reaction--diffusion equation in a domain consisting of two bulk regions connected via small channels periodically distributed within a thin layer. The height and the thickness of the channels are of order $\epsilon$, and the equation inside the layer depends on the parameter $\epsilon$ and an additional parameter $\gamma \in [-1,1)$, which describes the size of the diffusion in the layer. We derive effective models for the limit $\epsilon \to 0 $, when the channel-domain is replaced by an interface between the two bulk-domains.


title: On the worst-case error of least squares algorithms for $L_2$-approximation with high probability

abstract: It was recently shown in [4] that, for $L_2$-approximation of functions from a Hilbert space, function values are almost as powerful as arbitrary linear information, if the approximation numbers are square-summable. That is, we showed that \[ e_n \,\lesssim\, \sqrt{\frac{1}{k_n} \sum_{j\geq k_n} a_j^2} \qquad \text{ with }\quad k_n \asymp \frac{n}{\ln(n)}, \] where $e_n$ are the sampling numbers and $a_k$ are the approximation numbers. In particular, if $(a_k)\in\ell_2$, then $e_n$ and $a_n$ are of the same polynomial order. For this, we presented an explicit (weighted least squares) algorithm based on i.i.d. random points and proved that this works with positive probability. This implies the existence of a good deterministic sampling algorithm. Here, we present a modification of the proof in [4] that shows that the same algorithm works with probability at least $1-{n^{-c}}$ for all $c>0$.


title: Bag of biterms modeling for short texts

abstract: Analyzing texts from social media encounters many challenges due to their unique characteristics of shortness, massiveness, and dynamic. Short texts do not provide enough context information, causing the failure of the traditional statistical models. Furthermore, many applications often face with massive and dynamic short texts, causing various computational challenges to the current batch learning algorithms. This paper presents a novel framework, namely Bag of Biterms Modeling (BBM), for modeling massive, dynamic, and short text collections. BBM comprises of two main ingredients: (1) the concept of Bag of Biterms (BoB) for representing documents, and (2) a simple way to help statistical models to include BoB. Our framework can be easily deployed for a large class of probabilistic models, and we demonstrate its usefulness with two well-known models: Latent Dirichlet Allocation (LDA) and Hierarchical Dirichlet Process (HDP). By exploiting both terms (words) and biterms (pairs of words), the major advantages of BBM are: (1) it enhances the length of the documents and makes the context more coherent by emphasizing the word connotation and co-occurrence via Bag of Biterms, (2) it inherits inference and learning algorithms from the primitive to make it straightforward to design online and streaming algorithms for short texts. Extensive experiments suggest that BBM outperforms several state-of-the-art models. We also point out that the BoB representation performs better than the traditional representations (e.g, Bag of Words, tf-idf) even for normal texts.


title: Top Comment or Flop Comment? Predicting and Explaining User Engagement in Online News Discussions

abstract: Comment sections below online news articles enjoy growing popularity among readers. However, the overwhelming number of comments makes it infeasible for the average news consumer to read all of them and hinders engaging discussions. Most platforms display comments in chronological order, which neglects that some of them are more relevant to users and are better conversation starters. In this paper, we systematically analyze user engagement in the form of the upvotes and replies that a comment receives. Based on comment texts, we train a model to distinguish comments that have either a high or low chance of receiving many upvotes and replies. Our evaluation on user comments from TheGuardian.com compares recurrent and convolutional neural network models, and a traditional feature-based classifier. Further, we investigate what makes some comments more engaging than others. To this end, we identify engagement triggers and arrange them in a taxonomy. Explanation methods for neural networks reveal which input words have the strongest influence on our model's predictions. In addition, we evaluate on a dataset of product reviews, which exhibit similar properties as user comments, such as featuring upvotes for helpfulness.


title: Harder-Narasimhan theory

abstract: An abstract formulation of Harder-Narasimhan theory is stated without proof by L. Fargues, and I found it helpful to write it all out.


title: On the Complexity and Approximability of Optimal Sensor Selection and Attack for Kalman Filtering

abstract: Given a linear dynamical system affected by stochastic noise, we consider the problem of selecting an optimal set of sensors (at design-time) to minimize the trace of the steady state a priori or a posteriori error covariance of the Kalman filter, subject to certain selection budget constraints. We show the fundamental result that there is no polynomial-time constant-factor approximation algorithm for this problem. This contrasts with other classes of sensor selection problems studied in the literature, which typically pursue constant-factor approximations by leveraging greedy algorithms and submodularity (or supermodularity) of the cost function. Here, we provide a specific example showing that greedy algorithms can perform arbitrarily poorly for the problem of design-time sensor selection for Kalman filtering. We then study the problem of attacking (i.e., removing) a set of installed sensors, under predefined attack budget constraints, to maximize the trace of the steady state a priori or a posteriori error covariance of the Kalman filter. Again, we show that there is no polynomial-time constant-factor approximation algorithm for this problem, and show specifically that greedy algorithms can perform arbitrarily poorly.


title: Performance of Julia for High Energy Physics Analyses

abstract: We argue that the Julia programming language is a compelling alternative to implementations in Python and C++ for common data analysis workflows in high energy physics. We compare the speed of implementations of different workflows in Julia with those in Python and C++. Our studies show that the Julia implementations are competitive for tasks that are dominated by computational load rather than data access. For work that is dominated by data access, we demonstrate an application with concurrent file reading and parallel data processing.


title: Revisiting type-2 triangular norms on normal convex fuzzy truth values

abstract: This paper studies t-norms on the space $\mathbf{L}$ of all normal and convex fuzzy truth values. We first prove that the only non-convolution form type-2 t-norm constructed by Wu et al. satisfies the distributivity law for meet-convolution and show that t-norm in the sense of Walker and Walker is strictly stronger than t$_{r}$-norm on $\mathbf{L}$, which is strictly stronger than t-norm on $\mathbf{L}$. Furthermore, we characterize some restrictive axioms of t$_{r}$-norms for convolution operations on $\mathbf{L}$ and obtain some necessary conditions for t$_{r}$-(co)norm convolution operations on $\mathbf{L}$ .


title: Bounded State Estimation over Finite-State Channels: Relating Topological Entropy and Zero-Error Capacity

abstract: We investigate state estimation of linear systems over channels having a finite state not known by the transmitter or receiver. We show that similar to memoryless channels, zero-error capacity is the right figure of merit for achieving bounded estimation errors. We then consider finite-state, worst-case versions of the common erasure and additive noise channels models, in which the noise is governed by a finite-state machine without any statistical structure. Upper and lower bounds on their zero-error capacities are derived, revealing a connection with the {\em topological entropy} of the channel dynamics. Separate necessary and sufficient conditions for bounded linear state estimation errors via such channels are obtained. These estimation conditions bring together the topological entropies of the linear system and the discrete channel.


title: Local maximizers of adjoint Fourier restriction estimates for the cone, paraboloid and sphere

abstract: We show that, possibly after a compactification of spacetime, constant functions are local maximizers of the Tomas-Stein adjoint Fourier restriction inequality for the cone and paraboloid in every dimension, and for the sphere in dimension up to 60. For the cone and paraboloid we work from the PDE framework, which enables the use of the Penrose and the Lens transformations, which map the conjectured optimal functions into constants.


title: Search for long-lived, massive particles in events with a displaced vertex and a muon with large impact parameter in $pp$ collisions at $\sqrt{s} = 13$ TeV with the ATLAS detector

abstract: A search for long-lived particles decaying into hadrons and at least one muon is presented. The analysis selects events that pass a muon or missing-transverse-momentum trigger and contain a displaced muon track and a displaced vertex. The analyzed dataset of proton-proton collisions at $\sqrt{s} = 13$ TeV was collected with the ATLAS detector and corresponds to 136 fb$^{-1}$. The search employs dedicated reconstruction techniques that significantly increase the sensitivity to long-lived particle decays that occur in the ATLAS inner detector. Background estimates for Standard Model processes and instrumental effects are extracted from data. The observed event yields are compatible with those expected from background processes. The results are presented as limits at 95% confidence level on model-independent cross sections for processes beyond the Standard Model, and interpreted as exclusion limits in scenarios with pair-production of long-lived top squarks that decay via a small $R$-parity-violating coupling into a quark and a muon. Top squarks with masses up to 1.7 TeV are excluded for a lifetime of 0.1 ns, and masses below 1.3 TeV are excluded for lifetimes between 0.01 ns and 30 ns.


title: Large Nernst Power Factor over a Broad Temperature Range in Polycrystalline Weyl Semimetal NbP

abstract: The discovery of topological materials has provided new opportunities to exploit advanced materials for heat-to-electricity energy conversion as they share many common characteristics with thermoelectric materials. In this work, we report the magneto-thermoelectric properties and Nernst effect of the topological Weyl semimetal NbP. We find that polycrystalline, bulk NbP shows a significantly larger Nernst thermopower than its conventional thermopower under magnetic field. As a result, a maximum Nernst power factor of ~ 35*10-4 Wm-1K-2 is achieved at 9 T and 136 K, which is 4 times higher than its conventional power factor and is also comparable to that of state-of-the-art thermoelectrics. Moreover, the Nernst power factor maintains relatively large value over a broad temperature range. These results highlight that the enhancement of thermoelectric performance can be achieved in topological semimetals based on the Nernst effect and transverse transport.


title: StrokeCoder: Path-Based Image Generation from Single Examples using Transformers

abstract: This paper demonstrates how a Transformer Neural Network can be used to learn a Generative Model from a single path-based example image. We further show how a data set can be generated from the example image and how the model can be used to generate a large set of deviated images, which still represent the original image's style and concept.


title: Pedestrian Models for Autonomous Driving Part II: High-Level Models of Human Behavior

abstract: Autonomous vehicles (AVs) must share space with pedestrians, both in carriageway cases such as cars at pedestrian crossings and off-carriageway cases such as delivery vehicles navigating through crowds on pedestrianized high-streets. Unlike static obstacles, pedestrians are active agents with complex, interactive motions. Planning AV actions in the presence of pedestrians thus requires modelling of their probable future behaviour as well as detecting and tracking them. This narrative review article is Part II of a pair, together surveying the current technology stack involved in this process, organising recent research into a hierarchical taxonomy ranging from low-level image detection to high-level psychological models, from the perspective of an AV designer. This self-contained Part II covers the higher levels of this stack, consisting of models of pedestrian behaviour, from prediction of individual pedestrians' likely destinations and paths, to game-theoretic models of interactions between pedestrians and autonomous vehicles. This survey clearly shows that, although there are good models for optimal walking behaviour, high-level psychological and social modelling of pedestrian behaviour still remains an open research question that requires many conceptual issues to be clarified. Early work has been done on descriptive and qualitative models of behaviour, but much work is still needed to translate them into quantitative algorithms for practical AV control.


title: Measurements of the production cross-section for a $Z$ boson in association with $b$-jets in proton-proton collisions at $\sqrt{s} = 13$ TeV with the ATLAS detector

abstract: This paper presents a measurement of the production cross-section of a $Z$ boson in association with $b$-jets, in proton-proton collisions at $\sqrt{s} = 13$ TeV with the ATLAS experiment at the Large Hadron Collider using data corresponding to an integrated luminosity of 35.6 fb$^{-1}$. Inclusive and differential cross-sections are measured for events containing a $Z$ boson decaying into electrons or muons and produced in association with at least one or at least two $b$-jets with transverse momentum $p_\textrm{T}>$ 20 GeV and rapidity $|y| < 2.5$. Predictions from several Monte Carlo generators based on leading-order (LO) or next-to-leading-order (NLO) matrix elements interfaced with a parton-shower simulation and testing different flavour schemes for the choice of initial-state partons are compared with measured cross-sections. The 5-flavour number scheme predictions at NLO accuracy agree better with data than 4-flavour number scheme ones. The 4-flavour number scheme predictions underestimate data in events with at least one b-jet.


title: PROSPECT-PRO for estimating content of nitrogen-containing leaf proteins and other carbon-based constituents

abstract: Models of radiative transfer (RT) are important tools for remote sensing of vegetation, as they facilitate forward simulations of remotely sensed data as well as inverse estimation of biophysical and biochemical properties from vegetation optical properties. The remote sensing estimation of foliar protein content is a key to monitoring the nitrogen cycle in terrestrial ecosystems in particular to better understand photosynthetic capacity of plants and improve nitrogen management in agriculture. However, no physically based leaf RT model currently allows for proper decomposition of leaf dry matter into nitrogen-based proteins and carbon-based constituents (CBC), estimated from optical properties of fresh or dry foliage. We developed a new version of the PROSPECT model, named PROSPECT-PRO, which separates nitrogen-based constituents (proteins) from CBC (including cellulose, lignin, hemicellulose and starch). PROSPECT-PRO was calibrated and validated on subsets of the LOPEX dataset, accounting for both fresh and dry broadleaf and grass samples. We applied an iterative model inversion optimization algorithm to identify optimal spectral subdomains for retrieval of leaf protein and CBC contents, with 2125-2174 nm optimal for proteins and 2025-2349 nm optimal for CBCs. PROSPECT-PRO inversions revealed a better performance in estimating proteins from optical properties of fresh than dry leaves. We further tested the ability of PROSPECT-PRO to estimate leaf mass per area (LMA) as the sum of proteins and CBC using independent datasets acquired for numerous plant species. Results showed that PROSPECT-PRO is fully compatible and comparable with its predecessor PROSPECT-D in indirect estimation of LMA. We can conclude from findings of this study that PROSPECT-PRO has a high potential in establishing the carbon-to-nitrogen ratio based on the retrieved CBC-to-proteins ratio.


title: A micro-macro Markov chain Monte Carlo method for molecular dynamics using reaction coordinate proposals II: indirect reconstruction

abstract: We introduce a new micro-macro Markov chain Monte Carlo method (mM-MCMC) with indirect reconstruction to sample invariant distributions of molecular dynamics systems that exhibit a time-scale separation between the microscopic (fast) dynamics, and the macroscopic (slow) dynamics of some low-dimensional set of reaction coordinates. The algorithm enhances exploration of the state space in the presence of metastability by allowing larger proposal moves at the macroscopic level, on which a conditional accept-reject procedure is applied. Only when the macroscopic proposal is accepted, the full microscopic state is reconstructed from the newly sampled reaction coordinate value and is subjected to a second accept/reject procedure. The computational gain stems from the fact that most proposals are rejected at the macroscopic level, at low computational cost, while microscopic states, once reconstructed, are almost always accepted. This paper discusses an indirect method to reconstruct microscopic samples from macroscopic reaction coordinate values, that can also be applied in cases where direct reconstruction is cumbersome. The indirect reconstruction method generates a microscopic sample by performing a biased microscopic simulation, starting from the previous microscopic sample and driving the microscopic state towards the proposed reaction coordinate value. We show numerically that the mM-MCMC scheme with indirect reconstruction can significantly extend the range of applicability of the mM-MCMC method.


title: TLDR: Token Loss Dynamic Reweighting for Reducing Repetitive Utterance Generation

abstract: Natural Language Generation (NLG) models are prone to generating repetitive utterances. In this work, we study the repetition problem for encoder-decoder models, using both recurrent neural network (RNN) and transformer architectures. To this end, we consider the chit-chat task, where the problem is more prominent than in other tasks that need encoder-decoder architectures. We first study the influence of model architectures. By using pre-attention and highway connections for RNNs, we manage to achieve lower repetition rates. However, this method does not generalize to other models such as transformers. We hypothesize that the deeper reason is that in the training corpora, there are hard tokens that are more difficult for a generative model to learn than others and, once learning has finished, hard tokens are still under-learned, so that repetitive generations are more likely to happen. Based on this hypothesis, we propose token loss dynamic reweighting (TLDR) that applies differentiable weights to individual token losses. By using higher weights for hard tokens and lower weights for easy tokens, NLG models are able to learn individual tokens at different paces. Experiments on chit-chat benchmark datasets show that TLDR is more effective in repetition reduction for both RNN and transformer architectures than baselines using different weighting functions.


title: Rigorous State Evolution Analysis for Approximate Message Passing with Side Information

abstract: A common goal in many research areas is to reconstruct an unknown signal x from noisy linear measurements. Approximate message passing (AMP) is a class of low-complexity algorithms that can be used for efficiently solving such high-dimensional regression tasks. Often, it is the case that side information (SI) is available during reconstruction. For this reason, a novel algorithmic framework that incorporates SI into AMP, referred to as approximate message passing with side information (AMP-SI), has been recently introduced. In this work, we provide rigorous performance guarantees for AMP-SI when there are statistical dependencies between the signal and SI pairs and the entries of the measurement matrix are independent and identically distributed Gaussian. The AMP-SI performance is shown to be provably tracked by a scalar iteration referred to as state evolution. Moreover, we provide numerical examples that demonstrate empirically that the SE can predict the AMP-SI mean square error accurately.


title: Surface homogenization of an array of Helmholtz resonators for a viscoacoustic model using two-scale convergence

abstract: We derive the weak limit of a linear viscoacoustic model in an acoustic liner that is a chamber connected to a periodic repetition of elongated chambers -- the Helmholtz resonators. As model we consider the time-harmonic and linearized compressible Navier-Stokes equations for the acoustic velocity and pressure. Following the approach in Schmidt et al., J. Math. Ind 8:15, 2018 for the viscoacoustic transmission problem of multiperforated plates the viscosity is scaled as $\delta^{4}$ with the period $\delta$ of the array of chambers and the size of the necks as well as the wall thickness like $\delta^2$ such that the viscous boundary layers are of the order of the size of the necks. Applying the method of two-scale convergence we obtain with a stability assumption in the limit $\delta \to 0$ that the acoustic pressure fulfills the Helmholz equation with impedance boundary conditions. These boundary conditions depend on the frequency, the length of the resonators and through the effective Rayleigh conductivity -- that can be computed numerically -- on the shape of their necks. We compare the limit model to semi-analytical models in the literature.


title: Two-Dimensional Materials for Energy-Efficient Spin-Orbit Torque Devices

abstract: Spin-orbit torques (SOTs), which rely on spin current generation from charge current in a nonmagnetic material, promise an energy-efficient scheme for manipulating magnetization in magnetic devices. A critical topic for spintronic devices using SOTs is to enhance the charge to spin conversion efficiency. Besides, the current-induced spin polarization is usually limited to in-plane, whereas out-of-plane spin polarization could be favored for efficient perpendicular magnetization switching. Recent advances in utilizing two important classes of van der Waals materials$-$topological insulators and transition-metal dichalcogenides$-$as spin sources to generate SOT shed light on addressing these challenges. Topological insulators such as bismuth selenide have shown a giant SOT efficiency, which is larger than those from three-dimensional heavy metals by at least one order of magnitude. Transition-metal dichalcogenides such as tungsten telluride have shown a current-induced out-of-plane spin polarization, which is allowed by the reduced symmetry. In this review, we use symmetry arguments to predict and analyze SOTs in van der Waal material-based heterostructures. We summarize the recent progress of SOT studies based on topological insulators and transition-metal dichalcogenides and show how these results are in line with the symmetry arguments. At last, we identify unsolved issues in the current studies and suggest three potential research directions in this field.


title: XBlock-EOS: Extracting and Exploring Blockchain Data From EOSIO

abstract: Blockchain-based cryptocurrencies and applications have flourished in blockchain research community. Massive data generated from diverse blockchain systems bring not only huge business values but also technological challenges in data analytics of heterogeneous blockchain data. Different from Bitcoin and Ethereum, EOSIO has richer diversity and a higher volume of blockchain data due to its unique architectural design in resource management, consensus scheme and high throughput. Despite its popularity (e.g., 89,800,000 blocks generated till November 14, 2019 since its launch on June 8, 2018), few studies have been made on data analysis of EOSIO. To fill this gap, we collect and process the up-to-date on-chain data from EOSIO. We name these well-processed EOSIO datasets as XBlock-EOS, which consists of 7 well-processed datasets: 1) Block, Transaction and Action, 2) Internal and External EOS Transfer Action, 3) Contract Information, 4) Contract Invocation, 5) Token Action, 6) Account Creation, 7) Resource Management. It is challenging to process and analyze a high volume of raw EOSIO data and establish the mapping from original raw data to the well-grained datasets since it requires substantial efforts in extracting various types of data as well as sophisticated knowledge on software engineering and data analytics. Meanwhile, we present statistics and exploration on these datasets. Moreover, we also outline the possible research opportunities based on XBlock-EOS.


title: Efficient Calculation of Excitonic Effects in Solids Including Approximated Quasiparticle Energies

abstract: In this work we present a new procedure to compute optical spectra including excitonic effects and approximated quasiparticle corrections with reduced computational effort. The excitonic effects on optical spectra are included by solving the Bethe-Salpeter equation, considering quasiparticle eigenenergies and respective wavefunctions obtained within DFT-1/2 method. The electron-hole ladder diagrams are approximated by the screened exchange. To prove the capability of the procedure, we compare the calculated imaginary part of the dielectric functions of Si, Ge, GaAs, GaP, GaSb, InAs, InP, and InSb with experimental data. The energy position of the absorption peaks are correctly described. The good agreement with experimental results together with the very significant reduction of computational effort makes the procedure suitable on the investigation of optical spectra of more complex systems.


title: Ultrafast strain-induced charge transport in semiconductor superlattices

abstract: We investigate the effect of hypersonic (> 1 GHz) acoustic phonon wavepackets on electron transport in a semiconductor superlattice. Our quantum mechanical simulations demonstrate that a GHz train of picosecond deformation strain pulses propagating through a superlattice can generate current oscillations whose frequency is several times higher than that of the strain pulse train. The shape and polarity of the calculated current pulses agree well with experimentally measured electric signals. The calculations also explain and accurately reproduce the measured variation of the induced current pulse magnitude with the strain pulse amplitude and applied bias voltage. Our results open a route to developing acoustically-driven semiconductor superlattices as sources of millimetre and sub-millimetre electromagnetic waves.


title: Efficient search of optimal Flower Constellations

abstract: We derive an analytical closed expression to compute the minimum distance (quantified by the angle of separation measured from the center of the Earth) between any two satellites located at the same altitude and in circular orbits. We also exploit several properties of Flower Constellations (FCs) that, combined with our formula for the distance, give an efficient method to compute the minimum angular distance between satellites, for all possible FCs with up to a given number of satellites.


title: Enabling Pulse-level Programming, Compilation, and Execution in XACC

abstract: Noisy gate-model quantum processing units (QPUs) are currently available from vendors over the cloud, and digital quantum programming approaches exist to run low-depth circuits on physical hardware. These digital representations are ultimately lowered to pulse-level instructions by vendor quantum control systems to affect unitary evolution representative of the submitted digital circuit. Vendors are beginning to open this pulse-level control system to the public via specified interfaces. Robust programming methodologies, software frameworks, and backend simulation technologies for this analog model of quantum computation will prove critical to advancing pulse-level control research and development. Prototypical use cases for this include error mitigation, optimal pulse control, and physics-inspired pulse construction. Here we present an extension to the XACC quantum-classical software framework that enables pulse-level programming for superconducting, gate-model quantum computers, and a novel, general, and extensible pulse-level simulation backend for XACC that scales on classical compute clusters via MPI. Our work enables custom backend Hamiltonian definitions and gate-level compilation to available pulses with a focus on performance and scalability. We end with a demonstration of this capability, and show how to use XACC for pertinent pulse-level programming tasks.


title: Hybrid Precoding For Millimeter Wave MIMO Systems: A Matrix Factorization Approach

abstract: This paper investigates the hybrid precoding design for millimeter wave (mmWave) multiple-input multiple-output (MIMO) systems with finite-alphabet inputs. The precoding problem is a joint optimization of analog and digital precoders, and we treat it as a matrix factorization problem with power and constant modulus constraints. Our work presents three main contributions: First, we present a sufficient condition and a necessary condition for hybrid precoding schemes to realize unconstrained optimal precoders exactly when the number of data streams Ns satisfies Ns = minfrank(H);Nrfg, where H represents the channel matrix and Nrf is the number of radio frequency (RF) chains. Second, we show that the coupled power constraint in our matrix factorization problem can be removed without loss of optimality. Third, we propose a Broyden-Fletcher-Goldfarb-Shanno (BFGS)-based algorithm to solve our matrix factorization problem using gradient and Hessian information. Several numerical results are provided to show that our proposed algorithm outperforms existing hybrid precoding algorithms.


title: GISNet: Graph-Based Information Sharing Network For Vehicle Trajectory Prediction

abstract: The trajectory prediction is a critical and challenging problem in the design of an autonomous driving system. Many AI-oriented companies, such as Google Waymo, Uber and DiDi, are investigating more accurate vehicle trajectory prediction algorithms. However, the prediction performance is governed by lots of entangled factors, such as the stochastic behaviors of surrounding vehicles, historical information of self-trajectory, and relative positions of neighbors, etc. In this paper, we propose a novel graph-based information sharing network (GISNet) that allows the information sharing between the target vehicle and its surrounding vehicles. Meanwhile, the model encodes the historical trajectory information of all the vehicles in the scene. Experiments are carried out on the public NGSIM US-101 and I-80 Dataset and the prediction performance is measured by the Root Mean Square Error (RMSE). The quantitative and qualitative experimental results show that our model significantly improves the trajectory prediction accuracy, by up to 50.00%, compared to existing models.


title: A Topological Characterization of Modulo-$p$ Arguments and Implications for Necklace Splitting

abstract: The classes PPA-$p$ have attracted attention lately, because they are the main candidates for capturing the complexity of Necklace Splitting with $p$ thieves, for prime $p$. However, these classes were not known to have complete problems of a topological nature, which impedes any progress towards settling the complexity of the Necklace Splitting problem. On the contrary, topological problems have been pivotal in obtaining completeness results for PPAD and PPA, such as the PPAD-completeness of finding a Nash equilibrium [Daskalakis et al., 2009, Chen et al., 2009b] and the PPA-completeness of Necklace Splitting with 2 thieves [Filos-Ratsikas and Goldberg, 2019]. In this paper, we provide the first topological characterization of the classes PPA-$p$. First, we show that the computational problem associated with a simple generalization of Tucker's Lemma, termed $p$-polygon-Tucker, as well as the associated Borsuk-Ulam-type theorem, $p$-polygon-Borsuk-Ulam, are PPA-$p$-complete. Then, we show that the computational version of the well-known BSS Theorem [Barany et al., 1981], as well as the associated BSS-Tucker problem are PPA-$p$-complete. Finally, using a different generalization of Tucker's Lemma (termed $\mathbb{Z}_p$-star-Tucker), which we prove to be PPA-$p$-complete, we prove that $p$-thief Necklace Splitting is in PPA-$p$. This latter result gives a new combinatorial proof for the Necklace Splitting theorem, the only proof of this nature other than that of Meunier [2014]. All of our containment results are obtained through a new combinatorial proof for $\mathbb{Z}_p$-versions of Tucker's lemma that is a natural generalization of the standard combinatorial proof of Tucker's lemma by Freund and Todd [1981]. We believe that this new proof technique is of independent interest.


title: Purim: a rapid method with reduced cost for massive detection of CoVid-19

abstract: The CoVid-19 is spreading pandemically all over the world. A rapid defeat of the pandemic requires carrying out on the population a mass screening, able to separate positive from negative cases. Such a cleaning will free a flow of productive population. The current rate and cost of testing, performed with the common PCR (polymerase chain reaction) method and with the available resources, is forcing a selection of the subjects to be tested. Indeed, each one must be examined individually at the cost of precious time. Moreover, the exclusion of potentially positive individuals from screening induces health risks, a broad slowdown in the effort to curb the viral spread, and the consequent mortality rates. We present a new procedure, the Purified by Unified Resampling of Infected Multitudes, in short Purim, able to untangle any massive candidate sample with inexpensive screening, through the cross-correlated analysis of the joint speciments. This procedure can reveal and detect most negative patients and in most cases discover the identity of the few positives already in the first or few secondary tests. We investigate the the two-dimensional correlation case in function of the infection probability. The multi-dimensional topology, the scaled Purim procedure are also considered. Extensive Purim tests may measure and weight the degree of epidemic: their outcome may identify focal regions in the early stages. Assuming hundreds or thousand subjects, the saving both in time and in cost will be remarkable. Purim may be able to filter scheduled flights, scholar acceptance, popular international event participants. The optimal extension of Purim outcome is growing as the inverse of the epidemia expansion. Therefore, the earlier, the better.


title: Does a solar filament barb always correspond to a prominence foot?

abstract: Solar filaments are dark structures on the solar disk, with an elongated spine and several barbs extending out from the spine. When appearing above the solar limb, a filament is called a prominence, with several feet extending down to the solar surface. It was generally thought that filament barbs are simply the prominence feet veering away from the spine and down to the solar surface. However, it was recently noticed that there might be another dynamic type of barbs, which were proposed to be due to filament thread longitudinal oscillation. If this is the case, the dynamic barbs would not extend down to the solar surface. With the quadrature observations of a filament barb on 2011 June 5 from the {\it Solar Dynamics Observatory} and the {\it STEREO} satellites, we confirm that the filament barb is due to filament thread longitudinal oscillations. Viewed from the side, the filament barb looks like an appendix along the spine of the prominence, and does not extend down to the solar surface as a foot.


title: High frequency mechanical excitation of a silicon nanostring with piezoelectric aluminum nitride layers

abstract: A strong trend for quantum based technologies and applications follows the avenue of combining different platforms to exploit their complementary technological and functional advantages. Micro and nano-mechanical devices are particularly suitable for hybrid integration due to the easiness of fabrication at multi-scales and their pervasive coupling with electrons and photons. Here, we report on a nanomechanical technological platform where a silicon chip is combined with an aluminum nitride layer. Exploiting the AlN piezoelectricity, Surface Acoustic Waves are injected in the Si layer where the material has been localy patterned and etched to form a suspended nanostring. Characterizing the nanostring vertical displacement induced by the SAW, we found an external excitation peak efficiency in excess of 500 pm/V at 1 GHz mechanical frequency. Exploiting the long term expertise in silicon photonic and electronic devices as well as the SAW robustness and versatility, our technological platform represents a strong candidate for hybrid quantum systems.


title: Higher derivative scalar-tensor monomials and their classification

abstract: We make a full classification of scalar monomials built of the Riemann curvature tensor up to the quadratic order and of the covariant derivatives of the scalar field up to the third order. From the point of view of the effective field theory, the third or even higher order covariant derivatives of the scalar field are of the same order as the higher curvature terms, and thus should be taken into account. Moreover, higher curvature terms and higher order derivatives of the scalar field are complementary to each other, of which novel ghost-free combinations may exist. We make a systematic classification of all the possible monomials, according to the numbers of Riemann tensor and higher derivatives of the scalar field in each monomial. Complete basis of monomials at each order are derived, of which linear combinations may yield novel ghost-free Lagrangians. We also develop a diagrammatic representation of the monomials, which may help to simplify the analysis.


title: Minimising Good-for-Games automata is NP complete

abstract: This paper discusses the hardness of finding minimal good-for-games (GFG) Buchi, Co-Buchi, and parity automata with state based acceptance. The problem appears to sit between finding small deterministic and finding small nondeterministic automata, where minimality is NP-complete and PSPACE-complete, respectively. However, recent work of Radi and Kupferman has shown that minimising Co-Buchi automata with transition based acceptance is tractable, which suggests that the complexity of minimising GFG automata might be cheaper than minimising deterministic automata. We show for the standard state based acceptance that the minimality of a GFG automaton is NP-complete for Buchi, Co-Buchi, and parity GFG automata. The proofs are a surprisingly straight forward generalisation of the proofs from deterministic Buchi automata: they use a similar reductions, and the same hard class of languages.


title: The Planetary Luminosity Problem: "Missing Planets" and the Observational Consequences of Episodic Accretion

abstract: The high occurrence rates of spiral arms and large central clearings in protoplanetary disks, if interpreted as signposts of giant planets, indicate that gas giants form commonly as companions to young stars ($<$ few Myr) at orbital separations of 10--300\,au. However, attempts to directly image this giant planet population as companions to more mature stars ($> 10$\, Myr) have yielded few successes. This discrepancy could be explained if most giant planets form "cold start," i.e., by radiating away much of their formation energy as they assemble their mass, rendering them faint enough to elude detection at later times. In that case, giant planets should be bright at early times, during their accretion phase, and yet forming planets are detected only rarely through direct imaging techniques. Here we explore the possibility that the low detection rate of accreting planets is the result of episodic accretion through a circumplanetary disk. We also explore the possibility that the companion orbiting the Herbig Ae star HD~142527 may be a giant planet undergoing such an accretion outburst.


title: Tomography of time-bin quantum states using time-resolved detection

abstract: We present a method for measuring quantum states encoded in the temporal modes of photons. The basis for the multilevel quantum states is defined by the use of modes propagating in a dispersive medium, which is a fiber in this case. The propagation and time-resolved single photon detection allow us to define a positive-operator valued measure (POVM). The POVM depends on the amount of dispersion and the characteristics of a detector. This framework is numerically tested by performing quantum state tomography on a large number of states for a set of realistic experimental settings. Finally, the average fidelity between the expected and reconstructed states is computed for qubits, qutrits and entangled qubits.


title: In defence of metric learning for speaker recognition

abstract: The objective of this paper is 'open-set' speaker recognition of unseen speakers, where ideal embeddings should be able to condense information into a compact utterance-level representation that has small intra-speaker and large inter-speaker distance. A popular belief in speaker recognition is that networks trained with classification objectives outperform metric learning methods. In this paper, we present an extensive evaluation of most popular loss functions for speaker recognition on the VoxCeleb dataset. We demonstrate that the vanilla triplet loss shows competitive performance compared to classification-based losses, and those trained with our proposed metric learning objective outperform state-of-the-art methods.


title: A vignetting advantage for thin-film filter arrays in hyperspectral cameras

abstract: Vignetting in camera lenses is generally seen as something to avoid. For spectral cameras with thin-film interference filters, however, we argue that vignetting can be an advantage. When illuminated by focused light, the bandwidth of interference filters increases with the chief-ray angle, causing position-dependent smoothing of the spectra. We show that vignetting can be used to reduce smoothing and preserve important spectral features. Furthermore, we demonstrate that by adding additional vignetting to a lens, measurements can be made more consistent across the scene. This makes vignetting a useful parameter during spectral camera design.


title: Completeness in affine and statistical geometry

abstract: We begin the study of completeness of affine connections, especially those on statistical manifolds as well as on affine hypersurfaces. We collect basic facts, prove new theorems and provide examples with remarkable properties.


title: Is the Juice Worth the Squeeze? Machine Learning (ML) In and For Agent-Based Modelling (ABM)

abstract: In recent years, many scholars praised the seemingly endless possibilities of using machine learning (ML) techniques in and for agent-based simulation models (ABM). To get a more comprehensive understanding of these possibilities, we conduct a systematic literature review (SLR) and classify the literature on the application of ML in and for ABM according to a theoretically derived classification scheme. We do so to investigate how exactly machine learning has been utilized in and for agent-based models so far and to critically discuss the combination of these two promising methods. We find that, indeed, there is a broad range of possible applications of ML to support and complement ABMs in many different ways, already applied in many different disciplines. We see that, so far, ML is mainly used in ABM for two broad cases: First, the modelling of adaptive agents equipped with experience learning and, second, the analysis of outcomes produced by a given ABM. While these are the most frequent, there also exist a variety of many more interesting applications. This being the case, researchers should dive deeper into the analysis of when and how which kinds of ML techniques can support ABM, e.g. by conducting a more in-depth analysis and comparison of different use cases. Nonetheless, as the application of ML in and for ABM comes at certain costs, researchers should not use ML for ABMs just for the sake of doing it.


title: Demonstration of Muon-Beam Transverse Phase-Space Compression

abstract: We demonstrate efficient transverse compression of a 12.5 MeV/c muon beam stopped in a helium gas target featuring a vertical density gradient and crossed electric and magnetic fields. The muon stop distribution extending vertically over 14 mm was reduced to a 0.25 mm size (RMS) within 3.5 $\mu$s. The simulation including cross sections for low-energy $\mu^+$-$\text{He}$ elastic and charge exchange ($\mu^+\leftrightarrow $ muonium) collisions describes the measurements well. By combining the transverse compression stage with a previously demonstrated longitudinal compression stage, we can improve the phase space density of a $\mu^+ $ beam by a factor of $ 10^{10} $ with $ 10^{-3} $ efficiency.


title: Partially Observed Discrete-Time Risk-Sensitive Mean Field Games

abstract: In this paper, we consider discrete-time partially observed mean-field games with the risk-sensitive optimality criterion. We introduce risk-sensitivity behaviour for each agent via an exponential utility function. In the game model, each agent is weakly coupled with the rest of the population through its individual cost and state dynamics via the empirical distribution of states. We establish the mean-field equilibrium in the infinite-population limit using the technique of converting the underlying original partially observed stochastic control problem to a fully observed one on the belief space and the dynamic programming principle. Then, we show that the mean-field equilibrium policy, when adopted by each agent, forms an approximate Nash equilibrium for games with sufficiently many agents. We first consider finite-horizon cost function, and then, discuss extension of the result to infinite-horizon cost in the next-to-last section of the paper.


title: Severity Assessment of Coronavirus Disease 2019 (COVID-19) Using Quantitative Features from Chest CT Images

abstract: Background: Chest computed tomography (CT) is recognized as an important tool for COVID-19 severity assessment. As the number of affected patients increase rapidly, manual severity assessment becomes a labor-intensive task, and may lead to delayed treatment. Purpose: Using machine learning method to realize automatic severity assessment (non-severe or severe) of COVID-19 based on chest CT images, and to explore the severity-related features from the resulting assessment model. Materials and Method: Chest CT images of 176 patients (age 45.3$\pm$16.5 years, 96 male and 80 female) with confirmed COVID-19 are used, from which 63 quantitative features, e.g., the infection volume/ratio of the whole lung and the volume of ground-glass opacity (GGO) regions, are calculated. A random forest (RF) model is trained to assess the severity (non-severe or severe) based on quantitative features. Importance of each quantitative feature, which reflects the correlation to the severity of COVID-19, is calculated from the RF model. Results: Using three-fold cross validation, the RF model shows promising results, i.e., 0.933 of true positive rate, 0.745 of true negative rate, 0.875 of accuracy, and 0.91 of area under receiver operating characteristic curve (AUC). The resulting importance of quantitative features shows that the volume and its ratio (with respect to the whole lung volume) of ground glass opacity (GGO) regions are highly related to the severity of COVID-19, and the quantitative features calculated from the right lung are more related to the severity assessment than those of the left lung. Conclusion: The RF based model can achieve automatic severity assessment (non-severe or severe) of COVID-19 infection, and the performance is promising. Several quantitative features, which have the potential to reflect the severity of COVID-19, were revealed.


title: Superdeterministic hidden-variables models I: nonequilibrium and signalling

abstract: This is the first of two papers which attempt to comprehensively analyse superdeterministic hidden-variables models of Bell correlations. We first give an overview of superdeterminism and discuss various criticisms of it raised in the literature. We argue that the most common criticism, the violation of `free-will', is incorrect. We take up Bell's intuitive criticism that these models are `conspiratorial'. To develop this further, we introduce nonequilibrium extensions of superdeterministic models. We show that the measurement statistics of these extended models depend on the physical system used to determine the measurement settings. This suggests a fine-tuning in order to eliminate this dependence from experimental observation. We also study the signalling properties of these extended models. We show that although they generally violate the formal no-signalling constraints, this violation cannot be equated to an actual signal. We therefore suggest that the so-called no-signalling constraints be more appropriately named the marginal-independence constraints. We discuss the mechanism by which marginal-independence is violated in superdeterministic models. Lastly, we consider a hypothetical scenario where two experimenters use the apparent-signalling of a superdeterministic model to communicate with each other. This scenario suggests another conspiratorial feature peculiar to superdeterminism. These suggestions are quantitatively developed in the second paper.


title: Supervisory model predictive control for PV battery and heat pump system with phase change slurry thermal storage

abstract: We present the design, implementation and experimental validation of a supervisory predictive control approach for an electrical heating system featuring a phase change slurry as heat storage and transfer medium. The controller optimizes the energy flows that are used as set points for the heat generation and energy distribution components. The optimization handles the thermal and electrical subsystems simultaneously and is able to switch between different objectives. We show the control can be implemented on low-cost embedded hardware and validate it with an experimental test bed comprising an installation of the complete heating system, including all hydraulic and all electrical components. Experimental results demonstrate the feasibility of both, a heat pump heating system with a phase change slurry, and the optimal control approach. The main control objectives, i.e., thermal comfort and maximum self-consumption of solar energy, can be met. In addition, the system and its controller provide a load shifting potential.


title: Estimating Treatment Effects with Observed Confounders and Mediators

abstract: Given a causal graph, the do-calculus can express treatment effects as functionals of the observational joint distribution that can be estimated empirically. Sometimes the do-calculus identifies multiple valid formulae, prompting us to compare the statistical properties of the corresponding estimators. For example, the backdoor formula applies when all confounders are observed and the frontdoor formula applies when an observed mediator transmits the causal effect. In this paper, we investigate the over-identified scenario where both confounders and mediators are observed, rendering both estimators valid. Addressing the linear Gaussian causal model, we demonstrate that either estimator can dominate the other by an unbounded constant factor. Next, we derive an optimal estimator, which leverages all observed variables, and bound its finite-sample variance. We show that it strictly outperforms the backdoor and frontdoor estimators and that this improvement can be unbounded. We also present a procedure for combining two datasets, one with observed confounders and another with observed mediators. Finally, we evaluate our methods on both simulated data and the IHDP and JTPA datasets.


title: Georg Mohr's "Euclides Danicus" -- Preliminary Version

abstract: We present here a preliminary version of a translation with comments of Georg Mohr's book "Euclides Danicus", where the first proof of Mohr-Mascheroni Theorem appeared in 1672, 125 years before Mascheroni's book.


title: Several extremal problems on graphs involving the circumference, girth, and hyperbolicity constant

abstract: To compute the hyperbolicity constant is an almost intractable problem, thus it is natural to try to bound it in terms of some parameters of the graph. Let $\mathcal{G}(g,c,n)$ be the set of graphs $G$ with girth $g(G)=g$, circumference $c(G)=c$, and $n$ vertices; and let $\mathcal{H}(g,c,m)$ be the set of graphs with girth $g$, circumference $c$, and $m$ edges. In this work, we study the four following extremal problems on graphs: $A(g,c,n)=\min\{\delta(G)\,|\; G \in \mathcal{G}(g,c,n) \}$, $B(g,c,n)=\max\{\delta(G)\,|\; G \in \mathcal{G}(g,c,n) \}$, $\alpha(g,c,m)=\min\{\delta(G)\,|\; \in \mathcal{H}(g,c,m) \}$ and $\beta(g,c,m)=\max\{\delta(G)\,|\; G \in \mathcal{H}(g,c,m) \}$. In particular, we obtain bounds for $A(g,c,n)$ and $\alpha(g,c,m)$, and we compute the precise value of $B(g,c,n)$ and $\beta(g,c,m)$ for all values of $g$, $c$, $n$ and $m$.


title: Single-Shot 3D Widefield Fluorescence Imaging with a Computational Miniature Mesoscope

abstract: Fluorescence imaging is indispensable to biology and neuroscience. The need for large-scale imaging in freely behaving animals has further driven the development in miniaturized microscopes (miniscopes). However, conventional microscopes / miniscopes are inherently constrained by their limited space-bandwidth-product, shallow depth-of-field, and the inability to resolve 3D distributed emitters. Here, we present a Computational Miniature Mesoscope (CM$^2$) that overcomes these bottlenecks and enables single-shot 3D imaging across an 8 $\times$ 7-mm$^2$ field-of-view and 2.5-mm depth-of-field, achieving 7-$\mu$m lateral resolution and better than 200-$\mu$m axial resolution. Notably, the CM$^2$ has a compact lightweight design that integrates a microlens array for imaging and an LED array for excitation in a single platform. Its expanded imaging capability is enabled by computational imaging that augments the optics by algorithms. We experimentally validate the mesoscopic 3D imaging capability on volumetrically distributed fluorescent beads and fibers. We further quantify the effects of bulk scattering and background fluorescence on phantom experiments.


title: Secure Groupcast with Shared Keys

abstract: We consider a transmitter and $K$ receivers, each of which shares a key variable with the transmitter. Through a noiseless broadcast channel, the transmitter wishes to send a common message $W$ securely to $N$ out of the $K$ receivers while the remaining $K-N$ receivers learn no information about $W$. We are interested in the maximum message rate, i.e., the maximum number of bits of $W$ that can be securely groupcast to the legitimate receivers per key block and the minimum broadcast bandwidth, i.e., the minimum number of bits of the broadcast information required to securely groupcast the message bits. We focus on the setting of combinatorial keys, where every subset of the $K$ receivers share an independent key of arbitrary size. Under this combinatorial key setting, the maximum message rate is characterized for the following scenarios - 1) $N=1$ or $N=K-1$, i.e., secure unicast to 1 receiver with $K-1$ eavesdroppers or secure groupcast to $K-1$ receivers with $1$ eavesdropper, 2) $N=2, K=4$, i.e., secure groupcast to $2$ out of 4 receivers, and 3) the symmetric setting where the key size for any subset of the same cardinality is equal for any $N,K$. Further, for the latter two cases, the minimum broadcast bandwidth for the maximum message rate is characterized.


title: Accelerated Analog Neuromorphic Computing

abstract: This paper presents the concepts behind the BrainScales (BSS) accelerated analog neuromorphic computing architecture. It describes the second-generation BrainScales-2 (BSS-2) version and its most recent in-silico realization, the HICANN-X Application Specific Integrated Circuit (ASIC), as it has been developed as part of the neuromorphic computing activities within the European Human Brain Project (HBP). While the first generation is implemented in an 180nm process, the second generation uses 65nm technology. This allows the integration of a digital plasticity processing unit, a highly-parallel micro processor specially built for the computational needs of learning in an accelerated analog neuromorphic systems. The presented architecture is based upon a continuous-time, analog, physical model implementation of neurons and synapses, resembling an analog neuromorphic accelerator attached to build-in digital compute cores. While the analog part emulates the spike-based dynamics of the neural network in continuous-time, the latter simulates biological processes happening on a slower time-scale, like structural and parameter changes. Compared to biological time-scales, the emulation is highly accelerated, i.e. all time-constants are several orders of magnitude smaller than in biology. Programmable ion channel emulation and inter-compartmental conductances allow the modeling of nonlinear dendrites, back-propagating action-potentials as well as NMDA and Calcium plateau potentials. To extend the usability of the analog accelerator, it also supports vector-matrix multiplication. Thereby, BSS-2 supports inference of deep convolutional networks as well as local-learning with complex ensembles of spiking neurons within the same substrate.


title: Short-term predictions of country-specific Covid-19 infection rates based on power law scaling exponents

abstract: The number of corona virus (COVID-19) infections grows worldwide. In order to create short term predictions to prepare for the extent of the global pandemic we analyze infection data from the top 25 affected countries. It is shown that all country-specific infection rates follow a power law growth behavior and the scaling exponents per country are calculated. We find two different growth patterns: either steady power law growth from the very beginning with moderate scaling exponents of 3-5 or explosive power law growth with dramatic scaling exponents of 8-11. In the case of the USA we even find an exponent of 16.59. By means of data analysis we confirm that instituting strict measures of lock-downs combined with a strong adherence by the population are effective means to bring the growth rates down. While many countries have instituted measures there are only three countries (Denmark, Norway, and South Korea) so far where such lock-downs led to a significant reduction of the growth rate. In the case of Denmark we calculate the reduction of the scaling exponents to move from 6.82 to 1.47.


title: A Blind Permutation Similarity Algorithm

abstract: This paper introduces a polynomial blind algorithm that determines when two square matrices, $A$ and $B$, are permutation similar. The shifted and translated matrices $(A+\beta I+\gamma J)$ and $(B+\beta I+\gamma J)$ are used to color the vertices of two square, edge weighted, rook's graphs. Then the orbits are found by repeated symbolic squaring of the vertex colored and edge weighted adjacency matrices. Multisets of the diagonal symbols from non-permutation similar matrices are distinct within a few iterations, typically four or less.