title: Spectral flow of Dirac operators with magnetic cable knot

abstract: We study the spectral flow of Dirac operators with magnetic links on $\mathbb{S}^3$. These are generalisations of Aharonov-Bohm solenoids where the magnetic fields contain finitely many field lines coinciding with the components of a link, the flux of each exhibiting the same $2\pi$-periodicity as A-B solenoids. We study the spectral flow of the loop obtained as tuning the flux from $0$ to $2\pi$ in the case of only one field line: we relate the spectral flows obtained for one given knot and its cable knots, and obtain that torus knots have trivial spectral flow. The operators are studied in their Coulomb gauge in $\mathbb{R}^3$ (seen as a chart of $\mathbb{S}^3$ through the stereographic projection), which is simply given by the Biot and Savart formula.


title: A Spiking Neural Network with Local Learning Rules Derived From Nonnegative Similarity Matching

abstract: The design and analysis of spiking neural network algorithms will be accelerated by the advent of new theoretical approaches. In an attempt at such approach, we provide a principled derivation of a spiking algorithm for unsupervised learning, starting from the nonnegative similarity matching cost function. The resulting network consists of integrate-and-fire units and exhibits local learning rules, making it biologically plausible and also suitable for neuromorphic hardware. We show in simulations that the algorithm can perform sparse feature extraction and manifold learning, two tasks which can be formulated as nonnegative similarity matching problems.


title: Lower bounds for estimates of the Schr\"odinger maximal function

abstract: We give new lower bounds for $L^p$ estimates of the Schr\"odinger maximal function by generalizing an example of Bourgain.


title: A review on non-relativistic fully numerical electronic structure calculations on atoms and diatomic molecules

abstract: The need for accurate calculations on atoms and diatomic molecules is motivated by the opportunities and challenges of such studies. The most commonly-used approach for all-electron electronic structure calculations in general - the linear combination of atomic orbitals (LCAO) method - is discussed in combination with Gaussian, Slater a.k.a. exponential, and numerical radial functions. Even though LCAO calculations have major benefits, their shortcomings motivate the need for fully numerical approaches based on, e.g. finite differences, finite elements, or the discrete variable representation, which are also briefly introduced. Applications of fully numerical approaches for general molecules are briefly reviewed, and their challenges are discussed. It is pointed out that the high level of symmetry present in atoms and diatomic molecules can be exploited to fashion more efficient fully numerical approaches for these special cases, after which it is possible to routinely perform all-electron Hartree-Fock and density functional calculations directly at the basis set limit on such systems. Applications of fully numerical approaches to calculations on atoms as well as diatomic molecules are reviewed. Finally, a summary and outlook is given.


title: Quantum affine algebras and cluster algebras

abstract: This article is an extended version of the minicourse given by the second author at the summer school of the conference "Interactions of quantum affine algebras with cluster algebras, current algebras and categorification", held in June 2018 in Washington. The aim of the minicourse, consisting of three lectures, was to present a number of results and conjectures on certain monoidal categories of finite-dimensional representations of quantum affine algebras, obtained by exploiting the fact that their Grothendieck rings have the natural structure of a cluster algebra.


title: Think Global, Act Local: The Influence of Environment Age and Host Mass on Type Ia Supernova Light Curves

abstract: The reliability of Type Ia supernovae (SNIa) may be limited by the imprint of their galactic origins. To investigate the connection between supernovae and their host characteristics, we developed an improved method to estimate the stellar population age of the host as well as the local environment around the site of the supernova. We use a Bayesian method to estimate the star formation history and mass weighted age of a supernova's environment by matching observed spectral energy distributions to a synthesized stellar population. Applying this age estimator to both the photometrically and spectroscopically classified Sloan Digital Sky Survey II supernovae (N=103) we find a $0.114 \pm 0.039~{\rm mag}$ `step' in the average Hubble residual at a stellar age of $\sim 8~\text{Gyr}$; it is nearly twice the size of the currently popular mass step. We then apply a principal component analysis on the SALT2 parameters, host stellar mass, and local environment age. We find that a new parameter, PC$_1$, consisting of a linear combination of stretch, host stellar mass, and local age, shows a very significant ($4.7\sigma$) correlation with Hubble residuals. There is a much broader range of PC$_1$ values found in the Hubble flow sample when compared with the Cepheid calibration galaxies. These samples have mildly statistically different average PC$_1$ values, at $\sim 2.5\sigma$, resulting in at most a 1.3% reduction in the evaluation of H$_0$. Despite accounting for the highly significant trend in SNIa Hubble residuals, there remains a 9% discrepancy between the most recent precision estimates of H$_0$ using SNIa and the CMB.


title: On correlation functions in $J\bar T$-deformed CFTs

abstract: The $J\bar T$ deformation, built from the components of the stress tensor and of a $U(1)$ current, is a universal irrelevant deformation of two-dimensional CFTs that preserves the left-moving conformal symmetry, while breaking locality on the right-moving side. Operators in the $J\bar T$-deformed CFT are naturally labeled by the left-moving position and right-moving momentum and transform in representations of the one-dimensional extended conformal group. We derive an all-orders formula for the spectrum of conformal dimensions and charges of the deformed CFT, which we cross-check at leading order using conformal perturbation theory. We also compute the linear corrections to the one-dimensional OPE coefficients and comment on the extent to which the correlation functions in $J\bar T$-deformed CFTs can be obtained from field-dependent coordinate transformations.


title: A Forest from the Trees: Generation through Neighborhoods

abstract: In this work, we propose to learn a generative model using both learned features (through a latent space) and memories (through neighbors). Although human learning makes seamless use of both learned perceptual features and instance recall, current generative learning paradigms only make use of one of these two components. Take, for instance, flow models, which learn a latent space of invertible features that follow a simple distribution. Conversely, kernel density techniques use instances to shift a simple distribution into an aggregate mixture model. Here we propose multiple methods to enhance the latent space of a flow model with neighborhood information. Not only does our proposed framework represent a more human-like approach by leveraging both learned features and memories, but it may also be viewed as a step forward in non-parametric methods. The efficacy of our model is shown empirically with standard image datasets. We observe compelling results and a significant improvement over baselines.


title: Visualization tools for parameter selection in cluster analysis

abstract: We propose an algorithm, HPREF (Hierarchical Partitioning by Repeated Features), that produces a hierarchical partition of a set of clusterings of a fixed dataset, such as sets of clusterings produced by running a clustering algorithm with a range of parameters. This gives geometric structure to such sets of clustering, and can be used to visualize the set of results one obtains by running a clustering algorithm with a range of parameters.


title: Blaze: Simplified High Performance Cluster Computing

abstract: MapReduce and its variants have significantly simplified and accelerated the process of developing parallel programs. However, most MapReduce implementations focus on data-intensive tasks while many real-world tasks are compute intensive and their data can fit distributedly into the memory. For these tasks, the speed of MapReduce programs can be much slower than those hand-optimized ones. We present Blaze, a C++ library that makes it easy to develop high performance parallel programs for such compute intensive tasks. At the core of Blaze is a highly-optimized in-memory MapReduce function, which has three main improvements over conventional MapReduce implementations: eager reduction, fast serialization, and special treatment for a small fixed key range. We also offer additional conveniences that make developing parallel programs similar to developing serial programs. These improvements make Blaze an easy-to-use cluster computing library that approaches the speed of hand-optimized parallel code. We apply Blaze to some common data mining tasks, including word frequency count, PageRank, k-means, expectation maximization (Gaussian mixture model), and k-nearest neighbors. Blaze outperforms Apache Spark by more than 10 times on average for these tasks, and the speed of Blaze scales almost linearly with the number of nodes. In addition, Blaze uses only the MapReduce function and 3 utility functions in its implementation while Spark uses almost 30 different parallel primitives in its official implementation.


title: Outer automorphism groups of graph products: subgroups and quotients

abstract: We show that the outer automorphism groups of graph products of finitely generated abelian groups satisfy the Tits alternative, are residually finite, their so-called Torelli subgroups are finitely generated, and they satisfy a dichotomy between being virtually nilpotent and containing a non-abelian free subgroup that is determined by a graphical condition on the underlying labelled graph. Graph products of finitely generated abelian groups simultaneously generalize right-angled Artin groups (RAAGs) and right-angled Coxter groups (RACGs), providing a common framework for studying these groups. Our results extend a number of known results for the outer automorphism groups of RAAGs and/or RACGs by a variety of authors, including Caprace, Charney, Day, Ferov, Guirardel, Horbez, Minasyan, Vogtmann, Wade, and the current authors.


title: Very Long Term Field of View Prediction for 360-degree Video Streaming

abstract: 360-degree videos have gained increasing popularity in recent years with the developments and advances in Virtual Reality (VR) and Augmented Reality (AR) technologies. In such applications, a user only watches a video scene within a field of view (FoV) centered in a certain direction. Predicting the future FoV in a long time horizon (more than seconds ahead) can help save bandwidth resources in on-demand video streaming while minimizing video freezing in networks with significant bandwidth variations. In this work, we treat the FoV prediction as a sequence learning problem, and propose to predict the target user's future FoV not only based on the user's own past FoV center trajectory but also other users' future FoV locations. We propose multiple prediction models based on two different FoV representations: one using FoV center trajectories and another using equirectangular heatmaps that represent the FoV center distributions. Extensive evaluations with two public datasets demonstrate that the proposed models can significantly outperform benchmark models, and other users' FoVs are very helpful for improving long-term predictions.


title: Radio Source Extraction with ProFound

abstract: In the current era of radio astronomy, continuum surveys observe a multitude of objects with complex morphologies and sizes, and are not limited to observing point sources. Typical radio source extraction software generates catalogues by using Gaussian components to form a model of the emission. This may not be well suited to complicated jet structures and extended emission, particularly in the era of interferometers with a high density of short baselines, that are sensitive to extended emission. In this paper, we investigate how the optically-motivated source detection package ProFound (Robotham et al. 2018) may be used to model radio emission of both complicated and point-like radio sources. We use a combination of observations and simulations to investigate how ProFound compares to other source extractor packages used for radio surveys. We find that ProFound can accurately recover both the flux densities of simulated Gaussian sources as well as extended radio galaxies. ProFound can create models that trace the complicated nature of these extended galaxies, which we show is not necessarily the case with other source extraction software. Our work suggests that our knowledge of the emission from extended radio objects may be both over or under-estimated using traditional software. We suggest that ProFound offers a useful alternative to the fitting of Gaussian components for generating catalogues from current and future radio surveys. Furthermore, ProFound's multi-wavelength capabilities will be useful in investigating radio sources in combination with multiwavelength data.


title: Exponential quasi-ergodicity for processes with discontinuous trajectories

abstract: This paper tackles the issue of establishing a lower-bound on the asymptotic ratio of survival probabilities between two different initial conditions, asymptotically in time for a given Markov process with extinction. Such a comparison is a crucial step in recent techniques for proving exponential convergence to a quasi-stationary distribution. We introduce a weak form of the Harnack inequality as the essential ingredient for such a comparison. This property is actually a consequence of the convergence property we intend to prove. Its complexity appears as the price to pay for the level of flexibility required by our applications. We show in our illustrations how simply and efficiently it can be used nonetheless. As illustrations, we consider two continuous-time processes on R^d that do not satisfy the classical Harnack inequalities, even in a local version. The first one is a piecewise deterministic process while the second is a pure jump process with restrictions on the directions of its jumps.


title: Random Lindblad Dynamics

abstract: We study the mixing behavior of random Lindblad generators with no symmetries, using the dynamical map or propagator of the dissipative evolution. In particular, we determine the long-time behavior of a dissipative form factor, which is the trace of the propagator, and use this as a diagnostic for the existence or absence of a spectral gap in the distribution of eigenvalues of the Lindblad generator. We find that simple generators with a single jump operator are slowly mixing, and relax algebraically in time, due to the closing of the spectral gap in the thermodynamic limit. Introducing additional jump operators or a Hamiltonian opens up a spectral gap which remains finite in the thermodynamic limit, leading to exponential relaxation and thus rapid mixing. We use the method of moments and introduce a novel diagrammatic expansion to determine exactly the form factor to leading order in Hilbert space dimension $N$. We also present numerical support for our main results.

categories: quant-ph cond-mat.mes-hall cond-mat.stat-mech math-ph math.MP


title: Identification and Estimation of Causal Effects from Dependent Data

abstract: The assumption that data samples are independent and identically distributed (iid) is standard in many areas of statistics and machine learning. Nevertheless, in some settings, such as social networks, infectious disease modeling, and reasoning with spatial and temporal data, this assumption is false. An extensive literature exists on making causal inferences under the iid assumption [18, 12, 28, 22], even when unobserved confounding bias may be present. But, as pointed out in [20], causal inference in non-iid contexts is challenging due to the presence of both unobserved confounding and data dependence. In this paper we develop a general theory describing when causal inferences are possible in such scenarios. We use segregated graphs [21], a generalization of latent projection mixed graphs [30], to represent causal models of this type and provide a complete algorithm for non-parametric identification in these models. We then demonstrate how statistical inference may be performed on causal parameters identified by this algorithm. In particular, we consider cases where only a single sample is available for parts of the model due to full interference, i.e., all units are pathwise dependent and neighbors' treatments affect each others' outcomes [26]. We apply these techniques to a synthetic data set which considers users sharing fake news articles given the structure of their social network, user activity levels, and baseline demographics and socioeconomic covariates.


title: "Choix de Bruxelles": A New Operation on Positive Integers

abstract: The "Choix de Bruxelles" operation replaces a positive integer n by any of the numbers that can be obtained by halving or doubling a substring of the decimal representation of n. For example, 16 can become any of 16, 26, 13, 112, 8, or 32. We investigate the properties of this interesting operation and its iterates.


title: Residual Energy Based Cluster-head Selection in WSNs for IoT Application

abstract: Wireless sensor networks (WSN) groups specialized transducers that provide sensing services to Internet of Things (IoT) devices with limited energy and storage resources. Since replacement or recharging of batteries in sensor nodes is almost impossible, power consumption becomes one of the crucial design issues in WSN. Clustering algorithm plays an important role in power conservation for the energy constrained network. Choosing a cluster head can appropriately balance the load in the network thereby reducing energy consumption and enhancing lifetime. The paper focuses on an efficient cluster head election scheme that rotates the cluster head position among the nodes with higher energy level as compared to other. The algorithm considers initial energy, residual energy and an optimum value of cluster heads to elect the next group of cluster heads for the network that suits for IoT applications such as environmental monitoring, smart cities, and systems. Simulation analysis shows the modified version performs better than the LEACH protocol by enhancing the throughput by 60%, lifetime by 66%, and residual energy by 64%.


title: Non-uniqueness of entropy-conserving solutions to the ideal compressible MHD equations

abstract: In this note we consider the ideal compressible magneto-hydrodynamics (MHD) equations in a special two dimensional setting. We show that there exist particular initial data for which one obtains infinitely many entropy-conserving weak solutions by using the convex integration technique. Our result is also true for the isentropic case.


title: Search for rare decay $J/ \psi \to \phi e^+ e^-$

abstract: Using a data sample of $448.1\times10^6$ $\psi(3686)$ events collected at $\sqrt{s}=$ 3.686 GeV with the BESIII detector at the BEPCII, we search for the rare decay $J/\psi \to \phi e^+ e^-$ via $\psi(3686) \to \pi^+\pi^- J/\psi $. No signal events are observed and the upper limit on the branching fraction is set to be $\mathcal{B}(J/\psi \to \phi e^+ e^-) < 1.2 \times 10^{-7}$ at the 90\% confidence level, which is still about one order of magnitude higher than the Standard Model prediction.


title: Factoring semi-primes with (quantum) SAT-solvers

abstract: The assumed computationally difficulty of factoring large integers forms the basis of security for RSA public-key cryptography, which specifically relies on products of two large primes or semi-primes. The best-known factoring algorithms for classical computers run in sub-exponential time. Since integer factorization is in NP, one can reduce this problem to any NP-hard problem, such as Boolean Satisfiability (SAT). While reducing factoring to SAT has proved to be useful for studying SAT solvers, attempting to factor large integers via such a reduction has not been found to be successful. Shor's quantum factoring algorithm factors any integer in polynomial time, although large-scale fault-tolerant quantum computers capable of implementing Shor's algorithm are not yet available, so relevant benchmarking experiments for factoring via Shor's algorithm are not yet possible. In recent years, however, several authors have attempted factorizations with the help of quantum processors via reductions to NP-hard problems. While this approach may shed some light on some algorithmic approaches for quantum solutions to NP-hard problems, in this paper we study and question the practical effectiveness of this approach for factoring large numbers. We find no evidence that this is a viable path toward factoring large numbers, even for scalable fault-tolerant quantum computers, as well as for various quantum annealing or other special purpose quantum hardware.


title: Generalization Bounds For Unsupervised and Semi-Supervised Learning With Autoencoders

abstract: Autoencoders are widely used for unsupervised learning and as a regularization scheme in semi-supervised learning. However, theoretical understanding of their generalization properties and of the manner in which they can assist supervised learning has been lacking. We utilize recent advances in the theory of deep learning generalization, together with a novel reconstruction loss, to provide generalization bounds for autoencoders. To the best of our knowledge, this is the first such bound. We further show that, under appropriate assumptions, an autoencoder with good generalization properties can improve any semi-supervised learning scheme. We support our theoretical results with empirical demonstrations.


title: Linking Scalar Dark Matter and Neutrino Masses with IceCube 170922A

abstract: Two of the key unresolved issues facing Standard Model physics are (i) the appearance of a small but non-zero neutrino mass, and, (ii) the missing mass problem in the Universe. The focus of this paper is a previously proposed low energy effective theory that couples a dark scalar to Standard Model neutrinos. This provides a stable dark matter candidate as well as radiatively generating a neutrino mass. Within this framework we will then construct an entirely new bound from the IceCube-170922A event which takes into account (i) the possible neutrino mass hierarchies, (ii) the effect of cosmological redshift on e.g. the number density of cosmic neutrino background neutrinos, and, (iii) the non-degeneracy of neutrino mass and flavour eigenstates. This builds on work by Kelly and Machado (2018), where the authors placed new constraints on neutrinophilic and axion dark matter models. At low mediator masses, we find an improvement of an order of magnitude on current constraints from kaon decays. The constraint is complimentary (and slightly weaker) than current constraints from Big Bang Nucleosynthesis and the Cosmic Microwave Background. We explore how future higher energy events could improve this bound.


title: A census of symbiotic stars in the 2MASS, WISE and Gaia surveys

abstract: We present a new census of Galactic and extragalactic symbiotic stars (SySts). This compilation contains 323 known and 87 candidate SySts. Of the confirmed SySts, 257 are Galactic and 66 extragalactic. The spectral energy distributions (SEDs) of 348 sources have been constructed using 2MASS and AllWISE data. Regarding the Galactic SySts, 74% are S-types, 13% D and 3.5% D$^{\prime}$. S-types show an SED peak between 0.8 and 1.7 $\mu$m, whereas D-type show a peak at longer wavelengths between 2 and 4 $\mu$m. D$^{\prime}$-type, on the other hand, display a nearly flat profile. Gaia distances and effective temperatures are also presented. According to their Gaia distances, S-type are found to be members of both thin and thick Galactic disk populations, while S+IR- and D-types are mainly thin disk sources. Gaia temperatures show a reasonable agreement with the temperatures derived from SEDs within their uncertainties. A new census of the OVI $\lambda$6830 Raman-scattered line in SySts is also presented. From a sample of 298 SySts with available optical spectra, 55% are found to emit the line. No significant preference is found among the different types. The report of the OVI $\lambda$6830 Raman-scattered line in non-SySts is also discussed as well as the correlation between the Raman-scattered OVI line and X-ray emission. We conclude that the presence of the OVI Raman-scattered line still provides a strong criterion for identifying a source as a SySt.


title: Efficient description of experimental effects in amplitude analyses

abstract: Amplitude analysis is a powerful technique to study hadron decays. A significant complication in these analyses is the treatment of instrumental effects, such as background and selection efficiency variations, in the multidimensional kinematic phase space. This paper reviews conventional methods to estimate efficiency and background distributions and outlines the methods of density estimation using Gaussian processes and artificial neural networks. Such techniques see widespread use elsewhere, but have not gained popularity in use for amplitude analyses. Finally, novel applications of these models are proposed, to estimate background density in the signal region from the sidebands in multiple dimensions, and a more general method for model-assisted density estimation using artificial neural networks.


title: PVNet: A LRCN Architecture for Spatio-Temporal Photovoltaic PowerForecasting from Numerical Weather Prediction

abstract: Photovoltaic (PV) power generation has emerged as one of the lead renewable energy sources. Yet, its production is characterized by high uncertainty, being dependent on weather conditions like solar irradiance and temperature. Predicting PV production, even in the 24-hour forecast, remains a challenge and leads energy providers to left idling - often carbon emitting - plants. In this paper, we introduce a Long-Term Recurrent Convolutional Network using Numerical Weather Predictions (NWP) to predict, in turn, PV production in the 24-hour and 48-hour forecast horizons. This network architecture fully leverages both temporal and spatial weather data, sampled over the whole geographical area of interest. We train our model on an NWP dataset from the National Oceanic and Atmospheric Administration (NOAA) to predict spatially aggregated PV production in Germany. We compare its performance to the persistence model and state-of-the-art methods.


title: External Labeling Techniques: A Taxonomy and Survey

abstract: External labeling is frequently used for annotating features in graphical displays and visualizations, such as technical illustrations, anatomical drawings, or maps, with textual information. Such a labeling connects features within an illustration by thin leader lines with their labels, which are placed in the empty space surrounding the image. Over the last twenty years, a large body of literature in diverse areas of computer science has been published that investigates many different aspects, models, and algorithms for automatically placing external labels for a given set of features. This state-of-the-art report introduces a first unified taxonomy for categorizing the different results in the literature and then presents a comprehensive survey of the state of the art, a sketch of the most relevant algorithmic techniques for external labeling algorithms, as well as a list of open research challenges in this multidisciplinary research field.


title: COME TOGETHER: Multi-Agent Geometric Consensus (Gathering, Rendezvous, Clustering, Aggregation)

abstract: This report surveys results on distributed systems comprising mobile agents that are identical and anonymous, oblivious and interact solely by adjusting their motion according to the relative location of their neighbours. The agents are assumed capable of sensing the presence of other agents within a given sensing range and able to implement rules of motion based on full or partial information on the geometric constellation of their neighbouring agents. Eight different problems that cover assumptions of finite vs infinite sensing range, direction and distance vs direction only sensing and discrete vs continuous motion, are analyzed in the context of geometric consensus, clustering or gathering tasks.


title: A Sieve-SMM Estimator for Dynamic Models

abstract: This paper proposes a Sieve Simulated Method of Moments (Sieve-SMM) estimator for the parameters and the distribution of the shocks in nonlinear dynamic models where the likelihood and the moments are not tractable. An important concern with SMM, which matches sample with simulated moments, is that a parametric distribution is required. However, economic quantities that depend on this distribution, such as welfare and asset-prices, can be sensitive to misspecification. The Sieve-SMM estimator addresses this issue by flexibly approximating the distribution of the shocks with a Gaussian and tails mixture sieve. The asymptotic framework provides consistency, rate of convergence and asymptotic normality results, extending existing results to a new framework with more general dynamics and latent variables. An application to asset pricing in a production economy shows a large decline in the estimates of relative risk-aversion, highlighting the empirical relevance of misspecification bias.


title: ParBlockchain: Leveraging Transaction Parallelism in Permissioned Blockchain Systems

abstract: Many existing blockchains do not adequately address all the characteristics of distributed system applications and suffer from serious architectural limitations resulting in performance and confidentiality issues. While recent permissioned blockchain systems, have tried to overcome these limitations, their focus has mainly been on workloads with no-contention, i.e., no conflicting transactions. In this paper, we introduce OXII, a new paradigm for permissioned blockchains to support distributed applications that execute concurrently. OXII is designed for workloads with (different degrees of) contention. We then present ParBlockchain, a permissioned blockchain designed specifically in the OXII paradigm. The evaluation of ParBlockchain using a series of benchmarks reveals that its performance in workloads with any degree of contention is better than the state of the art permissioned blockchain systems.


title: Randomized resolvent analysis

abstract: Performing global resolvent analysis for high-Reynolds-number turbulent flow calls for the handling of a large discrete operator. Even though such large operator is required in the analysis, most applications of resolvent analysis extracts only a few dominant resolvent response and forcing modes. Here, we consider the use of randomized numerical linear algebra to reduce the dimension of the resolvent operator for achieving computational speed up and memory saving compared to the standard resolvent analysis. To accomplish this goal, we utilize sketching of the linear operator with random test matrices with a Gaussian distribution and with insights from the base flow incorporated to perform singular value decomposition on a low-rank matrix holding dominant characteristics of the full resolvent operator. The strength of the randomized resolvent analysis is demonstrated on a turbulent separated flow over an airfoil. This randomized approach clears the path towards tackling resolvent analysis for higher-Reynolds number bi- and tri-global base flows.


title: Autonomous Tissue Manipulation via Surgical Robot Using Learning Based Model Predictive Control

abstract: Tissue manipulation is a frequently used fundamental subtask of any surgical procedures, and in some cases it may require the involvement of a surgeon's assistant. The complex dynamics of soft tissue as an unstructured environment is one of the main challenges in any attempt to automate the manipulation of it via a surgical robotic system. Two AI learning based model predictive control algorithms using vision strategies are proposed and studied: (1) reinforcement learning and (2) learning from demonstration. Comparison of the performance of these AI algorithms in a simulation setting indicated that the learning from demonstration algorithm can boost the learning policy by initializing the predicted dynamics with given demonstrations. Furthermore, the learning from demonstration algorithm is implemented on a Raven IV surgical robotic system and successfully demonstrated feasibility of the proposed algorithm using an experimental approach. This study is part of a profound vision in which the role of a surgeon will be redefined as a pure decision maker whereas the vast majority of the manipulation will be conducted autonomously by a surgical robotic system. A supplementary video can be found at: http://bionics.seas.ucla.edu/research/surgeryproject17.html


title: Efficient representation of Gaussian states for multi-mode non-Gaussian quantum state engineering via subtraction of arbitrary number of photons

abstract: We introduce a complete description of a multi-mode bosonic quantum state in the coherent-state basis (which in this work is denoted as "$K$" function ), which---up to a phase---is the square root of the well-known Husimi "$Q$" representation. We express the $K$ function of any $N$-mode Gaussian state as a function of its covariance matrix and displacement vector, and also that of a general continuous-variable cluster state in terms of the modal squeezing and graph topology of the cluster. This formalism lets us characterize the non Gaussian state left over when one measures a subset of modes of a Gaussian state using photon number resolving detection, the fidelity of the obtained non-Gaussian state with any target state, and the associated heralding probability, all analytically. We show that this probability can be expressed as a Hafnian, re-interpreting the output state of a circuit claimed to demonstrate quantum supremacy termed Gaussian boson sampling. As an example-application of our formalism, we propose a method to prepare a two-mode coherent-cat-basis Bell state with fidelity close to unity and success probability that is fundamentally higher than that of a well-known scheme that splits an approximate single-mode cat state---obtained by photon number subtraction on a squeezed vacuum mode---on a balanced beam splitter. This formalism could enable exploration of efficient generation of cat-basis entangled states, which are known to be useful for quantum error correction against photon loss.


title: (Near) Optimal Adaptivity Gaps for Stochastic Multi-Value Probing

abstract: Consider a kidney-exchange application where we want to find a max-matching in a random graph. To find whether an edge $e$ exists, we need to perform an expensive test, in which case the edge $e$ appears independently with a \emph{known} probability $p_e$. Given a budget on the total cost of the tests, our goal is to find a testing strategy that maximizes the expected maximum matching size. The above application is an example of the stochastic probing problem. In general the optimal stochastic probing strategy is difficult to find because it is \emph{adaptive}---decides on the next edge to probe based on the outcomes of the probed edges. An alternate approach is to show the \emph{adaptivity gap} is small, i.e., the best \emph{non-adaptive} strategy always has a value close to the best adaptive strategy. This allows us to focus on designing non-adaptive strategies that are much simpler. Previous works, however, have focused on Bernoulli random variables that can only capture whether an edge appears or not. In this work we introduce a multi-value stochastic probing problem, which can also model situations where the weight of an edge has a probability distribution over multiple values. Our main technical contribution is to obtain (near) optimal bounds for the (worst-case) adaptivity gaps for multi-value stochastic probing over prefix-closed constraints. For a monotone submodular function, we show the adaptivity gap is at most $2$ and provide a matching lower bound. For a weighted rank function of a $k$-extendible system (a generalization of intersection of $k$ matroids), we show the adaptivity gap is between $O(k\log k)$ and $k$. None of these results were known even in the Bernoulli case where both our upper and lower bounds also apply, thereby resolving an open question of Gupta et al.


title: Modeling and Analysis of Non-unique Behaviors in Multiple Frictional Impacts

abstract: Many fundamental challenges in robotics, based in manipulation or locomotion, require making and breaking contact with the environment. To represent the complexity of frictional contact events, impulsive impact models are especially popular, as they often lead to mathematically and computationally tractable approaches. However, when two or more impacts occur simultaneously, the precise sequencing of impact forces is generally unknown, leading to the potential for multiple possible outcomes. This simultaneity is far from pathological, and occurs in many common robotics applications. In this work, we propose an approach for resolving simultaneous frictional impacts, represented as a differential inclusion. Solutions to our model, an extension to multiple contacts of Routh's method, naturally capture the set of potential post-impact velocities.We prove that solutions to the presented model must terminate. This is, to the best of our knowledge, the first such guarantee for set-valued outcomes to simultaneous frictional impacts.


title: Auxiliary Field Quantum Monte Carlo for Multiband Hubbard Models: Controlling the Sign and Phase Problems to Capture Hund's Physics

abstract: In the study of strongly-correlated, many-electron systems, the Hubbard Kanamori (HK) model has emerged as one of the prototypes for transition metal oxide physics. The model is multi-band in nature and contains Hund's coupling terms, which have pronounced effects on metal-insulator transitions, high-temperature superconductivity, and other physical properties. In the following, we present a complete theoretical framework for treating the HK model using the ground state Auxiliary Field Quantum Monte Carlo (AFQMC) method and analyze its performance on few-band models whose parameters approximate those observed in ruthenate, rhodates, and other materials exhibiting Hund's physics. Unlike previous studies, the constrained path and phaseless approximations are used to respectively control the sign and phase problems, which enables high accuracy modeling of the HK model's ground state properties within parameter regimes of experimental interest. We demonstrate that, after careful consideration of the Hubbard-Stratonovich transformations and trial wave functions employed, relative errors in the energy of less than 1% can routinely be achieved for moderate to large values of the Hund's coupling constant. Crucially, our methodology also accurately predicts magnetic ordering and phase transitions. The results presented open the door to more predictive modeling of Hund's physics within a wide range of strongly-correlated materials using AFQMC.


title: Optical and Opto-Mechanical Design of a Novel "Macro" Image Slicer for the MIRADAS Instrument

abstract: We present the innovative macro-slicer optical and opto-mechanical designs for the third-generation Mid-resolution InfraReD Astronomical Spectrograph (MIRADAS) instrument for the 10.4m Gran Telescopio Canarias (GTC) in the 1-2.5 $\mu$m bandpass. MIRADAS uses up to 12 cryogenic, fully steerable probes to select simultaneous targets in a 5 arcminute field of view. The spectrograph module is a cross-dispersed echelle spectrograph. The macro-slicer is effectively a stack of six advanced image slicer Integral Field Units (IFUs) such as FRIDA or FISICA, and like other IFUs designed and built at the University of Florida by our group, uses a `bolt-and-go' approach to minimize the difficulty in alignment and maximize robustness. Like other advanced image slicer IFUs, there are three sets of mirrors that work together to geometrically rearrange the loosely packed inputs from the probe arms into a tightly packed pseudo-slit. The macro-slicer also passively keeps the spectral resolution of MIRADAS fixed at $R>20,000$ in seeing from 1.2 arcseconds down to 0.4 arcseconds, (typical observing conditions at GTC).


title: Quad-rotor Flight Simulation in Realistic Atmospheric Conditions

abstract: In trajectory planning and control design for unmanned air vehicles, highly simplified models are typically used to represent the vehicle dynamics and the operating environment. The goal of this work is to perform real-time, but realistic flight simulations and trajectory planning for quad-copters in low altitude (<500m) atmospheric conditions. The aerodynamic model for rotor performance is adapted from blade element momentum theory and validated against experimental data. Large-eddy simulations of the atmospheric boundary layer are used to accurately represent the operating environment of unmanned air vehicles. A reduced-order version of the atmospheric boundary layer data as well as the popular Dryden model are used to assess the impact of accuracy of the wind field model on the predicted vehicle performance and trajectory. The wind model, aerodynamics and control modules are integrated into a six-degree-of-freedom flight simulation environment with a fully nonlinear flight controller. Simulations are performed for two representative flight paths, namely, straight and circular paths. Results for different wind models are compared and the impact of simplifying assumptions in representing rotor aerodynamics is discussed. The simulation framework and codes are open-sourced for use by the community.


title: TrackNet: Simultaneous Object Detection and Tracking and Its Application in Traffic Video Analysis

abstract: Object detection and object tracking are usually treated as two separate processes. Significant progress has been made for object detection in 2D images using deep learning networks. The usual tracking-by-detection pipeline for object tracking requires that the object is successfully detected in the first frame and all subsequent frames, and tracking is done by associating detection results. Performing object detection and object tracking through a single network remains a challenging open question. We propose a novel network structure named trackNet that can directly detect a 3D tube enclosing a moving object in a video segment by extending the faster R-CNN framework. A Tube Proposal Network (TPN) inside the trackNet is proposed to predict the objectness of each candidate tube and location parameters specifying the bounding tube. The proposed framework is applicable for detecting and tracking any object and in this paper, we focus on its application for traffic video analysis. The proposed model is trained and tested on UA-DETRAC, a large traffic video dataset available for multi-vehicle detection and tracking, and obtained very promising results.


title: Circles in self dual symmetric R-spaces

abstract: Self dual symmetric R-spaces have special curves, called circles, introduced by Burstall, Donaldson, Pedit and Pinkall in 2011, whose definition does not involve the choice of any Riemannian metric. We characterize the elements of the big transformation group G of a self dual symmetric R-space M as those diffeomorphisms of M sending circles in circles. Besides, although these curves belong to the realm of the invariants by G, we manage to describe them in Riemannian geometric terms: Given a circle c in M, there is a maximal compact subgroup K of G such that c, except for a projective reparametrization, is a diametrical geodesic in M (or equivalently, a diagonal geodesic in a maximal totally geodesic flat torus of M), provided that M carries the canonical symmetric K-invariant metric. We include examples for the complex quadric and the split standard or isotropic Grassmannians.


title: Hyperbolicity in Spherical Gravitational Collapse in a Horndeski Theory

abstract: We numerically study spherical gravitational collapse in shift symmetric Einstein dilaton Gauss Bonnet (EdGB) gravity. We find evidence that there are open sets of initial data for which the character of the system of equations changes from hyperbolic to elliptic type in a compact region of the spacetime. In these cases evolution of the system, treated as a hyperbolic initial boundary value problem, leads to the equations of motion becoming ill-posed when the elliptic region forms. No singularities or discontinuities are encountered on the corresponding effective "Cauchy horizon". Therefore it is conceivable that a well-posed formulation of EdGB gravity (at least within spherical symmetry) may be possible if the equations are appropriately treated as mixed-type.


title: Balleans, hyperballeans and ideals

abstract: A ballean $\mathcal{B}$ (or a coarse structure) on a set $X$ is a family of subsets of $X$ called balls (or entourages of the diagonal in $X\times X$) defined in such a way that $\mathcal{B}$ can be considered as the asymptotic counterpart of a uniform topological space. The aim of this paper is to study two concrete balleans defined by the ideals in the Boolean algebra of all subsets of $X$ and their hyperballeans, with particular emphasis on their connectedness structure, more specifically the number of their connected components.


title: Recursive projection-aggregation decoding of Reed-Muller codes

abstract: We propose a new class of efficient decoding algorithms for Reed-Muller (RM) codes over binary-input memoryless channels. The algorithms are based on projecting the code on its cosets, recursively decoding the projected codes (which are lower-order RM codes), and aggregating the reconstructions (e.g., using majority votes). We further provide extensions of the algorithms using list-decoding. We run our algorithm for AWGN channels and Binary Symmetric Channels at the short code length ($\le 1024$) regime for a wide range of code rates. Simulation results show that in both low code rate and high code rate regimes, the new algorithm outperforms the widely used decoder for polar codes (SCL+CRC) with the same parameters. The performance of the new algorithm for RM codes in those regimes is in fact close to that of the maximal likelihood decoder. Finally, the new decoder naturally allows for parallel implementations.


title: Strong convergence rates for Markovian representations of fractional processes

abstract: Many fractional processes can be represented as an integral over a family of Ornstein-Uhlenbeck processes. This representation naturally lends itself to numerical discretizations, which are shown in this paper to have strong convergence rates of arbitrarily high polynomial order. This explains the potential, but also some limitations of such representations as the basis of Monte Carlo schemes for fractional volatility models such as the rough Bergomi model.


title: Hyperballeans of groups

abstract: In this paper we define some ballean structure on the power set of a group and, in particular, we study the subballean with support the lattice of all its subgroups. If $G$ is a group, we denote by $L(G)$ the family of all subgroups of $G$. For two groups $G$ and $H$, we relate their algebraic structure via the ballean structure of $L(G)$ and $L(H)$.


title: Weyl-Schr\"odinger representations of Heisenberg groups in infinite dimensions

abstract: A complexified Heisenberg matrix group $\mathrm{H}_\mathbb{C}$ with entries from an infinite-dimensional Hilbert space $H$ is investigated. The Weyl--Schr\"odinger type irreducible representations of $\mathrm{H}_\mathbb{C}$ on the space $L^2_\chi$ of square-integrable scalar functions is described. The integrability is understood under the invariant probability measure $\chi$ which satisfies an abstract Kolmogorov consistency conditions over the infinite-dimensional unitary group $U(\infty)$ irreducible acted on $H$. The space $L^2_\chi$ is generated by Schur polynomials in variables on Paley--Wiener maps over $U(\infty)$. Therewith, the Fourier-image of $L^2_\chi$ coincides with a space of Hilbert--Schmidt entire analytic functions on $H$ generated by suitable Fock space. Applications to linear and nonlinear heat equations over the group $\mathrm{H}_\mathbb{C}$ are considered.


title: Optimized Compilation of Aggregated Instructions for Realistic Quantum Computers

abstract: Recent developments in engineering and algorithms have made real-world applications in quantum computing possible in the near future. Existing quantum programming languages and compilers use a quantum assembly language composed of 1- and 2-qubit (quantum bit) gates. Quantum compiler frameworks translate this quantum assembly to electric signals (called control pulses) that implement the specified computation on specific physical devices. However, there is a mismatch between the operations defined by the 1- and 2-qubit logical ISA and their underlying physical implementation, so the current practice of directly translating logical instructions into control pulses results in inefficient, high-latency programs. To address this inefficiency, we propose a universal quantum compilation methodology that aggregates multiple logical operations into larger units that manipulate up to 10 qubits at a time. Our methodology then optimizes these aggregates by (1) finding commutative intermediate operations that result in more efficient schedules and (2) creating custom control pulses optimized for the aggregate (instead of individual 1- and 2-qubit operations). Compared to the standard gate-based compilation, the proposed approach realizes a deeper vertical integration of high-level quantum software and low-level, physical quantum hardware. We evaluate our approach on important near-term quantum applications on simulations of superconducting quantum architectures. Our proposed approach provides a mean speedup of $5\times$, with a maximum of $10\times$. Because latency directly affects the feasibility of quantum computation, our results not only improve performance but also have the potential to enable quantum computation sooner than otherwise possible.


title: HAHE: Hierarchical Attentive Heterogeneous Information Network Embedding

abstract: Heterogeneous information network (HIN) embedding has recently attracted much attention due to its effectiveness in dealing with the complex heterogeneous data. Meta path, which connects different object types with various semantic meanings, is widely used by existing HIN embedding works. However, several challenges have not been addressed so far. First, different meta paths convey different semantic meanings, while existing works assume that all nodes share same weights for meta paths and ignore the personalized preferences of different nodes on different meta paths. Second, given a meta path, nodes in HIN are connected by path instances while existing works fail to fully explore the differences between path instances that reflect nodes' preferences in the semantic space. rTo tackle the above challenges, we propose aHierarchical Attentive Heterogeneous information network Embedding (HAHE) model to capture the personalized preferences on meta paths and path instances in each semantic space. As path instances are based on a particular meta path, a hierarchical attention mechanism is naturally utilized to model the personalized preference on meta paths and path instances. Extensive experiments on several real-world datasets show that our proposed \model model significantly outperforms the state-of-the-art methods in terms of various data mining tasks.


title: Electron energy partition across interplanetary shocks: I. Methodology and Data Product

abstract: Analysis of 15314 electron velocity distribution functions (VDFs) within $\pm$2 hours of 52 interplanetary (IP) shocks observed by the \emph{Wind} spacecraft near 1 AU are introduced. The electron VDFs are fit to the sum of three model functions for the cold dense core, hot tenuous halo, and field-aligned beam/strahl component. The best results were found by modeling the core as either a bi-kappa or a symmetric (or asymmetric) bi-self-similar velocity distribution function, while both the halo and beam/strahl components were best fit to bi-kappa velocity distribution function. This is the first statistical study to show that the core electron distribution is better fit to a self-similar velocity distribution function than a bi-Maxwellian under all conditions. The self-similar distribution deviation from a Maxwellian is a measure of inelasticity in particle scattering from waves and/or turbulence. The range of values defined by the lower and upper quartiles for the kappa exponents are $\kappa{\scriptstyle_{ec}}$ $\sim$ 5.40--10.2 for the core, $\kappa{\scriptstyle_{eh}}$ $\sim$ 3.58--5.34 for the halo, and $\kappa{\scriptstyle_{eb}}$ $\sim$ 3.40--5.16 for the beam/strahl. The lower-to-upper quartile range of symmetric bi-self-similar core exponents are $s{\scriptstyle_{ec}}$ $\sim$ 2.00--2.04, and asymmetric bi-self-similar core exponents are $p{\scriptstyle_{ec}}$ $\sim$ 2.20--4.00 for the parallel exponent, and $q{\scriptstyle_{ec}}$ $\sim$ 2.00--2.46 for the perpendicular exponent. The nuanced details of the fit procedure and description of resulting data product are also presented. The statistics and detailed analysis of the results are presented in Paper II and Paper III of this three-part study.


title: Faster way to agony: Discovering hierarchies in directed graphs

abstract: Many real-world phenomena exhibit strong hierarchical structure. Consequently, in many real-world directed social networks vertices do not play equal role. Instead, vertices form a hierarchy such that the edges appear mainly from upper levels to lower levels. Discovering hierarchies from such graphs is a challenging problem that has gained attention. Formally, given a directed graph, we want to partition vertices into levels such that ideally there are only edges from upper levels to lower levels. From computational point of view, the ideal case is when the underlying directed graph is acyclic. In such case, we can partition the vertices into a hierarchy such that there are only edges from upper levels to lower edges. In practice, graphs are rarely acyclic, hence we need to penalize the edges that violate the hierarchy. One practical approach is agony, where each violating edge is penalized based on the severity of the violation. The fastest algorithm for computing agony requires $O(nm^2)$ time. In the paper we present an algorithm for computing agony that has better theoretical bound, namely $O(m^2)$. We also show that in practice the obtained bound is pessimistic and that we can use our algorithm to compute agony for large datasets. Moreover, our algorithm can be used as any-time algorithm.


title: Accurate and robust segmentation of neuroanatomy in T1-weighted MRI by combining spatial priors with deep convolutional neural networks

abstract: Neuroanatomical segmentation in magnetic resonance imaging (MRI) of the brain is a prerequisite for volume, thickness and shape measurements. This work introduces a new highly accurate and versatile method based on 3D convolutional neural networks for the automatic segmentation of neuroanatomy in T1-weighted MRI. In combination with a deep 3D fully convolutional architecture, efficient linear registration-derived spatial priors are used to incorporate additional spatial context into the network. An aggressive data augmentation scheme using random elastic deformations is also used to regularize the networks, allowing for excellent performance even in cases where only limited labelled training data are available. Applied to hippocampus segmentation in an elderly population (mean Dice coefficient = 92.1%) and sub-cortical segmentation in a healthy adult population (mean Dice coefficient = 89.5%), we demonstrate new state-of-the-art accuracies and a high robustness to outliers with the same architecture. Further validation on a multi-structure segmentation task in a scan-rescan dataset demonstrates accuracy (mean Dice coefficient = 86.6%) similar to the scan-rescan reliability of expert manual segmentations (mean Dice coefficient = 86.9%), and improved reliability compared to both expert manual segmentations and automated segmentations using FIRST. Furthermore, our method maintains a highly competitive runtime performance (e.g. requiring only 10 seconds for left/right hippocampal segmentation in 1x1x1 MNI stereotaxic space), orders of magnitude faster than conventional multi-atlas segmentation methods.


title: Asymptotics of quasi-classical localized states in 2D system of charged hard-core bosons

abstract: The continuous quasi-classical two-sublattice approximation is constructed for the 2D system of charged hard-core bosons to explore metastable inhomogeneous states analogous to inhomogeneous localized excitations in magnetic systems. The types of localized excitations are determined by asymptotic analysis and compared with numerical results. Depending on the homogeneous ground state, the excitations are the ferro and antiferro type vortices, the skyrmion-like topological excitations or linear domain walls.


title: What is the dimension of your binary data?

abstract: Many 0/1 datasets have a very large number of variables; on the other hand, they are sparse and the dependency structure of the variables is simpler than the number of variables would suggest. Defining the effective dimensionality of such a dataset is a nontrivial problem. We consider the problem of defining a robust measure of dimension for 0/1 datasets, and show that the basic idea of fractal dimension can be adapted for binary data. However, as such the fractal dimension is difficult to interpret. Hence we introduce the concept of normalized fractal dimension. For a dataset $D$, its normalized fractal dimension is the number of columns in a dataset $D'$ with independent columns and having the same (unnormalized) fractal dimension as $D$. The normalized fractal dimension measures the degree of dependency structure of the data. We study the properties of the normalized fractal dimension and discuss its computation. We give empirical results on the normalized fractal dimension, comparing it against baseline measures such as PCA. We also study the relationship of the dimension of the whole dataset and the dimensions of subgroups formed by clustering. The results indicate interesting differences between and within datasets.


title: Numerical analysis of a mechanotransduction dynamical model reveals homoclinic bifurcations of extracellular matrix mediated oscillations of the mesenchymal stem cell fate

abstract: We perform one and two-parameter numerical bifurcation analysis of a mechanotransduction model approximating the dynamics of mesenchymal stem cell differentiation into neurons, adipocytes, myocytes and osteoblasts. For our analysis, we use as bifurcation parameters the stiffness of the extracellular matrix and parameters linked with the positive feedback mechanisms that up-regulate the production of the YAP/TAZ transcriptional regulators (TRs) and the cell adhesion area. Our analysis reveals a rich nonlinear behaviour of the cell differentiation including regimes of hysteresis and multistability, stable oscillations of the effective adhesion area, the YAP/TAZ TRs and the PPAR$\gamma$ receptors associated with the adipogenic fate, as well as homoclinic bifurcations that interrupt relatively high-amplitude oscillations abruptly. The two-parameter bifurcation analysis of the Andronov-Hopf points that give birth to the oscillating patterns predicts their existence for soft extracellular substrates ($<1kPa$), a regime that favours the neurogenic and the adipogenic cell fate. Furthermore, in these regimes, the analysis reveals the presence of homoclinic bifurcations that result in the sudden loss of the stable oscillations of the cell-substrate adhesion towards weaker adhesion and high expression levels of the gene encoding Tubulin beta-3 chain, thus favouring the phase transition from the adipogenic to the neurogenic fate.


title: Bootstrapped Coordinate Search for Multidimensional Scaling

abstract: In this work, a unified framework for gradient-free Multidimensional Scaling (MDS) based on Coordinate Search (CS) is proposed. This family of algorithms is an instance of General Pattern Search (GPS) methods which avoid the explicit computation of derivatives but instead evaluate the objective function while searching on coordinate steps of the embedding space. The backbone element of CSMDS framework is the corresponding probability matrix that correspond to how likely is each corresponding coordinate to be evaluated. We propose a Bootstrapped instance of CSMDS (BS CSMDS) which enhances the probability of the direction that decreases the most the objective function while also reducing the corresponding probability of all the other coordinates. BS CSMDS manages to avoid unnecessary function evaluations and result to significant speedup over other CSMDS alternatives while also obtaining the same error rate. Experiments on both synthetic and real data reveal that BS CSMDS performs consistently better than other CSMDS alternatives under various experimental setups.


title: Discovering Nested Communities

abstract: Finding communities in graphs is one of the most well-studied problems in data mining and social-network analysis. In many real applications, the underlying graph does not have a clear community structure. In those cases, selecting a single community turns out to be a fairly ill-posed problem, as the optimization criterion has to make a difficult choice between selecting a tight but small community or a more inclusive but sparser community. In order to avoid the problem of selecting only a single community we propose discovering a sequence of nested communities. More formally, given a graph and a starting set, our goal is to discover a sequence of communities all containing the starting set, and each community forming a denser subgraph than the next. Discovering an optimal sequence of communities is a complex optimization problem, and hence we divide it into two subproblems: 1) discover the optimal sequence for a fixed order of graph vertices, a subproblem that we can solve efficiently, and 2) find a good order. We employ a simple heuristic for discovering an order and we provide empirical and theoretical evidence that our order is good.


title: On the Complexity of Toric Ideals

abstract: We investigate the computational complexity of problems on toric ideals such as normal forms, Gr\"obner bases, and Graver bases. We show that all these problems are strongly NP-hard in the general case. Nonetheless, we can derive efficient algorithms by taking advantage of the sparsity pattern of the matrix. We describe this sparsity pattern with a graph, and study the parameterized complexity of toric ideals in terms of graph parameters such as treewidth and treedepth. In particular, we show that the normal form problem can be solved in parameter-tractable time in terms of the treedepth. An important application of this result is in multiway ideals arising in algebraic statistics. We also give a parameter-tractable membership test to the reduced Gr\"obner basis. This test leads to an efficient procedure for computing the reduced Gr\"obner basis. Similar results hold for Graver bases computation.


title: Conversion Measure of Faraday Rotation-Conversion with application to Fast Radio Bursts

abstract: Faraday Rotation-Conversion is the simultaneous rotation of all three Stokes polarization parameters $Q$, $U$, $V$ as an electromagnetic wave propagates through a magnetized plasma. In this regime the Faraday plasma screen is characterized by more than just a Rotation Measure. We define the Conversion Measure that characterizes the wavelength-dependent conversion between the linear and circular polarization. In a cold plasma, the conversion occurs at the localized regions along the wave's path, where the large-scale magnetic field is perpendicular to the propagation direction. The number of these regions along the line of sight through the screen, and their individual contributions to the conversion measure, can be inferred from the polarization measurements. We argue that the simultaneous measurement of wavelength-dependent linear and circular polarization might give an important insight into the magnetic-field geometry of the Faraday screen in FRB121102 and other repeating Fast Radio Bursts.


title: Stiefel manifolds and polygons

abstract: Polygons are compound geometric objects, but when trying to understand the expected behavior of a large collection of random polygons -- or even to formalize what a random polygon is -- it is convenient to interpret each polygon as a point in some parameter space, essentially trading the complexity of the object for the complexity of the space. In this paper I describe such an interpretation where the parameter space is a Stiefel manifold and show how to exploit the geometry of the Stiefel manifold both to generate random polygons and to morph one polygon into another.


title: Confusion matrices and rough set data analysis

abstract: A widespread approach in machine learning to evaluate the quality of a classifier is to cross -- classify predicted and actual decision classes in a confusion matrix, also called error matrix. A classification tool which does not assume distributional parameters but only information contained in the data is based on the rough set data model which assumes that knowledge is given only up to a certain granularity. Using this assumption and the technique of confusion matrices, we define various indices and classifiers based on rough confusion matrices.


title: Anapole Dark Matter via Vector Boson Fusion Processes at the LHC

abstract: Dark matter that is electrically neutral but couples to the electromagnetic current through higher-dimensional operators constitutes an interesting class of models. We investigate this class of models at the Large Hadron Collider, focusing on the anapole moment operator in an effective field theory (EFT) framework, and utilizing the vector boson fusion (VBF) topology. Assuming proton-proton collisions at $\sqrt{s} = 13$ TeV, we present the VBF anapole dark matter (ADM) cross sections and kinematic distributions as functions of the free parameters of the EFT, the cutoff scale $\Lambda$ and the ADM mass $m_{\chi}$. We find that the distinctive VBF topology of two forward jets and large dijet pseudorapidity gap is effective at reducing SM backgrounds, leading to a $5\sigma $ discovery reach for all kinematically allowed ADM masses with $\Lambda \le 1.65$ (1.15) TeV, assuming an integrated luminosity of 3000 (100) fb$^{-1}$.


title: Global Stability of a Class of Difference Equations on Solvable Lie Algebras

abstract: Motivated by the ubiquitous sampled-data setup in applied control, we examine the stability of a class of difference equations that arises by sampling a right- or left-invariant flow on a matrix Lie group. The map defining such a difference equation has three key properties that facilitate our analysis: 1) its power series expansion enjoys a type of strong convergence; 2) the origin is an equilibrium; 3) the algebraic ideals enumerated in the lower central series of the Lie algebra are dynamically invariant. We show that certain global stability properties are implied by stability of the Jacobian linearization of dynamics at the origin. In particular global asymptotic stability. If the Lie algebra is nilpotent, then the origin enjoys semiglobal exponential stability.


title: Higher symplectic capacities

abstract: We construct new families of symplectic capacities indexed by certain symmetric polynomials, defined using rational symplectic field theory. In particular, we introduce a sequence of capacities based on an L-infinity structure on linearized contact homology and rational curve counts with local tangency constraints. We prove various structural properties of these capacities and give some preliminary computations which show that they give state of the art symplectic embedding obstructions in basic examples.


title: Massive Argon Space Telescope (MAST): a concept of heavy time projection chamber for gamma-ray astronomy in the 100 MeV --- 1 TeV energy range

abstract: We explore the concept of liquid Argon time projection chamber (TPC) for gamma-ray astronomy in the 100 MeV --- 1 TeV energy range. We propose a basic layout for such a telescope called MAST. Using a last-generation rocket such as Falcon Heavy, it is possible to launch a detector with the effective area and the differential sensitivity about one order of magnitude better than the Fermi-LAT ones. At the same time, the MAST concept allows for an excellent angular resolution, 3-10 times better than the Fermi-LAT one depending on the energy, and good energy resolution ($\approx$ 20 % at 100 MeV and 6-10 % for the 10 GeV --- 1 TeV energy range). We show that such a telescope would be instrumental in a broad range of long-standing astrophysical problems.


title: Optimally Scheduling CNN Convolutions for Efficient Memory Access

abstract: Embedded inference engines for convolutional networks must be parsimonious in memory bandwidth and buffer sizing to meet power and cost constraints. We present an analytical memory bandwidth model for loop-nest optimization targeting architectures with application managed buffers. We applied this model to optimize the CNN convolution loop-nest. We show that our model is more accurate than previously published models. Using this model we can identify non-trivial dataflow schedules that result in lowest communication bandwidth given tight local buffering constraints. We show that optimal dataflow schedules are implementable in practice and that our model is accurate with respect to a real implementation; moreover, we introduce an accelerator architecture, named Hardware Convolution Block (HWC), which implements the optimal schedules, and we show it achieves up to 14x memory bandwidth reduction compared to a previously published accelerator with a similar memory interface, but implementing a non-optimal schedule.


title: Determination of pressure-viscosity relation of 2,2,4-trimethylhexane by all-atom molecular dynamics simulations

abstract: The Newtonian viscosity of 2,2,4-trimethylhexane at 293K is determined at pressures from 0.1MPa to 1000MPa. Non-equilibrium molecular dynamics simulations are performed using AIREBO-M, an all-atom potential for hydrocarbons especially parameterized for high pressures. The steady-state shear stress and viscosity are determined from simple shear simulations at rates between $10^7$ and $5\cdot 10^9\ \textrm{s}^{-1}$. At low pressures, simulation rates are low enough to reach the Newtonian regime. At high pressures, results are extrapolated to the Newtonian limit by fitting rate-dependent viscosities to Eyring theory. The resulting pressure dependent viscosity is typical of small molecules and fits to a common model are discussed.


title: Symbolic QED Pre-silicon Verification for Automotive Microcontroller Cores: Industrial Case Study

abstract: We present an industrial case study that demonstrates the practicality and effectiveness of Symbolic Quick Error Detection (Symbolic QED) in detecting logic design flaws (logic bugs) during pre-silicon verification. Our study focuses on several microcontroller core designs (~1,800 flip-flops, ~70,000 logic gates) that have been extensively verified using an industrial verification flow and used for various commercial automotive products. The results of our study are as follows: 1. Symbolic QED detected all logic bugs in the designs that were detected by the industrial verification flow (which includes various flavors of simulation-based verification and formal verification). 2. Symbolic QED detected additional logic bugs that were not recorded as detected by the industrial verification flow. (These additional bugs were also perhaps detected by the industrial verification flow.) 3. Symbolic QED enables significant design productivity improvements: (a) 8X improved (i.e., reduced) verification effort for a new design (8 person-weeks for Symbolic QED vs. 17 person-months using the industrial verification flow). (b) 60X improved verification effort for subsequent designs (2 person-days for Symbolic QED vs. 4-7 person-months using the industrial verification flow). (c) Quick bug detection (runtime of 20 seconds or less), together with short counterexamples (10 or fewer instructions) for quick debug, using Symbolic QED.


title: Existence and regularity of minimizers for nonlocal energy functionals

abstract: In this paper we consider minimizers for nonlocal energy functionals generalizing elastic energies that are connected with the theory of peridynamics \cite{Silling2000} or nonlocal diffusion models \cite{Rossi}. We derive nonlocal versions of the Euler-Lagrange equations under two sets of growth assumptions for the integrand. Existence of minimizers is shown for integrands with joint convexity (in the function and nonlocal gradient components). By using the convolution structure we show regularity of solutions for certain Euler-Lagrange equations. No growth assumptions are needed for the existence and regularity of minimizers results, in contrast with the classical theory.


title: A Two-Stream Siamese Neural Network for Vehicle Re-Identification by Using Non-Overlapping Cameras

abstract: We describe in this paper a Two-Stream Siamese Neural Network for vehicle re-identification. The proposed network is fed simultaneously with small coarse patches of the vehicle shape's, with 96 x 96 pixels, in one stream, and fine features extracted from license plate patches, easily readable by humans, with 96 x 48 pixels, in the other one. Then, we combined the strengths of both streams by merging the Siamese distance descriptors with a sequence of fully connected layers, as an attempt to tackle a major problem in the field, false alarms caused by a huge number of car design and models with nearly the same appearance or by similar license plate strings. In our experiments, with 2 hours of videos containing 2982 vehicles, extracted from two low-cost cameras in the same roadway, 546 ft away, we achieved a F-measure and accuracy of 92.6% and 98.7%, respectively. We show that the proposed network, available at https://github.com/icarofua/siamese-two-stream, outperforms other One-Stream architectures, even if they use higher resolution image features.


title: Asymptotic Theory for Clustered Samples

abstract: We provide a complete asymptotic distribution theory for clustered data with a large number of independent groups, generalizing the classic laws of large numbers, uniform laws, central limit theory, and clustered covariance matrix estimation. Our theory allows for clustered observations with heterogeneous and unbounded cluster sizes. Our conditions cleanly nest the classical results for i.n.i.d. observations, in the sense that our conditions specialize to the classical conditions under independent sampling. We use this theory to develop a full asymptotic distribution theory for estimation based on linear least-squares, 2SLS, nonlinear MLE, and nonlinear GMM.


title: Explicit Formula for Preimages of Relaxed One-Sided Lipschitz Mappings with Negative Lipschitz Constants

abstract: This paper addresses Lipschitzian stability issues that play an important role in both theoretical and numerical aspects of variational analysis, optimization, and their applications. We particularly concentrate on the so-called relaxed one-sided Lipschitz property of set-valued mappings with negative Lipschitz constants. This property has been much less investigated than more conventional Lipschitzian behavior while being well recognized in a variety of applications. Recent work has revealed that set-valued mappings satisfying the relaxed one-sided Lipschitz condition with negative Lipschitz constant possess a localization property that is stronger than uniform metric regularity. The present paper complements this fact by providing a characterization not only of one specific single point of a preimage, but of entire preimages of such mappings. Developing a geometric approach, we derive an explicit formula to calculate preimages of relaxed one-sided Lipschitz mappings between finite-dimensional spaces and obtain a further specification of this formula via extreme points of image sets.


title: Differentially Private Release of High-Dimensional Datasets using the Gaussian Copula

abstract: We propose a generic mechanism to efficiently release differentially private synthetic versions of high-dimensional datasets with high utility. The core technique in our mechanism is the use of copulas. Specifically, we use the Gaussian copula to define dependencies of attributes in the input dataset, whose rows are modelled as samples from an unknown multivariate distribution, and then sample synthetic records through this copula. Despite the inherently numerical nature of Gaussian correlations we construct a method that is applicable to both numerical and categorical attributes alike. Our mechanism is efficient in that it only takes time proportional to the square of the number of attributes in the dataset. We propose a differentially private way of constructing the Gaussian copula without compromising computational efficiency. Through experiments on three real-world datasets, we show that we can obtain highly accurate answers to the set of all one-way marginal, and two-and three-way positive conjunction queries, with 99\% of the query answers having absolute (fractional) error rates between 0.01 to 3\%. Furthermore, for a majority of two-way and three-way queries, we outperform independent noise addition through the well-known Laplace mechanism. In terms of computational time we demonstrate that our mechanism can output synthetic datasets in around 6 minutes 47 seconds on average with an input dataset of about 200 binary attributes and more than 32,000 rows, and about 2 hours 30 mins to execute a much larger dataset of about 700 binary attributes and more than 5 million rows. To further demonstrate scalability, we ran the mechanism on larger (artificial) datasets with 1,000 and 2,000 binary attributes (and 5 million rows) obtaining synthetic outputs in approximately 6 and 19 hours, respectively.


title: Parameter-Free Online Convex Optimization with Sub-Exponential Noise

abstract: We consider the problem of unconstrained online convex optimization (OCO) with sub-exponential noise, a strictly more general problem than the standard OCO. In this setting, the learner receives a subgradient of the loss functions corrupted by sub-exponential noise and strives to achieve optimal regret guarantee, without knowledge of the competitor norm, i.e., in a parameter-free way. Recently, Cutkosky and Boahen (COLT 2017) proved that, given unbounded subgradients, it is impossible to guarantee a sublinear regret due to an exponential penalty. This paper shows that it is possible to go around the lower bound by allowing the observed subgradients to be unbounded via stochastic noise. However, the presence of unbounded noise in unconstrained OCO is challenging; existing algorithms do not provide near-optimal regret bounds or fail to have a guarantee. So, we design a novel parameter-free OCO algorithm for Banach space, which we call BANCO, via a reduction to betting on noisy coins. We show that BANCO achieves the optimal regret rate in our problem. Finally, we show the application of our results to obtain a parameter-free locally private stochastic subgradient descent algorithm, and the connection to the law of iterated logarithms.


title: Measurement of the scintillation efficiency for nuclear recoils in liquid argon under electric fields up to 3 kV/cm

abstract: We present a measurement of scintillation efficiency for a few tens of keV nuclear recoils (NR) with a liquid argon time projection chamber under electric fields ranging from $0$ to $3\ {\rm kV/cm}$. The calibration data are taken with $^{252}{\rm Cf}$ radioactive source. Observed scintillation and electroluminescence spectra are simultaneously fit with spectra derived from Geant4-based Monte Carlo simulation and an NR model. The scintillation efficiency extracted from the fit is reported as a function of recoil energy and electric field. This result can be used for designing the detector and for the interpretation of experimental data in searching for scintillation and ionization signals induced by WIMP dark matter.


title: An analysis of a mathematical model describing the growth of a tumor treated with chemotherapy

abstract: We present a mathematical analysis of a mixed ODE-PDE model describing the spatial distribution and temporal evolution of tumor and normal cells within a tissue subject to the effects of a chemotherapeutic drug. The model assumes that the influx of chemotherapy is restricted to a limited region of the tissue, mimicking a blood vessel passing transversely. We provide results on the existence and uniqueness of the model solution and numerical simulations illustrating different model behaviors.


title: Environmental effects on halo abundance and weak lensing peak statistics toward large underdense regions

abstract: The cosmic microwave background (CMB) contains an anomalous cold spot with a surrounding hot ring, known as the Cold Spot. Inoue & Silk (2006) proposed that this feature could be explained by postulating a supervoid: if such a large underdense region exists, then the growth of matter perturbing around the spot might differ from the average value in the Universe and the differences might affect weak lensing analysis of peak statistics. To investigate environmental effects on halo number count and peak statistics, we used a publicly available ray-tracing simulation for a box size of 2250$h^{-1}$Mpc on a side (Takahashi et al. 2017). We found that the number counts for massive haloes toward the largest underdense region in the simulation decreases and the corresponding significance of the difference, based on a cosmic average, is $\geq3\sigma$. On the basis of the results of peak statistics analysis, the number of high peaks decreases with the decrement of massive haloes, but the number of low peaks increases with the lack of matter in the line of sight. The highest significance of the decrement in peak counts in large underdense regions is $5\sigma$ in the total signal-to-noise ratio. Our result implies that environmental effects on halo abundance and weak lensing peak statistic can be used to probe the presence and properties of supervoids.


title: Weak gravitational quantum effects in boson particles

abstract: We rewrite the Klein-Gordon (KG) equation in an arbitrary space-time transforming it into a generalized Schr\"odinger equation. Then we take the weak field limit and show that this equation has some differences with the traditional Schr\"odinger equation plus a gravitational field. Thus, this procedure shows that the Schr\"odinger equation derived in a covariant manner is different from the traditional one. With this new Schr\"odinger equation, we study the KG equation in a Newtonian space-time to describe the behavior of a scalar particle in an inertial system. We give some examples where it is possible to study the energy levels, effective potential and the wave function of the systems, these results contain the gravitational effects due to the curvature of space-time. We show that it is possible to verify experimentally these effects in a laboratory using non-inertial reference frames.


title: Optimal control of a thermistor problem with vanishing conductivity

abstract: An optimal control of a steady state thermistor problem is considered, where the convective boundary coefficient is taken as the control variable. A distinctive feature of this paper is that the problem is considered in arbitrary dimensions, and the electrical conductivity is allowed to vanish above a threshold temperature value. The existence of a steady state is proved, an objective functional is introduced, the existence of the optimal control is proved, and the optimality system is derived.


title: Learning to Prescribe Interventions for Tuberculosis Patients Using Digital Adherence Data

abstract: Digital Adherence Technologies (DATs) are an increasingly popular method for verifying patient adherence to many medications. We analyze data from one city served by 99DOTS, a phone-call-based DAT deployed for Tuberculosis (TB) treatment in India where nearly 3 million people are afflicted with the disease each year. The data contains nearly 17,000 patients and 2.1M dose records. We lay the groundwork for learning from this real-world data, including a method for avoiding the effects of unobserved interventions in training data used for machine learning. We then construct a deep learning model, demonstrate its interpretability, and show how it can be adapted and trained in different clinical scenarios to better target and improve patient care. In the real-time risk prediction setting our model could be used to proactively intervene with 21% more patients and before 76% more missed doses than current heuristic baselines. For outcome prediction, our model performs 40% better than baseline methods, allowing cities to target more resources to clinics with a heavier burden of patients at risk of failure. Finally, we present a case study demonstrating how our model can be trained in an end-to-end decision focused learning setting to achieve 15% better solution quality in an example decision problem faced by health workers.


title: Cyclic Symmetry on Complex Tori and Bagnera-De Franchis Manifolds

abstract: We describe the possible linear actions of a cyclic group $G = \mathbb{Z} /n$ on a complex torus, using the cyclotomic exact sequence for the group algebra $\mathbb{Z} [G]$. The main application is devoted to a structure theorem for Bagnera-De Franchis Manifolds, but we also give an application to hypergeometric integrals.


title: The covering number of the strong measure zero ideal can be above almost everything else

abstract: We show that certain type of tree forcings, including Sacks forcing, increases the covering of the strong measure zero ideal $\mathcal{SN}$. As a consequence, in Sacks model, such covering number is equal to the size of the continuum, which indicates that this covering number is consistently larger than any other classical cardinal invariant of the continuum. Even more, Sacks forcing can be used to force that $\mathrm{non}(\mathcal{SN})<\mathrm{cov}(\mathcal{SN})<\mathrm{cof}(\mathcal{SN})$, which is the first consistency result where more than two cardinal invariants associated with $\mathcal{SN}$ are pairwise different. Another consequence is that $\mathcal{SN}\subseteq s^0$ in ZFC where $s^0$ denotes the Marczewski's ideal.


title: Training on Synthetic Noise Improves Robustness to Natural Noise in Machine Translation

abstract: We consider the problem of making machine translation more robust to character-level variation at the source side, such as typos. Existing methods achieve greater coverage by applying subword models such as byte-pair encoding (BPE) and character-level encoders, but these methods are highly sensitive to spelling mistakes. We show how training on a mild amount of random synthetic noise can dramatically improve robustness to these variations, without diminishing performance on clean text. We focus on translation performance on natural noise, as captured by frequent corrections in Wikipedia edit logs, and show that robustness to such noise can be achieved using a balanced diet of simple synthetic noises at training time, without access to the natural noise data or distribution.


title: Proceedings Tenth International Workshop on Computing with Terms and Graphs

abstract: This volume contains a selection of the papers presented at TERMGRAPH 2018, the tenth edition of the international workshop on computing with terms and graphs. Graphs, and graph transformation systems, are used in many areas within Computer Science: to represent data structures and algorithms, to define computation models, as a general modelling tool to study complex systems, etc. Research in this area addresses a range of theoretical and practical issues, including the modelling of first- and higher-order term rewriting by (acyclic or cyclic) graph rewriting, graphical frameworks such as interaction nets and sharing graphs (optimal reduction), rewrite calculi for the analysis of functional programs, graph reduction implementations of programming languages, graphical calculi modelling concurrent and mobile computations, object-oriented systems, graphs as a model of biological or chemical systems, and automated reasoning and symbolic computation systems working on shared structures. Previous editions of TERMGRAPH took place in Barcelona (2002), Rome (2004), Vienna (2006), Braga (2007), York (2009), Saarbrucken (2011), Rome (2013), Vienna (2014) and Eindhoven (2016). TERMGRAPH 2018 is affiliated with FSCD, which is part of FLOC.


title: Observing monomer - dimer transitions of neurotensin receptors 1 in single SMALPs by homoFRET and in an ABELtrap

abstract: G protein-coupled receptors (GPCRs) are a large superfamily of membrane proteins that are activated by extracellular small molecules or photons. Neurotensin receptor 1 (NTSR1) is a GPCR that is activated by neurotensin, i.e. a 13 amino acid peptide. Binding of neurotensin induces conformational changes in the receptor that trigger the intracellular signaling processes. While recent single-molecule studies have reported a dynamic monomer - dimer equilibrium of NTSR1 in vitro, a biophysical characterization of the oligomerization status of NTSR1 in living mammalian cells is complicated. Here we report on the oligomerization state of the human NTSR1 tagged with mRuby3 by dissolving the plasma membranes of living HEK293T cells into 10 nm-sized soluble lipid nanoparticles by addition of styrene-maleic acid copolymers (SMALPs). Single SMALPs were analyzed one after another in solution by multi-parameter single molecule spectroscopy including brightness, fluorescence lifetime and anisotropy for homoFRET. Brightness analysis was improved using single SMALP detection in a confocal ABELtrap for extended observation times in solution. A bimodal brightness distribution indicated a significant fraction of dimeric NTSR1 in SMALPs or in the plasma membrane, respectively, before addition of neurotensin.


title: Generalized solutions to the Dirichlet problem of translating mean curvature equations

abstract: In this paper we study the Dirichlet problem of translating mean curvature equations over domains in Riemannian manifolds with dimension $n$. Imitating the generalized solution theory of Miranda-Giusti, we define a new conformal area functional and a generalized solution to this Dirichlet problem. The existence of generalized solutions to this problem on bounded Lipschitz domains is established. If the domain is mean convex and bounded with $C^2$ boundary, its closure does not contain any closed minimal hypersurface except a singular set with its Hausdorff dimension at most $n-7$ and the boundary data is continuous, the generalized solution is the desirable classical smooth solution. The non-minimal condition could not be removed in general.


title: Bitcoin Selfish Mining and Dyck Words

abstract: We give a straightforward proof for the formula giving the long-term apparent hashrate of the Selfish Mining strategy in Bitcoin using only elementary probabilities and combinatorics, and more precisely, Dyck words. There is no need to compute stationary probabilities on Markov chain nor stopping times for Poisson processes as it was previously done. We consider also several other block withholding strategies.


title: Perturbative GAN: GAN with Perturbation Layers

abstract: Perturbative GAN, which replaces convolution layers of existing convolutional GANs (DCGAN, WGAN-GP, BIGGAN, etc.) with perturbation layers that adds a fixed noise mask, is proposed. Compared with the convolu-tional GANs, the number of parameters to be trained is smaller, the convergence of training is faster, the incep-tion score of generated images is higher, and the overall training cost is reduced. Algorithmic generation of the noise masks is also proposed, with which the training, as well as the generation, can be boosted with hardware acceleration. Perturbative GAN is evaluated using con-ventional datasets (CIFAR10, LSUN, ImageNet), both in the cases when a perturbation layer is adopted only for Generators and when it is introduced to both Generator and Discriminator.


title: Model for a Noise Matched Phased Array Feed

abstract: We present a model for a Noise Matched Phased Array Feed (PAF) system and compare model predictions with the measurement results. The PAF system consists of an array feed, a receiver, a beamformer and a parabolic reflector. The novel aspect of our model is the characterization of the {\em PAF system} by a single matrix. This characteristic matrix is constructed from the open-circuit voltage covariance at the output of the PAF due to signal from the observing source, ground spillover noise, sky background noise and (low-noise) amplifier (LNA) noise. The best signal-to-noise ratio on the source achievable with the PAF system will be the maximum eigenvalue of the characteristic matrix. The voltage covariance due to signal and spillover noise are derived by applying the Lorentz reciprocity theorem. The receiver noise covariance and noise temperature are obtained in terms of Lange invariants such that they are suitable for noise matching the array feed. The model predictions are compared with the measured performance of a 1.4 GHz, 19-element, dual-polarized PAF on the Robert C. Byrd Green Bank Telescope. We show that the model predictions, obtained with an additional noise contribution due to the measured losses ahead of the low-noise amplifier, compare well with the measured ratio of system temperature to aperture efficiency as a function of frequency and as a function of offset from the boresight. Further, our modeling indicates that the bandwidth over which this ratio is optimum can be improved by a factor of at least two by noise matching the PAF with the LNA.


title: A Model for Phased Array Feed

abstract: In this report we present a model for phased array feed (PAF) and compare the model predictions with measurements. A theory for loss-less PAF is presented first. To develop the theory we ask the question -- what is the best $T_{sys}/\eta_{ap}$ that can be achieved when a PAF is used on a telescope to observe a source at an angle $\theta_s, \phi_s$ from the boresight direction ? We show that a characteristic matrix for the {\em system} (i.e. PAF+telescope+receiver) can be constructed starting from the signal-to-noise ratio of the observations and the best $T_{sys}/\eta_{ap}$ can be obtained from the maximum eigenvalue of the characteristic matrix. For constructing the characteristic matrix, we derive the open-circuit voltage at the output of the antenna elements in the PAF due to (a) radiation from source, (b) radiation from ground (spillover), (c) radiation from sky background and (d) noise due to the receiver. The characteristic matrix is then obtained from the correlation matrices of these voltages. We then describe a modeling program developed to implement the theory presented here. Finally the model predictions are compared with results from test observations made toward Virgo A with a prototype PAF (Kite array) on the GBT (Roshi et al. 2015).


title: Phased Array Feed Model Equations corresponding to two definitions of embedded beam pattern

abstract: In this report, we present the phased array feed (PAF) model equations for two definitions of embedded beam patterns. In Roshi \& Fisher (2016), we presented the PAF model by defining the embedded beam pattern as the beam pattern due to a 1 V excitation to one port and all other ports short circuited. This embedded beam pattern is referred to as voltage-embedded-beam (VEB). The embedded beam pattern can also be defined as the beam pattern due to a 1 A excitation to one port and all other ports open circuited. This definition is usually used in engineering literature and we refer to the pattern as current-embedded-beam (CEB). Here we derive the relationship between the two embedded beam patterns and present the corresponding model equations.


title: Anomalous magnetic moment of the muon with dynamical QCD+QED

abstract: The current $3.5\sigma$ discrepancy between experimental and Standard Model determinations of the anomalous magnetic moment of the muon $a_\mu=(g-2)/2$ can only be extended to the discovery $5\sigma$ regime through a reduction of both experimental and theoretical uncertainties. On the theory side, this means a determination of the hadronic vacuum polarisation (HVP) contribution to better than 0.5%, a level of precision that demands the inclusion of QCD + QED effects to properly understand how the behaviour of quarks are modified when their electric charges are turned on. The QCDSF collaboration has generated an ensemble of configurations with dynamical QCD and QED fields with the specific aim of studying flavour breaking effects arising from differences in the quark masses and charges in physical quantities. Here we study these effects in a calculation of HVP around the SU(3) symmetric point. Furthermore, by performing partially-quenched simulations we are able to cover a larger range of quark masses and charges on these configurations and then fit the results to an SU(3) flavour breaking expansion. Subsequently, this allows for an extrapolation to the physical point.


title: A new approach to norm inequalities on weighted and variable Hardy spaces

abstract: We give new proofs of Hardy space estimates for fractional and singular integral operators on weighted and variable exponent Hardy spaces. Our proofs consist of several interlocking ideas: finite atomic decompositions in terms of $L^\infty$ atoms, vector-valued inequalities for maximal and other operators, and Rubio de Francia extrapolation. Many of these estimates are not new, but we give new and substantially simpler proofs, which in turn significantly simplifies the proofs of the Hardy spaces inequalities.