title: Risk-Aware Reasoning for Autonomous Vehicles

abstract: A significant barrier to deploying autonomous vehicles (AVs) on a massive scale is safety assurance. Several technical challenges arise due to the uncertain environment in which AVs operate such as road and weather conditions, errors in perception and sensory data, and also model inaccuracy. In this paper, we propose a system architecture for risk-aware AVs capable of reasoning about uncertainty and deliberately bounding the risk of collision below a given threshold. We discuss key challenges in the area, highlight recent research developments, and propose future research directions in three subsystems. First, a perception subsystem that detects objects within a scene while quantifying the uncertainty that arises from different sensing and communication modalities. Second, an intention recognition subsystem that predicts the driving-style and the intention of agent vehicles (and pedestrians). Third, a planning subsystem that takes into account the uncertainty, from perception and intention recognition subsystems, and propagates all the way to control policies that explicitly bound the risk of collision. We believe that such a white-box approach is crucial for future adoption of AVs on a large scale.


title: The $O(N)$ Model in $4<d<6$: Instantons and Complex CFTs

abstract: We revisit the scalar $O(N)$ model in the dimension range $4<d<6$ and study the effects caused by its metastability. As shown in previous work, this model formally possesses a fixed point where, perturbatively in the $1/N$ expansion, the operator scaling dimensions are real and above the unitarity bound. Here, we further show that these scaling dimensions do acquire small imaginary parts due to the instanton effects. In $d$ dimensions and for large $N$, we find that they are of order $e^{-N f(d)}$, where, remarkably, the function $f(d)$ equals the sphere free energy of a conformal scalar in $d-2$ dimensions. The non-perturbatively small imaginary parts also appear in other observables, such as the sphere free energy and two and three-point function coefficients, and we present some of their calculations. Therefore, at sufficiently large $N$, the $O(N)$ models in $4<d<6$ may be thought of as complex CFTs. When $N$ is large enough for the imaginary parts to be numerically negligible, the five-dimensional $O(N)$ models may be studied using the techniques of numerical bootstrap.


title: Weyl group characters afforded by zero weight spaces

abstract: Let $G$ be a simple complex Lie group with Weyl group $W$. We give a formula for the character of $W$ on the zero weight space of any finite dimensional representation of $G$. The formula involves partition functions, generalizing Kostant's partition function. On the elliptic set of $W$ the partition functions are trivial. On the elliptic regular set, the character formula is a monomial product of certain co-roots, up to a constant equal to $0$ or $\pm 1$. For a Coxeter element we recover Kostant's formula for this trace. If the long element $w_0=-1$, our formula leads to a method for determining all representations of $G$ for which the zero weight space is irreducible.


title: Resource Preservability

abstract: Resource theory is a general, model-independent approach aiming to understand the qualitative notion of resource quantitatively. In a given resource theory, free operations are physical processes that do not create the resource and are considered zero-cost. This brings the following natural question: For a given free operation, what is its ability to preserve a resource? We axiomatically formulate this ability as the resource preservability, which is constructed as a channel resource theory induced by a state resource theory. We provide two general classes of resource preservability monotones: One is based on state resource monotones, and another is based on channel distance measures. Specifically, the latter gives the robustness monotone, which has been recently found to have an operational interpretation. As examples, we show that athermality preservability of a Gibbs-preserving channel can be related to the smallest bath size needed to thermalize all its outputs, and it also bounds the capacity of a classical communication scenario under certain thermodynamic constraints. We further apply our theory to the study of entanglement preserving local thermalization (EPLT) and provide a new family of EPLT which admits arbitrarily small nonzero entanglement preservability and free entanglement preservation at the same time. Our results give the first systematic and general formulation of the resource preservation character of free operations.


title: On the Probabilistic Degrees of Symmetric Boolean functions

abstract: The probabilistic degree of a Boolean function $f:\{0,1\}^n\rightarrow \{0,1\}$ is defined to be the smallest $d$ such that there is a random polynomial $\mathbf{P}$ of degree at most $d$ that agrees with $f$ at each point with high probability. Introduced by Razborov (1987), upper and lower bounds on probabilistic degrees of Boolean functions --- specifically symmetric Boolean functions --- have been used to prove explicit lower bounds, design pseudorandom generators, and devise algorithms for combinatorial problems. In this paper, we characterize the probabilistic degrees of all symmetric Boolean functions up to polylogarithmic factors over all fields of fixed characteristic (positive or zero).


title: Resolving asset pricing puzzles using price-impact

abstract: We solve in closed-form an equilibrium model in which a finite number of exponential investors continuously consume and trade with price-impact. Compared to the analogous Pareto-efficient equilibrium model, price-impact has an amplification effect on risk-sharing distortions that helps resolve the interest rate puzzle and the stock-price volatility puzzle and, to a lesser extent, affects the equity premium puzzle.


title: A Physical Framework for the Earth System in the Anthropocene: Towards an Accountancy System

abstract: At a time when humanity has achieved global dominance at a scale that was previously thought impossible, it might also face an existential threat due to the consequences of that overwhelming influence on our common home, the Earth System (ES). In this work we explore how Physics may help us to understand the transitions that the ES is going through and lead us to a physically motivated accounting system that allows for setting boundaries to our negative influence on the ecosystems.


title: Electromagnetic Energy Balance Equations and Poynting Theorem

abstract: Poynting theorem plays a very important role in analyzing electromagnetic phenomena. The electromagnetic power flux density is usually expressed with the Poynting vector. However, since Poynting theorem basically focuses on the power balance in a system, it is not so efficient in some situations to use it for evaluating the electromagnetic energies. The energy balance issue for time varying fields is revisited in this paper. Energy balance equations are introduced and a modified method for evaluating power flux is proposed.


title: On the Existence of Block-Diagonal Solutions to Lyapunov and $\mathcal{H}_{\infty}$ Riccati Inequalities

abstract: In this paper, we describe sufficient conditions when block-diagonal solutions to Lyapunov and $\mathcal{H}_{\infty}$ Riccati inequalities exist. In order to derive our results, we define a new type of comparison systems, which are positive and are computed using the state-space matrices of the original (possibly nonpositive) systems. Computing the comparison system involves only the calculation of $\mathcal{H}_{\infty}$ norms of its subsystems. We show that the stability of this comparison system implies the existence of block-diagonal solutions to Lyapunov and Riccati inequalities. Furthermore, our proof is constructive and the overall framework allows the computation of block-diagonal solutions to these matrix inequalities with linear algebra and linear programming. Numerical examples illustrate our theoretical results.


title: Approximation algorithms for maximally balanced connected graph partition

abstract: Given a simple connected graph $G = (V, E)$, we seek to partition the vertex set $V$ into $k$ non-empty parts such that the subgraph induced by each part is connected, and the partition is maximally balanced in the way that the maximum cardinality of these $k$ parts is minimized. We refer this problem to as {\em min-max balanced connected graph partition} into $k$ parts and denote it as {\sc $k$-BGP}. The general vertex-weighted version of this problem on trees has been studied since about four decades ago, which admits a linear time exact algorithm; the vertex-weighted {\sc $2$-BGP} and {\sc $3$-BGP} admit a $5/4$-approximation and a $3/2$-approximation, respectively; but no approximability result exists for {\sc $k$-BGP} when $k \ge 4$, except a trivial $k$-approximation. In this paper, we present another $3/2$-approximation for our cardinality {\sc $3$-BGP} and then extend it to become a $k/2$-approximation for {\sc $k$-BGP}, for any constant $k \ge 3$. Furthermore, for {\sc $4$-BGP}, we propose an improved $24/13$-approximation. To these purposes, we have designed several local improvement operations, which could be useful for related graph partition problems.


title: Spiked Laplacian Graphs: Bayesian Community Detection in Heterogeneous Networks

abstract: In network data analysis, it is becoming common to work with a collection of graphs that exhibit \emph{heterogeneity}. For example, neuroimaging data from patient cohorts are increasingly available. A critical analytical task is to identify communities, and graph Laplacian-based methods are routinely used. However, these methods are currently limited to a single network and do not provide measures of uncertainty on the community assignment. In this work, we propose a probabilistic network model called the ``Spiked Laplacian Graph'' that considers each network as an invertible transform of the Laplacian, with its eigenvalues modeled by a modified spiked structure. This effectively reduces the number of parameters in the eigenvectors, and their sign patterns allow efficient estimation of the community structure. Further, the posterior distribution of the eigenvectors provides uncertainty quantification for the community estimates. Subsequently, we introduce a Bayesian non-parametric approach to address the issue of heterogeneity in a collection of graphs. Theoretical results are established on the posterior consistency of the procedure and provide insights on the trade-off between model resolution and accuracy. We illustrate the performance of the methodology on synthetic data sets, as well as a neuroscience study related to brain activity in working memory. Keywords: Hierarchical Community Detection, Isoperimetric Constant, Mixed-Effect Eigendecomposition, Normalized Graph Cut, Stiefel Manifold


title: C*-algebras of self-similar actions of groupoids on higher-rank graphs and their equilibrium states

abstract: We introduce the notion of a self-similar action of a groupoid G on a finite higher-rank graph. To these actions we associate a compactly aligned product system of Hilbert bimodules, and we show that the corresponding Nica-Toeplitz and Cuntz-Pimsner algebras are universal for generators and relations. We consider natural actions of the real numbers on both algebras and study the KMS states of the associated dynamics. For large inverse temperatures, we describe the simplex of KMS states on the Nica-Toeplitz algebra in terms of traces on the full C*-algebra of G. To study the KMS structure of the Cuntz-Pimsner algebra, we restrict to strongly-connected finite higher-rank graphs and a preferred dynamics. We generalise the G-periodicity group of Li and Yang, which is built to encode the periodicity of the underlying graph in the presence of the action of G. We prove that the KMS states of the Cuntz-Pimsner algebra are parametrised by states on the C*-algebra of the G-periodicity group, and we show that if the graph is G-aperiodic and the action satisfies a finite-state condition, then there is at most one KMS state on the Cuntz-Pimsner algebra. We illustrate our results by introducing the notion of a coloured-graph automaton, which we use to construct examples of self-similar actions. We compute the unique KMS states of the Cuntz-Pimsner algebra for some concrete examples.


title: Generalized period-index problem with an application to quadratic forms

abstract: Let $F$ be the function field of a curve over a complete discretely valued field. Let $\ell$ be a prime not equal to the characteristic of the residue field. Given a finite subgroup $B$ in the $\ell$ torsion part of the Brauer group ${}_{\ell}Br(F)$, we define the index of $B$ as the minimum of the degrees of field extensions which split all elements in $B$. In this manuscript, we give an upper bound for the index of any finite subgroup $B$ in terms of arithmetic invariants of $F$. As a simple application of our result, given a quadratic form $q/F$, where $F$ is the function field of a curve over an $n$-local field, we provide an upper bound to the minimum of degrees of field extensions $L/F$ so that the Witt index of $q\otimes L$ becomes the largest possible.


title: Neither Presentism nor Eternalism

abstract: Is reality three-dimensional and becoming real (Presentism), or is reality four-dimensional and becoming illusory (Eternalism)? Both options raise difficulties. I argue that we do not need to be trapped by this dilemma. There is a third possibility: reality has a more complex temporal structure than either of these two naive options. Fundamental becoming is real, but local and unoriented. A notion of present is well defined, but only locally and in the context of approximations.


title: Dust-acoustic rogue waves in an electron-positron-ion-dust plasma medium

abstract: A precise theoretical investigation has been made on dust-acoustic (DA) waves (DAWs) in a four components dusty plasma medium having inertial warm adiabatic dust grains and inertialess $q$-distributed electrons as well as isothermal ions and positrons. The nonlinear and dispersive parameters of the nonlinear Schr\"{o}dinger equation (NLSE), which develops by using reductive perturbations technique, have been used to recognize the stable and unstable parametric regions of the DAWs as well as associated DA rogue waves (DARWs) in the unstable parametric regime of the DAWs. The effects of the light positrons and massive dust grains in determining the amplitude and width of the DARWs associated with DAWs are examined. The findings of our present investigation may be useful for understanding different nonlinear electrostatic phenomena in both space and laboratory plasmas.


title: Selection Games on Continuous Functions

abstract: In this paper we study the selection principle of closed discrete selection, first researched by Tkachuk in [13] and strengthened by Clontz, Holshouser in [3], in set-open topologies on the space of continuous real-valued functions. Adapting the techniques involving point-picking games on \(X\) and \(C_p(X)\), the current authors showed similar equivalences in [1] involving the compact subsets of \(X\) and \(C_k(X)\). By pursuing a bitopological setting, we have touched upon a unifying framework which involves three basic techniques: general game duality via reflections (Clontz), general game equivalence via topological connections, and strengthening of strategies (Pawlikowski and Tkachuk). Moreover, we develop a framework which identifies topological notions to match with generalized versions of the point-open game.


title: Fourier-Mukai transforms and stable sheaves on Weierstrass elliptic surfaces

abstract: On a Weierstra{\ss} elliptic surface $X$, we define a `limit' of Bridgeland stability conditions, denoted as $Z^l$-stability, by moving the polarisation towards the fiber direction in the ample cone while keeping the volume of the polarisation fixed. We describe conditions under which a slope stable torsion-free sheaf is taken by a Fourier-Mukai transform to a $Z^l$-stable object, and describe a modification upon which a $Z^l$-semistable object is taken by the inverse Fourier-Mukai transform to a slope semistable torsion-free sheaf. We also study wall-crossing for Bridgeland stability, and show that 1-dimensional twisted Gieseker semistable sheaves are taken by a Fourier-Mukai transform to Bridgeland semistable objects.


title: Exact and/or Fast Nearest Neighbors

abstract: Prior methods for retrieval of nearest neighbors in high dimensions are fast and approximate--providing probabilistic guarantees of returning the correct answer--or slow and exact performing an exhaustive search. We present Certified Cosine, a novel approach to nearest-neighbors which takes advantage of structure present in the cosine similarity distance metric to offer certificates. When a certificate is constructed, it guarantees that the nearest neighbor set is correct, possibly avoiding an exhaustive search. Certified Cosine's certificates work with high dimensional data and outperform previous exact nearest neighbor methods on these datasets.


title: Convergence of the likelihood ratio method for linear response of non-equilibrium stationary states

abstract: We consider numerical schemes for computing the linear response of steady-state averages of stochastic dynamics with respect to a perturbation of the drift part of the stochastic differential equation. The schemes are based on Girsanov's change-of-measure theory to reweight trajectories with factors derived from a linearization of the Girsanov weights. We investigate both the discretization error and the finite time approximation error. The designed numerical schemes are shown to be of bounded variance with respect to the integration time, which is a desirable feature for long time simulation. We also show how the discretization error can be improved to second order accuracy in the time step by modifying the weight process in an appropriate way.


title: Deep Radiance Caching: Convolutional Autoencoders Deeper in Ray Tracing

abstract: Rendering realistic images with global illumination is a computationally demanding task and often requires dedicated hardware for feasible runtime. Recent research uses Deep Neural Networks to predict indirect lighting on image level, but such methods are commonly limited to diffuse materials and require training on each scene.We present Deep Radiance Caching (DRC), an efficient variant of Radiance Caching utilizing Convolutional Autoencoders for rendering global illumination. DRC employs a denoising neural network with Radiance Caching to support a wide range of material types, without the requirement of offline pre-computation or training for each scene.This offers high performance CPU rendering for maximum accessibility. Our method has been evaluated on interior scenes, and is able to produce high-quality images within 180 seconds on a single CPU.


title: Learn to Explain Efficiently via Neural Logic Inductive Learning

abstract: The capability of making interpretable and self-explanatory decisions is essential for developing responsible machine learning systems. In this work, we study the learning to explain problem in the scope of inductive logic programming (ILP). We propose Neural Logic Inductive Learning (NLIL), an efficient differentiable ILP framework that learns first-order logic rules that can explain the patterns in the data. In experiments, compared with the state-of-the-art methods, we find NLIL can search for rules that are x10 times longer while remaining x3 times faster. We also show that NLIL can scale to large image datasets, i.e. Visual Genome, with 1M entities.


title: Two-band model for magnetism and superconductivity in nickelates

abstract: The recently discovered superconductivity in Nd$_{1-x}$Sr$_x$NiO$_2$ provides a new opportunity for studying strongly correlated unconventional superconductivity. The single-hole Ni$^+$ ($3d^9$) configuration in the parent compound NdNiO$_2$ is similar to that of Cu$^{2+}$ in cuprates. We suggest that after doping, the intra-orbital spin-singlet and inter-orbital spin-triplet double-hole (doublon) configurations of Ni$^{2+}$ are competing, and we construct a two-band Hubbard model by including both the $3d_{x^2-y^2}$ and $3d_{xy}$-orbitals. The effective spin-orbital super-exchange model in the undoped case is a variant of the $SU(4)$ Kugel-Khomskii model augmented by symmetry breaking terms. Upon doping, the effective exchange interactions between spin-$\frac{1}{2}$ single-holes, spin-1 (triplet) doublons, and singlet doublons are derived. Possible superconducting pairing symmetries are classified in accordance to the $D_{4h}$ crystalline symmetry, and their connections to the superexchange interactions are analyzed.


title: Auto-Rotating Perceptrons

abstract: This paper proposes an improved design of the perceptron unit to mitigate the vanishing gradient problem. This nuisance appears when training deep multilayer perceptron networks with bounded activation functions. The new neuron design, named auto-rotating perceptron (ARP), has a mechanism to ensure that the node always operates in the dynamic region of the activation function, by avoiding saturation of the perceptron. The proposed method does not change the inference structure learned at each neuron. We test the effect of using ARP units in some network architectures which use the sigmoid activation function. The results support our hypothesis that neural networks with ARP units can achieve better learning performance than equivalent models with classic perceptrons.


title: Achieving higher taxi outflows from a congested drop-off lane: a simulation-based policy study

abstract: We examine special lanes used by taxis and other shared-ride services to drop-off patrons at airport and rail terminals. Vehicles are prohibited from overtaking each other within the lane. They must therefore wait in a first-in-first-out queue during busy periods. Patrons are often discharged from vehicles only upon reaching a desired drop-off area near the terminal entrance. When wait times grow long, however, some vehicles discharge their patrons in advance of that desired area. A train station in Eastern China is selected as a case study. Its FIFO drop-off lane is presently managed by policemen who allow taxis to enter the lane in batched fashion. Inefficiencies are observed because curb space near the terminal often goes unused. This is true even when supplemental batches of taxis are released into the lane in efforts to fill those spaces. A microscopic simulation model of a FIFO drop-off lane is developed in-house, and is painstakingly calibrated to data measured at the study site. Simulation experiments indicate that rescinding the FIFO lane's present batching strategy can increase taxi outflow by more than 26 percent. Further experiments show that even greater gains can be achieved by batching taxis, but requiring them to discharge patrons when forced by downstream queues to stop a prescribed distance in advance of a desired drop-off area. Further gains were predicted by requiring the lead taxi in each batch to discharge its patron only after travelling a prescribed distance beyond a desired location. Practical implications are discussed in light of the present boom in shared-ride services.


title: Sums of Cubes in Quaternion Rings

abstract: We investigate a version of Waring's Problem over quaternion rings, focusing on cubes in quaternion rings with integer coefficients. We determine the global upper and lower bounds for the number of cubes necessary to represent all such quaternions.


title: Interpretable neural networks based on continuous-valued logic and multicriteria decision operators

abstract: Combining neural networks with continuous logic and multicriteria decision making tools can reduce the black box nature of neural models. In this study, we show that nilpotent logical systems offer an appropriate mathematical framework for a hybridization of continuous nilpotent logic and neural models, helping to improve the interpretability and safety of machine learning. In our concept, perceptrons model soft inequalities; namely membership functions and continuous logical operators. We design the network architecture before training, using continuous logical operators and multicriteria decision tools with given weights working in the hidden layers. Designing the structure appropriately leads to a drastic reduction in the number of parameters to be learned. The theoretical basis offers a straightforward choice of activation functions (the cutting function or its differentiable approximation, the squashing function), and also suggests an explanation to the great success of the rectified linear unit (ReLU). In this study, we focus on the architecture of a hybrid model and introduce the building blocks for future application in deep neural networks. The concept is illustrated with some toy examples taken from an extended version of the tensorflow playground.


title: Optimality of feedback control for qubit purification under inefficient measurement

abstract: A quantum system may be purified, i.e., projected into a pure state, faster if one applies feedback operations during the measurement process. However existing results suggest that such an enhancement is only possible when the measurement efficiency exceeds 0.5, which is difficult to achieve experimentally. We address the task of finding the global optimal feedback control for purifying a single qubit in the presence of measurement inefficiency. We use the Bloch vector length, a more physical and practical quantity than purity, to assess the quality of the state, and employ a backward iteration algorithm to find the globally optimal strategy. Our results show that a speedup is available for quantum efficiencies well below 0.5, which opens the possibility of experimental implementation in existing systems.


title: Statistical Analysis of Stationary Solutions of Coupled Nonconvex Nonsmooth Empirical Risk Minimization

abstract: This paper has two main goals: (a) establish several statistical properties---consistency, asymptotic distributions, and convergence rates---of stationary solutions and values of a class of coupled nonconvex and nonsmoothempirical risk minimization problems, and (b) validate these properties by a noisy amplitude-based phase retrieval problem, the latter being of much topical interest.Derived from available data via sampling, these empirical risk minimization problems are the computational workhorse of a population risk model which involves the minimization of an expected value of a random functional. When these minimization problems are nonconvex, the computation of their globally optimal solutions is elusive. Together with the fact that the expectation operator cannot be evaluated for general probability distributions, it becomes necessary to justify whether the stationary solutions of the empirical problems are practical approximations of the stationary solution of the population problem. When these two features, general distribution and nonconvexity, are coupled with nondifferentiability that often renders the problems "non-Clarke regular", the task of the justification becomes challenging. Our work aims to address such a challenge within an algorithm-free setting. The resulting analysis is therefore different from the much of the analysis in the recent literature that is based on local search algorithms. Furthermore, supplementing the classical minimizer-centric analysis, our results offer a first step to close the gap between computational optimization and asymptotic analysis of coupled nonconvex nonsmooth statistical estimation problems, expanding the former with statistical properties of the practically obtained solution and providing the latter with a more practical focus pertaining to computational tractability.


title: Open sets in computability theory and Reverse Mathematics

abstract: To enable the study of open sets in computational approaches to mathematics, lots of extra data and structure on these sets is assumed. For both foundational and mathematical reasons, it is then a natural question, and the subject of this paper, what the influence of this extra data and structure is on the logical and computational properties of basic theorems pertaining to open sets. To answer this question, we study various basic theorems of analysis, like the Baire category, Heine, Heine-Borel, Urysohn, and Tietze theorems, all for open sets given by their (third-order) characteristic functions. Regarding computability theory, the objects claimed to exist by the aforementioned theorems undergo a shift from `computable' to `not computable in any type two functional', following Kleene's S1-S9. Regarding Reverse Mathematics, the latter's so-called Main Question, namely which set existence axioms are necessary for proving a given theorem, does not have a unique or unambiguous answer for the aforementioned theorems, working in Kohlenbach's higher-order framework. A finer study of representations of open sets leads to the new `$\Delta$-functional' which has unique (computational) properties.


title: Kimera: an Open-Source Library for Real-Time Metric-Semantic Localization and Mapping

abstract: We provide an open-source C++ library for real-time metric-semantic visual-inertial Simultaneous Localization And Mapping (SLAM). The library goes beyond existing visual and visual-inertial SLAM libraries (e.g., ORB-SLAM, VINS- Mono, OKVIS, ROVIO) by enabling mesh reconstruction and semantic labeling in 3D. Kimera is designed with modularity in mind and has four key components: a visual-inertial odometry (VIO) module for fast and accurate state estimation, a robust pose graph optimizer for global trajectory estimation, a lightweight 3D mesher module for fast mesh reconstruction, and a dense 3D metric-semantic reconstruction module. The modules can be run in isolation or in combination, hence Kimera can easily fall back to a state-of-the-art VIO or a full SLAM system. Kimera runs in real-time on a CPU and produces a 3D metric-semantic mesh from semantically labeled images, which can be obtained by modern deep learning methods. We hope that the flexibility, computational efficiency, robustness, and accuracy afforded by Kimera will build a solid basis for future metric-semantic SLAM and perception research, and will allow researchers across multiple areas (e.g., VIO, SLAM, 3D reconstruction, segmentation) to benchmark and prototype their own efforts without having to start from scratch.


title: Different single photon response of wide and narrow superconducting MoSi strips

abstract: The photon count rate of superconducting single photon detectors made of MoSi films shaped as a $2\,\mu$m-wide strip and a 115 nm-wide meander stripline is studied experimentally as a function of the dc biasing current at different values of the perpendicular magnetic field. For the wide strip a crossover current $I_\textrm{cross}$ is observed, below which the photon count rate increases with increasing magnetic field and above which it decreases. This behavior contrasts with the narrow MoSi meander for which no crossover current is observed, thus suggesting different photon detection mechanisms in the wide and narrow strips. Namely, we argue that in the wide strip the absorbed photon destroys superconductivity locally via the vortex-antivortex mechanism for the emergence of resistance while in the narrow meander superconductivity is destroyed across the whole stripline, forming a hot belt. Accordingly, the different photon detection mechanisms associated with vortices and the hot belt stipulate the qualitative difference in the dependence of the photon count rate on the magnetic field.


title: Hunt for dark subhalos in the galactic stellar field using computer vision

abstract: The lack of tangible evidence for non-gravitational interactions between dark and visible sectors drives the need for exploring new avenues of inferring dark matter properties through purely gravitational probes. In particular, addressing small-scale distribution of dark matter could lead to valuable new insights into its particle nature, either confirming predictions of the collisionless cold dark matter hypothesis or favouring models with suppressed small-scale matter power spectrum. In this work we propose a novel machine learning approach for constraining the abundance of galactic dark matter subhaloes through the analysis of Milky Way's stellar field that has been only recently mapped with sufficient coverage thanks to the Gaia mission. Our method is based on convolutional neural networks which represent a powerful tool for identifying characteristic perturbations in spatial maps of stellar number density and velocity distribution moments. For generating the training data we develop a robust and computationally efficient algorithm, capable of generating mock stellar fields from an arbitrary underlying phase-space distribution of stars. By preforming a preliminary study of the outlined approach on synthetic datasets we demonstrate that sensitivities down to (or even below) $10^8 M_\odot$ could be reached. Furthermore, our results show that the accuracy of the advocated technique crucially depends on the kinematic properties of mapped stars and could be further improved by applying it to abundant stellar populations with particularly low velocity dispersion, such as the galactic thin disc stars.


title: A Riemann-Hilbert approach to the lower tail of the KPZ equation

abstract: Fredholm determinants associated to deformations of the Airy kernel are closely connected to the solution to the Kardar-Parisi-Zhang (KPZ) equation with narrow wedge initial data, and they also appear as largest particle distribution in models of positive-temperature free fermions. We show that logarithmic derivatives of the Fredholm determinants can be expressed in terms of a 2x2 Riemann-Hilbert problem, and we use this to derive asymptotics for the Fredholm determinants. As an application of our result, we derive precise lower tail asymptotics for the solution of the KPZ equation with narrow wedge initial data, refining recent results by Corwin and Ghosal.


title: Pomeron/glueball and odderon/oddball trajectories

abstract: We predict glueball/oddball resonances lying on the pomeron/odderon trajectories. A simple new form of the trajectories, with threshold and asymptotic behaviour required by analyticity and unitarity, is proposed. The parameters of these (pomeron and odderon) trajectories are fitted to the data on high-energy elastic proton-proton scattering. The fitted trajectories are extrapolated to the resonance region to predict masses and widths of glueballs and oddballs. The (pomeron and odderon) trajectories may be used to calculate processes of central exclusive diffraction (CED).


title: Confinement and activity regulate bacterial motion in porous media

abstract: Understanding how bacteria move in porous media is critical to applications in healthcare, agriculture, environmental remediation, and chemical sensing. Recent work has demonstrated that E. coli, which moves by run-and-tumble dynamics in a homogeneous medium, exhibits a new form of motility when confined in a disordered porous medium: hopping-and-trapping motility, in which cells perform rapid, directed hops punctuated by intervals of slow, undirected trapping. Here, we use direct visualization to shed light on how these processes depend on pore-scale confinement and cellular activity. We find that hopping is determined by pore-scale confinement, and is independent of cellular activity; by contrast, trapping is determined by the competition between pore-scale confinement and cellular activity, as predicted by an entropic trapping model. These results thus help to elucidate the factors that regulate bacterial motion in porous media, and could help aid the development of new models of motility in heterogeneous environments.

categories: cond-mat.soft cond-mat.dis-nn cond-mat.mes-hall cond-mat.stat-mech physics.bio-ph


title: Machine learning design of a trapped-ion quantum spin simulator

abstract: Trapped ions have emerged as one of the highest quality platforms for the quantum simulation of interacting spin models of interest to various fields of physics. In such simulators, two effective spins can be made to interact with arbitrary strengths by coupling to the collective vibrational or phonon states of ions, controlled by precisely tuned laser beams. However, the task of determining laser control parameters required for a given spin-spin interaction graph is a type of inverse problem, which can be highly mathematically complex. In this paper, we adapt a modern machine learning technique developed for similar inverse problems to the task of finding the laser control parameters for a number of interaction graphs. We demonstrate that typical graphs, forming regular lattices of interest to physicists, can easily be produced for up to 50 ions using a single GPU workstation. The scaling of the machine learning method suggests that this can be expanded to hundreds of ions with moderate additional computational effort.


title: mfEGRA: Multifidelity Efficient Global Reliability Analysis through Active Learning for Failure Boundary Location

abstract: This paper develops mfEGRA, a multifidelity active learning method using data-driven adaptively refined surrogates for failure boundary location in reliability analysis. This work addresses the issue of prohibitive cost of reliability analysis using Monte Carlo sampling for expensive-to-evaluate high-fidelity models by using cheaper-to-evaluate approximations of the high-fidelity model. The method builds on the Efficient Global Reliability Analysis (EGRA) method, which is a surrogate-based method that uses adaptive sampling for refining Gaussian process surrogates for failure boundary location using a single-fidelity model. Our method introduces a two-stage adaptive sampling criterion that uses a multifidelity Gaussian process surrogate to leverage multiple information sources with different fidelities. The method combines expected feasibility criterion from EGRA with one-step lookahead information gain to refine the surrogate around the failure boundary. The computational savings from mfEGRA depends on the discrepancy between the different models, and the relative cost of evaluating the different models as compared to the high-fidelity model. We show that accurate estimation of reliability using mfEGRA leads to computational savings of $\sim$46% for an analytic multimodal test problem and 24% for a three-dimensional acoustic horn problem, when compared to single-fidelity EGRA. We also show the effect of using a priori drawn Monte Carlo samples in the implementation for the acoustic horn problem, where mfEGRA leads to computational savings of 45% for the three-dimensional case and 48% for a rarer event four-dimensional case as compared to single-fidelity EGRA.


title: Predicting popularity of EV charging infrastructure from GIS data

abstract: The availability of charging infrastructure is essential for large-scale adoption of electric vehicles (EV). Charging patterns and the utilization of infrastructure have consequences not only for the energy demand, loading local power grids but influence the economic returns, parking policies and further adoption of EVs. We develop a data-driven approach that is exploiting predictors compiled from GIS data describing the urban context and urban activities near charging infrastructure to explore correlations with a comprehensive set of indicators measuring the performance of charging infrastructure. The best fit was identified for the size of the unique group of visitors (popularity) attracted by the charging infrastructure. Consecutively, charging infrastructure is ranked by popularity. The question of whether or not a given charging spot belongs to the top tier is posed as a binary classification problem and predictive performance of logistic regression regularized with an l-1 penalty, random forests and gradient boosted regression trees is evaluated. Obtained results indicate that the collected predictors contain information that can be used to predict the popularity of charging infrastructure. The significance of predictors and how they are linked with the popularity are explored as well. The proposed methodology can be used to inform charging infrastructure deployment strategies.


title: The ground state solution for Kirchhoff-Schrodinger equations with singular exponential nonlinearities in R^4

abstract: In this paper, we consider Kirchhoff-Schrodinger equations with singular exponential nonlinearities in R^4,using singular Adams inequality and variational techniques, we get the existence of ground state solutions. Moreover, we also get the same results without the Ambrosetti-Rabinowitz (AR) condition.


title: Data-Driven Reachable Set Computation using Adaptive Gaussian Process Classification and Monte Carlo Methods

abstract: We present two data-driven methods for estimating reachable sets with probabilistic guarantees. Both methods make use of a probabilistic formulation allowing for a formal definition of a data-driven reachable set approximation that is correct in a probabilistic sense. The first method recasts the reachability problem as a binary classification problem, using a Gaussian process classifier to represent the reachable set. The quantified uncertainty of the Gaussian process model allows for an adaptive approach to the selection of new sample points. The second method uses a Monte Carlo sampling approach to compute an interval-based approximation of the reachable set. This method comes with a guarantee of probabilistic correctness, and an explicit bound on the number of sample points needed to achieve a desired accuracy and confidence. Each method is illustrated with a numerical example.


title: Bayesian analysis of dynamic binary data: A simulation study and application to economic index SP

abstract: It is proposed in the literature that in some complicated problems maximum likelihood estimates (MLE) are not suitable or even do not exist. An alternative to MLE for estimation of the parameters is the Bayesian method. The Markov chain Monte Carlo (MCMC) simulation procedure is designed to fit Bayesian models. Bayesian method like classical method (MLE) has advantages and disadvantages. One of the advantages of Bayesian method over MLE method is the ability of saving the information included in past data through the posterior distributions of the model parameters to be used for modelling future data. In this article we investigate the performance of Bayesian method in modelling dynamic binary data when the data are growing over time and individuals.


title: The Hartree and Hartree-Fock equations in Lebesgue $L^p$ and Fourier-Lebesgue $\hat{L}^p$ spaces

abstract: We establish some local and global well-posedness for Hartree-Fock equations of $N$ particles (HFP) with Cauchy data in Lebesgue spaces $L^p \cap L^2 $ for $1\leq p \leq \infty$. Similar results are proven for fractional HFP in Fourier-Lebesgue spaces $ \hat{L}^p \cap L^2 \ (1\leq p \leq \infty).$ On the other hand, we show that the Cauchy problem for HFP is ill-posed if we simply work in $\hat{L}^p \ (2<p\leq \infty).$ Analogue results hold for reduced HFP. In the process, we prove the boundedeness of various trilinear estimates for Hartree type non linearity in these spaces which may be of independent interest. As a consequence, we get natural $L^p$ and $\hat{L}^p$ extension of classical well-posedness theories of Hartree and Hartree-Fock equations with Cauchy data in just $L^2-$based Sobolev spaces.


title: Pomeron-pomeron scattering

abstract: Central exclusive diffractive (CED) production of meson resonances potentially is a factory producing new particles, in particular a glueball. The produced resonances lie in trajectories with vacuum quantum numbers, essentially on the pomeron trajectory. A tower of resonance recurrences, the production cross section and the resonances widths are predicted. A new feature is the form of the non-linear pomeron trajectory, producing resonances (glueballs) with increasing widths. At LHC energies, in the nearly forward direction the $t$-channel both in elastic, single or double diffraction dissociation as well as in CED is dominated by pomeron exchange (the role of secondary trajectories is negligible, however a small contribution from the odderon may be present).


title: Detecting Magnetic Fields in Exoplanets with Spectropolarimetry of the Helium Line at 1083 nm

abstract: The magnetic fields of the solar system planets provide valuable insights into the planets' interiors and can have dramatic consequences for the evolution of their atmospheres and interaction with the solar wind. However, we have little direct knowledge of magnetic fields in exoplanets. Here we present a method for detecting magnetic fields in the atmospheres of close-in exoplanets based on spectropolarimetric transit observations at the wavelength of the helium line at 1083 nm. This methodology has been successfully applied before for exploring magnetic fields in solar coronal filaments. Strong absorption signatures (transit depths on the order of a few percent) in the 1083 nm line have recently been observed for several close-in exoplanets. We show that in the conditions in these escaping atmospheres, metastable helium atoms should be optically pumped by the starlight and, for field strengths more than a few $\times 10^{-4}$ G, should align with the magnetic field. This results in linearly polarized absorption at 1083 nm that traces the field direction (the Hanle effect), which we explore by both analytic computation and with the Hazel numerical code. The linear polarization $\sqrt{Q^2+U^2}/I$ ranges from $\sim 10^{-3}$ in optimistic cases down to a few $\times 10^{-5}$ for particularly unfavorable cases, with very weak dependence on field strength. The line-of-sight component of the field results in a slight circular polarization (the Zeeman effect), also reaching $V/I\sim {\rm few}\times 10^{-5}(B_\parallel/10\,{\rm G})$. We discuss the detectability of these signals with current (SPIRou) and future (extremely large telescope) high-resolution infrared spectropolarimeters, and we briefly comment on possible sources of astrophysical contamination.


title: Boosting Local Causal Discovery in High-Dimensional Expression Data

abstract: We study the performance of Local Causal Discovery (LCD), a simple and efficient constraint-based method for causal discovery, in predicting causal effects in large-scale gene expression data. We construct practical estimators specific to the high-dimensional regime. Inspired by the ICP algorithm, we use an optional preselection method and two different statistical tests. Empirically, the resulting LCD estimator is seen to closely approach the accuracy of ICP, the state-of-the-art method, while it is algorithmically simpler and computationally more efficient.


title: Predicting Phenotypes from Brain Connection Structure

abstract: This article focuses on the problem of predicting a response variable based on a network-valued predictor. Our motivation is the development of interpretable and accurate predictive models for cognitive traits and neuro-psychiatric disorders based on an individual's brain connection network (connectome). Current methods reduce the complex, high dimensional brain network into low-dimensional pre-specified features prior to applying standard predictive algorithms. These methods are sensitive to feature choice and inevitably discard important information. Instead, we propose a nonparametric Bayes class of models that utilize the entire adjacency matrix defining brain region connections to adaptively detect predictive algorithms, while maintaining interpretability. The Bayesian Connectomics (BaCon) model class utilizes Poisson-Dirichlet processes to find a lower-dimensional, bidirectional (covariate, subject) pattern in the adjacency matrix. The small n, large p problem is transformed into a "small n, small q" problem, facilitating an effective stochastic search of the predictors. A spike-and-slab prior for the cluster predictors strikes a balance between regression model parsimony and flexibility, resulting in improved inferences and test case predictions. We describe basic properties of the BaCon model and develop efficient algorithms for posterior computation. The resulting methods are found to outperform existing approaches and applied to a creative reasoning data set.


title: Predicting publication productivity for researchers: A latent variable model

abstract: This study provided a model for the publication dynamics of researchers, which is based on the relationship between the publication productivity of researchers and two covariates: time and historical publication quantity. The relationship allows to estimate the latent variable the publication creativity of researchers. The variable is applied to the prediction of publication productivity for researchers. The statistical significance of the relationship is validated by the high quality dblp dataset. The effectiveness of the model is testified on the dataset by the fine fittings on the quantitative distribution of researchers' publications, the evolutionary trend of their publication productivity, and the occurrence of publication events. Due to its nature of regression, the model has the potential to be extended for assessing the confidence level of prediction results, and thus has clear applicability to empirical research.


title: Mullins-Sekerka as the Wasserstein flow of the perimeter

abstract: We prove the convergence of an implicit time discretization for the one-phase Mullins-Sekerka equation, possibly with additional non-local repulsion, proposed in [F. Otto, Arch. Rational Mech. Anal. 141 (1998) 63--103]. Our simple argument shows that the limit satisfies the equation in a distributional sense as well as an optimal energy-dissipation relation. The proof combines arguments from optimal transport, gradient flows & minimizing movements, and basic geometric measure theory.


title: REMIND Your Neural Network to Prevent Catastrophic Forgetting

abstract: People learn throughout life. However, incrementally updating conventional neural networks leads to catastrophic forgetting. A common remedy is replay, which is inspired by how the brain consolidates memory. Replay involves fine-tuning a network on a mixture of new and old instances. While there is neuroscientific evidence that the brain replays compressed memories, existing methods for convolutional networks replay raw images. Here, we propose REMIND, a brain-inspired approach that enables efficient replay with compressed representations. REMIND is trained in an online manner, meaning it learns one example at a time, which is closer to how humans learn. Under the same constraints, REMIND outperforms other methods for incremental class learning on the ImageNet ILSVRC-2012 dataset. We probe REMIND's robustness to data ordering schemes known to induce catastrophic forgetting. We demonstrate REMIND's generality by pioneering online learning for Visual Question Answering (VQA).


title: Strichartz estimates and Strauss conjecture on non-trapping asymptotically hyperbolic manifolds

abstract: We prove global-in-time Strichartz estimates for the shifted wave equations on non-trapping asymptotically hyperbolic manifolds. The key tools are the spectral measure estimates from \cite{CH2} and arguments borrowed from \cite{HZ, Zhang}. As an application, we prove the small data global existence for any power $p\in(1, 1+\frac{4}{n-1})$ for the shifted wave equation in this setting, involving nonlinearities of the form $\pm|u|^p$ or $\pm|u|^{p-1}u$, which answers partially an open question raised in \cite{SSW}.


title: Patterning Sn-based EUV resists with low-energy electrons

abstract: Extreme Ultraviolet (EUV) lithography is the newest technology that will be used in the semiconductor industry for printing circuitry in the sub-20 nm scale. Low-energy electrons (LEEs) produced upon illumination of resist materials with EUV photons (92 eV) play a central role in the formation of the nanopatterns. However, up to now the details of this process are not well understood. In this work, a novel experimental approach that combines Low-Energy Electron Microscopy (LEEM), Electron Energy Loss Spectroscopy (EELS), and Atomic Force Microscopy (AFM) is used to study changes induced by electrons in the 0-40 eV range in thin films of molecular organometallic EUV resists known as tin-oxo cages. LEEM-EELS spectroscopic experiments were used to detect surface charging upon electron exposure and to estimate the electron landing energy. AFM post-exposure analyses revealed that irradiation of the resist with LEEs leads to the densification of the resist layer associated to carbon loss. The same chemical processes that yield densification render the solubility change responsible for the pattern formation in the lithographic application. Remarkably, electrons as low as 1.2 eV are able to induce chemical reactions in the Sn-based resist. Based on the thickness profiles resulting from LEE exposures in the 3-48 mC/cm 2 dose range, a simplified reaction model is proposed where the resist undergoes sequential chemical reactions, yielding first a sparsely cross-linked network, followed by the formation of a denser cross-linked network. This model allows us to estimate a maximum reaction volume on the initial material of 0.15 nm 3 per incident electron in the energy range studied, which means that less than 10 LEEs per molecule on average are needed to turn the material insoluble and thus render a pattern.


title: Numerical investigation on the Hill's type lunar problem with homogeneous potential

abstract: We consider the planar Hill's lunar problem with a homogeneous gravitational potential. The investigation of the system is twofold. First, the starting conditions of the trajectories are classified into three classes, that is bounded, escaping, and collisional. Second, we study the no-return property of the Lagrange point $L_2$ and we observe that the escaping trajectories are scattered exponentially. Moreover, it is seen that in the supercritical case, with $\alpha \geq 2$, the basin boundaries are smooth. On the other hand, in the subcritical case, with $\alpha < 2$ the boundaries between the different types of basins exhibit fractal properties.


title: Automated Isolation for White-box Test Generation

abstract: Context. White-box test generation is a technique used for automatically selecting test inputs using only the source or binary code. However, such techniques encounter challenges when applying them to complex programs. One of the main challenges is handling the dependencies of the unit under test. Objective. Without proper actions, generated tests cannot cover all parts of the source code, or calling the dependencies may cause unexpected side effects (e.g., file system or network access). These issues should be tackled while maintaining the advantages of white-box test generation. Method. In this paper, we present an automated source code transformation approach tackling the dependency issue for white-box test generation. This technique isolates the test execution by creating a parameterized sandbox wrapped around the transformed unit. We implemented the approach in a ready-to-use tool using Microsoft Pex as a test generator, and evaluated it on 10 open-source projects from GitHub having more than 38.000 lines of code in total. Results. The results from the evaluation indicate that if the lack of isolation hinders white-box test generation, then our approach is able to help: it increases the code coverage reached by the automatically generated test, while it reduces unwanted side effects. Also, our results act as a unique baseline for the test generation performance of Microsoft Pex on open-source projects. Conclusion. Based on the results, our source code transformations might serve well for alleviating the isolation problem in white-box test generation as it increases the coverage reached in such situations, while maintaining the practical applicability of the tests generated on the isolated code.


title: Adaptive Krylov-Type Time Integration Methods

abstract: The Rosenbrock-Krylov family of time integration schemes is an extension of Rosenbrock-W methods that employs a specific Krylov based approximation of the linear system solutions arising within each stage of the integrator. This work proposes an extension of Rosenbrock-Krylov methods to address stability questions which arise for methods making use of inexact linear system solution strategies. Two approaches for improving the stability and efficiency of Rosenbrock-Krylov methods are proposed, one through direct control of linear system residuals and the second through a novel extension of the underlying Krylov space to include stage right hand side vectors. Rosenbrock-Krylov methods employing the new approaches show a substantial improvement in computational efficiency relative to prior implementations.


title: The lost boarding pass, and other practical problems

abstract: The reader is reminded of several puzzles involving randomness. These may be ill-posed, and if well-posed there is sometimes a solution that uses probabilistic intuition in a special way. Various examples are presented including the well known problem of the lost boarding pass: what is the probability that the last passenger boarding a fully booked plane sits in the assigned seat if the first passenger has occupied a randomly chosen seat? This problem, and its striking answer of $\frac12$, has attracted a good deal of attention since around 2000. We review elementary solutions to this, and to the more general problem of finding the probability the $m$th passenger sits in the assigned seat when in the presence of some number $k$ of passengers with lost boarding passes. A simple proof is presented of the independence of the occupancy status of different seats, and a connection to the Poisson--Dirichlet distribution is mentioned.


title: Optimising energy and overhead for large parameter space simulations

abstract: Many systems require optimisation over multiple objectives, where objectives are characteristics of the system such as energy consumed or increase in time to perform the work. Optimisation is performed by selecting the `best' set of input parameters to elicit the desired objectives. However, the parameter search space can often be far larger than can be searched in a reasonable time. Additionally, the objectives are often mutually exclusive -- leading to a decision being made as to which objective is more important or optimising over a combination of the objectives. This work is an application of a Genetic Algorithm to identify the Pareto frontier for finding the optimal parameter sets for all combinations of objectives. A Pareto frontier can be used to identify the sets of optimal parameters for which each is the `best' for a given combination of objectives -- thus allowing decisions to be made with full knowledge. We demonstrate this approach for the HTC-Sim simulation system in the case where a Reinforcement Learning scheduler is tuned for the two objectives of energy consumption and task overhead. Demonstrating that this approach can reduce the energy consumed by ~36% over previously published work without significantly increasing the overhead.


title: Fine-Grained Analysis of Propaganda in News Articles

abstract: Propaganda aims at influencing people's mindset with the purpose of advancing a specific agenda. Previous work has addressed propaganda detection at the document level, typically labelling all articles from a propagandistic news outlet as propaganda. Such noisy gold labels inevitably affect the quality of any learning system trained on them. A further issue with most existing systems is the lack of explainability. To overcome these limitations, we propose a novel task: performing fine-grained analysis of texts by detecting all fragments that contain propaganda techniques as well as their type. In particular, we create a corpus of news articles manually annotated at the fragment level with eighteen propaganda techniques and we propose a suitable evaluation measure. We further design a novel multi-granularity neural network, and we show that it outperforms several strong BERT-based baselines.


title: Self-similar dynamics of order parameter fluctuations in pump-probe experiments

abstract: Upon excitation by a laser pulse, broken-symmetry phases of a wide variety of solids demonstrate similar order parameter dynamics characterized by a dramatic slowing down of relaxation for stronger pump fluences. Motivated by this recurrent phenomenology, we develop a simple non-perturbative effective model of dynamics of collective bosonic excitations in pump-probe experiments. We find that as the system recovers after photoexcitation, it shows universal prethermalized dynamics manifesting a power-law, as opposed to exponential, relaxation, explaining the slowing down of the recovery process. For strong quenches, long-wavelength over-populated transverse modes dominate the long-time dynamics; their distribution function exhibits universal scaling in time and space, whose universal exponents can be computed analytically. Our model offers a unifying description of order parameter fluctuations in a regime far from equilibrium, and our predictions can be tested with available time-resolved techniques.


title: FIS-GAN: GAN with Flow-based Importance Sampling

abstract: Generative Adversarial Networks (GAN) training process, in most cases, apply Uniform or Gaussian sampling methods in the latent space, which probably spends most of the computation on examples that can be properly handled and easy to generate. Theoretically, importance sampling speeds up stochastic optimization in supervised learning by prioritizing training examples. In this paper, we explore the possibility of adapting importance sampling into adversarial learning. We use importance sampling to replace Uniform and Gaussian sampling methods in the latent space and employ normalizing flow to approximate latent space posterior distribution by density estimation. Empirically, results on MNIST and Fashion-MNIST demonstrate that our method significantly accelerates GAN's optimization while retaining visual fidelity in generated samples.


title: Analysis of tensor methods for stochastic models of gene regulatory networks

abstract: The tensor-structured parametric analysis (TPA) has been recently developed for simulating and analysing stochastic behaviours of gene regulatory networks [Liao et. al., 2015]. The method employs the Fokker-Planck approximation of the chemical master equation, and uses the Quantized Tensor Train (QTT) format, as a low-parametric tensor-structured representation of classical matrices and vectors, to approximate the high-dimensional stationary probability distribution. This paper presents a detailed error analysis of all approximation steps of the TPA regarding validity and accuracy, including modelling error, artificial boundary error, discretization error, tensor rounding error, and algebraic error. The error analysis is illustrated using computational examples, including the death-birth process and a 50-dimensional isomerization reaction chain.


title: Representing Dependencies in Event Structures

abstract: Event structures where the causality may explicitly change during a computation have recently gained the stage. In this kind of event structures the changes in the set of the causes of an event are triggered by modifiers that may add or remove dependencies, thus making the happening of an event contextual. Still the focus is always on the dependencies of the event. In this paper we promote the idea that the context determined by the modifiers plays a major role, and the context itself determines not only the causes but also what causality should be. Modifiers are then used to understand when an event (or a set of events) can be added to a configuration, together with a set of events modeling dependencies, which will play a less important role. We show that most of the notions of Event Structure presented in literature can be translated into this new kind of event structure, preserving the main notion, namely the one of configuration.


title: Drop fragmentation by laser-pulse impact

abstract: We study the fragmentation of a liquid drop that is hit by a laser pulse. The drop expands into a thin sheet that breaks by the radial expulsion of ligaments from its rim and the nucleation and growth of holes on the sheet. By combining experimental data from two liquid systems with vastly different time- and length scales we show how the early-time laser-matter interaction affects the late-time fragmentation. We identify two Rayleigh--Taylor instabilities of different origins as the prime cause of the fragmentation and derive scaling laws for the characteristic breakup time and wavenumber. The final web of ligaments results from a subtle interplay between these instabilities and deterministic modulations of the local sheet thickness, which originate from the drop deformation dynamics and spatial variations in the laser-beam profile.


title: Understanding how a falling ball chain can be speeded up by impact onto a surface

abstract: When a falling ball chain strikes a surface, a tension is created that pulls the chain downward. This causes a downward acceleration that is larger than free-fall, which has been observed by recent experiments. Here a theoretical description of this surprising phenomenon is developed. The equation of motion for the falling chain is derived, and then solved for a general form of the tension. the size of the tension needed to produce the observed motion is relatively small and is explained here as coming from the rotation of a link just above where the chain collides with the surface. This simple model is used to calculate the size of the tension in terms of physically measurable quantities: the length and width of link, the maximum bending angle at a junction, the inclination angle of the surface, and the coefficients of friction and restitution between the chain and the surface. The model's predictions agree with the results of current experiments. New experiments are proposed that can test the model.


title: Finite simple exceptional groups of Lie type in which all the subgroups of odd index are pronormal

abstract: A subgroup $H$ of a group $G$ is said to be pronormal in $G$ if $H$ and $H^g$ are conjugate in $\langle H, H^g \rangle$ for every $g \in G$. In this paper we classify finite simple groups $E_6(q)$ and ${}^2E_6(q)$ in which all the subgroups of odd index are pronormal. Thus, we complete a classification of finite simple exceptional groups of Lie type in which all the subgroups of odd index are pronormal.


title: Local Langlands correspondence for the twisted exterior and symmetric square $\epsilon$-factors of $\textrm{GL}_n$

abstract: Let $F$ be a non-Archimedean local field. Let $\mathcal{A}_n(F)$ be the set of equivalence classes of irreducible admissible representations of $\textrm{GL}_n(F)$, and $\mathcal{G}_n(F)$ be the set of equivalence classes of n-dimensional Frobenius semisimple Weil-Deligne representations of $W'_F$. The local Langlands correspondence(LLC) establishes the reciprocity maps $\textrm{Rec}_{n,F}: \mathcal{A}_n(F)\longrightarrow \mathcal{G}_n(F)$ , satisfying some nice properties. An important invariant under this correspondence is the L- and $\epsilon$-factors. This is also expected to be true under parallel compositions with a complex analytic representations of $\textrm{GL}_n(\mathbb{C})$. J.W. Cogdell, F. Shahidi, and T.-L. Tsai proved the equality of the symmetric and exterior square L- and $\epsilon$-factors [7] in 2017. But the twisted symmetric and exterior square L- and $\epsilon$-factor are new and very different from the untwisted case. In this paper we will define the twisted symmetric square L- and $\gamma$-factors using $\textrm{GSpin}_{2n+1}$, and establish the equality of the corresponding L- and $\epsilon$-factors. We will first reduce the problem to the analytic stability of their $\gamma$-factors for supercuspidal representations, then prove the supercuspidal stability by establishing general asymptotic expansions of partial Bessel function following the ideas in [7].


title: Joint Image and Depth Estimation with Mask-Based Lensless Cameras

abstract: Mask-based lensless cameras replace the lens of a conventional camera with a custom mask. These cameras can potentially be very thin and even flexible. Recently, it has been demonstrated that such mask-based cameras can recover light intensity and depth information of a scene. Existing depth recovery algorithms either assume that the scene consists of a small number of depth planes or solve a sparse recovery problem over a large 3D volume. Both these approaches fail to recover the scenes with large depth variations. In this paper, we propose a new approach for depth estimation based on an alternating gradient descent algorithm that jointly estimates a continuous depth map and light distribution of the unknown scene from its lensless measurements. We present simulation results on image and depth reconstruction for a variety of 3D test scenes. A comparison between the proposed algorithm and other method shows that our algorithm is more robust for natural scenes with a large range of depths. We built a prototype lensless camera and present experimental results for reconstruction of intensity and depth maps of different real objects.


title: 3D Scene Graph: A Structure for Unified Semantics, 3D Space, and Camera

abstract: A comprehensive semantic understanding of a scene is important for many applications - but in what space should diverse semantic information (e.g., objects, scene categories, material types, texture, etc.) be grounded and what should be its structure? Aspiring to have one unified structure that hosts diverse types of semantics, we follow the Scene Graph paradigm in 3D, generating a 3D Scene Graph. Given a 3D mesh and registered panoramic images, we construct a graph that spans the entire building and includes semantics on objects (e.g., class, material, and other attributes), rooms (e.g., scene category, volume, etc.) and cameras (e.g., location, etc.), as well as the relationships among these entities. However, this process is prohibitively labor heavy if done manually. To alleviate this we devise a semi-automatic framework that employs existing detection methods and enhances them using two main constraints: I. framing of query images sampled on panoramas to maximize the performance of 2D detectors, and II. multi-view consistency enforcement across 2D detections that originate in different camera locations.


title: Computational models for active matter

abstract: A variety of computational models have been developed to describe active matter at different length and time scales. The diversity of the methods and the challenges in modeling active matter---ranging from molecular motors and cytoskeletal filaments over artificial and biological swimmers on microscopic to groups of animals on macroscopic scales---mainly originate from their out-of-equilibrium character, multiscale nature, nonlinearity, and multibody interactions. In the present review, various modeling approaches and numerical techniques are addressed, compared, and differentiated to illuminate the innovations and current challenges in understanding active matter. The complexity increases from minimal microscopic models of dry active matter toward microscopic models of active matter in fluids. Complementary, coarse-grained descriptions and continuum models are elucidated. Microscopic details are often relevant and strongly affect collective behaviors, which implies that the selection of a proper level of modeling is a delicate choice, with simple models emphasizing universal properties and detailed models capturing specific features. Finally, current approaches to further advance the existing models and techniques to cope with real-world applications, such as complex media and biological environments, are discussed.


title: A Morawetz inequality for gravity-capillary water waves at low Bond number

abstract: This paper is devoted to the 2D gravity-capillary water waves equations in their Hamiltonian formulation, addressing the general question of proving Morawetz inequalities. We continue the analysis initiated in our previous work, where we have established local energy decay estimates for gravity waves. Here we add surface tension and prove a stronger estimate with a local regularity gain, akin to the smoothing effect for dispersive equations. Our main result holds globally in time and holds for genuinely nonlinear waves, since we are only assuming some very mild uniform Sobolev bounds for the solutions. Furthermore, it is uniform both in the infinite depth limit and the zero surface tension limit.


title: On the Hausdorff dimension of Riemann's non-differentiable function

abstract: Recent findings show that the classical Riemann's non-differentiable function has a physical and geometric nature as the irregular trajectory of a polygonal vortex filament driven by the binormal flow. In this article, we give an upper estimate of its Hausdorff dimension. We also adapt this result to the multifractal setting. To prove these results, we recalculate the asymptotic behavior of Riemann's function around rationals from a novel perspective, underlining its connections with the Talbot effect and Gauss sums, with the hope that it is useful to give a lower bound of its dimension and to answer further geometric questions.


title: K\"ahler manifolds with almost non-negative curvature

abstract: In this paper, we construct local and global solutions to the K\"ahler-Ricci flow from a non-collapsed K\"ahler manifold with curvature bounded from below. Combines with the mollification technique of McLeod-Simon-Topping, we show that the Gromov-Hausdorff limit of sequence of complete noncompact non-collapsed K\"ahler manifolds with orthogonal bisectional curvature and Ricci curvature bounded from below is homeomorphic to a complex manifold. We also use it to study the complex structure of complete K\"ahler manifolds with nonnegative orthogonal bisectional curvature, nonnegative Ricci curvature and maximal volume growth.


title: Probabilistic Successor Representations with Kalman Temporal Differences

abstract: The effectiveness of Reinforcement Learning (RL) depends on an animal's ability to assign credit for rewards to the appropriate preceding stimuli. One aspect of understanding the neural underpinnings of this process involves understanding what sorts of stimulus representations support generalisation. The Successor Representation (SR), which enforces generalisation over states that predict similar outcomes, has become an increasingly popular model in this space of inquiries. Another dimension of credit assignment involves understanding how animals handle uncertainty about learned associations, using probabilistic methods such as Kalman Temporal Differences (KTD). Combining these approaches, we propose using KTD to estimate a distribution over the SR. KTD-SR captures uncertainty about the estimated SR as well as covariances between different long-term predictions. We show that because of this, KTD-SR exhibits partial transition revaluation as humans do in this experiment without additional replay, unlike the standard TD-SR algorithm. We conclude by discussing future applications of the KTD-SR as a model of the interaction between predictive and probabilistic animal reasoning.


title: Compressed Video Action Recognition with Refined Motion Vector

abstract: Although CNN has reached satisfactory performance in image-related tasks, using CNN to process videos is much more challenging due to the enormous size of raw video streams. In this work, we propose to use motion vectors and residuals from modern video compression techniques to effectively learn the representation of the raw frames and greatly remove the temporal redundancy, giving a faster video processing model. Compressed Video Action Recognition(CoViAR) has explored to directly use compressed video to train the deep neural network, where the motion vectors were utilized to present temporal information. However, motion vector is designed for minimizing video size where precise motion information is not obligatory. Compared with optical flow, motion vectors contain noisy and unreliable motion information. Inspired by the mechanism of video compression codecs, we propose an approach to refine the motion vectors where unreliable movement will be removed while temporal information is largely reserved. We prove that replacing the original motion vector with refined one and using the same network as CoViAR has achieved state-of-art performance on the UCF-101 and HMDB-51 with negligible efficiency degrades comparing with original CoViAR.


title: The CEO problem with inter-block memory

abstract: An $n$-dimensional source with memory is observed by $K$ isolated encoders via parallel channels, who compress their observations to transmit to the decoder via noiseless rate-constrained links while leveraging their memory of the past. At each time instant, the decoder receives $K$ new codewords from the observers, combines them with the past received codewords, and produces a minimum-distortion estimate of the latest block of $n$ source symbols. This scenario extends the classical one-shot CEO problem to multiple rounds of communication with communicators maintaining the memory of the past. We extend the Berger-Tung inner and outer bounds to the scenario with inter-block memory, showing that the minimum asymptotically (as $n \to \infty$) achievable sum rate required to achieve a target distortion is bounded by minimal directed mutual information problems. For the Gauss-Markov source observed via $K$ parallel AWGN channels, we show that the inner bound is tight and solve the corresponding minimal directed mutual information problem, thereby establishing the minimum asymptotically achievable sum rate. Finally, we explicitly bound the rate loss due to a lack of communication among the observers; that bound is attained with equality in the case of identical observation channels. The general coding theorem is proved via a new nonasymptotic bound that uses stochastic likelihood coders and whose asymptotic analysis yields an extension of the Berger-Tung inner bound to the causal setting. The analysis of the Gaussian case is facilitated by reversing the channels of the observers.


title: Highly Embedded 8 micron Cores of Star Formation in the Spiral Arms and Filaments of 15 Nearby Disk Galaxies

abstract: Spitzer Space Telescope observations of 15 spiral galaxies show numerous dense cores at 8 microns that are revealed primarily in unsharp mask images. The cores are generally invisible in optical bands because of extinction, and they are also indistinct at 8 microns alone because of contamination by more widespread diffuse emission. Several hundred core positions, magnitudes, and colors from the four IRAC bands are measured and tabulated for each galaxy. The larger galaxies, which tend to have longer and more regular spiral arms, often have their infrared cores aligned along these arms, with additional cores in spiral arm spurs. Galaxies without regular spirals have their cores in more irregular spiral-like filaments, with typically only one or two cores per filament. Nearly every elongated emission feature has 8 micron cores strung out along its length. The occurrence of dense cores in long and thin filaments is reminiscent of filamentary star formation in the solar neighborhood, although on a scale 100 times larger in galaxies. The cores most likely form by gravitational instabilities and cloud agglomeration in the filaments. The simultaneous occurrence of several cores with regular spacings in some spiral arms suggests that in these cases, all of the cores formed at about the same time and the corresponding filaments are young. Total star formation rates for the galaxies correlate with the total embedded stellar masses in the cores with an average ratio corresponding to a possible age between 0.2 Myr and 2 Myr. This suggests that the identified cores are the earliest phase for most star formation.


title: Geometric differentiability of Riemann's non-differentiable function

abstract: Riemann's non-differentiable function is a classic example of a continuous function which is almost nowhere differentiable, and many results concerning its analytic regularity have been shown so far. However, it can also be given a geometric interpretation, so questions on its geometric regularity arise. This point of view is developed in the context of the evolution of vortex filaments, modelled by the Vortex Filament Equation or the binormal flow, in which a generalisation of Riemann's function to the complex plane can be regarded as the trajectory of a particle. The objective of this document is to show that the trajectory represented by its image does not have a tangent anywhere. For that, we discuss several concepts of tangent vectors in view of the set's irregularity.


title: A Note on Alberti's Luzin-Type Theorem for Gradients

abstract: We give a "soft" proof of Alberti's Luzin-type theorem in [1] (G. Alberti, A Lusintype theorem for gradients, J. Funct. Anal. 100 (1991)), using elementary geometric measure theory and topology. Applications to the $C^2$-rectifiability problem are also discussed.


title: Upper Triangularity for Unipotent Representations

abstract: Suppose $G$ is a real reductive group. The determination of the irreducible unitary representations of $G$ is one of the major unsolved problem in representation theory. There is evidence to suggest that every irreducible unitary representation of $G$ can be constructed through a sequence of well-understood operations from a finite set of building blocks, called the unipotent representations. These representations are `attached' (in a certain mysterious sense) to the nilpotent orbits of $G$ on the dual space of its Lie algebra. Inside this finite set is a still smaller set, consisting of the unipotent representations attached to non-induced nilpotent orbits. In this paper, we prove that in many cases this smaller set generates (through a suitable kind of induction) all unipotent representations.


abstract: We investigate R-optimal designs for multi-response regression models with multi-factors, where the random errors in these models are correlated. Several theoretical results are derived for Roptimal designs, including scale invariance, reflection symmetry, line and plane symmetry, and dependence on the covariance matrix of the errors. All the results can be applied to linear and nonlinear models. In addition, an efficient algorithm based on an interior point method is developed for finding R-optimal designs on discrete design spaces. The algorithm is very flexible, and can be applied to any multi-response regression model.


title: A Double Categorical View on Representations of Etendues

abstract: In this paper we introduce a description of ordered groupoids as a particular type of double categories. This enables us to turn Lawson's correspondence between ordered groupoids and left-cancellative categories into a biequivalence. We use this to identify which ordered functors are maps of sites in the sense that they give rise to geometric morphisms between the induced sheaf categories, and establish a Comparison Lemma for maps between Ehresmann sites.


title: Finsler metrics that are both Douglas and generalized Berwald in dimension two

abstract: We proof that in dimension two, a Finsler metric is Douglas and generalized Berwald, if and only if it is Berwald or a Randers metric $\alpha + \beta$, where $\beta$ is closed and is of constant length with respect to $\alpha$.


title: On Inference of Overlapping Coefficients in Two Inverse Lomax Populations

abstract: Overlapping coefficient is a direct measure of similarity between two distributions which is recently becoming very useful. This paper investigates estimation for some well-known measures of overlap, namely Matusita's measure $\rho$, Weitzman's measure $\Delta$ and $\Lambda$ based on Kullback-Leibler. Two estimation methods considered in this study are point estimation and Bayesian approach. Two Inverse Lomax populations with different shape parameters are considered. The bias and mean square error properties of the estimators are studied through a simulation study and a real data example.


title: Improving One-shot NAS by Suppressing the Posterior Fading

abstract: There is a growing interest in automated neural architecture search (NAS). To improve the efficiency of NAS, previous approaches adopt weight sharing method to force all models share the same set of weights. However, it has been observed that a model performing better with shared weights does not necessarily perform better when trained alone. In this paper, we analyse existing weight sharing one-shot NAS approaches from a Bayesian point of view and identify the posterior fading problem, which compromises the effectiveness of shared weights. To alleviate this problem, we present a practical approach to guide the parameter posterior towards its true distribution. Moreover, a hard latency constraint is introduced during the search so that the desired latency can be achieved. The resulted method, namely Posterior Convergent NAS (PC-NAS), achieves state-of-the-art performance under standard GPU latency constraint on ImageNet. In our small search space, our model PC-NAS-S attains 76.8 % top-1 accuracy, 2.1% higher than MobileNetV2 (1.4x) with the same latency. When adopted to the large search space, PC-NAS-L achieves 78.1 % top-1 accuracy within 11ms. The discovered architecture also transfers well to other computer vision applications such as object detection and person re-identification.


title: Using Deep Learning and Machine Learning to Detect Epileptic Seizure with Electroencephalography (EEG) Data

abstract: The prediction of epileptic seizure has always been extremely challenging in medical domain. However, as the development of computer technology, the application of machine learning introduced new ideas for seizure forecasting. Applying machine learning model onto the predication of epileptic seizure could help us obtain a better result and there have been plenty of scientists who have been doing such works so that there are sufficient medical data provided for researchers to do training of machine learning models.


title: Early Prediction of 30-day ICU Re-admissions Using Natural Language Processing and Machine Learning

abstract: ICU readmission is associated with longer hospitalization, mortality and adverse outcomes. An early recognition of ICU re-admission can help prevent patients from worse situation and lower treatment cost. As the abundance of Electronics Health Records (EHR), it is popular to design clinical decision tools with machine learning technique manipulating on healthcare large scale data. We designed data-driven predictive models to estimate the risk of ICU readmission. The discharge summary of each hospital admission was carefully represented by natural language processing techniques. Unified Medical Language System (UMLS) was further used to standardize inconsistency of discharge summaries. 5 machine learning classifiers were adopted to construct predictive models. The best configuration yielded a competitive AUC of 0.748. Our work suggests that natural language processing of discharge summaries is capable to send clinicians warning of unplanned 30-day readmission upon discharge.


title: A theorem of Kalman and minimal state-space realization of Vector Autoregressive Models

abstract: We introduce a concept of $autoregressive$ (AR)state-space realization that could be applied to all transfer functions $\boldsymbol{T}(L)$ with $\boldsymbol{T}(0)$ invertible. We show that a theorem of Kalman implies each Vector Autoregressive model (with exogenous variables) has a minimal $AR$-state-space realization of form $\boldsymbol{y}_t = \sum_{i=1}^p\boldsymbol{H}\boldsymbol{F}^{i-1}\boldsymbol{G}\boldsymbol{x}_{t-i}+\boldsymbol{\epsilon}_t$ where $\boldsymbol{F}$ is a nilpotent Jordan matrix and $\boldsymbol{H}, \boldsymbol{G}$ satisfy certain rank conditions. The case $VARX(1)$ corresponds to reduced-rank regression. Similar to that case, for a fixed Jordan form $\boldsymbol{F}$, $\boldsymbol{H}$ could be estimated by least square as a function of $\boldsymbol{G}$. The likelihood function is a determinant ratio generalizing the Rayleigh quotient. It is unchanged if $\boldsymbol{G}$ is replaced by $\boldsymbol{S}\boldsymbol{G}$ for an invertible matrix $\boldsymbol{S}$ commuting with $\boldsymbol{F}$. Using this invariant property, the search space for maximum likelihood estimate could be constrained to equivalent classes of matrices satisfying a number of orthogonal relations, extending the results in reduced-rank analysis. Our results could be considered a multi-lag canonical-correlation-analysis. The method considered here provides a solution in the general case to the polynomial product regression model of Velu et. al. We provide estimation examples. We also explore how the estimates vary with different Jordan matrix configurations and discuss methods to select a configuration. Our approach could provide an important dimensional reduction technique with potential applications in time series analysis and linear system identification. In the appendix, we link the reduced configuration space of $\boldsymbol{G}$ with a geometric object called a vector bundle.

categories: math.ST cs.SY eess.SY math.AG q-fin.ST stat.TH


title: Initial Characterization of Interstellar Comet 2I/2019 Q4 (Borisov)

abstract: We present initial observations of the interstellar body 2I/(2019 Q4) Borisov taken to determine its nature prior to the perihelion in 2019 December. Images from the Nordic Optical Telescope show a prominent, morphologically stable dust coma and tail. The dust cross-section within 15,000 km of the nucleus averages 130 sq. km (assuming geometric albedo 0.1) and increases by about 1 percent per day. If sustained, this rate indicates that the comet has been active for 100 days prior to the observations. Cometary activity thus started in 2019 June, at which time C/Borisov was at 4.5 AU from the Sun, a typical distance for the onset of water ice sublimation in comets. The dust optical colors, B-V = 0.80+/-0.05, V-R = 0.47+/-0.03 and R-I = 0.49+/-0.05 are identical to those of a sample of (solar system) long-period comets. The colors are similar to those of 1I/(2017 U1) 'Oumuamua, indicating a lack of the ultrared matter that is common in the Kuiper belt, on both interstellar objects. The effective size of the dust particles is estimated as 100 $\mu$m, based on the length of the dust tail and the 100 day lifetime. With this size, the ejected dust mass is of order 1.3e7 kg and the current dust mass loss rate 2 kg/s. We set an upper limit to the nucleus radius using photometry of <3.8 km (again for albedo 0.1). However, we use a statistical argument to show that the nucleus must be much smaller, likely a few hundred meters in radius.


title: Rethinking Kernel Methods for Node Representation Learning on Graphs

abstract: Graph kernels are kernel methods measuring graph similarity and serve as a standard tool for graph classification. However, the use of kernel methods for node classification, which is a related problem to graph representation learning, is still ill-posed and the state-of-the-art methods are heavily based on heuristics. Here, we present a novel theoretical kernel-based framework for node classification that can bridge the gap between these two representation learning problems on graphs. Our approach is motivated by graph kernel methodology but extended to learn the node representations capturing the structural information in a graph. We theoretically show that our formulation is as powerful as any positive semidefinite kernels. To efficiently learn the kernel, we propose a novel mechanism for node feature aggregation and a data-driven similarity metric employed during the training phase. More importantly, our framework is flexible and complementary to other graph-based deep learning models, e.g., Graph Convolutional Networks (GCNs). We empirically evaluate our approach on a number of standard node classification benchmarks, and demonstrate that our model sets the new state of the art.


title: On the Compactness of Hamiltonian Stationary Lagrangian Surfaces in K{\"a}hler Surfaces

abstract: We prove a bubble tree convergence theorem for a sequence of closed Hamiltonian Stationary Lagrangian surfaces with bounded areas and Willmore energies in a complete K{\"a}hler surface. We also prove two strong compactness theorems on the space of Hamiltonian stationary Lagrangian tori in $\mathbb C^2$ and $\mathbb{CP}^2$ respectively.


title: ClearGrasp: 3D Shape Estimation of Transparent Objects for Manipulation

abstract: Transparent objects are a common part of everyday life, yet they possess unique visual properties that make them incredibly difficult for standard 3D sensors to produce accurate depth estimates for. In many cases, they often appear as noisy or distorted approximations of the surfaces that lie behind them. To address these challenges, we present ClearGrasp -- a deep learning approach for estimating accurate 3D geometry of transparent objects from a single RGB-D image for robotic manipulation. Given a single RGB-D image of transparent objects, ClearGrasp uses deep convolutional networks to infer surface normals, masks of transparent surfaces, and occlusion boundaries. It then uses these outputs to refine the initial depth estimates for all transparent surfaces in the scene. To train and test ClearGrasp, we construct a large-scale synthetic dataset of over 50,000 RGB-D images, as well as a real-world test benchmark with 286 RGB-D images of transparent objects and their ground truth geometries. The experiments demonstrate that ClearGrasp is substantially better than monocular depth estimation baselines and is capable of generalizing to real-world images and novel objects. We also demonstrate that ClearGrasp can be applied out-of-the-box to improve grasping algorithms' performance on transparent objects. Code, data, and benchmarks will be released. Supplementary materials available on the project website: https://sites.google.com/view/cleargrasp


title: Soft-Label Dataset Distillation and Text Dataset Distillation

abstract: Dataset distillation is a method for reducing dataset sizes by learning a small number of synthetic samples containing all the information of a large dataset. This has several benefits like speeding up model training, reducing energy consumption, and reducing required storage space. Currently, each synthetic sample is assigned a single `hard' label, and also, dataset distillation can currently only be used with image data. We propose to simultaneously distill both images and their labels, thus assigning each synthetic sample a `soft' label (a distribution of labels). Our algorithm increases accuracy by 2-4% over the original algorithm for several image classification tasks. Using `soft' labels also enables distilled datasets to consist of fewer samples than there are classes as each sample can encode information for multiple classes. For example, training a LeNet model with 10 distilled images (one per class) results in over 96% accuracy on MNIST, and almost 92% accuracy when trained on just 5 distilled images. We also extend the dataset distillation algorithm to distill sequential datasets including texts. We demonstrate that text distillation outperforms other methods across multiple datasets. For example, models attain almost their original accuracy on the IMDB sentiment analysis task using just 20 distilled sentences. Our code can be found at $\href{https://github.com/ilia10000/dataset-distillation}{\text{https://github.com/ilia10000/dataset-distillation}}$.


title: Constraint-Preconditioned Krylov Solvers for Regularized Saddle-Point Systems

abstract: We consider the iterative solution of regularized saddle-point systems. When the leading block is symmetric and positive semi-definite on an appropriate subspace, Dollar, Gould, Schilders, and Wathen (2006) describe how to apply the conjugate gradient (CG) method coupled with a constraint preconditioner, a choice that has proved to be effective in optimization applications. We investigate the design of constraint-preconditioned variants of other Krylov methods for regularized systems by focusing on the underlying basis-generation process. We build upon principles laid out by Gould, Orban, and Rees (2014) to provide general guidelines that allow us to specialize any Krylov method to regularized saddle-point systems. In particular, we obtain constraint-preconditioned variants of Lanczos and Arnoldi-based methods, including the Lanczos version of CG, MINRES, SYMMLQ, GMRES(m) and DQGMRES. We also provide MATLAB implementations in hopes that they are useful as a basis for the development of more sophisticated software. Finally, we illustrate the numerical behavior of constraint-preconditioned Krylov solvers using symmetric and nonsymmetric systems arising from constrained optimization.