title: PM-LSH: a fast and accurate in-memory framework for high-dimensional approximate NN and closest pair search

abstract: Nearest neighbor (NN) search is inherently computationally expensive in high-dimensional spaces due to the curse of dimensionality. As a well-known solution, locality-sensitive hashing (LSH) is able to answer c-approximate NN (c-ANN) queries in sublinear time with constant probability. Existing LSH methods focus mainly on building hash bucket-based indexing such that the candidate points can be retrieved quickly. However, existing coarse-grained structures fail to offer accurate distance estimation for candidate points, which translates into additional computational overhead when having to examine unnecessary points. This in turn reduces the performance of query processing. In contrast, we propose a fast and accurate in-memory LSH framework, called PM-LSH, that aims to compute c-ANN queries on large-scale, high-dimensional datasets. First, we adopt a simple yet effective PM-tree to index the data points. Second, we develop a tunable confidence interval to achieve accurate distance estimation and guarantee high result quality. Third, we propose an efficient algorithm on top of the PM-tree to improve the performance of computing c-ANN queries. In addition, we extend PM-LSH to support closest pair (CP) search in high-dimensional spaces. We again adopt the PM-tree to organize the points in a lowdimensional space, and we propose a branch and bound algorithm together with a radius pruning technique to improve the performance of computing c-approximate closest pair (c-ACP) queries. Extensive experiments with real-world data offer evidence that PM-LSH is capable of outperforming existing proposals with respect to both efficiency and accuracy for both NN and CP search.


title: Rate-Exponent Region for a Class of Distributed Hypothesis Testing Against Conditional Independence Problems

abstract: We study a class of $K$-encoder hypothesis testing against conditional independence problems. Under the criterion that stipulates minimization of the Type II error subject to a (constant) upper bound $\epsilon$ on the Type I error, we characterize the set of encoding rates and exponent for both discrete memoryless and memoryless vector Gaussian settings. For the DM setting, we provide a converse proof and show that it is achieved using the Quantize-Bin-Test scheme of Rahman and Wagner. For the memoryless vector Gaussian setting, we develop a tight outer bound by means of a technique that relies on the de Bruijn identity and the properties of Fisher information. In particular, the result shows that for memoryless vector Gaussian sources the rate-exponent region is exhausted using the Quantize-Bin-Test scheme with \textit{Gaussian} test channels; and there is \textit{no} loss in performance caused by restricting the sensors' encoders not to employ time sharing. Furthermore, we also study a variant of the problem in which the source, not necessarily Gaussian, has finite differential entropy and the sensors' observations noises under the null hypothesis are Gaussian. For this model, our main result is an upper bound on the exponent-rate function. The bound is shown to mirror a corresponding explicit lower bound, except that the lower bound involves the source power (variance) whereas the upper bound has the source entropy power. Part of the utility of the established bound is for investigating asymptotic exponent/rates and losses incurred by distributed detection as function of the number of sensors.


title: Towards GaAs Thin-Film Tracking Detectors

abstract: Silicon-based tracking detectors have been used in several important applications, such as in cancer therapy using particle beams, and for the discovery of new elementary particles at the Large Hadron Collider at CERN. III-V semiconductor materials are an attractive alternative to silicon for this application, as they have some superior physical properties. They could meet the demands for fast timing detectors allowing time-of-flight measurements with ps resolution while being radiation tolerant and cost-efficient. As a material with a larger density, higher atomic number Z and much higher electron mobility than silicon, GaAs exhibits faster signal collection and a larger signal per {\mu}m of sensor thickness. In this work, we report on the fabrication of n-in-n GaAs thin-film devices intended to serve next-generation high-energy particle tracking detectors. Molecular beam epitaxy (MBE) was used to grow high-quality GaAs films with doping levels sufficiently low to achieve full depletion for detectors with an active thickness of 10 {\mu}m. The signal collection speed of the detector structures was assessed using the transient current technique (TCT). To elucidate the structural properties of the detector, Kelvin probe force microscopy (KPFM) was used, which confirmed the formation of the junction in the detector and revealed residual doping in the intrinsic layer. Our results suggest that GaAs thin films are suitable candidates to achieve thin and radiation-tolerant tracking detectors.


title: A new metaheuristic approach for the art gallery problem

abstract: In the problem "Localization and trilateration with the minimum number of landmarks", we faced the 3-Guard and classic Art Gallery Problems. The goal of the art gallery problem is to find the minimum number of guards within a simple polygon to observe and protect its entirety. It has many applications in robotics, telecommunications, etc. There are some approaches to handle the art gallery problem that is theoretically NP-hard. This paper offers an efficient method based on the Particle Filter algorithm which solves the most fundamental state of the problem in a nearly optimal manner. The experimental results on the random polygons generated by Bottino et al. \cite{bottino2011nearly} show that the new method is more accurate with fewer or equal guards. Furthermore, we discuss resampling and particle numbers to minimize the run time.


title: End-to-End Natural Language Understanding Pipeline for Bangla Conversational Agents

abstract: Chatbots are intelligent software built to be used as a replacement for human interaction. Existing studies typically do not provide enough support for low-resource languages like Bangla. Due to the increasing popularity of social media, we can also see the rise of interactions in Bangla transliteration (mostly in English) among the native Bangla speakers. In this paper, we propose a novel approach to build a Bangla chatbot aimed to be used as a business assistant which can communicate in low-resource languages like Bangla and Bangla Transliteration in English with high confidence consistently. Since annotated data was not available for this purpose, we had to work on the whole machine learning life cycle (data preparation, machine learning modeling, and model deployment) using Rasa Open Source Framework, fastText embeddings, Polyglot embeddings, Flask, and other systems as building blocks. While working with the skewed annotated dataset, we try out different components and pipelines to evaluate which works best and provide possible reasoning behind the observed results. Finally, we present a pipeline for intent classification and entity extraction which achieves reasonable performance (accuracy: 83.02%, precision: 80.82%, recall: 83.02%, F1-score: 80%).


title: Contractive inequalities for mixed norm spaces and the Beta function

abstract: For a wide range of pairs of mixed norm spaces such that one space is contained in another, we characterize all cases when contractive norm inequalities hold. In particular, this yields such results for many pairs of weighted Bergman spaces. Some inequalities of this type are motivated by their applications in Number Theory and in Mathematical Physics.


title: Conics meeting eight lines over perfect fields

abstract: Over the complex numbers, there are 92 plane conics meeting 8 general lines in projective 3-space. Using the Euler class and local degree from motivic homotopy theory, we give an enriched version of this result over any perfect field. This provides a weighted count of the number of plane conics meeting 8 general lines, where the weight of each conic is determined the geometry of its intersections with the 8 given lines. As a corollary, real conics meeting 8 general lines come in two families of equal size.


title: Meta-learning PINN loss functions

abstract: We propose a meta-learning technique for offline discovery of physics-informed neural network (PINN) loss functions. We extend earlier works on meta-learning, and develop a gradient-based meta-learning algorithm for addressing diverse task distributions based on parametrized partial differential equations (PDEs) that are solved with PINNs. Furthermore, based on new theory we identify two desirable properties of meta-learned losses in PINN problems, which we enforce by proposing a new regularization method or using a specific parametrization of the loss function. In the computational examples, the meta-learned losses are employed at test time for addressing regression and PDE task distributions. Our results indicate that significant performance improvement can be achieved by using a shared-among-tasks offline-learned loss function even for out-of-distribution meta-testing. In this case, we solve for test tasks that do not belong to the task distribution used in meta-training, and we also employ PINN architectures that are different from the PINN architecture used in meta-training. To better understand the capabilities and limitations of the proposed method, we consider various parametrizations of the loss function and describe different algorithm design options and how they may affect meta-learning performance.


title: Towards Better Laplacian Representation in Reinforcement Learning with Generalized Graph Drawing

abstract: The Laplacian representation recently gains increasing attention for reinforcement learning as it provides succinct and informative representation for states, by taking the eigenvectors of the Laplacian matrix of the state-transition graph as state embeddings. Such representation captures the geometry of the underlying state space and is beneficial to RL tasks such as option discovery and reward shaping. To approximate the Laplacian representation in large (or even continuous) state spaces, recent works propose to minimize a spectral graph drawing objective, which however has infinitely many global minimizers other than the eigenvectors. As a result, their learned Laplacian representation may differ from the ground truth. To solve this problem, we reformulate the graph drawing objective into a generalized form and derive a new learning objective, which is proved to have eigenvectors as its unique global minimizer. It enables learning high-quality Laplacian representations that faithfully approximate the ground truth. We validate this via comprehensive experiments on a set of gridworld and continuous control environments. Moreover, we show that our learned Laplacian representations lead to more exploratory options and better reward shaping.


title: Calliope -- A Polyphonic Music Transformer

abstract: The polyphonic nature of music makes the application of deep learning to music modelling a challenging task. On the other hand, the Transformer architecture seems to be a good fit for this kind of data. In this work, we present Calliope, a novel autoencoder model based on Transformers for the efficient modelling of multi-track sequences of polyphonic music. The experiments show that our model is able to improve the state of the art on musical sequence reconstruction and generation, with remarkably good results especially on long sequences.


title: Elekes-Szab\'o for groups, and approximate subgroups in weak general position

abstract: We show that with a suitable weak notion of general position, the Elekes-Szab\'o condition on the group operation of a connected complex algebraic group characterises nilpotence of the group. Along the way, we prove a Mordell-Lang result for generic finitely generated subgroups of commutative complex algebraic groups.


title: Multi-modality Deep Restoration of Extremely Compressed Face Videos

abstract: Arguably the most common and salient object in daily video communications is the talking head, as encountered in social media, virtual classrooms, teleconferences, news broadcasting, talk shows, etc. When communication bandwidth is limited by network congestions or cost effectiveness, compression artifacts in talking head videos are inevitable. The resulting video quality degradation is highly visible and objectionable due to high acuity of human visual system to faces. To solve this problem, we develop a multi-modality deep convolutional neural network method for restoring face videos that are aggressively compressed. The main innovation is a new DCNN architecture that incorporates known priors of multiple modalities: the video-synchronized speech signal and semantic elements of the compression code stream, including motion vectors, code partition map and quantization parameters. These priors strongly correlate with the latent video and hence they are able to enhance the capability of deep learning to remove compression artifacts. Ample empirical evidences are presented to validate the superior performance of the proposed DCNN method on face videos over the existing state-of-the-art methods.


title: BAMCAFE: A Bayesian Machine Learning Advanced Forecast Ensemble Method for Complex Turbulent Systems with Partial Observations

abstract: Ensemble forecast based on physics-informed models is one of the most widely used forecast algorithms for complex turbulent systems. A major difficulty in such a method is the model error that is ubiquitous in practice. Data-driven machine learning (ML) forecasts can reduce the model error but they often suffer from the partial and noisy observations. In this paper, a simple but effective Bayesian machine learning advanced forecast ensemble (BAMCAFE) method is developed, which combines an available imperfect physics-informed model with data assimilation (DA) to facilitate the ML ensemble forecast. In the BAMCAFE framework, a Bayesian ensemble DA is applied to create the training data of the ML model, which reduces the intrinsic error in the imperfect physics-informed model simulations and provides the training data of the unobserved variables. Then a generalized DA is employed for the initialization of the ML ensemble forecast. In addition to forecasting the optimal point-wise value, the BAMCAFE also provides an effective approach of quantifying the forecast uncertainty utilizing a non-Gaussian probability density function that characterizes the intermittency and extreme events. It is shown using a two-layer Lorenz 96 model that the BAMCAFE method can significant improve the forecasting skill compared to the typical reduced-order imperfect models with bare truncation or stochastic parameterization for both the observed and unobserved large-scale variables. It is also shown via a nonlinear conceptual model that the BAMCAFE leads to a comparable non-Gaussian forecast uncertainty as the perfect model while the associated imperfect physics-informed model suffers from large forecast biases.


title: Extending Text-to-Speech Synthesis with Articulatory Movement Prediction using Ultrasound Tongue Imaging

abstract: In this paper, we present our first experiments in text-to-articulation prediction, using ultrasound tongue image targets. We extend a traditional (vocoder-based) DNN-TTS framework with predicting PCA-compressed ultrasound images, of which the continuous tongue motion can be reconstructed in synchrony with synthesized speech. We use the data of eight speakers, train fully connected and recurrent neural networks, and show that FC-DNNs are more suitable for the prediction of sequential data than LSTMs, in case of limited training data. Objective experiments and visualized predictions show that the proposed solution is feasible and the generated ultrasound videos are close to natural tongue movement. Articulatory movement prediction from text input can be useful for audiovisual speech synthesis or computer-assisted pronunciation training.


title: On the Radon transform and the Dirac delta distribution in superspace

abstract: In this manuscript, we obtain a plane wave decomposition for the delta distribution in superspace, provided that the superdimension is not odd and negative. This decomposition allows for explicit inversion formulas for the super Radon transform in these cases. Moreover, we prove a more general Radon inversion formula valid for all possible integer values of the superdimension. The proof of this result comes along with the study of fractional powers of the super Laplacian, their fundamental solutions, and the plane wave decompositions of super Riesz kernels.


title: Ground State Cooling of an Ultracoherent Electromechanical System

abstract: Cavity electromechanics relies on parametric coupling between microwave and mechanical modes to manipulate the mechanical quantum state, and provide a coherent interface between different parts of hybrid quantum systems. High coherence of the mechanical mode is of key importance in such applications, in order to protect the quantum states it hosts from thermal decoherence. Here, we introduce an electromechanical system based around a soft-clamped mechanical resonator with an extremely high Q-factor (>$10^9$) held at very low (30 mK) temperatures. This ultracoherent mechanical resonator is capacitively coupled to a microwave mode, strong enough to enable ground-state-cooling of the mechanics ($\bar{n}_\mathrm{min}= 0.76\pm 0.16$). This paves the way towards exploiting the extremely long coherence times ($t_\mathrm{coh}>100 ms) offered by such systems for quantum information processing and state conversion.


title: Self-consistent dynamical maps for open quantum systems

abstract: In several cases, open quantum systems can be successfully described using master equations relying on Born-Markov approximations, but going beyond these approaches has become often necessary. In this work, we introduce the NCA and NCA-Markov dynamical maps for open quantum systems, which go beyond these master equations replacing the Born approximation with a self-consistent approximation, known as non-crossing approximation (NCA). These maps are formally similar to master equations, but allow to capture non-perturbative effects of the environment at a moderate extra numerical cost. To demonstrate their capabilities, we apply them to the spin-boson model at zero temperature for both a Ohmic and a sub-Ohmic environment, showing that they can both qualitatively capture its strong-coupling behaviour, and be quantitatively correct beyond standard master equations.


title: Quantile-Based Random Kaczmarz for corrupted linear systems of equations

abstract: We consider linear systems $Ax = b$ where $A \in \mathbb{R}^{m \times n}$ consists of normalized rows, $\|a_i\|_{\ell^2} = 1$, and where up to $\beta m$ entries of $b$ have been corrupted (possibly by arbitrarily large numbers). Haddock, Needell, Rebrova and Swartworth propose a quantile-based Random Kaczmarz method and show that for certain random matrices $A$ it converges with high likelihood to the true solution. We prove a deterministic version by constructing, for any matrix $A$, a number $\beta_A$ such that there is convergence for all perturbations with $\beta < \beta_A$. Assuming a random matrix heuristic, this proves convergence for tall Gaussian matrices with up to $\sim 0.5\%$ corruption (a number that can likely be improved).


title: A Thermodynamically-Consistent Phase Field Crystal Model of Solidification with Heat Flux

abstract: In this paper we describe a new model for solidification with heat flux using the phase field crystal (PFC) framework. The equations are thermodynamically consistent, in the sense that the time rate of change of the entropy density is positive in the bulk and at the boundaries of the domain of interest. The resulting model consists of two equations, a heat-like equation and a mass-conservation equation that describes how the atom density changes in time and space. The model is simple, yet it can properly capture the variation in the free energy landscape as the temperature is changed. We describe to construct a temperature-atom-density phase diagram using this energy landscape, and we give a simple demonstration of solidification using the model.


title: DebiasedDTA: A Framework for Improving the Generalizability of Drug-Target Affinity Prediction Models

abstract: Computational models that accurately predict the binding affinity of an input protein-chemical pair can accelerate drug discovery studies. These models are trained on available protein-chemical interaction datasets, which may contain dataset biases that may lead the model to learn dataset-specific patterns, instead of generalizable relationships. As a result, the prediction performance of models drops for previously unseen biomolecules, $\textit{i.e.}$ the prediction models cannot generalize to biomolecules outside of the dataset. The latest approaches that aim to improve model generalizability either have limited applicability or introduce the risk of degrading prediction performance. Here, we present DebiasedDTA, a novel drug-target affinity (DTA) prediction model training framework that addresses dataset biases to improve the generalizability of affinity prediction models. DebiasedDTA reweights the training samples to mitigate the effect of dataset biases and is applicable to most DTA prediction models. The results suggest that models trained in the DebiasedDTA framework can achieve improved generalizability in predicting the interactions of the previously unseen biomolecules, as well as performance improvements on those previously seen. Extensive experiments with different biomolecule representations, model architectures, and datasets demonstrate that DebiasedDTA can upgrade DTA prediction models irrespective of the biomolecule representation, model architecture, and training dataset. Last but not least, we release DebiasedDTA as an open-source python library to enable other researchers to debias their own predictors and/or develop their own debiasing methods. We believe that this python library will corroborate and foster research to develop more generalizable DTA prediction models.


title: Dynamical characterization of central sets along filter

abstract: Using the notions of Topological dynamics, H. Furstenberg defined central sets and proved the Central Sets Theorem. Later V. Bergelson and N. Hindman characterized central sets in terms of algebra of the Stone-\v{C}ech Compactification of discrete semigroup. They found that central sets are the members of the minimal idempotents of $\beta S$, the Stone-\v{C}ech Compactification of a semigroup $\left(S,\cdot\right)$. We know that any closed subsemigroup of $\beta S$ is generated by a filter. We call a set $A$ to be a $\mathcal{F}$-central set if it is a member of a minimal idempotent of a closed subsemigroup of $\beta S$, generated by the filter $\mathcal{F}$. In this article we will characterize the $\mathcal{F}$-central sets dynamically.


title: Research on Metro Service Quality Improvement Schemes Considering Feasibility

abstract: It is an important management task of metro agencies to formulate reasonable improvement schemes based on the result of service quality surveys. Considering scores, weights, and improvement feasibility of service quality attributes in a certain period, this paper integrates Decision Tree (DT) into Importance-Performance analysis (IPA) to build a DT-IPA model, which is used to determine the improvement priority of attributes, and to quantify the improvement degree. If-then rules extracted from the optimal decision tree and the improvement feasibility computed by analytic hierarchy process are two main items derived from the DT-IPA model. They are used to optimize the initial improvement priority of attributes determined by IPA and to quantify the degree of improvement of the adjusted attributes. Then, the overall service quality can reach a high score, realizing the operation goal. The effectiveness of the DT-IPA model was verified through an empirical study which was taken place in Changsha Metro, China. The proposed method can be a decision-making tool for metro agency managers to improve the quality of metro service.


title: Inference on Individual Treatment Effects in Nonseparable Triangular Models

abstract: In nonseparable triangular models with a binary endogenous treatment and a binary instrumental variable, Vuong and Xu (2017) established identification results for individual treatment effects (ITEs) under the rank invariance assumption. Using their approach, Feng, Vuong, and Xu (2019) proposed a uniformly consistent kernel estimator for the density of the ITE that utilizes estimated ITEs. In this paper, we establish the asymptotic normality of the density estimator of Feng, Vuong, and Xu (2019) and show that the ITE estimation errors have a non-negligible effect on the asymptotic distribution of the estimator. We propose asymptotically valid standard errors that account for ITEs estimation, as well as a bias correction. Furthermore, we develop uniform confidence bands for the density of the ITE using the jackknife multiplier or nonparametric bootstrap critical values.


title: Geometric pumping and dephasing at topological phase transition

abstract: A measure-preserving formalism (MPF) is constructed and applied to spin/band models, which yield observations about pumping. It occurs at topological phase transition (TPT), i.e., a $Z_2$-flip, suggesting that $Z_2$ can imply bulk effects. The model's asymptotic behavior is analytically solved via MPF. The pumping probability is geometric, fractional, and has a ceiling of $\frac{1}{2}$. Intriguingly, theorems are proved about occurrence conditions, which are linked to the system's dimension and the distinction between rational and irrational numbers. Experimental detection is discussed.


title: LATTE: LSTM Self-Attention based Anomaly Detection in Embedded Automotive Platforms

abstract: Modern vehicles can be thought of as complex distributed embedded systems that run a variety of automotive applications with real-time constraints. Recent advances in the automotive industry towards greater autonomy are driving vehicles to be increasingly connected with various external systems (e.g., roadside beacons, other vehicles), which makes emerging vehicles highly vulnerable to cyber-attacks. Additionally, the increased complexity of automotive applications and the in-vehicle networks results in poor attack visibility, which makes detecting such attacks particularly challenging in automotive systems. In this work, we present a novel anomaly detection framework called LATTE to detect cyber-attacks in Controller Area Network (CAN) based networks within automotive platforms. Our proposed LATTE framework uses a stacked Long Short Term Memory (LSTM) predictor network with novel attention mechanisms to learn the normal operating behavior at design time. Subsequently, a novel detection scheme (also trained at design time) is used to detect various cyber-attacks (as anomalies) at runtime. We evaluate our proposed LATTE framework under different automotive attack scenarios and present a detailed comparison with the best-known prior works in this area, to demonstrate the potential of our approach.


title: Effect of pseudo datasets for the classification-based engineering design

abstract: Machine learning classification techniques have been used widely to recognize the feasible design domain and discover hidden patterns in engineering design. An accurate classification model needs a large dataset; however, generating a large dataset is costly for complex simulation-based problems. After training by a small dataset, surrogate models can generate a large pseudo dataset efficiently. Errors, however, may be introduced by surrogate modeling. This paper investigates the mutual effect of a large pseudo dataset and surrogate modeling uncertainty. Four widely used methods, i.e., Naive Bayes classifier, support vector machine, random forest regression, and artificial neural network for classification, are studied on four benchmark problems. Kriging is used as the basic surrogate model method. The results show that a large pseudo dataset improves the classification accuracy, which depends on both design problems and classification algorithms. Except for the Naive Bayes, the other three methods are recommended for using pseudo data to improve classification performance. Also, a wind turbine design problem is used to illustrate the effect of the pseudo dataset on feasible subspace recognition. The large pseudo dataset improves the recognized subspace bound greatly, which can be reproduced by classification models well except for the Naive Bayes classifier. Under the uncertainty of surrogate modeling, the random forest presents high robustness to recognize the feasible design domain bound, while the artificial neural network demonstrates a high sensitivity to uncertainty with the recognized bound deteriorated.


title: Infrastructure-less Wireless Connectivity for Mobile Robotic Systems in Logistics: Why Bluetooth Mesh Networking is Important?

abstract: Mobile robots have disrupted the material handling industry which is witnessing radical changes. The requirement for enhanced automation across various industry segments often entails mobile robotic systems operating in logistics facilities with little/no infrastructure. In such environments, out-of-box low-cost robotic solutions are desirable. Wireless connectivity plays a crucial role in successful operation of such mobile robotic systems. A wireless mesh network of mobile robots is an attractive solution; however, a number of system-level challenges create unique and stringent service requirements. The focus of this paper is the role of Bluetooth mesh technology, which is the latest addition to the Internet-of-Things (IoT) connectivity landscape, in addressing the challenges of infrastructure-less connectivity for mobile robotic systems. It articulates the key system-level design challenges from communication, control, cooperation, coverage, security, and navigation/localization perspectives, and explores different capabilities of Bluetooth mesh technology for such challenges. It also provides performance insights through real-world experimental evaluation of Bluetooth mesh while investigating its differentiating features against competing solutions.


title: Enumerating conjugacy classes of graphical groups over finite fields

abstract: Each graph and choice of a commutative ring gives rise to an associated graphical group. In this article, we introduce and investigate graph polynomials that enumerate conjugacy classes of graphical groups over finite fields according to their sizes.


title: Multifractality in non-unitary random dynamics

abstract: We explore the multifractality of the steady state wave function in non-unitary random quantum dynamics in one dimension. We focus on two classes of random systems: the hybrid Clifford circuit model and the non-unitary free fermion dynamics. In the hybrid Clifford model, we map the measurement driven transition to an Anderson localization transition in an effective graph space by using properties of the stabilizer state. We show that the volume law phase with nonzero measurement rate is non-ergodic in the graph space and exhibits weak multifractal behavior. We apply the same method to the hybrid Clifford quantum automaton circuit and obtain similar multifractality in the volume law phase. For the non-unitary random free fermion system with a critical steady state, we compute the moments of the probability distribution of the single particle wave function and demonstrate that it is also weakly multifractal and has strong variations in real space.


title: ProGS: Property Graph Shapes Language (Extended Version)

abstract: Property graphs constitute data models for representing knowledge graphs. They allow for the convenient representation of facts, including facts about facts, represented by triples in subject or object position of other triples. Knowledge graphs such as Wikidata are created by a diversity of contributors and a range of sources leaving them prone to two types of errors. The first type of error, falsity of facts, is addressed by property graphs through the representation of provenance and validity, making triples occur as first-order objects in subject position of metadata triples. The second type of error, violation of domain constraints, has not been addressed with regard to property graphs so far. In RDF representations, this error can be addressed by shape languages such as SHACL or ShEx, which allow for checking whether graphs are valid with respect to a set of domain constraints. Borrowing ideas from the syntax and semantics definitions of SHACL, we design a shape language for property graphs, ProGS, which allows for formulating shape constraints on property graphs including their specific constructs, such as edges with identities and key-value annotations to both nodes and edges. We define a formal semantics of ProGS, investigate the resulting complexity of validating property graphs against sets of ProGS shapes, compare with corresponding results for SHACL, and implement a prototypical validator that utilizes answer set programming.


title: Strong recovery of geometric planted matchings

abstract: We study the problem of efficiently recovering the matching between an unlabelled collection of $n$ points in $\mathbb{R}^d$ and a small random perturbation of those points. We consider a model where the initial points are i.i.d. standard Gaussian vectors, perturbed by adding i.i.d. Gaussian vectors with variance $\sigma^2$. In this setting, the maximum likelihood estimator (MLE) can be found in polynomial time as the solution of a linear assignment problem. We establish thresholds on $\sigma^2$ for the MLE to perfectly recover the planted matching (making no errors) and to strongly recover the planted matching (making $o(n)$ errors) both for $d$ constant and $d = d(n)$ growing arbitrarily. Between these two thresholds, we show that the MLE makes $n^{\delta + o(1)}$ errors for an explicit $\delta \in (0, 1)$. These results extend to the geometric setting a recent line of work on recovering matchings planted in random graphs with independently-weighted edges. Our proof techniques rely on careful analysis of the combinatorial structure of partial matchings in large, weakly dependent random graphs using the first and second moment methods.


title: Asymmetric Multi-Stability from Relaxing the Rigid-Folding Conditions in a Stacked Miura-ori Cellular Solid

abstract: Traditionally, origami has been categorized into two groups according to their kinematics design: rigid and non-rigid origami. However, such categorization can be superficial, and rigid origami can obtain new mechanical properties by intentionally relaxing the rigid-folding kinematics. Based on numerical simulations using the bar-hinge approach and experiments, this study examines the multi-stability of a stacked Miura-origami cellular structure with different levels of facet compliance. The simulation and experiment results show that a unit cell in such cellular solid exhibits only two stable states if it follows the rigid origami kinematics; however, two more stable states are reachable if the origami facets become sufficiently compliant. Moreover, the switch between two certain stable states shows an asymmetric energy barrier, meaning that the unit cell follows fundamentally different deformation paths when it extends from one state to another compared to the opposite compression switch. As a result, the reaction force required for extending this unit cell between these two states can be higher than the compression switch. Such asymmetric multi-stability can be fine-tuned by tailoring the underlying origami design, and it can be extended into cellular solids with carefully placed voids. By showing the benefits of exploiting facet compliance, this study could foster multi-functional structures and material systems that traditional rigid origami cannot create.


title: A mathematical model to study the impact of intra-tumour heterogeneity on anti-tumour CD8+ T cell immune response

abstract: Intra-tumour heterogeneity (ITH) has a strong impact on the efficacy of the immune response against solid tumours. The number of sub-populations of cancer cells expressing different antigens and the percentage of immunogenic cells (i.e. tumour cells that are effectively targeted by immune cells) in a tumour are both expressions of ITH. Here, we present a spatially explicit stochastic individual-based model of the interaction dynamics between tumour cells and CD8+ T cells, which makes it possible to dissect out the specific impact of these two expressions of ITH on anti-tumour immune response. The set-up of numerical simulations of the model is defined so as to mimic scenarios considered in previous experimental studies. Moreover, the ability of the model to qualitatively reproduce experimental observations of successful and unsuccessful immune surveillance is demonstrated. First, the results of numerical simulations of this model indicate that the presence of a larger number of sub-populations of tumour cells that express different antigens is associated with a reduced ability of CD8+ T cells to mount an effective anti-tumour immune response. Second, the presence of a larger percentage of tumour cells that are not effectively targeted by CD8+ T cells may reduce the effectiveness of anti-tumour immunity. Ultimately, the mathematical model presented in this paper may provide a framework to help biologists and clinicians to better understand the mechanisms that are responsible for the emergence of different outcomes of immunotherapy.


title: On Mesh Deformation Techniques for Topology Optimization of Fluid-Structure Interaction Problems

abstract: Fluid-structure interactions are a widespread phenomenon in nature. Although their numerical modeling have come a long way, the application of numerical design tools to these multiphysics problems is still lagging behind. Gradient-based optimization is the most popular approach in topology optimization currently. Hence, it's a necessity to utilize mesh deformation techniques that have continuous, smooth derivatives. In this work, we address mesh deformation techniques for structured, quadrilateral meshes. We discuss and comment on two legacy mesh deformation techniques; namely the spring analogy model and the linear elasticity model. In addition, we propose a new technique based on the Yeoh hyperelasticity model. We focus on mesh quality as a gateway to mesh admissibility. We propose layered selective stiffening such that the elements adjacent to the fluid-structure interface - where the bulk of the mesh distortion occurs - are stiffened in consecutive layers. The legacy and the new models are able to sustain large deformations without deprecating the mesh quality, and the results are enhanced with using layered selective stiffening.


title: Connection between antennas, beam steering, and the moir\'{e} effect

abstract: The moir\'{e} effect provides an interpretation for the steering of antennas that form beams through internal spatial interferences. We make an explicit connection between such antennas and the moir\'{e} effect, and use it to model six planar antennas that steer by scaling, rotating, or translating operations. Three of the antennas illustrate how to use moir\'{e} patterns to generate antenna designs.


title: The Integer Sequence Transform $a \mapsto b$ where $b_n$ is the Number of Real Roots of the Polynomial $a_0 + a_1x + a_2x^2 + \cdots + a_nx^n$

abstract: We discuss the integer sequence transform $a \mapsto b$ where $b_n$ is the number of real roots of the polynomial $a_0 + a_1x + a_2x^2 + \cdots + a_nx^n$. It is shown that several sequences $a$ give the trivial sequence $b = (0,1,0,1, 0,1,\ldots)$, i.e., ${b_n = n \bmod 2}$, among them the Catalan numbers, central binomial coefficients, $n!$ and $\binom{n+k}{n}$ for a fixed $k$. We also look at some sequences $a$ for which $b$ is more interesting such as $a_n = (n+1)^k$ for $k \geq 3$. Further, general procedures are given for constructing real sequences $a_n$ for which $b_n$ is either always maximal or minimal.


title: Tropical cyclone intensity estimations over the Indian ocean using Machine Learning

abstract: Tropical cyclones are one of the most powerful and destructive natural phenomena on earth. Tropical storms and heavy rains can cause floods, which lead to human lives and economic loss. Devastating winds accompanying cyclones heavily affect not only the coastal regions, even distant areas. Our study focuses on the intensity estimation, particularly cyclone grade and maximum sustained surface wind speed (MSWS) of a tropical cyclone over the North Indian Ocean. We use various machine learning algorithms to estimate cyclone grade and MSWS. We have used the basin of origin, date, time, latitude, longitude, estimated central pressure, and pressure drop as attributes of our models. We use multi-class classification models for the categorical outcome variable, cyclone grade, and regression models for MSWS as it is a continuous variable. Using the best track data of 28 years over the North Indian Ocean, we estimate grade with an accuracy of 88% and MSWS with a root mean square error (RMSE) of 2.3. For higher grade categories (5-7), accuracy improves to an average of 98.84%. We tested our model with two recent tropical cyclones in the North Indian Ocean, Vayu and Fani. For grade, we obtained an accuracy of 93.22% and 95.23% respectively, while for MSWS, we obtained RMSE of 2.2 and 3.4 and $R^2$ of 0.99 and 0.99, respectively.


title: Non-binary Snow Index for Multi-Component Surfaces

abstract: A Non-Binary Snow Index for Multi-Component Surfaces (NBSI-MS) is proposed to map snow/ice cover. The NBSI-MS is based on the spectral characteristics of different Land Cover Types (LCTs) such as snow, water, vegetation, bare land, impervious, and shadow surfaces. This index can increase the separability between NBSI-MS values corresponding to snow from other LCTs and accurately delineate the snow/ice cover in non-binary maps. To test the robustness of the NBSI-MS, Greenland and France-Italy regions were examined where snow interacts with highly diversified geographical ecosystem. Data recorded by Landsat 5 TM, Landsat 8 OLI, and Sentinel-2A MSI satellites have been used. The NBSI-MS performance was also compared against the well-known NDSI, NDSII-1, S3, and SWI methods and evaluated based on Ground Reference Test Pixels (GRTPs) over non-binarized results. The results show that the NBSI-MS achieves overall accuracy (OA) ranging from 0.99 to 1 with kappa coefficient values in the same range as OA. The precision assessment confirms the performance superiority of the proposed NBSI-MS method for removing water and shadow surfaces over the compared relevant indices.


title: New Constraint of the Hubble Constant by Proper Motions of Radio Components Observed in AGN Twin-jets

abstract: As the advent of precision cosmology, the Hubble constant ($H_0$) inferred from the Lambda Cold Dark Matter fit to the Cosmic Microwave Background data is increasingly in tension with the measurements from the local distance ladder. To approach its real value, we need more independent methods to measure, or to make constraint of, the Hubble constant. In this paper, we apply a plain method, which is merely based on the Friedman-Lema\^itre-Robertson-Walker cosmology together with geometrical relations, to constrain the Hubble constant by proper motions of radio components observed in AGN twin-jets. Under the assumption that the ultimate ejection strengths in both sides of the twin-jet concerned are intrinsically the same, we obtain a lower limit of the $H_{\rm 0,min}=51.5\pm2.3\,\rm km\,s^{-1}\,Mpc^{-1}$ from the measured maximum proper motions of the radio components observed in the twin-jet of NGC 1052.


title: Real-Space Inversion and Super-Resolution of Ultrafast Scattering

abstract: Ultrafast scattering using X-rays or electrons is an emerging method to obtain structure dynamics at the atomic length and time scales. However, directly resolving in real-space atomic motions is inherently limited by the finite detector range and the probe energy. As a result, the time-resolved signal interpretation is mostly done in reciprocal space and relies on modeling and simulations of specific structures and processes. Here, we introduce a model-free approach to directly resolve scattering signals in real space, surpassing the diffraction limit, using scattering kernels and signal priors that naturally arise from the measurement constraints. We demonstrate the approach on simulated and experimental data, recover multiple atomic motions at sub-$\angstrom$ngstrom resolutions, and discuss the recovery accuracy and resolution limits vs signal fidelity. The approach offers a robust path to obtain high-resolution real-space information of atomic-scale structure dynamics using current time-resolved X-ray or electron scattering sources.


title: On Divergence- and Gradient-Preserving Coarse-Graining for Finite Volume Primitive Equation Ocean Models

abstract: We consider the problem of coarse-graining in the context of finite-volume fluid models. If a variable is defined on a high-resolution grid it may be coarse-grained so that it is defined on a grid of lower resolution. In general this will cause some information about the variable to be lost. In particular, horizontal divergences, gradients or other operators calculated on the coarse grid after projecting may differ from those calculated on the fine grid. In some cases we are able to choose averaging weights for coarse-graining such that the coarse-grid operators will give a result approximating that of the corresponding fine-grid operators applied on the fine grid. In this work we derive general conditions on the averaging weights that allow the divergence and gradient to be preserved. These conditions are applied to a regular triangular mesh with B-grid variable placement in which the fine-grid resolution is some integer multiple $N$ of the coarse-grid resolution. For this case we find particular values for the averaging weights that preserve the divergence, and a set of different averaging weights that preserve the gradient. We observe that the vertical component of the curl is also preserved by the same coarse-graining that preserves divergence. These coarse-grainings are applied to data from FESOM2 simulations and we demonstrate that using the particular coarse-grainings derived herein gives an overall reduction in $L^1$ error when compared with other methods.

categories: physics.ao-ph cs.NA math.NA physics.flu-dyn


title: Impact of Scene-Specific Enhancement Spectra on Matched Filter Greenhouse Gas Retrievals from Imaging Spectroscopy

abstract: Matched filter (MF) techniques have been widely used for retrieval of greenhouse gas enhancements (enh.) from imaging spectroscopy datasets. While multiple algorithmic techniques and refinements have been proposed, the greenhouse gas target spectrum used for concentration enh. estimation has remained largely unaltered since the introduction of quantitative MF retrievals. The magnitude of retrieved methane and carbon dioxide enh., and thereby integrated mass enh. (IME) and estimated flux of point-source emitters, is heavily dependent on this target spectrum. Current standard use of molecular absorption coefficients to create unit enh. target spectra does not account for absorption by background concentrations of greenhouse gases, solar and sensor geometry, or atmospheric water vapor absorption. We introduce geometric and atmospheric parameters into the generation of scene-specific (SS) unit enh. spectra to provide target spectra that are compatible with all greenhouse gas retrieval MF techniques. For methane plumes, IME resulting from use of standard, generic enh. spectra varied from -22 to +28.7% compared to SS enh. spectra. Due to differences in spectral shape between the generic and SS enh. spectra, differences in methane plume IME were linked to surface spectral characteristics in addition to geometric and atmospheric parameters. IME differences for carbon dioxide plumes, with generic enh. spectra producing integrated mass enh. -76.1 to -48.1% compared to SS enh. spectra. Fluxes calculated from these integrated enh. would vary by the same %s, assuming equivalent wind conditions. Methane and carbon dioxide IME were most sensitive to changes in solar zenith angle and ground elevation. SS target spectra can improve confidence in greenhouse gas retrievals and flux estimates across collections of scenes with diverse geometric and atmospheric conditions.

categories: physics.ao-ph astro-ph.EP astro-ph.IM eess.IV stat.AP


title: Unifying incidence and prevalence under a time-varying general branching process

abstract: Renewal equations are a popular approach used in modelling the number of new infections, i.e., incidence, in an outbreak. We develop a stochastic model of an outbreak based on a time-varying variant of the Crump-Mode-Jagers branching process. This model accommodates a time-varying reproduction number and a time-varying distribution for the generation interval. We then derive renewal-like integral equations for incidence, cumulative incidence and prevalence under this model. We show that the equations for incidence and prevalence are consistent with the so-called back-calculation relationship. We analyse two particular cases of these integral equations, one that arises from a Bellman-Harris process and one that arises from an inhomogeneous Poisson process model of transmission. We also show that the incidence integral equations that arise from both of these specific models agree with the renewal equation used ubiquitously in infectious disease modelling. We present a numerical discretisation scheme to solve these equations, and use this scheme to estimate rates of transmission from serological prevalence of SARS-CoV-2 in the UK and historical incidence data on Influenza, Measles, SARS and Smallpox.


title: Equivalent Laplacian and Adjacency Quantum Walks on Irregular Graphs

abstract: The continuous-time quantum walk is a particle evolving by Schr\"odinger's equation in discrete space. Encoding the space as a graph of vertices and edges, the Hamiltonian is proportional to the discrete Laplacian. In some physical systems, however, the Hamiltonian is proportional to the adjacency matrix instead. It is well-known that these quantum walks are equivalent when the graph is regular, i.e., when each vertex has the same number of neighbors. If the graph is irregular, however, the quantum walks evolve differently. In this paper, we show that for some irregular graphs, if the particle is initially localized at a certain vertex, the probability distributions of the two quantum walks are identical, even though the amplitudes differ. We analytically prove this for a graph with five vertices and a graph with six vertices. By simulating the walks on all 1,018,689,568 simple, connected, irregular graphs with eleven vertices or less, we found sixty-four graphs with this notion of equivalence. We also give eight infinite families of graphs supporting these equivalent walks.


title: 21cmVAE: A Very Accurate Emulator of the 21-cm Global Signal

abstract: Considerable observational efforts are being dedicated to measuring the sky-averaged (global) 21-cm signal of neutral hydrogen from Cosmic Dawn and the Epoch of Reionization. Deriving observational constraints on the astrophysics of this era requires modeling tools that can quickly and accurately generate theoretical signals across the wide astrophysical parameter space. For this purpose artificial neural networks were used to create the only two existing global signal emulators, 21cmGEM and globalemu. In this paper we introduce 21cmVAE, a neural network-based global signal emulator, trained on the same dataset of ~30,000 global signals as the other two emulators, but with a more direct prediction algorithm that prioritizes accuracy and simplicity. Using neural networks, we compute derivatives of the signals with respect to the astrophysical parameters and establish the most important astrophysical processes that drive the global 21-cm signal at different epochs. 21cmVAE has a relative rms error of only 0.34 - equivalently 0.54 mK - on average, which is a significant improvement compared to the existing emulators, and a run time of 0.04 seconds per parameter set. The emulator, the code, and the processed datasets are publicly available at https://github.com/christianhbye/21cmVAE and through https://zenodo.org/record/5904939.


title: Forster Decomposition and Learning Halfspaces with Noise

abstract: A Forster transform is an operation that turns a distribution into one with good anti-concentration properties. While a Forster transform does not always exist, we show that any distribution can be efficiently decomposed as a disjoint mixture of few distributions for which a Forster transform exists and can be computed efficiently. As the main application of this result, we obtain the first polynomial-time algorithm for distribution-independent PAC learning of halfspaces in the Massart noise model with strongly polynomial sample complexity, i.e., independent of the bit complexity of the examples. Previous algorithms for this learning problem incurred sample complexity scaling polynomially with the bit complexity, even though such a dependence is not information-theoretically necessary.


title: Few-shot Learning with Global Relatedness Decoupled-Distillation

abstract: Despite the success that metric learning based approaches have achieved in few-shot learning, recent works reveal the ineffectiveness of their episodic training mode. In this paper, we point out two potential reasons for this problem: 1) the random episodic labels can only provide limited supervision information, while the relatedness information between the query and support samples is not fully exploited; 2) the meta-learner is usually constrained by the limited contextual information of the local episode. To overcome these problems, we propose a new Global Relatedness Decoupled-Distillation (GRDD) method using the global category knowledge and the Relatedness Decoupled-Distillation (RDD) strategy. Our GRDD learns new visual concepts quickly by imitating the habit of humans, i.e. learning from the deep knowledge distilled from the teacher. More specifically, we first train a global learner on the entire base subset using category labels as supervision to leverage the global context information of the categories. Then, the well-trained global learner is used to simulate the query-support relatedness in global dependencies. Finally, the distilled global query-support relatedness is explicitly used to train the meta-learner using the RDD strategy, with the goal of making the meta-learner more discriminative. The RDD strategy aims to decouple the dense query-support relatedness into the groups of sparse decoupled relatedness. Moreover, only the relatedness of a single support sample with other query samples is considered in each group. By distilling the sparse decoupled relatedness group by group, sharper relatedness can be effectively distilled to the meta-learner, thereby facilitating the learning of a discriminative meta-learner. We conduct extensive experiments on the miniImagenet and CIFAR-FS datasets, which show the state-of-the-art performance of our GRDD method.


title: Implementing a LoRa Software-Defined Radio on a General-Purpose ULP Microcontroller

abstract: Emerging Internet-of-Things sensing applications rely on ultra low-power (ULP) microcontroller units (MCUs) that wirelessly transmit data to the cloud. Typical MCUs nowadays consist of generic blocks, except for the protocol-specific radios implemented in hardware. Hardware radios however slow down the evolution of wireless protocols due to retrocompatiblity concerns. In this work, we explore a software-defined radio architecture by demonstrating a LoRa transceiver running on custom ULP MCU codenamed SleepRider with an ARM Cortex-M4 CPU. In SleepRider MCU, we offload the generic baseband operations (e.g., low-pass filtering) to a reconfigurable digital front-end block and use the Cortex-M4 CPU to perform the protocol-specific computations. Our software implementation of the LoRa physical layer only uses the native SIMD instructions of the Cortex-M4 to achieve real-time transmission and reception of LoRa packets. SleepRider MCU has been fabricated in a 28nm FDSOI CMOS technology and is used in a testbed to experimentally validate the software implementation. Experimental results show that the proposed software-defined radio requires only a CPU frequency of 20 MHz to correctly receive a LoRa packet, with an ultra-low power consumption of 0.42 mW on average.


title: Differentially Private Stochastic Optimization: New Results in Convex and Non-Convex Settings

abstract: We study differentially private stochastic optimization in convex and non-convex settings. For the convex case, we focus on the family of non-smooth generalized linear losses (GLLs). Our algorithm for the $\ell_2$ setting achieves optimal excess population risk in near-linear time, while the best known differentially private algorithms for general convex losses run in super-linear time. Our algorithm for the $\ell_1$ setting has nearly-optimal excess population risk $\tilde{O}\big(\sqrt{\frac{\log{d}}{n\varepsilon}}\big)$, and circumvents the dimension dependent lower bound of \cite{Asi:2021} for general non-smooth convex losses. In the differentially private non-convex setting, we provide several new algorithms for approximating stationary points of the population risk. For the $\ell_1$-case with smooth losses and polyhedral constraint, we provide the first nearly dimension independent rate, $\tilde O\big(\frac{\log^{2/3}{d}}{{(n\varepsilon)^{1/3}}}\big)$ in linear time. For the constrained $\ell_2$-case with smooth losses, we obtain a linear-time algorithm with rate $\tilde O\big(\frac{1}{n^{1/3}}+\frac{d^{1/5}}{(n\varepsilon)^{2/5}}\big)$. Finally, for the $\ell_2$-case we provide the first method for {\em non-smooth weakly convex} stochastic optimization with rate $\tilde O\big(\frac{1}{n^{1/4}}+\frac{d^{1/6}}{(n\varepsilon)^{1/3}}\big)$ which matches the best existing non-private algorithm when $d= O(\sqrt{n})$. We also extend all our results above for the non-convex $\ell_2$ setting to the $\ell_p$ setting, where $1 < p \leq 2$, with only polylogarithmic (in the dimension) overhead in the rates.


title: Wide-field X-ray observations of the supernova remnant Puppis A with the SRG/ART-XC telescope

abstract: The Spectrum-Roentgen-Gamma (SRG) observatory is currently conducting its 4-year all-sky X-ray survey, started on December 12, 2019. The survey is periodically interrupted for technological operations with the spacecraft. These time intervals are usually used by the Mikhail Pavlinsky ART-XC telescope to perform calibrations. In this context, SRG carried out scanning observations of the Puppis A supernova remnant (SNR) with the aim to check the imaging performance of ART-XC and to optimize the technique of image reconstruction for extended objects. Using the unique imaging capabilities of ART-XC and its uniform coverage of the entire Puppis A region, we attempted to investigate the morphology of this SNR at energies >4 keV, and to search for previously unknown X-ray sources. Puppis A was observed in 2019-2020, conducting 1.5 x 1.5 degrees shallow surveys with an exposure of 36 hours. Additional deep pointed observations of the central part of Puppis A were carried out in 2021 lasted 31 hours to highlight the morphology of the extended emission. The X-ray emission of the Puppis A was significantly detected as an extended structure in the 4-6 keV energy band. The morphology of the emission is in general agreement with that observed in soft X-rays previously. The deep sky image of Puppis A obtained with the ART-XC telescope is characterized by a typical SNR shell rim morphology, an extended emission and a bright emission knot in the north-eastern part of the supernova shell. Also, four point X-ray sources have been detected, including three objects identified in catalogs, and one newly discovered X-ray emitter.


title: SRG/ART-XC discovery of SRGA J204318.2+443815: towards the complete population of faint X-ray pulsars

abstract: We report a discovery of a new long-period X-ray pulsar SRGA J204318.2+443815/SRGe J204319.0+443820 in the Be binary system. The source was found in the second all-sky survey by the Mikhail Pavlinsky telescope on board the SRG mission. The follow-up observations with XMM-Newton, NICER and NuSTAR observatories allowed us to discover a strong coherent signal in the source light curve with the period of $\sim742$ s. The pulsed fraction was found to depend on the energy increasing from $\sim20$% in soft X-rays to $>50$% at high energies, as it is typical for X-ray pulsars. The source demonstrate a quite hard spectrum with an exponential cutoff at high energies and bolometric luminosity of $L_X \simeq 4\times10^{35}$ erg/s. Dedicated optical and infrared observations with the RTT-150, NOT, Keck and Palomar telescopes revealed a number of emission lines (H$_{\alpha}$, HeI, Pashen and Braket series) with the strongly absorbed continuum. All of above suggests that SRGAJ204318.2+443815/ SRGeJ204319.0+443820 is a new persistent low luminosity X-ray pulsar in a distant binary system with a Be-star of the B0-B2e class. Thus the SRG observatory allow us to unveil the hidden population of faint persistent objects including the population of slowly rotating X-ray pulsars in Be systems.


title: Peculiar X-ray transient SRGA J043520.9+552226/AT2019wey discovered with SRG/ART-XC

abstract: Context: During the ongoing all-sky survey, the Mikhail Pavlinsky ART-XC telescope on board the SRG observatory should discover new X-ray sources, many of which can be transient. Here we report on the discovery and multiwavelength follow-up of a peculiar X-ray source SRGA J043520.9+552226=SRGe J043523.3+552234 - the high-energy counterpart of the optical transient AT2019wey. Aims: Thanks to its sensitivity and the survey strategy, the Mikhail Pavlinsky ART-XC telescope uncovers poorly studied weak transient populations. Using a synergy with current public optical surveys, we are aiming at revealing the nature of these transients to study its parent populations. The SRGA J043520.9+552226 is the first transient detected by ART-XC which has a bright optical counterpart suitable for further studies. Methods: We have used available public X-ray and optical data and observations with SRG, INTEGRAL, NuSTAR, NICER and ground-based telescopes to investigate the source spectral energy distributions at different phases of the outburst. Results: Based on X-ray spectral and timing properties derived from space observations, optical spectroscopy and photometry obtained with the 2.5-m and RC600 CMO SAI MSU telescopes, we propose the source to be a black hole in a low-mass close X-ray binary system.


title: Losing a minute every two years: SRG X-ray view on the rapidly accelerating X-ray pulsar SXP1323

abstract: SXP 1323 is a peculiar high-mass X-ray binary located in the Small Magellanic Cloud, renowned for its rapid spin-up. We investigate for the first time broadband X-ray properties of SXP1323, as observed by the Mikhail Pavlinsky ART-XC and eRosita telescopes on board the SRG observatory. Using ART-XC and eRosita, data we produced first broadband 1-20 keV X-ray spectrum and estimated pulsed fraction above 8 keV. With the addition of archival XMM-Newton observations we traced evolution of the SXP 1323 spin period over the last five years and found that after 2016 the source switched to a linear spin-up with rate of -29.9 s yr$^{-1}$. Broadband X-ray spectrum is typical for accreting X-ray pulsars, with steep powerlaw index ($\Gamma$=-0.15) and exponential cutoff energy of 5.1 keV. No significant difference between spectra obtained in states with and without pulsations were found.


title: An equation-of-state-meter for CBM using PointNet

abstract: A novel method for identifying the nature of QCD transitions in heavy-ion collision experiments is introduced. PointNet based Deep Learning (DL) models are developed to classify the equation of state (EoS) that drives the hydrodynamic evolution of the system created in Au-Au collisions at 10 AGeV. The DL models were trained and evaluated in different hypothetical experimental situations. A decreased performance is observed when more realistic experimental effects (acceptance cuts and decreased resolutions) are taken into account. It is shown that the performance can be improved by combining multiple events to make predictions. The PointNet based models trained on the reconstructed tracks of charged particles from the CBM detector simulation discriminate a crossover transition from a first order phase transition with an accuracy of up to 99.8%. The models were subjected to several tests to evaluate the dependence of its performance on the centrality of the collisions and physical parameters of fluid dynamic simulations. The models are shown to work in a broad range of centralities (b=0-7 fm). However, the performance is found to improve for central collisions (b=0-3 fm). There is a drop in the performance when the model parameters lead to reduced duration of the fluid dynamic evolution or when less fraction of the medium undergoes the transition. These effects are due to the limitations of the underlying physics and the DL models are shown to be superior in its discrimination performance in comparison to conventional mean observables.


title: Origin of Perpendicular Magnetic Anisotropy in Yttrium Iron Garnet Thin Films Grown on Si (100)

abstract: We report the magnetic properties of yttrium iron garnet (YIG) thin films grown by pulsed laser deposition technique. The films were deposited on Si (100) substrates in the range of 15-50 nm thickness. Magnetic characterizations were investigated by ferromagnetic resonance spectra. Perpendicular magnetic easy axis was achieved up to 50 nm thickness. We observed that the perpendicular anisotropy values decreased by increasing the film thickness. The origin of the perpendicular magnetic anisotropy (PMA) was attributed to the texture and the lattice distortion in the YIG thin films. We anticipate that perpendicularly magnetized YIG thin films on Si substrates pave the way for a cheaper and compatible fabrication process.


title: Investor Behavior Modeling by Analyzing Financial Advisor Notes: A Machine Learning Perspective

abstract: Modeling investor behavior is crucial to identifying behavioral coaching opportunities for financial advisors. With the help of natural language processing (NLP) we analyze an unstructured (textual) dataset of financial advisors' summary notes, taken after every investor conversation, to gain first ever insights into advisor-investor interactions. These insights are used to predict investor needs during adverse market conditions; thus allowing advisors to coach investors and help avoid inappropriate financial decision-making. First, we perform topic modeling to gain insight into the emerging topics and trends. Based on this insight, we construct a supervised classification model to predict the probability that an advised investor will require behavioral coaching during volatile market periods. To the best of our knowledge, ours is the first work on exploring the advisor-investor relationship using unstructured data. This work may have far-reaching implications for both traditional and emerging financial advisory service models like robo-advising.


title: Leveraging Explainability for Comprehending Referring Expressions in the Real World

abstract: For effective human-robot collaboration, it is crucial for robots to understand requests from users and ask reasonable follow-up questions when there are ambiguities. While comprehending the users' object descriptions in the requests, existing studies have focused on this challenge for limited object categories that can be detected or localized with existing object detection and localization modules. On the other hand, in the wild, it is impossible to limit the object categories that can be encountered during the interaction. To understand described objects and resolve ambiguities in the wild, for the first time, we suggest a method by leveraging explainability. Our method focuses on the active regions of a scene to find the described objects without putting the previous constraints on object categories and natural language instructions. We evaluate our method in varied real-world images and observe that the regions suggested by our method can help resolve ambiguities. When we compare our method with a state-of-the-art baseline, we show that our method performs better in scenes with ambiguous objects which cannot be recognized by existing object detectors.


title: Repulsive dynamics of strongly attractive one-dimensional quantum gases

abstract: We analyze the dynamics of one-dimensional quantum gases with strongly attractive contact interactions. We specify a class of initial states for which attractive forces effectively act as strongly repulsive ones during the time evolution. Our findings extend the theoretical results on the super-Tonks-Girardeau gas to a highly nonequilibrium dynamics. The novel mechanism is illustrated on the prototypical problem of the domain stability in a two-component Fermi gas. We also discuss nonlocal interactions and analyze universality of the presented results. Moreover, we use our conclusions to argue for the existence of metastable quantum droplets in the regime of strongly attractive contact and attractive dipolar interactions.


title: Coloring graphs with forbidden bipartite subgraphs

abstract: A conjecture of Alon, Krivelevich, and Sudakov states that, for any graph $F$, there is a constant $c_F > 0$ such that if $G$ is an $F$-free graph of maximum degree $\Delta$, then $\chi(G) \leq c_F \Delta / \log\Delta$. Alon, Krivelevich, and Sudakov verified this conjecture for a class of graphs $F$ that includes all bipartite graphs. Moreover, it follows from recent work by Davies, Kang, Pirot, and Sereni that if $G$ is $K_{t,t}$-free, then $\chi(G) \leq (t + o(1)) \Delta / \log\Delta$ as $\Delta \to \infty$. We improve this bound to $(1+o(1)) \Delta/\log \Delta$, making the constant factor independent of $t$. We further extend our result to the DP-coloring setting (also known as correspondence coloring), introduced by Dvo\v{r}\'ak and Postle.


title: Stellar intensity interferometry of Vega in photon counting mode

abstract: Stellar Intensity Interferometry is a technique based on the measurement of the second order spatial correlation of the light emitted from a star. The physical information provided by these measurements is the angular size and structure of the emitting source. A worldwide effort is presently under way to implement stellar intensity interferometry on telescopes separated by long baselines and on future arrays of Cherenkov telescopes. We describe an experiment of this type, realized at the Asiago Observatory (Italy), in which we performed for the first time measurements of the correlation counting photon coincidences in post-processing by means of a single photon software correlator and exploiting entirely the quantum properties of the light emitted from a star. We successfully detected the temporal correlation of Vega at zero baseline and performed a measurement of the correlation on a projected baseline of $\sim$2 km. The average discrete degree of coherence at zero baseline for Vega is $< g^{(2)} > \, = 1.0034 \pm 0.0008$, providing a detection with a signal-to-noise ratio $S/N \gtrsim 4$. No correlation is detected over the km baseline. The measurements are consistent with the expected degree of spatial coherence for a source with the 3.3 mas angular diameter of Vega. The experience gained with the Asiago experiment will serve for future implementations of stellar intensity interferometry on long-baseline arrays of Cherenkov telescopes.


title: A Five-Planet Resonant Chain: Reevaluation of the Kepler-80 System

abstract: Since the launch of the Kepler space telescope in 2009 and the subsequent K2 mission, hundreds of multi-planet systems have been discovered. The study of such systems, both as individual systems and as a population, leads to a better understanding of planetary formation and evolution. Kepler-80, a K-dwarf hosting six super-Earths, was the first system known to have four planets in a chain of resonances, a repeated geometric configuration. Transiting planets in resonant chains can enable us to estimate not only the planets' orbits and sizes but also their masses. Since the original resonance analysis and TTV fitting of Kepler-80, a new planet has been discovered whose signal likely altered the measured masses of the other planets. Here, we determine masses and orbits for all six planets hosted by Kepler-80 by direct forward photodynamical modeling of the lightcurve of this system. We then explore the resonant behaviour of the system. We find that the four middle planets are in a resonant chain, but that the outermost planet only dynamically interacts in $\sim14$\% of our solutions. We also find that the system and its dynamic behaviour are consistent with \emph{in situ} formation and compare our results to two other resonant chain systems, Kepler-60 and TRAPPIST-1.


title: Nonlinear Least Squares for Large-Scale Machine Learning using Stochastic Jacobian Estimates

abstract: For large nonlinear least squares loss functions in machine learning we exploit the property that the number of model parameters typically exceeds the data in one batch. This implies a low-rank structure in the Hessian of the loss, which enables effective means to compute search directions. Using this property, we develop two algorithms that estimate Jacobian matrices and perform well when compared to state-of-the-art methods.


title: Active Divergence with Generative Deep Learning -- A Survey and Taxonomy

abstract: Generative deep learning systems offer powerful tools for artefact generation, given their ability to model distributions of data and generate high-fidelity results. In the context of computational creativity, however, a major shortcoming is that they are unable to explicitly diverge from the training data in creative ways and are limited to fitting the target data distribution. To address these limitations, there have been a growing number of approaches for optimising, hacking and rewriting these models in order to actively diverge from the training data. We present a taxonomy and comprehensive survey of the state of the art of active divergence techniques, highlighting the potential for computational creativity researchers to advance these methods and use deep generative models in truly creative systems.


title: Double and Charge-Transfer Excitations in Time-Dependent Density Functional Theory

abstract: Time-dependent density functional theory has emerged as a method of choice for calculations of spectra and response properties in physics, chemistry, and biology, with its system-size scaling enabling computations on systems much larger than possible otherwise. While increasingly complex and interesting systems have been successfully tackled with relatively simple functional approximations, there has also been increasing awareness that these functionals tend to fail for certain classes of approximations. I review the fundamental challenges the approximate functionals have in describing double-excitations and charge-transfer excitations, which are two of the most common impediments for the theory to be applied in a black box way. At the same time, I describe the progress made in recent decades in developing functional approximations that give useful predictions for these excitations.


title: Statistics of galaxy mergers: bridging the gap between theory and observation

abstract: We present a study of galaxy mergers up to $z=10$ using the Planck Millennium cosmological dark matter simulation and the {\tt GALFORM} semi-analytical model of galaxy formation. Utilising the full ($800$ Mpc)$^3$ volume of the simulation, we studied the statistics of galaxy mergers in terms of merger rates and close pair fractions. We predict that merger rates begin to drop rapidly for high-mass galaxies ($M_*>10^{11.3}-10^{10.5}$ $M_\odot$ for $z=0-4$), as a result of the exponential decline in the galaxy stellar mass function. The predicted merger rates increase and then turn over with increasing redshift, by $z=3.5$, in disagreement with hydrodynamical simulations and semi-empirical models. In agreement with most other models and observations, we find that close pair fractions flatten or turn over at some redshift (dependent on the mass selection). We conduct an extensive comparison of close pair fractions, and highlight inconsistencies among models, but also between different observations. We provide a fitting formula for the major merger timescale for close galaxy pairs, in which the slope of the stellar mass dependence is redshift dependent. This is in disagreement with previous theoretical results that implied a constant slope. Instead we find a weak redshift dependence only for massive galaxies ($M_*>10^{10}$ M$_\odot$): in this case the merger timescale varies approximately as $M_*^{-0.55}$. We find that close pair fractions and merger timescales depend on the maximum projected separation as $r_\mathrm{max}^{1.32}$. This is in agreement with observations of small-scale clustering of galaxies, but is at odds with the linear dependence on projected separation that is often assumed.


title: Understanding the Communist Party of China's Information Operations

abstract: The Communist Party of China is known to engage in Information Operations to influence public opinion. In this paper, we seek to understand the tactics used by the Communist Party in a recent Information Operation - the one conducted to influence the narrative around the pro-democracy movement in Hong Kong. We use a Twitter dataset containing account information and tweets for the operation. Our research shows that the Hong Kong operation was (at least) partially conducted manually by humans rather than entirely by automated bots. We also show that the Communist Party mixed in personal attacks on Chinese dissidents and messages on COVID-19 with the party's views on the protests during the operation. Finally, we conclude that the Information Operation network in the Twitter dataset was set up to amplify content generated elsewhere rather than to influence the narrative with original content.


title: Positive mass theorems for asymptotically hyperbolic Riemannian manifolds with boundary

abstract: We prove positive mass theorems for asymptotically hyperbolic and asymptotically locally hyperbolic Riemannian manifolds with black-hole-type boundaries.


title: Direct speech-to-speech translation with discrete units

abstract: We present a direct speech-to-speech translation (S2ST) model that translates speech from one language to speech in another language without relying on intermediate text generation. We tackle the problem by first applying a self-supervised discrete speech encoder on the target speech and then training a sequence-to-sequence speech-to-unit translation (S2UT) model to predict the discrete representations of the target speech. When target text transcripts are available, we design a joint speech and text training framework that enables the model to generate dual modality output (speech and text) simultaneously in the same inference pass. Experiments on the Fisher Spanish-English dataset show that the proposed framework yields improvement of 6.7 BLEU compared with a baseline direct S2ST model that predicts spectrogram features. When trained without any text transcripts, our model performance is comparable to models that predict spectrograms and are trained with text supervision, showing the potential of our system for translation between unwritten languages. Audio samples are available at https://facebookresearch.github.io/speech_translation/direct_s2st_units/index.html .


title: Interpretable Mammographic Image Classification using Case-Based Reasoning and Deep Learning

abstract: When we deploy machine learning models in high-stakes medical settings, we must ensure these models make accurate predictions that are consistent with known medical science. Inherently interpretable networks address this need by explaining the rationale behind each decision while maintaining equal or higher accuracy compared to black-box models. In this work, we present a novel interpretable neural network algorithm that uses case-based reasoning for mammography. Designed to aid a radiologist in their decisions, our network presents both a prediction of malignancy and an explanation of that prediction using known medical features. In order to yield helpful explanations, the network is designed to mimic the reasoning processes of a radiologist: our network first detects the clinically relevant semantic features of each image by comparing each new image with a learned set of prototypical image parts from the training images, then uses those clinical features to predict malignancy. Compared to other methods, our model detects clinical features (mass margins) with equal or higher accuracy, provides a more detailed explanation of its prediction, and is better able to differentiate the classification-relevant parts of the image.


title: First direct measurement of $^{59}$Cu(p,$\alpha$)$^{56}$Ni: A step towards constraining the Ni-Cu cycle in the Cosmos

abstract: Reactions on the proton-rich nuclides drive the nucleosynthesis in Core-Collapse Supernovae (CCSNe) and in X-ray bursts (XRBs). CCSNe eject the nucleosynthesis products to the interstellar medium and hence are a potential inventory of p-nuclei, whereas in XRBs nucleosynthesis powers the light curves. In both astrophysical sites the Ni-Cu cycle, which features a competition between $^{59}$Cu(p,$\alpha$)$^{56}$Ni and $^{59}$Cu(p,$\gamma$)$^{60}$Zn, could potentially halt the production of heavier elements. Here, we report the first direct measurement of $^{59}$Cu(p,$\alpha$)$^{56}$Ni using a re-accelerated $^{59}$Cu beam and cryogenic solid hydrogen target. Our results show that the reaction proceeds predominantly to the ground state of $^{56}$Ni and the experimental rate has been found to be lower than Hauser-Feshbach-based statistical predictions. New results hint that the $\nu p$-process could operate at higher temperatures than previously inferred and therefore remains a viable site for synthesizing the heavier elements.


title: On strong integrability of the dressing cosets

abstract: We formulate sufficient conditions for the strong integrability of dressing cosets. We provide several sigma-model backgrounds solving those conditions, some of them are new and some of them were not so far formulated as the dressing cosets. The new models are based on the Drinfeld doubles having the structure of higher order jet bundles of quadratic Lie groups.


title: Non-abelian Mellin Transformations and Applications

abstract: We study non-abelian versions of the Mellin transformations, originally introduced by Gabber-Loeser on complex affine tori. Our main result is a generalization to the non-abelian context and with arbitrary coefficients of the t-exactness of Gabber-Loeser's Mellin transformation. As an intermediate step, we obtain vanishing results for the Sabbah specialization functors. Our main application is to construct new examples of duality spaces in the sense of Bieri-Eckmann, generalizing results of Denham-Suciu.


title: Analyzing black-hole ringdowns

abstract: A perturbed black hole rings down by emitting gravitational waves in tones with specific frequencies and durations. Such tones encode prized information about the geometry of the source spacetime and the fundamental nature of gravity, making the measurement of black hole ringdowns a key goal of gravitational wave astronomy. However, this task is plagued by technical challenges that invalidate the naive application of standard data analysis methods and complicate sensitivity projections. In this paper, we provide a comprehensive account of the formalism required to properly carry out ringdown analyses, examining in detail the foundations of recent observational results, and providing a framework for future measurements. We build on those insights to clarify the concepts of ringdown detectability and resolvability -- touching on the drawbacks of both Bayes factors and naive Fisher matrix approaches -- and find that overly pessimistic heuristics have led previous works to underestimate the role of ringdown overtones for black hole spectroscopy. We put our framework to work on the analysis of a variety of simulated signals in colored noise, including analytic injections and a numerical relativity simulation consistent with GW150914. We demonstrate that we can use tones of the quadrupolar angular harmonic to test the no-hair theorem at current sensitivity, with precision comparable to published constraints from real data. Finally, we assess the role of modeling systematics, and project measurements for future, louder signals. We release ringdown, a Python library for analyzing black hole ringdowns using the the methods discussed in this paper, under a permissive open-source license at https://github.com/maxisi/ringdown


title: The Fundamental Theorem of Natural Selection

abstract: Suppose we have $n$ different types of self-replicating entity, with the population $P_i$ of the $i$th type changing at a rate equal to $P_i$ times the fitness $f_i$ of that type. Suppose the fitness $f_i$ is any continuous function of all the populations $P_1, \dots, P_n$. Let $p_i$ be the fraction of replicators that are of the $i$th type. Then $p = (p_1, \dots, p_n)$ is a time-dependent probability distribution, and we prove that its speed as measured by the Fisher information metric equals the variance in fitness. In rough terms, this says that the speed at which information is updated through natural selection equals the variance in fitness. This result can be seen as a modified version of Fisher's fundamental theorem of natural selection. We compare it to Fisher's original result as interpreted by Price, Ewens and Edwards.


title: Identification of 3 cataclysmic variables detected by the ART-XC and eROSITA telescopes aboard SRG during the all-sky X-ray survey

abstract: We report the discovery of three previously unknown cataclysmic variables in the data of the first year of the all-sky X-ray survey by the SRG orbital observatory. The sources were selected due to their brightness in the 4--12 keV band in the data of the Mikhail Pavlinsky ART-XC telescope. They are also detected by the eROSITA telescope, which provided accurate localizations and spectral data for broad-band spectral analysis. All three objects had been previously known as X-ray sources from the ROSAT all-sky survey and XMM-Newton slew survey, but their nature remained unknown. The X-ray spectra obtained by eROSITA and ART-XC are consistent with optically thin thermal emission with a temperature kT>~15 keV for SRGAJ194638.9+704552 and SRGAJ225412.8+690658 and kT>~5 keV for SRGAJ204547.8+672642. This, together with the inferred high X-ray luminosities ($2\times 10^{32}$-$3\times 10^{33}$ erg s$^{-1}$), strongly suggests that all three sources are CVs. We have obtained optical photometry and spectroscopy for these objects using the AZT-33IK 1.6-m telescope of the Sayan Observatory. The optical properties confirm the CV nature of the objects. We conclude that SRGAJ194638.9+704552 is an intermediate polar, SRGAJ204547.8+672642 is most likely a polar or an intermediate polar, and SRGAJ225412.8+690658 can be either a magnetic or a non-magnetic CV. We also measured an orbital period of 2.98~hours for SRGAJ204547.8+672642, based on TESS data. Three out of the planned eight SRG all-sky surveys have now been completed. We expect to find plenty of new CVs during the survey and to continue our optical follow-up program.


title: A Persistent Spatial Semantic Representation for High-level Natural Language Instruction Execution

abstract: Natural language provides an accessible and expressive interface to specify long-term tasks for robotic agents. However, non-experts are likely to specify such tasks with high-level instructions, which abstract over specific robot actions through several layers of abstraction. We propose that key to bridging this gap between language and robot actions over long execution horizons are persistent representations. We propose a persistent spatial semantic representation method, and show how it enables building an agent that performs hierarchical reasoning to effectively execute long-term tasks. We evaluate our approach on the ALFRED benchmark and achieve state-of-the-art results, despite completely avoiding the commonly used step-by-step instructions.


title: Parallel Element-based Algebraic Multigrid for H(curl) and H(div) Problems Using the ParELAG Library

abstract: This paper presents the use of element-based algebraic multigrid (AMGe) hierarchies, implemented in the ParELAG (Parallel Element Agglomeration Algebraic Multigrid Upscaling and Solvers) library, to produce multilevel preconditioners and solvers for H(curl) and H(div) formulations. ParELAG constructs hierarchies of compatible nested spaces, forming an exact de Rham sequence on each level. This allows the application of hybrid smoothers on all levels and AMS (Auxiliary-space Maxwell Solver) or ADS (Auxiliary-space Divergence Solver) on the coarsest levels, obtaining complete multigrid cycles. Numerical results are presented, showing the parallel performance of the proposed methods. As a part of the exposition, this paper demonstrates some of the capabilities of ParELAG and outlines some of the components and procedures within the library.


title: Quantum noise and vacuum fluctuations in balanced homodyne detections through ideal multi-mode detectors

abstract: The balanced homodyne detection as a readout scheme of gravitational-wave detectors is carefully examined from the quantum field theoretical point of view. The readout scheme in gravitational-wave detectors specifies the directly measured quantum operator in the detection. This specification is necessary when we apply the recently developed quantum measurement theory to gravitational-wave detections. We examine the two models of measurement. One is the model in which the directly measured quantum operator at the photodetector is Glauber's photon number operator, and the other is the model in which the power operator of the optical field is directly measured. These two are regarded as ideal models of photodetectors. We first show these two models yield the same expectation value of the measurement. Since it is consensus in the gravitational-wave community that vacuum fluctuations contribute to the noises in the detectors, we also clarify the contributions of vacuum fluctuations to the quantum noise spectral density without using the two-photon formulation which is used in the gravitational-wave community. We found that the conventional noise spectral density in the two-photon formulation includes vacuum fluctuations from the main interferometer but does not include those from the local oscillator. Although the contribution of vacuum fluctuations from the local oscillator theoretically yields the difference between the above two models in the noise spectral densities, this difference is negligible in realistic situations.


title: Measurement of the Masses and Widths of the $\Sigma_c(2455)^+$ and $\Sigma_c(2520)^+$ Baryons

abstract: Using 980 ${\rm fb}^{-1}$ of data {collected} with the Belle detector operating at the KEKB asymmetric-energy $e^+e^-$ collider, we report the measurements of the masses, and the first measurements of the instrinsic widths, of the $\Sigma_c(2455)^+$ and $\Sigma_c(2520)^+$ charmed baryons. We find $M(\Sigma_c(2455)^+)-M(\Lambda_c^+) = 166.17\pm 0.05^{+0.16}_{-0.07}\ {\rm MeV}/c^2$, $\Gamma(\Sigma_c(2455)^+) = 2.3 \pm 0.3 \pm 0.3\ {\rm MeV/c^2}$, $M(\Sigma_c(2520)^+)-M(\Lambda_c^+) = 230.9 \pm 0.5 ^{+0.5}_{-0.1}\ {\rm MeV}/c^2$, and $\Gamma(\Sigma_c(2520)^+) = 17.2^{+2.3\ +3.1}_{-2.1\ -0.7}\ {\rm MeV}/c^2$, where the uncertainties are statistical and systematic, respectively. These measurements can be used to test models of the underlying quark structure of the $\Sigma_c$ states.


title: Fast Contact-Implicit Model-Predictive Control

abstract: We present a general approach for controlling robotic systems that make and break contact with their environments. Contact-implicit model predictive control (CI-MPC) generalizes linear MPC to contact-rich settings by utilizing a bi-level planning formulation with lower-level contact dynamics formulated as time-varying linear complementarity problems (LCPs) computed using strategic Taylor approximations about a reference trajectory. These dynamics enable the upper-level planning problem to reason about contact timing and forces, and generate entirely new contact-mode sequences online. To achieve reliable and fast numerical convergence, we devise a structure-exploiting interior-point solver for these LCP contact dynamics and a custom trajectory optimizer for the tracking problem. We demonstrate real-time solution rates for CI-MPC and the ability to generate and track non-periodic behaviours in hardware experiments on a quadrupedal robot. We also show that the controller is robust to model mismatch and can respond to disturbances by discovering and exploiting new contact modes across a variety of robotic systems in simulation, including a pushbot, planar hopper, planar quadruped, and planar biped.


title: Let's Play for Action: Recognizing Activities of Daily Living by Learning from Life Simulation Video Games

abstract: Recognizing Activities of Daily Living (ADL) is a vital process for intelligent assistive robots, but collecting large annotated datasets requires time-consuming temporal labeling and raises privacy concerns, e.g., if the data is collected in a real household. In this work, we explore the concept of constructing training examples for ADL recognition by playing life simulation video games and introduce the SIMS4ACTION dataset created with the popular commercial game THE SIMS 4. We build Sims4Action by specifically executing actions-of-interest in a "top-down" manner, while the gaming circumstances allow us to freely switch between environments, camera angles and subject appearances. While ADL recognition on gaming data is interesting from the theoretical perspective, the key challenge arises from transferring it to the real-world applications, such as smart-homes or assistive robotics. To meet this requirement, Sims4Action is accompanied with a GamingToReal benchmark, where the models are evaluated on real videos derived from an existing ADL dataset. We integrate two modern algorithms for video-based activity recognition in our framework, revealing the value of life simulation video games as an inexpensive and far less intrusive source of training data. However, our results also indicate that tasks involving a mixture of gaming and real data are challenging, opening a new research direction. We will make our dataset publicly available at https://github.com/aroitberg/sims4action.


title: Exponents of Diophantine approximation in dimension two for a general class of numbers

abstract: We study the Diophantine properties of a new class of transcendental real numbers which contains, among others, Roy's extremal numbers, Bugeaud-Laurent Sturmian continued fractions, and more generally the class of Sturmian type numbers. We compute, for each real number $\xi$ of this set, several exponents of Diophantine approximation to the pair $(\xi,\xi^2)$, together with $\omega_2^*(\xi)$ and $\widehat{\omega}_2^*(\xi)$, the so-called ordinary and uniform exponent of approximation to $\xi$ by algebraic numbers of degree $\leq 2$. As an application, we get new information on the set of values taken by $\widehat{\omega}_2^*$ at transcendental numbers, and we give a partial answer to a question of Fischler about his exponent $\beta_0$.


title: Statistical Modeling for Practical Pooled Testing During the COVID-19 Pandemic

abstract: Pooled testing offers an efficient solution to the unprecedented testing demands of the COVID-19 pandemic, although with potentially lower sensitivity and increased costs to implementation in some settings. Assessments of this trade-off typically assume pooled specimens are independent and identically distributed. Yet, in the context of COVID-19, these assumptions are often violated: testing done on networks (housemates, spouses, co-workers) captures correlated individuals, while infection risk varies substantially across time, place and individuals. Neglecting dependencies and heterogeneity may bias established optimality grids and induce a sub-optimal implementation of the procedure. As a lesson learned from this pandemic, this paper highlights the necessity of integrating field sampling information with statistical modeling to efficiently optimize pooled testing. Using real data, we show that (a) greater gains can be achieved at low logistical cost by exploiting natural correlations (non-independence) between samples -- allowing improvements in sensitivity and efficiency of up to 30% and 90% respectively; and (b) these gains are robust despite substantial heterogeneity across pools (non-identical). Our modeling results complement and extend the observations of Barak et al (2021) who report an empirical sensitivity well beyond expectations. Finally, we provide an interactive tool for selecting an optimal pool size using contextual information