title: Nonthermal radiation of evaporating black holes

abstract: Black hole (BH) evaporation is caused by creation of entangled particle-antiparticle pairs near the event horizon, with one carrying positive energy to infinity and the other carrying negative energy into the BH. Since under the event horizon, particles always move toward the BH center, they can only be absorbed but not emitted at the center. This breaks absorption-emission symmetry and, as a result, annihilation of the particle at the BH center is described by a non-Hermitian Hamiltonian. We show that due to entanglement between photons moving inside and outside the event horizon, nonunitary absorption of the negative energy photons near the BH center, alters the outgoing radiation. As a result, radiation of the evaporating BH is not thermal; it carries information about BH interior, and entropy is preserved during evaporation.


title: Synthetic-reflection self-injection-locked microcombs

abstract: Laser-driven microresonators have enabled chip-integrated light sources with unique properties, including the self-organized formation of ultrashort soliton pulses and frequency combs (microcombs). While poised to impact major photonic applications, such as spectroscopy, sensing and optical data processing, microcombs still necessitate complex scientific equipment to achieve and maintain suitable single-pulse operation. Here, to address this challenge, we demonstrate microresonators with programmable synthetic reflection providing an injection-feedback to the driving laser. When designed appropriately, synthetic reflection enables deterministic access to self-injection-locked microcombs operating exclusively in the single-soliton regime. These results provide a route to easily-operable microcombs for portable sensors, autonomous navigation, or extreme-bandwidth data processing. The novel concept of synthetic reflection may also be generalized to other integrated photonic systems.


title: Falsification of Internal and External Validity in Observational Studies via Conditional Moment Restrictions

abstract: Randomized Controlled Trials (RCT)s are relied upon to assess new treatments, but suffer from limited power to guide personalized treatment decisions. On the other hand, observational (i.e., non-experimental) studies have large and diverse populations, but are prone to various biases (e.g. residual confounding). To safely leverage the strengths of observational studies, we focus on the problem of falsification, whereby RCTs are used to validate causal effect estimates learned from observational data. In particular, we show that, given data from both an RCT and an observational study, assumptions on internal and external validity have an observable, testable implication in the form of a set of Conditional Moment Restrictions (CMRs). Further, we show that expressing these CMRs with respect to the causal effect, or "causal contrast", as opposed to individual counterfactual means, provides a more reliable falsification test. In addition to giving guarantees on the asymptotic properties of our test, we demonstrate superior power and type I error of our approach on semi-synthetic and real world datasets. Our approach is interpretable, allowing a practitioner to visualize which subgroups in the population lead to falsification of an observational study.


title: The fundamental theorem of calculus in differential rings

abstract: In this paper, we study the consequences of the fundamental theorem of calculus from an algebraic point of view. For functions with singularities, this leads to a generalized notion of evaluation. We investigate properties of such integro-differential rings and discuss many examples. We also construct corresponding integro-differential operators and provide normal forms via rewrite rules. They are then used to derive several identities and properties in a purely algebraic way, generalizing well-known results from analysis. In identities like shuffle relations for nested integrals and the Taylor formula, additional terms are obtained that take singularities into account. Another focus lies on treating basics of linear ODEs in this framework of integro-differential operators. These operators can have matrix coefficients, which allow to treat systems of arbitrary size in a unified way. In the appendix, using tensor reduction systems, we give the technical details of normal forms and prove them for operators including other functionals besides evaluation.


title: A Herschel study of the high-mass protostar IRAS20126+4104

abstract: We performed Herschel observations of the continuum and line emission from the high-mass star-forming region IRAS20126+4104, which hosts a well-studied B-type (proto)star powering a bipolar outflow and is associated with a Keplerian circumstellar disk. The continuum images at six wavelengths allowed us to derive an accurate estimate of the bolometric luminosity and mass of the molecular clump enshrouding the disk. The same region has been mapped in 12 rotational transitions of carbon monoxide, which were used in synergy with the continuum data to determine the temperature and density distribution inside the clump and improve upon the mass estimate. The maps of two fine structure oxygen far-IR lines were used to estimate the volume density of the shocked region at the surface of the southern lobe of the outflow and the mass-loss rate. Our findings lend further support to the scenario previously proposed by various authors, confirming that at the origin of the bolometric luminosity and bipolar outflow from IRAS20126+4104 is a B-type star located at the centre of the Keplerian disk.


title: Contrastive Meta-Learning for Partially Observable Few-Shot Learning

abstract: Many contrastive and meta-learning approaches learn representations by identifying common features in multiple views. However, the formalism for these approaches generally assumes features to be shared across views to be captured coherently. We consider the problem of learning a unified representation from partial observations, where useful features may be present in only some of the views. We approach this through a probabilistic formalism enabling views to map to representations with different levels of uncertainty in different components; these views can then be integrated with one another through marginalisation over that uncertainty. Our approach, Partial Observation Experts Modelling (POEM), then enables us to meta-learn consistent representations from partial observations. We evaluate our approach on an adaptation of a comprehensive few-shot learning benchmark, Meta-Dataset, and demonstrate the benefits of POEM over other meta-learning methods at representation learning from partial observations. We further demonstrate the utility of POEM by meta-learning to represent an environment from partial views observed by an agent exploring the environment.


title: Multiple pion pair production in a Regge-based model

abstract: Central diffractive event topologies at the LHC energies can be identified by two different approaches. First, the forward scattered protons can be measured in Roman pots. Second, a veto on hadronic activity away from midrapidity can be imposed to define a double-gap topology. Such a double-gap topology trigger has been implemented by the ALICE collaboration in Run 1 and Run 2 of the LHC. The analysis of these events allows to determine the charged-particle multiplicity within the acceptance. The excellent particle identification capabilities of ALICE allows to study two-track events both in the pion and kaon sector. Events with measured charged particle multiplicity larger than two can arise from multiple pair production. A Regge-based approach for modeling such multiple pair production is presented.


title: Towards the Linear Algebra Based Taxonomy of XAI Explanations

abstract: This paper proposes an alternative approach to the basic taxonomy of explanations produced by explainable artificial intelligence techniques. Methods of Explainable Artificial Intelligence (XAI) were developed to answer the question why a certain prediction or estimation was made, preferably in terms easy to understand by the human agent. XAI taxonomies proposed in the literature mainly concentrate their attention on distinguishing explanations with respect to involving the human agent, which makes it complicated to provide a more mathematical approach to distinguish and compare different explanations. This paper narrows its attention to the cases where the data set of interest belongs to $\mathbb{R} ^n$ and proposes a simple linear algebra-based taxonomy for local explanations.


title: A Novel Framework for Policy Mirror Descent with General Parameterization and Linear Convergence

abstract: Modern policy optimization methods in reinforcement learning, such as TRPO and PPO, owe their success to the use of parameterized policies. However, while theoretical guarantees have been established for this class of algorithms, especially in the tabular setting, the use of general parameterization schemes remains mostly unjustified. In this work, we introduce a novel framework for policy optimization based on mirror descent that naturally accommodates general parameterizations. The policy class induced by our scheme recovers known classes, e.g., softmax, and generates new ones depending on the choice of mirror map. Using our framework, we obtain the first result that guarantees linear convergence for a policy-gradient-based method involving general parameterization. To demonstrate the ability of our framework to accommodate general parameterization schemes, we provide its sample complexity when using shallow neural networks, show that it represents an improvement upon the previous best results, and empirically validate the effectiveness of our theoretical claims on classic control tasks.


title: Phase conjugation with spatially incoherent light in complex media

abstract: Shaping light deep inside complex media, such as biological tissue, is critical to many research fields. Although the coherent control of scattered light via wavefront shaping has made significant advances in addressing this challenge, controlling light over extended or multiple targets without physical access inside a medium remains elusive. Here we present a phase conjugation method for spatially incoherent light, which enables the non-invasive light control based on incoherent emission from multiple target positions. Our method characterizes the scattering responses of hidden sources by retrieving mutually incoherent scattered fields from speckle patterns. By time-reversing scattered fluorescence with digital phase conjugation, we experimentally demonstrate focusing of light on individual and multiple targets. We also demonstrate maximum energy delivery to an extended target through a scattering medium by exploiting transmission eigenchannels. This paves the way to control light propagation in complex media using incoherent contrasts mechanisms.


title: Consistency Regularisation in Varying Contexts and Feature Perturbations for Semi-Supervised Semantic Segmentation of Histology Images

abstract: Semantic segmentation of various tissue and nuclei types in histology images is fundamental to many downstream tasks in the area of computational pathology (CPath). In recent years, Deep Learning (DL) methods have been shown to perform well on segmentation tasks but DL methods generally require a large amount of pixel-wise annotated data. Pixel-wise annotation sometimes requires expert's knowledge and time which is laborious and costly to obtain. In this paper, we present a consistency based semi-supervised learning (SSL) approach that can help mitigate this challenge by exploiting a large amount of unlabelled data for model training thus alleviating the need for a large annotated dataset. However, SSL models might also be susceptible to changing context and features perturbations exhibiting poor generalisation due to the limited training data. We propose an SSL method that learns robust features from both labelled and unlabelled images by enforcing consistency against varying contexts and feature perturbations. The proposed method incorporates context-aware consistency by contrasting pairs of overlapping images in a pixel-wise manner from changing contexts resulting in robust and context invariant features. We show that cross-consistency training makes the encoder features invariant to different perturbations and improves the prediction confidence. Finally, entropy minimisation is employed to further boost the confidence of the final prediction maps from unlabelled data. We conduct an extensive set of experiments on two publicly available large datasets (BCSS and MoNuSeg) and show superior performance compared to the state-of-the-art methods.


title: Self-Compressing Neural Networks

abstract: This work focuses on reducing neural network size, which is a major driver of neural network execution time, power consumption, bandwidth, and memory footprint. A key challenge is to reduce size in a manner that can be exploited readily for efficient training and inference without the need for specialized hardware. We propose Self-Compression: a simple, general method that simultaneously achieves two goals: (1) removing redundant weights, and (2) reducing the number of bits required to represent the remaining weights. This is achieved using a generalized loss function to minimize overall network size. In our experiments we demonstrate floating point accuracy with as few as 3% of the bits and 18% of the weights remaining in the network.


title: RRT Guided Model Predictive Path Integral Method

abstract: This work presents an optimal sampling-based method to solve the real-time motion planning problem in static and dynamic environments, exploiting the Rapid-exploring Random Trees (RRT) algorithm and the Model Predictive Path Integral (MPPI) algorithm. The RRT algorithm provides a nominal mean value of the random control distribution in the MPPI algorithm, resulting in satisfactory control performance in static and dynamic environments without a need for fine parameter tuning. We also discuss the importance of choosing the right mean of the MPPI algorithm, which balances exploration and optimality gap, given a fixed sample size. In particular, a sufficiently large mean is required to explore the state space enough, and a sufficiently small mean is required to guarantee that the samples reconstruct the optimal controls. The proposed methodology automates the procedure of choosing the right mean by incorporating the RRT algorithm. The simulations demonstrate that the proposed algorithm can solve the motion planning problem in real-time for static or dynamic environments.


title: Missing Data in Discrete Time State-Space Modeling of Ecological Momentary Assessment Data: A Monte-Carlo Study of Imputation Methods

abstract: When using ecological momentary assessment data (EMA), missing data is pervasive as participant attrition is a common issue. Thus, any EMA study must have a missing data plan. In this paper, we discuss missingness in time series analysis and the appropriate way to handle missing data when the data is modeled as a discrete time continuous measure state-space model. We found that Missing Completely At Random, Missing At Random, and Time-dependent Missing At Random data have less bias and variability than Autoregressive Time-dependent Missing At Random and Missing Not At Random. The Kalman filter excelled at handling missing data. Contrary to the literature, we found that, with either default package settings or a lag-1 imputation model, multiple imputation struggled to recover the parameters.


title: Accurate and efficient multiscale simulation of a heterogeneous elastic beam via computation on small sparse patches

abstract: Modern `smart' materials have complex microscale structure, often with unknown macroscale closure. The Equation-Free Patch Scheme empowers us to non-intrusively, efficiently, and accurately simulate over large scales through computations on only small well-separated patches of the microscale system. Here the microscale system is a solid beam of random heterogeneous elasticity. The continuing challenge is to compute the given physics on just the microscale patches, and couple the patches across un-simulated macroscale space, in order to establish efficiency, accuracy, consistency, and stability on the macroscale. Dynamical systems theory supports the scheme. This research program is to develop a systematic non-intrusive approach, both computationally and analytically proven, to model and compute accurately macroscale system levels of general complex physical and engineering systems.


title: Enhancing Neural Network Differential Equation Solvers

abstract: We motivate the use of neural networks for the construction of numerical solutions to differential equations. We prove that there exists a feed-forward neural network that can arbitrarily minimise an objective function that is zero at the solution of Poisson's equation, allowing us to guarantee that neural network solution estimates can get arbitrarily close to the exact solutions. We also show how these estimates can be appreciably enhanced through various strategies, in particular through the construction of error correction networks, for which we propose a general method. We conclude by providing numerical experiments that attest to the validity of all such strategies for variants of Poisson's equation.


title: Neargeodesics in Gromov hyperbolic John domains in Banach spaces

abstract: In this paper, we prove that neargeodesics in Gromov hyperbolic John domains in Banach space are cone arcs. This result gives an improvement of a result of Li [Theorem 1, Int. J. Math. 25 (2014)].


title: Realistic pattern formations on surfaces by adding arbitrary roughness

abstract: We are interested in generating surfaces with arbitrary roughness and forming patterns on the surfaces. Two methods are applied to construct rough surfaces. In the first method, some superposition of wave functions with random frequencies and angles of propagation are used to get periodic rough surfaces with analytic parametric equations. The amplitude of such surfaces is also an important variable in the provided eigenvalue analysis for the Laplace-Beltrami operator and in the generation of pattern formation. Numerical experiments show that the patterns become irregular as the amplitude and frequency of the rough surface increase. For the sake of easy generalization to closed manifolds, we propose a second construction method for rough surfaces, which uses random nodal values and discretized heat filters. We provide numerical evidence that both surface {construction methods} yield comparable patterns to those {observed} in real-life animals.


title: Recovering Dantzig-Wolfe Bounds by Cutting Planes

abstract: Dantzig-Wolfe (DW) decomposition is a well-known technique in mixed-integer programming (MIP) for decomposing and convexifying constraints to obtain potentially strong dual bounds. We investigate cutting planes that can be derived using the DW decomposition algorithm and show that these cuts can provide the same dual bounds as DW decomposition. More precisely, we generate one cut for each DW block, and when combined with the constraints in the original formulation, these cuts imply the objective function cut one can simply write using the DW bound. This approach typically leads to a formulation with lower dual degeneracy that consequently has a better computational performance when solved by standard MIP solvers in the original space. We also discuss how to strengthen these cuts to improve the computational performance further. We test our approach on the Multiple Knapsack Assignment Problem and the Temporal Knapsack Problem, and show that the proposed cuts are helpful in accelerating the solution time without the need to implement branch and price.


title: Nearly Forgotten Cosmological Concept of E. B. Gliner

abstract: E. B. Gliner started his scientific career in 1963 at the age of 40. In 1965, when the existence of the cosmological constant $\lambda$ seemed unnecessary to most cosmologists, he renewed interest in the problem by emphasizing a material interpretation of de Sitter space (i.e., the space curved in the presence of $\lambda$). According to that interpretation, the curvature is produced by a cosmological vacuum (now identified as dark energy of the universe). In 1970, Gliner proposed a description of exponential expansion (or contraction) of the universe at the early (or late) evolution stage dominated by cosmological vacuum. In 1975, Gliner (with I.G. Dyminikova) suggested a model of the early universe free of Big Bang singularity, and developed a scenario of nonsingular Friedmann cosmology. Many of these findings were used in the modern inflation scenarios of the universe, first proposed by A.A. Starobinsky (1979) and A. Guth (1981) and greatly multiplied later. However, these inflation scenarios differ from the scenario of Gliner and Dymnikova, and Gliner's contribution to cosmology is nearly forgotten. The history and the essence of this contribution are outlined, as well the difference from the inflation theories.


title: Convolutional Neural Network-Based Automatic Classification of Colorectal and Prostate Tumor Biopsies Using Multispectral Imagery: System Development Study

abstract: Colorectal and prostate cancers are the most common types of cancer in men worldwide. To diagnose colorectal and prostate cancer, a pathologist performs a histological analysis on needle biopsy samples. This manual process is time-consuming and error-prone, resulting in high intra and interobserver variability, which affects diagnosis reliability. This study aims to develop an automatic computerized system for diagnosing colorectal and prostate tumors by using images of biopsy samples to reduce time and diagnosis error rates associated with human analysis. We propose a CNN model for classifying colorectal and prostate tumors from multispectral images of biopsy samples. The key idea was to remove the last block of the convolutional layers and halve the number of filters per layer. Our results showed excellent performance, with an average test accuracy of 99.8% and 99.5% for the prostate and colorectal data sets, respectively. The system showed excellent performance when compared with pretrained CNNs and other classification methods, as it avoids the preprocessing phase while using a single CNN model for classification. Overall, the proposed CNN architecture was globally the best-performing system for classifying colorectal and prostate tumor images. The proposed CNN was detailed and compared with previously trained network models used as feature extractors. These CNNs were also compared with other classification techniques. As opposed to pretrained CNNs and other classification approaches, the proposed CNN yielded excellent results. The computational complexity of the CNNs was also investigated, it was shown that the proposed CNN is better at classifying images than pretrained networks because it does not require preprocessing. Thus, the overall analysis was that the proposed CNN architecture was globally the best-performing system for classifying colorectal and prostate tumor images.


title: STEEL: Singularity-aware Reinforcement Learning

abstract: Batch reinforcement learning (RL) aims at leveraging pre-collected data to find an optimal policy that maximizes the expected total rewards in a dynamic environment. The existing methods require absolutely continuous assumption (e.g., there do not exist non-overlapping regions) on the distribution induced by target policies with respect to the data distribution over either the state or action or both. We propose a new batch RL algorithm that allows for singularity for both state and action spaces (e.g., existence of non-overlapping regions between offline data distribution and the distribution induced by the target policies) in the setting of an infinite-horizon Markov decision process with continuous states and actions. We call our algorithm STEEL: SingulariTy-awarE rEinforcement Learning. Our algorithm is motivated by a new error analysis on off-policy evaluation, where we use maximum mean discrepancy, together with distributionally robust optimization, to characterize the error of off-policy evaluation caused by the possible singularity and to enable model extrapolation. By leveraging the idea of pessimism and under some technical conditions, we derive a first finite-sample regret guarantee for our proposed algorithm under singularity. Compared with existing algorithms,by requiring only minimal data-coverage assumption, STEEL improves the applicability and robustness of batch RL. In addition, a two-step adaptive STEEL, which is nearly tuning-free, is proposed. Extensive simulation studies and one (semi)-real experiment on personalized pricing demonstrate the superior performance of our methods in dealing with possible singularity in batch RL.


title: $H_{3}S$ is a high-$\kappa$ superconductor with columnar pinning defects

abstract: Recently, the existence of superconductivity in $\mathrm{H_{3}S}$ and other high-pressure hydrides is called into question, because of some flawed magnetic measurements. In Ref.[1] the SCPC model is proposed, where strong pinning of vortices is due to long columnar defects in $\mathrm{H_{3}S}$ - with lengths $\mathrm{L}$ of the order of vortex lengths $\mathrm{L_{v}}$. Two relevant type of experiments in magnetic fields are explained by this model: (1) Reduction (with respect to standard superconductors) of the thermal broadening of resistance ($\mathit{TBR}$) in magnetic field $\mathrm{h=H/H_{c2}}$, $\mathrm{\mathrm{\delta t_{c}(h)}}$, is governed by the small parameter $\mathrm{C=}$$\mathrm{\xi_{0}}/\mathrm{L}$. The $\mathit{TBR}$ measurements give, that in $\mathrm{H_{3}S}$ there is a large shift of the irreversible line $\mathrm{B_{irr}(T)}$ towards the $\mathrm{B_{c2}(T)}$ line. (2) In ZFC (zero field cooled) experiments on penetration of the magnetic field ($\mathit{PMF}$) in $\mathrm{H_{3}S}$ [3], the latter reaches the center of a superconducting disk at much larger external fields, i. e. for $\mathrm{B_{0}\gg B_{c1}}$. The later is due to the pronounced pinning of vortices, but not to the Meissner effect. The calculated $\mathrm{T}$- dependence of the penetrated, perpendicular $\mathrm{B_{\perp}(T)}$ and parallel $\mathrm{B_{\parallel}}(\mathrm{T})$, magnetic field into the sample is in satisfactory agreement with the experimental results for $\mathrm{H_{3}S}$ in [3], where $\mathrm{B_{\perp}>B_{\parallel}}$. Problems related to measurements of the Meissner effect in $\mathrm{H_{3}S}$ are also discussed. The SCPC model, applied on the bulk $\mathrm{H_{3}S}$ sample, predicts, that the latter is a high-$\kappa$ superconductor.


title: Protein Representation Learning via Knowledge Enhanced Primary Structure Modeling

abstract: Protein representation learning has primarily benefited from the remarkable development of language models (LMs). Accordingly, pre-trained protein models also suffer from a problem in LMs: a lack of factual knowledge. The recent solution models the relationships between protein and associated knowledge terms as the knowledge encoding objective. However, it fails to explore the relationships at a more granular level, i.e., the token level. To mitigate this, we propose Knowledge-exploited Auto-encoder for Protein (KeAP), which performs token-level knowledge graph exploration for protein representation learning. In practice, non-masked amino acids iteratively query the associated knowledge tokens to extract and integrate helpful information for restoring masked amino acids via attention. We show that KeAP can consistently outperform the previous counterpart on 9 representative downstream applications, sometimes surpassing it by large margins. These results suggest that KeAP provides an alternative yet effective way to perform knowledge enhanced protein representation learning.


title: Advancing Radiograph Representation Learning with Masked Record Modeling

abstract: Modern studies in radiograph representation learning rely on either self-supervision to encode invariant semantics or associated radiology reports to incorporate medical expertise, while the complementarity between them is barely noticed. To explore this, we formulate the self- and report-completion as two complementary objectives and present a unified framework based on masked record modeling (MRM). In practice, MRM reconstructs masked image patches and masked report tokens following a multi-task scheme to learn knowledge-enhanced semantic representations. With MRM pre-training, we obtain pre-trained models that can be well transferred to various radiography tasks. Specifically, we find that MRM offers superior performance in label-efficient fine-tuning. For instance, MRM achieves 88.5% mean AUC on CheXpert using 1% labeled data, outperforming previous R$^2$L methods with 100% labels. On NIH ChestX-ray, MRM outperforms the best performing counterpart by about 3% under small labeling ratios. Besides, MRM surpasses self- and report-supervised pre-training in identifying the pneumonia type and the pneumothorax area, sometimes by large margins.


title: SeaFormer++: Squeeze-enhanced Axial Transformer for Mobile Visual Recognition

abstract: Since the introduction of Vision Transformers, the landscape of many computer vision tasks (e.g., semantic segmentation), which has been overwhelmingly dominated by CNNs, recently has significantly revolutionized. However, the computational cost and memory requirement renders these methods unsuitable on the mobile device. In this paper, we introduce a new method squeeze-enhanced Axial Transformer (SeaFormer) for mobile visual recognition. Specifically, we design a generic attention block characterized by the formulation of squeeze Axial and detail enhancement. It can be further used to create a family of backbone architectures with superior cost-effectiveness. Coupled with a light segmentation head, we achieve the best trade-off between segmentation accuracy and latency on the ARM-based mobile devices on the ADE20K, Cityscapes, Pascal Context and COCO-Stuff datasets. Critically, we beat both the mobilefriendly rivals and Transformer-based counterparts with better performance and lower latency without bells and whistles. Furthermore, we incorporate a feature upsampling-based multi-resolution distillation technique, further reducing the inference latency of the proposed framework. Beyond semantic segmentation, we further apply the proposed SeaFormer architecture to image classification and object detection problems, demonstrating the potential of serving as a versatile mobile-friendly backbone. Our code and models are made publicly available at https://github.com/fudan-zvg/SeaFormer.


title: Rank 2 $\ell$-adic local systems and Higgs bundles over a curve

abstract: Let $X$ be a smooth, projective, and geometrically connected curve defined over a finite field $\mathbb{F}_q$ and $S\subseteq X$ a subset of closed points. Let $\bar{X}$ and $\bar{S}$ be their base changes to an algebraic closure of $\mathbb{F}_q$. We study the number of $\ell$-adic local systems in rank $2$ over $\bar{X}-\bar{S}$ with prescribed tame local monodromies fixed by $k$-fold iterated action of Frobenius endomorphism for every $k\geq 1$. We confirm some conjectures of Deligne predicting that these numbers behave as if they were obtained from a Lefschetz fixed point formula. In fact, in all cases, our counting results are expressed in terms of the numbers of some Higgs bundles.


title: Selection principles and proofs from the Book

abstract: I provide simplified proofs for each of the following fundamental theorems regarding selection principles: 1. The Quasinormal Convergence Theorem, due to the author and Zdomskyy, asserting that a certain, important property of the space of continuous functions on a space is actually preserved by Borel images of that space. 2. The Scheepers Diagram Last Theorem, due to Peng, completing all provable implications in the diagram. 3. The Menger Game Theorem, due to Telg\'arsky, determining when Bob has a winning strategy in the game version of Menger's covering property. 4. A lower bound on the additivity of Rothberger's covering property, due to Carlson. The simplified proofs lead to several new results.


title: Spectral properties of the Laplacian of temporal networks following a constant block Jacobi model

abstract: We study the behavior of the eigenvectors associated with the smallest eigenvalues of the Laplacian matrix of temporal networks. We consider the multilayer representation of temporal networks, i.e. a set of networks linked through ordinal interconnected layers. We analyze the Laplacian matrix, known as supra-Laplacian, constructed through the supra-adjacency matrix associated with the multilayer formulation of temporal networks, using a constant block Jacobi model which has closed-form solution. To do this, we assume that the inter-layer weights are perturbations of the Kronecker sum of the separate adjacency matrices forming the temporal network. Thus we investigate the properties of the eigenvectors associated with the smallest eigenvalues (close to zero) of the supra-Laplacian matrix. Using arguments of perturbation theory, we show that these eigenvectors can be approximated by linear combinations of the zero eigenvectors of the individual time layers. This finding is crucial in reconsidering and generalizing the role of the Fielder vector in supra-Laplacian matrices.


title: Mathematical modelling and numerical simulation of reverse-osmosis desalination

abstract: The reverse osmosis membrane module is an integral element of a desalination system as it determines the overall performance of the desalination plant. The fraction of clean water that can be recovered via this process is often limited by salt precipitation which plays a critical role in its sustainability. In this work, we present a model to study the complex interplay between flow, transport and precipitation processes in reverse osmosis membranes, which together influence recovery and in turn process sustainability. A reactive porous interface model describes the membrane with a dynamic evolving porosity and permeability to capture the scaling and clogging of the membrane. An open-source finite-volume numerical solver is implemented within the OpenFOAM library and numerical tests are presented here showing the effect of the various parameters of the model and the robustness of the model to describe a wide range of operating conditions.

categories: math.NA cs.NA physics.chem-ph physics.flu-dyn


title: Echoes of the hexagon: remnants of hexagonal packing inside regular polygons

abstract: Based on numerical simulations that we have carried out, we provide evidence that for regular polygons with $\sigma= 6j$ sides (with $j=2,3,\dots$), $N(k)=3 k (k+1)+1$ (with $k=1,2,\dots$) congruent disks of appropriate size can be nicely packed inside these polygons in highly symmetrical configurations which apparently have maximal density for $N$ sufficiently small. These configurations are invariant under rotations of $\pi/3$ and are closely related to the configurations with perfect hexagonal packing in the regular hexagon and to the configurations with {\sl curved hexagonal packing} (CHP) in the circle found long time ago by Graham and Lubachevsky. At the basis of our explorations are the algorithms that we have devised, which are very efficient in producing the CHP and more general configurations inside regular polygons. We have used these algorithms to generate a large number of CHP configurations for different regular polygons and numbers of disks; a careful study of these results has made possible to fully characterize the general properties of the CHP configurations and to devise a {\sl deterministic} algorithm that completely ensembles a given CHP configuration once an appropriate input ("DNA") is specified. Our analysis shows that the number of CHP configurations for a given $N$ is highly degenerate in the packing fraction and it can be explicitly calculated in terms of $k$ (number of shells), of the building block of the DNA itself and of the number of vertices in the fundamental domain (because of the symmetry we work in $1/6$ of the whole domain). With the help of our deterministic algorithm we are able to build {\sl all} the CHP configurations for a polygon with $k$ shells.


title: A deep learning approach for adaptive zoning

abstract: We propose a supervised deep learning (DL) approach to perform adaptive zoning on time dependent partial differential equations that model the propagation of 1D shock waves in a compressible medium. We train a neural network on a dataset composed of different static shock profiles associated with the corresponding adapted meshes computed with standard adaptive zoning techniques. We show that the trained DL model learns how to capture the presence of shocks in the domain and generates at each time step an adapted non-uniform mesh that relocates the grid nodes to improve the accuracy of Lax-Wendroff and fifth order weighted essentially non-oscillatory (WENO5) space discretization schemes. We also show that the surrogate DL model reduces the computational time to perform adaptive zoning by at least a 2x factor with respect to standard techniques without compromising the accuracy of the reconstruction of the physical quantities of interest.


title: Randomized GCUR decompositions

abstract: By exploiting the random sampling techniques, this paper derives an efficient randomized algorithm for computing a generalized CUR decomposition, which provides low-rank approximations of both matrices simultaneously in terms of some of their rows and columns. For large-scale data sets that are expensive to store and manipulate, a new variant of the discrete empirical interpolation method known as L-DEIM, which needs much lower cost and provides a significant acceleration in practice, is also combined with the random sampling approach to further improve the efficiency of our algorithm. Moreover, adopting the randomized algorithm to implement the truncation process of restricted singular value decomposition (RSVD), combined with the L-DEIM procedure, we propose a fast algorithm for computing an RSVD based CUR decomposition, which provides a coordinated low-rank approximation of the three matrices in a CUR-type format simultaneously and provides advantages over the standard CUR approximation for some applications. We establish detailed probabilistic error analysis for the algorithms and provide numerical results that show the promise of our approaches.


title: A Rellich's result revisited and sensitivity of solutions of parametrized linear systems

abstract: In this paper we revisit a result due to Franz Rellich on smoothness of solutions of parametrized linear systems. With this result as a starting point, we obtain finer smoothness results in an elementary fashion and propose an efficient adjoint algorithm for computing sensitivities of $(n-1)$-deficient systems, being $n$ the order of the system.


title: Label-free learning of elliptic partial differential equation solvers with generalizability across boundary value problems

abstract: Traditional, numerical discretization-based solvers of partial differential equations (PDEs) are fundamentally agnostic to domains, boundary conditions and coefficients. In contrast, machine learnt solvers have a limited generalizability across these elements of boundary value problems. This is strongly true in the case of surrogate models that are typically trained on direct numerical simulations of PDEs applied to one specific boundary value problem. In a departure from this direct approach, the label-free machine learning of solvers is centered on a loss function that incorporates the PDE and boundary conditions in residual form. However, their generalization across boundary conditions is limited and they remain strongly domain-dependent. Here, we present a framework that generalizes across domains, boundary conditions and coefficients simultaneously with learning the PDE in weak form. Our work explores the ability of simple, convolutional neural network (CNN)-based encoder-decoder architectures to learn to solve a PDE in greater generality than its restriction to a particular boundary value problem. In this first communication, we consider the elliptic PDEs of Fickean diffusion, linear and nonlinear elasticity. Importantly, the learning happens independently of any labelled field data from either experiments or direct numerical solutions. Extensive results for these problem classes demonstrate the framework's ability to learn PDE solvers that generalize across hundreds of thousands of domains, boundary conditions and coefficients, including extrapolation beyond the learning regime. Once trained, the machine learning solvers are orders of magnitude faster than discretization-based solvers. We place our work in the context of recent continuous operator learning frameworks, and note extensions to transfer learning, active learning and reinforcement learning.


title: ESC: Exploration with Soft Commonsense Constraints for Zero-shot Object Navigation

abstract: The ability to accurately locate and navigate to a specific object is a crucial capability for embodied agents that operate in the real world and interact with objects to complete tasks. Such object navigation tasks usually require large-scale training in visual environments with labeled objects, which generalizes poorly to novel objects in unknown environments. In this work, we present a novel zero-shot object navigation method, Exploration with Soft Commonsense constraints (ESC), that transfers commonsense knowledge in pre-trained models to open-world object navigation without any navigation experience nor any other training on the visual environments. First, ESC leverages a pre-trained vision and language model for open-world prompt-based grounding and a pre-trained commonsense language model for room and object reasoning. Then ESC converts commonsense knowledge into navigation actions by modeling it as soft logic predicates for efficient exploration. Extensive experiments on MP3D, HM3D, and RoboTHOR benchmarks show that our ESC method improves significantly over baselines, and achieves new state-of-the-art results for zero-shot object navigation (e.g., 288% relative Success Rate improvement than CoW on MP3D).


title: A complete classification of categoricity spectra of accessible categories with directed colimits

abstract: We provide a complete classification of all the possible categoricity spectra, in terms of internal size, that can appear in a large accessible category with directed colimits, assuming the Singular Cardinal Hypothesis ($SCH$), and providing as well explicit threshold cardinals for eventual categoricity. This includes as a particular case the first complete classification of categoricity spectra of abstract elementary classes (AEC's) entirely in $ZFC$. More specifically, we have the following theorem: Let $\mathcal{K}$ be a large $\kappa$-accessible category with directed colimits. Assume the Singular Cardinal Hypothesis $SCH$ (only if the restriction to monomorphisms is not an AEC). Then the categoricity spectrum $\mathcal{C}at(\mathcal{K})=\{\lambda\geq \kappa: \mathcal{K} \text{ is $\lambda$-categorical}\}$ is one of the following: 1) $\mathcal{C}at(\mathcal{K})=\emptyset$. 2) $\mathcal{C}at(\mathcal{K})=[\alpha, \beta]$ for some $\alpha, \beta \in [\kappa, \beth_{\omega}(\kappa))$. 3) $\mathcal{C}at(\mathcal{K})=[\chi, \infty)$ for some $\chi \in [\kappa, \beth_{(2^{\kappa})^+})$. This solves in particular Shelah categoricity conjecture for AEC's. There are examples of each of the three cases of the classification, showing that they indeed occur.


title: The noncommutative minimal model program

abstract: This note aims to clarify the deep relationship between birational modifications of a variety and semiorthogonal decompositions of its derived category of coherent sheaves. The result is a conjecture on the existence and properties of canonical semiorthogonal decompositions, which is a noncommutative analog of the minimal model program. We identify a mechanism for constructing semiorthogonal decompositions using Bridgeland stability conditions, and we propose that through this mechanism the quantum differential equation of the variety controls the conjectured semiorthogonal decompositions. We establish several implications of the conjectures: one direction of Dubrovin's conjecture on the existence of full exceptional collections; the $D$-equivalence conjecture; the existence of new categorical birational invariants for varieties of positive genus; and the existence of minimal noncommutative resolutions of singular varieties. Finally, we verify the conjectures for smooth projective curves by establishing a previously conjectured description of the stability manifold of $\mathbb{P}^1$.


title: Improved machine learning algorithm for predicting ground state properties

abstract: Finding the ground state of a quantum many-body system is a fundamental problem in quantum physics. In this work, we give a classical machine learning (ML) algorithm for predicting ground state properties with an inductive bias encoding geometric locality. The proposed ML model can efficiently predict ground state properties of an $n$-qubit gapped local Hamiltonian after learning from only $\mathcal{O}(\log(n))$ data about other Hamiltonians in the same quantum phase of matter. This improves substantially upon previous results that require $\mathcal{O}(n^c)$ data for a large constant $c$. Furthermore, the training and prediction time of the proposed ML model scale as $\mathcal{O}(n \log n)$ in the number of qubits $n$. Numerical experiments on physical systems with up to 45 qubits confirm the favorable scaling in predicting ground state properties using a small training dataset.


title: Hamiltonian-Oriented Homotopy QAOA

abstract: The classical homotopy optimization approach has the potential to deal with highly nonlinear landscape, such as the energy landscape of QAOA problems. Following this motivation, we introduce Hamiltonian-Oriented Homotopy QAOA (HOHo-QAOA), that is a heuristic method for combinatorial optimization using QAOA, based on classical homotopy optimization. The method consists of a homotopy map that produces an optimization problem for each value of interpolating parameter. Therefore, HOHo-QAOA decomposes the optimization of QAOA into several loops, each using a mixture of the mixer and the objective Hamiltonian for cost function evaluation. Furthermore, we conclude that the HOHo-QAOA improves the search for low energy states in the nonlinear energy landscape and outperforms other variants of QAOA.


title: Lensing constraints on ultradense dark matter halos

abstract: Cosmological observations precisely measure primordial variations in the density of the Universe at megaparsec and larger scales, but much smaller scales remain poorly constrained. However, sufficiently large initial perturbations at small scales can lead to an abundance of ultradense dark matter minihalos that form during the radiation epoch and survive into the late-time Universe. Because of their early formation, these objects can be compact enough to produce detectable microlensing signatures. We investigate whether the EROS, OGLE, and HSC surveys can probe these halos by fully accounting for finite source size and extended lens effects. We find that current data may already constrain the amplitudes of primordial curvature perturbations in a new region of parameter space, but this conclusion is strongly sensitive to yet undetermined details about the internal structures of these ultradense halos. Under optimistic assumptions, current and future HSC data would constrain a power spectrum that features an enhancement at scales $k \sim 10^7/{\rm Mpc}$, and an amplitude as low as $\mathcal{P}_\zeta\simeq 10^{-4}$ may be accessible. This is a particularly interesting regime because it connects to primordial black hole formation in a portion of the LIGO/Virgo/Kagra mass range and the production of scalar-induced gravitational waves in the nanohertz frequency range reachable by pulsar timing arrays. These prospects motivate further study of the ultradense halo formation scenario to clarify their internal structures.


title: Cell Systems for $\overline{\operatorname{Rep}(U_q(\mathfrak{sl}_N))}$ Module Categories

abstract: In this paper, we define the KW cell system on a graph $\Gamma$, depending on parameters $N\in \mathbb{N}$, $q$ a root of unity, and $\omega$ an $N$-th root of unity. This is a polynomial system of equations depending on $\Gamma$ and the parameters. Using the graph planar algebra embedding theorem, we prove that when $q = e^{2\pi i \frac{1}{2(N+k)}}$, solutions to the KW cell system on $\Gamma$ classify module categories over $\overline{\mathrm{Rep}(U_q(sl_N))^\omega}$ whose action graph for the object $\Lambda_1$ is $\Gamma$. The KW cell system is a generalisation of the Etingof-Ostrik and the De Commer-Yamashita classifying data for $\overline{\mathrm{Rep}(U_q(sl_2))}$ module categories, and Ocneanu's cell calculus for $\overline{\mathrm{Rep}(U_q(sl_3))}$ module categories. To demonstrate the effectiveness of this cell calculus, we solve the KW cell systems corresponding to the exceptional module categories over $\overline{\mathrm{Rep}(U_q(sl_4))}$ when $q= e^{2\pi i \frac{1}{2(4+k)}}$, as well as for all three infinite families of charge conjugation modules. Building on the work of the second author, this explicitly constructs and classifies all irreducible module categories over $\mathcal{C}(sl_4, k)$ for all $k\in \mathbb{N}$. These results prove claims made by Ocneanu on the quantum subgroups of $SU(4)$. We also construct exceptional module categories over $\overline{\mathrm{Rep}(U_q(sl_4))^\omega}$ where $\omega\in \{-1, i, -i\}$. Two of these module categories have no analogue when $\omega=1$. The main technical contributions of this paper are a proof of the graph planar algebra embedding theorem for oriented planar algebras, and a refinement of Kazhdan and Wenzl's skein theory presentation of the category $\overline{\mathrm{Rep}(U_q(sl_N))^\omega}$. We also explicitly describe the subfactors coming from a solution to a KW cell system.


title: Shape-aware Text-driven Layered Video Editing

abstract: Temporal consistency is essential for video editing applications. Existing work on layered representation of videos allows propagating edits consistently to each frame. These methods, however, can only edit object appearance rather than object shape changes due to the limitation of using a fixed UV mapping field for texture atlas. We present a shape-aware, text-driven video editing method to tackle this challenge. To handle shape changes in video editing, we first propagate the deformation field between the input and edited keyframe to all frames. We then leverage a pre-trained text-conditioned diffusion model as guidance for refining shape distortion and completing unseen regions. The experimental results demonstrate that our method can achieve shape-aware consistent video editing and compare favorably with the state-of-the-art.


title: Future of Supply Chain: Challenges, Trends, and Prospects

abstract: This paper discusses the broad challenges shared by e-commerce and the process industries operating global supply chains. Specifically, we discuss how process industries and e-commerce differ in many aspects but have similar challenges ahead of them in order to remain competitive, keep up with the always increasing requirements of the customers and stakeholders, and gain efficiency. While both industries have been early adopters of decision support tools based on machine intelligence, both share unresolved challenges related to scalability, integration of decision-making over different time horizons (e.g. strategic, tactical and execution-level decisions) and across internal business units, and orchestration of human and computer-based decision-makers. We discuss future trends and research opportunities in the area of supply chain, and suggest that the methods of multi-agent systems supported by rigorous treatment of human decision-making in combination with machine intelligence is a great contender to address these critical challenges.


title: Cops and robbers on $P_5$-free graphs

abstract: We prove that every connected $P_5$-free graph has cop number at most two, solving a conjecture of Sivaraman. In order to do so, we first prove that every connected $P_5$-free graph $G$ with independence number at least three contains a three-vertex induced path with vertices $a \hbox{-} b \hbox{-} c$ in order, such that every neighbour of $c$ is also adjacent to one of $a,b$.


title: Graphene Oxide Photoreduction Recovers Graphene Hot Electron Cooling Dynamics

abstract: Reduced graphene oxide (rGO) is a bulk-processable quasi-amorphous 2D material with broad spectral coverage and fast electronic response. rGO sheets are suspended in a polymer matrix and sequentially photoreduced while measuring the evolving optical spectra and ultrafast electron relaxation dynamics. Photoreduced rGO yields optical absorption spectra that fit with the same Fano lineshape parameters as monolayer graphene. With increasing photoreduction time, rGO transient absorption kinetics accelerate monotonically, reaching an optimal point that matches the hot electron cooling in graphene. All stages of rGO ultrafast kinetics are simulated with a hot-electron cooling model mediated by disorder-assisted supercollisions. While the rGO room temperature 0.31 ps$^{-1}$ electronic cooling rate matches monolayer graphene, subsequent photoreduction can rapidly increase the rate by ~10-12$\times$. Such accelerated supercollision rates imply a reduced mean-free scattering length caused by photoionized point-defects on the rGO sp$^2$ sub-lattice. For visible range excitations of rGO, photoreduction shows three increasing spectral peaks that match graphene quantum dot (GQD) transitions, while a broad peak from oxygenated defect edge states shrinks. These three confined GQD states donate their hot carriers to the graphene sub-lattice with a 0.17 ps rise-time that accelerates with photoreduction. Collectively, many desirable photophysical properties of 2D graphene are replicated through selectively reducing rGO scaffolded within a 3D bulk polymeric network.


title: Infinite-Variate $L^2$-Approximation with Nested Subspace Sampling

abstract: We consider $L^2$-approximation on weighted reproducing kernel Hilbert spaces of functions depending on infinitely many variables. We focus on unrestricted linear information, admitting evaluations of arbitrary continuous linear functionals. We distinguish between ANOVA and non-ANOVA spaces, where, by ANOVA spaces, we refer to function spaces whose norms are induced by an underlying ANOVA function decomposition. In ANOVA spaces, we provide an optimal algorithm to solve the approximation problem using linear information. We determine the upper and lower error bounds on the polynomial convergence rate of $n$-th minimal worst-case errors, which match if the weights decay regularly. For non-ANOVA spaces, we also establish upper and lower error bounds. Our analysis reveals that for weights with a regular and moderate decay behavior, the convergence rate of $n$-th minimal errors is strictly higher in ANOVA than in non-ANOVA spaces.


title: Accurate modelling of the Lyman-$\alpha$ coupling for the 21-cm signal, observability with NenuFAR and SKA

abstract: The measurement of the $21$ cm signal from the Cosmic Dawn is a major goal for several existing and upcoming radio interferometers such as NenuFAR and the SKA. During this era before the beginning of the Epoch of Reionization, the signal is more difficult to observe due to brighter foregrounds but reveals additional information on the underlying astrophysical processes encoded in the spatial fluctuations of the spin temperature of hydrogen. To interpret future measurements, controlling the level of accuracy of the Lyman-$\alpha$ flux modelling is mandatory. In this work, we evaluate the impact of various approximations that exist in the main fast modelling approach compared to the results of a costly full radiative transfer simulation. The fast SPINTER code, presented in this work, computes the Lyman-$\alpha$ flux including the effect of wing scatterings for an inhomogeneous emissivity field, but assuming an otherwise homogeneous expanding universe. The LICORICE code computes the full radiative transfer in the Lyman-$\alpha$ line without any substantial approximation. We find that the difference between homogeneous and inhomogeneous gas density and temperature is very small for the computed flux. On the contrary, neglecting the effect of gas velocities produces a significant change in the computed flux. We identify the causes (mainly Doppler shifts due to velocity gradients) and quantify the magnitude of the effect in both an idealised setup and a realistic cosmological situation. We find that the amplitude of the effect, up to a factor of $\sim 2$ on the $21$ cm signal power spectrum on some scales (depending on both other model parameters and the redshift), can be easily discriminated with an SKA-like survey and already be approached, particularly for exotic signals, by the ongoing NenuFAR Cosmic Dawn Key Science Program.


title: Viscoelastic response of an active particle under the action of magnetic field

abstract: We consider the dynamics of a charged inertial active Ornstein-Uhlenbeck particle in a viscoelastic suspension under the action of an uniform magnetic field. With the help of both numerical simulation and analytical framework, we exactly investigate the viscoelastic response of the particle to the magnetic field by means of particle trajectories, mean displacement, and mean square displacement (MSD) calculations. The simulated particle trajectories and MSD calculations reveal that the steady state response of a confined harmonic particle to the magnetic field show interesting features due to the complex interplay of underlying physical processes such as elastic dissipation and active fluctuations in the medium. When the activity time scale of the dynamics is larger than the elastic dissipation timescale (or vice-versa), the steady state MSD shows an enhancement (or suppression) with increase in the field strength. In both the cases, for very strong magnetic field the MSD interestingly approaches the value at the equilibrium limit where the generalized fluctuation dissipation relation is satisfied and the particle behaves like an inertial passive Brownian particle. Thus, for a finite magnetic field, the system exhibits a re-entrant type transition from active to passive and then to active behaviour with persistence of elastic dissipation in the medium. However, for a free particle, the overall displacement covered always gets suppressed with increase in the strength of the magnetic field and finally approaches zero for a very strong magnetic field.


title: Massive Gravity and Lorentz Symmetry

abstract: We consider Lorentz-symmetry properties of the ghost-free massive gravity theory proposed by de Rham, Gabadadze, and Tolley. In particular, we present potentially observable effects in gravitational-wave propagation and in Newton's law, including Lorentz-violating signals.


title: Partitioned Matching Games for International Kidney Exchange

abstract: We introduce partitioned matching games as a suitable model for international kidney exchange programmes, where in each round the total number of available kidney transplants needs to be distributed amongst the participating countries in a "fair" way. A partitioned matching game $(N,v)$ is defined on a graph $G=(V,E)$ with an edge weighting $w$ and a partition $V=V_1 \cup \dots \cup V_n$. The player set is $N = \{1, \dots, n\}$, and player $p \in N$ owns the vertices in $V_p$. The value $v(S)$ of a coalition $S \subseteq N$ is the maximum weight of a matching in the subgraph of $G$ induced by the vertices owned by the players in $S$. If $|V_p|=1$ for all $p\in N$, then we obtain the classical matching game. Let $c=\max\{|V_p| \; |\; 1\leq p\leq n\}$ be the width of $(N,v)$. We prove that checking core non-emptiness is polynomial-time solvable if $c\leq 2$ but co-NP-hard if $c\leq 3$. We do this via pinpointing a relationship with the known class of $b$-matching games and completing the complexity classification on testing core non-emptiness for $b$-matching games. With respect to our application, we prove a number of complexity results on choosing, out of possibly many optimal solutions, one that leads to a kidney transplant distribution that is as close as possible to some prescribed fair distribution.


title: Open string field theory with stubs

abstract: There are various reasons why adding stubs to the vertices of open string field theory (OSFT) is interesting: Not only the stubs can tame certain singularities and make the theory more well-behaved, but also the new theory shares a lot of similarities with closed string field theory, which helps to improve our understanding of its structure and possible solutions. In this paper we explore two natural ways of implementing stubs into the framework of OSFT, resulting in an A-infinity-algebra giving rise to infinitely many vertices. We find two distinct consistent actions, both generated by a field redefinition, interestingly sharing the same equations of motion. In the last section we illustrate their relationship and physical meaning by applying our construction to nearly marginal solutions.


title: Learning Control from Raw Position Measurements

abstract: We propose a Model-Based Reinforcement Learning (MBRL) algorithm named VF-MC-PILCO, specifically designed for application to mechanical systems where velocities cannot be directly measured. This circumstance, if not adequately considered, can compromise the success of MBRL approaches. To cope with this problem, we define a velocity-free state formulation which consists of the collection of past positions and inputs. Then, VF-MC-PILCO uses Gaussian Process Regression to model the dynamics of the velocity-free state and optimizes the control policy through a particle-based policy gradient approach. We compare VF-MC-PILCO with our previous MBRL algorithm, MC-PILCO4PMS, which handles the lack of direct velocity measurements by modeling the presence of velocity estimators. Results on both simulated (cart-pole and UR5 robot) and real mechanical systems (Furuta pendulum and a ball-and-plate rig) show that the two algorithms achieve similar results. Conveniently, VF-MC-PILCO does not require the design and implementation of state estimators, which can be a challenging and time-consuming activity to be performed by an expert user.


title: Evolution of the Thermodynamic Properties of a Coronal Mass Ejection in the Inner Corona

abstract: The thermodynamic evolution of Coronal Mass Ejections (CMEs) in the inner corona (< 1.5 R$_{sun}$) is not yet completely understood. In this work, we study the evolution of thermodynamic properties of a CME core observed in the inner corona on July 20, 2017, by combining the MLSO/K-Cor white-light and the MLSO/CoMP Fe XIII 10747 {\AA} line spectroscopic data. We also estimate the emission measure weighted temperature (T$_{EM}$) of the CME core by applying the Differential Emission Measure (DEM) inversion technique on the SDO/AIA six EUV channels data and compare it with the effective temperature (T$_{eff}$) obtained using Fe XIII line width measurements. We find that the T$_{eff}$ and T$_{EM}$ of the CME core show similar variation and remain almost constant as the CME propagates from ~1.05 to 1.35 R$_{sun}$. The temperature of the CME core is of the order of million-degree kelvin, indicating that it is not associated with a prominence. Further, we estimate the electron density of this CME core using K-Cor polarized brightness (pB) data and found it decreasing by a factor of ~ 3.6 as the core evolves. An interesting finding is that the temperature of the CME core remains almost constant despite expected adiabatic cooling due to the expansion of the CME core, which suggests that the CME core plasma must be heated as it propagates. We conclude that the expansion of this CME core behaves more like an isothermal than an adiabatic process.


title: Optimal Decision Tree Policies for Markov Decision Processes

abstract: Interpretability of reinforcement learning policies is essential for many real-world tasks but learning such interpretable policies is a hard problem. Particularly rule-based policies such as decision trees and rules lists are difficult to optimize due to their non-differentiability. While existing techniques can learn verifiable decision tree policies there is no guarantee that the learners generate a decision that performs optimally. In this work, we study the optimization of size-limited decision trees for Markov Decision Processes (MPDs) and propose OMDTs: Optimal MDP Decision Trees. Given a user-defined size limit and MDP formulation OMDT directly maximizes the expected discounted return for the decision tree using Mixed-Integer Linear Programming. By training optimal decision tree policies for different MDPs we empirically study the optimality gap for existing imitation learning techniques and find that they perform sub-optimally. We show that this is due to an inherent shortcoming of imitation learning, namely that complex policies cannot be represented using size-limited trees. In such cases, it is better to directly optimize the tree for expected return. While there is generally a trade-off between the performance and interpretability of machine learning models, we find that OMDTs limited to a depth of 3 often perform close to the optimal limit.


title: Accurate Gaze Estimation using an Active-gaze Morphable Model

abstract: Rather than regressing gaze direction directly from images, we show that adding a 3D shape model can: i) improve gaze estimation accuracy, ii) perform well with lower resolution inputs and iii) provide a richer understanding of the eye-region and its constituent gaze system. Specifically, we use an `eyes and nose' 3D morphable model (3DMM) to capture the eye-region 3D facial geometry and appearance and we equip this with a geometric vergence model of gaze to give an `active-gaze 3DMM'. We show that our approach achieves state-of-the-art results on the Eyediap dataset and we present an ablation study. Our method can learn with only the ground truth gaze target point and the camera parameters, without access to the ground truth gaze origin points, thus widening the applicability of our approach compared to other methods.


title: Weighted Flow Diffusion for Local Graph Clustering with Node Attributes: an Algorithm and Statistical Guarantees

abstract: Local graph clustering methods aim to detect small clusters in very large graphs without the need to process the whole graph. They are fundamental and scalable tools for a wide range of tasks such as local community detection, node ranking and node embedding. While prior work on local graph clustering mainly focuses on graphs without node attributes, modern real-world graph datasets typically come with node attributes that provide valuable additional information. We present a simple local graph clustering algorithm for graphs with node attributes, based on the idea of diffusing mass locally in the graph while accounting for both structural and attribute proximities. Using high-dimensional concentration results, we provide statistical guarantees on the performance of the algorithm for the recovery of a target cluster with a single seed node. We give conditions under which a target cluster generated from a fairly general contextual random graph model, which includes both the stochastic block model and the planted cluster model as special cases, can be fully recovered with bounded false positives. Empirically, we validate all theoretical claims using synthetic data, and we show that incorporating node attributes leads to superior local clustering performances using real-world graph datasets.


title: Extracting Training Data from Diffusion Models

abstract: Image diffusion models such as DALL-E 2, Imagen, and Stable Diffusion have attracted significant attention due to their ability to generate high-quality synthetic images. In this work, we show that diffusion models memorize individual images from their training data and emit them at generation time. With a generate-and-filter pipeline, we extract over a thousand training examples from state-of-the-art models, ranging from photographs of individual people to trademarked company logos. We also train hundreds of diffusion models in various settings to analyze how different modeling and data decisions affect privacy. Overall, our results show that diffusion models are much less private than prior generative models such as GANs, and that mitigating these vulnerabilities may require new advances in privacy-preserving training.


title: Localised graph Maclaurin inequalities

abstract: The Maclaurin inequalities for graphs are a broad generalisation of the classical theorems of Tur\'an and Zykov. In a nutshell they provide an asymptotically sharp answer to the following question: what is the maximum number of cliques of size $q$ in a $K_{r+1}$-free graph with a given number of cliques of size $s$? We prove an extensions of the graph Maclaurin inequalities with a weight function that captures the local structure of the graph. As a corollary, we settle a recent conjecture of Kirsch and Nir, which simultaneously encompass the previous localised results of Brada\v{c}, Malec and Tompkins and of Kirsch and Nir.


title: Audio-Visual Segmentation with Semantics

abstract: We propose a new problem called audio-visual segmentation (AVS), in which the goal is to output a pixel-level map of the object(s) that produce sound at the time of the image frame. To facilitate this research, we construct the first audio-visual segmentation benchmark, i.e., AVSBench, providing pixel-wise annotations for sounding objects in audible videos. It contains three subsets: AVSBench-object (Single-source subset, Multi-sources subset) and AVSBench-semantic (Semantic-labels subset). Accordingly, three settings are studied: 1) semi-supervised audio-visual segmentation with a single sound source; 2) fully-supervised audio-visual segmentation with multiple sound sources, and 3) fully-supervised audio-visual semantic segmentation. The first two settings need to generate binary masks of sounding objects indicating pixels corresponding to the audio, while the third setting further requires generating semantic maps indicating the object category. To deal with these problems, we propose a new baseline method that uses a temporal pixel-wise audio-visual interaction module to inject audio semantics as guidance for the visual segmentation process. We also design a regularization loss to encourage audio-visual mapping during training. Quantitative and qualitative experiments on AVSBench compare our approach to several existing methods for related tasks, demonstrating that the proposed method is promising for building a bridge between the audio and pixel-wise visual semantics. Code is available at https://github.com/OpenNLPLab/AVSBench. Online benchmark is available at http://www.avlbench.opennlplab.cn.


title: Minkowski Functionals in $SO(3)$ for the spin-2 CMB polarisation field

abstract: The study of the angular power spectrum of Cosmic Microwave Background (CMB) anisotropies, both in intensity and in polarisation, has led to the tightest constraints on cosmological parameters. However, this statistical quantity is not sensitive to any deviation from Gaussianity and statistical isotropy in the CMB data. Minkowski Functionals (MFs) have been adopted as one of the most powerful statistical tools to study such deviations, since they characterise the topology and geometry of the field of interest. In this paper, we extend the application of MFs to CMB polarisation data by introducing a new formalism, where we lift the spin $2$ polarisation field to a scalar function in a higher-dimensional manifold: the group of rotations of the sphere, $SO(3)$. Such a function is defined as $f = Q \cos(2\psi) - U \sin(2\psi)$. We analytically obtain the expected values for the MFs of $f$ in the case of Gaussian isotropic polarisation maps. Furthermore, we present a new pipeline which estimates these MFs from input HEALPix polarisation maps. We apply it to CMB simulations in order to validate the theoretical results and the methodology. The pipeline is to be included in the publicly available Python package $\texttt{Pynkowski}$ available at https://github.com/javicarron/pynkowski.


title: Robust empirical risk minimization via Newton's method

abstract: A new variant of Newton's method for empirical risk minimization is studied, where at each iteration of the optimization algorithm, the gradient and Hessian of the objective function are replaced by robust estimators taken from existing literature on robust mean estimation for multivariate data. After proving a general theorem about the convergence of successive iterates to a small ball around the population-level minimizer, consequences of the theory in generalized linear models are studied when data are generated from Huber's epsilon-contamination model and/or heavytailed distributions. An algorithm for obtaining robust Newton directions based on the conjugate gradient method is also proposed, which may be more appropriate for high-dimensional settings, and conjectures about the convergence of the resulting algorithm are offered. Compared to robust gradient descent, the proposed algorithm enjoys the faster rates of convergence for successive iterates often achieved by second-order algorithms for convex problems, i.e., quadratic convergence in a neighborhood of the optimum, with a stepsize that may be chosen adaptively via backtracking linesearch.


title: Effect of crystal field engineering and Fermi level optimization on thermoelectric properties of Ge$_{1.01}$Te: Experimental investigation and theoretical insight

abstract: This study shows a method of enhancing the thermoelectric properties of GeTe-based materials by Ti and Bi co-doping on cation sites along with self-doping with Ge via simultaneous optimization of electronic (via crystal field engineering, and precise Fermi level optimization) and thermal (via point-defect scattering) transport properties. The pristine GeTe possesses high carrier concentration ($n$) due to intrinsic Ge vacancies, low Seebeck coefficient ($\alpha$), and high thermal conductivity ($\kappa$). The Ge vacancy optimization and crystal field engineering results in an enhanced $\alpha$ via excess Ge and Ti doping, which is further improved by band structure engineering through Bi doping. As a result of improved $\alpha$ and optimized Fermi level (carrier concentration), an enhanced power factor ($\alpha^2\sigma$) is obtained for Ti--Bi co-doped Ge$_{1.01}$Te. These experimental results are also evidenced by theoretical calculations of band structure, and thermoelectric parameters using density functional theory and Boltztrap calculations, respectively. A significant reduction in the phonon thermal conductivity ($\kappa_{ph}$) from $\sim$ 3.5 W.m$^{-1}$.K$^{-1}$ to $\sim$ 1.06 W.m$^{-1}$.K$^{-1}$ at 300\,K for Ti--Bi co-doping in GeTe, attributed to point-defect scattering due to mass and strain field fluctuation, in line with the Debye-Callaway model. The phonon dispersion calculations show a decreasing group velocity in Ti--Bi co-doped GeTe, supporting the obtained reduced $\kappa_{ph}$. The strategies used in the present study can significantly increase the effective mass, optimize the carrier concentration, and decrease phonon thermal conductivity while achieving an impressive maximum zT value of 1.75 at 773\,K and average zT (zT$_{av}$) of 1.03 for Ge$_{0.91}$Ti$_{0.02}$Bi$_{0.08}$Te over a temperature range of 300-773\,K.


title: Polynomial Preconditioning for Gradient Methods

abstract: We study first-order methods with preconditioning for solving structured nonlinear convex optimization problems. We propose a new family of preconditioners generated by symmetric polynomials. They provide first-order optimization methods with a provable improvement of the condition number, cutting the gaps between highest eigenvalues, without explicit knowledge of the actual spectrum. We give a stochastic interpretation of this preconditioning in terms of coordinate volume sampling and compare it with other classical approaches, including the Chebyshev polynomials. We show how to incorporate a polynomial preconditioning into the Gradient and Fast Gradient Methods and establish the corresponding global complexity bounds. Finally, we propose a simple adaptive search procedure that automatically chooses the best possible polynomial preconditioning for the Gradient Method, minimizing the objective along a low-dimensional Krylov subspace. Numerical experiments confirm the efficiency of our preconditioning strategies for solving various machine learning problems.


title: Adaptive Computation with Elastic Input Sequence

abstract: Humans have the ability to adapt the type of information they use, the procedure they employ, and the amount of time they spend when solving problems. However, most standard neural networks have a fixed function type and computation budget regardless of the sample's nature or difficulty. Adaptivity is a powerful paradigm as it not only imbues practitioners with flexibility pertaining to the downstream usage of these models but can also serve as a powerful inductive bias for solving certain challenging classes of problems. In this work, we introduce a new approach called AdaTape, which allows for dynamic computation in neural networks through adaptive tape tokens. AdaTape utilizes an elastic input sequence by equipping an architecture with a dynamic read-and-write tape. Specifically, we adaptively generate input sequences using tape tokens obtained from a tape bank which can be either trainable or derived from input data. We examine the challenges and requirements to obtain dynamic sequence content and length, and propose the Adaptive Tape Reading (ATR) algorithm to achieve both goals. Through extensive experiments on image recognition tasks, we show that AdaTape can achieve better performance while maintaining the computational cost. To facilitate further research, we have released code at https://github.com/google-research/scenic.


title: Looped Transformers as Programmable Computers

abstract: We present a framework for using transformer networks as universal computers by programming them with specific weights and placing them in a loop. Our input sequence acts as a punchcard, consisting of instructions and memory for data read/writes. We demonstrate that a constant number of encoder layers can emulate basic computing blocks, including embedding edit operations, non-linear functions, function calls, program counters, and conditional branches. Using these building blocks, we emulate a small instruction-set computer. This allows us to map iterative algorithms to programs that can be executed by a looped, 13-layer transformer. We show how this transformer, instructed by its input, can emulate a basic calculator, a basic linear algebra library, and in-context learning algorithms that employ backpropagation. Our work highlights the versatility of the attention mechanism, and demonstrates that even shallow transformers can execute full-fledged, general-purpose programs.


title: Unlocking Slot Attention by Changing Optimal Transport Costs

abstract: Slot attention is a powerful method for object-centric modeling in images and videos. However, its set-equivariance limits its ability to handle videos with a dynamic number of objects because it cannot break ties. To overcome this limitation, we first establish a connection between slot attention and optimal transport. Based on this new perspective we propose MESH (Minimize Entropy of Sinkhorn): a cross-attention module that combines the tiebreaking properties of unregularized optimal transport with the speed of regularized optimal transport. We evaluate slot attention using MESH on multiple object-centric learning benchmarks and find significant improvements over slot attention in every setting.


title: Transport across interfaces in symmetric orbifolds

abstract: We examine how conformal boundaries encode energy transport coefficients -- namely transmission and reflection probabilities -- of corresponding conformal interfaces in symmetric orbifold theories. These constitute a large class of irrational theories and are closely related to holographic setups. Our central goal is to compare such coefficients at the orbifold point (a field theory calculation) against their values when the orbifold is highly deformed (a gravity calculation) -- an approach akin to past AdS/CFT-guided comparisons of physical quantities at strong versus weak coupling. At the orbifold point, we find that the (weighted-average) transport coefficients are simply averages of coefficients in the underlying seed theory. We then focus on the symmetric orbifold of the $\mathbb{T}^4$ sigma model interface CFT dual to type IIB supergravity on the 3d Janus solution. We compare the holographic transmission coefficient, which was found by [1], to that of the orbifold point. We find that the profile of the transmission coefficient substantially increases with the coupling, in contrast to boundary entropy. We also present some related ideas about twisted-sector data encoded by boundary states.


title: Streaming Anomaly Detection

abstract: Anomaly detection is critical for finding suspicious behavior in innumerable systems. We need to detect anomalies in real-time, i.e. determine if an incoming entity is anomalous or not, as soon as we receive it, to minimize the effects of malicious activities and start recovery as soon as possible. Therefore, online algorithms that can detect anomalies in a streaming manner are essential. We first propose MIDAS which uses a count-min sketch to detect anomalous edges in dynamic graphs in an online manner, using constant time and memory. We then propose two variants, MIDAS-R which incorporates temporal and spatial relations, and MIDAS-F which aims to filter away anomalous edges to prevent them from negatively affecting the internal data structures. We then extend the count-min sketch to a Higher-Order sketch to capture complex relations in graph data, and to reduce detecting suspicious dense subgraph problem to finding a dense submatrix in constant time. Using this sketch, we propose four streaming methods to detect edge and subgraph anomalies. Next, we broaden the graph setting to multi-aspect data. We propose MStream which detects explainable anomalies in multi-aspect data streams. We further propose MStream-PCA, MStream-IB, and MStream-AE to incorporate correlation between features. Finally, we consider multi-dimensional data streams with concept drift and propose MemStream. MemStream leverages the power of a denoising autoencoder to learn representations and a memory module to learn the dynamically changing trend in data without the need for labels. We prove a theoretical bound on the size of memory for effective drift handling. In addition, we allow quick retraining when the arriving stream becomes sufficiently different from the training data. Furthermore, MemStream makes use of two architecture design choices to be robust to memory poisoning.


title: Liouville conformal field theory and the quantum zipper

abstract: Sheffield showed that conformally welding a $\gamma$-Liouville quantum gravity (LQG) surface to itself gives a Schramm-Loewner evolution (SLE) curve with parameter $\kappa = \gamma^2$ as the interface, and Duplantier-Miller-Sheffield proved similar stories for $\kappa = {16}/{\gamma^2}$ for $\gamma$-LQG surfaces with boundaries decorated by looptrees of disks or by continuum random trees. We study these dynamics for LQG surfaces coming from Liouville conformal field theory (LCFT). At stopping times depending only on the curve, we give an explicit description of the surface and curve in terms of LCFT and SLE. This has applications to both LCFT and SLE. We prove the boundary BPZ equation for LCFT, which is crucial to solving boundary LCFT. With Yu we will prove the reversibility of whole-plane $\mathrm{SLE}_\kappa$ for $\kappa > 8$ via a novel radial mating-of-trees, and show the space of LCFT surfaces is closed under conformal welding.


title: Numerical Issues for a Non-autonomous Logistic Model

abstract: The logistic equation has been extensively used to model biological phenomena across a variety of disciplines and has provided valuable insight into how our universe operates. Incorporating time-dependent parameters into the logistic equation allows the modeling of more complex behavior than its autonomous analog, such as a tumor's varying growth rate under treatment, or the expansion of bacterial colonies under varying resource conditions. Some of the most commonly used numerical solvers produce vastly different approximations for a non-autonomous logistic model with a periodically-varying growth rate changing signum. Incorrect, inconsistent, or even unstable approximate solutions for this non-autonomous problem can occur from some of the most frequently used numerical methods, including the lsoda, implicit backwards difference, and Runge-Kutta methods, all of which employ a black-box framework. Meanwhile, a simple, manually-programmed Runge-Kutta method is robust enough to accurately capture the analytical solution for biologically reasonable parameters and consistently produce reliable simulations. Consistency and reliability of numerical methods are fundamental for simulating non-autonomous differential equations and dynamical systems, particularly when applications are physically or biologically informed.


title: Multipoint Turbulence Analysis with Helioswarm

abstract: Exploration of plasma dynamics in space, including turbulence, is entering a new era of multi-satellite constellation measurements that will determine fundamental properties with unprecedented precision. Familiar but imprecise approximations will need to be abandoned and replaced with more advanced approaches. We present a preparatory study of the evaluation of second- and third-order statistics, using simultaneous measurements at many points. Here, for specificity, the orbital configuration of the NASA Helioswarm mission is employed in conjunction with three-dimensional magnetohydrodynamics numerical simulations of turbulence. The Helioswarm 9-spacecraft constellation flies virtually through the turbulence to compare results with the exact numerical statistics. We demonstrate novel increment-based techniques for the computation of (1) the multidimensional spectra and (2) the turbulent energy flux. This latter increment-space estimate of the cascade rate, based on the third-order Yaglom-Politano-Pouquet theory, uses numerous increment-space tetrahedra. Our investigation reveals that Helioswarm will provide crucial information on the nature of astrophysical turbulence.


title: The moment map for the variety of Leibniz algebras

abstract: We consider the moment map $m:\mathbb{P}V_n\rightarrow \text{i}\mathfrak{u}(n)$ for the action of $\text{GL}(n)$ on $V_n=\otimes^{2}(\mathbb{C}^{n})^{*}\otimes\mathbb{C}^{n}$, and study the functional $F_n=\|m\|^{2}$ restricted to the projectivizations of the algebraic varieties of all $n$-dimensional Leibniz algebras $L_n$ and all $n$-dimensional symmetric Leibniz algebras $S_n$, respectively. Firstly, we give a description of the maxima and minima of the functional $F_n: L_n \rightarrow \mathbb{R}$, proving that they are actually attained at the symmetric Leibniz algebras. Then, for an arbitrary critical point $[\mu]$ of $F_n: S_n \rightarrow \mathbb{R}$, we characterize the structure of $[\mu]$ by virtue of the nonnegative rationality. Finally, we classify the critical points of $F_n: S_n \rightarrow \mathbb{R}$ for $n=2$, $3$, respectively.


title: Forex Trading Strategy That Might Be Executed Due to the Popularity of Gotobi Anomaly

abstract: Our previous research has confirmed that the USD/JPY rate tends to rise toward 9:55 every morning in the Gotobi days, which are divisible by five. This is called the Gotobi anomaly. In the present study, we verify the possible trading strategy and its validity under the condition that investors recognize the existence of the anomaly. Moreover, we illustrate the possibility that the wealth of Japanese companies might leak to FX traders due to the arbitrage opportunity if Japanese companies blindly keep making payments in the Gotobi days as a business custom.


title: Representations and identities of Baxter monoids with involution

abstract: Let $(\mathsf{baxt}_n,~^\sharp)$ be the Baxter monoid of finite rank $n$ with Sch\"{u}tzenberger's involution $^{\sharp}$. In this paper, it is shown that $(\mathsf{baxt}_n,~^\sharp)$ admits a faithful representation by an involution monoid of upper triangular matrices over any semiring from a large class including the tropical semiring under the skew transposition. Then a transparent combinatorial characterization of the word identities satisfied by $(\mathsf{baxt}_n,~^\sharp)$ is given. Further, it is proved that $(\mathsf{baxt}_n,~^\sharp)$ is finitely based if and only if $n\neq 3$, and shown that the identity checking problem for $(\mathsf{baxt}_n,~^\sharp)$ can be done in polynomial time.


title: Purity and quasi-split torsors over Pr\"ufer bases

abstract: We establish an analogue of the Zariski--Nagata purity theorem for finite \'etale covers on smooth schemes over Pr\"ufer rings by demonstrating Auslander's flatness criterion in this non-Noetherian context. We derive an Auslander--Buchsbaum formula for general local rings, which provides a useful tool for studying the algebraic structures involved in our work. Through analysis of reflexive sheaves, we prove various purity theorems for torsors under certain group algebraic spaces, such as the reductive ones. Specifically, using results from EGAIV4 on parafactoriality on smooth schemes over normal bases, we prove the purity for cohomology groups of multiplicative type groups at this level of generality. Subsequently, we leverage the aforementioned purity results to resolve the Grothendieck--Serre conjecture for torsors under a quasi-split reductive group scheme over schemes smooth over Pr\"ufer rings. Along the way, we also prove a version of the Nisnevich purity conjecture for quasi-split reductive group schemes in our Pr\"uferian context, inspired by the recent work of Cesnavicius.


title: A quantum trajectory analysis of singular wave functions

abstract: The Schr\"{o}dinger equation admits smooth and finite solutions that spontaneously evolve into a singularity, even for a free particle. This blowup is generally ascribed to the intrinsic dispersive character of the associated time evolution. We resort to the notion of quantum trajectories to reinterpret this singular behavior. We show that the blowup can be directly related to local phase variations, which generate an underlying velocity field responsible for driving the quantum flux toward the singular region.


title: The impact of magnetic fields on cosmological galaxy mergers -- II. Modified angular momentum transport and feedback

abstract: The role of magnetic fields in galaxy evolution is still an unsolved question in astrophysics. We have previously shown that magnetic fields play a crucial role in major mergers between disc galaxies; in hydrodynamic simulations of such mergers, the Auriga model produces compact remnants with a distinctive bar and ring morphology. In contrast, in magnetohydrodynamic (MHD) simulations, remnants form radially-extended discs with prominent spiral arm structure. In this paper, we analyse a series of cosmological "zoom-in" simulations of major mergers and identify exactly \textit{how} magnetic fields are able to alter the outcome of the merger. We find that magnetic fields modify the transport of angular momentum, systematically hastening the merger progress. The impact of this altered transport depends on the orientation of the field, with a predominantly non-azimuthal (azimuthal) orientation increasing the central baryonic concentration (providing support against collapse). Both effects act to suppress an otherwise existent bar-instability, which in turn leads to a fundamentally different morphology and manifestation of feedback. We note, in particular, that stellar feedback is substantially less influential in MHD simulations, which allows for the later accretion of higher angular momentum gas and the subsequent rapid radial growth of the remnant disc. A corollary of the increased baryonic concentration in MHD simulations is that black holes are able to grow twice as large, although this turns out to have little impact on the remnant's development. Our results show that galaxy evolution cannot be modelled correctly without including magnetic fields.


title: Finite key performance of satellite quantum key distribution under practical constraints

abstract: Global-scale quantum communication networks will require efficient long-distance distribution of quantum signals. Optical fibre communication channels have range constraints due to exponential losses in the absence of quantum memories and repeaters. Satellites enable intercontinental quantum communication by exploiting more benign inverse square free-space attenuation and long sight lines. However, the design and engineering of satellite quantum key distribution (QKD) systems are difficult and characteristic differences to terrestrial QKD networks and operations pose additional challenges. The typical approach to modelling satellite QKD (SatQKD) has been to estimate performances with a fully optimised protocol parameter space and with few payload and platform resource limitations. Here, we analyse how practical constraints affect the performance of SatQKD for the Bennett-Brassard 1984 (BB84) weak coherent pulse decoy state protocol with finite key size effects. We consider engineering limitations and trade-offs in mission design including limited in-orbit tunability, quantum random number generation rates and storage, and source intensity uncertainty. We quantify practical SatQKD performance limits to determine the long-term key generation capacity and provide important performance benchmarks to support the design of upcoming missions.


title: Observables for moving, stupendously charged and massive primordial black holes

abstract: Stupendously large black holes exceeding $10^{11} M_\odot$ could exist, supported by recent observations of unexpectedly massive black holes at high redshifts. These objects may constitute a part of dark matter or even dark energy. One possibility to explain the cosmic accelerated expansion could be to consider charged black holes whose mutual repulsion overcomes their gravitational attraction. However, the extreme charge required turns these black holes into naked singularities, whose existence is questioned by the cosmic censorship hypothesis. Since the latter is driven by theoretical assumptions, we work out the most promising observables which are least cosmology-dependent to test their existence. We derive the electro-magnetic and gravitational lensing effects caused by such extreme objects at distances much larger than their extent to investigate possible ways for a discovery. Restricting searches to black holes between $10^{12}$ to $10^{14} M_\odot$, we show that such objects do not cause totally disruptive catastrophes, like dissociation of neutral hydrogen clouds or proton decay induced by strong electro-magnetic fields. Einstein rings of the order of 10" and rotation measures of plasma clouds subject to the magnetic fields induced by the moving black holes are identified as optimum observable signatures for now. Future space-based black-hole telescopes will follow up on these candidates and finally check the cosmic censorship hypothesis by their strong-field strong-lensing signatures, like an additional sub-arcsecond inner Einstein ring. Observable effects are so surprisingly moderate that a violation of cosmic censorship is hard to detect and even explaining cosmic expansion with moving naked singularities might be possible.


title: Earth mover's distance as a measure of CP violation

abstract: We introduce a new unbinned two sample test statistic sensitive to CP violation utilizing the optimal transport plan associated with the Wasserstein (earth mover's) distance. The efficacy of the test statistic is shown via two examples of CP asymmetric distributions with varying sample sizes: the Dalitz distributions of $B^0 \rightarrow K^+\pi^-\pi^0$ and of $D^0 \rightarrow \pi^+\pi^-\pi^0$ decays. The windowed version of the Wasserstein distance test statistic is shown to have comparable sensitivity to CP violation as the commonly used energy test statistic, but also retains information about the localized distributions of CP asymmetry over the Dalitz plot. For large statistic datasets we introduce two modified Wasserstein distance based test statistics -- the binned and the sliced Wasserstein distance statistics, which show comparable sensitivity to CP violation, but improved computing time and memory scalings. Finally, general extensions and applications of the introduced statistics are discussed.


title: Entanglement is better teleported than transmitted

abstract: We show that, for the purpose of quantum communication via a quantum field, it is essential to view the field not only as a medium for transmission but also as a source of entanglement that can aid in the communication task. To this end, we consider the quantum communication scenario where Alice is initially entangled with an ancilla and intends to communicate with Bob through a quantum field, so as to make Bob entangled with the ancilla. We find that if Alice and Bob communicate by directly coupling to the quantum field, then they can generate negativity between Bob and the ancilla only at orders that are higher than second perturbative order. We then present a protocol based on quantum teleportation in which Alice and Bob consume entanglement that they obtained from the field via interaction or harvesting. We show that this protocol can transfer negativity already to second perturbative order.