title: Nearest matrix polynomials with a specified elementary divisor

abstract: The problem of finding the distance from a given $n \times n$ matrix polynomial of degree $k$ to the set of matrix polynomials having the elementary divisor $(\lambda-\lambda_0)^j, \, j \geqslant r,$ for a fixed scalar $\lambda_0$ and $2 \leqslant r \leqslant kn$ is considered. It is established that polynomials that are not regular are arbitrarily close to a regular matrix polynomial with the desired elementary divisor. For regular matrix polynomials the problem is shown to be equivalent to finding minimal structure preserving perturbations such that a certain block Toeplitz matrix becomes suitably rank deficient. This is then used to characterize the distance via two different optimizations. The first one shows that if $\lambda_0$ is not already an eigenvalue of the matrix polynomial, then the problem is equivalent to computing a generalized notion of a structured singular value. The distance is computed via algorithms like BFGS and Matlab's globalsearch algorithm from the second optimization. Upper and lower bounds of the distance are also derived and numerical experiments are performed to compare them with the computed values of the distance.


title: Locally interacting diffusions as space-time Markov random fields

abstract: We consider a countable system of interacting (possibly non-Markovian) stochastic differential equations driven by independent Brownian motions and indexed by the vertices of a locally finite graph $G = (V,E)$. The drift of the process at each vertex is influenced by the states of that vertex and its neighbors, and the diffusion coefficient depends on the state of only that vertex. Such processes arise in a variety of applications including statistical physics, neuroscience, engineering and math finance. Under general conditions on the coefficients, we show that if the initial conditions form a second-order Markov random field on $d$-dimensional Euclidean space, then at any positive time, the collection of histories of the processes at different vertices forms a second-order Markov random field on path space. We also establish a bijection between (second-order) Gibbs measures on $(\R^d)^V$ (with finite second moments) and a set of space-time (second-order) Gibbs measures on path space, corresponding respectively to the initial law and the law of the solution to the stochastic differential equation. As a corollary, we establish a Gibbs uniqueness property that shows that for infinite graphs the joint distribution of the paths is completely determined by the initial condition and the specifications, namely the family of conditional distributions on finite vertex sets given the configuration on the complement. Along the way, we establish various approximation and projection results for Markov random fields on locally finite graphs that may be of independent interest.


title: Poisson and Gaussian fluctuations for the $\mathbf{f}$-vector of high-dimensional random simplicial complexes

abstract: We investigate the high-dimensional asymptotic distributional behavior of the $\mathbf{f}$-vector of a random Vietoris-Rips complex, that is generated over a stationary Poisson point process in $[-\frac{1}{2}, \frac{1}{2}]^d$ as the space dimension and the intensity tend to infinity while the radius parameter tends to zero simultaneously.


title: On some quasi-analytic classes

abstract: Using the so called monotonicity property, we prove that the Borel mapping restricted to some quasi-anlytic classes is never onto.


title: Parameter uncertainties in weighted unbinned maximum likelihood fits

abstract: Parameter estimation via unbinned maximum likelihood fits is central for many analyses performed in high energy physics. Unbinned maximum likelihood fits using event weights, for example to statistically subtract background contributions via the sPlot formalism, or to correct for acceptance effects, have recently seen increasing use in the community. However, it is well known that the naive approach to the estimation of parameter uncertainties via the second derivative of the logarithmic likelihood does not yield confidence intervals with the correct coverage in the presence of event weights. This paper derives the asymptotically correct expressions and compares them with several commonly used approaches for the determination of parameter uncertainties, some of which are shown to not generally be asymptotically correct. In addition, the effect of uncertainties on event weights is discussed, including uncertainties that can arise from the presence of nuisance parameters in the determination of sWeights.


title: Fetal cardiovascular decompensation during labor predicted from the individual heart rate: a prospective study in fetal sheep near term and the impact of low sampling rate

abstract: We present a novel computerized fetal heart rate intrapartum algorithm for early and individualized prediction of fetal cardiovascular decompensation, a key event in the causal chain leading to brain injury. This real-time machine learning algorithm performs well on noisy fetal heart rate data and requires ~2 hours to train on the individual fetal heart rate tracings in the first stage of labor; once trained, the algorithm predicts the event of fetal cardiovascular decompensation with 92% sensitivity. We show that the algorithm's performance suffers reducing sensitivity to 67% when the fetal heart rate is acquired at the sampling rate of 4 Hz used in ultrasound cardiotocographic monitors compared to the electrocardiogram(ECG)-derived signals as can be acquired from maternal abdominal ECG.


title: On braided co-Frobenius Hopf algebras

abstract: We present characterizations of braided co-Frobenius Hopf algebras in the braided tensor category of Yetter-Drinfeld modules over a Hopf algebra extending those already known for co-Frobenius Hopf algebras.


abstract: This thesis is entitled Dynamic programming systems for modeling and control of the traffic in transportation networks. Two parts are distinguished in this dissertation: 1) methods and approaches based on min-plus or max-plus algebra, where the dynamics are deterministic dynamic programming systems; 2) methods and approaches whose dynamic systems are non-linear but are interpreted as stochastic dynamic programming systems. Each of the two parts includes a chapter of necessary reviews, two main chapters and a chapter summarizing other works related to the concerned part. Part 1 includes a first chapter containing an introduction and some necessary reviews; two main chapters, one on the max-plus algebra model for the train dynamics on a metro line, the other one on the network calculus approach for modeling and calculating performance bounds on road networks; and a final chapter summarizing my other contributions on the topic of this part. Part 2 includes a first chapter containing an introduction and some necessary reviews; two main chapters, one on the microscopic modeling of traffic taking into account anticipation in driving, the other one on the modeling of the train dynamics on a metro line taking into account the passenger travel demand; and a final chapter summarizing my other contributions on the topic of this part.


title: Lie Theory suffices to understand, and Locally Resolve, the Problem of Time

abstract: The Lie claw digraph controls Background Independence and thus the Problem of Time and indeed the Fundamental Nature of Physical Law. This has been established in the realms of Flat and Differential Geometry with varying amounts of extra mathematical structure. This Lie claw digraph has Generator Closure at its centre (Lie brackets), Relationalism at its root (implemented by Lie derivatives), and, as its leaves, Assignment of Observables (zero commutants under Lie brackets) and Constructability from Less Structure Assumed (working if generator Deformation leads to Lie brackets algebraic Rigidity). This centre is enabled by automorphisms and powered by the Generalized Lie Algorithm extension of the Dirac Algorithm (itself sufficing for the canonical subcase, for which generators are constraints). The Problem of Time's facet ordering problem is resolved.


title: Nonequilibrium Phase Transition in Constrained Adsorption

abstract: We study the adsorption-desorption of fluid molecules on a solid substrate by introducing a schematic model in which the adsorption/desorption transition probabilities are given by irreversible kinetic constraints with a tunable violation of local detailed balance condition. Numerical simulations show that in one spatial dimension the model undergoes a continuous nonequilibrium phase transition whose location depends on the irreversibility strength. We show that the hierarchy of equations obeyed by multi-point correlation functions can be closed to the second order by means of a simple decoupling approximation, and that the approximated solution for the steady state yields a very good description of the overall phase diagram.


title: Jet quenching test of the QCD matter created at RHIC and the LHC needs opacity-resummed medium induced radiation

abstract: After almost two decades of investigation, jet quenching has become a fundamental tool to study the properties of the QCD matter produced in high-energy nuclear collisions. Despite the large progress in both experimental and theoretical tools, several unknowns remain to be solved. In particular the systematics of the QCD matter parameters extracted at different energies and centralities is still puzzling. On the other hand, a lot of attention has been put in the last few years to a correct resummation of the multiple scatterings that the particles of the jet suffer when traversing the QCD matter. Present formalisms have substantially improved the multiple-soft or single hard scattering approximations, common in the last twenty years. We show here that medium parameters extracted with full resummed medium-induced gluon radiation spectra significantly differ from these two extreme approximations, providing a natural systematics for different collision energies and centralities, with no inconsistencies. For the case of the full resummation this ensures that the density of scattering centers as a function of the temperature qualitatively agree with the equation of state computed in lattice QCD.


title: Tustin neural networks: a class of recurrent nets for adaptive MPC of mechanical systems

abstract: The use of recurrent neural networks to represent the dynamics of unstable systems is difficult due to the need to properly initialize their internal states, which in most of the cases do not have any physical meaning, consequent to the non-smoothness of the optimization problem. For this reason, in this paper focus is placed on mechanical systems characterized by a number of degrees of freedom, each one represented by two states, namely position and velocity. For these systems, a new recurrent neural network is proposed: Tustin-Net. Inspired by second-order dynamics, the network hidden states can be straightforwardly estimated, as their differential relationships with the measured states are hardcoded in the forward pass. The proposed structure is used to model a double inverted pendulum and for model-based Reinforcement Learning, where an adaptive Model Predictive Controller scheme using the Unscented Kalman Filter is proposed to deal with parameter changes in the system.


title: The shape of the cosmic ray proton spectrum

abstract: Recent observations of cosmic ray protons in the energy range $10^2$--$10^5$~GeV have revealed that the spectrum cannot be described by a simple power law. A hardening of the spectrum around an energy of order few hundred~GeV, first observed by the magnetic spectrometers PAMELA and AMS02, has now been confirmed by several calorimeter detectors (ATIC, CREAM, CALET, NUCLEON and DAMPE). These new measurements reach higher energy and indicate that the hardening corresponds to a larger step in spectral index than what estimated by the magnetic spectrometers. Data at still higher energy (by CREAM, NUCLEON and DAMPE) show that the proton spectrum undergoes a marked softening at $E \approx 10^4$~GeV. Understanding the origin of these unexpected spectral features is a significant challenge for models of the Galactic cosmic rays. An important open question is whether additional features are present in the proton spectrum between the softening and the "Knee". Extensive Air Shower detectors, using unfolding procedures that require the modeling of cosmic ray showers in the atmosphere, estimated the proton flux below and around the Knee (at $E \simeq 3$~PeV). These results however have large systematic uncertainties and are in poor agreement with each other. The measurement in the PeV energy range, recently presented by IceTop/IceCube, indicates a proton flux higher than extrapolations of the direct measurements calculated assuming a constant slope, and therefore requires the existence of an additional spectral hardening below the Knee. A clarification of this point is very important for an understanding of the origin of the Galactic cosmic rays, and is also essential for a precise calculation of the spectra of atmospheric neutrinos in the energy range ($E \gtrsim 10$~TeV) where they constitute the foreground for the emerging astrophysical $\nu$ signal.


title: Advancing Fourier: space-time concepts in ultrafast optics, imaging and photonic neural networks

abstract: The concepts of Fourier optics were established in France in the 1940s by Pierre-Michel Duffieux, and laid the foundations of an extensive series of activities in the French research community that have touched on nearly every aspect of contemporary optics and photonics. In this paper, we review a selection of results where applications of the Fourier transform and transfer functions in optics have been applied to yield significant advances in unexpected areas of optics, including the spatial shaping of complex laser beams in amplitude and in phase, real-time ultrafast measurements, novel ghost imaging techniques, and the development of parallel processing methodologies for photonic artificial intelligence.


title: Rigidly rotating gravitationally bound systems of point particles, compared to polytropes

abstract: In order to simulate rigidly rotating polytropes we have simulated systems of $N$ point particles, with $N$ up to 1800. Two particles at a distance $r$ interact by an attractive potential $-1/r$ and a repulsive potential $1/r^2$. The repulsion simulates the pressure in a polytropic gas of polytropic index $3/2$. We take the total angular momentum $L$ to be conserved, but not the total energy $E$. The particles are stationary in the rotating coordinate system. The rotational energy is $L^2/(2I)$ where $I$ is the moment of inertia. Configurations where the energy $E$ has a local minimum are stable. In the continuum limit $N\to\infty$ the particles become more and more tightly packed in a finite volume, with the interparticle distances decreasing as $N^{-1/3}$. We argue that $N^{-1/3}$ is a good parameter for describing the continuum limit. We argue further that the continuum limit is the polytropic gas of index $3/2$. For example, the density profile of the nonrotating gas approaches that computed from the Lane--Emden equation describing the nonrotating polytropic gas. In the case of maximum rotation the instability occurs by the loss of particles from the equator, which becomes a sharp edge, as predicted by Jeans in his study of rotating polytropes. We describe the minimum energy nonrotating configurations for a number of small values of $N$.


title: Understanding the plasma and magnetic field evolution of a filament using observations and Nonlinear force-free field modelling

abstract: We present observations and magnetic field models of an intermediate filament present on the Sun in August 2012, associated with a polarity inversion line that extends from AR 11541 in the east into the quiet sun at its western end. A combination of SDO/AIA, SDO/HMI, and GONG H alpha data allow us to analyse the structure and evolution of the filament from 2012 August 4 23:00 UT to 2012 August 6 08:00 UT when the filament was in equilibrium. By applying the flux rope insertion method, nonlinear force-free field models of the filament are constructed using SDO/HMI line-of-sight magnetograms as the boundary condition at the two times given above. Guided by observed filament barbs, both modelled flux ropes are split into three sections each with a different value of axial flux to represent the non-uniform photospheric field distribution. The flux in the eastern section of the rope increases by 4$\times$10$^{20}$ Mx between the two models, which is in good agreement with the amount of flux cancelled along the internal PIL of AR 11541, calculated to be 3.2$\times$10$^{20}$ Mx. This suggests that flux cancellation builds flux into the filament's magnetic structure. Additionally, the number of field line dips increases between the two models in the locations where flux cancellation, the formation of new filament threads and growth of the filament is observed. This suggests that flux cancellation associated with magnetic reconnection forms concave-up magnetic field that lifts plasma into the filament. During this time, the free magnetic energy in the models increases by 0.2$\times$10$^{31}$ ergs.


title: Forecasting Megaelectron-Volt Electrons inside Earth's Outer Radiation Belt: PreMevE 2.0 Based on Supervised Machine Learning Algorithms

abstract: Here we present the recent progress in upgrading a predictive model for Megaelectron-Volt (MeV) electrons inside the Earth's outer Van Allen belt. This updated model, called PreMevE 2.0, is demonstrated to make much improved forecasts, particularly at outer Lshells, by including upstream solar wind speeds to the model's input parameter list. Furthermore, based on several kinds of linear and artificial machine learning algorithms, a list of models were constructed, trained, validated and tested with 42-month MeV electron observations from Van Allen Probes. Out-of-sample test results from these models show that, with optimized model hyperparameters and input parameter combinations, the top performer from each category of models has the similar capability of making reliable 1-day (2-day) forecasts with Lshell-averaged performance efficiency values ~ 0.87 (~0.82). Interestingly, the linear regression model is often the most successful one when compared to other models, which indicates the relationship between 1 MeV electron dynamics and precipitating electrons is dominated by linear components. It is also shown that PreMevE 2.0 can reasonably predict the onsets of MeV electron events in 2-day forecasts. This improved PreMevE model is driven by observations from longstanding space infrastructure (a NOAA satellite in low-Earth-orbit, the solar wind monitor at the L1 point, and one LANL satellite in geosynchronous orbit) to make high-fidelity forecasts for MeV electrons, and thus can be an invaluable space weather forecasting tool for the future.


title: Robust low-delay Streaming PIR using convolutional codes

abstract: In this paper we investigate the design of a low-delay robust streaming PIR scheme on coded data that is resilient to unresponsive or slow servers and can privately retrieve streaming data in a sequential fashion subject to a fixed decoding delay. We present a scheme based on convolutional codes and the star product and assume no collusion between servers. In particular we propose the use of convolutional codes that have the maximum distance increase, called Maximum Distance Profile (MDP). We show that the proposed scheme can deal with many different erasure patterns.


title: Optimization of photoluminescence from W centers in silicon-on-insulator

abstract: W centers are trigonal defects generated by self-ion implantation in silicon that exhibit photoluminescence at 1.218 $\mu$m. We have shown previously that they can be used in waveguide-integrated all-silicon light-emitting diodes (LEDs). Here we optimize the implant energy, fluence and anneal conditions to maximize the photoluminescence intensity for W centers implanted in silicon-on-insulator, a substrate suitable for waveguide-integrated devices. After optimization, we observe near two orders of magnitude improvement in photoluminescence intensity relative to the conditions with the stopping range of the implanted ions at the center of the silicon device layer. The previously demonstrated waveguide-integrated LED used implant conditions with the stopping range at the center of this layer. We further show that such light sources can be manufactured at the 300-mm scale by demonstrating photoluminescence of similar intensity from 300 mm silicon-on-insulator wafers. The luminescence uniformity across the entire wafer is within the measurement error.


title: Sequential/Spatial, a Survey of Interactive Information Retrieval Methods for Controlled Experimentation and Evaluation

abstract: This survey presents studies that investigated non-spatial (sequential) and spatial information retrieval systems in parallel during a battery of information-seeking tasks with respect to user navigational behaviors, incidental learning, retrieval performance, cognitive abilities & load, direct manipulation of 2D & 3D interfaces, and satisfaction. I consider how information theory has contributed to the concepts of foraging, sense-making, exploration, and how the applied areas of interactive information retrieval (IIR) and cognitive/behavioral psychology have implemented these concepts into architecture, interface design, experimental design, user study, and evaluation methodology.


title: Fast sampling and counting k-SAT solutions in the local lemma regime

abstract: We give new algorithms based on Markov chains to sample and approximately count satisfying assignments to $k$-uniform CNF formulas where each variable appears at most $d$ times. For any $k$ and $d$ satisfying $kd<n^{o(1)}$ and $k\ge 20\log k + 20\log d + 60$, the new sampling algorithm runs in close to linear time, and the counting algorithm runs in close to quadratic time. Our approach is inspired by Moitra (JACM, 2019) which remarkably utilizes the Lov\'{a}sz local lemma in approximate counting. Our main technical contribution is to use the local lemma to bypass the connectivity barrier in traditional Markov chain approaches, which makes the well developed MCMC method applicable on disconnected state spaces such as SAT solutions. The benefit of our approach is to avoid the enumeration of local structures and obtain fixed polynomial running times, even if $k=\omega(1)$ or $d=\omega(1)$.


title: Synthetic Video Generation for Robust Hand Gesture Recognition in Augmented Reality Applications

abstract: Hand gestures are a natural means of interaction in Augmented Reality and Virtual Reality (AR/VR) applications. Recently, there has been an increased focus on removing the dependence of accurate hand gesture recognition on complex sensor setup found in expensive proprietary devices such as the Microsoft HoloLens, Daqri and Meta Glasses. Most such solutions either rely on multi-modal sensor data or deep neural networks that can benefit greatly from abundance of labelled data. Datasets are an integral part of any deep learning based research. They have been the principal reason for the substantial progress in this field, both, in terms of providing enough data for the training of these models, and, for benchmarking competing algorithms. However, it is becoming increasingly difficult to generate enough labelled data for complex tasks such as hand gesture recognition. The goal of this work is to introduce a framework capable of generating photo-realistic videos that have labelled hand bounding box and fingertip that can help in designing, training, and benchmarking models for hand-gesture recognition in AR/VR applications. We demonstrate the efficacy of our framework in generating videos with diverse backgrounds.


title: Decentralized physiology and the molecular basis of social life in eusocial insects

abstract: The traditional focus of physiological and functional genomic research is on molecular processes that play out within a single body. In contrast, when social interactions occur, molecular and behavioral responses in interacting individuals can lead to physiological processes that are distributed across multiple individuals. In eusocial insect colonies, such multi-body processes are tightly integrated, involving social communication mechanisms that regulate the physiology of colony members. As a result, conserved physiological mechanisms, for example related to pheromone detection and neural signaling pathways, are deployed in novel contexts and regulate emergent colony traits during the evolutionary origin and elaboration of social complexity. Here we review conceptual frameworks for organismal and colony physiology, and highlight functional genomic, physiological, and behavioral research exploring how colony-level traits arise from physical and chemical interactions among nestmates. We highlight mechanistic work exploring how colony traits arise from physical and chemical interactions among physiologically-specialized nestmates of various developmental stages. We consider similarities and differences between organismal and colony physiology, and make specific predictions based on a decentralized perspective on the function and evolution of colony traits. Integrated models of colony physiological function will be useful to address fundamental questions related to the evolution and ecology of collective behavior in natural systems.


title: The matching condition for larger size Riemann-Hilbert problems

abstract: In a larger size Riemann-Hilbert problem matching the local parametrices with the global parametrix is often a major issue. In this article we present a result that should tackle this problem in natural situations. We prove that, in a general setting, it is possible to obtain a double matching, that is, a matching condition on two circles instead of one circle. We discuss how this matching approach can be used to obtain local scaling limits of correlation kernels and apply our result to several examples from the existing literature.


title: Characterizing multiphoton excitation using time-resolved X-ray scattering

abstract: Molecular iodine was photoexcited by a strong 800 nm laser, driving several channels of multiphoton excitation. The motion following photoexcitation was probed using time-resolved X-ray scattering, which produces a scattering map $S(Q,\tau)$. Temporal Fourier transform methods were employed to obtain a frequency-resolved X-ray scattering signal $\tilde{S}(Q,\omega)$. Taken together, $S(Q,\tau)$ and $\tilde{S}(Q,\omega)$ separate different modes of motion, so that mode-specific nuclear oscillatory positions, oscillation amplitudes, directions of motions, and times may be measured accurately. Molecular dissociations likewise have a distinct signature, which may be used to identify both velocities and dissociation time shifts, and also can reveal laser-induced couplings among the molecular potentials.


title: Examining UK drill music through sentiment trajectory analysis

abstract: This paper presents how techniques from natural language processing can be used to examine the sentiment trajectories of gang-related drill music in the United Kingdom (UK). This work is important because key public figures are loosely making controversial linkages between drill music and recent escalations in youth violence in London. Thus, this paper examines the dynamic use of sentiment in gang-related drill music lyrics. The findings suggest two distinct sentiment use patterns and statistical analyses revealed that lyrics with a markedly positive tone attract more views and engagement on YouTube than negative ones. Our work provides the first empirical insights into the language use of London drill music, and it can, therefore, be used in future studies and by policymakers to help understand the alleged drill-gang nexus.


title: Optimal Transport Based Change Point Detection and Time Series Segment Clustering

abstract: Two common problems in time series analysis are the decomposition of the data stream into disjoint segments that are each in some sense "homogeneous" - a problem known as Change Point Detection (CPD) - and the grouping of similar nonadjacent segments, a problem that we call Time Series Segment Clustering (TSSC). Building upon recent theoretical advances characterizing the limiting distribution-free behavior of the Wasserstein two-sample test (Ramdas et al. 2015), we propose a novel algorithm for unsupervised, distribution-free CPD which is amenable to both offline and online settings. We also introduce a method to mitigate false positives in CPD and address TSSC by using the Wasserstein distance between the detected segments to build an affinity matrix to which we apply spectral clustering. Results on both synthetic and real data sets show the benefits of the approach.


title: Quantum thermal transistor in superconducting circuits

abstract: Logical devices based on electrical currents are ubiquitous in modern society. However, digital logic does have some drawbacks such as a relatively high power consumption. It is therefore of great interest to seek alternative means to build logical circuits that can either work as stand-alone devices or in conjunction with more traditional electronic circuits. One direction that holds great promise is the use of heat currents for logical components. In the present paper, we discuss a recent abstract proposal for a quantum thermal transistor and provide a concrete design of such a device using superconducting circuits. Using a circuit quantum electrodynamics Jaynes-Cummings model, we propose a three-terminal device that allows heat transfer from source to drain, depending on the temperature of a bath coupled at the gate modulator, and show that it provides similar properties to a conventional semiconductor transistor.


title: Noncoercive Lyapunov functions for input-to-state stability of infinite-dimensional systems

abstract: We consider an abstract class of infinite-dimensional dynamical systems with inputs. For this class, the significance of noncoercive Lyapunov functions is analyzed. It is shown that the existence of such Lyapunov functions implies norm-to-integral input-to-state stability. This property in turn is equivalent to input-to-state stability, if the system satisfies certain mild regularity assumptions. For a particular class of linear systems with unbounded admissible input operators, explicit constructions of noncoercive Lyapunov functions are provided. The theory is applied to a heat equation with Dirichlet boundary conditions.


title: Theory Summary at Strangeness in Quark Matter 2019

abstract: This is the theory summary of Strangeness in Quark Matter 2019 conference. Results include the state-of-the-art updates to the Quantum Chromodynamics (QCD) phase diagram with contributions both from heavy-ion collisions and nuclear astrophysics, studies on the QCD freeze-out lines, and several aspects regarding small systems including collectivity, heavy flavor dynamics, strangeness, and hard probes.


title: Constraints on Running of Non-Gaussianity from Large Scale Structure Probes

abstract: In this letter we present constraints on the scale-dependent "local" type primordial non-Gaussianity, which is described by non-Gaussianity's spectral index $n_{\rm NG}$, from the NRAO VLA Sky Survey and the quasar catalog of the Sloan Digital Sky Survey (SDSS) Data Release 6, together with the SDSS Data Release 12 photo-z sample. Here, we use the auto-correlation analyses of these three probes and their cross-correlation analyses with the cosmic microwave background (CMB) temperature map, and obtain the tight constraint on the spectral index: $n_{\rm NG}=0.2 ^{+0.7}_{-1.0}$ ($1\,\sigma$ C.L.), which shows the first competitive constraint on the running of non-Gaussianity from current large-scale structure clustering data. Furthermore, we also perform the forecast calculations and improve the limit of $n_{\rm NG}$ using the future Euclid mission, and obtain the standard deviation at 68\% confidence level: $\Delta n_{\rm NG}=1.74$ when considering the fiducial value $f_{\rm NL}=3$, which provides the complementary constraining power to those from the CMB bispectrum information.


title: Bitcoin Coin Selection with Leverage

abstract: We present a new Bitcoin coin selection algorithm, "coin selection with leverage", which aims to improve upon cost savings than that of standard knapsack like approaches. Parameters to the new algorithm are available to be tuned at the users discretion to address other goals of coin selection. Our approach naturally fits as a replacement for the standard knapsack ingredient of full coin selection procedures.


title: Annotated Hypergraphs: Models and Applications

abstract: Hypergraphs offer a natural modeling language for studying polyadic interactions between sets of entities. Many polyadic interactions are asymmetric, with nodes playing distinctive roles. In an academic collaboration network, for example, the order of authors on a paper often reflects the nature of their contributions to the completed work. To model these networks, we introduce \emph{annotated hypergraphs} as natural polyadic generalizations of directed graphs. Annotated hypergraphs form a highly general framework for incorporating metadata into polyadic graph models. To facilitate data analysis with annotated hypergraphs, we construct a role-aware configuration null model for these structures and prove an efficient Markov Chain Monte Carlo scheme for sampling from it. We proceed to formulate several metrics and algorithms for the analysis of annotated hypergraphs. Several of these, such as assortativity and modularity, naturally generalize dyadic counterparts. Other metrics, such as local role densities, are unique to the setting of annotated hypergraphs. We illustrate our techniques on six digital social networks, and present a detailed case-study of the Enron email data set.


title: On the structure of dg categories of relative singularities

abstract: In this paper we show that every object in the dg category of relative singularities Sing$(B,\underline{f})$ associated to a pair $(B,\underline{f})$, where $B$ is a ring and $\underline{f}\in B^n$, is equivalent to a retract of a $K(B,\underline{f})$-dg module concentrated in $n+1$ degrees. When $n=1$, we show that Orlov's comparison theorem, which relates the dg category of relative singularities and that of matrix factorizations of an LG-model, holds true without any regularity assumption on the potential.


title: Using image-extracted features to determine heart rate and blink duration for driver sleepiness detection

abstract: Heart rate and blink duration are two vital physiological signals which give information about cardiac activity and consciousness. Monitoring these two signals is crucial for various applications such as driver drowsiness detection. As there are several problems posed by the conventional systems to be used for continuous, long-term monitoring, a remote blink and ECG monitoring system can be used as an alternative. For estimating the blink duration, two strategies are used. In the first approach, pictures of open and closed eyes are fed into an Artificial Neural Network (ANN) to decide whether the eyes are open or close. In the second approach, they are classified and labeled using Linear Discriminant Analysis (LDA). The labeled images are then be used to determine the blink duration. For heart rate variability, two strategies are used to evaluate the passing blood volume: Independent Component Analysis (ICA); and a chrominance based method. Eye recognition yielded 78-92% accuracy in classifying open/closed eyes with ANN and 71-91% accuracy with LDA. Heart rate evaluations had a mean loss of around 16 Beats Per Minute (BPM) for the ICA strategy and 13 BPM for the chrominance based technique.


title: Toward the Pole

abstract: We study the analytic structure of semiclassical conformal blocks, namely of the 1-point conformal block on the torus and of the 4-point conformal block on the sphere, as functions of the intermediate dimension. We interpret their discontinuities, which can be revealed with the use of a particular resummation procedure, holographically in terms of configurations of geodesics in AdS3. In other words, we study the behavior of the $\Omega$-deformed $\mathcal N=2$ SYM theory in the Nekrasov-Shatashvili limit near those singular points, where naively the W-boson with non-vanishing angular momentum becomes massless. Upon a proper resummation of instanton contributions, these singularities disappear, which is similar to the Seiberg-Witten solution in the undeformed case where there is no massless W-boson. It is shown that the order parameter undergoes a non-analytic behavior near positions of the poles. The jump of the order parameter is interpreted holographically in terms of geodesic networks


title: Remarks on the complete integrability of quantum and classical dynamical systems

abstract: It is noted that the Schrodinger equation with any self-adjoint Hamiltonian is unitary equivalent to a set of non-interacting classical harmonic oscillators and in this sense any quantum dynamics is completely integrable. Higher order integrals of motion are presented. That does not mean that we can explicitly compute the time dependence for expectation value of any quantum observable. A similar result is indicated for classical dynamical systems in terms of Koopman's approach. Explicit transformations of quantum and classical dynamics to the free evolution by using direct methods of scattering theory and wave operators are considered. Examples from classical and quantum mechanics, and also from nonlinear partial differential equations and quantum field theory are discussed. Higher order integrals of motion for the multi-dimensional nonlinear Klein-Gordon and Schrodinger equations are mentioned.


title: Rosenthal families, pavings and generic cardinal invariants

abstract: Following D. Sobota we call a family $\mathcal F$ of infinite subsets of $\mathbb N$ a Rosenthal family if it can replace the family of all infinite subsets of $\mathbb N$ in classical Rosenthal's Lemma concerning sequences of measures on pairwise disjoint sets. We resolve two problems on Rosenthal families: every ultrafilter is a Rosenthal family and the minimal size of a Rosenthal family is exactly equal to the reaping cardinal $\mathfrak r$. This is achieved through analyzing nowhere reaping families of subsets of $\mathbb N$ and through applying a paving lemma which is a consequence of a paving lemma concerning linear operators on $\ell_1^n$ due to Bourgain. We use connections of the above results with free set results for functions on $\mathbb N$ and with linear operators on $c_0$ to determine the values of several other derived cardinal invariants.


title: Active Contact Forces Drive Non-Equilibrium Fluctuations in Membrane Vesicles

abstract: We analyze the non-equilibrium shape fluctuations of giant unilamellar vesicles encapsulating motile bacteria. Owing to bacteria--membrane collisions, we experimentally observe a significant increase in the magnitude of membrane fluctuations at low wave numbers, compared to the well-known thermal fluctuation spectrum. We interrogate these results by numerically simulating membrane height fluctuations via a modified Langevin equation, which includes bacteria--membrane contact forces. Taking advantage of the length and time scale separation of these contact forces and thermal noise, we further corroborate our results with an approximate theoretical solution to the dynamical membrane equations. Our theory and simulations demonstrate excellent agreement with non-equilibrium fluctuations observed in experiments. Moreover, our theory reveals that the fluctuation--dissipation theorem is not broken by the bacteria; rather, membrane fluctuations can be decomposed into thermal and active components.


title: Periodic coherent states decomposition and quantum dynamics on the flat torus

abstract: We provide a result on the coherent states decomposition for functions in $L^2 (\Bbb T^n)$ where $\Bbb T^n := (\Bbb R / 2\pi \Bbb Z)^n$. Suddenly, we study such a decomposition with respect to the quantum dynamics related to semiclassical elliptic Pseudodifferential operators, and we prove a related invariance result.


title: Optimizing the LO Distribution Architecture of mm-Wave Massive MIMO Receivers

abstract: Wireless networks at millimeter wavelengths have significant implementation difficulties. The path loss at these frequencies naturally leads us to consider antenna arrays with many elements. In these arrays, local oscillator (LO) generation is particularly challenging since the LO specifications affect the system architecture, signal processing design, and circuit implementation. We thoroughly analyze the effect of LO architecture design choices on the performance of a mm-wave massive MIMO uplink. This investigation focuses on the tradeoffs involved in centralized and distributed LO generation, correlated and uncorrelated phase noise sources, and the bandwidths of PLLs and carrier recovery loops. We show that, from both a performance and implementation complexity standpoint, the optimal LO architecture uses several distributed subarrays locked to a single intermediate-frequency reference in the low GHz range. Additionally, we show that the choice of PLL and carrier recovery loop bandwidths strongly affects the performance; for typical system parameters, loop bandwidths on the order of tens of MHz achieve SINRs suitable for high-order constellations. Finally, we present system simulations incorporating a complete model of the LO generation system and consider the case of a 128-element array with 16x-spatial multiplexing and a 2 GHz channel bandwidth at 75 GHz carrier. Using our optimization procedure we show that the system can support 16-way spatial multiplexing with 64-QAM modulation.


title: Non-reversible jump algorithms for Bayesian nested model selection

abstract: Non-reversible Markov chain Monte Carlo methods often outperform their reversible counterparts in terms of asymptotic variance of ergodic averages and mixing properties. Lifting the state-space (Chen et al., 1999; Diaconis et al., 2000) is a generic technique for constructing such samplers. The idea is to think of the random variables we want to generate as position variables and to associate to them direction variables so as to design Markov chains which do not have the diffusive behaviour often exhibited by reversible schemes. In this paper, we explore the benefits of using such ideas in the context of Bayesian model choice for nested models, a class of models for which the model indicator variable is an ordinal random variable. By lifting this model indicator variable, we obtain non-reversible jump algorithms, a non-reversible version of the popular reversible jump algorithms introduced by Green (1995). This simple algorithmic modification provides samplers which can empirically outperform their reversible counterparts at no extra computational cost. The code to reproduce all experiments is available online.


title: THH and traces of enriched categories

abstract: We prove that topological Hochschild homology (THH) arises from a presheaf of circles on a certain combinatorial category, which gives a universal construction of THH for any enriched infinity category. Our results rely crucially on an elementary, model-independent framework for enriched higher category theory, which may be of independent interest.


title: Improved Bounds for Burning Fence Graphs

abstract: Graph burning studies how fast a contagion, modeled as a set of fires, spreads in a graph. The burning process takes place in synchronous, discrete rounds. In each round, a fire breaks out at a vertex, and the fire spreads to all vertices that are adjacent to a burning vertex. The burning number of a graph $G$ is the minimum number of rounds necessary for each vertex of $G$ to burn. We consider the burning number of the $m \times n$ Cartesian grid graphs, written $G_{m,n}$.\ For $m = \omega(\sqrt{n})$, the asymptotic value of the burning number of $G_{m,n}$ was determined, but only the growth rate of the burning number was investigated in the case $m = O(\sqrt{n})$, which we refer to as fence graphs. We provide new explicit bounds on the burning number of fence graphs $G_{c\sqrt{n},n}$, where $c > 0$.


title: Jet Feedback and the Photon Underproduction Crisis in Simba

abstract: We examine the impact of black hole jet feedback on the properties of the low-redshift intergalactic medium (IGM) in the SIMBA simulation, with a focus on the Ly$\alpha$ forest mean flux decrement $D_A$. Without jet feedback, we confirm the Photon Underproduction Crisis (PUC) in which $\Gamma_{\rm HI}$ at $z=0$ must be increased by $\times6$ over the Haardt & Madau value in order to match the observed $D_{A}$. Turning on jet feedback lowers this discrepancy to $\sim\times 2.5$, and additionally using the recent Faucher-Gigu\`ere background mostly resolves the PUC, along with producing a flux probability distribution function in accord with observations. The PUC becomes apparent at late epochs ($z \lesssim 1$) where the jet and no-jet simulations diverge; at higher redshifts SIMBA reproduces the observed $D_{A}$ with no adjustment, with or without jets. The main impact of jet feedback is to lower the cosmic baryon fraction in the diffuse IGM from 39% to 16% at $z=0$, while increasing the warm-hot intergalactic medium (WHIM) baryon fraction from 30% to 70%; the lowering of the diffuse IGM content directly translates into a lowering of $D_{A}$ by a similar factor. Comparing to the older MUFASA simulation that employs different quenching feedback but is otherwise similar to SIMBA, MUFASA matches $D_{A}$ less well than SIMBA, suggesting that low-redshift measurements of $D_{A}$ and $\Gamma_{\rm HI}$ could provide constraints on feedback mechanisms. Our results suggest that widespread IGM heating at late times is a plausible solution to the PUC, and that SIMBA's jet AGN feedback model, included to quench massive galaxies, approximately yields this required heating.


title: Minkowski Symmetry Sets for 1-Parameter Families of Plane Curves

abstract: In this paper the generic bifurcations of the Minkowski symmetry set for 1-parameter families of plane curves are classified and the necessary and sufficient geometric criteria for each type are given. The Minkowski symmetry set is an analogue of the standard Euclidean symmetry set, and is defined to be the locus of centres of all its bitangent pseudo-circles. It is shown that the list of possible bifurcation types are different to those that occur in the list of possible types for the Euclidean symmetry set.


title: High-precision numerical solution of the Fermi polaron problem and large-order behavior of its diagrammatic series

abstract: We introduce a simple determinant diagrammatic Monte Carlo algorithm to compute the ground-state properties of a particle interacting with a Fermi sea through a zero-range interaction. The fermionic sign does not cause any fundamental problem when going to high diagram orders, and we reach order $N=30$. The data reveal that the diagrammatic series diverges exponentially as $(-1/R)^{N}$ with a radius of convergence $R<1$. Furthermore, on the polaron side of the polaron-dimeron transition, the value of $R$ is determined by a special class of three-body diagrams, corresponding to repeated scattering of the impurity between two particles of the Fermi sea. A power-counting argument explains why finite $R$ is possible for zero-range interactions in three dimensions. Resumming the divergent series through a conformal mapping yields the polaron energy with record accuracy.


title: CloudifierNet -- Deep Vision Models for Artificial Image Processing

abstract: Today, more and more, it is necessary that most applications and documents developed in previous or current technologies to be accessible online on cloud-based infrastructures. That is why the migration of legacy systems including their hosts of documents to new technologies and online infrastructures, using modern Artificial Intelligence techniques, is absolutely necessary. With the advancement of Artificial Intelligence and Deep Learning with its multitude of applications, a new area of research is emerging - that of automated systems development and maintenance. The underlying work objective that led to this paper aims to research and develop truly intelligent systems able to analyze user interfaces from various sources and generate real and usable inferences ranging from architecture analysis to actual code generation. One key element of such systems is that of artificial scene detection and analysis based on deep learning computer vision systems. Computer vision models and particularly deep directed acyclic graphs based on convolutional modules are generally constructed and trained based on natural images datasets. Due to this fact, the models will develop during the training process natural image feature detectors apart from the base graph modules that will learn basic primitive features. In the current paper, we will present the base principles of a deep neural pipeline for computer vision applied to artificial scenes (scenes generated by user interfaces or similar). Finally, we will present the conclusions based on experimental development and benchmarking against state-of-the-art transfer-learning implemented deep vision models.


title: Flat bands and gaps in twisted double bilayer graphene

abstract: We present electronic structure calculations of twisted double bilayer graphene (TDBG): A tetralayer graphene structure composed of two AB-stacked graphene bilayers with a relative rotation angle between them. Using first-principles calculations, we find that TDBG is semiconducting with a band gap that depends on the twist angle, that can be tuned by an external electric field. The gap is consistent with TDBG symmetry and its magnitude is related to surface effects, driving electron transfer from outer to inner layers. The surface effect competes with an energy upshift of localized states at inner layers, giving rise to the peculiar angle dependence of the band gap, which reduces at low angles. For these low twist angles, the TDBG develops flat bands, in which electrons in the inner layers are localized at the AA regions, as in twisted bilayer graphene.


title: Nearly Optimal Static Las Vegas Succinct Dictionary

abstract: Given a set $S$ of $n$ (distinct) keys from key space $[U]$, each associated with a value from $\Sigma$, the \emph{static dictionary} problem asks to preprocess these (key, value) pairs into a data structure, supporting value-retrieval queries: for any given $x\in [U]$, $\mathtt{valRet}(x)$ must return the value associated with $x$ if $x\in S$, or return $\bot$ if $x\notin S$. The special case where $|\Sigma|=1$ is called the \emph{membership} problem. The "textbook" solution is to use a hash table, which occupies linear space and answers each query in constant time. On the other hand, the minimum possible space to encode all (key, value) pairs is only $\mathtt{OPT}:= \lceil\lg_2\binom{U}{n}+n\lg_2|\Sigma|\rceil$ bits, which could be much less. In this paper, we design a randomized dictionary data structure using $\mathtt{OPT}+\mathrm{poly}\lg n+O(\lg\lg\lg\lg\lg U)$ bits of space, and it has \emph{expected constant} query time, assuming the query algorithm can access an external lookup table of size $n^{0.001}$. The lookup table depends only on $U$, $n$ and $|\Sigma|$, and not the input. Previously, even for membership queries and $U\leq n^{O(1)}$, the best known data structure with constant query time requires $\mathtt{OPT}+n/\mathrm{poly}\lg n$ bits of space (Pagh [Pag01] and P\v{a}tra\c{s}cu [Pat08]); the best-known using $\mathtt{OPT}+n^{0.999}$ space has query time $O(\lg n)$; the only known non-trivial data structure with $\mathtt{OPT}+n^{0.001}$ space has $O(\lg n)$ query time and requires a lookup table of size $\geq n^{2.99}$ (!). Our new data structure answers open questions by P\v{a}tra\c{s}cu and Thorup [Pat08,Tho13]. We also present a scheme that compresses a sequence $X\in\Sigma^n$ to its zeroth order (empirical) entropy up to $|\Sigma|\cdot\mathrm{poly}\lg n$ extra bits, supporting decoding each $X_i$ in $O(\lg |\Sigma|)$ expected time.


title: Mice and Covert Channels

abstract: Any secure network is only as secure as its weakest component. With overt channels tightly secured and attackers have started focusing on optical, audible, magnetic, and thermal covert channels to access sensitive systems. In this paper, we present a novel, reliable and bidirectional optical covert channel which uses optical mice. In this channel, the photocell in the mouse is used as a receiver while the LED is used as a transmitter. Our multiple experiments, which use mouse to mouse, mouse to camera and torch to mouse, show that the transmission rate can go as high as 10 bits per second. Additionally, we study the effects of infrared, distance and brightness on mouse input. We also show that infrared mice are susceptible to a similar kind of attack.


title: Invariants of models of genus one curves and modular forms

abstract: An invariant of a model of genus one curve is a polynomial in the coefficients of the model that is stable under certain linear transformations. The classical example of an invariant is the discriminant, which characterizes the singularity of models. The ring of invariants of genus one models over a field is generated by two elements. Fisher normalized these invariants for models of degree n=2,3,4 in such a way that these invariants are moreover defined over the integers. We provide an alternative way to express these normalized invariants using modular forms. This method relies on a direct computation for the discriminants based on their own geometric properties.


title: Faster Update Time for Turnstile Streaming Algorithms

abstract: In this paper, we present a new algorithm for maintaining linear sketches in turnstile streams with faster update time. As an application, we show that $\log n$ \texttt{Count} sketches or \texttt{CountMin} sketches with a constant number of columns (i.e., buckets) can be implicitly maintained in \emph{worst-case} $O(\log^{0.582} n)$ update time using $O(\log n)$ words of space, on a standard word RAM with word-size $w=\Theta(\log n)$. The exponent $0.582\approx 2\omega/3-1$, where $\omega$ is the current matrix multiplication exponent. Due to the numerous applications of linear sketches, our algorithm improves the update time for many streaming problems in turnstile streams, in the high success probability setting, without using more space, including $\ell_2$ norm estimation, $\ell_2$ heavy hitters, point query with $\ell_1$ or $\ell_2$ error, etc. Our algorithm generalizes, with the same update time and space, to maintaining $\log n$ linear sketches, where each sketch partitions the coordinates into $k<\log^{o(1)} n$ buckets using a $c$-wise independent hash function for constant $c$, and maintains the sum of coordinates for each bucket. Moreover, if arbitrary word operations are allowed, the update time can be further improved to $O(\log^{0.187} n)$, where $0.187\approx \omega/2-1$. Our update algorithm is adaptive, and it circumvents the non-adaptive cell-probe lower bounds for turnstile streaming algorithms by Larsen, Nelson and Nguy{\^{e}}n (STOC'15). On the other hand, our result also shows that proving unconditional cell-probe lower bound for the update time seems very difficult, even if the space is restricted to be (nearly) the optimum. If $\omega=2$, the cell-probe update time of our algorithm would be $\log^{o(1)} n$. Hence, proving any higher lower bound would imply $\omega>2$.


title: Learning from Explanations with Neural Execution Tree

abstract: While deep neural networks have achieved impressive performance on a range of NLP tasks, these data-hungry models heavily rely on labeled data, which restricts their applications in scenarios where data annotation is expensive. Natural language (NL) explanations have been demonstrated very useful additional supervision, which can provide sufficient domain knowledge for generating more labeled data over new instances, while the annotation time only doubles. However, directly applying them for augmenting model learning encounters two challenges: (1) NL explanations are unstructured and inherently compositional, which asks for a modularized model to represent their semantics, (2) NL explanations often have large numbers of linguistic variants, resulting in low recall and limited generalization ability. In this paper, we propose a novel Neural Execution Tree (NExT) framework to augment training data for text classification using NL explanations. After transforming NL explanations into executable logical forms by semantic parsing, NExT generalizes different types of actions specified by the logical forms for labeling data instances, which substantially increases the coverage of each NL explanation. Experiments on two NLP tasks (relation extraction and sentiment analysis) demonstrate its superiority over baseline methods. Its extension to multi-hop question answering achieves performance gain with light annotation effort.


title: Searching for pulsars in extreme orbits -- GPU acceleration of the Fourier domain 'jerk' search

abstract: Binary pulsars are an important target for radio surveys because they present a natural laboratory for a wide range of astrophysics for example testing general relativity, including detection of gravitational waves. The orbital motion of a pulsar which is locked in a binary system causes a frequency shift (a Doppler shift) in their normally very periodic pulse emissions. These shifts cause a reduction in the sensitivity of traditional periodicity searches. To correct this smearing Ransom [2001], Ransom et al. [2002] developed the Fourier domain acceleration search (FDAS) which uses a matched filtering technique. This method is however limited to a constant pulsar acceleration. Therefore, Andersen and Ransom [2018] broadened the Fourier domain acceleration search to account also for a linear change in the acceleration by implementing the Fourier domain "jerk" search into the PRESTO software package. This extension increases the number of matched filters used significantly. We have implemented the Fourier domain "jerk" search (JERK) on GPUs using CUDA. We have achieved 90x performance increase when compared to the parallel implementation of JERK in PRESTO. This work is part of the AstroAccelerate project Armour et al. [2019], a many-core accelerated time-domain signal processing library for radio astronomy.


title: Robust universal Hamiltonian quantum computing using two-body interactions

abstract: We present a new scheme to perform noise resilient universal adiabatic quantum computation using two-body interactions. To achieve this, we introduce a new family of error detecting subsystem codes whose gauge generators and a set of their logical operators - capable of encoding universal Hamiltonian computations - can be implemented using two-body interactions. Logical operators of the code are used to encode any given computational Hamiltonian, and the gauge operators are used to construct a penalty Hamiltonian whose ground subspace is protected against local-errors. In contrast to previous approaches, the constructed penalty Hamiltonian does not necessarily commute with the encoded computational Hamiltonians, but for our construction, the undesirable effect of the penalty Hamiltonian on the computation can be compensated by a simple modification of the implemented Hamiltonians. We also investigate whether a similar scheme can be constructed by encoding the computational Hamiltonian using only bare-logical operators of subsystem codes, to guarantee that computational Hamiltonian commutes with the penalty Hamiltonian. We prove a no-go theorem showing that restricting to two-body interactions and using any general CSS-type subsystem codes, such a construction cannot encode systems beyond an Ising chain in a transverse field. We point out that such a chain is universal for Hamiltonian-based quantum computation, but it is not universal for ground-state quantum computation.


title: Video-based compression for plenoptic point clouds

abstract: The plenoptic point cloud that has multiple colors from various directions, is a more complete representation than the general point cloud that usually has only one color. It is more realistic but also brings a larger volume of data that needs to be compressed efficiently. The state-of-the-art method to compress the plenoptic point cloud is an extension of the region-based adaptive hierarchical transform (RAHT). As far as we can see, in addition to RAHT, the video-based point cloud compression (V-PCC) is also an efficient point cloud compression method. However, to the best of our knowledge, no works have used a video-based solution to compress the plenoptic point cloud yet. In this paper, we first extend the V-PCC to support the plenoptic point cloud compression by generating multiple attribute videos. Then based on the observation that these videos from multiple views have very high correlations, we propose encoding them using multiview high efficiency video coding. We further propose a block-based padding method that unifies the unoccupied attribute pixels from different views to reduce their bit cost. The proposed algorithms are implemented in the V-PCC reference software. The experimental results show that the proposed algorithms can bring significant bitrate savings compared with the state-of-the-art method for plenoptic point cloud compression.


title: Fine properties of functions of bounded deformation -- an approach via linear PDEs

abstract: In this survey we collect some recent results obtained by the authors and collaborators concerning the fine structure of functions of bounded deformation (BD). These maps are $\mathrm{L}^1$-functions with the property that the symmetric part of their distributional derivative is representable as a bounded (matrix-valued) Radon measure. It has been known for a long time that for a (matrix-valued) Radon measure the property of being a symmetrized gradient can be characterized by an under-determined second-order PDE system, the Saint-Venant compatibility conditions. This observation gives rise to a new approach to the fine properties of BD-maps via the theory of PDEs for measures, which complements and partially replaces classical arguments. Starting from elementary observations, here we elucidate the ellipticity arguments underlying this recent progress and give an overview of the state of the art. We also present some open problems.


title: Time-evolution of magnetic field in hot nuclear matter with fluctuating topological charge

abstract: The time-evolution of the magnetic field in hot homogeneous nuclear matter has two qualitatively different stages separated by the sphaleron transition time $\tau_c$. At early times the axial chemical potential and the corresponding chiral conductivity $\sigma_\chi$ are slow functions of time. The soft chiral modes $k<\sigma_\chi$ of the magnetic field grow exponentially with time, which is known as the chiral instability. At later times $\sigma_\chi$ fluctuates due to the sphaleron transitions and can be regarded as a random process. It is argued that the average magnetic field is exponentially damped at later times. The time-evolution of the average magnetic energy is more complicated and depends on the electrical conductivity of the chiral matter but does not depend on chirality. It exhibits instability only if the matter is a poor electrical conductor, such as the quark-gluon plasma near the critical temperature. The precise conditions for the instability and the growth rate of the unstable modes are derived.


title: Random Forest Models for Accurate Identification of Coordination Environments from X-ray Absorption Near-Edge Structure

abstract: Analyzing coordination environments using X-ray absorption spectroscopy has broad applications ranging from solid-state physics to material chemistry. Here, we show that random forest models can identify the main coordination environment from K-edge X-ray absorption near edge structure (XANES) with a high accuracy of 85.4% and all associated coordination environments with a high Jaccard score of 81.8% for 33 cation elements in oxides, significantly outperforming other machine learning (ML) models. In a departure from prior works, we used a robust description of the coordination environment as a distribution over 25 distinct coordination motifs with coordination numbers ranging from 1-12. The random forest models were trained on the world's largest database of ~190,000 computed K-edge XANES spectra. Furthermore, the random forest models can be applied to predict the coordination environment from experimental K-edge XANES with minimal loss in accuracy (82.1%) due to the use of data augmentation. A drop-out feature importance analysis highlights the key roles that the pre-edge and main-peak regions play in coordination environment identification, with the post-peak region becoming increasingly important at higher coordination numbers. This work provides a general strategy to identify the coordination environment from K-edge XANES across broad chemistries, paving the way for future advancements in the application of ML to spectroscopy.


title: Isospin-Violating Contributions to $\epsilon'/\epsilon$

abstract: The known isospin-breaking contributions to the $K\rightarrow \pi\pi$ amplitudes are reanalyzed, taking into account our current understanding of the quark masses and the relevant non-perturbative inputs. We present a complete numerical reappraisal of the direct CP-violating ratio $\epsilon'/\epsilon$, where these corrections play a quite significant role. We obtain the Standard Model prediction $\text{Re}\left(\epsilon'/\epsilon\right)\, =\,\left(14\,\pm\,5\right)\cdot 10^{-4}$, which is in very good agreement with the measured ratio. The uncertainty, which has been estimated conservatively, is dominated by our current ignorance about $1/N_C$-suppressed contributions to some relevant chiral-perturbation-theory low-energy constants.


title: The possible new force in nature

abstract: The {\it zitterbewegung region} is studied in this paper and we investigate the possible presence of a new force in this region.This force was first proposed by Sidharth in \cite{cu-bgs}. This article intends to study the existence of a fifth force. Here, essentially, equations that characterize new particles and their antiparticles are obtained. In this manner a pedestal is built for a field theoretic approach that will be the subject of the second part of this theory.


title: The next decade of black hole spectroscopy

abstract: Gravitational wave observations of the ringdown of the remnant black hole in a binary black hole coalescence provide a unique opportunity of confronting the black hole no-hair theorem in general relativity with observational data. The most robust tests are possible if multiple ringdown modes can be observed. In this paper, using state-of-the-art Bayesian inference methods and the most up-to-date knowledge of binary black hole population parameters and ringdown mode amplitudes, we evaluate the prospects for black hole spectroscopy with current and future ground based gravitational wave detectors over the next 10 years. For different population models, we estimate the likely number of events for which the subdominant mode can be detected and distinguished from the dominant mode. We show that black hole spectroscopy could significantly test general relativity for events seen by the proposed LIGO Voyager detectors.


title: What's happening in traditional and inquiry-based introductory labs? An integrative analysis at a large research university

abstract: There is a growing recognition of the need to replace "cookbook"-style introductory labs with more-meaningful learning experiences. To identify the strengths and weaknesses of a mix of cookbook-style and inquiry-based labs, an introductory lab course currently being reformed was observed following a reflexive ethnographic protocol and pre and post E-CLASS surveys were administered. We analyzed data to identify shortcomings of the current labs and to determine areas for improvement.


title: All aboard! Challenges and successes in professional development for physics lab TAs

abstract: At large research universities in the USA, introductory physics labs are often run by graduate student teaching assistants (TAs). Thus, efforts to reform introductory labs should address the need for effective and relevant TA professional development. We developed and implemented a research-based professional development program that focuses on preparing TAs to effectively support inquiry-based learning in the lab. We identify positive effects by examining three possible ways in which the professional development might have impacted TAs and their work. First, we examine lab TAs' written reflections to understand the effect of the program on TAs' ways of thinking about student learning. Second, we observe and categorize TA-student interactions in the lab in order to investigate whether TA behaviors are changing after the professional development. Third, we examine students' attitudes toward experimental science and present one example case in which students' attitudes improve for those TAs who `buy in' to the professional development. Our results suggest lab TA professional development may have a tangible positive impact on TA performance and student learning.


title: Modeling of the drag and thermal conductivity coefficients in tropical cyclones taking into account foam cover

abstract: A formal averaging procedure over the air-sea interface is developed for both momentum and enthalpy surface-transfer coefficients, C_D and C_K, in hurricane conditions. This leads to splitting of both the transfer coefficients across the total area of the sea surface, C_(D,K), into the sums of their partial values over the foam-free, C_(D,Kw), and foam-covered, C_(D,Kf), fractions weighted with foam- and water-coverage coefficients {\alpha}f and 1-{\alpha}f . The transfer coefficients, C_(D,K) and C_(D,Kw), are estimated by measurements in open-sea and in laboratory conditions, respectively, while the transfer coefficients across the foam-covered fraction, C_(D,Kf), which cannot be measured directly, are estimated from the splitting relations. Applying the Monin-Obukhov similarity theory at the neutral stability atmospheric conditions to the transfer coefficients, separately to the foam-free, the foam-covered, and the total sea surfaces, yields the roughness lengths, Z_(D,Kw),Z_(D,Kf),Z_(D,K). The study is aimed at explaining an anomalous behavior of the momentum and enthalpy transfer coefficients, C_D and C_K, with wind speed U10 under hurricane conditions, by the effect of the foam slipping layer sandwiched between the atmosphere and the ocean.


title: Interatomic Potential in a Simple Dense Neural Network Representation

abstract: Simulations at the atomic scale provide a direct and effective way to understand the mechanical properties of materials. In the regime of classical mechanics, simulations for the thermodynamic properties of metals and alloys can be done by either solving the equations of motion or performing Monte Carlo sampling. The key component for an accurate simulation of such physical systems to produce faithful physical quantities is the use of an appropriate potential or a force field. In this paper, we explore the use of methods from the realm of machine learning to overcome and bypass difficulties encountered when fitting potentials for atomic systems. Particularly, we will show that classical potentials can be represented by a dense neural network with good accuracy.


title: Framework for Inferring Following Strategies from Time Series of Movement Data

abstract: How do groups of individuals achieve consensus in movement decisions? Do individuals follow their friends, the one predetermined leader, or whomever just happens to be nearby? To address these questions computationally, we formalize "Coordination Strategy Inference Problem". In this setting, a group of multiple individuals moves in a coordinated manner towards a target path. Each individual uses a specific strategy to follow others (e.g. nearest neighbors, pre-defined leaders, preferred friends). Given a set of time series that includes coordinated movement and a set of candidate strategies as inputs, we provide the first methodology (to the best of our knowledge) to infer whether each individual uses local-agreement-system or dictatorship-like strategy to achieve movement coordination at the group level. We evaluate and demonstrate the performance of the proposed framework by predicting the direction of movement of an individual in a group in both simulated datasets as well as two real-world datasets: a school of fish and a troop of baboons. Moreover, since there is no prior methodology for inferring individual-level strategies, we compare our framework with the state-of-the-art approach for the task of classification of group-level-coordination models. The results show that our approach is highly accurate in inferring the correct strategy in simulated datasets even in complicated mixed strategy settings, which no existing method can infer. In the task of classification of group-level-coordination models, our framework performs better than the state-of-the-art approach in all datasets. Animal data experiments show that fish, as expected, follow their neighbors, while baboons have a preference to follow specific individuals. Our methodology generalizes to arbitrary time series data of real numbers, beyond movement data.

categories: stat.ML cs.AI cs.LG cs.MA physics.data-an


title: Stochastic model for the CheY-P molarity in the neighbourhood of E. coli flagella motors

abstract: E.coli serves as prototype for the study of peritrichous enteric bacteria that perform runs and tumbles alternately. Bacteria run forward as a result of the counterclockwise (CCW) rotation of their flagella bundle and perform tumbles when at least one of their flagella rotates clockwise (CW), moving away from the bundle. The flagella are hooked to molecular rotary motors of nanometric diameter able to make transitions between CCW and CW rotations that last up to one hundredth of a second. At the same time, flagella move or rotate the bacteria's body microscopically during lapses that range between a tenth and ten seconds. We assume that the transitions between CCW and CW rotations occur solely by fluctuations of CheY-P molarity in the presence of two threshold values, and that a veto rule selects the run or tumble motions. We present Langevin eqs for the CheY-P molarity in the vicinity of each molecular motor. This model allows to obtain the run- or tumble-time distribution as a linear combination of decreasing exponentials that is a function of the steady molarity of CheY-P in the neighbourhood of the molecular motor, which fits experimental data. In turn, if the internal signaling system is unstimulated, we show that the runtime distributions reach power-law behaviour, a characteristic of self-organized systems, in some time range and, afterwards, exponential cutoff. In addition, our model explains without any fitting parameters the ultrasensitivity of the flagella motors as a function of the steady state of CheY-P molarity. In addition, we show that the tumble bias for peritrichous bacterium has a similar sigmoid-shape to the CW bias, although shifted to lower concentrations when the flagella number increases. Thus, the increment in the flagella number allows lower operational values for each motor increasing amplification and robustness of the chemotatic pathway.


title: Instrumentation for solar spectropolarimetry: state of the art and prospects

abstract: Given its unchallenged capabilities in terms of sensitivity and spatial resolution, the combination of imaging spectropolarimetry and numeric Stokes inversion represents the dominant technique currently used to remotely sense the physical properties of the solar atmosphere and, in particular, its important driving magnetic field. Solar magnetism manifests itself in a wide range of spatial, temporal, and energetic scales. The ubiquitous but relatively small and weak fields of the so-called quiet Sun are believed today to be crucial for answering many open questions in solar physics, some of which have substantial practical relevance due to the strong Sun-Earth connection. However, such fields are very challenging to detect because they require spectropolarimetric measurements with high spatial (sub-arcsec), spectral (<100 mA), and temporal (<10 s) resolution along with high polarimetric sensitivity (<0.001 of the intensity). We collect and discuss both well-established and upcoming instrumental solutions developed during the last decades to push solar observations toward the above-mentioned parameter regime. This typically involves design trade-offs due to the high dimensionality of the data and signal-to-noise-ratio considerations, among others. We focus on the main three components that form a spectro-polarimeter, namely, wavelength discriminators, the devices employed to encode the incoming polarization state into intensity images (polarization modulators), and the sensor technologies used to register them. We consider the instrumental solutions introduced to perform this kind of measurements at different optical wavelengths and from various observing locations, i.e., ground-based, from the stratosphere or near space.


title: Quantum fluctuations of the center-of-mass and relative parameters of NLS breathers

abstract: We study quantum fluctuations of macroscopic parameters of an NLS breathers, i.e., the second-order soliton solution of the nonlinear Schr\"odinger equation. Uncertainty relations for the parameters are derived and compared to similar relations for fundamental solitons. We compare two models for the state of the quantum field of fluctuations surrounding the classical field of the Bose-Einstein condensate: a conventionally used, computationally convenient "white noise", and a correlated noise which assumes that the breather has been created from a fundamental soliton, by means of the application of the factor-of-four quench of the nonlinearity strength. Theoretical methods used in the work are well suited for a large number of particles, $N$. We thus confirm the possibility of experimental observation of macroscopic quantum fluctuations, which is suggested by an extrapolation to large $N$ of recently reported low-$N$ Bethe-ansatz results [Phys. Rev. Lett. 119 220401 (2017)].


title: Self-Supervised Difference Detection for Weakly-Supervised Semantic Segmentation

abstract: To minimize the annotation costs associated with the training of semantic segmentation models, researchers have extensively investigated weakly-supervised segmentation approaches. In the current weakly-supervised segmentation methods, the most widely adopted approach is based on visualization. However, the visualization results are not generally equal to semantic segmentation. Therefore, to perform accurate semantic segmentation under the weakly supervised condition, it is necessary to consider the mapping functions that convert the visualization results into semantic segmentation. For such mapping functions, the conditional random field and iterative re-training using the outputs of a segmentation model are usually used. However, these methods do not always guarantee improvements in accuracy; therefore, if we apply these mapping functions iteratively multiple times, eventually the accuracy will not improve or will decrease. In this paper, to make the most of such mapping functions, we assume that the results of the mapping function include noise, and we improve the accuracy by removing noise. To achieve our aim, we propose the self-supervised difference detection module, which estimates noise from the results of the mapping functions by predicting the difference between the segmentation masks before and after the mapping. We verified the effectiveness of the proposed method by performing experiments on the PASCAL Visual Object Classes 2012 dataset, and we achieved 64.9\% in the val set and 65.5\% in the test set. Both of the results become new state-of-the-art under the same setting of weakly supervised semantic segmentation.


title: Affective Behaviour Analysis of On-line User Interactions: Are On-line Support Groups more Therapeutic than Twitter?

abstract: The increase in the prevalence of mental health problems has coincided with a growing popularity of health related social networking sites. Regardless of their therapeutic potential, On-line Support Groups (OSGs) can also have negative effects on patients. In this work we propose a novel methodology to automatically verify the presence of therapeutic factors in social networking websites by using Natural Language Processing (NLP) techniques. The methodology is evaluated on On-line asynchronous multi-party conversations collected from an OSG and Twitter. The results of the analysis indicate that therapeutic factors occur more frequently in OSG conversations than in Twitter conversations. Moreover, the analysis of OSG conversations reveals that the users of that platform are supportive, and interactions are likely to lead to the improvement of their emotional state. We believe that our method provides a stepping stone towards automatic analysis of emotional states of users of online platforms. Possible applications of the method include provision of guidelines that highlight potential implications of using such platforms on users' mental health, and/or support in the analysis of their impact on specific individuals.


title: A new simple proof for Lum-Chua's conjecture

abstract: The already proved Lum-Chua's conjecture says that a continuous planar piecewise linear differential system with two zones separated by a straight line has at most one limit cycle. In this paper, we provide a new proof by using a novel characterization for Poincar\'e half-maps in planar linear systems. This proof is very short and straightforward, because this characterization avoids the inherent flaws of the usual methods to study piecewise linear systems (the appearance of large case-by-case analysis due to the unnecessary discrimination between the different spectra of the involved matrices, to deal with transcendental equations forced by the implicit occurrence of flight time, ...). In addition, the application of the characterization allow us to prove that if a limit cycle exists, then it is hyperbolic and its stability is determined by a simple relationship between the parameters. To the best of our knowledge, the hyperbolicity of the limit cycle and this simple expression for its stability have not been pointed out before.


title: Gradient-based Adaptive Markov Chain Monte Carlo

abstract: We introduce a gradient-based learning method to automatically adapt Markov chain Monte Carlo (MCMC) proposal distributions to intractable targets. We define a maximum entropy regularised objective function, referred to as generalised speed measure, which can be robustly optimised over the parameters of the proposal distribution by applying stochastic gradient optimisation. An advantage of our method compared to traditional adaptive MCMC methods is that the adaptation occurs even when candidate state values are rejected. This is a highly desirable property of any adaptation strategy because the adaptation starts in early iterations even if the initial proposal distribution is far from optimum. We apply the framework for learning multivariate random walk Metropolis and Metropolis-adjusted Langevin proposals with full covariance matrices, and provide empirical evidence that our method can outperform other MCMC algorithms, including Hamiltonian Monte Carlo schemes.


title: Algorithms for Intersection Graphs of Multiple Intervals and Pseudo Disks

abstract: Intersection graphs of planar geometric objects such as intervals, disks, rectangles and pseudo-disks are well studied. Motivated by various applications, Butman et al. in SODA 2007 considered algorithmic questions in intersection graphs of $t$-intervals. A $t$-interval is a union of at most $t$ distinct intervals (here $t$ is a parameter) -- these graphs are referred to as Multiple-Interval Graphs. Subsequent work by Kammer et al. in Approx 2010 also considered $t$-disks and other geometric shapes. In this paper we revisit some of these algorithmic questions via more recent developments in computational geometry. For the minimum weight dominating set problem, we give a simple $O(t \log t)$ approximation for Multiple-Interval Graphs, improving on the previously known bound of $t^2$ . We also show that it is NP-hard to obtain an $o(t)$-approximation in this case. In fact, our results hold for the intersection graph of a set of t-pseudo-disks which is a much larger class. We obtain an ${\Omega}(1/t)$-approximation for the maximum weight independent set in the intersection graph of $t$-pseudo-disks. Our results are based on simple reductions to existing algorithms by appropriately bounding the union complexity of the objects under consideration.


title: Highly nonlinear ionization of atoms induced by intense high-harmonic pulses

abstract: Intense extreme-ultraviolet (XUV) pulses enable the investigation of XUV-induced nonlinear processes and are a prerequisite for the development of attosecond pump - attosecond probe experiments. While highly nonlinear processes in the XUV range have been studied at free-electron lasers (FELs), high-harmonic generation (HHG) has allowed the investigation of low-order nonlinear processes. Here we suggest a concept to optimize the HHG intensity, which surprisingly requires a scaling of the experimental parameters that differs substantially from optimizing the HHG pulse energy. As a result, we are able to study highly nonlinear processes in the XUV range using a driving laser with a modest ($\approx 10$ mJ) pulse energy. We demonstrate our approach by ionizing Ar atoms up to Ar$^{5+}$, requiring the absorption of at least 10 XUV photons.


title: CANet: Cross-disease Attention Network for Joint Diabetic Retinopathy and Diabetic Macular Edema Grading

abstract: Diabetic retinopathy (DR) and diabetic macular edema (DME) are the leading causes of permanent blindness in the working-age population. Automatic grading of DR and DME helps ophthalmologists design tailored treatments to patients, thus is of vital importance in the clinical practice. However, prior works either grade DR or DME, and ignore the correlation between DR and its complication, i.e., DME. Moreover, the location information, e.g., macula and soft hard exhaust annotations, are widely used as a prior for grading. Such annotations are costly to obtain, hence it is desirable to develop automatic grading methods with only image-level supervision. In this paper, we present a novel cross-disease attention network (CANet) to jointly grade DR and DME by exploring the internal relationship between the diseases with only image-level supervision. Our key contributions include the disease-specific attention module to selectively learn useful features for individual diseases, and the disease-dependent attention module to further capture the internal relationship between the two diseases. We integrate these two attention modules in a deep network to produce disease-specific and disease-dependent features, and to maximize the overall performance jointly for grading DR and DME. We evaluate our network on two public benchmark datasets, i.e., ISBI 2018 IDRiD challenge dataset and Messidor dataset. Our method achieves the best result on the ISBI 2018 IDRiD challenge dataset and outperforms other methods on the Messidor dataset. Our code is publicly available at https://github.com/xmengli999/CANet.


title: The $\textit{u}$-series: A separable decomposition for electrostatics computation with improved accuracy

abstract: The evaluation of electrostatic energy for a set of point charges in a periodic lattice is a computationally expensive part of molecular dynamics simulations (and other applications) because of the long-range nature of the Coulomb interaction. A standard approach is to decompose the Coulomb potential into a near part, typically evaluated by direct summation up to a cutoff radius, and a far part, typically evaluated in Fourier space. In practice, all decomposition approaches involve approximations---such as cutting off the near-part direct sum---but it may be possible to find new decompositions with improved tradeoffs between accuracy and performance. Here we present the $\textit{u-series}$, a new decomposition of the Coulomb potential that is more accurate than the standard (Ewald) decomposition for a given amount of computational effort, and achieves the same accuracy as the Ewald decomposition with approximately half the computational effort. These improvements, which we demonstrate numerically using a lipid membrane system, arise because the $\textit{u}$-series is smooth on the entire real axis and exact up to the cutoff radius. Additional performance improvements over the Ewald decomposition may be possible in certain situations because the far part of the $\textit{u}$-series is a sum of Gaussians, and can thus be evaluated using algorithms that require a separable convolution kernel; we describe one such algorithm that reduces communication latency at the expense of communication bandwidth and computation, a tradeoff that may be advantageous on modern massively parallel supercomputers.


title: The liquid argon detector and measurement of SiPM array at liquid argon temperature

abstract: Particle detectors based on liquid argon (LAr) have recently become recognized as an extremely attractive technology for the direct detection of dark matter as well as the measurement of coherent elastic neutrino-nucleus scattering (CE$\nu$NS). The Chinese argon group at Institute of High Energy Physics has been studying the LAr detector technology and a LAr detector has been operating steadily. A program of using a dual phase LAr detector to measure the CE$\nu$NS at Taishang Nuclear Power Plant has been proposed and the R\&D work is ongoing. Considering the requirements of ultra-low radio-purity and high photon collection efficiency, SiPMs will be a good choice and will be used in the detector. In this proceeding, an introduction of the LAr detector and the measurement results of SiPM array at LAr temperature will be presented.


title: Constraining short gamma-ray burst jet properties with gravitational waves and gamma rays

abstract: Gamma-ray burst (GRB) prompt emission is highly beamed, and understanding the jet geometry and beaming configuration can provide information on the poorly understood central engine and circum-burst environment. Prior to the advent of gravitational-wave astronomy, astronomers relied on observations of jet breaks in the multi-wavelength afterglow to determine the GRB opening angle, since the observer's viewing angle relative to the system cannot be determined from the electromagnetic data alone. Gravitational-wave observations, however, provide an independent measurement of the viewing angle. We describe a Bayesian method for determining the geometry of short GRBs using coincident electromagnetic and gravitational-wave observations. We demonstrate how an ensemble of multi-messenger detections can be used to measure the distributions of the jet energy, opening angle, Lorentz factor, and angular profile of short GRBs; we find that for a population of 100 such observations, we can constrain the mean of the opening angle distribution to within 10 degrees regardless of the angular emission profile. Conversely, the constraint on the energy distribution depends on the shape of the profile, which can be distinguished.


title: A Decentralized Time- and Energy-Optimal Control Framework for Connected Automated Vehicles: From Simulation to Field Test

abstract: The implementation of connected and automated vehicle (CAV) technologies enables a novel computational framework for real-time control aimed at optimizing energy consumption with associated benefits. In this paper, we implement an optimal control framework, developed previously, in an Audi A3 etron plug-in hybrid electric vehicle, and demonstrate that we can improve the vehicle's efficiency and travel time in a corridor including an on-ramp merging, a speed reduction zone, and a roundabout. Our exposition includes the development, integration, implementation and validation of the proposed framework in (1) simulation, (2) hardware-in-the-loop (HIL) testing, (3) connectivity enabled virtual reality based bench-test, and (4) field test in Mcity. We show that by adopting such inexpensive, yet effective process, we can efficiently integrate and test the controller framework, ensure proper connectivity and data transmission between different modules of the system, and reduce uncertainty. We evaluate the performance and effectiveness of the control framework and observe significant improvement in terms of energy and travel time compared to the baseline scenario.


title: Bare quantum simultaneity versus classical interactivity in communication complexity

abstract: A relational bipartite communication problem is presented that has an efficient quantum simultaneous-messages protocol, but no efficient classical two-way protocol.


title: Amortized Population Gibbs Samplers with Neural Sufficient Statistics

abstract: We develop amortized population Gibbs (APG) samplers, a class of scalable methods that frames structured variational inference as adaptive importance sampling. APG samplers construct high-dimensional proposals by iterating over updates to lower-dimensional blocks of variables. We train each conditional proposal by minimizing the inclusive KL divergence with respect to the conditional posterior. To appropriately account for the size of the input data, we develop a new parameterization in terms of neural sufficient statistics. Experiments show that APG samplers can train highly structured deep generative models in an unsupervised manner, and achieve substantial improvements in inference accuracy relative to standard autoencoding variational methods.


title: On the performance of particle filters with adaptive number of particles

abstract: We investigate the performance of a class of particle filters (PFs) that can automatically tune their computational complexity by evaluating online certain predictive statistics which are invariant for a broad class of state-space models. To be specific, we propose a family of block-adaptive PFs based on the methodology of Elvira et al (2017). In this class of algorithms, the number of Monte Carlo samples (known as particles) is adjusted periodically, and we prove that the theoretical error bounds of the PF actually adapt to the updates in the number of particles. The evaluation of the predictive statistics that lies at the core of the methodology is done by generating fictitious observations, i.e., particles in the observation space. We study, both analytically and numerically, the impact of the number $K$ of these particles on the performance of the algorithm. In particular, we prove that if the predictive statistics with $K$ fictitious observations converged exactly, then the particle approximation of the filtering distribution would match the first $K$ elements in a series of moments of the true filter. This result can be understood as a converse to some convergence theorems for PFs. From this analysis, we deduce an alternative predictive statistic that can be computed (for some models) without sampling any fictitious observations at all. Finally, we conduct an extensive simulation study that illustrates the theoretical results and provides further insights into the complexity, performance and behavior of the new class of algorithms.


title: Statistical properties of a granular gas fluidized by turbulent air wakes

abstract: We perform experiments with a granular system that consists of a collection of identical hollow spheres (ping-pong balls). Particles rest on a horizontal metallic grid and are confined within a circular region. Fluidization is achieved by means of a turbulent air current coming from below. Air flow is adjusted so that the balls do not elevate over the grid, as an approach to 2D dynamics. With a high-speed camera, we take images of the system. From these images we can infer horizontal particle positions and velocities by means of particle-tracking algorithms. With the obtained data we analyze: a) the systematic measurement error in the determination of positions and velocities from our digital images; b) the degree of homogeneity achieved in our experiments (which depends on possible deviations of the grid from the horizontal and on the homogeneity of turbulent air wakes). Interestingly, we have observed evidences of crystallization at high enough densities.


title: "Predicting" after peeking into the future: Correcting a fundamental flaw in the SAOM -- TERGM comparison of Leifeld and Cranmer (2019)

abstract: We review the empirical comparison of SAOMs and TERGMs by Leifeld and Cranmer (2019) in Network Science. We note that their model specification uses nodal covariates calculated from observed degrees instead of using structural effects, thus turning endogeneity into circularity. In consequence, their out-of-sample predictions using TERGMs are based on out-of-sample information and thereby predict the future using observations from the future. We conclude that their analysis rest on erroneous model specifications that render the article's conclusions meaningless. Consequently, researchers should disregard recommendations from the criticized paper when making informed modelling choices.


title: On the derivatives of the powers of trigonometric and hyperbolic sine and cosine

abstract: This work contains different expressions for the k'th derivative of the n'th power of the trigonometric and hyperbolic sine and cosine. The first set of expressions follow from the complex definitions of the trigonometric and hyperbolic sine and cosine, and the binomial theorem. The other expressions are polynomial-based. They are perhaps less obvious, and use only polynomials in sin(x) and cos(x), or in sinh(x) and cosh(x). No sines or cosines of arguments other than x appear in these polynomial-based expressions. The final expressions are dependent only on sin(x), cos(x), sinh(x), or cosh(x) respectively when k is even; and they only have a single additional factor cos(x), sin(x), cosh(x), or sinh(x) respectively when k is odd.