title: A Quantum Solution for Efficient Use of Symmetries in the Simulation of Many-Body Systems

abstract: A many-body Hamiltonian can be block-diagonalized by expressing it in terms of symmetry-adapted basis states. Finding the group orbit representatives of these basis states and their corresponding symmetries is currently a memory/computational bottleneck on classical computers during exact diagonalization. We apply Grover's search in the form of a minimization procedure to solve this problem. Our quantum solution provides an exponential reduction in memory, and a quadratic speedup in time over classical methods. We discuss explicitly the full circuit implementation of Grover minimization as applied to this problem, finding that the oracle only scales as polylog in the size of the group, which acts as the search space. Further, we design an error mitigation scheme that, with no additional qubits, reduces the impact of bit-flip errors on the computation, with the magnitude of mitigation directly correlated with the error rate, improving the utility of the algorithm in the Noisy Intermediate Scale Quantum era.


title: Hybrid-Vehfog: A Robust Approach for Reliable Dissemination of Critical Messages in Connected Vehicles

abstract: Vehicular Ad-hoc Networks (VANET) enable efficient communication between vehicles with the aim of improving road safety. However, the growing number of vehicles in dense regions and obstacle shadowing regions like Manhattan and other downtown areas leads to frequent disconnection problems resulting in disrupted radio wave propagation between vehicles. To address this issue and to transmit critical messages between vehicles and drones deployed from service vehicles to overcome road incidents and obstacles, we proposed a hybrid technique based on fog computing called Hybrid-Vehfog to disseminate messages in obstacle shadowing regions, and multi-hop technique to disseminate messages in non-obstacle shadowing regions. Our proposed algorithm dynamically adapts to changes in an environment and benefits in efficiency with robust drone deployment capability as needed. Performance of Hybrid-Vehfog is carried out in Network Simulator (NS-2) and Simulation of Urban Mobility (SUMO) simulators. The results showed that Hybrid-Vehfog outperformed Cloud-assisted Message Downlink Dissemination Scheme (CMDS), Cross-Layer Broadcast Protocol (CLBP), PEer-to-Peer protocol for Allocated REsource (PrEPARE), Fog-Named Data Networking (NDN) with mobility, and flooding schemes at all vehicle densities and simulation times.


title: Bayesian Anomaly Detection and Classification

abstract: Statistical uncertainties are rarely incorporated in machine learning algorithms, especially for anomaly detection. Here we present the Bayesian Anomaly Detection And Classification (BADAC) formalism, which provides a unified statistical approach to classification and anomaly detection within a hierarchical Bayesian framework. BADAC deals with uncertainties by marginalising over the unknown, true, value of the data. Using simulated data with Gaussian noise, BADAC is shown to be superior to standard algorithms in both classification and anomaly detection performance in the presence of uncertainties, though with significantly increased computational cost. Additionally, BADAC provides well-calibrated classification probabilities, valuable for use in scientific pipelines. We show that BADAC can work in online mode and is fairly robust to model errors, which can be diagnosed through model-selection methods. In addition it can perform unsupervised new class detection and can naturally be extended to search for anomalous subsets of data. BADAC is therefore ideal where computational cost is not a limiting factor and statistical rigour is important. We discuss approximations to speed up BADAC, such as the use of Gaussian processes, and finally introduce a new metric, the Rank-Weighted Score (RWS), that is particularly suited to evaluating the ability of algorithms to detect anomalies.


title: Trajectories of Blocked Community Members: Redemption, Recidivism and Departure

abstract: Community norm violations can impair constructive communication and collaboration online. As a defense mechanism, community moderators often address such transgressions by temporarily blocking the perpetrator. Such actions, however, come with the cost of potentially alienating community members. Given this tradeoff, it is essential to understand to what extent, and in which situations, this common moderation practice is effective in reinforcing community rules. In this work, we introduce a computational framework for studying the future behavior of blocked users on Wikipedia. After their block expires, they can take several distinct paths: they can reform and adhere to the rules, but they can also recidivate, or straight-out abandon the community. We reveal that these trajectories are tied to factors rooted both in the characteristics of the blocked individual and in whether they perceived the block to be fair and justified. Based on these insights, we formulate a series of prediction tasks aiming to determine which of these paths a user is likely to take after being blocked for their first offense, and demonstrate the feasibility of these new tasks. Overall, this work builds towards a more nuanced approach to moderation by highlighting the tradeoffs that are in play.


title: The Resolved Distributions of Dust Mass and Temperature in Local Group Galaxies

abstract: We utilize archival far-infrared maps from the Herschel Space Observatory in four Local Group galaxies (Small and Large Magellanic Clouds, M31, and M33). We model their Spectral Energy Distribution (SED) from 100 to 500 $\mu$m using a single-temperature modified blackbody emission with a fixed emissivity index of $\beta = 1.8$. From the best-fit model, we derive the dust temperature, $T_{\rm d}$, and the dust mass surface density, $\Sigma_{\rm d}$, at 13 parsec resolution for SMC and LMC, and at 167 parsec resolution for all targets. This measurement allows us to build the distribution of dust mass and luminosity as functions of dust temperature and mass surface density. We compare those distribution functions among galaxies and between regions in a galaxy. We find that LMC has the highest mass-weighted average $T_{\rm d}$, while M31 and M33 have the lowest mass-weighted average $T_{\rm d}$. Within a galaxy, star forming regions have higher $T_{\rm d}$ and $\Sigma_{\rm d}$ relative to the overall distribution function, due to more intense heating by young stars and higher gas mass surface density. When we degrade the resolutions to mimic distant galaxies, the mass-weighted mean temperature gets warmer as the resolution gets coarser, meaning the temperature derived from unresolved observation is systematically higher than that in highly resolved observation. As an implication, the total dust mass is lower (underestimated) in coarser resolutions. This resolution-dependent effect is more prominent in clumpy star-forming galaxies (SMC, LMC, and M33), and less prominent in more quiescent massive spiral (M31).


title: The flavor composition of astrophysical neutrinos after 8 years of IceCube: an indication of neutron decay scenario?

abstract: In this work we present an updated study of the flavor composition suggested by astrophysical neutrinos observed by IceCube. The main novelties compared to previous studies are the following: 1) we use the most recent measurements, namely 8 years of throughgoing muons and 7.5 years of High Energy Starting Events (HESE); 2) we consider a broken power law spectrum, in order to be consistent with the observations between 30 TeV and few PeV; 3) we use the throughgoing muon flux to predict the number of astrophysical HESE tracks. We show that accounting for the three previous elements, the result favors surprisingly the hypothesis of neutrinos produced by neutron decay, disfavoring the standard picture of neutrinos from pion decay at 2.0$\sigma$ and the damped muons regime at $2.6 \sigma$, once the atmospheric background is considered. Although the conventional scenario is not yet completely ruled out in the statistically and alternative interpretations are also plausible, such as an energy spectrum characterized by a non trivial shape, this intriguing result may suggest new directions for both theoretical interpretation and experimental search strategies.


title: Constraining the Thermal Properties of Planetary Surfaces using Machine Learning: Application to Airless Bodies

abstract: We present a new method for the determination of the surface properties of airless bodies from measurements of the emitted infrared flux. Our approach uses machine learning techniques to train, validate, and test a neural network representation of the thermophysical behavior of the atmosphereless body given shape model, illumination and observational geometry of the remote sensors. The networks are trained on a dataset of thermal simulations of the emitted infrared flux for different values of surface rock abundance, roughness, and values of the thermal inertia of the regolith and of the rock components. These surrogate models are then employed to retrieve the surface thermal properties by Markov Chain Monte Carlo Bayesian inversion of observed infrared fluxes. We apply the method to the inversion of simulated infrared fluxes of asteroid (101195) Bennu -- according to a geometry of observations similar to those planned for NASA's OSIRIS-REx mission -- and infrared observations of asteroid (25143) Itokawa. In both cases, the surface properties of the asteroid -- such as surface roughness, thermal inertia of the regolith and rock component, and relative rock abundance -- are retrieved; the contribution from the regolith and rock components are well separated. For the case of Itokawa, we retrieve a rock abundance of about 85% for pebbles larger than the diurnal skin depth, which is about 2 cm. The thermal inertia of the rock is found to be lower than the expected value for LL chondrites, indicating that the rocks on Itokawa could be fractured. The average thermal inertia of the surface is around 750 $J s^{-1/2} K^{-1} m^{-2}$ and the measurement of thermal inertia of the regolith corresponds to an average regolith particle diameter of about 10 mm, consistently with in situ measurements as well as results from previous studies.


title: Optimal regularity in time and space for the porous medium equation

abstract: Regularity estimates in time and space for solutions to the porous medium equation are shown in the scale of Sobolev spaces. In addition, higher spatial regularity for powers of the solutions is obtained. Scaling arguments indicate that these estimates are optimal. In the linear limit, the proven regularity estimates are consistent with the optimal regularity of the linear case.


title: All-in-one Relaxion, a unified solution to five BSM puzzles

abstract: We present a unified relaxion solution to the five major outstanding issues in particle physics: the hierarchy problem, dark matter, matter-antimatter asymmetry, neutrino masses and the strong CP problem. The only additional field content in our construction with respect to standard relaxion models is an up-type vector-like fermion pair and three right-handed neutrinos charged under the relaxion shift symmetry. The observed dark matter abundance is generated automatically by oscillations of the relaxion field that begin once it is misaligned from its original stopping point after reheating. The matter-antimatter asymmetry arises from spontaneous baryogenesis induced by the CPT violation due to the rolling of the relaxion after reheating. The CPT violation is communicated to the baryons and leptons via an operator, $\partial_\mu \phi J^\mu$, where $J^\mu$ consists of right-handed neutrino currents arising naturally from a simple neutrino mass model. Finally, the strong CP problem is solved via the Nelson-Barr mechanism, i.e. by imposing CP as a symmetry of the Lagrangian that is broken only spontaneously by the relaxion. The CP breaking is such that although an ${\cal O}(1)$ strong CKM phase is generated, the induced strong CP phase is much smaller, i.e., within experimental bounds.


title: Simultaneous calibration of spectro-photometric distances and the Gaia DR2 parallax zero-point offset with deep learning

abstract: Gaia measures the five astrometric parameters for stars in the Milky Way, but only four of them (positions and proper motion, but not parallax) are well measured beyond a few kpc from the Sun. Modern spectroscopic surveys such as APOGEE cover a large area of the Milky Way disk and we can use the relation between spectra and luminosity to determine distances to stars beyond Gaia's parallax reach. Here, we design a deep neural network trained on stars in common between Gaia and APOGEE that determines spectro-photometric distances to APOGEE stars, while including a flexible model to calibrate parallax zero-point biases in Gaia DR2. We determine the zero-point offset to be $-52.3 \pm 2.0uas$ when modeling it as a global constant, but also train a multivariate zero-point offset model that depends on $G$, $G_{BP} - G_{RP}$ color, and $T_\mathrm{eff}$ and that can be applied to all 139 million stars in Gaia DR2 within APOGEE's color--magnitude range. Our spectro-photometric distances are more precise than Gaia at distances $\approx 2kpc$ from the Sun. We release a catalog of spectro-photometric distances for the entire APOGEE DR14 data set which covers Galactocentric radii $2kpc\lesssim R \lesssim19kpc$; $\approx 150,000$ stars have <10% uncertainty, making this a powerful sample to study the chemo-dynamical structure of the disk. We use this sample to map the mean [Fe/H] and 15 abundance ratios [X/Fe] from the Galactic center to the edge of the disk. Among many interesting trends, we find that the bulge and bar region at $R \lesssim 5kpc$ clearly stands out in [Fe/H] and most abundance ratios.


title: Multiwavelength survey of X-ray sources in the Sculptor Dwarf Spheroidal Galaxy

abstract: We present an unprecedented, deep study of the primordial low-mass X-ray binary population in an isolated, lower-metallicity environment. We perform followup observations of previously-identified X-ray binary candidates in the Sculptor Dwarf Galaxy by combining a second Chandra observation with Spitzer and Gemini photometry, as well as Gemini spectroscopy of selected targets. Of the original nine bright X-ray sources identified, we are able to classify all but one as quasars, active galactic nuclei, or background galaxies. We further discover four new X-ray sources in the second-epoch Chandra observation. Three of these new sources are background sources and one is a foreground flaring star. We have found that Sculptor is effectively devoid of X-ray sources above a few 1e34 erg/s. If Sculptor is able to retain primordial binaries at a similar rate to globular clusters, this implies that bright X-ray binaries observed in globular clusters in the present epoch are all formed dynamically.


title: Effects of spin on magnetized binary neutron star mergers and jet launching

abstract: Events GW170817 and GRB 170817A provide the best confirmation so far that compact binary mergers where at least one of the companions is a neutron star (NS) can be the progenitors of short gamma-ray bursts (sGRBs). An open question for GW170817 remains the values and impact of the initial NS spins. The initial spins could possibly affect the remnant black hole (BH) mass and spin, the remnant disk and the formation and lifetime of a jet and its luminosity. Here we summarize our general relativistic magnetohydrodynamic simulations of spinning, NS binaries undergoing merger and delayed collapse to a BH. The binaries consist of two identical NSs, modeled as $\Gamma=2$ polytropes, in quasicircular orbit, each with spins $\chi_{\rm{NS}}=-0.053,\,0,\,0.24$, or $0.36$. The stars are endowed initially with a dipolar magnetic field extending from the interior into the exterior, as in a radio pulsar. Following merger, the redistribution of angular momentum by magnetic braking and magnetic turbulent viscosity in the hypermassive neutron star (HMNS) remnant, along with the loss of angular momentum due to gravitational radiation, induce the formation of a massive, nearly uniformly rotating inner core surrounded by a magnetized Keplerian disk-like envelope. The HMNS eventually collapses to a BH, with spin $a/M_{\rm BH} \simeq 0.78$ independent of the initial spin of the NSs, surrounded by a magnetized accretion disk. The larger the initial NS spin the heavier the disk. After $\Delta t\sim 3000-4000 M \sim 45-60(M_{\rm NS}/1.625M_\odot)\rm ms$ following merger, a mildly relativistic jet is launched. The lifetime of the jet [$\Delta t\sim 100-140(M_{\rm NS}/1.625M_\odot)\rm ms$] and its outgoing Poynting luminosity [$L_{\rm EM}\sim 10^{51.5\pm 1}\rm erg/s$] are consistent with typical sGRBs, as well as with the Blandford--Znajek mechanism for launching jets and their associated Poynting luminosities.


title: On some properties of multidimensional Bochner-Phillips functional calculus

abstract: The multidimensional functional calculus of semigroup generators, based on the class of Bernstein functions in several variables is developed, the spectral mapping theorems for joint spectra have been stated, the condition for holomorphy of semigroups, generated by operators which arises in the calculus is given, and the moment inequality for such operators is proved.


title: MPP: Model Performance Predictor

abstract: Operations is a key challenge in the domain of machine learning pipeline deployments involving monitoring and management of real-time prediction quality. Typically, metrics like accuracy, RMSE etc., are used to track the performance of models in deployment. However, these metrics cannot be calculated in production due to the absence of labels. We propose using an ML algorithm, Model Performance Predictor (MPP), to track the performance of the models in deployment. We argue that an ensemble of such metrics can be used to create a score representing the prediction quality in production. This in turn facilitates formulation and customization of ML alerts, that can be escalated by an operations team to the data science team. Such a score automates monitoring and enables ML deployments at scale.


title: Learning Hash Function through Codewords

abstract: In this paper, we propose a novel hash learning approach that has the following main distinguishing features, when compared to past frameworks. First, the codewords are utilized in the Hamming space as ancillary techniques to accomplish its hash learning task. These codewords, which are inferred from the data, attempt to capture grouping aspects of the data's hash codes. Furthermore, the proposed framework is capable of addressing supervised, unsupervised and, even, semi-supervised hash learning scenarios. Additionally, the framework adopts a regularization term over the codewords, which automatically chooses the codewords for the problem. To efficiently solve the problem, one Block Coordinate Descent algorithm is showcased in the paper. We also show that one step of the algorithms can be casted into several Support Vector Machine problems which enables our algorithms to utilize efficient software package. For the regularization term, a closed form solution of the proximal operator is provided in the paper. A series of comparative experiments focused on content-based image retrieval highlights its performance advantages.


title: The least prime ideal in the Chebotarev Density Theorem

abstract: In this article, we prove a new bound for the least prime ideal in the Chebotarev density theorem, which improves the main theorem of Zaman [Funct. Approx. Comment. Math. 57 (2017), no.1, 115-142] by a factor of $5/2$. Our main improvement comes from a new version of Tur\'an's power sum method. The key new idea is to use Harnack's inequality for harmonic functions to derive a superior lower bound for the generalised Fej\'er kernel.


title: Seeing beyond the light: Vison and photon electrodynamics in quantum spin ice

abstract: Understanding the nature and behaviour of excitations in quantum spin liquids, and in topological phases of matter in general, is of fundamental importance and has proven crucial for experimental detection and characterisation of candidate materials. Current theoretical and numerical techniques, however, have limited capabilities, especially when it comes to studying gapped excitations. Here, we propose a semiclassical numerical method to study systems whose spin liquid behaviour is underpinned by perturbative ring-exchange Hamiltonians. Our method can readily access both thermodynamic and spectral properties. We focus in particular on quantum spin ice and its photon and vison excitations. After benchmarking the method against existing results on photons, we use it to characterise visons and their thermodynamic behaviour, which remained hitherto largely unexplored. We find that visons, in contrast to spinons in classical spin ice, form a weak electrolyte: vison pairs are the dominant population at low temperatures. This is reflected in the behaviour of thermodynamic quantities, such as pinch point motifs in the relevant correlators. Visons also appear to strongly hybridise with the photon background, a phenomenon that affects the way these quasiparticles may show up in inelastic response measurements. Our results demonstrate that the method, and generalisations thereof, can substantially help our understanding of quasiparticles and their interplay in quantum spin ice and other quantum spin liquids, quantum dimer models, and lattice gauge theories in general.


title: The asymptotic analysis of a Darcy-Stokes system coupled through a curved interface

abstract: The asymptotic analysis of a Darcy-Stokes system modeling the fluid exchange between a narrow channel (Stokes) and a porous medium (Darcy) coupled through a $ C^{2} $ curved interface, is presented. The channel is a cylindrical domain between the interface ($ \Gamma $) and a parallel translation of it ($ \Gamma + \epsilon \, \boldsymbol{\widehat{e}}_{N} $). The introduction of a change variable to fix the domain's geometry and the introduction of two systems of coordinates: the Cartesian and a local one (consistent with the geometry of the surface), permit to find a Darcy-Brinkman lower dimensional coupled system as the limiting form, when the width of the channel tends to zero ($ \epsilon \rightarrow 0 $).


title: Experiments made before the instituted society in confirmation of the theory of the pressures which the walls of the channel sustain due to the flowing water

abstract: This is an annotated translation from Latin of "Experimenta coram societate instituta in confirmationem theoriae pressionum quas latera canalis ab aqua tranfluente sustinet" in which Daniel Bernoulli describes six experiments conducted before the St. Petersburg Academy in 1729, in confirmation of his theory dedicated to the motion of water through ducts, which contains the developments that led to the relation between pressure and velocity in fluid flows, known as Bernoulli's principle.


title: Groups normalized by the odd unitary group

abstract: We will give a definition of quadratic forms on bimodules and prove the sandwich classification theorem for subgroups of the general linear group $\mathrm{GL}(P)$ normalized by the elementary unitary group $\mathrm{EU}(P)$ if $P$ is a nondegenerate bimodule with large enough hyperbolic part.


title: Realizing ergodic properties in zero entropy subshifts

abstract: A subshift with linear block complexity has at most countably many ergodic measures, and we continue of the study of the relation between such complexity and the invariant measures. By constructing minimal subshifts whose block complexity is arbitrarily close to linear but has uncountably many ergodic measures, we show that this behavior fails as soon as the block complexity is superlinear. With a different construction, we show that there exists a minimal subshift with an ergodic measure whose slow entropy grows slower than any given rate tending to infinitely but faster than any other rate majorizing this one yet still growing subexponentially. These constructions lead to obstructions in using subshifts in applications to properties of the prime numbers and in finding a measurable version of the complexity gap that arises for shifts of sublinear complexity.


title: OpenKiwi: An Open Source Framework for Quality Estimation

abstract: We introduce OpenKiwi, a PyTorch-based open source framework for translation quality estimation. OpenKiwi supports training and testing of word-level and sentence-level quality estimation systems, implementing the winning systems of the WMT 2015-18 quality estimation campaigns. We benchmark OpenKiwi on two datasets from WMT 2018 (English-German SMT and NMT), yielding state-of-the-art performance on the word-level tasks and near state-of-the-art in the sentence-level tasks.


title: Better Algorithms for Stochastic Bandits with Adversarial Corruptions

abstract: We study the stochastic multi-armed bandits problem in the presence of adversarial corruption. We present a new algorithm for this problem whose regret is nearly optimal, substantially improving upon previous work. Our algorithm is agnostic to the level of adversarial contamination and can tolerate a significant amount of corruption with virtually no degradation in performance.


title: Scalable Hyperbolic Recommender Systems

abstract: We present a large scale hyperbolic recommender system. We discuss why hyperbolic geometry is a more suitable underlying geometry for many recommendation systems and cover the fundamental milestones and insights that we have gained from its development. In doing so, we demonstrate the viability of hyperbolic geometry for recommender systems, showing that they significantly outperform Euclidean models on datasets with the properties of complex networks. Key to the success of our approach are the novel choice of underlying hyperbolic model and the use of the Einstein midpoint to define an asymmetric recommender system in hyperbolic space. These choices allow us to scale to millions of users and hundreds of thousands of items.


title: Saliency Learning: Teaching the Model Where to Pay Attention

abstract: Deep learning has emerged as a compelling solution to many NLP tasks with remarkable performances. However, due to their opacity, such models are hard to interpret and trust. Recent work on explaining deep models has introduced approaches to provide insights toward the model's behaviour and predictions, which are helpful for assessing the reliability of the model's predictions. However, such methods do not improve the model's reliability. In this paper, we aim to teach the model to make the right prediction for the right reason by providing explanation training and ensuring the alignment of the model's explanation with the ground truth explanation. Our experimental results on multiple tasks and datasets demonstrate the effectiveness of the proposed method, which produces more reliable predictions while delivering better results compared to traditionally trained models.


title: Functions of bounded mean oscillation and Hankel operators on compact abelian groups

abstract: Generalization of functions of bounded mean oscillation and Hankel operators to the case of compact abelian groups with linearly ordered dual is considered. Spaces of functions of bounded mean oscillation and of bounded mean oscillation of analytic type on such groups are described in terms of boundedness of corresponding Hankel operators under the assumption that the dual group contains a minimal positive element.


title: Streaming Low-Rank Matrix Approximation with an Application to Scientific Simulation

abstract: This paper argues that randomized linear sketching is a natural tool for on-the-fly compression of data matrices that arise from large-scale scientific simulations and data collection. The technical contribution consists in a new algorithm for constructing an accurate low-rank approximation of a matrix from streaming data. This method is accompanied by an a priori analysis that allows the user to set algorithm parameters with confidence and an a posteriori error estimator that allows the user to validate the quality of the reconstructed matrix. In comparison to previous techniques, the new method achieves smaller relative approximation errors and is less sensitive to parameter choices. As concrete applications, the paper outlines how the algorithm can be used to compress a Navier--Stokes simulation and a sea surface temperature dataset.


title: Quantum Field Theory and Functional Integrals

abstract: These notes were inspired by the course ''Quantum Field Theory from a Functional Integral Point of View'' given at the University of Zurich in Spring 2017 by Santosh Kandel. We describe Feynman's path integral approach to quantum mechanics and quantum field theory from a functional integral point of view, where the main focus lies in Euclidean field theory. The notion of Gaussian measure and the construction of the Wiener measure are covered. Moreover, we recall the notion of classical mechanics and the Schr\"odinger picture of quantum mechanics, where it shows the equivalence to the path integral formalism, by deriving the quantum mechanical propagator out of it. Additionally, we give an introduction to elements of constructive quantum field theory.


title: Decentralized Coordinate-Descent Data Detection and Precoding for Massive MU-MIMO

abstract: Massive multiuser (MU) multiple-input multiple-output (MIMO) promises significant improvements in spectral efficiency compared to small-scale MIMO. Typical massive MU-MIMO base-station (BS) designs rely on centralized linear data detectors and precoders which entail excessively high complexity, interconnect data rates, and chip input/output (I/O) bandwidth when executed on a single computing fabric. To resolve these complexity and bandwidth bottlenecks, we propose new decentralized algorithms for data detection and precoding that use coordinate descent. Our methods parallelize computations across multiple computing fabrics, while minimizing interconnect and I/O bandwidth. The proposed decentralized algorithms achieve near-optimal error-rate performance and multi-Gbps throughput at sub-1 ms latency when implemented on a multi-GPU cluster with half-precision floating-point arithmetic.


title: What makes a good conversation? How controllable attributes affect human judgments

abstract: A good conversation requires balance -- between simplicity and detail; staying on topic and changing it; asking questions and answering them. Although dialogue agents are commonly evaluated via human judgments of overall quality, the relationship between quality and these individual factors is less well-studied. In this work, we examine two controllable neural text generation methods, conditional training and weighted decoding, in order to control four important attributes for chitchat dialogue: repetition, specificity, response-relatedness and question-asking. We conduct a large-scale human evaluation to measure the effect of these control parameters on multi-turn interactive conversations on the PersonaChat task. We provide a detailed analysis of their relationship to high-level aspects of conversation, and show that by controlling combinations of these variables our models obtain clear improvements in human quality judgments.


title: On the essential spectrum of $\lambda$-Toeplitz operators over compact Abelian groups

abstract: In the recent paper by Mark C. Ho (2014) the notion of a $\lambda$-Toeplitz operator on the Hardy space $H^2(\mathbb{T})$ over the one-dimensional torus $\mathbb{T}$ was introduced and it was shown (under the supplementary condition) that for $\lambda\in \mathbb{T}$ the essential spectrum of such an operator is invariant with respect to the rotation $z\mapsto \lambda z$; if in addition $\lambda$ is not of finite order the essential spectrum is circular. In this paper, we generalize these results to the case when $\mathbb{T}$ is replaced by an arbitrary compact Abelian group whose dual is totally ordered.


title: Principal Component Projection with Low-Degree Polynomials

abstract: In this paper, we consider approximations of principal component projection (PCP) without explicitly computing principal components. This problem has been studied in several recent works. The main feature of existing approaches is viewing the PCP matrix as a matrix function. This underlying function is the composition of a step function with a rational function. To find an approximate PCP, the step function is approximated by a polynomial while the rational function is evaluated by a fast ridge regression solver. In this work, we further improve this process by replacing the rational function with carefully constructed polynomials of low degree. We characterize the properties of polynomials that are suitable for approximating PCP, and establish an optimization problem to select the optimal one from those polynomials. We show theoretically and confirm numerically that the resulting approximate PCP approach with optimal polynomials is indeed effective for approximations of principal component projection.


title: Two-Multicast Channel with Confidential Messages

abstract: Motivated in part by the problem of secure multicast distributed storage, we analyze secrecy rates for a channel in which two transmitters simultaneously multicast to two receivers in the presence of an eavesdropper. Achievable rates are calculated via extensions of a technique due to Chia and El Gamal and the method of output statistics of random binning. Outer bounds are derived for both the degraded and non-degraded versions of the channel, and examples are provided in which the inner and outer bounds meet. The inner bounds recover known results for the multiple-access wiretap channel, broadcast channel with confidential messages, and the compound MAC channel. An auxiliary result is also produced that derives an inner bound on the minimal randomness necessary to achieve secrecy in multiple-access wiretap channels.


title: An SDN-Based Transmission Protocol with In-Path Packet Caching and Retransmission

abstract: In this paper, a comprehensive software-defined networking (SDN) based transmission protocol (SDTP) is presented for fifth generation (5G) communication networks, where an SDN controller gathers network state information from the physical network to improve data transmission efficiency between end hosts, with in-path packet retransmission. In the SDTP, we first develop a new two-way handshake mechanism for connection establishment between a pair of end host. With the aid of SDN control module, signaling exchanges for establishing E2E connections are migrated to the control plane to improve resource utilization in the data plane. A new SDTP packet header format is designed to support efficient data transmission with in-path packet caching and packet retransmission. Based on the new data packet format, a novel in-path receiver-based packet loss detection and caching-based packet retransmission scheme is proposed to achieve in-path fast recovery of lost packets. Extensive simulation results are presented to validate the effectiveness of the proposed protocol in terms of low connection establishment delay and low end-to-end packet transmission delay.


title: Shear-Driven Flow of Athermal, Frictionless, Spherocylinder Suspensions in Two Dimensions: Stress, Jamming, and Contacts

abstract: We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two dimensional spherocylinders driven in uniform steady state shear. Energy dissipation is via a viscous drag with respect to a uniformly sheared host fluid, giving a model for a non-Brownian suspension with a Newtonian rheology. We study pressure $p$ and deviatoric shear stress $\sigma$ as a function of packing fraction $\phi$, strain rate $\dot\gamma$, and a parameter $\alpha$ that measures the asphericity of the particles. We consider the anisotropy of the stress tensor, the macroscopic friction $\mu=\sigma/p$, and the divergence of the transport coefficient $\eta_p=p/\dot\gamma$ as $\phi$ is increased to the jamming $\phi_J$. From an analysis of Herschel-Bulkley rheology above jamming, we estimate $\phi_J$ as a function of $\alpha$ and show that the variation of $\phi_J$ with $\alpha$ is the main cause for differences in rheology as $\alpha$ is varied. However a detailed scaling analysis of the divergence of $\eta_p$ for our most elongated particles suggests that the jamming transition of spherocylinders may be in a different universality class than that of circular disks. We compute the number of contacts per particle $Z$ in the system and show that at jamming $Z_J$ is a non-monotonic function of $\alpha$ that is always smaller than the isostatic value. We measure the probability distribution of contacts per unit surface length $\mathcal{P}(\vartheta)$ at polar angle $\vartheta$ with respect to the spherocylinder spine, and find that as $\alpha\to 0$ this distribution seems to diverge at $\vartheta=\pi/2$, giving a finite limiting probability for contacts on the vanishingly small flat sides of the spherocylinder. Finally we consider the variation of the average contact force as a function of location on the particle surface.


title: A Computational Comparison of Optimization Methods for the Golomb Ruler Problem

abstract: The Golomb ruler problem is defined as follows: Given a positive integer n, locate n marks on a ruler such that the distance between any two distinct pair of marks are different from each other and the total length of the ruler is minimized. The Golomb ruler problem has applications in information theory, astronomy and communications, and it can be seen as a challenge for combinatorial optimization algorithms. Although constructing high quality rulers is well-studied, proving optimality is a far more challenging task. In this paper, we provide a computational comparison of different optimization paradigms, each using a different model (linear integer, constraint programming and quadratic integer) to certify that a given Golomb ruler is optimal. We propose several enhancements to improve the computational performance of each method by exploring bound tightening, valid inequalities, cutting planes and branching strategies. We conclude that a certain quadratic integer programming model solved through a Benders decomposition and strengthened by two types of valid inequalities performs the best in terms of solution time for small-sized Golomb ruler problem instances. On the other hand, a constraint programming model improved by range reduction and a particular branching strategy could have more potential to solve larger size instances due to its promising parallelization features.


title: Learning protein sequence embeddings using information from structure

abstract: Inferring the structural properties of a protein from its amino acid sequence is a challenging yet important problem in biology. Structures are not known for the vast majority of protein sequences, but structure is critical for understanding function. Existing approaches for detecting structural similarity between proteins from sequence are unable to recognize and exploit structural patterns when sequences have diverged too far, limiting our ability to transfer knowledge between structurally related proteins. We newly approach this problem through the lens of representation learning. We introduce a framework that maps any protein sequence to a sequence of vector embeddings --- one per amino acid position --- that encode structural information. We train bidirectional long short-term memory (LSTM) models on protein sequences with a two-part feedback mechanism that incorporates information from (i) global structural similarity between proteins and (ii) pairwise residue contact maps for individual proteins. To enable learning from structural similarity information, we define a novel similarity measure between arbitrary-length sequences of vector embeddings based on a soft symmetric alignment (SSA) between them. Our method is able to learn useful position-specific embeddings despite lacking direct observations of position-level correspondence between sequences. We show empirically that our multi-task framework outperforms other sequence-based methods and even a top-performing structure-based alignment method when predicting structural similarity, our goal. Finally, we demonstrate that our learned embeddings can be transferred to other protein sequence problems, improving the state-of-the-art in transmembrane domain prediction.


title: Sharp solvability criteria for Dirichlet problems of mean curvature type in Riemannian manifolds: non-existence results

abstract: It is well known that the Serrin condition is a necessary condition for the solvability of the Dirichlet problem for the prescribed mean curvature equation in bounded domains of $\mathbb{R}^n$ with certain regularity. In this paper we investigate the sharpness of the Serrin condition for the vertical mean curvature equation in the product $ M^n \times \mathbb{R} $. Precisely, given a $\mathscr{C}^2$ bounded domain $\Omega$ in $M$ and a function $ H = H (x, z) $ continuous in $\overline{\Omega}\times\mathbb{R}$ and non-decreasing in the variable $z$, we prove that the strong Serrin condition $(n-1)\mathcal{H}_{\partial\Omega}(y)\geq n\sup\limits_{z\in\mathbb{R}}|H(y,z)| \ \forall \ y\in\partial\Omega $, is a necessary condition for the solvability of the Dirichlet problem in a large class of Riemannian manifolds within which are the Hadamard manifolds and manifolds whose sectional curvatures are bounded above by a positive constant. As a consequence of our results we deduce Jenkins-Serrin and Serrin type sharp solvability criteria.


title: Dressed Power-law Inflation with Cuscuton

abstract: We study dressed inflation with a cuscuton and find a novel exact power-law solution. It is well known that the conventional power-law inflation is inconsistent with the Planck data. In contrast to this standard lore, we find that power-law inflation with a cuscuton can be reconciled with the Planck data. Moreover, we argue that the cuscuton generally ameliorates inflation models so that predictions are consistent with observations.


title: Quantum structures from association schemes

abstract: Starting from an association scheme induced by a finite group and the corresponding Bose-Mesner algebra we construct quantum Markov chains (QMC), their entangled versions, and interacting Fock spaces (IFS) using the quantum probabilistic approach. Our constructions are based on the intersection numbers and their duals Krien parameters of the schemes with examples focused on regular (distance-regular and distance-transitive) graphs.


title: Application of deconvolution to recover frequency-domain multiplexed detector pulses

abstract: Multiplexing of radiation detectors reduces the number of readout channels, which in turn reduces the number of digitizer input channels for data acquisition. We recently demonstrated frequency domain multiplexing (FDM) of pulse mode radiation detectors using a resonator that converts the detector signal into a damped sinusoid by convolution. The detectors were given unique "tags" by the oscillation frequency of each resonator. The charge collected and the time-of-arrival of the detector pulse were estimated from the corresponding resonator output in the frequency domain. In this paper, we demonstrate a new method to recover the detector pulse from the damped sinusoidal output by deconvolution. Deconvolution converts the frequency-encoded detector signal back to the original detector pulse. We have developed a new prototype FDM system to multiplex organic scintillators based on convolution and deconvolution. Using the new prototype, the charge collected under the anode pulse can be estimated from the recovered pulse with an uncertainty of about 4.4 keVee (keV electron equivalent). The time-of-arrival can be estimated from the recovered pulse with an uncertainty of about 102 ps. We also used a CeBr3 inorganic scintillator to measure the Cs-137 gamma spectrum using the recovered pulses and found a standard deviation of 13.8 keV at 662 keV compared to a standard deviation of 13.5 keV when the original pulses were used. Coincidence measurements with Na-22 using the deconvolved pulses resulted in a timing uncertainty of 617 ps compared to an uncertainty of 603 ps using the original pulses. Pulse shape discrimination was also performed using Cf-252 source and EJ-309 organic scintillator pulses recovered by deconvolution. A figure of merit value of 1.08 was observed when the recovered pulses were used compared to 1.2 for the original pulses.


title: From open learners to open games

abstract: The categories of open learners (due to Fong, Spivak and Tuy\'eras) and open games (due to the present author, Ghani, Winschel and Zahn) bear a very striking and unexpected similarity. The purpose of this short note is to prove that there is a faithful symmetric monoidal functor from the former to the latter, which means that any supervised neural network (without feedback or other complicating features) can be seen as an open game in a canonical way. Roughly, each parameter is controlled by a different player, and the game's best response relation encodes the dynamics of gradient descent. We suggest paths for further work exploiting the link.


title: Pressure-temperature phase diagram of EuRbFe$_4$As$_4$ superconductor

abstract: The pressure dependencies of the magnetic and superconducting transitions, as well as that of the superconducting upper critical field are reported for single crystalline EuRbFe$_4$As$_4$. Resistance measurements were performed under hydrostatic pressures up to 6.21 GPa and in magnetic fields up to 9 T. Zero-field-cool magnetization measurements were performed under hydrostatic pressures up to 1.24 GPa under 20 mT applied field. Superconducting transition temperature, $T_\text c$, up to 6.21 GPa and magnetic transition temperature, $T_\text M$, up to 1.24 GPa were obtained and a pressure-temperature phase diagram was constructed. Our results show that $T_\text c$ is monotonically suppressed upon increasing pressure. $T_\text M$ is linearly increased up to 1.24 GPa. For the studied pressure range, no signs of the crossing of $T_\text M$ and $T_\text c$ lines are observed. The normalized slope of the superconducting upper critical field is gradually suppressed with increasing pressure, which may be due to the continuous change of Fermi-velocity $v_F$ with pressure.


title: Beating SGD Saturation with Tail-Averaging and Minibatching

abstract: While stochastic gradient descent (SGD) is one of the major workhorses in machine learning, the learning properties of many practically used variants are poorly understood. In this paper, we consider least squares learning in a nonparametric setting and contribute to filling this gap by focusing on the effect and interplay of multiple passes, mini-batching and averaging, and in particular tail averaging. Our results show how these different variants of SGD can be combined to achieve optimal learning errors, hence providing practical insights. In particular, we show for the first time in the literature that tail averaging allows faster convergence rates than uniform averaging in the nonparametric setting. Finally, we show that a combination of tail-averaging and minibatching allows more aggressive step-size choices than using any one of said components.


title: A Non-Equilibrium Approach To Holographic Superconductors Using Gradient Flow

abstract: We study a charged scalar field in a bulk 3+1 dimensional anti-deSitter spacetime with a planar black hole background metric. Through the AdS/CFT correspondence this is equivalent to a strongly coupled field theory in 2+1 dimensions describing a superconductor. We use the gradient flow method and solve the flow equations numerically between two fixed points: a vacuum solution and a hairy black hole solution. We study the corresponding flow on the boundary between a normal metal phase and a superconducting phase. We show how the gradient flow moves fields between two fixed points in a way that minimizes the free energy of the system. At the fixed points of the flow the AdS/CFT correspondence provides an equivalence between the Euclidean on-shell action in the bulk and the free energy of the boundary, but it does not tell us about fields away from equilibrium. However, we can formally link static off-shell configurations in the bulk and in the boundary at the same point along the flow. For quasi-static evolution at least, it may be reasonable to think of this link as an extension of the AdS/CFT correspondance.


title: Discriminative Pattern Mining for Breast Cancer Histopathology Image Classification via Fully Convolutional Autoencoder

abstract: Accurate diagnosis of breast cancer in histopathology images is challenging due to the heterogeneity of cancer cell growth as well as of a variety of benign breast tissue proliferative lesions. In this paper, we propose a practical and self-interpretable invasive cancer diagnosis solution. With minimum annotation information, the proposed method mines contrast patterns between normal and malignant images in unsupervised manner and generates a probability map of abnormalities to verify its reasoning. Particularly, a fully convolutional autoencoder is used to learn the dominant structural patterns among normal image patches. Patches that do not share the characteristics of this normal population are detected and analyzed by one-class support vector machine and 1-layer neural network. We apply the proposed method to a public breast cancer image set. Our results, in consultation with a senior pathologist, demonstrate that the proposed method outperforms existing methods. The obtained probability map could benefit the pathology practice by providing visualized verification data and potentially leads to a better understanding of data-driven diagnosis solutions.


title: Optical observations of the nearby galaxy NGC 2366 through narrowband H$\alpha$ and SII filters. Supernova remnants status

abstract: We present detection of 67 HII regions and two optical supernova remnant (SNR) candidates in the nearby irregular galaxy NGC 2366. The SNR candidates were detected by applying [SII]/H$\alpha$ ratio criterion to observations made with the 2-m RCC telescope at Rozhen National Astronomical Observatory in Bulgaria. In this paper we report coordinates, diameters, H$\alpha$ and [SII] fluxes for detected objects across the two fields of view in NGC 2366 galaxy. Using archival XMM-Newton observations we suggest possible X-ray counterparts of two optical SNR candidates. Also, we discard classification of two previous radio SNR candidates in this galaxy, since they appear to be background galaxies.


title: Temperature-dependent Saturation of Weibel-type Instabilities in Counter-streaming Plasmas

abstract: We present the first 2X2V continuum Vlasov-Maxwell simulations of interpenetrating, unmagnetized plasmas to study the competition between two-stream, Oblique, and filamentation modes in the weakly relativistic regime. We find that after nonlinear saturation of the fastest-growing two-stream and Oblique modes, the effective temperature anisotropy, which drives current filament formation via the secular Weibel instability, has a strong dependence on the internal temperature of the counter-streaming plasmas. The effective temperature anisotropy is significantly more reduced in colder than in hotter plasmas, leading to orders of magnitude lower magnetization for colder plasmas. A strong dependence of the energy conversion efficiency of Weibel-type instabilities on internal beam temperature has implications for determining their contribution to the observed magnetization of many astrophysical and laboratory plasmas.


title: Ising-Dropout: A Regularization Method for Training and Compression of Deep Neural Networks

abstract: Overfitting is a major problem in training machine learning models, specifically deep neural networks. This problem may be caused by imbalanced datasets and initialization of the model parameters, which conforms the model too closely to the training data and negatively affects the generalization performance of the model for unseen data. The original dropout is a regularization technique to drop hidden units randomly during training. In this paper, we propose an adaptive technique to wisely drop the visible and hidden units in a deep neural network using Ising energy of the network. The preliminary results show that the proposed approach can keep the classification performance competitive to the original network while eliminating optimization of unnecessary network parameters in each training cycle. The dropout state of units can also be applied to the trained (inference) model. This technique could compress the network in terms of number of parameters up to 41.18% and 55.86% for the classification task on the MNIST and Fashion-MNIST datasets, respectively.


title: Partial Fermionization---Spectral Universality in 1D Repulsive Bose Gases

abstract: Due to the vast growth of the many-body level density with excitation energy, its smoothed form is of central relevance for spectral and thermodynamic properties of interacting quantum systems. We compute the cumulative of this level density for confined one-dimensional continuous systems with repulsive short-range interactions. We show that the crossover from an ideal Bose gas to the strongly correlated, fermionized gas, i.e., partial fermionization, exhibits universal behavior: Systems with very few up to many particles share the same underlying spectral features. In our derivation we supplement quantum cluster expansions with short-time dynamical information. Our nonperturbative analytical results are in excellent agreement with numerics for systems of experimental relevance in cold atom physics, such as interacting bosons on a ring (Lieb-Liniger model) or subject to harmonic confinement. Our method provides predictions for excitation spectra that enable access to finite-temperature thermodynamics in large parameter ranges.


title: Drug-drug interaction prediction based on co-medication patterns and graph matching

abstract: Background: The problem of predicting whether a drug combination of arbitrary orders is likely to induce adverse drug reactions is considered in this manuscript. Methods: Novel kernels over drug combinations of arbitrary orders are developed within support vector machines for the prediction. Graph matching methods are used in the novel kernels to measure the similarities among drug combinations, in which drug co-medication patterns are leveraged to measure single drug similarities. Results: The experimental results on a real-world dataset demonstrated that the new kernels achieve an area under the curve (AUC) value 0.912 for the prediction problem. Conclusions: The new methods with drug co-medication based single drug similarities can accurately predict whether a drug combination is likely to induce adverse drug reactions of interest. Keywords: drug-drug interaction prediction; drug combination similarity; co-medication; graph matching


title: Radar and Communication Co-existence: an Overview

abstract: Increased amounts of bandwidth are required to guarantee both high-quality/high-rate wireless services (4G and 5G) and reliable sensing capabilities such as automotive radar, air traffic control, earth geophysical monitoring and security applications. Therefore, co-existence between radar and communication systems using overlapping bandwidths has been a primary investigation field in recent years. Various signal processing techniques such as interference mitigation, pre-coding or spatial separation, and waveform design allow both radar and communications to share the spectrum. This article reviews recent work on co-existence between radar and communication systems, including signal models, waveform design and signal processing techniques. Our goal is to survey contributions in this area in order to provide a primary starting point for new researchers interested in these problems.


title: Ideals on countable sets: a survey with questions

abstract: An ideal on a set $X$ is a collection of subsets of $X$ closed under the operations of taking finite unions and subsets of its elements. Ideals are a very useful notion in topology and set theory and have been studied for a long time. We present a survey of results about ideals on countable sets and include many open questions.


title: On the validation and use of high-fidelity numerical simulations for gust response analysis

abstract: Specific gust response is considered as one of the most important loads encountered by an aircraft. The Certification Specification (CS) 25, defined by the European Aviation Safety Agency (EASA), and the Federal Aviation Regulations (FAR) 25, defined by the Federal Aviation Administration (FAA), describe the critical gusts that an aircraft must withstand. They must be analyzed for a large range of flight points (Altitude and Equivalent Air speed) and mass configurations. For some load cases, the standard tools could not be accurate enough to correctly predict the gust response and the use of high-fidelity computation could be required. Therefore, ONERA has implemented in its in-house Computational Fluid Dynamics (CFD) code elsA (ONERA-Airbus-Safran property) the capability to compute the high-fidelity aeroelastic gust response, directly in the time-domain, for different discrete gust shapes.This paper presents some recent work achieved at ONERA concerning high-fidelity simulations for gust response. First, a physical validation of the gust response simulation is performed by comparing the results to those obtained experimentally on a scaled model. Second, numerical comparisons are performed using various techniques, in order to model the gust. Finally, an application for generic regional aircraft is shown.


title: Spatial Analysis Made Easy with Linear Regression and Kernels

abstract: Kernel methods are an incredibly popular technique for extending linear models to non-linear problems via a mapping to an implicit, high-dimensional feature space. While kernel methods are computationally cheaper than an explicit feature mapping, they are still subject to cubic cost on the number of points. Given only a few thousand locations, this computational cost rapidly outstrips the currently available computational power. This paper aims to provide an overview of kernel methods from first-principals (with a focus on ridge regression), before progressing to a review of random Fourier features (RFF), a set of methods that enable the scaling of kernel methods to big datasets. At each stage, the associated R code is provided. We begin by illustrating how the dual representation of ridge regression relies solely on inner products and permits the use of kernels to map the data into high-dimensional spaces. We progress to RFFs, showing how only a few lines of code provides a significant computational speed-up for a negligible cost to accuracy. We provide an example of the implementation of RFFs on a simulated spatial data set to illustrate these properties. Lastly, we summarise the main issues with RFFs and highlight some of the advanced techniques aimed at alleviating them.


title: Optical diffraction from isolated nanoparticles

abstract: When subjected to monochromatic incident light a nanoparticle will emit light which then interferes with the incident beam. With sufficient contrast and sufficiently close to the particle this interference pattern may be recorded with a pointed optical fiber in collection mode. It is shown that the analytic dipole model accurately reproduces the observed interference pattern. Using this model and measuring only the lengths of the first two major axes of the observed elliptical fringes we are able to reproduce and quantify the fringe pattern. Importantly, we are able to locate the nanoparticle, with respect to the fibre, using only visible light in a simple experimental setup. For the case described where the image plane is of the order microns above the substrate, hence the fringe number is large, it is shown that the prediction for the particle location and fringe number is insensitive to measurement errors. The phase shift of the scattered wave, a quantity that is notoriously difficult to measure, is easily determined from the theory however it is very sensitive to errors.


title: Influencing factors that determine the usage of the crowd-shipping services

abstract: The objective of this study is to understand how senders choose shipping services for different products, given the availability of both emerging crowd-shipping (CS) and traditional carriers in a logistics market. Using data collected from a US survey, Random Utility Maximization (RUM) and Random Regret Minimization (RRM) models have been employed to reveal factors that influence the diversity of decisions made by senders. Shipping costs, along with additional real-time services such as courier reputations, tracking info, e-notifications, and customized delivery time and location, have been found to have remarkable impacts on senders' choices. Interestingly, potential senders were willing to pay more to ship grocery items such as food, beverages, and medicines by CS services. Moreover, the real-time services have low elasticities, meaning that only a slight change in those services will lead to a change in sender-behavior. Finally, data-science techniques were used to assess the performance of the RUM and RRM models and found to have similar accuracies. The findings from this research will help logistics firms address potential market segments, prepare service configurations to fulfill senders' expectations, and develop effective business operations strategies.


title: The Kalman condition for the boundary controllability of coupled 1-d wave equations

abstract: This paper is devoted to prove the exact controllability of a system of N one-dimensional coupled wave equations when the control is exerted on a part of the boundary by means of one control. We consider the case where the coupling matrix A has distinct eigenvalues. We give a Kalman condition (necessary and sufficient) and give a description, non-optimal in general, of the attainable set.


title: Closed orbit correction at synchrotrons for symmetric and near-symmetric lattices

abstract: This contribution compiles the benefits of lattice symmetry in the context of closed orbit correction. A symmetric arrangement of BPMs and correctors results in structured orbit response matrices of Circulant or block Circulant type. These forms of matrices provide favorable properties in terms of computational complexity, information compression and interpretation of mathematical vector spaces of BPMs and correctors. For broken symmetries, a nearest-Circulant approximation is introduced and the practical advantages of symmetry exploitation are demonstrated with the help of simulations and experiments in the context of FAIR synchrotrons.


title: Discovering Language of the Stocks

abstract: Stock prediction has always been attractive area for researchers and investors since the financial gains can be substantial. However, stock prediction can be a challenging task since stocks are influenced by a multitude of factors whose influence vary rapidly through time. This paper proposes a novel approach (Word2Vec) for stock trend prediction combining NLP and Japanese candlesticks. First, we create a simple language of Japanese candlesticks from the source OHLC data. Then, sentences of words are used to train the NLP Word2Vec model where training data classification also takes into account trading commissions. Finally, the model is used to predict trading actions. The proposed approach was compared to three trading models Buy & Hold, MA and MACD according to the yield achieved. We first evaluated Word2Vec on three shares of Apple, Microsoft and Coca-Cola where it outperformed the comparative models. Next we evaluated Word2Vec on stocks from Russell Top 50 Index where our Word2Vec method was also very successful in test phase and only fall behind the Buy & Hold method in validation phase. Word2Vec achieved positive results in all scenarios while the average yields of MA and MACD were still lower compared to Word2Vec.


title: Nonnegative Weak Solutions of Thin Film Equations Related to Viscous Flows in Cylindrical Geometries

abstract: Motivated by models for thin films coating cylinders in two physical cases proposed by V.I. Kerchman and A.L. Frenkel, we analyze the dynamics of corresponding thin film models. The models are governed by nonlinear, fourth-order, degenerate, parabolic PDEs. We prove, given positive and suitably regular initial data, the existence of weak solutions in all length scales of the cylinder, where all solutions are only local in time. We also prove that given a length constraint on the cylinder, long-time and global in time weak solutions exist. This analytical result is motivated by numerical work on related models in the Ph.D. Thesis of R. Ogrosky in conjunction with multiple further works jointly worked on by combinations of Camassa, Forest, Lee, the first author, Ogrosky, Olander, and Vaughn.


title: RAMHU: A New Robust Lightweight Scheme for Mutual Users Authentication in Healthcare Applications

abstract: Providing a mechanism to authenticate users in healthcare applications is an essential security requirement to prevent both external and internal attackers from penetrating patients' identities and revealing their health data. Many schemes have been developed to provide authentication mechanisms to ensure that only legitimate users are authorized to connect, but these schemes still suffer from vulnerable security. Various attacks expose patients' data for malicious tampering or destruction. Transferring health-related data and information between users and the health centre makes them exposed to penetration by adversaries as they may move through an insecure channel. In addition, previous mechanisms have suffered from the poor protection of users' authentication information. To ensure the protection of patients' information and data, we propose a scheme that authenticates users based on the information of both the device and the legitimate user. In this paper, we propose a Robust Authentication Model for Healthcare Users (RAMHU) that provides mutual authentication between server and clients. This model utilizes an Elliptic Curve Integrated Encryption Scheme (ECIES) and PHOTON to achieve strong security and a good overall performance. RAMHU relies on multi pseudonyms, physical address, and one-time password mechanisms to authenticate legitimate users. Moreover, extensive informal and formal security analysis with the automated validation of Internet security protocols and applications (AVISPA) tool demonstrates that our model offers a high level of security in repelling a wide variety of possible attacks.


title: Weighted integral solvers for elastic scattering by open arcs in two dimensions

abstract: We present a novel approach for the numerical solution of problems of elastic scattering by open arcs in two dimensions. Our methodology relies on the composition of weighted versions of the classical operators associated with Dirichlet and Neumann boundary conditions in conjunction with a certain "open-arc elastic Calder\'on relation" whose validity is demonstrated in this paper on the basis of numerical experiments, but whose rigorous mathematical proof is left for future work. Using this Calder\'on relation in conjunction with spectrally accurate quadrature rules and the Krylov-subspace linear algebra solver GMRES, the proposed overall open-arc elastic solver produces results of high accuracy in small number of iterations---for low and high frequencies alike. A variety of numerical examples in this paper demonstrate the accuracy and efficiency of the proposed methodology.


title: Acting Is Seeing: Navigating Tight Space Using Flapping Wings

abstract: Wings of flying animals can not only generate lift and control torques but also can sense their surroundings. Such dual functions of sensing and actuation coupled in one element are particularly useful for small sized bio-inspired robotic flyers, whose weight, size, and power are under stringent constraint. In this work, we present the first flapping-wing robot using its flapping wings for environmental perception and navigation in tight space, without the need for any visual feedback. As the test platform, we introduce the Purdue Hummingbird, a flapping-wing robot with 17cm wingspan and 12 grams weight, with a pair of 30-40Hz flapping wings driven by only two actuators. By interpreting the wing loading feedback and its variations, the vehicle can detect the presence of environmental changes such as grounds, walls, stairs, obstacles and wind gust. The instantaneous wing loading can be obtained through the measurements and interpretation of the current feedback by the motors that actuate the wings. The effectiveness of the proposed approach is experimentally demonstrated on several challenging flight tasks without vision: terrain following, wall following and going through a narrow corridor. To ensure flight stability, a robust controller was designed for handling unforeseen disturbances during the flight. Sensing and navigating one's environment through actuator loading is a promising method for mobile robots, and it can serve as an alternative or complementary method to visual perception.


title: On specific factors in graphs

abstract: It is well known that if $G = (V, E)$} is a multigraph and $X\subset V$ is a subset of even order, then $G$ contains a spanning forest $H$ such that each vertex from $X$ has an odd degree in $H$ and all the other vertices have an even degree in $H$. This spanning forest may have isolated vertices. If this is not allowed in $H$, then the situation is much more complicated. In this paper, we study this problem and generalize the concepts of even-factors and odd-factors in a unified form.


title: Explicit unconditionally stable methods for the heat equation via potential theory

abstract: We study the stability properties of explicit marching schemes for second-kind Volterra integral equations that arise when solving boundary value problems for the heat equation by means of potential theory. It is well known that explicit finite difference or finite element schemes for the heat equation are stable only if the time step $\Delta t$ is of the order $O(\Delta x^2)$, where $\Delta x$ is the finest spatial grid spacing. In contrast, for the Dirichlet and Neumann problems on the unit ball in all dimensions $d\ge 1$, we show that the simplest Volterra marching scheme, i.e., the forward Euler scheme, is unconditionally stable. Our proof is based on an explicit spectral radius bound of the marching matrix, leading to an estimate that an $L^2$-norm of the solution to the integral equation is bounded by $c_dT^{d/2}$ times the norm of the right hand side. For the Robin problem on the half space in any dimension, with constant Robin (heat transfer) coefficient $\kappa$, we exhibit a constant $C$ such that the forward Euler scheme is stable if $\Delta t < C/\kappa^2$, independent of any spatial discretization. This relies on new lower bounds on the spectrum of real symmetric Toeplitz matrices defined by convex sequences. Finally, we show that the forward Euler scheme is unconditionally stable for the Dirichlet problem on any smooth convex domain in any dimension, in $L^\infty$-norm.


title: Enhancing Clinical Concept Extraction with Contextual Embeddings

abstract: Neural network-based representations ("embeddings") have dramatically advanced natural language processing (NLP) tasks, including clinical NLP tasks such as concept extraction. Recently, however, more advanced embedding methods and representations (e.g., ELMo, BERT) have further pushed the state-of-the-art in NLP, yet there are no common best practices for how to integrate these representations into clinical tasks. The purpose of this study, then, is to explore the space of possible options in utilizing these new models for clinical concept extraction, including comparing these to traditional word embedding methods (word2vec, GloVe, fastText). Both off-the-shelf open-domain embeddings and pre-trained clinical embeddings from MIMIC-III are evaluated. We explore a battery of embedding methods consisting of traditional word embeddings and contextual embeddings, and compare these on four concept extraction corpora: i2b2 2010, i2b2 2012, SemEval 2014, and SemEval 2015. We also analyze the impact of the pre-training time of a large language model like ELMo or BERT on the extraction performance. Last, we present an intuitive way to understand the semantic information encoded by contextual embeddings. Contextual embeddings pre-trained on a large clinical corpus achieves new state-of-the-art performances across all concept extraction tasks. The best-performing model outperforms all state-of-the-art methods with respective F1-measures of 90.25, 93.18 (partial), 80.74, and 81.65. We demonstrate the potential of contextual embeddings through the state-of-the-art performance these methods achieve on clinical concept extraction. Additionally, we demonstrate contextual embeddings encode valuable semantic information not accounted for in traditional word representations.


title: The Generalized Complex Kernel Least-Mean-Square Algorithm

abstract: We propose a novel adaptive kernel based regression method for complex-valued signals: the generalized complex-valued kernel least-mean-square (gCKLMS). We borrow from the new results on widely linear reproducing kernel Hilbert space (WL-RKHS) for nonlinear regression and complex-valued signals, recently proposed by the authors. This paper shows that in the adaptive version of the kernel regression for complex-valued signals we need to include another kernel term, the so-called pseudo-kernel. This new solution is endowed with better representation capabilities in complex-valued fields, since it can efficiently decouple the learning of the real and the imaginary part. Also, we review previous realizations of the complex KLMS algorithm and its augmented version to prove that they can be rewritten as particular cases of the gCKLMS. Furthermore, important conclusions on the kernels design are drawn that help to greatly improve the convergence of the algorithms. In the experiments, we revisit the nonlinear channel equalization problem to highlight the better convergence of the gCKLMS compared to previous solutions. Also, the flexibility of the proposed generalized approach is tested in a second experiment with non-independent real and imaginary parts. The results illustrate the significant performance improvements of the gCKLMS approach when the complex-valued signals have different properties for the real and imaginary parts.


title: Attacking Hardware AES with DFA

abstract: We present the first practical attack on a hardware AES accelerator with 256 bit embedded keys using DFA. We identify the challenges of adapting well-known theoretical AES DFA models to hardware under attack from voltage fault injection and present solutions to those challenges. As a result, we managed to recover 278 real-world AES-256 keys from a secure computing system in a matter of hours with minimal cost.


title: Formation of planetary systems by pebble accretion and migration: How the radial pebble flux determines a terrestrial-planet or super-Earth growth mode

abstract: Super-Earths are found in tighter orbits than the Earth's around more than one third of main sequence stars. It has been proposed that super-Earths are scaled-up terrestrial planets that formed similarly, through mutual accretion of planetary embryos, but in discs much denser than the solar protoplanetary disc. We argue instead that terrestrial planets and super-Earths have two distinct formation pathways that are regulated by the disc's pebble reservoir. Through numerical integrations, which combine pebble accretion and N-body gravity between embryos, we show that a difference of a factor of two in the pebble mass-flux is enough to change the evolution from the terrestrial to the super-Earth growth mode. If the pebble mass-flux is small, then the initial embryos within the ice line grow slowly and do not migrate substantially, resulting in a widely spaced population of Mars-mass embryos when the gas disc dissipates. Without gas being present, the embryos become unstable and a small number of terrestrial planets are formed by mutual collisions. The final terrestrial planets are at most 5 Earth masses. Instead, if the pebble mass-flux is high, then the initial embryos within the ice line rapidly become sufficiently massive to migrate through the gas disc. Embryos concentrate at the inner edge of the disc and growth accelerates through mutual merging. This leads to the formation of a system of closely spaced super-Earths in the 5 to 20 Earth-mass range, bounded by the pebble isolation mass. Generally, instabilities of these super-Earth systems after the disappearance of the gas disc trigger additional merging events and dislodge the system from resonant chains. The pebble flux - which controls the transition between the two growth modes - may be regulated by the initial reservoir of solids in the disc or the presence of more distant giant planets that can halt the radial flow of pebbles.


title: The Donaldson-Thomas partition function of the banana manifold

abstract: A banana manifold is a compact Calabi-Yau threefold, fibered by Abelian surfaces, whose singular fibers have a singular locus given by a "banana configuration of curves". A basic example is given by $X_{ban}$, the blowup along the diagonal of the fibered product of a generic rational elliptic surface $S\to \mathbb{P}^{1}$ with itself. In this paper we give a closed formula for the Donaldson-Thomas partition function of the banana manifold $X_{ban }$ restricted to the 3-dimensional lattice $\Gamma$ of curve classes supported in the fibers of $X_{ban}\to \mathbb{P}^{1}$. It is given by \[ Z_{\Gamma}(X_{ban}) = \prod_{d_{1},d_{2},d_{3}\geq 0} \prod_{k} \left(1-p^{k}Q_{1}^{d_{1}}Q_{2}^{d_{2}}Q_{3}^{d_{3}}\right)^{-12c(||\mathbf{d} ||,k)} \] where $||\mathbf{d} || = 2d_{1}d_{2}+ 2d_{2}d_{3}+ 2d_{3}d_{1}-d_{1}^{2}-d_{2}^{2}-d_{3}^{2}$, and the coefficients $c(a,k)$ have a generating function given by an explicit ratio of theta functions. This formula has interesting properties and is closely realated to the equivariant elliptic genera of $\operatorname{Hilb} (\mathbb{C}^{2})$. In an appendix with S. Pietromonaco, it is shown that the corresponding genus $g$ Gromov-Witten potential $F_{g}$ is a genus 2 Siegel modular form of weight $2g-2$ for $g\geq 2$; namely it is the Skoruppa-Maass lift of a multiple of an Eisenstein series: $\frac{6|B_{2g}|}{g(2g-2)!} E_{2g}(\tau )$.


title: Observational constraints on inflationary potentials within the quantum collapse framework

abstract: The physical mechanism responsible for the emergence of primordial cosmic seeds from a perfect isotropic and homogeneous Universe has not been fully addressed in standard cosmic inflation. To handle this shortcoming, D. Sudarsky et al have developed a proposal: the self-induced collapse hypothesis. In this scheme, the objective collapse of the inflaton's wave function generates the inhomogeneity and anisotropy at all scales. In this paper we analyze the viability of a set of inflationary potentials in both the context of the collapse proposal and within the standard inflationary framework. For this, we perform a statistical analysis using recent CMB and BAO data to obtain the prediction for the scalar spectral index $n_s$ in the context of a particular collapse model: the Wigner scheme. The predicted $n_s$ and the tensor-to-scalar ratio $r$ in terms of the slow roll parameters is different between the collapse scheme and the standard inflationary model. For each potential considered we compare the prediction of $n_s$ and $r$ with the limits established by observational data in both pictures. The result of our analysis shows in most cases a difference in the inflationary potentials allowed by the observational limits in both frameworks. In particular, in the standard approach the more concave a potential is, the more is favored by the data. On the other hand, in the Wigner scheme, the data favors equally all type of concave potentials, including those at the border between convex and concave families.


title: Pragmatic classification of movement primitives for stroke rehabilitation

abstract: Rehabilitation training is the primary intervention to improve motor recovery after stroke, but a tool to measure functional training does not currently exist. To bridge this gap, we previously developed an approach to classify functional movement primitives using wearable sensors and a machine learning (ML) algorithm. We found that this approach had encouraging classification performance but had computational and practical limitations, such as training time, sensor cost, and magnetic drift. Here, we sought to refine this approach and determine the algorithm, sensor configurations, and data requirements needed to maximize computational and practical performance. Motion data had been previously collected from 6 stroke patients wearing 11 inertial measurement units (IMUs) as they moved objects on a target array. To identify optimal ML performance, we evaluated 4 algorithms that are commonly used in activity recognition (linear discriminant analysis (LDA), na\"ive Bayes, support vector machine, and k-nearest neighbors). We compared their classification accuracy, computational complexity, and tuning requirements. To identify optimal sensor configuration, we progressively sampled fewer sensors and compared classification accuracy. To identify optimal data requirements, we compared accuracy using data from IMUs versus accelerometers. We found that LDA had the highest classification accuracy (92%) of the algorithms tested. It also was the most pragmatic, with low training and testing times and modest tuning requirements. We found that 7 sensors on the paretic arm and back resulted in the best accuracy. Using this array, accelerometers had a lower accuracy (84%). We refined strategies to accurately and pragmatically quantify functional movement primitives in stroke patients. We propose that this optimized ML-sensor approach could be a means to quantify training dose after stroke.


title: $\ell_1$-sparsity Approximation Bounds for Packing Integer Programs

abstract: We consider approximation algorithms for packing integer programs (PIPs) of the form $\max\{\langle c, x\rangle : Ax \le b, x \in \{0,1\}^n\}$ where $c$, $A$, and $b$ are nonnegative. We let $W = \min_{i,j} b_i / A_{i,j}$ denote the width of $A$ which is at least $1$. Previous work by Bansal et al. \cite{bansal-sparse} obtained an $\Omega(\frac{1}{\Delta_0^{1/\lfloor W \rfloor}})$-approximation ratio where $\Delta_0$ is the maximum number of nonzeroes in any column of $A$ (in other words the $\ell_0$-column sparsity of $A$). They raised the question of obtaining approximation ratios based on the $\ell_1$-column sparsity of $A$ (denoted by $\Delta_1$) which can be much smaller than $\Delta_0$. Motivated by recent work on covering integer programs (CIPs) \cite{cq,chs-16} we show that simple algorithms based on randomized rounding followed by alteration, similar to those of Bansal et al. \cite{bansal-sparse} (but with a twist), yield approximation ratios for PIPs based on $\Delta_1$. First, following an integrality gap example from \cite{bansal-sparse}, we observe that the case of $W=1$ is as hard as maximum independent set even when $\Delta_1 \le 2$. In sharp contrast to this negative result, as soon as width is strictly larger than one, we obtain positive results via the natural LP relaxation. For PIPs with width $W = 1 + \epsilon$ where $\epsilon \in (0,1]$, we obtain an $\Omega(\epsilon^2/\Delta_1)$-approximation. In the large width regime, when $W \ge 2$, we obtain an $\Omega((\frac{1}{1 + \Delta_1/W})^{1/(W-1)})$-approximation. We also obtain a $(1-\epsilon)$-approximation when $W = \Omega(\frac{\log (\Delta_1/\epsilon)}{\epsilon^2})$.


title: Two-Stage Dual Dynamic Programming with Application to Nonlinear Hydro Scheduling

abstract: We present an approximate method for solving nonlinear control problems over long time horizons, in which the full nonlinear model is preserved over an initial part of the horizon, while the remainder of the horizon is modeled using a linear relaxation. As this approximate problem may still be too large to solve directly, we present a Benders decomposition-based solution algorithm that iterates between solving the nonlinear and linear parts of the horizon. This extends the Dual Dynamic Programming approach commonly employed for optimization of linearized hydro power systems. We prove that the proposed algorithm converges after a finite number of iterations, even when the nonlinear initial stage problems are solved inexactly. We also bound the suboptimality of the split-horizon method with respect to the original nonlinear problem, in terms of the properties of a map between the linear and nonlinear state-input trajectories. We then apply this method to a case study concerning a multiple reservoir hydro system, approximating the nonlinear head effects in the second stage using McCormick envelopes. We demonstrate that near-optimal solutions can be obtained in a shrinking horizon setting when the full nonlinear model is used for only a short initial section of the horizon. For this example, the approach is shown to be more practical than both conventional dynamic programming and a multi-cell McCormick envelope approximation from literature.


title: Strongly Enhanced Gilbert Damping in 3d Transition Metal Ferromagnet Monolayers in Contact with Topological Insulator Bi2Se3

abstract: Engineering Gilbert damping of ferromagnetic metal films is of great importance to exploit and design spintronic devices that are operated with an ultrahigh speed. Based on scattering theory of Gilbert damping, we extend the torque method originally used in studies of magnetocrystalline anisotropy to theoretically determine Gilbert dampings of ferromagnetic metals. This method is utilized to investigate Gilbert dampings of 3d transition metal ferromagnet iron, cobalt and nickel monolayers that are contacted by the prototypical topological insulator Bi2Se3. Amazingly, we find that their Gilbert dampings are strongly enhanced by about one order in magnitude, compared with dampings of their bulks and free-standing monolayers, owing to the strong spin-orbit coupling of Bi2Se3. Our work provides an attractive route to tailoring Gilbert damping of ferromagnetic metallic films by putting them in contact with topological insulators.


title: Systematic uncertainty due to background-gas collisions in trapped-ion optical clocks

abstract: We describe a framework for calculating the frequency shift and uncertainty of trapped-ion optical atomic clocks caused by background-gas collisions, and apply this framework to an $^{27}$Al$^+$ clock to enable a total fractional systematic uncertainty below $10^{-18}$. For this clock, with 38(19) nPa of room temperature H$_2$ background gas, we find that collisional heating generates a non-thermal distribution of motional states with a mean time-dilation shift of order $10^{-16}$ at the end of a 150 ms probe, which is not detected by sideband thermometry energy measurements. However, the contribution of collisional heating to the spectroscopy signal is highly suppressed and we calculate the BGC shift to be $-0.6(2.4)\times 10^{-19}$, where the shift is due to collisional heating time-dilation and the uncertainty is dominated by the worst case $\pm \pi/2$ bound used for collisional phase shift of the $^{27}$Al$^+$ superposition state. We experimentally validate the framework and determine the background-gas pressure in situ using measurements of the rate of collisions that cause reordering of mixed-species ion pairs.


title: pyro: a framework for hydrodynamics explorations and prototyping

abstract: pyro is a Python-based simulation framework designed for ease of implementation and exploration of hydrodynamics methods. It is built in a object-oriented fashion, allowing for the reuse of the core components and fast prototyping of new methods.


title: Retrieval Approach for Determining Surface Susceptibilities and Surface Porosities of a Symmetric Metascreen from Reflection and Transmission Coefficients

abstract: Recently we derived generalized sheet transition conditions (GSTCs) for electromagnetic fields at the surface of a metascreen (a metasurface with a ``fishnet'' structure, i.~e., a periodic array of arbitrary spaced apertures in a relatively impenetrable surface). The parameters in these GSTCs are interpreted as effective surface susceptibilities and surface porosities, which themselves are related to the geometry of the apertures that constitute the metascreen. In this paper, we use these GSTCs to derive the plane-wave reflection ($R$) and transmission ($T$) coefficients of a symmetric metascreen, expressed in terms of these surface parameters. From these equations, we develop a retrieval approach for determining the uniquely defined effective surface susceptibilities and surface porosities that characterize the metascreen from measured or simulated data for the $R$ and $T$ coefficients. We present the retrieved surface parameters for metascreens composed of five different types of apertures (circular holes, square holes, crosses, slots, and a square aperture filled with a high-contrast dielectric). The last example exhibits interesting resonances at frequencies where no resonances exist when the aperture is not filled, which opens up the possibility of designing metasurfaces with unique filtering properties. The retrieved surface parameters are validated by comparing them to other approaches.


title: A `Numbers' Approach to Astronomical Correlations I: Introduction and Application to galaxy Scaling Relations

abstract: We propose a new systematic method of studying correlations between parameters that describe an astronomical (or any) physical system. We recall that behind Dimensionless scaling laws in complex, self-interacting physical objects lies a rigorous theorem of Dimensional analysis, known widely as the Buckingham theorem. Once a {\it catalogue} of properties and forces that define an object or physical system is established, the theorem allows one to select a complete set of Dimensionless quantities or {\it Numbers} on which structure must depend. The internal structure takes the form of a functionally defined manifold in the space of these Numbers. Simple and familiar examples are discussed by way of introduction. Correlations in properties of astronomical objects can be sought either through the constancy of these Numbers or between pairs of the Numbers. In either case, within errors, the functional dependences take on an absolute numerical character. As our principal application, we study a well defined sample of galaxies in order to reveal the implied Tully Fisher and Baryonic Tully Fisher relations. We find that $L\,\propto\,v_{rot}^4$ for the former and $M_b\,\propto\,v_{rot}^3$ for the latter, suggesting that these relations may have different causal origins.


title: A General Framework for Structured Learning of Mechanical Systems

abstract: Learning accurate dynamics models is necessary for optimal, compliant control of robotic systems. Current approaches to white-box modeling using analytic parameterizations, or black-box modeling using neural networks, can suffer from high bias or high variance. We address the need for a flexible, gray-box model of mechanical systems that can seamlessly incorporate prior knowledge where it is available, and train expressive function approximators where it is not. We propose to parameterize a mechanical system using neural networks to model its Lagrangian and the generalized forces that act on it. We test our method on a simulated, actuated double pendulum. We show that our method outperforms a naive, black-box model in terms of data-efficiency, as well as performance in model-based reinforcement learning. We also conduct a systematic study of our method's ability to incorporate available prior knowledge about the system to improve data efficiency.


title: The impact of non-Gaussianity on the error covariance for observations of the Epoch of Reionization 21-cm power spectrum

abstract: Recent simulations show the Epoch of Reionization (EoR) 21-cm signal to be inherently non-Gaussian whereby the error covariance matrix $\mathbf{C}_{ij}$ of the 21-cm power spectrum (PS) contains a trispectrum contribution that would be absent if the signal were Gaussian. Using the binned power spectrum and trispectrum from simulations, here we present a methodology for incorporating these with the baseline distribution and system noise to make error predictions for observations with any radio-interferometric array. Here we consider the upcoming SKA-Low. Non-Gaussianity enhances the errors introducing a positive deviation $\Delta$ relative to the Gaussian predictions. $\Delta$ increases with observation time $t_{\rm obs}$ and saturates as the errors approach the cosmic variance. Considering $t_{\rm obs}=1024$ hours where a $5 \sigma$ detection is possible at all redshifts $7 \le z \le 13$, in the absence of foregrounds we find that the deviations are important at small $k$ where we have $\Delta \sim 40-100 \%$ at $k~\sim 0.04 ~{\rm Mpc}^{-1}$ for some of the redshifts and also at intermediate $k \, (\sim 0.4 ~{\rm Mpc}^{-1})$ where we have $\Delta \sim 200 \%$ at $z=7$. Non-Gaussianity also introduces correlations between the errors in different $k$ bins, and we find both correlations and anticorrelations with the correlation coefficient value spanning $-0.4 \le r_{ij} \le 0.8$. Incorporating the foreground wedge, $\Delta$ continues to be important ($> 50\%$) at $z=7$. We conclude that non-Gaussianity makes a significant contribution to the errors and this is important in the context of the future instruments that aim to achieve high-sensitivity measurements of the EoR 21-cm PS.


title: ABM of osteoblast's mechanotransduction pathway: time patterns of critical events

abstract: Background: Mechanotransduction in bone cells plays a pivotal role in osteoblast differentiation and bone remodelling. Mechanotransduction provides the link between modulation of the extracellular matrix and intracellular actions. By controlling the balance between the intracellular and extracellular domains, the mechanotransduction process determines optimal functionality of the skeletal dynamics, and it is one of the possible causes of osteophatological diseases. Results: Mechanotransduction in a single osteoblast under external mechanical perturbations has been modelled in the agent based framework to reproduce the dynamics of the stochastic reaction diffusion process among molecules in the cytoplasm, nuclear and extracellular domains. The amount of molecules and fluctuations of each molecular class has been analysed in terms of recurrences of critical events. A numerical approach has been developed to invert subordination processes and to extract the directing processes from the molecular signals in order to derive the distribution of recurrence of events. Conclusions: We observed large fluctuations enclose information hidden in the noise which is beyond the dynamic variations of molecular baselines. Studying the system under different parametric conditions and stimuli, the results have shown that the waiting time distributions of each molecule are a signature of the system's dynamics. The behaviours of the molecular waiting times change with the changing of parameters presenting the same variation of patterns for similar interacting molecules and identifying specific alterations for key molecules in the mechanotransduction pathway.


title: Distributionally Robust Reinforcement Learning

abstract: Real-world applications require RL algorithms to act safely. During learning process, it is likely that the agent executes sub-optimal actions that may lead to unsafe/poor states of the system. Exploration is particularly brittle in high-dimensional state/action space due to increased number of low-performing actions. In this work, we consider risk-averse exploration in approximate RL setting. To ensure safety during learning, we propose the distributionally robust policy iteration scheme that provides lower bound guarantee on state-values. Our approach induces a dynamic level of risk to prevent poor decisions and yet preserves the convergence to the optimal policy. Our formulation results in a efficient algorithm that accounts for a simple re-weighting of policy actions in the standard policy iteration scheme. We extend our approach to continuous state/action space and present a practical algorithm, distributionally robust soft actor-critic, that implements a different exploration strategy: it acts conservatively at short-term and it explores optimistically in a long-run. We provide promising experimental results on continuous control tasks.


title: Quantum security of hash functions and property-preservation of iterated hashing

abstract: This work contains two major parts: comprehensively studying the security notions of cryptographic hash functions against quantum attacks and the relationships between them; and revisiting whether Merkle-Damgard and related iterated hash constructions preserve the security properties of the compression function in the quantum setting. Specifically, we adapt the seven notions in Rogaway and Shrimpton (FSE'04) to the quantum setting and prove that the seemingly stronger attack model where an adversary accesses a challenger in quantum superposition does not make a difference. We confirm the implications and separations between the seven properties in the quantum setting, and in addition we construct explicit examples separating an inherently quantum notion called collapsing from several proposed properties. Finally, we pin down the properties that are preserved under several iterated hash schemes. In particular, we prove that the ROX construction in Andreeva et al. (Asiacrypt'07) preserves the seven properties in the quantum random oracle model.