title: Don't Patronize Me! An Annotated Dataset with Patronizing and Condescending Language towards Vulnerable Communities

abstract: In this paper, we introduce a new annotated dataset which is aimed at supporting the development of NLP models to identify and categorize language that is patronizing or condescending towards vulnerable communities (e.g. refugees, homeless people, poor families). While the prevalence of such language in the general media has long been shown to have harmful effects, it differs from other types of harmful language, in that it is generally used unconsciously and with good intentions. We furthermore believe that the often subtle nature of patronizing and condescending language (PCL) presents an interesting technical challenge for the NLP community. Our analysis of the proposed dataset shows that identifying PCL is hard for standard NLP models, with language models such as BERT achieving the best results.


title: Nonparametric Bayesian volatility estimation for gamma-driven stochastic differential equations

abstract: We study a nonparametric Bayesian approach to estimation of the volatility function of a stochastic differential equation driven by a gamma process. The volatility function is modelled a priori as piecewise constant, and we specify a gamma prior on its values. This leads to a straightforward procedure for posterior inference via an MCMC procedure. We give theoretical performance guarantees (contraction rates for the posterior) for the Bayesian estimate in terms of the regularity of the unknown volatility function. We illustrate the method on synthetic and real data examples.


title: $p$-Summable commutators in Bergman spaces on egg domains

abstract: The exact range of the positive real parameter $p$ is determined so that the multiplications by polynomial functions acting on Bergman spaces over complex ellipsoids are Schatten $p$-summable.


title: First-principles study of the electronic, magnetic, and crystal structure of perovskite molybdates

abstract: The molybdate oxides SrMoO$_3$, PbMoO$_3$, and LaMoO$_3$ are a class of metallic perovskites that exhibit interesting properties including high mobility, and unusual resistivity behavior. We use first-principles methods based on density functional theory to explore the electronic, crystal, and magnetic structure of these materials. In order to account for the electron correlations in the partially-filled Mo $4d$ shell, a local Hubbard $U$ interaction is included. The value of $U$ is estimated via the constrained random-phase approximation approach, and the dependence of the results on the choice of $U$ are explored. For all materials, GGA+$U$ predicts a metal with an orthorhombic, antiferromagnetic structure. For LaMoO$_3$, the $Pnma$ space group is the most stable, while for SrMoO$_3$ and PbMoO$_3$, the $Imma$ and $Pnma$ structures are close in energy. The $R_4^+$ octahedral rotations for SrMoO$_3$ and PbMoO$_3$ are found to be overestimated compared to the experimental low-temperature structure.


title: Examining the Feasibility of Off-the-Shelf Algorithms for Masking Directly Identifiable Information in Social Media Data

abstract: The identification and removal/replacement of protected information from social media data is an understudied problem, despite being desirable from an ethical and legal perspective. This paper identifies types of potentially directly identifiable information (inspired by protected health information in clinical texts) contained in tweets that may be readily removed using off-the-shelf algorithms, introduces an English dataset of tweets annotated for identifiable information, and compiles these off-the-shelf algorithms into a tool (Nightjar) to evaluate the feasibility of using Nightjar to remove directly identifiable information from the tweets. Nightjar as well as the annotated data can be retrieved from https://bitbucket.org/mdredze/nightjar.


title: A New Similarity Space Tailored for Supervised Deep Metric Learning

abstract: We propose a novel deep metric learning method. Differently from many works on this area, we defined a novel latent space obtained through an autoencoder. The new space, namely S-space, is divided into different regions that describe the positions where pairs of objects are similar/dissimilar. We locate makers to identify these regions. We estimate the similarities between objects through a kernel-based t-student distribution to measure the markers' distance and the new data representation. In our approach, we simultaneously estimate the markers' position in the S-space and represent the objects in the same space. Moreover, we propose a new regularization function to avoid similar markers to collapse altogether. We present evidences that our proposal can represent complex spaces, for instance, when groups of similar objects are located in disjoint regions. We compare our proposal to 9 different distance metric learning approaches (four of them are based on deep-learning) on 28 real-world heterogeneous datasets. According to the four quantitative metrics used, our method overcomes all the nine strategies from the literature.


title: Cryptanalysis of a code-based full-time signature

abstract: We present an attack against a code-based signature scheme based on the Lyubashevsky protocol that was recently proposed by Song, Huang, Mu, Wu and Wang (SHMWW). The private key in the SHMWW scheme contains columns coming in part from an identity matrix and in part from a random matrix. The existence of two types of columns leads to a strong bias in the distribution of set bits in produced signatures. Our attack exploits such a bias to recover the private key from a bunch of collected signatures. We provide a theoretical analysis of the attack along with experimental evaluations, and we show that as few as 10 signatures are enough to be collected for successfully recovering the private key. As for previous attempts of adapting Lyubashevsky's protocol to the case of code-based cryptography, the SHMWW scheme is thus proved unable to provide acceptable security. This confirms that devising secure code-based signature schemes with efficiency comparable to that of other post-quantum solutions (e.g., based on lattices) is still a challenging task.


title: On the Capacity of Intensity-Modulation Direct-Detection Gaussian Optical Wireless Communication Channels: A Tutorial

abstract: Optical wireless communication (OWC) using intensity-modulation and direct-detection (IM/DD) has a channel model which possesses unique features, due to the constraints imposed on the channel input. The aim of this tutorial is to overview results on the capacity of IM/DD channels with input-independent Gaussian noise as a model of OWC channels. It provides the reader with an entry point to the topic, and highlights some major contributions in this area. It begins with a discussion on channel models and how this IM/DD Gaussian channel model comes about, in addition to an explanation of input constraints. Then, it discusses the capacity of the single-input single-output channel, its computation, and capacity bounds and asymptotic capacity results. Then, it extends the discussion to the multiple-input multiple-output setup, and reviews capacity bounds for this channel model. Finally, it discusses multi-user channels modelled as a broadcast channel (downlink) or a multiple-access channel (uplink), with their associated capacity bounds.


title: The LOFAR Two Meter Sky Survey: Deep Fields, I -- Direction-dependent calibration and imaging

abstract: The Low Frequency Array (LOFAR) is an ideal instrument to conduct deep extragalactic surveys. It has a large field of view and is sensitive to large scale and compact emission. It is, however, very challenging to synthesize thermal noise limited maps at full resolution, mainly because of the complexity of the low-frequency sky and the direction dependent effects (phased array beams and ionosphere). In this first paper of a series we present a new calibration and imaging pipeline that aims at producing high fidelity, high dynamic range images with LOFAR High Band Antenna data, while being computationally efficient and robust against the absorption of unmodeled radio emission. We apply this calibration and imaging strategy to synthesize deep images of the Bootes and LH fields at 150 MHz, totaling $\sim80$ and $\sim100$ hours of integration respectively and reaching unprecedented noise levels at these low frequencies of $\lesssim30$ and $\lesssim23$ $\mu$Jy/beam in the inner $\sim3$ deg$^2$. This approach is also being used to reduce the LoTSS-wide data for the second data release.


title: Magnetic Field Line Random Walk and Solar Energetic Particle Path Lengths: Stochastic Theory and PSP/ISoIS Observation

abstract: Context:In 2020 May-June, six solar energetic ion events were observed by the Parker Solar Probe/ISoIS instrument suite at 0.35 AU from the Sun. From standard velocity-dispersion analysis, the apparent ion path length is 0.625 AU at the onset of each event. Aims:We develop a formalism for estimating the path length of random-walking magnetic field lines, to explain why the apparent ion pathlength at event onset greatly exceeds the radial distance from the Sun for these events. Methods:We developed analytical estimates of the average increase in pathlength of random-walking magnetic field lines, relative to the unperturbed mean field. Monte Carlo simulations of fieldline and particle trajectories in a model of solar wind turbulence are used to validate the formalism and study the path lengths of particle guiding-center and full-orbital trajectories. The formalism is implemented in a global solar wind model, and results are compared with ion pathlengths inferred from ISoIS observations. Results:Both a simple estimate and a rigorous theoretical formulation are obtained for fieldlines' pathlength increase as a function of pathlength along the large-scale field. From simulated fieldline and particle trajectories, we find that particle guiding centers can have pathlengths somewhat shorter than the average fieldline pathlength, while particle orbits can have substantially larger pathlengths due to their gyromotion with a nonzero effective pitch angle. Conclusions:The long apparent path length during these solar energetic ion events can be explained by 1) a magnetic field line path length increase due to the field line random walk, and 2) particle transport about the guiding center with a nonzero effective pitch angle. Our formalism for computing the magnetic field line path length, accounting for turbulent fluctuations, may be useful for application to solar particle transport in general.


title: Laserless quantum gates for electric dipoles in thermal motion

abstract: Internal states of polar molecules can be controlled by microwave-frequency electric dipole transitions. If the applied microwave electric field has a spatial gradient, these transitions also affect the motion of these dipolar particles. This capability can be used to engineer phonon-mediated quantum gates between e.g. trapped polar molecular ion qubits without laser illumination and without the need for cooling near the motional ground state. The result is a high-speed quantum processing toolbox for dipoles in thermal motion that combines the precision microwave control of solid-state qubits with the long coherence times of trapped ion qubits.


title: Predicting high rate granular transition and fragment statistics at the onset of granular flow for brittle ceramics

abstract: Brittle materials under impact loading exhibit a transition from a cracked solid to a granular medium. Appropriate representation of this transition to granular mechanics and the resulting initial fragment size and shape distribution in computational models is not well understood. The current work provides a numerical model to analyze competitive crack coalescence in the transition regime and provides insight into the onset of comminution and the initial conditions for subsequent granular flow. Crack statistics obtained from initial flaws using a wing crack growth based damage model have been used to discretely model elliptical cracks in three dimensions, with and without a minimal intersection constraint. These cracks are then allowed to coalesce with nearby cracks along favourable directions and the output fragment statistics have been predicted. The evolving fragmentation offers insight into the onset of comminution as well as the final transition to granular mechanics and the resulting initial fragment statistics. A simple phenomenological model has been proposed that suggests a transition criterion resembling the one obtained from the numerical model.


title: EffiScene: Efficient Per-Pixel Rigidity Inference for Unsupervised Joint Learning of Optical Flow, Depth, Camera Pose and Motion Segmentation

abstract: This paper addresses the challenging unsupervised scene flow estimation problem by jointly learning four low-level vision sub-tasks: optical flow $\textbf{F}$, stereo-depth $\textbf{D}$, camera pose $\textbf{P}$ and motion segmentation $\textbf{S}$. Our key insight is that the rigidity of the scene shares the same inherent geometrical structure with object movements and scene depth. Hence, rigidity from $\textbf{S}$ can be inferred by jointly coupling $\textbf{F}$, $\textbf{D}$ and $\textbf{P}$ to achieve more robust estimation. To this end, we propose a novel scene flow framework named EffiScene with efficient joint rigidity learning, going beyond the existing pipeline with independent auxiliary structures. In EffiScene, we first estimate optical flow and depth at the coarse level and then compute camera pose by Perspective-$n$-Points method. To jointly learn local rigidity, we design a novel Rigidity From Motion (RfM) layer with three principal components: \emph{}{(i)} correlation extraction; \emph{}{(ii)} boundary learning; and \emph{}{(iii)} outlier exclusion. Final outputs are fused based on the rigid map $M_R$ from RfM at finer levels. To efficiently train EffiScene, two new losses $\mathcal{L}_{bnd}$ and $\mathcal{L}_{unc}$ are designed to prevent trivial solutions and to regularize the flow boundary discontinuity. Extensive experiments on scene flow benchmark KITTI show that our method is effective and significantly improves the state-of-the-art approaches for all sub-tasks, i.e. optical flow ($5.19 \rightarrow 4.20$), depth estimation ($3.78 \rightarrow 3.46$), visual odometry ($0.012 \rightarrow 0.011$) and motion segmentation ($0.57 \rightarrow 0.62$).


title: 2D+3D Facial Expression Recognition via Discriminative Dynamic Range Enhancement and Multi-Scale Learning

abstract: In 2D+3D facial expression recognition (FER), existing methods generate multi-view geometry maps to enhance the depth feature representation. However, this may introduce false estimations due to local plane fitting from incomplete point clouds. In this paper, we propose a novel Map Generation technique from the viewpoint of information theory, to boost the slight 3D expression differences from strong personality variations. First, we examine the HDR depth data to extract the discriminative dynamic range $r_{dis}$, and maximize the entropy of $r_{dis}$ to a global optimum. Then, to prevent the large deformation caused by over-enhancement, we introduce a depth distortion constraint and reduce the complexity from $O(KN^2)$ to $O(KN\tau)$. Furthermore, the constrained optimization is modeled as a $K$-edges maximum weight path problem in a directed acyclic graph, and we solve it efficiently via dynamic programming. Finally, we also design an efficient Facial Attention structure to automatically locate subtle discriminative facial parts for multi-scale learning, and train it with a proposed loss function $\mathcal{L}_{FA}$ without any facial landmarks. Experimental results on different datasets show that the proposed method is effective and outperforms the state-of-the-art 2D+3D FER methods in both FER accuracy and the output entropy of the generated maps.


title: Widening the Dialogue Workflow Modeling Bottleneck in Ontology-Based Personal Assistants

abstract: We present a new approach to dialogue specification for Virtual Personal Assistants (VPAs) based on so-called dialogue workflow graphs, with several demonstrated advantages over current ontology-based methods. Our new dialogue specification language (DSL) enables customers to more easily participate in the VPA modeling process due to a user-friendly modeling framework. Resulting models are also significantly more compact. VPAs can be developed much more rapidly. The DSL is a new modeling layer on top of our ontology-based Dialogue Management (DM) framework OntoVPA. We explain the rationale and benefits behind the new language and support our claims with concrete reduced Level-of-Effort (LOE) numbers from two recent OntoVPA projects.


title: Ab initio electronic density in solids by many-body plane-wave auxiliary-field quantum Monte Carlo calculations

abstract: We present accurate many-body results of the electronic densities in several solid materials, including Si, NaCl, and Cu. These results are obtained using the ab initio auxiliary-field quantum Monte Carlo (AFQMC) method working in a plane-wave basis with norm-conserving, multiple-projector pseudopotentials. AFQMC has been shown to be an excellent many-body total energy method. Computation of observables and correlation functions other than the ground-state energy requires back-propagation, whose adaption and implementation in the plane-wave basis AFQMC framework are discussed in the present paper. This development allows us to compute correlation functions, electronic densities and interatomic forces, paving the way for geometry optimizations and calculations of thermodynamic properties in solids. Finite supercell size effects are considerably more subtle in the many-body framework than in independent-electron calculations. We analyze the convergence of the electronic density, and obtain best estimates for the thermodynamic limit. The densities from several typical density functionals are benchmarked against our near-exact results. The electronic densities we have obtained can also be used to help construct improved density functionals.


title: Statistical Dependencies Beyond Linear Correlations in Light Scattered by Disordered Media

abstract: Imaging through scattering and random media is an outstanding problem that to date has been tackled by either measuring the medium transmission matrix or exploiting linear correlations in the transmitted speckle patterns. However, transmission matrix techniques require interferometric stability and linear correlations, such as the memory effect, can be exploited only in thin scattering media. Here we show the existence of a statistical dependency in strongly scattered optical fields in a case where first-order correlations are not expected. We also show that this statistical dependence and the related information transport is directly linked to artificial neural network imaging in strongly scattering, dynamic media. These non-trivial dependencies provide a key to imaging through dynamic and thick scattering media with applications for deep-tissue imaging or imaging through smoke or fog


title: Optimal Transport-based Coverage Control for Swarm Robot Systems: Generalization of the Voronoi Tessellation-based Method

abstract: Swarm robot systems, which consist of many cooperating mobile robots, have attracted attention for their environmental adaptability and fault tolerance advantages. One of the most important tasks for such systems is coverage control, in which robots autonomously deploy to approximate a given spatial distribution. In this study, we formulate a coverage control paradigm using the concept of optimal transport and propose a novel control technique, which we have termed the optimal transport-based coverage control (OTCC) method. The proposed OTCC, derived via the gradient flow of the cost function in the Kantorovich dual problem, is shown to covers a widely used existing control method as a special case. We also perform a Lyapunov stability analysis of the controlled system, and provide numerical calculations to show that the OTCC reproduces target distributions with better performance than the existing control method.


title: Noise and Pulse Dynamics in Backward Stimulated Brillouin Scattering

abstract: We theoretically and numerically study the effects of thermal noise on pulses in backwards Stimulated Brillouin Scattering (SBS). Using a combination of stochastic calculus and numerical methods, we derive a theoretical model that can be used to quantitatively predict noise measurements. We study how the optical pulse configuration, including the input powers of the pump and Stokes fields, pulse durations and interaction time, affects the noise in the output Stokes field. We investigate the effects on the noise of the optical loss and waveguide length, and we find that the signal-to-noise ratio can be significantly improved, or reduced, for specific combinations of waveguide properties and pulse parameters.


title: Vector valued nonuniform multiresolution analysis associated with linear canonical transform

abstract: A multiresolution analysis associated with linear canonical transform was defined by Shah and Waseem for which the translation set is a discrete set which is not a group. In this paper, we continue the study based on this nonstandard setting and introduce vector-valued nonuniform multiresolution analysis associated with linear canonical transform (LCT-VNUMRA). We establish a necessary and sufficient condition for the existence of associated wavelets and derive an algorithm for the construction of vector-valued nonuniform multiresolution analysis on local fields starting from a vector refinement mask with appropriate conditions.


title: Better Automatic Program Repair by Using Bug Reports and Tests Together

abstract: Automated program repair is already deployed in industry, but concerns remain about repair quality. Recent research has shown that one of the main reasons repair tools produce incorrect (but seemingly correct) patches is imperfect fault localization (FL). This paper demonstrates that combining information from natural-language bug reports and test executions when localizing faults can have a significant positive impact on repair quality. For example, existing repair tools with such FL are able to correctly repair 7 defects in the Defects4J benchmark that no prior tools have repaired correctly. We develop, Blues, the first information-retrieval-based, statement-level FL technique that requires no training data. We further develop RAFL, the first unsupervised method for combining multiple FL techniques, which outperforms a supervised method. Using RAFL, we create SBIR by combining Blues with a spectrum-based (SBFL) technique. Evaluated on 815 real-world defects, SBIR consistently ranks buggy statements higher than its underlying techniques. We then modify three state-of-the-art repair tools, Arja, SequenceR, and SimFix, to use SBIR, SBFL, and Blues as their internal FL. We evaluate the quality of the produced patches on 689 real-world defects. Arja and SequenceR significantly benefit from SBIR: Arja using SBIR correctly repairs 28 defects, but only 21 using SBFL, and only 15 using Blues; SequenceR using SBIR correctly repairs 12 defects, but only 10 using SBFL, and only 4 using Blues. SimFix, (which has internal mechanisms to overcome poor FL), correctly repairs 30 defects using SBIR and SBFL, but only 13 using Blues. Our work is the first investigation of simultaneously using multiple software artifacts for automated program repair, and our promising findings suggest future research in this directions is likely to be fruitful.


title: Robust Deep Learning with Active Noise Cancellation for Spatial Computing

abstract: This paper proposes CANC, a Co-teaching Active Noise Cancellation method, applied in spatial computing to address deep learning trained with extreme noisy labels. Deep learning algorithms have been successful in spatial computing for land or building footprint recognition. However a lot of noise exists in ground truth labels due to how labels are collected in spatial computing and satellite imagery. Existing methods to deal with extreme label noise conduct clean sample selection and do not utilize the remaining samples. Such techniques can be wasteful due to the cost of data retrieval. Our proposed CANC algorithm not only conserves high-cost training samples but also provides active label correction to better improve robust deep learning with extreme noisy labels. We demonstrate the effectiveness of CANC for building footprint recognition for spatial computing.


title: Community Detection through Vector-label Propagation Algorithms

abstract: Community detection is a fundamental and important problem in network science, as community structures often reveal both topological and functional relationships between different components of the complex system. In this paper, we first propose a gradient descent framework of modularity optimization called vector-label propagation algorithm (VLPA), where a node is associated with a vector of continuous community labels instead of one label. Retaining weak structural information in vector-label, VLPA outperforms some well-known community detection methods, and particularly improves the performance in networks with weak community structures. Further, we incorporate stochastic gradient strategies into VLPA to avoid stuck in the local optima, leading to the stochastic vector-label propagation algorithm (sVLPA). We show that sVLPA performs better than Louvain Method, a widely used community detection algorithm, on both artificial benchmarks and real-world networks. Our theoretical scheme based on vector-label propagation can be directly applied to high-dimensional networks where each node has multiple features, and can also be used for optimizing other partition measures such as modularity with resolution parameters.


title: Option Pricing Incorporating Factor Dynamics in Complete Markets

abstract: Using the Donsker-Prokhorov invariance principle we extend the Kim-Stoyanov-Rachev-Fabozzi option pricing model to allow for variably-spaced trading instances, an important consideration for short-sellers of options. Applying the Cherny-Shiryaev-Yor invariance principles, we formulate a new binomial path-dependent pricing model for discrete- and continuous-time complete markets where the stock price dynamics depends on the log-return dynamics of a market influencing factor. In the discrete case, we extend the results of this new approach to a financial market with informed traders employing a statistical arbitrage strategy involving trading of forward contracts. Our findings are illustrated with numerical examples employing US financial market data. Our work provides further support for the conclusion that any option pricing model must preserve valuable information on the instantaneous mean log-return, the probability of the stock's upturn movement (per trading interval), and other market microstructure features.


title: A Stochastic Approach to Unitary Coupled Cluster

abstract: Unitary coupled cluster (UCC), originally developed as a variational alternative to the popular traditional coupled cluster method, has seen a resurgence as a functional form for use on quantum computers. However, the number of excitors present in the ansatz often presents a barrier to implementation on quantum computers. Given the natural sparsity of wavefunctions obtained from Quantum Monte Carlo methods, we consider here a stochastic solution to the UCC problem. Using the Coupled Cluster Monte Carlo framework, we develop cluster selection schemes that capture the structure of the UCC wavefunction, as well as its Trotterized approximation, and use these to solve the corresponding projected equations. Due to the fast convergence of the equations with order in the cluster expansion, this approach scales polynomially with the size of the system. Unlike traditional UCC implementations, our approach naturally produces a non-variational estimator for the energy in the form of the projected energy. For UCCSD in small systems, we find this agrees well with the expectation value of the energy and, in the case of two electrons, with full configuration interaction results. For the larger N$_2$ system, the two estimators diverge, with the projected energy approaching the coupled cluster result, while the expectation value is close to results from traditional UCCSD.


title: Distilling a Hierarchical Policy for Planning and Control via Representation and Reinforcement Learning

abstract: We present a hierarchical planning and control framework that enables an agent to perform various tasks and adapt to a new task flexibly. Rather than learning an individual policy for each particular task, the proposed framework, DISH, distills a hierarchical policy from a set of tasks by representation and reinforcement learning. The framework is based on the idea of latent variable models that represent high-dimensional observations using low-dimensional latent variables. The resulting policy consists of two levels of hierarchy: (i) a planning module that reasons a sequence of latent intentions that would lead to an optimistic future and (ii) a feedback control policy, shared across the tasks, that executes the inferred intention. Because the planning is performed in low-dimensional latent space, the learned policy can immediately be used to solve or adapt to new tasks without additional training. We demonstrate the proposed framework can learn compact representations (3- and 1-dimensional latent states and commands for a humanoid with 197- and 36-dimensional state features and actions) while solving a small number of imitation tasks, and the resulting policy is directly applicable to other types of tasks, i.e., navigation in cluttered environments. Video: https://youtu.be/HQsQysUWOhg


title: Refining Automatic Speech Recognition System for older adults

abstract: Building a high quality automatic speech recognition (ASR) system with limited training data has been a challenging task particularly for a narrow target population. Open-sourced ASR systems, trained on sufficient data from adults, are susceptible on seniors' speech due to acoustic mismatch between adults and seniors. With 12 hours of training data, we attempt to develop an ASR system for socially isolated seniors (80+ years old) with possible cognitive impairments. We experimentally identify that ASR for the adult population performs poorly on our target population and transfer learning (TL) can boost the system's performance. Standing on the fundamental idea of TL, tuning model parameters, we further improve the system by leveraging an attention mechanism to utilize the model's intermediate information. Our approach achieves 1.58% absolute improvements over the TL model.


title: Complexity, Exactness, and Rationality in Polynomial Optimization

abstract: We focus on rational solutions or nearly-feasible rational solutions that serve as certificates of feasibility for polynomial optimization problems. We show that, under some separability conditions, certain cubic polynomially constrained sets admit rational solutions. However, we show in other cases that it is NP Hard to detect if rational solutions exist or if they exist of any reasonable size. We extend this idea to various settings including near feasible, but super optimal solutions and detecting rational rays on which a cubic function is unbounded. Lastly, we show that in fixed dimension, the feasibility problem over a set defined by polynomial inequalities is in NP by providing a simple certificate to verify feasibility. We conclude with several related examples of irrationality and encoding size issues in QCQPs and SOCPs.


title: Uniform effective estimates for $\vert L(1,\chi)\vert$

abstract: Let $L(s,\chi)$ be the Dirichlet $L$-function associated to a non-principal primitive Dirichlet character $\chi$ defined modulo $q$, where $q\ge 3$. We prove, under the assumption of the Generalised Riemann Hypothesis, the validity of estimates given by Lamzouri, Li, and Soundararajan on $\vert L(1,\chi) \vert$. As a corollary, we have that similar estimates hold for the class number of the imaginary quadratic field $\mathbb{Q}(\sqrt{-q})$, $q\ge 5$.


title: First Search for Low-Frequency CH with a Square Kilometre Array Precursor Telescope

abstract: The diatomic free radical methylidyne (CH) is an important tracer of the interstellar medium and the study of it was critical to our earliest understanding of star formation. Although it is detectable across the electromagnetic spectrum, observations at radio frequencies allow for a study of the kinematics of the diffuse and dense gas in regions of new star formation. There are only two published (single-dish) detections of the low-frequency hyperfine transitions between 700 and 725 MHz, despite the precise frequencies being known. These low-frequency transitions are of particular interest as they are shown in laboratory experiments to be more sensitive to magnetic fields than their high-frequency counterparts (with more pronounced Zeeman splitting). In this work we take advantage of the radio quiet environment and increased resolution of the Australian Square Kilometre Array Pathfinder (ASKAP) over previous searches to make a pilot interferometric search for CH at 724.7883 MHz (the strongest of the hyperfine transitions) in RCW 38. We found the band is clean of radio frequency interference, but we did not detect the signal from this transition to a five sigma sensitivity limit of 0.09 Jy, which corresponds to a total column density upper limit of 1.9x10^18cm^-2 for emission and 1.3x10^14cm^-2 for absorption with an optical depth limit of 0.95. Achieved within 5 hrs of integration, this column density sensitivity should have been adequate to detect the emission or absorption in RCW 38, if it had similar properties to the only previous reported detections in W51.


title: Defying the Circadian Rhythm: Clustering Participant Telemetry in the UK Biobank Data

abstract: The UK Biobank dataset follows over 500,000 volunteers and contains a diverse set of information related to societal outcomes. Among this vast collection, a large quantity of telemetry collected from wrist-worn accelerometers provides a snapshot of participant activity. Using this data, a population of shift workers, subjected to disrupted circadian rhythms, is analysed using a mixture model-based approach to yield protective effects from physical activity on survival outcomes. In this paper, we develop a scalable, standardized, and unique methodology that efficiently clusters a vast quantity of participant telemetry. By building upon the work of Doherty et al. (2017), we introduce a standardized, low-dimensional feature for clustering purposes. Participants are clustered using a matrix variate mixture model-based approach. Once clustered, survival analysis is performed to demonstrate distinct lifetime outcomes for individuals within each cluster. In summary, we process, cluster, and analyse a subset of UK Biobank participants to show the protective effects from physical activity on circadian disrupted individuals.


title: GEMS embeddings of Schwarzschild and RN black holes in Painlev\'e-Gullstrand spacetimes

abstract: Making use of the higher dimensional global embedding Minkowski spacetime (GEMS), we embed (3+1)-dimensional Schwarzschild and Reissner-Nordstr\"om (RN) black holes written by the Painlev\'e-Gullstrand (PG) spacetimes, which have off-diagonal components in metrics, into (5+1)- and (5+2)-dimensional flat ones, respectively. As a result, we have shown the equivalence of the GEMS embeddings of the spacetimes with the diagonal and off-diagonal terms in metrics. Moreover, with the aid of their geodesic equations satisfying various boundary conditions in the flat embedded spacetimes, we directly obtain freely falling temperatures. We also show that freely falling temperatures in the PG spacetimes are well-defined beyond the event horizons, while they are equivalent to the Hawking temperatures, which are obtained in the original curved ones in the ranges between the horizon and the infinity. These will be helpful to study GEMS embeddings of more realistic Kerr, or rotating BTZ black holes.


title: Brillouin-Mandelstam Light Scattering Spectroscopy: Applications in Phononics and Spintronics

abstract: Recent years witnessed much broader use of Brillouin inelastic light scattering spectroscopy for the investigation of phonons and magnons in novel materials, nanostructures, and devices. Driven by developments in instrumentation and the strong need for accurate knowledge of energies of elemental excitations, the Brillouin - Mandelstam spectroscopy is rapidly becoming an essential technique, complementary to the Raman inelastic light scattering spectroscopy. We provide an overview of recent progress in the Brillouin light scattering technique, focusing on the use of this photonic method for the investigation of confined acoustic phonons, phononic metamaterials, magnon propagation and scattering. The Review emphasizes emerging applications of the Brillouin - Mandelstam spectroscopy for phonon engineered structures and spintronic devices and concludes with a perspective for future directions.


title: AXES: Approximation Manager for Emerging Memory Architectures

abstract: Memory approximation techniques are commonly limited in scope, targeting individual levels of the memory hierarchy. Existing approximation techniques for a full memory hierarchy determine optimal configurations at design-time provided a goal and application. Such policies are rigid: they cannot adapt to unknown workloads and must be redesigned for different memory configurations and technologies. We propose AXES: the first self-optimizing runtime manager for coordinating configurable approximation knobs across all levels of the memory hierarchy. AXES continuously updates and optimizes its approximation management policy throughout runtime for diverse workloads. AXES optimizes the approximate memory configuration to minimize power consumption without compromising the quality threshold specified by application developers. AXES can (1) learn a policy at runtime to manage variable application quality of service (QoS) constraints, (2) automatically optimize for a target metric within those constraints, and (3) coordinate runtime decisions for interdependent knobs and subsystems. We demonstrate AXES' ability to efficiently provide functions 1-3 on a RISC-V Linux platform with approximate memory segments in the on-chip cache and main memory. We demonstrate AXES' ability to save up to 37% energy in the memory subsystem without any design-time overhead. We show AXES' ability to reduce QoS violations by 75% with $<5\%$ additional energy.


title: Graded Lie algebras of maximal class of type $p$

abstract: The algebras of the title are infinite-dimensional graded Lie algebras $L= \bigoplus_{i=1}^{\infty}L_i$, over a field of positive characteristic $p$, which are generated by an element of degree $1$ and an element of degree $p$, and satisfy $[L_i,L_1]=L_{i+1}$ for $i\ge p$. %of maximal class in the sense that $L/L^i$ has dimension $i$ for all $i>1$. In case $p=2$ such algebras were classified by Caranti and Vaughan-Lee in 2003. We announce an extension of that classification to arbitrary prime characteristic, and prove several major steps in its proof.


title: On a Reaction-Diffusion System Modeling Infectious Diseases Without Life-time Immunity

abstract: In this paper we study a mathematical model for an infectious disease such as Cholera without life-time immunity. Due to the different mobility for susceptible, infected human and recovered human hosts, the diffusion coefficients are assumed to be different. The resulting system is governed by a strongly coupled reaction-diffusion system with different diffusion coefficients. Global existence and uniqueness are established under certain assumptions on known data. Moreover, global asymptotic behavior of the solution is obtained when some parameters satisfy certain conditions. These results extend the existing results in the literature. The main tool used in this paper comes from the delicate theory of elliptic and parabolic equations. Moreover, the energy method and Sobolev embedding are used in deriving {\em apriori} estimates. The analysis developed in this paper can be employed to study other epidemic models in biological and ecological systems.


title: Phenotyping Clusters of Patient Trajectories suffering from Chronic Complex Disease

abstract: Recent years have seen an increased focus into the tasks of predicting hospital inpatient risk of deterioration and trajectory evolution due to the availability of electronic patient data. A common approach to these problems involves clustering patients time-series information such as vital sign observations) to determine dissimilar subgroups of the patient population. Most clustering methods assume time-invariance of vital-signs and are unable to provide interpretability in clusters that is clinically relevant, for instance, event or outcome information. In this work, we evaluate three different clustering models on a large hospital dataset of vital-sign observations from patients suffering from Chronic Obstructive Pulmonary Disease. We further propose novel modifications to deal with unevenly sampled time-series data and unbalanced class distribution to improve phenotype separation. Lastly, we discuss further avenues of investigation for models to learn patient subgroups with distinct behaviour and phenotype.


title: Explicit formulas for Eigenvalues of Capelli operators for the Lie superalgebra $\frak{osp}(1|2n)$

abstract: We define a natural basis for the algebra of $\frak{gosp}(1|2n)$-invariant differential operators on the affine superspace $\mathbb{C}^{1|2n}$. We prove that these operators lie in the image of the centre of the enveloping algebra of $\frak{gosp}(1|2n)$. Using this result, we compute explicit formulas for the eigenvalues of these operators on irreducible summands of $\mathcal{P}(\mathbb{C}^{1|2n})$. This settles the Capelli eigenvalue problem for orthosymplectic Lie superalgebras in the cases that were not addressed in Sahi-Salmasian-Serganova. Our main technique relies on an explicit calculation of a certain determinant with polynomial entries.


title: Analytic Geometry and Hodge-Frobenius Structure

abstract: In this paper, we study Frobenius structures in higher dimensional $p$-adic analytic geometry and the corresponding $p$-adic functional analysis. This will build up foundations for further study on some generalized cohomology of Frobenius modules and the corresponding generalized Iwasawa theory and generalized noncommutative Tamagawa number conjectures in the spirit of Burns-Flach-Fukaya-Kato and Nakamura (as well as certainly the original noncommutative Tamagawa number conjectures as observed by Pal-Z\'abr\'adi). We will work in the program proposed by Carter-Kedlaya-Z\'abr\'adi and after Pal-Z\'abr\'adi, and we will follow closely the approach from Kedlaya-Pottharst-Xiao to investigate the corresponding deformation of the generalized $p$-adic Hodge structures.


title: Characterization of the Chemical and Electrical Properties of Defects at the Niobium-Silicon Interface

abstract: The nature and concentration of defects near niobium-silicon interfaces prepared with different silicon surface treatments were characterized using current-voltage (I-V), deep level transient spectroscopy (DLTS), and secondary ion mass spectroscopy (SIMS). All samples have H, C, O, F, and Cl chemical contamination in the Si within 50 nm of the interface and electrically active defects with activation energies of 0.147, 0.247, 0.339, and 0.556 eV above the valence band maximum (E$_{vbm}$). In all cases, the deep level defect concentration is dominated by the hole trap at E$_{vbm}$ + 0.556eV, which we assign to a Nb point defect in Si, presumably Nb$_\textrm{Si}$. This defect is present with concentrations ranging from $7\times10^{13}$ to $5\times10^{14}$ cm$^{-3}$ and depends on the final surface clean process. The optimum surface treatment used in this study is an HF etch followed by an in-situ 100 eV Ar-gas ion milling process. Higher energy ion milling is found to increase the electrically active Nb defect concentration in the Si, and increase the concentration of defects. The HF etch alone removes O from the interface, but results in significant H and F contamination, electrically-active point defect concentrations, and levels of Shockley-Reed-Hall recombination (i.e. Nb/Si Schottky diodes with an ideality factor, n, of $\approx$ 1.6). The RCA clean increases the depth and concentration of H, F, C, and Nb contamination.


title: Recursive Importance Sketching for Rank Constrained Least Squares: Algorithms and High-order Convergence

abstract: In this paper, we propose {\it \underline{R}ecursive} {\it \underline{I}mportance} {\it \underline{S}ketching} algorithm for {\it \underline{R}ank} constrained least squares {\it \underline{O}ptimization} (RISRO). The key step of RISRO is recursive importance sketching, a new sketching framework based on deterministically designed recursive projections, which significantly differs from the randomized sketching in the literature \citep{mahoney2011randomized,woodruff2014sketching}. Several existing algorithms in the literature can be reinterpreted under this new sketching framework and RISRO offers clear advantages over them. RISRO is easy to implement and computationally efficient, where the core procedure in each iteration is to solve a dimension-reduced least squares problem. We establish the local quadratic-linear and quadratic rate of convergence for RISRO under some mild conditions. We also discover a deep connection of RISRO to the Riemannian Gauss-Newton algorithm on fixed rank matrices. The effectiveness of RISRO is demonstrated in two applications in machine learning and statistics: low-rank matrix trace regression and phase retrieval. Simulation studies demonstrate the superior numerical performance of RISRO.

categories: math.OC cs.LG cs.NA math.NA stat.CO stat.ML


title: Knowledge-Augmented Dexterous Grasping with Incomplete Sensing

abstract: Humans can determine a proper strategy to grasp an object according to the measured physical attributes or the prior knowledge of the object. This paper proposes an approach to determining the strategy of dexterous grasping by using an anthropomorphic robotic hand simply based on a label or a description of an object. Object attributes are parsed from natural-language descriptions and augmented with an object knowledge base that is scraped from retailer websites. A novel metric named joint probability distance is defined to measure distance between object attributes. The probability distribution of grasp types for the given object is learned using a deep neural network which takes in object features as input. The action of the multi-fingered hand with redundant degrees of freedom (DoF) is controlled by a linear inverse-kinematics model of grasp topology and scales. The grasping strategy generated by the proposed approach is evaluated both by simulation and execution on a Sawyer robot with an AR10 robotic hand.


title: CG-Net: Conditional GIS-aware Network for Individual Building Segmentation in VHR SAR Images

abstract: Object retrieval and reconstruction from very high resolution (VHR) synthetic aperture radar (SAR) images are of great importance for urban SAR applications, yet highly challenging owing to the complexity of SAR data. This paper addresses the issue of individual building segmentation from a single VHR SAR image in large-scale urban areas. To achieve this, we introduce building footprints from GIS data as complementary information and propose a novel conditional GIS-aware network (CG-Net). The proposed model learns multi-level visual features and employs building footprints to normalize the features for predicting building masks in the SAR image. We validate our method using a high resolution spotlight TerraSAR-X image collected over Berlin. Experimental results show that the proposed CG-Net effectively brings improvements with variant backbones. We further compare two representations of building footprints, namely complete building footprints and sensor-visible footprint segments, for our task, and conclude that the use of the former leads to better segmentation results. Moreover, we investigate the impact of inaccurate GIS data on our CG-Net, and this study shows that CG-Net is robust against positioning errors in GIS data. In addition, we propose an approach of ground truth generation of buildings from an accurate digital elevation model (DEM), which can be used to generate large-scale SAR image datasets. The segmentation results can be applied to reconstruct 3D building models at level-of-detail (LoD) 1, which is demonstrated in our experiments.


title: Vis-CRF, A Classical Receptive Field Model for VISION

abstract: Over the last decade, a variety of new neurophysiological experiments have led to new insights as to how, when and where retinal processing takes place, and the nature of the retinal representation encoding sent to the cortex for further processing. Based on these neurobiological discoveries, in our previous work, we provided computer simulation evidence to suggest that Geometrical illusions are explained in part, by the interaction of multiscale visual processing performed in the retina. The output of our retinal stage model, named Vis-CRF, is presented here for a sample of natural image and for several types of Tilt Illusion, in which the final tilt percept arises from multiple scale processing of Difference of Gaussians (DoG) and the perceptual interaction of foreground and background elements (Nematzadeh and Powers, 2019; Nematzadeh, 2018; Nematzadeh, Powers and Lewis, 2017; Nematzadeh, Lewis and Powers, 2015).


title: On Integer Balancing of Digraphs

abstract: A weighted digraph is balanced if the sums of the weights of the incoming and of the outgoing edges are equal at each vertex. We show that if these sums are integers, then the edge weights can be integers as well.


title: The K2 M67 Study: Precise Mass for a Turnoff Star in the Old Open Cluster M67

abstract: We present a study of the bright detached eclipsing main sequence binary WOCS 11028 (Sanders 617) in the open cluster M67. Although the binary has only one eclipse per orbital cycle, we show that the masses of the stars can be derived very precisely thanks to a strong constraint on the orbital inclination: $M_A = 1.222\pm0.006 M_\odot$ and $M_B = 0.909\pm0.004 M_\odot$. We use a spectral energy distribution fitting method to derive characteristics of the component stars in lieu of the precise radii that would normally be derived from a doubly-eclipsing binary. The deconvolution of the SEDs reveals that the brighter component of the binary is at the faint turnoff point for the cluster -- a distinct evolutionary point that occurs after the convective core has been established and while the star is in the middle of its movement toward lower surface temperature, before the so-called hook at the end of main sequence. The measurements are in distinct disagreement with evolution models at solar metallicity: higher metal abundances are needed to reproduce the characteristics of WOCS 11028 A. We discuss the changes to model physics that are likely to be needed to address the discrepancies. The clearest conclusions are that diffusion is probably necessary to reconcile spectroscopic abundances of M67 stars with the need for higher metallicity models, and that reduced strength convective overshooting is occurring for stars at the turnoff. At super-solar bulk metallicity, various indicators agree on a cluster age between about 3.5 and 4.0 Gyr.


title: Uniform Bipartition in the Population Protocol Model with Arbitrary Communication Graphs

abstract: In this paper, we focus on the uniform bipartition problem in the population protocol model. This problem aims to divide a population into two groups of equal size. In particular, we consider the problem in the context of \emph{arbitrary} communication graphs. As a result, we clarify the solvability of the uniform bipartition problem with arbitrary communication graphs when agents in the population have designated initial states, under various assumptions such as the existence of a base station, symmetry of the protocol, and fairness of the execution. When the problem is solvable, we present protocols for uniform bipartition. When global fairness is assumed, the space complexity of our solutions is tight.


title: Extreme Value Preserving Networks

abstract: Recent evidence shows that convolutional neural networks (CNNs) are biased towards textures so that CNNs are non-robust to adversarial perturbations over textures, while traditional robust visual features like SIFT (scale-invariant feature transforms) are designed to be robust across a substantial range of affine distortion, addition of noise, etc with the mimic of human perception nature. This paper aims to leverage good properties of SIFT to renovate CNN architectures towards better accuracy and robustness. We borrow the scale-space extreme value idea from SIFT, and propose extreme value preserving networks (EVPNets). Experiments demonstrate that EVPNets can achieve similar or better accuracy than conventional CNNs, while achieving much better robustness on a set of adversarial attacks (FGSM,PGD,etc) even without adversarial training.


title: Acceleration radiation of an atom freely falling into a Kerr black hole and near-horizon conformal quantum mechanics

abstract: An atom falling freely into a Kerr black hole in a Boulware-like vacuum is shown to emit radiation with a Planck spectrum at the Hawking temperature. For a cloud of falling atoms with random initial times, the radiation is thermal. The existence of this radiation is due to the acceleration of the vacuum field modes with respect to the falling atom. Its properties can be traced to the dominant role of conformal quantum mechanics (CQM) in the neighborhood of the event horizon. We display this effect for a scalar field, though the acceleration radiation has a universal conformal behavior that is exhibited by all fields in the background of generic black holes.


title: Fredholm property and essential spectrum of $3-D$ Dirac operators with regular and singular potentials

abstract: We consider the $3-D$ Dirac operator $\mathfrak{D}_{\boldsymbol{A},\Phi ,Q_{\sin }}$ with variable regular magnetic and electrostatic potentials $ \boldsymbol{A}$,$\Phi $ and with singular potentials $Q_{\sin }$ with support on a smooth unbounded surface $\Sigma \subset \mathbb{R}^{3}$ which divides $\mathbb{R}^{3}$ on two open domains $\Omega_{\pm }$. We associate with the formal Dirac operator $\mathfrak{D}_{\boldsymbol{A},\Phi ,Q_{\sin }} $ an unbounded operator $\mathcal{D}_{\boldsymbol{A},\Phi ,Q_{\sin }}$ in $ L^{2}(\mathbb{R}^{3},\mathbb{C}^{4})$ generated by the regular part of $ \mathfrak{D}_{\boldsymbol{A},\Phi ,Q_{\sin }}$ with domain in $H^{1}(\Omega_{+},\mathbb{C}^{4})\oplus H^{1}(\Omega_{-},\mathbb{C}^{4})$ consisting of functions satisfying transmission conditions on $\Sigma .$ We consider the self-adjointness of operator $\mathcal{D}_{\boldsymbol{A},\Phi ,Q_{\sin }}$ for unbounded $C^{2}-$uniformly regular surfaces $\Sigma ,$ and the essential spectrum of $\mathcal{D}_{\boldsymbol{A},\Phi ,Q_{\sin }}$ if $ \Sigma $ is a $C^{2}$-surfaces with conic exits to infinity. As application we consider the electrostatic and Lorentz scalar $\delta_{\Sigma }-$shell interactions on unbounded surfaces $\Sigma .$


title: Change the coefficients of conditional entropies in extensivity

abstract: The Boltzmann--Gibbs entropy is a functional on the space of probability measures. When a state space is countable, one characterization of the Boltzmann--Gibbs entropy is given by the Shannon--Khinchin axioms, which consist of continuity, maximality, expandability and extensivity. Among these four properties, the extensivity is generalized in various ways. The extensivity of a functional is interpreted as the property that, for any random variables $(X,Y)$ taking finitely many values in $\mathbb{N}$, the difference between the value of the functional at the joint law of $(X,Y)$ and that at the law of $X$ coincides with the linear combinations of the values at the conditional laws of $Y$ given $X=n$ with coefficients given by the probabilities of each event $X=n$. A generalization of the extensivity obtained by replacing the coefficients with a power of the probabilities of the events $X=n$ provides a characterization of the Tsallis entropy. In this paper, we first prove the impossibility to replace the coefficients with a non-power function of the probabilities of the events $X=n$. Then we estimate the difference between the value at the joint law of $(X,Y)$ and that at the law of $X$ for a general functional.


title: The Crossing Number of Circulant Graph C(3k+1;{1,k}) on the Projective Plane

abstract: In this paper, we prove that the crossing number of circulant graph $C(3k+1;\{1,k\})$ on the projective plane is $k$ for $k \geq 3$.


title: The asymptotic spectrum of flipped multilevel Toeplitz matrices and of certain preconditionings

abstract: In this work, we perform a spectral analysis of flipped multilevel Toeplitz sequences, i.e., we study the asymptotic spectral behaviour of $\{Y_{\boldsymbol{n}}T_{\boldsymbol{n}}(f)\}_{\boldsymbol{n}}$, where $T_{\boldsymbol{n}}(f)$ is a real, square multilevel Toeplitz matrix generated by a function $f\in L^1([-\pi,\pi]^d)$ and $Y_{\boldsymbol{n}}$ is the exchange matrix, which has $1$s on the main anti-diagonal. In line with what we have shown for unilevel flipped Toeplitz matrix sequences, the asymptotic spectrum is determined by a $2\times 2$ matrix-valued function whose eigenvalues are $\pm |f|$. Furthermore, we characterize the eigenvalue distribution of certain preconditioned flipped multilevel Toeplitz sequences with an analysis that covers both multilevel Toeplitz and circulant preconditioners. Finally, all our findings are illustrated by several numerical experiments.


title: GPURepair: Automated Repair of GPU Kernels

abstract: This paper presents a tool for repairing errors in GPU kernels written in CUDA or OpenCL due to data races and barrier divergence. Our novel extension to prior work can also remove barriers that are deemed unnecessary for correctness. We implement these ideas in our tool called GPURepair, which uses GPUVerify as the verification oracle for GPU kernels. We also extend GPUVerify to support CUDA Cooperative Groups, allowing GPURepair to perform inter-block synchronization for CUDA kernels. To the best of our knowledge, GPURepair is the only tool that can propose a fix for intra-block data races and barrier divergence errors for both CUDA and OpenCL kernels and the only tool that fixes inter-block data races for CUDA kernels. We perform extensive experiments on about 750 kernels and provide a comparison with prior work. We demonstrate the superiority of GPURepair through its capability to fix more kernels and its unique ability to remove redundant barriers and handle inter-block data races.


title: Symmetric functions and Springer representations

abstract: The characters of the (total) Springer representations are identified with the Green functions by Kazhdan [Israel J. Math. {\bf 28} (1977)], and the latter are identified with Hall-Littlewood's $Q$-functions by Green [Trans. Amer. Math. Soc. (1955)]. In this paper, we present a purely algebraic proof that the (total) Springer representations of $\mathop{GL} ( n )$ are $\mathrm{Ext}$-orthogonal to each other, and show that it is compatible with the natural categorification of the ring of symmetric functions.


title: Arbitrary high-order linearly implicit energy-preserving algorithms for Hamiltonian PDEs

abstract: In this paper, we present a novel strategy to systematically construct linearly implicit energy-preserving schemes with arbitrary order of accuracy for Hamiltonian PDEs. Such novel strategy is based on the newly developed exponential scalar variable (ESAV) approach that can remove the bounded-from-blew restriction of nonlinear terms in the Hamiltonian functional and provides a totally explicit discretization of the auxiliary variable without computing extra inner products, which make it more effective and applicable than the traditional scalar auxiliary variable (SAV) approach. To achieve arbitrary high-order accuracy and energy preservation, we utilize the symplectic Runge-Kutta method for both solution variables and the auxiliary variable, where the values of internal stages in nonlinear terms are explicitly derived via an extrapolation from numerical solutions already obtained in the preceding calculation. A prediction-correction strategy is proposed to further improve the accuracy. Fourier pseudo-spectral method is then employed to obtain fully discrete schemes. Compared with the SAV schemes, the solution variables and the auxiliary variable in these ESAV schemes are now decoupled. Moreover, when the linear terms are of constant coefficients, the solution variables can be explicitly solved by using the fast Fourier transform. Numerical experiments are carried out for three Hamiltonian PDEs to demonstrate the efficiency and conservation of the ESAV schemes.


title: Stochastic Decomposition Method for Two-Stage Distributionally Robust Optimization

abstract: In this paper, we present a sequential sampling-based algorithm for the two-stage distributionally robust linear programming (2-DRLP) models. The 2-DRLP models are defined over a general class of ambiguity sets with discrete or continuous probability distributions. The algorithm is a distributionally robust version of the well-known stochastic decomposition algorithm of Higle and Sen (Math. of OR 16(3), 650-669, 1991) for a two-stage stochastic linear program. We refer to the algorithm as the distributionally robust stochastic decomposition (DRSD) method. The key features of the algorithm include (1) it works with data-driven approximations of ambiguity sets that are constructed using samples of increasing size and (2) efficient construction of approximations of the worst-case expectation function that solves only two second-stage subproblems in every iteration. We identify conditions under which the ambiguity set approximations converge to the true ambiguity sets and show that the DRSD method asymptotically identifies an optimal solution, with probability one. We also computationally evaluate the performance of the DRSD method for solving distributionally robust versions of instances considered in stochastic programming literature. The numerical results corroborate the analytical behavior of the DRSD method and illustrate the computational advantage over an external sampling-based decomposition approach (distributionally robust L-shaped method).


title: Cosmological dynamics of dark energy in scalar-torsion $f(T,\phi)$ gravity

abstract: It is investigated the cosmological dynamics of scalar-torsion $f(T,\phi)$ gravity as a dark energy model, where $T$ is the torsion scalar of teleparallel gravity and $\phi$ is a canonical scalar field. In this context, we are concerned with the phenomenology of the class of models with non-linear coupling to gravity and exponential potential. We obtain the critical points of the autonomous system, along with the stability conditions of each one of them and their cosmological properties. Particularly, we show the existence of new attractors with accelerated expansion, as well as, new scaling solutions in which the energy density of dark energy scales as the background fluid density, thus, defining the so-called scaling radiation and scaling matter epochs. The scaling solutions are saddle points, and therefore, the system exits these solutions to the current epoch of cosmic acceleration, towards an attractor point describing the dark energy-dominated era.


title: Numerical twist-even SU(1,1)-singlet solutions in open string field theory around the identity-based solution

abstract: Using the level truncation method, we construct numerical solutions, which are twist even and SU(1,1) singlet, in the theory around the Takahashi-Tanimoto identity-based solution (TT solution) with a real parameter $a$ in the framework of bosonic open string field theory. We find solutions corresponding to "double brane" and "ghost brane" solutions which were constructed by Kudrna and Schnabl in the conventional theory around the perturbative vacuum. Our solutions show somewhat similar $a$-dependence to tachyon vacuum and single brane solutions, which we found in the earlier works. In this sense, we might be able to expect that they are consistent with the conventional interpretation of $a$-dependence of the TT solution. We observe that numerical complex solutions at low levels become real ones at higher levels for some region of the parameter $a$. However, these real solutions do not so improve interpretation for double brane.


title: 5G Air-to-Ground Network Design and Optimization: A Deep Learning Approach

abstract: Direct air-to-ground (A2G) communications leveraging the fifth-generation (5G) new radio (NR) can provide high-speed broadband in-flight connectivity to aircraft in the sky. A2G network deployment entails optimizing various design parameters such as inter-site distances, number of sectors per site, and the up-tilt angles of sector antennas. The system-level design guidelines in the existing work on A2G network are rather limited. In this paper, a novel deep learning-based framework is proposed for efficient design and optimization of a 5G A2G network. The devised architecture comprises two deep neural networks (DNNs): the first DNN is used for approximating the 5G A2G network behavior in terms of user throughput, and the second DNN is developed as a function optimizer to find the throughput-optimal deployment parameters including antenna up-tilt angles and inter-site distances. Simulation results are provided to validate the proposed model and reveal system-level design insights.


title: Eigenvalue inequalities for the buckling problem of the drifting Laplacian of arbitrary order

abstract: In this paper, we investigate the buckling problem of the drifting Laplacian of arbitrary order on a bounded connected domain in complete smooth metric measure spaces (SMMSs) supporting a special function, and successfully get a general inequality for its eigenvalues. By applying this general inequality, if the complete SMMSs considered satisfy some curvature constraints, we can obtain a universal inequalities for eigenvalues of this buckling problem.


title: Optimal Accuracy-Time Trade-off for Deep Learning Services in Edge Computing Systems

abstract: With the increasing demand for computationally intensive services like deep learning tasks, emerging distributed computing platforms such as edge computing (EC) systems are becoming more popular. Edge computing systems have shown promising results in terms of latency reduction compared to the traditional cloud systems. However, their limited processing capacity imposes a trade-off between the potential latency reduction and the achieved accuracy in computationally-intensive services such as deep learning-based services. In this paper, we focus on finding the optimal accuracy-time trade-off for running deep learning services in a three-tier EC platform where several deep learning models with different accuracy levels are available. Specifically, we cast the problem as an Integer Linear Program, where optimal task scheduling decisions are made to maximize overall user satisfaction in terms of accuracy-time trade-off. We prove that our problem is NP-hard and then provide a polynomial constant-time greedy algorithm, called GUS, that is shown to attain near-optimal results. Finally, upon vetting our algorithmic solution through numerical experiments and comparison with a set of heuristics, we deploy it on a test-bed implemented to measure for real-world results. The results of both numerical analysis and real-world implementation show that GUS can outperform the baseline heuristics in terms of the average percentage of satisfied users by a factor of at least 50%.


title: Learning Efficient GANs for Image Translation via Differentiable Masks and co-Attention Distillation

abstract: Generative Adversarial Networks (GANs) have been widely-used in image translation, but their high computation and storage costs impede the deployment on mobile devices. Prevalent methods for CNN compression cannot be directly applied to GANs due to the peculiarties of GAN tasks and the unstable adversarial training. To solve these, in this paper, we introduce a novel GAN compression method, termed DMAD, by proposing a Differentiable Mask and a co-Attention Distillation. The former searches for a light-weight generator architecture in a training-adaptive manner. To overcome channel inconsistency when pruning the residual connections, an adaptive cross-block group sparsity is further incorporated. The latter simultaneously distills informative attention maps from both the generator and discriminator of a pre-trained model to the searched generator, effectively stabilizing the adversarial training of our light-weight model. Experiments show that DMAD can reduce the Multiply Accumulate Operations (MACs) of CycleGAN by 13x and that of Pix2Pix by 4x while retaining a comparable performance against the full model. Our code can be available at https://github.com/SJLeo/DMAD.


title: Boundary regularity and stability for spaces with Ricci bounded below

abstract: This paper studies the structure and stability of boundaries in noncollapsed $\text{RCD}(K,N)$ spaces, that is, metric-measure spaces $(X,\mathsf{d},\mathscr{H}^N)$ with lower Ricci curvature bounded below. Our main structural result is that the boundary $\partial X$ is homeomorphic to a manifold away from a set of codimension 2, and is $N-1$ rectifiable. Along the way we show effective measure bounds on the boundary and its tubular neighborhoods. These results are new even for Gromov-Hausdorff limits $(M_i^N,\mathsf{d}_{g_i},p_i) \rightarrow (X,\mathsf{d},p)$ of smooth manifolds with boundary, and require new techniques beyond those needed to prove the analogous statements for the regular set, in particular when it comes to the manifold structure of the boundary $\partial X$. The key local result is an $\epsilon$-regularity theorem, which tells us that if a ball $B_{2}(p)\subset X$ is sufficiently close to a half space $B_{2}(0)\subset \mathbb{R}^N_+$ in the Gromov-Hausdorff sense, then $B_1(p)$ is biH\"older to an open set of $\mathbb{R}^N_+$. In particular, $\partial X$ is itself homeomorphic to $B_1(0^{N-1})$ near $B_1(p)$. Further, the boundary $\partial X$ is $N-1$ rectifiable and the boundary measure $\mathscr{H}^{N-1}_{\partial X}$ is Ahlfors regular on $B_1(p)$ with volume close to the Euclidean volume. Our second collection of results involve the stability of the boundary with respect to noncollapsed mGH convergence $X_i\to X$. Specifically, we show a boundary volume convergence which tells us that the $N-1$ Hausdorff measures on the boundaries converge $\mathscr{H}^{N-1}_{\partial X_i}\to \mathscr{H}^{N-1}_{\partial X}$ to the limit Hausdorff measure on $\partial X$. We will see that a consequence of this is that if the $X_i$ are boundary free then so is $X$.


title: Optimal Sub-Gaussian Mean Estimation in $\mathbb{R}$

abstract: We revisit the problem of estimating the mean of a real-valued distribution, presenting a novel estimator with sub-Gaussian convergence: intuitively, "our estimator, on any distribution, is as accurate as the sample mean is for the Gaussian distribution of matching variance." Crucially, in contrast to prior works, our estimator does not require prior knowledge of the variance, and works across the entire gamut of distributions with bounded variance, including those without any higher moments. Parameterized by the sample size $n$, the failure probability $\delta$, and the variance $\sigma^2$, our estimator is accurate to within $\sigma\cdot(1+o(1))\sqrt{\frac{2\log\frac{1}{\delta}}{n}}$, tight up to the $1+o(1)$ factor. Our estimator construction and analysis gives a framework generalizable to other problems, tightly analyzing a sum of dependent random variables by viewing the sum implicitly as a 2-parameter $\psi$-estimator, and constructing bounds using mathematical programming and duality techniques.

categories: math.ST cs.DS cs.IT cs.LG math.IT stat.ML stat.TH


title: First constraints on the intrinsic CMB dipole and our velocity with Doppler and aberration

abstract: We test the usual hypothesis that the Cosmic Microwave Background (CMB) dipole, its largest anisotropy, is due to our peculiar velocity with respect to the Hubble flow by measuring independently the Doppler and aberration effects on the CMB using Planck 2018 data. We remove the spurious contributions from the conversion of intensity into temperature and arrive at measurements which are independent from the CMB dipole itself for both temperature and polarization maps and both SMICA and NILC component-separation methods. Combining these new measurements with the dipole one we get the first constraints on the intrinsic CMB dipole. Assuming a standard dipolar lensing contribution we can put an upper limit on the intrinsic amplitude: 3.7 mK (95% CI). We estimate the peculiar velocity of the solar system without assuming a negligible intrinsic dipole contribution: $v = (300^{+111}_{-93})$ km/s with $(l,b) = (276 \pm 33, 51 \pm 19)^\circ$ [SMICA], and $v = (296^{+111}_{-88})$ km/s with $(l,b) = (280 \pm 33, 50 \pm 20)^\circ$ [NILC] with negligible systematic contributions. These values are consistent with the peculiar velocity hypothesis of the dipole.


title: Partition-theoretic Frobenius-type limit formulas

abstract: Using partition generating function techniques, we prove $q$-series analogues of a formula of Frobenius generalizing Abel's convergence theorem for complex power series. Frobenius' result states that for $|q|<1$, $\lim_{q\to 1}(1-q)\sum_{n\geq 1} f(n) q^n $ is equal to the average value $\lim_{N\to \infty}$ $\frac{1}{N}\sum_{k=1}^{N}f(k)$ of the sequence $\{f(n)\}$ as $n\to \infty$, if the average value exists.


title: Field dependence of the Yukawa coupling in the three flavor quark-meson model

abstract: We investigate the renormalization group flow of the field-dependent Yukawa coupling in the framework of the three flavor quark-meson model. In a conventional perturbative calculation, given that the field rescaling is trivial, the Yukawa coupling does not get renormalized at the one-loop level if it is coupled to an equal number of scalar and pseudoscalar fields. Its field-dependent version, however, does flow with respect to the scale. Using the functional renormalization group technique, we show that it is highly nontrivial how to extract the actual flow of the Yukawa coupling as there are several new chirally invariant operators that get generated by quantum fluctuations in the effective action, which need to be distinguished from that of the Yukawa interaction.


title: Domain Adaptation based Interpretable Image Emotion Recognition using Facial Expression Recognition

abstract: A domain adaptation technique has been proposed in this paper to identify the emotions in generic images containing facial & non-facial objects and non-human components. It addresses the challenge of the insufficient availability of pre-trained models and well-annotated datasets for image emotion recognition (IER). It starts with proposing a facial emotion recognition (FER) system and then moves on to adapting it for image emotion recognition. First, a deep-learning-based FER system has been proposed that classifies a given facial image into discrete emotion classes. Further, an image recognition system has been proposed that adapts the proposed FER system to recognize the emotions portrayed by images using domain adaptation. It classifies the generic images into 'happy,' 'sad,' 'hate,' and 'anger' classes. A novel interpretability approach, Divide and Conquer based Shap (DnCShap), has also been proposed to interpret the highly relevant visual features for emotion recognition. The proposed system's architecture has been decided through ablation studies, and the experiments are conducted on four FER and four IER datasets. The proposed IER system has shown an emotion classification accuracy of 59.61% for the IAPSa dataset, 57.83% for the ArtPhoto dataset, 67.93% for the FI dataset, and 55.13% for the EMOTIC dataset. The important visual features leading to a particular emotion class have been identified, and the embedding plots for various emotion classes have been analyzed to explain the proposed system's predictions.


title: Model selection for count timeseries with applications in forecasting number of trips in bike-sharing systems and its volatility

abstract: Forecasting the number of trips in bike-sharing systems and its volatility over time is crucial for planning and optimizing such systems. This paper develops timeseries models to forecast hourly count timeseries data, and estimate its volatility. Such models need to take into account the complex patterns over various temporal scales including hourly, daily, weekly and annual as well as the temporal correlation. To capture this complex structure, a large number of parameters are needed. Here a structural model selection approach is utilized to choose the parameters. This method explores the parameter space for a group of covariates at each step. These groups of covariate are constructed to represent a particular structure in the model. The statistical models utilized are extensions of Generalized Linear Models to timeseries data. One challenge in using such models is the explosive behavior of the simulated values. To address this issue, we develop a technique which relies on damping the simulated value, if it falls outside of an admissible interval. The admissible interval is defined using measures of variability of the left and right tails. A new definition of outliers is proposed based on these variability measures. This new definition is shown to be useful in the context of asymmetric distributions.


title: Direct growth of germanene at interfaces between van der Waals materials and Ag(111)

abstract: Germanene, a two-dimensional honeycomb germanium crystal, is grown at graphene/Ag(111) and hexagonal boron nitride (h-BN)/Ag(111) interfaces by segregating germanium atoms. A simple annealing process in N2 or H2/Ar at ambient pressure leads to the formation of germanene, indicating that an ultrahigh-vacuum condition is not necessary. The grown germanene is stable in air and uniform over the entire area covered with a van der Waals (vdW) material. As an important finding, it is necessary to use a vdW material as a cap layer for the present germanene growth method since the use of an Al2O3 cap layer resulted in no germanene formation. The present study also proved that Raman spectroscopy in air is a powerful tool for characterizing germanene at the interfaces, which is concluded by multiple analyses including first-principles density functional theory calculations. The direct growth of h-BN-capped germanene on Ag(111), which is demonstrated in the present study, is considered to be a promising technique for the fabrication of future germanene-based electronic devices.


title: Ultra-efficient magnetism modulation in a Weyl ferromagnet by current-assisted domain wall motion

abstract: Flexible and efficient manipulation of magnetic configurations can be challenging. In the design of practical devices, achieving a high effective magnetic field with a low working current is under tight demand. Here, we report a unique method for efficient magnetism modulation by direct current injection in magnetic Weyl semimetal Co3Sn2S2. We demonstrate that the modulation process stems from current-assisted domain wall motion. Through two independent methods, we reveal that the spin-transfer torque efficiency of Co3Sn2S2 reaches as high as 2.4-5.6 kOe MA^(-1) cm^2, and the threshold current density for driving the magnetic domain walls is as low as <5.1*10^5 A/cm^2 without an external field, and <1.5*10^5 A/cm^2 with a moderate external field. Our findings manifest a new and powerful approach for sub-micron magnetism manipulation, and also open the door towards a new paradigm of spintronics that combines magnetism, topology, and metallicity for low-energy consumption memory and computing.


title: Laplace Green's Functions for Infinite Ground Planes with Local Roughness

abstract: The Green's functions for the Laplace equation respectively satisfying the Dirichlet and Neumann boundary conditions on the upper side of an infinite plane with a circular hole are introduced and constructed. These functions enables solution of the boundary value problems in domains where the hole is closed by any surface. This approach enables accounting for arbitrary positive and negative ground elevations inside the domain of interest, which, generally, is not possible to achieve using the regular method of images. Such problems appear in electrostatics, however, the methods developed apply to other domains where the Laplace or Poisson equations govern. Integral and series representations of the Green's functions are provided. An efficient computational technique based on the boundary element method with fast multipole acceleration is developed. A numerical study of some benchmark problems is presented.

categories: math.NA cs.NA math-ph math.MP physics.comp-ph


title: Data-aided Sensing for Distributed Detection

abstract: In this paper, we study data-aided sensing (DAS) for distributed detection in wireless sensor networks (WSNs) when sensors' measurements are correlated. In particular, we derive a node selection criterion based on the J-divergence in DAS for reliable decision subject to a decision delay constraint. Based on the proposed J-divergence based DAS, the nodes can be selected to rapidly increase the log-likelihood ratio (LLR), which leads to a reliable decision with a smaller number of the sensors that upload measurements for a shorter decision delay. From simulation results, it is confirmed that the J-divergence based DAS can provide a reliable decision with a smaller number of sensors compared to other approaches.


title: Each Thurston geometry admits a Heegaard genus two example

abstract: We collect examples of genus two Heegaard diagrams for compact 3-manifolds admitting each of the eight Thurston geometries.


title: Joint Resource Optimization for IRS-Assisted mmWave MIMO under QoS Constraints

abstract: This letter focuses on the non-convex joint optimization with a dynamic resource of multi-user for an intelligent reflecting surface-enhanced mmWave system, where all users are concentrated on the unique cluster beam. Firstly, the objective function of the above non-linear problem is converted into a quadratic programming form under the quality of service constraints. Further, a multi-blocks alternating optimization framework with dynamic power allocation is proposed to obtain the maximum sum-rate, where the relaxed ADMM algorithm is adopted to tackle the optimal fulldigital precoder and the corresponding passive reflecting matrix is obtained by the gradient-projection. The numerical results verify that beam optimization should be emphasized in high SINR, but joint dynamic resource allocation can further improve system performance even if the hardware dimensions reaches the limit.


title: Three-dimensional matter-wave interferometry of a trapped single ion

abstract: We report on a demonstration of Ramsey interferometry by three-dimensional motion with a trapped $^{171}Yb^+$ ion. We applied a momentum kick to the ion in a direction diagonal to the trap axes to initiate three-dimensional motion using a mode-locked pulse laser. The interference signal was analyzed theoretically to demonstrate three-dimensional matter-wave interference. This work paves the way to realizing matter-wave interferometry using trapped ions.


title: Ultra-Lightweight Speech Separation via Group Communication

abstract: Model size and complexity remain the biggest challenges in the deployment of speech enhancement and separation systems on low-resource devices such as earphones and hearing aids. Although methods such as compression, distillation and quantization can be applied to large models, they often come with a cost on the model performance. In this paper, we provide a simple model design paradigm that explicitly designs ultra-lightweight models without sacrificing the performance. Motivated by the sub-band frequency-LSTM (F-LSTM) architectures, we introduce the group communication (GroupComm), where a feature vector is split into smaller groups and a small processing block is used to perform inter-group communication. Unlike standard F-LSTM models where the sub-band outputs are concatenated, an ultra-small module is applied on all the groups in parallel, which allows a significant decrease on the model size. Experiment results show that comparing with a strong baseline model which is already lightweight, GroupComm can achieve on par performance with 35.6 times fewer parameters and 2.3 times fewer operations.


title: Augmented Fairness: An Interpretable Model Augmenting Decision-Makers' Fairness

abstract: We propose a model-agnostic approach for mitigating the prediction bias of a black-box decision-maker, and in particular, a human decision-maker. Our method detects in the feature space where the black-box decision-maker is biased and replaces it with a few short decision rules, acting as a "fair surrogate". The rule-based surrogate model is trained under two objectives, predictive performance and fairness. Our model focuses on a setting that is common in practice but distinct from other literature on fairness. We only have black-box access to the model, and only a limited set of true labels can be queried under a budget constraint. We formulate a multi-objective optimization for building a surrogate model, where we simultaneously optimize for both predictive performance and bias. To train the model, we propose a novel training algorithm that combines a nondominated sorting genetic algorithm with active learning. We test our model on public datasets where we simulate various biased "black-box" classifiers (decision-makers) and apply our approach for interpretable augmented fairness.


title: Efficient and Effective Community Search on Large-scale Bipartite Graphs

abstract: Bipartite graphs are widely used to model relationships between two types of entities. Community search retrieves densely connected subgraphs containing a query vertex, which has been extensively studied on unipartite graphs. However, community search on bipartite graphs remains largely unexplored. Moreover, all existing cohesive subgraph models on bipartite graphs can only be applied to measure the structure cohesiveness between two sets of vertices while overlooking the edge weight in forming the community. In this paper, we study the significant (alpha, beta)-community search problem on weighted bipartite graphs. Given a query vertex q, we aim to find the significant (alpha, beta)-community R of q which adopts (alpha, beta)-core to characterize the engagement level of vertices, and maximizes the minimum edge weight (significance) within R. To support fast retrieval of R, we first retrieve the maximal connected subgraph of (alpha, beta)-core containing the query vertex (the (alpha, beta)-community), and the search space is limited to this subgraph with a much smaller size than the original graph. A novel index structure is presented which can be built in O(delta * m) time and takes O(delta * m) space where m is the number of edges in G, delta is bounded by the square root of m and is much smaller in practice. Utilizing the index, the (alpha, beta)-community can be retrieved in optimal time. To further obtain R, we develop peeling and expansion algorithms to conduct searches by shrinking from the (alpha, beta)-community and expanding from the query vertex, respectively. The experimental results on real graphs not only demonstrate the effectiveness of the significant (alpha, beta)-community model but also validate the efficiency of our query processing and indexing techniques.


title: Rethinking the Separation Layers in Speech Separation Networks

abstract: Modules in all existing speech separation networks can be categorized into single-input-multi-output (SIMO) modules and single-input-single-output (SISO) modules. SIMO modules generate more outputs than input, and SISO modules keep the numbers of input and output the same. While the majority of separation models only contain SIMO architectures, it has also been shown that certain two-stage separation systems integrated with a post-enhancement SISO module can improve the separation quality. Why performance improvements can be achieved by incorporating the SISO modules? Are SIMO modules always necessary? In this paper, we empirically examine those questions by designing models with varying configurations in the SIMO and SISO modules. We show that comparing with the standard SIMO-only design, a mixed SIMO-SISO design with a same model size is able to improve the separation performance especially under low-overlap conditions. We further validate the necessity of SIMO modules and show that SISO-only models are still able to perform separation without sacrificing the performance. The observations allow us to rethink the model design paradigm and present different views on how the separation is performed.


title: Implicit Filter-and-sum Network for Multi-channel Speech Separation

abstract: Various neural network architectures have been proposed in recent years for the task of multi-channel speech separation. Among them, the filter-and-sum network (FaSNet) performs end-to-end time-domain filter-and-sum beamforming and has shown effective in both ad-hoc and fixed microphone array geometries. In this paper, we investigate multiple ways to improve the performance of FaSNet. From the problem formulation perspective, we change the explicit time-domain filter-and-sum operation which involves all the microphones into an implicit filter-and-sum operation in the latent space of only the reference microphone. The filter-and-sum operation is applied on a context around the frame to be separated. This allows the problem formulation to better match the objective of end-to-end separation. From the feature extraction perspective, we modify the calculation of sample-level normalized cross correlation (NCC) features into feature-level NCC (fNCC) features. This makes the model better matches the implicit filter-and-sum formulation. Experiment results on both ad-hoc and fixed microphone array geometries show that the proposed modification to the FaSNet, which we refer to as iFaSNet, is able to significantly outperform the benchmark FaSNet across all conditions with an on par model complexity.


title: Time-Resolved Focused Ion Beam Microscopy: Modeling, Estimation Methods, and Analyses

abstract: In a focused ion beam (FIB) microscope, source particles interact with a small volume of a sample to generate secondary electrons that are detected, pixel by pixel, to produce a micrograph. Randomness of the number of incident particles causes excess variation in the micrograph, beyond the variation in the underlying particle-sample interaction. We recently demonstrated that joint processing of multiple time-resolved measurements from a single pixel can mitigate this effect of source shot noise in helium ion microscopy. This paper is focused on establishing a rigorous framework for understanding the potential for this approach. It introduces idealized continuous- and discrete-time abstractions of FIB microscopy with direct electron detection and estimation-theoretic limits of imaging performance under these measurement models. Novel estimators for use with continuous-time measurements are introduced and analyzed, and estimators for use with discrete-time measurements are analyzed and shown to approach their continuous-time counterparts as time resolution is increased. Simulated FIB microscopy results are consistent with theoretical analyses and demonstrate that substantial improvements over conventional FIB microscopy image formation are made possible by time-resolved measurement.


title: Large and moderate deviation principles for McKean-Vlasov SDEs with jumps

abstract: In this paper, we consider McKean-Vlasov stochastic differential equations (MVSDEs) driven by L\'evy noise. By identifying the right equations satisfied by the solutions of the MVSDEs with shifted driving L\'evy noise, we build up a framework to fully apply the weak convergence method to establish large and moderate deviation principles for MVSDEs. In the case of ordinary SDEs, the rate function is calculated by using the solutions of the corresponding skeleton equations simply replacing the noise by the elements of the Cameron-Martin space. It turns out that the correct rate function for MVSDEs is defined through the solutions of skeleton equations replacing the noise by smooth functions and replacing the distributions involved in the equation by the distribution of the solution of the corresponding deterministic equation (without the noise). This is somehow surprising. With this approach, we obtain large and moderate deviation principles for much wider classes of MVSDEs in comparison with the existing literature.


title: Regularity via Links and Stein Factorization

abstract: Here, we introduce a new definition of regular point for piecewise-linear (PL) functions on combinatorial (PL triangulated) manifolds. This definition is given in terms of the restriction of the function to the link of the point. We show that our definition of regularity is distinct from other definitions that exist in the combinatorial topology literature. Next, we stratify the Jacobi set/critical locus of such a map as a poset stratified space. As an application, we consider the Reeb space of a PL function, stratify the Reeb space as well as the target of the function, and show that the Stein factorization is a map of stratified spaces.


title: Peer groups for organisational learning: clustering with practical constraints

abstract: Peer-grouping is used in many sectors for organisational learning, policy implementation, and benchmarking. Clustering provides a statistical, data-driven method for constructing meaningful peer groups, but peer groups must be compatible with business constraints such as size and stability considerations. Additionally, statistical peer groups are constructed from many different variables, and can be difficult to understand, especially for non-statistical audiences. We developed methodology to apply business constraints to clustering solutions and allow the decision-maker to choose the balance between statistical goodness-of-fit and conformity to business constraints. Several tools were utilised to identify complex distinguishing features in peer groups, and a number of visualisations are developed to explain high-dimensional clusters for non-statistical audiences. In a case study where peer group size was required to be small ($\leq 100$ members), we applied constrained clustering to a noisy high-dimensional data-set over two subsequent years, ensuring that the clusters were sufficiently stable between years. Our approach not only satisfied clustering constraints on the test data, but maintained an almost monotonic negative relationship between goodness-of-fit and stability between subsequent years. We demonstrated in the context of the case study how distinguishing features between clusters can be communicated clearly to different stakeholders with substantial and limited statistical knowledge.


title: Reinforcement Learning of Graph Neural Networks for Service Function Chaining

abstract: In the management of computer network systems, the service function chaining (SFC) modules play an important role by generating efficient paths for network traffic through physical servers with virtualized network functions (VNF). To provide the highest quality of services, the SFC module should generate a valid path quickly even in various network topology situations including dynamic VNF resources, various requests, and changes of topologies. The previous supervised learning method demonstrated that the network features can be represented by graph neural networks (GNNs) for the SFC task. However, the performance was limited to only the fixed topology with labeled data. In this paper, we apply reinforcement learning methods for training models on various network topologies with unlabeled data. In the experiments, compared to the previous supervised learning method, the proposed methods demonstrated remarkable flexibility in new topologies without re-designing and re-training, while preserving a similar level of performance.


title: A statistical machine learning approach for benchmarking in the presence of complex contextual factors and peer groups

abstract: The ability to compare between individuals or organisations fairly is important for the development of robust and meaningful quantitative benchmarks. To make fair comparisons, contextual factors must be taken into account, and comparisons should only be made between similar organisations such as peer groups. Previous benchmarking methods have used linear regression to adjust for contextual factors, however linear regression is known to be sub-optimal when nonlinear relationships exist between the comparative measure and covariates. In this paper we propose a random forest model for benchmarking that can adjust for these potential nonlinear relationships, and validate the approach in a case-study of high noise data. We provide new visualisations and numerical summaries of the fitted models and comparative measures to facilitate interpretation by both analysts and non-technical audiences. Comparisons can be made across the cohort or within peer groups, and bootstrapping provides a means of estimating uncertainty in both adjusted measures and rankings. We conclude that random forest models can facilitate fair comparisons between organisations for quantitative measures including in cases on complex contextual factor relationships, and that the models and outputs are readily interpreted by stakeholders.