title: Two Resolutions of the Margin Loan Pricing Puzzle

abstract: This paper supplies two possible resolutions of Fortune's (2000) margin-loan pricing puzzle. Fortune (2000) noted that the margin loan interest rates charged by stock brokers are very high in relation to the actual (low) credit risk and the cost of funds. If we live in the Black-Scholes world, the brokers are presumably making arbitrage profits by shorting dynamically precise amounts of their clients' portfolios. First, we extend Fortune's (2000) application of Merton's (1974) no-arbitrage approach to allow for brokers that can only revise their hedges finitely many times during the term of the loan. We show that extremely small differences in the revision frequency can easily explain the observed variation in margin loan pricing. In fact, four additional revisions per three-day period serve to explain all of the currently observed heterogeneity. Second, we study monopolistic (or oligopolistic) margin loan pricing by brokers whose clients are continuous-time Kelly gamblers. The broker solves a general stochastic control problem that yields simple and pleasant formulas for the optimal interest rate and the net interest margin. If the author owned a brokerage, he would charge an interest rate of $(r+\nu)/2-\sigma^2/4$, where $r$ is the cost of funds, $\nu$ is the compound-annual growth rate of the S&P 500 index, and $\sigma$ is the volatility.

categories: econ.GN econ.TH q-fin.EC q-fin.GN q-fin.PM q-fin.RM


title: NodeDrop: A Condition for Reducing Network Size without Effect on Output

abstract: Determining an appropriate number of features for each layer in a neural network is an important and difficult task. This task is especially important in applications on systems with limited memory or processing power. Many current approaches to reduce network size either utilize iterative procedures, which can extend training time significantly, or require very careful tuning of algorithm parameters to achieve reasonable results. In this paper we propose NodeDrop, a new method for eliminating features in a network. With NodeDrop, we define a condition to identify and guarantee which nodes carry no information, and then use regularization to encourage nodes to meet this condition. We find that NodeDrop drastically reduces the number of features in a network while maintaining high performance, reducing the number of parameters by a factor of 114x for a VGG like network on CIFAR10 without a drop in accuracy.


title: Wave breaking for shallow water models with time decaying solutions

abstract: A family of Camassa-Holm type equations with a linear term and cubic and quartic nonlinearities is considered. Local well-posedness results are established via Kato's approach. Conserved quantities for the equation are determined and from them we prove that the energy functional of the solutions is a time-dependent, monotonically decreasing function of time, and bounded from above by the Sobolev norm of the initial data under some conditions. The existence of wave breaking phenomenon is investigated and necessary conditions for its existence are obtained. In our framework the wave breaking is guaranteed, among other conditions, when the coefficient of the linear term is sufficiently small, which allows us to interpret the equation as a linear perturbation of some recent Camassa-Holm type equations considered in the literature.


title: A Hybrid RNN-HMM Approach for Weakly Supervised Temporal Action Segmentation

abstract: Action recognition has become a rapidly developing research field within the last decade. But with the increasing demand for large scale data, the need of hand annotated data for the training becomes more and more impractical. One way to avoid frame-based human annotation is the use of action order information to learn the respective action classes. In this context, we propose a hierarchical approach to address the problem of weakly supervised learning of human actions from ordered action labels by structuring recognition in a coarse-to-fine manner. Given a set of videos and an ordered list of the occurring actions, the task is to infer start and end frames of the related action classes within the video and to train the respective action classifiers without any need for hand labeled frame boundaries. We address this problem by combining a framewise RNN model with a coarse probabilistic inference. This combination allows for the temporal alignment of long sequences and thus, for an iterative training of both elements. While this system alone already generates good results, we show that the performance can be further improved by approximating the number of subactions to the characteristics of the different action classes as well as by the introduction of a regularizing length prior. The proposed system is evaluated on two benchmark datasets, the Breakfast and the Hollywood extended dataset, showing a competitive performance on various weak learning tasks such as temporal action segmentation and action alignment.


title: Dealing With $\delta$-Scuti Variables: Transit Light Curve Analysis of Planets Orbiting Rapidly-Rotating, Seismically Active A/F Stars

abstract: We measure the bulk system parameters of the seismically active, rapidly-rotating $\delta$-Scuti KOI-976 and constrain the orbit geometry of its transiting binary companion using a combined approach of asteroseismology and gravity-darkening light curve analysis. KOI-976 is a $1.62\pm0.2~\mathrm{M_\odot}$ star with a measured $v\sin(i)$ of $120\pm2$ km/s and seismically-induced variable signal that varies by $\sim$ 0.6\% of the star's total photometric brightness. We take advantage of the star's oblate shape and seismic activity to perform three measurements of its obliquity angle relative to the plane of the sky. We first apply rotational splitting theory to the star's variable signal observed in short-cadence \emph{Kepler} photometry to constrain KOI-976's obliquity angle, and then subtract off variability from that dataset using the linear algorithm for significance reduction software {\tt LASR}. We perform gravity-darkened fits to \emph{Kepler} variability-subtracted short-cadence photometry and to \emph{Kepler's} phase-folded long-cadence photometry to obtain two more measurements of the star's obliquity. We find that the binary system transits in a grazing configuration with measured obliquity values of $36^\circ\pm17^\circ$, $46^\circ\pm16^\circ$, and $43^\circ\pm20^\circ$ respectively for the three measurements. We perform these analyses as a way to demonstrate overcoming the challenges high-mass stars can present to transit light curve fitting and to prepare for the large number of exoplanets \emph{TESS} will discover orbiting A/F stars.


title: Correctness Verification of Neural Networks

abstract: We present a novel framework for specifying and verifying correctness globally for neural networks on perception tasks. Most previous works on neural network verification for perception tasks focus on robustness verification. Unlike robustness verification, which aims to verify that the prediction of a network is stable in some local regions around labelled points, our framework provides a way to specify correctness globally in the whole target input space and verify that the network is correct for all target inputs (or find the regions where the network is not correct). We provide a specification through 1) a state space consisting of all relevant states of the world and 2) an observation process that produces neural network inputs from the states of the world. Tiling the state and input spaces with a finite number of tiles, obtaining ground truth bounds from the state tiles and network output bounds from the input tiles, then comparing the ground truth and network output bounds delivers an upper bound on the network output error for any inputs of interest. The presented framework also enables detecting illegal inputs -- inputs that are not contained in (or close to) the target input space as defined by the state space and observation process (the neural network is not designed to work on them), so that we can flag when we don't have guarantees. Results from two case studies highlight the ability of our technique to verify error bounds over the whole target input space and show how the error bounds vary over the state and input spaces.


title: Probabilistic Existence Results for Parent-Identifying Schemes

abstract: Parent-identifying schemes provide a way to identify causes from effects for some information systems such as digital fingerprinting and group testing. In this paper, we consider combinatorial structures for parent-identifying schemes. First, we establish an equivalent relationship between parent-identifying schemes and forbidden configurations. Based on this relationship, we derive probabilistic existence lower bounds for two related combinatorial structures, that is, $t$-parent-identifying set systems ($t$-IPPS) and $t$-multimedia parent-identifying codes ($t$-MIPPC), which are used in broadcast encryption and multimedia fingerprinting respectively. The probabilistic lower bound for the maximum size of a $t$-IPPS has the asymptotically optimal order of magnitude in many cases, and that for $t$-MIPPC provides the asymptotically optimal code rate when $t=2$ and the best known asymptotic code rate when $t\geq 3$. Furthermore, we analyze the structure of $2$-IPPS and prove some bounds for certain cases.


title: A Language-Agnostic Model for Semantic Source Code Labeling

abstract: Code search and comprehension have become more difficult in recent years due to the rapid expansion of available source code. Current tools lack a way to label arbitrary code at scale while maintaining up-to-date representations of new programming languages, libraries, and functionalities. Comprehensive labeling of source code enables users to search for documents of interest and obtain a high-level understanding of their contents. We use Stack Overflow code snippets and their tags to train a language-agnostic, deep convolutional neural network to automatically predict semantic labels for source code documents. On Stack Overflow code snippets, we demonstrate a mean area under ROC of 0.957 over a long-tailed list of 4,508 tags. We also manually validate the model outputs on a diverse set of unlabeled source code documents retrieved from Github, and we obtain a top-1 accuracy of 86.6%. This strongly indicates that the model successfully transfers its knowledge from Stack Overflow snippets to arbitrary source code documents.


title: The Mechanism of the Origin and Development of Lightning from Initiating Event to Initial Breakdown Pulses (v.2)

abstract: Based on experimental results of recent years, this article presents a qualitative description of a possible mechanism (termed the Mechanism) covering the main stages of lightning initiation, starting before and including the initiating event, followed by the initial electric field change (IEC), followed by the first few initial breakdown pulses (IBPs). The Mechanism assumes initiation occurs in a region of ~1 km3 with average electric field E>0.3 MV/(m atm), which contains, because of turbulence, numerous small 'Eth-volumes' of 10^-4-10^-3 m3 with E>3 MV/(m atm). The Mechanism allows for lightning initiation by either of two observed types of events: a high power VHF event such as a Narrow Bipolar Event, or a weak VHF event. According to the Mechanism, both types of initiating events are caused by a group of relativistic runaway electron avalanche particles (where the initial electrons are secondary particles of an extensive air shower) passing through many Eth-volumes, thereby causing the nearly simultaneous launching of many positive streamer flashes. Due to ionization-heating instability, unusual plasma formations (UPFs) appear along the streamers' trajectories of the streamers. These UPFs combine into three-dimensional (3D) networks of hot plasma channels during the IEC, resulting in its observed weak current flow. The subsequent development and combination of two (or more) of these 3D networks of hot plasma channels then causes the first IBP. Each subsequent IBP is caused when another 3D network of hot plasma channels combines with the chain of networks caused by earlier IBPs.


title: Modelling double emulsion formation in planar flow-focusing microchannels

abstract: Double emulsion formation in a hierarchical flow-focusing channel is systematically investigated using a free energy ternary lattice Boltzmann model. A three dimensional formation regime diagram is constructed based on the capillary numbers of the inner ($Ca_i$), middle ($Ca_m$) and outer ($Ca_o$) phase fluids. The results show that the formation diagram can be classified into periodic two-step region, periodic one-step region, and non-periodic region. By varying $Ca_i$ and $Ca_m$ in the two-step formation region, different morphologies are obtained, including the regular double emulsions, decussate regimes with one or two alternate empty droplets, and structures with multiple inner droplets contained in the continuous middle phase thread. Bidisperse behaviors are also frequently encountered in the two-step formation region. In the periodic one-step formation region, scaling laws are proposed for the double emulsion size and for the size ratio between the inner droplet and the overall double emulsion. Furthermore, we show that the interfacial tension ratio can greatly change the morphologies of the obtained emulsion droplets, and the channel geometry plays an important role in changing the formation regimes and the double emulsion sizes. In particular, narrowing the side inlets or the distance between the two side inlets promotes the conversion from the two-step formation regime to the one-step formation regime.


title: A Perspective on Objects and Systematic Generalization in Model-Based RL

abstract: In order to meet the diverse challenges in solving many real-world problems, an intelligent agent has to be able to dynamically construct a model of its environment. Objects facilitate the modular reuse of prior knowledge and the combinatorial construction of such models. In this work, we argue that dynamically bound features (objects) do not simply emerge in connectionist models of the world. We identify several requirements that need to be fulfilled in overcoming this limitation and highlight corresponding inductive biases.


title: Transient amplifiers of selection and reducers of fixation for death-Birth updating on graphs

abstract: The spatial structure of an evolving population affects which mutations become fixed. Some structures amplify selection, increasing the likelihood that beneficial mutations become fixed while deleterious mutations do not. Other structures suppress selection, reducing the effect of fitness differences and increasing the role of random chance. This phenomenon can be modeled by representing spatial structure as a graph, with individuals occupying vertices. Births and deaths occur stochastically, according to a specified update rule. We study death-Birth updating: An individual is chosen to die and then its neighbors compete to reproduce into the vacant spot. Previous numerical experiments suggested that amplifiers of selection for this process are either rare or nonexistent. We introduce a perturbative method for this problem for weak selection regime, meaning that mutations have small fitness effects. We show that fixation probability under weak selection can be calculated in terms of the coalescence times of random walks. This result leads naturally to a new definition of effective population size. Using this and other methods, we uncover the first known examples of transient amplifiers of selection (graphs that amplify selection for a particular range of fitness values) for the death-Birth process. We also exhibit new families of "reducers of fixation", which decrease the fixation probability of all mutations, whether beneficial or deleterious.


title: Better Character Language Modeling Through Morphology

abstract: We incorporate morphological supervision into character language models (CLMs) via multitasking and show that this addition improves bits-per-character (BPC) performance across 24 languages, even when the morphology data and language modeling data are disjoint. Analyzing the CLMs shows that inflected words benefit more from explicitly modeling morphology than uninflected words, and that morphological supervision improves performance even as the amount of language modeling data grows. We then transfer morphological supervision across languages to improve language modeling performance in the low-resource setting.


title: Transforming Complex Sentences into a Semantic Hierarchy

abstract: We present an approach for recursively splitting and rephrasing complex English sentences into a novel semantic hierarchy of simplified sentences, with each of them presenting a more regular structure that may facilitate a wide variety of artificial intelligence tasks, such as machine translation (MT) or information extraction (IE). Using a set of hand-crafted transformation rules, input sentences are recursively transformed into a two-layered hierarchical representation in the form of core sentences and accompanying contexts that are linked via rhetorical relations. In this way, the semantic relationship of the decomposed constituents is preserved in the output, maintaining its interpretability for downstream applications. Both a thorough manual analysis and automatic evaluation across three datasets from two different domains demonstrate that the proposed syntactic simplification approach outperforms the state of the art in structural text simplification. Moreover, an extrinsic evaluation shows that when applying our framework as a preprocessing step the performance of state-of-the-art Open IE systems can be improved by up to 346% in precision and 52% in recall. To enable reproducible research, all code is provided online.


title: Neural networks grown and self-organized by noise

abstract: Living neural networks emerge through a process of growth and self-organization that begins with a single cell and results in a brain, an organized and functional computational device. Artificial neural networks, however, rely on human-designed, hand-programmed architectures for their remarkable performance. Can we develop artificial computational devices that can grow and self-organize without human intervention? In this paper, we propose a biologically inspired developmental algorithm that can 'grow' a functional, layered neural network from a single initial cell. The algorithm organizes inter-layer connections to construct a convolutional pooling layer, a key constituent of convolutional neural networks (CNN's). Our approach is inspired by the mechanisms employed by the early visual system to wire the retina to the lateral geniculate nucleus (LGN), days before animals open their eyes. The key ingredients for robust self-organization are an emergent spontaneous spatiotemporal activity wave in the first layer and a local learning rule in the second layer that 'learns' the underlying activity pattern in the first layer. The algorithm is adaptable to a wide-range of input-layer geometries, robust to malfunctioning units in the first layer, and so can be used to successfully grow and self-organize pooling architectures of different pool-sizes and shapes. The algorithm provides a primitive procedure for constructing layered neural networks through growth and self-organization. Broadly, our work shows that biologically inspired developmental algorithms can be applied to autonomously grow functional 'brains' in-silico.


title: A Surprising Density of Illusionable Natural Speech

abstract: Recent work on adversarial examples has demonstrated that most natural inputs can be perturbed to fool even state-of-the-art machine learning systems. But does this happen for humans as well? In this work, we investigate: what fraction of natural instances of speech can be turned into "illusions" which either alter humans' perception or result in different people having significantly different perceptions? We first consider the McGurk effect, the phenomenon by which adding a carefully chosen video clip to the audio channel affects the viewer's perception of what is said (McGurk and MacDonald, 1976). We obtain empirical estimates that a significant fraction of both words and sentences occurring in natural speech have some susceptibility to this effect. We also learn models for predicting McGurk illusionability. Finally we demonstrate that the Yanny or Laurel auditory illusion (Pressnitzer et al., 2018) is not an isolated occurrence by generating several very different new instances. We believe that the surprising density of illusionable natural speech warrants further investigation, from the perspectives of both security and cognitive science. Supplementary videos are available at: https://www.youtube.com/playlist?list=PLaX7t1K-e_fF2iaenoKznCatm0RC37B_k.


title: Dihedral Invariant Polynomials in the effective Lagrangian of QED

abstract: We present a new group-theoretical technique to calculate weak field expansions for some Feynman diagrams using invariant polynomials of the dihedral group. In particular we show results obtained for the first coefficients of the three loop effective lagrangian of 1+1 QED in an external constant field, where the dihedral symmetry appears. Our results suggest that a closed form involving rational numbers and the Riemann zeta function might exist for these coefficients.


title: Magnon-phonon interactions in magnetic insulators

abstract: We address the theory of magnon-phonon interactions and compute the corresponding quasi-particle and transport lifetimes in magnetic insulators with focus on yttrium iron garnet at intermediate temperatures from anisotropy- and exchange-mediated magnon-phonon interactions, the latter being derived from the volume dependence of the Curie temperature. We find in general weak effects of phonon scattering on magnon transport and the Gilbert damping of the macrospin Kittel mode. The magnon transport lifetime differs from the quasi-particle lifetime at shorter wavelengths.


title: Classification of exterior and proper fibrations

abstract: We classify exterior fibrations in the exterior homotopy category. As a result we also classify proper fibrations between CW-complexes.


title: Weakly Supervised Disentanglement by Pairwise Similarities

abstract: Recently, researches related to unsupervised disentanglement learning with deep generative models have gained substantial popularity. However, without introducing supervision, there is no guarantee that the factors of interest can be successfully recovered. Motivated by a real-world problem, we propose a setting where the user introduces weak supervision by providing similarities between instances based on a factor to be disentangled. The similarity is provided as either a binary (yes/no) or a real-valued label describing whether a pair of instances are similar or not. We propose a new method for weakly supervised disentanglement of latent variables within the framework of Variational Autoencoder. Experimental results demonstrate that utilizing weak supervision improves the performance of the disentanglement method substantially.


title: Fault-Tolerant Quantum Gates with Defects in Topological Stabiliser Codes

abstract: Braiding defects in topological stabiliser codes has been widely studied as a promising approach to fault-tolerant quantum computing. Here, we explore the potential and limitations of such schemes in codes of all spatial dimensions. We prove that a universal gate set for quantum computing cannot be realised by supplementing locality-preserving logical operators with defect braiding, even in more than two dimensions. However, notwithstanding this no-go theorem, we demonstrate that higher dimensional defect-braiding schemes have the potential to play an important role in realising fault-tolerant quantum computing. Specifically, we present an approach to implement the full Clifford group via braiding in any code possessing twist defects on which a fermion can condense. We explore three such examples in higher dimensional codes, specifically: in self-dual surface codes; the three dimensional Levin-Wen fermion mode; and the checkerboard model. Finally, we show how our no-go theorems can be circumvented to provide a universal scheme in three-dimensional surface codes without magic state distillation. Specifically, our scheme employs adaptive implementation of logical operators conditional on logical measurement outcomes to lift a combination of locality-preserving and braiding logical operators to universality.


title: Adiabatic projection method with Euclidean time subspace projection

abstract: Euclidean time projection is a powerful tool that uses exponential decay to extract the low-energy information of quantum systems. The adiabatic projection method, which is based on Euclidean time projection, is a procedure for studying scattering and reactions on the lattice. The method constructs the adiabatic Hamiltonian that gives the low-lying energies and wave functions of two-cluster systems. In this paper we seek the answer to the question whether an adiabatic Hamiltonian constructed in a smaller subspace of the two-cluster state space can still provide information on the low-lying spectrum and the corresponding wave functions. We present the results from our investigations on constructing the adiabatic Hamiltonian using Euclidean time projection and extracting details of the low-energy spectrum and wave functions by diagonalizing it. In our analyses we consider systems of fermion-fermion and fermion-dimer interacting via a zero-range attractive potential in one dimension, and fermion-fermion interacting via an attractive Gaussian potential in three dimensions. The results presented here provide a guide for improving the adiabatic projection method and for reducing the computational costs of large-scale calculations of \emph{ab initio} nuclear scattering and reactions using Monte Carlo methods.


title: A constraint for twist equivalence of cusp forms on GL$(n)$

abstract: This Note answers, and generalizes, a question of Kaisa Matom\"aki. We show that give two cuspidal automorphic representations $\pi_1$ and $\pi_2$ of $GL_n$ over a number field $F$ of respective conductors $N_1,$ $N_2,$ every character $\chi$ such that $\pi_1\otimes\chi\simeq\pi_2$ of conductor $Q,$ satisfies the bound: $Q^n\mid N_1N_2.$ If at every finite place $v,$ $\pi_{1,v}$ is a discrete series whenever it is ramified, then $Q^n$ divides the least common multiple $[N_1, N_2].$


title: Understanding nebular spectra of Type Ia supernovae

abstract: In this study, we present one-dimensional, non-local-thermodynamic-equilibrium, radiative transfer simulations (using CMFGEN) in which we introduce micro-clumping at nebular times into two Type Ia supernova ejecta models. We use one sub-Chandrasekhar (sub-M$_{\rm Ch}$) ejecta with 1.02 M$_\odot$ and one M$_{\rm Ch}$ ejecta model with 1.40 M$_\odot$. We introduce clumping factors $f=$0.33,0.25, and 0.10 which are constant throughout the ejecta and compared to the unclumped $f=1.0$ case. We find that clumping is a natural mechanism to reduce the ionization of the ejecta, reducing emission from [Fe III], [Ar III], and [S III] by a factor of a few. For decreasing values of the clumping factor $f$, the [Ca II] $\lambda\lambda$7291,7324 doublet became a dominant cooling line for our M$_{\rm Ch}$ model but still weak in our sub-M$_{\rm Ch}$ model. Strong [Ca II] $\lambda\lambda$7291,7324 indicates non-thermal heating in that region and may constrain explosion modelling. Due to the low abundance of stable nickel, our sub-M$_{\rm Ch}$ model never showed the [Ni II] 1.939 micron diagnostic feature for all clumping values.


title: Resolving Overlapping Convex Objects in Silhouette Images by Concavity Analysis and Gaussian Process

abstract: Segmentation of overlapping convex objects has various applications, for example, in nanoparticles and cell imaging. Often the segmentation method has to rely purely on edges between the background and foreground making the analyzed images essentially silhouette images. Therefore, to segment the objects, the method needs to be able to resolve the overlaps between multiple objects by utilizing prior information about the shape of the objects. This paper introduces a novel method for segmentation of clustered partially overlapping convex objects in silhouette images. The proposed method involves three main steps: pre-processing, contour evidence extraction, and contour estimation. Contour evidence extraction starts by recovering contour segments from a binarized image by detecting concave points. After this, the contour segments which belong to the same objects are grouped. The grouping is formulated as a combinatorial optimization problem and solved using the branch and bound algorithm. Finally, the full contours of the objects are estimated by a Gaussian process regression method. The experiments on a challenging dataset consisting of nanoparticles demonstrate that the proposed method outperforms three current state-of-art approaches in overlapping convex objects segmentation. The method relies only on edge information and can be applied to any segmentation problems where the objects are partially overlapping and have a convex shape.


title: Cores and Other Dense Structures in Complex Networks

abstract: Complex networks are a powerful paradigm to model complex systems. Specific network models, e.g., multilayer networks, temporal networks, and signed networks, enrich the standard network representation with additional information to better capture real-world phenomena. Despite the keen interest in a variety of problems, algorithms, and analysis methods for these types of network, the problem of extracting cores and dense structures still has unexplored facets. In this work, we present advancements to the state of the art by the introduction of novel definitions and algorithms for the extraction of dense structures from complex networks, mainly cores. At first, we define core decomposition in multilayer networks together with a series of applications built on top of it, i.e., the extraction of maximal multilayer cores only, densest subgraph in multilayer networks, the speed-up of the extraction of frequent cross-graph quasi-cliques, and the generalization of community search to the multilayer setting. Then, we introduce the concept of core decomposition in temporal networks; also in this case, we are interested in the extraction of maximal temporal cores only. Finally, in the context of discovering polarization in large-scale online data, we study the problem of identifying polarized communities in signed networks. The proposed methodologies are evaluated on a large variety of real-world networks against na\"{\i}ve approaches, non-trivial baselines, and competing methods. In all cases, they show effectiveness, efficiency, and scalability. Moreover, we showcase the usefulness of our definitions in concrete applications and case studies, i.e., the temporal analysis of contact networks, and the identification of polarization in debate networks.


title: Quantitative Propagation of Chaos in the bimolecular chemical reaction-diffusion model

abstract: We study a stochastic system of $N$ interacting particles which models bimolecular chemical reaction-diffusion. In this model, each particle $i$ carries two attributes: the spatial location $X_t^i\in \mathbb{T}^d$, and the type $\Xi_t^i\in \{1,\cdots,n\}$. While $X_t^i$ is a standard (independent) diffusion process, the evolution of the type $\Xi_t^i$ is described by pairwise interactions between different particles under a series of chemical reactions described by a chemical reaction network. We prove that in the large particle limit the stochastic dynamics converges to a mean field limit which is described by a nonlocal reaction-diffusion partial differential equation. In particular, we obtain a quantitative propagation of chaos result for the interacting particle system. Our proof is based on the relative entropy method used recently by Jabin and Wang \cite{JW18}. The key ingredient of the relative entropy method is a large deviation estimate for a special partition function, which was proved previously by technical combinatorial estimates. We give a simple probabilistic proof based on a novel martingale argument.


title: Galaxy rotation curves and preferred reference frame effects

abstract: As an alternative to dark matter models we use generalized Jordan-Brans-Dicke scalar-vector-tensor (JBD-SVT) gravity model to study the behavior of the rotational velocities of test particles moving around galaxies. To do so we consider an interaction potential $U(\phi,N_{\mu})$ between the Brans-Dicke scalar field $\phi$ and time like dynamical four-vector field $N_\mu$ which plays as four velocity of a preferred reference frame. We show that at in weak field limits metric solution of the galaxy under consideration reaches to a modified Schwarzschild-de Sitter space in which mass of the vector field plays as an effective cosmological constant. In fact the present work proposes modification on the formulation of Newton's gravitational acceleration. This is used to explain circular velocity of galaxies without postulating dark matter. We also check our theoretical results with empirical baryonic Tully Fisher relation which states a linear relations between the rotational speed of galaxies and their mass. Mathematical calculation predict a good correspondence between our theoretical results and experimental observations for a set of 12 spiral galaxies.


title: Multi-time distribution in discrete polynuclear growth

abstract: We study the multi-time distribution in a discrete polynuclear growth model or, equivalently, in directed last-passage percolation with geometric weights. A formula for the joint multi-time distribution function is derived in the discrete setting. It takes the form of a multiple contour integral of a block Fredholm determinant. The asymptotic multi-time distribution is then computed by taking the appropriate KPZ-scaling limit of this formula. This distribution is expected to be universal for models in the Kardar-Parisi-Zhang universality class.

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


title: Deep 3D Convolutional Neural Network for Automated Lung Cancer Diagnosis

abstract: Computer Aided Diagnosis has emerged as an indispensible technique for validating the opinion of radiologists in CT interpretation. This paper presents a deep 3D Convolutional Neural Network (CNN) architecture for automated CT scan-based lung cancer detection system. It utilizes three dimensional spatial information to learn highly discriminative 3 dimensional features instead of 2D features like texture or geometric shape whick need to be generated manually. The proposed deep learning method automatically extracts the 3D features on the basis of spatio-temporal statistics.The developed model is end-to-end and is able to predict malignancy of each voxel for given input scan. Simulation results demonstrate the effectiveness of proposed 3D CNN network for classification of lung nodule in-spite of limited computational capabilities.


title: A Systemic IoT-Fog-Cloud Architecture for Big-Data Analytics and Cyber Security Systems: A Review of Fog Computing

abstract: Abstract--- With the rapid growth of the Internet of Things (IoT), current Cloud systems face various drawbacks such as lack of mobility support, location-awareness, geo-distribution, high latency, as well as cyber threats. Fog/Edge computing has been proposed for addressing some of the drawbacks, as it enables computing resources at the network's edges and it locally offers big-data analytics rather than transmitting them to the Cloud. The Fog is defined as a Cloud-like system having similar functions, including software-, platform- and infrastructure-as services. The deployment of Fog applications faces various security issues related to virtualisation, network monitoring, data protection and attack detection. This paper proposes a systemic IoT-Fog-Cloud architecture that clarifies the interactions between the three layers of IoT, Fog and Cloud for effectively implementing big-data analytics and cyber security applications. It also reviews security challenges, solutions and future research directions in the architecture.


title: A separator-based method for generating weakly chordal graphs

abstract: We propose a scheme for generating a weakly chordal graph on n vertices with m edges. In this method, we first construct a tree and then generate an orthogonal layout (which is a weakly chordal graph on the n vertices) based on this tree. In the next and final step, we insert additional edges to give us a weakly chordal graph on m edges. Our algorithm ensures that the graph remains weakly chordal after each edge is inserted. The time complexity of an insertion query is O(n^3) time and an insertion takes constant time. On the other hand, a generation algorithm based on finding a 2-pair takes O(nm) time using the algorithm of Arikati and Rangan [1].


title: Semi-parametric Bayesian variable selection for gene-environment interactions

abstract: Many complex diseases are known to be affected by the interactions between genetic variants and environmental exposures beyond the main genetic and environmental effects. Study of gene-environment (G$\times$E) interactions is important for elucidating the disease etiology. Existing Bayesian methods for G$\times$E interaction studies are challenged by the high-dimensional nature of the study and the complexity of environmental influences. Many studies have shown the advantages of penalization methods in detecting G$\times$E interactions in "large p, small n" settings. However, Bayesian variable selection, which can provide fresh insight into G$\times$E study, has not been widely examined. We propose a novel and powerful semi-parametric Bayesian variable selection model that can investigate linear and nonlinear G$\times$E interactions simultaneously. Furthermore, the proposed method can conduct structural identification by distinguishing nonlinear interactions from main-effects-only case within the Bayesian framework. Spike and slab priors are incorporated on both individual and group levels to identify the sparse main and interaction effects. The proposed method conducts Bayesian variable selection more efficiently than existing methods. Simulation shows that the proposed model outperforms competing alternatives in terms of both identification and prediction. The proposed Bayesian method leads to the identification of main and interaction effects with important implications in a high-throughput profiling study with high-dimensional SNP data.


title: Spinor-helicity variables for cosmological horizons in de Sitter space

abstract: We consider massless fields of arbitrary spin in de Sitter space. We introduce a spinor-helicity formalism, which encodes the field data on a cosmological horizon. These variables reduce the free S-matrix in an observer's causal patch, i.e. the evolution of free fields from one horizon to another, to a simple Fourier transform. We show how this result arises via twistor theory, by decomposing the horizon<->horizon problem into a pair of (more symmetric) horizon<->twistor problems.


title: An entire function connected with the approximation of the golden ratio

abstract: In 1987, R. B. Paris uses the analytic function \[\label{main} g(w)=\lim_{n\to\infty}(2\varphi)^n\biggl(\underbrace{\sqrt{1+\sqrt{1+...\sqrt{1+w}}}}_n-\varphi\biggr),\ \ \ \varphi=\frac{1+\sqrt{5}}2, \] to estimate the convergence of nested squares to the golden ratio. The function $g$ is non-entire and, perhaps, can not be expressed in terms of some standard known functions. We show that $f(z):=g^{-1}(z)$ is an entire function satisfying Poincare equality. While $f$ has zeros of various multiplicities, it can be expressed in terms of its simple zeros, forming fractal structures similar to Julia sets.


title: Robust stability of moving horizon estimation for nonlinear systems with bounded disturbances using adaptive arrival cost

abstract: In this paper, the robust stability and convergence to the true state of moving horizon estimator based on an adaptive arrival cost are established for nonlinear detectable systems. Robust global asymptotic stability is shown for the case of non-vanishing bounded disturbances whereas the convergence to the true state is proved for the case of vanishing disturbances. Several simulations were made in order to show the estimator behaviour under different operational conditions and to compare it with the state of the art estimation methods.


title: Localization Requirements for Autonomous Vehicles

abstract: Autonomous vehicles require precise knowledge of their position and orientation in all weather and traffic conditions for path planning, perception, control, and general safe operation. Here we derive these requirements for autonomous vehicles based on first principles. We begin with the safety integrity level, defining the allowable probability of failure per hour of operation based on desired improvements on road safety today. This draws comparisons with the localization integrity levels required in aviation and rail where similar numbers are derived at 10^-8 probability of failure per hour of operation. We then define the geometry of the problem, where the aim is to maintain knowledge that the vehicle is within its lane and to determine what road level it is on. Longitudinal, lateral, and vertical localization error bounds (alert limits) and 95% accuracy requirements are derived based on US road geometry standards (lane width, curvature, and vertical clearance) and allowable vehicle dimensions. For passenger vehicles operating on freeway roads, the result is a required lateral error bound of 0.57 m (0.20 m, 95%), a longitudinal bound of 1.40 m (0.48 m, 95%), a vertical bound of 1.30 m (0.43 m, 95%), and an attitude bound in each direction of 1.50 deg (0.51 deg, 95%). On local streets, the road geometry makes requirements more stringent where lateral and longitudinal error bounds of 0.29 m (0.10 m, 95%) are needed with an orientation requirement of 0.50 deg (0.17 deg, 95%).


title: A Curated Image Parameter Dataset from Solar Dynamics Observatory Mission

abstract: We provide a large image parameter dataset extracted from the Solar Dynamics Observatory (SDO) mission's AIA instrument, for the period of January 2011 through the current date, with the cadence of six minutes, for nine wavelength channels. The volume of the dataset for each year is just short of 1 TiB. Towards achieving better results in the region classification of active regions and coronal holes, we improve upon the performance of a set of ten image parameters, through an in depth evaluation of various assumptions that are necessary for calculation of these image parameters. Then, where possible, a method for finding an appropriate settings for the parameter calculations was devised, as well as a validation task to show our improved results. In addition, we include comparisons of JP2 and FITS image formats using supervised classification models, by tuning the parameters specific to the format of the images from which they are extracted, and specific to each wavelength. The results of these comparisons show that utilizing JP2 images, which are significantly smaller files, is not detrimental to the region classification task that these parameters were originally intended for. Finally, we compute the tuned parameters on the AIA images and provide a public API (http://dmlab.cs.gsu.edu/dmlabapi) to access the dataset. This dataset can be used in a range of studies on AIA images, such as content-based image retrieval or tracking of solar events, where dimensionality reduction on the images is necessary for feasibility of the tasks.


title: Constraining R-parity-violating couplings in $\tau$-processes at the LHC and in electroweak precision measurements

abstract: We find new limits on the $\lambda'_{3jk}$ R-parity violating (RPV) couplings of the minimal supersymmetric standard model, using Drell-Yan differential cross sections at the LHC and electroweak precision measurements from LEP and SLC. Specifically, limits on six out of the nine $\lambda'_{3jk}$-couplings, with $j=1\text{ or }2$, are obtained using Drell-Yan data, with the remaining three (for $j=3$) bounded by precision electroweak data. We also update the limits on all $\lambda'_{ijk}$-couplings using electroweak data and find new bounds on $\lambda'_{132}$ and $\lambda'_{232}$ that are stronger than obtained elsewhere. A table of all current bounds on $\lambda'_{ijk}$ is given in an appendix.


title: Photometric observations of the asteroid 3200 Phaethon using small and middle telescopes

abstract: The main aim of photometrical observations of the asteroid 3200 Phaethon was searching for its low-level cometary activity (possible coma and/or dust tail) in the pre-perihelion passage. We performed observational runs with telescopes ranging from 0.61-m to 2-m and BVR color imaging. Three longer photometric series were used for modeling of the 3D shape of Phaethon. The color indices and size of the asteroid were estimated.


title: Partial Phase Cohesiveness in Networks of Kuramoto Oscillator Networks

abstract: Partial, instead of complete, synchronization has been widely observed in various networks including, in particular, brain networks. Motivated by data from human brain functional networks, in this technical note, we analytically show that partial synchronization can be induced by strong regional connections in coupled subnetworks of Kuramoto oscillators. To quantify the required strength of regional connections, we first obtain a critical value for the algebraic connectivity of the corresponding subnetwork using the incremental 2 norm. We then introduce the concept of the generalized complement graph, and obtain another condition on the weighted nodal degree by using the incremental infinity norm. Under these two conditions, regions of attraction for partial phase cohesiveness are estimated in the forms of the incremental 2 and infinity norms, respectively. Our result based on the incremental infinity norm is the first known criterion that is applicable to non-complete graphs. Numerical simulations are performed on a two-level network to illustrate our theoretical results; more importantly, we use real anatomical brain network data to show how our results may reveal the interplay between anatomical structure and empirical patterns of synchrony.


title: Clues on the Origin and Evolution of Massive Contact Binaries: Atmosphere Analysis of VFTS 352

abstract: The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment.


title: Dynamics of geodesics, and Maass cusp forms

abstract: The correspondence principle in physics between quantum mechanics and classical mechanics suggests deep relations between spectral and geometric entities of Riemannian manifolds. We survey---in a way intended to be accessible to a wide audience of mathematicians---a mathematically rigorous instance of such a relation that emerged in recent years, showing a dynamical interpretation of certain Laplace eigenfunctions of hyperbolic surfaces.


title: Layer-dependent spin-orbit torques generated by the centrosymmetric transition metal dichalcogenide $\beta$-MoTe$_2$

abstract: Single-crystal materials with sufficiently low crystal symmetry and strong spin-orbit interactions can be used to generate novel forms of spin-orbit torques on adjacent ferromagnets, such as the out-of-plane antidamping torque previously observed in WTe$_2$/ferromagnet heterostructures. Here, we present measurements of spin-orbit torques produced by the low-symmetry material $\beta$-MoTe$_2$, which unlike WTe$_2$ retains bulk inversion symmetry. We measure spin-orbit torques on $\beta$-MoTe$_2$/Permalloy heterostructures using spin-torque ferromagnetic resonance as a function of crystallographic alignment and MoTe$_2$ thickness down to the monolayer limit. We observe an out-of-plane antidamping torque with a spin torque conductivity as strong as 1/3 of that of WTe$_2$, demonstrating that the breaking of bulk inversion symmetry in the spin-generation material is not a necessary requirement for producing an out-of-plane antidamping torque. We also measure an unexpected dependence on the thickness of the $\beta$-MoTe$_2$ -- the out-of-plane antidamping torque is present in MoTe$_2$/Permalloy heterostructures when the $\beta$-MoTe$_2$ is a monolayer or trilayer thick, but goes to zero for devices with bilayer $\beta$-MoTe$_2$.


title: Selling Demand Response Using Options

abstract: Wholesale electricity markets in many jurisdictions use a two-settlement structure: a day-ahead market for bulk power transactions and a real-time market for fine-grain supply-demand balancing. This paper explores trading demand response assets within this two-settlement market structure. We consider two approaches for trading demand response assets: (a) an intermediate spot market with contingent pricing, and (b) an over-the-counter options contract. In the first case, we characterize the competitive equilibrium of the spot market, and show that it is socially optimal. Economic orthodoxy advocates spot markets, but these require expensive infrastructure and regulatory blessing. In the second case, we characterize competitive equilibria and compare its efficiency with the idealized spot market. Options contract are private bilateral over-the-counter transactions and do not require regulatory approval. We show that the optimal social welfare is, in general, not supported. We then design optimal option prices that minimize the social welfare gap. This optimal design serves to approximate the ideal spot market for demand response using options with modest loss of efficiency. Our results are validated through numerical simulations.


title: Attracting and repelling 2-body problems on a family of surfaces of constant curvature

abstract: We first provide a classification of the pure rotational motion of 2 particles on a sphere interacting via a repelling potential. This is achieved by providing a simple geometric equivalence between repelling particles and attracting particles, and relying on previous work on the similar classification for attracting particles. The second theme of the paper is to study the 2-body problem on a surface of constant curvature treating the curvature as a parameter, and with particular interest in how families of relative equilibria and their stability behave as the curvature passes through zero and changes sign. We consider two cases: firstly one where the particles are always attracting throughout the family, and secondly where they are attracting for negative curvature and repelling for positive curvature, interpolated by no interaction when the curvature vanishes. Our analysis clarifies the role of curvature in the existence and stability of relative equilibria.


title: Sea of Genes: Combining Animation and Narrative Strategies to Visualize Metagenomic Data for Museums

abstract: We examine the application of narrative strategies to present a complex and unfamiliar metagenomics dataset to the public in a science museum. Our dataset contains information about microbial gene expressions that scientists use to infer the behavior of microbes. This exhibit had three goals: to inform (the) public about microbes' behavior, cycles, and patterns; to link their behavior to the concept of gene expression; and to highlight scientists' use of gene expression data to understand the role of microbes. To address these three goals, we created a visualization with three narrative layers, each layer corresponding to a goal. This study presented us with an opportunity to assess existing frameworks for narrative visualization in a naturalistic setting. We present three successive rounds of design and evaluation of our attempts to engage visitors with complex data through narrative visualization. We highlight our design choices and their underlying rationale based on extant theories. We conclude that a central animation based on a curated dataset could successfully achieve our first goal, i.e., to communicate the aggregate behavior and interactions of microbes. We failed to achieve our second goal and had limited success with the third goal. Overall, this study highlights the challenges of telling multi-layered stories and the need for new frameworks for communicating layered stories in public settings.


title: The host galaxies of double compact objects across cosmic time

abstract: We explore the host galaxies of compact-object binaries (black hole--black hole binaries, BHBs; neutron star--black hole binaries, NSBHs; double-neutron stars; DNSs) across cosmic time, by means of population-synthesis simulations combined with the Illustris cosmological simulation. At high redshift ($z\gtrsim{}4$) the host galaxies of BHBs, NSBHs and DNSs are very similar and are predominantly low-mass galaxies (stellar mass $M<10^{11}$ M$_\odot$). If $z\gtrsim{}4$ most compact objects form and merge in the same galaxy, with a short delay time. At low redshift ($z\leq{}2$), the host galaxy populations of DNSs differ significantly from the host galaxies of both BHBs and NSBHs. DNSs merging at low redshift tend to form and merge in the same galaxy, with relatively short delay time. The stellar mass of DNS hosts peaks around $\sim{}10^{10}-10^{11}$ M$_\odot$. In contrast, BHBs and NSBHs merging at low redshift tend to form in rather small galaxies at high redshift and then to merge in larger galaxies with long delay times. This difference between DNSs and black hole binaries is a consequence of their profoundly different metallicity dependence.


title: Three waves of chemical dynamics

abstract: Three epochs in development of chemical dynamics are presented. We try to understand the modern research programs in the light of classical works. Three eras (or waves) of chemical dynamics can be revealed in the flux of research and publications. These waves may be associated with leaders: the first is the van't Hoof wave, the second may be called the Semenov--Hinshelwood wave and the third is definitely the Aris wave. Of course, the whole building was impossible without efforts of hundreds of other researchers. Some of of them are mentioned in our brief review.


title: Bruggeman Approach for Isotropic Chiral Mixtures Revisited

abstract: Two interpretations of the Bruggeman approach for the homogenization of isotropic chiral mixtures are shown to lead to different results. Whereas the standard interpretation is shown to yield the average-polarizability-density approach, a recent interpretation turns out to deliver a null-excess-polarization approach. The difference between the two interpretations arises from differing treatments of the local field


title: Using your ADC and Backend Capacitors to Authenticate for Free: (Virtually) Free from Database, Enrollment, and Excessive Area/Power Overhead

abstract: Detection of counterfeit chips has emerged as a crucial concern. Physically-unclonable-function (PUF)-based techniques are widely used for authentication, however, require dedicated hardware and large signature database. In this work, we show intrinsic & database-free authentication using back-end capacitors. The discussed technique simplifies authentication setup and reduces the test cost. We show that an analog-to-digital converter (ADC) can be modified for back-end capacitor-based authentication in addition to its regular functionality; hence, a dedicated authentication module is not necessary. Moreover, since back-end capacitors are quite insensitive to temperature and aging-induced variations than transistors, the discussed technique result in a more reliable authentication than transistor PUF-based authentication. The modifications to conventional ADC incur 3.2% power overhead and 75% active-area overhead; however, arguably, the advantages of the discussed intrinsic & database-free authentication outweigh the overheads. Full version of this article is published at IEEE TVLSI.


title: Episodic Memory in Lifelong Language Learning

abstract: We introduce a lifelong language learning setup where a model needs to learn from a stream of text examples without any dataset identifier. We propose an episodic memory model that performs sparse experience replay and local adaptation to mitigate catastrophic forgetting in this setup. Experiments on text classification and question answering demonstrate the complementary benefits of sparse experience replay and local adaptation to allow the model to continuously learn from new datasets. We also show that the space complexity of the episodic memory module can be reduced significantly (~50-90%) by randomly choosing which examples to store in memory with a minimal decrease in performance. We consider an episodic memory component as a crucial building block of general linguistic intelligence and see our model as a first step in that direction.


title: Discussion about Different Methods for Introducing the Turbulent Boundary Layer Excitation in Vibroacoustic Models

abstract: For controlling the noise radiated from vibrating structures excited by turbulent boundary layer (TBL) it is relevant to develop numerical tools for understanding how the structure reacts to TBL excitation. Usually, the wall pressure fluctuations of the TBL are described through statistical quantities (i.e. space-frequency or wavenumber-frequency spectra) which depend on the TBL parameters. On the other hand, the vibro-acoustic models (i.e. Finite Elements, Boundary Elements, Transfer Matrix Methods, Analytical models, etc) evaluate deterministic transfer functions which characterise the response of the considered structures. The first part of this paper focuses on the coupling between the stochastic TBL and the deterministic vibro-acoustic models. Five techniques are presented. Numerical applications on an academic marine test case are proposed in order to discuss the calculation parameters and the interests / drawbacks of each technique. In the second part of the paper, the high frequency modelling with the Statistical Energy Analysis (SEA) method is considered. The focus is placed on the estimation of an important input of this method: the injected power by the TBL into the structure for each third octave band.

categories: physics.app-ph math.NA physics.class-ph physics.comp-ph physics.flu-dyn


title: Increasing Compactness Of Deep Learning Based Speech Enhancement Models With Parameter Pruning And Quantization Techniques

abstract: Most recent studies on deep learning based speech enhancement (SE) focused on improving denoising performance. However, successful SE applications require striking a desirable balance between denoising performance and computational cost in real scenarios. In this study, we propose a novel parameter pruning (PP) technique, which removes redundant channels in a neural network. In addition, a parameter quantization (PQ) technique was applied to reduce the size of a neural network by representing weights with fewer cluster centroids. Because the techniques are derived based on different concepts, the PP and PQ can be integrated to provide even more compact SE models. The experimental results show that the PP and PQ techniques produce a compacted SE model with a size of only 10.03% compared to that of the original model, resulting in minor performance losses of 1.43% (from 0.70 to 0.69) for STOI and 3.24% (from 1.85 to 1.79) for PESQ. The promising results suggest that the PP and PQ techniques can be used in a SE system in devices with limited storage and computation resources.


title: Single T gate in a Clifford circuit drives transition to universal entanglement spectrum statistics

abstract: Clifford circuits are insufficient for universal quantum computation or creating $t$-designs with $t\ge 4$. While the entanglement entropy is not a telltale of this insufficiency, the entanglement spectrum is: the entanglement levels are Poisson-distributed for circuits restricted to the Clifford gate-set, while the levels follow Wigner-Dyson statistics when universal gates are used. In this paper we show, using finite-size scaling analysis of different measures of level spacing statistics, that in the thermodynamic limit, inserting a single T $(\pi/8)$ gate in the middle of a random Clifford circuit is sufficient to alter the entanglement spectrum from a Poisson to a Wigner-Dyson distribution.


title: Building a Maxey--Riley framework for surface ocean inertial particle dynamics

abstract: A framework for the study of surface ocean inertial particle motion is built from the Maxey--Riley set. A new set is obtained by vertically averaging each term of the original set, adapted to account for Earth's rotation effects, across the extent of a sufficiently small spherical particle that floats at an assumed unperturbed air--sea interface with unsteady nonuniform winds and ocean currents above and below, respectively. The inertial particle velocity is shown to exponentially decay in time to a velocity that lies close to an average of seawater and air velocities, weighted by a function of the seawater-to-particle density ratio. Such a weighted average velocity turns out to fortuitously be of the type commonly discussed in the search-and-rescue literature, which alone cannot explain the observed role of anticyclonic mesoscale eddies as traps for marine debris or the formation of great garbage patches in the subtropical gyres, phenomena dominated by finite-size effects. A heuristic extension of the theory is proposed to describe the motion of nonspherical particles by means of a simple shape factor correction, and recommendations are made for incorporating wave-induced Stokes drift, and allowing for inhomogeneities of the carrying fluid density. The new Maxey--Riley set outperforms an ocean adaptation that ignored wind drag effects and the first reported adaption that attempted to incorporate them.


title: Handling Divergent Reference Texts when Evaluating Table-to-Text Generation

abstract: Automatically constructed datasets for generating text from semi-structured data (tables), such as WikiBio, often contain reference texts that diverge from the information in the corresponding semi-structured data. We show that metrics which rely solely on the reference texts, such as BLEU and ROUGE, show poor correlation with human judgments when those references diverge. We propose a new metric, PARENT, which aligns n-grams from the reference and generated texts to the semi-structured data before computing their precision and recall. Through a large scale human evaluation study of table-to-text models for WikiBio, we show that PARENT correlates with human judgments better than existing text generation metrics. We also adapt and evaluate the information extraction based evaluation proposed by Wiseman et al (2017), and show that PARENT has comparable correlation to it, while being easier to use. We show that PARENT is also applicable when the reference texts are elicited from humans using the data from the WebNLG challenge.


title: Representation Theoretic Patterns in Multi-Frequency Class Averaging for Three-Dimensional Cryo-Electron Microscopy

abstract: We develop in this paper a novel intrinsic classification algorithm -- multi-frequency class averaging (MFCA) -- for classifying noisy projection images obtained from three-dimensional cryo-electron microscopy (cryo-EM) by the similarity among their viewing directions. This new algorithm leverages multiple irreducible representations of the unitary group to introduce additional redundancy into the representation of the optimal in-plane rotational alignment, extending and outperforming the existing class averaging algorithm that uses only a single representation. The formal algebraic model and representation theoretic patterns of the proposed MFCA algorithm extend the framework of Hadani and Singer to arbitrary irreducible representations of the unitary group. We conceptually establish the consistency and stability of MFCA by inspecting the spectral properties of a generalized local parallel transport operator through the lens of Wigner $D$-matrices. We demonstrate the efficacy of the proposed algorithm with numerical experiments.


title: MelNet: A Generative Model for Audio in the Frequency Domain

abstract: Capturing high-level structure in audio waveforms is challenging because a single second of audio spans tens of thousands of timesteps. While long-range dependencies are difficult to model directly in the time domain, we show that they can be more tractably modelled in two-dimensional time-frequency representations such as spectrograms. By leveraging this representational advantage, in conjunction with a highly expressive probabilistic model and a multiscale generation procedure, we design a model capable of generating high-fidelity audio samples which capture structure at timescales that time-domain models have yet to achieve. We apply our model to a variety of audio generation tasks, including unconditional speech generation, music generation, and text-to-speech synthesis---showing improvements over previous approaches in both density estimates and human judgments.


title: First Results on the Search for Chameleons with the KWISP Detector at CAST

abstract: We report on a first measurement with a sensitive opto-mechanical force sensor designed for the direct detection of coupling of real chameleons to matter. These dark energy candidates could be produced in the Sun and stream unimpeded to Earth. The KWISP detector installed on the CAST axion search experiment at CERN looks for tiny displacements of a thin membrane caused by the mechanical effect of solar chameleons. The displacements are detected by a Michelson interferometer with a homodyne readout scheme. The sensor benefits from the focusing action of the ABRIXAS X-ray telescope installed at CAST, which increases the chameleon flux on the membrane. A mechanical chopper placed between the telescope output and the detector modulates the incoming chameleon stream. We present the results of the solar chameleon measurements taken at CAST in July 2017, setting an upper bound on the force acting on the membrane of $80$~pN at 95\% confidence level. The detector is sensitive for direct coupling to matter $10^4 \leq\beta_m \leq 10^8$, where the coupling to photons is locally bound to $\beta_\gamma \leq 10^{11}$.

categories: hep-ex astro-ph.CO astro-ph.SR physics.ins-det


title: Unavoidable disorder and entropy in multi-component systems

abstract: The need for improved functionalities is driving the search for more complicated multi-component materials. Despite the factorially increasing composition space, ordered compounds with 4 or more species are rare. Here, we unveil the competition between the gain in enthalpy and entropy with increasing number of species by statistical analysis of the AFLOW data repositories. A threshold in the number of species is found where entropy gain exceeds enthalpy gain. Beyond that, enthalpy can be neglected, and disorder - complete or partial - is unavoidable.


title: Coherent electrical control of a single high-spin nucleus in silicon

abstract: Nuclear spins are highly coherent quantum objects. In large ensembles, their control and detection via magnetic resonance is widely exploited, e.g. in chemistry, medicine, materials science and mining. Nuclear spins also featured in early ideas and demonstrations of quantum information processing. Scaling up these ideas requires controlling individual nuclei, which can be detected when coupled to an electron. However, the need to address the nuclei via oscillating magnetic fields complicates their integration in multi-spin nanoscale devices, because the field cannot be localized or screened. Control via electric fields would resolve this problem, but previous methods relied upon transducing electric signals into magnetic fields via the electron-nuclear hyperfine interaction, which severely affects the nuclear coherence. Here we demonstrate the coherent quantum control of a single antimony (spin-7/2) nucleus, using localized electric fields produced within a silicon nanoelectronic device. The method exploits an idea first proposed in 1961 but never realized experimentally with a single nucleus. Our results are quantitatively supported by a microscopic theoretical model that reveals how the purely electrical modulation of the nuclear electric quadrupole interaction, in the presence of lattice strain, results in coherent nuclear spin transitions. The spin dephasing time, 0.1 seconds, surpasses by orders of magnitude those obtained via methods that require a coupled electron spin for electrical drive. These results show that high-spin quadrupolar nuclei could be deployed as chaotic models, strain sensors and hybrid spin-mechanical quantum systems using all-electrical controls. Integrating electrically controllable nuclei with quantum dots could pave the way to scalable nuclear- and electron-spin-based quantum computers in silicon that operate without the need for oscillating magnetic fields.


title: Graph Sampling for Matrix Completion Using Recurrent Gershgorin Disc Shift

abstract: Matrix completion algorithms fill missing entries in a large matrix given a subset of observed samples. However, how to best pre-select informative matrix entries given a sampling budget is largely unaddressed. In this paper, we propose a fast sample selection strategy for matrix completion from a graph signal processing perspective. Specifically, we first regularize the matrix reconstruction objective using a dual graph signal smoothness prior, resulting in a system of linear equations for solution. We then select appropriate samples to maximize the smallest eigenvalue $\lambda_{\min}$ of the coefficient matrix, thus maximizing the stability of the linear system. To efficiently solve this combinatorial problem, we derive a greedy sampling strategy, leveraging on Gershgorin circle theorem, that iteratively selects one sample (equivalent to shifting one Gershgorin disc) at a time corresponding to the largest magnitude entry in the first eigenvector of a modified graph Laplacian matrix. Our algorithm benefits computationally from warm start as the first eigenvectors of incremented Laplacian matrices are computed recurrently for more samples. To achieve computation scalability when sampling large matrices, we further rewrite the coefficient matrix as a sum of two separate components, each of which exhibits block-diagonal structure that we exploit for alternating block-wise sampling. Extensive experiments on both synthetic and real-world datasets show that our graph sampling algorithm substantially outperforms existing sampling schemes for matrix completion and reduces the completion error, when combined with a range of modern matrix completion algorithms.


title: hankel: A Python library for performing simple and accurate Hankel transformations

abstract: This paper presents \textsc{hankel}, a pure-python code for solving Hankel-type integrals and transforms. Such transforms are common in the physical sciences, especially appearing as the radial solution to angularly symmetric Fourier Transforms in arbitrary dimensions. The code harnesses the advantages of solving such transforms via the one-dimensional Hankel transform -- an increase in conceptual simplicity and efficiency -- and implements them in the user-friendly and flexible Python language. We discuss several limitations of the adopted method, and point to the code's extensive documentation for further examples.


title: Screen-printed and spray coated graphene-based RFID transponders

abstract: We report Ultra-High-Frequency (UHF, 800MHz-1GHz) Radio Frequency Identification (RFID) transponders consisting of printed dipole antennas combined with RFID microchips. These are fabricated on Kapton via screen printing and on paper via spray coating, using inks obtained via microfluidization of graphite. We introduce a hybrid antenna structure, where an Al loop (small compared to the overall size of the antenna) is connected to a microchip with the double function of matching the impedances of antenna and microchip and avoiding bonding between exfoliated graphite and chip. The transponders have reading distance~11m at UHF RFID frequencies, larger than previously reported for graphene-based RFID and comparable with commercial transponders based on metallic antennas.


title: Well-posedness for Photoacoustic Tomography with Fabry-Perot Sensors

abstract: In the mathematical analysis of photoacoustic imaging, it is usually assumed that the acoustic pressure (Dirichlet data) is measured on a detection surface. However, actual ultrasound detectors gather data of a different type. In this paper, we propose a more realistic mathematical model of ultrasound measurements acquired by the Fabry--Perot sensor. This modeling incorporates directional response of such sensors. We study the solvability of the resulting photoacoustic tomography problem, concluding that the problem is well-posed under certain assumptions. Numerical reconstructions are implemented using the Landweber iterations, after discretization of the governing equations using the finite element method.


title: Molecule-surface interaction from van der Waals-corrected semilocal density functionals: the example of thiophene on transition-metal surfaces

abstract: Semi-local density functional approximations are widely used. None of them can capture the long-range van der Waals (vdW) attraction between separated subsystems, but they differ remarkably in the extent to which they capture intermediate-range vdW effects responsible for equilibrium bonds between neighboring small closed-shell subsystems. The local density approximation (LDA) often overestimates this effect, while the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) underestimates it. The strongly-constrained and appropriately normed (SCAN) meta-GGA often estimates it well. All of these semi-local functionals require an additive non-local correction such as the revised Vydrov-Van Voorhis 2010 (rVV10) to capture the long-range part. This work reports adsorption energies and the corresponding geometry of the aromatic thiophene (C$_4$H$_4$S) bound to transition metal surfaces. The adsorption process requires a genuine interplay of covalent and weak binding and requires a simultaneously accurate description of surface and adsorption energies with the correct prediction of the adsorption site. All these quantities must come from well balanced short and long-range correlation effects for a universally applicable method for weak interactions with chemical accuracy.


title: Ramsey, Paper, Scissors

abstract: We introduce a graph Ramsey game called Ramsey, Paper, Scissors. This game has two players, Proposer and Decider. Starting from an empty graph on $n$ vertices, on each turn Proposer proposes a potential edge and Decider simultaneously decides (without knowing Proposer's choice) whether to add it to the graph. Proposer cannot propose an edge which would create a triangle in the graph. The game ends when Proposer has no legal moves remaining, and Proposer wins if the final graph has independence number at least $s$. We prove a threshold phenomenon exists for this game by exhibiting randomized strategies for both players that are optimal up to constants. Namely, there exist constants $0<A<B$ such that (under optimal play) Proposer wins with high probability if $s<A\sqrt{n}\log{n}$, while Decider wins with high probability if $s>B\sqrt{n}\log{n}$. This is a factor of $\Theta(\sqrt{\log{n}})$ larger than the lower bound coming from the off-diagonal Ramsey number $r(3,s)$.


title: Absolving the SSINS of Precision Interferometric Radio Data: A New Technique for Mitigating Faint Radio Frequency Interference

abstract: We introduce a new pipeline for analyzing and mitigating radio frequency interference (RFI), which we call Sky-Subtracted Incoherent Noise Spectra (SSINS). SSINS is designed to identify and remove faint RFI below the single baseline thermal noise by employing a frequency-matched detection algorithm on baseline-averaged amplitudes of time-differenced visibilities. We demonstrate the capabilities of SSINS using the Murchison Widefield Array (MWA) in Western Australia. We successfully image aircraft flying over the array via digital television (DTV) reflection detected using SSINS and summarize an RFI occupancy survey of MWA Epoch of Reionization data. We describe how to use SSINS with new data using a documented, publicly available implementation with comprehensive usage tutorials.


title: Exploring weight initialization, diversity of solutions, and degradation in recurrent neural networks trained for temporal and decision-making tasks

abstract: Recurrent Neural Networks (RNNs) are frequently used to model aspects of brain function and structure. In this work, we trained small fully-connected RNNs to perform temporal and flow control tasks with time-varying stimuli. Our results show that different RNNs can solve the same task by converging to different underlying dynamics and also how the performance gracefully degrades as either network size is decreased, interval duration is increased, or connectivity damage is increased. For the considered tasks, we explored how robust the network obtained after training can be according to task parameterization. In the process, we developed a framework that can be useful to parameterize other tasks of interest in computational neuroscience. Our results are useful to quantify different aspects of the models, which are normally used as black boxes and need to be understood in order to model the biological response of cerebral cortex areas.


title: Robust Gaussian Process Regression for Real-Time High Precision GPS Signal Enhancement

abstract: Satellite-based positioning system such as GPS often suffers from large amount of noise that degrades the positioning accuracy dramatically especially in real-time applications. In this work, we consider a data-mining approach to enhance the GPS signal. We build a large-scale high precision GPS receiver grid system to collect real-time GPS signals for training. The Gaussian Process (GP) regression is chosen to model the vertical Total Electron Content (vTEC) distribution of the ionosphere of the Earth. Our experiments show that the noise in the real-time GPS signals often exceeds the breakdown point of the conventional robust regression methods resulting in sub-optimal system performance. We propose a three-step approach to address this challenge. In the first step we perform a set of signal validity tests to separate the signals into clean and dirty groups. In the second step, we train an initial model on the clean signals and then reweigting the dirty signals based on the residual error. A final model is retrained on both the clean signals and the reweighted dirty signals. In the theoretical analysis, we prove that the proposed three-step approach is able to tolerate much higher noise level than the vanilla robust regression methods if two reweighting rules are followed. We validate the superiority of the proposed method in our real-time high precision positioning system against several popular state-of-the-art robust regression methods. Our method achieves centimeter positioning accuracy in the benchmark region with probability $78.4\%$ , outperforming the second best baseline method by a margin of $8.3\%$. The benchmark takes 6 hours on 20,000 CPU cores or 14 years on a single CPU.


title: On the divergence of Birkhoff Normal Forms

abstract: It is well known that a real analytic symplectic diffeomorphism of the $2d$-dimensional disk ($d\geq 1$) admitting the origin as a non-resonant elliptic fixed can be {\it formally} conjugated to its Birkhoff Normal Form, a formal power series defining a {\it formal integrable} symplectic diffeomorphism at the origin. We prove in this paper that this Birkhoff Normal Form is in general divergent. This solves, in any dimension, the question of determining which of the two alternatives of Perez-Marco's theorem \cite{PM} is true and answers a question by H. Eliasson. Our result is a consequence of the fact that when $d=1$ the convergence of the formal object that is the BNF has strong dynamical consequences on the Lebesgue measure of the set of invariant circles in arbitrarily small neighborhoods of the origin. Our proof, as well as our results, extend to the case of real-analytic diffeomorphisms of the annulus admitting a Diophantine invariant torus.


title: Explicit $L^2$ bounds for the Riemann $\zeta$ function

abstract: Explicit bounds on the tails of the zeta function $\zeta$ are needed for applications, notably for integrals involving $\zeta$ on vertical lines or other paths going to infinity. Here we bound weighted $L^2$ norms of tails of $\zeta$. Two approaches are followed, each giving the better result on a different range. The first one is inspired by the proof of the standard mean value theorem for Dirichlet polynomials. The second approach, superior for large $T$, is based on classical lines, starting with an approximation to $\zeta$ via Euler-Maclaurin. Both bounds give main terms of the correct order for $0<\sigma\leq 1$ and are strong enough to be of practical use for the rigorous computation of improper integrals. We also present bounds for the $L^{2}$ norm of $\zeta$ in $[1,T]$ for $0\leq\sigma\leq 1$.


title: Three-dimensional modeling of chromospheric spectral lines in a simulated active region

abstract: Because of the complex physics that governs the formation of chromospheric lines, interpretation of solar chromospheric observations is difficult. The origin and characteristics of many chromospheric features are, because of this, unresolved. We focus here on studying two prominent features: long fibrils and flare ribbons. To model them, we use a 3D MHD simulation of an active region which self-consistently reproduces both of them. We model the H$\alpha$, Mg II k, Ca II K, and Ca II 8542 {\AA} lines using the 3D non-LTE radiative transfer code Multi3D. This simulation reproduces long fibrils that span between the opposite-polarity sunspots and go up to 4 Mm in height. They can be traced in all lines due to density corrugation. Opposite to previous studies, H$\alpha$, Mg II h&k, and Ca II H&K, are formed at similar height in this model. Magnetic field lines are aligned with the H$\alpha$ fibrils, but the latter holds to a lesser extent for the Ca II 8542 {\AA} line. The simulation shows structures in the H$\alpha$ line core that look like flare ribbons. The emission in the ribbons is caused by a dense chromosphere and a transition region at high column mass. The ribbons are visible in all chromospheric lines, but least prominent in Ca II 8542 {\AA} line. In some pixels, broad asymmetric profiles with a single emission peak are produced, similar to the profiles observed in flare ribbons. They are caused by a deep onset of the chromospheric temperature rise and large velocity gradients. The simulation produces long fibrils similar to what is seen in observations. It also produces structures similar to flare ribbons despite the lack of non-thermal electrons in the simulation. The latter suggests that thermal conduction might be a significant agent in transporting flare energy to the chromosphere in addition to non-thermal electrons.


title: Integrated Access and Backhaul in 5G mmWave Networks: Potentials and Challenges

abstract: Integrated Access and Backhaul (IAB) is being investigated as a means to overcome deployment costs of ultra-dense 5G millimeter wave (mmWave) networks by realizing wireless backhaul links to relay the access traffic. For the development of these systems, however, it is fundamental to validate the performance of IAB in realistic scenarios through end-to-end system level simulations. In this paper, we shed light on the most recent standardization activities on IAB, and compare architectures with and without IAB in mmWave deployments. While it is well understood that IAB networks reduce deployment costs by obviating the need to provide wired backhaul to each cellular base-station, in this paper we demonstrate the cell-edge throughput advantage offered by IAB using end-to-end system level simulations. We further highlight some research challenges associated with this architecture that will require further investigations.


title: A Dyadic IRT Model

abstract: We propose a dyadic Item Response Theory (dIRT) model for measuring interactions of pairs of individuals when the responses to items represent the actions (or behaviors, perceptions, etc.) of each individual (actor) made within the context of a dyad formed with another individual (partner). Examples of its use include the assessment of collaborative problem solving, or the evaluation of intra-team dynamics. The dIRT model generalizes both Item Response Theory (IRT) models for measurement and the Social Relations Model (SRM) for dyadic data. The responses of an actor when paired with a partner are modeled as a function of not only the actor's inclination to act and the partner's tendency to elicit that action, but also the unique relationship of the pair, represented by two directional, possibly correlated, interaction latent variables. Generalizations are discussed, such as accommodating triads or larger groups. Estimation is performed using Markov-chain Monte Carlo implemented in Stan, making it straightforward to extend the dIRT model in various ways. Specifically, we show how the basic dIRT model can be extended to accommodate latent regressions, multilevel settings with cluster-level random effects, as well as joint modeling of dyadic data and a distal outcome. A simulation study demonstrates that estimation performs well. We apply our proposed approach to speed-dating data and find new evidence of pairwise interactions between participants, describing a mutual attraction that is inadequately characterized by individual properties alone.


title: MEMe: An Accurate Maximum Entropy Method for Efficient Approximations in Large-Scale Machine Learning

abstract: Efficient approximation lies at the heart of large-scale machine learning problems. In this paper, we propose a novel, robust maximum entropy algorithm, which is capable of dealing with hundreds of moments and allows for computationally efficient approximations. We showcase the usefulness of the proposed method, its equivalence to constrained Bayesian variational inference and demonstrate its superiority over existing approaches in two applications, namely, fast log determinant estimation and information-theoretic Bayesian optimisation.


title: Do place cells dream of conditional probabilities? Learning Neural Nystr\"om representations

abstract: We posit that hippocampal place cells encode information about future locations under a transition distribution observed as an agent explores a given (physical or conceptual) space. The encoding of information about the current location, usually associated with place cells, then emerges as a necessary step to achieve this broader goal. We formally derive a biologically-inspired neural network from Nystr\"om kernel approximations and empirically demonstrate that the network successfully approximates transition distributions. The proposed network yields representations that, just like place cells, soft-tile the input space with highly sparse and localized receptive fields. Additionally, we show that the proposed computational motif can be extended to handle supervised problems, creating class-specific place cells while exhibiting low sample complexity.


title: Attractive versus truncated repulsive supercooled liquids: The dynamics is encoded in the pair correlation function

abstract: We compare glassy dynamics in two liquids that differ in the form of their interaction potentials. Both systems have the same repulsive interactions but one has also an attractive part in the potential. These two systems exhibit very different dynamics despite having nearly identical pair correlation functions. We demonstrate that a properly weighted integral of the pair correlation function, which amplifies the subtle differences between the two systems, correctly captures their dynamical differences. The weights are obtained from a standard machine learning algorithm.

categories: cond-mat.stat-mech cond-mat.dis-nn cond-mat.soft cs.LG


title: Rauzy induction of polygon partitions and toral $\mathbb{Z}^2$-rotations

abstract: We extend the notion of Rauzy induction of interval exchange transformations to the case of toral $\mathbb{Z}^2$-rotation, i.e., $\mathbb{Z}^2$-action defined by rotations on a 2-torus. If $\mathcal{X}_{\mathcal{P},R}$ denotes the symbolic dynamical system corresponding to a partition $\mathcal{P}$ and $\mathbb{Z}^2$-action $R$ such that $R$ is Cartesian on a sub-domain $W$, we express the 2-dimensional configurations in $\mathcal{X}_{\mathcal{P},R}$ as the image under a $2$-dimensional morphism (up to a shift) of a configuration in $\mathcal{X}_{\widehat{\mathcal{P}}|_W,\widehat{R}|_W}$ where $\widehat{\mathcal{P}}|_W$ is the induced partition and $\widehat{R}|_W$ is the induced $\mathbb{Z}^2$-action on $W$. We focus on one example $\mathcal{X}_{\mathcal{P}_0,R_0}$ for which we obtain an eventually periodic sequence of 2-dimensional morphisms. We prove that it is the same as the substitutive structure of the minimal subshift $X_0$ of the Jeandel-Rao Wang shift computed in an earlier work by the author. As a consequence, $\mathcal{P}_0$ is a Markov partition for the associated toral $\mathbb{Z}^2$-rotation $R_0$. It also implies that the subshift $X_0$ is uniquely ergodic and is isomorphic to the toral $\mathbb{Z}^2$-rotation $R_0$ which can be seen as a generalization for 2-dimensional subshifts of the relation between Sturmian sequences and irrational rotations on a circle. Batteries included: the algorithms and code to reproduce the proofs are provided.


title: Training Neural Machine Translation To Apply Terminology Constraints

abstract: This paper proposes a novel method to inject custom terminology into neural machine translation at run time. Previous works have mainly proposed modifications to the decoding algorithm in order to constrain the output to include run-time-provided target terms. While being effective, these constrained decoding methods add, however, significant computational overhead to the inference step, and, as we show in this paper, can be brittle when tested in realistic conditions. In this paper we approach the problem by training a neural MT system to learn how to use custom terminology when provided with the input. Comparative experiments show that our method is not only more effective than a state-of-the-art implementation of constrained decoding, but is also as fast as constraint-free decoding.


title: Polyelectrolyte solution under spatial and dielectric confinement

abstract: Polyelectrolytes under confinement are crucial for energy storage and for understanding biomolecular functions. Using molecular dynamics simulations, we analyze a polyelectrolyte solution confined between two oppositely charged planar dielectric surfaces and include surface polarization effects due to dielectric mismatch at the two electrodes. Although the effect of polarization on the charge distribution seems minor, we find that surface polarization enhances energy storage and also leads to the emergence of negative differential capacitance in confined polyelectrolyte solutions.


title: A stochastic variant of the abelian sandpile model

abstract: We introduce a natural stochastic extension, called SSP, of the abelian sandpile model(ASM), which shares many mathematical properties with ASM, yet radically differs in its physical behavior, for example in terms of the shape of the steady state and of the avalanche size distribution. We establish a basic theory of SSP analogous to that of ASM, and present a brief numerical study of its behavior. Our original motivation for studying SSP stems from its connection to the LLL algorithm established in another work by the authors [5]. The importance of understanding how LLL works cannot be stressed more, especially from the point of view of lattice-based cryptography. We believe SSP serves as a tractable toy model of LLL that would help further our understanding of it.


title: Socially Inspired Communication in Swarm Robotics

abstract: Localized communication in swarms has been shown to increase swarm effectiveness in some situations by allowing for additional opportunities for cooperation. However, communication and utilization of potentially outdated information is also a concern. We present an explicit non-directional goal-based communication model and message accept/reject scheme, and test our model in a set of object gathering experiments with a swarm of robots. The results of the experiments indicate that even low levels of communication regarding the swarm's goal outperform high levels of random information communication.


title: Elliott-Yafet Spin-Phonon Relaxation Times from First Principles

abstract: We present a first-principles approach for computing the phonon-limited $T_1$ spin relaxation time due to the Elliot-Yafet mechanism. Our scheme combines fully-relativistic spin-flip electron-phonon interactions with an approach to compute the effective spin of band electrons in materials with inversion symmetry. We apply our method to silicon and diamond, for which we compute the temperature dependence of the spin relaxation times and analyze the contributions to spin relaxation from different phonons and valley processes. The computed spin relaxation times in silicon are in excellent agreement with experiment in the 50$-$300 K temperature range. In diamond, we predict intrinsic spin relaxation times of 540 $\mu$s at 77 K and 2.3 $\mu$s at 300 K. Our work enables precise predictions of spin-phonon relaxation times in a wide range of materials, providing microscopic insight into spin relaxation and guiding the development of spin-based quantum technologies.


title: RL-Based Method for Benchmarking the Adversarial Resilience and Robustness of Deep Reinforcement Learning Policies

abstract: This paper investigates the resilience and robustness of Deep Reinforcement Learning (DRL) policies to adversarial perturbations in the state space. We first present an approach for the disentanglement of vulnerabilities caused by representation learning of DRL agents from those that stem from the sensitivity of the DRL policies to distributional shifts in state transitions. Building on this approach, we propose two RL-based techniques for quantitative benchmarking of adversarial resilience and robustness in DRL policies against perturbations of state transitions. We demonstrate the feasibility of our proposals through experimental evaluation of resilience and robustness in DQN, A2C, and PPO2 policies trained in the Cartpole environment.


title: Rotation and pulsation in Ap stars: first light results from TESS sectors 1 and 2

abstract: We present the first results from the Transiting Exoplanet Survey Satellite (TESS) on the rotational and pulsational variability of magnetic chemically peculiar A-type stars. We analyse TESS 2-min cadence data from sectors 1 and 2 on a sample of 83 stars. Five new rapidly oscillating Ap (roAp) stars are announced. One of these pulsates with periods around 4.7 min, making it the shortest period roAp star known to date. Four out of the five new roAp stars are multiperiodic. Three of these, and the singly-periodic one show the presence of rotational mode splitting. Individual frequencies are provided in all cases. In addition, seven previously known roAp stars are analysed. Additional modes of oscillation are found in some stars, while in others we are able to distinguish the true pulsations from possible aliases present in the ground-based data. We find that the pulsation amplitude in the TESS filter is typically a factor 6 smaller than that in the $B$ filter which is usually used for ground-based observations. For four roAp stars we set constraints on the inclination angle and magnetic obliquity, through the application of the oblique pulsator model. We also confirm the absence of roAp-type pulsations down to amplitude limits of 6 and 13 micromag, respectively, in two of the best characterised non-oscillating Ap (noAp) stars. We announce 27 new rotational variables along with their rotation periods, and provide different rotation periods for seven other stars. Finally, we discuss how these results challenge state-of-the-art pulsation models for roAp stars.