title: Comment Ranking Diversification in Forum Discussions

abstract: Viewing consumption of discussion forums with hundreds or more comments depends on ranking because most users only view top-ranked comments. When comments are ranked by an ordered score (e.g. number of replies or up-votes) without adjusting for semantic similarity of near-ranked comments, top-ranked comments are more likely to emphasize the majority opinion and incur redundancy. In this paper, we propose a top K comment diversification re-ranking model using Maximal Marginal Relevance (MMR) and evaluate its impact in three categories: (1) semantic diversity, (2) inclusion of the semantics of lower-ranked comments, and (3) redundancy, within the context of a HarvardX course discussion forum. We conducted a double-blind, small-scale evaluation experiment requiring subjects to select between the top 5 comments of a diversified ranking and a baseline ranking ordered by score. For three subjects, across 100 trials, subjects selected the diversified (75% score, 25% diversification) ranking as significantly (1) more diverse, (2) more inclusive, and (3) less redundant. Within each category, inter-rater reliability showed moderate consistency, with typical Cohen-Kappa scores near 0.2. Our findings suggest that our model improves (1) diversification, (2) inclusion, and (3) redundancy, among top K ranked comments in online discussion forums.


title: Massive Spin-2 Scattering Amplitudes in Extra-Dimensional Theories

abstract: In this paper we describe in detail the computation of the scattering amplitudes of massive spin-2 Kaluza-Klein excitations in a gravitational theory with a single compact extra dimension, whether flat or warped. These scattering amplitudes are characterized by intricate cancellations between different contributions: although individual contributions may grow as fast as ${\cal O}(s^5)$, the full results grow only as ${\cal O}(s)$. We demonstrate that the cancellations persist for all incoming and outgoing particle helicities and examine how truncating the computation to only include a finite number of intermediate states impacts the accuracy of the results. We also carefully assess the range of validity of the low energy effective Kaluza-Klein theory. In particular, for the warped case we demonstrate directly how an emergent low energy scale controls the size of the scattering amplitude, as conjectured by the AdS/CFT correspondence.


title: Fast Join Project Query Evaluation using Matrix Multiplication

abstract: In the last few years, much effort has been devoted to developing join algorithms in order to achieve worst-case optimality for join queries over relational databases. Towards this end, the database community has had considerable success in developing succinct algorithms that achieve worst-case optimal runtime for full join queries, i.e the join is over all variables present in the input database. However, not much is known about join evaluation with {\em projections} beyond some simple techniques of pushing down the projection operator in the query execution plan. Such queries have a large number of applications in entity matching, graph analytics and searching over compressed graphs. In this paper, we study how a class of join queries with projections can be evaluated faster using worst-case optimal algorithms together with matrix multiplication. Crucially, our algorithms are parameterized by the output size of the final result, allowing for choice of the best execution strategy. We implement our algorithms as a subroutine and compare the performance with state-of-the-art techniques to show they can be improved upon by as much as 50x. More importantly, our experiments indicate that matrix multiplication is a useful operation that can help speed up join processing owing to highly optimized open source libraries that are also highly parallelizable.


title: Correlated Feature Selection with Extended Exclusive Group Lasso

abstract: In many high dimensional classification or regression problems set in a biological context, the complete identification of the set of informative features is often as important as predictive accuracy, since this can provide mechanistic insight and conceptual understanding. Lasso and related algorithms have been widely used since their sparse solutions naturally identify a set of informative features. However, Lasso performs erratically when features are correlated. This limits the use of such algorithms in biological problems, where features such as genes often work together in pathways, leading to sets of highly correlated features. In this paper, we examine the performance of a Lasso derivative, the exclusive group Lasso, in this setting. We propose fast algorithms to solve the exclusive group Lasso, and introduce a solution to the case when the underlying group structure is unknown. The solution combines stability selection with random group allocation and introduction of artificial features. Experiments with both synthetic and real-world data highlight the advantages of this proposed methodology over Lasso in comprehensive selection of informative features.


title: Target Detection, Tracking and Avoidance System for Low-cost UAVs using AI-Based Approaches

abstract: An onboard target detection, tracking and avoidance system has been developed in this paper, for low-cost UAV flight controllers using AI-Based approaches. The aim of the proposed system is that an ally UAV can either avoid or track an unexpected enemy UAV with a net to protect itself. In this point of view, a simple and robust target detection, tracking and avoidance system is designed. Two open-source tools were used for the aim: a state-of-the-art object detection technique called SSD and an API for MAVLink compatible systems called MAVSDK. The MAVSDK performs velocity control when a UAV is detected so that the manoeuvre is done simply and efficiently. The proposed system was verified with Software in the loop (SITL) and Hardware in the loop (HITL) simulators. The simplicity of this algorithm makes it innovative, and therefore it should be used in future applications needing robust performances with low-cost hardware such as delivery drone applications.


title: LEEP: A New Measure to Evaluate Transferability of Learned Representations

abstract: We introduce a new measure to evaluate the transferability of representations learned by classifiers. Our measure, the Log Expected Empirical Prediction (LEEP), is simple and easy to compute: when given a classifier trained on a source data set, it only requires running the target data set through this classifier once. We analyze the properties of LEEP theoretically and demonstrate its effectiveness empirically. Our analysis shows that LEEP can predict the performance and convergence speed of both transfer and meta-transfer learning methods, even for small or imbalanced data. Moreover, LEEP outperforms recently proposed transferability measures such as negative conditional entropy and H scores. Notably, when transferring from ImageNet to CIFAR100, LEEP can achieve up to 30% improvement compared to the best competing method in terms of the correlations with actual transfer accuracy.


title: Certified Defense to Image Transformations via Randomized Smoothing

abstract: We extend randomized smoothing to cover parameterized transformations (e.g., rotations, translations) and certify robustness in the parameter space (e.g., rotation angle). This is particularly challenging as interpolation and rounding effects mean that image transformations do not compose, in turn preventing direct certification of the perturbed image (unlike certification with $\ell^p$ norms). We address this challenge by introducing three different kinds of defenses, each with a different guarantee (heuristic, distributional and individual) stemming from the method used to bound the interpolation error. Importantly, we show how individual certificates can be obtained via either statistical error bounds or efficient online inverse computation of the image transformation. We provide an implementation of all methods at https://github.com/eth-sri/transformation-smoothing.


title: "Do the Right Thing" for Whom? An Experiment on Ingroup Favouritism, Group Assorting and Moral Suasion

abstract: In this paper we investigate the effect of moral suasion on ingroup favouritism. We report a well-powered, pre-registered, two-stage 2x2 mixed-design experiment. In the first stage, groups are formed on the basis of how participants answer to a set of questions, concerning non-morally relevant issues in one treatment (assorting on non-moral preferences), and morally relevant issues in another treatment (assorting on moral preferences). In the second stage, participants choose how to split a given amount of money between participants of their own group and participants of the other group, first in the baseline setting and then in a setting where they are told to do what they believe to be morally right (moral suasion). Our main results are: (i) in the baseline, participants tend to favour their own group to a greater extent when groups are assorted according to moral preferences, compared to when they are assorted according to non-moral preferences; (ii) the net effect of moral suasion is to decrease ingroup favouritism, but there is also a non-negligible proportion of participants for whom moral suasion increases ingroup favouritism; (iii) the effect of moral suasion is substantially stable across group assorting and four pre-registered individual characteristics (gender, political orientation, religiosity, pro-life vs pro-choice ethical convictions).


title: Optical projection and spatial separation of spin entangled triplet-pairs from the S1 (21Ag-) state of pi-conjugated systems

abstract: The S1 (21Ag-) state is an optically dark state of natural and synthetic pi-conjugated materials that can play a critical role in optoelectronic processes such as, energy harvesting, photoprotection and singlet fission. Despite this widespread importance, direct experimental characterisations of the electronic structure of the S1 (21Ag-) wavefunction have remained scarce and uncertain, although advanced theory predicts it to have a rich multi-excitonic character. Here, studying an archetypal polymer, polydiacetylene, and carotenoids, we experimentally demonstrate that S1 (21Ag-) is a superposition state with strong contributions from spin-entangled pairs of triplet excitons (1(TT)). We further show that optical manipulation of the S1 (21Ag-) wavefunction using triplet absorption transitions allows selective projection of the 1(TT) component into a manifold of spatially separated triplet-pairs with lifetimes enhanced by up to one order of magnitude and whose yield is strongly dependent on the level of inter-chromophore coupling. Our results provide a unified picture of 21Ag-states in pi-conjugated materials and open new routes to exploit their dynamics in singlet fission, photobiology and for the generation of entangled (spin-1) particles for molecular quantum technologies.

categories: physics.chem-ph cond-mat.mes-hall cond-mat.mtrl-sci physics.optics quant-ph


title: Piecewise linear regressions for approximating distance metrics

abstract: This paper presents a data structure that summarizes distances between configurations across a robot configuration space, using a binary space partition whose cells contain parameters used for a locally linear approximation of the distance function. Querying the data structure is extremely fast, particularly when compared to the graph search required for querying Probabilistic Roadmaps, and memory requirements are promising. The paper explores the use of the data structure constructed for a single robot to provide a heuristic for challenging multi-robot motion planning problems. Potential applications also include the use of remote computation to analyze the space of robot motions, which then might be transmitted on-demand to robots with fewer computational resources.


title: Improving Learning Effectiveness For Object Detection and Classification in Cluttered Backgrounds

abstract: Usually, Neural Networks models are trained with a large dataset of images in homogeneous backgrounds. The issue is that the performance of the network models trained could be significantly degraded in a complex and heterogeneous environment. To mitigate the issue, this paper develops a framework that permits to autonomously generate a training dataset in heterogeneous cluttered backgrounds. It is clear that the learning effectiveness of the proposed framework should be improved in complex and heterogeneous environments, compared with the ones with the typical dataset. In our framework, a state-of-the-art image segmentation technique called DeepLab is used to extract objects of interest from a picture and Chroma-key technique is then used to merge the extracted objects of interest into specific heterogeneous backgrounds. The performance of the proposed framework is investigated through empirical tests and compared with that of the model trained with the COCO dataset. The results show that the proposed framework outperforms the model compared. This implies that the learning effectiveness of the framework developed is superior to the models with the typical dataset.


title: Usual stochastic ordering results for series and parallel systems with components having Exponentiated Chen distribution

abstract: In this paper, we have discussed the usual stochastic ordering relations between two systems. Each system consists of n mutually independent components. The components follow Exponentiated (Extended) Chen distribution with three parameters $\alpha, \beta, \lambda$. Two popular systems are taken into consideration, one is the series system and another is the parallel system. The results in this paper are obtained by varying one parameter and the other parameters are kept constant. The hazard rate ordering or reversed hazard rate ordering relations that are not possible for series or parallel systems have been demonstrated with the help of counterexamples.


title: Spin structure relation to phase contrast imaging of isolated magnetic Bloch and Neel skyrmions

abstract: Magnetic skyrmions are promising candidates for future storage devices with a large data density. A great variety of materials have been found that host skyrmions up to the room-temperature regime. Lorentz microscopy, usually performed in a transmission electron microscope (TEM), is one of the most important tools for characterizing skyrmion samples in real space. Using numerical calculations, this work relates the phase contrast in a TEM to the actual magnetization profile of an isolated Neel or Bloch skyrmion, the two most common skyrmion types. Within the framework of the used skyrmion model, the results are independent of skyrmion size and wall width and scale with sample thickness for purely magnetic specimens. Simple rules are provided to extract the actual skyrmion configuration of pure Bloch or Neel skyrmions without the need of simulations. Furthermore, first differential phase contrast (DPC) measurements on Neel skyrmions that meet experimental expectations are presented and showcase the described principles. The work is relevant for material sciences where it enables the engineering of skyrmion profiles via convenient characterization.


title: RSANet: Recurrent Slice-wise Attention Network for Multiple Sclerosis Lesion Segmentation

abstract: Brain lesion volume measured on T2 weighted MRI images is a clinically important disease marker in multiple sclerosis (MS). Manual delineation of MS lesions is a time-consuming and highly operator-dependent task, which is influenced by lesion size, shape and conspicuity. Recently, automated lesion segmentation algorithms based on deep neural networks have been developed with promising results. In this paper, we propose a novel recurrent slice-wise attention network (RSANet), which models 3D MRI images as sequences of slices and captures long-range dependencies through a recurrent manner to utilize contextual information of MS lesions. Experiments on a dataset with 43 patients show that the proposed method outperforms the state-of-the-art approaches. Our implementation is available online at https://github.com/tinymilky/RSANet.


title: A list-mode OSEM-based attenuation and scatter compensation method for SPECT

abstract: Reliable attenuation and scatter compensation (ASC) is a prerequisite for quantification and beneficial for visual interpretation tasks in SPECT. In this paper, we develop a reconstruction method that uses the entire SPECT emission data, i.e. data in both the photopeak and scatter windows, acquired in list-mode format and including the energy attribute of the detected photon, to perform ASC. We implemented a GPU-based version of this method using an ordered subsets expectation maximization (OSEM) algorithm. The method was objectively evaluated using realistic simulation studies on the task of estimating uptake in the striatal regions of the brain in a 2-D dopamine transporter (DaT)-scan SPECT study. We observed that inclusion of data from the scatter window and using list-mode data yielded improved quantification compared to using data only from the photopeak window or using binned data. These results motivate further development of list-mode-based ASC methods that include scatter-window data for SPECT.


title: Global stabilization and destabilization by the state dependent noise with particular distributions

abstract: Under natural assumptions, an unstable equilibrium of a difference equation can be stabilized by a bounded multiplicative noise, identically distributed at each step. This includes stabilization of an otherwise unstable positive equilibrium of Ricker, logistic, and Beverton-Holt maps. Introduction of a multiplicative noise also allows to destabilize a stable equilibrium in a sense that all solutions stay away from this point, almost surely. In our examples a noise has symmetric, discrete or continuous, distribution with support $[-1,1]$, including Bernoulli and uniform continuous distribution. We obtain conditions on the noise amplitudes in each case that allow to either stabilize or destabilize an equilibrium. Computer simulations illustrate our results.


title: Beyond the Trees: Resilient Multipath for Last-mile WISP Networks

abstract: Expanding the reach of the Internet is a topic of widespread interest today. Google and Facebook, among others, have begun investing substantial research efforts toward expanding Internet access at the edge. Compared to data center networks, which are relatively over-engineered, last-mile networks are highly constrained and end up being ultimately responsible for the performance issues that impact the user experience. The most viable and cost-effective approach for providing last-mile connectivity has proved to be Wireless ISPs (WISPs), which rely on point-to-point wireless backhaul infrastructure to provide connectivity using cheap commodity wireless hardware. However, individual WISP network links are known to have poor reliability and the networks as a whole are highly cost constrained. Motivated by these observations, we propose Wireless ISPs with Redundancy (WISPR), which leverages the cost-performance tradeoff inherent in commodity wireless hardware to move toward a greater number of inexpensive links in WISP networks thereby lowering costs. To take advantage of this new path diversity, we introduce a new, general protocol that provides increased performance, reliability, or a combination of the two.


title: A stochastic comparison study for the smallest and largest ordered statistic from Weibull-G and Gompertz Makeham distribution

abstract: In this paper, we have discussed the stochastic comparison of the smallest and largest ordered statistic from independent heterogeneous Weibull-G random variables and Gompertz Makeham random variables. We compare systems arising from taking different model parameters and obtain stochastic ordering results under the condition of multivariate chain majorization. Using the notion of vector majorization, we compare different systems and obtain stochastic ordering results.


title: Cautious Reinforcement Learning via Distributional Risk in the Dual Domain

abstract: We study the estimation of risk-sensitive policies in reinforcement learning problems defined by a Markov Decision Process (MDPs) whose state and action spaces are countably finite. Prior efforts are predominately afflicted by computational challenges associated with the fact that risk-sensitive MDPs are time-inconsistent. To ameliorate this issue, we propose a new definition of risk, which we call caution, as a penalty function added to the dual objective of the linear programming (LP) formulation of reinforcement learning. The caution measures the distributional risk of a policy, which is a function of the policy's long-term state occupancy distribution. To solve this problem in an online model-free manner, we propose a stochastic variant of primal-dual method that uses Kullback-Lieber (KL) divergence as its proximal term. We establish that the number of iterations/samples required to attain approximately optimal solutions of this scheme matches tight dependencies on the cardinality of the state and action spaces, but differs in its dependence on the infinity norm of the gradient of the risk measure. Experiments demonstrate the merits of this approach for improving the reliability of reward accumulation without additional computational burdens.

categories: stat.ML cs.AI cs.LG cs.SY eess.SY math.OC


title: Probing quantum processor performance with pyGSTi

abstract: PyGSTi is a Python software package for assessing and characterizing the performance of quantum computing processors. It can be used as a standalone application, or as a library, to perform a wide variety of quantum characterization, verification, and validation (QCVV) protocols on as-built quantum processors. We outline pyGSTi's structure, and what it can do, using multiple examples. We cover its main characterization protocols with end-to-end implementations. These include gate set tomography, randomized benchmarking on one or many qubits, and several specialized techniques. We also discuss and demonstrate how power users can customize pyGSTi and leverage its components to create specialized QCVV protocols and solve user-specific problems.


title: ISR and IFI in Precision AFB Studies with KKMC-hh

abstract: KKMC-hh is a hadronic event generator for Z boson production and decays, which includes exponentiated multi-photon radiation and first-order electroweak corrections. We have used KKMC-hh to investigate the role of initial sate radiation (ISR) and initial-final interference (IFI) in precision electroweak analyses at the LHC. We compare the effect of this radiation on angular distributions and forward-backward asymmetry, which are particularly important for the measurement of the weak mixing angle. We discuss the relation of the ISR implementation in KKMC-hh to ISR from parton distribution functions with QED corrections.


title: Time Series Data Augmentation for Deep Learning: A Survey

abstract: Deep learning performs remarkably well on many time series analysis tasks recently. The superior performance of deep neural networks relies heavily on a large number of training data to avoid overfitting. However, the labeled data of many real-world time series applications may be limited such as classification in medical time series and anomaly detection in AIOps. As an effective way to enhance the size and quality of the training data, data augmentation is crucial to the successful application of deep learning models on time series data. In this paper, we systematically review different data augmentation methods for time series. We propose a taxonomy for the reviewed methods, and then provide a structured review for these methods by highlighting their strengths and limitations. We also empirically compare different data augmentation methods for different tasks including time series classification, anomaly detection, and forecasting. Finally, we discuss and highlight five future directions to provide useful research guidance.


title: Stochastic Pore Collapse Models in Granular Materials

abstract: Stochastic models for pore collapse in granular materials are developed. First, a general fluctuating stress-strain relation for a plastic flow rule is derived. The fluctuations account for non-associativity in plastic deformations typically observed in heterogeneous materials. Second, an axisymmetric spherical shell compaction model is extended to account for fluctuations in the material microstructure due to granular interactions at the pore scale. This changes the stress-strain constitutive equation determining the dynamics of pore collapse. Results show that stochastic differential equations can account for multiscale interactions in a statistical sense.


title: A new approach to the generalized Springer correspondence

abstract: The Springer resolution of the nilpotent cone is used to give a geometric construction of the irreducible representations of Weyl groups. Borho and MacPherson obtain the Springer correspondence by applying the decomposition theorem to the Springer resolution, establishing an injective map from the set of irreducible Weyl group representations to simple equivariant perverse sheaves on the nilpotent cone. In this manuscript, we consider a generalization of the Springer resolution using a variety defined by the first author. Our main result shows that in the type A case, applying the decomposition theorem to this map yields all simple perverse sheaves on the nilpotent cone with multiplicity as predicted by Lusztig's generalized Springer correspondence.


title: The UTMOST pulsar timing programme II: Timing noise across the pulsar population

abstract: While pulsars possess exceptional rotational stability, large scale timing studies have revealed at least two distinct types of irregularities in their rotation: red timing noise and glitches. Using modern Bayesian techniques, we investigated the timing noise properties of 300 bright southern-sky radio pulsars that have been observed over 1.0-4.8 years by the upgraded Molonglo Observatory Synthesis Telescope (MOST). We reanalysed the spin and spin-down changes associated with nine previously reported pulsar glitches, report the discovery of three new glitches and four unusual glitch-like events in the rotational evolution of PSR J1825$-$0935. We develop a refined Bayesian framework for determining how red noise strength scales with pulsar spin frequency ($\nu$) and spin-down frequency ($\dot{\nu}$), which we apply to a sample of 280 non-recycled pulsars. With this new method and a simple power-law scaling relation, we show that red noise strength scales across the non-recycled pulsar population as $\nu^{a} |\dot{\nu}|^{b}$, where $a = -0.84^{+0.47}_{-0.49}$ and $b = 0.97^{+0.16}_{-0.19}$. This method can be easily adapted to utilise more complex, astrophysically motivated red noise models. Lastly, we highlight our timing of the double neutron star PSR J0737$-$3039, and the rediscovery of a bright radio pulsar originally found during the first Molonglo pulsar surveys with an incorrectly catalogued position.


title: An Improved Generic ER Schema for Conceptual Modeling of Information Systems

abstract: The Entity-Relationship (ER) model is widely used for creating ER schemas for modeling application domains in the field of Information Systems development. However, when an ER schema is transformed to a Relational Database Schema (RDS), some important information on the ER schema may not be represented meaningfully on the RDS. This causes a loss of information during the transformation process. Although, several previous researches have proposed solutions to remedy the situation, the problem still exists. Thus, in this on-going research, we wish to improve the proposed solutions and maximize information preservation in the ER to relational transformation process. Cardinality ratio constraints, role names, composite attributes, and certain relationship types are among the information frequently lost in the transformation process. Deficiencies in the ER model and the transformation method seems to cause this situation. We take the view that if the information lost is resolved; a one-to-one mapping should exist from the ER schema to its RDS. We modified the ER model and the transformation algorithm following a heuristic research method with a view to eliminating the deficiencies and thereby achieving a one-to-one mapping. We should show that the mapping exists for any real-world application. We create a generic ER schema - an ER schema that represents any phenomena in symbolic form - and use it to show that a one-to-one mapping exists for any real-world application. In this paper, we explore our generic ER schema and its advantages over its predecessors in view of representing any real-world application.


title: Atmospheric Circulation on Black Widow Companions

abstract: We present a model for atmospheric wind circulation in binary millisecond pulsar (MSP) companions, showing how the optical light curve (LC) and radial velocities (RV) are sensitive to the wind flow, causing LC orbital phase shifts and asymmetries, as observed for several `spider' MSP. Velocity widths of spectral lines offer additional opportunities for measuring surface wind speed. As examples, we fit optical data for a black widow pulsar J1959+2048 and a redback pulsar J2215+5135; the wind heating models (WH) are statistically strongly preferred over direct heating (DH) for both objects, although the latter is even better fit with a heated spot. In general, WH effects tend to increase the inferred orbital inclination $i$ and decrease the inferred companion center-of-mass radial velocity amplitude $K_\mathrm{c}$; both effects decrease the inferred neutron star mass. Even with such a decrease, we find large masses for the two neutron stars: $2.18\pm0.09M_\odot$ and $2.28^{+0.10}_{-0.09}M_\odot$, respectively (for the modest surface speeds fit from the bulk heat flow; supersonic photospheric winds can slightly change these values). These are among the highest masses known, and our improved modeling increases confidence that the results are important for understanding the dense matter equation of state.


title: The Simulation, Fabrication Technology and Characteristic Research of Micro-Pressure Sensor with Isosceles Trapezoidal Beam-Membrane

abstract: A micro-pressure sensor with an isosceles trapezoidal beam-membrane (ITBM) is proposed in this paper, consisting of a square silicon membrane, four isosceles trapezoidal beams and four piezoresistors.To investigate how the elastic silicon membrane affects pressure sensitive characteristics, a simulation models based on ANSYS 15.0 software were used to analyze the effect of structural dimension on characteristics of pressure sensor. According to that, the chips of micro-pressure sensors were fabricated by micro-electro-mechanical system (MEMS) technology on a silicon wafer with <100> orientation.The experimental results show that the proposed sensor achieves a better sensitivity of 9.64 mV/kPa and an excellent linearity of 0.09%F.S. in the range of 0~3.0 kPa at room temperature and a supply voltage of 5.0 V,with a super temperature coefficient of sensitivity(TCS) about - 684 ppm/K from 235.15 K to 360.15 K and low pressure measurement less than 3.0 kPa.


title: Generalized Self-Adapting Particle Swarm Optimization algorithm with archive of samples

abstract: In this paper we enhance Generalized Self-Adapting Particle Swarm Optimization algorithm (GAPSO), initially introduced at the Parallel Problem Solving from Nature 2018 conference, and to investigate its properties. The research on GAPSO is underlined by the two following assumptions: (1) it is possible to achieve good performance of an optimization algorithm through utilization of all of the gathered samples, (2) the best performance can be accomplished by means of a combination of specialized sampling behaviors (Particle Swarm Optimization, Differential Evolution, and locally fitted square functions). From a software engineering point of view, GAPSO considers a standard Particle Swarm Optimization algorithm as an ideal starting point for creating a generalpurpose global optimization framework. Within this framework hybrid optimization algorithms are developed, and various additional techniques (like algorithm restart management or adaptation schemes) are tested. The paper introduces a new version of the algorithm, abbreviated as M-GAPSO. In comparison with the original GAPSO formulation it includes the following four features: a global restart management scheme, samples gathering within an R-Tree based index (archive/memory of samples), adaptation of a sampling behavior based on a global particle performance, and a specific approach to local search. The above-mentioned enhancements resulted in improved performance of M-GAPSO over GAPSO, observed on both COCO BBOB testbed and in the black-box optimization competition BBComp. Also, for lower dimensionality functions (up to 5D) results of M-GAPSO are better or comparable to the state-of-the art version of CMA-ES (namely the KL-BIPOP-CMA-ES algorithm presented at the GECCO 2017 conference).


title: Distributionally Robust Chance Constrained Programming with Generative Adversarial Networks (GANs)

abstract: This paper presents a novel deep learning based data-driven optimization method. A novel generative adversarial network (GAN) based data-driven distributionally robust chance constrained programming framework is proposed. GAN is applied to fully extract distributional information from historical data in a nonparametric and unsupervised way without a priori approximation or assumption. Since GAN utilizes deep neural networks, complicated data distributions and modes can be learned, and it can model uncertainty efficiently and accurately. Distributionally robust chance constrained programming takes into consideration ambiguous probability distributions of uncertain parameters. To tackle the computational challenges, sample average approximation method is adopted, and the required data samples are generated by GAN in an end-to-end way through the differentiable networks. The proposed framework is then applied to supply chain optimization under demand uncertainty. The applicability of the proposed approach is illustrated through a county-level case study of a spatially explicit biofuel supply chain in Illinois.


title: Accelerating and Stopping Resistance Drift in Phase Change Memory Cells via High Electric Field Stress

abstract: We observed resistance drift in 125 K - 300 K temperature range in melt quenched amorphous Ge2Sb2Te5 line-cells with length x width x thickness = ~500 nm x ~100 nm x ~ 50 nm. Drift coefficients measured using small voltage sweeps appear to decrease from 0.12 +/- 0.029 at 300 K to 0.075 +/- 0.006 at 125 K. The current-voltage characteristics of the amorphized cells measured in the 85 K - 300 K using high-voltage sweeps (0 to ~25 V) show a combination of a linear, low-field exponential and high-field exponential conduction mechanisms, all of which are strong functions of temperature. The very first high-voltage sweep after amorphization (with electric fields up to ~70% of the breakdown field) shows clear hysteresis in the current-voltage characteristics due to accelerated drift, while the consecutive sweeps show stable characteristics. Stabilization was achieved with 50 nA compliance current (current densities ~104 A/cm^2), preventing appreciable self-heating in the cells. The observed acceleration and stoppage of the resistance drift with the application of high electric fields is attributed to changes in the electrostatic potential profile within amorphous Ge2Sb2Te5 due to trapped charges, reducing tunneling current. Stable current-voltage characteristics are used to extract carrier activation energies for the conduction mechanisms in 85 K - 300 K temperature range. The carrier activation energy associated with linear current-voltage response is extracted to be 331 +/- 5 meV in 200 - 300 K range, while carrier activation energies of 233 +/- 2 meV and 109 +/- 5 meV are extracted in 85 K to 300 K range for the mechanisms that give exponential current-voltage responses.


title: The almost optimal multilinear restriction estimate for hypersurfaces with curvature: the case of $n-1$ hypersurfaces in R^n

abstract: In this paper we establish the optimal multilinear restriction estimate for n-1 hypersurfaces with some curvature, where $n$ is the dimension of the underlying space. The result is sharp up to the endpoint and the role of curvature is made precise in terms of the shape operator.


title: Cross-modality Person re-identification with Shared-Specific Feature Transfer

abstract: Cross-modality person re-identification (cm-ReID) is a challenging but key technology for intelligent video analysis. Existing works mainly focus on learning common representation by embedding different modalities into a same feature space. However, only learning the common characteristics means great information loss, lowering the upper bound of feature distinctiveness. In this paper, we tackle the above limitation by proposing a novel cross-modality shared-specific feature transfer algorithm (termed cm-SSFT) to explore the potential of both the modality-shared information and the modality-specific characteristics to boost the re-identification performance. We model the affinities of different modality samples according to the shared features and then transfer both shared and specific features among and across modalities. We also propose a complementary feature learning strategy including modality adaption, project adversarial learning and reconstruction enhancement to learn discriminative and complementary shared and specific features of each modality, respectively. The entire cm-SSFT algorithm can be trained in an end-to-end manner. We conducted comprehensive experiments to validate the superiority of the overall algorithm and the effectiveness of each component. The proposed algorithm significantly outperforms state-of-the-arts by 22.5% and 19.3% mAP on the two mainstream benchmark datasets SYSU-MM01 and RegDB, respectively.


title: Resonances as Viscosity Limits for Exterior Dilation Analytic Potentials

abstract: For exterior dilation analytic potential, $V$, we use the method of complex scaling to show that the resonances of $ - \Delta + V $, in a conic neighbourhood of the real axis, are limits of eigenvalues of $ - \Delta + V - i \epsilon x^2 $ as $ \epsilon \to 0+ $, if $V$ can be analytically extended from $\mathbb{R}^n$ to a truncated cone in $\mathbb{C}^n$.


title: Qualitative properties for solutions to conformally invariant fourth order critical systems

abstract: We study qualitative properties for nonnegative solutions to a conformally invariant coupled system of fourth order equations involving critical exponents. For solutions defined in the punctured space, there exist essentially two cases to analyze. If the origin is a removable singularity, we prove that non-singular solutions are rotationally invariant and weakly positive. More precisely, they are the product of a fourth order spherical solution by a unit vector with nonnegative coordinates. If the origin is a non-removable singularity, we show that the solutions are radially symmetric and strongly positive. Furthermore, using a Pohozaev-type invariant, we prove the non-existence of semi-singular solutions, that is, all components equally blow-up in the neighborhood of origin. Namely, they are classified as multiples of the Emden--Fowler solution. Our results are natural generalizations of the famous classification due to [L. A. Caffarelli, B. Gidas and J. Spruck, Comm. Pure Appl. Math. (1989)] on the classical singular Yamabe equation.


title: Road Curb Detection and Localization with Monocular Forward-view Vehicle Camera

abstract: We propose a robust method for estimating road curb 3D parameters (size, location, orientation) using a calibrated monocular camera equipped with a fisheye lens. Automatic curb detection and localization is particularly important in the context of Advanced Driver Assistance System (ADAS), i.e. to prevent possible collision and damage of the vehicle's bumper during perpendicular and diagonal parking maneuvers. Combining 3D geometric reasoning with advanced vision-based detection methods, our approach is able to estimate the vehicle to curb distance in real time with mean accuracy of more than 90%, as well as its orientation, height and depth. Our approach consists of two distinct components - curb detection in each individual video frame and temporal analysis. The first part comprises of sophisticated curb edges extraction and parametrized 3D curb template fitting. Using a few assumptions regarding the real world geometry, we can thus retrieve the curb's height and its relative position w.r.t. the moving vehicle on which the camera is mounted. Support Vector Machine (SVM) classifier fed with Histograms of Oriented Gradients (HOG) is used for appearance-based filtering out outliers. In the second part, the detected curb regions are tracked in the temporal domain, so as to perform a second pass of false positives rejection. We have validated our approach on a newly collected database of 11 videos under different conditions. We have used point-wise LIDAR measurements and manual exhaustive labels as a ground truth.


title: Optimization with Momentum: Dynamical, Control-Theoretic, and Symplectic Perspectives

abstract: We analyze the convergence rate of various momentum-based optimization algorithms from a dynamical systems point of view. Our analysis exploits fundamental topological properties, such as the continuous dependence of iterates on their initial conditions, to provide a simple characterization of convergence rates. In many cases, closed-form expressions are obtained that relate algorithm parameters to the convergence rate. The analysis encompasses discrete time and continuous time, as well as time-invariant and time-variant formulations, and is not limited to a convex or Euclidean setting. In addition, the article rigorously establishes why symplectic discretization schemes are important for momentum-based optimization algorithms, and provides a characterization of algorithms that exhibit accelerated convergence.


title: Rationally Inattentive Path-Planning via RRT*

abstract: We consider a path-planning scenario for a mobile robot traveling in a configuration space with obstacles under the presence of stochastic disturbances. A novel path length metric is proposed on the uncertain configuration space and then integrated with the existing RRT* algorithm. The metric is a weighted sum of two terms which capture both the Euclidean distance traveled by the robot and the perception cost, i.e., the amount of information the robot must perceive about the environment to follow the path safely. The continuity of the path length function with respect to the topology of the total variation metric is shown and the optimality of the Rationally Inattentive RRT* algorithm is discussed. Three numerical studies are presented which display the utility of the new algorithm.


title: Spectral convergence in geometric quantization --- the case of toric symplectic manifolds

abstract: In this paper, we show the spectral convergence result of $\overline{\partial}$-Laplacians when $(X,\omega)$ is a compact toric symplectic manifold equipped with the natural prequantum line bundle $L$. We consider a family $\{ J_s\}_s$ of $\omega$-compatible complex structures tending to the large complex structure limit, and obtain the spectral convergence of $\overline{\partial}$-Laplacians acting on $L^k$.


title: Double-Differential Inclusive Charged-Current $\nu_\mu$ Cross Sections on Hydrocarbon in MINERvA at $\langle E_{\nu} \rangle \sim$ 3.5 GeV

abstract: MINERvA reports inclusive charged-current cross sections for muon neutrinos on hydrocarbon in the NuMI beamline. We measured the double-differential cross section in terms of the longitudinal and transverse muon momenta, as well as the single-differential cross sections in those variables. The data used in this analysis correspond to an exposure of $3.34 \times 10^{20}$ protons on target with a peak neutrino energy of approximately 3.5 GeV. Measurements are compared to the GENIE, NuWro and GiBUU neutrino cross-section predictions, as well as a version of GENIE modified to produce better agreement with prior exclusive MINERvA measurements. None of the models or variants were able to successfully reproduce the data across the entire phase space, which includes areas dominated by each interaction channel.


title: Digital quantum simulation of molecular dynamics and control

abstract: Optimally-shaped electromagnetic fields have the capacity to coherently control the dynamics of quantum systems and thus offer a promising means for controlling molecular transformations relevant to chemical, biological, and materials applications. Currently, advances in this area are hindered by the prohibitive cost of the quantum dynamics simulations needed to explore the principles and possibilities of molecular control. However, the emergence of nascent quantum-computing devices suggests that efficient simulations of quantum dynamics may be on the horizon. In this article, we study how quantum computers could be employed to design optimally-shaped fields to control molecular systems. We introduce a hybrid algorithm that utilizes a quantum computer for simulating the field-induced quantum dynamics of a molecular system in polynomial time, in combination with a classical optimization approach for updating the field. Qubit encoding methods relevant for molecular control problems are described, and procedures for simulating the quantum dynamics and obtaining the simulation results are discussed. Numerical illustrations are then presented that explicitly treat paradigmatic vibrational and rotational control problems, and also consider how optimally-shaped fields could be used to elucidate the mechanisms of energy transfer in light-harvesting complexes. Resource estimates, as well as a numerical assessment of the impact of hardware noise and the prospects of near-term hardware implementations, are provided for the latter task.


title: Lie Biderivations on Triangular Algebras

abstract: Let $\mathcal{T}$ be a triangular algebra over a commutative ring $\mathcal{R}$ and $\varphi: \mathcal{T} \times \mathcal{T}\longrightarrow \mathcal{T}$ be an arbitrary Lie biderivation of $\mathcal{T}$. We will address the question of describing the form of $\varphi$ in the current work. It is shown that under certain mild assumptions, $\varphi$ is the sum of an inner biderivation and an extremal biderivation and a some central bilinear mapping. Our results is immediately applied to block upper triangular algebras and Hilbert space nest algebras .


title: On Catastrophic Interference in Atari 2600 Games

abstract: Model-free deep reinforcement learning is sample inefficient. One hypothesis -- speculated, but not confirmed -- is that catastrophic interference within an environment inhibits learning. We test this hypothesis through a large-scale empirical study in the Arcade Learning Environment (ALE) and, indeed, find supporting evidence. We show that interference causes performance to plateau; the network cannot train on segments beyond the plateau without degrading the policy used to reach there. By synthetically controlling for interference, we demonstrate performance boosts across architectures, learning algorithms and environments. A more refined analysis shows that learning one segment of a game often increases prediction errors elsewhere. Our study provides a clear empirical link between catastrophic interference and sample efficiency in reinforcement learning.


title: Efficiently Guiding Imitation Learning Agents with Human Gaze

abstract: Human gaze is known to be an intention-revealing signal in human demonstrations of tasks. In this work, we use gaze cues from human demonstrators to enhance the performance of agents trained via three popular imitation learning methods -- behavioral cloning (BC), behavioral cloning from observation (BCO), and Trajectory-ranked Reward EXtrapolation (T-REX). Based on similarities between the attention of reinforcement learning agents and human gaze, we propose a novel approach for utilizing gaze data in a computationally efficient manner, as part of an auxiliary loss function, which guides a network to have higher activations in image regions where the human's gaze fixated. This work is a step towards augmenting any existing convolutional imitation learning agent's training with auxiliary gaze data. Our auxiliary coverage-based gaze loss (CGL) guides learning toward a better reward function or policy, without adding any additional learnable parameters and without requiring gaze data at test time. We find that our proposed approach improves the performance by 95% for BC, 343% for BCO, and 390% for T-REX, averaged over 20 different Atari games. We also find that compared to a prior state-of-the-art imitation learning method assisted by human gaze (AGIL), our method achieves better performance, and is more efficient in terms of learning with fewer demonstrations. We further interpret trained CGL agents with a saliency map visualization method to explain their performance. At last, we show that CGL can help alleviate a well-known causal confusion problem in imitation learning.


title: Learning Multivariate Hawkes Processes at Scale

abstract: Multivariate Hawkes Processes (MHPs) are an important class of temporal point processes that have enabled key advances in understanding and predicting social information systems. However, due to their complex modeling of temporal dependencies, MHPs have proven to be notoriously difficult to scale, what has limited their applications to relatively small domains. In this work, we propose a novel model and computational approach to overcome this important limitation. By exploiting a characteristic sparsity pattern in real-world diffusion processes, we show that our approach allows to compute the exact likelihood and gradients of an MHP -- independently of the ambient dimensions of the underlying network. We show on synthetic and real-world datasets that our model does not only achieve state-of-the-art predictive results, but also improves runtime performance by multiple orders of magnitude compared to standard methods on sparse event sequences. In combination with easily interpretable latent variables and influence structures, this allows us to analyze diffusion processes at previously unattainable scale.


title: Spontaneous Pulse Formation in Edge-Less Photonic Crystal Resonators

abstract: Complex systems are a proving ground for fundamental interactions between components and their collective emergent phenomena. Through intricate design, integrated photonics offers intriguing nonlinear interactions that create new patterns of light. In particular, the canonical Kerr-nonlinear resonator becomes unstable with a sufficiently intense traveling-wave excitation, yielding instead a Turing pattern composed of a few interfering waves. These resonators also support the localized soliton pulse as a separate nonlinear stationary state. Kerr solitons are remarkably versatile for applications, but they cannot emerge from constant excitation. Here, we explore an edge-less photonic-crystal resonator (PhCR) that enables spontaneous formation of a soliton pulse in place of the Turing pattern. We design a PhCR in the regime of single-azimuthal-mode engineering to re-balance Kerr-nonlinear frequency shifts in favor of the soliton state, commensurate with how group-velocity dispersion balances nonlinearity. Our experiments establish PhCR solitons as mode-locked pulses by way of ultraprecise optical-frequency measurements, and we characterize their fundamental properties. Our work shows that sub-wavelength nanophotonic design expands the palette for nonlinear engineering of light.


title: Impact of epitaxial strain on the topological-nontopological phase diagram and semimetallic behavior of InAs/GaSb composite quantum wells

abstract: We study the influence of epitaxial strain on the electronic properties of InAs/GaSb composite quantum wells (CQWs), host structures for quantum spin Hall insulators, by transport measurements and eight-band $\mathbf{k\cdot p}$ calculations. Using different substrates and buffer layer structures for crystal growth, we prepare two types of samples with vastly different strain conditions. CQWs with a nearly strain-free GaSb layer exhibit a resistance peak at the charge neutrality point that reflects the opening of a topological gap in the band-inverted regime. In contrast, for CQWs with 0.50\% biaxial tensile strain in the GaSb layer, semimetallic behavior indicating a gap closure is found for the same degree of band inversion. Additionally, with the tensile strain, the boundary between the topological and nontopological regimes is located at a larger InAs thickness. Eight-band $\mathbf{k\cdot p}$ calculations reveal that tensile strain in GaSb not only shifts the phase boundary but also significantly modifies the band structure, which can result in the closure of an indirect gap and make the system semimetallic even in the topological regime. Our results thus provide a global picture of the topological-nontopological phase diagram as a function of layer thicknesses and strain.


title: Detecting Patch Adversarial Attacks with Image Residuals

abstract: We introduce an adversarial sample detection algorithm based on image residuals, specifically designed to guard against patch-based attacks. The image residual is obtained as the difference between an input image and a denoised version of it, and a discriminator is trained to distinguish between clean and adversarial samples. More precisely, we use a wavelet domain algorithm for denoising images and demonstrate that the obtained residuals act as a digital fingerprint for adversarial attacks. To emulate the limitations of a physical adversary, we evaluate the performance of our approach against localized (patch-based) adversarial attacks, including in settings where the adversary has complete knowledge about the detection scheme. Our results show that the proposed detection method generalizes to previously unseen, stronger attacks and that it is able to reduce the success rate (conversely, increase the computational effort) of an adaptive attacker.


title: Measurement of the Sixth-Order Cumulant of Net-Proton Distributions in Au+Au Collisions from the STAR Experiment

abstract: Higher-order cumulants of the net-baryon multiplicity distributions are predicted to be sensitive to the properties of the nuclear matter created in high-energy nuclear collisions. In this talk, we present the collision centrality and acceptance (rapidity and transverse momentum) dependence of the ratio of the 6$^{{\rm th}}$- to the 2$^{\rm nd}$-order cumulant ratio ($C_{6}/C_{2}$) of net-proton in Au+Au collisions at $\sqrt{s_{\rm NN}}=54.4$ and 200 GeV measured by the STAR detector at RHIC. The new results are compared to hadron transport model and lattice QCD calculations.


title: Forensic analysis of the Windows telemetry for diagnostics

abstract: Telemetry is the automated sensing and collection of data from a remote device. It is often used to provide better services for users. Microsoft uses telemetry to periodically collect information about Windows systems and to help improve user experience and fix potential issues. Windows telemetry service functions by creating RBS files on the local system to reliably transfer and manage the telemetry data, and these files can provide useful information in a digital forensic investigation. Combined with the information derived from traditional Windows forensics, investigators can have greater confidence in the evidence derived from various artifacts. It is possible to acquire information that can be confirmed only for live systems, such as the computer hardware serial number, the connection records for external storage devices, and traces of executed processes. This information is included in the RBS files that are created for use in Windows telemetry. In this paper, we introduced how to acquire RBS files telemetry and analyzed the data structure of these RBS files, which are able to determine the types of information that Windows OS have been collected. We also discussed the reliability and the novelty by comparing the conventional artifacts with the RBS files, which could be useful in digital forensic investigation.


title: Decentralized Federated Learning via SGD over Wireless D2D Networks

abstract: Federated Learning (FL), an emerging paradigm for fast intelligent acquisition at the network edge, enables joint training of a machine learning model over distributed data sets and computing resources with limited disclosure of local data. Communication is a critical enabler of large-scale FL due to significant amount of model information exchanged among edge devices. In this paper, we consider a network of wireless devices sharing a common fading wireless channel for the deployment of FL. Each device holds a generally distinct training set, and communication typically takes place in a Device-to-Device (D2D) manner. In the ideal case in which all devices within communication range can communicate simultaneously and noiselessly, a standard protocol that is guaranteed to converge to an optimal solution of the global empirical risk minimization problem under convexity and connectivity assumptions is Decentralized Stochastic Gradient Descent (DSGD). DSGD integrates local SGD steps with periodic consensus averages that require communication between neighboring devices. In this paper, wireless protocols are proposed that implement DSGD by accounting for the presence of path loss, fading, blockages, and mutual interference. The proposed protocols are based on graph coloring for scheduling and on both digital and analog transmission strategies at the physical layer, with the latter leveraging over-the-air computing via sparsity-based recovery.


title: Near-optimal ground state preparation

abstract: Preparing the ground state of a given Hamiltonian and estimating its ground energy are important but computationally hard tasks. However, given some additional information, these problems can be solved efficiently on a quantum computer. We assume that an initial state with non-trivial overlap with the ground state can be efficiently prepared, and the spectral gap between the ground energy and the first excited energy is bounded from below. With these assumptions we design an algorithm that prepares the ground state when an upper bound of the ground energy is known, whose runtime has a logarithmic dependence on the inverse error. When such an upper bound is not known, we propose a hybrid quantum-classical algorithm to estimate the ground energy, where the dependence of the number of queries to the initial state on the desired precision is exponentially improved compared to the current state-of-the-art algorithm proposed in [Ge et al. 2019]. These two algorithms can then be combined to prepare a ground state without knowing an upper bound of the ground energy. We also prove that our algorithms reach the complexity lower bounds by applying it to the unstructured search problem and the quantum approximate counting problem.


title: DGST : Discriminator Guided Scene Text detector

abstract: Scene text detection task has attracted considerable attention in computer vision because of its wide application. In recent years, many researchers have introduced methods of semantic segmentation into the task of scene text detection, and achieved promising results. This paper proposes a detector framework based on the conditional generative adversarial networks to improve the segmentation effect of scene text detection, called DGST (Discriminator Guided Scene Text detector). Instead of binary text score maps generated by some existing semantic segmentation based methods, we generate a multi-scale soft text score map with more information to represent the text position more reasonably, and solve the problem of text pixel adhesion in the process of text extraction. Experiments on standard datasets demonstrate that the proposed DGST brings noticeable gain and outperforms state-of-the-art methods. Specifically, it achieves an F-measure of 87% on ICDAR 2015 dataset.


title: The Complexity of Possible Winners on Partial Chains

abstract: The Possible Winner (PW) problem, a fundamental algorithmic problem in computational social choice, concerns elections where voters express only partial preferences between candidates. Via a sequence of investigations, a complete classification of the complexity of the PW problem was established for all pure positional scoring rules: the PW problem is in P for the plurality and veto rules, and NP-complete for all other such rules. More recently, the PW problem was studied on classes of restricted partial orders that arise in natural settings, such as partitioned partial orders and truncated partial orders; in particular, it was shown that there are rules for which the PW problem drops from NP-complete to P on such restricted partial orders. Here, we investigate the PW problem on partial chains, i.e., partial orders that are a total order on a subset of their domains. Such orders arise naturally in a variety of settings, including rankings of movies or restaurants. We classify the complexity of the PW problem on partial chains by establishing that, perhaps surprisingly, this restriction does not change the complexity of the problem, namely, the PW problem is NP-complete for all pure positional scoring rules other than the plurality and veto rules. As a byproduct, we obtain a new and more principled proof of the complexity of the PW problem on arbitrary partial orders.


title: Localization with Deep Neural Networks using mmWave Ray Tracing Simulations

abstract: The world is moving towards faster data transformation with more efficient localization of a user being the preliminary requirement. This work investigates the use of a deep learning technique for wireless localization, considering both millimeter-wave (mmWave) and sub-6 GHz frequencies. The capability of learning a new neural network model makes the localization process easier and faster. In this study, a Deep Neural Network (DNN) was used to localize User Equipment (UE) in two static scenarios. We propose two different methods to train a neural network, one using channel parameters (features) and another using a channel response vector and compare their performances using preliminary computer simulations. We observe that the former approach produces high localization accuracy considering that all of the users have a fixed number of multipath components (MPCs), this method is reliant on the number of MPCs. On the other hand, the latter approach is independent of the MPCs, but it performs relatively poorly compared to the first approach.


title: Distortion of Magnetic Fields in BHR 71

abstract: The magnetic field structure of a star-forming Bok globule BHR 71 was determined based on near-infrared polarimetric observations of background stars. The magnetic field in BHR 71 was mapped from 25 stars. By using a simple 2D parabolic function, the plane-of-sky magnetic axis of the core was found to be $\theta_{\rm mag} = 125^{\circ} \pm 11^{\circ}$. The plane-of-sky mean magnetic field strength of BHR 71 was found to be $B_{\rm pos} = 8.8 - 15.0$ $\mu$G, indicating that the BHR 71 core is magnetically supercritical with $\lambda = 1.44 - 2.43$. Taking into account the effect of thermal/turbulent pressure and the plane-of-sky magnetic field component, the critical mass of BHR 71 was $M_{\rm cr} = 14.5-18.7$ M$_{\odot}$, which is consistent with the observed core mass of $M_{\rm core} \approx 14.7$ M$_{\odot}$ (Yang et al. 2017). We conclude that BHR 71 is in a condition close to a kinematically critical state, and the magnetic field direction lies close to the plane of sky. Since BHR 71 is a star-forming core, a significantly subcritical condition (i.e., the magnetic field direction deviating from the plane of sky) is unlikely, and collapsed from a condition close to a kinematically critical state. There are two possible scenarios to explain the curved magnetic fields of BHR 71, one is an hourglass-like field structure due to mass accumulation and the other is the Inoue \& Fukui (2013) mechanism, which proposes the interaction of the core with a shock wave to create curved magnetic fields wrapping around the core.


title: Proper Motions of the Radio Source Orion MR, Formerly Known as Orion n, and New Sources with Large Proper Motions in Orion BN/KL

abstract: The infrared source known as Orion n was detected in 1980 with observations made with the 3.8-m United Kingdom Infrared Telescope. About two decades later, sensitive observations made with the Very Large Array revealed the presence of a mJy double radio source apparently coincident in position with the infrared source n. The radio source was assumed to be the counterpart of the infrared source. However, over the years it has been concluded that the radio source shows large proper motions to the south while the infrared source n is stationary. Here we reanalyze the proper motions of the radio source adding both older and newer VLA observations than previously used. We confirm the proper motions of the radio source that at present no longer coincides positionally with the infrared source. The solution to this problem is, most probably, that the infrared source n and the radio source are not the same object: the infrared source is a stationary object in the region while the radio counterpart is moving as a result of the explosion that took place in this region some 500 years ago and that expelled large amounts of molecular gas as well as several compact sources. Considering the paper where it was first reported, we refer to this double radio source as Orion MR. In addition, we use these new observations to fully confirm the large proper motions of the sources IRc23 and Zapata 11. Together with sources BN, I, Orion MR, and x, there are at least six compact sources that recede from a point in common in Orion BN/KL. However, IRc23 is peculiar in that its ejection age appears to be only $\sim$300 years. The relatively large number of sources rules out as a possible mechanism the classic three-body scenario since then only two escaping bodies are expected: a tight binary plus the third star involved in the encounter.


title: Critical length for the spreading-vanishing dichotomy in higher dimensions

abstract: We consider an extension of the classical Fisher-Kolmogorov equation, called the \textit{Fisher-Stefan} model, which is a moving boundary problem on $0 < x < L(t)$. A key property of the Fisher-Stefan model is the \textit{spreading-vanishing dichotomy}, where solutions with $L(t) > L_{\textrm{c}}$ will eventually spread as $t \to \infty$, whereas solutions where $L(t) \ngtr L_{\textrm{c}}$ will vanish as $t \to \infty$. In one dimension is it well-known that the critical length is $L_{\textrm{c}} = \pi/2$. In this work we re-formulate the Fisher-Stefan model in higher dimensions and calculate $L_{\textrm{c}}$ as a function of spatial dimensions in a radially symmetric coordinate system. Our results show how $L_{\textrm{c}}$ depends upon the dimension of the problem and numerical solutions of the governing partial differential equation are consistent with our calculations.


title: A Prospective Randomized Clinical Trial for Measuring Radiology Study Reporting Time on Artificial Intelligence-Based Detection of Intracranial Hemorrhage in Emergent Care Head CT

abstract: We propose Artificial Intelligence Prospective Randomized Observer Blinding Evaluation (AI-PROBE) for quantitative clinical performance evaluation of radiology AI systems within prospective randomized clinical trials. AI-PROBE encompasses a study design and a matching radiology IT infrastructure that randomly blinds radiologists for results provided by AI-based image analysis. To demonstrate the applicability of our evaluation framework, we present a first prospective randomized clinical trial on the effect of Intra-Cranial Hemorrhage (ICH) detection in emergent care head CT on radiology study Turn-Around Time (TAT). Here, we acquired 620 non-contrast head CT scans from inpatient and emergency room patients at a large academic hospital. Following acquisition, scans were automatically analyzed for the presence of ICH using commercially available software (Aidoc, Tel Aviv, Israel). Cases identified positive for ICH by AI (ICH-AI+) were flagged in radiologists' reading worklists, where flagging was randomly switched off with probability 50%. TAT was measured as time difference between study completion and first clinically communicated reporting, with time stamps automatically retrieved from various IT systems. TATs for flagged cases (73+/-143 min) were significantly lower than TATs for non-flagged (132+/-193 min) cases (p<0.05, one-sided t-test), where 105 of 122 ICH-AI+ cases were true positive. Total sensitivity, specificity, and accuracy over all analyzed cases were 95.0%, 96.7%, and 96.4%, respectively. We conclude that automatic identification of ICH reduces TAT for ICH in emergent care head CT, which carries the potential for improving timely clinical management of ICH. Our results suggest that AI-PROBE can contribute to systematic quantitative evaluation of AI systems in clinical practice using clinically meaningful quantities, such as TAT or diagnostic accuracy.


title: Bringing Inter-Thread Cache Benefits to Federated Scheduling -- Extended Results & Technical Report

abstract: Multiprocessor scheduling of hard real-time tasks modeled by directed acyclic graphs (DAGs) exploits the inherent parallelism presented by the model. For DAG tasks, a node represents a request to execute an object on one of the available processors. In one DAG task, there may be multiple execution requests for one object, each represented by a distinct node. These distinct execution requests offer an opportunity to reduce their combined cache overhead through coordinated scheduling of objects as threads within a parallel task. The goal of this work is to realize this opportunity by incorporating the cache-aware BUNDLE-scheduling algorithm into federated scheduling of sporadic DAG task sets. This is the first work to incorporate instruction cache sharing into federated scheduling. The result is a modification of the DAG model named the DAG with objects and threads (DAG-OT). Under the DAG-OT model, descriptions of nodes explicitly include their underlying executable object and number of threads. When possible, nodes assigned the same executable object are collapsed into a single node; joining their threads when BUNDLE-scheduled. Compared to the DAG model, the DAG-OT model with cache-aware scheduling reduces the number of cores allocated to individual tasks by approximately 20 percent in the synthetic evaluation and up to 50 percent on a novel parallel computing platform implementation. By reducing the number of allocated cores, the DAG-OT model is able to schedule a subset of previously infeasible task sets.


title: Time-Resolved and Energy-Resolved Polarizations of GRB Prompt Emission

abstract: Besides light curves and spectra, polarization provides a different powerful tool of studying the $\gamma-$ray burst (GRB) prompt phase. Compared with the time-integrated and energy-integrated polarization, time-resolved and energy-resolved polarization can deliver more physical information about the emitting region. Here we investigate time-resolved and energy-resolved polarization of GRB prompt emission using the synchrotron models. We find that the equal arrival time surface effect is very important in shaping the PD curves when the physical conditions of emitting region changes violently with radius. Polarization properties are neither correlated with the spectral lag nor the peak energy evolution patterns. Polarization properties with a mixed magnetic field are very similar to those for a corresponding ordered magnetic field but the former has a smaller polarization degree. The emission at the MeV peak can be highly polarized for a synchrotron model while it is unpolarized as predicted by a dissipative photosphere model. Future energy-resolved polarization observations can distinguish between these two models.


title: Multistage Distributionally Robust Mixed-Integer Programming with Decision-Dependent Moment-Based Ambiguity Sets

abstract: We study multistage distributionally robust mixed-integer programs under endogenous uncertainty, where the probability distribution of stage-wise uncertainty depends on the decisions made in previous stages. We first consider two ambiguity sets defined by decision-dependent bounds on the first and second moments of uncertain parameters and by mean and covariance matrix that exactly match decision-dependent empirical ones, respectively. For both sets, we show that the subproblem in each stage can be recast as a mixed-integer linear program (MILP). Moreover, we extend the general moment-based ambiguity set in (Delage and Ye, 2010) to the multistage decision-dependent setting, and derive mixed-integer semidefinite programming (MISDP) reformulations of stage-wise subproblems. We develop methods for attaining lower and upper bounds of the optimal objective value of the multistage MISDPs, and approximate them using a series of MILPs. We deploy the Stochastic Dual Dynamic integer Programming (SDDiP) method for solving the problem under the three ambiguity sets with risk-neutral or risk-averse objective functions, and conduct numerical studies on multistage facility-location instances having diverse sizes under different parameter and uncertainty settings. Our results show that the SDDiP quickly finds optimal solutions for moderate-sized instances under the first two ambiguity sets, and also finds good approximate bounds for the multistage MISDPs derived under the third ambiguity set. We also demonstrate the efficacy of incorporating decision-dependent distributional ambiguity in multistage decision-making processes.


title: Flavor physics and Jarlskog determinant

abstract: The flavor problem is reviewed starting with the chiral symmetry, and the $A_4$ symmetry derivation and its realization in GUTs are presented.


title: Utilizing Network Properties to Detect Erroneous Inputs

abstract: Neural networks are vulnerable to a wide range of erroneous inputs such as adversarial, corrupted, out-of-distribution, and misclassified examples. In this work, we train a linear SVM classifier to detect these four types of erroneous data using hidden and softmax feature vectors of pre-trained neural networks. Our results indicate that these faulty data types generally exhibit linearly separable activation properties from correct examples, giving us the ability to reject bad inputs with no extra training or overhead. We experimentally validate our findings across a diverse range of datasets, domains, pre-trained models, and adversarial attacks.


title: Improved Image Coding Autoencoder With Deep Learning

abstract: In this paper, we build autoencoder based pipelines for extreme end-to-end image compression based on Ball\'e's approach, which is the state-of-the-art open source implementation in image compression using deep learning. We deepened the network by adding one more hidden layer before each strided convolutional layer with exactly the same number of down-samplings and up-samplings. Our approach outperformed Ball\'e's approach, and achieved around 4.0% reduction in bits per pixel (bpp), 0.03% increase in multi-scale structural similarity (MS-SSIM), and only 0.47% decrease in peak signal-to-noise ratio (PSNR), It also outperforms all traditional image compression methods including JPEG2000 and HEIC by at least 20% in terms of compression efficiency at similar reconstruction image quality. Regarding encoding and decoding time, our approach takes similar amount of time compared with traditional methods with the support of GPU, which means it's almost ready for industrial applications.


title: Sylvester rank functions for amenable normal extensions

abstract: We introduce a notion of amenable normal extension S of a unital ring R with a finite approximation system F, encompassing the amenable algebras over a field of Gromov and Elek, the twisted crossed product by an amenable group, and the tensor product with a field extension. It is shown that every Sylvester matrix rank function rk of R preserved by S has a canonical extension to a Sylvester matrix rank function rk_F for S. In the case of twisted crossed product by an amenable group, and the tensor product with a field extension, it is also shown that rk_F depends on rk continuously. When an amenable group has a twisted action on a unital C*-algebra preserving a tracial state, we also show that two natural Sylvester matrix rank functions on the algebraic twisted crossed product constructed out of the tracial state coincide.


title: Well-Defined Quantum Soliton Masses Without Supersymmetry

abstract: 22 years ago, Rebhan and van Nieuwenhuizen showed that loop corrections to the mass of a quantum soliton depend on a choice of matching condition for the regulators of the vacuum and one-soliton sector Hamiltonians. In supersymmetric theories, regulators which preserve supersymmetry yield the correct quantum corrections, as these are dictated by supersymmetry. However, in a general theory it is not known which matching condition yields the correct mass. We use the leading term in the operator that creates the soliton to construct the regulated one-soliton sector Hamiltonian from that of the vacuum sector, providing the correct matching condition. As an application, we derive a simple formula for the one-loop mass of a kink in a large class of 1+1 dimensional scalar field theories and also, at one loop, we diagonalize the Hamiltonian which describes the kink excitations.


title: Enhanced quantum oscillations in Kondo insulators

abstract: Quantum oscillations have long been regarded as the manifestation of the Fermi surface in metals. However, they were recently observed in Kondo insulators. We examine the Kondo screening due to Landau levels in Kondo insulators. It is shown that even for large Kondo insulating gaps, appreciable amplitudes of quantum oscillations that are consistent with experimental observations are present both in magnetization and resistivity. Specifically, we show that due to the periodic alignment between the Landau levels in the conduction and the f-orbit electrons, the Kondo screening itself undergoes oscillations so that the electronic structure oscillates with the magnetic field. Our results explain main features of quantum oscillations observed in experiments. They indicate that the non-rigidity of the electronic structure results in observable quantum oscillations Kondo insulators. This new effect provides a new way to probe the Fermi surface geometry of insulators.


title: An ALMA Survey of H$_2$CO in Protoplanetary Disks

abstract: H$_2$CO is one of the most abundant organic molecules in protoplanetary disks and can serve as a precursor to more complex organic chemistry. We present an ALMA survey of H$_2$CO towards 15 disks covering a range of stellar spectral types, stellar ages, and dust continuum morphologies. H$_2$CO is detected towards 13 disks and tentatively detected towards a 14th. We find both centrally-peaked and centrally-depressed emission morphologies, and half of the disks show ring-like structures at or beyond expected CO snowline locations. Together these morphologies suggest that H$_2$CO in disks is commonly produced through both gas-phase and CO-ice-regulated grain-surface chemistry. We extract disk-averaged and azimuthally-averaged H$_2$CO excitation temperatures and column densities for four disks with multiple H$_2$CO line detections. The temperatures are between 20-50K, with the exception of colder temperatures in the DM Tau disk. These temperatures suggest that H$_2$CO emission in disks is generally emerging from the warm molecular layer, with some contributions from the colder midplane. Applying the same H$_2$CO excitation temperatures to all disks in the survey, we find that H$_2$CO column densities span almost three orders of magnitude ($\sim 5 \times 10^{11} - 5 \times 10^{14} \mathrm{cm}^{-2}$). The column densities appear uncorrelated with disk size and stellar age, but Herbig Ae disks may have less H$_2$CO compared to T Tauri disks, possibly because of less CO freeze-out. More H$_2$CO observations towards Herbig Ae disks are needed to confirm this tentative trend, and to better constrain under which disk conditions H$_2$CO and other oxygen-bearing organics efficiently form during planet formation.


title: Future global stability for relativistic perfect fluids with linear equations of state $p=K\rho$ where $1/3<K<1/2$

abstract: We establish the future stability of nonlinear perturbations of a class of homogeneous solutions to the relativistic Euler equations with a linear equation of state $p=K\rho$ on exponentially expanding FLRW spacetimes for the equation of state parameter values $1/3 < K < 1/2$.


title: Are L2 adversarial examples intrinsically different?

abstract: Deep Neural Network (DDN) has achieved notable success in various tasks, including many security concerning scenarios. However, a considerable amount of work has proved its vulnerability to adversaries. We unravel the properties that can intrinsically differentiate adversarial examples and normal inputs through theoretical analysis. That is, adversarial examples generated by $L_2$ attacks usually have larger input sensitivity which can be used to identify them efficiently. We also found that those generated by $L_\infty$ attacks will be different enough in the pixel domain to be detected empirically. To verify our analysis, we proposed a \textbf{G}uided \textbf{C}omplementary \textbf{D}efense module (\textbf{GCD}) integrating detection and recovery processes. When compared with adversarial detection methods, our detector achieves a detection AUC of over 0.98 against most of the attacks. When comparing our guided rectifier with commonly used adversarial training methods and other rectification methods, our rectifier outperforms them by a large margin. We achieve a recovered classification accuracy of up to 99\% on MNIST, 89\% on CIFAR-10, and 87\% on ImageNet subsets against $L_2$ attacks. Furthermore, under the white-box setting, our holistic defensive module shows a promising degree of robustness. Thus, we confirm that at least $L_2$ adversarial examples are intrinsically different enough from normal inputs both theoretically and empirically. And we shed light upon designing simple yet effective defensive methods with these properties.


title: Handling Position Bias for Unbiased Learning to Rank in Hotels Search

abstract: Nowadays, search ranking and recommendation systems rely on a lot of data to train machine learning models such as Learning-to-Rank (LTR) models to rank results for a given query, and implicit user feedbacks (e.g. click data) have become the dominant source of data collection due to its abundance and low cost, especially for major Internet companies. However, a drawback of this data collection approach is the data could be highly biased, and one of the most significant biases is the position bias, where users are biased towards clicking on higher ranked results. In this work, we will investigate the marginal importance of properly handling the position bias in an online test environment in Tripadvisor Hotels search. We propose an empirically effective method of handling the position bias that fully leverages the user action data. We take advantage of the fact that when user clicks a result, he has almost certainly observed all the results above, and the propensities of the results below the clicked result will be estimated by a simple but effective position bias model. The online A/B test results show that this method leads to an improved search ranking model.


title: Remarks on the speeds of a class of random walks on the integers

abstract: In recent years, there has been an interest in deriving certain important probabilistic results as consequences of deterministic ones; see for instance \cite{beig} and \cite{acc}. In this work, we continue on this path by deducing a well known equivalence between the speed of random walks on the integers and the growth of the size of their ranges. This result is an immediate consequence of the Kesten-Spitzer-Whitman theorem, and by appearances is probabilistic in nature, but we will show that it follows easily from an elementary deterministic result. We also investigate the common property of recurrent random walks of having speed zero, and show by example that this property need not be shared by deterministic sequences. However, if we consider the inter-arrival times (times at which the sequence is equal to 0) then we find a sufficient deterministic condition for a sequence to have zero speed, and show that this can be used to derive several probabilistic results.


title: Temporal Convolutional Attention-based Network For Sequence Modeling

abstract: With the development of feed-forward models, the default model for sequence modeling has gradually evolved to replace recurrent networks. Many powerful feed-forward models based on convolutional networks and attention mechanism were proposed and show more potential to handle sequence modeling tasks. We wonder that is there an architecture that can not only achieve an approximate substitution of recurrent network, but also absorb the advantages of feed-forward models. So we propose an exploratory architecture referred to Temporal Convolutional Attention-based Network (TCAN) which combines temporal convolutional network and attention mechanism. TCAN includes two parts, one is Temporal Attention (TA) which captures relevant features inside the sequence, the other is Enhanced Residual (ER) which extracts shallow layer's important information and transfers to deep layers. We improve the state-of-the-art results of bpc/perplexity to 30.28 on word-level PTB, 1.092 on character-level PTB, and 9.20 on WikiText-2.


title: Recent developments in the PySCF program package

abstract: PYSCF is a Python-based general-purpose electronic structure platform that both supports first-principles simulations of molecules and solids, as well as accelerates the development of new methodology and complex computational workflows. The present paper explains the design and philosophy behind PYSCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PYSCF as a development environment. We then summarize the capabilities of PYSCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PYSCF across the domains of quantum chemistry, materials science, machine learning and quantum information science.


title: Fermi-surface topology and renormalization of bare ellipticity in an interacting anisotropic electron gas

abstract: We investigate effects of electron-electron interactions on the shape of the Fermi surface in an anisotropic two-dimensional electron gas using the `RPA-GW' self-energy approximation. We find that the interacting Fermi surface deviates from an ellipse, but not in an arbitrary way. The interacting Fermi surface has only two qualitatively distinct shapes for most values of $r_s$. The Fermi surface undergoes two distinct transitions between these two shapes as $r_s$ increases. For larger $r_s$, the degree of the deviation from an ellipse rapidly increases, but, in general, our theory provides a justification for the widely used elliptical Fermi surface approximation even for the interacting system since the non-elliptic corrections are quantitatively rather small except for very large $r_s$.


title: First Constraints on Compact Dark Matter from Fast Radio Burst Microstructure

abstract: Despite existing constraints, it remains possible that up to $35\%$ of all dark matter is comprised of compact objects, such as the black holes in the 10-100\,M$_\odot$ range whose existence has been confirmed by LIGO. The strong gravitational lensing of transients such as FRBs and GRBs has been suggested as a more sensitive probe for compact dark matter than intensity fluctuations observed in microlensing experiments. Recently ASKAP has reported burst substructure down to $15\mu$s timescales in FRBs in the redshift range $0.3-0.5$. We investigate here the implications of this for the detectability of compact dark matter by FRBs. We find that a sample size of $\sim130$ FRBs would be required to constrain compact dark matter to less than the existing 35$\%$ limit with 95$\%$ confidence, if it were distributed along $\gtrsim 1\,$Gpc-long FRB sightlines through the cosmic web. Conversely, existing constraints on the fraction of compact dark matter permit as many as 1 in $\approx 40$ of all $z \lesssim 0.4$ FRBs to exhibit micro-lensed burst structure. Approximately $170$ FRBs intercepting halos within $\sim 50\,$kpc would be required to exclude the fraction of compact dark matter in each intercepted halo to a similar level. Furthermore, we consider the cumulative effects of lensing of the FRB signal by a macroscopic dark matter distribution. We conclude that lensing from a uniform distribution of compact objects is likely not observable, but suggest that FRBs may set meaningful limits on power-law distributions of dark matter.


title: Current status and muon $g-2$ explanation of lepton portal dark matter

abstract: In this paper, we summarize phenomenology in lepton portal dark matter (DM) models, where DM couples to leptons and extra leptons/sleptons. There are several possible setups: complex/real scalar DM and Dirac/Majorana fermion DM. In addition, there are choices for the lepton chirality that couples to DM. We discuss the prediction of each model and compare it with the latest experimental constraints from the DM, the LHC, and the flavor experiments. We also propose a simple setup to achieve the discrepancy in the anomalous magnetic moment of muon.


title: A Video Analysis Method on Wanfang Dataset via Deep Neural Network

abstract: The topic of object detection has been largely improved recently, especially with the development of convolutional neural network. However, there still exist a lot of challenging cases, such as small object, compact and dense or highly overlapping object. Existing methods can detect multiple objects wonderfully, but because of the slight changes between frames, the detection effect of the model will become unstable, the detection results may result in dropping or increasing the object. In the pedestrian flow detection task, such phenomenon can not accurately calculate the flow. To solve this problem, in this paper, we describe the new function for real-time multi-object detection in sports competition and pedestrians flow detection in public based on deep learning. Our work is to extract a video clip and solve this frame of clips efficiently. More specfically, our algorithm includes two stages: judge method and optimization method. The judge can set a maximum threshold for better results under the model, the threshold value corresponds to the upper limit of the algorithm with better detection results. The optimization method to solve detection jitter problem. Because of the occurrence of frame hopping in the video, and it will result in the generation of video fragments discontinuity. We use optimization algorithm to get the key value, and then the detection result value of index is replaced by key value to stabilize the change of detection result sequence. Based on the proposed algorithm, we adopt wanfang sports competition dataset as the main test dataset and our own test dataset for YOLOv3-Abnormal Number Version(YOLOv3-ANV), which is 5.4% average improvement compared with existing methods. Also, video above the threshold value can be obtained for further analysis. Spontaneously, our work also can used for pedestrians flow detection and pedestrian alarm tasks.


title: Automated classification of stems and leaves of potted plants based on point cloud data

abstract: The accurate classification of plant organs is a key step in monitoring the growing status and physiology of plants. A classification method was proposed to classify the leaves and stems of potted plants automatically based on the point cloud data of the plants, which is a nondestructive acquisition. The leaf point training samples were automatically extracted by using the three-dimensional convex hull algorithm, while stem point training samples were extracted by using the point density of a two-dimensional projection. The two training sets were used to classify all the points into leaf points and stem points by utilizing the support vector machine (SVM) algorithm. The proposed method was tested by using the point cloud data of three potted plants and compared with two other methods, which showed that the proposed method can classify leaf and stem points accurately and efficiently.


title: Generalized Sliced Distances for Probability Distributions

abstract: Probability metrics have become an indispensable part of modern statistics and machine learning, and they play a quintessential role in various applications, including statistical hypothesis testing and generative modeling. However, in a practical setting, the convergence behavior of the algorithms built upon these distances have not been well established, except for a few specific cases. In this paper, we introduce a broad family of probability metrics, coined as Generalized Sliced Probability Metrics (GSPMs), that are deeply rooted in the generalized Radon transform. We first verify that GSPMs are metrics. Then, we identify a subset of GSPMs that are equivalent to maximum mean discrepancy (MMD) with novel positive definite kernels, which come with a unique geometric interpretation. Finally, by exploiting this connection, we consider GSPM-based gradient flows for generative modeling applications and show that under mild assumptions, the gradient flow converges to the global optimum. We illustrate the utility of our approach on both real and synthetic problems.


title: Explainable $k$-Means and $k$-Medians Clustering

abstract: Clustering is a popular form of unsupervised learning for geometric data. Unfortunately, many clustering algorithms lead to cluster assignments that are hard to explain, partially because they depend on all the features of the data in a complicated way. To improve interpretability, we consider using a small decision tree to partition a data set into clusters, so that clusters can be characterized in a straightforward manner. We study this problem from a theoretical viewpoint, measuring cluster quality by the $k$-means and $k$-medians objectives: Must there exist a tree-induced clustering whose cost is comparable to that of the best unconstrained clustering, and if so, how can it be found? In terms of negative results, we show, first, that popular top-down decision tree algorithms may lead to clusterings with arbitrarily large cost, and second, that any tree-induced clustering must in general incur an $\Omega(\log k)$ approximation factor compared to the optimal clustering. On the positive side, we design an efficient algorithm that produces explainable clusters using a tree with $k$ leaves. For two means/medians, we show that a single threshold cut suffices to achieve a constant factor approximation, and we give nearly-matching lower bounds. For general $k \geq 2$, our algorithm is an $O(k)$ approximation to the optimal $k$-medians and an $O(k^2)$ approximation to the optimal $k$-means. Prior to our work, no algorithms were known with provable guarantees independent of dimension and input size.


title: Spectropolarimetric analysis of FRB 181112 at microsecond resolution: Implications for Fast Radio Burst emission mechanism

abstract: We have developed a new coherent dedispersion mode to study the emission of Fast Radio Bursts that trigger the voltage capture capability of the Australian SKA Pathfinder (ASKAP) interferometer. In principle the mode can probe emission timescales down to 3 ns with full polarimetric information preserved. Enabled by the new capability, here we present a spectropolarimetric analysis of FRB 181112 detected by ASKAP, localized to a galaxy at redshift 0.47. At microsecond time resolution the burst is resolved into four narrow pulses with a rise time of just $15 \mu$s for the brightest. The pulses have a diversity of morphology, but do not show evidence for temporal broadening by turbulent plasma along the line of sight, nor is there any evidence for periodicity in their arrival times. The pulses are highly polarized (up to 95%), with the polarization position angle varying both between and within pulses. The pulses have apparent rotation measures that vary by $15\pm 2\, {\rm rad \,m^{-2}}$ and apparent dispersion measures that vary by $0.041\pm 0.004\,{\rm pc\,cm^{-3}}$. Conversion between linear and circular polarization is observed across the brightest pulse. We conclude that the FRB 181112 pulses are most consistent with being a direct manifestation of the emission process or the result of propagation through a relativistic plasma close to the source. This demonstrates that our method, which facilitates high-time-resolution polarimetric observations of FRBs, can be used to study not only burst emission processes, but also a diversity of propagation effects present on the gigaparsec paths they traverse.


title: UKARA 1.0 Challenge Track 1: Automatic Short-Answer Scoring in Bahasa Indonesia

abstract: We describe our third-place solution to the UKARA 1.0 challenge on automated essay scoring. The task consists of a binary classification problem on two datasets | answers from two different questions. We ended up using two different models for the two datasets. For task A, we applied a random forest algorithm on features extracted using unigram with latent semantic analysis (LSA). On the other hand, for task B, we only used logistic regression on TF-IDF features. Our model results in F1 score of 0.812.