title: Oscillatory, non-progressing flows induce directed cell motion

abstract: We present a deformation-dependent propulsion phenomenon for soft particles such as cells in microchannels. It is based on a broken time reversal symmetry generated by a fast forward and slow backward motion of a fluid which does not progress on average. In both sections, soft particles deform differently and thus progress relatively to the liquid. We demonstrate this by using Lattice-Boltzmann simulations of ubiquitous red blood cells in microchannels, as well as simulations for capsules and minimal soft tissue models in unbounded Poiseuille flows. The propulsion of the soft particles depends besides the oscillation asymmetry on their size, deformation type and elasticity. This is also demonstrated by analytical calculations for a minimal model. Our findings may stimulate a rethinking of particle sorting methods. For example, healthy and malignant cells often differ in their elasticity. With the proposed method, several cell types with different deformability can be separated simultaneously without labeling or obstacles in a microfluidic device.


title: Deep Optimization of Parametric IIR Filters for Audio Equalization

abstract: This paper describes a novel Deep Learning method for the design of IIR parametric filters for automatic audio equalization. A simple and effective neural architecture, named BiasNet, is proposed to determine the IIR equalizer parameters. An output denormalization technique is used to obtain accurate tuning of the IIR filters center frequency, quality factor and gain. All layers involved in the proposed method are shown to be differentiable, allowing backpropagation to optimize the network weights and achieve, after a number of training iterations, the optimal output. The parameters are optimized with respect to a loss function based on a spectral distance between the measured and desired magnitude response, and a regularization term used to achieve a spatialization of the acoustc scene. Two scenarios with different characteristics were considered for the experimental evaluation: a room and a car cabin. The performance of the proposed method improves over the baseline techniques and achieves an almost flat band. Moreover IIR filters provide a consistently lower computational cost during runtime with respect to FIR filters.


title: Something is wrong in the state of QED

abstract: Quantum electrodynamics (QED) is considered the most accurate theory in the history of science. However, this precision is based on a single experimental value: the anomalous magnetic moment of the electron (g-factor). An examination of the history of QED reveals that this value was obtained in a very suspicious way. These suspicions include the case of Karplus & Kroll, who admitted to having lied in their presentation of the most relevant calculation in the history of QED. As we will demonstrate in this paper, the Karplus & Kroll affair was not an isolated case, but one in a long series of errors, suspicious coincidences, mathematical inconsistencies and renormalized infinities swept under the rug.


title: Mixed-effects location scale models for joint modelling school value-added effects on the mean and variance of student achievement

abstract: School value-added models are widely applied to study, monitor, and hold schools to account for school differences in student learning. The traditional model is a mixed-effects linear regression of student current achievement on student prior achievement, background characteristics, and a school random intercept effect. The latter is referred to as the school value-added score and measures the mean student covariate-adjusted achievement in each school. In this article, we argue that further insights may be gained by additionally studying the variance in this quantity in each school. These include the ability to identify both individual schools and school types that exhibit unusually high or low variability in student achievement, even after accounting for differences in student intakes. We explore and illustrate how this can be done via fitting mixed-effects location scale versions of the traditional school value-added model. We discuss the implications of our work for research and school accountability systems.


title: A Methodology to Identify Cognition Gaps in Visual Recognition Applications Based on Convolutional Neural Networks

abstract: Developing consistently well performing visual recognition applications based on convolutional neural networks, e.g. for autonomous driving, is very challenging. One of the obstacles during the development is the opaqueness of their cognitive behaviour. A considerable amount of literature has been published which describes irrational behaviour of trained CNNs showcasing gaps in their cognition. In this paper, a methodology is presented that creates worstcase images using image augmentation techniques. If the CNN's cognitive performance on such images is weak while the augmentation techniques are supposedly harmless, a potential gap in the cognition has been found. The presented worst-case image generator is using adversarial search approaches to efficiently identify the most challenging image. This is evaluated with the well-known AlexNet CNN using images depicting a typical driving scenario.


title: Proper motions of OB stars in the far Carina Arm

abstract: In large scale maps of the Galactic disc, the Carina Arm stands out as a clear spiral feature, hosting prominent star clusters and associations rich in massive stars. We study the proper motions of 4199 O and early B most likely in the far Carina Arm, at distances mainly in excess of 4 kpc from the Sun, within the sky region, $282^{\circ} < \ell < 294^{\circ}$ and $-3^{\circ} < b < +1^{\circ}$ (Galactic coordinates). The sample is constructed by extending an existing blue-selected catalogue, and cross-matching with Gaia EDR3 astrometry. The observed pattern of proper motions is modulated into a saw-tooth pattern, with full amplitude approaching 1 mas yr$^{-1}$, recurring roughly every 2--3 degrees of longitude (200--300 pc at the median OB-star distance of 5.8 kpc). Kinematic perturbation of underlying circular rotation is most likely present. The data also reveal a moving group containing $>50$ OB stars at $\ell \sim 286^{\circ}$, $b \sim -1^{\circ}.4$ behind the main run of the far arm. An analysis of relative proper motions is performed that yields an incidence of runaway O stars of at least 10\% (potentially $>20$\% when full space motions become available). To map where runaways have run away from, we set up simulations for the region that assume linear trajectories and test for trajectory impact parameter in order to identify likely ejection hot spots. We find the method currently gives good results for times of flight of up to $\sim$4 Myr. It shows convincingly that only NGC 3603 and Westerlund 2 have ejected OB stars in significant numbers. Indeed, both clusters have experienced intense spells of ejection between 0.6--0.9 and 0.5--0.8 Myr ago, respectively.


title: Overview and Design of the sPHENIX TPC

abstract: The proposed sPHENIX detector at RHIC is being designed to study with unprecedented precision the nature of the Quark Gluon Plasma, or QGP, at RHIC through measurement of jets, jet correlations, and Y particles. To reach the goal of producing a high mass-resolution reconstruction of the Y(nS) states, a compact TPC covering pseudorapidities < |1.1| is under construction to provide sPHENIX with the necessary invariant mass resolution. Similar to the ALICE-TPC upgrade, the sPHENIX TPC design presently utilizes quadruple- GEM detectors. Utilizing novel zigzag patterned readout pads, a modified SAMPA readout chip, and an optimized Ne-CF4 gas mixture immersed in a 1.4 T magnetic field to reduce ion backflow, sPHENIX will have a unique lens into the QGP. The motivation behind the design and innovations of the TPC will be presented, along with the current status of the project.


title: Applying Machine Learning to Study Fluid Mechanics

abstract: This paper provides a short overview of how to use machine learning to build data-driven models in fluid mechanics. The process of machine learning is broken down into five stages: (1) formulating a problem to model, (2) collecting and curating training data to inform the model, (3) choosing an architecture with which to represent the model, (4) designing a loss function to assess the performance of the model, and (5) selecting and implementing an optimization algorithm to train the model. At each stage, we discuss how prior physical knowledge may be embedding into the process, with specific examples from the field of fluid mechanics.


title: Models for Narrative Information: A Study

abstract: The major objective of this work is to study and report the existing ontology-driven models for narrative information. The paper aims to analyze these models across various domains. The goal of this work is to bring the relevant literature, and ontology models under one umbrella, and perform a parametric comparative study. A systematic literature review methodology was adopted for an extensive literature selection. A random stratified sampling technique was used to select the models from the literature. The findings explicate a comparative view of the narrative models across domains. The differences and similarities of knowledge representation across domains, in case of narrative information models based on ontology was identified. There are significantly fewer studies that reviewed the ontology-based narrative models. This work goes a step further by evaluating the ontologies using the parameters from narrative components. This paper will explore the basic concepts and top-level concepts in the models. Besides, this study provides a comprehensive study of the narrative theories in the context of ongoing research. The findings of this work demonstrate the similarities and differences among the elements of the ontology across domains. It also identifies the state of the art literature for ontology-based narrative information.


title: Enhancing Computational Fluid Dynamics with Machine Learning

abstract: Machine learning is rapidly becoming a core technology for scientific computing, with numerous opportunities to advance the field of computational fluid dynamics. In this Perspective, we highlight some of the areas of highest potential impact, including to accelerate direct numerical simulations, to improve turbulence closure modeling, and to develop enhanced reduced-order models. We also discuss emerging areas of machine learning that are promising for computational fluid dynamics, as well as some potential limitations that should be taken into account.


title: Two simple criterion to obtain exact controllability and stabilization of a linear family of Dispersive PDE's on a periodic domain

abstract: In this work, we use the classical moment method to find a practical and simple criterion to determine if a family of linearized Dispersive equations on a periodic domain is exactly controllable and exponentially stabilizable with any given decay rate in $H^{s}_{p}(\mathbb{T})$ with $s\in \mathbb{R}$. We apply these results to prove that the linearized Smith equation, the linearized dispersion-generalized Benjamin-Ono equation, the linearized fourth-order Schr\"odinger equation, and the Higher-order Schr\"odinger equations are exactly controllable and exponentially stabilizable with any given decay rate in $H^{s}_{p}(\mathbb{T})$ with $s\in \mathbb{R}$.


title: Collective magnetic and plasma excitations in Josephson $\psi$ junctions

abstract: We show that Josephson $\psi$ junctions with the half-metallic (HM) weak link coupled to the superconducting (S) electrodes through the ferromagnetic (F) layers host collective excitations of magnetic moment and the Josephson phase. This results in the shift of the ferromagnetic resonance frequency, anomalies in the current-voltage characteristics and the appearance of additional magnetic anisotropy in the F layers. In contrast to the previously studied S/F/S junctions, the coupling between magnetic and plasma modes emerges even in the long-wavelength limit. Such coupling is shown to enable the controllable magnetization reversal in the F layer governed by the d.c. current pulse which provides the effective mechanism for magnetic moment manipulation in the devices of superconducting spintronics.


title: Sets of fractional operators and numerical estimation of the order of convergence of a family of fractional fixed point methods

abstract: Considering the large number of fractional operators that exist, and since it does not seem that their number will stop increasing soon at the time of writing this paper, it is presented for the first time, as far as the authors know, a simplified and compact way to work the fractional calculus through the classification of fractional operators using sets. This new way of working with fractional operators, which may be called as fractional calculus of sets, allows to generalize objects of the conventional calculus such as tensor operators, the diffusion equation, the heat equation, the Taylor series of a vector-valued function, and the fixed point method in several variables which allows to generate the method known as the fractional fixed point method. It is also shown that each fractional fixed point method that generates a convergent sequence has the ability to generate an uncountable family of fractional fixed point methods that generate convergent sequences. So, it is shown one way to estimate numerically the mean order of convergence of any fractional fixed point method in a region $\Omega$ through the problem of determining the critical points of a scalar function, and it is shown how to construct a hybrid fractional iterative method to determine the critical points of a scalar function.

categories: math.NA cs.NA math.CV math.FA math.OA physics.comp-ph


title: Extending the Hong-Ou-Mandel Effect: the power of nonclassicality

abstract: We show that the parity (evenness or oddness) of a nonclassical state of light has a dominant influence on the interference effects at a balanced beam splitter, irrespective of the state initially occupying the other input mode. Specifically, the parity of the nonclassical state gives rise to destructive interference effects that result in deep valleys in the output joint number distribution of which the Hong-Ou-Mandel (HOM) effect is a limiting case. The counter-intuitive influence of even a single photon to control the output of a beam splitter illuminated by any field, be it a coherent or even a noisy thermal field, demonstrates the extraordinary power of non-classicality. The canonical example of total destructive interference of quantum amplitudes leading to the absence of coincidence counts from a 50/50 beam splitter is the celebrated HOM effect, characterized by the vanishing of the joint probability of detecting singe photons in each of the output beams. We show that this is a limiting case of more general input states upon which a 50/50 BS can create total, or near total, destructive interference of quantum amplitudes. For odd photon number input Fock states of arbitrary value n>0 we show that the joint photon number probabilities vanish when detecting identical photon numbers in each output beams. We examine the mixing of photon number states of n = 1, 2, and 3 with a CV state, such as a coherent state of arbitrary amplitude, and a thermal state. These vanishing joint probabilities form what we call a central nodal line -- a contiguous set of zeros representing complete destructive interference of quantum amplitudes. For odd or even photon number Fock states with n>1 there will be additional off-diagonal pseudo-nodal curves along which the joint photon number probabilities are either zero, or near zero, which constitute a near, but not complete, destructive interference.


title: A set with no Riesz basis of exponentials

abstract: We show that there exists a bounded subset of R such that no system of exponentials can be a Riesz basis for the corresponding Hilbert space. An additional result gives a lower bound for the Riesz constant of any putative Riesz basis of the two dimensional disk.


title: A Lie algebra of Grassmannian Dirac operators and vector variables

abstract: The Lie algebra generated by $m\ $ $p$-dimensional Grassmannian Dirac operators and $m\ $ $p$-dimensional vector variables is identified as the orthogonal Lie algebra $\mathfrak{so}(2m+1)$. In this paper, we study the space $\mathcal{P}$ of polynomials in these vector variables, corresponding to an irreducible $\mathfrak{so}(2m+1)$ representation. In particular, a basis of $\mathcal{P}$ is constructed, using various Young tableaux techniques. Throughout the paper, we also indicate the relation to the theory of parafermions.


title: Universal deterministic quantum operations in microwave quantum links

abstract: We propose a realistic setup, inspired by already existing experiments, within which we develop a general formalism for the implementation of distributed quantum gates. Mediated by a quantum link that establishes a bidirectional quantum channel between distant nodes, our proposal works both for inter- and intra node communication and handles scenarios ranging from the few to the many modes limit of the quantum link. We are able to design fast and reliable state transfer protocols in every regime of operation, which, together with a detailed description of the scattering process, allows us to engineer two sets of deterministic universal distributed quantum gates. Gates whose implementation in quantum networks does not need entanglement distribution nor measurements. By employing a realistic description of the physical setup we identify the most relevant imperfections in the quantum links as well as optimal points of operation with resulting infidelities of $1-F \approx 10^{-2}-10^{-3}$.


title: The Generic Isogeny Decomposition of the Prym Variety of a Cyclic Branched Covering

abstract: Let $f\colon S'\longrightarrow S$ be a cyclic branched covering of smooth projective surfaces over $\mathbb{C}$ whose branch locus $\Delta\subset S$ is a smooth ample divisor. Pick a very ample complete linear system $|\mathcal{H}|$ on $S$, such that the polarized surface $(S, |\mathcal{H}|)$ is not a scroll nor has rational hyperplane sections. For the general member $[C]\in|\mathcal{H}|$ consider the $\mu_{n}$-equivariant isogeny decomposition of the Prym variety $Prym(C'/C)$ of the induced covering $f\colon C'= f^{-1}(C)\longrightarrow C$:\[Prym(C'/C)\sim\prod_{d|n,\ d\neq1}\mathcal{P}_{d}(C'/C).\] We show that for the very general member $[C]\in|\mathcal{H}|$ the isogeny component $\mathcal{P}_{d}(C'/C)$ is $\mu_{d}$-simple with $End_{\mu_{d}}(\mathcal{P}_{d}(C'/C))\cong\mathbb{Z}[\zeta_{d}]$. In addition, for the non-ample case we reformulate the result by considering the identity component of the kernel of the map $\mathcal{P}_{d}(C'/C)\subset Jac(C')\longrightarrow Alb(S')$.


title: Deep Chandra observations of merging galaxy cluster ZwCl 2341+0000

abstract: Knowledge of X-ray shock and radio relic connection in merging galaxy clusters has been greatly extended in terms of both observation and theory over the last decade. ZwCl 2341+0000 is a double-relic merging galaxy cluster; previous studies have shown that half of the S relic is associated with an X-ray surface brightness discontinuity, while the other half not. The discontinuity was believed to be a shock front. Therefore, it is a mysterious case of an only partial shock-relic connection. By using the 206.5 ks deep Chandra observations, we aim to investigate the nature of the S discontinuity. Meanwhile, we aim to explore new morphological and thermodynamical features. In addition, we utilize the GMRT and JVLA images to compute radio spectral index (SI) maps. In the deep observations, the previously reported S surface brightness discontinuity is better described as a sharp change in slope or as a kink, which is likely contributed by the disrupted core of the S subcluster. The radio SI maps show spectral flattening at the SE edge of the S relic, suggesting that the location of the shock front is 640 kpc away from the kink. We update the radio shock Mach number to be $2.2\pm0.1$ and $2.4\pm0.4$ for the S and N radio relics, respectively, based on the injection SI. We also put a 3 sigma lower limit on the X-ray Mach number of the S shock to be >1.6. Meanwhile, the deep observations reveal that the N subcluster is in a perfect cone shape, with a ~400 kpc linear cold front on each side. This type of conic subcluster has been predicted by simulations but is observed here for the first time. It represents a transition stage between a blunt-body cold front and a slingshot cold front. Strikingly, we found a 400 kpc long gas trail attached to the apex of the cone, which could be due to the gas stripping. In addition, an over-pressured hot region is found in the SW flank of the cluster.


title: Exploring the Limits of Large Scale Pre-training

abstract: Recent developments in large-scale machine learning suggest that by scaling up data, model size and training time properly, one might observe that improvements in pre-training would transfer favorably to most downstream tasks. In this work, we systematically study this phenomena and establish that, as we increase the upstream accuracy, the performance of downstream tasks saturates. In particular, we investigate more than 4800 experiments on Vision Transformers, MLP-Mixers and ResNets with number of parameters ranging from ten million to ten billion, trained on the largest scale of available image data (JFT, ImageNet21K) and evaluated on more than 20 downstream image recognition tasks. We propose a model for downstream performance that reflects the saturation phenomena and captures the nonlinear relationship in performance of upstream and downstream tasks. Delving deeper to understand the reasons that give rise to these phenomena, we show that the saturation behavior we observe is closely related to the way that representations evolve through the layers of the models. We showcase an even more extreme scenario where performance on upstream and downstream are at odds with each other. That is, to have a better downstream performance, we need to hurt upstream accuracy.


title: Top-N: Equivariant set and graph generation without exchangeability

abstract: This work addresses one-shot set and graph generation, and, more specifically, the parametrization of probabilistic decoders that map a vector-shaped prior to a distribution over sets or graphs. Sets and graphs are most commonly generated by first sampling points i.i.d. from a normal distribution, and then processing these points along with the prior vector using Transformer layers or Graph Neural Networks. This architecture is designed to generate exchangeable distributions, i.e., all permutations of the generated outputs are equally likely. We however show that it only optimizes a proxy to the evidence lower bound, which makes it hard to train. We then study equivariance in generative settings and show that non-exchangeable methods can still achieve permutation equivariance. Using this result, we introduce Top-n creation, a differentiable generation mechanism that uses the latent vector to select the most relevant points from a trainable reference set. Top-n can replace i.i.d. generation in any Variational Autoencoder or Generative Adversarial Network. Experimentally, our method outperforms i.i.d. generation by 15% at SetMNIST reconstruction, by 33% at object detection on CLEVR, generates sets that are 74% closer to the true distribution on a synthetic molecule-like dataset, and generates more valid molecules on QM9.


title: LCDM halo substructure properties revealed with high resolution and large volume cosmological simulations

abstract: We investigate the structural properties, distribution and abundance of LCDM dark matter subhaloes using the Phi-4096 and Uchuu suite of N-body cosmological simulations. Thanks to the combination of their large volume, high mass resolution and superb statistics, we are able to quantify -- for the first time consistently over more than seven decades in ratio of subhalo-to-host-halo mass -- dependencies of subhalo properties with mass, maximum circular velocity, Vmax, host halo mass and distance to host halo centre. We also dissect the evolution of these dependencies over cosmic time. We provide accurate fits for the subhalo mass and velocity functions, both exhibiting decreasing power-law slopes in the expected range of values and with no significant dependence on redshift. We also find subhalo abundance to depend weakly on host halo mass. We explore the distribution of subhaloes within their hosts and its evolution over cosmic time for subhaloes located as deep as ~0.1 per cent of the host virial radius. Subhalo structural properties are codified via a concentration parameter, cV, that does not depend on any specific, pre-defined density profile and relies only on Vmax. We derive the cV-Vmax relation in the range 7-1500 km/s and find an important dependence on distance of the subhalo to the host halo centre, as already described in Molin\'e et al. (2017). Interestingly, we also find subhaloes of the same mass to be significantly more concentrated into more massive hosts. Finally, we investigate the redshift evolution of cV, and provide accurate fits that take into account all mentioned dependencies. Our results offer an unprecedented detailed characterization of the subhalo population, consistent over a wide range of subhalo and host halo masses, as well as cosmic times. Our work enables precision work in any future research involving dark matter halo substructure.


title: Thermometry of Gaussian quantum systems using Gaussian measurements

abstract: We study the problem of estimating the temperature of Gaussian systems with feasible measurements, namely Gaussian and photo-detection-like measurements. For Gaussian measurements, we develop a general method to identify the optimal measurement numerically, and derive the analytical solutions in some relevant cases. For a class of single-mode states that includes thermal ones, the optimal Gaussian measurement is either Heterodyne or Homodyne, depending on the temperature regime. This is in contrast to the general setting, in which a projective measurement in the eigenbasis of the Hamiltonian is optimal regardless of temperature. In the general multi-mode case, and unlike the general unrestricted scenario where joint measurements are not helpful for thermometry (nor for any parameter estimation task), it is open whether joint Gaussian measurements provide an advantage over local ones. We conjecture that they are not useful for thermal systems, supported by partial analytical and numerical evidence. We further show that Gaussian measurements become optimal in the limit of large temperatures, while on/off photo-detection-like measurements do it for when the temperature tends to zero. Our results therefore pave the way for effective thermometry of Gaussian quantum systems using experimentally realizable measurements.


title: Complexity, Information Geometry, and Loschmidt Echo near Quantum Criticality

abstract: We consider the Nielsen complexity ${\mathcal C}_N$, the Loschmidt echo ${\mathcal L}$, and the Fubini-Study complexity $\tau$ in the transverse XY model, following a sudden quantum quench, in the thermodynamic limit. At small times, the first two are related by ${\mathcal L} \sim e^{-{\mathcal C}_N}$. By computing a novel time-dependent quantum information metric, we show that in this regime, ${\mathcal C}_N \sim d\tau^2$, up to lowest order in perturbation. The former relation continues to hold in the same limit at large times, whereas the latter does not. Our results indicate that in the thermodynamic limit, the Nielsen complexity and the Loschmidt echo show enhanced temporal oscillations when one quenches from a close neighbourhood of the critical line, while such oscillations are notably absent when the quench is on such a line. We explain this behaviour by studying the nature of quasi-particle excitations in the vicinity of criticality. Finally, we argue that the triangle inequality for the Nielsen complexity might be violated in certain regions of the parameter space, and point out why one should be careful about the nature of the interaction Hamiltonian, while using this measure.


title: Peaks of sound velocity in two color dense QCD: quark saturation effects and semishort range correlations

abstract: We discuss stiffening of dense matter in two color QCD (QC$_2$D) where hadrons are mesons and diquark baryons. We study two models which describe a transition of matter from the Bose-Einstein-Condensation regime at low density to the Bardeen-Cooper-Schrieffer regime at high density. The first model is based on coherent states of diquarks, and the second is the Nambu-Jona-Lasinio model with diquark pairing terms. We particularly focus on how quark states are occupied as baryon density increases. We find that, due to the occupied quark levels, the ideal gas picture of diquarks breaks down at density significantly less than the density where baryon cores overlap. The saturation of quark states at low momenta stiffens equations of state. We also study the effects of interactions which depend on the quark occupation probability. We argue that equations of state become very stiff when the bulk part of the quark Fermi sea has the effective repulsion but the Fermi surface enjoys the attractive correlations. This disparity for different momentum domains is possible due to the strong channel dependence in gluon exchanges with momentum transfer of $0.2-1$ GeV. These concepts can be transferred from QC$_2$D to QCD in any numbers of colors.


title: The transformation of edge-regular and pseudo strongly regular graphs under graph operations

abstract: The graph $G$ is said to be strongly regular with parameters $(n,k,\lambda,\mu)$ if the following conditions hold: (1) each vertex has $k$ neighbours; (2) any two adjacent vertices of $G$ have $\lambda$ common neighbours; (3) any two non-adjacent vertices of $G$ have $\mu$ common neighbours. In this paper we study two weaker notions of strongly regular graphs. A graph satisfying the conditions $(1)$ and $(2)$ is called an edge-regular graph with parameters $(n,k,\lambda)$. We call a graph satisfying the conditions $(1)$ and $(3)$ a pseudo strongly regular graph with parameters $(n,k,\mu)$. In this paper we study the impact of various graph operations on edge regular graphs and pseudo strongly regular graphs.


title: CARL: A Benchmark for Contextual and Adaptive Reinforcement Learning

abstract: While Reinforcement Learning has made great strides towards solving ever more complicated tasks, many algorithms are still brittle to even slight changes in their environment. This is a limiting factor for real-world applications of RL. Although the research community continuously aims at improving both robustness and generalization of RL algorithms, unfortunately it still lacks an open-source set of well-defined benchmark problems based on a consistent theoretical framework, which allows comparing different approaches in a fair, reliable and reproducibleway. To fill this gap, we propose CARL, a collection of well-known RL environments extended to contextual RL problems to study generalization. We show the urgent need of such benchmarks by demonstrating that even simple toy environments become challenging for commonly used approaches if different contextual instances of this task have to be considered. Furthermore, CARL allows us to provide first evidence that disentangling representation learning of the states from the policy learning with the context facilitates better generalization. By providing variations of diverse benchmarks from classic control, physical simulations, games and a real-world application of RNA design, CARL will allow the community to derive many more such insights on a solid empirical foundation.


title: Notarial timestamps savings in logs management via Merkle trees and Key Derivation Functions

abstract: Nowadays log files handling imposes to ISPs (intended in their widest scope) strict normative duties apart from common technological issues. This work analyses how retention time policies and timestamping are deeply interlinked from the point of view of service providers, possibly leading to costs rise. A new schema is proposed trying to mitigate the need for third-party suppliers, enforcing cryptographic primitives well established in other fields of Information Technology but perhaps not yet widespread in logs management. The foundations of timestamping are recapped, and properties of cryptographic primitives introduced as a natural way to bypass legacy schema inefficiency and as an extra level of protection: these choices are justified by savings estimation (with regard to different ISP magnitudes) and by some basic security considerations.


title: Anomalies and phases of strongly-coupled chiral gauge theories: recent developments

abstract: After many years of investigations, our understanding of the dynamics of strongly-coupled chiral gauge theories is still quite unsatisfactory today. Conventional wisdom about strongly-coupled gauge theories, successfully applied to QCD, is not always as useful in chiral gauge theories. Recently some new ideas and techniques have been developed, which involve concepts of generalized symmetries, of gauging a discrete center symmetry, and of generalizing the 't Hooft anomaly matching constraints to include certain mixed symmetries. This new development has been applied to chiral gauge theories, leading to many interesting, sometimes quite unexpected, results. For instance, in the context of generalized Bars-Yankielowicz and generalized Georgi-Glashow models, these new types of anomalies give a rather clear indication in favor of the dynamical Higgs phase, against confining, flavor symmetric vacua. Another closely related topics is strong anomaly and the effective low-energy action representing it. It turns out that they have significant implications on the phase of chiral gauge theories, giving indications consistent with the findings based on the generalized anomalies. Some striking analogies and contrasts between the massless QCD and chiral gauge theories seem to emerge from these discussions. The aim of this work is to review these developments.


title: A novel formulation of the unintegrated gluon distribution for DIS

abstract: In this talk, we revisit inclusive DIS in the small $x$ limit and derive a new factorization formula that accounts for leading powers in both Bjorken and Regge limits. In this semi-classical description, we obtain a new unintegrated gluon distribution which encompasses both the dipole operator and the gluon Parton Distribution Function with an explicit dependence on the longitudinal momentum fraction $x$.


title: Unified Theory of Characteristic Modes: Part II -- Tracking, Losses, and FEM Evaluation

abstract: This is the second component of a two-part paper dealing with a unification of characteristic mode decomposition. This second part addresses modal tracking and losses and presents several numerical examples for both surface- and volume-based method-of-moment formulations. A new tracking algorithm based on algebraic properties of the transition matrix is developed, achieving excellent precision and requiring a very low number of frequency samples as compared to procedures previously reported in the literature. The transition matrix is further utilized to show that characteristic mode decomposition of lossy objects fails to deliver orthogonal far fields and to demonstrate how characteristic modes can be evaluated using the finite element method.


title: Revisiting closed asymptotic couples

abstract: Every discrete definable subset of a closed asymptotic couple with ordered scalar field $\boldsymbol k$ is shown to be contained in a finite-dimensional $\boldsymbol k$-linear subspace of that couple. It follows that the differential-valued field $\mathbb T$ of transseries induces more structure on its value group than what is definable in its asymptotic couple equipped with its scalar multiplication by real numbers, where this asymptotic couple is construed as a two-sorted structure with $\mathbb R$ as the underlying set for the second sort.


title: Spectral density reconstruction with Chebyshev polynomials

abstract: Accurate calculations of the spectral density in a strongly correlated quantum many-body system are of fundamental importance to study its dynamics in the linear response regime. Typical examples are the calculation of inclusive and semi-exclusive scattering cross sections in atomic nuclei and transport properties of nuclear and neutron star matter. Integral transform techniques play an important role in accessing the spectral density in a variety of nuclear systems. However, their accuracy is in practice limited by the need to perform a numerical inversion which is often ill-conditioned. In the present work we extend a recently proposed quantum algorithm which circumvents this problem. We show how to perform controllable reconstructions of the spectral density over a finite energy resolution with rigorous error estimates. An appropriate expansion in Chebyshev polynomials allows for efficient simulations also on classical computers. We apply our idea to reconstruct a simple model -- response function as a proof of principle. This paves the way for future applications in nuclear and condensed matter physics.


title: Deep learning enabled strategies for modelling of complex aperiodic plasmonic metasurfaces of arbitrary size

abstract: Optical interactions have an important impact on the optical response of nanostructures in complex environments. Accounting for interactions in large ensembles of structures requires computationally demanding numerical calculations. In particular if no periodicity can be exploited, full field simulations can become prohibitively expensive. Here we propose a method for the numerical description of aperiodic assemblies of plasmonic nanostructures. Our approach is based on dressed polarizabilities, which are conventionally very expensive to calculate, a problem which we alleviate using a deep convolutional neural network as surrogate model. We demonstrate that the method offers high accuracy with errors in the order of a percent. In cases where the interactions are predominantly short-range, e.g. for out-of-plane illumination of planar metasurfaces, it can be used to describe aperiodic metasurfaces of basically unlimited size, containing many thousands of unordered plasmonic nanostructures. We furthermore show that the model is capable to spectrally resolve coupling effects. The approach is therefore of highest interest for the field of metasurfaces. It provides significant advantages in applications like homogenization of large aperiodic planar metastructures or the design of sophisticated wavefronts at the micrometer scale, where optical interactions play a crucial role.


title: OCCASO IV. Radial Velocities and Open Cluster Kinematics

abstract: Context: Open clusters (OCs) are widely used as test particles to investigate a variety of astrophysical phenomena, from stellar evolution to Galactic evolution. Gaia and the complementary massive spectroscopic surveys are providing an unprecedented wealth of information about these systems. Aims: The Open Cluster Chemical Abundances from Spanish Observatories (OCCASO) survey aims to complement all this work by determining OCs accurate radial velocities and chemical abundances from high-resolution, R$\geq$60\,000, spectra. Methods: Radial velocities have been obtained by cross-correlating the observed spectra with a library of synthetic spectra which covers from early M to A spectral types. Results: We provide radial velocities for 336 stars, including several Gaia Benchmark Stars and objects belonging to 51 open clusters. The internal uncertainties of the derived radial velocities go from 10 m/s to 21 m/s as a function of the instrumental configuration used. The derived radial velocities, together with the Gaia proper motions, have been used to investigate the cluster membership of the observed stars. After this careful membership analysis, we obtain average velocities for 47 open clusters. To our knowledge, this is the first radial velocity determination for 5 of these clusters. Finally, the radial velocities, proper motions, distances and ages have been used to investigate the kinematics of the observed clusters and in the integration of their orbits.


title: Open bibliographic data and the Italian National Scientific Qualification: measuring coverage of academic fields

abstract: The importance of open bibliographic repositories is widely accepted by the scientific community. For evaluation processes, however, there is still some skepticism: even if large repositories of open access articles and free publication indexes exist and are continuously growing, assessment procedures still rely on proprietary databases, mainly due to the richness of the data available in these proprietary databases and the services provided by the companies they are offered by. This paper investigates the status of open bibliographic data of three of the most used open resources, namely Microsoft Academic Graph, Crossref and OpenAIRE, evaluating their potentialities as substitutes of proprietary databases for academic evaluation processes. We focused on the Italian National Scientific Qualification (NSQ), the Italian process for University Professor qualification, which uses data from commercial indexes, and investigated similarities and differences between research areas, disciplines and application roles. The main conclusion is that open datasets are ready to be used for some disciplines, among which mathematics, natural sciences, economics and statistics, even if there is still room for improvement; but there is still a large gap to fill in others - like history, philosophy, pedagogy and psychology - and a stronger effort is required from researchers and institutions.


title: Deep Learning for the Approximation of a Shape Functional

abstract: Artificial Neuronal Networks are models widely used for many scientific tasks. One of the well-known field of application is the approximation of high-dimensional problems via Deep Learning. In the present paper we investigate the Deep Learning techniques applied to Shape Functionals, and we start from the so--called Torsional Rigidity. Our aim is to feed the Neuronal Network with digital approximations of the planar domains where the Torsion problem (a partial differential equation problem) is defined, and look for a prediction of the value of Torsion. Dealing with images, our choice fell on Convolutional Neural Network (CNN), and we train such a network using reference solutions obtained via Finite Element Method. Then, we tested the network against some well-known properties involving the Torsion as well as an old standing conjecture. In all cases, good approximation properties and accuracies occurred.


title: Halos and undecidability of tensor stable positive maps

abstract: A map $\mathcal{P}$ is tensor stable positive (tsp) if $\mathcal{P}^{\otimes n}$ is positive for all $n$, and essential tsp if it is not completely positive or completely co-positive. Are there essential tsp maps? Here we prove that there exist essential tsp maps on the hypercomplex numbers. It follows that there exist bound entangled states with a negative partial transpose (NPT) on the hypercomplex, that is, there exists NPT bound entanglement in the halo of quantum states. We also prove that tensor stable positivity on the matrix multiplication tensor is undecidable, and conjecture that tensor stable positivity is undecidable. Proving this conjecture would imply existence of essential tsp maps, and hence of NPT bound entangled states.


title: A geometric multiscale model for the numerical simulation of blood flow in the human left heart

abstract: We present a new computational model for the numerical simulation of blood flow in the human left heart. To this aim, we use the Navier-Stokes equations in an Arbitrary Lagrangian Eulerian formulation to account for the endocardium motion and we model the cardiac valves by means of the Resistive Immersed Implicit Surface method. To impose a physiological displacement of the domain boundary, we use a 3D cardiac electromechanical model of the left ventricle coupled to a lumped-parameter (0D) closed-loop model of the remaining circulation. We thus obtain a one-way coupled electromechanics-fluid dynamics model in the left ventricle. To extend the left ventricle motion to the endocardium of the left atrium and to that of the ascending aorta, we introduce a preprocessing procedure according to which an harmonic extension of the left ventricle displacement is combined with the motion of the left atrium based on the 0D model. To better match the 3D cardiac fluid flow with the external blood circulation, we couple the 3D Navier-Stokes equations to the 0D circulation model, obtaining a multiscale coupled 3D-0D fluid dynamics model that we solve via a segregated numerical scheme. We carry out numerical simulations for a healthy left heart and we validate our model by showing that meaningful hemodynamic indicators are correctly reproduced.

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


title: Wasserstein distance and metric trees

abstract: We study the Wasserstein (or earthmover) metric on the space $P(X)$ of probability measures on a metric space $X$. We show that, if a finite metric space $X$ embeds stochastically with distortion $D$ in a family of finite metric trees, then $P(X)$ embeds bi-Lipschitz into $\ell^1$ with distortion $D$. Next, we re-visit the closed formula for the Wasserstein metric on finite metric trees due to Evans-Matsen \cite{EvMat}. We advocate that the right framework for this formula is real trees, and we give two proofs of extensions of this formula: one making the link with Lipschitz-free spaces from Banach space theory, the other one algorithmic (after reduction to finite metric trees).


title: Local stability of McKean-Vlasov equations arising from heterogeneous Gibbs systems using limit of relative entropies

abstract: A family of heterogeneous mean-field systems with jumps is analyzed. These systems are constructed as a Gibbs measure on block graphs. When the total number of particles goes to infinity, a law of large numbers is shown to hold in a multi-class context resulting in the weak convergence of the empirical vector towards the solution of a McKean-Vlasov system of equations. We then investigate the local stability of the limiting McKean-Vlasov system through the construction of a local Lyapunov function. We first compute the limit of adequately scaled relative entropy functions associated with the explicit stationary distribution of the N-particles system. Using a Laplace principle for empirical vectors we show that the limit takes an explicit form. Then we demonstrate that this limit satisfies a descent property which, combined with some mild assumptions shows that it is indeed a local Lyapunov function.


title: Self-Supervised Generative Style Transfer for One-Shot Medical Image Segmentation

abstract: In medical image segmentation, supervised deep networks' success comes at the cost of requiring abundant labeled data. While asking domain experts to annotate only one or a few of the cohort's images is feasible, annotating all available images is impractical. This issue is further exacerbated when pre-trained deep networks are exposed to a new image dataset from an unfamiliar distribution. Using available open-source data for ad-hoc transfer learning or hand-tuned techniques for data augmentation only provides suboptimal solutions. Motivated by atlas-based segmentation, we propose a novel volumetric self-supervised learning for data augmentation capable of synthesizing volumetric image-segmentation pairs via learning transformations from a single labeled atlas to the unlabeled data. Our work's central tenet benefits from a combined view of one-shot generative learning and the proposed self-supervised training strategy that cluster unlabeled volumetric images with similar styles together. Unlike previous methods, our method does not require input volumes at inference time to synthesize new images. Instead, it can generate diversified volumetric image-segmentation pairs from a prior distribution given a single or multi-site dataset. Augmented data generated by our method used to train the segmentation network provide significant improvements over state-of-the-art deep one-shot learning methods on the task of brain MRI segmentation. Ablation studies further exemplified that the proposed appearance model and joint training are crucial to synthesize realistic examples compared to existing medical registration methods. The code, data, and models are available at https://github.com/devavratTomar/SST.


title: Electrostatic treatment of charged interfaces in classical atomistic simulations

abstract: Artificial electrostatic potentials can be present in supercells constructed for atomistic simulations of surfaces and interfaces in ionic crystals. Treating the ions as point charges, we systematically derive an electrostatic formalism for model systems of increasing complexity, both neutral and charged, and with either open or periodic boundary conditions. This allows to correctly interpret results of classical atomistic simulations which are directly affected by the appearance of these potentials. We demonstrate our approach at the example of a strontium titanite (SrTiO$_3$) supercell containing an asymmetric tilt grain boundary. The formation energies of charged oxygen vacancies and the relaxed interface structure are calculated based on an interatomic rigid-ion potential, and the results are analyzed in consideration of the electrostatic effects.


title: Unveiling the coupling of single metallic nanoparticles to whispering-gallery microcavities

abstract: Whispering-gallery mode resonators host multiple trapped narrowband circulating optical resonances that find applications in quantum electrodynamics, optomechanics, and sensing. However, the spherical symmetry and low field leakage of dielectric microspheres make it difficult to probe their high-quality optical modes using far-field radiation. Even so, local field enhancement from metallic nanoparticles (MNPs) coupled to the resonators can interface the optical far-field and the bounded cavity modes. In this work, we study the interaction between whispering-gallery modes and MNP surface plasmons with nanometric spatial resolution by using electron-beam spectroscopies in a scanning transmission electron microscope. We show that gallery modes are induced over a selective spectral range of the nanoparticle plasmons, and additionally, their polarization can be controlled by the induced dipole moment of the MNP. Our study demonstrates a viable mechanism to effectively excite high-quality-factor whispering-gallery modes and holds potential for applications in optical sensing and light manipulation.


title: Efficient Modelling Across Time of Human Actions and Interactions

abstract: This thesis focuses on video understanding for human action and interaction recognition. We start by identifying the main challenges related to action recognition from videos and review how they have been addressed by current methods. Based on these challenges, and by focusing on the temporal aspect of actions, we argue that current fixed-sized spatio-temporal kernels in 3D convolutional neural networks (CNNs) can be improved to better deal with temporal variations in the input. Our contributions are based on the enlargement of the convolutional receptive fields through the introduction of spatio-temporal size-varying segments of videos, as well as the discovery of the local feature relevance over the entire video sequence. The resulting extracted features encapsulate information that includes the importance of local features across multiple temporal durations, as well as the entire video sequence. Subsequently, we study how we can better handle variations between classes of actions, by enhancing their feature differences over different layers of the architecture. The hierarchical extraction of features models variations of relatively similar classes the same as very dissimilar classes. Therefore, distinctions between similar classes are less likely to be modelled. The proposed approach regularises feature maps by amplifying features that correspond to the class of the video that is processed. We move away from class-agnostic networks and make early predictions based on feature amplification mechanism. The proposed approaches are evaluated on several benchmark action recognition datasets and show competitive results. In terms of performance, we compete with the state-of-the-art while being more efficient in terms of GFLOPs. Finally, we present a human-understandable approach aimed at providing visual explanations for features learned over spatio-temporal networks.


title: Optimization with Constraint Learning: A Framework and Survey

abstract: Many real-life optimization problems frequently contain one or more constraints or objectives for which there are no explicit formulas. If data is however available, these data can be used to learn the constraints. The benefits of this approach are clearly seen, however there is a need for this process to be carried out in a structured manner. This paper therefore provides a framework for Optimization with Constraint Learning (OCL) which we believe will help to formalize and direct the process of learning constraints from data. This framework includes the following steps: (i) setup of the conceptual optimization model, (ii) data gathering and preprocessing, (iii) selection and training of predictive models, (iv) resolution of the optimization model, and (v) verification and improvement of the optimization model. We then review the recent OCL literature in light of this framework, and highlight current trends, as well as areas for future research.


title: The generalized Floquet-Bloch spectrum for periodic thermodiffusive layered materials

abstract: The dynamic behaviour of periodic thermodiffusive multi-layered media excited by harmonic oscillations is studied. In the framework of linear thermodiffusive elasticity, periodic laminates, whose elementary cell is composed by an arbitrary number of layers, are considered. The generalized Floquet-Bloch conditions are imposed, and the universal dispersion relation of the composite is obtained by means of an approach based on the formal solution for a single layer together with the transfer matrix method. The eigenvalue problem associated with the dispersion equation is solved by means of an analytical procedure based on the symplecticity properties of the transfer matrix to which corresponds a palindromic characteristic polynomial, and the frequency band structure associated to wave propagating inside the medium are finally derived. The proposed approach is tested through illustrative examples where thermodiffusive multilayered structures of interest for renewable energy devices fabrication are analyzed. The effects of thermodiffusion coupling on both the propagation and attenuation of Bloch waves in these systems are investigated in detail.


title: Physics of complex systems: Discovery in the age of Goedel

abstract: From the perspective of the physics of complex systems (1) we deal with the current state of mod-ern physics including the crisis in physics demonstrated through its epistemological, psychological, economical as well as the social context; (2) considering the strength of the Goedel s Incompleteness Theorems we point out the following open questions in physics: (i) the limits of the precision of certainty, (ii) the limits of making decisions on information and (iii) the limitations of reasoning that sometimes affect progress, (3) since further advances will necessarily require synergy between physics and seemingly distinct fields - mathematics, information science, chemistry, biology, medicine, psychology, and art. We illustrate this relation by providing examples based on our research. Point (1) is discussed in Ch1 while Ch2-Ch5 encompasses the point (2) and point (3) is covered in Ch6-Ch10.


title: FacialFilmroll: High-resolution multi-shot video editing

abstract: We present FacialFilmroll, a solution for spatially and temporally consistent editing of faces in one or multiple shots. We build upon unwrap mosaic [Rav-Acha et al. 2008] by specializing it to faces. We leverage recent techniques to fit a 3D face model on monocular videos to (i) improve the quality of the mosaic for edition and (ii) permit the automatic transfer of edits from one shot to other shots of the same actor. We explain how FacialFilmroll is integrated in post-production facility. Finally, we present video editing results using FacialFilmroll on high resolution videos.


title: Adversarial Robustness Verification and Attack Synthesis in Stochastic Systems

abstract: Probabilistic model checking is a useful technique for specifying and verifying properties of stochastic systems including randomized protocols and reinforcement learning models. Existing methods rely on the assumed structure and probabilities of certain system transitions. These assumptions may be incorrect, and may even be violated by an adversary who gains control of system components. In this paper, we develop a formal framework for adversarial robustness in systems modeled as discrete time Markov chains (DTMCs). We base our framework on existing methods for verifying probabilistic temporal logic properties and extend it to include deterministic, memoryless policies acting in Markov decision processes (MDPs). Our framework includes a flexible approach for specifying structure-preserving and non structure-preserving adversarial models. We outline a class of threat models under which adversaries can perturb system transitions, constrained by an $\varepsilon$ ball around the original transition probabilities. We define three main DTMC adversarial robustness problems: adversarial robustness verification, maximal $\delta$ synthesis, and worst case attack synthesis. We present two optimization-based solutions to these three problems, leveraging traditional and parametric probabilistic model checking techniques. We then evaluate our solutions on two stochastic protocols and a collection of Grid World case studies, which model an agent acting in an environment described as an MDP. We find that the parametric solution results in fast computation for small parameter spaces. In the case of less restrictive (stronger) adversaries, the number of parameters increases, and directly computing property satisfaction probabilities is more scalable. We demonstrate the usefulness of our definitions and solutions by comparing system outcomes over various properties, threat models, and case studies.


title: Ultrahigh nitrogen-vacancy center concentration in diamond

abstract: High concentration of negatively charged nitrogen-vacancy ($\text{NV}^{-}$) centers was created in diamond single crystals containing approximately 100 ppm nitrogen using electron and neutron irradiation and subsequent thermal annealing in a stepwise manner. Continuous wave electron paramagnetic resonance (EPR) was used to determine the transformation efficiency from isolated N atoms to $\text{NV}^{-}$ centers in each production step and its highest value was as high as 17.5 %. Charged vacancies are formed after electron irradiation as shown by EPR spectra, but the thermal annealing restores the sample quality as the defect signal diminishes. We find that about 25 % of the vacancies form NVs during the annealing process. The large $\text{NV}^{-}$ concentration allows to observe orientation dependent spin-relaxation times and also the determination of the hyperfine and quadrupole coupling constants with high precision using electron spin echo (ESE) and electron-nuclear double resonance (ENDOR). We also observed the EPR signal associated with the so-called W16 centers, whose spectroscopic properties might imply a nitrogen dimer-vacancy center for its origin.


title: Rate Splitting Multiple Access for Semi-Grant-Free Transmissions

abstract: Enabled by hybrid grant-based (GB) and grant-free (GF) transmission techniques, GF users of internet of things (IoT) devices and massive machine-type communications (mMTC) meet opportunities to share wireless resources with GB users. In this paper, we propose a rate splitting multiple access (RSMA) strategy for an emerging semi-grant-free (SGF) transmission system to increase connectivity and reliability. In the proposed RSMA assisted SGF (RSMA-SGF) scheme, the GF users apply the rate splitting principle to realize distributed contentions and utilize transmit power most effectively for robust transmissions, meanwhile keeping themselves transparent to the GB user. Compared to existing non-orthogonal multiple access (NOMA) assisted SGF schemes, the RSMA-SGF scheme significantly decreases outage probability and achieves full multiuser diversity gain without restricting the GB and GF users' target rates to a limited value region. Exact expressions and asymptotic analysis for the outage probability are provided to facilitate the system performance evaluation of the proposed RSMA-SGF scheme. Computer simulation results clarify the superior outage performance of the RSMA-SGF scheme and verify the accuracy of the developed analytical results.


title: NeurWIN: Neural Whittle Index Network For Restless Bandits Via Deep RL

abstract: Whittle index policy is a powerful tool to obtain asymptotically optimal solutions for the notoriously intractable problem of restless bandits. However, finding the Whittle indices remains a difficult problem for many practical restless bandits with convoluted transition kernels. This paper proposes NeurWIN, a neural Whittle index network that seeks to learn the Whittle indices for any restless bandits by leveraging mathematical properties of the Whittle indices. We show that a neural network that produces the Whittle index is also one that produces the optimal control for a set of Markov decision problems. This property motivates using deep reinforcement learning for the training of NeurWIN. We demonstrate the utility of NeurWIN by evaluating its performance for three recently studied restless bandit problems. Our experiment results show that the performance of NeurWIN is significantly better than other RL algorithms.


title: A study of first-passage time minimization via Q-learning in heated gridworlds

abstract: Optimization of first-passage times is required in applications ranging from nanobots navigation to market trading. In such settings, one often encounters unevenly distributed noise levels across the environment. We extensively study how a learning agent fares in 1- and 2- dimensional heated gridworlds with an uneven temperature distribution. The results show certain bias effects in agents trained via simple tabular Q-learning, SARSA, Expected SARSA and Double Q-learning. While high learning rate prevents exploration of regions with higher temperature, low enough rate increases the presence of agents in such regions. The discovered peculiarities and biases of temporal-difference-based reinforcement learning methods should be taken into account in real-world physical applications and agent design.

categories: math.OC cond-mat.stat-mech cs.AI cs.LG math.DS


title: Fractional spin excitations and conductance in the spiral staircase Heisenberg ladder

abstract: We investigate theoretically the spiral staircase Heisenberg spin-$1/2$ ladder in the presence of antiferromagnetic long-range spin interactions and a uniform magnetic field. As a special case we also consider the Kondo necklace model. If the magnetizations of the two chains forming the ladder satisfy a certain resonance condition, involving interchain couplings as perturbations, the system is in a partially gapped magnetic phase hosting excitations characterized by fractional spins, whose values can be changed by the magnetic field. We show that these fractional spin excitations can be probed by spin currents in a transport setup with a spin conductance that reveals the fractionalized spin. In some special cases, the spin conductance reaches universal values in units of $(g\mu_B)^2/h$, where $g$ is the $g$-factor, $\mu_B$ the Bohr magneton, and $h$ the Planck constant. We obtain our results with the help of bosonization and numerical density matrix renormalization group methods.


title: Electrochemistry meets polymer physics: polymerized ionic liquids on an electrified electrode

abstract: Polymeric ionic liquids are emerging polyelectrolyte materials for modern electrochemical applications. In this paper, we propose a self-consistent field theory of the polymeric ionic liquid on a charged conductive electrode. Taking into account the conformation entropy of rather long polymerized cations within the Lifshitz theory and electrostatic and excluded volume interactions of ionic species within the mean-field approximation, we obtain a system of self-consistent field equations for the local electrostatic potential and average concentrations of monomeric units and counterions. We solve these equations in the linear approximation for the cases of a point-like charge and a flat infinite uniformly charged electrode immersed in a polymeric ionic liquid and derive analytical expressions for local ionic concentrations and electrostatic potential, and derive an analytical expression for the linear differential capacitance of the electric double layer. We also find a numerical solution to the self-consistent field equations for two types of boundary conditions for the local polymer concentration on the electrode, corresponding to the cases of the specific adsorption absence (indifferent surface) and strong short-range repulsion of the monomeric units near the charged surface (hard wall case). For both cases, we investigate the behavior of differential capacitance as a function of applied voltage for a pure polymeric ionic liquid and a polymeric ionic liquid dissolved in a polar organic solvent. We observe that the differential capacitance profile shape is strongly sensitive to the adopted boundary condition for the local polymer concentration on the electrode.


title: A space-time multiscale mortar mixed finite element method for parabolic equations

abstract: We develop a space-time mortar mixed finite element method for parabolic problems. The domain is decomposed into a union of subdomains discretized with non-matching spatial grids and asynchronous time steps. The method is based on a space-time variational formulation that couples mixed finite elements in space with discontinuous Galerkin in time. Continuity of flux (mass conservation) across space-time interfaces is imposed via a coarse-scale space-time mortar variable that approximates the primary variable. Uniqueness, existence, and stability, as well as a priori error estimates for the spatial and temporal errors are established. A space-time non-overlapping domain decomposition method is developed that reduces the global problem to a space-time coarse-scale mortar interface problem. Each interface iteration involves solving in parallel space-time subdomain problems. The spectral properties of the interface operator and the convergence of the interface iteration are analyzed. Numerical experiments are provided that illustrate the theoretical results and the flexibility of the method for modeling problems with features that are localized in space and time.


title: Strict deformation quantization of abelian lattice gauge fields

abstract: This paper shows how to construct classical and quantum field C*-algebras modeling a $U(1)^n$-gauge theory in any dimension using a novel approach to lattice gauge theory, while simultaneously constructing a strict deformation quantization between the respective field algebras. The construction starts with quantization maps defined on operator systems (instead of C*-algebras) associated to the lattices, in a way that quantization commutes with all lattice refinements, therefore giving rise to a quantization map on the continuum (meaning ultraviolet and infrared) limit. Although working with operator systems at the finite level, in the continuum limit we obtain genuine C*-algebras. We also prove that the C*-algebras (classical and quantum) are invariant under time evolutions related to the electric part of abelian Yang--Mills. Our classical and quantum systems at the finite level are essentially the ones of [van Nuland and Stienstra, 2020], which admit completely general dynamics, and we briefly discuss ways to extend this powerful result to the continuum limit. We also briefly discuss reduction, and how the current set-up should be generalized to the non-abelian case.


title: Stochastic Multiplicative Weights Updates in Zero-Sum Games

abstract: We study agents competing against each other in a repeated network zero-sum game while applying the multiplicative weights update (MWU) algorithm with fixed learning rates. In our implementation, agents select their strategies probabilistically in each iteration and update their weights/strategies using the realized vector payoff of all strategies, i.e., stochastic MWU with full information. We show that the system results in an irreducible Markov chain where agent strategies diverge from the set of Nash equilibria. Further, we show that agents will play pure strategies with probability 1 in the limit.


title: The Quality of the 2020 Census: An Independent Assessment of Census Bureau Activities Critical to Data Quality

abstract: This report summarizes major findings from an independent evaluation of 2020 census operations. The American Statistical Association 2020 Census Quality Indicators Task Force selected the authors to conduct the evaluation using nonpublic operations data provided by the Census Bureau. The evaluation focused on the quality of state-level population counts released by the Census Bureau for congressional apportionment. The authors first partitioned the census enumeration process into five operation phases. Within each phase, one or more activities considered to be critical to census data quality were identified. Operational data from each activity were then analyzed to assess the risk of error, particularly as they related to similar activities in the 2010 census. Overall, the evaluation found that census operations relied on higher risk activities at a higher rate in 2020 than in 2010, suggesting that the risk of error may be higher in 2020 than in 2010. However, the available data were insufficient to determine whether the apportionment counts are of lower quality in 2020 than in 2010.


title: Learned Uncertainty Calibration for Visual Inertial Localization

abstract: The widely-used Extended Kalman Filter (EKF) provides a straightforward recipe to estimate the mean and covariance of the state given all past measurements in a causal and recursive fashion. For a wide variety of applications, the EKF is known to produce accurate estimates of the mean and typically inaccurate estimates of the covariance. For applications in visual inertial localization, we show that inaccuracies in the covariance estimates are \emph{systematic}, i.e. it is possible to learn a nonlinear map from the empirical ground truth to the estimated one. This is demonstrated on both a standard EKF in simulation and a Visual Inertial Odometry system on real-world data.


title: Polarization selectivity of aloof-beam electron energy-loss spectroscopy in one-dimensional ZnO nanorods

abstract: Orientation dependent electronic properties of wurtzite zinc oxide nanorods are characterized by aloof beam electron energy-loss spectroscopy (EELS) carried out in a scanning transmission electron microscope (STEM). The two key crystal orientation differentiating transitions specific to the in-plane (13.0 eV) and out-of-plane (11.2 eV) directions with respect to the wurtzite structure are examined by first principles density-functional theory calculations. We note some degree of orientation dependence at the onset of direct band gap transition near 3.4 eV. We demonstrate that good polarization selectivity can be achieved by placing the electron probe at different locations around the specimen with increasing impact parameter while keeping the beam-specimen orientation fixed. The observed results are qualitatively elucidated in terms of the perpendicular electric fields generated by the fast electron (60 kV) used in the microscope. The fact that good polarization selectivity can be achieved by aloof beam EELS without the requirement of sample reorientation is an attractive aspect from the characterization method point of view in the STEM-EELS community.


title: Fast barrier-free switching in synthetic antiferromagnets

abstract: We analytically solve the Landau-Lifshitz equations for the collective magnetization dynamics in a synthetic antiferromagnet (SAF) nanoparticle and uncover a regime of barrier-free switching under a short small-amplitude magnetic field pulse applied perpendicular to the SAF plane. We give examples of specific implementations for forming such low-power and ultra-fast switching pulses. For fully optical, resonant, barrier-free SAF switching we estimate the power per write operation to be $ \sim 100 $ pJ, 10-100 times smaller than for conventional quasi-static rotation, which should be attractive for memory applications.


title: Investigating the co-evolution of massive black holes in dual active galactic nuclei and their host galaxies via galaxy merger simulations

abstract: Major galaxy mergers can trigger nuclear activities and are responsible for high-luminosity quasi-stellar objects /active galactic nuclei (QSOs/AGNs). In certain circumstances, such mergers may cause dual active galactic nuclei (dAGN) phenomenon. This study investigates dAGN triggering and evolution of massive black holes (MBHs) during the merging processes using hydrodynamic code GADGET-2 to simulate several gas-rich major mergers at redshift $z=2$ and $3$, respectively. Results reveal that gas-rich major mergers can trigger significant nuclear activities after the second and third pericentric passages and the formation of dAGN with significant time duration ($\sim 10 - 390$ Myr). During the merging processes, galactic bulge evolves with time because of the rapid star formation in each (or both) galactic centers and initial mixing of stars in galactic disks due to violent relaxation. MBHs grow substantially due to accretion and finally merge into a bigger black hole. The growth of galactic bulges and corresponding increases of its velocity dispersions predate the growth of MBHs in the dAGN stages. The MBHs in these stages deviate below the relation between MBH mass and bulge mass (or velocity dispersion), and they revert to the relation after the final mergers due to the significant accretion that occurs mostly at a separation less than a few kpc. Then, the two MBHs merge with each other.


title: S2 Reducer: High-Performance Sparse Communication to Accelerate Distributed Deep Learning

abstract: Distributed stochastic gradient descent (SGD) approach has been widely used in large-scale deep learning, and the gradient collective method is vital to ensure the training scalability of the distributed deep learning system. Collective communication such as AllReduce has been widely adopted for the distributed SGD process to reduce the communication time. However, AllReduce incurs large bandwidth resources while most gradients are sparse in many cases since many gradient values are zeros and should be efficiently compressed for bandwidth saving. To reduce the sparse gradient communication overhead, we propose Sparse-Sketch Reducer (S2 Reducer), a novel sketch-based sparse gradient aggregation method with convergence guarantees. S2 Reducer reduces the communication cost by only compressing the non-zero gradients with count-sketch and bitmap, and enables the efficient AllReduce operators for parallel SGD training. We perform extensive evaluation against four state-of-the-art methods over five training models. Our results show that S2 Reducer converges to the same accuracy, reduces 81\% sparse communication overhead, and achieves 1.8$ \times $ speedup compared to state-of-the-art approaches.


title: Short-and-Sparse Deconvolution Via Rank-One Constrained Optimization (ROCO)

abstract: Short-and-sparse deconvolution (SaSD) aims to recover a short kernel and a long and sparse signal from their convolution. In the literature, formulations of blind deconvolution is either a convex programming via a matrix lifting of convolution, or a bilinear Lasso. Optimization solvers are typically based on bilinear factorizations. In this paper, we formulate SaSD as a non-convex optimization with a rank-one matrix constraint, hence referred to as Rank-One Constrained Optimization (ROCO). The solver is based on alternating direction method of multipliers (ADMM). It operates on the full rank-one matrix rather than bilinear factorizations. Closed form updates are derived for the efficiency of ADMM. Simulations include both synthetic data and real images. Results show substantial improvements in recovery accuracy (at least 19dB in PSNR for real images) and comparable runtime compared with benchmark algorithms based on bilinear factorization.


title: Lossy compression of statistical data using quantum annealer

abstract: We present a new lossy compression algorithm for statistical floating-point data through a representation learning with binary variables. The algorithm finds a set of basis vectors and their binary coefficients that precisely reconstruct the original data. The optimization for the basis vectors is performed classically, while binary coefficients are retrieved through both simulated and quantum annealing for comparison. A bias correction procedure is also presented to estimate and eliminate the error and bias introduced from the inexact reconstruction of the lossy compression for statistical data analyses. The compression algorithm is demonstrated on two different datasets of lattice quantum chromodynamics simulations. The results obtained using simulated annealing show 3.5 times better compression performance than the algorithms based on a neural-network autoencoder and principal component analysis. Calculations using quantum annealing also show promising results, but performance is limited by the integrated control error of the quantum processing unit, which yields large uncertainties in the biases and coupling parameters. Hardware comparison is further studied between the previous generation D-Wave 2000Q and the current D-Wave Advantage system. Our study shows that the Advantage system is more likely to obtain low-energy solutions for the problems than the 2000Q.


abstract: Let $\mathbb{F}_q$ be the finite field of order $q$, and $\mathbb P^1(\mathbb{F}_q) = \mathbb F_q\cup \{\infty\}$. Write $(x+\sqrt y)^m$ as $N(x,y)+D(x,y)\sqrt{y}$. For $m\in\mathbb N$ and $a \in \mathbb{F}_q$, the R\'edei function $R_{m,a}\colon \mathbb P^1(\mathbb F_q) \to \mathbb P^1(\mathbb F_q)$ is defined by $N(x,a)/D(x,a)$ if $D(x,a)\neq 0$ and $x\neq\infty$, and $\infty$, otherwise. In this paper we give a complete characterization of all pairs $(m,n)\in\mathbb N^2$ such that the R\'edei permutations $R_{m,a}$ and $R_{n,b}$ have the same cycle structure when $a$ and $b$ have the same quadratic character and $q$ is odd. We explore some relationships between such pairs $(m,n)$, and provide explicit families of R\'edei permutations with the same cycle structure. When a R\'edei permutation has a unique cycle structure that is not shared by any other R\'edei permutation, we call it isolated. We show that the only isolated R\'edei permutations are the isolated R\'edei involutions. Moreover, all our results can be transferred to bijections of the form $mx$ and $x^m$ on certain domains.


title: Late reverberation suppression using U-nets

abstract: In real-world settings, speech signals are almost always affected by reverberation produced by the working environment; these corrupted signals need to be \emph{dereverberated} prior to performing, e.g., speech recognition, speech-to-text conversion, compression, or general audio enhancement. In this paper, we propose a supervised dereverberation technique using \emph{U-nets with skip connections}, which are fully-convolutional encoder-decoder networks with layers arranged in the form of an "U" and connections that "skip" some layers. Building on this architecture, we address speech dereverberation through the lens of Late Reverberation Suppression (LS). Via experiments on synthetic and real-world data with different noise levels and reverberation settings, we show that our proposed method termed "LS U-net" improves quality, intelligibility and other performance metrics compared to the original U-net method and it is on par with the state-of-the-art GAN-based approaches.


title: Analyzing the Intrinsic Magnetic Field in the Galactic Center Radio Arc

abstract: The Radio Arc is a system of organized non-thermal filaments (NTFs) located within the Galactic Center (GC) region of the Milky Way. Recent observations of the Radio Arc NTFs revealed a magnetic field which alternates between being parallel and rotated with respect to the orientation of the filaments. This pattern is in stark contrast to the predominantly parallel magnetic field orientations observed in other GC NTFs. To help elucidate the origin of this pattern, we analyze spectro-polarimetric data of the Radio Arc NTFs using an Australian Telescope Compact Array data set covering the continuous frequency range from $\sim$4 to 11 GHz at a spectral resolution of 2 MHz. We fit depolarization models to the spectral polarization data to characterize Faraday effects along the line-of-sight. We assess whether structures local to the Radio Arc NTFs may contribute to the unusual magnetic field orientation. External Faraday effects are identified as the most likely origin of the rotation observed for the Radio Arc NTFs; however, internal Faraday effects are also found to be likely in regions of parallel magnetic field. The increased likelihood of internal Faraday effects in parallel magnetic field regions may be attributed to the effects of structures local to the GC. One such structure could be the Radio Shell local to the Radio Arc NTFs. Future studies are needed to determine whether this alternating magnetic field pattern is present in other multi-stranded NTFs, or is a unique property resulting from the complex interstellar region local to the Radio Arc NTFs.


title: Carrier transport and performance limit of semi-transparent photovoltaics: CuIn$_{1-x}$Ga$_x$Se$_2$ as a case study

abstract: Semi-transparent photovoltaic devices for building integrated applications have the potential to provide simultaneous power generation and natural light penetration. CuIn$_{1-x}$Ga$_x$Se$_2$ (CIGS) has been established as a mature technology for thin-film photovoltaics, however, its potential for Semi-Transparent Photovoltaics (STPV) is yet to be explored. In this paper, we present its carrier transport physics explaining the trend seen in recently published experiments. STPV requires deposition of films of only a few hundred nanometers to make them transparent and manifests several unique properties compared to a conventional thin-film solar cell. Our analysis shows that the short-circuit current, Jsc is dominated by carriers generated in the depletion region, making it nearly independent of bulk and back-surface recombination. The bulk recombination, which limits the open-circuit voltage Voc, appears to be higher than usual attributable to numerous grain boundaries. When the absorber layer is reduced below 500 nm, grain size reduces resulting in more grain boundaries and higher resistance. This produces an inverse relationship between series resistance and absorber thickness. We also present a thickness-dependent model of shunt resistance showing its impact in these ultra-thin devices. For various scenarios of bulk and interface recombinations, shunt and series resistances, AVT and composition of CuIn$_{1-x}$Ga$_x$Se$_2$, we project the efficiency limit which - for most practical cases - is found to be $\leq$10% for AVT $\geq$25%.


title: Drift and Matrix coefficients for discrete group extensions of countable Markov shifts

abstract: There has been much interest in generalizing Kesten's criterion for amenability in terms of a random walk to other contexts, such as determining amenability of a deck covering group by the bottom of the spectrum of the Laplacian or entropy of the geodesic flow. One outcome of this work is to generalise the results to so-called discrete group extensions of countable Markov shifts that satisfy a strong positive recurrence hypothesis. The other outcome is to further develop the language of unitary representation theory in this problem, and to bring some of the machinery developed by Coulon-Dougall-Schapira-Tapie [Twisted Patterson-Sullivan measures and applications to amenability and coverings, arXiv:1809.10881, 2018] to the countable Markov shift setting. In particular we recast the problem of determining a drop in Gurevi\v{c} pressure in terms of eventual almost sure decay for matrix coefficients, and explain that a so-called twisted measure "finds points with the worst decay." We are also able locate the results of Dougall-Sharp [Anosov flows, growth rates on covers and group extensions of subshifts, {\em Inventiones Mathematicae}, 223, 445-483, 2021] within this framework.


title: NaRLE: Natural Language Models using Reinforcement Learning with Emotion Feedback

abstract: Current research in dialogue systems is focused on conversational assistants working on short conversations in either task-oriented or open domain settings. In this paper, we focus on improving task-based conversational assistants online, primarily those working on document-type conversations (e.g., emails) whose contents may or may not be completely related to the assistant's task. We propose "NARLE" a deep reinforcement learning (RL) framework for improving the natural language understanding (NLU) component of dialogue systems online without the need to collect human labels for customer data. The proposed solution associates user emotion with the assistant's action and uses that to improve NLU models using policy gradients. For two intent classification problems, we empirically show that using reinforcement learning to fine tune the pre-trained supervised learning models improves performance up to 43%. Furthermore, we demonstrate the robustness of the method to partial and noisy implicit feedback.


title: Bar-driven Leading Spiral Arms in a Counter-rotating Dark Matter Halo

abstract: An overwhelming majority of galactic spiral arms trail with respect to the rotation of the galaxy, though a small sample of leading spiral arms has been observed. The formation of these leading spirals is not well understood. Here we show, using collisionless $N$-body simulations, that a barred disc galaxy in a retrograde dark matter halo can produce long-lived ($\sim3$ Gyr) leading spiral arms. Due to the strong resonant coupling of the disc to the halo, the bar slows rapidly and spiral perturbations are forced ahead of the bar. We predict that such a system, if observed, will also host a dark matter wake oriented perpendicular to the stellar bar. More generally, we propose that any mechanism that rapidly decelerates the stellar bar will allow leading spiral arms to flourish.


title: Online Application Guidance for Heterogeneous Memory Systems

abstract: Many high end and next generation computing systems to incorporated alternative memory technologies to meet performance goals. Since these technologies present distinct advantages and tradeoffs compared to conventional DDR* SDRAM, such as higher bandwidth with lower capacity or vice versa, they are typically packaged alongside conventional SDRAM in a heterogeneous memory architecture. To utilize the different types of memory efficiently, new data management strategies are needed to match application usage to the best available memory technology. However, current proposals for managing heterogeneous memories are limited because they either: 1) do not consider high-level application behavior when assigning data to different types of memory, or 2) require separate program execution (with a representative input) to collect information about how the application uses memory resources. This work presents a toolset for addressing the limitations of existing approaches for managing complex memories. It extends the application runtime layer with automated monitoring and management routines that assign application data to the best tier of memory based on previous usage, without any need for source code modification or a separate profiling run. It evaluates this approach on a state-of-the-art server platform with both conventional DDR4 SDRAM and non-volatile Intel Optane DC memory, using both memory-intensive high performance computing (HPC) applications as well as standard benchmarks. Overall, the results show that this approach improves program performance significantly compared to a standard unguided approach across a variety of workloads and system configurations. Additionally, we show that this approach achieves similar performance as a comparable offline profiling-based approach after a short startup period, without requiring separate program execution or offline analysis steps.


title: Detection of blue whale vocalisations using a temporal-domain convolutional neural network

abstract: We present a framework for detecting blue whale vocalisations from acoustic submarine recordings. The proposed methodology comprises three stages: i) a preprocessing step where the audio recordings are conditioned through normalisation, filtering, and denoising; ii) a label-propagation mechanism to ensure the consistency of the annotations of the whale vocalisations, and iii) a convolutional neural network that receives audio samples. Based on 34 real-world submarine recordings (28 for training and 6 for testing) we obtained promising performance indicators including an Accuracy of 85.4\% and a Recall of 93.5\%. Furthermore, even for the cases where our detector did not match the ground-truth labels, a visual inspection validates the ability of our approach to detect possible parts of whale calls unlabelled as such due to not being complete calls.


title: Operation-Adversarial Scenario Generation

abstract: This paper proposes a modified conditional generative adversarial network (cGAN) model to generate net load scenarios for power systems that are statistically credible, conditioned by given labels (e.g., seasons), and, at the same time, "stressful" to the system operations and dispatch decisions. The measure of stress used in this paper is based on the operating cost increases due to net load changes. The proposed operation-adversarial cGAN (OA-cGAN) internalizes a DC optimal power flow model and seeks to maximize the operating cost and achieve a worst-case data generation. The training and testing stages employed in the proposed OA-cGAN use historical day-ahead net load forecast errors and has been implemented for the realistic NYISO 11-zone system. Our numerical experiments demonstrate that the generated operation-adversarial forecast errors lead to more cost-effective and reliable dispatch decisions.


title: Inference and De-Noising of Non-Gaussian Particle Distribution Functions: A Generative Modeling Approach

abstract: The particle-in-cell numerical method of plasma physics balances a trade-off between computational cost and intrinsic noise. Inference on data produced by these simulations generally consists of binning the data to recover the particle distribution function, from which physical processes may be investigated. In addition to containing noise, the distribution function is temporally dynamic and can be non-gaussian and multi-modal, making the task of modeling it difficult. Here we demonstrate the use of normalizing flows to learn a smooth, tractable approximation to the noisy particle distribution function. We demonstrate that the resulting data driven likelihood conserves relevant physics and may be extended to encapsulate the temporal evolution of the distribution function.


title: Weighted inequalities for discrete iterated kernel operators

abstract: We develop a new method that enables us to solve the open problem of characterizing discrete inequalities for kernel operators involving suprema. More precisely, we establish necessary and sufficient conditions under which there exists a positive constant $C$ such that \begin{equation*} \Bigg(\sum_{n\in\mathbb{Z}}\Bigg(\sum_{i=-\infty}^n {U}(i,n) a_i\Bigg)^{q} {w}_n\Bigg)^{\frac{1}{q}} \le C \Bigg(\sum_{n\in\mathbb{Z}}a_n^p {v}_n\Bigg)^{\frac{1}{p}} \end{equation*} holds for every sequence of nonnegative numbers $\{a_n\}_{n\in\mathbb{Z}}$ where $U$ is a kernel satisfying certain regularity condition, $0 < p,q \leq \infty$ and $\{v_n\}_{n\in\mathbb{Z}}$ and $\{w_n\}_{n\in\mathbb{Z}}$ are fixed weight sequences. We do the same for the inequality \begin{equation*} \Bigg( \sum_{n\in\mathbb{Z}} w_n \Big[ \sup_{-\infty<i\le n} U(i,n) \sum_{j=-\infty}^{i} a_j \Big]^q \Bigg)^{\frac{1}{q}} \le C \Bigg( \sum_{n\in\mathbb{Z}} a_n^p v_n \Bigg)^{\frac{1}{p}}. \end{equation*} We characterize these inequalities by conditions of both discrete and continuous nature.


title: Urns & Tubes

abstract: Urn models play an important role to express various basic ideas in probability theory. Here we extend this urn model with tubes. An urn contains coloured balls, which can be drawn with probabilities proportional to the numbers of balls of each colour. For each colour a tube is assumed. These tubes have different sizes (lengths). The idea is that after drawing a ball from the urn it is dropped in the urn of the corresponding colour. We consider two associated probability distributions. The first-full distribution on colours gives for each colour the probability that the corresponding tube is full first, before any of the other tubes. The negative distribution on natural numbers captures for a number k the probability that all tubes are full for the first time after k draws. This paper uses multisets to systematically describe these first-full and negative distributions in the urns & tubes setting, in fully multivariate form, for all three standard drawing modes (multinomial, hypergeometric, and Polya).


abstract: Autoregressive (AR) time series models are widely used in parametric spectral estimation (SE), where the power spectral density (PSD) of the time series is approximated by that of the \emph{best-fit} AR model, which is available in closed form. Since AR parameters are usually found via maximum-likelihood, least squares or the method of moments, AR-based SE fails to account for the uncertainty of the approximate PSD, and thus only yields point estimates. We propose to handle the uncertainty related to the AR approximation by finding the full posterior distribution of the AR parameters to then propagate this uncertainty to the PSD approximation by \emph{integrating out the AR parameters}; we implement this concept by assuming two different priors over the model noise. Through practical experiments, we show that the proposed Bayesian autoregressive spectral estimation (BASE) provides point estimates that follow closely those of standard autoregressive spectral estimation (ASE), while also providing error bars. BASE is validated against ASE and the Periodogram on both synthetic and real-world signals.


title: Double Inozemtsev Limits of the Quantum DELL System

abstract: In this letter we study various Inozemtsev-type limits of the quantum double elliptic (DELL) system when both elliptic parameters are sent to zero at different rates, while the coupling constant is sent to infinity, such that a certain combination of the three parameters is kept fixed. We find a regime in which such double Inozemtsev limit of DELL produces the elliptic Ruijsenaars-Schneider (eRS) Hamiltonians albeit in an unconventional normalization. We discuss other double scaling limits and anisotropic scaling of coordinates and momenta. In addition, we provide a formal expression for the eigenvalues of the eRS Hamiltonians solely in terms of their eigenfunctions.


title: Exploring usability of Reddit in data science and knowledge processing

abstract: This contribution argues that Reddit, as a massive, categorized, open-access dataset, is a useful data source, for "almost any topic". Hence, it can be used in data science, e.g. for knowledge exploration. This statement is backed-up with presented analysis, based on 180 manually annotated papers, related to Reddit itself, and data acquired from popular databases of scientific papers. Finally, an open source tool is introduced, which provides an easy access to Reddit resources, and an exploratory data analysis of how Reddit covers selected topics. These functions can be used as a prelude analysis to a broader exploration of Reddit's applicability.


title: Label differential privacy via clustering

abstract: We present new mechanisms for \emph{label differential privacy}, a relaxation of differentially private machine learning that only protects the privacy of the labels in the training set. Our mechanisms cluster the examples in the training set using their (non-private) feature vectors, randomly re-sample each label from examples in the same cluster, and output a training set with noisy labels as well as a modified version of the true loss function. We prove that when the clusters are both large and high-quality, the model that minimizes the modified loss on the noisy training set converges to small excess risk at a rate that is comparable to the rate for non-private learning. We describe both a centralized mechanism in which the entire training set is stored by a trusted curator, and a distributed mechanism where each user stores a single labeled example and replaces her label with the label of a randomly selected user from the same cluster. We also describe a learning problem in which large clusters are necessary to achieve both strong privacy and either good precision or good recall. Our experiments show that randomizing the labels within each cluster significantly improves the privacy vs. accuracy trade-off compared to applying uniform randomized response to the labels, and also compared to learning a model via DP-SGD.