title: Measuring $H_0$ with pulsar timing arrays

abstract: Pulsar Timing Arrays have yet to convincingly observe gravitational waves. Some time ago it was pointed out by one of the authors that a dramatic enhancement of the signal would take place for particular values of the angle subtended by the source and the observed pulsar. This enhancement is due to the fact that waves propagate in a Friedmann-Lemaitre-Robertson-Walker metric where, contrary to some wide-spread belief, a simple harmonic function with a red-shifted frequency is not a solution of the equation of motion. At the first non-trivial order, proper solutions have an effective wave number that differs from the frequency. This leads to some interesting effects in Pulsar Timing Arrays whose most visible manifestation is the enhancement of the signal that, all other parameters kept fixed, is related in a simple manner to the value of $H_0$. In this work, we rederive in an alternative way the main results, extend the formalism to a more realistic setting where all components in the cosmological budget are included, investigate in detail the dependence of the signal on the various parameters involved and propose an observational set-up to hopefully detect this very relevant effect.


title: Point disclinations in the Chern-Simons geometric theory of defects

abstract: We use the Chern-Simons action for a SO(3)-connection for the description of point disclinations in the geometric theory of defects. The most general spherically symmetric SO(3)-connection with zero curvature is found. The corresponding orthogonal spherically symmetric SO(3) matrix and n-field are computed. Two examples of point disclinations are described.


title: The many Shapley values for model explanation

abstract: The Shapley value has become a popular method to attribute the prediction of a machine-learning model on an input to its base features. The use of the Shapley value is justified by citing [16] showing that it is the \emph{unique} method that satisfies certain good properties (\emph{axioms}). There are, however, a multiplicity of ways in which the Shapley value is operationalized in the attribution problem. These differ in how they reference the model, the training data, and the explanation context. These give very different results, rendering the uniqueness result meaningless. Furthermore, we find that previously proposed approaches can produce counterintuitive attributions in theory and in practice---for instance, they can assign non-zero attributions to features that are not even referenced by the model. In this paper, we use the axiomatic approach to study the differences between some of the many operationalizations of the Shapley value for attribution, and propose a technique called Baseline Shapley (BShap) that is backed by a proper uniqueness result. We also contrast BShap with Integrated Gradients, another extension of Shapley value to the continuous setting.


title: Quantum Spin Liquid in a depleted triangular lattice Iridate K$_x$Ir$_y$O$_2$

abstract: We report discovery of a new iridate family K$_x$Ir$_y$O$_2$ with depleted triangular lattice planes made up of edge sharing IrO$_6$ octahedra separated by K planes. Such a material interpolates between the triangular and honeycomb lattices and is a new playground for Kitaev physics. The materials are Mott insulators with $y = 1 - x/4$. Physical property measurements for the $x \approx 0.85$ material are reported. Using magnetic susceptibility $\chi$ versus temperature $T$ measurements we find $S_{eff} = 1/2$ moments interacting strongly with a Weiss temperature $\theta \approx - 180$~K and no magnetic order or spin freezing down to $T = 1.8$~K\@. Heat capacity shows a broad maximum around $30$~K which is insensitive to magnetic fields and a $T$-linear low temperature behaviour with $\gamma \sim 10$~mJ/mol~K$^2$. These results are consistent with a gapless QSL state in K$_{0.85}$Ir$_{0.79}$O$_2$.


title: WiFi Motion Detection: A Study into Efficacy and Classification

abstract: WiFi and security pose both an issue and act as a growing presence in everyday life. Today's motions detection implementations are severely lacking in the areas of secrecy, scope, and cost. To combat this problem, we aim to develop a motion detection system that utilizes WiFi Channel State Information (CSI), which describes how a wireless signal propagates from the transmitter to the receiver. The goal of this study is to develop a real-time motion detection and classification system that is discreet, cost-effective, and easily implementable. The system would only require an Ubuntu laptop with an Intel Ultimate N WiFi Link 5300 and a standard router. The system will be developed in two parts: (1) a robust system to track CSI variations in real-time, and (2) an algorithm to classify the motion. The system used to track CSI variance in real-time was completed in August 2018. Initial results show that introduction of motion to a previously motionless area is detected with high confidence. We present the development of (1) anomaly detection, utilizing the moving average filter implemented in the initial program and/or unsupervised machine learning, and (2) supervised machine learning algorithms to classify a set of simple motions using a proposed feature extraction methods. Lastly, classification methods such as Decision Tree, Naive Bayes, and Long Short-Term Memory can be used to classify basic actions regardless of speed, location, or orientation.


title: Using SPTpol, Planck 2015, and non-CMB data to constrain tilted spatially-flat and untilted non-flat $\Lambda$CDM, XCDM, and $\phi$CDM dark energy inflation cosmologies

abstract: We use six tilted spatially-flat and untilted non-flat dark energy cosmological models in analyses of South Pole Telescope polarization (SPTpol) cosmic microwave background (CMB) data, alone and in combination with Planck 2015 CMB data and non-CMB data, namely, the Pantheon Type Ia supernovae apparent magnitudes, a collection of baryon acoustic oscillation data points, Hubble parameter measurements, and growth rates. Although the cosmological models that best-fit the Planck CMB and non-CMB data do not provide good fits to the SPTpol data, with the $\chi^2$'s exceeding the expected value, given the uncertainties, in each model the cosmological parameter constraints from the SPTpol data and from the Planck CMB and non-CMB data are largely mutually consistent. When the smaller angular scale SPTpol data are used jointly with either the Planck data alone or with the Planck CMB and the non-CMB data to constrain untilted non-flat models, spatially-closed models remain favored over their corresponding flat limits. When used in conjunction with Planck data, non-CMB data (baryon acoustic oscillation measurements in particular, from six experiments) have significantly more constraining power than the SPTpol data.


title: Two Decades of Network Science as seen through the co-authorship network of network scientists

abstract: Complex networks have attracted a great deal of research interest in the last two decades since Watts & Strogatz, Barab\'asi & Albert and Girvan & Newman published their highly-cited seminal papers on small-world networks, on scale-free networks and on the community structure of complex networks, respectively. These fundamental papers initiated a new era of research establishing an interdisciplinary field called network science. Due to the multidisciplinary nature of the field, a diverse but not divided network science community has emerged in the past 20 years. This paper honors the contributions of network science by exploring the evolution of this community as seen through the growing co-authorship network of network scientists (here the notion refers to a scholar with at least one paper citing at least one of the three aforementioned milestone papers). After investigating various characteristics of 29,528 network science papers, we construct the co-authorship network of 52,406 network scientists and we analyze its topology and dynamics. We shed light on the collaboration patterns of the last 20 years of network science by investigating numerous structural properties of the co-authorship network and by using enhanced data visualization techniques. We also identify the most central authors, the largest communities, investigate the spatiotemporal changes, and compare the properties of the network to scientometric indicators.


title: Iterative Hard Thresholding for Low CP-rank Tensor Models

abstract: Recovery of low-rank matrices from a small number of linear measurements is now well-known to be possible under various model assumptions on the measurements. Such results demonstrate robustness and are backed with provable theoretical guarantees. However, extensions to tensor recovery have only recently began to be studied and developed, despite an abundance of practical tensor applications. Recently, a tensor variant of the Iterative Hard Thresholding method was proposed and theoretical results were obtained that guarantee exact recovery of tensors with low Tucker rank. In this paper, we utilize the same tensor version of the Restricted Isometry Property (RIP) to extend these results for tensors with low CANDECOMP/PARAFAC (CP) rank. In doing so, we leverage recent results on efficient approximations of CP decompositions that remove the need for challenging assumptions in prior works. We complement our theoretical findings with empirical results that showcase the potential of the approach.


title: On the Close Correspondence between Storm-time ULF Wave Power and the POES VLF Chorus Wave Amplitude Proxy

abstract: Ground-based Pc5 ULF wave power in multiple ground-based meridians is compared to the VLF wave amplitude proxy, derived from POES precipitation, for the 33 storms studied by Li et al. [2015]. The results reveal common L-shell and time profiles for the ULF waves and VLF proxy for every single storm, especially at $L\leq 6$, and identical discrimination between efficient and inefficient radiation belt electron acceleration. The observations imply either ULF waves play a role in driving precipitation which is falsely interpreted as VLF wave power in the proxy, ULF waves drive VLF waves (the reverse being energetically unfeasible), or both have a common driver with nearly identical L-shell and time-dependence. Global ground-based ULF wave power coherence implies a small number of meridians can be used to estimate storm-time radial diffusion coefficients. However, the strong correspondence between ULF wave power and VLF wave proxy complicates causative assessments of electron acceleration.


title: Rayleigh Random Flights on the Poisson line SIRSN

abstract: We study scale-invariant Rayleigh Random Flights ("RRF") in random environments given by planar Scale-Invariant Random Spatial Networks ("SIRSN") based on speed-marked Poisson line processes. A natural one-parameter family of such RRF (with scale-invariant dynamics) can be viewed as producing "randomly-broken local geodesics" on the SIRSN; we aim to shed some light on a conjecture that a (non-broken) geodesic on such a SIRSN will never come to a complete stop en route. (If true, then all such geodesics can be represented as doubly-infinite sequences of sequentially connected line segments. This would justify a natural procedure for computing geodesics.) The family of these RRF ("SIRSNRRF"), is introduced via a novel axiomatic theory of abstract scattering representations for Markov chains (itself of independent interest). Palm conditioning (specifically the Mecke-Slivnyak theorem for Palm probabilities of Poisson point processes) and ideas from the ergodic theory of random walks in random environments are used to show that at a critical value of the parameter the speed of the scale-invariant SIRSNRRF neither diverges to infinity nor tends to zero, thus supporting the conjecture.


title: A mathematical model describing the localization and spread of influenza A virus infection within the human respiratory tract

abstract: Within the human respiratory tract (HRT), viruses diffuse through the periciliary fluid (PCF) bathing the epithelium, and travel upwards via advection towards the nose and mouth, as the mucus escalator entrains the PCF. While many mathematical models (MMs) to date have described the course of influenza A virus (IAV) infections in vivo, none have considered the impact of both diffusion and advection on the kinetics and localization of the infection. The MM herein represents the HRT as a one-dimensional track extending from the nose down to a depth of 30 cm, wherein stationary cells interact with the concentration of IAV which move along within the PCF. When IAV advection and diffusion are both considered, the former is found to dominate infection kinetics, and a 10-fold increase in the virus production rate is required to counter its effects. The MM predicts that advection prevents infection from disseminating below the depth at which virus first deposits. Because virus is entrained upwards, the upper HRT sees the most virus, whereas the lower HRT sees far less. As such, infection peaks and resolves faster in the upper than in the lower HRT, making it appear as though infection progresses from the upper towards the lower HRT. When the spatial MM is expanded to include cellular regeneration and an immune response, it can capture the time course of infection with a seasonal and an avian IAV strain by shifting parameters in a manner consistent with what is expected to differ between these two types of infection. The impact of antiviral therapy with neuraminidase inhibitors was also investigated. This new MM offers a convenient and unique platform from which to study the localization and spread of respiratory viral infections within the HRT.


title: Improved bounds for the sunflower lemma

abstract: A sunflower with $r$ petals is a collection of $r$ sets so that the intersection of each pair is equal to the intersection of all of them. Erd\H{o}s and Rado proved the sunflower lemma: for any fixed $r$, any family of sets of size $w$, with at least about $w^w$ sets, must contain a sunflower with $r$ petals. The famous sunflower conjecture states that the bound on the number of sets can be improved to $c^w$ for some constant $c$. In this paper, we improve the bound to about $(\log w)^w$. In fact, we prove the result for a robust notion of sunflowers, for which the bound we obtain is sharp up to lower order terms.


title: Minimum Description Length Revisited

abstract: This is an up-to-date introduction to and overview of the Minimum Description Length (MDL) Principle, a theory of inductive inference that can be applied to general problems in statistics, machine learning and pattern recognition. While MDL was originally based on data compression ideas, this introduction can be read without any knowledge thereof. It takes into account all major developments since 2007, the last time an extensive overview was written. These include new methods for model selection and averaging and hypothesis testing, as well as the first completely general definition of {\em MDL estimators}. Incorporating these developments, MDL can be seen as a powerful extension of both penalized likelihood and Bayesian approaches, in which penalization functions and prior distributions are replaced by more general luckiness functions, average-case methodology is replaced by a more robust worst-case approach, and in which methods classically viewed as highly distinct, such as AIC vs BIC and cross-validation vs Bayes can, to a large extent, be viewed from a unified perspective.


title: Simulation Model of Two-Robot Cooperation in Common Operating Environment

abstract: The article considers a simulation modelling problem related to the chess game process occurring between two three-tier manipulators. The objective of the game construction lies in developing the procedure of effective control of the autonomous manipulator robots located in a common operating environment. The simulation model is a preliminary stage of building a natural complex that would provide cooperation of several manipulator robots within a common operating environment. The article addresses issues of training and research.


title: Dialogue Coherence Assessment Without Explicit Dialogue Act Labels

abstract: Recent dialogue coherence models use the coherence features designed for monologue texts, e.g. nominal entities, to represent utterances and then explicitly augment them with dialogue-relevant features, e.g., dialogue act labels. It indicates two drawbacks, (a) semantics of utterances is limited to entity mentions, and (b) the performance of coherence models strongly relies on the quality of the input dialogue act labels. We address these issues by introducing a novel approach to dialogue coherence assessment. We use dialogue act prediction as an auxiliary task in a multi-task learning scenario to obtain informative utterance representations for coherence assessment. Our approach alleviates the need for explicit dialogue act labels during evaluation. The results of our experiments show that our model substantially (more than 20 accuracy points) outperforms its strong competitors on the DailyDialogue corpus, and performs on par with them on the SwitchBoard corpus for ranking dialogues concerning their coherence.


title: Lower bounds and fixed points for the centered Hardy--Littlewood maximal operator

abstract: For all $p>1$ and all centrally symmetric convex bodies $K\subset \mathbb{R}^d$ define $Mf$ as the centered maximal function associated to $K$. We show that when $d=1$ or $d=2$, we have $||Mf||_p\ge (1+\epsilon(p,K))||f||_p$. For $d\ge 3$, let $q_0(K)$ be the infimum value of $p$ for which $M$ has a fixed point. We show that for generic shapes $K$, we have $q_0(K)>q_0(B(0,1))$.


title: On the dependent product in toposes

abstract: We give an explicit construction of the dependent product in an elementary topos, and a site-theoretic description for it in the case of a Grothendieck topos.


title: Time series model selection with a meta-learning approach; evidence from a pool of forecasting algorithms

abstract: One of the challenging questions in time series forecasting is how to find the best algorithm. In recent years, a recommender system scheme has been developed for time series analysis using a meta-learning approach. This system selects the best forecasting method with consideration of the time series characteristics. In this paper, we propose a novel approach to focusing on some of the unanswered questions resulting from the use of meta-learning in time series forecasting. Therefore, three main gaps in previous works are addressed including, analyzing various subsets of top forecasters as inputs for meta-learners; evaluating the effect of forecasting error measures; and assessing the role of the dimensionality of the feature space on the forecasting errors of meta-learners. All of these objectives are achieved with the help of a diverse state-of-the-art pool of forecasters and meta-learners. For this purpose, first, a pool of forecasting algorithms is implemented on the NN5 competition dataset and ranked based on the two error measures. Then, six machine-learning classifiers known as meta-learners, are trained on the extracted features of the time series in order to assign the most suitable forecasting method for the various subsets of the pool of forecasters. Furthermore, two-dimensionality reduction methods are implemented in order to investigate the role of feature space dimension on the performance of meta-learners. In general, it was found that meta-learners were able to defeat all of the individual benchmark forecasters; this performance was improved even after applying the feature selection method.


title: Constraining the Sea Quark Distributions Through W$^\pm$ Cross Section Ratio Measurements at STAR

abstract: Over the past several years, parton distribution functions (PDFs) have become more precise. However there are still kinematic regions where more data are needed to help constrain global PDF extractions, such as the ratio of the sea quark distributions $\bar{d}$/$\bar{u}$ near the valence region. Furthermore, current measurements appear to suggest different high-$x$ behaviors of this ratio. The $W$ cross section ratio ($W^+$/$W^-$) is sensitive to the unpolarized quark distributions at large $Q^2$ set by the $W$ mass. Such a measurement can be used to help constrain the $\bar{d}$/$\bar{u}$ ratio. The STAR experiment at RHIC is well equipped to measure the leptonic decays of $W$ bosons, in the mid-pseudorapdity range $\left(|\eta| \leq 1 \right)$, produced in proton-proton collisions at $\sqrt{s}$ = 500/510 GeV. At these kinematics STAR is sensitive to quark distributions near $x$ of 0.16. STAR can also measure $W^+$/$W^-$ in a more forward region ranging from 1.0 $< \eta <$1.5, which extends the sea quark sensitivity to higher $x$. RHIC runs from 2011 through 2013 have collected about 350 pb$^{-1}$ of integrated luminosity, and an additional 350 pb$^{-1}$ from the 2017 run. These proceedings will present preliminary results of the 2011-2013 $W^+$/$W^-$ cross section ratio measurements. Additionally, the $W/Z$ cross section ratio, differential and total $W$ and $Z$ cross sections are presented.


title: On spectral curves and complexified boundaries of the phase-lock areas in a model of Josephson junction

abstract: The paper deals with a three-parameter family of special double confluent Heun equations that was introduced and studied by V.M.Buchstaber and S.I.Tertychnyi as an equivalent presentation of a model of overdamped Josephson junction in superconductivity. The parameters are $l,\lambda,\mu\in\mathbb R$. Buchstaber and Tertychnyi described those parameter values, for which the corresponding equation has a polynomial solution. They have shown that for $\mu\neq0$ this happens exactly when $l\in\mathbb N$ and the parameters $(\lambda,\mu)$ lie on an algebraic curve $\Gamma_l\subset\mathbb C^2_{(\lambda,\mu)}$ called the $l$-th spectral curve and defined as zero locus of determinant of a remarkable three-diagonal $l\times l$-matrix. They studied the real spectral curves and obtained important results with applications to phase-lock areas in model of Josephson junction, which is a family of dynamical systems on 2-torus. In the present paper we prove irreducibility of complex spectral curves. We also calculate their genera for $l\leqslant20$ and present a conjecture on general genus formula. We apply the irreducibility result to the phase-lock areas, which are those level sets of the rotation number function $\rho$ on the parameter space of the above-mentioned family of dynamical systems that have non-empty interiors. The family of their boundaries is a countable union of analytic surfaces. We show that, unexpectedly, its complexification is a complex analytic subset consisting of just four irreducible components, and we describe them. We present a Monotonicity Conjecture on the evolution of the phase-lock area portraits and a partial positive result towards its confirmation.


title: Some remarks on the asymptotic profile of solutions to structurally damped $\sigma$-evolution equations

abstract: In this paper, we are interested in analyzing the asymptotic profiles of solutions to the Cauchy problem for linear structurally damped $\sigma$-evolution equations in $L^2$-sense. Depending on the parameters $\sigma$ and $\delta$ we would like to not only indicate approximation formula of solutions but also recognize the optimality of their decay rates as well in the distinct cases of parabolic like damping and $\sigma$-evolution like damping. Moreover, such results are also discussed when we mix these two kinds of damping terms in a $\sigma$-evolution equation to investigate how each of them affects the asymptotic profile of solutions.


title: A Hybrid Method for Online Trajectory Planning of Mobile Robots in Cluttered Environments

abstract: This paper presents a method for online trajectory planning in known environments. The proposed algorithm is a fusion of sampling-based techniques and model-based optimization via quadratic programming. The former is used to efficiently generate an obstacle-free path while the latter takes into account the robot dynamical constraints to generate a time-dependent trajectory. The main contribution of this work lies on the formulation of a convex optimization problem over the generated obstacle-free path that is guaranteed to be feasible. Thus, in contrast with previously proposed methods, iterative formulations are not required. The proposed method has been compared with state-of-the-art approaches showing a significant improvement in success rate and computation time. To illustrate the effectiveness of this approach for online planning, the proposed method was applied to the fluid autonomous navigation of a quadcopter in multiple environments consisting of up to two hundred obstacles. The scenarios hereinafter presented are some of the most densely cluttered experiments for online planning and navigation reported to date. See video at https://youtu.be/DJ1IZRL5t1Q


title: The double traveling salesman problem with partial last-in-first-out loading constraints

abstract: In this paper, we introduce the Double Traveling Salesman Problem with Partial Last-In-First-Out Loading Constraints (DTSPPL). It is a pickup-and-delivery single-vehicle routing problem, where all pickup operations must be performed before any delivery one because the pickup and delivery areas are geographically separated. The vehicle collects items in the pickup area and loads them into its container, a horizontal stack. After performing all pickup operations, the vehicle begins delivering the items in the delivery area. Loading and unloading operations must obey a partial Last-In-First-Out (LIFO) policy, i.e., a version of the LIFO policy that may be violated within a given reloading depth. The objective of the DTSPPL is to minimize the total cost, which involves the total distance traveled by the vehicle and the number of items that are unloaded and then reloaded due to violations of the standard LIFO policy. We formally describe the DTSPPL through two Integer Linear Programming (ILP) formulations and propose a heuristic algorithm based on the Biased Random-Key Genetic Algorithm (BRKGA) to find high-quality solutions. The performance of the proposed solution approaches is assessed over a broad set of instances. Computational results have shown that both ILP formulations have been able to solve only the smaller instances, whereas the BRKGA obtained good quality solutions for almost all instances, requiring short computational times.


title: Applying machine learning optimization methods to the production of a quantum gas

abstract: We apply three machine learning strategies to optimize the atomic cooling processes utilized in the production of a Bose-Einstein condensate (BEC). For the first time, we optimize both laser cooling and evaporative cooling mechanisms simultaneously. We present the results of an evolutionary optimization method (Differential Evolution), a method based on non-parametric inference (Gaussian Process regression) and a gradient-based function approximator (Artificial Neural Network). Online optimization is performed using no prior knowledge of the apparatus, and the learner succeeds in creating a BEC from completely randomized initial parameters. Optimizing these cooling processes results in a factor of four increase in BEC atom number compared to our manually-optimized parameters. This automated approach can maintain close-to-optimal performance in long-term operation. Furthermore, we show that machine learning techniques can be used to identify the main sources of instability within the apparatus.


title: Cavity electromechanics with parametric mechanical driving

abstract: Microwave optomechanical circuits have been demonstrated in the past years to be extremely powerfool tools for both, exploring fundamental physics of macroscopic mechanical oscillators as well as being promising candidates for novel on-chip quantum limited microwave devices. In most experiments so far, the mechanical oscillator is either used as a passive device element and its displacement is detected using the superconducting cavity or manipulated by intracavity fields. Here, we explore the possibility to directly and parametrically manipulate the mechanical nanobeam resonator of a cavity electromechanical system, which provides additional functionality to the toolbox of microwave optomechanical devices. In addition to using the cavity as an interferometer to detect parametrically modulated mechanical displacement and squeezed thermomechanical motion, we demonstrate that parametric modulation of the nanobeam resonance frequency can realize a phase-sensitive parametric amplifier for intracavity microwave photons. In contrast to many other microwave amplification schemes using electromechanical circuits, the presented technique allows for simultaneous cooling of the mechanical element, which potentially enables this type of optomechanical microwave amplifier to be quantum-limited.


title: DynGraph2Seq: Dynamic-Graph-to-Sequence Interpretable Learning for Health Stage Prediction in Online Health Forums

abstract: Online health communities such as the online breast cancer forum enable patients (i.e., users) to interact and help each other within various subforums, which are subsections of the main forum devoted to specific health topics. The changing nature of the users' activities in different subforums can be strong indicators of their health status changes. This additional information could allow health-care organizations to respond promptly and provide additional help for the patient. However, modeling complex transitions of an individual user's activities among different subforums over time and learning how these correspond to his/her health stage are extremely challenging. In this paper, we first formulate the transition of user activities as a dynamic graph with multi-attributed nodes, then formalize the health stage inference task as a dynamic graph-to-sequence learning problem, and hence propose a novel dynamic graph-to-sequence neural networks architecture (DynGraph2Seq) to address all the challenges. Our proposed DynGraph2Seq model consists of a novel dynamic graph encoder and an interpretable sequence decoder that learn the mapping between a sequence of time-evolving user activity graphs and a sequence of target health stages. We go on to propose dynamic graph hierarchical attention mechanisms to facilitate the necessary multi-level interpretability. A comprehensive experimental analysis of its use for a health stage prediction task demonstrates both the effectiveness and the interpretability of the proposed models.


title: EPIC-Fusion: Audio-Visual Temporal Binding for Egocentric Action Recognition

abstract: We focus on multi-modal fusion for egocentric action recognition, and propose a novel architecture for multi-modal temporal-binding, i.e. the combination of modalities within a range of temporal offsets. We train the architecture with three modalities -- RGB, Flow and Audio -- and combine them with mid-level fusion alongside sparse temporal sampling of fused representations. In contrast with previous works, modalities are fused before temporal aggregation, with shared modality and fusion weights over time. Our proposed architecture is trained end-to-end, outperforming individual modalities as well as late-fusion of modalities. We demonstrate the importance of audio in egocentric vision, on per-class basis, for identifying actions as well as interacting objects. Our method achieves state of the art results on both the seen and unseen test sets of the largest egocentric dataset: EPIC-Kitchens, on all metrics using the public leaderboard.


title: General Composite Non-Abelian Strings and Flag Manifold Sigma Models

abstract: We fully investigate the symmetry breaking patterns occurring upon creation of composite non-Abelian strings: vortex strings in non-Abelian theories where different sets of colours have different amounts of flux. After spontaneous symmetry breaking, there remains some internal colour degrees of freedom attached to these objects, which we argue must exist in a Flag manifold, a more general kind of projective space than both $\mathbb{CP}(N)$ and the Grassmannian manifold. These strings are expected to be BPS, since its constituents are. We demonstrate that this is true and construct a low-energy effective action for the fluctuations of the internal Flag moduli, which we then re-write it in two different ways for the dynamics of these degrees of freedom: a gauged linear sigma model with auxiliary fields and a non-linear sigma model with an explicit target space metric for the Flag Manifolds, both of which $\mathcal{N}=(2,2)$ supersymmetric. We finish by performing some groundwork analysis of the resulting theory.


title: Search for supersymmetry using Higgs boson to diphoton decays at $\sqrt{s} =$ 13 TeV

abstract: A search for supersymmetry (SUSY) is presented where at least one Higgs boson is produced and decays to two photons in the decay chains of pair-produced SUSY particles. Two analysis strategies are pursued: one focused on strong SUSY production and the other focused on electroweak SUSY production. The presence of charged leptons, additional Higgs boson candidates, and various kinematic variables are used to categorize events into search regions that are sensitive to different SUSY scenarios. The results are based on data from proton-proton collisions at the Large Hadron Collider at a center-of-mass energy of 13 TeV collected by the CMS experiment, corresponding to an integrated luminosity of 77.5 fb$^{-1}$. No statistically significant excess of events is observed relative to the standard model expectations. We exclude bottom squark pair production for bottom squark masses below 530 GeV and a lightest SUSY particle mass of 1 GeV; wino-like chargino-neutralino production in gauge-mediated SUSY breaking (GMSB) for chargino and neutralino masses below 235 GeV with a gravitino mass of 1 GeV; and higgsino-like chargino-neutralino production in GMSB, where the neutralino decays exclusively to a Higgs boson and a gravitino for neutralino masses below 290 GeV.


title: IMPRESSION -- Prediction of NMR Parameters for 3-dimensional chemical structures using Machine Learning with near quantum chemical accuracy

abstract: The IMPRESSION (Intelligent Machine PREdiction of Shift and Scalar Information Of Nuclei) machine learning system provides an efficient and accurate route to the prediction of NMR parameters from 3-dimensional chemical structures. Here we demonstrate that machine learning predictions, trained on quantum chemical computed values for NMR parameters, are essentially as accurate but computationally much more efficient (tens of milliseconds per molecule) than quantum chemical calculations (hours/days per molecule). Training the machine learning systems on quantum chemical, rather than experimental, data circumvents the need for existence of large, structurally diverse, error-free experimental databases and makes IMPRESSION applicable to solving 3-dimensional problems such as molecular conformation and isomerism


title: A Pieri rule for Demazure characters of the general linear group

abstract: The Pieri rule is a nonnegative, multiplicity-free formula for the Schur function expansion of the product of an arbitrary Schur function with a single row Schur function. Key polynomials are characters of Demazure modules for the general linear group that generalize the Schur function basis of symmetric functions to a basis of the full polynomial ring. We prove a nonsymmetric generalization of the Pieri rule by giving a cancellation-free, multiplicity-free formula for the key polynomial expansion of the product of an arbitrary key polynomial with a single part key polynomial. Our proof is combinatorial, generalizing the Robinson--Schensted--Knuth insertion algorithm on tableaux to an insertion algorithm on Kohnert diagrams.


title: Blockchain based access control systems: State of the art and challenges

abstract: Access to the system resources. The current access control systems face many problems, such as the presence of the third-party, inefficiency, and lack of privacy. These problems can be addressed by blockchain, the technology that received major attention in recent years and has many potentials. In this study, we overview the problems of the current access control systems, and then, we explain how blockchain can help to solve them. We also present an overview of access control studies and proposed platforms in different domains. This paper presents the state of the art and the challenges of blockchain-based access control systems.


title: In-phase and anti-phase synchronization in a laser frequency comb

abstract: Coupled clocks are a classic example of a synchronization system leading to periodic collective oscillations. This phenomenon already attracted the attention of Christian Huygens back in 1665,who described it as a kind of "sympathy" among oscillators. In this work we describe the formation of two types of laser frequency combs as a system of oscillators coupled through the beating of the lasing modes. We experimentally show two completely different types of synchronizations in a quantum dot laser { in-phase and splay states. Both states can be generated in the same device, just by varying the damping losses of the system. This effectively modifes the coupling among the oscillators. The temporal output of the laser is characterized using both linear and quadratic autocorrelation techniques. Our results show that both pulses and frequency-modulated states can be generated on demand. These findings allow to connect laser frequency combs produced by amplitude-modulated and frequency-modulated lasers, and link these to pattern formation in coupled systems such as Josephson-junction arrays.


title: ColorNet -- Estimating Colorfulness in Natural Images

abstract: Measuring the colorfulness of a natural or virtual scene is critical for many applications in image processing field ranging from capturing to display. In this paper, we propose the first deep learning-based colorfulness estimation metric. For this purpose, we develop a color rating model which simultaneously learns to extracts the pertinent characteristic color features and the mapping from feature space to the ideal colorfulness scores for a variety of natural colored images. Additionally, we propose to overcome the lack of adequate annotated dataset problem by combining/aligning two publicly available colorfulness databases using the results of a new subjective test which employs a common subset of both databases. Using the obtained subjectively annotated dataset with 180 colored images, we finally demonstrate the efficacy of our proposed model over the traditional methods, both quantitatively and qualitatively.


title: Predicting Animation Skeletons for 3D Articulated Models via Volumetric Nets

abstract: We present a learning method for predicting animation skeletons for input 3D models of articulated characters. In contrast to previous approaches that fit pre-defined skeleton templates or predict fixed sets of joints, our method produces an animation skeleton tailored for the structure and geometry of the input 3D model. Our architecture is based on a stack of hourglass modules trained on a large dataset of 3D rigged characters mined from the web. It operates on the volumetric representation of the input 3D shapes augmented with geometric shape features that provide additional cues for joint and bone locations. Our method also enables intuitive user control of the level-of-detail for the output skeleton. Our evaluation demonstrates that our approach predicts animation skeletons that are much more similar to the ones created by humans compared to several alternatives and baselines.


title: Transfer Learning for Relation Extraction via Relation-Gated Adversarial Learning

abstract: Relation extraction aims to extract relational facts from sentences. Previous models mainly rely on manually labeled datasets, seed instances or human-crafted patterns, and distant supervision. However, the human annotation is expensive, while human-crafted patterns suffer from semantic drift and distant supervision samples are usually noisy. Domain adaptation methods enable leveraging labeled data from a different but related domain. However, different domains usually have various textual relation descriptions and different label space (the source label space is usually a superset of the target label space). To solve these problems, we propose a novel model of relation-gated adversarial learning for relation extraction, which extends the adversarial based domain adaptation. Experimental results have shown that the proposed approach outperforms previous domain adaptation methods regarding partial domain adaptation and can improve the accuracy of distance supervised relation extraction through fine-tuning.


title: Scaling of wetting and pre-wetting transitions on nano-patterned walls

abstract: We consider a nano-patterned planar wall consisting of a periodic array of stripes of width $L$, which are completely wet by liquid (contact angle $\theta=0$), separated by regions of width $D$ which are completely dry (contact angle $\theta=\pi)$. Using microscopic Density Functional Theory we show that in the presence of long-ranged dispersion forces, the wall-gas interface undergoes a first-order wetting transition, at bulk coexistence, as the separation $D$ is reduced to a value $D_w\propto\ln L$, induced by the bridging between neighboring liquid droplets. Associated with this is a line of pre-wetting transitions occurring off coexistence. By varying the stripe width $L$ we show that the pre-wetting line shows universal scaling behaviour and data collapse. This verifies predictions based on mesoscopic models for the scaling properties associated with finite-size effects at complete wetting including the logarithmic singular contribution to the surface free-energy.


title: Navigation of a Quadratic Potential with Ellipsoidal Obstacles

abstract: Given a convex quadratic potential of which its minimum is the agent's goal and a Euclidean space populated with ellipsoidal obstacles, one can construct a Rimon-Koditschek (RK) artificial potential to navigate. Its negative gradient attracts the agent toward the goal and repels the agent away from the boundary of the obstacles. This is a popular approach to navigation problems since it can be implemented with local spatial information that is acquired during operation time. However, navigation is only successful in situations where the obstacles are not too eccentric (flat). This paper proposes a modification to gradient dynamics that allows successful navigation of an environment with a quadratic cost and ellipsoidal obstacles regardless of their eccentricity. This is accomplished by altering gradient dynamics with a Hessian correction that is intended to imitate worlds with spherical obstacles in which RK potentials are known to work. The resulting dynamics simplify by the quadratic form of the obstacles. Convergence to the goal and obstacle avoidance is established from almost every initial position (up to a set of measure one) in the free space, with mild conditions on the location of the target. Results are corroborated empirically with numerical simulations.


title: Speeding-up Ab Initio Molecular Dynamics with Hybrid Functionals using Adaptively Compressed Exchange Operator based Multiple Timestepping

abstract: Ab initio molecular dynamics (AIMD) simulations using hybrid density functionals and plane waves are of great interest owing to the accuracy of this approach in treating condensed matter systems. On the other hand, such AIMD calculations are not routinely carried out since the computational cost involved in applying the Hartree Fock exchange operator is very high. In this work, we make use of a strategy that combines adaptively compressed exchange operator formulation and multiple time step integration to significantly reduce the computational cost of these simulations. We demonstrate the efficiency of this approach for a realistic condensed matter system.


title: Spatially resolved spectroscopy of the debris disk HD 32297: Further evidence of small dust grains

abstract: Spectro-photometry of debris disks in total intensity and polarimetry can provide new insight into the properties of the dust grains therein (size distribution and optical properties). We aim to constrain the morphology of the highly inclined debris disk HD 32297. We also intend to obtain spectroscopic and polarimetric measurements to retrieve information on the particle size distribution within the disk for certain grain compositions. We observed HD 32297 with SPHERE in Y, J, and H bands in total intensity and in J band in polarimetry. The observations are compared to synthetic models of debris disks and we developed methods to extract the photometry in total intensity overcoming the data-reduction artifacts, namely the self-subtraction. The spectro-photometric measurements averaged along the disk mid-plane are then compared to model spectra of various grain compositions. These new images reveal the very inner part of the system as close as 0.15". The disk image is mostly dominated by the forward scattering making one side (half-ellipse) of the disk more visible, but observations in total intensity are deep enough to also detect the back side for the very first time. The images as well as the surface brightness profiles of the disk rule out the presence of a gap as previously proposed. We do not detect any significant asymmetry between the northeast and southwest sides of the disk. The spectral reflectance features a "gray to blue" color which is interpreted as the presence of grains far below the blowout size. The presence of sub-micron grains in the disk is suspected to be the result of gas drag and/or "avalanche mechanisms". The blue color of the disk could be further investigated with additional total intensity and polarimetric observations in K and H bands respectively to confirm the spectral slope and the fraction of polarization.


title: Dynamics of Bubbling Wave Maps with Prescribed Radiation

abstract: We study energy critical one-equivariant wave maps taking values in the two-sphere. It is known that any finite energy wave map that develops a singularity does so by concentrating the energy of (possibly) several copies of the ground state harmonic map at the origin. If only a single bubble of energy is concentrated, the solution decomposes into a dynamically rescaled harmonic map plus a term that accounts for the energy that radiates away from the singularity. In this paper, we construct blow up solutions by prescribing the radiative component of the map. In addition, we give a sharp classification of the dynamical blow up rate for every solution with this prescribed radiation.


title: A Decomposition and Metric-Based Evaluation Framework for Microservices

abstract: Migrating from monolithic systems into microservice is a very complex task. Companies are commonly decomposing the monolithic system manually, analyzing dependencies of the monolith and then assessing different decomposition options. The goal of our work is two-folded: 1) we provide a microservice measurement framework to objectively evaluate and compare the quality of microservices-based systems; 2) we propose a decomposition system based on business process mining. The microservice measurement framework can be applied independently from the decomposition process adopted, but is also useful to continuously evaluate the architectural evolution of a system. Results show that the decomposition framework helps companies to easily identify the different decomposition options. The measurement framework can help to decrease the subjectivity of the decision between different decomposition options and to evaluate architectural erosion in existing systems.


title: Reflections on Virasoro circuit complexity and Berry phase

abstract: Recently, the notion of circuit complexity defined in symmetry group manifolds has been related to geometric actions which generally arise in the coadjoint orbit method in representation theory and play an important role in geometric quantization. On the other hand, it is known that there exists a precise relation between geometric actions and Berry phases defined in group representations. Motivated by these connections, we elaborate on a relation between circuit complexity and the group theoretic Berry phase. As the simplest setup relevant for holography, we discuss the case of two dimensional conformal field theories. In the large central charge limit, we identify the computational cost function with the Berry connection in the unitary representation of the Virasoro group. We then use the latter identification to express the Berry phase in terms of the Virasoro circuit complexity. The former can be seen as the holonomy of the Berry connection along the path in the group manifold which defines the protocol. In addition, we derive a proportionality relation between Virasoro circuit complexity and the logarithm of the inner product between a particularly chosen reference state and the prepared target state. In this sense, the logarithmic formula turns out to be approximating the complexity up to some additive constant if the building blocks of the circuit are taken to be the underlying symmetry gates. Predictions based on this formula have recently been shown to coincide with the holographic complexity proposals and the path integral optimization procedure. The found connections may therefore help to better understand such coincidences. We also discuss that our findings, put together with earlier observations, may suggest a connection between the Virasoro Berry phase and the complexity measure in the path integral optimization proposal.


title: Entropy-cooled nonequilibrium states of the Hubbard model

abstract: We show that the recently proposed cooling-by-doping mechanism allows to efficiently prepare interesting nonequilibrium states of the Hubbard model. Using nonequilibrium dynamical mean field theory and a particle-hole symmetric setup with dipolar excitations to full and empty bands we produce cold photo-doped Mott insulating states with a sharp Drude peak in the optical conductivity, a superconducting state in the repulsive Hubbard model with an inverted population, and $\eta$-paired states in systems with a large density of doublons and holons. The reshuffling of entropy into full and empty bands not only provides an efficient cooling mechanism, it also allows to overcome thermalization bottlenecks and slow dynamics that have been observed in systems cooled by the coupling to boson baths.


title: Experimental Demonstration of Exceptional Points of Degeneracy in Linear Time Periodic Systems and Exceptional Sensitivity

abstract: We present the experimental demonstration of the occurrence of exceptional points of degeneracy (EPDs) in a single resonator by introducing a linear time-periodic variation of one of its components, in contrast to EPDs in parity time (PT)-symmetric systems that require two coupled resonators with precise values of gain and loss. In the proposed scheme, only the tuning of the modulation frequency is required that is easily achieved in electronic systems. The EPD is a point in a system parameters' space at which two or more eigenstates coalesce, and this leads to unique properties not occurring at other non-degenerate operating points. We show theoretically and experimentally the existence of a second order EPD in a time-varying single resonator. Furthermore, we measure the sensitivity of the proposed system to a small structural perturbation and show that the operation of the system at an EPD dramatically boosts its sensitivity performance to very small perturbations. Also, we show experimentally how this unique sensitivity induced by an EPD can be used to devise new exceptionally-sensitive sensors based on a single resonator by simply applying time modulation.


title: Six-dimensional gauge theories and (twisted) generalized cohomology

abstract: We consider the global aspects of the 6-dimensional $\mathcal{N}=(1, 0)$ theory arising from the coupling of the vector multiplet to the tensor multiplet. We show that the Yang-Mills field and its dual, when both are abelianized, combine to define a class in twisted cohomology with the twist arising from the class of the $B$-field, in a duality-symmetric manner. We then show that this lifts naturally to a class in twisted (differential) K-theory. Alternatively, viewing the B-field in both $\mathcal{N}=(1,0)$ and $\mathcal{N}=(2,0)$ theories, not as a twist but as an invertible element, leads to a description within untwisted chromatic level two generalized cohomology theories, including forms of elliptic cohomology and Morava K-theory.


title: Aspects of AdS$_2$ classification in M-theory: Solutions with mesonic and baryonic charges

abstract: We construct necessary and sufficient geometric conditions for a class of AdS$_2$ solutions of M-theory with, at least, minimal supersymmetry to exist. We generalize previous results in the literature for ${\cal N}=(2,0)$ supersymmetry in AdS$_2$ to ${\cal N}=(1,0)$. When the solution can be locally described as AdS$_2\times \Sigma_g \times \,$SE$_7$ with $\Sigma_g$ a Riemann surface of genus $g$ and SE$_7$ a seven-dimensional Sasaki-Einstein manifold, we clarify and unify various solutions present in the literature. In the case of SE$_7=Q^{1,1,1}$ we find a new solution with baryonic and mesonic charges turned on simultaneously.


title: Unfolding as Quantum Annealing

abstract: High-energy physics is replete with hard computational problems and it is one of the areas where quantum computing could be used to speed up calculations. We present an implementation of likelihood-based regularized unfolding on a quantum computer. The inverse problem is recast in terms of quadratic unconstrained binary optimization (QUBO), which has the same form of the Ising hamiltonian and hence it is solvable on a programmable quantum annealer. We tested the method using a model that captures the essence of the problem, and compared the results with a baseline method commonly used in precision measurements at the Large Hadron Collider (LHC) at CERN. The unfolded distribution is in very good agreement with the original one. We also show how the method can be extended to include the effect of nuisance parameters representing sources of systematic uncertainties affecting the measurement.


title: Adversarial-Based Knowledge Distillation for Multi-Model Ensemble and Noisy Data Refinement

abstract: Generic Image recognition is a fundamental and fairly important visual problem in computer vision. One of the major challenges of this task lies in the fact that single image usually has multiple objects inside while the labels are still one-hot, another one is noisy and sometimes missing labels when annotated by humans. In this paper, we focus on tackling these challenges accompanying with two different image recognition problems: multi-model ensemble and noisy data recognition with a unified framework. As is well-known, usually the best performing deep neural models are ensembles of multiple base-level networks, as it can mitigate the variation or noise containing in the dataset. Unfortunately, the space required to store these many networks, and the time required to execute them at runtime, prohibit their use in applications where test sets are large (e.g., ImageNet). In this paper, we present a method for compressing large, complex trained ensembles into a single network, where the knowledge from a variety of trained deep neural networks (DNNs) is distilled and transferred to a single DNN. In order to distill diverse knowledge from different trained (teacher) models, we propose to use adversarial-based learning strategy where we define a block-wise training loss to guide and optimize the predefined student network to recover the knowledge in teacher models, and to promote the discriminator network to distinguish teacher vs. student features simultaneously. Extensive experiments on CIFAR-10/100, SVHN, ImageNet and iMaterialist Challenge Dataset demonstrate the effectiveness of our MEAL method. On ImageNet, our ResNet-50 based MEAL achieves top-1/5 21.79%/5.99% val error, which outperforms the original model by 2.06%/1.14%. On iMaterialist Challenge Dataset, our MEAL obtains a remarkable improvement of top-3 1.15% (official evaluation metric) on a strong baseline model of ResNet-101.


title: On the dichotomy of a locally compact semitopological monoid of order isomorphisms between principal filters of $\mathbb{N}^n$ with adjoined zero

abstract: Let $n$ be any positive integer and $\mathscr{I\!P\!F}(\mathbb{N}^n)$ be the semigroup of all order isomorphisms between principal filters of the $n$-th power of the set of positive integers $\mathbb{N}$ with the product order. We prove that a Hausdorff locally compact semitopological semigroup ${\mathscr{I\!P\!F}(\mathbb{N}^n)}$ with an adjoined zero is either compact or discrete.


title: Compositional Video Prediction

abstract: We present an approach for pixel-level future prediction given an input image of a scene. We observe that a scene is comprised of distinct entities that undergo motion and present an approach that operationalizes this insight. We implicitly predict future states of independent entities while reasoning about their interactions, and compose future video frames using these predicted states. We overcome the inherent multi-modality of the task using a global trajectory-level latent random variable, and show that this allows us to sample diverse and plausible futures. We empirically validate our approach against alternate representations and ways of incorporating multi-modality. We examine two datasets, one comprising of stacked objects that may fall, and the other containing videos of humans performing activities in a gym, and show that our approach allows realistic stochastic video prediction across these diverse settings. See https://judyye.github.io/CVP/ for video predictions.


title: Jackiw-Teitelboim Model Coupled to Conformal Matter in the Semi-Classical Limit

abstract: We analyse the Jackiw-Teitelboim model of 2D gravity coupled to $N$ massless free scalar fields in the semi-classical limit. Two systems are studied which essentially differ in the boundary conditions that are imposed. We find that the thermodynamics has interesting differences. We also analyse the response to additional infalling matter which satisfies the null energy condition. The second law is shown to be valid in both systems for the generalised entropy which takes into account the entanglement across the event horizon due to the matter fields. Similarly we find that the generalised entropy increases along future Q-screens in both systems.


title: Dynamical resolution scale in transverse momentum distributions at the LHC

abstract: The QCD evolution of transverse momentum dependent (TMD) distribution functions has recently been formulated in a parton branching (PB) formalism. In this approach, soft-gluon coherence effects are taken into account by introducing the soft-gluon resolution scale and exploiting the relation between transverse-momentum recoils and branching scales. In this work we investigate the implications of dynamical, i.e., branching scale dependent, resolution scales. We present both analytical studies and numerical solution of PB evolution equations in the presence of dynamical resolution scales. We use this to compare PB results with other approaches in the literature, and to analyze predictions for transverse momentum distributions in $Z$-boson production at the Large Hadron Collider (LHC).


title: Some inequalities for reversible Markov chains and branching random walks via spectral optimization

abstract: We present results relating mixing times to the intersection time of branching random walk (BRW) in which the logarithm of the expected number of particles grows at rate of the spectral-gap $\mathrm{gap}$ . This is a finite state space analog of a critical branching process. Namely, we show that the maximal expected hitting time of a state by such a BRW is up to a universal constant larger than the $L_{\infty}$ mixing-time, whereas under transitivity the same is true for the intersection time of two independent such BRWs. Using the same methodology, we show that for a sequence of reversible Markov chains, the $L_{\infty}$ mixing-times $t_{\mathrm{mix}}^{(\infty)} $ are of smaller order than the maximal hitting times $t_{\mathrm{hit}}$ iff the product of the spectral-gap and $t_{\mathrm{hit}}$ diverges, by establishing the inequality $t_{\mathrm{mix}}^{(\infty)} \le \frac{1}{\mathrm{gap}}\log(et_{\mathrm{hit}} \cdot \mathrm{gap}) $. This resolves a conjecture of Aldous and Fill (Reversible Markov chains and random walks on graphs, Open Problem 14.12) asserting that under transitivity the condition that $ t_{\mathrm{hit}} \gg \frac{1}{\mathrm{gap}} $ implies mean-field behavior for the coalescing time of coalescing random walks.


title: Double Reinforcement Learning for Efficient Off-Policy Evaluation in Markov Decision Processes

abstract: Off-policy evaluation (OPE) in reinforcement learning allows one to evaluate novel decision policies without needing to conduct exploration, which is often costly or otherwise infeasible. We consider for the first time the semiparametric efficiency limits of OPE in Markov decision processes (MDPs), where actions, rewards, and states are memoryless. We show existing OPE estimators may fail to be efficient in this setting. We develop a new estimator based on cross-fold estimation of $q$-functions and marginalized density ratios, which we term double reinforcement learning (DRL). We show that DRL is efficient when both components are estimated at fourth-root rates and is also doubly robust when only one component is consistent. We investigate these properties empirically and demonstrate the performance benefits due to harnessing memorylessness.


title: ViCo: Word Embeddings from Visual Co-occurrences

abstract: We propose to learn word embeddings from visual co-occurrences. Two words co-occur visually if both words apply to the same image or image region. Specifically, we extract four types of visual co-occurrences between object and attribute words from large-scale, textually-annotated visual databases like VisualGenome and ImageNet. We then train a multi-task log-bilinear model that compactly encodes word "meanings" represented by each co-occurrence type into a single visual word-vector. Through unsupervised clustering, supervised partitioning, and a zero-shot-like generalization analysis we show that our word embeddings complement text-only embeddings like GloVe by better representing similarities and differences between visual concepts that are difficult to obtain from text corpora alone. We further evaluate our embeddings on five downstream applications, four of which are vision-language tasks. Augmenting GloVe with our embeddings yields gains on all tasks. We also find that random embeddings perform comparably to learned embeddings on all supervised vision-language tasks, contrary to conventional wisdom.


title: Unsupervised Lemmatization as Embeddings-Based Word Clustering

abstract: We focus on the task of unsupervised lemmatization, i.e. grouping together inflected forms of one word under one label (a lemma) without the use of annotated training data. We propose to perform agglomerative clustering of word forms with a novel distance measure. Our distance measure is based on the observation that inflections of the same word tend to be similar both string-wise and in meaning. We therefore combine word embedding cosine similarity, serving as a proxy to the meaning similarity, with Jaro-Winkler edit distance. Our experiments on 23 languages show our approach to be promising, surpassing the baseline on 23 of the 28 evaluation datasets.


title: Sequential Latent Spaces for Modeling the Intention During Diverse Image Captioning

abstract: Diverse and accurate vision+language modeling is an important goal to retain creative freedom and maintain user engagement. However, adequately capturing the intricacies of diversity in language models is challenging. Recent works commonly resort to latent variable models augmented with more or less supervision from object detectors or part-of-speech tags. Common to all those methods is the fact that the latent variable either only initializes the sentence generation process or is identical across the steps of generation. Both methods offer no fine-grained control. To address this concern, we propose Seq-CVAE which learns a latent space for every word position. We encourage this temporal latent space to capture the 'intention' about how to complete the sentence by mimicking a representation which summarizes the future. We illustrate the efficacy of the proposed approach to anticipate the sentence continuation on the challenging MSCOCO dataset, significantly improving diversity metrics compared to baselines while performing on par w.r.t sentence quality.


title: VL-BERT: Pre-training of Generic Visual-Linguistic Representations

abstract: We introduce a new pre-trainable generic representation for visual-linguistic tasks, called Visual-Linguistic BERT (VL-BERT for short). VL-BERT adopts the simple yet powerful Transformer model as the backbone, and extends it to take both visual and linguistic embedded features as input. In it, each element of the input is either of a word from the input sentence, or a region-of-interest (RoI) from the input image. It is designed to fit for most of the visual-linguistic downstream tasks. To better exploit the generic representation, we pre-train VL-BERT on the massive-scale Conceptual Captions dataset, together with text-only corpus. Extensive empirical analysis demonstrates that the pre-training procedure can better align the visual-linguistic clues and benefit the downstream tasks, such as visual commonsense reasoning, visual question answering and referring expression comprehension. It is worth noting that VL-BERT achieved the first place of single model on the leaderboard of the VCR benchmark. Code is released at \url{https://github.com/jackroos/VL-BERT}.


title: Semi-Lagrangian implicit Bhatnagar-Gross-Krook collision model for the finite-volume discrete Boltzmann method

abstract: A new implicit BGK collision model using a semi-Lagrangian approach is proposed in this paper. Unlike existing models, in which the implicit BGK collision is resolved either by a temporal extrapolation or by a variable transformation, the new model removes the implicitness by tracing the particle distribution functions (PDFs) back in time along their characteristic paths during the collision process. An interpolation scheme is needed to evaluate the PDFs at the traced-back locations. By using the first-order interpolation, the resulting model allows for the straightforward replacement of ${f_{\alpha}}^{eq,n+1}$ by ${f_{\alpha}}^{eq,n}$ no matter where it appears. After comparing the new model with the existing models under different numerical conditions (e.g. different flux schemes and time marching schemes) and using the new model to successfully modify the variable transformation technique, three conclusions can be drawn. First, the new model can improve the accuracy by almost an order of magnitude. Second, it can slightly reduce the computational cost. Therefore, the new scheme improves accuracy without extra cost. Finally, the new model can significantly improve the ${\Delta}t/{\tau}$ limit compared to the temporal interpolation model while having the same ${\Delta}t/{\tau}$ limit as the variable transformation approach. The new scheme with a second-order interpolation is also developed and tested; however, that technique displays no advantage over the simple first-order interpolation approach. Both numerical and theoretical analyses are also provided to explain why the new implicit scheme with simple first-order interpolation can outperform the same scheme with second-order interpolation, as well as the existing temporal extrapolation and variable transformation schemes.


title: Design, Assembly, Calibration, and Measurement of an Augmented Reality Haploscope

abstract: A haploscope is an optical system which produces a carefully controlled virtual image. Since the development of Wheatstone's original stereoscope in 1838, haploscopes have been used to measure perceptual properties of human stereoscopic vision. This paper presents an augmented reality (AR) haploscope, which allows the viewing of virtual objects superimposed against the real world. Our lab has used generations of this device to make a careful series of perceptual measurements of AR phenomena, which have been described in publications over the previous 8 years. This paper systematically describes the design, assembly, calibration, and measurement of our AR haploscope. These methods have been developed and improved in our lab over the past 10 years. Despite the fact that 180 years have elapsed since the original report of Wheatstone's stereoscope, we have not previously found a paper that describes these kinds of details.


title: First i-TED demonstrator: a Compton imager with Dynamic Electronic Collimation

abstract: i-TED consists of both a total energy detector and a Compton camera primarily intended for the measurement of neutron capture cross sections by means of the simultaneous combination of neutron time-of-flight (TOF) and $\gamma$-ray imaging techniques. TOF allows one to obtain a neutron-energy differential capture yield, whereas the imaging capability is intended for the discrimination of radiative background sources, that have a spatial origin different from that of the capture sample under investigation. A distinctive feature of i-TED is the embedded Dynamic Electronic Collimation (DEC) concept, which allows for a trade-off between efficiency and image resolution. Here we report on some general design considerations and first performance characterization measurements made with an i-TED demonstrator in order to explore its $\gamma$-ray detection and imaging capabilities.


title: A Proof of Brouwer's Conjecture

abstract: Building up from recent results by other authors, we present proof for Brouwer's Conjecture.


title: A Class D Power Amplifier for Multi-Frequency Eddy Current Testing Based on Multi-Simultaneous-Frequency Selective Harmonic Elimination Pulse Width Modulation

abstract: Efficiency and multisimultaneous-frequency (MSF) output capability are two major criteria characterizing the performance of a power amplifier in the application of multifrequency eddy current testing (MECT). Switch-mode power amplifiers are known to have a very high efficiency, yet they have rarely been adopted in the instrumental development of MECT. In addition, switch-mode power amplifiers themselves are lacking in the research literature for MSF capability. In this article, a Class D power amplifier is designed so as to address the two issues. An MSF selective harmonic elimination pulsewidth modulation method is proposed to generate alternating magnetic fields, which are rich in selected harmonics. A field-programmable-gate-array-based experimental system has been developed to verify the design. Results show that the proposed methodology is capable of generating high MSF currents in the transmitting coil with a low distortion of signal.


title: Electric Circuit Realizations of Fracton Physics

abstract: We design a set of classical macroscopic electric circuits in which charge exhibits the mobility restrictions of fracton quasiparticles. The crucial ingredient in these circuits is a transformer, which induces currents between pairs of adjacent wires. For an appropriately designed geometry, this induction serves to enforce conservation of dipole moment. We show that a network of capacitors connected via ideal transformers will forever remember the dipole moment of its initial charge configuration. Relaxation of the dipole moment in realistic systems can only occur via flux leakage in the transformers, which will lead to violations of fracton physics at the longest times. We propose a concrete diagnostic for these "fractolectric" circuits in the form of their characteristic equilibrium charge configurations, which we verify using simple circuit simulation software. These circuits not only provide an experimental testing ground for fracton physics, but also serve as DC filters. We outline extensions of these ideas to circuits featuring other types of higher moment conservation laws, as well as to higher-dimensional circuits which act as fracton "current-ice." While our focus is on classical circuits, we discuss how these ideas can be straightforwardly extended to realize quantized fractons in superconducting circuits.


title: The tidal evolution of dark matter substructure -- I. Subhalo density profiles

abstract: Accurately predicting the abundance and structural evolution of dark matter subhaloes is crucial for understanding galaxy formation, modeling galaxy clustering, and constraining the nature of dark matter. Due to the nonlinear nature of subhalo evolution, cosmological $N$-body simulations remain its primary method of investigation. However, it has recently been demonstrated that such simulations are still heavily impacted by artificial disruption, diminishing the information content on small scales and reducing the reliability of all simulation-calibrated semi-analytical models. In this paper, we utilize the recently released DASH library of high-resolution, idealized simulations of the tidal evolution of subhaloes, which are unhindered by numerical overmerging due to discreteness noise or force softening, to calibrate an improved, more-accurate model of the evolution of the density profiles of subhaloes that undergo tidal heating and stripping within their host halo. By testing previous findings that the structural evolution of a tidally truncated subhalo depends solely on the fraction of mass stripped, independent of the details of the stripping, we identify an additional dependence on the initial subhalo concentration. We provide significantly improved fitting functions for the subhalo density profiles and structural parameters ($V_\mathrm{max}$ and $r_\mathrm{max}$) that are unimpeded by numerical systematics and applicable to a wide range of parameter space. This model will be an integral component of a future semi-analytical treatment of substructure evolution, which can be used to predict key quantities, such as the evolved subhalo mass function and annihilation boost factors, and validate such calculations performed with cosmological simulations.


title: On the black hole content and initial mass function of 47 Tuc

abstract: The globular cluster (GC) 47 Tuc has recently been proposed to host an intermediate-mass black hole (IMBH) or a population of stellar-mass black holes (BHs). To shed light on its dark content, we present an application of self-consistent multimass models with a varying mass function and content of stellar remnants, which we fit to various observational constraints. Our best-fitting model successfully matches the observables and correctly predicts the radial distribution of millisecond pulsars and their gravitational accelerations inferred from long-term timing observations. The data favours a population of BHs with a total mass of $430^{+386}_{-301}$ $M_{\odot}$, but the most likely model has very few BHs. Since our models do not include a central IMBH and accurately reproduce the observations, we conclude that there is currently no need to invoke the presence of an IMBH in 47 Tuc. The global present-day mass function inferred is significantly depleted in low-mass stars (power-law slope $\alpha=-0.52^{+0.17}_{-0.16}$). Given the orbit and predicted mass-loss history of this massive GC, the dearth of low-mass stars is difficult to explain with a standard initial mass function (IMF) followed by long-term preferential escape of low-mass stars driven by two-body relaxation, and instead suggests that 47 Tuc may have formed with a bottom-light IMF. We discuss alternative evolutionary origins for the flat mass function and ways to reconcile this with the low BH retention fraction. Finally, by capturing the effect of dark remnants, our method offers a new way to probe the IMF in a GC above the current main-sequence turn-off mass, for which we find a slope of $-2.49\pm0.08$.


title: Hydro+ in Action: Understanding the Out-of-Equilibrium Dynamics Near a Critical Point in the QCD Phase Diagram

abstract: Upcoming experimental programs will look for signatures of a possible critical point in the QCD phase diagram in fluctuation observables. To understand and predict these signatures, one must account for the fact that the dynamics of any critical fluctuations must be out-of-equilibrium: because of critical slowing down, the fluctuations cannot stay in equilibrium as the droplet of QGP produced in a collision expands and cools. Furthermore, their out-of-equilibrium dynamics must also influence the hydrodynamic evolution of the cooling droplet. The recently developed Hydro+ formalism allows for a consistent description of both the hydrodynamics and the out-of-equilibrium fluctuations, including the feedback between them. We shall explicitly demonstrate how this works, setting up a Hydro+ simulation in a simplified setting: a rapidity-independent fireball undergoing radial flow with an equation of state in which we imagine a critical point close to the $\mu_B=0$ axis of the phase diagram. Within this setup, we show that we can quantitatively capture non-equilibrium phenomena, including critical fluctuations over a range of scales and memory effects. Furthermore, we illustrate the interplay between the dynamics of the fluctuations and the hydrodynamic flow of the fireball: as the fluid cools and flows, the dynamical fluctuations lag relative to how they would evolve if they stayed in equilibrium; there is then a backreaction on the flow itself due to the out-of-equilibrium fluctuations; and, in addition, the radial flow transports fluctuations outwards by advection. Within our model, we find that the backreaction from the out-of-equilibrium fluctuations does not yield dramatically large effects in the hydrodynamic variables. Further work will be needed in order to check this quantitative conclusion in other settings but, if it persists, this will considerably simplify future modelling.


title: Emergent fractons and algebraic quantum liquid from plaquette melting transitions

abstract: Paramagnetic spin systems with spontaneously broken spatial symmetries, such as valence bond solid (VBS) phases, can host topological defects carrying non-trivial quantum numbers, which enables the paradigm of deconfined quantum criticality. In this work, we study the properties of topological defects in valence plaquette solid (VPS) phases on square and cubic lattices. We show that the defects of the VPS order parameter, in addition to possessing non-trivial quantum numbers, have fracton mobility constraints deep in the VPS phase, which has been overlooked previously. The spinon inside a single vortex cannot move freely in any direction, while a dipolar pair of vortices with spinon pairs can only move perpendicular to its dipole moment. These mobility constraints, while they persist, can potentially inhibit the condensation of vortices and preclude a continuous transition from the VPS to the N\'eel antiferromagnet. Instead, the VPS melting transition can be driven by proliferation of spinon dipoles. For example, we argue that a $2d$ VPS can melt into a stable gapless phase in the form of an algebraic bond liquid with algebraic correlations and long range entanglement. Such a bond liquid phase yields a concrete example of the elusive $2d$ Bose metal with symmetry fractionalization. We also study $3d$ valence plaquette and valence cube ordered phase, and demonstrate that the topological defects therein also have fractonic dynamics. Possible nearby phases after melting the valence plaquettes or cubes are also discussed.


title: Crystallographic splitting theorem for band representations and fragile topological photonic crystals

abstract: The fundamental building blocks in band theory are band representations (BRs): bands whose infinitely-numbered Wannier functions are generated (by action of a space group) from a finite number of symmetric Wannier functions centered on a point in space. This work aims to simplify questions on a multi-rank BR by splitting it into unit-rank bands, via the following crystallographic splitting theorem: being a rank-$N$ BR is equivalent to being splittable into a finite sum of bands indexed by $\{1,2,\ldots,N\}$, such that each band is spanned by a single, analytic Bloch function of $k$, and any symmetry in the space group acts by permuting $\{1,2,\ldots,N\}$. Applying this theorem, we develop computationally efficient methods to determine whether a given energy band (of a tight-binding or Schr\"odinger Hamiltonian) is a BR, and, if so, how to numerically construct the corresponding symmetric Wannier functions. Thus we prove that rotation-symmetric topological insulators in class AI are fragile, meaning that the obstruction to symmetric Wannier functions is removable by addition of BRs. An implication of fragility is that its boundary states, while robustly covering the bulk energy gap in finite-rank tight-binding models, are unstable if the Hilbert space is expanded to include all symmetry-allowed representations. These fragile insulators have photonic analogs that we identify; in particular, we prove that an existing photonic crystal built by Yang et al. [Nature 565, 622 (2019)] is fragile topological with removable surface states, which disproves a widespread perception of 'topologically-protected' surface states in time-reversal-invariant, gapped photonic/phononic crystals. Our theorem is finally applied to derive various symmetry obstructions on the Wannier functions of topological insulators, and to prove their equivalence with the nontrivial holonomy of Bloch functions.


title: Exploring the Space of Jets with CMS Open Data

abstract: We explore the metric space of jets using public collider data from the CMS experiment. Starting from 2.3/fb of 7 TeV proton-proton collisions collected at the Large Hadron Collider in 2011, we isolate a sample of 1,690,984 central jets with transverse momentum above 375 GeV. To validate the performance of the CMS detector in reconstructing the energy flow of jets, we compare the CMS Open Data to corresponding simulated data samples for a variety of jet kinematic and substructure observables. Even without detector unfolding, we find very good agreement for track-based observables after using charged hadron subtraction to mitigate the impact of pileup. We perform a range of novel analyses, using the "energy mover's distance" (EMD) to measure the pairwise difference between jet energy flows. The EMD allows us to quantify the impact of detector effects, visualize the metric space of jets, extract correlation dimensions, and identify the most and least typical jet configurations. To facilitate future jet studies with CMS Open Data, we make our datasets and analysis code available, amounting to around two gigabytes of distilled data and one hundred gigabytes of simulation files.


title: A Real-Time Search for Interstellar Impacts on the Moon

abstract: The discovery of `Oumuamua and CNEOS 2014-01-08 allowed for a calibration of the impact rate of interstellar objects. We propose a new telescope in lunar orbit to study in real-time interstellar meteoroid impacts and to serve as a laboratory for hypervelocity collisions. We show that a telescope with diameter $D_a \gtrsim 2 \mathrm{\; m}$ would be capable of detecting $\gtrsim 1$ interstellar meteoroid impacts (among hundreds of Solar System meteoroid impacts) per year. For each meteoroid, measurements of the reflected sunlight and shadow, as well as the impact's optical flash and crater, would allow for the determination of the 3D velocity, mass, density, and composition, as well as the radiative efficiency.


title: Structure of Brightest Cluster Galaxies and Intracluster Light

abstract: Observations of 170 local ($z\lesssim0.08$) galaxy clusters in the northern hemisphere have been obtained with the Wendelstein Telescope Wide Field Imager (WWFI). We correct for systematic effects such as point-spread function broadening, foreground star contamination, relative bias offsets, and charge persistence. Background inhomogeneities induced by scattered light are reduced down to $\Delta {\rm SB} > 31~g'$ mag arcsec$^{-2}$ by large dithering and subtraction of night-sky flats. Residual background inhomogeneities brighter than ${\rm SB}_{\sigma}< 27.6~g'$ mag arcsec$^{-2}$ caused by galactic cirrus are detected in front of 23% of the clusters. However, the large field of view allows discrimination between accretion signatures and galactic cirrus. We detect accretion signatures in the form of tidal streams in 22%, shells in 9.4%, and multiple nuclei in 47% of the Brightest Cluster Galaxies (BCGs) and find two BCGs in 7% of the clusters. We measure semimajor-axis surface brightness profiles of the BCGs and their surrounding Intracluster Light (ICL) down to a limiting surface brightness of ${\rm SB} = 30~g'$ mag arcsec$^{-2}$. The spatial resolution in the inner regions is increased by combining the WWFI light profiles with those that we measured from archival \textit{Hubble Space Telescope} images or deconvolved WWFI images. We find that 71% of the BCG+ICL systems have surface brightness (SB) profiles that are well described by a single S\'ersic (SS) function, whereas 29% require a double S\'ersic (DS) function to obtain a good fit. We find that BCGs have scaling relations that differ markedly from those of normal ellipticals, likely due to their indistinguishable embedding in the ICL.


title: The VMC Survey -- XXXIV. Morphology of Stellar Populations in the Magellanic Clouds

abstract: The Magellanic Clouds are nearby dwarf irregular galaxies whose morphologies show different properties when traced by different stellar populations, making them an important laboratory for studying galaxy morphologies. We study the morphology of the Magellanic Clouds using data from the VISTA survey of the Magellanic Clouds system (VMC). We used about $10$ and $2.5$ million sources across an area of $\sim105$ deg$^2$ and $\sim42$ deg$^2$ towards the Large and Small Magellanic Cloud (LMC and SMC), respectively. We estimated median ages of stellar populations occupying different regions of the near-infrared ($J-K_\mathrm{s}, K_\mathrm{s}$) colour-magnitude diagram. Morphological maps were produced and detailed features in the central regions were characterised for the first time with bins corresponding to a spatial resolution of $0.13$ kpc (LMC) and $0.16$ kpc (SMC). In the LMC, we find that main sequence stars show coherent structures that grow with age and trace the multiple spiral arms of the galaxy, star forming regions become dimmer as we progress in age, while supergiant stars are centrally concentrated. Intermediate-age stars, despite tracing a regular and symmetrical morphology, show central clumps and hints of spiral arms. In the SMC, young main sequence stars depict a broken bar. Intermediate-age populations show signatures of elongation towards the Magellanic Bridge that can be attributed to the LMC-SMC interaction $\sim200$ Myr ago. They also show irregular central features suggesting that the inner SMC has also been influenced by tidal interactions.


title: IceCube Search for Galactic Neutrino Sources based on HAWC Observations of the Galactic Plane

abstract: We present a search in IceCube data for neutrino emission from Galactic TeV gamma-ray sources detected by the HAWC gamma-ray observatory. HAWC serves as the excellent instrument to complement IceCube with its energy range extending to very high energies. Assuming that the highest energy photons originate from the decay of pions, rather than from accelerated leptons, the very high energy gamma-rays observed by HAWC are expected to be correlated with neutrinos. Using eight years of IceCube data, we report on two analyses that investigate a possible neutrino--gamma ray correlation. The first is a stacked analysis of identified HAWC point sources and the second is a template method which accounts for the full morphology of HAWC sources, including their measured extension.


title: Search for Neutrinos from Populations of Optical Transients

abstract: Since the detection of high-energy cosmic neutrinos at the IceCube Neutrino Observatory in 2013, there has been an on-going search to find the origins of this flux. Despite recent evidence identifying a flaring blazar as a possible neutrino source, the vast majority of the diffuse neutrino flux measured by IceCube remains unexplained. Here, the latest IceCube results testing time-dependent correlation between neutrinos and Tidal Disruption Events (TDEs) are presented, limiting the contribution of jetted and non-jetted TDEs to the diffuse astrophysical neutrino flux to be less than 1.3% and 26% respectively at 90% confidence level. In addition, a dedicated search for neutrinos from the extraordinary transient AT2018cow are presented, and upper limits on the integrated neutrino emission are derived. Expected improvements from new and upcoming time domain optical surveys (such as ZTF and LSST) are also introduced.


title: Estimation of singly-transiting K2 planet periods with Gaia parallaxes

abstract: When a planet is only observed to transit once, direct measurement of its period is impossible. It is possible, however, to constrain the periods of single transiters, and this is desirable as they are likely to represent the cold and far extremes of the planet population observed by any particular survey. Improving the accuracy with which the period of single transiters can be constrained is therefore critical to enhance the long-period planet yield of surveys. Here, we combine Gaia parallaxes with stellar models and broad-band photometry to estimate the stellar densities of K2 planet host stars, then use that stellar density information to model individual planet transits and infer the posterior period distribution. We show that the densities we infer are reliable by comparing with densities derived through asteroseismology, and apply our method to 27 validation planets of known (directly measured) period, treating each transit as if it were the only one, as well as to 12 true single transiters. When we treat eccentricity as a free parameter, we achieve a fractional period uncertainty over the true single transits of $94^{+87}_{-58}\%$, and when we fix $e=0$, we achieve fractional period uncertainty $15^{+30}_{-6}\%$, a roughly threefold improvement over typical period uncertainties of previous studies.


title: Evolution of Neutral Oxygen During the Epoch of Reionization and its Use in Estimating the Neutral Hydrogen Fraction

abstract: We use synthetic sightlines drawn through snapshots of the Technicolor Dawn simulations to explore how the statistics of neutral oxygen OI absorbers respond to hydrogen reionization. The ionization state of the circumgalactic medium (CGM) initially roughly tracks that of the intergalactic medium, but beginning at $z=8$ the CGM grows systematically more neutral owing to self-shielding. Weak absorbers trace diffuse gas that lies farther from halos, hence they are ionized first, whereas stronger systems are less sensitive to reionization. The overall OI covering fraction decreases slowly with time owing to competition between ongoing enrichment and gradual encroachment of ionization fronts into increasingly overdense gas. While the declining covering fraction is partially offset by continued formation of new halos, the ionization of the diffuse gas causes the predicted line-of-sight incidence rate of OI absorbers to decline abruptly at the overlap epoch, in qualitative agreement with observations. In comparison to the recently-observed equivalent width (EW) distribution at $z\approx6$, the simulations underproduce systems with $EW \geq 0.1 \unicode{x212B}$, although they reproduce weaker systems with $EW \geq 0.05 \unicode{x212B}$. By $z\approx5$, the incidence of $EW < 0.1 \unicode{x212B}$ systems are overproduced, consistent with previous indications that the simulated ionizing background is too weak at $z<6$. The summed column densities of SiII and SiIV trace the total oxygen column, and hence the ratio of the OI and SiII + SiIV comoving mass densities traces the progress of reionization. This probe may prove particularly useful in the regime where $x_{HI} > 10\%$


title: Decay of correlations for certain isometric extensions of Anosov flows

abstract: We establish exponential decay of correlations of all orders for locally $G$-accessible isometric extensions of transitive Anosov flows, under the assumption that the strong stable and strong unstable foliations of the base Anosov flow are jointly $C^1$. This is accomplished by translating accessibility properties of the extension into local non-integrability estimates measured by Dolgopyat's infinitesimal transitivity group, from which we obtain contraction properties for a class of 'twisted' symbolic transfer operators.


title: Clustered Hierarchical Entropy-Scaling Search of Astronomical and Biological Data

abstract: Both astronomy and biology are experiencing explosive growth of data, resulting in a "big data" problem that stands in the way of a "big data" opportunity for discovery. One common question asked of such data is that of approximate search ($\rho-$nearest neighbors search). We present a hierarchical search algorithm for such data sets that takes advantage of particular geometric properties apparent in both astronomical and biological data sets, namely the metric entropy and fractal dimensionality of the data. We present CHESS (Clustered Hierarchical Entropy-Scaling Search), a search tool with virtually no loss in specificity or sensitivity, demonstrating a $13.6\times$ speedup over linear search on the Sloan Digital Sky Survey's APOGEE data set and a $68\times$ speedup on the GreenGenes 16S metagenomic data set, as well as asymptotically fewer distance comparisons on APOGEE when compared to the FALCONN locality-sensitive hashing library. CHESS demonstrates an asymptotic complexity not directly dependent on data set size, and is in practice at least an order of magnitude faster than linear search by performing fewer distance comparisons. Unlike locality-sensitive hashing approaches, CHESS can work with any user-defined distance function. CHESS also allows for implicit data compression, which we demonstrate on the APOGEE data set. We also discuss an extension allowing for efficient k-nearest neighbors search.


title: An extremely X--ray weak blazar at z=5

abstract: We present the discovery and properties of DESJ014132.4-542749.9 (DES0141-54), a new powerful radio-loud active galactic nucleus (AGN) in the early Universe (z=5.0). It was discovered by cross-matching the first data release of the Dark Energy Survey (DES DR1) with the Sidney University Molonglo Survey (SUMSS) radio catalog at 0.843 GHz. This object is the first radio-loud AGN at high redshift discovered in the DES. The radio properties of DES0141-54, namely its very large radio-loudness (R>10$^{4}$), the high radio luminosity (L$_{0.8 GHz}$=1.73$\times$10$^{28}$ W Hz$^{-1}$), and the flatness of the radio spectrum ($\alpha$=0.35) up to very high frequencies (120 GHz in the source's rest frame), classify this object as a blazar, meaning, a radio-loud AGN observed along the relativistic jet axis. However, the X--ray luminosity of DESJ0141-54 is much lower compared to those of the high redshift (z$\geq$4.5) blazars discovered so far. Moreover its X-ray-to-radio luminosity ratio (log($\frac{L_{[0.5-10]keV}}{L_{1.4GHz}}$)=9.96$\pm$0.30 Hz) is small also when compared to lower redshift blazars: only 2\% of the low-z population has a similar ratio. By modeling the spectral energy distribution we found that this peculiar X--ray weakness and the powerful radio emission could be related to a particularly high value of the magnetic field. Finally, the mass of the central black hole is relatively small (M$_{BH}$ = 3-8 $\times$10$^8$ M$_{\odot}$) compared to other confirmed blazars at similar redshift, making DES0141-54 the radio-loud AGN that host the smallest supermassive black hole ever discovered at z$\geq$5.


title: Simulation of Quantum Many-Body Systems on Amazon Cloud

abstract: Quantum many-body systems (QMBs) are some of the most challenging physical systems to simulate numerically. Methods involving approximations for tensor network (TN) contractions have proven to be viable alternatives to algorithms such as quantum Monte Carlo or simulated annealing. However, these methods are cumbersome, difficult to implement, and often have significant limitations in their accuracy and efficiency when considering systems in more than one dimension. In this paper, we explore the exact computation of TN contractions on two-dimensional geometries and present a heuristic improvement of TN contraction that reduces the computing time, the amount of memory, and the communication time. We run our algorithm for the Ising model using memory optimized x1.32x large instances on Amazon Web Services (AWS) Elastic Compute Cloud (EC2). Our results show that cloud computing is a viable alternative to supercomputers for this class of scientific applications.