title: Simulation-Based Inference for Global Health Decisions

abstract: The COVID-19 pandemic has highlighted the importance of in-silico epidemiological modelling in predicting the dynamics of infectious diseases to inform health policy and decision makers about suitable prevention and containment strategies. Work in this setting involves solving challenging inference and control problems in individual-based models of ever increasing complexity. Here we discuss recent breakthroughs in machine learning, specifically in simulation-based inference, and explore its potential as a novel venue for model calibration to support the design and evaluation of public health interventions. To further stimulate research, we are developing software interfaces that turn two cornerstone COVID-19 and malaria epidemiology models COVID-sim, (https://github.com/mrc-ide/covid-sim/) and OpenMalaria (https://github.com/SwissTPH/openmalaria) into probabilistic programs, enabling efficient interpretable Bayesian inference within those simulators.


title: New model of chromite and magnesiochromite solubility in silicate melts

abstract: Chromite is a key magmatic mineral frequently used as petrogenetic indicator of physico-chemical conditions of mafic magma crystallization. In this work, magnesiochromite and chromite solubility in a natural basalt and an iron-free haplobasalt was investigated at 1440{\deg}C and atmospheric pressure under controlled CO-CO2 gas mixtures corresponding to the range two log units below to two log units above the fayalite-magnetite-quartz buffer. The source of chromium was either natural chromite or synthetic Cr2O3, the latter reacting with the basaltic liquids to form a magnesiochromite. The highest concentrations of Cr in haplobasaltic melts are found in equilibrium with magnesiochromite, depending on redox conditions. In detail, at low fO2, liquids have high Cr contents, but the variation of log [Cr, ppm] is not a linear function of log fO2. Using our new data and data from the literature a model for Cr concentrations at chromite/magnesiochromite saturation in silicate melts has been developed based upon a thermodynamic formalism. Our model may be used to assess the effect of melt composition on chromite and magnesiochromite solubility in silicate melts during peridotite melting and assimilation. At moderate and high oxygen fugacities, concentration levels of Cr at chromite saturation are higher for ultramafic than for felsic rocks. Our data imply that assimilation of magnesiochromite-bearing serpentinite lithosphere could result in high Cr contents in mafic melts, triggering massive crystallization of chromite, especially upon the system hydridization, cooling, oxidation and magma degassing. Our model may be applied for quantitative prediction of chromite crystallization and formation of the terrestrial mantle-crust transition zones and mantle represented by chromitites and dunites.


title: Multi-agent Communication meets Natural Language: Synergies between Functional and Structural Language Learning

abstract: We present a method for combining multi-agent communication and traditional data-driven approaches to natural language learning, with an end goal of teaching agents to communicate with humans in natural language. Our starting point is a language model that has been trained on generic, not task-specific language data. We then place this model in a multi-agent self-play environment that generates task-specific rewards used to adapt or modulate the model, turning it into a task-conditional language model. We introduce a new way for combining the two types of learning based on the idea of reranking language model samples, and show that this method outperforms others in communicating with humans in a visual referential communication task. Finally, we present a taxonomy of different types of language drift that can occur alongside a set of measures to detect them.


title: Optical Response and Drift Matrix of Quadratic Optomechanical System

abstract: Nonlinear interactions in optomechanical systems play a crucial role in many emerging number of interesting studies and phenomena such as existence of optomechanical chaos introduced by F. Monifi et al. [Nature Photonics 10, 399405 (2016)] and optomechanical symmetry breaking proposed by Zhong-Peng Liu et al. [Phys. Rev. Lett.117, 110802 (2016)]. In this article we have theoretically examined quadratically coupled optomechanical system containing two atomic levels. We have first studied the solution of various modes of the system at steady state and later we have observed the variation of Transmission Intensity (T) with several parameters of the system. Further we have extended our analyzation to find Drift matrix of the quadratic optomechanical system and stability conditions by adiabetically eliminating atomic degree of freedom.


title: Quantifying the Redshift Space Distortion of the Bispectrum II: Induced Non-Gaussianity at Second Order Perturbation

abstract: The anisotrpy of the redshift space bispectrum $B^s(\mathbf{k_1},\mathbf{k_2},\mathbf{k_3})$, which contains a wealth of cosmological information, is completely quantified using multipole moments $\bar{B}^m_{\ell}(k_1,\mu,t)$ where $k_1$, the length of the largest side, and $(\mu,t)$ respectively quantify the size and shape of the triangle $(\mathbf{k_1},\mathbf{k_2},\mathbf{k_3})$. We present analytical expressions for all the multipoles which are predicted to be non-zero ($\ell \le 8, m \le 6$ ) at second order perturbation theory. The multipoles also depend on $\beta_1,b_1$ and $\gamma_2$, which quantify the linear redshift distortion parameter, linear bias and quadratic bias respectively. Considering triangles of all possible shapes, we analyse the shape dependence of all of the multipoles holding $k_1=0.2 \, {\rm Mpc}^{-1}, \beta_1=1, b_1=1$ and $\gamma_2=0$ fixed. The monopole $\bar{B}^0_0$, which is positive everywhere, is minimum for equilateral triangles. $\bar{B}_0^0$ increases towards linear triangles, and is maximum for linear triangles close to the squeezed limit. Both $\bar{B}^0_{2}$ and $\bar{B}^0_4$ are similar to $\bar{B}^0_0$, however the quadrupole $\bar{B}^0_2$ exceeds $\bar{B}^0_0$ over a significant range of shapes. The other multipoles, many of which become negative, have magnitudes smaller than $\bar{B}^0_0$. In most cases the maxima or minima, or both, occur very close to the squeezed limit. $\mid \bar{B}^m_{\ell} \mid $ is found to decrease rapidly if $\ell$ or $m$ are increased. The shape dependence shown here is characteristic of non-linear gravitational clustering. Non-linear bias, if present, will lead to a different shape dependence.


title: Existence and Uniqueness of Recursive Utility Models in $L_p$

abstract: Recursive preferences, of the sort developed by Epstein and Zin (1989), play an integral role in modern macroeconomics and asset pricing theory. Unfortunately, it is non-trivial to establish the unique existence of a solution to recursive utility models. We show that the tightest known existence and uniqueness conditions can be extended to (i) Schorfheide, Song and Yaron (2018) recursive utilities and (ii) recursive utilities with `narrow framing'. Further, we sharpen the solution space of Borovicka and Stachurski (2019) from $L_1$ to $L_p$ so that the results apply to a broader class of modern asset pricing models. For example, using $L_2$ Hilbert space theory, we find the class of parameters which generate a unique $L_2$ solution to the Bansal and Yaron (2004) and Schorfheide, Song and Yaron (2018) models.


title: Recognition of 26 Degrees of Freedom of Hands Using Model-based approach and Depth-Color Images

abstract: In this study, we present an model-based approach to recognize full 26 degrees of freedom of a human hand. Input data include RGB-D images acquired from a Kinect camera and a 3D model of the hand constructed from its anatomy and graphical matrices. A cost function is then defined so that its minimum value is achieved when the model and observation images are matched. To solve the optimization problem in 26 dimensional space, the particle swarm optimization algorimth with improvements are used. In addition, parallel computation in graphical processing units (GPU) is utilized to handle computationally expensive tasks. Simulation and experimental results show that the system can recognize 26 degrees of freedom of hands with the processing time of 0.8 seconds per frame. The algorithm is robust to noise and the hardware requirement is simple with a single camera.


title: On Learned Operator Correction in Inverse Problems

abstract: We discuss the possibility to learn a data-driven explicit model correction for inverse problems and whether such a model correction can be used within a variational framework to obtain regularised reconstructions. This paper discusses the conceptual difficulty to learn such a forward model correction and proceeds to present a possible solution as forward-adjoint correction that explicitly corrects in both data and solution spaces. We then derive conditions under which solutions to the variational problem with a learned correction converge to solutions obtained with the correct operator. The proposed approach is evaluated on an application to limited view photoacoustic tomography and compared to the established framework of Bayesian approximation error method.

categories: math.NA cs.CV cs.LG cs.NA eess.IV math.OC


title: On complex dynamics in a Purkinje and a ventricular cardiac cell model

abstract: Cardiac muscle cells can exhibit complex patterns including irregular behaviour such as chaos or (chaotic) early afterdepolarisations (EADs), which can lead to sudden cardiac death. Suitable mathematical models and their analysis help to predict the occurrence of such phenomena and to decode their mechanisms. The focus of this paper is the investigation of dynamics of cardiac muscle cells described by systems of ordinary differential equations. This is generically performed by studying a Purkinje cell model and a modified ventricular cell model. We find chaotic dynamics with respect to the leak current in the Purkinje cell model, and EADs and chaos with respect to a reduced fast potassium current and an enhanced calcium current in the ventricular cell model -- features that have been experimentally observed and are known to exist in some models, but are new to the models under present consideration. We also investigate the related monodomain models of both systems to study synchronisation and the behaviour of the cells on macro-scale in connection with the discovered features. The models show qualitatively the same behaviour to what has been experimentally observed. However, for certain parameter settings the dynamics occur within a non-physiological range.


title: Modeling and forecasting of the COVID-19 pandemic in India

abstract: Background: The ongoing COVID-19 epidemic dilated rapidly throughout India. To end the global COVID-19 pandemic major behavioral, social distancing, contact tracing, and state interventions has been undertaken to reduce the outbreak and avert the persistence of the coronavirus in humans in India and worldwide. In absence of any vaccine or therapeutics, forecasting is of utmost priority for health care planning and control the transmission of COVID-19. Methods: We have proposed a mathematical model that explain the transmission dynamics of COVID-19 in India. Based on the estimated data our model predicts the evolution of epidemics and the end of SARS-CoV-2 and aids to evaluate the influence of different policies to control the spread of the diseases. Findings: With the real data for infected individuals, we find the basic reproduction number $R_0$ for 17 states of India and overall India. A complete figure is given to demonstrate the estimated pandemic life cycle along with the real data or history to date. Our study reveals that the strict control measures implemented in India substantially mitigated the disseminate of SARS-CoV-2. Importantly, model simulations predict that 95% reduction of outbreak on June 26, 2020 and 99% reduction of outbreak on July 26, 2020 in India. Interpretation: Our model simulation demonstrates that the community-wide elimination of SARS-CoV-2 is possible by mitigating the social distancing and use essential precautions. Lockdown can be implemented strictly to prevent the human-to-human transmission of COVID-19. The model-based and parameter estimation of epidemic life cycle, and end dates, if can be done precisely, may decrease distress and over optimism and develop the mentality for all of us for the next stages of the outbreak evolution.


title: Brane-jet stability of non-supersymmetric AdS vacua

abstract: We classify the non-supersymmetric, and perturbatively stable within $D=4$, AdS vacua of maximal $D=4$ supergravity with a dyonic ISO(7) gauging in a large sector of the supergravity. Seven such vacua are established within this sector, all of them giving rise to non-supersymmetric $\textrm{AdS}_{4} \times \textrm{S}^{6}$ type IIA backgrounds with and without non-trivial warpings and with internal fluxes. Then, we analyse the dynamics of various probe D$p$-branes in these backgrounds searching for potential brane-jet instabilities. In all these cases, such instabilities are absent. Finally, an alternative decay channel through tunnelling is investigated, focusing on one of the seven backgrounds. We do not find instabilities either, but the analysis remains inconclusive.


title: Probabilistic Guarantees for Safe Deep Reinforcement Learning

abstract: Deep reinforcement learning has been successfully applied to many control tasks, but the application of such agents in safety-critical scenarios has been limited due to safety concerns. Rigorous testing of these controllers is challenging, particularly when they operate in probabilistic environments due to, for example, hardware faults or noisy sensors. We propose MOSAIC, an algorithm for measuring the safety of deep reinforcement learning agents in stochastic settings. Our approach is based on the iterative construction of a formal abstraction of a controller's execution in an environment, and leverages probabilistic model checking of Markov decision processes to produce probabilistic guarantees on safe behaviour over a finite time horizon. It produces bounds on the probability of safe operation of the controller for different initial configurations and identifies regions where correct behaviour can be guaranteed. We implement and evaluate our approach on agents trained for several benchmark control problems.


title: FaceFilter: Audio-visual speech separation using still images

abstract: The objective of this paper is to separate a target speaker's speech from a mixture of two speakers using a deep audio-visual speech separation network. Unlike previous works that used lip movement on video clips or pre-enrolled speaker information as an auxiliary conditional feature, we use a single face image of the target speaker. In this task, the conditional feature is obtained from facial appearance in cross-modal biometric task, where audio and visual identity representations are shared in latent space. Learnt identities from facial images enforce the network to isolate matched speakers and extract the voices from mixed speech. It solves the permutation problem caused by swapped channel outputs, frequently occurred in speech separation tasks. The proposed method is far more practical than video-based speech separation since user profile images are readily available on many platforms. Also, unlike speaker-aware separation methods, it is applicable on separation with unseen speakers who have never been enrolled before. We show strong qualitative and quantitative results on challenging real-world examples.


title: You Only Search Once: A Fast Automation Framework for Single-Stage DNN/Accelerator Co-design

abstract: DNN/Accelerator co-design has shown great potential in improving QoR and performance. Typical approaches separate the design flow into two-stage: (1) designing an application-specific DNN model with high accuracy; (2) building an accelerator considering the DNN specific characteristics. However, it may fail in promising the highest composite score which combines the goals of accuracy and other hardware-related constraints (e.g., latency, energy efficiency) when building a specific neural-network-based system. In this work, we present a single-stage automated framework, YOSO, aiming to generate the optimal solution of software-and-hardware that flexibly balances between the goal of accuracy, power, and QoS. Compared with the two-stage method on the baseline systolic array accelerator and Cifar10 dataset, we achieve 1.42x~2.29x energy or 1.79x~3.07x latency reduction at the same level of precision, for different user-specified energy and latency optimization constraints, respectively.


title: Evolutionary Simplicial Learning as a Generative and Compact Sparse Framework for Classification

abstract: Dictionary learning for sparse representations has been successful in many reconstruction tasks. Simplicial learning is an adaptation of dictionary learning, where subspaces become clipped and acquire arbitrary offsets, taking the form of simplices. Such adaptation is achieved through additional constraints on sparse codes. Furthermore, an evolutionary approach can be chosen to determine the number and the dimensionality of simplices composing the simplicial, in which most generative and compact simplicials are favored. This paper proposes an evolutionary simplicial learning method as a generative and compact sparse framework for classification. The proposed approach is first applied on a one-class classification task and it appears as the most reliable method within the considered benchmark. Most surprising results are observed when evolutionary simplicial learning is considered within a multi-class classification task. Since sparse representations are generative in nature, they bear a fundamental problem of not being capable of distinguishing two classes lying on the same subspace. This claim is validated through synthetic experiments and superiority of simplicial learning even as a generative-only approach is demonstrated. Simplicial learning loses its superiority over discriminative methods in high-dimensional cases but can further be modified with discriminative elements to achieve state-of-the-art performance in classification tasks.


title: Stress-energy of the quantized fields in the spacetime of the Damour-Solodukhin wormhole

abstract: The static traversable wormhole should made out of some type of exotic matter which satisfies the Morris-Thorne conditions. Although the characteristic size of the region with the exotic matter can be made arbitrary small, the calculations performed so far suggest that the Morris-Thorne conditions are quite restrictive and it is hard to find matter with the desired properties. Traditionally, the quantized fields are considered as the best candidates because they can violate the weak-energy condition. In this paper we employ the Schwinger-DeWitt expansion to construct and examine the approximate stress-energy tensor of the quantized massive scalar (with an arbitrary curvature coupling), spinor and vector field in the spacetime of the Damour-Solodukhin wormhole. We find that for the scalar field there is a region in a parameter space in which the stress-energy tensor has the desired properties. That means that of the twenty-one cases considered so far (the seven types of the wormhole geometries and the three types of the massive fields) only in the four cases (for certain values of the parameters) the stress-energy tensor does satisfy the Morris-Thorne conditions.


title: Anomaly Detection And Classification In Time Series With Kervolutional Neural Networks

abstract: Recently, with the development of deep learning, end-to-end neural network architectures have been increasingly applied to condition monitoring signals. They have demonstrated superior performance for fault detection and classification, in particular using convolutional neural networks. Even more recently, an extension of the concept of convolution to the concept of kervolution has been proposed with some promising results in image classification tasks. In this paper, we explore the potential of kervolutional neural networks applied to time series data. We demonstrate that using a mixture of convolutional and kervolutional layers improves the model performance. The mixed model is first applied to a classification task in time series, as a benchmark dataset. Subsequently, the proposed mixed architecture is used to detect anomalies in time series data recorded by accelerometers on helicopters. We propose a residual-based anomaly detection approach using a temporal auto-encoder. We demonstrate that mixing kervolutional with convolutional layers in the encoder is more sensitive to variations in the input data and is able to detect anomalous time series in a better way.


title: Should Classical Physics Be Interpreted Indeterministically?

abstract: An indeterministic interpretation of classical physics has been proposed recently, in which the argument relies on attacking an alleged unwarranted metaphysical hidden assumption of the standard deterministic interpretation. This short paper aims at showing that it is arguably a strawman attack.


title: Optimal long-term investment in illiquid markets when prices have negative memory

abstract: In a discrete-time financial market model with instantaneous price impact, we find an asymptotically optimal strategy for an investor maximizing her expected wealth. The asset price is assumed to follow a process with negative memory. We determine how the optimal growth rate depends on the impact parameter and on the covariance decay rate of the price.


title: A Berkeley View of Teaching CS at Scale

abstract: Over the past decade, undergraduate Computer Science (CS) programs across the nation have experienced an explosive growth in enrollment as computational skills have proven increasingly important across many domains and in the workforce at large. Motivated by this unprecedented student demand, the CS program at the University of California, Berkeley has tripled the size of its graduating class in five years. The first two introductory courses for majors, each taught by one faculty instructor and several hundred student teachers, combine to serve nearly 2,900 students per term. This report presents three strategies that have enabled the effective teaching, delivery, and management of large-scale CS courses: (1) the development of autograder infrastructure and online platforms to provide instant feedback with minimal instructor intervention and deliver the course at scale; (2) the expansion of academic and social student support networks resulting from changes in teaching assistant responsibilities and the development of several near-peer mentoring communities; and (3) the expansion of undergraduate teacher preparation programs to meet the increased demand for qualified student teachers. These interventions have helped both introductory and advanced courses address capacity challenges and expand enrollments while receiving among the highest student evaluations of teaching in department history. Implications for inclusivity and diversity are discussed.


title: Theoretical evidence of a significant modification of the electronic structure of double-walled carbon nanotubes due to the interlayer interaction

abstract: The recently reported experimental optical spectra of double-walled carbon nanotubes exhibit more peaks than it could be expected based on the layers alone. The appearance of excess peaks has been attributed to the interlayer interaction. In order to elucidate the origin of the excess peaks, we perform calculations of the optical absorption of a particular nanotube using the recursion method with non-orthogonal tight-binding basis functions. Our study shows that the interlayer interaction can give rise to major changes in the electronic structure of this nanotube, manifesting themselves with shifts of the optical transitions and appearance of new optical transitions. The derived absorption spectrum is found to be in excellent agreement with the available experimental data, which justifies the use of the proposed approach for double-walled carbon nanotubes.


title: Connectivity estimation of high dimensional data recorded from neuronal cells

abstract: The main result of this thesis is the development of a novel connectivity estimation method, called Total Spiking Probability Edges (TSPE). Based on cross-correlation and edge filtering at different time scales this method is proposed and the theoretical framework is outlined in this work. TSPE enables the classification between inhibitory and excitatory connections by using recorded action potentials. To compare this method learning about state of the art algorithms to estimate connectivity is necessary. After a research, promising algorithms are implemented and evaluated for further research topics, among others in the biomems lab of UAS Aschaffenburg. To evaluate these algorithms in silico networks are used, because of their known connectivity. This makes it possible to validate the correctness of our algorithm results. Therefore, a biophysically representative neuronal network simulation is needed first. Datasets were simulated in different ways and analysed in order to develop an evaluation framework. After a successful evaluation with in silico networks, in vitro experiments and their analyses complete this project.


title: Consensus-Based Global Optimization with Personal Best

abstract: In this paper we propose a variant of a consensus-based global optimization (CBO) method that uses personal best information in order to compute the global minimum of a non-convex, locally Lipschitz continuous function. The proposed approach is motivated by the original particle swarming algorithms, in which particles adjust their position with respect to the personal best, the current global best, and some additive noise. The personal best information along an individual trajectory is included with the help of a weighted mean. This weighted mean can be computed very efficiently due to its accumulative structure. It enters the dynamics via an additional drift term. We illustrate the performance with a toy example, analyze the respective memory-dependent stochastic system and compare the performance with the original CBO with component-wise noise for several benchmark problems. The proposed method has a higher success rate for computational experiments with a small particle number and where the initial particle distribution is disadvantageous with respect to the global minimum.


title: Towards Lagrangian construction for infinite half-integer spin field

abstract: We formulate the conditions for the generalized fields in the space with additional commuting Weyl spinor coordinates which define the infinite half-integer spin representation of the four-dimensional Poincar\'e group. Using this formulation we develop the BRST approach and derive the Lagrangian for the half-integer infinite spin fields.


title: A possible transit of a disintegrating exoplanet in the nearby multiplanet system DMPP-1

abstract: We analyse TESS photometry of DMPP-1 (HD 38677; TIC 66560666), a nearby F8V star hosting hot super-Earth planets and a warm Neptune. Using the Transit Least Squares algorithm and other methods we find a transit signal at $P=3.2854^{+0.0032}_{-0.0025}$ d with depth 87$^{+25}_{-30}$ ppm and false alarm probability 1.6%. This is shallower than hitherto published TESS discoveries. The 3.285 d signal is recovered for several, but not all, methods for detrending stellar astrophysical variability. Further observations are needed to improve the significance of the detection. If this transit were due to an Earth-like rocky planet it would have been detected in the RV data, but it is not. The TESS data cover seven individual transits, one of which is consistent with zero depth. The insolation of the putative planet is $990 {\rm S_{\oplus}}$, typical of fluxes experienced by the three known catastrophically disintegrating exoplanets (CDEs). The transits can be self-consistently attributed to a CDE with a mass below the RV detection threshold. We searched for transits of the known RV planets, finding null results and detection thresholds of $< 100$ ppm, which we quantify for each. The DMPP-1 planetary system was discovered as a consequence of circumstellar gas attributed to ablation of hot planets. The RV planets may have been ablated to near-pure iron cores. We place limits on the orbital inclinations of the RV planets where the expected transit depth exceeds the detection threshold. If the 3.2854 day transit detection is confirmed, e.g. with CHEOPS photometry, DMPP-1 would be a first-rate target for JWST spectroscopy.


title: A new technique for preserving conservation laws

abstract: This paper introduces a new symbolic-numeric strategy for finding semidiscretizations of a given PDE that preserve multiple local conservation laws. We prove that for one spatial dimension, various one-step time integrators from the literature preserve fully discrete local conservation laws whose densities are either quadratic or a Hamiltonian. The approach generalizes to time integrators with more steps and conservation laws of other kinds; higher-dimensional PDEs can be treated by iterating the new strategy. We use the Boussinesq equation as a benchmark and introduce new families of schemes of order two and four that preserve three conservation laws. We show that the new technique is practicable for PDEs with three dependent variables, introducing as an example new families of second-order schemes for the potential Kadomtsev-Petviashvili equation.


title: Structural origin of plasticity in strained high-entropy alloy

abstract: High-entropy alloys (HEAs) are solid solutions of multiple elements with equal atomic ratios which present an innovative pathway for de novo alloy engineering. While there exist extensive studies to ascertain the important structural aspects governing their mechanical behaviors, elucidating the underlying deformation mechanisms still remains a challenge. Using atomistic simulations, we probe the particle rearrangements in a yielding, model HEA system to understand the structural origin of its plasticity. We find the plastic deformation is initiated by irreversible topological fluctuations which tend to spatially localize in regions termed as soft spots which consist of particles actively participating in slow vibrational motions, an observation strikingly reminiscent of nonlinear glassy rheology. Due to the varying local elastic moduli resulting from the loss of compositional periodicity, these plastic responses exhibit significant spatial heterogeneity and are found to be inversely correlated with the distribution of local electronegativity. Further mechanical loading promotes the cooperativity among these local plastic events and triggers the formation of dislocation loops. As in strained crystalline solids, different dislocation loops can further merge together and propagate as the main carrier of large-scale plastic deformation. However, the energy barriers located at the spatial regions with higher local electronegativity severely hinders the motion of dislocations. By delineating the transient mechanical response in terms of atomic configuration, our computational findings shed new light on understanding the nature of plasticity of single-phase HEA.


title: ECoPANN: A Framework for Estimating Cosmological Parameters using Artificial Neural Networks

abstract: In this work, we present a new method to estimate cosmological parameters accurately based on the artificial neural network (ANN), and a code called ECoPANN (Estimating Cosmological Parameters with ANN) is developed to achieve parameter inference. We test the ANN method by estimating the basic parameters of the concordance cosmological model using the simulated temperature power spectrum of the cosmic microwave background (CMB). The results show that the ANN performs excellently on best-fit values and errors of parameters, as well as correlations between parameters when compared with that of the Markov Chain Monte Carlo (MCMC) method. Besides, for a well-trained ANN model, it is capable of estimating parameters for multiple experiments that have different precisions, which can greatly reduce the consumption of time and computing resources for parameter inference. Furthermore, we extend the ANN to a multibranch network to achieve a joint constraint on parameters. We test the multibranch network using the simulated temperature and polarization power spectra of the CMB, Type Ia supernovae, and baryon acoustic oscillations, and almost obtain the same results as the MCMC method. Therefore, we propose that the ANN can provide an alternative way to accurately and quickly estimate cosmological parameters, and ECoPANN can be applied to the research of cosmology and even other broader scientific fields.


title: SND electromagnetic calorimeter time measurement and its applications

abstract: The SND is a non-magnetic detector deployed at the VEPP-2000 $e^+e^-$ collider (BINP, Novosibirsk) for hadronic cross-section measurements in the center of mass energy region below 2 GeV. The important part of the detector is a three-layer hodoscopic electromagnetic calorimeter (EMC) based on NaI(Tl) counters. Until the recent EMC spectrometric channel upgrade, only the energy deposition measurement in counters was possible. A new EMC signal shaping and digitizing electronics based on FADC allows us to obtain also the event time structure. The new electronics and supporting software, including digital signal processing algorithms, are used for data taking in the ongoing experiment. We discuss the amplitude and time extraction algorithms, the new system performance on experimental events and physical analysis applications.


title: Semi-supervised Neural Chord Estimation Based on a Variational Autoencoder with Latent Chord Labels and Features

abstract: This paper describes a statistically-principled semi-supervised method of automatic chord estimation (ACE) that can make effective use of music signals regardless of the availability of chord annotations. The typical approach to ACE is to train a deep classification model (neural chord estimator) in a supervised manner by using only annotated music signals. In this discriminative approach, prior knowledge about chord label sequences (model output) has scarcely been taken into account. In contrast, we propose a unified generative and discriminative approach in the framework of amortized variational inference. More specifically, we formulate a deep generative model that represents the generative process of chroma vectors (observed variables) from discrete labels and continuous features (latent variables), which are assumed to follow a Markov model favoring self-transitions and a standard Gaussian distribution, respectively. Given chroma vectors as observed data, the posterior distributions of the latent labels and features are computed approximately by using deep classification and recognition models, respectively. These three models form a variational autoencoder and can be trained jointly in a semi-supervised manner. The experimental results show that the regularization of the classification model based on the Markov prior of chord labels and the generative model of chroma vectors improved the performance of ACE even under the supervised condition. The semi-supervised learning using additional non-annotated data can further improve the performance.


title: A Many-Body Theory of the Optical Conductivity of Excitons and Trions in Two-Dimensional Materials

abstract: The optical spectra of two dimensional (2D) materials exhibit sharp absorption peaks that are commonly identified with exciton and trions (or charged excitons). In this paper, we show that excitons and trions in doped 2D materials can be described by two coupled Schrodinger-like equations - one two-body equation for excitons and another four-body equation for trions. In electron doped 2D materials, a bound trion state is identified with a four-body bound state of an exciton and an excited conduction band electron-hole pair. In doped 2D materials, the exciton and the trions states are the not the eigenstates of the full Hamiltonian and their respective Schrodinger equations are coupled due to Coulomb interactions. The strength of this coupling increases with the doping density. Solutions of these two coupled equations can quantitatively explain all the prominent features experimentally observed in the optical absorption spectra of 2D materials including the observation of two prominent absorption peaks and the variation of their energy splittings and spectral shapes and strengths with the electron density. The optical conductivity obtained in our work satisfies the optical conductivity sum rule exactly. A superposition of exciton and trion states can be used to construct a solution of the two coupled Schrodinger equations and this solution resembles the variational exciton-polaron state, thereby establishing the relationship between our approach and Fermi polaron physics.


title: Bayesian Bits: Unifying Quantization and Pruning

abstract: We introduce Bayesian Bits, a practical method for joint mixed precision quantization and pruning through gradient based optimization. Bayesian Bits employs a novel decomposition of the quantization operation, which sequentially considers doubling the bit width. At each new bit width, the residual error between the full precision value and the previously rounded value is quantized. We then decide whether or not to add this quantized residual error for a higher effective bit width and lower quantization noise. By starting with a power-of-two bit width, this decomposition will always produce hardware-friendly configurations, and through an additional 0-bit option, serves as a unified view of pruning and quantization. Bayesian Bits then introduces learnable stochastic gates, which collectively control the bit width of the given tensor. As a result, we can obtain low bit solutions by performing approximate inference over the gates, with prior distributions that encourage most of them to be switched off. We experimentally validate our proposed method on several benchmark datasets and show that we can learn pruned, mixed precision networks that provide a better trade-off between accuracy and efficiency than their static bit width equivalents.


title: Approval-Based Shortlisting

abstract: Shortlisting is the task of reducing a long list of alternatives to a (smaller) set of best or most suitable alternatives. Shortlisting is often used in the nomination process of awards or in recommender systems to display featured objects. In this paper, we analyze shortlisting methods that are based on approval data, a common type of preferences. Furthermore, we assume that the size of the shortlist, i.e., the number of best or most suitable alternatives, is not fixed but determined by the shortlisting method. We axiomatically analyze established and new shortlisting methods and complement this analysis with an experimental evaluation based on synthetic and real-world data. Our results lead to recommendations which shortlisting methods to use, depending on the desired properties.


title: Pressure corrections in decoupling SU(2) Yang-Mills Theory: The case of dihedral diagrams involving both massive and massless modes

abstract: In this work we show the step by step calculations needed to quantify the contribution of a three-loop order diagram with dihedral symmetry to the radiative corrections of the pressure in SU(2) thermal Yang-Mills theory in deconfining phase. We surveyed past developments, and performed computations for separate channel combinations, defined by Mandelstam variables which are constrained by two 4-vertices. An analytically integrable approximation for high-temperature conditions was found, to verify the relevance of the corrections for this diagram. A numerical analysis with Monte Carlo methods was carried out to check the validity of such approximation, to compare it with the full integral. A Dyson-Schwinger resummation had to be performed to all dihedral loop orders in order to control the temperature dependency found.


title: WISE J135501.90-825838.9 is a Nearby, Young, Extremely Low-mass Substellar Binary

abstract: We present a parallax solution for WISE J135501.90-825838.9, a spectral binary with spectral types L7+T7.5 and candidate AB Doradus member. Using $WISE$ astrometry, we obtain a distance of $d = 16.7\pm5.3$ pc. This preliminary parallax solution provides further evidence that WISE J135501.90-825838.9 is a member of AB Doradus (130-200 Myr), and when combined with evolutionary models predicts masses of 11 $M_\mathrm{Jup}$ and 9 $M_\mathrm{Jup}$ for both components.


title: Ambient Sound Helps: Audiovisual Crowd Counting in Extreme Conditions

abstract: Visual crowd counting has been recently studied as a way to enable people counting in crowd scenes from images. Albeit successful, vision-based crowd counting approaches could fail to capture informative features in extreme conditions, e.g., imaging at night and occlusion. In this work, we introduce a novel task of audiovisual crowd counting, in which visual and auditory information are integrated for counting purposes. We collect a large-scale benchmark, named auDiovISual Crowd cOunting (DISCO) dataset, consisting of 1,935 images and the corresponding audio clips, and 170,270 annotated instances. In order to fuse the two modalities, we make use of a linear feature-wise fusion module that carries out an affine transformation on visual and auditory features. Finally, we conduct extensive experiments using the proposed dataset and approach. Experimental results show that introducing auditory information can benefit crowd counting under different illumination, noise, and occlusion conditions. The dataset and code will be released. Code and data have been made available


title: Seiberg-Witten-Casson Invariant of Homology $S^1 \times S^3$ with Circle Action

abstract: In this paper we shall compute the Mrowka-Ruberman-Saveliev invariant introduced in [17] for the case when the manifold admits a free circle action.


title: Stealthy and Efficient Adversarial Attacks against Deep Reinforcement Learning

abstract: Adversarial attacks against conventional Deep Learning (DL) systems and algorithms have been widely studied, and various defenses were proposed. However, the possibility and feasibility of such attacks against Deep Reinforcement Learning (DRL) are less explored. As DRL has achieved great success in various complex tasks, designing effective adversarial attacks is an indispensable prerequisite towards building robust DRL algorithms. In this paper, we introduce two novel adversarial attack techniques to \emph{stealthily} and \emph{efficiently} attack the DRL agents. These two techniques enable an adversary to inject adversarial samples in a minimal set of critical moments while causing the most severe damage to the agent. The first technique is the \emph{critical point attack}: the adversary builds a model to predict the future environmental states and agent's actions, assesses the damage of each possible attack strategy, and selects the optimal one. The second technique is the \emph{antagonist attack}: the adversary automatically learns a domain-agnostic model to discover the critical moments of attacking the agent in an episode. Experimental results demonstrate the effectiveness of our techniques. Specifically, to successfully attack the DRL agent, our critical point technique only requires 1 (TORCS) or 2 (Atari Pong and Breakout) steps, and the antagonist technique needs fewer than 5 steps (4 Mujoco tasks), which are significant improvements over state-of-the-art methods.


title: Spectrum of the Hidden-Bottom and the Hidden-Charm/Strange Exotics in the Dynamical Diquark Model

abstract: The lightest hidden-bottom tetraquarks in the dynamical diquark model fill an $S$-wave multiplet consisting of 12 isomultiplets. We predict their masses and dominant bottomonium decay channels using a simple 3-parameter Hamiltonian that captures the core fine-structure features of the model, including isospin dependence. The only experimental inputs needed are the corresponding observables for $Z_b(10610)$ and $Z_b(10650)$. The mass of $X_b$, the bottom analogue to $X(3872)$, is highly constrained in this scheme. In addition, using lattice-calculated potentials we predict the location of the center of mass of the $P$-wave multiplet and find that $Y(10860)$ fits well but the newly discovered $Y(10750)$ does not, more plausibly being a $D$-wave bottomonium state. Using similar methods, we also examine the lowest $S$-wave multiplet of 6 $c\bar c s\bar s$ states, assuming as in earlier work that $X(3915)$ and $Y(4140)$ are members, and predict the masses and dominant charmonium decay modes of the other states. We again use lattice potentials to compute the centers of mass of higher multiplets, and find them to be compatible with the masses of $Y(4626)$ ($1P$) and $X(4700)$ ($2S$), respectively.


title: Scalar Active Mixtures: The Non-Reciprocal Cahn-Hilliard Model

abstract: Pair interactions between active particles need not follow Newton's third law. In this work we propose a continuum model of pattern formation due to non-reciprocal interaction between multiple species of scalar active matter. The classical Cahn-Hilliard model is minimally modified by supplementing the equilibrium Ginzburg-Landau dynamics with particle number conserving currents which cannot be derived from a free energy, reflecting the microscopic departure from action-reaction symmetry. The strength of the asymmetry in the interaction determines whether the steady state exhibits a macroscopic phase separation or a traveling density wave displaying global polar order. The latter structure, which is equivalent to an active self-propelled smectic phase, coarsens via annihilation of defects, whereas the former structure undergoes Ostwald ripening. The emergence of traveling density waves, which is a clear signature of broken time-reversal symmetry in this active system, is a generic feature of any multi-component mixture with microscopic non-reciprocal interactions.

categories: cond-mat.stat-mech cond-mat.soft physics.bio-ph physics.chem-ph


title: Pair production of dark particles in meson decays

abstract: Rare decays of $K$ and $B$ mesons provide a powerful probe of dark sectors with light new particles. We show that the pair production of $O(100\,{\rm MeV})$ dark states can be probed with the decays of $K_L$ mesons, owing to the enhanced two-body kinematics, $K_L\to X_1X_2$ or $X_2X_2$. If either or these two particles is unstable, e.g. $X_2\to X_1\pi^0$, $X_2\to X_1\gamma$ or $X_{1,2}\to \gamma\gamma$, such decays could easily mimic $K_L\to \pi^0 \nu\overline{\nu}$ signatures, while not being ruled out by the decays of charged kaons. We construct explicit models that have enhanced $K_L$ decay signatures, and are constrained by the results of the KOTO experiment. We note that recently reported excess events can also be accommodated while satisfying all other constraints ($B$ decays, colliders, beam dumps). These models are based on the extensions of the gauge and/or scalar sector of the theory. The lightest of $X_{1,2}$ particles, if stable, could constitute the entirety of dark matter.


title: Phase transition in cohomology groups of non-uniform random simplicial complexes

abstract: We consider a generalised model of a random simplicial complex, which arises from a random hypergraph. Our model is generated by taking the downward-closure of a non-uniform binomial random hypergraph, in which for each $k$, each set of $k+1$ vertices forms an edge with some probability $p_k$ independently. As a special case, this contains an extensively studied model of a (uniform) random simplicial complex, introduced by Meshulam and Wallach [Random Structures & Algorithms 34 (2009), no. 3, pp. 408-417]. We consider a higher-dimensional notion of connectedness on this new model according to the vanishing of cohomology groups over an arbitrary abelian group $R$. We prove that this notion of connectedness displays a phase transition and determine the threshold. We also prove a hitting time result for a natural process interpretation, in which simplices and their downward-closure are added one by one. In addition, we determine the asymptotic behaviour of cohomology groups inside the critical window around the time of the phase transition.


title: Some computations on the characteristic variety of a line arrangement

abstract: We find monodromy formulas for line arrangements which are fibered with respect to the projection from one point. We use them to find $0$-dimensional translated components in the first characteristic variety of the arrangement $\mathcal R(2n)$ determined by a regular $n$-polygon and its diagonals.


title: ZeroShotCeres: Zero-Shot Relation Extraction from Semi-Structured Webpages

abstract: In many documents, such as semi-structured webpages, textual semantics are augmented with additional information conveyed using visual elements including layout, font size, and color. Prior work on information extraction from semi-structured websites has required learning an extraction model specific to a given template via either manually labeled or distantly supervised data from that template. In this work, we propose a solution for "zero-shot" open-domain relation extraction from webpages with a previously unseen template, including from websites with little overlap with existing sources of knowledge for distant supervision and websites in entirely new subject verticals. Our model uses a graph neural network-based approach to build a rich representation of text fields on a webpage and the relationships between them, enabling generalization to new templates. Experiments show this approach provides a 31% F1 gain over a baseline for zero-shot extraction in a new subject vertical.


title: Electrically Tunable Room Temperature Hysteresis Crossover in Underlap MoS$_2$ FETs

abstract: Clockwise to anti-clockwise hysteresis crossover in current-voltage transfer characteristics of field effect transistors (FETs) with graphene and MoS$_2$ channels holds significant promise for non-volatile memory applications. However such crossovers have been shown to manifest only at high temperature. In this work, for the first time, we demonstrate room temperature hysteresis crossover in few-layer MoS$_2$ FETs by using a gate-drain underlap design to induce a differential response from traps at the MoS$_2$-HfO$_2$ channel-gate dielectric interface to applied gate bias. The appearance of interface trap-driven anti-clockwise hysteresis at high gate voltages in underlap FETs can be unambiguously attributed to the presence of an underlap since transistors with and without the underlap region were fabricated on the same MoS$_2$ channel flake. The underlap design also enables room temperature tuning of the anti-clockwise hysteresis window (by 140$\times$) as well as the crossover gate voltage (by 2.6$\times$) with applied drain bias and underlap length. Comprehensive measurements of the transfer curves in ambient and vacuum conditions at varying sweep rates and temperatures (RT, 45 $^\circ$C and 65 $^\circ$C) help segregate the quantitative contributions of adsorbates, interface traps, and bulk HfO$_2$ traps to the clockwise and anti-clockwise hysteresis.


title: Natural Way to Overcome the Catastrophic Forgetting in Neural Networks

abstract: Not so long ago, a method was discovered that successfully overcomes the catastrophic forgetting in neural networks. Although we know about the cases of using this method to preserve skills when adapting pre-trained networks to particular tasks, it has not obtained widespread distribution yet. In this paper, we would like to propose an alternative method of overcoming catastrophic forgetting based on the total absolute signal passing through each connection in the network. This method has a simple implementation and seems to us essentially close to the processes occurring in the brain of animals to preserve previously learned skills during subsequent learning. We hope that the ease of implementation of this method will serve its wide application.


title: Effective Model of Loop Extrusion Predicts Chromosomal Domains

abstract: An active loop-extrusion mechanism is regarded as the main out--of--equilibrium mechanism responsible for the structuring of megabase-sized domains in chromosomes. We developed a model to study the dynamics of the chromosome fibre by solving the kinetic equations associated with the motion of the extruder. By averaging out the position of the extruder along the chain, we build an effective equilibrium model capable of reproducing experimental contact maps based solely on the positions of extrusion--blocking proteins. We assessed the quality of the effective model using numerical simulations of chromosomal segments and comparing the results with explicit-extruder models and experimental data.


title: The $\alpha$-$\eta$-$\mathcal{F}$ and $\alpha$-$\kappa$-$\mathcal{F}$ Composite Fading Distributions

abstract: In this paper, we present the $\alpha$-$\eta$-$\mathcal{F}$ and $\alpha$-$\kappa$-$\mathcal{F}$ composite fading distributions. The two distributions generalize the two well-known composite fading distributions, namely the $\eta$-$\mu$/inverse gamma and the $\kappa$-$\mu$/inverse gamma distributions. For both distributions, we derive new exact closed-form expressions for the probability density function (PDF) and cumulative density function (CDF) of the instantaneous signal-to-noise ratio (SNR). Additionally, the outage probabilities (OPs) of a wireless communication system operating over the new composite fading channels are obtained. Moreover, we derive asymptotic OPs, in the high SNR regime, to gain more insight into the influence of fading parameters. Some numerical results on the envelope PDFs and OPs are provided and compared with Monte-Carlo simulation results to validate the analyses.


title: Robust On-Manifold Optimization for Uncooperative Space Relative Navigation with a Single Camera

abstract: Optical cameras are gaining popularity as the suitable sensor for relative navigation in space due to their attractive sizing, power and cost properties when compared to conventional flight hardware or costly laser-based systems. However, a camera cannot infer depth information on its own, which is often solved by introducing complementary sensors or a second camera. In this paper, an innovative model-based approach is instead demonstrated to estimate the six-dimensional pose of a target object relative to the chaser spacecraft using solely a monocular setup. The observed facet of the target is tackled as a classification problem, where the three-dimensional shape is learned offline using Gaussian mixture modeling. The estimate is refined by minimizing two different robust loss functions based on local feature correspondences. The resulting pseudo-measurements are then processed and fused with an extended Kalman filter. The entire optimization framework is designed to operate directly on the $SE\text{(3)}$ manifold, uncoupling the process and measurement models from the global attitude state representation. It is validated on realistic synthetic and laboratory datasets of a rendezvous trajectory with the complex spacecraft Envisat. It is demonstrated how it achieves an estimate of the relative pose with high accuracy over its full tumbling motion.


title: Distilling neural networks into skipgram-level decision lists

abstract: Several previous studies on explanation for recurrent neural networks focus on approaches that find the most important input segments for a network as its explanations. In that case, the manner in which these input segments combine with each other to form an explanatory pattern remains unknown. To overcome this, some previous work tries to find patterns (called rules) in the data that explain neural outputs. However, their explanations are often insensitive to model parameters, which limits the scalability of text explanations. To overcome these limitations, we propose a pipeline to explain RNNs by means of decision lists (also called rules) over skipgrams. For evaluation of explanations, we create a synthetic sepsis-identification dataset, as well as apply our technique on additional clinical and sentiment analysis datasets. We find that our technique persistently achieves high explanation fidelity and qualitatively interpretable rules.


title: See Change: VLT spectroscopy of a sample of high-redshift Type Ia supernova host galaxies

abstract: The Supernova Cosmology Project has conducted the `See Change' programme, aimed at discovering and observing high-redshift (1.13 $\leq$ z $\leq$ 1.75) Type Ia supernovae (SNe Ia). We used multi-filter Hubble Space Telescope (HST) observations of massive galaxy clusters with sufficient cadence to make the observed SN Ia light curves suitable for a cosmological probe of dark energy at z > 0.5. This See Change sample of SNe Ia with multi-colour light curves will be the largest to date at these redshifts. As part of the See Change programme, we obtained ground-based spectroscopy of each discovered transient and/or its host galaxy. Here we present Very Large Telescope (VLT) spectra of See Change transient host galaxies, deriving their redshifts, and host parameters such as stellar mass and star formation rate. Of the 39 See Change transients/hosts that were observed with the VLT, we successfully determined the redshift for 26, including 15 SNe Ia at z > 0.97. We show that even in passive environments, it is possible to recover secure redshifts for the majority of SN hosts out to z = 1.5. We find that with typical exposure times of 3 - 4 hrs on an 8m-class telescope we can recover ~75% of SN Ia redshifts in the range of 0.97 < z < 1.5. Furthermore, we show that the combination of HST photometry and VLT spectroscopy is able to provide estimates of host galaxy stellar mass that are sufficiently accurate for use in a mass-step correction in the cosmological analysis.


title: Staggered Diagonal Embedding Based Linear Field Size Streaming Codes

abstract: An $(a,b,\tau)$ streaming code is a packet-level erasure code that can recover under a strict delay constraint of $\tau$ time units, from either a burst of $b$ erasures or else of $a$ random erasures, occurring within a sliding window of time duration $w$. While rate-optimal constructions of such streaming codes are available for all parameters $\{a,b,\tau,w\}$ in the literature, they require in most instances, a quadratic, $O(\tau^2)$ field size. In this work, we make further progress towards field size reduction and present rate-optimal $O(\tau)$ field size streaming codes for two regimes: (i) $gcd(b,\tau+1-a)\ge a$ (ii) $\tau+1 \ge a+b$ and $b \mod \ a \in \{0,a-1\}$.


title: A Deeper Look at the Unsupervised Learning of Disentangled Representations in $\beta$-VAE from the Perspective of Core Object Recognition

abstract: The ability to recognize objects despite there being differences in appearance, known as Core Object Recognition, forms a critical part of human perception. While it is understood that the brain accomplishes Core Object Recognition through feedforward, hierarchical computations through the visual stream, the underlying algorithms that allow for invariant representations to form downstream is still not well understood. (DiCarlo et al., 2012) Various computational perceptual models have been built to attempt and tackle the object identification task in an artificial perceptual setting. Artificial Neural Networks, computational graphs consisting of weighted edges and mathematical operations at vertices, are loosely inspired by neural networks in the brain and have proven effective at various visual perceptual tasks, including object characterization and identification. (Pinto et al., 2008) (DiCarlo et al., 2012) For many data analysis tasks, learning representations where each dimension is statistically independent and thus disentangled from the others is useful. If the underlying generative factors of the data are also statistically independent, Bayesian inference of latent variables can form disentangled representations. This thesis constitutes a research project exploring a generalization of the Variational Autoencoder (VAE), $\beta$-VAE, that aims to learn disentangled representations using variational inference. $\beta$-VAE incorporates the hyperparameter $\beta$, and enforces conditional independence of its bottleneck neurons, which is in general not compatible with the statistical independence of latent variables. This text examines this architecture, and provides analytical and numerical arguments, with the goal of demonstrating that this incompatibility leads to a non-monotonic inference performance in $\beta$-VAE with a finite optimal $\beta$.


title: CoSimGNN: Towards Large-scale Graph Similarity Computation

abstract: The ability to compute similarity scores between graphs based on metrics such as Graph Edit Distance (GED) is important in many real-world applications. Computing exact GED values is typically an NP-hard problem and traditional algorithms usually achieve an unsatisfactory trade-off between accuracy and efficiency. Recently, Graph Neural Networks (GNNs) provide a data-driven solution for this task, which is more efficient while maintaining prediction accuracy in small graph (around 10 nodes per graph) similarity computation. Existing GNN-based methods, which either respectively embeds two graphs (lack of low-level cross-graph interactions) or deploy cross-graph interactions for whole graph pairs (redundant and time-consuming), are still not able to achieve competitive results when the number of nodes in graphs increases. In this paper, we focus on similarity computation for large-scale graphs and propose the "embedding-coarsening-matching" framework CoSimGNN, which first embeds and coarsens large graphs with adaptive pooling operation and then deploys fine-grained interactions on the coarsened graphs for final similarity scores. Furthermore, we create several synthetic datasets which provide new benchmarks for graph similarity computation. Detailed experiments on both synthetic and real-world datasets have been conducted and CoSimGNN achieves the best performance while the inference time is at most 1/3 of that of previous state-of-the-art.


title: Detecting a target with quantum entanglement

abstract: In the last decade a lot of research activity focused on the use of quantum entanglement as a resource for remote target detection, i.e. on the design of a quantum radar. The literature on this subject uses tools of quantum optics and quantum information theory, and therefore it often results obscure to radar engineers. This review has been written with the purpose of removing this obscurity.As such, it contains a review of the main advances in the quantum radar literature accompanied by a thorough introduction of the quantum optics background necessary for its understanding.


title: Using Deep Learning to Enhance Event Geometry Reconstruction for the Telescope Array Surface Detector

abstract: The extremely low flux of ultra-high energy cosmic rays (UHECR) makes their direct observation by orbital experiments practically impossible. For this reason all current and planned UHECR experiments detect cosmic rays indirectly by observing the extensive air showers (EAS) initiated by cosmic ray particles in the atmosphere. The world largest statistics of the ultra-high energy EAS events is recorded by the networks of surface stations. In this paper we consider a novel approach for reconstruction of the arrival direction of the primary particle based on the deep convolutional neural network. The latter is using raw time-resolved signals of the set of the adjacent trigger stations as an input. The Telescope Array (TA) Surface Detector (SD) is an array of 507 stations, each containing two layers plastic scintillator with an area of $3$ m$^2$. The training of the model is performed with the Monte-Carlo dataset. It is shown that within the Monte-Carlo simulations, the new approach yields better resolution than the traditional reconstruction method based on the fitting of the EAS front. The details of the network architecture and its optimization for this particular task are discussed.


title: Data-driven Dynamic Multi-objective Optimal Control: An Aspiration-satisfying Reinforcement Learning Approach

abstract: This paper presents an iterative data-driven algorithm for solving dynamic multi-objective (MO) optimal control problems arising in control of nonlinear continuous-time systems. It is first shown that the Hamiltonian functional corresponding to each objective can be leveraged to compare the performance of admissible policies. Hamiltonian-inequalities are then used for which their satisfaction guarantees satisfying the objectives' aspirations. An aspiration-satisfying dynamic optimization framework is then presented to optimize the main objective while satisfying the aspiration of other objectives. Relation to satisficing (good enough) decision-making framework is shown. A Sum-of-Square (SOS) based iterative algorithm is developed to solve the formulated aspiration-satisfying MO optimization. To obviate the requirement of complete knowledge of the system dynamics, a data-driven satisficing reinforcement learning approach is proposed to solve the SOS optimization problem in real-time using only the information of the system trajectories measured during a time interval without having full knowledge of the system dynamics. Finally, two simulation examples are provided to show the effectiveness of the proposed algorithm.


title: Efficacy of the FDA nozzle benchmark and the lattice Boltzmann method for the analysis of biomedical flows in transitional regime

abstract: Flows through medical devices as well as in anatomical vessels despite being at moderate Reynolds number may exhibit transitional or even turbulent character. In order to validate numerical methods and codes used for biomedical flow computations, the U.S. food and drug administration (FDA) established an experimental benchmark, which was a pipe with gradual contraction and sudden expansion representing a nozzle. The experimental results for various Reynolds numbers ranging from 500 to 6500 were publicly released. Previous and recent computational investigations of flow in the FDA nozzle found limitations in various CFD approaches and some even questioned the adequacy of the benchmark itself. This communication reports the results of a lattice Boltzmann method (LBM) based direct numerical simulation (DNS) approach applied to the FDA nozzle benchmark for transitional cases of Reynolds numbers 2000 and 3500. The goal is to evaluate if a simple off the shelf LBM would predict the experimental results without the use of complex models or synthetic turbulence at the inflow. LBM computations with various spatial and temporal resolutions are performed - in the extremities of 44 million to 2.8 billion lattice cells - conducted respectively on 32 CPU cores of a desktop to more than 300'000 cores of a modern supercomputer to explore and characterize miniscule flow details and quantify Kolmogorov scales. The LBM simulations transition to turbulence at a Reynolds number 2000 like the FDA's experiments and acceptable agreement in jet breakdown locations, average velocity, shear stress and pressure is found for both the Reynolds numbers.

categories: physics.flu-dyn cs.DC nlin.CG physics.bio-ph


title: A model theoretic solution to a problem of L\'{a}szl\'{o} Fuchs

abstract: Problem 5.1 in page 181 of [Fuc15] asks to find the cardinals $\lambda$ such that there is a universal abelian $p$-group for purity of cardinality $\lambda$, i.e., an abelian $p$-group $U_\lambda$ of cardinality $\lambda$ such that every abelian $p$-group of cardinality $\leq \lambda$ purely embeds in $U_\lambda$. In this paper we use ideas from the theory of abstract elementary classes to show: $\textbf{Theorem.}$ Let $p$ be a prime number. If $\lambda^{\aleph_0}=\lambda$ or $\forall \mu < \lambda( \mu^{\aleph_0} < \lambda)$, then there is a universal abelian $p$-group for purity of cardinality $\lambda$. Moreover for $n\geq 2$, there is a universal abelian $p$-group for purity of cardinality $\aleph_n$ if and only if $2^{\aleph_0} \leq \aleph_n$. As the theory of abstract elementary classes has barely been used to tackle algebraic questions, an effort was made to introduce this theory from an algebraic perspective.


title: Using Bayesian Optimization to Accelerate Virtual Screening for the Discovery of Therapeutics Appropriate for Repurposing for COVID-19

abstract: The novel Wuhan coronavirus known as SARS-CoV-2 has brought almost unprecedented effects for a non-wartime setting, hitting social, economic and health systems hard.~ Being able to bring to bear pharmaceutical interventions to counteract its effects will represent a major turning point in the fight to turn the tides in this ongoing battle.~ Recently, the World's most powerful supercomputer, SUMMIT, was used to identify existing small molecule pharmaceuticals which may have the desired activity against SARS-CoV-2 through a high throughput virtual screening approach. In this communication, we demonstrate how the use of Bayesian optimization can provide a valuable service for the prioritisation of these calculations, leading to the accelerated identification of high-performing candidates, and thus expanding the scope of the utility of HPC systems for time critical screening


title: Velocity Dispersions of Massive Quiescent Galaxies from Weak Lensing and Spectroscopy

abstract: We use MMT spectroscopy and deep Subaru Hyper Suprime-Cam (HSC) imaging to compare the spectroscopic central stellar velocity dispersion of quiescent galaxies with the effective dispersion of the dark matter halo derived from the stacked lensing signal. The spectroscopic survey (the Smithsonian Hectospec Lensing Survey) provides a sample of 4585 quiescent galaxy lenses with measured line-of-sight central stellar velocity dispersion ($\sigma_{\rm SHELS}$) that is more than 85% complete for $R < 20.6$, $D_{n}4000> 1.5$ and $M_{\star} > 10^{9.5}{\rm M}_{\odot}$. The median redshift of the sample of lenses is 0.32. We measure the stacked lensing signal from the HSC deep imaging. The central stellar velocity dispersion is directly proportional to the velocity dispersion derived from the lensing $\sigma_{\rm Lens}$, $\sigma_{\rm Lens} = (1.05\pm0.15)\sigma_{\rm SHELS}+(-21.17\pm35.19)$. The independent spectroscopic and weak lensing velocity dispersions probe different scales, $\sim3$kpc and $\gtrsim$ 100 kpc, respectively, and strongly indicate that the observable central stellar velocity dispersion for quiescent galaxies is a good proxy for the velocity dispersion of the dark matter halo. We thus demonstrate the power of combining high-quality imaging and spectroscopy to shed light on the connection between galaxies and their dark matter halos.


title: Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

abstract: When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.


title: A convex programming approach to solve posynomial systems

abstract: We exhibit a class of classical or tropical posynomial systems which can be solved by reduction to linear or convex programming problems. This relies on a notion of colorful vectors with respect to a collection of Newton polytopes. This extends the convex programming approach of one player stochastic games.


title: Vaporization frequency response to pressure oscillations: an approximate analytical solution for mixed injection regimes

abstract: This work is devoted to a theoretical analysis of mass frequency response to pressure oscillations of a spray of repetitively injected drops into a combustion chamber. A single stationary spherical droplet continuously fed with the same liquid fuel so that its volume remains constant despite the evaporation, the so-called 'mean droplet' in the Heidmann analogy, represents this vaporizing spray. The feeding is realized with a liquid-liquid heat transfer coefficient by using a source point placed at the mean droplet centre, in such a way that only radial thermal convection and conduction effects are allowed inside the droplet during the process. This feeding procedure is now viewed as a proper boundary condition that is a mixed or a generalized feeding regime controlling the liquid fuel injection process into the combustion chamber. Drawing upon a linear analysis based on the Rayleigh criterion, the evaporating mass response factor is evaluated. Effects due to the variation of the heat transfer coefficient and that of the process characteristic times are analysed. An abrupt increase appears in the response when a fuel thermodynamic coefficient approaches a particular value.


title: Calibration Requirements for Epoch of Reionization 21-cm signal observations -- I. Effect of time-correlated gains

abstract: The residual gain errors add to the systematics of the radio interferometric observations. In case of the high dynamic range observations, these systematic effects dominates over the thermal noise of the observation. In this work, we investigate the effect of time-correlated residual gain errors in the estimation of the power spectrum of the sky brightness distribution in high dynamic range observations. Particularly, we discuss a methodology to estimate the bias in the power spectrum estimator of the redshifted 21-cm signal from neutral hydrogen in the presence of bright extragalactic compact sources. We find, that for the visibility based power spectrum estimators, particularly those use nearby baseline correlations to avoid noise bias, the bias in the power spectrum arises mainly from the time correlation in the residual gain error. The bias also depends on the baseline distribution for a particular observation. Analytical calculations show that the bias is dominant for certain types of baseline pairs used for the visibility correlation. We perform simulated observation of extragalactic compact sources in the presence of residual gain errors with the GMRT like array and estimate the bias in the power spectrum. Our results indicate that in order to estimate the redshifted 21-cm power spectrum, better calibration techniques and estimator development are required.


title: Minimum Race-Time Planning-Strategy for an Autonomous Electric Racecar

abstract: Increasing attention to autonomous passenger vehicles has also attracted interest in an autonomous racing series. Because of this, platforms such as Roborace and the Indy Autonomous Challenge are currently evolving. Electric racecars face the challenge of a limited amount of stored energy within their batteries. Furthermore, the thermodynamical influence of an all-electric powertrain on the race performance is crucial. Severe damage can occur to the powertrain components when thermally overstressed. In this work we present a race-time minimal control strategy deduced from an Optimal Control Problem (OCP) that is transcribed into a Nonlinear Problem (NLP). Its optimization variables stem from the driving dynamics as well as from a thermodynamical description of the electric powertrain. We deduce the necessary first-order Ordinary Differential Equations (ODE)s and form simplified loss models for the implementation within the numerical optimization. The significant influence of the powertrain behavior on the race strategy is shown.


title: Formation of massive globular clusters with dark matter and its implication on dark matter annihilation

abstract: Recent observational studies of $\gamma$-ray emission from massive globular clusters (GCs) have revealed possible evidence of dark matter (DM) annihilation within GCs. It is, however, still controversial whether the emission comes from DM or from milli-second pulsars. We here present the new results of numerical simulations, which demonstrate that GCs with DM can originate from nucleated dwarfs orbiting the ancient MW. The simulated stripped nuclei (i.e., GCs) have the central DM densities ranging from 0.1 to several ${\rm M_\odot pc^{-3}}$, depending on the orbits and the masses of the host dwarf galaxies. However, GCs born outside the central regions of their hosts can have no/little DM after their hosts are destroyed and the GCs become the Galactic halo GCs. These results suggest that only GCs originating from stellar nuclei of dwarfs can possibly have DM. We further calculate the expected $\gamma$-ray emission from these simulated GCs and compare them to observations of $\omega$ Cen. Given the large range of DM densities in the simulated GCs, we suggest that the recent possible detection of DM annihilation from GCs should be more carefully interpreted.


title: Renormalization of Higher-Dimensional Operators from On-shell Amplitudes

abstract: On-shell amplitude methods allow to derive one-loop renormalization effects from just tree-level amplitudes, with no need of loop calculations. We derive a simple formula to obtain the anomalous dimensions of higher-dimensional operators from a product of tree-level amplitudes. We show how this works for dimension-6 operators of the Standard Model, providing explicit examples of the simplicity, elegance and efficiency of the method. Many anomalous dimensions can be calculated from the same Standard Model tree-level amplitude, displaying the attractive recycling aspect of the on-shell method. With this method, it is possible to relate anomalous dimensions that in the Feynman approach arise from very different diagrams, and obtain non-trivial checks of their relative coefficients. We compare our results to those in the literature, where ordinary methods have been applied.


title: AdamMC: A Model Checker for Petri Nets with Transits against Flow-LTL (Full Version)

abstract: The correctness of networks is often described in terms of the individual data flow of components instead of their global behavior. In software-defined networks, it is far more convenient to specify the correct behavior of packets than the global behavior of the entire network. Petri nets with transits extend Petri nets and Flow-LTL extends LTL such that the data flows of tokens can be tracked. We present the tool AdamMC as the first model checker for Petri nets with transits against Flow-LTL. We describe how AdamMC can automatically encode concurrent updates of software-defined networks as Petri nets with transits and how common network specifications can be expressed in Flow-LTL. Underlying AdamMC is a reduction to a circuit model checking problem. We introduce a new reduction method that results in tremendous performance improvements compared to a previous prototype. Thereby, AdamMC can handle software-defined networks with up to 82 switches.


title: Hierarchical Eclipses

abstract: The obscuration of a celestial body that covers another one in the background will be called a ``hierarchical eclipse''. The most obvious case is that a star or a planet will be hidden from sight by the moon during a lunar eclipse. Four objects of the solar system will line up then. We investigate this phenomenon with respect to the region of visibility and periodicity. There exists a parallax field constraining the chances for observation. A historic account from the Middle Ages is preserved that we analyse from different viewing angles. Furthermore, we provide a list of events from 0 to 4000 AD. From this, it is apparent that Jupiter is most often involved in such spectacles because its orbit inclination is small. High-inclination orbits reduce the probability to have a coincidence of an occultation of that object with a lunar eclipse.

categories: physics.pop-ph astro-ph.EP physics.geo-ph physics.hist-ph


title: Real-time and high-throughput Raman signal extraction and processing in CARS hyperspectral imaging

abstract: We present a new collection of processing techniques, collectively "factorized Kramers--Kronig and error correction" (fKK-EC), for (a) Raman signal extraction, (b) denoising, and (c) phase- and scale-error correction in coherent anti-Stokes Raman scattering (CARS) hyperspectral imaging and spectroscopy. These new methods are orders-of-magnitude faster than conventional methods and are capable of real-time performance, owing to the unique core concept: performing all processing on a small basis vector set and using matrix/vector multiplication afterwards for direct and fast transformation of the entire dataset. Experimentally, we demonstrate that a 703026 spectra image of chicken cartilage can be processed in 70 s (approximately 0.1 ms / spectrum), which is approximately 70 times faster than with the conventional workflow (approximately 7.0 ms / spectrum). Additionally, we discuss how this method may be used for machine learning (ML) by re-using the transformed basis vector sets with new data. Using this ML paradigm, the same tissue image was processed (post-training) in approximately 33 s, which is a speed-up of approximately 150 times when compared with the conventional workflow.


title: PENNI: Pruned Kernel Sharing for Efficient CNN Inference

abstract: Although state-of-the-art (SOTA) CNNs achieve outstanding performance on various tasks, their high computation demand and massive number of parameters make it difficult to deploy these SOTA CNNs onto resource-constrained devices. Previous works on CNN acceleration utilize low-rank approximation of the original convolution layers to reduce computation cost. However, these methods are very difficult to conduct upon sparse models, which limits execution speedup since redundancies within the CNN model are not fully exploited. We argue that kernel granularity decomposition can be conducted with low-rank assumption while exploiting the redundancy within the remaining compact coefficients. Based on this observation, we propose PENNI, a CNN model compression framework that is able to achieve model compactness and hardware efficiency simultaneously by (1) implementing kernel sharing in convolution layers via a small number of basis kernels and (2) alternately adjusting bases and coefficients with sparse constraints. Experiments show that we can prune 97% parameters and 92% FLOPs on ResNet18 CIFAR10 with no accuracy loss, and achieve 44% reduction in run-time memory consumption and a 53% reduction in inference latency.


title: Early soft and flexible fusion of EEG and fMRI via tensor decompositions

abstract: Data fusion refers to the joint analysis of multiple datasets which provide complementary views of the same task. In this preprint, the problem of jointly analyzing electroencephalography (EEG) and functional Magnetic Resonance Imaging (fMRI) data is considered. Jointly analyzing EEG and fMRI measurements is highly beneficial for studying brain function because these modalities have complementary spatiotemporal resolution: EEG offers good temporal resolution while fMRI is better in its spatial resolution. The fusion methods reported so far ignore the underlying multi-way nature of the data in at least one of the modalities and/or rely on very strong assumptions about the relation of the two datasets. In this preprint, these two points are addressed by adopting for the first time tensor models in the two modalities while also exploring double coupled tensor decompositions and by following soft and flexible coupling approaches to implement the multi-modal analysis. To cope with the Event Related Potential (ERP) variability in EEG, the PARAFAC2 model is adopted. The results obtained are compared against those of parallel Independent Component Analysis (ICA) and hard coupling alternatives in both simulated and real data. Our results confirm the superiority of tensorial methods over methods based on ICA. In scenarios that do not meet the assumptions underlying hard coupling, the advantage of soft and flexible coupled decompositions is clearly demonstrated.


title: Introducing Coherent MIMO Sensing, a fading-resilient, polarization-independent approach to phase-OTDR

abstract: Nowadays, long distance optical fibre transmission systems use polarization diversity multiplexed signals to enhance transmission performance. Distributed acoustic sensors (DAS) use the same propagation medium ie. single mode optical fibre, and aims at comparable targets such as covering the highest distance with the best signal quality. In the case of sensors, a noiseless transmission enables to monitor a large quantity of mechanical events along the fibre. This paper aims at extending the perspectives of DAS systems with regard to technology breakthroughs introduced in long haul transmission systems over the last decade. We recently developed a sensor interrogation method based on coherent phase-sensitive optical time domain reflectometry ($\phi$-OTDR), with dual polarization multiplexing at the transmitter and polarization diversity at the receiver. We name this technique Coherent-MIMO sensing. A study is performed from a dual-polarization numerical model to compare several sensor interrogation techniques, including Coherent-MIMO. We demonstrate that dual-polarization probing of a fibre sensor makes it insensitive to polarization effects, decreases the risks of false alarms and thus strongly enhances its sensitivity. The simulations results are validated with an experiment, and finally quantitative data are given on the performance increase enabled by Coherent-MIMO sensing.

categories: eess.SP physics.app-ph physics.ins-det physics.optics


title: Latency-Sensitive Web Service Workflows: A Case for a Software-Defined Internet

abstract: The Internet, at large, remains under the control of service providers and autonomous systems. The Internet of Things (IoT) and edge computing provide an increasing demand and potential for more user control for their web service workflows. Network Softwarization revolutionizes the network landscape in various stages, from building, incrementally deploying, and maintaining the environment. Software-Defined Networking (SDN) and Network Functions Virtualization (NFV) are two core tenets of network softwarization. SDN offers a logically centralized control plane by abstracting away the control of the network devices in the data plane. NFV virtualizes dedicated hardware middleboxes and deploys them on top of servers and data centers as network functions. Thus, network softwarization enables efficient management of the system by enhancing its control and improving the reusability of the network services. In this work, we propose our vision for a Software-Defined Internet (SDI) for latency-sensitive web service workflows. SDI extends network softwarization to the Internet-scale, to enable a latency-aware user workflow execution on the Internet.


title: ChewBaccaNN: A Flexible 223 TOPS/W BNN Accelerator

abstract: Binary Neural Networks enable smart IoT devices, as they significantly reduce the required memory footprint and computational complexity while retaining a high network performance and flexibility. This paper presents ChewBaccaNN, a 0.7 mm$^2$ sized binary convolutional neural network (CNN) accelerator designed in GlobalFoundries 22 nm technology. By exploiting efficient data re-use, data buffering, latch-based memories, and voltage scaling, a throughput of 241 GOPS is achieved while consuming just 1.1 mW at 0.4V/154MHz during inference of binary CNNs with up to 7x7 kernels, leading to a peak core energy efficiency of 223 TOPS/W. ChewBaccaNN's flexibility allows to run a much wider range of binary CNNs than other accelerators, drastically improving the accuracy-energy trade-off beyond what can be captured by the TOPS/W metric. In fact, it can perform CIFAR-10 inference at 86.8% accuracy with merely 1.3 $\mu J$, thus exceeding the accuracy while at the same time lowering the energy cost by 2.8x compared to even the most efficient and much larger analog processing-in-memory devices, while keeping the flexibility of running larger CNNs for higher accuracy when needed. It also runs a binary ResNet-18 trained on the 1000-class ILSVRC dataset and improves the energy efficiency by 4.4x over accelerators of similar flexibility. Furthermore, it can perform inference on a binarized ResNet-18 trained with 8-bases Group-Net to achieve a 67.5% Top-1 accuracy with only 3.0 mJ/frame -- at an accuracy drop of merely 1.8% from the full-precision ResNet-18.


title: Analysis and Design Considerations for Achieving the Fundamental Limits of Phase Noise in mmWave Oscillators with On-Chip MEMS Resonator

abstract: Very small electromechanical coupling coefficient in micro-electromechanical systems (MEMS) or acoustic resonators is quite of a concern for oscillator performance, specially at mmWave frequencies. This small coefficient is the manifestation of the small ratio of motional capacitance to static capacitance in the resonators. This work provides a general solution to overcome the problem of relatively high static capacitance at mmWave frequencies and presents analysis and design techniques for achieving extremely low phase noise and a very high figure-of-merit (FoM) in an on-chip MEMS resonator based mmWave oscillator. The proposed analysis and techniques are validated with design and simulation of a 30 GHz oscillator with MEMS resonator having quality factor of 10,000 in 14 nm GF technology. Post layout simulation results show that it achieves a phase noise of -132 dBc/Hz and FoM of 217 dBc/Hz at offset of 1 MHz.


title: On a Subclass of Meromorphic Univalent Functions Involving Hypergeometric Function

abstract: The main object of the present paper is to, introduce the. class of meromorphic univalent functions Involving! hypergeomatrc function .We obtain~ some interesting geometric properties according to coefficient inequality , growth and distortion theorems , radii of starlikeness and convexity for the" functions in our subclass.


title: Certain Classes of Univalent Functions With Negative Coefficients Defined By General Linear Operator

abstract: In this study, a subclass of an univalent function with negative coefficients which is defined by a new general Linear operator have been introduced. The sharp results for coefficients estimators, distortion and closure bounds, Hadamard product, and Neighborhood, and this paper deals with the utilizing of many of the results for classical hypergeometric function, where there can be generalized to m-hypergeometric functions. A subclasses of univalent functions are presented, and it has involving operator which generalizes many well-known. Denote A the class of functions f and we have other results have been studied.


title: Muons in supernovae: implications for the axion-muon coupling

abstract: The high temperature and electron degeneracy attained during a supernova allow for the formation of a large muon abundance within the core of the resulting proto-neutron star. If new pseudoscalar degrees of freedom have large couplings to the muon, they can be produced by this muon abundance and contribute to the cooling of the star. By generating the largest collection of supernova simulations with muons to date, we show that observations of the cooling rate of SN 1987A place strong constraints on the coupling of axion-like particles to muons, limiting the coupling to $g_{a\mu} < 10^{-7.5}~\text{GeV}^{-1}$.


title: Moving from two- to multi-way interactions among binary risk factors on the additive scale

abstract: Many studies have focused on investigating deviations from additive interaction of two dichotomous risk factors on a binary outcome. There is, however, a gap in the literature with respect to interactions on the additive scale of >2 risk factors. In this paper, we present an approach for examining deviations from additive interaction among three on more binary exposures. The relative excess risk due to interaction (RERI) is used as measure of additive interaction. First, we concentrate on three risk factors - we propose to decompose the total RERI to: the RERI owned to the joint presence of all 3 risk factors and the RERI of any two risk factors, given that the third is absent. We then extend this approach, to >3 binary risk factors. For illustration, we use a sample from data from the Greek EPIC cohort and we investigate the association with overall mortality of Mediterranean diet, body mass index (BMI), and, smoking. Our formulae enable better interpretability of any evidence for deviations from additivity owned to more than two risk factors and provide simple ways of communicating such results from a public health perspective by attributing any excess relative risk to specific combinations of these factors.


title: ECAPA-TDNN: Emphasized Channel Attention, Propagation and Aggregation in TDNN Based Speaker Verification

abstract: Current speaker verification techniques rely on a neural network to extract speaker representations. The successful x-vector architecture is a Time Delay Neural Network (TDNN) that applies statistics pooling to project variable-length utterances into fixed-length speaker characterizing embeddings. In this paper, we propose multiple enhancements to this architecture based on recent trends in the related fields of face verification and computer vision. Firstly, the initial frame layers can be restructured into 1-dimensional Res2Net modules with impactful skip connections. Similarly to SE-ResNet, we introduce Squeeze-and-Excitation blocks in these modules to explicitly model channel interdependencies. The SE block expands the temporal context of the frame layer by rescaling the channels according to global properties of the recording. Secondly, neural networks are known to learn hierarchical features, with each layer operating on a different level of complexity. To leverage this complementary information, we aggregate and propagate features of different hierarchical levels. Finally, we improve the statistics pooling module with channel-dependent frame attention. This enables the network to focus on different subsets of frames during each of the channel's statistics estimation. The proposed ECAPA-TDNN architecture significantly outperforms state-of-the-art TDNN based systems on the VoxCeleb test sets and the 2019 VoxCeleb Speaker Recognition Challenge.


title: Eyes-Closed Safety Kernels: Safety for Autonomous Systems Under Loss of Observability

abstract: A framework is presented for handling a potential loss of observability of a dynamical system in a provably-safe way. Inspired by the fragility of data-driven perception systems used by autonomous vehicles, we formulate the problem that arises when a sensing modality fails or is found to be untrustworthy during autonomous operation. We cast this problem as a differential game played between the dynamical system being controlled and the external system factor(s) for which observations are lost. The game is a zero-sum Stackelberg game in which the controlled system (leader) is trying to find a trajectory which maximizes a function representing the safety of the system, and the unobserved factor (follower) is trying to minimize the same function. The set of winning initial configurations of this game for the controlled system represent the set of all states in which safety can be maintained with respect to the external factor, even if observability of that factor is lost. This is the set we refer to as the Eyes-Closed Safety Kernel. In practical use, the policy defined by the winning strategy of the controlled system is only needed to be executed whenever observability of the external system is lost or the system deviates from the Eyes-Closed Safety Kernel due to other, non-safety oriented control schemes. We present a means for solving this game offline, such that the resulting winning strategy can be used for computationally efficient, provably-safe, online control when needed. The solution approach presented is based on representing the game using the solutions of two Hamilton-Jacobi partial differential equations. We illustrate the applicability of our framework by working through a realistic example in which an autonomous car must avoid a dynamic obstacle despite potentially losing observability.