title: Sherman-Morrison-Woodbury Identity for Tensors

abstract: In linear algebra, the sherman-morrison-woodbury identity says that the inverse of a rank-$k$ correction of some matrix can be computed by doing a rank-k correction to the inverse of the original matrix. This identity is crucial to accelerate the matrix inverse computation when the matrix involves correction. Many scientific and engineering applications have to deal with this matrix inverse problem after updating the matrix, e.g., sensitivity analysis of linear systems, covariance matrix update in kalman filter, etc. However, there is no similar identity in tensors. In this work, we will derive the sherman-morrison-woodbury identity for invertible tensors first. Since not all tensors are invertible, we further generalize the sherman-morrison-woodbury identity for tensors with moore-penrose generalized inverse by utilizing orthogonal projection of the correction tensor part into the original tensor and its Hermitian tensor. According to this new established the sherman-morrison-woodbury identity for tensors, we can perform sensitivity analysis for multi-linear systems by deriving the normalized upper bound for the solution of a multilinear system. Several numerical examples are also presented to demonstrate how the normalized error upper bounds are affected by perturbation degree of tensor coefficients.


title: Serre functors and graded categories

abstract: We study Serre structures on categories enriched in pivotal monoidal categories, and apply this to study Serre structures on two types of graded k-linear categories: categories with group actions and categories with graded hom spaces. We check that Serre structures are preserved by taking orbit categories and skew group categories, and describe the relationship with graded Frobenius algebras. Using a formal version of Auslander-Reiten translations, we show that the derived category of a d-representation finite algebra is fractionally Calabi-Yau if and only if its preprojective algebra has a graded Nakayama automorphism of finite order. This connects various results in the literature and gives new examples of fractional Calabi-Yau algebras.


title: Image-based Vehicle Re-identification Model with Adaptive Attention Modules and Metadata Re-ranking

abstract: Vehicle Re-identification is a challenging task due to intra-class variability and inter-class similarity across non-overlapping cameras. To tackle these problems, recently proposed methods require additional annotation to extract more features for false positive image exclusion. In this paper, we propose a model powered by adaptive attention modules that requires fewer label annotations but still out-performs the previous models. We also include a re-ranking method that takes account of the importance of metadata feature embeddings in our paper. The proposed method is evaluated on CVPR AI City Challenge 2020 dataset and achieves mAP of 37.25% in Track 2.


title: Addressing the interpretability problem for deep learning using many valued quantum logic

abstract: Deep learning models are widely used for various industrial and scientific applications. Even though these models have achieved considerable success in recent years, there exists a lack of understanding of the rationale behind decisions made by such systems in the machine learning community. This problem of interpretability is further aggravated by the increasing complexity of such models. This paper utilizes concepts from machine learning, quantum computation and quantum field theory to demonstrate how a many valued quantum logic system naturally arises in a specific class of generative deep learning models called Convolutional Deep Belief Networks. It provides a robust theoretical framework for constructing deep learning models equipped with the interpretability of many valued quantum logic systems without compromising their computing efficiency.


title: A Machine Learning Pipeline Stage for Adaptive Frequency Adjustment

abstract: A machine learning (ML) design framework is proposed for adaptively adjusting clock frequency based on propagation delay of individual instructions. A random forest model is trained to classify propagation delays in real time, utilizing current operation type, current operands, and computation history as ML features. The trained model is implemented in Verilog as an additional pipeline stage within a baseline processor. The modified system is experimentally tested at the gate level in 45 nm CMOS technology, exhibiting a speedup of 70% and energy reduction of 30% with coarse-grained ML classification. A speedup of 89% is demonstrated with finer granularities with 15.5% reduction in energy consumption.


title: Regularization of the movement of a material point along a flat trajectory: application to robotics problems

abstract: The control problem of the working tool movement along a predefined trajectory is considered. The integral of kinetic energy and weighted inertia forces for the whole period of motion is considered as a cost functional. The trajectory is assumed to be planar and defined in advance. The problem is reduced to a system of ordinary differential equations of the fourth order. Numerical examples of solving the problem for movement along straight, circular and elliptical trajectories are given.


title: Construction of weak solutions to compressible Navier--Stokes equations with general inflow/outflow boundary conditions via a numerical approximation

abstract: The construction of weak solutions to compressible Navier-Stokes equations via a numerical method (including a rigorous proof of the convergence) is in a short supply, and so far, available only for one sole numerical scheme suggested in Karper [{\em Numer. Math.}, 125(3) : 441--510, 2013] for the no slip boundary conditions and the isentropic pressure with adiabatic coefficient $\gamma>3$. Here we consider the same problem for the general non zero inflow-outflow boundary conditions, which is definitely more appropriate setting from the point of view of applications, but which is essentially more involved as far as the existence of weak solutions is concerned. There is a few recent proofs of existence of weak solutions in this setting, but none of them is performed via a numerical method. The goal of this paper is to fill this gap. The existence of weak solutions on the continuous level requires several tools of functional and harmonic analysis and differential geometry whose numerical counterparts are not known. Our main strategy therefore consists in rewriting of the numerical scheme in its variational form modulo remainders and to apply and/or to adapt to the new variational formulation the tools developed in the theoretical analysis. In addition to the result, which is new, the synergy between numerical and theoretical analysis is the main originality of the present paper.


title: WordPress on AWS: a Communication Framework

abstract: Every organization needs to communicate with its audience, and social media is an attractive and inexpensive way to maintain dialogic communication. About 1/3 of the Internet web pages are powered by WordPress, and about a million companies have moved their IT infrastructure to the AWS cloud. Together, AWS and WordPress offer an attractive, effective and inexpensive way for companies, both large and small, to maintain their presence on the web.


title: Star versions of Lindel\"of spaces

abstract: A space $ X $ is said to be set star-Lindel\"{o}f (resp., set strongly star-Lindel\"{o}f) if for each nonempty subset $ A $ of $ X $ and each collection $ \mathcal{U} $ of open sets in $ X $ such that $ \overline{A} \subseteq \bigcup \mathcal{U} $, there is a countable subset $ \mathcal{V}$ of $ \mathcal{U} $ (resp., countable subset $ F $ of $ \overline{A} $) such that $ A \subseteq {\rm St}( \bigcup \mathcal{V}, \mathcal{U})$ (resp., $ A \subseteq {\rm St}( F, \mathcal{U})$). The classes of set star-Lindel\"{o}f spaces and set strongly star-Lindel\"{o}f spaces lie between the class of Lindel\"{o}f spaces and the class of star-Lindel\"{o}f spaces. In this paper, we investigate the relationship among set star-Lindel\"{o}f spaces, set strongly star-Lindel\"{o}f spaces, and other related spaces by providing some suitable examples and study the topological properties of set star-Lindel\"{o}f and set strongly star-Lindel\"{o}f spaces.


title: Jamming Distance Dictates Colloidal Shear Thickening

abstract: We report experimental and computational observations of dynamic contact networks for colloidal suspensions undergoing shear thickening. The dense suspensions are comprised of sterically stabilized poly(methyl methacrylate) hard sphere colloids that are spherically symmetric and have varied surface roughness. Confocal rheometry and dissipative particle dynamics simulations show that the shear thickening strength scales exponentially with the scaled deficit contact number and the scaled jamming distance. Rough colloids, which experience additional tangential and rolling constraints, require an average of 1.5 - 2 fewer particle contacts as compared to smooth colloids, in order to generate the same shear thickening strength. This is because the surface roughness enhances geometric friction in a way that the rough colloids do not experience a large change in the free volume near the jamming point. In contrast, smooth colloids must undergo significant reduction in the free volume to support an equivalent shear stress. The available free volume for different colloid roughness is related to the deficiency from the maximum number of nearest neighbors at jamming under shear. Our results further suggest that the force per contact is different for particles with different morphologies.


title: Diffractive $\gamma \gamma$ production in $pp$ collisions at the LHC

abstract: In this letter we estimate the contribution of the double diffractive processes for the diphoton production in $pp$ collisions at the Large Hadron Collider (LHC). The acceptance of the central and forward LHC detectors is taken into account and predictions for the invariant mass, rapidity and, transverse momentum distributions are presented. A comparison with the predictions for the Light -- by -- Light (LbL) scattering and exclusive diphoton production is performed. We demonstrate that the events associated to double diffractive processes can be separated and its study can be used to constrain the behavior of the diffractive parton distribution functions.


title: Forbidding $K_{2,t}$ traces in triple systems

abstract: Let $H$ and $F$ be hypergraphs. We say $H$ contains $F$ as a trace if there exists some set $S \subseteq V(H)$ such that $H|_S:=\{E\cap S: E \in E(H)\}$ contains a subhypergraph isomorphic to $F$. In this paper we give an upper bound on the number of edges in a $3$-uniform hypergraph that does not contain $K_{2,t}$ as a trace when $t$ is large. In particular, we show that $ \lim_{t\to \infty}\lim_{n\to \infty} \frac{\mathrm{ex}(n, \mathrm{Tr}_3(K_{2,t}))}{t^{3/2}n^{3/2}} = \frac{1}{6}.$ Moreover, we show $\frac{1}{2} n^{3/2} + o(n^{3/2}) \leq \mathrm{ex}(n, \mathrm{Tr}_3(C_4)) \leq \frac{5}{6} n^{3/2} + o(n^{3/2})$.


title: Gravitational Wave Test of the Strong Equivalence Principle

abstract: The Strong Equivalence Principle (SEP) holds the full essence and meaning of the General Theory of Relativity as the nonlinear relativistic theory of gravitation. It asserts the universal coupling of gravity to all matter and its interactions including the gravitational interaction and the gravitational self energy. We point out that confirming the gravitational coupling to gravitons, and hence to the gravitational waves, is the direct test of the SEP. We show that the near simultaneous detection of gravitational waves and gamma rays from the merger of binary neutron stars provides a unique and the most precise test of the SEP, better than a part in $10^{9}$, which is also the only test of the SEP in the radiation sector.


title: The geometric classification of nilpotent $\mathfrak{CD}$-algebras

abstract: We give a geometric classification of complex $4$-dimensional nilpotent $\mathfrak{CD}$-algebras. The corresponding geometric variety has dimension 18 and decomposes into $2$ irreducible components determined by the Zariski closures of a two-parameter family of algebras and a four-parameter family of algebras (see Theorem 2). In particular, there are no rigid $4$-dimensional complex nilpotent $\mathfrak{CD}$-algebras.


title: Fractional Covers of Hypergraphs with Bounded Multi-Intersection

abstract: Fractional (hyper-)graph theory is concerned with the specific problems that arise when fractional analogues of otherwise integer-valued (hyper-)graph invariants are considered. The focus of this paper is on fractional edge covers of hypergraphs. Our main technical result generalizes and unifies previous conditions under which the size of the support of fractional edge covers is bounded independently of the size of the hypergraph itself. This allows us to extend previous tractability results for checking if the fractional hypertree width of a given hypergraph is $\leq k$ for some constant $k$. We also show how our results translate to fractional vertex covers.


title: The Joint Automated Repository for Various Integrated Simulations (JARVIS) for data-driven materials design

abstract: The Joint Automated Repository for Various Integrated Simulations (JARVIS) is an integrated infrastructure to accelerate materials discovery and design using density functional theory (DFT), classical force-fields (FF), and machine learning (ML) techniques. JARVIS is motivated by the Materials Genome Initiative (MGI) principles of developing open-access databases and tools to reduce the cost and development time of materials discovery, optimization, and deployment. The major features of JARVIS are: JARVIS-DFT, JARVIS-FF, JARVIS-ML, and JARVIS-Tools. To date, JARVIS consists of 40,000 materials and 1 million calculated properties in JARVIS-DFT, 1,500 materials and 110 force-fields in JARVIS-FF, and 25 ML models for material-property predictions in JARVIS-ML, all of which are continuously expanding. JARVIS-Tools provides scripts and workflows for running and analyzing various simulations. We compare our computational data to experiments or high-fidelity computational methods wherever applicable to evaluate error/uncertainty in predictions. In addition to the existing workflows, the infrastructure can support a wide variety of other technologically important applications as part of the data-driven materials design paradigm. The JARVIS datasets and tools are publicly available at the website: https://jarvis.nist.gov .


title: A Dynamic Stability and Performance Analysis of Automatic Generation Control

abstract: Automatic generation control (AGC) is one of the most important coordinated control systems present in modern interconnected power systems. Despite being heavily studied, no interconnected dynamic stability and performance analysis of AGC is available in the literature. This paper presents such an analysis for a class of multi-area interconnected nonlinear power systems, providing a nonlinear stability proof and examining dynamic performance from the perspectives of non-interactive tuning, response speed, and sensitivity to disturbances. The key insight is that dynamic stability and performance can be rigorously captured by a reduced dynamic model, which depends only on the AGC controller parameters and on the area frequency characteristic constants. Our analytical results clarify some of the historical controversy concerning the tuning of frequency bias constants in AGC for dynamic stability and performance, and are validated by simulations on a detailed test system.


title: PsychFM: Predicting your next gamble

abstract: There is a sudden surge to model human behavior due to its vast and diverse applications which includes modeling public policies, economic behavior and consumer behavior. Most of the human behavior itself can be modeled into a choice prediction problem. Prospect theory is a theoretical model that tries to explain the anomalies in choice prediction. These theories perform well in terms of explaining the anomalies but they lack precision. Since the behavior is person dependent, there is a need to build a model that predicts choices on a per-person basis. Looking on at the average persons choice may not necessarily throw light on a particular person's choice. Modeling the gambling problem on a per person basis will help in recommendation systems and related areas. A novel hybrid model namely psychological factorisation machine ( PsychFM ) has been proposed that involves concepts from machine learning as well as psychological theories. It outperforms the popular existing models namely random forest and factorisation machines for the benchmark dataset CPC-18. Finally,the efficacy of the proposed hybrid model has been verified by comparing with the existing models.


title: Universal Distances for Extended Persistence

abstract: The extended persistence diagram is an invariant of piecewise linear functions, which is known to be stable under perturbations of functions with respect to the bottleneck distance as introduced by Cohen-Steiner, Edelsbrunner, and Harer. We address the question of universality, which asks for the largest possible stable distance on extended persistence diagrams, showing that a more discriminative variant of the bottleneck distance is universal. Our result applies more generally to settings where persistence diagrams are considered only up to a certain degree. We achieve our results by establishing a functorial construction and several characteristic properties of relative interlevel set homology, which mirror the classical Eilenberg--Steenrod axioms. Finally, we contrast the bottleneck distance with the interleaving distance of sheaves on the real line by showing that the latter is not intrinsic, let alone universal. This particular result has the further implication that the interleaving distance of Reeb graphs is not intrinsic either.


title: A scanning tunneling microscope capable of electron spin resonance and pump-probe spectroscopy at mK temperature and in vector magnetic field

abstract: In the last decade, detecting spin dynamics at the atomic scale has been enabled by combining techniques like electron spin resonance (ESR) or pump-probe spectroscopy with scanning tunneling microscopy (STM). Here, we demonstrate an ultra-high vacuum (UHV) STM operational at milliKelvin (mK) and in a vector magnetic field capable of both ESR and pump-probe spectroscopy. By implementing GHz compatible cabling, we achieve appreciable RF amplitudes at the junction while maintaining mK base temperature. We demonstrate the successful operation of our setup by utilizing two experimental ESR modes (frequency sweep and magnetic field sweep) on an individual TiH molecule on MgO/Ag(100) and extract the effective g-factor. We trace the ESR transitions down to MHz into an unprecedented low frequency band enabled by the mK base temperature. We also implement an all-electrical pump-probe scheme based on waveform sequencing suited for studying dynamics down to the nanoseconds range. We benchmark our system by detecting the spin relaxation time T1 of individual Fe atoms on MgO/Ag(100) and note a field strength and orientation dependent relaxation time.


title: Pretrained Semantic Speech Embeddings for End-to-End Spoken Language Understanding via Cross-Modal Teacher-Student Learning

abstract: Spoken language understanding is typically based on pipeline architectures including speech recognition and natural language understanding steps. These components are optimized independently to allow usage of available data, but the overall system suffers from error propagation. In this paper, we propose a novel training method that enables pretrained contextual embeddings to process acoustic features. In particular, we extend it with an encoder of pretrained speech recognition systems in order to construct end-to-end spoken language understanding systems. Our proposed method is based on the teacher-student framework across speech and text modalities that aligns the acoustic and the semantic latent spaces. Experimental results in three benchmarks show that our system reaches the performance comparable to the pipeline architecture without using any training data and outperforms it after fine-tuning with ten examples per class on two out of three benchmarks.


title: LOOC: Localize Overlapping Objects with Count Supervision

abstract: Acquiring count annotations generally requires less human effort than point-level and bounding box annotations. Thus, we propose the novel problem setup of localizing objects in dense scenes under this weaker supervision. We propose LOOC, a method to Localize Overlapping Objects with Count supervision. We train LOOC by alternating between two stages. In the first stage, LOOC learns to generate pseudo point-level annotations in a semi-supervised manner. In the second stage, LOOC uses a fully-supervised localization method that trains on these pseudo labels. The localization method is used to progressively improve the quality of the pseudo labels. We conducted experiments on popular counting datasets. For localization, LOOC achieves a strong new baseline in the novel problem setup where only count supervision is available. For counting, LOOC outperforms current state-of-the-art methods that only use count as their supervision. Code is available at: https://github.com/ElementAI/looc.


title: Minimal Modification to Nos\'e-Hoover Barostat Enables Correct NPT Sampling

abstract: The Nos\'e-Hoover dynamics for isothermal-isobaric (NPT) computer simulations do not generate the appropriate partition function for ergodic systems. The present paper points out that this can be corrected with a simple addition of a constant term to only one of the equations of motion. The solution proposed is much simpler than previous modifications done towards the same goal. The present modification is motivated by the work virial theorem, which has been derived for the special case of an infinitely periodic system in the first part of this paper.


title: Expected Eligibility Traces

abstract: The question of how to determine which states and actions are responsible for a certain outcome is known as the credit assignment problem and remains a central research question in reinforcement learning and artificial intelligence. Eligibility traces enable efficient credit assignment to the recent sequence of states and actions experienced by the agent, but not to counterfactual sequences that could also have led to the current state. In this work, we introduce expected eligibility traces. Expected traces allow, with a single update, to update states and actions that could have preceded the current state, even if they did not do so on this occasion. We discuss when expected traces provide benefits over classic (instantaneous) traces in temporal-difference learning, and show that sometimes substantial improvements can be attained. We provide a way to smoothly interpolate between instantaneous and expected traces by a mechanism similar to bootstrapping, which ensures that the resulting algorithm is a strict generalisation of TD($\lambda$). Finally, we discuss possible extensions and connections to related ideas, such as successor features.


title: Rejection criteria based on outliers in the KiDS photometric redshifts and PDF distributions derived by machine learning

abstract: The Probability Density Function (PDF) provides an estimate of the photometric redshift (zphot) prediction error. It is crucial for current and future sky surveys, characterized by strict requirements on the zphot precision, reliability and completeness. The present work stands on the assumption that properly defined rejection criteria, capable of identifying and rejecting potential outliers, can increase the precision of zphot estimates and of their cumulative PDF, without sacrificing much in terms of completeness of the sample. We provide a way to assess rejection through proper cuts on the shape descriptors of a PDF, such as the width and the height of the maximum PDF's peak. In this work we tested these rejection criteria to galaxies with photometry extracted from the Kilo Degree Survey (KiDS) ESO Data Release 4, proving that such approach could lead to significant improvements to the zphot quality: e.g., for the clipped sample showing the best trade-off between precision and completeness, we achieve a reduction in outliers fraction of $\simeq 75\%$ and an improvement of $\simeq 6\%$ for NMAD, with respect to the original data set, preserving the $\simeq 93\%$ of its content.


title: Galois representations on the cohomology of hyper-K\"{a}hler varieties

abstract: We show that the Andr\'{e} motive of a hyper-K\"{a}hler variety $X$ over a field $K \subset \mathbb{C}$ with $b_2(X)>6$ is governed by its component in degree $2$. More precisely, we prove that if $X_1$ and $X_2$ are deformation equivalent hyper-K\"{a}hler varieties with $b_2(X_i)>6$ and if there exists a Hodge isometry $f\colon H^2(X_1,\mathbb{Q})\to H^2(X_2,\mathbb{Q})$, then the Andr\'e motives of $X_1$ and $X_2$ are isomorphic after a finite extension of $K$, up to an additional technical assumption in presence of non-trivial odd cohomology. As a consequence, the Galois representations on the \'{e}tale cohomology of $X_1$ and $X_2$ are isomorphic as well. We prove a similar result for varieties over a finite field which can be lifted to hyper-K\"{a}hler varieties for which the Mumford--Tate conjecture is true.


title: Incidence Hypergraphs: Box Products & the Laplacian

abstract: The box product and its associated box exponential are characterized for the categories of quivers (directed graphs), multigraphs, set system hypergraphs, and incidence hypergraphs. It is shown that only the quiver case of the box exponential can be characterized via homs entirely within their own category. An asymmetry in the incidence hypergraphic box product is rectified via an incidence dual-closed generalization that effectively treats vertices and edges as real and imaginary parts of a complex number, respectively. This new hypergraphic box product is shown to have a natural interpretation as the canonical box product for graphs via the bipartite representation functor, and its associated box exponential is represented as homs entirely in the category of incidence hypergraphs; with incidences determined by incidence-prism mapping. The evaluation of the box exponential at paths is shown to correspond to the entries in half-powers of the oriented hypergraphic signless Laplacian matrix.


title: Sharp discontinuous traveling waves in a hyperbolic Keller--Segel equation

abstract: In this work we describe a hyperbolic model with cell-cell repulsion with a dynamics in the population of cells. More precisely, we consider a population of cells producing a field (which we call "pressure") which induces a motion of the cells following the opposite of the gradient. The field indicates the local density of population and we assume that cells try to avoid crowded areas and prefer locally empty spaces which are far away from the carrying capacity. We analyze the well-posedness property of the associated Cauchy problem on the real line. We start from bounded initial conditions and we consider some invariant properties of the initial conditions such as the continuity, smoothness and monotonicity. We also describe in detail the behavior of the level sets near the propagating boundary of the solution and we find that an asymptotic jump is formed on the solution for a natural class of initial conditions. Finally, we prove the existence of sharp traveling waves for this model, which are particular solutions traveling at a constant speed, and argue that sharp traveling waves are necessarily discontinuous. This analysis is confirmed by numerical simulations of the PDE problem.


title: KiDS-1000 Methodology: Modelling and inference for joint weak gravitational lensing and spectroscopic galaxy clustering analysis

abstract: We present the methodology for a joint cosmological analysis of weak gravitational lensing from the fourth data release of the ESO Kilo-Degree Survey (KiDS-1000) and galaxy clustering from the partially overlapping BOSS and 2dFLenS surveys. Cross-correlations between galaxy positions and ellipticities have been incorporated into the analysis, necessitating a hybrid model of non-linear scales that blends perturbative and non-perturbative approaches, and an assessment of contributions by astrophysical effects. All weak lensing signals are measured consistently via Fourier-space statistics that are insensitive to the survey mask and display low levels of mode mixing. The calibration of photometric redshift distributions and multiplicative gravitational shear bias has been updated, and a more complete tally of residual calibration uncertainties is propagated into the likelihood. A dedicated suite of more than 20000 mocks is used to assess the performance of covariance models and to quantify the impact of survey geometry and spatial variations of survey depth on signals and their errors. The sampling distributions for the likelihood and the $\chi^2$ goodness-of-fit statistic have been validated, with proposed changes to the number of degrees of freedom. Standard weak lensing point estimates on $S_8=\sigma_8\,(\Omega_{\rm m}/0.3)^{1/2}$ derived from its marginal posterior are easily misinterpreted to be biased low, and an alternative estimator and associated credible interval have been proposed. Known systematic effects pertaining to weak lensing modelling and inference are shown to bias $S_8$ by no more than 0.1 standard deviations, with the caveat that no conclusive validation data exist for models of intrinsic galaxy alignments. Compared to the previous KiDS analyses, $S_8$ constraints are expected to improve by 20% for weak lensing alone and by 29% for the joint analysis. [abridged]


title: KiDS-1000 catalogue: Weak gravitational lensing shear measurements

abstract: We present weak lensing shear catalogues from the fourth data release of the Kilo-Degree Survey, KiDS-1000, spanning 1006 square degrees of deep and high-resolution imaging. Our `gold-sample' of galaxies, with well-calibrated photometric redshift distributions, consists of 21 million galaxies with an effective number density of $6.17$ galaxies per square arcminute. We quantify the accuracy of the spatial, temporal, and flux-dependent point-spread function (PSF) model, verifying that the model meets our requirements to induce less than a $0.1\sigma$ change in the inferred cosmic shear constraints on the clustering cosmological parameter $S_8 = \sigma_8\sqrt{\Omega_{\rm m}/0.3}$. Through a series of two-point null-tests, we validate the shear estimates, finding no evidence for significant non-lensing B-mode distortions in the data. The PSF residuals are detected in the highest-redshift bins, originating from object selection and/or weight bias. The amplitude is, however, shown to be sufficiently low and within our stringent requirements. With a shear-ratio null-test, we verify the expected redshift scaling of the galaxy-galaxy lensing signal around luminous red galaxies. We conclude that the joint KiDS-1000 shear and photometric redshift calibration is sufficiently robust for combined-probe gravitational lensing and spectroscopic clustering analyses.


title: Testing KiDS cross-correlation redshifts with simulations

abstract: Measuring cosmic shear in wide-field imaging surveys requires accurate knowledge of the redshift distribution of all sources. The clustering-redshift technique exploits the angular cross-correlation of a target galaxy sample with unknown redshifts and a reference sample with known redshifts, and is an attractive alternative to colour-based methods of redshift calibration. We test the performance of such clustering redshift measurements using mock catalogues that resemble the Kilo-Degree Survey (KiDS). These mocks are created from the MICE simulation and closely mimic the properties of the KiDS source sample and the overlapping spectroscopic reference samples. We quantify the performance of the clustering redshifts by comparing the cross-correlation results with the true redshift distributions in each of the five KiDS photometric redshift bins. Such a comparison to an informative model is necessary due to the incompleteness of the reference samples at high redshifts. Clustering mean redshifts are unbiased at $|\Delta z|<0.006$ under these conditions. The redshift evolution of the galaxy bias can be reliably mitigated at this level of precision using auto-correlation measurements and self-consistency relations, and will not become a dominant source of systematic error until the arrival of Stage-IV cosmic shear surveys. Using redshift distributions from a direct colour-based estimate instead of the true redshift distributions as a model for comparison with the clustering redshifts increases the biases in the mean to up to $|\Delta z|\sim0.04$. This indicates that the interpretation of clustering redshifts in real-world applications will require more sophisticated (parameterised) models of the redshift distribution in the future. If such better models are available, the clustering-redshift technique promises to be a highly complementary alternative to other methods of redshift calibration.


title: Causality in Curved Spacetimes: The Speed of Light & Gravity

abstract: Within the low-energy effective field theories of QED and gravity, the low-energy speed of light or that of gravitational waves can typically be mildly superluminal in curved spacetimes. Related to this, small scattering time advances relative to the curved background can emerge from known effective field theory coefficients for photons or gravitons. We clarify why these results are not in contradiction with causality, analyticity or Lorentz invariance, and highlight various subtleties that arise when dealing with superluminalities and time advances in the gravitational context. Consistent low-energy effective theories are shown to self-protect by ensuring that any time advance and superluminality calculated within the regime of validity of the effective theory is necessarily unresolvable, and cannot be argued to lead to a macroscopically larger lightcone. Such considerations are particularly relevant for putting constraints on cosmological and gravitational effective field theories and we provide explicit criteria to be satisfied so as to ensure causality.


title: A double $(\infty,1)$-categorical nerve for double categories

abstract: We construct a nerve from double categories into double $(\infty,1)$-categories and show that it gives a right Quillen and homotopically fully faithful functor between the model structure for weakly horizontally invariant double categories and the model structure on bisimplicial spaces for double $(\infty,1)$-categories seen as double Segal objects in spaces complete in the horizontal direction. We then restrict the nerve along a homotopical horizontal embedding of 2-categories into double categories, and show that it gives a right Quillen and homotopically fully faithful functor between Lack's model structure for 2-categories and the model structure for 2-fold complete Segal spaces. We further show that Lack's model structure is right-induced along this nerve from the model structure for 2-fold complete Segal spaces.


title: Effective birationality for sub-pairs with real coefficients

abstract: For $\epsilon$-lc Fano type varieties $X$ of dimension $d$ and a given finite set $\Gamma$, we show that there exists a positive integer $m_0$ which only depends on $\epsilon,d$ and $\Gamma$, such that both $|-mK_X-\sum_i\lceil mb_i\rceil B_i|$ and $|-mK_X-\sum_i\lfloor mb_i\rfloor B_i|$ define birational maps for any $m\ge m_0$ provided that $B_i$ are pseudo-effective Weil divisors, $b_i\in\Gamma$, and $-(K_X+\sum_ib_iB_i)$ is big. When $\Gamma\subset[0,1]$ satisfies the DCC but is not finite, we construct an example to show that the effective birationality may fail even if $X$ is fixed, $B_i$ are fixed prime divisors, and $(X,B)$ is $\epsilon'$-lc for some $\epsilon'>0$.


title: Variational Autoencoders for Anomalous Jet Tagging

abstract: We present a detailed study on Variational Autoencoders (VAEs) for anomalous jet tagging at the Large Hadron Collider. By taking in low-level jet constituents' information, and training with background QCD jets in an unsupervised manner, the VAE is able to encode important information for reconstructing jets, while learning an expressive posterior distribution in the latent space. When using the VAE as an anomaly detector, we present different approaches to detect anomalies: directly comparing in the input space or, instead, working in the latent space. In order to facilitate general search approaches such as bump-hunt, mass-decorrelated VAEs based on distance correlation regularization are also studied. We find that the naive mass-decorrelated VAEs fail at maintaining proper detection performance, by assigning higher probabilities to some anomalous samples. To build a performant mass-decorrelated anomalous jet tagger, we propose the Outlier Exposed VAE (OE-VAE), for which some outlier samples are introduced in the training process to guide the learned information. OE-VAEs are employed to achieve two goals at the same time: increasing sensitivity of outlier detection and decorrelating jet mass from the anomaly score. We succeed in reaching excellent results from both aspects. Code implementation of this work can be found at https://github.com/taolicheng/VAE-Jet


title: Swoosh! Rattle! Thump! -- Actions that Sound

abstract: Truly intelligent agents need to capture the interplay of all their senses to build a rich physical understanding of their world. In robotics, we have seen tremendous progress in using visual and tactile perception; however, we have often ignored a key sense: sound. This is primarily due to the lack of data that captures the interplay of action and sound. In this work, we perform the first large-scale study of the interactions between sound and robotic action. To do this, we create the largest available sound-action-vision dataset with 15,000 interactions on 60 objects using our robotic platform Tilt-Bot. By tilting objects and allowing them to crash into the walls of a robotic tray, we collect rich four-channel audio information. Using this data, we explore the synergies between sound and action and present three key insights. First, sound is indicative of fine-grained object class information, e.g., sound can differentiate a metal screwdriver from a metal wrench. Second, sound also contains information about the causal effects of an action, i.e. given the sound produced, we can predict what action was applied to the object. Finally, object representations derived from audio embeddings are indicative of implicit physical properties. We demonstrate that on previously unseen objects, audio embeddings generated through interactions can predict forward models 24% better than passive visual embeddings. Project videos and data are at https://dhiraj100892.github.io/swoosh/


title: Language-agnostic BERT Sentence Embedding

abstract: While BERT is an effective method for learning monolingual sentence embeddings for semantic similarity and embedding based transfer learning (Reimers and Gurevych, 2019), BERT based cross-lingual sentence embeddings have yet to be explored. We systematically investigate methods for learning multilingual sentence embeddings by combining the best methods for learning monolingual and cross-lingual representations including: masked language modeling (MLM), translation language modeling (TLM) (Conneau and Lample, 2019), dual encoder translation ranking (Guo et al., 2018), and additive margin softmax (Yang et al., 2019a). We show that introducing a pre-trained multilingual language model dramatically reduces the amount of parallel training data required to achieve good performance by 80%. Composing the best of these methods produces a model that achieves 83.7% bi-text retrieval accuracy over 112 languages on Tatoeba, well above the 65.5% achieved by Artetxe and Schwenk (2019b), while still performing competitively on monolingual transfer learning benchmarks (Conneau and Kiela, 2018). Parallel data mined from CommonCrawl using our best model is shown to train competitive NMT models for en-zh and en-de. We publicly release our best multilingual sentence embedding model for 109+ languages at https://tfhub.dev/google/LaBSE.


title: From Weak Antilocalization to Kondo Scattering in a Magnetic Complex Oxide Interface

abstract: Quantum corrections to electrical resistance can serve as sensitive probes of the magnetic landscape of a material. For example, interference between time-reversed electron paths gives rise to weak localization effects, which can provide information about the coupling between spins and orbital motion, while the Kondo effect is sensitive to the presence of spin impurities. Here we use low-temperature magnetotransport measurements to reveal a transition from weak antilocalization (WAL) to Kondo scattering in the quasi-two-dimensional electron gas formed at the interface between SrTiO$_3$ and the Mott insulator NdTiO$_3$. This transition occurs as the thickness of the NdTiO$_3$ layer is increased. Analysis of the Kondo scattering and WAL points to the presence of atomic-scale magnetic impurities coexisting with extended magnetic regions that affect transport via a strong magnetic exchange interaction. This leads to distinct magnetoresistance behaviors that can serve as a sensitive probe of magnetic properties in two dimensions.


title: Bernoulli Randomness and Biased Normality

abstract: One can consider $\mu$-Martin-L\"of randomness for a probability measure $\mu$ on $2^{\omega}$, such as the Bernoulli measure $\mu_p$ given $p \in (0, 1)$. We study Bernoulli randomness of sequences in $n^{\omega}$ with parameters $p_0, p_1, \dotsc, p_{n-1}$, and we introduce a biased version of normality. We prove that every Bernoulli random real is normal in the biased sense, and this has the corollary that the set of biased normal reals has full Bernoulli measure in $n^{\omega}$. We give an algorithm for computing biased normal sequences from normal sequences, so that we can give explicit examples of biased normal reals. We investigate an application of randomness to iterated function systems. Finally, we list a few further questions relating to Bernoulli randomness and biased normality.


title: Trace-Norm Adversarial Examples

abstract: White box adversarial perturbations are sought via iterative optimization algorithms most often minimizing an adversarial loss on a $l_p$ neighborhood of the original image, the so-called distortion set. Constraining the adversarial search with different norms results in disparately structured adversarial examples. Here we explore several distortion sets with structure-enhancing algorithms. These new structures for adversarial examples, yet pervasive in optimization, are for instance a challenge for adversarial theoretical certification which again provides only $l_p$ certificates. Because adversarial robustness is still an empirical field, defense mechanisms should also reasonably be evaluated against differently structured attacks. Besides, these structured adversarial perturbations may allow for larger distortions size than their $l_p$ counter-part while remaining imperceptible or perceptible as natural slight distortions of the image. Finally, they allow some control on the generation of the adversarial perturbation, like (localized) bluriness.


title: $^{26}$Al-$^{26}$Mg isotopic, mineralogy, petrography of a Hibonite-Pyroxene Spherule in Allan Hills 77307 (CO3.03): Implications for the origin and evolution of these objects

abstract: 10 Hibonite-pyroxene/glass spherules discovered hitherto are a rare suite of refractory inclusions that show the largest range of exotic isotopic properties (anomalies in neutron rich isotopes (e.g., $^{48}$Ca, $^{50}$Ti), abundance of $^{26}$Al) despite their defining simple spherical morphology and mineralogy consisting predominantly of few hibonites nestled within/with glassy or crystallised calcium, aluminium-rich pyroxene. $^{26}$Al-$^{26}$Mg chronological studies along with petrography and mineralogy of a relatively large (~120 micron diameter), found in Allan Hills 77307 (CO3.03) has been performed. Uniquely, both hibonite and pyroxene show discordant abundance of short-lived now-extinct radionuclide $^{26}$Al that suggest disparate and distinct regions of origin of hibonite and pyroxene. The pristine petrography and mineralogy of this inclusion allow discernment of their genesis and trend of alteration in hibonite-pyroxene/glass spherules.


title: Deep Learning Models for Visual Inspection on Automotive Assembling Line

abstract: Automotive manufacturing assembly tasks are built upon visual inspections such as scratch identification on machined surfaces, part identification and selection, etc, which guarantee product and process quality. These tasks can be related to more than one type of vehicle that is produced within the same manufacturing line. Visual inspection was essentially human-led but has recently been supplemented by the artificial perception provided by computer vision systems (CVSs). Despite their relevance, the accuracy of CVSs varies accordingly to environmental settings such as lighting, enclosure and quality of image acquisition. These issues entail costly solutions and override part of the benefits introduced by computer vision systems, mainly when it interferes with the operating cycle time of the factory. In this sense, this paper proposes the use of deep learning-based methodologies to assist in visual inspection tasks while leaving very little footprints in the manufacturing environment and exploring it as an end-to-end tool to ease CVSs setup. The proposed approach is illustrated by four proofs of concept in a real automotive assembly line based on models for object detection, semantic segmentation, and anomaly detection.


title: On Cauchy dual operator and duality for Banach spaces of analytic functions

abstract: In this paper, two related types of dualities are investigated. The first is the duality between left-invertible operators and the second is the duality between Banach spaces of vector-valued analytic functions. We will examine a pair ($\mathcal{B},\Psi)$ consisting of a reflexive Banach spaces $\mathcal{B}$ of vector-valued analytic functions on which a left-invertible multiplication operator acts and an operator-valued holomorphic function $\Psi$. We prove that there exist a dual pair ($\mathcal{B}^\prime,\Psi^\prime)$ such that the space $\mathcal{B}^\prime$ is unitarily equivalent to the space $\mathcal{B}^*$ and the following intertwining relations hold \begin{equation*} \mathscr{L} \mathcal{U} = \mathcal{U}\mathscr{M}_z^* \quad\text{and}\quad \mathscr{M}_z\mathcal{U} = \mathcal{U} \mathscr{L}^*, \end{equation*} where $\mathcal{U}$ is the unitary operator between $\mathcal{B}^\prime$ and $\mathcal{B}^*$. In addition we show that $\Psi$ and $\Psi^\prime$ are connected through the relation\begin{equation*} \langle(\Psi^\prime( \bar{z}) e_1) (\lambda),e_2\rangle= \langle e_1,(\Psi( \bar{ \lambda}) e_2)(z)\rangle \end{equation*} for every $e_1,e_2\in E$, $z\in \varOmega$, $\lambda\in \varOmega^\prime$. If a left-invertible operator $T$ satisfies certain conditions, then both $T$ and the Cauchy dual operator $T^\prime$ can be modelled as a multiplication operator on reproducing kernel Hilbert spaces of vector-valued analytic functions $\mathscr{H}$ and $\mathscr{H}^\prime$, respectively. We prove that Hilbert space of the dual pair of $(\mathscr{H},\Psi)$ coincide with $\mathscr{H}^\prime$, where $\Psi$ is a certain operator-valued holomorphic function. Moreover, we characterize when the duality between spaces $\mathscr{H}$ and $\mathscr{H}^\prime$ obtained by identifying them with $\mathcal{H}$ is the same as the duality obtained from the Cauchy pairing.


title: Improved flat-back 3D gadgets in origami extrusions completely downward compatible with the conventional pyramid-supported 3D gadgets

abstract: An origami extrusion is a folding of a 3D object in the middle of a flat piece of paper, using 3D gadgets which create faces with solid angles. In this paper we focus on 3D gadgets which create a top face parallel to the ambient paper and two side faces sharing a ridge, with two outgoing simple pleats, where a simple pleat is a pair of a mountain fold and a valley fold. There are two such types of 3D gadgets. One is the conventional type of 3D gadgets with a triangular pyramid supporting the two side faces from inside. The other is the newer type of 3D gadgets presented in our previous paper, which improve the conventional ones in several respects: They have flat back sides above the ambient paper and no gap between the side faces; they are less interfering with adjacent gadgets so that we can make the extrusion higher at one time; they are downward compatible with conventional ones if constructible; they have a modified flat-back gadget used for repetition which does not interfere with adjacent gadgets; the angles of their outgoing pleats can be changed under certain conditions. However, there are cases where we can apply the conventional gadgets while we cannot our previous ones. The purpose of this paper is to improve our previous 3D gadgets to be completely downward compatible with the conventional ones, in the sense that any conventional gadget can be replaced by our improved one with the same outgoing pleats, but the converse is not always possible. To be more precise, we prove that for any given conventional 3D gadget there are an infinite number of improved 3D gadgets which are compatible with it, and the conventional 3D gadget can be replaced with any of these 3D gadgets without affecting any other conventional 3D gadget. Also, we see that our improved 3D gadget keep all of the above advantages over the conventional ones.


title: Sensitivity Analysis for the 2D Navier-Stokes Equations with Applications to Continuous Data Assimilation

abstract: We rigorously prove the well-posedness of the formal sensitivity equations with respect to the Reynolds number corresponding to the 2D incompressible Navier-Stokes equations. Moreover, we do so by showing a sequence of difference quotients converges to the unique solution of the sensitivity equations for both the 2D Navier-Stokes equations and the related data assimilation equations, which utilize the continuous data assimilation algorithm proposed by Azouani, Olson, and Titi. As a result, this method of proof provides uniform bounds on difference quotients, demonstrating parameter recovery algorithms that change parameters as the system evolves will not blow-up. We also note that this appears to be the first such rigorous proof of global existence and uniqueness to strong or weak solutions to the sensitivity equations for the 2D Navier-Stokes equations (in the natural case of zero initial data), and that they can be obtained as a limit of difference quotients with respect to the Reynolds number.


title: IPO: Iterative Physical Optics Image Approximation

abstract: An improved Iterative Physical Optics (IPO) image approximation method has been presented to dramatically increase the accuracy of the approximation and extend its applicability to PEC surfaces with smaller radii or larger curvatures. Starting from the first-order conventional PO image approximation, the IPO image approximation method iteratively correct the surface current to compensate the deviation of the electric field boundary condition on the PEC surfaces, making use of the local plane wave approximation. Numerical validations with two popular PEC surfaces, i.e., the parabolic dish antennas and the PEC spheres, are carried out and the results show that the IPO approximation method increases the surface current accuracy by more than two orders of magnitude, compared to the conventional PO image approximation method.


title: An Analysis of Data Driven, Decision-Making Capabilities of Managers in Banks

abstract: Organizations are adopting data analytics and Business Intelligence (BI) tools to gain insights from the past data, forecast future events, and to get timely and reliable information for decision making. While the tools are becoming mature, affordable, and more comfortable to use, it is also essential to understand whether the contemporary managers and leaders are ready for Data-Driven Decision Making (DDDM). We explore the extent the Decision Makers (DMs) utilize data and tools, as well as their ability to interpret various forms of outputs from tools and to apply those insights to gain competitive advantage. Our methodology was based on a qualitative survey, where we interviewed 12 DMs of six commercial banks in Sri Lanka at the branch, regional, and CTO, CIO, and Head of IT levels. We identified that on many occasions, DMs' intuition overrules the DDDM due to uncertainty, lack of trust, knowledge, and risk-taking. Moreover, it was identified that the quality of visualizations has a significant impact on the use of intuition by overruling DDDM. We further provide a set of recommendations on the adoption of BI tools and how to overcome the struggles faced while performing DDDM.


title: An optimal design problem with non-standard growth and no concentration effects

abstract: We obtain an integral representation for certain functionals arising in the context of optimal design and damage evolution problems under non-standard growth conditions and perimeter penalisation. Under our hypotheses, the integral representation includes a term which is absolutely continuous with respect to the Lebesgue measure and a perimeter term, but no additional singular term. We also study some dimension reduction problems providing results for the optimal design of thin films.


title: Accurate Bounding-box Regression with Distance-IoU Loss for Visual Tracking

abstract: Most existing trackers are based on using a classifier and multi-scale estimation to estimate the target state. Consequently, and as expected, trackers have become more stable while tracking accuracy has stagnated. While trackers adopt a maximum overlap method based on an intersection-over-union (IoU) loss to mitigate this problem, there are defects in the IoU loss itself, that make it impossible to continue to optimize the objective function when a given bounding box is completely contained within/without another bounding box; this makes it very challenging to accurately estimate the target state. Accordingly, in this paper, we address the above-mentioned problem by proposing a novel tracking method based on a distance-IoU (DIoU) loss, such that the proposed tracker consists of target estimation and target classification. The target estimation part is trained to predict the DIoU score between the target ground-truth bounding-box and the estimated bounding-box. The DIoU loss can maintain the advantage provided by the IoU loss while minimizing the distance between the center points of two bounding boxes, thereby making the target estimation more accurate. Moreover, we introduce a classification part that is trained online and optimized with a Conjugate-Gradient-based strategy to guarantee real-time tracking speed. Comprehensive experimental results demonstrate that the proposed method achieves competitive tracking accuracy when compared to state-of-the-art trackers while with a real-time tracking speed.


title: A New Linear Inversion Formula for a class of Hypergeometric polynomials

abstract: Given complex parameters $x$, $\nu$, $\alpha$, $\beta$ and $\gamma \notin -\mathbb{N}$, consider the infinite lower triangular matrix $\mathbf{A}(x,\nu;\alpha, \beta,\gamma)$ with elements $$ A_{n,k}(x,\nu;\alpha,\beta,\gamma) = \displaystyle (-1)^k\binom{n+\alpha}{k+\alpha} \cdot F(k-n,-(\beta+n)\nu;-(\gamma+n);x) $$ for $1 \leqslant k \leqslant n$, depending on the Hypergeometric polynomials $F(-n,\cdot;\cdot;x)$, $n \in \mathbb{N}^*$. After stating a general criterion for the inversion of infinite matrices in terms of associated generating functions, we prove that the inverse matrix $\mathbf{B}(x,\nu;\alpha, \beta,\gamma) = \mathbf{A}(x,\nu;\alpha, \beta,\gamma)^{-1}$ is given by \begin{align} B_{n,k}(x,\nu;\alpha, \beta,\gamma) = & \; \displaystyle (-1)^k\binom{n+\alpha}{k+\alpha} \; \cdot \nonumber \\ & \; \biggl [ \; \frac{\gamma+k}{\beta+k} \, F(k-n,(\beta+k)\nu;\gamma+k;x) \; + \nonumber \\ & \; \; \; \frac{\beta-\gamma}{\beta+k} \, F(k-n,(\beta+k)\nu;1+\gamma+k;x) \; \biggr ] \nonumber \end{align} for $1 \leqslant k \leqslant n$, thus providing a new class of linear inversion formulas. Functional relations for the generating functions of related sequences $S$ and $T$, that is, $T = \mathbf{A}(x,\nu;\alpha, \beta,\gamma) \, S \Longleftrightarrow S = \mathbf{B}(x,\nu;\alpha, \beta,\gamma) \, T$, are also provided.


title: Selecting Regions of Interest in Large Multi-Scale Images for Cancer Pathology

abstract: Recent breakthroughs in object detection and image classification using Convolutional Neural Networks (CNNs) are revolutionizing the state of the art in medical imaging, and microscopy in particular presents abundant opportunities for computer vision algorithms to assist medical professionals in diagnosis of diseases ranging from malaria to cancer. High resolution scans of microscopy slides called Whole Slide Images (WSIs) offer enough information for a cancer pathologist to come to a conclusion regarding cancer presence, subtype, and severity based on measurements of features within the slide image at multiple scales and resolutions. WSIs' extremely high resolutions and feature scales ranging from gross anatomical structures down to cell nuclei preclude the use of standard CNN models for object detection and classification, which have typically been designed for images with dimensions in the hundreds of pixels and with objects on the order of the size of the image itself. We explore parallel approaches based on Reinforcement Learning and Beam Search to learn to progressively zoom into the WSI to detect Regions of Interest (ROIs) in liver pathology slides containing one of two types of liver cancer, namely Hepatocellular Carcinoma (HCC) and Cholangiocarcinoma (CC). These ROIs can then be presented directly to the pathologist to aid in measurement and diagnosis or be used for automated classification of tumor subtype.


title: TLIO: Tight Learned Inertial Odometry

abstract: In this work we propose a tightly-coupled Extended Kalman Filter framework for IMU-only state estimation. Strap-down IMU measurements provide relative state estimates based on IMU kinematic motion model. However the integration of measurements is sensitive to sensor bias and noise, causing significant drift within seconds. Recent research by Yan et al. (RoNIN) and Chen et al. (IONet) showed the capability of using trained neural networks to obtain accurate 2D displacement estimates from segments of IMU data and obtained good position estimates from concatenating them. This paper demonstrates a network that regresses 3D displacement estimates and its uncertainty, giving us the ability to tightly fuse the relative state measurement into a stochastic cloning EKF to solve for pose, velocity and sensor biases. We show that our network, trained with pedestrian data from a headset, can produce statistically consistent measurement and uncertainty to be used as the update step in the filter, and the tightly-coupled system outperforms velocity integration approaches in position estimates, and AHRS attitude filter in orientation estimates.


title: Dalek -- a deep-learning emulator for TARDIS

abstract: Supernova spectral time series contain a wealth of information about the progenitor and explosion process of these energetic events. The modeling of these data requires the exploration of very high dimensional posterior probabilities with expensive radiative transfer codes. Even modest parametrizations of supernovae contain more than ten parameters and a detailed exploration demands at least several million function evaluations. Physically realistic models require at least tens of CPU minutes per evaluation putting a detailed reconstruction of the explosion out of reach of traditional methodology. The advent of widely available libraries for the training of neural networks combined with their ability to approximate almost arbitrary functions with high precision allows for a new approach to this problem. Instead of evaluating the radiative transfer model itself, one can build a neural network proxy trained on the simulations but evaluating orders of magnitude faster. Such a framework is called an emulator or surrogate model. In this work, we present an emulator for the TARDIS supernova radiative transfer code applied to Type Ia supernova spectra. We show that we can train an emulator for this problem given a modest training set of a hundred thousand spectra (easily calculable on modern supercomputers). The results show an accuracy on the percent level (that are dominated by the Monte Carlo nature of TARDIS and not the emulator) with a speedup of several orders of magnitude. This method has a much broader set of applications and is not limited to the presented problem.


title: New Recipes for Brownian Loop Soups

abstract: We define a large new class of conformal primary operators in the ensemble of Brownian loops in two dimensions known as the ``Brownian loop soup,'' and compute their correlation functions analytically and in closed form. The loop soup is a conformally invariant statistical ensemble with central charge $c = 2 \lambda$, where $\lambda > 0$ is the intensity of the soup. Previous work identified exponentials of the layering operator $e^{i \beta N(z)}$ as primary operators. Each Brownian loop was assigned $\pm 1$ randomly, and $N(z)$ was defined to be the sum of these numbers over all loops that encircle the point $z$. These exponential operators then have conformal dimension ${\frac{\lambda}{10}}(1 - \cos \beta)$. Here we generalize this procedure by assigning a more general random value to each loop. The operator $e^{i \beta N(z)}$ remains primary with conformal dimension $\frac {\lambda}{10}(1 - \phi(\beta))$, where $\phi(\beta)$ is the characteristic function of the probability distribution used to assign random values to each loop. Using recent results we compute in closed form the exact two-point functions in the upper half-plane and four-point functions in the full plane of this very general class of operators. These correlation functions depend analytically on the parameters $\lambda, \beta_i, z_i$, and on the characteristic function $\phi(\beta)$. They satisfy the conformal Ward identities and are crossing symmetric. As in previous work, the conformal block expansion of the four-point function reveals the existence of additional and as-yet uncharacterized conformal primary operators.

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


title: NGC 474 as viewed with KCWI: diagnosing a shell galaxy

abstract: We present new spectra obtained using Keck/KCWI and perform kinematics and stellar population analyses of the shell galaxy NGC 474, from both the galaxy centre and a region from the outer shell. We show that both regions have similarly extended star formation histories although with different stellar population properties. The central region of NGC 474 is dominated by intermediate-aged stars (8.3 \pm 0.3 Gyr) with subsolar metallicity ([Z/H]= -0.24 \pm 0.07 dex) while the observed shell region, which hosts a substantial population of younger stars, has a mean luminosity-weighted age of 4.0 \pm 0.5 Gyr with solar metallicities ([Z/H]=-0.03 \pm 0.09 dex). Our results are consistent with a scenario in which NGC 474 experienced a major to intermediate merger with a log\((M_*/M_\odot)\sim10 \) mass satellite galaxy at least \sim 2 Gyr ago which produced its shell system. This work shows that the direct spectroscopic study of low-surface brightness stellar features, such as shells, is now feasible and opens up a new window to understanding galaxy formation and evolution.


title: Evidence for a high mutual inclination between the cold Jupiter and transiting super Earth orbiting $\pi$ Men

abstract: $\pi$ Men hosts a transiting super Earth ($P\approx6.27$ d, $m\approx4.82$ $M_{\oplus}$, $R\approx2.04$ $R_{\oplus}$) discovered by TESS and a cold Jupiter ($P\approx2093$ d, $m \sin I\approx10.02$ $M_{\rm{Jup}}$, $e\approx0.64$) discovered from radial velocity. We use Gaia DR2 and Hipparcos astrometry to derive the star's velocity caused by the orbiting planets and constrain the cold Jupiter's sky-projected inclination ($I_b=41-65^{\circ}$). From this we derive the mutual inclination ($\Delta I$) between the two planets, and find that $49^{\circ}< \Delta I < 131^{\circ}$ (1$\sigma$), and $28^{\circ} < \Delta I < 152^{\circ}$ (2$\sigma$). We examine the dynamics of the system using $N$-body simulations, and find that potentially large oscillations in the super Earth's eccentricity and inclination are suppressed by general relativistic precession. However, nodal precession of the inner orbit around the invariable plane causes the super Earth to only transit between 7-22 per cent of the time, and to usually be observed as misaligned with the stellar spin axis. We repeat our analysis for HAT-P-11, finding a large $\Delta I$ between its close-in Neptune and cold Jupiter and similar dynamics. $\pi$ Men and HAT-P-11 are prime examples of systems where dynamically hot outer planets excite their inner planets, with the effects of increasing planet eccentricities, planet-star misalignments, and potentially reducing the transit multiplicity. Formation of such systems likely involves scattering between multiple giant planets or misaligned protoplanetary discs. Future imaging of the faint debris disc in $\pi$ Men and precise constraints on the stellar spin orientation would provide strong tests for these formation scenarios.


title: Discovery of diffuse optical emission lines from the inner Galaxy: Evidence for LI(N)ER-like gas

abstract: Optical emission lines are used to categorize galaxies into three groups according to their dominant central radiation source: active galactic nuclei, star formation, or low-ionization (nuclear) emission regions [LI(N)ERs] that may trace ionizing radiation from older stellar populations. Using the Wisconsin H-Alpha Mapper, we detect optical line emission in low-extinction windows within eight degrees of Galactic Center. The emission is associated with the 1.5-kiloparsec-radius "Tilted Disk" of neutral gas. We modify a model of this disk and find that the hydrogen gas observed is at least $48\%$ ionized. The ratio [NII] $\lambda$6584 $\overset{\lower.5em\circ}{\mathrm{A}}$/H$\alpha$ $\lambda$6563 $\overset{\lower.5em\circ}{\mathrm{A}}$ increases from 0.3 to 2.5 with Galactocentric radius; [OIII] $\lambda$5007 $\overset{\lower.5em\circ}{\mathrm{A}}$ and H$\beta$ $\lambda$4861 $\overset{\lower.5em\circ}{\mathrm{A}}$ are also sometimes detected. The line ratios for most Tilted Disk sightlines are characteristic of LI(N)ER galaxies.


title: Constraints on Axion-Lepton coupling from Big Bang Nucleosynthesis

abstract: In this article, we study the implications of the coupling between Axion-Like-Particles (ALPs) and Leptons to cosmology, in particular, the Big Bang Nucleosynthesis (BBN). We show that the BBN, through the constraint on the effective number of relativistic neutrino species, provides the most stringent bound on the ALP-electron interaction strength for the mass of axion between 20 keV and 1 MeV. For other values of the mass, the BBN bound complements the stellar-evolution and laboratory bounds.


title: The Effect of Bars on the Ionized ISM: Optical Emission Lines from Milky Way Analogs

abstract: Gas interior to the bar of the Milky Way has recently been shown as the closest example of a Low Ionization (Nuclear) Emission Region--LI(N)ER--in the universe. To better understand the nature of this gas, a sample of face-on galaxies with integral field spectroscopy are used to study the ionized gas conditions of 240 barred and 250 nonbarred galaxies, focusing on those that are most similar to the Milky Way. Strong optical line emission of $[NII]$ $\lambda 6584$, H$\alpha$, $[OIII]$ $\lambda 5007$, and H$\beta$ are used to diagnose the dominant ionization mechanisms of gas across galaxies and the Galaxy via Baldwin-Phillips-Terlevich (BPT) Diagrams. Barred galaxies show a strong suppression of star formation and an increase in composite and LI(N)ER like spectra in their inner regions when compared with similar nonbarred counterparts. This effect is lessened in galaxies of very low ($\log_{10}(M_\star/M_\odot) \lesssim 10.4$) or very high ($\log_{10}(M_\star/M_\odot) \gtrsim 11.1$) total stellar mass. Bar masks from Galaxy Zoo:3D show the bar's non-axisymmetric effect on the ionized gas and help predict the face-on distribution of ionized gas conditions near the bar of the Milky Way.


title: Implementation of stellar heating feedback in simulations of star cluster formation: effects on the initial mass function

abstract: Explaining the initial mass function (IMF) of stars is a long-standing problem in astrophysics. The number of complex mechanisms involved in the process of star cluster formation, such as turbulence, magnetic fields and stellar feedback, make understanding and modeling the IMF a challenging task. In this paper, we aim to assert the importance of stellar heating feedback in the star cluster formation process and its effect on the shape of the IMF. We use an analytical sub-grid model to implement the radiative feedback in fully three-dimensional magnetohydrodynamical (MHD) simulations of star cluster formation, with the ultimate objective of obtaining numerical convergence on the IMF. We compare a set of MHD adaptive-mesh-refinement (AMR) simulations with three different implementations of the heating of the gas: 1) a polytropic equation of state (EOS), 2) a spherically symmetric stellar heating feedback, and 3) our newly developed polar heating model that takes into account the geometry of the accretion disc and the resulting shielding of stellar radiation by dust. For each of the three heating models, we analyse the distribution of stellar masses formed in ten molecular cloud simulations with different realizations of the turbulence to obtain a statistically representative IMF. We conclude that stellar heating feedback has a profound influence on the number of stars formed and plays a crucial role in controlling the IMF. We find that the simulations with the polar heating model achieve the best convergence on the observed IMF.


title: Selective branching, quenching, and converting of topological modes

abstract: A salient feature of topological phases are surface states and many of the widely studied physical properties are directly tied to their existence. Although less explored, a variety of topological phases can however similarly be distinguished by their response to localized flux defects, resulting in the binding of modes whose stability can be traced back to that of convectional edge states. The reduced dimensionality of these objects renders the possibility of arranging them in distinct geometries, such as arrays that branch or terminate in the bulk. We show that the prospect of hybridizing the modes in such new kinds of channels poses profound opportunities in a dynamical context. In particular, we find that creating junctions of $\pi$-flux chains or extending them as function of time can induce transistor and stop-and-go effects. Pending controllable initial conditions certain branches of the extended defect array can be actively biased. Discussing these physical effects within a generally applicable framework that relates to a variety of established artificial topological materials, such as mass-spring setups and LC circuits, our results offer an avenue to explore and manipulate new transport effects that are rooted in the topological characterization of the underlying system.


title: Ubiquitous velocity fluctuations throughout the molecular interstellar medium

abstract: The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.


title: The MUSE Hubble Ultra Deep Field Survey. XV. The mean rest-UV spectra of Ly-alpha emitters at z>3

abstract: We investigate the ultraviolet (UV) spectral properties of faint Lyman-$\alpha$ emitters (LAEs) in the redshift range 2.9<z<4.6 and provide material to prepare future observations of the faint Universe. We use data from the MUSE Hubble Ultra Deep Survey to construct mean rest-frame spectra of continuum-faint (median M$_{UV}$ of -18 and down to M$_{UV}$ of -16), low stellar mass (median value of $10^{8.4}$ and down to $10^{7}M_{\odot}$) LAEs at redshift z>3. We compute various averaged spectra of LAEs sub-sampled on the basis of their observational (e.g., Ly$\alpha$ strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We search for UV spectral features other than Ly$\alpha$, such as higher-ionization nebular emission lines and absorption features. We successfully observe the OIII]1666 and CIII]909 collisionally excited emission lines and the HeII1640 recombination feature, as well as the resonant CIV1550 doublet either in emission or P-Cygni. We compare the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Ly$\alpha$ equivalent widths, fainter UV magnitudes, bluer UV spectral slopes and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z>2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass and faint LAEs. We publicly release the stacked spectra as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope.


title: Gravitational Waves, Holography, and Black Hole Microstates

abstract: Gravitational wave observations of the near-horizon region of black holes lend insight into the quantum nature of gravity. In particular, gravitational wave echoes have been identified as a potential signature of quantum gravity-inspired structure near the horizon. In this paper, we connect such observables to the language of black hole microstates in string theory and holography. To that end, we propose a toy model describing the AdS$_3$ near-horizon region of five-dimensional black holes, inspired by earlier work of Solodukhin. This model captures key features of recently constructed microstate geometries, and allows us to make three observations. First, we relate the language of AdS/CFT, in particular the holographic retarded two-point correlator, to effective parameters modeling the structure that are used in flat space gravitational wave literature. Second, we find that for a typical microstate, the `cap' of the microstructure is exponentially close to the horizon, making it an effective sub-Planckian correction to the black hole geometry, although the microstate geometry itself is classical. Third, using a microcanonical ensemble average over geometries, we find support for the claim that the gravitational wave echo amplitude in a typical quantum microstate of the black hole is exponentially suppressed by the black hole entropy.


title: Measuring the total infrared light from galaxy clusters at z=0.5-1.6: connecting stellar populations to dusty star formation

abstract: Massive galaxy clusters undergo strong evolution from z~1.6 to z~0.5, with overdense environments at high-z characterized by abundant dust-obscured star formation and stellar mass growth which rapidly give way to widespread quenching. Data spanning the near- to far-infrared (IR) spectrum can directly trace this transformation; however, such studies have largely been limited to the massive galaxy end of cluster populations. In this work, we present ``total light" stacking techniques spanning 3.4-500{\mu}m aimed at revealing the total cluster IR emission, including low mass members and potential intracluster dust. We detail our procedures for WISE, Spitzer, and Herschel imaging, including corrections to recover the total stacked emission in the case of high fractions of detected galaxies. We apply our stacking techniques to 232 well-studied massive (log M200/Msun~13.8) clusters across multiple z bins, recovering extended cluster emission at all wavelengths, typically at >5sigma. We measure the averaged near- to far-IR radial profiles and SEDs, quantifying the total stellar and dust content. The near-IR radial profiles are well described by an NFW model with a high (c~7) concentration parameter. Dust emission is similarly concentrated, albeit suppressed at small radii (r<0.2Mpc). The measured SEDs lack warm dust, consistent with the colder SEDs expected for low mass galaxies. We derive total stellar masses consistent with the theoretical Mhalo-M_star relation and specific-star formation rates that evolve strongly with redshift, echoing that of massive (log Mstar/Msun>10) cluster galaxies. Separating out the massive galaxy population reveals that the majority of cluster far-IR emission (~70-80%) is provided by the low mass constituents, which differs from field galaxies. This effect may be a combination of mass-dependent quenching and excess dust in low mass cluster galaxies.


title: Classical-quantum correspondence for two-level pseudo-Hermitian systems

abstract: In this work, a classical-quantum correspondence for two-level pseudo-Hermitian systems is proposed and analyzed. We show that the presence of a complex external field can be described by a pseudo-Hermitian Hamiltonian if there is a suitable canonical transformation that links it to a real field. We construct a covariant quantization scheme which maps canonically related pseudoclassical theories to unitarily equivalent quantum realizations, such that there is a unique metric-inducing isometry between the distinct Hilbert spaces. In this setting, the pseudo-Hermiticity condition for the operators induces an involution which guarantees the reality of the corresponding symbols, even for the complex field case. We assign a physical meaning for the dynamics in the presence of a complex field by constructing a classical correspondence. As an application of our theoretical framework, we propose a damped version of the Rabi problem and determine the configuration of the parameters of the setup for which damping is completely suppressed. The experimental viability of the proposal is studied within a specific context. We suggest that the main theoretical results developed in the present work could be experimentally verified.

categories: quant-ph hep-th math-ph math.MP physics.atom-ph


title: Quantum fluctuations hinder finite-time information erasure near the Landauer limit

abstract: Information is physical but information is also processed in finite time. Where computing protocols are concerned, finite-time processing in the quantum regime can dynamically generate coherence. Here we show that this can have significant thermodynamic implications. We demonstrate that quantum coherence generated in the energy eigenbasis of a system undergoing a finite-time information erasure protocol yields rare events with extreme dissipation. These fluctuations are of purely quantum origin. By studying the full statistics of the dissipated heat in the slow driving limit, we prove that coherence provides a non-negative contribution to all statistical cumulants. Using the simple and paradigmatic example of single bit erasure, we show that these extreme dissipation events yield distinct, experimentally distinguishable signatures.


title: Egocentric Action Recognition by Video Attention and Temporal Context

abstract: We present the submission of Samsung AI Centre Cambridge to the CVPR2020 EPIC-Kitchens Action Recognition Challenge. In this challenge, action recognition is posed as the problem of simultaneously predicting a single `verb' and `noun' class label given an input trimmed video clip. That is, a `verb' and a `noun' together define a compositional `action' class. The challenging aspects of this real-life action recognition task include small fast moving objects, complex hand-object interactions, and occlusions. At the core of our submission is a recently-proposed spatial-temporal video attention model, called `W3' (`What-Where-When') attention~\cite{perez2020knowing}. We further introduce a simple yet effective contextual learning mechanism to model `action' class scores directly from long-term temporal behaviour based on the `verb' and `noun' prediction scores. Our solution achieves strong performance on the challenge metrics without using object-specific reasoning nor extra training data. In particular, our best solution with multimodal ensemble achieves the 2$^{nd}$ best position for `verb', and 3$^{rd}$ best for `noun' and `action' on the Seen Kitchens test set.


title: High-recall causal discovery for autocorrelated time series with latent confounders

abstract: We present a new method for linear and nonlinear, lagged and contemporaneous constraint-based causal discovery from observational time series in the presence of latent confounders. We show that existing causal discovery methods such as FCI and variants suffer from low recall in the autocorrelated time series case and identify low effect size of conditional independence tests as the main reason. Information-theoretical arguments show that effect size can often be increased if causal parents are included in the conditioning sets. To identify parents early on, we suggest an iterative procedure that utilizes novel orientation rules to determine ancestral relationships already during the edge removal phase. We prove that the method is order-independent, and sound and complete in the oracle case. Extensive simulation studies for different numbers of variables, time lags, sample sizes, and further cases demonstrate that our method indeed achieves much higher recall than existing methods for the case of autocorrelated continuous variables while keeping false positives at the desired level. This performance gain grows with stronger autocorrelation. At https://github.com/jakobrunge/tigramite we provide Python code for all methods involved in the simulation studies.


title: Submillimetre galaxies in cosmological hydrodynamical simulations -- an opportunity for constraining feedback models

abstract: Submillimetre galaxies (SMGs) have long posed a challenge for theorists, and self-consistently reproducing the properties of the SMG population in a large-volume cosmological hydrodynamical simulation has not yet been achieved. We use a scaling relation derived from previous simulations plus radiative transfer calculations to predict the submm flux densities of simulated SMGs drawn from cosmological simulations from the Illustris and IllustrisTNG projects based on the simulated galaxies' star formation rates (SFRs) and dust masses and compare the predicted number counts with observations. We find that the predicted SMG number counts based on IllustrisTNG are significantly less than observed (more than 1 dex at $S_{850} \gtrsim 4$ mJy). The simulation from the original Illustris project yields more SMGs than IllustrisTNG: the predicted counts are consistent with those observed at both $S_{850} \lesssim 5$ mJy and $S_{850} \gtrsim 9$ mJy and only a factor of $\sim 2$ lower than observed at intermediate flux densities. The redshift distribution of SMGs with $S_{850} > 3$ mJy in IllustrisTNG is consistent with the observed distribution, whereas the Illustris redshift distribution peaks at significantly lower redshift (1.5 vs. 2.8). We demonstrate that IllustrisTNG hosts fewer SMGs than Illustris because in the former, high-mass ($M_{\star} \sim 10^{11} \, \text{M}_{\odot}$) $z \sim 2-3$ galaxies have lower dust masses and SFRs than in Illustris owing to differences in the sub-grid models for stellar or/and active galactic nucleus (AGN) feedback between the two simulations (we unfortunately cannot isolate the specific cause(s) post hoc). Our results demonstrate that because our method enables predicting SMG number counts in post-processing with a negligible computational expense, SMGs can provide useful constraints for tuning sub-grid models in future large-volume cosmological simulations.


title: Nonrelativistic Open String and Yang-Mills Theory

abstract: The classical and quantum worldsheet theory describing nonrelativistic open string theory in an arbitrary nonrelativistic open and closed string background is constructed. We show that the low energy dynamics of open strings ending on n coincident D-branes in flat spacetime is described by a Galilean invariant U(n) Yang-Mills theory. We also study nonrelativistic open string excitations with winding number and demonstrate that their dynamics can be encoded into a local gauge theory in one higher dimension. By demanding conformal invariance of the boundary couplings, the nonlinear equations of motion that govern the consistent open string backgrounds coupled to an arbitrary closed background (described by a string Newton-Cartan geometry, Kalb-Ramond, and dilaton field) are derived and shown to emerge from a Galilean invariant Dirac-Born-Infeld type action.


title: Supernova Triggers for End-Devonian Extinctions

abstract: The Late Devonian was a protracted period of low speciation resulting in biodiversity decline, culminating in extinction events near the Devonian-Carboniferous boundary. Recent evidence indicates that the final extinction event may have coincided with a dramatic drop in stratospheric ozone, possibly due to a global temperature rise. Here we study an alternative possible cause for the postulated ozone drop: a nearby supernova explosion that could inflict damage by accelerating cosmic rays that can deliver ionizing radiation for up to $\sim 100$ kyr. We therefore propose that the end-Devonian extinctions were triggered by supernova explosions at $\sim 20$ pc, somewhat beyond the "kill distance" that would have precipitated a full mass extinction. Such nearby supernovae are likely due to core-collapses of massive stars; these are concentrated in the thin Galactic disk where the Sun resides. Detecting either of the long-lived radioisotopes Sm-146 or Pu-244 in one or more end-Devonian extinction strata would confirm a supernova origin, point to the core-collapse explosion of a massive star, and probe supernova nucleosythesis. Other possible tests of the supernova hypothesis are discussed.


title: Inference on the change point in high dimensional time series models via plug in least squares

abstract: We study a plug in least squares estimator for the change point parameter where change is in the mean of a high dimensional random vector under subgaussian or subexponential distributions. We obtain sufficient conditions under which this estimator possesses sufficient adaptivity against plug in estimates of mean parameters in order to yield an optimal rate of convergence $O_p(\xi^{-2})$ in the integer scale. This rate is preserved while allowing high dimensionality as well as a potentially diminishing jump size $\xi,$ provided $s\log (p\vee T)=o(\surd(Tl_T))$ or $s\log^{3/2}(p\vee T)=o(\surd(Tl_T))$ in the subgaussian and subexponential cases, respectively. Here $s,p,T$ and $l_T$ represent a sparsity parameter, model dimension, sampling period and the separation of the change point from its parametric boundary. Moreover, since the rate of convergence is free of $s,p$ and logarithmic terms of $T,$ it allows the existence of limiting distributions. These distributions are then derived as the {\it argmax} of a two sided negative drift Brownian motion or a two sided negative drift random walk under vanishing and non-vanishing jump size regimes, respectively. Thereby allowing inference of the change point parameter in the high dimensional setting. Feasible algorithms for implementation of the proposed methodology are provided. Theoretical results are supported with monte-carlo simulations.


title: Joint galaxy-galaxy lensing and clustering constraints on galaxy formation

abstract: We compare predictions for galaxy-galaxy lensing profiles and clustering from the Henriques et al. (2015) public version of the Munich semi-analytical model of galaxy formation (SAM) and the IllustrisTNG suite, primarily TNG300, with observations from KiDS+GAMA and SDSS-DR7 using four different selection functions for the lenses (stellar mass, stellar mass and group membership, stellar mass and isolation criteria, stellar mass and colour). We find that this version of the SAM does not agree well with the current data for stellar mass-only lenses with $M_\ast > 10^{11}\,M_\odot$. By decreasing the merger time for satellite galaxies as well as reducing the radio-mode AGN accretion efficiency in the SAM, we obtain better agreement, both for the lensing and the clustering, at the high mass end. We show that the new model is consistent with the signals for central galaxies presented in Velliscig et al. (2017). Turning to the hydrodynamical simulation, TNG300 produces good lensing predictions, both for stellar mass-only ($\chi^2 = 1.81$ compared to $\chi^2 = 7.79$ for the SAM), and locally brightest galaxies samples ($\chi^2 = 3.80$ compared to $\chi^2 = 5.01$). With added dust corrections to the colours it matches the SDSS clustering signal well for red low mass galaxies. We find that both the SAMs and TNG300 predict $\sim 50\,\%$ excessive lensing signals for intermediate mass red galaxies with $10.2 < \log_{10} M_\ast [ M_\odot ] < 11.2$ at $r \approx 0.6\,h^{-1}\,\mathrm{Mpc}$, which require further theoretical development.


title: Bypassing the energy functional in density functional theory: Direct calculation of electronic energies from conditional probability densities

abstract: Density functional calculations can fail for want of an accurate exchange-correlation approximation. The energy can instead be extracted from a sequence of density functional calculations of conditional probabilities (CP-DFT). Simple CP approximations yield usefully accurate results for two-electron ions, the hydrogen dimer, and the uniform gas at all temperatures. CP-DFT has no self-interaction error for one electron, and correctly dissociates H2, both major challenges. For warm dense matter, classical CP-DFT calculations can overcome the convergence problems of Kohn-Sham DFT.


title: A Unifying View of Optimism in Episodic Reinforcement Learning

abstract: The principle of optimism in the face of uncertainty underpins many theoretically successful reinforcement learning algorithms. In this paper we provide a general framework for designing, analyzing and implementing such algorithms in the episodic reinforcement learning problem. This framework is built upon Lagrangian duality, and demonstrates that every model-optimistic algorithm that constructs an optimistic MDP has an equivalent representation as a value-optimistic dynamic programming algorithm. Typically, it was thought that these two classes of algorithms were distinct, with model-optimistic algorithms benefiting from a cleaner probabilistic analysis while value-optimistic algorithms are easier to implement and thus more practical. With the framework developed in this paper, we show that it is possible to get the best of both worlds by providing a class of algorithms which have a computationally efficient dynamic-programming implementation and also a simple probabilistic analysis. Besides being able to capture many existing algorithms in the tabular setting, our framework can also address largescale problems under realizable function approximation, where it enables a simple model-based analysis of some recently proposed methods.


title: Generalized Path Pairs and Fuss-Catalan Triangles

abstract: Path pairs are a modification of parallelogram polyominoes that provide yet another combinatorial interpretation of the Catalan numbers. More generally, the number of path pairs of length $n$ and distance $\delta$ corresponds to the $(n-1,\delta-1)$ entry of Shapiro's so-called Catalan triangle. In this paper, we widen the notion of path pairs $(\gamma_1,\gamma_2)$ to the situation where $\gamma_1$ and $\gamma_2$ may have different lengths, and then enforce divisibility conditions on runs of vertical steps in $\gamma_2$. This creates a two-parameter family of integer triangles that generalize the Catalan triangle and qualify as proper Riordan arrays for many choices of parameters. In particular, we use generalized path pairs to provide a new combinatorial interpretation for all entries in every proper Riordan array $\mathcal{R}(d(t),h(t))$ of the form $d(t) = C_k(t)^i$, $h(t) = t \kern+1pt C_k(t)^k$, where $1 \leq i \leq k$ and $C_k(t)$ is the generating function for some sequence of Fuss-Catalan numbers (some $k \geq 2$). Closed formulas are then provided for the number of generalized path pairs across an even broader range of parameters, as well as for the number of weak path pairs with a fixed number of non-initial intersections.


title: The $\nu_{R}$-philic scalar: its loop-induced interactions and Yukawa forces in LIGO observations

abstract: Right-handed neutrinos ($\nu_{R}$) are often considered as a portal to new hidden physics. It is tempting to consider a gauge singlet scalar $(\phi)$ that exclusively couples to $\nu_{R}$ via a $\nu_{R}\nu_{R}\phi$ term. Such a $\nu_{R}$-philic scalar does not interact with charged fermions at tree level but loop-induced effective interactions are inevitable, which are systematically investigated in this work. The magnitude of the loop-induced couplings coincidentally meets the current sensitivity of fifth-force searches. In particular, the loop-induced coupling to muons could be tested in the recent LIGO observations of neutron star mergers as there might be a sizable Yukawa force in the binary system mediated by the $\nu_{R}$-philic scalar.


title: Integral representation for Euler sums of hyperharmonic numbers

abstract: In this short paper, we derive an integral representation for Euler sums of hyperharmonic numbers. We use results established by other authors to then show that the integral has a closed-form in terms of zeta values and Stirling numbers of the first kind. Specifically, the integral has the form of $$\int_0^\infty \frac{t^{m-1}\ln(1-e^{-t})}{(1-e^{-t})^r} \ dt$$ where $m, r \in \mathbb{N}$, $m > r$ and $r\ge1$.


title: On 3-distance spherical 5-designs

abstract: Inspired by a recently formulated conjecture by Bannai et al. we investigate spherical codes which admit exactly three different distances and are spherical 5-designs. Computing and analyzing distance distributions we provide new proof of the fact (due to Levenshtein) that such codes are maximal and rule out certain cases towards a proof of the conjecture.


title: Spatial Iterated Prisoner's Dilemma as a Transformation Semigroup

abstract: The prisoner's dilemma (PD) is a game-theoretic model studied in a wide array of fields to understand the emergence of cooperation between rational self-interested agents. In this work, we formulate a spatial iterated PD as a discrete-event dynamical system where agents play the game in each time-step and analyse it algebraically using Krohn-Rhodes algebraic automata theory using a computational implementation of the holonomy decomposition of transformation semigroups. In each iteration all players adopt the most profitable strategy in their immediate neighbourhood. Perturbations resetting the strategy of a given player provide additional generating events for the dynamics. Our initial study shows that the algebraic structure, including how natural subsystems comprising permutation groups acting on the spatial distributions of strategies, arise in certain parameter regimes for the pay-off matrix, and are absent for other parameter regimes. Differences in the number of group levels in the holonomy decomposition (an upper bound for Krohn-Rhodes complexity) are revealed as more pools of reversibility appear when the temptation to defect is at an intermediate level. Algebraic structure uncovered by this analysis can be interpreted to shed light on the dynamics of the spatial iterated PD.


title: b-articulation points and b-bridges in strongly biconnected directed graphs

abstract: A directed graph $G=(V,E)$ is called strongly biconnected if $G$ is strongly connected and the underlying graph of $G$ is biconnected. This class of directed graphs was first introduced by Wu and Grumbach. Let $G=(V,E)$ be a strongly biconnected directed graph. An edge $e\in E$ is a b-bridge if the subgraph $G\setminus \left\lbrace e\right\rbrace =(V,E\setminus \left\lbrace e\right\rbrace) $ is not strongly biconnected. A vertex $w\in V$ is a b-articulation point if $G\setminus \left\lbrace w\right\rbrace$ is not strongly biconnected, where $G\setminus \left\lbrace w\right\rbrace$ is the subgraph obtained from $G$ by removing $w$. In this paper we study b-articulation points and b-bridges.