title: Non-homogeneous wave equation on a cone

abstract: The wave equation $\left(\partial_{tt} - c^2 \Delta_x\right) u(x,t) = e^{-t} f(x,t)$ is shown to have a unique solution if $u$ and its partial derivatives in $x$ are in $L^2(e^{-t})$ on the cone, and the solution can be explicit given in the Fourier series of orthogonal polynomials on the cone. This provides a particular solution for the boundary value problems of the non-homogeneous wave equation on the cone, which can be combined with a solution to the homogeneous wave equation in the cone to obtain the full solution.


title: Scalar Fields and the FLRW Singularity

abstract: The dynamics of multiple scalar fields on a flat FLRW spacetime can be described entirely as a relational system in terms of the matter alone. The matter dynamics is an autonomous system from which the geometrical dynamics can be inferred, and this autonomous system remains deterministic at the point corresponding to the singularity of the cosmology. We show the continuation of this system corresponds to a parity inversion at the singularity, and that the singularity itself is a surface on which the space-time manifold becomes non-orientable.


title: Exact Recovery of Tensor Robust Principal Component Analysis under Linear Transforms

abstract: This work studies the Tensor Robust Principal Component Analysis (TRPCA) problem, which aims to exactly recover the low-rank and sparse components from their sum. Our model is motivated by the recently proposed linear transforms based tensor-tensor product and tensor SVD. We define a new transforms depended tensor rank and the corresponding tensor nuclear norm. Then we solve the TRPCA problem by convex optimization whose objective is a weighted combination of the new tensor nuclear norm and the $\ell_1$-norm. In theory, we show that under certain incoherence conditions, the convex program exactly recovers the underlying low-rank and sparse components. It is of great interest that our new TRPCA model generalizes existing works. In particular, if the studied tensor reduces to a matrix, our TRPCA model reduces to the known matrix RPCA. Our new TRPCA which is allowed to use general linear transforms can be regarded as an extension of our former TRPCA work which uses the discrete Fourier transform. But their proof of the recovery guarantee is different. Numerical experiments verify our results and the application on image recovery demonstrates the superiority of our method.


title: Data-Driven Wide-Area Control Design of Power System Using the Passivity Shortage Framework

abstract: A novel wide-area control design is presented to mitigate inter-area power frequency oscillations. A large-scale power system is decomposed into a network of passivity-short subsystems whose nonlinear interconnections have a state-dependent affine form, and by utilizing the passivity shortage framework a two-level design procedure is developed. At the lower level, any generator control can be viewed as one that makes the generator passivity-short and $L_2$ stable, and the stability impact of the lower-level control on the overall system can be characterized in terms of two parameters. While the system is nonlinear, the impact parameters can be optimized by solving a data-driven matrix inequality (DMI), and the high-level wide-area control is then designed by solving another Lyapunov matrix inequality in terms of the design parameters. The proposed methodology makes the design modular, and the resulting control is adaptive with respect to operating conditions of the power system. A test system is used to illustrate the proposed design, including DMI and the wide-area control, and simulation results demonstrate effectiveness in damping out inter-area oscillations.


title: Evidence of non collisional femtosecond laser electron heating in dielectric materials

abstract: Electron dynamics in the bulk of large band gap dielectric crystals induced by intense femtosecond laser pulses at 800 nm is studied. With laser intensities under the ablation threshold (a few 10 TW/cm\textsuperscript{2}), electrons with unexpected energies in excess of 40-50 eV are observed by using the photoemission spectroscopy. A theoretical approach based on the Boltzmann kinetic equation including state-of-the-art modeling for various particles interactions is developed to interpret these experimental observations. A direct comparison shows that both electron heating in the bulk and a further laser field acceleration after ejection from the material contribute equivalently to the final electron energy gain. The electron heating in the bulk is shown to be significantly driven by a non-collisional process, i.e. direct multiphoton transitions between sub-bands of the conduction band. This work also sheds light on the contribution of the standard electron excitation/relaxation collisional processes, providing a new baseline to study the electron dynamics in dielectric materials and associated applications as laser material structuring.


title: Non-diffusive angular momentum transport in rotating $\bf z$-pinches

abstract: The stability of conducting Taylor-Couette flows under the presence of toroidal magnetic background fields is considered. For strong enough magnetic amplitudes such magnetohydrodynamic flows are unstable against nonaxisymmetric perturbations which may also transport angular momentum. In accordance with the often used diffusion approximation one expects the angular momentum transport vanishing for rigid rotation. In the sense of a nondiffusive $\Lambda$ effect, however, even for {\em rigidly} rotating $z$-pinches an axisymmetric angular momentum flux appears which is directed outward (inward) for large (small) magnetic Mach numbers. The internal rotation in a magnetized rotating tank can thus never be uniform. Those particular rotation laws are used to estimate the value of the instability-induced eddy viscosity for which the nondiffusive $\Lambda$ effect and the diffusive shear-induced transport compensate each other. The results provide the Shakura-Sunyaev viscosity ansatz leading to numerical values linearly growing with the Reynolds number of rotation.


title: Compositional Deep Learning

abstract: Neural networks have become an increasingly popular tool for solving many real-world problems. They are a general framework for differentiable optimization which includes many other machine learning approaches as special cases. In this thesis we build a category-theoretic formalism around a class of neural networks exemplified by CycleGAN. CycleGAN is a collection of neural networks, closed under composition, whose inductive bias is increased by enforcing composition invariants, i.e. cycle-consistencies. Inspired by Functorial Data Migration, we specify the interconnection of these networks using a categorical schema, and network instances as set-valued functors on this schema. We also frame neural network architectures, datasets, models, and a number of other concepts in a categorical setting and thus show a special class of functors, rather than functions, can be learned using gradient descent. We use the category-theoretic framework to conceive a novel neural network architecture whose goal is to learn the task of object insertion and object deletion in images with unpaired data. We test the architecture on three different datasets and obtain promising results.


title: Investigating Target Set Reduction for End-to-End Speech Recognition of Hindi-English Code-Switching Data

abstract: End-to-end (E2E) systems are fast replacing the conventional systems in the domain of automatic speech recognition. As the target labels are learned directly from speech data, the E2E systems need a bigger corpus for effective training. In the context of code-switching task, the E2E systems face two challenges: (i) the expansion of the target set due to multiple languages involved, and (ii) the lack of availability of sufficiently large domain-specific corpus. Towards addressing those challenges, we propose an approach for reducing the number of target labels for reliable training of the E2E systems on limited data. The efficacy of the proposed approach has been demonstrated on two prominent architectures, namely CTC-based and attention-based E2E networks. The experimental validations are performed on a recently created Hindi-English code-switching corpus. For contrast purpose, the results for the full target set based E2E system and a hybrid DNN-HMM system are also reported.


title: DNN-based Speaker Embedding Using Subjective Inter-speaker Similarity for Multi-speaker Modeling in Speech Synthesis

abstract: This paper proposes novel algorithms for speaker embedding using subjective inter-speaker similarity based on deep neural networks (DNNs). Although conventional DNN-based speaker embedding such as a $d$-vector can be applied to multi-speaker modeling in speech synthesis, it does not correlate with the subjective inter-speaker similarity and is not necessarily appropriate speaker representation for open speakers whose speech utterances are not included in the training data. We propose two training algorithms for DNN-based speaker embedding model using an inter-speaker similarity matrix obtained by large-scale subjective scoring. One is based on similarity vector embedding and trains the model to predict a vector of the similarity matrix as speaker representation. The other is based on similarity matrix embedding and trains the model to minimize the squared Frobenius norm between the similarity matrix and the Gram matrix of $d$-vectors, i.e., the inter-speaker similarity derived from the $d$-vectors. We crowdsourced the inter-speaker similarity scores of 153 Japanese female speakers, and the experimental results demonstrate that our algorithms learn speaker embedding that is highly correlated with the subjective similarity. We also apply the proposed speaker embedding to multi-speaker modeling in DNN-based speech synthesis and reveal that the proposed similarity vector embedding improves synthetic speech quality for open speakers whose speech utterances are unseen during the training.


title: High-pressure hybrid materials that can store hydrogen in table salt

abstract: We demonstrate in this paper that high pressure can promote the reactions between the ionic compounds and H2 molecules and form thermodynamically stable hybrid compounds. Using crystal structure search method based on first principles calculations and particle swarm optimization algorithm, we show that many alkali halides XY (X=Na, K, Rb, Cs; Y=Cl, Br, I) can form stable hybrid compounds with various H2 compositions under different pressures that could be as low as 0.8 GPa (KIH2). Especially, the commonly known table salt, NaCl, can form stable NaClH2 and NaCl(H2)4 hybrid compounds under pressures of 19.7 and 37.9 GPa, respectively. Our results demonstrate that pressure can promote the formation of solid-molecule hybrid compounds, which reveals a new unique way of discovering novel materials that can combine the advantages of inorganic compounds and small molecules. A direct application of our work is the proof of turning the common table salt and other similar ionic compounds into hydrogen storage materials under moderate pressures.


title: Changes: The Past, Present, and Future of the Nearby Dwarf Galaxy NGC 55

abstract: Spectra that cover wavelengths from 0.6 to 1.1um are used to examine the behavior of emission and absorption features in a contiguous 22 x 300 arcsec region centered on the nearby dwarf galaxy NGC 55. Based on the relative strengths of various emission features measured over spatial scales of many tens of parsecs, it is concluded that the ionization states and sulphur abundances in most of the star-forming regions near the center of NGC 55 are similar. A large star-forming region is identified in the north west part of the disk at a projected distance of ~1 kpc from the galaxy center that has distinct ionization properties. Fossil star-forming regions are also identified using the depth of the near-infrared Ca triplet. One such area is identified near the intersection of the major and minor axes, and it is suggested that this area is a proto-nucleus. The spectra of bright unresolved sources that are blended stellar asterisms, compact HII regions, and star clusters are also discussed. The spectra of some of the HII regions contain Ca triplet absorption lines, signalling a concentration of stars in the resolution element that span many Myr. Six of the unresolved sources have spectroscopic characteristics that are indicative of C stars embedded in intermediate age clusters. The peculiar properties of NGC 55 have been well documented in the literature, and it is argued that these may indicate that NGC 55 is transforming into a dwarf lenticular galaxy.


title: Double neutron star mergers from hierarchical triple-star systems

abstract: The isolated binary evolution model for merging neutron stars (NSs) involves processes such as mass transfer, common-envelope evolution, and natal kicks, all of which are poorly understood. Also, the predicted NS-NS merger rates are typically lower than the rates inferred from the LIGO GW170817 event. Here, we investigate merger rates of NS and black hole (BH)-NS binaries in hierarchical triple-star systems. In such systems, the tertiary can induce Lidov-Kozai (LK) oscillations in the inner binary, accelerating its coalescence, and potentially enhancing compact object merger rates. However, since compact objects originate from massive stars, the prior evolution should also be taken into account. Natal kicks, in particular, could significantly reduce the rates by unbinding the tertiary before it can affect the inner binary through LK evolution. We carry out simulations of massive triples taking into account stellar evolution starting from the main sequence, secular and tidal evolution, and the effects of supernovae. For large NS birth kicks (\sigma_k = 265 km/s), we find that the triple NS-NS merger rate (several hundred Gpc^{-3} yr^{-1}) is lower by a factor of ~2-3 than the binary rate, but for no kicks (\sigma_k = 0 km/s), the triple rate (several thousand Gpc^{-3} yr^{-1}) is comparable to the binary rate. Our results indicate that a significant fraction of NS-NS mergers could originate from triples if a substantial portion of the NS population is born with low kick velocities, as indicated by other work. However, uncertainties and open questions remain because of our simplifying assumption of dynamical decoupling after inner binary interaction has been triggered.


title: MaxiMask and MaxiTrack: two new tools for identifying contaminants in astronomical images using convolutional neural networks

abstract: In this work, we propose two convolutional neural network classifiers for detecting contaminants in astronomical images. Once trained, our classifiers are able to identify various contaminants, such as cosmic rays, hot and bad pixels, persistence effects, satellite or plane trails, residual fringe patterns, nebulous features, saturated pixels, diffraction spikes, and tracking errors in images. They encompass a broad range of ambient conditions, such as seeing, image sampling, detector type, optics, and stellar density. The first classifier, MaxiM ask , performs semantic segmentation and generates bad pixel maps for each contaminant, based on the probability that each pixel belongs to a given contaminant class. The second classifier, MaxiTrack , classifies entire images and mosaics, by computing the probability for the focal plane to be affected by tracking errors. We gathered training and testing data from real data originating from various modern charged-coupled devices and near-infrared cameras, that are augmented with image simulations. We quantified the performance of both classifiers and show that M axi M ask achieves state-of-the-art performance for the identification of cosmic ray hits. Thanks to a built-in Bayesian update mechanism, both classifiers can be tuned to meet specific science goals in various observational contexts.


title: Thermal evolution and quiescent emission of transiently accreting neutron stars

abstract: We study long-term thermal evolution of neutron stars in soft X-ray transients (SXTs), taking the deep crustal heating into account consistently with the changes of the composition of the crust. We collect observational estimates of average accretion rates and thermal luminosities of such neutron stars and compare the theory with observations. We perform simulations of thermal evolution of accreting neutron stars, considering the gradual replacement of the original nonaccreted crust by the reprocessed accreted matter, the neutrino and photon energy losses, and the deep crustal heating due to nuclear reactions in the accreted crust. We test and compare results for different modern theoretical models. We update a compilation of the observational estimates of the thermal luminosities in quiescence and average accretion rates in the SXTs and compare the observational estimates with the theoretical results. Long-term thermal evolution of transiently accreting neutron stars is nonmonotonic. The quasi-equilibrium temperature in quiescence reaches a minimum and then increases toward the final steady state. The quasi-equilibrium thermal luminosity of a neutron star in an SXT can be substantially lower at the minimum than in the final state. This enlarges the range of possibilities for theoretical interpretation of observations of such neutron stars. The updates of the theory and observations leave unchanged the previous conclusions that the direct Urca process operates in relatively cold neutron stars and that an accreted heat-blanketing envelope is likely present in relatively hot neutron stars in the SXTs in quiescence. The results of the comparison of theory with observations favor suppression of the triplet pairing type of nucleon superfluidity in the neutron-star matter.


title: Approximation by group invariant subspaces

abstract: In this article we study the structure of $\Gamma$-invariant spaces of $L^2(\bf R)$. Here $\bf R$ is a second countable LCA group. The invariance is with respect to the action of $\Gamma$, a non commutative group in the form of a semidirect product of a discrete cocompact subgroup of $\bf R$ and a group of automorphisms. This class includes in particular most of the crystallographic groups. We obtain a complete characterization of $\Gamma$-invariant subspaces in terms of range functions associated to shift-invariant spaces. We also define a new notion of range function adapted to the $\Gamma$-invariance and construct Parseval frames of orbits of some elements in the subspace, under the group action. These results are then applied to prove the existence and construction of a $\Gamma$-invariant subspace that best approximates a set of functional data in $L^2(\bf R)$. This is very relevant in applications since in the euclidean case, $\Gamma$-invariant subspaces are invariant under rigid movements, a very sought feature in models for signal processing.


title: Forward-Backward Drell-Yan Asymmetry and PDF Determination

abstract: We investigate the impact of high-statistics Drell-Yan (DY) measurements at the LHC on the study of non-perturbative QCD effects from parton distribution functions (PDF). We present the results of a PDF profiling analysis based on the neutral-current DY forward-backward asymmetry, using the open source fit platform xFitter.


title: Quantitative Impact Evaluation of an Abstraction Layer for Data Stream Processing Systems

abstract: With the demand to process ever-growing data volumes, a variety of new data stream processing frameworks have been developed. Moving an implementation from one such system to another, e.g., for performance reasons, requires adapting existing applications to new interfaces. Apache Beam addresses these high substitution costs by providing an abstraction layer that enables executing programs on any of the supported streaming frameworks. In this paper, we present a novel benchmark architecture for comparing the performance impact of using Apache Beam on three streaming frameworks: Apache Spark Streaming, Apache Flink, and Apache Apex. We find significant performance penalties when using Apache Beam for application development in the surveyed systems. Overall, usage of Apache Beam for the examined streaming applications caused a high variance of query execution times with a slowdown of up to a factor of 58 compared to queries developed without the abstraction layer. All developed benchmark artifacts are publicly available to ensure reproducible results.


title: Fully-automated deep learning-powered system for DCE-MRI analysis of brain tumors

abstract: Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) plays an important role in diagnosis and grading of brain tumor. Although manual DCE biomarker extraction algorithms boost the diagnostic yield of DCE-MRI by providing quantitative information on tumor prognosis and prediction, they are time-consuming and prone to human error. In this paper, we propose a fully-automated, end-to-end system for DCE-MRI analysis of brain tumors. Our deep learning-powered technique does not require any user interaction, it yields reproducible results, and it is rigorously validated against benchmark (BraTS'17 for tumor segmentation, and a test dataset released by the Quantitative Imaging Biomarkers Alliance for the contrast-concentration fitting) and clinical (44 low-grade glioma patients) data. Also, we introduce a cubic model of the vascular input function used for pharmacokinetic modeling which significantly decreases the fitting error when compared with the state of the art, alongside a real-time algorithm for determination of the vascular input region. An extensive experimental study, backed up with statistical tests, showed that our system delivers state-of-the-art results (in terms of segmentation accuracy and contrast-concentration fitting) while requiring less than 3 minutes to process an entire input DCE-MRI study using a single GPU.


title: A Constant Factor Approximation for Capacitated Min-Max Tree Cover

abstract: Given a graph $G=(V,E)$ with non-negative real edge lengths and an integer parameter $k$, the Min-Max k-Tree Cover problem seeks to find a set of at most $k$ subtrees of $G$, such that the union of the trees is the vertex set $V$. The objective is to minimize the maximum length among all the trees. We give the first constant factor approximation for the hard uniform capacitated version of this problem, where, an input parameter $\lambda$ upper bounds the number of vertices that can be covered by any of the trees. Our result extends to the rooted version of the problem, where we are given a set of $k$ root vertices, $R$ and each of the covering trees is required to include a distinct vertex in $R$ as the root. Prior to our work, the only result known was a $(2k-1)$-approximation algorithm for the special case when the total number of vertices in the graph is $k\lambda$ [Guttmann-Beck and Hassin, J. of Algorithms, 1997]. Our technique circumvents the difficulty of using the minimum spanning tree of the graph as a lower bound, which is standard for the uncapacitated version of the problem [Even et al., OR Letters 2004] [Khani et al., Algorithmica 2010]. Instead, we use Steiner trees that cover $\lambda$ vertices along with an iterative refinement procedure that ensures that the output trees have low cost and the vertices are well distributed among the trees.


title: A variational finite volume scheme for Wasserstein gradient flows

abstract: We propose a variational finite volume scheme to approximate the solutions to Wasserstein gradient flows. The time discretization is based on an implicit linearization of the Wasserstein distance expressed thanks to Benamou-Brenier formula, whereas space discretization relies on upstream mobility two-point flux approximation finite volumes. Our scheme is based on a first discretize then optimize approach in order to preserve the variational structure of the continuous model at the discrete level. Our scheme can be applied to a wide range of energies, guarantees non-negativity of the discrete solutions as well as decay of the energy. We show that our scheme admits a unique solution whatever the convex energy involved in the continuous problem, and we prove its convergence in the case of the linear Fokker-Planck equation with positive initial density. Numerical illustrations show that it is first order accurate in both time and space, and robust with respect to both the energy and the initial profile.


title: A Polynomial Time Algorithm for Log-Concave Maximum Likelihood via Locally Exponential Families

abstract: We consider the problem of computing the maximum likelihood multivariate log-concave distribution for a set of points. Specifically, we present an algorithm which, given $n$ points in $\mathbb{R}^d$ and an accuracy parameter $\epsilon>0$, runs in time $poly(n,d,1/\epsilon),$ and returns a log-concave distribution which, with high probability, has the property that the likelihood of the $n$ points under the returned distribution is at most an additive $\epsilon$ less than the maximum likelihood that could be achieved via any log-concave distribution. This is the first computationally efficient (polynomial time) algorithm for this fundamental and practically important task. Our algorithm rests on a novel connection with exponential families: the maximum likelihood log-concave distribution belongs to a class of structured distributions which, while not an exponential family, "locally" possesses key properties of exponential families. This connection then allows the problem of computing the log-concave maximum likelihood distribution to be formulated as a convex optimization problem, and solved via an approximate first-order method. Efficiently approximating the (sub) gradients of the objective function of this optimization problem is quite delicate, and is the main technical challenge in this work.


title: XferNAS: Transfer Neural Architecture Search

abstract: The term Neural Architecture Search (NAS) refers to the automatic optimization of network architectures for a new, previously unknown task. Since testing an architecture is computationally very expensive, many optimizers need days or even weeks to find suitable architectures. However, this search time can be significantly reduced if knowledge from previous searches on different tasks is reused. In this work, we propose a generally applicable framework that introduces only minor changes to existing optimizers to leverage this feature. As an example, we select an existing optimizer and demonstrate the complexity of the integration of the framework as well as its impact. In experiments on CIFAR-10 and CIFAR-100, we observe a reduction in the search time from 200 to only 6 GPU days, a speed up by a factor of 33. In addition, we observe new records of 1.99 and 14.06 for NAS optimizers on the CIFAR benchmarks, respectively. In a separate study, we analyze the impact of the amount of source and target data. Empirically, we demonstrate that the proposed framework generally gives better results and, in the worst case, is just as good as the unmodified optimizer.


title: Magnetic effects in general relativity: quadrupoles, shear and memory

abstract: Gravitational effect of "magnetic type" -- those having a curl-like character over large spheres -- are investigated, for isolated systems. The Bondi-Sachs-Newman-Penrose formalism clarifies a number of points, especially related to radiation memory. It is shown that the "memory tensor" is equivalent to the change in Bondi shear, from before to after the emission of radiation. This means that if magnetic radiation memory is present, at least one of the intervals bracketing the radiation must have non-zero magnetic shear but vanishing radiation. Such intervals, called here CPMS regimes, are shown to be necessarily non-stationary, however, raising a variety of technical and interpretative issues. In linearized general relativity, the gravitational fields due to point magnetic quadrupoles with arbitrary time-dependence are computed, and some of their physics studied. In the far zone, there is a red-shift effect which could be searched for astrophysically: light coming from behind a source generating magnetic shear would be red-shifted by an amount varying with the angle around the source of shear, and in the far-field limit this red-shift goes inversely with the impact parameter. Induction-zone effects are also considered. An induction-zone memory effect should exist which could possibly be within the reach of laboratory experiments, but no good candidates for astrophysically detectable effects are found. Also a quadrupole will induce test particles to move in such a fashion as to create an opposing quadrupole, an effect reminiscent of Lenz's law.


title: A roadmap for Generalized Plane Waves and their interpolation properties

abstract: This work focuses on the study of partial differential equation (PDE) based basis function for Discontinuous Galerkin methods to solve numerically wave-related boundary value problems with variable coefficients. To tackle problems with constant coefficients, wave-based methods have been widely studied in the literature: they rely on the concept of Trefftz functions, i.e. local solutions to the governing PDE, using oscillating basis functions rather than polynomial functions to represent the numerical solution. Generalized Plane Waves (GPWs) are an alternative developed to tackle problems with variable coefficients, in which case Trefftz functions are not available. In a similar way, they incorporate information on the PDE, however they are only approximate Trefftz functions since they don't solve the governing PDE exactly, but only an approximated PDE. Considering a new set of PDEs beyond the Helmholtz equation, we propose to set a roadmap for the construction and study of local interpolation properties of GPWs. Identifying carefully the various steps of the process, we provide an algorithm to summarize the construction of these functions, and establish necessary conditions to obtain high order interpolation properties of the corresponding basis.


title: Deep Perceptual Compression

abstract: Several deep learned lossy compression techniques have been proposed in the recent literature. Most of these are optimized by using either MS-SSIM (multi-scale structural similarity) or MSE (mean squared error) as a loss function. Unfortunately, neither of these correlate well with human perception and this is clearly visible from the resulting compressed images. In several cases, the MS-SSIM for deep learned techniques is higher than say a conventional, non-deep learned codec such as JPEG-2000 or BPG. However, the images produced by these deep learned techniques are in many cases clearly worse to human eyes than those produced by JPEG-2000 or BPG. We propose the use of an alternative, deep perceptual metric, which has been shown to align better with human perceptual similarity. We then propose Deep Perceptual Compression (DPC) which makes use of an encoder-decoder based image compression model to jointly optimize on the deep perceptual metric and MS-SSIM. Via extensive human evaluations, we show that the proposed method generates visually better results than previous learning based compression methods and JPEG-2000, and is comparable to BPG. Furthermore, we demonstrate that for tasks like object-detection, images compressed with DPC give better accuracy.


title: White Paper on New Opportunities at the Next-Generation Neutrino Experiments (Part 1: BSM Neutrino Physics and Dark Matter)

abstract: With the advent of a new generation of neutrino experiments which leverage high-intensity neutrino beams for precision measurements, it is timely to explore physics topics beyond the standard neutrino-related physics. Given that the realm of beyond the standard model (BSM) physics has been mostly sought at high-energy regimes at colliders, such as the LHC at CERN, the exploration of BSM physics in neutrino experiments will enable complementary measurements at the energy regimes that balance that of the LHC. This is in concert with new ideas for high-intensity beams for fixed target and beam-dump experiments world-wide, e.g., those at CERN. The combination of the high intensity proton beam facilities and massive detectors for precision neutrino oscillation parameter measurements and for CP violation phase measurements will help make BSM physics reachable even in low energy regimes in accelerator based experiments. Large mass detectors with highly precise tracking and energy measurements, excellent timing resolution, and low energy thresholds will enable searches for BSM phenomena from cosmogenic origin, as well. Therefore, it is conceivable that BSM topics in the next generation neutrino experiments could be the dominant physics topics in the foreseeable future, as the precision of the neutrino oscillation parameter and CPV measurements continues to improve. In this spirit, this white paper provides a review of the current landscape of BSM theory in neutrino experiments in two selected areas of the BSM topics - dark matter and neutrino related BSM - and summarizes the current results from existing neutrino experiments to set benchmarks for both theory and experiment. This paper then provides a review of upcoming neutrino experiments throughout the next 10 - 15 year time scale and their capabilities to set the foundation for potential reach in BSM physics in the two aforementioned themes.


title: Production of electromagnetic and sound radiation from gravitational waves in core-collapse supernovae

abstract: This paper presents the scenario that gravitational waves, generated in core-collapse of a pre-supernova star, can produce both electromagnetic radiation and sound radiation as the gravitational wave propagates outward from the collapsing core. While the energy of this co-produced electromagnetic and/or sound radiation is orders of magnitude smaller than the initiating gravitational radiation, the electromagnetic/sonic power may be sufficient to re-ignite fusion outside the collapsing core. The non-equilibrium re-ignition of fusion, in roughly the same time frame as the strongest neutrino emissions, would change the configuration of the pre-supernova star. Although the co-produced electromagnetic/sonic radiation is not powerful enough to contribute directly to the supernova, the associated non-equilibrium re-ignition of fusion could alter the state outside the core. The aim of this paper is to argue that the effect of this hypothesized co-produced electromagnetic/sonic radiation should be included in computational models of core collapse supernovae.


title: Learning High-Level Planning Symbols from Intrinsically Motivated Experience

abstract: In symbolic planning systems, the knowledge on the domain is commonly provided by an expert. Recently, an automatic abstraction procedure has been proposed in the literature to create a Planning Domain Definition Language (PDDL) representation, which is the most widely used input format for most off-the-shelf automated planners, starting from `options', a data structure used to represent actions within the hierarchical reinforcement learning framework. We propose an architecture that potentially removes the need for human intervention. In particular, the architecture first acquires options in a fully autonomous fashion on the basis of open-ended learning, then builds a PDDL domain based on symbols and operators that can be used to accomplish user-defined goals through a standard PDDL planner. We start from an implementation of the above mentioned procedure tested on a set of benchmark domains in which a humanoid robot can change the state of some objects through direct interaction with the environment. We then investigate some critical aspects of the information abstraction process that have been observed, and propose an extension that mitigates such criticalities, in particular by analysing the type of classifiers that allow a suitable grounding of symbols.


title: Is there a non-stationary {\gamma}-ray emission zone 42 pc away from the 3C 279 core?

abstract: We investigate the relationship between the variable gamma-ray emission and jet properties in the blazar 3C 279, by combining the Fermi/LAT data spanning a period of eight years with concurrent radio measurements made at multiple epochs with VLBA at 15 and 43 GHz within the MOJAVE and VLBA-BU monitoring programs. The aim of this paper is to compare the flux variability of the different components found in the VLBA observations, to the variability in the gamma-rays. This analysis helps to investigate whether any of the jet components can be associated with the gamma-ray variability. Through Spearman rank correlation we found that the gamma-ray variability is correlated with a particular region (feature B, in the MOJAVE images) downstream from the observed base (core) of the jet. This jet component is, therefore, a likely location at which an important fraction of the variable gamma-ray emission is produced. We also calculated the average proper motion of the component with respect to the VLBA core and found that it moves at an apparent superluminal velocity of $(3.70 \pm 0.35)c$, implying that one of the gamma-ray emission zones is not stationary. This jet component is also found between 6.86 mas and 8.68 mas, which translates to a distance from the radio core of at least 42 pc.


title: Turbulence statistics in a negatively buoyant multiphase plume

abstract: We investigate the turbulence statistics in a {multiphase plume made of heavy particles (particle Reynolds number at terminal velocity is 450)}. Using refractive-index-matched stereoscopic particle image velocimetry, we measure the locations of particles {whose buoyancy drives the formation of a multiphase plume,} {together with the local velocity of the induced flow in the ambient salt-water}. {Measurements in the plume centerplane exhibit self-similarity in mean flow characteristics consistent with classic integral plume theories.} The turbulence characteristics resemble those measured in a bubble plume, {including strong anisotropy in the normal Reynolds stresses. However, we observe structural differences between the two multiphase plumes. First, the skewness of the probability density function (PDF) of the axial velocity fluctuations is not that which would be predicted by simply reversing the direction of a bubble plume. Second, in contrast to a bubble plume, the particle plume has a non-negligible fluid-shear production term in the turbulent kinetic energy (TKE) budget. Third, the radial decay of all measured terms in the TKE budget is slower than those in a bubble plume.} Despite these dissimilarities, a bigger picture emerges that applies to both flows. The TKE production by particles (or bubbles) roughly balances the viscous dissipation, except near the plume centerline. The one-dimensional power-spectra of the velocity fluctuations show a -3 power-law that puts both the particle and bubble plume in a category different from single-phase shear-flow turbulence.


title: Semi-Lagrangian Vlasov simulation on GPUs

abstract: In this paper, our goal is to efficiently solve the Vlasov equation on GPUs. A semi-Lagrangian discontinuous Galerkin scheme is used for the discretization. Such kinetic computations are extremely expensive due to the high-dimensional phase space. The SLDG code, which is publicly available under the MIT license abstracts the number of dimensions and uses a shared codebase for both GPU and CPU based simulations. We investigate the performance of the implementation on a range of both Tesla (V100, Titan V, K80) and consumer (GTX 1080 Ti) GPUs. Our implementation is typically able to achieve a performance of approximately 470 GB/s on a single GPU and 1600 GB/s on four V100 GPUs connected via NVLink. This results in a speedup of about a factor of ten (comparing a single GPU with a dual socket Intel Xeon Gold node) and approximately a factor of 35 (comparing a single node with and without GPUs). In addition, we investigate the effect of single precision computation on the performance of the SLDG code and demonstrate that a template based dimension independent implementation can achieve good performance regardless of the dimensionality of the problem.


title: Characterization and decomposition of the natural van der Waals heterostructure SnSb2Te4 under compression

abstract: This joint experimental and theoretical study of the structural, vibrational and electrical properties of rhombohedral SnSb2Te4 at high pressure unveils the internal mechanisms of its compression. The equation of state and the internal polyhedral compressibility, the symmetry and behavior of the Raman-active modes and the electrical behavior of this topological insulator under compression have been discussed and compared with the parent binary alpha-Sb2Te3 and SnTe compounds and with related ternary compounds. Our X-ray diffraction and Raman measurements together with theoretical calculations, which include topological electron density and electronic localization function analysis, evidence the presence of an isostructural phase transition around 2 GPa and a Fermi resonance around 3.5 GPa. The Raman spectrum of SnSb2Te4 shows vibrational modes that are forbidden in rocksalt SnTe; thus showing a novel way to experimentally observe the forbidden vibrational modes of some compounds. Additionally, since SnSb2Te4 is an incipient metal, like its parent binary compounds, we establish a new criterion to identify the recently proposed metavalent bonding in complex materials when different bond characters coexist in the system. Finally, SnSb2Te4 exhibits a pressure-induced decomposition into the high-pressure phases of its parent binary compounds above 7 GPa, which is supported by an analysis of their formation enthalpies. We have framed the behavior of SnSb2Te4 within the extended orbital radii map of BA2Te4 compounds, which paves the way to understand the pressure behavior and stability ranges of other layered van der Waals-type compounds with similar stoichiometry.


title: A study of convex convex-composite functions via infimal convolution with applications

abstract: In this note we provide a full conjugacy and subdifferential calculus for convex convex-composite functions in finite-dimensional space. Our approach, based on infimal convolution and cone-convexity, is straightforward and yields the desired results under a verifiable Slater-type condition, with relaxed monotonicity and without lower semicontinuity assumptions on the functions in play. The versatility of our findings is illustrated by a series of applications in optimization and matrix analysis, including conic programming, matrix-fractional, variational Gram, and spectral functions.


title: ALMA Observations of the Terahertz Spectrum of Sagittarius A*

abstract: We present ALMA observations at 233, 678, and 870 GHz of the Galactic Center black hole, Sagittarius A*. These observations reveal a flat spectrum over this frequency range with spectral index $\alpha \approx -0.3$, where the flux density $S \propto \nu^\alpha$. We model the submm and far infrared spectrum with a one zone synchrotron model of thermal electrons. We infer electron densities $n = (2-5) \times 10^6$ cm$^{-3}$, electron temperatures $T_e = (1-3) \times 10^{11}$ K, and magnetic field strength $B = 10-50$ G. The parameter range can be further constrained using the observed quiescent X-ray luminosity. The flat submm spectrum results in a high electron temperature and implies that the emitting electrons are efficiently heated. We also find that the emission is most likely optically thin at 233 GHz. These results indicate that millimeter and submillimeter wavelength very long baseline interferometry of Sgr A* including those of the Event Horizon Telescope should see a transparent emission region down to event horizon scales.


title: Multi-Task Regression-based Learning for Autonomous Unmanned Aerial Vehicle Flight Control within Unstructured Outdoor Environments

abstract: Increased growth in the global Unmanned Aerial Vehicles (UAV) (drone) industry has expanded possibilities for fully autonomous UAV applications. A particular application which has in part motivated this research is the use of UAV in wide area search and surveillance operations in unstructured outdoor environments. The critical issue with such environments is the lack of structured features that could aid in autonomous flight, such as road lines or paths. In this paper, we propose an End-to-End Multi-Task Regression-based Learning approach capable of defining flight commands for navigation and exploration under the forest canopy, regardless of the presence of trails or additional sensors (i.e. GPS). Training and testing are performed using a software in the loop pipeline which allows for a detailed evaluation against state-of-the-art pose estimation techniques. Our extensive experiments demonstrate that our approach excels in performing dense exploration within the required search perimeter, is capable of covering wider search regions, generalises to previously unseen and unexplored environments and outperforms contemporary state-of-the-art techniques.


title: SentiMATE: Learning to play Chess through Natural Language Processing

abstract: We present SentiMATE, a novel end-to-end Deep Learning model for Chess, employing Natural Language Processing that aims to learn an effective evaluation function assessing move quality. This function is pre-trained on the sentiment of commentary associated with the training moves and is used to guide and optimize the agent's game-playing decision making. The contributions of this research are three-fold: we build and put forward both a classifier which extracts commentary describing the quality of Chess moves in vast commentary datasets, and a Sentiment Analysis model trained on Chess commentary to accurately predict the quality of said moves, to then use those predictions to evaluate the optimal next move of a Chess agent. Both classifiers achieve over 90 % classification accuracy. Lastly, we present a Chess engine, SentiMATE, which evaluates Chess moves based on a pre-trained sentiment evaluation function. Our results exhibit strong evidence to support our initial hypothesis - "Can Natural Language Processing be used to train a novel and sample efficient evaluation function in Chess Engines?" - as we integrate our evaluation function into modern Chess engines and play against agents with traditional Chess move evaluation functions, beating both random agents and a DeepChess implementation at a level-one search depth - representing the number of moves a traditional Chess agent (employing the alpha-beta search algorithm) looks ahead in order to evaluate a given chess state.


title: MIMIC-Extract: A Data Extraction, Preprocessing, and Representation Pipeline for MIMIC-III

abstract: Robust machine learning relies on access to data that can be used with standardized frameworks in important tasks and the ability to develop models whose performance can be reasonably reproduced. In machine learning for healthcare, the community faces reproducibility challenges due to a lack of publicly accessible data and a lack of standardized data processing frameworks. We present MIMIC-Extract, an open-source pipeline for transforming raw electronic health record (EHR) data for critical care patients contained in the publicly-available MIMIC-III database into dataframes that are directly usable in common machine learning pipelines. MIMIC-Extract addresses three primary challenges in making complex health records data accessible to the broader machine learning community. First, it provides standardized data processing functions, including unit conversion, outlier detection, and aggregating semantically equivalent features, thus accounting for duplication and reducing missingness. Second, it preserves the time series nature of clinical data and can be easily integrated into clinically actionable prediction tasks in machine learning for health. Finally, it is highly extensible so that other researchers with related questions can easily use the same pipeline. We demonstrate the utility of this pipeline by showcasing several benchmark tasks and baseline results.


title: Universal sets for ideals

abstract: In this paper we consider a notion of universal sets for ideals. We show that there exist universal sets of minimal Borel complexity for classic ideals like null subsets of $2^\omega$ and meager subsets of any Polish space, and demonstrate that the existence of such sets is helpful in establishing some facts about the real line in generic extensions. We also construct universal sets for $\mathcal{E}$ - the $\sigma$-ideal generated by closed null subsets of $2^\omega$, and for some ideals connected with forcing notions: $\mathcal{K}_\sigma$ subsets of $\omega^{\omega}$ and the Laver ideal. We also consider Fubini products of ideals and show that there are $\Sigma^0_3$ universal sets for $\mathcal{N}\otimes\mathcal{M}$ and $\mathcal{M}\otimes\mathcal{N}$.


title: Dynamics of relaxation and dressing of a quenched Bose polaron

abstract: We study the non-equilibrium dynamics of relaxation and dressing of a mobile impurity suddenly immersed--or quenched-- into a zero temperature homogeneous Bose Einstein condensate (BEC) with velocity $v$. A many body generalization of Weisskopf-Wigner theory is implemented to obtain the impurity fidelity, reduced density matrix and entanglement entropy. The dynamics depend crucially on the Mach number $\beta =v/c$, with $c$ the speed of sound of superfluid phonons and features many different time scales. Quantum Zeno behavior at early time is followed by relaxational and dressing dynamics determined by Cerenkov emission of long-wavelength phonons for $\beta >1$ with a decay rate $\Gamma_p \propto (\beta-1)^3$. The polaron dressing dynamics \emph{slows-down} as $\beta \rightarrow 1$ and is characterized by power laws $t^{-\alpha}$ with different exponents for $\beta \lessgtr 1$. The asymptotic entanglement entropy features a sharp discontinuity and the residue features a cusp at $\beta =1$. These non-equilibrium features suggest \emph{universal} dynamical critical phenomena near $\beta \simeq 1$, and are a direct consequence of the linear dispersion relation of long wavelength superfluid phonons. We conjecture on the emergence of an asymptotic dynamical attractor with $\beta \leq 1$.


title: Scalable Topological Data Analysis and Visualization for Evaluating Data-Driven Models in Scientific Applications

abstract: With the rapid adoption of machine learning techniques for large-scale applications in science and engineering comes the convergence of two grand challenges in visualization. First, the utilization of black box models (e.g., deep neural networks) calls for advanced techniques in exploring and interpreting model behaviors. Second, the rapid growth in computing has produced enormous datasets that require techniques that can handle millions or more samples. Although some solutions to these interpretability challenges have been proposed, they typically do not scale beyond thousands of samples, nor do they provide the high-level intuition scientists are looking for. Here, we present the first scalable solution to explore and analyze high-dimensional functions often encountered in the scientific data analysis pipeline. By combining a new streaming neighborhood graph construction, the corresponding topology computation, and a novel data aggregation scheme, namely topology aware datacubes, we enable interactive exploration of both the topological and the geometric aspect of high-dimensional data. Following two use cases from high-energy-density (HED) physics and computational biology, we demonstrate how these capabilities have led to crucial new insights in both applications.


title: I Stand With You: Using Emojis to Study Solidarity in Crisis Events

abstract: We study how emojis are used to express solidarity in social media in the context of two major crisis events - a natural disaster, Hurricane Irma in 2017 and terrorist attacks that occurred on November 2015 in Paris. Using annotated corpora, we first train a recurrent neural network model to classify expressions of solidarity in text. Next, we use these expressions of solidarity to characterize human behavior in online social networks, through the temporal and geospatial diffusion of emojis. Our analysis reveals that emojis are a powerful indicator of sociolinguistic behaviors (solidarity) that are exhibited on social media as the crisis events unfold.


title: Explicit bounds on exceptional zeroes of Dirichlet L-function II

abstract: This paper improves the upper bound for the exceptional zeroes of Dirichlet L-functions with even characters. The result is obtained by improving on explicit estimate for $L'(\sigma;\chi)$ for $\sigma$ close to unity, using a result on the average of Dirichlet characters, and on the lower bound for $L(1;\chi)$, with computational aid.


title: A multiscale Laplacian of Gaussian (LoG) filtering approach to pulmonary nodule detection from whole-lung CT scans

abstract: Candidate generation, the first stage for most computer aided detection (CAD) systems, rapidly scans the entire image data for any possible abnormality locations, while the subsequent stages of the CAD system refine the candidates list to determine the most probable or significant of these candidates. The candidate generator creates a list of the locations and provides a size estimate for each candidate. A multiscale scale-normalized Laplacian of Gaussian (LoG) filtering method for detecting pulmonary nodules in whole-lung CT scans, presented in this paper, achieves a high sensitivity for both solid and nonsolid pulmonary nodules. The pulmonary nodule LoG filtering method was validated on a size-enriched database of 706 whole-lung low-dose CT scans containing 499 solid (>= 4 mm) and 107 nonsolid (>= 6 mm) pulmonary nodules. The method achieved a sensitivity of 0.998 (498/499) for solid nodules and a sensitivity of 1.000 (107/107) for nonsolid nodules. Furthermore, compared to radiologist measurements, the method provided low average nodule size estimation error of 0.12 mm for solid and 1.27 mm for nonsolid nodules. The average distance between automatically and manually determined nodule centroids were 1.41 mm and 1.43 mm, respectively.


title: Control of LED Emission with Functional Dielectric Metasurfaces

abstract: The improvement of light-emitting diodes (LEDs) is one of the major goals of optoelectronics and photonics research. While emission rate enhancement is certainly one of the targets, in this regard, for LED integration to complex photonic devices, one would require to have, additionally, precise control of the wavefront of the emitted light. Metasurfaces are spatial arrangements of engineered scatters that may enable this light manipulation capability with unprecedented resolution. Most of these devices, however, are only able to function properly under irradiation of light with a large spatial coherence, typically normally incident lasers. LEDs, on the other hand, have angularly broad, Lambertian-like emission patterns characterized by a low spatial coherence, which makes the integration of metasurface devices on LED architectures extremely challenging. A novel concept for metasurface integration on LED is proposed, using a cavity to increase the LED spatial coherence through an angular collimation. Due to the resonant character of the cavity, extending the spatial coherence of the emitted light does not come at the price of any reduction in the total emitted power. The experimental demonstration of the proposed concept is implemented on a GaP LED architecture including a hybrid metallic-Bragg cavity. By integrating a silicon metasurface on top we demonstrate two different functionalities of these compact devices: directional LED emission at a desired angle and LED emission of a vortex beam with an orbital angular momentum. The presented concept is general, being applicable to other incoherent light sources and enabling metasurfaces designed for plane waves to work with incoherent light emitters.


title: Measurement of differential cross sections and charge ratios for $t$-channel single top quark production in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A measurement is presented of differential cross sections for $t$-channel single top quark and antiquark production in proton-proton collisions at a centre-of-mass energy of 13 TeV by the CMS experiment at the LHC. From a data set corresponding to an integrated luminosity of 35.9 fb$^{-1}$, events containing one muon or electron and two or three jets are analysed. The cross section is measured as a function of the top quark transverse momentum ($p_\mathrm{T}$), rapidity, and polarisation angle, the charged lepton $p_\mathrm{T}$ and rapidity, and the $p_\mathrm{T}$ of the W boson from the top quark decay. In addition, the charge ratio is measured differentially as a function of the top quark, charged lepton, and W boson kinematic observables. The results are found to be in agreement with standard model predictions using various next-to-leading-order event generators and sets of parton distribution functions. Additionally, the spin asymmetry, sensitive to the top quark polarisation, is determined from the differential distribution of the polarisation angle at parton level to be 0.440 $\pm$ 0.070, in agreement with the standard model prediction.


title: Recasting the Proof of Parseval's Identity

abstract: We generalize aspects of Fourier Analysis from intervals on $\mathbb{R}$ to bounded and measurable subsets of $\mathbb{R}^n$. In doing so, we obtain a few interesting results. The first is a new proof of the famous Integral Cauchy-Schwarz Inequality. The second is a restatement of Parseval's Identity that doubles as a representation of integrating bounded and measurable functions over bounded and measurable subsets of $\mathbb{R}^n$. Finally, we apply these first two results to develop some sufficient criteria for additional integral inequalities that are elementary in nature.


title: Hour-scale persistent negative anomaly of atmospheric electrostatic field near the epicenter before earthquake

abstract: Although earthquake prediction is a big challenge in the world, some simple observational tools can capture many physical signals and demonstrate that an earthquake (EQ) may be forthcoming in short period. Many researchers have studied the significant variation of atmospheric electrostatic field related to the forthcoming earthquake. However, until now, there is not a compelling physical mechanism which can explain why atmospheric electrostatic abnormal signal could appear just before an earthquake. Here we present a precursor signal and propose a brief physical interpretation. Under fair air conditions, if the near-surface atmospheric electrostatic field Ez (oriented down when it is positive) presents a very stable negative anomaly (from -100 V/m to -5000 V/m), it will forebode that an earthquake (seismicity from 3-8) would take place in the next several to tens of hours within a distance less than 100 km. We name this prediction technique as "DC Ez Determination"(DED). In addition, the mechanism of such abnormal quasi-static electric field before a forthcoming earthquake has been proposed here: (1) Radon gas releases from the rock clefts near the epicenter during seismogenic process. (2) {\alpha} particles are produced due to the radioactive decay of Radon gas in the air. (3) {\alpha} particle ionizing radiation creates more positive particles in the air, which is much more effective than that {\beta} and {\gamma} particles produced during Radon radioactive decay. (4) The new positive particles change formal positive atmospheric electric field (Ez) into stably negative near the earth-surface. (5) The closer the instrument is to the epicenter, the more likely it is to observe the negative Ez signal related to the earthquake. It is recommended to establish an instrument network to capture the reliable precursor and make some warnings before the earthquake disaster.


title: Toxicity Prediction by Multimodal Deep Learning

abstract: Prediction of toxicity levels of chemical compounds is an important issue in Quantitative Structure-Activity Relationship (QSAR) modeling. Although toxicity prediction has achieved significant progress in recent times through deep learning, prediction accuracy levels obtained by even very recent methods are not yet very high. We propose a multimodal deep learning method using multiple heterogeneous neural network types and data representations. We represent chemical compounds by strings, images, and numerical features. We train fully connected, convolutional, and recurrent neural networks and their ensembles. Each data representation or neural network type has its own strengths and weaknesses. Our motivation is to obtain a collective performance that could go beyond individual performance of each data representation or each neural network type. On a standard toxicity benchmark, our proposed method obtains significantly better accuracy levels than that by the state-of-the-art toxicity prediction methods.


title: When can we improve on sample average approximation for stochastic optimization?

abstract: We explore the performance of sample average approximation in comparison with several other methods for stochastic optimization when there is information available on the underlying true probability distribution. The methods we evaluate are (a) bagging; (b) kernel smoothing; (c) maximum likelihood estimation (MLE); and (d) a Bayesian approach. We use two test sets, the first has a quadratic objective function allowing for very different types of interaction between the random component and the univariate decision variable. Here the sample average approximation is remarkably effective and only consistently outperformed by a Bayesian approach. The second test set is a portfolio optimization problem in which we use different covariance structures for a set of 5 stocks. Here bagging, MLE and a Bayesian approach all do well.


title: Proceedings Seventeenth Conference on Theoretical Aspects of Rationality and Knowledge

abstract: This is the proceedings of the Seventeenth conference on Theoretical Aspects of Rationality and Knowledge, 17-19 July 2019, Institut de Recherche en Informatique de Toulouse (IRIT), Toulouse University Toulouse, France. The mission of the TARK conferences is to bring together researchers from a wide variety of fields, including Artificial Intelligence, Cryptography, Distributed Computing, Economics and Game Theory, Linguistics, Philosophy, and Psychology, in order to further our understanding of interdisciplinary issues involving reasoning about rationality and knowledge.


title: Override and update

abstract: Override and update are natural constructions for combining partial functions, which arise in various program specification contexts. We use an unexpected connection with combinatorial geometry to provide a complete finite system of equational axioms for the first order theory of the override and update constructions on partial functions, resolving the main unsolved problem in the area.


title: Online Set-Based Dynamic Analysis for Sound Predictive Race Detection

abstract: Predictive data race detectors find data races that exist in executions other than the observed execution. Smaragdakis et al. introduced the causally-precedes (CP) relation and a polynomial-time analysis for sound (no false races) predictive data race detection. However, their analysis cannot scale beyond analyzing bounded windows of execution traces. This work introduces a novel dynamic analysis called Raptor that computes CP soundly and completely. Raptor is inherently an online analysis that analyzes and finds all CP-races of an execution trace in its entirety. An evaluation of a prototype implementation of Raptor shows that it scales to program executions that the prior CP analysis cannot handle, finding data races that the prior CP analysis cannot find.


title: Batch Uniformization for Minimizing Maximum Anomaly Score of DNN-based Anomaly Detection in Sounds

abstract: Use of an autoencoder (AE) as a normal model is a state-of-the-art technique for unsupervised-anomaly detection in sounds (ADS). The AE is trained to minimize the sample mean of the anomaly score of normal sounds in a mini-batch. One problem with this approach is that the anomaly score of rare-normal sounds becomes higher than that of frequent-normal sounds, because the sample mean is strongly affected by frequent-normal samples, resulting in preferentially decreasing the anomaly score of frequent-normal samples. To decrease anomaly scores for both frequent- and rare-normal sounds, we propose batch uniformization, a training method for unsupervised-ADS for minimizing a weighted average of the anomaly score on each sample in a mini-batch. We used the reciprocal of the probabilistic density of each sample as the weight, more intuitively, a large weight is given for rare-normal sounds. Such a weight works to give a constant anomaly score for both frequent- and rare-normal sounds. Since the probabilistic density is unknown, we estimate it by using the kernel density estimation on each training mini-batch. Verification- and objective-experiments show that the proposed batch uniformization improves the performance of unsupervised-ADS.


title: Estimating the Moon to Earth radius ratio with a smartphone, a telescope and an eclipse

abstract: On January 20th, 2019, a total lunar eclipse was possible to be observed in Santiago, Chile. Using a smartphone attached to a telescope, photographs of the phenomenon were taken. With Earth's shadow on those images, and using textbook geometry, a simple open-source software and analytical procedures, we were allowed to calculate the ratio between the radii of the Moon and the Earth. The results are in very good agreement with the correct value for such ratio. This shows the strength of the smartphone technology to get powerful astronomical results in a very simple way and in a very short amount of time.


title: Only Time Can Tell: Discovering Temporal Data for Temporal Modeling

abstract: Understanding temporal information and how the visual world changes over time is a fundamental ability of intelligent systems. In video understanding, temporal information is at the core of many current challenges, including compression, efficient inference, motion estimation or summarization. However, in current video datasets it has been observed that action classes can often be recognized without any temporal information from a single frame of video. As a result, both benchmarking and training in these datasets may give an unintentional advantage to models with strong image understanding capabilities, as opposed to those with strong temporal understanding. In this paper we address this problem head on by identifying action classes where temporal information is actually necessary to recognize them and call these "temporal classes". Selecting temporal classes using a computational method would bias the process. Instead, we propose a methodology based on a simple and effective human annotation experiment. We remove just the temporal information by shuffling frames in time and measure if the action can still be recognized. Classes that cannot be recognized when frames are not in order are included in the temporal Dataset. We observe that this set is statistically different from other static classes, and that performance in it correlates with a network's ability to capture temporal information. Thus we use it as a benchmark on current popular networks, which reveals a series of interesting facts. We also explore the effect of training on the temporal dataset, and observe that this leads to better generalization in unseen classes, demonstrating the need for more temporal data. We hope that the proposed dataset of temporal categories will help guide future research in temporal modeling for better video understanding.


title: Efficient algorithm based on Liechtenstein method for computing exchange coupling constants using localized basis set

abstract: For large-scale computation of the exchange coupling constants $J_{ij}$, we reconstruct the Liechtenstein formula for localized orbital representation and simplify the energy integrations by adopting the finite pole approximation of the Fermi function proposed by Ozaki [Phys. Rev. B 75, 035123 (2007)]. We calculate the exchange coupling constant $J_{\mathrm{1NN}}$ of the first-nearest-neighbor sites in body-centered-cubic Fe systems of various sizes to estimate the optimal computational parameters that yield appropriate values at the lowest computational cost. It is shown that the number of poles needed for a computational accuracy of 0.05 meV is determined as $\sim$ 60, whereas the number of necessary Matsubara poles needed to obtain similar accuracy, which was determined in previous studies, is on the order of 1000. Finally, we show $J_{ij}$ as a function of atomic distance, and compared it with one derived from Korringa-Kohn-Rostoker Green's function formalism. The distance profile of $J_{ij}$ derived by KKR formalism agrees well with that derived by our study, and this agreement supports the reliability of our newly derived formalism.


title: The Luminosity Function and Formation Rate of A Complete Sample of Long Gamma-Ray Bursts

abstract: We study the luminosity function and formation rate of long gamma-ray bursts (GRBs) by using a maximum likelihood method. This is the first time this method is applied to a well-defined sample of GRBs that is complete in redshift. The sample is composed of 99 bursts detected by the $Swift$ satellite, 81 of them with measured redshift and luminosity for a completeness level of $82\%$. We confirm that a strong redshift evolution in luminosity (with an evolution index of $\delta=2.22^{+0.32}_{-0.31}$) or in density ($\delta=1.92^{+0.20}_{-0.21}$) is needed in order to reproduce the observations well. But since the predicted redshift and luminosity distributions in the two scenarios are very similar, it is difficult to distinguish between these two kinds of evolutions only on the basis of the current sample. Furthermore, we also consider an empirical density case in which the GRB rate density is directly described as a broken power-law function and the luminosity function is taken to be non-evolving. In this case, we find that the GRB formation rate rises like $(1+z)^{3.85^{+0.48}_{-0.45}}$ for $z\leq2$ and is proportional to $(1+z)^{-1.07^{+0.98}_{-1.12}}$ for $z\geq2$. The local GRB rate is $1.49^{+0.63}_{-0.64}$ Gpc$^{-3}$ yr$^{-1}$. The GRB rate may be consistent with the cosmic star formation rate (SFR) at $z\leq2$, but shows an enhancement compared to the SFR at $z\geq2$.


title: Hamiltonian of new general relativity using differential forms

abstract: In a recent work we derived the kinematic Hamiltonian and primary constraints of the new general relativity class of teleparallel gravity theories and showed that these theories can be grouped in 9 classes, based on the presence or absence of primary constraints in their Hamiltonian. Here we demonstrate an alternative approach towards this result, by using differential forms instead of tensor components throughout the calculation. We prove that also this alternative derivation yields the same results and show how they are related to each other.


title: Asymptotic Lech's inequality

abstract: We explore the classical Lech's inequality relating the Hilbert--Samuel multiplicity and colength of an $\mathfrak{m}$-primary ideal in a Noetherian local ring $(R,\mathfrak{m})$. We prove optimal versions of Lech's inequality for sufficiently deep ideals in characteristic $p>0$, and we conjecture that they hold in all characteristics. Our main technical result shows that if $(R,\mathfrak{m})$ has characteristic $p>0$ and $\widehat{R}$ is reduced, equidimensional, and has an isolated singularity, then for any sufficiently deep $\mathfrak{m}$-primary ideal $I$, the colength and Hilbert--Kunz multiplicity of $I$ are sufficiently close to each other. More precisely, for all $\varepsilon>0$, there exists $N\gg0$ such that for any $I\subseteq R$ with $l(R/I)>N$, we have $(1-\varepsilon)l(R/I)\leq e_{HK}(I)\leq(1+\varepsilon)l(R/I)$.


title: Liger: Combining Interaction Paradigms for Visual Analysis

abstract: Visualization tools usually leverage a single interaction paradigm (e.g., manual view specification, visualization by demonstration, etc.), which fosters the process of visualization construction. A large body of work has investigated the effectiveness of individual interaction paradigms, building an understanding of advantages and disadvantages of each in isolation. However, how can we leverage the benefits of multiple interaction paradigms by combining them into a single tool? We currently lack a holistic view of how interaction paradigms that use the same input modality (e.g., mouse) can be combined into a single tool and how people use such tools. To investigate opportunities and challenges in combining paradigms, we first created a multi-paradigm prototype (Liger) that combines two mouse-based interaction paradigms (manual view specification and visualization by demonstration) in a unified tool. We then conducted an exploratory study with Liger, providing initial evidence that people 1) use both paradigms interchangeably, 2) seamlessly switch between paradigms based on the operation at hand, and 3) choose to successfully complete a single operation using a combination of both paradigms.


title: Greedy Optimized Multileaving for Personalization

abstract: Personalization plays an important role in many services. To evaluate personalized rankings, online evaluation, such as A/B testing, is widely used today. Recently, multileaving has been found to be an efficient method for evaluating rankings in information retrieval fields. This paper describes the first attempt to optimize the multileaving method for personalization settings. We clarify the challenges of applying this method to personalized rankings. Then, to solve these challenges, we propose greedy optimized multileaving (GOM) with a new credit feedback function. The empirical results showed that GOM was stable for increasing ranking lengths and the number of rankers. We implemented GOM on our actual news recommender systems, and compared its online performance. The results showed that GOM evaluated the personalized rankings precisely, with significantly smaller sample sizes (< 1/10) than A/B testing.


title: Homoepitaxial growth of SrTiO$_3$ by Pulsed Laser Deposition: energetic vs thermal growth

abstract: Pulsed Laser Deposition (PLD) is widely used to grow epitaxial thin films of quantum materials such as complex oxides. Here, we use in-situ X-ray scattering to study homoepitaxy of SrTiO$_3$ by energetic (e-) and thermalized (th-) PLD. We find that e-PLD suppresses the lateral growth of two-dimensional islands, which suggests that energetic particles break up smaller islands. Fast interlayer transport occurs for both e-PLD and th-PLD, implying a process operating on sub-microsecond timescales that doesn't depend strongly on the kinetic energy of the incident particles.


title: Asymptotic resolvents of a product of two marginals of a random tensor

abstract: Random tensors can be used to produce random matrices. This idea is, for instance, very natural when one studies random quantum states with the aim of exploring properties that are generically true, or true with some probability. We hereby study the moments generating function, in the sense of the Stieltjes transform - i.e. the resolvent -, of a random matrix defined as a product of two different marginals of the same random tensor. We study the resolvent in two different asymptotical regimes. In the first regime, the resolvent is easily computed thanks to freeness results for the two different marginals and straightforward application of free harmonic analysis. In the second regime, we show that the resolvent satisfies an algebraic equation of degree six. This algebraic equation possesses roots whose expressions can be given explicitly in terms of radicals. We obtain this result by using an enumerative combinatorics approach. One of the interesting aspects of the second regime is that the corresponding probability density function interpolates between the square of a Marchenko-Pastur and the free multiplicative square of a Marchenko-Pastur law.


title: Convergence of Edge Computing and Deep Learning: A Comprehensive Survey

abstract: Ubiquitous sensors and smart devices from factories and communities are generating massive amounts of data, and ever-increasing computing power is driving the core of computation and services from the cloud to the edge of the network. As an important enabler broadly changing people's lives, from face recognition to ambitious smart factories and cities, developments of artificial intelligence (especially deep learning, DL) based applications and services are thriving. However, due to efficiency and latency issues, the current cloud computing service architecture hinders the vision of "providing artificial intelligence for every person and every organization at everywhere". Thus, unleashing DL services using resources at the network edge near the data sources has emerged as a desirable solution. Therefore, edge intelligence, aiming to facilitate the deployment of DL services by edge computing, has received significant attention. In addition, DL, as the representative technique of artificial intelligence, can be integrated into edge computing frameworks to build intelligent edge for dynamic, adaptive edge maintenance and management. With regard to mutually beneficial edge intelligence and intelligent edge, this paper introduces and discusses: 1) the application scenarios of both; 2) the practical implementation methods and enabling technologies, namely DL training and inference in the customized edge computing framework; 3) challenges and future trends of more pervasive and fine-grained intelligence. We believe that by consolidating information scattered across the communication, networking, and DL areas, this survey can help readers to understand the connections between enabling technologies while promoting further discussions on the fusion of edge intelligence and intelligent edge, i.e., Edge DL.


title: Spatially Aggregated Gaussian Processes with Multivariate Areal Outputs

abstract: We propose a probabilistic model for inferring the multivariate function from multiple areal data sets with various granularities. Here, the areal data are observed not at location points but at regions. Existing regression-based models can only utilize the sufficiently fine-grained auxiliary data sets on the same domain (e.g., a city). With the proposed model, the functions for respective areal data sets are assumed to be a multivariate dependent Gaussian process (GP) that is modeled as a linear mixing of independent latent GPs. Sharing of latent GPs across multiple areal data sets allows us to effectively estimate the spatial correlation for each areal data set; moreover it can easily be extended to transfer learning across multiple domains. To handle the multivariate areal data, we design an observation model with a spatial aggregation process for each areal data set, which is an integral of the mixed GP over the corresponding region. By deriving the posterior GP, we can predict the data value at any location point by considering the spatial correlations and the dependences between areal data sets, simultaneously. Our experiments on real-world data sets demonstrate that our model can 1) accurately refine coarse-grained areal data, and 2) offer performance improvements by using the areal data sets from multiple domains.


title: Heteroclinic solutions for a generalized Frenkel-Kontorova model by minimization methods of Rabinowitz and Stredulinsky

abstract: We study heteroclinic solutions of a generalized Frenkel-Kontorova model. Using the methods of Rabinowitz and Stredulinsky, we prove that if the rotation vector of the configuration is rational and if there is an adjacent pair of periodic configurations, then there is a solution that is heteroclinic in one fixed direction and periodic in other directions. Furthermore, if the above heteroclinic solutions have an adjacent pair, then there is a solution that is heteroclinic in two directions and periodic in other directions. The procedure can be repeated to produce more complex solutions. Thus we obtain a variational construction for these minimal and Birkhoff solutions.


title: Representation Learning for Classical Planning from Partially Observed Traces

abstract: Specifying a complete domain model is time-consuming, which has been a bottleneck of AI planning technique application in many real-world scenarios. Most classical domain-model learning approaches output a domain model in the form of the declarative planning language, such as STRIPS or PDDL, and solve new planning instances by invoking an existing planner. However, planning in such a representation is sensitive to the accuracy of the learned domain model which probably cannot be used to solve real planning problems. In this paper, to represent domain models in a vectorization representation way, we propose a novel framework based on graph neural network (GNN) integrating model-free learning and model-based planning, called LP-GNN. By embedding propositions and actions in a graph, the latent relationship between them is explored to form a domain-specific heuristics. We evaluate our approach on five classical planning domains, comparing with the classical domain-model learner ARMS. The experimental results show that the domain models learned by our approach are much more effective on solving real planning problems.


title: An extension of the Andrews-Warnaar partial theta function identity

abstract: In this paper, by applying a range of classic summation and transformation formulas for basic hypergeometric series, we obtain a three-term identity for partial theta functions. It extends the Andrews-Warnaar partial theta function identity, and also unifies several results on partial theta functions due to Ramanujan, Lovejoy and Kim. We also establish a two-term version of the extension, which can be used to derive identities for partial and false theta functions. Finally, we present a relation between the big $q$-Jacobi polynomials and the Andrews-Warnaar partial theta function identity.


title: Search for anomalous triple gauge couplings in WW and WZ production in lepton + jet events in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A search is presented for three additional operators that would lead to anomalous WW$\gamma$ or WWZ couplings with respect to those in the standard model. They are constrained by studying events with two vector bosons; a W boson decaying to e$\nu$ or $\mu\nu$, and a W or Z boson decaying hadronically, reconstructed as a single, massive, large-radius jet. The search uses a data set of proton-proton collisions at a centre-of-mass energy of 13 TeV, recorded by the CMS experiment at the CERN LHC in 2016, and corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Using the reconstructed diboson invariant mass, 95% confidence intervals are obtained for the anomalous coupling parameters of $-$1.58 $<$ $c_\mathrm{WWW}/\Lambda^2$ $<$ 1.59 TeV$^{-2}$, $-$2.00 $<$ $c_\mathrm{W}/\Lambda^2$ $<$ 2.65 TeV$^{-2}$, and $-$8.78 $<$ $c_\mathrm{B}/\Lambda^2$ $<$ 8.54 TeV$^{-2}$, in agreement with standard model expectations of zero for each parameter. These are the strictest bounds on these parameters to date.


title: Data Structures Meet Cryptography: 3SUM with Preprocessing

abstract: This paper shows several connections between data structure problems and cryptography against preprocessing attacks. Our results span data structure upper bounds, cryptographic applications, and data structure lower bounds, as summarized next. First, we apply Fiat--Naor inversion, a technique with cryptographic origins, to obtain a data structure upper bound. In particular, our technique yields a suite of algorithms with space $S$ and (online) time $T$ for a preprocessing version of the $N$-input 3SUM problem where $S^3\cdot T = \widetilde{O}(N^6)$. This disproves a strong conjecture (Goldstein et al., WADS 2017) that there is no data structure that solves this problem for $S=N^{2-\delta}$ and $T = N^{1-\delta}$ for any constant $\delta>0$. Secondly, we show equivalence between lower bounds for a broad class of (static) data structure problems and one-way functions in the random oracle model that resist a very strong form of preprocessing attack. Concretely, given a random function $F: [N] \to [N]$ (accessed as an oracle) we show how to compile it into a function $G^F: [N^2] \to [N^2]$ which resists $S$-bit preprocessing attacks that run in query time $T$ where $ST=O(N^{2-\varepsilon})$ (assuming a corresponding data structure lower bound on 3SUM). In contrast, a classical result of Hellman tells us that $F$ itself can be more easily inverted, say with $N^{2/3}$-bit preprocessing in $N^{2/3}$ time. We also show that much stronger lower bounds follow from the hardness of kSUM. Our results can be equivalently interpreted as security against adversaries that are very non-uniform, or have large auxiliary input, or as security in the face of a powerfully backdoored random oracle. Thirdly, we give non-adaptive lower bounds for 3SUM and a range of geometric problems which match the best known lower bounds for static data structure problems.


title: New Era of Deeplearning-Based Malware Intrusion Detection: The Malware Detection and Prediction Based On Deep Learning

abstract: With the development of artificial intelligence algorithms like deep learning models and the successful applications in many different fields, further similar trails of deep learning technology have been made in cyber security area. It shows the preferable performance not only in academic security research but also in industry practices when dealing with part of cyber security issues by deep learning methods compared to those conventional rules. Especially for the malware detection and classification tasks, it saves generous time cost and promotes the accuracy for a total pipeline of malware detection system. In this paper, we construct special deep neural network, ie, MalDeepNet (TB-Malnet and IB-Malnet) for malware dynamic behavior classification tasks. Then we build the family clustering algorithm based on deep learning and fulfil related testing. Except that, we also design a novel malware prediction model which could detect the malware coming in future through the Mal Generative Adversarial Network (Mal-GAN) implementation. All those algorithms present fairly considerable value in related datasets afterwards.


title: Charged current $b \to c \tau \bar{\nu}_\tau$ anomalies in a general $W^\prime$ boson scenario

abstract: Very recent experimental information obtained from the Belle experiment, along with that accumulated by the BABAR and LHCb experiments, has shown the existence of anomalies in the ratios $R(D)$ and $R(D^{*})$ associated with the charged-current transition $b \to c \tau \bar{\nu}_\tau$. Although the Belle measurements are in agreement with the standard model (SM) predictions, the new experimental world averages still exhibit a tension. In addition, the $D^\ast$ longitudinal polarization $F_L(D^\ast)$ related with the channel $B \to D^\ast \tau \bar{\nu}_\tau$ observed by the Belle Collaboration and the ratio $R(J/\psi)$ measured by the LHCb Collaboration also show discrepancies with their corresponding SM estimations. We present a model-independent study based on the most general effective Lagrangian that yields a tree-level effective contribution to the transition $b \to c \tau \bar{\nu}_\tau$ induced by a general $W^\prime$ boson. Instead of considering any specific new physics (NP) realization, we perform an analysis by considering all of the different chiral charges to the charm-bottom and $\tau$-$\nu_{\tau}$ interaction terms with a charged $W^\prime$ boson that explain the anomalies. We present a phenomenological study of parameter space allowed by the new experimental $b \to c \tau \bar{\nu}_\tau$ data and with the mono-tau signature $pp \to \tau_h X + \rm{MET}$ at the LHC. For comparison, we include some of the $W^\prime$ boson NP realizations that have already been studied in the literature.


title: Dark Matter Spin Characterisation in Mono-$Z$ Channels

abstract: The $B-L$ Supersymmetric Standard Model (BLSSM) is an ideal testing ground of the spin nature of Dark Matter (DM) as it offers amongst its candidates both a spin-1/2 (the lightest neutralino) and spin-0 (the lightest right-handed sneutrino) state. We show that the mono-$Z$ channel can be used at the Large Hadron Collider (LHC) to diagnose whether a DM signal is characterised within the BLSSM by a fermionic or (pseudo)scalar DM particle. Sensitivity to either hypothesis can be obtained after only 100 fb$^{-1}$ of luminosity following Runs 2 and 3 of the LHC.


title: First-principles prediction of a new family of layered topological insulators

abstract: Dozens of layered V$_2$IV$_2$VI$_6$ (V=P, As, Sb, Bi, IV=Si, Ge, Sn, Pb, VI=S, Se, Te) materials are investigated, several of which have been successfully synthesized in experiment. Among them, we predict nine strong topological insulators (TIs), two strong topological metals (TMs) and nearly twenty trivial insulators at their equilibrium structures. The TIs are in the (1;111) topological class, with energy gaps ranging from 0.04 to 0.2 eV. The strong TMs and the trivial insulators belong to the (1;111) and (0;000) topological classes, respectively. Small compressive strains easily turn some of the trivial insulators into strong TIs. This study enriches not only the family of topological materials but also the family of van der Waals layered materials, providing promising candidates for the future spintronic devices.


title: Matrix Completion for Survey Data Prediction with Multivariate Missingness

abstract: The National Health and Nutrition Examination Survey (NHANES) studies the nutritional and health status over the whole U.S. population with comprehensive physical examinations and questionnaires. However, survey data analyses become challenging due to inevitable missingness in almost all variables. In this paper, we develop a new imputation method to deal with multivariate missingness at random using matrix completion. In contrast to existing imputation schemes either conducting row-wise or column-wise imputation, we treat the data matrix as a whole which allows exploiting both row and column patterns to impute the missing values in the whole data matrix at one time. We adopt a column-space-decomposition model for the population data matrix with easy-to-obtain demographic data as covariates and a low-rank structured residual matrix. A unique challenge arises due to lack of identification of parameters in the sample data matrix. We propose a projection strategy to uniquely identify the parameters and corresponding penalized estimators, which are computationally efficient and possess desired statistical properties. The simulation study shows that the doubly robust estimator using the proposed matrix completion for imputation has smaller mean squared error than other competitors. To demonstrate practical relevance, we apply the proposed method to the 2015-2016 NHANES Questionnaire Data.


title: Surface Anchoring Energy of Cholesteric Liquid Crystals

abstract: In this paper, we propose a suitable surface energy expression for cholesteric liquid crystals. We show that there exists a symmetry allowed term for chiral nematics that doesn't appear in the traditional Rapini-Papoular surface energy form. We discuss some consequences of this new surface anchoring term.


title: Sparse Recovery for Orthogonal Polynomial Transforms

abstract: In this paper we consider the following sparse recovery problem. We have query access to a vector $\vx \in \R^N$ such that $\vhx = \vF \vx$ is $k$-sparse (or nearly $k$-sparse) for some orthogonal transform $\vF$. The goal is to output an approximation (in an $\ell_2$ sense) to $\vhx$ in sublinear time. This problem has been well-studied in the special case that $\vF$ is the Discrete Fourier Transform (DFT), and a long line of work has resulted in sparse Fast Fourier Transforms that run in time $O(k \cdot \mathrm{polylog} N)$. However, for transforms $\vF$ other than the DFT (or closely related transforms like the Discrete Cosine Transform), the question is much less settled. In this paper we give sublinear-time algorithms---running in time $\poly(k \log(N))$---for solving the sparse recovery problem for orthogonal transforms $\vF$ that arise from orthogonal polynomials. More precisely, our algorithm works for any $\vF$ that is an orthogonal polynomial transform derived from Jacobi polynomials. The Jacobi polynomials are a large class of classical orthogonal polynomials (and include Chebyshev and Legendre polynomials as special cases), and show up extensively in applications like numerical analysis and signal processing. One caveat of our work is that we require an assumption on the sparsity structure of the sparse vector, although we note that vectors with random support have this property with high probability. Our approach is to give a very general reduction from the $k$-sparse sparse recovery problem to the $1$-sparse sparse recovery problem that holds for any flat orthogonal polynomial transform; then we solve this one-sparse recovery problem for transforms derived from Jacobi polynomials.


title: Joint Coverage and Power Control in Highly Dynamic and Massive UAV Networks: An Aggregative Game-theoretic Learning Approach

abstract: Unmanned aerial vehicles (UAV) ad-hoc network is a significant contingency plan for communication after a natural disaster, such as typhoon and earthquake. To achieve efficient and rapid networks deployment, we employ noncooperative game theory and amended binary log-linear algorithm (BLLA) seeking for the Nash equilibrium which achieves the optimal network performance. We not only take channel overlap and power control into account but also consider coverage and the complexity of interference. However, extensive UAV game theoretical models show limitations in post-disaster scenarios which require large-scale UAV network deployments. Besides, the highly dynamic post-disaster scenarios cause strategies updating constraint and strategy-deciding error on UAV ad-hoc networks. To handle these problems, we employ aggregative game which could capture and cover those characteristics. Moreover, we propose a novel synchronous payoff-based binary log-linear learning algorithm (SPBLLA) to lessen information exchange and reduce time consumption. Ultimately, the experiments indicate that, under the same strategy-deciding error rate, SPBLLA's learning rate is manifestly faster than that of the revised BLLA. Hence, the new model and algorithm are more suitable and promising for large-scale highly dynamic scenarios.


title: Interlayer correlated fractional quantum Hall state in the {\nu} = 4/5 bilayer system

abstract: We perform exact diagonalization studies for fractional quantum Hall states at filing factor 4/5 in a bilayer system, on a torus with various aspect ratios and angles. We find that in the absence of tunneling, two weakly coupled 2/5-layers undergo a phase transition into an interlayer-correlated regime, which is also Abelian with the five-fold degeneracy on the torus. In the limit of zero layer separation, this phase becomes a singlet in the pseudospin variable describing the layer degree-of-freedom. By studying the Chern-number matrix, we show that the K-matrix describing the interlayer-correlated regime requires matrix dimension larger than two and this regime is in particular not described by a Halperin state. A detailed analysis of possible 4 X 4 K-matrices having the requisite symmetries and quantum numbers shows that there is only one equivalence class of such matrices. A model wave function representing this universality class is constructed. The role of separate particle number conservation in both layers is discussed, and it is argued that this additional symmetry allows for the further distinction of two different symmetry protected Abelian phases in the interlayer correlated regime. Interlayer tunneling breaks this symmetry, and can drive the system into a single-layer regime when strong enough. A qualitative phase diagram in the tunneling-layer separation parameter space is proposed based on our numerical results.


title: Existence of a conjugate point in the incompressible Euler flow on an ellipsoid

abstract: Existence of a conjugate point in the incompressible Euler flow on a sphere and an ellipsoid is considered. Misiolek (1996) formulated a differential-geometric criterion (we call M-criterion) for the existence of a conjugate point in a fluid flow. In this paper, it is shown that no zonal flow (stationary Euler flow) satisfies M-criterion if the background manifold is a sphere, on the other hand, there are zonal flows satisfy M-criterion if the background manifold is an ellipsoid (even it is sufficiently close to the sphere). The conjugate point is created by the fully nonlinear effect of the inviscid fluid flow with differential geometric mechanism.

categories: math.DG math.AP physics.flu-dyn physics.geo-ph


title: Components of Brill-Noether Loci for Curves with Fixed Gonality

abstract: We describe a conjectural stratification of the Brill-Noether variety for general curves of fixed genus and gonality. As evidence for this conjecture, we show that this Brill-Noether variety has at least as many irreducible components as predicted by the conjecture, and that each of these components has the expected dimension. Our proof uses combinatorial and tropical techniques. Specifically, we analyze containment relations between the various strata of tropical Brill-Noether loci identified by Pflueger in his classification of special divisors on chains of loops.


title: Optimizing Validation Phase of Hyperledger Fabric

abstract: Blockchain technologies are on the rise, and Hyperledger Fabric is one of the most popular permissioned blockchain platforms. In this paper, we re-architect the validation phase of Fabric based on our analysis from fine-grained breakdown of the validation phase's latency. Our optimized validation phase uses a chaincode cache during validation of transactions, initiates state database reads in parallel with validation of transactions, and writes to the ledger and databases in parallel. Our experiments reveal performance improvements of 2x for CouchDB and 1.3x for LevelDB. Notably, our optimizations can be adopted in a future release of Hyperledger Fabric.