title: Reduction Theory of Algebraic Modules and their Successive Minima

abstract: Lattices defined as modules over algebraic rings or orders have garnered interest recently, particularly in the fields of cryptography and coding theory. Whilst there exist many attempts to generalise the conditions for LLL reduction to such lattices, there do not seem to be any attempts so far to generalise stronger notions of reduction such as Minkowski, HKZ and BKZ reduction. Moreover, most lattice reduction methods for modules over algebraic rings involve applying traditional techniques to the embedding of the module into real space, which distorts the structure of the algebra. In this paper, we generalise some classical notions of reduction theory to that of free modules defined over an order. Moreover, we extend the definitions of Minkowski, HKZ and BKZ reduction to that of such modules and show that bases reduced in this manner have vector lengths that can be bounded above by the successive minima of the lattice multiplied by a constant that depends on the algebra and the dimension of the module. In particular, we show that HKZ reduced bases are polynomially close to the successive minima of the lattice in terms of the module dimension. None of our definitions require the module to be embedded and thus preserve the structure of the module.


title: Crystal Structure of the $\tau_{11}$ Al$_4$Fe$_{1.7}$Si Phase from Neutron Diffraction and Ab Initio Calculations

abstract: The intermetallic $\tau_{11}$ Al$_4$Fe$_{1.7}$Si phase is of interest for high-temperature structural application due to its combination of low density and high strength. We determine the crystal structure of the $\tau_{11}$ phase through a combination of powder neutron diffraction and density functional theory calculations. Using Pawley and Rietveld refinements of the neutron diffraction data provides an initial crystal structure model. Since Al and Si have nearly identical neutron scattering lengths, we use density-functional calculations to determine their preferred site occupations. The $\tau_{11}$ phase exhibits a hexagonal crystal structure with space group $P6_3/mmc$ and lattice parameters of $a=7.478$~\AA\ and $c=7.472$~\AA. The structure comprises five Wyckoff positions; Al occupies the $6h$ and $12k$ sites, Fe the $2a$ and $6h$ sites, and Si the $2a$ sites. We observe site disorder and partial occupancies on all sites with a large fraction of 80\% Fe vacancies on the $2d$ sites, indicating an entropic stabilization of the $\tau_{11}$ phase at high temperature.


title: Committing to Interdependence: Implications from Game Theory for Human-Robot Trust

abstract: Human-robot interaction and game theory have developed distinct theories of trust for over three decades in relative isolation from one another. Human-robot interaction has focused on the underlying dimensions, layers, correlates, and antecedents of trust models, while game theory has concentrated on the psychology and strategies behind singular trust decisions. Both fields have grappled to understand over-trust and trust calibration, as well as how to measure trust expectations, risk, and vulnerability. This paper presents initial steps in closing the gap between these fields. Using insights and experimental findings from interdependence theory and social psychology, this work starts by analyzing a large game theory competition data set to demonstrate that the strongest predictors for a wide variety of human-human trust interactions are the interdependence-derived variables for commitment and trust that we have developed. It then presents a second study with human subject results for more realistic trust scenarios, involving both human-human and human-machine trust. In both the competition data and our experimental data, we demonstrate that the interdependence metrics better capture social `overtrust' than either rational or normative psychological reasoning, as proposed by game theory. This work further explores how interdependence theory--with its focus on commitment, coercion, and cooperation--addresses many of the proposed underlying constructs and antecedents within human-robot trust, shedding new light on key similarities and differences that arise when robots replace humans in trust interactions.


title: MeerKAT view of the diffuse radio sources in Abell 3667 and their interactions with the thermal plasma

abstract: During their lifetime, galaxy clusters grow through the accretion of matter from the filaments of the large scale structure and from mergers with other clusters. These mergers release a large amount of energy into the intracluster medium (ICM) through merger shocks and turbulence. These phenomena are associated with the formation of radio sources known as radio relics and radio halos, respectively. Radio relics and halos are unique proxies to study the complex properties of these dynamically active regions of clusters and in general the micro physics of the ICM. Abell 3667 is a spectacular examples of a merging system hosting a large pair of radio relics. Due to its proximity (z=0.0553) and large mass, the system enables the study of these sources to a uniquely high level of detail. We observed Abell 3667 with MeerKAT as part of the MeerKAT Galaxy Cluster Legacy Survey. We used these data to study the large scale emission of the cluster, including its polarisation and spectral properties. We present the most detailed view of the radio relic system in Abell 3667 to date, with a resolution reaching 3 kpc. The relics are filled with a network of filaments with different spectral and polarisation properties that are likely associated with multiple regions of particle acceleration and local enhancements of the magnetic field. Conversely, the magnetic field in the space between filaments has strengths close to that expected in unperturbed regions at the same cluster-centric distance. Comparisons with MHD simulations supports the idea of filaments as multiple acceleration sites. Our observations also confirm the presence of an elongated radio halo, developed in the wake of the bullet-like sub-cluster that merged from the South-East. Finally, we associate the process of magnetic draping to a thin polarised radio source surrounding the remnant of the bullet's cool core.


title: Absolute and Relative Bias in Eight Common Observational Study Designs: Evidence from a Meta-analysis

abstract: Observational studies are needed when experiments are not possible. Within study comparisons (WSC) compare observational and experimental estimates that test the same hypothesis using the same treatment group, outcome, and estimand. Meta-analyzing 39 of them, we compare mean bias and its variance for the eight observational designs that result from combining whether there is a pretest measure of the outcome or not, whether the comparison group is local to the treatment group or not, and whether there is a relatively rich set of other covariates or not. Of these eight designs, one combines all three design elements, another has none, and the remainder include any one or two. We found that both the mean and variance of bias decline as design elements are added, with the lowest mean and smallest variance in a design with all three elements. The probability of bias falling within 0.10 standard deviations of the experimental estimate varied from 59 to 83 percent in Bayesian analyses and from 86 to 100 percent in non-Bayesian ones -- the ranges depending on the level of data aggregation. But confounding remains possible due to each of the eight observational study design cells including a different set of WSC studies.


title: Predictive coding, precision and natural gradients

abstract: There is an increasing convergence between biologically plausible computational models of inference and learning with local update rules and the global gradient-based optimization of neural network models employed in machine learning. One particularly exciting connection is the correspondence between the locally informed optimization in predictive coding networks and the error backpropagation algorithm that is used to train state-of-the-art deep artificial neural networks. Here we focus on the related, but still largely under-explored connection between precision weighting in predictive coding networks and the Natural Gradient Descent algorithm for deep neural networks. Precision-weighted predictive coding is an interesting candidate for scaling up uncertainty-aware optimization -- particularly for models with large parameter spaces -- due to its distributed nature of the optimization process and the underlying local approximation of the Fisher information metric, the adaptive learning rate that is central to Natural Gradient Descent. Here, we show that hierarchical predictive coding networks with learnable precision indeed are able to solve various supervised and unsupervised learning tasks with performance comparable to global backpropagation with natural gradients and outperform their classical gradient descent counterpart on tasks where high amounts of noise are embedded in data or label inputs. When applied to unsupervised auto-encoding of image inputs, the deterministic network produces hierarchically organized and disentangled embeddings, hinting at the close connections between predictive coding and hierarchical variational inference.


title: Associative algebras and intertwining operators

abstract: Let $V$ be a vertex operator algebra and $A^{\infty}(V)$ and $A^{N}(V)$ for $N\in \mathbb{N}$ the associative algebras introduced by the author in [H5]. For a lower-bounded generalized $V$-module $W$, we give $W$ a structure of graded $A^{\infty}(V)$-module and we introduce an $A^{\infty}(V)$-bimodule $A^{\infty}(W)$ and an $A^{N}(V)$-bimodule $A^{N}(W)$. We prove that the space of (logarithmic) intertwining operators of type $\binom{W_{3}}{W_{1}W_{2}}$ for lower-bounded generalized $V$-modules $W_{1}$, $W_{2}$ and $W_{3}$ is isomorphic to the space $\hom_{A^{\infty}(V)}(A^{\infty}(W_{1})\otimes_{A^{\infty}(V)}W_{2}, W_{3})$. Assuming that $W_{2}$ and $W_{3}'$ are equivalent to certain universal lower-bounded generalized $V$-modules generated by their $A^{N}(V)$-submodules consisting of elements of levels less than or equal to $N\in \mathbb{N}$, we also prove that the space of (logarithmic) intertwining operators of type $\binom{W_{3}}{W_{1}W_{2}}$ is isomorphic to the space of $\hom_{A^{N}(V)}(A^{N}(W_{1})\otimes_{A^{N}(V)}\Omega_{N}^{0}(W_{2}), \Omega_{N}^{0}(W_{3}))$.


title: New results and open problems on subgraph centrality

abstract: Subgraph centrality, introduced by Estrada and Rodr\'iguez-Vel\'azquez in [12], has become a widely used centrality measure in the analysis of networks, with applications in biology, neuroscience, economics and many other fields. It is also worthy of study from a strictly mathematical point of view, in view of its connections to topics in spectral graph theory, number theory, analytic matrix functions, and combinatorics. In this paper we present some new results and a list of open questions about subgraph centrality and other node centrality measures based on graph walks.


title: Learning Interpretation with Explainable Knowledge Distillation

abstract: Knowledge Distillation (KD) has been considered as a key solution in model compression and acceleration in recent years. In KD, a small student model is generally trained from a large teacher model by minimizing the divergence between the probabilistic outputs of the two. However, as demonstrated in our experiments, existing KD methods might not transfer critical explainable knowledge of the teacher to the student, i.e. the explanations of predictions made by the two models are not consistent. In this paper, we propose a novel explainable knowledge distillation model, called XDistillation, through which both the performance the explanations' information are transferred from the teacher model to the student model. The XDistillation model leverages the idea of convolutional autoencoders to approximate the teacher explanations. Our experiments shows that models trained by XDistillation outperform those trained by conventional KD methods not only in term of predictive accuracy but also faithfulness to the teacher models.


title: Leptonic non-thermal emission from supernova remnants evolving in the circumstellar magnetic field

abstract: The very-high-energy (VHE; E > 100 GeV) gamma-ray emission observed from a number of Supernova remnants (SNRs) indicates particle acceleration to high energies at the shock of the remnants and a potentially significant contribution to Galactic cosmic rays. It is extremely difficult to determine whether protons (through hadronic interactions and subsequent pion decay) or electrons (through inverse Compton scattering on ambient photon fields) are responsible for this emission. For a successful diagnostic, a good understanding of the spatial and energy distribution of the underlying particle population is crucial. Most SNRs are created in core-collapse explosions and expand into the wind bubble of their progenitor stars. This circumstellar medium features a complex spatial distribution of gas and magnetic field which naturally strongly affects the resulting particle population. In this work, we conduct a detailed study of the spectro-spatial evolution of the electrons accelerated at the forward shock of core-collapse SNRs and their non-thermal radiation, using the RATPaC code that is designed for the time- and spatially dependent treatment of particle acceleration at SNR shocks. We focus on the impact of the spatially inhomogeneous magnetic field through the efficiency of diffusion and synchrotron cooling. It is demonstrated that the structure of the circumstellar magnetic field can leave strong signatures in the spectrum and morphology of the resulting non-thermal emission.


title: KdV Equation Model in Open Cylindrical Channel under Precession

abstract: A new model for Korteweg and de-Vries equation (KdV) is derived. The system under study is an open channel consisting of two concentric cylinders, rotating about their vertical axis, which is tilted by slope {\tau} from the inertial vertical z, in uniform rate {\Omega}_1={\tau}{\Omega}, and the whole tank is elevated over other table rotating at rate {\Omega}. Under these conditions, a set of Kelvin waves is formed on the free surface depending on the angle of tilt, characterized by the slope {\tau}, volume of water, and rotation rate. The resonant mode in the system appears in the form of a single Kelvin solitary wave, whose amplitude satisfies the Korteweg-de Vries equation with forced term. The equation was derived following classical perturbation methods, the additional term made the equation a non-integrable one, that cannot be solved without the help of numerical methods. Invoking the simple finite difference scheme method, it was found that the numerical results are in a good agreement with the experiment.


title: Light-cone PDFs and GPDs from Lattice QCD

abstract: In this article, we review recent lattice calculations on the $x$-dependence of PDFs and GPDs from lattice QCD.


title: Fractional resonances and prethermal states in Floquet systems

abstract: In periodically-driven quantum systems, resonances can induce exotic nonequilibrium behavior and new phases of matter without static analog. We report on the emergence of fractional and integer resonances in a broad class of many-body Hamiltonians with a modulated hopping with a frequency that is either a fraction or an integer of the on-site interaction. We contend that there is a fundamental difference between these resonances when interactions bring the system to a Floquet prethermal state. Second-order processes dominate the dynamics in the fractional resonance case, leading to less entanglement and more localized quantum states than in the integer resonance case dominated by first-order processes. We demonstrate the dominating emergence of fractional resonances using the Magnus expansion of the effective Hamiltonian and quantify their effects on the many-body dynamics via quantum states' von Neumann entropy and Loschmidt echo. Our findings reveal novel features of the nonequilibrium quantum many-body system, such as the coexistence of Floquet prethermalization and localization, that may allow to development of quantum memories for quantum technologies and quantum information processing.


title: A note on the multiplicities of the determinantal thickenings of maximal minors

abstract: Let $S=\mathbb{C}[x_{ij}]$ be a polynomial ring of $m\times n$ variables over $\mathbb{C}$ and let $I$ be the determinantal ideal of maximal minors of $S$. Using the representation theoretic techniques introduced in arXiv:1305.1719, arXiv:1309.0617 and arXiv:1611.00415, we prove the existence of the generalized $j$-multiplicities $\epsilon^j(I)$ defined by Dao and Monta\~{n}o in arXiv:1705.05033. We will also give a closed formula of $\epsilon^j(I)$, which generalized the results in arXiv:1308.0582 and arXiv:1912.02917 in the maximal minors case.


title: Diffraction integral computation using sinc approximation

abstract: We propose a method based on sinc series approximations for computing the Rayleigh-Sommerfeld and Fresnel diffraction integrals of optics. The diffraction integrals are given in terms of a convolution, and our proposed numerical approach is not only super-algebraically convergent, but it also satisfies an important property of the convolution -- namely, the preservation of bandwidth. Furthermore, the accuracy of the proposed method depends only on how well the source field is approximated; it is independent of wavelength, propagation distance, and observation plane discretization. In contrast, methods based on the fast Fourier transform (FFT), such as the angular spectrum method (ASM) and its variants, approximate the optical fields in the source and observation planes using Fourier series. We will show that the ASM introduces artificial periodic boundary conditions and violates the preservation of bandwidth property, resulting in limited accuracy which decreases for longer propagation distances. The sinc-based approach avoids both of these problems. Numerical results are presented for Gaussian beam propagation and circular aperture diffraction to demonstrate the high-order accuracy of the sinc method for both short-range and long-range propagation. For comparison, we also present numerical results obtained with the angular spectrum method.


title: Micro-X Sounding Rocket Payload Re-flight Progress

abstract: Micro-X is an X-ray sounding rocket payload that had its first flight on July 22, 2018. The goals of the first flight were to operate a transition edge sensor (TES) X-ray microcalorimeter array in space and take a high-resolution spectrum of the Cassiopeia A supernova remnant. The first flight was considered a partial success. The array and its time-division multiplexing readout system were successfully operated in space, but due to a failure in the attitude control system, no time on-target was acquired. A re-flight has been scheduled for summer 2022. Since the first flight, modifications have been made to the detector systems to improve noise and reduce the susceptibility to magnetic fields. The three-stage SQUID circuit, NIST MUX06a, has been replaced by a two-stage SQUID circuit, NIST MUX18b. The initial laboratory results for the new detector system will be presented in this paper.


title: Morphological computation and decentralized learning in a swarm of sterically interacting robots

abstract: Whereas naturally occurring swarms thrive when crowded, physical interactions in robotic swarms are either avoided or carefully controlled, thus limiting their operational density. Here we present a mechanical design rule that allows robots to act in a collision-dominated environment. We introduce the Morphobots -- a robotic swarm platform developed to implement embodied computation through a morpho-functional design. By engineering a 3D-printed exoskeleton we encode a re-orientation response to an external body force (such as gravity) or a surface force (such as a collision). We show that the force-orientation response is generic, and can augment existing swarm-robotic platforms (e.g Kilobots) as well as custom robots even 10 times larger. At the individual level, the exoskeleton improves the motility and stability, and also allows to encode two contrasting dynamical behaviors in response to an external force or a collision (including collision with a wall or a movable obstacle, and on a dynamically tilting plane). This force-orientation response adds a mechanical layer to the robot's sense-act cycle at the swarm level, leveraging steric interactions for collective phototaxis when crowded. Enabling collisions also promotes information flow, facilitating online distributed learning. Each robot runs an embedded algorithm that ultimately optimizes collective performance. We identify an effective parameter that controls the force-orientation response and explore its implications in swarms that transition from dilute to crowded. \green{Experimenting with both physical swarms (of up to 64 robots), and simulated swarms (of up to 8192 agents) show that the effect of morphological-computation increases with growing swarm size.


title: Homotopy invariance of tame homotopy groups of regular schemes

abstract: The \'etale homotopy groups of schemes as defined by Artin and Mazur have the disadvantage of being homotopy invariant only in characteristic zero. This and other related problems led to the definition of the tame topology which is coarser than the \'etale topology by disallowing wild ramification along the boundary of compactifications. The object of this paper is to show that the associated tame homotopy groups are indeed ($\mathbb{A}^1$-)homotopy invariant, at least for regular schemes.


title: Nanoscale Optical Imaging of 2D Semiconductor Stacking Orders by Exciton-Enhanced Second Harmonic Generation

abstract: Second harmonic generation (SHG) is a nonlinear optical response arising exclusively from broken inversion symmetry in the electric-dipole limit. Recently, SHG has attracted widespread interest as a versatile and noninvasive tool for characterization of crystal symmetry and emerging ferroic or topological orders in quantum materials. However, conventional far-field optics is unable to probe local symmetry at the deep subwavelength scale. Here, we demonstrate near-field SHG imaging of 2D semiconductors and heterostructures with the spatial resolution down to 20 nm using a scattering-type nano-optical apparatus. We show that near-field SHG efficiency is greatly enhanced by excitons in atomically thin transition metal dichalcogenides. Furthermore, by correlating nonlinear and linear scattering-type nano-imaging, we resolve nanoscale variations of interlayer stacking order in bilayer WSe2, and reveal the stacking-tuned excitonic light-matter-interactions. Our work demonstrates nonlinear optical interrogation of crystal symmetry and structure-property relationships at the nanometer length scales relevant to emerging properties in quantum materials.

categories: physics.optics cond-mat.mes-hall cond-mat.mtrl-sci physics.app-ph


title: Human irrationality: both bad and good for reward inference

abstract: Assuming humans are (approximately) rational enables robots to infer reward functions by observing human behavior. But people exhibit a wide array of irrationalities, and our goal with this work is to better understand the effect they can have on reward inference. The challenge with studying this effect is that there are many types of irrationality, with varying degrees of mathematical formalization. We thus operationalize irrationality in the language of MDPs, by altering the Bellman optimality equation, and use this framework to study how these alterations would affect inference. We find that wrongly modeling a systematically irrational human as noisy-rational performs a lot worse than correctly capturing these biases -- so much so that it can be better to skip inference altogether and stick to the prior! More importantly, we show that an irrational human, when correctly modelled, can communicate more information about the reward than a perfectly rational human can. That is, if a robot has the correct model of a human's irrationality, it can make an even stronger inference than it ever could if the human were rational. Irrationality fundamentally helps rather than hinder reward inference, but it needs to be correctly accounted for.


title: Fabrication of voltage gated spin Hall nano-oscillators

abstract: We demonstrate an optimized fabrication process for electric field (voltage gate) controlled nano-constriction spin Hall nano-oscillators (SHNOs), achieving feature sizes of <30 nm with easy to handle ma-N 2401 e-beam lithography negative tone resist. For the nanoscopic voltage gates, we utilize a two-step tilted ion beam etching approach and through-hole encapsulation using 30 nm HfOx. The optimized tilted etching process reduces sidewalls by 75% compared to no tilting. Moreover, the HfOx encapsulation avoids any sidewall shunting and improves gate breakdown. Our experimental results on W/CoFeB/MgO/SiO2 SHNOs show significant frequency tunability (6 MHz/V) even for moderate perpendicular magnetic anisotropy. Circular patterns with diameter of 45 nm are achieved with an aspect ratio better than 0.85 for 80% of the population. The optimized fabrication process allows incorporating a large number of individual gates to interface to SHNO arrays for unconventional computing and densely packed spintronic neural networks.


abstract: This paper examines the interdisciplinary research question of how to integrate Computational Argumentation, as studied in AI, with Cognition, as can be found in Cognitive Science, Linguistics, and Philosophy. It stems from the work of the 1st Workshop on Computational Argumentation and Cognition (COGNITAR), which was organized as part of the 24th European Conference on Artificial Intelligence (ECAI), and took place virtually on September 8th, 2020. The paper begins with a brief presentation of the scientific motivation for the integration of Computational Argumentation and Cognition, arguing that within the context of Human-Centric AI the use of theory and methods from Computational Argumentation for the study of Cognition can be a promising avenue to pursue. A short summary of each of the workshop presentations is given showing the wide spectrum of problems where the synthesis of the theory and methods of Computational Argumentation with other approaches that study Cognition can be applied. The paper presents the main problems and challenges in the area that would need to be addressed, both at the scientific level but also at the epistemological level, particularly in relation to the synthesis of ideas and approaches from the various disciplines involved.


title: Through-Foliage Tracking with Airborne Optical Sectioning

abstract: Detecting and tracking moving targets through foliage is difficult, and for many cases even impossible in regular aerial images and videos. We present an initial light-weight and drone-operated 1D camera array that supports parallel synthetic aperture aerial imaging. Our main finding is that color anomaly detection benefits significantly from image integration when compared to conventional raw images or video frames (on average 97% vs. 42% in precision in our field experiments). We demonstrate, that these two contributions can lead to the detection and tracking of moving people through densely occluding forest.


title: A new proof of reversibility of $SLE_{\kappa}$ for $\kappa \leq 4$

abstract: We give a new proof of the reversibility of the Schramm Loewner evolution for $\kappa \leq 4$. The main ideas used in the proof are similar to those used in the original proof of this result, given by Zhan.


title: Adversarially Robust Learning for Security-Constrained Optimal Power Flow

abstract: In recent years, the ML community has seen surges of interest in both adversarially robust learning and implicit layers, but connections between these two areas have seldom been explored. In this work, we combine innovations from these areas to tackle the problem of N-k security-constrained optimal power flow (SCOPF). N-k SCOPF is a core problem for the operation of electrical grids, and aims to schedule power generation in a manner that is robust to potentially k simultaneous equipment outages. Inspired by methods in adversarially robust training, we frame N-k SCOPF as a minimax optimization problem - viewing power generation settings as adjustable parameters and equipment outages as (adversarial) attacks - and solve this problem via gradient-based techniques. The loss function of this minimax problem involves resolving implicit equations representing grid physics and operational decisions, which we differentiate through via the implicit function theorem. We demonstrate the efficacy of our framework in solving N-3 SCOPF, which has traditionally been considered as prohibitively expensive to solve given that the problem size depends combinatorially on the number of potential outages.


title: Accounting for data heterogeneity in integrative analysis and prediction methods: An application to Chronic Obstructive Pulmonary Disease

abstract: Epidemiologic and genetic studies in chronic obstructive pulmonary disease (COPD) and many complex diseases suggest subgroup disparities (e.g., by sex). We consider this problem from the standpoint of integrative analysis where we combine information from different views (e.g., genomics, proteomics, clinical data). Existing integrative analysis methods ignore the heterogeneity in subgroups, and stacking the views and accounting for subgroup heterogeneity does not model the association among the views. To address analytical challenges in the problem of our interest, we propose a statistical approach for joint association and prediction that leverages the strengths in each view to identify molecular signatures that are shared by and specific to males and females and that contribute to the variation in COPD, measured by airway wall thickness. HIP (Heterogeneity in Integration and Prediction) accounts for subgroup heterogeneity, allows for sparsity in variable selection, is applicable to multi-class and to univariate or multivariate continuous outcomes, and incorporates covariate adjustment. We develop efficient algorithms in PyTorch. Our COPD findings have identified several proteins, genes, and pathways that are common and specific to males and females, some of which have been implicated in COPD, while others could lead to new insights into sex differences in COPD mechanisms.


title: Bertrand's Postulate for Carmichael Numbers

abstract: Alford, Granville, and Pomerance proved that there are infinitely many Carmichael numbers. In the same paper, they ask if a statement analogous to Bertrand's postulate could be proven for Carmichael numbers. In this paper, we answer this question, proving the stronger statement that for all $\delta>0$ and $x$ sufficiently large in terms of $\delta$, there exist at least $e^{\frac{\log x}{(\log\log x)^{2+\delta}}}$ Carmichael numbers between $x$ and $x+\frac{x}{(\log x)^{\frac{1}{2+\delta}}}$.


title: Synchronization of networks of piecewise-smooth systems

abstract: We study convergence in networks of piecewise-smooth (PWS) systems that commonly arise in applications to model dynamical systems whose evolution is affected by macroscopic events such as switches and impacts. Existing approaches were typically oriented toward guaranteeing global bounded synchronizability, local stability of the synchronization manifold, or achieving synchronization by exerting a control action on each node. Here we start by generalizing existing results on QUAD systems to the case of PWS systems, accounting for a large variety of nonlinear coupling laws. Then, we propose that a discontinuous coupling can be used to guarantee global synchronizability of a network of N PWS agents under mild assumptions on the individual dynamics. We provide extensive numerical simulations to gain insights on larger networks.


title: On Krull-Schmidt bicategories

abstract: We study the existence and uniqueness of direct sum decompositions in additive bicategories. We find a simple definition of Krull-Schmidt bicategories, for which we prove the uniqueness of decompositions into indecomposable objects as well as a characterization in terms of splitting of idempotents and properties of 2-cell endomorphism rings. Examples of Krull-Schmidt bicategories abound, with many arising from the various flavors of 2-dimensional linear representation theory.


title: Axion-fermion Coupling and Dyon Charge as Physical Signatures of a Space-time Inner Symmetry

abstract: In this paper we intend to complement the identification given in Kuerten & Fernandes-Silva [1] (Mod. Phys. Lett. A. Vol. 33, No. 16) which relates the axion to a metric spinor phase by means of Maxwell's theory in the Infeld-van der Waerden's {\gamma}-formalism. Thus, we obtain two alternative identifications: the first focuses on Dirac's theory so that when obtaining an axion-like phase-fermion coupling, we achieve the first identification, and the last one investigates the phase behavior under Peccei-Quinn rotations in order to show that the phase changes as an axion pseudoparticle. With the formal aspects established, we also study the semiclassical fermion-photon system to demonstrate that the magnetic monopole current defined in [1] has dyon charge in flat universe and acquires a Witten effect form when there is a demand for chiral symmetry.


title: On the complexity of SAT

abstract: In this paper, we prove that no deterministic algorithm can solve SAT in polynomial time in the number of boolean variables.


title: Hierarchical clustering by aggregating representatives in sub-minimum-spanning-trees

abstract: One of the main challenges for hierarchical clustering is how to appropriately identify the representative points in the lower level of the cluster tree, which are going to be utilized as the roots in the higher level of the cluster tree for further aggregation. However, conventional hierarchical clustering approaches have adopted some simple tricks to select the "representative" points which might not be as representative as enough. Thus, the constructed cluster tree is less attractive in terms of its poor robustness and weak reliability. Aiming at this issue, we propose a novel hierarchical clustering algorithm, in which, while building the clustering dendrogram, we can effectively detect the representative point based on scoring the reciprocal nearest data points in each sub-minimum-spanning-tree. Extensive experiments on UCI datasets show that the proposed algorithm is more accurate than other benchmarks. Meanwhile, under our analysis, the proposed algorithm has O(nlogn) time-complexity and O(logn) space-complexity, indicating that it has the scalability in handling massive data with less time and storage consumptions.


title: Computing with Cognitive States

abstract: Basic experimental findings about human working memory can be described by an algebra built on high-dimensional binary states, representing information items, and two operations: multiplication for binding and addition for bundling. In contrast to common VSA algebras, bundling is not associative. Consequently bundling a sequence of items preserves their sequential ordering. The cognitive states representing a memorised list exhibit a primacy as well as a recency gradient. The typical concave-up and asymmetrically shaped serial position curve is derived as a linear combination of those gradients. Quantitative implications of the algebra are shown to agree well with empirical data from basic cognitive tasks including storage and retrieval of information in human working memory.


title: Real topological Hochschild homology via the norm and Real Witt vectors

abstract: We prove that Real topological Hochschild homology can be characterized as the norm from the cyclic group of order $2$ to the orthogonal group $O(2)$. From this perspective, we then prove a multiplicative double coset formula for the restriction of this norm to dihedral groups of order $2m$. This informs our new definition of Real Hochschild homology of rings with anti-involution, which we show is the algebraic analogue of Real topological Hochschild homology. Using extra structure on Real Hochschild homology, we define a new theory of $p$-typical Witt vectors of rings with anti-involution. We end with an explicit computation of the degree zero $D_{2m}$-Mackey functor homotopy groups of $\operatorname{THR}(\underline{\mathbb{Z}})$ for $m$ odd. This uses a Tambara reciprocity formula for sums for general finite groups, which may be of independent interest.


title: Exploiting all samples in low-resource sentence classification: early stopping and initialization parameters

abstract: In low resource settings, deep neural models have often shown lower performance due to overfitting. The primary method to solve the overfitting problem is to generalize model parameters. To this end, many researchers have depended on large external resources with various manipulation techniques. In this study, we discuss how to exploit all available samples in low resource settings, without external datasets and model manipulation. This study focuses on natural language processing task. We propose a simple algorithm to find out good initialization parameters that improve robustness to a small sample set. We apply early stopping techniques that enable the use of all samples for training. Finally, the proposed learning strategy is to train all samples with the good initialization parameters and stop the model with the early stopping techniques. Extensive experiments are conducted on seven public sentence classification datasets, and the results demonstrate that the proposed learning strategy achieves better performance than several state-of-the-art works across the seven datasets.


title: Adaptive and quasi-sliding control of shimmy in landing gears

abstract: Shimmy is a dangerous phenomenon that occurs when aircraft's nose landing gears oscillate in a rapid and uncontrollable fashion. In this paper, we propose the use of two nonlinear control approaches (zero average control and model reference adaptive control based on minimal control synthesis) as simple yet effective strategies to suppress undesired oscillations, even in the presence of uncertainties and partial state measurements. Numerical results are presented to validate the proposed control approaches.


title: Universal Families of Eulerian Multiple Zeta Values in Positive Characteristics

abstract: We study positive characteristic multiple zeta values associated to general curves over $\mathbb F_q$ together with an $\mathbb F_q$-rational point $\infty$ as introduced by Thakur. For the case of the projective line these values were defined as analogues of classical multiple zeta values. In the present paper we first establish a general non-commutative factorization of exponential series associated to certain lattices of rank one. Next we introduce universal families of multiple zeta values of Thakur and show that they are Eulerian in full generality. In particular, we prove a conjecture of Lara Rodriguez and Thakur arXiv:2003.12910. One of the main ingredients of the proofs is the notion of L-series in Tate algebras introduced by the third author arXiv:1107.4511 in 2012.


title: Control Barrier Function Augmentation in Sampling-based Control Algorithm for Sample Efficiency

abstract: For a nonlinear stochastic path planning problem, sampling-based algorithms generate thousands of random sample trajectories to find the optimal path while guaranteeing safety by Lagrangian penalty methods. However, the sampling-based algorithm can perform poorly in obstacle-rich environments because most samples might violate safety constraints, invalidating the corresponding samples. To improve the sample efficiency of sampling-based algorithms in cluttered environments, we propose an algorithm based on model predictive path integral control and control barrier functions. The proposed algorithm needs fewer samples and time-steps and has a better performance in cluttered environments compared to the original model predictive path integral control algorithm.


title: Meshfree implementation of the cardiac monodomain model through the Fragile Points Method

abstract: Meshfree methods for in silico modelling and simulation of cardiac electrophysiology are gaining more and more popularity. These methods do not require a mesh and are more suitable than the Finite Element Method (FEM) to simulate the activity of complex geometrical structures like the human heart. However, challenges such as numerical integration accuracy and time efficiency remain and limit their applicability. Recently, the Fragile Points Method (FPM) has been introduced in the meshfree methods family. It uses local, simple, polynomial, discontinuous functions to construct trial and test functions in the Galerkin weak form. This allows for accurate integration and improved efficiency while enabling the imposition of essential and natural boundary conditions as in the FEM. In this work, we consider the application of FPM for cardiac electrophysiology simulation. We derive the cardiac monodomain model using the FPM formulation and we solve several benchmark problems in 2D and 3D. We show that FPM leads to solutions of similar accuracy and efficiency with FEM while alleviating the need for a mesh. Additionally, FPM demonstrates better convergence than FEM in the considered benchmarks.


title: Reflectance of Rhenium as a Function of Pressure in a Diamond Anvil Cell

abstract: We have measured the reflectance of rhenium in the visible region to pressures up to 100 GPa in a diamond anvil cell (DAC). By photographing the reflecting surface, we visually show that there are challenges to obtaining accurate values in a DAC for several reasons and propose some useful procedures. We also show that knowledge of the reflectance of rhenium can overcome the problem of absorption of light by diamonds when studying the reflectance of materials at high pressure in a DAC.


title: Scalable Diverse Model Selection for Accessible Transfer Learning

abstract: With the preponderance of pretrained deep learning models available off-the-shelf from model banks today, finding the best weights to fine-tune to your use-case can be a daunting task. Several methods have recently been proposed to find good models for transfer learning, but they either don't scale well to large model banks or don't perform well on the diversity of off-the-shelf models. Ideally the question we want to answer is, "given some data and a source model, can you quickly predict the model's accuracy after fine-tuning?" In this paper, we formalize this setting as "Scalable Diverse Model Selection" and propose several benchmarks for evaluating on this task. We find that existing model selection and transferability estimation methods perform poorly here and analyze why this is the case. We then introduce simple techniques to improve the performance and speed of these algorithms. Finally, we iterate on existing methods to create PARC, which outperforms all other methods on diverse model selection. We have released the benchmarks and method code in hope to inspire future work in model selection for accessible transfer learning.


title: RLOps: Development Life-cycle of Reinforcement Learning Aided Open RAN

abstract: Radio access network (RAN) technologies continue to evolve, with Open RAN gaining the most recent momentum. In the O-RAN specifications, the RAN intelligent controllers (RICs) are software-defined orchestration and automation functions for the intelligent management of RAN. This article introduces principles for machine learning (ML), in particular, reinforcement learning (RL) applications in the O-RAN stack. Furthermore, we review the state-of-the-art research in wireless networks and cast it onto the RAN framework and the hierarchy of the O-RAN architecture. We provide a taxonomy for the challenges faced by ML/RL models throughout the development life-cycle: from the system specification to production deployment (data acquisition, model design, testing and management, etc.). To address the challenges, we integrate a set of existing MLOps principles with unique characteristics when RL agents are considered. This paper discusses a systematic model development, testing and validation life-cycle, termed: RLOps. We discuss fundamental parts of RLOps, which include: model specification, development, production environment serving, operations monitoring and safety/security. Based on these principles, we propose the best practices for RLOps to achieve an automated and reproducible model development process. At last, a holistic data analytics platform rooted in the O-RAN deployment is designed and implemented, aiming to embrace and fulfil the aforementioned principles and best practices of RLOps.


title: Neural Population Geometry Reveals the Role of Stochasticity in Robust Perception

abstract: Adversarial examples are often cited by neuroscientists and machine learning researchers as an example of how computational models diverge from biological sensory systems. Recent work has proposed adding biologically-inspired components to visual neural networks as a way to improve their adversarial robustness. One surprisingly effective component for reducing adversarial vulnerability is response stochasticity, like that exhibited by biological neurons. Here, using recently developed geometrical techniques from computational neuroscience, we investigate how adversarial perturbations influence the internal representations of standard, adversarially trained, and biologically-inspired stochastic networks. We find distinct geometric signatures for each type of network, revealing different mechanisms for achieving robust representations. Next, we generalize these results to the auditory domain, showing that neural stochasticity also makes auditory models more robust to adversarial perturbations. Geometric analysis of the stochastic networks reveals overlap between representations of clean and adversarially perturbed stimuli, and quantitatively demonstrates that competing geometric effects of stochasticity mediate a tradeoff between adversarial and clean performance. Our results shed light on the strategies of robust perception utilized by adversarially trained and stochastic networks, and help explain how stochasticity may be beneficial to machine and biological computation.


title: GraSSNet: Graph Soft Sensing Neural Networks

abstract: In the era of big data, data-driven based classification has become an essential method in smart manufacturing to guide production and optimize inspection. The industrial data obtained in practice is usually time-series data collected by soft sensors, which are highly nonlinear, nonstationary, imbalanced, and noisy. Most existing soft-sensing machine learning models focus on capturing either intra-series temporal dependencies or pre-defined inter-series correlations, while ignoring the correlation between labels as each instance is associated with multiple labels simultaneously. In this paper, we propose a novel graph based soft-sensing neural network (GraSSNet) for multivariate time-series classification of noisy and highly-imbalanced soft-sensing data. The proposed GraSSNet is able to 1) capture the inter-series and intra-series dependencies jointly in the spectral domain; 2) exploit the label correlations by superimposing label graph that built from statistical co-occurrence information; 3) learn features with attention mechanism from both textual and numerical domain; and 4) leverage unlabeled data and mitigate data imbalance by semi-supervised learning. Comparative studies with other commonly used classifiers are carried out on Seagate soft sensing data, and the experimental results validate the competitive performance of our proposed method.


title: Soft-Sensing ConFormer: A Curriculum Learning-based Convolutional Transformer

abstract: Over the last few decades, modern industrial processes have investigated several cost-effective methodologies to improve the productivity and yield of semiconductor manufacturing. While playing an essential role in facilitating real-time monitoring and control, the data-driven soft-sensors in industries have provided a competitive edge when augmented with deep learning approaches for wafer fault-diagnostics. Despite the success of deep learning methods across various domains, they tend to suffer from bad performance on multi-variate soft-sensing data domains. To mitigate this, we propose a soft-sensing ConFormer (CONvolutional transFORMER) for wafer fault-diagnostic classification task which primarily consists of multi-head convolution modules that reap the benefits of fast and light-weight operations of convolutions, and also the ability to learn the robust representations through multi-head design alike transformers. Another key issue is that traditional learning paradigms tend to suffer from low performance on noisy and highly-imbalanced soft-sensing data. To address this, we augment our soft-sensing ConFormer model with a curriculum learning-based loss function, which effectively learns easy samples in the early phase of training and difficult ones later. To further demonstrate the utility of our proposed architecture, we performed extensive experiments on various toolsets of Seagate Technology's wafer manufacturing process which are shared openly along with this work. To the best of our knowledge, this is the first time that curriculum learning-based soft-sensing ConFormer architecture has been proposed for soft-sensing data and our results show strong promise for future use in soft-sensing research domain.


title: Soft Sensing Model Visualization: Fine-tuning Neural Network from What Model Learned

abstract: The growing availability of the data collected from smart manufacturing is changing the paradigms of production monitoring and control. The increasing complexity and content of the wafer manufacturing process in addition to the time-varying unexpected disturbances and uncertainties, make it infeasible to do the control process with model-based approaches. As a result, data-driven soft-sensing modeling has become more prevalent in wafer process diagnostics. Recently, deep learning has been utilized in soft sensing system with promising performance on highly nonlinear and dynamic time-series data. Despite its successes in soft-sensing systems, however, the underlying logic of the deep learning framework is hard to understand. In this paper, we propose a deep learning-based model for defective wafer detection using a highly imbalanced dataset. To understand how the proposed model works, the deep visualization approach is applied. Additionally, the model is then fine-tuned guided by the deep visualization. Extensive experiments are performed to validate the effectiveness of the proposed system. The results provide an interpretation of how the model works and an instructive fine-tuning method based on the interpretation.


title: Positroids are 3-colorable

abstract: We show that every positroid of rank $r \geq 3$ has a good coline. Using the definition of the chromatic number of oriented matroid introduced by J.\ Ne\v{s}et\v{r}il, R.\ Nickel, and W.~Hochst\"{a}ttler, this shows that every orientation of a positroid is 3-colorable.


title: Exact energy-eigenstates of the Coulomb-Stark Hamiltonian

abstract: An approximation-free, numerically efficient algorithm is presented for the Hamiltonian eigen-states of the Stark-Hydrogen problem describing a quantum particle exposed to the central Coulomb force and a homogeneous external field. As an example of application in a state-expansion with continuous energy, we calculate the time-dependent wavefunction of an electron tunneling from a hydrogen atom suddenly exposed to an external electric field.


title: Bayesian Knockoff Generators for Robust Inference Under Complex Data Structure

abstract: The recent proliferation of medical data, such as genetics and electronic health records (EHR), offers new opportunities to find novel predictors of health outcomes. Presented with a large set of candidate features, interest often lies in selecting the ones most likely to be predictive of an outcome for further study such that the goal is to control the false discovery rate (FDR) at a specified level. Knockoff filtering is an innovative strategy for FDR-controlled feature selection. But, existing knockoff methods make strong distributional assumptions that hinder their applicability to real world data. We propose Bayesian models for generating high quality knockoff copies that utilize available knowledge about the data structure, thus improving the resolution of prognostic features. Applications to two feature sets are considered: those with categorical and/or continuous variables possibly having a population substructure, such as in EHR; and those with microbiome features having a compositional constraint and phylogenetic relatedness. Through simulations and real data applications, these methods are shown to identify important features with good FDR control and power.


title: Expanding the Reach of Gravitational Wave Astronomy to the Edge of the Universe: The Gravitational-Wave International Committee Study Reports on Next Generation Ground-based Gravitational-Wave Observatories

abstract: The first direct detection of gravitational waves emitted from a pair of merging black holes in 2015 has been heralded as one of most significant scientific breakthroughs in physics and astronomy of the 21st century. Motivated by the tremendous scientific opportunities now opened by gravitational-wave observatories and recognizing that to fully exploit the new field will require new observatories that may take 15 to 20 years from conception until operations begin, the Gravitational Wave International Committee (GWIC) convened a subcommittee to examine the path to build and operate a network of future ground-based observatories, capable of extending the observational GW horizon well beyond that currently attainable with the current generation of detectors. This report is the first in a six part series of reports by the GWIC 3G Subcommittee: i) Expanding the Reach of Gravitational Wave Observatories to the Edge of the Universe (this report), ii) The Next Generation Global Gravitational Wave Observatory: The Science Book, iii) 3G R&D: R&D for the Next Generation of Ground-based Gravitational Wave Detectors, iv) Gravitational Wave Data Analysis: Computing Challenges in the 3G Era, v) Future Ground-based Gravitational-wave Observatories: Synergies with Other Scientific Communities, and vi) An Exploration of Possible Governance Models for the Future Global Gravitational-Wave Observatory Network.


abstract: Cyber infrastructure will be a critical consideration in the development of next generation gravitational-wave detectors. The demand for data analysis computing in the 3G era will be driven by the high number of detections as well as the expanded search parameter space for compact astrophysical objects and the subsequent parameter estimation follow-up required to extract the nature of the sources. Additionally, there will be an increased need to develop appropriate and scalable computing cyberinfrastructure, including data access and transfer protocols, and storage and management of software tools, that have sustainable development, support, and management processes. This report identifies the major challenges and opportunities facing 3G gravitational-wave observatories and presents recommendations for addressing them. This report is the fourth in a six part series of reports by the GWIC 3G Subcommittee: i) Expanding the Reach of Gravitational Wave Observatories to the Edge of the Universe, ii) The Next Generation Global Gravitational Wave Observatory: The Science Book, iii) 3G R&D: R&D for the Next Generation of Ground-based Gravitational Wave Detectors, iv) Gravitational Wave Data Analysis: Computing Challenges in the 3G Era (this report), v) Future Ground-based Gravitational-wave Observatories: Synergies with Other Scientific Communities, and vi) An Exploration of Possible Governance Models for the Future Global Gravitational-Wave Observatory Network.


title: Future Ground-Based Gravitational-Wave Observatories: Synergies with Other Scientific Communities

abstract: Planning for the development of a 3rd generation global gravitational-wave detector array is a multifaceted and complex effort that will necessarily need a high level of community input. Interfacing to extant and new stakeholders in the broader scientific constituencies is necessary to keep them aware of the activities taking place in the ground-based gravitational-wave community and receive input to inform and evolve the planning. In this report, we present the approaches GWIC and gravitational-wave collaborations and projects should consider taking to engage with broader community. This report is the fifth in a six part series of reports by the GWIC 3G Subcommittee: i) Expanding the Reach of Gravitational Wave Observatories to the Edge of the Universe, ii) The Next Generation Global Gravitational Wave Observatory: The Science Book, iii) 3G R&D: R&D for the Next Generation of Ground-based Gravitational Wave Detectors, iv) Gravitational Wave Data Analysis: Computing Challenges in the 3G Era, v) Future Ground-based Gravitational-wave Observatories: Synergies with Other Scientific Communities (this report), and vi) An Exploration of Possible Governance Models for the Future Global Gravitational-Wave Observatory Network.


abstract: The construction of a global network of detectors is the cornerstone to scientific success for 3rd generation gravitational wave astronomy. If carried out with a vision to the future, the third generation implementation, in its infrastructure, technology base and governance can provide the point of departure for subsequent developments beyond the third generation, a path to the long and revolutionary future of gravitational wave astrophysics. The goal of this document is to lay out considerations that influence the optimal choice of governance and to lay out a possible path that can lead the community to an optimal governance model. This report is the six in a six part series of reports by the GWIC 3G Subcommittee: i) Expanding the Reach of Gravitational Wave Observatories to the Edge of the Universe, ii) The Next Generation Global Gravitational Wave Observatory: The Science Book, iii) 3G R&D: R&D for the Next Generation of Ground-based Gravitational Wave Detectors, iv) Gravitational Wave Data Analysis: Computing Challenges in the 3G Era, v) Future Ground-based Gravitational-wave Observatories: Synergies with Other Scientific Communities, and vi) An Exploration of Possible Governance Models for the Future Global Gravitational-Wave Observatory Network (this report).


abstract: The next generation of ground-based gravitational-wave detectors will observe coalescences of black holes and neutron stars throughout the cosmos, thousands of them with exceptional fidelity. The Science Book is the result of a 3-year effort to study the science capabilities of networks of next generation detectors. Such networks would make it possible to address unsolved problems in numerous areas of physics and astronomy, from Cosmology to Beyond the Standard Model of particle physics, and how they could provide insights into workings of strongly gravitating systems, astrophysics of compact objects and the nature of dense matter. It is inevitable that observatories of such depth and finesse will make new discoveries inaccessible to other windows of observation. In addition to laying out the rich science potential of the next generation of detectors, this report provides specific science targets in five different areas in physics and astronomy and the sensitivity requirements to accomplish those science goals. This report is the second in a six part series of reports by the GWIC 3G Subcommittee: i) Expanding the Reach of Gravitational Wave Observatories to the Edge of the Universe, ii) The Next Generation Global Gravitational Wave Observatory: The Science Book (this report), iii) 3G R&D: R&D for the Next Generation of Ground-based Gravitational Wave Detectors, iv) Gravitational Wave Data Analysis: Computing Challenges in the 3G Era, v) Future Ground-based Gravitational-wave Observatories: Synergies with Other Scientific Communities, and vi) An Exploration of Possible Governance Models for the Future Global Gravitational-Wave Observatory Network.


title: 3G R&D: R&D for the Next Generation of Ground-based Gravitational-wave Detectors

abstract: To deliver on the promise of next generation gravitational-wave observatories, a sustained and coordinated detector research and development program is required. This report examines in detail the wide range of nearer- and longer-term detector R&D programs needed for next generation GW detectors commensurate with the key science targets presented in "The Next Generation Global Gravitational Wave Observatory: The Science Book", including considerations of site selection and large-scale vacuum infrastructure. The report makes a series of detailed recommendations on the needed advances in detector technology and the timescales needed to achieve those advances. It also identifies areas where larger-scale globally coordinated R&D efforts will be critical to ensuring success while minimizing costs. This report is the third in a six part series of reports by the GWIC 3G Subcommittee: i) Expanding the Reach of Gravitational Wave Observatories to the Edge of the Universe, ii) The Next Generation Global Gravitational Wave Observatory: The Science Book, iii) 3G R&D: R&D for the Next Generation of Ground-based Gravitational Wave Detectors (this report), iv) Gravitational Wave Data Analysis: Computing Challenges in the 3G Era, v) Future Ground-based Gravitational-wave Observatories: Synergies with Other Scientific Communities, and vi) An Exploration of Possible Governance Models for the Future Global Gravitational-Wave Observatory Network.


title: Drying-induced stresses in poroelastic drops on rigid substrates

abstract: We develop a theory for drying-induced stresses in sessile, poroelastic drops undergoing evaporation on rigid surfaces. Using a lubrication-like approximation, the governing equations of three-dimensional nonlinear poroelasticity are reduced to a single thin-film equation for the drop thickness. We find that thin drops experience compressive elastic stresses but the total in-plane stresses are tensile. The mechanical response of the drop is dictated by the initial profile of the solid skeleton, which controls the in-plane deformation, the dominant components of elastic stress, and sets a limit on the depth of delamination that can potentially occur. Our theory suggests that the alignment of desiccation fractures in colloidal drops is selected by the shape of the drop at the point of gelation. We propose that the emergence of three distinct fracture patterns in dried blood drops is a consequence of a non-monotonic drop profile at gelation. We also show that depletion fronts, which separate wet and dry solid, can invade the drop from the contact line and localise the generation of mechanical stress during drying. Finally, the finite element method is used to explore the stress profiles in drops with large contact angles.


title: On Affine Hilbert Functions of Unions of Layers in Finite Grids

abstract: The affine Hilbert function is a classical algebraic object that has been central, among other tools, to the development of the polynomial method in combinatorics. Owing to its concrete connections with Gr\"obner basis theory, as well as its applicability in several areas like computational complexity, combinatorial geometry, and coding theory, an important line of enquiry is to understand the affine Hilbert function of structured sets of points in the affine space. In this work, we determine the affine Hilbert function (over the reals) of arbitrary unions of layers of points in a uniform grid (a finite grid with the component sets having equispaced points), where each layer of points is determined by a fixed sum of components for all the points. This extends a result of Bernasconi and Egidi (Inf. Comput. 1999) from the Boolean cube setting to the uniform grid setting. Our proofs follow a similar outline as that of Bernasconi and Egidi. However, there are two bottlenecks that arise in the uniform grid setting. We resolve these by using (i) a classical fact that a symmetric Jordan basis exists for the function space on a uniform grid, and (ii) an extension to multisets of an algebraic interpretation by Friedl and R\'onyai (Discrete Math. 2003) of the notion of order shattering. The affine Hilbert function is, in fact, a stronger notion than the finite-degree Zariski closure, which is yet another important tool in the polynomial method toolkit. We conclude by giving an alternative proof of a combinatorial characterization of a variant of finite-degree Zariski closures, for unions of layers in uniform grids, obtained in an earlier work of the author (arXiv, 2021).


title: Learning Online for Unified Segmentation and Tracking Models

abstract: Tracking requires building a discriminative model for the target in the inference stage. An effective way to achieve this is online learning, which can comfortably outperform models that are only trained offline. Recent research shows that visual tracking benefits significantly from the unification of visual tracking and segmentation due to its pixel-level discrimination. However, it imposes a great challenge to perform online learning for such a unified model. A segmentation model cannot easily learn from prior information given in the visual tracking scenario. In this paper, we propose TrackMLP: a novel meta-learning method optimized to learn from only partial information to resolve the imposed challenge. Our model is capable of extensively exploiting limited prior information hence possesses much stronger target-background discriminability than other online learning methods. Empirically, we show that our model achieves state-of-the-art performance and tangible improvement over competing models. Our model achieves improved average overlaps of66.0%,67.1%, and68.5% in VOT2019, VOT2018, and VOT2016 datasets, which are 6.4%,7.3%, and6.4% higher than our baseline. Code will be made publicly available.


title: A Central Difference Graph Convolutional Operator for Skeleton-Based Action Recognition

abstract: This paper proposes a new graph convolutional operator called central difference graph convolution (CDGC) for skeleton based action recognition. It is not only able to aggregate node information like a vanilla graph convolutional operation but also gradient information. Without introducing any additional parameters, CDGC can replace vanilla graph convolution in any existing Graph Convolutional Networks (GCNs). In addition, an accelerated version of the CDGC is developed which greatly improves the speed of training. Experiments on two popular large-scale datasets NTU RGB+D 60 & 120 have demonstrated the efficacy of the proposed CDGC. Code is available at https://github.com/iesymiao/CD-GCN.


title: Dynamics of photosynthetic light harvesting systems interacting with N-photon Fock states

abstract: We develop a method to simulate the excitonic dynamics of realistic photosynthetic light harvesting systems including non-Markovian coupling to phonon degrees of freedom, under excitation by N-photon Fock state pulses. This method combines the input-output formalism and the hierarchical equations of motion (HEOM) formalism into a double hierarchy of coupled linear equations in density matrices. We show analytically that, under weak field excitation relevant to natural photosynthesis conditions, an N-photon Fock state input and a corresponding coherent state input give rise to equal density matrices in the excited manifold. However, an important difference is that an N-photon Fock state input has no off-diagonal coherence between the ground and excited subspaces, in contrast with the coherences created by a coherent state input. We derive expressions for the probability to absorb a single Fock state photon, with or without the influence of phonons. For short pulses (or equivalently, wide bandwidth pulses), we show that the absorption probability has a universal behavior that depends only upon a system-dependent effective energy spread parameter {\Delta} and an exciton-light coupling constant {\Gamma}. This holds for a broad range of chromophore systems and for a variety of pulse shapes. We also analyse the absorption probability in the opposite long pulses (narrow bandwidth) regime. We also derive an expression for the long time emission rate in the presence of phonons and use it to study the difference between collective versus independent emission. Finally, we present a numerical simulation for the LHCII monomer (14-mer) system under single photon excitation that illustrates the use of the double hierarchy for calculation of Fock state excitation of a light harvesting system including chromophore coupling to a non-Markovian phonon bath.


title: A discrete complement of Lyapunov's inequality and its information theoretic consequences

abstract: We establish a reversal of Lyapunov's inequality for monotone log-concave sequences, settling a conjecture of Havrilla-Tkocz and Melbourne-Tkocz. A strengthened version of the same conjecture is disproved through counter example. We also derive several information theoretic inequalities as consequences. In particular sharp bounds are derived for the varentropy, R\'enyi entropies, and the concentration of information of monotone log-concave random variables. Moreover, the majorization approach utilized in the proof of the main theorem, is applied to derive analogous information theoretic results in the symmetric setting, where the Lyapunov reversal is known to fail.


title: A Hybrid EM Algorithm for Linear Two-Way Interactions with Missing Data

abstract: We study an EM algorithm for estimating product-term regression models with missing data. The study of such problems in the likelihood tradition has thus far been restricted to an EM algorithm method using full numerical integration. However, under most missing data patterns, we show that this problem can be solved analytically, and numerical approximations are only needed under specific conditions. Thus we propose a hybrid EM algorithm, which uses analytic solutions when available and approximate solutions only when needed. The theoretical framework of our algorithm is described herein, along with two numerical experiments using both simulated and real data. We show that our algorithm confers higher accuracy to the estimation process, relative to the existing full numerical integration method. We conclude with a discussion of applications, extensions, and topics of further research.


title: Quasifibrations of Graphs to Find Symmetries in Biological Networks

abstract: A fibration of graphs is an homomorphism that is a local isomorphism of in-neighbourhoods, much in the same way a covering projection is a local isomorphism of neighbourhoods. Recently, it has been shown that graph fibrations are useful tools to uncover symmetries and synchronization patterns in biological networks ranging from gene, protein,and metabolic networks to the brain. However, the inherent incompleteness and disordered nature of biological data precludes the application of the definition of fibration as it is; as a consequence, also the currently known algorithms to identify fibrations fail in these domains. In this paper, we introduce and develop systematically the theory of quasifibrations which attempts to capture more realistic patterns of almost-synchronization of units in biological networks. We provide an algorithmic solution to the problem of finding quasifibrations in networks where the existence of missing links and variability across samples preclude the identification of perfect symmetries in the connectivity structure. We test the algorithm against other strategies to repair missing links in incomplete networks using real connectome data and synthetic networks. Quasifibrations can be applied to reconstruct any incomplete network structure characterized by underlying symmetries and almost synchronized clusters.

categories: q-bio.QM cs.IT math.IT math.OC physics.bio-ph physics.data-an


title: Equilibrium conformations and surface charge regulation of spherical polymer brushes in stretched regimes

abstract: In the present work, we study the equilibrium conformations of linear polyelectrolytes tethered onto a spherical, oppositely charged core in equilibrium with an ionic reservoir of fixed concentration. Particular focus is placed on the situation of stretched chains, where the monomer concentration is known to display an inverse square-law decay far away from the spherical surface, which is then further extrapolated all the way down to the grafting core. While the equilibrium distributions of mobile ions are computed in the framework of a classical Density Functional Theory (cDFT) that incorporates both their size and electrostatic correlations within the grafted polyelectrolyte, the equilibrium configuration of the latter is described by its averaged radius of gyration, which is taken as a variational parameter that guarantees mechanical equilibrium across the polymer-solvent interface. The averaged particle size is then analyzed over a wide range of polymerization degrees, ionic concentrations and functionality of the polymer backbones. Two distinct regimes can be identified: at high ionic strengths, swelling of the grafted polymers is dominated by ionic entropic contribution as well as polymer size effects, whereas at low ionic concentrations a balance between electrostatic and entropic effects is the main driven mechanism for particle stretching. Using Monte Carlo simulations, we then proceed to investigate the effects of charge regulation when the brush core is further decorated with active functional sites randomly distributed over its surface, which act as receptors onto which dissolved acidic ions can be adsorbed.


title: LoMEF: A Framework to Produce Local Explanations for Global Model Time Series Forecasts

abstract: Global Forecasting Models (GFM) that are trained across a set of multiple time series have shown superior results in many forecasting competitions and real-world applications compared with univariate forecasting approaches. One aspect of the popularity of statistical forecasting models such as ETS and ARIMA is their relative simplicity and interpretability (in terms of relevant lags, trend, seasonality, and others), while GFMs typically lack interpretability, especially towards particular time series. This reduces the trust and confidence of the stakeholders when making decisions based on the forecasts without being able to understand the predictions. To mitigate this problem, in this work, we propose a novel local model-agnostic interpretability approach to explain the forecasts from GFMs. We train simpler univariate surrogate models that are considered interpretable (e.g., ETS) on the predictions of the GFM on samples within a neighbourhood that we obtain through bootstrapping or straightforwardly as the one-step-ahead global black-box model forecasts of the time series which needs to be explained. After, we evaluate the explanations for the forecasts of the global models in both qualitative and quantitative aspects such as accuracy, fidelity, stability and comprehensibility, and are able to show the benefits of our approach.


title: Refactoring for Reuse: An Empirical Study

abstract: Refactoring is the de-facto practice to optimize software health. While several studies propose refactoring strategies to optimize software design through applying design patterns and removing design defects, little is known about how developers actually refactor their code to improve its reuse. Therefore, we extract, from 1,828 open-source projects, a set of refactorings that were intended to improve the software reusability. We analyze the impact of reusability refactorings on the state-of-the-art reusability metrics, and we compare the distribution of reusability refactoring types, with the distribution of the remaining mainstream refactorings. Overall, we found that the distribution of refactoring types, applied in the context of reusability, is different from the distribution of refactoring types in mainstream development. In the refactorings performed to improve reusability, source files are subject to more design-level types of refactorings. Reusability refactorings significantly impact, high-level code elements, such as packages, classes, and methods, while typical refactorings, impact all code elements, including identifiers, and parameters. These findings provide practical insights into the current practice of refactoring in the context of code reuse involving the act of refactoring.


title: A hybrid-mixed finite element method for single-phase Darcy flow in fractured porous media

abstract: We present a hybrid-mixed finite element method for a novel hybrid-dimensional model of single-phase Darcy flow in a fractured porous media. In this model, the fracture is treated as an $(d-1)$-dimensional interface within the $d$-dimensional fractured porous domain, for $d=2, 3$. Two classes of fracture are distinguished based on the permeability magnitude ratio between the fracture and its surrounding medium: when the permeability in the fracture is (significantly) larger than in its surrounding medium, it is considered as a {\it conductive} fracture; when the permeability in the fracture is (significantly) smaller than in its surrounding medium, it is considered as a {\it blocking} fracture. The conductive fractures are treated using the classical hybrid-dimensional approach of the interface model where pressure is assumed to be continuous across the fracture interfaces, while the blocking fractures are treated using the recent Dirac-$\delta$ function approach where normal component of Darcy velocity is assumed to be continuous across the interface. Due to the use of Dirac-$\delta$ function approach for the blocking fractures, our numerical scheme allows for nonconforming meshes with respect to the blocking fractures. This is the major novelty of our model and numerical discretization. Moreover, our numerical scheme produces locally conservative velocity approximations and leads to a symmetric positive definite linear system involving pressure degrees of freedom on the mesh skeleton only. The performance of the proposed method is demonstrated by various benchmark test cases in both two- and three-dimensions. Numerical results indicate that the proposed scheme is highly competitive with existing methods in the literature.


title: Fault Diagnosis of Nonlinear Systems Using a Hybrid-Degree Dual Cubature-based Estimation Scheme

abstract: In this paper, a novel hybrid-degree dual estimation approach based on cubature rules and cubature-based nonlinear filters is proposed for fault diagnosis of nonlinear systems through simultaneous state and time-varying parameter estimation. Our proposed dual nonlinear filtering scheme is developed based on case-dependent cubature rules that are motivated by the following observations and facts, namely (i) dynamic characteristics of nonlinear system states and parameters generally are distinct and posses different degrees of complexities, and (ii) performance of cubature rules depend on the system dynamics and vary due to handling of high-dimensional integrations approximations. For improving the robustness capability of our proposed methodologies modified cubature point propagation method is incorporated. The performance of our proposed dual estimation strategy is demonstrated and evaluated by application to a nonlinear gas turbine engine for addressing the component fault diagnosis problem within an integrated fault detection, isolation and identification framework. Robustness analysis is implemented to verify the capability of our proposed approaches to deal with parametric uncertainties and unmodeled dynamics. Extensive simulation case studies and discussions with respect to component fouling, erosion or abrupt faults are provided to substantiate and justify the superiority of our proposed fault diagnosis methodology when compared to other well-known alternative diagnostic techniques such as the Unscented Kalman Filters (UKF) and Particle Filters (PF) that are commonly available in the literature.


title: Understanding and Assessment of Mission-Centric Key Cyber Terrains for joint Military Operations

abstract: Since the cyberspace consolidated as fifth warfare dimension, the different actors of the defense sector began an arms race toward achieving cyber superiority, on which research, academic and industrial stakeholders contribute from a dual vision, mostly linked to a large and heterogeneous heritage of developments and adoption of civilian cybersecurity capabilities. In this context, augmenting the conscious of the context and warfare environment, risks and impacts of cyber threats on kinetic actuations became a critical rule-changer that military decision-makers are considering. A major challenge on acquiring mission-centric Cyber Situational Awareness (CSA) is the dynamic inference and assessment of the vertical propagations from situations that occurred at the mission supportive Information and Communications Technologies (ICT), up to their relevance at military tactical, operational and strategical views. In order to contribute on acquiring CSA, this paper addresses a major gap in the cyber defence state-of-the-art: the dynamic identification of Key Cyber Terrains (KCT) on a mission-centric context. Accordingly, the proposed KCT identification approach explores the dependency degrees among tasks and assets defined by commanders as part of the assessment criteria. These are correlated with the discoveries on the operational network and the asset vulnerabilities identified thorough the supported mission development. The proposal is presented as a reference model that reveals key aspects for mission-centric KCT analysis and supports its enforcement and further enforcement by including an illustrative application case.


title: Optimization Framework for Splitting DNN Inference Jobs over Computing Networks

abstract: Ubiquitous artificial intelligence (AI) is considered one of the key services in 6G systems. AI services typically rely on deep neural network (DNN) requiring heavy computation. Hence, in order to support ubiquitous AI, it is crucial to provide a solution for offloading or distributing computational burden due to DNN, especially at end devices with limited resources. We develop an optimization framework for assigning the computation tasks of DNN inference jobs to computing resources in the network, so as to reduce the inference latency. To this end, we propose a layered graph model with which simple conventional routing jointly solves the problem of selecting nodes for computation and paths for data transfer between nodes. We show that using our model, the existing approaches to splitting DNN inference jobs can be equivalently reformulated as a routing problem that possesses better numerical properties. We also apply the proposed framework to derive algorithms for minimizing the end-to-end inference latency. We show through numerical evaluations that our new formulation can find a solution for DNN inference job distribution much faster than the existing formulation, and that our algorithms can select computing nodes and data paths adaptively to the computational attributes of given DNN inference jobs, so as to reduce the end-to-end latency.


title: GPMR: An Iterative Method for Unsymmetric Partitioned Linear Systems

abstract: We introduce an iterative method named GPMR for solving 2x2 block unsymmetric linear systems. GPMR is based on a new process that reduces simultaneously two rectangular matrices to upper Hessenberg form and that is closely related to the block-Arnoldi process. GPMR is tantamount to Block-GMRES with two right-hand sides in which the two approximate solutions are summed at each iteration, but requires less storage and work per iteration. We compare the performance of GPMR with GMRES and Block-GMRES on linear systems from the SuiteSparse Matrix Collection. In our experiments, GPMR terminates significantly earlier than GMRES on a residual-based stopping condition with an improvement ranging from around 10% up to 50% in terms of number of iterations. We also illustrate by experiment that GPMR appears more resilient to loss of orthogonality than Block-GMRES.


title: Nii: a Bayesian orbit retrieval code applied to differential astrometry

abstract: Here we present an open source Python-based Bayesian orbit retrieval code (Nii) that implements an automatic parallel tempering Markov chain Monte Carlo (APT-MCMC) strategy. Nii provides a module to simulate the observations of a space-based astrometry mission in the search for exoplanets, a signal extraction process for differential astrometric measurements using multiple reference stars, and an orbital parameter retrieval framework using APT-MCMC. We further verify the orbit retrieval ability of the code through two examples corresponding to a single-planet system and a dual-planet system. In both cases, efficient convergence on the posterior probability distribution can be achieved. Although this code specifically focuses on the orbital parameter retrieval problem of differential astrometry, Nii can also be widely used in other Bayesian analysis applications.


title: Leveraging Unsupervised Image Registration for Discovery of Landmark Shape Descriptor

abstract: In current biological and medical research, statistical shape modeling (SSM) provides an essential framework for the characterization of anatomy/morphology. Such analysis is often driven by the identification of a relatively small number of geometrically consistent features found across the samples of a population. These features can subsequently provide information about the population shape variation. Dense correspondence models can provide ease of computation and yield an interpretable low-dimensional shape descriptor when followed by dimensionality reduction. However, automatic methods for obtaining such correspondences usually require image segmentation followed by significant preprocessing, which is taxing in terms of both computation as well as human resources. In many cases, the segmentation and subsequent processing require manual guidance and anatomy specific domain expertise. This paper proposes a self-supervised deep learning approach for discovering landmarks from images that can directly be used as a shape descriptor for subsequent analysis. We use landmark-driven image registration as the primary task to force the neural network to discover landmarks that register the images well. We also propose a regularization term that allows for robust optimization of the neural network and ensures that the landmarks uniformly span the image domain. The proposed method circumvents segmentation and preprocessing and directly produces a usable shape descriptor using just 2D or 3D images. In addition, we also propose two variants on the training loss function that allows for prior shape information to be integrated into the model. We apply this framework on several 2D and 3D datasets to obtain their shape descriptors, and analyze their utility for various applications.


title: Nonperturbative electromagnetic nonlinearities, n-photon reflectors, and Fock-state lasers based on deep-strong coupling of light and matter

abstract: Light and matter can now interact in a regime where their coupling is stronger than their bare energies. This deep-strong coupling (DSC) regime of quantum electrodynamics promises to challenge many conventional assumptions about the physics of light and matter. Here, we show how light and matter interactions in this regime give rise to electromagnetic nonlinearities dramatically different from those of naturally existing materials. Excitations in the DSC regime act as photons with a linear energy spectrum up to a critical excitation number, after which, the system suddenly becomes strongly anharmonic, thus acting as an effective intensity-dependent nonlinearity of an extremely high order. We show that this behavior allows for N-photon blockade (with $N \gg 1$), enabling qualitatively new kinds of quantum light sources. For example, this nonlinearity forms the basis for a new type of gain medium, which when integrated into a laser or maser, produces large Fock states (rather than coherent states). Such Fock states could in principle have photon numbers orders of magnitude larger than any realized previously, and would be protected from dissipation by a new type of equilibrium between nonlinear gain and linear loss. We discuss paths to experimental realization of the effects described here.


title: An explicit predictor/multicorrector time marching with automatic adaptivity for finite-strain elastodynamics

abstract: We propose a time-adaptive predictor/multi-corrector method to solve hyperbolic partial differential equations, based on the generalized-$\alpha$ scheme that provides user-control on the numerical dissipation and second-order accuracy in time. Our time adaptivity uses an error estimation that exploits the recursive structure of the variable updates. The predictor/multicorrector method explicitly updates the equation system but computes the residual of the system implicitly. We analyze the method's stability and describe how to determine the parameters that ensure high-frequency dissipation and accurate low-frequency approximation. Subsequently, we solve a linear wave equation, followed by non-linear finite strain deformation problems with different boundary conditions. Thus, our method is a straightforward, stable and computationally efficient approach to simulate real-world engineering problems. Finally, to show the performance of our method, we provide several numerical examples in two and three dimensions. These challenging tests demonstrate that our predictor/multicorrector scheme dynamically adapts to sudden energy releases in the system, capturing impacts and boundary shocks. The method efficiently and stably solves dynamic equations with consistent and under-integrated mass matrices conserving the linear and angular momenta as well as the system's energy for long-integration times.


title: Ab initio calculations for void swelling bias in $\alpha$ and $\delta$-plutonium

abstract: Void swelling can develop in materials under persistent irradiation when non-equilibrium vacancy and self-interstitial populations migrate under sufficiently asymmetric interaction biases. In conventional metals, the propensity is determined to first approximation by comparing point-defect relaxation strains. We thus present DFT-based calculations of structures and formation energies and volumes of point defects in the $\alpha$ and the $\delta$-phases of plutonium. Our results show that lattice defects in $\delta$-Pu, in contrast to most fcc metals, have surprisingly small formation volumes. Equally unexpected are the large defect formation volumes found in the low-symmetry $\alpha$-Pu phase. Both these unusual properties can be satisfactorily explained from defect-induced spin/orbital moment formation and destruction in the Pu phases. When we use the calculated defect properties to estimate the classic void swelling bias in each of the phases, we find it to be unusually small in $\delta$-Pu, but likely much larger in $\alpha$-Pu. Hence, swelling rates and mechanisms can diverge dramatically between the different phases of Pu. Especially in the transient regime before formation of large defect clusters, the swelling rate of $\alpha$-Pu can reliably be expected to be much larger than $\delta$-Pu. However, accurate forecasts over longer times will require the conventional void-swelling theory to be modified to handle the complexities presented by the different Pu phases. As a case in point, we show possible anomalous temperature dependence of vacancy properties in $\delta$-Pu, caused by entropic contributions from defect-induced spin-lattice fluctuations. Such complications may affect defect-defect interactions and thus alter the void swelling bias.


title: Identities for a parametric Weyl algebra over a ring

abstract: In 2013 Benkart, Lopes and Ondrus introduced and studied in a series of papers the infinite-dimensional unital associative algebra $\A_h$ generated by elements $x,y,$ which satisfy the relation $yx-xy=h$ for some $0\neq h\in \FF[x]$. We generalize this construction to $\A_h(\B)$ by working over the fixed $\FF$-algebra $\B$ instead of $\FF$. We describe the polynomial identities for $\A_h(\B)$ over the infinite field $\FF$ in case $h\in\B[x]$ satisfies certain restrictions.


title: Milnor's triple linking number and Gauss diagram formulas of 3-bouquet graphs

abstract: In this paper, we introduce two functions such that the subtraction corresponds to the Milnor's triple linking number; the addition obtains a new integer-valued link homotopy invariant of $3$-component links. We also have found a series of integer-valued invariants derived from four terms whose sum equals the Milnor's triple linking number. We apply this structure to give invariants of $3$-bouquet graphs.


title: HydraGAN A Multi-head, Multi-objective Approach to Synthetic Data Generation

abstract: Synthetic data generation overcomes limitations of real-world machine learning. Traditional methods are valuable for augmenting costly datasets but only optimize one criterion: realism. In this paper, we tackle the problem of generating synthetic data that optimize multiple criteria. This goal is necessary when real data are replaced by synthetic for privacy preservation. We introduce HydraGAN, a new approach to synthetic data generation that introduces multiple generator and discriminator agents into the system. The multi-agent GAN optimizes the goal of privacy-preservation as well as data realism. To facilitate multi-agent training, we adapt game-theoretic principles to offer equilibrium guarantees. We observe that HydraGAN outperforms baseline methods for three datasets for multiple criteria of maximizing data realism, maximizing model accuracy, and minimizing re-identification risk.


title: Robust Multi-Robot Trajectory Optimization Using Alternating Direction Method of Multiplier

abstract: We propose a variant of alternating direction method of multiplier (ADMM) to solve constrained trajectory optimization problems. Our ADMM framework breaks a joint optimization into small sub-problems, leading to a low iteration cost and decentralized parameter updates. Starting from a collision-free initial trajectory, our method inherits the theoretical properties of primal interior point method (P-IPM), i.e., guaranteed collision avoidance and homotopy preservation throughout optimization, while being orders of magnitude faster. We have analyzed the convergence and evaluated our method for time-optimal multi-UAV trajectory optimizations and simultaneous goal-reaching of multiple robot arms, where we take into consider kinematics-, dynamics-limits, and homotopy-preserving collision constraints. Our method highlights an order of magnitude's speedup, while generating trajectories of comparable qualities as state-of-the-art P-IPM solver.


title: Delay-Oriented Distributed Scheduling Using Graph Neural Networks

abstract: In wireless multi-hop networks, delay is an important metric for many applications. However, the max-weight scheduling algorithms in the literature typically focus on instantaneous optimality, in which the schedule is selected by solving a maximum weighted independent set (MWIS) problem on the interference graph at each time slot. These myopic policies perform poorly in delay-oriented scheduling, in which the dependency between the current backlogs of the network and the schedule of the previous time slot needs to be considered. To address this issue, we propose a delay-oriented distributed scheduler based on graph convolutional networks (GCNs). In a nutshell, a trainable GCN module generates node embeddings that capture the network topology as well as multi-step lookahead backlogs, before calling a distributed greedy MWIS solver. In small- to medium-sized wireless networks with heterogeneous transmit power, where a few central links have many interfering neighbors, our proposed distributed scheduler can outperform the myopic schedulers based on greedy and instantaneously optimal MWIS solvers, with good generalizability across graph models and minimal increase in communication complexity.


title: Identification and Adaptive Control of Markov Jump Systems: Sample Complexity and Regret Bounds

abstract: Learning how to effectively control unknown dynamical systems is crucial for intelligent autonomous systems. This task becomes a significant challenge when the underlying dynamics are changing with time. Motivated by this challenge, this paper considers the problem of controlling an unknown Markov jump linear system (MJS) to optimize a quadratic objective. By taking a model-based perspective, we consider identification-based adaptive control for MJSs. We first provide a system identification algorithm for MJS to learn the dynamics in each mode as well as the Markov transition matrix, underlying the evolution of the mode switches, from a single trajectory of the system states, inputs, and modes. Through mixing-time arguments, sample complexity of this algorithm is shown to be $\mathcal{O}(1/\sqrt{T})$. We then propose an adaptive control scheme that performs system identification together with certainty equivalent control to adapt the controllers in an episodic fashion. Combining our sample complexity results with recent perturbation results for certainty equivalent control, we prove that when the episode lengths are appropriately chosen, the proposed adaptive control scheme achieves $\mathcal{O}(\sqrt{T})$ regret, which can be improved to $\mathcal{O}(polylog(T))$ with partial knowledge of the system. Our proof strategy introduces innovations to handle Markovian jumps and a weaker notion of stability common in MJSs. Our analysis provides insights into system theoretic quantities that affect learning accuracy and control performance. Numerical simulations are presented to further reinforce these insights.