title: Design Considerations for Data Daemons: Co-creating Design Futures to Explore Ethical Personal Data Management

abstract: Mobile applications and online service providers track our virtual and physical behaviour more actively and with a broader scope than ever before. This has given rise to growing concerns about ethical personal data management. Even though regulation and awareness around data ethics are increasing, end-users are seldom engaged when defining and designing what a future with ethical personal data management should look like. We explore a participatory process that uses design futures, the Future workshop method and design fictions to envision ethical personal data management with end-users and designers. To engage participants effectively, we needed to bridge their differential expertise and make the abstract concepts of data and ethics tangible. By concretely presenting personal data management and control as fictitious entities called Data Daemons, we created a shared understanding of these abstract concepts, and empowered non-expert end-users and designers to become actively engaged in the design process.


title: Federated Dynamic Spectrum Access

abstract: Due to the growing volume of data traffic produced by the surge of Internet of Things (IoT) devices, the demand for radio spectrum resources is approaching their limitation defined by Federal Communications Commission (FCC). To this end, Dynamic Spectrum Access (DSA) is considered as a promising technology to handle this spectrum scarcity. However, standard DSA techniques often rely on analytical modeling wireless networks, making its application intractable in under-measured network environments. Therefore, utilizing neural networks to approximate the network dynamics is an alternative approach. In this article, we introduce a Federated Learning (FL) based framework for the task of DSA, where FL is a distributive machine learning framework that can reserve the privacy of network terminals under heterogeneous data distributions. We discuss the opportunities, challenges, and opening problems of this framework. To evaluate its feasibility, we implement a Multi-Agent Reinforcement Learning (MARL)-based FL as a realization associated with its initial evaluation results.


title: The Food Recognition Benchmark: Using DeepLearning to Recognize Food on Images

abstract: The automatic recognition of food on images has numerous interesting applications, including nutritional tracking in medical cohorts. The problem has received significant research attention, but an ongoing public benchmark to develop open and reproducible algorithms has been missing. Here, we report on the setup of such a benchmark using publicly available food images sourced through the mobile MyFoodRepo app. Through four rounds, the benchmark released the MyFoodRepo-273 dataset constituting 24,119 images and a total of 39,325 segmented polygons categorized in 273 different classes. Models were evaluated on private tests sets from the same platform with 5,000 images and 7,865 annotations in the final round. Top-performing models on the 273 food categories reached a mean average precision of 0.568 (round 4) and a mean average recall of 0.885 (round 3). We present experimental validation of round 4 results, and discuss implications of the benchmark setup designed to increase the size and diversity of the dataset for future rounds.


title: Topology of contact points in Lieb-kagom\'e model

abstract: We analyse Lieb-kagom\'e model, a three-band model with contact points showing particular examples of the merging of Dirac contact points. We prove that eigenstates can be parametrized in a classification surface, which is a hypersurface of a 4-dimension space. This classification surface is a powerful device giving topological properties of the energy band structure; the analysis of its fundamental group proves that all singularities of the band structure can be characterized by four independent winding (integer) numbers. Lieb case separates: its classification surface differs and there is only one winding number.


title: On component interactions in two-stage recommender systems

abstract: Thanks to their scalability, two-stage recommenders are used by many of today's largest online platforms, including YouTube, LinkedIn, and Pinterest. These systems produce recommendations in two steps: (i) multiple nominators, tuned for low prediction latency, preselect a small subset of candidates from the whole item pool; (ii) a slower but more accurate ranker further narrows down the nominated items, and serves to the user. Despite their popularity, the literature on two-stage recommenders is relatively scarce, and the algorithms are often treated as mere sums of their parts. Such treatment presupposes that the two-stage performance is explained by the behavior of the individual components in isolation. This is not the case: using synthetic and real-world data, we demonstrate that interactions between the ranker and the nominators substantially affect the overall performance. Motivated by these findings, we derive a generalization lower bound which shows that independent nominator training can lead to performance on par with uniformly random recommendations. We find that careful design of item pools, each assigned to a different nominator, alleviates these issues. As manual search for a good pool allocation is difficult, we propose to learn one instead using a Mixture-of-Experts based approach. This significantly improves both precision and recall at K.


title: Notes on $\{a,b,c\}$-Modular Matrices

abstract: Let $A \in \mathbb{Z}^{m \times n}$ be an integral matrix and $a$, $b$, $c \in \mathbb{Z}$ satisfy $a \geq b \geq c \geq 0$. The question is to recognize whether $A$ is $\{a,b,c\}$-modular, i.e., whether the set of $n \times n$ subdeterminants of $A$ in absolute value is $\{a,b,c\}$. We will succeed in solving this problem in polynomial time unless $A$ possesses a duplicative relation, that is, $A$ has nonzero $n \times n$ subdeterminants $k_1$ and $k_2$ satisfying $2 \cdot |k_1| = |k_2|$. This is an extension of the well-known recognition algorithm for totally unimodular matrices. As a consequence of our analysis, we present a polynomial time algorithm to solve integer programs in standard form over $\{a,b,c\}$-modular constraint matrices for any constants $a$, $b$ and $c$.


title: Fast Bayesian Variable Selection in Binomial and Negative Binomial Regression

abstract: Bayesian variable selection is a powerful tool for data analysis, as it offers a principled method for variable selection that accounts for prior information and uncertainty. However, wider adoption of Bayesian variable selection has been hampered by computational challenges, especially in difficult regimes with a large number of covariates or non-conjugate likelihoods. Generalized linear models for count data, which are prevalent in biology, ecology, economics, and beyond, represent an important special case. Here we introduce an efficient MCMC scheme for variable selection in binomial and negative binomial regression that exploits Tempered Gibbs Sampling (Zanella and Roberts, 2019) and that includes logistic regression as a special case. In experiments we demonstrate the effectiveness of our approach, including on cancer data with seventeen thousand covariates.


title: Epidemiological theory of virus variants

abstract: We propose a physical theory underlying the temporal evolution of competing virus variants that relies on the existence of (quasi) fixed points capturing the large time scale invariance of the dynamics. To motivate our result we first modify the time-honoured compartmental models of the SIR type to account for the existence of competing variants and then show how their evolution can be naturally re-phrased in terms of flow equations ending at quasi fixed points. As the natural next step we employ (near) scale invariance to organise the time evolution of the competing variants within the effective description of the epidemic Renormalization Group framework. We test the resulting theory against the time evolution of COVID-19 virus variants that validate the theory empirically.


title: Creating Quantum Emitters in Hexagonal Boron Nitride Deterministically on Chip-Compatible Substrates

abstract: Two-dimensional hexagonal boron nitride (hBN) that hosts bright room-temperature single-photon emitters (SPEs) is a promising material platform for quantum information applications. An important step towards the practical application of hBN is the on-demand, position-controlled generation of SPEs. Several strategies have been reported to achieve the deterministic creation of hBN SPEs. However, they either rely on a substrate nanopatterning procedure that is not compatible with integrated photonic devices or utilize a radiation source that might cause unpredictable damage to hBN and underlying substrates. Here, we report a radiation- and lithography-free route to deterministically activate hBN SPEs by nanoindentation with an atomic force microscope (AFM) tip. The method is applied to thin hBN flakes (less than 25 nm in thickness) on flat silicon-dioxide-silicon substrates that can be readily integrated into on-chip photonic devices. The achieved SPEs yields are above 30% by utilizing multiple indent sizes, and a maximum yield of 36% is demonstrated for the indent size of around 400 nm. Our results mark an important step towards the deterministic creation and integration of hBN SPEs with photonic and plasmonic on-chip devices.


title: Online Estimation and Coverage Control with Heterogeneous Sensing Information

abstract: Heterogeneous multi-robot sensing systems are able to characterize physical processes more comprehensively than homogeneous systems. Access to multiple modalities of sensory data allow such systems to fuse information between complementary sources and learn richer representations of a phenomenon of interest. Often, these data are correlated but vary in fidelity, i.e., accuracy (bias) and precision (noise). Low-fidelity data may be more plentiful, while high-fidelity data may be more trustworthy. In this paper, we address the problem of multi-robot online estimation and coverage control by combining low- and high-fidelity data to learn and cover a sensory function of interest. We propose two algorithms for this task of heterogeneous learning and coverage -- namely Stochastic Sequencing of Multi-fidelity Learning and Coverage (SMLC) and Deterministic Sequencing of Multi-fidelity Learning and Coverage (DMLC) -- and prove that they converge asymptotically. In addition, we demonstrate the empirical efficacy of SMLC and DMLC through numerical simulations.


title: Sensitivity of airborne transmission of enveloped viruses to seasonal variation in indoor relative humidity

abstract: In temperate climates, infection rates of enveloped viruses peak during the winter. While these seasonal trends are established in influenza and human coronaviruses, the mechanisms driving the variation remain poorly understood and thus difficult to extend to similar viruses like SARS-CoV-2. In this study, we use the Quadrature-based model of Respiratory Aerosol and Droplets (QuaRAD) to explore the sensitivity of airborne transmission to the seasonal variation in indoor relative humidity across the wide range of relevant conditions, using SARS-CoV-2 as an example. Relative humidity impacts the evaporation rate and equilibrium size of airborne particles, which in turn may impact particle removal rates and virion viability. Across a large ensemble of scenarios, we found that the dry indoor conditions typical of the winter season lead to slower inactivation than in the more humid summer season; in poorly ventilated spaces, this reduction in inactivation rates increases the concentration of active virions, but this effect was important when the susceptible person was farther than 2 m downwind of the infectious person. On the other hand, changes in particle settling velocity with relative humidity did not significantly affect the removal or travel distance of virus-laden scenarios.


title: Multitask Learning for Scalable and Dense Multilayer Bayesian Map Inference

abstract: This article presents a novel and flexible multitask multilayer Bayesian mapping framework with readily extendable attribute layers. The proposed framework goes beyond modern metric-semantic maps to provide even richer environmental information for robots in a single mapping formalism while exploiting intralayer and interlayer correlations. It removes the need for a robot to access and process information from many separate maps when performing a complex task, advancing the way robots interact with their environments. To this end, we design a multitask deep neural network with attention mechanisms as our front-end to provide heterogeneous observations for multiple map layers simultaneously. Our back-end runs a scalable closed-form Bayesian inference with only logarithmic time complexity. We apply the framework to build a dense robotic map including metric-semantic occupancy and traversability layers. Traversability ground truth labels are automatically generated from exteroceptive sensory data in a self-supervised manner. We present extensive experimental results on publicly available datasets and data collected by a 3D bipedal robot platform and show reliable mapping performance in different environments. Finally, we also discuss how the current framework can be extended to incorporate more information such as friction, signal strength, temperature, and physical quantity concentration using Gaussian map layers. The software for reproducing the presented results or running on customized data is made publicly available.


title: Derived homotopy algebras

abstract: We develop a theory of minimal models for algebras over an operad defined over a commutative ring, not necessarily a field, extending and supplementing the work of Sagave in the associative case.


title: Dynamics of flags over wide ranges of mass and bending stiffness

abstract: There have been many studies of the instability of a flexible plate or flag to flapping motions, and of large-amplitude flapping. Here we use inviscid simulations and a linearized model to study more generally how key quantities -- mode number (or wavenumber), frequency, and amplitude -- depend on the two dimensionless parameters, flag mass and bending stiffness. In the limit of small flag mass, flags perform traveling wave motions that move at nearly the speed of the oncoming flow. The flag mode number scales as the -1/4 power of bending stiffness. The flapping frequency has the same scaling, with an additional slight increase with flag mass in the small-mass regime. The flapping amplitude scales approximately as flag mass to the 1/2 power. For large flag mass, the dominant mode number is low (0 or 1), the flapping frequency tends to zero, and the amplitude saturates in the neighborhood of its upper limit (the flag length). In a linearized model, the fastest growing modes have somewhat different power law scalings for wavenumber and frequency. We discuss how the numerical scalings are consistent with a weakly nonlinear model.


title: Object Detection Based Handwriting Localization

abstract: We present an object detection based approach to localize handwritten regions from documents, which initially aims to enhance the anonymization during the data transmission. The concatenated fusion of original and preprocessed images containing both printed texts and handwritten notes or signatures are fed into the convolutional neural network, where the bounding boxes are learned to detect the handwriting. Afterwards, the handwritten regions can be processed (e.g. replaced with redacted signatures) to conceal the personally identifiable information (PII). This processing pipeline based on the deep learning network Cascade R-CNN works at 10 fps on a GPU during the inference, which ensures the enhanced anonymization with minimal computational overheads. Furthermore, the impressive generalizability has been empirically showcased: the trained model based on the English-dominant dataset works well on the fictitious unseen invoices, even in Chinese. The proposed approach is also expected to facilitate other tasks such as handwriting recognition and signature verification.


title: Committor Functions for Climate Phenomena at the Predictability Margin: The example of El Ni\~no Southern Oscillation in the Jin and Timmerman model

abstract: Many phenomena in the climate system lie in the gray zone between weather and climate: they are not amenable to deterministic forecast, but they still depend on the initial condition. A natural example is medium-range forecasting, which is inherently probabilistic because it lies beyond the deterministic predictability time of the atmosphere, but for which statistically significant prediction can be made which depend on the current state of the system. Similarly, one may ask the probability of occurrence of an El Ni\~no event several months ahead of time. In this paper, we introduce a quantity which corresponds precisely to this type of prediction problem: the committor function is the probability that an event takes place within a given time window, as a function of the initial condition. We explain the main mathematical properties of this probabilistic concept, and compute it in the case of a low-dimensional stochastic model for El-Ni\~no, the Jin and Timmerman model. In this context, we show that the ability to predict the probability of occurrence of the event of interest may differ strongly depending on the initial state. The main result is the new distinction between intrinsic probabilistic predictability (when the committor function is smooth and probability can be computed which does not depend sensitively on the initial condition) and intrinsic probabilistic unpredictability (when the committor function depends sensitively on the initial condition). We also demonstrate that the Jin and Timmerman model might be the first example of a stochastic differential equation with weak noise for which transition between attractors do not follow the Arrhenius law, which is expected based on large deviation theory and generic hypothesis.

categories: physics.ao-ph cond-mat.stat-mech nlin.CD physics.data-an


title: Discovery of two TNO-like bodies in the asteroid belt

abstract: Two extremely red main-belt asteroids: 203 Pompeja and 269 Justitia, were identified from combined visible and near-infrared spectroscopic observations collected at the IRTF and SAO observatories. These two asteroids have a redder spectral slope than any other D-type body, which are the reddest objects in the asteroid belt, and similar to RR and IR-class objects found in the outer Solar System among trans-Neptunian objects and Centaurs. Spectroscopic results suggest the presence of complex organic materials on the surface layer of these asteroids, implying that they could have formed in the vicinity of Neptune and been transplanted to the main belt region during a phase of planetary migration. 203 Pompeia is the only very red asteroid known so far among the ~250 bodies with diameter larger than 110 km (i.e. presumably structurally intact) found in the asteroid belt. These discoveries add another piece of evidence that the main asteroid belt hosts a population of bodies that were formed in the outskirt of the Solar System.


title: Connectedness percolation of fractal liquids

abstract: We apply connectedness percolation theory to fractal liquids of hard particles, and make use of a Percus-Yevick liquid state theory combined with a geometric connectivity criterion. We find that in fractal dimensions the percolation threshold interpolates continuously between integer-dimensional values, and that it decreases monotonically with increasing (fractal) dimension. The influence of hard-core interactions is only significant for dimensions below three. Finally, our theory incorrectly suggests that a percolation threshold is absent below about two dimensions, which we attribute to the breakdown of the connectedness Percus-Yevick closure.


title: Modularity in Reinforcement Learning via Algorithmic Independence in Credit Assignment

abstract: Many transfer problems require re-using previously optimal decisions for solving new tasks, which suggests the need for learning algorithms that can modify the mechanisms for choosing certain actions independently of those for choosing others. However, there is currently no formalism nor theory for how to achieve this kind of modular credit assignment. To answer this question, we define modular credit assignment as a constraint on minimizing the algorithmic mutual information among feedback signals for different decisions. We introduce what we call the modularity criterion for testing whether a learning algorithm satisfies this constraint by performing causal analysis on the algorithm itself. We generalize the recently proposed societal decision-making framework as a more granular formalism than the Markov decision process to prove that for decision sequences that do not contain cycles, certain single-step temporal difference action-value methods meet this criterion while all policy-gradient methods do not. Empirical evidence suggests that such action-value methods are more sample efficient than policy-gradient methods on transfer problems that require only sparse changes to a sequence of previously optimal decisions.


title: A Pilot Radio Search for Magnetic Activity in Directly Imaged Exoplanets

abstract: We present the first systematic search for GHz frequency radio emission from directly imaged exoplanets using Very Large Array (VLA) observations of sufficient angular resolution to separate the planets from their host stars. We obtained results for five systems and eight exoplanets located at $\lesssim 50$ pc, through new observations (Ross 458, GU Psc, and 51 Eri) and archival data (GJ 504 and HR 8799). We do not detect radio emission from any of the exoplanets, with $3\sigma$ luminosity upper limits of $(0.9-23)\times10^{21}$ erg s$^{-1}$. These limits are comparable to the level of radio emission detected in several ultracool dwarfs, including T dwarfs, whose masses are only a factor of two times higher than those of the directly-imaged exoplanets. Despite the lack of detections in this pilot study, we highlight the need for continued GHz frequency radio observations of nearby exoplanets at $\mu$Jy-level sensitivity.


title: Leveraging GPU batching for scalable nonlinear programming through massive Lagrangian decomposition

abstract: We present the implementation of a trust-region Newton algorithm ExaTron for bound-constrained nonlinear programming problems, fully running on multiple GPUs. Without data transfers between CPU and GPU, our implementation has achieved the elimination of a major performance bottleneck under a memory-bound situation, particularly when solving many small problems in batch. We discuss the design principles and implementation details for our kernel function and core operations. Different design choices are justified by numerical experiments. By using the application of distributed control of alternating current optimal power flow, where a large problem is decomposed into many smaller nonlinear programs using a Lagrangian approach, we demonstrate computational performance of ExaTron on the Summit supercomputer at Oak RidgeNational Laboratory. Our numerical results show the linear scaling with respect to the batch size and the number of GPUs and more than 35 times speedup on 6 GPUs than on 40 CPUs available on a single node.


title: Massey products for algebras over operads

abstract: We define a generalization of Massey products for algebras over a Koszul operad in characteristic zero, extending Massey's and Allday's and Retah's in the associative and Lie cases, respectively. We establish connections with minimal models and with Dimitrova's universal operadic cohomology class. We compute a Gerstenhaber algebra example and a hypercommutative algebra example related to the Chevalley-Eilenberg complex of the Heisenberg Lie algebra.


title: Sharp Lower Bounds on the Approximation Rate of Shallow Neural Networks

abstract: We consider the approximation rates of shallow neural networks with respect to the variation norm. Upper bounds on these rates have been established for sigmoidal and ReLU activation functions, but it has remained an important open problem whether these rates are sharp. In this article, we provide a solution to this problem by proving sharp lower bounds on the approximation rates for shallow neural networks, which are obtained by lower bounding the $L^2$-metric entropy of the convex hull of the neural network basis functions. In addition, our methods also give sharp lower bounds on the Kolmogorov $n$-widths of this convex hull, which show that the variation spaces corresponding to shallow neural networks cannot be efficiently approximated by linear methods. These lower bounds apply to both sigmoidal activation functions with bounded variation and to activation functions which are a power of the ReLU. Our results also quantify how much stronger the Barron spectral norm is than the variation norm and, combined with previous results, give the asymptotics of the $L^\infty$-metric entropy up to logarithmic factors in the case of the ReLU activation function.


title: Finite Element Approximations of a Class of Nonlinear Stochastic Wave Equation with Multiplicative Noise

abstract: Wave propagation problems have many applications in physics and engineering, and the stochastic effects are important in accurately modeling them due to the uncertainty of the media. This paper considers and analyzes a fully discrete finite element method for a class of nonlinear stochastic wave equations, where the diffusion term is globally Lipschitz continuous while the drift term is only assumed to satisfy weaker conditions as in [11]. The novelties of this paper are threefold. First, the error estimates cannot not be directly obtained if the numerical scheme in primal form is used. The numerical scheme in mixed form is introduced and several H\"{o}lder continuity results of the strong solution are proved, which are used to establish the error estimates in both $L^2$ norm and energy norms. Second, two types of discretization of the nonlinear term are proposed to establish the $L^2$ stability and energy stability results of the discrete solutions. These two types of discretization and proper test functions are designed to overcome the challenges arising from the stochastic scaling in time issues and the nonlinear interaction. These stability results play key roles in proving the probability of the set on which the error estimates hold approaches one. Third, higher order moment stability results of the discrete solutions are proved based on an energy argument and the underlying energy decaying property of the method. Numerical experiments are also presented to show the stability results of the discrete solutions and the convergence rates in various norms.


title: Test-Time Adaptation to Distribution Shift by Confidence Maximization and Input Transformation

abstract: Deep neural networks often exhibit poor performance on data that is unlikely under the train-time data distribution, for instance data affected by corruptions. Previous works demonstrate that test-time adaptation to data shift, for instance using entropy minimization, effectively improves performance on such shifted distributions. This paper focuses on the fully test-time adaptation setting, where only unlabeled data from the target distribution is required. This allows adapting arbitrary pretrained networks. Specifically, we propose a novel loss that improves test-time adaptation by addressing both premature convergence and instability of entropy minimization. This is achieved by replacing the entropy by a non-saturating surrogate and adding a diversity regularizer based on batch-wise entropy maximization that prevents convergence to trivial collapsed solutions. Moreover, we propose to prepend an input transformation module to the network that can partially undo test-time distribution shifts. Surprisingly, this preprocessing can be learned solely using the fully test-time adaptation loss in an end-to-end fashion without any target domain labels or source domain data. We show that our approach outperforms previous work in improving the robustness of publicly available pretrained image classifiers to common corruptions on such challenging benchmarks as ImageNet-C.


title: Optimal Convergence Rates for the Orthogonal Greedy Algorithm

abstract: We analyze the orthogonal greedy algorithm when applied to dictionaries $\mathbb{D}$ whose convex hull has small entropy. We show that if the metric entropy of the convex hull of $\mathbb{D}$ decays at a rate of $O(n^{-\frac{1}{2}-\alpha})$ for $\alpha > 0$, then the orthogonal greedy algorithm converges at the same rate on the variation space of $\mathbb{D}$. This improves upon the well-known $O(n^{-\frac{1}{2}})$ convergence rate of the orthogonal greedy algorithm in many cases, most notably for dictionaries corresponding to shallow neural networks. These results hold under no additional assumptions on the dictionary beyond the decay rate of the entropy of its convex hull. In addition, they are robust to noise in the target function and can be extended to convergence rates on the interpolation spaces of the variation norm. We show empirically that the predicted rates are obtained for the dictionary corresponding to shallow neural networks with Heaviside activation function in two dimensions. Finally, we show that these improved rates are sharp and prove a negative result showing that the iterates generated by the orthogonal greedy algorithm cannot in general be bounded in the variation norm of $\mathbb{D}$.


title: Generalized cohomology theories for algebraic stacks

abstract: We extend the stable motivic homotopy category of Voevodsky to the class of scalloped algebraic stacks, and show that it admits the formalism of Grothendieck's six operations. Objects in this category represent generalized cohomology theories for stacks like algebraic K-theory, as well as new examples like genuine motivic cohomology and algebraic cobordism. These cohomology theories admit Gysin maps and satisfy homotopy invariance, localization, and Mayer-Vietoris. We also prove a fixed point localization formula for torus actions. Finally, the construction is contrasted with a "lisse-extended" stable motivic homotopy category, defined for arbitrary stacks: we show for example that lisse-extended motivic cohomology of quotient stacks is computed by the equivariant higher Chow groups of Edidin-Graham, and we also get a good new theory of Borel-equivariant algebraic cobordism. Moreover, the lisse-extended motivic homotopy type is shown to recover all previous constructions of motives of stacks.


title: Characterization of the Variation Spaces Corresponding to Shallow Neural Networks

abstract: We study the variation space corresponding to a dictionary of functions in $L^2(\Omega)$ for a bounded domain $\Omega\subset \mathbb{R}^d$. Specifically, we compare the variation space, which is defined in terms of a convex hull with related notions based on integral representations. This allows us to show that three important notions relating to the approximation theory of shallow neural networks, the Barron space, the spectral Barron space, and the Radon BV space, are actually variation spaces with respect to certain natural dictionaries.


title: A Note on the Topology of the First Stage of 2SLS with Many Instruments

abstract: The finite sample properties of estimators are usually understood or approximated using asymptotic theories. Two main asymptotic constructions have been used to characterize the presence of many instruments. The first assumes that the number of instruments increases with the sample size. I demonstrate that in this case, one of the key assumptions used in the asymptotic construction may imply that the number of ``effective" instruments should be finite, resulting in an internal contradiction. The second asymptotic representation considers that the number of instrumental variables (IVs) may be finite, infinite, or even a continuum. The number does not change with the sample size. In this scenario, the regularized estimator obtained depends on the topology imposed on the set of instruments as well as on a regularization parameter. These restrictions may induce a bias or restrict the set of admissible instruments. However, the assumptions are internally coherent. The limitations of many IVs asymptotic assumptions provide support for finite sample distributional studies to better understand the behavior of many IV estimators.


title: Multimodal Trajectory Prediction Conditioned on Lane-Graph Traversals

abstract: Accurately predicting the future motion of surrounding vehicles requires reasoning about the inherent uncertainty in driving behavior. This uncertainty can be loosely decoupled into lateral (e.g., keeping lane, turning) and longitudinal (e.g., accelerating, braking). We present a novel method that combines learned discrete policy rollouts with a focused decoder on subsets of the lane graph. The policy rollouts explore different goals given current observations, ensuring that the model captures lateral variability. Longitudinal variability is captured by our latent variable model decoder that is conditioned on various subsets of the lane graph. Our model achieves state-of-the-art performance on the nuScenes motion prediction dataset, and qualitatively demonstrates excellent scene compliance. Detailed ablations highlight the importance of the policy rollouts and the decoder architecture.


title: Untidy Data: The Unreasonable Effectiveness of Tables

abstract: Working with data in table form is usually considered a preparatory and tedious step in the sensemaking pipeline; a way of getting the data ready for more sophisticated visualization and analytical tools. But for many people, spreadsheets -- the quintessential table tool -- remain a critical part of their information ecosystem, allowing them to interact with their data in ways that are hidden or abstracted in more complex tools. This is particularly true for data workers: people who work with data as part of their job but do not identify as professional analysts or data scientists. We report on a qualitative study of how these workers interact with and reason about their data. Our findings show that data tables serve a broader purpose beyond data cleanup at the initial stage of a linear analytic flow: users want to see and "get their hands on" the underlying data throughout the analytics process, reshaping and augmenting it to support sensemaking. They reorganize, mark up, layer on levels of detail, and spawn alternatives within the context of the base data. These direct interactions and human-readable table representations form a rich and cognitively important part of building understanding of what the data mean and what they can do with it. We argue that interactive tables are an important visualization idiom in their own right; that the direct data interaction they afford offers a fertile design space for visual analytics; and that sense making can be enriched by more flexible human-data interaction than is currently supported in visual analytics tools.


title: Hamilton-Jacobi Equations for Two Classes of State-Constrained Zero-Sum Games

abstract: This paper presents Hamilton-Jacobi (HJ) formulations for two classes of two-player zero-sum games: one with a maximum cost value over time, and one with a minimum cost value over time. In the zero-sum game setting, player A minimizes the given cost while satisfying state constraints, and player B wants to prevent player A's success. For each class of problems, this paper presents two HJ equations: one for time-varying dynamics, cost, and state constraint; the other for time-invariant dynamics, cost, and state constraint. Utilizing the HJ equations, the optimal control for each player is analyzed, and a numerical algorithm is presented to compute the solution to the HJ equations. A two-dimensional water system is introduced as an example to demonstrate the proposed HJ framework.


title: An Uncertainty Estimation Framework for Probabilistic Object Detection

abstract: In this paper, we introduce a new technique that combines two popular methods to estimate uncertainty in object detection. Quantifying uncertainty is critical in real-world robotic applications. Traditional detection models can be ambiguous even when they provide a high-probability output. Robot actions based on high-confidence, yet unreliable predictions, may result in serious repercussions. Our framework employs deep ensembles and Monte Carlo dropout for approximating predictive uncertainty, and it improves upon the uncertainty estimation quality of the baseline method. The proposed approach is evaluated on publicly available synthetic image datasets captured from sequences of video.


title: Comaximal Factorization Lattices

abstract: Brewer and Heinzer studied the (integral) domains D having the property that each proper ideal A of D has a comaximal ideal factorization with some additional property. They proved that for a domain D, the following are equivalent: (1) Each proper ideal A of D has a comaximal factorization where the factors have prime radical (resp. are primary, resp. are prime powers). (2) The prime spectrum of D is a tree under inclusion and each ideal of D has only finitely many minimal primes (resp. D is one dimensional and each ideal of D has only finitely many minimal primes, resp. D is a Dedekind domain). The aim of this paper is to show that most of the results can be obtained in the setup of multiplicative lattices.


title: Understanding Cognitive Fatigue from fMRI Scans with Self-supervised Learning

abstract: Functional magnetic resonance imaging (fMRI) is a neuroimaging technique that records neural activations in the brain by capturing the blood oxygen level in different regions based on the task performed by a subject. Given fMRI data, the problem of predicting the state of cognitive fatigue in a person has not been investigated to its full extent. This paper proposes tackling this issue as a multi-class classification problem by dividing the state of cognitive fatigue into six different levels, ranging from no-fatigue to extreme fatigue conditions. We built a spatio-temporal model that uses convolutional neural networks (CNN) for spatial feature extraction and a long short-term memory (LSTM) network for temporal modeling of 4D fMRI scans. We also applied a self-supervised method called MoCo (Momentum Contrast) to pre-train our model on a public dataset BOLD5000 and fine-tuned it on our labeled dataset to predict cognitive fatigue. Our novel dataset contains fMRI scans from Traumatic Brain Injury (TBI) patients and healthy controls (HCs) while performing a series of N-back cognitive tasks. This method establishes a state-of-the-art technique to analyze cognitive fatigue from fMRI data and beats previous approaches to solve this problem.


title: Penrose process for a charged black hole in a uniform magnetic field

abstract: Spinning black holes create electromagnetic storms when immersed in ambient magnetic fields, illuminating the otherwise epically dark terrain. In an electromagnetic extension of the Penrose process, tremendous energy can be extracted, boosting the energy of radiating particles far more efficiently than the mechanical Penrose process. We locate the regions from which energy can be mined and demonstrate explicitly that they are no longer restricted to the ergosphere. We also show that there can be toroidal regions that trap negative energy particles in orbit around the black hole. We find that the effective charge coupling between the black hole and the super-radiant particles decreases as energy is extracted, much like the spin of a black hole decreases in the mechanical analogue. While the effective coupling decreases, the actual charge of the black hole increases in magnitude reaching the energetically-favored Wald value, at which point energy extraction is impeded. We demonstrate the array of orbits for products from the electromagnetic Penrose process.


title: Are conditional GANs explicitly conditional?

abstract: This paper proposes two important contributions for conditional Generative Adversarial Networks (cGANs) to improve the wide variety of applications that exploit this architecture. The first main contribution is an analysis of cGANs to show that they are not explicitly conditional. In particular, it will be shown that the discriminator and subsequently the cGAN does not automatically learn the conditionality between inputs. The second contribution is a new method, called a contrario cGAN, that explicitly models conditionality for both parts of the adversarial architecture via a novel a contrario loss that involves training the discriminator to learn unconditional (adverse) examples. This leads to a novel type of data augmentation approach for GANs (a contrario learning) which allows to restrict the search space of the generator to conditional outputs using adverse examples. Extensive experimentation is carried out to evaluate the conditionality of the discriminator by proposing a probability distribution analysis. Comparisons with the cGAN architecture for different applications show significant improvements in performance on well known datasets including, semantic image synthesis, image segmentation, monocular depth prediction and "single label"-to-image using different metrics including Fr\'echet Inception Distance (FID), mean Intersection over Union (mIoU), Root Mean Square Error log (RMSE log) and Number of statistically-Different Bins (NDB).


title: Level-rank duality of knot and link invariants

abstract: A number of results for the level-rank duality of $G(N)_K$ $\leftrightarrow$ $G(K)_N$ Chern-Simons theory are summarized, with emphasis on the applications to knot and link invariants. Explicit examples for $SU(2)_K$ $\leftrightarrow$ $SU(K)_2$ illustrate general results. A criterion to distinguish torus knots and links from hyperbolic knots and links, based on tables constructed by Kaul for one and two strand invariants, is presented. Possible symmetries of hyperbolic knot and link invariants are discussed. The level-rank duality of torus knot and link invariants of minimal models is examined


title: Small random initialization is akin to spectral learning: Optimization and generalization guarantees for overparameterized low-rank matrix reconstruction

abstract: Recently there has been significant theoretical progress on understanding the convergence and generalization of gradient-based methods on nonconvex losses with overparameterized models. Nevertheless, many aspects of optimization and generalization and in particular the critical role of small random initialization are not fully understood. In this paper, we take a step towards demystifying this role by proving that small random initialization followed by a few iterations of gradient descent behaves akin to popular spectral methods. We also show that this implicit spectral bias from small random initialization, which is provably more prominent for overparameterized models, also puts the gradient descent iterations on a particular trajectory towards solutions that are not only globally optimal but also generalize well. Concretely, we focus on the problem of reconstructing a low-rank matrix from a few measurements via a natural nonconvex formulation. In this setting, we show that the trajectory of the gradient descent iterations from small random initialization can be approximately decomposed into three phases: (I) a spectral or alignment phase where we show that that the iterates have an implicit spectral bias akin to spectral initialization allowing us to show that at the end of this phase the column space of the iterates and the underlying low-rank matrix are sufficiently aligned, (II) a saddle avoidance/refinement phase where we show that the trajectory of the gradient iterates moves away from certain degenerate saddle points, and (III) a local refinement phase where we show that after avoiding the saddles the iterates converge quickly to the underlying low-rank matrix. Underlying our analysis are insights for the analysis of overparameterized nonconvex optimization schemes that may have implications for computational problems beyond low-rank reconstruction.

categories: cs.LG cs.IT math.IT math.OC math.ST stat.ML stat.TH


title: Ray Velocity Derivatives in Anisotropic Elastic Media. Part I -- General Anisotropy

abstract: We present an original, generic, and efficient approach for computing the first and second partial derivatives of ray velocities along ray paths in general anisotropic elastic media. These derivatives are used in solving kinematic problems, like two-point ray bending methods and seismic tomography, and they are essential for evaluating the dynamic properties along the rays (amplitudes and phases). The traveltime is delivered through an integral over a given Lagrangian defined at each point along the ray. Although the Lagrangian cannot be explicitly expressed in terms of the medium properties and the ray direction components, its derivatives can still be formulated analytically using the corresponding arclength-related Hamiltonian that can be explicitly expressed in terms of the medium properties and the slowness vector components. This requires first to invert for the slowness vector components, given the ray direction components. Computation of the slowness vector and the ray velocity derivatives is considerably simplified by using an auxiliary scaled-time-related Hamiltonian obtained directly from the Christoffel equation and connected to the arclength-related Hamiltonian by a simple scale factor. This study consists of two parts. In Part I, we consider general anisotropic (triclinic) media, and provide the derivatives (gradients and Hessians) of the ray velocity, with respect to (1) the spatial/directional vectors and (2) the elastic model parameters. In Part II, we apply the theory of Part I explicitly to polar anisotropic media (transverse isotropy with tilted axis of symmetry, TTI), and obtain the explicit ray velocity derivatives for the coupled qP and qSV waves and for SH waves.


title: A Survey on Trust Metrics for Autonomous Robotic Systems

abstract: This paper surveys the area of Trust Metrics related to security for autonomous robotic systems. As the robotics industry undergoes a transformation from programmed, task oriented, systems to Artificial Intelligence-enabled learning, these autonomous systems become vulnerable to several security risks, making a security assessment of these systems of critical importance. Therefore, our focus is on a holistic approach for assessing system trust which requires incorporating system, hardware, software, cognitive robustness, and supplier level trust metrics into a unified model of trust. We set out to determine if there were already trust metrics that defined such a holistic system approach. While there are extensive writings related to various aspects of robotic systems such as, risk management, safety, security assurance and so on, each source only covered subsets of an overall system and did not consistently incorporate the relevant costs in their metrics. This paper attempts to put this prior work into perspective, and to show how it might be extended to develop useful system-level trust metrics for evaluating complex robotic (and other) systems.


title: Pytearcat: PYthon TEnsor AlgebRa calCulATor

abstract: This paper introduces the first release of Pytearcat, a Python package developed to compute tensor algebra operations in the context of theoretical physics, for instance, in general relativity. Given that working with tensors can become a complex task, people often rely on computational tools to perform tensor calculations. We aim to build a tensor calculator based on Python, which benefits from being free and easy to use. Pytearcat syntax resembles the usual physics notation for tensor calculus, such as the Einstein notation for index contraction. This version allows the user to perform many tensor operations, including derivatives and series expansions, along with routines to obtain the typical General Relativity tensors. A particular concern was put in the execution times, leading to incorporate an alternative core for the symbolic calculations, enabling to reach much faster execution times. The syntax and the versatility of Pytearcat are the most important features of this package, where the latter can be used to extend Pytearcat to other areas of theoretical physics.


title: Early Mobility Recognition for Intensive Care Unit Patients Using Accelerometers

abstract: With the development of the Internet of Things(IoT) and Artificial Intelligence(AI) technologies, human activity recognition has enabled various applications, such as smart homes and assisted living. In this paper, we target a new healthcare application of human activity recognition, early mobility recognition for Intensive Care Unit(ICU) patients. Early mobility is essential for ICU patients who suffer from long-time immobilization. Our system includes accelerometer-based data collection from ICU patients and an AI model to recognize patients' early mobility. To improve the model accuracy and stability, we identify features that are insensitive to sensor orientations and propose a segment voting process that leverages a majority voting strategy to recognize each segment's activity. Our results show that our system improves model accuracy from 77.78\% to 81.86\% and reduces the model instability (standard deviation) from 16.69\% to 6.92\%, compared to the same AI model without our feature engineering and segment voting process.


title: Representing polynomial of ST-CONNECTIVITY

abstract: We show that the coefficients of the representing polynomial of any monotone Boolean function are the values of the M\"obius function of an atomistic lattice related to this function. Using this we determine the representing polynomial of any Boolean function corresponding to a ST-CONNECTIVITY problem in acyclic quivers (directed acyclic multigraphs). Only monomials corresponding to unions of paths have non-zero coefficients which are $(-1)^D$ where $D$ is an easily computable function of the quiver corresponding to the monomial (it is the number of plane regions in the case of planar graphs). We determine that the number of monomials with non-zero coefficients for the two-dimensional $n \times n$ grid connectivity problem is $2^{\Omega(n^2)}$.


title: On a toroidalization for klt singularities

abstract: In this article, we prove a toroidalization principle for finite actions on klt singularities. As an application, we prove that the Jordan property for the regional fundamental group of klt singularities can be realized geometrically: by extracting a toric singularity over the klt germ. In the course of the proof, we will prove statements about finite actions on dual complexes and almost fixed points in the fibers of equivariant Fano type morphisms. Furthermore, we will prove that the rank of a fundamental group of the klt singularity is bounded above by its regularity.


title: On the potential of sequential and non-sequential regression models for Sentinel-1-based biomass prediction in Tanzanian miombo forests

abstract: This study derives regression models for above-ground biomass (AGB) estimation in miombo woodlands of Tanzania that utilise the high availability and low cost of Sentinel-1 data. The limited forest canopy penetration of C-band SAR sensors along with the sparseness of available ground truth restrict their usefulness in traditional AGB regression models. Therefore, we propose to use AGB predictions based on airborne laser scanning (ALS) data as a surrogate response variable for SAR data. This dramatically increases the available training data and opens for flexible regression models that capture fine-scale AGB dynamics. This becomes a sequential modelling approach, where the first regression stage has linked in situ data to ALS data and produced the AGB prediction map; We perform the subsequent stage, where this map is related to Sentinel-1 data. We develop a traditional, parametric regression model and alternative non-parametric models for this stage. The latter uses a conditional generative adversarial network (cGAN) to translate Sentinel-1 images into ALS-based AGB prediction maps. The convolution filters in the neural networks make them contextual. We compare the sequential models to traditional, non-sequential regression models, all trained on limited AGB ground reference data. Results show that our newly proposed non-sequential Sentinel-1-based regression model performs better quantitatively than the sequential models, but achieves less sensitivity to fine-scale AGB dynamics. The contextual cGAN-based sequential models best reproduce the distribution of ALS-based AGB predictions. They also reach a lower RMSE against in situ AGB data than the parametric sequential model, indicating a potential for further development.


title: How to Reach Real-Time AI on Consumer Devices? Solutions for Programmable and Custom Architectures

abstract: The unprecedented performance of deep neural networks (DNNs) has led to large strides in various Artificial Intelligence (AI) inference tasks, such as object and speech recognition. Nevertheless, deploying such AI models across commodity devices faces significant challenges: large computational cost, multiple performance objectives, hardware heterogeneity and a common need for high accuracy, together pose critical problems to the deployment of DNNs across the various embedded and mobile devices in the wild. As such, we have yet to witness the mainstream usage of state-of-the-art deep learning algorithms across consumer devices. In this paper, we provide preliminary answers to this potentially game-changing question by presenting an array of design techniques for efficient AI systems. We start by examining the major roadblocks when targeting both programmable processors and custom accelerators. Then, we present diverse methods for achieving real-time performance following a cross-stack approach. These span model-, system- and hardware-level techniques, and their combination. Our findings provide illustrative examples of AI systems that do not overburden mobile hardware, while also indicating how they can improve inference accuracy. Moreover, we showcase how custom ASIC- and FPGA-based accelerators can be an enabling factor for next-generation AI applications, such as multi-DNN systems. Collectively, these results highlight the critical need for further exploration as to how the various cross-stack solutions can be best combined in order to bring the latest advances in deep learning close to users, in a robust and efficient manner.


title: Towards a theory of coarse geometry of operator spaces

abstract: We introduce two notions of coarse embeddability between operator spaces: almost complete coarse embeddability of bounded subsets and spherically-complete coarse embeddability. We provide examples showing that these notions are strictly weaker than complete isomorphic embeddability - in fact, they do not even imply isomorphic embeddability. On the other hand, we show that, despite their nonlinearity, the existence of such embeddings provides restrictions on the linear operator space structures of the spaces. Examples of nonlinear equivalences between operator spaces are also provided.


title: Evading Adversarial Example Detection Defenses with Orthogonal Projected Gradient Descent

abstract: Evading adversarial example detection defenses requires finding adversarial examples that must simultaneously (a) be misclassified by the model and (b) be detected as non-adversarial. We find that existing attacks that attempt to satisfy multiple simultaneous constraints often over-optimize against one constraint at the cost of satisfying another. We introduce Orthogonal Projected Gradient Descent, an improved attack technique to generate adversarial examples that avoids this problem by orthogonalizing the gradients when running standard gradient-based attacks. We use our technique to evade four state-of-the-art detection defenses, reducing their accuracy to 0% while maintaining a 0% detection rate.


title: Birkhoff Averages and the Breakdown of Invariant Tori in Volume-Preserving Maps

abstract: In this paper, we develop numerical methods based on the weighted Birkhoff average for studying two-dimensional invariant tori for volume-preserving maps. The methods do not rely on symmetries, such as time-reversal symmetry, nor on approximating tori by periodic orbits. The rate of convergence of the average gives a sharp distinction between chaotic and regular dynamics and allows accurate computation of rotation vectors for regular orbits. Resonant and rotational tori are distinguished by computing the resonance order of the rotation vector to a given precision. Critical parameter values, where tori are destroyed, are computed by a sharp decrease in convergence rate of the Birkhoff average. We apply these methods for a three-dimensional generalization of Chirikov's standard map: an angle-action map with two angle variables. Computations on grids in frequency and perturbation amplitude allow estimates of the critical set. We also use continuation to follow tori with fixed rotation vectors. We test three conjectures for cubic fields that have been proposed to give locally robust invariant tori.


title: A Low Power Threshold, Ultrathin Optical Limiter Based on a Nonlinear Zone Plate

abstract: Ultrathin optical limiters are needed to protect light sensitive components in miniaturized optical systems. However, it has proven challenging to achieve a sufficiently low optical limiting threshold. In this work, we theoretically show that an ultrathin optical limiter with low threshold intensity can be realized using a nonlinear zone plate. The zone plate is embedded with nonlinear saturable absorbing materials that allow the device to focus low intensity light, while high intensity light is transmitted as a plane wave without a focal spot. Based on this proposed mechanism, we use the finite-difference time-domain method to computationally design a zone plate embedded with InAs quantum dots as the saturable absorbing material. The device has a thickness of just 0.5 $\mu m$ and exhibits good optical limiting behavior with a threshold intensity as low as 0.45 kW/$cm^2$, which is several orders of magnitude lower than current ultrathin flat-optics-based optical limiters. This design can be optimized for different operating wavelengths and threshold intensities by using different saturable absorbing materials. Additionally, the diameter and focal length of the nonlinear zone plate can be easily adjusted to fit different systems and applications. Due to its flexible design, low power threshold, and ultrathin thickness, this optical limiting concept may be promising for application in miniaturized optical systems.


title: Nonlinear dynamics and phase space transport by chorus emission

abstract: Chorus emission in planetary magnetospheres is taken as working paradigm to motivate a short tutorial trip through theoretical plasma physics methods and their applications. Starting from basic linear theory, readers are first made comfortable with whistler wave packets and their propagation in slowly varying weakly nonuniform media, such as the Earth's magnetosphere, where they can be amplified by a population of supra-thermal electrons. The nonlinear dynamic description of energetic electrons in the phase space in the presence of self-consistently evolving whistler fluctuation spectrum is progressively introduced by addressing renormalization of the electron response and spectrum evolution equations. Analytical and numerical results on chorus frequency chirping are obtained and compared with existing observations and particle in cell simulations. Finally, the general theoretical framework constructed during this short trip through chorus physics is used to draw analogies with condensed matter and laser physics as well as magnetic confinement fusion research. Discussing these analogies ultimately presents plasma physics as an exciting cross-disciplinary field to study.


title: A Phase Space Approach to the Conformal Construction of Non-Vacuum Initial Data Sets in General Relativity

abstract: We present a uniform (and unambiguous) procedure for scaling the matter fields in implementing the conformal method to parameterize and construct solutions of Einstein constraint equations with coupled matter sources. The approach is based on a phase space representation of the space-time matter fields after a careful $n+1$ decomposition into spatial fields $B$ and conjugate momenta $\Pi_B$, which are specified directly and are conformally invariant quantities. We show that if the Einstein-matter field theory is specified by a Lagrangian which is diffeomorphism invariant and involves no dependence on derivatives of the space-time metric in the matter portion of the Lagrangian, then fixing $B$ and $\Pi_B$ results in conformal constraint equations that, for constant-mean curvature initial data, semi-decouple just as they do for the vacuum Einstein conformal constraint equations. We prove this result by establishing a structural property of the Einstein momentum constraint that is independent of the conformal method: For an Einstein-matter field theory which satisfies the conditions just stated, if $B$ and $\Pi_B$ satisfy the matter Euler-Lagrange equations, then (in suitable form) the right-hand side of the momentum constraint on each spatial slice depends only on $B$ and $\Pi_B$ and is independent of the space-time metric. We discuss the details of our construction in the special cases of the following models: Einstein-Maxwell-charged scalar field, Einstein-Proca, Einstein-perfect fluid, and Einstein-Maxwell-charged dust. In these examples we find that our technique gives a theoretical basis for scaling rules, such as those for electromagnetism, that have worked pragmatically in the past, but also generates new equations with advantageous features for perfect fluids that allow direct specification of total rest mass and total charge in any spatial region.


title: BO Ceti: Dwarf Nova Showing Both IW And and SU UMa-Type Features

abstract: IW And stars are a recently recognized subgroup of dwarf novae which are characterized by (often repetitive) slowly rising standstills terminated by brightening, but the exact mechanism for this variation is not yet identified. We have identified BO Cet, which had been considered as a novalike cataclysmic variable, as a new member of IW And stars based on the behavior in 2019-2020. In addition to this, the object showed dwarf nova-type outbursts in 2020-2021, and superhumps having a period 7.8% longer than the orbital one developed at least during one long outburst. This object has been confirmed as an SU UMa-type dwarf nova with an exceptionally long orbital period (0.1398 d). BO Cet is thus the first cataclysmic variable showing both SU UMa-type and IW And-type features. We obtained a mass ratio (q) of 0.31-0.34 from the superhumps in the growing phase (stage A superhumps). At this q, the radius of the 3:1 resonance, responsible for tidal instability and superhumps, and the tidal truncation radius are very similar. We interpret that in some occasions this object showed IW And-type variation when the disk size was not large enough, but that the radius of the 3:1 resonance could be reached as the result of thermal instability. We also discuss that there are SU UMa-type dwarf novae above q=0.30, which is above the previously considered limit (q~0.25) derived from numerical simulations and that this is possible since the radius of the 3:1 resonance is inside the tidal truncation radius. We constrained the mass of the white dwarf larger than 1.0Msol, which may be responsible for the IW And-type behavior and the observed strength of the He II emission. The exact reason, however, why this object is unique in that it shows both SU UMa-type and IW And-type features is still unsolved.


title: Multi-Sensor Fusion based Robust Row Following for Compact Agricultural Robots

abstract: This paper presents a state-of-the-art LiDAR based autonomous navigation system for under-canopy agricultural robots. Under-canopy agricultural navigation has been a challenging problem because GNSS and other positioning sensors are prone to significant errors due to attentuation and multi-path caused by crop leaves and stems. Reactive navigation by detecting crop rows using LiDAR measurements is a better alternative to GPS but suffers from challenges due to occlusion from leaves under the canopy. Our system addresses this challenge by fusing IMU and LiDAR measurements using an Extended Kalman Filter framework on low-cost hardwware. In addition, a local goal generator is introduced to provide locally optimal reference trajectories to the onboard controller. Our system is validated extensively in real-world field environments over a distance of 50.88~km on multiple robots in different field conditions across different locations. We report state-of-the-art distance between intervention results, showing that our system is able to safely navigate without interventions for 386.9~m on average in fields without significant gaps in the crop rows, 56.1~m in production fields and 47.5~m in fields with gaps (space of 1~m without plants in both sides of the row).


title: Fast Radio Bursts as crustal dynamical events induced by magnetic field evolution in young magnetars

abstract: We revisit in this work a model for repeating Fast Radio Bursts based of the release of energy provoked by the magnetic field dynamics affecting a magnetar's crust. We address the basics of such a model by solving the propagation of the perturbation approximately, and quantify the energetics and the radiation by bunches of charges in the so-called {\it charge starved} region in the magnetosphere. The (almost) simultaneous emission of newly detected X-rays from SGR 1935+2154 is tentatively associated to a reconnection behind the propagation. The strength of $f$-mode gravitational radiation excited by the event is quantified, and more detailed studies of the non-linear (spiky) soliton solutions suggested.


title: Regression analysis of distributional data through Multi-Marginal Optimal transport

abstract: We formulate and solve a regression problem with time-stamped distributional data. Distributions are considered as points in the Wasserstein space of probability measures, metrized by the 2-Wasserstein metric, and may represent images, power spectra, point clouds of particles, and so on. The regression seeks a curve in the Wasserstein space that passes closest to the dataset. Our regression problem allows utilizing general curves in a Euclidean setting (linear, quadratic, sinusoidal, and so on), lifted to corresponding measure-valued curves in the Wasserstein space. It can be cast as a multi-marginal optimal transport problem that allows efficient computation. Illustrative academic examples are presented.


title: Magnetic Resonance Force Microscopy with Overlapping Frequencies of Cantilever and Spin

abstract: We have studied theoretically magnetic resonance force microscopy (MRFM) with a high frequency nanomechanical cantilever when the cantilever frequency matches the resonant frequency of a single electron spin. Our estimations show that in this scenario the relative frequency shift of the cantilever can be much greater than the record MRFM frequency shift achieved in experiments with a single spin detection. Experimental realization of our proposal could open the way for fast detection of a single electron spin and even for detection of a single nuclear spin.


title: Quasiparticle scattering in a superconductor near a nematic critical point: resonance mode and multiple attractive channels

abstract: We analyze the scattering rate for 2D fermions interacting via soft nematic fluctuations. The ground state is an s-wave superconductor, but other pairing channels are almost equally attractive. This strongly alters the scattering rate: At energies beyond the pairing gap $\Delta$, it is renormalized by contributions from all pairing channels. At energies of order $\Delta$, it is determined by the competition between scattering into a gapped continuum and dispersing nematic resonance. The outcome is a "peak-peak-dip-hump" spectrum, similar, but not identical, to the "peak-dip-hump" structure in the cuprates.


title: Approximation Schemes for Capacitated Vehicle Routing on Graphs of Bounded Treewidth, Bounded Doubling, or Highway Dimension

abstract: In this paper, we present Approximation Schemes for Capacitated Vehicle Routing Problem (CVRP) on several classes of graphs. In CVRP, introduced by Dantzig and Ramser (1959), we are given a graph $G=(V,E)$ with metric edges costs, a depot $r\in V$, and a vehicle of bounded capacity $Q$. The goal is to find minimum cost collection of tours for the vehicle that returns to the depot, each visiting at most $Q$ nodes, such that they cover all the nodes. This generalizes classic TSP and has been studied extensively. In the more general setting, each node $v$ has a demand $d_v$ and the total demand of each tour must be no more than $Q$. Either the demand of each node must be served by one tour (unsplittable) or can be served by multiple tour (splittable). The best known approximation algorithm for general graphs has ratio $\alpha+2(1-\epsilon)$ (for the unsplittable) and $\alpha+1-\epsilon$ (for the splittable) for some fixed $\epsilon>\frac{1}{3000}$, where $\alpha$ is the best approximation for TSP. Even for the case of trees, the best approximation ratio is $4/3$ by Becker (2018) and it has been an open question if there is an approximation scheme for this simple class of graphs. Das and Mathieu (2015) presented an approximation scheme with time $n^{\log^{O(1/\epsilon)}n}$ for Euclidean plane $\mathbb{R}^2$. No other approximation scheme is known for any other class of metrics (without further restrictions on $Q$). In this paper, we make significant progress on this classic problem by presenting Quasi-Polynomial Time Approximation Schemes (QPTAS) for graphs of bounded treewidth, graphs of bounded highway dimensions, and graphs of bounded doubling dimensions. For comparison, our result implies an approximation scheme for Euclidean plane with run time $n^{O(\log^{10}n/\epsilon^{9})}$.


title: Empirical Framework for Cournot Oligopoly with Private Information

abstract: We propose an empirical framework for asymmetric Cournot oligopoly with private information about variable costs. First, considering a linear demand for a homogenous product with a random intercept, we characterize the Bayesian Cournot-Nash equilibrium. Then we establish the identification of the joint distribution of demand and firm-specific cost distributions. Following the identification steps, we propose a likelihood-based estimation method and apply it to the global market for crude-oil and quantify the welfare effect of private information. We also consider extensions of the model to include either product differentiation, conduct parameters, nonlinear demand, or selective entry.


title: The first photometric and spectroscopic analysis of the extremely low mass ratio contact binary NSVS 5029961

abstract: We performed photometric and spectroscopic investigations of NSVS 5029961 for the first time. The new BV(RI)$_c$-band light curves were obtained with the 1.0-m telescope at Weihai Observatory of Shandong University. Applying the Wilson-Devinney program, we found that NSVS 5029961 is an A-subtype shallow contact binary with extremely low mass ratio (q = 0.1515, f = 19.1\%). Six spectra have been obtained by LAMOST, and many chromospheric activity emission line indicators were detected in the spectra, revealing that the target exhibits strong chromospheric activity. We calculated the absolute parameters with the photometric solutions and Gaia distance, and estimated the initial masses of the two components and the age of the binary. The evolutionary status was discussed by using the mass-radius and mass-luminosity diagrams. The result shows the primary component is a little evolved star and the secondary component has evolved away from the main sequence. The formation and evolution investigations of NSVS 5029661 indicate that it may have evolved from a detached binary with short period and low mass ratio by angular momentum loss via magnetic braking and case A mass transfer, and is in a stable contact stage at present.


title: Directional asymptotics of Fej\'er monotone sequences

abstract: The notion of Fej\'er monotonicity is instrumental in unifying the convergence proofs of many iterative methods, such as the Krasnoselskii-Mann iteration, the proximal point method, the Douglas-Rachford splitting algorithm, and many others. In this paper, we present directionally asymptotical results of strongly convergent subsequences of Fej\'er monotone sequences. We also provide examples to show that the sets of directionally asymptotic cluster points can be large and that weak convergence is needed in infinite-dimensional spaces.


title: On the arithmetic Siegel--Weil formula for GSpin Shimura varieties

abstract: We formulate and prove a local arithmetic Siegel--Weil formula for GSpin Rapoport--Zink spaces, which is a precise identity between the arithmetic intersection numbers of special cycles on GSpin Rapoport--Zink spaces and the derivatives of local representation densities of quadratic forms. As a first application, we prove a semi-global arithmetic Siegel--Weil formula as conjectured by Kudla, which relates the arithmetic intersection numbers of special cycles on GSpin Shimura varieties at a place of good reduction and the central derivatives of nonsingular Fourier coefficients of incoherent Siegel Eisenstein series.


title: Breakup of the Synchronous State of Binary Asteroid Systems

abstract: This paper continues the authors' previous work and presents a coplanar averaged ellipsoid-ellipsoid model of synchronous binary asteroid system (BAS) plus thermal and tidal effects. Using this model, we analyze the breakup mechanism of the synchronous BAS. Different from the classical spin-orbit coupling model which neglects the rotational motion's influence on the orbital motion, our model considers simultaneously the orbital motion and the rotational motions. Our findings are following. (1) Stable region of the secondary's synchronous state is mainly up to the secondary's shape. The primary's shape has little influence on it. (2) The stable region shrinks continuously with the increasing value of the secondary's shape parameter $a_B/b_B$. Beyond the value of $a_B/b_B=\sqrt{2}$, the planar stable region for the secondary's synchronous rotation is small but not zero. (3) Considering the BYORP torque, our model shows agreement with the 1-degree of freedom adiabatic invariance theory in the outwards migration process, but an obvious difference in the inwards migration process. In particular, our studies show that the so-called 'long-term' stable equilibrium between the BYORP torque and the tidal torque is never a real equilibrium state, although the binary asteroid system can be captured in this state for quite a long time. (4) In case that the primary's angular velocity gradually reduces due to the YORP effect, the secondary's synchronous state may be broken when the primary's rotational motion crosses some major spin-orbit resonances.


title: Probing the superfluid-insulator phase transition by a non-Hermitian external field

abstract: We study the response of a thermal state of the Hubbard-like system to either global or local non-Hermitian perturbation, which coalesces the degenerate ground state within the $U(1)$ symmetry breaking phase. We show that the dynamical response of the system is strongly sensitive to the underlying quantum phase transition (QPT) from a Mott insulator to a superfluid state. The Uhlmann fidelity in the superfluid phase decays to a steady value determined by the order of the exceptional point (EP) within the subspace spanned by the degenerate ground states but remains almost unchanged in the Mott insulating phase. It demonstrates that the phase diagram at zero temperature is preserved even though a local probing field is applied. Specifically, two celebrated models including the Bose-Hubbard model and the Jaynes-Cummings-Hubbard model are employed to demonstrate this property in the finite-size system, wherein fluctuations of the boson and polariton number are observed based on EP dynamics. This work presents an alternative approach to probe the superfluid-insulator QPT at non-zero temperature.


title: NOEMA Observations of CO Emission in Arp 142 and Arp 238

abstract: Previous studies have shown significant differences in the enhancement of the star-formation rate (SFR) and the star-formation efficiency (SFE=SFR/M_mol) between spiral-spiral and spiral-elliptical mergers. In order to shed light on the physical mechanism of these differences, we present NOEMA observations of the molecular gas distribution and kinematics (linear resolutions of ~ 2kpc) in two representative close major-merger star-forming pairs: the spiral-elliptical pair Arp142 and the spiral-spiral pair Arp238. The CO in Arp142 is widely distributed over a highly distorted disk without any nuclear concentration, and an off-centric ring-like structure is discovered in channel maps. The SFE varies significantly within Arp142, with a starburst region (Region 1) near the eastern tip of the distorted disk showing an SFE ~0.3 dex above the mean of the control sample of isolated galaxies, and the SFE of the main disk (Region 4) 0.43 dex lower than the mean of the control sample. In contrast, the CO emission in Arp238 is detected only in two compact sources at the galactic centers. Compared to the control sample, Arp238-E shows an SFE enhancement of more than 1 dex whereas Arp238-W has an enhancement of ~0.7 dex. We suggest that the extended CO distribution and the large SFE variation in Arp142 are due to an expanding large-scale ring triggered by a recent high-speed head-on collision between the spiral galaxy and the elliptical galaxy, and the compact CO sources with high SFEs in Arp238 are associated with nuclear starbursts induced by gravitational tidal torques in a low-speed coplanar interaction.


title: LNL polycategories and doctrines of linear logic

abstract: We define and study LNL polycategories, which abstract the judgmental structure of classical linear logic with exponentials. Many existing structures can be represented as LNL polycategories, including LNL adjunctions, linear exponential comonads, LNL multicategories, IL-indexed categories, linearly distributive categories with storage, commutative and strong monads, CBPV-structures, models of polarized calculi, Freyd-categories, and skew multicategories, as well as ordinary cartesian, symmetric, and planar multicategories and monoidal categories, symmetric polycategories, and linearly distributive and *-autonomous categories. To study such classes of structures uniformly, we define a notion of LNL doctrine, such that each of these classes of structures can be identified with the algebras for some such doctrine. We show that free algebras for LNL doctrines can be presented by a sequent calculus, and that every morphism of doctrines induces an adjunction between their 2-categories of algebras.


title: Stability of isoperimetric inequalities for Laplace eigenvalues on surfaces

abstract: We prove stability estimates for the isoperimetric inequalities for the first and the second nonzero Laplace eigenvalues on surfaces, both globally and in a fixed conformal class. We employ the notion of eigenvalues of measures and show that if a normalized eigenvalue is close to its maximal value, the corresponding measure must be close in the Sobolev space $W^{-1,2}$ to the set of maximizing measures. In particular, this implies a qualitative stability result: metrics almost maximizing the normalized eigenvalue must be $W^{-1,2}$-close to a maximal metric. Following this approach, we prove sharp quantitative stability of the celebrated Hersch's inequality for the first eigenvalue on the sphere, as well as of its counterpart for the second eigenvalue. Similar results are also obtained for the precise isoperimetric eigenvalue inequalities on the projective plane, torus, and Klein bottle. The square of the $W^{-1,2}$ distance to a maximizing measure in these stability estimates is controlled by the difference between the normalized eigenvalue and its maximal value, indicating that the maxima are in a sense nondegenerate. We construct examples showing that the power of the distance can not be improved, and that the choice of the Sobolev space $W^{-1,2}$ is optimal.


title: A novel augmented Lagrangian method of multipliers for optimization with general inequality constraints

abstract: We introduce a twice differentiable augmented Lagrangian for nonlinear optimization with general inequality constraints and show that a strict local minimizer of the original problem is an approximate strict local solution of the augmented Lagrangian. A novel augmented Lagrangian method of multipliers (ALM) is then presented. Our method is originated from a generalization of the Hetenes-Powell augmented Lagrangian, and is a combination of the augmented Lagrangian and the interior-point technique. It shares a similar algorithmic framework with existing ALMs for optimization with inequality constraints, but it can use the second derivatives and does not depend on projections on the set of inequality constraints. In each iteration, our method solves a twice continuously differentiable unconstrained optimization subproblem on primal variables. The dual iterates, penalty and smoothing parameters are updated adaptively. The global and local convergence are analyzed. Without assuming any constraint qualification, it is proved that the proposed method has strong global convergence. The method may converge to either a Kurash-Kuhn-Tucker (KKT) point or a singular stationary point when the converging point is a minimizer. It may also converge to an infeasible stationary point of nonlinear program when the problem is infeasible. Furthermore, our method is capable of rapidly detecting the possible infeasibility of the solved problem. Under suitable conditions, it is locally linearly convergent to the KKT point, which is consistent with ALMs for optimization with equality constraints. The preliminary numerical experiments on some small benchmark test problems demonstrate our theoretical results.


title: EVPropNet: Detecting Drones By Finding Propellers For Mid-Air Landing And Following

abstract: The rapid rise of accessibility of unmanned aerial vehicles or drones pose a threat to general security and confidentiality. Most of the commercially available or custom-built drones are multi-rotors and are comprised of multiple propellers. Since these propellers rotate at a high-speed, they are generally the fastest moving parts of an image and cannot be directly "seen" by a classical camera without severe motion blur. We utilize a class of sensors that are particularly suitable for such scenarios called event cameras, which have a high temporal resolution, low-latency, and high dynamic range. In this paper, we model the geometry of a propeller and use it to generate simulated events which are used to train a deep neural network called EVPropNet to detect propellers from the data of an event camera. EVPropNet directly transfers to the real world without any fine-tuning or retraining. We present two applications of our network: (a) tracking and following an unmarked drone and (b) landing on a near-hover drone. We successfully evaluate and demonstrate the proposed approach in many real-world experiments with different propeller shapes and sizes. Our network can detect propellers at a rate of 85.1% even when 60% of the propeller is occluded and can run at upto 35Hz on a 2W power budget. To our knowledge, this is the first deep learning-based solution for detecting propellers (to detect drones). Finally, our applications also show an impressive success rate of 92% and 90% for the tracking and landing tasks respectively.


title: Self-organized structures of two-component laser fields and their active controls in a cold Rydberg atomic gas

abstract: We investigate the formation and control of stationary optical patterns in a cold Rydberg atomic gas via double electromagnetically induced transparency. We show that, through the modulational instability of plane-wave state of a laser field with two polarization components, the system undergoes a spontaneous symmetry breaking and hence the emergence of plentiful self-organized spatial optical structures, which can be manipulated by the ratio between the cross- and self-Kerr nonlinearities, the nonlocality degree of the Kerr nonlinearities, and the populations initially prepared in the two atomic ground states. Interestingly, a crossover from mixture to separation in space (optical phase separation) of the two polarization components occurs when the ratio between the cross- and self-Kerr nonlinearities exceeds a critical value. We also show that the system supports nonlocal two-component spatial optical solitons and vortices when the parameters of the system are selected suitably. The rich diversity and active controllability of the self-organized optical structures reported here provide a way for realizing novel optical patterns and solitons and their structural phase transitions based on Rydberg atomic gases.

categories: cond-mat.quant-gas nlin.AO nlin.PS physics.atom-ph physics.optics


title: Benchmarking Knowledge-driven Zero-shot Learning

abstract: External knowledge (a.k.a. side information) plays a critical role in zero-shot learning (ZSL) which aims to predict with unseen classes that have never appeared in training data. Several kinds of external knowledge, such as text and attribute, have been widely investigated, but they alone are limited with incomplete semantics. Some very recent studies thus propose to use Knowledge Graph (KG) due to its high expressivity and compatibility for representing kinds of knowledge. However, the ZSL community is still in short of standard benchmarks for studying and comparing different external knowledge settings and different KG-based ZSL methods. In this paper, we proposed six resources covering three tasks, i.e., zero-shot image classification (ZS-IMGC), zero-shot relation extraction (ZS-RE), and zero-shot KG completion (ZS-KGC). Each resource has a normal ZSL benchmark and a KG containing semantics ranging from text to attribute, from relational knowledge to logical expressions. We have clearly presented these resources including their construction, statistics, data formats and usage cases w.r.t. different ZSL methods. More importantly, we have conducted a comprehensive benchmarking study, with two general and state-of-the-art methods, two setting-specific methods and one interpretable method. We discussed and compared different ZSL paradigms w.r.t. different external knowledge settings, and found that our resources have great potential for developing more advanced ZSL methods and more solutions for applying KGs for augmenting machine learning. All the resources are available at https://github.com/China-UK-ZSL/Resources_for_KZSL.


title: Bottleneck Time Minimization for Distributed Iterative Processes: Speeding Up Gossip-Based Federated Learning on Networked Computers

abstract: We present a novel task scheduling scheme for accelerating computational applications involving distributed iterative processes that are executed on networked computing resources. Such an application consists of multiple tasks, each of which outputs data at each iteration to be processed by neighboring tasks; these dependencies between the tasks can be represented as a directed graph. We first mathematically formulate the problem as a Binary Quadratic Program (BQP), accounting for both computation and communication costs. We show that the problem is NP-hard. We then relax the problem as a Semi-Definite Program (SDP) and utilize a randomized rounding technique based on sampling from a suitably-formulated multi-variate Gaussian distribution. Furthermore, we derive the expected value of bottleneck time. Finally, we apply our proposed scheme on gossip-based federated learning as an application of iterative processes. Through numerical evaluations on the MNIST and CIFAR-10 datasets, we show that our proposed approach outperforms well-known scheduling techniques from distributed computing. In particular, for arbitrary settings, we show that it reduces bottleneck time by $91\%$ compared to HEFT and $84\%$ compared to throughput HEFT.


title: FallDeF5: A Fall Detection Framework Using 5G-based Deep Gated Recurrent Unit Networks

abstract: Fall prevalence is high among elderly people, which is challenging due to the severe consequences of falling. This is why rapid assistance is a critical task. Ambient assisted living (AAL) uses recent technologies such as 5G networks and the internet of medical things (IoMT) to address this research area. Edge computing can reduce the cost of cloud communication, including high latency and bandwidth use, by moving conventional healthcare services and applications closer to end-users. Artificial intelligence (AI) techniques such as deep learning (DL) have been used recently for automatic fall detection, as well as supporting healthcare services. However, DL requires a vast amount of data and substantial processing power to improve its performance for the IoMT linked to the traditional edge computing environment. This research proposes an effective fall detection framework based on DL algorithms and mobile edge computing (MEC) within 5G wireless networks, the aim being to empower IoMT-based healthcare applications. We also propose the use of a deep gated recurrent unit (DGRU) neural network to improve the accuracy of existing DL-based fall detection methods. DGRU has the advantage of dealing with time-series IoMT data, and it can reduce the number of parameters and avoid the vanishing gradient problem. The experimental results on two public datasets show that the DGRU model of the proposed framework achieves higher accuracy rates compared to the current related works on the same datasets.


title: Decentralized Blockchain-based model for Edge Computing

abstract: Blockchain technology is among the fastest-growing technologies in the world today. It has been adopted in diverse areas but mostly in financial systems, such as Bitcoin cryptocurrency. Therefore, it is a niche that has attracted interest from researchers from various fields, including computer science. Other areas where Blockchain is being embraced are the Smart Grid and Internet of Things (IoT) technologies, among others. While it is all good and improving many areas of applications, Blockchain still has some shortcomings. For example, it is not designed for high scalability when accommodating normal transactions. On the other hand, a parallel technology that has diverse applications in distributed networks better known as edge computing has emerged. Its main advantage is that it increases the speed of pf processes within those networks. However, like Blockchain, edge computing has its shortcomings. Its security systems and management systems have been found to be wanting. Hence the idea to integrate the two technologies and take advantage of their strengths. A blend of the two would lead to advanced network servers, huge data storage, and heightened security in transactions. However, this integration will best happen when some measures are taken. For example, there is a need to address scalability, resource management satisfactorily, and the security of the systems. To solve the integration problem, a decentralized Blockchain-based model of Edge computing is proposed in this paper.


title: Microbiome compositional analysis with logistic-tree normal models

abstract: Modern microbiome compositional data are often high-dimensional and exhibit complex dependency among the microbial taxa. However, existing statistical models for such data either do not adequately account for the dependency among the microbial taxa or lack computational scalability with respect to the number of taxa. This presents challenges in important applications such as association analysis between microbiome compositions and disease risk in which valid statistical analysis requires appropriately incorporating the "variance components" or "random effects" in the microbiome composition. We introduce a generative model, called the "logistic-tree normal" (LTN) model, that addresses this need. LTN marries two popular classes of models-namely, log-ratio normal (LN) and Dirichlet-tree (DT)-and inherits the key benefits of each. LTN incorporates the tree-based binomial decomposition as the DT does, but it jointly models the corresponding binomial probabilities using a (multivariate) logistic-normal distribution as in LN models. It therefore allows rich covariance structures as LN, along with computational efficiency realized through a Polya-Gamma augmentation on the binomial models associated with the tree splits. Accordingly, Bayesian inference on LTN can readily proceed by Gibbs sampling. LTN also allows common techniques for effective inference on high-dimensional data to be readily incorporated. We construct a general mixed-effects model using LTN to characterize compositional random effects, which allows flexible taxa covariance. We demonstrate its use in testing association between microbiome composition and disease risk as well as in estimating the covariance among taxa. We carry out an extensive case study using this LTN-enriched compositional mixed-effects model to analyze a longitudinal dataset from the T1D cohort of the DIABIMMUNE project.


title: Intrinsic Metallicity Variation in the Intermediate Mass Type II Globular Cluster NGC 1261

abstract: Globular Clusters (GCs) are now well known to almost universally show multiple popu-lations (MPs). The HST UV Legacy Survey of a large number of Galactic GCs in UV filters optimized to explore MPs finds that a small fraction of GCs, termed Type II, also display more complex, anomalous behavior. Several well-studied Type II GCs show intrinsic Fe abundance variations, suggesting that the other, less well-studied, Type II GCs should also exhibit similar behavior. Our aim is to perform the first detailed metallicity analysis of NGC 1261, an intermediate mass Type II GC, in order to determine if this object shows an intrinsic Fe variation. We determined the Fe abundance in eight red giant members using Magellan-MIKE and UVES-FLAMES high-resolution, high S/N spectroscopy. The full range of [Fe/H] for the entire sample from the spectra is from -1.05 to -1.43 dexwith an observed spread sigma_obs=0.133 dex. Compared with the total internal error of Sigma_tot=0.06,this indicates a significant intrinsic metallicity spread of Sigma_int=0.119 dex. We found a very similar variation in [Fe/H] using an independent method to derive the atmospheric parameters based on near-IR photometry. More importantly, the mean metallicity of the five presumed normal metallicity stars is -1.37+/-0.02, while that of the three presumed anomalous/highmetallicity stars is -1.18+/-0.09. This difference is significant at the $\pm$2.4Sigma level. We find indications from existing data of other Type II GCs that several of them presumedto have real metallicity spreads may in fact posses none. The minimum mass required for a GC to acquire an intrinsic Fe spread appears to be $\pm$10^5 Msun. We find no strong correlation betwee nmass and metallicity variation for Type II GCs. The metallicity spread is also independent of the fraction of anomalous stars within the Type II GCs and of GC origin.


title: Magnetic moments and transition magnetic moments of $P_c$ and $P_{cs}$ states

abstract: We study the magnetic moments and transition magnetic moments of $P_c$ and $P_{cs}$ states in the molecular picture. We first revisit the magnetic moments of $P_c$ states as the $S$ wave molecular states without coupled channel effects. The coupled channel effects and the $D$ wave contributions are then investigated carefully. The coupled channel effects contribute to the change of $0.1\sim 0.4$ nuclear magneton $\mu_N$ for most cases while the $D$ wave only induces the variation of less than $0.03 ~\mu_N$. In addition, we obtain the transition magnetic moments between different $P_c$ states and the related electromagnetic decay widths of $P_c'\to P_c\gamma$. The magnetic moments of $P_{cs}$ are much different for the assumption of spin being 1/2 and 3/2. The study of electromagnetic properties will help us disclose further the structure of these unconventional states.


title: Time-Domain Doppler Biomotion Detections Immune to Unavoidable DC Offsets

abstract: In the past decades, continuous Doppler radar sensor-based bio-signal detections have attracted many research interests. A typical example is the Doppler heartbeat detection. While significant progresses have been achieved, reliable, time-domain accurate demodulation of bio-signals in the presence of unavoidable DC offsets remains a technical challenge. Aiming to overcome this difficulty, we propose in this paper a novel demodulation algorithm that does not need to trace and eliminate dynamic DC offsets based on approximating segmented arcs in a quadrature constellation of sampling data to directional chords. Assisted by the principal component analysis, such chords and their directions can be deterministically determined. Simulations and experimental validations showed fully recovery of micron-level pendulum movements and strongly noised human heartbeats, verifying the effectiveness and accuracy of the proposed approach.