title: Learning Inverse Kinodynamics for Accurate High-Speed Off-Road Navigation on Unstructured Terrain

abstract: This paper presents a learning-based approach to consider the effect of unobservable world states in kinodynamic motion planning in order to enable accurate high-speed off-road navigation on unstructured terrain. Existing kinodynamic motion planners either operate in structured and homogeneous environments and thus do not need to explicitly account for terrain-vehicle interaction, or assume a set of discrete terrain classes. However, when operating on unstructured terrain, especially at high speeds, even small variations in the environment will be magnified and cause inaccurate plan execution. In this paper, to capture the complex kinodynamic model and mathematically unknown world state, we learn a kinodynamic planner in a data-driven manner with onboard inertial observations. Our approach is tested on a physical robot in different indoor and outdoor environments, enables fast and accurate off-road navigation, and outperforms environment-independent alternatives, demonstrating 52.4% to 86.9% improvement in terms of plan execution success rate while traveling at high speeds.


title: Learning-based Robust Motion Planning with Guaranteed Stability: A Contraction Theory Approach

abstract: This paper presents Learning-based Autonomous Guidance with RObustness and Stability guarantees (LAG-ROS), which provides machine learning-based nonlinear motion planners with formal robustness and stability guarantees, by designing a differential Lyapunov function using contraction theory. LAG-ROS utilizes a neural network to model a robust tracking controller independently of a target trajectory, for which we show that the Euclidean distance between the target and controlled trajectories is exponentially bounded linearly in the learning error, even under the existence of bounded external disturbances. We also present a convex optimization approach that minimizes the steady-state bound of the tracking error to construct the robust control law for neural network training. In numerical simulations, it is demonstrated that the proposed method indeed possesses superior properties of robustness and nonlinear stability resulting from contraction theory, whilst retaining the computational efficiency of existing learning-based motion planners.


title: ISALT: Inference-based schemes adaptive to large time-stepping for locally Lipschitz ergodic systems

abstract: Efficient simulation of SDEs is essential in many applications, particularly for ergodic systems that demand efficient simulation of both short-time dynamics and large-time statistics. However, locally Lipschitz SDEs often require special treatments such as implicit schemes with small time-steps to accurately simulate the ergodic measure. We introduce a framework to construct inference-based schemes adaptive to large time-steps (ISALT) from data, achieving a reduction in time by several orders of magnitudes. The key is the statistical learning of an approximation to the infinite-dimensional discrete-time flow map. We explore the use of numerical schemes (such as the Euler-Maruyama, a hybrid RK4, and an implicit scheme) to derive informed basis functions, leading to a parameter inference problem. We introduce a scalable algorithm to estimate the parameters by least squares, and we prove the convergence of the estimators as data size increases. We test the ISALT on three non-globally Lipschitz SDEs: the 1D double-well potential, a 2D multi-scale gradient system, and the 3D stochastic Lorenz equation with degenerate noise. Numerical results show that ISALT can tolerate time-step magnitudes larger than plain numerical schemes. It reaches optimal accuracy in reproducing the invariant measure when the time-step is medium-large.


title: Real-Time Ellipse Detection for Robotics Applications

abstract: We propose a new algorithm for real-time detection and tracking of elliptic patterns suitable for real-world robotics applications. The method fits ellipses to each contour in the image frame and rejects ellipses that do not yield a good fit. The resulting detection and tracking method is lightweight enough to be used on robots' resource-limited onboard computers, can deal with lighting variations and detect the pattern even when the view is partial. The method is tested on an example application of an autonomous UAV landing on a fast-moving vehicle to show its performance indoors, outdoors, and in simulation on a real-world robotics task. The comparison with other well-known ellipse detection methods shows that our proposed algorithm outperforms other methods with the F1 score of 0.981 on a dataset with over 1500 frames. The videos of experiments, the source codes, and the collected dataset are provided with the paper at https://theairlab.org/landing-on-vehicle .


title: LET: Linguistic Knowledge Enhanced Graph Transformer for Chinese Short Text Matching

abstract: Chinese short text matching is a fundamental task in natural language processing. Existing approaches usually take Chinese characters or words as input tokens. They have two limitations: 1) Some Chinese words are polysemous, and semantic information is not fully utilized. 2) Some models suffer potential issues caused by word segmentation. Here we introduce HowNet as an external knowledge base and propose a Linguistic knowledge Enhanced graph Transformer (LET) to deal with word ambiguity. Additionally, we adopt the word lattice graph as input to maintain multi-granularity information. Our model is also complementary to pre-trained language models. Experimental results on two Chinese datasets show that our models outperform various typical text matching approaches. Ablation study also indicates that both semantic information and multi-granularity information are important for text matching modeling.


title: A Layered Grouping Random Access Scheme Based on Dynamic Preamble Selection for Massive Machine Type Communications

abstract: Massive machine type communication (mMTC) has been identified as an important use case in Beyond 5G networks and future massive Internet of Things (IoT). However, for the massive multiple access in mMTC, there is a serious access preamble collision problem if the conventional 4-step random access (RA) scheme is employed. Consequently, a range of grantfree (GF) RA schemes were proposed. Nevertheless, if the number of cellular users (devices) significantly increases, both the energy and spectrum efficiency of the existing GF schemes still rapidly degrade owing to the much longer preambles required. In order to overcome this dilemma, a layered grouping strategy is proposed, where the cellular users are firstly divided into clusters based on their geographical locations, and then the users of the same cluster autonomously join in different groups by using optimum energy consumption (Opt-EC) based K-means algorithm. With this new layered cellular architecture, the RA process is divided into cluster load estimation phase and active group detection phase. Based on the state evolution theory of approximated message passing algorithm, a tight lower bound on the minimum preamble length for achieving a certain detection accuracy is derived. Benefiting from the cluster load estimation, a dynamic preamble selection (DPS) strategy is invoked in the second phase, resulting the required preambles with minimum length. As evidenced in our simulation results, this two-phase DPS aided RA strategy results in a significant performance improvement


title: Phonon heat conduction in Al1-xScxN thin films

abstract: Aluminum scandium nitride alloy (Al1-xScxN) is regarded as a promising material for high-performance acoustic devices used in wireless communication systems. Phonon scattering and heat conduction processes govern the energy dissipation in acoustic resonators, ultimately determining their performance quality. This work reports, for the first time, on phonon scattering processes and thermal conductivity in Al1-xScxN alloys with the Sc content (x) up to 0.26. The thermal conductivity measured presents a descending trend with increasing x. Temperature-dependent measurements show an increase in thermal conductivity as the temperature increases at temperatures below 200K, followed by a plateau at higher temperatures (T> 200K). Application of a virtual crystal phonon conduction model allows us to elucidate the effects of boundary and alloy scattering on the observed thermal conductivity behaviors. We further demonstrate that the alloy scattering is caused mainly by strain-field difference, and less by the atomic mass difference between ScN and AlN, which is in contrast to the well-studied Al1-xGaxN and SixGe1-x alloy systems where atomic mass difference dominates the alloy scattering. This work studies and provides the quantitative knowledge for phonon scattering and the thermal conductivity in Al1-xScxN, paving the way for future investigation of materials and design of acoustic devices.


title: Gemini/Phoenix H-band analysis of the globular cluster AL3

abstract: The globular cluster AL~3 is old and located in the inner bulge. Three individual stars were observed with the Phoenix spectrograph at the Gemini South telescope. The wavelength region contains prominent lines of CN, OH, and CO, allowing the derivation of C, N, and O abundances of cool stars. We aim to derive C, N, O abundances of three stars in the bulge globular cluster AL3, and additionally in stars of NGC 6558 and HP1. The spectra of AL3 allows us to derive the cluster's radial velocity. For AL3, we applied a new code to analyse its colour-magnitude diagram. Synthetic spectra were computed and compared to observed spectra for the three clusters. We present a detailed identification of lines in the spectral region centred at 15555 A, covering the wavelength range 15525-15590 A. C, N, and O abundances are tentatively derived for the sample stars.


title: Computing Accurate Probabilistic Estimates of One-D Entropy from Equiprobable Random Samples

abstract: We develop a simple Quantile Spacing (QS) method for accurate probabilistic estimation of one-dimensional entropy from equiprobable random samples, and compare it with the popular Bin-Counting (BC) method. In contrast to BC, which uses equal-width bins with varying probability mass, the QS method uses estimates of the quantiles that divide the support of the data generating probability density function (pdf) into equal-probability-mass intervals. Whereas BC requires optimal tuning of a bin-width hyper-parameter whose value varies with sample size and shape of the pdf, QS requires specification of the number of quantiles to be used. Results indicate, for the class of distributions tested, that the optimal number of quantile-spacings is a fixed fraction of the sample size (empirically determined to be ~0.25-0.35), and that this value is relatively insensitive to distributional form or sample size, providing a clear advantage over BC since hyperparameter tuning is not required. Bootstrapping is used to approximate the sampling variability distribution of the resulting entropy estimate, and is shown to accurately reflect the true uncertainty. For the four distributional forms studied (Gaussian, Log-Normal, Exponential and Bimodal Gaussian Mixture), expected estimation bias is less than 1% and uncertainty is relatively low even for very small sample sizes. We speculate that estimating quantile locations, rather than bin-probabilities, results in more efficient use of the information in the data to approximate the underlying shape of an unknown data generating pdf.


title: Efficient computational algorithms for approximate optimal designs

abstract: In this paper, we propose two simple yet efficient computational algorithms to obtain approximate optimal designs for multi-dimensional linear regression on a large variety of design spaces. We focus on the two commonly used optimal criteria, $D$- and $A$-optimal criteria. For $D$-optimality, we provide an alternative proof for the monotonic convergence for $D$-optimal criterion and propose an efficient computational algorithm to obtain the approximate $D$-optimal design. We further show that the proposed algorithm converges to the $D$-optimal design, and then prove that the approximate $D$-optimal design converges to the continuous $D$-optimal design under certain conditions. For $A$-optimality, we provide an efficient algorithm to obtain approximate $A$-optimal design and conjecture the monotonicity of the proposed algorithm. Numerical comparisons suggest that the proposed algorithms perform well and they are comparable or superior to some existing algorithms.


title: Do Not Let Privacy Overbill Utility: Gradient Embedding Perturbation for Private Learning

abstract: The privacy leakage of the model about the training data can be bounded in the differential privacy mechanism. However, for meaningful privacy parameters, a differentially private model degrades the utility drastically when the model comprises a large number of trainable parameters. In this paper, we propose an algorithm \emph{Gradient Embedding Perturbation (GEP)} towards training differentially private deep models with decent accuracy. Specifically, in each gradient descent step, GEP first projects individual private gradient into a non-sensitive anchor subspace, producing a low-dimensional gradient embedding and a small-norm residual gradient. Then, GEP perturbs the low-dimensional embedding and the residual gradient separately according to the privacy budget. Such a decomposition permits a small perturbation variance, which greatly helps to break the dimensional barrier of private learning. With GEP, we achieve decent accuracy with reasonable computational cost and modest privacy guarantee for deep models. Especially, with privacy bound $\epsilon=8$, we achieve $74.9\%$ test accuracy on CIFAR10 and $95.1\%$ test accuracy on SVHN, significantly improving over existing results.


title: Deep learning polarization distributions in ferroelectrics from STEM data: with and without atom finding

abstract: Over the last decade, scanning transmission electron microscopy (STEM) has emerged as a powerful tool for probing atomic structures of complex materials with picometer precision, opening the pathway toward exploring ferroelectric, ferroelastic, and chemical phenomena on the atomic-scale. Analyses to date extracting a polarization signal from lattice coupled distortions in STEM imaging rely on discovery of atomic positions from intensity maxima/minima and subsequent calculation of polarization and other order parameter fields from the atomic displacements. Here, we explore the feasibility of polarization mapping directly from the analysis of STEM images using deep convolutional neural networks (DCNNs). In this approach, the DCNN is trained on the labeled part of the image (i.e., for human labelling), and the trained network is subsequently applied to other images. We explore the effects of the choice of the descriptors (centered on atomic columns and grid-based), the effects of observational bias, and whether the network trained on one composition can be applied to a different one. This analysis demonstrates the tremendous potential of the DCNN for the analysis of high-resolution STEM imaging and spectral data and highlights the associated limitations.


title: Variational Selective Autoencoder: Learning from Partially-Observed Heterogeneous Data

abstract: Learning from heterogeneous data poses challenges such as combining data from various sources and of different types. Meanwhile, heterogeneous data are often associated with missingness in real-world applications due to heterogeneity and noise of input sources. In this work, we propose the variational selective autoencoder (VSAE), a general framework to learn representations from partially-observed heterogeneous data. VSAE learns the latent dependencies in heterogeneous data by modeling the joint distribution of observed data, unobserved data, and the imputation mask which represents how the data are missing. It results in a unified model for various downstream tasks including data generation and imputation. Evaluation on both low-dimensional and high-dimensional heterogeneous datasets for these two tasks shows improvement over state-of-the-art models.


title: Confidence Calibration with Bounded Error Using Transformations

abstract: As machine learning techniques become widely adopted in new domains, especially in safety-critical systems such as autonomous vehicles, it is crucial to provide accurate output uncertainty estimation. As a result, many approaches have been proposed to calibrate neural networks to accurately estimate the likelihood of misclassification. However, while these methods achieve low calibration error, there is space for further improvement, especially in large-dimensional settings such as ImageNet. In this paper, we introduce a calibration algorithm, named Hoki, that works by applying random transformations to the neural network logits. We provide a sufficient condition for perfect calibration based on the number of label prediction changes observed after applying the transformations. We perform experiments on multiple datasets and show that the proposed approach generally outperforms state-of-the-art calibration algorithms across multiple datasets and models, especially on the challenging ImageNet dataset. Finally, Hoki is scalable as well, as it requires comparable execution time to that of temperature scaling.


title: Transition Density of an Infinite-dimensional diffusion with the Jack Parameter

abstract: From the Poisson-Dirichlet diffusions to the $Z$-measure diffusions, they all have explicit transition densities. In this paper, we will show that the transition densities of the $Z$-measure diffusions can also be expressed as a mixture of a sequence of probability measures on the Thoma simplex. The coefficients are still the transition probabilities of the Kingman coalescent stopped at state $1$. This fact will be uncovered by a dual process method in a special case where the $Z$-measure diffusions is established through up-down chain in the Young graph.


title: Aximorphic Perspective Projection Model for Immersive Imagery

abstract: A wide choice of cinematic lenses enables motion-picture creators to adapt image visual-appearance to their creative vision. Such choice does not exist in the realm of real-time computer graphics, where only one type of perspective projection is widely used, a linear perspective. This paper presents an extended perspective imaging model, which can represent distortion and FoV parameters of entire variety of film and photographic lenses (e.g., wide-angle, fisheye, anamorphic), while preserving parametrization in an artistically convincing manner. Self-experimentation with the model revealed that each projection type provides accurate perception of a different aspect of depicted space (e.g., speed, distance, shape). Presented model, enables combination of multiple projections, each on a different axis of the image, to achieve optimal perception for a given scenario. This new projection, named aximorphic, was made available here, under an open license (CC BY-SA 3.0), for a wide and easy adoption.


title: Convex geometric reasoning for crystalline energies

abstract: The present work revisits the classical Wulff problem restricted to crystalline integrands, a class of surface energies that gives rise to finitely faceted crystals. The general proof of the Wulff theorem was given by J.E. Taylor (1978) by methods of Geometric Measure Theory. This work follows a simpler and direct way through Minkowski Theory by taking advantage of the convex properties of the considered Wulff shapes.


title: Ensuring Progress for Multiple Mobile Robots via Space Partitioning, Motion Rules, and Adaptively Centralized Conflict Resolution

abstract: In environments where multiple robots must coordinate in a shared space, decentralized approaches allow for decoupled planning at the cost of global guarantees, while centralized approaches make the opposite trade-off. These solutions make a range of assumptions - commonly, that all the robots share the same planning strategies. In this work, we present a framework that ensures progress for all robots without assumptions on any robot's planning strategy by (1) generating a partition of the environment into "flow", "open", and "passage" regions and (2) imposing a set of rules for robot motion in these regions. These rules for robot motion prevent deadlock through an adaptively centralized protocol for resolving spatial conflicts between robots. Our proposed framework ensures progress for all robots without a grid-like discretization of the environment or strong requirements on robot communication, coordination, or cooperation. Each robot can freely choose how to plan and coordinate for itself, without being vulnerable to other robots or groups of robots blocking them from their goals, as long as they follow the rules when necessary. We describe our space partition and motion rules, prove that the motion rules suffice to guarantee progress in partitioned environments, and demonstrate several cases in simulated polygonal environments. This work strikes a balance between each robot's planning independence and a guarantee that each robot can always reach any goal in finite time.


title: A Local Method for Identifying Causal Relations under Markov Equivalence

abstract: Causality is important for designing interpretable and robust methods in artificial intelligence research. We propose a local approach to identify whether a variable is a cause of a given target under the framework of causal graphical models of directed acyclic graphs (DAGs). In general, the causal relation between two variables may not be identifiable from observational data as many causal DAGs encoding different causal relations are Markov equivalent. In this paper, we first introduce a sufficient and necessary graphical condition to check the existence of a causal path from a variable to a target in every Markov equivalent DAG. Next, we provide local criteria for identifying whether a variable is a cause/non-cause of a target based only on the local structure instead of the entire graph. Finally, we propose a local learning algorithm for this causal query via learning the local structure of the variable and some additional statistical independence tests related to the target. Simulation studies show that our local algorithm is efficient and effective, compared with other state-of-art methods.


title: Radiative cooling of colored paint based on Fe3+ doped Y2Ce2O7

abstract: Materials with both low absorption of incoming solar radiation and high emittance in mid-infrared band can be applied for daytime radiative cooling. Current state-of-the-art materials for passive radiative cooling often utilize a combination of solar reflector and infrared emitter by different structures, or even by expensive nanofabricated photonic structures, which limits the applications in practice. In this study, possessing these two specified radiative properties, pure Y2Ce2O7 is demonstrated with a performance of passive radiative cooling. With a bandgap at 375.7 nm, the prepared Y2Ce2O7 shows a high solar reflectance of 91%, while with lattice strain and distortion of various bonds (e.g., Y-O, Ce-O), it also shows a high emittance of 0.96 in MIR band. More attracting, the aesthetics performance of Y2Ce2O7 can be modified by doping Fe3+ ions to change its color from ivory white to light yellow or red with high NIR reflection and MIR emission, indicating that the Y2Ce2-xFexO7 shows a better cooling performance than a common paint with a similar color. According to the field demonstration of cooling performance at noon time, the Y2Ce2O7 and Y2Ce1.9Fe0.1O7 paints are 2.2 K and 1.8 K lower than the common white and umber paints, respectively, while at night, all paints are 2.3 K lower than the ambient air. If applied on the envelop of a building, the simulation shows that the Y2Ce2O7 and Y2Ce1.9Fe0.1O7 paints save 54.45% and 21.14% energy consumption compared with a common white and umber paints, respectively, in a hot season. The demonstrated Y2Ce2-xFexO7 holds potentials for energy-saving applications in hot climates.


title: A flapping feathered wing-powered aerial vehicle

abstract: An aerial vehicle powered by flapping feathered wings was designed, developed and fabricated. Different from legacy flapping-wing aerial vehicles with membrane wings, the new design uses authentic bird feathers to fabricate wings. In field tests, a radio-controlled electric-powered aerial vehicle with flapping feathered wings successfully took off, flew up to 63.88 s and landed safely. It was found that flapping feathered wings can generate sufficient thrust and lift to make a man-made aerial vehicle accomplish takeoff, sustainable flight and a safe landing.


title: Summation formulas of $q$-hyperharmonic numbers

abstract: In this paper, several weighted summation formulas of $q$-hyperharmonic numbers are derived. As special cases, several formulas of hyperharmonic numbers of type $\sum_{\ell=1}^{n} {\ell}^{p} H_{\ell}^{(r)}$ and $\sum_{\ell=0}^{n} {\ell}^{p} H_{n-\ell}^{(r)}$ are obtained.


title: Campana points, Height zeta functions, and log Manin's conjecture

abstract: This is a report of the author's talk at RIMS workshop 2020 Problems and Prospects in Analytic Number Theory held online on Zoom. We discuss a recent formulation of log Manin's conjecture for klt Campana points and an approach to this conjecture using the height zeta function method.


title: The Control of the Elementary Quantum Systems Radiation Using Metamaterials and Nanometaparticles

abstract: The most important direction in the development of fundamental and applied physics is the study of the properties of optical systems at the nanoscale in order to create optical and quantum computers, biosensors, single-photon sources for quantum informatics, devices for DNA sequencing, sensors of various fields, etc. In all these cases, nanoscale light sources - dye molecules, quantum dots (epitaxial or colloidal), color centers in crystals, and nanocontacts in metals - are of key importance. In the nanoenvironment, the characteristics of these elementary quantum systems - pumping rates, radiative and non-radiative decay rates, the local density of states, lifetimes, level shifts - experience changes that can be used intentionally to create nanoscale light sources with desired properties. This review presents an analysis of actual theoretical and experimental works in the field of elementary quantum systems radiation control using plasmonic and dielectric nanostructures, metamaterials, and nanoparticles made from metamaterials.


title: Weyl Semimetal Path to Valley Filtering in Graphene

abstract: We propose a device in which a sheet of graphene is coupled to a Weyl semimetal, allowing for the physical access to the study of tunneling from two-dimensional to three dimensional massless Dirac fermions. Due to the reconstructed band structure, we find that this device acts as a robust valley filter for electrons in the graphene sheet. We show that, by appropriate alignment, the Weyl semimetal draws away current in one of the two graphene valleys while allowing current in the other to pass unimpeded. In contrast to other proposed valley filters, the mechanism of our proposed device occurs in the bulk of the graphene sheet, obviating the need for carefully shaped edges or dimensions.


title: New identities for the Shannon function and applications

abstract: We show how the Shannon entropy function H(p,q)is expressible as a linear combination of other Shannon entropy functions involving quotients of polynomials in p,q of degree n for any given positive integer n. An application to cryptographic keys is presented.


title: Decay properties of Axially Symmetric D-solutions to the Steady Incompressible Magnetohydrodynamic Equations

abstract: In this paper, we investigate the decay properties of axially symmetric solutions to the steady incompressible magnetohydrodynamic equations in $\mathbb{R}^3$ with finite Dirichlet integral. We first derive the decay rates of general D-solutions to the axisymmetric MHD equations. In the special case where the magnetic field only has the swirl component ${\textbf {h}}(r,z)= h_{\theta}(r,z) {\mathbf e_{\theta}}$, we obtain better decay rates. The last result examines the decay rates along the axis $Oz$ also within the special class of D-solutions with only swirl magnetic field. The main tool in this paper is the combination of the scaling argument, the \textit{Brezis-Gallouet inequality} and the weighted energy estimate.


title: Deep Equal Risk Pricing of Financial Derivatives with Multiple Hedging Instruments

abstract: This paper studies the equal risk pricing (ERP) framework for the valuation of European financial derivatives. This option pricing approach is consistent with global trading strategies by setting the premium as the value such that the residual hedging risk of the long and short positions in the option are equal under optimal hedging. The ERP setup of Marzban et al. (2020) is considered where residual hedging risk is quantified with convex risk measures. The main objective of this paper is to assess through extensive numerical experiments the impact of including options as hedging instruments within the ERP framework. The reinforcement learning procedure developed in Carbonneau and Godin (2020), which relies on the deep hedging algorithm of Buehler et al. (2019b), is applied to numerically solve the global hedging problems by representing trading policies with neural networks. Among other findings, numerical results indicate that in the presence of jump risk, hedging long-term puts with shorter-term options entails a significant decrease of both equal risk prices and market incompleteness as compared to trading only the stock. Monte Carlo experiments demonstrate the potential of ERP as a fair valuation approach providing prices consistent with observable market prices. Analyses exhibit the ability of ERP to span a large interval of prices through the choice of convex risk measures which is close to encompass the variance-optimal premium.


title: Learning orbital dynamics of binary black hole systems from gravitational wave measurements

abstract: We introduce a gravitational waveform inversion strategy that discovers mechanical models of binary black hole (BBH) systems. We show that only a single time series of (possibly noisy) waveform data is necessary to construct the equations of motion for a BBH system. Starting with a class of universal differential equations parameterized by feed-forward neural networks, our strategy involves the construction of a space of plausible mechanical models and a physics-informed constrained optimization within that space to minimize the waveform error. We apply our method to various BBH systems including extreme and comparable mass ratio systems in eccentric and non-eccentric orbits. We show the resulting differential equations apply to time durations longer than the training interval, and relativistic effects, such as perihelion precession, radiation reaction, and orbital plunge, are automatically accounted for. The methods outlined here provide a new, data-driven approach to studying the dynamics of binary black hole systems.


title: Resolving ultrahigh-contrast ultrashort pulses with single-shot cross-correlator at the photon noise limit

abstract: In strong-field physics experiments with intense lasers, it is of paramount importance to single-shot diagnose the temporal contrast between laser pulse peak and its noise pedestal. This allows fast optimization of pulse contrast and meaningful comparison with theory for each pulse shot, and it can help new outcomes from clean laser-plasma interactions. Thus far, high contrast ratios up to ~10^10, required by present petawatt (PW) class lasers, have been accessible in both generation and single-shot characterization. However, ultrahigh contrast ~10^13, required by the planned 200-PW lasers, challenges intense laser technology and remains an open question. This paper reports on the first demonstration of such an ultrahigh-contrast measurement by adapting single-shot cross-correlator (SSCC). We introduce an ultrafast method that enables to determine the SSCC detection limit. Our strategy mimics the test laser having known ultrahigh contrast in the measurement frame of time-to-space mapping. The ultimate contrast-measurement limit of 10^13 is achieved, which corresponds to the highest pulse intensity set by SSCC damage threshold and the lowest noise pedestal set by single-photon detection. As a consequence, photon noise in the detection is observed and increases as the noise pedestal reduces. The demonstrated measurement ability at the photon noise limit is applied to a high-contrast laser system based on second-harmonic generation and optical parametric chirped-pulse amplification, suggesting accessible of ultrahigh contrast pulses.


title: Strapdown Inertial Navigation System Initial Alignment based on Group of Double Direct Spatial Isometries

abstract: The task of strapdown inertial navigation system (SINS) initial alignment is to calculate the attitude transformation matrix from body frame to navigation frame. In this paper, such attitude transformation matrix is divided into two parts through introducing the initial inertially fixed navigation frame as inertial frame. The attitude changes of the navigation frame corresponding to the defined inertial frame can be exactly calculated with known velocity and position provided by GNSS. The attitude from body frame to the defined inertial frame is estimated based on the SINS mechanization in inertial frame. The attitude, velocity and position in inertial frame are formulated together as element of the group of double direct spatial isometries.It is proven that the group state model in inertial frame satisfies a particular "group affine" property and the corresponding error model satisfies a "log linear" autonomous differential equation on the Lie algebra. Based on such striking property, the attitude from body frame to the defined inertial frame can be estimated based on the linear error model with even extreme large misalignments. Two different error state vectors are extracted based on right and left matrix multiplications and the detailed linear state space models are derived based on the right and left errors for SINS mechanization in inertial frame. With the derived linear state space models, the explicit initial alignment procedures have been presented. Extensive simulation and field tests indicate that the initial alignment based on the left error model can perform quite well within a wide range of initial attitude errors, although the used filter is still a type of linear Kalman filter. This method is promising in practical products abandoning the traditional coarse alignment stage.


title: Improving the Hosmer-Lemeshow Goodness-of-Fit Test in Large Models with Replicated Trials

abstract: The Hosmer-Lemeshow (HL) test is a commonly used global goodness-of-fit (GOF) test that assesses the quality of the overall fit of a logistic regression model. In this paper, we give results from simulations showing that the type 1 error rate (and hence power) of the HL test decreases as model complexity grows, provided that the sample size remains fixed and binary replicates are present in the data. We demonstrate that the generalized version of the HL test by Surjanovic et al. (2020) can offer some protection against this power loss. We conclude with a brief discussion explaining the behaviour of the HL test, along with some guidance on how to choose between the two tests.


title: File fragment recognition based on content and statistical features

abstract: Nowadays, the speed up development and use of digital devices such as smartphones have put people at risk of internet crimes. The evidence of present crimes in a computer file can be easily unreachable by changing the prefix of a file or other algorithms. In more complex cases, either file divided into different parts or the parts of a file that has information about the file type are deleted, where the file fragment recognition issue is discussed. The known files are divided into different fragments, and different classification algorithms are used to solve the problems of file fragment recognition. The issue of identifying the type of file fragment due to its importance in cybercrime issues as well as antivirus has been highly emphasized and has been addressed in many articles. Increasing the accuracy in this field on the types of widely used files due to the sensitivity of the subject of recognizing the type of file under study is the main goal of researchers in this field. Failure to identify the correct type of file will lead to deviations of the results and evidence from the main issue or failure to conclude. In this paper, first, the file is divided into different fragments. Then, the file fragment features, which are obtained from Binary Frequency Distribution, are reduced by 2 feature reduction algorithms; Sequential Forward Selection algorithm as well as Sequential Floating Forward Selection algorithm to delete sparse features that result in increased accuracy and speed. Finally, the reduced features are given to 3 Multiclass classifier algorithms, Multilayer Perceptron, Support Vector Machines, and K-Nearest Neighbor for classification and comparison of the results. The proposed recognition algorithm can recognize 6 types of useful files and may distinguish a type of file fragments with higher accuracy than the similar works done.


title: Sentiment Analysis of Persian-English Code-mixed Texts

abstract: The rapid production of data on the internet and the need to understand how users are feeling from a business and research perspective has prompted the creation of numerous automatic monolingual sentiment detection systems. More recently however, due to the unstructured nature of data on social media, we are observing more instances of multilingual and code-mixed texts. This development in content type has created a new demand for code-mixed sentiment analysis systems. In this study we collect, label and thus create a dataset of Persian-English code-mixed tweets. We then proceed to introduce a model which uses BERT pretrained embeddings as well as translation models to automatically learn the polarity scores of these Tweets. Our model outperforms the baseline models that use Na\"ive Bayes and Random Forest methods.


title: Dimension of divergence set of the wave equation

abstract: We consider the Hausdorff dimension of the divergence set on which the pointwise convergence $\lim_{t\rightarrow 0} e^{it\sqrt{-\Delta}} f(x) = f(x)$ fails when $f \in H^s(\mathbb R^d)$. We especially prove the conjecture raised by Barcel\'o, Bennett, Carbery and Rogers \cite{BBCR} for $d=3$, and improve the previous results in higher dimensions $d\ge4$. We also show that a Strichartz type estimate for $f\to e^{it\sqrt{-\Delta}} f$ with the measure $ dt\,d\mu(x)$ is essentially equivalent to the estimate for the spherical average of $\widehat \mu$ which has been extensively studied for the Falconer distance set problem. The equivalence provides shortcuts to the recent results due to B. Liu and K. Rogers.


title: LazyFormer: Self Attention with Lazy Update

abstract: Improving the efficiency of Transformer-based language pre-training is an important task in NLP, especially for the self-attention module, which is computationally expensive. In this paper, we propose a simple but effective solution, called \emph{LazyFormer}, which computes the self-attention distribution infrequently. LazyFormer composes of multiple lazy blocks, each of which contains multiple Transformer layers. In each lazy block, the self-attention distribution is only computed once in the first layer and then is reused in all upper layers. In this way, the cost of computation could be largely saved. We also provide several training tricks for LazyFormer. Extensive experiments demonstrate the effectiveness of the proposed method.


title: Externally corrected CCSD with renormalized perturbative triples (R-ecCCSD(T)) and density matrix renormalization group and selected configuration interaction external sources

abstract: We investigate the renormalized perturbative triples correction together with the externally corrected coupled-cluster singles and doubles (ecCCSD) method. We take the density matrix renormalization group (DMRG) and heatbath CI (HCI) as external sources for the ecCCSD equations. The accuracy is assessed for the potential energy surfaces of H2O, N2, and F2. We find that the triples correction significantly improves on ecCCSD and we do not see any instability of the renormalized triples with respect to dissociation. We explore how to balance the cost of computing the external source amplitudes with respect to the accuracy of the subsequent CC calculation. In this context, we find that very approximate wavefunctions (and their large amplitudes) serve as an efficient and accurate external source. Finally, we characterize the domain of correlation treatable using the externally corrected method and renormalized triples combination studied in this work via a well-known wavefunction diagnostic.


title: Collective Bias Models in Two-Tier Voting Systems and the Democracy Deficit

abstract: We analyse optimal voting weights in two-tier voting systems. In our model, the overall population (or union) is split in groups (or member states) of different sizes. The individuals comprising the overall population constitute the first tier, and the council is the second tier. Each group has a representative in the council that casts votes on their behalf. By "optimal weights", we mean voting weights in the council which minimise the democracy deficit, i.e. the expected deviation of the council vote from a (hypothetical) popular vote. We assume that the voters within each group interact via what we call a local collective bias or common belief (through tradition, common values, strong religious beliefs, etc.). We allow in addition an interaction across group borders via a global bias. Thus, the voting behaviour of each voter depends on the behaviour of all other voters. This correlation may be stronger between voters in the same group, but is in general not zero for voters in different groups. We call the respective voting measure a Collective Bias Model (CBM). The "simple CBM" introduced in [12] and in particular the Impartial Culture and the Impartial Anonymous Culture are special cases of our general model. We compute the optimal weights in the large population limit. Those optimal weights are unique as long as there is no "complete" correlation between the groups. In this case, we obtain optimal weights which are the sum of a common constant equal for all groups and a summand which is proportional to the population of each group. We also analyse the conditions under which the optimal weights are negative, thus making it impossible to reach the theoretical minimum of the democracy deficit. This is a new aspect of the model owed to the correlation between votes belonging to different groups.


title: Ab initio MCDHF calculations of the In and Tl electron affinities and their isotope shifts

abstract: We report multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction calculations on the Thallium (Tl) electron affinity, as well as on the excited energy levels arising from the ground configuration of Tl$^-$. The results are compared with the available experimental values and further validated by extending the study to its homologous, lighter element, Indium (In), belonging to Group 13 (III.A) of the periodic table. The calculated electron affinities of In and Tl, 383.4 and 322.8 meV, agree with the latest measurements by within 1\%. Three bound states $^3P_{0,1,2}$ are confirmed in the $5s^25p^2$ configuration of In$^-$ while only the ground state $^3P_{0}$ is bound in the $6s^26p^2$ configuration of Tl$^-$. The isotope shifts on the In and Tl electron affinities are also estimated. The E2/M1 intraconfiguration radiative transition rates within $5s^25p^2 \; ^3P_{0,1,2}$ of In$^-$ are used to calculate the radiative lifetimes of the metastable $^3P_{1,2}$ levels.


title: Schr\"odinger operators with complex sparse potentials

abstract: We establish quantitative upper and lower bounds for Schr\"odinger operators with complex potentials that satisfy some weak form of sparsity. Our first result is a quantitative version of an example, due to S.\ Boegli (Comm. Math. Phys., 2017, 352, 629-639), of a Schr\"odinger operator with eigenvalues accumulating to every point of the essential spectrum. The second result shows that the eigenvalue bounds of Frank (Bull. Lond. Math. Soc., 2011, 43, 745-750 and Trans. Amer. Math. Soc., 2018, 370, 219-240) can be improved for sparse potentials. The third result generalizes a theorem of Klaus (Ann. Inst. H. Poincar\'e Sect. A (N.S.), 1983, 38, 7-13) on the characterization of the essential spectrum to the multidimensional non-selfadjoint case. The fourth result shows that, in one dimension, the purely imaginary (non-sparse) step potential has unexpectedly many eigenvalues, comparable to the number of resonances. Our examples show that several known upper bounds are sharp.


title: Enhanced displacement of phase separating liquid mixtures in 2D confined spaces

abstract: Displacing liquid in a confined space is important for technological processes, ranging from porous membrane separation to CO$_{2}$ sequestration. The liquid to be displaced usually consists of multiple components with different solubilities in the displacing liquid. Phase separation and chemical composition gradients in the liquids can influence the displacement rate. In this work, we investigate the effects of liquid composition on the displacement process of ternary liquid mixtures in a quasi-2D microchannel where liquid-liquid phase separation occurs concurrently. We focused on model ternary mixtures containing 1-octanol, ethanol, and water. These mixtures are displaced with water or with ethanol aqueous solution. The spatial distribution of subphases arising from phase separation and the displacement rates of the solution are impacted by the initial ternary solution composition. The boundary between the solution and displacing liquid changes from a defined interface to a diffusive interface as the initial 1-octanol composition in the solution is reduced. The displacement rate also varies non-linearly with the initial 1-octanol composition. The slowest displacement rate arises in the intermediate 1-octanol concentration, where a stable three-zone configuration forms at the boundary. At very low 1-octanol concentration, the displacement rate is fast, associated with droplet formation and motion driven by the chemical concentration gradients formed during phase separation. The excessive energy provided from phase separation may contribute to the enhanced displacement at intermediate to high 1-octanol concentrations, but not at the low 1-octanol concentration with enhancement from induced flow in confinement. The knowledge gained from this study highlights the importance of manipulating phase separation to enhance mass transport in confinement for a wide range of separation processes.


title: Control of Scanning Quantum Dot Microscopy

abstract: Scanning quantum dot microscopy is a recently developed high-resolution microscopy technique that is based on atomic force microscopy and is capable of imaging the electrostatic potential of nanostructures like molecules or single atoms. Recently, it could be shown that it not only yields qualitatively but also quantitatively cutting edge images even on an atomic level. In this paper we present how control is a key enabling element to this. The developed control approach consists of a two-degree-of-freedom control framework that comprises a feedforward and a feedback part. For the latter we design two tailored feedback controllers. The feedforward part generates a reference for the current scanned line based on the previously scanned one. We discuss in detail various aspects of the presented control approach and its implications for scanning quantum dot microscopy. We evaluate the influence of the feedforward part and compare the two proposed feedback controllers. The proposed control algorithms speed up scanning quantum dot microscopy by more than a magnitude and enable to scan large sample areas.


title: A Wells type exact sequence for non-degenerate unitary solutions of the Yang--Baxter equation

abstract: Cycle sets are known to give non-degenerate unitary solutions of the Yang--Baxter equation and linear cycle sets are enriched versions of these algebraic systems. The paper explores the recently developed cohomology and extension theory for linear cycle sets. We derive a four term exact sequence relating 1-cocycles, second cohomology and certain groups of automorphisms arising from central extensions of linear cycle sets. This is an analogue of a similar exact sequence for group extensions known due to Wells. We also compare the exact sequence for linear cycle sets with that for their underlying abelian groups via the forgetful functor and discuss generalities on dynamical 2-cocycles.


title: The classification of left-invariant para-K\"ahler structures on simply connected four-dimensional Lie groups

abstract: We give a complete classification of left invariant para-K\"ahler structures on four-dimensional simply connected Lie groups up to an automorphism. As an application we discuss some curvatures properties of the canonical connection associated to these structures as flat, Ricci flat and existence of Ricci solitons.


title: Quantum states of the Kapitza pendulum

abstract: The quantum states of the Kapitza pendulum are described within the effective potential obtained by the method of averaging over the fast oscillations. An analytical estimate of the energy spectrum of stabilized states is given using approximate model potential. For the lowest states of an inverted pendulum, the spectrum is repeduced by the energies of a harmonic oscillator with perturbation theory corrections. Tunneling effect contribution to the energies of resonance states in the double-well effective potential is estimated. The results of numerical simulations of vibrational and rotational spectra of the Kapitza pendulum by the semiclassical method and by the Numerov algorithm are compared.


title: Machine Learning for Scientific Discovery

abstract: Machine Learning algorithms are good tools for both classification and prediction purposes. These algorithms can further be used for scientific discoveries from the enormous data being collected in our era. We present ways of discovering and understanding astronomical phenomena by applying machine learning algorithms to data collected with radio telescopes. We discuss the use of supervised machine learning algorithms to predict the free parameters of star formation histories and also better understand the relations between the different input and output parameters. We made use of Deep Learning to capture the non-linearity in the parameters. Our models are able to predict with low error rates and give the advantage of predicting in real time once the model has been trained. The other class of machine learning algorithms viz. unsupervised learning can prove to be very useful in finding patterns in the data. We explore how we use such unsupervised techniques on solar radio data to identify patterns and variations, and also link such findings to theories, which help to better understand the nature of the system being studied. We highlight the challenges faced in terms of data size, availability, features, processing ability and importantly, the interpretability of results. As our ability to capture and store data increases, increased use of machine learning to understand the underlying physics in the information captured seems inevitable.


title: Quasi feedback forms for differential-algebraic systems

abstract: We investigate feedback forms for linear time-invariant systems described by differential-algebraic equations. Feedback forms are representatives of certain equivalence classes. For example state space transformations, invertible transformations from the left, and proportional state feedback constitute an equivalence relation. The representative of such an equivalence class, which we call proportional feedback form for the above example, allows to read off relevant system theoretic properties. Our main contribution is to derive a quasi proportional feedback form. This form is advantageous since it provides some geometric insight and is simple to compute, but still allows to read off the relevant structural properties of the control system. We also derive a quasi proportional and derivative feedback form. Similar advantages hold.


title: A continuation method for computing the multilinear Pagerank

abstract: The multilinear Pagerank model [Gleich, Lim and Yu, 2015] is a tensor-based generalization of the Pagerank model. Its computation requires solving a system of polynomial equations that contains a parameter $\alpha \in [0,1)$. For $\alpha \approx 1$, this computation remains a challenging problem, especially since the solution may be non-unique. Extrapolation strategies that start from smaller values of $\alpha$ and `follow' the solution by slowly increasing this parameter have been suggested; however, there are known cases where these strategies fail, because a globally continuous solution curve cannot be defined as a function of $\alpha$. In this paper, we improve on this idea, by employing a predictor-corrector continuation algorithm based on a more general representation of the solutions as a curve in $\mathbb{R}^{n+1}$. We prove several global properties of this curve that ensure the good behavior of the algorithm, and we show in our numerical experiments that this method is significantly more reliable than the existing alternatives.


title: Distributional robustness in minimax linear quadratic control with Wasserstein distance

abstract: To address the issue of inaccurate distributions in practical stochastic systems, a minimax linear-quadratic control method is proposed using the Wasserstein metric. Our method aims to construct a control policy that is robust against errors in an empirical distribution of underlying uncertainty, by adopting an adversary that selects the worst-case distribution. The opponent receives a Wasserstein penalty proportional to the amount of deviation from the empirical distribution. A closed-form expression of the finite-horizon optimal policy pair is derived using a Riccati equation. The result is then extended to the infinite-horizon average cost setting by identifying conditions under which the Riccati recursion converges to the unique positive semi-definite solution to an algebraic Riccati equation. Our method is shown to possess several salient features including closed-loop stability, and an out-of-sample performance guarantee. We also discuss how to optimize the penalty parameter for enhancing the distributional robustness of our control policy. Last but not least, a theoretical connection to the classical $H_\infty$-method is identified from the perspective of distributional robustness.


title: Critical points of the random cluster model with Newman-Ziff sampling

abstract: We present a method for computing transition points of the random cluster model using a generalization of the Newman-Ziff algorithm, a celebrated technique in numerical percolation, to the random cluster model. The new method is straightforward to implement and works for real cluster weight $q>0$. Furthermore, results for an arbitrary number of values of $q$ can be found at once within a single simulation. Because the algorithm used to sweep through bond configurations is identical to that of Newman and Ziff, which was conceived for percolation, the method loses accuracy for large lattices when $q>1$. However, by sampling the critical polynomial, accurate estimates of critical points in two dimensions can be found using relatively small lattice sizes, which we demonstrate here by computing critical points for non-integer values of $q$ on the square lattice, to compare with the exact solution, and on the unsolved non-planar square matching lattice. The latter results would be much more difficult to obtain using other techniques.


title: Cloud Broker: A Systematic Mapping Study

abstract: In a cloud environment, a cloud broker is an important entity that works as an independent middleware between cloud customers and providers to address issues and conduct negotiations related to satisfying both customer preferences and service provider profits. In recent years, researchers have published many articles which directly or indirectly address this research area. A systematic method is vital for extracting all search spaces (journals, conferences, and workshops) and primary studies (articles) conducted in the cloud broker field and then selecting some of the highest quality studies. The proposed systematic review includes a comprehensive three-tier search strategy (manual search, backward snowballing, and database search). The detailed explanation of the reviewing process is inserted in Appendix A. In the search methodology, qualitative criteria have been defined to select studies with the highest quality and the most relevance among all search spaces. In the present study, out of 1,928 extracted search spaces, 171 search spaces have been selected based on the defined quality criteria. Then, 1,298 articles have been extracted from these 171 selected search spaces. As a result, 496 high-quality papers have been selected among the mentioned papers. The chosen papers were published in prestigious journals, conferences, and workshops from 2009 through 2019. In the current Systematic Mapping Study (SMS), eight research questions have been designed for the purpose of identifying information that is significant to the cloud broker field, such as the most critical and debated topics, existing trends and issues, active researchers and countries, commonly used techniques in building cloud brokers, evaluation methods, the amount of research conducted by year and the place of publication, and the most important active search spaces.


title: A Simulation-based End-to-End Learning Framework for Evidential Occupancy Grid Mapping

abstract: Evidential occupancy grid maps (OGMs) are a popular representation of the environment of automated vehicles. Inverse sensor models (ISMs) are used to compute OGMs from sensor data such as lidar point clouds. Geometric ISMs show a limited performance when estimating states in unobserved but inferable areas and have difficulties dealing with ambiguous input. Deep learning-based ISMs face the challenge of limited training data and they often cannot handle uncertainty quantification yet. We propose a deep learning-based framework for learning an OGM algorithm which is both capable of quantifying first- and second-order uncertainty and which does not rely on manually labeled data. Results on synthetic and on real-world data show superiority over other approaches. Source code and datasets are available at https://github.com/ika-rwth-aachen/EviLOG


title: Defining Preferred and Natural Robot Motions in Immersive Telepresence from a First-Person Perspective

abstract: This paper presents some early work and future plans regarding how the autonomous motions of a telepresence robot affect a person embodied in the robot through a head-mounted display. We consider the preferences, comfort, and the perceived naturalness of aspects of piecewise linear paths compared to the same aspects on a smooth path. In a user study, thirty-six subjects (eighteen females) watched panoramic videos of three different paths through a simulated museum in virtual reality and responded to questionnaires regarding each path. We found that comfort had a strong effect on path preference, and that the subjective feeling of naturalness also had a strong effect on path preference, even though people consider different things as natural. We describe a categorization of the responses regarding the naturalness of the robot's motion and provide a recommendation on how this can be applied more broadly. Although immersive robotic telepresence is increasingly being used for remote education, clinical care, and to assist people with disabilities or mobility complications, the full potential of this technology is limited by issues related to user experience. Our work addresses these shortcomings and will enable the future personalization of telepresence experiences for the improvement of overall remote communication and the enhancement of the feeling of presence in a remote location.


title: Compact freeform illumination design by deblurring the extended sources

abstract: Illumination is the deliberate utilization of light to realize practical or aesthetic effects. The designers combine with the environmental considerations, energy-saving goals, and technology advances with fundamental physics to develop lighting solutions to satisfy all of our ever-changing needs. Achieving highly efficient and precise control of the energy output of light sources while maintaining compact optical structures is the ultimate goal of illumination design. To realize miniaturized and lightweight luminaires, the design process must consider the extents of light sources. However, the illumination design for extended sources is still a challenging and unsolved problem. Here, we propose a method to design ultra-performance illumination optics enabled by freeform optical surfaces. The proposed method is very general with no limitations of far-field approximation and Lambertian luminescent property. We demonstrate the feasibility and efficiency of the proposed method by designing several freeform lenses realizing accurate and highly efficient illumination control as well as ultra-compact structures.


title: The Distance and Peculiar Velocity of the Norma cluster (ACO 3627) using the Near-Infrared $J$ and $K_s$-band Fundamental Plane Relations

abstract: We report distance measurements for the Norma cluster based on the near-infrared $J$- and $K_s$-band Fundamental Plane (FP) relations. Our simultaneous $J$ and $K_s$-band photometry analyses were performed using 31 early-type galaxies in the nearby Norma cluster obtained using the 1.4 m InfraRed Survey Facility (IRSF) at the South African Astronomical Observatory. Our final $K_s$-band FP sample consists of 41 early-type galaxies from the Norma cluster observed using the IRSF and the New Technology Telescope (NTT) at the European Southern Observatory. This is the largest cluster sample used for peculiar velocity studies in the Great Attractor region to date. From the $K_s$-band FP, we find a distance to the Norma cluster of $4915 \pm 121$ km s$^{-1}$. The implied peculiar velocity for Norma is $44 \pm 151$ km s$^{-1}$ which further supports a small peculiar velocity for the Norma cluster.


title: Combinatorial Bandits under Strategic Manipulations

abstract: Strategic behavior against sequential learning methods, such as "click framing" in real recommendation systems, have been widely observed. Motivated by such behavior we study the problem of combinatorial multi-armed bandits (CMAB) under strategic manipulations of rewards, where each arm can modify the emitted reward signals for its own interest. This characterization of the adversarial behavior is a relaxation of previously well-studied settings such as adversarial attacks and adversarial corruption. We propose a strategic variant of the combinatorial UCB algorithm, which has a regret of at most $O(m\log T + m B_{max})$ under strategic manipulations, where $T$ is the time horizon, $m$ is the number of arms, and $B_{max}$ is the maximum budget of an arm. We provide lower bounds on the budget for arms to incur certain regret of the bandit algorithm. Extensive experiments on online worker selection for crowdsourcing systems, online influence maximization and online recommendations with both synthetic and real datasets corroborate our theoretical findings on robustness and regret bounds, in a variety of regimes of manipulation budgets.


title: On Interpretability and Similarity in Concept-Based Machine Learning

abstract: Machine Learning (ML) provides important techniques for classification and predictions. Most of these are black-box models for users and do not provide decision-makers with an explanation. For the sake of transparency or more validity of decisions, the need to develop explainable/interpretable ML-methods is gaining more and more importance. Certain questions need to be addressed: How does an ML procedure derive the class for a particular entity? Why does a particular clustering emerge from a particular unsupervised ML procedure? What can we do if the number of attributes is very large? What are the possible reasons for the mistakes for concrete cases and models? For binary attributes, Formal Concept Analysis (FCA) offers techniques in terms of intents of formal concepts, and thus provides plausible reasons for model prediction. However, from the interpretable machine learning viewpoint, we still need to provide decision-makers with the importance of individual attributes to the classification of a particular object, which may facilitate explanations by experts in various domains with high-cost errors like medicine or finance. We discuss how notions from cooperative game theory can be used to assess the contribution of individual attributes in classification and clustering processes in concept-based machine learning. To address the 3rd question, we present some ideas on how to reduce the number of attributes using similarities in large contexts.


title: Joint value distribution of $L$-functions on the critical line

abstract: In this paper, we discuss the joint value distribution of $L$-functions in a suitable class. We obtain joint large deviations results in the central limit theorem for these $L$-functions and some mean value theorems, which give evidence that different $L$-functions are "statistically independent".


title: The hollow Gaussian beam propagation on curved surface based on matrix optics method

abstract: In this paper, ABCD matrix is introduced to study the paraxial transmission of light on a constant gaussian curvature surface (CGCS), which is the first time to our knowledge. It is also proved that the curved surface can be used as the implementation of fractional Fourier transform, we further generalize that we can obtain the transfer matrix of an arbitrary surface with gently varying curvature by matrix optics. As a beam propagation example, based on the Collins integral, an analytical propagation formula for the hollow Gaussian beams (HGBs) on the CGCS is derived. The propagation characteristics of HGBs on one CGCS are illustrated graphically in detail, including the change of dark spots size and splitting rays. Besides its propagation periodicity and diffraction properties, a criterion for convergence and divergence of the spot size is proposed. The area of the dark region of the HGBs can easily be controlled by proper choice of the beam parameters and the shape of CGCS. In addition, we also study the special propagation properties of the hollow beam with a fractional order. Compared with HGBs in flat space, these novel characteristics of HGBs propagation on curved surface may further expand the application range of hollow beam.


title: Design and Control of a Highly Redundant Rigid-Flexible Coupling Robot to Assist the COVID-19 Oropharyngeal-Swab Sampling

abstract: The outbreak of novel coronavirus pneumonia (COVID-19) has caused mortality and morbidity worldwide. Oropharyngeal-swab (OP-swab) sampling is widely used for the diagnosis of COVID-19 in the world. To avoid the clinical staff from being affected by the virus, we developed a 9-degree-of-freedom (DOF) rigid-flexible coupling (RFC) robot to assist the COVID-19 OP-swab sampling. This robot is composed of a visual system, UR5 robot arm, micro-pneumatic actuator and force-sensing system. The robot is expected to reduce risk and free up the clinical staff from the long-term repetitive sampling work. Compared with a rigid sampling robot, the developed force-sensing RFC robot can facilitate OP-swab sampling procedures in a safer and softer way. In addition, a varying-parameter zeroing neural network-based optimization method is also proposed for motion planning of the 9-DOF redundant manipulator. The developed robot system is validated by OP-swab sampling on both oral cavity phantoms and volunteers.


title: What's in a GitHub Repository? -- A Software Documentation Perspective

abstract: Developers use and contribute to repositories on GitHub. Documentation present in the repositories serves as an important source by helping developers to understand, maintain and contribute to the project. Currently, documentation in a repository is diversified, among various files, with most of it present in ReadMe files. However, other software artifacts in the repository, such as issue reports and pull requests could also contribute to documentation, without documentation being explicitly specified. Hence, in this paper, we propose a taxonomy of documentation sources by analyzing different software artifacts, developer interviews and card-sorting approach. We inspected multiple artifacts of 950 public GitHub repositories, written in four different programming languages, C++, C#, Python and Java, and analyzed the type and amount of documentation that could be extracted from these artifacts. To this end, we observe that, about 25.93% of information extracted from all sources proposed in the taxonomy contains error-related documentation, and that pull requests contribute to around 18.21% of extracted information.


title: Scene Retrieval for Contextual Visual Mapping

abstract: Visual navigation localizes a query place image against a reference database of place images, also known as a `visual map'. Localization accuracy requirements for specific areas of the visual map, `scene classes', vary according to the context of the environment and task. State-of-the-art visual mapping is unable to reflect these requirements by explicitly targetting scene classes for inclusion in the map. Four different scene classes, including pedestrian crossings and stations, are identified in each of the Nordland and St. Lucia datasets. Instead of re-training separate scene classifiers which struggle with these overlapping scene classes we make our first contribution: defining the problem of `scene retrieval'. Scene retrieval extends image retrieval to classification of scenes defined at test time by associating a single query image to reference images of scene classes. Our second contribution is a triplet-trained convolutional neural network (CNN) to address this problem which increases scene classification accuracy by up to 7% against state-of-the-art networks pre-trained for scene recognition. The second contribution is an algorithm `DMC' that combines our scene classification with distance and memorability for visual mapping. Our analysis shows that DMC includes 64% more images of our chosen scene classes in a visual map than just using distance interval mapping. State-of-the-art visual place descriptors AMOS-Net, Hybrid-Net and NetVLAD are finally used to show that DMC improves scene class localization accuracy by a mean of 3% and localization accuracy of the remaining map images by a mean of 10% across both datasets.


title: Circularly polarized states and propagating bound states in the continuum in a periodic array of cylinders

abstract: Bound states in the continuum (BICs) in a periodic structure sandwiched between two homogeneous media have interesting properties and useful applications in photonics. The topological nature of BICs was previously revealed based on a topological charge related to the far-field polarization vector of the surrounding resonant states. Recently, it was established that when a symmetry-protected BIC (with a nonzero topological charge) is destroyed by a generic symmetry-breaking perturbation, a pair of circularly polarized resonant states (CPSs) emerge and the net topological charge is conserved. A periodic structure can also support propagating BICs with a nonzero wavevector. These BICs are not protected by symmetry in the sense of symmetry mismatch, but they need symmetry for their robust existence. Based on a highly accurate computational method for a periodic array of slightly noncircular cylinders, we show that a propagating BIC is typically destroyed by a structural perturbation that breaks only the in-plane inversion symmetry, and when this happens, a pair of CPSs of opposite handedness emerge so that the net topological charge is conserved. We also study the generation and annihilation of CPSs when a structural parameter is varied. It is shown that two CPSs with opposite topological charge and same handedness, connected to two BICs or in a continuous branch from one BIC, may collapse and become a CPS with a zero charge. Our study clarifies the important connection between symmetry and topological charge conservation.


title: Swivel: Hardening WebAssembly against Spectre

abstract: We describe Swivel, a new compiler framework for hardening WebAssembly (Wasm) against Spectre attacks. Outside the browser, Wasm has become a popular lightweight, in-process sandbox and is, for example, used in production to isolate different clients on edge clouds and function-as-a-service platforms. Unfortunately, Spectre attacks can bypass Wasm's isolation guarantees. Swivel hardens Wasm against this class of attacks by ensuring that potentially malicious code can neither use Spectre attacks to break out of the Wasm sandbox nor coerce victim code-another Wasm client or the embedding process-to leak secret data. We describe two Swivel designs, a software-only approach that can be used on existing CPUs, and a hardware-assisted approach that uses extension available in Intel 11th generation CPUs. For both, we evaluate a randomized approach that mitigates Spectre and a deterministic approach that eliminates Spectre altogether. Our randomized implementations impose under 10.3% overhead on the Wasm-compatible subset of SPEC 2006, while our deterministic implementations impose overheads between 3.3% and 240.2%. Though high on some benchmarks, Swivel's overhead is still between 9x and 36.3x smaller than existing defenses that rely on pipeline fences.


title: Improving Approximate Optimal Transport Distances using Quantization

abstract: Optimal transport (OT) is a popular tool in machine learning to compare probability measures geometrically, but it comes with substantial computational burden. Linear programming algorithms for computing OT distances scale cubically in the size of the input, making OT impractical in the large-sample regime. We introduce a practical algorithm, which relies on a quantization step, to estimate OT distances between measures given cheap sample access. We also provide a variant of our algorithm to improve the performance of approximate solvers, focusing on those for entropy-regularized transport. We give theoretical guarantees on the benefits of this quantization step and display experiments showing that it behaves well in practice, providing a practical approximation algorithm that can be used as a drop-in replacement for existing OT estimators.


title: Energy Decay of some boundary coupled systems involving wave$\backslash$ Euler-Bernoulli beam with one locally singular fractional Kelvin-Voigt damping

abstract: In this paper, we investigate the energy decay of hyperbolic systems of wave-wave, wave-Euler- Bernoulli beam and beam-beam types. The two equations are coupled through boundary connection with only one localized non-smooth fractional Kelvin-Voigt damping. First, we reformulate each system into an augmented model and using a general criteria of Arendt-Batty, we prove that our models are strongly stable. Next, by using frequency domain approach, combined with multiplier technique and some interpolation inequalities, we establish different types of polynomial energy decay rate which depends on the order of the fractional derivative and the type of the damped equation in the system.


title: Distributed Online Learning with Multiple Kernels

abstract: We consider the problem of learning a nonlinear function over a network of learners in a fully decentralized fashion. Online learning is additionally assumed, where every learner receives continuous streaming data locally. This learning model is called a fully distributed online learning (or a fully decentralized online federated learning). For this model, we propose a novel learning framework with multiple kernels, which is named DOMKL. The proposed DOMKL is devised by harnessing the principles of an online alternating direction method of multipliers and a distributed Hedge algorithm. We theoretically prove that DOMKL over T time slots can achieve an optimal sublinear regret, implying that every learner in the network can learn a common function which has a diminishing gap from the best function in hindsight. Our analysis also reveals that DOMKL yields the same asymptotic performance of the state-of-the-art centralized approach while keeping local data at edge learners. Via numerical tests with real datasets, we demonstrate the effectiveness of the proposed DOMKL on various online regression and time-series prediction tasks.


title: Synthesis of Hybrid Automata with Affine Dynamics from Time-Series Data

abstract: Formal design of embedded and cyber-physical systems relies on mathematical modeling. In this paper, we consider the model class of hybrid automata whose dynamics are defined by affine differential equations. Given a set of time-series data, we present an algorithmic approach to synthesize a hybrid automaton exhibiting behavior that is close to the data, up to a specified precision, and changes in synchrony with the data. A fundamental problem in our synthesis algorithm is to check membership of a time series in a hybrid automaton. Our solution integrates reachability and optimization techniques for affine dynamical systems to obtain both a sufficient and a necessary condition for membership, combined in a refinement framework. The algorithm processes one time series at a time and hence can be interrupted, provide an intermediate result, and be resumed. We report experimental results demonstrating the applicability of our synthesis approach.


title: Random Forest based Qantile Oriented Sensitivity Analysis indices estimation

abstract: We propose a random forest based estimation procedure for Quantile Oriented Sensitivity Analysis-QOSA. In order to be efficient, a cross validation step on the leaf size of trees is required. Our full estimation procedure is tested on both simulated data and a real dataset.


title: Time-Series Imputation with Wasserstein Interpolation for Optimal Look-Ahead-Bias and Variance Tradeoff

abstract: Missing time-series data is a prevalent practical problem. Imputation methods in time-series data often are applied to the full panel data with the purpose of training a model for a downstream out-of-sample task. For example, in finance, imputation of missing returns may be applied prior to training a portfolio optimization model. Unfortunately, this practice may result in a look-ahead-bias in the future performance on the downstream task. There is an inherent trade-off between the look-ahead-bias of using the full data set for imputation and the larger variance in the imputation from using only the training data. By connecting layers of information revealed in time, we propose a Bayesian posterior consensus distribution which optimally controls the variance and look-ahead-bias trade-off in the imputation. We demonstrate the benefit of our methodology both in synthetic and real financial data.


title: IPFS and Friends: A Qualitative Comparison of Next Generation Peer-to-Peer Data Networks

abstract: Decentralized, distributed storage offers a way to reduce the impact of data silos as often fostered by centralized cloud storage. While the intentions of this trend are not new, the topic gained traction due to technological advancements, most notably blockchain networks. As a consequence, we observe that a new generation of peer-to-peer data networks emerges. In this survey paper, we therefore provide a technical overview of the next generation data networks. We use select data networks to introduce general concepts and to emphasize new developments. Specifically, we provide a deeper outline of the Interplanetary File System and a general overview of Swarm, the Hypercore Protocol, SAFE, Storj, and Arweave. We identify common building blocks and provide a qualitative comparison. From the overview, we derive future challenges and research goals concerning data networks.


title: Models of random subtrees of a graph

abstract: Consider a connected graph $G=(E,V)$ with $N=|V|$ vertices. The main purpose of this paper is to explore the question of uniform sampling of a subtree of $G$ with $n$ nodes, for some $n\leq N$ (the spanning tree case correspond to $n=N$, and is already deeply studied in the literature). We provide new asymptotically exact simulation methods using Markov chains for general connected graphs $G$, and any $n\leq N$. We highlight the case of the uniform subtree of $\mathbb{Z}^2$ with $n$ nodes, containing the origin $(0,0)$ for which Schramm asked several questions. We produce pictures, statistics, and some conjectures. A second aim of the paper is devoted to surveying other models of random subtrees of a graph, among them, DLA models, the first passage percolation, the uniform spanning tree and the minimum spanning tree. We also provide new models, some statistics, and some conjectures.


title: Implementation of 3D degridding algorithm on the NVIDIA GPUs using CUDA

abstract: Practical aperture synthesis imaging algorithms work by iterating between estimating the sky brightness distribution and a comparison of a prediction based on this estimate with the measured data ("visibilities"). Accuracy in the latter step is crucial but is made difficult by irregular and non-planar sampling of data by the telescope. In this work we present a GPU implementation of 3d de-gridding which accurately deals with these two difficulties and is designed for distributed operation. We address the load balancing issues caused by large variation in visibilities that need to be computed. Using CUDA and NVidia GPUs we measure performance up to 1.2 billion visibilities per second.


title: Performance Comparison for Scientific Computations on the Edge via Relative Performance

abstract: In a typical Internet-of-Things setting that involves scientific applications, a target computation can be evaluated in many different ways depending on the split of computations among various devices. On the one hand, different implementations (or algorithms)--equivalent from a mathematical perspective--might exhibit significant difference in terms of performance. On the other hand, some of the implementations are likely to show similar performance characteristics. In this paper, we focus on analyzing the performance of a given set of algorithms by clustering them into performance classes. To this end, we use a measurement-based approach to evaluate and score algorithms based on pair-wise comparisons; we refer to this approach as"Relative performance analysis". Each comparison yields one of three outcomes: one algorithm can be "better", "worse", or "equivalent" to another; those algorithms evaluating to have equivalent performance are merged into the same performance class. We show that our clustering methodology facilitates algorithm selection with respect to more than one metric; for instance, from the subset of equivalently fast algorithms, one could then select an algorithm that consumes the least energy on a certain device.


title: Spiraling of sub-Riemannian geodesics around the Reeb flow in the 3D contact case

abstract: We consider a closed three-dimensional contact sub-Riemannian manifold. The objective of this note is to provide a precise description of the sub-Riemannian geodesics with large initial momenta: we prove that they "spiral around the Reeb orbits", not only in the phase space but also in the configuration space. Our analysis is based on a normal form along any Reeb orbit due to Melrose.


title: The derived sequence of a pre-Jaffard family

abstract: We introduce the concept of \emph{pre-Jaffard family}, a generalization of Jaffard families obtained by substituting the locally finite hypothesis with a much weaker compactness hypothesis. From any such family, we construct a sequence of overrings of the starting domain that allows to decompose stable semistar operations and singular length functions in more cases than what is allowed by Jaffard families. We also apply the concept to one-dimensional domains, unifying the treatment of sharp and dull degree of a Pr\"ufer domain.


title: Bragg-Curve Simulation of Carbon-Ion Beams for Particle-therapy Applications: a study with the GEANT4 toolkit

abstract: We used the GEANT4 Monte Carlo MC Toolkit to simulate carbon ion beams incident on water, tissue, and bone, taking into account nuclear fragmentation reactions. Upon increasing the energy of the primary beam, the position of the Bragg-Peak transfers to a location deeper inside the phantom. For different materials, the peak is located at a shallower depth along the beam direction and becomes sharper with increasing electron density NZ. Subsequently, the generated depth dose of the Bragg curve is then benchmarked with experimental data from GSI in Germany. The results exhibit a reasonable correlation with GSI experimental data with an accuracy of between 0.02 and 0.08 cm, thus establishing the basis to adopt MC in heavy-ion treatment planning. The Kolmogorov-Smirnov K-S test further ascertained from a statistical point of view that the simulation data matched the experimentally measured data very well. The two-dimensional isodose contours at the entrance were compared to those around the peak position and in the tail region beyond the peak, showing that bone produces more dose, in comparison to both water and tissue, due to secondary doses. In the water, the results show that the maximum energy deposited per fragment is mainly attributed to secondary carbon ions, followed by secondary boron and beryllium. Furthermore, the number of protons produced is the highest, thus making the maximum contribution to the total dose deposition in the tail region. Finally, the associated spectra of neutrons and photons were analyzed. The mean neutron energy value was found to be 16.29 MeV, and 1.03 MeV for the secondary gamma. However, the neutron dose was found to be negligible as compared to the total dose due to their longer range.


title: On the viscosity approach to a class of fully nonlinear elliptic equations

abstract: In this paper, we study some properties of viscosity sub/super-solutions of a class of fully nonlinear elliptic equations relative to the eigenvalues of the complex Hessian. We show that every viscosity subsolution is approximated by a decreasing sequence of smooth subsolutions. When the equations satisfy some conditions on the limit at infinity, we verify that the comparison principle holds, and as a sequence, we obtain a result about the existence of solution of the Dirichlet problem. Using the comparison principle, we show that, under suitable conditions, a Perron-Bremermann envelope can be approximated by a decreasing sequence of viscosity solutions.


title: Quantum invariants of knotoids

abstract: In this paper, we construct quantum invariants for knotoid diagrams in $\mathbb{R}^2$. The diagrams are arranged with respect to a given direction in the plane ({\it Morse knotoids}). A Morse knotoid diagram can be decomposed into basic elementary diagrams each of which is associated to a matrix that yields solutions of the quantum Yang-Baxter equation. We recover the bracket polynomial, and define the rotational bracket polynomial, the binary bracket polynomial, the Alexander polynomial, the generalized Alexander polynomial and an infinity of specializations of the Homflypt polynomial for Morse knotoids via quantum state sum models.


title: Blockchained Federated Learning for Threat Defense

abstract: Given the increasing complexity of threats in smart cities, the changing environment, and the weakness of traditional security systems, which in most cases fail to detect serious threats such as zero-day attacks, the need for alternative more active and more effective security methods keeps increasing. Such approaches are the adoption of intelligent solutions to prevent, detect and deal with threats or anomalies under the conditions and the operating parameters of the infrastructure in question. This research paper introduces the development of an intelligent Threat Defense system, employing Blockchain Federated Learning, which seeks to fully upgrade the way passive intelligent systems operate, aiming at implementing an Advanced Adaptive Cooperative Learning (AACL) mechanism for smart cities networks. The AACL is based on the most advanced methods of computational intelligence while ensuring privacy and anonymity for participants and stakeholders. The proposed framework combines Federated Learning for the distributed and continuously validated learning of the tracing algorithms. Learning is achieved through encrypted smart contracts within the blockchain technology, for unambiguous validation and control of the process. The aim of the proposed Framework is to intelligently classify smart cities networks traffic derived from Industrial IoT (IIoT) by Deep Content Inspection (DCI) methods, in order to identify anomalies that are usually due to Advanced Persistent Threat (APT) attacks.


title: Self-consistent adiabatic inspiral and transition motion

abstract: The transition motion of a point particle around the last stable orbit of Kerr is described at leading order in the transition-timescale expansion. Taking systematically into account all self-force effects, we prove that the transition motion is still described by the Painlev\'e transcendent equation of the first kind. Using an asymptotically matched expansions scheme, we consistently match the quasi-circular adiabatic inspiral with the transition motion. The matching requires us to take into account the secular change of angular velocity due to radiation reaction during the adiabatic inspiral.


title: Active Modular Environment for Robot Navigation

abstract: This paper presents a novel robot-environment interaction in navigation tasks such that robots have neither a representation of their working space nor planning function, instead, an active environment takes charge of these aspects. This is realized by spatially deploying computing units, called cells, and making cells manage traffic in their respective physical region. Different from stigmegic approaches, cells interact with each other to manage environmental information and to construct instructions on how robots move. As a proof-of-concept, we present an architecture called AFADA and its prototype, consisting of modular cells and robots moving on the cells. The instructions from cells are based on a distributed routing algorithm and a reservation protocol. We demonstrate that AFADA achieves efficient robot moves for single-robot navigation in a dynamic environment changing its topology with a stochastic model, comparing to self-navigation by a robot itself. This is followed by several demos, including multi-robot navigation, highlighting the power of offloading both representation and planning from robots to the environment. We expect that the concept of AFADA contributes to developing the infrastructure for multiple robots because it can engage online and lifelong planning and execution.


title: Investigation of spatiotemporal output beam profile instabilities from differentially pumped capillaries

abstract: Differentially pumped capillaries, i.e., capillaries operated in a pressure gradient environment, are widely used for nonlinear pulse compression. In this work, we show that strong pressure gradients and high gas throughputs can cause spatiotemporal instabilities of the output beam profile. The instabilities occur with a sudden onset as the flow evolves from laminar to turbulent. Based on the experimental and numerical results, we derive guidelines to predict the onset of those instabilities and discuss possible applications in the context of nonlinear flow dynamics.


title: Observable quality assessment of broadband very long baseline interferometry system

abstract: The next-generation, broadband geodetic very long baseline interferometry system, named VGOS, is developing its global network, and VGOS networks with a small size of 3--7 stations have already made broadband observations from 2017 to 2019. We made quality assessments for two kinds of observables in the 21 VGOS sessions currently available: group delay and differential total electron content ($\delta$TEC). Our study reveals that the random measurement noise of VGOS group delays is at the level of less than 2 ps (1 ps = 10$^{-12}$ s), while the contributions from systematic error sources, mainly source structure related, are at the level of 20 ps. Due to the significant improvement in measurement noise, source structure effects with relatively small magnitudes that are not overwhelming in the S/X VLBI system, for instance 10 ps, are clearly visible in VGOS observations. Another critical error source in VGOS observations is discrete delay jumps, for instance, a systematic offset of about 310 ps or integer multiples of that. The predominant causative factor is found to be related to source structure. The measurement noise level of $\delta$TEC observables is about 0.07 TECU, but the systematic effects are five times larger than that. A strong correlation between group delay and $\delta$TEC observables is discovered with a trend of 40 ps/TECU for observations with large structure effects; there is a second trend in the range 60 ps/TECU to 70 ps/TECU when the measurement noise is dominant.


title: Migration of CMSWEB Cluster at CERN to Kubernetes

abstract: The CMS experiment heavily relies on the CMSWEB cluster to host critical services for its operational needs. The cluster is deployed on virtual machines (VMs) from the CERN OpenStack cloud and is manually maintained by operators and developers. The release cycle is composed of several steps, from building RPMs, their deployment to perform validation, and integration tests. To enhance the sustainability of the CMSWEB cluster, CMS decided to migrate its cluster to a containerized solution such as Docker, orchestrated with Kubernetes (k8s). This allows us to significantly reduce the release upgrade cycle, follow the end-to-end deployment procedure, and reduce operational cost. This paper gives an overview of the current CMSWEB cluster and its issues. We describe the new architecture of the CMSWEB cluster in Kubernetes. We also provide a comparison of VM and Kubernetes deployment approaches and report on lessons learned during the migration process.


title: Adaptive Robust Large Volatility Matrix Estimation Based on High-Frequency Financial Data

abstract: Several novel statistical methods have been developed to estimate large integrated volatility matrices based on high-frequency financial data. To investigate their asymptotic behaviors, they require a sub-Gaussian or finite high-order moment assumption for observed log-returns, which cannot account for the heavy-tail phenomenon of stock-returns. Recently, a robust estimator was developed to handle heavy-tailed distributions with some bounded fourth-moment assumption. However, we often observe that log-returns have heavier tail distribution than the finite fourth-moment and that the degrees of heaviness of tails are heterogeneous across asset and over time. In this paper, to deal with the heterogeneous heavy-tailed distributions, we develop an adaptive robust integrated volatility estimator that employs pre-averaging and truncation schemes based on jump-diffusion processes. We call this an adaptive robust pre-averaging realized volatility (ARP) estimator. We show that the ARP estimator has a sub-Weibull tail concentration with only finite 2$\alpha$-th moments for any $\alpha>1$. In addition, we establish matching upper and lower bounds to show that the ARP estimation procedure is optimal. To estimate large integrated volatility matrices using the approximate factor model, the ARP estimator is further regularized using the principal orthogonal complement thresholding (POET) method. The numerical study is conducted to check the finite sample performance of the ARP estimator.


title: A crank for bipartitions with designated summands

abstract: Andrews, Lewis and Lovejoy introduced the partition function $PD(n)$ as the number of partitions of $n$ with designated summands. A bipartition of $n$ is an ordered pair of partitions $(\pi_1, \pi_2)$ with the sum of all of the parts being $n$. In this paper, we introduce a generalized crank named the $pd$-crank for bipartitions with designated summands and give some inequalities for the $pd$-crank of bipartitions with designated summands modulo 2 and 3. We also define the $pd$-crank moments weighted by the parity of $pd$-cranks $\mu_{2k,bd}(-1,n)$ and show the positivity of $(-1)^n\mu_{2k,bd}(-1,n)$. Let $M_{bd}(m,n)$ denote the number of bipartitions of $n$ with designated summands with $pd$-crank $m$. We prove a monotonicity property of $pd$-cranks of bipartitions with designated summands and find that the sequence $\{M_{bd}(m,n)\}_{|m|\leq n}$ is unimodal for $n\not= 1,5,7$.