title: Slow Accretion of Helium Rich Matter onto C-O White Dwarf: Does Composition of WD Matter?

abstract: In binary systems, the helium accretion onto carbon-oxygen (CO) white dwarfs (WDs) plays a vital role in many astrophysical scenarios, especially in supernovae type Ia. Moreover, ignition density for accretion rate $\dot{M} \lesssim 10^{-9} M_{\odot}$ $yr^{-1}$ in helium accreting CO white dwarfs decides the triggering mechanism of a supernova explosion which could be either off-centre helium flash or central carbon flash. We aim to study the accretion of helium with a slow accretion rate $5 \times 10^{-10} M_{\odot}$ $yr^{-1}$ onto relatively cool and hot white dwarfs of different abundances of carbon and oxygen. The simulation code Modules for Experiments in Stellar Astrophysics (MESA) has been used for our study. We analyze the variation in several properties like surface gravity (g), helium luminosity ($L_{He}$), and effective temperature ($T_{eff}$) during the accretion phase of the white dwarfs. We also calculate the ignition density (${\rho}_{He}$) and ignition temperature ($T_{He}$) of helium burning. As expected, the size of WD decreases and g increases during the accretion. However, a red-giant-like expansion is observed after the rapid ignition towards the end. The dependence of helium accreting WD evolution on its composition has also been explored in this study. We find that white dwarfs of the lower abundance of carbon accrete slightly longer before the onset of helium ignition.


title: Ontology-Based Skill Description Learning for Flexible Production Systems

abstract: The increasing importance of resource-efficient production entails that manufacturing companies have to create a more dynamic production environment, with flexible manufacturing machines and processes. To fully utilize this potential of dynamic manufacturing through automatic production planning, formal skill descriptions of the machines are essential. However, generating those skill descriptions in a manual fashion is labor-intensive and requires extensive domain-knowledge. In this contribution an ontology-based semi-automatic skill description system that utilizes production logs and industrial ontologies through inductive logic programming is introduced and benefits and drawbacks of the proposed solution are evaluated.


title: Casimir preserving stochastic Lie-Poisson integrators

abstract: Casimir preserving integrators for stochastic Lie-Poisson equations with Stratonovich noise are developed extending Runge-Kutta Munthe-Kaas methods. The underlying Lie-Poisson structure is preserved along stochastic trajectories. A related stochastic differential equation on the Lie algebra is derived. The solution of this differential equation updates the evolution of the Lie-Poisson dynamics by means of the exponential map. The constructed numerical method conserves Casimir-invariants exactly, which is important for long time integration. This is illustrated numerically for the case of the stochastic heavy top and the stochastic sine-Euler equations.


title: Learning to Search in Task and Motion Planning with Streams

abstract: Task and motion planning problems in robotics combine symbolic planning over discrete task variables with motion optimization over continuous state and action variables. Recent works such as PDDLStream have focused on optimistic planning with an incrementally growing set of objects until a feasible trajectory is found. However, this set is exhaustively expanded in a breadth-first manner, regardless of the logical and geometric structure of the problem at hand, which makes long-horizon reasoning with large numbers of objects prohibitively time-consuming. To address this issue, we propose a geometrically informed symbolic planner that expands the set of objects and facts in a best-first manner, prioritized by a Graph Neural Network that is learned from prior search computations. We evaluate our approach on a diverse set of problems and demonstrate an improved ability to plan in difficult scenarios. We also apply our algorithm on a 7DOF robotic arm in block-stacking manipulation tasks.


title: Unravelling multi-agent ranked delegations

abstract: We introduce a voting model with multi-agent ranked delegations. This model generalises liquid democracy in two aspects: first, an agent's delegation can use the votes of multiple other agents to determine their own -- for instance, an agent's vote may correspond to the majority outcome of the votes of a trusted group of agents; second, agents can submit a ranking over multiple delegations, so that a backup delegation can be used when their preferred delegations are involved in cycles. The main focus of this paper is the study of unravelling procedures that transform the delegation ballots received from the agents into a profile of direct votes, from which a winning alternative can then be determined by using a standard voting rule. We propose and study six such unravelling procedures, two based on optimisation and four using a greedy approach. We study both algorithmic and axiomatic properties, as well as related computational complexity problems of our unravelling procedures for different restrictions on the types of ballots that the agents can submit.


title: Gradient expansions for the large-coupling strength limit of the M{\o}ller-Plesset adiabatic connection

abstract: The adiabatic connection that has as weak-interaction expansion the M{\o}ller-Plesset perturbation series has been recently shown to have a large coupling-strength expansion in terms of functionals of the Hartree-Fock density with a clear physical meaning. In this work we accurately evaluate these density functionals and we extract second-order gradient coefficients from the data for neutral atoms, following ideas similar to the ones used in the literature for exchange, with some modifications. These new gradient expansions will be the key ingredient for performing interpolations that have already been shown to reduce dramatically MP2 errors for large non-covalent complexes. As a byproduct, our investigation of neutral atoms with large number of electrons $N$ indicates that the second-order gradient expansion for exchange grows as $N\log(N)$ rather than as $N$ as often reported in the literature.


title: Fundamental groups of 3-dimensional small covers

abstract: Small covers arising from 3-dimensional simple polytopes are an interesting class of 3-manifolds. The fundamental group is a rigid invariant for wide classes of 3-manifolds, particularly for orientable Haken manifolds, which include orientable small covers. By using Morse-theoretic approach we give a procedure to get an explicit, balanced presentation of the fundamental group of a closed, orientable 3-dimensional simple cover with minimal number of generators. Beside that the minimal Heegaard splitting is determined by this presentation.


title: Well-posedness of singular-degenerate porous medium type equations and application to biofilm models

abstract: We show the well-posedness for a large class of degenerate parabolic equations with an additional singularity and mixed Dirichlet-Neumann boundary conditions on bounded Lipschitz domains. The proof is based on an $L^1$-contraction result. In addition, we analyze systems where degenerate equations are coupled to semilinear reaction diffusion equations. This setting includes mathematical models for biofilm growth which are the motivation for our analysis.


title: A Comparative Analysis of Machine Learning Techniques for IoT Intrusion Detection

abstract: The digital transformation faces tremendous security challenges. In particular, the growing number of cyber-attacks targeting Internet of Things (IoT) systems restates the need for a reliable detection of malicious network activity. This paper presents a comparative analysis of supervised, unsupervised and reinforcement learning techniques on nine malware captures of the IoT-23 dataset, considering both binary and multi-class classification scenarios. The developed models consisted of Support Vector Machine (SVM), Extreme Gradient Boosting (XGBoost), Light Gradient Boosting Machine (LightGBM), Isolation Forest (iForest), Local Outlier Factor (LOF) and a Deep Reinforcement Learning (DRL) model based on a Double Deep Q-Network (DDQN), adapted to the intrusion detection context. The most reliable performance was achieved by LightGBM. Nonetheless, iForest displayed good anomaly detection results and the DRL model demonstrated the possible benefits of employing this methodology to continuously improve the detection. Overall, the obtained results indicate that the analyzed techniques are well suited for IoT intrusion detection.


title: Measuring the Hannay geomatric phase

abstract: An electric network whose dynamics is entirely equivalent to the Foucault pendulum is presented. This circuit allows one to measure in a simulation of the circuit, and eventually in a built circuit in the laboratory, the Hannay geometric phase. This phase corresponds, via the equivalence discussed in the text, exactly to the latitude effect on the rotation of the oscillation plane of the pendulum. The circuit includes a gyrator, that introduces a non-reciprocal character to the network.


title: Computing weakly singular and near-singular integrals over curved boundary elements

abstract: We present algorithms for computing weakly singular and near-singular integrals arising when solving the 3D Helmholtz equation with curved boundary elements. These are based on the computation of the preimage of the singularity in the reference element's space using Newton's method, singularity subtraction, the continuation approach, and transplanted Gauss quadrature. We demonstrate the accuracy of our method for quadratic basis functions and quadratic triangles with several numerical experiments, including the scattering by two half-spheres.


title: Scene Representation Transformer: Geometry-Free Novel View Synthesis Through Set-Latent Scene Representations

abstract: A classical problem in computer vision is to infer a 3D scene representation from few images that can be used to render novel views at interactive rates. Previous work focuses on reconstructing pre-defined 3D representations, e.g. textured meshes, or implicit representations, e.g. radiance fields, and often requires input images with precise camera poses and long processing times for each novel scene. In this work, we propose the Scene Representation Transformer (SRT), a method which processes posed or unposed RGB images of a new area, infers a "set-latent scene representation", and synthesises novel views, all in a single feed-forward pass. To calculate the scene representation, we propose a generalization of the Vision Transformer to sets of images, enabling global information integration, and hence 3D reasoning. An efficient decoder transformer parameterizes the light field by attending into the scene representation to render novel views. Learning is supervised end-to-end by minimizing a novel-view reconstruction error. We show that this method outperforms recent baselines in terms of PSNR and speed on synthetic datasets, including a new dataset created for the paper. Further, we demonstrate that SRT scales to support interactive visualization and semantic segmentation of real-world outdoor environments using Street View imagery.


title: Non-Markovian systems out of equilibrium: Exact results for two routes of coarse graining

abstract: Generalized Langevin equations (GLEs) can be systematically derived via dimensional reduction from high-dimensional microscopic systems. For linear models the derivation can either be based on projection operator techniques such as the Mori-Zwanzig (MZ) formalism or by "integrating out" the bath degrees of freedom. Based on exact analytical results we show that both routes can lead to fundamentally different GLEs and that the origin of these differences is based inherently on the non-equilibrium nature of the microscopic stochastic model. The most important conceptional difference between the two routes is that the MZ result intrinsically fulfills the generalized second fluctuation-dissipation theorem while the integration result can lead to its violation. We supplement our theoretical findings with numerical and simulation results for two popular non-equilibrium systems: Time-delayed feedback control and the active Ornstein-Uhlenbeck process.


title: Country-wide Retrieval of Forest Structure From Optical and SAR Satellite Imagery With Deep Ensembles

abstract: Monitoring and managing Earth's forests in an informed manner is an important requirement for addressing challenges like biodiversity loss and climate change. While traditional in situ or aerial campaigns for forest assessments provide accurate data for analysis at regional level, scaling them to entire countries and beyond with high temporal resolution is hardly possible. In this work, we propose a method based on deep ensembles that densely estimates forest structure variables at country-scale with 10-meter resolution, using freely available satellite imagery as input. Our method jointly transforms Sentinel-2 optical images and Sentinel-1 synthetic-aperture radar images into maps of five different forest structure variables: 95th height percentile, mean height, density, Gini coefficient, and fractional cover. We train and test our model on reference data from 41 airborne laser scanning missions across Norway and demonstrate that it is able to generalize to unseen test regions, achieving normalized mean absolute errors between 11% and 15%, depending on the variable. Our work is also the first to propose a variant of so-called Bayesian deep learning to densely predict multiple forest structure variables with well-calibrated uncertainty estimates from satellite imagery. The uncertainty information increases the trustworthiness of the model and its suitability for downstream tasks that require reliable confidence estimates as a basis for decision making. We present an extensive set of experiments to validate the accuracy of the predicted maps as well as the quality of the predicted uncertainties. To demonstrate scalability, we provide Norway-wide maps for the five forest structure variables.


title: Spectral functions and localization landscape theory in speckle potentials

abstract: Spectral function is a key tool for understanding the behavior of Bose-Einstein condensates of cold atoms in random potentials generated by a laser speckle. In this paper we introduce a new method for computing the spectral functions in disordered potentials. Using a combination of the Wigner-Weyl approach with the landscape theory, we build an approximation for the Wigner distributions of the eigenstates in the phase space and show its accuracy in all regimes, from the deep quantum regime to the intermediate and semiclassical. Based on this approximation, we devise a method to compute the spectral functions using only the landscape-based effective potential. The paper demonstrates the efficiency of the proposed approach for disordered potentials with various statistical properties without requiring any adjustable parameters.


title: Global Interaction Modelling in Vision Transformer via Super Tokens

abstract: With the popularity of Transformer architectures in computer vision, the research focus has shifted towards developing computationally efficient designs. Window-based local attention is one of the major techniques being adopted in recent works. These methods begin with very small patch size and small embedding dimensions and then perform strided convolution (patch merging) in order to reduce the feature map size and increase embedding dimensions, hence, forming a pyramidal Convolutional Neural Network (CNN) like design. In this work, we investigate local and global information modelling in transformers by presenting a novel isotropic architecture that adopts local windows and special tokens, called Super tokens, for self-attention. Specifically, a single Super token is assigned to each image window which captures the rich local details for that window. These tokens are then employed for cross-window communication and global representation learning. Hence, most of the learning is independent of the image patches $(N)$ in the higher layers, and the class embedding is learned solely based on the Super tokens $(N/M^2)$ where $M^2$ is the window size. In standard image classification on Imagenet-1K, the proposed Super tokens based transformer (STT-S25) achieves 83.5\% accuracy which is equivalent to Swin transformer (Swin-B) with circa half the number of parameters (49M) and double the inference time throughput. The proposed Super token transformer offers a lightweight and promising backbone for visual recognition tasks.


title: DA$^{\textbf{2}}$-Net : Diverse & Adaptive Attention Convolutional Neural Network

abstract: Standard Convolutional Neural Network (CNN) designs rarely focus on the importance of explicitly capturing diverse features to enhance the network's performance. Instead, most existing methods follow an indirect approach of increasing or tuning the networks' depth and width, which in many cases significantly increases the computational cost. Inspired by a biological visual system, we propose a Diverse and Adaptive Attention Convolutional Network (DA$^{2}$-Net), which enables any feed-forward CNNs to explicitly capture diverse features and adaptively select and emphasize the most informative features to efficiently boost the network's performance. DA$^{2}$-Net incurs negligible computational overhead and it is designed to be easily integrated with any CNN architecture. We extensively evaluated DA$^{2}$-Net on benchmark datasets, including CIFAR100, SVHN, and ImageNet, with various CNN architectures. The experimental results show DA$^{2}$-Net provides a significant performance improvement with very minimal computational overhead.


title: The Southern Wide-field Gamma-ray Observatory: Status and Prospects

abstract: The Southern Wide-field Gamma-ray Observatory (SWGO) Collaboration is currently engaged in design and prototyping work towards the realisation of this future gamma-ray facility. SWGO will complement CTA and the existing ground-particle based-detectors of the Northern Hemisphere (HAWC and LHAASO) with a very wide field and high duty cycle view of the southern sky. Here I summarise the status of the project and plans for the future, including expectations for sensitivity and science targets as well as the status of the site search and technological developments.


title: Formation of a supergiant quantum vortex in a relativistic Bose-Einstein condensate driven by rotation and a parallel magnetic field

abstract: Analysis based on the energy spectrum of noninteracting bosons shows that, under the circumstance of parallel rotation and magnetic field, charged bosons form a Bose-Einstein condensate because of the lift of the Landau level degeneracy by rotation [\textcolor{blue}{Y. Liu and I. Zahed, Phys. Rev. Lett. {\bf120}, 032001 (2018)}]. In this work, we study the interaction effect on the ground state of this Bose-Einstein condensate of charged bosons from the viewpoint of spontaneous symmetry breaking. We employ a minimal model for charged bosons with repulsive self-interaction. We find that the ground state of such a Bose-Einstein condensate is a supergiant quantum vortex, i.e., a quantized vortex with a large circulation. The size of the vortex is as large as the system size. The low-energy dispersion of the excitation spectra exhibits quadratic behavior, which is an anisotropic realization of the type-II Goldstone boson. Our study may give some implications to off-central relativistic heavy ion collisions, where large vorticity and magnetic fields can be generated.

categories: nucl-th cond-mat.quant-gas cond-mat.supr-con hep-ph


title: Introduction to SPDEs from Probability and PDE

abstract: Lecture notes accompanying a 6 hour mini-course on SPDE given at the University of Oxford in April/May 2021. The main focus of these notes is on an exposition of the variational method for monotone and coercive SPDE. A recap of the necessary functional analysis, operator theory and stochastic analysis on Hilbert spaces is included, with additional references. The final chapter contains a discussion of the pathwise approach.


title: Confined self-propulsion of an isotropic active colloid

abstract: To spontaneously break their intrinsic symmetry and self-propel at the micron scale, isotropic active colloidal particles and droplets exploit the non-linear convective transport of chemical solutes emitted/consumed at their surface by the surface-driven fluid flows generated by these solutes. Significant progress was recently made to understand the onset of self-propulsion and non-linear dynamics. Yet, most models ignore a fundamental experimental feature, namely the spatial confinement of the colloid, and its effect on propulsion. In this work, the self-propulsion of an isotropic colloid inside a capillary tube is investigated numerically. A flexible computational framework is proposed based on a finite-volume approach on adaptative octree-grids and embedded boundary methods. This method is able to account for complex geometric confinement, the nonlinear coupling of chemical transport and flow fields, and the precise resolution of the surface boundary conditions, that drive the system's dynamics. Somewhat counter-intuitively, spatial confinement promotes the colloid's spontaneous motion by reducing the minimum advection-to-diffusion ratio or P\'clet number, Pe, required to self-propel; furthermore, self-propulsion velocities are significantly modified as the colloid-to-capillary size ratio $\kappa$ is increased, reaching a maximum at fixed Pe for an optimal confinement $0<\kappa<1$. These properties stem from a fundamental change in the dominant chemical transport mechanism with respect to the unbounded problem : with diffusion now restricted in most directions by the confining walls, the excess solute is predominantly convected away downstream from the colloid, enhancing front-back concentration contrasts. These results are confirmed quantitatively using conservation arguments and lubrication analysis of the tightly-confined limit, $\kappa\rightarrow 1$.


title: Randomized Stochastic Gradient Descent Ascent

abstract: An increasing number of machine learning problems, such as robust or adversarial variants of existing algorithms, require minimizing a loss function that is itself defined as a maximum. Carrying a loop of stochastic gradient ascent (SGA) steps on the (inner) maximization problem, followed by an SGD step on the (outer) minimization, is known as Epoch Stochastic Gradient \textit{Descent Ascent} (ESGDA). While successful in practice, the theoretical analysis of ESGDA remains challenging, with no clear guidance on choices for the inner loop size nor on the interplay between inner/outer step sizes. We propose RSGDA (Randomized SGDA), a variant of ESGDA with stochastic loop size with a simpler theoretical analysis. RSGDA comes with the first (among SGDA algorithms) almost sure convergence rates when used on nonconvex min/strongly-concave max settings. RSGDA can be parameterized using optimal loop sizes that guarantee the best convergence rates known to hold for SGDA. We test RSGDA on toy and larger scale problems, using distributionally robust optimization and single-cell data matching using optimal transport as a testbed.


title: Semantic-Aware Generation for Self-Supervised Visual Representation Learning

abstract: In this paper, we propose a self-supervised visual representation learning approach which involves both generative and discriminative proxies, where we focus on the former part by requiring the target network to recover the original image based on the mid-level features. Different from prior work that mostly focuses on pixel-level similarity between the original and generated images, we advocate for Semantic-aware Generation (SaGe) to facilitate richer semantics rather than details to be preserved in the generated image. The core idea of implementing SaGe is to use an evaluator, a deep network that is pre-trained without labels, for extracting semantic-aware features. SaGe complements the target network with view-specific features and thus alleviates the semantic degradation brought by intensive data augmentations. We execute SaGe on ImageNet-1K and evaluate the pre-trained models on five downstream tasks including nearest neighbor test, linear classification, and fine-scaled image recognition, demonstrating its ability to learn stronger visual representations.


title: Neural network stochastic differential equation models with applications to financial data forecasting

abstract: In this article, we employ a collection of stochastic differential equations with drift and diffusion coefficients approximated by neural networks to predict the trend of chaotic time series which has big jump properties. Our contributions are, first, we propose a model called L\'evy induced stochastic differential equation network, which explores compounded stochastic differential equations with $\alpha$-stable L\'evy motion to model complex time series data and solve the problem through neural network approximation. Second, we theoretically prove that the numerical solution through our algorithm converges in probability to the solution of corresponding stochastic differential equation, without curse of dimensionality. Finally, we illustrate our method by applying it to real financial time series data and find the accuracy increases through the use of non-Gaussian L\'evy processes. We also present detailed comparisons in terms of data patterns, various models, different shapes of L\'evy motion and the prediction lengths.


title: Bianchi type I cosmology with a Kalb-Ramond background field

abstract: An effect of the Lorentz symmetry breaking is pointed out in the cosmological context. Using a Bianchi I geometry coupled to the Kalb-Ramond field, a consequence of the Lorentz symmetry violation is indicated by a different rate of expansion in a given spatial direction. This article focuses on the coupling constant $\xi_1$, which generates, from the Kalb-Ramond field, all three coefficients that give rise to the Lorentz violation in the gravity sector of the minimal Standard Model Extension. The coupling constant $\xi_1$ increases the rate of expansion of the universe in a given direction during a dark energy era. As a consequence, a range of validity of that coupling constant is also obtained.


title: Towards a Precision Measurement of Binary Black Holes Formation Channels Using Gravitational Waves and Emission Lines

abstract: The formation of compact objects-neutron stars, black holes, and supermassive black holes-and its connection to the chemical composition of the galaxies is one of the central questions in astrophysics. We propose a novel data-driven, multi-messenger technique to address this question by exploiting the inevitable correlation between gravitational waves and atomic/molecular emission line signals. For a fiducial probability distribution function $p(t_d)\propto t_d^{-\kappa}$ of time delays, this method can provide a measurement of the minimum delay time of $0.5$ Gyr and the power-law index $\kappa=1$ with a standard deviation $0.12$ (and $0.45$) and $0.06$ (and $0.34$), respectively from five years of LIGO-Virgo-KAGRA observation in synergy with SPHEREx line intensity mapping (and DESI emission-line galaxies). Such measurements will provide data-driven multi-messenger constraints on the delay time distribution which is currently not well known.

categories: astro-ph.GA astro-ph.CO astro-ph.HE gr-qc


title: An efficient IMEX-DG solver for the compressible Navier-Stokes equations for non-ideal gases

abstract: We propose an efficient, accurate and robust IMEX solver for the compressible Navier-Stokes equation describing non-ideal gases with general equations of state. The method, which is based on an $h-$adaptive Discontinuos Galerkin spatial discretization and on an Additive Runge Kutta IMEX method for time discretization, is tailored for low Mach number applications and allows to simulate low Mach regimes at a significantly reduced computational cost, while maintaining full second order accuracy also for higher Mach number regimes. The method has been implemented in the framework of the $deal.II$ numerical library, whose adaptive mesh refinement capabilities are employed to enhance efficiency. Refinement indicators appropriate for real gas phenomena have been introduced. A number of numerical experiments on classical benchmarks for compressible flows and their extension to real gases demonstrate the properties of the proposed method.


title: A Markov chain for metapopulations of small sizes with attraction landscape

abstract: Mathematical models represent one of the fundamental ways of studying nature. In special, epidemic models have shown to be particularly useful in the understanding of the course of diseases and in the planning effective control policies. A particular type of epidemic model considers the individuals divided into populations. When studied in graphs, it is already known that the graph topology can play an important role in the evolution of the disease. At the same time, one may want to study the effect of the presence of an underlying \emph{attraction landscape} of the vertices, apart from the respectively underlying topology. In this work, we study metapopulations with small number of individuals in the presence of an attraction landscape. Individuals move across populations and get infected according to the SIS compartmental model. By using a Markov chain approach, we provide a numerical approximation to the prediction of the long-term prevalence of the disease. More specifically, an approach that combines two binomial distributions for mobility, with appropriate assumptions, is proposed to approximate the model. The problem setting is simulated through Monte-Carlo experiments and the obtained results are compared to the mathematic-analytical approach. Substantial agreement is observed between both approaches, which corroborates the effectiveness of the reported numerical approach. In addition, we also study the impact of different levels of attraction landscapes, as well as propagation on the local scale of the entire population. All in all, this study proposes a potentially effective approach to a mostly unexplored setting of disease transmission.


title: Bosonic $\eta$-deformations of Non-integrable Backgrounds

abstract: We consider the non-integrable bosonic backgrounds $W_{2,4}\times T^{1,1}$ and $AdS_5\times T^{1,1}$ and derive their bosonic $\eta$-deformed versions using an $r$-matrix that solves the modified Yang-Baxter equation obtaining new integrable deformed backgrounds.


title: Universal equations for maximal isotropic Grassmannians

abstract: The isotropic Grassmannian parametrizes isotropic subspaces of a vector space equipped with a quadratic form. In this paper, we show that any maximal isotropic Grassmannian in its Pl\"ucker embedding can be defined by pulling back the equations of $Gr_{\operatorname{iso}}(3,7)$ or $Gr_{\operatorname{iso}}(4,8)$.


title: Intrinsic Dimension, Persistent Homology and Generalization in Neural Networks

abstract: Disobeying the classical wisdom of statistical learning theory, modern deep neural networks generalize well even though they typically contain millions of parameters. Recently, it has been shown that the trajectories of iterative optimization algorithms can possess fractal structures, and their generalization error can be formally linked to the complexity of such fractals. This complexity is measured by the fractal's intrinsic dimension, a quantity usually much smaller than the number of parameters in the network. Even though this perspective provides an explanation for why overparametrized networks would not overfit, computing the intrinsic dimension (e.g., for monitoring generalization during training) is a notoriously difficult task, where existing methods typically fail even in moderate ambient dimensions. In this study, we consider this problem from the lens of topological data analysis (TDA) and develop a generic computational tool that is built on rigorous mathematical foundations. By making a novel connection between learning theory and TDA, we first illustrate that the generalization error can be equivalently bounded in terms of a notion called the 'persistent homology dimension' (PHD), where, compared with prior work, our approach does not require any additional geometrical or statistical assumptions on the training dynamics. Then, by utilizing recently established theoretical results and TDA tools, we develop an efficient algorithm to estimate PHD in the scale of modern deep neural networks and further provide visualization tools to help understand generalization in deep learning. Our experiments show that the proposed approach can efficiently compute a network's intrinsic dimension in a variety of settings, which is predictive of the generalization error.


title: Security Threats and Cellular Network Procedures for Unmanned Aircraft Systems

abstract: This paper discusses cellular network security for unmanned aircraft systems (UASs) and provides insights into the ongoing Third Generation Partnership Project (3GPP) standardization efforts with respect to authentication and authorization, location information privacy, and command and control signaling. We introduce the 3GPP reference architecture for network connected UAS and the new network functions as part of the 5G core network, discuss introduce the three security contexts, potential threats, and the 3GPP procedures. The paper identifies research opportunities for UAS communications security and recommends critical security features and processes to be considered for standardization.


title: A new lower bound on the pebbling number of the grid

abstract: A pebbling move on a graph consists of removing $2$ pebbles from a vertex and adding $1$ pebble to one of the neighbouring vertices. A vertex is called reachable if we can put $1$ pebble on it after a sequence of moves. The optimal pebbling number of a graph is the minimum number $m$ such that there exists a distribution of $m$ pebbles so that each vertex is reachable. For the case of a square grid $n \times m$, Gy\H{o}ri, Katona and Papp recently showed that its optimal pebbling number is at least $\frac{2}{13}nm \approx 0.1538nm$ and at most $\frac{2}{7}nm +O(n+m) \approx 0.2857nm$. We improve the lower bound to $\frac{5092}{28593}nm +O(m+n) \approx 0.1781nm$.


title: Partial Eruption, Confinement, and Twist Buildup and Release of a Double-decker Filament

abstract: We investigate the failed partial eruption of a filament system in NOAA AR 12104 on 2014 July 5, using multiwavelength EUV, magnetogram, and H$\alpha$ observations, as well as magnetic field modeling. The filament system consists of two almost co-spatial segments with different end points, both resembling a C shape. Following an ejection and a precursor flare related to flux cancellation, only the upper segment rises and then displays a prominent twisted structure, while rolling over toward its footpoints. The lower segment remains undisturbed, indicating that the system possesses a double-decker structure. The erupted segment ends up with a reverse-C shape, with material draining toward its footpoints, while losing its twist. Using the flux rope insertion method, we construct a model of the source region that qualitatively reproduces key elements of the observed evolution. At the eruption onset, the model consists of a flux rope atop a flux bundle with negligible twist, which is consistent with the observational interpretation that the filament possesses a double-decker structure. The flux rope reaches the critical height of the torus instability during its initial relaxation, while the lower flux bundle remains in stable equilibrium. The eruption terminates when the flux rope reaches a dome-shaped quasi-separatrix layer that is reminiscent of a magnetic fan surface, although no magnetic null is found. The flux rope is destroyed by reconnection with the confining overlying flux above the dome, transferring its twist in the process.


title: Homogeneous Low-Resolution Face Recognition Method based Correlation Features

abstract: Face recognition technology has been widely adopted in many mission-critical scenarios like means of human identification, controlled admission, and mobile device access, etc. Security surveillance is a typical scenario of face recognition technology. Because the low-resolution feature of surveillance video and images makes it difficult for high-resolution face recognition algorithms to extract effective feature information, Algorithms applied to high-resolution face recognition are difficult to migrate directly to low-resolution situations. As face recognition in security surveillance becomes more important in the era of dense urbanization, it is essential to develop algorithms that are able to provide satisfactory performance in processing the video frames generated by low-resolution surveillance cameras. This paper study on the Correlation Features-based Face Recognition (CoFFaR) method which using for homogeneous low-resolution surveillance videos, the theory, experimental details, and experimental results are elaborated in detail. The experimental results validate the effectiveness of the correlation features method that improves the accuracy of homogeneous face recognition in surveillance security scenarios.


title: Using Color To Identify Insider Threats

abstract: Insider threats are costly, hard to detect, and unfortunately rising in occurrence. Seeking to improve detection of such threats, we develop novel techniques to enable us to extract powerful features and augment attack vectors for greater classification power. Most importantly, we generate high quality color image encodings of user behavior that do not have the downsides of traditional greyscale image encodings. Combined, they form Computer Vision User and Entity Behavior Analytics, a detection system designed from the ground up to improve upon advancements in academia and mitigate the issues that prevent the usage of advanced models in industry. The proposed system beats state-of-art methods used in academia and as well as in industry on a gold standard benchmarking dataset.


title: Finite-density-induced motility and turbulence of chimera solitons

abstract: We consider a one-dimensional oscillatory medium with a coupling through a diffusive linear field. In the limit of fast diffusion this setup reduces to the classical Kuramoto-Battogtokh model. We demonstrate that for a finite diffusion stable chimera solitons, namely localized synchronous domain in an infinite asynchronous environment, are possible. The solitons are stable also for finite density of oscillators, but in this case they sway with a nearly constant speed. This finite-density-induced motility disappears in the continuum limit, as the velocity of the solitons is inverse proportional to the density. A long-wave instability of the homogeneous asynchronous state causes soliton turbulence, which appears as a sequence of soliton mergings and creations. As the instability of the asynchronous state becomes stronger, this turbulence develops into a spatio-temporal intermittency.


title: A Versatile Optimization Framework For Sustainable Post-Disaster Building Reconstruction

abstract: This paper proposes an optimization framework for sustainable post-disaster building reconstruction. Based on mathematical optimization, it is intended to provide decision-makers with a versatile tool to optimize building designs and to explore the trade-off between costs and environmental impact (represented by embodied energy) of alternative building materials. The mixed-integer nonlinear optimization model includes an analytical building model that considers structural and safety constraints and incorporates regional building codes. Using multi-objective optimization concepts, Pareto-optimal designs are computed that represent the best trade-off designs from which a decision-maker can choose when they take additional criteria into consideration. As a case study, we consider the design of a multi-room one-story masonry building in Nepal. We demonstrate how the framework can be employed to answer a variety of questions, such as the optimal building design and material selection, the sensitivity of the decision to the prices, and the impact of regional safety regulation thresholds.


title: A Sphere Packing Bound for Vector Gaussian Fading Channels under Peak Amplitude Constraints

abstract: An upper bound on the capacity of multiple-input multiple-output (MIMO) Gaussian fading channels is derived under peak amplitude constraints. The upper bound is obtained borrowing concepts from convex geometry and it extends to MIMO channels notable results from the geometric analysis on the capacity of scalar Gaussian channels. Relying on a sphere packing argument and on the renowned Steiner's formula, the proposed upper bound depends on the intrinsic volumes of the constraint region, i.e., functionals defining a measure of the geometric features of a convex body. The tightness of the bound is investigated at high signal-to-noise ratio (SNR) for any arbitrary convex amplitude constraint region, for any channel matrix realization, and any dimension of the MIMO system. In addition, two variants of the upper bound are proposed: one is useful to ensure the feasibility in the evaluation of the bound and the other to improve the bound's performance in the low SNR regime. Finally, the upper bound is specialized for two practical transmitter configurations, either employing a single power amplifier for all transmitting antennas or a power amplifier for each antenna.


title: Variational Gibbs Inference for Statistical Model Estimation from Incomplete Data

abstract: Statistical models are central to machine learning with broad applicability across a range of downstream tasks. The models are controlled by free parameters that are typically estimated from data by maximum-likelihood estimation or approximations thereof. However, when faced with real-world data sets many of the models run into a critical issue: they are formulated in terms of fully-observed data, whereas in practice the data sets are plagued with missing data. The theory of statistical model estimation from incomplete data is conceptually similar to the estimation of latent-variable models, where powerful tools such as variational inference (VI) exist. However, in contrast to standard latent-variable models, parameter estimation with incomplete data often requires estimating exponentially-many conditional distributions of the missing variables, hence making standard VI methods intractable. We address this gap by introducing variational Gibbs inference (VGI), a new general-purpose method to estimate the parameters of statistical models from incomplete data. We validate VGI on a set of synthetic and real-world estimation tasks, estimating important machine learning models such as variational autoencoders and normalising flows from incomplete data. The proposed method, whilst general-purpose, achieves competitive or better performance than existing model-specific estimation methods.


title: Embeddings and Integrable Charges for Extended Corner Symmetry

abstract: We revisit the problem of extending the phase space of diffeomorphism-invariant theories to account for embeddings associated with the boundary of sub-regions. We do so by emphasizing the importance of a careful treatment of embeddings in all aspects of the covariant phase space formalism. In so doing we introduce a new notion of the extension of field space associated with the embeddings which has the important feature that the Noether charges associated with all extended corner symmetries are in fact integrable, but not necessarily conserved. We give an intuitive understanding of this description. We then show that the charges give a representation of the extended corner symmetry via the Poisson bracket, without central extension.


title: Time-resolved reversible optical switching of the ultralow-loss phase change material Sb2Se3

abstract: The antimony-based chalcogenide Sb2Se3 is a rapidly emerging material for photonic phase change applications owing to its ultra-low optical losses at telecommunication wavelengths in both crystalline and amorphous phases. Here, we investigate the dynamical response of these materials from nanoseconds to milliseconds under optical pumping conditions. We apply bichromatic pump-probe transient reflectance spectroscopy which is a widely used method to study the optical performance of optical phase change materials. Amorphous regions of several hundreds of nanometers in diameter are induced by pulsed excitation of the material using a wavelength of 488 nm above the absorption edge, while the transient reflectance is probed using a continuous wave 980 nm laser, well below the absorption edge of the material. We find vitrification dynamics in the nanosecond range and observe crystallization on millisecond time scales. These results show a large five-orders of magnitude difference in time scales between crystallization and vitrification dynamics in this material. The insights provided in this work are fundamental for the optimisation of the material family and its employment in photonic applications.


title: Effect of laser surface hardening on a low carbon steel substrate wear and hardness

abstract: Laser surface hardening (LSH) is an effective process to produce hard and wear resistant steel surfaces. This work produced hardened surfaces using a fiber laser in a bare (B. uncoated) and carbon-coated (C) conditions. The obtained hardness and microstructure were compatible to martensite (B) and cementite (C) major constituents. A low coefficient of friction was obtained using C-condition. The treated surface also presented little damage during a reciprocating wear against a hard metal sphere.


title: Multiple target tracking with interaction using an MCMC MRF Particle Filter

abstract: This paper presents and discusses an implementation of a multiple target tracking method, which is able to deal with target interactions and prevent tracker failures due to hijacking. The referenced approach uses a Markov Chain Monte Carlo (MCMC) sampling step to evaluate the filter and constructs an efficient proposal density to generate new samples. This density integrates target interaction terms based on Markov Random Fields (MRFs) generated per time step. The MRFs model the interactions between targets in an attempt to reduce tracking ambiguity that typical particle filters suffer from when tracking multiple targets. A test sequence of 662 grayscale frames containing 20 interacting ants in a confined space was used to test both the proposed approach and a set of importance sampling based independent particle filters, to establish a performance comparison. It is shown that the implemented approach of modeling target interactions using MRF successfully corrects many of the tracking errors made by the independent, interaction unaware, particle filters.


title: Learning Conditional Invariance through Cycle Consistency

abstract: Identifying meaningful and independent factors of variation in a dataset is a challenging learning task frequently addressed by means of deep latent variable models. This task can be viewed as learning symmetry transformations preserving the value of a chosen property along latent dimensions. However, existing approaches exhibit severe drawbacks in enforcing the invariance property in the latent space. We address these shortcomings with a novel approach to cycle consistency. Our method involves two separate latent subspaces for the target property and the remaining input information, respectively. In order to enforce invariance as well as sparsity in the latent space, we incorporate semantic knowledge by using cycle consistency constraints relying on property side information. The proposed method is based on the deep information bottleneck and, in contrast to other approaches, allows using continuous target properties and provides inherent model selection capabilities. We demonstrate on synthetic and molecular data that our approach identifies more meaningful factors which lead to sparser and more interpretable models with improved invariance properties.


title: Federated Data Science to Break Down Silos [Vision]

abstract: Similar to Open Data initiatives, data science as a community has launched initiatives for sharing not only data but entire pipelines, derivatives, artifacts, etc. (Open Data Science). However, the few efforts that exist focus on the technical part on how to facilitate sharing, conversion, etc. This vision paper goes a step further and proposes KEK, an open federated data science platform that does not only allow for sharing data science pipelines and their (meta)data but also provides methods for efficient search and, in the ideal case, even allows for combining and defining pipelines across platforms in a federated manner. In doing so, KEK addresses the so far neglected challenge of actually finding artifacts that are semantically related and that can be combined to achieve a certain goal.


title: Information Bottleneck-Based Hebbian Learning Rule Naturally Ties Working Memory and Synaptic Updates

abstract: Artificial neural networks have successfully tackled a large variety of problems by training extremely deep networks via back-propagation. A direct application of back-propagation to spiking neural networks contains biologically implausible components, like the weight transport problem or separate inference and learning phases. Various methods address different components individually, but a complete solution remains intangible. Here, we take an alternate approach that avoids back-propagation and its associated issues entirely. Recent work in deep learning proposed independently training each layer of a network via the information bottleneck (IB). Subsequent studies noted that this layer-wise approach circumvents error propagation across layers, leading to a biologically plausible paradigm. Unfortunately, the IB is computed using a batch of samples. The prior work addresses this with a weight update that only uses two samples (the current and previous sample). Our work takes a different approach by decomposing the weight update into a local and global component. The local component is Hebbian and only depends on the current sample. The global component computes a layer-wise modulatory signal that depends on a batch of samples. We show that this modulatory signal can be learned by an auxiliary circuit with working memory (WM) like a reservoir. Thus, we can use batch sizes greater than two, and the batch size determines the required capacity of the WM. To the best of our knowledge, our rule is the first biologically plausible mechanism to directly couple synaptic updates with a WM of the task. We evaluate our rule on synthetic datasets and image classification datasets like MNIST, and we explore the effect of the WM capacity on learning performance. We hope our work is a first-step towards understanding the mechanistic role of memory in learning.


title: BioLeaF: A Bio-plausible Learning Framework for Training of Spiking Neural Networks

abstract: Our brain consists of biological neurons encoding information through accurate spike timing, yet both the architecture and learning rules of our brain remain largely unknown. Comparing to the recent development of backpropagation-based (BP-based) methods that are able to train spiking neural networks (SNNs) with high accuracy, biologically plausible methods are still in their infancy. In this work, we wish to answer the question of whether it is possible to attain comparable accuracy of SNNs trained by BP-based rules with bio-plausible mechanisms. We propose a new bio-plausible learning framework, consisting of two components: a new architecture, and its supporting learning rules. With two types of cells and four types of synaptic connections, the proposed local microcircuit architecture can compute and propagate error signals through local feedback connections and support training of multi-layers SNNs with a globally defined spiking error function. Under our microcircuit architecture, we employ the Spike-Timing-Dependent-Plasticity (STDP) rule operating in local compartments to update synaptic weights and achieve supervised learning in a biologically plausible manner. Finally, We interpret the proposed framework from an optimization point of view and show the equivalence between it and the BP-based rules under a special circumstance. Our experiments show that the proposed framework demonstrates learning accuracy comparable to BP-based rules and may provide new insights on how learning is orchestrated in biological systems.


title: Humanode Whitepaper: You are [not] a bot

abstract: The advent of blockchain technology has led to a massive wave of different decentralized ledger technology (DLT) solutions. Such projects as Bitcoin and Ethereum have shifted the paradigm of how to transact value in a decentralized manner, but their various core technologies have their own advantages and disadvantages. This paper aims to describe an alternative to modern decentralized financial networks by introducing the Humanode network. Humanode is a network safeguarded by cryptographically secure bio-authorized nodes. Users will be able to deploy nodes by staking their encrypted biometric data. This approach can potentially lead to the creation of a public, permissionless financial network based on consensus between equal human nodes with algorithm-based emission mechanisms targeting real value growth and proportional emission. Humanode combines different technological stacks to achieve a decentralized, secure, scalable, efficient, consistent, immutable, and sustainable financial network: 1) a bio-authorization module based on cryptographically secure neural networks for the private classification of 3D templates of users' faces 2) a private Liveness detection mechanism for identification of real human beings 3) a Substrate module as a blockchain layer 4) a cost-based fee system 5) a Vortex decentralized autonomous organization (DAO) governing system 6) a monetary policy and algorithm, Fath, where monetary supply reacts to real value growth and emission is proportional. All of these implemented technologies have nuances that are crucial for the integrity of the network. In this paper we address these details, describing problems that might occur and their possible solutions. The main goal of Humanode is to create a stable and just financial network that relies on the existence of human life.


title: Observations on holographic aspects of four-dimensional asymptotically flat ${\cal N}=2$ black holes

abstract: In this note, we explore holographic attributes of four-dimensional near-extremal Reissner-Nordstrom black hole solutions in ungauged ${\cal N}=2$ supergravity theories at the two-derivative level by recasting them as a specific first-order deformation in solution space, associated with an infinitesimal Harrison transformation, of black holes in an $AdS_2$ space-time. Specifically, we use this link to exhibit how bulk properties, such as mass and entropy, of four-dimensional near-extremal black holes are holographically encoded in the one-dimensional boundary theory dual to gravity in an infinitesimally deformed $AdS_2$ space-time. We do so for the case of four-dimensional near-extremal black holes that arise as deformations in solution space of BPS black holes by changing the non-extremality parameter. For these near-extremal black holes, we further show that the nAdS$_2$ attractor mechanism can be recast as a specific deformation of the BPS flow equations in four dimensions. Additionally, we also discuss time-dependent perturbations of the four-dimensional near-extremal Reissner-Nordstrom solutions from a two-dimensional point of view.


title: Analysis of emergent patterns in crossing flows of pedestrians reveals an invariant of `stripe' formation in human data

abstract: When two streams of pedestrians cross at an angle, striped patterns spontaneously emerge as a result of local pedestrian interactions. This clear case of self-organized pattern formation remains to be elucidated. In counterflows, with a crossing angle of 180{\deg}, alternating lanes of traffic are commonly observed moving in opposite directions, whereas in crossing flows at an angle of 90{\deg} diagonal stripes have been reported. Naka (1977) hypothesized that stripe orientation is perpendicular to the bisector of the crossing angle. However, studies of crossing flows at acute and obtuse angles remain underdeveloped. We tested the bisector hypothesis in experiments on small groups (18-19 participants each) crossing at seven angles (30{\deg} intervals), and analyzed the geometric properties of stripes. We present two novel computational methods for analyzing striped patterns in pedestrian data: (i) an edge-cutting algorithm, which detects the dynamic formation of stripes and allows us to measure local properties of individual stripes; and (ii) a pattern-matching technique, based on the Gabor function, which allows us to estimate global properties (orientation and wavelength) of the striped pattern at a time T. We find an invariant property: stripes in the two groups are parallel and perpendicular to the bisector at all crossing angles. In contrast, other properties depend on the crossing angle: stripe spacing (wavelength), stripe size (number of pedestrians per stripe), and crossing time all decrease as the crossing angle increases from 30{\deg} to 180{\deg}, whereas the number of stripes increases with crossing angle. We also observe that the width of individual stripes is dynamically squeezed as the two groups cross each other. The findings thus support the bisector hypothesis at a wide range of crossing angles, although the theoretical reasons for this invariant remain unclear.


title: On a class of Cheeger inequalities

abstract: We study a general version of the Cheeger inequality by considering the shape functional $\mathcal{F}_{p,q}(\Omega)=\lambda_p^{1/p}(\Omega)/\lambda_q(\Omega)^{1/q}$. The infimum and the supremum of $\mathcal{F}_{p,q}$ are studied in the class of all domains $\Omega$ of $\mathbb{R}^d$ and in the subclass of convex domains. In the latter case the issue concerning the existence of an optimal domain for $\mathcal{F}_{p,q}$ is discussed.


title: Attitude determination for nano-satellites -- II. Dead reckoning with a multiplicative extended Kalman filter

abstract: This paper is the second part of a series of studies discussing a novel attitude determination method for nano-satellites. Our approach is based on the utilization of thermal imaging sensors to determine the direction of the Sun and the nadir with respect to the satellite with sub-degree accuracy. The proposed method is planned to be applied during the Cubesats Applied for MEasuring and LOcalising Transients (CAMELOT) mission aimed at detecting and localizing gamma-ray bursts with an efficiency and accuracy comparable to large gamma-ray space observatories. In this paper we introduce a simulation model aimed at testing the applicability of our attitude determination approach. Its first part simulates the orbit and rotation of a satellite with arbitrary initial conditions while its second part applies our attitude determination algorithm which is based on a multiplicative extended Kalman filter. The simulated satellite is assumed to be equipped with a GPS system, MEMS gyroscopes and the infrasensors. These instruments provide the required data input for the Kalman filter. We demonstrate the applicability of our attitude determination algorithm by simulating the motion of a nano-satellite on Low Earth Orbit. Our results show that the attitude determination may have a 1$\sigma$ error of $\sim30'$ even with a large gyroscope drift during the orbital periods when the infrasensors provide both the direction of the Sun and the Earth (the nadir). This accuracy is an improvement on the point source detection accuracy of the infrasensors. However, the attitude determination error can get as high as 25$^{\circ}$ during periods when the Sun is occulted by the Earth. We show that following an occultation period the attitude information is immediately recovered by the Kalman filter once the Sun is observed again.

categories: astro-ph.IM astro-ph.EP astro-ph.SR physics.space-ph


title: Laser frequency offset locking at 10-Hz-level instability using hybrid electronic filters

abstract: Lasers with well controlled relative frequencies are indispensable for many applications in science and technology. We present a frequency offset locking method for lasers based on beat frequency discrimination utilizing hybrid electronic LC filters. The method is specifically designed for decoupling the tightness of the lock from the broadness of its capture range. The presented demonstration locks two free running diode lasers at 780 nm with a 5.5 GHz offset. It displays an offset frequency instability below 55 Hz for timescales in excess of 1000 s and a minimum of 12 Hz at 10 s averaging, outperforming the best reported instabilities of methods based on beat frequency discrimination. The performance is complemented with a 190 MHz lock capture range, a tuning range of up to 1 GHz, and a frequency ramp agility of 200 kHz/$\mu$s.


title: A categorification of the colored Jones polynomial at a root of unity

abstract: There is a $p$-differential on the triply-graded Khovanov--Rozansky homology of knots and links over a field of positive characteristic $p$ that gives rise to an invariant in the homotopy category finite-dimensional $p$-complexes. A differential on triply-graded homology discovered by Cautis is compatible with the $p$-differential structure. As a consequence we get a categorification of the colored Jones polynomial evaluated at a $2p$th root of unity.


title: SwinBERT: End-to-End Transformers with Sparse Attention for Video Captioning

abstract: The canonical approach to video captioning dictates a caption generation model to learn from offline-extracted dense video features. These feature extractors usually operate on video frames sampled at a fixed frame rate and are often trained on image/video understanding tasks, without adaption to video captioning data. In this work, we present SwinBERT, an end-to-end transformer-based model for video captioning, which takes video frame patches directly as inputs, and outputs a natural language description. Instead of leveraging multiple 2D/3D feature extractors, our method adopts a video transformer to encode spatial-temporal representations that can adapt to variable lengths of video input without dedicated design for different frame rates. Based on this model architecture, we show that video captioning can benefit significantly from more densely sampled video frames as opposed to previous successes with sparsely sampled video frames for video-and-language understanding tasks (e.g., video question answering). Moreover, to avoid the inherent redundancy in consecutive video frames, we propose adaptively learning a sparse attention mask and optimizing it for task-specific performance improvement through better long-range video sequence modeling. Through extensive experiments on 5 video captioning datasets, we show that SwinBERT achieves across-the-board performance improvements over previous methods, often by a large margin. The learned sparse attention masks in addition push the limit to new state of the arts, and can be transferred between different video lengths and between different datasets. Code is available at https://github.com/microsoft/SwinBERT


title: Quantum state transfer on the complete bipartite graph

abstract: Previously it was shown that (almost) perfect state transfer can be achieved on the complete bipartite graph by a discrete-time coined quantum walk based algorithm when both the sender and receiver vertices are in the same partition of the graph and when the sender and receiver are in opposite partitions of the same size. By changing the coin operator, we analyze the state transfer problem and we show that it is still possible to achieve state transfer with high fidelity even when the sender and receiver are in different partitions with different sizes. Moreover, it is also possible to use an active switch approach using lackadaisical quantum walks where the marked vertex is switched between the sender and receiver during the algorithm.


title: The Implicit Graph Conjecture is False

abstract: An efficient implicit representation of an $n$-vertex graph $G$ in a family $\mathcal{F}$ of graphs assigns to each vertex of $G$ a binary code of length $O(\log n)$ so that the adjacency between every pair of vertices can be determined only as a function of their codes. This function can depend on the family but not on the individual graph. Every family of graphs admitting such a representation contains at most $2^{O(n\log(n))}$ graphs on $n$ vertices, and thus has at most factorial speed of growth. The Implicit Graph Conjecture states that, conversely, every hereditary graph family with at most factorial speed of growth admits an efficient implicit representation. We refute this conjecture by establishing the existence of hereditary graph families with factorial speed of growth that require codes of length $n^{\Omega(1)}$.


title: The concentration-compactness principle for Orlicz spaces and applications

abstract: In this paper we extend the well-known concentration -- compactness principle of P.L. Lions to Orlicz spaces. As an application we show an existence result to some critical elliptic problem with nonstandard growth.


title: A large-deviations principle for all the components in a sparse inhomogeneous random graph

abstract: We study an inhomogeneous sparse random graph on [N] = {1, . . . , N } as introduced in a seminal paper by Bollobas, Janson and Riordan (2007): vertices have a type (here in a compact metric space S), and edges between different vertices occur randomly and independently over all vertex pairs, with a probability depending on the two vertex types. In the limit N to infinity, we consider the sparse regime, where the average degree is O(1). We prove a large-deviations principle with explicit rate function for the statistics of the collection of all the connected components, registered according to their vertex type sets, and distinguished according to being microscopic (of finite size) or macroscopic (of size proportional to N). In doing so, we derive explicit logarithmic asymptotics for the probability that GN is connected. We present a full analysis of the rate function including its minimizers. From this analysis we deduce a number of limit laws, conditional and unconditional, which provide comprehensive information about all the microscopic and macroscopic components of the graph. In particular, we recover the criterion for the existence of the phase transition given in [5].


title: Fracture mechanical behavior of polymers: 1. Amorphous glassy state

abstract: Theoretical analyses and experiments have been carried out to investigate fracture behavior of glassy polymers. Our birefringence measurements quantify the local stress buildup at cut tip during different stages of drawing. Based on polymethyl methacrylate (PMMA), bisphenol A polycarbonate (PC) and polyethylene terephthalate (PET), we find several key results beyond the existing knowledge base. (1) The inherent fracture and yield strengths sigma_F(inh) and sigma_Y(inh) differ little in magnitude from the breaking and yield stress (sigma_b and sigma_y). (2) Stress intensification (SI) near a pre-through-cut builds up deviates from the theoretical description of linear elastic fracture mechanics (LEFM) upon approaching the cut tip. (3) SI meets a natural cutoff below which stress ceases to increase. (4) The stress stip at cut tip shows a trend of approximate linear increase with the far-field load s0 for all three polymers and different cut size a. (5) A characteristic length scale P emerges from the linear relation between stip and KI. For these glassy polymers, P is on the order of 0.1 mm. (6) Fracture toughness of brittle polymers is characterized by critical stress intensity factor K_Ic = sigma_F(inh)(2*pi*P)1/2, revealing relevance of the two crucial quantities. (7) The critical energy release rate GIc for brittle glass polymers such as PMMA is determined by the product of its work of fracture wF (of uncut specimen) and P. (8) The elusive fractocohesive length Lfc defined in the literature as G_Ic/w_F naturally arises from the new expression for G_Ic as stated in (7), i.e., it is essentially P. These results suggest that a great deal of future work is required to acquire additional understanding with regards to fracture and failure behaviors of plastics.


title: Killing spinors and hypersurfaces

abstract: We consider spin manifolds with an Einstein metric, either Riemannian or indefinite, for which there exists a Killing spinor. We describe the intrinsic geometry of nondegenerate hypersurfaces in terms of a PDE satisfied by a pair of induced spinors, akin to the generalized Killing spinor equation. Conversely, we prove an embedding result for real analytic pseudo-Riemannian manifolds carrying a pair of spinors satisfying this condition.


title: Online Sampling and Decision Making with Low Entropy

abstract: Consider the problem: we are given $n$ boxes, labeled $\{1,2,\ldots, n\}$ by an adversary, each containing a single number chosen from an unknown distribution; these $n$ distributions are not necessarily identical. We are also given an integer $k \leq n$. We have to choose an order in which we will sequentially open these boxes, and each time we open the next box in this order, we learn the number in the box. Once we reject a number in a box, the box cannot be recalled. Our goal is to accept $k$ of these numbers, without necessarily opening all boxes, such that the accepted numbers are the $k$ largest numbers in the boxes (thus their sum is maximized). A natural approach to solve such problems is to use randomness to sample randomly ordered elements, however, as indicated in several sources, e.g., Turan et al. NIST'15, Bierhorst et al. Nature'18, pure randomness is hard to get in reality. We present an algorithm for this problem, which is provably and simultaneously near-optimal with respect to the achieved competitive ratio and the used amount of randomness. In particular, we construct a distribution on the orders with entropy $\Theta(\log\log n)$ such that a deterministic multiple-threshold algorithm gives a competitive ratio $1-O(\sqrt{\log k/k})$, for $k < \log n/\log \log n$. Our competitive ratio is simultaneously optimal and uses optimal entropy $\Theta(\log\log n)$, improving in three ways the previous best known algorithm, whose competitive ratio is $1 - O(1/k^{1/3}) - o(1)$.


title: BaLeNAS: Differentiable Architecture Search via the Bayesian Learning Rule

abstract: Differentiable Architecture Search (DARTS) has received massive attention in recent years, mainly because it significantly reduces the computational cost through weight sharing and continuous relaxation. However, more recent works find that existing differentiable NAS techniques struggle to outperform naive baselines, yielding deteriorative architectures as the search proceeds. Rather than directly optimizing the architecture parameters, this paper formulates the neural architecture search as a distribution learning problem through relaxing the architecture weights into Gaussian distributions. By leveraging the natural-gradient variational inference (NGVI), the architecture distribution can be easily optimized based on existing codebases without incurring more memory and computational consumption. We demonstrate how the differentiable NAS benefits from Bayesian principles, enhancing exploration and improving stability. The experimental results on NAS-Bench-201 and NAS-Bench-1shot1 benchmark datasets confirm the significant improvements the proposed framework can make. In addition, instead of simply applying the argmax on the learned parameters, we further leverage the recently-proposed training-free proxies in NAS to select the optimal architecture from a group architectures drawn from the optimized distribution, where we achieve state-of-the-art results on the NAS-Bench-201 and NAS-Bench-1shot1 benchmarks. Our best architecture in the DARTS search space also obtains competitive test errors with 2.37\%, 15.72\%, and 24.2\% on CIFAR-10, CIFAR-100, and ImageNet datasets, respectively.


title: Fractionalized holes in one-dimensional $\mathbb{Z}_2$ gauge theory coupled to fermion matter -- deconfined dynamics and emergent integrability

abstract: We investigate the interplay of quantum one-dimensional discrete $\mathbb{Z}_2$ gauge fields and fermion matter near full filling in terms of deconfined fractionalized hole excitations that constitute mobile domain walls between vacua that break spontaneously translation symmetry. In the limit of strong string tension, we uncover emergent integrable correlated hopping dynamics of holes which is complementary to the constrained XXZ description in terms of bosonic dimers. We analyze numerically quantum dynamics of spreading of an isolated hole together with the associated time evolution of entanglement and provide analytical understanding of its salient features. We also study the model enriched with a short-range interaction and clarify the nature of the resulting ground state at low filling of holes and identify deconfined hole excitations near the hole filling $\nu^h=1/3$.

categories: cond-mat.str-el cond-mat.quant-gas hep-lat hep-th


title: Length of stationary Gaussian excursions

abstract: Given that a stationary Gaussian process is above a high threshold, the length of time it spends before going below that threshold is studied. The asymptotic order is determined by the smoothness of the sample paths, which in turn is a function of the tails of the spectral measure. Two disjoint regimes are studied - one in which the second spectral moment is finite and the other in which the tails of the spectral measure are regularly varying and the second moment is infinite.


title: Characteristic Neural Ordinary Differential Equations

abstract: We propose Characteristic-Neural Ordinary Differential Equations (C-NODEs), a framework for extending Neural Ordinary Differential Equations (NODEs) beyond ODEs. While NODEs model the evolution of a latent variables as the solution to an ODE, C-NODE models the evolution of the latent variables as the solution of a family of first-order quasi-linear partial differential equations (PDEs) along curves on which the PDEs reduce to ODEs, referred to as characteristic curves. This in turn allows the application of the standard frameworks for solving ODEs, namely the adjoint method. Learning optimal characteristic curves for given tasks improves the performance and computational efficiency, compared to state of the art NODE models. We prove that the C-NODE framework extends the classical NODE on classification tasks by demonstrating explicit C-NODE representable functions not expressible by NODEs. Additionally, we present C-NODE-based continuous normalizing flows, which describe the density evolution of latent variables along multiple dimensions. Empirical results demonstrate the improvements provided by the proposed method for classification and density estimation on CIFAR-10, SVHN, and MNIST datasets under a similar computational budget as the existing NODE methods. The results also provide empirical evidence that the learned curves improve the efficiency of the system through a lower number of parameters and function evaluations compared with baselines.


title: Evaluation of Interpretability for Deep Learning algorithms in EEG Emotion Recognition: A case study in Autism

abstract: Current models on Explainable Artificial Intelligence (XAI) have shown an evident and quantified lack of reliability for measuring feature-relevance when statistically entangled features are proposed for training deep classifiers. There has been an increase in the application of Deep Learning in clinical trials to predict early diagnosis of neuro-developmental disorders, such as Autism Spectrum Disorder (ASD). However, the inclusion of more reliable saliency-maps to obtain more trustworthy and interpretable metrics using neural activity features is still insufficiently mature for practical applications in diagnostics or clinical trials. Moreover, in ASD research the inclusion of deep classifiers that use neural measures to predict viewed facial emotions is relatively unexplored. Therefore, in this study we propose the evaluation of a Convolutional Neural Network (CNN) for electroencephalography (EEG)-based facial emotion recognition decoding complemented with a novel RemOve-And-Retrain (ROAR) methodology to recover highly relevant features used in the classifier. Specifically, we compare well-known relevance maps such as Layer-Wise Relevance Propagation (LRP), PatternNet, Pattern-Attribution, and Smooth-Grad Squared. This study is the first to consolidate a more transparent feature-relevance calculation for a successful EEG-based facial emotion recognition using a within-subject-trained CNN in typically-developed and ASD individuals.


title: Mitigating Noise-Induced Gradient Vanishing in Variational Quantum Algorithm Training

abstract: Variational quantum algorithms are expected to demonstrate the advantage of quantum computing on near-term noisy quantum computers. However, training such variational quantum algorithms suffers from gradient vanishing as the size of the algorithm increases. Previous work cannot handle the gradient vanishing induced by the inevitable noise effects on realistic quantum hardware. In this paper, we propose a novel training scheme to mitigate such noise-induced gradient vanishing. We first introduce a new cost function of which the gradients are significantly augmented by employing traceless observables in truncated subspace. We then prove that the same minimum can be reached by optimizing the original cost function with the gradients from the new cost function. Experiments show that our new training scheme is highly effective for major variational quantum algorithms of various tasks.


title: Testing the quasicentroid molecular dynamics method on gas-phase ammonia

abstract: Quasicentroid molecular dynamics (QCMD) is a path-integral method for approximating nuclear quantum effects in dynamics simulations, which has given promising results for gas- and condensed-phase water. Here, by simulating the infrared spectrum of gas-phase ammonia, we test the feasibility of extending QCMD beyond water. Overall, QCMD works as well for ammonia as for water, reducing or eliminating blue shifts from the classical spectrum without introducing the artificial red-shifts or broadening associated with other imaginary-time path-integral methods. However, QCMD gives only a modest improvement over the classical spectrum for the position of the symmetric bend mode, which is highly anharmonic (since it correlates with the inversion pathway). We expect QCMD to have similar problems with large-amplitude degrees of freedom in other molecules, but otherwise to work as well as for water.


title: A Family of Complex Kleinian Split Solvable Groups

abstract: It is shown that lattices of a family of split solvable subgroups of PSL(N + 1, C) are complex Kleinian using techniques of Lie groups and dynamical systems, also that there exists a minimal limit set for the action of these lattices on the N dimensional complex projective space and that there are exactly two maximal discontinuity regions.


title: Searching for stop LSP at the LHC

abstract: We analyze relevant signals expected at the LHC for a stop as the lightest supersymmetric particle (LSP). The discussion is carried out in the framework of the $\mu\nu$SSM, where the presence of $R$-parity violating couplings involving right-handed neutrinos solves the $\mu$-problem and reproduces neutrino data. The stops are pair produced at the LHC, decaying with displaced vertices to a lepton and a bottom quark or a neutrino and a top quark. We compare the predictions of this scenario with ATLAS and CMS searches for long-lived particles. To analyze the parameter space we sample the $\mu\nu$SSM for a stop LSP using a likelihood data-driven method, and paying special attention to reproduce the current experimental data on neutrino and Higgs physics, as well as flavor observables. Our results translate into a lower limit on the mass of the left (right) stop LSP of 1068 GeV (1341 GeV), corresponding to a decay length of 1.86 mm (6.61 mm).


title: OTB-morph: One-Time Biometrics via Morphing applied to Face Templates

abstract: Cancelable biometrics refers to a group of techniques in which the biometric inputs are transformed intentionally using a key before processing or storage. This transformation is repeatable enabling subsequent biometric comparisons. This paper introduces a new scheme for cancelable biometrics aimed at protecting the templates against potential attacks, applicable to any biometric-based recognition system. Our proposed scheme is based on time-varying keys obtained from morphing random biometric information. An experimental implementation of the proposed scheme is given for face biometrics. The results confirm that the proposed approach is able to withstand against leakage attacks while improving the recognition performance.


title: Toric and tropical Bertini theorems in positive characteristic

abstract: We generalize the toric Bertini theorem of Fuchs, Mantova, and Zannier to positive characteristic. A key part of the proof is a new algebraically closed field containing the field \kk(t_1,\dots,t_d) of rational functions over an algebraically closed field \kk of prime characteristic. As a corollary, we extend the tropical Bertini theorem of Maclagan and Yu to arbitrary characteristic, which removes the characteristic dependence from the d-connectivity result for tropical varieties from that paper.


title: Bootstrapping Closed Hyperbolic Surfaces

abstract: The eigenvalues of the Laplace-Beltrami operator and the integrals of products of eigenfunctions and holomorphic $s$-differentials satisfy certain consistency conditions on closed hyperbolic surfaces. These consistency conditions can be derived by using spectral decompositions to write quadruple overlap integrals in terms of triple overlap integrals in different ways. We show how to efficiently construct these consistency conditions and use them to derive upper bounds on eigenvalues, following the approach of the conformal bootstrap. As an example of such a bootstrap bound, we find a numerical upper bound on the spectral gap of closed orientable hyperbolic surfaces that is nearly saturated by the Bolza surface.


title: Cross-Domain Adaptive Teacher for Object Detection

abstract: We address the task of domain adaptation in object detection, where there is a domain gap between a domain with annotations (source) and a domain of interest without annotations (target). As an effective semi-supervised learning method, the teacher-student framework (a student model is supervised by the pseudo labels from a teacher model) has also yielded a large accuracy gain in cross-domain object detection. However, it suffers from the domain shift and generates many low-quality pseudo labels (\textit{e.g.,} false positives), which leads to sub-optimal performance. To mitigate this problem, we propose a teacher-student framework named Adaptive Teacher (AT) which leverages domain adversarial learning and weak-strong data augmentation to address the domain gap. Specifically, we employ feature-level adversarial training in the student model, allowing features derived from the source and target domains to share similar distributions. This process ensures the student model produces domain-invariant features. Furthermore, we apply weak-strong augmentation and mutual learning between the teacher model (taking data from the target domain) and the student model (taking data from both domains). This enables the teacher model to learn the knowledge from the student model without being biased to the source domain. We show that AT demonstrates superiority over existing approaches and even Oracle (fully-supervised) models by a large margin. For example, we achieve 50.9% (49.3%) mAP on Foggy Cityscape (Clipart1K), which is 9.2% (5.2%) and 8.2% (11.0%) higher than previous state-of-the-art and Oracle, respectively.


title: A Uniformly Selected, Southern-Sky 6dF, Optical AGN Catalog

abstract: We have constructed a catalog of active galactic nuclei (AGNs) with $z < 0.13$, based on optical spectroscopy, from the parent sample of galaxies in the 6dF galaxy survey (Final Release of 6dFGS), a census of the Southern hemisphere. This work is an extension of our all sky AGN catalog \citet [ZCF, here after]{ZCF19}. The ZCF is based on 43,533 galaxies with $\rm K_s \leq$ 11.75 ($z \leq 0.09$) in the 2MASS Redshift Survey (2MRS). The parent catalog of this work, the 6dF catalog, consists of 136,304 publicly available digital spectra for 125,071 galaxies with $\rm Dec \leq 0^\circ$ and $\rm K_s \leq 12.65$ (median $z = 0.053$). Our AGN catalog consists of 3109 broad line AGNs and 12,156 narrow line AGNs which satisfy the \citet{Kauffmann03} criteria, of which 3865 also satisfy the \citet{Kewley01} criteria. We also provide emission line widths, fluxes, flux errors, and signal-to-noise ratios of all the galaxies in our spectroscopic sample, allowing users to customize the selection criteria. In addition, we provide AGN likelihood for the rest of galaxies based on the availability and quality of their spectra. These likelihood values can be used for rigorous statistical analyses.


title: Contextuality and Wigner negativity are equivalent for continuous-variable quantum measurements

abstract: Quantum computers will provide considerable speedups with respect to their classical counterparts. However, the identification of the innately quantum features that enable these speedups is challenging. In the continuous-variable setting - a promising paradigm for the realisation of universal, scalable, and fault-tolerant quantum computing - contextuality and Wigner negativity have been perceived as two such distinct resources. Here we show that they are in fact equivalent for the standard models of continuous-variable quantum computing. While our results provide a unifying picture of continuous-variable resources for quantum speedup, they also pave the way towards practical demonstrations of continuous-variable contextuality, and shed light on the significance of negative probabilities in phase-space descriptions of quantum mechanics.


title: Differentially private stochastic expectation propagation (DP-SEP)

abstract: We are interested in privatizing an approximate posterior inference algorithm called Expectation Propagation (EP). EP approximates the posterior by iteratively refining approximations to the local likelihoods, and is known to provide better posterior uncertainties than those by variational inference (VI). However, EP needs a large memory to maintain all local approximates associated with each datapoint in the training data. To overcome this challenge, stochastic expectation propagation (SEP) considers a single unique local factor that captures the average effect of each likelihood term to the posterior and refines it in a way analogous to EP. In terms of privacy, SEP is more tractable than EP because at each refining step of a factor, the remaining factors are fixed and do not depend on other datapoints as in EP, which makes the sensitivity analysis straightforward. We provide a theoretical analysis of the privacy-accuracy trade-off in the posterior estimates under our method, called differentially private stochastic expectation propagation (DP-SEP). Furthermore, we demonstrate the performance of our DP-SEP algorithm evaluated on both synthetic and real-world datasets in terms of the quality of posterior estimates at different levels of guaranteed privacy.


title: Mono-chromatic single photon events at the muon collider

abstract: The cross section for lepton pair annihilation into a photon and a dark photon or an axion-like particle is constant for large center-of-mass energies because some of the portal operators coupling Standard Model and dark sector are proportional to the energy. Feebly coupled though they are, these portal operators will be enhanced by the large center-of-mass energy made available by a muon collider and thus provide the ideal example of possible physics beyond the Standard Model to be studied with such a machine. We discuss the characteristic signature of the presence of these operators: mono-chromatic single photon events for the two benchmarks of having center-of-mass energies of 3 and 10 TeV and integrated luminosity of, respectively, 1 and 10 ab$^{-1}$. We find that an effective scale of the portal operator as large as $\Lambda=112$ TeV for an axion-like particle and $\Lambda=141$ TeV for a dark photon can be separated from the background with a confidence level of 95% in the first benchmark; these interaction scales can be raised to $\Lambda=375$ TeV and $\Lambda=459$ TeV in the case of the second benchmark. The signal for the pseudo scalar particle can be distinguished from that of the spin-1 with about 500 events. The response of the detector to high-energy photons is examined.


title: Families of solutions of the heterotic G$_2$ system

abstract: We construct new families of solutions of the heterotic G$_2$ system on squashed homogeneous 3-Sasakian manifolds, that is, using squashed metrics on either the 7-sphere or the Aloff-Wallach space $N_{1,1}$. We obtain AdS$_3$ solutions for all values of the squashing parameter $s$ except for the nearly-parallel G$_2$ value $s=1/\sqrt{5}$, for which we don't find any solutions. Along the process, we construct different G$_2$-instanton connections on bundles over these squashed manifolds.


title: Quantum Motif Clustering

abstract: We present three quantum algorithms for clustering graphs based on higher-order patterns, known as motif clustering. One uses a straightforward application of Grover search, the other two make use of quantum approximate counting, and all of them obtain square-root like speedups over the fastest classical algorithms in various settings. In order to use approximate counting in the context of clustering, we show that for general weighted graphs the performance of spectral clustering is mostly left unchanged by the presence of constant (relative) errors on the edge weights. Finally, we extend the original analysis of motif clustering in order to better understand the role of multiple `anchor nodes' in motifs and the types of relationships that this method of clustering can and cannot capture.


title: The coupling flow of ${\cal N}=4$ super Yang-Mills theory

abstract: We offer a novel perspective on ${\cal N}=4$ supersymmetric Yang-Mills (SYM) theory through the framework of the Nicolai map, a transformation of the bosonic fields that allows one to compute quantum correlators in terms of a free, purely bosonic functional measure. Generally, any Nicolai map is obtained through a path-ordered exponential of the so-called coupling flow operator. The latter can be canonically constructed in any gauge using an ${\cal N}=1$ off-shell superfield formulation of ${\cal N}=4$ SYM, or alternatively through dimensional reduction of the result from ${\cal N}=1$ $D=10$ SYM, in which case we need to restrict to the Landau gauge. We propose a general theory of the ${\cal N}=4$ coupling flow operator, arguing that it exhibits an ambiguity in form of an R-symmetry freedom given by the Lie algebra $\mathfrak{su}(4)$. This theory incorporates our two construction approaches as special points in $\mathfrak{su}(4)$ and defines a broad class of Nicolai maps for ${\cal N}=4$ SYM.


title: Solving systems of Boolean multivariate equations with quantum annealing

abstract: Polynomial systems over the binary field have important applications, especially in symmetric and asymmetric cryptanalysis, multivariate-based post-quantum cryptography, coding theory, and computer algebra. In this work, we study the quantum annealing model for solving Boolean systems of multivariate equations of degree 2, usually referred to as the Multivariate Quadratic problem. We present different methodologies to embed the problem into a Hamiltonian that can be solved by available quantum annealing platforms. In particular, we provide three embedding options, and we highlight their differences in terms of quantum resources. Moreover, we design a machine-agnostic algorithm that adopts an iterative approach to better solve the problem Hamiltonian by repeatedly reducing the search space. Finally, we use D-Wave devices to successfully implement our methodologies on several instances of the Multivariate Quadratic problem.


title: Boundary Heisenberg Algebras and Their Deformations

abstract: We investigate the deformations and rigidity of boundary Heisenberg-like algebras. In particular, we focus on the Heisenberg and $\text{Heisenberg}\oplus\mathfrak{witt}$ algebras which arise as symmetry algebras in three-dimensional gravity theories. As a result of the deformation procedure we find a large class of algebras. While some of these algebras are new, some of them have already been obtained as asymptotic and boundary symmetry algebras, supporting the idea that symmetry algebras associated to diverse boundary conditions and spacetime loci are algebraically interconnected through deformation of algebras. The deformation/contraction relationships between the new algebras are investigated. In addition, it is also shown that the deformation procedure reaches new algebras inaccessible to the Sugawara construction. As a byproduct of our analysis, we obtain that $\text{Heisenberg}\oplus\mathfrak{witt}$ and the asymptotic symmetry algebra Weyl-$\mathfrak{bms}_3$ are not connected via single deformation but in a more subtle way.


title: A Kernel Test for Causal Association via Noise Contrastive Backdoor Adjustment

abstract: Causal inference grows increasingly complex as the number of confounders increases. Given treatments $X$, confounders $Z$ and outcomes $Y$, we develop a non-parametric method to test the \textit{do-null} hypothesis $H_0:\; p(y|\text{\it do}(X=x))=p(y)$ against the general alternative. Building on the Hilbert Schmidt Independence Criterion (HSIC) for marginal independence testing, we propose backdoor-HSIC (bd-HSIC) and demonstrate that it is calibrated and has power for both binary and continuous treatments under a large number of confounders. Additionally, we establish convergence properties of the estimators of covariance operators used in bd-HSIC. We investigate the advantages and disadvantages of bd-HSIC against parametric tests as well as the importance of using the do-null testing in contrast to marginal independence testing or conditional independence testing. A complete implementation can be found at \hyperlink{https://github.com/MrHuff/kgformula}{\texttt{https://github.com/MrHuff/kgformula}}.


title: Spectral analysis and stabilization of the dissipative Schr\"{o}dinger operator on the tadpole graph

abstract: We consider the damped Schr\"odinger semigroup $e^{-it \frac{d^2}{dx^2}}$ on the tadpole graph ${\mathcal R}$. We first give a careful spectral analysis and an appropriate decomposition of the kernel of the resolvent. As a consequence and by showing that the generalized eigenfunctions form a Riesz basis of some subspace of $L^2({\mathcal R})$, we prove that the corresponding energy decay exponentially.


title: Securities Lending Haircuts and Indemnification Pricing

abstract: Securities borrowing and lending are critical to proper functioning of securities markets. To alleviate securities owners' exposure to borrower default risk, overcollateralization and indemnification are provided by the borrower and the lending agent respectively. Haircuts as the level of overcollateralization and the cost of indemnification are naturally interrelated: the higher haircut is, the lower cost shall become. This article presents a method of quantifying their relationship. Borrower dependent haircuts satisfying the lender's credit risk appetite are computed for US Treasuries and main equities by applying a repo haircut model to bilateral securities lending transactions. Indemnification is designed to fulfill a triple-A risk appetite when the transaction haircut fails to deliver. The cost of indemnification consists of a risk charge, a capital charge, and a funding charge, each corresponding to the expected loss, the economic capital, and the redundant fund needed to arrive at the triple-A haircut.


title: Generalizing Clinical Trials with Convex Hulls

abstract: Randomized clinical trials eliminate confounding but impose strict exclusion criteria that limit recruitment to a subset of the population. Observational datasets are more inclusive but suffer from confounding -- often providing overly optimistic estimates of treatment response over time due to partially optimized physician prescribing patterns. We therefore assume that the unconfounded treatment response lies somewhere in-between the observational estimate before and the observational estimate after treatment assignment. This assumption allows us to extrapolate results from exclusive trials to the broader population by analyzing observational and trial data simultaneously using an algorithm called Optimum in Convex Hulls (OCH). OCH represents the treatment effect either in terms of convex hulls of conditional expectations or convex hulls (also known as mixtures) of conditional densities. The algorithm first learns the component expectations or densities using the observational data and then learns the linear mixing coefficients using trial data in order to approximate the true treatment effect; theory importantly explains why this linear combination should hold. OCH estimates the treatment effect in terms both expectations and densities with state of the art accuracy.