title: Classes of critical avalanche dynamics in complex networks

abstract: Dynamical processes exhibiting absorbing states are essential in the modeling of a large variety of situations from material science to epidemiology and social sciences. Such processes exhibit the possibility of avalanching behavior upon slow driving. Here, we study the distribution of sizes and durations of avalanches for well-known dynamical processes on complex networks. We find that all analyzed models display a similar critical behavior, characterized by the presence of two distinct regimes. At small scales, sizes and durations of avalanches exhibit distributions that are dependent on the network topology and the model dynamics. At asymptotically large scales instead -- irrespective of the type of dynamics and of the topology of the underlying network -- sizes and durations of avalanches are characterized by power-law distributions with the exponents of the standard mean-field critical branching process.


title: On the quasicompactness of the moduli stack of logarithmic G-connections over a curve

abstract: Fix a smooth projective curve over a field of characteristic zero and a finite set of punctures. Let G be a connected linear algebraic group. We prove that the moduli of G-bundles with logarithmic connections having fixed residue classes at the punctures is an algebraic stack of finite type.


title: Policy Evaluation Networks

abstract: Many reinforcement learning algorithms use value functions to guide the search for better policies. These methods estimate the value of a single policy while generalizing across many states. The core idea of this paper is to flip this convention and estimate the value of many policies, for a single set of states. This approach opens up the possibility of performing direct gradient ascent in policy space without seeing any new data. The main challenge for this approach is finding a way to represent complex policies that facilitates learning and generalization. To address this problem, we introduce a scalable, differentiable fingerprinting mechanism that retains essential policy information in a concise embedding. Our empirical results demonstrate that combining these three elements (learned Policy Evaluation Network, policy fingerprints, gradient ascent) can produce policies that outperform those that generated the training data, in zero-shot manner.


title: Understanding How and Why University Students Use Virtual Private Networks

abstract: We study how and why university students chose and use VPNs, and whether they are aware of the security and privacy risks that VPNs pose. To answer these questions, we conducted 32 in-person interviews and a survey with 349 respondents, all university students in the United States. We find students are mostly concerned with access to content and privacy concerns were often secondary. They made tradeoffs to achieve a particular goal, such as using a free commercial VPN that may collect their online activities to access an online service in a geographic area. Many users expected that their VPNs were collecting data about them, although they did not understand how VPNs work. We conclude with a discussion of ways to help users make choices about VPNs.


title: Tensor Decompositions in Deep Learning

abstract: The paper surveys the topic of tensor decompositions in modern machine learning applications. It focuses on three active research topics of significant relevance for the community. After a brief review of consolidated works on multi-way data analysis, we consider the use of tensor decompositions in compressing the parameter space of deep learning models. Lastly, we discuss how tensor methods can be leveraged to yield richer adaptive representations of complex data, including structured information. The paper concludes with a discussion on interesting open research challenges.


title: No computation without representation: Avoiding data and algorithm biases through diversity

abstract: The emergence and growth of research on issues of ethics in AI, and in particular algorithmic fairness, has roots in an essential observation that structural inequalities in society are reflected in the data used to train predictive models and in the design of objective functions. While research aiming to mitigate these issues is inherently interdisciplinary, the design of unbiased algorithms and fair socio-technical systems are key desired outcomes which depend on practitioners from the fields of data science and computing. However, these computing fields broadly also suffer from the same under-representation issues that are found in the datasets we analyze. This disconnect affects the design of both the desired outcomes and metrics by which we measure success. If the ethical AI research community accepts this, we tacitly endorse the status quo and contradict the goals of non-discrimination and equity which work on algorithmic fairness, accountability, and transparency seeks to address. Therefore, we advocate in this work for diversifying computing as a core priority of the field and our efforts to achieve ethical AI practices. We draw connections between the lack of diversity within academic and professional computing fields and the type and breadth of the biases encountered in datasets, machine learning models, problem formulations, and interpretation of results. Examining the current fairness/ethics in AI literature, we highlight cases where this lack of diverse perspectives has been foundational to the inequity in treatment of underrepresented and protected group data. We also look to other professional communities, such as in law and health, where disparities have been reduced both in the educational diversity of trainees and among professional practices. We use these lessons to develop recommendations that provide concrete steps for the computing community to increase diversity.


title: Full Duplex Hybrid A/D Beamforming with Reduced Complexity Multi-Tap Analog Cancellation

abstract: Although the hardware complexity of the analog self-interference canceller in full duplex Multiple Input Multiple Output (MIMO) designs does not necessarily scale with the number of transceiver antennas, exploiting the benefits of analog cancellation in massive MIMO systems with hundreds of antenna elements is still quite impractical. Hybrid Analog and Digital (A/D) beamforming architectures have been lately considered as a candidate technology for realizing massive MIMO transceivers with very large number of antenna elements, but with much fewer numbers of Radio Frequency (RF) chains. In this paper, we present a novel architecture for full duplex hybrid A/D beamforming transceivers including multi-tap analog cancellation with reduced number of taps and simple multiplexers for efficient signal routing among the transceiver RF chains. Capitalizing on the proposed transceiver architecture, we present a joint design of analog cancellation and A/D beamforming with the objective to maximize the achievable full duplex rate performance. Representative millimeter wave simulation results demonstrate the effectiveness of the proposed architecture and algorithmic framework for enabling simultaneous uplink and downlink communications with reduced complexity analog self-interference cancellation.


title: Presence and absence of itinerant gapless excitations in the quantum spin liquid candidate EtMe$_3$Sb[Pd(dmit)$_2$]$_2$

abstract: EtMe$_3$Sb[Pd(dmit)$_2$]$_2$, an organic Mott insulator with nearly isotropic triangular lattice, is a candidate material for a quantum spin liquid, in which the zero-point fluctuations do not allow the spins to order. The itinerant gapless excitations inferred from the thermal transport measurements in this system have been a hotly debated issue recently. While the presence of a finite linear residual thermal conductivity, $\kappa_0/T \equiv \kappa/T (T \rightarrow 0)$, has been shown [M. Yamashita {\it et al.} Science {\bf 328}, 1246 (2010)], recent experiments [P. Bourgeois-Hope {\it et al.}, Phys. Rev. X {\bf 9}, 041051 (2019); J. M. Ni {\it et al.}, Phys. Rev. Lett. {\bf 123}, 247204 (2019)] have reported the absence of $\kappa_0/T$. Here we show that the low-temperature thermal conductivity strongly depends on the cooling process of the sample. When cooling down very slowly, a sizable $\kappa_0/T$ is observed. In contrast, when cooling down rapidly, $\kappa_0/T$ vanishes and, in addition, the phonon thermal conductivity is strongly suppressed. These results suggest that possible random scatterers introduced during the cooling process are responsible for the apparent discrepancy of the thermal conductivity data in this organic system. The present results provide evidence that the true ground state of EtMe$_3$Sb[Pd(dmit)$_2$]$_2$ is likely to be a quantum spin liquid with itinerant gapless excitations.


title: Two-neutron knockout as a probe of the composition of states in $^{22}$Mg, $^{23}$Al, and $^{24}$Si

abstract: Simpson and Tostevin proposed that the width and shape of exclusive parallel momentum distributions of the A-2 residue in direct two-nucleon knockout reactions carry a measurable sensitivity to the nucleon single-particle configurations and their couplings within the wave functions of exotic nuclei. We report here on the first benchmarks and use of this new spectroscopic tool. Exclusive parallel momentum distributions for states in the neutron-deficient nuclei $^{22}$Mg, $^{23}$Al, and $^{24}$Si populated in such direct two-neutron removal reactions were extracted and compared to predictions combining eikonal reaction theory and shell-model calculations. For the well-known $^{22}$Mg and $^{23}$Al nuclei, measurements and calculations were found to agree, supporting the dependence of the parallel momentum distribution width on the angular momentum composition of the shell-model two-neutron amplitudes. In $^{24}$Si, a level at 3439(9) keV, of relevance for the important $^{23}$Al(p,$\gamma$)$^{24}$Si astrophysical reaction rate, was confirmed to be the $2^+_2$ state, while the $4^+_1$ state, expected to be strongly populated in two-neutron knockout, was not observed. This puzzle is resolved by theoretical considerations of the Thomas-Ehrman shift, which also suggest that a previously reported 3471-keV state in $^{24}$Si is in fact the ($0^+_2$) level with one of the largest experimental mirror-energy shifts ever observed.


title: Light curves from highly compact neutron stars with spot size effect

abstract: We systematically examine light curves from a single hot spot on a slowly rotating neutron star with very high compactness, where the so-called invisible zone does not exist. In particular, we adopt three different shapes of hot spot and take into account the finite size effect of hot spot on the light curves. Then, we find that the brightening of flux occurs when the hot spot with small area crosses the opposite side to the observer, where the brightening becomes stronger as the area of hot spot decreases. Since this brightening happens only when a part of hot spot crosses the opposite side to the observer, one may constrain the neutron star geometry, i.e., the combination of an inclination angle, angle between the rotational axis and normal vector at the center of hot spot, and the opening angle of hot spot, if the brightening would be observed. In addition, by counting such a brightening in the light curves, one may know how many narrow bands of hot spot and/or spots with small area cross the opposite side to the observer.


title: Towards Universal Representation Learning for Deep Face Recognition

abstract: Recognizing wild faces is extremely hard as they appear with all kinds of variations. Traditional methods either train with specifically annotated variation data from target domains, or by introducing unlabeled target variation data to adapt from the training data. Instead, we propose a universal representation learning framework that can deal with larger variation unseen in the given training data without leveraging target domain knowledge. We firstly synthesize training data alongside some semantically meaningful variations, such as low resolution, occlusion and head pose. However, directly feeding the augmented data for training will not converge well as the newly introduced samples are mostly hard examples. We propose to split the feature embedding into multiple sub-embeddings, and associate different confidence values for each sub-embedding to smooth the training procedure. The sub-embeddings are further decorrelated by regularizing variation classification loss and variation adversarial loss on different partitions of them. Experiments show that our method achieves top performance on general face recognition datasets such as LFW and MegaFace, while significantly better on extreme benchmarks such as TinyFace and IJB-S.


title: Edge-Matching Graph Contractions and their Interlacing Properties

abstract: For a given graph $\mathcal{G}$ of order $n$ with $m$ edges, and a real symmetric matrix associated to the graph, $M\left(\mathcal{G}\right)\in\mathbb{R}^{n\times n}$, the interlacing graph reduction problem is to find a graph $\mathcal{G}_{r}$ of order $r<n$ such that the eigenvalues of $M\left(\mathcal{G}_{r}\right)$ interlace the eigenvalues of $M\left(\mathcal{G}\right)$. Graph contractions over partitions of the vertices are widely used as a combinatorial graph reduction tool. In this study, we derive a graph reduction interlacing theorem based on subspace mappings and the minmax theory. We then define a class of edge-matching graph contractions and show how two types of edge-matching contractions provide Laplacian and normalized Laplacian interlacing. An $\mathcal{O}\left(mn\right)$ algorithm is provided for finding a normalized Laplacian interlacing contraction and an $\mathcal{O}\left(n^{2}+nm\right)$ algorithm is provided for finding a Laplacian interlacing contraction.


title: A Deep Unsupervised Feature Learning Spiking Neural Network with Binarized Classification Layers for EMNIST Classification using SpykeFlow

abstract: End user AI is trained on large server farms with data collected from the users. With ever increasing demand for IOT devices, there is a need for deep learning approaches that can be implemented (at the edge) in an energy efficient manner. In this work we approach this using spiking neural networks. The unsupervised learning technique of spike timing dependent plasticity (STDP) using binary activations are used to extract features from spiking input data. Gradient descent (backpropagation) is used only on the output layer to perform the training for classification. The accuracies obtained for the balanced EMNIST data set compare favorably with other approaches. The effect of stochastic gradient descent (SGD) approximations on learning capabilities of our network are also explored.


title: The hypergeometric test performs comparably to TF-IDF on standard text analysis tasks

abstract: Term frequency-inverse document frequency, or TF-IDF for short, and its many variants form a class of term weighting functions the members of which are widely used in text analysis applications. While TF-IDF was originally proposed as a heuristic, theoretical justifications grounded in information theory, probability, and the divergence from randomness paradigm have been advanced. In this work, we present an empirical study showing that TF-IDF corresponds very nearly with the hypergeometric test of statistical significance on selected real-data document retrieval, summarization, and classification tasks. These findings suggest that a fundamental mathematical connection between TF-IDF and the negative logarithm of the hypergeometric test P-value (i.e., a hypergeometric distribution tail probability) remains to be elucidated. We advance the empirical analyses herein as a first step toward explaining the long-standing effectiveness of TF-IDF from a statistical significance testing lens. It is our aspiration that these results will open the door to the systematic evaluation of significance testing derived term weighting functions in text analysis applications.


title: Demonstration of NV-detected ESR spectroscopy at 115 GHz and 4.2 Tesla

abstract: High frequency electron spin resonance (ESR) spectroscopy is an invaluable tool for identification and characterization of spin systems. Nanoscale ESR using the nitrogen-vacancy (NV) center has been demonstrated down to the level of a single spin. However, NV-detected ESR has exclusively been studied at low magnetic fields, where spectral overlap prevents clear identification of spectral features. Within this work, we demonstrate NV-detected ESR measurements of single-substitutional nitrogen impurities in diamond at a NV Larmor frequency of 115 GHz and the corresponding magnetic field of 4.2 Tesla. The NV-ESR measurements utilize a double electron-electron resonance sequence and are performed using both ensemble and single NV spin systems. In the single NV experiment, chirp pulses are used to improve the population transfer and for NV-ESR measurements. This work provides the basis for NV-based ESR measurements of external spins at high magnetic fields.


title: Causal inference with limited resources: proportionally-representative interventions

abstract: Investigators often evaluate treatment effects by considering settings in which all individuals are assigned a treatment of interest, assuming that an unlimited number of treatment units are available. However, many real-life treatments are of limited supply and cannot be provided to all individuals in the population. For example, patients on the liver transplant waiting list cannot be assigned a liver transplant immediately at the time they reach highest priority because a suitable organ is not likely to be immediately available. In these cases, investigators may still be interested in the effects of treatment strategies in which a finite number of organs are available at a given time, that is, treatment regimes that satisfy resource constraints. Here, we describe an estimand that can be used to define causal effects of treatment strategies that satisfy resource constraints: proportionally-representative interventions for limited resources. We derive a simple class of inverse probability weighted estimators, and apply one such estimator to evaluate the effect of restricting or expanding utilization of "increased risk" liver organs to treat patients with end-stage liver disease. Our method is designed to evaluate policy-relevant interventions in the setting of finite treatment resources.


title: Echo State Neural Machine Translation

abstract: We present neural machine translation (NMT) models inspired by echo state network (ESN), named Echo State NMT (ESNMT), in which the encoder and decoder layer weights are randomly generated then fixed throughout training. We show that even with this extremely simple model construction and training procedure, ESNMT can already reach 70-80% quality of fully trainable baselines. We examine how spectral radius of the reservoir, a key quantity that characterizes the model, determines the model behavior. Our findings indicate that randomized networks can work well even for complicated sequence-to-sequence prediction NLP tasks.


title: Analysis of diversity-accuracy tradeoff in image captioning

abstract: We investigate the effect of different model architectures, training objectives, hyperparameter settings and decoding procedures on the diversity of automatically generated image captions. Our results show that 1) simple decoding by naive sampling, coupled with low temperature is a competitive and fast method to produce diverse and accurate caption sets; 2) training with CIDEr-based reward using Reinforcement learning harms the diversity properties of the resulting generator, which cannot be mitigated by manipulating decoding parameters. In addition, we propose a new metric AllSPICE for evaluating both accuracy and diversity of a set of captions by a single value.


title: Comparison of Multi-Class and Binary Classification Machine Learning Models in Identifying Strong Gravitational Lenses

abstract: Typically, binary classification lens-finding schemes are used to discriminate between lens candidates and non-lenses. However, these models often suffer from substantial false-positive classifications. Such false positives frequently occur due to images containing objects such as crowded sources, galaxies with arms, and also images with a central source and smaller surrounding sources. Therefore, a model might confuse the stated circumstances with an Einstein ring. It has been proposed that by allowing such commonly misclassified image types to constitute their own classes, machine learning models will more easily be able to learn the difference between images that contain real lenses, and images that contain lens imposters. Using Hubble Space Telescope (HST) images, in the F814W filter, we compare the usage of binary and multi-class classification models applied to the lens finding task. From our findings, we conclude there is not a significant benefit to using the multi-class model over a binary model. We will also present the results of a simple lens search using a multi-class machine learning model, and potential new lens candidates.


title: Joint Optimization of Signal Design and Resource Allocation in Wireless D2D Edge Computing

abstract: In this paper, we study the distributed computational capabilities of device-to-device (D2D) networks. A key characteristic of D2D networks is that their topologies are reconfigurable to cope with network demands. For distributed computing, resource management is challenging due to limited network and communication resources, leading to inter-channel interference. To overcome this, recent research has addressed the problems of wireless scheduling, subchannel allocation, power allocation, and multiple-input multiple-output (MIMO) signal design, but has not considered them jointly. In this paper, unlike previous mobile edge computing (MEC) approaches, we propose a joint optimization of wireless MIMO signal design and network resource allocation to maximize energy efficiency. Given that the resulting problem is a non-convex mixed integer program (MIP) which is prohibitive to solve at scale, we decompose its solution into two parts: (i) a resource allocation subproblem, which optimizes the link selection and subchannel allocations, and (ii) MIMO signal design subproblem, which optimizes the transmit beamformer, transmit power, and receive combiner. Simulation results using wireless edge topologies show that our method yields substantial improvements in energy efficiency compared with cases of no offloading and partially optimized methods and that the efficiency scales well with the size of the network.


title: Quantum geometry, logic and probability

abstract: Quantum geometry on a discrete set means a directed graph with a weight associated to each arrow defining the quantum metric. However, these `lattice spacing' weights do not have to be independent of the direction of the arrow. We use this greater freedom to give a quantum geometric interpretation of discrete Markov processes with transition probabilities as arrow weights, namely taking the diffusion form $\partial_+ f=(-\Delta_\theta+ q-p)f$ for the graph Laplacian $\Delta_\theta$, potential functions $q,p$ built from the probabilities, and finite difference $\partial_+$ in the time direction. Motivated by this new point of view, we introduce a `discrete Schroedinger process' as $\partial_+\psi=\imath(-\Delta+V)\psi$ for the Laplacian associated to a bimodule connection such that the discrete evolution is unitary. We solve this explicitly for the 2-state graph, finding a 1-parameter family of such connections and an induced `generalised Markov process' for $f=|\psi|^2$ in which there is an additional source current built from $\psi$. We also discuss our recent work on the quantum geometry of logic in `digital' form over the field $\Bbb F_2=\{0,1\}$, including de Morgan duality and its possible generalisations.


title: A multiscale scheme accurately simulates macroscale shocks in an equation-free framework

abstract: Scientists and engineers often create accurate, trustworthy, computational simulation schemes - but all too often these are too computationally expensive to execute over the time or spatial domain of interest. The equation-free approach is to marry such trusted simulations to a framework for numerical macroscale reduction - the patch dynamics scheme. This article extends the patch scheme to scenarios in which the trusted simulation resolves abrupt state changes on the microscale that appear as shocks on the macroscale. Accurate simulation for problems in these scenarios requires extending the patch scheme by capturing the shock within a novel patch, and also modifying the patch coupling rules in the vicinity in order to maintain accuracy. With these two extensions to the patch scheme, straightforward arguments derive consistency conditions that match the usual order of accuracy for patch schemes. The new scheme is successfully tested on four archetypal problems. This technique will empower scientists and engineers to accurately and efficiently simulate, over large spatial domains, multiscale multiphysics systems that have rapid transition layers on the microscale.


title: Posterior Sampling for Anytime Motion Planning on Graphs with Expensive-to-Evaluate Edges

abstract: Collision checking is a computational bottleneck in motion planning, requiring lazy algorithms that explicitly reason about when to perform this computation. Optimism in the face of collision uncertainty minimizes the number of checks before finding the shortest path. However, this may take a prohibitively long time to compute, with no other feasible paths discovered during this period. For many real-time applications, we instead demand strong anytime performance, defined as minimizing the cumulative lengths of the feasible paths yielded over time. We introduce Posterior Sampling for Motion Planning (PSMP), an anytime lazy motion planning algorithm that leverages learned posteriors on edge collisions to quickly discover an initial feasible path and progressively yield shorter paths. PSMP obtains an expected regret bound of $\tilde{O}(\sqrt{\mathcal{S} \mathcal{A} T})$ and outperforms comparative baselines on a set of 2D and 7D planning problems.


title: Violation of generalized fluctuation theorems in adaptively driven steady states: Applications to hair cell oscillations

abstract: The spontaneously oscillating hair bundle of sensory cells in the inner ear is an example of a stochastic, nonlinear oscillator driven by internal active processes. Moreover, this internal activity is adaptive -- its power input depends on the current state of the system. We study fluctuation dissipation relations in such adaptively-driven, nonequilibrium limit-cycle oscillators. We observe the expected violation of the well-known, equilibrium fluctuation-dissipation theorem (FDT), and verify the existence of a generalized fluctuation-dissipation theorem (GFDT) in the non-adaptively driven model of the hair cell oscillator. This generalized fluctuation theorem requires the system to be analyzed in the co-moving frame associated with the mean limit cycle of the stochastic oscillator. We then demonstrate, via numerical simulations and analytic calculations, that the adaptively-driven dynamical hair cell model violates both the FDT and the GFDT. We go on to show, using stochastic, finite-state, dynamical models, that such a feedback-controlled drive in stochastic limit cycle oscillators generically violates both the FDT and GFDT. We propose that one may in fact use the breakdown of the GFDT as a tool to more broadly look for and quantify the effect of adaptive, feedback mechanisms associated with driven (nonequilibrium) biological dynamics.


title: Computational molecular field theory for nematic liquid crystals

abstract: Nematic liquid crystals exhibit configurations in which the underlying ordering changes markedly on macroscopic length scales. Such structures include topological defects in the nematic phase and tactoids within nematic-isotropic coexistence. We discuss a computational study of inhomogeneous configurations that is based on a field theory extension of the Maier-Saupe molecular model of a uniaxial, nematic liquid crystal. A tensor order parameter is defined as the second moment of an orientational probability distribution, leading to a free energy that is not convex within the isotropic-nematic coexistence region, and that goes to infinity if the eigenvalues of the order parameter become non-physical. Computations of the spatial profile of the order parameter are presented for an isotropic-nematic interface in one dimension, a tactoid in two dimensions, and a nematic disclination in two dimensions. We compare our results to those given by the Landau de-Gennes free energy for the same configurations and discuss the advantages of such a model over the latter.


title: A canonical embedding of $\textbf{Aut}_{\textbf{hol}}({\bf \mathbb C^n})$

abstract: The group $\text{Aut}_{\text{hol}}(\mathbb C^n)$ of self-biholomorphisms of $\mathbb C^n$ consists of affine maps if $n=1$, but in higher dimensions it is a large object that has not been described explicitly. Despite the intricacies involved when $n>1$, surprisingly every $F\in \text{Aut}_{\text{hol}}(\mathbb C^n)$ is uniquely determined inside the group by only two data, of infinitesimal and global nature: the $1$-jet of $F$ at $0$, and the complex Hessian of a certain plurisubharmonic function associated to $F$. If $n=1$ this global datum is zero for all $F$, which is then determined solely by its $1$-jet at $0$, and one recovers $\text{Aut}_{\text{hol}}(\mathbb C)= \text{Aff}(\mathbb C)\cong \mathbb C \times \mathbb C^{*}$. Our main result, formulated as the existence of a canonical embedding of $ \text{Aut}_{\text{hol}} ( \mathbb C^n)$, also singles out a natural candidate for moduli space of $ \text{Aut}_{\text{hol}} ( \mathbb C^n)$, for all $n>1$.


title: HI Asymmetries in LVHIS, VIVA and HALOGAS Galaxies

abstract: We present an analysis of morphological, kinematic and spectral asymmetries in observations of atomic neutral hydrogen (HI) gas from the Local Volume HI Survey (LVHIS), the VLA Imaging of Virgo in Atomic Gas (VIVA) survey and the Hydrogen Accretion in Local Galaxies Survey (HALOGAS). With the aim of investigating the impact of the local environment density and stellar mass on the measured HI asymmetries in future large HI surveys, we provide recommendations for the most meaningful measures of asymmetry for use in future analysis. After controlling for stellar mass, we find signs of statistically significant trends of increasing asymmetries with local density. The most significant trend we measure is for the normalised flipped spectrum residual ($A_{\mathrm{spec}}$), with mean LVHIS and VIVA values of $0.204\pm0.011$ and $0.615\pm0.068$ at average weighted $10^{\mathrm{th}}$ nearest-neighbour galaxy number densities of $\log(\rho_{10}/\mathrm{Mpc}^{-3})=-1.64$ and 0.88, respectively. Looking ahead to the WALLABY survey on the Australian Square Kilometre Array Pathfinder (ASKAP), we estimate that the number of detections will be sufficient to provide coverage over 5 orders of magnitude in both local density and stellar mass increasing the dynamic range and accuracy with which we can probe the effect of these properties on the asymmetry in the distribution of atomic gas in galaxies.


title: Restriction of Eisenstein series and Stark-Heegner points

abstract: In a recent work of Darmon, Pozzi and Vonk, the authors consider a particular $p$-adic family of Hilbert Eisenstein series $E_k(1,\brch)$ associated with an odd character $\brch$ of the narrow ideal class group of a real quadratic field $F$ and compute the first derivative of a certain one-variable twisted triple product $p$-adic $L$-series attached to $E_k(1,\brch)$ and an elliptic newform $f$ of weight $2$ on $\Gamma_0(p)$. In this paper, we generalize their construction to include the cyclotomic variable and thus obtain a two-variable twisted triple product $p$-adic $L$-series. Moreover, when $f$ is associated with an elliptic curve $E$ over $\Q$, we prove that the first derivative of this $p$-adic $L$-series along the weight direction is a product of the $p$-adic logarithm of a Stark-Heegner point of $E$ over $F$ introduced by Darmon and the cyclotomic $p$-adic $L$-function for $E$.


title: Combining Laue diffraction with Bragg coherent diffraction imaging at 34-ID-C

abstract: Measurement modalities in Bragg coherent diffraction imaging (BCDI) rely on finding signal from a single nanoscale crystal object, which satisfies the Bragg condition among a large number of arbitrarily oriented nanocrystals. However, even when the signal from a single Bragg reflection with (hkl) Miller indices is found, the crystallographic axes on the retrieved three-dimensional (3D) image of the crystal remain unknown, and thus, localizing in reciprocal space other Bragg reflections becomes in reality impossible or requires good knowledge of the orientation of the crystal. We report the commissioning of a movable double-bounce Si (111) monochromator at the 34-ID-C end station of the Advanced Photon Source, which aims at delivering multi-reflection BCDI as a standard tool in a single beamline instrument. The new instrument enables this through rapid switching from monochromatic to broadband (pink) beam permitting the use of Laue diffraction to determine crystal orientation. With a proper orientation matrix determined for the lattice, one can measure coherent diffraction near multiple Bragg peaks, thus providing sufficient information to image the full strain tensor in 3D. We discuss the design, concept of operation, the developed procedures for indexing Laue patterns, and automated measuring of Bragg coherent diffraction data from multiple reflections of the same nanocrystal.


title: Stochastic Frank-Wolfe for Constrained Finite-Sum Minimization

abstract: We propose a novel Stochastic Frank-Wolfe (a.k.a. conditional gradient) algorithm for constrained smooth finite-sum minimization with a generalized linear prediction/structure. This class of problems includes empirical risk minimization with sparse, low-rank, or other structured constraints. The proposed method is simple to implement, does not require step-size tuning, and has a constant per-iteration cost that is independent of the dataset size. Furthermore, as a byproduct of the method we obtain a stochastic estimator of the Frank-Wolfe gap that can be used as a stopping criterion. Depending on the setting, the proposed method matches or improves on the best computational guarantees for Stochastic Frank-Wolfe algorithms. Benchmarks on several datasets highlight different regimes in which the proposed method exhibits a faster empirical convergence than related methods. Finally, we provide an implementation of all considered methods in an open-source package.


title: Simulation of Real-time Routing for UAS traffic Management with Communication and Airspace Safety Considerations

abstract: Small Unmanned Aircraft Systems (sUAS) will be an important component of the smart city and intelligent transportation environments of the near future. The demand for sUAS related applications, such as commercial delivery and land surveying, is expected to grow rapidly in next few years. In general, sUAS traffic routing and management functions are needed to coordinate the launching of sUAS from different launch sites and determine their trajectories to avoid conflict while considering several other constraints such as expected arrival time, minimum flight energy, and availability of communication resources. However, as the airborne sUAS density grows in a certain area, it is difficult to foresee the potential airspace and communications resource conflicts and make immediate decisions to avoid them. To address this challenge, we present a temporal and spatial routing algorithm and simulation platform for sUAS trajectory management in a high density urban area that plans sUAS movements in a spatial and temporal maze taking into account obstacles that are either static or dynamic in time. The routing allows the sUAS to avoid static no-fly areas (i.e. static obstacles) or other in-flight sUAS and areas that have congested communication resources (i.e. dynamic obstacles). The algorithm is evaluated using an agent-based simulation platform. The simulation results show that the proposed algorithm outperforms other route management algorithms in many areas, especially in processing speed and memory efficiency. Detailed comparisons are provided for the sUAS flight time, the overall throughput, conflict rate and communication resource utilization. The results demonstrate that our proposed algorithm can be used to address the airspace and communication resource utilization needs for a next generation smart city and smart transportation.


title: SWARM: Adaptive Load Balancing in Distributed Streaming Systems for Big Spatial Data

abstract: The proliferation of GPS-enabled devices has led to the development of numerous location-based services. These services need to process massive amounts of spatial data in real-time. The current scale of spatial data cannot be handled using centralized systems. This has led to the development of distributed spatial streaming systems. Existing systems are using static spatial partitioning to distribute the workload. In contrast, the real-time streamed spatial data follows non-uniform spatial distributions that are continuously changing over time. Distributed spatial streaming systems need to react to the changes in the distribution of spatial data and queries. This paper introduces SWARM, a light-weight adaptivity protocol that continuously monitors the data and query workloads across the distributed processes of the spatial data streaming system, and redistribute and rebalance the workloads soon as performance bottlenecks get detected. SWARM is able to handle multiple query-execution and data-persistence models. A distributed streaming system can directly use SWARM to adaptively rebalance the system's workload among its machines with minimal changes to the original code of the underlying spatial application. Extensive experimental evaluation using real and synthetic datasets illustrate that, on average, SWARM achieves 200% improvement over a static grid partitioning that is determined based on observing a limited history of the data and query workloads. Moreover, SWARM reduces execution latency on average 4x compared with the other technique.


title: GATCluster: Self-Supervised Gaussian-Attention Network for Image Clustering

abstract: We propose a self-supervised Gaussian ATtention network for image Clustering (GATCluster). Rather than extracting intermediate features first and then performing the traditional clustering algorithm, GATCluster directly outputs semantic cluster labels without further post-processing. Theoretically, we give a Label Feature Theorem to guarantee the learned features are one-hot encoded vectors, and the trivial solutions are avoided. To train the GATCluster in a completely unsupervised manner, we design four self-learning tasks with the constraints of transformation invariance, separability maximization, entropy analysis, and attention mapping. Specifically, the transformation invariance and separability maximization tasks learn the relationships between sample pairs. The entropy analysis task aims to avoid trivial solutions. To capture the object-oriented semantics, we design a self-supervised attention mechanism that includes a parameterized attention module and a soft-attention loss. All the guiding signals for clustering are self-generated during the training process. Moreover, we develop a two-step learning algorithm that is memory-efficient for clustering large-size images. Extensive experiments demonstrate the superiority of our proposed method in comparison with the state-of-the-art image clustering benchmarks. Our code has been made publicly available at https://github.com/niuchuangnn/GATCluster.


title: The maximal excess charge in reduced Hartree-Fock molecule

abstract: We consider a molecule described by the Hartree-Fock model without the exchange term. We prove that nuclei of total charge $Z$ can bind at most $Z+C$ electrons, where $C$ is a constant independent of $Z$.


title: On the extension property of dilatation monotone risk measures

abstract: Let $\mathcal{X}$ be a subset of $L^1$ that contains the space of simple random variables $\mathcal{L}$ and $\rho: \mathcal{X} \rightarrow (-\infty,\infty]$ a dilatation monotone functional with the Fatou property. In this note, we show that $\rho$ extends uniquely to a $\sigma(L^1,\mathcal{L})$ lower semicontinuous and dilatation monotone functional $\overline{\rho}: L^1 \rightarrow (-\infty,\infty]$. Moreover, $\overline{\rho}$ preserves monotonicity, (quasi)convexity, and cash-additivity of $\rho$. Our findings complement recent extension results for quasiconvex law-invariant functionals proved in [17,20]. As an application of our results, we show that transformed norm risk measures on Orlicz hearts admit a natural extension to $L^1$ that retains the robust representations obtained in [4,6].


title: Structure and morphology of low mechanical loss TiO$_2$-doped Ta$_2$O$_5$

abstract: Amorphous oxide thin films play a fundamental role in state-of-the art interferometry experiments, such as gravitational wave detectors where these films compose the high reflectance mirrors of end and input masses. The sensitivity of these detectors is affected by thermal noise in the mirrors with its main source being the mechanical loss of the high index layers. These thermally driven fluctuations are a fundamental limit to optical interferometry experiments and there is a pressing need to understand the underlying processes that lead to mechanical dissipation in materials at room temperature. Two strategies are known to lower the mechanical loss: employing a mixture of Ta$_2$O$_5$ with $\approx$ 20\% of TiO$_2$ and post-deposition annealing, but the reasons behind this are not completely understood. In this work, we present a systematic study of the structural and optical properties of ion beam sputtered TiO$_2$-doped Ta$_2$O$_5$ films as a function of the annealing temperature. We show for the first time that low mechanical loss is associated with a material morphology that consists of nanometer sized Ar-rich bubbles embedded into an atomically homogeneous mixed titanium-tantalum oxide. When the Ti cation ratio is high, however, phase separation occurs in the film which leads to increased mechanical loss. These results indicate that for designing low mechanical loss mixed oxide coatings for interferometry applications it would be beneficial to identify materials with the ability to form ternary compounds while the dopant ratio needs to be kept low to avoid phase separation.


title: Bridging the Gap between Spatial and Spectral Domains: A Survey on Graph Neural Networks

abstract: Deep learning's success has been widely recognized in a variety of machine learning tasks, including image classification, audio recognition, and natural language processing. As an extension of deep learning beyond these domains, graph neural networks (GNNs) are designed to handle the non-Euclidean graph-structure which is intractable to previous deep learning techniques. Existing GNNs are presented using various techniques, making direct comparison and cross-reference more complex. Although existing studies categorize GNNs into spatial-based and spectral-based techniques, there hasn't been a thorough examination of their relationship. To close this gap, this study presents a single framework that systematically incorporates most GNNs. We organize existing GNNs into spatial and spectral domains, as well as expose the connections within each domain. A review of spectral graph theory and approximation theory builds a strong relationship across the spatial and spectral domains in further investigation.


title: Grothendieck meeting [Wess & Bagger]: [Supersymmetry and supergravity: IV, V, VI, VII, XXII] (2nd ed.) reconstructed in complexified ${\Bbb Z}/2$-graded $C^\infty$-Algebraic Geometry, I. Construction under trivialization of spinor bundle

abstract: Forty-six years after the birth of supersymmetry in 1973 from works of Julius Wess and Bruno Zumino, the standard quantum-field-theorists and particle physicists' language of `superspaces', `supersymmetry', and `supersymmetric action functionals in superspace formulation' as given in Chapters IV, V, VI, VII, XXII of the classic on supersymmetry and supergravity: Julius Wess & Jonathan Bagger: Supersymmetry and Supergravity (2nd ed.), is finally polished, with only minimal mathematical patches added for consistency and accuracy in dealing with nilpotent objects from the Grassmann algebra involved, to a precise setting in the language of complexified ${\Bbb Z}/2$-graded $C^\infty$-Algebraic Geometry. This is completed after the lesson learned from D(14.1) (arXiv:1808.05011 [math.DG]) and the notion of `$d=3+1$, $N=1$ towered superspaces' as complexified ${\Bbb Z}/2$-graded $C^\infty$-schemes, their distinguished sectors, and purge-evaluation maps first developed in SUSY(1) (= D(14.1.Supp.1)) (arXiv:1902.06246 [hep-th]) and further polished in the current work. While the construction depends on a choice of a trivialization of the spinor bundle by covariantly constant sections, as long as the transformation law and the induced isomorphism under a change of trivialization of the spinor bundle by covariantly constant sections are understood, any object or structure thus defined or constructed is mathematically well-defined. The construction can be generalized to all other space-time dimensions with simple or extended supersymmetries. This is part of the mathematical foundation required to study fermionic D-branes in the Ramond-Neveu-Schwarz formulation.


title: Controllable Level Blending between Games using Variational Autoencoders

abstract: Previous work explored blending levels from existing games to create levels for a new game that mixes properties of the original games. In this paper, we use Variational Autoencoders (VAEs) for improving upon such techniques. VAEs are artificial neural networks that learn and use latent representations of datasets to generate novel outputs. We train a VAE on level data from Super Mario Bros. and Kid Icarus, enabling it to capture the latent space spanning both games. We then use this space to generate level segments that combine properties of levels from both games. Moreover, by applying evolutionary search in the latent space, we evolve level segments satisfying specific constraints. We argue that these affordances make the VAE-based approach especially suitable for co-creative level design and compare its performance with similar generative models like the GAN and the VAE-GAN.


title: A Structure-based Memory Maintenance Model for Neural Tracking of Linguistic Structures

abstract: It is recently demonstrated that cortical activity can track the time courses of phrases and sentences during speech listening. Here, we propose a plausible neural processing framework to explain this phenomenon. It is argued that the brain maintains the neural representation of a linguistic unit, i.e., a word or a phrase, in a processing buffer until the unit is integrated into a higher-level structure. After being integrated, the unit is removed from the buffer and becomes activated long-term memory. In this model, the duration each unit is maintained in the processing buffer depends on the linguistic structure of the speech input. It is shown that the number of items retained in the processing buffer follows the time courses of phrases and sentences, in line with neurophysiological data, whether the syntactic structure of a sentence is mentally parsed using a bottom-up or top-down predictive model. This model generates a range of testable predictions about the link between linguistic structures, their dynamic psychological representations and their neural underpinnings.


title: Direct Multipixel Imaging and Spectroscopy of an Exoplanet with a Solar Gravity Lens Mission

abstract: We examined the solar gravitational lens (SGL) as the means to produce direct high-resolution, multipixel images of exoplanets. The properties of the SGL are remarkable: it offers maximum light amplification of ~1e11 and angular resolution of ~1e-10 arcsec. A probe with a 1-m telescope in the SGL focal region can image an exoplanet at 30 pc with 10-kilometer resolution on its surface, sufficient to observe seasonal changes, oceans, continents, surface topography. We reached and exceeded all objectives set for our study: We developed a new wave-optical approach to study the imaging of exoplanets while treating them as extended, resolved, faint sources at large but finite distances. We properly accounted for the solar corona brightness. We developed deconvolution algorithms and demonstrated the feasibility of high-quality image reconstruction under realistic conditions. We have proven that multipixel imaging and spectroscopy of exoplanets with the SGL are feasible. We have developed a new mission concept that delivers an array of optical telescopes to the SGL focal region relying on three innovations: i) a new way to enable direct exoplanet imaging, ii) use of smallsats solar sails fast transit through the solar system and beyond, iii) an open architecture to take advantage of swarm technology. This approach enables entirely new missions, providing a great leap in capabilities for NASA and the greater aerospace community. Our results are encouraging as they lead to a realistic design for a mission that will be able to make direct resolved images of exoplanets in our stellar neighborhood. It could allow exploration of exoplanets relying on the SGL capabilities decades, if not centuries, earlier than possible with other extant technologies. The architecture and mission concepts for a mission to the strong interference region of the SGL are promising and should be explored further.

categories: astro-ph.IM astro-ph.EP gr-qc physics.optics


title: Purely singular splittings of cyclic groups

abstract: Let $G$ be a finite abelian group. We say that $M$ and $S$ form a \textsl{splitting} of $G$ if every nonzero element $g$ of $G$ has a unique representation of the form $g=ms$ with $m\in M$ and $s\in S$, while $0$ has no such representation. The splitting is called \textit{purely singular} if for each prime divisor $p$ of $|G|$, there is at least one element of $M$ is divisible by $p$. In this paper, we mainly study the purely singular splittings of cyclic groups. We first prove that if $k\ge3$ is a positive integer such that $[-k+1, \,k]^*$ splits a cyclic group $\mathbb{Z}_m$, then $m=2k$. Next, we have the following general result. Suppose $M=[-k_1, \,k_2]^*$ splits $\mathbb{Z}_{n(k_1+k_2)+1}$ with $1\leq k_1< k_2$. If $n\geq 2$, then $k_1\leq n-2$ and $k_2\leq 2n-5$. Applying this result, we prove that if $M=[-k_1, \,k_2]^*$ splits $\mathbb{Z}_m$ purely singularly, and either $(i)$ $\gcd(s, \,m)=1$ for all $s\in S$ or $(ii)$ $m=2^{\alpha}p^{\beta}$ or $2^{\alpha}p_1p_2$ with $\alpha\geq 0$, $\beta\geq 1$ and $p$, $p_1$, $p_2$ odd primes, then $m=k_1+k_2+1$ or $k_1=0$ and $m=k_2+1$ or $2k_2+1$.


title: Evidence for Returning Disk Radiation in the Black Hole X-ray Binary XTE J1550-564

abstract: We explore the accretion properties of the black hole X-ray binary \j1550\ during its outbursts in 1998/99 and 2000. We model the disk, corona, and reflection components of X-ray spectra taken with the {\it Rossi X-ray Timing Explorer} (\rxte), using the {\tt relxill} suite of reflection models. The key result of our modeling is that the reflection spectrum in the very soft state is best explained by disk self-irradiation, i.e., photons from the inner disk are bent by the strong gravity of the black hole, and reflected off the disk surface. This is the first known detection of thermal disk radiation reflecting off the inner disk. There is also an apparent absorption line at $\sim6.9$ keV which may be evidence of an ionized disk wind. The coronal electron temperature ($kT_{\rm e}$) is, as expected, lower in the brighter outburst of 1998/99, explained qualitatively by more efficient coronal cooling due to irradiating disk photons. The disk inner radius is consistent with being within a few times the innermost stable circular orbit (ISCO) throughout the bright-hard-to-soft states (10s of $r_{\rm g}$ in gravitational units). The disk inclination is low during the hard state, disagreeing with the binary inclination value, and very close to $90^{\circ}$ in the soft state, recovering to a lower value when adopting a blackbody spectrum as the irradiating continuum.


title: Learning Scalable Multi-Agent Coordination by Spatial Differentiation for Traffic Signal Control

abstract: The intelligent control of the traffic signal is critical to the optimization of transportation systems. To achieve global optimal traffic efficiency in large-scale road networks, recent works have focused on coordination among intersections, which have shown promising results. However, existing studies paid more attention to observations sharing among intersections (both explicit and implicit) and did not care about the consequences after decisions. In this paper, we design a multiagent coordination framework based on Deep Reinforcement Learning methods for traffic signal control, defined as {\gamma}-Reward that includes both original {\gamma}-Reward and {\gamma}-Attention-Reward. Specifically, we propose the Spatial Differentiation method for coordination which uses the temporal-spatial information in the replay buffer to amend the reward of each action. A concise theoretical analysis that proves the proposed model can converge to Nash equilibrium is given. By extending the idea of Markov Chain to the dimension of space-time, this truly decentralized coordination mechanism replaces the graph attention method and realizes the decoupling of the road network, which is more scalable and more in line with practice. The simulation results show that the proposed model remains a state-of-the-art performance even not use a centralized setting. Code is available in https://github.com/Skylark0924/Gamma Reward.


title: Optimality and Stability in Non-Convex Smooth Games

abstract: Convergence to a saddle point for convex-concave functions has been studied for decades, while recent years has seen a surge of interest in non-convex (zero-sum) smooth games, motivated by their recent wide applications. It remains an intriguing research challenge how local optimal points are defined and which algorithm can converge to such points. An interesting concept is known as the local minimax point, which strongly correlates with the widely-known gradient descent ascent algorithm. This paper aims to provide a comprehensive analysis of local minimax points, such as their relation with other solution concepts and their optimality conditions. We find that local saddle points can be regarded as a special type of local minimax points, called uniformly local minimax points, under mild continuity assumptions. In (non-convex) quadratic games, we show that local minimax points are (in some sense) equivalent to global minimax points. Finally, we study the stability of gradient algorithms near local minimax points. Although gradient algorithms can converge to local/global minimax points in the non-degenerate case, they would often fail in general cases. This implies the necessity of either novel algorithms or concepts beyond saddle points and minimax points in non-convex smooth games.


title: Quantitative estimate of the continuum approximations of interacting particle systems in one dimension

abstract: We consider a large class of interacting particle systems in 1D described by an energy whose interaction potential is singular and non-local. This class covers Riesz gases (in particular, log gases) and applications to plasticity and approximation theory of functions. While it is well established that the minimisers of such interaction energies converge to a certain particle density profile as the number of particles tends to infinity, any bound on the rate of this convergence is only known in special cases by means of quantitative estimates. The main result of this paper extends these quantitative estimates to a large class of interaction energies by a different proof. The proof relies on one-dimensional features such as the convexity of the interaction potential and the ordering of the particles. The main novelty of the proof is the treatment of the singularity of the interaction potential by means of a carefully chosen renormalisation.


title: Anapole mediated giant photothermal nonlinearity in nanostructured silicon

abstract: Featured with a plethora of electric and magnetic Mie resonances, high index dielectric nanostructures offer a versatile platform to concentrate light-matter interactions at the nanoscale. By integrating unique features of far-field scattering control and near-field concentration from radiationless anapole states, here, we demonstrate a giant photothermal nonlinearity in single subwavelength-sized silicon nanodisks. The nanoscale energy concentration and consequent near-field enhancements mediated by the anapole mode yield a reversible nonlinear scattering with a large modulation depth and a broad dynamic range, unveiling a record-high nonlinear index change up to 0.5 at mild incident light intensities on the order of MW/cm2. The observed photothermal nonlinearity showcases three orders of magnitude enhancement compared with that of unstructured bulk silicon, as well as nearly one order of magnitude higher than that through the radiative electric dipolar mode. Such nonlinear scattering can empower distinctive point spread functions in confocal reflectance imaging, offering the potential for far-field localization of nanostructured Si with an accuracy approaching 40 nm. Our findings shed new light on active silicon photonics based on optical anapoles.


title: Subsets of rectifiable curves in Banach spaces I: sharp exponents in traveling salesman theorems

abstract: The Analyst's Traveling Salesman Problem is to find a characterization of subsets of rectifiable curves in a metric space. This problem was introduced and solved in the plane by Jones in 1990 and subsequently solved in higher-dimensional Euclidean spaces by Okikiolu in 1992 and in the infinite-dimensional Hilbert space $\ell_2$ by Schul in 2007. In this paper, we establish sharp extensions of Schul's necessary and sufficient conditions for a bounded set $E\subset\ell_p$ to be contained in a rectifiable curve from $p=2$ to $1<p<\infty$. While the necessary and sufficient conditions coincide when $p=2$, we demonstrate that there is a strict gap between the necessary condition and sufficient condition when $p\neq 2$. We also identify and correct technical errors in the proof by Schul. This investigation is partly motivated by recent work of Edelen, Naber, and Valtorta on Reifenberg-type theorems in Banach spaces and complements work of Hahlomaa and recent work of David and Schul on the Analyst's TSP in general metric spaces.


title: State-only Imitation with Transition Dynamics Mismatch

abstract: Imitation Learning (IL) is a popular paradigm for training agents to achieve complicated goals by leveraging expert behavior, rather than dealing with the hardships of designing a correct reward function. With the environment modeled as a Markov Decision Process (MDP), most of the existing IL algorithms are contingent on the availability of expert demonstrations in the same MDP as the one in which a new imitator policy is to be learned. This is uncharacteristic of many real-life scenarios where discrepancies between the expert and the imitator MDPs are common, especially in the transition dynamics function. Furthermore, obtaining expert actions may be costly or infeasible, making the recent trend towards state-only IL (where expert demonstrations constitute only states or observations) ever so promising. Building on recent adversarial imitation approaches that are motivated by the idea of divergence minimization, we present a new state-only IL algorithm in this paper. It divides the overall optimization objective into two subproblems by introducing an indirection step and solves the subproblems iteratively. We show that our algorithm is particularly effective when there is a transition dynamics mismatch between the expert and imitator MDPs, while the baseline IL methods suffer from performance degradation. To analyze this, we construct several interesting MDPs by modifying the configuration parameters for the MuJoCo locomotion tasks from OpenAI Gym.


title: Stochastic Matching with Few Queries: $(1-\varepsilon)$ Approximation

abstract: Suppose that we are given an arbitrary graph $G=(V, E)$ and know that each edge in $E$ is going to be realized independently with some probability $p$. The goal in the stochastic matching problem is to pick a sparse subgraph $Q$ of $G$ such that the realized edges in $Q$, in expectation, include a matching that is approximately as large as the maximum matching among the realized edges of $G$. The maximum degree of $Q$ can depend on $p$, but not on the size of $G$. This problem has been subject to extensive studies over the years and the approximation factor has been improved from $0.5$ to $0.5001$ to $0.6568$ and eventually to $2/3$. In this work, we analyze a natural sampling-based algorithm and show that it can obtain all the way up to $(1-\epsilon)$ approximation, for any constant $\epsilon > 0$. A key and of possible independent interest component of our analysis is an algorithm that constructs a matching on a stochastic graph, which among some other important properties, guarantees that each vertex is matched independently from the vertices that are sufficiently far. This allows us to bypass a previously known barrier towards achieving $(1-\epsilon)$ approximation based on existence of dense Ruzsa-Szemer\'edi graphs.


title: Defense-PointNet: Protecting PointNet Against Adversarial Attacks

abstract: Despite remarkable performance across a broad range of tasks, neural networks have been shown to be vulnerable to adversarial attacks. Many works focus on adversarial attacks and defenses on 2D images, but few focus on 3D point clouds. In this paper, our goal is to enhance the adversarial robustness of PointNet, which is one of the most widely used models for 3D point clouds. We apply the fast gradient sign attack method (FGSM) on 3D point clouds and find that FGSM can be used to generate not only adversarial images but also adversarial point clouds. To minimize the vulnerability of PointNet to adversarial attacks, we propose Defense-PointNet. We compare our model with two baseline approaches and show that Defense-PointNet significantly improves the robustness of the network against adversarial samples.


title: A Visual Communication Map for Multi-Agent Deep Reinforcement Learning

abstract: Deep reinforcement learning has been applied successfully to solve various real-world problems and the number of its applications in the multi-agent settings has been increasing. Multi-agent learning distinctly poses significant challenges in the effort to allocate a concealed communication medium. Agents receive thorough knowledge from the medium to determine subsequent actions in a distributed nature. Apparently, the goal is to leverage the cooperation of multiple agents to achieve a designated objective efficiently. Recent studies typically combine a specialized neural network with reinforcement learning to enable communication between agents. This approach, however, limits the number of agents or necessitates the homogeneity of the system. In this paper, we have proposed a more scalable approach that not only deals with a great number of agents but also enables collaboration between dissimilar functional agents and compatibly combined with any deep reinforcement learning methods. Specifically, we create a global communication map to represent the status of each agent in the system visually. The visual map and the environmental state are fed to a shared-parameter network to train multiple agents concurrently. Finally, we select the Asynchronous Advantage Actor-Critic (A3C) algorithm to demonstrate our proposed scheme, namely Visual communication map for Multi-agent A3C (VMA3C). Simulation results show that the use of visual communication map improves the performance of A3C regarding learning speed, reward achievement, and robustness in multi-agent problems.


title: Review, Analysis and Design of a Comprehensive Deep Reinforcement Learning Framework

abstract: The integration of deep learning to reinforcement learning (RL) has enabled RL to perform efficiently in high-dimensional environments. Deep RL methods have been applied to solve many complex real-world problems in recent years. However, development of a deep RL-based system is challenging because of various issues such as the selection of a suitable deep RL algorithm, its network configuration, training time, training methods, and so on. This paper proposes a comprehensive software framework that not only plays a vital role in designing a connect-the-dots deep RL architecture but also provides a guideline to develop a realistic RL application in a short time span. We have designed and developed a deep RL-based software framework that strictly ensures flexibility, robustness, and scalability. By inheriting the proposed architecture, software managers can foresee any challenges when designing a deep RL-based system. As a result, they can expedite the design process and actively control every stage of software development, which is especially critical in agile development environments. To enforce generalization, the proposed architecture does not depend on a specific RL algorithm, a network configuration, the number of agents, or the type of agents. Using our framework, software developers can develop and integrate new RL algorithms or new types of agents, and can flexibly change network configuration or the number of agents.


title: Tighter uncertainty relations based on Wigner-Yanase skew information for observables and channels

abstract: Uncertainty principle is the basis of quantum mechanics. It reflects the basic law of the movement of microscopic particles. Wigner-Yanase skew information, as a measure of quantum uncertainties, is used to characterize the intrinsic features of the state and the observable. In this paper, we mainly investigate the sum uncertainty relations for both quantum mechanical observables and quantum channels based on skew information. We establish a new uncertainty relation in terms of Wigner-Yanase skew information for $n$ observables, which is saturated (thus it holds as equality) for two incompatible observables. We also present two uncertainty relations for arbitrary finite $N$ quantum channels by using skew information. Our uncertainty relations have tighter lower bounds than the existing ones. Detailed examples are provided.


title: Kernel Bi-Linear Modeling for Reconstructing Data on Manifolds: The Dynamic-MRI Case

abstract: This paper establishes a kernel-based framework for reconstructing data on manifolds, tailored to fit the dynamic-(d)MRI-data recovery problem. The proposed methodology exploits simple tangent-space geometries of manifolds in reproducing kernel Hilbert spaces and follows classical kernel-approximation arguments to form the data-recovery task as a bi-linear inverse problem. Departing from mainstream approaches, the proposed methodology uses no training data, employs no graph Laplacian matrix to penalize the optimization task, uses no costly (kernel) pre-imaging step to map feature points back to the input space, and utilizes complex-valued kernel functions to account for k-space data. The framework is validated on synthetically generated dMRI data, where comparisons against state-of-the-art schemes highlight the rich potential of the proposed approach in data-recovery problems.


title: Hierarchical Memory Decoding for Video Captioning

abstract: Recent advances of video captioning often employ a recurrent neural network (RNN) as the decoder. However, RNN is prone to diluting long-term information. Recent works have demonstrated memory network (MemNet) has the advantage of storing long-term information. However, as the decoder, it has not been well exploited for video captioning. The reason partially comes from the difficulty of sequence decoding with MemNet. Instead of the common practice, i.e., sequence decoding with RNN, in this paper, we devise a novel memory decoder for video captioning. Concretely, after obtaining representation of each frame through a pre-trained network, we first fuse the visual and lexical information. Then, at each time step, we construct a multi-layer MemNet-based decoder, i.e., in each layer, we employ a memory set to store previous information and an attention mechanism to select the information related to the current input. Thus, this decoder avoids the dilution of long-term information. And the multi-layer architecture is helpful for capturing dependencies between frames and word sequences. Experimental results show that even without the encoding network, our decoder still could obtain competitive performance and outperform the performance of RNN decoder. Furthermore, compared with one-layer RNN decoder, our decoder has fewer parameters.


title: Using a thousand optimization tasks to learn hyperparameter search strategies

abstract: We present TaskSet, a dataset of tasks for use in training and evaluating optimizers. TaskSet is unique in its size and diversity, containing over a thousand tasks ranging from image classification with fully connected or convolutional neural networks, to variational autoencoders, to non-volume preserving flows on a variety of datasets. As an example application of such a dataset we explore meta-learning an ordered list of hyperparameters to try sequentially. By learning this hyperparameter list from data generated using TaskSet we achieve large speedups in sample efficiency over random search. Next we use the diversity of the TaskSet and our method for learning hyperparameter lists to empirically explore the generalization of these lists to new optimization tasks in a variety of settings including ImageNet classification with Resnet50 and LM1B language modeling with transformers. As part of this work we have opensourced code for all tasks, as well as ~29 million training curves for these problems and the corresponding hyperparameters.


title: Magnetic effects on the solvability of 2D MHD boundary layer equations without resistivity in Sobolev spaces

abstract: In this paper, we are concerned with the magnetic effect on the Sobolev solvability of boundary layer equations for the 2D incompressible MHD system without resistivity. The MHD boundary layer is described by the Prandtl type equations derived from the incompressible viscous MHD system without resistivity under the no-slip boundary condition on the velocity. Assuming that the initial tangential magnetic field does not degenerate, a local-in-time well-posedness in Sobolev spaces is proved without the monotonicity condition on the velocity field. Moreover, we show that if the tangential magnetic field shear layer is degenerate at one point, then the linearized MHD boundary layer system around the shear layer profile is ill-posed in the Sobolev settings provided that the initial velocity shear flow is non-degenerately critical at the same point.


title: Long-range collectivity in small collision-systems with two- and four-particle correlations @ STAR

abstract: New STAR differential and integral v{2,3} measurements that explicitly account for non-flow contributions are reported for p/d/3He+Au, collisions at Roots=200 GeV. The measurements, which leverage the two-particle correlators for p/d/3He+Au and minimum-bias p+p collisions in tandem with three well-established methods of non-flow subtraction, are observed to be method-independent. For comparable multiplicities, they further indicate system-independent v2{2} and v3{2} values that are consistent with the critical role of both size (N{ch}) and the subnucleonic-fluctuations-driven eccentricities e{2,3} but are inconsistent with the notion of shape engineering in p/d/3He+Au collisions. The scaling function derived from the measurements, confirm the important role of final-state effects across a broad spectrum of collision-system sizes and energies and suggests an increase in eta/s(T,muB) for small collision-systems.


title: HAM: Hybrid Associations Models for Sequential Recommendation

abstract: Sequential recommendation aims to identify and recommend the next few items for a user that the user is most likely to purchase/review, given the user's purchase/rating trajectories. It becomes an effective tool to help users select favorite items from a variety of options. In this manuscript, we developed hybrid associations models (HAM) to generate sequential recommendations using three factors: 1) users' long-term preferences, 2) sequential, high-order and low-order association patterns in the users' most recent purchases/ratings, and 3) synergies among those items. HAM uses simplistic pooling to represent a set of items in the associations, and element-wise product to represent item synergies of arbitrary orders. We compared HAM models with the most recent, state-of-the-art methods on six public benchmark datasets in three different experimental settings. Our experimental results demonstrate that HAM models significantly outperform the state of the art in all the experimental settings, with an improvement as much as 46.6%. In addition, our run-time performance comparison in testing demonstrates that HAM models are much more efficient than the state-of-the-art methods, and are able to achieve significant speedup as much as 139.7 folds.


title: BBAND Index: A No-Reference Banding Artifact Predictor

abstract: Banding artifact, or false contouring, is a common video compression impairment that tends to appear on large flat regions in encoded videos. These staircase-shaped color bands can be very noticeable in high-definition videos. Here we study this artifact, and propose a new distortion-specific no-reference video quality model for predicting banding artifacts, called the Blind BANding Detector (BBAND index). BBAND is inspired by human visual models. The proposed detector can generate a pixel-wise banding visibility map and output a banding severity score at both the frame and video levels. Experimental results show that our proposed method outperforms state-of-the-art banding detection algorithms and delivers better consistency with subjective evaluations.


title: Multi-Robot Path Planning Using Medial-Axis-Based Pebble-Graph Embedding

abstract: We present a centralized algorithm for labeled, disk-shaped Multi-Robot Path Planning (MPP) in a continuous planar workspace with polygonal boundaries. Our method automatically transform the continuous problem into a discrete, graph-based variant termed the pebble motion problem, which can be solved efficiently. To construct the underlying pebble graph, we identify inscribed circles in the workspace via a medial axis transform and organize robots into layers within each inscribed circle. We show that our layered pebble-graph enables collision-free motions, allowing all graph-restricted MPP instances to be feasible. MPP instances with continuous start and goal positions can then be solved via local navigations that route robots from and to graph vertices. We tested our method on several environments with high robot-packing densities (up to $61.6\%$ of the workspace). For environments with narrow passages, such density violates the well-separated assumptions made by state-of-the-art MPP planners, while our method achieves an average success rate of $83\%$.


title: CrossWOZ: A Large-Scale Chinese Cross-Domain Task-Oriented Dialogue Dataset

abstract: To advance multi-domain (cross-domain) dialogue modeling as well as alleviate the shortage of Chinese task-oriented datasets, we propose CrossWOZ, the first large-scale Chinese Cross-Domain Wizard-of-Oz task-oriented dataset. It contains 6K dialogue sessions and 102K utterances for 5 domains, including hotel, restaurant, attraction, metro, and taxi. Moreover, the corpus contains rich annotation of dialogue states and dialogue acts at both user and system sides. About 60% of the dialogues have cross-domain user goals that favor inter-domain dependency and encourage natural transition across domains in conversation. We also provide a user simulator and several benchmark models for pipelined task-oriented dialogue systems, which will facilitate researchers to compare and evaluate their models on this corpus. The large size and rich annotation of CrossWOZ make it suitable to investigate a variety of tasks in cross-domain dialogue modeling, such as dialogue state tracking, policy learning, user simulation, etc.


title: Unshuffling Data for Improved Generalization

abstract: Generalization beyond the training distribution is a core challenge in machine learning. The common practice of mixing and shuffling examples when training neural networks may not be optimal in this regard. We show that partitioning the data into well-chosen, non-i.i.d. subsets treated as multiple training environments can guide the learning of models with better out-of-distribution generalization. We describe a training procedure to capture the patterns that are stable across environments while discarding spurious ones. The method makes a step beyond correlation-based learning: the choice of the partitioning allows injecting information about the task that cannot be otherwise recovered from the joint distribution of the training data. We demonstrate multiple use cases with the task of visual question answering, which is notorious for dataset biases. We obtain significant improvements on VQA-CP, using environments built from prior knowledge, existing meta data, or unsupervised clustering. We also get improvements on GQA using annotations of "equivalent questions", and on multi-dataset training (VQA v2 / Visual Genome) by treating them as distinct environments.


title: Proceedings 8th International Workshop on Theorem Proving Components for Educational Software

abstract: This EPTCS volume contains the proceedings of the ThEdu'19 workshop, promoted on August 25, 2019, as a satellite event of CADE-27, in Natal, Brazil. Representing the eighth installment of the ThEdu series, ThEdu'19 was a vibrant workshop, with an invited talk by Sarah Winkler, four contributions, and the first edition of a Geometry Automated Provers Competition. After the workshop an open call for papers was issued and attracted seven submissions, six of which have been accepted by the reviewers, and collected in the present post-proceedings volume. The ThEdu series pursues the smooth transition from an intuitive way of doing mathematics at secondary school to a more formal approach to the subject in STEM education, while favoring software support for this transition by exploiting the power of theorem-proving technologies. The volume editors hope that this collection of papers will further promote the development of theorem-proving-based software, and that it will collaborate on improving mutual understanding between computer mathematicians and stakeholders in education.


title: Gradient Boosted Normalizing Flows

abstract: By chaining a sequence of differentiable invertible transformations, normalizing flows (NF) provide an expressive method of posterior approximation, exact density evaluation, and sampling. The trend in normalizing flow literature has been to devise deeper, more complex transformations to achieve greater flexibility. We propose an alternative: Gradient Boosted Normalizing Flows (GBNF) model a density by successively adding new NF components with gradient boosting. Under the boosting framework, each new NF component optimizes a sample weighted likelihood objective, resulting in new components that are fit to the residuals of the previously trained components. The GBNF formulation results in a mixture model structure, whose flexibility increases as more components are added. Moreover, GBNFs offer a wider, as opposed to strictly deeper, approach that improves existing NFs at the cost of additional training---not more complex transformations. We demonstrate the effectiveness of this technique for density estimation and, by coupling GBNF with a variational autoencoder, generative modeling of images. Our results show that GBNFs outperform their non-boosted analog, and, in some cases, produce better results with smaller, simpler flows.


title: Missing water in Class I protostellar disks

abstract: Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determine the water abundance in embedded Keplerian disks through spatially-resolved observations of H$_2^{18}$O lines to understand the evolution of water during star and planet formation. We present H$_2^{18}$O line observations with ALMA and NOEMA millimeter interferometers toward five young stellar objects. NOEMA observed the 3$_{1,3}$ - $2_{2,0}$ line (E$_{\rm up}$ = 203.7 K) while ALMA targeted the $4_{1,4}$ - $3_{2,1}$ line (E$_{\rm up}$ = 322.0 K). Water column densities are derived considering optically thin and thermalized emission. Our observations are sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample. No H$_2^{18}$O emission is detected toward any of our five Class I disks. We report upper limits to the integrated line intensities. The inferred water column densities in Class I disks are N < 10$^{15}$ cm$^{-2}$ on 100 au scales which include both disk and envelope. The upper limits imply a disk-averaged water abundance of $\lesssim 10^{-6}$ with respect to H$_2$ for Class I objects. After taking into account the physical structure of the disk, the upper limit to the water abundance averaged over the inner warm disk with $T>$ 100 K is between 10$^{-7}$ up to 10$^{-5}$. Water vapor is not abundant in warm protostellar envelopes around Class I protostars. Upper limits to the water vapor column densities in Class I disks are at least two orders magnitude lower than values found in Class 0 disk-like structures.


title: A Neuromorphic Proto-Object Based Dynamic Visual Saliency Model with an FPGA Implementation

abstract: The ability to attend to salient regions of a visual scene is an innate and necessary preprocessing step for both biological and engineered systems performing high-level visual tasks (e.g. object detection, tracking, and classification). Computational efficiency, in regard to processing bandwidth and speed, is improved by only devoting computational resources to salient regions of the visual stimuli. In this paper, we first present a neuromorphic, bottom-up, dynamic visual saliency model based on the notion of proto-objects. This is achieved by incorporating the temporal characteristics of the visual stimulus into the model, similarly to the manner in which early stages of the human visual system extracts temporal information. This neuromorphic model outperforms state-of-the-art dynamic visual saliency models in predicting human eye fixations on a commonly used video dataset with associated eye tracking data. Secondly, for this model to have practical applications, it must be capable of performing its computations in real-time under low-power, small-size, and lightweight constraints. To address this, we introduce a Field-Programmable Gate Array implementation of the model on an Opal Kelly 7350 Kintex-7 board. This novel hardware implementation allows for processing of up to 23.35 frames per second running on a 100 MHz clock - better than 26x speedup from the software implementation.


title: Moments and Non-vanishing of central values of Quadratic Hecke $L$-functions in the Gaussian Field

abstract: We evaluate the first three moments of central values of a family of qudratic Hecke $L$-functions in the Gaussian field with power saving error terms. In particular, we obtain asymptotic formulas for the first two moments with error terms of size $O(X^{1/2+\varepsilon})$. We also study the first and second mollified moments of the same family of $L$-functions to show that at least $87.5\%$ of the members of this family have non-vanishing central values.


title: Scaling exponents saturate in three-dimensional isotropic turbulence

abstract: From a database of direct numerical simulations of homogeneous and isotropic turbulence, generated in periodic boxes of various sizes, we extract the spherically symmetric part of moments of velocity increments and first verify the following (somewhat contested) results: the $4/5$-ths law holds in an intermediate range of scales and that the second order exponent over the same range of scales is {\it{anomalous}}, departing from the self-similar value of $2/3$ and approaching a constant of $0.72$ at high Reynolds numbers. We compare with some typical theories the dependence of longitudinal exponents as well as their derivatives with respect to the moment order $n$, and estimate the most probable value of the H\"older exponent. We demonstrate that the transverse scaling exponents saturate for large $n$, and trace this trend to the presence of large localized jumps in the signal. The saturation value of about $2$ at the highest Reynolds number suggests, when interpreted in the spirit of fractals, the presence of vortex sheets rather than more complex singularities. In general, the scaling concept in hydrodynamic turbulence appears to be more complex than even the multifractal description.

categories: physics.flu-dyn physics.comp-ph physics.data-an stat.CO


title: A Purely-Reactive Manipulability-Maximising Motion Controller

abstract: We present a novel approach to controlling the instantaneous velocity of a robot end-effector that is able to simultaneously maximise manipulability and avoid joint limits. It operates on non-redundant and redundant robots, which is achieved by adding artificial redundancy in the form of controlled path deviation. We formulate the problem as a quadratic programme and provide an open-source Python implementation that provides solutions in just a few milliseconds. It accepts a robot model expressed using URDF or Denavit-Hartenberg parameterisation. We compare our method to previous work in simulation and on a physical robot.


title: Relations among $k$-ME concurrence, negativity, polynomial invariants, and tangle

abstract: The $k$-ME concurrence as a measure of multipartite entanglement (ME) unambiguously detects all $k$-nonseparable states in arbitrary dimensions, and satisfies many important properties of an entanglement measure. Negativity is a simple computable bipartite entanglement measure. Invariant and tangle are useful tools to study the properties of the quantum states. In this paper we mainly investigate the internal relations among the $k$-ME concurrence, negativity, polynomial invariants, and tangle. Strong links between $k$-ME concurrence and negativity as well as between $k$-ME concurrence and polynomial invariants are derived. We obtain the quantitative relation between $k$-ME ($k$=$n$) concurrence and negativity for all $n$-qubit states, give a exact value of the $n$-ME concurrence for the mixture of $n$-qubit GHZ states and white noise, and derive an connection between $k$-ME concurrence and tangle for $n$-qubit W state. Moreover, we find that for any $3$-qubit pure state the $k$-ME concurrence ($k$=2, 3) is related to negativity, tangle and polynomial invariants, while for $4$-qubit states the relations between $k$-ME concurrence (for $k$=2, 4) and negativity, and between $k$-ME concurrence and polynomial invariants also exist. Our work provides clear quantitative connections between $k$-ME concurrence and negativity, and between $k$-ME concurrence and polynomial invariants.


title: Acceleration of Actor-Critic Deep Reinforcement Learning for Visual Grasping in Clutter by State Representation Learning Based on Disentanglement of a Raw Input Image

abstract: For a robotic grasping task in which diverse unseen target objects exist in a cluttered environment, some deep learning-based methods have achieved state-of-the-art results using visual input directly. In contrast, actor-critic deep reinforcement learning (RL) methods typically perform very poorly when grasping diverse objects, especially when learning from raw images and sparse rewards. To make these RL techniques feasible for vision-based grasping tasks, we employ state representation learning (SRL), where we encode essential information first for subsequent use in RL. However, typical representation learning procedures are unsuitable for extracting pertinent information for learning the grasping skill, because the visual inputs for representation learning, where a robot attempts to grasp a target object in clutter, are extremely complex. We found that preprocessing based on the disentanglement of a raw input image is the key to effectively capturing a compact representation. This enables deep RL to learn robotic grasping skills from highly varied and diverse visual inputs. We demonstrate the effectiveness of this approach with varying levels of disentanglement in a realistic simulated environment.


title: Layered Sampling for Robust Optimization Problems

abstract: In real world, our datasets often contain outliers. Moreover, the outliers can seriously affect the final machine learning result. Most existing algorithms for handling outliers take high time complexities (e.g. quadratic or cubic complexity). {\em Coreset} is a popular approach for compressing data so as to speed up the optimization algorithms. However, the current coreset methods cannot be easily extended to handle the case with outliers. In this paper, we propose a new variant of coreset technique, {\em layered sampling}, to deal with two fundamental robust optimization problems: {\em $k$-median/means clustering with outliers} and {\em linear regression with outliers}. This new coreset method is in particular suitable to speed up the iterative algorithms (which often improve the solution within a local range) for those robust optimization problems. Moreover, our method is easy to be implemented in practice. We expect that our framework of layered sampling will be applicable to other robust optimization problems.


title: Globally optimal consensus maximization for robust visual inertial localization in point and line map

abstract: Map based visual inertial localization is a crucial step to reduce the drift in state estimation of mobile robots. The underlying problem for localization is to estimate the pose from a set of 3D-2D feature correspondences, of which the main challenge is the presence of outliers, especially in changing environment. In this paper, we propose a robust solution based on efficient global optimization of the consensus maximization problem, which is insensitive to high percentage of outliers. We first introduce translation invariant measurements (TIMs) for both points and lines to decouple the consensus maximization problem into rotation and translation subproblems, allowing for a two-stage solver with reduced solution dimensions. Then we show that (i) the rotation can be calculated by minimizing TIMs using only 1-dimensional branch-and-bound (BnB), (ii) the translation can be found by running 1-dimensional search for three times with prioritized progressive voting. Compared with the popular randomized solver, our solver achieves deterministic global convergence without depending on an initial value. While compared with existing BnB based methods, ours is exponentially faster. Finally, by evaluating the performance on both simulation and real-world datasets, our approach gives accurate pose even when there are 90\% outliers (only 2 inliers).


title: High-resolution imaging of Rydberg atoms in optical lattices using an aspheric-lens objective in vacuum

abstract: We present a high-resolution, simple and versatile system for imaging ultracold Rydberg atoms in optical lattices. The imaging objective is a single aspheric lens (with a working distance of 20.6 mm and a numerical aperture (NA) of 0.51) placed inside the vacuum chamber. Adopting a large-working-distance lens leaves room for electrodes and electrostatic shields to control electric fields around Rydberg atoms. With this setup, we achieve an Rayleigh resolution of 1.10 $\mu$m or $1.41\lambda$ ($\lambda=780$ nm), limited by the NA of the aspheric lens. For systems of highly excited Rydberg states with blockade radii greater than a few $\mu$m, the resolution achieved is sufficient for studying many physical processes of interest.


title: Cosmological Constraint on $\Omega_m$ and $\sigma_8$ from Cluster Abundances using the $\mathtt{GalWCat19}$ Optical-Spectroscopic SDSS Catalog

abstract: We derive cosmological constraints on the matter density, \om, and the amplitude of fluctuations, \sig, using $\mathtt{GalWCat19}$, a catalog of 1800 galaxy clusters we identified in the Sloan Digital Sky Survey-DR13 spectroscopic data set using our GalWeight technique to determine cluster membership \citep{Abdullah18,Abdullah19}. By analyzing a subsample of 756 clusters in a redshift range of $0.045\leq z \leq 0.125$ and virial masses of $M\geq 0.8\times10^{14}$ \hm ~with mean redshift of $z = 0.085$, we obtain \om ~$=0.310^{+0.023}_{-0.027} \pm 0.041$ (systematic) and \sig ~$=0.810^{+0.031}_{-0.036}\pm 0.035$ (systematic), with a cluster normalization relation of $\sigma_8= 0.43 \Omega_m^{-0.55}$. There are several unique aspects to our approach: we use the largest spectroscopic data set currently available, and we assign membership using the GalWeight technique which we have shown to be very effective at simultaneously maximizing the number of {\it{bona fide}} cluster members while minimizing the number of contaminating interlopers. Moreover, rather than employing scaling relations, we calculate cluster masses individually using the virial mass estimator. Since $\mathtt{GalWCat19}$ is a low-redshift cluster catalog we do not need to make any assumptions about evolution either in cosmological parameters or in the properties of the clusters themselves. Our constraints on \om ~and \sig ~are consistent and very competitive with those obtained from non-cluster abundance cosmological probes such as Cosmic Microwave Background (CMB), Baryonic Acoustic Oscillation (BAO), and supernovae (SNe). The joint analysis of our cluster data with Planck18+BAO+Pantheon gives \om ~$=0.315^{+0.013}_{-0.011}$ and \sig ~$=0.810^{+0.011}_{-0.010}$.


title: A variable service rate queue model for hub median problem

abstract: Hub location problems have multiple applications in logistic systems, airways industry, supply chain network design, and telecommunication. In the hub location problem, a number of nodes should be selected as the hub nodes to act as the main distributors and other nodes are connected together by these hubs. The input flow to the hub nodes is very large so more often we face the congestion in the hub nodes that causes disturbances in the whole system. Also, we have different service rates which is another cause of disturbance and should be addressed by the models. This paper addresses these issues by providing a model that prevents congestion in the system. We incorporated queuing system in the p-median hub location problem by considering multiple server options and different service rates. CAB dataset (contains 25 US cities) was used in the implementation and our findings show the big impact of considering congestion on the hub location network design.


title: Improving the Performance of Stochastic Local Search for Maximum Vertex Weight Clique Problem Using Programming by Optimization

abstract: The maximum vertex weight clique problem (MVWCP) is an important generalization of the maximum clique problem (MCP) that has a wide range of real-world applications. In situations where rigorous guarantees regarding the optimality of solutions are not required, MVWCP is usually solved using stochastic local search (SLS) algorithms, which also define the state of the art for solving this problem. However, there is no single SLS algorithm which gives the best performance across all classes of MVWCP instances, and it is challenging to effectively identify the most suitable algorithm for each class of MVWCP instances. In this work, we follow the paradigm of Programming by Optimization (PbO) to develop a new, flexible and highly parametric SLS framework for solving MVWCP, combining, for the first time, a broad range of effective heuristic mechanisms. By automatically configuring this PbO-MWC framework, we achieve substantial advances in the state-of-the-art in solving MVWCP over a broad range of prominent benchmarks, including two derived from real-world applications in transplantation medicine (kidney exchange) and assessment of research excellence.


title: Integrating Boundary Assembling into a DNN Framework for Named Entity Recognition in Chinese Social Media Text

abstract: Named entity recognition is a challenging task in Natural Language Processing, especially for informal and noisy social media text. Chinese word boundaries are also entity boundaries, therefore, named entity recognition for Chinese text can benefit from word boundary detection, outputted by Chinese word segmentation. Yet Chinese word segmentation poses its own difficulty because it is influenced by several factors, e.g., segmentation criteria, employed algorithm, etc. Dealt improperly, it may generate a cascading failure to the quality of named entity recognition followed. In this paper we integrate a boundary assembling method with the state-of-the-art deep neural network model, and incorporate the updated word boundary information into a conditional random field model for named entity recognition. Our method shows a 2% absolute improvement over previous state-of-the-art results.