title: A comparison of particle and fluid models for positive streamer discharges in air

abstract: Both fluid and particle models are commonly used to simulate streamer discharges. In this paper, we quantitatively study the agreement between these approaches for axisymmetric and 3D simulations of positive streamers in air. We use a drift-diffusion-reaction fluid model with the local field approximation and a PIC-MCC (particle-in-cell, Monte Carlo collision) particle model. The simulations are performed at 300 K and 1 bar in a 10 mm plate-plate gap with a 2 mm needle electrode. Applied voltages between 11.7 and 15.6 kV are used, which correspond to background fields of about 15 to 20 kV/cm. Streamer properties like maximal electric field, head position and velocity are compared as a function of time or space. Our results show good agreement between the particle and fluid simulations, in contrast to some earlier comparisons that were carried out in 1D or for negative streamers. To quantify discrepancies between the models, we mainly look at streamer velocities as a function of streamer length. For the test cases considered here, the mean deviation in streamer velocity between the particle and fluid simulations is less than 4%. We study the effect of different types of transport data for the fluid model, and find that flux coefficients lead to good agreement whereas bulk coefficients do not. Furthermore, we find that with a two-term Boltzmann solver, data should be computed using a temporal growth model for the best agreement. The numerical convergence of the particle and fluid models is also studied. In fluid simulations the streamer velocity increases somewhat using finer grids, whereas the particle simulations are less sensitive to the grid. Photoionization is the dominant source of stochastic fluctuations in our simulations. When the same stochastic photoionization model is used, particle and fluid simulations exhibit similar fluctuations.


title: Convergence of Restricted Additive Schwarz with impedance transmission conditions for discretised Helmholtz problems

abstract: The Restricted Additive Schwarz method with impedance transmission conditions, also known as the Optimised Restricted Additive Schwarz (ORAS) method, is a simple overlapping one-level parallel domain decomposition method, which has been successfully used as an iterative solver and as a preconditioner for discretized Helmholtz boundary-value problems. In this paper, we give, for the first time, a convergence analysis for ORAS as an iterative solver -- and also as a preconditioner -- for nodal finite element Helmholtz systems of any polynomial order. The analysis starts by showing (for general domain decompositions) that ORAS as an unconventional finite element approximation of a classical parallel iterative Schwarz method, formulated at the PDE (non-discrete) level. This non-discrete Schwarz method was recently analysed in [Gong, Gander, Graham, Lafontaine, Spence, arXiv 2106.05218], and the present paper gives a corresponding discrete version of this analysis. In particular, for domain decompositions in strips in 2-d, we show that, when the mesh size is small enough, ORAS inherits the convergence properties of the Schwarz method, independent of polynomial order. The proof relies on characterising the ORAS iteration in terms of discrete `impedance-to-impedance maps', which we prove (via a novel weighted finite-element error analysis) converge as $h\rightarrow 0$ in the operator norm to their non-discrete counterparts.


title: Static, spherically symmetric objects in Type-II minimally modified gravity

abstract: Static, spherically symmetric solutions representing stars made of barotropic perfect fluid are studied in the context of two theories of type-II minimally modified gravity, VCDM and VCCDM. Both of these theories share the property that no additional degree of freedom is introduced in the gravity sector, and propagate only two gravitational waves besides matter fields, as in General Relativity (GR). We find that, on imposing physical boundary conditions on the Misner-Sharp mass of the system, the solutions in V(C)CDM exactly coincide with the ones in GR, namely they also satisfy the Tolman-Oppenheimer-Volkoff equation.


title: Protein induced lipid demixing in homogeneous membranes

abstract: Specific lipid environments are necessary for the establishment of protein signalling platforms in membranes, yet their origin has been highly debated. We present a continuum, exactly solvable model of protein induced local demixing of lipid membranes. The coupling between a local composition and a local thickness of the membrane induces lipid domains around inclusions with hydrophobic mismatch, even for temperatures above the miscibility critical point of the membrane. The model qualitatively explains the experimentally observed formation of lipid domains induced by anchoring of reconstituted actin in flat supported lipid bilayers.


title: Structural Parameterizations of Budgeted Graph Coloring

abstract: We introduce a variant of the graph coloring problem, which we denote as {\sc Budgeted Coloring Problem} (\bcp). Given a graph $G$, an integer $c$ and an ordered list of integers $\{b_1, b_2, \ldots, b_c\}$, \bcp asks whether there exists a proper coloring of $G$ where the $i$-th color is used to color at most $b_i$ many vertices. This problem generalizes two well-studied graph coloring problems, {\sc Bounded Coloring Problem} (\bocp) and {\sc Equitable Coloring Problem} (\ecp) and as in the case of other coloring problems, it is \nph even for constant values of $c$. So we study \bcp under the paradigm of parameterized complexity. \begin{itemize} \item We show that \bcp is \fpt (fixed-parameter tractable) parameterized by the vertex cover size. This generalizes a similar result for \ecp and immediately extends to the \bocp, which was earlier not known. \item We show that \bcp is polynomial time solvable for cluster graphs generalizing a similar result for \ecp. However, we show that \bcp is \fpt, but unlikely to have polynomial kernel, when parameterized by the deletion distance to clique, contrasting the linear kernel for \ecp for the same parameter. \item While the \bocp is known to be polynomial time solvable on split graphs, we show that \bcp is \nph on split graphs. As \bocp is hard on bipartite graphs when $c>3$, the result follows for \bcp as well. We provide a dichotomy result by showing that \bcp is polynomial time solvable on bipartite graphs when $c=2$. We also show that \bcp is \nph on co-cluster graphs, contrasting the polynomial time algorithm for \ecp and \bocp. \end{itemize} Finally we present an $\mathcal{O}^*(2^{|V(G)|})$ algorithm for the \bcp, generalizing the known algorithm with a similar bound for the standard chromatic number.


title: Robustness of ARS Leptogenesis in Scalar Extensions

abstract: Extensions of the Standard Model (SM) with sterile neutrinos are well motivated from the observed oscillations of the light neutrinos and they have shown to successfully explain the Baryon Asymmetry of the Universe (BAU) through, for instance, the so-called ARS leptogenesis. Sterile neutrinos can be added in minimal ways to the SM, but many theories exist where sterile neutrinos are not the only new fields. Such theories often include scalar bosons, which brings about the possibility of further interactions between the sterile neutrinos and the SM. In this paper we consider an extension of the SM with two sterile neutrinos and one scalar singlet particle and investigate the effect that an additional, thermalised, scalar has on the ARS leptogenesis mechanism. We show that in general the created asymmetry is reduced due to additional sterile neutrino production from scalar decays. When sterile neutrinos and scalars are discovered in the laboratory, our results will provide information on the applicability of the ARS leptogenesis mechanism.


title: Production of Chern-Simons bosons in decays of mesons

abstract: We consider the effective interaction of quarks with a new GeV-scale vector particle that couples to electroweak gauge bosons by the so-called effective Chern-Simons interaction. We call this particle the Chern-Simons (CS) boson. We construct an effective Lagrangian of the CS boson interaction with quarks of two different flavors. This interaction is given by a divergent loop diagram, however, it turns out that the divergent part is equal to zero as a consequence of the CKM matrix unitarity in the SM. Therefore, we are able to predict the effective interaction of the CS boson with quarks of different flavors without introducing new unknown parameters to the model, using only parameters of the initial effective Lagrangian. Our result shows that the effective interaction of the CS boson with down-type quarks is sufficiently stronger compared with up-type quarks. Based on our results, we give a prediction for the production of CS bosons in mesons decays. Branching fractions were obtained for the main reactions of the CS production in meson decays. The results obtained will be useful for searching for the long-lived GeV-scale CS boson in intensity frontier experiments.


title: Size-dependent enhancement of passive microwave rectification in magnetic tunnel junctions with perpendicular magnetic anisotropy

abstract: Spintronic rf detectors were demonstrated, recently, for energy harvesting and wireless communication at low input power. Here we report on the optimization of the rectified output dc voltage using magnetic tunnel junctions (MTJ) with strong perpendicular anisotropy (PMA) of both the polarizing and the free layer. The magnetization of the polarizing layer is fixed out of plane, while the free layer thickness is adjusted so that its magnetization orientation changes from in plane to out of plane. The rectification dc output voltage lies in the mV range for moderate rf powers, with a signal to noise ratio of 10 to 100 for Prf = -25dBm. It shows a strong dependence on the dimensions of the MTJ: it increases by a factor of 5 to 6 when reducing the diameter from 150nm to 20nm. This enhancement can be doubled when reducing the FL thickness from 1.8nm to 1.6nm. This dimensional enhancement is attributed to the change of the effective anisotropy of the excited free layer, and the MTJ resistance. The results are of interest for the design of spintronic based rf detectors with optimized sensitivity.


title: Closing the "Quantum Supremacy" Gap: Achieving Real-Time Simulation of a Random Quantum Circuit Using a New Sunway Supercomputer

abstract: We develop a high-performance tensor-based simulator for random quantum circuits(RQCs) on the new Sunway supercomputer. Our major innovations include: (1) a near-optimal slicing scheme, and a path-optimization strategy that considers both complexity and compute density; (2) a three-level parallelization scheme that scales to about 42 million cores; (3) a fused permutation and multiplication design that improves the compute efficiency for a wide range of tensor contraction scenarios; and (4) a mixed-precision scheme to further improve the performance. Our simulator effectively expands the scope of simulatable RQCs to include the 10*10(qubits)*(1+40+1)(depth) circuit, with a sustained performance of 1.2 Eflops (single-precision), or 4.4 Eflops (mixed-precision)as a new milestone for classical simulation of quantum circuits; and reduces the simulation sampling time of Google Sycamore to 304 seconds, from the previously claimed 10,000 years.


title: Simple data balancing achieves competitive worst-group-accuracy

abstract: We study the problem of learning classifiers that perform well across (known or unknown) groups of data. After observing that common worst-group-accuracy datasets suffer from substantial imbalances, we set out to compare state-of-the-art methods to simple balancing of classes and groups by either subsampling or reweighting data. Our results show that these data balancing baselines achieve state-of-the-art-accuracy, while being faster to train and requiring no additional hyper-parameters. In addition, we highlight that access to group information is most critical for model selection purposes, and not so much during training. All in all, our findings beg closer examination of benchmarks and methods for research in worst-group-accuracy optimization.


title: An introduction to Goldstone boson physics and to the coset construction

abstract: These lecture notes are based on a six-hour series of lectures given at the XVII Modave summer school in mathematical physics, aimed at Ph.D. students in high-energy theoretical physics. The manuscript starts by briefly stating Goldstone's theorem and emphasises the motivations behind Goldstone physics; the main asset being the universality of spontaneous symmetry breaking (SSB) which is the fundamental hypothesis of Goldstone's theorem. Once the different notions of SSB will be clarified/reviewed, Goldstone's theorem will be stated and proved. A prediction of this theorem is the existence of gapless particles, called Nambu-Goldstone modes (NG modes). From the discussion on Goldstone's results, some aspects of the NG modes will emerge. Besides to be gapless, they are systematically weakly coupled at low energy. Therefore, an effective field theory (EFT) building tool called ``coset construction'' will be presented to explicitly display these specific features of the NG modes. The coset construction suits our goal since it is mainly based on the symmetry realisation of the perturbed theory around the background inducing SSB. From the general obtained EFT, a counting rule for the NG modes will be derived. The limitations of this rule as well as the still ongoing generalisation will be discussed (e.g. spacetime symmetry breaking). The tools developed during this course will be illustrated with a concrete example in physics: ferromagnetism. The notes end with a brief state of the art of Goldstone physics. This provides some directions into which the interested reader could investigate to expand his knowledge on the subject. N.B.: No prerequisites are required beside the standard courses of a Master in theoretical physics.

categories: hep-th cond-mat.other hep-ph math-ph math.MP


title: Dimer-dimer correlations at the rough-smooth boundary

abstract: Three phases of macroscopic domains have been seen for large but finite periodic dimer models; these are known as the frozen, rough and smooth phases. The transition region between the frozen and rough region has received a lot of attention for the last twenty years and recently work has been underway to understand the rough-smooth transition region in the case of the two-periodic Aztec diamond. We compute uniform asymptotics for dimer-dimer correlations of the two-periodic Aztec diamond when the dimers lie in the rough-smooth transition region. These asymptotics rely on a formula found in [5] for the inverse Kasteleyn matrix, they also apply to a related infinite dimer model.


title: Frequency multiplexed entanglement for continuous-variable quantum key distribution

abstract: Quantum key distribution with continuous variables already uses advantageous high-speed single-mode homodyne detection with low electronic noise at room temperature. Together with continuous-variable information encoding to nonclassical states, the distance for secure key transmission through lossy channels can approach 300 km in current optical fibers. Such protocols tolerate higher channel noise and also limited data processing efficiency compared to coherent-state protocols. The secret key rate can be further increased by increasing the system clock rates, and, further, by a suitable frequency-mode-multiplexing of optical transmission channels. However, the multiplexed modes couple together in the source or any other part of the protocol. Therefore, multiplexed communication will experience crosstalk and the gain can be minuscule. Advantageously, homodyne detectors allow solving this crosstalk problem by proper data processing. It is a potential advantage over protocols with single-photon detectors, which do not enable similar data processing techniques. We demonstrate the positive outcome of this methodology on the experimentally characterized frequency-multiplexed entangled source of femtosecond optical pulses with natural crosstalk between eight entangled pairs of modes. As the main result, we predict almost 15-fold higher secret key rate. This experimental test and analysis of frequency-multiplexed entanglement source opens the way for the field implementation of high-capacity quantum key distribution with continuous variables.


title: Detecting the ultra low dimensionality of real networks

abstract: Reducing dimension redundancy to find simplifying patterns in high-dimensional datasets and complex networks has become a major endeavor in many scientific fields. However, detecting the dimensionality of their latent space is challenging but necessary to generate efficient embeddings to be used in a multitude of downstream tasks. Here, we propose a method to infer the dimensionality of networks without the need for any a priori spatial embedding. Due to the ability of hyperbolic geometry to capture the complex connectivity of real networks, we detect ultra low dimensionality far below values reported using other approaches. We applied our method to real networks from different domains and found unexpected regularities, including: tissue-specific biomolecular networks being extremely low dimensional; brain connectomes being close to the three dimensions of their anatomical embedding; and social networks and the Internet requiring slightly higher dimensionality. Beyond paving the way towards an ultra efficient dimensional reduction, our findings help address fundamental issues that hinge on dimensionality, such as universality in critical behavior.


title: Finding Regions of Heterogeneity in Decision-Making via Expected Conditional Covariance

abstract: Individuals often make different decisions when faced with the same context, due to personal preferences and background. For instance, judges may vary in their leniency towards certain drug-related offenses, and doctors may vary in their preference for how to start treatment for certain types of patients. With these examples in mind, we present an algorithm for identifying types of contexts (e.g., types of cases or patients) with high inter-decision-maker disagreement. We formalize this as a causal inference problem, seeking a region where the assignment of decision-maker has a large causal effect on the decision. Our algorithm finds such a region by maximizing an empirical objective, and we give a generalization bound for its performance. In a semi-synthetic experiment, we show that our algorithm recovers the correct region of heterogeneity accurately compared to baselines. Finally, we apply our algorithm to real-world healthcare datasets, recovering variation that aligns with existing clinical knowledge.


title: Deep Transfer Learning for Multi-source Entity Linkage via Domain Adaptation

abstract: Multi-source entity linkage focuses on integrating knowledge from multiple sources by linking the records that represent the same real world entity. This is critical in high-impact applications such as data cleaning and user stitching. The state-of-the-art entity linkage pipelines mainly depend on supervised learning that requires abundant amounts of training data. However, collecting well-labeled training data becomes expensive when the data from many sources arrives incrementally over time. Moreover, the trained models can easily overfit to specific data sources, and thus fail to generalize to new sources due to significant differences in data and label distributions. To address these challenges, we present AdaMEL, a deep transfer learning framework that learns generic high-level knowledge to perform multi-source entity linkage. AdaMEL models the attribute importance that is used to match entities through an attribute-level self-attention mechanism, and leverages the massive unlabeled data from new data sources through domain adaptation to make it generic and data-source agnostic. In addition, AdaMEL is capable of incorporating an additional set of labeled data to more accurately integrate data sources with different attribute importance. Extensive experiments show that our framework achieves state-of-the-art results with 8.21% improvement on average over methods based on supervised learning. Besides, it is more stable in handling different sets of data sources in less runtime.


title: Proton deflectometry with in situ x-ray reference for absolute measurement of electromagnetic fields in high-energy-density plasmas

abstract: We report a technique of proton deflectometry which uses a grid and an in situ reference x-ray grid image for precise measurements of magnetic fields in high-energy density plasmas. A D$^3$He fusion implosion provides a bright point-source of both protons and x-rays, which is split into beamlets by a mesh grid. The protons undergo deflections as they propagate through the plasma region of interest, whereas the x-rays travel along straight lines. The x-ray image therefore provides a zero-deflection reference image. The line-integrated magnetic fields are inferred from the shifts of beamlets between the deflected (proton) and reference (x-ray) images. We developed a system for analysis of this data, including automatic algorithms to find beamlet locations and calculate their deflections from the reference image. The technique is verified in an experiment performed at OMEGA to measure a non-uniform magnetic field in vacuum, then applied to observe the interaction of an expanding plasma plume with the magnetic field.


title: Standard meta-analysis methods are not robust

abstract: P values or risk ratios from multiple, independent studies, observational or randomized, can be computationally combined to provide an overall assessment of a research question in meta-analysis. There is a need to examine the reliability of these methods of combination. It is typical in observational studies to statistically test many questions and not correct the analysis results for multiple testing or multiple modeling, MTMM. The same problem can happen for randomized, experimental trials. There is the additional problem that some of the base studies may be using fabricated or fraudulent data. If there is no attention to MTMM or fraud in the base studies, there is no guarantee that the results to be combined are unbiased, the key requirement for the valid combining of results. We note that methods of combination are not robust; even one extreme base study value can overwhelm standard methods of combination. It is possible that multiple, extreme (MTMM or fraudulent) results can feed from the base studies to bias the combined result. A meta-analysis of observational (or even randomized studies) may not be reliable. Examples are given along with some methods to evaluate existing base studies and meta-analysis studies.


title: Learning algebraic structures with the help of Borel equivalence relations

abstract: We study algorithmic learning of algebraic structures. In our framework, a learner receives larger and larger pieces of an arbitrary copy of a computable structure and, at each stage, is required to output a conjecture about the isomorphism type of such a structure. The learning is successful if the conjectures eventually stabilize to a correct guess. We prove that a family of structures is learnable if and only if its learning domain is continuously reducible to the relation $E_0$ of eventual agreement on reals. This motivates a novel research program, that is, using descriptive set theoretic tools to calibrate the (learning) complexity of nonlearnable families. Here, we focus on the learning power of well-known benchmark Borel equivalence relations (i.e., $E_1$, $E_2$, $E_3$, $Z_0$, and $E_{set}$).


title: Neural Analysis and Synthesis: Reconstructing Speech from Self-Supervised Representations

abstract: We present a neural analysis and synthesis (NANSY) framework that can manipulate voice, pitch, and speed of an arbitrary speech signal. Most of the previous works have focused on using information bottleneck to disentangle analysis features for controllable synthesis, which usually results in poor reconstruction quality. We address this issue by proposing a novel training strategy based on information perturbation. The idea is to perturb information in the original input signal (e.g., formant, pitch, and frequency response), thereby letting synthesis networks selectively take essential attributes to reconstruct the input signal. Because NANSY does not need any bottleneck structures, it enjoys both high reconstruction quality and controllability. Furthermore, NANSY does not require any labels associated with speech data such as text and speaker information, but rather uses a new set of analysis features, i.e., wav2vec feature and newly proposed pitch feature, Yingram, which allows for fully self-supervised training. Taking advantage of fully self-supervised training, NANSY can be easily extended to a multilingual setting by simply training it with a multilingual dataset. The experiments show that NANSY can achieve significant improvement in performance in several applications such as zero-shot voice conversion, pitch shift, and time-scale modification.


title: Streaming Generalized Canonical Polyadic Tensor Decompositions

abstract: In this paper, we develop a method which we call OnlineGCP for computing the Generalized Canonical Polyadic (GCP) tensor decomposition of streaming data. GCP differs from traditional canonical polyadic (CP) tensor decompositions as it allows for arbitrary objective functions which the CP model attempts to minimize. This approach can provide better fits and more interpretable models when the observed tensor data is strongly non-Gaussian. In the streaming case, tensor data is gradually observed over time and the algorithm must incrementally update a GCP factorization with limited access to prior data. In this work, we extend the GCP formalism to the streaming context by deriving a GCP optimization problem to be solved as new tensor data is observed, formulate a tunable history term to balance reconstruction of recently observed data with data observed in the past, develop a scalable solution strategy based on segregated solves using stochastic gradient descent methods, describe a software implementation that provides performance and portability to contemporary CPU and GPU architectures and integrates with Matlab for enhanced useability, and demonstrate the utility and performance of the approach and software on several synthetic and real tensor data sets.


title: Automatic extraction of one-loop Wilson coefficients in general BSM scenarios using MARTY-1.4

abstract: We present a fully automated procedure providing an easy way to perform, systematically, phenomenological analyses in flavor physics for general BSM scenarios. This procedure relies on MARTY-1.4, is model independent and requires as input only the Lagrangian of the theory. Once the Lagrangian has been defined, tree-level and one-loop Wilson coefficients can be calculated symbolically by MARTY, from which flavor observables can be computed numerically by available software programs. We focus in particular on $b\to s\gamma$ and the recently updated $b\to s\ell^+\ell^-$ observables which are in tension with the SM, and present a general procedure to extract the relevant one-loop coefficients, such as $C_7$, $C_7^\prime$, $C_9$ and $C_{10}$.


title: Self-Contained Kinematic Calibration of a Novel Whole-Body Artificial Skin for Human-Robot Collaboration

abstract: In this paper, we present an accelerometer-based kinematic calibration algorithm to accurately estimate the pose of multiple sensor units distributed along a robot body. Our approach is self-contained, can be used on any robot provided with a Denavit-Hartenberg kinematic model, and on any skin equipped with Inertial Measurement Units (IMUs). To validate the proposed method, we first conduct extensive experimentation in simulation and demonstrate a sub-cm positional error from ground truth data --an improvement of six times with respect to prior work; subsequently, we then perform a real-world evaluation on a seven degrees-of-freedom collaborative platform. For this purpose, we additionally introduce a novel design for a stand-alone artificial skin equipped with an IMU for use with the proposed algorithm and a proximity sensor for sensing distance to nearby objects. In conclusion, in this work, we demonstrate seamless integration between a novel hardware design, an accurate calibration method, and preliminary work on applications: the high positional accuracy effectively enables to locate distributed proximity data and allows for a distributed avoidance controller to safely avoid obstacles and people without the need of additional sensing.


title: Geometric phase in a dissipative Jaynes-Cummings model: theoretical explanation for resonance robustness

abstract: We follow a generalized kinematic approach to compute the geometric phases acquired in both unitary and dissipative Jaynes-Cummings models, which provide a fully quantum description for a two-level system interacting with a single mode of the (cavity) electromagnetic field, in a perfect or dissipative cavity respectively. In the dissipative model, the non-unitary effects arise from the outflow of photons through the cavity walls and the incoherent pumping of the two-level system. Our approach allows to compare the geometric phases acquired in these models, leading to an exhaustive characterization of the corrections introduced by the presence of the environment. We also provide geometric interpretations for the observed behaviors. When the resonance condition is satisfied, we show the geometric phase is robust, exhibiting a vanishing correction under a non-unitary evolution. This fact is supported with a geometrical explanation as well.


title: NIDA-CLIFGAN: Natural Infrastructure Damage Assessment through Efficient Classification Combining Contrastive Learning, Information Fusion and Generative Adversarial Networks

abstract: During natural disasters, aircraft and satellites are used to survey the impacted regions. Usually human experts are needed to manually label the degrees of the building damage so that proper humanitarian assistance and disaster response (HADR) can be achieved, which is labor-intensive and time-consuming. Expecting human labeling of major disasters over a wide area gravely slows down the HADR efforts. It is thus of crucial interest to take advantage of the cutting-edge Artificial Intelligence and Machine Learning techniques to speed up the natural infrastructure damage assessment process to achieve effective HADR. Accordingly, the paper demonstrates a systematic effort to achieve efficient building damage classification. First, two novel generative adversarial nets (GANs) are designed to augment data used to train the deep-learning-based classifier. Second, a contrastive learning based method using novel data structures is developed to achieve great performance. Third, by using information fusion, the classifier is effectively trained with very few training data samples for transfer learning. All the classifiers are small enough to be loaded in a smart phone or simple laptop for first responders. Based on the available overhead imagery dataset, results demonstrate data and computational efficiency with 10% of the collected data combined with a GAN reducing the time of computation from roughly half a day to about 1 hour with roughly similar classification performances.


title: Cauchy Pairs

abstract: The notion of pairable functions is introduced and some of its properties are developed. In this connection the famous Euler identity is interpreted as a property of certain pairable functions and finite cyclic groups.


title: Conditional Invertible Neural Networks for Medical Imaging

abstract: Over the last years, deep learning methods have become an increasingly popular choice to solve tasks from the field of inverse problems. Many of these new data-driven methods have produced impressive results, although most only give point estimates for the reconstruction. However, especially in the analysis of ill-posed inverse problems, the study of uncertainties is essential. In our work, we apply generative flow-based models based on invertible neural networks to two challenging medical imaging tasks, i.e. low-dose computed tomography and accelerated medical resonance imaging. We test different architectures of invertible neural networks and provide extensive ablation studies. In most applications, a standard Gaussian is used as the base distribution for a flow-based model. Our results show that the choice of a radial distribution can improve the quality of reconstructions.


title: Active clustering for labeling training data

abstract: Gathering training data is a key step of any supervised learning task, and it is both critical and expensive. Critical, because the quantity and quality of the training data has a high impact on the performance of the learned function. Expensive, because most practical cases rely on humans-in-the-loop to label the data. The process of determining the correct labels is much more expensive than comparing two items to see whether they belong to the same class. Thus motivated, we propose a setting for training data gathering where the human experts perform the comparatively cheap task of answering pairwise queries, and the computer groups the items into classes (which can be labeled cheaply at the very end of the process). Given the items, we consider two random models for the classes: one where the set partition they form is drawn uniformly, the other one where each item chooses its class independently following a fixed distribution. In the first model, we characterize the algorithms that minimize the average number of queries required to cluster the items and analyze their complexity. In the second model, we analyze a specific algorithm family, propose as a conjecture that they reach the minimum average number of queries and compare their performance to a random approach. We also propose solutions to handle errors or inconsistencies in the experts' answers.


title: New findings for the old problem: Exact solutions for domain walls in coupled real Ginzburg-Landau equations

abstract: This work reports new exact solutions for domain-wall (DW) states produced by a system of coupled real Ginzburg-Landau (GL) equations which model patterns in thermal convection, optics, and Bose-Einstein condensates (BECs). An exact solution for symmetric DW was known for a single value of the cross-interaction coefficient, G = 3 (defined so that its self-interaction counterpart is 1). Here an exact asymmetric DW is obtained for the system in which the diffusion term is absent in one component. It exists for all G > 1. Also produced is an exact solution for DW in the symmetric real-GL system which includes linear coupling. In addition, an effect of a trapping potential on the DW is considered, which is relevant to the case of BEC. In a system of three GL equations, an exact solution is obtained for a composite state including a two-component DW and a localized state in the third component. Bifurcations which create two lowest composite states are identified too. Lastly, exact solutions are found for the system of real GL equations for counterpropagating waves, which represent a sink or source of the waves, as well as for a system of three equations which includes a standing localized component.


title: Solving eigenvalue PDEs of metastable diffusion processes using artificial neural networks

abstract: In this paper, we consider the eigenvalue PDE problem of the infinitesimal generators of metastable diffusion processes. We propose a numerical algorithm based on training artificial neural networks for solving the leading eigenvalues and eigenfunctions of such high-dimensional eigenvalue problem. The algorithm is able to find multiple leading eigenpairs by solving a single training task. It is useful in understanding the dynamical behaviors of metastable processes on large timescales. We demonstrate the capability of our algorithm on a high-dimensional model problem, and on the simple molecular system alanine dipeptide.


title: Model based Multi-agent Reinforcement Learning with Tensor Decompositions

abstract: A challenge in multi-agent reinforcement learning is to be able to generalize over intractable state-action spaces. Inspired from Tesseract [Mahajan et al., 2021], this position paper investigates generalisation in state-action space over unexplored state-action pairs by modelling the transition and reward functions as tensors of low CP-rank. Initial experiments on synthetic MDPs show that using tensor decompositions in a model-based reinforcement learning algorithm can lead to much faster convergence if the true transition and reward functions are indeed of low rank.


title: Doubly Robust Criterion for Causal Inference

abstract: The semiparametric estimation approach, which includes inverse-probability-weighted and doubly robust estimation using propensity scores, is a standard tool in causal inference, and it is rapidly being extended in various directions. On the other hand, although model selection is indispensable in statistical analysis, an information criterion for selecting an appropriate regression structure has just started to be developed. In this paper, based on the original definition of Akaike information criterion (AIC; \citealt{Aka73}), we derive an AIC-type criterion for propensity score analysis. Here, we define a risk function based on the Kullback-Leibler divergence as the cornerstone of the information criterion and treat a general causal inference model that is not necessarily a linear one. The causal effects to be estimated are those in the general population, such as the average treatment effect on the treated or the average treatment effect on the untreated. In light of the fact that this field attaches importance to doubly robust estimation, which allows either the model of the assignment variable or the model of the outcome variable to be wrong, we make the information criterion itself doubly robust so that either one can be wrong and it will still be an asymptotically unbiased estimator of the risk function. In simulation studies, we compare the derived criterion with an existing criterion obtained from a formal argument and confirm that the former outperforms the latter. Specifically, we check that the divergence between the estimated structure from the derived criterion and the true structure is clearly small in all simulation settings and that the probability of selecting the true or nearly true model is clearly higher. Real data analyses confirm that the results of variable selection using the two criteria differ significantly.


title: An ubiquitous three-term recurrence relation

abstract: We solve an eigenvalue equation that appears in several papers about a wide range of physical problems. The Frobenius method leads to a three-term recurrence relation for the coefficients of the power series that, under suitable truncation, yields exact analytical eigenvalues and eigenfunctions for particular values of a model parameter. From these solutions some researchers have derived a variety of predictions like allowed angular frequencies, allowed field intensities and the like. We also solve the eigenvalue equation numerically by means of the variational Rayleigh-Ritz method and compare the resulting eigenvalues with those provided by the truncation condition. In this way we prove that those physical predictions are merely artifacts of the truncation condition.


title: The Small-Sized Telescopes for the Southern Site of the Cherenkov Telescope Array

abstract: The Cherenkov Telescope Array (CTA) will use three telescope sizes to efficiently detect cosmic gamma rays in the energy range from several tens of GeV to hundreds of TeV. The Small-Sized Telescopes (SSTs) will form the largest section of the array, covering an area of many square kilometres on the CTA southern site in Paranal, Chile. Up to 70 SSTs will be implemented by an international consortium of institutes and teams as an in-kind contribution to the CTA Observatory. The SSTs will provide unprecedented sensitivity to gamma rays above 1 TeV and the highest angular resolution of any instrument above the hard X-ray band. CTA has recently finalised the technology that will be used for the SSTs: the telescopes will be a dual-reflector design with a primary reflector of ~4 m diameter, equipped with an SiPM-based camera with full waveform readout from $\sim$2000 channels covering a $\sim$9$^\circ$ field of view. The Schwarzschild-Couder optical configuration leads to a small plate-scale, and consequently a compact, cost-efficient camera. In this contribution, we describe the experience gained operating telescope and camera prototypes during the CTA preparatory phase, and the development of the final SST design.


title: On Competitive Permutations for Set Cover by Intervals

abstract: We revisit the problem of computing an optimal partial cover of points by intervals. We show that the greedy algorithm computes a permutation $\Pi = \pi_1, \pi_2,\ldots$ of the intervals that is $3/4$-competitive for any prefix of $k$ intervals. That is, for any $k$, the intervals $\pi_1 \cup \cdots \cup \pi_k$ covers at least $3/4$-fraction of the points covered by the optimal solution using $k$ intervals. We also provide an approximation algorithm that, in $O(n + m/\varepsilon)$ time, computes a cover by $(1+\varepsilon)k$ intervals that is as good as the optimal solution using $k$ intervals, where $n$ is the number of input points, and $m$ is the number of intervals (we assume here the input is presorted). Finally, we show a counter example illustrating that the optimal solutions for set cover do not have the diminishing return property -- that is, the marginal benefit from using more sets is not monotonically decreasing. Fortunately, the diminishing returns does hold for intervals.


title: HSVI for zs-POSGs using Concavity, Convexity and Lipschitz Properties

abstract: Dynamic programming and heuristic search are at the core of state-of-the-art solvers for sequential decision-making problems. In partially observable or collaborative settings (\eg, POMDPs and Dec-POMDPs), this requires introducing an appropriate statistic that induces a fully observable problem as well as bounding (convex) approximators of the optimal value function. This approach has succeeded in some subclasses of 2-player zero-sum partially observable stochastic games (zs-POSGs) as well, but failed in the general case despite known concavity and convexity properties, which only led to heuristic algorithms with poor convergence guarantees. We overcome this issue, leveraging on these properties to derive bounding approximators and efficient update and selection operators, before deriving a prototypical solver inspired by HSVI that provably converges to an $\epsilon$-optimal solution in finite time, and which we empirically evaluate. This opens the door to a novel family of promising approaches complementing those relying on linear programming or iterative methods.


title: Fully device-independent quantum key distribution using synchronous correlations

abstract: We derive a device-independent quantum key distribution protocol based on synchronous correlations and their Bell inequalities. This protocol offers several advantages over other device-independent schemes including symmetry between the two users and no need for preshared randomness. We close a "synchronicity" loophole by showing that an almost synchronous correlation inherits the self-testing property of the associated synchronous correlation. We also pose a new security assumption that closes the "locality" (or "causality") loophole: an unbounded adversary with even a small uncertainty about the users' choice of measurement bases cannot produce any almost synchronous correlation that approximately maximally violates a synchronous Bell inequality.


title: Our Fundamental Physical Space: An Essay on the Metaphysics of the Wave Function

abstract: The mathematical structure of realist quantum theories has given rise to a debate about how our ordinary 3-dimensional space is related to the 3N-dimensional configuration space on which the wave function is defined. Which of the two spaces is our (more) fundamental physical space? I review the debate between 3N-Fundamentalists and 3D-Fundamentalists and evaluate it based on three criteria. I argue that when we consider which view leads to a deeper understanding of the physical world, especially given the deeper topological explanation from the unordered configurations to the Symmetrization Postulate, we have strong reasons in favor of 3D-Fundamentalism. I conclude that our evidence favors the view that our fundamental physical space in a quantum world is 3-dimensional rather than 3N-dimensional. I outline lines of future research where the evidential balance can be restored or reversed. Finally, I draw lessons from this case study to the debate about theoretical equivalence.


title: FacTeR-Check: Semi-automated fact-checking through Semantic Similarity and Natural Language Inference

abstract: Our society produces and shares overwhelming amounts of information through Online Social Networks (OSNs). Within this environment, misinformation and disinformation have proliferated, becoming a public safety concern in most countries. Allowing the public and professionals to efficiently find reliable evidences about the factual veracity of a claim is a crucial step to mitigate this harmful spread. To this end, we propose FacTeR-Check, a multilingual architecture for semi-automated fact-checking that can be used for either applications designed for the general public and by fact-checking organisations. FacTeR-Check enables retrieving fact-checked information, unchecked claims verification and tracking dangerous information over social media. This architectures involves several modules developed to evaluate semantic similarity, to calculate natural language inference and to retrieve information from Online Social Networks. The union of all these components builds a semi-automated fact-checking tool able of verifying new claims, to extract related evidence, and to track the evolution of a hoax on a OSN. While individual modules are validated on related benchmarks (mainly MSTS and SICK), the complete architecture is validated using a new dataset called NLI19-SP that is publicly released with COVID-19 related hoaxes and tweets from Spanish social media. Our results show state-of-the-art performance on the individual benchmarks, as well as producing a useful analysis of the evolution over time of 61 different hoaxes.


title: Addressing delayed case reporting in infectious disease forecast modeling

abstract: Infectious disease forecasting is of great interest to the public health community and policymakers, since forecasts can provide insight into disease dynamics in the near future and inform interventions. Due to delays in case reporting, however, forecasting models may often underestimate the current and future disease burden. In this paper, we propose a general framework for addressing reporting delay in disease forecasting efforts with the goal of improving forecasts. We propose strategies for leveraging either historical data on case reporting or external internet-based data to estimate the amount of reporting error. We then describe several approaches for adapting general forecasting pipelines to account for under- or over-reporting of cases. We apply these methods to address reporting delay in data on dengue fever cases in Puerto Rico from 1990 to 2009 and to reports of influenza-like illness (ILI) in the United States between 2010 and 2019. Through a simulation study, we compare method performance and evaluate robustness to assumption violations. Our results show that forecasting accuracy and prediction coverage almost always increase when correction methods are implemented to address reporting delay. Some of these methods required knowledge about the reporting error or high quality external data, which may not always be available. Provided alternatives include excluding recently-reported data and performing sensitivity analysis. This work provides intuition and guidance for handling delay in disease case reporting and may serve as a useful resource to inform practical infectious disease forecasting efforts.


title: Differential Deep Detection in Massive MIMO With One-Bit ADC

abstract: This article presents a differential detection scheme for the uplink of a massive MIMO system that employs one-bit quantizers on each receive antenna. We focus on the detection of differential amplitude and phase shift keying symbols and we use the Bussgang theorem to express the quantized received signal in terms of quantized signals received during previous channel uses. Subsequently, we derive the maximum likelihood detector for the differentially encoded amplitude and phase information symbols. We note that while the one-bit detector can decode the differentially encoded phase information symbols, it fails to decode the differentially encoded amplitude information. To decode the amplitude information, we present a one-bit variable quantization level (VQL) system and train a deep neural network to perform two-symbol differential amplitude detection. Through Monte-Carlo simulations, we empirically validate the performance of the proposed amplitude and phase detectors. The presented numerical results show that the spectral efficiency attained in one-bit differential systems is better than the spectral efficiency attained in one-bit coherent systems


title: Comparing Heuristics, Constraint Optimization, and Reinforcement Learning for an Industrial 2D Packing Problem

abstract: Cutting and Packing problems are occurring in different industries with a direct impact on the revenue of businesses. Generally, the goal in Cutting and Packing is to assign a set of smaller objects to a set of larger objects. To solve Cutting and Packing problems, practitioners can resort to heuristic and exact methodologies. Lately, machine learning is increasingly used for solving such problems. This paper considers a 2D packing problem from the furniture industry, where a set of wooden workpieces must be assigned to different modules of a trolley in the most space-saving way. We present an experimental setup to compare heuristics, constraint optimization, and deep reinforcement learning for the given problem. The used methodologies and their results get collated in terms of their solution quality and runtime. In the given use case a greedy heuristic produces optimal results and outperforms the other approaches in terms of runtime. Constraint optimization also produces optimal results but requires more time to perform. The deep reinforcement learning approach did not always produce optimal or even feasible solutions. While we assume this could be remedied with more training, considering the good results with the heuristic, deep reinforcement learning seems to be a bad fit for the given use case.


title: Connecting density fluctuations and Kirkwood-Buff integrals for finite-size systems

abstract: Kirkwood-Buff integrals (KBI) connect the microscopic structure and thermodynamic properties of liquid solutions. KBI are defined in the grand canonical ensemble and evaluated assuming the thermodynamic limit (TL). In order to reconcile analytical and numerical approaches, finite-size KBI have been proposed in the literature, resulting in two strategies to obtain their TL values from computer simulations. (i) The spatial block-analysis method in which the simulation box is divided into subdomains of volume $V$ to compute fluctuations of the number of particles. (ii) A direct integration method where a corrected radial distribution function and a kernel that accounts for the geometry of the integration subvolumes are combined to obtain KBI as a function of $V$. In this work, we propose a method that connects both strategies into a single framework. We start from the definition of finite-size KBI, including the integration subdomain and an asymptotic correction to the radial distribution function, and solve them in Fourier space where periodic boundary conditions are trivially introduced. The limit $q\to 0$, equivalent to the value of the KBI in the TL, is obtained via the spatial block-analysis method. When compared to the latter, our approach gives nearly identical results for all values of $V$. Moreover, all finite-size effect contributions (ensemble, finite-integration domains and periodic boundary conditions) are easily identifiable in the calculation. This feature allows us to analyse finite-size effects independently and extrapolate the results of a single simulation to different box sizes. To validate our approach, we investigate prototypical systems, including SPC/E water and aqueous urea mixtures.


title: The contact process with fitness on random trees

abstract: The contact process is a simple model for the spread of an infection in a structured population. We consider a variant of this process on Bienaym\'e-Galton-Watson trees, where vertices are equipped with a random fitness representing inhomogeneous transmission rates among individuals. In this paper, we establish conditions under which this inhomogeneous contact process exhibits a phase transition. We first prove that if certain mixed moments of the joint offspring and fitness distribution are finite, then the survival threshold is strictly positive. Further, we show that, if slightly different mixed moments are infinite, then this implies that there is no phase transition and the process survives with positive probability for any choice of the infection parameter. A similar dichotomy is known for the contact process on a Bienaym\'e-Galton-Watson tree. However, we show that the introduction of fitness means that we have to take into account the combined effect of fitness and offspring distribution to decide which scenario occurs.


title: Reinforcement Learning in Factored Action Spaces using Tensor Decompositions

abstract: We present an extended abstract for the previously published work TESSERACT [Mahajan et al., 2021], which proposes a novel solution for Reinforcement Learning (RL) in large, factored action spaces using tensor decompositions. The goal of this abstract is twofold: (1) To garner greater interest amongst the tensor research community for creating methods and analysis for approximate RL, (2) To elucidate the generalised setting of factored action spaces where tensor decompositions can be used. We use cooperative multi-agent reinforcement learning scenario as the exemplary setting where the action space is naturally factored across agents and learning becomes intractable without resorting to approximation on the underlying hypothesis space for candidate solutions.


title: Decoupling of the many-body effects from the electron mass in GaAs by means of reduced dimensionality

abstract: Determining the (bare) electron mass $m_0$ in crystals is often hindered by many-body effects since Fermi-liquid physics renormalises the band mass, making the observed effective mass $m^*$ depend on density. Here, we use a one-dimensional (1D) geometry to amplify the effect of interactions, forcing the electrons to form a nonlinear Luttinger liquid with separate holon and spinon bands, therefore separating the interaction effects from $m_0$. Measuring the spectral function of gated quantum wires formed in GaAs by means of magnetotunnelling spectroscopy and interpreting them using the 1D Fermi-Hubbard model, we obtain $m_0=(0.0525\pm0.0015)m_\textrm{e}$ in this material, where $m_\textrm{e}$ is the free-electron mass. By varying the density in the wires, we change the interaction parameter $r_\textrm{s}$ in the range from $\sim$1-4 and show that $m_0$ remains constant. The determined value of $m_0$ is $\sim 22$% lighter than observed in GaAs in geometries of higher dimensionality $D$ ($D>1$), consistent with the quasi-particle picture of a Fermi liquid that makes electrons heavier in the presence of interactions.


title: Skewness and Kurtosis in Stochastic Thermodynamics

abstract: The thermodynamic uncertainty relation is a prominent result in stochastic thermodynamics that provides a bound on the fluctuations of any thermodynamic flux, also known as current, in terms of the average rate of entropy production. Such fluctuations are quantified by the second moment of the probability distribution of the current. The role of higher order standardized moments such as skewness and kurtosis remains largely unexplored. We analyze the skewness and kurtosis associated with the first passage time of thermodynamic currents within the framework of stochastic thermodynamics. We develop a method to evaluate higher order standardized moments associated with the first passage time of any current. For systems with a unicyclic network of states, we conjecture upper and lower bounds on skewness and kurtosis associated with entropy production. These bounds depend on the number of states and the thermodynamic force that drives the system out of equilibrium. We show that these bounds for skewness and kurtosis do not hold for multicyclic networks. We discuss the application of our results to infer an underlying network of states.


title: Deep Reinforcement Learning for Simultaneous Sensing and Channel Access in Cognitive Networks

abstract: We consider the problem of dynamic spectrum access (DSA) in cognitive wireless networks, where only partial observations are available to the users due to narrowband sensing and transmissions. The cognitive network consists of primary users (PUs) and a secondary user (SU), which operate in a time duplexing regime. The traffic pattern for each PU is assumed to be unknown to the SU and is modeled as a finite-memory Markov chain. Since observations are partial, then both channel sensing and access actions affect the throughput. The objective is to maximize the SU's long-term throughput. To achieve this goal, we develop a novel algorithm that learns both access and sensing policies via deep Q-learning, dubbed Double Deep Q-network for Sensing and Access (DDQSA). To the best of our knowledge, this is the first paper that solves both sensing and access policies for DSA via deep Q-learning. Second, we analyze the optimal policy theoretically to validate the performance of DDQSA. Although the general DSA problem is P-SPACE hard, we derive the optimal policy explicitly for a common model of a cyclic user dynamics. Our results show that DDQSA learns a policy that implements both sensing and channel access, and significantly outperforms existing approaches.


title: Heat transfer process with solid-solid interface: Analytical and numerical solutions

abstract: This work is aimed at the study and analysis of the heat transport on a metal bar of length $L$ with a solid-solid interface. The process is assumed to be developed along one direction, across two homogeneous and isotropic materials. Analytical and numerical solutions are obtained under continuity conditions at the interface, that is a perfect assembly. The lateral side is assumed to be isolated and a constant thermal source is located at the left-boundary while the right-end stays free allowing the heat to transfer to the surrounding fluid by a convective process. The differences between the analytic solution and temperature measurements at any point on the right would indicate the presence of discontinuities. The greater these differences, the greater the discontinuity in the interface due to thermal resistances, providing a measure of its propagation from the interface and they could be modeled as temperature perturbations. The problem of interest may be described by a parabolic equation with initial, interface and boundary conditions, where the thermal properties, the conductivity and diffusivity coefficients, are piecewise constant functions. The analytic solution is derived by using Fourier methods. Special attention is given to the Sturm-Liouville problem that arises when deriving the solution, since a complicated eigenvalue equation must to be solved. Numerical simulations are conducted by using finite difference schemes where its convergence and stability properties are discussed along with physical interpretations of the results.


title: Solving the Yamabe Problem by an Iterative Method on a Small Riemannian Domain

abstract: We introduce an iterative scheme to solve the Yamabe equation $ - a\Delta_{g} u + S u = \lambda u^{p-1} $ on small domains $(\Omega,g)\subset {\mathbb R}^n$ equipped with a Riemannian metric $g$. Thus $g$ admits a conformal change to a constant scalar curvature metric. The proof does not use the traditional functional minimization. Applications to the Yamabe problem on closed manifolds, manifolds with boundary, and noncompact manifolds are given in forthcoming papers.


title: Power Minimization of Downlink Spectrum Slicing for eMBB and URLLC Users

abstract: 5G technology allows heterogeneous services to share the wireless spectrum within the same radio access network. In this context, spectrum slicing of the shared radio resources is a critical task to guarantee the performance of each service. We analyze a downlink communication serving two types of traffic: enhanced mobile broadband (eMBB) and ultra-reliable low-latency communication (URLLC). Due to the nature of low-latency traffic, the base station knows the channel state information (CSI) of the eMBB users while having statistical CSI for the URLLC users. We study the power minimization problem employing orthogonal multiple access (OMA) and non-orthogonal multiple access (NOMA) schemes. Based on this analysis, we propose a lookup table-based approach and a block coordinated descent (BCD) algorithm. We show that the BCD is optimal for the URLLC power allocation. The numerical results show that NOMA leads to lower power consumption than OMA, except when the average channel gain of the URLLC user is very high. For the latter case, the optimal approach depends on the channel condition of the eMBB user. Even when OMA attains the best performance, the gap with NOMA is negligible, showing the capability of NOMA to reduce power consumption in practically every condition.


title: Performance prediction of massively parallel computation by Bayesian inference

abstract: A performance prediction method for massively parallel computation is proposed. The method is based on performance modeling and Bayesian inference to predict elapsed time T as a function of the number of used nodes P (T=T(P)). The focus is on extrapolation for larger values of P from the perspective of application researchers. The proposed method has several improvements over the method developed in a previous paper, and application to real-symmetric generalized eigenvalue problem shows promising prediction results. The method is generalizable and applicable to many other computations.


title: Stacked polarimeters with twisted black phosphorus

abstract: The real-time, in-line analysis of light polarization is critical in optical communication networks, which suffers from the complex systems with numerous bulky opto-electro-mechanical elements tandemly arranged along optical path. Here, we propose a fiber-integrated polarimeter with nano-thickness by vertically stacking three two-dimensional (2D) materials based photodetection units. We demonstrate a self-power-calibrated, ultrafast, unambiguous detection of linear (LP) and circular polarized (CP) light according to the symmetry broken induced linear photogalvanic effects (LPGE) and circular photogalvanic effects (CPGE) in black phosphorous (BP) units, which are twistedly stacked to substitute traditional mechanical rotation of polarizers. As a demonstration, we achieve Hadamard single-pixel polarimetric imaging by the polarimeter to recognize the polarization distributions, showing potential in high-speed polarization-division-multiplexed imaging and real-time polarized endoscopy. This work provides a new strategy for next-generation ultracompact optical and optoelectronic systems, and guides a way for developing high-resolution arrayed devices with multifunctional pixels.


title: On sufficient conditions for spanning structures in dense graphs

abstract: We study structural conditions in dense graphs that guarantee the existence of vertex-spanning substructures such as Hamilton cycles. It is easy to see that every Hamiltonian graph is connected, has a perfect fractional matching and, excluding the bipartite case, contains an odd cycle. A simple consequence of the Robust Expander Theorem of K\"uhn, Osthus and Treglown tells us that any large enough graph that robustly satisfies these properties must already be Hamiltonian. Our main result generalises this phenomenon to powers of cycles and graphs of sublinear bandwidth subject to natural generalisations of connectivity, matchings and odd cycles. This answers a question of Ebsen, Maesaka, Reiher, Schacht and Sch\"ulke and solves the embedding problem that underlies multiple lines of research on sufficient conditions for spanning structures in dense graphs. As applications, we recover and establish Bandwidth Theorems in a variety of settings including Ore-type degree conditions, P\'osa-type degree conditions, deficiency-type conditions, locally dense and inseparable graphs, multipartite graphs as well as robust expanders.


title: Stability estimates for radial basis function methods applied to time-dependent hyperbolic PDEs

abstract: We derive stability estimates for three commonly used radial basis function (RBF) methods to solve hyperbolic time-dependent PDEs: the RBF generated finite difference (RBF-FD) method, the RBF partition of unity method (RBF-PUM) and Kansa's (global) RBF method. We give the estimates in the discrete $\ell_2$-norm intrinsic to each of the three methods.The results show that Kansa's method and RBF-PUM can be $\ell_2$-stable in time under a sufficiently large oversampling of the discretized system of equations. On the other hand, the RBF-FD method is not $\ell_2$-stable by construction, no matter how large the oversampling is. We show that this is due to the jumps (discontinuities) in the RBF-FD cardinal basis functions. We also provide a stabilization of the RBF-FD method that penalizes the spurious jumps. Numerical experiments show an agreement with our theoretical observations.


title: Latent Equilibrium: A unified learning theory for arbitrarily fast computation with arbitrarily slow neurons

abstract: The response time of physical computational elements is finite, and neurons are no exception. In hierarchical models of cortical networks each layer thus introduces a response lag. This inherent property of physical dynamical systems results in delayed processing of stimuli and causes a timing mismatch between network output and instructive signals, thus afflicting not only inference, but also learning. We introduce Latent Equilibrium, a new framework for inference and learning in networks of slow components which avoids these issues by harnessing the ability of biological neurons to phase-advance their output with respect to their membrane potential. This principle enables quasi-instantaneous inference independent of network depth and avoids the need for phased plasticity or computationally expensive network relaxation phases. We jointly derive disentangled neuron and synapse dynamics from a prospective energy function that depends on a network's generalized position and momentum. The resulting model can be interpreted as a biologically plausible approximation of error backpropagation in deep cortical networks with continuous-time, leaky neuronal dynamics and continuously active, local plasticity. We demonstrate successful learning of standard benchmark datasets, achieving competitive performance using both fully-connected and convolutional architectures, and show how our principle can be applied to detailed models of cortical microcircuitry. Furthermore, we study the robustness of our model to spatio-temporal substrate imperfections to demonstrate its feasibility for physical realization, be it in vivo or in silico.


title: Testing and Estimating Structural Breaks in Time Series and Panel Data in Stata

abstract: Identifying structural change is a crucial step in analysis of time series and panel data. The longer the time span, the higher the likelihood that the model parameters have changed as a result of major disruptive events, such as the 2007-2008 financial crisis and the 2020 COVID-19 outbreak. Detecting the existence of breaks, and dating them is therefore necessary not only for estimation purposes but also for understanding drivers of change and their effect on relationships. This article introduces a new community contributed command called xtbreak, which provides researchers with a complete toolbox for analysing multiple structural breaks in time series and panel data. xtbreak can detect the existence of breaks, determine their number and location, and provide break date confidence intervals. The new command is used to explore changes in the relationship between COVID-19 cases and deaths in the US using both country-level time series data and state-level panel data.


title: Thermalization without chaos in harmonic systems

abstract: Recent numerical results showed that thermalization of Fourier modes is achieved in short time-scales in the Toda model, despite its integrability and the absence of chaos. Here we provide numerical evidence that the scenario according to which chaos is irrelevant for thermalization is realized even in the simplest of all classical integrable system: the harmonic chain. We study relaxation from an atypical condition given with respect to "random" modes, showing that a thermal state with equilibrium properties is attained in short times. Such a result is independent from the orthonormal base used to represent the chain state, provided it is random.


title: On plethysms and Sylow branching coefficients

abstract: We prove a recursive formula for plethysm coefficients of the form $a^\mu_{\lambda,(m)}$, generalising results on plethysms due to Bruns--Conca--Varbaro and de Boeck--Paget--Wildon. From this we deduce a stability result and resolve two conjectures of de Boeck concerning plethysms, as well as obtain new results on Sylow branching coefficients for symmetric groups for the prime 2. Further, letting $P_n$ denote a Sylow 2-subgroup of $S_n$, we show that almost all Sylow branching coefficients of $S_n$ corresponding to the trivial character of $P_n$ are positive.


title: GenURL: A General Framework for Unsupervised Representation Learning

abstract: Unsupervised representation learning (URL), which learns compact embeddings of high-dimensional data without supervision, has made remarkable progress recently. However, the development of URLs for different requirements is independent, which limits the generalization of the algorithms, especially prohibitive as the number of tasks grows. For example, dimension reduction methods, t-SNE, and UMAP optimize pair-wise data relationships by preserving the global geometric structure, while self-supervised learning, SimCLR, and BYOL focus on mining the local statistics of instances under specific augmentations. To address this dilemma, we summarize and propose a unified similarity-based URL framework, GenURL, which can smoothly adapt to various URL tasks. In this paper, we regard URL tasks as different implicit constraints on the data geometric structure that help to seek optimal low-dimensional representations that boil down to data structural modeling (DSM) and low-dimensional transformation (LDT). Specifically, DMS provides a structure-based submodule to describe the global structures, and LDT learns compact low-dimensional embeddings with given pretext tasks. Moreover, an objective function, General Kullback-Leibler divergence (GKL), is proposed to connect DMS and LDT naturally. Comprehensive experiments demonstrate that GenURL achieves consistent state-of-the-art performance in self-supervised visual learning, unsupervised knowledge distillation (KD), graph embeddings (GE), and dimension reduction.


title: Systematic errors induced by the elliptical power-law model in galaxy-galaxy strong lens modeling

abstract: The elliptical power-law (EPL) model of the mass in a galaxy is widely used in strong gravitational lensing analyses. However, the distribution of mass in real galaxies is more complex. We quantify the biases due to this model mismatch by simulating and then analysing mock {\it Hubble Space Telescope} imaging of lenses with mass distributions inferred from SDSS-MaNGA stellar dynamics data. We find accurate recovery of source galaxy morphology, except for a slight tendency to infer sources to be more compact than their true size. The Einstein radius of the lens is also robustly recovered with 0.1% accuracy, as is the global density slope, with 2.5% relative systematic error, compared to the 3.4% intrinsic dispersion. However, asymmetry in real lenses also leads to a spurious fitted `external shear' with typical strength, $\gamma_{\rm ext}=0.015$. Furthermore, time delays inferred from lens modelling without measurements of stellar dynamics are typically underestimated by $\sim$5%. Using such measurements from a sub-sample of 37 lenses would bias measurements of the Hubble constant $H_0$ by $\sim$9%. Although this work is based on a particular set of MaNGA galaxies, and the specific value of the detected biases may change for another set of strong lenses, our results strongly suggest the next generation cosmography needs to use more complex lens mass models.


title: V-Learning -- A Simple, Efficient, Decentralized Algorithm for Multiagent RL

abstract: A major challenge of multiagent reinforcement learning (MARL) is the curse of multiagents, where the size of the joint action space scales exponentially with the number of agents. This remains to be a bottleneck for designing efficient MARL algorithms even in a basic scenario with finitely many states and actions. This paper resolves this challenge for the model of episodic Markov games. We design a new class of fully decentralized algorithms -- V-learning, which provably learns Nash equilibria (in the two-player zero-sum setting), correlated equilibria and coarse correlated equilibria (in the multiplayer general-sum setting) in a number of samples that only scales with $\max_{i\in[m]} A_i$, where $A_i$ is the number of actions for the $i^{\rm th}$ player. This is in sharp contrast to the size of the joint action space which is $\prod_{i=1}^m A_i$. V-learning (in its basic form) is a new class of single-agent RL algorithms that convert any adversarial bandit algorithm with suitable regret guarantees into a RL algorithm. Similar to the classical Q-learning algorithm, it performs incremental updates to the value functions. Different from Q-learning, it only maintains the estimates of V-values instead of Q-values. This key difference allows V-learning to achieve the claimed guarantees in the MARL setting by simply letting all agents run V-learning independently.


title: Minimal Integral Models for Principal Series Weil Characters

abstract: We prove a conjecture of Udo Riese about the minimal ring of definition for principal series Weil characters of $\text{SL}_2(\mathbb{F}_p)$, for $p$ an odd prime. More precisely, we show that the $(p+1)/2$-dimensional Weil characters can be realized over the ring of integers of $\mathbb{Q}(\sqrt{\varepsilon p})$, where $\varepsilon =(-1)^{(p-1)/2}$, and we provide explicit integral models over these quadratic rings. We do so by studying the Galois action on the integral models of Weil characters recently discovered by Yilong Wang.


title: Semi-Lagrangian 4d, 5d, and 6d kinetic plasma simulation on large scale GPU equipped supercomputer

abstract: Running kinetic plasma physics simulations using grid-based solvers is very demanding both in terms of memory as well as computational cost. This is primarily due to the up to six-dimensional phase space and the associated unfavorable scaling of the computational cost as a function of grid spacing (often termed the curse of dimensionality). In this paper, we present 4d, 5d, and 6d simulations of the Vlasov--Poisson equation with a split-step semi-Lagrangian discontinuous Galerkin scheme on graphic processing units (GPUs). The local communication pattern of this method allows an efficient implementation on large-scale GPU-based systems and emphasizes the importance of considering algorithmic and high-performance computing aspects in unison. We demonstrate a single node performance above 2 TB/s effective memory bandwidth (on a node with 4 A100 GPUs) and show excellent scaling (parallel efficiency between 30% and 67%) for up to 1536 A100 GPUs on JUWELS Booster.


title: Improved Nuclear Physics Near $A=61$ Refines Urca Neutrino Luminosities in Accreted Neutron Star Crusts

abstract: We performed a Penning trap mass measurement of $^{61}{\rm Zn}$ at the National Superconducting Cyclotron Laboratory and NuShellX calculations of the $^{61}{\rm Zn}$ and $^{62}{\rm Ga}$ structure using the GXPF1A Hamiltonian to obtain improved estimates of the $^{61}{\rm Zn}(p,\gamma)^{62}{\rm Ga}$ and $^{60}{\rm Cu}(p,\gamma)^{61}{\rm Zn}$ reaction rates. Surveying astrophysical conditions for type-I X-ray bursts with the code MESA, implementing our improved reaction rates, and taking into account updated nuclear masses for $^{61}{\rm V}$ and $^{61}{\rm Cr}$ from the recent literature, we refine the neutrino luminosity from the important mass number $A=61$ urca cooling source in accreted neutron star crusts. This improves our understanding of the thermal barrier between deep heating in the crust and the shallow depths where extra heat is needed to explain X-ray superbursts, as well as the expected signature of crust urca neutrino emission in light curves of cooling transients.


title: Stochastic transport equations with unbounded divergence

abstract: We study in this article the existence and uniqueness of solutions to a class of stochastic transport equations with irregular coefficients and unbounded divergence. In the first result we assume the drift is $L^{2}([0,T] \times \R^{d})\cap L^{\infty}([0,T] \times \R^{d})$ and the divergence is the locally integrable. In the second result we show that the smoothing acts as a selection criterion when the drift is in $L^{2}([0,T] \times \R^{d})\cap L^{\infty}([0,T] \times \R^{d})$ without any condition on the divergence.


title: Bugs in Quantum Computing Platforms: An Empirical Study

abstract: The interest in quantum computing is growing, and with it, the importance of software platforms to develop quantum programs. Ensuring the correctness of such platforms is important, and it requires a thorough understanding of the bugs they typically suffer from. To address this need, this paper presents the first in-depth study of bugs in quantum computing platforms. We gather and inspect a set of 223 real-world bugs from 18 open-source quantum computing platforms. Our study shows that a significant fraction of these bugs (39.9%) are quantum-specific, calling for dedicated approaches to prevent and find them. The bugs are spread across various components, but quantum-specific bugs occur particularly often in components that represent, compile, and optimize quantum programming abstractions. Many quantum-specific bugs manifest through unexpected outputs, rather than more obvious signs of misbehavior, such as crashes. Finally, we present a hierarchy of recurrent bug patterns, including ten novel, quantum-specific patterns. Our findings not only show the importance and prevalence bugs in quantum computing platforms, but they help developers to avoid common mistakes and tool builders to tackle the challenge of preventing, finding, and fixing these bugs.


title: Relative stable pairs and a non-Calabi-Yau wall crossing

abstract: Let $Y$ be a smooth projective threefold and let $f:Y\to X$ be a birational map with $Rf_*\mathcal{O}_Y=\mathcal{O}_X$. When $Y$ is Calabi-Yau, Bryan-Steinberg defined enumerative invariants associated to such maps called $f$-relative stable (or Bryan-Steinberg) invariants. When $X$ has Gorenstein singularities and $f$ has relative dimension one, they compared these invariants to the Donaldson-Thomas, or equivalently the Pandharipande-Thomas invariants of $Y$. We define Bryan-Steinberg invariants for maps $f$ as above without assuming that $Y$ is Calabi-Yau. For $X$ with Gorenstein and rational singularities, $f$ of relative dimension one, and for insertions from $X$ and arbitrary descendant levels, we conjecture a relation between the generating functions of Bryan-Steinberg and Pandharipande-Thomas invariants of $Y$. We check the conjecture for the contraction $f: Y\to X$ of a rational curve $C$ with normal bundle $N_{C/Y}\cong \mathcal{O}_C(-1)^{\oplus 2}$ using degeneration and localization techniques to reduce to a Calabi-Yau situation, which we then treat using Joyce's motivic Hall algebra.


title: An expanded ultraluminous X-ray source catalogue

abstract: Context. ULXs ($L_\textrm{X}>10^{39}$ erg s$^{-1}$) are excellent probes of accretion physics, star formation and IMBHs searches. As the sample size of X-ray data from modern observatories increases, producing extensive catalogues of ULXs and studying their collective properties is a possibility and a priority. Aims. We build a ULX catalogue based on one of the latest XMM-Newton releases, 4XMM-DR9, and the galaxy catalogue HECATE and study whether the properties of the expanded XMM-Newton ULX population are consistent with previous findings. Methods. We perform cross-matching between XMM-Newton sources and HECATE objects to identify host galaxies, and flag known interlopers in external catalogues and databases. Manual inspection of image data from PanSTARRS1 and the NASA/IPAC database is occasionally performed. We use distance and luminosity arguments to identify ULX candidates. Spectral, abundance and variability properties of the candidates are studied from 4XMM-DR9, HECATE and 4XMM-DR9s parameters. Results. We identify 779 ULX candidates. In spiral galaxies the number of ULX candidates per star forming rate is consistent with previous studies, and a significant ULX population in elliptical and lenticular is also found. Candidates hosted by late-type galaxies are more abundant, present harder spectra and undergo more and more extreme inter-observation variability than the early-type hosted ones. Around 30 candidates with $L_\textrm{X}>10^{41}$ erg s$^{-1}$ are also identified. Conclusions. Our results on the spectral and abundance properties of ULXs confirm the findings made by previous studies based on XMM-Newton and Chandra data, while our population-scale study on variability properties is unprecedented. Our study provides limited insight on the properties of the brightest ULX candidates due to the small sample size.


title: The role of representational conventions in assessing the empirical significance of symmetries

abstract: This paper explicates the direct empirical significance (DES) of symmetries in gauge theory, with comparisons to classical mechanics. Given a physical system composed of subsystems, such significance is to be awarded to physical differences of the composite system that arise from symmetries acting solely on its subsystems. So my overarching main question is: can DES be associated to the local gauge symmetries, acting solely on subsystems? In local gauge theories, any quantity with physical significance must be a gauge-invariant quantity. To attack the question of DES from this gauge-invariant angle, we require a split of the state into its physical and its representational content: a split that is relative to a representational convention, or a gauge-fixing. Using this method, we propose a rigorous definition of DES, valid for any state. This definition fills the gaps in influential previous construals of DES. In particular, Wallace's need to specialize to `generic' states is explained and dispensed with.


title: A general framework for unbiased tests of gravity using galaxy clusters

abstract: We present a Markov chain Monte Carlo pipeline which can be used for unbiased large-scale tests of gravity using galaxy cluster observations. The pipeline, which currently uses cluster number counts to constrain the present-day background scalar field $f_{R0}$ of Hu-Sawicki $f(R)$ gravity, fully accounts for the effects of the fifth force on cluster properties including the dynamical mass, the halo concentration and the observable-mass scaling relations. This is achieved using general models which have been calibrated over a wide and continuous mass range ($10^{11}M_{\odot}\lesssim M\lesssim10^{15}M_{\odot}$) using a large suite of cosmological simulations, including the first to simultaneously incorporate both screened modified gravity and full baryonic physics. We show, using mock cluster catalogues, that an incomplete treatment of the observable-mass scaling relations in $f(R)$ gravity, which do not necessarily follow the usual power-law behaviour, can lead to unbiased and imprecise constraints. It is therefore essential to fully account for these effects in future cosmological tests of gravity that will make use of vast cluster catalogues from ongoing and upcoming galaxy surveys. Our constraint framework can be easily extended to other gravity models; to demonstrate this, we have carried out a similar modelling of cluster properties in the normal-branch Dvali-Gabadadze-Porrati model (nDGP), which features a very different screening mechanism. Using our full-physics simulations, we also study the angular power spectra of the thermal and kinetic Sunyaev-Zel'dovich effects in $f(R)$ gravity and nDGP, and demonstrate the potential for precise constraints of gravity using data from upcoming CMB experiments. Finally, we present a retuned baryonic physics model, based on the IllustrisTNG model, which can be used for full-physics simulations within large cosmological volumes.


title: DreamerPro: Reconstruction-Free Model-Based Reinforcement Learning with Prototypical Representations

abstract: Top-performing Model-Based Reinforcement Learning (MBRL) agents, such as Dreamer, learn the world model by reconstructing the image observations. Hence, they often fail to discard task-irrelevant details and struggle to handle visual distractions. To address this issue, previous work has proposed to contrastively learn the world model, but the performance tends to be inferior in the absence of distractions. In this paper, we seek to enhance robustness to distractions for MBRL agents. Specifically, we consider incorporating prototypical representations, which have yielded more accurate and robust results than contrastive approaches in computer vision. However, it remains elusive how prototypical representations can benefit temporal dynamics learning in MBRL, since they treat each image independently without capturing temporal structures. To this end, we propose to learn the prototypes from the recurrent states of the world model, thereby distilling temporal structures from past observations and actions into the prototypes. The resulting model, DreamerPro, successfully combines Dreamer with prototypes, making large performance gains on the DeepMind Control suite both in the standard setting and when there are complex background distractions. Code available at https://github.com/fdeng18/dreamer-pro .


title: IndoNLI: A Natural Language Inference Dataset for Indonesian

abstract: We present IndoNLI, the first human-elicited NLI dataset for Indonesian. We adapt the data collection protocol for MNLI and collect nearly 18K sentence pairs annotated by crowd workers and experts. The expert-annotated data is used exclusively as a test set. It is designed to provide a challenging test-bed for Indonesian NLI by explicitly incorporating various linguistic phenomena such as numerical reasoning, structural changes, idioms, or temporal and spatial reasoning. Experiment results show that XLM-R outperforms other pre-trained models in our data. The best performance on the expert-annotated data is still far below human performance (13.4% accuracy gap), suggesting that this test set is especially challenging. Furthermore, our analysis shows that our expert-annotated data is more diverse and contains fewer annotation artifacts than the crowd-annotated data. We hope this dataset can help accelerate progress in Indonesian NLP research.


title: Ostrogradsky-Hamilton approach to geodetic brane gravity

abstract: We develop the Ostrogradsky-Hamilton formalism for geodetic brane gravity, described by the Regge-Teitelboim geometric model in higher codimension. We treat this gravity theory as a second-order derivative theory, based on the extrinsic geometric structure of the model. As opposed to previous treatments of geodetic brane gravity, our Lagrangian is linearly dependent on second-order time derivatives of the field variables, the embedding functions. The difference resides in a boundary term in the action, usually discarded. Certainly, this suggests applying an appropriate Ostrogradsky-Hamiltonian approach to this type of theories. The price to pay for this choice is the appearance of second class constraints. We determine the full set of phase space constraints, as well as the gauge transformations they generate in the reduced phase space. Additionally, we compute the algebra of constraints and explain its physical content. In the same spirit, we deduce the counting of the physical degrees of freedom. We comment briefly on the naive formal canonical quantization emerging from our development.


title: Tetraquark channels with $\bar b b$ pair in the static limit

abstract: Belle experiment discovered two hadrons with exotic quark content $Z_b^+\simeq \bar bb \bar du$. We present a lattice study of the $\bar bb\bar du$ systems with various quantum numbers using static bottom quarks. Only one set of quantum numbers that couples to $Z_b$ and $\Upsilon\;\pi$ was explored on the lattice before: these studies found an attractive potential between $B$ and $\bar B^*$ which leads to a bound state below the threshold. In the present study, we consider the other three sets of quantum numbers. Eigen-energies of the $\bar bb \bar du$ system are extracted as a function of separation between $b$ and $\bar b$. The resulting eigen-energies do not show any sizable deviation from non-interacting energies of the systems $\bar bb+\bar du$ and $\bar bu+\bar db$, so no significant attraction or repulsion is found. A slight exception is a small attraction between $B$ and $\bar B^*$ at small distance for the quantum number that couples to $Z_b$ and $\eta_b\;\rho$.


title: The first-principles phase diagram of monolayer nanoconfined water

abstract: Water in nanoscale cavities is ubiquitous and of central importance to everyday phenomena in geology and biology. However, the properties of nanoscale water can be remarkably different from bulk, as shown e.g., by the anomalously low dielectric constant of water in nanochannels [1], near frictionless water flow [2], or the possible existence of a square ice phase [3]. Such properties suggest that nanoconfined water could be engineered for technological applications in nanouidics [4], electrolyte materials [5], and water desalination [6]. Unfortunately, challenges in experimentally characterising water on the nanoscale and the high cost of first-principles simulations have prevented the molecular level understanding required to control the behavior of water. Here we combine a range of computational approaches to enable a first-principles level investigation of a single layer of water within a graphene-like channel. We find that monolayer water exhibits surprisingly rich and diverse phase behavior that is highly sensitive to temperature and the van der Waals pressure acting within the nanochannel. In addition to multiple molecular phases with melting temperatures varying non-monotonically by over 400 degrees with pressure, we predict a hexatic phase, which is an intermediate between a solid and a liquid, and a superionic phase with a high electrical conductivity exceeding that of battery materials. Notably, this suggests that nanoconfinement could be a promising route towards superionic behavior at easily accessible conditions.


title: Staggered Pseudo Magnetic Field in Twisted Transition Metal Dichalcogenides: Physical Origin and Experimental Consequences

abstract: Strong magnetic fields profoundly affect the quantum physics of charged particles, as seen for example by the integer and fractionally quantized Hall effects, and the fractal `Hofstadter butterfly' spectrum of electrons in the presence of a periodic potential and a magnetic field. Intrinsic physics can lead to effects equivalent to those produced by an externally applied magnetic field. Examples include the `staggered flux' phases emerging in some theories of quantum spin liquids and the Chern insulator behavior of twisted bilayer graphene when valley symmetry is broken. In this paper we show that when two layers of the transition metal dichalcogenide material WSe2 are stacked at a small relative twist angle to form a Moire bilayer, the resulting low energy physics can be understood in terms of electrons moving in a strong and tunable staggered flux. We predict experimental consequences including sign reversals of the Hall coefficient on application of an interlayer potential and spin currents appearing at sample edges and interfaces.


title: Globular Cluster Intrinsic Iron Abundance Spreads: II. Protocluster Metallicities and the Age-Metallicity Relations of Milky Way Progenitors

abstract: Intrinsic iron abundance spreads in globular clusters, although usually small, are very common, and are signatures of self enrichment: some stars within the cluster have been enriched by supernova ejecta from other stars within the same cluster. We use the Bailin (2018) self enrichment model to predict the relationship between properties of the protocluster -- its mass and the metallicity of the protocluster gas cloud -- and the final observable properties today -- its current metallicity and the internal iron abundance spread. We apply this model to an updated catalog of Milky Way globular clusters where the initial mass and/or the iron abundance spread is known to reconstruct their initial metallicities. We find that with the exception of the known anomalous bulge cluster Terzan 5 and three clusters strongly suspected to be nuclear star clusters from stripped dwarf galaxies, the model provides a good lens for understanding their iron spreads and initial metallicities. We then use these initial metallicities to construct age-metallicity relations for kinematically-identified major accretion events in the Milky Way's history. We find that using the initial metallicity instead of the current metallicity does not alter the overall picture of the Milky Way's history, since the difference is usually small, but does provide information that can help distinguish which accretion event some individual globular clusters with ambiguous kinematics should be associated with, and points to potential complexity within the accretion events themselves.


title: Boundary metrics on soliton moduli spaces

abstract: The geodesic approximation is a powerful method for studying the dynamics of BPS solitons. However, there are systems, such as BPS monopoles in three-dimensional hyperbolic space, where this approach is not applicable because the moduli space metric defined by the kinetic energy is not finite. In the case of hyperbolic monopoles, an alternative metric has been defined using the abelian connection on the sphere at infinity, but its relation to the dynamics of hyperbolic monopoles is unclear. Here this metric is placed in a more general context of boundary metrics on soliton moduli spaces. Examples are studied in systems in one and two space dimensions, where it is much easier to compare the results with simulations of the full nonlinear field theory dynamics. It is found that geodesics of the boundary metric provide a reasonable description of soliton dynamics.


title: Autonomous Exploration Development Environment and the Planning Algorithms

abstract: Autonomous Exploration Development Environment is an open-source repository released to facilitate the development of high-level planning algorithms and integration of complete autonomous navigation systems. The repository contains representative simulation environment models, fundamental navigation modules, e.g., local planner, terrain traversability analysis, waypoint following, and visualization tools. Together with two of our high-level planner releases -- TARE planner for exploration and FAR planner for route planning, we detail usage of the three open-source repositories and share experiences in the integration of autonomous navigation systems. We use DARPA Subterranean Challenge as a use case where the repositories together form the main navigation system of the CMU-OSU Team. In the end, we discuss a few potential use cases in extended applications.


title: Generalized Hydrodynamics of the attractive Non-Linear Schroedinger Equation

abstract: We study the generalized hydrodynamics of the one-dimensional classical Non Linear Schroedinger equation in the attractive phase. We thereby show that the thermodynamic limit is entirely captured by solitonic modes and radiation is absent. Our results are derived by considering the semiclassical limit of the quantum Bose gas, where the Planck constant has a key role as a regulator of the classical soliton gas. We use our result to study adiabatic interaction changes from the repulsive to the attractive phase, observing soliton production and obtaining exact analytical results which are in excellent agreement with Monte Carlo simulations.

categories: cond-mat.stat-mech cond-mat.quant-gas cond-mat.str-el nlin.SI


title: Growing uniform planar maps face by face

abstract: We provide "growth schemes" for inductively generating uniform random $2p$-angulations of the sphere with $n$ faces, as well as uniform random simple triangulations of the sphere with $2n$ faces. In the case of $2p$-angulations, we provide a way to insert a new face at a random location in a uniform $2p$-angulation with $n$ faces in such a way that the new map is precisely a uniform $2p$-angulation with $n+1$ faces. Similarly, given a uniform simple triangulation of the sphere with $2n$ faces, we describe a way to insert two new adjacent triangles so as to obtain a uniform simple triangulation of the sphere with $2n+2$ faces. The latter is based on a new bijective presentation of simple triangulations that relies on a construction by Poulalhon and Schaeffer.