title: Voxelmorph++ Going beyond the cranial vault with keypoint supervision and multi-channel instance optimisation

abstract: The majority of current research in deep learning based image registration addresses inter-patient brain registration with moderate deformation magnitudes. The recent Learn2Reg medical registration benchmark has demonstrated that single-scale U-Net architectures, such as VoxelMorph that directly employ a spatial transformer loss, often do not generalise well beyond the cranial vault and fall short of state-of-the-art performance for abdominal or intra-patient lung registration. Here, we propose two straightforward steps that greatly reduce this gap in accuracy. First, we employ keypoint self-supervision with a novel network head that predicts a discretised heatmap and robustly reduces large deformations for better robustness. Second, we replace multiple learned fine-tuning steps by a single instance optimisation with hand-crafted features and the Adam optimiser. Different to other related work, including FlowNet or PDD-Net, our approach does not require a fully discretised architecture with correlation layer. Our ablation study demonstrates the importance of keypoints in both self-supervised and unsupervised (using only a MIND metric) settings. On a multi-centric inspiration-exhale lung CT dataset, including very challenging COPD scans, our method outperforms VoxelMorph by improving nonlinear alignment by 77% compared to 19% - reaching target registration errors of 2 mm that outperform all but one learning methods published to date. Extending the method to semantic features sets new stat-of-the-art performance on inter-subject abdominal CT registration.


title: Local and Global GANs with Semantic-Aware Upsampling for Image Generation

abstract: In this paper, we address the task of semantic-guided image generation. One challenge common to most existing image-level generation methods is the difficulty in generating small objects and detailed local textures. To address this, in this work we consider generating images using local context. As such, we design a local class-specific generative network using semantic maps as guidance, which separately constructs and learns subgenerators for different classes, enabling it to capture finer details. To learn more discriminative class-specific feature representations for the local generation, we also propose a novel classification module. To combine the advantages of both global image-level and local class-specific generation, a joint generation network is designed with an attention fusion module and a dual-discriminator structure embedded. Lastly, we propose a novel semantic-aware upsampling method, which has a larger receptive field and can take far-away pixels that are semantically related for feature upsampling, enabling it to better preserve semantic consistency for instances with the same semantic labels. Extensive experiments on two image generation tasks show the superior performance of the proposed method. State-of-the-art results are established by large margins on both tasks and on nine challenging public benchmarks. The source code and trained models are available at https://github.com/Ha0Tang/LGGAN.


title: Multi-modal Alignment using Representation Codebook

abstract: Aligning signals from different modalities is an important step in vision-language representation learning as it affects the performance of later stages such as cross-modality fusion. Since image and text typically reside in different regions of the feature space, directly aligning them at instance level is challenging especially when features are still evolving during training. In this paper, we propose to align at a higher and more stable level using cluster representation. Specifically, we treat image and text as two "views" of the same entity, and encode them into a joint vision-language coding space spanned by a dictionary of cluster centers (codebook). We contrast positive and negative samples via their cluster assignments while simultaneously optimizing the cluster centers. To further smooth out the learning process, we adopt a teacher-student distillation paradigm, where the momentum teacher of one view guides the student learning of the other. We evaluated our approach on common vision language benchmarks and obtain new SoTA on zero-shot cross modality retrieval while being competitive on various other transfer tasks.


title: Towards Targeted Change Detection with Heterogeneous Remote Sensing Images for Forest Mortality Mapping

abstract: Several generic methods have recently been developed for change detection in heterogeneous remote sensing data, such as images from synthetic aperture radar (SAR) and multispectral radiometers. However, these are not well suited to detect weak signatures of certain disturbances of ecological systems. To resolve this problem we propose a new approach based on image-to-image translation and one-class classification (OCC). We aim to map forest mortality caused by an outbreak of geometrid moths in a sparsely forested forest-tundra ecotone using multisource satellite images. The images preceding and following the event are collected by Landsat-5 and RADARSAT-2, respectively. Using a recent deep learning method for change-aware image translation, we compute difference images in both satellites' respective domains. These differences are stacked with the original pre- and post-event images and passed to an OCC trained on a small sample from the targeted change class. The classifier produces a credible map of the complex pattern of forest mortality.


title: Magnetising the ${\cal N}=4$ Super Yang-Mills plasma

abstract: We investigate the thermodynamics of the anisotropic magnetic $\rm AdS_5$ black brane solution found by D'Hoker and Kraus arXiv:0908.3875. This solution is the gravity dual of a strongly coupled ${\cal N}=4$ Super Yang-Mills plasma in ${\mathbb R}^{3,1}$, with temperature $T$, in the presence of a magnetic field ${\cal B}$. Following the procedure of holographic renormalisation we calculate the Gibbs free energy and the holographic stress tensor of the conformal plasma. We evaluate several thermodynamic quantities including the magnetisation, the anisotropic pressures and the speeds of sound. Our results are consistent with an RG flow from a perturbed $\rm AdS_5$ black brane at small ${\cal B}/T^2$ to a $\rm \bf BTZ \times {\mathbb R}^2$ black brane at large ${\cal B}/T^2$. We also perform a phenomenological analysis where we compare the thermodynamics of a magnetised conformal plasma against the lattice QCD results for the thermodynamics of the magnetised quark-gluon plasma.


title: ERF: Explicit Radiance Field Reconstruction From Scratch

abstract: We propose a novel explicit dense 3D reconstruction approach that processes a set of images of a scene with sensor poses and calibrations and estimates a photo-real digital model. One of the key innovations is that the underlying volumetric representation is completely explicit in contrast to neural network-based (implicit) alternatives. We encode scenes explicitly using clear and understandable mappings of optimization variables to scene geometry and their outgoing surface radiance. We represent them using hierarchical volumetric fields stored in a sparse voxel octree. Robustly reconstructing such a volumetric scene model with millions of unknown variables from registered scene images only is a highly non-convex and complex optimization problem. To this end, we employ stochastic gradient descent (Adam) which is steered by an inverse differentiable renderer. We demonstrate that our method can reconstruct models of high quality that are comparable to state-of-the-art implicit methods. Importantly, we do not use a sequential reconstruction pipeline where individual steps suffer from incomplete or unreliable information from previous stages, but start our optimizations from uniformed initial solutions with scene geometry and radiance that is far off from the ground truth. We show that our method is general and practical. It does not require a highly controlled lab setup for capturing, but allows for reconstructing scenes with a vast variety of objects, including challenging ones, such as outdoor plants or furry toys. Finally, our reconstructed scene models are versatile thanks to their explicit design. They can be edited interactively which is computationally too costly for implicit alternatives.


title: Gaussian quantum estimation of the lossy parameter in a thermal environment

abstract: Lossy bosonic channels play an important role in a number of quantum information tasks, since they well approximate thermal dissipation in an experiment. Here, we characterize their metrological power in the idler-free and entanglement-assisted cases, using respectively single- and two-mode Gaussian states as probes. In the problem of estimating the lossy parameter, we study the energy-constrained quantum Fisher information (QFI) for generic temperature and lossy parameter regimes, showing qualitative behaviours of the optimal probes. We show semi-analytically that the two-mode squeezed-vacuum state optimizes the QFI for any value of the lossy parameter and temperature. We discuss the optimization of the {\it total} QFI, where the number of probes is allowed to vary by keeping the total energy-constrained. In this context, we elucidate the role of the "shadow-effect" for reaching a quantum advantage. We also consider a photon-number normalization for the environment, widely used in the analysis of quantum illumination and quantum reading protocols. With this normalization, we prove that the large bandwidth TMSV state is the optimal probe for any parameter value. Here, the quantum advantage is of at most a factor of $2$, and is reached in the bright environment case for {\it any} lossy parameter values. Finally, we discuss the implications of our results for quantum illumination and quantum reading applications.


title: Fast Bayesian estimation of brain activation with cortical surface and subcortical fMRI data using EM

abstract: Analysis of brain imaging scans is critical to understanding the way the human brain functions, which can be leveraged to treat injuries and conditions that affect the quality of life for a significant portion of the human population. In particular, functional magnetic resonance imaging (fMRI) scans give detailed data on a living subject at high spatial and temporal resolutions. Due to the high cost involved in the collection of these scans, robust methods of analysis are of critical importance in order to produce meaningful inference. Bayesian methods in particular allow for the inclusion of expected behavior from prior study into an analysis, increasing the power of the results while circumventing problems that arise in classical analyses, including the effects of smoothing results and sensitivity to multiple comparison testing corrections. Recent development of a surface-based spatial Bayesian general linear model for cortical surface fMRI (cs-fMRI) data provides the desired power increase in task fMRI data using stochastic partial differential equation (SPDE) priors. This model relies on the computational efficiencies of the integrated nested Laplace approximation (INLA) to perform powerful analyses that have been validated to outperform classical analyses. In this article, we develop an exact Bayesian analysis method for the GLM, employing an expectation-maximization (EM) algorithm to find maximum a posteriori (MAP) estimates of task-based regressors on cs-fMRI and subcortical fMRI data while using minimal computational resources. Our proposed method is compared to the INLA implementation of the Bayesian GLM, as well as a classical GLM on simulated data. A validation of the method on data from the Human Connectome Project is also provided.


title: LISA: Learning Interpretable Skill Abstractions from Language

abstract: Learning policies that effectively utilize language instructions in complex, multi-task environments is an important problem in sequential decision-making. While it is possible to condition on the entire language instruction directly, such an approach could suffer from generalization issues. In our work, we propose \emph{Learning Interpretable Skill Abstractions (LISA)}, a hierarchical imitation learning framework that can learn diverse, interpretable primitive behaviors or skills from language-conditioned demonstrations to better generalize to unseen instructions. LISA uses vector quantization to learn discrete skill codes that are highly correlated with language instructions and the behavior of the learned policy. In navigation and robotic manipulation environments, LISA outperforms a strong non-hierarchical Decision Transformer baseline in the low data regime and is able to compose learned skills to solve tasks containing unseen long-range instructions. Our method demonstrates a more natural way to condition on language in sequential decision-making problems and achieve interpretable and controllable behavior with the learned skills.


title: Risk-averse controller design against data injection attacks on actuators for uncertain control systems

abstract: In this paper, we consider the optimal controller design problem against data injection attacks on actuators for an uncertain control system. We consider attacks that aim at maximizing the attack impact while remaining stealthy in the finite horizon. To this end, we use the Conditional Value-at-Risk to characterize the risk associated with the impact of attacks. The worst-case attack impact is characterized using the recently proposed output-to-output $\ell_2$-gain (OOG). We formulate the design problem and observe that it is non-convex and hard to solve. Using the framework of scenario-based optimization and a convex proxy for the OOG, we propose a convex optimization problem that approximately solves the design problem with probabilistic certificates. Finally, we illustrate the results through a numerical example.


title: An Empirical Study on Explanations in Out-of-Domain Settings

abstract: Recent work in Natural Language Processing has focused on developing approaches that extract faithful explanations, either via identifying the most important tokens in the input (i.e. post-hoc explanations) or by designing inherently faithful models that first select the most important tokens and then use them to predict the correct label (i.e. select-then-predict models). Currently, these approaches are largely evaluated on in-domain settings. Yet, little is known about how post-hoc explanations and inherently faithful models perform in out-of-domain settings. In this paper, we conduct an extensive empirical study that examines: (1) the out-of-domain faithfulness of post-hoc explanations, generated by five feature attribution methods; and (2) the out-of-domain performance of two inherently faithful models over six datasets. Contrary to our expectations, results show that in many cases out-of-domain post-hoc explanation faithfulness measured by sufficiency and comprehensiveness is higher compared to in-domain. We find this misleading and suggest using a random baseline as a yardstick for evaluating post-hoc explanation faithfulness. Our findings also show that select-then predict models demonstrate comparable predictive performance in out-of-domain settings to full-text trained models.


title: Differentiable Matrix Elements with MadJax

abstract: MadJax is a tool for generating and evaluating differentiable matrix elements of high energy scattering processes. As such, it is a step towards a differentiable programming paradigm in high energy physics that facilitates the incorporation of high energy physics domain knowledge, encoded in simulation software, into gradient based learning and optimization pipelines. MadJax comprises two components: (a) a plugin to the general purpose matrix element generator MadGraph that integrates matrix element and phase space sampling code with the JAX differentiable programming framework, and (b) a standalone wrapping API for accessing the matrix element code and its gradients, which are computed with automatic differentiation. The MadJax implementation and example applications of simulation based inference and normalizing flow based matrix element modeling, with capabilities enabled uniquely with differentiable matrix elements, are presented.

categories: hep-ph cs.LG physics.comp-ph physics.data-an


title: The r-Camassa-Holm equation: smooth and singular solutions

abstract: This paper introduces the r-Camassa-Holm (r-CH) equation, which describes a geodesic flow on the manifold of diffeomorphisms acting on the real line induced by the W1,r metric. The conserved energy is for the problem is given by the full W1,r norm and the for r = 2, we recover the Camassa-Holm equation. We compute the Lie symmetries for r-CH and study various symmetry reductions. We introduce singular weak solutions of the r-CH equation for r >= 2 and demonstrates their robustness in numerical simulations of their nonlinear interactions in both overtaking and head-on collisions. Several open questions are formulated about the unexplored properties of the r-CH weak singular solutions, including the question of whether they would emerge from smooth initial conditions.


title: Large-charge conformal dimensions at the $O(N)$ Wilson-Fisher fixed point

abstract: Recent work using a large-charge expansion for the $O(N)$ Wilson-Fisher conformal field theory has shown that the anomalous dimensions of large-charge operators can be expressed in terms of a few low-energy constants (LECs) of a large-charge effective field theory (EFT). By performing lattice Monte Carlo computations at the $O(N)$ Wilson-Fisher fixed point, we compute the anomalous dimensions of large-charge operators up to $N=8$ and charge $Q=10$, and extract the leading and subleading LECs of the $O(N)$ large-charge EFT. To alleviate the signal-to-noise ratio problem present in the large-charge sector of conventional lattice formulations of the $O(N)$ theory, we employ a recently developed qubit formulation of the $O(N)$ nonlinear sigma models with a worm algorithm. This enables us to test the validity of the large-charge expansion and the recent large-$N$ predictions for the coefficients of the large-charge EFT.


title: Strongly anisotropic electronic and magnetic structures in oxide dichlorides RuOCl$_2$ and OsOCl$_2$

abstract: Here, using density functional theory and density matrix renormalization group methods, we investigate the electronic and magnetic properties of RuOCl$_2$ and OsOCl$_2$ with $d^4$ electronic configurations. Different from a previous study using VOI$_2$ with $d^1$ configuration, these systems with $4d^4$ or $5d^4$ do not exhibit a ferroelectric instability along the $a$-axis. Due to the fully-occupied $d_{xy}$ orbital in RuOCl$_2$ and OsOCl$_2$, the Peierls instability distortion disappears along the $b$-axis, leading to an undistorted I${\rm mmm}$ phase (No. 71). Furthermore, we observe strongly anisotropic electronic and magnetic structures along the $a$-axis. The large crystal-field splitting energy (between $d_{xz/yz}$ and $d_{xy}$ orbitals) and large hopping between nearest-neighbor Ru and Os atoms suppresses the spin-orbital effect in $M$OCl$_2$ ($M$ = Ru or Os) with electronic density $n = 4$, resulting in a spin-1 system instead of a $J = 0$ singlet ground state. Moreover, we find staggered antiferromagnetic order with $\pi$ wavevector along the $M$-O chain direction ($a$-axis) while the magnetic coupling along the $b$-axis is weak. Based on Wannier functions from first-principles calculations, we calculated the relevant hopping amplitudes and crystal-field splitting energies of the $t_{2g}$ orbitals for the Os atoms to construct a multi-orbital Hubbard model for the $M$-O chains. Staggered AFM with $\uparrow$-$\downarrow$-$\uparrow$-$\downarrow$ spin structure dominates in our DMRG calculations, in agreement with DFT calculations.


title: Coefficients of Catalan States of Lattice Crossing I: $\Theta_{A}$-state Expansion

abstract: Plucking polynomial for plane rooted trees was introduced by J.H. Przytycki in 2014. As it was shown later, this polynomial can be used to find coefficients $C(A)$ of Catalan states $C$ of $m \times n$-lattice crossing $L(m,n)$ without returns on one side. In this paper, we show that $C(A)$ for any $C$ can be found by using $\Theta_{A}$-state expansion which represents $C(A)$ as a linear combination of coefficients of Catalan states with no top returns over $\mathbb{Q}(A)$. We also provide an algorithm for finding $\Theta_{A}$-state expansions and examples of its applications. Finally, an example of a Catalan state with non-unimodal coefficient is given.


title: Robust Causal Inference of Drug-drug Interactions

abstract: There is growing interest in developing causal inference methods for multi-valued treatments with a focus on pairwise average treatment effects. Here we focus on a clinically important, yet less-studied estimand: causal drug-drug interactions (DDIs), which quantifies the degree to which the causal effect of drug A is altered by the presence versus the absence of drug B. Confounding adjustment when studying the effects of DDIs can be accomplished via inverse probability of treatment weighting (IPTW), a standard approach originally developed for binary treatments and later generalized to multi-valued treatments. However, this approach generally results in biased results when the propensity score model is misspecified. Motivated by the need for more robust techniques, we propose two empirical likelihood-based weighting approaches that allow for specifying a set of propensity score models, with the second method balancing user-specified covariates directly, by incorporating additional, nonparametric constraints. The resulting estimators from both methods are consistent when the postulated set of propensity score models contains a correct one; this property has been termed multiple robustness. We then evaluate their finite sample performance through simulation. The results demonstrate that the proposed estimators outperform the standard IPTW method in terms of both robustness and efficiency. Finally, we apply the proposed methods to evaluate the impact of renin-angiotensin system inhibitors (RAS-I) on the comparative nephrotoxicity of nonsteroidal anti-inflammatory drugs (NSAID) and opioids, using data derived from electronic medical records from a large multi-hospital health system.


title: Structure from Voltage

abstract: Effective resistance (ER) is an attractive way to interrogate the structure of graphs. It is an alternative to computing the eigen-vectors of the graph Laplacian. Graph laplacians are used to find low dimensional structures in high dimensional data. Here too, ER based analysis has advantages over eign-vector based methods. Unfortunately Von Luxburg et al. (2010) show that, when vertices correspond to a sample from a distribution over a metric space, the limit of the ER between distant points converges to a trivial quantity that holds no information about the structure of the graph. We show that by using scaling resistances in a graph with $n$ vertices by $n^2$, one gets a meaningful limit of the voltages and of effective resistances. We also show that by adding a "ground" node to a metric graph one gets a simple and natural way to compute all of the distances from a chosen point to all other points.


title: Piezoelectric Strain FET (PeFET) based Non-Volatile Memories

abstract: We propose non-volatile memory (NVM) designs based on Piezoelectric Strain FET (PeFET) utilizing a piezoelectric/ferroelectric (PE/FE such as PZT) coupled with 2D Transition Metal Dichalcogenide (2D-TMD such as MoS2) transistor. The proposed NVMs store bit information in the form of polarization (P) of the FE/PE, use electric-field driven P-switching for write and employ piezoelectricity induced dynamic bandgap modulation of 2D-TMD channel for bit sensing. We analyze PeFET with COMSOL based 3D modeling showing that the circuit-driven optimization of PeFET geometry is essential to achieve effective hammer-and-nail effect and adequate bandgap modulation for NVM read. Our results show that distinguishability of binary states to up to 11X is achieved in PeFETs.We propose various flavors of PeFET NVMs, namely (a) high density (HD) NVM featuring a compact access-transistor-less bit-cell, (b) 1T-1PeFET NVM with segmented architecture, targeted for optimized write energy and latency and (c) cross-coupled (CC) NVM offering a trade-off between area and latency.PeFET NVMs offer up to 7X smaller cell area, 66% lower write energy, 87% lower read energy and 44% faster read compared to 2D-FET SRAM. This comes at the cost of high write latency in PeFET NVMs, which can be minimized by virtue of optimized PE geometry.


title: Self-Repelling Elastic Manifolds with Low Dimensional Range

abstract: We consider self-repelling elastic manifolds with a domain $[-N,N]^d \cap \mathbb{Z}^d$, that take values in $\mathbb{R}^D$. Our main result states that when the dimension of the domain is $d=2$ and the dimension of the range is $D=1$, the effective radius $R_N$ of the manifold is approximately $N^{4/3}$. This verifies the conjecture of Kantor, Kardar and Nelson [7]. Our results for the case where $d \geq 3$ and $D <d$ give a lower bound on $R_N$ of order $N^{\frac{1}{D} \left(d-\frac{2(d-D)}{D+2} \right)}$ and an upper bound proportional to $N^{\frac{d}{2}+\frac{d-D}{D+2}}$. These results imply that self-repelling elastic manifolds with a low dimensional range undergo a significantly stronger stretching than in the case where $d=D$, which was studied by the authors in [10].


title: Magnetic Interactions of the Centrosymmetric Skyrmion Material Gd2PdSi3

abstract: The experimental realization of magnetic skyrmions in centrosymmetric materials has been driven by theoretical understanding of how a delicate balance of anisotropy and frustration can stabilize topological spin structures in applied magnetic fields. Recently, the centrosymmetric material Gd$_{2}$PdSi$_{3}$ was shown to host a field-induced skyrmion phase, but the skyrmion stabilization mechanism remains unclear. Here, we employ neutron-scattering measurements on an isotopically-enriched polycrystalline Gd$_{2}$PdSi$_{3}$ sample to quantify the interactions that drive skyrmion formation. Our analysis reveals spatially-extended interactions in triangular planes that are consistent with an RKKY mechanism, and large ferromagnetic inter-planar magnetic interactions that are modulated by the Pd/Si superstructure. The skyrmion phase emerges from a zero-field helical magnetic order with magnetic moments perpendicular to the magnetic propagation vector, indicating that the magnetic dipolar interaction plays a significant role. Our experimental results establish an interaction space that can promote skyrmion formation, facilitating identification and design of centrosymmetric skyrmion materials.


title: Spin dynamics in the square-lattice cupola system Ba(TiO)Cu$_4$(PO$_4$)$_4$

abstract: We report high-resolution single-crystal inelastic neutron scattering measurements on the spin-1/2 antiferromagnet Ba(TiO)Cu$_4$(PO$_4$)$_4$. This material is formed from layers of four-site \cupola" structures, oriented alternately upwards and downwards, which constitute a rather special realization of two-dimensional (2D) square-lattice magnetism. The strong Dzyaloshinskii-Moriya (DM) interaction within each cupola, or plaquette, unit has a geometry largely unexplored among the numerous studies of magnetic properties in 2D Heisenberg models with spin and spatial anisotropies. We have measured the magnetic excitations at zero field and in fields up to 5 T, finding a complex mode structure with multiple characteristic features that allow us to extract all the relevant magnetic interactions by modelling within the linear spin-wave approximation. We demonstrate that Ba(TiO)Cu$_4$(PO$_4$)$_4$ is a checkerboard system with almost equal intra- and inter-plaquette couplings, in which the intra-plaquette DM interaction is instrumental both in enforcing robust magnetic order and in opening a large gap at the Brillouin-zone center. We place our observations in the perspective of generalized phase diagrams for spin-1/2 square-lattice models and materials, where exploring anisotropies and frustration as routes to quantum disorder remains a frontier research problem.


title: Perturbation of invariant subspaces for ill-conditioned eigensystem

abstract: Given a diagonalizable matrix $A$, we study the stability of its invariant subspaces when its matrix of eigenvectors is ill-conditioned. Let $\mathcal{X}_1$ be some invariant subspace of $A$ and $X_1$ be the matrix storing the right eigenvectors that spanned $\mathcal{X}_1$. It is generally believed that when the condition number $\kappa_2(X_1)$ gets large, the corresponding invariant subspace $\mathcal{X}_1$ will become unstable to perturbation. This paper proves that this is not always the case. Specifically, we show that the growth of $\kappa_2(X_1)$ alone is not enough to destroy the stability. As a direct application, our result ensures that when $A$ gets closer to a Jordan form, one may still estimate its invariant subspaces from the noisy data stably.


title: Optimal Transport-based Graph Matching for 3D retinal OCT image registration

abstract: Registration of longitudinal optical coherence tomography (OCT) images assists disease monitoring and is essential in image fusion applications. Mouse retinal OCT images are often collected for longitudinal study of eye disease models such as uveitis, but their quality is often poor compared with human imaging. This paper presents a novel but efficient framework involving an optimal transport based graph matching (OT-GM) method for 3D mouse OCT image registration. We first perform registration of fundus-like images obtained by projecting all b-scans of a volume on a plane orthogonal to them, hereafter referred to as the x-y plane. We introduce Adaptive Weighted Vessel Graph Descriptors (AWVGD) and 3D Cube Descriptors (CD) to identify the correspondence between nodes of graphs extracted from segmented vessels within the OCT projection images. The AWVGD comprises scaling, translation and rotation, which are computationally efficient, whereas CD exploits 3D spatial and frequency domain information. The OT-GM method subsequently performs the correct alignment in the x-y plane. Finally, registration along the direction orthogonal to the x-y plane (the z-direction) is guided by the segmentation of two important anatomical features peculiar to mouse b-scans, the Internal Limiting Membrane (ILM) and the hyaloid remnant (HR). Both subjective and objective evaluation results demonstrate that our framework outperforms other well-established methods on mouse OCT images within a reasonable execution time.


title: Decisions over Sequences

abstract: This paper introduces a class of objects called decision rules that map infinite sequences of alternatives to a decision space. These objects can be used to model situations where a decision maker encounters alternatives in a sequence such as receiving recommendations. Within the class of decision rules, we study natural subclasses: stopping and uniform stopping rules. Our main result establishes the equivalence of these two subclasses of decision rules. Next, we introduce the notion of computability of decision rules using Turing machines and show that computable rules can be implemented using a simpler computational device: a finite automaton. We further show that computability of choice rules -- an important subclass of decision rules -- is implied by their continuity with respect to a natural topology. Finally, we introduce some natural heuristics in this framework and provide their behavioral characterization.


title: A graph related to the sum of element orders of a finite group

abstract: A finite group is called $\psi$-divisible iff $\psi(H)|\psi(G)$ for any subgroup $H$ of a finite group $G$. Here, $\psi(G)$ is the sum of element orders of $G$. For now, the only known examples of such groups are the cyclic ones of square-free order. The existence of non-abelian $\psi$-divisible groups still constitutes an open question. The aim of this paper is to make a connection between the $\psi$-divisibility property and graph theory. Hence, for a finite group $G$, we introduce a simple undirected graph called the $\psi$-divisibility graph of $G$. We denote it by $\psi_G$. Its vertices are the non-trivial subgroups of $G$, while two distinct vertices $H$ and $K$ are adjacent iff $H\subset K$ and $\psi(H)|\psi(K)$ or $K\subset H$ and $\psi(K)|\psi(H)$. We prove that $G$ is $\psi$-divisible iff $\psi_G$ has a universal (dominating) vertex. Also, we study various properties of $\psi_G$, when $G$ is a finite cyclic group. The choice of restricting our study to this specific class of groups is motivated in the paper.


title: Parametrized and equivariant higher algebra

abstract: We develop the rudiments of a theory of parametrized $\infty$-operads, including parametrized generalizations of monoidal envelopes, Day convolution, operadic left Kan extensions, results on limits and colimits of algebras, and the symmetric monoidal Yoneda embedding.


title: Structure Extraction in Task-Oriented Dialogues with Slot Clustering

abstract: Extracting structure information from dialogue data can help us better understand user and system behaviors. In task-oriented dialogues, dialogue structure has often been considered as transition graphs among dialogue states. However, annotating dialogue states manually is expensive and time-consuming. In this paper, we propose a simple yet effective approach for structure extraction in task-oriented dialogues. We first detect and cluster possible slot tokens with a pre-trained model to approximate dialogue ontology for a target domain. Then we track the status of each identified token group and derive a state transition structure. Empirical results show that our approach outperforms unsupervised baseline models by far in dialogue structure extraction. In addition, we show that data augmentation based on extracted structures enriches the surface formats of training data and can achieve a significant performance boost in dialogue response generation.


title: How Well Can We Measure Galaxy Dust Attenuation Curves? The Impact of the Assumed Star-Dust Geometry Model in SED Fitting

abstract: One of the most common methods for inferring galaxy attenuation curves is via spectral energy distribution (SED) modeling, where the dust attenuation properties are modeled simultaneously with other galaxy physical properties. In this paper, we assess the ability of SED modeling to infer these dust attenuation curves from broadband photometry, and suggest a new flexible model that greatly improves the accuracy of attenuation curve derivations. To do this, we fit mock SEDs generated from the Simba cosmological simulation with the Prospector SED fitting code. We consider the impact of the commonly-assumed uniform screen model and introduce a new non-uniform screen model parameterized by the fraction of unobscured stellar light. This non-uniform screen model allows for a non-zero fraction of stellar light to remain unattenuated, resulting in a more flexible attenuation curve shape by decoupling the shape of the UV attenuation curve from the optical attenuation curve. The ability to constrain the dust attenuation curve is significantly improved with the use of a non-uniform screen model, with the median offset in UV attenuation decreasing from $-0.30$ dex with a uniform screen model to $-0.17$ dex with the non-uniform screen model. With this increase in dust attenuation modeling accuracy, we also improve the star formation rates (SFRs) inferred with the non-uniform screen model, decreasing the SFR offset on average by $0.12$ dex. We discuss the efficacy of this new model, focusing on caveats with modeling star-dust geometries and the constraining power of available SED observations.


title: Long-time asymptotics and regularity estimates for weak solutions to a doubly degenerate thin-film equation in the Taylor-Couette setting

abstract: We study the dynamic behaviour of solutions to a fourth-order quasilinear degenerate parabolic equation for large times arising in fluid dynamical applications. The degeneracy occurs both with respect to the unknown and with respect to the sum of its first and third spatial derivative. The modelling equation appears as a thin-film limit for the interface separating two immiscible viscous fluid films confined between two cylinders rotating at small relative angular velocity. More precisely, the fluid occupying the layer next to the outer cylinder is considered to be Newtonian, i.e. it has constant viscosity, while we assume that the layer next to the inner cylinder is filled by a shear-thinning power-law fluid. Using energy methods, Fourier analysis and suitable regularity estimates for higher-order parabolic equations, we prove global existence of positive weak solutions in the case of low initial energy. Moreover, these global solutions are polynomially stable, in the sense that interfaces which are initially close to a circle, converge at rate $1/t^{1/\beta}$ for some $\beta > 0$ to a circle, as time tends to infinity. In addition, we provide regularity estimates for general nonlinear degenerate parabolic equations of fourth order.


title: Robust Multi-Agent Bandits Over Undirected Graphs

abstract: We consider a multi-agent multi-armed bandit setting in which $n$ honest agents collaborate over a network to minimize regret but $m$ malicious agents can disrupt learning arbitrarily. Assuming the network is the complete graph, existing algorithms incur $O( (m + K/n) \log (T) / \Delta )$ regret in this setting, where $K$ is the number of arms and $\Delta$ is the arm gap. For $m \ll K$, this improves over the single-agent baseline regret of $O(K\log(T)/\Delta)$. In this work, we show the situation is murkier beyond the case of a complete graph. In particular, we prove that if the state-of-the-art algorithm is used on the undirected line graph, honest agents can suffer (nearly) linear regret until time is doubly exponential in $K$ and $n$. In light of this negative result, we propose a new algorithm for which the $i$-th agent has regret $O( ( d_{\text{mal}}(i) + K/n) \log(T)/\Delta)$ on any connected and undirected graph, where $d_{\text{mal}}(i)$ is the number of $i$'s neighbors who are malicious. Thus, we generalize existing regret bounds beyond the complete graph (where $d_{\text{mal}}(i) = m$), and show the effect of malicious agents is entirely local (in the sense that only the $d_{\text{mal}}(i)$ malicious agents directly connected to $i$ affect its long-term regret).


title: One Model is All You Need: Multi-Task Learning Enables Simultaneous Histology Image Segmentation and Classification

abstract: The recent surge in performance for image analysis of digitised pathology slides can largely be attributed to the advances in deep learning. Deep models can be used to initially localise various structures in the tissue and hence facilitate the extraction of interpretable features for biomarker discovery. However, these models are typically trained for a single task and therefore scale poorly as we wish to adapt the model for an increasing number of different tasks. Also, supervised deep learning models are very data hungry and therefore rely on large amounts of training data to perform well. In this paper, we present a multi-task learning approach for segmentation and classification of nuclei, glands, lumina and different tissue regions that leverages data from multiple independent data sources. While ensuring that our tasks are aligned by the same tissue type and resolution, we enable meaningful simultaneous prediction with a single network. As a result of feature sharing, we also show that the learned representation can be used to improve the performance of additional tasks via transfer learning, including nuclear classification and signet ring cell detection. As part of this work, we train our developed Cerberus model on a huge amount of data, consisting of over 600K objects for segmentation and 440K patches for classification. We use our approach to process 599 colorectal whole-slide images from TCGA, where we localise 377 million, 900K and 2.1 million nuclei, glands and lumina, respectively and make the results available to the community for downstream analysis.


title: Elliptical Slice Sampling for Probabilistic Verification of Stochastic Systems with Signal Temporal Logic Specifications

abstract: Autonomous robots typically incorporate complex sensors in their decision-making and control loops. These sensors, such as cameras and Lidars, have imperfections in their sensing and are influenced by environmental conditions. In this paper, we present a method for probabilistic verification of linearizable systems with Gaussian and Gaussian mixture noise models (e.g. from perception modules, machine learning components). We compute the probabilities of task satisfaction under Signal Temporal Logic (STL) specifications, using its robustness semantics, with a Markov Chain Monte-Carlo slice sampler. As opposed to other techniques, our method avoids over-approximations and double-counting of failure events. Central to our approach is a method for efficient and rejection-free sampling of signals from a Gaussian distribution such that satisfy or violate a given STL formula. We show illustrative examples from applications in robot motion planning.


title: Direct observation of discommensurate charge density wave modulation in the quasi-1D Weyl semimetal candidate NbTe$_4$

abstract: The transition-metal tetrachalcogenides are a model system to explore the conjunction of correlated electronic states such as charge density waves (CDW), with topological phases of matter. Understanding the connection between these phases requires a thorough understanding of the individual states, which for the case of the CDW in this system, is still missing. In this paper we combine phonon-structure calculations and scanning tunneling microscopy measurements of NbTe$_4$ in order to provide a full characterization of the CDW state. We find that, at short range, the superstructure formed by the CDW is fully commensurate with the lattice parameters. Moreover, our data reveals the presence of phase-slip domain-walls separating regions of commensurate-CDW in the nanoscale, indicating that the CDW in this compound is discommensurate at long-range. Our results solve a long-standing discussion about the nature of the CDW in these materials, and provide a strong basis for the study of the interplay between this state and other novel quantum electronic states.

categories: cond-mat.mtrl-sci cond-mat.mes-hall cond-mat.other cond-mat.str-el


title: Deep Camera Pose Regression Using Pseudo-LiDAR

abstract: An accurate and robust large-scale localization system is an integral component for active areas of research such as autonomous vehicles and augmented reality. To this end, many learning algorithms have been proposed that predict 6DOF camera pose from RGB or RGB-D images. However, previous methods that incorporate depth typically treat the data the same way as RGB images, often adding depth maps as additional channels to RGB images and passing them through convolutional neural networks (CNNs). In this paper, we show that converting depth maps into pseudo-LiDAR signals, previously shown to be useful for 3D object detection, is a better representation for camera localization tasks by projecting point clouds that can accurately determine 6DOF camera pose. This is demonstrated by first comparing localization accuracies of a network operating exclusively on pseudo-LiDAR representations, with networks operating exclusively on depth maps. We then propose FusionLoc, a novel architecture that uses pseudo-LiDAR to regress a 6DOF camera pose. FusionLoc is a dual stream neural network, which aims to remedy common issues with typical 2D CNNs operating on RGB-D images. The results from this architecture are compared against various other state-of-the-art deep pose regression implementations using the 7 Scenes dataset. The findings are that FusionLoc performs better than a number of other camera localization methods, with a notable improvement being, on average, 0.33m and 4.35{\deg} more accurate than RGB-D PoseNet. By proving the validity of using pseudo-LiDAR signals over depth maps for localization, there are new considerations when implementing large-scale localization systems.


title: Asymptotic Normality of Gini Correlation in High Dimension with Applications to the K-sample Problem

abstract: The categorical Gini correlation proposed by Dang et al. is a dependence measure to characterize independence between categorical and numerical variables. The asymptotic distributions of the sample correlation under dependence and independence have been established when the dimension of the numerical variable is fixed. However, its asymptotic behavior for high dimensional data has not been explored. In this paper, we develop the central limit theorem for the Gini correlation in the more realistic setting where the dimensionality of the numerical variable is diverging. We then construct a powerful and consistent test for the $K$-sample problem based on the asymptotic normality. The proposed test not only avoids computation burden but also gains power over the permutation procedure. Simulation studies and real data illustrations show that the proposed test is more competitive to existing methods across a broad range of realistic situations, especially in unbalanced cases.


title: High-precision real-space simulation of electrostatically-confined few-electron states

abstract: In this paper we present a computational procedure that utilizes real-space grids to obtain high precision approximations of electrostatically confined few-electron states such as those that arise in gated semiconductor quantum dots. We use the Full Configuration Interaction (FCI) method with a continuously adapted orthonormal orbital basis to approximate the ground and excited states of such systems. We also introduce a benchmark problem based on a realistic analytical electrostatic potential for quantum dot devices. We show that our approach leads to highly precise computed energies and energy differences over a wide range of model parameters. The analytic definition of the benchmark allows for a collection of tests that are easily replicated, thus facilitating comparisons with other computational approaches.


title: Sampling-Based Winner Prediction in District-Based Elections

abstract: In a district-based election, we apply a voting rule $r$ to decide the winners in each district, and a candidate who wins in a maximum number of districts is the winner of the election. We present efficient sampling-based algorithms to predict the winner of such district-based election systems in this paper. When $r$ is plurality and the margin of victory is known to be at least $\varepsilon$ fraction of the total population, we present an algorithm to predict the winner. The sample complexity of our algorithm is $\mathcal{O}\left(\frac{1}{\varepsilon^4}\log \frac{1}{\varepsilon}\log\frac{1}{\delta}\right)$. We complement this result by proving that any algorithm, from a natural class of algorithms, for predicting the winner in a district-based election when $r$ is plurality, must sample at least $\Omega\left(\frac{1}{\varepsilon^4}\log\frac{1}{\delta}\right)$ votes. We then extend this result to any voting rule $r$. Loosely speaking, we show that we can predict the winner of a district-based election with an extra overhead of $\mathcal{O}\left(\frac{1}{\varepsilon^2}\log\frac{1}{\delta}\right)$ over the sample complexity of predicting the single-district winner under $r$. We further extend our algorithm for the case when the margin of victory is unknown, but we have only two candidates. We then consider the median voting rule when the set of preferences in each district is single-peaked. We show that the winner of a district-based election can be predicted with $\mathcal{O}\left(\frac{1}{\varepsilon^4}\log\frac{1}{\varepsilon}\log\frac{1}{\delta}\right)$ samples even when the harmonious order in different districts can be different and even unknown. Finally, we also show some results for estimating the margin of victory of a district-based election within both additive and multiplicative error bounds.


title: Competitors-Aware Stochastic Lap Strategy Optimisation for Race Hybrid Vehicles

abstract: World Endurance Championship (WEC) racing events are characterised by a relevant performance gap among competitors. The fastest vehicles category, consisting in hybrid vehicles, has to respect energy usage constraints set by the technical regulation. Considering absence of competitors, i.e. traffic conditions, the optimal energy usage strategy for lap time minimisation is typically computed through a constrained optimisation problem. To the best of our knowledge, the majority of state-of-the-art works neglects competitors. This leads to a mismatch with the real world, where traffic generates considerable time losses. To bridge this gap, we propose a new framework to offline compute optimal strategies for the powertrain energy management considering competitors. Through analysis of the available data from previous events, statistics on the sector times and overtaking probabilities are extracted to encode the competitors' behaviour. Adopting a multi-agent model, the statistics are then used to generate realistic Monte Carlo (MC) simulation of their position along the track. The simulator is then adopted to identify the optimal strategy as follows. We develop a longitudinal vehicle model for the ego-vehicle and implement an optimisation problem for lap time minimisation in absence of traffic, based on Genetic Algorithms. Solving the optimisation problem for a variety of constraints generates a set of candidate optimal strategies. Stochastic Dynamic Programming is finally implemented to choose the best strategy considering competitors, whose motion is generated by the MC simulator. Our approach, validated on data from a real stint of race, allows to significantly reduce the lap time.


title: On energy for accelerating observers in black hole spacetimes

abstract: Quasi-local energies for constantly accelerating observers in Ba\~nados, Teitelboim, and Zanelli (BTZ), Schwarzschild and Schwarzschild-de Sitter spacetimes are derived. The energies are expressed in terms of acceleration, cosmological constant, and area, quantities measurable by the observers. Based on results from quantum fields in curved spacetime for the redshifted Hawking temperature, entropy and thermodynamic-like laws are briefly explored in the three spacetimes.


title: Pippi: Practical Protocol Instantiation

abstract: A protocol specifies interactions between roles, which together constitute a multiagent system (MAS). Enacting a protocol presupposes that agents are bound to the its roles. Existing protocol-based approaches, however, do not adequately treat the practical aspects of how roles bindings come about. Pippi addresses this problem of MAS instantiation. It proposes the notion of a metaprotocol, enacting which instantiates a MAS suitable for enacting a given protocol. Pippi demonstrates the subtleties involved in instantiating MAS arising from protocol composition, correlation, and decentralization. To address these subtleties and further support practical application patterns, we introduce an enhanced protocol language, with support for parameter types (including role and protocol typed parameters, for metaprotocols), interface flexibility, and binding constraints. We discuss the realization of our approach through an extended agent architecture, including the novel concept of a MAS adapter for contact management. We evaluate Pippi's expressiveness by demonstrating common patterns for agent discovery.


title: Using Multivariate Imputation by Chained Equations to Predict Redshifts of Active Galactic Nuclei

abstract: Redshift measurement of active galactic nuclei (AGNs) remains a time-consuming and challenging task, as it requires follow up spectroscopic observations and detailed analysis. Hence, there exists an urgent requirement for alternative redshift estimation techniques. The use of machine learning (ML) for this purpose has been growing over the last few years, primarily due to the availability of large-scale galactic surveys. However, due to observational errors, a significant fraction of these data sets often have missing entries, rendering that fraction unusable for ML regression applications. In this study, we demonstrate the performance of an imputation technique called Multivariate Imputation by Chained Equations (MICE), which rectifies the issue of missing data entries by imputing them using the available information in the catalog. We use the Fermi-LAT Fourth Data Release Catalog (4LAC) and impute 24% of the catalog. Subsequently, we follow the methodology described in Dainotti et al. (2021) and create an ML model for estimating the redshift of 4LAC AGNs. We present results which highlight positive impact of MICE imputation technique on the machine learning models performance and obtained redshift estimation accuracy.


title: Virtual Reference Feedback Tuning for linear discrete-time systems with robust stability guarantees based on Set Membership

abstract: In this paper we propose a novel methodology that allows to design, in a purely data-based fashion and for linear single-input and single-output systems, both robustly stable and performing control systems for tracking piecewise constant reference signals. The approach uses both (i) Virtual Reference Feedback Tuning for enforcing suitable performances and (ii) the Set Membership framework for providing a-priori robust stability guarantees. Indeed, an uncertainty set for the system parameters is obtained through Set Membership identification, where an algorithm based on the scenario approach is proposed to estimate the inflation parameter in a probabilistic way. Based on this set, robust stability conditions are enforced as Linear Matrix Inequality constraints within an optimization problem whose linear cost function relies on Virtual Reference Feedback Tuning. To show the generality and effectiveness of our approach, we apply it to two of the most widely used yet simple control schemes, i.e., where tracking is achieved thanks to (i) a static feedforward action and (ii) an integrator in closed-loop. The proposed method is not fully direct due to the Set Membership identification. However, the uncertainty set is used with the only objective of providing robust stability guarantees for the closed-loop system and it is not directly used for the performances optimization, which instead is totally based on data. The effectiveness of the developed method is demonstrated with reference to two simulation examples. A comparison with other data-driven methods is also carried out.


title: Amortized Proximal Optimization

abstract: We propose a framework for online meta-optimization of parameters that govern optimization, called Amortized Proximal Optimization (APO). We first interpret various existing neural network optimizers as approximate stochastic proximal point methods which trade off the current-batch loss with proximity terms in both function space and weight space. The idea behind APO is to amortize the minimization of the proximal point objective by meta-learning the parameters of an update rule. We show how APO can be used to adapt a learning rate or a structured preconditioning matrix. Under appropriate assumptions, APO can recover existing optimizers such as natural gradient descent and KFAC. It enjoys low computational overhead and avoids expensive and numerically sensitive operations required by some second-order optimizers, such as matrix inverses. We empirically test APO for online adaptation of learning rates and structured preconditioning matrices for regression, image reconstruction, image classification, and natural language translation tasks. Empirically, the learning rate schedules found by APO generally outperform optimal fixed learning rates and are competitive with manually tuned decay schedules. Using APO to adapt a structured preconditioning matrix generally results in optimization performance competitive with second-order methods. Moreover, the absence of matrix inversion provides numerical stability, making it effective for low precision training.


title: Assigned rational functions of a rooted tree

abstract: We investigate the spectral properties of rooted trees with the intention of improving the currently existing results that deal with this matter. The concept of an assigned rational function is recursively defined for each vertex of a rooted tree. Afterwards, two mathematical formulas are given which show how the characteristic polynomials of the adjacency and Laplacian matrix can be represented as products of the aforementioned rational functions. In order to demonstrate their general use case scenario, the obtained formulas are subsequently implemented on balanced trees, with a special focus on the Bethe trees. In the end, some of the previously derived results are used in order to construct a tree merging procedure which preserves the spectra of all of the starting trees.


title: Dynamic N:M Fine-grained Structured Sparse Attention Mechanism

abstract: Transformers are becoming the mainstream solutions for various tasks like NLP and Computer vision. Despite their success, the high complexity of the attention mechanism hinders them from being applied to latency-sensitive tasks. Tremendous efforts have been made to alleviate this problem, and many of them successfully reduce the asymptotic complexity to linear. Nevertheless, most of them fail to achieve practical speedup over the original full attention under moderate sequence lengths and are unfriendly to finetuning. In this paper, we present DFSS, an attention mechanism that dynamically prunes the full attention weight matrix to N:M fine-grained structured sparse pattern. We provide both theoretical and empirical evidence that demonstrates DFSS is a good approximation of the full attention mechanism. We propose a dedicated CUDA kernel design that completely eliminates the dynamic pruning overhead and achieves speedups under arbitrary sequence length. We evaluate the 1:2 and 2:4 sparsity under different configurations and achieve 1.27~ 1.89x speedups over the full-attention mechanism. It only takes a couple of finetuning epochs from the pretrained model to achieve on par accuracy with full attention mechanism on tasks from various domains under different sequence lengths from 384 to 4096.


title: Explicit Non-Abelian Gerbes with Connections

abstract: We define the notion of adjustment for strict Lie 2-groups and provide the complete cocycle description for non-Abelian gerbes with connections whose structure 2-group is an adjusted 2-group. Most importantly, we depart from the common fake-flat connections and employ adjusted connections. This is an important generalisation that is needed for physical applications especially in the context of supergravity. We give a number of explicit examples; in particular, we lift the spin structure on $S^4$, corresponding to an instanton-anti-instanton pair, to a string structure, a 2-group bundle with connection. We also outline how categorified forms of Bogomolny monopoles known as self-dual strings can be obtained via a Penrose-Ward transform of string bundles over twistor space.


title: Dynamic Control of Service Systems with Returns: Application to Design of Post-Discharge Hospital Readmission Prevention Programs

abstract: We study a control problem for queueing systems where customers may return for additional episodes of service after their initial service completion. At each service completion epoch, the decision maker can choose to reduce the probability of return for the departing customer but at a cost that is convex increasing in the amount of reduction in the return probability. Other costs are incurred as customers wait in the queue and every time they return for service. Our primary motivation comes from post-discharge Quality Improvement (QI) interventions (e.g., follow up phone-calls, appointments) frequently used in a variety of healthcare settings to reduce unplanned hospital readmissions. Our objective is to understand how the cost of interventions should be balanced with the reductions in congestion and service costs. To this end, we consider a fluid approximation of the queueing system and characterize the structure of optimal long-run average and bias-optimal transient control policies for the fluid model. Our structural results motivate the design of intuitive surge protocols whereby different intensities of interventions (corresponding to different levels of reduction in the return probability) are provided based on the congestion in the system. Through extensive simulation experiments, we study the performance of the fluid policy for the stochastic system and identify parameter regimes where it leads to significant cost savings compared to a fixed long-run average optimal policy that ignores holding costs and a simple policy that uses the highest level of intervention whenever the queue is non-empty. In particular, we find that in a parameter regime relevant to our motivating application, dynamically adjusting the intensity of interventions could result in up to 25.4% reduction in long-run average cost and 33.7% in finite-horizon costs compared to the simple aggressive policy.


title: On the decategorification of some higher actions in Heegaard Floer homology

abstract: We decategorify the higher actions on bordered Heegaard Floer strands algebras from recent work of Rouquier and the author and identify the decategorifications with certain actions on exterior powers of homology groups of surfaces. We also suggest an interpretation for these actions in the language of open-closed TQFT, and we prove a corresponding gluing formula.


title: Distributed randomized Kaczmarz for the adversarial workers

abstract: Developing large-scale distributed methods that are robust to the presence of adversarial or corrupted workers is an important part of making such methods practical for real-world problems. Here, we propose an iterative approach that is adversary-tolerant for least-squares problems. The algorithm utilizes simple statistics to guarantee convergence and is capable of learning the adversarial distributions. Additionally, the efficiency of the proposed method is shown in simulations in the presence of adversaries. The results demonstrate the great capability of such methods to tolerate different levels of adversary rates and to identify the erroneous workers with high accuracy.


title: On quantitative hypocoercivity estimates based on Harris-type theorems

abstract: This review concerns recent results on the quantitative study of convergence towards the stationary state for spatially inhomogeneous kinetic equations. We focus on analytical results obtained by means of certain probabilistic techniques from the ergodic theory of Markov processes. These techniques are sometimes referred to as Harris-type theorems. They provide constructive proofs for convergence results in the $L^1$ (or total variation) setting for a large class of initial data. The convergence rates can be made explicit (both for geometric and sub-geometric rates) by tracking the constants appearing in the hypotheses. Harris-type theorems are particularly well-adapted for equations exhibiting non-explicit and non-equilibrium steady states since they do not require prior information on the existence of stationary states. This allows for significant improvements of some already-existing results by relaxing assumptions and providing explicit convergence rates. We aim to present Harris-type theorems, providing a guideline on how to apply these techniques to the kinetic equations at hand. We discuss recent quantitative results obtained for kinetic equations in gas theory and mathematical biology, giving some perspectives on potential extensions to nonlinear equations.


title: Estimating causal effects with optimization-based methods: A review and empirical comparison

abstract: In the absence of randomized controlled and natural experiments, it is necessary to balance the distributions of (observable) covariates of the treated and control groups in order to obtain an unbiased estimate of a causal effect of interest; otherwise, a different effect size may be estimated, and incorrect recommendations may be given. To achieve this balance, there exist a wide variety of methods. In particular, several methods based on optimization models have been recently proposed in the causal inference literature. While these optimization-based methods empirically showed an improvement over a limited number of other causal inference methods in their relative ability to balance the distributions of covariates and to estimate causal effects, they have not been thoroughly compared to each other and to other noteworthy causal inference methods. In addition, we believe that there exist several unaddressed opportunities that operational researchers could contribute with their advanced knowledge of optimization, for the benefits of the applied researchers that use causal inference tools. In this review paper, we present an overview of the causal inference literature and describe in more detail the optimization-based causal inference methods, provide a comparative analysis of the prevailing optimization-based methods, and discuss opportunities for new methods.


title: Non-linear Smoothing for the Periodic Generalized Non-linear Schr\"odinger Equation

abstract: We consider the periodic non-linear Schr\"odinger equation with non-linearity given by $|u|^{p-1}u$ for odd $p > 1$ in dimension $1$. We first establish that the difference between the non-linear evolution and a phase rotation of the the linear evolution is in a smoother space. We then study forced and damped defocusing non-linear Schr\"odinger equations of the above type and establish an analogous smoothing statement that extends globally in time. As a corollary we establish both existence and smootheness for global attractors in the energy space.


title: A computational model of self-organized shape dynamics of active surfaces in fluids

abstract: Mechanochemical processes on surfaces such as the cellular cortex or epithelial sheets, play a key role in determining patterns and shape changes of biological systems. To understand the complex interplay of hydrodynamics and material flows on such active surfaces requires novel numerical tools. Here, we present a finite-element method for an active deformable surface interacting with the surrounding fluids. The underlying model couples surface and bulk hydrodynamics to surface flow of a diffusible species which generates active contractile forces. The method is validated with previous results based on linear stability analysis and shows almost perfect agreement regarding predicted patterning. Away from the linear regime we find rich non-linear behavior, such as the presence of multiple stationary states. We study the formation of a contractile ring on the surface and the corresponding shape changes. Finally, we explore mechanochemical pattern formation on various surface geometries and find that patterning strongly adapts to local surface curvature. The developed method provides a basis to analyze a variety of systems that involve mechanochemical pattern formation on active surfaces interacting with surrounding fluids.


title: Making sense of complex systems through resolution, relevance, and mapping entropy

abstract: Complex systems are characterised by a tight, nontrivial interplay of their constituents, which gives rise to a multi-scale spectrum of emergent properties. In this scenario, it is practically and conceptually difficult to identify those degrees of freedom that mostly determine the behaviour of the system and separate them from less prominent players. Here, we tackle this problem making use of three measures of statistical information: resolution, relevance, and mapping entropy. We address the links existing among them, taking the moves from the established relation between resolution and relevance and further developing novel connections between resolution and mapping entropy; by these means we can identify, in a quantitative manner, the number and selection of degrees of freedom of the system that preserve the largest information content about the generative process that underlies an empirical dataset. The method, which is implemented in a freely available software, is fully general, as it is shown through the application to three very diverse systems, namely a toy model of independent binary spins, a coarse-grained representation of the financial stock market, and a fully atomistic simulation of a protein.


title: ApacheJIT: A Large Dataset for Just-In-Time Defect Prediction

abstract: In this paper, we present ApacheJIT, a large dataset for Just-In-Time defect prediction. ApacheJIT consists of clean and bug-inducing software changes in popular Apache projects. ApacheJIT has a total of 106,674 commits (28,239 bug-inducing and 78,435 clean commits). Having a large number of commits makes ApacheJIT a suitable dataset for machine learning models, especially deep learning models that require large training sets to effectively generalize the patterns present in the historical data to future data.


title: First-principle validation of Fourier's law in d=1,2,3 classical systems

abstract: We numerically study the thermal transport in the classical inertial nearest-neighbor XY ferromagnet in $d=1,2,3$, the total number of sites being given by $N=L^d$, where $L$ is the linear size of the system. For the thermal conductance $\sigma$, we obtain $\sigma(T,L)\, L^{\delta(d)} = A(d)\, e_{q(d)}^{- B(d)\,[L^{\gamma(d)}T]^{\eta(d)}}$ (with $e_q^z \equiv [1+(1-q)z]^{1/(1-q)};\,e_1^z=e^z;\,A(d)>0;\,B(d)>0;\,q(d)>1;\,\eta(d)>2;\,\delta \ge 0; \,\gamma(d)>0)$, for all values of $L^{\gamma(d)}T$ for $d=1,2,3$. In the $L\to\infty$ limit, we have $\sigma \propto 1/L^{\rho_\sigma(d)}$ with $\rho_\sigma(d)= \delta(d)+ \gamma(d) \eta(d)/[q(d)-1]$. The material conductivity is given by $\kappa=\sigma L^d \propto 1/L^{\rho_\kappa(d)}$ ($L\to\infty$) with $\rho_\kappa(d)=\rho_\sigma(d)-d$. Our numerical results are consistent with 'conspiratory' $d$-dependences of $(q,\eta,\delta,\gamma)$, which comply with normal thermal conductivity (Fourier law) for all dimensions.


title: The XP Stabiliser Formalism: a Generalisation of the Pauli Stabiliser Formalism with Arbitrary Phases

abstract: We propose an extension to the Pauli stabiliser formalism that includes fractional $2\pi/N$ rotations around the $Z$ axis for some integer $N$. The resulting generalised stabiliser formalism - denoted the XP stabiliser formalism - allows for a wider range of states and codespaces to be represented. We describe the states which arise in the formalism, and demonstrate an equivalence between XP stabiliser states and 'weighted hypergraph states' - a generalisation of both hypergraph and weighted graph states. Given an arbitrary set of XP operators, we present algorithms for determining the codespace and logical operators for an XP code. Finally, we consider whether measurements of XP operators on XP codes can be classically simulated.


title: Entanglement between quasiparticles in superconducting islands mediated by a single spin

abstract: Condensed matter is composed of a small set of identical units, yet it shows an immense range of behaviour. Recently, an array of cold atoms was used to generate long-range quantum entanglement, a property of topological matter. Another approach to strong non-local correlations employs the macroscopic coherence of superconductors. Impurity spins in superconductors are thought to be unamenable to the formation of long-range spin entanglement because each spin tends to be screened by binding to a quasiparticle from the superconductor to form a local singlet. Here we demonstrate that it is possible to attach a second quasiparticle to the spin, overscreening it into a doublet state carrying ferromagnetic correlations between two quasiparticles over a micrometer distance. To demonstrate this effect, which is strongest for equal binding, we symmetrically couple the spin of a quantum dot to two ultrasmall superconducting islands. The overscreened state requires sufficiently large Coulomb repulsion and exchange binding to become well defined. We predict that this state will carry long-range correlations for an alternating chain of quantum dots and superconducting islands, opening a new route to controllable large-scale entanglement in the solid state.


title: Parameterizations of black-hole spacetimes beyond circularity

abstract: We discuss parameterizations of black-hole spacetimes in and beyond General Relativity in view of their symmetry constraints: within the class of axisymmetric, stationary spacetimes, we propose a parameterization that includes non-circular spacetimes, both in Boyer-Lindquist as well as in horizon-penetrating coordinates. We show how existing parameterizations, which make additional symmetry assumptions (first, circularity; second, a hidden constant of motion), are included in the new parameterization. Further, we explain why horizon-penetrating coordinates may be more suitable to parameterize non-circular deviations from the Kerr geometry. Our investigation is motivated by our result that the regular, spinning black-hole spacetimes proposed in [1,2] are non-circular. This particular deviation from circularity can result in cusps, a dent and an asymmetry in the photon rings surrounding the black-hole shadow. Finally, we explore a new class of non-circular deviations from Kerr black holes, which promote the spin parameter to a function, and find indications that regularity cannot be achieved in this setting. This result strengthens the case for regular black holes based on a promotion of the mass parameter to a function.


title: Touching multifunctions on a Hilbert space

abstract: We introduce the concept of the touching of two multifunctions on a real Hilbert space, and deduce that certain multifunctions on the space have a unique fixed point. These result are applied to the theory of genaralized cycles and generalized gap vectors for the composition of the projections onto a finite number of closed convex space in a real Hilbert space.


title: Tunable superconducting coupling of quantum dots via Andreev bound states in semiconductor-superconductor nanowires

abstract: Semiconductor quantum dots have proven to be a useful platform for quantum simulation in the solid state. However, implementing a superconducting coupling between quantum dots mediated by a Cooper pair has so far suffered from limited tunability and strong suppression. This has limited applications such as Cooper pair splitting and quantum dot simulation of topological Kitaev chains. In this work, we propose how to mediate tunable effective couplings via Andreev bound states in a semiconductor-superconductor nanowire connecting two quantum dots. We show that in this way it is possible to individually control both the coupling mediated by Cooper pairs and by single electrons by changing the properties of the Andreev bound states with easily accessible experimental parameters. In addition, the problem of coupling suppression is greatly mitigated. We also propose how to experimentally extract the coupling strengths from resonant current in a three-terminal junction. Our proposal will enable future experiments that have not been possible so far.


title: MRI-GAN: A Generalized Approach to Detect DeepFakes using Perceptual Image Assessment

abstract: DeepFakes are synthetic videos generated by swapping a face of an original image with the face of somebody else. In this paper, we describe our work to develop general, deep learning-based models to classify DeepFake content. We propose a novel framework for using Generative Adversarial Network (GAN)-based models, we call MRI-GAN, that utilizes perceptual differences in images to detect synthesized videos. We test our MRI-GAN approach and a plain-frames-based model using the DeepFake Detection Challenge Dataset. Our plain frames-based-model achieves 91% test accuracy and a model which uses our MRI-GAN framework with Structural Similarity Index Measurement (SSIM) for the perceptual differences achieves 74% test accuracy. The results of MRI-GAN are preliminary and may be improved further by modifying the choice of loss function, tuning hyper-parameters, or by using a more advanced perceptual similarity metric.


title: The first large-scale shell-model calculation of the two-neutrino double beta decay of $^{76}$Ge to the excited states in $^{76}$Se

abstract: Large-scale shell-model calculations were carried out for the half-lives and branching ratios of the $2\nu\beta\beta$ decay of $^{76}$Ge to the ground state and the lowest three excited states $2_1^+$, $0_2^+$ and $2_2^+$ in $^{76}$Se. In total, the wave functions of more than 10,000 intermediate $1^+$ states in $^{76}$As were calculated in a three-step procedure allowing an efficient use of the available computer resources. In the first step, 250 lowest states, below some 5 MeV of excitation energy, were calculated without truncations within a full major shell $0f_{5/2}-1p-0g_{9/2}$ for both protons and neutrons. The wave functions of the rest of the states, up to some 30 MeV, were computed in two more steps by introducing two consecutive stages of truncation. The computed magnitudes of the $2\nu\beta\beta$ nuclear matrix elements (including the value of the axial-vector coupling $g_{\rm A}$), $\vert M_{2\nu}\vert g_{\rm A}^2$, converged to the values 0.168$g_{\rm A}^2$, $1.2\times10^{-3}$$g_{\rm A}^2$, 0.121$g_{\rm A}^2$, and $3.1\times10^{-3}$$g_{\rm A}^2$ for the $0^+_{\rm g.s.}$, $2^+_1$, $0^+_2$, and $2^+_2$ states, respectively. Using up-to-date phase-space integrals, the corresponding branching ratios were derived to be 99.926\%, 4.4$\times10^{-5}$\%, 0.074\% and 2.5$\times10^{-7}$\%. The experimental half-life $(1.926\pm0.094)\times10^{21}$ yr of the ground-state transition was used to derive the value $g_{\rm A}=0.80\pm0.01$ for the axial-vector coupling, which is consistent with other shell-model calculations suggesting a quenched value of $g_{\rm A}$. Using this value of $g_{\rm A}$, predictions for the transition half-lives were derived.


title: An Achievable Rate Region for $3-$User Classical-Quantum Broadcast Channels

abstract: We consider the scenario of communicating on a $3\mhyphen$user classical-quantum broadcast channel. We undertake an information theoretic study and focus on the problem of characterizing an inner bound to its capacity region. We design a new coding scheme based \textit{partitioned coset codes} - an ensemble of codes possessing algebraic properties. Analyzing its information-theoretic performance, we characterize a new inner bound. We identify examples for which the derived inner bound is strictly larger than that achievable using IID random codes.


title: Pedagogical Demonstrations and Pragmatic Learning in Artificial Tutor-Learner Interactions

abstract: When demonstrating a task, human tutors pedagogically modify their behavior by either "showing" the task rather than just "doing" it (exaggerating on relevant parts of the demonstration) or by giving demonstrations that best disambiguate the communicated goal. Analogously, human learners pragmatically infer the communicative intent of the tutor: they interpret what the tutor is trying to teach them and deduce relevant information for learning. Without such mechanisms, traditional Learning from Demonstration (LfD) algorithms will consider such demonstrations as sub-optimal. In this paper, we investigate the implementation of such mechanisms in a tutor-learner setup where both participants are artificial agents in an environment with multiple goals. Using pedagogy from the tutor and pragmatism from the learner, we show substantial improvements over standard learning from demonstrations.


title: GraphWorld: Fake Graphs Bring Real Insights for GNNs

abstract: Despite advances in the field of Graph Neural Networks (GNNs), only a small number (~5) of datasets are currently used to evaluate new models. This continued reliance on a handful of datasets provides minimal insight into the performance differences between models, and is especially challenging for industrial practitioners who are likely to have datasets which look very different from those used as academic benchmarks. In the course of our work on GNN infrastructure and open-source software at Google, we have sought to develop improved benchmarks that are robust, tunable, scalable,and generalizable. In this work we introduce GraphWorld, a novel methodology and system for benchmarking GNN models on an arbitrarily-large population of synthetic graphs for any conceivable GNN task. GraphWorld allows a user to efficiently generate a world with millions of statistically diverse datasets. It is accessible, scalable, and easy to use. GraphWorld can be run on a single machine without specialized hardware, or it can be easily scaled up to run on arbitrary clusters or cloud frameworks. Using GraphWorld, a user has fine-grained control over graph generator parameters, and can benchmark arbitrary GNN models with built-in hyperparameter tuning. We present insights from GraphWorld experiments regarding the performance characteristics of tens of thousands of GNN models over millions of benchmark datasets. We further show that GraphWorld efficiently explores regions of benchmark dataset space uncovered by standard benchmarks, revealing comparisons between models that have not been historically obtainable. Using GraphWorld, we also are able to study in-detail the relationship between graph properties and task performance metrics, which is nearly impossible with the classic collection of real-world benchmarks.


title: Understanding Digits in Identifier Names: An Exploratory Study

abstract: Before any software maintenance can occur, developers must read the identifier names found in the code to be maintained. Thus, high-quality identifier names are essential for productive program comprehension and maintenance activities. With developers free to construct identifier names to their liking, it can be difficult to automatically reason about the quality and semantics behind an identifier name. Studying the structure of identifier names can help alleviate this problem. Existing research focuses on studying words within identifiers, but there are other symbols that appear in identifier names -- such as digits. This paper explores the presence and purpose of digits in identifier names through an empirical study of 800 open-source Java systems. We study how digits contribute to the semantics of identifier names and how identifier names that contain digits evolve over time through renaming. We envision our findings improving the efficiency of name appraisal and recommendation tools and techniques.


title: Tuning the Topological Band Gap of Bismuthene with Silicon-based Substrate

abstract: Some meta-stable polymorphs of bismuth monolayer (bismuthene) can host topologically nontrivial phases. However, it remains unclear if these polymorphs can become stable through interaction with a substrate, whether their topological properties are preserved, and how to design an optimal substrate to make the topological phase more robust. Using first-principles techniques we demonstrate that bismuthene polymorphs can become stable over silicon carbide (SiC), silicon (Si), silicon dioxide (SiO2) and that the proximity interaction in the heterostructures has a significant effect on the electronic structure of the monolayer, even when bonding is weak. We show that the van der Waals interactions and the breaking of the sublattice symmetry are the main factors driving changes in the electronic structure. Our work demonstrates that substrate interaction can strengthen the topological properties of bismuthene polymorphs and make them accessible for experimental investigation and technological applications.


title: The Right Spin: Learning Object Motion from Rotation-Compensated Flow Fields

abstract: Both a good understanding of geometrical concepts and a broad familiarity with objects lead to our excellent perception of moving objects. The human ability to detect and segment moving objects works in the presence of multiple objects, complex background geometry, motion of the observer and even camouflage. How humans perceive moving objects so reliably is a longstanding research question in computer vision and borrows findings from related areas such as psychology, cognitive science and physics. One approach to the problem is to teach a deep network to model all of these effects. This contrasts with the strategy used by human vision, where cognitive processes and body design are tightly coupled and each is responsible for certain aspects of correctly identifying moving objects. Similarly from the computer vision perspective, there is evidence that classical, geometry-based techniques are better suited to the "motion-based" parts of the problem, while deep networks are more suitable for modeling appearance. In this work, we argue that the coupling of camera rotation and camera translation can create complex motion fields that are difficult for a deep network to untangle directly. We present a novel probabilistic model to estimate the camera's rotation given the motion field. We then rectify the flow field to obtain a rotation-compensated motion field for subsequent segmentation. This strategy of first estimating camera motion, and then allowing a network to learn the remaining parts of the problem, yields improved results on the widely used DAVIS benchmark as well as the recently published motion segmentation data set MoCA (Moving Camouflaged Animals).


title: Enhancing Satellite Imagery using Deep Learning for the Sensor To Shooter Timeline

abstract: The sensor to shooter timeline is affected by two main variables: satellite positioning and asset positioning. Speeding up satellite positioning by adding more sensors or by decreasing processing time is important only if there is a prepared shooter, otherwise the main source of time is getting the shooter into position. However, the intelligence community should work towards the exploitation of sensors to the highest speed and effectiveness possible. Achieving a high effectiveness while keeping speed high is a tradeoff that must be considered in the sensor to shooter timeline. In this paper we investigate two main ideas, increasing the effectiveness of satellite imagery through image manipulation and how on-board image manipulation would affect the sensor to shooter timeline. We cover these ideas in four scenarios: Discrete Event Simulation of onboard processing versus ground station processing, quality of information with cloud cover removal, information improvement with super resolution, and data reduction with image to caption. This paper will show how image manipulation techniques such as Super Resolution, Cloud Removal, and Image to Caption will improve the quality of delivered information in addition to showing how those processes effect the sensor to shooter timeline.


title: The Visible Spectro-Polarimeter of the Daniel K. Inouye Solar Telescope

abstract: The Daniel K. Inouye Solar Telescope (DKIST) Visible Spectro-Polarimeter (ViSP) is a traditional slit-scanning spectrograph, with the ability to observe solar regions up to a $120\times78~\mathrm{arcsec}^2$ area. The design implements dual-beam polarimetry, a polychromatic polarization modulator, a high-dispersion echelle grating, and three spectral channels that can be automatically positioned. A defining feature of the instrument is its capability to tune anywhere within the 380-900~nm range of the solar spectrum, allowing for a virtually infinite number of combinations of three wavelengths to be observed simultaneously. This enables the ViSP user to pursue well-established spectro-polarimetric studies of the magnetic structure and plasma dynamics of the solar atmosphere, as well as completely novel investigations of the solar spectrum. Within the suite of first-generation instruments at the DKIST, ViSP is the only wavelength-versatile spectro-polarimeter available to the scientific community. It was specifically designed to be a discovery instrument, for the exploration of new spectroscopic and polarimetric diagnostics, and to test improved models of polarized line formation, through high spatial-, spectral-, and temporal-resolution observations of the Sun's polarized spectrum. In this instrument article, we describe the science requirements and design drivers of ViSP, and we present preliminary science data collected during the commissioning of the instrument.


title: A Gelfand Transform for Spinor Fields on Embedded Riemannian Manifolds

abstract: A classical result of Gelfand shows that the topologized spectrum of characters on commutative Banach algebra is homeomorphic to the underlying space. This fact is used in solving the Calder\'on problem in dimension 2 via the boundary control (BC) method. To apply the BC method in dimension 3, the algebra of complex holomorphic functions can be replaced by the space of harmonic quaternion fields, but this space is no longer an algebra and is not commutative. Nonetheless, it has been shown that a suitable notion of a spectrum exists and in the case when the underlying space is convex, the spectrum is homeomorphic to the ball. Our goal is to generalize this result to more general manifolds in arbitrary dimension. To do so, we use Clifford algebras of multivector fields and the set of functions we care about is the space monogenic spinor fields. The spectrum consists of spinor valued functionals that respect the module and subalgebra structure of the space of monogenic spinor fields. We prove that for compact regions in $\mathbb{R}^n$, the spectrum is homeomorphic to the manifold itself. Furthermore, some essential lemmas and propositions are proven for arbitrary compact Riemannian $n$-manifolds. Finally, we end with a Stone--Weierstrass theorem which shows the algebra generated by the monogenic spinor fields on arbitrary compact $M$ with boundary is dense in the algebra of continuous spinor fields.


title: Efficient Task Allocation in Smart Warehouses with Multi-delivery Stations and Heterogeneous Robots

abstract: The task allocation problem in multi-robot systems (MRTA) is an NP-hard problem whose viable solutions are usually found by heuristic algorithms. Considering the increasing need of improvement on logistics, the use of robots for increasing the efficiency of logistics warehouses is becoming a requirement. In a smart warehouse the main tasks consist of employing a fleet of automated picking and mobile robots that coordinate by picking up items from a set of orders from the shelves and dropping them at the delivery stations. Two aspects generally justify multi-robot task allocation complexity: (i) environmental aspects, such as multi-delivery stations and dispersed robots (since they remain in constant motion) and (ii) fleet heterogeneity, where robots' traffic speed and capacity loads are different from each other. Despite these properties have been widely researched in the literature, they usually are investigated separately. This work proposes a scalable and efficient task allocation algorithm for smart warehouses with multi-delivery stations and heterogeneous fleets. Our strategy employs a novel cost estimator, which computes costs as a function of the robots' variable characteristics and capacity while they receive new tasks. For validating the strategy a series of experiments is performed simulating the operation of smart warehouses with multiple delivery stations and heteregenous fleets. The results show that our strategy generates routes costing up to $33$\% less than the routes generated by a state-of-the-art task allocation algorithm and $96$\% faster in test instances representing our target scenario. Considering single-delivery stations and non-dispersed robots, we reduced the number of robots by up to $18$\%, allocating tasks $92$\% faster, and generating routes whose costs are statistically similar to the routes generated by the state-of-the-art algorithm.


title: Neural Ordinary Differential Equations for Nonlinear System Identification

abstract: Neural ordinary differential equations (NODE) have been recently proposed as a promising approach for nonlinear system identification tasks. In this work, we systematically compare their predictive performance with current state-of-the-art nonlinear and classical linear methods. In particular, we present a quantitative study comparing NODE's performance against neural state-space models and classical linear system identification methods. We evaluate the inference speed and prediction performance of each method on open-loop errors across eight different dynamical systems. The experiments show that NODEs can consistently improve the prediction accuracy by an order of magnitude compared to benchmark methods. Besides improved accuracy, we also observed that NODEs are less sensitive to hyperparameters compared to neural state-space models. On the other hand, these performance gains come with a slight increase of computation at the inference time.


title: Dense Gas History of the Universe: from ASPECS to the ngVLA

abstract: We review the evolution of the cosmic average molecular gas density to large look-back times, using observations of rotational transitions of CO. Molecular gas is the fuel for star formation in galaxies. Deep searches for CO emission from distant galaxies have delineated the density of molecular gas back to $z \sim 5$, or within 1~Gyr of the Big Bang. The results show a rise and fall in the gas density that parallels, and likely drives, the rise and fall of the cosmic star formation rate density. We present the potential for the next generation Very Large Array to image the distribution and dynamics of the molecular gas in early galaxies, and to make a precise measurement of the dense gas history of the Universe.


title: Uniqueness for nonlinear Fokker-Planck equations and for McKean-Vlasov SDEs: The degenerate case

abstract: This work is concerned with the existence and uniqueness of generalized (mild or distributional) solutions to (possibly degenerate) Fokker-Planck equations $\rho_t-\Delta\beta(\rho)+{\rm div}(Db(\rho)\rho)=0$ in $(0,\infty)\times\mathbb{R}^d,$ $\rho(0,x) \equiv \rho_0(x)$. Under suitable assumptions on $\beta:\mathbb{R}\to\mathbb{R},\,b:\mathbb{R}\to\mathbb{R}$ and $D:\mathbb{R}^d\to\mathbb{R}^d$, $d\ge1$, this equation generates a unique flow $\rho(t)=S(t)\rho_0:[0,\infty)\to L^1(\mathbb{R}^d)$ as a mild solution in the sense of nonlinear semigroup theory. This flow is also unique in the class of $L^\infty((0,T)\times\mathbb{R}^d)\cap L^1((0,T)\times\mathbb{R}^d),$ $\forall T>0$, Schwartz distributional solutions on $(0,\infty)\times\mathbb{R}^d$. Moreover, for $\rho_0\in L^1(\mathbb{R}^d)\cap H^{-1}(\mathbb{R}^d)$, $t\to S(t)\rho_0$ is differentiable from the right on $[0,\infty)$ in $H^{-1}(\mathbb{R}^d)$-norm. As a main application, the weak uniqueness of the corresponding McKean-Vlasov SDEs is proven.


title: Rectifying homographies for stereo vision: analytical solution for minimal distortion

abstract: Stereo rectification is the determination of two image transformations (or homographies) that map corresponding points on the two images, projections of the same point in the 3D space, onto the same horizontal line in the transformed images. Rectification is used to simplify the subsequent stereo correspondence problem and speeding up the matching process. Rectifying transformations, in general, introduce perspective distortion on the obtained images, which shall be minimised to improve the accuracy of the following algorithm dealing with the stereo correspondence problem. The search for the optimal transformations is usually carried out relying on numerical optimisation. This work proposes a closed-form solution for the rectifying homographies that minimise perspective distortion. The experimental comparison confirms its capability to solve the convergence issues of the previous formulation. Its Python implementation is provided.


title: On classification of strategic agents who can both game and improve

abstract: In this work, we consider classification of agents who can both game and improve. For example, people wishing to get a loan may be able to take some actions that increase their perceived credit-worthiness and others that also increase their true credit-worthiness. A decision-maker would like to define a classification rule with few false-positives (does not give out many bad loans) while yielding many true positives (giving out many good loans), which includes encouraging agents to improve to become true positives if possible. We consider two models for this problem, a general discrete model and a linear model, and prove algorithmic, learning, and hardness results for each. For the general discrete model, we give an efficient algorithm for the problem of maximizing the number of true positives subject to no false positives, and show how to extend this to a partial-information learning setting. We also show hardness for the problem of maximizing the number of true positives subject to a nonzero bound on the number of false positives, and that this hardness holds even for a finite-point version of our linear model. We also show that maximizing the number of true positives subject to no false positive is NP-hard in our full linear model. We additionally provide an algorithm that determines whether there exists a linear classifier that classifies all agents accurately and causes all improvable agents to become qualified, and give additional results for low-dimensional data.