title: New examples and partial classification of 15-vertex triangulations of the quaternionic projective plane

abstract: Brehm and K\"uhnel (1992) constructed three 15-vertex combinatorial 8-manifolds `like the quaternionic projective plane' with symmetry groups $\mathrm{A}_5$, $\mathrm{A}_4$, and $\mathrm{S}_3$, respectively. Gorodkov (2016) proved that these three manifolds are in fact PL homeomorphic to $\mathbb{HP}^2$. Note that 15 is the minimal number of vertices of a combinatorial 8-manifold that is not PL homeomorphic to $S^8$. In the present paper we construct a lot of new 15-vertex triangulations of $\mathbb{HP}^2$. A surprising fact is that such examples are found for very different symmetry groups, including those not in any way related to the group $\mathrm{A}_5$. Namely, we find 19 triangulations with symmetry group $\mathrm{C}_7$, one triangulation with symmetry group $\mathrm{C}_6\times\mathrm{C}_2$, 14 triangulations with symmetry group $\mathrm{C}_6$, 26 triangulations with symmetry group $\mathrm{C}_5$, one new triangulation with symmetry group $\mathrm{A}_4$, and 11 new triangulations with symmetry group $\mathrm{S}_3$. Further, we obtain the following classification result. We prove that, up to isomorphism, there are exactly 75 triangulations of $\mathbb{HP}^2$ with 15 vertices and symmetry group of order at least 4: the three Brehm-K\"uhnel triangulations and the 72 new triangulations listed above. On the other hand, we show that there are plenty of triangulations with symmetry groups $\mathrm{C}_3$ and $\mathrm{C}_2$, as well as the trivial symmetry group.


title: The Global Existence for the Derivative Nonlinear Schr\"odinger Equation with solitons

abstract: The global existence of the solution for the second-type derivative nonlinear Schr\"odinger (DNLSII) equation with solitons is presented for the first time on the line with weighted Sobolev initial data in $H^2( \mathbb{R}) \cap H^{1,1}(\mathbb{R} )$. The modified Darboux transformation is applied as a main research tool to analyze the global existence for the Cauchy problem by adding or subtracting a finite number of zeros of the scattering data. For this purpose, it is necessary to study the invertibility and uniqueness of the Darboux transformation first acting on the weighted Sobolev space. The analysis of time evolution comprehensively demonstrates the existence of a unique global solution to the DNLSII equation with Cauchy problem in $H^2( \mathbb{R}) \cap H^{1,1}(\mathbb{R} )$.


title: Actively Learning Numerical Relativity

abstract: Data analysis of gravitational waves detected by the Ligo-Virgo-Kagra collaboration and future observatories relies on precise modelling of the sources. In order to build, calibrate and validate current models, we resort to expensive simulations in Numerical Relativity (NR), the fully-fledged simulation of Einstein's Equations. Since simulation costs and the dimensionality of parameter space are prohibitive to perform a dense coverage, approximate models interpolate among the available simulation data. We put forward the technique of Gaussian Process Active Learning (GPAL), an adaptive, data-driven protocol, for parameter space exploration and training of gravitational wave approximants. We evaluate this proposal by studying a computationally inexpensive scenario, in which we calibrate the approximant TEOBResumS using the NR-informed model as a proxy for NR. In this case study, we find that GPAL reduces the computational cost of training by a factor of 4 with respect to uniform or randomly distributed simulations. Moreover, we consider a parallel implementation which reduces computational time, and hybrid strategies which improve pre-calibrated models. The Gaussian Process regression employed in this approach naturally endows the algorithm with notion of model uncertainty. We comment on the implications of this feature for data analysis.


title: Optimizing rgb-d semantic segmentation through multi-modal interaction and pooling attention

abstract: Semantic segmentation of RGB-D images involves understanding the appearance and spatial relationships of objects within a scene, which requires careful consideration of various factors. However, in indoor environments, the simple input of RGB and depth images often results in a relatively limited acquisition of semantic and spatial information, leading to suboptimal segmentation outcomes. To address this, we propose the Multi-modal Interaction and Pooling Attention Network (MIPANet), a novel approach designed to harness the interactive synergy between RGB and depth modalities, optimizing the utilization of complementary information. Specifically, we incorporate a Multi-modal Interaction Fusion Module (MIM) into the deepest layers of the network. This module is engineered to facilitate the fusion of RGB and depth information, allowing for mutual enhancement and correction. Additionally, we introduce a Pooling Attention Module (PAM) at various stages of the encoder. This module serves to amplify the features extracted by the network and integrates the module's output into the decoder in a targeted manner, significantly improving semantic segmentation performance. Our experimental results demonstrate that MIPANet outperforms existing methods on two indoor scene datasets, NYUDv2 and SUN-RGBD, underscoring its effectiveness in enhancing RGB-D semantic segmentation.


title: Symbolic Execution for Quantum Error Correction Programs

abstract: We define QSE, a symbolic execution framework for quantum programs by integrating symbolic variables into quantum states and the outcomes of quantum measurements. The soundness of QSE is established through a theorem that ensures the correctness of symbolic execution within operational semantics. We further introduce symbolic stabilizer states, which symbolize the phases of stabilizer generators, for the efficient analysis of quantum error correction (QEC) programs. Within the QSE framework, we can use symbolic expressions to characterize the possible discrete Pauli errors in QEC, providing a significant improvement over existing methods that rely on sampling with simulators. We implement QSE with the support of symbolic stabilizer states in a prototype tool named QuantumSE.jl. Our experiments on representative QEC codes, including quantum repetition codes, Kitaev's toric codes, and quantum Tanner codes, demonstrate the efficiency of QuantumSE.jl for debugging QEC programs with over 1000 qubits. In addition, by substituting concrete values in symbolic expressions of measurement results, QuantumSE.jl is also equipped with a sampling feature for stabilizer circuits. Despite a longer initialization time than the state-of-the-art stabilizer simulator, Google's Stim, QuantumSE.jl offers a quicker sampling rate in the experiments.


title: Simulation of Matrix Product States to Unveil the Initial State Dependency of non-Gaussian Dynamics of Kerr Nonlinearity

abstract: We simulate a free dissipative and coherent-driven Kerr nonlinear system using a time-evolving block decimation (TEBD) algorithm, to study the impact of the initial state on the exact quantum dynamics of the system. The superposition of two coherent branches results in non-classical time dynamics. The Wigner state representation confirms that the system ends up saturating to two different branches, through evolving different trajectories, resulting in de-Gaussification throughout evolution. Furthermore, we also see that the time evolution suffers a residual effect of the initial state.


title: TPTU-v2: Boosting Task Planning and Tool Usage of Large Language Model-based Agents in Real-world Systems

abstract: Large Language Models (LLMs) have demonstrated proficiency in addressing tasks that necessitate a combination of task planning and the usage of external tools that require a blend of task planning and the utilization of external tools, such as APIs. However, real-world complex systems present three prevalent challenges concerning task planning and tool usage: (1) The real system usually has a vast array of APIs, so it is impossible to feed the descriptions of all APIs to the prompt of LLMs as the token length is limited; (2) the real system is designed for handling complex tasks, and the base LLMs can hardly plan a correct sub-task order and API-calling order for such tasks; (3) Similar semantics and functionalities among APIs in real systems create challenges for both LLMs and even humans in distinguishing between them. In response, this paper introduces a comprehensive framework aimed at enhancing the Task Planning and Tool Usage (TPTU) abilities of LLM-based agents operating within real-world systems. Our framework comprises three key components designed to address these challenges: (1) the API Retriever selects the most pertinent APIs for the user task among the extensive array available; (2) LLM Finetuner tunes a base LLM so that the finetuned LLM can be more capable for task planning and API calling; (3) the Demo Selector adaptively retrieves different demonstrations related to hard-to-distinguish APIs, which is further used for in-context learning to boost the final performance. We validate our methods using a real-world commercial system as well as an open-sourced academic dataset, and the outcomes clearly showcase the efficacy of each individual component as well as the integrated framework.


title: Word Measures on Wreath Products II

abstract: Every word $w$ in $F_r$, the free group of rank $r$, induces a probability measure (the $w$-measure) on every finite group $G$, by substitution of random $G$-elements in the letters. This measure is determined by its Fourier coefficients: the $w$-expectations $E_w[\chi]$ of the irreducible characters of $G$. For every finite group $G$, every stable character $\chi$ of $G\wr S_n$ (trace of a finitely generated $FI_G$-module), and every word $w\in F_r$, we approximate $E_w[\chi]$ up to an error term of $O(n^{-\pi(w)})$, where $\pi(w)$ is the primitivity rank of $w$. This generalizes previous works by Puder, Hanany, Magee and the author. As an application we show that random Schreier graphs of representation-stable actions of $G\wr S_n$ are close-to-optimal expanders. The paper reveals a surprising relation between stable representation theory of wreath products and not-necessarily connected Stallings core graphs.


title: Discrete approximations of Gaussian smoothing and Gaussian derivatives

abstract: This paper develops an in-depth treatment concerning the problem of approximating the Gaussian smoothing and Gaussian derivative computations in scale-space theory for application on discrete data. With close connections to previous axiomatic treatments of continuous and discrete scale-space theory, we consider three main ways discretizing these scale-space operations in terms of explicit discrete convolutions, based on either (i) sampling the Gaussian kernels and the Gaussian derivative kernels, (ii) locally integrating the Gaussian kernels and the Gaussian derivative kernels over each pixel support region and (iii) basing the scale-space analysis on the discrete analogue of the Gaussian kernel, and then computing derivative approximations by applying small-support central difference operators to the spatially smoothed image data. We study the properties of these three main discretization methods both theoretically and experimentally, and characterize their performance by quantitative measures, including the results they give rise to with respect to the task of scale selection, investigated for four different use cases, and with emphasis on the behaviour at fine scales. The results show that the sampled Gaussian kernels and derivatives as well as the integrated Gaussian kernels and derivatives perform very poorly at very fine scales. At very fine scales, the discrete analogue of the Gaussian kernel with its corresponding discrete derivative approximations performs substantially better. The sampled Gaussian kernel and the sampled Gaussian derivatives do, on the other hand, lead to numerically very good approximations of the corresponding continuous results, when the scale parameter is sufficiently large, in the experiments presented in the paper, when the scale parameter is greater than a value of about 1, in units of the grid spacing.


title: Dynamics of a two-strain epidemic model with waning immunity -- a perturbative approach

abstract: Many infectious diseases are comprised of multiple strains with examples including Influenza, tuberculosis, and Dengue virus. The time evolution of such systems is linked to a complex landscape shaped by interactions between competing strains. Possible long-term dynamics include the extinction of less competitive strains, convergence to multi-strain steady-states, or self-sustained oscillations. This work considers a two-strain epidemic model in which the strains can interact indirectly via the immunity response generated following infections, and in which this immune response wanes with time. In particular, we focus on scenarios where the rate of waning immunity is significantly faster than the rate of demographic turnover. The first key result of this study is the explicit computation of the steady states of the nonlinear system of seven equations. Following this result, we take advantage of the separation of time scales in the problem and use perturbation methods to analyze the stability of the fixed points. In particular, we establish the conditions under which the system gives rise to the coexistence of the two strains and whether coexistence is attained via convergence to an endemic steady-state or via self-sustained oscillations. Our study unveils two parameter regimes of distinct qualitative behavior of the system and characterizes the separatrix between them. Within the first regime, the system gives rise to oscillatory coexistence for all feasible conditions. In the second regime, the system's behavior is governed by a solution to a quadratic equation, potentially resulting in the convergence to a multi-strain endemic equilibrium or the persistence of oscillatory coexistence.


title: GeoSAM: Fine-tuning SAM with Sparse and Dense Visual Prompting for Automated Segmentation of Mobility Infrastructure

abstract: The Segment Anything Model (SAM) has shown impressive performance when applied to natural image segmentation. However, it struggles with geographical images like aerial and satellite imagery, especially when segmenting mobility infrastructure including roads, sidewalks, and crosswalks. This inferior performance stems from the narrow features of these objects, their textures blending into the surroundings, and interference from objects like trees, buildings, vehicles, and pedestrians - all of which can disorient the model to produce inaccurate segmentation maps. To address these challenges, we propose Geographical SAM (GeoSAM), a novel SAM-based framework that implements a fine-tuning strategy using the dense visual prompt from zero-shot learning, and the sparse visual prompt from a pre-trained CNN segmentation model. The proposed GeoSAM outperforms existing approaches for geographical image segmentation, specifically by 26%, 7%, and 17% for road infrastructure, pedestrian infrastructure, and on average, respectively, representing a momentous leap in leveraging foundation models to segment mobility infrastructure including both road and pedestrian infrastructure in geographical images. The source code can be found on this GitHub repository: https://github.com/rafiibnsultan/GeoSAM/tree/main.


title: Equivalence of lattice operators and graph matrices

abstract: We explore the relationship between lattice field theory and graph theory, placing special emphasis on the interplay between Dirac and scalar lattice operators and matrices within the realm of spectral graph theory. Beyond delving into fundamental concepts of spectral graph theory, such as adjacency and Laplacian matrices, we introduce a novel matrix named as "anti-symmetrized adjacency matrix", specifically tailored for cycle digraphs ($T^1$ lattice) and simple directed paths ($B^1$ lattice). The nontrivial relation between graph theory matrices and lattice operators shows that the graph Laplacian matrix mirrors the lattice scalar operator and the Wilson term in lattice fermions, while the anti-symmetrized adjacency matrix, along with its extensions to higher dimensions, are equivalent to naive lattice Dirac operators. Building upon these connections, we provide rigorous proofs for two key assertions: (i) The count of zero-modes in a free lattice scalar operator coincides with the zeroth Betti number of the underlying graph (lattice). (ii) The maximum count of Dirac zero-modes in a free lattice fermion operator is equivalent to the cumulative sum of all Betti numbers when the $D$-dimensional graph results from a cartesian product of cycle digraphs ($T^1$ lattice) and simple directed paths ($B^1$ lattice).


title: Bounds on Representation-Induced Confounding Bias for Treatment Effect Estimation

abstract: State-of-the-art methods for conditional average treatment effect (CATE) estimation make widespread use of representation learning. Here, the idea is to reduce the variance of the low-sample CATE estimation by a (potentially constrained) low-dimensional representation. However, low-dimensional representations can lose information about the observed confounders and thus lead to bias, because of which the validity of representation learning for CATE estimation is typically violated. In this paper, we propose a new, representation-agnostic refutation framework for estimating bounds on the representation-induced confounding bias that comes from dimensionality reduction (or other constraints on the representations) in CATE estimation. First, we establish theoretically under which conditions CATE is non-identifiable given low-dimensional (constrained) representations. Second, as our remedy, we propose a neural refutation framework which performs partial identification of CATE or, equivalently, aims at estimating lower and upper bounds of the representation-induced confounding bias. We demonstrate the effectiveness of our bounds in a series of experiments. In sum, our refutation framework is of direct relevance in practice where the validity of CATE estimation is of importance.


title: An efficient numerical approach to modeling the effects of particle shape on rubble-pile dynamics

abstract: We present an approach for the inclusion of non-spherical constituents in high-resolution N-body discrete element method (DEM) simulations. We use aggregates composed of bonded spheres to model non-spherical components. Though the method may be applied more generally, we detail our implementation in the existing N-body code pkdgrav. It has long been acknowledged that non-spherical grains confer additional shear strength and resistance to flow when compared with spheres. As a result, we expect that rubble-pile asteroids will also exhibit these properties and may behave differently than comparable rubble piles composed of idealized spheres. Since spherical particles avoid some significant technical challenges, most DEM gravity codes have used only spherical particles, or have been confined to relatively low resolutions. We also discuss the work that has gone into improving performance with non-spherical grains, building on pkdgrav's existing leading-edge computational efficiency among DEM gravity codes. This allows for the addition of non-spherical shapes while maintaining the efficiencies afforded by pkdgrav's tree implementation and parallelization. As a test, we simulated the gravitational collapse of 25,000 non-spherical bodies in parallel. In this case, the efficiency improvements allowed for an increase in speed by nearly a factor of three when compared with the naive implementation. Without these enhancements, large runs with non-spherical components would remain prohibitively expensive. Finally, we present the results of several small-scale tests: spinup due to the YORP effect, tidal encounters, and the Brazil-nut Effect. In all cases, we find that the inclusion of non-spherical constituents has a measurable impact on simulation outcomes.


title: Structure and substructure connectivity of folded divide-and-swap cube

abstract: Let $ H $ be a connected subgraph of a graph $ G $. The structure connectivity of $ G $, denoted by $ \kappa(G;H) $, is the minimum number of a set of connected subgraphs in $ G $, whose removal disconnects $ G $ and each element in the set is isomorphic to $ H $. The substructure connectivity of $ G $, denoted by $ \kappa^s(G;H) $, is the minimum number of a set of connected subgraphs in $ G $, whose removal disconnects $ G $ and each element in the set is isomorphic to a connected subgraph of $ H $. In this paper, we determine $ H $-structure connectivity and $ H $-substructure connectivity of folded divide-and-swap cube $ FDSC_n $ for $ H\in\{K_1, K_{1,1}, K_{1,m} (2\leq m \leq d+1) \} $ where $ n=2^d $. We show that $\kappa(FDSC_n;K_1)=\kappa^s(FDSC_n;K_1)=d+2$, $\kappa(FDSC_n;K_{1,1})=\kappa^s(FDSC_n;K_{1,1})=d+1 $ for $ d\geq1 $ and $\kappa(FDSC_n;K_{1,m})=\kappa^s(FDSC_n;K_{1,m})=\lfloor\frac{d}{2}\rfloor+1$ for $d\geq1 $ and $ 2\leq m \leq d+1$.


title: The total number of descents and levels in tensor words and cyclic tensor words

abstract: We obtain an explicit formula for the total number of descents and levels in tensor words and cyclic tensor words of arbitrary dimension. We also establish the exact maximal number of descents in words and use it to obtain an upper bound on the maximal number of descents in tensor words.


title: MoVideo: Motion-Aware Video Generation with Diffusion Models

abstract: While recent years have witnessed great progress on using diffusion models for video generation, most of them are simple extensions of image generation frameworks, which fail to explicitly consider one of the key differences between videos and images, i.e., motion. In this paper, we propose a novel motion-aware video generation (MoVideo) framework that takes motion into consideration from two aspects: video depth and optical flow. The former regulates motion by per-frame object distances and spatial layouts, while the later describes motion by cross-frame correspondences that help in preserving fine details and improving temporal consistency. More specifically, given a key frame that exists or generated from text prompts, we first design a diffusion model with spatio-temporal modules to generate the video depth and the corresponding optical flows. Then, the video is generated in the latent space by another spatio-temporal diffusion model under the guidance of depth, optical flow-based warped latent video and the calculated occlusion mask. Lastly, we use optical flows again to align and refine different frames for better video decoding from the latent space to the pixel space. In experiments, MoVideo achieves state-of-the-art results in both text-to-video and image-to-video generation, showing promising prompt consistency, frame consistency and visual quality.


title: On P\'olya's random walk constants

abstract: A celebrated result in probability theory is that a simple symmetric random walk on the $d$-dimensional lattice $\mathbb{Z}^d$ is recurrent for $d=1,2$ and transient for $d\geq 3$. In this note, we derive a closed-form expression, in terms of the Lauricella function $F_C$, for the return probability for all $d\geq3$. Previously, a closed-form formula had only been available for $d=3$.


title: Analysis of the $Z_b(10650)$ state based on electromagnetic properties

abstract: In this study, the magnetic and quadrupole moments of the $Z_b(10650)$ state are determined using the compact diquark-antidiquark interpolating current through the QCD light-cone sum rule. The values that are obtained as a result of the analysis are as follows: $\mu_{Z_b} = 2.35^{+0.34}_{-0.33}~\mu_N$ and $\mathcal{D}_{Z_b} =(1.82^{+0.35}_{-0.31})\times 10^{-2}~\mbox{fm}^2$. Examining the results obtained, it can be seen that the magnetic moments are large enough to be measured experimentally, while the quadrupole moment is obtained as a small but non-zero value, corresponding to a prolate charge distribution. The magnetic moment is the leading-order response of a bound system to a weak external magnetic field. It therefore provides an excellent platform to probe the internal structures of hadrons governed by the quark-gluon dynamics of QCD.


title: LABCAT: Locally adaptive Bayesian optimization using principal-component-aligned trust regions

abstract: Bayesian optimization (BO) is a popular method for optimizing expensive black-box functions. BO has several well-documented shortcomings, including computational slowdown with longer optimization runs, poor suitability for non-stationary or ill-conditioned objective functions, and poor convergence characteristics. Several algorithms have been proposed that incorporate local strategies, such as trust regions, into BO to mitigate these limitations; however, none address all of them satisfactorily. To address these shortcomings, we propose the LABCAT algorithm, which extends trust-region-based BO by adding a rotation aligning the trust region with the weighted principal components and an adaptive rescaling strategy based on the length-scales of a local Gaussian process surrogate model with automatic relevance determination. Through extensive numerical experiments using a set of synthetic test functions and the well-known COCO benchmarking software, we show that the LABCAT algorithm outperforms several state-of-the-art BO and other black-box optimization algorithms.


title: Matrix manipulations via unitary transformations and ancilla-state measurements

abstract: We propose protocols for calculating inner product, matrix addition and matrix multiplication based on multiqubit Toffoli-type and the simplest one-qubit operations and employ ancilla measurements to remove all garbage of calculations. The depth (runtime) of the addition protocol is $O(1)$ and that of other protocols logarithmically increases with the dimensionality of the considered matrices.


title: Generalized Ricci surfaces

abstract: We consider smooth Riemannian surfaces whose curvature $K$ satisfies the relation $\Delta\log|K-c|=aK+b$ away from points where $K=c$ for some $(a,b,c)\in\mathbb{R}^3$, which we call generalized Ricci surfaces. We prove some isometric immersion theorems allowing points where $K=c$ using properties of log-harmonic functions. For instance, we obtain a characterization of Riemannian surfaces that locally admit minimal isometric immersions, possibly with umbilical points, into a $3$-dimensional Riemannian manifold of constant sectional curvature. We also give an application to convex affine spheres. Finally, we study compact generalized Ricci surfaces: we obtain topological obstructions and construct examples.


title: Portuguese FAQ for Financial Services

abstract: Scarcity of domain-specific data in the Portuguese financial domain has disfavored the development of Natural Language Processing (NLP) applications. To address this limitation, the present study advocates for the utilization of synthetic data generated through data augmentation techniques. The investigation focuses on the augmentation of a dataset sourced from the Central Bank of Brazil FAQ, employing techniques that vary in semantic similarity. Supervised and unsupervised tasks are conducted to evaluate the impact of augmented data on both low and high semantic similarity scenarios. Additionally, the resultant dataset will be publicly disseminated on the Hugging Face Datasets platform, thereby enhancing accessibility and fostering broader engagement within the NLP research community.


title: Approximation Algorithms for Packing Cycles and Paths in Complete Graphs

abstract: Given an edge-weighted (metric/general) complete graph with $n$ vertices, the maximum weight (metric/general) $k$-cycle/path packing problem is to find a set of $\frac{n}{k}$ vertex-disjoint $k$-cycles/paths such that the total weight is maximized. In this paper, we consider approximation algorithms. For metric $k$-cycle packing, we improve the previous approximation ratio from $3/5$ to $7/10$ for $k=5$, and from $7/8\cdot(1-1/k)^2$ for $k>5$ to $(7/8-0.125/k)(1-1/k)$ for constant odd $k>5$ and to $7/8\cdot (1-1/k+\frac{1}{k(k-1)})$ for even $k>5$. For metric $k$-path packing, we improve the approximation ratio from $7/8\cdot (1-1/k)$ to $\frac{27k^2-48k+16}{32k^2-36k-24}$ for even $10\geq k\geq 6$. For the case of $k=4$, we improve the approximation ratio from $3/4$ to $5/6$ for metric 4-cycle packing, from $2/3$ to $3/4$ for general 4-cycle packing, and from $3/4$ to $14/17$ for metric 4-path packing.


title: Stable $(r+1)$-th capillary hypersurfaces

abstract: In this paper, we propose a new definition of stable $(r+1)$-th capillary hypersurfaces from variational perspective for any $1\leq r\leq n-1$. More precisely, we define stable $(r+1)$-th capillary hypersurfaces to be smooth local minimizers of a new energy functional under volume-preserving and contact angle-preserving variations. Using the new concept of the stable $(r+1)$-th capillary hypersurfaces, we generalize the stability results of Souam \cite{Souam} in an Euclidean half-space and Guo-Wang-Xia \cite{GWX} in a horoball in hyperbolic space for capillary hypersurface to $(r+1)$-th capillary hypersurface case.


title: Using Causal Threads to Explain Changes in a Dynamic System

abstract: We explore developing rich semantic models of systems. Specifically, we consider structured causal explanations about state changes in those systems. Essentially, we are developing process-based dynamic knowledge graphs. As an example, we construct a model of the causal threads for geological changes proposed by the Snowball Earth theory. Further, we describe an early prototype of a graphical interface to present the explanations. Unlike statistical approaches to summarization and explanation such as Large Language Models (LLMs), our approach of direct representation can be inspected and verified directly.


title: Self-Distilled Representation Learning for Time Series

abstract: Self-supervised learning for time-series data holds potential similar to that recently unleashed in Natural Language Processing and Computer Vision. While most existing works in this area focus on contrastive learning, we propose a conceptually simple yet powerful non-contrastive approach, based on the data2vec self-distillation framework. The core of our method is a student-teacher scheme that predicts the latent representation of an input time series from masked views of the same time series. This strategy avoids strong modality-specific assumptions and biases typically introduced by the design of contrastive sample pairs. We demonstrate the competitiveness of our approach for classification and forecasting as downstream tasks, comparing with state-of-the-art self-supervised learning methods on the UCR and UEA archives as well as the ETT and Electricity datasets.


title: Error estimate of a quasi-Monte Carlo time-splitting pseudospectral method for nonlinear Schrodinger equation with random potentials

abstract: In this paper, we consider the numerical solution of a nonlinear Schrodinger equation with spatial random potential. The randomly shifted quasi-Monte Carlo (QMC) lattice rule combined with the time-splitting pseudospectral discretization is applied and analyzed. The nonlinearity in the equation induces difficulties in estimating the regularity of the solution in random space. By the technique of weighted Sobolev space, we identify the possible weights and show the existence of QMC that converges optimally at the almost-linear rate without dependence on dimensions. The full error estimate of the scheme is established. We present numerical results to verify the accuracy and investigate the wave propagation.


title: H2 suboptimal containment control of homogeneous and heterogeneous multi-agent systems

abstract: This paper deals with the H2 suboptimal state containment control problem for homogeneous linear multi-agent systems and the H2 suboptimal output containment control problem for heterogeneous linear multi-agent systems. For both problems, given multiple autonomous leaders and a number of followers, we introduce suitable performance outputs and an associated H2 cost functional, respectively. The aim is to design a distributed protocol by dynamic output feedback that achieves state/output containment control while the associated H2 cost is smaller than an a priori given upper bound. To this end, we first show that the H2 suboptimal state/output containment control problem can be equivalently transformed into H2 suboptimal control problems for a set of independent systems. Based on this, design methods are then provided to compute such distributed dynamic output feedback protocols. Simulation examples are provided to illustrate the performance of our proposed protocols.


title: Synchronization rates and limit laws for random dynamical systems

abstract: We study general random dynamical systems of continuous maps on some compact metric space. Assuming a local contraction condition and uniqueness of the stationary measure, we establish probabilistic limit laws such as the central limit theorem, the strong law of large numbers, and the law of the iterated logarithm. Moreover, we study exponential synchronization and synchronization on average. In the particular case of iterated function systems on $\mathbb S^1$, we analyze synchronization rates and describe their large deviations. In the case of $C^{1+\beta}$-diffeomorphisms, these deviations on random orbits are obtained from the large deviations of the expected Lyapunov exponent.


title: Assessment of Transmission-level Fault Impacts on 3-phase and 1-phase Distribution IBR Operation

abstract: The widespread deployment of inverter-based resources (IBRs) renders distribution systems susceptible to transmission-level faults. This paper presents a comprehensive analysis of the impact of transmission-level faults on 3-phase and 1-phase distribution IBR operation. To evaluate distributed IBR tripping across various phases and locations on a distribution feeder, we conduct simulations of both symmetrical and unsymmetrical transmission faults at progressively greater electrical distances on a real-time transmission and distribution (T&D) co-simulation platform. The IBR power-to-load ratios (PLRs) at 50%, 100%, and 300% are considered to emulate low, medium, and high IBR conditions. Our results indicate that, while 1-phase and 2-phase faults typically trigger fewer IBR trips when compared to 3-phase faults, a significant power imbalance arises from the tripping of 1-phase IBRs on the affected phases. The imbalance can result in significant power quality problems and unintended equipment tripping. It may be necessary to design fault-ride-through mechanisms specifically tailored to 1-phase IBRs to help mitigate the power imbalances caused by unbalanced faults.


title: RflyMAD: A Dataset for Multicopter Fault Detection and Health Assessment

abstract: This paper presents an open-source dataset RflyMAD, a Multicopter Abnomal Dataset developed by Reliable Flight Control (Rfly) Group aiming to promote the development of research fields like fault detection and isolation (FDI) or health assessment (HA). The entire 114 GB dataset includes 11 types of faults under 6 flight statuses which are adapted from ADS-33 file to cover more occasions in which the multicopters have different mobility levels when faults occur. In the total 5629 flight cases, the fault time is up to 3283 minutes, and there are 2566 cases for software-in-the-loop (SIL) simulation, 2566 cases for hardware-in-the-loop (HIL) simulation and 497 cases for real flight. As it contains simulation data based on RflySim and real flight data, it is possible to improve the quantity while increasing the data quality. In each case, there are ULog, Telemetry log, Flight information and processed files for researchers to use and check. The RflyMAD dataset could be used as a benchmark for fault diagnosis methods and the support relationship between simulation data and real flight is verified through transfer learning methods. More methods as a baseline will be presented in the future, and RflyMAD will be updated with more data and types. In addition, the dataset and related toolkit can be accessed through https://rfly-openha.github.io/documents/4_resources/dataset.html.


title: Fermi surface topology and electronic transport properties of a chiral crystal NbGe$_2$ with strong electron-phonon interaction

abstract: We report the electronic structures and transport properties of a chiral crystal NbGe$_2$, which is a candidate for a coupled electron-phonon liquid. The electrical resistivity and thermoelectric power of NbGe$_2$ exhibit clear differences compared to those of NbSi2 even though both niobium ditetrelides are isostructural and isoelectronic. We discuss the intriguing transport properties of NbGe$_2$ based on a van Hove-type singularity in the density of states. The analysis of de Haas-van Alphen oscillations measured by the field modulation and magnetic torque methods reveals the detailed shape of the Fermi surface of NbGe$_2$ by comparison with the results of energy band structure calculations using a local density approximation. The electron and hole Fermi surfaces of NbGe$_2$ split into two because of the anti-symmetric spin-orbit interaction. The temperature dependence of quantum oscillations indicates that the effective mass is isotropically enhanced in NbGe$_2$ due to strong electron-phonon interaction.


title: On the Communication Complexity of Decentralized Bilevel Optimization

abstract: Stochastic bilevel optimization finds widespread applications in machine learning, including meta-learning, hyperparameter optimization, and neural architecture search. To extend stochastic bilevel optimization to distributed data, several decentralized stochastic bilevel optimization algorithms have been developed. However, existing methods often suffer from slow convergence rates and high communication costs in heterogeneous settings, limiting their applicability to real-world tasks. To address these issues, we propose two novel decentralized stochastic bilevel gradient descent algorithms based on simultaneous and alternating update strategies. Our algorithms can achieve faster convergence rates and lower communication costs than existing methods. Importantly, our convergence analyses do not rely on strong assumptions regarding heterogeneity. More importantly, our theoretical analysis clearly discloses how the additional communication required for estimating hypergradient under the heterogeneous setting affects the convergence rate. To the best of our knowledge, this is the first time such favorable theoretical results have been achieved with mild assumptions in the heterogeneous setting. Furthermore, we demonstrate how to establish the convergence rate for the alternating update strategy when combined with the variance-reduced gradient. Finally, experimental results confirm the efficacy of our algorithms.


title: A Generative Model for Accelerated Inverse Modelling Using a Novel Embedding for Continuous Variables

abstract: In materials science, the challenge of rapid prototyping materials with desired properties often involves extensive experimentation to find suitable microstructures. Additionally, finding microstructures for given properties is typically an ill-posed problem where multiple solutions may exist. Using generative machine learning models can be a viable solution which also reduces the computational cost. This comes with new challenges because, e.g., a continuous property variable as conditioning input to the model is required. We investigate the shortcomings of an existing method and compare this to a novel embedding strategy for generative models that is based on the binary representation of floating point numbers. This eliminates the need for normalization, preserves information, and creates a versatile embedding space for conditioning the generative model. This technique can be applied to condition a network on any number, to provide fine control over generated microstructure images, thereby contributing to accelerated materials design.


title: Unveiling Deception: Establishing a Taxonomic Framework for Disinformation within Scientific Discourse

abstract: Disinformation spreads among the public and in scientific discourse through the actions of individuals, organizations, and governments that distort scholarly communications, media narratives, and institutional trust. This taxonomy introduces a structured framework and specialized set of definitions to elucidate the key participants, platforms, and strategies employed in the propagation of disinformation. Enhanced comprehension of the mechanisms and pathways of scientific disinformation equips journalists and policymakers with the tools necessary to more effectively recognize and address these issues. The authors developed this taxonomy of disinformation through a multi-faceted approach, encompassing a literature review, expert review, and case study analysis. The literature review revealed a scarcity of taxonomical models amidst prevalent algorithmic detection studies. Subsequently, an expert review process refined our taxonomy through collaborative analysis of twenty-two cases of identified disinformation, categorized by their methods, motives, and impacts. Finally, we validated and fine-tuned our taxonomy through detailed case studies of twelve diverse disinformation instances, assessing the taxonomy's effectiveness in capturing the essential characteristics of each case and making necessary adjustments to ensure its relevance and accuracy in real-world applications.


title: Celestial conformal blocks of massless scalars and analytic continuation of the Appell function $F_1$

abstract: In celestial conformal field theory (CCFT), the 4d massless scalars are represented by 2d conformal operators with conformal dimensions $h=\bar{h}=(1+i\lambda)/2$. The Mellin transform of 4d massless scalar amplitudes gives the conformal correlators of CCFT. We study the conformal block decomposition of these celestial correlators in CCFT and obtain the explicit blocks. This conformal block decomposition is highly nontrivial, even for the simplest 4d massless scalar amplitude. We use the analytic continuation of the Appell hypergeometric function $F_1$ and the method of monodromy projection of conformal blocks, to achieve this block decomposition. This procedure is consistent with the crossing symmetry, in both the correlator-level and each explicit block-level. We also investigate its behavior in the conformal soft limit and find that the Appell hypergeometric function $F_1$ does not reduce to the Gauss hypergeometric function. This is different from the block decomposition of celestial gluons we studied before, where the Appell hypergeometric function $F_1$ reduces to the Gauss hypergeometric function. This difference comes from the shift of conformal dimensions and is the reason why we adopt the new method here for the block decomposition of celestial massless scalars.


title: Multicritical dissipative phase transitions in the anisotropic open quantum Rabi model

abstract: We investigate the nonequilibrium steady state of the anisotropic open quantum Rabi model, which exhibits first-order and second-order dissipative phase transitions upon varying the degree of anisotropy between the coupling strengths of rotating and counterrotating terms. Using both semiclassical and quantum approaches, we find a rich phase diagram resulting from the interplay between the anisotropy and the dissipation. First, there exists a bistable phase where both the normal and superradiant phases are stable. Second, there are multicritical points where the phase boundaries for the first- and second-order phase transitions meet. We show that a new set of critical exponents governs the scaling of the multicritical points. Finally, we discuss the feasibility of observing the multicritical transitions and bistability using a pair of trapped ions where the anisotropy can be tuned by the controlling the intensity of the Raman transitions. Our study enlarges the scope of critical phenomena that may occur in finite-component quantum systems, which could be useful for the applications in the critical quantum sensing.


title: Large-scale Mixed Traffic Control Using Dynamic Vehicle Routing and Privacy-Preserving Crowdsourcing

abstract: Controlling and coordinating urban traffic flow through robot vehicles is emerging as a novel transportation paradigm for the future. While this approach garners growing attention from researchers and practitioners, effectively managing and coordinating large-scale mixed traffic remains a challenge. We introduce an effective framework for large-scale mixed traffic control via privacy-preserving crowdsourcing and dynamic vehicle routing. Our framework consists of three modules: a privacy-protecting crowdsensing method, a graph propagation-based traffic forecasting method, and a privacy-preserving route selection mechanism. We evaluate our framework using a real-world road network. The results show that our framework accurately forecasts traffic flow, efficiently mitigates network-wide RV shortage issue, and coordinates large-scale mixed traffic. Compared to other baseline methods, our framework not only reduces the RV shortage issue up to 69.4% but also reduces the average waiting time of all vehicles in the network up to 27%.


title: p-adaptive discontinuous Galerkin method for the shallow water equations on heterogeneous computing architectures

abstract: Heterogeneous computing and exploiting integrated CPU-GPU architectures has become a clear current trend since the flattening of Moore's Law. In this work, we propose a numerical and algorithmic re-design of a p-adaptive quadrature-free discontinuous Galerkin method (DG) for the shallow water equations (SWE). Our new approach separates the computations of the non-adaptive (lower-order) and adaptive (higher-order) parts of the discretization form each other. Thereby, we can overlap computations of the lower-order and the higher-order DG solution components. Furthermore, we investigate execution times of main computational kernels and use automatic code generation to optimize their distribution between the CPU and GPU. Several setups, including a prototype of a tsunami simulation in a tide-driven flow scenario, are investigated, and the results show that significant performance improvements can be achieved in suitable setups.


title: Coverage-Validity-Aware Algorithmic Recourse

abstract: Algorithmic recourse emerges as a prominent technique to promote the explainability, transparency and hence ethics of machine learning models. Existing algorithmic recourse approaches often assume an invariant predictive model; however, the predictive model is usually updated upon the arrival of new data. Thus, a recourse that is valid respective to the present model may become invalid for the future model. To resolve this issue, we propose a novel framework to generate a model-agnostic recourse that exhibits robustness to model shifts. Our framework first builds a coverage-validity-aware linear surrogate of the nonlinear (black-box) model; then, the recourse is generated with respect to the linear surrogate. We establish a theoretical connection between our coverage-validity-aware linear surrogate and the minimax probability machines (MPM). We then prove that by prescribing different covariance robustness, the proposed framework recovers popular regularizations for MPM, including the $\ell_2$-regularization and class-reweighting. Furthermore, we show that our surrogate pushes the approximate hyperplane intuitively, facilitating not only robust but also interpretable recourses. The numerical results demonstrate the usefulness and robustness of our framework.


title: An Alternative to Regulation: The Case for Public AI

abstract: Can governments build AI? In this paper, we describe an ongoing effort to develop ``public AI'' -- publicly accessible AI models funded, provisioned, and governed by governments or other public bodies. Public AI presents both an alternative and a complement to standard regulatory approaches to AI, but it also suggests new technical and policy challenges. We present a roadmap for how the ML research community can help shape this initiative and support its implementation, and how public AI can complement other responsible AI initiatives.


title: Scaling Law of Large Sequential Recommendation Models

abstract: Scaling of neural networks has recently shown great potential to improve the model capacity in various fields. Specifically, model performance has a power-law relationship with model size or data size, which provides important guidance for the development of large-scale models. However, there is still limited understanding on the scaling effect of user behavior models in recommender systems, where the unique data characteristics (e.g. data scarcity and sparsity) pose new challenges to explore the scaling effect in recommendation tasks. In this work, we focus on investigating the scaling laws in large sequential recommendation models. Specially, we consider a pure ID-based task formulation, where the interaction history of a user is formatted as a chronological sequence of item IDs. We don't incorporate any side information (e.g. item text), because we would like to explore how scaling law holds from the perspective of user behavior. With specially improved strategies, we scale up the model size to 0.8B parameters, making it feasible to explore the scaling effect in a diverse range of model sizes. As the major findings, we empirically show that scaling law still holds for these trained models, even in data-constrained scenarios. We then fit the curve for scaling law, and successfully predict the test loss of the two largest tested model scales. Furthermore, we examine the performance advantage of scaling effect on five challenging recommendation tasks, considering the unique issues (e.g. cold start, robustness, long-term preference) in recommender systems. We find that scaling up the model size can greatly boost the performance on these challenging tasks, which again verifies the benefits of large recommendation models.


title: Bell-INGARCH Model

abstract: Integer-valued time series exist widely in economics, finance, biology, computer science, medicine, insurance, and many other fields. In recent years, many types of models have been proposed to model integer-valued time series data, in which the integer autoregressive model and integer-valued GARCH model are the most representative. Although there have been many results of integer-valued time series data, the parameters of integer-valued time series model structure are more complicated. This paper is dedicated to proposing a new simple integer-valued GARCH model. First, the Bell integer-valued GARCH model is given based on Bell distribution. Then, the conditional maximum likelihood estimation method is used to obtain the estimators of parameters. Later, numerical simulations confirm the finite sample properties of the estimation of unknown parameters. Finally, the model is applied in the two real examples. Compared with the existing models, the proposed model is more simple and applicable.


title: Label-Synchronous Neural Transducer for Adaptable Online E2E Speech Recognition

abstract: Although end-to-end (E2E) automatic speech recognition (ASR) has shown state-of-the-art recognition accuracy, it tends to be implicitly biased towards the training data distribution which can degrade generalisation. This paper proposes a label-synchronous neural transducer (LS-Transducer), which provides a natural approach to domain adaptation based on text-only data. The LS-Transducer extracts a label-level encoder representation before combining it with the prediction network output. Since blank tokens are no longer needed, the prediction network performs as a standard language model, which can be easily adapted using text-only data. An Auto-regressive Integrate-and-Fire (AIF) mechanism is proposed to generate the label-level encoder representation while retaining low latency operation that can be used for streaming. In addition, a streaming joint decoding method is designed to improve ASR accuracy while retaining synchronisation with AIF. Experiments show that compared to standard neural transducers, the proposed LS-Transducer gave a 12.9% relative WER reduction (WERR) for intra-domain LibriSpeech data, as well as 21.4% and 24.6% relative WERRs on cross-domain TED-LIUM 2 and AESRC2020 data with an adapted prediction network.


title: Scale-aware competition network for palmprint recognition

abstract: Palmprint biometrics garner heightened attention in palm-scanning payment and social security due to their distinctive attributes. However, prevailing methodologies singularly prioritize texture orientation, neglecting the significant texture scale dimension. We design an innovative network for concurrently extracting intra-scale and inter-scale features to redress this limitation. This paper proposes a scale-aware competitive network (SAC-Net), which includes the Inner-Scale Competition Module (ISCM) and the Across-Scale Competition Module (ASCM) to capture texture characteristics related to orientation and scale. ISCM efficiently integrates learnable Gabor filters and a self-attention mechanism to extract rich orientation data and discern textures with long-range discriminative properties. Subsequently, ASCM leverages a competitive strategy across various scales to effectively encapsulate the competitive texture scale elements. By synergizing ISCM and ASCM, our method adeptly characterizes palmprint features. Rigorous experimentation across three benchmark datasets unequivocally demonstrates our proposed approach's exceptional recognition performance and resilience relative to state-of-the-art alternatives.


title: The Gamma Ray Origin in RXJ0852.0-4622 Quantifying the Hadronic and Leptonic Components: Further Evidence for the Cosmic Ray Acceleration in Young Shell-type SNRs

abstract: Fukui et al. (2021) quantified the hadronic and leptonic gamma rays in the young TeV gamma ray shell-type supernova remnant (SNR) RXJ1713.7-3946 (RXJ1713), and demonstrated that the gamma rays are a combination of the hadronic and leptonic gamma ray components with a ratio of $\sim 6:4$ in gamma ray counts $N_\mathrm{g}$. This discovery, which adopted a new methodology of multiple-linear gamma-ray decomposition, was the first quantification of the two gamma ray components. In the present work, we applied the same methodology to another TeV gamma ray shell-type SNR RX~J0852.0$-$4622 (RXJ0852) in the 3D space characterized by [the interstellar proton column density $N_{\mathrm{p}}$]-[the nonthermal X-ray count $N_{\mathrm{x}}$]-[$N_{\mathrm{g}}$], and quantified the hadronic and leptonic gamma ray components to have a ratio of $\sim 5:5$ in $N_{\mathrm{g}}$. The present work adopted fitting of two/three flat planes in the 3D space instead of a single flat plane, which allowed to suppress fitting errors. The quantification indicates that the hadronic and leptonic gamma rays are in the same order of magnitude in these two core-collapse SNRs, verifying the significant hadronic gamma ray components. We argue that the target interstellar protons, in particular their spatial distribution, are essential in any attempts to identify type of particles responsible for the gamma-ray emission. The present results confirm that the CR energy $\lesssim 100$\,TeV is compatible with a scheme that SNRs are the dominant source of these Galactic CRs.


title: Charged particles in the background of the Kiselev solution in power-Maxwell electrodynamics

abstract: In this work we analyze the motion of charged particles in the background of the Kiselev geometry, which is considered here as an exact solution in the context of power-Maxwell electrodynamics. As it is well known, one can use either an electric ansatz or a magnetic one for the nonlinear electromagnetic field. We study the motion of an electrically charged particle for an electrically charged black hole and also for a magnetically charged black hole. In the second case the motion is restricted to Poincar\'e cones of various angles, as expected.


title: Diamagnetic susceptibility of neon and argon including leading relativistic effects

abstract: We report theoretical calculations of the static diamagnetic susceptibility, $\chi_0$, of neon and argon atoms. The calculations were performed using a hierarchy coupled-cluster methods combined with application of both the Gaussian and Slater orbital basis sets. We included the complete relativistic correction of order of $\alpha^4$, where $\alpha$ is the fine structure constant, and obtained an estimate of the quantum electrodynamics (QED) contributions. The finite nuclear mass and size corrections were also considered but are found to be small. The final results are $\chi_0=-8.4786(7)\cdot10^{-5}$ $a_0^3$ for the neon and $\chi_0=-22.9545(32)\cdot10^{-5}$ $a_0^3$ for the argon atom, where $a_0$ is the Bohr radius. The uncertainties in the last digits, shown in the parentheses, are primarily due to the errors in the non-relativistic electronic wavefunction, as well as due to the neglected quantum electrodynamics corrections.


title: Quasi-invariant theorem on the Gaussian path space

abstract: In this article, we will first introduce a class of Gaussian processes, and prove the quasi-invariant theorem with respect to the Gaussian Wiener measure, which is the law of the associated Gaussian process. In particular, it includes the case of the fractional Brownian motion. As applications, we will establish the integration by parts formula and Bismut-Elworthy-Li formula on the Gaussian path space, and by which some logarithmic Sobolev inequalities will be presented. Moreover, we will also provides some applications in the field of financial mathematics.


title: On the co-rotation of Milky Way satellites: LMC-mass satellites induce apparent motions in outer halo tracers

abstract: Understanding the physical mechanism behind the formation of a co-rotating thin plane of satellite galaxies, like the one observed around the Milky Way (MW), has been challenging. The perturbations induced by a massive satellite galaxy, like the Large Magellanic Cloud (LMC) provide valuable insight into this problem. The LMC induces an apparent co-rotating motion in the outer halo by displacing the inner regions of the halo with respect to the outer halo. Using the Latte suite of FIRE-2 cosmological simulations of MW-mass galaxies, we confirm that the apparent motion of the outer halo induced by the infall of a massive satellite changes the observed distribution of orbital poles of outer-halo tracers, including satellites. We quantify the changes in the distribution of orbital poles using the two-point angular correlation function and find that all satellites induce changes. However, the most massive satellites with pericentric passages between 30-100kpc induce the largest changes. The best LMC-like satellite analog shows the largest change in orbital pole distribution. The dispersion of orbital poles decreases by 20{\deg} during the first two pericentric passages. Even when excluding the satellites brought in with the LMC-like satellite, there is clustering of orbital poles. These results suggest that in the MW, the recent pericentric passage of the LMC should have changed the observed distribution of orbital poles of all other satellites. Therefore, studies of kinematically-coherent planes of satellites that seek to place the MW in a cosmological context should account for the existence of a massive satellite like the LMC.


title: Locally compact modules over abelian groups and compactly generated metabelian groups

abstract: We perform a general study of the structure of locally compact modules over compactly generated abelian groups. We obtain a devissage result for such modules of the form "compact-by-sheer-by-discrete", and then study more specifically the sheer part. The main typical example of a sheer module is a polycontractable module, i.e., a finite direct product of modules, each of which is contracted by some group element. We show that every sheer module has a "large" polycontractable submodule, in a suitable sense. We apply this to the study of compactly generated metabelian groups. For instance, we prove that they always have a maximal compact normal subgroup, and we extend the Bieri-Strebel characterization of compactly presentable metabelian groups from the discrete case to this more general setting.


title: Quantum flux operators in higher spin theories

abstract: We construct Carrollian higher spin field theories by reducing the bosonic Fronsdal theories in flat spacetime to future null infinity. We extend the Poincar\'e fluxes to quantum flux operators which generate Carrollian diffeomorphism, namely supertranslation and superrotation. These flux operators form a closed symmetry algebra once including a helicity flux operator which follows from higher spin super-duality transformation. The super-duality transformation is an angle-dependent transformation at future null infinity which generalizes the usual electro-magnetic duality transformation. The results agree with the lower spin cases when restricting to $s=0,1,2$.


title: Symmetry-invariant quantum machine learning force fields

abstract: Machine learning techniques are essential tools to compute efficient, yet accurate, force fields for atomistic simulations. This approach has recently been extended to incorporate quantum computational methods, making use of variational quantum learning models to predict potential energy surfaces and atomic forces from ab initio training data. However, the trainability and scalability of such models are still limited, due to both theoretical and practical barriers. Inspired by recent developments in geometric classical and quantum machine learning, here we design quantum neural networks that explicitly incorporate, as a data-inspired prior, an extensive set of physically relevant symmetries. We find that our invariant quantum learning models outperform their more generic counterparts on individual molecules of growing complexity. Furthermore, we study a water dimer as a minimal example of a system with multiple components, showcasing the versatility of our proposed approach and opening the way towards larger simulations. Our results suggest that molecular force fields generation can significantly profit from leveraging the framework of geometric quantum machine learning, and that chemical systems represent, in fact, an interesting and rich playground for the development and application of advanced quantum machine learning tools.

categories: quant-ph cs.LG physics.chem-ph physics.comp-ph


title: Encoding Performance Data in MEI with the Automatic Music Performance Analysis and Comparison Toolkit (AMPACT)

abstract: This paper presents a new method of encoding performance data in MEI using the recently added \texttt{<extData>} element. Performance data was extracted using the Automatic Music Performance Analysis and Comparison Toolkit (AMPACT) and encoded as a JSON object within an \texttt{<extData>} element linked to a specific musical note. A set of pop music vocals has was encoded to demonstrate both the range of descriptors that can be encoded in <extData> and how AMPACT can be used for extracting performance data in the absence of a fully specified musical score.


title: Local environment-based machine learning for molecular adsorption energy prediction

abstract: Most machine learning (ML) models in Materials Science are developed by global geometric features, often falling short in describing localized characteristics, like molecular adsorption on materials. In this study, we introduce a local environment framework that extracts local features from crystal structures to portray the environment surrounding specific adsorption sites. Upon OC20 database (~20,000 3D entries), we apply our local environment framework on several ML models, such as random forest, convolutional neural network, and graph neural network. It is found that our framework achieves remarkable prediction accuracy in predicting molecular adsorption energy, significantly outperforming other examined global-environment-based models. Moreover, the employment of this framework reduces data requirements and augments computational speed, specifically for deep learning algorithms. Finally, we directly apply our Local Environment ResNet (LERN) on a small 2DMatPedia database (~2,000 2D entries), which also achieves highly accurate prediction, demonstrating the model transferability and remarkable data efficiency. Overall, the prediction accuracy, data-utilization efficiency, and transferability of our local-environment-based ML framework hold a promising high applicability across a broad molecular adsorption field, such as catalysis and sensor technologies.


title: Circuit complexity of quantum access models for encoding classical data

abstract: Classical data encoding is usually treated as a black-box in the oracle-based quantum algorithms. On the other hand, their constructions are crucial for practical algorithm implementations. Here, we open the black-boxes of data encoding and study the Clifford$+T$ complexity of constructing some typical quantum access models. For general matrices, we show that both sparse-access input models and block-encoding require nearly linear circuit complexities relative to the matrix dimension, even if matrices are sparse. We also gives construction protocols achieving near-optimal gate complexities. On the other hand, the construction becomes efficient with respect to the data qubit when the matrix is the linear combination polynomial terms of efficient unitaries. As a typical example, we propose improved block encoding when these unitaries are Pauli strings. Our protocols are built upon improved quantum state preparation and a selective oracle for Pauli strings, which hold independent value. Our access model constructions offer considerable flexibility, allowing for tunable ancillary qubit number and offers corresponding space-time trade-offs.


title: Ambiguity aversion as a route to randomness in a duopoly game

abstract: The global dynamics is investigated for a duopoly game where the perfect foresight hypothesis is relaxed and firms are worst-case maximizers. Overlooking the degree of product substitutability as well as the sensitivity of price to quantity, the unique and globally stable Cournot-Nash equilibrium of the complete-information duopoly game, loses stability when firms are not aware if they are playing a duopoly game, as it is, or an oligopoly game with more than two competitors. This finding resembles Theocharis' condition for the stability of the Cournot-Nash equilibrium in oligopolies without uncertainty. As opposed to complete-information oligopoly games, coexisting attractors, disconnected basins of attractions and chaotic dynamics emerge when the Cournot-Nash equilibrium loses stability. This difference in the global dynamics is due to the nonlinearities introduced by the worst-case approach to uncertainty, which mirror in bimodal best-reply functions. Conducted with techniques that require a symmetric setting of the game, the investigation of the dynamics reveals that a chaotic regime prevents firms from being ambiguity averse, that is, firms are worst-case maximizers only in the quantity-expectation space. Therefore, chaotic dynamics are the result and at the same time the source of profit uncertainty.


title: Evidential Uncertainty Quantification: A Variance-Based Perspective

abstract: Uncertainty quantification of deep neural networks has become an active field of research and plays a crucial role in various downstream tasks such as active learning. Recent advances in evidential deep learning shed light on the direct quantification of aleatoric and epistemic uncertainties with a single forward pass of the model. Most traditional approaches adopt an entropy-based method to derive evidential uncertainty in classification, quantifying uncertainty at the sample level. However, the variance-based method that has been widely applied in regression problems is seldom used in the classification setting. In this work, we adapt the variance-based approach from regression to classification, quantifying classification uncertainty at the class level. The variance decomposition technique in regression is extended to class covariance decomposition in classification based on the law of total covariance, and the class correlation is also derived from the covariance. Experiments on cross-domain datasets are conducted to illustrate that the variance-based approach not only results in similar accuracy as the entropy-based one in active domain adaptation but also brings information about class-wise uncertainties as well as between-class correlations. The code is available at https://github.com/KerryDRX/EvidentialADA. This alternative means of evidential uncertainty quantification will give researchers more options when class uncertainties and correlations are important in their applications.


title: Self-Supervised Pretraining for Heterogeneous Hypergraph Neural Networks

abstract: Recently, pretraining methods for the Graph Neural Networks (GNNs) have been successful at learning effective representations from unlabeled graph data. However, most of these methods rely on pairwise relations in the graph and do not capture the underling higher-order relations between entities. Hypergraphs are versatile and expressive structures that can effectively model higher-order relationships among entities in the data. Despite the efforts to adapt GNNs to hypergraphs (HyperGNN), there are currently no fully self-supervised pretraining methods for HyperGNN on heterogeneous hypergraphs. In this paper, we present SPHH, a novel self-supervised pretraining framework for heterogeneous HyperGNNs. Our method is able to effectively capture higher-order relations among entities in the data in a self-supervised manner. SPHH is consist of two self-supervised pretraining tasks that aim to simultaneously learn both local and global representations of the entities in the hypergraph by using informative representations derived from the hypergraph structure. Overall, our work presents a significant advancement in the field of self-supervised pretraining of HyperGNNs, and has the potential to improve the performance of various graph-based downstream tasks such as node classification and link prediction tasks which are mapped to hypergraph configuration. Our experiments on two real-world benchmarks using four different HyperGNN models show that our proposed SPHH framework consistently outperforms state-of-the-art baselines in various downstream tasks. The results demonstrate that SPHH is able to improve the performance of various HyperGNN models in various downstream tasks, regardless of their architecture or complexity, which highlights the robustness of our framework.


title: Optimal Locally Private Nonparametric Classification with Public Data

abstract: In this work, we investigate the problem of public data assisted non-interactive Local Differentially Private (LDP) learning with a focus on non-parametric classification. Under the posterior drift assumption, we for the first time derive the mini-max optimal convergence rate with LDP constraint. Then, we present a novel approach, the locally differentially private classification tree, which attains the mini-max optimal convergence rate. Furthermore, we design a data-driven pruning procedure that avoids parameter tuning and provides a fast converging estimator. Comprehensive experiments conducted on synthetic and real data sets show the superior performance of our proposed methods. Both our theoretical and experimental findings demonstrate the effectiveness of public data compared to private data, which leads to practical suggestions for prioritizing non-private data collection.


title: Higher-order moments of the elliptic flow distribution in PbPb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV

abstract: The hydrodynamic flow-like behavior of charged hadrons in high-energy lead-lead collisions is studied through multiparticle correlations. The elliptic anisotropy values based on different orders of multiparticle cumulants, $v_{2}\{2k\}$, are measured up to the tenth order ($k$ = 5) as functions of the collision centrality at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The data were recorded by the CMS experiment at the LHC and correspond to an integrated luminosity of 0.607 nb$^{-1}$. A hierarchy is observed between the coefficients, with $v_{2}\{2\} > v_{2}\{4\} \gtrsim v_{2}\{6\} \gtrsim v_{2}\{8\} \gtrsim v_{2}\{10\}$. Based on these results, centrality-dependent moments for the fluctuation-driven event-by-event $v_{2}$ distribution are determined, including the skewness, kurtosis and, for the first time, superskewness. Assuming a hydrodynamic expansion of the produced medium, these moments directly probe the initial-state geometry in high-energy nucleus-nucleus collisions.


title: SOccDPT: Semi-Supervised 3D Semantic Occupancy from Dense Prediction Transformers trained under memory constraints

abstract: We present SOccDPT, a memory-efficient approach for 3D semantic occupancy prediction from monocular image input using dense prediction transformers. To address the limitations of existing methods trained on structured traffic datasets, we train our model on unstructured datasets including the Indian Driving Dataset and Bengaluru Driving Dataset. Our semi-supervised training pipeline allows SOccDPT to learn from datasets with limited labels by reducing the requirement for manual labelling by substituting it with pseudo-ground truth labels to produce our Bengaluru Semantic Occupancy Dataset. This broader training enhances our model's ability to handle unstructured traffic scenarios effectively. To overcome memory limitations during training, we introduce patch-wise training where we select a subset of parameters to train each epoch, reducing memory usage during auto-grad graph construction. In the context of unstructured traffic and memory-constrained training and inference, SOccDPT outperforms existing disparity estimation approaches as shown by the RMSE score of 9.1473, achieves a semantic segmentation IoU score of 46.02% and operates at a competitive frequency of 69.47 Hz. We make our code and semantic occupancy dataset public.


title: Dynamic System Stability Verification Using Numerical Simulator

abstract: There are recent shifts in demand for design controllers from simplified to complex model-based. Although simplification approaches are successful in many areas of engineering control systems, high-fidelity simulation-based control design, for example, reinforcement learning, has been rising in robotics areas. On the other hand, the lack of assurances about the stability and robustness of simulation-based control design restricts its applications to safety-critical systems. We develop computational methods to verify the stability and robustness of safety-critical systems. By extending the inverse Lyapunov theorem, we present a practical method to compute the constants required to check the exponential stability conditions of dynamic systems implemented in a numerical simulator. It is shown that the norm-bound of the propagated states is a function of the numerical integration steps, where the numerical simulator may include discontinuous jumps of states. The energy bounds for the transition states are obtained based on the exponential stability assumption of the inverse Lyapunov theorem. Finally, a finite sampling algorithm provides the deterministic stability guarantee for the continuous state space.


title: Ion Dynamics Across a Low Mach Number Bow Shock

abstract: A thorough understanding of collisionless shocks requires knowledge of how different ion species are accelerated across the shock. We investigate a bow shock crossing using the Magnetospheric Multiscale spacecraft after a coronal mass ejection crossed Earth, which led to solar wind consisting of protons, alpha particles, and singly charge helium ions. The low Mach number of the bow shock enabled the ions to be distinguished upstream and sometimes downstream of the shock. Some of the protons are specularly reflected and produce quasi-periodic fine structures in the velocity distribution functions downstream of the shock. Heavier ions are shown to transit the shock without reflection. However, the gyromotion of the heavier ions partially obscures the fine structure of proton distributions. Additionally, the calculated proton moments are unreliable when the different ion species are not distinguished by the particle detector. The need to high time-resolution mass-resolving ion detectors when investigating collisionless shocks is discussed.


title: Rapidity scan approach for net-baryon cumulants with a statistical thermal model

abstract: Utilizing rapidity-dependent measurements to map the QCD phase diagram provides a complementary approach to traditional beam-energy-dependent measurements around midrapidity. The changing nature of thermodynamic properties of QCD matter along the beam axis in heavy-ion collisions at low collision energies both motivates and poses challenges for this method. In this study, we derive the analytical cumulant-generating function for subsystems within distinct rapidity windows, while accounting for global net-baryon charge conservation of the full system. Rapidity-dependent net-baryon cumulants are then calculated for a system exhibiting inhomogeneity along the beam axis, and their sensitivity to finite acceptances through changing rapidity bin widths is explored. We highlight the nontrivial behaviors exhibited by these cumulants, underscoring their importance in establishing a noncritical baseline for interpreting net-proton cumulants in the search for the QCD critical point. Finally, we discuss the implications of the rapidity scan for mapping the QCD phase diagram within the current context.


title: ML-LMCL: Mutual Learning and Large-Margin Contrastive Learning for Improving ASR Robustness in Spoken Language Understanding

abstract: Spoken language understanding (SLU) is a fundamental task in the task-oriented dialogue systems. However, the inevitable errors from automatic speech recognition (ASR) usually impair the understanding performance and lead to error propagation. Although there are some attempts to address this problem through contrastive learning, they (1) treat clean manual transcripts and ASR transcripts equally without discrimination in fine-tuning; (2) neglect the fact that the semantically similar pairs are still pushed away when applying contrastive learning; (3) suffer from the problem of Kullback-Leibler (KL) vanishing. In this paper, we propose Mutual Learning and Large-Margin Contrastive Learning (ML-LMCL), a novel framework for improving ASR robustness in SLU. Specifically, in fine-tuning, we apply mutual learning and train two SLU models on the manual transcripts and the ASR transcripts, respectively, aiming to iteratively share knowledge between these two models. We also introduce a distance polarization regularizer to avoid pushing away the intra-cluster pairs as much as possible. Moreover, we use a cyclical annealing schedule to mitigate KL vanishing issue. Experiments on three datasets show that ML-LMCL outperforms existing models and achieves new state-of-the-art performance.


title: Post-Giant Impact Planetesimals Sustaining Extreme Debris Disks

abstract: Extreme debris disks can show short term behaviour through the evolution and clearing of small grains produced in giant impacts, and potentially a longer period of variability caused by a planetesimal population formed from giant impact ejecta. In this paper, we present results of numerical simulations to explain how a planetesimal populated disk can supply an observed extreme debris disk with small grains. We simulated a sample of giant impacts from which we form a planetesimal population. We then use the $N$-body code {\sc Rebound} to evolve the planetesimals spatially and collisionally. We adopt a simplistic collision criteria in which we define destructive collisions to be between planetesimals with a mutual impact velocity that exceeds two times the catastrophic disruption threshold, $V^*$. We find that for some configurations, a planetesimal populated disk can produce a substantial amount of dust to sustain an observable disk. The semi-major axis at which the giant impact occurs changes the mass added to the observed disk substantially while the orientation of the impact has less of an effect. We determine how the collision rate at the collision point changes over time and show that changes in semi-major axis and orientation only change the initial collision rate of the disk. Collision rates across all disks evolve at a similar rate.


title: The fluid mechanics of splat painting

abstract: In splat painting, a collection of liquid droplets is projected onto the substrate by imposing a controlled acceleration to a paint-loaded brush. To unravel the physical phenomena at play in this artistic technique, we perform a series of experiments where the amount of expelled liquid and the resulting patterns on the substrate are systematically characterized as a function of the liquid viscosity and brush acceleration. Experimental trends and orders of magnitude are rationalized by simple physical models, revealing the existence of an inertia-dominated flow in the anisotropic, porous tip of the brush. We argue that splat painting artists intuitively tune their parameters to work in this regime, which may also play a role in other pulsed flows, like violent expiratory events or sudden geophysical processes.


title: Inspecting Explainability of Transformer Models with Additional Statistical Information

abstract: Transformer becomes more popular in the vision domain in recent years so there is a need for finding an effective way to interpret the Transformer model by visualizing it. In recent work, Chefer et al. can visualize the Transformer on vision and multi-modal tasks effectively by combining attention layers to show the importance of each image patch. However, when applying to other variants of Transformer such as the Swin Transformer, this method can not focus on the predicted object. Our method, by considering the statistics of tokens in layer normalization layers, shows a great ability to interpret the explainability of Swin Transformer and ViT.


title: On the zero set of holomorphic sectional curvature

abstract: A notable example due to Heier, Lu, Wong, and Zheng shows that there exist compact complex K\"ahler manifolds with ample canonical line bundle such that the holomorphic sectional curvature is negative semi-definite and vanishes along high-dimensional linear subspaces in every tangent space. The main result of this note is an upper bound for the dimensions of these subspaces. Due to the holomorphic sectional curvature being a real-valued bihomogeneous polynomial of bidegree $(2,2)$ on every tangent space, the proof is based on making a connection with the work of D'Angelo on complex subvarieties of real algebraic varieties and the decomposition of polynomials into differences of squares. Our bound involves an invariant that we call the holomorphic sectional curvature square decomposition length, and our arguments work as long as the holomorphic sectional curvature is semi-definite, be it negative or positive.


title: Equilibrate Parametrization: Optimal Metric Selection with Provable One-iteration Convergence for $ l_1 $-minimization

abstract: Incorporating a non-Euclidean variable metric to first-order algorithms is known to bring enhancement. However, due to the lack of an optimal choice, such an enhancement appears significantly underestimated. In this work, we establish a metric selection principle via optimizing a convergence rate upper-bound. For general l1-minimization, we propose an optimal metric choice with closed-form expressions guaranteed. Equipping such a variable metric, we prove that the optimal solution to the l1 problem will be obtained via a one-time proximal operator evaluation. Our technique applies to a large class of fixed-point algorithms, particularly the ADMM, which is popular, general, and requires minimum assumptions. The key to our success is the employment of an unscaled/equilibrate upper-bound. We show that there exists an implicit scaling that poses a hidden obstacle to optimizing parameters. This turns out to be a fundamental issue induced by the classical parametrization. We note that the conventional way always associates the parameter to the range of a function/operator. This turns out not a natural way, causing certain symmetry losses, definition inconsistencies, and unnecessary complications, with the well-known Moreau identity being the best example. We propose equilibrate parametrization, which associates the parameter to the domain of a function, and to both the domain and range of a monotone operator. A series of powerful results are obtained owing to the new parametrization. Quite remarkably, the preconditioning technique can be shown as equivalent to the metric selection issue.


title: Algorithms for Gromov-Witten Invariants of Elliptic Curves

abstract: We present an enhanced algorithm for exploring mirror symmetry for elliptic curves through the correspondence of algebraic and tropical geometry, focusing on Gromov-Witten invariants of elliptic curves and, in particular, Hurwitz numbers. We present a new highly efficient algorithm for computing generating series for these numbers. We have implemented the algorithm both using Singular and OSCAR. The implementations outperform by far the current method provided in Singular. The OSCAR implementation, benefiting in particular from just-in-time compilation, again by far outperforms the implementation of the new algorithm in Singular. This advancement in computing the Gromov-Witten invariants facilitates a study of number theoretic and geometric properties of the generating series, including quasi-modularity and homogeneity.


title: Systematic Analysis of Security and Vulnerabilities in Miniapps

abstract: The past few years have witnessed a boom of miniapps, as lightweight applications, miniapps are of great importance in the mobile internet sector. Consequently, the security of miniapps can directly impact compromising the integrity of sensitive data, posing a potential threat to user privacy. However, after a thorough review of the various research efforts in miniapp security, we found that their actions in researching the safety of miniapp web interfaces are limited. This paper proposes a triad threat model focusing on users, servers and attackers to mitigate the security risk of miniapps. By following the principle of least privilege and the direction of permission consistency, we design a novel analysis framework for the security risk assessment of miniapps by this model. Then, we analyzed the correlation between the security risk assessment and the threat model associated with the miniapp. This analysis led to identifying potential scopes and categorisations with security risks. In the case study, we identify nine major categories of vulnerability issues, such as SQL injection, logical vulnerabilities and cross-site scripting. We also assessed a total of 50,628 security risk hazards and provided specific examples.


title: A Survey of Emerging Applications of Diffusion Probabilistic Models in MRI

abstract: Diffusion probabilistic models (DPMs) which employ explicit likelihood characterization and a gradual sampling process to synthesize data, have gained increasing research interest. Despite their huge computational burdens due to the large number of steps involved during sampling, DPMs are widely appreciated in various medical imaging tasks for their high-quality and diversity of generation. Magnetic resonance imaging (MRI) is an important medical imaging modality with excellent soft tissue contrast and superb spatial resolution, which possesses unique opportunities for DPMs. Although there is a recent surge of studies exploring DPMs in MRI, a survey paper of DPMs specifically designed for MRI applications is still lacking. This review article aims to help researchers in the MRI community to grasp the advances of DPMs in different applications. We first introduce the theory of two dominant kinds of DPMs, categorized according to whether the diffusion time step is discrete or continuous, and then provide a comprehensive review of emerging DPMs in MRI, including reconstruction, image generation, image translation, segmentation, anomaly detection, and further research topics. Finally, we discuss the general limitations as well as limitations specific to the MRI tasks of DPMs and point out potential areas that are worth further exploration.


title: PIMSAB: A Processing-In-Memory System with Spatially-Aware Communication and Bit-Serial-Aware Computation

abstract: Bit-serial Processing-In-Memory (PIM) is an attractive paradigm for accelerator architectures, for parallel workloads such as Deep Learning (DL), because of its capability to achieve massive data parallelism at a low area overhead and provide orders-of-magnitude data movement savings by moving computational resources closer to the data. While many PIM architectures have been proposed, improvements are needed in communicating intermediate results to consumer kernels, for communication between tiles at scale, for reduction operations, and for efficiently performing bit-serial operations with constants. We present PIMSAB, a scalable architecture that provides spatially aware communication network for efficient intra-tile and inter-tile data movement and provides efficient computation support for generally inefficient bit-serial compute patterns. Our architecture consists of a massive hierarchical array of compute-enabled SRAMs (CRAMs) and is codesigned with a compiler to achieve high utilization. The key novelties of our architecture are: (1) providing efficient support for spatially-aware communication by providing local H-tree network for reductions, by adding explicit hardware for shuffling operands, and by deploying systolic broadcasting, and (2) taking advantage of the divisible nature of bit-serial computations through adaptive precision, bit-slicing and efficient handling of constant operations. When compared against a similarly provisioned modern Tensor Core GPU (NVIDIA A100), across common DL kernels and an end-to-end DL network (Resnet18), PIMSAB outperforms the GPU by 3x, and reduces energy by 4.2x. We compare PIMSAB with similarly provisioned state-of-the-art SRAM PIM (Duality Cache) and DRAM PIM (SIMDRAM) and observe a speedup of 3.7x and 3.88x respectively.


title: Multi-Task Reinforcement Learning with Mixture of Orthogonal Experts

abstract: Multi-Task Reinforcement Learning (MTRL) tackles the long-standing problem of endowing agents with skills that generalize across a variety of problems. To this end, sharing representations plays a fundamental role in capturing both unique and common characteristics of the tasks. Tasks may exhibit similarities in terms of skills, objects, or physical properties while leveraging their representations eases the achievement of a universal policy. Nevertheless, the pursuit of learning a shared set of diverse representations is still an open challenge. In this paper, we introduce a novel approach for representation learning in MTRL that encapsulates common structures among the tasks using orthogonal representations to promote diversity. Our method, named Mixture Of Orthogonal Experts (MOORE), leverages a Gram-Schmidt process to shape a shared subspace of representations generated by a mixture of experts. When task-specific information is provided, MOORE generates relevant representations from this shared subspace. We assess the effectiveness of our approach on two MTRL benchmarks, namely MiniGrid and MetaWorld, showing that MOORE surpasses related baselines and establishes a new state-of-the-art result on MetaWorld.


title: Blowups of smooth Fano hypersurfaces, their birational geometry and divisorial stability

abstract: Let $\Gamma$ be a smooth $k$-dimensional hypersurface in $\mathbb P^{k+1}$ and $X \supset \Gamma$ a smooth $n$-dimensional Fano hypersurface in $\mathbb P^{n+1}$ where $n\geq 3$ and $k\geq 1$. Let $Y \rightarrow X$ be the blowup of $X$ along $\Gamma$. We give a constructive proof that $Y$ is a Mori dream space. In particular, we describe its Mori chamber decomposition and the associated birational models of $Y$. We classify for which $X$ and $\Gamma$ the variety $Y$ is a Fano manifold and we initiate the study of K-stability of $Y$.


title: A Simple and Transparent Method for Improving the Energetics and Thermodynamics of Seawater Approximations: Static Energy Asymptotics (SEA)

abstract: The static energy encodes all possible information about the thermodynamics and potential energy (and all related forces) of stratified geophysical fluids. In this paper, we develop a systematic methodology, called static energy asymptotics, that exploits this property for constructing energetically and thermodynamically consistent sound-proof approximations of the equations of motion. By approximating the static energy to various orders of accuracy, two main families of approximations are (re-)derived and discussed: the pseudo-incompressible (PI) approximation and the anelastic (AN) approximation. For all approximations, the background and available potential energies (in Lorenz sense) can be constructed to match their exact counterparts as closely as feasible and to be expressible in terms of the exact (as opposed to ad-hoc) thermodynamic potentials. For hydrostatic motions, the AN approximation (of which the Boussinesq approximation is a special case) has the same structure as that of legacy Seawater Boussinesq primitive equations. The energetics of such models could therefore be made transparently traceable to that of the full Navier-Stokes equations at little to no additional cost, thus allowing them to take full advantage of the Gibbs Sea Water (GSW) library developed as part of the new thermodynamic standard for seawater TEOS-10


title: Machine Culture

abstract: The ability of humans to create and disseminate culture is often credited as the single most important factor of our success as a species. In this Perspective, we explore the notion of machine culture, culture mediated or generated by machines. We argue that intelligent machines simultaneously transform the cultural evolutionary processes of variation, transmission, and selection. Recommender algorithms are altering social learning dynamics. Chatbots are forming a new mode of cultural transmission, serving as cultural models. Furthermore, intelligent machines are evolving as contributors in generating cultural traits--from game strategies and visual art to scientific results. We provide a conceptual framework for studying the present and anticipated future impact of machines on cultural evolution, and present a research agenda for the study of machine culture.


title: Analytic continuation of multipoint correlation functions

abstract: Conceptually, the Matsubara formalism (MF), using imaginary frequencies, and the Keldysh formalism (KF), formulated in real frequencies, give equivalent results for systems in thermal equilibrium. The MF has less complexity and is thus more convenient than the KF. However, computing dynamical observables in the MF requires the analytic continuation from imaginary to real frequencies. The analytic continuation is well-known for two-point correlation functions (having one frequency argument), but, for multipoint correlators, a straightforward recipe for deducing all Keldysh components from the MF correlator had not been formulated yet. Recently, a representation of MF and KF correlators in terms of formalism-independent partial spectral functions and formalism-specific kernels was introduced by Kugler, Lee, and von Delft [Phys. Rev. X 11, 041006 (2021)]. We use this representation to formally elucidate the connection between both formalisms. We show how a multipoint MF correlator can be analytically continued to recover all partial spectral functions and yield all Keldysh components of its KF counterpart. The procedure is illustrated for various correlators of the Hubbard atom.


title: Addressing the speed-accuracy simulation trade-off for adaptive spiking neurons

abstract: The adaptive leaky integrate-and-fire (ALIF) model is fundamental within computational neuroscience and has been instrumental in studying our brains $\textit{in silico}$. Due to the sequential nature of simulating these neural models, a commonly faced issue is the speed-accuracy trade-off: either accurately simulate a neuron using a small discretisation time-step (DT), which is slow, or more quickly simulate a neuron using a larger DT and incur a loss in simulation accuracy. Here we provide a solution to this dilemma, by algorithmically reinterpreting the ALIF model, reducing the sequential simulation complexity and permitting a more efficient parallelisation on GPUs. We computationally validate our implementation to obtain over a $50\times$ training speedup using small DTs on synthetic benchmarks. We also obtained a comparable performance to the standard ALIF implementation on different supervised classification tasks - yet in a fraction of the training time. Lastly, we showcase how our model makes it possible to quickly and accurately fit real electrophysiological recordings of cortical neurons, where very fine sub-millisecond DTs are crucial for capturing exact spike timing.


title: A study guide for "Trilinear smoothing inequalities and a variant of the triangular Hilbert transform"

abstract: This article is a study guide for "Trilinear smoothing inequalities and a variant of the triangular Hilbert transform" by Christ, Durcik, and Roos. We first present the standard techniques in the study of oscillatory integrals with the simpler toy model of a Hilbert transform along a parabola. These standard techniques prove to be insufficient in the study of the triangular Hilbert transform with curvature. The central and novel idea in their proof of the $L^p$-boundedness of the triangular Hilbert transform with curvature is a trilinear smoothing inequality which we also examine in this article.