title: Simultaneous Preference and Metric Learning from Paired Comparisons

abstract: A popular model of preference in the context of recommendation systems is the so-called \emph{ideal point} model. In this model, a user is represented as a vector $\mathbf{u}$ together with a collection of items $\mathbf{x_1}, \ldots, \mathbf{x_N}$ in a common low-dimensional space. The vector $\mathbf{u}$ represents the user's "ideal point," or the ideal combination of features that represents a hypothesized most preferred item. The underlying assumption in this model is that a smaller distance between $\mathbf{u}$ and an item $\mathbf{x_j}$ indicates a stronger preference for $\mathbf{x_j}$. In the vast majority of the existing work on learning ideal point models, the underlying distance has been assumed to be Euclidean. However, this eliminates any possibility of interactions between features and a user's underlying preferences. In this paper, we consider the problem of learning an ideal point representation of a user's preferences when the distance metric is an unknown Mahalanobis metric. Specifically, we present a novel approach to estimate the user's ideal point $\mathbf{u}$ and the Mahalanobis metric from paired comparisons of the form "item $\mathbf{x_i}$ is preferred to item $\mathbf{x_j}$." This can be viewed as a special case of a more general metric learning problem where the location of some points are unknown a priori. We conduct extensive experiments on synthetic and real-world datasets to exhibit the effectiveness of our algorithm.


title: Duality theory for enriched Priestley spaces

abstract: The term Stone-type duality often refers to a dual equivalence between a category of lattices or other partially ordered structures on one side and a category of topological structures on the other. This paper is part of a larger endeavour that aims to extend a web of Stone-type dualities from ordered to metric structures and, more generally, to quantale-enriched categories. In particular, we improve our previous work and show how certain duality results for categories of [0,1]-enriched Priestley spaces and [0,1]-enriched relations can be restricted to functions. In a broader context, we investigate the category of quantale-enriched Priestley spaces and continuous functors, with emphasis on those properties which identify the algebraic nature of the dual of this category.


title: Singular BPS boundary conditions in $\mathcal{N} = (2,2)$ supersymmetric gauge theories

abstract: We derive general BPS boundary conditions in two-dimensional $\mathcal{N}=(2,2)$ supersymmetric gauge theories. We analyze the solutions of these boundary conditions, and in particular those that allow the bulk fields to have poles at the boundary. We also present the brane configurations for the half- and quarter-BPS boundary conditions of the $\mathcal{N}=(2,2)$ supersymmetric gauge theories in terms of branes in Type IIA string theory. We find that both A-type and B-type brane configurations are lifted to M-theory as a system of M2-branes ending on an M5-brane wrapped on a product of a holomorphic curve in $\mathbb{C}^2$ with a special Lagrangian 3-cycle in $\mathbb{C}^3$.


title: Composite Estimation for Quantile Regression Kink Models with Longitudinal Data

abstract: Kink model is developed to analyze the data where the regression function is twostage linear but intersects at an unknown threshold. In quantile regression with longitudinal data, previous work assumed that the unknown threshold parameters or kink points are heterogeneous across different quantiles. However, the location where kink effect happens tend to be the same across different quantiles, especially in a region of neighboring quantile levels. Ignoring such homogeneity information may lead to efficiency loss for estimation. In view of this, we propose a composite estimator for the common kink point by absorbing information from multiple quantiles. In addition, we also develop a sup-likelihood-ratio test to check the kink effect at a given quantile level. A test-inversion confidence interval for the common kink point is also developed based on the quantile rank score test. The simulation study shows that the proposed composite kink estimator is more competitive with the least square estimator and the single quantile estimator. We illustrate the practical value of this work through the analysis of a body mass index and blood pressure data set.


title: Distributed Synchronous Visualization Design: Challenges and Strategies

abstract: We reflect on our experiences as designers of COVID-19 data visualizations working in a distributed synchronous design space during the pandemic. This is especially relevant as the pandemic posed new challenges to distributed collaboration amidst civic lockdown measures and an increased dependency on spatially distributed teamwork across almost all sectors. Working from home being 'the new normal', we explored potential solutions for collaborating and prototyping remotely from our own homes using the existing tools at our disposal. Since members of our cross-disciplinary team had different technical skills, we used a range of synchronous remote design tools and methods. We aimed to preserve the richness of co-located collaboration such as face-to-face physical presence, body gestures, facial expressions, and the making and sharing of physical artifacts. While meeting over Zoom, we sketched on paper and used digital collaboration tools, such as Miro and Google Docs. Using an auto-ethnographic approach, we articulate our challenges and strategies throughout the process, providing useful insights about synchronous distributed collaboration.


title: Adaptive preferential sampling in phylodynamics

abstract: Longitudinal molecular data of rapidly evolving viruses and pathogens provide information about disease spread and complement traditional surveillance approaches based on case count data. The coalescent is used to model the genealogy that represents the sample ancestral relationships. The basic assumption is that coalescent events occur at a rate inversely proportional to the effective population size $N_{e}(t)$, a time-varying measure of genetic diversity. When the sampling process (collection of samples over time) depends on $N_{e}(t)$, the coalescent and the sampling processes can be jointly modeled to improve estimation of $N_{e}(t)$. Failing to do so can lead to bias due to model misspecification. However, the way that the sampling process depends on the effective population size may vary over time. We introduce an approach where the sampling process is modeled as an inhomogeneous Poisson process with rate equal to the product of $N_{e}(t)$ and a time-varying coefficient, making minimal assumptions on their functional shapes via Markov random field priors. We provide scalable algorithms for inference, show the model performance vis-a-vis alternative methods in a simulation study, and apply our model to SARS-CoV-2 sequences from Los Angeles and Santa Clara counties. The methodology is implemented and available in the R package adapref.


title: Baryon Acoustic Oscillations in the projected cross-correlation function between the eBOSS DR16 quasars and photometric galaxies from the DESI Legacy Imaging Surveys

abstract: We search for the Baryon Acoustic Oscillations in the projected cross-correlation function binned into transverse comoving radius between the SDSS-IV DR16 eBOSS quasars and a dense photometric sample of galaxies selected from the DESI Legacy Imaging Surveys. We estimate the density of the photometric sample of galaxies in this redshift range to be about 2900 deg$^{-2}$, which is deeper than the official DESI ELG selection, and the density of the spectroscopic sample is about 20 deg$^{-2}$. In order to mitigate the systematics related to the use of different imaging surveys close to the detection limit, we use a neural network approach that accounts for complex dependencies between the imaging attributes and the observed galaxy density. We find that we are limited by the depth of the imaging surveys which affects the density and purity of the photometric sample and its overlap in redshift with the quasar sample, which thus affects the performance of the method. When cross-correlating the photometric galaxies with quasars in $0.6 \leq z \leq 1.2$, the cross-correlation function can provide better constraints on the comoving angular distance, $D_{\rm M}$ (6\% precision) compared to the constraint on the spherically-averaged distance $D_{\rm V}$ (9\% precision) obtained from the auto-correlation. Although not yet competitive, this technique will benefit from the arrival of deeper photometric data from upcoming surveys which will enable it to go beyond the current limitations we have identified in this work.


title: Pion condensation in the early Universe at nonvanishing lepton flavor asymmetry and its gravitational wave signatures

abstract: We investigate the possible formation of a Bose-Einstein condensed phase of pions in the early Universe at nonvanishing values of lepton flavor asymmetries. A hadron resonance gas model with pion interactions, based on first-principle lattice QCD simulations at nonzero isospin density, is used to evaluate cosmic trajectories at various values of electron, muon, and tau lepton asymmetries that satisfy the available constraints on the total lepton asymmetry. The cosmic trajectory can pass through the pion condensed phase if the combined electron and muon asymmetry is sufficiently large: $|l_e + l_{\mu}| \gtrsim 0.1$, with little sensitivity to the difference $l_e - l_\mu$ between the individual flavor asymmetries. Future constraints on the values of the individual lepton flavor asymmetries will thus be able to either confirm or rule out the condensation of pions during the cosmic QCD epoch. We demonstrate that the pion condensed phase leaves an imprint both on the spectrum of primordial gravitational waves and on the mass distribution of primordial black holes at the QCD scale e.g. the black hole binary of recent LIGO event GW190521 can be formed in that phase.


title: A Differential Topological Model for Olfactory Learning and Representation

abstract: This thesis is designed to be a self-contained exposition of the neurobiological and mathematical aspects of sensory perception, memory, and learning with a bias towards olfaction. The final chapters introduce a new approach to modeling focusing more on the geometry of the system as opposed to element wise dynamics. Additionally, we construct an organism independent model for olfactory processing: something which is currently missing from the literature.


title: Can scale-dependent cosmology alleviate the $H_0$ tension?

abstract: Scale-dependence is a common feature to all effective models of quantum gravity. In this paper, a cosmological model based on the scale-dependent scenario of gravity is presented. It is argued that such models, where the scale-dependence appears as a correction to the classical $\Lambda$CDM evolution, have the potential of addressing the tensions between early and late time measurements of $H_0$. After defining criteria to parametrize this tension, we perform a numerical scan over the parameter space of the scale-dependent model, subject to these criteria. In this analysis, it is found that, indeed, the tension can be released.


title: Symbolic Powers of Cover Ideals of Graphs and Koszul Property

abstract: We show that attaching a whisker (or a pendant) at the vertices of a cycle cover of a graph results in a new graph with the following property: all symbolic powers of its cover ideal are Koszul or, equivalently, componentwise linear. This extends previous work where the whiskers were added to all vertices or to the vertices of a vertex cover of the graph.


title: Investigation of the N\'eel phase of the frustrated Heisenberg antiferromagnet by differentiable symmetric tensor networks

abstract: The recent progress in the optimization of two-dimensional tensor networks [H.-J. Liao, J.-G. Liu, L. Wang, and T. Xiang, Phys. Rev. X ${\bf 9}$, 031041 (2019)] based on automatic differentiation opened the way towards precise and fast optimization of such states and, in particular, infinite projected entangled-pair states (iPEPS) that constitute a generic-purpose ${\it Ansatz}$ for lattice problems governed by local Hamiltonians. In this work, we perform an extensive study of a paradigmatic model of frustrated magnetism, the $J_1-J_2$ Heisenberg antiferromagnet on the square lattice. By using advances in both optimization and subsequent data analysis, through finite correlation-length scaling, we report accurate estimations of the magnetization curve in the N\'eel phase for $J_2/J_1 \le 0.45$. The unrestricted iPEPS simulations reveal an $U(1)$ symmetric structure, which we identify and impose on tensors, resulting in a clean and consistent picture of antiferromagnetic order vanishing at the phase transition with a quantum paramagnet at $J_2/J_1 \approx 0.46(1)$. The present methodology can be extended beyond this model to study generic order-to-disorder transitions in magnetic systems.


title: Heterogeneous Coefficients, Control Variables, and Identification of Multiple Treatment Effects

abstract: Multidimensional heterogeneity and endogeneity are important features of models with multiple treatments. We consider a heterogeneous coefficients model where the outcome is a linear combination of dummy treatment variables, with each variable representing a different kind of treatment. We use control variables to give necessary and sufficient conditions for identification of average treatment effects. With mutually exclusive treatments we find that, provided the heterogeneous coefficients are mean independent from treatments given the controls, a simple identification condition is that the generalized propensity scores (Imbens, 2000) be bounded away from zero and that their sum be bounded away from one, with probability one. Our analysis extends to distributional and quantile treatment effects, as well as corresponding treatment effects on the treated. These results generalize the classical identification result of Rosenbaum and Rubin (1983) for binary treatments.


title: Dark Fluxes from Accreting Black Holes and Direct Detections

abstract: We discuss the possibility that accreting black hole systems can be sources for dark matter flux through several different mechanisms. We firstly discuss two types of systems: coronal thermal plasmas around supermassive black holes in active galactic nuclei (AGNs), and accretion disks of stellar-mass X-ray black hole binaries (BHBs). We explore how these black hole systems may produce keV light dark matter fluxes and find that the dark fluxes from those sources might be too weak to account for the current XENON1T excess. On the other hand, black holes can be good accelerators to accrete and boost heavy dark matter particles. If considering collisions or dark electromagnetism, those particles can then escape and reach the benchmark speed of 0.1c at the detector. We also extend the black hole mass region to primordial black holes (PBHs) and discuss the possibility of contributing to keV light dark flux via superradiance of PBHs.


title: Proposing a two-step Decision Support System (TPIS) based on Stacked ensemble classifier for early and low cost (step-1) and final (step-2) differential diagnosis of Mycobacterium Tuberculosis from non-tuberculosis Pneumonia

abstract: Background: Mycobacterium Tuberculosis (TB) is an infectious bacterial disease presenting similar symptoms to pneumonia; therefore, differentiating between TB and pneumonia is challenging. Therefore, the main aim of this study is proposing an automatic method for differential diagnosis of TB from Pneumonia. Methods: In this study, a two-step decision support system named TPIS is proposed for differential diagnosis of TB from pneumonia based on stacked ensemble classifiers. The first step of our proposed model aims at early diagnosis based on low-cost features including demographic characteristics and patient symptoms (including 18 features). TPIS second step makes the final decision based on the meta features extracted in the first step, the laboratory tests and chest radiography reports. This retrospective study considers 199 patient medical records for patients suffering from TB or pneumonia, which has been registered in a hospital in Arak, Iran. Results: Experimental results show that TPIS outperforms the compared machine learning methods for early differential diagnosis of pulmonary tuberculosis from pneumonia with AUC of 90.26 and accuracy of 91.37 and final decision making with AUC of 92.81 and accuracy of 93.89. Conclusions: The main advantage of early diagnosis is beginning the treatment procedure for confidently diagnosed patients as soon as possible and preventing latency in treatment. Therefore, early diagnosis reduces the maturation of late treatment of both diseases.


title: Approximation, characterization, and continuity of multivariate monotonic regression functions

abstract: We deal with monotonic regression of multivariate functions $f: Q \to \mathbb{R}$ on a compact rectangular domain $Q$ in $\mathbb{R}^d$, where monotonicity is understood in a generalized sense: as isotonicity in some coordinate directions and antitonicity in some other coordinate directions. As usual, the monotonic regression of a given function $f$ is the monotonic function $f^*$ that has the smallest (weighted) mean-squared distance from $f$. We establish a simple general approach to compute monotonic regression functions: namely, we show that the monotonic regression $f^*$ of a given function $f$ can be approximated arbitrarily well -- with simple bounds on the approximation error in both the $2$-norm and the $\infty$-norm -- by the monotonic regression $f_n^*$ of grid-constant functions $f_n$. We also establish the continuity of the monotonic regression $f^*$ of a continuous function $f$ along with an explicit averaging formula for $f^*$. And finally, we deal with generalized monotonic regression where the mean-squared distance from standard monotonic regression is replaced by more complex distance measures which arise, for instance, in maximum smoothed likelihood estimation. We will see that the solution of such generalized monotonic regression problems is simply given by the standard monotonic regression $f^*$.


title: Threshold selection for extremal index estimation

abstract: We propose a new threshold selection method for the nonparametric estimation of the extremal index of stochastic processes. The so-called discrepancy method was proposed as a data-driven smoothing tool for estimation of a probability density function. Now it is modified to select a threshold parameter of an extremal index estimator. To this end, a specific normalization of the discrepancy statistic based on the Cram\'{e}r-von Mises-Smirnov statistic $\omega^2$ is calculated by the $k$ largest order statistics instead of an entire sample. Its asymptotic distribution as $k\to\infty$ is proved to be the same as the $\omega^2$-distribution. The quantiles of the latter distribution are used as discrepancy values. The rate of convergence of an extremal index estimate coupled with the discrepancy method is derived. The discrepancy method is used as an automatic threshold selection for the intervals and $K-$gaps estimators and it may be applied to other estimators of the extremal index.


title: The \'etale-open topology and the stable fields conjecture

abstract: For an arbitrary field $K$ and $K$-variety $V$, we introduce the \'etale-open topology on the set $V(K)$ of $K$-points of $V$. This topology agrees with the Zariski topology, Euclidean topology, or valuation topology when $K$ is separably closed, real closed, or $p$-adically closed, respectively. Topological properties of the \'etale-open topology corresponds to algebraic properties of $K$. For example, the \'etale-open topology on $\mathbb{A}^1(K)$ is not discrete if and only if $K$ is large. As an application, we show that a large stable field is separably closed.


title: Quantum Coulomb glass on the Bethe lattice

abstract: We study the Coulomb glass emerging from the interplay of strong interactions and disorder in a model of spinless fermions on the Bethe lattice. In the infinite coordination number limit, strong interactions induce a metallic Coulomb glass phase with a pseudogap structure at the Fermi energy. Quantum and thermal fluctuations both melt this glass and induce a disordered quantum liquid phase. We combine self-consistent diagrammatic perturbation theory with continuous time quantum Monte-Carlo simulations to obtain the complete phase diagram of the electron glass, and to characterize its dynamical properties in the quantum liquid, as well as in the replica symmetry broken glassy phase. Tunneling spectra display an Efros-Shklovskii pseudogap upon decreasing temperatures, but the density of states remains finite at the Fermi energy due to residual quantum fluctuations. Our results bear relevance to the metallic glass phase observed in Si inversion layers.


title: A Tale of Planet Formation: From Dust to Planets

abstract: The characterization of exoplanets and their birth protoplanetary disks has enormously advanced in the last decade. Benefitting from that, our global understanding of the planet formation processes has been substantially improved. In this review, we first summarize the cutting-edge states of the exoplanet and disk observations. We further present a comprehensive panoptic view of modern core accretion planet formation scenarios, including dust growth and radial drift, planetesimal formation by the streaming instability, core growth by planetesimal accretion and pebble accretion. We discuss the key concepts and physical processes in each growth stage and elaborate on the connections between theoretical studies and observational revelations. Finally, we point out the critical questions and future directions of planet formation studies.


title: Absolute irreducibility of the binomial polynomials

abstract: In this paper we investigate the factorization behaviour of the binomial polynomials $\binom{x}{n} = \frac{x(x-1)\cdots (x-n+1)}{n!}$ and their powers in the ring of integer-valued polynomials $\operatorname{Int}(\mathbb{Z})$. While it is well-known that the binomial polynomials are irreducible elements in $\operatorname{Int}(\mathbb{Z})$, the factorization behaviour of their powers has not yet been fully understood. We fill this gap and show that the binomial polynomials are absolutely irreducible in $\operatorname{Int}(\mathbb{Z})$, that is, $\binom{x}{n}^m$ factors uniquely into irreducible elements in $\operatorname{Int}(\mathbb{Z})$ for all $m\in \mathbb{N}$. By reformulating the problem in terms of linear algebra and number theory, we show that the question can be reduced to determining the rank of, what we call, the valuation matrix of $n$. A main ingredient in computing this rank is the following number-theoretical result for which we also provide a proof: If $n>10$ and $n$, $n-1$, \ldots, $n-(k-1)$ are composite integers, then there exists a prime number $p > 2k$ that divides one of these integers.


title: Uniform pointwise estimates for ultraspherical polynomials

abstract: We prove pointwise bounds for two-parameter families of Jacobi polynomials. Our bounds imply estimates for a class of functions arising from the spectral analysis of distinguished Laplacians and sub-Laplacians on the unit sphere in arbitrary dimension, and are instrumental in the proof of sharp multiplier theorems for those operators.


title: Intelligent Reflecting Surface Aided Multi-User Communication: Capacity Region and Deployment Strategy

abstract: Intelligent reflecting surface (IRS) is a new promising technology that is able to reconfigure the wireless propagation channel via smart and passive signal reflection. In this paper, we investigate the capacity region of a two-user communication network with one access point (AP) aided by $M$ IRS elements for enhancing the user-AP channels, where the IRS incurs negligible delay, thus the user-AP channels via the IRS follow the classic discrete memoryless channel model. In particular, we consider two practical IRS deployment strategies that lead to different effective channels between the users and AP, namely, the distributed deployment where the $M$ elements form two IRSs, each deployed in the vicinity of one user, versus the centralized deployment where all the $M$ elements are deployed in the vicinity of the AP. First, we consider the uplink multiple-access channel (MAC) and derive the capacity/achievable rate regions for both deployment strategies under different multiple access schemes. It is shown that the centralized deployment generally outperforms the distributed deployment under symmetric channel setups in terms of achievable user rates. Next, we extend the results to the downlink broadcast channel (BC) by leveraging the celebrated uplink-downlink (or MAC-BC) duality framework, and show that the superior rate performance of centralized over distributed deployment also holds. Numerical results are presented that validate our analysis, and reveal new and useful insights for optimal IRS deployment in wireless networks.


title: How to fairly share a watermelon

abstract: Geometry, calculus and in particular integrals, are too often seen by young students as technical tools with no link to the reality. This fact generates into the students a loss of interest with a consequent removal of motivation in the study of such topics and more widely in pursuing scientific curricula. With this note we put to the fore a simple example of practical interest where the above concepts prove central; our aim is thus to motivate students and to reverse the dropout trend by proposing an introduction to the theory starting from practical applications. More precisely, we will show how using a mixture of geometry, calculus and integrals one can easily share a watermelon into regular slices with equal volume.


title: CLEANN: Accelerated Trojan Shield for Embedded Neural Networks

abstract: We propose CLEANN, the first end-to-end framework that enables online mitigation of Trojans for embedded Deep Neural Network (DNN) applications. A Trojan attack works by injecting a backdoor in the DNN while training; during inference, the Trojan can be activated by the specific backdoor trigger. What differentiates CLEANN from the prior work is its lightweight methodology which recovers the ground-truth class of Trojan samples without the need for labeled data, model retraining, or prior assumptions on the trigger or the attack. We leverage dictionary learning and sparse approximation to characterize the statistical behavior of benign data and identify Trojan triggers. CLEANN is devised based on algorithm/hardware co-design and is equipped with specialized hardware to enable efficient real-time execution on resource-constrained embedded platforms. Proof of concept evaluations on CLEANN for the state-of-the-art Neural Trojan attacks on visual benchmarks demonstrate its competitive advantage in terms of attack resiliency and execution overhead.


title: OnsagerNet: Learning Stable and Interpretable Dynamics using a Generalized Onsager Principle

abstract: We propose a systematic method for learning stable and physically interpretable dynamical models using sampled trajectory data from physical processes based on a generalized Onsager principle. The learned dynamics are autonomous ordinary differential equations parameterized by neural networks that retain clear physical structure information, such as free energy, diffusion, conservative motion and external forces. For high dimensional problems with a low dimensional slow manifold, an autoencoder with metric preserving regularization is introduced to find the low dimensional generalized coordinates on which we learn the generalized Onsager dynamics. Our method exhibits clear advantages over existing methods on benchmark problems for learning ordinary differential equations. We further apply this method to study Rayleigh-Benard convection and learn Lorenz-like low dimensional autonomous reduced order models that capture both qualitative and quantitative properties of the underlying dynamics. This forms a general approach to building reduced order models for forced dissipative systems.


title: On the origin of a rotating metal-poor stellar population in the Milky Way Nuclear Cluster

abstract: We explore the origin of a population of stars recently detected in the inner parsec of the Milky Way Nuclear Cluster (NC), which exhibit sub-solar metallicity and a higher rotation compared to the dominant population. Using state-of-the-art $N$-body simulations, we model the infall of a massive stellar system into the Galactic center, both of Galactic and extra-galactic origin. We show that the newly discovered population can either be the remnant of a massive star cluster formed a few kpc away from the Galactic center (Galactic scenario) or be accreted from a dwarf galaxy originally located at 10-100 kpc (extragalactic scenario) and that reached the Galactic center 3-5 Gyr ago. A comparison between our models and characteristic Galactocentric distance and metallicity distributions of Milky Way satellites and globular clusters favours the Galactic scenario. A comparison with clusters associated with the Enceladus-Sausage, Sequoia, Sagittarius and Canis Major structures suggests that the progenitor of the observed metal-poor substructure formed in-situ rather than being accreted.


title: MOA-2007-BLG-400 A Super-Jupiter Mass Planet Orbiting a Galactic BulgeK-dwarf Revealed by Keck Adaptive Optics Imaging

abstract: We present Keck/NIRC2 adaptive optics imaging of planetary microlensing event MOA-2007-BLG-400 that resolves the lens star system from the source. We find that the MOA-2007-BLG-400L planetary system consists of a $1.71\pm 0.27 M_{\rm Jup}$ planet orbiting a $0.69\pm 0.04M_{\odot}$ K-dwarf host star at a distance of $6.89\pm 0.77\,$kpc from the Sun. So, this planetary system probably resides in the Galactic bulge. The planet-host star projected separation is only weakly constrained due to the close-wide light curve degeneracy; the 2$\sigma$ projected separation range is 0.6--$7.2\,$AU. This host mass is at the top end of the range of masses predicted by a standard Bayesian analysis that assumes that all stars have an equal chance of hosting a star of the observed mass ratio. This and the similar result for event MOA-2013-BLG-220 suggests that more massive stars may be more likely to host planets with a mass ratio in the $0.002 < q < 0.004$ range that orbit beyond the snow line. These results also indicate the importance of host star mass measurements for exoplanets found by microlensing. The microlensing survey imaging data from NASA's Nancy Grace Roman Space Telescope (formerly WFIRST) mission will be doing mass measurements like this for a huge number of planetary events. This host lens is the highest contrast lens-source detected in microlensing mass measurement analysis (the lens being 10$\times$ fainter than the source). We present an improved method of calculating photometry and astrometry uncertainties based on the Jackknife method, which produces more accurate errors that are $\sim$$2.5 \times$ larger than previous estimates.

categories: astro-ph.EP astro-ph.GA astro-ph.IM astro-ph.SR


title: A PSF-based Approach to TESS High quality data Of Stellar clusters (PATHOS) -- III. Exploring the properties of young associations through their variables, dippers, and candidate exoplanets

abstract: Young associations in star forming regions are stellar systems that allow us to understand the mechanisms that characterise the stars in their early life and what happens around them. In particular, the analysis of the disks and of the exoplanets around young stars allows us to know the key processes that prevail in their evolution and understand the properties of the exoplanets orbiting older stars. The TESS mission is giving us the opportunity to extract and analyse the light curves of association members with high accuracy, but the crowding that affects these regions makes difficult the light curve extraction. In the PATHOS project, cutting-edge tools are used to extract high-precision light curves and identify variable stars and transiting exoplanets in open clusters and associations. In this work, I analysed the light curves of stars in five young ($\lesssim 10$ Myr) associations, searching for variables and candidate exoplanets. By using the rotational periods of the association members, I constrained the ages of the five stellar systems ($\sim 2$-$10$ Myr). I searched for dippers and I investigated the properties of the dust that forms the circumstellar disks. Finally, I searched for transiting signals, finding 6 strong candidate exoplanets. No candidates with radius $R_{\rm P}\lesssim 0.9 R_J$ have been detected, in agreement with the expectations. The frequency of giant planets resulted to be $\sim$2-3 %, higher than that expected for field stars ($\lesssim 1$ %); the low statistic makes this conclusion not strong, and new investigations on young objects are mandatory to confirm this result.

categories: astro-ph.SR astro-ph.EP astro-ph.GA astro-ph.IM


title: History entanglement entropy

abstract: A formalism is proposed to describe entangled quantum histories, and their entanglement entropy. We define a history vector, living in a tensor space with basis elements corresponding to the allowed histories, i.e. histories with nonvanishing amplitudes. The amplitudes are the components of the history vector, and contain the dynamical information. Probabilities of measurement sequences, and resulting collapse, are given by generalized Born rules: they are all expressed by means of projections and scalar products involving the history vector. Entangled history states are introduced, and a history density matrix is defined in terms of ensembles of history vectors. The corresponding history entropies (and history entanglement entropies for composite systems) are explicitly computed in two examples taken from quantum computation circuits.


title: An unusually low density ultra-short period super-Earth and three mini-Neptunes around the old star TOI-561

abstract: Based on HARPS-N radial velocities (RVs) and TESS photometry, we present a full characterisation of the planetary system orbiting the late G dwarf TOI-561. After the identification of three transiting candidates by TESS, we discovered two additional external planets from RV analysis. RVs cannot confirm the outer TESS transiting candidate, which would also make the system dynamically unstable. We demonstrate that the two transits initially associated with this candidate are instead due to single transits of the two planets discovered using RVs. The four planets orbiting TOI-561 include an ultra-short period (USP) super-Earth (TOI-561 b) with period $P_{\rm b} = 0.45$ d, mass $M_{\rm b} =1.59 \pm 0.36$ M$_\oplus$ and radius $R_{\rm b}=1.42 \pm 0.07$ R$_\oplus$, and three mini-Neptunes: TOI-561 c, with $P_{\rm c} = 10.78$ d, $M_{\rm c} = 5.40 \pm 0.98$ M$_\oplus$, $R_{\rm c}= 2.88 \pm 0.09$ R$_\oplus$; TOI-561 d, with $P_{\rm d} = 25.6$ d, $M_{\rm d} = 11.9 \pm 1.3$ M$_\oplus$, $R_{\rm d} = 2.53 \pm 0.13$ R$_\oplus$; and TOI-561 e, with $P_{\rm e} = 77.2$ d, $M_{\rm e} = 16.0 \pm 2.3$ M$_\oplus$, $R_{\rm e} = 2.67 \pm 0.11$ R$_\oplus$. Having a density of $3.0 \pm 0.8$ g cm$^{-3}$, TOI-561 b is the lowest density USP planet known to date. Our N-body simulations confirm the stability of the system and predict a strong, anti-correlated, long-term transit time variation signal between planets d and e. The unusual density of the inner super-Earth and the dynamical interactions between the outer planets make TOI-561 an interesting follow-up target.


title: Radial variation of the stellar mass functions in the globular clusters M15 and M30: clues of a non-standard IMF?

abstract: We exploit a combination of high-resolution Hubble Space Telescope and wide-field ESO-VLT observations to study the slope of the global mass function (alphaG) and its radial variation (alpha(r)) in the two dense, massive and post core-collapse globular clusters M15 and M30. The available data-set samples the clusters' Main Sequence down to 0.2 Msun and the photometric completeness allows the study of the mass function between 0.40 Msun and 0.75 Msun from the central regions out to their tidal radii. We find that both clusters show a very similar variation in alpha(r) as a function of clustercentric distance. They both exhibit a very steep variation in alpha(r) in the central regions, which then attains almost constant values in the outskirts. Such a behavior can be interpreted as the result of long-term dynamical evolution of the systems driven by mass-segregation and mass-loss processes. We compare these results with a set of direct N-body simulations and find that they are only able to reproduce the observed values of alpha(r) and alphaG at dynamical ages (t/trh) significantly larger than those derived from the observed properties of both clusters. We investigate possible physical mechanisms responsible for such a discrepancy and argue that both clusters might be born with a non-standard (flatter/bottom-lighter) initial mass function.


title: First L band detection of hot exozodiacal dust with VLTI/MATISSE

abstract: For the first time we observed the emission of hot exozodiacal dust in L band. We used the new instrument MATISSE at the Very Large Telescope Interferometer to detect the hot dust around $\kappa$ Tuc with a significance of 3$\sigma$ to 6$\sigma$ at wavelengths between 3.37 and 3.85 $\mu$m and a dust-to-star flux ratio of 5 to 7 %. We modelled the spectral energy distribution based on the new L band data alone and in combination with H band data published previously. In all cases we find 0.58 $\mu$m grains of amorphous carbon to fit the $\kappa$ Tuc observations the best, however, also nanometre or micrometre grains and other carbons or silicates reproduce the observations well. Since the H band data revealed a temporal variability, while our L band data were taken at a different epoch, we combine them in different ways. Depending on the approach, the best fits are obtained for a narrow dust ring at a stellar distance in the 0.1 to 0.29 au range and thus with a temperature between 940 and 1430 K. Within the 1$\sigma$ uncertainty dust location and temperature are confined to 0.032 - 1.18 au and 600 - 2000 K.


title: Revealing the Formation of the Milky Way Nuclear Star Cluster via Chemo-Dynamical Modeling

abstract: The Milky Way nuclear star cluster (MW NSC) has been used as a template to understand the origin and evolution of galactic nuclei and the interaction of nuclear star clusters with supermassive black holes. It is the only nuclear star cluster with a supermassive black hole where we can resolve individual stars to measure their kinematics and metal abundance to reconstruct its formation history. Here, we present results of the first chemo-dynamical model of the inner 1 pc of the MW NSC using metallicity and radial velocity data from the KMOS spectrograph on the Very Large Telescope. We find evidence for two kinematically and chemically distinct components in this region. The majority of the stars belong to a previously known super-solar metallicity component with a rotation axis perpendicular to the Galactic plane. However, we identify a new kinematically distinct sub-solar metallicity component which contains about 7\% of the stars and appears to be rotating faster than the main component with a rotation axis that may be misaligned. This second component may be evidence for an infalling star cluster or remnants of a dwarf galaxy, merging with the MW NSC. These measurements show that the combination of chemical abundances with kinematics is a promising method to directly study the MW NSC's origin and evolution.


title: Non-Fermi Liquids in Conducting 2D Networks

abstract: We explore the physics of novel fermion liquids emerging from conducting networks, where 1D metallic wires form a periodic 2D superstructure. Such structure naturally appears in marginally-twisted bilayer graphenes, moire transition metal dichalcogenides, and also in some charge-density wave materials. For these network systems, we theoretically show that a remarkably wide variety of new non-Fermi liquids emerge and that these non-Fermi liquids can be classified by the characteristics of the junctions in networks. Using this, we calculate the electric conductivity of the non-Fermi liquids as a function of temperature, which show markedly different scaling behaviors than a regular 2D Fermi liquid.


title: Geometric entanglement in integer quantum Hall states

abstract: We study the quantum entanglement structure of integer quantum Hall states via the reduced density matrix of spatial subregions. In particular, we examine the eigenstates, spectrum and entanglement entropy (EE) of the density matrix for various ground and excited states, with or without mass anisotropy. We focus on an important class of regions that contain sharp corners or cusps, leading to a geometric angle-dependent contribution to the EE. We unravel surprising relations by comparing this corner term at different fillings. We further find that the corner term, when properly normalized, has nearly the same angle dependence as numerous conformal field theories (CFTs) in two spatial dimensions, which hints at a broader structure. In fact, the Hall corner term is found to obey bounds that were previously obtained for CFTs. In addition, the low-lying entanglement spectrum and the corresponding eigenfunctions reveal "excitations" localized near corners. Finally, we present an outlook for fractional quantum Hall states.

categories: cond-mat.str-el cond-mat.mes-hall hep-th quant-ph


title: On the construction of convolution-like operators associated with multidimensional diffusion processes

abstract: When is it possible to interpret a given Markov process as a L\'evy-like process? Since the class of L\'evy processes can be defined by the relation between transition probabilities and convolutions, the answer to this question lies in the existence of a convolution-like operator satisfying the same relation with the transition probabilities of the process. It is known that the so-called Sturm-Liouville convolutions have the desired properties and therefore the question above has a positive answer for a certain class of one-dimensional diffusions. However, more general processes have never been systematically treated in the literature. This study addresses this gap by considering the general problem of constructing a convolution-like operator for a given strong Feller process on a general locally compact metric space. Both necessary and sufficient conditions for the existence of such convolution-like structures are determined, which reveal a connection between the answer to the above question and certain analytical and geometrical properties of the eigenfunctions of the transition semigroup. The case of reflected Brownian motions on bounded domains of R d and compact Riemannian manifolds is considered in greater detail: various special cases are analysed, and a general discussion on the existence of appropriate convolution-like structures is presented.


title: The Complex Upper HR Diagram

abstract: Several decades of observations of the most massive and most luminous stars have revealed a complex upper HR Diagram, shaped by mass loss, and inhabited by a variety of evolved stars exhibiting the consequences of their mass loss histories. This review presents a brief historical overview of the HR Diagram for massive stars, highlighting some of the primary discoveries and results from their observation in nearby galaxies. The chapters in this volume include reviews of our current understanding of different groups of evolved massive stars, all losing mass and in different stages of their evolution; the Luminous Blue Variables (LBVs), B[e] supergiants, the warm hypergiants, Wolf-Rayet stars, and the population of OB stars and supergiants in the Magellanic Clouds.


title: Radiation-Driven Stellar Eruptions

abstract: Very massive stars occasionally expel material in colossal eruptions, driven by continuum radiation pressure rather than blast waves. Some of them rival supernovae in total radiative output, and the mass loss is crucial for subsequent evolution. Some are supernova impostors, including SN precursor outbursts, while others are true SN events shrouded by material that was ejected earlier. Luminous Blue Variable stars (LBV's) are traditionally cited in relation with giant eruptions, though this connection is not well established. After four decades of research, the fundamental causes of giant eruptions and LBV events remain elusive. This review outlines the basic relevant physics, with a brief summary of essential observational facts. Reasons are described for the spectrum and emergent radiation temperature of an opaque outflow. Many proposed mechanisms are noted for instabilities in the star's photosphere, in its iron opacity peak zones, and in its central region. Some of the remarks and conjectures here have not yet become familiar in the published literature.


title: Conservative interpolation of aeroacoustic sources in a hybrid workflow applied to fan

abstract: In low Mach number aeroacoustics, the well known disparity of scales makes it possible to apply efficient hybrid simulation models using different meshes for flow and acoustics, which leads to a powerful computational procedure. Our study applies the hybrid workflow to the computationally efficient perturbed convective wave equation with only one scalar unknown, the acoustic velocity potential. The workflow of this aeroacoustic approach is based on three steps: 1. perform unsteady incompressible flow computations on a sub-domain; 2. compute the acoustic sources; 3. simulate the acoustic field using a mesh specifically suited for computational aeroacoustics. In general, hybrid aeroacoustic methods seek for robust and conservative mesh-to-mesh transformation of the aeroacoustic sources while high computational efficiency is ensured. In this paper, we investigate the accuracy of a cell-centroid based conservative interpolation scheme compared to the more accurate cut-volume cell approach and their application to the computation of rotating systems, namely an axial fan. The capability and robustness of the cut-volume cell interpolation in a hybrid workflow on different meshes are investigated by a grid convergence study. The results of the acoustic simulation are in good agreement with measurements thus demonstrating the applicability of the conservative cut-volume cell interpolation to rotating systems.


title: Measuring the mixing scale of the ISM within nearby spiral galaxies

abstract: The spatial distribution of metals reflects, and can be used to constrain, the processes of chemical enrichment and mixing. Using PHANGS-MUSE optical integral field spectroscopy, we measure the gas phase oxygen abundances (metallicities) across 7,138 HII regions in a sample of eight nearby disc galaxies. In Paper I (Kreckel et al. 2019) we measure and report linear radial gradients in the metallicities of each galaxy, and qualitatively searched for azimuthal abundance variations. Here, we examine the two-dimensional variation in abundances once the radial gradient is subtracted, Delta(O/H), in order to quantify the homogeneity of the metal distribution and to measure the mixing scale over which HII region metallicities are correlated. We observe low (0.03--0.05 dex) scatter in Delta(O/H) globally in all galaxies, with significantly lower (0.02--0.03 dex) scatter on small (<600 pc) spatial scales. This is consistent with the measurement uncertainties, and implies the two-dimensional metallicity distribution is highly correlated on scales of <600 pc. We compute the two point correlation function for metals in the disc in order to quantify the scale lengths associated with the observed homogeneity. This mixing scale is observed to correlate better with the local gas velocity dispersion (of both cold and ionized gas) than with the star formation rate. Selecting only HII regions with enhanced abundances relative to a linear radial gradient, we do not observe increased homogeneity on small scales. This suggests that the observed homogeneity is driven by the mixing introducing material from large scales rather than by pollution from recent and on-going star formation.


title: Age dissection of the vertical breathing motions in Gaia DR2: evidence for spiral driving

abstract: Gaia DR2 has revealed breathing motions in the Milky Way, with stars on both sides of the Galactic mid-plane moving coherently towards or away from it. The generating mechanism of these breathing motions is thought to be spiral density waves. Here we test this hypothesis. Using a self-consistent, high-resolution simulation with star formation, and which hosts prominent spirals, we first study the signatures of breathing motions excited by spirals. In the model, the breathing motions induced by the spiral structure have an increasing amplitude with distance from the mid-plane, pointing to an internal cause for them. We then show that, at fixed height, the breathing motion amplitude decreases with age. Next, we investigate the signature of the breathing motions in the Gaia DR2 dataset. We demonstrate that, at the location with a consistently large breathing motion, the corresponding amplitude increases monotonically with distance from the mid-plane, in agreement with the model. Furthermore, we show that at the same location, the breathing motion amplitude decreases with age, again similar to what we find in the model. This strengthens the case that the observed breathing motions are driven by spiral density waves.


title: Anderson Molecules

abstract: Atoms can form molecules if they attract each other. Here, we show that atoms are also able to form bound states not due to the attractive interaction but because of destructive interference. If the interaction potential changes in a disordered way with a change of the distance between two atoms, Anderson localization can lead to the formation of exponentially localized bound states. While disordered interaction potentials do not exist in nature, we show that they can be created by means of random modulation in time of the strength of the original interaction potential between atoms and objects that we dub Anderson molecules can be realized in the laboratory.

categories: cond-mat.quant-gas cond-mat.mes-hall physics.atm-clus quant-ph


title: End-to-End Deep Learning Model for Cardiac Cycle Synchronization from Multi-View Angiographic Sequences

abstract: Dynamic reconstructions (3D+T) of coronary arteries could give important perfusion details to clinicians. Temporal matching of the different views, which may not be acquired simultaneously, is a prerequisite for an accurate stereo-matching of the coronary segments. In this paper, we show how a neural network can be trained from angiographic sequences to synchronize different views during the cardiac cycle using raw x-ray angiography videos exclusively. First, we train a neural network model with angiographic sequences to extract features describing the progression of the cardiac cycle. Then, we compute the distance between the feature vectors of every frame from the first view with those from the second view to generate distance maps that display stripe patterns. Using pathfinding, we extract the best temporally coherent associations between each frame of both videos. Finally, we compare the synchronized frames of an evaluation set with the ECG signals to show an alignment with 96.04% accuracy.


title: SlackQ : Approaching the Qubit Mapping Problem with A Slack-aware Swap Insertion Scheme

abstract: The rapid progress of physical implementation of quantum computers paved the way for the design of tools to help users write quantum programs for any given quantum device. The physical constraints inherent in current NISQ architectures prevent most quantum algorithms from being directly executed on quantum devices. To enable two-qubit gates in the algorithm, existing works focus on inserting SWAP gates to dynamically remap logical qubits to physical qubits. However, their schemes lack consideration of the execution time of generated quantum circuits. In this work, we propose a slack-aware SWAP insertion scheme for the qubit mapping problem in the NISQ era. Our experiments show performance improvement by up to 2.36X at maximum, by 1.62X on average, over 106 representative benchmarks from RevLib, IBM Qiskit , and ScaffCC.


title: The peculiar Ca-rich SN 2019ehk: Evidence for a Type IIb core-collapse supernova from a low mass stripped progenitor

abstract: The nature of the peculiar `Ca-rich' SN 2019ehk in the nearby galaxy M100 remains unclear. Its origin has been debated as either a stripped core-collapse supernova or a thermonuclear helium detonation event. Here, we present very late-time photometry of the transient obtained with the Keck I telescope at $\approx280$ days from peak light. Using the photometry to perform accurate flux calibration of a contemporaneous nebular phase spectrum, we measure an [O I] luminosity of $(0.19-1.08)\times10^{38}$ erg s$^{-1}$ and [Ca II] luminosity of $(2.7-15.6)\times10^{38}$ erg s$^{-1}$ over the range of the uncertain extinction along the line of sight. We use these measurements to derive lower limits on the synthesized oxygen mass of $\approx0.004-0.069$ M$_\odot$. The oxygen mass is a sensitive tracer of the progenitor mass for core-collapse supernovae, and our estimate is consistent with explosions of very low mass CO cores of $1.45-1.5$ M$_\odot$, corresponding to He core masses of $\approx1.8-2.0$ M$_\odot$. We present high quality peak light optical spectra of the transient and highlight features of hydrogen in both the early (`flash') and photospheric phase spectra, that suggest the presence of $\gtrsim0.02$ M$_\odot$ of hydrogen in the progenitor at the time of explosion. The presence of H, together with the large [Ca II]/[O I] ratio ($\approx10-15$) in the nebular phase is consistent with SN 2019ehk being a Type IIb core-collapse supernova from a stripped low mass ($\approx9-9.5$ M$_\odot$) progenitor, similar to the Ca-rich SN IIb iPTF 15eqv. These results provide evidence for a likely class of `Ca-rich' core-collapse supernovae from stripped low mass progenitors in star forming environments, distinct from the thermonuclear Ca-rich gap transients in old environments.


title: Pilaster: A Collection of Citation Metadata Extracted From Publications on Visualization for the Digital Humanities

abstract: In this paper, we present Pilaster (https://visusal.github.io/pilaster/), a collection of citation metadata extracted from publications in visualization for the digital humanities. The collection is generated from a seed set of relevant publications from which we extracted cited works, including journal and conference papers, books, theses, or blog posts, among other resources. The main aim of this work revolves around three main points: first, the collection may serve as an entry point to the discipline for digital humanists and visualization scholars without previous experience in the field. Second, Pilaster can be regarded as a meeting point for more established visualization or humanities scholars seeking to collaborate in the development of novel research ideas and related visualization design studies in the context of the humanities. Third, and given the large amount of visualization design spaces that were captured, we believe the dataset has the potential to become the starting point for future studies aimed at understanding the particularities of problem-driven visualization research in this and other contexts.


title: An optimal home energy management system for modulating heat pumps and photovoltaic systems

abstract: Efficient residential sector coupling plays a key role in supporting the energy transition. In this study, we analyze the structural properties associated with the optimal control of a home energy management system and the effects of common technological configurations and objectives. We conduct this study by modeling a representative building with a modulating air-sourced heat pump, a photovoltaic (PV) system, a battery, and thermal storage systems for floor heating and hot-water supply. In addition, we allow grid feed-in by assuming fixed feed-in tariffs and consider user comfort. In our numerical analysis, we find that the battery, naturally, is the essential building block for improving self-sufficiency. However, in order to use the PV surplus efficiently grid feed-in is necessary. The commonly considered objective of maximizing self-consumption is not economically viable under the given tariff structure; however, close-to-optimal performance and significant reduction in solution times can be achieved by maximizing self-sufficiency. Based on optimal control and considering seasonal effects, the dominant order of PV distribution and the target states of charge of the storage systems can be derived. Using a rolling horizon approach, the solution time can be reduced to less than 1 min (achieving a time resolution of 1 h per year). By evaluating the value of information, we find that the common value of 24 h for the prediction and control horizons results in unintended but avoidable end-of-horizon effects. Our input data and mixed-integer linear model developed using the Julia JuMP programming language are available in an open-source manner.


title: Whitham modulation theory for generalized Whitham equations and a general criterion for modulational instability

abstract: The Whitham equation was proposed as a model for surface water waves that combines the quadratic flux nonlinearity $f(u) = \tfrac{1}{2}u^2$ of the Korteweg-de Vries equation and the full linear dispersion relation $\Omega(k) = \sqrt{k\tanh k}$ of uni-directional gravity water waves in suitably scaled variables. This paper proposes and analyzes a generalization of Whitham's model to unidirectional nonlinear wave equations consisting of a general nonlinear flux function $f(u)$ and a general linear dispersion relation $\Omega(k)$. Assuming the existence of periodic traveling wave solutions to this generalized Whitham equation, their slow modulations are studied in the context of Whitham modulation theory. A multiple scales calculation yields the modulation equations, a system of three conservation laws that describe the slow evolution of the periodic traveling wave's wavenumber, amplitude, and mean. In the weakly nonlinear limit, explicit, simple criteria in terms of general $f(u)$ and $\Omega(k)$ establishing the strict hyperbolicity and genuine nonlinearity of the modulation equations are determined. This result is interpreted as a generalized Lighthill-Whitham criterion for modulational instability.


title: Risk-Averse Optimization for Resilience Enhancement of Complex Engineering Systems under Uncertainties

abstract: With the growth of complexity and extent, large scale interconnected network systems, e.g. transportation networks or infrastructure networks, become more vulnerable towards external disruptions. Hence, managing potential disruptive events during the design, operating, and recovery phase of an engineered system therefore improving the system's resilience is an important yet challenging task. In order to ensure system resilience after the occurrence of failure events, this study proposes a mixed-integer linear programming (MILP) based restoration framework using heterogeneous dispatchable agents. Scenario-based stochastic optimization (SO) technique is adopted to deal with the inherent uncertainties imposed on the recovery process from nature. Moreover, different from conventional SO using deterministic equivalent formulations, additional risk measure is implemented for this study because of the temporal sparsity of the decision making in applications such as the recovery from extreme events. The resulting restoration framework involves a large-scale MILP problem and thus an adequate decomposition technique, i.e. modified Lagrangian dual decomposition, is also employed in order to achieve tractable computational complexity. Case study results based on the IEEE 37-bus test feeder demonstrate the benefits of using the proposed framework for resilience improvement as well as the advantages of adopting SO formulations.


title: Grassmannian-parameterized solutions to direct-sum polygon and simplex equations

abstract: We consider polygon and simplex equations, of which the simplest nontrivial examples are pentagon (5-gon) and Yang--Baxter (2-simplex), respectively. We examine the general structure of (2n+1)-gon and 2n-simplex equations in direct sums of vector spaces. Then we provide a construction for their solutions, parameterized by elements of the Grassmannian Gr(n+1,2n+1).


title: Running Neural Networks on the NIC

abstract: In this paper we show that the data plane of commodity programmable (Network Interface Cards) NICs can run neural network inference tasks required by packet monitoring applications, with low overhead. This is particularly important as the data transfer costs to the host system and dedicated machine learning accelerators, e.g., GPUs, can be more expensive than the processing task itself. We design and implement our system -- N3IC -- on two different NICs and we show that it can greatly benefit three different network monitoring use cases that require machine learning inference as first-class-primitive. N3IC can perform inference for millions of network flows per second, while forwarding traffic at 40Gb/s. Compared to an equivalent solution implemented on a general purpose CPU, N3IC can provide 100x lower processing latency, with 1.5x increase in throughput.


title: Efficient simulation of moire materials using the density matrix renormalization group

abstract: We present an infinite density-matrix renormalization group (DMRG) study of an interacting continuum model of twisted bilayer graphene (tBLG) near the magic angle. Because of the long-range Coulomb interaction and the large number of orbital degrees of freedom, tBLG is difficult to study with standard DMRG techniques -- even constructing and storing the Hamiltonian already poses a major challenge. To overcome these difficulties, we use a recently developed compression procedure to obtain a matrix product operator representation of the interacting tBLG Hamiltonian which we show is both efficient and accurate even when including the spin, valley and orbital degrees of freedom. To benchmark our approach, we focus mainly on the spinless, single-valley version of the problem where, at half-filling, we find that the ground state is a nematic semimetal. Remarkably, we find that the ground state is essentially a k-space Slater determinant, so that Hartree-Fock and DMRG give virtually identical results for this problem. Our results show that the effects of long-range interactions in magic angle graphene can be efficiently simulated with DMRG, and opens up a new route for numerically studying strong correlation physics in spinful, two-valley tBLG, and other moire materials, in future work.


title: "Reduction of Monetary Cost in Cloud Storage System by Using Extended Strict Timed Causal Consistency"

abstract: Cloud storage systems have been introduced to provide a scalable, secure, reliable, and highly available data storage environment for the organizations and end-users. Therefore, the service provider should grow in a geographical extent. Consequently, extensive storage service provision requires a replication mechanism. Replication imposes many costs on the cloud storage, including the synchronization, communications, storage, etc., costs among the replicas. Moreover, the synchronization process among replicas is a major challenge in cloud storage. Therefore, consistency can be defined as the coordination among the replicas. In this paper, we propose an extension to the strict timed causal consistency by adding the considerations for the monetary costs and the number of violations in the cloud storage systems and call it the extended strict timed causal consistency. Our proposed supports monotonic read, read your write, monotonic write, and write follow read, models by taking into account the causal relations between users' operations, at the client-side. Besides, it supports timed causal at the server-side. We employed the Cassandra cloud database that supports various consistencies such as all, one, quorum, etc. Our method performs better in reducing staleness rate, the severity of violations, and monetary cost in comparison with all, one, quorum, and causal.


title: Symmetry of electron bands in graphene: (nearly) free electron vs. tight-binding

abstract: We present the symmetry labelling of all electron bands in graphene obtained by combining numerical band calculations and analytical analysis based on group theory. The latter was performed both in the framework of the (nearly) free electron model, or in the framework of the tight-binding model. The predictions about relative positions of the bands which can be made on the basis of each of the models just using the group theory (and additional simple qualitative arguments, if necessary) are complimentary.


title: Nature vs. nurture: a Bayesian framework for assessing apparent correlations between planetary orbital properties and stellar ages

abstract: Many exoplanets have orbital characteristics quite different from those seen in our own solar system, including planets locked in orbital resonances and planets on orbits that are elliptical or highly inclined from their host star's spin axis. It is debated whether the wide variety in system architecture is primarily due to differences in formation conditions (nature) or due to evolution over time (nurture). Identifying trends between planetary and stellar properties, including stellar age, can help distinguish between these competing theories and offer insights as to how planets form and evolve. However, it can be challenging to determine whether observed trends between planetary properties and stellar age are driven by the age of the system -- pointing to evolution over time being an important factor -- or other parameters to which the age may be related, such as stellar mass or stellar temperature. The situation is complicated further by the possibilities of selection biases, small number statistics, uncertainties in stellar age, and orbital evolution timescales that are typically much shorter than the range of observed ages. Here we develop a Bayesian statistical framework to assess the robustness of such observed correlations and to determine whether they are indeed due to evolutionary processes, are more likely to reflect different formation scenarios, or are merely coincidental. We apply this framework to reported trends between stellar age and 2:1 orbital resonances, spin-orbit misalignments, and hot Jupiters' orbital eccentricities. We find strong support for the nurture hypothesis only in the final case.


title: Recent Trends in the Use of Deep Learning Models for Grammar Error Handling

abstract: Grammar error handling (GEH) is an important topic in natural language processing (NLP). GEH includes both grammar error detection and grammar error correction. Recent advances in computation systems have promoted the use of deep learning (DL) models for NLP problems such as GEH. In this survey we focus on two main DL approaches for GEH: neural machine translation models and editor models. We describe the three main stages of the pipeline for these models: data preparation, training, and inference. Additionally, we discuss different techniques to improve the performance of these models at each stage of the pipeline. We compare the performance of different models and conclude with proposed future directions.


title: The boundary density profile of a Coulomb droplet. Freezing at the edge

abstract: We revisit the problem of computing the boundary density profile of a droplet of two-dimensional one-component plasma (2D OCP) with logarithmic interaction between particles in a confining harmonic potential. At a sufficiently low temperature but still in the liquid phase, the density exhibits oscillations as a function of the distance to the boundary of the droplet. We obtain the density profile numerically using Monte-Carlo simulations of the 2D OCP. We argue that the decay and period of those oscillations can be explained within a picture of the Wigner crystallization near the boundary, where the crystal is gradually melted with the increasing distance to the boundary.


title: Self-similar potentials in quantum mechanics and coherent states

abstract: A brief description of the relations between the factorization method in quantum mechanics, self-similar potentials, integrable systems and the theory of special functions is given. New coherent states of the harmonic oscillator related to the Fourier transformation are constructed.

categories: quant-ph cond-mat.stat-mech hep-th math-ph math.CA math.MP


title: Energy-Resolved Information Scrambling in Energy-Space Lattices

abstract: Weakly interacting Fermi gases simulate spin-lattices in energy-space, offering a rich platform for investigating information spreading and spin coherence in a large many-body quantum system. We show that the collective spin vector can be determined as a function of energy from the measured spin density, enabling general energy-space resolved protocols. We measure an out-of-time-order correlation function in this system and observe the energy dependence of the many-body coherence.


title: Bootstrap p-values reduce type 1 error of the robust rank-order test of difference in medians

abstract: The robust rank-order test (Fligner and Policello, 1981) was designed as an improvement of the non-parametric Wilcoxon-Mann-Whitney U-test to be more appropriate when the samples being compared have unequal variance. However, it tends to be excessively liberal when the samples are asymmetric. This is likely because the test statistic is assumed to have a standard normal distribution for sample sizes > 12. This work proposes an on-the-fly method to obtain the distribution of the test statistic from which the critical/p-value may be computed directly. The method of likelihood maximization is used to estimate the parameters of the parent distributions of the samples being compared. Using these estimated populations, the null distribution of the test statistic is obtained by the Monte-Carlo method. Simulations are performed to compare the proposed method with that of standard normal approximation of the test statistic. For small sample sizes (<= 20), the Monte-Carlo method outperforms the normal approximation method. This is especially true for low values of significance levels (< 5%). Additionally, when the smaller sample has the larger standard deviation, the Monte-Carlo method outperforms the normal approximation method even for large sample sizes (= 40/60). The two methods do not differ in power. Finally, a Monte-Carlo sample size of 10^4 is found to be sufficient to obtain the aforementioned relative improvements in performance. Thus, the results of this study pave the way for development of a toolbox to perform the robust rank-order test in a distribution-free manner.


title: The auroral radio emission of the magnetic B-type star rho OphC

abstract: The non-thermal radio emission of main-sequence early-type stars is a signature of stellar magnetism. We present multi-wavelength (1.6-16.7 GHz) ATCA measurements of the early-type magnetic star rho OphC, which is a flat-spectrum non-thermal radio source. The rho OphC radio emission is partially circularly polarized with a steep spectral dependence: the fraction of polarized emission is about 60% at the lowest frequency sub-band (1.6 GHz) while is undetected at 16.7 GHz. This is clear evidence of coherent Auroral Radio Emission (ARE) from the rho OphC magnetosphere. Interestingly, the detection of the rho OphC's ARE is not related to a peculiar rotational phase. This is a consequence of the stellar geometry, which makes the strongly anisotropic radiation beam of the amplified radiation always pointed towards Earth. The circular polarization sign evidences mainly amplification of the ordinary mode of the electromagnetic wave, consistent with a maser amplification occurring within dense regions. This is indirect evidence of the plasma evaporation from the polar caps, a phenomenon responsible for the thermal X-ray aurorae. rho OphC is not the first early-type magnetic star showing the O-mode dominated ARE but is the first star with the ARE always on view.


title: Hybrid reaction-diffusion and clock-and-wavefront model for the arrest of oscillations in the somitogenesis segmentation clock

abstract: The clock and wavefront paradigm is arguably the most widely accepted model for explaining the embryonic process of somitogenesis. According to this model, somitogenesis is based upon the interaction between a genetic oscillator, known as segmentation clock, and a differentiation wavefront, which provides the positional information indicating where each pair of somites is formed. Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation. Cotterell et al. called it a progressive oscillatory reaction-diffusion (PORD) model. In the Meinhardt- PORD model, somitogenesis is driven by short-range interactions and the posterior movement of the front is a local, emergent phenomenon, which is not controlled by global positional information. With this model, it is possible to explain some experimental observations that are incompatible with the clock and wavefront model. However the Meinhardt-PORD model has some important disadvantages of its own. Namely, it is quite sensitive to fluctuations and depends on very specific initial conditions (which are not biologically realistic). In this work, we propose an equivalent Meinhardt-PORD model, and then amend it to couple it with a wavefront consisting of a receding morphogen gradient. By doing so, we get a hybrid model between the Meinhardt-PORD and the clock-and-wavefront ones, which overcomes most of the deficiencies of the two originating models.


title: LFGCN: Levitating over Graphs with Levy Flights

abstract: Due to high utility in many applications, from social networks to blockchain to power grids, deep learning on non-Euclidean objects such as graphs and manifolds, coined Geometric Deep Learning (GDL), continues to gain an ever increasing interest. We propose a new L\'evy Flights Graph Convolutional Networks (LFGCN) method for semi-supervised learning, which casts the L\'evy Flights into random walks on graphs and, as a result, allows both to accurately account for the intrinsic graph topology and to substantially improve classification performance, especially for heterogeneous graphs. Furthermore, we propose a new preferential P-DropEdge method based on the Girvan-Newman argument. That is, in contrast to uniform removing of edges as in DropEdge, following the Girvan-Newman algorithm, we detect network periphery structures using information on edge betweenness and then remove edges according to their betweenness centrality. Our experimental results on semi-supervised node classification tasks demonstrate that the LFGCN coupled with P-DropEdge accelerates the training task, increases stability and further improves predictive accuracy of learned graph topology structure. Finally, in our case studies we bring the machinery of LFGCN and other deep networks tools to analysis of power grid networks - the area where the utility of GDL remains untapped.


title: Impact of the 50 Hz harmonics on the beam evolution of the Large Hadron Collider

abstract: Harmonics of the mains frequency (50 Hz) have been systematically observed in the transverse beam spectrum of the Large Hadron Collider (LHC) since the start of its operation in the form of dipolar excitations. In the presence of strong non-linearities such as beam-beam interactions, as many of these power supply ripple tones reside in the vicinity of the betatron tune they can increase the tune diffusion of the particles in the distribution, leading to proton losses and eventually to a significant reduction of the beam lifetime. The aim of this paper is to determine whether the 50~Hz harmonics have an impact on the beam performance of the LHC. A quantitative characterization of the ripple spectrum present in the operation of the accelerator, together with an understanding of its source is an essential ingredient to also evaluate the impact of the 50~Hz harmonics on the future upgrade of the LHC, the High Luminosity LHC (HL-LHC). To this end, simulations with the single-particle tracking code, SixTrack, are employed including a realistic ripple spectrum as extracted from experimental observations to quantify the impact of such effects in terms of tune diffusion, Dynamic Aperture and beam lifetime. The methods and results of the tracking studies are reported and discussed in this paper.


title: Scalar Perturbations of a Single-Horizon Regular Black Hole

abstract: We investigate the massless scalar field perturbations, including the quasinormal mode spectrum and the ringdown waveform, of a regular black hole spacetime that was derived via the Loop Quantum Gravity inspired polymer quantization of spherical $4$D black holes. In contrast to most, if not all, of the other regular black holes considered in the literature, the resulting nonsingular spacetime has a single bifurcative horizon and hence no mass inflation. In the interior, the areal radius decreases to a minimum given by the Polymerization constant, $k$, and then re-expands into a Kantowski-Sachs universe. We find indications that this black hole model is stable against small scalar perturbations. We also show that an increase in the magnitude of $k$ will decrease the height of the QNM potential and gives oscillations with lower frequency and less damping.


title: Spatially resolved BPT mapping of nearby Seyfert 2 galaxies

abstract: We present spatially resolved BPT mapping of the extended narrow line regions (ENLRs) of seven nearby Seyfert 2 galaxies, using HST narrow band filter imaging. We construct the BPT diagrams using $\leq$ 0.1" resolution emission line images of [O III]$\lambda$5007, H$\alpha$, [S II]$\lambda$$\lambda$6717,6731, and H$\beta$. By mapping these diagnostic lines according to the BPT classification, we dissect the ENLR into Seyfert, LINER, and star-forming regions. The nucleus and ionization cones are dominated by Seyfert-type emission, which can be interpreted as predominantly photoionization by the active galactic nucleus (AGN). The Seyfert nucleus and ionization cones transition to and are surrounded by a LINER cocoon, extending up to $\sim$ 250 pc in thickness. The ubiquity of the LINER cocoon in Seyfert 2 galaxies suggests that the circumnuclear regions are not necessarily Seyfert-type, and LINER activity plays an important role in Seyfert 2 galaxies. We demonstrate that spatially resolved diagnostics are crucial to understanding the excitation mechanisms in different regions and the AGN-host galaxy interactions.


title: Decomposing the Internal Faraday Rotation of Black Hole Accretion Flows

abstract: Faraday rotation has been seen at millimeter wavelengths in several low luminosity active galactic nuclei, including Event Horizon Telescope (EHT) targets M87* and Sgr A*. The observed rotation measure (RM) probes the density, magnetic field, and temperature of material integrated along the line of sight. To better understand how accretion disc conditions are reflected in the RM, we perform polarized radiative transfer calculations using a set of general relativistic magneto-hydrodynamic (GRMHD) simulations appropriate for M87*. We find that in spatially resolved millimetre wavelength images on event horizon scales, the RM can vary by orders of magnitude and even flip sign. The observational consequences of this spatial structure include significant time-variability, sign-flips, and non-$\lambda^2$ evolution of the polarization plane. For some models, we find that internal rotation measure can cause significant bandwidth depolarization even across the relatively narrow fractional bandwidths observed by the EHT. We decompose the linearly polarized emission in these models based on their RM and find that emission in front of the mid-plane can exhibit orders of magnitude less Faraday rotation than emission originating from behind the mid-plane or within the photon ring. We confirm that the spatially unresolved (i.e., image integrated) RM is a poor predictor of the accretion rate, with substantial scatter stemming from time variability and inclination effects. Models can be constrained with repeated observations to characterise time variability and the degree of non-$\lambda^2$ evolution of the polarization plane.


title: Intra- and Intermolecular Inhomogeneities in the Course of Glassy Densification of Polyalcohols

abstract: To test basic assumptions inherent to most theories of molecular liquids and glasses, Infrared spectroscopy is carried out on short polyalcohols at temperatures ranging from far above to far below their glass transition. By analyzing specific vibrations, the thermal expansion of covalent bonds and hydrogen bridges is determined revealing striking differences. A comparison with density verifies the negligibility of \textit{intra}molecular expansion but exposes severe inhomogeneities on \textit{inter}molecular scale. This signals distinct roles in glassy densification and promotes an atomistic understanding.


title: NV-Diamond Magnetic Microscopy using a Double Quantum 4-Ramsey Protocol

abstract: We introduce a double quantum (DQ) 4-Ramsey measurement protocol that enables wide-field magnetic imaging using nitrogen vacancy (NV) centers in diamond, with enhanced homogeneity of the magnetic sensitivity relative to conventional single quantum (SQ) techniques. The DQ 4-Ramsey protocol employs microwave-phase alternation across four consecutive Ramsey (4-Ramsey) measurements to isolate the desired DQ magnetic signal from any residual SQ signal induced by microwave pulse errors. In a demonstration experiment employing a 1-$\mu$m-thick NV layer in a macroscopic diamond chip, the DQ 4-Ramsey protocol provides volume-normalized DC magnetic sensitivity of $\eta^\text{V}=34\,$nTHz$^{-1/2} \mu$m$^{3/2}$ across a $125\,\mu$m$ \,\times\,125\,\mu $m field of view, with about 5$\times$ less spatial variation in sensitivity across the field of view compared to a SQ measurement. The improved robustness and magnetic sensitivity homogeneity of the DQ 4-Ramsey protocol enable imaging of dynamic, broadband magnetic sources such as integrated circuits and electrically-active cells.

categories: quant-ph cond-mat.mes-hall physics.app-ph physics.ins-det


title: Evaluating the relative contribution of data sources in a Bayesian analysis with the application of estimating the size of hard to reach populations

abstract: When using multiple data sources in an analysis, it is important to understand the influence of each data source on the analysis and the consistency of the data sources with each other and the model. We suggest the use of a retrospective value of information framework in order to address such concerns. Value of information methods can be computationally difficult. We illustrate the use of computational methods that allow these methods to be applied even in relatively complicated settings. In illustrating the proposed methods, we focus on an application in estimating the size of hard to reach populations. Specifically, we consider estimating the number of injection drug users in Ukraine by combining all available data sources spanning over half a decade and numerous sub-national areas in the Ukraine. This application is of interest to public health researchers as this hard to reach population that plays a large role in the spread of HIV. We apply a Bayesian hierarchical model and evaluate the contribution of each data source in terms of absolute influence, expected influence, and level of surprise. Finally we apply value of information methods to inform suggestions on future data collection.


title: Table Scraps: An Actionable Framework for Multi-Table Data Wrangling From An Artifact Study of Computational Journalism

abstract: For the many journalists who use data and computation to report the news, data wrangling is an integral part of their work.Despite an abundance of literature on data wrangling in the context of enterprise data analysis, little is known about the specific operations, processes, and pain points journalists encounter while performing this tedious, time-consuming task. To better understand the needs of this user group, we conduct a technical observation study of 50 public repositories of data and analysis code authored by 33 professional journalists at 26 news organizations. We develop two detailed and cross-cutting taxonomies of data wrangling in computational journalism, for actions and for processes. We observe the extensive use of multiple tables, a notable gap in previous wrangling analyses. We develop a concise, actionable framework for general multi-table data wrangling that includes wrangling operations documented in our taxonomy that are without clear parallels in other work. This framework, the first to incorporate tablesas first-class objects, will support future interactive wrangling tools for both computational journalism and general-purpose use. We assess the generative and descriptive power of our framework through discussion of its relationship to our set of taxonomies.


title: Literal Encoding: Text is a first-class data encoding

abstract: Digital humanities are rooted in text analysis. However, most visualization paradigms use only categoric, ordered or quantitative data. Literal text must be considered a base data type to encode into visualizations. Literal text offers functional, perceptual, cognitive, semantic and operational benefits. These are briefly illustrated with a subset of sample visualizations focused on semantic word sequences, indicating benefits over standard graphs, maps, treemaps, bar charts and narrative layouts.


title: Quantifying Inaccuracies in Modeling COVID-19 Pandemic within a Continuous Time Picture

abstract: Typically, mathematical simulation studies on COVID-19 pandemic forecasting are based on deterministic differential equations which assume that both the number ($n$) of individuals in various epidemiological classes and the time ($t$) on which they depend are quantities that vary continuous. This picture contrasts with the discrete representation of $n$ and $t$ underlying the real epidemiological data reported in terms daily numbers of infection cases, for which a description based on finite difference equations would be more adequate. Adopting a logistic growth framework, in this paper we present a quantitative analysis of the errors introduced by the continuous time description. This analysis reveals that, although the height of the epidemiological curve maximum is essentially unaffected, the position $T_{1/2}^{c}$ obtained within the continuous time representation is systematically shifted backwards in time with respect to the position $T_{1/2}^{d}$ predicted within the discrete time representation. Rather counterintuitively, the magnitude of this temporal shift $\tau \equiv T_{1/2}^{c} - T_{1/2}^{d} < 0$ is basically insensitive to changes in infection rate $\kappa$. For a broad range of $\kappa$ values deduced from COVID-19 data at extreme situations (exponential growth in time and complete lockdown), we found a rather robust estimate $\tau \simeq -2.65\,\mbox{day}^{-1}$. Being obtained without any particular assumption, the present mathematical results apply to logistic growth in general without any limitation to a specific real system.

categories: physics.soc-ph physics.bio-ph q-bio.PE q-bio.QM


title: Verifying Results of the IBM Qiskit Quantum Circuit Compilation Flow

abstract: Realizing a conceptual quantum algorithm on an actual physical device necessitates the algorithm's quantum circuit description to undergo certain transformations in order to adhere to all constraints imposed by the hardware. In this regard, the individual high-level circuit components are first synthesized to the supported low-level gate-set of the quantum computer, before being mapped to the target's architecture---utilizing several optimizations in order to improve the compilation result. Specialized tools for this complex task exist, e.g., IBM's Qiskit, Google's Cirq, Microsoft's QDK, or Rigetti's Forest. However, to date, the circuits resulting from these tools are hardly verified, which is mainly due to the immense complexity of checking if two quantum circuits indeed realize the same functionality. In this paper, we propose an efficient scheme for quantum circuit equivalence checking---specialized for verifying results of the IBM Qiskit quantum circuit compilation flow. To this end, we combine characteristics unique to quantum computing, e.g., its inherent reversibility, and certain knowledge about the compilation flow into a dedicated equivalence checking strategy. Experimental evaluations confirm that the proposed scheme allows to verify even large circuit instances with tens of thousands of operations within seconds or even less, whereas state-of-the-art techniques frequently time-out or require substantially more runtime. A corresponding open source implementation of the proposed method is publicly available at https://github.com/iic-jku/qcec.


title: Power Grid State Estimation under General Cyber-Physical Attacks

abstract: Effective defense against cyber-physical attacks in power grid requires the capability of accurate damage assessment within the attacked area. While some solutions have been proposed to recover the phase angles and the link status (i.e., breaker status) within the attacked area, existing solutions made the limiting assumption that the grid stays connected after the attack. To fill this gap, we study the problem of recovering the phase angles and the link status under a general cyber-physical attack that may partition the grid into islands. To this end, we (i) show that the existing solutions and recovery conditions still hold if the post-attack power injections in the attacked area are known, and (ii) propose a linear programming-based algorithm that can perfectly recover the link status under certain conditions even if the post-attack power injections are unknown. Our numerical evaluations based on the Polish power grid demonstrate that the proposed algorithm is highly accurate in localizing failed links once the phase angles are known.


title: Distributed Continuous-Time Optimization with Time-Varying Objective Functions and Inequality Constraints

abstract: This paper is devoted to the distributed continuous-time optimization problem with time-varying objective functions and time-varying nonlinear inequality constraints. Different from most studied distributed optimization problems with time-invariant objective functions and constraints, the optimal solution in this paper is time varying and forms a trajectory. To minimize the global time-varying objective function subject to time-varying local constraint functions using only local information and local interaction, we present a distributed control algorithm that consists of a sliding-mode part and a Hessian-based optimization part. The asymptotical convergence of the proposed algorithm to the optimal solution is studied under suitable assumptions. The effectiveness of the proposed scheme is demonstrated through a simulation example.


title: Observation and analysis of creation, decay, and regeneration of annular soliton clusters in a lossy cubic-quintic optical medium

abstract: We observe and analyze formation, decay, and subsequent regeneration of ring-shaped clusters of (2+1)-dimensional spatial solitons (filaments) in a medium with the cubic-quintic (focusing-defocusing) self-interaction and strong dissipative nonlinearity. The cluster of filaments, that remains stable over ~17.5 Rayleigh lengths, is produced by the azimuthal modulational instability from a parent ring-shaped beam with embedded vorticity l = 1. In the course of still longer propagation, the stability of the soliton cluster is lost under the action of nonlinear losses. The annular cluster is then spontaneously regenerated due to power transfer from the reservoir provided by the unsplit part of the parent vortex ring. A (secondary) interval of the robust propagation of the regenerated cluster is identified. The experiments use a laser beam (at wavelength 800 nm), built of pulses with temporal duration 150 fs, at the repetition rate of 1 kHz, propagating in a cell filled by liquid carbon disulfide. Numerical calculations, based on a modified nonlinear Schr\"odinger equation which includes the cubic-quintic refractive terms and nonlinear losses, provide results in close agreement with the experimental findings.


title: Bilayer artificial spin ice: magnetic force switching and basic thermodynamics

abstract: We study an artificial spin ice system consisting of two identical layers separated by a height offset $h$. For small separation, the layers are shown to attract each other, provided the whole system is in the ground state. Such an attraction comes about by means of a power-law force that we compare to van der Waals forces. When magnetic monopoles occur in one (or both) layers, the scenario becomes even more interesting and these layers may also repel each other. By tuning parameters like $h$ and monopole distance, switching between attraction and repulsion may be accomplished in a feasible way. Regarding its thermodynamics, the specific heat peak shifts to lower temperature as $h$ increases.


title: Sparse Systolic Tensor Array for Efficient CNN Hardware Acceleration

abstract: Convolutional neural network (CNN) inference on mobile devices demands efficient hardware acceleration of low-precision (INT8) general matrix multiplication (GEMM). Exploiting data sparsity is a common approach to further accelerate GEMM for CNN inference, and in particular, structural sparsity has the advantages of predictable load balancing and very low index overhead. In this paper, we address a key architectural challenge with structural sparsity: how to provide support for a range of sparsity levels while maintaining high utilization of the hardware. We describe a time unrolled formulation of variable density-bound block (VDBB) sparsity that allows for a configurable number of non-zero elements per block, at constant utilization. We then describe a systolic array microarchitecture that implements this scheme, with two data reuse optimizations. Firstly, we increase reuse in both operands and partial products by increasing the number of MACs per PE. Secondly, we introduce a novel approach of moving the IM2COL transform into the hardware, which allows us to achieve a 3x data bandwidth expansion just before the operands are consumed by the datapath, reducing the SRAM power consumption. The optimizations for weight sparsity, activation sparsity and data reuse are all interrelated and therefore the optimal combination is not obvious. Therefore, we perform an design space evaluation to find the pareto-optimal design characteristics. The resulting design achieves 16.8 TOPS/W in 16nm with modest 50% model sparsity and scales with model sparsity up to 55.7TOPS/W at 87.5%. As well as successfully demonstrating the variable DBB technique, this result significantly outperforms previously reported sparse CNN accelerators.


title: Scalable multiphoton generation from cavity-synchronized single-photon sources

abstract: We propose an efficient, scalable, and deterministic scheme to generate multiple indistinguishable photons over independent channels, on demand. Our design relies on multiple single-photon sources, each coupled to a waveguide, and all of them interact with a common cavity mode. The cavity synchronizes and triggers the simultaneous emission of one photon by each source, which are collected by the waveguides. For a state-of-the-art circuit QED implementation, this scheme supports the creation of single photons with purity, indistinguishability, and efficiency of $99\%$ at rates of $\sim $MHz. We also discuss conditions to produce up to 100 photons simultaneously with generation rates of hundreds of kHz. This is orders of magnitude more efficient than previous demultiplexed sources for boson sampling and enables the realization of deterministic multi-photon sources and scalable quantum information processing with photons.


title: Don't miss the Mismatch: Investigating the Objective Function Mismatch for Unsupervised Representation Learning

abstract: Finding general evaluation metrics for unsupervised representation learning techniques is a challenging open research question, which recently has become more and more necessary due to the increasing interest in unsupervised methods. Even though these methods promise beneficial representation characteristics, most approaches currently suffer from the objective function mismatch. This mismatch states that the performance on a desired target task can decrease when the unsupervised pretext task is learned too long - especially when both tasks are ill-posed. In this work, we build upon the widely used linear evaluation protocol and define new general evaluation metrics to quantitatively capture the objective function mismatch and the more generic metrics mismatch. We discuss the usability and stability of our protocols on a variety of pretext and target tasks and study mismatches in a wide range of experiments. Thereby we disclose dependencies of the objective function mismatch across several pretext and target tasks with respect to the pretext model's representation size, target model complexity, pretext and target augmentations as well as pretext and target task types. In our experiments, we find that the objective function mismatch reduces performance by ~0.1-5.0% for Cifar10, Cifar100 and PCam in many setups, and up to ~25-59% in extreme cases for the 3dshapes dataset.