title: String-to-String Interpretations with Polynomial-Size Output

abstract: String-to-string MSO interpretations are like Courcelle's MSO transductions, except that a single output position can be represented using a tuple of input positions instead of just a single input position. In particular, the output length is polynomial in the input length, as opposed to MSO transductions, which have output of linear length. We show that string-to-string MSO interpretations are exactly the polyregular functions. The latter class has various characterizations, one of which is that it consists of the string-to-string functions recognized by pebble transducers. Our main result implies the surprising fact that string-to-string MSO interpretations are closed under composition.


title: Ridesharing with Driver Location Preferences

abstract: We study revenue-optimal pricing and driver compensation in ridesharing platforms when drivers have heterogeneous preferences over locations. If a platform ignores drivers' location preferences, it may make inefficient trip dispatches; moreover, drivers may strategize so as to route towards their preferred locations. In a model with stationary and continuous demand and supply, we present a mechanism that incentivizes drivers to both (i) report their location preferences truthfully and (ii) always provide service. In settings with unconstrained driver supply or symmetric demand patterns, our mechanism achieves (full-information) first-best revenue. Under supply constraints and unbalanced demand, we show via simulation that our mechanism improves over existing mechanisms and has performance close to the first-best.


title: Graph Neural Tangent Kernel: Fusing Graph Neural Networks with Graph Kernels

abstract: While graph kernels (GKs) are easy to train and enjoy provable theoretical guarantees, their practical performances are limited by their expressive power, as the kernel function often depends on hand-crafted combinatorial features of graphs. Compared to graph kernels, graph neural networks (GNNs) usually achieve better practical performance, as GNNs use multi-layer architectures and non-linear activation functions to extract high-order information of graphs as features. However, due to the large number of hyper-parameters and the non-convex nature of the training procedure, GNNs are harder to train. Theoretical guarantees of GNNs are also not well-understood. Furthermore, the expressive power of GNNs scales with the number of parameters, and thus it is hard to exploit the full power of GNNs when computing resources are limited. The current paper presents a new class of graph kernels, Graph Neural Tangent Kernels (GNTKs), which correspond to infinitely wide multi-layer GNNs trained by gradient descent. GNTKs enjoy the full expressive power of GNNs and inherit advantages of GKs. Theoretically, we show GNTKs provably learn a class of smooth functions on graphs. Empirically, we test GNTKs on graph classification datasets and show they achieve strong performance.


title: De Giorgi type results for equations with nonlocal lower-order terms

abstract: It is known that the De Giorgi's conjecture does not hold in two dimensions for semilinear elliptic equations with a nonzero drift, in general, $$ \Delta u+ q\cdot \nabla u+f(u)=0 \ \ \text{in } \ \ \mathbb R^2, $$ when $q=(0,-c)$ for $c\neq 0$. This equation arises in the modeling of Bunsen burner flames. Bunsen flames are usually made of two flames: a diffusion flame and a premixed flame. In this article, we prove De Giorgi type results, and stability conjecture, for the following local-nonlocal counterpart of the above equation (with a nonlocal premixed flame) in two dimensions, $$\Delta u + c L[u] + f(u)=0 \quad \text{in} \ \ \mathbb R^n, $$ when $L$ is a nonlocal operator, $f\in C^1(\mathbb R)$ and $c\in\mathbb R^+$. In addition, we provide a priori estimates for the above equation, when $n\ge 1$, with various jumping kernels. The operator $\Delta+cL$ is an infinitesimal generator of jump-diffusion processes in the context of probability theory.


title: Sinkhorn Barycenters with Free Support via Frank-Wolfe Algorithm

abstract: We present a novel algorithm to estimate the barycenter of arbitrary probability distributions with respect to the Sinkhorn divergence. Based on a Frank-Wolfe optimization strategy, our approach proceeds by populating the support of the barycenter incrementally, without requiring any pre-allocation. We consider discrete as well as continuous distributions, proving convergence rates of the proposed algorithm in both settings. Key elements of our analysis are a new result showing that the Sinkhorn divergence on compact domains has Lipschitz continuous gradient with respect to the Total Variation and a characterization of the sample complexity of Sinkhorn potentials. Experiments validate the effectiveness of our method in practice.


title: Modeling Uncertainty by Learning a Hierarchy of Deep Neural Connections

abstract: Modeling uncertainty in deep neural networks, despite recent important advances, is still an open problem. Bayesian neural networks are a powerful solution, where the prior over network weights is a design choice, often a normal distribution or other distribution encouraging sparsity. However, this prior is agnostic to the generative process of the input data, which might lead to unwarranted generalization for out-of-distribution tested data. We suggest the presence of a confounder for the relation between the input data and the discriminative function given the target label. We propose an approach for modeling this confounder by sharing neural connectivity patterns between the generative and discriminative networks. This approach leads to a new deep architecture, where networks are sampled from the posterior of local causal structures, and coupled into a compact hierarchy. We demonstrate that sampling networks from this hierarchy, proportionally to their posterior, is efficient and enables estimating various types of uncertainties. Empirical evaluations of our method demonstrate significant improvement compared to state-of-the-art calibration and out-of-distribution detection methods.


title: Persistent homology detects curvature

abstract: In topological data analysis, persistent homology is used to study the "shape of data". Persistent homology computations are completely characterized by a set of intervals called a bar code. It is often said that the long intervals represent the "topological signal" and the short intervals represent "noise". We give evidence to dispute this thesis, showing that the short intervals encode geometric information. Specifically, we prove that persistent homology detects the curvature of disks from which points have been sampled. We describe a general computational framework for solving inverse problems using the average persistence landscape, a continuous mapping from metric spaces with a probability measure to a Hilbert space. In the present application, the average persistence landscapes of points sampled from disks of constant curvature results in a path in this Hilbert space which may be learned using standard tools from statistical and machine learning.


title: Supernovae SDSS Host Metallicity Catalog and Resulting Metallicity Distributions by SNe Type

abstract: In support of ongoing projects we have compiled a database of supernovae with host metallicities determined from spectroscopy in the Sloan Digital Sky Survey (SDSS) which we are providing as a catalog for the general use of the community. Here we provide metallicity values in the KK04 (Kobulnicky & Kewley 2004), T04 (Tremonti et al. 2004), and D16 (Dopita et al. 2016) diagnostics as well as a compilation of the relevant line fluxes to allow users to compute their own metallicities. The sample presented here has approximately 1700 objects with a metallicity in at least one diagnostic, with about 1200 objects in KK04 of which over 300 are core-collapse SNe or about twice the number of core-collapse objects present in the compilation of Kelly & Kirshner (2012) reflecting the recent considerable improvement in SNe detection capabilities. We find some differences in the metallicity distributions between the subtypes of both Type II and Type Ia SNe and speculate that metallicity may play a role in if Type II SNe events happen as IIb & IIn or as standard Type II.


title: Dispersion relations for $\gamma^*\gamma^*\to\pi\pi$: helicity amplitudes, subtractions, and anomalous thresholds

abstract: We present a comprehensive analysis of the dispersion relations for the doubly-virtual process $\gamma^*\gamma^*\to\pi\pi$. Starting from the Bardeen-Tung-Tarrach amplitudes, we first derive the kernel functions that define the system of Roy-Steiner equations for the partial-wave helicity amplitudes. We then formulate the solution of these partial-wave dispersion relations in terms of Omn\`es functions, with special attention paid to the role of subtraction constants as critical for the application to hadronic light-by-light scattering. In particular, we explain for the first time why for some amplitudes the standard Muskhelishvili-Omn\`es solution applies, while for others a modified approach based on their left-hand cut is required unless subtractions are introduced. In the doubly-virtual case, the analytic structure of the vector-resonance partial waves then gives rise to anomalous thresholds, even for space-like virtualities. We develop a strategy to account for these effects in the numerical solution, illustrated in terms of the $D$-waves in $\gamma^*\gamma^*\to\pi\pi$, which allows us to predict the doubly-virtual responses of the $f_2(1270)$ resonance. In general, our results form the basis for the incorporation of two-meson intermediate states into hadronic light-by-light scattering beyond the $S$-wave contribution.


title: Cross-section measurements of radiative proton-capture reactions in $^{112}$Cd at energies of astrophysical interest

abstract: Reactions involving the group of nuclei commonly known as p nuclei are part of the nucleosynthetic mechanisms at astrophysical sites. The $^{113}$In nucleus is such a case with several open questions regarding its origin at extreme stellar environments. In this work, the experimental study of the cross sections of the radiative proton-capture reaction $^{112}$Cd(p,$\gamma$)$^{113}$In is attempted for the first time at energies lying inside the Gamow window with an isotopically enriched $^{112}$Cd target. Two different techniques, the in-beam $\gamma$-ray spectroscopy and the activation method, have been applied. The latter method is required to account for the presence of a low-lying $^{113}$In isomer at 392 keV having a halflife of $\approx 100$ min. From the cross sections, the astrophysical S factors and the isomeric ratios have been additionally deduced. The experimental results are compared to detailed Hauser-Feshbach theoretical calculations using TALYS, and discussed in terms of their significance to the optical model potential involved.


title: Exploiting Uncertainty of Loss Landscape for Stochastic Optimization

abstract: We introduce novel variants of momentum by incorporating the variance of the stochastic loss function. The variance characterizes the confidence or uncertainty of the local features of the averaged loss surface across the i.i.d. subsets of the training data defined by the mini-batches. We show two applications of the gradient of the variance of the loss function. First, as a bias to the conventional momentum update to encourage conformity of the local features of the loss function (e.g. local minima) across mini-batches to improve generalization and the cumulative training progress made per epoch. Second, as an alternative direction for "exploration" in the parameter space, especially, for non-convex objectives, that exploits both the optimistic and pessimistic views of the loss function in the face of uncertainty. We also introduce a novel data-driven stochastic regularization technique through the parameter update rule that is model-agnostic and compatible with arbitrary architectures. We further establish connections to probability distributions over loss functions and the REINFORCE policy gradient update with baseline in RL. Finally, we incorporate the new variants of momentum proposed into Adam, and empirically show that our methods improve the rate of convergence of training based on our experiments on the MNIST and CIFAR-10 datasets.


title: Sr-induced dipole scatter in BST: Insights from MD simulations using a transferable bond valence-based interatomic potential

abstract: In order to design next-generation ferroelectrics, a microscopic understanding of their macroscopic properties is critical. One means to achieving an atomistic description of ferroelectric and dielectric phenomena is classical molecular dynamics simulations. Previously, we have shown that interatomic potentials based on the bond valence molecular dynamics (BVMD) method can be used to study structural phase transitions, ferroelectric domain nucleation, and domain wall migration in several perovskite oxides and fixed-composition binary and ternary alloys. Most modern devices, however, use variable-composition perovskite oxide alloys such as Ba$_x$Sr$_{1-x}$TiO$_3$ (BST). In this paper, we extend our bond valence approach to BST solid solutions and, in so doing, show that the potential parameters for each element are transferable between materials with different $x$. Using this potential, we perform BVMD simulations investigating the temperature and composition dependence of the lattice constants, Ti displacements, and ferroelectric polarization of BST and find that our predictions match experiments and first-principles theory. Additionally, based on a detailed analysis of local dipole distributions, we demonstrate that substitution of Sr for Ba scrambles dipoles, reduces global polarization, and enhances the order-disorder character of the ferroelectric-paraelectric phase transition.


title: Primordial black hole formation with non-Gaussian curvature perturbations

abstract: In the context of transient constant-roll inflation near a local maximum, we derive the non-perturbative field redefinition that relates a Gaussian random field with the true non-Gaussian curvature perturbation. Our analysis shows the emergence of a new critical amplitude $\zeta_*$, corresponding to perturbations that prevent the inflaton from overshooting the local maximum, thus becoming trapped in the false minimum of the potential. For potentials with a mild curvature at the local maximum (and thus small non-Gaussianity), we recover the known perturbative field redefinition. We apply these results to the formation of primordial black holes, and discuss the cases for which $\zeta_*$ is smaller or of the same order than the critical value for collapse of spherically symmetric overdensities. In the latter case, we present a simple potential for which the power spectrum needs an amplitude 10 times smaller that in the Gaussian case for producing a sizeable amount of primordial black holes.


title: Parametrised Algorithms for Directed Modular Width

abstract: Many well-known NP-hard algorithmic problems on directed graphs resist efficient parametrisations with most known width measures for directed graphs, such as directed treewidth, DAG-width, Kelly-width and many others. While these focus on measuring how close a digraph is to an oriented tree resp. a directed acyclic graph, in this paper, we investigate directed modular width as a parameter, which is closer to the concept of clique-width. We investigate applications of modular decompositions of directed graphs to a wide range of algorithmic problems and derive FPT-algorithms for several well-known digraph-specific NP-hard problems, namely minimum (weight) directed feedback vertex set, minimum (weight) directed dominating set, digraph colouring, directed Hamiltonian path/cycle, partitioning into paths, (capacitated) vertex-disjoint directed paths, and the directed subgraph homeomorphism problem. The latter yields a polynomial-time algorithm for detecting topological minors in digraphs of bounded directed modular width. Finally we illustrate that also other structural digraph parameters, such as the directed pathwidth and the cycle-rank can be computed efficiently using directed modular width as a parameter.


title: A qubit regularization of the $O(3)$ sigma model

abstract: We construct a qubit regularization of the $O(3)$ non-linear sigma model in two and three spatial dimensions using a quantum Hamiltonian with two qubits per lattice site. Using a worldline formulation and worm algorithms, we show that in two spatial dimensions our model has a quantum critical point where the well-known scale-invariant physics of the three-dimensional Wilson-Fisher fixed point is reproduced. In three spatial dimensions, we recover mean-field critical exponents at a similar quantum critical point. These results show that our qubit Hamiltonian is in the same universality class as the traditional classical lattice model close to the critical points. Simple modifications to our model also allow us to study the physics of traditional lattice models with $O(2)$ and $Z_2$ symmetries close to the corresponding critical points.


title: Near-Term Quantum-Classical Associative Adversarial Networks

abstract: We introduce a new hybrid quantum-classical adversarial machine learning architecture called a quantum-classical associative adversarial network (QAAN). This architecture consists of a classical generative adversarial network with a small auxiliary quantum Boltzmann machine that is simultaneously trained on an intermediate layer of the discriminator of the generative network. We numerically study the performance of QAANs compared to their classical counterparts on the MNIST and CIFAR-10 data sets, and show that QAANs attain a higher quality of learning when evaluated using the Inception score and the Fr\'{e}chet Inception distance. As the QAAN architecture only relies on sampling simple local observables of a small quantum Boltzmann machine, this model is particularly amenable for implementation on the current and next generations of quantum devices.


title: Phase transition in the Kolkata Paise Restaurant problem

abstract: A novel phase transition behaviour is observed in the Kolkata Paise Restaurant (KPR) problem where large number ($N$) of agents or customers collectively (and iteratively) learn to choose among the $N$ restaurants where she would expect to be alone that evening and would get the only dish available there (or may get randomly picked up if more than one agent arrive there that evening). The players are expected to evolve their strategy such that the publicly available information about past crowd in different restaurants can be utilized and each of them is able to make the best minority choice. For equally ranked restaurants we follow two crowd-avoiding strategies: Strategy I, where each of the $n_i(t)$ number of agents arriving at the $i$-th restaurant on the $t$-th evening goes back to the same restaurant on the next evening with probability $[n_i(t)]^{-\alpha}$, while in Strategy II, with probability $p$, when $n_i(t) > 1$. We study the steady state ($t$-independent) utilization fraction $f:(1-f)$ giving the steady state (wastage) fraction of restaurants going without any customer in any particular evening. With both the strategies we find, near $\alpha_c=0_+$ (in strategy I) or $p=1_-$ (in strategy II), the steady state wastage fraction $(1-f)\propto(\alpha - \alpha_c)^{\beta}$ or $(p_c - p)^\beta$ with $\beta \simeq 0.8, 0.87, 1.0$ and the convergence time $\tau$ (for $f(t)$ becoming independent of $t$) varies as $\tau\propto{(\alpha-\alpha_c)}^{-\gamma}$ or ${(p_c-p)}^{-\gamma}$, with $\gamma \simeq 1.18, 1.11, 1.05$ in infinite-dimension (rest of the $N-1$ neighboring restaurants), three-dimension ($6$ neighbors) and two-dimension ($4$ neighbors) respectively.


title: Convergence of uniform triangulations under the Cardy embedding

abstract: We consider an embedding of planar maps into an equilateral triangle $\Delta$ which we call the Cardy embedding. The embedding is a discrete approximation of a conformal map based on percolation observables that are used in Smirnov's proof of Cardy's formula. Under the Cardy embedding, the planar map induces a metric and an area measure on $\Delta$ and a boundary measure on $\partial \Delta$. We prove that for uniformly sampled triangulations, the metric and the measures converge jointly in the scaling limit to the Brownian disk conformally embedded into $\Delta$ (i.e., to the $\sqrt{8/3}$-Liouville quantum gravity disk). As part of our proof, we prove scaling limit results for critical site percolation on the uniform triangulations, in a quenched sense. In particular, we establish the scaling limit of the percolation crossing probability for a uniformly sampled triangulation with four boundary marked points.


title: An attention-based multi-resolution model for prostate whole slide imageclassification and localization

abstract: Histology review is often used as the `gold standard' for disease diagnosis. Computer aided diagnosis tools can potentially help improve current pathology workflows by reducing examination time and interobserver variability. Previous work in cancer grading has focused mainly on classifying pre-defined regions of interest (ROIs), or relied on large amounts of fine-grained labels. In this paper, we propose a two-stage attention-based multiple instance learning model for slide-level cancer grading and weakly-supervised ROI detection and demonstrate its use in prostate cancer. Compared with existing Gleason classification models, our model goes a step further by utilizing visualized saliency maps to select informative tiles for fine-grained grade classification. The model was primarily developed on a large-scale whole slide dataset consisting of 3,521 prostate biopsy slides with only slide-level labels from 718 patients. The model achieved state-of-the-art performance for prostate cancer grading with an accuracy of 85.11\% for classifying benign, low-grade (Gleason grade 3+3 or 3+4), and high-grade (Gleason grade 4+3 or higher) slides on an independent test set.


title: AssembleNet: Searching for Multi-Stream Neural Connectivity in Video Architectures

abstract: Learning to represent videos is a very challenging task both algorithmically and computationally. Standard video CNN architectures have been designed by directly extending architectures devised for image understanding to include the time dimension, using modules such as 3D convolutions, or by using two-stream design to capture both appearance and motion in videos. We interpret a video CNN as a collection of multi-stream convolutional blocks connected to each other, and propose the approach of automatically finding neural architectures with better connectivity and spatio-temporal interactions for video understanding. This is done by evolving a population of overly-connected architectures guided by connection weight learning. Architectures combining representations that abstract different input types (i.e., RGB and optical flow) at multiple temporal resolutions are searched for, allowing different types or sources of information to interact with each other. Our method, referred to as AssembleNet, outperforms prior approaches on public video datasets, in some cases by a great margin. We obtain 58.6% mAP on Charades and 34.27% accuracy on Moments-in-Time.


title: Generalization Bounds of Stochastic Gradient Descent for Wide and Deep Neural Networks

abstract: We study the training and generalization of deep neural networks (DNNs) in the over-parameterized regime, where the network width (i.e., number of hidden nodes per layer) is much larger than the number of training data points. We show that, the expected $0$-$1$ loss of a wide enough ReLU network trained with stochastic gradient descent (SGD) and random initialization can be bounded by the training loss of a random feature model induced by the network gradient at initialization, which we call a neural tangent random feature (NTRF) model. For data distributions that can be classified by NTRF model with sufficiently small error, our result yields a generalization error bound in the order of $\tilde{\mathcal{O}}(n^{-1/2})$ that is independent of the network width. Our result is more general and sharper than many existing generalization error bounds for over-parameterized neural networks. In addition, we establish a strong connection between our generalization error bound and the neural tangent kernel (NTK) proposed in recent work.


title: What Can Neural Networks Reason About?

abstract: Neural networks have succeeded in many reasoning tasks. Empirically, these tasks require specialized network structures, e.g., Graph Neural Networks (GNNs) perform well on many such tasks, but less structured networks fail. Theoretically, there is limited understanding of why and when a network structure generalizes better than others, although they have equal expressive power. In this paper, we develop a framework to characterize which reasoning tasks a network can learn well, by studying how well its computation structure aligns with the algorithmic structure of the relevant reasoning process. We formally define this algorithmic alignment and derive a sample complexity bound that decreases with better alignment. This framework offers an explanation for the empirical success of popular reasoning models, and suggests their limitations. As an example, we unify seemingly different reasoning tasks, such as intuitive physics, visual question answering, and shortest paths, via the lens of a powerful algorithmic paradigm, dynamic programming (DP). We show that GNNs align with DP and thus are expected to solve these tasks. On several reasoning tasks, our theory is supported by empirical results.


title: Deep Learning for Direct Hybrid Precoding in Millimeter Wave Massive MIMO Systems

abstract: This paper proposes a novel neural network architecture, that we call an auto-precoder, and a deep-learning based approach that jointly senses the millimeter wave (mmWave) channel and designs the hybrid precoding matrices with only a few training pilots. More specifically, the proposed machine learning model leverages the prior observations of the channel to achieve two objectives. First, it optimizes the compressive channel sensing vectors based on the surrounding environment in an unsupervised manner to focus the sensing power on the most promising spatial directions. This is enabled by a novel neural network architecture that accounts for the constraints on the RF chains and models the transmitter/receiver measurement matrices as two complex-valued convolutional layers. Second, the proposed model learns how to construct the RF beamforming vectors of the hybrid architectures directly from the projected channel vector (the received sensing vector). The auto-precoder neural network that incorporates both the channel sensing and beam prediction is trained end-to-end as a multi-task classification problem. Thanks to this design methodology that leverages the prior channel observations and the implicit awareness about the surrounding environment/user distributions, the proposed approach significantly reduces the training overhead compared to classical (non-machine learning) solutions. For example, for a system of 64 transmit and 64 receive antennas, with 3 RF chains at both sides, the proposed solution needs only 8 or 16 channel training pilots to directly predict the RF beamforming/combining vectors of the hybrid architectures and achieve near-optimal achievable rates. This highlights a promising solution for the channel estimation and hybrid precoding design problem in mmWave and massive MIMO systems.


title: Swampland Bounds on the Abelian Gauge Sector

abstract: We derive bounds on the number of abelian gauge group factors in six-dimensional gravitational theories with minimal supersymmetry and in their F-theoretic realisations. These bounds follow by requiring consistency of certain BPS strings in the spectrum of the theory, as recently proposed in the literature. Under certain assumptions this approach constrains the number of abelian gauge group factors in six-dimensional supergravity theories with at least one tensor multiplet to be $N \leq 20$ (or $N \leq 22$ in absence of charged matter). For any geometric F-theory realisation with at least one tensor multiplet we establish the bound $N \leq 16$ by demanding unitarity of a heterotic solitonic string which exists even in absence of a perturbative heterotic dual. This result extends to four-dimensional F-theory vacua on any blowup of a rational fibration. Our findings lead to universal bounds on the rank of the Mordell-Weil group of elliptically fibered Calabi-Yau threefolds.


title: On Network Design Spaces for Visual Recognition

abstract: Over the past several years progress in designing better neural network architectures for visual recognition has been substantial. To help sustain this rate of progress, in this work we propose to reexamine the methodology for comparing network architectures. In particular, we introduce a new comparison paradigm of distribution estimates, in which network design spaces are compared by applying statistical techniques to populations of sampled models, while controlling for confounding factors like network complexity. Compared to current methodologies of comparing point and curve estimates of model families, distribution estimates paint a more complete picture of the entire design landscape. As a case study, we examine design spaces used in neural architecture search (NAS). We find significant statistical differences between recent NAS design space variants that have been largely overlooked. Furthermore, our analysis reveals that the design spaces for standard model families like ResNeXt can be comparable to the more complex ones used in recent NAS work. We hope these insights into distribution analysis will enable more robust progress toward discovering better networks for visual recognition.


title: Image classification using quantum inference on the D-Wave 2X

abstract: We use a quantum annealing D-Wave 2X computer to obtain solutions to NP-hard sparse coding problems. To reduce the dimensionality of the sparse coding problem to fit on the quantum D-Wave 2X hardware, we passed downsampled MNIST images through a bottleneck autoencoder. To establish a benchmark for classification performance on this reduced dimensional data set, we used an AlexNet-like architecture implemented in TensorFlow, obtaining a classification score of $94.54 \pm 0.7 \%$. As a control, we showed that the same AlexNet-like architecture produced near-state-of-the-art classification performance $(\sim 99\%)$ on the original MNIST images. To obtain a set of optimized features for inferring sparse representations of the reduced dimensional MNIST dataset, we imprinted on a random set of $47$ image patches followed by an off-line unsupervised learning algorithm using stochastic gradient descent to optimize for sparse coding. Our single-layer of sparse coding matched the stride and patch size of the first convolutional layer of the AlexNet-like deep neural network and contained $47$ fully-connected features, $47$ being the maximum number of dictionary elements that could be embedded onto the D-Wave $2$X hardware. Recent work suggests that the optimal level of sparsity corresponds to a critical value of the trade-off parameter associated with a putative second order phase transition, an observation supported by a free energy analysis of D-Wave energy states. When the sparse representations inferred by the D-Wave $2$X were passed to a linear support vector machine, we obtained a classification score of $95.68\%$. Thus, on this problem, we find that a single-layer of quantum inference is able to outperform a standard deep neural network architecture.


title: A generalisation of the honeycomb dimer model to higher dimensions

abstract: Linde, Moore, and Nordahl introduced a generalisation of the honeycomb dimer model to higher dimensions. The purpose of this article is to describe a number of structural properties of this generalised model. First, it is shown that the samples of the model are in one-to-one correspondence with the perfect matchings of a hypergraph. This leads to a generalised Kasteleyn theory: the partition function of the model equals the Cayley hyperdeterminant of the adjacency hypermatrix of the hypergraph. Second, we prove an identity which relates the covariance matrix of the random height function directly to the random geometrical structure of the model. This identity is known in the planar case but is new for higher dimensions. It relies on a more explicit formulation of Sheffield's cluster swap which is made possible by the structure of the honeycomb dimer model. Finally, we use the special properties of this explicit cluster swap to give a new and simplified proof of strict convexity of the surface tension in this case.


title: Crystalline droplets with emergent topological color-charge in many-body systems with sign-changing interactions

abstract: We introduce a novel type of self-bound droplet which carries an emergent color charge. We consider a system of particles hopping on a lattice and interacting via a commensurately sign-changing potential which is attractive at a short range. The droplet formation is heralded by spontaneous crystallization into topologically distinct domains. This endows each droplet with an emergent color charge governing their mutual interactions: attractive for equal colors and repulsive otherwise. The number of allowed colors is fixed only by the discrete spatial symmetries of the sign-changing part of the interaction potential. With increasing interaction range, the droplets become progressively more mobile, with their color charge still being energetically protected, allowing for nontrivial viscous dynamics of the interacting droplet plasmas formed during cooling. Sign-changing potentials with a short-range attraction appear quite naturally for light-mediated interactions and we concretely propose a realization in state-of-the-art experiments with cold atoms in a multimode optical cavity.


title: A Model for Energy Buildup and Eruption Onset in Coronal Mass Ejections

abstract: Coronal mass ejections (CMEs) and eruptive flares (EFs) are the most energetic explosions in the solar system. Their underlying origin is the free energy that builds up slowly in the sheared magnetic field of a filament channel. We report the first end-to-end numerical simulation of a CME/EF, from zero-free-energy initial state through filament-channel formation to violent eruption, driven solely by the magnetic-helicity condensation process. Helicity is the topological measure of linkages between magnetic flux systems, and is conserved in the corona, building up inexorably until it is ejected into interplanetary space. Numerous investigations have demonstrated that helicity injected by small-scale vortical motions, such as those observed in the photosphere, undergoes an inverse cascade from small scales to large, `condensing' at magnetic-polarity boundaries. Our new results verify that this process forms a filament channel within a compact bipolar region embedded in a background dipole field, and show for the first time that a fast CME eventually occurs via the magnetic-breakout mechanism. We further show that the trigger for explosive eruption is reconnection onset in the flare current sheet that develops above the polarity inversion line: this reconnection forms flare loops below the sheet and a CME flux rope above, and initiates high-speed outward flow of the CME. Our findings have important implications for magnetic self-organization and explosive behavior in solar and other astrophysical plasmas, and for understanding and predicting explosive solar activity.


title: Evolution of Magnetic Order from the Localized to the Itinerant Limit

abstract: Quantum materials that feature magnetic long-range order often reveal complex phase diagrams when localized electrons become mobile. In many materials magnetism is rapidly suppressed as electronic charges dissolve into the conduction band. In materials where magnetism persists, it is unclear how the magnetic properties are affected. Here we study the evolution of the magnetic structure in Nd(1-x)Ce(x)CoIn(5) from the localized to the highly itinerant limit. We observe two magnetic ground states inside a heavy-fermion phase that are detached from unconventional superconductivity. The presence of two different magnetic phases provides evidence that increasing charge delocalization affects the magnetic interactions via anisotropic band hybridization.


title: Positivity of Multiparameter Persistence Diagrams and Bottleneck Stability

abstract: Persistent homology studies the birth and death of cycles in a parameterized family of spaces. In this paper, we study the birth and death of cycles in a multifiltration of a chain complex with the goal of producing a persistence diagram that satisfies bottleneck stability.


title: Video from Stills: Lensless Imaging with Rolling Shutter

abstract: Because image sensor chips have a finite bandwidth with which to read out pixels, recording video typically requires a trade-off between frame rate and pixel count. Compressed sensing techniques can circumvent this trade-off by assuming that the image is compressible. Here, we propose using multiplexing optics to spatially compress the scene, enabling information about the whole scene to be sampled from a row of sensor pixels, which can be read off quickly via a rolling shutter CMOS sensor. Conveniently, such multiplexing can be achieved with a simple lensless, diffuser-based imaging system. Using sparse recovery methods, we are able to recover 140 video frames at over 4,500 frames per second, all from a single captured image with a rolling shutter sensor. Our proof-of-concept system uses easily-fabricated diffusers paired with an off-the-shelf sensor. The resulting prototype enables compressive encoding of high frame rate video into a single rolling shutter exposure, and exceeds the sampling-limited performance of an equivalent global shutter system for sufficiently sparse objects.


title: Correlation Effects in Quadrupole Insulators: a Quantum Monte Carlo Study

abstract: The quadrupole insulator, a high-order topological insulator, with on-site Hubbard interaction is numerically studied by large-scale projector quantum Monte Carlo (PQMC) simulations. The Green's function formalism is successfully used to characterize topological properties in interacting quadrupole insulators for the first time. We find that the topological quadrupole insulator is stable against weak interactions and turns into a trivial antiferromagnetic (AFM) insulator by a continuous topological phase transition (TPT) for strong interactions. The critical exponents related to the TPT are estimated to be $\nu=0.67(4)$, $\beta=0.40(2)$, which are distinct from those of the known AFM transitions and suggest a new universality class.


title: Explaining the ANITA Anomaly with Inelastic Boosted Dark Matter

abstract: We propose a new physics scenario in which the decay of a very heavy dark-matter candidate which does not interact with the neutrino sector could explain the two anomalous events recently reported by ANITA. The model is composed of two components of dark matter, an unstable dark-sector state, and a massive dark gauge boson. We assume that the heavier dark-matter particle of EeV-range mass is distributed over the galactic halo and disintegrates into a pair of lighter -- highly boosted -- dark-matter states in the present universe which reach and penetrate the Earth. The latter scatters {\it in}elastically off a nucleon and produces a heavier dark-sector unstable state which subsequently decays back to the lighter dark matter along with hadrons, which induce Extensive Air Showers, via on-/off-shell dark gauge boson. Depending on the mass hierarchy within the dark sector, either the dark gauge boson or the unstable dark-sector particle can be long-lived, hence transmitted significantly through the Earth. We study the angular distribution of the signal and show that our model favors emergence angles in the range $\sim 25^\circ -35^\circ$ if the associated parameter choices bear the situation where the mean free path of the boosted incident particle is much larger than the Earth diameter while its long-lived decay product has a decay length of dimensions comparable to the Earth radius. Our model, in particular, avoids any constraints from complementary neutrino searches such as IceCube or the Auger observatory.


title: Inverse Seesaw and Portal Dark Matter

abstract: We study the phenomenology of the inverse seesaw mechanism in the scalar-Higgs portal dark matter model. The model is an extension of the Standard Model including two additional neutrinos, a singlet scalar and a fermionic dark matter. We consider the inverse seesaw mechanism where the mass of 2 additional neutrinos are made dynamic by the singlet scalar. We found that the natural scale for the scalar vacuum expectation value is naturally close to the weak scale. Motivating by this fact, we focus on the possibility of the singlet scalar connecting with dark matter, i.e., the scalar is also the mediator between dark sector and the Standard Model. We perform a numerical analysis over the parameter space subject to the indirect and direct detection constraints. The feasible region of the parameter space will be discussed.


title: Modeling Theory of Mind in Multi-Agent Games Using Adaptive Feedback Control

abstract: A major challenge in cognitive science and AI has been to understand how autonomous agents might acquire and predict behavioral and mental states of other agents in the course of complex social interactions. How does such an agent model the goals, beliefs, and actions of other agents it interacts with? What are the computational principles to model a Theory of Mind (ToM)? Deep learning approaches to address these questions fall short of a better understanding of the problem. In part, this is due to the black-box nature of deep networks, wherein computational mechanisms of ToM are not readily revealed. Here, we consider alternative hypotheses seeking to model how the brain might realize a ToM. In particular, we propose embodied and situated agent models based on distributed adaptive control theory to predict actions of other agents in five different game theoretic tasks (Harmony Game, Hawk-Dove, Stag-Hunt, Prisoner's Dilemma and Battle of the Exes). Our multi-layer control models implement top-down predictions from adaptive to reactive layers of control and bottom-up error feedback from reactive to adaptive layers. We test cooperative and competitive strategies among seven different agent models (cooperative, greedy, tit-for-tat, reinforcement-based, rational, predictive and other's-model agents). We show that, compared to pure reinforcement-based strategies, probabilistic learning agents modeled on rational, predictive and other's-model phenotypes perform better in game-theoretic metrics across tasks. Our autonomous multi-agent models capture systems-level processes underlying a ToM and highlight architectural principles of ToM from a control-theoretic perspective.


title: Enriching Bibliographic Data by Combining String Matching and the Wikidata Knowledge Graph to Improve the Measurement of International Research Collaboration

abstract: Measuring international research collaboration is necessary when evaluating, for example, the efficacy of policy meant to increase cooperation between countries, but is currently very difficult as bibliographic records contain only affiliation data from which there is no standard method to identify the relevant countries. In this paper we describe a method to address this difficulty, and evaluate it using both general and domain-specific data sets.


title: Equation of State of Colloidal Membranes

abstract: In the presence of a non-adsorbing polymer, monodisperse rod-like colloids assemble into one-rod-length thick liquid-like monolayers, called colloidal membranes. The density of the rods within a colloidal membrane is determined by a balance between the osmotic pressure exerted by the enveloping polymer suspension and the repulsion between the colloidal rods. We developed a microfluidic device for continuously observing an isolated membrane while dynamically controlling the osmotic pressure of the polymer suspension. Using this technology we measured the membrane rod density over a range of osmotic pressures than is wider that what is accessible in equilibrium samples. With increasing density we observed a first-order phase transition, in which the in-plane membrane order transforms from a 2D fluid into a 2D solid. In the limit of low osmotic pressures, we measured the rate at which individual rods evaporate from the membrane. The developed microfluidic technique could have wide applicability for in situ investigation of various soft materials and how their properties depend on the solvent composition.


title: Cospectral Bipartite Graphs with the Same Degree Sequences but with Different Number of Large Cycles

abstract: Finding the multiplicity of cycles in bipartite graphs is a fundamental problem of interest in many fields including the analysis and design of low-density parity-check (LDPC) codes. Recently, Blake and Lin computed the number of shortest cycles ($g$-cycles, where $g$ is the girth of the graph) in a bi-regular bipartite graph, in terms of the degree sequences and the spectrum (eigenvalues of the adjacency matrix) of the graph [{\em IEEE Trans. Inform. Theory 64(10):6526--6535, 2018}]. This result was subsequently extended in [{\em IEEE Trans. Inform. Theory, accepted for publication, Dec. 2018}] to cycles of length $g+2, \ldots, 2g-2$, in bi-regular bipartite graphs, as well as $4$-cycles and $6$-cycles in irregular and half-regular bipartite graphs, with $g \geq 4$ and $g \geq 6$, respectively. In this paper, we complement these positive results with negative results demonstrating that the information of the degree sequences and the spectrum of a bipartite graph is, in general, insufficient to count (a) the $i$-cycles, $i \geq 2g$, in bi-regular graphs, (b) the $i$-cycles for any $i > g$, regardless of the value of $g$, and $g$-cycles for $g \geq 6$, in irregular graphs, and (c) the $i$-cycles for any $i > g$, regardless of the value of $g$, and $g$-cycles for $g \geq 8$, in half-regular graphs. To obtain these results, we construct counter-examples using the Godsil-McKay switching.


title: Private Hypothesis Selection

abstract: We provide a differentially private algorithm for hypothesis selection. Given samples from an unknown probability distribution $P$ and a set of $m$ probability distributions $\mathcal{H}$, the goal is to output, in a $\varepsilon$-differentially private manner, a distribution from $\mathcal{H}$ whose total variation distance to $P$ is comparable to that of the best such distribution (which we denote by $\alpha$). The sample complexity of our basic algorithm is $O\left(\frac{\log m}{\alpha^2} + \frac{\log m}{\alpha \varepsilon}\right)$, representing a minimal cost for privacy when compared to the non-private algorithm. We also can handle infinite hypothesis classes $\mathcal{H}$ by relaxing to $(\varepsilon,\delta)$-differential privacy. We apply our hypothesis selection algorithm to give learning algorithms for a number of natural distribution classes, including Gaussians, product distributions, sums of independent random variables, piecewise polynomials, and mixture classes. Our hypothesis selection procedure allows us to generically convert a cover for a class to a learning algorithm, complementing known learning lower bounds which are in terms of the size of the packing number of the class. As the covering and packing numbers are often closely related, for constant $\alpha$, our algorithms achieve the optimal sample complexity for many classes of interest. Finally, we describe an application to private distribution-free PAC learning.


title: Photoevaporation of Jeans-unstable molecular clumps

abstract: We study the photoevaporation of Jeans-unstable molecular clumps by isotropic FUV (6 eV $< {\rm h}\nu$ < 13.6 eV) radiation, through 3D radiative transfer hydrodynamical simulations implementing a non-equilibrium chemical network that includes the formation and dissociation of H$_2$. We run a set of simulations considering different clump masses ($M=10-200$ M$_\odot$) and impinging fluxes ($G_0=2\times 10^3-8\times 10^4$ in Habing units). In the initial phase, the radiation sweeps the clump as an R-type dissociation front, reducing the H$_2$ mass by a factor 40-90%. Then, a weak ($\mathcal{M}\simeq 2$) shock develops and travels towards the centre of the clump, which collapses while loosing mass from its surface. All considered clumps remain gravitationally unstable even if radiation rips off most of the clump mass, showing that external FUV radiation is not able to stop clump collapse. However, the FUV intensity regulates the final H$_2$ mass available for star formation: for example, for $G_0 < 10^4$ more than 10% of the initial clump mass survives. Finally, for massive clumps ($\sim 100$ M$_\odot$) the H$_2$ mass increases by 25-50% during the collapse, mostly because of the rapid density growth that implies a more efficient H$_2$ self-shielding.


title: The kinematics of star clusters undergoing gas expulsion in Newtonian and Milgromian dynamics

abstract: We study the kinematics of stars in clusters undergoing gas expulsion in standard Newtonian dynamics and also in Milgromian dynamics (MOND). Gas expulsion can explain the observed line-of-sight (LoS) velocity dispersion profile of NGC 2419 in Newtonian dynamics. For a given star formation efficiency (SFE), the shapes of the velocity dispersion profiles, which are normalised by the velocity dispersion at the projected half-mass radius, are almost indistinguishable for different SFE models in Newtonian dynamics. The velocity dispersion of a star cluster in the outer halo of a galaxy can indeed have a strong radial anisotropy in Newtonian dynamics after gas expulsion. MOND displays several different properties from Newtonian dynamics. In particular, the slope of the central velocity dispersion profile is less steep in MOND for the same SFE. Moreover, for a given SFE, more massive embedded cluster models result in more rapidly declining central velocity dispersion profiles for the final star clusters, while less massive embedded cluster models lead to flatter velocity dispersion profiles for the final products. The onset of the radial-orbit instability in post-gas-expulsion MOND models is discussed. SFEs as low as a few percent, typical of molecular clouds, lead to surviving ultra-diffuse objects. Gas expulsion alone is unlikely the physical mechanism for the observed velocity dispersion profile of NGC 2419 in MOND.


title: Discrete Time Crystals

abstract: Experimental advances have allowed for the exploration of nearly isolated quantum many-body systems whose coupling to an external bath is very weak. A particularly interesting class of such systems is those which do not thermalize under their own isolated quantum dynamics. In this review, we highlight the possibility for such systems to exhibit new non-equilibrium phases of matter. In particular, we focus on "discrete time crystals", which are many-body phases of matter characterized by a spontaneously broken discrete time translation symmetry. We give a definition of discrete time crystals from several points of view, emphasizing that they are a non-equilibrium phenomenon, which is stabilized by many-body interactions, with no analog in non-interacting systems. We explain the theory behind several proposed models of discrete time crystals, and compare a number of recent realizations, in different experimental contexts.

categories: cond-mat.str-el cond-mat.dis-nn cond-mat.quant-gas cond-mat.stat-mech


title: Horizon-AGN virtual observatory -- 2: Template-free estimates of galaxy properties from colours

abstract: Using the Horizon-AGN hydrodynamical simulation and self-organising maps (SOMs), we show how to compress the complex data structure of a cosmological simulation into a 2-d grid which is much easier to analyse. We first verify the tight correlation between the observed 0.3$\!-\!5\mu$m broad-band colours of Horizon-AGN galaxies and their high-resolution spectra. The correlation is found to extend to physical properties such as redshift, stellar mass, and star formation rate (SFR). This direct mapping from colour to physical parameter space is shown to work also after including photometric uncertainties that mimic the COSMOS survey. We then label the SOM grid with a simulated calibration sample and estimate redshift and SFR for COSMOS-like galaxies up to $z\sim3$. In comparison to state-of-the-art techniques based on synthetic templates, our method is comparable in performance but less biased at estimating redshifts, and significantly better at predicting SFRs. In particular our "data-driven" approach, in contrast to model libraries, intrinsically allows for the complexity of galaxy formation and can handle sample biases. We advocate that obtaining the calibration for this method should be one of the goals of next-generation galaxy surveys.


title: Predictions for the spatial distribution of the dust continuum emission in 1<z<5 star-forming galaxies

abstract: We present the first detailed study of the spatially-resolved dust continuum emission of simulated galaxies at 1<z<5. We run the radiative transfer code SKIRT on a sample of submillimeter-bright galaxies drawn from the Feedback in Realistic Environments (FIRE) project. These simulated galaxies reach Milky Way masses by z=2. Our modelling provides predictions for the full rest-frame far-ultraviolet-to-far-infrared spectral energy distributions of these simulated galaxies, as well as 25-pc-resolution maps of their emission across the wavelength spectrum. The derived morphologies are notably different in different wavebands, with the same galaxy often appearing clumpy and extended in the far-ultraviolet yet an ordered spiral at far-infrared wavelengths. The observed-frame 870-$\mu$m half-light radii of our FIRE-2 galaxies are ~0.5-4kpc, consistent with existing ALMA observations of galaxies with similarly high redshifts and stellar masses. In both simulated and observed galaxies, the dust continuum emission is generally more compact than the cold gas and the dust mass, but more extended than the stellar component. The most extreme cases of compact dust emission seem to be driven by particularly compact recent star-formation, which generates steep dust temperature gradients. Our results confirm that the spatial extent of the dust continuum emission is sensitive to both the dust mass and SFR distributions.


title: Ultra-diffuse galaxies without dark matter

abstract: I develop a high velocity galaxy collision model to explain a rare but puzzling phenomenon, namely the apparent existence of ultra-diffuse galaxies with little dark matter. Predictions include simultaneous triggering of overpressured dense clouds to form luminous old globular clusters, a protogroup environment to generate high relative velocities of the initially gas-rich galaxies in the early universe, and spatially separated dark halos, possibly detectable via gravitational lensing and containing relic low metallicity stars with enhanced [alpha/Fe] at ultralow surface brightness.


title: Enhanced cluster lensing models with measured galaxy kinematics

abstract: We present an improved determination of the total mass distribution of three CLASH/HFF massive clusters, MACS J1206.2-0847 (z=0.44), MACS J0416.1-2403 (z=0.40), Abell S1063 (z=0.35). We specifically reconstruct the sub-halo mass component with robust stellar kinematics information of cluster galaxies, in combination with precise strong lensing models based on large samples of spectroscopically identified multiple images. We use VLT/MUSE integral-field spectroscopy in the cluster cores to measure the stellar velocity dispersion, $\sigma$, of 40-60 member galaxies per cluster, covering 4-5 magnitudes to $m_{F160W}\simeq 21.5$. We verify the robustness and quantify the accuracy of the velocity dispersion measurements with extensive spectral simulations. With these data, we determine the normalization and slope of the galaxy $L\mbox{-}\sigma$ Faber-Jackson relation in each cluster and use these parameters as a prior for the scaling relations of the sub-halo population in the mass distribution modeling. When compared to our previous lens models, the inclusion of member galaxies' kinematics provides a similar precision in reproducing the positions of the multiple images. However, the inherent degeneracy between the central effective velocity dispersion, $\sigma_0$, and truncation radius, $r_{cut}$, of sub-halos is strongly reduced, thus significantly alleviating possible systematics in the measurements of sub-halo masses. The three independent determinations of the $\sigma_0\mbox{-}r_{cut}$ scaling relation in each cluster are found to be fully consistent, enabling a statistical determination of sub-halo sizes as a function of $\sigma_0$, or halo masses. We derive galaxy central velocity dispersion functions of the three clusters and found them in agreement with each other. Sub-halo mass functions determined with this method can be compared with those obtained from cosmological simulations.


title: Repeated faint quasinormal bursts in extreme-mass-ratio inspiral waveforms: Evidence from frequency-domain scalar self-force calculations on generic Kerr orbits

abstract: We report development of a code to calculate the scalar self-force on a scalar-charged particle moving on generic bound orbits in the Kerr spacetime. The scalar self-force model allows rapid development of computational techniques relevant to generic gravitational extreme-mass-ratio inspirals (EMRIs). Our frequency-domain calculations are made with arbitrary numerical precision code written in \textsc{Mathematica}. We extend spectral source integration techniques to the Kerr spacetime, increasing computational efficiency. We model orbits with nearly arbitrary inclinations $0\leq\iota<\pi/2$ and eccentricities up to $e \lesssim 0.8$. This effort extends earlier work by Warburton and Barack where motion was restricted to the equatorial plane or to inclined spherical orbits. Consistent with a recent discovery by Thornburg and Wardell \cite{ThorWard17} in time-domain calculations, we observe self-force oscillations during the radially-outbound portion of highly eccentric orbits around a rapidly rotating black hole. As noted previously, these oscillations reflect coupling into the self-force by quasinormal modes excited during pericenter passage. Our results confirm the effect with a frequency-domain code. \emph{More importantly, we find that quasinormal bursts (QNBs) appear directly in the waveform following each periastron passage.} These faint bursts are shown to be a superposition of the least-damped overtone (i.e., fundamental) of at least four ($l=m \le 4$) quasinormal modes. Our results suggest that QNBs should appear in gravitational waveforms, and thus provide a gauge-invariant signal. Potentially observable in high signal-to-noise ratio EMRIs, QNBs would provide high-frequency components to the parameter estimation problem that would complement low-frequency elements of the waveform.


title: Astrophotonic Spectrographs

abstract: Astrophotonics is the application of photonic technologies to channel, manipulate, and disperse light from one or more telescopes to achieve scientific objectives in astronomy in an efficient and cost-effective way. Utilizing photonic advantage for astronomical spectroscopy is a promising approach to miniaturizing the next generation of spectrometers for large telescopes. It can be primarily attained by leveraging the two-dimensional nature of photonic structures on a chip or a set of fibers, thus reducing the size of spectroscopic instrumentation to a few centimeters and the weight to a few hundred grams. A wide variety of astrophotonic spectrometers is currently being developed, including arrayed waveguide gratings (AWGs), photonic echelle gratings (PEGs), and Fourier-transform spectrometer (FTS). These astrophotonic devices are flexible, cheaper to mass produce, easier to control, and much less susceptible to vibrations and flexure than conventional astronomical spectrographs. The applications of these spectrographs range from astronomy to biomedical analysis. This paper provides a brief review of this new class of astronomical spectrographs.


title: Study of the Mass-Ratio Distribution of Spectroscopic Binaries. II. The Boundaries of the Brown-Dwarf Desert as Seen with the APOGEE Spectroscopic Binaries

abstract: Analysis of APOGEE DR12 stellar radial-velocities by Troup et al. (2016) affirmed the existence of the well-known Brown-Dwarf Desert (BDD). They detected a dearth of spectroscopic binaries (SB) with periods shorter than $\sim 10$ - $30$ days and secondaries with masses in the range of $\sim0.01$ - $0.1\, M_{\odot}$. We reconsider here their sample of binaries, focusing on 116 systems on the main sequence of the Gaia color-magnitude diagram, with mostly K-dwarf primaries. Using our recently devised algorithm to analyze the mass-ratio distribution of a sample of SBs we confirm the BDD existence and delineate its boundaries. For the K-dwarf APOGEE $1$ - $25$ days binaries, the companion-mass range of the BDD is $\sim0.02$ - $0.2\, M_{\odot}$. The mass ratio distribution of the long-period ($25$ - $500$ days) binaries does not show any dearth at the $q$-range studied. Instead, their distribution displays a linear increase in $\log q$, implying a tendency towards low-$q$ values. The limits of the BDD do not coincide with the frequently used mass limits of the brown-dwarf population, sometimes defined as $0.013$ and $0.08\, M_{\odot}$, based on theoretically derived stellar minimum masses for burning deuterium and hydrogen in their cores. Trying to draw the boundaries of the desert, we suggest either a wedged or trapezoidal shape. We discuss briefly different scenarios that can account for the formation of the BDD, in terms of differentiating between stellar secondaries and planets in particular, and compare this desert to the Neptunian desert that can distinguish between Jovian planets and super Earths of short periods.


title: High order regularised symplectic integrator for collisional planetary systems

abstract: We present a new mixed variable symplectic (MVS) integrator for planetary systems, that fully resolve close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators, while also reducing the effective step size whenever two planets encounter. We use a high order MVS scheme such that it is possible to integrate with large time steps far away from close encounters. We show that this algorithm is able to resolve almost exact collisions (i.e with a mutual separation of a fraction of the physical radius) while using the same time-step as in weakly perturbed problem such as the Solar System. We demonstrate the long term behaviour on systems of six super-Earths experiencing strong scattering for 50 kyr. We compare our algorithm to hybrid methods such as MERCURY and show that for an equivalent cost we obtain much better energy conservation.


title: Towards a multi-input astrophotonic AWG spectrograph

abstract: Astrophotonics is the new frontier technology to develop diffraction-limited, miniaturized, and cost-effective instruments for the next generation of large telescopes. For various astronomical studies such as probing the early universe, observing in near infrared (NIR) is crucial. To address this, we are developing moderate resolution (R = 1500) on-chip astrophotonic spectrographs in the NIR bands (J Band: 1.1-1.4 $\mu m$; H band: 1.45-1.7 $\mu m$) using the concept of arrayed waveguide gratings (AWGs). We fabricate the AWGs using a silica-on-silicon substrate. The waveguides on these AWGs are 2 $\mu m$ wide and 0.1 $\mu m$ high Si3N4 core buried inside a 15 $\mu m$ thick SiO2 cladding. To make the maximal use of astrophotonic integration such as coupling the AWGs with multiple single-mode fibers coming from photonic lanterns or fiber Bragg gratings (FBGs), we require a multi-input AWG design. In a multi-input AWG, the output spectrum due to each individual input channel overlaps to produce a combined spectrum from all inputs. This on-chip combination of light effectively improves the signal-to-noise ratio as compared to spreading the photons to several AWGs with single inputs. In this paper, we present the design and simulation results of an AWG in the H band with 3 input waveguides (channels). The resolving power of individual input channels is 1500, while the overall resolving power with three inputs together is 500, 600, 750 in three different configurations simulated here. The free spectral range of the device is 9.5 nm around a central wavelength of 1600 nm. For the standard multi-input AWG, the relative shift between the output spectra due to adjacent input channels is about 1.6 nm, which roughly equals one spectral channel spacing. In this paper, we discuss ways to increase the resolving power and the number of inputs without compromising the free spectral range or throughput.


title: A Do-it-yourself Spectrograph Kit for Educational Outreach in Optics and Photonics

abstract: We designed and built a do-it-yourself spectrograph assembly to demonstrate the concept of spectroscopy, an indispensable tool for exploring the cosmos. This spectrograph is designed for optical band (400-750 nm). It uses a transmission grating to disperse the light and a webcam to measure the spectrum. This spectrograph provides a resolving power ($\lambda/\delta\lambda$) of about 1000. This demonstration involves off-the-shelf materials costing less than \$500, thus making it an easy to build demonstration kit for a school or public setting. The kit is well-suited for performing various science experiments and acquiring hands-on experience for students to learn the concepts such as coherence, spectral orders, resolving power, absorption and emission spectra. All of these concepts are an integral part of modern astronomical observations as well as various other fields in STEM such as biomedical engineering, chemical analysis, food and water quality, etc. This kit is portable and fully modular, making it apt for outreach purposes.


title: Evidence for UHECR origin in starburst galaxies

abstract: The quest for the origin(s) of ultra-high-energy cosmic rays (UHECRs) continues to be a far-reaching pillar of high energy astrophysics. The source scrutiny is mostly based on three observables: the energy spectrum, the nuclear composition, and the distribution of arrival directions. We show that each of these three observables can be well reproduced with UHECRs originating in starburst galaxies.


title: Composite 2HDM with singlets: a viable dark matter scenario

abstract: We study the non-minimal composite Higgs model with global symmetry $\text{SO}(7)$ broken to $\text{SO}(5)\times \text{SO}(2)$. The model results in a composite Two-Higgs doublet model (2HDM) equipped with two extra singlets, the lightest of which can be a viable dark matter candidate. The model is able to reproduce the correct dark matter relic density both via the usual thermal freeze-out and through late time decay of the heavier singlet. In the case of thermal freeze-out, it is possible to evade current experimental constraints even with the minimum fine tuning allowed by electroweak precision tests.


title: The geometry of graded cotangent bundles

abstract: Given a vector bundle $A\to M$ we study the geometry of the graded manifolds $T^*[k]A[1]$, including their canonical symplectic structures, compatible Q-structures and Lagrangian Q-submanifolds. We relate these graded objects to classical structures, such as higher Courant algebroids on $A\oplus\bigwedge^{k-1}A^*$ and higher Dirac structures therein, semi-direct products of Lie algebroid structures on $A$ with their coadjoint representations up to homotopy, and branes on certain AKSZ $\sigma$-models.


title: Largest Inscribed Rectangles in Geometric Convex Sets

abstract: This paper considers the problem of finding maximum volume (axis-aligned) inscribed boxes in a compact convex set, defined by a finite number of convex inequalities, and presents optimization and geometric approaches for solving them. Several optimization models are developed that can be easily generalized to find other inscribed geometric shapes such as triangles, rhombi, and squares. To find the largest axis-aligned inscribed rectangles in the higher dimensions, an interior-point method algorithm is presented and analyzed. For 2-dimensional space, a parametrized optimization approach is developed to find the largest (axis-aligned) inscribed rectangles in convex sets. The optimization approach provides a uniform framework for solving a wide variety of relevant problems. Finally, two computational geometric $(1-\varepsilon)$--approximation algorithms with sub-linear running times are presented that improve the previous results.


title: Universality of the critical point mapping between Ising model and QCD at small quark mass

abstract: The universality of the QCD equation of state near the critical point is expressed by mapping pressure as a function of temperature $T$ and baryon chemical potential $\mu$ in QCD to Gibbs free energy as a function of reduced temperature $r$ and magnetic field $h$ in the Ising model. The mapping parameters are, in general, not universal, i.e., determined by details of the microscopic dynamics, rather than by symmetries and long-distance dynamics. In this paper we point out that in the limit of small quark masses, when the critical point is close to the tricritical point, the mapping parameters show universal dependence on the quark mass $m_q$. In particular, the angle between the $r=0$ and $h=0$ lines in the $(\mu,T)$ plane vanishes as $m_q^{2/5}$. We discuss possible phenomenological consequences of these findings.

categories: hep-ph cond-mat.stat-mech hep-lat nucl-th


title: The hidden assumptions of Frauchiger and Renner

abstract: This note is a critical examination of the argument of Frauchiger and Renner (Nature Communications 9:3711 (2018)), in which they claim to show that three reasonable assumptions about the use of quantum mechanics jointly lead to a contradiction. It is shown that further assumptions are needed to establish the contradiction, and that each of these assumptions is invalid in some version of quantum mechanics.


title: ALMA unveils widespread molecular gas clumps in the ram pressure stripped tail of the Norma jellyfish galaxy

abstract: We present the first high-resolution map of the cold molecular gas distribution, as traced by CO(2-1) emission with ALMA, in a prominent ram pressure stripped tail. The Norma cluster galaxy ESO 137-001 is undergoing a strong interaction with the surrounding intra-cluster medium and is one of the nearest jellyfish galaxies with a long multi-phase tail. We have mapped the full extent of the tail at 1" (350 pc) angular resolution and found a rich distribution of mostly compact CO regions extending to nearly 60 kpc in length and 25 kpc in width. In total about 10^9 M_sun of molecular gas was detected. The CO features are found predominantly at the heads of numerous small-scale (~ 1.5 kpc) fireballs (i.e., star-forming clouds with linear streams of young stars extending toward the galaxy) but also of large-scale (~ 8 kpc) super-fireballs, and double-sided fireballs that have additional diffuse ionized gas tails extending in the direction opposite to the stellar tails. The new data help to shed light on the origin of the molecular tail - CO filaments oriented in the direction of the tail with only diffuse associated Halpha emission are likely young molecular features formed in situ, whereas other large CO features tilted with respect to the tail may have originated from the densest gas complexes that were pushed gradually away from the disk. The ALMA observations of ESO 137-001, together with observations from HST, Chandra and VLT/MUSE, offer the most complete view of a spectacular ram pressure stripped tail to date.


title: A New Two-Molecule Combination Band as Diagnostic of Carbon Monoxide Diluted in Nitrogen Ice On Triton

abstract: A combination band due to a mechanism whereby a photon excites two or more vibrational modes ({\it e.g.} a bend and a stretch) of an individual molecule is commonly seen in laboratory and astronomical spectroscopy. Here, we present evidence of a much less commonly seen combination band $-$ one where a photon simultaneously excites two adjacent molecules in an ice. In particular, we present near-infrared spectra of laboratory CO/N$_2$ ice samples where we identify a band at 4467.5 cm$^{-1}$ (2.239 $\mu$m) that results from single photons exciting adjacent pairs of CO and N$_2$ molecules. We also present a near-infrared spectrum of Neptune's largest satellite Triton taken with the Gemini-South 8.1 meter telescope and the Immersion Grating Infrared Spectrograph (IGRINS) that shows this 4467.5 cm$^{-1}$ (2.239 $\mu$m) CO-N$_2$ combination band. The existence of the band in a spectrum of Triton indicates that CO and N$_2$ molecules are intimately mixed in the ice rather than existing as separate regions of pure CO and pure N$_2$ deposits. Our finding is important because CO and N$_2$ are the most volatile species on Triton and so dominate seasonal volatile transport across its surface. Our result will place constraints on the interaction between the surface and atmosphere of Triton. 1


title: Clustered Gaussian Graphical Model via Symmetric Convex Clustering

abstract: Knowledge of functional groupings of neurons can shed light on structures of neural circuits and is valuable in many types of neuroimaging studies. However, accurately determining which neurons carry out similar neurological tasks via controlled experiments is both labor-intensive and prohibitively expensive on a large scale. Thus, it is of great interest to cluster neurons that have similar connectivity profiles into functionally coherent groups in a data-driven manner. In this work, we propose the clustered Gaussian graphical model (GGM) and a novel symmetric convex clustering penalty in an unified convex optimization framework for inferring functional clusters among neurons from neural activity data. A parallelizable multi-block Alternating Direction Method of Multipliers (ADMM) algorithm is used to solve the corresponding convex optimization problem. In addition, we establish convergence guarantees for the proposed ADMM algorithm. Experimental results on both synthetic data and real-world neuroscientific data demonstrate the effectiveness of our approach.


title: Radio Frequency Magnet-free Circulators Based on Spatiotemporal Modulation of Surface Acoustic Wave Filters

abstract: In this paper, a new generation of magnet-free circulators with high performance is proposed. Circulators are crucial devices in modern communication systems due to their ability to enable full-duplexing and double the spectral efficiency directly in the physical layer of the radio-frequency (RF) front-end. Traditionally, Lorentz reciprocity is broken by applying magnetic bias to ferrite materials, therefore conventional circulators are bulky and expensive. In this paper, this problem is addressed by replacing the magnetic bias with periodic spatiotemporal modulation. Compared to previous works, the proposed circulator is constructed using surface acoustic wave (SAW) filters instead of transmission lines (TL), which reduces the modulation frequency by at least a factor of 20 and ensures ultra-low power consumption and high linearity. The miniaturized high quality (Q) factor SAW filters also lead to a low-loss non-reciprocal band with strong isolation (IX) and broad bandwidth (BW) on a chip scale, therefore addressing such limitations in previous magnet-free demonstrations. Furthermore, compared to the conventional differential circuit configuration, a novel quad configuration is developed, which doubles the intermodulation-free bandwidth.


title: Interface moment dynamics and its contribution to spin-transfer torque switching process in magnetic tunnel junctions

abstract: A practical problem for memory applications involving perpendicularly magnetized magnetic tunnel junctions is the reliability of switching characteristics at high-bias voltage. Often it has been observed that at high-bias, additional error processes are present that cause a decrease in switching probability upon further increase of bias voltage. We identify the main cause of such error-rise process through examination of switching statistics as a function of bias voltage and applied field, and the junction switching dynamics in real time. These experiments show a coincidental onset of error-rise and the presence of a new low-frequency microwave emission well below that dictated by the anisotropy field. We show that in a few-macrospin coupled numerical model, this is consistent with an interface region with concentrated perpendicular anisotropy, and where the magnetic moment has limited exchange coupling to the rest of the layers. These results point to the important role high-frequency interface magnetic moment dynamics play in determining the switching characteristics of these tunnel junction devices.


title: SDN-based In-network Honeypot: Preemptively Disrupt and Mislead Attacks in IoT Networks

abstract: Detecting cyber attacks in the network environments used by Internet-of-things (IoT) and preventing them from causing physical perturbations play an important role in delivering dependable services. To achieve this goal, we propose in-network Honeypot based on Software-Defined Networking (SDN) to disrupt and mislead adversaries into exposures while they are in an early stage of preparing an attack. Different from traditional Honeypot requiring dedicated hardware setup, the in-network Honeypot directly reroutes traffic from suspicious nodes and intelligently spoofs the network traffic to them by adding misleading information into normal traffic. Preliminary evaluations on real networks demonstrate that the in-network Honeypot can have little impacts on the performance of IoT networks.


title: Study of counterintuitive transport properties in the Aubry-Andr\'e-Harper model via entanglement entropy and persistent current

abstract: The single particle eigenstates of the Aubry-Andr\'e-Harper model are known to show a delocalization-localization transition at a finite strength of the quasi-periodic disorder. In this work, we point out that an intimate relationship exists between the sub-band structure of the spectrum and transport properties of the model. To capture the transport properties we have not only used a variety of single-particle measures like inverse participation ratio, and von Neumann entropy, but also many-particle measures such as persistent current and its variance, and many body entanglement entropy. The many-particle measures are very sensitive to the sub-band structure of the spectrum. Even in the delocalized phase, surprisingly the entanglement entropy is substantially suppressed when the Fermi level is in the band gaps whereas the persistent current is vanishingly small for the same locations of the Fermi level. The entanglement entropy seems to follow area-law exclusively for these special locations of Fermi level or filling fractions of free fermions. A study of the standard deviation of persistent current offers further distinguishing features for the special fillings. In the delocalized phase, the standard deviation vs. mean persistent current curves are discontinuous for the non-special values of filling fractions and continuous (closed) for the special values of filling fractions whereas in the localized phase, these curves become straight lines for both types of filling fractions. We have also discussed how the results depend on the system size. Our results, specially on the persistent current, can potentially be tested experimentally using the present day set-ups based on ultra-cold atoms.


title: Rigidity in dynamics and M\"obius disjointness

abstract: Let $(X, T)$ be a topological dynamical system. We show that if each invariant measure of $(X, T)$ gives rise to a measure-theoretic dynamical system that is either: a. rigid along a sequence of "bounded prime volume" or b. admits a polynomial rate of rigidity on a linearly dense subset in $C(X)$, then $(X, T)$ satisfies Sarnak's conjecture on M\"obius disjointness. We show that the same conclusion also holds if there are countably many invariant ergodic measures, and each of them satisfies a. or b. This recovers some earlier results and implies Sarnak's conjecture in the following new cases: for almost every interval exchange map of $d$ intervals with $d \geq 2$, for $C^{2+\epsilon}$-smooth skew products over rotations and $C^{2+\epsilon}$-smooth flows (without fixed points) on the torus. In particular, these are improvements of earlier results of respectively Chaika-Eskin, Wang and Huang-Wang-Ye. We also discuss some purely arithmetic consequences for the Liouville function.


title: Pass-efficient methods for compression of high-dimensional turbulent flow data

abstract: The future of high-performance computing, specifically on future Exascale computers, will presumably see memory capacity and bandwidth fail to keep pace with data generated, for instance, from massively parallel partial differential equation (PDE) systems. Current strategies proposed to address this bottleneck entail the omission of large fractions of data, as well as the incorporation of $\textit{in situ}$ compression algorithms to avoid overuse of memory. To ensure that post-processing operations are successful, this must be done in a way that a sufficiently accurate representation of the solution is stored. Moreover, in situations where the input/output system becomes a bottleneck in analysis, visualization, etc., or the execution of the PDE solver is expensive, the the number of passes made over the data must be minimized. In the interest of addressing this problem, this work focuses on the utility of pass-efficient, parallelizable, low-rank, matrix decomposition methods in compressing high-dimensional simulation data from turbulent flows. A particular emphasis is placed on using coarse representation of the data -- compatible with the PDE discretization grid -- to accelerate the construction of the low-rank factorization. This includes the presentation of a novel single-pass matrix decomposition algorithm for computing the so-called interpolative decomposition. The methods are described extensively and numerical experiments on two turbulent channel flow data are performed. In the first (unladen) channel flow case, compression factors exceeding $400$ are achieved while maintaining accuracy with respect to first- and second-order flow statistics. In the particle-laden case, compression factors of 100 are achieved and the compressed data is used to recover particle velocities.


title: A survey of advances in vision-based vehicle re-identification

abstract: Vehicle re-identification (V-reID) has become significantly popular in the community due to its applications and research significance. In particular, the V-reID is an important problem that still faces numerous open challenges. This paper reviews different V-reID methods including sensor based methods, hybrid methods, and vision based methods which are further categorized into hand-crafted feature based methods and deep feature based methods. The vision based methods make the V-reID problem particularly interesting, and our review systematically addresses and evaluates these methods for the first time. We conduct experiments on four comprehensive benchmark datasets and compare the performances of recent hand-crafted feature based methods and deep feature based methods. We present the detail analysis of these methods in terms of mean average precision (mAP) and cumulative matching curve (CMC). These analyses provide objective insight into the strengths and weaknesses of these methods. We also provide the details of different V-reID datasets and critically discuss the challenges and future trends of V-reID methods.


title: L\'evy bridges with random length

abstract: In this paper our first goal is to give precise definition of the L\'evy bridges with random length. Our second task is to establish the Markov property of this process with respect to its completed natural filtration and thus with respect to the usual augmentation of this one. This property will be crucial for the right-continuity of completed natural filtration.


title: Large Scale Incremental Learning

abstract: Modern machine learning suffers from catastrophic forgetting when learning new classes incrementally. The performance dramatically degrades due to the missing data of old classes. Incremental learning methods have been proposed to retain the knowledge acquired from the old classes, by using knowledge distilling and keeping a few exemplars from the old classes. However, these methods struggle to scale up to a large number of classes. We believe this is because of the combination of two factors: (a) the data imbalance between the old and new classes, and (b) the increasing number of visually similar classes. Distinguishing between an increasing number of visually similar classes is particularly challenging, when the training data is unbalanced. We propose a simple and effective method to address this data imbalance issue. We found that the last fully connected layer has a strong bias towards the new classes, and this bias can be corrected by a linear model. With two bias parameters, our method performs remarkably well on two large datasets: ImageNet (1000 classes) and MS-Celeb-1M (10000 classes), outperforming the state-of-the-art algorithms by 11.1% and 13.2% respectively.


title: Effects of coplanar satellite bands on galactic disc evolution

abstract: Small dwarf companions have been long thought to have minimal influence on their host galaxy's evolution without undergoing direct impacts to the host's disc. However, in light of recent discoveries of coplanar, corotating satellite structures around the Milky Way, Andromeda, and Centaurus A, we use an N-body/test particle simulation to show that low-mass dwarf satellites within such structures are able to exert significant influence on their host's disc, driving spiral waves and inducing stellar scattering. This is accomplished through quasi-periodic alignments of multiple small satellites within the structure that emulate the gravitational influence of a single, larger satellite such as Sagittarius Dwarf or the Large Magellanic Cloud. We find that the coplanar, corotational nature of such structures allows for repeated alignments on short enough timescales to overcome damping within the disc, and in a consistent enough fashion to continually drive spiral waves over the course of 2 Gyr of simulation time. The spirals driven by this phenomenon tend to be flocculent and many-armed due to the irregular intervals over which alignments occur. We additionally find that while the aligned satellites are able to induce noticeable thickening of the disc, their ability to drive surface density profile evolution is secondary to other effects.


title: Sub-nanosecond switching in a cryogenic spin-torque spin-valve memory element with a dilute permalloy free layer

abstract: We present a study of the pulsed current switching characteristics of spin-valve nanopillars with in-plane magnetized dilute permalloy and undiluted permalloy free layers in the ballistic regime at low temperature. The dilute permalloy free layer device switches much faster: the characteristic switching time for a permalloy free (Ni0.83Fe0.17) layer device is 1.18 ns, while that for a dilute permalloy ([Ni0.83Fe0.17]0.6Cu0.4) free layer device is 0.475 ns. A ballistic macrospin model can capture the data trends with a reduced spin torque asymmetry parameter, reduced spin polarization and increased Gilbert damping for the dilute permalloy free layer relative to the permalloy devices. Our study demonstrates that reducing the magnetization of the free layer increases the switching speed while greatly reducing the switching energy and shows a promising route toward even lower power magnetic memory devices compatible with superconducting electronics.


title: Singularity formation in fractional Burgers' equations

abstract: The formation of singularities in finite time in non-local Burgers' equations, with time-fractional derivative, is studied in detail. The occurrence of finite time singularity is proved, revealing the underlying mechanism, and precise estimates on the blow-up time are provided. The employment of the present equation to model a problem arising in job market is also analyzed.


title: Using Metrics Suites to Improve the Measurement of Privacy in Graphs

abstract: Social graphs are widely used in research (e.g., epidemiology) and business (e.g., recommender systems). However, sharing these graphs poses privacy risks because they contain sensitive information about individuals. Graph anonymization techniques aim to protect individual users in a graph, while graph de-anonymization aims to re-identify users. The effectiveness of anonymization and de-anonymization algorithms is usually evaluated with privacy metrics. However, it is unclear how strong existing privacy metrics are when they are used in graph privacy. In this paper, we study 26 privacy metrics for graph anonymization and de-anonymization and evaluate their strength in terms of three criteria: monotonicity indicates whether the metric indicates lower privacy for stronger adversaries; for within-scenario comparisons, evenness indicates whether metric values are spread evenly; and for between-scenario comparisons, shared value range indicates whether metrics use a consistent value range across scenarios. Our extensive experiments indicate that no single metric fulfills all three criteria perfectly. We therefore use methods from multi-criteria decision analysis to aggregate multiple metrics in a metrics suite, and we show that these metrics suites improve monotonicity compared to the best individual metric. This important result enables more monotonic, and thus more accurate, evaluations of new graph anonymization and de-anonymization algorithms.


title: Ternary mappings of triangular algebras

abstract: We take a categorical approach to describe ternary derivations and ternary automorphisms of triangular algebras. New classes of automorphisms and derivations of triangular algebras are also introduced and studied.


title: Epsilon-Lexicase Selection for Regression

abstract: Lexicase selection is a parent selection method that considers test cases separately, rather than in aggregate, when performing parent selection. It performs well in discrete error spaces but not on the continuous-valued problems that compose most system identification tasks. In this paper, we develop a new form of lexicase selection for symbolic regression, named epsilon-lexicase selection, that redefines the pass condition for individuals on each test case in a more effective way. We run a series of experiments on real-world and synthetic problems with several treatments of epsilon and quantify how epsilon affects parent selection and model performance. epsilon-lexicase selection is shown to be effective for regression, producing better fit models compared to other techniques such as tournament selection and age-fitness Pareto optimization. We demonstrate that epsilon can be adapted automatically for individual test cases based on the population performance distribution. Our experiments show that epsilon-lexicase selection with automatic epsilon produces the most accurate models across tested problems with negligible computational overhead. We show that behavioral diversity is exceptionally high in lexicase selection treatments, and that epsilon-lexicase selection makes use of more fitness cases when selecting parents than lexicase selection, which helps explain the performance improvement.


title: Learning Nearest Neighbor Graphs from Noisy Distance Samples

abstract: We consider the problem of learning the nearest neighbor graph of a dataset of n items. The metric is unknown, but we can query an oracle to obtain a noisy estimate of the distance between any pair of items. This framework applies to problem domains where one wants to learn people's preferences from responses commonly modeled as noisy distance judgments. In this paper, we propose an active algorithm to find the graph with high probability and analyze its query complexity. In contrast to existing work that forces Euclidean structure, our method is valid for general metrics, assuming only symmetry and the triangle inequality. Furthermore, we demonstrate efficiency of our method empirically and theoretically, needing only O(n log(n)Delta^-2) queries in favorable settings, where Delta^-2 accounts for the effect of noise. Using crowd-sourced data collected for a subset of the UT Zappos50K dataset, we apply our algorithm to learn which shoes people believe are most similar and show that it beats both an active baseline and ordinal embedding.


title: Interpretable PID Parameter Tuning for Control Engineering using General Dynamic Neural Networks: An Extensive Comparison

abstract: Modern automation systems rely on closed loop control, wherein a controller interacts with a controlled process, based on observations. These systems are increasingly complex, yet most controllers are linear Proportional-Integral-Derivative (PID) controllers. PID controllers perform well on linear and near-linear systems but their simplicity is at odds with the robustness required to reliably control complex processes. Modern machine learning offers a way to extend PID controllers beyond their linear capabilities by using neural networks. However, such an extension comes at the cost of losing stability guarantees and controller interpretability. In this paper, we examine the utility of extending PID controllers with recurrent neural networks-namely, General Dynamic Neural Networks (GDNN); we show that GDNN (neural) PID controllers perform well on a range of control systems and highlight how they can be a scalable and interpretable option for control systems. To do so, we provide an extensive study using four benchmark systems that represent the most common control engineering benchmarks. All control benchmarks are evaluated with and without noise as well as with and without disturbances. The neural PID controller performs better than standard PID control in 15 of 16 tasks and better than model-based control in 13 of 16 tasks. As a second contribution, we address the lack of interpretability that prevents neural networks from being used in real-world control processes. We use bounded-input bounded-output stability analysis to evaluate the parameters suggested by the neural network, thus making them understandable. This combination of rigorous evaluation paired with better interpretability is an important step towards the acceptance of neural-network-based control approaches. It is furthermore an important step towards interpretable and safely applied artificial intelligence.


title: Quantum matrix diagonalization visualized

abstract: We show how to visualize the process of diagonalizing the Hamiltonian matrix to find the energy eigenvalues and eigenvectors of a generic one-dimensional quantum system. Starting in the familiar sine-wave basis of an embedding infinite square well, we display the Hamiltonian matrix graphically with the basis functions alongside. Each step in the diagonalization process consists of selecting a nonzero off-diagonal matrix element, then rotating the two corresponding basis vectors in their own subspace until this element is zero. We provide Mathematica code to display the effects of these rotations on both the matrix and the basis functions. As an electronic supplement we also provide a JavaScript web app to interactively carry out this process.


title: Global smooth solutions with large data for a system modeling aurora type phenomena in the 2-torus

abstract: We prove existence and uniqueness of smooth solutions with large initial data for a system of equations modeling the interaction of short waves, governed by a nonlinear Schr\"odinger equation, and long waves, described by the equations of magnetohydrodynamics. In the model, the short waves propagate along the streamlines of the fluid flow. This is translated in the system by setting up the nonlinear Schr\"odinger equation in the Lagrangian coordinates of the fluid. Besides, the equations are coupled by nonlinear terms accounting for the strong interaction of the dynamics. The system provides a simplified mathematical model for studying aurora type phenomena. We focus on the 2-dimensional case with periodic boundary conditions. This is the first result on existence of smooth solutions with large data for the multidimensional case of the model under consideration.


title: Exploring Computational User Models for Agent Policy Summarization

abstract: AI agents are being developed to support high stakes decision-making processes from driving cars to prescribing drugs, making it increasingly important for human users to understand their behavior. Policy summarization methods aim to convey strengths and weaknesses of such agents by demonstrating their behavior in a subset of informative states. Some policy summarization methods extract a summary that optimizes the ability to reconstruct the agent's policy under the assumption that users will deploy inverse reinforcement learning. In this paper, we explore the use of different models for extracting summaries. We introduce an imitation learning-based approach to policy summarization; we demonstrate through computational simulations that a mismatch between the model used to extract a summary and the model used to reconstruct the policy results in worse reconstruction quality; and we demonstrate through a human-subject study that people use different models to reconstruct policies in different contexts, and that matching the summary extraction model to these can improve performance. Together, our results suggest that it is important to carefully consider user models in policy summarization.