title: Access Structure Hiding Secret Sharing from Novel Set Systems and Vector Families

abstract: Secret sharing provides a means to distribute shares of a secret such that any authorized subset of shares, specified by an access structure, can be pooled together to recompute the secret. The standard secret sharing model requires public access structures, which violates privacy and facilitates the adversary by revealing high-value targets. In this paper, we address this shortcoming by introducing \emph{hidden access structures}, which remain secret until some authorized subset of parties collaborate. The central piece of this work is the construction of a set-system $\mathcal{H}$ with strictly greater than $\exp\left(c \dfrac{1.5 (\log h)^2}{\log \log h}\right)$ subsets of a set of $h$ elements. Our set-system $\mathcal{H}$ is defined over $\mathbb{Z}_m$, where $m$ is a non-prime-power, such that the size of each set in $\mathcal{H}$ is divisible by $m$ but the sizes of their pairwise intersections are not divisible by $m$, unless one set is a subset of another. We derive a vector family $\mathcal{V}$ from $\mathcal{H}$ such that superset-subset relationships in $\mathcal{H}$ are represented by inner products in $\mathcal{V}$. We use $\mathcal{V}$ to "encode" the access structures and thereby develop the first \emph{access structure hiding} secret sharing scheme. For a setting with $\ell$ parties, our scheme supports $2^{\binom{\ell}{\ell/2+1}}$ out of the $2^{2^{\ell - O(\log \ell)}}$ total monotone access structures, and its maximum share size for any access structures is $(1+ o(1)) \dfrac{2^{\ell+1}}{\sqrt{\pi \ell/2}}$. The scheme assumes semi-honest polynomial-time parties, and its security relies on the Generalized Diffie-Hellman assumption.


title: Training Deep Neural Networks Without Batch Normalization

abstract: Training neural networks is an optimization problem, and finding a decent set of parameters through gradient descent can be a difficult task. A host of techniques has been developed to aid this process before and during the training phase. One of the most important and widely used class of method is normalization. It is generally favorable for neurons to receive inputs that are distributed with zero mean and unit variance, so we use statistics about dataset to normalize them before the first layer. However, this property cannot be guaranteed for the intermediate activations inside the network. A widely used method to enforce this property inside the network is batch normalization. It was developed to combat covariate shift inside networks. Empirically it is known to work, but there is a lack of theoretical understanding about its effectiveness and potential drawbacks it might have when used in practice. This work studies batch normalization in detail, while comparing it with other methods such as weight normalization, gradient clipping and dropout. The main purpose of this work is to determine if it is possible to train networks effectively when batch normalization is removed through adaption of the training process.


title: Super-Human Performance in Gran Turismo Sport Using Deep Reinforcement Learning

abstract: Autonomous car racing is a major challenge in robotics. It raises fundamental problems for classical approaches such as planning minimum-time trajectories under uncertain dynamics and controlling the car at the limits of its handling. Besides, the requirement of minimizing the lap time, which is a sparse objective, and the difficulty of collecting training data from human experts have also hindered researchers from directly applying learning-based approaches to solve the problem. In the present work, we propose a learning-based system for autonomous car racing by leveraging a high-fidelity physical car simulation, a course-progress proxy reward, and deep reinforcement learning. We deploy our system in Gran Turismo Sport, a world-leading car simulator known for its realistic physics simulation of different race cars and tracks, which is even used to recruit human race car drivers. Our trained policy achieves autonomous racing performance that goes beyond what had been achieved so far by the built-in AI, and, at the same time, outperforms the fastest driver in a dataset of over 50,000 human players.


title: Fast algorithms for robust principal component analysis with an upper bound on the rank

abstract: The robust principal component analysis (RPCA) decomposes a data matrix into a low-rank part and a sparse part. There are mainly two types of algorithms for RPCA. The first type of algorithm applies regularization terms on the singular values of a matrix to obtain a low-rank matrix. However, calculating singular values can be very expensive for large matrices. The second type of algorithm replaces the low-rank matrix as the multiplication of two small matrices. They are faster than the first type because no singular value decomposition (SVD) is required. However, the rank of the low-rank matrix is required, and an accurate rank estimation is needed to obtain a reasonable solution. In this paper, we propose algorithms that combine both types. Our proposed algorithms require an upper bound of the rank and SVD on small matrices. First, they are faster than the first type because the cost of SVD on small matrices is negligible. Second, they are more robust than the second type because an upper bound of the rank instead of the exact rank is required. Furthermore, we apply the Gauss-Newton method to increase the speed of our algorithms. Numerical experiments show the better performance of our proposed algorithms.


title: Tailoring the Chirality of Microlaser with Topological Structured Chiral Droplets

abstract: The generation of chiral laser emission offers promising opportunities for modern photonic applications and the study of chiral light-mater interactions. Despite the great process made in recent years, the direct generation of chiral lasers with controllable chirality remains challenging in a microcavity. This study reports a strategy to control the emission chirality of whispering-gallery-mode organic microlasers with topological-structured chiral liquid crystal droplets. The findings suggest that the topological transformations in a microdroplet can induce different optical chirality strength and lasing polarization characteristics. In particular, the role of optical rotatory power was also investigated under linear and circular polarized excitation, where a vast lasing intensity difference between left/right circular polarized excitation was revealed. Theoretical analysis and simulation were carried out to support the significant findings. This study provides a facile tool for manipulating chiral laser emissions in microcavity, giving inspiration for topological-controlled photonics, chiral optical devices, and chiral molecular detection.


title: Enhanced anomalous Nernst effect in disordered Dirac and Weyl materials

abstract: We analyze the thermoelectric response of Dirac and Weyl semimetals using the semiclassical approach, focusing on the extrinsic contributions due to skew-scattering and side jump. Our results apply to linear response Nernst effect in ferromagnetic Dirac materials such as Fe$_3$Sn$_2$, Weyl semimetal Co$_3$Sn$_2$S$_2$ and to second order response of monolayer graphene on hBN with trigonal warping. Our analysis indicates that the extrinsic contributions can be a significant component of anomalous Nernst response and used to explain an enhanced thermoelectric response.


title: Revisiting $^{129}$Xe electric dipole moment measurements applying a new global phase fitting approach

abstract: By measuring the nuclear magnetic spin precession frequencies of polarized $^{129}$Xe and $^{3}$He, a new upper limit on the $^{129}$Xe atomic electric dipole moment (EDM) $ d_\mathrm{A} (^{129}\mathrm{Xe})$ was reported in Phys. Rev. Lett. 123, 143003 (2019). Here, we propose a new evaluation method based on global phase fitting (GPF) for analyzing the continuous phase development of the $^{3}$He-$^{129}$Xe comagnetometer signal. The Cramer-Rao Lower Bound on the $^{129}$Xe EDM for the GPF method is theoretically derived and shows the potential benefit of our new approach. The robustness of the GPF method is verified with Monte-Carlo studies. By optimizing the analysis parameters and adding data that could not be analyzed with the former method, we obtain a result of $d_\mathrm{A} (^{129}\mathrm{Xe}) = 1.1 \pm 3.6~\mathrm{(stat)} \pm 2.0~\mathrm{(syst)} \times 10^{-28}~ e~\mathrm{cm}$ in an unblinded analysis. For the systematic uncertainty analyses, we adopted all methods from the aforementioned PRL publication except the comagnetometer phase drift, which can be omitted using the GPF method. The updated null result can be interpreted as a new upper limit of $| d_\mathrm{A} (^{129}\mathrm{Xe}) | < 8.3 \times 10^{-28}~e~\mathrm{cm}$ at the 95\% C.L.


title: Integration of singular subalgebroids by diffeological groupoids

abstract: We establish an integration theory for singular subalgebroids, by diffeological groupoids. To do so, we single out a class of diffeological groupoids satisfying specific properties, and we introduce a differentiation-integration procedure under which they correspond to singular subalgebroids. Our definition of integration distinguishes the holonomy groupoid from the graph, although both differentiate to the original singular subalgebroid: the holonomy groupoid satisfies a certain submersive property, while the graph does not.


title: Frobenius nilHecke algebras

abstract: To any Frobenius superalgebra $A$ we associate towers of Frobenius nilCoxeter algebras and Frobenius nilHecke algebras. These act naturally, via Frobenius divided difference operators, on Frobenius polynomial algebras. When $A$ is the ground ring, our algebras recover the classical nilCoxeter and nilHecke algebras. When $A$ is the two-dimensional Clifford algebra, they are Morita equivalent to the odd nilCoxeter and odd nilHecke algebras.


title: Stochastic Adaptive Line Search for Differentially Private Optimization

abstract: The performance of private gradient-based optimization algorithms is highly dependent on the choice of step size (or learning rate) which often requires non-trivial amount of tuning. In this paper, we introduce a stochastic variant of classic backtracking line search algorithm that satisfies R\'enyi differential privacy. Specifically, the proposed algorithm adaptively chooses the step size satsisfying the the Armijo condition (with high probability) using noisy gradients and function estimates. Furthermore, to improve the probability with which the chosen step size satisfies the condition, it adjusts per-iteration privacy budget during runtime according to the reliability of noisy gradient. A naive implementation of the backtracking search algorithm may end up using unacceptably large privacy budget as the ability of adaptive step size selection comes at the cost of extra function evaluations. The proposed algorithm avoids this problem by using the sparse vector technique combined with the recent privacy amplification lemma. We also introduce a privacy budget adaptation strategy in which the algorithm adaptively increases the budget when it detects that directions pointed by consecutive gradients are drastically different. Extensive experiments on both convex and non-convex problems show that the adaptively chosen step sizes allow the proposed algorithm to efficiently use the privacy budget and show competitive performance against existing private optimizers.


title: Embedding a Heavy-Ball type of Momentum into the Estimating Sequences

abstract: We present a new accelerated gradient-based method for solving smooth unconstrained optimization problems. The goal is to embed a heavy-ball type of momentum into the Fast Gradient Method (FGM). For this purpose, we devise a generalization of the estimating sequences, which allows for encoding any form of information about the cost function that can aid in further accelerating the minimization process. In the black box framework, we propose a construction for the generalized estimating sequences, which is obtained by exploiting the history of the previously constructed estimating functions. From the viewpoint of efficiency estimates, we prove that the lower bound on the number of iterations for the proposed method is $\mathcal{O} \left(\sqrt{\frac{\kappa}{2}}\right)$. Our theoretical results are further corroborated by extensive numerical experiments on various types of optimization problems, often dealt within signal processing. Both synthetic and real-world datasets are utilized to demonstrate the efficiency of our proposed method in terms of decreasing the distance to the optimal solution, as well as in terms of decreasing the norm of the gradient.


title: Entanglement harvesting for Unruh-DeWitt detectors in circular motion

abstract: We study the properties of the transition probability and entanglement harvesting phenomenon for circularly accelerated detectors locally interacting with massless scalar fields. The dependence of the transition probability on the parameters associated with the circular motion is first analyzed in detail. By a cross-comparison with the situation of the uniformly accelerated motion, we obtain that the transition probability and the possible thermalization behavior for detectors rotating with an extremely large circular radius are analogous to that for uniformly accelerated detectors, but for a very small linear speed and a large acceleration, the effective temperature which characterizes the detectors' thermalization in a finite duration is much lower than that for uniformly accelerated detectors. We then focus on the phenomenon of entanglement harvesting in two special situations of circular trajectories, i.e., the coaxial rotation and the mutually perpendicular axial rotation by examining the concurrence as the entanglement measure in detail. We find that when two circularly accelerated detectors have equivalent acceleration and size of circular trajectory, the harvested entanglement rapidly decays with increasing acceleration or separation between two detectors. In contrast with the situation of uniform acceleration, the angular velocity would have significant impacts on entanglement harvesting. Especially for those detectors circularly moving in different directions, both the acceleration and trajectory radius play an important inhibiting role in entanglement harvesting. When two circularly accelerated detectors have different values of acceleration or angular velocity, we find that the entanglement can still be extracted by such detectors, even in the situation that one detector is at rest and the other is in a circular motion.


title: Comparison of Convolutional neural network training parameters for detecting Alzheimers disease and effect on visualization

abstract: Convolutional neural networks (CNN) have become a powerful tool for detecting patterns in image data. Recent papers report promising results in the domain of disease detection using brain MRI data. Despite the high accuracy obtained from CNN models for MRI data so far, almost no papers provided information on the features or image regions driving this accuracy as adequate methods were missing or challenging to apply. Recently, the toolbox iNNvestigate has become available, implementing various state of the art methods for deep learning visualizations. Currently, there is a great demand for a comparison of visualization algorithms to provide an overview of the practical usefulness and capability of these algorithms. Therefore, this thesis has two goals: 1. To systematically evaluate the influence of CNN hyper-parameters on model accuracy. 2. To compare various visualization methods with respect to the quality (i.e. randomness/focus, soundness).


title: A linearised inverse conductivity problem for the Maxwell system at a high frequency

abstract: We consider a linearised inverse conductivity problem for electromagnetic waves in a three dimensional bounded domain at a high time-harmonic frequency. Increasing stability bounds for the conductivity coefficient in the full Maxwell system and in a simplified transverse electric mode are derived. These bounds contain a Lipschitz term with a factor growing polynomially in terms of the frequency, a Holder term, and a logarithmic term which decays with respect to the frequency as a power. To validate this increasing stability numerically, we propose a reconstruction algorithm aiming at the recovery of sufficiently many Fourier modes of the conductivity. A numerical evidence sheds light on the influence of the growing frequency and confirms the improved resolution at higher frequencies.


title: Reinforcement Learning Evaluation and Solution for the Feedback Capacity of the Ising Channel with Large Alphabet

abstract: We propose a new method to compute the feedback capacity of unifilar finite state channels (FSCs) with memory using reinforcement learning (RL). The feedback capacity was previously estimated using its formulation as a Markov decision process (MDP) with dynamic programming (DP) algorithms. However, their computational complexity grows exponentially with the channel alphabet size. Therefore, we use RL, and specifically its ability to parameterize value functions and policies with neural networks, to evaluate numerically the feedback capacity of channels with a large alphabet size. The outcome of the RL algorithm is a numerical lower bound on the feedback capacity, which is used to reveal the structure of the optimal solution. The structure is modeled by a graph-based auxiliary random variable that is utilized to derive an analytic upper bound on the feedback capacity with the duality bound. The capacity computation is concluded by verifying the tightness of the upper bound by testing whether it is BCJR invariant. We demonstrate this method on the Ising channel with an arbitrary alphabet size. For an alphabet size smaller than or equal to 8, we derive the analytic solution of the capacity. Next, the structure of the numerical solution is used to deduce a simple coding scheme that achieves the feedback capacity and serves as a lower bound for larger alphabets. For an alphabet size greater than 8, we present an upper bound on the feedback capacity. For an asymptotically large alphabet size, we present an asymptotic optimal coding scheme.


title: A plane wave study on the localized-extended transitions in the one-dimensional incommensurate systems

abstract: Based on our recently proposed plane wave framework, we theoretically study the localized-extended transition in the one dimensional incommensurate systems with cosine type of potentials, which are in close connection to many recent experiments in the ultracold atom and photonic crystal. We formulate a propagator based scattering picture for the transition at the ground state and single particle mobility edge, in which the deeper connection between the incommensurate potentials, eigenstate compositions and transition mechanism is revealed. We further show that there exists a upper limit of localization length for all localized eigenstates, leading to an fundamental difference to the Anderson localization. Numerical calculations are presented alongside the analysis to justify our statements. The theoretical analysis and numerical methods can also be generalized to systems in higher dimensions, with different potentials or beyond the single particle regime, which would benefit the future studies in the related fields.


title: A port-Hamiltonian approach to modeling the structural dynamics of complex systems

abstract: With this contribution, we give a complete and comprehensive framework for modeling the dynamics of complex mechanical structures as port-Hamiltonian systems. This is motivated by research on the potential of lightweight construction using active load-bearing elements integrated into the structure. Such adaptive structures are of high complexity and very heterogeneous in nature. Port-Hamiltonian systems theory provides a promising approach for their modeling and control. Subsystem dynamics can be formulated in a domain-independent way and interconnected by means of power flows. The modular approach is also suitable for robust decentralized control schemes. Starting from a distributed-parameter port-Hamiltonian formulation of beam dynamics, we show the application of an existing structure-preserving mixed finite element method to arrive at finite-dimensional approximations. In contrast to the modeling of single bodies with a single boundary, we consider complex structures composed of many simple elements interconnected at the boundary. This is analogous to the usual way of modeling civil engineering structures which has not been transferred to port-Hamiltonian systems before. A block diagram representation of the interconnected systems is used to generate coupling constraints which leads to differential algebraic equations of index one. After the elimination of algebraic constraints, systems in input-state-output(ISO) port-Hamiltonian form are obtained. Port-Hamiltonian system models for the considered class of systems can also be constructed from the mass and stiffness matrices obtained via conventional finite element methods. We show how this relates to the presented approach and discuss the differences, promoting a better understanding across engineering disciplines.


title: Password Guessers Under a Microscope: An In-Depth Analysis to Inform Deployments

abstract: Password guessers are instrumental for assessing the strength of passwords. Despite their diversity and abundance, little is known about how different guessers compare to each other. We perform in-depth analyses and comparisons of the guessing abilities and behavior of password guessers. To extend analyses beyond number of passwords cracked, we devise an analytical framework to compare the types of passwords that guessers generate under various conditions (e.g., limited training data, limited number of guesses, and dissimilar training and target data). Our results show that guessers often produce dissimilar guesses, even when trained on the same data. We leverage this result to show that combinations of computationally-cheap guessers are as effective as computationally intensive guessers, but more efficient. Our insights allow us to provide a concrete set of recommendations for system administrators when performing password checking.


title: A Finite Graph Approach to the Probabilistic Hadwiger-Nelson Problem

abstract: We advance a probabilistic approach to the Hadwiger-Nelson problem initially developed by the Polymath16 project, in particular relating the approach to finite unit-distance graphs. We define the numerical \textit{badness} of a given $k$-coloring of the plane to be the probability that a randomly chosen unit-distance edge is monochromatic under the coloring, and we provide lower bounds on the badness of arbitrary $k$-colorings using a probabilistic technique relating to finite graphs. The contrapositive of the resulting bounds lets us compute lower bounds on the order of non $k$-colorable unit-distance graphs, improving bounds produced by Pritikin and the Polymath16 project in the $k = 4$ and $k = 5$ cases. Additionally, we make partial progress on a probabilistic analog of the de Bruijn-Erd\H{o}s compactness theorem.


title: Overdetermined boundary problems with nonconstant Dirichlet and Neumann data

abstract: In this paper we consider the overdetermined boundary problem for a general second order semilinear elliptic equation on bounded domains of $\mathbf{R}^n$, where one prescribes both the Dirichlet and Neumann data of the solution. We are interested in the case where the data are not necessarily constant and where the coefficients of the equation can depend on the position, so that the overdetermined problem does not generally admit a radial solution. Our main result is that, nevertheless, under minor technical hypotheses nontrivial solutions to the overdetermined boundary problem always exist.


title: Anomaly Detection with Convolutional Autoencoders for Fingerprint Presentation Attack Detection

abstract: In recent years, the popularity of fingerprint-based biometric authentication systems significantly increased. However, together with many advantages, biometric systems are still vulnerable to presentation attacks (PAs). In particular, this applies for unsupervised applications, where new attacks unknown to the system operator may occur. Therefore, presentation attack detection (PAD) methods are used to determine whether samples stem from a bona fide subject or from a presentation attack instrument (PAI). In this context, most works are dedicated to solve PAD as a two-class classification problem, which includes training a model on both bona fide and PA samples. In spite of the good detection rates reported, these methods still face difficulties detecting PAIs from unknown materials. To address this issue, we propose a new PAD technique based on autoencoders (AEs) trained only on bona fide samples (i.e. one-class), which are captured in the short wave infrared domain. On the experimental evaluation over a database of 19,711 bona fide and 4,339 PA images including 45 different PAI species, a detection equal error rate (D-EER) of 2.00% was achieved. Additionally, our best performing AE model is compared to further one-class classifiers (support vector machine, Gaussian mixture model). The results show the effectiveness of the AE model as it significantly outperforms the previously proposed methods.


title: Inducing a many-body topological state of matter through Coulomb-engineered local interactions

abstract: The engineering of artificial systems hosting topological excitations is at the heart of current condensed matter research. Most of these efforts focus on single-particle properties neglecting possible engineering routes via the modifications of the fundamental many-body interactions. Interestingly, recent experimental breakthroughs have shown that Coulomb interactions can be efficiently controlled by substrate screening engineering. Inspired by this success } we propose a simple platform in which topologically non-trivial many-body excitations emerge solely from dielectrically-engineered Coulomb interactions in an otherwise topologically trivial single-particle band structure. Furthermore, by performing a realistic microscopic modeling of screening engineering, we demonstrate how our proposal can be realized in one-dimensional systems such as quantum-dot chains. Our results put forward Coulomb engineering as a powerful tool to create topological excitations, with potential applications in a variety of solid-state platforms.


title: Generators of the plane Cremona group over the field with two elements

abstract: The plane Cremona group over the finite field $\mathbb{F}_2$ is generated by three infinite families and finitely many birational maps with small base orbits. One family preserves the pencil of lines through a point, the other two preserve the pencil of conics through four points that form either one Galois orbit of size 4, or two Galois orbits of size 2. For each family, we give a generating set that is parametrized by the rational functions over $\mathbb{F}_2$. Moreover, we describe the finitely many remaining maps and give an upper bound on the number needed to generate the Cremona group. Finally, we prove that the plane Cremona group over $\mathbb{F}_2$ is generated by involutions.


title: High performance Wannier interpolation of Berry curvature and related quantities: WannierBerri code

abstract: Wannier interpolation is a powerful tool for performing Brillouin zone integrals over dense grids of $\mathbf{k}$ points, which are essential to evaluate such quantities as the intrinsic anomalous Hall conductivity or Boltzmann transport coefficients. However, new physical problems and new materials create new numerical challenges, and computations with the existing codes become very expensive, which often prevents reaching the desired accuracy. In this article, I present a series of methods that boost the speed of Wannier interpolation by several orders of magnitude. They include a combination of fast and slow Fourier transforms, explicit use of symmetries and recursive adaptive grid refinement among others. The proposed methodology has been implemented in the new python code WannierBerri, which also aims to serve as a convenient platform for the future development of interpolation schemes for novel phenomena.


title: XNAP: Making LSTM-based Next Activity Predictions Explainable by Using LRP

abstract: Predictive business process monitoring (PBPM) is a class of techniques designed to predict behaviour, such as next activities, in running traces. PBPM techniques aim to improve process performance by providing predictions to process analysts, supporting them in their decision making. However, the PBPM techniques` limited predictive quality was considered as the essential obstacle for establishing such techniques in practice. With the use of deep neural networks (DNNs), the techniques` predictive quality could be improved for tasks like the next activity prediction. While DNNs achieve a promising predictive quality, they still lack comprehensibility due to their hierarchical approach of learning representations. Nevertheless, process analysts need to comprehend the cause of a prediction to identify intervention mechanisms that might affect the decision making to secure process performance. In this paper, we propose XNAP, the first explainable, DNN-based PBPM technique for the next activity prediction. XNAP integrates a layer-wise relevance propagation method from the field of explainable artificial intelligence to make predictions of a long short-term memory DNN explainable by providing relevance values for activities. We show the benefit of our approach through two real-life event logs.


title: Pure de Sitter space and the island moving back in time

abstract: Observers in de Sitter space can only access the space up to their cosmological horizon. Assuming thermal equilibrium, we use the quantum Ryu-Takayanagi or island formula to compute the entanglement entropy between the states inside the cosmological horizon and states outside, as a function of time. We obtain a Page curve that is bound at a value corresponding to the Gibbons-Hawking entropy. At this transition an 'island' forms, which is in a significantly different location as compared to when considering black hole horizons and even moves back in time. These differences turn out to be essential for non-violation of the no-cloning theorem in combination with entanglement wedge reconstruction. This consideration furthermore introduces the need for a scrambling time, the entropy dependence of which turns out to coincide with what is expected for black holes. The model we employ has pure three-dimensional de Sitter space as a solution. We dimensionally reduce to two dimensions in order to take into account semi-classical effects. Nevertheless, we expect the aforementioned qualitative features of the island to persist in higher dimensions.


title: How Archimedes showed that $\pi$ is approximately equal to 22/7

abstract: The ratio of the circumference, C, of a circle to its diameter, D, is a constant number denoted by $\pi$ and is independent of the size of the circle. It is known that $\pi$ is an irrational number and therefore cannot be expressed as a common fraction. Its value is approximately equal to 3.141592. Since Archimedes was one of the first persons to suggest a rational approximation of 22/7 for $\pi$, it is sometimes referred to as Archimedes' constant. In this article, we discuss how Archimedes came up with his formula. Archimedes in fact proved that 223/71 < $\pi$ < 22/7. Here we provide an improved lower bound.


title: Mining Large Quasi-cliques with Quality Guarantees from Vertex Neighborhoods

abstract: Mining dense subgraphs is an important primitive across a spectrum of graph-mining tasks. In this work, we formally establish that two recurring characteristics of real-world graphs, namely heavy-tailed degree distributions and large clustering coefficients, imply the existence of substantially large vertex neighborhoods with high edge-density. This observation suggests a very simple approach for extracting large quasi-cliques: simply scan the vertex neighborhoods, compute the clustering coefficient of each vertex, and output the best such subgraph. The implementation of such a method requires counting the triangles in a graph, which is a well-studied problem in graph mining. When empirically tested across a number of real-world graphs, this approach reveals a surprise: vertex neighborhoods include maximal cliques of non-trivial sizes, and the density of the best neighborhood often compares favorably to subgraphs produced by dedicated algorithms for maximizing subgraph density. For graphs with small clustering coefficients, we demonstrate that small vertex neighborhoods can be refined using a local-search method to ``grow'' larger cliques and near-cliques. Our results indicate that contrary to worst-case theoretical results, mining cliques and quasi-cliques of non-trivial sizes from real-world graphs is often not a difficult problem, and provides motivation for further work geared towards a better explanation of these empirical successes.


title: Deepcode and Modulo-SK are Designed for Different Settings

abstract: We respond to [1] which claimed that "Modulo-SK scheme outperforms Deepcode [2]". We demonstrate that this statement is not true: the two schemes are designed and evaluated for entirely different settings. DeepCode is designed and evaluated for the AWGN channel with (potentially delayed) uncoded output feedback. Modulo-SK is evaluated on the AWGN channel with coded feedback and unit delay. [1] also claimed an implementation of Schalkwijk and Kailath (SK) [3] which was numerically stable for any number of information bits and iterations. However, we observe that while their implementation does marginally improve over ours, it also suffers from a fundamental issue with precision. Finally, we show that Deepcode dominates the optimized performance of SK, over a natural choice of parameterizations when the feedback is noisy.


title: Search for Nearby Earth Analogs. III. Detection of ten new planets, three planet candidates, and confirmation of three planets around eleven nearby M dwarfs

abstract: Earth-sized planets in the habitable zones of M dwarfs are good candidates for the study of habitability and detection of biosignatures. To search for these planets, we analyze all available radial velocity data and apply four signal detection criteria to select the optimal candidates. We find ten strong candidates satisfying these criteria and three weak candidates showing inconsistency over time due to data samplings. We also confirm three previous planet candidates and improve their orbital solutions through combined analyses of updated data sets. Among the strong planet candidates, HIP 38594 b is a temperate super-Earth with a mass of $8.2 \pm 1.7$ $M_\oplus$ and an orbital period of $60.7\pm0.1$ days, orbiting around an early-type M dwarf. Early-type M dwarfs are less active and thus are better hosts for habitable planets than mid-type and late-type M dwarfs. Moreover, we report the detection of five two-planet systems, including two systems made up of a warm or cold Neptune and a cold Jupiter, consistent with a positive correlation between these two types of planets. We also detect three temperate Neptunes, four cold Neptunes, and four cold Jupiters, contributing to a rarely explored planet population. Due to their proximity to the Sun, these planets on wide orbits are appropriate targets for direct imaging by future facilities such as HabEx and ELT.


title: Classification of generic spherical quadrilaterals

abstract: Generic spherical quadrilaterals are classified up to isometry. Condition of genericity consists in the requirement that the images of the sides under the developing map belong to four distinct circles which have no triple intersections. Under this condition, it is shown that the space of quadrilaterals with prescribed angles consists of finitely many open curves. Degeneration at the endpoints of these curves is also determined.


title: The time-dependent angular analysis of $B_d\to K_S\ell\ell$, a new benchmark for new physics

abstract: We consider the time-dependent analysis of $B_d\to K_S\ell\ell$ taking into account the time-evolution of the $B_d$ meson and its mixing into $\bar{B}_d$. We discuss the angular conventions required to define the angular observables in a transparent way with respect to CP conjugation. The inclusion of time evolution allows us to identify six new observables, out of which three could be accessed from a time-dependent tagged analysis. We also show that these observables could be obtained by time-integrated measurements in a hadronic environment if flavour tagging is available. We provide simple and precise predictions for these observables in the SM and in NP models with real contributions to SM and chirally flipped operators, which are independent of form factors and charm-loop contributions. As such, these observables provide robust and powerful cross-checks of the New Physics scenarios currently favoured by global fits to $b\to s\ell\ell$ data. In addition, we discuss the sensitivity of these observables with respect to NP scenarios involving scalar and tensor operators, or CP-violating phases. We illustrate how these new observables can provide a benchmark to discriminate among the various NP scenarios in $b\to s\mu\mu$. We discuss the extension of these results for $B_s$ decays into $f_0$, $\eta$ or $\eta'$.


title: Offloading Optimization in Edge Computing for Deep Learning Enabled Target Tracking by Internet-of-UAVs

abstract: The empowering unmanned aerial vehicles (UAVs) have been extensively used in providing intelligence such as target tracking. In our field experiments, a pre-trained convolutional neural network (CNN) is deployed at the UAV to identify a target (a vehicle) from the captured video frames and enable the UAV to keep tracking. However, this kind of visual target tracking demands a lot of computational resources due to the desired high inference accuracy and stringent delay requirement. This motivates us to consider offloading this type of deep learning (DL) tasks to a mobile edge computing (MEC) server due to limited computational resource and energy budget of the UAV, and further improve the inference accuracy. Specifically, we propose a novel hierarchical DL tasks distribution framework, where the UAV is embedded with lower layers of the pre-trained CNN model, while the MEC server with rich computing resources will handle the higher layers of the CNN model. An optimization problem is formulated to minimize the weighted-sum cost including the tracking delay and energy consumption introduced by communication and computing of the UAVs, while taking into account the quality of data (e.g., video frames) input to the DL model and the inference errors. Analytical results are obtained and insights are provided to understand the tradeoff between the weighted-sum cost and inference error rate in the proposed framework. Numerical results demonstrate the effectiveness of the proposed offloading framework.


title: Two-pole structures in a relativistic Friedrichs-Lee-QPC scheme

abstract: A general appearance of two-pole structures is exhibited in a relativistic Friedrichs-Lee model combined with a relativistic quark pair creation model in a consistent manner. This kind of two-pole structure could be found when a $q\bar q$ state couples to the open-flavor continuum state in the $S$ partial wave. We found that many enigmatic states, such as $f_0(500)/\sigma$, $K_0^*(700)/\kappa$, $a_0(980)$, $f_0(980)$, $D_0^*(2300)$, $D_{s0}^*(2317)$, and $X(3872)$, together with another higher state for each, all result from this kind of two-pole structures. Furthermore, an interesting observation is that this kind of two-pole structure will contribute roughly a total of 180$^\circ$ phase shift for the scattering process in a single channel approximation. This relativistic scheme may provide more insights into the understanding of the properties of non-$q\bar q$ state. It is also suggested that such two-pole structure could be a common phenomenon which deserves studying both from theoretical and experimental perspectives.


title: Three-dimensional seismic characterization and imaging of the Soda Lake geothermal field

abstract: Accurate characterization of subsurface geophysical properties and detection of the fault system are essential for geothermal energy exploration and production. The Soda Lake geothermal field is in western Nevada with a complex fault system. Previous seismic characterization only produced a low-resolution, smooth velocity model along with a simple, conceptual fault model. Using optimized correlation-based full-waveform inversion, wavefield-separation-based reverse-time migration, and automatic fault detection techniques, we present 3D seismic characterization for the Soda Lake geothermal field using 3D surface seismic data acquired with Vibroseis sources. We obtain 3D high-resolution velocity, density, and acoustic impedance models, 3D seismic images with different grid spacings, and a high-resolution fault system. Consistency check between the constructed faults and currently active injection and production geothermal wells verifies that our seismic inversion and imaging results and detected faults are reliable. These results can provide valuable information for optimizing well placement and geothermal energy production at the Soda Lake geothermal field.


title: Forecasting day-ahead electricity prices: A review of state-of-the-art algorithms, best practices and an open-access benchmark

abstract: While the field of electricity price forecasting has benefited from plenty of contributions in the last two decades, it arguably lacks a rigorous approach to evaluating new predictive algorithms. The latter are often compared using unique, not publicly available datasets and across too short and limited to one market test samples. The proposed new methods are rarely benchmarked against well established and well performing simpler models, the accuracy metrics are sometimes inadequate and testing the significance of differences in predictive performance is seldom conducted. Consequently, it is not clear which methods perform well nor what are the best practices when forecasting electricity prices. In this paper, we tackle these issues by performing a literature survey of state-of-the-art models, comparing state-of-the-art statistical and deep learning methods across multiple years and markets, and by putting forward a set of best practices. In addition, we make available the considered datasets, forecasts of the state-of-the-art models, and a specifically designed python toolbox, so that new algorithms can be rigorously evaluated in future studies.


title: Reinforcement Learning for Improving Object Detection

abstract: The performance of a trained object detection neural network depends a lot on the image quality. Generally, images are pre-processed before feeding them into the neural network and domain knowledge about the image dataset is used to choose the pre-processing techniques. In this paper, we introduce an algorithm called ObjectRL to choose the amount of a particular pre-processing to be applied to improve the object detection performances of pre-trained networks. The main motivation for ObjectRL is that an image which looks good to a human eye may not necessarily be the optimal one for a pre-trained object detector to detect objects.


title: Neural networks in day-ahead electricity price forecasting: Single vs. multiple outputs

abstract: Recent advancements in the fields of artificial intelligence and machine learning methods resulted in a significant increase of their popularity in the literature, including electricity price forecasting. Said methods cover a very broad spectrum, from decision trees, through random forests to various artificial neural network models and hybrid approaches. In electricity price forecasting, neural networks are the most popular machine learning method as they provide a non-linear counterpart for well-tested linear regression models. Their application, however, is not straightforward, with multiple implementation factors to consider. One of such factors is the network's structure. This paper provides a comprehensive comparison of two most common structures when using the deep neural networks -- one that focuses on each hour of the day separately, and one that reflects the daily auction structure and models vectors of the prices. The results show a significant accuracy advantage of using the latter, confirmed on data from five distinct power exchanges.


title: Differentially Private Clustering: Tight Approximation Ratios

abstract: We study the task of differentially private clustering. For several basic clustering problems, including Euclidean DensestBall, 1-Cluster, k-means, and k-median, we give efficient differentially private algorithms that achieve essentially the same approximation ratios as those that can be obtained by any non-private algorithm, while incurring only small additive errors. This improves upon existing efficient algorithms that only achieve some large constant approximation factors. Our results also imply an improved algorithm for the Sample and Aggregate privacy framework. Furthermore, we show that one of the tools used in our 1-Cluster algorithm can be employed to get a faster quantum algorithm for ClosestPair in a moderate number of dimensions.


title: On the cavity size in circumbinary discs

abstract: How does the cavity size in circumbinary discs depend on disc and binary properties? We investigate by simulating disc cavities carved by binary companions using smoothed particle hydrodynamics (SPH). We find that a cavity is quickly opened on the dynamical time, while the cavity size is set on the viscous time. In agreement with previous findings, we find long term cavity sizes of 2-5 times the binary semi-major axis, increasing with eccentricity and decreasing with disc aspect ratio. When considering binaries inclined with respect to the disc we find three regimes: i) discs that evolve towards a coplanar orbit have a large cavity, slightly smaller than that of an initially coplanar disc; ii) discs that evolve towards a polar orbit by breaking have a small cavity, equal in size to that of an initially polar disc; iii) discs that evolve towards a polar orbit via warping have an intermediate-sized cavity. We find typical gas depletions inside the cavity of $\gtrsim 2$ orders of magnitude in surface density.


title: On the base point locus of surface parametrizations: formulas and consequences

abstract: This paper shows that the multiplicity of the base points locus of a projective rational surface parametrization can be expressed as the degree of the content of a univariate resultant. As a consequence, we get a new proof of the degree formula relating the degree of the surface, the degree of the parametrization, the base points multiplicity, and the degree of the rational map induced by the parametrization. In addition, we extend both formulas to the case of dominant rational maps of the projective plane and describe how the base point loci of a parametrization and its reparametrizations are related. As an application of these results, we explore how the degree of a surface reparametrization is affected by the presence of base points.


title: The role of rare events in the pinning problem

abstract: Type II superconductors exhibit a fascinating phenomenology that is determined by the dynamical properties of the vortex matter hosted by the material. A crucial element in this phenomenology is vortex pinning by material defects, e.g., immobilizing vortices at small drives and thereby guaranteeing dissipation-free current flow. Pinning models for vortices and other topological defects, such as domain walls in magnets or dislocations in crystals, come in two standard variants: i) weak collective pinning, where individual weak defects are unable to pin, while the random accumulation of many force centers within a collective pinning volume combines into an effective pin, and ii) strong pinning, where strong defects produce large vortex displacements and bistabilities that lead to pinning on the level of individual defects. The transition between strong and weak pinning is quantified by the Labusch criterion $\kappa \approx f_p/\bar{C}\xi = 1$, where $f_p$ and $\bar{C}$ are the force of one defect and the effective elasticity of the vortex lattice, respectively ($\xi$ is the coherence length). Here, we show that a third generic type of pinning becomes dominant when the pinning force $f_p$ enters the weak regime, the pinning by rare events. We find that within an intermediate regime $1/2 < \kappa < 1$, compact pairs of weak defects define strong pinning clusters that extend the mechanism of strong pinning into the weak regime. We present a detailed analysis of this cluster-pinning mechanism and show that its pinning-force density parametrically dominates over the weak pinning result. The present work is a first attempt to include correlations between defects into the discussion of strong pinning.


title: A Computer-Assisted Study of Red Coral Population Dynamics

abstract: We consider a 13-dimensional age-structured discrete red coral population model varying with respect to a fitness parameter. Our numerical results give a bifurcation diagram of both equilibria and stable invariant curves of orbits. We observe that not only for low levels of fitness, but also for high levels of fitness, populations are extremely vulnerable, in that they spend long time periods near extinction. We then use computer-assisted proofs techniques to rigorously validate the set of regular and bifurcation fixed points that have been found numerically.


title: Linguistically-aware Attention for Reducing the Semantic-Gap in Vision-Language Tasks

abstract: Attention models are widely used in Vision-language (V-L) tasks to perform the visual-textual correlation. Humans perform such a correlation with a strong linguistic understanding of the visual world. However, even the best performing attention model in V-L tasks lacks such a high-level linguistic understanding, thus creating a semantic gap between the modalities. In this paper, we propose an attention mechanism - Linguistically-aware Attention (LAT) - that leverages object attributes obtained from generic object detectors along with pre-trained language models to reduce this semantic gap. LAT represents visual and textual modalities in a common linguistically-rich space, thus providing linguistic awareness to the attention process. We apply and demonstrate the effectiveness of LAT in three V-L tasks: Counting-VQA, VQA, and Image captioning. In Counting-VQA, we propose a novel counting-specific VQA model to predict an intuitive count and achieve state-of-the-art results on five datasets. In VQA and Captioning, we show the generic nature and effectiveness of LAT by adapting it into various baselines and consistently improving their performance.


title: Stability of a point charge for the Vlasov-Poisson system: the radial case

abstract: We consider the Vlasov-Poisson system with initial data a small, radial, absolutely continuous perturbation of a point charge. We show that the solution is global and disperses to infinity via a modified scattering along trajectories of the linearized flow. This is done by an exact integration of the linearized equation, followed by the analysis of the perturbed Hamiltonian equation in action-angle coordinates.


title: Deconstructing the Decentralization Trilemma

abstract: The vast majority of applications at this moment rely on centralized servers to relay messages between clients, where these servers are considered trusted third-parties. With the rise of blockchain technologies over the last few years, there has been a move away from both centralized servers and traditional federated models to more decentralized peer-to-peer alternatives. However, there appears to be a trilemma between security, scalability, and decentralization in blockchain-based systems. Deconstructing this trilemma using well-known threat models, we define a typology of centralized, federated, and decentralized architectures. Each of the different architectures has this trilemma play out differently. Facing a possible decentralized future, we outline seven hard problems facing decentralization and theorize that the differences between centralized, federated, and decentralized architectures depend on differing social interpretations of trust.


title: Enhancing the Erd\H{o}s-Lov\'asz Tihany Conjecture for line graphs of multigraphs

abstract: In this paper, we prove an enhanced version of the Erd\H{o}s-Lov\'asz Tihany Conjecture for line graphs of multigraphs. That is, for every graph $G$ whose chromatic number $\chi(G)$ is more than its clique number $\omega(G)$ and for nonnegative integer $\ell$, any two integers $s,t \geq 3.5\ell+2$ with $s+t = \chi(G)+1$, there is a partition $(S,T)$ of the vertex set $V(G)$ such that $\chi(G[S])\geq s$ and $\chi(G[T])\geq t+\ell$. In particular, when $\ell=1$, we can obtain the same result just for any $s,t\geq4$. The Erd\H{o}s-Lov\'asz Tihany conjecture is a special case when $\ell=0$.


title: MaskedFace-Net -- A Dataset of Correctly/Incorrectly Masked Face Images in the Context of COVID-19

abstract: The wearing of the face masks appears as a solution for limiting the spread of COVID-19. In this context, efficient recognition systems are expected for checking that people faces are masked in regulated areas. To perform this task, a large dataset of masked faces is necessary for training deep learning models towards detecting people wearing masks and those not wearing masks. Some large datasets of masked faces are available in the literature. However, at the moment, there are no available large dataset of masked face images that permits to check if detected masked faces are correctly worn or not. Indeed, many people are not correctly wearing their masks due to bad practices, bad behaviors or vulnerability of individuals (e.g., children, old people). For these reasons, several mask wearing campaigns intend to sensitize people about this problem and good practices. In this sense, this work proposes three types of masked face detection dataset; namely, the Correctly Masked Face Dataset (CMFD), the Incorrectly Masked Face Dataset (IMFD) and their combination for the global masked face detection (MaskedFace-Net). Realistic masked face datasets are proposed with a twofold objective: i) to detect people having their faces masked or not masked, ii) to detect faces having their masks correctly worn or incorrectly worn (e.g.; at airport portals or in crowds). To the best of our knowledge, no large dataset of masked faces provides such a granularity of classification towards permitting mask wearing analysis. Moreover, this work globally presents the applied mask-to-face deformable model for permitting the generation of other masked face images, notably with specific masks. Our datasets of masked face images (137,016 images) are available at https://github.com/cabani/MaskedFace-Net.


title: Enhancing the Erd\H{o}s-Lov\'asz Tihany Conjecture for graphs with independence number two

abstract: Let $s\ge2$ and $t\ge2$ be integers. A graph $G$ is $(s,t)$-\emph{splittable} if $V(G)$ can be partitioned into two sets $S$ and $T$ such that $\chi(G[S])\geq s$ and $\chi(G[T])\geq t$. The well-known Erd\H{o}s-Lov\'asz Tihany Conjecture from 1968 states that every graph $G$ whose chromatic number $\chi(G)=s+t-1$ is more than its clique number $\omega(G)$ is $(s,t)$-splittable. In this paper, we prove an enhanced version of the Erd\H{o}s-Lov\'asz Tihany Conjecture for graphs with independence number two. That is, for every graph $G$ with $\chi(G)=s+t-1>\omega(G)+1$ is $(s,t+1)$-splittable. There are examples showing that this result is best possible.


title: Multichannel quantum-defect theory for anisotropic interactions

abstract: We present a general formulation of multichannel quantum-defect theory (MQDT) for anisotropic long-range potentials. The theory unifies the treatment of atomic and molecular interactions of all types, and greatly expands the set of interactions that can be treated and understood systematically, including complex interactions involving molecules. In one exemplary manifestation, the theory provides a methodology to make the classification of atomic interactions based on the Periodic Table quantitative, instead of qualitative, and to generalize the Table to include molecular classes. Through the concept of effective potential, the theory further establishes a foundation for new classes of quantum theories for chemistry and for a broad range of quantum systems made of either a few or many atoms and/or molecules.


title: Multi-critical topological transition at quantum criticality

abstract: The investigation and characterization of topological quantum phase transition between gapless phases is one of the recent interest of research in topological states of matter. We consider transverse field Ising model with three spin interaction in one dimension and observe a topological transition between gapless phases on one of the critical lines of this model. We study the distinct nature of these gapless phases and show that they belong to different universality classes. The topological invariant number (winding number) characterize different topological phases for the different regime of parameter space. We observe the evidence of two multi-critical points, one is topologically trivial and the other one is topologically active. Topological quantum phase transition between the gapless phases on the critical line occurs through the non-trivial multi-critical point in the Lifshitz universality class. We calculate and analyze the behavior of Wannier state correlation function close to the multi-critical point and confirm the topological transition between gapless phases. We show the breakdown of Lorentz invariance at this multi-critical point through the energy dispersion analysis. We also show that the scaling theories and curvature function renormalization group can also be effectively used to understand the topological quantum phase transitions between gapless phases. The model Hamiltonian which we study is more applicable for the system with gapless excitations, where the conventional concept of topological quantum phase transition fails.


title: V Tree -- Continued Fraction Expansion, Stern-Brocot Tree, Minkowski's $?(x)$ Function In Binary: Exponentially Faster

abstract: The Stern-Brocot tree and Minkowki's question mark function $?(x)$ (or Conway's box function) are related to the continued fraction expansion of numbers from Q with unary encoding of the partial denominators. We first define binary encodings $C_I, C_{II}$ of the natural numbers, adapted to the Gau\ss-Kuz'min measure for the distribution of partial denominators. We then define the V$_1$ tree as analogue to the Stern-Brocot tree, using the binary encondings $C_I, C_{II}$. We shall see that all numbers with denominator $q$ are present in the first $3.44\log_2(q)$ levels, instead of $1/q$ appearing in level $q$ in the Stern-Brocot tree. The extension of the V$_1$ tree, the V tree, covers all numbers from Q exactly once. We also define the binary version of Minkowski's question mark function, $?_V$, and conjecture that it has no derivative at rational points (for the original, $?'(x)=0, x\in Q$).


title: Sumsets of Semiconvex sets

abstract: We investigate additive properties of sets $A,$ where $A=\{a_1,a_2,\ldots ,a_k\}$ is a monotone increasing set of real numbers, and the differences of consecutive elements are all distinct. It is known that $|A+B|\geq c|A||B|^{1/2}$ for any finite set of numbers $B.$ The bound is tight up to the constant multiplier. We give a new proof to this result using bounds on crossing numbers of geometric graphs. We construct examples showing the limits of possible improvements. In particular, we show that there are arbitrarily large sets with different consecutive differences and sub-quadratic sumset sizes.


title: Quantum backflow in a ring

abstract: Free motion of a quantum particle with the wave function entirely comprised of plane waves with non-negative momenta may be accompanied by negative probability current, an effect called quantum backflow. The effect is weak and fragile, and has not yet been observed experimentally. Here we show that quantum backflow becomes significantly more pronounced and more amenable to experimental observation if, instead of letting the particle move along a straight line, one forces it to move in a circular ring.


title: Multilanguage Number Plate Detection using Convolutional Neural Networks

abstract: Object Detection is a popular field of research for recent technologies. In recent years, profound learning performance attracts the researchers to use it in many applications. Number plate (NP) detection and classification is analyzed over decades however, it needs approaches which are more precise and state, language and design independent since cars are now moving from state to another easily. In this paperwe suggest a new strategy to detect NP and comprehend the nation, language and layout of NPs. YOLOv2 sensor with ResNet attribute extractor heart is proposed for NP detection and a brand new convolutional neural network architecture is suggested to classify NPs. The detector achieves average precision of 99.57% and country, language and layout classification precision of 99.33%. The results outperforms the majority of the previous works and can move the area forward toward international NP detection and recognition.


title: Self-supervised Denoising via Diffeomorphic Template Estimation: Application to Optical Coherence Tomography

abstract: Optical Coherence Tomography (OCT) is pervasive in both the research and clinical practice of Ophthalmology. However, OCT images are strongly corrupted by noise, limiting their interpretation. Current OCT denoisers leverage assumptions on noise distributions or generate targets for training deep supervised denoisers via averaging of repeat acquisitions. However, recent self-supervised advances allow the training of deep denoising networks using only repeat acquisitions without clean targets as ground truth, reducing the burden of supervised learning. Despite the clear advantages of self-supervised methods, their use is precluded as OCT shows strong structural deformations even between sequential scans of the same subject due to involuntary eye motion. Further, direct nonlinear alignment of repeats induces correlation of the noise between images. In this paper, we propose a joint diffeomorphic template estimation and denoising framework which enables the use of self-supervised denoising for motion deformed repeat acquisitions, without empirically registering their noise realizations. Strong qualitative and quantitative improvements are achieved in denoising OCT images, with generic utility in any imaging modality amenable to multiple exposures.


title: How to organize a hackathon -- A planning kit

abstract: Hackathons and similar time-bounded events have become a global phenomenon. Their proliferation in various domains and their usefulness for a variety of goals has subsequently led to the emergence of different formats. While there are a multitude of guidelines available on how to prepare and run a hackathon, most of them focus on a particular format that was created for a specific purpose within a domain for a certain type of participants. This makes it difficult in particular for novice organizers to decide how to run an event that fits their needs. To address this gap we developed a planning kit that is organized around 12 key decision that organizers need to make when preparing and running a hackathon, and the tradeoffs that drive decision-making. The main planning kit is available online while this report is meant as a downloadable and citable resource.


title: Discovery potentials of double-charm tetraquarks

abstract: In this study, we investigate the discovery potential of double-charm tetraquarks $T^{\{cc\}}_{[\bar{q}\bar{q}']}$. We find that their production cross sections at the LHCb with $\sqrt{s} = 13$ TeV reach $\mathcal{O}(10^4)$ pb, which indicates that the LHCb has collected $\mathcal{O}(10^8)$ such particles. Through the decay channels of $T^{\{cc\}}_{[\bar{u}\bar{d}]}\to D^{+}K^{-}\pi^{+}$ or $D^0D^+\gamma$ (if stable) or $T^{\{cc\}}_{[\bar{u}\bar{d}]}\to D^0D^{*+}\to D^0D^0\pi^+$ (if unstable), it is highly expected that they get discovered at the LHCb in the near future. We also discuss the productions and decays of the double-charm tetraquarks at future Tera-$Z$ factories.


title: Uncovering the primary X-ray emission and possible starburst component in the polarized NLS1 Mrk 1239

abstract: X-ray observations of the unique NLS1 galaxy Mrk 1239 spanning 18-years are presented. Data from XMM-Newton, Suzaku, Swift and NuSTAR are combined to obtain a broad-band, multi-epoch view of the source. There is spectral variability in the 3-10 keV band over the 18-years. An analysis of the NuSTAR and Suzaku light curves also suggests rapid variability in the 3-10 keV band, which is consistent with the NLS1 definition of the source. However, no variability is seen below 3 keV on any timescale. Two distinct physical models are adopted to describe the data above and below ~3 keV. The low energies are dominated by a hot, diffuse gas likely associated with a starburst component at large physical scales. The higher energy spectrum is dominated by emission from the central region. Ionised partial covering and relativistic blurred reflection are considered for the central region emission. In both cases, the underlying power law has a photon index of $\Gamma\sim2.3-2.4$. A distant reflector, a neutral partial covering component with a covering fraction near $\sim1$, and contributions from starburst emission are always required. The blurred reflection model requires a reflection dominated spectrum, which may be at odds with the low emissivity index and radio properties of the source. By contrast, the two absorption components required in the ionised partial covering model may correspond to the two distinct regions of polarization observed in the optical. Regardless of the physical model, spectral changes between epochs are driven by the absorption components and on short time scales by intrinsic AGN variability.


title: A Harnack inequality for weak solutions of the Finsler $\gamma$-Laplacian

abstract: We study regularity of the Finsler $\gamma$-Laplacian, a general class of degenerate elliptic PDEs which naturally appear in anisotropic geometric problems. Precisely, given any strictly convex family of $C^{1}$-norms $\{ \rho_{x}\}$ on $\mathbb{R}^{n}$ and $\gamma > 1$, we consider the $W^{1,\gamma}(\Omega)$ solutions of the anisotropic PDE $$ \displaystyle \int_{\Omega} \left \langle \rho_{x}(Du)^{\gamma-1} (D \rho_{x})(Du), D \varphi \right \rangle = \int_{\Omega} \vec{F} \cdot D \varphi + f \varphi \qquad \forall \varphi \in W^{1,\gamma^{\prime}}_{0}(\Omega). $$ Under the mild assumption $|\xi|^{-1} \rho_{x}( \xi) \in [\nu, \Lambda]$ for all $(x,\xi) \in \Omega \times \mathbb{R}^{n}$ and some $0 < \nu \le \Lambda < \infty$ we perform a Moser iteration, verifying that sub- and super-solutions satisfy one-sided $\| \cdot \|_{\infty}$ bounds, which together imply solutions are locally bounded. When $u$ is non-negative this also implies a (weak) Harnack inequality. If $f, \vec{F} \equiv 0$ weak solutions also benefit from a strong maximum principle, and a Liouville-type theorem.


title: Overconvergent modular forms are highest weight vectors in the Hodge-Tate weight zero part of completed cohomology

abstract: We construct a $(\mathfrak{gl}_2, B(\mathbb{Q}_p))$ and Hecke-equivariant cup product pairing between overconvergent modular forms and the local cohomology at $0$ of a sheaf on $\mathbb{P}^1$, landing in the compactly supported completed $\mathbb{C}_p$-cohomology of the modular curve. The local cohomology group is a highest-weight Verma module, and the cup product is non-trivial on a highest weight vector for any overconvergent modular form of infinitesimal weight not equal to $1$. For classical weight $k\geq 2$, the Verma has an algebraic quotient $H^1(\mathbb{P}^1, \mathcal{O}(-k))$, and on classical forms the pairing factors through this quotient, giving a geometric description of "half" of the locally algebraic vectors in completed cohomology; the other half is described by a pairing with the roles of $H^1$ and $H^0$ reversed between the modular curve and $\mathbb{P}^1$. Under minor assumptions, we deduce a conjecture of Gouvea on the Hodge-Tate-Sen weights of Galois representations attached to overconvergent modular forms. Our main results are essentially a strict subset of those obtained independently by Lue Pan in arXiv:2008.07099, but the perspective here is different and the proofs are short and use simple tools: a Mayer-Vietoris cover, a cup product, and a boundary map in group cohomology.


title: Gradients as a Measure of Uncertainty in Neural Networks

abstract: Despite tremendous success of modern neural networks, they are known to be overconfident even when the model encounters inputs with unfamiliar conditions. Detecting such inputs is vital to preventing models from making naive predictions that may jeopardize real-world applications of neural networks. In this paper, we address the challenging problem of devising a simple yet effective measure of uncertainty in deep neural networks. Specifically, we propose to utilize backpropagated gradients to quantify the uncertainty of trained models. Gradients depict the required amount of change for a model to properly represent given inputs, thus providing a valuable insight into how familiar and certain the model is regarding the inputs. We demonstrate the effectiveness of gradients as a measure of model uncertainty in applications of detecting unfamiliar inputs, including out-of-distribution and corrupted samples. We show that our gradient-based method outperforms state-of-the-art methods by up to 4.8% of AUROC score in out-of-distribution detection and 35.7% in corrupted input detection.


title: A Two Stage Generalized Block Orthogonal Matching Pursuit (TSGBOMP) Algorithm

abstract: Recovery of an unknown sparse signal from a few of its projections is the key objective of compressed sensing. Often one comes across signals that are not ordinarily sparse but are sparse blockwise. Existing block sparse recovery algorithms like BOMP make the assumption of uniform block size and known block boundaries, which are, however, not very practical in many applications. This paper addresses this problem and proposes a two step procedure, where the first stage is a coarse block location identification stage while the second stage carries out finer localization of a non-zero cluster within the window selected in the first stage. A detailed convergence analysis of the proposed algorithm is carried out by first defining the so-called pseudoblock-interleaved block RIP of the given generalized block sparse signal and then imposing upper bounds on the corresponding RIC. We also extend the analysis for complex vector as well as matrix entries where it turns out that the extension is non-trivial and requires special care. Furthermore, assuming real Gaussian sensing matrix entries, we find a lower bound on the probability that the derived recovery bounds are satisfied. The lower bound suggests that there are sets of parameters such that the derived bound is satisfied with high probability. Simulation results confirm significantly improved performance of the proposed algorithm as compared to BOMP.


title: Sampling Multiple Edges Efficiently

abstract: We present a sublinear time algorithm that allows one to sample multiple edges from a distribution that is pointwise $\epsilon$-close to the uniform distribution, in an \emph{amortized-efficient} fashion. We consider the adjacency list query model, where access to a graph $G$ is given via degree and neighbor queries. The problem of sampling a single edge in this model has been raised by Eden and Rosenbaum (SOSA 18). Let $n$ and $m$ denote the number of vertices and edges of $G$, respectively. Eden and Rosenbaum provided upper and lower bounds of $\Theta^*(n/\sqrt m)$ for sampling a single edge in general graphs (where $O^*(\cdot)$ suppresses $\textrm{poly}(1/\epsilon)$ and $\textrm{poly}(\log n)$ dependencies). We ask whether the query complexity lower bound for sampling a single edge can be circumvented when multiple samples are required. That is, can we get an improved amortized per-sample cost if we allow a preprocessing phase? We answer in the affirmative. We present an algorithm that, if one knows the number of required samples $q$ in advance, has an overall cost that is sublinear in $q$, namely, $O^*(\sqrt q \cdot(n/\sqrt m))$, which is strictly preferable to $O^*(q\cdot (n/\sqrt m))$ cost resulting from $q$ invocations of the algorithm by Eden and Rosenbaum. Subsequent to a preliminary version of this work, T\v{e}tek and Thorup (arXiv, preprint) proved that this bound is essentially optimal.


title: Probing Nuclear and Circumnuclear Properties of NGC 6300 using X-ray Observations

abstract: We present the results obtained from a detailed X-ray timing and spectral analysis of Seyfert 2 galaxy NGC 6300 by using observations with the {\it Suzaku}, {\it Chandra} and {\it NuSTAR} observatories between 2007 and 2016. We calculate variance, rms fractional variability of the source in different energy bands and find variabilities in various energy bands. Spectral properties of the source are studied by using various phenomenological and physical models. The properties of the Compton clouds, reflection, Fe K$\alpha$ line emission and soft X-ray excess are studied in detail. Several physical parameters of the source are extracted and investigated to establish the presence/absence of any correlation between them. We also investigate the nature of the circumnuclear `torus' and find that the torus is not uniform, rather clumpy. The observed changes in the line-of-sight column density can be explained in terms of transiting clouds. The iron line emitting region is found to be different in the different epoch of observations. We also observe that the torus and the nucleus independently evolve over the years.


title: Experimental Analysis on Variations and Accuracy of Crosstalk in Trench-Assisted Multi-core Fibers

abstract: Space division multiplexing using multi-core fiber (MCF) is a promising solution to cope with the capacity crunch in standard single-mode fiber based optical communication systems. Nevertheless, the achievable capacity of MCF is limited by inter-core crosstalk (IC-XT). Many existing researches treat IC-XT as a static interference, however, recent research shows that IC-XT varies with time, wavelength and baud rate. This inherent stochastic feature requires a comprehensive characterization of the behaviour of MCF to its application in practical transmission systems and the theoretical understanding of IC-XT phenomenon. In this paper, we experimentally investigate the IC-XT behaviour of an 8-core trench-assisted MCF in a temperature-controlled environment, using popular modulation formats. We compare the measured results with the theoretical prediction to validate the analytical IC-XT models previously developed. Moreover, we explore the effects of the measurement configurations on the IC-XT accuracy and present an analysis on the IC-XT step distribution. Our results indicate that a number of transmission parameters have significant influence on the strength and volatility of IC-XT. Moreover, the averaging time of the power meter and the observation time window can affect the value of the observed IC-XT, the degrees of the effects vary with the type of the source signals.


title: Predicting Coordinated Actuated Traffic Signal Change Times using LSTM Neural Networks

abstract: Vehicle acceleration and deceleration maneuvers at traffic signals results in significant fuel and energy consumption levels. Green light optimal speed advisory systems require reliable estimates of signal switching times to improve vehicle fuel efficiency. Obtaining these estimates is difficult for actuated signals where the length of each green indication changes to accommodate varying traffic conditions. This study details a four-step Long Short-Term Memory deep learning-based methodology that can be used to provide reasonable switching time estimates from green to red and vice versa while being robust to missing data. The four steps are data gathering, data preparation, machine learning model tuning, and model testing and evaluation. The input to the models included controller logic, signal timing parameters, time of day, traffic state from detectors, vehicle actuation data, and pedestrian actuation data. The methodology is applied and evaluated on data from an intersection in Northern Virginia. A comparative analysis is conducted between different loss functions including the mean squared error, mean absolute error, and mean relative error used in LSTM and a new loss function is proposed. The results show that while the proposed loss function outperforms conventional loss functions in terms of overall absolute error values, the choice of the loss function is dependent on the prediction horizon. In particular, the proposed loss function is outperformed by the mean relative error for very short prediction horizons and mean squared error for very long prediction horizons.


title: Short-term origin-destination demand prediction in urban rail transit systems: A channel-wise attentive split-convolutional neural network method

abstract: Short-term origin-destination (OD) flow prediction in urban rail transit (URT) plays a crucial role in smart and real-time URT operation and management. Different from other short-term traffic forecasting methods, the short-term OD flow prediction possesses three unique characteristics: (1) data availability: real-time OD flow is not available during the prediction; (2) data dimensionality: the dimension of the OD flow is much higher than the cardinality of transportation networks; (3) data sparsity: URT OD flow is spatiotemporally sparse. There is a great need to develop novel OD flow forecasting method that explicitly considers the unique characteristics of the URT system. To this end, a channel-wise attentive split-convolutional neural network (CAS-CNN) is proposed. The proposed model consists of many novel components such as the channel-wise attention mechanism and split CNN. In particular, an inflow/outflow-gated mechanism is innovatively introduced to address the data availability issue. We further originally propose a masked loss function to solve the data dimensionality and data sparsity issues. The model interpretability is also discussed in detail. The CAS-CNN model is tested on two large-scale real-world datasets from Beijing Subway, and it outperforms the rest of benchmarking methods. The proposed model contributes to the development of short-term OD flow prediction, and it also lays the foundations of real-time URT operation and management.


title: Moment Multicalibration for Uncertainty Estimation

abstract: We show how to achieve the notion of "multicalibration" from H\'ebert-Johnson et al. [2018] not just for means, but also for variances and other higher moments. Informally, it means that we can find regression functions which, given a data point, can make point predictions not just for the expectation of its label, but for higher moments of its label distribution as well-and those predictions match the true distribution quantities when averaged not just over the population as a whole, but also when averaged over an enormous number of finely defined subgroups. It yields a principled way to estimate the uncertainty of predictions on many different subgroups-and to diagnose potential sources of unfairness in the predictive power of features across subgroups. As an application, we show that our moment estimates can be used to derive marginal prediction intervals that are simultaneously valid as averaged over all of the (sufficiently large) subgroups for which moment multicalibration has been obtained.


title: Fast Approximate Bayesian Contextual Cold Start Learning (FAB-COST)

abstract: Cold-start is a notoriously difficult problem which can occur in recommendation systems, and arises when there is insufficient information to draw inferences for users or items. To address this challenge, a contextual bandit algorithm -- the Fast Approximate Bayesian Contextual Cold Start Learning algorithm (FAB-COST) -- is proposed, which is designed to provide improved accuracy compared to the traditionally used Laplace approximation in the logistic contextual bandit, while controlling both algorithmic complexity and computational cost. To this end, FAB-COST uses a combination of two moment projection variational methods: Expectation Propagation (EP), which performs well at the cold start, but becomes slow as the amount of data increases; and Assumed Density Filtering (ADF), which has slower growth of computational cost with data size but requires more data to obtain an acceptable level of accuracy. By switching from EP to ADF when the dataset becomes large, it is able to exploit their complementary strengths. The empirical justification for FAB-COST is presented, and systematically compared to other approaches on simulated data. In a benchmark against the Laplace approximation on real data consisting of over $670,000$ impressions from autotrader.co.uk, FAB-COST demonstrates at one point increase of over $16\%$ in user clicks. On the basis of these results, it is argued that FAB-COST is likely to be an attractive approach to cold-start recommendation systems in a variety of contexts.


title: Spectral Processing of COVID-19 Time-Series Data

abstract: The presence of oscillations in aggregated COVID-19 data not only raises questions about the data's accuracy, it hinders understanding of the pandemic. Spectral analysis is used to reveal additional properties of the data, and the oscillations are replicated using sinusoidal resynthesis. The precise behavior of the seven-day moving average is also discussed, specifically, the cause of its jaggedness and the phase error it introduces. In comparison, other filtering techniques and Fourier processing produce superior smoothing and have zero phase error. Both of these are presented, and they are extended to isolate several frequency ranges. This extracts some of the same short-term variability that is resynthesized, and it shows that fluctuations with periods between 8 and 21 days are present in U.S. mortality data. These methods have applications that include modeling epidemiological time-series data as well as identifying less obvious properties of the data.


title: The agnostic sampling transceiver

abstract: Increasing capacity demands in the access networks require inventive concepts for the transmission and distribution of digital as well as analog signals over the same network. Here a new transceiver system, which is completely agnostic for the signals to be transmitted is presented. Nyquist sampling and time multiplexing of N phase and intensity modulated digital and analog channels with one single modulator, as well as the transmission and demultiplexing with another modulator have been demonstrated. The aggregate symbol rate corresponds to the modulator bandwidth and can be further increased by a modification of the setup. No high-speed electronic signal processing or high bandwidth photonics is required. Apart from its simplicity and the possibility to process high bandwidth signals with low bandwidth electronics and photonics, the method has the potential to be easily integrated into any platform and thus, might be a solution for the increasing data rates in future access networks.


title: Deep Learning Based on Generative Adversarial and Convolutional Neural Networks for Financial Time Series Predictions

abstract: In the big data era, deep learning and intelligent data mining technique solutions have been applied by researchers in various areas. Forecast and analysis of stock market data have represented an essential role in today's economy, and a significant challenge to the specialist since the market's tendencies are immensely complex, chaotic and are developed within a highly dynamic environment. There are numerous researches from multiple areas intending to take on that challenge, and Machine Learning approaches have been the focus of many of them. There are multiple models of Machine Learning algorithms been able to obtain competent outcomes doing that class of foresight. This paper proposes the implementation of a generative adversarial network (GAN), which is composed by a bi-directional Long short-term memory (LSTM) and convolutional neural network(CNN) referred as Bi-LSTM-CNN to generate synthetic data that agree with existing real financial data so the features of stocks with positive or negative trends can be retained to predict future trends of a stock. The novelty of this proposed solution that distinct from previous solutions is that this paper introduced the concept of a hybrid system (Bi-LSTM-CNN) rather than a sole LSTM model. It was collected data from multiple stock markets such as TSX, SHCOMP, KOSPI 200 and the S&P 500, proposing an adaptative-hybrid system for trends prediction on stock market prices, and carried a comprehensive evaluation on several commonly utilized machine learning prototypes, and it is concluded that the proposed solution approach outperforms preceding models. Additionally, during the research stage from preceding works, gaps were found between investors and researchers who dedicated to the technical domain.


title: Just another quantum assembly language (Jaqal)

abstract: The Quantum Scientific Computing Open User Testbed (QSCOUT) is a trapped-ion quantum computer testbed realized at Sandia National Laboratories on behalf of the Department of Energy's Office of Science and its Advanced Scientific Computing (ASCR) program. Here we describe Jaqal, for Just another quantum assembly language, the programming language we invented to specify programs executed on QSCOUT. Jaqal is useful beyond QSCOUT---it can support mutliple hardware targets because it offloads gate names and their pulse-sequence definitions to external files. We describe the capabilities of the Jaqal language, our approach in designing it, and the reasons for its creation. To learn more about QSCOUT, Jaqal, or JaqalPaq, the metaprogramming Python package we developed for Jaqal, please visit https://qscout.sandia.gov, https://gitlab.com/jaqal, or send an e-mail to qscout@sandia.gov.


title: A class of Finite difference Methods for solving inhomogeneous damped wave equations

abstract: In this paper, a class of finite difference numerical techniques is presented to solve the second-order linear inhomogeneous damped wave equation. The consistency, stability, and convergences of these numerical schemes are discussed. The results obtained are compared to the exact solution, ordinary explicit, implicit finite difference methods, and the fourth-order compact method (FOCM). The general idea of these methods is developed by using the C0-semigroups operator theory. We also showed that the stability region for the explicit finite difference scheme depends on the damping coefficient.


title: Bayesian neural networks and dimensionality reduction

abstract: In conducting non-linear dimensionality reduction and feature learning, it is common to suppose that the data lie near a lower-dimensional manifold. A class of model-based approaches for such problems includes latent variables in an unknown non-linear regression function; this includes Gaussian process latent variable models and variational auto-encoders (VAEs) as special cases. VAEs are artificial neural networks (ANNs) that employ approximations to make computation tractable; however, current implementations lack adequate uncertainty quantification in estimating the parameters, predictive densities, and lower-dimensional subspace, and can be unstable and lack interpretability in practice. We attempt to solve these problems by deploying Markov chain Monte Carlo sampling algorithms (MCMC) for Bayesian inference in ANN models with latent variables. We address issues of identifiability by imposing constraints on the ANN parameters as well as by using anchor points. This is demonstrated on simulated and real data examples. We find that current MCMC sampling schemes face fundamental challenges in neural networks involving latent variables, motivating new research directions.


title: Algorithm Based on One Monocular Video Delivers Highly Valid and Reliable Gait Parameters

abstract: Despite its paramount importance for manifold use cases (e.g., in the health care industry, sports, rehabilitation and fitness assessment), sufficiently valid and reliable gait parameter measurement is still limited to high-tech gait laboratories mostly. Here, we demonstrate the excellent validity and test-retest repeatability of a novel gait assessment system which is built upon modern convolutional neural networks to extract three-dimensional skeleton joints from monocular frontal-view videos of walking humans. The validity study is based on a comparison to the GAITRite pressure-sensitive walkway system. All measured gait parameters (gait speed, cadence, step length and step time) showed excellent concurrent validity for multiple walk trials at normal and fast gait speeds. The test-retest-repeatability is on the same level as the GAITRite system. In conclusion, we are convinced that our results can pave the way for cost, space and operationally effective gait analysis in broad mainstream applications. Most sensor-based systems are costly, must be operated by extensively trained personnel (e.g., motion capture systems) or - even if not quite as costly - still possess considerable complexity (e.g., wearable sensors). In contrast, a video sufficient for the assessment method presented here can be obtained by anyone, without much training, via a smartphone camera.


title: TactileSGNet: A Spiking Graph Neural Network for Event-based Tactile Object Recognition

abstract: Tactile perception is crucial for a variety of robot tasks including grasping and in-hand manipulation. New advances in flexible, event-driven, electronic skins may soon endow robots with touch perception capabilities similar to humans. These electronic skins respond asynchronously to changes (e.g., in pressure, temperature), and can be laid out irregularly on the robot's body or end-effector. However, these unique features may render current deep learning approaches such as convolutional feature extractors unsuitable for tactile learning. In this paper, we propose a novel spiking graph neural network for event-based tactile object recognition. To make use of local connectivity of taxels, we present several methods for organizing the tactile data in a graph structure. Based on the constructed graphs, we develop a spiking graph convolutional network. The event-driven nature of spiking neural network makes it arguably more suitable for processing the event-based data. Experimental results on two tactile datasets show that the proposed method outperforms other state-of-the-art spiking methods, achieving high accuracies of approximately 90\% when classifying a variety of different household objects.


title: A geodesic interior-point method for linear optimization over symmetric cones

abstract: We develop a new interior-point method (IPM) for symmetric-cone optimization, a common generalization of linear, second-order-cone, and semidefinite programming. In contrast to classical IPMs, we update iterates with a geodesic of the cone instead of the kernel of the linear constraints. This approach yields a primal-dual-symmetric, scale-invariant, and line-search-free algorithm that uses just half the variables of a standard primal-dual IPM. With elementary arguments, we establish polynomial-time convergence matching the standard square-root-n bound. Finally, we prove global convergence of a long-step variant and provide an implementation that supports all symmetric cones. For linear programming, our algorithms reduce to central-path tracking in the log domain.


title: Learning Structure in Nested Logit Models

abstract: This paper introduces a new data-driven methodology for nested logit structure discovery. Nested logit models allow the modeling of positive correlations between the error terms of the utility specifications of the different alternatives in a discrete choice scenario through the specification of a nesting structure. Current nested logit model estimation practices require an a priori specification of a nesting structure by the modeler. In this we work we optimize over all possible specifications of the nested logit model that are consistent with rational utility maximization. We formulate the problem of learning an optimal nesting structure from the data as a mixed integer nonlinear programming (MINLP) optimization problem and solve it using a variant of the linear outer approximation algorithm. We exploit the tree structure of the problem and utilize the latest advances in integer optimization to bring practical tractability to the optimization problem we introduce. We demonstrate the ability of our algorithm to correctly recover the true nesting structure from synthetic data in a Monte Carlo experiment. In an empirical illustration using a stated preference survey on modes of transportation in the U.S. state of Massachusetts, we use our algorithm to obtain an optimal nesting tree representing the correlations between the unobserved effects of the different travel mode choices. We provide our implementation as a customizable and open-source code base written in the Julia programming language.


title: Efficient planning of peen-forming patterns via artificial neural networks

abstract: Robust automation of the shot peen forming process demands a closed-loop feedback in which a suitable treatment pattern needs to be found in real-time for each treatment iteration. In this work, we present a method for finding the peen-forming patterns, based on a neural network (NN), which learns the nonlinear function that relates a given target shape (input) to its optimal peening pattern (output), from data generated by finite element simulations. The trained NN yields patterns with an average binary accuracy of 98.8\% with respect to the ground truth in microseconds.


title: The MRS UAV System: Pushing the Frontiers of Reproducible Research, Real-world Deployment, and Education with Autonomous Unmanned Aerial Vehicles

abstract: We present a multirotor Unmanned Aerial Vehicle control (UAV) and estimation system for supporting replicable research through realistic simulations and real-world experiments. We propose a unique multi-frame localization paradigm for estimating the states of a UAV in various frames of reference using multiple sensors simultaneously. The system enables complex missions in GNSS and GNSS-denied environments, including outdoor-indoor transitions and the execution of redundant estimators for backing up unreliable localization sources. Two feedback control designs are presented: one for precise and aggressive maneuvers, and the other for stable and smooth flight with a noisy state estimate. The proposed control and estimation pipeline are constructed without using the Euler/Tait-Bryan angle representation of orientation in 3D. Instead, we rely on rotation matrices and a novel heading-based convention to represent the one free rotational degree-of-freedom in 3D of a standard multirotor helicopter. We provide an actively maintained and well-documented open-source implementation, including realistic simulation of UAV, sensors, and localization systems. The proposed system is the product of years of applied research on multi-robot systems, aerial swarms, aerial manipulation, motion planning, and remote sensing. All our results have been supported by real-world system deployment that shaped the system into the form presented here. In addition, the system was utilized during the participation of our team from the CTU in Prague in the prestigious MBZIRC 2017 and 2020 robotics competitions, and also in the DARPA SubT challenge. Each time, our team was able to secure top places among the best competitors from all over the world. On each occasion, the challenges has motivated the team to improve the system and to gain a great amount of high-quality experience within tight deadlines.


title: On dropping the first Sobol' point

abstract: Quasi-Monte Carlo (QMC) points are a substitute for plain Monte Carlo (MC) points that greatly improve integration accuracy under mild assumptions on the problem. Because QMC can give errors that are $o(1/n)$ as $n\to\infty$, changing even one point can change the estimate by an amount much larger than the error would have been and worsen the convergence rate. As a result, certain practices that fit quite naturally and intuitively with MC points are very detrimental to QMC performance. These include thinning, burn-in, and taking sample sizes such as powers of $10$, other than the ones for which the QMC points were designed. This article looks at the effects of a common practice in which one skips the first point of a Sobol' sequence. The retained points ordinarily fail to be a digital net and when scrambling is applied, skipping over the first point can increase the numerical error by a factor proportional to $\sqrt{n}$ where $n$ is the number of function evaluations used.