title: Early detection of knee osteoarthritis using deep learning on knee magnetic resonance images

abstract: The aim of this study was to investigate the influence of MRI and patient data on the prediction of knee osteoarthritis (OA) incidence using different deep learning architectures. Knee OA incidence within 24 months was predicted using the intermediate-weighted turbo spin-echo (IW-TSE) sequence of 593 patients from the Osteoarthritis Initiative. To extract a region of interest containing the knee joint from the IW-TSE sequence, a U-Net model was trained and used to segment bone on a dual echo steady state (DESS) sequence. Subsequently, IW-TSE and DESS sequences were registered and the DESS segmentations were transformed to the corresponding IW-TSE scans. The performance of MRI-based features in the prediction of knee OA incidence was tested using three different deep learning architectures: a residual network (ResNet), a densely connected convolutional network (DenseNet), and a convolutional variational autoencoder (CVAE). To evaluate the predictive performance of MRI-based features alone, the outputs of ResNet, DenseNet, and CVAE were coupled with patient data (i.e., age, gender, BMI) and used as input to a Logistic Regression (LR) Classifier. Knee OA was defined based on visual MRI and X-ray-based OA features. The ResNet and DenseNet showed similar results, with both methods having the area under the receiver operating characteristic curve (AUC) values up to 0.6269. The best performing OA detection model was CVAE with an AUC of 0.6699 when combined with patient data and an AUC of 0.6689 when used alone as input to the LR classifier. The results showed that three deep learning algorithms have similar metrics when using IW-TSE MRIs and their performance increased with the inclusion of patient data, which shows the strong influence of variables such as age, gender, and BMI on the detection of knee OA.


title: Oxygen dissociation on the C3N monolayer: A first-principles study

abstract: The oxygen dissociation and the oxidized structure on the pristine C3N monolayer in exposure to air are the inevitably critical issues for the C3N engineering and surface functionalization yet have not been revealed in detail. Using the first-principles calculations, we have systematically investigated the possible O2 adsorption sites, various O2 dissociation pathways and the oxidized structures. It is demonstrated that the pristine C3N monolayer shows more O2 physisorption sites and exhibits stronger O2 adsorption than the pristine graphene. Among various dissociation pathways, the most preferable one is a two-step process involving an intermediate state with the chemisorbed O2 and the barrier is lower than that on the pristine graphene, indicating that the pristine C3N monolayer is more susceptible to oxidation than the pristine graphene. Furthermore, we found that the most stable oxidized structure is not produced by the most preferable dissociation pathway but generated from a direct dissociation process. These results can be generalized into a wide range of temperatures and pressures using ab initio atomistic thermodynamics. Our findings deepen the understanding of the chemical stability of 2D crystalline carbon nitrides under ambient conditions, and could provide insights into the tailoring of the surface chemical structures via doping and oxidation.


title: CLONeR: Camera-Lidar Fusion for Occupancy Grid-aided Neural Representations

abstract: Recent advances in neural radiance fields (NeRFs) achieve state-of-the-art novel view synthesis and facilitate dense estimation of scene properties. However, NeRFs often fail for large, unbounded scenes that are captured under very sparse views with the scene content concentrated far away from the camera, as is typical for field robotics applications. In particular, NeRF-style algorithms perform poorly: (1) when there are insufficient views with little pose diversity, (2) when scenes contain saturation and shadows, and (3) when finely sampling large unbounded scenes with fine structures becomes computationally intensive. This paper proposes CLONeR, which significantly improves upon NeRF by allowing it to model large outdoor driving scenes that are observed from sparse input sensor views. This is achieved by decoupling occupancy and color learning within the NeRF framework into separate Multi-Layer Perceptrons (MLPs) trained using LiDAR and camera data, respectively. In addition, this paper proposes a novel method to build differentiable 3D Occupancy Grid Maps (OGM) alongside the NeRF model, and leverage this occupancy grid for improved sampling of points along a ray for volumetric rendering in metric space. Through extensive quantitative and qualitative experiments on scenes from the KITTI dataset, this paper demonstrates that the proposed method outperforms state-of-the-art NeRF models on both novel view synthesis and dense depth prediction tasks when trained on sparse input data.


title: Decohering Tensor Network Quantum Machine Learning Models

abstract: Tensor network quantum machine learning (QML) models are promising applications on near-term quantum hardware. While decoherence of qubits is expected to decrease the performance of QML models, it is unclear to what extent the diminished performance can be compensated for by adding ancillas to the models and accordingly increasing the virtual bond dimension of the models. We investigate here the competition between decoherence and adding ancillas on the classification performance of two models, with an analysis of the decoherence effect from the perspective of regression. We present numerical evidence that the fully-decohered unitary tree tensor network (TTN) with two ancillas performs at least as well as the non-decohered unitary TTN, suggesting that it is beneficial to add at least two ancillas to the unitary TTN regardless of the amount of decoherence may be consequently introduced.


title: A Simple Method for Determining Center Shift and Spectral Area of a M\"ossbauer Spectrum

abstract: A very simple method for determining the center (isomer) shift, CS, and the spectral area, A, of a M\"ossbauer spectrum is outlined. Its applicability is demonstrated in two examples viz. pyrite and a ternary sigma-phase Fe-Cr-Ni compound. Sets of the spectra recorded in the temperature interval of 78-295 K for the former, and 78-295 K for the latter were analyzed with the simple and a standard method. CS(T)- and A(T)-data obtained with both ways of the spectra analysis were analyzed in terms of the Debye model. The determined values of the Debye temperatures agree well within the error limit with each other proving that this very simple method is suitable for a quick determining of CS and A spectral parameters.


title: Models of $VTC^0$ as exponential integer parts

abstract: We prove that (additive) ordered group reducts of nonstandard models of the bounded arithmetical theory $\mathsf{VTC^0}$ are recursively saturated in a rich language with predicates expressing the integers, rationals, and logarithmically bounded numbers. Combined with our previous results on the construction of the real exponential function on completions of models of $\mathsf{VTC^0}$, we show that every countable model of $\mathsf{VTC^0}$ is an exponential integer part of a real-closed exponential field.


title: Estimation of Correlation Matrices from Limited time series Data using Machine Learning

abstract: Correlation matrices contain a wide variety of spatio-temporal information about a dynamical system. Predicting correlation matrices from partial time series information of a few nodes characterizes the spatio-temporal dynamics of the entire underlying system. This information can help to predict the underlying network structure, e.g., inferring neuronal connections from spiking data, deducing causal dependencies between genes from expression data, and discovering long spatial range influences in climate variations. Traditional methods of predicting correlation matrices utilize time series data of all the nodes of the underlying networks. Here, we use a supervised machine learning technique to predict the correlation matrix of entire systems from finite time series information of a few randomly selected nodes. The accuracy of the prediction validates that only a limited time series of a subset of the entire system is enough to make good correlation matrix predictions. Furthermore, using an unsupervised learning algorithm, we furnish insights into the success of the predictions from our model. Finally, we employ the machine learning model developed here to real-world data sets.


title: Revisiting Outer Optimization in Adversarial Training

abstract: Despite the fundamental distinction between adversarial and natural training (AT and NT), AT methods generally adopt momentum SGD (MSGD) for the outer optimization. This paper aims to analyze this choice by investigating the overlooked role of outer optimization in AT. Our exploratory evaluations reveal that AT induces higher gradient norm and variance compared to NT. This phenomenon hinders the outer optimization in AT since the convergence rate of MSGD is highly dependent on the variance of the gradients. To this end, we propose an optimization method called ENGM which regularizes the contribution of each input example to the average mini-batch gradients. We prove that the convergence rate of ENGM is independent of the variance of the gradients, and thus, it is suitable for AT. We introduce a trick to reduce the computational cost of ENGM using empirical observations on the correlation between the norm of gradients w.r.t. the network parameters and input examples. Our extensive evaluations and ablation studies on CIFAR-10, CIFAR-100, and TinyImageNet demonstrate that ENGM and its variants consistently improve the performance of a wide range of AT methods. Furthermore, ENGM alleviates major shortcomings of AT including robust overfitting and high sensitivity to hyperparameter settings.


title: $h \to \Upsilon \gamma$ Decay: Smoking Gun Signature of Wrong-Sign $hb\bar{b}$ Coupling

abstract: We analyze in a model-independent way the potential to probe new physics using Higgs decay to $\Upsilon \gamma$. The $h \to \Upsilon \gamma$ decay width is unusually small in the Standard Model because of an accidental cancellation among the direct and indirect decay diagrams. Thus, any new physics that can modify the direct or the indirect decay amplitudes disrupts the accidental Standard Model cancellation and can potentially lead to a relatively large decay width for $h \to \Upsilon \gamma$. Here, we carry out a detailed model-independent analysis of the possible new physics that can disrupt this cancellation. We demonstrate that after taking into account all possible constraints on Higgs production and decay processes from experimental measurements, the wrong-sign $h b \bar{b}$ coupling is the only scenario in which the $h \to \Upsilon \gamma$ decay width can be changed by almost two orders of magnitude. We conclude that an observation of a significantly enhanced $h \to \Upsilon \gamma$ decay width at the LHC or any future collider will be a conclusive evidence of a wrong-sign $h b \bar{b}$ coupling.


title: Recent results on homological properties of binomial edge ideal of graphs

abstract: In this article, we give a comprehensive survey of the recent progress of research on binomial edge ideal of a graph since 2018.


title: Radiation forces and confinement of neutral particles into the pulse envelope. New regime of collision ionization

abstract: The aim of this work is to find analytical expressions of the radiation force and potential densities arising from a laser pulse propagating in a dielectric media. In the femtosecond region the force becomes strong enough to confine neutral particles into the pulse envelope and translate them with the group velocity in gases. Additionally, if the trapped particles are with high density, the probability of collision with the free atoms and molecules in air become significant. The collision energies are in the range of 12 - 24 GeV and high enough to ionize the neutral atoms. Thus, a new type of collision ionization can be observed, when powerful femtosecond pulses propagate in gaseous media.


title: Role of electronic thermal transport in amorphous metal recrystallization: a molecular dynamics study

abstract: Recrystallization of glasses is important in a wide range of applications including electronics and reactive materials. Molecular dynamics (MD) has been used to provide an atomic picture of this process, but prior work has neglected the thermal transport role of electrons, the dominant thermal carrier in metallic systems. We characterize the role of electronic thermal conductivity on the velocity of recrystallization in Ni using MD coupled to a continuum description of electronic thermal transport via a two-temperature model. Our simulations show that for strong enough coupling between electrons and ions, the increased thermal conductivity removes the heat from the exothermic recrystallization process more efficiently, leading to a lower effective temperature at the recrystallization front and, consequently, lower propagation velocity. We characterize how electron-phonon coupling strength and system size affects front propagation velocity. Interestingly, we find that initial recrystallization velocity increases with decreasing in system size due to higher overall temperatures. Overall, we show that a more accurate description of thermal transport due to the incorporation of electrons results in better agreement with experiments.


title: Spinning correlators in $\mathcal{N} = 2$ SCFTs: Superspace and AdS amplitudes

abstract: We study four-point functions of spinning operators in the flavor current multiplet in four dimensional $\mathcal{N}=2$ SCFTs, using superspace techniques. In particular we explicitly construct the differential operators relating the different components of the super-correlator. As a byproduct of our analysis, we report the computation of the four-point amplitudes of gluons in bosonic Yang-Mills theories on $\mathrm{AdS}_5$ and we give evidence of an AdS double copy relation between the gluon amplitude and its gravitational counterpart.


title: Hierarchical Relational Learning for Few-Shot Knowledge Graph Completion

abstract: Knowledge graphs (KGs) are known for their large scale and knowledge inference ability, but are also notorious for the incompleteness associated with them. Due to the long-tail distribution of the relations in KGs, few-shot KG completion has been proposed as a solution to alleviate incompleteness and expand the coverage of KGs. It aims to make predictions for triplets involving novel relations when only a few training triplets are provided as reference. Previous methods have mostly focused on designing local neighbor aggregators to learn entity-level information and/or imposing sequential dependency assumption at the triplet level to learn meta relation information. However, valuable pairwise triplet-level interactions and context-level relational information have been largely overlooked for learning meta representations of few-shot relations. In this paper, we propose a hierarchical relational learning method (HiRe) for few-shot KG completion. By jointly capturing three levels of relational information (entity-level, triplet-level and context-level), HiRe can effectively learn and refine the meta representation of few-shot relations, and consequently generalize very well to new unseen relations. Extensive experiments on two benchmark datasets validate the superiority of HiRe against other state-of-the-art methods.


title: Transformers in Remote Sensing: A Survey

abstract: Deep learning-based algorithms have seen a massive popularity in different areas of remote sensing image analysis over the past decade. Recently, transformers-based architectures, originally introduced in natural language processing, have pervaded computer vision field where the self-attention mechanism has been utilized as a replacement to the popular convolution operator for capturing long-range dependencies. Inspired by recent advances in computer vision, remote sensing community has also witnessed an increased exploration of vision transformers for a diverse set of tasks. Although a number of surveys have focused on transformers in computer vision in general, to the best of our knowledge we are the first to present a systematic review of recent advances based on transformers in remote sensing. Our survey covers more than 60 recent transformers-based methods for different remote sensing problems in sub-areas of remote sensing: very high-resolution (VHR), hyperspectral (HSI) and synthetic aperture radar (SAR) imagery. We conclude the survey by discussing different challenges and open issues of transformers in remote sensing. Additionally, we intend to frequently update and maintain the latest transformers in remote sensing papers with their respective code at: https://github.com/VIROBO-15/Transformer-in-Remote-Sensing


title: Co-Imitation: Learning Design and Behaviour by Imitation

abstract: The co-adaptation of robots has been a long-standing research endeavour with the goal of adapting both body and behaviour of a system for a given task, inspired by the natural evolution of animals. Co-adaptation has the potential to eliminate costly manual hardware engineering as well as improve the performance of systems. The standard approach to co-adaptation is to use a reward function for optimizing behaviour and morphology. However, defining and constructing such reward functions is notoriously difficult and often a significant engineering effort. This paper introduces a new viewpoint on the co-adaptation problem, which we call co-imitation: finding a morphology and a policy that allow an imitator to closely match the behaviour of a demonstrator. To this end we propose a co-imitation methodology for adapting behaviour and morphology by matching state distributions of the demonstrator. Specifically, we focus on the challenging scenario with mismatched state- and action-spaces between both agents. We find that co-imitation increases behaviour similarity across a variety of tasks and settings, and demonstrate co-imitation by transferring human walking, jogging and kicking skills onto a simulated humanoid.


title: Linear response quantum transport through interacting multi-orbital nanostructures

abstract: Nanoelectronics devices, such as quantum dot systems or single-molecule transistors, consist of a quantum nanostructure coupled to a macroscopic external electronic circuit. Thermoelectric transport between source and drain leads is controlled by the quantum dynamics of the lead-coupled nanostructure, through which a current must pass. Strong electron interactions due to quantum confinement on the nanostructure produce nontrivial conductance signatures such as Coulomb blockade and Kondo effects, which become especially pronounced at low temperatures. In this work we first provide a modern review of standard quantum transport techniques, focusing on the linear response regime, and highlight the strengths and limitations of each. In the second part, we develop an improved numerical scheme for calculation of the ac linear electrical conductance through generic interacting nanostructures, based on the numerical renormalization group (NRG) method, and explicitly demonstrate its utility in terms of accuracy and efficiency. In the third part we derive low-energy effective models valid in various commonly-encountered situations, and from them we obtain simple analytical expressions for the low-temperature conductance. This indirect route via effective models, although approximate, allows certain limitations of conventional methodologies to be overcome, and provides physical insights into transport mechanisms. Finally, we apply and compare the various techniques, taking the two-terminal triple quantum dot and the serial multi-level double dot devices as nontrivial benchmark systems.


title: Strategies for quantum-optimized construction of interpolating operators in classical simulations of lattice quantum field theories

abstract: It has recently been argued that noisy intermediate-scale quantum computers may be used to optimize interpolating operator constructions for lattice quantum field theory (LQFT) calculations on classical computers. Here, two concrete realizations of the method are developed and implemented. The first approach is to maximize the overlap, or fidelity, of the state created by an interpolating operator acting on the vacuum state to the target eigenstate. The second is to instead minimize the energy expectation value of the interpolated state. These approaches are implemented in a proof-of-concept calculation in (1+1)-dimensions for a single-flavor massive Schwinger model to obtain quantum-optimized interpolating operator constructions for a vector meson state in the theory. Although fidelity maximization is preferable in the absence of noise due to quantum gate errors, it is found that energy minimization is more robust to these effects in the proof-of-concept calculation. This work serves as a concrete demonstration of how quantum computers in the intermediate term might be used to accelerate classical LQFT calculations.


title: GOALS-JWST: Resolving the Circumnuclear Gas Dynamics in NGC 7469 in the Mid-Infrared

abstract: The nearby, luminous infrared galaxy (LIRG) NGC 7469 hosts a Seyfert nucleus with a circumnuclear star-forming ring and is thus the ideal local laboratory for investigating the starburst--AGN connection in detail. We present integral-field observations of the central 1.3 kpc region in NGC 7469 obtained with the JWST Mid-InfraRed Instrument. Molecular and ionized gas distributions and kinematics at a resolution of {\sim}100 pc over the 4.9 - 7.6{\mu}m region are examined to study gas dynamics influenced by the central AGN. The low-ionization [Fe II] {\lambda}5.34{\mu}m and [Ar II] {\lambda}6.99{\mu}m lines are bright on the nucleus and in the starburst ring, as opposed to H2 S(5) {\lambda}6.91{\mu}m which is strongly peaked at the center and surrounding ISM. The high-ionization [Mg V] line is resolved and shows a broad, blueshifted component associated with the outflow. It has a nearly face-on geometry that is strongly peaked on the nucleus, where it reaches a maximum velocity of -650 km/s, and extends about 400 pc to the East. Regions of enhanced velocity dispersion in H2 and [Fe II] {\sim}180 pc from the AGN that also show high L(H2)/L(PAH) and L([Fe II])/L(Pf{\alpha}) ratios to the W and N of the nucleus pinpoint regions where the ionized outflow is depositing energy, via shocks, into the dense interstellar medium between the nucleus and the starburst ring. These resolved mid-infrared observations of the nuclear gas dynamics demonstrate the power of JWST and its high-sensitivity integral-field spectroscopic capability to resolve feedback processes around supermassive black holes in the dusty cores of nearby LIRGs.


title: Cross-Camera Deep Colorization

abstract: In this paper, we consider the color-plus-mono dual-camera system and propose an end-to-end convolutional neural network to align and fuse images from it in an efficient and cost-effective way. Our method takes cross-domain and cross-scale images as input, and consequently synthesizes HR colorization results to facilitate the trade-off between spatial-temporal resolution and color depth in the single-camera imaging system. In contrast to the previous colorization methods, ours can adapt to color and monochrome cameras with distinctive spatial-temporal resolutions, rendering the flexibility and robustness in practical applications. The key ingredient of our method is a cross-camera alignment module that generates multi-scale correspondences for cross-domain image alignment. Through extensive experiments on various datasets and multiple settings, we validate the flexibility and effectiveness of our approach. Remarkably, our method consistently achieves substantial improvements, i.e., around 10dB PSNR gain, upon the state-of-the-art methods. Code is at: https://github.com/IndigoPurple/CCDC


title: AutoPET Challenge 2022: Automatic Segmentation of Whole-body Tumor Lesion Based on Deep Learning and FDG PET/CT

abstract: Automatic segmentation of tumor lesions is a critical initial processing step for quantitative PET/CT analysis. However, numerous tumor lesion with different shapes, sizes, and uptake intensity may be distributed in different anatomical contexts throughout the body, and there is also significant uptake in healthy organs. Therefore, building a systemic PET/CT tumor lesion segmentation model is a challenging task. In this paper, we propose a novel training strategy to build deep learning models capable of systemic tumor segmentation. Our method is validated on the training set of the AutoPET 2022 Challenge. We achieved 0.7574 Dice score, 0.0299 false positive volume and 0.2538 false negative volume on preliminary test set.The code of our work is available on the following link: https://github.com/ZZZsn/MICCAI2022-autopet.


title: The current orbit of Atlas (SXV)

abstract: With the success of the Cassini-Huygens mission, the dynamic complexity surrounding natural satellites of Saturn began to be elucidated. New ephemeris could be calculated with a higher level of precision, which made it possible to study in detail the resonant phenomena and, in particular, the 54:53 near mean-motion resonance between Prometheus and Atlas. For this task, we have mapped in details the domains of the resonance with dense sets of initial conditions and distinct ranges of parameters. Our initial goal was to identify possible regions in the phase space of Atlas for which some critical angles, associated with the 54:53 mean motion have a stable libration. Our investigations revealed that there is no possibility for the current Atlas orbital configuration to have any regular behavior since it is in a chaotic region located at the boundary of the 54:53 mean-motion resonance phase space. This result is in accordance with previous works (Cooper et al. 2015, Renner et al. 2016). In this work, we generalize such investigations by showing detailed aspects of the Atlas-Prometheus 54:53 mean-motion resonance, like the extension of the chaotic layers, the thin domain of the center of the 54:53 resonance, the proximity of other neighborhood resonances, among other secondary conclusions. In particular, we have also shown that even in the deep interior of the resonance, it is difficult to map periodic motion of the resonant pair for very long time spans.


title: General theory of emergent elasticity for second-order topological phase transitions

abstract: The raise of the symmetry breaking mechanism by Landau[1] is a landmark in the studies of phase transitions. The Kosterlitz-Thouless phase transition[2-3] and the fractional quantum Hall effect[4], however, are believed to be induced by another mechanism: topology change. Despite rapid development of the theory of topological orders[5-7], a unified theoretical framework describing this new paradigm of phase transition and its relation to the Landau paradigm, is not seen. Here, we establish such a framework based on variational principle, and show that the critical condition for any second-order topological phase transitions is loss of positive-definiteness of a topologically protected second-order variation of the free energy. A topologically protected variation of a field solution is performed with respect to its emergent displacements, which is introduced by constructing an emergent elasticity problem for the solution. The Landau paradigm of phase transitions studies global property changes induced instability, while topological phase transitions study local property changes induced topological instability. The general effectiveness of this criterion is shown through analyzing the topological stability of several prototype solutions in both real space (two-kink solutions of the sine-Gordon model; an isolated skyrmion in chiral magnets) and reciprocal space (phonon spectrum of a monoatomic chain and a diatomic chain; band structure of a monoatomic chain within the Kronig-Penny model). Every field is emergent elastic with spatially modulated stiffness, and changes of topological property occur at its softened points. We anticipate our work to be a starting point for a general study of topological phase transitions in all field theories.


title: Exploiting Fairness to Enhance Sensitive Attributes Reconstruction

abstract: In recent years, a growing body of work has emerged on how to learn machine learning models under fairness constraints, often expressed with respect to some sensitive attributes. In this work, we consider the setting in which an adversary has black-box access to a target model and show that information about this model's fairness can be exploited by the adversary to enhance his reconstruction of the sensitive attributes of the training data. More precisely, we propose a generic reconstruction correction method, which takes as input an initial guess made by the adversary and corrects it to comply with some user-defined constraints (such as the fairness information) while minimizing the changes in the adversary's guess. The proposed method is agnostic to the type of target model, the fairness-aware learning method as well as the auxiliary knowledge of the adversary. To assess the applicability of our approach, we have conducted a thorough experimental evaluation on two state-of-the-art fair learning methods, using four different fairness metrics with a wide range of tolerances and with three datasets of diverse sizes and sensitive attributes. The experimental results demonstrate the effectiveness of the proposed approach to improve the reconstruction of the sensitive attributes of the training set.


title: Flexible transition probability model for assessing cost-effectiveness of breast cancer screening

abstract: Breast cancer is the most common cancer among Western women. Fortunately, organized screening has reduced breast cancer mortality and, consequently, the European Union has recommended screening with mammography for 50-69-year-old women. This recommendation is followed well in Europe. Widening the screening target age further is supported by conditional recommendations for 45-49- and 70-74-year-old women. However, before extending screening to new age groups, it's essential to carefully consider the benefits and costs locally as circumstances vary between different regions and/or countries. We propose a new approach to assess cost-effectiveness of breast cancer screening for a long-ongoing program with incomplete historical screening data. The new model is called flexible stage distribution model. It is based on estimating the stage distributions of breast cancer cases under different screening strategies. In the model, an ongoing screening strategy may be used as a baseline and other screening strategies may be incorporated by changes in the incidence rates. The model is flexible, as it enables to apply different approaches for estimating the altered stage distributions. Thus, if randomized data is available, one may rely on that. On the other hand, if randomized data is not available, altered stage distributions may be estimated by extrapolating the stage distributions of the youngest and oldest screened/non-screened age groups. We apply the proposed flexible stage distribution model for assessing incremental cost of extending the current biennial breast cancer screening to younger and older target ages in Finland.


title: A Method for Discovering Novel Classes in Tabular Data

abstract: In Novel Class Discovery (NCD), the goal is to find new classes in an unlabeled set given a labeled set of known but different classes. While NCD has recently gained attention from the community, no framework has yet been proposed for heterogeneous tabular data, despite being a very common representation of data. In this paper, we propose TabularNCD, a new method for discovering novel classes in tabular data. We show a way to extract knowledge from already known classes to guide the discovery process of novel classes in the context of tabular data which contains heterogeneous variables. A part of this process is done by a new method for defining pseudo labels, and we follow recent findings in Multi-Task Learning to optimize a joint objective function. Our method demonstrates that NCD is not only applicable to images but also to heterogeneous tabular data. Extensive experiments are conducted to evaluate our method and demonstrate its effectiveness against 3 competitors on 7 diverse public classification datasets.


title: Connections and loops intertwinning

abstract: On a finite graph, we prove that trace of holonomies determine an intertwining relation between merge-and-split generators on collections of geodesic loops ensembles and Casimir operators on unitary connections. By adding a deformation part to the generator on loops, this result is extended to the Casimir operator modified in order to be self adjoint with respect to Yang- Mills measure.


title: A Framework for Extracting and Encoding Features from Object-Centric Event Data

abstract: Traditional process mining techniques take event data as input where each event is associated with exactly one object. An object represents the instantiation of a process. Object-centric event data contain events associated with multiple objects expressing the interaction of multiple processes. As traditional process mining techniques assume events associated with exactly one object, these techniques cannot be applied to object-centric event data. To use traditional process mining techniques, the object-centric event data are flattened by removing all object references but one. The flattening process is lossy, leading to inaccurate features extracted from flattened data. Furthermore, the graph-like structure of object-centric event data is lost when flattening. In this paper, we introduce a general framework for extracting and encoding features from object-centric event data. We calculate features natively on the object-centric event data, leading to accurate measures. Furthermore, we provide three encodings for these features: tabular, sequential, and graph-based. While tabular and sequential encodings have been heavily used in process mining, the graph-based encoding is a new technique preserving the structure of the object-centric event data. We provide six use cases: a visualization and a prediction use case for each of the three encodings. We use explainable AI in the prediction use cases to show the utility of both the object-centric features and the structure of the sequential and graph-based encoding for a predictive model.


title: TESS search for substellar companions through pulsation timing of $\delta$ Scuti stars. I. Discovery of companions around Chang 134 and V393 Car

abstract: Early-type main-sequence pulsating stars such as \ds{} variables are one of the least explored class of targets in the search for exoplanets. Pulsation timing (PT) is an alternative technique to the most effective search methods. It exploits the light-travel-time effect (LTTE) to infer the presence of additional massive bodies around a pulsating star by measuring a periodic phase modulation of its signal. PT has been extremely successful in discovering and characterizing stellar binaries when it was applied to high-precision light curves over long temporal baselines, such as those delivered by the Kepler mission. In favorable conditions, the sensitivity of PT can reach the planetary-mass regime; one such candidate has already been claimed. The advent of TESS, with its nearly full-sky coverage and the availability of full-frame images, opens a great opportunity to expand this field of research. In this work, we present a pilot study aimed to understand the potential of PT as applied to TESS data, which are considerably different from Kepler data in terms of photometric noise, sampling cadence, and temporal baseline. We focused on the most favorable class of \ds{}, that is, those showing strong pulsations and very simple frequency spectra. After the development of a customized pipeline, we were able to detect candidate companions for two targets in the (sub-)stellar mass regime: Chang~134 ($43\pm 5$~$M_\mathrm{jup}$, $P\simeq 82$~d) and V393 Car ($\gtrsim 100$~$M_\mathrm{jup}$, $P\gtrsim 700$~d). Our results also highlight the limiting factors of this technique and the importance of an accurate absolute time calibration for future missions such as PLATO.


title: Compact Binary Foreground Subtraction in Next-Generation Ground-Based Observatories

abstract: The stochastic gravitational-wave backgrounds (SGWBs) for current detectors are dominated by binary black-hole (BBH) and binary neutron-star (BNS) coalescences. The sensitivity of current networks of gravitational-wave (GW) detectors allows only a small fraction of BBHs and BNSs to be resolved and subtracted, but previous work indicated that the situation should significantly improve with next-generation (XG) observatories. We revisit these conclusions by taking into account waveform-modeling uncertainties, updated astrophysical models, and (crucially) the full set of parameters that must be estimated to remove the resolved sources. Compared to previous studies, we find that the residual background from BBHs and BNSs is large even with XG detector networks. New data analysis methods will thus be required to observe the SGWB from cosmic supernovae or contributions from early-Universe phenomena like cosmic strings, stiff post-inflation fluids, or axion inflation.


title: Accessing CKM suppressed top decays at the LHC

abstract: We propose an extension of the existing experimental strategy for measuring branching fractions of top quark decays, targeting specifically $t\to j_q W$, where $j_q$ is a light quark jet. The improved strategy uses orthogonal $b$- and $q$-taggers, and adds a new observable, the number of light-quark-tagged jets, to the already commonly used observable, the fraction of $b$-tagged jets in an event. Careful inclusion of the additional complementary observable significantly increases the expected statistical power of the analysis, with the possibility of excluding $|V_{tb}|=1$ at $95\%$ C.L. at the HL-LHC, and accessing directly the standard model value of $|V_{td}|^2+|V_{ts}|^2$.


title: Breaking the intensity-bias degeneracy in line intensity mapping

abstract: Line intensity mapping is rapidly gaining prominence as a tool to understand galaxy evolution and cosmology at high redshift. However, the standard power spectrum analysis for these surveys suffers from a tight degeneracy between the overall mean line intensity and the bias of the emitting galaxies. We outline a new formalism for the Voxel Intensity Distribution (VID), a one-point statistic which is known to contain information beyond the power spectrum. We use this new calculation to show for the first time that the VID can be used to break this key degeneracy in intensity mapping data.


title: A minimal model for adaptive SIS epidemics

abstract: The interplay between disease spreading and personal risk perception is of key importance for modelling the spread of infectious diseases. We propose a planar system of ordinary differential equations (ODEs) to describe the co-evolution of a spreading phenomenon and the average link density in the personal contact network. Contrary to standard epidemic models,we assume that the contact network changes based on the current prevalence of the disease in the population, i.e.\ it adapts to the current state of the epidemic. We assume that personal risk perception is described using two functional responses: one for link-breaking and one for link-creation. The focus is on applying the model to epidemics, but we highlight other possible fields of application. We derive an explicit form for the basic reproduction number and guarantee the existence of at least one endemic equilibrium, for all possible functional responses. Moreover, we show that for all functional responses, limit cycles do not exist.


title: Suppressing variance in 21-cm signal simulations during reionization

abstract: Current best limits on the 21-cm signal during reionization are provided at large scales ($\gtrsim$100 Mpc). To model these scales, enormous simulation volumes are required which are computationally expensive. We find that the primary source of uncertainty at these large scales is sample variance, which decides the minimum size of simulations required to analyse current and upcoming observations. In large-scale structure simulations, the method of `fixing' the initial conditions (ICs) to exactly follow the initial power spectrum and `pairing' two simulations with exactly out-of-phase ICs has been shown to significantly reduce sample variance. Here we apply this `fixing and pairing' (F\&P) approach to reionization simulations whose clustering signal originates from both density fluctuations and reionization bubbles. Using a semi-numerical code, we show that with the traditional method, simulation boxes of $L\simeq 500$ (300) Mpc are required to model the large-scale clustering signal at $k$=0.1 Mpc$^{-1}$ with a precision of 5 (10) per cent. Using F\&P, the simulation boxes can be reduced by a factor of 2 to obtain the same precision level. We conclude that the computing costs can be reduced by at least a factor of 4 when using the F\&P approach.


title: Nanoscale crack propagation in clay with water adsorption through reactive MD modeling

abstract: The atomic-scale cracking mechanism in clay is vital in discovering the cracking mechanism of clay at the continuum scale in that clay is a nanomaterial. In this article, we investigate mechanisms of mode I and mode II crack propagations in pyrophyllite and Ca-montmorillonite with water adsorption through reactive molecular dynamics with a bond-order force field. Clay water adsorption is considered by adding water molecules to the clay surface. During the equilibration stage, water adsorption could cause bending deformation of the pre-defined edge crack region. The relatively small orientating angle of water molecules indicates the formation of hydrogen bonds in the crack propagation process. The peak number density of adsorbed water decreases with the increasing strains. The atomistic structure evolution of the crack tip under loading is analyzed to interpret the nanoscale crack propagation mechanism. The numerical results show that the crack tip first gets blunted with a significant increase in the radius of the curvature of the crack tip and a slight change in crack length. The crack tip blunting is studied by tracking the crack tip opening distance and O-Si-O angle in the tetrahedral Si-O cell in mode I and II cracks. We compare bond-breaking behaviors between Al-O and Si-O. It is found that Si-O bond breaking is primarily responsible for crack propagation. The critical stress intensity factor and critical energy release rate are determined from MD simulation results.


title: Experience With Machine Protection Systems At PIP2IT

abstract: The PIP-II Injector Test 1 PIP2IT facility accelerator was assembled in multiple stages in 2014 $-$ 2021 to test concepts and components of the future PIPII linac that is being constructed at Fermilab. In its final configuration, PIP2IT accelerated a 0.55 ms x 20 Hz x 2 mA H$-$ beam to 16 MeV. To protect elements of the beam line, a Machine Protection System MPS was implemented and commissioned. The beam was interrupted faster than 10$\mu$s when excessive beam loss was detected. The paper describes the MPS architecture, methods of the loss detection, procedure of the beam interruption, and operational experience at PIP2IT.


title: The Cosmic Ultraviolet Baryon Survey (CUBS) V: On the Thermodynamic Properties of the Cool Circumgalactic Medium at $z < 1$

abstract: This paper presents a systematic study of the photoionization and thermodynamic properties of the cool circumgalactic medium (CGM) as traced by rest-frame ultraviolet absorption lines around 26 galaxies at redshift $z\lesssim1$. The study utilizes both high-quality far-ultraviolet and optical spectra of background QSOs and deep galaxy redshift surveys to characterize the gas density, temperature, and pressure of individual absorbing components and to resolve their internal non-thermal motions. The derived gas density spans more than three decades, from $\log (n_{\rm H}/{\rm cm^{-3}}) \approx -4$ to $-1$, while the temperature of the gas is confined in a narrow range of $\log (T/{\rm K})\approx 4.3\pm 0.3$. In addition, a weak anti-correlation between gas density and temperature is observed, consistent with the expectation of the gas being in photoionization equilibrium. Furthermore, decomposing the observed line widths into thermal and non-thermal contributions reveals that more than 30% of the components at $z\lesssim 1$ exhibit line widths driven by non-thermal motions, in comparison to $<20$% found at $z\approx 2$-3. Attributing the observed non-thermal line widths to intra-clump turbulence, we find that massive quenched galaxies on average exhibit higher non-thermal broadening/turbulent energy in their CGM compared to star-forming galaxies at $z\lesssim 1$. Finally, strong absorption features from multiple ions covering a wide range of ionization energy (e.g., from Mg II to O IV) can be present simultaneously in a single absorption system with kinematically aligned component structure, but the inferred pressure in different phases may differ by a factor of $\approx 10$.


title: Cubic-Regularized Newton for Spectral Constrained Matrix Optimization and its Application to Fairness

abstract: Matrix functions are utilized to rewrite smooth spectral constrained matrix optimization problems as smooth unconstrained problems over the set of symmetric matrices which are then solved via the cubic-regularized Newton method. A second-order chain rule identity for matrix functions is proven to compute the higher-order derivatives to implement cubic-regularized Newton, and a new convergence analysis is provided for cubic-regularized Newton for matrix vector spaces. We demonstrate the applicability of our approach by conducting numerical experiments on both synthetic and real datasets. In our experiments, we formulate a new model for estimating fair and robust covariance matrices in the spirit of the Tyler's M-estimator (TME) model and demonstrate its advantage.


title: Efficient Adiabatic Preparation of Tensor Network States

abstract: We propose and study a specific adiabatic path to prepare those tensor network states that are unique ground states of few-body parent Hamiltonians in finite lattices, which include normal tensor network states, as well as other relevant nonnormal states. This path guarantees a gap for finite systems and allows for efficient numerical simulation. In one dimension, we numerically investigate the preparation of a family of states with varying correlation lengths and the one-dimensional Affleck-Kennedy-Lieb-Tasaki (AKLT) state and show that adiabatic preparation can be much faster than standard methods based on sequential preparation. We also apply the method to the two-dimensional AKLT state on the hexagonal lattice, for which no method based on sequential preparation is known, and show that it can be prepared very efficiently for relatively large lattices.


title: How Descriptive are GMRES Convergence Bounds?

abstract: GMRES is a popular Krylov subspace method for solving linear systems of equations involving a general non-Hermitian coefficient matrix. The conventional bounds on GMRES convergence involve polynomial approximation problems in the complex plane. Three popular approaches pose this approximation problem on the spectrum, the field of values, or pseudospectra of the coefficient matrix. We analyze and compare these bounds, illustrating with six examples the success and failure of each. When the matrix departs from normality due only to a low-dimensional invariant subspace, we discuss how these bounds can be adapted to exploit this structure. Since the Arnoldi process that underpins GMRES provides approximations to the pseudospectra, one can estimate the GMRES convergence bounds as an iteration proceeds.


title: Elaboration-Generating Commonsense Question Answering at Scale

abstract: In question answering requiring common sense, language models (e.g., GPT-3) have been used to generate text expressing background knowledge that helps improve performance. Yet the cost of working with such models is very high; in this work, we finetune smaller language models to generate useful intermediate context, referred to here as elaborations. Our framework alternates between updating two language models -- an elaboration generator and an answer predictor -- allowing each to influence the other. Using less than 0.5% of the parameters of GPT-3, our model outperforms alternatives with similar sizes and closes the gap on GPT-3 on four commonsense question answering benchmarks. Human evaluations show that the quality of the generated elaborations is high.


title: Gravitational form factors of the delta resonance in chiral EFT

abstract: The leading one-loop corrections to the gravitational form factors of the delta resonance are calculated in the framework of chiral effective field theory. Various contributions to the energy-momentum tensor and the renormalization of the low-energy constants are worked out. Using the small scale expansion, expressions for static quantities are obtained and the real and imaginary parts of the gravitational form factors are calculated numerically.


title: Asymptotically Hyperbolic Einstein Constraint Equations with Apparent Horizon Boundary and the Penrose Inequality for Perturbations of Schwarzschild-AdS

abstract: We prove the existence of asymptotically hyperbolic solutions to the vacuum Einstein constraint equations with a marginally outer trapped boundary of positive mean curvature, using the constant mean curvature conformal method. As an application of this result, we verify the Penrose inequality for certain perturbations of Schwarzschild Anti-de Sitter black hole initial data.


title: Smiles in Profiles: Improving Fairness and Efficiency Using Estimates of User Preferences in Online Marketplaces

abstract: Online platforms often face challenges being both fair (i.e., non-discriminatory) and efficient (i.e., maximizing revenue). Using computer vision algorithms and observational data from a micro-lending marketplace, we find that choices made by borrowers creating online profiles impact both of these objectives. We further support this conclusion with a web-based randomized survey experiment. In the experiment, we create profile images using Generative Adversarial Networks that differ in a specific feature and estimate its impact on lender demand. We then counterfactually evaluate alternative platform policies and identify particular approaches to influencing the changeable profile photo features that can ameliorate the fairness-efficiency tension.


title: Resumming Quark's Longitudinal Momentum Logarithms in LaMET Expansion of Lattice PDFs

abstract: In the large-momentum expansion for parton distribution functions (PDFs), the natural physics scale is the longitudinal momentum ($p_z$) of the quarks (or gluons) in a large-momentum hadron. We show how to expose this scale dependence through resumming logarithms of the type $\ln^n p_z/\mu$ in the matching coefficient, where $\mu$ is a fixed renormalization scale. The result enhances the accuracy of the expansion at moderate $p_z>1$ GeV, and at the same time, clearly shows that the partons cannot be approximated from quarks with $p_z\sim \Lambda_{\rm QCD}$ which are not predominantly collinear with the parent hadron momentum, consistent with power counting of the large-momentum effective theory. The same physics mechanism constrains the coordinate space expansion at large distances $z$, the conjugate of $p_z$, as illustrated in the example of fitting the moments of the PDFs.


title: A critical analysis of `Relative facts do not exist. Relational quantum mechanics is incompatible with quantum mechanics' by Jay Lawrence, Marcin Markiewicz and Marek \'{Z}ukowski

abstract: We discuss a recent work by J.~Lawrence et al.[arxiv.org/abs/2208.11793] criticizing relational quantum mechanics (RQM) and based on a famous nonlocality theorem Going back to Greenberger Horne and Zeilinger (GHZ). Here, we show that the claims presented in this recent work are unjustified and we debunk the analysis.


title: Noise In Intense Electron Bunches

abstract: We report on our investigations into density fluctuations in electron bunches. Noise and density fluctuations in relativistic electron bunches, accelerated in a linac, are of critical importance to various Coherent Electron Cooling (CEC) [1-5] concepts as well as to free-electron lasers (FELs). For CEC, the beam noise results in additional diffusion that counteracts cooling. In SASE FELs, a micro-wave instability starts from the initial noise in the beam and eventually leads to the beam microbunching yielding coherent radiation, and the initial noise in the FEL bandwidth plays a useful role. In seeded FELs, in contrast, such noise interferes with the seed signal, so that reducing noise at the initial seed wavelength would lower the seed laser power requirement [6-8]. Status of the project will be presented.


title: The static energy of a quark-antiquark pair from Laplacian eigenmodes

abstract: We test a method for computing the static quark-antiquark potential in lattice QCD, which is not based on Wilson loops, but where the trial states are formed by eigenvector components of the covariant lattice Laplace operator. The runtime of this method is significantly smaller than the standard Wilson loop calculation, when computing the static potential not only for on-axis, but also for many off-axis quark-antiquark separations, i.e., when a fine spatial resolution is required. We further improve the signal by using multiple eigenvector pairs, weighted with Gaussian profile functions of the eigenvalues, providing a basis for a generalized eigenvalue problem (GEVP), as it was recently introduced to improve distillation in meson spectroscopy. We show results with the new method for the static potential with dynamical fermions and demonstrate its efficiency compared to traditional Wilson loop calculations. The method presented here can also be applied to compute hybrid or tetra-quark potentials and to static-light systems.


title: Constraining the origin of giant exoplanets via elemental abundance measurements

abstract: The origin of close-in giant planets is a key open question in planet formation theory. The two leading models are (i) formation at the outer disk followed by migration and (ii) in situ formation. In this work we determine the atmospheric composition of warm Jupiters for both formation scenarios. We perform N-body simulations of planetesimal accretion interior and exterior to the water ice-line for various planetary formation locations, planetary masses, and planetesimal sizes to estimate the accreted heavy-element mass and final planetary composition. We find that the two models differ significantly: migrating giant planets have 2-14 times higher metallicities than planets that form in situ. The ratio between refractories and volatiles is found to be above one for migrating planets but below 0.4 for planets that form in situ. We also identify very different trends between heavy-element enrichment and planetary mass for these two formation mechanisms. While the metallicity of migrating planets is found to increase with decreasing planetary mass, it is about constant for in situ formation. Our study highlights the importance of measuring the atmospheric composition of warm Jupiters and its connection to their formation and evolutionary paths.


title: Weighted $p(\cdot)$-Poincar\'e and Sobolev inequality]{Weighted $% p(\cdot )$-Poincar\'{e} and Sobolev inequalities for vector fields satisfiying H\"{o}rmander's condition and applications

abstract: In this paper we will establish different weighted Poincar\'{e} inequalities with variable exponents on Carnot-Carath\'{e}odory spaces or Carnot groups. We will use different techniques to obtain these inequalities. For vector fields satisfying H\"{o}rmander's condition in variable non-isotropic Sobolev spaces, we consider a weight in the variable Muckenhoupt class $% A_{p(\cdot ),p^{\ast }(\cdot )}$, where the exponent $p(\cdot )$ satisfies appropriate hypotheses, and in this case we obtain the first order weighted Poincar\'{e} inequalities with variable exponents. In the case of Carnot groups we also set up the higher order weighted Poincar\'{e} inequalities with variable exponents. For these results the crucial part is proving the boundedness of the fractional integral operator on Lebesgue spaces with weighted and variable exponents on spaces of homogeneous type. Moreover, using other techniques, we extend some of these results when the exponent satisfies a jump condition and the weight is in a smaller Muckenhoupt class. Finally, we will use these weighted Poincar\'{e} inequalities to establish the existence and uniqueness of a minimizer to the Dirichlet energy integral for a problem involving a degenerate $p(\cdot )$-Laplacian with zero boundary values in Carnot groups.


title: Better Peer Grading through Bayesian Inference

abstract: Peer grading systems aggregate noisy reports from multiple students to approximate a true grade as closely as possible. Most current systems either take the mean or median of reported grades; others aim to estimate students' grading accuracy under a probabilistic model. This paper extends the state of the art in the latter approach in three key ways: (1) recognizing that students can behave strategically (e.g., reporting grades close to the class average without doing the work); (2) appropriately handling censored data that arises from discrete-valued grading rubrics; and (3) using mixed integer programming to improve the interpretability of the grades assigned to students. We show how to make Bayesian inference practical in this model and evaluate our approach on both synthetic and real-world data obtained by using our implemented system in four large classes. These extensive experiments show that grade aggregation using our model accurately estimates true grades, students' likelihood of submitting uninformative grades, and the variation in their inherent grading error; we also characterize our models' robustness.


title: Local vanishing mean oscillation

abstract: We consider various notions of vanishing mean oscillation on a (possibly unbounded) domain $\Omega \subset \mathbb{R}^n$, and prove an analogue of Sarason's theorem, giving sufficient conditions for the density of bounded Lipschitz functions in the nonhomogeneous space $\rm{vmo}(\Omega)$. We also study $\rm{cmo}(\Omega)$, the closure in $\rm{bmo}(\Omega)$ of the continuous functions with compact support in $\Omega$. Using these approximation results, we prove that there is a bounded extension from $\rm{vmo}(\Omega)$ and $\rm{cmo}(\Omega)$ to the corresponding spaces on $\mathbb{R}^n$, if and only if $\Omega$ is a locally uniform domain.


title: FuzzerAid: Grouping Fuzzed Crashes Based On Fault Signatures

abstract: Fuzzing has been an important approach for finding bugs and vulnerabilities in programs. Many fuzzers deployed in industry run daily and can generate an overwhelming number of crashes. Diagnosing such crashes can be very challenging and time-consuming. Existing fuzzers typically employ heuristics such as code coverage or call stack hashes to weed out duplicate reporting of bugs. While these heuristics are cheap, they are often imprecise and end up still reporting many "unique" crashes corresponding to the same bug. In this paper, we present FuzzerAid that uses fault signatures to group crashes reported by the fuzzers. Fault signature is a small executable program and consists of a selection of necessary statements from the original program that can reproduce a bug. In our approach, we first generate a fault signature using a given crash. We then execute the fault signature with other crash inducing inputs. If the failure is reproduced, we classify the crashes into the group labeled with the fault signature; if not, we generate a new fault signature. After all the crash inducing inputs are classified, we further merge the fault signatures of the same root cause into a group. We implemented our approach in a tool called FuzzerAid and evaluated it on 3020 crashes generated from 15 real-world bugs and 4 large open source projects. Our evaluation shows that we are able to correctly group 99.1% of the crashes and reported only 17 (+2) "unique" bugs, outperforming the state-of-the-art fuzzers.


title: Modeling frequency shifts of collective bubble resonances with the boundary element method

abstract: Increasing the number of closely-packed air bubbles immersed in water changes the frequency of the Minnaert resonance. The collective interactions between bubbles in a small ensemble are primarily in the same phase, causing them to radiate a spherically-symmetric field that peaks at a frequency lower than the Minnaert resonance for a single bubble. In contrast, large periodic arrays include bubbles that are further apart than half the wavelength, so that collective resonances have bubbles oscillating in opposite phases, ultimately creating a fundamental resonance at a frequency higher than the single-bubble Minnaert resonance. This work investigates the transition in resonance behavior using a modal analysis of a mass-spring system and a boundary element method. We significantly reduce the computational complexity of the full-wave solver to a linear dependence on the number of bubbles in a rectangular array. The simulated acoustic fields confirm the initial downshift in resonance frequency and the strong influence of collective resonances when the array has hundreds of bubbles covering more than half the wavelength. These results are essential in understanding the low-frequency resonance characteristics of bubble ensembles, which have important applications in diverse fields such as underwater acoustics, quantum physics, and metamaterial design.


title: Spectral Asymptotics at Thresholds for a Dirac-type Operator on $\mathbb{Z}^2$

abstract: In this article, we provide the spectral analysis of a Dirac-type operator on $\mathbb{Z}^2$ by describing the behavior of the spectral shift function associated with a sign-definite trace-class perturbation by a multiplication operator. We prove that it remains bounded outside a single threshold and obtain its main asymptotic term in the unbounded case. Interestingly, we show that the constant in the main asymptotic term encodes the interaction between a flat band and whole non-constant bands. The strategy used is the reduction of the spectral shift function to the eigenvalue counting function of some compact operator which can be studied as a toroidal pseudo-differential operator.


title: Population level information combined parameter estimation from complex survey datasets

abstract: We consider an empirical likelihood framework for inference for a statistical model based on an informative sampling design and population-level information. The population-level information is summarized in the form of estimating equations and incorporated into the inference through additional constraints. Covariate information is incorporated both through the weights and the estimating equations. The estimator is based on conditional weights. We show that under usual conditions, with population size increasing unbounded, the estimates are strongly consistent, asymptotically unbiased, and normally distributed. Moreover, they are more efficient than other probability-weighted analogs. Our framework provides additional justification for inverse probability weighted score estimators in terms of conditional empirical likelihood. We give an application to demographic hazard modeling by combining birth registration data with panel survey data to estimate annual first birth probabilities.


title: Supersymmetric Quantum Potentials Analogs of Classical Electrostatic Fields

abstract: A relation between classical electrostatic fields and Schr\"odinger-like Hamiltonians is evidenced. Hence, supersymmetric quantum potentials analogous to classical electrostatic fields can be constructed. Proposing an ansatz for the electrostatic potential as the natural logarithm of a nodeless function, it is demonstrated that the electrostatic fields fulfil the Bernoulli equation associated to a second-order confluent supersymmetric transformation. By using the so-called confluent algorithm, it is possible, given a charge density, to find the corresponding electrostatic field as well as the supersymmetric potentials. Furthermore, the associated charge density and the electrostatic field profile of Schr\"odinger-like solvable potentials can be determined.

categories: math-ph math.MP physics.class-ph quant-ph


title: Ultramassive black holes formed by triple quasar mergers at $z\sim 2$

abstract: The origin of rare and elusive ultramassive black holes (UMBH, with MBH > 1e10 Msun) is an open question. Using the large volume cosmological hydrodynamic simulation ASTRID, we report on the formation of an extremely massive UMBH with MBH ~ 1e11 Msun at z~2. The UMBH is assembled as a result of two successive mergers of massive galaxies each with stellar mass M* > 3e11 Msun that also produces a bright, rare triple quasar system powered by three ~10^9 Msun black holes. The second merger of supermassive black holes (SMBHs) follows the first after 150 Myrs. The merger events lead to sustained Eddington accretion onto the central SMBH, forming an UMBH in the center of a massive compact stellar core with M* > 2e12 Msun. The strong feedback of the UMBH quenches the surrounding star formation to < 10 Msun/yr in the inner 50 kpc/h region. There are two more UMBHs with MBH > 5e10 Msun at z>2 in ASTRID which are also produced by major mergers of galaxies, and their progenitors can be observed as quasar triplets of lower luminosity. The rarely observed quasar multiples can be the cradle of UMBHs at high redshift, and likely end up in the center of the most massive clusters.


title: Improving Contextual Recognition of Rare Words with an Alternate Spelling Prediction Model

abstract: Contextual ASR, which takes a list of bias terms as input along with audio, has drawn recent interest as ASR use becomes more widespread. We are releasing contextual biasing lists to accompany the Earnings21 dataset, creating a public benchmark for this task. We present baseline results on this benchmark using a pretrained end-to-end ASR model from the WeNet toolkit. We show results for shallow fusion contextual biasing applied to two different decoding algorithms. Our baseline results confirm observations that end-to-end models struggle in particular with words that are rarely or never seen during training, and that existing shallow fusion techniques do not adequately address this problem. We propose an alternate spelling prediction model that improves recall of rare words by 34.7% relative and of out-of-vocabulary words by 97.2% relative, compared to contextual biasing without alternate spellings. This model is conceptually similar to ones used in prior work, but is simpler to implement as it does not rely on either a pronunciation dictionary or an existing text-to-speech system.


title: Recent progress with observations and models to characterize the magnetic fields from star-forming cores to protostellar disks

abstract: In this review paper, we aim at providing a global outlook on the progresses made in the recent years to characterize the role of magnetic fields during the embedded phases of the star formation process. Thanks to the development of observational capabilities and the parallel progress in numerical models capturing most of the important physics at work during star formation, it has recently become possible to confront detailed predictions of magnetized models to observational properties of the youngest protostars. We provide an overview of the most important consequences when adding magnetic fields to state-of-the-art models of protostellar formation, emphasizing their role to shape the resulting star(s) and their disk(s). We discuss the importance of magnetic field coupling to set the efficiency of magnetic processes, and provide a review of observational works putting constraints on the two main agents responsible for the coupling in star-forming cores: dust grains and ionized gas. We recall the physical processes and observational methods allowing to trace the magnetic field topology and its intensity in embedded protostars, and review the main steps, success and limitations in comparing real observations to synthetic observations from the non-ideal MHD models. Finally, we discuss the main threads of observational evidence that suggest a key role of magnetic fields for star and disk formation, and propose a scenario solving the angular momentum for star formation, also highlighting the remaining tensions that exist between models and observations.


title: Mixed Reality for Mechanical Design and Assembly Planning

abstract: Design for Manufacturing and Assembly (DFMA) is a crucial design stage within the heavy vehicle manufacturing process that involves optimising the order and feasibility of the parts assembly process to reduce manufacturing complexity and overall cost. Existing work has focused on conducting DFMA within virtual environments to reduce manufacturing costs, but users are less able to relate and compare physical characteristics of a virtual component with real physical objects. Therefore, a Mixed Reality (MR) application is developed for engineers to visualise and manipulate assembly parts virtually, conduct and plan out an assembly within its intended physical environment. Two pilot evaluations were conducted with both engineering professionals and non-engineers to assess effectiveness of the software for assembly planning. Usability results suggest that the application is overall usable (M=56.1, SD=7.89), and participants felt a sense of involvement in the activity (M=13.1, SD=3.3). Engineering professionals see the application as a useful and cost-effective tool for optimising their mechanical assembly designs.


title: On rigidity properties of time-changes of unipotent flows

abstract: We study time-changes of unipotent flows on finite volume quotients of semisimple linear groups, generalising previous work by Ratner on time-changes of horocycle flows. Any measurable isomorphism between time-changes of unipotent flows gives rise to a non-trivial joining supported on its graph. Under a spectral gap assumption on the groups, we show the following rigidity result: either the only limit point of this graph joining under the action of a one-parameter renormalising subgroup is the trivial joining, or the isomorphism is "affine", namely it is obtained composing an algebraic isomorphism with a (non-constant) translation along the centraliser.


title: A Fredholm Alternative for Elliptic with Interior and Boundary Nonlinear Reactions

abstract: In this paper we study the existence of solutions to the following generalized nonlinear two-parameter problem \begin{equation*} a(u, v) \; =\; \lambda\, b(u, m) + \mu\, m(u, v) + \varepsilon\, F(u, v), \end{equation*} for a triple $(a, b, m)$ of continuous, symmetric bilinear forms on a real separable Hilbert space $V$ and nonlinear form $F$. This problem is a natural abstraction of nonlinear problems that occur for a large class of differential operators, various elliptic pde's with nonlinearities in either the differential equation and/or the boundary conditions being a special subclass. First, a Fredholm alternative for the associated linear two-parameter eigenvalue problem is developed, and then this is used to construct a nonlinear version of the Fredholm alternative. Lastly, the Steklov-Robin Fredholm equation is used to exemplify the abstract results.


title: Spontaneous symmetry breaking in the BFSS model: Analytical results using the Gaussian expansion method

abstract: We apply the Gaussian expansion method to the BFSS matrix model in the high temperature limit. When the (Euclidean) BFSS action is expanded about a Gaussian ansatz, it is shown that the SO(9) symmetry is spontaneously broken, analogous to what happens in the IKKT model. The analysis of the free energy, using the set of gap equations which determines the width of the Gaussian terms, is sufficient to show that this symmetry breaking happens only when the fermionic terms are included and is absent in the bosonic case.


title: A PDE approach for regret bounds under partial monitoring

abstract: In this paper, we study a learning problem in which a forecaster only observes partial information. By properly rescaling the problem, we heuristically derive a limiting PDE on Wasserstein space which characterizes the asymptotic behavior of the regret of the forecaster. Using a verification type argument, we show that the problem of obtaining regret bounds and efficient algorithms can be tackled by finding appropriate smooth sub/supersolutions of this parabolic PDE.


title: Decentralized Eigendecomposition for Online Learning over Graphs with Applications

abstract: In this paper, the problem of decentralized eigenvalue decomposition of a general symmetric matrix that is important, e.g., in Principal Component Analysis, is studied, and a decentralized online learning algorithm is proposed. Instead of collecting all information in a fusion center, the proposed algorithm involves only local interactions among adjacent agents. It benefits from the representation of the matrix as a sum of rank-one components which makes the algorithm attractive for online eigenvalue and eigenvector tracking applications. We examine the performance of the proposed algorithm in two types of important application examples: First, we consider the online eigendecomposition of a sample covariance matrix over the network, with application in decentralized Direction-of-Arrival (DoA) estimation and DoA tracking applications. Then, we investigate the online computation of the spectra of the graph Laplacian that is important in, e.g., Graph Fourier Analysis and graph dependent filter design. We apply our proposed algorithm to track the spectra of the graph Laplacian in static and dynamic networks. Simulation results reveal that the proposed algorithm outperforms existing decentralized algorithms both in terms of estimation accuracy as well as communication cost.


title: Object-based active inference

abstract: The world consists of objects: distinct entities possessing independent properties and dynamics. For agents to interact with the world intelligently, they must translate sensory inputs into the bound-together features that describe each object. These object-based representations form a natural basis for planning behavior. Active inference (AIF) is an influential unifying account of perception and action, but existing AIF models have not leveraged this important inductive bias. To remedy this, we introduce 'object-based active inference' (OBAI), marrying AIF with recent deep object-based neural networks. OBAI represents distinct objects with separate variational beliefs, and uses selective attention to route inputs to their corresponding object slots. Object representations are endowed with independent action-based dynamics. The dynamics and generative model are learned from experience with a simple environment (active multi-dSprites). We show that OBAI learns to correctly segment the action-perturbed objects from video input, and to manipulate these objects towards arbitrary goals.


title: Category Theory for Programming

abstract: In these lecture notes, we give a brief introduction to some elements of category theory. The choice of topics is guided by applications to functional programming. Firstly, we study initial algebras, which provide a mathematical characterization of datatypes and recursive functions on them. Secondly, we study monads, which give a mathematical framework for effects in functional languages. The notes include many problems and solutions.


title: A Failure Identification and Recovery Framework for a Planar Reconfigurable Cable Driven Parallel Robot

abstract: In cable driven parallel robots (CDPRs), a single cable malfunction usually induces complete failure of the entire robot. However, the lost static workspace (due to failure) can often be recovered through reconfiguration of the cable attachment points on the frame. This capability is introduced by adding kinematic redundancies to the robot in the form of moving linear sliders that are manipulated in a real-time redundancy resolution controller. The presented work combines this controller with an online failure detection framework to develop a complete fault tolerant control scheme for automatic task recovery. This solution provides robustness by combining pose estimation of the end-effector with the failure detection through the application of an Interactive Multiple Model (IMM) algorithm relying only on end-effector information. The failure and pose estimation scheme is then tied into the redundancy resolution approach to produce a seamless automatic task (trajectory) recovery approach for cable failures.


title: Splashes in isotropic media

abstract: The response of a medium to a sudden localized perturbation (a "splash") will be explained for isotropic media within the framework of linear response theory. In this theory splashes result from the interference of the collective excitations of the medium, with the outcome determined by the interplay between their phase and group velocities as well as the sign of the latter. The salient features of splashes are controlled by the existence of extremal values of the phase and the group velocities: the group velocity gives the expansion rate of the locus of the points where new wavefronts nucleate or existing ones disappear, while the phase velocity determines the large-time expansion rate of a group of wavefronts. If the group velocity is negative in a spectral range and takes on a minimal value within it, then converging wavefronts will be present in the splash. These results are relevant to the studies of several experimentally viable setups, such as a splash on the surface of deep water due to a small pebble or a raindrop, a splash in the two-dimensional electron gas caused by a short voltage pulse applied with the tip of a scanning tunneling microscope, or a bulk splash in superfluid $^{4}He$ due to formation of an electron bubble. Specifically, the gross features of the latter are determined by five extremal velocities. Additionally, due to the existence of a negative group velocity spectral range, some of the wavefronts in superfluid $^{4}He$ splash are converging.


title: On metric approximate subgroups

abstract: Let $G$ be a group with a metric $\mathrm{d}$ invariant under left and right translations, and let $\bar{\mathbb{D}}_r$ be the ball of radius $r$ around the identity. A $(k,r)$-metric approximate subgroup is a symmetric subset $X$ of $G$ such that the pairwise product set $XX$ is covered by at most $k$ translates of $X\bar{\mathbb{D}}_r$. This notion was introduced in arXiv:math/0601431 along with the version for discrete groups (approximate subgroups). In arXiv:0909.2190, it was shown for the discrete case that, at the asymptotic limit of $X$ finite but large, the "approximateness" (or need for more than one translate) can be attributed to a canonically associated Lie group. Here we prove an analogous result in the metric setting, under a certain finite covering assumption on $X$ replacing finiteness. In particular, if $G$ has bounded exponent, we show that any $(k,r)$-metric approximate subgroup is close to a $(1,r')$-metric approximate subgroup for an appropriate $r'$.


title: On static manifolds satisfying an overdetermined Robin type condition on the boundary

abstract: In this work, we consider static manifolds $M$ with nonempty boundary $\partial M$. In this case, we suppose that the potential $V$ also satisfies an overdetermined Robin type condition on $\partial M$. We prove a rigidity theorem for the Euclidean closed unit ball $B^3$ in $\mathbb{R}^3$. More precisely, we give a sharp upper bound for the area of the zero set $\Sigma=V^{-1}(0)$ of the potential $V$, when $\Sigma$ is connected and intersects $\partial M$. We also consider the case where $\Sigma=V^{-1}(0)$ does not intersect $\partial M$.


title: NGC 6240 Supermassive Black Hole Binary dynamical evolution based on Chandra data

abstract: The main idea of our research is to estimate the physical coalescence time of the double supermassive black hole (SMBH) system in the centre of NGC 6240 based on the X-ray observations from the Chandra space observatory. The spectra of the Northern and Southern nuclei were fitted by spectral models from Sherpa and both presented the narrow component of the Fe K$\alpha$ emission line. It enabled us to apply the spectral model to these lines and to find relative offset $\approx0.02$ keV. The enclosed dynamical mass of the central region of NGC 6240 with radius 1 kpc was estimated $\approx 2.04\times 10^{11} \rm\; M_{\odot}$. These data allowed us to carry on the high resolution direct N-body simulations with Newtonian and post-Newtonian (up to $2.5\mathcal{PN}$ correction) dynamics for this particular double SMBH system. As a result, from our numerical models we approximated the central SMBH binary merging time for the different binary eccentricities. In our numerical parameters range the upper limit for the merging time, even for the very small eccentricities, is still below $\approx70$ Myr. Gravitational waveforms and amplitude-frequency pictures from such events can be detected using Pulsar Timing Array (PTA) projects at the last merging phase.


title: Langevin dynamics of generalized spins as SU($N$) coherent states

abstract: Classical models of spin systems traditionally retain only the dipole moments, but a quantum spin state will frequently have additional structure. Spins of magnitude $S$ have $N=2S+1$ levels. Alternatively, the spin state is fully characterized by a set of $N^{2}-1$ local physical observables, which we interpret as generalized spin components. For example, a spin with $S=1$ has three dipole components and five quadrupole components. These components evolve under a generalization of the classical Landau-Lifshitz dynamics, which can be extended with noise and damping terms. In this paper, we reformulate the dynamical equations of motion as a Langevin dynamics of SU($N$) coherent states in the Schr\"odinger picture. This viewpoint is especially useful as the basis for an efficient numerical method to sample spin configurations in thermal equilibrium and to simulate the relaxation and driven motion of topological solitons. To illustrate the approach, we simulate a non-equilibrium relaxation process that creates CP$^{2}$ Skyrmions, which are topological defects with both dipole and quadrupole character.


title: Neuromorphic implementation of ECG anomaly detection using delay chains

abstract: Real-time analysis and classification of bio-signals measured using wearable devices is computationally costly and requires dedicated low-power hardware. One promising approach is to use spiking neural networks implemented using in-memory computing architectures and neuromorphic electronic circuits. However, as these circuits process data in streaming mode without the possibility of storing it in external buffers, a major challenge lies in the processing of spatio-temporal signals that last longer than the time constants present in the network synapses and neurons. Here we propose to extend the memory capacity of a spiking neural network by using parallel delay chains. We show that it is possible to map temporal signals of multiple seconds into spiking activity distributed across multiple neurons which have time constants of few milliseconds. We validate this approach on an ECG anomaly detection task and present experimental results that demonstrate how temporal information is properly preserved in the network activity.


title: Electron Cloud Simulations In The Fermilab Recycler

abstract: We present a simulation study to characterize the stability region of the Fermilab (FNAL) Recycler in the context of secondary emission yield (SEY). Interactions between electrons and beam pipe material can produce electron clouds that jeopardize beam stability in certain focusing configurations. Such an instability was documented in the Recycler, and the work presented here reflects improvements to better understand that finding. We incorporated the Furman-Pivi Model into a PyECLOUD analysis, and we determined the instability threshold given various bunch lengths, intensities, SEY magnitudes, and model parameters.


title: Deep-PANTHER: Learning-Based Perception-Aware Trajectory Planner in Dynamic Environments

abstract: This paper presents Deep-PANTHER, a learning-based perception-aware trajectory planner for unmanned aerial vehicles (UAVs) in dynamic environments. Given the current state of the UAV, and the predicted trajectory and size of the obstacle, Deep-PANTHER generates multiple trajectories to avoid a dynamic obstacle while simultaneously maximizing its presence in the field of view (FOV) of the onboard camera. To obtain a computationally tractable real-time solution, imitation learning is leveraged to train a Deep-PANTHER policy using demonstrations provided by a multimodal optimization-based expert. Extensive simulations show replanning times that are two orders of magnitude faster than the optimization-based expert, while achieving a similar cost. By ensuring that each expert trajectory is assigned to one distinct student trajectory in the loss function, Deep-PANTHER can also capture the multimodality of the problem and achieve a mean squared error (MSE) loss with respect to the expert that is up to 18 times smaller than state-of-the-art (Relaxed) Winner-Takes-All approaches. Deep-PANTHER is also shown to generalize well to obstacle trajectories that differ from the ones used in training.


title: A Two-step Metropolis Hastings Method for Bayesian Empirical Likelihood Computation with Application to Bayesian Model Selection

abstract: In recent times empirical likelihood has been widely applied under Bayesian framework. Markov chain Monte Carlo (MCMC) methods are frequently employed to sample from the posterior distribution of the parameters of interest. However, complex, especially non-convex nature of the likelihood support erects enormous hindrances in choosing an appropriate MCMC algorithm. Such difficulties have restricted the use of Bayesian empirical likelihood (BayesEL) based methods in many applications. In this article, we propose a two-step Metropolis Hastings algorithm to sample from the BayesEL posteriors. Our proposal is specified hierarchically, where the estimating equations determining the empirical likelihood are used to propose values of a set of parameters depending on the proposed values of the remaining parameters. Furthermore, we discuss Bayesian model selection using empirical likelihood and extend our two-step Metropolis Hastings algorithm to a reversible jump Markov chain Monte Carlo procedure to sample from the resulting posterior. Finally, several applications of our proposed methods are presented.


title: On the arboreal jump number of a poset

abstract: A jump is a pair of consecutive elements in an extension of a poset which are incomparable in the original poset. The arboreal jump number is an NP-hard problem that aims to find an arboreal extension of a given poset with minimum number of jumps. The contribution of this paper is twofold: (i)~a characterization that reveals a relation between the number of jumps of an arboreal order extension and the size of a partition of its elements that satisfy some structural properties of the covering graph; (ii)~a compact integer programming model and a heuristic to solve the arboreal jump number problem along with computational results comparing both strategies. The exact method provides an optimality certificate for 18 out of 41 instances with execution time limited to two hours. Furthermore, our heuristic was able to find good feasible solutions for all instances in less than three minutes.


title: Geometry-controlled phase transition in vibrated granular media

abstract: We report experiments on the dynamics of vibrated particles constrained in a two-dimensional vertical container, motivated by the following question: how to get the most out of a given external vibration to maximize internal disorder (e.g. to blend particles) and agitation (e.g. to absorb vibrations)? Granular media are analogs to classical thermodynamic systems, where the injection of energy can be achieved by shaking them: fluidization arises by tuning either the amplitude or the frequency of the oscillations. Alternatively, we explore what happens when another feature, the container geometry, is modified while keeping constant the energy injection. Our method consists in modifying the container base into a V-shape to break the symmetries of the inner particulate arrangement. The lattice contains a compact hexagonal solid-like crystalline phase coexisting with a loose amorphous fluid-like phase, at any thermal agitation. We show that both the solid-to-fluid volume fraction and the granular temperature depend not only on the external vibration but also on the number of topological defects triggered by the asymmetry of the container. The former relies on the statistics of the energy fluctuations and the latter is consistent with a two-dimensional melting transition described by the KTHNY theory.


title: Convergence rate analysis of the gradient descent-ascent method for convex-concave saddle-point problems

abstract: In this paper, we study the gradient descent-ascent method for convex-concave saddle-point problems. We derive a new non-asymptotic global convergence rate in terms of distance to the solution set by using the semidefinite programming performance estimation method. The given convergence rate incorporates most parameters of the problem and it is exact for a large class of strongly convex-strongly concave saddle-point problems for one iteration. We also investigate the algorithm without strong convexity and we provide some necessary and sufficient conditions under which the gradient descent-ascent enjoys linear convergence.


title: Activists Want Better, Safer Technology

abstract: This work surveys 50 Black Lives Matter activists in the United States about the role of technology in organizing protests and other actions. Broad questions about the overall ease and safety of existing technology allowed them to spontaneously offer features and concerns that were most important to them. While the overwhelming majority (84%) of activists reported using social media to find information, 64% had concerns about data privacy or surveillance on social media, and 40% had concerns about the credibility or reliability of information online. Community played a big role in activists' interactions with technology, with 72% reporting that personal networks helped them find information or feel safe attending protests in general. We hope this study will provide a framework for deeper analysis in areas that are important to activists, so technologists can focus on building solutions that fully serve activists' needs and interests.


title: Person Monitoring by Full Body Tracking in Uniform Crowd Environment

abstract: Full body trackers are utilized for surveillance and security purposes, such as person-tracking robots. In the Middle East, uniform crowd environments are the norm which challenges state-of-the-art trackers. Despite tremendous improvements in tracker technology documented in the past literature, these trackers have not been trained using a dataset that captures these environments. In this work, we develop an annotated dataset with one specific target per video in a uniform crowd environment. The dataset was generated in four different scenarios where mainly the target was moving alongside the crowd, sometimes occluding with them, and other times the camera's view of the target is blocked by the crowd for a short period. After the annotations, it was used in evaluating and fine-tuning a state-of-the-art tracker. Our results have shown that the fine-tuned tracker performed better on the evaluation dataset based on two quantitative evaluation metrics, compared to the initial pre-trained tracker.


title: Feature diversity in self-supervised learning

abstract: Many studies on scaling laws consider basic factors such as model size, model shape, dataset size, and compute power. These factors are easily tunable and represent the fundamental elements of any machine learning setup. But researchers have also employed more complex factors to estimate the test error and generalization performance with high predictability. These factors are generally specific to the domain or application. For example, feature diversity was primarily used for promoting syn-to-real transfer by Chen et al. (2021). With numerous scaling factors defined in previous works, it would be interesting to investigate how these factors may affect overall generalization performance in the context of self-supervised learning with CNN models. How do individual factors promote generalization, which includes varying depth, width, or the number of training epochs with early stopping? For example, does higher feature diversity result in higher accuracy held in complex settings other than a syn-to-real transfer? How do these factors depend on each other? We found that the last layer is the most diversified throughout the training. However, while the model's test error decreases with increasing epochs, its diversity drops. We also discovered that diversity is directly related to model width.


title: Distributed Primal-dual Optimization for Heterogeneous Multi-agent Systems

abstract: Heterogeneous networks comprise agents with varying capabilities in terms of computation, storage, and communication. In such settings, it is crucial to factor in the operating characteristics in allowing agents to choose appropriate updating schemes, so as to better distribute computational tasks and utilize the network more efficiently. We consider the multi-agent optimization problem of cooperatively minimizing the sum of local strongly convex objectives. We propose an asynchronous distributed primal-dual protocol, which allows for the primal update steps to be agent-dependent (an agent can opt between first-order or Newton updates). Our analysis introduces a unifying framework for such hybrid optimization scheme and establishes global linear convergence in expectation, under strongly convex objectives and general agent activation schemes. Numerical experiments on real life datasets attest to the merits of the proposed algorithm.


title: IRS-Enabled Backscattering in a Downlink Non-Orthogonal Multiple Access System

abstract: Intelligent reflecting surface (IRS)-enabled backscatter communications can be enabled by an access point (AP) that splits its transmit signal into modulated and unmodulated parts. This letter integrates non-orthogonal multiple access (NOMA) with this method to create a two-user primary system and a secondary system of IRS data. Considering the decoding order, we maximize the rate of the strongest primary user by jointly optimizing the IRS phase shifts, power splitting (PS) factor at the AP, and NOMA power coefficients while guaranteeing the quality of service (QoS) for both weak user and IRS data in the primary and secondary systems, respectively. The resulting optimization problem is non-convex. Thus, we split it into three parts and develop an alternating optimization (AO) algorithm. The advantage is that we derive closed-form solutions for the PS factor and NOMA power coefficients in the first two parts. In the third part, we optimize the phase shifts by exploiting semi-definite relaxation (SDR) and penalty techniques to handle the unit-modulus constraints. This algorithm achieves substantial gains (e.g., 40--68%) compared to relevant baseline schemes.


title: A 4~Kpc Molecular Gas Lane in Cygnus A

abstract: We present the discovery of a 4 kpc molecular gas lane in the Cygnus A host galaxy, using ALMA CO 2-1 observations. The gas lane is oriented roughly perpendicular to the projected radio jet axis. The CO emission generally follows the clumpy dust lanes seen in HST I-band images. The total molecular gas mass is $30\times 10^8$ M$_\odot$ for Milky Way type clouds, and $3.6 \times 10^8$ M$_\odot$ for starburst conditions. There is a velocity change from the northern to southern CO peaks of about $\pm 175$~km~s$^{-1}$, and an apparently smooth velocity gradient between the peaks, although the emission in the central region is weak. In the inner $\sim 0.5"$ projected distance from the radio core, comparison of the CO velocities to those observed for H$_2$ 2.1218 $\mu$m emission shows higher velocities for the vibrationally excited warm molecular gas than the cooler CO 2-1 line emitting gas at similar projected radii. A possible explanation for these different projected velocities at a given radius is that the cooler CO gas is distributed in a clumpy ring at radius $\sim 1.5"$ to $2"$, while the warm H$_2$ 2.12$\mu$m emitting gas is interior to this ring. Of course, the current data cannot rule-out a clumpy, amorphous molecular gas distribution linearly distributed perpendicular to the radio jet axis. We consider surface brightness properties on scales down to $\sim 265$~pc, and discuss the Cygnus A results in the context of other radio galaxies with CO emission.