title: GRMHD beyond Kerr: An extension of the HARM code for thin disks to non-Kerr spacetimes

abstract: Black hole based tests of general relativity have proliferated in recent times with new and improved detectors and telescopes. Modelling of the black hole neighborhood, where most of the radiation carrying strong-field signature originates, is of utmost importance for robust and accurate constraints on possible violations of general relativity. As a first step, this paper presents the extension of general relativistic magnetohydrodynamic simulations of thin accretion disks to parametrically deformed black holes that generalize the Kerr solution. The extension is based on \textsc{harmpi}, a publicly available member of the \textsc{harm} family of codes, and uses a phenomenological metric to study parametric deviations away from Kerr. The extended model is used to study the disk structure, stability, and radiative efficiency. We also compute the Fe K$\alpha$ profiles in simplified scenarios and present an outlook for the future.


title: New possibilities in identification of binary choice models with fixed effects

abstract: We study the identification of binary choice models with fixed effects. We provide a condition called sign saturation and show that this condition is sufficient for the identification of the model. In particular, we can guarantee identification even with bounded regressors. We also show that without this condition, the model is not identified unless the error distribution belongs to a small class. The same sign saturation condition is also essential for identifying the sign of treatment effects. A test is provided to check the sign saturation condition and can be implemented using existing algorithms for the maximum score estimator.


title: Action of Hecke algebra on the double flag variety of type AIII

abstract: Consider a connected reductive algebraic group $ G $ and a symmetric subgroup $ K $. Let $ \mathfrak{X} = K/B_K \times G/P $ be a double flag variety of finite type, where $ B_K $ is a Borel subgroup of $ K $, and $ P $ a parabolic subgroup of $ G $. A general argument shows that the orbit space $ \mathbb{C}\,\mathfrak{X}/K $ inherits a natural action of the Hecke algebra $ \mathscr{H} = \mathscr{H}(K, B_K) $ of double cosets via convolutions. However, to find out the explicit structure of the Hecke module is a quite different problem. In this paper, we determine the explicit action of $ \mathscr{H} $ on $ \mathbb{C}\,\mathfrak{X}/K $ in a combinatorial way using graphs for the double flag variety of type AIII. As a by-product, we also get the description of the representation of the Weyl group on $ \mathbb{C}\,\mathfrak{X}/K $ as a direct sum of induced representations.


title: Survival Kernets: Scalable and Interpretable Deep Kernel Survival Analysis with an Accuracy Guarantee

abstract: Kernel survival analysis models estimate individual survival distributions with the help of a kernel function, which measures the similarity between any two data points. Such a kernel function can be learned using deep kernel survival models. In this paper, we present a new deep kernel survival model called a survival kernet, which scales to large datasets in a manner that is amenable to model interpretation and also theoretical analysis. Specifically, the training data are partitioned into clusters based on a recently developed training set compression scheme for classification and regression called kernel netting that we extend to the survival analysis setting. At test time, each data point is represented as a weighted combination of these clusters, and each such cluster can be visualized. For a special case of survival kernets, we establish a finite-sample error bound on predicted survival distributions that is, up to a log factor, optimal. Whereas scalability at test time is achieved using the aforementioned kernel netting compression strategy, scalability during training is achieved by a warm-start procedure based on tree ensembles such as XGBoost and a heuristic approach to accelerating neural architecture search. On four standard survival analysis datasets of varying sizes (up to roughly 3 million data points), we show that survival kernets are highly competitive compared to various baselines tested in terms of time-dependent concordance index. Our code is available at: https://github.com/georgehc/survival-kernets


title: De-biasing particle filtering for a continuous time hidden Markov model with a Cox process observation model

abstract: We develop a (nearly) unbiased particle filtering algorithm for a specific class of continuous-time state-space models, such that (a) the latent process $X_t$ is a linear Gaussian diffusion; and (b) the observations arise from a Poisson process with intensity $\lambda(X_t)$. The likelihood of the posterior probability density function of the latent process includes an intractable path integral. Our algorithm relies on Poisson estimates which approximate unbiasedly this integral. We show how we can tune these Poisson estimates to ensure that, with large probability, all but a few of the estimates generated by the algorithm are positive. Then replacing the negative estimates by zero leads to a much smaller bias than what would obtain through discretisation. We quantify the probability of negative estimates for certain special cases and show that our particle filter is effectively unbiased. We apply our method to a challenging 3D single molecule tracking example with a Born and Wolf observation model.


title: Policy Learning with Asymmetric Counterfactual Utilities

abstract: Data-driven decision making plays an important role even in high stakes settings like medicine and public policy. Learning optimal policies from observed data requires a careful formulation of the utility function whose expected value is maximized across a population. Although researchers typically use utilities that depend on observed outcomes alone, in many settings the decision maker's utility function is more properly characterized by the joint set of potential outcomes under all actions. For example, the Hippocratic principle to "do no harm" implies that the cost of causing death to a patient who would otherwise survive without treatment is greater than the cost of forgoing life-saving treatment. We consider optimal policy learning with asymmetric counterfactual utility functions of this form that consider the joint set of potential outcomes. We show that asymmetric counterfactual utilities lead to an unidentifiable expected utility function, and so we first partially identify it. Drawing on statistical decision theory, we then derive minimax decision rules by minimizing the maximum expected utility loss relative to different alternative policies. We show that one can learn minimax loss decision rules from observed data by solving intermediate classification problems, and establish that the finite sample excess expected utility loss of this procedure is bounded by the regret of these intermediate classifiers. We apply this conceptual framework and methodology to the decision about whether or not to use right heart catheterization for patients with possible pulmonary hypertension.


title: Learning to Estimate and Refine Fluid Motion with Physical Dynamics

abstract: Extracting information on fluid motion directly from images is challenging. Fluid flow represents a complex dynamic system governed by the Navier-Stokes equations. General optical flow methods are typically designed for rigid body motion, and thus struggle if applied to fluid motion estimation directly. Further, optical flow methods only focus on two consecutive frames without utilising historical temporal information, while the fluid motion (velocity field) can be considered a continuous trajectory constrained by time-dependent partial differential equations (PDEs). This discrepancy has the potential to induce physically inconsistent estimations. Here we propose an unsupervised learning based prediction-correction scheme for fluid flow estimation. An estimate is first given by a PDE-constrained optical flow predictor, which is then refined by a physical based corrector. The proposed approach outperforms optical flow methods and shows competitive results compared to existing supervised learning based methods on a benchmark dataset. Furthermore, the proposed approach can generalize to complex real-world fluid scenarios where ground truth information is effectively unknowable. Finally, experiments demonstrate that the physical corrector can refine flow estimates by mimicking the operator splitting method commonly utilised in fluid dynamical simulation.


title: Metallicity profiles of Ultra Diffuse Galaxies in NIHAO simulations

abstract: Supernovae feedback driven expansion has proven to be a viable mechanism to explain the average properties of Ultra Diffuse Galaxies (UDGs) such as the sizes, colors, mass and internal kinematics. Here, we explore the origin of stellar metallicity gradients in feedback driven simulated UDGs from the NIHAO project and compare them with the observed distribution of metallicity gradients of both Local Group dwarfs as well as of the recently observed UDG DF44. Simulated UDGs display a large variety of metallicity profiles, showing flat to negative gradients, similarly to what is observed in LG dwarfs, while DF44 data suggest a flat to positive gradient. The variety of metallicity gradients in simulations is set by the interplay between the radius at which star formation occurs and the subsequent supernovae feedback driven stellar redistribution: rotation supported systems tend to have flat metallicity profiles while dispersion supported galaxies show negative and steep profiles. Our results suggest that UDGs are not peculiar in what regards their metallicity gradients, when compared to regular dwarfs. Desirably, a larger observational sample of UDGs' gradients shall be available in the future, in order to test our predictions.


title: Random tensor networks with nontrivial links

abstract: Random tensor networks are a powerful toy model for understanding the entanglement structure of holographic quantum gravity. However, unlike holographic quantum gravity, their entanglement spectra are flat. It has therefore been argued that a better model consists of random tensor networks with link states that are not maximally entangled, i.e., have nontrivial spectra. In this work, we initiate a systematic study of the entanglement properties of these networks. We employ tools from free probability, random matrix theory, and one-shot quantum information theory to study random tensor networks with bounded and unbounded variation in link spectra, and in cases where a subsystem has one or multiple minimal cuts. If the link states have bounded spectral variation, the limiting entanglement spectrum of a subsystem with two minimal cuts can be expressed as a free product of the entanglement spectra of each cut, along with a Marchenko-Pastur distribution. For a class of states with unbounded spectral variation, analogous to semiclassical states in quantum gravity, we relate the limiting entanglement spectrum of a subsystem with two minimal cuts to the distribution of the minimal entanglement across the two cuts. In doing so, we draw connections to previous work on split transfer protocols, entanglement negativity in random tensor networks, and Euclidean path integrals in quantum gravity.


title: The common envelope origins of the fast jet in the planetary nebula M 3-38

abstract: We present the analysis of Multi-Espectr\'ografo en GTC de Alta Resoluci\'on para Astronom\'ia (MEGARA) high-dispersion integral field spectroscopic observations of the bipolar planetary nebula (PN) M 3-38. These observations unveil the presence of a fast outflow aligned with the symmetry axis of M 3-38 that expands with a velocity up to $\pm$225 km s$^{-1}$. The deprojected space velocity of this feature can be estimated to be $\approx$320$^{+130}_{-60}$ km s$^{-1}$, which together with its highly collimated morphology suggests that it is one of the fastest jet detected in a PN. We have also used Kepler observations of the central star of M 3-38 to unveil variability associated with a dominant period of 17.7 days. We attribute this to the presence of a low-mass star with an orbital separation of $\approx$0.12-0.16 au. The fast and collimated ejection and the close binary system point towards a common envelope formation scenario for M 3-38.


title: The weave product of two conics

abstract: Tensor diagrams are a handy way to depict complicated relationships between objects in projective geometry. One of the simpler ones takes two copies of a $3\times 3$ matrix and computes its adjugate. In this paper, we give a geometric interpretation of this construction when two different matrices are used. To do so, we interpret them as coefficient matrices of conic sections in the projective plane and relate the diagram construction to their common tangents.


title: Tight bounds for the learning of homotopy \`a la Niyogi, Smale, and Weinberger for subsets of Euclidean spaces and of Riemannian manifolds

abstract: In this article we extend and strengthen the seminal work by Niyogi, Smale, and Weinberger on the learning of the homotopy type from a sample of an underlying space. In their work, Niyogi, Smale, and Weinberger studied samples of $C^2$ manifolds with positive reach embedded in $\mathbb{R}^d$. We extend their results in the following ways: In the first part of our paper we consider both manifolds of positive reach -- a more general setting than $C^2$ manifolds -- and sets of positive reach embedded in $\mathbb{R}^d$. The sample $P$ of such a set $\mathcal{S}$ does not have to lie directly on it. Instead, we assume that the two one-sided Hausdorff distances -- $\varepsilon$ and $\delta$ -- between $P$ and $\mathcal{S}$ are bounded. We provide explicit bounds in terms of $\varepsilon$ and $ \delta$, that guarantee that there exists a parameter $r$ such that the union of balls of radius $r$ centred at the sample $P$ deformation-retracts to $\mathcal{S}$. In the second part of our paper we study homotopy learning in a significantly more general setting -- we investigate sets of positive reach and submanifolds of positive reach embedded in a \emph{Riemannian manifold with bounded sectional curvature}. To this end we introduce a new version of the reach in the Riemannian setting inspired by the cut locus. Yet again, we provide tight bounds on $\varepsilon$ and $\delta$ for both cases (submanifolds as well as sets of positive reach), exhibiting the tightness by an explicit construction.


title: Radio Nebul\ae\ from Hyper-Accreting X-ray Binaries as Common Envelope Precursors and Persistent Counterparts of Fast Radio Bursts

abstract: Roche lobe overflow from a donor star onto a black hole or neutron star binary companion can evolve to a phase of unstable runaway mass-transfer, lasting as short as hundreds of orbits ($\lesssim 10^{2}$ yr for a giant donor), and eventually culminating in a common envelope event. The highly super-Eddington accretion rates achieved during this brief phase ($\dot{M} \gtrsim 10^{5}\dot{M}_{\rm Edd})$ are accompanied by intense mass-loss in disk winds, analogous but even more extreme than ultra-luminous X-ray (ULX) sources in the nearby universe. Also in analogy with observed ULX, this expanding outflow will inflate an energetic `bubble' of plasma into the circumbinary medium. Embedded within this bubble is a nebula of relativistic electrons heated at the termination shock of the faster $v \gtrsim 0.1 c$ wind/jet from the inner accretion flow. We present a time-dependent, one-zone model for the synchrotron radio emission and other observable properties of such ULX `hyper-nebulae'. If ULX jets are sources of repeating fast radio bursts (FRB), as recently proposed, such hyper-nebulae could generate persistent radio emission and contribute large and time-variable rotation measure to the bursts, consistent with those seen from FRB 20121102 and FRB 190520B. ULX hyper-nebulae can be discovered independent of an FRB association in radio surveys such as VLASS, as off-nuclear point-sources whose fluxes can evolve significantly on timescales as short as years, possibly presaging energetic transients from common envelope mergers.


title: Reconfigurable Elastic Metamaterials: Engineering Dispersion with Meccano$^{\text{TM}}$

abstract: We design, simulate and experimentally characterise a reconfigurable elastic metamaterial with beyond-nearest-neighbour (BNN) coupling. The structure is composed from the popular British model construction system, Meccano$^{\text{TM}}$. The Meccano$^{\text{TM}}$ metamaterial supports backwards waves with opposite directions of phase and group velocities. We experimentally verify three distinct configurations and acoustically infer their spatial vibration spectra.


title: Theory of rotationally resolved two-dimensional infrared spectroscopy including polarization dependence and rotational coherence dynamics

abstract: Two-dimensional infrared (2DIR) spectroscopy is widely used to study molecular dynamics but it is typically restricted to solid and liquid phase samples and modest spectral resolution. Only recently, its potential to study gas-phase dynamics is beginning to be realized. Moreover, the recently proposed technique of cavity-enhanced 2D spectroscopy using frequency combs and developments in multi-comb spectroscopy are expected to dramatically advance capabilities for acquisition of rotationally-resolved 2DIR spectra. This demonstrates the need for rigorous and quantitative treatment of rotationally-resolved, polarization-dependent third-order response of gas-phase samples. In this article, we provide a rigorous and quantitative description of rotationally-resolved 2DIR spectroscopy using density-matrix, time-dependent perturbation theory and angular momentum algebra techniques. We describe the band and branch structure of 2D spectra, decompose the molecular response into polarization-dependence classes, use this decomposition to derive and explain special polarization conditions and relate the liquid-phase polarization conditions to gas-phase ones. Furthermore, we discuss the rotational coherence dynamics during waiting time.


title: Optimal Investment and Equilibrium Pricing under Ambiguity

abstract: We consider portfolio selection under nonparametric $\alpha$-maxmin ambiguity in the neighbourhood of a reference distribution. We show strict concavity of the portfolio problem under ambiguity aversion. Implied demand functions are nondifferentiable, resemble observed bid-ask spreads, and are consistent with existing parametric limiting participation results under ambiguity. Ambiguity seekers exhibit a discontinuous demand function, implying an empty set of reservation prices. If agents have identical, or sufficiently similar prior beliefs, the first-best equilibrium is no trade. Simple conditions yield the existence of a Pareto-efficient second-best equilibrium, implying that heterogeneity in ambiguity preferences is sufficient for mutually beneficial transactions among all else homogeneous traders. These equilibria reconcile many observed phenomena in liquid high-information financial markets, such as liquidity dry-ups, portfolio inertia, and negative risk premia.


title: Next-to-leading-order QCD Corrections to Higgs Production in association with a Jet

abstract: We compute the next-to-leading-order (NLO) QCD corrections to the Higgs pT distribution in Higgs production in association with a jet via gluon fusion at the LHC, with exact dependence on the mass of the quark circulating in the heavy-quark loops. The NLO corrections are presented including the top-quark mass, and for the first time, the bottom-quark mass as well. Further, besides the on-shell mass scheme, we consider for the first time a running mass renormalisation scheme. The computation is based on amplitudes which are valid for arbitrary heavy-quark masses.


title: Bi-Calibration Networks for Weakly-Supervised Video Representation Learning

abstract: The leverage of large volumes of web videos paired with the searched queries or surrounding texts (e.g., title) offers an economic and extensible alternative to supervised video representation learning. Nevertheless, modeling such weakly visual-textual connection is not trivial due to query polysemy (i.e., many possible meanings for a query) and text isomorphism (i.e., same syntactic structure of different text). In this paper, we introduce a new design of mutual calibration between query and text to boost weakly-supervised video representation learning. Specifically, we present Bi-Calibration Networks (BCN) that novelly couples two calibrations to learn the amendment from text to query and vice versa. Technically, BCN executes clustering on all the titles of the videos searched by an identical query and takes the centroid of each cluster as a text prototype. The query vocabulary is built directly on query words. The video-to-text/video-to-query projections over text prototypes/query vocabulary then start the text-to-query or query-to-text calibration to estimate the amendment to query or text. We also devise a selection scheme to balance the two corrections. Two large-scale web video datasets paired with query and title for each video are newly collected for weakly-supervised video representation learning, which are named as YOVO-3M and YOVO-10M, respectively. The video features of BCN learnt on 3M web videos obtain superior results under linear model protocol on downstream tasks. More remarkably, BCN trained on the larger set of 10M web videos with further fine-tuning leads to 1.6%, and 1.8% gains in top-1 accuracy on Kinetics-400, and Something-Something V2 datasets over the state-of-the-art TDN, and ACTION-Net methods with ImageNet pre-training. Source code and datasets are available at \url{https://github.com/FuchenUSTC/BCN}.


title: Controlling rotationally-resolved two-dimensional infrared spectra with polarization

abstract: Recent advancements in infrared frequency combs will enable recording of coherent two-dimensional infrared spectra of gas-phase molecules with rotational resolution (RR2DIR). We describe how RR2DIR spectra can be controlled with polarization to enhance or suppress certain features of the spectrum, with new polarization conditions unique to freely rotating molecules and absent in the condensed-phase. With the polarization control methods described here, RR2DIR spectroscopy can be a powerful tool for studying complex gas mixtures of polyatomic molecules.


title: Theoretical results for hadronic masses and their widths in the framework of the SO(4) model

abstract: The prediction of properties of the low energy portion of the hadronic spectrum is a challenging task which, up to day, is still tentatively given due to the non-perturbative nature of QCD at low energies. In this paper we are exploring the validity of the SO(4) scheme, as representative of the fundamental QCD structure of meson-like states in the region below 2.5 GeV. We have focussed the attention in the calculation of the energy and width of states of various spin, isospin and parities. The theoretical results are compared, sistematically, to the available experimental information.


title: Characterizing Distances Between Points in the Level Sets of a Class of Continuous Functions on a Closed Interval

abstract: Given a continuous function $f:[a,b]\to\mathbb{R}$ such that $f(a)=f(b)$, we investigate the set of distances $|x-y|$ where $f(x)=f(y)$. In particular, we show that the only distances this set must contain are ones which evenly divide $[a,b]$. Additionally, we show that it must contain at least one third of the interval $[0,b-a]$. Lastly, we explore some higher dimensional generalizations.


title: Online Coordination: Methods and Comparative Case Studies of Coordinated Groups across Four Events in the United States

abstract: Coordinated groups of user accounts working together in online social media can be used to manipulate the online discourse and thus is an important area of study. In this study, we work towards a general theory of coordination. There are many ways to coordinate groups online: semantic, social, referral and many more. Each represents a coordination dimension, where the more dimensions of coordination are present for one event, the stronger the coordination present. We build on existing approaches that detect coordinated groups by identifying high levels of synchronized actions within a specified time window. A key concern with this approach is the selection of the time window. We propose a method that selects the optimal window size to accurately capture local coordination while avoiding the capture of coincidental synchronicity. With this enhanced method of coordination detection, we perform a comparative study across four events: US Elections Primaries 2020, Reopen America 2020, Capitol Riots 2021 and COVID Vaccine Release 2021. Herein, we explore the following three dimensions of coordination for each event -- semantic, referral and social coordination -- and perform group and user analysis within and among the events. This allows us to expose different user coordination behavior patterns and identify narratives and user support themes, hence estimating the degree and theme of coordination.


title: Convex subsets of non-convex Lorentz balls

abstract: Many star bodies have convex subsets with approximately the same Gaussian measure (of the complement). Inspired by this phenomenon, and in connection with the randomized Dvoretzky theorem for Lorentz spaces, we derive bounds on the distribution of certain functions of a Gaussian random vector by approximating their sub-level sets by convex subsets.


title: A fixed point index approach to Krasnosel'skii-Precup fixed point theorem in cones and applications

abstract: We present an alternative approach to the vector version of Krasnosel'skii compression-expansion fixed point theorem due to Precup, which is based on the fixed point index. It allows us to obtain new general versions of this fixed point theorem and also multiplicity results. We emphasize that all of them are coexistence fixed point theorems for operator systems, that means that every component of the fixed points obtained is non-trivial. Finally, these coexistence fixed point theorems are applied to obtain results concerning the existence of positive solutions for systems of Hammerstein integral equations and radially symmetric solutions of $(p_1,p_2)$-Laplacian systems.


title: PlanBench: An Extensible Benchmark for Evaluating Large Language Models on Planning and Reasoning about Change

abstract: Generating plans of action, and reasoning about change have long been considered a core competence of intelligent agents. It is thus no surprise that evaluating the planning and reasoning capabilities of large language models (LLMs) has become a hot topic of research. Most claims about LLM planning capabilities are however based on common sense tasks-where it becomes hard to tell whether LLMs are planning or merely retrieving from their vast world knowledge. There is a strong need for systematic and extensible planning benchmarks with sufficient diversity to evaluate whether LLMs have innate planning capabilities. Motivated by this, we propose PlanBench, an extensible benchmark suite based on the kinds of domains used in the automated planning community, especially in the International Planning Competition, to test the capabilities of LLMs in planning or reasoning about actions and change. PlanBench provides sufficient diversity in both the task domains and the specific planning capabilities. Our studies also show that on many critical capabilities-including plan generation-LLM performance falls quite short, even with the SOTA models. PlanBench can thus function as a useful marker of progress of LLMs in planning and reasoning.


title: A Learning Aided Flexible Gradient Descent Approach to MISO Beamforming

abstract: This paper proposes a learning aided gradient descent (LAGD) algorithm to solve the weighted sum rate (WSR) maximization problem for multiple-input single-output (MISO) beamforming. The proposed LAGD algorithm directly optimizes the transmit precoder through implicit gradient descent based iterations, at each of which the optimization strategy is determined by a neural network, and thus, is dynamic and adaptive. At each instance of the problem, this network is initialized randomly, and updated throughout the iterative solution process. Therefore, the LAGD algorithm can be implemented at any signal-to-noise ratio (SNR) and for arbitrary antenna/user numbers, does not require labelled data or training prior to deployment. Numerical results show that the LAGD algorithm can outperform of the well-known WMMSE algorithm as well as other learning-based solutions with a modest computational complexity. Our code is available at https://github.com/XiaGroup/LAGD.


title: A survey on association schemes on triples

abstract: Association schemes on triples (ASTs) are ternary analogues of classical association schemes, whose relations and adjacency algebras are ternary instead of binary. We provide a survey of the current progress in the study of ASTs, highlighting open questions, suggesting research directions, and producing some related results. We review properties of the ternary adjacency algebras of ASTs, ASTs whose relations are invariant under some group action, and ASTs obtained from 2-designs and two-graphs. We also provide a notion of fusion and fission ASTs, using the AST obtained from the affine special linear group $ASL(2,q)$ as an example.


title: A Quantitative Analysis of Dynamic Mechanisms Regulating HIV Latency and Activation

abstract: Objective: The reservoir of human immunodeficiency virus (HIV) latently infected cells is the major obstacle for eradication of acquired immunodeficiency syndrome (AIDS). Due to the noisy environment and multiple influencing factors in the organism, current dynamical models cannot reach a common understanding of the molecular mechanism of HIV latency. In this work, through a new dynamical structure decomposition, the deterministic part of the equation can be separated from the stochastic noise. Thus, the fixed-point analysis of ordinary differential equation is enough to obtain the different steady states of the system. Methods: We established a dynamical model of HIV transcription process by using continuous stochastic differential equations, which simplifies the dimensions of equations needed to describe the system and increases the explorable space of the model. Different states between latency and activation of virus and their relations can be intuitively represented by potential functions and probability distribution functions. Results: Based on our model, the influence of different dynamical parameters on stability is quantitatively analyzed, the parameter ranges of the system in bistable and monostable states are obtained respectively. The theoretical basis of this work is verified by comparing the effects of different factors on dynamic bifurcation with the results of biological experiments. Conclusion: This paper goes beyond previous discrete stochastic methods, and can quantitatively analyze the dynamic mechanism of HIV transcriptional regulation through ordinary differential equations, which is beneficial to the promotion to deal with the high-dimensional situation, and further study the occurrence and development of AIDS in vivo, so as to guide the design of experiments and search for clinical treatment.


title: Quantum self-consistent equation-of-motion method for computing molecular excitation energies, ionization potentials, and electron affinities on a quantum computer

abstract: Near-term quantum computers are expected to facilitate material and chemical research through accurate molecular simulations. Several developments have already shown that accurate ground-state energies for small molecules can be evaluated on present-day quantum devices. Although electronically excited states play a vital role in chemical processes and applications, the search for a reliable and practical approach for routine excited-state calculations on near-term quantum devices is ongoing. Inspired by excited-state methods developed for the unitary coupled-cluster theory in quantum chemistry, we present an equation-of-motion-based method to compute excitation energies following the variational quantum eigensolver algorithm for ground-state calculations on a quantum computer. We perform numerical simulations on H$_2$, H$_4$, H$_2$O, and LiH molecules to test our quantum self-consistent equation-of-motion (q-sc-EOM) method and compare it to other current state-of-the-art methods. q-sc-EOM makes use of self-consistent operators to satisfy the vacuum annihilation condition, a critical property for accurate calculations. It provides real and size-intensive energy differences corresponding to vertical excitation energies, ionization potentials and electron affinities. We also find that q-sc-EOM is more suitable for implementation on NISQ devices as it is expected to be more resilient to noise compared with the currently available methods.


title: On a certified VMS-Smagorinsky Reduced Basis model with LPS pressure stabilisation

abstract: In this work we introduce a certified Reduced Basis VMS-Smagorinsky turbulence model with local projection stabilisation (LPS) on the pressure. We prove its stability for Taylor-Hood discretisations of velocity-pressure. We construct an \textit{a posteriori} error estimator for the snapshot selection through a Greedy algorithm, based on the Brezzi-Rappaz-Raviart theory of approximation of non-singular branches of non-linear PDEs. The Empirical Interpolation Method (EIM) is used for the approximation of the non-linear terms. We present some numerical tests in which we show an improved speedup on the computation of the reduced basis problem with the LPS pressure stabilisation, with respect to the method of using pressure supremizers.


title: A Theory of Sub-Barcodes

abstract: From the work of Bauer and Lesnick, it is known that there is no functor from the category of pointwise finite-dimensional persistence modules to the category of barcodes and overlap matchings. In this work, we introduce sub-barcodes and show that there is a functor from the category of factorizations of persistence module homomorphisms to a poset of barcodes ordered by the sub-barcode relation. Sub-barcodes and factorizations provide a looser alternative to bottleneck matchings and interleavings that can give strong guarantees in a number of settings that arise naturally in topological data analysis. The main use of sub-barcodes is to make strong claims about an unknown barcode in the absence of an interleaving. For example, given only upper and lower bounds $g\geq f\geq \ell$ of an unknown real-valued function $f$, a sub-barcode associated with $f$ can be constructed from $\ell$ and $g$ alone. We propose a theory of sub-barcodes and observe that the subobjects in the category of functors from intervals to matchings naturally correspond to sub-barcodes.


title: Reconstructing the Disrupted Dwarf Galaxy $Gaia$-Sausage/Enceladus Using its Stars and Globular Clusters

abstract: We combine spectroscopic, photometric, and astrometric information from APOGEE data release 17 and $Gaia$ early data release 3 to perform a self-consistent characterization of $Gaia$-Sausage/Enceladus (GSE), the remnant of the last major merger experienced by the Milky Way, considering stars and globular clusters (GCs) altogether. Our novel set of chemodynamical criteria to select genuine stars of GSE yields a metallicity distribution function with a median [Fe/H] of $-1.22$ dex and $0.23$ dex dispersion. Stars from GSE present an excess of [Al/Fe] and [Mg/Mn] (also [Mg/Fe]) in comparison to surviving Milky Way dwarf satellites, which can be explained by differences in star-formation efficiencies and timescales between these systems. However, stars from Sequoia, another proposed accreted halo substructure, essentially overlap the GSE footprint in all analyzed chemical-abundance spaces, but present lower metallicities. Among probable GCs of GSE with APOGEE observations available, we find no evidence for atypical [Fe/H] spreads with the exception of $\omega$ Centauri ($\omega$Cen). Under the assumption that $\omega$Cen is a stripped nuclear star cluster, we estimate the stellar mass of its progenitor to be $M_\star \approx 1.3 \times 10^9 M_\odot$, well-within literature expectations for GSE. This leads us to envision GSE as the best available candidate for the original host galaxy of $\omega$Cen. We also take advantage of $Gaia$'s photometry and APOGEE metallicities as priors to determine fundamental parameters for eight high-probability ($>$70%) GC members of GSE via statistical isochrone fitting. Finally, the newly determined ages and APOGEE [Fe/H] values are utilized to model the age-metallicity relation of GSE.


title: Finite Groups With Two Irredundant Covers

abstract: An irredundant cover of a finite group $G$ is a collection of proper subgroups whose union is $G$ and which contains no smaller subcover. We classify finite groups which possess exactly two irredundant covers, thereby initiating an answer to a question of Brodie, who classified finite groups with one irredundant cover.


title: SuperFGD prototype time resolution studies

abstract: The SuperFGD will be a part of the ND280 near detector of the T2K and Hyper Kamiokande projects, that will help to reduce systematic uncertainties related with neutrino flux and cross-section modeling. The upgraded ND280 will be able to perform a full exclusive reconstruction of the final state from neutrino-nucleus interactions, including measurements of low momentum protons, pions and, for the first time, event-by event measurements of neutron kinematics. The time resolution defines the neutron energy resolution. We present the results of time resolution measurements made with the SuperFGD prototype that consists of 9216 plastic scintillator cubes (cube size is 1 cm$^3$) readout with 1728 wavelength-shifting fibers going along three orthogonal directions. We use data from the muon beam exposure at CERN. The time resolution of 0.97 ns was obtained for one readout channel after implementing the time calibration with a correction for the time-walk effect. The time resolution improves with energy deposited in a scintillator cube. Averaging two readout channels for one scintillator cube improves the time resolution to 0.68 ns which means that signals in different channels are not synchronous. Therefore the contribution from the time recording step of 2.5 ns is averaged as well. Averaging time values from N channels improves the time resolution by $\sim 1/\sqrt{N}$. Therefore a very good time resolution should be achievable for neutrons since neutron recoils hit typically several scintillator cubes and in addition produce larger amplitudes than muons. Measurements performed with a laser and a wide-bandwidth oscilloscope demonstrated that the time resolution obtained with the muon beam is not far from its expected limit. The intrinsic time resolution of one channel is 0.67 ns for signals of 56 photo-electron typical for minimum ionizing particles.


title: Topological mean dimension of induced systems

abstract: For a topological system with positive topological entropy, we show that the induced transformation on the set of probability measures endowed with the weak-$*$ topology has infinite topological mean dimension. We also estimate the rate of divergence of the entropy with respect to the Wasserstein distance when the scale goes to zero.


title: Bayesian modeling and clustering for spatio-temporal areal data: An application to Italian unemployment

abstract: Spatio-temporal areal data can be seen as a collection of time series which are spatially correlated according to a specific neighboring structure. Incorporating the temporal and spatial dimension into a statistical model poses challenges regarding the underlying theoretical framework as well as the implementation of efficient computational methods. We propose to include spatio-temporal random effects using a conditional autoregressive prior, where the temporal correlation is modeled through an autoregressive mean decomposition and the spatial correlation by the precision matrix inheriting the neighboring structure. Their joint distribution constitutes a Gaussian Markov random field, whose sparse precision matrix enables the usage of efficient sampling algorithms. We cluster the areal units using a nonparametric prior, thereby learning latent partitions of the areal units. The performance of the model is assessed via an application to study regional unemployment patterns in Italy. When compared to other spatial and spatio-temporal competitors, the proposed model shows more precise estimates and the additional information obtained from the clustering allows for an extended economic interpretation of the unemployment rates of the Italian provinces.


title: A nation-wide experiment: fuel tax cuts and almost free public transport for three months in Germany -- Report 2 First wave results

abstract: In spring 2022, the German federal government agreed on a set of measures that aim at reducing households' financial burden resulting from a recent price increase, especially in energy and mobility. These measures include among others, a nation-wide public transport ticket for 9 EUR per month and a fuel tax cut that reduces fuel prices by more than 15%. In transportation research this is an almost unprecedented behavioral experiment. It allows to study not only behavioral responses in mode choice and induced demand but also to assess the effectiveness of transport policy instruments. We observe this natural experiment with a three-wave survey and an app-based travel diary on a sample of hundreds of participants as well as an analysis of traffic counts. In this second report, we update the information on study participation, provide first insights on the smartphone app usage as well as insights on the first wave results, particularly on the 9 EUR-ticket purchase intention.


title: Shift limits of a non-autonomous system

abstract: Let $t=t_1t_2\cdots$ be an element of the full shift with shift map $\tau$ on a finite set of characters $\mathcal{A}$ and let $ \Sigma=\text{ closure} \{\tau^i(t):\;i\in\N\cup\{0\}\}$. Let $f_t=f_{t_1,\,\infty}=\cdots\circ f_{t_2}\circ f_{t_1} $ be a non-autonomous system over a compact metric space $ X $ where $t_i\in \mathcal A $. The set $\F_t^+=\{f_{\tau^i(t)}:\; i\in\N\}$ is called the shifted family of $f_t$. If $t$ is a transitive point of the full shift on $\mathcal A$, then by introducing a natural topology, $\overline{\F_t^+}$ is a classical IFS; otherwise, $\overline{\F_t^+}=\{f_\sigma=f_{\sigma_1,\,\infty}:\; \sigma\in\Sigma\}$ is a generalized IFS. We will show that if $ f_t$ has some various shadowing and specification properties, then this is true for $f_{\sigma}\in\overline{\F^+_t}$; however, this claim is not true for other properties such as transitivity, mixing and exactness. Also, if $ \Sigma $ is sofic and $x\in X$ is periodic point for some $f_\sigma\in\overline{\F^+_t}$, then there is a periodic $\sigma'\in\Sigma$ such that $x$ is periodic for $f_{\sigma'}\in\overline{\F^+_t}$.


title: Effects of reactive, dissipative and rate-limited nonlinearity on the behaviour of superconducting resonator parametric amplifiers

abstract: We present a formalism for modelling parametric amplification by resonators subject to rate-limited nonlinearity of mixed reactive/dissipative character, with particular relevance to superconducting devices. The non-linearity is assumed to be characterised by a single state parameter, which responds to changes in the energy stored in the resonator with finite response time. We show how the operating point and small signal amplification behaviour of the pumped resonator can be calculated, characterised and optimised in terms of a set of three dimensionless parameters. The formalism is then illustrated with a simple, first-order, model nonlinearity and the implications for amplification via quasiparticle generation in a superconductor discussed. Throughout we describe how the parameters needed to characterise the device can be determined experimentally from steady-state measurements. A key result of this paper is that rate-limiting of a nonlinear mechanism does not preclude amplification, although it does limit the bandwidth over which it may be achieved.


title: Fundamental-mode eigenfrequencies of neutral and charged twin neutron stars

abstract: We investigate the effects of rapid and slow conversions on the fundamental-mode eigenfrequencies of hybrid neutron stars having highly discontinuous transitions between hadronic and quark matter, the so-called twin stars. We analyze some characteristic cases in the available parameter space of the equations of state for the hadronic and quark phases. Furthermore, we also consider the possibility that these stellar configurations are electrically charged. Our results indicate that for neutral configurations under rapid conversions the stability window coincides with the usual stability criterion, i.e. $\partial M/\partial p_c > 0$ and that the two branches are disconnected. This discontinuity remains when electric charge is considered, but the usual criterion is not sufficient to determine the star's stability. On the other hand, slow conversions connect initially disconnected branches and increase the stability window of the hybrid configurations. For both conversion speeds, the presence of electric charge diminishes the magnitude of the eigenfrequencies and its stability window.


title: How to quantise probabilities while preserving their convex order

abstract: We introduce an algorithm which, given probabilities $\mu \leq_{\text{cx}} \nu$ in convex order and defined on a separable Banach space $B$, constructs finitely-supported approximations $\mu_n \to \mu, \nu_n\to \nu$ which are in convex order $\mu_n \leq_{\text{cx}} \nu_n$. We provide upper-bounds for the speed of convergence, in terms of the Wasserstein distance. We discuss the (dis)advantages of our algorithm and its link with the discretisation of the Martingale Optimal Transport problem, and we illustrate its implementation with numerical examples. We study the operation which, given $\mu$/$\nu$ and some (finite) partition of $B$, outputs $\mu_n$/$\nu_n$, showing that applied to a probability $\gamma$ and to all partitions it outputs the set of all probabilities $\zeta \leq_{\text{cx}} \gamma$.


title: Root-$T \overline{T}$ Deformations in Two-Dimensional Quantum Field Theories

abstract: In this letter we introduce a one-parameter deformation of two-dimensional quantum field theories generated by a non-analytic operator which we call Root-$T \overline{T}$. For a conformal field theory, the operator coincides with the square-root of the $T \overline{T}$ operator. More generally, the operator is defined so that classically it is marginal and generates a flow which commutes with the $T \overline{T}$-flow. Intriguingly, the Root-$T \overline{T}$ flow is closely related to the ModMax theory recently constructed by Bandos, Lechner, Sorokin and Townsend.

categories: hep-th cond-mat.stat-mech math-ph math.MP nlin.SI


title: Roadmap to Autonomous Surgery -- A Framework to Surgical Autonomy

abstract: Robotic surgery has increased the domain of surgeries possible. Several examples of partial surgical automation have been seen in the past decade. We break down the path of automation tasks into features required and provide a checklist that can help reach higher levels of surgical automation. Finally, we discuss the current challenges and advances required to make this happen.


title: What Goes Around: Leveraging a Constant-curvature Motion Constraint in Radar Odometry

abstract: This paper presents a method that leverages vehicle motion constraints to refine data associations in a point-based radar odometry system. By using the strong prior on how a non-holonomic robot is constrained to move smoothly through its environment, we develop the necessary framework to estimate ego-motion from a single landmark association rather than considering all of these correspondences at once. This allows for informed outlier detection of poor matches that are a dominant source of pose estimate error. By refining the subset of matched landmarks, we see an absolute decrease of 2.15% (from 4.68% to 2.53%) in translational error, approximately halving the error in odometry (reducing by 45.94%) than when using the full set of correspondences. This contribution is relevant to other point-based odometry implementations that rely on a range sensor and provides a lightweight and interpretable means of incorporating vehicle dynamics for ego-motion estimation.


title: Overcoming High Frequency Limitations of Current-Mode Control Using a Control Conditioning Approach -- Part I: Modeling and Analysis

abstract: Current-mode control is one of the most popular controller strategies for power converters. With the advent of wide bandgap devices including GaN and SiC, higher switching frequencies have become more viable at higher power because of lower switching losses. However, the advantage of higher switching frequency for faster, higher bandwidth control is squandered because of current sensor interference. We present a framework for characterizing and analyzing this interference as uncertainties to the controller model. These uncertainties introduce additional dynamics and nonlinearity that can result in instability and poor transient performance of the current control loop. In this paper, we provide a model framework based on a new control conditioning approach that guarantees global stability and a strategy for optimizing transient performance. In Part II of this paper series, we present the analysis, design, and hardware validation of three effective solutions.


title: Long-range Order and Directional Defect Propagation in the Nonreciprocal XY Model with Vision Cone Interactions

abstract: We study a two-dimensional, nonreciprocal XY model, where each spin interacts only with its nearest neighbours in a certain angle around its current orientation, i.e., its `vision cone'. Using energetic arguments and Monte-Carlo simulations we show that a true long-range ordered phase emerges. A necessary ingredient is a configuration-dependent bond dilution entailed by the vision cones. Strikingly, defects propagate in a directional manner, thereby breaking the parity and time-reversal symmetry of the spin dynamics. This is detectable by a non-zero entropy production rate.


title: SFace: Privacy-friendly and Accurate Face Recognition using Synthetic Data

abstract: Recent deep face recognition models proposed in the literature utilized large-scale public datasets such as MS-Celeb-1M and VGGFace2 for training very deep neural networks, achieving state-of-the-art performance on mainstream benchmarks. Recently, many of these datasets, e.g., MS-Celeb-1M and VGGFace2, are retracted due to credible privacy and ethical concerns. This motivates this work to propose and investigate the feasibility of using a privacy-friendly synthetically generated face dataset to train face recognition models. Towards this end, we utilize a class-conditional generative adversarial network to generate class-labeled synthetic face images, namely SFace. To address the privacy aspect of using such data to train a face recognition model, we provide extensive evaluation experiments on the identity relation between the synthetic dataset and the original authentic dataset used to train the generative model. Our reported evaluation proved that associating an identity of the authentic dataset to one with the same class label in the synthetic dataset is hardly possible. We also propose to train face recognition on our privacy-friendly dataset, SFace, using three different learning strategies, multi-class classification, label-free knowledge transfer, and combined learning of multi-class classification and knowledge transfer. The reported evaluation results on five authentic face benchmarks demonstrated that the privacy-friendly synthetic dataset has high potential to be used for training face recognition models, achieving, for example, a verification accuracy of 91.87\% on LFW using multi-class classification and 99.13\% using the combined learning strategy.


title: Robustness against data loss with Algebraic Statistics

abstract: The paper describes an algorithm that, given an initial design $\mathcal{F}_n$ of size $n$ and a linear model with $p$ parameters, provides a sequence $\mathcal{F}_n \supset \ldots \supset \mathcal{F}_{n-k} \supset \ldots \supset \mathcal{F}_p$ of nested \emph{robust} designs. The sequence is obtained by the removal, one by one, of the runs of $\mathcal{F}_n$ till a $p$-run \emph{saturated} design $\mathcal{F}_p$ is obtained. The potential impact of the algorithm on real applications is high. The initial fraction $\mathcal{F}_n$ can be of any type and the output sequence can be used to organize the experimental activity. The experiments can start with the runs corresponding to $\mathcal{F}_p$ and continue adding one run after the other (from $\mathcal{F}_{n-k}$ to $\mathcal{F}_{n-k+1}$) till the initial design $\mathcal{F}_n$ is obtained. In this way, if for some unexpected reasons the experimental activity must be stopped before the end when only $n-k$ runs are completed, the corresponding $\mathcal{F}_{n-k}$ has a high value of robustness for $k \in \{1, \ldots, n-p\}$. The algorithm uses the circuit basis, a special representation of the kernel of a matrix with integer entries. The effectiveness of the algorithm is demonstrated through the use of simulations.


title: Ideal Polytopes for Representations of $GL_n(\mathbb{C})$

abstract: In this paper we use the superpotential for the flag variety $GL_n/B$ and particular coordinate systems that we call ideal coordinates for $\mathbf{i}$, to construct polytopes $\mathcal{P}^{\mathbf{i}}_\lambda$ inside $\mathbb{R}^{R_+}$, associated to highest weight representations $V_\lambda$ of $GL_n(\mathbb{C})$. Here $\mathbf{i}$ is a reduced expression of the longest element of the Weyl group and $R_+$ is the set of positive roots of $GL_n$. The lattice points of $\mathcal{P}^{\mathbf{i}}_\lambda$ can be used to encode a basis of the representation $V_\lambda$. In particular, for a specific choice of $\mathbf{i}$, the polytope $\mathcal{P}^{\mathbf{i}}_\lambda$ is unimodularly equivalent to a Gelfand-Tsetlin polytope. The construction of the polytopes involves tropicalisation of the superpotential. Using work of Judd (arXiv:1606.06883) we have that there is a unique positive critical point of the superpotential $\mathcal{W}_{t^\lambda}$ over the field of Puiseux series. Its coordinates, in terms of the ideal coordinates for $\mathbf{i}$, are positive in the sense of having positive leading term. The remarkable property of our new polytopes relates to the tropical version of this critical point, which, for every choice of $\mathbf{i}$, gives a point in $\mathbb{R}^{R_+}$ that lies in the interior of the polytope $\mathcal{P}^{\mathbf{i}}_\lambda$. We prove that this tropical critical point is independent of the reduced expression $\mathbf{i}$, and that it is given by a pattern called the ideal filling for $\lambda$ that was introduced by Judd. Finally, combining these results with work of Rietsch (arXiv:math/0511124) relating critical points of the superpotential with Toeplitz matrices, we show that for a totally positive lower-triangular Toeplitz matrix over the field of Puiseux series factorized into simple root subgroups, the valuations of the factors give an ideal filling.


title: Overcoming High Frequency Limitations of Current-Mode Control Using a Control Conditioning Approach -- Part II: Implementation and Hardware

abstract: This article is the second part of a paper series about interference in extremum (i.e., peak or valley) current-mode control, which applies to both fixed and variable switching frequency power converters. Specifically, this part presents three control conditioning methods that mitigate the adverse effect of interference. These methods are new ways to use: (i) slope compensation; (ii) low-pass filtering; and (iii) the phenomenon of comparator-overdrive-delay, for control conditioning. The stability criterion, closed-loop dynamics, and transient performance are derived with mathematical rigor for each method. The design tradeoffs are illustrated, discussed, and compared. The effectiveness of all three methods are demonstrated and validated in hardware using a power converter operating at multi-MHz switching frequencies.


title: Lyapunov Density Models: Constraining Distribution Shift in Learning-Based Control

abstract: Learned models and policies can generalize effectively when evaluated within the distribution of the training data, but can produce unpredictable and erroneous outputs on out-of-distribution inputs. In order to avoid distribution shift when deploying learning-based control algorithms, we seek a mechanism to constrain the agent to states and actions that resemble those that it was trained on. In control theory, Lyapunov stability and control-invariant sets allow us to make guarantees about controllers that stabilize the system around specific states, while in machine learning, density models allow us to estimate the training data distribution. Can we combine these two concepts, producing learning-based control algorithms that constrain the system to in-distribution states using only in-distribution actions? In this work, we propose to do this by combining concepts from Lyapunov stability and density estimation, introducing Lyapunov density models: a generalization of control Lyapunov functions and density models that provides guarantees on an agent's ability to stay in-distribution over its entire trajectory.


title: PRIVIC: A privacy-preserving method for incremental collection of location data

abstract: With recent advancements in technology, the threats of privacy violations of individuals' sensitive data are surging. Location data, in particular, have been shown to carry a substantial amount of sensitive information. A standard method to mitigate the privacy risks for location data consists in adding noise to the true values to achieve geo-indistinguishability (geo-ind). However, geo-ind alone is not sufficient to cover all privacy concerns. In particular, isolated locations are not sufficiently protected by the state-of-the-art Laplace mechanism (LAP) for geo-ind. In this paper, we focus on a mechanism based on the Blahut-Arimoto algorithm (BA) from the rate-distortion theory. We show that BA, in addition to providing geo-ind, enforces an elastic metric that mitigates the problem of isolation. Furthermore, BA provides an optimal trade-off between information leakage and quality of service. We then proceed to study the utility of BA in terms of the statistics that can be derived from the reported data, focusing on the inference of the original distribution. To this purpose, we de-noise the reported data by applying the iterative Bayesian update (IBU), an instance of the expectation-maximization method. It turns out that BA and IBU are dual to each other, and as a result, they work well together, in the sense that the statistical utility of BA is quite good and better than LAP for high privacy levels. Exploiting these properties of BA and IBU, we propose an iterative method, PRIVIC, for a privacy-friendly incremental collection of location data from users by service providers. We illustrate the soundness and functionality of our method both analytically and with experiments.


title: QuantFace: Towards Lightweight Face Recognition by Synthetic Data Low-bit Quantization

abstract: Deep learning-based face recognition models follow the common trend in deep neural networks by utilizing full-precision floating-point networks with high computational costs. Deploying such networks in use-cases constrained by computational requirements is often infeasible due to the large memory required by the full-precision model. Previous compact face recognition approaches proposed to design special compact architectures and train them from scratch using real training data, which may not be available in a real-world scenario due to privacy concerns. We present in this work the QuantFace solution based on low-bit precision format model quantization. QuantFace reduces the required computational cost of the existing face recognition models without the need for designing a particular architecture or accessing real training data. QuantFace introduces privacy-friendly synthetic face data to the quantization process to mitigate potential privacy concerns and issues related to the accessibility to real training data. Through extensive evaluation experiments on seven benchmarks and four network architectures, we demonstrate that QuantFace can successfully reduce the model size up to 5x while maintaining, to a large degree, the verification performance of the full-precision model without accessing real training datasets.


title: Note on Connes spectral distances of qubits

abstract: By virtue of the mathematical tools in noncommutative geometry, we study the Connes spectral distances between one- and two-qubit states. We construct a spectral triple corresponding to the $2D$ fermionic phase spaces, and calculate the Connes spectral distance between one qubits. Based on the Connes spectral distance, we define a coherence measure of quantum states, and calculate the coherence of one-qubit states. We also study some simple cases about two-qubit states, and the corresponding spectral distances satisfy the Pythagoras theorem. We find that the Connes spectral distances are different from quantum trace distances. The Connes spectral distance can be considered as a significant supplement to the trace distances in quantum information sciences. These results are significant for the study of physical relations and geometric structures of qubits and other quantum states.


title: A quantum information perspective on meson melting

abstract: We propose to use quantum information notions to characterize thermally induced melting of nonperturbative bound states at high temperatures. We apply tensor networks to investigate this idea in static and dynamical settings within the Ising quantum field theory, where bound states are confined fermion pairs - mesons. An equilibrium signature of meson melting is identified in the temperature dependence of the thermal-state second R\'enyi entropy, which varies from exponential to power-law scaling. Out of equilibrium, we identify as the relevant signature the transition from an oscillatory to a linear growing behavior of reflected entropy after a thermal quench. These analyses apply more broadly, which brings new ways of describing in-medium meson phenomena in quantum many-body and high-energy physics.

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


title: Magnetic Forces in Paramagnetic Fluids

abstract: An overview of the effect of a magnetic field gradient on fluids with linear magnetic susceptibilities is given. It is shown that two commonly encountered expressions, the magnetic field gradient force and the concentration gradient force for paramagnetic species in solution are equivalent for incompressible fluids. The magnetic field gradient and concentration gradient forces are approximations of the Kelvin force and Korteweg-Helmholtz force densities, respectively. The criterion for the appearance of magnetically induced convection is derived. Experimental work in which magnetically induced convection plays a role is reviewed.


title: Non-Abelian covering and new recursion operators for the 4D Mart\'inez Alonso-Shabat equation

abstract: We present new recursion operators for (shadows of nonlocal) symmetries of the 4D Mart\'inez Alonso-Shabat equation $u_{ty} = u_z u_{xy} - u_y u_{xz}$, and we show that their actions can produce new symmetries which are not contained in the Lie algebra of nonlocal symmetries presented in [I.S.Krasil'shchik, P.Voj\v{c}\'{a}k, On the algebra of nonlocal symmetries for the 4D Mart\'inez Alonso-Shabat equation. J. of Geom. and Phys. 163, (2021), 104122, (arXiv:2008.10281v1)]. To this end, we construct a non-Abelian covering of the equation in question using the Lax pair with two non-removable parameters.


title: Neural Transformers for Intraductal Papillary Mucosal Neoplasms (IPMN) Classification in MRI images

abstract: Early detection of precancerous cysts or neoplasms, i.e., Intraductal Papillary Mucosal Neoplasms (IPMN), in pancreas is a challenging and complex task, and it may lead to a more favourable outcome. Once detected, grading IPMNs accurately is also necessary, since low-risk IPMNs can be under surveillance program, while high-risk IPMNs have to be surgically resected before they turn into cancer. Current standards (Fukuoka and others) for IPMN classification show significant intra- and inter-operator variability, beside being error-prone, making a proper diagnosis unreliable. The established progress in artificial intelligence, through the deep learning paradigm, may provide a key tool for an effective support to medical decision for pancreatic cancer. In this work, we follow this trend, by proposing a novel AI-based IPMN classifier that leverages the recent success of transformer networks in generalizing across a wide variety of tasks, including vision ones. We specifically show that our transformer-based model exploits pre-training better than standard convolutional neural networks, thus supporting the sought architectural universalism of transformers in vision, including the medical image domain and it allows for a better interpretation of the obtained results.


title: Terabit Indoor Laser-Based Wireless Communications: LiFi 2.0 for 6G

abstract: This paper provides a summary of available technologies required for implementing indoor laser-based wireless networks capable of achieving aggregate data-rates of terabits per second as widely accepted as a sixth generation (6G) key performance indicator. The main focus of this paper is on the technologies supporting the near infrared region of the optical spectrum. The main challenges in the design of the transmitter and receiver systems and communication/networking schemes are identified and new insights are provided. This paper also covers the previous and recent standards as well as industrial applications for optical wireless communications (OWC) and LiFi.


title: RRT and RRT* Using Vehicle Dynamics

abstract: The trajectory derived from RRT and RRT* is zagged. A holonomic drive is able to follow this trajectory. But real-life vehicle which has dynamical constraints cannot follow this trajectory. In this work, we are going to modify the RRT and RRT* algorithm to generate a trajectory that a vehicle with dynamical constraint can follow. The continuous nature of steering control and acceleration control in a real-world vehicle introduces the complexity in its model. To introduce constraint in the vehicle's motion, while reducing the number of control and hence complexity, we are modeling our vehicle as a Dubins car. A Dubins car has only three controls (turning left, turning right, and moving forward) with a fixed velocity which makes our model simple. We use dubins curve (path that dubins car can follow) to trace the trajectory in RRT and RRT* algorithm.


title: Chiral flow in a binary mixture of two-dimensional active disks

abstract: We study, experimentally, the dynamics of a binary mixture of air-fluidized disks. The disks are chiral since they incorporate a set of blades with constant tilt. Both species are identical except for their blades tilt angle, which is rotated by 180o in the second species. We analyze the phase behavior of the system. Our analysis reveals a wide range of different fluid dynamics, including chiral flow. This chiral flow features in its base state a large vortex. We report, for certain ranges of relative particle density of each species, inversion of the vorticity of this vortex. We discuss on the possible mechanisms behind these chiral flow transitions.


title: EpiGRAF: Rethinking training of 3D GANs

abstract: A very recent trend in generative modeling is building 3D-aware generators from 2D image collections. To induce the 3D bias, such models typically rely on volumetric rendering, which is expensive to employ at high resolutions. During the past months, there appeared more than 10 works that address this scaling issue by training a separate 2D decoder to upsample a low-resolution image (or a feature tensor) produced from a pure 3D generator. But this solution comes at a cost: not only does it break multi-view consistency (i.e. shape and texture change when the camera moves), but it also learns the geometry in a low fidelity. In this work, we show that it is possible to obtain a high-resolution 3D generator with SotA image quality by following a completely different route of simply training the model patch-wise. We revisit and improve this optimization scheme in two ways. First, we design a location- and scale-aware discriminator to work on patches of different proportions and spatial positions. Second, we modify the patch sampling strategy based on an annealed beta distribution to stabilize training and accelerate the convergence. The resulted model, named EpiGRAF, is an efficient, high-resolution, pure 3D generator, and we test it on four datasets (two introduced in this work) at $256^2$ and $512^2$ resolutions. It obtains state-of-the-art image quality, high-fidelity geometry and trains ${\approx} 2.5 \times$ faster than the upsampler-based counterparts. Project website: https://universome.github.io/epigraf.


title: HealNet -- Self-Supervised Acute Wound Heal-Stage Classification

abstract: Identifying, tracking, and predicting wound heal-stage progression is a fundamental task towards proper diagnosis, effective treatment, facilitating healing, and reducing pain. Traditionally, a medical expert might observe a wound to determine the current healing state and recommend treatment. However, sourcing experts who can produce such a diagnosis solely from visual indicators can be difficult, time-consuming and expensive. In addition, lesions may take several weeks to undergo the healing process, demanding resources to monitor and diagnose continually. Automating this task can be challenging; datasets that follow wound progression from onset to maturation are small, rare, and often collected without computer vision in mind. To tackle these challenges, we introduce a self-supervised learning scheme composed of (a) learning embeddings of wound's temporal dynamics, (b) clustering for automatic stage discovery, and (c) fine-tuned classification. The proposed self-supervised and flexible learning framework is biologically inspired and trained on a small dataset with zero human labeling. The HealNet framework achieved high pre-text and downstream classification accuracy; when evaluated on held-out test data, HealNet achieved 97.7% pre-text accuracy and 90.62% heal-stage classification accuracy.


title: Neural network analysis of X-ray polarimeter data

abstract: This chapter presents deep neural network based methods for enhancing the sensitivity of X-ray telescopic observations with imaging polarimeters. Deep neural networks can be used to determine photoelectron emission directions, photon absorptions points, and photon energies from 2D photoelectron track images, with estimates for both the statistical and model uncertainties. Deep neural network predictive uncertainties can be incorporated into a weighted maximum likelihood to estimate source polarization parameters. Events converting outside of the fiducial gas volume, whose tracks have little polarization sensitivity, complicate polarization estimation. Deep neural network based classifiers can be used to select against these events to improve energy resolution and polarization sensitivity. The performance of deep neural network methods is compared against standard data analysis methods, revealing a < 0.75x improvement in minimum detectable polarization for IXPE-specific simulations. Potential future developments and improvements to these methods are discussed.


title: Coherent noise enables probabilistic sequence replay in spiking neuronal networks

abstract: Animals rely on different decision strategies when faced with ambiguous or uncertain cues. Depending on the context, decisions may be biased towards events that were most frequently experienced in the past, or be more explorative. A particular type of decision making central to cognition is sequential memory recall in response to ambiguous cues. A previously developed spiking neuronal network implementation of sequence prediction and recall learns complex, high-order sequences in an unsupervised manner by local, biologically inspired plasticity rules. In response to an ambiguous cue, the model deterministically recalls the sequence shown most frequently during training. Here, we present an extension of the model enabling a range of different decision strategies. In this model, explorative behavior is generated by supplying neurons with noise. As the model relies on population encoding, uncorrelated noise averages out, and the recall dynamics remain effectively deterministic. In the presence of locally correlated noise, the averaging effect is avoided without impairing the model performance, and without the need for large noise amplitudes. We investigate two forms of correlated noise occurring in nature: shared synaptic background inputs, and random locking of the stimulus to spatiotemporal oscillations in the network activity. Depending on the noise characteristics, the network adopts various replay strategies. This study thereby provides potential mechanisms explaining how the statistics of learned sequences affect decision making, and how decision strategies can be adjusted after learning.


title: Cascade of transitions in twisted and non-twisted graphene layers within the van Hove scenario

abstract: Motivated by measurements of compressibility and STM spectra in twisted bilayer graphene, we analyze the pattern of symmetry breaking for itinerant fermions near a van Hove singularity. Making use of an approximate SU(4) symmetry of the Landau functional, we show that the structure of the spin/isospin order parameter changes with increasing filling via a cascade of transitions. We compute the feedback from different spin/isospin orders on fermions and argue that each order splits the initially 4-fold degenerate van Hove peak in a particular fashion, consistent with the STM data and compressibility measurements, providing a unified interpretation of the cascade of transitions in twisted bilayer graphene. Our results follow from a generic analysis of an SU(4)-symmetric Landau functional and are valid beyond a specific underlying fermionic model. We argue that an analogous van Hove scenario explains the cascade of phase transitions in non-twisted Bernal bilayer and rhombohedral trilayer graphene.


title: Rethinking Symbolic Regression Datasets and Benchmarks for Scientific Discovery

abstract: This paper revisits datasets and evaluation criteria for Symbolic Regression (SR), specifically focused on its potential for scientific discovery. Focused on a set of formulas used in the existing datasets based on Feynman Lectures on Physics, we recreate 120 datasets to discuss the performance of symbolic regression for scientific discovery (SRSD). For each of the 120 SRSD datasets, we carefully review the properties of the formula and its variables to design reasonably realistic sampling ranges of values so that our new SRSD datasets can be used for evaluating the potential of SRSD such as whether or not an SR method can (re)discover physical laws from such datasets. We also create another 120 datasets that contain dummy variables to examine whether SR methods can choose necessary variables only. Besides, we propose to use normalized edit distances (NED) between a predicted equation and the true equation trees for addressing a critical issue that existing SR metrics are either binary or errors between the target values and an SR model's predicted values for a given input. We conduct benchmark experiments on our new SRSD datasets using various representative SR methods. The experimental results show that we provide a more realistic performance evaluation, and our user study shows that the NED correlates with human judges significantly more than an existing SR metric. We publish repositories of our code and 240 SRSD datasets.


title: Two-component parton fractional quantum Hall state in graphene

abstract: We study the $\nu={\pm}1/2$ fractional quantum Hall states in graphene observed by Zibrov et al. [Nat. Phys. 14, 930 (2018)]. The parton construction is employed to provide a valley unpolarized trial wave function for these states. The lattice scale valley anisotropic terms in the Hamiltonian soften the repulsion between electrons in different valleys to favor the valley unpolarized state. The validity of our proposal is corroborated by numerical calculations.


title: One-axis twisting as a method of generating many-body Bell correlations

abstract: We demonstrate that the one-axis twisting (OAT), a versatile method of creating non-classical states of bosonic qubits, is a powerful source of many-body Bell correlations. We develop a fully analytical and universal treatment of the process, which allows us to identify the critical time at which the Bell correlations emerge, and predict the depth of Bell correlations at all subsequent times. Our findings are illustrated with a highly non-trivial example of the OAT dynamics generated using the Bose-Hubbard model.


title: Faster Diffusion Cardiac MRI with Deep Learning-based breath hold reduction

abstract: Diffusion Tensor Cardiac Magnetic Resonance (DT-CMR) enables us to probe the microstructural arrangement of cardiomyocytes within the myocardium in vivo and non-invasively, which no other imaging modality allows. This innovative technology could revolutionise the ability to perform cardiac clinical diagnosis, risk stratification, prognosis and therapy follow-up. However, DT-CMR is currently inefficient with over six minutes needed to acquire a single 2D static image. Therefore, DT-CMR is currently confined to research but not used clinically. We propose to reduce the number of repetitions needed to produce DT-CMR datasets and subsequently de-noise them, decreasing the acquisition time by a linear factor while maintaining acceptable image quality. Our proposed approach, based on Generative Adversarial Networks, Vision Transformers, and Ensemble Learning, performs significantly and considerably better than previous proposed approaches, bringing single breath-hold DT-CMR closer to reality.


title: Multi-UAV Planning for Cooperative Wildfire Coverage and Tracking with Quality-of-Service Guarantees

abstract: In recent years, teams of robot and Unmanned Aerial Vehicles (UAVs) have been commissioned by researchers to enable accurate, online wildfire coverage and tracking. While the majority of prior work focuses on the coordination and control of such multi-robot systems, to date, these UAV teams have not been given the ability to reason about a fire's track (i.e., location and propagation dynamics) to provide performance guarantee over a time horizon. Motivated by the problem of aerial wildfire monitoring, we propose a predictive framework which enables cooperation in multi-UAV teams towards collaborative field coverage and fire tracking with probabilistic performance guarantee. Our approach enables UAVs to infer the latent fire propagation dynamics for time-extended coordination in safety-critical conditions. We derive a set of novel, analytical temporal, and tracking-error bounds to enable the UAV-team to distribute their limited resources and cover the entire fire area according to the case-specific estimated states and provide a probabilistic performance guarantee. Our results are not limited to the aerial wildfire monitoring case-study and are generally applicable to problems, such as search-and-rescue, target tracking and border patrol. We evaluate our approach in simulation and provide demonstrations of the proposed framework on a physical multi-robot testbed to account for real robot dynamics and restrictions. Our quantitative evaluations validate the performance of our method accumulating 7.5x and 9.0x smaller tracking-error than state-of-the-art model-based and reinforcement learning benchmarks, respectively.


title: The Impact of Visibility on the Right to Opt-out of Sale under CCPA

abstract: The California Consumer Protection Act (CCPA) gives users the right to opt-out of sale of their personal information, but prior work has found that opt-out mechanisms provided under this law result in very low opt-out rates. Privacy signals offer a solution for users who are willing to proactively take steps to enable privacy-enhancing tools, but many users are not aware of their rights under CCPA. We therefore explore an alternative approach to enhancing privacy under CCPA: increasing the visibility of opt-out of sale mechanisms. We conduct an user study with 54 participants and find that visible, standardized banners significantly increase opt-out of sale rates in the wild. Participants also report less difficulty opting out and more satisfaction with opt-out mechanisms compared to the native mechanisms currently provided by websites. Our results suggest that effective privacy regulation depends on imposing clear, enforceable visibility standards, and that CCPA's requirements for opt-out of sale mechanisms fall short.


title: FedHiSyn: A Hierarchical Synchronous Federated Learning Framework for Resource and Data Heterogeneity

abstract: Federated Learning (FL) enables training a global model without sharing the decentralized raw data stored on multiple devices to protect data privacy. Due to the diverse capacity of the devices, FL frameworks struggle to tackle the problems of straggler effects and outdated models. In addition, the data heterogeneity incurs severe accuracy degradation of the global model in the FL training process. To address aforementioned issues, we propose a hierarchical synchronous FL framework, i.e., FedHiSyn. FedHiSyn first clusters all available devices into a small number of categories based on their computing capacity. After a certain interval of local training, the models trained in different categories are simultaneously uploaded to a central server. Within a single category, the devices communicate the local updated model weights to each other based on a ring topology. As the efficiency of training in the ring topology prefers devices with homogeneous resources, the classification based on the computing capacity mitigates the impact of straggler effects. Besides, the combination of the synchronous update of multiple categories and the device communication within a single category help address the data heterogeneity issue while achieving high accuracy. We evaluate the proposed framework based on MNIST, EMNIST, CIFAR10 and CIFAR100 datasets and diverse heterogeneous settings of devices. Experimental results show that FedHiSyn outperforms six baseline methods, e.g., FedAvg, SCAFFOLD, and FedAT, in terms of training accuracy and efficiency.


title: Eikonal Approximation in Celestial CFT

abstract: We identify an eikonal regime in celestial CFT$_2$ in which massless 2-2 scattering is dominated by t-channel exchange. We derive a formula for the celestial amplitude that resums exchanges of arbitrary integer spin to all orders in the coupling. The resulting eikonal phase takes the same form as in flat space with the powers of center-of-mass energy replaced by weight-shifting operators on the celestial sphere. We independently compute the celestial two-point function for a scalar propagating in a shockwave background and show that to leading order in the gravitational coupling and for a suitable choice of the source, the result agrees with the prediction from the celestial eikonal formula for graviton exchange. We demonstrate that this two-point function can be directly obtained from the corresponding formula in AdS$_4$ in a flat space limit. We finally establish a general relation between scalar celestial amplitudes in celestial CFT$_{d-1}$ and the flat space limit of scalar AdS$_{d+1}$ Witten diagrams.


title: Nonlinear Dynamic analysis of vector-host model for Zika infection with predatory fish Gambusia Affinis

abstract: In the present paper, we study the dynamics of a nine compartmental vector-host model for Zika virus infection where the predatory fish Gambusia Affinis is introduced into the system to control the zika infection by preying on the vector. The system has six practically feasible equilibrium points where four of them are disease-free, and the rest are endemic. We discuss the existence and stability conditions for the equilibria. We find that when sexual transmission of zika comes to a halt then in absence of mosquitoes infection cannot persist. Hence, one needs to eradicate mosquitoes to eradicate infection. Moreover, we deduce that in the case of zika infection pushing the basic reproduction number below unity is next to impossible. Therefore, O_0, the mosquito survival threshold parameter, and O, the mosquito survival threshold parameter with predation play a crucial role in getting rid of the infection in respective cases since mosquitoes cannot survive when these are less than unity. Sensitivity analysis shows the importance of reducing mosquito biting rate and mutual contact rates between vector and host. It exhibits the importance of increasing the natural mortality rate of vectors to reduce the basic reproduction number. Numerical simulation shows that when the basic reproduction number is close but greater than unity, the introduction of a small amount of predatory fish Gambusia Affinis can completely swipe off the infection. In case of high transmission or high basic reproduction number, this fish increases the susceptible human population and keeps the infection under control, hence, prohibiting it from becoming an epidemic.


title: Strong Simulation of Linear Optical Processes

abstract: In this paper, we provide an algorithm and general framework for the simulation of photons passing through linear optical interferometers. Given $n$ photons at the input of an $m$-mode interferometer, our algorithm computes the probabilities of all possible output states with time complexity $O\left({n\binom{n+m-1}{m-1}}\right)$, linear in the number of output states $\binom{n+m-1}{m-1}$. It outperforms the permanent-based method by an exponential factor, and for the restricted problem of computing the probability for one given output it improves the time complexity over the state-of-the-art for the permanent of matrices with multiple rows or columns, with a tradeoff in the memory usage. Our algorithm also has additional versatility by virtue of its use of memorisation -- the storing of intermediate results -- which is advantageous in situations where several input states may be of interest. Additionally it allows for hybrid simulations, in which outputs are sampled from output states whose probability exceeds a given threshold, or from a restricted set of states. We consider a concrete, optimised implementation, and we benchmark the efficiency of our approach compared to existing tools.


title: (Certified!!) Adversarial Robustness for Free!

abstract: In this paper we show how to achieve state-of-the-art certified adversarial robustness to 2-norm bounded perturbations by relying exclusively on off-the-shelf pretrained models. To do so, we instantiate the denoised smoothing approach of Salman et al. 2020 by combining a pretrained denoising diffusion probabilistic model and a standard high-accuracy classifier. This allows us to certify 71% accuracy on ImageNet under adversarial perturbations constrained to be within an 2-norm of 0.5, an improvement of 14 percentage points over the prior certified SoTA using any approach, or an improvement of 30 percentage points over denoised smoothing. We obtain these results using only pretrained diffusion models and image classifiers, without requiring any fine tuning or retraining of model parameters.


title: On the effectiveness of persistent homology

abstract: Persistent homology (PH) is one of the most popular methods in Topological Data Analysis. Even though PH has been used in many different types of applications, the reasons behind its success remain elusive; in particular, it is not known for which classes of problems it is most effective, or to what extent it can detect geometric or topological features. The goal of this work is to identify some types of problems where PH performs well or even better than other methods in data analysis. We consider three fundamental shape analysis tasks: the detection of the number of holes, curvature and convexity from 2D and 3D point clouds sampled from shapes. Experiments demonstrate that PH is successful in these tasks, outperforming several baselines, including PointNet, an architecture inspired precisely by the properties of point clouds. In addition, we observe that PH remains effective for limited computational resources and limited training data, as well as out-of-distribution test data, including various data transformations and noise. For convexity detection, we provide a theoretical guarantee that PH is effective for this task in $\mathbb{R}^d$, and demonstrate the detection of a convexity measure on the FLAVIA data set of plant leaf images. Due to the crucial role of shape classification in understanding mathematical and physical structures and objects, and in many applications, the findings of this work will provide some knowledge about the types of problems that are appropriate for PH, so that it can - to borrow the words from Wigner 1960 - ``remain valid in future research, and extend, to our pleasure", but to our lesser bafflement, to a variety of applications.


title: Vicinity Vision Transformer

abstract: Vision transformers have shown great success on numerous computer vision tasks. However, its central component, softmax attention, prohibits vision transformers from scaling up to high-resolution images, due to both the computational complexity and memory footprint being quadratic. Although linear attention was introduced in natural language processing (NLP) tasks to mitigate a similar issue, directly applying existing linear attention to vision transformers may not lead to satisfactory results. We investigate this problem and find that computer vision tasks focus more on local information compared with NLP tasks. Based on this observation, we present a Vicinity Attention that introduces a locality bias to vision transformers with linear complexity. Specifically, for each image patch, we adjust its attention weight based on its 2D Manhattan distance measured by its neighbouring patches. In this case, the neighbouring patches will receive stronger attention than far-away patches. Moreover, since our Vicinity Attention requires the token length to be much larger than the feature dimension to show its efficiency advantages, we further propose a new Vicinity Vision Transformer (VVT) structure to reduce the feature dimension without degenerating the accuracy. We perform extensive experiments on the CIFAR100, ImageNet1K, and ADE20K datasets to validate the effectiveness of our method. Our method has a slower growth rate of GFlops than previous transformer-based and convolution-based networks when the input resolution increases. In particular, our approach achieves state-of-the-art image classification accuracy with 50% fewer parameters than previous methods.


title: Uncertainty Quantification for Competency Assessment of Autonomous Agents

abstract: For safe and reliable deployment in the real world, autonomous agents must elicit appropriate levels of trust from human users. One method to build trust is to have agents assess and communicate their own competencies for performing given tasks. Competency depends on the uncertainties affecting the agent, making accurate uncertainty quantification vital for competency assessment. In this work, we show how ensembles of deep generative models can be used to quantify the agent's aleatoric and epistemic uncertainties when forecasting task outcomes as part of competency assessment.


title: On-Site Potential Creates Complexity in Systems with Disordered Coupling

abstract: We calculate the average number of critical points $\overline{\mathcal{N}}$ of the energy landscape of a many-body system with disordered two-body interactions and a weak on-site potential. We find that introducing a weak nonlinear on-site potential dramatically increases $\overline{\mathcal{N}}$ to exponential in system size and give a complete picture of the organization of critical points. Our results extend solvable spin-glass models to physically more realistic models and are of relevance to glassy systems, nonlinear oscillator networks and many-body interacting systems.

categories: cond-mat.dis-nn cond-mat.stat-mech math-ph math.MP


title: LargeKernel3D: Scaling up Kernels in 3D Sparse CNNs

abstract: Recent advance in 2D CNNs has revealed that large kernels are important. However, when directly applying large convolutional kernels in 3D CNNs, severe difficulties are met, where those successful module designs in 2D become surprisingly ineffective on 3D networks, including the popular depth-wise convolution. To address this vital challenge, we instead propose the spatial-wise partition convolution and its large-kernel module. As a result, it avoids the optimization and efficiency issues of naive 3D large kernels. Our large-kernel 3D CNN network, LargeKernel3D, yields notable improvement in 3D tasks of semantic segmentation and object detection. It achieves 73.9% mIoU on the ScanNetv2 semantic segmentation and 72.8% NDS nuScenes object detection benchmarks, ranking 1st on the nuScenes LIDAR leaderboard. The performance further boosts to 74.2% NDS with a simple multi-modal fusion. In addition, LargeKernel3D can be scaled to 17x17x17 kernel size on Waymo 3D object detection. For the first time, we show that large kernels are feasible and essential for 3D visual tasks.


title: Weak approximation on Ch\^{a}telet surfaces

abstract: We study weak approximation for Ch\^{a}telet surfaces over number fields when all singular fibers are defined over rational points. We consider Ch\^{a}telet surfaces which satisfy weak approximation over every finite extension of the ground field. We prove many of these results by showing that the Brauer-Manin obstruction vanishes, then apply results of Colliot-Th\'el\`ene, Sansuc, and Swinnerton-Dyer.


title: A Study of Weisfeiler-Leman Colorings on Planar Graphs

abstract: The Weisfeiler-Leman (WL) algorithm is a combinatorial procedure that computes colorings on graphs, which can often be used to detect their (non-)isomorphism. Particularly the 1- and 2-dimensional versions 1-WL and 2-WL have received much attention, due to their numerous links to other areas of computer science. Knowing the expressive power of a certain dimension of the algorithm usually amounts to understanding the computed colorings. An increase in the dimension leads to finer computed colorings and, thus, more graphs can be distinguished. For example, on the class of planar graphs, 3-WL solves the isomorphism problem. However, the expressive power of 2-WL on the class is poorly understood (and, in particular, it may even well be that it decides isomorphism). In this paper, we investigate the colorings computed by 2-WL on planar graphs. Towards this end, we analyze the graphs induced by edge color classes in the graph. Based on the obtained classification, we show that for every 3-connected planar graph, it holds that: a) after coloring all pairs with their 2-WL color, the graph has fixing number 1 with respect to 1-WL, or b) there is a 2-WL-definable matching that can be used to transform the graph into a smaller one, or c) 2-WL detects a connected subgraph that is essentially the graph of a Platonic or Archimedean solid, a prism, a cycle, or a bipartite graph K_{2,\ell}. In particular, the graphs from case (a) are identified by 2-WL.


title: EnvPool: A Highly Parallel Reinforcement Learning Environment Execution Engine

abstract: There has been significant progress in developing reinforcement learning (RL) training systems. Past works such as IMPALA, Apex, Seed RL, Sample Factory, and others, aim to improve the system's overall throughput. In this paper, we aim to address a common bottleneck in the RL training system, i.e., parallel environment execution, which is often the slowest part of the whole system but receives little attention. With a curated design for paralleling RL environments, we have improved the RL environment simulation speed across different hardware setups, ranging from a laptop and a modest workstation, to a high-end machine such as NVIDIA DGX-A100. On a high-end machine, EnvPool achieves one million frames per second for the environment execution on Atari environments and three million frames per second on MuJoCo environments. When running EnvPool on a laptop, the speed is 2.8x that of the Python subprocess. Moreover, great compatibility with existing RL training libraries has been demonstrated in the open-sourced community, including CleanRL, rl_games, DeepMind Acme, etc. Finally, EnvPool allows researchers to iterate their ideas at a much faster pace and has great potential to become the de facto RL environment execution engine. Example runs show that it only takes five minutes to train agents to play Atari Pong and MuJoCo Ant on a laptop. EnvPool is open-sourced at https://github.com/sail-sg/envpool.