title: Edge states and topological pumping in elastic lattices with periodically modulated coupling

abstract: We investigate the dispersion topology of elastic lattices characterized by spatial stiffness modulation. The modulation is defined by the sampling of a two-dimensional surface, which provides the lattices with topological properties that are usually attributed to two-dimensional crystals. We show that the cyclic variation of the phase of the stiffness modulation leads to a Berry phase accumulation for the Bloch eigenmodes, which is characterized by integer valued Chern numbers for the bands and associated gap labels. The resulting non-trivial gaps are spanned by edge modes localized at one of the boundaries of the considered 1D lattices. The edge mode location is governed by the phase of the stiffness distribution, whose spatial modulation drives these modes to transition from one edge to the other, through a bulk state that occurs when the corresponding dispersion branch touches the bulk bands. This property enables the implementation of a topological pump that is obtained by stacking and coupling a family of modulated 1D lattices along a second spatial dimension. The gradual variation of the stiffness phase modulation drives the adiabatic transition of the edge states, which transition from their localized state at one boundary, to a bulk mode and, finally, to another localized state at the opposite boundary. Similar effects were illustrated under the assumption of quasiperiodic modulation of the lattice interactions. We here demonstrate that a topological pump can be achieved also in periodic media, and illustrate this for the first time in elastic discrete lattices.


title: Characterizing Task Completion Latencies in Fog Computing

abstract: Fog computing, which distributes computing resources to multiple locations between the Internet of Things (IoT) devices and the cloud, is attracting considerable attention from academia and industry. Yet, despite the excitement about the potential of fog computing, few comprehensive quantitative characteristics of the properties of fog computing architectures have been conducted. In this paper we examine the properties of task completion latencies in fog computing. First, we present the results of our empirical benchmarking-based study of task completion latencies. The study covered a range of settings, and uniquely considered both traditional and serverless fog computing execution points. It demonstrated the range of execution point characteristics in different locations and the relative stability of latency characteristics for a given location. It also highlighted properties of serverless execution that are not incorporated in existing fog computing algorithms. Second, we present a framework we developed for co-optimizing task completion quality and latency, which was inspired by the insights of our empirical study. We describe fog computing task assignment problems we formulated under this framework, and present the algorithms we developed for solving them.


title: Demystifying Neural Network Filter Pruning

abstract: Based on filter magnitude ranking (e.g. L1 norm), conventional filter pruning methods for Convolutional Neural Networks (CNNs) have been proved with great effectiveness in computation load reduction. Although effective, these methods are rarely analyzed in a perspective of filter functionality. In this work, we explore the filter pruning and the retraining through qualitative filter functionality interpretation. We find that the filter magnitude based method fails to eliminate the filters with repetitive functionality. And the retraining phase is actually used to reconstruct the remained filters for functionality compensation for the wrongly-pruned critical filters. With a proposed functionality-oriented pruning method, we further testify that, by precisely addressing the filter functionality redundancy, a CNN can be pruned without considerable accuracy drop, and the retraining phase is unnecessary.


title: Deep Probabilistic Ensembles: Approximate Variational Inference through KL Regularization

abstract: In this paper, we introduce Deep Probabilistic Ensembles (DPEs), a scalable technique that uses a regularized ensemble to approximate a deep Bayesian Neural Network (BNN). We do so by incorporating a KL divergence penalty term into the training objective of an ensemble, derived from the evidence lower bound used in variational inference. We evaluate the uncertainty estimates obtained from our models for active learning on visual classification. Our approach steadily improves upon active learning baselines as the annotation budget is increased.


title: Building Corpora for Single-Channel Speech Separation Across Multiple Domains

abstract: To date, the bulk of research on single-channel speech separation has been conducted using clean, near-field, read speech, which is not representative of many modern applications. In this work, we develop a procedure for constructing high-quality synthetic overlap datasets, necessary for most deep learning-based separation frameworks. We produced datasets that are more representative of realistic applications using the CHiME-5 and Mixer 6 corpora and evaluate standard methods on this data to demonstrate the shortcomings of current source-separation performance. We also demonstrate the value of a wide variety of data in training robust models that generalize well to multiple conditions.


title: Computational Histological Staining and Destaining of Prostate Core Biopsy RGB Images with Generative Adversarial Neural Networks

abstract: Histopathology tissue samples are widely available in two states: paraffin-embedded unstained and non-paraffin-embedded stained whole slide RGB images (WSRI). Hematoxylin and eosin stain (H&E) is one of the principal stains in histology but suffers from several shortcomings related to tissue preparation, staining protocols, slowness and human error. We report two novel approaches for training machine learning models for the computational H&E staining and destaining of prostate core biopsy RGB images. The staining model uses a conditional generative adversarial network that learns hierarchical non-linear mappings between whole slide RGB image (WSRI) pairs of prostate core biopsy before and after H&E staining. The trained staining model can then generate computationally H&E-stained prostate core WSRIs using previously unseen non-stained biopsy images as input. The destaining model, by learning mappings between an H&E stained WSRI and a non-stained WSRI of the same biopsy, can computationally destain previously unseen H&E-stained images. Structural and anatomical details of prostate tissue and colors, shapes, geometries, locations of nuclei, stroma, vessels, glands and other cellular components were generated by both models with structural similarity indices of 0.68 (staining) and 0.84 (destaining). The proposed staining and destaining models can engender computational H&E staining and destaining of WSRI biopsies without additional equipment and devices.


title: Distilling Critical Paths in Convolutional Neural Networks

abstract: Neural network compression and acceleration are widely demanded currently due to the resource constraints on most deployment targets. In this paper, through analyzing the filter activation, gradients, and visualizing the filters' functionality in convolutional neural networks, we show that the filters in higher layers learn extremely task-specific features, which are exclusive for only a small subset of the overall tasks, or even a single class. Based on such findings, we reveal the critical paths of information flow for different classes. And by their intrinsic property of exclusiveness, we propose a critical path distillation method, which can effectively customize the convolutional neural networks to small ones with much smaller model size and less computation.


title: DeepDPM: Dynamic Population Mapping via Deep Neural Network

abstract: Dynamic high resolution data on human population distribution is of great importance for a wide spectrum of activities and real-life applications, but is too difficult and expensive to obtain directly. Therefore, generating fine-scaled population distributions from coarse population data is of great significance. However, there are three major challenges: 1) the complexity in spatial relations between high and low resolution population; 2) the dependence of population distributions on other external information; 3) the difficulty in retrieving temporal distribution patterns. In this paper, we first propose the idea to generate dynamic population distributions in full-time series, then we design dynamic population mapping via deep neural network(DeepDPM), a model that describes both spatial and temporal patterns using coarse data and point of interest information. In DeepDPM, we utilize super-resolution convolutional neural network(SRCNN) based model to directly map coarse data into higher resolution data, and a time-embedded long short-term memory model to effectively capture the periodicity nature to smooth the finer-scaled results from the previous static SRCNN model. We perform extensive experiments on a real-life mobile dataset collected from Shanghai. Our results demonstrate that DeepDPM outperforms previous state-of-the-art methods and a suite of frequent data-mining approaches. Moreover, DeepDPM breaks through the limitation from previous works in time dimension so that dynamic predictions in all-day time slots can be obtained.


title: Nonlinear Dynamics Semi-classical Model of Quantum Spin

abstract: A nonlinear dynamics semi-classical model is used to show that standard quantum spin analysis can be obtained. The model includes a classically driven nonlinear differential equation with dissipation and a semi-classical interpretation of the torque on a spin magnetic moment in the presence of a realistic magnetic field, which will represent two equilibrium positions. The highly complicated driven nonlinear dissipative semi-classical model is used to introduce chaos, which is necessary to produce the correct statistical quantum results. The resemblance between this semi-classical spin model and the thoroughly studied classical driven-damped nonlinear pendulum are shown and discussed.


title: A Particle-Hole-Symmetric Model for a Paired Fractional Quantum Hall State in a Half-filled Landau Level

abstract: The fractional quantum Hall effect (FQHE) observed at half filling of the second Landau level is believed to be caused by a pairing of composite fermions captured by the Moore-Read Pfaffian wave function. The generating Hamiltonian for the Moore-Read Pfaffian is a purely three-body model that breaks particle-hole symmetry and lacks other properties, such as dominate two-body repulsive interactions, expected from a physical model of the FQHE. We use exact diagonalization to study the low energy states of a more physical two-body generator model derived from the three-body model. We find that the two-body model exhibits the essential features expected from the Moore-Read Pfaffian: pairing, non-Abelian anyon excitations, and a neutral fermion mode. The model also satisfies constraints expected for a physical model of the FQHE at half-filling because it is: short range, spatially decaying, particle-hole symmetric, and supports a roton mode with a robust spectral gap in the thermodynamic limit. Hence, this two-body model offers a bridge between artificial three-body generator models for paired states and the physical Coulomb interaction and can be used to further explore properties of non-Abelian physics in the FQHE.


title: A random cocycle with non H\"older Lyapunov exponent

abstract: We provide an example of a Schr\"odinger cocycle over a mixing Markov shift for which the integrated density of states has a very weak modulus of continuity, close to the log-H\"older lower bound established by W. Craig and B. Simon. This model is based upon a classical example due to Y. Kifer of a random Bernoulli cocycle with zero Lyapunov exponents which is not strongly irreducible. It follows that the Lyapunov exponent of a Bernoulli cocycle near this Kifer example cannot be H\"older or weak-H\"older continuous, thus providing a limitation on the modulus of continuity of the Lyapunov exponent of random cocycles.


title: Spherical photon orbits in the field of Kerr naked singularities

abstract: For the Kerr naked singularity (KNS) spacetimes, we study properties of spherical photon orbits (SPOs) confined to constant Boyer-Lindquist radius r. Some new features of the SPOs are found, having no counterparts in the Kerr black hole (KBH) spacetimes, especially stable orbits that could be pure prograde/retrograde, or with turning point in the azimuthal direction. At r > 1 (r < 1) the covariant photon energy E > 0 (E < 0), at r = 1 there is E = 0. All unstable orbits must have E > 0. It is shown that the polar SPOs can exist only in the spacetimes with dimensionless spin a < 1.7996. Existence of closed SPOs with vanishing total change of the azimuth is demonstrated. Classification of the KNS and KBH spacetimes in dependence on their dimensionless spin a is proposed, considering the properties of the SPOs. For selected types of the KNS spacetimes, typical SPOs are constructed, including the closed paths. It is shown that the stable SPOs intersect the equatorial plane in a region of stable circular orbits of test particles, depending on the spin a. Relevance of this intersection for the Keplerian accretion discs is outlined and observational effects are estimated.


title: Non-minimally coupled Scalar Dark Matter from Inflationary Fluctuations

abstract: It is well known that light scalar fields present during inflation are coherently excited. We show that if the field couples to gravity in a non-minimal way, the fluctuations at large scales are suppressed with respect to the small scales ones. This fact allows for the field excitations to make a sizeable contribution to the energy density of the universe without generating too large isocurvature fluctuations at observable scales. We show that this mechanism could generate all the observed dark matter and study the main cosmological implications of this setup.


title: Visual Attention is Beyond One Single Saliency Map

abstract: Of later years, numerous bottom-up attention models have been proposed on different assumptions. However, the produced saliency maps may be different from each other even from the same input image. We also observe that human fixation map varies across time greatly. When people freely view an image, they tend to allocate attention at salient regions of large scale at first, and then search more and more detailed regions. In this paper, we argue that, for one input image visual attention cannot be described by only one single saliency map, and this mechanism should be modeled as a dynamic process. Under the frequency domain paradigm, we proposed a global inhibition model to mimic this process by suppressing the {\it non-saliency} in the input image; we also show that the dynamic process is influenced by one parameter in the frequency domain. Experiments illustrate that the proposed model is capable of predicting human dynamic fixation distribution.


title: Finding and Following of Honeycombing Regions in Computed Tomography Lung Images by Deep Learning

abstract: In recent years, besides the medical treatment methods in medical field, Computer Aided Diagnosis (CAD) systems which can facilitate the decision making phase of the physician and can detect the disease at an early stage have started to be used frequently. The diagnosis of Idiopathic Pulmonary Fibrosis (IPF) disease by using CAD systems is very important in that it can be followed by doctors and radiologists. It has become possible to diagnose and follow up the disease with the help of CAD systems by the development of high resolution computed imaging scanners and increasing size of computation power. The purpose of this project is to design a tool that will help specialists diagnose and follow up the IPF disease by identifying areas of honeycombing and ground glass patterns in High Resolution Computed Tomography (HRCT) lung images. Creating a program module that segments the lung pair and creating a self-learner deep learning model from given Computed Tomography (CT) images for the specific diseased regions thanks to doctors are the main purposes of this work. Through the created model, program module will be able to find special regions in given new CT images. In this study, the performance of lung segmentation was tested by the S{\o}rensen-Dice coefficient method and the mean performance was measured as 90.7%, testing of the created model was performed with data not used in the training stage of the CNN network, and the average performance was measured as 87.8% for healthy regions, 73.3% for ground-glass areas and 69.1% for honeycombing zones.


title: Planning Low-Carbon Campus Energy Hubs

abstract: Multi-energy systems can provide a constant level of service to end-use energy demands, while deriving delivered energy from a variety of primary/secondary energy sources. This fuel-switching capability can be used to reduce operating expenses, reduce environmental impacts, improve flexibility to accommodate renewable energy, and improve reliability. This paper presents four frameworks for incentivizing energy hub equipment investments for low-carbon operation targets. These frameworks vary in the measures taken to achieve low-carbon operation (explicit constraint vs. carbon pricing) and in the relationship between the hub builder and operator (cooperative vs. uncoordinated). The underlying energy hub model upon which these frameworks are built is an enhanced greenfield model, introducing 'energy buses' to reduce dimensionality. A case study is conducted for a campus being designed in Beijing, and results from each framework are compared to illustrate their relative costs. When the operator cannot be trusted to cooperate in controlling emissions, the system must be 'overbuilt' with more expensive equipment to ensure emissions target are met. A taxation-based approach increases overall costs at moderate emissions targets, but this effect decreases at aggressive targets. This paper also compares the cost of less efficient institutional frameworks with the most efficient approach, i.e. cooperation between builder and operator with constraints on emissions.


title: OverSketch: Approximate Matrix Multiplication for the Cloud

abstract: We propose OverSketch, an approximate algorithm for distributed matrix multiplication in serverless computing. OverSketch leverages ideas from matrix sketching and high-performance computing to enable cost-efficient multiplication that is resilient to faults and straggling nodes pervasive in low-cost serverless architectures. We establish statistical guarantees on the accuracy of OverSketch and empirically validate our results by solving a large-scale linear program using interior-point methods and demonstrate a 34% reduction in compute time on AWS Lambda.


title: A Volumetric Convolutional Neural Network for Brain Tumor Segmentation

abstract: Brain cancer can be very fatal, but chances of survival increase through early detection and treatment. Doctors use Magnetic Resonance Imaging (MRI) to detect and locate tumors in the brain, and very carefully analyze scans to segment brain tumors. Manual segmentation is time consuming and tiring for doctors, and it can be difficult for them to notice extremely small abnormalities. Automated segmentations performed by computers offer quicker diagnoses, the ability to notice small details, and more accurate segmentations. Advances in deep learning and computer hardware have allowed for high-performing automated segmentation approaches. However, several problems persist in practice: increased training time, class imbalance, and low performance. In this paper, I propose applying V-Net, a volumetric, fully convolutional neural network, to segment brain tumors in MRI scans from the BraTS Challenges. With this approach, I achieve a whole tumor dice score of 0.89 and train the network in a short time while addressing class imbalance with the use of a dice loss layer. Then, I propose applying an existing technique to improve automated segmentation performance in practice.


title: Sparse and Smooth Signal Estimation: Convexification of L0 Formulations

abstract: Signal estimation problems with smoothness and sparsity priors can be naturally modeled as quadratic optimization with $\ell_0$-"norm" constraints. Since such problems are non-convex and hard-to-solve, the standard approach is, instead, to tackle their convex surrogates based on $\ell_1$-norm relaxations. In this paper, we propose a new iterative (convex) conic quadratic relaxations that exploit not only the $\ell_0$-"norm" terms, but also the fitness and smoothness functions. The iterative convexification approach substantially closes the gap between the $\ell_0$-"norm" and its $\ell_1$ surrogate. These stronger relaxations lead to significantly better estimators than $\ell_1$-norm approaches and also allow one to utilize affine sparsity priors. In addition, the parameters of the model and the resulting estimators are easily interpretable. Experiments with a tailored Lagrangian decomposition method indicate that the proposed iterative convex relaxations \rev{yield solutions within 1\% of the exact $\ell_0$ approach, and can tackle instances with up to 100,000 variables under one minute.


title: Quaternion Convolutional Neural Networks for Heterogeneous Image Processing

abstract: Convolutional neural networks (CNN) have recently achieved state-of-the-art results in various applications. In the case of image recognition, an ideal model has to learn independently of the training data, both local dependencies between the three components (R,G,B) of a pixel, and the global relations describing edges or shapes, making it efficient with small or heterogeneous datasets. Quaternion-valued convolutional neural networks (QCNN) solved this problematic by introducing multidimensional algebra to CNN. This paper proposes to explore the fundamental reason of the success of QCNN over CNN, by investigating the impact of the Hamilton product on a color image reconstruction task performed from a gray-scale only training. By learning independently both internal and external relations and with less parameters than real valued convolutional encoder-decoder (CAE), quaternion convolutional encoder-decoders (QCAE) perfectly reconstructed unseen color images while CAE produced worst and gray-scale versions.


title: A Bayesian Perspective of Convolutional Neural Networks through a Deconvolutional Generative Model

abstract: Inspired by the success of Convolutional Neural Networks (CNNs) for supervised prediction in images, we design the Deconvolutional Generative Model (DGM), a new probabilistic generative model whose inference calculations correspond to those in a given CNN architecture. The DGM uses a CNN to design the prior distribution in the probabilistic model. Furthermore, the DGM generates images from coarse to finer scales. It introduces a small set of latent variables at each scale, and enforces dependencies among all the latent variables via a conjugate prior distribution. This conjugate prior yields a new regularizer based on paths rendered in the generative model for training CNNs-the Rendering Path Normalization (RPN). We demonstrate that this regularizer improves generalization, both in theory and in practice. In addition, likelihood estimation in the DGM yields training losses for CNNs, and inspired by this, we design a new loss termed as the Max-Min cross entropy which outperforms the traditional cross-entropy loss for object classification. The Max-Min cross entropy suggests a new deep network architecture, namely the Max-Min network, which can learn from less labeled data while maintaining good prediction performance. Our experiments demonstrate that the DGM with the RPN and the Max-Min architecture exceeds or matches the-state-of-art on benchmarks including SVHN, CIFAR10, and CIFAR100 for semi-supervised and supervised learning tasks.


title: When Not to Classify: Detection of Reverse Engineering Attacks on DNN Image Classifiers

abstract: This paper addresses detection of a reverse engineering (RE) attack targeting a deep neural network (DNN) image classifier; by querying, RE's aim is to discover the classifier's decision rule. RE can enable test-time evasion attacks, which require knowledge of the classifier. Recently, we proposed a quite effective approach (ADA) to detect test-time evasion attacks. In this paper, we extend ADA to detect RE attacks (ADA-RE). We demonstrate our method is successful in detecting "stealthy" RE attacks before they learn enough to launch effective test-time evasion attacks.


title: Machine Learning Algorithms for Classification of Microcirculation Images from Septic and Non-Septic Patients

abstract: Sepsis is a life-threatening disease and one of the major causes of death in hospitals. Imaging of microcirculatory dysfunction is a promising approach for automated diagnosis of sepsis. We report a machine learning classifier capable of distinguishing non-septic and septic images from dark field microcirculation videos of patients. The classifier achieves an accuracy of 89.45%. The area under the receiver operating characteristics of the classifier was 0.92, the precision was 0.92 and the recall was 0.84. Codes representing the learned feature space of trained classifier were visualized using t-SNE embedding and were separable and distinguished between images from critically ill and non-septic patients. Using an unsupervised convolutional autoencoder, independent of the clinical diagnosis, we also report clustering of learned features from a compressed representation associated with healthy images and those with microcirculatory dysfunction. The feature space used by our trained classifier to distinguish between images from septic and non-septic patients has potential diagnostic application.


title: Super-ideal diodes at the Schottky-Mott limit in gated graphene-WSe$_2$ heterojunctions

abstract: Metal-semiconductor interfaces, known as Schottky junctions, have long been hindered by defects and impurities. Such imperfections dominate the electrical characteristics of the junction by pinning the metal Fermi energy. We report measurements on a boron nitride encapsulated graphene-tungsten diselenide (WSe$_2$) Schottky junction which exhibits ideal diode characteristics and a complete lack of Fermi-level pinning. The Schottky barrier height of the device is rigidly tuned by electrostatic gating of the WSe$_2$, enabling experimental verification of the Schottky-Mott limit in a single device. Utilizing this exceptional gate control, we realize a super-ideal gated-Schottky diode which surpasses the ideal diode limit. Our results provide a pathway for defect-free electrical contact to two-dimensional semiconductors and open up possibilities for circuits with efficient switching characteristics and higher efficiency optoelectronic devices.


title: MAMMO: A Deep Learning Solution for Facilitating Radiologist-Machine Collaboration in Breast Cancer Diagnosis

abstract: With an aging and growing population, the number of women requiring either screening or symptomatic mammograms is increasing. To reduce the number of mammograms that need to be read by a radiologist while keeping the diagnostic accuracy the same or better than current clinical practice, we develop Man and Machine Mammography Oracle (MAMMO) - a clinical decision support system capable of triaging mammograms into those that can be confidently classified by a machine and those that cannot be, thus requiring the reading of a radiologist. The first component of MAMMO is a novel multi-view convolutional neural network (CNN) with multi-task learning (MTL). MTL enables the CNN to learn the radiological assessments known to be associated with cancer, such as breast density, conspicuity, suspicion, etc., in addition to learning the primary task of cancer diagnosis. We show that MTL has two advantages: 1) learning refined feature representations associated with cancer improves the classification performance of the diagnosis task and 2) issuing radiological assessments provides an additional layer of model interpretability that a radiologist can use to debug and scrutinize the diagnoses provided by the CNN. The second component of MAMMO is a triage network, which takes as input the radiological assessment and diagnostic predictions of the first network's MTL outputs and determines which mammograms can be correctly and confidently diagnosed by the CNN and which mammograms cannot, thus needing to be read by a radiologist. Results obtained on a private dataset of 8,162 patients show that MAMMO reduced the number of radiologist readings by 42.8% while improving the overall diagnostic accuracy in comparison to readings done by radiologists alone. We analyze the triage of patients decided by MAMMO to gain a better understanding of what unique mammogram characteristics require radiologists' expertise.


title: Similarity Learning with Higher-Order Graph Convolutions for Brain Network Analysis

abstract: Learning a similarity metric has gained much attention recently, where the goal is to learn a function that maps input patterns to a target space while preserving the semantic distance in the input space. While most related work focused on images, we focus instead on learning a similarity metric for neuroimages, such as fMRI and DTI images. We propose an end-to-end similarity learning framework called Higher-order Siamese GCN for multi-subject fMRI data analysis. The proposed framework learns the brain network representations via a supervised metric-based approach with siamese neural networks using two graph convolutional networks as the twin networks. Our proposed framework performs higher-order convolutions by incorporating higher-order proximity in graph convolutional networks to characterize and learn the community structure in brain connectivity networks. To the best of our knowledge, this is the first community-preserving similarity learning framework for multi-subject brain network analysis. Experimental results on four real fMRI datasets demonstrate the potential use cases of the proposed framework for multi-subject brain analysis in health and neuropsychiatric disorders. Our proposed approach achieves an average AUC gain of 75% compared to PCA, an average AUC gain of 65.5% compared to Spectral Embedding, and an average AUC gain of 24.3% compared to S-GCN across the four datasets, indicating promising application in clinical investigation and brain disease diagnosis.


title: Strong consistency of kernel estimator in a semiparametric regression model

abstract: Estimating the effective dimension reduction (EDR) space, related to the semiparametric regression model introduced by Li \cite{sir}, is based on the estimation of the covariance matrix $\Lambda$ of the conditional expectation of the vector of predictors given the response. An estimator $\widehat{\Lambda}_n$ of $\Lambda $ based on kernel method was introduced by Zhu and Fang \cite{Asymptotics} who then derived, under some conditions, the asymptotic distribution of $\sqrt{n}\left(\widehat{\Lambda}_n-\Lambda\right)$, as $n\rightarrow +\infty$. In this paper, we obtain, under specified conditions, the almost sure convergence of $\widehat{\Lambda}_n$ to $\Lambda$, as $n\rightarrow +\infty$.


title: An upper bound on the Wiener Index of a k-connected graph

abstract: The Wiener index of a connected graph is the summation of all distances between unordered pairs of vertices of the graph. In this paper, we give an upper bound on the Wiener index of a $k$-connected graph $G$ of order $n$ for integers $n-1>k \ge 1$: \[W(G) \le \frac{1}{4} n \lfloor \frac{n+k-2}{k} \rfloor (2n+k-2-k\lfloor \frac{n+k-2}{k} \rfloor).\] Moreover, we show that this upper bound is sharp when $k \ge 2$ is even, and can be obtained by the Wiener index of Harary graph $H_{k,n}$.


title: Quasi-incompressible Multi-species Ionic Fluid Models

abstract: In traditional hydrodynamic theories for ionic fluids, conservation of the mass and linear momentum is not properly taken care of. In this paper, we develop hydrodynamic theories for a viscous, ionic fluid of $N$ ionic species enforcing mass and momentum conservation as well as considering the size effect of the ionic particles. The theories developed are quasi-incompressible in that the mass-average velocity is no longer divergence-free whenever there exists variability in densities of the fluid components, and the models are dissipative. We present several ways to derive the transport equations for the ions, which lead to different rates of energy dissipation. The theories can be formulated in either number densities, volume fractions or mass densities of the ionic fluid components. We show that the theory with the Cahn-Hilliard transport equation for ionic species reduces to the classical Poisson-Nernst-Planck (PNP) model with the size effect for ionic fluids when the densities of the fluid components are equal and the entropy of the solvent is neglected. It further reduces to the PNP model when the size effect is neglected. A linear stability analysis of the model together with two of its limits, which is the extended PNP model (EPNP defined in the text) and the classical PNP model (CPNP) with the finite size effect, on a constant state and a comparison among the three models in 1D space are presented to highlight the similarity and the departure of this model from the EPNP and the CPNP model.


title: Topological Semantics for Lumped Parameter Systems Modeling

abstract: Behaviors of many engineering systems are described by lumped parameter models that encapsulate the spatially distributed nature of the system into networks of lumped elements; the dynamics of such a network is governed by a system of ordinary differential and algebraic equations. Languages and simulation tools for modeling such systems differ in syntax, informal semantics, and in the methods by which such systems of equations are generated and simulated, leading to numerous interoperability challenges. Logical extensions of SysML aim specially at unifying a subset of the underlying concepts in such languages. We propose to unify semantics of all such systems using standard notions from algebraic topology. In particular, Tonti diagrams classify all physical theories in terms of physical laws (topological and constitutive) defined over a pair of dual cochain complexes and may be used to describe different types of lumped parameter systems. We show that all possible methods for generating the corresponding state equations within each physical domain correspond to paths over Tonti diagrams. We further propose a generalization of Tonti diagram that captures the behavior and supports canonical generation of state equations for multi-domain lumped parameter systems. The unified semantics provides a basis for greater interoperability in systems modeling, supporting automated translation, integration, reuse, and numerical simulation of models created in different authoring systems and applications. Notably, the proposed algebraic topological semantics is also compatible with spatially and temporally distributed models that are at the core of modern CAD and CAE systems.


title: Band Selection from Hyperspectral Images Using Attention-based Convolutional Neural Networks

abstract: This paper introduces new attention-based convolutional neural networks for selecting bands from hyperspectral images. The proposed approach re-uses convolutional activations at different depths, identifying the most informative regions of the spectrum with the help of gating mechanisms. Our attention techniques are modular and easy to implement, and they can be seamlessly trained end-to-end using gradient descent. Our rigorous experiments showed that deep models equipped with the attention mechanism deliver high-quality classification, and repeatedly identify significant bands in the training data, permitting the creation of refined and extremely compact sets that retain the most meaningful features.


title: Automated Diagnosis of Lymphoma with Digital Pathology Images Using Deep Learning

abstract: Recent studies have shown promising results in using Deep Learning to detect malignancy in whole slide imaging. However, they were limited to just predicting positive or negative finding for a specific neoplasm. We attempted to use Deep Learning with a convolutional neural network algorithm to build a lymphoma diagnostic model for four diagnostic categories: benign lymph node, diffuse large B cell lymphoma, Burkitt lymphoma, and small lymphocytic lymphoma. Our software was written in Python language. We obtained digital whole slide images of Hematoxylin and Eosin stained slides of 128 cases including 32 cases for each diagnostic category. Four sets of 5 representative images, 40x40 pixels in dimension, were taken for each case. A total of 2,560 images were obtained from which 1,856 were used for training, 464 for validation and 240 for testing. For each test set of 5 images, the predicted diagnosis was combined from prediction of 5 images. The test results showed excellent diagnostic accuracy at 95% for image-by-image prediction and at 10% for set-by-set prediction. This preliminary study provided a proof of concept for incorporating automated lymphoma diagnostic screen into future pathology workflow to augment the pathologists' productivity.


title: Robust multiple-set linear canonical analysis based on minimum covariance determinant estimator

abstract: By deriving influence functions related to multiple-set linear canonical analysis (MSLCA) we show that the classical version of this analysis, based on empirical covariance operators, is not robust. Then, we introduce a robust version of MSLCA by using the MCD estimator of the covariance operator of the involved random vector. The related influence functions are then derived and are shown to be bounded. Asymptotic properties of the introduced robust MSLCA are obtained and permit to propose a robust test for mutual non-correlation. This test is shown to be robust by studying the related second order influence function under the null hypothesis.


title: Hausdorff closed limits and the causal boundary of globally hyperbolic spacetimes with timelike boundary

abstract: We show that when a spacetime $\mathcal{M}(=M \cup \partial M)$ is globally hyperbolic with (possibly empty) smooth timelike boundary $\partial M$, a metrizable topology, the closed limit topology (CLT) introduced by F. Hausdorff himself in the 1950's in set theory, can be advantageously adopted on the Geroch-Kronheimer-Penrose causal completion of M, retaining essentially all the good properties of previous topologies in this ambient. In particular, we show that if the globally hyperbolic spacetime $M$ admits a conformal boundary, defined in such broad terms as to include all the standard examples in the literature, then the latter is homeomorphic to the causal boundary endowed with the CLT. We also discuss how our recent proposal arXiv:1807.00152 for a definition of null infinity using only causal boundaries can be translated when using the CLT, simplifying a number of technical aspects in the pertinent definitions. In a more technical vein, in the appendix we discuss the relationship of the CLT with the more generally applicable (but not always Hausdorff) chronological topology, and show that they coincide exactly in those cases when the latter happens to be Hausdorff.


title: Angular reduction in multiparticle matrix elements

abstract: A general method for the reduction of coupled spherical harmonic products is presented. When the total angular coupling is zero, the reduction leads to an explicitly real expression in the scalar products within the unit vector arguments of the spherical harmonics. For non-scalar couplings, the reduction gives Cartesian tensor forms for the spherical harmonic products, with tensors built from the physical vectors in the original expression. The reduction for arbitrary couplings is given in closed form, making it amenable to symbolic manipulation on a computer. The final expressions do not depend on a special choice of coordinate axes, nor do they contain azimuthal quantum number summations, nor do they have complex tensor terms for couplings to a scalar. Consequently, they are easily interpretable from the properties of the physical vectors they contain.


title: CAB: Continuous Adaptive Blending Estimator for Policy Evaluation and Learning

abstract: The ability to perform offline A/B-testing and off-policy learning using logged contextual bandit feedback is highly desirable in a broad range of applications, including recommender systems, search engines, ad placement, and personalized health care. Both offline A/B-testing and off-policy learning require a counterfactual estimator that evaluates how some new policy would have performed, if it had been used instead of the logging policy. In this paper, we identify a family of counterfactual estimators which subsumes most such estimators proposed to date. Our analysis of this family identifies a new estimator - called Continuous Adaptive Blending (CAB) - which enjoys many advantageous theoretical and practical properties. In particular, it can be substantially less biased than clipped Inverse Propensity Score (IPS) weighting and the Direct Method, and it can have less variance than Doubly Robust and IPS estimators. In addition, it is sub-differentiable such that it can be used for learning, unlike the SWITCH estimator. Experimental results show that CAB provides excellent evaluation accuracy and outperforms other counterfactual estimators in terms of learning performance.


title: Gramian-Based Optimization for the Analysis and Control of Traffic Networks

abstract: This paper proposes a simplified version of classical models for urban transportation networks, and studies the problem of controlling intersections with the goal of optimizing network-wide congestion. Differently from traditional approaches to control traffic signaling, a simplified framework allows for a more tractable analysis of the network overall dynamics, and enables the design of critical parameters while considering network-wide measures of efficiency. Motivated by the increasing availability of real-time high-resolution traffic data, we cast an optimization problem that formalizes the goal of minimizing the overall network congestion by optimally controlling the durations of green lights at intersections. Our formulation allows us to relate congestion objectives with the problem of optimizing a metric of controllability of an associated dynamical network. We then provide a technique to efficiently solve the optimization by parallelizing the computation among a group of distributed agents. Lastly, we assess the benefits of the proposed modeling and optimization framework through microscopic simulations on typical traffic commute scenarios for the area of Manhattan. The optimization framework proposed in this study is made available online on a Sumo microscopic simulator based interface [1].


title: Electrically-Tunable High Curie Temperature Two-Dimensional Ferromagnetism in Van der Waals Layered Crystals

abstract: Identifying intrinsic low-dimensional ferromagnets with high transition temperature and electrically tunable magnetism is crucial for the development of miniaturized spintronics and magnetoelectrics. Recently long-range 2D ferromagnetism was observed in van der Waals crystals CrI$_3$ and Cr$_2$Ge$_2$Te$_6$, however their Curie temperature is significantly lowered when reducing down to monolayer/few layers. Herein, using renormalized spin-wave theory and first-principles electronic structure theory, we present a theoretical study of electrically tunable 2D ferromagnetism in van der Waals layered CrSBr and CrSeBr semiconductors with high Curie temperature of ~150K and sizable band gap. High transition temperature is attributed to strong anion-mediated superexchange interaction and a sizable spin-wave excitation gap due to large exchange and single-ion anisotropy. Remarkably, hole and electron doping can switch magnetization easy axis from in-plane to out-of-plane direction. These unique characteristics establish monolayer CrSBr and CrSeBr as promising platform for realizing 2D spintronics and magnetoelectrics such as 2D spin field effect transistor.


title: Poisson type operators on the Fock space of type B and in the Blitvi{\'c} model

abstract: In \cite{Bia97} Biane proposed a new statistic on set partitions which he called \emph{restricted crossings}. In a series of papers \cite{Ans01,Ans04,Ans04b,Ans05} Anshelevich showed that this statistic is an essential tool to investigate stochastic processes on $q$-Fock space. In particular, Anshelevich constructed operators whose moments count restricted crossings and used these operators to develop a beautiful theory of noncommutative $q$-L\'{e}vy processes. In the present paper following Anshelevich we define gauge operators on $(\alpha,q)$-Fock and cumulants which are governed by statistics on partitions of type B. In addition we investigate this construction in the context of a model of Blitvi{\'c} model \cite{B12}, where some related but different combinatorial structures appear, and we explain their relation with $t$-free probability.


title: Oblivious Set-maxima for Intersection of Convex Polygons

abstract: In this paper we revisit the well known set-maxima problem in the oblivious setting. Let $X=\{x_1,\ldots, x_n\}$ be a set of $n$ elements with an underlying total order. Let $\mathcal{S}=\{S_1,\ldots,S_m\}$ be a collection of $m$ distinct subsets of $X$. The set-maxima problem asks to determine the maxima of all the sets in the collection. In the comparison tree model we are interested in determining the number of comparisons necessary and sufficient to solve the problem. We present an oblivious algorithm based on the lattice structure of the input set system. Our algorithm is simple and yet for many set systems gives a non-trivial improvement over known deterministic algorithms. We apply our algorithm to a special $\cal S$ which is determined by an intersection structure of convex polygons and show that $O(n)$ comparisons suffice.


title: Ongoing Resurfacing of KBO Eris by Volatile Transport in Local, Collisional, Sublimation Atmosphere Regime

abstract: Kuiper belt object (KBO) Eris is exceptionally bright with a greater visible geometric albedo than any other known KBO. Its infrared reflectance spectrum is dominated by methane, which should form tholins that darken the surface on timescales much shorter than the age of the Solar System. Thus one or more ongoing processes probably maintain its brightness. Eris is predicted to have a primarily nitrogen atmosphere that is in vapor pressure equilibrium with nitrogen-ice and is collisional (not ballistic). Eris's eccentric orbit is expected to result in two atmospheric regimes: (1) a period near perihelion when the atmosphere is global (analogous to the atmospheres of Mars, Triton, and Pluto) and (2) a period near aphelion when only a local atmosphere exists near the warmest region (analogous to the atmosphere of Io). A numerical model developed to simulate Eris's thermal and volatile evolution in the local atmosphere regime is presented. The model conserves energy, mass, and momentum while maintaining vapor pressure equilibrium. It is adaptable to other local, collisional, sublimation atmospheres, which in addition to Io and Eris, may occur on several volatile-bearing KBOs. The model was applied for a limiting case where Eris is fixed at aphelion and has an initial nitrogen-ice mass everywhere equal to the precipitable column of nitrogen in Pluto's atmosphere during the New Horizons encounter (the resultant mass if the Pluto atmosphere collapsed uniformly onto the surface). The model results indicate that (1) transport of nitrogen in the local, collisional, sublimation atmosphere regime is significant, (2) changes of Eris's albedo or color from nitrogen transport may be observable, and (3) uniform collapse of a global, nitrogen atmosphere likely cannot explain Eris's anomalous albedo in the present epoch. Seasonal volatile transport remains a plausible hypothesis to explain...


title: Protein Hydration Waters are Susceptible to Unfavorable Perturbations

abstract: The interactions of a protein, its phase behavior, and ultimately, its ability to function, are all influenced by the interactions between the protein and its hydration waters. Here we study proteins with a variety of sizes, shapes, chemistries, and biological functions, and characterize their interactions with their hydration waters using molecular simulation and enhanced sampling techniques. We find that akin to extended hydrophobic surfaces, proteins situate their hydration waters at the edge of a dewetting transition, making them susceptible to unfavorable perturbations. We also find that the strength of the unfavorable potential needed to trigger dewetting is roughly the same, regardless of the protein being studied, and depends only the width of the hydration shell being perturbed. Our findings establish a framework for systematically classifying protein patches according to how favorably they interact with water.


title: Molecular ion abundances in the diffuse ISM : CF+, HCO+, HOC+, and C3H+

abstract: The transition between atomic and molecular hydrogen is associated with important changes in the structure of interstellar clouds, and marks the beginning of interstellar chemistry. Because of the relatively simple networks controlling their abundances, molecular ions are usually good probes of the underlying physical conditions including for instance the fraction of gas in molecular form or the fractional ionization. In this paper we focus on three possible probes of the molecular hydrogen column density, HCO+, HOC+, and CF+. We presented high sensitivity ALMA absorption data toward a sample of compact HII regions and bright QSOs with prominent foreground absorption, in the ground state transitions of the molecular ions HCO+, HOC+, and CF+ and the neutral species HCN and HNC, and from the excited state transitions of C3H+(4-3) and 13CS(2-1). These data are compared with Herschel absorption spectra of the ground state transition of HF and p-H2O. We show that the HCO+, HOC+, and CF+ column densities are well correlated with each other. HCO+ and HOC+ are tightly correlated with p-H2O, while they exhibit a different correlation pattern with HF depending on whether the absorbing matter is located in the Galactic disk or in the central molecular zone. We report new detections of C3H+ confirming that this ion is ubiquitous in the diffuse matter, with an abundance relative to H2 of ~7E-11. We confirm that the CF+ abundance is lower than predicted by simple chemical models and propose that the rate of the main formation reaction is lower by a factor of about 3 than usually assumed. In the absence of CH or HF data, we recommend to use the ground state transitions of HCO+, CCH, and HOC+ to trace diffuse molecular hydrogen, with mean abundances relative to H2 of 3E-9, 4E-8 and 4E-11.


title: Development of a SiPM-based readout-module for the characterization of various scintillator materials

abstract: In this paper we discuss a general approach to use SiPMs in different electrical configurations in combination with numerous scintillator materials for applications in nuclear instrumentation, such as calorimetry and timing. The "hybrid" configuration used in the PANDA detector was found to show the fastest timing properties with several hundred pico second response by combining the advantages of a large active area while avoiding high operation voltages and dark noise. A "parallel" configuration shows the largest light collection efficiency and is therefore suitable for energy measurements. Several calibration spectra were taken to measure the energy deposit of minimal ionizing cosmic muons in a two-centimeter-thick PbWO$_4$ crystal.


title: Proactive On-call Scheduling during a Seasonal Epidemic

abstract: Overcrowding in Emergency Departments (EDs) is particularly problematic during seasonal epidemic crises. Each year during this period, EDs set off recourse actions to cope with the increase in workload. Uncertainty in the length and amplitude of epidemics make managerial decisions difficult. We propose in this study a staff allocation model to manage the situation using on-calls. An on-call scheduling policy is proposed to best balance between demand coverage and labor cost under legal constraints of working time. The problem is modelled as a two-stage stochastic Integer Linear Program (ILP) and solved using a Sample Average Approximation (SAA) method. Several epidemic scenarios are defined with data from an ED in Lille, France.


title: Attention-Mechanism-based Tracking Method for Intelligent Internet of Vehicles

abstract: Vehicle tracking task plays an important role on the internet of vehicles and intelligent transportation system. Beyond the traditional GPS sensor, the image sensor can capture different kinds of vehicles, analyze their driving situation and can interact with them. Aiming at the problem that the traditional convolutional neural network is vulnerable to background interference, this paper proposes vehicle tracking method based on human attention mechanism for self-selection of deep features with an inter-channel fully connected layer. It mainly includes the following contents: 1) A fully convolutional neural network fused attention mechanism with the selection of the deep features for convolution. 2) A separation method for template and semantic background region to separate target vehicles from the background in the initial frame adaptively. 3) A two-stage method for model training using our traffic dataset. The experimental results show that the proposed method improves the tracking accuracy without an increase in tracking time. Meanwhile, it strengthens the robustness of algorithm under the condition of the complex background region. The success rate of the proposed method in overall traffic datasets is higher than Siamese network by about 10 percent, and the overall precision is higher than Siamese network by 8 percent.


title: Design and simplified calibration of a Mueller imaging polarimeter for material classification

abstract: This study is concerned with the design of a Mueller imaging polarimeter for the visualization of spatially-varying Mueller matrix fields. A simplified calibration procedure is advocated, where all the optical elements are calibrated simultaneously rather than independently as in the state-of-the-art. This is shown to significantly reduce the bias inherent to sequential calibration methods. In addition, this procedure requires no reference sample, it allows calibration both in transmission or in reflection modes, and it relies on ready-to-use cameras. Put together, these novelties should help non-specialists in optics designing and calibrating a Mueller imaging polarimeter for applications such as material classification.


title: Comb-resolved spectroscopy with immersion grating in long-wave infrared

abstract: We have developed a dispersive spectrometer using a compact immersion grating for direct frequency comb spectroscopy in the long-wave infrared region of 8-10 {\mu}m. A frequency resolution of 463 MHz is achieved, which is the highest reported in this wavelength region with a dispersive direct frequency comb spectrometer. We also demonstrate individual mode-resolved imaging of the frequency comb spectrum by Vernier cavity filtering and apply this to obtain both simple and complex molecular spectra. These results indicate that the immersion grating spectrometer offers the next advancement for sensitive, high-resolution spectroscopy of transient and large/complex molecules when combined with cavity enhancement and cooling techniques.


title: Flow-Cut Gaps and Face Covers in Planar Graphs

abstract: The relationship between the sparsest cut and the maximum concurrent multi-flow in graphs has been studied extensively. For general graphs with $k$ terminal pairs, the flow-cut gap is $O(\log k)$, and this is tight. But when topological restrictions are placed on the flow network, the situation is far less clear. In particular, it has been conjectured that the flow-cut gap in planar networks is $O(1)$, while the known bounds place the gap somewhere between $2$ (Lee and Raghavendra, 2003) and $O(\sqrt{\log k})$ (Rao, 1999). A seminal result of Okamura and Seymour (1981) shows that when all the terminals of a planar network lie on a single face, the flow-cut gap is exactly $1$. This setting can be generalized by considering planar networks where the terminals lie on $\gamma>1$ faces in some fixed planar drawing. Lee and Sidiropoulos (2009) proved that the flow-cut gap is bounded by a function of $\gamma$, and Chekuri, Shepherd, and Weibel (2013) showed that the gap is at most $3\gamma$. We prove that the flow-cut gap is $O(\log\gamma)$, by showing that the edge-weighted shortest-path metric induced on the terminals admits a stochastic embedding into trees with distortion $O(\log\gamma)$, which is tight. The preceding results refer to the setting of edge-capacitated networks. For vertex-capacitated networks, it can be significantly more challenging to control flow-cut gaps. While there is no exact vertex-capacitated version of the Okamura-Seymour Theorem, an approximate version holds; Lee, Mendel, and Moharrami (2015) showed that the vertex-capacitated flow-cut gap is $O(1)$ on planar networks whose terminals lie on a single face. We prove that the flow-cut gap is $O(\gamma)$ for vertex-capacitated instances when the terminals lie on at most $\gamma$ faces. In fact, this result holds in the more general setting of submodular vertex capacities.


title: Support-Area Dependence of Vibration-Insensitive Optical Cavities

abstract: The vibration sensitivities of optical cavities depending on the support-area were investigated both numerically and experimentally. We performed the numerical simulation with two models; one with total constraint over the support area, and the other with only vertical constraint. A support-area-size insensitive optimal support condition could be found by the numerical simulation. The support-area was determined in the experiment by a Viton rubber pad. The vertical, transverse, and longitudinal vibration sensitivities were measured experimentally. The experimental result agreed with the numerical simulation with a sliding model (only vertical constraint).


title: Finding Independent Transversals Efficiently

abstract: We give an efficient algorithm that, given a graph $G$ and a partition $V_1,\ldots,V_m$ of its vertex set, finds either an independent transversal (an independent set $\{v_1,\ldots,v_m\}$ in $G$ such that $v_i\in V_i$ for each $i$), or a subset $\mathcal B$ of vertex classes such that the subgraph of $G$ induced by $\bigcup\mathcal B$ has a small dominating set. A non-algorithmic proof of this result has been known for a number of years and has been applied to solve many other problems. Thus we are able to give algorithmic versions of many of these applications, a few of which we describe explicitly here.


title: Automatic Assessment of Full Left Ventricular Coverage in Cardiac Cine Magnetic Resonance Imaging with Fisher-Discriminative 3D CNN

abstract: Cardiac magnetic resonance (CMR) images play a growing role in the diagnostic imaging of cardiovascular diseases. Full coverage of the left ventricle (LV), from base to apex, is a basic criterion for CMR image quality and necessary for accurate measurement of cardiac volume and functional assessment. Incomplete coverage of the LV is identified through visual inspection, which is time-consuming and usually done retrospectively in the assessment of large imaging cohorts. This paper proposes a novel automatic method for determining LV coverage from CMR images by using Fisher-discriminative three-dimensional (FD3D) convolutional neural networks (CNNs). In contrast to our previous method employing 2D CNNs, this approach utilizes spatial contextual information in CMR volumes, extracts more representative high-level features and enhances the discriminative capacity of the baseline 2D CNN learning framework, thus achieving superior detection accuracy. A two-stage framework is proposed to identify missing basal and apical slices in measurements of CMR volume. First, the FD3D CNN extracts high-level features from the CMR stacks. These image representations are then used to detect the missing basal and apical slices. Compared to the traditional 3D CNN strategy, the proposed FD3D CNN minimizes within-class scatter and maximizes between-class scatter. We performed extensive experiments to validate the proposed method on more than 5,000 independent volumetric CMR scans from the UK Biobank study, achieving low error rates for missing basal/apical slice detection (4.9\%/4.6\%). The proposed method can also be adopted for assessing LV coverage for other types of CMR image data.


title: Training Domain Specific Models for Energy-Efficient Object Detection

abstract: We propose an end-to-end framework for training domain specific models (DSMs) to obtain both high accuracy and computational efficiency for object detection tasks. DSMs are trained with distillation \cite{hinton2015distilling} and focus on achieving high accuracy at a limited domain (e.g. fixed view of an intersection). We argue that DSMs can capture essential features well even with a small model size, enabling higher accuracy and efficiency than traditional techniques. In addition, we improve the training efficiency by reducing the dataset size by culling easy to classify images from the training set. For the limited domain, we observed that compact DSMs significantly surpass the accuracy of COCO trained models of the same size. By training on a compact dataset, we show that with an accuracy drop of only 3.6\%, the training time can be reduced by 93\%. The codes are uploaded in https://github.com/kentaroy47/training-domain-specific-models.


title: Evaluating Methods for End-User Creation of Robot Task Plans

abstract: How can we enable users to create effective, perception-driven task plans for collaborative robots? We conducted a 35-person user study with the Behavior Tree-based CoSTAR system to determine which strategies for end user creation of generalizable robot task plans are most usable and effective. CoSTAR allows domain experts to author complex, perceptually grounded task plans for collaborative robots. As a part of CoSTAR's wide range of capabilities, it allows users to specify SmartMoves: abstract goals such as "pick up component A from the right side of the table." Users were asked to perform pick-and-place assembly tasks with either SmartMoves or one of three simpler baseline versions of CoSTAR. Overall, participants found CoSTAR to be highly usable, with an average System Usability Scale score of 73.4 out of 100. SmartMove also helped users perform tasks faster and more effectively; all SmartMove users completed the first two tasks, while not all users completed the tasks using the other strategies. SmartMove users showed better performance for incorporating perception across all three tasks.


title: Optimal embeddings into Lorentz spaces for some vector differential operators via Gagliardo's lemma

abstract: We prove a family of Sobolev inequalities of the form $$ \Vert u \Vert_{L^{\frac{n}{n-1}, 1} (\mathbb{R}^n,V)} \le \Vert A (D) u \Vert_{L^1 (\mathbb{R}^n,E)} $$ where $A (D) : C^\infty_c (\mathbb{R}^n, V) \to C^\infty_c (\mathbb{R}^n, E)$ is a vector first-order homogeneous linear differential operator with constant coefficients, $u$ is a vector field on $\mathbb{R}^n$ and $L^{\frac{n}{n - 1}, 1} (\mathbb{R}^{n})$ is a Lorentz space. These new inequalities imply in particular the extension of the classical Gagliardo-Nirenberg inequality to Lorentz spaces originally due to Alvino and a sharpening of an inequality in terms of the deformation operator by Strauss (Korn-Sobolev inequality) on the Lorentz scale. The proof relies on a nonorthogonal application of the Loomis--Whitney inequality and Gagliardo's lemma.


title: A second countable locally compact transitive groupoid without open range map

abstract: Dana P. Williams raised in [Proc. Am. Math. Soc., Ser. B, 2016] the following question: Must a second countable, locally compact, transitive groupoid have open range map? This paper gives a negative answer to that question. Although a second countable, locally compact transitive groupoid G may fail to have open range map, we prove that we can replace its topology with a topology which is also second countable, locally compact, and with respect to which G is a topological groupoid whose range map is open. Moreover, the two topologies generate the same Borel structure and coincide on the fibres of G.


title: Deep Reinforcement Learning via L-BFGS Optimization

abstract: Reinforcement Learning (RL) algorithms allow artificial agents to improve their action selections so as to increase rewarding experiences in their environments. Deep Reinforcement Learning algorithms require solving a nonconvex and nonlinear unconstrained optimization problem. Methods for solving the optimization problems in deep RL are restricted to the class of first-order algorithms, such as stochastic gradient descent (SGD). The major drawback of the SGD methods is that they have the undesirable effect of not escaping saddle points and their performance can be seriously obstructed by ill-conditioning. Furthermore, SGD methods require exhaustive trial and error to fine-tune many learning parameters. Using second derivative information can result in improved convergence properties, but computing the Hessian matrix for large-scale problems is not practical. Quasi-Newton methods require only first-order gradient information, like SGD, but they can construct a low rank approximation of the Hessian matrix and result in superlinear convergence. The limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) approach is one of the most popular quasi-Newton methods that construct positive definite Hessian approximations. In this paper, we introduce an efficient optimization method, based on the limited memory BFGS quasi-Newton method using line search strategy -- as an alternative to SGD methods. Our method bridges the disparity between first order methods and second order methods by continuing to use gradient information to calculate a low-rank Hessian approximations. We provide formal convergence analysis as well as empirical results on a subset of the classic ATARI 2600 games. Our results show a robust convergence with preferred generalization characteristics, as well as fast training time and no need for the experience replaying mechanism.


title: Reconstructing Speech Stimuli From Human Auditory Cortex Activity Using a WaveNet Approach

abstract: The superior temporal gyrus (STG) region of cortex critically contributes to speech recognition. In this work, we show that a proposed WaveNet, with limited available data, is able to reconstruct speech stimuli from STG intracranial recordings. We further investigate the impulse response of the fitted model for each recording electrode and observe phoneme level temporospectral tuning properties for the recorded area of cortex. This discovery is consistent with previous studies implicating the posterior STG (pSTG) in a phonetic representation of speech and provides detailed acoustic features that certain electrode sites possibly extract during speech recognition.


title: Seismic Signal Denoising and Decomposition Using Deep Neural Networks

abstract: Denoising and filtering are widely used in routine seismic-data-processing to improve the signal-to-noise ratio (SNR) of recorded signals and by doing so to improve subsequent analyses. In this paper we develop a new denoising/decomposition method, DeepDenoiser, based on a deep neural network. This network is able to learn simultaneously a sparse representation of data in the time-frequency domain and a non-linear function that maps this representation into masks that decompose input data into a signal of interest and noise (defined as any non-seismic signal). We show that DeepDenoiser achieves impressive denoising of seismic signals even when the signal and noise share a common frequency band. Our method properly handles a variety of colored noise and non-earthquake signals. DeepDenoiser can significantly improve the SNR with minimal changes in the waveform shape of interest, even in presence of high noise levels. We demonstrate the effect of our method on improving earthquake detection. There are clear applications of DeepDenoiser to seismic imaging, micro-seismic monitoring, and preprocessing of ambient noise data. We also note that potential applications of our approach are not limited to these applications or even to earthquake data, and that our approach can be adapted to diverse signals and applications in other settings.


title: ACE: An Actor Ensemble Algorithm for Continuous Control with Tree Search

abstract: In this paper, we propose an actor ensemble algorithm, named ACE, for continuous control with a deterministic policy in reinforcement learning. In ACE, we use actor ensemble (i.e., multiple actors) to search the global maxima of the critic. Besides the ensemble perspective, we also formulate ACE in the option framework by extending the option-critic architecture with deterministic intra-option policies, revealing a relationship between ensemble and options. Furthermore, we perform a look-ahead tree search with those actors and a learned value prediction model, resulting in a refined value estimation. We demonstrate a significant performance boost of ACE over DDPG and its variants in challenging physical robot simulators.


title: Quantum Particle Production Effects on Cosmic Expansion

abstract: Quantum fields in cosmological spacetimes can experience particle production due to their interaction with the expanding background. This effect is particularly relevant for models of the very early Universe, when the energy density generated through this process may back-react on the cosmological expansion. Yet, these scenarios have not been fully explored due to the several technical hurdles imposed by the back-reaction calculations. In this work we review the basics of cosmological quantum particle production and demonstrate a numerical algorithm to solve the back-reaction problem in regimes dominated by particle production. As an illustration, we compute the effects of a massive quantized scalar field on a cosmological bounce scenario, explicitly showing that quantum particle production can cause the contracting phase to end in a radiation crunch, or can delay the bounce. Finally, we discuss the relevance of quantum particle production/annihilation to bounce and inflationary models of the early Universe.


title: Solutions of nonlinear Sobolev-Burgers PDEs

abstract: Nonlinear Sobolev-Burgers PDEs are considered. Their solutions are investigated. A technique of noncommutative line integration is utilized for their description. A new method of PDEs solution with the help of Cayley-Dickson algebras is developed in the article. Moreover, random operator valued measures are studied and applied to solutions of PDEs.


title: The Sitting Closer to Friends than Enemies Problem in the Circumference

abstract: The Sitting Closer to Friends than Enemies (SCFE) problem is to find an embedding in a metric space for the vertices of a given signed graph so that, for every pair of incident edges with different sign, the positive edge is shorter (in the metric of the space) than the negative edge. In this document, we present new results regarding the SCFE problem when the metric space in consideration is the circumference. Our main results say that, given a signed graph, it is NP-complete to decide whether such an embedding exists in the circumference or not. Nevertheless, if the given signed graph is complete, then such decision can be made in polynomial time. In particular, we prove that, given a complete signed graph, it has such an embedding if and only if its positive part is a proper circular arc graph.


title: First Measurement of $\nu_{\mu}$ Charged-Current $\pi^{0}$ Production on Argon with a LArTPC

abstract: We report the first measurement of the flux-integrated cross section of $\nu_{\mu}$ charged-current single $\pi^{0}$ production on argon. This measurement is performed with the MicroBooNE detector, an 85 ton active mass liquid argon time projection chamber exposed to the Booster Neutrino Beam at Fermilab. This result on argon is compared to past measurements on lighter nuclei to investigate the scaling assumptions used in models of the production and transport of pions in neutrino-nucleus scattering. The techniques used are an important demonstration of the successful reconstruction and analysis of neutrino interactions producing electromagnetic final states using a liquid argon time projection chamber operating at the earth's surface.


title: Proceedings of the 2018 Workshop on Compositional Approaches in Physics, NLP, and Social Sciences

abstract: The ability to compose parts to form a more complex whole, and to analyze a whole as a combination of elements, is desirable across disciplines. This workshop bring together researchers applying compositional approaches to physics, NLP, cognitive science, and game theory. Within NLP, a long-standing aim is to represent how words can combine to form phrases and sentences. Within the framework of distributional semantics, words are represented as vectors in vector spaces. The categorical model of Coecke et al. [2010], inspired by quantum protocols, has provided a convincing account of compositionality in vector space models of NLP. There is furthermore a history of vector space models in cognitive science. Theories of categorization such as those developed by Nosofsky [1986] and Smith et al. [1988] utilise notions of distance between feature vectors. More recently G\"ardenfors [2004, 2014] has developed a model of concepts in which conceptual spaces provide geometric structures, and information is represented by points, vectors and regions in vector spaces. The same compositional approach has been applied to this formalism, giving conceptual spaces theory a richer model of compositionality than previously [Bolt et al., 2018]. Compositional approaches have also been applied in the study of strategic games and Nash equilibria. In contrast to classical game theory, where games are studied monolithically as one global object, compositional game theory works bottom-up by building large and complex games from smaller components. Such an approach is inherently difficult since the interaction between games has to be considered. Research into categorical compositional methods for this field have recently begun [Ghani et al., 2018]. Moreover, the interaction between the three disciplines of cognitive science, linguistics and game theory is a fertile ground for research. Game theory in cognitive science is a well-established area [Camerer, 2011]. Similarly game theoretic approaches have been applied in linguistics [J\"ager, 2008]. Lastly, the study of linguistics and cognitive science is intimately intertwined [Smolensky and Legendre, 2006, Jackendoff, 2007]. Physics supplies compositional approaches via vector spaces and categorical quantum theory, allowing the interplay between the three disciplines to be examined.


title: Greedy Frank-Wolfe Algorithm for Exemplar Selection

abstract: In this paper, we consider the problem of selecting representatives from a data set for arbitrary supervised/unsupervised learning tasks. We identify a subset $S$ of a data set $A$ such that 1) the size of $S$ is much smaller than $A$ and 2) $S$ efficiently describes the entire data set, in a way formalized via convex optimization. In order to generate $|S| = k$ exemplars, our kernelizable algorithm, Frank-Wolfe Sparse Representation (FWSR), only needs to execute $\approx k$ iterations with a per-iteration cost that is quadratic in the size of $A$. This is in contrast to other state of the art methods which need to execute until convergence with each iteration costing an extra factor of $d$ (dimension of the data). Moreover, we also provide a proof of linear convergence for our method. We support our results with empirical experiments; we test our algorithm against current methods in three different experimental setups on four different data sets. FWSR outperforms other exemplar finding methods both in speed and accuracy in almost all scenarios.


abstract: We developed and experimentally tested a Symmetric Traveling Wave Parametric Amplifier (STWPA) based on Three-Wave Mixing, using the new concept of a Symmetric rf- SQUID. This allows to fully control the second and third order nonlinearities of the STWPA by applying external currents. In this way, the optimal bias point can be reached, taking into account both phase mismatch and pump depletion minimization. The structure was tested at 4.2K, showing a 4GHz bandwidth and a maximum estimated gain of 17dB. STWPA showed also great flexibility, allowing up-down conversion mixer operations and rf-controlled switch.


title: The Euler Equations in Planar Domains with Corners

abstract: When the velocity field is not a priori known to be globally almost Lipschitz, global uniqueness of solutions to the two-dimensional Euler equations has been established only in some special cases, and the solutions to which these results apply share the property that the diffuse part of the vorticity is constant near the points where the velocity is insufficiently regular. Assuming that the latter holds initially, the challenge is then to propagate this property along the Euler dynamic via an appropriate control of the Lagrangian trajectories. In domains with obtuse corners and sufficiently smooth elsewhere, Yudovich solutions fail to be almost Lipschitz only near these corners, and we investigate necessary and sufficient conditions for the vorticity to remain constant there. We show that if the vorticity is initially constant near the whole boundary, then it remains such forever (and global weak solutions are unique), provided no corner has angle greater than $\pi$. We also show that this fails in general for domains that do have such corners.


title: Heuristic Recurrent Algorithms for Photonic Ising Machines

abstract: The inability of conventional electronic architectures to efficiently solve large combinatorial problems motivates the development of novel computational hardware. There has been much effort recently toward developing novel, application-specific hardware, across many different fields of engineering, such as integrated circuits, memristors, and photonics. However, unleashing the true potential of such novel architectures requires the development of featured algorithms which optimally exploit their fundamental properties. We here present the Photonic Recurrent Ising Sampler (PRIS), a heuristic method tailored for parallel architectures that allows for fast and efficient sampling from distributions of combinatorially hard Ising problems. Since the PRIS relies essentially on vector-to-fixed matrix multiplications, we suggest the implementation of the PRIS in photonic parallel networks, which realize these operations at an unprecedented speed. The PRIS provides sample solutions to the ground state of arbitrary Ising models, by converging in probability to their associated Gibbs distribution. By running the PRIS at various noise levels, we probe the critical behavior of universality classes and their critical exponents. In addition to the attractive features of photonic networks, the PRIS relies on intrinsic dynamic noise and eigenvalue dropout to find ground states more efficiently. Our work suggests speedups in heuristic methods via photonic implementations of the PRIS. We also hint at a broader class of (meta)heuristic algorithms derived from the PRIS, such as combined simulated annealing on the noise and eigenvalue dropout levels. Our algorithm can also be implemented in a competitive manner on fast parallel electronic hardware, such as FPGAs and ASICs.


title: The planning problem in Mean Field Games as regularized mass transport

abstract: In this paper, using variational approaches, we investigate the first order planning problem arising in the theory of mean field games. We show the existence and uniqueness of weak solutions of the problem in the case of a large class of Hamiltonians with arbitrary superlinear order of growth at infinity and local coupling functions. We require the initial and final measures to be merely summable. At the same time (relying on the techniques developed recently by Graber and M\'esz\'aros), under stronger monotonicity and convexity conditions on the data, we obtain Sobolev estimates on the solutions of the planning problem both for space and time derivatives.


title: Investigating First Returns: The Effect of Multicolored Vectors

abstract: By definition, a first return is the immediate moment that a path, using vectors in the Cartesian plane, touches the $x$-axis after leaving it previously from a given point; the initial point is often the origin. In this case, using certain diagonal and horizontal vectors while restricting the movements to the first quadrant will cause almost every first return to end at the point $(2n,0)$, where $2n$ counts the equal number of up and down steps in a path. The exception will be explained further in the sections below. Using the first returns of Catalan, Schr\"{o}der, and Motzkin numbers, which resulted from the lattice paths formed using a combination of diagonal and/or horizontal vectors, we then investigated the effect that coloring select vectors will have on each of the original generating functions.


title: The Relationship Between Pascal's Triangle and Random Walks

abstract: Random walks are a series of up, down, and level steps that enumerate distinct paths from $(0,0)$ to $(2n,0)$, where $n$ is the semi-length of the path. We used these paths to analyze Catalan, Schr\"{o}der, and Motzkin number sequences through a combination of matrix operations, quadratic functions, and inductive reasoning. Our results revealed a number of distinct patterns, some unnamed, between these number sequences and Pascal's triangle that can be explained through generating functions, first returns, group theory, and the Riordan matrix. Various proofs and properties of these number sequences are provided, including each generating function, their respective first returns, and matrix properties. These findings lead to a deeper understanding of combinatorics and graph theory.


title: Global well-posedness and asymptotic behavior in Besov-Morrey spaces for chemotaxis-Navier-Stokes fluids

abstract: In this work we consider the Keller-Segel system coupled with Navier-Stokes equations in $\mathbb{R}^{N}$ for $N\geq2$. We prove the global well-posedness with small initial data in Besov-Morrey spaces. Our initial data class extends previous ones found in the literature such as that obtained by Kozono-Miura-Sugiyama (J. Funct. Anal. 2016). It allows to consider initial cell density and fluid velocity concentrated on smooth curves or at points depending on the spatial dimension. Self-similar solutions are obtained depending on the homogeneity of the initial data and considering the case of chemical attractant without degradation rate. Moreover, we analyze the asymptotic stability of solutions at infinity and obtain a class of asymptotically self-similar ones.


title: Abstract hypernormalisation, and normalisation-by-trace-evaluation for generative systems

abstract: Jacobs' hypernormalisation is a construction on finitely supported discrete probability distributions, obtained by generalising certain patterns occurring in quantitative information theory. In this paper, we generalise Jacobs' notion in turn, by describing a notion of hypernormalisation in the abstract setting of a symmetric monoidal category endowed with a linear exponential monad -- a structure arising in the categorical semantics of linear logic. We show that Jacobs' hypernormalisation arises in this fashion from the finitely supported probability measure monad on the category of sets, which can be seen as a linear exponential monad with respect to a non-standard monoidal structure on sets which we term the convex monoidal structure. We give the construction of this monoidal structure in terms of a quantum-algebraic notion known as a tricocycloid. Besides the motivating example, and its natural generalisations to the continuous context, we give a range of other instances of our abstract hypernormalisation, which swap out the side-effect of probabilistic choice for other important side-effects such as non-deterministic choice, logical choice via tests in a Boolean algebra, and input from a stream of values. Finally, we exploit our framework to describe a normalisation-by-trace-evaluation process for behaviours of various kinds of coalgebraic generative systems, including labelled transition systems, probabilistic generative systems, and stream processors.


title: Resource-efficient analyzer of Bell and Greenberger-Horne-Zeilinger states of multiphoton systems

abstract: We propose a resource-efficient error-rejecting entangled-state analyzer for polarization-encoded multiphoton systems. Our analyzer is based on two single-photon quantum-nondemolition detectors, where each of them is implemented with a four-level emitter (e.g., a quantum dot) coupled to a one-dimensional system (such as a micropillar cavity or a photonic nanocrystal waveguide). The analyzer works in a passive way and can completely distinguish $2^n$ Greenberger-Horne-Zeilinger~(GHZ) states of $n$ photons without using any active operation or fast switching. The efficiency and fidelity of the GHZ-state analysis can, in principle, be close to unity, when an ideal single-photon scattering condition is fulfilled. For a nonideal scattering, which typically reduces the fidelity of a GHZ-state analysis, we introduce a passively error-rejecting circuit to enable a near-perfect fidelity at the expense of a slight decrease of its efficiency. Furthermore, the protocol can be directly used to perform a two-photon Bell-state analysis. This passive, resource-efficient, and error-rejecting protocol can, therefore, be useful for practical quantum networks.


title: Polyadic Hopf algebras and quantum groups

abstract: This article continues the study of concrete algebra-like structures in our polyadic approach, where the arities of all operations are initially taken as arbitrary, but the relations between them, the arity shapes, are to be found from some natural conditions ("arity freedom principle"). In this way, generalized associative algebras, coassociative coalgebras, bialgebras and Hopf algebras are defined and investigated. They have many unusual features in comparison with the binary case. For instance, both the algebra and its underlying field can be zeroless and nonunital, the existence of the unit and counit is not obligatory, and the dimension of the algebra is not arbitrary, but "quantized". The polyadic convolution product and bialgebra can be defined, and when the algebra and coalgebra have unequal arities, the polyadic version of the antipode, the querantipode, has different properties. As a possible application to quantum group theory, we introduce the polyadic version of braidings, almost co-commutativity, quasitriangularity and the equations for the $R$-matrix (which can be treated as a polyadic analog of the Yang-Baxter equation). Finally, we propose another concept of deformation which is governed not by the twist map, but by the medial map, where only the latter is unique in the polyadic case. We present the corresponding braidings, almost co-mediality and $M$-matrix, for which the compatibility equations are found.


title: Ferroelectricity in [111]-oriented epitaxially strained SrTiO$_3$ from first principles

abstract: We use first principles density functional theory calculations to investigate the effect of biaxial strain in the low-temperature structural and ferroelectric properties of [111]-oriented SrTiO$_3$. We find that [111] biaxial strain, achievable by coherent epitaxial growth along the [111] direction, induces structural distortions in SrTiO$_3$ that are not present in either bulk or [001]-oriented SrTiO$_3$. Under [111] biaxial strain, SrTiO$_3$ displays ferroelectricity at tensile strain, and paraelectricity at compressive strain. We compute the phonon spectrum and macroscopic polarization of SrTiO$_3$ as a function of [111] biaxial strain, and relate our results to the predictions of the free energy phenomenological model of Pertsev, Tagantsev and Setter [Phys. Rev. B 61, 825 (2000); Phys. Rev. B 65, 219901 (2002)].


title: The RLLChatbot: a solution to the ConvAI challenge

abstract: Current conversational systems can follow simple commands and answer basic questions, but they have difficulty maintaining coherent and open-ended conversations about specific topics. Competitions like the Conversational Intelligence (ConvAI) challenge are being organized to push the research development towards that goal. This article presents in detail the RLLChatbot that participated in the 2017 ConvAI challenge. The goal of this research is to better understand how current deep learning and reinforcement learning tools can be used to build a robust yet flexible open domain conversational agent. We provide a thorough description of how a dialog system can be built and trained from mostly public-domain datasets using an ensemble model. The first contribution of this work is a detailed description and analysis of different text generation models in addition to novel message ranking and selection methods. Moreover, a new open-source conversational dataset is presented. Training on this data significantly improves the Recall@k score of the ranking and selection mechanisms compared to our baseline model responsible for selecting the message returned at each interaction.


title: First assessment of the binary lens OGLE-2015_BLG-0232

abstract: We present an analysis of the microlensing event OGLE-2015-BLG-0232. This event is challenging to characterize for two reasons. First, the light curve is not well sampled during the caustic crossing due to the proximity of the full Moon impacting the photometry quality. Moreover, the source brightness is difficult to estimate because this event is blended with a nearby K dwarf star. We found that the light curve deviations are likely due to a close brown dwarf companion (i.e., s = 0.55 and q = 0.06), but the exact nature of the lens is still unknown. We finally discuss the potential of follow-up observations to estimate the lens mass and distance in the future.


title: Dark Matter microphysics and 21 cm observations

abstract: Dark matter interactions with massless or very light Standard Model particles, as photons or neutrinos, may lead to a suppression of the matter power spectrum at small scales and of the number of low mass haloes. Bounds on the dark matter scattering cross section with light degrees of freedom in such interacting dark matter (IDM) scenarios have been obtained from e.g. early time cosmic microwave background physics and large scale structure observations. Here we scrutinize dark matter microphysics in light of the claimed 21 cm EDGES 78 MHz absorption signal. IDM is expected to delay the 21 cm absorption features due to collisional damping effects. We identify the astrophysical conditions under which the existing constraints on the dark matter scattering cross section could be largely improved due to the IDM imprint on the 21 cm signal, providing also an explicit comparison to the WDM scenario.


title: Leibniz type rule: $\Psi-$Hilfer fractional derivative

abstract: In this paper, we present the Leibniz rule for the $\Psi-$Hilfer ($\Psi-$H) fractional derivative in two versions, the first in relation to $\Psi-$RL fractional derivative and the second in relation to the $\Psi-$H fractional derivative. In this sense, we present some particular cases of Leibniz rules and Leibniz type rules from the investigated case.


title: Two-dimensional composite solitons in a spin-orbit-coupled Fermi gas in free space

abstract: We address a possibility of creating soliton states in oblate binary-fermionic clouds in the framework of the density-functional theory, which includes the spin-orbit coupling (SOC) and nonlinear attraction between spin-up and down-polarized components of the spinor wave function. In the limit when the inter-component attraction is much stronger than the effective intra-component Pauli repulsion, the resulting model also represents a system of Gross-Pitaevskii equations for a binary Bose-Einstein condensate including the SOC effect. We show that the model gives rise to two-dimensional quiescent composite solitons in free space. A stability region is identified for solitons of the mixed-mode type (which feature mixtures of zero-vorticity and vortical terms in both components), while solitons of the other type, semi-vortices (with the vorticity carried by one component) are unstable. Due to breaking of the Galilean invariance by SOC, the system supports moving solitons with velocities up to a specific critical value. Collisions between moving solitons are briefly considered too. The collisions lead, in particular, to a quasi-elastic rebound, or an inelastic outcome, which features partial merger of the solitons.


title: Expansion dynamics of a spherical Bose-Einstein condensate

abstract: We experimentally and theoretically observe the expansion behaviors of a spherical Bose-Einstein condensate. A rubidium condensate is produced in an isotropic optical dipole trap with an asphericity of 0.037. We measure the variation of the condensate size during the expansion process. The free expansion of the condensate is isotropic, which is different from that of the condensate usually produced in the anisotropic trap. The expansion in the short time is speeding and then after a long time the expansion velocity asymptotically approaches a constant value. We derive an analytic solution of the expansion behavior based on the spherical symmetry, allowing a quantitative comparison with the experimental measurement. The interaction energy of the condensate is gradually converted into the kinetic energy at the beginning of the expansion and the kinetic energy dominates after a long-time expansion. We obtain the interaction energy of the condensate in the trap by probing the expansion velocity, which is consistent with the theoretical prediction.


title: Zernike Polynomials: Evaluation, Quadrature, and Interpolation

abstract: Zernike polynomials are a basis of orthogonal polynomials on the unit disk that are a natural basis for representing smooth functions. They arise in a number of applications including optics and atmospheric sciences. In this paper, we provide a self-contained reference on Zernike polynomials, algorithms for evaluating them, and what appear to be new numerical schemes for quadrature and interpolation. We also introduce new properties of Zernike polynomials in higher dimensions. The quadrature rule and interpolation scheme use a tensor product of equispaced nodes in the angular direction and roots of certain Jacobi polynomials in the radial direction. An algorithm for finding the roots of these Jacobi polynomials is also described. The performance of the interpolation and quadrature schemes is illustrated through numerical experiments. Discussions of higher dimensional Zernike polynomials are included in appendices.


title: Performant TCP for Low-Power Wireless Networks

abstract: Low-power and lossy networks (LLNs) enable diverse applications integrating many resource-constrained embedded devices, often requiring interconnectivity with existing TCP/IP networks as part of the Internet of Things. But TCP has received little attention in LLNs due to concerns about its overhead and performance, leading to LLN-specific protocols that require specialized gateways for interoperability. We present a systematic study of a well-designed TCP stack in IEEE 802.15.4-based LLNs, based on the TCP protocol logic in FreeBSD. Through careful implementation and extensive experiments, we show that modern low-power sensor platforms are capable of running full-scale TCP and that TCP, counter to common belief, performs well despite the lossy nature of LLNs. By carefully studying the interaction between the transport and link layers, we identify subtle but important modifications to both, achieving TCP goodput within 25% of an upper bound (5-40x higher than prior results) and low-power operation commensurate to CoAP in a practical LLN application scenario. This suggests that a TCP-based transport layer, seamlessly interoperable with existing TCP/IP networks, is viable and performant in LLNs.


title: Scalable Bottom-up Subspace Clustering using FP-Trees for High Dimensional Data

abstract: Subspace clustering aims to find groups of similar objects (clusters) that exist in lower dimensional subspaces from a high dimensional dataset. It has a wide range of applications, such as analysing high dimensional sensor data or DNA sequences. However, existing algorithms have limitations in finding clusters in non-disjoint subspaces and scaling to large data, which impinge their applicability in areas such as bioinformatics and the Internet of Things. We aim to address such limitations by proposing a subspace clustering algorithm using a bottom-up strategy. Our algorithm first searches for base clusters in low dimensional subspaces. It then forms clusters in higher-dimensional subspaces using these base clusters, which we formulate as a frequent pattern mining problem. This formulation enables efficient search for clusters in higher-dimensional subspaces, which is done using FP-trees. The proposed algorithm is evaluated against traditional bottom-up clustering algorithms and state-of-the-art subspace clustering algorithms. The experimental results show that the proposed algorithm produces clusters with high accuracy, and scales well to large volumes of data. We also demonstrate the algorithm's performance using real-life data, including ten genomic datasets and a car parking occupancy dataset.