title: Online Bin Packing with Known T

abstract: In the online bin packing problem, a sequence of items is revealed one at a time, and each item must be packed into an available bin instantly upon its arrival. In this paper, we revisit the problem under a setting where the total number of items T is known in advance, also known as the closed online bin packing problem. Specifically, we study both the stochastic model and the random permutation model. We develop and analyze an adaptive algorithm that solves an offline bin packing problem at geometric time intervals and uses the offline optimal solution to guide online packing decisions. Under both models, we show that the algorithm achieves C\sqrt{T} regret (in terms of the number of used bins) compared to the hindsight optimal solution, where C is a universal constant (<= 13) that bears no dependence on the underlying distribution or the item sizes. The result shows the lower bound barrier of \Omega(\sqrt{T \log T}) in Shor (1986) can be broken with solely the knowledge of the horizon T, but without a need of knowing the underlying distribution. As to the algorithm analysis, we develop a new approach to analyzing the packing dynamic using the notion of exchangeable random variables. The approach creates a symmetrization between the offline solution and the online solution, and it is used to analyze both the algorithm performance and various benchmarks related to the bin packing problem. For the latter one, our analysis provides an alternative (probably simpler) treatment and tightens the analysis of the asymptotic benchmark of the stochastic bin packing problem in Rhee and Talagrand (1989a,b). As the analysis only relies on a symmetry between the offline and online problems, the algorithm and benchmark analyses can be naturally extended from the stochastic model to the random permutation model.


title: Maxwell-scalar device based on the electric dipole

abstract: In this work we study the electric field of a dipole immersed in a medium with permittivity controlled by a real scalar field which is non-minimally coupled to the Maxwell field. We model the system with an interesting function, which allows the presence of exact solutions, describing the possibility of the permittivity to encapsulate the charges at very high values, giving rise to an effect that is not present in the standard situation. The results are of direct interest to applications for emission and absorption of radiation, and may motivate new studies concerning binary stars and black holes in gravity scenarios of current interest.


title: Collective variable discovery in the age of machine learning: reality, hype and everything in between

abstract: Understanding kinetics and thermodynamics profile of biomolecules is necessary to understand their functional roles which has a major impact in mechanism driven drug discovery. Molecular dynamics simulation has been routinely used to understand conformational dynamics and molecular recognition in biomolecules. Statistical analysis of high-dimensional spatiotemporal data generated from molecular dynamics simulation requires identification of few low-dimensional variables which can describe essential dynamics of a system without significant loss of informations. In physical chemistry, these low-dimensional variables often called collective variables. Collective variables are used to generated reduced representation of free energy surface and calculate transition probabilities between different metastable basins. However the choice of collective variables is not trivial for complex systems. Collective variables ranges from geometric criteria's such as distances, dihedral angles to abstract ones such as weighted linear combinations of multiple geometric variables. Advent of machine learning algorithms led to increasing use of abstract collective variables to represent biomolecular dynamics. In this review, I will highlight several nuances of commonly used collective variables ranging from geometric to abstract ones. Further, I will put forward some cases where machine learning based collective variables were used to describe simple systems which in principle could have been described by geometric ones. Finally, I will put forward my thoughts on artificial general intelligence and how it can be used to discover and predict collective variables from spatiotemporal data generated by molecular dynamics simulations.


title: A New Sentence Extraction Strategy for Unsupervised Extractive Summarization Methods

abstract: In recent years, text summarization methods have attracted much attention again thanks to the researches on neural network models. Most of the current text summarization methods based on neural network models are supervised methods which need large-scale datasets. However, large-scale datasets are difficult to obtain in practical applications. In this paper, we model the task of extractive text summarization methods from the perspective of Information Theory, and then describe the unsupervised extractive methods with a uniform framework. To improve the feature distribution and to decrease the mutual information of summarization sentences, we propose a new sentence extraction strategy which can be applied to existing unsupervised extractive methods. Experiments are carried out on different datasets, and results show that our strategy is indeed effective and in line with expectations.


title: Quantifying Adaptability in Pre-trained Language Models with 500 Tasks

abstract: When a neural language model (LM) is adapted to perform a new task, what aspects of the task predict the eventual performance of the model? In NLP, systematic features of LM generalization to individual examples are well characterized, but systematic aspects of LM adaptability to new tasks are not nearly as well understood. We present a large-scale empirical study of the features and limits of LM adaptability using a new benchmark, TaskBench500, built from 500 procedurally generated sequence modeling tasks. These tasks combine core aspects of language processing, including lexical semantics, sequence processing, memorization, logical reasoning, and world knowledge. Using TaskBench500, we evaluate three facets of adaptability, finding that: (1) adaptation procedures differ dramatically in their ability to memorize small datasets; (2) within a subset of task types, adaptation procedures exhibit compositional adaptability to complex tasks; and (3) failure to match training label distributions is explained by mismatches in the intrinsic difficulty of predicting individual labels. Our experiments show that adaptability to new tasks, like generalization to new examples, can be systematically described and understood, and we conclude with a discussion of additional aspects of adaptability that could be studied using the new benchmark.


title: Simultaneously Predicting Multiple Plant Traits from Multiple Sensors via Deformable CNN Regression

abstract: Trait measurement is critical for the plant breeding and agricultural production pipeline. Typically, a suite of plant traits is measured using laborious manual measurements and then used to train and/or validate higher throughput trait estimation techniques. Here, we introduce a relatively simple convolutional neural network (CNN) model that accepts multiple sensor inputs and predicts multiple continuous trait outputs - i.e. a multi-input, multi-output CNN (MIMO-CNN). Further, we introduce deformable convolutional layers into this network architecture (MIMO-DCNN) to enable the model to adaptively adjust its receptive field, model complex variable geometric transformations in the data, and fine-tune the continuous trait outputs. We examine how the MIMO-CNN and MIMO-DCNN models perform on a multi-input (i.e. RGB and depth images), multi-trait output lettuce dataset from the 2021 Autonomous Greenhouse Challenge. Ablation studies were conducted to examine the effect of using single versus multiple inputs, and single versus multiple outputs. The MIMO-DCNN model resulted in a normalized mean squared error (NMSE) of 0.068 - a substantial improvement over the top 2021 leaderboard score of 0.081. Open-source code is provided.


title: Distributed Interactive Proofs for the Recognition of Some Geometric Intersection Graph Classes

abstract: A graph $G=(V,E)$ is a geometric intersection graph if every node $v \in V$ is identified with a geometric object of some particular type, and two nodes are adjacent if the corresponding objects intersect. Geometric intersection graph classes have been studied from both the theoretical and practical point of view. On the one hand, many hard problems can be efficiently solved or approximated when the input graph is restricted to a geometric intersection class of graphs. On the other hand, these graphs appear naturally in many applications such as sensor networks, scheduling problems, and others. Recently, in the context of distributed certification and distributed interactive proofs, the recognition of graph classes has started to be intensively studied. Different results related to the recognition of trees, bipartite graphs, bounded diameter graphs, triangle-free graphs, planar graphs, bounded genus graphs, $H$-minor free graphs, etc., have been obtained. The goal of the present work is to design efficient distributed protocols for the recognition of relevant geometric intersection graph classes, namely permutation graphs, trapezoid graphs, circle graphs, and polygon-circle graphs. More precisely, for the two first classes, we give proof labeling schemes recognizing them with logarithmic-sized certificates. For the other two classes, we give three-round distributed interactive protocols that use messages and certificates of size $\mathcal{O}(\log n)$. Finally, we provide logarithmic lower-bounds on the size of the certificates on the proof labeling schemes for the recognition of any of the aforementioned geometric intersection graph classes.


title: On the questions of asymptotic recoverability of information and subsystems in quantum gravity

abstract: A longstanding question in quantum gravity regards the localization of quantum information; one way to formulate this question is to ask how subsystems can be defined in quantum-gravitational systems. The gauge symmetry and necessity of solving the constraints appear to imply that the answers to this question here are different than in finite quantum systems, or in local quantum field theory. Specifically, the constraints can be solved by providing a "gravitational dressing" for the underlying field-theory operators, but this modifies their locality properties. It has been argued that holography itself may be explained through this role of the gauge symmetry and constraints, at the nonperturbative level, but there are also subtleties in constructing a holographic map in this approach. There are also claims that holography is implied even by perturbative solution of the constraints. This short note provides further examination of these questions, and in particular investigates to what extent perturbative or nonperturbative solution of the constraints implies that information naively thought to be localized can be recovered by asymptotic measurements, and the relevance of this in defining subsystems. In the leading perturbative case, the relevant effects are seen to be exponentially suppressed. These questions are, for example, important in sharply characterizing the unitarity problem for black holes.


title: Anchoring to Exemplars for Training Mixture-of-Expert Cell Embeddings

abstract: Analyzing the morphology of cells in microscopy images can provide insights into the mechanism of compounds or the function of genes. Addressing this task requires methods that can not only extract biological information from the images, but also ignore technical variations, ie, changes in experimental procedure or differences between equipments used to collect microscopy images. We propose Treatment ExemplArs with Mixture-of-experts (TEAMs), an embedding learning approach that learns a set of experts that are specialized in capturing technical variations in our training set and then aggregates specialist's predictions at test time. Thus, TEAMs can learn powerful embeddings with less technical variation bias by minimizing the noise from every expert. To train our model, we leverage Treatment Exemplars that enable our approach to capture the distribution of the entire dataset in every minibatch while still fitting into GPU memory. We evaluate our approach on three datasets for tasks like drug discovery, boosting performance on identifying the true mechanism of action of cell treatments by 5.5-11% over the state-of-the-art.


title: Locomotion of magnetoelastic membranes in viscous fluids

abstract: The development of multifunctional and biocompatible microrobots for biomedical applications relies on achieving locomotion through viscous fluids. Here, we describe a framework for swimming in homogeneous magnetoelastic membranes composed of superparamagnetic particles. By solving the equations of motion, we find the dynamical modes of circular membranes in precessing magnetic fields, which are found to actuate in or out of synchronization with a magnetic field precessing above or below a critical precession frequency, $\omega_c$, respectively. For frequencies larger than $\omega_c$, synchronized rotational and radial waves propagate on the membrane. These waves give rise to locomotion in an incompressible fluid at low Reynolds number using the lattice Boltzmann approach. Non-reciprocal motion resulting in swimming is achieved by breaking the morphological symmetry of the membrane, attained via truncation of a circular segment. The membrane translation can be adapted to a predetermined path by programming the external magnetic field. Our results lay the foundation for achieving directed motion in thin, homogeneous magnetoelastic membranes with a diverse array of geometries.


title: Contextual Bandit Applications in Customer Support Bot

abstract: Virtual support agents have grown in popularity as a way for businesses to provide better and more accessible customer service. Some challenges in this domain include ambiguous user queries as well as changing support topics and user behavior (non-stationarity). We do, however, have access to partial feedback provided by the user (clicks, surveys, and other events) which can be leveraged to improve the user experience. Adaptable learning techniques, like contextual bandits, are a natural fit for this problem setting. In this paper, we discuss real-world implementations of contextual bandits (CB) for the Microsoft virtual agent. It includes intent disambiguation based on neural-linear bandits (NLB) and contextual recommendations based on a collection of multi-armed bandits (MAB). Our solutions have been deployed to production and have improved key business metrics of the Microsoft virtual agent, as confirmed by A/B experiments. Results include a relative increase of over 12% in problem resolution rate and relative decrease of over 4% in escalations to a human operator. While our current use cases focus on intent disambiguation and contextual recommendation for support bots, we believe our methods can be extended to other domains.


title: Quantum advantage in biometric authentication with single photons

abstract: It was recently proposed to use the human visual system's ability to perform efficient photon counting in order to devise a new biometric methodology. The relevant biometric fingerprint is represented by the optical losses light suffers along several different paths from the cornea to the retina. The fingerprint is accessed by interrogating a subject on perceiving or not weak light flashes, containing few tens of photons, so that the subject's visual system works at the threshold of perception, at which regime optical losses play a significant role. Here we show that if instead of weak laser light pulses we use quantum light sources, in particular single photon sources, we obtain a quantum advantage, which translates into a reduction of the resources required to achieve a desired performance. Besides the particular application on biometrics, our work further demonstrates that quantum light sources can provide deeper insights when studying human vision.


title: Physically Consistent Neural Networks for building thermal modeling: theory and analysis

abstract: Due to their high energy intensity, buildings play a major role in the current worldwide energy transition. Building models are ubiquitous since they are needed at each stage of the life of buildings, i.e. for design, retrofitting, and control operations. Classical white-box models, based on physical equations, are bound to follow the laws of physics but the specific design of their underlying structure might hinder their expressiveness and hence their accuracy. On the other hand, black-box models are better suited to capture nonlinear building dynamics and thus can often achieve better accuracy, but they require a lot of data and might not follow the laws of physics, a problem that is particularly common for neural network (NN) models. To counter this known generalization issue, physics-informed NNs have recently been introduced, where researchers introduce prior knowledge in the structure of NNs to ground them in known underlying physical laws and avoid classical NN generalization issues. In this work, we present a novel physics-informed NN architecture, dubbed Physically Consistent NN (PCNN), which only requires past operational data and no engineering overhead, including prior knowledge in a linear module running in parallel to a classical NN. We formally prove that such networks are physically consistent - by design and even on unseen data - with respect to different control inputs and temperatures outside and in neighboring zones. We demonstrate their performance on a case study, where the PCNN attains an accuracy up to 40% better than a classical physics-based resistance-capacitance model on 3-day long prediction horizons. Furthermore, despite their constrained structure, PCNNs attain similar performance to classical NNs on the validation data, overfitting the training data less and retaining high expressiveness to tackle the generalization issue.


title: Zero-shot hashtag segmentation for multilingual sentiment analysis

abstract: Hashtag segmentation, also known as hashtag decomposition, is a common step in preprocessing pipelines for social media datasets. It usually precedes tasks such as sentiment analysis and hate speech detection. For sentiment analysis in medium to low-resourced languages, previous research has demonstrated that a multilingual approach that resorts to machine translation can be competitive or superior to previous approaches to the task. We develop a zero-shot hashtag segmentation framework and demonstrate how it can be used to improve the accuracy of multilingual sentiment analysis pipelines. Our zero-shot framework establishes a new state-of-the-art for hashtag segmentation datasets, surpassing even previous approaches that relied on feature engineering and language models trained on in-domain data.


title: Piano Timbre Development Analysis using Machine Learning

abstract: A data set of recorded single played tones of a concert grand piano is investigated using Machine Learning (ML) on psychoacoustic timbre features. The examined instrument has been recorded at two stages: firstly right after manufacture and secondly after being played in a concert hall for one year. A previous study [Plath2019] revealed that listeners clearly distinguished both stages but no clear correlation with acoustics, signal processing tools or verbalizations of perceived differences could be found. Using a Self-Organizing Map (SOM), training single as well as double feature sets, it can be shown that spectral flux is able to perfectly cluster the two stages. Sound Pressure Level (SPL), roughness, and fractal correlation dimension (as a measure for initial transient chaoticity) are furthermore able to order the keys with respect to high and low notes. Combining spectral flux with the three other features in double-feature training sets maintains stage clustering only for SPL and fractal dimension, showing sub-clusters for both stages. These sub-clusters point to a homogenization of SPL for stage 2 with respect to stage 1 and a pronounced ordering and sub-clustering of key regions with respect to initial transient chaoticity.


title: Multi-scale Feature Learning Dynamics: Insights for Double Descent

abstract: A key challenge in building theoretical foundations for deep learning is the complex optimization dynamics of neural networks, resulting from the high-dimensional interactions between the large number of network parameters. Such non-trivial dynamics lead to intriguing behaviors such as the phenomenon of "double descent" of the generalization error. The more commonly studied aspect of this phenomenon corresponds to model-wise double descent where the test error exhibits a second descent with increasing model complexity, beyond the classical U-shaped error curve. In this work, we investigate the origins of the less studied epoch-wise double descent in which the test error undergoes two non-monotonous transitions, or descents as the training time increases. By leveraging tools from statistical physics, we study a linear teacher-student setup exhibiting epoch-wise double descent similar to that in deep neural networks. In this setting, we derive closed-form analytical expressions for the evolution of generalization error over training. We find that double descent can be attributed to distinct features being learned at different scales: as fast-learning features overfit, slower-learning features start to fit, resulting in a second descent in test error. We validate our findings through numerical experiments where our theory accurately predicts empirical findings and remains consistent with observations in deep neural networks.


title: Almost Everything About the Unitary Almost Mathieu Operator

abstract: We introduce a unitary almost-Mathieu operator, which is obtained from a two-dimensional quantum walk in a uniform magnetic field. We exhibit a version of Aubry--Andr\'{e} duality for this model, which partitions the parameter space into three regions: a supercritical region and a subcritical region that are dual to one another, and a critical regime that is self-dual. In each parameter region, we characterize the cocycle dynamics of the transfer matrix cocycle generated by the associated generalized eigenvalue equation. In particular, we show that supercritical, critical, and subcritical behavior all occur in this model. Using Avila's global theory of one-frequency cocycles, we exactly compute the Lyapunov exponent on the spectrum in terms of the given parameters. We also characterize the spectral type for each value of the coupling constant, almost every frequency, and almost every phase. Namely, we show that for almost every frequency and every phase the spectral type is purely absolutely continuous in the subcritical region, pure point in the supercritical region, and purely singular continuous in the critical region. In some parameter regions, we refine the almost-sure results. In the critical case for instance, we show that the spectrum is a Cantor set of zero Lebesgue measure for arbitrary irrational frequency and that the spectrum is purely singular continuous for all but countably many phases.

categories: math.SP math-ph math.FA math.MP quant-ph


title: Minimax properties of Dirichlet kernel density estimators

abstract: This paper considers the asymptotic behavior in $\beta$-H\"older spaces, and under $L^p$ losses, of a Dirichlet kernel density estimator proposed by Aitchison and Lauder (1985) for the analysis of compositional data. In recent work, Ouimet and Tolosana-Delgado (2022) established the uniform strong consistency and asymptotic normality of this estimator. As a complement, it is shown here that the Aitchison-Lauder estimator can achieve the minimax rate asymptotically for a suitable choice of bandwidth whenever $(p,\beta) \in [1, 3) \times (0, 2]$ or $(p, \beta) \in \mathcal{A}_d$, where $\mathcal{A}_d$ is a specific subset of $[3, 4) \times (0, 2]$ that depends on the dimension $d$ of the Dirichlet kernel. It is also shown that this estimator cannot be minimax when either $p \in [4, \infty)$ or $\beta \in (2, \infty)$. These results extend to the multivariate case, and also rectify in a minor way, earlier findings of Bertin and Klutchnikoff (2011) concerning the minimax properties of Beta kernel estimators.


title: Modeling synchronization in human musical rhythms using Impulse Pattern Formulation (IPF)

abstract: When musicians perform in an ensemble, synchronizing to a mutual pace is the foundation of their musical interaction. Clock generators, e.g., metronomes, or drum machines, might assist such synchronization, but these means, in general, will also distort this natural, self-organized, inter-human synchronization process. In this work, the synchronization of musicians to an external rhythm is modeled using the Impulse Pattern Formulation (IPF), an analytical modeling approach for synergetic systems motivated by research on musical instruments. Nonlinear coupling of system components is described as the interaction of individually propagating and exponentially damped impulse trains. The derived model is systematically examined by analyzing its behavior when coupled to numerical designed and carefully controlled rhythmical beat sequences. The results are evaluated by comparison in the light of other publications on tapping. Finally, the IPF model can be applied to analyze the personal rhythmical signature of specific musicians or to replace drum machines and click tracks with more musical and creative solutions.


title: On some locally symmetric embedded spaces with non-negative scalar curvature and their characterization

abstract: In this work we perform a general study of properties of a class of locally symmetric embedded hypersurfaces in spacetimes admitting a $1+1+2$ spacetime decomposition. The hypersurfaces are given by specifying the form of the Ricci tensor with respect to the induced metric. These are slices of constant time in the spacetime. Firstly, the form of the Ricci tensor for general hypersurfaces is obtained and the conditions under which the general case reduces to those of constant time slices are specified. We provide a characterization of these hypersurfaces, with key physical quantities in the spacetime playing a role in specifying the local geometry of these hypersurfaces. Furthermore, we investigate the case where these hypersurfaces admit a Ricci soliton structure. The particular cases where the vector fields associated to the solitons are Killing or conformal Killing vector fields are analyzed. Finally, in the context of spacetimes with local rotational symmetry it is shown that, only spacetimes in this class with vanishing rotation and spatial twist can admit the hypersurface types considered, and that the hypersurfaces are necessarily flat. And if such hypersurface do admit a Ricci soliton structure, the soliton is steady, with the components of the soliton field being constants.


title: Cross-validation for change-point regression: pitfalls and solutions

abstract: Cross-validation is the standard approach for tuning parameter selection in many non-parametric regression problems. However its use is less common in change-point regression, perhaps as its prediction error-based criterion may appear to permit small spurious changes and hence be less well-suited to estimation of the number and location of change-points. We show that in fact the problems of cross-validation with squared error loss are more severe and can lead to systematic under- or over-estimation of the number of change-points, and highly suboptimal estimation of the mean function in simple settings where changes are easily detectable. We propose two simple approaches to remedy these issues, the first involving the use of absolute error rather than squared error loss, and the second involving modifying the holdout sets used. For the latter, we provide conditions that permit consistent estimation of the number of change-points for a general change-point estimation procedure. We show these conditions are satisfied for least squares estimation using new results on its performance when supplied with the incorrect number of change-points. Numerical experiments show that our new approaches are competitive with common change-point methods using classical tuning parameter choices when error distributions are well-specified, but can substantially outperform these in misspecified models. An implementation of our methodology is available in the R package crossvalidationCP on CRAN.


title: Text2Mesh: Text-Driven Neural Stylization for Meshes

abstract: In this work, we develop intuitive controls for editing the style of 3D objects. Our framework, Text2Mesh, stylizes a 3D mesh by predicting color and local geometric details which conform to a target text prompt. We consider a disentangled representation of a 3D object using a fixed mesh input (content) coupled with a learned neural network, which we term neural style field network. In order to modify style, we obtain a similarity score between a text prompt (describing style) and a stylized mesh by harnessing the representational power of CLIP. Text2Mesh requires neither a pre-trained generative model nor a specialized 3D mesh dataset. It can handle low-quality meshes (non-manifold, boundaries, etc.) with arbitrary genus, and does not require UV parameterization. We demonstrate the ability of our technique to synthesize a myriad of styles over a wide variety of 3D meshes.


title: On Complexity of 1-Center in Various Metrics

abstract: We consider the classic 1-center problem: Given a set $P$ of $n$ points in a metric space find the point in $P$ that minimizes the maximum distance to the other points of $P$. We study the complexity of this problem in $d$-dimensional $\ell_p$-metrics and in edit and Ulam metrics over strings of length $d$. Our results for the 1-center problem may be classified based on $d$ as follows. $\bullet$ Small $d$: Assuming the hitting set conjecture (HSC), we show that when $d=\omega(\log n)$, no subquadratic algorithm can solve 1-center problem in any of the $\ell_p$-metrics, or in edit or Ulam metrics. $\bullet$ Large $d$: When $d=\Omega(n)$, we extend our conditional lower bound to rule out subquartic algorithms for 1-center problem in edit metric (assuming Quantified SETH). On the other hand, we give a $(1+\epsilon)$-approximation for 1-center in Ulam metric with running time $\tilde{O_{\varepsilon}}(nd+n^2\sqrt{d})$. We also strengthen some of the above lower bounds by allowing approximations or by reducing the dimension $d$, but only against a weaker class of algorithms which list all requisite solutions. Moreover, we extend one of our hardness results to rule out subquartic algorithms for the well-studied 1-median problem in the edit metric, where given a set of $n$ strings each of length $n$, the goal is to find a string in the set that minimizes the sum of the edit distances to the rest of the strings in the set.


title: Context-Aware Transfer Attacks for Object Detection

abstract: Blackbox transfer attacks for image classifiers have been extensively studied in recent years. In contrast, little progress has been made on transfer attacks for object detectors. Object detectors take a holistic view of the image and the detection of one object (or lack thereof) often depends on other objects in the scene. This makes such detectors inherently context-aware and adversarial attacks in this space are more challenging than those targeting image classifiers. In this paper, we present a new approach to generate context-aware attacks for object detectors. We show that by using co-occurrence of objects and their relative locations and sizes as context information, we can successfully generate targeted mis-categorization attacks that achieve higher transfer success rates on blackbox object detectors than the state-of-the-art. We test our approach on a variety of object detectors with images from PASCAL VOC and MS COCO datasets and demonstrate up to $20$ percentage points improvement in performance compared to the other state-of-the-art methods.


title: Extensions of quasidiagonal $C^*$-algebras and controlling the $K_0$-map of embeddings

abstract: We study the validity of the Blackadar-Kirchberg conjecture for extensions of separable, nuclear, quasidiagonal $C^*$-algebras that satisfy the UCT. More specifically, we show that the conjecture for the extension has an affirmative answer if the ideal lies in a class of $C^*$-algebras that is closed under local approximations and contains all separable ASH algebras, as well as certain classes of simple, unital $C^*$-algebras and crossed products of unital $C^*$-algebras with the integers.


title: Interior Structure Models Of Venus

abstract: Venus' mass and radius are similar to those of Earth. However, dissimilarities in atmospheric properties, geophysical activity and magnetic field generation could hint towards significant differences in the chemical composition and interior evolution of the two planets. Although various explanations for the differences between Venus and Earth have been proposed, the currently available data are insufficient to discriminate among the different solutions. Here we investigate the possible range of Venus structure models. We assume that core segregation happened as a single-stage event. The mantle composition is inferred from the core composition using a prescription for metal-silicate partitioning. We consider three different cases for the composition of Venus defined via the bulk Si and Mg content, and the core's S content. Permissible ranges for the core size, mantle and core composition as well as the normalized moment of inertia (MoI) are presented for these compositions. A solid inner core could exist for all compositions. We estimate that Venus' MoI is 0.317-0.351 and its core size 2930-4350 km for all assumed compositions. Higher MoI values correspond to more oxidizing conditions during core segregation. A determination of the abundance of FeO in Venus' mantle by future missions could further constrain its composition and internal structure. This can reveal important information on Venus' formation and evolution, and possibly, the reasons for the differences between Venus and our home planet.


title: Quadratic points on bielliptic modular curves

abstract: Bruin and Najman, Ozman and Siksek, and Box described all the quadratic points on the modular curves of genus $2\leq g(X_0(n)) \leq 5$. Since all the hyperelliptic curves $X_0(n)$ are of genus $\leq 5$ and as a curve can have infinitely many quadratic points only if it is either of genus $\leq 1$, hyperelliptic or bielliptic, the question of describing the quadratic points on the bielliptic modular curves $X_0(n)$ naturally arises; this question has recently also been posed by Mazur. We answer Mazur's question completely and describe the quadratic points on all the bielliptic modular curves $X_0(n)$ for which this has not been done already. The values of $n$ that we deal with are $n=60,62,69,79,83,89,92,94,95,101,119$ and $131$; the curves $X_0(n)$ are of genus up to $11$. We find all the exceptional points on these curves and show that they all correspond to CM elliptic curves. The two main methods we use are Box's relative symmetric Chabauty method and an application of a moduli description of $\Q$-curves of degree $d$ with an independent isogeny of degree $m$, which reduces the problem to finding the rational points on several quotients of modular curves.


title: CALVIN: A Benchmark for Language-Conditioned Policy Learning for Long-Horizon Robot Manipulation Tasks

abstract: General-purpose robots coexisting with humans in their environment must learn to relate human language to their perceptions and actions to be useful in a range of daily tasks. Moreover, they need to acquire a diverse repertoire of general-purpose skills that allow composing long-horizon tasks by following unconstrained language instructions. In this paper, we present CALVIN (Composing Actions from Language and Vision), an open-source simulated benchmark to learn long-horizon language-conditioned tasks. Our aim is to make it possible to develop agents that can solve many robotic manipulation tasks over a long horizon, from onboard sensors, and specified only via human language. CALVIN tasks are more complex in terms of sequence length, action space, and language than existing vision-and-language task datasets and supports flexible specification of sensor suites. We evaluate the agents in zero-shot to novel language instructions and to novel environments and objects. We show that a baseline model based on multi-context imitation learning performs poorly on CALVIN, suggesting that there is significant room for developing innovative agents that learn to relate human language to their world models with this benchmark.


title: Uniform even subgraphs and graphical representations of Ising as factors of i.i.d

abstract: We prove that the Loop O(1) model, a well-known graphical expansion of the Ising model, is a factor of i.i.d. on unimodular random rooted graphs under various conditions, including in the presence of a non-negative external field. As an application we show that the gradient of the free Ising model is a factor of i.i.d. on unimodular planar maps having a locally finite dual. The key idea is to develop an appropriate theory of local limits of uniform even subgraphs with various boundary conditions and prove that they can be sampled as a factor of i.i.d. Another key tool we prove and exploit is that the wired uniform spanning tree on a unimodular transient graph is a factor of i.i.d. This partially answers some questions posed by Hutchcroft.


title: SIMD-Optimized Search Over Sorted Data

abstract: Applications often require a fast, single-threaded search algorithm over sorted data, typical in table-lookup operations. We explore various search algorithms for a large number of search candidates over a relatively small array of logarithmically-distributed sorted data. These include an innovative hash-based search that takes advantage of floating point representation to bin data by the exponent. Algorithms that can be optimized to take advantage of SIMD vector instructions are of particular interest. We then conduct a case study applying our results and analyzing algorithmic performance with the EOSPAC package. EOSPAC is a table look-up library for manipulation and interpolation of SESAME equation-of-state data. Our investigation results in a couple of algorithms with better performance with a best case 8x speedup over the original EOSPAC Hunt-and-Locate implementation. Our techniques are generalizable to other instances of search algorithms seeking to get a performance boost from vectorization.


title: Traversing Time with Multi-Resolution Gaussian Process State-Space Models

abstract: Gaussian Process state-space models capture complex temporal dependencies in a principled manner by placing a Gaussian Process prior on the transition function. These models have a natural interpretation as discretized stochastic differential equations, but inference for long sequences with fast and slow transitions is difficult. Fast transitions need tight discretizations whereas slow transitions require backpropagating the gradients over long subtrajectories. We propose a novel Gaussian process state-space architecture composed of multiple components, each trained on a different resolution, to model effects on different timescales. The combined model allows traversing time on adaptive scales, providing efficient inference for arbitrarily long sequences with complex dynamics. We benchmark our novel method on semi-synthetic data and on an engine modeling task. In both experiments, our approach compares favorably against its state-of-the-art alternatives that operate on a single time-scale only.


title: Pearson's goodness-of-fit tests for sparse distributions

abstract: Pearson's chi-squared test is widely used to test the goodness of fit between categorical data and a given discrete distribution function. When the number of sets of the categorical data, say $k$, is a fixed integer, Pearson's chi-squared test statistic converges in distribution to a chi-squared distribution with $k-1$ degrees of freedom when the sample size $n$ goes to infinity. In real applications, the number $k$ often changes with $n$ and may be even much larger than $n$. By using the martingale techniques, we prove that Pearson's chi-squared test statistic converges to the normal under quite general conditions. We also propose a new test statistic which is more powerful than chi-squared test statistic based on our simulation study. A real application to lottery data is provided to illustrate our methodology.


title: A Risk-Averse Preview-based $Q$-Learning Algorithm: Application to Highway Driving of Autonomous Vehicles

abstract: A risk-averse preview-based $Q$-learning planner is presented for navigation of autonomous vehicles. To this end, the multi-lane road ahead of a vehicle is represented by a finite-state non-stationary Markov decision process (MDP). A risk assessment unit module is then presented that leverages the preview information provided by sensors along with a stochastic reachability module to assign reward values to the MDP states and update them as scenarios develop. A sampling-based risk-averse preview-based Q-learning algorithm is finally developed that generates samples using the preview information and reward function to learn risk-averse optimal planning strategies without actual interaction with the environment. The risk factor is imposed on the objective function to avoid fluctuation of the Q values, which can jeopardize the vehicle's safety and/or performance. The overall hybrid automaton model of the system is leveraged to develop a feasibility check unit module that detects unfeasible plans and enables the planner system to react proactively to the changes of the environment. Finally, to verify the efficiency of the presented algorithm, its implementation on two highway driving scenarios of an autonomous vehicle in a varying traffic density is considered.


title: EPR Pairs from Indefinite Causal Order

abstract: We develop a novel protocol to achieve entanglement creation between qubit pairs by making a qubit pair experience an indefinite causal order process. In particular, our protocol does not need two qubits to interact with each other, meaning that only local operations are applied on individual qubits. In order to demonstrate the usage of the protocol, we first identify the conditions under which, EPR pairs, the maximally entangled states emerge from our protocol. Furthermore, we show that, although our protocol is not deterministic, the success probability of getting an EPR pair could be non-vanishing, even comparably high enough. Our result not only sheds new light on entanglement creation, but revealing that a quantum switch has further potential to be explored.


title: The InSight HP$^3$ mole on Mars: Lessons learned from attempts to penetrate to depth in the Martian soil

abstract: The NASA InSight mission payload includes the Heat Flow and Physical Properties Package HP$^3$ to measure the surface heat flow. The package was designed to use a small penetrator -- nicknamed the mole -- to implement a string of temperature sensors in the soil to a depth of 5m. The mole itself is equipped with sensors to measure a thermal conductivity as it proceeds to depth. The heat flow would be calculated from the product of the temperature gradient and the thermal conductivity. To avoid the perturbation caused by annual surface temperature variations, the measurements would be taken at a depth between 3 m and 5 m. The mole was designed to penetrate cohesionless soil similar to Quartz sand which was expected to provide a good analogue material for Martian sand. The sand would provide friction to the buried mole hull to balance the remaining recoil of the mole hammer mechanism that drives the mole forward. Unfortunately, the mole did not penetrate more than a mole length of 40 cm. The failure to penetrate deeper was largely due to a few tens of centimeter thick cohesive duricrust that failed to provide the required friction. Although a suppressor mass and spring in the hammer mechanism absorbed much of the recoil, the available mass did not allow a system that would have eliminated the recoil. The mole penetrated to 40 cm depth benefiting from friction provided by springs in the support structure from which it was deployed. It was found in addition that the Martian soil provided unexpected levels of penetration resistance that would have motivated to designing a more powerful mole. It is concluded that more mass would have allowed to design a more robust system with little or no recoil, more energy of the mole hammer mechanism and a more massive support structure.


title: Simulation Intelligence: Towards a New Generation of Scientific Methods

abstract: The original "Seven Motifs" set forth a roadmap of essential methods for the field of scientific computing, where a motif is an algorithmic method that captures a pattern of computation and data movement. We present the "Nine Motifs of Simulation Intelligence", a roadmap for the development and integration of the essential algorithms necessary for a merger of scientific computing, scientific simulation, and artificial intelligence. We call this merger simulation intelligence (SI), for short. We argue the motifs of simulation intelligence are interconnected and interdependent, much like the components within the layers of an operating system. Using this metaphor, we explore the nature of each layer of the simulation intelligence operating system stack (SI-stack) and the motifs therein: (1) Multi-physics and multi-scale modeling; (2) Surrogate modeling and emulation; (3) Simulation-based inference; (4) Causal modeling and inference; (5) Agent-based modeling; (6) Probabilistic programming; (7) Differentiable programming; (8) Open-ended optimization; (9) Machine programming. We believe coordinated efforts between motifs offers immense opportunity to accelerate scientific discovery, from solving inverse problems in synthetic biology and climate science, to directing nuclear energy experiments and predicting emergent behavior in socioeconomic settings. We elaborate on each layer of the SI-stack, detailing the state-of-art methods, presenting examples to highlight challenges and opportunities, and advocating for specific ways to advance the motifs and the synergies from their combinations. Advancing and integrating these technologies can enable a robust and efficient hypothesis-simulation-analysis type of scientific method, which we introduce with several use-cases for human-machine teaming and automated science.


title: A criterion of quasi-infinite divisibility for discrete laws

abstract: We consider arbitrary discrete probability laws on the real line. We obtain a criterion of their belonging to a new class of quasi-infinitely divisible laws, which is a wide natural extension of the class of well known infinitely divisible laws through the L\'evy type representations.


title: From Coarse to Fine-grained Concept based Discrimination for Phrase Detection

abstract: Phrase detection requires methods to identify if a phrase is relevant to an image and localize it, if applicable. A key challenge for training more discriminative detection models is sampling negatives. Sampling techniques from prior work focus primarily on hard, often noisy, negatives disregarding the broader distribution of negative samples. Our proposed CFCD-Net addresses this through two novels methods. First, we generate groups of semantically similar words we call concepts (\eg, \{dog, cat, horse\} and \ \{car, truck, SUV\}), and then train our CFCD-Net to discriminate between a region of interest and its unrelated concepts. Second, for phrases containing fine-grained mutually-exclusive words (\eg, colors), we force the model to select only one applicable phrase for each region using our novel fine-grained module (FGM). We evaluate our approach on Flickr30K Entities and RefCOCO+, where we improve mAP over the state-of-the-art by 1.5-2 points. When considering only the phrases affected by our FGM module, we improve by 3-4 points on both datasets.


title: Deep brain stimulation with a computational model for the cortex-thalamus-basal-ganglia system and network dynamics of neurological disorders

abstract: Deep brain stimulation (DBS) can alleviate the movement disorders like Parkinson's disease (PD). Indeed, it is known that aberrant beta (13-30Hz)oscillations and the loss of dopaminergic neurons in the basal ganglia-thalamus (BGTH) and cortex characterize the akinesia symptoms of PD. However, the relevant biophysical mechanism behind this process still remains unclear. Based on the prior striatal inhibitory model, we propose an extended BGTH model incorporating medium spine neurons (MSNs) and fast-spiking interneurons (FSIs) along with the effect of DBS. We are focusing in this paper on an open-loop DBS mode, where the stimulation parameters stay constant independent of variations in the disease state, and modifications of parameters rely mainly on trial and error of medical experts. Additionally, we propose a novel combined model of the cerebellar-basal-ganglia thalamocortical network, MSNs, and FSIs, and show new results that indicate that Parkinsonian oscillations in the beta-band frequency range emerge from the dynamics of such a network. Our model predicts that DBS can be used to suppress beta oscillations in globus pallidus pars interna (GPi) neurons. This research will help our better understanding of the changes in brain activity caused by DBS, providing new insight for studying PD in the future.


title: Approximations for STERGMs Based on Cross-Sectional Data

abstract: Temporal exponential-family random graph models (TERGMs) are a flexible class of network models for the dynamics of tie formation and dissolution. In practice, separable TERGMs (STERGMs) are the subclass most often used, as these permit estimation from inexpensive cross-sectional study designs, and benefit from approximations designed to reduce the computational burden. Improving the approximations are the focus of this paper. We extend the work of Carnegie et al., which addressed the problem of constructing a STERGM with two specific equilibrium properties: a cross-sectional distribution defined by a given exponential-family random graph model (ERGM), and tie durations defined by given constant hazards of dissolution. We start with Carnegie et al.'s observation that the exact result is tractable in the dyad-independent case, and then show that taking the sparse limit of the exact result leads to a different approximation than the one they presented. We show that the new approximation outperforms theirs for sparse, dyad-independent models, and that for dyad-dependent models the errors tend to increase with the level of dependence for both approximations. We then extend the theoretical results of Carnegie et al. to the dyad-dependent case, proving that both the old and new approximations are asymptotically exact as the STERGM time step size goes to zero, for arbitrary dyad-dependent terms and some dyad-dependent constraints. We also show that the continuous-time limit of the discrete-time approximations has exactly the combination of cross-sectional and durational equilibrium behavior that we seek.


title: Analyzing a Carceral Algorithm used by the Pennsylvania Department of Corrections

abstract: Scholars have focused on algorithms used during sentencing, bail, and parole, but little work explores what we call carceral algorithms that are used during incarceration. This paper is focused on the Pennsylvania Additive Classification Tool (PACT) used to classify prisoners' custody levels while they are incarcerated. Algorithms that are used during incarceration warrant deeper attention by scholars because they have the power to enact the lived reality of the prisoner. The algorithm in this case determines the likelihood a person would endure additional disciplinary actions, can complete required programming, and gain experiences that, among other things, are distilled into variables feeding into the parole algorithm. Given such power, examining algorithms used on people currently incarcerated offers a unique analytic view to think about the dialectic relationship between data and algorithms. Our examination of the PACT is two-fold and complementary. First, our qualitative overview of the historical context surrounding PACT reveals that it is designed to prioritize incapacitation and control over rehabilitation. While it closely informs prisoner rehabilitation plans and parole considerations, it is rooted in population management for prison securitization. Second, on analyzing data for 146,793 incarcerated people in PA, along with associated metadata related to the PACT, we find it is replete with racial bias as well as errors, omissions, and inaccuracies. Our findings to date further caution against data-driven criminal justice reforms that rely on pre-existing data infrastructures and expansive, uncritical, data-collection routines.


title: Unsupervised Domain Adaptation for Semantic Image Segmentation: a Comprehensive Survey

abstract: Semantic segmentation plays a fundamental role in a broad variety of computer vision applications, providing key information for the global understanding of an image. Yet, the state-of-the-art models rely on large amount of annotated samples, which are more expensive to obtain than in tasks such as image classification. Since unlabelled data is instead significantly cheaper to obtain, it is not surprising that Unsupervised Domain Adaptation reached a broad success within the semantic segmentation community. This survey is an effort to summarize five years of this incredibly rapidly growing field, which embraces the importance of semantic segmentation itself and a critical need of adapting segmentation models to new environments. We present the most important semantic segmentation methods; we provide a comprehensive survey on domain adaptation techniques for semantic segmentation; we unveil newer trends such as multi-domain learning, domain generalization, test-time adaptation or source-free domain adaptation; we conclude this survey by describing datasets and benchmarks most widely used in semantic segmentation research. We hope that this survey will provide researchers across academia and industry with a comprehensive reference guide and will help them in fostering new research directions in the field.


title: Aspect Ratio Universal Rectangular Layouts

abstract: A \emph{generic rectangular layout} (for short, \emph{layout}) is a subdivision of an axis-aligned rectangle into axis-aligned rectangles, no four of which have a point in common. Such layouts are used in data visualization and in cartography. The contacts between the rectangles represent semantic or geographic relations. A layout is weakly (strongly) \emph{aspect ratio universal} if any assignment of aspect ratios to rectangles can be realized by a weakly (strongly) equivalent layout. We give combinatorial characterizations for weakly and strongly aspect ratio universal layouts. Furthermore, we describe a quadratic-time algorithm that decides whether a given graph is the dual graph of a strongly aspect ratio universal layout, and finds such a layout if one exists.


title: Input-level Inductive Biases for 3D Reconstruction

abstract: Much of the recent progress in 3D vision has been driven by the development of specialized architectures that incorporate geometrical inductive biases. In this paper we tackle 3D reconstruction using a domain agnostic architecture and study how instead to inject the same type of inductive biases directly as extra inputs to the model. This approach makes it possible to apply existing general models, such as Perceivers, on this rich domain, without the need for architectural changes, while simultaneously maintaining data efficiency of bespoke models. In particular we study how to encode cameras, projective ray incidence and epipolar geometry as model inputs, and demonstrate competitive multi-view depth estimation performance on multiple benchmarks.


title: Projection methods for Neural Field equations

abstract: Neural field models are nonlinear integro-differential equations for the evolution of neuronal activity, and they are a prototypical large-scale, coarse-grained neuronal model in continuum cortices. Neural fields are often simulated heuristically and, in spite of their popularity in mathematical neuroscience, their numerical analysis is not yet fully established. We introduce generic projection methods for neural fields, and derive a-priori error bounds for these schemes. We extend an existing framework for stationary integral equations to the time-dependent case, which is relevant for neuroscience applications. We find that the convergence rate of a projection scheme for a neural field is determined to a great extent by the convergence rate of the projection operator. This abstract analysis, which unifies the treatment of collocation and Galerkin schemes, is carried out in operator form, without resorting to quadrature rules for the integral term, which are introduced only at a later stage, and whose choice is enslaved by the choice of the projector. Using an elementary timestepper as an example, we demonstrate that the error in a time stepper has two separate contributions: one from the projector, and one from the time discretisation. We give examples of concrete projection methods: two collocation schemes (piecewise-linear and spectral collocation) and two Galerkin schemes (finite elements and spectral Galerkin); for each of them we derive error bounds from the general theory, introduce several discrete variants, provide implementation details, and present reproducible convergence tests.


title: GAM Changer: Editing Generalized Additive Models with Interactive Visualization

abstract: Recent strides in interpretable machine learning (ML) research reveal that models exploit undesirable patterns in the data to make predictions, which potentially causes harms in deployment. However, it is unclear how we can fix these models. We present our ongoing work, GAM Changer, an open-source interactive system to help data scientists and domain experts easily and responsibly edit their Generalized Additive Models (GAMs). With novel visualization techniques, our tool puts interpretability into action -- empowering human users to analyze, validate, and align model behaviors with their knowledge and values. Built using modern web technologies, our tool runs locally in users' computational notebooks or web browsers without requiring extra compute resources, lowering the barrier to creating more responsible ML models. GAM Changer is available at https://interpret.ml/gam-changer.


title: Using Evolutionary Algorithms to Design Antennas with Greater Sensitivity to Ultra High Energy Neutrinos

abstract: The Genetically Evolved NEutrino Telescopes for Improved Sensitivity, or GENETIS, project seeks to optimize detectors in physics for science outcomes in high dimensional parameter spaces. In this project, we designed an antenna using a genetic algorithm with a science outcome directly as the sole figure of merit. This paper presents initial results on the improvement of an antenna design for in ice neutrino detectors using the current Askaryan Radio Array, or ARA, experiment as a baseline. By optimizing for the effective volume using the evolved antenna design in ARA, we improve upon ARAs simulated sensitivity to ultra high energy neutrinos by 22 percent, despite using limited parameters in this initial investigation. Future improvements will continue to increase the computational efficiency of the genetic algorithm and the complexity and fitness of the antenna designs. This work lays the foundation for continued research and development of methods to increase the sensitivity of detectors in physics and other fields in parameter spaces of high dimensionality.


title: A Comparison of Multi-Phase Magnetic Field Tracers in a High-Galactic Latitude Region of the Filamentary Interstellar Medium

abstract: Understanding how the Galactic magnetic field threads the multi-phase interstellar medium (ISM) remains a considerable challenge, as different magnetic field tracers probe dissimilar phases and field components. We search for evidence of a common magnetic field shared between the ionized and neutral ISM by comparing 1.4 GHz radio continuum polarization and HI line emission from the Galactic Arecibo L-Band Feed Array Continuum Transit Survey (GALFACTS) and Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey, respectively. We compute the polarization gradient of the continuum emission and search for associations with diffuse/translucent HI structures. The polarization gradient is sensitive to changes in the integrated product of the thermal electron density and line-of-sight field strength ($B_\parallel$) in warm ionized gas, while narrow HI structures highlight the plane-of-sky field orientation in cold neutral gas. We identified one region in the high-Galactic latitude Arecibo sky, G216+26 centered on $(\ell,b)\sim(216\deg,+26\deg)$, containing filaments in the polarization gradient that are aligned with narrow HI structures roughly parallel to the Galactic plane. We present a comparison of multi-phase observations and magnetic field tracers of this region, demonstrating that the warm ionized and cold neutral media are connected likely via a common magnetic field. We quantify the physical properties of a polarization gradient filament associated with H$\alpha$ emission, measuring a line-of-sight field strength $B_\parallel=6{\pm}4 \mu$G and a plasma beta $\beta=2.1^{+3.1}_{-2.1}$. We discuss the lack of widespread multi-phase magnetic field alignments and consider whether this region is associated with a short-timescale or physically rare phenomenon. This work highlights the utility of multi-tracer analyses for understanding the magnetized ISM.


title: On the Stability of Additional Moons Orbiting Kepler-1625 b

abstract: Since it was proposed the exomoon candidate Kepler-1625 b-I changed the way we see satellite systems. Because of its unusual physical characteristics, many questions about the stability and origin of this candidate were raised. Currently, we have enough theoretical studies to assure that if Kepler-1625 b-I is indeed confirmed, it will be stable. The origin of this candidate was also explored. Previous works indicated that the most likely scenario is capture, even though conditions for in situ formation were also investigated. In this work, we assume that Kepler-1625 b-I is an exomoon and studied the possibility of an additional, massive exomoon being stable in the same system. To model this scenario we perform N-body simulations of a system including the planet, Kepler-1625 b-I and one extra Earth-like satellite. Based on previous results, the satellites in our system will be exposed to tidal interactions with the planet and gravitation effects due to the rotation of the planet. We found that the satellite system around Kepler-1625 b is capable of harbouring two massive satellites. The extra Earth-like satellite would be stable in different locations between the planet and Kepler-1625 b-I, with a preference for regions inside $25$ $R_p$. Our results suggest that the strong tidal interactions between the planet and the satellites is an important mechanism to assure the stability of satellites in circular orbits closer to the planet, while the 2:1 mean motion resonance between the Earth-like satellite and Kepler-1625 b-I would provide stability for satellites in wider orbits.


title: Assessing models of force-dependent unbinding rates via infrequent metadynamics

abstract: Protein-ligand interactions are crucial for a wide range of physiological processes. Many cellular functions result in these non-covalent `bonds' being mechanically strained, and this can be integral to proper cellular function. Broadly, two classes of force dependence have been observed -- slip bonds, where unbinding rate increases, and catch bonds where unbinding rate decreases. Despite much theoretical work, we cannot we predict for which protein-ligand pairs, pulling coordinates, and forces a particular rate dependence will appear. Here, we assess the ability of MD simulations combined with enhanced sampling techniques to probe the force dependence of unbinding rates. We show that the infrequent metadynamics technique correctly produces both catch and slip bonding kinetics for model potentials. We then apply it to the well-studied case of a buckyball in a hydrophobic cavity, which appears to exhibit an ideal slip bond. Finally, we compute the force-dependent unbinding rate of biotin-streptavidin. Here, the complex nature of the unbinding process causes the infrequent metadynamics method to begin to break down due to the presence of unbinding intermediates, despite use of a previously optimized sampling coordinate. Allowing for this limitation, a combination of kinetic and free energy computations predict an overall slip bond for larger forces consistent with prior experimental results, although there are substantial deviations at small forces that require further investigation. This work demonstrates the promise of predicting force-dependent unbinding rates using enhanced sampling MD techniques, while also revealing the methodological barriers that must be overcome to tackle more complex targets in the future.


title: Design of Planar and Conformal, Passive, Loss-less Metasurfaces that Beamform

abstract: A general technique for synthesizing both planar and conformal beamforming metasurfaces is presented that utilizes full-wave modeling techniques and rapid optimization methods. The synthesized metasurfaces consist of a patterned metallic cladding supported by a finite-size grounded dielectric substrate. The metasurfaces are modeled using integral equations which accurately account for mutual coupling and the metasurface's finite dimensions. The synthesis technique consists of three phases: a direct solve phase to obtain an initial metasurface design with complex-valued impedances satisfying the desired far-field beam specifications, a subsequent optimization phase that converts the complex-valued impedances to purely reactive ones, and a final patterning phase to realize the purely reactive impedances as a patterned metallic cladding. The optimization phase introduces surface waves which facilitate passivity. The metasurface is optimized using gradient descent with a semi-analytic gradient obtained using the adjoint variable method. Three examples are presented: a low-profile directly-fed metasurface antenna with near perfect aperture efficiency, a scanned-beam reflectarray design with controlled sidelobes, and a conformal metasurface reflectarray. The far-field and near-field performance of the metasurfaces are verified and the bandwidth and loss tolerance of the metasurfaces are investigated.


title: Computing a Link Diagram from its Exterior

abstract: A knot is a circle piecewise-linearly embedded into the 3-sphere. The topology of a knot is intimately related to that of its exterior, which is the complement of an open regular neighborhood of the knot. Knots are typically encoded by planar diagrams, whereas their exteriors, which are compact 3-manifolds with torus boundary, are encoded by triangulations. Here, we give the first practical algorithm for finding a diagram of a knot given a triangulation of its exterior. Our method applies to links as well as knots, allows us to recover links with hundreds of crossings. We use it to find the first diagrams known for 23 principal congruence arithmetic link exteriors; the largest has over 2,500 crossings. Other applications include finding pairs of knots with the same 0-surgery, which relates to questions about slice knots and the smooth 4D Poincar\'e conjecture.


title: CSG0: Continual Urban Scene Generation with Zero Forgetting

abstract: With the rapid advances in generative adversarial networks (GANs), the visual quality of synthesised scenes keeps improving, including for complex urban scenes with applications to automated driving. We address in this work a continual scene generation setup in which GANs are trained on a stream of distinct domains; ideally, the learned models should eventually be able to generate new scenes in all seen domains. This setup reflects the real-life scenario where data are continuously acquired in different places at different times. In such a continual setup, we aim for learning with zero forgetting, \IE, with no degradation in synthesis quality over earlier domains due to catastrophic forgetting. To this end, we introduce a novel framework that not only (i) enables seamless knowledge transfer in continual training but also (ii) guarantees zero forgetting with a small overhead cost. While being more memory efficient, thanks to continual learning, our model obtains better synthesis quality as compared against the brute-force solution that trains one full model for each domain. Especially, under extreme low-data regimes, our approach outperforms the brute-force one by a large margin.


title: The local PNG bias of neutral Hydrogen, ${\rm H_I}$

abstract: We use separate universe simulations with the IllustrisTNG galaxy formation model to predict the local PNG bias parameters $b_\phi$ and $b_{\phi\delta}$ of atomic neutral hydrogen, ${\rm H_I}$. These parameters and their relation to the linear density bias parameter $b_1$ play a key role in observational constraints of the local PNG parameter $f_{\rm NL}$ using the ${\rm H_I}$ power spectrum and bispectrum. Our results show that the popular calculation based on the universality of the halo mass function overpredicts the $b_\phi(b_1)$ and $b_{\phi\delta}(b_1)$ relations measured in the simulations. In particular, our results show that at $z \lesssim 1$ the ${\rm H_I}$ power spectrum is more sensitive to $f_{\rm NL}$ compared to previously thought ($b_\phi$ is more negative), but is less sensitive at other epochs ($b_\phi$ is less positive). We discuss how this can be explained by the competition of physical effects such as that large-scale gravitational potentials with local PNG (i) accelerate the conversion of hydrogen to heavy elements by star formation, (ii) enhance the effects of baryonic feedback that eject the gas to regions more exposed to ionizing radiation, and (iii) promote the formation of denser structures that shield the ${\rm H_I}$ more efficiently. Our numerical results can be used to revise existing forecast studies on $f_{\rm NL}$ using 21cm line-intensity mapping data. Despite this first step towards predictions for the local PNG bias parameters of ${\rm H_I}$, we emphasize that more work is needed to assess their sensitivity on the assumed galaxy formation physics and ${\rm H_I}$ modeling strategy.


title: Human Parity on CommonsenseQA: Augmenting Self-Attention with External Attention

abstract: Most of today's AI systems focus on using self-attention mechanisms and transformer architectures on large amounts of diverse data to achieve impressive performance gains. In this paper, we propose to augment the transformer architecture with an external attention mechanism to bring external knowledge and context to bear. By integrating external information into the prediction process, we hope to reduce the need for ever-larger models and increase the democratization of AI systems. We find that the proposed external attention mechanism can significantly improve the performance of existing AI systems, allowing practitioners to easily customize foundation AI models to many diverse downstream applications. In particular, we focus on the task of Commonsense Reasoning, demonstrating that the proposed external attention mechanism can augment existing transformer models and significantly improve the model's reasoning capabilities. The proposed system, Knowledgeable External Attention for commonsense Reasoning (KEAR), reaches human parity on the open CommonsenseQA research benchmark with an accuracy of 89.4\% in comparison to the human accuracy of 88.9\%.


title: NP-completeness of the Active Time Scheduling Problem

abstract: In this paper, we study the active time scheduling problem. We are given n jobs with integral processing times each of which has an integral release time and deadline. The goal is to schedule all the jobs on a machine that can work on b jobs simultaneously, and the objective is to minimize the number of time slots for which the machine is active. The active time scheduling model was introduced by Chang et al. in the context of energy-efficient scheduling. Surprisingly, despite the development of a number of constant factor approximation algorithms for the problem, the complexity of this fundamental problem had remained open. In this paper, we resolve this open problem and show that the active time scheduling problem is indeed NP-complete.


title: Impact of Target Word and Context on End-to-End Metonymy Detection

abstract: Metonymy is a figure of speech in which an entity is referred to by another related entity. The task of metonymy detection aims to distinguish metonymic tokens from literal ones. Until now, metonymy detection methods attempt to disambiguate only a single noun phrase in a sentence, typically location names or organization names. In this paper, we disambiguate every word in a sentence by reformulating metonymy detection as a sequence labeling task. We also investigate the impact of target word and context on metonymy detection. We show that the target word is less useful for detecting metonymy in our dataset. On the other hand, the entity types that are associated with domain-specific words in their context are easier to solve. This shows that the context words are much more relevant for detecting metonymy.


title: Functional Regularization for Reinforcement Learning via Learned Fourier Features

abstract: We propose a simple architecture for deep reinforcement learning by embedding inputs into a learned Fourier basis and show that it improves the sample efficiency of both state-based and image-based RL. We perform infinite-width analysis of our architecture using the Neural Tangent Kernel and theoretically show that tuning the initial variance of the Fourier basis is equivalent to functional regularization of the learned deep network. That is, these learned Fourier features allow for adjusting the degree to which networks underfit or overfit different frequencies in the training data, and hence provide a controlled mechanism to improve the stability and performance of RL optimization. Empirically, this allows us to prioritize learning low-frequency functions and speed up learning by reducing networks' susceptibility to noise in the optimization process, such as during Bellman updates. Experiments on standard state-based and image-based RL benchmarks show clear benefits of our architecture over the baselines. Website at https://alexanderli.com/learned-fourier-features


title: DoodleFormer: Creative Sketch Drawing with Transformers

abstract: Creative sketching or doodling is an expressive activity, where imaginative and previously unseen depictions of everyday visual objects are drawn. Creative sketch image generation is a challenging vision problem, where the task is to generate diverse, yet realistic creative sketches possessing the unseen composition of the visual-world objects. Here, we propose a novel coarse-to-fine two-stage framework, DoodleFormer, that decomposes the creative sketch generation problem into the creation of coarse sketch composition followed by the incorporation of fine-details in the sketch. We introduce graph-aware transformer encoders that effectively capture global dynamic as well as local static structural relations among different body parts. To ensure diversity of the generated creative sketches, we introduce a probabilistic coarse sketch decoder that explicitly models the variations of each sketch body part to be drawn. Experiments are performed on two creative sketch datasets: Creative Birds and Creative Creatures. Our qualitative, quantitative and human-based evaluations show that DoodleFormer outperforms the state-of-the-art on both datasets, yielding realistic and diverse creative sketches. On Creative Creatures, DoodleFormer achieves an absolute gain of 25 in terms of Fr`echet inception distance (FID) over the state-of-the-art. We also demonstrate the effectiveness of DoodleFormer for related applications of text to creative sketch generation and sketch completion.


title: Novel Local Radiomic Bayesian Classifiers for Non-Invasive Prediction of MGMT Methylation Status in Glioblastoma

abstract: Glioblastoma, an aggressive brain cancer, is amongst the most lethal of all cancers. Expression of the O6-methylguanine-DNA-methyltransferase (MGMT) gene in glioblastoma tumor tissue is of clinical importance as it has a significant effect on the efficacy of Temozolomide, the primary chemotherapy treatment administered to glioblastoma patients. Currently, MGMT methylation is determined through an invasive brain biopsy and subsequent genetic analysis of the extracted tumor tissue. In this work, we present novel Bayesian classifiers that make probabilistic predictions of MGMT methylation status based on radiomic features extracted from FLAIR-sequence magnetic resonance imagery (MRIs). We implement local radiomic techniques to produce radiomic activation maps and analyze MRIs for the MGMT biomarker based on statistical features of raw voxel-intensities. We demonstrate the ability for simple Bayesian classifiers to provide a boost in predictive performance when modelling local radiomic data rather than global features. The presented techniques provide a non-invasive MRI-based approach to determining MGMT methylation status in glioblastoma patients.


title: Estimating nonlinear stability from time series data

abstract: Basin stability (BS) is a measure of nonlinear stability in multi-stable dynamical systems. BS has previously been estimated using Monte-Carlo simulations, which requires the explicit knowledge of a dynamical model. We discuss the requirements for estimating BS from time series data in the presence of strong perturbations, and illustrate our approach for two simple models of climate tipping elements: the Amazon rain forest and the thermohaline ocean circulation. We discuss the applicability of our method to observational data as constrained by the relevant time scales of total observation time, typical return time of perturbations and internal convergence time scale of the system of interest and other factors.

categories: physics.data-an nlin.CD physics.ao-ph physics.geo-ph


title: Modeling Demand Flexibility of RES-based Virtual Power Plants

abstract: In this paper, an approach to evaluate the benefits of demand flexibility for Virtual Power Plants (VPPs) is presented. The flexible demands chosen in this study are part of a renewable energy source-based VPP that participates in Day-Ahead Market (DAM) and Intra-Day Market (IDM) and has dispatchable and non-dispatchable assets. A demand model with bi-level flexibility is proposed: the first level is associated with DAM, whereas the second level is related to IDM sessions. Simulations are carried out considering a 12-node network to ascertain the eventual impacts of modeling demand flexibility on VPP operation. The market structure considered in the case study resembles the different trading floors in the Spanish electricity market. Results obtained show that the proposed demand flexibility scheme increases the overall profit of the VPP, as well as the revenues of the demand owners without disrupting consumers' comfort.


title: Multi-scale Graph Convolutional Networks with Self-Attention

abstract: Graph convolutional networks (GCNs) have achieved remarkable learning ability for dealing with various graph structural data recently. In general, deep GCNs do not work well since graph convolution in conventional GCNs is a special form of Laplacian smoothing, which makes the representation of different nodes indistinguishable. In the literature, multi-scale information was employed in GCNs to enhance the expressive power of GCNs. However, over-smoothing phenomenon as a crucial issue of GCNs remains to be solved and investigated. In this paper, we propose two novel multi-scale GCN frameworks by incorporating self-attention mechanism and multi-scale information into the design of GCNs. Our methods greatly improve the computational efficiency and prediction accuracy of the GCNs model. Extensive experiments on both node classification and graph classification demonstrate the effectiveness over several state-of-the-art GCNs. Notably, the proposed two architectures can efficiently mitigate the over-smoothing problem of GCNs, and the layer of our model can even be increased to $64$.


title: Effect of dust-ion collision on superthermal plasmas in cylindrical and spherical geometry

abstract: Dust-ion collisional effect in dust ion acoustic waves is investigated in the framework of spherical and cylindrical geometry. For such study the fluid equation are takes and collisional effect is consider in the equation of momentum balance. Using reductive perturbation technique the evolution equation is derived with additional damping term for non planner geometry. Using the conservation principle an analytical time dependent solitary wave solution is obtained for the evolution equation. It is seen that dust ion collisional frequency has an significant effect on the width and amplitude of the solitary wave solution.


title: Computational study of non-unitary partitions

abstract: Following Cayley, MacMahon, and Sylvester, define a non-unitary partition to be an integer partition with no part equal to one, and let $\nu(n)$ denote the number of non-unitary partitions of size $n$. In a 2021 paper, the sixth author proved a formula to compute $p(n)$ by enumerating only non-unitary partitions of size $n$, and recorded a number of conjectures regarding the growth of $\nu(n)$ as $n\to \infty$. Here we refine and prove some of these conjectures. For example, we prove $p(n) \sim \nu(n)\sqrt{n/\zeta(2)}$ as $n\to \infty$, and give Ramanujan-like congruences between $p(n)$ and $\nu(n)$ such as $p(5n)\equiv \nu(5n)\ (\operatorname{mod} 5)$.


title: A Deep-Learning Intelligent System Incorporating Data Augmentation for Short-Term Voltage Stability Assessment of Power Systems

abstract: Facing the difficulty of expensive and trivial data collection and annotation, how to make a deep learning-based short-term voltage stability assessment (STVSA) model work well on a small training dataset is a challenging and urgent problem. Although a big enough dataset can be directly generated by contingency simulation, this data generation process is usually cumbersome and inefficient; while data augmentation provides a low-cost and efficient way to artificially inflate the representative and diversified training datasets with label preserving transformations. In this respect, this paper proposes a novel deep-learning intelligent system incorporating data augmentation for STVSA of power systems. First, due to the unavailability of reliable quantitative criteria to judge the stability status for a specific power system, semi-supervised cluster learning is leveraged to obtain labeled samples in an original small dataset. Second, to make deep learning applicable to the small dataset, conditional least squares generative adversarial networks (LSGAN)-based data augmentation is introduced to expand the original dataset via artificially creating additional valid samples. Third, to extract temporal dependencies from the post-disturbance dynamic trajectories of a system, a bi-directional gated recurrent unit with attention mechanism based assessment model is established, which bi-directionally learns the significant time dependencies and automatically allocates attention weights. The test results demonstrate the presented approach manages to achieve better accuracy and a faster response time with original small datasets. Besides classification accuracy, this work employs statistical measures to comprehensively examine the performance of the proposal.


title: Contrastive Cycle Adversarial Autoencoders for Single-cell Multi-omics Alignment and Integration

abstract: Muilti-modality data are ubiquitous in biology, especially that we have entered the multi-omics era, when we can measure the same biological object (cell) from different aspects (omics) to provide a more comprehensive insight into the cellular system. When dealing with such multi-omics data, the first step is to determine the correspondence among different modalities. In other words, we should match data from different spaces corresponding to the same object. This problem is particularly challenging in the single-cell multi-omics scenario because such data are very sparse with extremely high dimensions. Secondly, matched single-cell multi-omics data are rare and hard to collect. Furthermore, due to the limitations of the experimental environment, the data are usually highly noisy. To promote the single-cell multi-omics research, we overcome the above challenges, proposing a novel framework to align and integrate single-cell RNA-seq data and single-cell ATAC-seq data. Our approach can efficiently map the above data with high sparsity and noise from different spaces to a low-dimensional manifold in a unified space, making the downstream alignment and integration straightforward. Compared with the other state-of-the-art methods, our method performs better in both simulated and real single-cell data. The proposed method is helpful for the single-cell multi-omics research. The improvement for integration on the simulated data is significant.


title: Communication and Energy Efficient Slimmable Federated Learning via Superposition Coding and Successive Decoding

abstract: Mobile devices are indispensable sources of big data. Federated learning (FL) has a great potential in exploiting these private data by exchanging locally trained models instead of their raw data. However, mobile devices are often energy limited and wirelessly connected, and FL cannot cope flexibly with their heterogeneous and time-varying energy capacity and communication throughput, limiting the adoption. Motivated by these issues, we propose a novel energy and communication efficient FL framework, coined SlimFL. To resolve the heterogeneous energy capacity problem, each device in SlimFL runs a width-adjustable slimmable neural network (SNN). To address the heterogeneous communication throughput problem, each full-width (1.0x) SNN model and its half-width ($0.5$x) model are superposition-coded before transmission, and successively decoded after reception as the 0.5x or $1.0$x model depending on the channel quality. Simulation results show that SlimFL can simultaneously train both $0.5$x and $1.0$x models with reasonable accuracy and convergence speed, compared to its vanilla FL counterpart separately training the two models using $2$x more communication resources. Surprisingly, SlimFL achieves even higher accuracy with lower energy footprints than vanilla FL for poor channels and non-IID data distributions, under which vanilla FL converges slowly.


title: Synthetic ECG Signal Generation Using Generative Neural Networks

abstract: Electrocardiogram (ECG) datasets tend to be highly imbalanced due to the scarcity of abnormal cases. Additionally, the use of real patients' ECGs is highly regulated due to privacy issues. Therefore, there is always a need for more ECG data, especially for the training of automatic diagnosis machine learning models, which perform better when trained on a balanced dataset. We studied the synthetic ECG generation capability of 5 different models from the generative adversarial network (GAN) family and compared their performances, the focus being only on Normal cardiac cycles. Dynamic Time Warping (DTW), Fr\'echet, and Euclidean distance functions were employed to quantitatively measure performance. Five different methods for evaluating generated beats were proposed and applied. We also proposed 3 new concepts (threshold, accepted beat and productivity rate) and employed them along with the aforementioned methods as a systematic way for comparison between models. The results show that all the tested models can, to an extent, successfully mass-generate acceptable heartbeats with high similarity in morphological features, and potentially all of them can be used to augment imbalanced datasets. However, visual inspections of generated beats favors BiLSTM-DC GAN and WGAN, as they produce statistically more acceptable beats. Also, with regards to productivity rate, the Classic GAN is superior with a 72% productivity rate. We also designed a simple experiment with the state-of-the-art classifier (ECGResNet34) to show empirically that the augmentation of the imbalanced dataset by synthetic ECG signals could improve the performance of classification significantly.


title: DIY Graphics Tab: A Cost-Effective Alternative to Graphics Tablet for Educators

abstract: Everyday, more and more people are turning to online learning, which has altered our traditional classroom method. Recording lectures has always been a normal task for online educators, and it has lately become even more important during the epidemic because actual lessons are still being postponed in several countries. When recording lectures, a graphics tablet is a great substitute for a whiteboard because of its portability and ability to interface with computers. This graphic tablet, however, is too expensive for the majority of instructors. In this paper, we propose a computer vision-based alternative to the graphics tablet for instructors and educators, which functions largely in the same way as a graphic tablet but just requires a pen, paper, and a laptop's webcam. We call it "Do-It-Yourself Graphics Tab" or "DIY Graphics Tab". Our system receives a sequence of images of a person's writing on paper acquired by a camera as input and outputs the screen containing the contents of the writing from the paper. The task is not straightforward since there are many obstacles such as occlusion due to the person's hand, random movement of the paper, poor lighting condition, perspective distortion due to the angle of view, etc. A pipeline is used to route the input recording through our system, which conducts instance segmentation and preprocessing before generating the appropriate output. We also conducted user experience evaluations from the teachers and students, and their responses are examined in this paper.


title: Toward a Taxonomy of Trust for Probabilistic Machine Learning

abstract: Probabilistic machine learning increasingly informs critical decisions in medicine, economics, politics, and beyond. We need evidence to support that the resulting decisions are well-founded. To aid development of trust in these decisions, we develop a taxonomy delineating where trust in an analysis can break down: (1) in the translation of real-world goals to goals on a particular set of available training data, (2) in the translation of abstract goals on the training data to a concrete mathematical problem, (3) in the use of an algorithm to solve the stated mathematical problem, and (4) in the use of a particular code implementation of the chosen algorithm. We detail how trust can fail at each step and illustrate our taxonomy with two case studies: an analysis of the efficacy of microcredit and The Economist's predictions of the 2020 US presidential election. Finally, we describe a wide variety of methods that can be used to increase trust at each step of our taxonomy. The use of our taxonomy highlights steps where existing research work on trust tends to concentrate and also steps where establishing trust is particularly challenging.


title: Adapting BERT for Continual Learning of a Sequence of Aspect Sentiment Classification Tasks

abstract: This paper studies continual learning (CL) of a sequence of aspect sentiment classification (ASC) tasks. Although some CL techniques have been proposed for document sentiment classification, we are not aware of any CL work on ASC. A CL system that incrementally learns a sequence of ASC tasks should address the following two issues: (1) transfer knowledge learned from previous tasks to the new task to help it learn a better model, and (2) maintain the performance of the models for previous tasks so that they are not forgotten. This paper proposes a novel capsule network based model called B-CL to address these issues. B-CL markedly improves the ASC performance on both the new task and the old tasks via forward and backward knowledge transfer. The effectiveness of B-CL is demonstrated through extensive experiments.


title: Clarification of Basic Concepts for Electromagnetic Interference Shielding Effectiveness

abstract: There exists serious miscomprehension in the open literature about the electromagnetic interference shielding effectiveness (EMI SE) as a critical index to evaluate the shielding performance, which is misleading to the graduates and newcomers embarking on the field of electromagnetic shielding materials. EMI SE is defined as the sum of three terms including reflection loss, absorption loss and multiple reflection loss in the classical Schelkunoff theory, while it is decomposed into two terms named reflection loss and absorption loss in practice, which is called Calculation theory here. In this paper, we elucidate the widely-seen misconceptions connected with EMI SE via theoretical derivation and instance analysis. Firstly, the terms in Calculation theory are often mistakenly regarded as the approximation of the terms with the same names in Schelkunoff theory when multiple reflection loss is negligible. Secondly, it is insufficient and unreasonable to determine the absorption-dominant shielding performance in the case that absorption loss is higher than reflection loss since reflection loss and absorption loss cannot represent the actual levels of reflected and absorbed power. Power coefficients are recommended to compare the contribution of reflection and absorption to shielding performance. Thirdly, multiple reflection effect is included in the definitions of reflection loss and absorption loss in Calculation theory, and the effect of multiple reflections on shielding property is clarified as against the commonly wrong understandings. These clarifications offer correct comprehension about the shielding mechanism and assessment of reflection and absorption contribution to the total shielding.


title: Dynamic Graph Learning-Neural Network for Multivariate Time Series Modeling

abstract: Multivariate time series forecasting is a challenging task because the data involves a mixture of long- and short-term patterns, with dynamic spatio-temporal dependencies among variables. Existing graph neural networks (GNN) typically model multivariate relationships with a pre-defined spatial graph or learned fixed adjacency graph. It limits the application of GNN and fails to handle the above challenges. In this paper, we propose a novel framework, namely static- and dynamic-graph learning-neural network (SDGL). The model acquires static and dynamic graph matrices from data to model long- and short-term patterns respectively. Static matric is developed to capture the fixed long-term association pattern via node embeddings, and we leverage graph regularity for controlling the quality of the learned static graph. To capture dynamic dependencies among variables, we propose dynamic graphs learning method to generate time-varying matrices based on changing node features and static node embeddings. And in the method, we integrate the learned static graph information as inductive bias to construct dynamic graphs and local spatio-temporal patterns better. Extensive experiments are conducted on two traffic datasets with extra structural information and four time series datasets, which show that our approach achieves state-of-the-art performance on almost all datasets. If the paper is accepted, I will open the source code on github.


title: Fast verification of spent nuclear fuel dry casks using cosmic ray muons: Monte Carlo simulation study

abstract: An application of muon scattering tomography is investigated for the non-destructive verification of spent nuclear fuel dry storage casks in context of international treaty declarations in proving that the diversion of plutonium has not occurred during long term storage in the spent fuel storage facility. It was shown in the previous studies that the method of muon scattering tomography can provide a detailed 3-D image of the spent fuel container but with the measurement time required for that is in the range of a day or even longer. In this paper we investigate the feasibility of a fast verification of dry storage casks in order to obtain an immediate answer without a detailed 3-D image reconstruction of container content if at least one fuel assembly is missing. We have performed preliminary modelling using Geant4 of safeguards verification of spent fuel dry storage cask HI-STORM 190 with cosmic ray muons. The simulated data were analysed using the Point of Closest Approach method and Kolmogorov statistical test. To infer the state of content of dry casks we have performed the comparison and further classification of muon observable for the scenarios of one missing fuel assembly with respect to fully loaded casks. The results of this study indicate that one missing fuel assembly can be detected with measurement times of about one hour.