title: Improved memory in recurrent neural networks with sequential non-normal dynamics

abstract: Training recurrent neural networks (RNNs) is a hard problem due to degeneracies in the optimization landscape, a problem also known as vanishing/exploding gradients. Short of designing new RNN architectures, previous methods for dealing with this problem usually boil down to orthogonalization of the recurrent dynamics, either at initialization or during the entire training period. The basic motivation behind these methods is that orthogonal transformations are isometries of the Euclidean space, hence they preserve (Euclidean) norms and effectively deal with vanishing/exploding gradients. However, this ignores the crucial effects of non-linearity and noise. In the presence of a non-linearity, orthogonal transformations no longer preserve norms, suggesting that alternative transformations might be better suited to non-linear networks. Moreover, in the presence of noise, norm preservation itself ceases to be the ideal objective. A more sensible objective is maximizing the signal-to-noise ratio (SNR) of the propagated signal instead. Previous work has shown that in the linear case, recurrent networks that maximize the SNR display strongly non-normal, sequential dynamics and orthogonal networks are highly suboptimal by this measure. Motivated by this finding, here we investigate the potential of non-normal RNNs, i.e. RNNs with a non-normal recurrent connectivity matrix, in sequential processing tasks. Our experimental results show that non-normal RNNs outperform their orthogonal counterparts in a diverse range of benchmarks. We also find evidence for increased non-normality and hidden chain-like feedforward motifs in trained RNNs initialized with orthogonal recurrent connectivity matrices.


title: Reference Capabilities for Safe Parallel Array Programming

abstract: The array is a fundamental data structure that provides an efficient way to store and retrieve non-sparse data contiguous in memory. Arrays are important for the performance of many memory-intensive applications due to the design of modern memory hierarchies: contiguous storage facilitates spatial locality and predictive access patterns which enables prefetching. Operations on large arrays often lend themselves well to parallelisation, such as a fork-join style divide-and-conquer algorithm for sorting. For parallel operations on arrays to be deterministic, data-race freedom must be guaranteed. For operations on arrays of primitive data, data-race freedom is obtained by coordinating accesses so that no two threads operate on the same array indices. This is however not enough for arrays of non-primitives due to aliasing: accesses of separate array elements may return pointers to the same object, or overlapping structures. Reference capabilities have been used successfully in the past to statically guarantee the absence of data-races in object-oriented programs. This paper presents the first extension of reference capabilities -- called array capabilities -- that support concurrent and parallel operations on arrays of both primitive and non-primitive values. In addition to element access, array capabilities support the abstract manipulation of arrays, logical splitting of arrays into subarrays, and merging subarrays. These operations allow expressing a wide range of array use cases. This paper presents the array capability design space and show how it applies to a number of array use cases. The core ideas are formalised and proven sound in a simple calculus, along with a proof that shows that well-typed programs with array capabilities are free from data-races.


title: Behavior of quantum correlations under nondissipative decoherence by means of the correlation matrix

abstract: In this paper we use the Fano representation of two-qubit states from which we can identify a correlation matrix containing the information about the classical and quantum correlations present in the bipartite quantum state. To illustrate the use of this matrix, we analyze the behavior of the correlations under non-dissipative decoherence in two-qubit states with maximally mixed marginals. From the behavior of the elements of the correlation matrix before and after making measurements on one of the subsystems, we identify the classical and quantum correlations present in the Bell-diagonal states. In addition, we use the correlation matrix to study the phenomenon known as freezing of quantum discord. We find that under some initial conditions where freezing of quantum discord takes place, quantum correlation instead may remain not constant. In order to further explore into these results we also compute a non-commutativity measure of quantum correlations to analyze the behavior of quantum correlations under non-dissipative decoherence. We conclude from our study that freezing of quantum discord may not always be identified as equivalent to the freezing of the actual quantum correlations.


title: Periodic projections of alternating knots

abstract: This paper is devoted to prove the existence of $q$-periodic alternating projections of prime alternating $q$-periodic knots. The main tool is the Menasco-Thistlethwaite's Flyping theorem. Let $K$ be an oriented prime alternating knot that is $q$-periodic with $q\geq 3$, i.e. $K$ admits a symmetry that is a rotation of order $q$. Then $K$ has an alternating $q$-periodic projection. As applications, we obtain the crossing number of a $q$ -periodic alternating knot with $q\geq 3$ is a multiple of $q$ and we give an elementary proof that the knot $12_{a634} $ is not 3-periodic; this proof does not depend on computer computations as in "Periodic knots and Heegaard Floer correction terms" by Stanilav Jabuka and Swatee Naik (arXiv:1307.5116 [math.GT]).


title: Reinforcement Learning Experience Reuse with Policy Residual Representation

abstract: Experience reuse is key to sample-efficient reinforcement learning. One of the critical issues is how the experience is represented and stored. Previously, the experience can be stored in the forms of features, individual models, and the average model, each lying at a different granularity. However, new tasks may require experience across multiple granularities. In this paper, we propose the policy residual representation (PRR) network, which can extract and store multiple levels of experience. PRR network is trained on a set of tasks with a multi-level architecture, where a module in each level corresponds to a subset of the tasks. Therefore, the PRR network represents the experience in a spectrum-like way. When training on a new task, PRR can provide different levels of experience for accelerating the learning. We experiment with the PRR network on a set of grid world navigation tasks, locomotion tasks, and fighting tasks in a video game. The results show that the PRR network leads to better reuse of experience and thus outperforms some state-of-the-art approaches.


title: Review on tests of General Relativity And Modified Gravity Using Pulsar Timing

abstract: In this review we aim to describe the most stringent tests of the strong equivalence principle, the fundamental principle of General Relativity, using pulsar timing. For this purpose, we first construct the parametrized post-Newtonian and post-Keplerian formalism together with their parameters. Then we constrain the post-Newtonian parameters associated to the violation of the strong equivalence principle with the most stringent tests using pulsar timing techniques. We will in particular see that these bounds are the most constraining limits on a violation of the principle. Finally we will discuss the implication of these results on scalar-tensor theories as well as massive gravity. We will find that the set of possible scalar-tensor theories can be tightly reduce, however in the case of massive gravity, the level of precision is not sufficient to probe the predicted mass of the graviton.


title: A generalization of analytic torsion via differential forms on spaces of metrics

abstract: We introduce multi-torsion, a spectral invariant generalizing Ray-Singer analytic torsion. We define multi-torsion for compact manifolds with a certain local geometric product structure that gives a bigrading on differential forms. We prove that multi-torsion is metric-independent in a suitable sense. Our definition of multi-torsion is inspired by an interpretation of each of analytic torsion and the eta invariant as a regularized integral of a closed differential form on a space of metrics on a vector bundle or on a space of elliptic operators. We generalize the Stokes' theorem argument explaining the dependence of torsion and eta on the geometric data used to define them to the local product setting to prove our metric-independence theorem for multi-torsion.


title: On approximate solutions of the equations of incompressible magnetohydrodynamics

abstract: Inspired by an approach proposed previously for the incompressible Navier-Stokes (NS) equations, we present a general framework for the a posteriori analysis of the equations of incompressible magnetohydrodynamics (MHD) on a torus of arbitrary dimension d; this setting involves a Sobolev space of infinite order, made of C^infinity vector fields (with vanishing divergence and mean) on the torus. Given any approximate solution of the MHD Cauchy problem, its a posteriori analysis with the method of the present work allows to infer a lower bound on the time of existence of the exact solution, and to bound from above the Sobolev distance of any order between the exact and the approximate solution. In certain cases the above mentioned lower bound on the time of existence is found to be infinite, so one infers the global existence of the exact MHD solution. We present some applications of this general scheme; the most sophisticated one lives in dimension d=3, with the ABC flow (perturbed magnetically) as an initial datum, and uses for the Cauchy problem a Galerkin approximate solution in 124 Fourier modes. We illustrate the conclusions arising in this case from the a posteriori analysis of the Galerkin approximant; these include the derivation of global existence of the exact MHD solution with the ABC datum, when the dimensionless viscosity and resistivity are equal and stay above an explicitly given threshold value.


title: Non-relativistic expansion of the Einstein-Hilbert Lagrangian

abstract: We present a systematic technique to expand the Einstein-Hilbert Lagrangian in inverse powers of the speed of light squared. The corresponding result for the non-relativistic gravity Lagrangian is given up to next-to-next-to-leading order. The techniques are universal and can be used to expand any Lagrangian theory whose fields are a function of a given parameter.


title: Singular quasilinear elliptic systems involving gradient terms

abstract: In this paper we establish existence of smooth positive solutions for a singular quasilinear elliptic system involving gradient terms. The approach combines sub-supersolutions method and Schauder's fixed point theorem.


title: Are Labels Required for Improving Adversarial Robustness?

abstract: Recent work has uncovered the interesting (and somewhat surprising) finding that training models to be invariant to adversarial perturbations requires substantially larger datasets than those required for standard classification. This result is a key hurdle in the deployment of robust machine learning models in many real world applications where labeled data is expensive. Our main insight is that unlabeled data can be a competitive alternative to labeled data for training adversarially robust models. Theoretically, we show that in a simple statistical setting, the sample complexity for learning an adversarially robust model from unlabeled data matches the fully supervised case up to constant factors. On standard datasets like CIFAR-10, a simple Unsupervised Adversarial Training (UAT) approach using unlabeled data improves robust accuracy by 21.7% over using 4K supervised examples alone, and captures over 95% of the improvement from the same number of labeled examples. Finally, we report an improvement of 4% over the previous state-of-the-art on CIFAR-10 against the strongest known attack by using additional unlabeled data from the uncurated 80 Million Tiny Images dataset. This demonstrates that our finding extends as well to the more realistic case where unlabeled data is also uncurated, therefore opening a new avenue for improving adversarial training.


title: Minimal submanifolds from the abelian Higgs model

abstract: Given a Hermitian line bundle $L\to M$ over a closed, oriented Riemannian manifold $M$, we study the asymptotic behavior, as $\epsilon\to 0$, of couples $(u_\epsilon,\nabla_\epsilon)$ critical for the rescalings \begin{align*} &E_\epsilon(u,\nabla)=\int_M\Big(|\nabla u|^2+\epsilon^2|F_\nabla|^2+\frac{1}{4\epsilon^2}(1-|u|^2)^2\Big) \end{align*} of the self-dual Yang-Mills-Higgs energy, where $u$ is a section of $L$ and $\nabla$ is a Hermitian connection on $L$ with curvature $F_{\nabla}$. Under the natural assumption $\limsup_{\epsilon\to 0}E_\epsilon(u_\epsilon,\nabla_\epsilon)<\infty$, we show that the energy measures converge subsequentially to (the weight measure $\mu$ of) a stationary integral $(n-2)$-varifold. Also, we show that the $(n-2)$-currents dual to the curvature forms converge subsequentially to $2\pi\Gamma$, for an integral $(n-2)$-cycle $\Gamma$ with $|\Gamma|\le\mu$. Finally, we provide a variational construction of nontrivial critical points $(u_\epsilon,\nabla_\epsilon)$ on arbitrary line bundles, satisfying a uniform energy bound. As a byproduct, we obtain a PDE proof, in codimension two, of Almgren's existence result of (nontrivial) stationary integral $(n-2)$-varifolds in an arbitrary closed Riemannian manifold.


title: PowerSGD: Practical Low-Rank Gradient Compression for Distributed Optimization

abstract: We study gradient compression methods to alleviate the communication bottleneck in data-parallel distributed optimization. Despite the significant attention received, current compression schemes either do not scale well or fail to achieve the target test accuracy. We propose a new low-rank gradient compressor based on power iteration that can i) compress gradients rapidly, ii) efficiently aggregate the compressed gradients using all-reduce, and iii) achieve test performance on par with SGD. The proposed algorithm is the only method evaluated that achieves consistent wall-clock speedups when benchmarked against regular SGD with an optimized communication backend. We demonstrate reduced training times for convolutional networks as well as LSTMs on common datasets. Our code is available at https://github.com/epfml/powersgd.


title: Pre-Training Graph Neural Networks for Generic Structural Feature Extraction

abstract: Graph neural networks (GNNs) are shown to be successful in modeling applications with graph structures. However, training an accurate GNN model requires a large collection of labeled data and expressive features, which might be inaccessible for some applications. To tackle this problem, we propose a pre-training framework that captures generic graph structural information that is transferable across tasks. Our framework can leverage the following three tasks: 1) denoising link reconstruction, 2) centrality score ranking, and 3) cluster preserving. The pre-training procedure can be conducted purely on the synthetic graphs, and the pre-trained GNN is then adapted for downstream applications. With the proposed pre-training procedure, the generic structural information is learned and preserved, thus the pre-trained GNN requires less amount of labeled data and fewer domain-specific features to achieve high performance on different downstream tasks. Comprehensive experiments demonstrate that our proposed framework can significantly enhance the performance of various tasks at the level of node, link, and graph.


title: General solution to the U(1) anomaly equations

abstract: The anomaly cancellation equations for the $U(1)$ gauge group can be written as a cubic equation in $n-1$ integer variables, where $n$ is the number of Weyl fermions carrying the $U(1)$ charge. We solve this Diophantine cubic equation by providing a parametrization of the charges in terms of $n-2$ integers, and prove that this is the most general solution.


title: The pebbling threshold spectrum and paths

abstract: Given a distribution of pebbles on the vertices of a graph, say that we can pebble a vertex if a pebble is left on it after some sequence of moves, each of which takes two pebbles from some vertex and places one on an adjacent vertex. A distribution is solvable if all vertices are pebblable; the pebbling threshold of a sequence of graphs is, roughly speaking, the total number of pebbles for which random distributions with that number of pebbles on a graph in the sequence change from being almost never solvable to being almost always solvable. We show that any sequence of connected graphs with strictly increasing orders always has some pebbling threshold which is $\Omega(\sqrt{n})$ and $O(2^{\sqrt{2 \log_2 n}} n/\sqrt{\log_2 n})$, and that it is possible to construct such a sequence of connected graphs which has any desired pebbling threshold between these bounds. (Here, $n$ is the order of a graph in the sequence.) It follows that the sequence of paths, which, improving earlier estimates, we show has pebbling threshold $\Theta(2^{\sqrt{\log_2 n}} n/\sqrt{\log_2 n})$, does not have the greatest possible pebbling threshold.


title: Accretion into black hole, and formation of magnetically arrested accretion disks

abstract: The exact time-dependent solution is obtained for a magnetic field growth during a spherically symmetric accretion into a black hole (BH) with a Schwarzschild metric. Magnetic field is increasing with time, changing from the initially uniform into a quasi-radial field. Equipartition between magnetic and kinetic energies in the falling gas is established in the developed stages of the flow. Estimates of the synchrotron radiation intensity are presented for the stationary flow. The main part of the radiation is formed in the region $r \leq 7 r_g$, here $r_g$ is a BH gravitational radius. The two-dimensional stationary self-similar magnetohydrodynamic solution is obtained for the matter accretion into BH, in a presence of a large-scale magnetic field, when the magnetic field far from the BH is homogeneous and does not influence the flow. At the symmetry plane perpendicular to the direction of the distant magnetic field, the quasi-stationary disk is formed around BH, which structure is determined by dissipation processes. Parameters of the shock forming due to matter infall onto the disk are obtained. The radiation spectrum of the disk and the shock are obtained for the $10\,\, M_\odot$ BH. The luminosity of such object is about the solar one, for a characteristic galactic gas density, with possibility of observation at distances less than 1 kpc. The spectra of a laminar and a turbulent disk structure around BH are very different. The turbulent disk emits a large part of its flux in the infrared. It may occur that some of the galactic infrared star-like sources are a single BH in the turbulent accretion state. The radiative efficiency of the magnetized disk is very high, reaching $\sim 0.5\,\dot M\,c^2$ so it was called recently as a magnetically arrested disk (MAD). Numerical simulations of MAD, and its appearance during accretion into neutron stars are considered and discussed.


title: End to end learning and optimization on graphs

abstract: Real-world applications often combine learning and optimization problems on graphs. For instance, our objective may be to cluster the graph in order to detect meaningful communities (or solve other common graph optimization problems such as facility location, maxcut, and so on). However, graphs or related attributes are often only partially observed, introducing learning problems such as link prediction which must be solved prior to optimization. Standard approaches treat learning and optimization entirely separately, while recent machine learning work aims to predict the optimal solution directly from the inputs. Here, we propose an alternative decision-focused learning approach that integrates a differentiable proxy for common graph optimization problems as a layer in learned systems. The main idea is to learn a representation that maps the original optimization problem onto a simpler proxy problem that can be efficiently differentiated through. Experimental results show that our ClusterNet system outperforms both pure end-to-end approaches (that directly predict the optimal solution) and standard approaches that entirely separate learning and optimization. Code for our system is available at https://github.com/bwilder0/clusternet.


title: Data assimilation in price formation

abstract: We consider the problem of estimating the density of buyers and vendors in a nonlinear parabolic price formation model using measurements of the price and the transaction rate. Our approach is based on a work by Puel et al., see \cite{Puel2002}, and results in a optimal control problem. We analyse this problems and provide stability estimates for the controls as well as the unknown density in the presence of measurement errors. Our analytic findings are supported with numerical experiments.


title: Disorder-immune photonics based on Mie-resonant dielectric metamaterials

abstract: When the feature size of photonic structures becomes comparable or even smaller than the wavelength of light, the fabrication imperfections inevitably introduce disorder that may eliminate many functionalities of subwavelength photonic devices. Here we suggest a novel concept to achieve a robust bandgap which can endure disorder beyond 30 as a result of the transition from photonic crystals to Mie-resonant metamaterials. By utilizing Mie-resonant metamaterials with high refractive index, we demonstrate photonic waveguides and cavities with strong robustness to position disorder, thus providing a novel approach to the bandgap-based nanophotonic devices with new properties and functionalities.


title: Three-dimensional instability of flame fronts in type I X-ray bursts

abstract: We present the first realistic 3D simulations of flame front instabilities during type I X-ray bursts. The unperturbed front is characterised by the balance between the pressure gradient and the Coriolis force of a spinning neutron star ({\nu} = 450 Hz in our case). This balance leads to a fast horizontal velocity field parallel to the flame front. This flow is strongly sheared in the vertical direction. When we perturb the front an instability quickly corrugates the front. We identify this instability as the baroclinic instability. Most importantly, the flame is not disrupted by the instability and there are two major consequences: the overall flame propagation speed is {\sim} 10 times faster than in the unperturbed case and distinct flame vortices appear. The speedup is due to the corrugation of the front and the dynamics of the vortices. These vortices may also be linked to the oscillations observed in the lightcurves of the bursts.


title: Unlabeled Data Improves Adversarial Robustness

abstract: We demonstrate, theoretically and empirically, that adversarial robustness can significantly benefit from semisupervised learning. Theoretically, we revisit the simple Gaussian model of Schmidt et al. that shows a sample complexity gap between standard and robust classification. We prove that unlabeled data bridges this gap: a simple semisupervised learning procedure (self-training) achieves high robust accuracy using the same number of labels required for achieving high standard accuracy. Empirically, we augment CIFAR-10 with 500K unlabeled images sourced from 80 Million Tiny Images and use robust self-training to outperform state-of-the-art robust accuracies by over 5 points in (i) $\ell_\infty$ robustness against several strong attacks via adversarial training and (ii) certified $\ell_2$ and $\ell_\infty$ robustness via randomized smoothing. On SVHN, adding the dataset's own extra training set with the labels removed provides gains of 4 to 10 points, within 1 point of the gain from using the extra labels.


title: Protocols for Checking Compromised Credentials

abstract: To prevent credential stuffing attacks, industry best practice now proactively checks if user credentials are present in known data breaches. Recently, some web services, such as HaveIBeenPwned (HIBP) and Google Password Checkup (GPC), have started providing APIs to check for breached passwords. We refer to such services as compromised credential checking (C3) services. We give the first formal description of C3 services, detailing different settings and operational requirements, and we give relevant threat models. One key security requirement is the secrecy of a user's passwords that are being checked. Current widely deployed C3 services have the user share a small prefix of a hash computed over the user's password. We provide a framework for empirically analyzing the leakage of such protocols, showing that in some contexts knowing the hash prefixes leads to a 12x increase in the efficacy of remote guessing attacks. We propose two new protocols that provide stronger protection for users' passwords, implement them, and show experimentally that they remain practical to deploy.


title: Fiducial cross sections for the four-lepton decay mode in Higgs-plus-jet production up to NNLO QCD

abstract: The four-lepton decay mode of the Higgs boson allows for a clean kinematic reconstruction, thereby enabling precision studies of the Higgs boson properties and of its production dynamics. We compute the NNLO QCD corrections to fiducial cross sections relevant to this decay mode in the gluon-fusion production of a Higgs boson in association with a hadronic jet, and study the impact of the QCD corrections on the fiducial acceptance factors in inclusive Higgs and Higgs-plus-jet production. We investigate in detail the different definitions used in the ATLAS and CMS measurements to define the fiducial cross sections. Differences in the lepton isolation prescription are found to have a sizeable impact on the higher order corrections to the fiducial acceptance factors.


title: Threshold for blowup for the supercritical cubic wave equation

abstract: We consider the focusing cubic wave equation in the energy supercritical case, i.e., in dimensions $d \geq 5$. For this model an explicit nontrivial self-similar blowup solution was recently found by the first and third author in \cite{GlogicSchoerkhuber}. Furthermore, the solution is proven to be co-dimension one stable in $d=7$. In this paper, we study the equation from a numerical point of view. For $d=5$ and $d=7$ in the radial case, we provide evidence that this solution is at the threshold between generic ODE blowup and dispersion. In addition, we investigate the spectral problem that underlies the stability analysis and compute the spectrum in general supercritical dimensions.

categories: math.AP cs.NA math-ph math.MP math.NA math.SP


title: Budget Minimization with Precedence Constraints

abstract: Budget Minimization is a scheduling problem with precedence constraints, i.e., a scheduling problem on a partially ordered set of jobs $(N, \unlhd)$. A job $j \in N$ is available for scheduling, if all jobs $i \in N$ with $i \unlhd j$ are completed. Further, each job $j \in N$ is assigned real valued costs $c_{j}$, which can be negative or positive. A schedule is an ordering $j_{1}, \dots, j_{\vert N \vert}$ of all jobs in $N$. The budget of a schedule is the external investment needed to complete all jobs, i.e., it is $\max_{l \in \{0, \dots, \vert N \vert \} } \sum_{1 \le k \le l} c_{j_{k}}$. The goal is to find a schedule with minimum budget. Rafiey et al. (2015) showed that Budget Minimization is NP-hard following from a reduction from a molecular folding problem. We extend this result and prove that it is NP-hard to $\alpha(N)$-approximate the minimum budget even on bipartite partial orders. We present structural insights that lead to arguably simpler algorithms and extensions of the results by Rafiey et al. (2015). In particular, we show that there always exists an optimal solution that partitions the set of jobs and schedules each subset independently of the other jobs. We use this structural insight to derive polynomial-time algorithms that solve the problem to optimality on series-parallel and convex bipartite partial orders.


title: Self-Referencing Embedded Strings (SELFIES): A 100% robust molecular string representation

abstract: The discovery of novel materials and functional molecules can help to solve some of society's most urgent challenges, ranging from efficient energy harvesting and storage to uncovering novel pharmaceutical drug candidates. Traditionally matter engineering -- generally denoted as inverse design -- was based massively on human intuition and high-throughput virtual screening. The last few years have seen the emergence of significant interest in computer-inspired designs based on evolutionary or deep learning methods. The major challenge here is that the standard strings molecular representation SMILES shows substantial weaknesses in that task because large fractions of strings do not correspond to valid molecules. Here, we solve this problem at a fundamental level and introduce SELFIES (SELF-referencIng Embedded Strings), a string-based representation of molecules which is 100\% robust. Every SELFIES string corresponds to a valid molecule, and SELFIES can represent every molecule. SELFIES can be directly applied in arbitrary machine learning models without the adaptation of the models; each of the generated molecule candidates is valid. In our experiments, the model's internal memory stores two orders of magnitude more diverse molecules than a similar test with SMILES. Furthermore, as all molecules are valid, it allows for explanation and interpretation of the internal working of the generative models.


title: High Dimensional Classification via Regularized and Unregularized Empirical Risk Minimization: Precise Error and Optimal Loss

abstract: This article provides, through theoretical analysis, an in-depth understanding of the classification performance of the empirical risk minimization framework, in both ridge-regularized and unregularized cases, when high dimensional data are considered. Focusing on the fundamental problem of separating a two-class Gaussian mixture, the proposed analysis allows for a precise prediction of the classification error for a set of numerous data vectors $\mathbf{x} \in \mathbb R^p$ of sufficiently large dimension $p$. This precise error depends on the loss function, the number of training samples, and the statistics of the mixture data model. It is shown to hold beyond Gaussian distribution under some additional non-sparsity condition of the data statistics. Building upon this quantitative error analysis, we identify the simple square loss as the optimal choice for high dimensional classification in both ridge-regularized and unregularized cases, regardless of the number of training samples.


title: Age of Information in G/G/1/1 Systems: Age Expressions, Bounds, Special Cases, and Optimization

abstract: We consider the age of information in G/G/1/1 systems under two service discipline models. In the first model, if a new update arrives when the service is busy, it is blocked; in the second model, a new update preempts the current update in service. For the blocking model, we first derive an exact age expression, then we propose two simple to calculate upper bounds for the average age. The first upper bound assumes the interarrival times to have log-concave distribution. The second upper bound assumes both the interarrivals and service times to have log-concave distribution. Both upper bounds are tight in the case of M/M/1/1 systems. We show that deterministic interarrivals and service times are optimum for the blocking service model. In addition, using the age expression for G/G/1/1 systems, we calculate average age expressions for special cases, i.e., M/G/1/1 and G/M/1/1 systems. Next, for the preemption in service model, we first derive an exact average age expression for G/G/1/1 systems. Then, we propose a simple to calculate upper bound for the average age. In addition, similar to blocking discipline, using the age expression for G/G/1/1 systems, we calculate average age expressions for special cases, i.e., M/G/1/1 and G/M/1/1 systems. Average age for G/M/1/1 can be written as a summation of two terms, the first of which depends only on the first and second moments of interarrival times and the second of which depends only on the service rate. In other words, interarrival and service times are decoupled. We show that deterministic interarrivals are optimum for G/M/1/1 systems. On the other hand, we observe for non-exponential service times that the optimal distribution of interarrival times depends on the relative values of the mean interarrival time and the mean service time.


title: Detecting Light Dark Matter with Magnons

abstract: Scattering of light dark matter with sub-eV energy deposition can be detected with collective excitations in condensed matter systems. When dark matter has spin-independent couplings to atoms or ions, it has been shown to efficiently excite phonons. Here we show that, if dark matter couples to the electron spin, magnon excitations in materials with magnetic dipole order offer a promising detection path. We derive general formulae for single magnon excitation rates from dark matter scattering, and demonstrate as a proof of principle the projected reach of a yttrium iron garnet target for several dark matter models with spin-dependent interactions. This highlights the complementarity of various collective excitations in probing different dark matter interactions.

categories: hep-ph astro-ph.CO cond-mat.mtrl-sci hep-ex


title: The stability-limit conjecture revisited

abstract: The stability-limit conjecture (SLC) proposes that the liquid spinodal of water returns to positive pressure in the supercooled region, and that the apparent divergence of water's thermodynamic response functions as temperature decreases are explained by the approach to this reentrant spinodal. Subsequently, it has been argued that the predictions of the SLC are inconsistent with general thermodynamic principles. Here we reconsider the thermodynamic viability of the SLC by examining a model equation of state for water first studied to clarify the relationship of the SLC to the proposed liquid-liquid phase transition in supercooled water. By demonstrating that a binodal may terminate on a spinodal at a point that is not a critical point, we show that the SLC is thermodynamically permissible in a system that has both a liquid-gas and a liquid-liquid phase transition. We also describe and clarify other unusual thermodynamic behavior that may arise in such a system, particularly that associated with the so-called "critical-point-free" scenario for a liquid-liquid phase transition, which may apply to the case of liquid Si.


title: An Efficient Detection of Malware by Naive Bayes Classifier Using GPGPU

abstract: Due to continuous increase in the number of malware (according to AV-Test institute total ~8 x 10^8 malware are already known, and every day they register ~2.5 x 10^4 malware) and files in the computational devices, it is very important to design a system which not only effectively but can also efficiently detect the new or previously unseen malware to prevent/minimize the damages. Therefore, this paper presents a novel group-wise approach for the efficient detection of malware by parallelizing the classification using the power of GPGPU and shown that by using the Naive Bayes classifier the detection speed-up can be boosted up to 200x. The investigation also shows that the classification time increases significantly with the number of features.


title: A Survey on the Detection of Android Malicious Apps

abstract: Android-based smart devices are exponentially growing, and due to the ubiquity of the Internet, these devices are globally connected to the different devices/networks. Its popularity, attractive features, and mobility make malware creator to put a number of malicious apps in the market to disrupt and annoy the victims. Although to identify the malicious apps, time-to-time various techniques are proposed. However, it appears that malware developers are always ahead of the anti-malware group, and the proposed techniques by the anti-malware groups are not sufficient to counter the advanced malicious apps. Therefore, to understand the various techniques proposed/used for the identification of Android malicious apps, in this paper, we present a survey conducted by us on the work done by the researchers in this field.


title: Maxwell Electrodynamics in terms of Physical Potentials

abstract: A fully relativistically covariant and manifestly gauge invariant formulation of classical Maxwell electrodynamics is presented, purely in terms of gauge invariant potentials without entailing any gauge fixing. We show that the inhomogeneous equations satisfied by the physical scalar and vector potentials (originally discovered by Maxwell) have the same symmetry as the isometry of Minkowski spacetime, thereby reproducing Einstein's incipient approach leading to his discovery of special relativity as a spacetime symmetry. To arrive at this conclusion, we show how the Maxwell equations for the potentials follow from stationary electromagnetism by replacing the Laplacian operator by the d'Alembertian operator, while making all variables dependent on space and time. We also establish consistency of these equations by deriving them from the standard Maxwell equations for the field strengths, showing that there is a unique projection operator which projects onto the physical potentials. Properties of the physical potentials are elaborated through their iterative N\"other coupling to a charged scalar field leading to the Abelian Higgs model, and through a sketch of the Aharonov-Bohm effect, where dependence of the Aharonov-Bohm phase on the physical vector potential is highlighted.


title: Deep Quantitative Susceptibility Mapping for Background Field Removal and Total Field Inversion

abstract: Quantitative susceptibility mapping (QSM) utilizes MRI signal phase to estimate local tissue susceptibility, which has been shown useful to provide novel image contrast and as biomarkers of abnormal tissue. QSM requires addressing a challenging post-processing problem: filtering of image phase estimates and inversion of the phase to susceptibility relationship. A wide variety of quantification errors, robustness limitations, and artifacts constraints QSM clinical translation. To overcome these limitations, a robust deep-learning-based QSM reconstruction approach is proposed to perform background field removal and susceptibility inversion simultaneously from input MRI phase images. Synthetic training data based on in-vivo data sources and physics simulations were used for training. The network was quantitatively tested using gold-standard in-silico labeled dataset against established background field removal and QSM inversion approaches. In addition, the algorithm was applied to a QSM challenge data and clinical susceptibility-weighted imaging (SWI) data. When quantitatively compared against gold-standard in-silico labels, the proposed algorithm outperformed the existing comparable background field removal approaches and QSM reconstruction algorithms. The QSM challenge data and clinical SWI data demonstrated that the proposed approach was able to robustly generate high quality local field and QSM with improved accuracy.


title: Sketch2code: Generating a website from a paper mockup

abstract: An early stage of developing user-facing applications is creating a wireframe to layout the interface. Once a wireframe has been created it is given to a developer to implement in code. Developing boiler plate user interface code is time consuming work but still requires an experienced developer. In this dissertation we present two approaches which automates this process, one using classical computer vision techniques, and another using a novel application of deep semantic segmentation networks. We release a dataset of websites which can be used to train and evaluate these approaches. Further, we have designed a novel evaluation framework which allows empirical evaluation by creating synthetic sketches. Our evaluation illustrates that our deep learning approach outperforms our classical computer vision approach and we conclude that deep learning is the most promising direction for future research.


title: Functional Adversarial Attacks

abstract: We propose functional adversarial attacks, a novel class of threat models for crafting adversarial examples to fool machine learning models. Unlike a standard $\ell_p$-ball threat model, a functional adversarial threat model allows only a single function to be used to perturb input features to produce an adversarial example. For example, a functional adversarial attack applied on colors of an image can change all red pixels simultaneously to light red. Such global uniform changes in images can be less perceptible than perturbing pixels of the image individually. For simplicity, we refer to functional adversarial attacks on image colors as ReColorAdv, which is the main focus of our experiments. We show that functional threat models can be combined with existing additive ($\ell_p$) threat models to generate stronger threat models that allow both small, individual perturbations and large, uniform changes to an input. Moreover, we prove that such combinations encompass perturbations that would not be allowed in either constituent threat model. In practice, ReColorAdv can significantly reduce the accuracy of a ResNet-32 trained on CIFAR-10. Furthermore, to the best of our knowledge, combining ReColorAdv with other attacks leads to the strongest existing attack even after adversarial training. An implementation of ReColorAdv is available at https://github.com/cassidylaidlaw/ReColorAdv .


title: Generality of Lieb's Concavity Theorem

abstract: We show that Lieb's concavity theorem holds more generally for any unitarily invariant matrix function $\phi:\mathbf{H}^n_+\rightarrow \mathbb{R}$ that is monotone and concave. Concretely, we prove the joint concavity of the function $(A,B) \mapsto\phi\big[(B^\frac{qs}{2}K^*A^{ps}KB^\frac{qs}{2})^{\frac{1}{s}}\big] $ on $\mathbf{H}_+^m\times\mathbf{H}_+^n$, for any $K\in \mathbb{C}^{m\times n},s\in(0,1],p,q\in[0,1], p+q\leq 1$.


title: Counterfactual Analysis under Partial Identification Using Locally Robust Refinement

abstract: Structural models that admit multiple reduced forms, such as game-theoretic models with multiple equilibria, pose challenges in practice, especially when parameters are set-identified and the identified set is large. In such cases, researchers often choose to focus on a particular subset of equilibria for counterfactual analysis, but this choice can be hard to justify. This paper shows that some parameter values can be more "desirable" than others for counterfactual analysis, even if they are empirically equivalent given the data. In particular, within the identified set, some counterfactual predictions can exhibit more robustness than others, against local perturbations of the reduced forms (e.g. the equilibrium selection rule). We provide a representation of this subset which can be used to simplify the implementation. We illustrate our message using moment inequality models, and provide an empirical application based on a model with top-coded data.


title: Periodically Driven Many-Body Systems: A Floquet Density Matrix Renormalization Group Study

abstract: Driving a quantum system periodically in time can profoundly alter its long-time correlations and give rise to exotic quantum states of matter. The complexity of the combination of many-body correlations and dynamic manipulations has the potential to uncover a whole field of new phenomena, but the theoretical and numerical understanding becomes extremely difficult. We now propose a promising numerical method by generalizing the density matrix renormalization group to a superposition of Fourier components of periodically driven many-body systems using Floquet theory. With this method we can study the full time-dependent quantum solution in a large parameter range for all evolution times, beyond the commonly used high-frequency approximations. Numerical results are presented for the isotropic Heisenberg antiferromagnetic spin-1/2 chain under both local(edge) and global driving for spin-spin correlations and temporal fluctuations. As the frequency is lowered, we demonstrate that more and more Fourier components become relevant and determine strong length- and frequency-dependent changes of the quantum correlations that cannot be described by effective static models.


title: A Giant Intragroup Nebula Hosting a Damped Lya Absorber at z=0.313

abstract: This paper reports the discovery of spatially-extended line-emitting nebula, reaching to ~100 physical kpc (pkpc) from a damped Ly{\alpha} absorber (DLA) at z_DLA=0.313 along the sightline toward QSO PKS1127-145 (z_QSO=1.188). This DLA was known to be associated with a galaxy group of dynamical mass M_group ~3e12 M_sun, but its physical origin remained ambiguous. New wide-field integral field observations revealed a giant nebula detected in [OII], H{\beta}, [OIII], H{\alpha}, and [NII] emission, with the line-emitting gas following closely the motions of group galaxies. One of the denser streams passes directly in front of the QSO with kinematics consistent with the absorption profiles recorded in the QSO echelle spectra. The emission morphology, kinematics, and line ratios of the nebula suggest that shocks and turbulent mixing layers, produced as a result of stripped gaseous streams moving at supersonic speed across the ambient hot medium, contribute significantly to the ionization of the gas. While the DLA may not be associated with any specific detected member of the group, both the kinematic and dust properties are consistent with the DLA originating in streams of gas stripped from sub-L* group members at <~25 pkpc from the QSO sightline. This study demonstrates that gas stripping in low-mass galaxy groups is effective in releasing metal-enriched gas from star-forming regions, producing absorption systems in QSO spectra, and that combining absorption and emission-line observations provides an exciting new opportunity for studying gas and galaxy co-evolution.


title: A New Approach to Free-Form Cluster Lens Modeling Inspired by the JPEG Image Compression Method

abstract: I propose a new approach to free-form cluster lens modeling that is inspired by the JPEG image compression method. This approach is motivated specifically by the need for accurate modeling of high-magnification regions in galaxy clusters. Existing modeling methods may struggle in these regions due to their limited flexibility in the parametrization of the lens, even for a wide variety of free-form methods. This limitation especially hinders the characterization of faint galaxies at high redshifts, which have important implications for the formation of the first galaxies and even for the nature of dark matter. JPEG images are extremely accurate representations of their original, uncompressed counterparts but use only a fraction of number of parameters to represent that information. Its relevance is immediately obvious to cluster lens modeling. Using this technique, it is possible to construct flexible models that are capable of accurately reproducing the true mass distribution using only a small number of free parameters. Transferring this well-proven technology to cluster lens modeling, I demonstrate that this `JPEG parametrization' is indeed flexible enough to accurately approximate an N-body simulated cluster.


title: Star formation quenching imprinted on the internal structure of naked red nuggets

abstract: The formation and assembly process of massive galaxies is a combination of two phases: an initial in-situ-dominated one followed by an ex-situ-dominated evolution. Separating these two contributions is therefore crucial to understand the baryonic cycle within massive halos. A recently discovered population of so-called naked red nuggets, galaxies that shortcut the ex-situ stage preserving their pristine properties, presents a unique opportunity to study in detail star formation in massive galaxies without the confounding effect of later accretion. We investigate the spatially resolved star formation histories of a sample of 12 naked red nuggets. We measure how their radial light distributions, star formation rates and central densities evolved in time. We find that, while forming stars, red nuggets become gradually more concentrated, reaching a maximum concentration at quenching. After being quenched, they kept forming stars in a more disky-like configuration. Our measurements suggest that super-massive black holes and host galaxies grow their mass in a self-regulated way until a characteristic Mbh/Mhalo is reached. Once black holes are massive enough, red nuggets get quenched and depart from the star formation main sequence. While in the main sequence, red nuggets evolve at roughly constant star formation rate. This can explain up to ~0.3 dex of the scatter of the star formation main sequence, as well as its higher normalization observed in the early Universe. Hence, our results suggest that the main sequence is composed of populations of galaxies at different evolutionary stages, and that the scatter is therefore due to secular processes.


title: The Grism Lens-Amplified Survey from Space (GLASS). XIII. G800L optical spectra from the parallel fields

abstract: We present a catalogue of 22755 objects with slitless, optical, Hubble Space Telescope (HST) spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS). The data cover $\sim$220 sq. arcmin to 7-orbit ($\sim$10 ks) depth in 20 parallel pointings of the Advanced Camera for Survey's G800L grism. The fields are located 6' away from 10 massive galaxy clusters in the HFF and CLASH footprints. Thirteen of the fields have ancillary HST imaging from these or other programs to facilitate a large number of applications, from studying metal distributions at $z\sim0.5$, to quasars at $z\sim4$, to the star formation histories of hundreds of galaxies in between. The spectroscopic catalogue has a median redshift of $\langle z\rangle=0.6$ with a median uncertainty of $\Delta z / (1+z)\lesssim2\%$ at $\rm F814W\lesssim23$ AB. Robust continuum detections reach a magnitude fainter. The 5 $\sigma$ limiting line flux is $f_{\rm lim}\approx5\times10^{-17}\rm~erg~s^{-1}~cm^{-2}$ and half of all sources have 50% of pixels contaminated at $\lesssim$1%. All sources have 1- and 2-D spectra, line fluxes/uncertainties and identifications, redshift probability distributions, spectral models, and derived narrow-band emission line maps from the Grism Redshift and Line Analysis tool (GRIZLI). We provide other basic sample characterisations, show data examples, and describe sources and potential investigations of interest. All data and products will be available online along with software to facilitate their use.


title: Modelling the chemical enrichment of Population III supernovae: The origin of the metals in near-pristine gas clouds

abstract: The most metal-poor, high redshift damped Lyman-alpha systems (DLAs) provide a window to study some of the first few generations of stars. In this paper, we present a novel model to investigate the chemical enrichment of the near-pristine DLA population. This model accounts for the mass distribution of the enriching stellar population, the typical explosion energy of their supernovae, and the average number of stars that contribute to the enrichment of these DLAs. We conduct a maximum likelihood analysis of these model parameters using the observed relative element abundances ([C/O], [Si/O], and [Fe/O]) of the 11 most metal-poor DLAs currently known. We find that the mass distribution of the stars that have enriched this sample of metal-poor DLAs can be well-described by a Salpeter-like IMF slope at M > 10 M_sun and that a typical metal-poor DLA has been enriched by < 72 massive stars (95 per cent confidence), with masses < 40 M_sun. The inferred typical explosion energy (E_exp = 1.8^{+0.3}_{-0.2}x10^51 erg) is somewhat lower than that found by recent works that model the enrichment of metal-poor halo stars. These constraints suggest that some of the metal-poor DLAs in our sample may have been enriched by Population II stars. Using our enrichment model, we also infer some of the typical physical properties of the most metal-poor DLAs. We estimate that the total stellar mass content is log10(M_*/M_sun) = 3.5^{+0.3}_{-0.4} and the total gas mass is log10(M_gas/M_sun) = 7.0^{+0.3}_{-0.4} for systems with a relative oxygen abundance [O/H] ~ -3.0.


title: Annihilation Signatures of Hidden Sector Dark Matter Within Early-Forming Microhalos

abstract: If the dark matter is part of a hidden sector with only very feeble couplings to the Standard Model, the lightest particle in the hidden sector will generically be long-lived and could come to dominate the energy density of the universe prior to the onset of nucleosynthesis. During this early matter-dominated era, density perturbations will grow more quickly than otherwise predicted, leading to a large abundance of sub-earth-mass dark matter microhalos. Since the dark matter does not couple directly to the Standard Model, the minimum halo mass is much smaller than expected for weakly interacting dark matter, and the smallest halos could form during the radiation-dominated era. In this paper, we calculate the evolution of density perturbations within the context of such hidden sector models and use a series of $N$-body simulations to determine the outcome of nonlinear collapse during radiation domination. The resulting microhalos are extremely dense, which leads to very high rates of dark matter annihilation and to large indirect detection signals that resemble those ordinarily predicted for decaying dark matter. We find that the Fermi Collaboration's measurement of the high-latitude gamma-ray background rules out a wide range of parameter space within this class of models. The scenarios that are most difficult to constrain are those that feature a very long early matter-dominated era; if microhalos form prior to the decay of the unstable hidden sector matter, the destruction of these microhalos effectively heats the dark matter, suppressing the later formation of microhalos.


title: Effect of spin on the inspiral of binary neutron stars

abstract: We perform long-term simulations of spinning binary neutron stars, with our highest dimensionless spin being $\chi \sim 0.32$. To assess the importance of spin during the inspiral we vary the spin, and also use two equations of state, one that consists of plain nuclear matter and produces compact stars (SLy), and a hybrid one that contains both nuclear and quark matter and leads to larger stars (ALF2). Using high resolution that has grid spacing $\Delta x\sim 98$ m on the finest refinement level, we find that the effects of spin in the phase evolution of a binary system can be larger than the one that comes from tidal forces. Our calculations demonstrate explicitly that although tidal effects are dominant for small spins ($\lesssim 0.1$), this is no longer true when the spins are larger, but still much smaller than the Keplerian limit.


title: X-ray imaging polarimetry with a 2.5-$\mathrm{\mu}$m pixel CMOS sensor for visible light at room temperature

abstract: X-ray polarimetry in astronomy has not been exploited well, despite its importance. The recent innovation of instruments is changing this situation. We focus on a complementary MOS (CMOS) pixel detector with small pixel size and employ it as an x-ray photoelectron tracking polarimeter. The CMOS detector we employ is developed by GPixel Inc., and has a pixel size of 2.5 $\mathrm{\mu}$m $\times$ 2.5 $\mathrm{\mu}$m. Although it is designed for visible light, we succeed in detecting x-ray photons with an energy resolution of 176 eV (FWHM) at 5.9 keV at room temperature and the atmospheric condition. We measure the x-ray detection efficiency and polarimetry sensitivity by irradiating polarized monochromatic x-rays at BL20B2 in SPring-8, the synchrotron radiation facility in Japan. We obtain modulation factors of 7.63% $\pm$ 0.07% and 15.5% $\pm$ 0.4% at 12.4 keV and 24.8 keV, respectively. It demonstrates that this sensor can be used as an x-ray imaging spectrometer and polarimeter with the highest spatial resolution ever tested.


title: Parameterization of tensor network contraction

abstract: We present a conceptually clear and algorithmically useful framework for parameterizing the costs of tensor network contraction. Our framework is completely general, applying to tensor networks with arbitrary bond dimensions, open legs, and hyperedges. The fundamental objects of our framework are rooted and unrooted contraction trees, which represent classes of contraction orders. Properties of a contraction tree correspond directly and precisely to the time and space costs of tensor network contraction. The properties of rooted contraction trees give the costs of parallelized contraction algorithms. We show how contraction trees relate to existing tree-like objects in the graph theory literature, bringing to bear a wide range of graph algorithms and tools to tensor network contraction. Independent of tensor networks, we show that the edge congestion of a graph is almost equal to the branchwidth of its line graph.


title: Detecting Dark Matter with Aharonov-Bohm

abstract: While the evidence for dark matter continues to grow, the nature of the dark matter remains a mystery. A dark $U(1)_D$ gauge theory can have a small kinetic mixing with the visible photon which provides a portal to the dark sector. Magnetic monopoles of the dark $U(1)_D$ can obtain small magnetic couplings to our photon through this kinetic mixing. This coupling is only manifest below the mass of the dark photon; at these scales the monopoles are bound together by tubes of dark magnetic flux. These flux tubes can produce phase shifts in Aharonov-Bohm type experiments. We outline how this scenario might be realized, examine the existing constraints, and quantify the experimental sensitivity required to detect magnetic dipole dark matter in this novel way.


title: How Pulses in Short Gamma-Ray Bursts Constrain HMXRB Evolution

abstract: We demonstrate how pulse structures in Short gamma-ray bursts (SGRBs), coupled with observations of GRB/GW 170817A, constrain the geometries of dying HMXRB systems composed of merging neutron stars.


title: Many-body Chaos in Thermalised Fluids

abstract: Linking thermodynamic variables like temperature $T$ and the measure of chaos, the Lyapunov exponents $\lambda$, is a question of fundamental importance in many-body systems. By using nonlinear fluid equations in one and three dimensions, we prove that in thermalised flows $\lambda \propto \sqrt{T}$, in agreement with results from frustrated spin systems. This reveals an underlying universality and provides evidence for recent conjectures on the thermal scaling of $\lambda$. We also reconcile seemingly disparate effects -- equilibration on one hand and pushing systems out-of-equilibrium on the other -- of many-body chaos by relating $\lambda$ to $T$ through the dynamical structures of the flow.


title: ITO-based Electro-absorption Modulator for Photonic Neural Activation Function

abstract: Recently integrated optics has become an intriguing platform for implementing machine learning algorithms and inparticular neural networks. Integrated photonic circuits can straightforwardly perform vector-matrix multiplicationswith high efficiency and low power consumption by using weighting mechanism through linear optics. Although,this can not be said for the activation function which requires either nonlinear optics or an electro-optic module withan appropriate dynamic range. Even though all-optical nonlinear optics is potentially faster, its current integrationis challenging and is rather inefficient. Here we demonstrate an electro-absorption modulator based on an IndiumTin Oxide layer, whose dynamic range is used as nonlinear activation function of a photonic neuron. The nonlinearactivation mechanism is based on a photodiode, which integrates the weighed products, and whose photovoltage drivesthe elecro-absorption modulator. The synapse and neuron circuit is then constructed to execute a 200-node MNISTclassification neural network used for benchmarking the nonlinear activation function and compared with an equivalentelectronic module.


title: Graphene epsilon-near-zero plasmonic crystals

abstract: Plasmonic crystals are a class of optical metamaterials that consist of engineered structures at the sub-wavelength scale. They exhibit optical properties that are not found under normal circumstances in nature, such as negative-refractive-index and epsilon-near-zero (ENZ) behavior. Graphene-based plasmonic crystals present linear, elliptical, or hyperbolic dispersion relations that exhibit ENZ behavior, normal or negative-index diffraction. The optical properties can be dynamically tuned by controlling the operating frequency and the doping level of graphene. We propose a construction approach to expand the frequency range of the ENZ behavior. We demonstrate how the combination of a host material with an optical Lorentzian response in combination with a graphene conductivity that follows a Drude model leads to an ENZ condition spanning a large frequency range.


title: Magnetic inhibition of convection in O star envelopes

abstract: It has been suggested that the absence of macroturbulence in the atmosphere of NGC 1624 - 2 is due its strong magnetic field (the strongest known for a massive O star) suppressing convection in its outer layers, removing the mechanism thought responsible for the observed macroturbulence in stars with lower field strengths. Here, we develop and apply a criterion for a uniform magnetic field to suppress convection in stellar envelopes in which radiation pressure is a significant contributor to hydrostatic balance. We find upper mass limits of ~55 Msun and ~30 Msun for magnetic suppression to be possible in Zero Age Main Sequence and Terminal Age Main Sequence stars, respectively. For evolved stars, magnetic suppression of convection can significantly alter the stars' evolution. For NGC 1624 - 2, we find a polar dipole strength of 16.5 +/- 5.9 kG is required to suppress convection, in good agreement with the value ~20 kG measured by spectropolarimetry.


title: Predicatively unprovable termination of the Ackermannian Goodstein process

abstract: The classical Goodstein process gives rise to long but finite sequences of natural numbers whose termination is not provable in Peano arithmetic. In this manuscript we consider a variant based on the Ackermann function. We show that Ackermannian Goodstein sequences eventually terminate, but this fact is not provable using predicative means.


title: Exact Recovery in Block Spin Ising Models at the Critical Line

abstract: We show how to exactly reconstruct the block structure at the critical line in the so-called Ising block model. This model was re-introduced by Berthet, Rigollet and Srivastava in a recent paper. There the authors show how to exactly reconstruct blocks away from the critical line and they give an upper and a lower bound on the number of observations one needs; thereby they establish a minimax optimal rate (up to constants). Our technique relies on a combination of their methods with fluctuation results for block spin Ising models. The latter are extended to the full critical regime. We find that the number of necessary observations depends on whether the interaction parameter between two blocks is positive or negative: In the first case, there are about $N \log N$ observations required to exactly recover the block structure, while in the latter $\sqrt N \log N$ observations suffice.


title: Autonomous multipartite entanglement engines

abstract: The generation of genuine multipartite entangled states is challenging in practice. Here we explore a new route to this task, via autonomous entanglement engines which use only incoherent coupling to thermal baths and time-independent interactions. We present a general machine architecture, which allows for the generation of a broad range of multipartite entangled states in a heralded manner. Specifically, given a target multiple-qubit state, we give a sufficient condition ensuring that it can be generated by our machine. We discuss the cases of Greenberger-Horne-Zeilinger, Dicke and cluster states in detail. These results demonstrate the potential of purely thermal resources for creating multipartite entangled states useful for quantum information processing.


title: Not a simple relationship between Neptune's migration speed and Kuiper belt inclination excitation

abstract: We present numerical simulations of giant planet migration in our solar system and examine how the speed of planetary migration affects inclinations in the resulting population of small bodies (test particles) scattered outward and subsequently captured into Neptune's 3:2 mean motion resonance (the Plutinos) as well as the hot classical Kuiper belt population. We do not find a consistent relationship between the degree of test particle inclination excitation and e-folding planet migration timescales in the range 5-50 Myr. Our results present a counter-example to Nesvorny 2015's finding that the Plutino and hot classical inclinations showed a marked increase with increasing e-folding timescales for Neptune's migration. We argue that these differing results are likely due to differing secular architectures of the giant planets during and after migration. Small changes in the planets' initial conditions and differences in the numerical implementation of planet migration can result in different amplitudes of the planets' inclination secular modes, and this can lead to different final inclination distributions for test particles in the simulations. We conclude that the observed large inclination dispersion of Kuiper belt objects does not require Neptune's migration to be slow; planetary migration with e-folding timescales of 5, 10, 30, and 50 Myr can all yield inclination dispersions similar to the observed Plutino and hot classical populations, with no correlation between the degree of inclination excitation and migration speed.


title: Highly efficient optical add-drop filter with an angle-polished fiber coupler

abstract: Microbubble whispering-gallery resonators have shown great promise in fiber-optic communications because of their low confinement loss and hollow cores, which allow for facile stress-based tunability. Usually, the transmission spectrum of taper-coupled microbubbles contains closely spaced modes due to the relatively large radii and oblate geometry of microbubbles. In this letter, we develop an optical add-drop filter using a microbubble coupled to fiber taper and angle-polished fiber waveguides. Because of the extra degree of freedom in the angle of its polish surface, the angle-polished fiber can be used for the discriminatory excitation of certain radial-order modes in the optical microcavity, reducing the high modal density and enhancing add-drop selectivity. Our robust and tunable add-drop filter demonstrated a drop efficiency of 85.9% and quality factor of 2 x 10^7, corresponding to a linewidth of 9.68 MHz. As a proof of concept, the drop frequency was tuned using internal aerostatic pressure at a rate of 7.3 +/- 0.2 GHz/bar with no diminishing effects on the add-drop filter performance.


title: Minimum-Margin Active Learning

abstract: We present a new active sampling method we call min-margin which trains multiple learners on bootstrap samples and then chooses the examples to label based on the candidates' minimum margin amongst the bootstrapped models. This extends standard margin sampling in a way that increases its diversity in a supervised manner as it arises from the model uncertainty. We focus on the one-shot batch active learning setting, and show theoretically and through extensive experiments on a broad set of problems that min-margin outperforms other methods, particularly as batch size grows.


title: IoT Skullfort: Exploring the Impact of Internet Connected Cosplay

abstract: In this paper, we explore the potential impact of Internet of Things (IoT) technology may have on the cosplay community. We developed a costume (an IoT Skullfort) and embedded IoT technology to enhance its capabilities and user interactions. Sensing technologies are widely used in many different wearable domains including cosplay scenarios. However, in most of these scenarios, typical interaction pattern is that the costume responds to its environment or the player's behaviour (e.g., colour of lights may get changed when player moves hands). In contrast, our research focuses on exploring scenarios where the audience (third party) get to manipulate the costume behaviour (e.g., the audience get to change the colour of the Skullfort using a mobile application). We believe such an audience (third party) influenced cosplay brings new opportunities for enhanced entertainment. However, it also creates significant challenges. We report the results gathered through a focus group conducted in collaboration with cosplay community experts.


title: Dynamical system analysis of generalized energy-momentum-squared gravity

abstract: In this work we have investigated the dynamics of a recent modification to the general theory of relativity, the energy-momentum squared gravity model $f(R,\mathbf{T^2})$, where $R$ represents the scalar curvature and $\mathbf{T^2}$ the square of the energy-momentum tensor. By using dynamical system analysis for various types of gravity functions $f(R,\mathbf{T^2})$, we have studied the structure of the phase space and the physical implications of the energy-momentum squared coupling. In the first case of functional where $f(R,\mathbf{T^2})=f_0 R^n(\mathbf{T^2})^m$, with $f_0$ constant, we have shown that the phase space structure has a reduced complexity, with a high sensitivity to the values of the $m$ and $n$ parameters. Depending on the values of the $m$ and $n$ parameters, the model exhibits various cosmological epochs, corresponding to matter eras, solutions associated with an accelerated expansion, or decelerated periods. The second model studied corresponds to the $f(R,\mathbf{T^2})=\alpha R^n+\beta\mathbf{(T^2)}^m$ form with $\alpha, \beta$ constant parameters. In this case, a richer phase space structure is obtained which can recover different cosmological scenarios, associated to matter eras, de--Sitter solutions, and dark energy epochs. Hence, this model represents an interesting cosmological model which can explain the current evolution of the Universe and the emergence of the accelerated expansion as a geometrical consequence.


title: Independent Component Analysis based on multiple data-weighting

abstract: Independent Component Analysis (ICA) - one of the basic tools in data analysis - aims to find a coordinate system in which the components of the data are independent. In this paper we present Multiple-weighted Independent Component Analysis (MWeICA) algorithm, a new ICA method which is based on approximate diagonalization of weighted covariance matrices. Our idea is based on theoretical result, which says that linear independence of weighted data (for gaussian weights) guarantees independence. Experiments show that MWeICA achieves better results to most state-of-the-art ICA methods, with similar computational time.


title: Human-Usable Password Schemas: Beyond Information-Theoretic Security

abstract: Password users frequently employ passwords that are too simple, or they just reuse passwords for multiple websites. A common complaint is that utilizing secure passwords is too difficult. One possible solution to this problem is to use a password schema. Password schemas are deterministic functions which map challenges (typically the website name) to responses (passwords). Previous work has been done on developing and analyzing publishable schemas, but these analyses have been information-theoretic, not complexity-theoretic; they consider an adversary with infinite computing power. We perform an analysis with respect to adversaries having currently achievable computing capabilities, assessing the realistic practical security of such schemas. We prove for several specific schemas that a computer is no worse off than an infinite adversary and that it can successfully extract all information from leaked challenges and their respective responses, known as challenge-response pairs. We also show that any schema that hopes to be secure against adversaries with bounded computation should obscure information in a very specific way, by introducing many possible constraints with each challenge-response pair. These surprising results put the analyses of password schemas on a more solid and practical footing.


title: Pseudo-Riemannian geometry embeds information geometry in optimal transport

abstract: Optimal transport and information geometry both study geometric structures on spaces of probability distributions. Optimal transport characterizes the cost-minimizing movement from one distribution to another, while information geometry originates from coordinate-invariant properties of statistical inference. Their connections and applications in statistics and machine learning have started to gain more attention. In this paper we give a new differential geometric connection between the two fields. Namely, the pseudo-Riemannian framework of Kim and McCann, a geometric perspective on the fundamental Ma-Trudinger-Wang (MTW) condition in the regularity theory of optimal transport maps, encodes the dualistic structure of statistical manifold. This general relation is described using the natural framework of $c$-divergence, a divergence defined by an optimal transport map. As a by-product, we obtain a new information-geometric interpretation of the MTW tensor. This connection sheds light on old and new aspects of information geometry. The dually flat geometry of Bregman divergence corresponds to the quadratic cost and the pseudo-Euclidean space, and the $L^{(\alpha)}$-divergence introduced by Pal and the first author has constant sectional curvature in a sense to be made precise. In these cases we give a geometric interpretation of the information-geometric curvature in terms of the divergence between a primal-dual pair of geodesics.


title: Greedy inference with structure-exploiting lazy maps

abstract: We propose a framework for solving high-dimensional Bayesian inference problems using \emph{structure-exploiting} low-dimensional transport maps or flows. These maps are confined to a low-dimensional subspace (hence, lazy), and the subspace is identified by minimizing an upper bound on the Kullback--Leibler divergence (hence, structured). Our framework provides a principled way of identifying and exploiting low-dimensional structure in an inference problem. It focuses the expressiveness of a transport map along the directions of most significant discrepancy from the posterior, and can be used to build deep compositions of lazy maps, where low-dimensional projections of the parameters are iteratively transformed to match the posterior. We prove weak convergence of the generated sequence of distributions to the posterior, and we demonstrate the benefits of the framework on challenging inference problems in machine learning and differential equations, using inverse autoregressive flows and polynomial maps as examples of the underlying density estimators.


title: Optimal Control of Fractional Elliptic PDEs with State Constraints and Characterization of the dual of Fractional Order Sobolev Spaces

abstract: This paper introduces the notion of state constraints for optimal control problems governed by fractional elliptic PDEs of order $s \in (0,1)$. There are several mathematical tools that are developed during the process to study this problem, for instance, the characterization of the dual of the fractional order Sobolev spaces and well-posedness of fractional PDEs with measure-valued datum. These tools are widely applicable. We show well-posedness of the optimal control problem and derive the first order optimality conditions. Notice that the adjoint equation is a fractional PDE with measure as the right-hand-side datum. We use the characterization of the fractional order dual spaces to study the regularity of the state and adjoint equations. We emphasize that the classical case ($s=1$) was considered by E. Casas in \cite{ECasas_1986a} but almost none of the existing results are applicable to our fractional case.


title: K-stability of continuous C(X)-algebras

abstract: A C*-algebra is said to be K-stable if its nonstable K-groups are naturally isomorphic to the usual K-theory groups. We study continuous $C(X)$-algebras, each of whose fibers are K-stable. We show that such an algebra is itself K-stable under the assumption that the underlying space $X$ is compact, metrizable, and of finite covering dimension.


title: A variation of Broyden Class methods using Householder adaptive transforms

abstract: In this work we introduce and study novel Quasi Newton minimization methods based on a Hessian approximation Broyden Class-\textit{type} updating scheme, where a suitable matrix $\tilde{B}_k$ is updated instead of the current Hessian approximation $B_k$. We identify conditions which imply the convergence of the algorithm and, if exact line search is chosen, its quadratic termination. By a remarkable connection between the projection operation and Krylov spaces, such conditions can be ensured using low complexity matrices $\tilde{B}_k$ obtained projecting $B_k$ onto algebras of matrices diagonalized by products of two or three Householder matrices adaptively chosen step by step. Extended experimental tests show that the introduction of the adaptive criterion, which theoretically guarantees the convergence, considerably improves the robustness of the minimization schemes when compared with a non-adaptive choice; moreover, they show that the proposed methods could be particularly suitable to solve large scale problems where $L$-$BFGS$ performs poorly.


title: Improved Gauss law model and in-medium heavy quarkonium at finite density and velocity

abstract: We explore the in-medium properties of heavy-quarkonium states at finite baryo-chemical potential and finite transverse momentum based on a modern complex-valued potential model. Our starting point is a novel, rigorous derivation of the generalized Gauss law for in-medium quarkonium, combining the non-perturbative physics of the vacuum bound state with a weak coupling description of the medium degrees of freedom. Its relation to previous models in the literature is discussed. We show that our approach is able to reproduce the complex lattice QCD heavy quark potential even in the non-perturbative regime, using a single temperature dependent parameter, the Debye mass $m_D$. After vetting the Gauss-law potential with state-of-the-art lattice QCD data, we extend it to the regime of finite baryon density and finite velocity, currently inaccessible to first principles simulations. In-medium spectral functions computed from the Gauss-law potential are subsequently used to estimate the $\psi'/J/\psi$ ratio in heavy-ion collisions at different beam energies and transverse momenta. We find qualitative agreement with the predictions from the statistical model of hadronization for the $\sqrt{s_{NN}}$ dependence and a mild dependence on the transverse momentum.


title: Whitney Numbers for Poset Cones

abstract: Hyperplane arrangements dissect $\mathbb{R}^n$ into connected components called chambers, and a well-known theorem of Zaslavsky counts chambers as a sum of nonnegative integers called Whitney numbers of the first kind. His theorem generalizes to count chambers within any cone defined as the intersection of a collection of halfspaces from the arrangement, leading to a notion of Whitney numbers for each cone. This paper focuses on cones within the braid arrangement, consisting of the reflecting hyperplanes $x_i=x_j$ inside $\mathbb{R}^n$ for the symmetric group, thought of as the type $A_{n-1}$ reflection group. Here cones correspond to posets, chambers within the cone correspond to linear extensions of the poset, and the Whitney numbers of the cone interestingly refine the number of linear extensions of the poset. We interpret this refinement for all posets as counting linear extensions according to a statistic that generalizes the number of left-to-right maxima of a permutation. When the poset is a disjoint union of chains, we interpret this refinement differently, using Foata's theory of cycle decomposition for multiset permutations, leading to a simple generating function compiling these Whitney numbers.


title: Long-Step Path-Following Algorithm for Quantum Information Theory: Some Numerical Aspects and Applications

abstract: We consider some important computational aspects of the long-step path-following algorithm developed in our previous work and show that a broad class of complicated optimization problems arising in quantum information theory can be solved using this approach. In particular, we consider one difficult and important optimization problem in quantum key distribution and show that our method can solve problems of this type much faster in comparison with (very few) available options.


title: Visual Understanding and Narration: A Deeper Understanding and Explanation of Visual Scenes

abstract: We describe the task of Visual Understanding and Narration, in which a robot (or agent) generates text for the images that it collects when navigating its environment, by answering open-ended questions, such as 'what happens, or might have happened, here?'


title: On the Parallelization of Triangular Decomposition of Polynomial Systems

abstract: We discuss the parallelization of algorithms for solving polynomial systems symbolically by way of triangular decomposition. Algorithms for solving polynomial systems combine low-level routines for performing arithmetic operations on polynomials and high-level procedures which produce the different components (points, curves, surfaces) of the solution set. The latter "component-level" parallelization of triangular decompositions, our focus here, belongs to the class of dynamic irregular parallel applications. Possible speedup factors depend on geometrical properties of the solution set (number of components, their dimensions and degrees); these algorithms do not scale with the number of processors. In this paper we combine two different concurrency schemes, the fork-join model and producer-consumer patterns, to better capture opportunities for component-level parallelization. We report on our implementation with the publicly available BPAS library. Our experimentation with 340 systems yields promising results.


title: Discovery of four apparently cold dusty galaxies at z=3.62-5.85 in the COSMOS field: direct evidence of CMB impact on high-redshift galaxy observables

abstract: We report Atacama Large Millimetre Array (ALMA) observations of four high-redshift dusty star-forming galaxy candidates selected from far-Infrared (FIR)/submm observations in the COSMOS field. We securely detect all galaxies in the continuum and spectroscopically confirm them at z=3.62--5.85 using ALMA 3mm line scans, detecting multiple CO and/or [CI] transitions. This includes the most distant dusty galaxy currently known in the COSMOS field, ID85001929 at z=5.847. These redshifts are lower than we had expected as these galaxies have substantially colder dust temperatures (i.e., their SEDs peak at longer rest frame wavelengths) than most literature sources at z>4. The observed cold dust temperatures are best understood as evidence for optically thick dust continuum in the FIR, rather than the result of low star formation efficiency with rapid metal enrichment. We provide direct evidence that, given their cold spectral energy distributions, CMB plays a significant role biasing their observed Rayleigh-Jeans (RJ) slopes to unlikely steep values and, possibly, reducing their CO fluxes by a factor of two. We recover standard RJ slopes when the CMB contribution is taken into account. High resolution ALMA imaging shows compact morphology and evidence for mergers. This work reveals a population of cold dusty star-forming galaxies that were under-represented in current surveys, and are even colder than typical Main Sequence galaxies at the same redshift. High FIR dust optical depth might be a widespread feature of compact starbursts at any redshift.


title: Table2Vec: Neural Word and Entity Embeddings for Table Population and Retrieval

abstract: Tables contain valuable knowledge in a structured form. We employ neural language modeling approaches to embed tabular data into vector spaces. Specifically, we consider different table elements, such caption, column headings, and cells, for training word and entity embeddings. These embeddings are then utilized in three particular table-related tasks, row population, column population, and table retrieval, by incorporating them into existing retrieval models as additional semantic similarity signals. Evaluation results show that table embeddings can significantly improve upon the performance of state-of-the-art baselines.


title: Bayesian Profiling Multiple Imputation for Missing Electronic Health Records

abstract: Electronic health records (EHRs) are increasingly used for clinical and comparative effectiveness research, but suffer from missing data. Motivated by health services research on diabetes care, we seek to increase the quality of EHRs by focusing on missing values of longitudinal glycosylated hemoglobin (A1c), a key risk factor for diabetes complications and adverse events. Under the framework of multiple imputation (MI), we propose an individualized Bayesian latent profiling approach to capture A1c measurement trajectories subject to missingness. The proposed method is applied to EHRs of adult patients with diabetes in a large academic Midwestern health system between 2003 and 2013 and had Medicare A and B coverage. We combine MI inferences to evaluate the association of A1c levels with the incidence of acute adverse health events and examine patient heterogeneity across identified patient profiles. We investigate different missingness mechanisms and perform imputation diagnostics. Our approach is computationally efficient and fits flexible models that provide useful clinical insights.


title: Analyticity of the free energy for quantum Airy structures

abstract: It is shown that the free energy associated to a finite dimensional Airy structure is an analytic function at each finite order of the $\hbar$ expansion. Semiclassical series itself is in general divergent. Calculations are facilitated by putting the topological recursion equations into a form exhibiting more explicitly the semiclassical geometry. This formulation involves certain differential operators on the characteristic variety, which are found to satisfy a Lie algebra cocycle condition. It is proven that this cocycle is a coboundary. Developed formalism is applied in specific examples. In the case of a divergent $\hbar$ series, a simple resummation is performed. Analytic properties of the obtained partition functions are investigated.


title: Experimental observation of exceptional surface in synthetic dimensions with magnon polaritons

abstract: Exceptional points (EPs) are singularities of energy levels in non-Hermitian systems. In this Letter, we demonstrate the surface of EPs on a magnon polariton platform composed of coupled magnons and microwave photons. Our experiments show that EPs form a three-dimensional exceptional surface (ES) when the system is tuned in a four-dimensional synthetic space. We demonstrated that there exists an exceptional saddle point (ESP) in the ES which originates from the unique couplings between magnons and microwave photons. Such an ESP exhibits unique anisotropic behaviors in both the real and imaginary part of the eigenfrequencies. To the best of our knowledge, this is the first experimental observation of ES, opening up new opportunities for high-dimensional control of non-Hermitian systems.


title: Where Are We With Light Sterile Neutrinos?

abstract: We review the status of searches for sterile neutrinos in the $\sim 1$ eV range, with an emphasis on the latest results from short baseline oscillation experiments and how they fit within sterile neutrino oscillation models. We present global fit results to a three-active-flavor plus one-sterile-flavor model (3+1), where we find an improvement of $\Delta \chi^2=35$ for 3 additional parameters compared to a model with no sterile neutrino. This is a 5$\sigma$ improvement, indicating that an effect that is like that of a sterile neutrino is highly preferred by the data. However we note that separate fits to the appearance and disappearance oscillation data sets within a 3+1 model do not show the expected overlapping allowed regions in parameter space. This "tension" leads us to explore two options: 3+2, where a second additional mass state is introduced, and a 3+1+decay model, where the $\nu_4$ state can decay to invisible particles. The 3+1+decay model, which is also motivated by improving compatibility with cosmological observations, yields the larger improvement, with a $\Delta \chi^2=8$ for 1 additional parameter beyond the 3+1 model, which is a $2.6\sigma$ improvement. Moreover the tension between appearance and disappearance experiments is reduced compared to 3+1, although disagreement remains. In these studies, we use a frequentist approach and also a Bayesean method of finding credible regions. With respect to this tension, we review possible problems with the global fitting method. We note multiple issues, including problems with reproducing the experimental results, especially in the case of experiments that do not provide adequate data releases. We discuss an unexpected 5 MeV excess, observed in the reactor flux energy spectrum, that may be affecting the oscillation interpretation of the short baseline reactor data. We emphasize the care that must be taken in mapping to the true neutrino energy in the case of oscillation experiments that are subject to multiple interaction modes and nuclear effects. We point to problems with the "Parameter-Goodness-of-Fit test" that is used to quantify the tension. Lastly, we point out that analyses presenting limits often receive less scrutiny that signals. While we provide a snapshot of the status of sterile neutrino searches today and global fits to their interpretation, we emphasize that this is a fast-moving field. We briefly review experiments that are expected to report new data in the immediate future. Lastly, we consider the 5-year horizon, where we propose that decay-at-rest neutrino sources are the best method of finally resolving the confusing situation.


title: Interaction Trees: Representing Recursive and Impure Programs in Coq

abstract: "Interaction trees" (ITrees) are a general-purpose data structure for representing the behaviors of recursive programs that interact with their environments. A coinductive variant of "free monads," ITrees are built out of uninterpreted events and their continuations. They support compositional construction of interpreters from "event handlers", which give meaning to events by defining their semantics as monadic actions. ITrees are expressive enough to represent impure and potentially nonterminating, mutually recursive computations, while admitting a rich equational theory of equivalence up to weak bisimulation. In contrast to other approaches such as relationally specified operational semantics, ITrees are executable via code extraction, making them suitable for debugging, testing, and implementing software artifacts that are amenable to formal verification. We have implemented ITrees and their associated theory as a Coq library, mechanizing classic domain- and category-theoretic results about program semantics, iteration, monadic structures, and equational reasoning. Although the internals of the library rely heavily on coinductive proofs, the interface hides these details so that clients can use and reason about ITrees without explicit use of Coq's coinduction tactics. To showcase the utility of our theory, we prove the termination-sensitive correctness of a compiler from a simple imperative source language to an assembly-like target whose meanings are given in an ITree-based denotational semantics. Unlike previous results using operational techniques, our bisimulation proof follows straightforwardly by structural induction and elementary rewriting via an equational theory of combinators for control-flow graphs.


title: \emph{Ab initio} theory of nitrogen-vacancy center in diamond

abstract: Nitrogen-vacancy center in diamond is a solid state defect qubit with favorable coherence time up to room temperature which could be harnessed in several quantum enhanced sensor and quantum communication applications, and has a potential in quantum simulation and computing. The quantum control largely depends on the intricate details about the electronic structure and states of the nitrogen-vacancy center, radiative and non-radiative rates between these states and the coupling of these states to external spins, electrical, magnetic and strain fields and temperature. In this review paper it is shown how first principles calculations contributed to understanding the properties of nitrogen-vacancy center, and will be briefly discussed the issues to be solved towards full \emph{ab initio} description of solid state defect qubits.

categories: quant-ph cond-mat.str-el physics.atom-ph physics.comp-ph


title: Thinking Slow about Latency Evaluation for Simultaneous Machine Translation

abstract: Simultaneous machine translation attempts to translate a source sentence before it is finished being spoken, with applications to translation of spoken language for live streaming and conversation. Since simultaneous systems trade quality to reduce latency, having an effective and interpretable latency metric is crucial. We introduce a variant of the recently proposed Average Lagging (AL) metric, which we call Differentiable Average Lagging (DAL). It distinguishes itself by being differentiable and internally consistent to its underlying mathematical model.


title: Online Convex Optimization with Perturbed Constraints

abstract: This paper addresses Online Convex Optimization (OCO) problems where the constraints have additive perturbations that (i) vary over time and (ii) are not known at the time to make a decision. Perturbations may not be i.i.d. generated and can be used to model a time-varying budget or commodity in resource allocation problems. The problem is to design a policy that obtains sublinear regret while ensuring that the constraints are satisfied on average. To solve this problem, we present a primal-dual proximal gradient algorithm that has $O(T^\epsilon \vee T^{1-\epsilon})$ regret and $O(T^\epsilon)$ constraint violation, where $\epsilon \in [0,1)$ is a parameter in the learning rate. Our results match the bounds of previous work on OCO with time-varying constraints when $\epsilon = 1/2$; however, we (i) define the regret using a time-varying set of best fixed decisions; (ii) can balance between regret and constraint violation; and (iii) use an adaptive learning rate that allows us to run the algorithm for any time horizon.


title: DISCO: Depth Inference from Stereo using Context

abstract: Recent deep learning based approaches have outperformed classical stereo matching methods. However, current deep learning based end-to-end stereo matching methods adopt a generic encoder-decoder style network with skip connections. To limit computational requirement, many networks perform excessive down sampling, which results in significant loss of useful low-level information. Additionally, many network designs do not exploit the rich multi-scale contextual information. In this work, we address these aforementioned problems by carefully designing the network architecture to preserve required spatial information throughout the network, while at the same time achieve large effective receptive field to extract multiscale contextual information. For the first time, we create a synthetic disparity dataset reflecting real life images captured using a smartphone; this enables us to obtain state-of-the-art results on common real life images. The proposed model DISCO is pre-trained on the synthetic Scene Flow dataset and evaluated on popular benchmarks and our in-house dataset of challenging real life images. The proposed model outperforms existing state-of-the-art methods in terms of quality as well as quantitative metrics.


abstract: We consider the problem of decomposing a large covariance matrix into the sum of a low-rank matrix and a diagonally dominant matrix, and we call this problem the "Diagonally-Dominant Principal Component Analysis (DD-PCA)". DD-PCA is an effective tool for designing statistical methods for strongly correlated data. We showcase the use of DD-PCA in two statistical problems: covariance matrix estimation, and global detection in multiple testing. Using the output of DD-PCA, we propose a new estimator for estimating a large covariance matrix with factor structure. Thanks to a nice property of diagonally dominant matrices, this estimator enjoys the advantage of simultaneous good estimation of the covariance matrix and the precision matrix (by a plain inversion). A plug-in of this estimator to linear discriminant analysis and portfolio optimization yields appealing performance in real data. We also propose two new tests for testing the global null hypothesis in multiple testing when the $z$-scores have a factor covariance structure. Both tests first use DD-PCA to adjust the individual $p$-values and then plug in the adjusted $p$-values to the Higher Criticism (HC) test. These new tests significantly improve over the HC test and compare favorably with other existing tests. For computation of DD-PCA, we propose an iterative projection algorithm and an ADMM algorithm.