title: A predictive safety filter for learning-based racing control

abstract: The growing need for high-performance controllers in safety-critical applications like autonomous driving has been motivating the development of formal safety verification techniques. In this paper, we design and implement a predictive safety filter that is able to maintain vehicle safety with respect to track boundaries when paired alongside any potentially unsafe control signal, such as those found in learning-based methods. A model predictive control (MPC) framework is used to create a minimally invasive algorithm that certifies whether a desired control input is safe and can be applied to the vehicle, or that provides an alternate input to keep the vehicle in bounds. To this end, we provide a principled procedure to compute a safe and invariant set for nonlinear dynamic bicycle models using efficient convex approximation techniques. To fully support an aggressive racing performance without conservative safety interventions, the safe set is extended in real-time through predictive control backup trajectories. Applications for assisted manual driving and deep imitation learning on a miniature remote-controlled vehicle demonstrate the safety filter's ability to ensure vehicle safety during aggressive maneuvers.


title: Combined Iodine- and Sulfur-based Treatments for an Effective Passivation of GeSn Surface

abstract: GeSn alloys are metastable semiconductors that have been proposed as building blocks for silicon-integrated short-wave and mid-wave infrared photonic and sensing platforms. Exploiting these semiconductors requires, however, the control of their epitaxy and their surface chemistry to reduce non-radiative recombination that hinders the efficiency of optoelectronic devices. Herein, we demonstrate that a combined sulfur- and iodine-based treatments yields effective passivation of Ge and Ge0.9Sn0.1 surfaces. X-ray photoemission spectroscopy and in situ spectroscopic ellipsometry measurements were used to investigate the dynamics of surface stability and track the reoxidation mechanisms. Our analysis shows that the largest reduction in oxide after HI treatment, while HF+(NH4)2S results in a lower re-oxidation rate. A combined HI+(NH4)2S treatment preserves the lowest oxide ratio <10 % up to 1 hour of air exposure, while less than half of the initial oxide coverage is reached after 4 hours. These results highlight the potential of S- and I-based treatments in stabilizing the GeSn surface chemistry thus enabling a passivation method that is compatible with materials and device processing.


title: Unit Tangent Bundles, CR Leaf Spaces, and Hypercomplex Structures

abstract: Unit tangent bundles $UM$ of semi-Riemannian manifolds $M$ are shown to be examples of dynamical Legendrian contact structures, which were defined in recent work [25] of Sykes-Zelenko to generalize leaf spaces of 2-nondegenerate CR manifolds. In doing so, Sykes-Zelenko extended the classification in Porter-Zelenko [20] of regular, 2-nondegenerate CR structures to those that can be recovered from their leaf space. The present paper treats dynamical Legendrian contact structures associated with 2-nondegenerate CR structures which were called "strongly regular" in Porter-Zelenko, named "L-contact structures." Closely related to Lie-contact structures, L-contact manifolds have homogeneous models given by isotropic Grassmannians of complex 2-planes whose algebra of infinitesimal symmetries is one of $\mathfrak{so}(p+2,q+2)$ or $\mathfrak{so}^*(2p+4)$ for $p\geq1$, $q\geq0$. Each 2-plane in the homogeneous model is a split-quaternionic or quaternionic line, respectively, and more general L-contact structures arise on contact manifolds with hypercomplex structures, unit tangent bundles being a prime example. The Ricci curvature tensor of $M$ is used to define the "Ricci-shifted" L-contact structure on $UM$, whose Nijenhuis tensor vanishes when $M$ is conformally flat. In the language of Sykes-Zelenko (for $M$ analytic), such $UM$ is the leaf space of a 2-nondegenerate CR manifold which is "recoverable" from $UM$, providing a new source of examples of 2-nondegenerate CR structures.


title: Planar potential flow on Cartesian grids

abstract: Potential flow has many applications, including the modelling of unsteady flows in aerodynamics. For these models to work efficiently, it is best to avoid Biot-Savart interactions. This work presents a grid-based treatment of potential flows in two dimensions and its use in a vortex model for simulating unsteady aerodynamic flows. For flows consisting of vortex elements, the treatment follows the vortex-in-cell approach and solves the streamfunction-vorticity Poisson equation on a Cartesian grid after transferring the circulation from the vortices onto the grid. For sources and sinks, an analogous approach can be followed using the scalar potential. The combined velocity field due to vortices, sinks, and sources can then be obtained using the Helmholtz decomposition. In this work, we use several key tools that ensure the approach works on arbitrary geometries, with and without sharp edges. Firstly, the immersed boundary projection method is used to account for bodies in the flow and the resulting body-forcing Lagrange multiplier is identified as the bound vortex sheet strength. Secondly, sharp edges are treated by decomposing the vortex sheet strength into a singular and smooth part. To enforce the Kutta condition, the smooth part can then be constrained to remove the singularity introduced by the sharp edge. These constraints and the Poisson equation are formulated as a saddle-point system and solved using the Schur complement method. The lattice Green's function is used to efficiently solve the discrete Poisson equation with unbounded boundary conditions. The method and its accuracy are demonstrated for several problems.


title: Instance Specific Approximations for Submodular Maximization

abstract: For many optimization problems in machine learning, finding an optimal solution is computationally intractable and we seek algorithms that perform well in practice. Since computational intractability often results from pathological instances, we look for methods to benchmark the performance of algorithms against optimal solutions on real-world instances. The main challenge is that an optimal solution cannot be efficiently computed for intractable problems, and we therefore often do not know how far a solution is from being optimal. A major question is therefore how to measure the performance of an algorithm in comparison to an optimal solution on instances we encounter in practice. In this paper, we address this question in the context of submodular optimization problems. For the canonical problem of submodular maximization under a cardinality constraint, it is intractable to compute a solution that is better than a $1-1/e \approx 0.63$ fraction of the optimum. Algorithms like the celebrated greedy algorithm are guaranteed to achieve this $1-1/e$ bound on any instance and are used in practice. Our main contribution is not a new algorithm for submodular maximization but an analytical method that measures how close an algorithm for submodular maximization is to optimal on a given problem instance. We use this method to show that on a wide variety of real-world datasets and objectives, the approximation of the solution found by greedy goes well beyond $1-1/e$ and is often at least 0.95. We develop this method using a novel technique that lower bounds the objective of a dual minimization problem to obtain an upper bound on the value of an optimal solution to the primal maximization problem.


title: On Big Pieces approximations of parabolic hypersurfaces

abstract: Let $\Sigma$ be a closed subset of $\mathbb{R}^ {n+1}$ which is parabolic Ahlfors-David regular and assume that $\Sigma$ satisfies a 2-sided corkscrew condition. Assume, in addition, that $\Sigma$ is either time-forwards Ahlfors-David regular, time-backwards Ahlfors-David regular, or parabolic uniform rectifiable. We then first prove that $\Sigma$ satisfies a {\it weak synchronized two cube condition}. Based on this we are able to revisit the argument in \cite{NS} and prove that $\Sigma$ contains {\it uniform big pieces of Lip(1,1/2) graphs}. When $\Sigma$ is parabolic uniformly rectifiable the construction can be refined and in this case we prove that $\Sigma$ contains {\it uniform big pieces of regular parabolic Lip(1,1/2) graphs}. Similar results hold if $\Omega\subset\mathbb R^{n+1}$ is a connected component of $\mathbb R^{n+1}\setminus\Sigma$ and in this context we also give a parabolic counterpart of the main result in \cite{AHMNT} by proving that if $\Omega$ is a one-sided parabolic chord arc domain, and if $\Sigma$ is parabolic uniformly rectifiable, then $\Omega$ is in fact a parabolic chord arc domain. Our results give a flexible parabolic version of the classical (elliptic) result of G. David and D. Jerison concerning the existence of uniform big pieces of Lipschitz graphs for sets satisfying a two disc condition.


title: Are locally finite MV-algebras a variety?

abstract: We answer Mundici's problem number 3 (D. Mundici. Advanced {\L}ukasiewicz calculus. Trends in Logic Vol. 35. Springer 2011, p. 235): Is the category of locally finite MV-algebras equivalent to an equational class? We prove: (i) The category of locally finite MV-algebras is not equivalent to any finitary variety. (ii) More is true: the category of locally finite MV-algebras is not equivalent to any finitely-sorted finitary quasi-variety. (iii) The category of locally finite MV-algebras is equivalent to an infinitary variety; with operations of at most countable arity. (iv) The category of locally finite MV-algebras is equivalent to a countably-sorted finitary variety. Our proofs rest upon the duality between locally finite MV-algebras and the category of multisets by R. Cignoli, E. J. Dubuc and D. Mundici, and categorical characterisations of varieties and quasi-varieties proved by J. Duskin, J. R. Isbell, F. W. Lawvere and others. In fact no knowledge on MV-algebras is needed, apart from the aforementioned duality.


title: A Robotic Model of Hippocampal Reverse Replay for Reinforcement Learning

abstract: Hippocampal reverse replay is thought to contribute to learning, and particularly reinforcement learning, in animals. We present a computational model of learning in the hippocampus that builds on a previous model of the hippocampal-striatal network viewed as implementing a three-factor reinforcement learning rule. To augment this model with hippocampal reverse replay, a novel policy gradient learning rule is derived that associates place cell activity with responses in cells representing actions. This new model is evaluated using a simulated robot spatial navigation task inspired by the Morris water maze. Results show that reverse replay can accelerate learning from reinforcement, whilst improving stability and robustness over multiple trials. As implied by the neurobiological data, our study implies that reverse replay can make a significant positive contribution to reinforcement learning, although learning that is less efficient and less stable is possible in its absence. We conclude that reverse replay may enhance reinforcement learning in the mammalian hippocampal-striatal system rather than provide its core mechanism.


title: Transfer function interpolation remainder formula of rational Krylov subspace methods

abstract: Rational Krylov subspace projection methods are one of successful methods in MOR, mainly because some order derivatives of the approximate and original transfer functions are the same. This is the well known moments matching result. However, the properties of points which are far from the interpolating points are little known. In this paper, we obtain the error's explicit expression which involves shifts and Ritz values. The advantage of our result over than the known moments matches theory is, to some extent, similar to the one of Lagrange type remainder formula over than Peano Type remainder formula in Taylor theorem. Expect for the proof, we also provide three explanations for the error formula. One explanation shows that in the Gauss-Christoffel quadrature sense, the error is the Gauss quadrature remainder, when the Gauss quadrature formula is applied onto the resolvent function. By using the error formula, we propose some greedy algorithms for the interpolatory $H_{\infty}$ norm MOR.


title: Event Camera Based Real-Time Detection and Tracking of Indoor Ground Robots

abstract: This paper presents a real-time method to detect and track multiple mobile ground robots using event cameras. The method uses density-based spatial clustering of applications with noise (DBSCAN) to detect the robots and a single k-dimensional ($k - d$) tree to accurately keep track of them as they move in an indoor arena. Robust detections and tracks are maintained in the face of event camera noise and lack of events (due to robots moving slowly or stopping). An off-the-shelf RGB camera-based tracking system was used to provide ground truth. Experiments including up to 4 robots are performed to study the effect of i) varying DBSCAN parameters, ii) the event accumulation time, iii) the number of robots in the arena, iv) the speed of the robots, and v) variation in ambient light conditions on the detection and tracking performance. The experimental results showed 100% detection and tracking fidelity in the face of event camera noise and robots stopping for tests involving up to 3 robots (and upwards of 93% for 4 robots). When the lighting conditions were varied, a graceful degradation in detection and tracking fidelity was observed.


title: The Sensitivity of Word Embeddings-based Author Detection Models to Semantic-preserving Adversarial Perturbations

abstract: Authorship analysis is an important subject in the field of natural language processing. It allows the detection of the most likely writer of articles, news, books, or messages. This technique has multiple uses in tasks related to authorship attribution, detection of plagiarism, style analysis, sources of misinformation, etc. The focus of this paper is to explore the limitations and sensitiveness of established approaches to adversarial manipulations of inputs. To this end, and using those established techniques, we first developed an experimental frame-work for author detection and input perturbations. Next, we experimentally evaluated the performance of the authorship detection model to a collection of semantic-preserving adversarial perturbations of input narratives. Finally, we compare and analyze the effects of different perturbation strategies, input and model configurations, and the effects of these on the author detection model.


title: A Data-Driven Convex Programming Approach to Worst-Case Robust Tracking Controller Design

abstract: This paper studies finite-horizon robust tracking control for discrete-time linear systems, based on input-output data. We leverage behavioral theory to represent system trajectories through a set of noiseless historical data, instead of using an explicit system model. By assuming that recent output data available to the controller are affected by noise terms verifying a quadratic bound, we formulate an optimization problem with a linear cost and LMI constraints for solving the robust tracking problem without any approximations. Our approach hinges on a parameterization of noise trajectories compatible with the data-dependent system representation and on a reformulation of the tracking cost, which enables the application of the S-lemma. In addition, we propose a method for reducing the computational complexity and demonstrate that the size of the resulting LMIs does not scale with the number of historical data. Finally, we show that the proposed formulation can easily incorporate actuator disturbances as well as constraints on inputs and outputs. The performance of the new controllers is discussed through simulations.


title: {\Lambda} polarization from thermalized jet energy

abstract: We examine the formation of vortical "smoke rings" as a result of thermalization of energy lost by a jet. We simulate the formation and evolution of these rings using hydrodynamics and define an observable that allows to probe this phenomenon experimentally. We argue that observation of vorticity associated with jets would be an experimental confirmation of the thermalization of the energy lost by quenched jets, and also a probe of shear viscosity.


title: Dynamic Games among Teams with Delayed Intra-Team Information Sharing

abstract: We analyze a class of stochastic dynamic games among teams with asymmetric information, where members of a team share their observations internally with a delay of $d$. Each team is associated with a controlled Markov Chain, whose dynamics are coupled through the players' actions. These games exhibit challenges in both theory and practice due to the presence of signaling and the increasing domain of information over time. We develop a general approach to characterize a subset of Nash Equilibria where the agents can use a compressed version of their information, instead of the full information, to choose their actions. We identify two subclasses of strategies: Sufficient Private Information Based (SPIB) strategies, which only compress private information, and Compressed Information Based (CIB) strategies, which compress both common and private information. We show that while SPIB-strategy-based equilibria always exist, the same is not true for CIB-strategy-based equilibria. We develop a backward inductive sequential procedure, whose solution (if it exists) provides a CIB strategy-based equilibrium. We identify some instances where we can guarantee the existence of a solution to the above procedure. Our results highlight the tension among compression of information, existence of (compression based) equilibria, and backward inductive sequential computation of such equilibria in stochastic dynamic games with asymmetric information.


title: Nonthermal processes in hot accretion flows onto supermassive black holes: An inhomogeneous model

abstract: Context. Many low-redshift active galactic nuclei harbor a supermassive black hole accreting matter at low or medium rates. At such rates, the accretion flow usually consists of a cold optically thick disk, plus a hot, low density, collisionless corona. In the latter component, charged particles can be accelerated to high energies by various mechanisms. Aims. We aim to investigate, in detail, nonthermal processes in hot accretion flows onto supermassive black holes, covering a wide range of accretion rates and luminosities. Methods. We developed a model consisting of a thin Shakura-Sunyaev disk plus an inner hot accretion flow or corona, modeled as a radiatively inefficient accretion flow, where nonthermal processes take place. We solved the transport equations for relativistic particles and estimated the spectral energy distributions resulting from nonthermal interactions between the various particle species and the fields in the source. Results. We covered a variety of scenarios, from low accretion rates up to 10% of the Eddington limit, and identified the relevant cooling mechanisms in each case. The presence of hadrons in the hot flow is decisive for the spectral shape, giving rise to secondary particles and gamma-ray cascades. We applied our model to the source IC 4329A, confirming earlier results that showed evidence of nonthermal particles in the corona.


title: Dynamic Graph Modeling of Simultaneous EEG and Eye-tracking Data for Reading Task Identification

abstract: We present a new approach, that we call AdaGTCN, for identifying human reader intent from Electroencephalogram~(EEG) and Eye movement~(EM) data in order to help differentiate between normal reading and task-oriented reading. Understanding the physiological aspects of the reading process~(the cognitive load and the reading intent) can help improve the quality of crowd-sourced annotated data. Our method, Adaptive Graph Temporal Convolution Network (AdaGTCN), uses an Adaptive Graph Learning Layer and Deep Neighborhood Graph Convolution Layer for identifying the reading activities using time-locked EEG sequences recorded during word-level eye-movement fixations. Adaptive Graph Learning Layer dynamically learns the spatial correlations between the EEG electrode signals while the Deep Neighborhood Graph Convolution Layer exploits temporal features from a dense graph neighborhood to establish the state of the art in reading task identification over other contemporary approaches. We compare our approach with several baselines to report an improvement of 6.29% on the ZuCo 2.0 dataset, along with extensive ablation experiments


title: KAM Theory Meets Statistical Learning Theory: Hamiltonian Neural Networks with Non-Zero Training Loss

abstract: Many physical phenomena are described by Hamiltonian mechanics using an energy function (the Hamiltonian). Recently, the Hamiltonian neural network, which approximates the Hamiltonian as a neural network, and its extensions have attracted much attention. This is a very powerful method, but its use in theoretical studies remains limited. In this study, by combining the statistical learning theory and Kolmogorov-Arnold-Moser (KAM) theory, we provide a theoretical analysis of the behavior of Hamiltonian neural networks when the learning error is not completely zero. A Hamiltonian neural network with non-zero errors can be considered as a perturbation from the true dynamics, and the perturbation theory of the Hamilton equation is widely known as the KAM theory. To apply the KAM theory, we provide a generalization error bound for Hamiltonian neural networks by deriving an estimate of the covering number of the gradient of the multi-layer perceptron, which is the key ingredient of the model. This error bound gives an $L^\infty$ bound on the Hamiltonian that is required in the application of the KAM theory.


title: Finite Confluences and Closed Pattern Mining

abstract: The purpose of this article is to propose and investigate a partial order structure weaker than the lattice structure and which have nice properties regarding closure operators. We extend accordingly closed pattern mining and formal concept analysis to such structures we further call confluences. The primary motivation for investigating these structures is that it allows to reduce a lattice to a part whose elements are connected, as in some graph, still preserving a useful characterization of closure operators. Our investigation also considers how reducing one of the lattice involved in a Galois connection affects the structure of the closure operators ranges. When extending this way formal concept analysis we will focus on the intensional space, i.e. in reducing the pattern language, while recent investigations rather explored the reduction of the extensional space to connected elements.


title: Chasm in Hegemony: Explaining and Reproducing Disparities in Homophilous Networks

abstract: In networks with a minority and a majority community, it is well-studied that minorities are under-represented at the top of the social hierarchy. However, researchers are less clear about the representation of minorities from the lower levels of the hierarchy, where other disadvantages or vulnerabilities may exist. We offer a more complete picture of social disparities at each social level with empirical evidence that the minority representation exhibits two opposite phases: at the higher rungs of the social ladder, the representation of the minority community decreases; but, lower in the ladder, which is more populous, as you ascend, the representation of the minority community improves. We refer to this opposing phenomenon between the upper-level and lower-level as the \emph{chasm effect}. Previous models of network growth with homophily fail to detect and explain the presence of this chasm effect. We analyze the interactions among a few well-observed network-growing mechanisms with a simple model to reveal the sufficient and necessary conditions for both phases in the chasm effect to occur. By generalizing the simple model naturally, we present a complete bi-affiliation bipartite network-growth model that could successfully capture disparities at all social levels and reproduce real social networks. Finally, we illustrate that addressing the chasm effect can create fairer systems with two applications in advertisement and fact-checks, thereby demonstrating the potential impact of the chasm effect on the future research of minority-majority disparities and fair algorithms.


title: Estimating Average Treatment Effects with Support Vector Machines

abstract: Support vector machine (SVM) is one of the most popular classification algorithms in the machine learning literature. We demonstrate that SVM can be used to balance covariates and estimate average causal effects under the unconfoundedness assumption. Specifically, we adapt the SVM classifier as a kernel-based weighting procedure that minimizes the maximum mean discrepancy between the treatment and control groups while simultaneously maximizing effective sample size. We also show that SVM is a continuous relaxation of the quadratic integer program for computing the largest balanced subset, establishing its direct relation to the cardinality matching method. Another important feature of SVM is that the regularization parameter controls the trade-off between covariate balance and effective sample size. As a result, the existing SVM path algorithm can be used to compute the balance-sample size frontier. We characterize the bias of causal effect estimation arising from this trade-off, connecting the proposed SVM procedure to the existing kernel balancing methods. Finally, we conduct simulation and empirical studies to evaluate the performance of the proposed methodology and find that SVM is competitive with the state-of-the-art covariate balancing methods.


title: Forcing the escape: Partial control of escaping orbits from a transient chaotic region

abstract: A new control algorithm based on the partial control method has been developed. The general situation we are considering is an orbit starting in a certain phase space region Q having a chaotic transient behavior affected by noise, so that the orbit will definitely escape from Q in an unpredictable number of iterations. Thus, the goal of the algorithm is to control in a predictable manner when to escape. While partial control has been used as a way to avoid escapes, here we want to adapt it to force the escape in a controlled manner. We have introduced new tools such as escape functions and escape sets that once computed makes the control of the orbit straightforward. We have applied the new idea to three different cases in order to illustrate the various application possibilities of this new algorithm.


title: Investigating Moral Foundations from Web Trending Topics

abstract: Moral foundations theory helps understand differences in morality across cultures. In this paper, we propose a model to predict moral foundations (MF) from social media trending topics. We also investigate whether differences in MF influence emotional traits. Our results are promising and leave room for future research avenues.


title: PolicySpace2: modeling markets and endogenous public policies

abstract: Policymakers decide on alternative policies facing restricted budgets and uncertain, ever-changing future. Designing public policies is further difficult due to the need to decide on priorities and handle effects across policies. Housing policies, specifically, involve heterogeneous characteristics of properties themselves and the intricacy of housing markets and the spatial context of cities. We propose PolicySpace2 (PS2) as an adapted and extended version of the open source PolicySpace agent-based model. PS2 is a computer simulation that relies on empirically detailed spatial data to model real estate, along with labor, credit, and goods and services markets. Interaction among workers, firms, a bank, households and municipalities follow the literature benchmarks to integrate economic, spatial and transport scholarship. PS2 is applied to a comparison among three competing public policies aimed at reducing inequality and alleviating poverty: (a) house acquisition by the government and distribution to lower income households, (b) rental vouchers, and (c) monetary aid. Within the model context, the monetary aid, that is, smaller amounts of help for a larger number of households, makes the economy perform better in terms of production, consumption, reduction of inequality, and maintenance of financial duties. PS2 as such is also a framework that may be further adapted to a number of related research questions.


title: Coherent Suppression of Molecular Bremsstrahlung Radiation at GHz Frequencies in the Ionization Trail of Extensive Air Showers

abstract: Several attempts to detect extensive air showers (EAS) induced by ultrahigh-energy cosmic rays have been conducted in the last decade based on the molecular Bremsstrahlung radiation (MBR) at GHz frequencies from quasi-elastic collisions of ionisation electrons left in the atmosphere after the passage of the cascade of particles. These attempts have led to the detection of a handful of signals only, all of them forward-directed along the shower axis and hence suggestive of originating from geomagnetic and Askaryan emissions extending into GHz frequencies close to the Cherenkov angle. In this paper, the lack of detection of events is explained by the coherent suppression of the MBR in frequency ranges below the collision rate due to the destructive interferences impacting the emission amplitude of photons between the successive collisions of the electrons. The spectral intensity at the ground level is shown to be several orders of magnitude below the sensitivity of experimental setups. In particular, the spectral intensity at 10~km from the shower core for a vertical shower induced by a proton of $10^{17.5}$ eV is 7-to-8 orders of magnitude below the reference value anticipated from a scaling law converting a laboratory measurement to EAS expectations. Consequently, the MBR cannot be seen as the basis of a new detection technique of EAS for the next decades.


title: Boosting background suppression in the NEXT experiment through Richardson-Lucy deconvolution

abstract: Next-generation neutrinoless double beta decay experiments aim for half-life sensitivities of ~$10^{27}$ yr, requiring suppressing backgrounds to <1 count/tonne/yr. For this, any extra background rejection handle, beyond excellent energy resolution and the use of extremely radiopure materials, is of utmost importance. The NEXT experiment exploits differences in the spatial ionization patterns of double beta decay and single-electron events to discriminate signal from background. While the former display two Bragg peak dense ionization regions at the opposite ends of the track, the latter typically have only one such feature. Thus, comparing the energies at the track extremes provides an additional rejection tool. The unique combination of the topology-based background discrimination and excellent energy resolution (1% FWHM at the Q-value of the decay) is the distinguishing feature of NEXT. Previous studies demonstrated a topological background rejection factor of ~5 when reconstructing electron-positron pairs in the $^{208}$Tl 1.6 MeV double escape peak (with Compton events as background), recorded in the NEXT-White demonstrator at the Laboratorio Subterr\'aneo de Canfranc, with 72% signal efficiency. This was recently improved through the use of a deep convolutional neural network to yield a background rejection factor of ~10 with 65% signal efficiency. Here, we present a new reconstruction method, based on the Richardson-Lucy deconvolution algorithm, which allows reversing the blurring induced by electron diffusion and electroluminescence light production in the NEXT TPC. The new method yields highly refined 3D images of reconstructed events, and, as a result, significantly improves the topological background discrimination. When applied to real-data 1.6 MeV $e^-e^+$ pairs, it leads to a background rejection factor of 27 at 57% signal efficiency.


title: On Meritocracy in Optimal Set Selection

abstract: Typically, merit is defined with respect to some intrinsic measure of worth. We instead consider a setting where an individual's worth is \emph{relative}: when a Decision Maker (DM) selects a set of individuals from a population to maximise expected utility, it is natural to consider the \emph{Expected Marginal Contribution} (EMC) of each person to the utility. We show that this notion satisfies an axiomatic definition of fairness for this setting. We also show that for certain policy structures, this notion of fairness is aligned with maximising expected utility, while for linear utility functions it is identical to the Shapley value. However, for certain natural policies, such as those that select individuals with a specific set of attributes (e.g. high enough test scores for college admissions), there is a trade-off between meritocracy and utility maximisation. We analyse the effect of constraints on the policy on both utility and fairness in extensive experiments based on college admissions and outcomes in Norwegian universities.


title: Polytypism and Superconductivity in the NbS$_2$ System

abstract: We report on the phase formation and the superconducting properties in the NbS$_2$ system. Specifically, we have performed a series of standardized solid-state syntheses in this system, which allow us to establish a comprehensive synthesis map for the formation of the two polytypes 2H-NbS$_2$ and 3R-NbS$_2$, respectively. We show that the identification of two polytypes by means of X-ray diffraction is not always unambiguous. Our physical property measurements on a phase-pure sample of 3R-NbS$_2$, on a phase-pure sample of 2H-NbS$_2$, and a mixed phase sample confirm earlier reports that 2H-NbS$_2$ is a bulk superconductor and that 3R-NbS$_2$ is not a superconductor above $T =$ 1.75 K. Our results clearly show that specific heat measurements, as true bulk measurements, are crucial for the identification of superconducting materials in this and related systems. Our results indicate that for the investigation of van-der-Waals materials great care has to be taken on choosing the synthesis conditions for obtaining phase pure samples.


title: Feature Importance Explanations for Temporal Black-Box Models

abstract: Models in the supervised learning framework may capture rich and complex representations over the features that are hard for humans to interpret. Existing methods to explain such models are often specific to architectures and data where the features do not have a time-varying component. In this work, we propose TIME, a method to explain models that are inherently temporal in nature. Our approach (i) uses a model-agnostic permutation-based approach to analyze global feature importance, (ii) identifies the importance of salient features with respect to their temporal ordering as well as localized windows of influence, and (iii) uses hypothesis testing to provide statistical rigor.


title: Non-Singular Adversarial Robustness of Neural Networks

abstract: Adversarial robustness has become an emerging challenge for neural network owing to its over-sensitivity to small input perturbations. While being critical, we argue that solving this singular issue alone fails to provide a comprehensive robustness assessment. Even worse, the conclusions drawn from singular robustness may give a false sense of overall model robustness. Specifically, our findings show that adversarially trained models that are robust to input perturbations are still (or even more) vulnerable to weight perturbations when compared to standard models. In this paper, we formalize the notion of non-singular adversarial robustness for neural networks through the lens of joint perturbations to data inputs as well as model weights. To our best knowledge, this study is the first work considering simultaneous input-weight adversarial perturbations. Based on a multi-layer feed-forward neural network model with ReLU activation functions and standard classification loss, we establish error analysis for quantifying the loss sensitivity subject to $\ell_\infty$-norm bounded perturbations on data inputs and model weights. Based on the error analysis, we propose novel regularization functions for robust training and demonstrate improved non-singular robustness against joint input-weight adversarial perturbations.


title: AGN and star formation across cosmic time

abstract: We investigate the balance of power between stars and AGN across cosmic history, based on the comparison between the infrared (IR) galaxy luminosity function (LF) and the IR AGN LF. The former corresponds to emission from dust heated by stars and AGN, whereas the latter includes emission from AGN-heated dust only. We find that at all redshifts (at least up to z~2.5), the high luminosity tails of the two LFs converge, indicating that the most infrared-luminous galaxies are AGN-powered. Our results shed light to the decades-old conundrum regarding the flatter high-luminosity slope seen in the IR galaxy LF compared to that in the UV and optical. We attribute this difference to the increasing fraction of AGN-dominated galaxies with increasing total infrared luminosity (L_IR). We partition the L_IR-z parameter space into a star-formation and an AGN-dominated region, finding that the most luminous galaxies at all epochs lie in the AGN-dominated region. This sets a potential `limit' to attainable star formation rates, casting doubt on the abundance of `extreme starbursts': if AGN did not exist, L_IR>10^13 Lsun galaxies would be significantly rarer than they currently are in our observable Universe. We also find that AGN affect the average dust temperatures (T_dust) of galaxies and hence the shape of the well-known L_IR-T_dust relation. We propose that the reason why local ULIRGs are hotter than their high redshift counterparts is because of a higher fraction of AGN-dominated galaxies amongst the former group.


title: On the number of hyperedges in the hypergraph of lines and pseudo-discs

abstract: Consider the hypergraph whose vertex set is a family of $n$ lines in general position in the plane, and whose hyperedges are induced by intersections with a family of pseudo-discs. We prove that the number of $t$-hyperedges is bounded by $O_t(n^2)$ and that the total number of hyperedges is bounded by $O(n^3)$. Both bounds are tight.


title: Baby Intuitions Benchmark (BIB): Discerning the goals, preferences, and actions of others

abstract: To achieve human-like common sense about everyday life, machine learning systems must understand and reason about the goals, preferences, and actions of other agents in the environment. By the end of their first year of life, human infants intuitively achieve such common sense, and these cognitive achievements lay the foundation for humans' rich and complex understanding of the mental states of others. Can machines achieve generalizable, commonsense reasoning about other agents like human infants? The Baby Intuitions Benchmark (BIB) challenges machines to predict the plausibility of an agent's behavior based on the underlying causes of its actions. Because BIB's content and paradigm are adopted from developmental cognitive science, BIB allows for direct comparison between human and machine performance. Nevertheless, recently proposed, deep-learning-based agency reasoning models fail to show infant-like reasoning, leaving BIB an open challenge.


title: High order strong stability preserving multi-derivative implicit and IMEX Runge--Kutta methods with asymptotic preserving properties

abstract: In this work we present a class of high order unconditionally strong stability preserving (SSP) implicit multi-derivative Runge--Kutta schemes, and SSP implicit-explicit (IMEX) multi-derivative Runge--Kutta schemes where the time-step restriction is independent of the stiff term. The unconditional SSP property for a method of order $p>2$ is unique among SSP methods, and depends on a backward-in-time assumption on the derivative of the operator. We show that this backward derivative condition is satisfied in many relevant cases where SSP IMEX schemes are desired. We devise unconditionally SSP implicit Runge--Kutta schemes of order up to $p=4$, and IMEX Runge--Kutta schemes of order up to $p=3$. For the multi-derivative IMEX schemes, we also derive and present the order conditions, which have not appeared previously. The unconditional SSP condition ensures that these methods are positivity preserving, and we present sufficient conditions under which such methods are also asymptotic preserving when applied to a range of problems, including a hyperbolic relaxation system, the Broadwell model, and the Bhatnagar-Gross-Krook (BGK) kinetic equation. We present numerical results to support the theoretical results, on a variety of problems.


title: Geometric invariant decomposition of SU(3)

abstract: A novel invariant decomposition of diagonalizable $n \times n$ matrices into $n$ commuting matrices is presented. This decomposition is subsequently used to split the fundamental representation of $\mathfrak{su}(3)$ Lie algebra elements into at most three commuting elements of $\mathfrak{u}(3)$. As a result, the exponential of an $\mathfrak{su}(3)$ Lie algebra element can be split into three commuting generalized Euler's formulas, or conversely, a Lie group element can be factorized into at most three generalized Euler's formulas. After the factorization has been performed, the logarithm follows immediately.


title: State Augmented Constrained Reinforcement Learning: Overcoming the Limitations of Learning with Rewards

abstract: A common formulation of constrained reinforcement learning involves multiple rewards that must individually accumulate to given thresholds. In this class of problems, we show a simple example in which the desired optimal policy cannot be induced by any weighted linear combination of rewards. Hence, there exist constrained reinforcement learning problems for which neither regularized nor classical primal-dual methods yield optimal policies. This work addresses this shortcoming by augmenting the state with Lagrange multipliers and reinterpreting primal-dual methods as the portion of the dynamics that drives the multipliers evolution. This approach provides a systematic state augmentation procedure that is guaranteed to solve reinforcement learning problems with constraints. Thus, as we illustrate by an example, while previous methods can fail at finding optimal policies, running the dual dynamics while executing the augmented policy yields an algorithm that provably samples actions from the optimal policy.


title: Multi-Feature Multi-Scale CNN-Derived COVID-19 Classification from Lung Ultrasound Data

abstract: The global pandemic of the novel coronavirus disease 2019 (COVID-19) has put tremendous pressure on the medical system. Imaging plays a complementary role in the management of patients with COVID-19. Computed tomography (CT) and chest X-ray (CXR) are the two dominant screening tools. However, difficulty in eliminating the risk of disease transmission, radiation exposure and not being costeffective are some of the challenges for CT and CXR imaging. This fact induces the implementation of lung ultrasound (LUS) for evaluating COVID-19 due to its practical advantages of noninvasiveness, repeatability, and sensitive bedside property. In this paper, we utilize a deep learning model to perform the classification of COVID-19 from LUS data, which could produce objective diagnostic information for clinicians. Specifically, all LUS images are processed to obtain their corresponding local phase filtered images and radial symmetry transformed images before fed into the multi-scale residual convolutional neural network (CNN). Secondly, image combination as the input of the network is used to explore rich and reliable features. Feature fusion strategy at different levels is adopted to investigate the relationship between the depth of feature aggregation and the classification accuracy. Our proposed method is evaluated on the point-of-care US (POCUS) dataset together with the Italian COVID-19 Lung US database (ICLUS-DB) and shows promising performance for COVID-19 prediction.


title: Quantum disordered state of magnetic charges in nanoengineered honeycomb lattice

abstract: A quantum magnetic state due to magnetic charges is never observed, even though they are treated as quantum mechanical variable in theoretical calculations. Here, we demonstrate the occurrence of a novel quantum disordered state of magnetic charges in nanoengineered magnetic honeycomb lattice of ultra-small connecting elements. The experimental research, performed using spin resolved neutron scattering, reveals a massively degenerate ground state, comprised of low integer and energetically forbidden high integer magnetic charges, that manifests cooperative paramagnetism at low temperature. The system tends to preserve the degenerate configuration even under large magnetic field application. It exemplifies the robustness of disordered correlation of magnetic charges in 2D honeycomb lattice. The realization of quantum disordered ground state elucidates the dominance of exchange energy, which is enabled due to the nanoscopic magnetic element size in nanoengineered honeycomb. Consequently, an archetypal platform is envisaged to study quantum mechanical phenomena due to emergent magnetic charges.


title: Arguments for the Unsuitability of Convolutional Neural Networks for Non--Local Tasks

abstract: Convolutional neural networks have established themselves over the past years as the state of the art method for image classification, and for many datasets, they even surpass humans in categorizing images. Unfortunately, the same architectures perform much worse when they have to compare parts of an image to each other to correctly classify this image. Until now, no well-formed theoretical argument has been presented to explain this deficiency. In this paper, we will argue that convolutional layers are of little use for such problems, since comparison tasks are global by nature, but convolutional layers are local by design. We will use this insight to reformulate a comparison task into a sorting task and use findings on sorting networks to propose a lower bound for the number of parameters a neural network needs to solve comparison tasks in a generalizable way. We will use this lower bound to argue that attention, as well as iterative/recurrent processing, is needed to prevent a combinatorial explosion.


title: Quantum Algorithm for Quantum State Discrimination via Partial Negation and Weak Measurement

abstract: The quantum state discrimination problem is to distinguish between non-orthogonal quantum states. This problem has many applications in quantum information theory, quantum communication and quantum cryptography. In this paper a quantum algorithm using weak measurement and partial negation will be proposed to solve the quantum state discrimination problem using a single copy of an unknown qubit. The usage of weak measurement makes it possible to reconstruct the qubit after measurement since the superposition will not be destroyed due to measurement. The proposed algorithm will be able to determine, with high probability of success, the state of the unknown qubit and whether it is encoded in the Hadamard or the computational basis by counting the outcome of the successive measurements on an auxiliary qubit.


title: Conditions for Exact Convex Relaxation and No Spurious Local Optima

abstract: Non-convex optimization problems can be approximately solved via relaxation or local algorithms. For many practical problems such as optimal power flow (OPF) problems, both approaches tend to succeed in the sense that relaxation is usually exact and local algorithms usually converge to a global optimum. In this paper, we study conditions which are sufficient or necessary for such non-convex problems to simultaneously have exact relaxation and no spurious local optima. Those conditions help us explain the widespread empirical experience that local algorithms for OPF problems often work extremely well.


title: Multi-Group Multicast Beamforming by Superiorized Projections onto Convex Sets

abstract: In this paper, we propose an iterative algorithm to address the nonconvex multi-group multicast beamforming problem with quality-of-service constraints and per-antenna power constraints. We formulate a convex relaxation of the problem as a semidefinite program in a real Hilbert space, which allows us to approximate a point in the feasible set by iteratively applying a bounded perturbation resilient fixed-point mapping. Inspired by the superiorization methodology, we use this mapping as a basic algorithm, and we add in each iteration a small perturbation with the intent to reduce the objective value and the distance to nonconvex rank-constraint sets. We prove that the sequence of perturbations is bounded, so the algorithm is guaranteed to converge to a feasible point of the relaxed semidefinite program. Simulations show that the proposed approach outperforms existing algorithms in terms of both computation time and approximation gap in many cases.


title: Inferring the Type of Phase Transitions Undergone in Epileptic Seizures Using Random Graph Hidden Markov Models for Percolation in Noisy Dynamic Networks

abstract: In clinical neuroscience, epileptic seizures have been associated with the sudden emergence of coupled activity across the brain. The resulting functional networks - in which edges indicate strong enough coupling between brain regions - are consistent with the notion of percolation, which is a phenomenon in complex networks corresponding to the sudden emergence of a giant connected component. Traditionally, work has concentrated on noise-free percolation with a monotonic process of network growth, but real-world networks are more complex. We develop a class of random graph hidden Markov models (RG-HMMs) for characterizing percolation regimes in noisy, dynamically evolving networks in the presence of edge birth and edge death, as well as noise. This class is used to understand the type of phase transitions undergone in a seizure, and in particular, distinguishing between different percolation regimes in epileptic seizures. We develop a hypothesis testing framework for inferring putative percolation mechanisms. As a necessary precursor, we present an EM algorithm for estimating parameters from a sequence of noisy networks only observed at a longitudinal subsampling of time points. Our results suggest that different types of percolation can occur in human seizures. The type inferred may suggest tailored treatment strategies and provide new insights into the fundamental science of epilepsy.


title: Magnetic charge's relaxation propelled electricity in two-dimensional magnetic honeycomb lattice

abstract: Emerging new concepts, such as magnetic charge dynamics in two-dimensional magnetic material, can provide novel mechanism for spin based electrical transport at macroscopic length. In artificial spin ice of single domain elements, magnetic charge's relaxation can create an efficient electrical pathway for conduction by generating fluctuations in local magnetic field that couple with conduction electrons spins. In a first demonstration, we show that the electrical conductivity is propelled by more than an order of magnitude at room temperature due to magnetic charge defects sub-picosecond relaxation in artificial magnetic honeycomb lattice. The direct evidence to the proposed electrical conduction mechanism in two-dimensional frustrated magnet points to the untapped potential for spintronic applications in this system.


title: Structure and hardness of in situ synthesized nano-oxide strengthened CoCrFeNi high entropy alloy thin films

abstract: In this study, we report on face-centered cubic structured CoCrFeNi high-entropy alloy thin films with finely dispersed nano-oxide particles which are formed by internal oxidation. Analytical scanning transmission electron microscopy imaging found that the particles are Cr2O3. The oxide particles contribute to the hardening of the film increasing its hardness by 14% compared to that of the film without precipitates, through the Orowan-type strengthening mechanism. Our novel approach paves the way to design medium- and high-entropy alloys with high strength by making use of oxide phases.


title: Operator preconditioning: the simplest case

abstract: Using the framework of operator or Calder\'on preconditioning, uniform preconditioners are constructed for elliptic operators discretized with continuous finite (or boundary) elements. The preconditioners are constructed as the composition of an opposite order operator, discretized on the same ansatz space, and two diagonal scaling operators.


title: Learning to Drop Points for LiDAR Scan Synthesis

abstract: 3D laser scanning by LiDAR sensors plays an important role for mobile robots to understand their surroundings. Nevertheless, not all systems have high resolution and accuracy due to hardware limitations, weather conditions, and so on. Generative modeling of LiDAR data as scene priors is one of the promising solutions to compensate for unreliable or incomplete observations. In this paper, we propose a novel generative model for learning LiDAR data based on generative adversarial networks. As in the related studies, we process LiDAR data as a compact yet lossless representation, a cylindrical depth map. However, despite the smoothness of real-world objects, many points on the depth map are dropped out through the laser measurement, which causes learning difficulty on generative models. To circumvent this issue, we introduce measurement uncertainty into the generation process, which allows the model to learn a disentangled representation of the underlying shape and the dropout noises from a collection of real LiDAR data. To simulate the lossy measurement, we adopt a differentiable sampling framework to drop points based on the learned uncertainty. We demonstrate the effectiveness of our method on synthesis and reconstruction tasks using two datasets. We further showcase potential applications by restoring LiDAR data with various types of corruption.


title: The subspace $c_0$ is not complemented in $ac_0$

abstract: We prove that the subspace $c_0$ of sequences that converge to zero is not complemented in the space $ac_0$ of sequences that almost converge to zero. We proceed with applying the same approach to inclusion chain $c_0\subset A_0 \subset \ell_\infty$.


title: Radar Cross Section Based Statistical Recognition of UAVs at Microwave Frequencies

abstract: This paper presents a radar cross-section (RCS)-based statistical recognition system for identifying/ classifying unmanned aerial vehicles (UAVs) at microwave frequencies. First, the paper presents the results of the vertical (VV) and horizontal (HH) polarization RCS measurement of six commercial UAVs at 15 GHz and 25 GHz in a compact range anechoic chamber. The measurement results show that the average RCS of the UAVs depends on shape, size, material composition of the target UAV as well as the azimuth angle, frequency, and polarization of the illuminating radar. Afterward, radar characterization of the target UAVs is achieved by fitting the RCS measurement data to 11 different statistical models. From the model selection analysis, we observe that the lognormal, generalized extreme value, and gamma distributions are most suitable for modeling the RCS of the commercial UAVs while the Gaussian distribution performed relatively poorly. The best UAV radar statistics forms the class conditional probability densities for the proposed UAV statistical recognition system. The performance of the UAV statistical recognition system is evaluated at different signal noise ratio (SNR) with the aid of Monte Carlo analysis. At an SNR of 10 dB, the average classification accuracy of 97.43% or better is achievable.


title: Location Trace Privacy Under Conditional Priors

abstract: Providing meaningful privacy to users of location based services is particularly challenging when multiple locations are revealed in a short period of time. This is primarily due to the tremendous degree of dependence that can be anticipated between points. We propose a R\'enyi divergence based privacy framework for bounding expected privacy loss for conditionally dependent data. Additionally, we demonstrate an algorithm for achieving this privacy under Gaussian process conditional priors. This framework both exemplifies why conditionally dependent data is so challenging to protect and offers a strategy for preserving privacy to within a fixed radius for sensitive locations in a user's trace.


title: Integrated photonics on thin-film lithium niobate

abstract: Lithium niobate (LN), an outstanding and versatile material, has influenced our daily life for decades: from enabling high-speed optical communications that form the backbone of the Internet to realizing radio-frequency filtering used in our cell phones. This half-century-old material is currently embracing a revolution in thin-film LN integrated photonics. The success of manufacturing wafer-scale, high-quality, thin films of LN on insulator (LNOI), accompanied with breakthroughs in nanofabrication techniques, have made high-performance integrated nanophotonic components possible. With rapid development in the past few years, some of these thin-film LN devices, such as optical modulators and nonlinear wavelength converters, have already outperformed their legacy counterparts realized in bulk LN crystals. Furthermore, the nanophotonic integration enabled ultra-low-loss resonators in LN, which unlocked many novel applications such as optical frequency combs and quantum transducers. In this Review, we cover -- from basic principles to the state of the art -- the diverse aspects of integrated thin-film LN photonics, including the materials, basic passive components, and various active devices based on electro-optics, all-optical nonlinearities, and acousto-optics. We also identify challenges that this platform is currently facing and point out future opportunities. The field of integrated LNOI photonics is advancing rapidly and poised to make critical impacts on a broad range of applications in communication, signal processing, and quantum information.


title: Quantifying Confounding Bias in Generative Art: A Case Study

abstract: In recent years, AI generated art has become very popular. From generating art works in the style of famous artists like Paul Cezanne and Claude Monet to simulating styles of art movements like Ukiyo-e, a variety of creative applications have been explored using AI. Looking from an art historical perspective, these applications raise some ethical questions. Can AI model artists' styles without stereotyping them? Does AI do justice to the socio-cultural nuances of art movements? In this work, we take a first step towards analyzing these issues. Leveraging directed acyclic graphs to represent potential process of art creation, we propose a simple metric to quantify confounding bias due to the lack of modeling the influence of art movements in learning artists' styles. As a case study, we consider the popular cycleGAN model and analyze confounding bias across various genres. The proposed metric is more effective than state-of-the-art outlier detection method in understanding the influence of art movements in artworks. We hope our work will elucidate important shortcomings of computationally modeling artists' styles and trigger discussions related to accountability of AI generated art.


title: The SpaceNet Multi-Temporal Urban Development Challenge

abstract: Building footprints provide a useful proxy for a great many humanitarian applications. For example, building footprints are useful for high fidelity population estimates, and quantifying population statistics is fundamental to ~1/4 of the United Nations Sustainable Development Goals Indicators. In this paper we (the SpaceNet Partners) discuss efforts to develop techniques for precise building footprint localization, tracking, and change detection via the SpaceNet Multi-Temporal Urban Development Challenge (also known as SpaceNet 7). In this NeurIPS 2020 competition, participants were asked identify and track buildings in satellite imagery time series collected over rapidly urbanizing areas. The competition centered around a brand new open source dataset of Planet Labs satellite imagery mosaics at 4m resolution, which includes 24 images (one per month) covering ~100 unique geographies. Tracking individual buildings at this resolution is quite challenging, yet the winning participants demonstrated impressive performance with the newly developed SpaceNet Change and Object Tracking (SCOT) metric. This paper details the top-5 winning approaches, as well as analysis of results that yielded a handful of interesting anecdotes such as decreasing performance with latitude.


title: Exact site frequency spectra of neutrally evolving tumors: a transition between power laws reveals a signature of cell viability

abstract: The site frequency spectrum (SFS) is a popular summary statistic of genomic data. While the SFS of a constant-sized population undergoing neutral mutations has been extensively studied in population genetics, the rapidly growing amount of cancer genomic data has attracted interest in the spectrum of an exponentially growing population. Recent theoretical results have generally dealt with special or limiting cases, such as considering only cells with an infinite line of descent, assuming deterministic tumor growth, or taking large-time or large-population limits. In this work, we derive exact expressions for the expected SFS of a cell population that evolves according to a stochastic branching process, first for cells with an infinite line of descent and then for the total population, evaluated either at a fixed time (fixed-time spectrum) or at the stochastic time at which the population reaches a certain size (fixed-size spectrum). We find that while the rate of mutation scales the SFS of the total population linearly, the rates of cell birth and cell death change the shape of the spectrum at the small-frequency end, inducing a transition between a $1/j^2$ power-law spectrum and a $1/j$ spectrum as cell viability decreases. We show that this insight can in principle be used to estimate the ratio between the rate of cell death and cell birth, as well as the mutation rate, using the site frequency spectrum alone. Although the discussion is framed in terms of tumor dynamics, our results apply to any exponentially growing population of individuals undergoing neutral mutations.


title: Design and Analysis of a Logless Dynamic Reconfiguration Protocol

abstract: Distributed replication systems based on the replicated state machine model have become ubiquitous as the foundation of modern database systems. To ensure availability in the presence of faults, these systems must be able to dynamically replace failed nodes with healthy ones via dynamic reconfiguration. MongoDB is a document oriented database with a distributed replication mechanism derived from the Raft protocol. In this paper, we present MongoRaftReconfig, a novel dynamic reconfiguration protocol for the MongoDB replication system. MongoRaftReconfig utilizes a logless approach to managing configuration state and decouples the processing of configuration changes from the main database operation log. The protocol's design was influenced by engineering constraints faced when attempting to redesign an unsafe, legacy reconfiguration mechanism that existed previously in MongoDB. We provide a safety proof of MongoRaftReconfig, along with a formal specification in TLA+. To our knowledge, this is the first published safety proof and formal specification of a reconfiguration protocol for a Raft-based system. We also present results from model checking its safety properties on finite protocol instances. Finally, we discuss the conceptual novelties of MongoRaftReconfig, how it can be understood as an optimized and generalized version of the single server reconfiguration algorithm of Raft, and present an experimental evaluation of how its optimizations can provide performance benefits for reconfigurations.


title: Existence of disjoint frequently hypercyclic operators which fail to be disjoint weakly mixing

abstract: In this short note, we answer a question of Martin and Sanders [Integr. Equ. Oper. Theory, 85 (2) (2016), 191-220] by showing the existence of disjoint frequently hypercyclic operators which fail to be disjoint weakly mixing and, therefore, fail to satisfy the Disjoint Hypercyclicity Criterion. We also show that given an operator $T$ such that $T \oplus T$ is frequently hypercyclic, the set of operators $S$ such that $T, S$ are disjoint frequently hypercyclic but fail to satisfy the Disjoint Hypercyclicity Criterion is SOT dense in the algebra of bounded linear operators.


title: The Talbot effect as the fundamental solution to the free Schr\"odinger equation

abstract: The Talbot effect is usually modeled using the Helmholtz equation, but its main experimental features are captured by the solution to the free Schr\"odinger equation with the Dirac comb as initial datum. This simplified description is a consequence of the paraxial approximation in geometric optics. However, it is a heuristic approximation that is not mathematically well justified, so K. I. Oskolkov raised the problem of "mathematizing" it. We show that it holds exactly in the sense of distributions.


title: An experimental demonstration of the memristor test

abstract: A simple and unambiguous test has been recently suggested [J. Phys. D: Applied Physics, 52, 01LT01 (2018)] to check experimentally if a resistor with memory is indeed a memristor, namely a resistor whose resistance depends only on the charge that flows through it, or on the history of the voltage across it. However, although such a test would represent the litmus test for claims about memristors (in the ideal sense), it has yet to be applied widely to actual physical devices. In this paper, we experimentally apply it to a current-carrying wire interacting with a magnetic core, which was recently claimed to be a memristor (so-called `$\Phi$ memristor') [J. Appl. Phys. 125, 054504 (2019)]. The results of our experiment demonstrate unambiguously that this `$\Phi$ memristor' is not a memristor: it is simply an inductor with memory. This demonstration casts further doubts that ideal memristors do actually exist in nature or may be easily created in the lab.


title: Expanded view of electron-hole recollisions in solid-state high-order harmonic generation: Full-Brillouin-zone tunneling and imperfect recollisions

abstract: We theoretically investigate electron-hole recollisions in high-harmonic generation (HHG) in band-gap solids irradiated by linearly and elliptically polarized drivers. We find that in many cases the emitted harmonics do not originate in electron-hole pairs created at the minimum band gap, where the tunneling probability is maximized, but rather in pairs created across an extended region of the Brillouin zone (BZ). In these situations, the analogy to gas-phase HHG in terms of the short- and long-trajectory categorizations is inadequate. Our analysis methodology comprises three complementary levels of theory: the numerical solutions to the semiconductor Bloch equations, an extended semiclassical recollision model, and a quantum wave packet approach. We apply this methodology to two general material types with representative band structures: a bulk system and a hexagonal monolayer system. In the bulk, the interband harmonics generated using elliptically-polarized drivers are found to originate not from tunneling at the minimum band gap $\Gamma$, but from regions away from it. In the monolayer system driven by linearly-polarized pulses, tunneling regions near different symmetry points in the BZ lead to distinct harmonic energies and emission profiles. We show that the imperfect recollisions, where an electron-hole pair recollide while being spatially separated, are important in both bulk and monolayer materials. The excellent agreement between our three levels of theory highlights and characterizes the complexity behind the HHG emission dynamics in solids, and expands on the notion of interband HHG as always originating in trajectories tunnelled at the minimum band gap. Our work furthers the fundamental understanding of HHG in periodic systems and will benefit the future design of experiments.


title: Modular Design Patterns for Hybrid Learning and Reasoning Systems: a taxonomy, patterns and use cases

abstract: The unification of statistical (data-driven) and symbolic (knowledge-driven) methods is widely recognised as one of the key challenges of modern AI. Recent years have seen large number of publications on such hybrid neuro-symbolic AI systems. That rapidly growing literature is highly diverse and mostly empirical, and is lacking a unifying view of the large variety of these hybrid systems. In this paper we analyse a large body of recent literature and we propose a set of modular design patterns for such hybrid, neuro-symbolic systems. We are able to describe the architecture of a very large number of hybrid systems by composing only a small set of elementary patterns as building blocks. The main contributions of this paper are: 1) a taxonomically organised vocabulary to describe both processes and data structures used in hybrid systems; 2) a set of 15+ design patterns for hybrid AI systems, organised in a set of elementary patterns and a set of compositional patterns; 3) an application of these design patterns in two realistic use-cases for hybrid AI systems. Our patterns reveal similarities between systems that were not recognised until now. Finally, our design patterns extend and refine Kautz' earlier attempt at categorising neuro-symbolic architectures.


title: Magic squares, the symmetric group and M\"obius randomness

abstract: Diaconis and Gamburd computed moments of secular coefficients in the CUE ensemble. We use the characteristic map to give a new combinatorial proof of their result. We also extend their computation to moments of traces of symmetric powers, where the same result holds but in a wider range. Our combinatorial proof is inspired by gcd matrices, as used by Vaughan and Wooley and by Granville and Soundararajan. We use these CUE computations to suggest a conjecture about moments of characters sums twisted by the Liouville (or by the M\"obius) function, and establish a version of it in function fields. The moral of our conjecture (and its verification in function fields) is that the Steinhaus random multiplicative function is a good model for the Liouville (or for the M\"obius) function twisted by a random Dirichlet character. We also evaluate moments of secular coefficients and traces of symmetric powers, without any condition on the size of the matrix. As an application we give a new formula for a matrix integral that was considered by Keating, Rodgers, Roditty-Gershon and Rudnick in their study of the $k$-fold divisor function.


title: Sloppy Gear Mechanism for Coupled Stochastic Transportation: from anti-equilibrium flow to infinite selectivity

abstract: Non-equilibrium transportation of particles through a restricted space (such as porous media or narrow channels) significantly differs from free space. With a simple model of two types of particles competing to transport via a passive single-lane channel connecting two chemical baths, we find two dynamical modes of transportation -- the dud mode and the ratchet mode. At the ratchet mode, the gradient flow of one type of particle forces the other into a transient anti-equilibrium flow against its gradient. At the dud mode, the two types of particles both flow according to their gradients. This counter-intuitive non-equilibrium effect can be explained by a sloppy gear mechanism. In addition to the anti-gradient flow observed in the ratchet mode, we find that the non-equilibrium sloppy gear mechanism suppresses the flow of one particle while allowing the flow for the other. This mechanism provides a general theoretical framework to explain and design the selective transportation of particles via passive channels.


title: A Scaling Limit for Utility Indifference Prices in the Discretized Bachelier Model

abstract: We consider the discretized Bachelier model where hedging is done on an equidistant set of times. Exponential utility indifference prices are studied for path-dependent European options and we compute their non-trivial scaling limit for a large number of trading times $n$ and when risk aversion is scaled like $n\ell$ for some constant $\ell>0$. Our analysis is purely probabilistic. We first use a duality argument to transform the problem into an optimal drift control problem with a penalty term. We further use martingale techniques and strong invariance principles and get that the limiting problem takes the form of a volatility control problem.


title: Experimental measurement of the isolated magnetic susceptibility

abstract: The isolated susceptibility $\chi_{\rm I}$ may be defined as a (non-thermodynamic) average over the canonical ensemble, but while it has often been discussed in the literature, it has not been clearly measured. Here, we demonstrate an unambiguous measurement of $\chi_{\rm I}$ at avoided nuclear-electronic level crossings in a dilute spin ice system, containing well-separated holmium ions. We show that $\chi_{\rm I}$ quantifies the superposition of quasi-classical spin states at these points, and is a direct measure of state concurrence and populations.


title: Rotor-routing reachability is easy, chip-firing reachability is hard

abstract: Chip-firing and rotor-routing are two well-studied examples of abelian networks. We study the complexity of their respective reachability problems. We show that the rotor-routing reachability problem is decidable in polynomial time, and we give a simple characterization of when a chip-and-rotor configuration is reachable from another one. For chip-firing, it has been known that the reachability problem is in P if we have a class of graphs whose period length is polynomial (for example, Eulerian digraphs). Here we show that in the general case, chip-firing reachability is hard in the sense that if the chip-firing reachability problem were in P for general digraphs, then the polynomial hierarchy would collapse to NP. We encode graphs by their adjacency matrix, and we encode ribbon structures "succinctly", only remembering the number of consecutive parallel edges.


title: Brain Waves Analysis Via a Non-parametric Bayesian Mixture of Autoregressive Kernels

abstract: The standard approach to analyzing brain electrical activity is to examine the spectral density function (SDF) and identify predefined frequency bands that have the most substantial relative contributions to the overall variance of the signal. However, a limitation of this approach is that the precise frequency and bandwidth of oscillations vary with cognitive demands. Thus they should not be arbitrarily defined a priori in an experiment. In this paper, we develop a data-driven approach that identifies (i) the number of prominent peaks, (ii) the frequency peak locations, and (iii) their corresponding bandwidths (or spread of power around the peaks). We propose a Bayesian mixture auto-regressive decomposition method (BMARD), which represents the standardized SDFas a Dirichlet process mixture based on a kernel derived from second-order auto-regressive processes which completely characterize the location (peak)and scale (bandwidth) parameters. We present a Metropolis-Hastings within Gibbs algorithm to sample from the posterior distribution of the mixture parameters. Simulation studies demonstrate the robustness and performance of the BMARD method. Finally, we use the proposed BMARD method to analyze local field potential (LFP) activity from the hippocampus of laboratory rats across different conditions in a non-spatial sequence memory experiment to identify the most interesting frequency bands and examine the link between specific patterns of activity and trial-specific cognitive demands.


title: Do Transformer Modifications Transfer Across Implementations and Applications?

abstract: The research community has proposed copious modifications to the Transformer architecture since it was introduced over three years ago, relatively few of which have seen widespread adoption. In this paper, we comprehensively evaluate many of these modifications in a shared experimental setting that covers most of the common uses of the Transformer in natural language processing. Surprisingly, we find that most modifications do not meaningfully improve performance. Furthermore, most of the Transformer variants we found beneficial were either developed in the same codebase that we used or are relatively minor changes. We conjecture that performance improvements may strongly depend on implementation details and correspondingly make some recommendations for improving the generality of experimental results.


title: On the Fokker-Planck approximation in the kinetic equation of multicomponent classical nucleation theory

abstract: We examine the validity of the Fokker-Planck equation with linear force coefficients as an approximation to the kinetic equation of nucleation in homogeneous isothermal multicomponent condensation. Starting with a discrete equation of balance governing the temporal evolution of the distribution function of an ensemble of multicomponent droplets and reducing it (by means of Taylor series expansions) to the differential form in the vicinity of the saddle point of the free energy surface, we have identified the parameters whereof the smallness is necessary for the resulting kinetic equation to have the form of the Fokker-Planck equation with linear (in droplet variables) force coefficients. The "non-smallness" of these parameters results either in the appearance of the third or higher order partial derivatives of the distribution function in the kinetic equation or in its force coefficients becoming non-linear functions of droplet variables, or both; this would render the conventional kinetic equation of multicomponent nucleation and its predictions inaccurate. As a numerical illustration, we carried out calculations for isothermal condensation in five binary systems of various non-ideality at T=293.15 K: 1-butanol--1-hexanol, water--methanol, water--ethanol, water--1-propanol, water--1-butanol. Our results suggest that under typical experimental conditions the kinetic equation of binary nucleation of classical nucleation theory may require a two-fold modification and, hence, the conventional expression for the steady-state binary nucleation rate may not be adequate for the consistent comparison of theoretical predictions with experimental data.


title: Hospital management in the COVID-19 emergency: Abelian Sandpile paradigm and beyond

abstract: In this article, we propose a mathematical model -- based on a cellular automaton -- for the redistribution of patients within a network of hospitals with limited available resources, in order to reduce the risks of a local/global collapse of the healthcare system. We attempt at developing a conceptual tool to support making rational decisions relevant to the optimisation of the allocation of patients into accessible medical facilities. The strategy is based on a version of the Abelian Sandpile model for the Self-Organised Criticality, with the idea of testing the paradigm for the management of patients among the COVID-19 hospitals in Italian regions. In particular, we compare the novel proposal to the standard management of connections between hospitals, showing a number of advantages at a local and global level, by means of a reliable indicator function introduced for measuring the effectiveness of the allocation strategies.


title: Fermi acceleration in relativistic collisionless plasma shocks correlates with anisotropic energy gains

abstract: Collisionless shocks generated by two colliding relativistic electron-positron plasma shells are studied using particle-in-cell (PIC) simulations. Shocks are mediated by the Weibel instability (WI), and the kinetic energy of the fastest accelerated particles is found to be anisotropically modified by WI-induced electric fields. Specifically, we show that all particles interacting with the shock bifurcate into two groups based on their final relativistic Lorentz factor $\gamma$: slow ($\gamma < \gamma_{bf}$) and fast ($\gamma > \gamma_{bf}$), where $\gamma_{bf}$ is the bifurcation Lorentz factor that was found to be approximately twice the initial (upstream) Lorentz factor $\gamma_0$. We have found that the energies of the slow particles are equally affected by the longitudinal and transverse components of the shock electric field, whereas the fast particles are primarily accelerated by the transverse field component.


title: Learner-Private Convex Optimization

abstract: Convex optimization with feedback is a framework where a learner relies on iterative queries and feedback to arrive at the minimizer of a convex function. It has gained considerable popularity thanks to its scalability in large-scale optimization and machine learning. The repeated interactions, however, expose the learner to privacy risks from eavesdropping adversaries that observe the submitted queries. In this paper, we study how to optimally obfuscate the learner's queries in convex optimization with first-order feedback, so that their learned optimal value is provably difficult to estimate for an eavesdropping adversary. We consider two formulations of learner privacy: a Bayesian formulation in which the convex function is drawn randomly, and a minimax formulation in which the function is fixed and the adversary's probability of error is measured with respect to a minimax criterion. Suppose that the learner wishes to ensure the adversary cannot estimate accurately with probability greater than $1/L$ for some $L>0$. Our main results show that the query complexity overhead is additive in $L$ in the minimax formulation, but multiplicative in $L$ in the Bayesian formulation. Compared to existing learner-private sequential learning models with binary feedback, our results apply to the significantly richer family of general convex functions with full-gradient feedback. Our proofs learn on tools from the theory of Dirichlet processes, as well as a novel strategy designed for measuring information leakage under a full-gradient oracle.

categories: stat.ML cs.CR cs.LG math.OC math.ST stat.TH


title: The Sherrington-Kirkpatrick model for spin glasses: A new approach for the solution

abstract: We discuss the Sherrington-Kirkpatrick mean-field version of a spin glass within the distributional zeta-function method (DZFM). In the DZFM, since the dominant contribution to the average free energy is written as a series of moments of the partition function of the model, the spin-glass multivalley structure is obtained. Also, an exact expression for the saddle points corresponding to each valley and a global critical temperature showing the existence of many stables or at least metastables equilibrium states is presented. Near the critical point we obtain analytical expressions of the order parameters that are in agreement with phenomenological results. We evaluate the linear and nonlinear susceptibility and we find the expected singular behavior at the spin-glass critical temperature. Furthermore, we obtain a positive definite expression for the entropy and we show that ground-state entropy tends to zero as the temperature goes to zero. We show that our solution is stable for each term in the expansion. Finally, we analyze the behavior of the overlap distribution, where we find a general expression for each moment of the partition function.


title: Mean curvature flow with generic low-entropy initial data

abstract: We prove that sufficiently low-entropy closed hypersurfaces can be perturbed so that their mean curvature flow encounters only spherical and cylindrical singularities. Our theorem applies to all closed surfaces in $\mathbb{R}^3$ with entropy $\leq 2$ and to all closed hypersurfaces in $\mathbb{R}^4$ with entropy $\leq \lambda(\mathbb{S}^1 \times \mathbb{R}^2)$. When combined with recent work of Daniels-Holgate, this strengthens Bernstein-Wang's low-entropy Schoenflies-type theorem by relaxing the entropy bound to $\lambda(\mathbb{S}^1 \times \mathbb{R}^2)$. Our techniques, based on a novel density drop argument, also lead to a new proof of generic regularity result for area-minimizing hypersurfaces in eight dimensions (due to Hardt-Simon and Smale).


title: Contact statistics in populations of noninteracting random walkers in two dimensions

abstract: The interaction between individuals in biological populations, dilute components of chemical systems, or particles transported by turbulent flows depends critically on their contact statistics. This work clarifies those statistics under the simplifying assumptions that the underlying motions approximate a Brownian random walk and that the particles can be treated as noninteracting. We measure the contact-interval (also called the waiting-time or inter-arrival-time), contact-count, and contact-duration distributions in populations of individuals undergoing noninteracting continuous-space-time random walks on a periodic two-dimensional plane (a torus), as functions of the population number density, walker radius, and random-walk step size. The contact-interval is exponentially distributed for times longer than the mean-free-collision time but not for times shorter than that, and the contact duration distribution is strongly peak at the ballistic-crossing time for head-on collisions, when the ballistic-crossing time is short compared to the mean step duration. While successive contacts between individuals are independent, the probability of repeat contact decreases with time after a previous contact. This leads to a negative duration dependence of the waiting-time interval and over dispersion of the contact-count probability density function for all time intervals. The paper demonstrates that for populations of small particles (walker radius small compared to the mean-separation or random-walk step size) the mean-free-collision interval, the ballistic-crossing time, and the random-walk-step duration can be used to construct temporal scalings which allow for common waiting-time, contact-count, and contact-duration distributions across different populations. Semi-analytic approximations for both the waiting-time and contact-duration distributions are also presented.

categories: math-ph cond-mat.stat-mech math.MP q-bio.PE


title: Decomposition Methods for Global Solutions of Mixed-Integer Linear Programs

abstract: This paper introduces two decomposition-based methods for two-block mixed-integer linear programs (MILPs), which aim to take advantage of separable structures of the original problem by solving a sequence of lower-dimensional MILPs. The first method is based on the $\ell_1$-augmented Lagrangian method (ALM), and the second one is based on a modified alternating direction method of multipliers (ADMM). In the presence of certain block-angular structures, both methods create parallel subproblems in one block of variables, and add nonconvex cuts to update the other block; they converge to globally optimal solutions of the original MILP under proper conditions. Numerical experiments on three classes of MILPs demonstrate the advantages of the proposed methods on structured problems over the state-of-the-art MILP solvers.


title: Resolving the LMA-dark NSI degeneracy with coherent neutrino-nucleus scattering

abstract: In the presence of non-standard neutrino interactions (NSI), a degeneracy exists in neutrino oscillation data, which involves the flipping of the octant of the mixing angle ${\theta_{12}}$ and the type of the neutrino mass ordering. In this article, we revisit the status of this degeneracy in the light of recent data on coherent elastic neutrino-nucleus scattering (CE${\nu}$NS) from the COHERENT experiment. For general relative couplings to up and down quarks, the degeneracy is disfavoured at the $2{\sigma}$ level by the latest data but remains at a higher confidence level. We investigate the requirements of future CE${\nu}$NS measurements to resolve the degeneracy with high significance. We find that a measurement involving both, electron and muon neutrino flavours and a target with a neutron-to-proton ratio close to 1 is required. For example, an experiment with a silicon target at the European Spallation Source can resolve the degeneracy at more than $4{\sigma}$ for arbitrary relative couplings to up and down quarks.


title: Observation of vacuum-induced collective quantum beats

abstract: We demonstrate collectively enhanced vacuum-induced quantum beat dynamics from a three-level V-type atomic system. Exciting a dilute atomic gas of magneto-optically trapped $^{85}$Rb atoms with a weak drive resonant on one of the transitions, we observe the forward-scattered field after a sudden shut-off of the laser. The subsequent radiative dynamics, measured for various optical depths of the atomic cloud, exhibits superradiant decay rates, as well as collectively enhanced quantum beats. Our work is also the first experimental illustration of quantum beats arising from atoms initially prepared in a single excited level as a result of the vacuum-induced coupling between excited levels.


title: Lotka's Law and Pattern of Author Productivity in the Field of Brain Concussion Research: A Scientometric Analysis

abstract: The present study deals a scientometric analysis of 8486 bibliometric publications retrieved from the Web of Science database during the period 2008 to 2017. Data is collected and analyzed using Bibexcel software. The study focuses on various aspect of the quantitative research such as growth of papers (year wise), Collaborative Index (CI), Degree of Collaboration (DC), Co-authorship Index (CAI), Collaborative Co-efficient (CC), Modified Collaborative Co-Efficient (MCC), Lotka's Exponent value, Kolmogorov-Smirnov test (K-S Test).


title: Strain-tuning of nematicity and superconductivity in single crystals of FeSe

abstract: Strain is a powerful experimental tool to explore new electronic states and understand unconventional superconductivity. Here, we investigate the effect of uniaxial strain on the nematic and superconducting phase of single crystal FeSe using magnetotransport measurements. We find that the resistivity response to the strain is strongly temperature dependent and it correlates with the sign change in the Hall coefficient being driven by scattering, coupling with the lattice and multiband phenomena. Band structure calculations suggest that under strain the electron pockets develop a large in-plane anisotropy as compared with the hole pocket. Magnetotransport studies at low temperatures indicate that the mobility of the dominant carriers increases with tensile strain. Close to the critical temperature, all resistivity curves at constant strain cross in a single point, indicating a universal critical exponent linked to a strain-induced phase transition. Our results indicate that the superconducting state is enhanced under compressive strain and suppressed under tensile strain, in agreement with the trends observed in FeSe thin films and overdoped pnictides, whereas the nematic phase seems to be affected in the opposite way by the uniaxial strain. By comparing the enhanced superconductivity under strain of different systems, our results suggest that strain on its own cannot account for the enhanced high $T_c$ superconductivity of FeSe systems.


title: Strong-coupling expansion of multi-band interacting models: mapping onto the transverse-field $J_1$-$J_2$ Ising model

abstract: We investigate a class of two-dimensional two-band microscopic models in which the inter-band repulsive interactions play the dominant role. We first demonstrate three different schemes of constraining the ratios between the three types of inter-band interactions -- density-density, spin exchange, and pair-hopping -- that render the model free of the fermionic sign-problem for any filling and, consequently, amenable to efficient Quantum Monte Carlo simulations. We then study the behavior of these sign-problem-free models in the strong-coupling regime. In the cases where spin-rotational invariance is preserved or lowered to a planar symmetry, the strong-coupling ground state is a quantum paramagnet. However, in the case where there is only a residual Ising symmetry, the strong-coupling expansion maps onto the transverse-field $J_1$-$J_2$ Ising model, whose pseudospins are associated with local inter-band magnetic order. We show that by varying the band structure parameters within a reasonable range of values, a variety of ground states and quantum critical points can be accessed in the strong-coupling regime, some of which are not realized in the weak-coupling regime. We compare these results with the case of the single-band Hubbard model, where only intra-band repulsion is present, and whose strong-coupling behavior is captured by a simple Heisenberg model.


title: Cone theorem and Mori hyperbolicity

abstract: We discuss the cone theorem for quasi-log schemes and the Mori hyperbolicity. In particular, we establish that the log canonical divisor of a Mori hyperbolic projective normal pair is nef if it is nef when restricted to the non-lc locus. This answers Svaldi's question completely. We also treat the uniruledness of the degenerate locus of an extremal contraction morphism for quasi-log schemes. Furthermore, we prove that every fiber of a relative quasi-log Fano scheme is rationally chain connected modulo the non-qlc locus.


title: Nonconvex integro-differential sweeping process with applications

abstract: In this paper, we analyze and discuss the well-posedness of a new variant of the so-called sweeping process, introduced by J.J. Moreau in the early 70's \cite{More71} with motivation in plasticity theory. In this variant, the normal cone to the (mildly non-convex) prox-regular moving set $C(t)$, supposed to have an absolutely continuous variation, is perturbed by a sum of a Carath\'{e}odory mapping and an integral forcing term. The integrand of the forcing term depends on two time-variables, that is, we study a general integro-differential sweeping process of Volterra type. By setting up an appropriate semi-discretization method combined with a new Gronwall-like inequality (differential inequality), we show that the integro-differential sweeping process has one and only one absolutely continuous solution. We also establish the continuity of the solution with respect to the initial value. The results of the paper are applied to the study of nonlinear integro-differential complementarity systems which are combination of Volterra integro-differential equations with nonlinear complementarity constraints. Another application is concerned with non-regular electrical circuits containing time-varying capacitors and nonsmooth electronic device like diodes. Both applications represent an additional novelty of our paper.


title: Dispersion and the Speed-Limited Particle-in-Cell Algorithm

abstract: This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. These results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.


title: Maximality of Seidel matrices and switching roots of graphs

abstract: In this paper, we discuss maximality of Seidel matrices with a fixed largest eigenvalue. We present a classification of maximal Seidel matrices of largest eigenvalue $3$, which gives a classification of maximal equiangular lines in a Euclidean space with angle $\arccos1/3$. Motivated by the maximality of the exceptional root system $E_8$, we define strong maximality of a Seidel matrix, and show that every Seidel matrix achieving the absolute bound is strongly maximal.


title: Surface Gravity of Rotating Dumbbell Shapes

abstract: We investigate the problem of determining the shape of a rotating celestial object - e.g., a comet or an asteroid - under its own gravitational field. More specifically, we consider an object symmetric with respect to one axis - such as a dumbbell - that rotates around a second axis perpendicular to the symmetry axis. We assume that the object can be modeled as an incompressible fluid of constant mass density, which is regarded as a first approximation of an aggregate of particles. In the literature, the gravitational field of a body is often described as a multipolar expansion involving spherical coordinates (Kaula, 1966). In this work we describe the shape in terms of cylindrical coordinates, which are most naturally adapted to the symmetry of the body, and we express the gravitational potential generated by the rotating body as a simple formula in terms of elliptic integrals. An equilibrium shape occurs when the gravitational potential energy and the rotational kinetic energy at the surface of the body balance each other out. Such an equilibrium shape can be derived as a solution of an optimization problem, which can be found via the variational method. We give an example where we apply this method to a two-parameter family of dumbbell shapes, and find approximate numerical solutions to the corresponding optimization problem.

categories: astro-ph.EP math-ph math.DS math.MP physics.space-ph