title: Joint spectra of spherical Aluthge transforms of commuting n-tuples of Hilbert space operators

abstract: Let $\mathbf{T} \equiv (T_1,\cdots,T_n)$ be a commuting $n$-tuple of operators on a Hilbert space $\mathcal{H}$, and let $T_i \equiv V_i P \; (1 \le i \le n)$ be its canonical joint polar decomposition (i.e., $P:=\sqrt{T_1^*T_1+\cdots+T_n^*T_n}$, $(V_1,\cdots,V_n)$ a joint partial isometry, and $\bigcap_{i=1}^n \ker T_i = \bigcap_{i=1}^n \ker V_i = \ker P)$. \ The spherical Aluthge transform of $\mathbf{T}$ is the (necessarily commuting) $n$-tuple $\hat{\mathbf{T}}:=(\sqrt{P}V_1\sqrt{P},\cdots,\sqrt{P}V_n\sqrt{P})$. \ We prove that $\sigma_T(\hat{\mathbf{T}})=\sigma_T(\mathbf{T})$, where $\sigma_T$ denotes Taylor spectrum. \ We do this in two stages: away from the origin we use tools and techniques from criss-cross commutativity; at the origin we show that the left invertibility of $\mathbf{T}$ or $\hat{\mathbf{T}}$ implies the invertibility of $P$. \ As a consequence, we can readily extend our main result to other spectral systems that rely on the Koszul complex for their definitions.


title: Demonstration of a Reconfigurable Entangled Radiofrequency-Photonic Sensor Network

abstract: Quantum metrology takes advantage of nonclassical resources such as entanglement to achieve a sensitivity level below the standard quantum limit. To date, almost all quantum-metrology demonstrations are restricted to improving the measurement performance at a single sensor, but a plethora of applications require multiple sensors that work jointly to tackle distributed sensing problems. Here, we propose and experimentally demonstrate a reconfigurable sensor network empowered by continuous-variable (CV) multipartite entanglement. Our experiment establishes a connection between the entanglement structure and the achievable quantum advantage in different distributed sensing problems. The demonstrated entangled sensor network is composed of three sensor nodes each equipped with an electro-optic transducer for the detection of radiofrequency (RF) signals. By properly tailoring the CV multipartite entangled states, the entangled sensor network can be reconfigured to maximize the quantum advantage in distributed RF sensing problems such as measuring the angle of arrival of an RF field. The rich physics of CV multipartite entanglement unveiled by our work would open a new avenue for distributed quantum sensing and would lead to applications in ultrasensitive positioning, navigation, and timing.


title: Formulae for the Derivative of the Poincar\'e Constant of Gibbs Measures

abstract: We establish formulae for the derivative of the Poincar\'e constant of Gibbs measures on both compact domains and all of $\R^d$. As an application, we show that if the (not necessarily convex) Hamiltonian is an increasing function, then the Poincar\'e constant is strictly decreasing in the inverse temperature, and vice versa. Applying this result to the $O(2)$ model allows us to give a sharpened upper bound on its Poincar\'e constant. We further show that this model exhibits a qualitatively different zero-temperature behavior of the Poincar\'e and Log-Sobolev constants.


title: Quasinormality of powers of commuting pairs of bounded operators

abstract: We study jointly quasinormal and spherically quasinormal pairs of commuting operators on Hilbert space, as well as their powers. We first prove that, up to a constant multiple, the only jointly quasinormal $2$-variable weighted shift is the Helton-Howe shift. Second, we show that a left invertible subnormal operator $T$ whose square $T^{2}$ is quasinormal must be quasinormal. Third, we generalize a characterization of quasinormality for subnormal operators in terms of their normal extensions to the case of commuting subnormal $n$-tuples. Fourth, we show that if a $2$-variable weighted shift $W_{\left(\alpha ,\beta \right) }$ and its powers $W_{\left(\alpha ,\beta \right)}^{(2,1)}$ and $W_{\left(\alpha ,\beta \right)}^{(1,2)}$ are all spherically quasinormal, then $W_{\left( \alpha ,\beta \right)}$ may not necessarily be jointly quasinormal. Moreover, it is possible for both $W_{\left(\alpha ,\beta \right)}^{(2,1)}$ and $W_{\left(\alpha ,\beta \right)}^{(1,2)}$ to be spherically quasinormal without $W_{\left(\alpha ,\beta \right)}$ being spherically quasinormal. Finally, we prove that, for $2$-variable weighted shifts, the common fixed points of the toral and spherical Aluthge transforms are jointly quasinormal.


title: Dictionary Learning with Almost Sure Error Constraints

abstract: A dictionary is a database of standard vectors, so that other vectors / signals are expressed as linear combinations of dictionary vectors, and the task of learning a dictionary for a given data is to find a good dictionary so that the representation of data points has desirable features. Dictionary learning and the related matrix factorization methods have gained significant prominence recently due to their applications in Wide variety of fields like machine learning, signal processing, statistics etc. In this article we study the dictionary learning problem for achieving desirable features in the representation of a given data with almost sure recovery constraints. We impose the constraint that every sample is reconstructed properly to within a predefined threshold. This problem formulation is more challenging than the conventional dictionary learning, which is done by minimizing a regularised cost function. We make use of the duality results for linear inverse problems to obtain an equivalent reformulation in the form of a convex-concave min-max problem. The resulting min-max problem is then solved using gradient descent-ascent like algorithms.


title: The Tessellation-Level-Tree: characterising the nested hierarchy of density peaks and their spatial distribution in cosmological N-body simulations

abstract: We use the Millennium and Millennium-II simulations to illustrate the Tessellation-Level-Tree (TLT), a hierarchical tree structure linking density peaks in a field constructed by voronoi tessellation of the particles in a cosmological N-body simulation. The TLT uniquely partitions the simulation particles into disjoint subsets, each associated with a local density peak. Each peak is a subpeak of a unique higher peak. The TLT can be persistence filtered to suppress peaks produced by discreteness noise. Thresholding a peak's particle list at $\sim 80\left<\rho\right>$ results in a structure similar to a standard friend-of-friends halo and its subhaloes. For thresholds below $\sim 7\left<\rho\right>$, the largest structure percolates and is much more massive than other objects. It may be considered as defining the cosmic web. For a threshold of $5\left<\rho\right>$, it contains about half of all cosmic mass and occupies $\sim 1\%$ of all cosmic volume; a typical external point is then $\sim 7h^{-1}\mathrm{Mpc}$ from the web. We investigate the internal structure and clustering of TLT peaks. Defining the saddle point density $\rho_{\mathrm{lim}}$ as the density at which a peak joins its parent peak, we show the median value of $\rho_{\mathrm{lim}}$ for FoF-like peaks to be similar to the density threshold at percolation. Assembly bias as a function of $\rho_{\mathrm{lim}}$ is stronger than for any known internal halo property. For peaks of group mass and below, the lowest quintile in $\rho_{\mathrm{lim}}$ has $b\approx 0$, and is thus uncorrelated with the mass distribution.


title: Unitriangular Shape of Decomposition Matrices of Unipotent Blocks

abstract: We show that the decomposition matrix of unipotent $\ell$-blocks of a finite reductive group $\mathbf{G}(\mathbb{F}_q)$ has a unitriangular shape, assuming $q$ is a power of a good prime and $\ell$ is very good for $\mathbf{G}$. This was conjectured by Geck in 1990 as part of his PhD thesis. We establish this result by constructing projective modules using a modification of generalised Gelfand--Graev characters introduced by Kawanaka. We prove that each such character has at most one unipotent constituent which occurs with multiplicity one. This establishes a 30 year old conjecture of Kawanaka.


title: Novel Screening with Two Bodies: Summing the ladder in disformal scalar-tensor theories

abstract: When augmenting our cosmological models or gravitational theories with an additional light scalar field, any coupling between matter and this scalar can affect the orbital motion of binary systems. Ordinarily, the new force mediated by the scalar can be naturally the same order of magnitude as the usual gravitational force and therefore is tightly constrained. We show that a disformal coupling between the scalar and matter can lead to a novel screening mechanism in which these fifth forces are suppressed by several orders of magnitude at sufficiently small separations and large relative velocities (such as solar system scales). This is a result of resumming a class of ladder diagrams, which suppresses the propagation of scalar signals between the two bodies. Moreover, we are able to relate potential ambiguities in this resummation to non-perturbative effects (which are invisible to perturbation theory). As a result, solar system tests and future gravitational wave observations can now be used to place meaningful constraints on scalar-tensor theories with disformal couplings. We exemplify this using observational bounds on the precession of planetary orbits.


title: Natural Question Generation with Reinforcement Learning Based Graph-to-Sequence Model

abstract: Natural question generation (QG) aims to generate questions from a passage and an answer. In this paper, we propose a novel reinforcement learning (RL) based graph-to-sequence (Graph2Seq) model for QG. Our model consists of a Graph2Seq generator where a novel Bidirectional Gated Graph Neural Network is proposed to embed the passage, and a hybrid evaluator with a mixed objective combining both cross-entropy and RL losses to ensure the generation of syntactically and semantically valid text. The proposed model outperforms previous state-of-the-art methods by a large margin on the SQuAD dataset.


title: Plutos Far Side

abstract: The New Horizons spacecraft provided near global observations of Pluto that far exceed the resolution of Earth-based data sets. Most Pluto New Horizons analysis hitherto has focused on the encounter hemisphere of Pluto (i.e., the antiCharon hemisphere containing Sputnik Planitia). In this work, we summarize and interpret data on the far side (i.e., the non-encounter hemisphere), providing the first integrated New Horizons overview of the far side terrains. We find strong evidence for widespread bladed deposits, evidence for an impact crater about as large as any on the near side hemisphere, evidence for complex lineations approximately antipodal to Sputnik Planitia that may be causally related, and evidence that the far side maculae are smaller and more structured than the encounter hemisphere maculae.


title: Spectromicroscopic measurement of surface and bulk band structure interplay in a disordered topological insulator

abstract: Topological insulators are bulk semiconductors that manifest in-gap massless Dirac surface states due to the topological bulk-boundary correspondence principle [1-3]. These surface states have been a subject of tremendous ongoing interest, due both to their intrinsic properties and to higher order emergence phenomena that can be achieved by manipulating the interface environment [4-11]. Here, angle resolved photoemission (ARPES) spectromicroscopy and supplementary scanning tunneling microscopy (STM) are performed on the model topological insulator Bi2Se3 to investigate the interplay of crystallographic inhomogeneity with the topologically ordered bulk and surface band structure. Quantitative analysis methods are developed to obtain key spectroscopic information in spite of a limited dwell time on each measured point. Band energies are found to vary on the scale of 50 meV across the sample surface, enabling single-sample measurements that are analogous to a multi-sample doping series (termed a "binning series"). Focusing separately on the surface and bulk electrons reveals a nontrivial hybridization-like interplay between fluctuations in the surface and bulk state energetics.


title: Quantum break in high intensity gravitational wave interactions

abstract: The lowest order amplitudes for [$ {\rm graviton+graviton \rightarrow photon +photon}]$ lead to cross-sections of order $G^2$, where $G$ is the gravitational constant. These are too small to be of any interest. However in dense clouds of pure gravitons there are collective effects utilizing these same amplitudes that under the right circumstances can lead to copious production of photons.


title: Superintegrable Systems on 3 Dimensional Conformally Flat Spaces

abstract: We consider Hamiltonians associated with 3 dimensional conformally flat spaces, possessing 2, 3 and 4 dimensional isometry algebras. We use the conformal algebra to build additional {\em quadratic} first integrals, thus constructing a large class of superintegrable systems and the complete Poisson algebra of first integrals. We then use the isometries to reduce our systems to 2 degrees of freedom. For each isometry algebra we give a {\em universal} reduction of the corresponding general Hamiltonian. The superintegrable specialisations reduce, in this way, to systems of Darboux-Koenigs type, whose integrals are reductions of those of the 3 dimensional system.


title: Trans-Planckian Censorship and Inflation in Grand Unified Theories

abstract: The recently proposed trans-Planckian censorship conjecture (TCC) amounts to the claim that inflation models with an inflationary energy scale larger than Lambda_inf^max ~ 10^9 GeV belong to the swampland, i.e., cannot be embedded into a consistent theory of quantum gravity. In this paper, we point out that this constraint can be readily satisfied in D-term hybrid inflation (DHI), which is a well-motivated inflation scenario in the context of supersymmetric grand unification. In DHI, the amplitude of the primordial scalar power spectrum originates from a Fayet-Iliopoulos term of the order of the unification scale, sqrt{xi} ~ 10^16 GeV. At the same time, the TCC results in an upper bound on the corresponding gauge coupling constant of g_max ~ 10^-14. We are able to show that this constraint translates into an upper bound on the gravitino mass of m_3/2^max ~ 10 MeV, which opens the possibility that dark matter is accounted for by thermally produced gravitinos, if the reheating temperature is close to T_reh ~ 100 TeV. Interestingly enough, a somewhat similar gravitino mass range has recently been derived in a model that aims at explaining dark energy in terms of axion quintessence and resolving the Hubble tension by means of decaying gravitino dark matter.


title: Frequency-Sliding Generalized Cross-Correlation: A Sub-band Time Delay Estimation Approach

abstract: The generalized cross correlation (GCC) is regarded as the most popular approach for estimating the time difference of arrival (TDOA) between the signals received at two sensors. Time delay estimates are obtained by maximizing the GCC output, where the direct-path delay is usually observed as a prominent peak. Moreover, GCCs play also an important role in steered response power (SRP) localization algorithms, where the SRP functional can be written as an accumulation of the GCCs computed from multiple sensor pairs. Unfortunately, the accuracy of TDOA estimates is affected by multiple factors, including noise, reverberation and signal bandwidth. In this paper, a sub-band approach for time delay estimation aimed at improving the performance of the conventional GCC is presented. The proposed method is based on the extraction of multiple GCCs corresponding to different frequency bands of the cross-power spectrum phase in a sliding-window fashion. The major contributions of this paper include: 1) a sub-band GCC representation of the cross-power spectrum phase that, despite having a reduced temporal resolution, provides a more suitable representation for estimating the true TDOA; 2) such matrix representation is shown to be rank one in the ideal noiseless case, a property that is exploited in more adverse scenarios to obtain a more robust and accurate GCC; 3) we propose a set of low-rank approximation alternatives for processing the sub-band GCC matrix, leading to better TDOA estimates and source localization performance. An extensive set of experiments is presented to demonstrate the validity of the proposed approach.


title: A statistical mechanical model for drug release: relations between release parameters and porosity

abstract: A lattice gas model is proposed for investigating the release of drug molecules on devices with semi-permeable, porous membranes in two and three dimensions. The kinetic of this model was obtained through the analytical solution of the three-dimension diffusion equation for systems without membrane and with Monte Carlo simulations. Pharmaceutical data from drug release is usually adjusted to the Weibull function, $\exp [-(t/\tau)^b ]$, also known as stretched exponential, and the dependence of adjusted parameters $b$ and $\tau$ is usually associated, in the pharmaceutical literature, with physical mechanisms dominating the drug dynamics inside the capsule. The relation of parameters $\tau$ and $b$ with porosity $\lambda$ are found to satisfy, a simple linear relation for between $\tau$ and $\lambda^{-1}$, which can be explained through simple physically based arguments, and a scaling relation between $b$ and $\lambda$, with the scaling coefficient proportional to the system dimension.


title: Keyphrase Extraction from Scholarly Articles as Sequence Labeling using Contextualized Embeddings

abstract: In this paper, we formulate keyphrase extraction from scholarly articles as a sequence labeling task solved using a BiLSTM-CRF, where the words in the input text are represented using deep contextualized embeddings. We evaluate the proposed architecture using both contextualized and fixed word embedding models on three different benchmark datasets (Inspec, SemEval 2010, SemEval 2017) and compare with existing popular unsupervised and supervised techniques. Our results quantify the benefits of (a) using contextualized embeddings (e.g. BERT) over fixed word embeddings (e.g. Glove); (b) using a BiLSTM-CRF architecture with contextualized word embeddings over fine-tuning the contextualized word embedding model directly, and (c) using genre-specific contextualized embeddings (SciBERT). Through error analysis, we also provide some insights into why particular models work better than others. Lastly, we present a case study where we analyze different self-attention layers of the two best models (BERT and SciBERT) to better understand the predictions made by each for the task of keyphrase extraction.


title: Resilient Distributed Recovery of Large Fields

abstract: This paper studies the resilient distributed recovery of large fields under measurement attacks, by a team of agents, where each measures a small subset of the components of a large spatially distributed field. An adversary corrupts some of the measurements. The agents collaborate to process their measurements, and each is interested in recovering only a fraction of the field. We present a field recovery consensus+innovations type distributed algorithm that is resilient to measurement attacks, where an agent maintains and updates a local state based on its neighbors states and its own measurement. Under sufficient conditions on the attacker and the connectivity of the communication network, each agent's state, even those with compromised measurements, converges to the true value of the field components that it is interested in recovering. Finally, we illustrate the performance of our algorithm through numerical examples.


title: Machine Learning for AC Optimal Power Flow

abstract: We explore machine learning methods for AC Optimal Powerflow (ACOPF) - the task of optimizing power generation in a transmission network according while respecting physical and engineering constraints. We present two formulations of ACOPF as a machine learning problem: 1) an end-to-end prediction task where we directly predict the optimal generator settings, and 2) a constraint prediction task where we predict the set of active constraints in the optimal solution. We validate these approaches on two benchmark grids.


title: String C-group representations of almost simple groups: a survey

abstract: This survey paper aims at giving the state of the art in the study of string C-group representations of almost simple groups. It also suggest a series of problems and conjectures to the interested reader.


title: UW-SVC: Scalable Video Coding Transmission for In-network Underwater Imagery Analysis

abstract: Underwater imagery has enabled numerous civilian applications in various domains, ranging from academia to industry, and from industrial surveillance and maintenance to environmental protection and behavior of marine creatures studies. The accumulation of litter and plastic debris at the seafloor and the bottom of rivers are extremely harmful for the aquatic life. We propose a solution for this problem using a team of Autonomous Underwater Vehicles (AUVs) to exchange the recorded video in order to reconstruct the seafloor regions of interest. However, underwater video transmission is a challenge in the harsh environment in which radio-frequency waves are absorbed for distances above a few tens of meters, optical waves require narrow laser beams and suffer from scattering and ocean wave motion, and acoustic waves, while long range, provide a very low bandwidth and unreliable channel for communication. In our solution, the scalable coded video of each vehicle is shared with a selected group of receiving videos, pseudo-multicasting, through the acoustic channel. Presented evaluations, including both simulations and experiments, confirm the efficiency and flexibility of the proposed solution using acoustic software-defined modems.


title: ProxIQA: A Proxy Approach to Perceptual Optimization of Learned Image Compression

abstract: The use of $\ell_p$ $(p=1,2)$ norms has largely dominated the measurement of loss in neural networks due to their simplicity and analytical properties. However, when used to assess the loss of visual information, these simple norms are not very consistent with human perception. Here, we describe a different "proximal" approach to optimize image analysis networks against quantitative perceptual models. Specifically, we construct a proxy network, broadly termed ProxIQA, which mimics the perceptual model while serving as a loss layer of the network. We experimentally demonstrate how this optimization framework can be applied to train an end-to-end optimized image compression network. By building on top of an existing deep image compression model, we are able to demonstrate a bitrate reduction of as much as $31\%$ over MSE optimization, given a specified perceptual quality (VMAF) level.


title: Efficient Emulation of Computer Models Utilising Multiple Known Boundaries of Differing Dimensions

abstract: Emulation has been successfully applied across a wide variety of scientific disciplines for efficiently analysing computationally intensive models. We develop known boundary emulation strategies which utilise the fact that, for many computer models, there exist hyperplanes in the input parameter space for which the model output can be evaluated far more efficiently, whether this be analytically or just significantly faster using a more efficient and simpler numerical solver. The information contained on these known hyperplanes, or boundaries, can be incorporated into the emulation process via analytical update, thus involving no additional computational cost. In this article, we show that such analytical updates are available for multiple boundaries of various dimensions. We subsequently demonstrate which configurations of boundaries such analytical updates are available for, in particular by presenting a set of conditions that such a set of boundaries must satisfy. We demonstrate the powerful computational advantages of the known boundary emulation techniques developed on both an illustrative low-dimensional simulated example and a scientifically relevant and high-dimensional systems biology model of hormonal crosstalk in the roots of an Arabidopsis plant.


title: BUT System Description for DIHARD Speech Diarization Challenge 2019

abstract: This paper describes the systems developed by the BUT team for the four tracks of the second DIHARD speech diarization challenge. For tracks 1 and 2 the systems were based on performing agglomerative hierarchical clustering (AHC) over x-vectors, followed by the Bayesian Hidden Markov Model (HMM) with eigenvoice priors applied at x-vector level followed by the same approach applied at frame level. For tracks 3 and 4, the systems were based on performing AHC using x-vectors extracted on all channels.


title: Stability of tangent bundles on smooth toric Picard-rank-2 varieties and surfaces

abstract: We give a combinatorial criterion for the tangent bundle on a smooth toric variety to be stable with respect to a given polarisation in terms of the corresponding lattice polytope. Furthermore, we show that for a smooth toric surface and a smooth toric variety of Picard rank 2, there exists an ample line bundle with respect to which the tangent bundle is stable if and only if it is an iterated blow-up of projective space.


title: An Erd\H{o}s-Ko-Rado Theorem for unions of length 2 paths

abstract: A family of sets is intersecting if any two sets in the family intersect. Given a graph $G$ and an integer $r\geq 1$, let $\mathcal{I}^{(r)}(G)$ denote the family of independent sets of size $r$ of $G$. For a vertex $v$ of $G$, the family of independent sets of size $r$ that contain $v$ is called an $r$-star. Then $G$ is said to be $r$-EKR if no intersecting subfamily of $ \mathcal{I}^{(r)}(G)$ is bigger than the largest $r$-star. Let $n$ be a positive integer, and let $G$ consist of the disjoint union of $n$ paths each of length 2. We prove that if $1 \leq r \leq n/2$, then $G$ is $r$-EKR. This affirms a longstanding conjecture of Holroyd and Talbot for this class of graphs and can be seen as an analogue of a well-known theorem on signed sets, proved using different methods, by Deza and Frankl and by Bollob\'as and Leader. Our main approach is a novel probabilistic extension of Katona's elegant cycle method, which might be of independent interest.


title: H\"older parameterization of iterated function systems and a self-affine phenomenon

abstract: We investigate the H\"older geometry of curves generated by iterated function systems (IFS) in a complete metric space. A theorem of Hata from 1985 asserts that every connected attractor of an IFS is locally connected and path-connected. We give a quantitative strengthening of Hata's theorem. First we prove that every connected attractor of an IFS is $(1/s)$-H\"older path-connected, where $s$ is the similarity dimension of the IFS. Then we show that every connected attractor of an IFS is parameterized by a $(1/\alpha)$-H\"older curve for all $\alpha>s$. At the endpoint, $\alpha=s$, a theorem of Remes from 1998 already established that connected self-similar sets in Euclidean space that satisfy the open set condition are parameterized by $(1/s)$-H\"older curves. In a secondary result, we show how to promote Remes' theorem to self-similar sets in complete metric spaces, but in this setting require the attractor to have positive $s$-dimensional Hausdorff measure in lieu of the open set condition. To close the paper, we determine sharp H\"older exponents of parameterizations in the class of connected self-affine Bedford-McMullen carpets and build parameterizations of self-affine sponges. An interesting phenomenon emerges in the self-affine setting. While the optimal parameter $s$ for a self-similar curve in $\mathbb{R}^n$ is always at most the ambient dimension $n$, the optimal parameter $s$ for a self-affine curve in $\mathbb{R}^n$ may be strictly greater than $n$.


title: Atmospheric Neutrinos

abstract: Atmospheric neutrinos produced by cosmic-ray interactions around the globe provide a beam for the study of neutrino properties. They are also a background in searches for neutrinos of astrophysical origin. Both aspects are addressed in this chapter, which begins with a brief introduction on neutrino oscillations in relation to the spectrum of atmospheric neutrinos. Section 2 describes the cascade equation for hadrons in the atmosphere and the main features of atmospheric leptons from their decays. Next, uncertainties in the fluxes that arise from limited knowledge of the primary spectrum and of particle production are discussed. The final section covers aspects specific to neutrino telescopes.


title: An Analysis of the Shapes of Interstellar Extinction Curves. VII. Milky Way Spectrophotometric Optical-through-Ultraviolet Extinction and Its R-Dependence

abstract: We produce a set of 72 NIR through UV extinction curves by combining new HST/STIS optical spectrophotometry with existing IUE spectrophotometry (yielding gapless coverage from 1150 to 10000 Angstroms) and NIR photometry. These curves are used to determine a new, internally consistent, NIR through UV Milky Way Mean Curve and to characterize how the shapes of the extinction curves depend on R(V). We emphasize that while this dependence captures much of the curve variability, there remains considerable variation which is independent of R(V). We use the optical spectrophotometry to verify the presence of structure at intermediate wavelength scales in the curves. The fact that the optical through UV portions of the curves are sampled at relatively high resolution makes them very useful for determining how extinction affects different broad band systems, and we provide several examples. Finally, we compare our results to previous investigations.


title: Image Restoration Using Deep Regulated Convolutional Networks

abstract: While the depth of convolutional neural networks has attracted substantial attention in the deep learning research, the width of these networks has recently received greater interest. The width of networks, defined as the size of the receptive fields and the density of the channels, has demonstrated crucial importance in low-level vision tasks such as image denoising and restoration. However, the limited generalization ability, due to the increased width of networks, creates a bottleneck in designing wider networks. In this paper, we propose the Deep Regulated Convolutional Network (RC-Net), a deep network composed of regulated sub-network blocks cascaded by skip-connections, to overcome this bottleneck. Specifically, the Regulated Convolution block (RC-block), featured by a combination of large and small convolution filters, balances the effectiveness of prominent feature extraction and the generalization ability of the network. RC-Nets have several compelling advantages: they embrace diversified features through large-small filter combinations, alleviate the hazy boundary and blurred details in image denoising and super-resolution problems, and stabilize the learning process. Our proposed RC-Nets outperform state-of-the-art approaches with significant performance gains in various image restoration tasks while demonstrating promising generalization ability. The code is available at https://github.com/cswin/RC-Nets.


title: The Deepfake Detection Challenge (DFDC) Preview Dataset

abstract: In this paper, we introduce a preview of the Deepfakes Detection Challenge (DFDC) dataset consisting of 5K videos featuring two facial modification algorithms. A data collection campaign has been carried out where participating actors have entered into an agreement to the use and manipulation of their likenesses in our creation of the dataset. Diversity in several axes (gender, skin-tone, age, etc.) has been considered and actors recorded videos with arbitrary backgrounds thus bringing visual variability. Finally, a set of specific metrics to evaluate the performance have been defined and two existing models for detecting deepfakes have been tested to provide a reference performance baseline. The DFDC dataset preview can be downloaded at: deepfakedetectionchallenge.ai


title: A Recursion for the FiboNarayana and the Generalized Narayana Numbers

abstract: The Lucas polynomials, $\{n\}$, are polynomials in $s$ and $t$ given by $\{ n \} = s \{ n-1 \} + t \{ n-2 \}$ for $n \geq 2$ with $ \{ 0 \} = 0$ and $\{ 1 \} = 1$. The lucanomial coefficients, an analogue of the binomial coefficients, are given by \[ \Bigl\{ \begin{array}{c} n\\k \end{array} \Bigr \} = \frac{ \{n\}! }{ \{k\}! \{n-k\}!}. \] When $s = t = 1$ then $\{ n \} = F_n$ and the lucanomial coefficient becomes the fibonomial coefficient \[ \binom{n}{k}_F = \frac{F_n!}{F_k! F_{n-k}!}. \] The well-known Narayana numbers, $N_{n,k}$ satisfy the equation \[ N_{n,k} = \frac{1}{n} \binom{n}{k} \binom{n}{k-1}. \] \[ %C_n = \sum_{k=1}^n N_{n,k}. %\] In 2018, Bennett, Carrillo, Machacek and Sagan defined the generalized Narayana numbers and conjectured that these numbers are positive integers for $n \geq 1$. In this paper we define the FiboNarayana number $N_{n,k,F}$ and give a new recurrence relation for both the FiboNarayana numbers and the generalized Narayana numbers, proving the conjecture that these are positive integers for $n \geq 1$.


title: Multiscale-Spectral GFEM and Optimal Oversampling

abstract: In this work we address the Multiscale Spectral Generalized Finite Element Method (MS-GFEM) developed in [I. Babu\v{s}ka and R. Lipton, Multiscale Modeling and Simulation 9 (2011), pp. 373--406]. We outline the numerical implementation of this method and present simulations that demonstrate contrast independent exponential convergence of MS-GFEM solutions. We introduce strategies to reduce the computational cost of generating the optimal oversampled local approximating spaces used here. These strategies retain accuracy while reducing the computational work necessary to generate local bases. Motivated by oversampling we develop a nearly optimal local basis based on a partition of unity on the boundary and the associated A-harmonic extensions.


title: LRP2020: Astrostatistics in Canada

abstract: (Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada.


title: Beating the House: Identifying Inefficiencies in Sports Betting Markets

abstract: Inefficient markets allow investors to consistently outperform the market. To demonstrate that inefficiencies exist in sports betting markets, we created a betting algorithm that generates above market returns for the NFL, NBA, NCAAF, NCAAB, and WNBA betting markets. To formulate our betting strategy, we collected and examined a novel dataset of bets, and created a non-parametric win probability model to find positive expected value situations. As the United States Supreme Court has recently repealed the federal ban on sports betting, research on sports betting markets is increasingly relevant for the growing sports betting industry.


title: Microwave Photonic Exciter Unit for Radar System

abstract: Modern radar systems require low phase noise and long term phase-stable radar carrier signals for high-resolution imaging, micro-Doppler signatures applications. High stable signal generation through conventional crystal followed by frequency multiplier techniques are limited phase noise performance under vibration conditions. Optoelectronic oscillators (OEOs) offer low phase noise and long-term phase stability compared to the conventional oscillators. In this design of the Radar Exciter unit based on microwave photonics components viz., OEOs, frequency dividers, optical filters, optical delay lines, optical arbitrary waveform generators are discussed. The microwave photonics offers frequency-independent Exciter design with tunable filters and tunable laser sources. In this paper the performance comparison of microwave photonic-based exciter unit also discussed.


title: Distributed High Speed Optical Network for Digital Radar Systems

abstract: Modern Digital radar systems with multiple digital beamforming capability are built of a large number of receivers and requires high-speed data interface links for transmission of receiver baseband data to processor units. High data throughput (>250Mbyte/sec) from typical eight-channel receivers will be transmitted to Digital beamformer over high-speed serial interface links over an optical channel. Currently for digital radar systems with sub-array level beam former distribution of receiver data is through point to point optical interface links. For the modern element-level digital beamforming radars the distribution of baseband data increases the design complexity. In this paper novel scheme of usage of distributed optical interface network is discussed using high-speed optical transport networks, FPGAs as well as distributed techniques to over the above problem. The recent advances in optical communication and feature of FPGA devices are utilized in the implementation of optical distribution networks and these schemes are covered in this paper.


title: A Novel Scheme of Digital Instantaneous Automatic Gain Control (DIAGC) for Pulse Radars

abstract: Several schemes for gain control are used for preventing saturation of receiver, and overloading of data processor, tracker or display in pulse radars. The use of digital processing techniques open the door to a variety of digital automatic gain control schemes for analyzing digitized return signals and controlling receiver gain only at saturating clutter zones without affecting the detection at other zones. In this paper, we present a novel scheme of Digital Instantaneous Automatic Gain Control (DIAGC) which is based on storing digitally the dwell based clutter returns and deriving the gain control. The returns corresponding to the first two PRTs in a dwell are used to analyze the presence of saturating clutter zones and the depth of saturation. Third PRT onwards proper gain control is applied at the IF stage to prevent saturation of the following stages. FPGA based scheme is used for digital data processing, storing, threshold calculation and gain control generation. The effect of DIAGC on pulse compression is also addressed in this paper.


title: LSTM-Assisted Evolutionary Self-Expressive Subspace Clustering

abstract: Massive volumes of high-dimensional data that evolves over time is continuously collected by contemporary information processing systems, which brings up the problem of organizing this data into clusters, i.e. achieve the purpose of dimensional deduction, and meanwhile learning its temporal evolution patterns. In this paper, a framework for evolutionary subspace clustering, referred to as LSTM-ESCM, is introduced, which aims at clustering a set of evolving high-dimensional data points that lie in a union of low-dimensional evolving subspaces. In order to obtain the parsimonious data representation at each time step, we propose to exploit the so-called self-expressive trait of the data at each time point. At the same time, LSTM networks are implemented to extract the inherited temporal patterns behind data in an overall time frame. An efficient algorithm has been proposed based on MATLAB. Next, experiments are carried out on real-world datasets to demonstrate the effectiveness of our proposed approach. And the results show that the suggested algorithm dramatically outperforms other known similar approaches in terms of both run time and accuracy.


title: Reproducibility of high-performance quantum dot single-photon sources

abstract: Single-photon sources based on semiconductor quantum dots have emerged as an excellent platform for high efficiency quantum light generation. However, scalability remains a challenge since quantum dots generally present inhomogeneous characteristics. Here we benchmark the performance of fifteen deterministically fabricated single-photon sources. They display an average indistinguishability of 90.6 +/- 2.8 % with a single-photon purity of 95.4 +/- 1.5 % and high homogeneity in operation wavelength and temporal profile. Each source also has state-of-the-art brightness with an average first lens brightness value of 13.6 +/- 4.4 %. Whilst the highest brightness is obtained with a charged quantum dot, the highest quantum purity is obtained with neutral ones. We also introduce various techniques to identify the nature of the emitting state. Our study sets the groundwork for large-scale fabrication of identical sources by identifying the remaining challenges and outlining solutions.


title: Identification of Interaction Clusters Using a Semi-supervised Hierarchical Clustering Method

abstract: Motivation: Identifying interaction clusters of large gene regulatory networks (GRNs) is critical for its further investigation, while this task is very challenging, attributed to data noise in experiment data, large scale of GRNs, and inconsistency between gene expression profiles and function modules, etc. It is promising to semi-supervise this process by prior information, but shortage of prior information sometimes make it very challenging. Meanwhile, it is also annoying, and sometimes impossible to discovery gold standard for evaluation of clustering results.\\ Results: With assistance of an online enrichment tool, this research proposes a semi-supervised hierarchical clustering method via deconvolved correlation matrix~(SHC-DC) to discover interaction clusters of large-scale GRNs. Three benchmark networks including a \emph{Ecoli} network and two \emph{Yeast} networks are employed to test semi-supervision scheme of the proposed method. Then, SHC-DC is utilized to cluster genes in sleep study. Results demonstrates it can find interaction modules that are generally enriched in various signal pathways. Besides the significant influence on blood level of interleukins, impact of sleep on important pathways mediated by them is also validated by the discovered interaction modules.


title: Deck.gl: Large-scale Web-based Visual Analytics Made Easy

abstract: In this paper, we demonstrate how deck.gl, an open-source project born out of data-heavy visual analytics applications, has grown into the robust visualization framework it is today. We begin by explaining why we built another data visualization framework in the first place. Then, we summarize our design goals (distilled from our interactions with users) and discuss how they guided the development of the framework's main features. We use two real-world applications of deck.gl to showcase how it can be applied to simplify the creation of data-heavy visualizations. We also discuss our lessons learned as we continue to improve the framework for the larger visualization community.


title: Fusion Rings and Modular Invariance for Affine Lie Algebras: A Double Affine Weyl Group Formulation

abstract: From a certain induced representation $\mathcal{P}_\ell$ of a double affine Weyl group, we construct a ring $\mathcal{F}_\ell$ that is isomorphic to the fusion ring, or Verlinde algebra, associated to affine Lie algebras at fixed positive integer level for both twisted and untwisted type. The induced representation, which also has a natural commutative associative algebra structure and is modular invariant with respect to certain congruence subgroups, contains $\mathcal{F}_\ell$ as an ideal and we show how it naturally inherits the modular invariance property from $\mathcal{P}_\ell.$ This construction directly shows how the action of the modular transformation $S:\tau\to -1/\tau$ determines the structure constants, with respect to a natural basis, of $\mathcal{F}_\ell,$ which are precisely the fusion rules of Verlinde algebras. Using this ideal, we also give a simple proof of a well-known epimorphism between the representation rings of simple Lie algebras and the fusion rings of the corresponding affine Lie algebras.


title: KRNET: Image Denoising with Kernel Regulation Network

abstract: One popular strategy for image denoising is to design a generalized regularization term that is capable of exploring the implicit prior underlying data observation. Convolutional neural networks (CNN) have shown the powerful capability to learn image prior information through a stack of layers defined by a combination of kernels (filters) on the input. However, existing CNN-based methods mainly focus on synthetic gray-scale images. These methods still exhibit low performance when tackling multi-channel color image denoising. In this paper, we optimize CNN regularization capability by developing a kernel regulation module. In particular, we propose a kernel regulation network-block, referred to as KR-block, by integrating the merits of both large and small kernels, that can effectively estimate features in solving image denoising. We build a deep CNN-based denoiser, referred to as KRNET, via concatenating multiple KR-blocks. We evaluate KRNET on additive white Gaussian noise (AWGN), multi-channel (MC) noise, and realistic noise, where KRNET obtains significant performance gains over state-of-the-art methods across a wide spectrum of noise levels.


title: Downlink Performance of Dense Antenna Deployment: To Distribute or Concentrate?

abstract: Massive multiple-input multiple-output (massive MIMO) and small cell densification are complementary key 5G enablers. Given a fixed number of the entire base-station antennas per unit area, this paper fairly compares (i) to deploy few base stations (BSs) and concentrate many antennas on each of them, i.e. massive MIMO, and (ii) to deploy more BSs equipped with few antennas, i.e. small cell densification. We observe that small cell densification always outperforms for both signal-to-interference ratio (SIR) coverage and energy efficiency (EE), when each BS serves multiple users via L number of sub-bands (multi-carrier transmission). Moreover, we also observe that larger L increases SIR coverage while decreasing EE, thus urging the necessity of optimal 5G network design. These two observations are based on our novel closed-form SIR coverage probability derivation using stochastic geometry, also validated via numerical simulations.


title: Eigenvalues and Spectral Dimension of Random Geometric Graphs in Thermodynamic Regime

abstract: Network geometries are typically characterized by having a finite spectral dimension (SD), $d_{s}$ that characterizes the return time distribution of a random walk on a graph. The main purpose of this work is to determine the SD of a variety of random graphs called random geometric graphs (RGGs) in the thermodynamic regime, in which the average vertex degree is constant. The spectral dimension depends on the eigenvalue density (ED) of the RGG normalized Laplacian in the neighborhood of the minimum eigenvalues. In fact, the behavior of the ED in such a neighborhood characterizes the random walk. Therefore, we first provide an analytical approximation for the eigenvalues of the regularized normalized Laplacian matrix of RGGs in the thermodynamic regime. Then, we show that the smallest non zero eigenvalue converges to zero in the large graph limit. Based on the analytical expression of the eigenvalues, we show that the eigenvalue distribution in a neighborhood of the minimum value follows a power-law tail. Using this result, we find that the SD of RGGs is approximated by the space dimension $d$ in the thermodynamic regime.


title: Critical behavior for a semilinear parabolic equation with forcing term depending of time and space

abstract: We investigate the large-time behavior of the sign-changing solution of the inhomogeneous semilinear heat equation with a forcing term depending of time and space. we identify the critical exponent for this problem, which separates the nonexistence/existence of global-in-time solutions, and show the discontinuity of this critical exponent at the point which the inhomogeneous term becomes independent of time variable.


title: Spectral Analysis of the Adjacency Matrix of Random Geometric Graphs

abstract: In this article, we analyze the limiting eigenvalue distribution (LED) of random geometric graphs (RGGs). The RGG is constructed by uniformly distributing $n$ nodes on the $d$-dimensional torus $\mathbb{T}^d \equiv [0, 1]^d$ and connecting two nodes if their $\ell_{p}$-distance, $p \in [1, \infty]$ is at most $r_{n}$. In particular, we study the LED of the adjacency matrix of RGGs in the connectivity regime, in which the average vertex degree scales as $\log\left( n\right)$ or faster, i.e., $\Omega \left(\log(n) \right)$. In the connectivity regime and under some conditions on the radius $r_{n}$, we show that the LED of the adjacency matrix of RGGs converges to the LED of the adjacency matrix of a deterministic geometric graph (DGG) with nodes in a grid as $n$ goes to infinity. Then, for $n$ finite, we use the structure of the DGG to approximate the eigenvalues of the adjacency matrix of the RGG and provide an upper bound for the approximation error.


title: Principal specializations of Schubert polynomials and pattern containment

abstract: We show that the principal specialization of the Schubert polynomial at $w$ is bounded below by $1+p_{132}(w)+p_{1432}(w)$ where $p_u(w)$ is the number of occurrences of the pattern $u$ in $w$, strengthening a previous result by A. Weigandt. We then make a conjecture relating the principal specialization of Schubert polynomials to pattern containment. Finally, we characterize permutations $w$ whose RC-graphs are connected by simple ladder moves via pattern avoidance.


title: Spectral bounds of the regularized normalized Laplacian for random geometric graphs

abstract: In this work, we study the spectrum of the regularized normalized Laplacian for random geometric graphs (RGGs) in both the connectivity and thermodynamic regimes. We prove that the limiting eigenvalue distribution (LED) of the normalized Laplacian matrix for an RGG converges to the Dirac measure in one in the full range of the connectivity regime. In the thermodynamic regime, we propose an approximation for the LED and we provide a bound on the Levy distance between the approximation and the actual distribution. In particular, we show that the LED of the regularized normalized Laplacian matrix for an RGG can be approximated by the LED of the regularized normalized Laplacian for a deterministic geometric graph with nodes in a grid (DGG). Thereby, we obtain an explicit approximation of the eigenvalues in the thermodynamic regime.


title: Speech Emotion Recognition with Dual-Sequence LSTM Architecture

abstract: Speech Emotion Recognition (SER) has emerged as a critical component of the next generation human-machine interfacing technologies. In this work, we propose a new dual-level model that predicts emotions based on both MFCC features and mel-spectrograms produced from raw audio signals. Each utterance is preprocessed into MFCC features and two mel-spectrograms at different time-frequency resolutions. A standard LSTM processes the MFCC features, while a novel LSTM architecture, denoted as Dual-Sequence LSTM (DS-LSTM), processes the two mel-spectrograms simultaneously. The outputs are later averaged to produce a final classification of the utterance. Our proposed model achieves, on average, a weighted accuracy of 72.7% and an unweighted accuracy of 73.3%---a 6% improvement over current state-of-the-art unimodal models---and is comparable with multimodal models that leverage textual information as well as audio signals.


title: Integrating DFT and DRBD Formalizations in HOL4

abstract: Dynamic Fault Trees (DFT) and Dynamic Reliability Block Diagrams (DRBD) are two modeling approaches that capture the dynamic failure behavior of engineering systems for their reliability analysis. Recently, two independent higher-order logic (HOL) formalizations of DFT and DRBD algebras have been developed in the HOL4 theorem prover. In this work, we propose to integrate these two modeling approaches for the efficient formal reliability analysis of complex systems by leveraging upon the advantages of each method. The soundness of this integration is provided through a formal proof of equivalence between the DFT and DRBD algebras. We show the efficiency of the proposed integrated formal reliability analysis on a drive-by-wire system as a case study.


title: Release Practices for Mobile Apps--What do Users and Developers Think?

abstract: Large software organizations such as Facebook or Netflix, who otherwise make daily or even hourly releases of their web applications using continuous delivery, have had to invest heavily into a customized release strategy for their mobile apps, because the vetting process of app stores introduces lag and uncertainty into the release process. Amidst these large, resourceful organizations, it is unknown how the average mobile app developer organizes her app's releases, even though an incorrect strategy might bring a premature app update to the market that drives away customers towards the heavy market competition. To understand the common release strategies used for mobile apps, the rationale behind them and their perceived impact on users, we performed two surveys with users and developers. We found that half of the developers have a clear strategy for their mobile app releases, since especially the more experienced developers believe that it affects user feedback. We also found that users are aware of new app updates, yet only half of the surveyed users enables automatic updating of apps. While the release date and frequency is not a decisive factor to install an app, users prefer to install apps that were updated more recently and less frequently. Our study suggests that an app's release strategy is a factor that affects the ongoing success of mobile apps.


title: Targeted Estimation of Heterogeneous Treatment Effect in Observational Survival Analysis

abstract: The aim of clinical effectiveness research using repositories of electronic health records is to identify what health interventions 'work best' in real-world settings. Since there are several reasons why the net benefit of intervention may differ across patients, current comparative effectiveness literature focuses on investigating heterogeneous treatment effect and predicting whether an individual might benefit from an intervention. The majority of this literature has concentrated on the estimation of the effect of treatment on binary outcomes. However, many medical interventions are evaluated in terms of their effect on future events, which are subject to loss to follow-up. In this study, we describe a framework for the estimation of heterogeneous treatment effect in terms of differences in time-to-event (survival) probabilities. We divide the problem into three phases: (1) estimation of treatment effect conditioned on unique sets of the covariate vector; (2) identification of features important for heterogeneity using an ensemble of non-parametric variable importance methods; and (3) estimation of treatment effect on the reference classes defined by the previously selected features, using one-step Targeted Maximum Likelihood Estimation. We conducted a series of simulation studies and found that this method performs well when either sample size or event rate is high enough and the number of covariates contributing to the effect heterogeneity is moderate. An application of this method to a clinical case study was conducted by estimating the effect of oral anticoagulants on newly diagnosed non-valvular atrial fibrillation patients using data from the UK Clinical Practice Research Datalink.


title: Probabilistic Radiomics: Ambiguous Diagnosis with Controllable Shape Analysis

abstract: Radiomics analysis has achieved great success in recent years. However, conventional Radiomics analysis suffers from insufficiently expressive hand-crafted features. Recently, emerging deep learning techniques, e.g., convolutional neural networks (CNNs), dominate recent research in Computer-Aided Diagnosis (CADx). Unfortunately, as black-box predictors, we argue that CNNs are "diagnosing" voxels (or pixels), rather than lesions; in other words, visual saliency from a trained CNN is not necessarily concentrated on the lesions. On the other hand, classification in clinical applications suffers from inherent ambiguities: radiologists may produce diverse diagnosis on challenging cases. To this end, we propose a controllable and explainable {\em Probabilistic Radiomics} framework, by combining the Radiomics analysis and probabilistic deep learning. In our framework, 3D CNN feature is extracted upon lesion region only, then encoded into lesion representation, by a controllable Non-local Shape Analysis Module (NSAM) based on self-attention. Inspired from variational auto-encoders (VAEs), an Ambiguity PriorNet is used to approximate the ambiguity distribution over human experts. The final diagnosis is obtained by combining the ambiguity prior sample and lesion representation, and the whole network named $DenseSharp^{+}$ is end-to-end trainable. We apply the proposed method on lung nodule diagnosis on LIDC-IDRI database to validate its effectiveness.


title: Classification of complex hyperbolic triangle groups by types

abstract: We give a complete classification of complex hyperbolic $(n_1, n_2, n_3)$-triangle groups by types defined according to the ellipticity of two particular words of short length. This improves the Schwartz conjecture proved by Grossi.


title: Improved error rates for sparse (group) learning with Lipschitz loss functions

abstract: We study a family of sparse estimators defined as minimizers of some empirical Lipschitz loss function -- which include the hinge loss, the logistic loss and the quantile regression loss -- with a convex, sparse or group-sparse regularization. In particular, we consider the L1 norm on the coefficients, its sorted Slope version, and the Group L1-L2 extension. We propose a new theoretical framework that uses common assumptions in the literature to simultaneously derive new high-dimensional L2 estimation upper bounds for all three regularization schemes. %, and to improve over existing results. For L1 and Slope regularizations, our bounds scale as $(k^*/n) \log(p/k^*)$ -- $n\times p$ is the size of the design matrix and $k^*$ the dimension of the theoretical loss minimizer $\B{\beta}^*$ -- and match the optimal minimax rate achieved for the least-squares case. For Group L1-L2 regularization, our bounds scale as $(s^*/n) \log\left( G / s^* \right) + m^* / n$ -- $G$ is the total number of groups and $m^*$ the number of coefficients in the $s^*$ groups which contain $\B{\beta}^*$ -- and improve over the least-squares case. We show that, when the signal is strongly group-sparse, Group L1-L2 is superior to L1 and Slope. In addition, we adapt our approach to the sub-Gaussian linear regression framework and reach the optimal minimax rate for Lasso, and an improved rate for Group-Lasso. Finally, we release an accelerated proximal algorithm that computes the nine main convex estimators of interest when the number of variables is of the order of $100,000s$.


title: EQSA: Earthquake Situational Analytics from Social Media

abstract: This paper introduces EQSA, an interactive exploratory tool for earthquake situational analytics using social media. EQSA is designed to support users to characterize the condition across the area around the earthquake zone, regarding related events, resources to be allocated, and responses from the community. On the general level, changes in the volume of messages from chosen categories are presented, assisting users in conveying a general idea of the condition. More in-depth analysis is provided with topic evolution, community visualization, and location representation. EQSA is developed with intuitive, interactive features and multiple linked views, visualizing social media data, and supporting users to gain a comprehensive insight into the situation. In this paper, we present the application of EQSA with the VAST Challenge 2019: Mini-Challenge 3 (MC3) dataset.


title: Classical skew orthogonal polynomials in a two-component log-gas with charges $+1$ and $+2$

abstract: There is a two-component log-gas system with Boltzmann factor which provides an interpolation between the eigenvalue PDF for $\beta = 1$ and $\beta = 4$ invariant random matrix ensembles. The solvability of this log-gas system relies on the construction of particular skew orthogonal polynomials, with the skew inner product a linear combination of the $\beta = 1$ and $\beta = 4$ inner products, each involving weight functions. For suitably related classical weight functions, we seek to express the skew orthogonal polynomials as linear combinations of the underlying orthogonal polynomials. It is found that in each case (Gaussian, Laguerre, Jacobi and generalised Cauchy) the coefficients can be expressed in terms of hypergeometric polynomials with argument relating to the fugacity. In the Jacobi case, for example, these are a special case of the Wilson polynomials.


title: High-dimensional and universally consistent k-sample tests

abstract: The k-sample testing problem involves determining whether $k$ groups of data points are each drawn from the same distribution. The standard method for k-sample testing in biomedicine is Multivariate analysis of variance (MANOVA), despite that it depends on strong, and often unsuitable, parametric assumptions. Moreover, independence testing and k-sample testing are closely related, and several universally consistent high-dimensional independence tests such as distance correlation (Dcorr) and Hilbert-Schmidt-Independence-Criterion (Hsic) enjoy solid theoretical and empirical properties. In this paper, we prove that independence tests achieve universally consistent k-sample testing and that k-sample statistics such as Energy and Maximum Mean Discrepancy (MMD) are precisely equivalent to Dcorr. An empirical evaluation of nonparametric independence tests showed that they generally perform better than the popular MANOVA test, even in Gaussian distributed scenarios. The evaluation included several popular independence statistics and covered a comprehensive set of simulations. Additionally, the testing approach was extended to perform multiway and multilevel tests, which were demonstrated in a simulated study as well as a real-world fMRI brain scans with a set of attributes.


title: A Truly Exact and Optimal Perfect Absorbing Layer for Time-harmonic Acoustic Wave Scattering Problems

abstract: In this paper, we design a truly exact and optimal perfect absorbing layer (PAL) for domain truncation of the two-dimensional Helmholtz equation in an unbounded domain with bounded scatterers. This technique is based on a complex compression coordinate transformation in polar coordinates, and a judicious substitution of the unknown field in the artificial layer. Compared with the widely-used perfectly matched layer (PML) methods, the distinctive features of PAL lie in that (i) it is truly exact in the sense that the PAL-solution is identical to the original solution in the bounded domain reduced by the truncation layer; (ii) with the substitution, the PAL-equation is free of singular coefficients and the substituted unknown field is essentially non-oscillatory in the layer; and (iii) the construction is valid for general star-shaped domain truncation. By formulating the variational formulation in Cartesian coordinates, the implementation of this technique using standard spectral-element or finite-element methods can be made easy as a usual coding practice. We provide ample numerical examples to demonstrate that this method is highly accurate, parameter-free and robust for very high wave-number and thin layer. It outperforms the classical PML and the recently advocated PML using unbounded absorbing functions. Moreover, it can fix some flaws of the PML approach.


title: Convex co-compact actions of relatively hyperbolic groups

abstract: In this paper we consider discrete groups in ${\rm PGL}_d(\mathbb{R})$ acting convex co-compactly on a properly convex domain in real projective space. For such groups, we establish necessary and sufficient conditions for the group to be relatively hyperbolic in terms of the geometry of the convex domain. This answers a question of Danciger-Gu\'eritaud-Kassel and is analogous to a result of Hruska-Kleiner for ${\rm CAT}(0)$ spaces.


title: Speech-Based Parameter Estimation of an Asymmetric Vocal Fold Oscillation Model and Its Application in Discriminating Vocal Fold Pathologies

abstract: So far, several physical models have been proposed for the study of vocal fold oscillations during phonation. The parameters of these models, such as vocal fold elasticity, resistance, etc. are traditionally determined through the observation and measurement of the vocal fold vibrations in the larynx. Since such direct measurements tend to be the most accurate, the traditional practice has been to set the parameter values of these models based on measurements that are averaged across an ensemble of human subjects. However, the direct measurement process is hard to revise outside of clinical settings. In many cases, especially in pathological ones, the properties of the vocal folds often deviate from their generic values---sometimes asymmetrically wherein the characteristics of the two vocal folds differ for the same individual. In such cases, it is desirable to find a more scalable way to adjust the model parameters on a case by case basis. In this paper, we present a novel and alternate way to determine vocal fold model parameters from the speech signal. We focus on an asymmetric model and show that for such models, differences in estimated parameters can be successfully used to discriminate between voices that are characteristic of different underlying vocal fold pathologies.


title: Personalized Graph Neural Networks with Attention Mechanism for Session-Aware Recommendation

abstract: The problem of session-aware recommendation aims to predict users' next click based on their current session and historical sessions. Existing session-aware recommendation methods have defects in capturing complex item transition relationships. Other than that, most of them fail to explicitly distinguish the effects of different historical sessions on the current session. To this end, we propose a novel method, named Personalized Graph Neural Networks with Attention Mechanism (A-PGNN) for brevity. A-PGNN mainly consists of two components: one is Personalized Graph Neural Network (PGNN), which is used to extract the personalized structural information in each user behavior graph, compared with the traditional Graph Neural Network (GNN) model, which considers the role of the user when the node embeddding is updated. The other is Dot-Product Attention mechanism, which draws on the Transformer net to explicitly model the effect of historical sessions on the current session. Extensive experiments conducted on two real-world data sets show that A-PGNN evidently outperforms the state-of-the-art personalized session-aware recommendation methods.


title: Monotonic Properties of Completed Aggregates in Recursive Queries

abstract: The use of aggregates in recursion enables efficient and scalable support for a wide range of BigData algorithms, including those used in graph applications, KDD applications, and ML applications, which have proven difficult to be expressed and supported efficiently in BigData systems supporting Datalog or SQL. The problem with these languages and systems is that, to avoid the semantic and computational issues created by non-monotonic constructs in recursion, they only allow programs that are stratified with respect to negation and aggregates. Now, while this crippling restriction is well-justified for negation, it is frequently unjustified for aggregates, since (i) aggregates are often monotonic in the standard lattice of set-containment, (ii) the PreM property guarantees that programs with extrema in recursion are equivalent to stratified programs where extrema are used as post-constraints, and (iii) any program computing any aggregates on sets of facts of predictable cardinality tantamounts to stratified programs where the precomputation of the cardinality of the set is followed by a stratum where recursive rules only use monotonic constructs. With (i) and (ii) covered in previous papers, this paper focuses on (iii) using examples of great practical interest. For such examples, we provide a formal semantics that is conducive to efficient and scalable implementations via well-known techniques such as semi-naive fixpoint currently supported by most Datalog and SQL3 systems.


title: Planted Models for $k$-way Edge and Vertex Expansion

abstract: Graph partitioning problems are a central topic of study in algorithms and complexity theory. Edge expansion and vertex expansion, two popular graph partitioning objectives, seek a $2$-partition of the vertex set of the graph that minimizes the considered objective. However, for many natural applications, one might require a graph to be partitioned into $k$ parts, for some $k \geq 2$. For a $k$-partition $S_1, \ldots, S_k$ of the vertex set of a graph $G = (V,E)$, the $k$-way edge expansion (resp. vertex expansion) of $\{S_1, \ldots, S_k\}$ is defined as $\max_{i \in [k]} \Phi(S_i)$, and the balanced $k$-way edge expansion (resp. vertex expansion) of $G$ is defined as \[ \min_{ \{S_1, \ldots, S_k\} \in \mathcal{P}_k} \max_{i \in [k]} \Phi(S_i) \, , \] where $\mathcal{P}_k$ is the set of all balanced $k$-partitions of $V$ (i.e each part of a $k$-partition in $\mathcal{P}_k$ should have cardinality $|V|/k$), and $\Phi(S)$ denotes the edge expansion (resp. vertex expansion) of $S \subset V$. We study a natural planted model for graphs where the vertex set of a graph has a $k$-partition $S_1, \ldots, S_k$ such that the graph induced on each $S_i$ has large expansion, but each $S_i$ has small edge expansion (resp. vertex expansion) in the graph. We give bi-criteria approximation algorithms for computing the balanced $k$-way edge expansion (resp. vertex expansion) of instances in this planted model.


title: Dissipation-facilitated molecules in a Fermi gas with non-Hermitian spin-orbit coupling

abstract: We study the impact of non-Hermiticity on the molecule formation in a two-component spin-orbit-coupled Fermi gas near a wide Feshbach resonance. Under an experimentally feasible configuration where the two-photon Raman process is dissipative, the Raman-induced synthetic spin-orbit coupling acquires a complex strength. Remarkably, dissipation of the system facilitates the formation and binding of molecules, which, despite their dissipative nature and finite lifetime, exist over a wider parameter regime than in the corresponding Hermitian system. These dissipation-facilitated molecules can be probed by the inverse radio-frequency (rf) spectroscopy, provided the Raman lasers are blue detuned to the excited state. The effects of dissipation manifest in the rf spectra as shifted peaks with broadened widths, which serve as a clear experimental signature. Our results, readily observable in current cold-atom experiments, shed light on the fascinating interplay of non-Hermiticity and interaction in few- and many-body open quantum systems.


title: The attoclock and the tunnelling time debate

abstract: Attosecond angular streaking, also known as the "attoclock", employs a short elliptically polarized laser pulse to tunnel ionize an electron from an atom or a molecule and to put a time stamp on this process by deflecting the photoelectron in the angular spatial direction. This deflection can be used to evaluate the time the tunneling electron spends under the classically inaccessible barrier and to determine whether this time is finite. In this review, we examine the latest experimental and theoretical findings and present a comprehensive set of evidence supporting the zero tunneling time scenario.


title: Bayesian Symbolic Regression

abstract: Interpretability is crucial for machine learning in many scenarios such as quantitative finance, banking, healthcare, etc. Symbolic regression (SR) is a classic interpretable machine learning method by bridging X and Y using mathematical expressions composed of some basic functions. However, the search space of all possible expressions grows exponentially with the length of the expression, making it infeasible for enumeration. Genetic programming (GP) has been traditionally and commonly used in SR to search for the optimal solution, but it suffers from several limitations, e.g. the difficulty in incorporating prior knowledge; overly-complicated output expression and reduced interpretability etc. To address these issues, we propose a new method to fit SR under a Bayesian framework. Firstly, Bayesian model can naturally incorporate prior knowledge (e.g., preference of basis functions, operators and raw features) to improve the efficiency of fitting SR. Secondly, to improve interpretability of expressions in SR, we aim to capture concise but informative signals. To this end, we assume the expected signal has an additive structure, i.e., a linear combination of several concise expressions, whose complexity is controlled by a well-designed prior distribution. In our setup, each expression is characterized by a symbolic tree, and the proposed SR model could be solved by sampling symbolic trees from the posterior distribution using an efficient Markov chain Monte Carlo (MCMC) algorithm. Finally, compared with GP, the proposed BSR(Bayesian Symbolic Regression) method saves computer memory with no need to keep an updated 'genome pool'. Numerical experiments show that, compared with GP, the solutions of BSR are closer to the ground truth and the expressions are more concise. Meanwhile we find the solution of BSR is robust to hyper-parameter specifications such as the number of trees.


title: Compressible Euler Equations on a Sphere and Elliptic-Hyperbolic Property

abstract: In this work we systematically derive the governing equations of supersonic conical flow by projecting the 3D Euler equations onto the unit sphere. These equations result from taking the assumption of conical invariance on the 3D flow field. Under this assumption, the compressible Euler equations reduce to a system defined on the surface of the unit sphere. This compressible flow problem has been successfully used to study steady supersonic flow past cones of arbitrary cross section by reducing the number of spatial dimensions from 3 down to 2 while still capturing many of the relevant 3D effects. In this paper the powerful machinery of tensor calculus is utilized to avoid reference to any particular coordinate system. With the flexibility to use any coordinate system on the surface of a sphere, the equations can be more readily solved numerically when a structured mesh is used by defining the mesh lines to be the coordinate lines. The type of the system of partial differential equations would be hyperbolic or elliptic based on whether the crossflow Mach number is supersonic or subsonic.


title: How to extract the "Abelian" part of double-winding Wilson loop

abstract: It is known that the naive Abelian Wilson loop defined by the Abelian projection cannot reproduce the correct behavior of the double-winding Wilson loop. It is also known that the naive Abelian Wilson loop cannot reproduce the correct behavior of the Wilson loops in higher representations, but this problem was recently solved by using the redefined "Abelian" Wilson loop. In this talk, we will give another reason why this redefined "Abelian" Wilson loop behaves correctly, and by following the same line of the argument, we will propose redefined "Abelian" double-winding Wilson loop which is considered to behave correctly.


title: Further Bijections to Pattern-Avoiding Valid Hook Configurations

abstract: Valid hook configurations are combinatorial objects used to understand West's stack-sorting map. We extend existing bijections corresponding valid hook configurations to intervals in partial orders on Motzkin paths. To enumerate valid hook configurations on $312$-avoiding permutations, we build off of an existing bijection into a Motzkin poset and construct a bijection to certain well-studied closed lattice walks in the first quadrant. We use existing results about these lattice paths to show that valid hook configurations on $312$-avoiding permutations are not counted by a $D$-finite generating function, resolving a question of Defant's, and additionally to compute asymptotics for the number of such configurations. We also extend a bijection of Defant's to a correspondence between valid hook configurations on $132$-avoiding permutations and intervals in the Motzkin-Tamari posets, providing a more elegant proof of Defant's enumeration thereof. To investigate this bijection, we present a number of lemmas about valid hook configurations that are generally applicable and further study the bijections of Defant's.


title: Numerical Solution of Compressible Euler and Magnetohydrodynamic flow past an infinite cone

abstract: A numerical scheme is developed for systems of conservation laws on manifolds which arise in high speed aerodynamics and magneto-aerodynamics. The systems are presented in an arbitrary coordinate system on the manifold and involve source terms which account for the curvature of the domain. In order for a numerical method to accurately capture the behavior of the system it is solving, the equations must be discretized in a way that is not only consistent in value, but also models the appropriate character of the system. Such a discretization is presented in this work which preserves the tensorial transformation relationships involved in formulating equations in a curved space. A numerical method is then developed and applied to the conical Euler and Ideal Magnetohydrodynamic equations. To the author's knowledge, this is the first demonstration of a numerical solver for the conical Ideal MHD equations.


title: Unsupervised High-Resolution Depth Learning From Videos With Dual Networks

abstract: Unsupervised depth learning takes the appearance difference between a target view and a view synthesized from its adjacent frame as supervisory signal. Since the supervisory signal only comes from images themselves, the resolution of training data significantly impacts the performance. High-resolution images contain more fine-grained details and provide more accurate supervisory signal. However, due to the limitation of memory and computation power, the original images are typically down-sampled during training, which suffers heavy loss of details and disparity accuracy. In order to fully explore the information contained in high-resolution data, we propose a simple yet effective dual networks architecture, which can directly take high-resolution images as input and generate high-resolution and high-accuracy depth map efficiently. We also propose a Self-assembled Attention (SA-Attention) module to handle low-texture region. The evaluation on the benchmark KITTI and Make3D datasets demonstrates that our method achieves state-of-the-art results in the monocular depth estimation task.


title: Moving Indoor: Unsupervised Video Depth Learning in Challenging Environments

abstract: Recently unsupervised learning of depth from videos has made remarkable progress and the results are comparable to fully supervised methods in outdoor scenes like KITTI. However, there still exist great challenges when directly applying this technology in indoor environments, e.g., large areas of non-texture regions like white wall, more complex ego-motion of handheld camera, transparent glasses and shiny objects. To overcome these problems, we propose a new optical-flow based training paradigm which reduces the difficulty of unsupervised learning by providing a clearer training target and handles the non-texture regions. Our experimental evaluation demonstrates that the result of our method is comparable to fully supervised methods on the NYU Depth V2 benchmark. To the best of our knowledge, this is the first quantitative result of purely unsupervised learning method reported on indoor datasets.


title: Matrix-Product Codes over Commutative Rings and Constructions Arising from $(\sigma,\delta)$-Codes

abstract: A well-known lower bound (over finite fields and some special finite commutative rings) on the Hamming distance of a matrix-product code (MPC) is shown to remain valid over any commutative ring $R$. A sufficient condition is given, as well, for such a bound to be sharp. It is also shown that an MPC is free when its input codes are all free, in which case a generating matrix is given. If $R$ is finite, a sufficient condition is provided for the dual of an MPC to be an MPC, a generating matrix for such a dual is given, and characterizations of LCD, self-dual, and self-orthogonal MPCs are presented. Finally, results of this paper are used along with previous results of the authors to construct novel MPCs arising from $(\sigma, \delta)$-codes. Some properties of such constructions are also studied.


title: On Self-Orthogonality and Self-Duality of Matrix-Product Codes over Commutative Rings

abstract: Let $R$ be a commutative ring with identity. The paper studies the problem of self-orthogonality and self-duality matrix-product codes (MPCs) over $R$. Some methods as well as special matrices are introduced for the construction of such MPCs. A characterization of such codes (in a special case) is also given. Some concrete examples are presented throughout the paper.


title: Endowing Deep 3D Models with Rotation Invariance Based on Principal Component Analysis

abstract: In this paper, we propose a simple yet effective method to endow deep 3D models with rotation invariance by expressing the coordinates in an intrinsic frame determined by the object shape itself. Key to our approach is to find such an intrinsic frame which should be unique to the identical object shape and consistent across different instances of the same category, e.g. the frame axes of desks should be all roughly along the edges. Interestingly, the principal component analysis exactly provides an effective way to define such a frame, i.e. setting the principal components as the frame axes. As the principal components have direction ambiguity caused by the sign-ambiguity of eigenvector computation, there exist several intrinsic frames for each object. In order to achieve absolute rotation invariance for a deep model, we adopt the coordinates expressed in all intrinsic frames as inputs to obtain multiple output features, which will be further aggregated as a final feature via a self-attention module. Our method is theoretically rotation-invariant and can be flexibly embedded into the current network architectures. Comprehensive experiments demonstrate that our approach can achieve near state-of-the-art performance on rotation-augmented dataset for ModelNet40 classification and outperform other models on SHREC'17 perturbed retrieval task.


title: Privacy- and Utility-Preserving Textual Analysis via Calibrated Multivariate Perturbations

abstract: Accurately learning from user data while providing quantifiable privacy guarantees provides an opportunity to build better ML models while maintaining user trust. This paper presents a formal approach to carrying out privacy preserving text perturbation using the notion of dx-privacy designed to achieve geo-indistinguishability in location data. Our approach applies carefully calibrated noise to vector representation of words in a high dimension space as defined by word embedding models. We present a privacy proof that satisfies dx-privacy where the privacy parameter epsilon provides guarantees with respect to a distance metric defined by the word embedding space. We demonstrate how epsilon can be selected by analyzing plausible deniability statistics backed up by large scale analysis on GloVe and fastText embeddings. We conduct privacy audit experiments against 2 baseline models and utility experiments on 3 datasets to demonstrate the tradeoff between privacy and utility for varying values of epsilon on different task types. Our results demonstrate practical utility (< 2% utility loss for training binary classifiers) while providing better privacy guarantees than baseline models.


title: Benchmark Dataset for Timetable Optimization of Bus Routes in the City of New Delhi

abstract: Public transport is one of the major forms of transportation in the world. This makes it vital to ensure that public transport is efficient. This research presents a novel real-time GPS bus transit data for over 500 routes of buses operating in New Delhi. The data can be used for modeling various timetable optimization tasks as well as in other domains such as traffic management, travel time estimation, etc. The paper also presents an approach to reduce the waiting time of Delhi buses by analyzing the traffic behavior and proposing a timetable. This algorithm serves as a benchmark for the dataset. The algorithm uses a constrained clustering algorithm for classification of trips. It further analyses the data statistically to provide a timetable which is efficient in learning the inter- and intra-month variations.


title: $hv$-Block Cross Validation is not a BIBD: a Note on the Paper by Jeff Racine (2000)

abstract: This note corrects a mistake in the paper "consistent cross-validatory model-selection for dependent data: $hv$-block cross-validation" by Racine (2000). In his paper, he implied that the therein proposed $hv$-block cross-validation is consistent in the sense of Shao (1993). To get this intuition, he relied on the speculation that $hv$-block is a balanced incomplete block design (BIBD). This note demonstrates that this is not the case, and thus the theoretical consistency of $hv$-block remains an open question. In addition, I also provide a Python program counting the number of occurrences of each sample and each pair of samples.

categories: stat.ML cs.LG math.ST q-bio.QM stat.ME stat.TH


title: Global existence in the critical and subcritical cases to the Fisher-KPP model with nonlocal nonlinear reaction

abstract: The Cauchy problem considered in this paper is the following \begin{align} \left\{ \begin{array}{ll} u_t=\Delta u+u^\alpha\left(M_0- \int_{\mathbb{R}^n} u(x,t)dx\right),\quad & x \in \mathbb{R}^n, t>0, u(x,0)=U_0(x)\geq 0,\quad & x \in \mathbb{R}^n. \end{array} \right. \end{align} where $M_0>0, \alpha >1, n \ge 3$. When the coefficient $M_0-\int_{\mathbb{R}^n} u(x,t) dx$ remains positive, \er{nkpp0} is analogous to \begin{align} \left\{ \begin{array}{ll} u_t=\Delta u+u^\alpha,\quad & x \in \mathbb{R}^n, t>0, u(x,0)=U_0(x)\geq 0,\quad & x \in \mathbb{R}^n. \end{array} \right. \end{align} It is well known that when $1<\alpha \le 1+2/n$, the local solution of \er{fujita} blows up in finite time as long as the initial value is nontrivial. The present paper forms a contrast to \er{fujita} and shows the global existence of solutions to \er{nkpp0} for $1<\alpha \le 1+2/n$ by dealing with the mathematical challenge which is from the nonlocal term $\int_{\mathbb{R}^n} u dx$. It's proved that when $1<\alpha <1+2/n$, such a global bound is obtained for any positive $M_0$ and any non-negative initial data. While if $\alpha=1+2/n$, then the global solution does exist for sufficiently small $M_0$ and any non-negative initial data. Furthermore, the large time behavior of the global solution is also discussed for $\alpha=1+2/n$. Besides, this paper establishes the hyper-contractivity of a global solution in $L^\infty(\mathbb{R}^n)$ with $U_0 \in L^1(\mathbb{R}^n)$ for the case $\alpha=1+2/n$.


title: Self-Adaptive Network Pruning

abstract: Deep convolutional neural networks have been proved successful on a wide range of tasks, yet they are still hindered by their large computation cost in many industrial scenarios. In this paper, we propose to reduce such cost for CNNs through a self-adaptive network pruning method (SANP). Our method introduces a general Saliency-and-Pruning Module (SPM) for each convolutional layer, which learns to predict saliency scores and applies pruning for each channel. Given a total computation budget, SANP adaptively determines the pruning strategy with respect to each layer and each sample, such that the average computation cost meets the budget. This design allows SANP to be more efficient in computation, as well as more robust to datasets and backbones. Extensive experiments on 2 datasets and 3 backbones show that SANP surpasses state-of-the-art methods in both classification accuracy and pruning rate.


title: Visually Exploring Software Maintenance Activities

abstract: Lehman's Laws teach us that a software system will become progressively less satisfying to its users over time, unless it is continually adapted to meet new needs. A line of previous works sought to better understand software maintenance by studying how commits can be classified into three main software maintenance activities. Corrective: fault fixing; Perfective: system improvements; Adaptive: new feature introduction. In this work we suggest visualizations for exploring software maintenance activities in both project and individual developer scopes. We demonstrate our approach using a prototype we have built using the Shiny R framework. In addition, we have also published our prototype as an online demo. This demo allows users to explore the maintenance activities of a number of popular open source projects. We believe that the visualizations we provide can assist practitioners in monitoring and maintaining the health of software projects. In particular, they can be useful for identifying general imbalances, peaks, deeps and other anomalies in projects' and developers' maintenance activities.


title: Processing Large Datasets of Fined Grained Source Code Changes

abstract: In the era of Big Code, when researchers seek to study an increasingly large number of repositories to support their findings, the data processing stage may require manipulating millions and more of records. In this work we focus on studies involving fine-grained AST level source code changes. We present how we extended the CodeDistillery source code mining framework with data manipulation capabilities, aimed to alleviate the processing of large datasets of fine grained source code changes. The capabilities we have introduced allow researchers to highly automate their repository mining process and streamline the data acquisition and processing phases. These capabilities have been successfully used to conduct a number of studies, in the course of which dozens of millions of fine-grained source code changes have been processed.


title: Sparse-Dense Subspace Clustering

abstract: Subspace clustering refers to the problem of clustering high-dimensional data into a union of low-dimensional subspaces. Current subspace clustering approaches are usually based on a two-stage framework. In the first stage, an affinity matrix is generated from data. In the second one, spectral clustering is applied on the affinity matrix. However, the affinity matrix produced by two-stage methods cannot fully reveal the similarity between data points from the same subspace (intra-subspace similarity), resulting in inaccurate clustering. Besides, most approaches fail to solve large-scale clustering problems due to poor efficiency. In this paper, we first propose a new scalable sparse method called Iterative Maximum Correlation (IMC) to learn the affinity matrix from data. Then we develop Piecewise Correlation Estimation (PCE) to densify the intra-subspace similarity produced by IMC. Finally we extend our work into a Sparse-Dense Subspace Clustering (SDSC) framework with a dense stage to optimize the affinity matrix for two-stage methods. We show that IMC is efficient when clustering large-scale data, and PCE ensures better performance for IMC. We show the universality of our SDSC framework as well. Experiments on several data sets demonstrate the effectiveness of our approaches. Moreover, we are the first one to apply densification on affinity matrix before spectral clustering, and SDSC constitutes the first attempt to build a universal three-stage subspace clustering framework.


title: Improving Sequence Modeling Ability of Recurrent Neural Networks via Sememes

abstract: Sememes, the minimum semantic units of human languages, have been successfully utilized in various natural language processing applications. However, most existing studies exploit sememes in specific tasks and few efforts are made to utilize sememes more fundamentally. In this paper, we propose to incorporate sememes into recurrent neural networks (RNNs) to improve their sequence modeling ability, which is beneficial to all kinds of downstream tasks. We design three different sememe incorporation methods and employ them in typical RNNs including LSTM, GRU and their bidirectional variants. In evaluation, we use several benchmark datasets involving PTB and WikiText-2 for language modeling, SNLI for natural language inference and another two datasets for sentiment analysis and paraphrase detection. Experimental results show evident and consistent improvement of our sememe-incorporated models compared with vanilla RNNs, which proves the effectiveness of our sememe incorporation methods. Moreover, we find the sememe-incorporated models have higher robustness and outperform adversarial training in defending adversarial attack. All the code and data of this work can be obtained at https://github.com/thunlp/SememeRNN.


title: New Regularity Criteria for the Navier-Stokes Equations in Terms of Pressure

abstract: In this paper, we generalize the main results of [1] and [31] to Lorentz spaces, using a simple procedure. The main results are the following. Let $n\geq 3$ and let $u$ be a Leray-Hopf solution to the $n$-dimensional Navier-Stokes equations with viscosity $\nu$ and divergence free initial condition $u_0\in L^2(\mathbb{R}^n)\cap L^{k}(\mathbb{R}^n)$ (where $k=k(s)$ is sufficiently large). Then there exists a constant $c>0$ such that if \begin{equation} \|p\|_{L^{r,\infty}(0,\infty;L^{s,\infty}(\mathbb{R}^n))}<c\hspace{10mm}\frac{n}{s}+\frac{2}{r}\leq 2,\hspace{5mm}s>\frac{n}{2} \end{equation} or \begin{equation} \|\nabla p\|_{L^{r,\infty}(0,\infty;L^{s,\infty}(\mathbb{R}^n))}<c\hspace{10mm}\frac{n}{s}+\frac{2}{r}\leq 3,\hspace{5mm}s>\frac{n}{3} \end{equation} then $u$ is smooth on $(0, \infty) \times \mathbb{R}^n$. Partial results in the case $n=3$ were obtained in [32], [33] and then recently extended to all appropriate pairs of $r,s$ in [14]. Our results present a unified proof which works for all dimensions $n\geq 3$ and the full range or admissible pairs, $(s,r)$.


title: Magnetic and transport properties in magnetic topological insulators MnBi2Te4(Bi2Te3)n (n=1,2)

abstract: The observation of quantized anomalous Hall conductance in the forced ferromagnetic state of MnBi2Te4 thin flakes has attracted much attentions. However, strong magnetic field is needed to fully polarize the magnetic moments due to the large antiferromagnetic interlayer exchange coupling. Here, we reported the magnetic and electrical transport properties of the magnetic van der Waals MnBi2Te4(Bi2Te3)n (n=1,2) single crystals, in which the interlayer antiferromagnetic exchange coupling is greatly suppressed with the increase of the separation layers Bi2Te3. MnBi4Te7 and MnBi6Te10 show weak antiferromagnetic transition at 12.3 and 10.5 K, respectively. The ferromagnetic hysteresis was observed at low temperature for both of the crystals, which is quite crucial for realizing the quantum anomalous Hall effect without external magnetic field. Our work indicates that MnBi2Te4(Bi2Te3)n (n=1,2) provide ideal platforms to investigate the rich topological phases with going to their 2D limits.


title: Bijections on Dyck tilings: DTS/DTR bijections, Dyck tableaux and tree-like tableaux

abstract: Dyck tilings are certain tilings in the region surrounded by two Dyck paths. We study bijections and combinatorial objects bijective to Dyck tilings, which include Dyck tiling strip (DTS) and Dyck tiling ribbon (DTR) bijections, increasing and decreasing trees, Hermite histories, Dyck tableaux and tree-like tableaux. Dyck tableaux and tree-like tableaux are originally defined for a zigzag path, or equivalently a permutation. We generalize them to the case of general Dyck paths. We show that the most properties of Dyck tableaux can be generalized to the generic case, and show some enumerative results on generalized tree-like tableaux. We also show connections among DTS and DTR bijections, Hermite histories, involutions on increasing and decreasing trees and the reflection of Dyck tilings.