title: Deformations of path algebras of quivers with relations

abstract: Let $A = \Bbbk Q / I$ be the path algebra of any finite quiver $Q$ modulo any two-sided ideal $I$ of relations and let $R$ be any reduction system satisfying the diamond condition for $I$. We introduce an intrinsic notion of deformation of reduction systems and show that there is an equivalence of deformation problems between deformations of the associative algebra $A$ and deformations of the reduction system $R$, the latter being controlled by a natural, explicit L$_\infty$ algebra. It follows in particular that any formal deformation of the associative multiplication on $A$ can, up to gauge equivalence, be given by a combinatorially defined star product, and the approach via reduction systems can be used to give a concrete and complete description of the deformation theory of $A$. For the polynomial algebra in a finite number of variables, this combinatorial star product can be described via bidifferential operators associated to graphs, which we compare to the graphs appearing in Kontsevich's universal quantization formula. Using the notion of admissible orders on the set of paths of the quiver $Q$, we give criteria for the existence of algebraizations of formal deformations, which we also interpret geometrically via algebraic varieties of reduction systems. In this context the Maurer-Cartan equation of the L$_\infty$ algebra can be viewed as a generalization of the Braverman-Gaitsgory criterion for Poincar\'e-Birkhoff-Witt deformations of Koszul algebras.


title: On Thompson Sampling with Langevin Algorithms

abstract: Thompson sampling for multi-armed bandit problems is known to enjoy favorable performance in both theory and practice. However, it suffers from a significant limitation computationally, arising from the need for samples from posterior distributions at every iteration. We propose two Markov Chain Monte Carlo (MCMC) methods tailored to Thompson sampling to address this issue. We construct quickly converging Langevin algorithms to generate approximate samples that have accuracy guarantees, and we leverage novel posterior concentration rates to analyze the regret of the resulting approximate Thompson sampling algorithm. Further, we specify the necessary hyperparameters for the MCMC procedure to guarantee optimal instance-dependent frequentist regret while having low computational complexity. In particular, our algorithms take advantage of both posterior concentration and a sample reuse mechanism to ensure that only a constant number of iterations and a constant amount of data is needed in each round. The resulting approximate Thompson sampling algorithm has logarithmic regret and its computational complexity does not scale with the time horizon of the algorithm.


title: NeurIPS 2019 Disentanglement Challenge: Improved Disentanglement through Aggregated Convolutional Feature Maps

abstract: This report to our stage 1 submission to the NeurIPS 2019 disentanglement challenge presents a simple image preprocessing method for training VAEs leading to improved disentanglement compared to directly using the images. In particular, we propose to use regionally aggregated feature maps extracted from CNNs pretrained on ImageNet. Our method achieved the 2nd place in stage 1 of the challenge. Code is available at https://github.com/mseitzer/neurips2019-disentanglement-challenge.


title: Nanoscale electronic inhomogeneity in FeSe$_{0.4}$Te$_{0.6}$ revealed through unsupervised machine learning

abstract: We report on an apparent low-energy nanoscale electronic inhomogeneity in FeSe$_{0.4}$Te$_{0.6}$ due to the distribution of selenium and tellurium atoms revealed through unsupervised machine learning. Through an unsupervised clustering algorithm, characteristic spectra of selenium- and tellurium-rich regions are identified. The inhomogeneity linked to these spectra can clearly be traced in the differential conductance and is detected both at energy scales of a few electron volts as well as within a few millielectronvolts of the Fermi energy. By comparison with ARPES, this inhomogeneity can be linked to an electron-like band just above the Fermi energy. It is directly correlated with the local distribution of selenium and tellurium. There is no clear correlation with the magnitude of the superconducting gap, however the height of the coherence peaks shows significant correlation with the intensity with which this band is detected, and hence with the local chemical composition.


title: Surprising Variability of the Planetary Nebula IC 4997 = QV Sge

abstract: We present the results of a new epoch of a long-term photometric monitoring of the variable planetary nebula IC~4997. The integral $UBV$ light curves display a continuing brightening of 0.$^{m}15$ in $V$, a slight rise (<0.$^{m}1$) in $B$, and constancy in $U$. The $B-V$ color has got redder from 0.$^{m}4$ in 2000 to 0.$^{m}7$ in 2019, whereas the $U-B$ color has not changed significantly. New $JHKL$ photometry was obtained in 2019. We found the source to be fainter by 0.$^{m}4$ in $L$ and bluer by 0.$^{m}2$ in the $K-L$ color if compared to the data of 1999-2006. The long-term brightness variations are due mostly to the changing input of emission lines to the integral light. Low-dispersion spectroscopic observations carried out in 2010-2019 revealed a continuing decrease in the $\lambda4363/\lambda4340$ ratio: it decreased by a factor of $\sim3$ in 30 years and reached the level of 1960-1970. The absolute intensities of [OIII] nebualr lines increased by a factor of $>2$ from 1990 to 2019, whereas the [OIII] $\lambda$4363 line had weakened by a factor of 2 comparing to the maximum value observed in 2000. The variation of H$\beta$ absolute intensity in 1960-2019 was shown to be similar to that of [OIII] nebular lines, but of smaller amplitude. The electron density in the outer part of the nebula was estimated from the [SII] and [ClIII] lines. Basing on the data on absolute intensities we propose a possible scenario describing the change of $N_e,T_e$ in 1970-2019. The spectral variability of IC~4997 could be explained by a variation of electron temperature in the nebula caused by not so much the change in ionizing flux from the central star as the variable stellar wind and related processes. The photometric and spectral changes observed for IC~4997 in 1960-2019 may be interpreted as an observable consequence of a single episode of enhanced mass loss from the variable central star.


title: On the Modularity of Hypernetworks

abstract: In the context of learning to map an input $I$ to a function $h_I:\mathcal{X}\to \mathbb{R}$, two alternative methods are compared: (i) an embedding-based method, which learns a fixed function in which $I$ is encoded as a conditioning signal $e(I)$ and the learned function takes the form $h_I(x) = q(x,e(I))$, and (ii) hypernetworks, in which the weights $\theta_I$ of the function $h_I(x) = g(x;\theta_I)$ are given by a hypernetwork $f$ as $\theta_I=f(I)$. In this paper, we define the property of modularity as the ability to effectively learn a different function for each input instance $I$. For this purpose, we adopt an expressivity perspective of this property and extend the theory of Devore et al. 1996 and provide a lower bound on the complexity (number of trainable parameters) of neural networks as function approximators, by eliminating the requirements for the approximation method to be robust. Our results are then used to compare the complexities of $q$ and $g$, showing that under certain conditions and when letting the functions $e$ and $f$ be as large as we wish, $g$ can be smaller than $q$ by orders of magnitude. This sheds light on the modularity of hypernetworks in comparison with the embedding-based method. Besides, we show that for a structured target function, the overall number of trainable parameters in a hypernetwork is smaller by orders of magnitude than the number of trainable parameters of a standard neural network and an embedding method.


title: A Critical View of the Structural Causal Model

abstract: In the univariate case, we show that by comparing the individual complexities of univariate cause and effect, one can identify the cause and the effect, without considering their interaction at all. In our framework, complexities are captured by the reconstruction error of an autoencoder that operates on the quantiles of the distribution. Comparing the reconstruction errors of the two autoencoders, one for each variable, is shown to perform surprisingly well on the accepted causality directionality benchmarks. Hence, the decision as to which of the two is the cause and which is the effect may not be based on causality but on complexity. In the multivariate case, where one can ensure that the complexities of the cause and effect are balanced, we propose a new adversarial training method that mimics the disentangled structure of the causal model. We prove that in the multidimensional case, such modeling is likely to fit the data only in the direction of causality. Furthermore, a uniqueness result shows that the learned model is able to identify the underlying causal and residual (noise) components. Our multidimensional method outperforms the literature methods on both synthetic and real world datasets.


title: Conditional regression for single-index models

abstract: The single-index model is a statistical model for intrinsic regression where responses are assumed to depend on a single yet unknown linear combination of the predictors, allowing to express the regression function as $ \mathbb{E} [ Y | X ] = f ( \langle v , X \rangle ) $ for some unknown \emph{index} vector $v$ and \emph{link} function $f$. Conditional methods provide a simple and effective approach to estimate $v$ by averaging moments of $X$ conditioned on $Y$, but depend on parameters whose optimal choice is unknown and do not provide generalization bounds on $f$. In this paper we propose a new conditional method converging at $\sqrt{n}$ rate under an explicit parameter characterization. Moreover, we prove that polynomial partitioning estimates achieve the $1$-dimensional min-max rate for regression of H\"older functions when combined to any $\sqrt{n}$-convergent index estimator. Overall this yields an estimator for dimension reduction and regression of single-index models that attains statistical optimality in quasilinear time.


title: Fighting Fire with Light: A Case for Defending DDoS Attacks Using the Optical Layer

abstract: The DDoS attack landscape is growing at an unprecedented pace. Inspired by the recent advances in optical networking, we make a case for optical layer-aware DDoS defense (O-LAD) in this paper. Our approach leverages the optical layer to isolate attack traffic rapidly via dynamic reconfiguration of (backup) wavelengths using ROADMs---bridging the gap between (a) evolution of the DDoS attack landscape and (b) innovations in the optical layer (e.g., reconfigurable optics). We show that the physical separation of traffic profiles allows finer-grained handling of suspicious flows and offers better performance for benign traffic in the face of an attack. We present preliminary results modeling throughput and latency for legitimate flows while scaling the strength of attacks. We also identify a number of open problems for the security, optical, and systems communities: modeling diverse DDoS attacks (e.g., fixed vs. variable rate, detectable vs. undetectable), building a full-fledged defense system with optical advancements (e.g., OpenConfig), and optical layer-aware defenses for a broader class of attacks (e.g., network reconnaissance).


title: Driving with Data in the Motor City: Mining and Modeling Vehicle Fleet Maintenance Data

abstract: The City of Detroit maintains an active fleet of over 2500 vehicles, spending an annual average of over \$5 million on purchases and over \$7.7 million on maintenance. Modeling patterns and trends in this data is of particular importance to a variety of stakeholders, particularly as Detroit emerges from Chapter 9 bankruptcy, but the structure in such data is complex, and the city lacks dedicated resources for in-depth analysis. The City of Detroit's Operations and Infrastructure Group and the University of Michigan initiated a collaboration which seeks to address this unmet need by analyzing data from the City of Detroit's vehicle fleet. This work presents a case study and provides the first data-driven benchmark, demonstrating a suite of methods to aid in data understanding and prediction for large vehicle maintenance datasets. We present analyses to address three key questions raised by the stakeholders, related to discovering multivariate maintenance patterns over time; predicting maintenance; and predicting vehicle- and fleet-level costs. We present a novel algorithm, PRISM, for automating multivariate sequential data analyses using tensor decomposition. This work is a first of its kind that presents both methodologies and insights to guide future civic data research.


title: Geometric Algebra Power Theory (GAPoT): Revisiting Apparent Power under Non-Sinusoidal Conditions

abstract: Traditional power theories and one of their most important concepts --apparent power-- are still a source of debate and, as shown in the literature, they present several flaws that misinterpret the power-transfer phenomena under distorted grid conditions. In recent years, advanced mathematical tools such as geometric algebra (GA) have been applied to address these issues. However, the application of GA to electrical circuits requires more consensus, improvements and refinement. In this paper, power theories based on GA are revisited. Several drawbacks and inconsistencies of previous works are identified and modifications to the so-called geometric algebra power theory (GAPoT) are presented. This theory takes into account power components generated by cross-products between current and voltage harmonics in the frequency domain. Compared to other theories based on GA, it is compatible with the traditional definition of apparent power calculated as the product of RMS voltage and current. Also, mathematical developments are done in a multi-dimensional Euclidean space where the energy conservation principle is satisfied. The paper includes a basic example and experimental results in which measurements from a utility supply are analysed. Finally, suggestions for the extension to three-phase systems are drawn.


title: The square root of the vacuum I. Equivariance for the $\kappa$-symmetry superdistribution

abstract: A complete and natural geometric and physical interpretation of the tangential gauge supersymmetry, also known as $\kappa$-symmetry, of a large class of Green-Schwarz(-type) super-$\sigma$-models for the super-$p$-brane in a homogeneous space of a (supersymmetry) Lie supergroup is established in the convenient setting of the topological Hughes-Polchinski formulation of the super-$\sigma$-model and illustrated on a number of physical examples. The supersymmetry is identified as an odd superdistribution in the tangent sheaf of the supertarget of the super-$\sigma$-model, generating - through its weak derived flag - the vacuum foliation of the supertarget. It is also demonstrated to canonically lift to the vacuum restriction of the extended Hughes-Polchinski $p$-gerbe associated with the superbackground of the field theory, and that in the form of a canonical linearised equivariant structure thereon, canonically compatible with the residual global supersymmetry of the vacuum.


title: Persistent homotopy theory

abstract: Vietoris-Rips and degree Rips complexes are represented as homotopy types by their underlying posets of simplices, and basic homotopy stability theorems are recast in these terms. These homotopy types are viewed as systems (or functors), which are defined on a parameter space. The category of systems of spaces admits a partial homotopy theory that is based on controlled equivalences, suitably defined, that are the output of homotopy stability results.


title: Transversals and bipancyclicity in bipartite graph families

abstract: A bipartite graph is called bipancyclic if it contains cycles of every even length from four up to the number of vertices in the graph. A theorem of Schmeichel and Mitchem states that for $n \geq 4$, every balanced bipartite graph on $2n$ vertices in which each vertex in one color class has degree greater than $\frac{n}{2}$ and each vertex in the other color class has degree at least $\frac{n}{2}$ is bipancyclic. We prove a generalization of this theorem in the setting of graph transversals. Namely, we show that given a family $\mathcal{G}$ of $2n$ bipartite graphs on a common set $X$ of $2n$ vertices with a common balanced bipartition, if each graph of $\mathcal G$ has minimum degree greater than $\frac{n}{2}$ in one color class and minimum degree at least $\frac{n}{2}$ in the other color class, then there exists a cycle on $X$ of each even length $4 \leq \ell \leq 2n$ that uses at most one edge from each graph of $\mathcal G$. We also show that given a family $\mathcal G$ of $n$ bipartite graphs on a common set $X$ of $2n$ vertices meeting the same degree conditions, there exists a perfect matching on $X$ that uses exactly one edge from each graph of $\mathcal G$.


title: One-loop order effects from one universal extra dimension on $\lambda\phi^{4}$ theory

abstract: The self-interacting $\lambda\phi^{4}$ scalar field theory is a warhorse in quantum field theory. Here we explore the one-loop order impact from one universal extra dimension, $S^{1}/\mathbb{Z}_{2}$, to the self-energy and four point vertex functions associated to this theory. Such effects come as an infinite number of UV divergences corresponding to an infinite superposition of excited KK particles around the loop. We show that dimensional regularisation is adequate enough to control them in terms of the product of the one dimensional inhomogenous Epstein zeta function times the gamma function. From the analytical properties of these functions, the UV divergences are extracted and the counterterms defined; the latter turn out to be of canonical dimension four at the Lagrangian level. We use both, the MS-scheme and a mass-dependent subtraction scheme to remove divergences. Only the latter manifestly satisfy the decoupling theorem.


title: Deep Multimodal Image-Text Embeddings for Automatic Cross-Media Retrieval

abstract: This paper considers the task of matching images and sentences by learning a visual-textual embedding space for cross-modal retrieval. Finding such a space is a challenging task since the features and representations of text and image are not comparable. In this work, we introduce an end-to-end deep multimodal convolutional-recurrent network for learning both vision and language representations simultaneously to infer image-text similarity. The model learns which pairs are a match (positive) and which ones are a mismatch (negative) using a hinge-based triplet ranking. To learn about the joint representations, we leverage our newly extracted collection of tweets from Twitter. The main characteristic of our dataset is that the images and tweets are not standardized the same as the benchmarks. Furthermore, there can be a higher semantic correlation between the pictures and tweets contrary to benchmarks in which the descriptions are well-organized. Experimental results on MS-COCO benchmark dataset show that our model outperforms certain methods presented previously and has competitive performance compared to the state-of-the-art. The code and dataset have been made available publicly.


title: EDEN: Sensitivity Analysis and Transiting Planet Detection Limits for Nearby Late Red Dwarfs

abstract: Small planets are common around late-M dwarfs and can be detected through highly precise photometry by the transit method. Planets orbiting nearby stars are particularly important as they are often the best-suited for future follow-up studies. We present observations of three nearby M-dwarfs referred to as EIC-1, EIC-2, and EIC-3, and use them to search for transits and set limits on the presence of planets. On most nights our observations are sensitive to Earth-sized transiting planets, and photometric precision is similar to or better than TESS for faint late-M dwarfs of the same magnitude (I=15 mag). We present our photometry and transit search pipeline, which utilizes simple median detrending in combination with transit least squares based transit detection (Hippke & Heller 2019).For these targets, and transiting planets between one and two Earth radii, we achieve an average transit detection probability of 60% between periods of 0.5 and 2 days, 30% between 2 and 5 days,and 10% between 5 and 10 days. These sensitivities are conservative compared to visual searches.


title: Distributed Quantum Proofs for Replicated Data

abstract: The paper tackles the issue of $\textit{checking}$ that all copies of a large data set replicated at several nodes of a network are identical. The fact that the replicas may be located at distant nodes prevents the system from verifying their equality locally, i.e., by having each node consult only nodes in its vicinity. On the other hand, it remains possible to assign $\textit{certificates}$ to the nodes, so that verifying the consistency of the replicas can be achieved locally. However, we show that, as the data set is large, classical certification mechanisms, including distributed Merlin-Arthur protocols, cannot guarantee good completeness and soundness simultaneously, unless they use very large certificates. The main result of this paper is a distributed $\textit{quantum}$ Merlin-Arthur protocol enabling the nodes to collectively check the consistency of the replicas, based on small certificates, and in a single round of message exchange between neighbors, with short messages. In particular, the certificate-size is logarithmic in the size of the data set, which gives an exponential advantage over classical certification mechanisms.


title: Averaging in the case of multiple invariant measures for the fast system

abstract: We consider the averaging principle for deterministic or stochastic systems with a fast stochastic component (family of continuous-time Markov chains depending on the state of the system as a parameter). We show that, due to bifurcations in the simplex of invariant probability measures of the chains, the limiting system should be considered on a graph or on an open book with certain gluing conditions in the vertices of the graph (or on the bifurcation surface).


title: Optimal Jammer Placement in UAV-assisted Relay Networks

abstract: We consider the relaying application of unmanned aerial vehicles (UAVs), in which UAVs are placed between two transceivers (TRs) to increase the throughput of the system. Instead of studying the placement of UAVs as pursued in existing literature, we focus on investigating the placement of a jammer or a major source of interference on the ground to effectively degrade the performance of the system, which is measured by the maximum achievable data rate of transmission between the TRs. We demonstrate that the optimal placement of the jammer is in general a non-convex optimization problem, for which obtaining the solution directly is intractable. Afterward, using the inherent characteristics of the signal-to-interference ratio (SIR) expressions, we propose a tractable approach to find the optimal position of the jammer. Based on the proposed approach, we investigate the optimal positioning of the jammer in both dual-hop and multi-hop UAV relaying settings. Numerical simulations are provided to evaluate the performance of our proposed method.


title: How Transferable are the Representations Learned by Deep Q Agents?

abstract: In this paper, we consider the source of Deep Reinforcement Learning (DRL)'s sample complexity, asking how much derives from the requirement of learning useful representations of environment states and how much is due to the sample complexity of learning a policy. While for DRL agents, the distinction between representation and policy may not be clear, we seek new insight through a set of transfer learning experiments. In each experiment, we retain some fraction of layers trained on either the same game or a related game, comparing the benefits of transfer learning to learning a policy from scratch. Interestingly, we find that benefits due to transfer are highly variable in general and non-symmetric across pairs of tasks. Our experiments suggest that perhaps transfer from simpler environments can boost performance on more complex downstream tasks and that the requirements of learning a useful representation can range from negligible to the majority of the sample complexity, based on the environment. Furthermore, we find that fine-tuning generally outperforms training with the transferred layers frozen, confirming an insight first noted in the classification setting.


title: Application of ERA5 and MENA simulations to predict offshore wind energy potential

abstract: This study explores wind energy resources in different locations through the Gulf of Oman and also their future variability due climate change impacts. In this regard, EC-EARTH near surface wind outputs obtained from CORDEX-MENA simulations are used for historical and future projection of the energy. The ERA5 wind data are employed to assess suitability of the climate model. Moreover, the ERA5 wave data over the study area are applied to compute sea surface roughness as an important variable for converting near surface wind speeds to those of wind speed at turbine hub-height. Considering the power distribution, bathymetry and distance from the coats, some spots as tentative energy hotspots to provide detailed assessment of directional and temporal variability and also to investigate climate change impact studies. RCP8.5 as a common climatic scenario is used to project and extract future variation of the energy in the selected sites. The results of this study demonstrate that the selected locations have a suitable potential for wind power turbine plan and constructions.


title: Suboptimal Control of Unknown Second-order Nonlinear Systems with Guaranteed Global Convergence

abstract: A suboptimal active disturbance rejection controller (S-ADRC) is proposed for second-order systems with unknown time-varying nonlinear dynamics. The output-feedback controller guarantees a global convergence to the vicinity of an optimal solution by means of dynamic control gains, based on the estimated main and extended state variables obtained through a high-gain observer. Three numerical examples compare the performance of the proposed control scheme applied to linear and nonlinear systems with that of a fixed-gain conventional ADRC as well as several model-based optimal and suboptimal controllers.


title: Spacecraft Tracking Applications of the Square Kilometre Array

abstract: The Square Kilometre Array (SKA) is the next generation radio telescope distinguished by a superb sensitivity due to its large aperture (about one square kilometre) and advanced instrumentation. It will cover a broad range of observing bands including those used for tracking of and communications to deep space missions. While spacecraft tracking is not a main application defining the technical specifications of the SKA, this facility might play a role in tracking deep space probes as a backup to the ``dedicated'' deep space tracking networks. This paper presents possible applications of the SKA as a deep space tracking facility and major related technical specifications of various concepts of the SKA. It was presented at the 3rd International Workshop on Tracking, Telemetry and Command Systems for Space Applications, ESA-ESOC, Darmstadt, Germany, 7-9 September 2004. Over the past years, the SKA concept has developed to a much higher level of detalisation and is currently at the implementation phase. A number of specific considerations in this presentation no longer correspond to the actual status of the SKA project. However, the overall concept of the SKA applications for communication and tracking of interplanetary spacecraft remain topical, and some approaches presented here remain of interest for prospective deep space missions.


title: Triple Wins: Boosting Accuracy, Robustness and Efficiency Together by Enabling Input-Adaptive Inference

abstract: Deep networks were recently suggested to face the odds between accuracy (on clean natural images) and robustness (on adversarially perturbed images) (Tsipras et al., 2019). Such a dilemma is shown to be rooted in the inherently higher sample complexity (Schmidt et al., 2018) and/or model capacity (Nakkiran, 2019), for learning a high-accuracy and robust classifier. In view of that, give a classification task, growing the model capacity appears to help draw a win-win between accuracy and robustness, yet at the expense of model size and latency, therefore posing challenges for resource-constrained applications. Is it possible to co-design model accuracy, robustness and efficiency to achieve their triple wins? This paper studies multi-exit networks associated with input-adaptive efficient inference, showing their strong promise in achieving a "sweet point" in cooptimizing model accuracy, robustness and efficiency. Our proposed solution, dubbed Robust Dynamic Inference Networks (RDI-Nets), allows for each input (either clean or adversarial) to adaptively choose one of the multiple output layers (early branches or the final one) to output its prediction. That multi-loss adaptivity adds new variations and flexibility to adversarial attacks and defenses, on which we present a systematical investigation. We show experimentally that by equipping existing backbones with such robust adaptive inference, the resulting RDI-Nets can achieve better accuracy and robustness, yet with over 30% computational savings, compared to the defended original models.


title: Measure bound for translation surfaces with short saddle connections

abstract: We prove that any ergodic $SL_2(R)$-invariant probability measure on a stratum of translation surfaces satisfies strong regularity: the measure of the set of surfaces with two non-parallel saddle connections of length at most $\epsilon_1, \epsilon_2$ is $O(\epsilon_1^2 \epsilon_2^2)$. We prove a more general theorem which works for any number of short saddle connections. The proof uses the multi-scale compactification of strata recently introduced by Bainbridge-Chen-Gendron-Grushevsky-M\"oller and the algebraicity result of Filip.


title: Renormalization of the band gap in 2D materials through the competition between electromagnetic and four-fermion interactions

abstract: Recently the renormalization of the band gap $m$, in both WSe$_2$ and MoS$_2$, has been experimentally measured as a function of the carrier concentration $n$. The main result establishes a decreasing of hundreds of meV, in comparison with the bare band gap, as the carrier concentration increases. These materials are known as transition metal dichalcogenides and their low-energy excitations are, approximately, described by the massive Dirac equation. Using Pseudo Quantum Electrodynamics (PQED) to describe the electromagnetic interaction between these quasiparticles and from renormalization group analysis, we obtain that the renormalized mass describes the band gap renormalization with a function given by $m(n)/m_0=(n/n_0)^{C_\lambda/2}$, where $m_0=m(n_0)$ and $C_\lambda$ is a function of the coupling constant $\lambda$. We compare our theoretical results with the experimental findings for WSe$_2$ and MoS$_2$, and we conclude that our approach is in agreement with these experimental results for reasonable values of $\lambda$. In addition we introduced a Gross-Neveu (GN) interaction which could simulate an disorder/impurity-like microscopic interaction. In this case, we show that there exists a critical coupling constant, namely, $\lambda_c \approx 0,66$ in which the beta function of the mass vanishes, providing a stable fixed point in the ultraviolet limit. For $\lambda>\lambda_c$, the renormalized mass decreases while for $\lambda<\lambda_c$ it increases with the carrier concentration.


title: Lattice Gauge Theory for a Quantum Computer

abstract: The quantum link~\cite{Brower:1997ha} Hamiltonian was introduced two decades ago as an alternative to Wilson's Euclidean lattice QCD with gauge fields represented by bi-linear fermion/anti-fermion operators. When generalized this new microscopic representation of lattice field theories is referred as {\tt D-theory}~\cite{Brower:2003vy}. Recast as a Hamiltonian in Minkowski space for real time evolution, D-theory leads naturally to quantum Qubit algorithms. Here to explore digital quantum computing for gauge theories, the simplest example of U(1) compact QED on triangular lattice is defined and gauge invariant kernels for the Suzuki-Trotter expansions are expressed as Qubit circuits capable of being tested on the IBM-Q and other existing Noisy Intermediate Scale Quantum (NISQ) hardware. This is a modest step in exploring the quantum complexity of D-theory to guide future applications to high energy physics and condensed matter quantum field theories.


title: Symbolic Querying of Vector Spaces: Probabilistic Databases Meets Relational Embeddings

abstract: We propose unifying techniques from probabilistic databases and relational embedding models with the goal of performing complex queries on incomplete and uncertain data. We formalize a probabilistic database model with respect to which all queries are done. This allows us to leverage the rich literature of theory and algorithms from probabilistic databases for solving problems. While this formalization can be used with any relational embedding model, the lack of a well-defined joint probability distribution causes simple query problems to become provably hard. With this in mind, we introduce \TO, a relational embedding model designed to be a tractable probabilistic database, by exploiting typical embedding assumptions within the probabilistic framework. Using a principled, efficient inference algorithm that can be derived from its definition, we empirically demonstrate that \TOs is an effective and general model for these querying tasks.


title: New Extremal binary self-dual codes of length 68 from a novel approach to neighbors

abstract: In this work, we introduce the concept of distance between self-dual codes, which generalizes the concept of a neighbor for self-dual codes. Using the k-neighbors, we are able to construct extremal binary self-dual codes of length 68 with new weight enumerators. We construct 143 extremal binary self-dual codes of length 68 with new weight enumerators including 42 codes with gamma=8 in their W_{68,2} and 40 with gamma=9 in their W_{68,2}. These examples are the first in the literature for these gamma values. This completes the theoretical list of possible values for gamma in W_{68,2}.


title: On Incompressible Vibrations of the Stratified Atmosphere on the Flat Earth

abstract: We consider the vibrations and waves in the atmosphere under the gravitation on the flat earth. The stratified density distribution of the back ground equilibrium is supposed to touch the vacuum at the finite height of the stratosphere. We show that incompressible motions are possible to create vibrations or progressive waves with countably many time frequencies of the vibrations or speeds of the wave propagation which accumulate to 0, say, of countably infinitely many slow and slow modes.


title: Generalized Octave Convolutions for Learned Multi-Frequency Image Compression

abstract: Learned image compression has recently shown the potential to outperform the standard codecs. State-of-the-art rate-distortion (R-D) performance has been achieved by context-adaptive entropy coding approaches in which hyperprior and autoregressive models are jointly utilized to effectively capture the spatial dependencies in the latent representations. However, the latents are feature maps of the same spatial resolution in previous works, which contain some redundancies that affect the R-D performance. In this paper, we propose the first learned multi-frequency image compression and entropy coding approach that is based on the recently developed octave convolutions to factorize the latents into high and low frequency (resolution) components, where the low frequency is represented by a lower resolution. Therefore, its spatial redundancy is reduced, which improves the R-D performance. Novel generalized octave convolution and octave transposed-convolution architectures with internal activation layers are also proposed to preserve more spatial structure of the information. Experimental results show that the proposed scheme not only outperforms all existing learned methods as well as standard codecs such as the next-generation video coding standard VVC (4:2:0) on the Kodak dataset in both PSNR and MS-SSIM. We also show that the proposed generalized octave convolution can improve the performance of other auto-encoder-based computer vision tasks such as semantic segmentation and image denoising.


title: Status drives how we cite: Evidence from thousands of authors

abstract: Researchers cite works for a variety of reasons, including some having nothing to do with acknowledging influence. The distribution of different citation types in the literature, and which papers attract which types, is poorly understood. We investigate high-influence and low-influence citations and the mechanisms producing them using 17,154 ground-truth citation types provided via survey by 9,380 authors systematically sampled across academic fields. Overall, 54% of citations denote little-to-no influence and these citations are concentrated among low status (lightly cited) papers. In contrast, high-influence citations are concentrated among high status (highly cited) papers through a number of steps that resemble a pipeline. Authors discover highly cited papers earlier in their projects, more often through social contacts, and read them more closely. Papers' status, above and beyond any quality differences, directly helps determine their pipeline: experimentally revealing or hiding citation counts during the survey shows that low counts cause lowered perceptions of quality. Accounting for citation types thus reveals a "double status effect": in addition to affecting how often a work is cited, status affects how meaningfully it is cited. Consequently, highly cited papers are even more influential than their raw citation counts suggest.


title: Predicting Rate of Cognitive Decline at Baseline Using a Deep Neural Network with Multidata Analysis

abstract: Purpose: This study investigates whether a machine-learning-based system can predict the rate of cognitive decline in mildly cognitively impaired patients by processing only the clinical and imaging data collected at the initial visit. Approach: We built a predictive model based on a supervised hybrid neural network utilizing a 3-Dimensional Convolutional Neural Network to perform volume analysis of Magnetic Resonance Imaging and integration of non-imaging clinical data at the fully connected layer of the architecture. The experiments are conducted on the Alzheimers Disease Neuroimaging Initiative dataset. Results: Experimental results confirm that there is a correlation between cognitive decline and the data obtained at the first visit. The system achieved an area under the receiver operator curve (AUC) of 0.70 for cognitive decline class prediction. Conclusion: To our knowledge, this is the first study that predicts slowly deteriorating/stable or rapidly deteriorating classes by processing routinely collected baseline clinical and demographic data (Baseline MRI, Baseline MMSE, Scalar Volumetric data, Age, Gender, Education, Ethnicity, and Race). The training data is built based on MMSE-rate values. Unlike the studies in the literature that focus on predicting Mild Cognitive Impairment-to-Alzheimer`s disease conversion and disease classification, we approach the problem as an early prediction of cognitive decline rate in MCI patients.


title: Improving the Linkage Construction with Echelon-Ferrers for Constant-Dimension Codes

abstract: Echelon-Ferrers is an important method to improve lower bounds for constant-dimension codes, which can be applied on various parameters. Fagang Li [12] combined the linkage construction and echelon-Ferrers to obtain some new lower bounds of constant-dimension codes. In this letter, we generalize this linkage construction to obtain new lower bounds.


title: Phase diagram of the interacting persistent spin-helix state

abstract: We study the phase diagram of the interacting two-dimensional electron gas (2DEG) with equal Rashba and Dresselhaus spin-orbit coupling, which for weak coupling gives rise to the well-known persistent spin-helix phase. We construct the full Hartree-Fock phase diagram using a classical Monte-Carlo method analogous to that used in Phys.Rev.B 96, 235425 (2017). For the 2DEG with only Rashba spin-orbit coupling, it was found that at intermediate values of the Wigner-Seitz radius rs the system is characterized by a single Fermi surface with an out-of-plane spin polarization, while at slightly larger values of rs it undergoes a transition to a state with a shifted Fermi surface and an in-plane spin polarization. The various phase transitions are first-order, and this shows up in discontinuities in the conductivity and the appearance of anisotropic resistance in the in-plane polarized phase. In this work, we show that the out-of-plane spin-polarized region shrinks as the strength of the Dresselhaus spin-orbit interaction increases, and entirely vanishes when the Rashba and Dresselhaus spin-orbit coupling strengths are equal. At this point, the system can be mapped onto a 2DEG without spin-orbit coupling, and this transformation reveals the existence of an in-plane spin-polarized phase with a single, displaced Fermi surface beyond rs > 2.01. This is confirmed by classical Monte-Carlo simulations. We discuss experimental observation and useful applications of the novel phase, as well as caveats of using the classical Monte-Carlo method.


title: Detecting electron-phonon couplings during photo-induced phase transition

abstract: Photo-induced phase transitions have been intensively studied owing to the ability to control a material of interest in the ultrafast manner, which can induce exotic phases unable to be attained at equilibrium. However, the key mechanisms are still under debate, and it has currently been a central issue how the couplings between the electron, lattice, and spin degrees of freedom are evolving during photo-induced phase transitions. Here, we develop a new analysis method, frequency-domain angle-resolved photoemission spectroscopy, to gain precise insight into electron-phonon couplings during photo-induced insulator-to-metal transitions for Ta$_2$NiSe$_5$. We demonstrate that multiple coherent phonons generated by displacive excitations show band-selective coupling to the electrons. Furthermore, we find that the lattice modulation corresponding to the 2 THz phonon mode, where Ta lattice is sheared along the a-axis, is the most relevant for the photo-induced semimetallic state.


title: Bayesian bandwidth estimation and semi-metric selection for a functional partial linear model with unknown error density

abstract: This study examines the optimal selections of bandwidth and semi-metric for a functional partial linear model. Our proposed method begins by estimating the unknown error density using a kernel density estimator of residuals, where the regression function, consisting of parametric and nonparametric components, can be estimated by functional principal component and functional Nadayara-Watson estimators. The estimation accuracy of the regression function and error density crucially depends on the optimal estimations of bandwidth and semi-metric. A Bayesian method is utilized to simultaneously estimate the bandwidths in the regression function and kernel error density by minimizing the Kullback-Leibler divergence. For estimating the regression function and error density, a series of simulation studies demonstrate that the functional partial linear model gives improved estimation and forecast accuracies compared with the functional principal component regression and functional nonparametric regression. Using a spectroscopy dataset, the functional partial linear model yields better forecast accuracy than some commonly used functional regression models. As a by-product of the Bayesian method, a pointwise prediction interval can be obtained, and marginal likelihood can be used to select the optimal semi-metric.


title: Computing Bi-Lipschitz Outlier Embeddings into the Line

abstract: The problem of computing a bi-Lipschitz embedding of a graphical metric into the line with minimum distortion has received a lot of attention. The best-known approximation algorithm computes an embedding with distortion $O(c^2)$, where $c$ denotes the optimal distortion [B\u{a}doiu \etal~2005]. We present a bi-criteria approximation algorithm that extends the above results to the setting of \emph{outliers}. Specifically, we say that a metric space $(X,\rho)$ admits a $(k,c)$-embedding if there exists $K\subset X$, with $|K|=k$, such that $(X\setminus K, \rho)$ admits an embedding into the line with distortion at most $c$. Given $k\geq 0$, and a metric space that admits a $(k,c)$-embedding, for some $c\geq 1$, our algorithm computes a $({\mathsf p}{\mathsf o}{\mathsf l}{\mathsf y}(k, c, \log n), {\mathsf p}{\mathsf o}{\mathsf l}{\mathsf y}(c))$-embedding in polynomial time. This is the first algorithmic result for outlier bi-Lipschitz embeddings. Prior to our work, comparable outlier embeddings where known only for the case of additive distortion.


title: Links in Surfaces and Laplacian Modules

abstract: Laplacian matrices of weighted graphs in surfaces $S$ are used to define module and polynomial invariants of $Z/2$-homologically trivial links in $S \times [0,1]$. Information about virtual genus is obtained.


title: Photothermally Induced Transparency

abstract: Induced transparency is a common but remarkable effect in optics. It occurs when a strong driving field is used to render an otherwise opaque material transparent. The effect is known as electromagnetically induced transparency in atomic media and optomechanically induced transparency in systems that consist of coupled optical and mechanical resonators. In this work, we introduce the concept of photothermally induced transparency (PTIT). It happens when an optical resonator exhibits non-linear behavior due to optical heating of the resonator or its mirrors. Similar to the established mechanisms for induced transparency, PTIT can suppress the coupling between an optical resonator and a traveling optical field. We further show that the dispersion of the resonator can be modified to exhibit slow or fast light. Because of the relatively slow thermal response, we observe the bandwidth of the PTIT to be $2\pi\times15.9$ Hz which theoretically suggests a group velocity of as low as $5$ m/s.


title: Trapping $^{171}$Yb Atoms into a One-Dimensional Optical Lattice with a Small Waist

abstract: In most experiments with atoms trapped in optical lattices, the transverse size of the optical lattice beams is on the order of tens of micrometers, and loading many atoms into smaller optical lattices has not been carefully investigated. We report trapping 1500 $^{171}$Yb atoms in a one-dimensional optical lattice generated by a narrow cavity mode at a distance of 0.14 mm from a mirror surface. The simplest approach of loading atoms from a mirror magneto-optical trap overlapped with the cavity mode allows the adjustment of the loading position by tuning a uniform bias magnetic field. The number of atoms trapped in the optical lattice exhibits two local maxima for different lattice depths, with a global maximum in the deeper lattice. These results open a way to quantum mechanical manipulation of atoms based on strong interaction with a tightly focused light field.


title: Complete Dictionary Learning via $\ell_p$-norm Maximization

abstract: Dictionary learning is a classic representation learning method that has been widely applied in signal processing and data analytics. In this paper, we investigate a family of $\ell_p$-norm ($p>2,p \in \mathbb{N}$) maximization approaches for the complete dictionary learning problem from theoretical and algorithmic aspects. Specifically, we prove that the global maximizers of these formulations are very close to the true dictionary with high probability, even when Gaussian noise is present. Based on the generalized power method (GPM), an efficient algorithm is then developed for the $\ell_p$-based formulations. We further show the efficacy of the developed algorithm: for the population GPM algorithm over the sphere constraint, it first quickly enters the neighborhood of a global maximizer, and then converges linearly in this region. Extensive experiments will demonstrate that the $\ell_p$-based approaches enjoy a higher computational efficiency and better robustness than conventional approaches and $p=3$ performs the best.


title: Using dark states to charge and stabilise open quantum batteries

abstract: We introduce an open quantum battery protocol using dark states to achieve both superextensive capacity and power density, with non-interacting spins coupled to a reservoir. Further, our power density actually scales with the of number of spins $N$ in the battery. We show that the enhanced capacity and power is correlated with entanglement. Whilst connected to the charger, the charged state of the battery is a steady state, stabilized through quantum interference in the open system.


title: Optimal Advertising for Information Products

abstract: When selling information products, the seller can provide some free partial information to change people's valuations so that the overall revenue can possibly be increased. We study the general problem of advertising information products by revealing partial information. We consider buyers who are decision-makers. The outcomes of the decision problems depend on the state of the world that is unknown to the buyers. The buyers can make their own observations and thus can hold different personal beliefs about the state of the world. There is an information seller who has access to the state of the world. The seller can promote the information by revealing some partial information. We assume that the seller chooses a long-term advertising strategy and then commits to it. The seller's goal is to maximize the expected revenue. We study the problem in two settings. (1) The seller targets buyers of a certain type. In this case, finding the optimal advertising strategy is equivalent to finding the concave closure of a simple function. The function is a product of two quantities, the likelihood ratio and the cost of uncertainty. Based on this observation, we prove some properties of the optimal mechanism, which allow us to solve for the optimal mechanism by a finite-size convex program. The convex program will have a polynomial-size if the state of the world has a constant number of possible realizations or the buyers face a decision problem with a constant number of options. For the general problem, we prove that it is NP-hard to find the optimal mechanism. (2) When the seller faces buyers of different types and only knows the distribution of their types, we provide an approximation algorithm when it is not too hard to predict the possible type of buyers who will make the purchase. For the general problem, we prove that it is NP-hard to find a constant-factor approximation.


title: Permutation Inference for Canonical Correlation Analysis

abstract: Canonical correlation analysis (CCA) has become a key tool for population neuroimaging, allowing investigation of associations between many imaging and non-imaging measurements. As other variables are often a source of variability not of direct interest, previous work has used CCA on residuals from a model that removes these effects, then proceeded directly to permutation inference. We show that such a simple permutation test leads to inflated error rates. The reason is that residualisation introduces dependencies among the observations that violate the exchangeability assumption. Even in the absence of nuisance variables, however, a simple permutation test for CCA also leads to excess error rates for all canonical correlations other than the first. The reason is that a simple permutation scheme does not ignore the variability already explained by previous canonical variables. Here we propose solutions for both problems: in the case of nuisance variables, we show that transforming the residuals to a lower dimensional basis where exchangeability holds results in a valid permutation test; for more general cases, with or without nuisance variables, we propose estimating the canonical correlations in a stepwise manner, removing at each iteration the variance already explained, while dealing with different number of variables in both sides. We also discuss how to address the multiplicity of tests, proposing an admissible test that is not conservative, and provide a complete algorithm for permutation inference for CCA.

categories: stat.ME math.ST stat.AP stat.CO stat.ML stat.TH


title: Parallel Clique Counting and Peeling Algorithms

abstract: We present a new parallel algorithm for $k$-clique counting/listing that has polylogarithmic span (parallel time) and is work-efficient (matches the work of the best sequential algorithm) for sparse graphs. Our algorithm is based on computing low out-degree orientations, which we present new linear-work and polylogarithmic-span algorithms for computing in parallel. We also present new parallel algorithms for producing unbiased estimations of clique counts using graph sparsification. Finally, we design two new parallel work-efficient algorithms for approximating the $k$-clique densest subgraph, the first of which is a $1/k$-approximation and the second of which is a $1/(k(1+\epsilon))$-approximation and has polylogarithmic span. Our first algorithm does not have polylogarithmic span, but we prove that it solves a P-complete problem. In addition to the theoretical results, we also implement the algorithms and propose various optimizations to improve their practical performance. On a 30-core machine with two-way hyper-threading, our algorithms achieve 13.23--38.99x and 1.19--13.76x self-relative parallel speedup for $k$-clique counting and $k$-clique densest subgraph, respectively. Compared to the state-of-the-art parallel $k$-clique counting algorithms, we achieve up to 9.88x speedup, and compared to existing implementations of $k$-clique densest subgraph, we achieve up to 11.83x speedup. We are able to compute the $4$-clique counts on the largest publicly-available graph with over two hundred billion edges for the first time.


title: Nuclear multipole responses from chiral effective field theory interaction

abstract: We probe nuclear multipole resonances in the framework of the random-phase approximation by using the interaction obtained from the chiral effective field theory. The three-nucleon force is included in a form of the in-medium two-nucleon interaction which was derived from the chiral three-nucleon force. The isoscalar monopole, isoscalar dipole, isovector dipole and isoscalar quadrupole resonances of the closed-shell $^{56,68,78}$Ni have been investigated. The calculations reasonably reproduce the experimental multipole resonances of $^{56,68}$Ni, and well describe the pygmy dipole resonance and dipole polarizability measured in $^{68}$Ni. The multipole resonances of $^{78}$Ni, including pygmy dipole resonance and dipole polarizability, are predicted. The detailed effects of the tensor force and three-body force are analyzed by dissecting the chiral interaction. We find that in general the tensor force effect on electric giant resonances is not as significant as the effect from the three-body force, although the tensor force provides more than half of the binding energy. The effect from three-body force is strong in light nuclei. Particularly, three-body force is crucial for the formation of the pygmy resonance in calculations.


title: Relativistic double-ionization equation-of-motion coupled-cluster method: Application to low-lying doubly ionized states

abstract: The article deals with the extension of the relativistic double-ionization equation-of-motion coupled-cluster (DI-EOMCC) method [H. Pathak et al. Phys. Rev. A 90, 010501(R) (2014)] for the molecular systems. The Dirac-Coulomb (DC) Hamiltonian with four-component spinors is considered to take care of the relativistic effects. The implemented method is employed to compute a few low-lying doubly ionized states of noble gas atoms (Ar, Kr, Xe, and Rn) and Cl 2 , Br 2 , HBr, and HI. Additionally, we presented results with two intermediate schemes in the four-component relativistic DI-EOMCC framework to understand the role of electron correlation. The computed double ionization spectra for the atomic systems are compared with the values from the non-relativistic DI-EOMCC method with spin-orbit coupling (SOC) [Z. Wang et al. J. Chem. Phys. 142, 144109 (2015)] and the values from the National Institute of Science and Technology (NIST) database. Our atomic results are found to be in good agreement with the NIST values. Further, the obtained results for the molecular systems agree well with the available experimental values.


title: Higher order Massey products and applications

abstract: In this survey, we discuss two research areas related to Massey's higher operations. The first direction is connected with the cohomology of Lie algebras and the theory of representations. The second main theme is at the intersection of toric topology, homotopy theory of polyhedral products, and the homology theory of local rings, Stanley-Reisner rings of simplicial complexes.


title: On Seminal HEDP Research Opportunities Enabled by Colocating Multi-Petawatt Laser with High-Density Electron Beams

abstract: The scientific community is currently witnessing an expensive and worldwide race to achieve the highest possible light intensity. Within the next decade this effort is expected to reach nearly $10^{24}\,\mathrm{W}/\mathrm{cm^2}$ in the lab frame by focusing of 100 PW, near-infrared lasers. A major driving force behind this effort is the possibility to study strong-field vacuum breakdown and an accompanying electron-positron pair plasma via a quantum electrodynamic (QED) cascade [Edwin Cartlidge, "The light fantastic", Science 359, 382 (2018)]. Whereas Europe is focusing on all-optical 10 PW-class laser facilities (e.g., Apollon and ELI), China is already planning on co-locating a 100 PW laser system with a 25 keV superconducting XFEL and thus implicitly also a high-quality electron beam [Station of Extreme Light (SEL) at the Shanghai Superintense-Ultrafast Lasers Facility (SULF)]. This white paper elucidates the seminal scientific opportunities facilitated by colliding dense, multi-GeV electron beams with multi-PW optical laser pulses. Such a multi-beam facility would enable the experimental exploration of extreme HEDP environments by generating electron-positron pair plasmas with unprecedented densities and temperatures, where the interplay between strong-field quantum and collective plasma effects becomes decisive.


title: Optomechanically induced carrier-envelope-phase dependent effects and their analytical solutions

abstract: To date, investigations of carrier-envelope-phase (CEP) dependent effects have been limited to optical pulses with few cycles and high intensity, and have not been reported for other types of pulses. Optomechanical systems are shown to have the potential to go beyond these limits. We present an approach using optomechanics to extend the concept of the traditional CEP in the few-cycle regime to mechanical pulses and develop a two-step model to give a physical insight. By adding an auxiliary continuous optical field, we show that a CEP-dependent effect appears even in the multi-cycle regime of mechanical pulses. We obtain the approximated analytical solutions providing full understanding for these optomechanically induced CEP-dependent effects. In addition, our findings show that one can draw on the optomechanical interaction to revive the CEP-dependent effects on optical pulses with an arbitrary number of cycles and without specific intensity requirements. The effects of CEP, broadly extended to encompass few- and multi-cycle optical and mechanical pulses, may stimulate a variety of applications in the preparation of a CEP-stabilized pulse, the generation of ultrasonic pulses with a desired shape, the linear manipulation of optical combs, and more.


title: Effective statistical fringe removal algorithm for high-sensitivity imaging of ultracold atoms

abstract: High-sensitivity imaging of ultracold atoms is often challenging when interference patterns are imprinted on the imaging light. Such image noises result in low signal-to-noise ratio and limit the capability to extract subtle physical quantities. Here we demonstrate an advanced fringe removal algorithm for absorption imaging of ultracold atoms, which efficiently suppresses unwanted fringe patterns using a small number of sample images without taking additional reference images. The protocol is based on an image decomposition and projection method with an extended image basis. We apply this scheme to raw absorption images of degenerate Fermi gases for the measurement of atomic density fluctuations and temperatures. The quantitative analysis shows that image noises can be efficiently removed with only tens of reference images, which manifests the efficiency of our protocol. Our algorithm would be of particular interest for the quantum emulation experiments in which several physical parameters need to be scanned within a limited time duration.

categories: cond-mat.quant-gas eess.IV physics.optics quant-ph


title: Metric topology on the moduli space

abstract: We define the smooth Lipschitz topology on the moduli space and show that each conformal class is dense in the moduli space endowed with Gromov-Hausdorff topology, which offers an answer to the Tuschmann's question.


title: Ensuring Privacy in Location-Based Services: A Model-based Approach

abstract: In recent years, the widespread of mobile devices equipped with GPS and communication chips has led to the growing use of location-based services (LBS) in which a user receives a service based on his current location. The disclosure of user's location, however, can raise serious concerns about user privacy in general, and location privacy in particular which led to the development of various location privacy-preserving mechanisms aiming to enhance the location privacy while using LBS applications. In this paper, we propose to model the user mobility pattern and utility of the LBS as a Markov decision process (MDP), and inspired by probabilistic current state opacity notation, we introduce a new location privacy metric, namely $\epsilon-$privacy, that quantifies the adversary belief over the user's current location. We exploit this dynamic model to design a LPPM that while it ensures the utility of service is being fully utilized, independent of the adversary prior knowledge about the user, it can guarantee a user-specified privacy level can be achieved for an infinite time horizon. The overall privacy-preserving framework, including the construction of the user mobility model as a MDP, and design of the proposed LPPM, are demonstrated and validated with real-world experimental data.


title: The distribution of the generalized greatest common divisor and visibility of lattice points

abstract: For a fixed $b\in \mathbb{N}=\{1,2,3,\dots\}$, Goins et al. \cite{Harris} defined the concept of $b$-visibility for a lattice point $(r,s)$ in $L=\mathbb{N}\times \mathbb{N}$ which states that $(r,s)$ is $b$-visible from the origin if it lies on the graph of $f(x)=ax^b$, for some positive $a\in \mathbb{Q}$, and no other lattice point in $L$ lies on this graph between $(0,0)$ and $(r,s)$. Furthermore, to study the density of $b$-visible points in $L$ Goins et al. defined a generalization of greatest common divisor, denoted by $\gcd_b$, and proved that the proportion of $b$-visible lattice points in $L$ is given by $1/\zeta(b+1)$, where $\zeta(s)$ is the Riemann zeta function. In this paper we study the mean values of arithmetic functions $\Lambda:L\to \mathbb{ C}$ defined using $\gcd_b$ and recover the main result of \cite{Harris} as a consequence of the more general results of this paper. We also investigate a generalization of a result in \cite{Harris} that asserts that there are arbitrarily large rectangular arrangements of $b$-visible points in the lattice $L$ for a fixed $b$, more specifically, we give necessary and sufficient conditions for an arbitrary rectangular arrangement containing $b$-visible and $b$-invisible points to be realizable in the lattice $L$. Our result is inspired by the work of Herzog and Stewart \cite{Herzog} who proved this in the case $b=1$.


title: Application of a Modified Spheromak Model to Simulations of Coronal Mass Ejection in the Inner Heliosphere

abstract: The magnetic fields of interplanetary coronal mass ejections (ICMEs), which originate close to the Sun in the form of a flux rope, determine their geoeffectiveness. Therefore, robust flux rope-based models of CMEs are required to perform magnetohydrodynamic (MHD) simulations aimed at space weather predictions. We propose a modified spheromak model and demonstrate its applicability to CME simulations. In this model, such properties of a simulated CME as the poloidal and toroidal magnetic fluxes, and the helicity sign can be controlled with a set of input parameters. We propose a robust technique for introducing CMEs with an appropriate speed into a background, MHD solution describing the solar wind in the inner heliosphere. Through a parametric study, we find that the speed of a CME is much more dependent on its poloidal flux than on the toroidal flux. We also show that the CME speed increases with its total energy, giving us control over its initial speed. We further demonstrate the applicability of this model to simulations of CME-CME collisions. Finally, we use this model to simulate the 12 July 2012 CME and compare the plasma properties at 1 AU with observations. The predicted CME properties agree reasonably with observational data.


title: A mixed version for a Fuchs' Lemma

abstract: We prove a version for mixed groups for a Fuchs' result about connections between the cancellation property of a group and the unit lifting property of its (Walk-)endomorphism rings.


title: Cooperative Adaptive Learning Control for A Group of Nonholonomic UGVs by Output Feedback

abstract: A high-gain observer-based cooperative deterministic learning (CDL) control algorithm is proposed in this chapter for a group of identical unicycle-type unmanned ground vehicles (UGVs) to track over desired reference trajectories. For the vehicle states, the positions of the vehicles can be measured, while the velocities are estimated using the high-gain observer. For the trajectory tracking controller, the radial basis function (RBF) neural network (NN) is used to online estimate the unknown dynamics of the vehicle, and the NN weight convergence and estimation accuracy is guaranteed by CDL. The major challenge and novelty of this chapter is to track the reference trajectory using this observer-based CDL algorithm without the full knowledge of the vehicle state and vehicle model. In addition, any vehicle in the system is able to learn the knowledge of unmodeled dynamics along the union of trajectories experienced by all vehicle agents, such that the learned knowledge can be re-used to follow any reference trajectory defined in the learning phase. The learning-based tracking convergence and consensus learning results, as well as using learned knowledge for tracking experienced trajectories, are shown using the Lyapunov method. Simulation is given to show the effectiveness of this algorithm.


title: Handling the Positive-Definite Constraint in the Bayesian Learning Rule

abstract: The Bayesian learning rule is a natural-gradient variational inference method, which not only contains many existing learning algorithms as special cases but also enables the design of new algorithms. Unfortunately, when variational parameters lie in an open constraint set, the rule may not satisfy the constraint and requires line-searches which could slow down the algorithm. In this work, we address this issue for positive-definite constraints by proposing an improved rule that naturally handles the constraints. Our modification is obtained by using Riemannian gradient methods, and is valid when the approximation attains a \emph{block-coordinate natural parameterization} (e.g., Gaussian distributions and their mixtures). We propose a principled way to derive Riemannian gradients and retractions from scratch. Our method outperforms existing methods without any significant increase in computation. Our work makes it easier to apply the rule in the presence of positive-definite constraints in parameter spaces.


title: Omni-Scale CNNs: a simple and effective kernel size configuration for time series classification

abstract: The Receptive Field (RF) size has been one of the most important factors for One Dimensional Convolutional Neural Networks (1D-CNNs) on time series classification tasks. Large efforts have been taken to choose the appropriate size because it has a huge influence on the performance and differs significantly for each dataset. In this paper, we propose an Omni-Scale block (OS-block) for 1D-CNNs, where the kernel sizes are decided by a simple and universal rule. Particularly, it is a set of kernel sizes that can efficiently cover the best RF size across different datasets via consisting of multiple prime numbers according to the length of the time series. The experiment result shows that models with the OS-block can achieve a similar performance as models with the searched optimal RF size and due to the strong optimal RF size capture ability, simple 1D-CNN models with OS-block achieves the state-of-the-art performance on four time series benchmarks, including both univariate and multivariate data from multiple domains. Comprehensive analysis and discussions shed light on why the OS-block can capture optimal RF sizes across different datasets. Code available [https://github.com/Wensi-Tang/OS-CNN]


title: Reduction of multisymplectic manifolds

abstract: We extend the Marsden-Weinstein-Meyer symplectic reduction theorem to the setting of multisymplectic manifolds. In this context, we investigate the dependence of the reduced space on the reduction parameters. With respect to a distinguished class of multisymplectic moment maps, an exact stationary phase approximation and nonabelian localization theorem are also obtained.


title: A Paley-Wiener theorem for spherical $p$-adic spaces and Bernstein morphisms

abstract: Let $G$ be (the rational points of) a connected reductive group over a local non-archimedean field $F$. In this article we formulate and prove a property of an $F$-spherical homogeneous $G$-space (which in addition satisfies the finite multiplicity property, which is expected to hold for all $F$-spherical homogeneous $G$-spaces) which we call the Paley-Wiener property. This is much more elementary, but also contains much less information, than the recent relevant work of Delorme, Harinck and Sakellaridis (however, it holds for a wider class of spaces). The property results from a parallel categorical property. We also discuss how to define Bernstein morphisms via this approach.


title: Exploring the Connection Between Binary and Spiking Neural Networks

abstract: On-chip edge intelligence has necessitated the exploration of algorithmic techniques to reduce the compute requirements of current machine learning frameworks. This work aims to bridge the recent algorithmic progress in training Binary Neural Networks and Spiking Neural Networks - both of which are driven by the same motivation and yet synergies between the two have not been fully explored. We show that training Spiking Neural Networks in the extreme quantization regime results in near full precision accuracies on large-scale datasets like CIFAR-$100$ and ImageNet. An important implication of this work is that Binary Spiking Neural Networks can be enabled by "In-Memory" hardware accelerators catered for Binary Neural Networks without suffering any accuracy degradation due to binarization. We utilize standard training techniques for non-spiking networks to generate our spiking networks by conversion process and also perform an extensive empirical analysis and explore simple design-time and run-time optimization techniques for reducing inference latency of spiking networks (both for binary and full-precision models) by an order of magnitude over prior work.


title: Stochastic Model Predictive Control for Central HVAC Plants

abstract: We present a stochastic model predictive control (MPC) framework for central heating, ventilation, and air conditioning (HVAC) plants. The framework uses real data to forecast and quantify uncertainty of disturbances affecting the system over multiple timescales (electrical loads, heating/cooling loads, and energy prices). We conduct detailed closed-loop simulations and systematic benchmarks for the central HVAC plant of a typical university campus. Results demonstrate that deterministic MPC fails to properly capture disturbances and that this translates into economic penalties associated with peak demand charges and constraint violations in thermal storage capacity (overflow and/or depletion). Our results also demonstrate that stochastic MPC provides a more systematic approach to mitigate uncertainties and that this ultimately leads to cost savings of up to 7.5% and to mitigation of storage constraint violations. Benchmark results also indicate that these savings are close to ideal savings (9.6%) obtained under MPC with perfect information.


title: Causal Strategic Linear Regression

abstract: In many predictive decision-making scenarios, such as credit scoring and academic testing, a decision-maker must construct a model that accounts for agents' propensity to "game" the decision rule by changing their features so as to receive better decisions. Whereas the strategic classification literature has previously assumed that agents' outcomes are not causally affected by their features (and thus that strategic agents' goal is deceiving the decision-maker), we join concurrent work in modeling agents' outcomes as a function of their changeable attributes. As our main contribution, we provide efficient algorithms for learning decision rules that optimize three distinct decision-maker objectives in a realizable linear setting: accurately predicting agents' post-gaming outcomes (prediction risk minimization), incentivizing agents to improve these outcomes (agent outcome maximization), and estimating the coefficients of the true underlying model (parameter estimation). Our algorithms circumvent a hardness result of Miller et al. (2020) by allowing the decision maker to test a sequence of decision rules and observe agents' responses, in effect performing causal interventions through the decision rules.


title: Spatially mixed moir\'e excitons in two-dimensional van der Waals superlattices

abstract: Moir\'e superlattices open an unprecedented opportunity for tailoring interactions between quantum particles and their coupling to electromagnetic fields. Strong superlattice potential generates moir\'e minibands of excitons -- bound pairs of electrons and holes that reside either in a single layer (intralayer excitons) or two separate layers (interlayer excitons). The twist-angle-controlled interlayer hybridization of carriers can also mix the two types of excitons to combine the strengths of both. Here, we report a direct observation of spatially mixed moir\'e excitons in angle-aligned WSe2/WS2 and MoSe2/WS2 superlattices by optical reflectance spectroscopy. The strongly interacting interlayer and intralayer moir\'e excitons in WSe2/WS2 manifest energy level anticrossing and oscillator strength redistribution under a vertical electric field. We also observe doping-dependent exciton miniband renormalization and mixing near half filling of the first electron miniband of WS2. Our findings have significant implications for emerging correlated states in two-dimensional semiconductors, such as exciton condensates and Bose-Hubbard models, and optoelectronic applications of these materials.


title: Real-time Exciton-Manipulation of Plexcitonic Coupling in a WS2-Ag nanocavity

abstract: We demonstrate the real-time exciton-manipulation of plexcitonic coupling in monolayer WS2 coupled to a plasmonic nanocavity by immersing into a mixed solution of dichloromethane (DCM) and ethanol. By adjusting the mixture ratio, a continuous tuning of the Rabi splitting energy ranged from 178 meV (in ethanol) to 266 meV (in DCM) is achieved. The results are mainly attributed to the remarkable increase of the proportion of neutral exciton in the monolayer WS2 (from 59% to 100%) as the concentration of DCM is increased. It offers an important stepping stone towards a further study of plexcitonic coupling in layered materials, along with potential applications in quantum information processing and nonlinear optical materials.


title: Robust Learning-Based Control via Bootstrapped Multiplicative Noise

abstract: Despite decades of research and recent progress in adaptive control and reinforcement learning, there remains a fundamental lack of understanding in designing controllers that provide robustness to inherent non-asymptotic uncertainties arising from models estimated with finite, noisy data. We propose a robust adaptive control algorithm that explicitly incorporates such non-asymptotic uncertainties into the control design. The algorithm has three components: (1) a least-squares nominal model estimator; (2) a bootstrap resampling method that quantifies non-asymptotic variance of the nominal model estimate; and (3) a non-conventional robust control design method using an optimal linear quadratic regulator (LQR) with multiplicative noise. A key advantage of the proposed approach is that the system identification and robust control design procedures both use stochastic uncertainty representations, so that the actual inherent statistical estimation uncertainty directly aligns with the uncertainty the robust controller is being designed against. We show through numerical experiments that the proposed robust adaptive controller can significantly outperform the certainty equivalent controller on both expected regret and measures of regret risk.

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


title: An Overlapping Domain Decomposition Framework without Dual Formulation for Variational Imaging Problems

abstract: In this paper, we propose a novel overlapping domain decomposition method that can be applied to various problems in variational imaging such as total variation minimization. Most of recent domain decomposition methods for total variation minimization adopt the Fenchel--Rockafellar duality, whereas the proposed method is based on the primal formulation. Thus, the proposed method can be applied not only to total variation minimization but also to those with complex dual problems such as higher order models. In the proposed method, an equivalent formulation of the model problem with parallel structure is constructed using a custom overlapping domain decomposition scheme with the notion of essential domains. As a solver for the constructed formulation, we propose a decoupled augmented Lagrangian method for untying the coupling of adjacent subdomains. Convergence analysis of the decoupled augmented Lagrangian method is provided. We present implementation details and numerical examples for various model problems including total variation minimizations and higher order models.


title: Black-box sampling for weakly smooth Langevin Monte Carlo using p-generalized Gaussian smoothing

abstract: Discretization of continuous-time diffusion processes is a widely recognized method for sampling. However, the canonical Euler-Maruyama discretization of the Langevin diffusion process, also named as Langevin Monte Carlo (LMC), studied mostly in the context of smooth (gradient-Lipschitz) and strongly log-concave densities, a significant constraint for its deployment in many sciences, including computational statistics and statistical learning. In this paper, we establish several theoretical contributions to the literature on such sampling methods. Particularly, we generalize the Gaussian smoothing, approximate the gradient using p-generalized Gaussian smoothing and take advantage of it in the context of black-box sampling. We first present a non-strongly concave and weakly smooth black-box LMC algorithm, ideal for practical applicability of sampling challenges in a general setting.


title: Reconfigurable Intelligent Surface Assisted Multiuser MISO Systems Exploiting Deep Reinforcement Learning

abstract: Recently, the reconfigurable intelligent surface (RIS), benefited from the breakthrough on the fabrication of programmable meta-material, has been speculated as one of the key enabling technologies for the future six generation (6G) wireless communication systems scaled up beyond massive multiple input multiple output (Massive-MIMO) technology to achieve smart radio environments. Employed as reflecting arrays, RIS is able to assist MIMO transmissions without the need of radio frequency chains resulting in considerable reduction in power consumption. In this paper, we investigate the joint design of transmit beamforming matrix at the base station and the phase shift matrix at the RIS, by leveraging recent advances in deep reinforcement learning (DRL). We first develop a DRL based algorithm, in which the joint design is obtained through trial-and-error interactions with the environment by observing predefined rewards, in the context of continuous state and action. Unlike the most reported works utilizing the alternating optimization techniques to alternatively obtain the transmit beamforming and phase shifts, the proposed DRL based algorithm obtains the joint design simultaneously as the output of the DRL neural network. Simulation results show that the proposed algorithm is not only able to learn from the environment and gradually improve its behavior, but also obtains the comparable performance compared with two state-of-the-art benchmarks. It is also observed that, appropriate neural network parameter settings will improve significantly the performance and convergence rate of the proposed algorithm.


title: Electron interactions in strain-induced zero-energy flat band in twisted bilayer graphene near the magic angle

abstract: In the vicinity of the magic angle in twisted bilayer graphene (TBG), the two low-energy van Hove singularities (VHSs) become exceedingly narrow1-10 and many exotic correlated states, such as superconductivity, ferromagnetism, and topological phases, are observed11-16. Heterostrain, which is almost unavoidable in the TBG, can modify its single-particle band structure and lead to novel properties of the TBG that have never been considered so far. Here, we show that heterostrain in a TBG near the magic angle generates a new zero-energy flat band between the two VHSs. Doping the TBG to partially fill the zero-energy flat band, we observe a correlation-induced gap of about 10 meV that splits the flat band. By applying perpendicular magnetic fields, a large and linear response of the gap to magnetic fields is observed, attributing to the emergence of large orbital magnetic moments in the TBG when valley degeneracy of the flat band is lifted by electron-electron interactions. The orbital magnetic moment per moire supercell is measured as about 15 uB in the TBG.


title: Radiative topological biphoton states in modulated qubit arrays

abstract: We study topological properties of bound pairs of photons in spatially-modulated qubit arrays (arrays of two-level atoms) coupled to a waveguide. While bound pairs behave like Bloch waves, they are topologically nontrivial in the parameter space formed by the center-of-mass momentum and the modulation phase, where the latter plays the role of a synthetic dimension. In a superlattice where each unit cell contains three two-level atoms (qubits), we calculate the Chern numbers for the bound-state photon bands, which are found to be $(1,-2,1)$. For open boundary condition, we find exotic topological bound-pair edge states with radiative losses. Unlike the conventional case of the bulk-edge correspondence, these novel edge modes not only exist in gaps separating the bound-pair bands, but they also may merge with and penetrate into the bands. By joining two structures with different spatial modulations, we find long-lived interface states which may have applications in storage and quantum information processing.


title: Singlet Fermionic Dark Matter with Dark $Z$

abstract: We present a fermionic dark matter model mediated by the hidden gauge boson. We assume the QED-like hidden sector which consists of a Dirac fermion and U(1)$_X$ gauge symmetry, and introduce an additional scalar electroweak doublet field with the U(1)$_X$ charge as a mediator. The hidden U(1)$_X$ symmetry is spontaneously broken by the electroweak symmetry breaking and there exists a massive extra neutral gauge boson in this model which is the mediator between the hidden and visible sectors. Due to the U(1)$_X$ charge, the additional scalar doublet does not couple to the Standard Model fermions, which leads to the Higgs sector of type I two Higgs doublet model. The new gauge boson couples to the Standard Model fermions with couplings proportional to those of the ordinary $Z$ boson but very suppressed, thus we call it the dark $Z$ boson. We study the phenomenology of the dark $Z$ boson and the Higgs sector, and show the hidden fermion can be the dark matter candidate.


title: MCG+07-20-052: Interacting dwarf pair in a group environment

abstract: We present an observational study of the interacting pair of dwarf galaxies, MCG+07-20-052, in the vicinity of Milky Way mass spiral galaxy NGC 2998. MCG+07-20-052 is located at a sky-projected distance of 105 kpc from NGC 2998 and the two have a relative line-of-sight velocity of 60 kms. We observed tidal tail-like extensions on both members (D1 and D2) of the interacting pair MCG+07-20-052. The interacting dwarf galaxies, D1 and D2, have B-band absolute magnitudes of $-$17.17 and $-$17.14 mag, respectively, and D2 is significantly bluer than D1. We obtained HI 21 cm line data of the NGC 2998 system using the Giant Metrewave Radio Telescope (GMRT) to get a more detailed view of the neutral hydrogen (HI) emission in the interacting dwarf galaxies and in the galaxy members of the NGC 2998 group. Evidence of a merger between the dwarf galaxies in the MCG+07-20-052 is also present in the HI kinematics and morphology where we find that the HI is mostly concentrated around D2, which also shows a higher level of star-forming activity and bluer $g-r$ color index compared to D1. In addition, we detect extended tenuous HI emission around another member galaxy, NGC 3006, located close to the MCG+07-20-052-pair at a sky-projected distance of 41 kpc. We compare here our results from the MCG+07-20-052 pair-NGC 2998 system with other known LMC-SMC-Milky Way type systems and discuss the possible origin of the dwarf-dwarf interaction.


title: Approximate Data Deletion from Machine Learning Models

abstract: Deleting data from a trained machine learning (ML) model is a critical task in many applications. For example, we may want to remove the influence of training points that might be out of date or outliers. Regulations such as EU's General Data Protection Regulation also stipulate that individuals can request to have their data deleted. The naive approach to data deletion is to retrain the ML model on the remaining data, but this is too time consuming. In this work, we propose a new approximate deletion method for linear and logistic models whose computational cost is linear in the the feature dimension $d$ and independent of the number of training data $n$. This is a significant gain over all existing methods, which all have superlinear time dependence on the dimension. We also develop a new feature-injection test to evaluate the thoroughness of data deletion from ML models.


title: On Hiding Neural Networks Inside Neural Networks

abstract: Modern neural networks often contain significantly more parameters than the size of their training data. We show that this excess capacity provides an opportunity for embedding secret machine learning models within a trained neural network. Our novel framework hides the existence of a secret neural network with arbitrary desired functionality within a carrier network. We prove theoretically that the secret network's detection is computationally infeasible and demonstrate empirically that the carrier network does not compromise the secret network's disguise. Our paper introduces a previously unknown steganographic technique that can be exploited by adversaries if left unchecked.


title: A Distributed Architecture for Real-time Hybrid Traffic Light Control in Urban Transportation Networks

abstract: A macroscopic model is proposed to depict the traffic dynamics involved in urban traffic systems. The link dynamics are described based on the cell-transmission model and bounded by the link capacities, while the flow dynamics are proposed based on the discharge headways and saturation flow at intersections. To fulfill the requirement of a closed-loop traffic light control strategy, an approach to estimate the branching ratios at intersections is proposed and simulations show that the convergence would be achieved under constant cyclic flow profiles. Furthermore, a system partitioning approach is proposed based congestion level identification, which is achieved via a machine learning method and a hybrid traffic network control strategy is proposed to integrate different traffic light control schemes together.


title: Sparse Optimization for Green Edge AI Inference

abstract: With the rapid upsurge of deep learning tasks at the network edge, effective edge artificial intelligence (AI) inference becomes critical to provide low-latency intelligent services for mobile users via leveraging the edge computing capability. In such scenarios, energy efficiency becomes a primary concern. In this paper, we present a joint inference task selection and downlink beamforming strategy to achieve energy-efficient edge AI inference through minimizing the overall power consumption consisting of both computation and transmission power consumption, yielding a mixed combinatorial optimization problem. By exploiting the inherent connections between the set of task selection and group sparsity structural transmit beamforming vector, we reformulate the optimization as a group sparse beamforming problem. To solve this challenging problem, we propose a log-sum function based three-stage approach. By adopting the log-sum function to enhance the group sparsity, a proximal iteratively reweighted algorithm is developed. Furthermore, we establish the global convergence analysis and provide the ergodic worst-case convergence rate for this algorithm. Simulation results will demonstrate the effectiveness of the proposed approach for improving energy efficiency in edge AI inference systems.


title: Toward a mathematical theory of the crystallographic phase retrieval problem

abstract: Motivated by the X-ray crystallography technology to determine the atomic structure of biological molecules, we study the crystallographic phase retrieval problem, arguably the leading and hardest phase retrieval setup. This problem entails recovering a K-sparse signal of length N from its Fourier magnitude or, equivalently, from its periodic auto-correlation. Specifically, this work focuses on the fundamental question of uniqueness: what is the maximal sparsity level K/N that allows unique mapping between a signal and its Fourier magnitude, up to intrinsic symmetries. We design a systemic computational technique to affirm uniqueness for any specific pair (K,N), and establish the following conjecture: the Fourier magnitude determines a generic signal uniquely, up to intrinsic symmetries, as long as K<=N/2. Based on group-theoretic considerations and an additional computational technique, we formulate a second conjecture: if K<N/2, then for any signal the set of solutions to the crystallographic phase retrieval problem has measure zero in the set of all signals with a given Fourier magnitude. Together, these conjectures constitute the first attempt to establish a mathematical theory for the crystallographic phase retrieval problem.


title: Signatures of the vortical quark-gluon plasma in hadron yields

abstract: We investigate the hadron production from the vortical quark-gluon plasma created in heavy-ion collisions. Based on the quark-coalescence and statistical hadronization models, we show that total hadron yields summed over the spin components are enhanced by the local vorticity with quadratic dependence. The enhancement factor amounts to be a few percent and may be detectable within current experimental sensitivities. We also show that the effect is stronger for hadrons with larger spin, and thus propose a new signature of the local vorticity, which may be detected by the yield ratio of distinct hadron species having different spins such as $\phi$ and $\eta'$. The vorticity dependence of hadron yields seems robust, with consistent predictions in both of the hadron production mechanisms for reasonable values of the vorticity strength estimated for heavy-ion collisions.


title: Optimizing High Performance Markov Clustering for Pre-Exascale Architectures

abstract: HipMCL is a high-performance distributed memory implementation of the popular Markov Cluster Algorithm (MCL) and can cluster large-scale networks within hours using a few thousand CPU-equipped nodes. It relies on sparse matrix computations and heavily makes use of the sparse matrix-sparse matrix multiplication kernel (SpGEMM). The existing parallel algorithms in HipMCL are not scalable to Exascale architectures, both due to their communication costs dominating the runtime at large concurrencies and also due to their inability to take advantage of accelerators that are increasingly popular. In this work, we systematically remove scalability and performance bottlenecks of HipMCL. We enable GPUs by performing the expensive expansion phase of the MCL algorithm on GPU. We propose a CPU-GPU joint distributed SpGEMM algorithm called pipelined Sparse SUMMA and integrate a probabilistic memory requirement estimator that is fast and accurate. We develop a new merging algorithm for the incremental processing of partial results produced by the GPUs, which improves the overlap efficiency and the peak memory usage. We also integrate a recent and faster algorithm for performing SpGEMM on CPUs. We validate our new algorithms and optimizations with extensive evaluations. With the enabling of the GPUs and integration of new algorithms, HipMCL is up to 12.4x faster, being able to cluster a network with 70 million proteins and 68 billion connections just under 15 minutes using 1024 nodes of ORNL's Summit supercomputer.


title: Utilizing a null class to restrict decision spaces and defend against neural network adversarial attacks

abstract: Despite recent progress, deep neural networks generally continue to be vulnerable to so-called adversarial examples--input images with small perturbations that can result in changes in the output classifications, despite no such change in the semantic meaning to human viewers. This is true even for seemingly simple challenges such as the MNIST digit classification task. In part, this suggests that these networks are not relying on the same set of object features as humans use to make these classifications. In this paper we examine an additional, and largely unexplored, cause behind this phenomenon--namely, the use of the conventional training paradigm in which the entire input space is parcellated among the training classes. Owing to this paradigm, learned decision spaces for individual classes span excessively large regions of the input space and include images that have no semantic similarity to images in the training set. In this study, we train models that include a null class. That is, models may "opt-out" of classifying an input image as one of the digit classes. During training, null images are created through a variety of methods, in an attempt to create tighter and more semantically meaningful decision spaces for the digit classes. The best performing models classify nearly all adversarial examples as nulls, rather than mistaking them as a member of an incorrect digit class, while simultaneously maintaining high accuracy on the unperturbed test set. The use of a null class and the training paradigm presented herein may provide an effective defense against adversarial attacks for some applications. Code for replicating this study will be made available at https://github.com/mattroos/null_class_adversarial_defense .


title: Temporal Spike Sequence Learning via Backpropagation for Deep Spiking Neural Networks

abstract: Spiking neural networks (SNNs) are well suited for spatio-temporal learning and implementations on energy-efficient event-driven neuromorphic processors. However, existing SNN error backpropagation (BP) methods lack proper handling of spiking discontinuities and suffer from low performance compared with the BP methods for traditional artificial neural networks. In addition, a large number of time steps are typically required to achieve decent performance, leading to high latency and rendering spike-based computation unscalable to deep architectures. We present a novel Temporal Spike Sequence Learning Backpropagation (TSSL-BP) method for training deep SNNs, which breaks down error backpropagation across two types of inter-neuron and intra-neuron dependencies and leads to improved temporal learning precision. It captures inter-neuron dependencies through presynaptic firing times by considering the all-or-none characteristics of firing activities and captures intra-neuron dependencies by handling the internal evolution of each neuronal state in time. TSSL-BP efficiently trains deep SNNs within a much shortened temporal window of a few steps while improving the accuracy for various image classification datasets including CIFAR10.