title: From H\"older Triangles to the Whole Plane

abstract: We show how to determine whether two given real polynomial functions of a single variable are Lipschitz equivalent by comparing the values and also the multiplicities of the given polynomial functions at their critical points. Then we show how to reduce the problem of ${\cal R}$-semialgebraic Lipschitz classification of $\beta$-quasihomogeneous polynomials of two real variables to the problem of Lipschitz classification of real polynomial functions of a single variable, under some fairly general conditions.


title: An Ode to an ODE

abstract: We present a new paradigm for Neural ODE algorithms, called ODEtoODE, where time-dependent parameters of the main flow evolve according to a matrix flow on the orthogonal group O(d). This nested system of two flows, where the parameter-flow is constrained to lie on the compact manifold, provides stability and effectiveness of training and provably solves the gradient vanishing-explosion problem which is intrinsically related to training deep neural network architectures such as Neural ODEs. Consequently, it leads to better downstream models, as we show on the example of training reinforcement learning policies with evolution strategies, and in the supervised learning setting, by comparing with previous SOTA baselines. We provide strong convergence results for our proposed mechanism that are independent of the depth of the network, supporting our empirical studies. Our results show an intriguing connection between the theory of deep neural networks and the field of matrix flows on compact manifolds.


title: Deterministic homogenization under optimal moment assumptions for fast-slow systems. Part 1

abstract: We consider deterministic homogenization (convergence to a stochastic differential equation) for multiscale systems of the form \[ x_{k+1} = x_k + n^{-1} a_n(x_k,y_k) + n^{-1/2} b_n(x_k,y_k), \quad y_{k+1} = T_n y_k, \] where the fast dynamics is given by a family $T_n$ of nonuniformly expanding maps. Part 1 builds on our recent work on martingale approximations for families of nonuniformly expanding maps. We prove an iterated weak invariance principle and establish optimal iterated moment bounds for such maps. (The iterated moment bounds are new even for a fixed nonuniformly expanding map T.) The homogenization results are a consequence of this together with parallel developments on rough path theory in Part 2 by Chevyrev, Friz, Korepanov, Melbourne & Zhang.


title: An Adaptive MMC Synchronous Stability Control Method Based on Local PMU measurements

abstract: Reducing the current is a common method to ensure the synchronous stability of a modular multilevel converter (MMC) when there is a short-circuit fault at its AC side. However, the uncertainty of the fault location of the AC system leads to a significant difference in the maximum allowable stable operating current during the fault. This paper proposes an adaptive MMC fault-current control method using local phasor measurement unit (PMU) measurements. Based on the estimated Thevenin equivalent (TE) parameters of the system, the current can be directly calculated to ensure the maximum output power of the MMC during the fault. This control method does not rely on off-line simulation and adapts itself to various fault conditions. The effective measurements are firstly selected by the voltage threshold and parameter constraints, which allow us to handle the error due to the change on the system-side. The proposed TE estimation method can fast track the change of the system impedance without depending on the initial value and can deal with the TE potential changes after a large disturbance. The simulation shows that the TE estimation can accurately track the TE parameters after the fault, and the current control instruction during an MMC fault can ensure the maximum output power of the MMC.


title: Capturing Video Frame Rate Variations via Entropic Differencing

abstract: High frame rate videos are increasingly getting popular in recent years, driven by the strong requirements of the entertainment and streaming industries to provide high quality of experiences to consumers. To achieve the best trade-offs between the bandwidth requirements and video quality in terms of frame rate adaptation, it is imperative to understand the effects of frame rate on video quality. In this direction, we devise a novel statistical entropic differencing method based on a Generalized Gaussian Distribution model expressed in the spatial and temporal band-pass domains, which measures the difference in quality between reference and distorted videos. The proposed design is highly generalizable and can be employed when the reference and distorted sequences have different frame rates. Our proposed model correlates very well with subjective scores in the recently proposed LIVE-YT-HFR database and achieves state of the art performance when compared with existing methodologies.


title: Parity-based, bias-free optical quantum random number generation with min-entropy estimation

abstract: We describe the generation of sequences of random bits from the parity of photon counts produced by polarization measurements on a polarization-entangled state. The resulting sequences are bias free, pass the applicable tests in the NIST battery of statistical randomness tests, and are shown to be Borel normal, without the need for experimental calibration stages or postprocessing of the output. Because the photon counts are produced in the course of a measurement of the violation of the Clauser-Horne-Shimony-Holt inequality, we are able to concurrently verify the nonclassical nature of the photon statistics and estimate a lower bound on the min-entropy of the bit-generating source. The rate of bit production in our experiment is around 13 bits/s.


title: A note on Almgren-Chriss optimal execution problem with geometric Brownian motion

abstract: We solve explicitly the Almgren-Chriss optimal liquidation problem where the stock price process follows a geometric Brownian motion. Our technique is to work in terms of cash and to use functional analysis tools. We show that this framework extends readily to the case of a stochastic drift for the price process and the liquidation of a portfolio.


title: Velocity limits in the thermonuclear supernova ejection scenario for hypervelocity stars and the origin of US 708

abstract: Hypervelocity stars (HVS) are a class of stars moving at high enough velocities to be gravitationally unbound from the Galaxy. Ejection from a close binary system in which one of the components undergoes a thermonuclear supernova (SN) has emerged as a promising candidate production mechanism for the least massive specimens of this class. This study presents a thorough theoretical analysis of candidate progenitor systems of thermonuclear SNe in the single degenerate helium donor scenario in the relevant parameter space leading to the ejection of HVS. The primary goal is investigation of the, previously unclear, characteristics of the velocity spectra of the ejected component. Presented are the results of 390 binary model sequences computed with the MESA framework, investigating the evolution of supernova progenitors composed of a helium-rich hot subdwarf and a accreting white dwarf. Results are then correlated with an idealized kinematic analysis of the observed object US 708. It is seen that the ejection velocity spectra reach a maximum in the range $0.19~M_\odot < M_{HVS} < 0.25~M_\odot$. Depending on the local Galactic potential, all donors below $0.4~\text{M}_\odot$ are expected to become HVS. This channel is able to account for runaway velocities up to $\sim1150~\text{km s}^{-1}$ with a Chandrasekhar mass accretor, exceeding $1200~\text{km s}^{-1}$ if super-Chandrasekhar mass detonations are taken into account. It is found that the previously assumed mass of $0.3~M_\odot$ for US 708, combined with more recently obtained proper motions, favor a sub-Chandrasekhar mass explosion with a terminal WD mass between $1.1~M_\odot$ and $1.2~M_\odot$. The presence of clear ejection velocity maxima provides constraints on the terminal state of a supernova progenitor. It is possible to discern certain types of explosion mechanisms from the inferred ejection velocities alone.


title: Frequently recurrent operators

abstract: Motivated by a recent investigation of Costakis et al. on the notion of recurrence in linear dynamics, we study various stronger forms of recurrence for linear operators, in particular that of frequent recurrence. We study, among other things, the relationship between a type of recurrence and the corresponding notion of hypercyclicity, the influence of power boundedness, and the interplay between recurrence and spectral properties. We obtain, in particular, Ansari- and L\'eon-M\"uller-type theorems for $\mathcal{F}$-recurrence under very weak assumptions on the Furstenberg family $\mathcal{F}$. This allows us, as a by-product, to deduce Ansari- and L\'eon-M\"uller-type theorems for $\mathcal{F}$-hypercyclicity.


title: Absence of renormalization group pathologies in some critical Dyson-Ising ferromagnets

abstract: The Dyson-Ising ferromagnet is a one-dimensional Ising model with a power law interaction. When the power is between -1 and -2, the model has a phase transition. Van Enter and Le Ny proved that at sufficiently low temperatures the decimation renormalization group transformation is not defined in the sense that the renormalized measure is not a Gibbs measure. We consider a modified model in which the nearest neighbor couplings are much larger than the other couplings. For a family of Hamiltonians which includes critical cases, we prove that the first step of the renormalization group transformation can be rigorously defined for majority rule and decimation.


title: Learning Minimax Estimators via Online Learning

abstract: We consider the problem of designing minimax estimators for estimating the parameters of a probability distribution. Unlike classical approaches such as the MLE and minimum distance estimators, we consider an algorithmic approach for constructing such estimators. We view the problem of designing minimax estimators as finding a mixed strategy Nash equilibrium of a zero-sum game. By leveraging recent results in online learning with non-convex losses, we provide a general algorithm for finding a mixed-strategy Nash equilibrium of general non-convex non-concave zero-sum games. Our algorithm requires access to two subroutines: (a) one which outputs a Bayes estimator corresponding to a given prior probability distribution, and (b) one which computes the worst-case risk of any given estimator. Given access to these two subroutines, we show that our algorithm outputs both a minimax estimator and a least favorable prior. To demonstrate the power of this approach, we use it to construct provably minimax estimators for classical problems such as estimation in the finite Gaussian sequence model, and linear regression.


title: Band-limited Soft Actor Critic Model

abstract: Soft Actor Critic (SAC) algorithms show remarkable performance in complex simulated environments. A key element of SAC networks is entropy regularization, which prevents the SAC actor from optimizing against fine grained features, oftentimes transient, of the state-action value function. This results in better sample efficiency during early training. We take this idea one step further by artificially bandlimiting the target critic spatial resolution through the addition of a convolutional filter. We derive the closed form solution in the linear case and show that bandlimiting reduces the interdependency between the low and high frequency components of the state-action value approximation, allowing the critic to learn faster. In experiments, the bandlimited SAC outperformed the classic twin-critic SAC in a number of Gym environments, and displayed more stability in returns. We derive novel insights about SAC by adding a stochastic noise disturbance, a technique that is increasingly being used to learn robust policies that transfer well to the real world counterparts.


title: Online Kernel based Generative Adversarial Networks

abstract: One of the major breakthroughs in deep learning over the past five years has been the Generative Adversarial Network (GAN), a neural network-based generative model which aims to mimic some underlying distribution given a dataset of samples. In contrast to many supervised problems, where one tries to minimize a simple objective function of the parameters, GAN training is formulated as a min-max problem over a pair of network parameters. While empirically GANs have shown impressive success in several domains, researchers have been puzzled by unusual training behavior, including cycling so-called mode collapse. In this paper, we begin by providing a quantitative method to explore some of the challenges in GAN training, and we show empirically how this relates fundamentally to the parametric nature of the discriminator network. We propose a novel approach that resolves many of these issues by relying on a kernel-based non-parametric discriminator that is highly amenable to online training---we call this the Online Kernel-based Generative Adversarial Networks (OKGAN). We show empirically that OKGANs mitigate a number of training issues, including mode collapse and cycling, and are much more amenable to theoretical guarantees. OKGANs empirically perform dramatically better, with respect to reverse KL-divergence, than other GAN formulations on synthetic data; on classical vision datasets such as MNIST, SVHN, and CelebA, show comparable performance.


title: Local Fourier analysis of multigrid for hybridized and embedded discontinuous Galerkin methods

abstract: In this paper we present a geometric multigrid method with Jacobi and Vanka relaxation for hybridized and embedded discontinuous Galerkin discretizations of the Laplacian. We present a local Fourier analysis (LFA) of the two-grid error-propagation operator and show that the multigrid method applied to an embedded discontinuous Galerkin (EDG) discretization is almost as efficient as when applied to a continuous Galerkin discretization. We furthermore show that multigrid applied to an EDG discretization outperforms multigrid applied to a hybridized discontinuous Galerkin (HDG) discretization. Numerical examples verify our LFA predictions.


title: ADDENDUM Details of the Derivation of the Probability of Coverage for the Relaying Scheme (Section IV in the paper: "A Poisson Line Process based Framework for Determining the Needed RSU Density and Relaying Hops in Vehicular Networks'')

abstract: This paper develops a framework to study multi-hop relaying in a vehicular network consisting of vehicles and Road Side Units (RSUs), and the effect of this relaying on the network coverage and the communication delay. We use a stochastic geometry model that consists of a combination of Poisson Line Process (PLP) and 1D Poisson Point Process (PPP) to reliably characterize the vehicular network layout and the locations of the vehicles and the RSUs. Using this model, we analyze the effect of the different network parameters on the coverage provided by the RSUs to the vehicles. Then, we investigate how the uncovered vehicles can receive their intended packets by relaying them through multiple hops that form connected paths to the RSUs. We also analyze the delay introduced to packet delivery due to multi-hop relaying. Namely, we present results that illustrate the coverage gains achieved through multi-hop relaying and the delays induced. Such results could be used by network planners and operators to decide on the different configurations and operational parameters of the vehicular network to suit particular scenarios and objectives.


title: Channel Estimation for One-Bit Multiuser Massive MIMO Using Conditional GAN

abstract: Channel estimation is a challenging task, especially in a massive multiple-input multiple-output (MIMO) system with one-bit analog-to-digital converters (ADC). Traditional deep learning (DL) methods, that learn the mapping from inputs to real channels, have significant difficulties in estimating accurate channels because their loss functions are not well designed and investigated. In this paper, a conditional generative adversarial networks (cGAN) is developed to predict more realistic channels by adversarially training two DL networks. cGANs not only learn the mapping from quantized observations to real channels but also learn an adaptive loss function to correctly train the networks. Numerical results show that the proposed cGAN based approach outperforms existing DL methods and achieves high robustness in massive MIMO systems.


title: BEV-Seg: Bird's Eye View Semantic Segmentation Using Geometry and Semantic Point Cloud

abstract: Bird's-eye-view (BEV) is a powerful and widely adopted representation for road scenes that captures surrounding objects and their spatial locations, along with overall context in the scene. In this work, we focus on bird's eye semantic segmentation, a task that predicts pixel-wise semantic segmentation in BEV from side RGB images. This task is made possible by simulators such as Carla, which allow for cheap data collection, arbitrary camera placements, and supervision in ways otherwise not possible in the real world. There are two main challenges to this task: the view transformation from side view to bird's eye view, as well as transfer learning to unseen domains. Existing work transforms between views through fully connected layers and transfer learns via GANs. This suffers from a lack of depth reasoning and performance degradation across domains. Our novel 2-staged perception pipeline explicitly predicts pixel depths and combines them with pixel semantics in an efficient manner, allowing the model to leverage depth information to infer objects' spatial locations in the BEV. In addition, we transfer learning by abstracting high-level geometric features and predicting an intermediate representation that is common across different domains. We publish a new dataset called BEVSEG-Carla and show that our approach improves state-of-the-art by 24% mIoU and performs well when transferred to a new domain.


title: Obliquity of an Earth-like planet from frequency modulation of its direct imaged lightcurve: mock analysis from general circulation model simulation

abstract: Direct-imaging techniques of exoplanets have made significant progress recently, and will eventually enable to monitor photometric and spectroscopic signals of earth-like habitable planets in the future. The presence of clouds, however, would remain as one of the most uncertain components in deciphering such direct-imaged signals of planets. We attempt to examine how the planetary obliquity produce different cloud patterns by performing a series of GCM (General Circulation Model) simulation runs using a set of parameters relevant for our Earth. Then we use the simulated photometric lightcurves to compute their frequency modulation due to the planetary spin-orbit coupling over an entire orbital period, and attempt to see to what extent one can estimate the obliquity of an Earth-twin. We find that it is possible to estimate the obliquity of an Earth-twin within the uncertainty of several degrees with a dedicated 4 m space telescope at 10 pc away from the system if the stellar flux is completely blocked. While our conclusion is based on several idealized assumptions, a frequency modulation of a directly-imaged earth-like planet offers a unique methodology to determine its obliquity.


title: Deep Implicit Coordination Graphs for Multi-agent Reinforcement Learning

abstract: Multi-agent reinforcement learning (MARL) requires coordination to efficiently solve certain tasks. Fully centralized control is often infeasible in such domains due to the size of joint action spaces. Coordination graph based formalization allows reasoning about the joint action based on the structure of interactions. However, they often require domain expertise in their design. This paper introduces the deep implicit coordination graph (DICG) architecture for such scenarios. DICG consists of a module for inferring the dynamic coordination graph structure which is then used by a graph neural network based module to learn to implicitly reason about the joint actions or values. DICG allows learning the tradeoff between full centralization and decentralization via standard actor-critic methods to significantly improve coordination for domains with large number of agents. We apply DICG to both centralized-training-centralized-execution and centralized-training-decentralized-execution regimes. We demonstrate that DICG solves the relative overgeneralization pathology in predatory-prey tasks as well as outperforms various MARL baselines on the challenging StarCraft II Multi-agent Challenge (SMAC) and traffic junction environments.


title: Two Simple Ways to Learn Individual Fairness Metrics from Data

abstract: Individual fairness is an intuitive definition of algorithmic fairness that addresses some of the drawbacks of group fairness. Despite its benefits, it depends on a task specific fair metric that encodes our intuition of what is fair and unfair for the ML task at hand, and the lack of a widely accepted fair metric for many ML tasks is the main barrier to broader adoption of individual fairness. In this paper, we present two simple ways to learn fair metrics from a variety of data types. We show empirically that fair training with the learned metrics leads to improved fairness on three machine learning tasks susceptible to gender and racial biases. We also provide theoretical guarantees on the statistical performance of both approaches.


title: Local Convolutions Cause an Implicit Bias towards High Frequency Adversarial Examples

abstract: Adversarial Attacks are still a significant challenge for neural networks. Recent work has shown that adversarial perturbations typically contain high-frequency features, but the root cause of this phenomenon remains unknown. Inspired by theoretical work on linear full-width convolutional models, we hypothesize that the local (i.e. bounded-width) convolutional operations commonly used in current neural networks are implicitly biased to learn high frequency features, and that this is one of the root causes of high frequency adversarial examples. To test this hypothesis, we analyzed the impact of different choices of linear and nonlinear architectures on the implicit bias of the learned features and the adversarial perturbations, in both spatial and frequency domains. We find that the high-frequency adversarial perturbations are critically dependent on the convolution operation because the spatially-limited nature of local convolutions induces an implicit bias towards high frequency features. The explanation for the latter involves the Fourier Uncertainty Principle: a spatially-limited (local in the space domain) filter cannot also be frequency-limited (local in the frequency domain). Furthermore, using larger convolution kernel sizes or avoiding convolutions (e.g. by using Vision Transformers architecture) significantly reduces this high frequency bias, but not the overall susceptibility to attacks. Looking forward, our work strongly suggests that understanding and controlling the implicit bias of architectures will be essential for achieving adversarial robustness.


title: Task-Agnostic Online Reinforcement Learning with an Infinite Mixture of Gaussian Processes

abstract: Continuously learning to solve unseen tasks with limited experience has been extensively pursued in meta-learning and continual learning, but with restricted assumptions such as accessible task distributions, independently and identically distributed tasks, and clear task delineations. However, real-world physical tasks frequently violate these assumptions, resulting in performance degradation. This paper proposes a continual online model-based reinforcement learning approach that does not require pre-training to solve task-agnostic problems with unknown task boundaries. We maintain a mixture of experts to handle nonstationarity, and represent each different type of dynamics with a Gaussian Process to efficiently leverage collected data and expressively model uncertainty. We propose a transition prior to account for the temporal dependencies in streaming data and update the mixture online via sequential variational inference. Our approach reliably handles the task distribution shift by generating new models for never-before-seen dynamics and reusing old models for previously seen dynamics. In experiments, our approach outperforms alternative methods in non-stationary tasks, including classic control with changing dynamics and decision making in different driving scenarios.


title: 2D Retrieval Frameworks for Hot Jupiter Phase Curves

abstract: Spectroscopic phase curves provide unique access to the three-dimensional properties of transiting exoplanet atmospheres. However, a modeling framework must be developed to deliver accurate inferences of atmospheric properties for these complex data sets. Here, we develop an approach to retrieve temperature structures and molecular abundances from phase curve spectra at any orbital phase. In the context of a representative hot Jupiter with a large day-night temperature contrast, we examine the biases in typical one-dimensional (1D) retrievals as a function of orbital phase/geometry, compared to two-dimensional (2D) models that appropriately capture the disk-integrated phase geometry. We guide our intuition by applying our new framework on a simulated HST+Spitzer phase curve data set in which the "truth" is known, followed by an application to the spectroscopic phase curve of the canonical hot Jupiter, WASP-43b. We also demonstrate the retrieval framework on simulated JWST phase curve observations. We apply our new geometric framework to a joint-fit of all spectroscopic phases, assuming longitudinal molecular abundance homogeneity, resulting in an a factor of 2 improvement in abundances precision when compared to individual phase constraints. With a 1D retrieval model on simulated HST+Spitzer data, we find strongly biased molecular abundances for CH$_4$ and CO$_2$ at most orbital phases. With 2D, the day and night profiles retrieved from WASP-43b remain consistent throughout the orbit. JWST retrievals show that a 2D model is strongly favored at all orbital phases. Based on our new 2D retrieval implementation, we provide recommendations on when 1D models are appropriate and when more complex phase geometries involving multiple TP profiles are required to obtain an unbiased view of tidally locked planetary atmospheres.


title: PCA-based Antenna Impedance Estimation in Rayleigh Fading Channels

abstract: Impedance matching between receive antenna and front-end significantly impacts channel capacity in wireless channels. To implement capacity-optimal matching, the receiver must know the antenna impedance. But oftentimes this impedance varies with loading conditions in the antenna near-field. To mitigate this variation, several authors have proposed antenna impedance estimation techniques. However, the optimality of these techniques remains unclear. In this paper, we consider antenna impedance estimation at MISO receivers over correlated Rayleigh fading channels. We derive in closed-form the optimal ML estimator under i.i.d. fading and then show it can be found via a scalar optimization in generally correlated fading channels. Numerical results suggest a computationally efficient, principal-components approach that estimates antenna impedance in real-time for all Rayleigh fading channels. Furthermore, ergodic capacity can be significantly boosted at originally poorly matched receivers with adaptive matching using our proposed approach.


title: Optimising Monotone Chance-Constrained Submodular Functions Using Evolutionary Multi-Objective Algorithms

abstract: Many real-world optimisation problems can be stated in terms of submodular functions. A lot of evolutionary multi-objective algorithms have recently been analyzed and applied to submodular problems with different types of constraints. We present a first runtime analysis of evolutionary multi-objective algorithms for chance-constrained submodular functions. Here, the constraint involves stochastic components and the constraint can only be violated with a small probability of alpha. We show that the GSEMO algorithm obtains the same worst case performance guarantees as recently analyzed greedy algorithms. Furthermore, we investigate the behavior of evolutionary multi-objective algorithms such as GSEMO and NSGA-II on different submodular chance constrained network problems. Our experimental results show that this leads to significant performance improvements compared to the greedy algorithm.


title: Vertex Partitions into an Independent Set and a Forest with Each Component Small

abstract: For each integer k >= 2, we determine a sharp bound on mad(G) such that V(G) can be partitioned into sets I and F_k, where I is an independent set and G[F_k] is a forest in which each component has at most k vertices. For each k we construct an infinite family of examples showing our result is best possible. Our results imply that every planar graph G of girth at least 9 (resp. 8, 7) has a partition of V(G) into an independent set I and a set F such that G[F] is a forest with each component of order at most 3 (resp. 4, 6). Hendrey, Norin, and Wood asked for the largest function g(a,b) such that if mad(G) < g(a,b) then V(G) has a partition into sets A and B such that mad(G[A]) < a and mad(G[B]) < b. They specifically asked for the value of g(1,b), i.e., the case when A is an independent set. Previously, the only values known were g(1,4/3) and g(1,2). We find g(1,b) whenever 4/3 < b < 2.


title: MALOnt: An Ontology for Malware Threat Intelligence

abstract: Malware threat intelligence uncovers deep information about malware, threat actors, and their tactics, Indicators of Compromise(IoC), and vulnerabilities in different platforms from scattered threat sources. This collective information can guide decision making in cyber defense applications utilized by security operation centers(SoCs). In this paper, we introduce an open-source malware ontology - MALOnt that allows the structured extraction of information and knowledge graph generation, especially for threat intelligence. The knowledge graph that uses MALOnt is instantiated from a corpus comprising hundreds of annotated malware threat reports. The knowledge graph enables the analysis, detection, classification, and attribution of cyber threats caused by malware. We also demonstrate the annotation process using MALOnt on exemplar threat intelligence reports. A work in progress, this research is part of a larger effort towards auto-generation of knowledge graphs (KGs)for gathering malware threat intelligence from heterogeneous online resources.


title: Exact Large Time Behavior of Spherically-Symmetric Plasmas

abstract: We consider the classical and relativistic Vlasov-Poisson systems with spherically-symmetric initial data and prove the optimal decay rates for all suitable $L^p$ norms of the charge density and electric field, as well as, the optimal growth rates for the largest particle position and momentum on the support of the distribution function. Though a previous work \cite{Horst} established upper bounds on the decay of the supremum of the charge density and electric field, we provide a slightly different proof, attain optimal rates, and extend this result to include all other norms. Additionally, we prove sharp lower bounds on each of the aforementioned quantities and establish the time-asymptotic behavior of all spatial and momentum characteristics. Finally, we investigate the limiting behavior of the spatial average of the particle distribution as $t \to \infty$. In particular, we show that it converges uniformly to a smooth, compactly-supported function that preserves the mass, angular momentum, and energy of the system and depends only upon limiting particle momenta.


title: Interlacing and scaling exponents for the geodesic watermelon in last passage percolation

abstract: In discrete planar last passage percolation (LPP), random values are assigned independently to each vertex in $\mathbb Z^2$, and each finite upright path in $\mathbb Z^2$ is ascribed the weight given by the sum of values of its vertices. The weight of a collection of disjoint paths is the sum of its members' weights. The notion of a geodesic, a maximum weight path between two vertices, has a natural generalization concerning several disjoint paths: a $k$-geodesic watermelon in $[1,n]^2\cap\mathbb Z^2$ is a collection of $k$ disjoint paths contained in this square that has maximum weight among all such collections. While the weights of such collections are known to be important objects, the maximizing paths have been largely unexplored beyond the $k=1$ case. For exactly solvable models, such as exponential and geometric LPP, it is well known that for $k=1$ the exponents that govern fluctuation in weight and transversal distance are $1/3$ and $2/3$; that is, typically, the weight of the geodesic on the route $(1,1) \to (n,n)$ fluctuates around a dominant linear growth of the form $\mu n$ by the order of $n^{1/3}$; and the maximum Euclidean distance of the geodesic from the diagonal has order $n^{2/3}$. Assuming a strong but local form of convexity and one-point moderate deviation bounds for the geodesic weight profile---which are available in all known exactly solvable models---we establish that, typically, the $k$-geodesic watermelon's weight falls below $\mu nk$ by order $k^{5/3}n^{1/3}$, and its transversal fluctuation is of order $k^{1/3}n^{2/3}$. Our arguments crucially rely on, and develop, a remarkable deterministic interlacing property that the watermelons admit. Our methods also yield sharp rigidity estimates for naturally associated point processes, which improve on estimates obtained via tools from the theory of determinantal point processes available in the integrable setting.


title: Bi-orthogonal harmonics for the decomposition of gravitational radiation I: angular modes, completeness, and the introduction of adjoint-spheroidal harmonics

abstract: The estimation of radiative modes is a central problem in gravitational wave theory, with essential applications in signal modeling and data analysis. This problem is complicated by most astrophysically relevant systems' not having modes that are analytically tractable. A ubiquitous workaround is to use not modes, but multipole moments defined by spin weighted spherical harmonics. However, spherical multipole moments are only related to the modes of systems without angular momentum. As a result, they can obscure the underlying physics of astrophysically relevant systems, such as binary black hole merger and ringdown. In such cases, spacetime angular momentum means that radiative modes are not spherical, but spheroidal in nature. Here, we work through various problems related to spheroidal harmonics. We show for the first time that spheroidal harmonics are not only capable of representing arbitrary gravitational wave signals, but that they also possess a kind of orthogonality not used before in general relativity theory. Along the way we present a new class of spin weighted harmonic functions dubbed ``adjoint-spheroidal" harmonics. These new functions may be used for the general estimation of spheroidal multipole moments via complete bi-orthogonal decomposition (in the angular domain). By construction, adjoint-spheroidal harmonics suppress mode-mixing effects known to plague spherical harmonic decomposition; as a result, they better approximate a system's true radiative modes. We discuss potential applications of these results. Lastly, we briefly comment on the challenges posed by the analogous problem with Teukolsky's radial functions


title: Demand Estimation from Sales Transaction Data -- Practical Extensions

abstract: In this paper we discuss practical limitations of the standard choice-based demand models used in the literature to estimate demand from sales transaction data. We present modifications and extensions of the models and discuss data preprocessing and solution techniques which are useful for practitioners dealing with sales transaction data. Among these, we present an algorithm to split sales transaction data observed under partial availability, we extend a popular Expectation Maximization (EM) algorithm for non-homogeneous product sets, and we develop two iterative optimization algorithms which can handle much of the extensions discussed in the paper.


title: Beyond a pale blue dot : how to search for possible bio-signatures on earth-like planets

abstract: The Earth viewed from outside the Solar system would be identified merely like a pale blue dot, as coined by Carl Sagan. In order to detect possible signatures of the presence of life on a second earth among several terrestrial planets discovered in a habit-able zone, one has to develop and establish a methodology to characterize the planet as something beyond a mere pale blue dot. We pay particular attention to the periodic change of the color of the dot according to the rotation of the planet. Because of the large-scale inhomogeneous distribution of the planetary surface, the reflected light of the dot comprises different color components corresponding to land, ocean, ice, and cloud that cover the surface of the planet. If we decompose the color of the dot into several principle components, in turn, one can identify the presence of the different surface components. Furthermore, the vegetation on the earth is known to share a remarkable reflection signature; the reflection becomes significantly enhanced at wave-lengths longer than 760nm, which is known as a red-edge of the vegetation. If one can identify the corresponding color signature in a pale blue dot, it can be used as a unique probe of the presence of life. I will describe the feasibility of the methodology for future space missions, and consider the direction towards astrobiology from an astrophysicist's point of view.


title: I-mode pedestal relaxation events at ASDEX Upgrade

abstract: The I-mode confinement regime can feature small edge temperature drops that can lead to an increase in the energy deposited onto the divertor targets. In this work, we show that these events are associated with a relaxation of both electron temperature and density edge profiles, with the largest drop found at the pedestal top position. Stability analysis of edge profiles reveals that the operational points are far from the ideal peeling-ballooning boundary. Also, we show that these events appear close to the H-mode transition in the typical I-mode operational space in ASDEX Upgrade, and that no further enhancement of energy confinement is found when they occur. Moreover, scrape-off layer transport during these events is found to be very similar to type-I ELMs, with regard to timescales ($\approx$ 800 $\mu$s), filament propagation, toroidally asymmetric energy effluxes at the midplane and asymmetry between inner and outer divertor deposited energy. In particular, the latter reveals that more energy reaches the outer divertor target. Lastly, first measurements of the divertor peak energy fluence are reported, and projections to ARC - a reactor designed to operate in I-mode - are drawn.


title: Towards Adaptive Benthic Habitat Mapping

abstract: Autonomous Underwater Vehicles (AUVs) are increasingly being used to support scientific research and monitoring studies. One such application is in benthic habitat mapping where these vehicles collect seafloor imagery that complements broadscale bathymetric data collected using sonar. Using these two data sources, the relationship between remotely-sensed acoustic data and the sampled imagery can be learned, creating a habitat model. As the areas to be mapped are often very large and AUV systems collecting seafloor imagery can only sample from a small portion of the survey area, the information gathered should be maximised for each deployment. This paper illustrates how the habitat models themselves can be used to plan more efficient AUV surveys by identifying where to collect further samples in order to most improve the habitat model. A Bayesian neural network is used to predict visually-derived habitat classes when given broad-scale bathymetric data. This network can also estimate the uncertainty associated with a prediction, which can be deconstructed into its aleatoric (data) and epistemic (model) components. We demonstrate how these structured uncertainty estimates can be utilised to improve the model with fewer samples. Such adaptive approaches to benthic surveys have the potential to reduce costs by prioritizing further sampling efforts. We illustrate the effectiveness of the proposed approach using data collected by an AUV on offshore reefs in Tasmania, Australia.


title: The Lernaean Hydra of Data Series Similarity Search: An Experimental Evaluation of the State of the Art

abstract: Increasingly large data series collections are becoming commonplace across many different domains and applications. A key operation in the analysis of data series collections is similarity search, which has attracted lots of attention and effort over the past two decades. Even though several relevant approaches have been proposed in the literature, none of the existing studies provides a detailed evaluation against the available alternatives. The lack of comparative results is further exacerbated by the non-standard use of terminology, which has led to confusion and misconceptions. In this paper, we provide definitions for the different flavors of similarity search that have been studied in the past, and present the first systematic experimental evaluation of the efficiency of data series similarity search techniques. Based on the experimental results, we describe the strengths and weaknesses of each approach and give recommendations for the best approach to use under typical use cases. Finally, by identifying the shortcomings of each method, our findings lay the ground for solid further developments in the field.


title: First-principles study of the effective Hamiltonian for Dirac fermions with spin-orbit coupling in two-dimensional molecular conductor $\alpha$-(BETS)$_2$I$_3$

abstract: We employed first-principles density-functional theory (DFT) calculations to characterize Dirac electrons in quasi-two-dimensional molecular conductor $\alpha$-(BETS)$_2$I$_3$ [= $\alpha$-(BEDT-TSeF)$_2$I$_3$] at a low temperature of 30K. We provide a tight-binding model with intermolecular transfer energies evaluated from maximally localized Wannier functions, where the number of relevant transfer integrals is relatively large due to the delocalized character of Se $p$ orbitals. The spin-orbit coupling gives rise to an exotic insulating state with an indirect band gap of about 2 meV. We analyzed the energy spectrum with a Dirac cone close to the Fermi level to develop an effective Hamiltonian with site-potentials, which reproduces the spectrum obtained by the DFT band structure.


title: Implicit Crowdsourcing for Identifying Abusive Behavior in Online Social Networks

abstract: The increased use of online social networks for the dissemination of information comes with the misuse of the internet for cyberbullying, cybercrime, spam, vandalism, amongst other things. To proactively identify abuse in the networks, we propose a model to identify abusive posts by crowdsourcing. The crowdsourcing part of the detection mechanism is implemented implicitly, by simply observing the natural interaction between users encountering the messages. We explore the node-to-node spread of information on Twitter and propose a model that predicts the abuse level (abusive, hate, spam, normal) associated with the tweet by observing the attributes of the message, along with those of the users interacting with it. We demonstrate that the difference in users' interactions with abusive posts can be leveraged in identifying posts of varying abuse levels.


title: Optimal Mutation Rates for the $(1+\lambda)$ EA on OneMax

abstract: The OneMax problem, alternatively known as the Hamming distance problem, is often referred to as the "drosophila of evolutionary computation (EC)", because of its high relevance in theoretical and empirical analyses of EC approaches. It is therefore surprising that even for the simplest of all mutation-based algorithms, Randomized Local Search and the (1+1) EA, the optimal mutation rates were determined only very recently, in a GECCO 2019 poster. In this work, we extend the analysis of optimal mutation rates to two variants of the $(1+\lambda)$ EA and to the $(1+\lambda)$ RLS. To do this, we use dynamic programming and, for the $(1+\lambda)$ EA, numeric optimization, both requiring $\Theta(n^3)$ time for problem dimension $n$. With this in hand, we compute for all population sizes $\lambda \in \{2^i \mid 0 \le i \le 18\}$ and for problem dimension $n \in \{1000, 2000, 5000\}$ which mutation rates minimize the expected running time and which ones maximize the expected progress. Our results do not only provide a lower bound against which we can measure common evolutionary approaches, but we also obtain insight into the structure of these optimal parameter choices. For example, we show that, for large population sizes, the best number of bits to flip is not monotone in the distance to the optimum. We also observe that the expected remaining running time are not necessarily unimodal for the $(1+\lambda)$ EA$_{0 \rightarrow 1}$ with shifted mutation.


title: Model family selection for classification using Neural Decision Trees

abstract: Model selection consists in comparing several candidate models according to a metric to be optimized. The process often involves a grid search, or such, and cross-validation, which can be time consuming, as well as not providing much information about the dataset itself. In this paper we propose a method to reduce the scope of exploration needed for the task. The idea is to quantify how much it would be necessary to depart from trained instances of a given family, reference models (RMs) carrying `rigid' decision boundaries (e.g. decision trees), so as to obtain an equivalent or better model. In our approach, this is realized by progressively relaxing the decision boundaries of the initial decision trees (the RMs) as long as this is beneficial in terms of performance measured on an analyzed dataset. More specifically, this relaxation is performed by making use of a neural decision tree, which is a neural network built from DTs. The final model produced by our method carries non-linear decision boundaries. Measuring the performance of the final model, and its agreement to its seeding RM can help the user to figure out on which family of models he should focus on.


title: Return of the Lernaean Hydra: Experimental Evaluation of Data Series Approximate Similarity Search

abstract: Data series are a special type of multidimensional data present in numerous domains, where similarity search is a key operation that has been extensively studied in the data series literature. In parallel, the multidimensional community has studied approximate similarity search techniques. We propose a taxonomy of similarity search techniques that reconciles the terminology used in these two domains, we describe modifications to data series indexing techniques enabling them to answer approximate similarity queries with quality guarantees, and we conduct a thorough experimental evaluation to compare approximate similarity search techniques under a unified framework, on synthetic and real datasets in memory and on disk. Although data series differ from generic multidimensional vectors (series usually exhibit correlation between neighboring values), our results show that data series techniques answer approximate %similarity queries with strong guarantees and an excellent empirical performance, on data series and vectors alike. These techniques outperform the state-of-the-art approximate techniques for vectors when operating on disk, and remain competitive in memory.


title: Optimization of Express Train Service Network: Under the Competition of Highway Transportation

abstract: In order to reduce the carbon emission, the related government departments encourage road freights to be transferred more by railway transportation. In China freight transport system, the road transportation is usually responsible for the freights that are in a short distance or the ones with high value-added. To transfer more high value-added freights from highway to railway, except the transportation expenses of railway have an advantage over the road, the transportation time is of certain competitive force as well. Therefore, it is very essential for railway to provide freight train products that are of competitive power. Under such circumstance, a multi-objective programming model of optimizing the rail express train network is devised in this work on the basis of taking both road and railway transportation modes into account. The aims of optimization are to minimize the operation costs of rail trains, and to maximize the railway transport revenue. In a network with a given set of express train services, either the all-or-nothing (AON) method or the logit model can be employed when assigning high value-added freights. These two flow assignment patterns are investigated in this work.


title: PDE-based Dynamic Density Estimation for Large-scale Agent Systems

abstract: Large-scale agent systems have foreseeable applications in the near future. Estimating their macroscopic density is critical for many density-based optimization and control tasks, such as sensor deployment and city traffic scheduling. In this paper, we study the problem of estimating their dynamically varying probability density, given the agents' individual dynamics (which can be nonlinear and time-varying) and their states observed in real-time. The density evolution is shown to satisfy a linear partial differential equation uniquely determined by the agents' dynamics. We present a density filter which takes advantage of the system dynamics to gradually improve its estimation and is scalable to the agents' population. Specifically, we use kernel density estimators (KDE) to construct a noisy measurement and show that, when the agents' population is large, the measurement noise is approximately ``Gaussian''. With this important property, infinite-dimensional Kalman filters are used to design density filters. It turns out that the covariance of measurement noise depends on the true density. This state-dependence makes it necessary to approximate the covariance in the associated operator Riccati equation, rendering the density filter suboptimal. The notion of input-to-state stability is used to prove that the performance of the suboptimal density filter remains close to the optimal one. Simulation results suggest that the proposed density filter is able to quickly recognize the underlying modes of the unknown density and automatically ignore outliers, and is robust to different choices of kernel bandwidth of KDE.


title: Transporting Robotic Swarms via Mean-Field Feedback Control

abstract: With the rapid development of AI and robotics, transporting a large swarm of networked robots has foreseeable applications in the near future. Existing research in swarm robotics has mainly followed a bottom-up philosophy with predefined local coordination and control rules. However, it is arduous to verify the global requirements and analyze their performance. This motivates us to pursue a top-down approach, and develop a provable control strategy for deploying a robotic swarm to achieve a desired global configuration. Specifically, we use mean-field partial differential equations (PDEs) to model the swarm and control its mean-field density (i.e., probability density) over a bounded spatial domain using mean-field feedback. The presented control law uses density estimates as feedback signals and generates corresponding velocity fields that, by acting locally on individual robots, guide their global distribution to a target profile. The design of the velocity field is therefore centralized, but the implementation of the controller can be fully distributed -- individual robots sense the velocity field and derive their own velocity control signals accordingly. The key contribution lies in applying the concept of input-to-state stability (ISS) to show that the perturbed closed-loop system (a nonlinear and time-varying PDE) is locally ISS with respect to density estimation errors. The effectiveness of the proposed control laws is verified using agent-based simulations.


title: Stellar Spectral Interpolation using Machine Learning

abstract: Theoretical stellar spectra rely on model stellar atmospheres computed based on our understanding of the physical laws at play in the stellar interiors. These models, coupled with atomic and molecular line databases, are used to generate theoretical stellar spectral libraries (SSLs) comprising of stellar spectra over a regular grid of atmospheric parameters (temperature, surface gravity, abundances) at any desired resolution. Another class of SSLs is referred to as empirical spectral libraries; these contain observed spectra at limited resolution. SSLs play an essential role in deriving the properties of stars and stellar populations. Both theoretical and empirical libraries suffer from limited coverage over the parameter space. This limitation is overcome to some extent by generating spectra for specific sets of atmospheric parameters by interpolating within the grid of available parameter space. In this work, we present a method for spectral interpolation in the optical region using machine learning algorithms that are generic, easily adaptable for any SSL without much change in the model parameters, and computationally inexpensive. We use two machine learning techniques, Random Forest (RF) and Artificial Neural Networks (ANN), and train the models on the MILES library. We apply the trained models to spectra from the CFLIB for testing and show that the performance of the two models is comparable. We show that both the models achieve better accuracy than the existing methods of polynomial based interpolation and the Gaussian radial basis function (RBF) interpolation.


title: A characterization of $\omega$-limit sets in subshifts of Baire space

abstract: In this paper we consider the structure of $\omega$-limit sets in subshifts of Baire space. We consider both subshifts of finite type and subshifts of bounded type and we demonstrate that many classical structure theorems for $\omega$-limit sets fail in this context. Nevertheless, we obtain characterizations of $\omega$-limit sets in subshift of finite types and of attracting $\omega$-limit sets in subshifts of bounded type.


title: Towards a self-organizing pre-symbolic neural model representing sensorimotor primitives

abstract: The acquisition of symbolic and linguistic representations of sensorimotor behavior is a cognitive process performed by an agent when it is executing and/or observing own and others' actions. According to Piaget's theory of cognitive development, these representations develop during the sensorimotor stage and the pre-operational stage. We propose a model that relates the conceptualization of the higher-level information from visual stimuli to the development of ventral/dorsal visual streams. This model employs neural network architecture incorporating a predictive sensory module based on an RNNPB (Recurrent Neural Network with Parametric Biases) and a horizontal product model. We exemplify this model through a robot passively observing an object to learn its features and movements. During the learning process of observing sensorimotor primitives, i.e. observing a set of trajectories of arm movements and its oriented object features, the pre-symbolic representation is self-organized in the parametric units. These representational units act as bifurcation parameters, guiding the robot to recognize and predict various learned sensorimotor primitives. The pre-symbolic representation also accounts for the learning of sensorimotor primitives in a latent learning context.


title: The existence of a strongly polynomial time simplex algorithm for linear programs

abstract: It is well known that the most challenging question in optimization and discrete geometry is whether there is a strongly polynomial time simplex algorithm for linear programs (LPs). This paper gives a positive answer to this question by using the parameter analysis technique presented by us (arXiv:2006.08104). We show that there is a simplex algorithm whose number of pivoting steps does not exceed the number of variables of a LP problem.


title: Dot product chains

abstract: We study a variant of Erd\H os' unit distance problem, concerning dot products between successive pairs of points chosen from a large finite point set. Specifically, given a large finite set of $n$ points $E$, and a sequence of nonzero dot products $(\alpha_1,\ldots,\alpha_k)$, we give upper and lower bounds on the maximum possible number of tuples of distinct points $(A_1,\dots, A_{k+1})\in E^{k+1}$ satisfying $A_j \cdot A_{j+1}=\alpha_j$ for every $1\leq j \leq k$.


title: Beyond Homophily in Graph Neural Networks: Current Limitations and Effective Designs

abstract: We investigate the representation power of graph neural networks in the semi-supervised node classification task under heterophily or low homophily, i.e., in networks where connected nodes may have different class labels and dissimilar features. Many popular GNNs fail to generalize to this setting, and are even outperformed by models that ignore the graph structure (e.g., multilayer perceptrons). Motivated by this limitation, we identify a set of key designs -- ego- and neighbor-embedding separation, higher-order neighborhoods, and combination of intermediate representations -- that boost learning from the graph structure under heterophily. We combine them into a graph neural network, H2GCN, which we use as the base method to empirically evaluate the effectiveness of the identified designs. Going beyond the traditional benchmarks with strong homophily, our empirical analysis shows that the identified designs increase the accuracy of GNNs by up to 40% and 27% over models without them on synthetic and real networks with heterophily, respectively, and yield competitive performance under homophily.


title: Coupling-based Invertible Neural Networks Are Universal Diffeomorphism Approximators

abstract: Invertible neural networks based on coupling flows (CF-INNs) have various machine learning applications such as image synthesis and representation learning. However, their desirable characteristics such as analytic invertibility come at the cost of restricting the functional forms. This poses a question on their representation power: are CF-INNs universal approximators for invertible functions? Without a universality, there could be a well-behaved invertible transformation that the CF-INN can never approximate, hence it would render the model class unreliable. We answer this question by showing a convenient criterion: a CF-INN is universal if its layers contain affine coupling and invertible linear functions as special cases. As its corollary, we can affirmatively resolve a previously unsolved problem: whether normalizing flow models based on affine coupling can be universal distributional approximators. In the course of proving the universality, we prove a general theorem to show the equivalence of the universality for certain diffeomorphism classes, a theoretical insight that is of interest by itself.


title: Direct Measurement of the Solar-Wind Taylor Microscale using MMS Turbulence Campaign Data

abstract: Using the novel Magnetospheric Multiscale (MMS) mission data accumulated during the 2019 MMS Solar Wind Turbulence Campaign, we calculate the Taylor microscale $(\lambda_{\mathrm{T}})$ of the turbulent magnetic field in the solar wind. The Taylor microscale represents the onset of dissipative processes in classical turbulence theory. An accurate estimation of Taylor scale from spacecraft data is, however, usually difficult due to low time cadence, the effect of time decorrelation, and other factors. Previous reports were based either entirely on the Taylor frozen-in approximation, which conflates time dependence, or that were obtained using multiple datasets, which introduces sample-to-sample variation of plasma parameters, or where inter-spacecraft distance were larger than the present study. The unique configuration of linear formation with logarithmic spacing of the 4 MMS spacecraft, during the campaign, enables a direct evaluation of the $\lambda_{\mathrm{T}}$ from a single dataset, independent of the Taylor frozen-in approximation. A value of $\lambda_{\mathrm{T}} \approx 7000 \, \mathrm{km}$ is obtained, which is about 3 times larger than the previous estimates.


title: Stratonovich-Ito integration scheme in ultrafast spin caloritronics

abstract: The magnonic spin Seebeck effect is a key element of spin caloritronic, a field that exploits thermal effects for spintronic applications. Early studies were focused on investigating the steady-state nonequilibrium magnonic spin Seebeck current, and the underlying physics of the magnonic spin Seebeck effect is now relatively well established. However, the initial steps of the formation of the spin Seebeck current are in the scope of recent interest. To address this dynamical aspect theoretically we propose here a new approach to the time-resolved spin Seebeck effect. Our method exploits the supersymmetric theory of stochastics and Ito - Stratonovich integration scheme. We found that in the early step the spin Seebeck current has both nonzero transversal and longitudinal components. As the magnetization dynamics approaches the steady-state, the transversal components decay through dephasing over the dipole-dipole reservoir. The time scale for this process is typically in the sub-nanoseconds pointing thus to the potential of an ultrafast control of the dynamical spin Seebeck during its buildup.


title: Energy of the Interacting Self-Avoiding Walk at the $\theta-$point

abstract: We perform a numerical study of a new microcanonical polymer model on the three dimensional cubic lattice, consisting of ideal chains whose range and number of nearest-neighbor contacts are fixed to given values. Our simulations suggest an interesting exact relation concerning the internal energy per monomer of the Interacting Self-Avoiding Walk at the $\theta-$point.


title: Characteristics and applications of interplanetary coronal mass ejection composition

abstract: In situ measurements of interplanetary coronal mass ejection (ICME) composition, including elemental abundances and charge states of heavy ions, open a new avenue to study coronal mass ejections (CMEs) besides remote-sensing observations. The ratios between different elemental abundances can diagnose the plasma origin of CMEs (e.g., from the corona or chromosphere/photosphere) due to the first ionization potential (FIP) effect, which means elements with different FIP get fractionated between the photosphere and corona. The ratios between different charge states of a specific element can provide the electron temperature of CMEs in the corona due to the freeze-in effect, which can be used to investigate their eruption process. In this review, we first give an overview of the ICME composition and then demonstrate their applications in investigating some important subjects related to CMEs, such as the origin of filament plasma and the eruption process of magnetic flux ropes. Finally, we point out several important questions that should be addressed further for better utilizing the ICME composition to study CMEs.


title: Coconut: sortable summarizations for scalable indexes over static and streaming data series

abstract: Many modern applications produce massive streams of data series that need to be analyzed, requiring efficient similarity search operations. However, the state-of-the-art data series indexes that are used for this purpose do not scale well for massive datasets in terms of performance, or storage costs. We pinpoint the problem to the fact that existing summarizations of data series used for indexing cannot be sorted while keeping similar data series close to each other in the sorted order. To address this problem, we present Coconut, the first data series index based on sortable summarizations and the first efficient solution for indexing and querying streaming series. The first innovation in Coconut is an inverted, sortable data series summarization that organizes data series based on a z-order curve, keeping similar series close to each other in the sorted order. As a result, Coconut is able to use bulk loading and updating techniques that rely on sorting to quickly build and maintain a contiguous index using large sequential disk I/Os. We then explore prefix-based and median-based splitting policies for bottom-up bulk loading, showing that median-based splitting outperforms the state of the art, ensuring that all nodes are densely populated. Finally, we explore the impact of sortable summarizations on variable-sized window queries, showing that they can be supported in the presence of updates through efficient merging of temporal partitions. Overall, we show analytically and empirically that Coconut dominates the state-of-the-art data series indexes in terms of construction speed, query speed, and storage costs.


title: Zeros of a binomial combination of Chebyshev polynomials

abstract: For $0<\alpha<1$, we study the zeros of the sequence of polynomials $\left\{ P_{m}(z)\right\} _{m=0}^{\infty}$ generated by the reciprocal of $(1-t)^{\alpha}(1-2zt+t^{2})$, expanded as a power series in $t$. Equivalently, this sequence is obtained from a linear combination of Chebyshev polynomials whose coefficients have a binomial form. We show that the number of zeros of $P_{m}(z)$ outside the interval $(-1,1)$ is bounded by a constant independent of $m$.


title: Video Playback Rate Perception for Self-supervisedSpatio-Temporal Representation Learning

abstract: In self-supervised spatio-temporal representation learning, the temporal resolution and long-short term characteristics are not yet fully explored, which limits representation capabilities of learned models. In this paper, we propose a novel self-supervised method, referred to as video Playback Rate Perception (PRP), to learn spatio-temporal representation in a simple-yet-effective way. PRP roots in a dilated sampling strategy, which produces self-supervision signals about video playback rates for representation model learning. PRP is implemented with a feature encoder, a classification module, and a reconstructing decoder, to achieve spatio-temporal semantic retention in a collaborative discrimination-generation manner. The discriminative perception model follows a feature encoder to prefer perceiving low temporal resolution and long-term representation by classifying fast-forward rates. The generative perception model acts as a feature decoder to focus on comprehending high temporal resolution and short-term representation by introducing a motion-attention mechanism. PRP is applied on typical video target tasks including action recognition and video retrieval. Experiments show that PRP outperforms state-of-the-art self-supervised models with significant margins. Code is available at github.com/yuanyao366/PRP


title: wav2vec 2.0: A Framework for Self-Supervised Learning of Speech Representations

abstract: We show for the first time that learning powerful representations from speech audio alone followed by fine-tuning on transcribed speech can outperform the best semi-supervised methods while being conceptually simpler. wav2vec 2.0 masks the speech input in the latent space and solves a contrastive task defined over a quantization of the latent representations which are jointly learned. Experiments using all labeled data of Librispeech achieve 1.8/3.3 WER on the clean/other test sets. When lowering the amount of labeled data to one hour, wav2vec 2.0 outperforms the previous state of the art on the 100 hour subset while using 100 times less labeled data. Using just ten minutes of labeled data and pre-training on 53k hours of unlabeled data still achieves 4.8/8.2 WER. This demonstrates the feasibility of speech recognition with limited amounts of labeled data.


title: Representation via Representations: Domain Generalization via Adversarially Learned Invariant Representations

abstract: We investigate the power of censoring techniques, first developed for learning {\em fair representations}, to address domain generalization. We examine {\em adversarial} censoring techniques for learning invariant representations from multiple "studies" (or domains), where each study is drawn according to a distribution on domains. The mapping is used at test time to classify instances from a new domain. In many contexts, such as medical forecasting, domain generalization from studies in populous areas (where data are plentiful), to geographically remote populations (for which no training data exist) provides fairness of a different flavor, not anticipated in previous work on algorithmic fairness. We study an adversarial loss function for $k$ domains and precisely characterize its limiting behavior as $k$ grows, formalizing and proving the intuition, backed by experiments, that observing data from a larger number of domains helps. The limiting results are accompanied by non-asymptotic learning-theoretic bounds. Furthermore, we obtain sufficient conditions for good worst-case prediction performance of our algorithm on previously unseen domains. Finally, we decompose our mappings into two components and provide a complete characterization of invariance in terms of this decomposition. To our knowledge, our results provide the first formal guarantees of these kinds for adversarial invariant domain generalization.


title: Compiler Directed Speculative Intermittent Computation

abstract: This paper presents CoSpec, a new architecture/compiler co-design scheme that works for commodity in-order processors used in energy-harvesting systems. To achieve crash consistency without requiring unconventional architectural support, CoSpec leverages speculation assuming that power failure is not going to occur and thus holds all committed stores in a store buffer (SB), as if they were speculative, in case of mispeculation. CoSpec compiler first partitions a given program into a series of recoverable code regions with the SB size in mind, so that no region overflows the SB. When the program control reaches the end of each region, the speculation turns out to be successful, thus releasing all the buffered stores of the region to NVM. If power failure occurs during the execution of a region, all its speculative stores disappear in the volatile SB, i.e., they never affect program states in NVM. Consequently, the interrupted region can be restarted with consistent program states in the wake of power failure. To hide the latency of the SB release, i.e., essentially NVM writes, at each region boundary, CoSpec overlaps the NVM writes of the current region with the speculative execution of the next region. Such instruction level parallelism gives an illusion of out-of-order execution on top of the in-order processor, achieving a speedup of more than 1.2X when there is no power outage. Our experiments on a set of real energy harvesting traces with frequent outages demonstrate that CoSpec outperforms the state-of-the-art scheme by 1.8~3X on average.


title: Unsupervised Image Classification for Deep Representation Learning

abstract: Deep clustering against self-supervised learning is a very important and promising direction for unsupervised visual representation learning since it requires little domain knowledge to design pretext tasks. However, the key component, embedding clustering, limits its extension to the extremely large-scale dataset due to its prerequisite to save the global latent embedding of the entire dataset. In this work, we aim to make this framework more simple and elegant without performance decline. We propose an unsupervised image classification framework without using embedding clustering, which is very similar to standard supervised training manner. For detailed interpretation, we further analyze its relation with deep clustering and contrastive learning. Extensive experiments on ImageNet dataset have been conducted to prove the effectiveness of our method. Furthermore, the experiments on transfer learning benchmarks have verified its generalization to other downstream tasks, including multi-label image classification, object detection, semantic segmentation and few-shot image classification.


title: Pseudo-LiDAR Point Cloud Interpolation Based on 3D Motion Representation and Spatial Supervision

abstract: Pseudo-LiDAR point cloud interpolation is a novel and challenging task in the field of autonomous driving, which aims to address the frequency mismatching problem between camera and LiDAR. Previous works represent the 3D spatial motion relationship induced by a coarse 2D optical flow, and the quality of interpolated point clouds only depends on the supervision of depth maps. As a result, the generated point clouds suffer from inferior global distributions and local appearances. To solve the above problems, we propose a Pseudo-LiDAR point cloud interpolation network to generates temporally and spatially high-quality point cloud sequences. By exploiting the scene flow between point clouds, the proposed network is able to learn a more accurate representation of the 3D spatial motion relationship. For the more comprehensive perception of the distribution of point cloud, we design a novel reconstruction loss function that implements the chamfer distance to supervise the generation of Pseudo-LiDAR point clouds in 3D space. In addition, we introduce a multi-modal deep aggregation module to facilitate the efficient fusion of texture and depth features. As the benefits of the improved motion representation, training loss function, and model structure, our approach gains significant improvements on the Pseudo-LiDAR point cloud interpolation task. The experimental results evaluated on KITTI dataset demonstrate the state-of-the-art performance of the proposed network, quantitatively and qualitatively.


title: A Bakry-\'Emery Almost Splitting Result With Applications to the Topology of Black Holes

abstract: The almost splitting theorem of Cheeger-Colding is established in the setting of almost nonnegative generalized $m$-Bakry-\'{E}mery Ricci curvature, in which $m$ is positive and the associated vector field is not necessarily required to be the gradient of a function. In this context it is shown that with a diameter upper bound and volume lower bound the fundamental group of such manifolds is almost abelian. Furthermore, extensions of well-known results concerning Ricci curvature lower bounds are given for generalized $m$-Bakry-\'{E}mery Ricci curvature. These include: the first Betti number bound of Gromov and Gallot, Anderson's finiteness of fundamental group isomorphism types, volume comparison, the Abresch-Gromoll inequality, and a Cheng-Yau gradient estimate. Finally, this analysis is applied to stationary vacuum black holes in higher dimensions to find that low temperature horizons must have limited topology, similar to the restrictions exhibited by (extreme) horizons of zero temperature.


title: Predicting Temporal Sets with Deep Neural Networks

abstract: Given a sequence of sets, where each set contains an arbitrary number of elements, the problem of temporal sets prediction aims to predict the elements in the subsequent set. In practice, temporal sets prediction is much more complex than predictive modelling of temporal events and time series, and is still an open problem. Many possible existing methods, if adapted for the problem of temporal sets prediction, usually follow a two-step strategy by first projecting temporal sets into latent representations and then learning a predictive model with the latent representations. The two-step approach often leads to information loss and unsatisfactory prediction performance. In this paper, we propose an integrated solution based on the deep neural networks for temporal sets prediction. A unique perspective of our approach is to learn element relationship by constructing set-level co-occurrence graph and then perform graph convolutions on the dynamic relationship graphs. Moreover, we design an attention-based module to adaptively learn the temporal dependency of elements and sets. Finally, we provide a gated updating mechanism to find the hidden shared patterns in different sequences and fuse both static and dynamic information to improve the prediction performance. Experiments on real-world data sets demonstrate that our approach can achieve competitive performances even with a portion of the training data and can outperform existing methods with a significant margin.


title: Distinguishing resonance symmetries with energy-resolved photoion angular distributions from ion-pair formation in O$_2$ following two-photon absorption of a 9.3 eV femtosecond pulse

abstract: We present a combined experimental and theoretical study on the photodissociation dynamics of ion-pair formation in O$_2$ following resonant two-photon absorption of a 9.3 eV femtosecond pulse, where the resulting O$^+$ ions are detected using 3-D momentum imaging. Ion-pair formation states of $^3\Sigma^-_g$ and $^3\Pi_g$ symmetry are accessed through predissociation of optically dark continuum Rydberg states converging to the B $^2\Sigma^-_g$ ionic state, which are resonantly populated via a mixture of both parallel-parallel and parallel-perpendicular two-photon transitions. This mixture is evident in the angular distribution of the dissociation relative to the light polarization, and varies with the kinetic energy release (KER) of the fragmenting ion-pair. The KER-dependent photoion angular distribution reveals the underlying two-photon absorption dynamics involved in the ion-pair production mechanism and indicates the existence of two nearly degenerate continuum resonances possessing different symmetries, which can both decay by coupling to ion-pair states of the same total symmetry through internal conversion.


title: Generating Adjacency-Constrained Subgoals in Hierarchical Reinforcement Learning

abstract: Goal-conditioned hierarchical reinforcement learning (HRL) is a promising approach for scaling up reinforcement learning (RL) techniques. However, it often suffers from training inefficiency as the action space of the high-level, i.e., the goal space, is often large. Searching in a large goal space poses difficulties for both high-level subgoal generation and low-level policy learning. In this paper, we show that this problem can be effectively alleviated by restricting the high-level action space from the whole goal space to a $k$-step adjacent region of the current state using an adjacency constraint. We theoretically prove that the proposed adjacency constraint preserves the optimal hierarchical policy in deterministic MDPs, and show that this constraint can be practically implemented by training an adjacency network that can discriminate between adjacent and non-adjacent subgoals. Experimental results on discrete and continuous control tasks show that incorporating the adjacency constraint improves the performance of state-of-the-art HRL approaches in both deterministic and stochastic environments.


title: Unsupervised Vehicle Re-identification with Progressive Adaptation

abstract: Vehicle re-identification (reID) aims at identifying vehicles across different non-overlapping cameras views. The existing methods heavily relied on well-labeled datasets for ideal performance, which inevitably causes fateful drop due to the severe domain bias between the training domain and the real-world scenes; worse still, these approaches required full annotations, which is labor-consuming. To tackle these challenges, we propose a novel progressive adaptation learning method for vehicle reID, named PAL, which infers from the abundant data without annotations. For PAL, a data adaptation module is employed for source domain, which generates the images with similar data distribution to unlabeled target domain as ``pseudo target samples''. These pseudo samples are combined with the unlabeled samples that are selected by a dynamic sampling strategy to make training faster. We further proposed a weighted label smoothing (WLS) loss, which considers the similarity between samples with different clusters to balance the confidence of pseudo labels. Comprehensive experimental results validate the advantages of PAL on both VehicleID and VeRi-776 dataset.


title: Paying more attention to snapshots of Iterative Pruning: Improving Model Compression via Ensemble Distillation

abstract: Network pruning is one of the most dominant methods for reducing the heavy inference cost of deep neural networks. Existing methods often iteratively prune networks to attain high compression ratio without incurring significant loss in performance. However, we argue that conventional methods for retraining pruned networks (i.e., using small, fixed learning rate) are inadequate as they completely ignore the benefits from snapshots of iterative pruning. In this work, we show that strong ensembles can be constructed from snapshots of iterative pruning, which achieve competitive performance and vary in network structure. Furthermore, we present simple, general and effective pipeline that generates strong ensembles of networks during pruning with large learning rate restarting, and utilizes knowledge distillation with those ensembles to improve the predictive power of compact models. In standard image classification benchmarks such as CIFAR and Tiny-Imagenet, we advance state-of-the-art pruning ratio of structured pruning by integrating simple l1-norm filters pruning into our pipeline. Specifically, we reduce 75-80% of total parameters and 65-70% MACs of numerous variants of ResNet architectures while having comparable or better performance than that of original networks. Code associate with this paper is made publicly available at https://github.com/lehduong/kesi.


title: Recovering Accurate Labeling Information from Partially Valid Data for Effective Multi-Label Learning

abstract: Partial Multi-label Learning (PML) aims to induce the multi-label predictor from datasets with noisy supervision, where each training instance is associated with several candidate labels but only partially valid. To address the noisy issue, the existing PML methods basically recover the ground-truth labels by leveraging the ground-truth confidence of the candidate label, \ie the likelihood of a candidate label being a ground-truth one. However, they neglect the information from non-candidate labels, which potentially contributes to the ground-truth label recovery. In this paper, we propose to recover the ground-truth labels, \ie estimating the ground-truth confidences, from the label enrichment, composed of the relevance degrees of candidate labels and irrelevance degrees of non-candidate labels. Upon this observation, we further develop a novel two-stage PML method, namely \emph{\underline{P}artial \underline{M}ulti-\underline{L}abel \underline{L}earning with \underline{L}abel \underline{E}nrichment-\underline{R}ecovery} (\baby), where in the first stage, it estimates the label enrichment with unconstrained label propagation, then jointly learns the ground-truth confidence and multi-label predictor given the label enrichment. Experimental results validate that \baby outperforms the state-of-the-art PML methods.


title: Federated Learning Meets Multi-objective Optimization

abstract: Federated learning has emerged as a promising, massively distributed way to train a joint deep model over large amounts of edge devices while keeping private user data strictly on device. In this work, motivated from ensuring fairness among users and robustness against malicious adversaries, we formulate federated learning as multi-objective optimization and propose a new algorithm FedMGDA+ that is guaranteed to converge to Pareto stationary solutions. FedMGDA+ is simple to implement, has fewer hyperparameters to tune, and refrains from sacrificing the performance of any participating user. We establish the convergence properties of FedMGDA+ and point out its connections to existing approaches. Extensive experiments on a variety of datasets confirm that FedMGDA+ compares favorably against state-of-the-art.


title: Coherent control of quantum and entanglement dynamics via periodic modulations in optomechanical semi-conductor resonator coupled to quantum-dot excitons

abstract: We systematically study the influence of simultaneously modulating the input laser intensity and quantum dot (QD) resonance frequecy on the mean-field dynamics, fluctuation energy transfer and entanglement in a optomechanical semi-conductor resonator embedded with a QD. We show that the modulation and the hybrid system can be engineered to attain the desired mean-field values, control the fluctuation energy transfer and the entanglement between the various degrees of freedom. A remarkably high degree of entanglement can be achieved by modulating only the QD frequency. The interplay between the two modulations leads to an entanglement which lies between that generated solely by modulating either the QD or the pump laser intensity. A transition from low stationary to large dynamical entanglement occurs as we switch on the modulation. This study opens up new possibilities for optimal control strategies and can be used for data signal transfer and storage in quantum communication platforms.


title: Differential operators on quantized flag manifolds at roots of unity III

abstract: We describe the cohomology of the sheaf of twisted differential operators on the quantized flag manifold at a root of unity whose order is a prime power. It follows from this and our previous results that for the De Concini-Kac type quantized enveloping algebra, where the parameter $q$ is specialized to a root of unity whose order is a prime power, the number of irreducible modules with a certain specified central character coincides with the dimension of the total cohomology group of the corresponding Springer fiber. This gives a weak version of a conjecture of Lusztig concerning non-restricted representations of the quantized enveloping algebra.


title: A Distributed Multi-Vehicle Coordination Algorithm for Navigation in Tight Environments

abstract: This work presents a distributed method for multi-vehicle coordination based on nonlinear model predictive control (NMPC) and dual decomposition. Our approach allows the vehicles to coordinate in tight spaces (e.g., busy highway lanes or parking lots) by using a polytopic description of each vehicle's shape and formulating collision avoidance as a dual optimization problem. Our method accommodates heterogeneous teams of vehicles (i.e., vehicles with different polytopic shapes and dynamic models can be part of the same team). Our method allows the vehicles to share their intentions in a distributed fashion without relying on a central coordinator and efficiently provides collision-free trajectories for the vehicles. In addition, our method decouples the individual-vehicles' trajectory optimization from their collision-avoidance objectives enhancing the scalability of the method and allowing one to exploit parallel hardware architectures. All these features are particularly important for vehicular applications, where the systems operate at high-frequency rates in dynamic environments. To validate our method, we apply it in a vehicular application, that is, the autonomous lane-merging of a team of connected vehicles to form a platoon. We compare our design with the centralized NMPC design to show the computational benefits of the proposed distributed algorithm.


title: Real-time LCC-HVDC Maximum Emergency Power Capacity Estimation Based on Local PMU Measurements

abstract: The adjustable capacity of a line-commutated-converter High Voltage Direct Current (LCC-HVDC) connected to a power system, called the LCC-HVDC maximum emergency power capability or HVDC-MC for short, plays an important role in determining the response of that system to a large disturbance. However, it is a challenging task to obtain an accurate HVDC-MC due to system model uncertainties as well as to contingencies. To address this problem, this paper proposes to estimate the HVDC-MC using a Thevenin equivalent (TE) of the system seen at the HVDC terminal bus of connection with the power system, whose parameters are estimated by processing positive-sequences voltages and currents of local synchrophasor measurements. The impacts of TE potential changes on the impedance estimation under large disturbance have been extensively investigated and an adaptive screening process of current measurements is developed to reduce the error of TE impedance estimation. The uncertainties of phasor measurements have been further taken into account by resorting to the total least square estimation method. The limitations of the HVDC control characteristics, the voltage-dependent current order limit, the converter capacity, and the AC voltage on HVDC-MC estimation are also considered. The simulations show that the proposed method can accurately track the dynamics of the TE parameters and the real-time HVDC-MC after the large disturbances.


title: AOT: Pushing the Efficiency Boundary of Main-memory Triangle Listing

abstract: Triangle listing is an important topic significant in many practical applications. Efficient algorithms exist for the task of triangle listing. Recent algorithms leverage an orientation framework, which can be thought of as mapping an undirected graph to a directed acylic graph, namely oriented graph, with respect to any global vertex order. In this paper, we propose an adaptive orientation technique that satisfies the orientation technique but refines it by traversing carefully based on the out-degree of the vertices in the oriented graph during the computation of triangles. Based on this adaptive orientation technique, we design a new algorithm, namely aot, to enhance the edge-iterator listing paradigm. We also make improvements to the performance of aot by exploiting the local order within the adjacent list of the vertices. We show that aot is the first work which can achieve best performance in terms of both practical performance and theoretical time complexity. Our comprehensive experiments over $16$ real-life large graphs show a superior performance of our \aot algorithm when compared against the state-of-the-art, especially for massive graphs with billions of edges. Theoretically, we show that our proposed algorithm has a time complexity of $\Theta(\sum_{ \langle u,v \rangle \in \vec{E} } \min\{ deg^{+}(u),deg^{+}(v)\}))$, where $\vec{E}$ and $deg^{+}(x)$ denote the set of directed edges in an oriented graph and the out-degree of vertex $x$ respectively. As to our best knowledge, this is the best time complexity among in-memory triangle listing algorithms.


title: Predictions in the eye of the beholder: an active inference account of Watt governors

abstract: Active inference introduces a theory describing action-perception loops via the minimisation of variational (and expected) free energy or, under simplifying assumptions, (weighted) prediction error. Recently, active inference has been proposed as part of a new and unifying framework in the cognitive sciences: predictive processing. Predictive processing is often associated with traditional computational theories of the mind, strongly relying on internal representations presented in the form of generative models thought to explain different functions of living and cognitive systems. In this work, we introduce an active inference formulation of the Watt centrifugal governor, a system often portrayed as the canonical "anti-representational" metaphor for cognition. We identify a generative model of a steam engine for the governor, and derive a set of equations describing "perception" and "action" processes as a form of prediction error minimisation. In doing so, we firstly challenge the idea of generative models as explicit internal representations for cognitive systems, suggesting that such models serve only as implicit descriptions for an observer. Secondly, we consider current proposals of predictive processing as a theory of cognition, focusing on some of its potential shortcomings and in particular on the idea that virtually any system admits a description in terms of prediction error minimisation, suggesting that this theory may offer limited explanatory power for cognitive systems. Finally, as a silver lining we emphasise the instrumental role this framework can nonetheless play as a mathematical tool for modelling cognitive architectures interpreted in terms of Bayesian (active) inference.


title: Measure-resend authenticated semi-quantum key distribution with single photons

abstract: Yu et al. and Li et al. have proposed the measure-resend protocols of authenticated semi-quantum key distribution (ASQKD). A new measure-resend ASQKD protocol is proposed in this paper, which requires a lower burden of quantum resource, needs fewer bits of the pre-shared key, and even provides better qubit efficiency than their protocols. The security proof shows the robustness of the proposed protocol under the collective attack.


title: Euler equations for Cosserat media

abstract: We consider Cosserat media as SO(3)-structures over a domain D in R3. Motions of such media are given by automorphisms of the SO(3)-bundle. We present Euler-type equations for such media and discuss their structure.


title: Ancilla-assisted frequency estimation under phase covariant noises with Greenberger-Horne-Zeilinger states

abstract: It has been demonstrated that the optimal sensitivity achievable with Greenberger-Horne-Zeilinger states is the same as that with uncorrelated probes in the frequency estimation in the presence of uncorrelated Markovian dephasing [S. F. Huelga, et al., Phys. Rev. Lett. 79, 3865 (1997)]. Here, we extend this issue by examining the optimal frequency sensitivities achievable by the use of ancilla-assisted strategy, which has been proposed recently for robust phase estimation. We present the ultimate frequency sensitivities bounded by the quantum Fisher information for a general case in the presence of Markovian covariant phase noises, and the optimal measurement observables that can saturate the theoretical sensitivity bounds. We also demonstrate the effectiveness of the ancilla-assisted strategy for preserving frequency sensitivities suffering from specific physically ground noises.


title: Algorithms to enumerate superspecial Howe curves of genus 4

abstract: A Howe curve is a curve of genus $4$ obtained as the fiber product of two genus-$1$ double covers of $\mathbf{P}^1$. In this paper, we present a simple algorithm for testing isomorphism of Howe curves, and we propose two main algorithms for finding and enumerating these curves: One involves solving multivariate systems coming from Cartier--Manin matrices, while the other uses Richelot isogenies of curves of genus $2$. Comparing the two algorithms by implementation and by complexity analyses, we conclude that the latter enumerates curves more efficiently. Using these algorithms, we show that there exist superspecial curves of genus $4$ in characteristic $p$ for every prime $p$ with $7 < p < 20000$.


title: On some enriched contractions in Banach spaces

abstract: In this paper, we introduce two new types of enriched contractions, viz., enriched $\mathcal{A}$-contraction and enriched $\mathcal{A}'$-contraction. Then we obtain fixed points of mappings satisfying such contractions using the fixed point property of the average operator of the mappings. Further, we study the well-posedness and limit shadowing property of the fixed point problem involving the contractions, and give some examples to validate the results proved. We frame an open question related to the existence of a fixed point of such contractions. We also show that Berinde and P\u{a}curar's recent results on different kinds enriched contractions and some well known classical fixed point results are particular cases of our results.


title: Dielectric constant of supercritical water in a large pressure-temperature range

abstract: A huge amount of water at supercritical conditions exists in Earth's interior, where its dielectric properties play a critical role in determining how it stores and transports materials. However, it is very challenging to obtain the static dielectric constant of water, $\epsilon_0$, in a wide pressure-temperature (P-T) range as found in deep Earth either experimentally or by first-principles simulations. Here, we introduce a neural network dipole model, which, combined with molecular dynamics, can be used to compute P-T dependent dielectric properties of water as accurately as first-principles methods but much more efficiently. We found that $\epsilon_0$ may vary by one order of magnitude in Earth's upper mantle, suggesting that the solvation properties of water change dramatically at different depths. There is a subtle interplay between the molecular dipole moment and the dipolar angular correlation in governing the change of $\epsilon_0$. We also calculated the frequency-dependent dielectric constant of water in the microwave range, which, to the best of our knowledge, has not been calculated from first principles, and found that temperature affects the dielectric absorption more than pressure. Our results are of great use in many areas, e.g., modelling water-rock interactions in geochemistry. The computational approach introduced here can be readily applied to other molecular fluids.


title: von Neumann's Minimax Theorem for Continuous Quantum Games

abstract: The concept of a classical player, corresponding to a classical random variable, is extended to include quantum random variables in the form of self adjoint operators on infinite dimensional Hilbert space. A quantum version of Von Neumann's Minimax theorem for infinite dimensional (or continuous) games is proved.