title: Average-case analysis of the Gaussian Elimination with Partial Pivoting

abstract: The Gaussian Elimination with Partial Pivoting (GEPP) is a classical algorithm for solving systems of linear equations. Although in specific cases the loss of precision in GEPP due to roundoff errors can be very significant, empirical evidence strongly suggests that for a {\it typical} square coefficient matrix, GEPP is numerically stable. We obtain a (partial) theoretical justification of this phenomenon by showing that, given the random $n\times n$ standard Gaussian coefficient matrix $A$, the {\it growth factor} of the Gaussian Elimination with Partial Pivoting is at most polynomially large in $n$ with probability close to one. This implies that with probability close to one the number of bits of precision sufficient to solve $Ax = b$ to $m$ bits of accuracy using GEPP is $m+O(\log n)$, which improves an earlier estimate $m+O(\log^2 n)$ of Sankar, and which we conjecture to be optimal by the order of magnitude. We further provide tail estimates of the growth factor which can be used to support the empirical observation that GEPP is more stable than the Gaussian Elimination with no pivoting.


title: New Progress in Classic Area: Polynomial Root-squaring and Root-finding

abstract: The DLG root-squaring iterations, due to Dandelin 1826 and rediscovered by Lobachevsky 1834 and Graeffe 1837, have been the main approach to root-finding for a univariate polynomial p(x) in the 19th century and beyond, but not so nowadays because these iterations are prone to severe numerical stability problems. Trying to avoid these problems we have found simple but novel reduction of the iterations applied for Newton's inverse ratio -p'(x)/p(x) to approximation of the power sums of the zeros of p(x) and its reverse polynomial. The resulting polynomial root-finders can be devised and performed independently of DLG iterations, based on Newton's identities or Cauchy integrals. In the former case the computation involve a set of leading or tailing coefficients of an input polynomial. In the latter case we must scale the variable and increase the arithmetic computational cost to ensure numerical stability. Nevertheless the cost is still manageable, at least for fast root-refinement, and the algorithms can be applied to a black box polynomial p(x)--given by a black box for the evaluation of the ratio p'(x)/p(x) rather than by its coefficients. This enables important computational benefits, including efficient recursive as well as concurrent approximation of a set of zeros of p(x) or even all of its zeros, acceleration where an input polynomial can be evaluated fast, and extension to approximation of the eigenvalues of a matrix or a polynomial matrix, being efficient if the matrix can be inverted fast, e.g., is data sparse. We also recall our recent fast algorithms for approximation of the root radii, that is, the distances to the roots from the origin or any complex value to the zeros of p(x), and apply it for fast black box initialization of polynomial root-finding by means of functional iterations, e.g., Newton's, Ehrlich's, and Weierstrass's.


title: A review of machine learning approaches, challenges and prospects for computational tumor pathology

abstract: Computational pathology is part of precision oncology medicine. The integration of high-throughput data including genomics, transcriptomics, proteomics, metabolomics, pathomics, and radiomics into clinical practice improves cancer treatment plans, treatment cycles, and cure rates, and helps doctors open up innovative approaches to patient prognosis. In the past decade, rapid advances in artificial intelligence, chip design and manufacturing, and mobile computing have facilitated research in computational pathology and have the potential to provide better-integrated solutions for whole-slide images, multi-omics data, and clinical informatics. However, tumor computational pathology now brings some challenges to the application of tumour screening, diagnosis and prognosis in terms of data integration, hardware processing, network sharing bandwidth and machine learning technology. This review investigates image preprocessing methods in computational pathology from a pathological and technical perspective, machine learning-based methods, and applications of computational pathology in breast, colon, prostate, lung, and various tumour disease scenarios. Finally, the challenges and prospects of machine learning in computational pathology applications are discussed.


title: Torsional Diffusion for Molecular Conformer Generation

abstract: Molecular conformer generation is a fundamental task in computational chemistry. Several machine learning approaches have been developed, but none have outperformed state-of-the-art cheminformatics methods. We propose torsional diffusion, a novel diffusion framework that operates on the space of torsion angles via a diffusion process on the hypertorus and an extrinsic-to-intrinsic score model. On a standard benchmark of drug-like molecules, torsional diffusion generates superior conformer ensembles compared to machine learning and cheminformatics methods in terms of both RMSD and chemical properties, and is orders of magnitude faster than previous diffusion-based models. Moreover, our model provides exact likelihoods, which we employ to build the first generalizable Boltzmann generator. Code is available at https://github.com/gcorso/torsional-diffusion.


title: On the complexity of nonsmooth automatic differentiation

abstract: Using the notion of conservative gradient, we provide a simple model to estimate the computational costs of the backward and forward modes of algorithmic differentiation for a wide class of nonsmooth programs. The overhead complexity of the backward mode turns out to be independent of the dimension when using programs with locally Lipschitz semi-algebraic or definable elementary functions. This considerably extends Baur-Strassen's smooth cheap gradient principle. We illustrate our results by establishing fast backpropagation results of conservative gradients through feedforward neural networks with standard activation and loss functions. Nonsmooth backpropagation's cheapness contrasts with concurrent forward approaches, which have, to this day, dimensional-dependent worst-case overhead estimates. We provide further results suggesting the superiority of backward propagation of conservative gradients. Indeed, we relate the complexity of computing a large number of directional derivatives to that of matrix multiplication, and we show that finding two subgradients in the Clarke subdifferential of a function is an NP-hard problem.


title: Empirical Study of Quality Image Assessment for Synthesis of Fetal Head Ultrasound Imaging with DCGANs

abstract: In this work, we present an empirical study of DCGANs, including hyperparameter heuristics and image quality assessment, as a way to address the scarcity of datasets to investigate fetal head ultrasound. We present experiments to show the impact of different image resolutions, epochs, dataset size input, and learning rates for quality image assessment on four metrics: mutual information (MI), Fr\'echet inception distance (FID), peak-signal-to-noise ratio (PSNR), and local binary pattern vector (LBPv). The results show that FID and LBPv have stronger relationship with clinical image quality scores. The resources to reproduce this work are available at \url{https://github.com/budai4medtech/miua2022}.


title: Linear-Quadratic Mean Field Games of Controls with Non-Monotone Data

abstract: In this paper, we study a class of linear-quadratic (LQ) mean field games of controls with common noises and their corresponding $N$-player games. The theory of mean field game of controls considers a class of mean field games where the interaction is via the joint law of both the state and control. By the stochastic maximum principle, we first analyze the limiting behavior of the representative player and obtain his/her optimal control in a feedback form with the given distributional flow of the population and its control. The mean field equilibrium is determined by the Nash certainty equivalence (NCE) system. Thanks to the common noise, we do not require any monotonicity conditions for the solvability of the NCE system. We also study the master equation arising from LQ mean field games of controls, which is a finite-dimensional second-order parabolic equation. It can be shown that the master equation admits a unique classical solution over an arbitrary time horizon without any monotonicity conditions. Beyond that, we can solve the $N$-player games directly by further assuming the non-degeneracy of the idiosyncratic noises. As byproducts, we prove the quantitative convergence results from the $N$-player game to the mean field game and the propagation of chaos property for the related optimal trajectories.


title: Adversarial RAW: Image-Scaling Attack Against Imaging Pipeline

abstract: Deep learning technologies have become the backbone for the development of computer vision. With further explorations, deep neural networks have been found vulnerable to well-designed adversarial attacks. Most of the vision devices are equipped with image signal processing (ISP) pipeline to implement RAW-to-RGB transformations and embedded into data preprocessing module for efficient image processing. Actually, ISP pipeline can introduce adversarial behaviors to post-capture images while data preprocessing may destroy attack patterns. However, none of the existing adversarial attacks takes into account the impacts of both ISP pipeline and data preprocessing. In this paper, we develop an image-scaling attack targeting on ISP pipeline, where the crafted adversarial RAW can be transformed into attack image that presents entirely different appearance once being scaled to a specific-size image. We first consider the gradient-available ISP pipeline, i.e., the gradient information can be directly used in the generation process of adversarial RAW to launch the attack. To make the adversarial attack more applicable, we further consider the gradient-unavailable ISP pipeline, in which a proxy model that well learns the RAW-to-RGB transformations is proposed as the gradient oracles. Extensive experiments show that the proposed adversarial attacks can craft adversarial RAW data against the target ISP pipelines with high attack rates.


title: Using UAS Imagery and Computer Vision to Support Site-Specific Weed Control in Corn

abstract: Currently, weed control in a corn field is performed by a blanket application of herbicides that do not consider spatial distribution information of weeds and also uses an extensive amount of chemical herbicides. To reduce the amount of chemicals, we used drone-based high-resolution imagery and computer-vision techniques to perform site-specific weed control in corn.


title: Examining the behaviour of state-of-the-art convolutional neural networks for brain tumor detection with and without transfer learning

abstract: Distinguishing normal from malignant and determining the tumor type are critical components of brain tumor diagnosis. Two different kinds of dataset are investigated using state-of-the-art CNN models in this research work. One dataset(binary) has images of normal and tumor types, while another(multi-class) provides all images of tumors classified as glioma, meningioma, or pituitary. The experiments were conducted in these dataset with transfer learning from pre-trained weights from ImageNet as well as initializing the weights randomly. The experimental environment is equivalent for all models in this study in order to make a fair comparison. For both of the dataset, the validation set are same for all the models where train data is 60% while the rest is 40% for validation. With the proposed techniques in this research, the EfficientNet-B5 architecture outperforms all the state-of-the-art models in the binary-classification dataset with the accuracy of 99.75% and 98.61% accuracy for the multi-class dataset. This research also demonstrates the behaviour of convergence of validation loss in different weight initialization techniques.


title: Adaptive Adversarial Training to Improve Adversarial Robustness of DNNs for Medical Image Segmentation and Detection

abstract: It is known that Deep Neural networks (DNNs) are vulnerable to adversarial attacks, and the adversarial robustness of DNNs could be improved by adding adversarial noises to training data (e.g., the standard adversarial training (SAT)). However, inappropriate noises added to training data may reduce a model's performance, which is termed the trade-off between accuracy and robustness. This problem has been sufficiently studied for the classification of whole images but has rarely been explored for image analysis tasks in the medical application domain, including image segmentation, landmark detection, and object detection tasks. In this study, we show that, for those medical image analysis tasks, the SAT method has a severe issue that limits its practical use: it generates a fixed and unified level of noise for all training samples for robust DNN training. A high noise level may lead to a large reduction in model performance and a low noise level may not be effective in improving robustness. To resolve this issue, we design an adaptive-margin adversarial training (AMAT) method that generates sample-wise adaptive adversarial noises for robust DNN training. In contrast to the existing, classification-oriented adversarial training methods, our AMAT method uses a loss-defined-margin strategy so that it can be applied to different tasks as long as the loss functions are well-defined. We successfully apply our AMAT method to state-of-the-art DNNs, using five publicly available datasets. The experimental results demonstrate that: (1) our AMAT method can be applied to the three seemingly different tasks in the medical image application domain; (2) AMAT outperforms the SAT method in adversarial robustness; (3) AMAT has a minimal reduction in prediction accuracy on clean data, compared with the SAT method; and (4) AMAT has almost the same training time cost as SAT.


title: MaxStyle: Adversarial Style Composition for Robust Medical Image Segmentation

abstract: Convolutional neural networks (CNNs) have achieved remarkable segmentation accuracy on benchmark datasets where training and test sets are from the same domain, yet their performance can degrade significantly on unseen domains, which hinders the deployment of CNNs in many clinical scenarios. Most existing works improve model out-of-domain (OOD) robustness by collecting multi-domain datasets for training, which is expensive and may not always be feasible due to privacy and logistical issues. In this work, we focus on improving model robustness using a single-domain dataset only. We propose a novel data augmentation framework called MaxStyle, which maximizes the effectiveness of style augmentation for model OOD performance. It attaches an auxiliary style-augmented image decoder to a segmentation network for robust feature learning and data augmentation. Importantly, MaxStyle augments data with improved image style diversity and hardness, by expanding the style space with noise and searching for the worst-case style composition of latent features via adversarial training. With extensive experiments on multiple public cardiac and prostate MR datasets, we demonstrate that MaxStyle leads to significantly improved out-of-distribution robustness against unseen corruptions as well as common distribution shifts across multiple, different, unseen sites and unknown image sequences under both low- and high-training data settings. The code can be found at https://github.com/cherise215/MaxStyle.


title: RIDDLE: Lidar Data Compression with Range Image Deep Delta Encoding

abstract: Lidars are depth measuring sensors widely used in autonomous driving and augmented reality. However, the large volume of data produced by lidars can lead to high costs in data storage and transmission. While lidar data can be represented as two interchangeable representations: 3D point clouds and range images, most previous work focus on compressing the generic 3D point clouds. In this work, we show that directly compressing the range images can leverage the lidar scanning pattern, compared to compressing the unprojected point clouds. We propose a novel data-driven range image compression algorithm, named RIDDLE (Range Image Deep DeLta Encoding). At its core is a deep model that predicts the next pixel value in a raster scanning order, based on contextual laser shots from both the current and past scans (represented as a 4D point cloud of spherical coordinates and time). The deltas between predictions and original values can then be compressed by entropy encoding. Evaluated on the Waymo Open Dataset and KITTI, our method demonstrates significant improvement in the compression rate (under the same distortion) compared to widely used point cloud and range image compression algorithms as well as recent deep methods.


title: Mutual- and Self- Prototype Alignment for Semi-supervised Medical Image Segmentation

abstract: Semi-supervised learning methods have been explored in medical image segmentation tasks due to the scarcity of pixel-level annotation in the real scenario. Proto-type alignment based consistency constraint is an intuitional and plausible solu-tion to explore the useful information in the unlabeled data. In this paper, we propose a mutual- and self- prototype alignment (MSPA) framework to better utilize the unlabeled data. In specific, mutual-prototype alignment enhances the information interaction between labeled and unlabeled data. The mutual-prototype alignment imposes two consistency constraints in reverse directions between the unlabeled and labeled data, which enables the consistent embedding and model discriminability on unlabeled data. The proposed self-prototype alignment learns more stable region-wise features within unlabeled images, which optimizes the classification margin in semi-supervised segmentation by boosting the intra-class compactness and inter-class separation on the feature space. Extensive experimental results on three medical datasets demonstrate that with a small amount of labeled data, MSPA achieves large improvements by leveraging the unlabeled data. Our method also outperforms seven state-of-the-art semi-supervised segmentation methods on all three datasets.


title: Denoising Fast X-Ray Fluorescence Raster Scans of Paintings

abstract: Macro x-ray fluorescence (XRF) imaging of cultural heritage objects, while a popular non-invasive technique for providing elemental distribution maps, is a slow acquisition process in acquiring high signal-to-noise ratio XRF volumes. Typically on the order of tenths of a second per pixel, a raster scanning probe counts the number of photons at different energies emitted by the object under x-ray illumination. In an effort to reduce the scan times without sacrificing elemental map and XRF volume quality, we propose using dictionary learning with a Poisson noise model as well as a color image-based prior to restore noisy, rapidly acquired XRF data.


title: Patcher: Patch Transformers with Mixture of Experts for Precise Medical Image Segmentation

abstract: We present a new encoder-decoder Vision Transformer architecture, Patcher, for medical image segmentation. Unlike standard Vision Transformers, it employs Patcher blocks that segment an image into large patches, each of which is further divided into small patches. Transformers are applied to the small patches within a large patch, which constrains the receptive field of each pixel. We intentionally make the large patches overlap to enhance intra-patch communication. The encoder employs a cascade of Patcher blocks with increasing receptive fields to extract features from local to global levels. This design allows Patcher to benefit from both the coarse-to-fine feature extraction common in CNNs and the superior spatial relationship modeling of Transformers. We also propose a new mixture-of-experts (MoE) based decoder, which treats the feature maps from the encoder as experts and selects a suitable set of expert features to predict the label for each pixel. The use of MoE enables better specializations of the expert features and reduces interference between them during inference. Extensive experiments demonstrate that Patcher outperforms state-of-the-art Transformer- and CNN-based approaches significantly on stroke lesion segmentation and polyp segmentation. Code for Patcher is released with publication to facilitate future research.


title: Learning Probabilistic Topological Representations Using Discrete Morse Theory

abstract: Accurate delineation of fine-scale structures is a very important yet challenging problem. Existing methods use topological information as an additional training loss, but are ultimately making pixel-wise predictions. In this paper, we propose the first deep learning based method to learn topological/structural representations. We use discrete Morse theory and persistent homology to construct an one-parameter family of structures as the topological/structural representation space. Furthermore, we learn a probabilistic model that can perform inference tasks in such a topological/structural representation space. Our method generates true structures rather than pixel-maps, leading to better topological integrity in automatic segmentation tasks. It also facilitates semi-automatic interactive annotation/proofreading via the sampling of structures and structure-aware uncertainty.


title: Orthogonal Transform based Generative Adversarial Network for Image Dehazing

abstract: Image dehazing has become one of the crucial preprocessing steps for any computer vision task. Most of the dehazing methods try to estimate the transmission map along with the atmospheric light to get the dehazed image in the image domain. In this paper, we propose a novel end-to-end architecture that directly estimates dehazed image in Krawtchouk transform domain. For this a customized Krawtchouk Convolution Layer (KCL) in the architecture is added. KCL is constructed using Krawtchouk basis functions which converts the image from the spatial domain to the Krawtchouk transform domain. Another convolution layer is added at the end of the architecture named as Inverse Krawtchouk Convolution Layer (IKCL) which converts the image back to the spatial domain from the transform domain. It has been observed that the haze is mainly present in lower frequencies of hazy images, wherein the Krawtchouk transform helps to analyze the high and low frequencies of the images separately. We have divided our architecture into two branches, the upper branch deals with the higher frequencies while the lower branch deals with the lower frequencies of the image. The lower branch is made deeper in terms of the layers as compared to the upper branch to address the haze present in the lower frequencies. Using the proposed Orthogonal Transform based Generative Adversarial Network (OTGAN) architecture for image dehazing, we were able to achieve competitive results when compared to the present state-of-the-art methods.


title: Complementary ADF-STEM: a Flexible Approach to Quantitative 4D-STEM

abstract: Scanning transmission electron microscopy (STEM) has a broad range of applications in materials characterization, including real-space imaging, spectroscopy, and diffraction, at length scales from the micron to sub-{\AA}ngstr\"om. The recent development and adoption of high-speed, direct electron STEM detectors has enabled diffraction patterns to be collected at each probe position, generating four-dimensional STEM (4D-STEM) datasets and opening new imaging modalities. However, the limited pixel numbers in these detectors enforce a tradeoff between angular resolution and maximum collection angle. In this paper, we describe a straightforward method for quantifying 4D-STEM data by utilizing the full flux of the electron beam, including electrons scattered beyond the limits of the detector. This enables significantly increased experimental flexibility, including the synthesis of quantitative, high-contrast complementary annular dark field (cADF) STEM images from low-angle diffraction patterns whilst maintaining high angular resolution; as well as the optimization of electron dose and the more effective use of low dynamic range detectors.


title: Detection of Fibrosis in Cine Magnetic Resonance Images Using Artificial Intelligence Techniques

abstract: Background: Artificial intelligence techniques have demonstrated great potential in cardiology, especially to detect imperceptible patterns for the human eye. In this sense, these techniques seem to be adequate to identify patterns in the myocardial texture which could lead to characterize and quantify fibrosis. Purpose: The aim of this study was to postulate a new artificial intelligence method to identify fibrosis in cine cardiac magnetic resonance (CMR) imaging. Methods: A retrospective observational study was carried out in a population of 75 subjects from a clinical center of San Carlos de Bariloche. The proposed method analyzes the myocardial texture in cine CMR images using a convolutional neural network to determine local myocardial tissue damage. Results: An accuracy of 89% for quantifying local tissue damage was observed for the validation data set and 70% for the test set. In addition, the qualitative analysis showed a high spatial correlation in lesion location. Conclusions: The postulated method enables to spatially identify fibrosis using only the information from cine nuclear magnetic resonance studies, demonstrating the potential of this technique to quantify myocardial viability in the future or to study the lesions etiology


title: Automatic Quantification of Volumes and Biventricular Function in Cardiac Resonance. Validation of a New Artificial Intelligence Approach

abstract: Background: Artificial intelligence techniques have shown great potential in cardiology, especially in quantifying cardiac biventricular function, volume, mass, and ejection fraction (EF). However, its use in clinical practice is not straightforward due to its poor reproducibility with cases from daily practice, among other reasons. Objectives: To validate a new artificial intelligence tool in order to quantify the cardiac biventricular function (volume, mass, and EF). To analyze its robustness in the clinical area, and the computational times compared with conventional methods. Methods: A total of 189 patients were analyzed: 89 from a regional center and 100 from a public center. The method proposes two convolutional networks that include anatomical information of the heart to reduce classification errors. Results: A high concordance (Pearson coefficient) was observed between manual quantification and the proposed quantification of cardiac function (0.98, 0.92, 0.96 and 0.8 for volumes and biventricular EF) in about 5 seconds per study. Conclusions: This method quantifies biventricular function and volumes in seconds with an accuracy equivalent to that of a specialist.

categories: eess.IV cs.AI cs.CV cs.LG physics.med-ph


title: A field equation for induction-transduction of activation-deactivation probability on measurable space

abstract: Induction-transduction of activating-deactivating points are fundamental mechanisms of action that underlie innumerable systems and phenomena, mathematical, natural, and anthropogenic, and can exhibit complex behaviors such as self-excitation, phase transitions, hysteresis, polarization, periodicity, chaos, wave behavior, geometry, and energy transfer. We describe a class of primitives for induction-transduction based on dynamics on images of marked random counting measures under graphical random transformations. We derive a field equation for the law of the activation-deactivation (Bernoulli) process on an arbitrary measurable space and describe some mechanisms of action on the unit interval.


title: Federated Deep Learning Meets Autonomous Vehicle Perception: Design and Verification

abstract: Realizing human-like perception is a challenge in open driving scenarios due to corner cases and visual occlusions. To gather knowledge of rare and occluded instances, federated learning assisted connected autonomous vehicle (FLCAV) has been proposed, which leverages vehicular networks to establish federated deep neural networks (DNNs) from distributed data captured by vehicles and road sensors. Without the need of data aggregation, FLCAV preserves privacy while reducing communication costs compared with conventional centralized learning. However, it is challenging to determine the network resources and road sensor placements for multi-stage training with multi-modal datasets in multi-variant scenarios. This article presents networking and training frameworks for FLCAV perception. Multi-layer graph resource allocation and vehicle-road contrastive sensor placement are proposed to address the network management and sensor deployment problems, respectively. We also develop CarlaFLCAV, a software platform that implements the above system and methods. Experimental results confirm the superiority of the proposed techniques compared with various benchmarks.


title: Uncertainty Estimation in Machine Learning

abstract: Most machine learning techniques are based upon statistical learning theory, often simplified for the sake of computing speed. This paper is focused on the uncertainty aspect of mathematical modeling in machine learning. Regression analysis is chosen to further investigate the evaluation aspect of uncertainty in model coefficients and, more importantly, in the output feature value predictions. A survey demonstrates major stages in the conventional least squares approach to the creation of the regression model, along with its uncertainty estimation. On the other hand, it is shown that in machine learning the model complexity and severe nonlinearity become serious obstacles to uncertainty evaluation. Furthermore, the process of machine model training demands high computing power, not available at the level of personal computers. This is why so-called pre-trained models are widely used in such areas of machine learning as natural language processing. The latest example of a pre-trained model is the Generative Pre-trained Transformer 3 with hundreds of billions of parameters and a half-terabyte training dataset. Similarly, mathematical models built from real data are growing in complexity which is accompanied by the growing amount of training data. However, when machine models and their predictions are used in decision-making, one needs to estimate uncertainty and evaluate accompanying risks. This problem could be resolved with non-parametric techniques at the expense of greater demand for computing power, which can be offered by modern supercomputers available, including those utilizing graphical and tensor processing units along with the conventional central processors.


title: Implications of a Temperature Dependent IMF III: Mass Growth and Quiescence

abstract: The stellar initial mass function (IMF) is predicted to depend upon the temperature of gas in star-forming molecular clouds. The introduction of an additional parameter, $T_{IMF}$ , into photometric template fitting, suggest most galaxies obey an IMF top-heavier than the Galactic IMF. The implications of these revised fits on mass functions, quiescence and turnoff are discussed. At all redshifts the highest mass galaxies become quiescent first with the turnoff mass decreasing towards the present. The synchronous turnoff mass across galaxies suggests quiescence is driven by universal mechanisms rather than by stochastic or environmental processes.


title: Grand Unification of continuous-variable codes

abstract: Quantum error correction codes in continuous variables (also called CV codes, or single-mode bosonic codes) have recently been identified to be a technologically viable option for building fault-tolerant quantum computers. The best-known examples are the GKP code and the cat-code, both of which were shown to have some advantageous properties over any discrete-variable, or qubit codes. It was recently shown that the cat-code, as well as other kinds of CV codes, belong to a set of codes with common properties called rotation-symmetric codes. We expand this result by giving a general description of sets of codes with common properties, and rules by which they can be mapped into one another, effectively creating a unified description of continuous-variable codes. We prove that the properties of all of these sets of codes can be obtained from the properties of the GKP code. We also show explicitly how this construction works in the case of rotation-symmetric codes, re-deriving known properties and finding new ones.


title: Double brane holographic model dual to 2d ICFTs

abstract: A minimal single-brane holographic model can be used as a dual to 2d conformal interfaces (ICFTs) to calculate the transmission coefficient $\mathcal{T}$ of energy transported across the defect as well as boundary entropy $\log g$, the additional entanglement entropy for some sub-region that encloses the defect. Both $\mathcal{T}$ and $\log g$ are uniquely determined by the tension characterizing the brane. In contrast, in field theory defects typically the transmission coefficient can be dialed from 0 to 1 independently for each allowed value of $\log g$. To address this discrepancy, we look at a double brane (3-region bulk) holographic model. Merger of two single brane interfaces creates genuinely new interfaces which indeed allow a range of accessible transmission coefficients for a fixed value of $\log g$. In particular, the $\mathcal{T}=0$ limit of two completely decoupled BCFTs can be achieved.


title: Water emission tracing active star formation from the Milky Way to high-$z$ galaxies

abstract: (Abridged) The question of how most stars in the Universe form remains open. While star formation predominantly occurs in young massive clusters, the current framework focuses on isolated star formation. One way to access the bulk of protostellar activity within star-forming clusters is to trace signposts of active star formation with emission from molecular outflows. These outflows are bright in water emission, providing a direct observational link between nearby and distant galaxies. We propose to utilize the knowledge of local star formation as seen with molecular tracers to explore the nature of star formation in the Universe. We present a large-scale statistical galactic model of emission from galactic active star-forming regions. Our model is built on observations of well-resolved nearby clusters. By simulating emission from molecular outflows, which is known to scale with mass, we create a proxy that can be used to predict the emission from clustered star formation at galactic scales. We evaluated the impact of the most important global-star formation parameters (i.e., initial stellar mass function (IMF), molecular cloud mass distribution, star formation efficiency (SFE), and free-fall time efficiency) on simulation results. We observe that for emission from the para-H2O 202 - 111 line, the IMF and molecular cloud mass distribution have a negligible impact on the emission, both locally and globally, whereas the opposite holds for the SFE and free-fall time efficiency. Moreover, this water transition proves to be a low-contrast tracer of star formation. The fine-tuning of the model and adaptation to morphologies of distant galaxies should result in realistic predictions of observed molecular emission and make the galaxy-in-a-box model a tool to analyze and better understand star formation throughout cosmological times.


title: Radio Emission from Binary Ultracool Dwarf Systems

abstract: Well-characterized binary systems will provide valuable opportunities to study the conditions that are necessary for the onset of both auroral and non-auroral magnetospheric radio emission in the ultracool dwarf regime. We present new detections of non-auroral "quiescent" radio emission at 4-8 GHz of the three ultracool dwarf binary systems GJ 564 BC, LP 415-20, and 2MASS J21402931+1625183. We also tentatively detect a highly circularly polarized pulse at 4-6 GHz that may indicate aurorae from GJ 564 BC. Finally, we show that the brightest binary ultracool dwarf systems may be more luminous than predictions from single-object systems.


title: Dynamical Instability in Multi-Orbiter Systems with Gas Friction

abstract: Closely-packed multi-planet systems are known to experience dynamical instability if the spacings between the planets are too small. Such instability can be tempered by the frictional forces acting on the planets from gaseous discs. A similar situation applies to stellar-mass black holes embedded in AGN discs around supermassive black holes. In this paper, we use $N$-body integrations to evaluate how the frictional damping of orbital eccentricity affects the growth of dynamical instability for a wide range of planetary spacing and planet-to-star mass ratios. We find that the stability of a system depends on the damping timescale $\tau$ relative to the zero-friction instability growth timescale $t_{\rm inst}$. In a two-planet system, the frictional damping can stabilise the dynamical evolution if $t_{\rm inst}\gtrsim\tau$. With three planets, $t_{\rm inst} \gtrsim 10\tau - 100\tau$ is needed for stabilisation. When the separations between the planetary orbits are sufficiently small, $t_{\rm inst}$ can be less than the synodic period between the planets, which makes frictional stabilisation unlikely to occur. As the orbital spacing increases, the instability timescale tends to grow exponentially on average, but it can vary by a few orders of magnitude depending on the initial orbital phases of the planets. In general, the stable region (at large orbital spacings) and unstable region (at small orbital spacings) are separated by a transition zone, in which the (in)stability of the system is not guaranteed. We also devise a linear map to analyse the dynamical instability of the "planet + test-mass" system, and we find qualitatively similar results to the $N$-body simulations.


title: Probing finite-temperature observables in quantum simulators of spin systems with short-time dynamics

abstract: Preparing finite temperature states in quantum simulators of spin systems, such as trapped ions or Rydberg atoms in optical tweezers, is challenging due to their almost perfect isolation from the environment. Here, we show how finite-temperature observables can be obtained with an algorithm motivated from the Jarzynski equality and equivalent to the one in Lu, Banuls and Cirac, PRX Quantum 2, 020321 (2021). It consists of classical importance sampling of initial states and a measurement of the Loschmidt echo with a quantum simulator. We use the method as a quantum-inspired classical algorithm and simulate the protocol with matrix product states to analyze the requirements on a quantum simulator. This way, we show that a finite temperature phase transition in the long-range transverse field Ising model can be characterized in trapped ion quantum simulators. We propose a concrete measurement protocol for the Loschmidt echo and discuss the influence of measurement noise, dephasing, as well as state preparation and measurement errors. We argue that the algorithm is robust against those imperfections under realistic conditions.

categories: quant-ph cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el


title: A Note on the Analytic Structure of Celestial Amplitudes

abstract: Celestial amplitudes, obtained by applying Mellin transform and analytic continuation on "ordinary" amplitudes, have interesting properties which may provide useful insights on the underlying theory. Their analytic structures are thus of great interest and need to be better understood. In this paper, we critically examine the analytic structure of celestial amplitudes in a massless low-energy effective field theory. We find that, fixed-order loop contributions, which generate multipoles on the negative $\beta$-plane, in general do not provide an accurate description of the analytic structure of celestial amplitudes. By resumming over the leading logarithmic contributions using renormalization group equations (RGEs), we observe much richer analytic structures, which generally contain branch cuts. It is also possible to generate multipoles or shifted single poles if the RGEs satisfy certain relations. Including sub-leading logarithmic contributions is expected to introduce additional corrections to the picture. However, without a new approach, it is difficult to make a general statement since the analytic form of the Mellin transform is challenging to obtain.


title: The frustration-free fully packed loop model

abstract: We consider a quantum fully packed loop model on the square lattice with a frustration-free projector Hamiltonian and ring-exchange interactions acting on plaquettes. A boundary Hamiltonian is added to favour domain-wall boundary conditions and link ground state properties to the combinatorics and six-vertex model literature. We discuss how the boundary term fractures the Hilbert space into Krylov subspaces, and we prove that the Hamiltonian is ergodic within each subspace, leading to a series of energy-equidistant exact eigenstates in the lower end of the spectrum. Among them we systematically classify both finitely entangled eigenstates and product eigenstates. Using a recursion relation for enumerating half-plane configurations, we compute numerically the exact entanglement entropy of the ground state, confirming area law scaling. Finally, the spectrum is shown to be gapless in the thermodynamic limit with a trial state constructed by adding a twist to the ground state superposition.

categories: cond-mat.stat-mech math-ph math.MP quant-ph


title: Following the muon track of hierarchical sectors at LHCb

abstract: This paper reports a study of the experimental signatures of hierarchical sectors beyond the Standard Model characterized by a flavor-violating heavy vector and a set of light pseudo-Goldstone bosons $a_{1,2}$, spanning a large range of lifetimes. The non-minimal scalar spectrum triggers novel $B$ decays into multiple leptons that would have escaped the reach of current searches. Novel displaced vertex analyses at LHCb are therefore discussed to probe the hierarchical new physics, extending the scope of the tracking system of the detector. Additionally, the reach of the proposed CODEX-b experiment is studied as well. By exploiting the use of tracks only reconstructed in subsystems of the detector at the high-level LHCb trigger, ${\mathcal{B}(B_s^0\to a_1 a_2) < 10^{-8}}$ and ${\mathcal{B}(B^+\to K^+ a_1 a_2) < 10^{-9}}$ could be reached in the muon channel across seven orders of magnitude in the lifetime of the Goldstone bosons. Correspondingly, heavy-light particle couplings of order $\lesssim 1$ could be tested, potentially ruling out composite Higgs scenarios where the heavy and light sectors couple strongly.


title: General Relativistic Implicit Monte Carlo Radiation-Hydrodynamics

abstract: We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++: an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of the radiation-hydrodynamics equations, but generalized for both Newtonian and relativistic regimes. A multiple reference frame approach is used to geodesically transport photon packets (and solve the hydrodynamics equations) in the coordinate frame, while radiation-matter interactions are handled either in the fluid or electron frames then communicated via Lorentz boosts and orthonormal tetrad bases attached to the fluid. We describe a method for constructing estimators of radiation moments using path-weighting that generalizes to arbitrary coordinate systems in flat or curved spacetime. Absorption, emission, scattering, and relativistic Comptonization are among the matter interactions considered in this report. We discuss our formulations and numerical methods, and validate our models against a suite of radiation and coupled radiation-hydrodynamics test problems in both flat and curved spacetimes.


title: The Charged Higgs from the Bottom-Up: Probing Flavor at the LHC

abstract: We systematically study model-independent constraints on the three generic charged Higgs couplings to $b$-quarks and up-type quarks. While existing LHC searches have focussed on the $tb$ coupling, we emphasize that the LHC plays a crucial role in probing also $ub$ and $cb$ couplings, since constraints from flavor physics are weak. In particular we propose various new searches that can significantly extend the present reach on the parameter space by: i) looking for light charged Higgses that decay into $ub$-quarks, ii) probing charged Higgs couplings to light and top quarks using multi-$b$-jet signatures, iii) looking for single $b$-quarks in low-mass dijet searches, iv) searching for charge asymmetries induced by charged Higgs production via $ub$ couplings.


title: Quantum Interference of Hydrodynamic Modes in a Dirty Marginal Fermi Liquid

abstract: We study the electrical transport of a two-dimensional non-Fermi liquid with disorder, and we determine the first quantum correction to the semiclassical dc conductivity due to quantum interference. We consider a system with $N$ flavors of fermions coupled to SU($N$) critical matrix bosons. Motivated by the SYK model, we employ the bilocal field formalism and derive a set of finite-temperature saddle-point equations governing the fermionic and bosonic self-energies in the large-$N$ limit. Interestingly, disorder smearing induces a marginal Fermi liquid (MFL) self-energy for the fermions. We next consider fluctuations around the saddle points and derive a MFL-Finkel'stein nonlinear sigma model. We find that the Altshuler-Aronov quantum conductance correction gives linear-$T$ resistivity that can dominate over the Drude result at low temperature. The strong temperature dependence of the quantum correction arises due to rapid relaxation of the mediating quantum-critical bosons. We verify that our calculations explicitly satisfy the Ward identity at the semiclassical and quantum levels. Our results establish that quantum interference persists in two-particle hydrodynamic modes, even when quasiparticles are subject to strong (Planckian) dissipation.


title: Short GRB Host Galaxies I: Photometric and Spectroscopic Catalogs, Host Associations, and Galactocentric Offsets

abstract: We present a comprehensive optical and near-infrared census of the fields of 90 short gamma-ray bursts (GRBs) discovered in 2005-2021, constituting all short GRBs for which host galaxy associations are feasible ($\approx$ 60% of the total Swift short GRB population). We contribute 245 new multi-band imaging observations across 49 distinct GRBs and 25 spectra of their host galaxies. Supplemented by literature and archival survey data, the catalog contains 335 photometric and 40 spectroscopic data sets. The photometric catalog reaches $3\sigma$ depths of $\gtrsim 24-27$ mag and $\gtrsim 23-26$ mag for the optical and near-infrared bands, respectively. We identify host galaxies for 84 bursts, in which the most robust associations make up 54% (49/90) of events, while only a small fraction, 6.7%, have inconclusive host associations. Based on new spectroscopy, we determine 17 host spectroscopic redshifts with a range of $z\approx 0.15-1.6$ and find that $\approx$ 25-44% of Swift short GRBs originate from $z>1$. We also present the galactocentric offset catalog for 83 short GRBs. Taking into account the large range of individual measurement uncertainties, we find a median of projected offset of $\approx 7.9$ kpc, for which the bursts with the most robust associations have a smaller median of $\approx 4.9$ kpc. Our catalog captures more high-redshift and low-luminosity hosts, and more highly-offset bursts than previously found, thereby diversifying the population of known short GRB hosts and properties. In terms of locations and host luminosities, the populations of short GRBs with and without detectable extended emission are statistically indistinguishable. This suggests that they arise from the same progenitors, or from multiple progenitors which form and evolve in similar environments. All of the data products are available on the BRIGHT website.


title: Short GRB Host Galaxies. II. A Legacy Sample of Redshifts, Stellar Population Properties, and Implications for their Neutron Star Merger Origins

abstract: We present the stellar population properties of 69 short gamma-ray burst (GRB) host galaxies, representing the largest uniformly-modeled sample to-date. Using the Prospector stellar population inference code, we jointly fit photometry and/or spectroscopy of each host galaxy. We find a population median redshift of $z=0.64^{+0.83}_{-0.32}$ ($68\%$ confidence), including 10 new or revised photometric redshifts at $z\gtrsim1$. We further find a median mass-weighted age of $t_m=0.8^{+2.71}_{-0.53}$Gyr, stellar mass of $\log(M_*/M_\odot)=9.69^{+0.75}_{-0.65}$, star formation rate of SFR=$1.44^{+9.37}_{-1.35}M_\odot$yr$^{-1}$, stellar metallicity of $\log(Z_*/Z_\odot)=-0.38^{+0.44}_{-0.42}$, and dust attenuation of $A_V=0.43^{+0.85}_{-0.36}$~mag (68\% confidence). Overall, the majority of short GRB hosts are star-forming ($\approx84\%$), with small fractions that are either transitioning ($\approx6\%$) or quiescent ($\approx10\%$); however, we observe a much larger fraction ($\approx40\%$) of quiescent and transitioning hosts at $z\lesssim0.25$, commensurate with galaxy evolution. We find that short GRB hosts populate the star-forming main sequence of normal field galaxies, but do not include as many high-mass galaxies, implying that their binary neutron star (BNS) merger progenitors are dependent on a combination of host star formation and stellar mass. The distribution of ages and redshifts implies a broad delay-time distribution, with a fast-merging channel at $z>1$ and a decreased BNS formation efficiency at lower redshifts. If short GRB hosts are representative of BNS merger hosts within the horizon of current gravitational wave detectors, these results can inform future searches for electromagnetic counterparts. All of the data and modeling products are available on the BRIGHT website.


title: Anomalous Bootstrap on the half line

abstract: We study carefully the problem of the bootstrap on the half line. We show why one needs the full set of constraints derived from the Stieltjes theorem on the moment problem by reexamining previous results on the hydrogen atom. We also study the hydrogen atom at continuous angular momentum. We show that the constraints on the moment problem alone do not fix the boundary conditions in all cases and at least one of the positive matrices needs to be slightly enlarged to remove unphysical branches. We explain how to solve the more general problem of the bootstrap for Robin boundary conditions. The recursion relations that are usually used receive additional anomalous contributions. These corrections are necessary to compute the moments of the measure. We apply these to the linear potential and we show how the bootstrap matches the analytical results, based on the Airy function, for this example.


title: Advantages and limitations of quantum routing

abstract: The Swap gate is a ubiquitous tool for moving information on quantum hardware, yet it can be considered a classical operation because it does not entangle product states. Genuinely quantum operations could outperform Swap for the task of permuting qubits within an architecture, which we call routing. We consider quantum routing in two models: (1) allowing arbitrary two-qubit unitaries, or (2) allowing Hamiltonians with norm-bounded interactions. We lower bound the circuit depth or time of quantum routing in terms of spectral properties of graphs representing the architecture interaction constraints, and give a generalized upper bound for all simple connected $n$-vertex graphs. In particular, we give conditions for a superpolynomial classical-quantum routing separation, which exclude graphs with a small spectral gap and graphs of bounded degree. Finally, we provide examples of a quadratic separation between gate-based and Hamiltonian routing models with a constant number of local ancillas per qubit and of an $\Omega(n)$ speedup if we also allow fast local interactions.


title: [Re] Badder Seeds: Reproducing the Evaluation of Lexical Methods for Bias Measurement

abstract: Combating bias in NLP requires bias measurement. Bias measurement is almost always achieved by using lexicons of seed terms, i.e. sets of words specifying stereotypes or dimensions of interest. This reproducibility study focuses on the original authors' main claim that the rationale for the construction of these lexicons needs thorough checking before usage, as the seeds used for bias measurement can themselves exhibit biases. The study aims to evaluate the reproducibility of the quantitative and qualitative results presented in the paper and the conclusions drawn thereof. We reproduce most of the results supporting the original authors' general claim: seed sets often suffer from biases that affect their performance as a baseline for bias metrics. Generally, our results mirror the original paper's. They are slightly different on select occasions, but not in ways that undermine the paper's general intent to show the fragility of seed sets.


title: Introduction to special issue: Chaos Indicators, Phase Space and Chemical Reaction Dynamics

abstract: The study of the phase space of multidimensional systems is one of the central open problems in dynamical systems. Being able to distinguish chaoticity from regularity in nonlinear dynamical systems, as well as to determine the subspace of the phase space in which instabilities are expected to occur, is also an important field. To investigate these, diverse chaos indicators have been developed. The information provided from these indicators is important to understand the dynamical behavior of many systems in celestial mechanics and, more recently, in chemical reaction dynamics. This special issue contains new results around the topics: phase space structure, chaos indicators, and chemical reaction dynamics.


title: Free energy dissipation enhances spatial accuracy and robustness of Turing pattern in small reaction-diffusion systems

abstract: Accurate and robust spatial orders are ubiquitous in living systems. In 1952, Alan Turing proposed an elegant mechanism for pattern formation based on spontaneous breaking of the spatial translational symmetry in the underlying reaction-diffusion system. Much is understood about dynamics and structure of Turing patterns. However, little is known about the energetic cost of Turing pattern. Here, we study nonequilibrium thermodynamics of a small spatially extended biochemical reaction-diffusion system by using analytical and numerical methods. We find that the onset of Turing pattern requires a minimum energy dissipation to drive the nonequilibrium chemical reactions. Above onset, only a small fraction of the total energy expenditure is used to overcome diffusion for maintaining the spatial pattern. We show that the positioning error decreases as energy dissipation increases following the same tradeoff relationship between timing error and energy cost in biochemical oscillatory systems. In a finite system, we find that a specific Turing pattern exists only within a finite range of total molecule number, and energy dissipation broadens the range, which enhances the robustness of the Turing pattern against molecule number fluctuations in living cells. These results are verified in a realistic model of the Muk system underlying DNA segregation in E. coli, and testable predictions are made for the dependence of the accuracy and robustness of the spatial pattern on the ATP/ADP ratio. In general, the theoretical framework developed here can be applied to study nonequilibrium thermodynamics of spatially extended biochemical systems.


title: A Unified Spectroscopic and Photometric Model to Infer Surface Inhomogeneity: Application to Luhman 16B

abstract: Extremely large telescopes (ELTs) provide an opportunity to observe surface inhomogeneities for ultracool objects including M dwarfs, brown dwarfs (BDs), and gas giant planets via Doppler imaging and spectro-photometry techniques. These inhomogeneities can be caused by star spots, clouds, and vortices. Star spots and associated stellar flares play a significant role in habitability, either stifling life or catalyzing abiogenesis depending on the emission frequency, magnitude, and orientation. Clouds and vortices may be the source of spectral and photometric variability observed at the L/T transition of BDs and are expected in gas giant exoplanets. We develop a versatile analytical framework to model and infer surface inhomogeneities which can be applied to both spectroscopic and photometric data. This model is validated against a slew of numerical simulations. Using archival spectroscopic and photometric data, we infer star spot parameters (location, size, and contrast) and generate global surface maps for Luhman 16B (an early T dwarf and one of our solar system's nearest neighbors at a distance of approximately 2 pc). We confirm previous findings that Luhman 16B's atmosphere is inhomogeneous with time-varying features. In addition, we provide tentative evidence of longer timescale atmospheric structures such as dark equatorial and bright mid-latitude to polar spots. These findings are discussed in the context of atmospheric circulation and dynamics for ultracool dwarfs. Our analytical model will be valuable in assessing the feasibility of using ELTs to study surface inhomogeneities of gas giant exoplanets and other ultracool objects.


title: Traces on locally compact groups

abstract: We conduct a systematic study of traces on locally compact groups, in particular traces on their universal and reduced C*-algebras. We introduce the trace kernel, and examine its relation to the von Neumann kernel and to small-invariant neighbourhood (SIN) quotients. In doing so, we introduce the class of residually-$SIN$ groups, which contains both $SIN$ and maximally almost periodic groups. We examine in detail the trace kernel for connected groups. We study traces on reduced C*-algebras, giving a simple proof for compactly generated groups that existence of such a trace is equivalent to having an open normal amenable subgroup, and we display non-discrete groups admitting unique trace. We finish by examining amenable traces and the factorization property. We show for property (T) groups that amenable trace kernels coincide with von Neumann kernels. We show for totally disconnected groups that amenable trace separation implies the factorization property. We use amenable traces to give a simple proof that amenability of the group is equivalent to simultaneous nuclearity and possessing a trace of its reduced C*-algebra. As a final application of the results obtained in the paper, we address the embeddability of group C*-algebras into simple AF algebras. As a consequence, if a locally compact group is amenable and tracially separated (trace kernel is trivial), then its reduced C*-algebra is quasi-diagonal.


title: Radar Guided Dynamic Visual Attention for Resource-Efficient RGB Object Detection

abstract: An autonomous system's perception engine must provide an accurate understanding of the environment for it to make decisions. Deep learning based object detection networks experience degradation in the performance and robustness for small and far away objects due to a reduction in object's feature map as we move to higher layers of the network. In this work, we propose a novel radar-guided spatial attention for RGB images to improve the perception quality of autonomous vehicles operating in a dynamic environment. In particular, our method improves the perception of small and long range objects, which are often not detected by the object detectors in RGB mode. The proposed method consists of two RGB object detectors, namely the Primary detector and a lightweight Secondary detector. The primary detector takes a full RGB image and generates primary detections. Next, the radar proposal framework creates regions of interest (ROIs) for object proposals by projecting the radar point cloud onto the 2D RGB image. These ROIs are cropped and fed to the secondary detector to generate secondary detections which are then fused with the primary detections via non-maximum suppression. This method helps in recovering the small objects by preserving the object's spatial features through an increase in their receptive field. We evaluate our fusion method on the challenging nuScenes dataset and show that our fusion method with SSD-lite as primary and secondary detector improves the baseline primary yolov3 detector's recall by 14% while requiring three times fewer computational resources.


title: An Assessment Of Full-Wave Effects On Maxwellian Lower-Hybrid Wave Damping

abstract: Lower-hybrid current drive (LHCD) actuators are important components of modern day fusion experiments as well as proposed fusion reactors. However, simulations of LHCD often differ substantially from experimental results, and from each other, especially in the inferred power deposition profile shape. Here we investigate some possible causes of this discrepancy; "full-wave" effects such as interference and diffraction, which are omitted from standard raytracing simulations and the breakdown of the raytracing near reflections and caustics. We compare raytracing simulations to state-of-the-art full-wave simulations using matched hot-plasma dielectric tensors in realistic tokamak scenarios for the first time. We show that differences between full-wave simulations and raytracing in previous work were primarily due to numerical and physical inconsistencies in the simulations, and we demonstrate that good agreement between raytracing and converged full-wave simulations can be obtained in reactor relevant-scenarios with large ray caustics and in situations with weak damping.


title: Monkeypox Image Data collection

abstract: This paper explains the initial Monkeypox Open image data collection procedure. It was created by assembling images collected from websites, newspapers, and online portals and currently contains around 1905 images after data augmentation.


title: Seamless Interaction Design with Coexistence and Cooperation Modes for Robust Human-Robot Collaboration

abstract: A robot needs multiple interaction modes to robustly collaborate with a human in complicated industrial tasks. We develop a Coexistence-and-Cooperation (CoCo) human-robot collaboration system. Coexistence mode enables the robot to work with the human on different sub-tasks independently in a shared space. Cooperation mode enables the robot to follow human guidance and recover failures. A human intention tracking algorithm takes in both human and robot motion measurements as input and provides a switch on the interaction modes. We demonstrate the effectiveness of CoCo system in a use case analogous to a real world multi-step assembly task.


title: Properties of a Ternary Infinite Word

abstract: We study the properties of the ternary infinite word p = 012102101021012101021012 ... , that is, the fixed point of the map h:0->01, 1->21, 2->0. We determine its factor complexity, critical exponent, and prove that it is 2-balanced. We compute its abelian complexity and determine the lengths of its bispecial factors. Finally, we give a characterization of p in terms of avoided factors.


title: Real-Time Super-Resolution for Real-World Images on Mobile Devices

abstract: Image Super-Resolution (ISR), which aims at recovering High-Resolution (HR) images from the corresponding Low-Resolution (LR) counterparts. Although recent progress in ISR has been remarkable. However, they are way too computationally intensive to be deployed on edge devices, since most of the recent approaches are deep learning-based. Besides, these methods always fail in real-world scenes, since most of them adopt a simple fixed "ideal" bicubic downsampling kernel from high-quality images to construct LR/HR training pairs which may lose track of frequency-related details. In this work, an approach for real-time ISR on mobile devices is presented, which is able to deal with a wide range of degradations in real-world scenarios. Extensive experiments on traditional super-resolution datasets (Set5, Set14, BSD100, Urban100, Manga109, DIV2K) and real-world images with a variety of degradations demonstrate that our method outperforms the state-of-art methods, resulting in higher PSNR and SSIM, lower noise and better visual quality. Most importantly, our method achieves real-time performance on mobile or edge devices.


title: A probabilistic approach to vanishing viscosity for PDEs on the Wasserstein space

abstract: In this work we prove an analogue, for partial differential equations on the space of probability measures, of the classical vanishing viscosity result known for equations on the Euclidean space. Our result allows in particular to show that the value function arising in various problems of classical mechanics and games can be obtained as the limiting case of second order PDEs. The method of proof builds on stochastic analysis arguments and allows for instance to prove a Freindlin-Wentzell large deviation theorem for McKean-Vlasov equations.


title: Bayesian and Frequentist Inference for Synthetic Controls

abstract: The synthetic control method has become a widely popular tool to estimate causal effects with observational data. Despite this, inference for synthetic control methods remains challenging. Often, inferential results rely on linear factor model data generating processes. In this paper, we characterize the conditions on the factor model primitives (the factor loadings) for which the statistical risk minimizers are synthetic controls (in the simplex). Then, we propose a Bayesian alternative to the synthetic control method that preserves the main features of the standard method and provides a new way of doing valid inference. We explore a Bernstein-von Mises style result to link our Bayesian inference to the frequentist inference. For linear factor model frameworks we show that a maximum likelihood estimator (MLE) of the synthetic control weights can consistently estimate the predictive function of the potential outcomes for the treated unit and that our Bayes estimator is asymptotically close to the MLE in the total variation sense. Through simulations, we show that there is convergence between the Bayes and frequentist approach even in sparse settings. Finally, we apply the method to re-visit the study of the economic costs of the German re-unification and the Catalan secession movement. The Bayesian synthetic control method is available in the bsynth R-package.


title: An asynchronous object-oriented approach to the automation of the 0.8-meter George Mason University campus telescope in Python

abstract: We present a unique implementation of Python coding in an asynchronous object-oriented programming (OOP) framework to fully automate the process of collecting data with the George Mason University (GMU) Observatory's 0.8-meter telescope. The goal of this project is to perform automated follow-up observations for the Transiting Exoplanet Survey Satellite (TESS) mission, while still allowing for human control, monitoring, and adjustments. Prior to our implementation, the facility was computer-controlled by a human observer through a combination of webcams, TheSkyX, ASCOM Dome, MaxIm DL, and a weather station. We have automated slews and dome movements, CCD exposures, saving FITS images and metadata, initial focusing, guiding on the target, using the ambient temperature to adjust the focus as the telescope cools through the rest of the night, taking calibration images (darks and flats), and monitoring local weather data. The automated weather monitor periodically checks various weather data from multiple sources to automate the decision to close the observatory during adverse conditions. We have organized the OOP code structure in such a way that each hardware device or important higher-level process is categorized as its own object class or "module" with associated attributes and methods, with inherited common methods across modules for code reusability. To allow actions to be performed simultaneously across different modules, we implemented a multithreaded approach where each module is given its own CPU thread on which to operate concurrently with all other threads. After the initial few modules (camera, telescope, dome, data I/O) were developed, further development of the code was carried out in tandem with testing on sky on clear nights. The code, in its current state, has been tested and used for observations on 171 nights, with more planned usage and feature additions.


title: The Before, During, and After of Multi-Robot Deadlock

abstract: Collision avoidance for multirobot systems is a well-studied problem. Recently, control barrier functions (CBFs) have been proposed for synthesizing controllers that guarantee collision avoidance and goal stabilization for multiple robots. However, it has been noted that reactive control synthesis methods (such as CBFs) are prone to \textit{deadlock}, an equilibrium of system dynamics that causes the robots to stall before reaching their goals. In this paper, we analyze the closed-loop dynamics of robots using CBFs, to characterize controller parameters, initial conditions, and goal locations that invariably lead the system to deadlock. Using tools from duality theory, we derive geometric properties of robot configurations of an $N$ robot system once it is in deadlock and we justify them using the mechanics interpretation of KKT conditions. Our key deductions are that 1) system deadlock is characterized by a force-equilibrium on robots and 2) deadlock occurs to ensure safety when safety is on the brink of being violated. These deductions allow us to interpret deadlock as a subset of the state space, and we show that this set is non-empty and located on the boundary of the safe set. By exploiting these properties, we analyze the number of admissible robot configurations in deadlock and develop a provably-correct decentralized algorithm for deadlock resolution to safely deliver the robots to their goals. This algorithm is validated in simulations as well as experimentally on Khepera-IV robots.


title: Optimal Competitive-Ratio Control

abstract: Inspired by competitive policy designs approaches in online learning, new control paradigms such as competitive-ratio and regret-optimal control have been recently proposed as alternatives to the classical $\mathcal{H}_2$ and $\mathcal{H}_\infty$ approaches. These competitive metrics compare the control cost of the designed controller against the cost of a clairvoyant controller, which has access to past, present, and future disturbances in terms of ratio and difference, respectively. While prior work provided the optimal solution for the regret-optimal control problem, in competitive-ratio control, the solution is only provided for the sub-optimal problem. In this work, we derive the optimal solution to the competitive-ratio control problem. We show that the optimal competitive ratio formula can be computed as the maximal eigenvalue of a simple matrix, and provide a state-space controller that achieves the optimal competitive ratio. We conduct an extensive numerical study to verify this analytical solution, and demonstrate that the optimal competitive-ratio controller outperforms other controllers on several large scale practical systems. The key techniques that underpin our explicit solution is a reduction of the control problem to a Nehari problem, along with a novel factorization of the clairvoyant controller's cost. We reveal an interesting relation between the explicit solutions that now exist for both competitive control paradigms by formulating a regret-optimal control framework with weight functions that can also be utilized for practical purposes.


title: Odd-frequency pairing and time-reversal symmetry breaking for repulsive interactions

abstract: We study the pairing of fermions by an interaction consisting of a Hubbard repulsion, mimicking a screened Coulomb potential, and a dynamical phonon-mediated attraction. For such interaction, the gap equation allows even- and odd-frequency solutions $\Delta_e$ and $\Delta_o$. We show that odd-frequency pairing does not develop within the Eliashberg approximation due to over-critical pair-breaking from the self-energy. When vertex corrections are included, the pairing interaction gets stronger, and $\Delta_o$ can develop. We argue that even in this case keeping the self-energy is still a must as it cancels out the thermal piece in the gap equation. We further argue that $\Delta_o$ is not affected by Hubbard repulsion and for strong repulsion is comparable to a reduced $\Delta_e$. The resulting superconducting state is a superposition $\Delta_e \pm i \Delta_o$, which spontaneously breaks the time-reversal symmetry, despite that the pairing symmetry is an ordinary $s$-wave.


title: Onesweep: A Faster Least Significant Digit Radix Sort for GPUs

abstract: We present Onesweep, a least-significant digit (LSD) radix sorting algorithm for large GPU sorting problems residing in global memory. Our parallel algorithm employs a method of single-pass prefix sum that only requires ~2n global read/write operations for each digit-binning iteration. This exhibits a significant reduction in last-level memory traffic versus contemporary GPU radix sorting implementations, where each iteration of digit binning requires two passes through the dataset totaling ~3n global memory operations. On the NVIDIA A100 GPU, our approach achieves 29.4 GKey/s when sorting 256M random 32-bit keys. Compared to CUB, the current state-of-the-art GPU LSD radix sort, our approach provides a speedup of ~1.5x. For 32-bit keys with varied distributions, our approach provides more consistent performance compared to HRS, the current state-of-the-art GPU MSD radix sort, and outperforms it in almost all cases.


title: Seismology in the Solar System

abstract: The NASA InSight mission has helped to measure the deep interior of Mars using observations of seismic waves excited by marsquakes. Currently, installation of seismometers on the moon is foreseen. We review the case for seismic experiments on all major planetary bodies of the solar system. We discuss scientific goals in accordance with the Decadal survey for planetary science and astrobiology and the ESA Voyage 2050 program as well as technical challenges and potential mission concepts, to answer the question: Where could we do seismology on other planets and why should we do it?


title: Infinite Sumsets in Sets with Positive Density

abstract: Motivated by questions asked by Erdos, we prove that any set $A\subset{\mathbb N}$ with positive upper density contains, for any $k\in{\mathbb N}$, a sumset $B_1+\cdots+B_k$, where $B_1,\dots,B_k\subset{\mathbb N}$ are infinite. Our proof uses ergodic theory and relies on structural results for measure preserving systems. Our techniques are new, even for the previously known case of $k=2$.


title: Turbulent magnetic field in the HII region Sh 2-27

abstract: Magnetic fields in the turbulent interstellar medium (ISM) are a key element in understanding Galactic dynamics, but there are many observational challenges. One useful probe for studying the magnetic field component parallel to the line of sight (LoS) is Faraday rotation of linearly polarized radio synchrotron emission, combined with H$\alpha$ observations. HII regions are the perfect laboratories to probe such magnetic fields as they are localized in space, and are well-defined sources often with known distances and measurable electron densities. We chose the HII region Sharpless 2-27 (Sh 2-27). By using a map of the magnetic field strength along the LoS ($B_{\parallel}$) for the first time, we investigate the basic statistical properties of the turbulent magnetic field inside Sh 2-27. We study the scaling of the magnetic field fluctuations, compare it to the Kolmogorov scaling, and attempt to find an outer scale of the turbulent magnetic field fluctuations. We estimate the median value of $n_e$ as $7.3\pm0.1$ cm$^{-3}$, and the median value of $B_{\parallel}$ as $-4.5\pm0.1$ $\mu$G, which is comparable to the magnetic field strength in diffuse ISM. The slope of the structure function of the estimated $B_{\parallel}$-map is found to be slightly steeper than Kolmogorov, consistent with our Gaussian-random-field $B_{\parallel}$ simulations revealing that an input Kolmogorov slope in the magnetic field results in a somewhat steeper slope in $B_{\parallel}$. These results suggest that the lower limit to the outer scale of turbulence is 10 pc in the HII region, which is comparable to the size of the computation domain. This may indicate that the turbulence probed here could actually be cascading from the larger scales in the ambient medium, associated with the interstellar turbulence in the general ISM, which is illuminated by the presence of Sh 2-27.


title: Linear maps preserving products equal to primitive idempotents of an incidence algebra

abstract: Let $A$, $B$ be algebras and $a\in A$, $b\in B$ a fixed pair of elements. We say that a map $\varphi:A\to B$ preserves products equal to $a$ and $b$ if for all $a_1,a_2\in A$ the equality $a_1a_2=a$ implies $\varphi(a_1)\varphi(a_2)=b$. In this paper we study bijective linear maps $\varphi:I(X,F)\to I(X,F)$ preserving products equal to primitive idempotents of $I(X,F)$, where $I(X,F)$ is the incidence algebra of a finite connected poset $X$ over a field $F$. We fully characterize the situation, when such a map $\varphi$ exists, and whenever it does, $\varphi$ is either an automorphism of $I(X,F)$ or the negative of an automorphism of $I(X,F)$.


title: Continuum theories of structured dielectrics

abstract: Aequous dielectrics are ubiquitous in soft- and bio-nano matter systems. The theoretical description of such systems in terms of continuum (`macroscopic') theory remains a serious challenge. In this perspective we first review the existing continuum phenomenological approaches that have been developed in the past decades. In order to describe a path to advance continuum theory beyond these approaches we then take recourse to the Onsager-Dupuis theory of the dielectric behaviour of ice, which, for the case of a solid dielectric, exemplified important conceptual issues we deem relevant for the development of a more fundamental continuum theory of liquid dielectrics. Subsequently, we discuss our recently proposed continuum field theory of structured dielectrics, which provides a generalized approach to the dielectric behavior of such systems.


title: A Quantitative Simulation-based Modeling Approach for College Counseling Centers

abstract: College counseling centers in various universities have been tasked with the important responsibility of attending to the mental health needs of their students. Owing to the unprecedented recent surge of demand for such services, college counseling centers are facing several crippling resource-level challenges. This is leading to longer wait times which limits access to critical mental health services. To address these challenges, we construct a discrete-event simulation model that captures several intricate details of their operations and provides a data-driven framework to quantify the effect of different policy changes. In contrast to existing work on this matter, which are primarily based on qualitative assessments, the considered quantitative approach has the potential to lead to key observations that can assist counseling directors in constructing a system with desirable performance. To demonstrate the benefit of the considered simulation model, we use data specific to Texas A&M's Counseling & Psychological Services to run a series of numerical experiments. Our results demonstrate the predictive power of the simulation model, highlight a number of key observations, and identify policy changes that result in desirable system performance.


title: Seamless GPU acceleration for C++ based physics with the Metal Shading Language on Apple's M series unified chips

abstract: The M series of chips produced by Apple have proven a capable and power-efficient alternative to mainstream Intel and AMD x86 processors for everyday tasks. Additionally, the unified design integrating the central processing and graphics processing unit, have allowed these M series chips to excel at many tasks with heavy graphical requirements without the need for a discrete graphical processing unit (GPU), and in some cases even outperforming discrete GPUs. In this work, we show how the M series chips can be leveraged using the Metal Shading Language (MSL) to accelerate typical array operations in C++. More importantly, we show how the usage of MSL avoids the typical complexity of CUDA or OpenACC memory management, by allowing the central processing unit (CPU) and GPU to work in unified memory. We demonstrate how performant the M series chips are on standard one-dimensional and two-dimensional array operations such as array addition, SAXPY and finite difference stencils, with respect to serial and OpenMP accelerated CPU code. The reduced complexity of implementing MSL also allows us to accelerate an existing elastic wave equation solver (originally based on OpenMP accelerated C++) using MSL, with minimal effort, while retaining all CPU and OpenMP functionality. The resulting performance gain of simulating the wave equation is near an order of magnitude for specific settings. This gain attained from using MSL is similar to other GPU-accelerated wave-propagation codes with respect to their CPU variants, but does not come at much increased programming complexity that prohibits the typical scientific programmer to leverage these accelerators. This result shows how unified processing units can be a valuable tool to seismologists and computational scientists in general, lowering the bar to writing performant codes that leverage modern GPUs.


title: Gradient-preserving hyper-reduction of nonlinear dynamical systems via discrete empirical interpolation

abstract: This work proposes a hyper-reduction method for nonlinear parametric dynamical systems characterized by gradient fields such as Hamiltonian systems and gradient flows. The gradient structure is associated with conservation of invariants or with dissipation and hence plays a crucial role in the description of the physical properties of the system. Traditional hyper-reduction of nonlinear gradient fields yields efficient approximations that, however, lack the gradient structure. We focus on Hamiltonian gradients and we propose to first decompose the nonlinear part of the Hamiltonian, mapped into a suitable reduced space, into the sum of d terms, each characterized by a sparse dependence on the system state. Then, the hyper-reduced approximation is obtained via discrete empirical interpolation (DEIM) of the Jacobian of the derived d-valued nonlinear function. The resulting hyper-reduced model retains the gradient structure and its computationally complexity is independent of the size of the full model. Moreover, a priori error estimates show that the hyper-reduced model converges to the reduced model and the Hamiltonian is asymptotically preserved. Whenever the nonlinear Hamiltonian gradient is not globally reducible, i.e. its evolution requires high-dimensional DEIM approximation spaces, an adaptive strategy is performed. This consists in updating the hyper-reduced Hamiltonian via a low-rank correction of the DEIM basis. Numerical tests demonstrate the applicability of the proposed approach to general nonlinear operators and runtime speedups compared to the full and the reduced models.


title: R2U++: A Multiscale Recurrent Residual U-Net with Dense Skip Connections for Medical Image Segmentation

abstract: U-Net is a widely adopted neural network in the domain of medical image segmentation. Despite its quick embracement by the medical imaging community, its performance suffers on complicated datasets. The problem can be ascribed to its simple feature extracting blocks: encoder/decoder, and the semantic gap between encoder and decoder. Variants of U-Net (such as R2U-Net) have been proposed to address the problem of simple feature extracting blocks by making the network deeper, but it does not deal with the semantic gap problem. On the other hand, another variant UNET++ deals with the semantic gap problem by introducing dense skip connections but has simple feature extraction blocks. To overcome these issues, we propose a new U-Net based medical image segmentation architecture R2U++. In the proposed architecture, the adapted changes from vanilla U-Net are: (1) the plain convolutional backbone is replaced by a deeper recurrent residual convolution block. The increased field of view with these blocks aids in extracting crucial features for segmentation which is proven by improvement in the overall performance of the network. (2) The semantic gap between encoder and decoder is reduced by dense skip pathways. These pathways accumulate features coming from multiple scales and apply concatenation accordingly. The modified architecture has embedded multi-depth models, and an ensemble of outputs taken from varying depths improves the performance on foreground objects appearing at various scales in the images. The performance of R2U++ is evaluated on four distinct medical imaging modalities: electron microscopy (EM), X-rays, fundus, and computed tomography (CT). The average gain achieved in IoU score is 1.5+-0.37% and in dice score is 0.9+-0.33% over UNET++, whereas, 4.21+-2.72 in IoU and 3.47+-1.89 in dice score over R2U-Net across different medical imaging segmentation datasets.


title: Additive MIL: Intrinsically Interpretable Multiple Instance Learning for Pathology

abstract: Multiple Instance Learning (MIL) has been widely applied in pathology towards solving critical problems such as automating cancer diagnosis and grading, predicting patient prognosis, and therapy response. Deploying these models in a clinical setting requires careful inspection of these black boxes during development and deployment to identify failures and maintain physician trust. In this work, we propose a simple formulation of MIL models, which enables interpretability while maintaining similar predictive performance. Our Additive MIL models enable spatial credit assignment such that the contribution of each region in the image can be exactly computed and visualized. We show that our spatial credit assignment coincides with regions used by pathologists during diagnosis and improves upon classical attention heatmaps from attention MIL models. We show that any existing MIL model can be made additive with a simple change in function composition. We also show how these models can debug model failures, identify spurious features, and highlight class-wise regions of interest, enabling their use in high-stakes environments such as clinical decision-making.


title: Robust Topological Inference in the Presence of Outliers

abstract: The distance function to a compact set plays a crucial role in the paradigm of topological data analysis. In particular, the sublevel sets of the distance function are used in the computation of persistent homology -- a backbone of the topological data analysis pipeline. Despite its stability to perturbations in the Hausdorff distance, persistent homology is highly sensitive to outliers. In this work, we develop a framework of statistical inference for persistent homology in the presence of outliers. Drawing inspiration from recent developments in robust statistics, we propose a $\textit{median-of-means}$ variant of the distance function ($\textsf{MoM Dist}$), and establish its statistical properties. In particular, we show that, even in the presence of outliers, the sublevel filtrations and weighted filtrations induced by $\textsf{MoM Dist}$ are both consistent estimators of the true underlying population counterpart, and their rates of convergence in the bottleneck metric are controlled by the fraction of outliers in the data. Finally, we demonstrate the advantages of the proposed methodology through simulations and applications.

categories: math.ST cs.CG cs.LG math.AT stat.ML stat.TH


title: Accelerating the heat diffusion: fast thermal relaxation of a microcantilever

abstract: In most systems, thermal diffusion is intrinsically slow with respect to mechanical relaxation. We devise here a generic approach to accelerate the relaxation of the temperature field of a 1D object, in order to beat the mechanical time scales. This approach is applied to a micro-meter sized silicon cantilever, locally heated by a laser beam. A tailored driving protocol for the laser power is derived to quickly reach the thermal stationary state. The model is implemented experimentally yielding a significant acceleration of the thermal relaxation, up to a factor 30. An excellent agreement with the theoretical predictions is reported. This strategy allows a thermal steady state to be reached significantly faster than the natural mechanical relaxation.


title: Prophecy Variables for Hyperproperty Verification

abstract: Temporal logics for hyperproperties like HyperLTL use trace quantifiers to express properties that relate multiple system runs. In practice, the verification of such specifications is mostly limited to formulas without quantifier alternation, where verification can be reduced to checking a trace property over the self-composition of the system. Quantifier alternations like $\forall \pi. \exists \pi'. \phi$, can either be solved by complementation or with an interpretation as a two-person game between a $\forall$-player, who incrementally constructs the trace $\pi$, and an $\exists$-player, who constructs $\pi'$ in such a way that $\pi$ and $\pi'$ together satisfy $\phi$. The game-based approach is significantly cheaper but incomplete because the $\exists$-player does not know the future moves of the $\forall$-player. In this paper, we establish that the game-based approach can be made complete by adding ($\omega$-regular) temporal prophecies. Our proof is constructive, yielding an effective algorithm for the generation of a complete set of prophecies.


title: Dynamically constraining the length of the Milky Way bar

abstract: We present a novel method for constraining the length of the Galactic bar using 6D phase space information to directly integrate orbits. We define a pseudo-length for the Galactic bar, named $R_{Freq}$, based on the maximal extent of trapped bar orbits. We find the $R_{Freq}$ measured from orbits is consistent with the $R_{Freq}$ of the assumed potential only when the length of the bar and pattern speed of said potential is similar to the model from which the initial phase-space coordinates of the orbits are derived. Therefore, one can measure the model's or the Milky Way's bar length from 6D phase-space coordinates by determining which assumed potential leads to a self-consistent measured $R_{Freq}$. When we apply this method to $\approx$210,000 stars in APOGEE DR17 and $Gaia$ eDR3 data, we find a consistent result only for potential models with a dynamical bar length of $\approx$3.5 kpc. We find the Milky Way's trapped bar orbits extend out to only $\approx$3.5 kpc, but there is also an overdensity of stars at the end of the bar out to 4.8 kpc which could be related to an attached spiral arm. We also find that the measured orbital structure of the bar is strongly dependent on the properties of the assumed potential.


title: Two-Dimensional Spectroscopy of Two-Dimensional Materials

abstract: In this work we provide an exact and efficient numerical approach to simulate multi-time correlation functions in the Mahan-Nozi\`{e}res-De Dominicis model, which crudely mimics the spectral properties of doped two-dimensional semiconductors such as monolayer transition metal dichalcogenides. We apply this approach to study the coherent two-dimensional electronic spectra of the model. We show that several experimentally observed phenomena, such as peak asymmetry and coherent oscillations in the waiting-time dependence of the trion-exciton cross peaks of the two-dimensional rephasing spectrum, emerge naturally in our approach. Additional features are also present which find no correspondence with experimentally expected behavior. We trace these features to the infinite hole mass property of the model. We use this understanding to construct an efficient approach which filters out configurations associated with the lack of exciton recoil, enabling the connection to previous work and providing a route to the construction of realistic two-dimensional spectra over a broad doping range in two-dimensional semiconductors.


title: Optimal entanglement enhancing via conditional measurements

abstract: Enhancing quantum entanglement is important for many quantum information processing applications. In this paper, we consider a protocol for entanglement enhancing in a two-mode squeezed vacuum state (TMSVS), attained based on photon subtraction, photon catalysis, and photon addition. Central to such an operation is the task of mixing and detecting number states with each mode of TMSVS. We analyze various settings and find an optimal setup for improving the entanglement of the state.


title: Integral-direct Hartree-Fock and M{\o}ller-Plesset Perturbation Theory for Periodic Systems with Density Fitting: Application to the Benzene Crystal

abstract: We present an algorithm and implementation of integral-direct, density-fitted Hartree-Fock (HF) and second-order M{\o}ller-Plesset perturbation theory (MP2) for periodic systems. The new code eliminates the formerly prohibitive storage requirements and allows us to study systems one order of magnitude larger than before at the periodic MP2 level. We demonstrate the significance of the development by studying the benzene crystal in both the thermodynamic limit and the complete basis set limit, for which we predict an MP2 cohesive energy of $-72.8$ kJ/mol, which is about $10$--$15$ kJ/mol larger in magnitude than all previously reported MP2 calculations. Compared to the best theoretical estimate from literature, several modified MP2 models approach chemical accuracy in the predicted cohesive energy of the benzene crystal and hence may be promising cost-effective choices for future applications on molecular crystals.


title: Do-Operation Guided Causal Representation Learning with Reduced Supervision Strength

abstract: Causal representation learning has been proposed to encode relationships between factors presented in the high dimensional data. However, existing methods suffer from merely using a large amount of labeled data and ignore the fact that samples generated by the same causal mechanism follow the same causal relationships. In this paper, we seek to explore such information by leveraging do-operation to reduce supervision strength. We propose a framework that implements do-operation by swapping latent cause and effect factors encoded from a pair of inputs. Moreover, we also identify the inadequacy of existing causal representation metrics empirically and theoretically and introduce new metrics for better evaluation. Experiments conducted on both synthetic and real datasets demonstrate the superiorities of our method compared with state-of-the-art methods.


title: Spread of SARS-CoV-2 in a SIS model with vaccination and breakthrough infection

abstract: Although previous infection and vaccination provide protection against SARS-CoV-2 infection, both reinfection and breakthrough infection are possible events whose occurrence would increase with time after first exposure to the antigen and with the emergence of new variants of the virus. Periodic vaccination could counteract this decline in protection. In the present work, our aim was to develop and explore a model of SARS-CoV-2 spread with vaccination, reinfection and breakthrough infection. A modified deterministic SIS (Susceptible-Infected-Susceptible) model represented by a system of differential equations was designed. As in any SIS model, the population was divided into susceptible and infected individuals. But in our design, susceptible individuals were, in turn, grouped into three consecutive categories whose susceptibility increases with time after infection or vaccination. The model was studied by means of computer simulations, which were analysed qualitatively. The results obtained show that the prevalence, after oscillating between peaks and valleys, reaches a plateau phase. Moreover, as might be expected, the magnitude of the peaks and plateaus increases as the infection rate rises, the vaccination rate decreases and the rate of decay of protection conferred by vaccination or previous infection increases. Therefore, the present study suggests that, at least under certain conditions, the spread of SARS-CoV-2, although it could experience fluctuations, would finally evolve into an endemic form, with a more or less stable prevalence that would depend on the levels of infection and vaccination, and on the kinetics of post-infection and post-vaccination protection. However, it should be kept in mind that our development is a theoretical scheme with many limitations. For this reason, its predictions should be considered with great care.


title: A tensor factorization model of multilayer network interdependence

abstract: Multilayer networks describe the rich ways in which nodes are related by accounting for different relationships in separate layers. These multiple relationships are naturally represented by an adjacency tensor. In this work we study the use of the nonnegative Tucker decomposition (NNTuck) of such tensors under a KL loss as an expressive factor model that naturally generalizes existing stochastic block models of multilayer networks. Quantifying interdependencies between layers can identify redundancies in the structure of a network, indicate relationships between disparate layers, and potentially inform survey instruments for collecting social network data. We propose definitions of layer independence, dependence, and redundancy based on likelihood ratio tests between nested nonnegative Tucker decompositions. Using both synthetic and real-world data, we evaluate the use and interpretation of the NNTuck as a model of multilayer networks. Algorithmically, we show that using expectation maximization (EM) to maximize the log-likelihood under the NNTuck is step-by-step equivalent to tensorial multiplicative updates for the NNTuck under a KL loss, extending a previously known equivalence from nonnegative matrices to nonnegative tensors.