title: Implementation of the Quantum Equivalence Principle

abstract: The quantum equivalence principle says that, for any given point, it is possible to find a quantum coordinate system with respect to which we have definite causal structure in the vicinity of that point. It is conjectured that this principle will play a similar role in the construction of a theory of Quantum Gravity to the role played by the equivalence principle in the construction of the theory of General Relativity. To actually implement the quantum equivalence principle we need a suitable notion of quantum coordinate systems - setting up a framework for these is the main purpose of the present paper. First we introduce a notion of extended states consisting of a superposition of terms (labeled by $u$) where each term corresponds to a manifold, $\mathcal{M}_u$, with fields defined on it. A quantum coordinate system consists of an identification of points between some subsets, $\mathcal{O}_u\subseteq \mathcal{M}_u$, of these manifolds along with a coordinate, $x$, that takes the same value on those points identified. We also introduce a notion of quantum coordinate transformations (which can break the identification map between the manifolds) and show how these can be used to attain definite causal structure in the vicinity of a point. We discuss in some detail how the quantum equivalence principle might form a starting point for an approach to constructing a theory of Quantum Gravity that is analogous to way the equivalence principle is used to construct General Relativity.


title: Traditional Machine Learning for Pitch Detection

abstract: Pitch detection is a fundamental problem in speech processing as F0 is used in a large number of applications. Recent articles have proposed deep learning for robust pitch tracking. In this paper, we consider voicing detection as a classification problem and F0 contour estimation as a regression problem. For both tasks, acoustic features from multiple domains and traditional machine learning methods are used. The discrimination power of existing and proposed features is assessed through mutual information. Multiple supervised and unsupervised approaches are compared. A significant relative reduction of voicing errors over the best baseline is obtained: 20% with the best clustering method (K-means) and 45% with a Multi-Layer Perceptron. For F0 contour estimation, the benefits of regression techniques are limited though. We investigate whether those objective gains translate in a parametric synthesis task. Clear perceptual preferences are observed for the proposed approach over two widely-used baselines (RAPT and DIO).


title: Decomposition of Map Graphs with Applications

abstract: Bidimensionality is the most common technique to design subexponential-time parameterized algorithms on special classes of graphs, particularly planar graphs. The core engine behind it is a combinatorial lemma of Robertson, Seymour and Thomas that states that every planar graph either has a $\sqrt{k}\times \sqrt{k}$-grid as a minor, or its treewidth is $O(\sqrt{k})$. However, bidimensionality theory cannot be extended directly to several well-known classes of geometric graphs. Nevertheless, a relaxation of this lemma has been proven useful for unit disk graphs. Inspired by this, we prove a new decomposition lemma for map graphs. Informally, our lemma states the following. For any map graph $G$, there exists a collection $(U_1,\ldots,U_t)$ of cliques of $G$ with the following property: $G$ either contains a $\sqrt{k}\times \sqrt{k}$-grid as a minor, or it admits a tree decomposition where every bag is the union of $O(\sqrt{k})$ of the cliques in the above collection. The new lemma appears to be a handy tool in the design of subexponential parameterized algorithms on map graphs. We demonstrate its usability by designing algorithms on map graphs with running time $2^{O({\sqrt{k}\log{k}})} \cdot n^{O(1)}$ for the Connected Planar $\cal F$-Deletion problem (that encompasses problems such as Feedback Vertex Set and Vertex Cover). Obtaining subexponential algorithms for Longest Cycle/Path and Cycle Packing is more challenging. We have to construct tree decompositions with more powerful properties and to prove sublinear bounds on the number of ways an optimum solution could "cross" bags in these decompositions. For Longest Cycle/Path, these are the first subexponential-time parameterized algorithms on map graphs. For Feedback Vertex Set and Cycle Packing, we improve upon known $2^{O({k^{0.75}\log{k}})} \cdot n^{O(1)}$-time algorithms on map graphs.


title: The StreetLearn Environment and Dataset

abstract: Navigation is a rich and well-grounded problem domain that drives progress in many different areas of research: perception, planning, memory, exploration, and optimisation in particular. Historically these challenges have been separately considered and solutions built that rely on stationary datasets - for example, recorded trajectories through an environment. These datasets cannot be used for decision-making and reinforcement learning, however, and in general the perspective of navigation as an interactive learning task, where the actions and behaviours of a learning agent are learned simultaneously with the perception and planning, is relatively unsupported. Thus, existing navigation benchmarks generally rely on static datasets (Geiger et al., 2013; Kendall et al., 2015) or simulators (Beattie et al., 2016; Shah et al., 2018). To support and validate research in end-to-end navigation, we present StreetLearn: an interactive, first-person, partially-observed visual environment that uses Google Street View for its photographic content and broad coverage, and give performance baselines for a challenging goal-driven navigation task. The environment code, baseline agent code, and the dataset are available at http://streetlearn.cc


title: Asymptotics of MAP Inference in Deep Networks

abstract: Deep generative priors are a powerful tool for reconstruction problems with complex data such as images and text. Inverse problems using such models require solving an inference problem of estimating the input and hidden units of the multi-layer network from its output. Maximum a priori (MAP) estimation is a widely-used inference method as it is straightforward to implement, and has been successful in practice. However, rigorous analysis of MAP inference in multi-layer networks is difficult. This work considers a recently-developed method, multi-layer vector approximate message passing (ML-VAMP), to study MAP inference in deep networks. It is shown that the mean squared error of the ML-VAMP estimate can be exactly and rigorously characterized in a certain high-dimensional random limit. The proposed method thus provides a tractable method for MAP inference with exact performance guarantees.


title: Neutral island statistics during reionization from 21-cm tomography

abstract: We present the prospects of extracting information about the Epoch of Reionization by identifying the remaining neutral regions, referred to as islands, in tomographic observations of the redshifted 21-cm signal. Using simulated data sets we show that at late times the 21-cm power spectrum is fairly insensitive to the details of the reionization process but that the properties of the neutral islands can distinguish between different reionization scenarios. We compare the properties of these islands with those of ionized bubbles. At equivalent volume filling fractions, neutral islands tend to be fewer in number but larger compared to the ionized bubbles. In addition, the evolution of the size distribution of neutral islands is found to be slower than that of the ionized bubbles and also their percolation behaviour differs substantially. Even though the neutral islands are relatively rare, they will be easier to identify in observations with the low-frequency component of the Square Kilometre Array (SKA-Low) due to their larger size and the lower noise levels at lower redshifts. The size distribution of neutral islands at the late stages of reionization is found to depend on the source properties, such as the ionizing efficiency of the sources and their minimum mass. We find the longest line of sight through a neutral region to be more than 100 comoving Mpc until very late stages (90-95 per cent reionized), which may have relevance for the long absorption trough at $z=5.6-5.8$ in the spectrum of quasar ULAS J0148+0600.


title: Cooling of hypernuclear compact stars: Hartree-Fock models and high-density pairing

abstract: The thermal evolution of hypernuclear compact stars is studied for stellar models constructed on the basis of covariant density functional theory in Hartree and Hartree-Fock approximation. Parametrizations of both types are consistent with the astrophysical mass constraints on compact stars and available hypernuclear data. We discuss the differences of these density functionals and highlight the effects they have on the composition and on the cooling of hypernuclear stars. It is shown that hypernuclear stars computed with density functional models that have a low symmetry energy slope, $L$, are fairly consistent with the cooling data of observed compact stars. The class of stellar models based on larger $L$ values gives rise to the direct Urca process at low densities, which leads to significantly faster cooling. We conjecture high-density pairing for protons and $\Lambda$'s in the $P$-wave channel and provide simple scaling arguments to obtain these gaps. As a consequence the most massive stellar models with masses $1.8 \le M/M_{\odot} \le2$ experience slower cooling by hyperonic dUrca processes which involve $\Lambda$'s and protons.


title: High resolution $\gamma$-ray spectroscopy using GALILEO array

abstract: The GALILEO $\gamma$-ray spectrometer has been constructed at the Legnaro National Laboratory of INFN (LNL-INFN). It can be coupled to advanced ancillary devices which allows nuclear structure studies employing the variety of in-beam $\gamma$-ray spectroscopy methods. Such studies benefit from reactions induced by the intense stable beams delivered by the Tandem-ALPI-PIAVE accelerator complex and by the radioactive beams which will be provided by the SPES facility. In this paper we outline two experiments performed within the experimental campaign at GALILEO coupled to the EUCLIDES Si-ball and the Neutron Wall array. The first one was aimed at spectroscopic studies in A=31 mirror nuclei and the second one at measurements of lifetimes of excited states in nuclei in the vicinity of $^{100}$Sn.


title: Nonlinear Craig Interpolant Generation

abstract: Interpolation-based techniques have become popularized in recent years because of their inherently modular and local reasoning, which can scale up existing formal verification techniques like theorem proving, model-checking, abstraction interpretation, and so on, while the scalability is the bottleneck of these techniques. Craig interpolant generation plays a central role in interpolation-based techniques, and therefore has drawn increasing attentions. In the literature, there are various works done on how to automatically synthesize interpolants for decidable fragments of first-order logic, linear arithmetic, array logic, equality logic with uninterpreted functions (EUF), etc., and their combinations. But Craig interpolant generation for non-linear theory and its combination with the aforementioned theories are still in infancy, although some attempts have been done. In this paper, we first prove that a polynomial interpolant of the form $h(\mathbf{x})>0$ exists for two mutually contradictory polynomial formulas $\phi(\mathbf{x},\mathbf{y})$ and $\psi(\mathbf{x},\mathbf{z})$, with the form $f_1\ge0\wedge\cdots\wedge f_n\ge0$, where $f_i$ are polynomials in $\mathbf{x},\mathbf{y}$ or $\mathbf{x},\mathbf{z}$, and the quadratic module generated by $f_i$ is Archimedean. Then, we show that synthesizing such interpolant can be reduced to solving a semi-definite programming problem (${\rm SDP}$). In addition, we propose a verification approach to assure the validity of the synthesized interpolant and consequently avoid the unsoundness caused by numerical error in ${\rm SDP}$ solving. Finally, we discuss how to generalize our approach to general semi-algebraic formulas.


title: Generalizing Graph Convolutional Neural Networks with Edge-Variant Recursions on Graphs

abstract: This paper reviews graph convolutional neural networks (GCNNs) through the lens of edge-variant graph filters. The edge-variant graph filter is a finite order, linear, and local recursion that allows each node, in each iteration, to weigh differently the information of its neighbors. By exploiting this recursion, we formulate a general framework for GCNNs which considers state-of-the-art solutions as particular cases. This framework results useful to i) understand the tradeoff between local detail and the number of parameters of each solution and ii) provide guidelines for developing a myriad of novel approaches that can be implemented locally in the vertex domain. One of such approaches is presented here showing superior performance w.r.t. current alternatives in graph signal classification problems.


title: Traveling wave solutions for the FPU chain: a constructive approach

abstract: Traveling waves for the FPU chain are constructed by solving the associated equation for the spatial profile $u$ of the wave. We consider solutions whose derivatives $u'$ need not be small, may change sign several times, but decrease at least exponentially. Our method of proof is computer-assisted. Unlike other methods, it does not require that the FPU potential has an attractive (positive) quadratic term. But we currently need to restrict the size of that term. In particular, our solutions in the attractive case are all supersonic.


title: 3D analysis of spatial resolution of MIRO/Rosetta measurements at 67P/CG

abstract: The MIRO instrument's remote sensing capability is integral in constraining water density, temperature and velocity fields in the coma of 67P/Churyumov-Gersimenko. Our aim is to quantify the contribution to the water density from all facets inside and outside the field-of-view (FOV) of MIRO, in both, nadir and limb geometries. This information is crucial for understanding the MIRO derived coma production rates and their relation to the nucleus characteristics, and inherent spatial resolution of the data. This study relies on a detailed 3D nucleus shape model, illumination conditions and the pointing information of the viewing geometry. With these parameters, we can evaluate the relative contribution of water density originating from facets directly inside the MIRO beam as well outside of the beam as a function of distance along the MIRO line-of-sight. We also calculate the ratio of in-beam versus out-of-beam number density. We demonstrate that despite the rather small MIRO field-of-view there is only a small fraction of molecules that originate from facets within the MIRO beam. This is true for nadir, but also translated into the limb observing geometry. The MIRO instrument cannot discriminate active from non-active regions directly from observations. This study also suggests that the beam averaged solar incidence angle, local time and mean normal vectors are not necessary related to molecules within the MIRO beam. These results also illustrate why the 1D spherical Haser model can be applied with relative success for analyzing the MIRO data (and generally any Rosetta measurements). The future possibilities of constraining gas activity distribution on the surface should use 3D codes extracting information from the MIRO spectral line shapes which contain additional information. The presented results are in fact applicable to all relevant instruments onboard Rosetta.


title: On genetic correlation estimation with summary statistics from genome-wide association studies

abstract: Genome-wide association studies (GWAS) have been widely used to examine the association between single nucleotide polymorphisms (SNPs) and complex traits, where both the sample size n and the number of SNPs p can be very large. Recently, cross-trait polygenic risk score (PRS) method has gained extremely popular for assessing genetic correlation of complex traits based on GWAS summary statistics (e.g., SNP effect size). However, empirical evidence has shown a common bias phenomenon that even highly significant cross-trait PRS can only account for a very small amount of genetic variance (R^2 often <1%). The aim of this paper is to develop a novel and powerful method to address the bias phenomenon of cross-trait PRS. We theoretically show that the estimated genetic correlation is asymptotically biased towards zero when complex traits are highly polygenic/omnigenic. When all p SNPs are used to construct PRS, we show that the asymptotic bias of PRS estimator is independent of the unknown number of causal SNPs m. We propose a consistent PRS estimator to correct such asymptotic bias. We also develop a novel estimator of genetic correlation which is solely based on two sets of GWAS summary statistics. In addition, we investigate whether or not SNP screening by GWAS p-values can lead to improved estimation and show the effect of overlapping samples among GWAS. Our results may help demystify and tackle the puzzling "missing genetic overlap" phenomenon of cross-trait PRS for dissecting the genetic similarity of closely related heritable traits. We illustrate the finite sample performance of our bias-corrected PRS estimator by using both numerical experiments and the UK Biobank data, in which we assess the genetic correlation between brain white matter tracts and neuropsychiatric disorders.


title: Higher-order corrections to jet quenching

abstract: We compute the inclusive jet spectrum in the presence of a dense QCD medium by going beyond the single parton energy loss approximation. We show that higher-order corrections are important yielding large logarithmic contributions that must be resummed to all orders. This reflects the fact that jet quenching is sensitive to fluctuations of the jet substructure.


title: Patient-Specific Cerebral Aneurysm Hemodynamics: Comparison of in vitro Volumetric Particle Velocimetry, Computational Fluid Dynamics (CFD), and in vivo 4D Flow MRI

abstract: Typical approaches to patient-specific hemodynamic studies of cerebral aneurysms use image based computational fluid dynamics (CFD) and seek to statistically correlate parameters such as wall shear stress (WSS) and oscillatory shear index (OSI) to risk of growth and rupture. However, such studies have reported contradictory results, emphasizing the need for in-depth comparisons of volumetric experiments and CFD. In this work, we conducted tomographic particle velocimetry experiments using two patient-specific cerebral aneurysm models (basilar tip and internal carotid artery) under pulsatile flow conditions and processed the particle images using Shake the Box (STB), a particle tracking method. The STB data was compared to that obtained from in vivo 4D flow MRI and CFD. Although qualitative agreement of flow pathlines across modalities was observed, each modality maintained notably unique spatiotemporal distributions of low normalized WSS regions. Analysis of time averaged WSS (TAWSS), OSI, and Relative Residence Time (RRT) demonstrated that non-dimensional parameters, such as OSI, may be more robust to the varying assumptions, limitations, and spatial resolutions of each subject and modality. These results suggest a need for further multi-modality analysis as well as development of non-dimensional hemodynamic parameters and correlation of such metrics to aneurysm risk of growth and rupture.


title: Geometric phases in 2D and 3D polarized fields: geometrical, dynamical, and topological aspects

abstract: Geometric phases are a universal concept that underpins numerous phenomena involving multi-component wave fields. These polarization-dependent phases are inherent in interference effects, spin-orbit interaction phenomena, and topological properties of vector wave fields. Geometric phases have been thoroughly studied in two-component fields, such as two-level quantum systems or paraxial optical waves. However, their description for fields with three or more components, such as generic nonparaxial optical fields routinely used in modern nano-optics, constitutes a nontrivial problem. Here we describe geometric, dynamical, and total phases calculated along a closed spatial contour in a multi-component complex field, with particular emphasis on 2D (paraxial) and 3D (nonparaxial) optical fields. We present several equivalent approaches: (i) an algebraic formalism, universal for any multi-component field; (ii) a dynamical approach using the Coriolis coupling between the spin angular momentum and reference-frame rotations; and (iii) a geometric representation, which unifies the Pancharatnam-Berry phase for the 2D polarization on the Poincar\'e sphere and the Majorana-sphere representation for the 3D polarized fields. Most importantly, we reveal close connections between geometric phases, angular-momentum properties of the field, and topological properties of polarization singularities in 2D and 3D fields, such as C-points and polarization M\"obius strips.


title: Filippov flows and mean-field limits in the kinetic singular Kuramoto model

abstract: The agent-based singular Kuramoto model was proposed in [60] as a singular version of the Kuramoto model of coupled oscillators that is consistent with Hebb's rule of neuroscience. In such paper, the authors studied its well-posedness via the concept of Filippov solutions. Interestingly, they found some new emergent phenomena in the paradigm of Kuramoto model: clustering into subgroups and emergence of global phase synchronization taking place at finite time. This paper aims at introducing the associated kinetic singular Kuramoto model along $\mathbb{T}\times\mathbb{R}$, that is reminiscent of the classical Kuramoto-Sakaguchi equation. Our main goal is to propose a well-posedness theory of measure-valued solutions that remains valid after eventual phase collisions. The results will depend upon the specific regime of singularity: subcritical, critical and supercritical. The cornerstone is the existence of Filippov characteristic flows for interaction kernels with jump discontinuities. Our second goal is to study stability with respect to initial data, that in particular will provide quantitative estimates for the mean-field limit of the agent-based model towards the kinetic equation in quadratic Wasserstein-type distances. Finally, we will recover global phase synchronization at the macroscopic scale under appropriate generic assumptions on the initial data. The most singular regime will be tackled separately with an alternative method, namely, a singular hyperbolic limit on a regularized second order model with inertia. Note that we will work within the manifold $\mathbb{T}\times\mathbb{R}$ and will avoid resorting on Euclidean approximations. Also, no gradient-type structure will be needed, as opposed to some preceding literature.


title: Long-tail Relation Extraction via Knowledge Graph Embeddings and Graph Convolution Networks

abstract: We propose a distance supervised relation extraction approach for long-tailed, imbalanced data which is prevalent in real-world settings. Here, the challenge is to learn accurate "few-shot" models for classes existing at the tail of the class distribution, for which little data is available. Inspired by the rich semantic correlations between classes at the long tail and those at the head, we take advantage of the knowledge from data-rich classes at the head of the distribution to boost the performance of the data-poor classes at the tail. First, we propose to leverage implicit relational knowledge among class labels from knowledge graph embeddings and learn explicit relational knowledge using graph convolution networks. Second, we integrate that relational knowledge into relation extraction model by coarse-to-fine knowledge-aware attention mechanism. We demonstrate our results for a large-scale benchmark dataset which show that our approach significantly outperforms other baselines, especially for long-tail relations.


title: Magnetic ordering and ferroelectricity in multiferroic 2H-AgFeO2: Comparison between hexagonal and rhombohedral polytypes

abstract: Magnetic and dielectric properties of the hexagonal triangular lattice antiferromagnet 2H-AgFeO2 have been studied by neutron diffraction, magnetic susceptibility, specific heat, pyroelectric current, and dielectric constant measurements. The ferroelectric polarization, P ~ 5 {\mu}C/m2, has been found to appear below 11 K due to a polar nature of the magnetic ground state of the system. In the temperature range of 11 K < T < 18 K, an incommensurate spin density wave (ICM1) with the nonpolar magnetic point group mmm1' and the k1 = (0,q1_b,0; q1_b = 0.390-0.405)propagation vector takes place. Below 14 K, a proper screw ordering (ICM2) and k2 = (0,q2_b,0; q2_b = 0.385-0.396) appears as a minor phase which coexists with ICM1 and the ground state down to the lowestmeasured temperature 5.5 K. No ferroelectric polarization associated with the ICM2 phase was observed in agreement with its nonpolar point group 2221'. Finally, a spiral order with cycloid and proper screw components (ICM3), and k3 = (q3_a,q3_b,0; q3_a = 0.0467, q3_b = 0.349) emerges below 11 K as the ground state of the system. Based on the deduced magnetic point group 21', we conclude that the ferroelectric polarization in ICM3 is parallel to the c axis and is caused by the inverse Dzyloshinskii-Moriya effect with p1 ~ rij x (Si x Sj ). Unlike the rhombohedral 3R-AgFeO2 polytype, the additional contribution to the macroscopic polarization p2 ~ Si x Sj is not allowed in the present case due to the symmetry constraints imposed by the hexagonal lattice of 2H-AgFeO2.


title: Ultra-high-speed Terahertz Imaging Using Atomic Vapour

abstract: Terahertz (THz) technologies, generally defined as operating in the 0.1-10THz range, bridge the gap between electronic and photonic devices. Because THz radiation passes readily through materials such as plastics, paper and cloth it can be employed in non-destructive testing, and as it is non-ionising it is considered safe for security and biomedical applications. There is significant demand for high speed THz imaging across a range of applications but, despite ongoing efforts, fast full-field imaging remains an unfulfilled goal. Here we demonstrate a THz imaging system based upon efficient THz-to-optical conversion in atomic vapour, where full-field images can be collected at ultra-high speeds using conventional optical camera technology. For a 0.55 THz field we show an effective 1 cm $^{2}$ sensor with near diffraction-limited spatial resolution and a minimum detectable power of 190 $\pm$ 30 fW s$^{-1/2}$ per 40x40${\mu}$m pixel capable of video capture at 3000 frames per second. This combination of speed and sensitivity represents a step change in the state of the art of THz imaging, and will likely lead to its uptake in wider industrial settings. With further improvements we expect that even higher frame rates of up to 1 MHz would be possible.


title: Phase Portraits of Hyperbolic Geometry

abstract: Phase plotting is a useful way of visualising functions on complex space. We reinvent the method in the context of hyperbolic geometry, and we use it to plot functions on various representative surfaces for hyperbolic space, illustrating with direct motions in particular. The reinvention is nontrivial, and we discuss the essential features for ensuring that such visualisations convey useful information. Our approach is to exploit conformal maps between representative surfaces, in order to uniquely represent the preimages of geodesics. Our considerations and methods are prototypical of what one might consider for adapting similar methods of visualisation in other contexts.


title: Time Series Source Separation using Dynamic Mode Decomposition

abstract: The Dynamic Mode Decomposition (DMD) extracted dynamic modes are the non-orthogonal eigenvectors of the matrix that best approximates the one-step temporal evolution of the multivariate samples. In the context of dynamical system analysis, the extracted dynamic modes are a generalization of global stability modes. We apply DMD to a data matrix whose rows are linearly independent, additive mixtures of latent time series. We show that when the latent time series are uncorrelated at a lag of one time-step then, in the large sample limit, the recovered dynamic modes will approximate, up to a column-wise normalization, the columns of the mixing matrix. Thus, DMD is a time series blind source separation algorithm in disguise, but is different from closely related second order algorithms such as the Second-Order Blind Identification (SOBI) method and the Algorithm for Multiple Unknown Signals Extraction (AMUSE). All can unmix mixed stationary, ergodic Gaussian time series in a way that kurtosis-based Independent Components Analysis (ICA) fundamentally cannot. We use our insights on single lag DMD to develop a higher-lag extension, analyze the finite sample performance with and without randomly missing data, and identify settings where the higher lag variant can outperform the conventional single lag variant. We validate our results with numerical simulations, and highlight how DMD can be used in change point detection.


title: Masses and electric charges: gauge anomalies and anomalous thresholds

abstract: We work out in the forward limit and up to order $e^6$ in perturbation theory the collinear divergences. In this kinematical regime we discover new collinear divergences that we argue can be only cancelled using quantum interference with processes contributing to the gauge anomaly. This rules out the possibility of a quantum consistent and anomaly free theory with massless charges and long range interactions. We use the anomalous threshold singularities to derive a gravitational lower bound on the mass of the lightest charged fermion.


title: On the spectrum of asymptotic entropies of random walks

abstract: Given a random walk on a free group, we study the random walks it induces on the group's quotients. We show that the spectrum of asymptotic entropies of the induced random walks has no isolated points, except perhaps its maximum.


title: Simson Identity of Generalized m-step Fibonacci Numbers

abstract: One of the best known and oldest identities for the Fibonacci sequence $F_n$ is $F_{n+1}F_{n-1}-F_{n}^2=(-1)^n$ which was derived first by R. Simson in 1753 and it is now called as Simson or Cassini Identity. In this paper, we generalize this result to generalized m-step Fibonacci numbers and give an attractive formula. Furthermore, we present some Simson's identities of particular generalized m-step Fibonacci sequences.


title: Denial-of-Service Attacks on Shared Cache in Multicore: Analysis and Prevention

abstract: In this paper we investigate the feasibility of denial-of-service (DoS) attacks on shared caches in multicore platforms. With carefully engineered attacker tasks, we are able to cause more than 300X execution time increases on a victim task running on a dedicated core on a popular embedded multicore platform, regardless of whether we partition its shared cache or not. Based on careful experimentation on real and simulated multicore platforms, we identify an internal hardware structure of a non-blocking cache, namely the cache writeback buffer, as a potential target of shared cache DoS attacks. We propose an OS-level solution to prevent such DoS attacks by extending a state-of-the-art memory bandwidth regulation mechanism. We implement the proposed mechanism in Linux on a real multicore platform and show its effectiveness in protecting against cache DoS attacks.


title: On the index of reducibility of parameter ideals: The stable and limit values

abstract: Let $(R, \frak m)$ be a Noetherian local ring and $M$ a finitely generated $R$-module of dimension $d$. A famous result of Northcott says that if $M$ is Cohen-Macaulay, then the index of reducibility of parameter ideals on $M$ is an invariant of the module. The aim of this paper is to extend Northcott's theorem for any finitely generated $R$-module. We call this invariant the stable value of the indices of reducibility of parameter ideals of $M$. We also introduce the limit value of the indices of reducibility of parameter ideals of $M$.


title: Stochastic nonlinear model for somatic cell population dynamics during ovarian follicle activation

abstract: In mammals, female germ cells are sheltered within somatic structures called ovarian follicles, which remain in a quiescent state until they get activated, all along reproductive life. We investigate the sequence of somatic cell events occurring just after follicle activation, starting by the awakening of precursor somatic cells, and their transformation into proliferative cells. We introduce a nonlinear stochastic model accounting for the joint dynamics of the two cell types, and allowing us to investigate the potential impact of a feedback from proliferative cells onto precursor cells. To tackle the key issue of whether cell proliferation is concomitant or posterior to cell awakening, we assess both the time needed for all precursor cells to awake, and the corresponding increase in the total cell number with respect to the initial cell number. Using the probabilistic theory of first passage times, we design a numerical scheme based on a rigorous Finite State Projection and coupling techniques to compute the mean extinction time and the cell number at extinction time. We find that the feedback term clearly lowers the number of proliferative cells at the extinction time. We calibrate the model parameters using an exact likelihood approach. We carry out a comprehensive comparison between the initial model and a series of submodels, which helps to select the critical cell events taking place during activation, and suggests that awakening is prominent over proliferation.


title: Theoretical Study on Anisotropic Magnetoresistance Effects of I//[100], I//[110], and I//[001] for Ferromagnets with A Crystal Field of Tetragonal Symmetry

abstract: Using the electron scattering theory, we obtain analytic expressions for anisotropic magnetoresistance (AMR) ratios for ferromagnets with a crystal field of tetragonal symmetry. Here, a tetragonal distortion exists in the [001] direction, the magnetization ${\mbox{\boldmath $M$}}$ lies in the (001) plane, and the current ${\mbox{\boldmath $I$}}$ flows in the [100], [010], or [001] direction. When the ${\mbox{\boldmath $I$}}$ direction is denoted by $i$, we obtain the AMR ratio as ${\rm AMR}^i (\phi_i)= C_0^i + C_2^i \cos 2\phi_i + C_4^i \cos 4 \phi_i \ldots = \sum_{j=0,2,4,\ldots} C_j^i \cos j\phi_i$, with $i=[100]$, $[110]$, and $[001]$, $\phi_{[100]} = \phi_{[001]}=\phi$, and $\phi_{[110]}=\phi'$. The quantity $\phi$ ($\phi'$) is the relative angle between ${\mbox{\boldmath $M$}}$ and the $[100]$ ($[110]$) direction, and $C_j^i$ is a coefficient composed of a spin--orbit coupling constant, an exchange field, the crystal field, and resistivities. We elucidate the origin of $C_j^i \cos j\phi_i$ and the features of $C_j^i$. In addition, we obtain the relation $C_4^{[100]} = -C_4^{[110]}$, which was experimentally observed for Ni, under a certain condition. We also qualitatively explain the experimental results of $C_2^{[100]}$, $C_4^{[100]}$, $C_2^{[110]}$, and $C_4^{[110]}$ at 293 K for Ni.


title: A New, Larger Sample of Supernova Remnants in NGC 6946

abstract: The relatively nearby spiral galaxy NGC~6946 is one of the most actively star forming galaxies in the local Universe. Ten supernovae (SNe) have been observed since 1917, and hence NGC6946 surely contains a large number of supernova remnants (SNRs). Here we report a new optical search for these SNRs using narrow-band images obtained with the WIYN telescope. We identify 147 emission nebulae as likely SNRs, based on elevated [SII]:Halpha ratios compared to HII regions. We have obtained spectra of 102 of these nebulae with Gemini North-GMOS; of these, 89 have [SII]:Halpha ratios greater than 0.4, the canonical optical criterion for identifying SNRs. There is very little overlap between our sample and the SNR candidates identified by Lacey et al. (2001) from radio data. Also, very few of our SNR candidates are known X-ray sources, unlike the situation in some other galaxies such as M33 and M83. The emission line ratios, e.g., [NII]:Halpha, of the candidates in NGC6946 are typical of those observed in SNR samples from other galaxies with comparable metallicity. None of the candidates observed in our low-resolution spectra show evidence of anomalous abundances or significant velocity broadening. A search for emission at the sites of all the historical SNe in NGC6946 resulted in detections for only two: SN1980K and SN2004et. Spectra of both show very broad, asymmetric line profiles, consistent with the interaction between SN ejecta and the progenitor star's circumstellar material, as seen in late spectra from other core-collapse SNe of similar age.


title: Magnetic ordering in pressure-induced phases with giant spin-driven ferroelectricity in multiferroic TbMnO3

abstract: In order to clarify the mechanism associated with pressure/magnetic-field-induced giant ferroelectric polarization in TbMnO3, this work investigated changes in magnetic ordering brought about by variations in temperature, magnetic field, and pressure. This was accomplished by means of neutron diffraction analyses under high pressures and high magnetic fields, employing a single crystal. The incommensurate magnetic ordering of a cycloid structure was found to be stable below the reported critical pressure of 4.5 GPa. In contrast, a commensurate E-type spin ordering of Mn spins and a noncollinear configuration of Tb spins with k=(0,1/2,0) appeared above 4.5 GPa. The application of a magnetic field along the a axis (H_{||a}) under pressure induces a k=(0,0,0)antiferromagnetic structure in the case of Tb spins above H_{||a}, enhancing the ferroelectric polarization, while the E-type ordering of Mn spins is stable even above the critical field. From the present experimental findings, we conclude that the E-type ordering of Mn spins induces giant ferroelectric polarization through an exchange striction mechanism. The H_{||a}-induced polarization enhancement can be understood by considering that the polarization, reduced by the polar ordering of Tb moments in a zero field, can be recovered through a field-induced change to nonpolar k=(0,0,0) ordering at H_{||a} ~ 2T.


title: Algorithms for Piecewise Constant Signal Approximations

abstract: We consider the problem of finding optimal piecewise constant approximations of one-dimensional signals. These approximations should consist of a specified number of segments (samples) and minimise the mean squared error to the original signal. We formalise this goal as a discrete nonconvex optimisation problem, for which we study two algorithms. First we reformulate a recent adaptive sampling method by Dar and Bruckstein in a compact and transparent way. This allows us to analyse its limitations when it comes to violations of its three key assumptions: signal smoothness, local linearity, and error balancing. As a remedy, we propose a direct optimisation approach which does not rely on any of these assumptions and employs a particle swarm optimisation algorithm. Our experiments show that for nonsmooth signals or low sample numbers, the direct optimisation approach offers substantial qualitative advantages over the Dar--Bruckstein method. As a more general contribution, we disprove the optimality of the principle of error balancing for optimising data in the l^2 norm.


title: Adaptive computation of the Symmetric Nonnegative Matrix Factorization (NMF)

abstract: Nonnegative Matrix Factorization (NMF), first proposed in 1994 for data analysis, has received successively much attention in a great variety of contexts such as data mining, text clustering, computer vision, bioinformatics, etc. In this paper the case of a symmetric matrix is considered and the symmetric nonnegative matrix factorization (SymNMF) is obtained by using a penalized nonsymmetric minimization problem. Instead of letting the penalizing parameter increase according to an a priori fixed rule, as suggested in literature, we propose a heuristic approach based on an adaptive technique. Extensive experimentation shows that the proposed algorithm is effective.


title: Automatic Handgun Detection in X-ray Images using Bag of Words Model with Selective Search

abstract: Baggage inspection systems using X-ray screening are crucial for security. Only 90% of threat objects are recognized from the X-ray system based in human inspection. Manual detection requires high concentration due to the images complexity and the challenges objects points of view. An algorithm based on Bag of Visual Word (BoVW) with Selective Search is proposed in this paper for handguns detection in single energy X-ray images from the public GDXray database. This approach is an adaptation of BoVW for X-ray baggage images context. In order to evaluate the proposed method the algorithm effectiveness recognition was tested on all bounding boxes returned by selective search algorithm in 200 images. The most relevant result is the precision and true positive rate (PPV = 80%, TPR= 92%). This approach achieves good performance for handgun recognition. In addition, it is the first time the Selective Search localization algorithm was tested in baggage X-ray images and showed possibilities with Bag of Visual Words.


title: Reply to Glos et al. [arXiv:1801.01294] on QSWalk.m package performance

abstract: We address comments in a recent paper by Glos et al.\ [CPC 235 (2018) 414, arXiv:1801.01294] regarding a slowdown in the performance of our QSWalk.m package for large matrix sizes. We show that the underlying issue has been rectified in the most recent version of Mathematica. We recommend all users of QSWalk.m to upgrade to version 11.3 or above.


title: Uniform Length Dominating Sequence Graphs

abstract: A sequence of vertices $(v_1,\, \dots , \,v_k)$ of a graph $G$ is called a {\it dominating closed neighborhood sequence} if $\{v_1,\, \dots , \,v_k\}$ is a dominating set of $G$ and $N[v_i]\nsubseteq \cup _{j=1}^{i-1} N[v_j]$ for every $i$. A graph $G$ is said to be {\it $k-$uniform} if all dominating closed neighborhood sequences have equal length $k$. Bre{\v s}ar et al. (2014) characterized $k$-uniform graphs with $k\leq 3$. In this article we extend their work by giving a complete characterization of all $k$-uniform graphs with $k\geq 4$.


title: Noncontrast free-breathing respiratory self-navigated coronary artery cardiovascular magnetic resonance angiography at 3 T using lipid insensitive binomial off-resonant excitation (LIBRE)

abstract: Robust and homogeneous lipid suppression is mandatory for coronary magnetic resonance angiography (MRA) since coronary arteries are commonly embedded in fat. However, effective large volume lipid suppression becomes challenging when performing radial whole-heart coronary MRA and the problem may even be exacerbated at increasing magnetic field strengths. Incomplete fat suppression also generates artifacts, and may affect advanced motion correction methods. The aim was to evaluate a recently reported lipid insensitive MRI method for self-navigated coronary MRA at 3T. Lipid insensitive binomial off resonant excitation (LIBRE) radiofrequency (RF) excitation pulses were included into a self-navigated 3D radial GRE coronary MRA sequence at 3T. LIBRE was compared against conventional fat saturation (FS) and binomial 1-180{\deg}-1 water excitation (WE). First, fat suppression of all techniques was numerically characterized using Matlab and experimentally validated in phantoms and in legs of human volunteers. Subsequently, free-breathing self-navigated coronary MRA was performed using the LIBRE pulse as well as FS and WE in ten volunteers. Results obtained in the simulations were confirmed by the experimental validations as LIBRE enabled near complete fat suppression for 3D radial imaging in vitro and in vivo. For self-navigated whole-heart coronary MRA at 3T, fat SNR was significantly attenuated using LIBRE compared with conventional FS. LIBRE increased the RCA vessel sharpness significantly (37 +/- 9% (LIBRE) vs. 29 +/- 8% (FS) and 30 +/- 8% (WE), both p<0.05) and led to a significant increase in the measured RCA vessel length to (83 +/- 31 mm (LIBRE) vs. 56 +/- 12 mm (FS) and 59 +/- 27 (WE) p<0.05). LIBRE enables robust large volume fat suppression and significantly improves coronary artery image quality at 3T compared to the use of conventional fat suppression and water excitation.


title: New lower bounds for the energy of matrices and graphs

abstract: Let $R$ be a Hermitian matrix. The energy of $R$, $\mathcal{E}(R)$, corresponds to the sum of the absolute values of its eigenvalues. In this work it is obtained two lower bounds for $\mathcal{E}(R).$ The first one generalizes a lower bound obtained by Mc Clellands for the energy of graphs in $1971$ to the case of Hermitian matrices and graphs with a given nullity. The second one generalizes a lower bound obtained by K. Das, S. A. Mojallal and I. Gutman in 2013 to symmetric non-negative matrices and graphs with a given nullity. The equality cases are discussed. These lower bounds are obtained for graphs with $m$ edges and some examples are provided showing that, some obtained bounds are incomparable with the known lower bound for the energy $2\sqrt{m}$. Another family of lower bounds are obtained from an increasing sequence of lower bounds for the spectral radius of a graph. The bounds are stated for singular and non-singular graphs.


title: The cyclic sieving phenomenon on circular Dyck paths

abstract: We give a $q$-enumeration of circular Dyck paths, which is a superset of the classical Dyck paths enumerated by the Catalan numbers. These objects have recently been studied by Alexandersson and Panova. Furthermore, we show that this $q$-analogue exhibits the cyclic sieving phenomenon under a natural action of the cyclic group. The enumeration and cyclic sieving is generalized to M\"obius paths. We also discuss properties of a generalization of cyclic sieving, which we call subset cyclic sieving. Finally, we also introduce the notion of Lyndon-like cyclic sieving that concerns special recursive properties of combinatorial objects exhibiting the cyclic sieving phenomenon.


title: Persistent Quasi-Periodic Pulsations During a Large X-Class Solar Flare

abstract: Solar flares often display pulsating and oscillatory signatures in the emission, known as quasi-periodic pulsations (QPP). QPP are typically identified during the impulsive phase of flares, yet in some cases, their presence is detected late into the decay phase. Here, we report extensive fine structure QPP that are detected throughout the large X8.2 flare from 2017 September 10. Following the analysis of the thermal pulsations observed in the GOES/XRS and the 131 A channel of SDO/AIA, we find a pulsation period of ~65 s during the impulsive phase followed by lower amplitude QPP with a period of ~150 s in the decay phase, up to three hours after the peak of the flare. We find that during the time of the impulsive QPP, the soft X-ray source observed with RHESSI rapidly rises at a velocity of approximately 17 km/s following the plasmoid/coronal mass ejection (CME) eruption. We interpret these QPP in terms of a manifestation of the reconnection dynamics in the eruptive event. During the long-duration decay phase lasting several hours, extended downward contractions of collapsing loops/plasmoids that reach the top of the flare arcade are observed in EUV. We note that the existence of persistent QPP into the decay phase of this flare are most likely related to these features. The QPP during this phase are discussed in terms of MHD wave modes triggered in the post-flaring loops.


title: Kardar-Parisi-Zhang physics in the quantum Heisenberg magnet

abstract: Equilibrium spatio-temporal correlation functions are central to understanding weak nonequilibrium physics. In certain local one-dimensional classical systems with three conservation laws they show universal features. Namely, fluctuations around ballistically propagating sound modes can be described by the celebrated Kardar-Parisi-Zhang (KPZ) universality class. Can such universality class be found also in quantum systems? By unambiguously demonstrating that the KPZ scaling function describes magnetization dynamics in the SU(2) symmetric Heisenberg spin chain we show, for the first time, that this is so. We achieve that by introducing new theoretical and numerical tools, and make a puzzling observation that the conservation of energy does not seem to matter for the KPZ physics.


title: Joint segmentation and classification of retinal arteries/veins from fundus images

abstract: Objective Automatic artery/vein (A/V) segmentation from fundus images is required to track blood vessel changes occurring with many pathologies including retinopathy and cardiovascular pathologies. One of the clinical measures that quantifies vessel changes is the arterio-venous ratio (AVR) which represents the ratio between artery and vein diameters. This measure significantly depends on the accuracy of vessel segmentation and classification into arteries and veins. This paper proposes a fast, novel method for semantic A/V segmentation combining deep learning and graph propagation. Methods A convolutional neural network (CNN) is proposed to jointly segment and classify vessels into arteries and veins. The initial CNN labeling is propagated through a graph representation of the retinal vasculature, whose nodes are defined as the vessel branches and edges are weighted by the cost of linking pairs of branches. To efficiently propagate the labels, the graph is simplified into its minimum spanning tree. Results The method achieves an accuracy of 94.8% for vessels segmentation. The A/V classification achieves a specificity of 92.9% with a sensitivity of 93.7% on the CT-DRIVE database compared to the state-of-the-art-specificity and sensitivity, both of 91.7%. Conclusion The results show that our method outperforms the leading previous works on a public dataset for A/V classification and is by far the fastest. Significance The proposed global AVR calculated on the whole fundus image using our automatic A/V segmentation method can better track vessel changes associated to diabetic retinopathy than the standard local AVR calculated only around the optic disc.


title: Estimation of the heat conducted by a cluster of small cavities and characterization of the equivalent heat conduction

abstract: We estimate the heat conducted by a cluster of many small cavities. We show that the dominating heat is a sum, over the number of the cavities, of the heats generated by each cavity after interacting with each other. This interaction is described through densities computable as solutions of a close, and invertible, system of time domain integral equations of a second kind. As an application of these expansions, we derive the effective heat conductivity which generates approximately the same heat as the cluster of cavities, distributed in a 3D bounded domain, with explicit error estimates in terms of that cluster. At the analysis level, we use time domain integral equations. Doing that, we have two choices. First, we can favor the space variable by reducing the heat potentials to the ones related to the Laplace operator (avoiding Laplace transform). Second, we can favor the time variable by reducing the representation to the Abel integral operator. As the model under investigation has time-independent parameters, we follow here the first approach.


title: Time-Dependent Surveillance-Evasion Games

abstract: Surveillance-Evasion (SE) games form an important class of adversarial trajectory-planning problems. We consider time-dependent SE games, in which an Evader is trying to reach its target while minimizing the cumulative exposure to a moving enemy Observer. That Observer is simultaneously aiming to maximize the same exposure by choosing how often to use each of its predefined patrol trajectories. Following the framework introduced in Gilles and Vladimirsky (arXiv:1812.10620), we develop efficient algorithms for finding Nash Equilibrium policies for both players by blending techniques from semi-infinite game theory, convex optimization, and multi-objective dynamic programming on continuous planning spaces. We illustrate our method on several examples with Observers using omnidirectional and angle-restricted sensors on a domain with occluding obstacles.


title: A Serious Game for Introducing Software Engineering Ethics to University Students

abstract: This paper presents a game based on storytelling, in which the players are faced with ethical dilemmas related to software engineering specific issues. The players' choices have consequences on how the story unfolds and could lead to various alternative endings. This Ethics Game was used as a tool to mediate the learning activity and it was evaluated by 144 students during a Software Engineering Course on the 2017-2018 academic year. This evaluation was based on a within-subject pre-post design methodology and provided insights on the students learning gain (academic performance), as well as on the students' perceived educational experience. In addition, it provided the results of the students' usability evaluation of the Ethics Game. The results indicated that the students did improve their knowledge about software engineering ethics by playing this game. Also, they considered this game to be a useful educational tool and of high usability. Female students had statistically significant higher knowledge gain and higher evaluation scores than male students, while no statistically significant differences were measured in groups based on the year of study.


title: Quantitative Robustness of Localized Support Vector Machines

abstract: The huge amount of available data nowadays is a challenge for kernel-based machine learning algorithms like SVMs with respect to runtime and storage capacities. Local approaches might help to relieve these issues and to improve statistical accuracy. It has already been shown that these local approaches are consistent and robust in a basic sense. This article refines the analysis of robustness properties towards the so-called influence function which expresses the differentiability of the learning method: We show that there is a differentiable dependency of our locally learned predictor on the underlying distribution. The assumptions of the proven theorems can be verified without knowing anything about this distribution. This makes the results interesting also from an applied point of view.


title: A boron-coated CCD camera for direct detection of Ultracold Neutrons (UCN)

abstract: A new boron-coated CCD camera is described for direct detection of ultracold neutrons (UCN) through the capture reactions $^{10}$B (n,$\alpha$0$\gamma$)$^7$Li (6%) and $^{10}$B(n,$\alpha$1$\gamma$)$^7$Li (94%). The experiments, which extend earlier works using a boron-coated ZnS:Ag scintillator, are based on direct detections of the neutron-capture byproducts in silicon. The high position resolution, energy resolution and particle ID performance of a scientific CCD allows for observation and identification of all the byproducts $\alpha$, $^7$Li and $\gamma$ (electron recoils). A signal-to-noise improvement on the order of 10$^4$ over the indirect method has been achieved. Sub-pixel position resolution of a few microns is demonstrated. The technology can also be used to build UCN detectors with an area on the order of 1 m$^2$. The combination of micrometer scale spatial resolution, few electrons ionization thresholds and large area paves the way to new research avenues including quantum physics of UCN and high-resolution neutron imaging and spectroscopy.


title: An Application of Jackson's $(p, q)$-Derivative to a Subclass of Starlike Functions with Negative Coefficients

abstract: In this paper, we introduce and investigate the subclass $\mathcal{P}_{p,q}^{\xi ,\kappa}(\tau, \eta)$ of starlike functions with negative coefficients by using the differential operator $\Upsilon_{\tau ,p,q}^{\xi ,\kappa}$. Coefficient inequalities, growth and distortion theorems, closure theorems, and some properties of several functions belonging to this class are obtained. We also determine the radii of close-to-convexity, starlikeness, and convexity for functions belonging to the class $\mathcal{P}_{p,q}^{\xi ,\kappa}(\tau, \eta)$. Furthermore, we obtain the integral means inequality and neighborhood results for functions belonging to the class $\mathcal{P}_{p,q}^{\xi ,\kappa}(\tau, \eta)$. The results presented in this paper generalize or improve those in related works of several earlier authors.


title: Role of geometrical cues in neuronal growth

abstract: Geometrical cues play an essential role in neuronal growth. Here, we quantify axonal growth on surfaces with controlled geometries and report a general stochastic approach that quantitatively describes the motion of growth cones. We show that axons display a strong directional alignment on micro-patterned surfaces when the periodicity of the patterns matches the dimension of the growth cone. The growth cone dynamics on surfaces with uniform geometry is described by a linear Langevin equation with both deterministic and stochastic contributions. In contrast, axonal growth on surfaces with periodic patterns is characterized by a system of two generalized Langevin equations with both linear and quadratic velocity dependence and stochastic noise. We combine experimental data with theoretical analysis to measure the key parameters of the growth cone motion: angular distributions, correlation functions, diffusion coefficients, characteristics speeds and damping coefficients. We demonstrate that axonal dynamics displays a cross-over from an Ornstein-Uhlenbeck process to a non-linear stochastic regime when the geometrical periodicity of the pattern approaches the linear dimension of the growth cone. Growth alignment is determined by surface geometry, which is fully quantified by the deterministic part of the Langevin equation. These results provide new insight into the role of curvature sensing proteins and their interactions with geometrical cues.


title: An application of the Gy\'{a}rf\'{a}s path argument

abstract: We adapt the Gy\'{a}rf\'{a}s path argument to prove that $t-2$ cops can capture a robber, in at most $t-1$ moves, in the game of cops and robbers played in a graph that does not contain the $t$-vertex path as an induced subgraph.


title: A solid-state entangled photon pair source with high brightness and indistinguishability

abstract: The generation of high-quality entangled photon pairs has been being a long-sought goal in modern quantum communication and computation. To date, the most widely-used entangled photon pairs are generated from spontaneous parametric downconversion, a process that is intrinsically probabilistic and thus relegated to a regime of low pair-generation rates. In contrast, semiconductor quantum dots can generate triggered entangled photon pairs via a cascaded radiative decay process, and do not suffer from any fundamental trade-off between source brightness and multi-pair generation. However, a source featuring simultaneously high photon-extraction efficiency, high-degree of entanglement fidelity and photon indistinguishability has not yet been reported. Here, we present an entangled photon pair source with high brightness and indistinguishability by deterministically embedding GaAs quantum dots in broadband photonic nanostructures that enable Purcell-enhanced emission. Our source produces entangled photon pairs with a record pair collection probability of up to 0.65(4) (single-photon extraction efficiency of 0.85(3)), entanglement fidelity of 0.88(2), and indistinguishabilities of 0.901(3) and 0.903(3), which immediately creates opportunities for advancing quantum photonic technologies.


title: Beam Squint and Channel Estimation for Wideband mmWave Massive MIMO-OFDM Systems

abstract: With the increasing scale of antenna arrays in wideband millimeter-wave (mmWave) communications, the physical propagation delays of electromagnetic waves traveling across the whole array will become large and comparable to the time-domain sample period, which is known as the spatial-wideband effect. In this case, different subcarriers in an orthogonal frequency division multiplexing (OFDM) system will "see" distinct angles of arrival (AoAs) for the same path. This effect is known as beam squint, resulting from the spatial-wideband effect, and makes the approaches based on the conventional multiple-input multiple-output (MIMO) model, such as channel estimation and precoding, inapplicable. After discussing the relationship between beam squint and the spatial-wideband effect, we propose a channel estimation scheme for frequency-division duplex (FDD) mmWave massive MIMO-OFDM systems with hybrid analog/digital precoding, which takes the beam squint effect into consideration. A super-resolution compressed sensing approach is developed to extract the frequency-insensitive parameters of each uplink channel path, i.e., the AoA and the time delay, and the frequency-sensitive parameter, i.e., the complex channel gain. With the help of the reciprocity of these frequency-insensitive parameters in FDD systems, the downlink channel estimation can be greatly simplified, where only limited pilots are needed to obtain downlink complex gains and reconstruct downlink channels. Furthermore, the uplink and downlink channel covariance matrices can be constructed from these frequency-insensitive channel parameters rather than through a long-term average, which enables the minimum mean-squared error (MMSE) channel estimation to further enhance performance. Numerical results demonstrate the superiority of the proposed scheme over the conventional methods in mmWave communications.


title: Using stigmergy as a computational memory in the design of recurrent neural networks

abstract: In this paper, a novel architecture of Recurrent Neural Network (RNN) is designed and experimented. The proposed RNN adopts a computational memory based on the concept of stigmergy. The basic principle of a Stigmergic Memory (SM) is that the activity of deposit/removal of a quantity in the SM stimulates the next activities of deposit/removal. Accordingly, subsequent SM activities tend to reinforce/weaken each other, generating a coherent coordination between the SM activities and the input temporal stimulus. We show that, in a problem of supervised classification, the SM encodes the temporal input in an emergent representational model, by coordinating the deposit, removal and classification activities. This study lays down a basic framework for the derivation of a SM-RNN. A formal ontology of SM is discussed, and the SM-RNN architecture is detailed. To appreciate the computational power of an SM-RNN, comparative NNs have been selected and trained to solve the MNIST handwritten digits recognition benchmark in its two variants: spatial (sequences of bitmap rows) and temporal (sequences of pen strokes).


title: Random Walks on Dynamic Graphs: Mixing Times, HittingTimes, and Return Probabilities

abstract: We establish and generalise several bounds for various random walk quantities including the mixing time and the maximum hitting time. Unlike previous analyses, our derivations are based on rather intuitive notions of local expansion properties which allows us to capture the progress the random walk makes through $t$-step probabilities. We apply our framework to dynamically changing graphs, where the set of vertices is fixed while the set of edges changes in each round. For random walks on dynamic connected graphs for which the stationary distribution does not change over time, we show that their behaviour is in a certain sense similar to static graphs. For example, we show that the mixing and hitting times of any sequence of $d$-regular connected graphs is $O(n^2)$, generalising a well-known result for static graphs. We also provide refined bounds depending on the isoperimetric dimension of the graph, matching again known results for static graphs. Finally, we investigate properties of random walks on dynamic graphs that are not always connected: we relate their convergence to stationarity to the spectral properties of an average of transition matrices and provide some examples that demonstrate strong discrepancies between static and dynamic graphs.


title: Dynamics of a yoyo using a smartphone gyroscope sensor

abstract: The dynamics of a traditional toy, the yoyo, is investigated theoretically and experimentally using smartphone' sensors. In particular, using the gyroscope the angular velocity is measured. The experimental results are complemented thanks to a digital video analysis. The concordance between theoretical and experimental results is discussed. As the yoyo is a ubiquitous, simple and traditional toy this simple proposal could encourage students to experiment with everyday objects and modern technologies.


title: Hybrid Actor-Critic Reinforcement Learning in Parameterized Action Space

abstract: In this paper we propose a hybrid architecture of actor-critic algorithms for reinforcement learning in parameterized action space, which consists of multiple parallel sub-actor networks to decompose the structured action space into simpler action spaces along with a critic network to guide the training of all sub-actor networks. While this paper is mainly focused on parameterized action space, the proposed architecture, which we call hybrid actor-critic, can be extended for more general action spaces which has a hierarchical structure. We present an instance of the hybrid actor-critic architecture based on proximal policy optimization (PPO), which we refer to as hybrid proximal policy optimization (H-PPO). Our experiments test H-PPO on a collection of tasks with parameterized action space, where H-PPO demonstrates superior performance over previous methods of parameterized action reinforcement learning.


title: Teaching HCI Design in a Flipped Learning M.Sc. Course Using Eye-Tracking Peer Evaluation Data

abstract: This paper presents experiences from a flipped classroom M.Sc. course on Human-Computer Interaction (HCI). The students that finished successfully this course participated in twelve short workshops, based on a flipped classroom model. Each workshop focused on a specific HCI activity, while before the workshops, a two-hour lecture was used to introduce the students in the flipped learning concept. All the rest of the educational material was offered to the students online before each workshop. Such material was mainly short lectures from the professor, in the form of videos uploaded in the course's YouTube channel and documents delivered using the university learning management system (LMS). For each workshop the students had to be prepared to participate, which was tested using brief quizzes before the start of specific workshops. The activity presented in this paper was the design and evaluation of an interactive system. The students were asked to form six groups comprising of three to four students each. Then a system's description, vague enough to stimulate creativity, was randomly assigned to each group. This activity presented in this paper was the longest activity of the entire course and it was conducted in four consequent workshops. The paper presents the setting of this experiment, the peer assessment method and the use of eye-tracking data collected and analysed to aid the students towards improving their design. The students created a working model of the system with limited functionality and improved this model using eye-tracking data from the peer evaluation of this model. The use of these data offered them the insight to improve their models and to undergo design changes. The paper presents samples of the progress made between various versions of the models and concludes presenting the preliminary positive results of the students qualitative evaluation of this experiment.


title: Structure Constants from Modularity in Warped CFT

abstract: We derive a universal formula for the asymptotic growth of the mean value of three-point coefficient for Warped Conformal Field Theories (WCFTs), and provide a holographic calculation in BTZ black holes. WCFTs are two dimensional quantum field theories featuring a chiral Virasoro and U(1) Kac-Moody algebra, and are conjectured to be holographically dual to quantum gravity on asymptotically AdS$_3$ spacetime with Comp$\grave{\mathrm{e}}$re-Song-Strominger boundary conditions. The WCFT calculation amounts to the calculation of one-point functions on torus, whose high temperature limit can be approximated by using modular covariance of WCFT, similar to the derivation of Cardy formula. The bulk process is given by a tadpole diagram, with a massive spinning particle propagates from the infinity to the horizon, and splits into particle and antiparticle which annihilate after going around the horizon of BTZ black holes. The agreement between the bulk and WCFT calculations indicates that the black hole geometries in asymptotically AdS$_3$ spacetimes can emerge upon coarse-graining over microstates in WCFTs, similar to the results of Kraus and Maloney in the context of AdS/CFT[arXiv:1608.03284].


title: Reduced Focal Loss: 1st Place Solution to xView object detection in Satellite Imagery

abstract: This paper describes our approach to the DIUx xView 2018 Detection Challenge [1]. This challenge focuses on a new satellite imagery dataset. The dataset contains 60 object classes that are highly imbalanced. Due to the imbalanced nature of the dataset, the training process becomes significantly more challenging. To address this problem, we introduce a novel Reduced Focal Loss function, which brought us 1st place in the DIUx xView 2018 Detection Challenge.


title: Systematics of toroidal dipole modes in Ca, Ni, Zr, and Sn isotopes

abstract: We analyze the relation between isoscalar toroidal modes and so-called pygmy dipole resonance (PDR) which both appear in the same region of low-energy dipole excitations. To this end, we use a theoretical description within the fully self-consistent Skyrme quasiparticle random-phase approximation (QRPA). Test cases are spherical nuclei $^{40,48}$Ca, $^{58,72}$Ni, $^{90,100}$Zr, and $^{100,120,132}$Sn which cover four different elements and for each element at least two isotopes with different neutron excess, one small and another large. The structure of the modes is investigated in terms of strength functions, transition densities (TD) and transition currents (TC). For all considered nuclei, we see that, independently on whether PDR strength exists or not, the flow pattern in the lower part of the "PDR energy region" is basically the isoscalar vortical toroidal motion with a minor irrotational fraction. A one-to-one correspondence between calculated TD and TC can be established. The toroidal flow appears already in the uncoupled two-quasiparticle (2qp) excitations and becomes definitively strong for the QRPA modes. Altogether, we find that low-lying dipole strength often denoted as isoscalar PDR is actually an oversimplified imitation of the basically toroidal motion in nuclei with a sufficient neutron excess.


title: Knowledge of Quantum Hidden Variables Enables Backwards-In-Time Signaling

abstract: Bell's theorem implies that any completion of quantum mechanics which uses hidden variables (that is, preexisting values of all observables) must be nonlocal in the Einstein sense. This customarily indicates that knowledge of the hidden variables would permit superluminal communication. Such superluminal signaling, akin to the existence of a preferred reference frame, is to be expected. However, here we provide a protocol that allows an observer with knowledge of the hidden variables to communicate with her own causal past, without superluminal signaling. That is, such knowledge would contradict causality, irrespectively of the validity of relativity theory. Among the ways we propose for bypassing the paradox there is the possibility of hidden variables that change their values even when the state does not, and that means that signaling backwards in time is prohibited in Bohmian mechanics.


title: Toward Achieving Formal Guarantees for Human-Aware Controllers in Human-Robot Interactions

abstract: With the primary objective of human-robot interaction being to support humans' goals, there exists a need to formally synthesize robot controllers that can provide the desired service. Synthesis techniques have the benefit of providing formal guarantees for specification satisfaction. There is potential to apply these techniques for devising robot controllers whose specifications are coupled with human needs. This paper explores the use of formal methods to construct human-aware robot controllers to support the productivity requirements of humans. We tackle these types of scenarios via human workload-informed models and reactive synthesis. This strategy allows us to synthesize controllers that fulfill formal specifications that are expressed as linear temporal logic formulas. We present a case study in which we reason about a work delivery and pickup task such that the robot increases worker productivity, but not stress induced by high work backlog. We demonstrate our controller using the Toyota HSR, a mobile manipulator robot. The results demonstrate the realization of a robust robot controller that is guaranteed to properly reason and react in collaborative tasks with human partners.


title: Homological Berglund-H\"ubsch mirror symmetry for curve singularities

abstract: Given a two-variable invertible polynomial, we show that its category of maximally-graded matrix factorisations is quasi-equivalent to the Fukaya-Seidel category of its Berglund-H\"ubsch transpose. This was previously shown for Brieskorn-Pham and $D$-type singularities by Futaki-Ueda. The proof involves explicit construction of a tilting object on the B-side, and comparison with a specific basis of Lefschetz thimbles on the A-side.


title: Learning Sensory-Motor Associations from Demonstration

abstract: We propose a method which generates reactive robot behavior learned from human demonstration. In order to do so, we use the Playful programming language which is based on the reactive programming paradigm. This allows us to represent the learned behavior as a set of associations between sensor and motor primitives in a human readable script. Distinguishing between sensor and motor primitives introduces a supplementary level of granularity and more importantly enforces feedback, increasing adaptability and robustness. As the experimental section shows, useful behaviors may be learned from a single demonstration covering a very limited portion of the task space.


title: The nature of GRB host galaxies from chemical abundances

abstract: We try to identify the nature of high redshift long Gamma-Ray Bursts (LGRBs) host galaxies by comparing the observed abundance ratios in the interstellar medium with detailed chemical evolution models accounting for the presence of dust. We compared measured abundance data from LGRB afterglow spectra to abundance patterns as predicted by our models for different galaxy types. We analysed in particular [X/Fe] abundance ratios (where X is C, N, O, Mg, Si, S, Ni, Zn) as functions of [Fe/H]. Different galaxies (irregulars, spirals, ellipticals) are, in fact, characterised by different star formation histories, which produce different [X/Fe] ratios ("time-delay model"). This allows us to identify the morphology of the hosts and to infer their age (i.e. the time elapsed from the beginning of star formation) at the time of the GRB events, as well as other important parameters. Relative to previous works, we use newer models in which we adopt updated stellar yields and prescriptions for dust production, accretion and destruction. We have considered a sample of seven LGRB host galaxies. Our results have suggested that two of them (GRB 050820, GRB 120815A) are ellipticals, two (GRB 081008, GRB 161023A) are spirals and three (GRB 050730, GRB 090926A, GRB 120327A) are irregulars. We also found that in some cases changing the initial mass function can give better agreement with the observed data. The calculated ages of the host galaxies span from the order of 10 Myr to little more than 1 Gyr.


title: Spectral analysis and parameter estimation in levitated optomechanics

abstract: Optical levitation of nanoscale particles has emerged as a platform for precision measurement. Extremely low damping, together with optical interferometric position detection, makes possible exquisite force measurement and promises low-energy tests of fundamental physics. Essential to such measurement is an understanding of the confidence with which parameters can be inferred from spectra estimated from the indirect measurement provided by interferometry. We present an apparatus optimized for sensitivity along one motional degree of freedom, a theoretical model of the spectrum, and maximum likelihood estimation. The treatment accounts for the sinusoidal dependence of interferometric signal on particle position, and we use the technique to extract thermodynamic quantities in a regime where simpler treatments are confounded.


title: Secure Data Offloading Strategy for Connected and Autonomous Vehicles

abstract: Connected and Automated Vehicles (CAVs) are expected to constantly interact with a network of processing nodes installed in secure cabinets located at the side of the road -- thus, forming Fog Computing-based infrastructure for Intelligent Transportation Systems (ITSs). Future city-scale ITS services will heavily rely upon the sensor data regularly off-loaded by each CAV on the Fog Computing network. Due to the broadcast nature of the medium, CAVs' communications can be vulnerable to eavesdropping. This paper proposes a novel data offloading approach where the Random Linear Network Coding (RLNC) principle is used to ensure the probability of an eavesdropper to recover relevant portions of sensor data is minimized. Our preliminary results confirm the effectiveness of our approach when operated in a large-scale ITS networks.


title: A sequence of polynomials with optimal condition number

abstract: We find an explicit sequence of univariate polynomials of arbitrary degree with optimal condition number. This solves a problem posed by Michael Shub and Stephen Smale in 1993.


title: Quasioptical modeling of wave beams with and without mode conversion: II. Numerical simulations of single-mode beams

abstract: This work continues a series of papers where we propose an algorithm for quasioptical modeling of electromagnetic beams with and without mode conversion. The general theory was reported in the first paper of this series, where a parabolic partial differential equation was derived for the field envelope that may contain one or multiple modes with close group velocities. Here, we present a corresponding code PARADE (PAraxial RAy DEscription) and its test applications to single-mode beams in vacuum and also in inhomogeneous magnetized plasma. The numerical results are compared, respectively, with analytic formulas from Gaussian-beam optics and also with cold-plasma ray tracing. Quasioptical simulations of mode-converting beams are reported in the next, third paper of this series.


title: Extremal total distance of graphs of given radius

abstract: In 1984, Plesn\'{i}k determined the minimum total distance for given order and diameter and characterized the extremal graphs and digraphs. We prove the analog for given order and radius, when the order is sufficiently large compared to the radius. This confirms asymptotically a conjecture of Chen et al. We show the connection between minimizing the total distance and maximizing the size under the same conditions. We also prove some asymptotically optimal bounds for the maximum total distance.


title: Realizing Quantum Boltzmann Machines Through Eigenstate Thermalization

abstract: Quantum Boltzmann machines are natural quantum generalizations of Boltzmann machines that are expected to be more expressive than their classical counterparts, as evidenced both numerically for small systems and asymptotically under various complexity theoretic assumptions. However, training quantum Boltzmann machines using gradient-based methods requires sampling observables in quantum thermal distributions, a problem that is NP-hard. In this work, we find that the locality of the gradient observables gives rise to an efficient sampling method based on the Eigenstate Thermalization Hypothesis, and thus through Hamiltonian simulation an efficient method for training quantum Boltzmann machines on near-term quantum devices. Furthermore, under realistic assumptions on the moments of the data distribution to be modeled, the distribution sampled using our algorithm is approximately the same as that of an ideal quantum Boltzmann machine. We demonstrate numerically that under the proposed training scheme, quantum Boltzmann machines capture multimodal Bernoulli distributions better than classical restricted Boltzmann machines with the same connectivity structure. We also provide numerical results on the robustness of our training scheme with respect to noise.


title: A comprehensive real-time analysis model at the LHCb experiment

abstract: An evolved real-time data processing strategy is proposed for high-energy physics experiments, and its implementation at the LHCb experiment is presented. The reduced event model allows not only the signal candidate firing the trigger to be persisted, as previously available, but also an arbitrary set of other reconstructed or raw objects from the event. This allows for higher trigger rates for a given output data bandwidth, when compared to the traditional model of saving the full raw detector data for each trigger, whilst accommodating inclusive triggers and preserving data mining capabilities. The gains in physics reach and savings in computing resources already made possible by the model are discussed, along with the prospects of employing it more widely for Run 3 of the Large Hadron Collider.


title: Broadening and saturation effects in dijet azimuthal correlations in p-p and p-Pb collisions at $\mathbf{\sqrt{s_{NN}}} = $ 5.02 TeV

abstract: We demonstrate that the recent forward-forward dijet correlation data measured by the ATLAS collaboration for proton-proton and proton-lead collisions are consistent with the broadening effects due to both the gluon saturation and the resummation of large logarithms of the hard scale (the so-called Sudakov logarithms). We find that both effects are necessary to describe the experimental results.


title: Simulation study of estimating between-study variance and overall effect in meta-analysis of standardized mean difference

abstract: Methods for random-effects meta-analysis require an estimate of the between-study variance, $\tau^2$. The performance of estimators of $\tau^2$ (measured by bias and coverage) affects their usefulness in assessing heterogeneity of study-level effects, and also the performance of related estimators of the overall effect. For the effect measure standardized mean difference (SMD), we provide the results from extensive simulations on five point estimators of $\tau^2$ (the popular methods of DerSimonian-Laird, restricted maximum likelihood, and Mandel and Paule (MP); the less-familiar method of Jackson; the new method (KDB) based on the improved approximation to the distribution of the Q statistic by Kulinskaya, Dollinger and Bj{\o}rkest{\o}l (2011) ), five interval estimators for $\tau^2$ (profile likelihood, Q-profile, Biggerstaff and Jackson, Jackson, and the new KDB method), six point estimators of the overall effect (the five related to the point estimators of $\tau^2$ and an estimator whose weights use only study-level sample sizes), and eight interval estimators for the overall effect (five based on the point estimators for $\tau^2$; the Hartung-Knapp-Sidik-Jonkman (HKSJ) interval; a modification of HKSJ; and an interval based on the sample-size-weighted estimator).


title: Opportunistic View Materialization with Deep Reinforcement Learning

abstract: Carefully selected materialized views can greatly improve the performance of OLAP workloads. We study using deep reinforcement learning to learn adaptive view materialization and eviction policies. Our insight is that such selection policies can be effectively trained with an asynchronous RL algorithm, that runs paired counter-factual experiments during system idle times to evaluate the incremental value of persisting certain views. Such a strategy obviates the need for accurate cardinality estimation or hand-designed scoring heuristics. We focus on inner-join views and modeling effects in a main-memory, OLAP system. Our research prototype system, called DQM, is implemented in SparkSQL and we experiment on several workloads including the Join Order Benchmark and the TPC-DS workload. Results suggest that: (1) DQM can outperform heuristic when their assumptions are not satisfied by the workload or there are temporal effects like period maintenance, (2) even with the cost of learning, DQM is more adaptive to changes in the workload, and (3) DQM is broadly applicable to different workloads and skews.


title: Quasioptical modeling of wave beams with and without mode conversion: III. Numerical simulations of mode-converting beams

abstract: This work continues a series of papers where we propose an algorithm for quasioptical modeling of electromagnetic beams with and without mode conversion. The general theory was reported in the first paper of this series, where a parabolic partial differential equation was derived for the field envelope that may contain one or multiple modes with close group velocities. In the second paper, we presented a corresponding code PARADE (PAraxial RAy DEscription) and its test applications to single-mode beams. Here, we report quasioptical simulations of mode-converting beams for the first time. We also demonstrate that PARADE can model splitting of two-mode beams. The numerical results produced by PARADE show good agreement with those of one-dimensional full-wave simulations and also with conventional ray tracing (to the extent that one-dimensional and ray-tracing simulations are applicable).


title: Microscopic Traffic Simulation by Cooperative Multi-agent Deep Reinforcement Learning

abstract: Expert human drivers perform actions relying on traffic laws and their previous experience. While traffic laws are easily embedded into an artificial brain, modeling human complex behaviors which come from past experience is a more challenging task. One of these behaviors is the capability of communicating intentions and negotiating the right of way through driving actions, as when a driver is entering a crowded roundabout and observes other cars movements to guess the best time to merge in. In addition, each driver has its own unique driving style, which is conditioned by both its personal characteristics, such as age and quality of sight, and external factors, such as being late or in a bad mood. For these reasons, the interaction between different drivers is not trivial to simulate in a realistic manner. In this paper, this problem is addressed by developing a microscopic simulator using a Deep Reinforcement Learning Algorithm based on a combination of visual frames, representing the perception around the vehicle, and a vector of numerical parameters. In particular, the algorithm called Asynchronous Advantage Actor-Critic has been extended to a multi-agent scenario in which every agent needs to learn to interact with other similar agents. Moreover, the model includes a novel architecture such that the driving style of each vehicle is adjustable by tuning some of its input parameters, permitting to simulate drivers with different levels of aggressiveness and desired cruising speeds.


title: Graphical Calculus for products and convolutions

abstract: Graphical calculus is an intuitive visual notation for manipulating tensors and index contractions. Using graphical calculus leads to simple and memorable derivations, and with a bit of practice one can learn to prove complex identities even without the need for pen and paper. This manuscript is meant as a demonstration of the power and flexibility of graphical notation and we advocate exploring the use of graphical calculus in undergraduate courses. In the first part we define the following matrix products in graphical language: dot, tensor, Kronecker, Hadamard, Kathri-Rao and Tracy-Singh. We then use our definitions to prove several known identities in an entirely graphical way. Despite ordinary proofs consist in several lines of quite involved mathematical expressions, graphical calculus is so expressive that after writing an identity in graphical form we can realise by visual inspection that it is in fact true. As an example of the intuitiveness of graphical proofs, we derive two new identities. In the second part we develop a graphical description of convolutions, which is a central ingredient of convolutional neural networks and signal processing. Our single definition includes as special cases the circular discrete convolution and the cross-correlation. We illustrate how convolution can be seen as another type of product and we derive a generalised convolution theorem. We conclude with a quick guide on implementing tensor contractions in python.


title: Measuring valley polarization in two-dimensional materials with second-harmonic spectroscopy

abstract: A population imbalance at different valleys of an electronic system lowers its effective rotational symmetry. We introduce a technique to measure such imbalance - a valley polarization - that exploits the unique fingerprints of this symmetry reduction in the polarization-dependent second-harmonic generation (SHG). We present the principle and detection scheme in the context of hexagonal two-dimensional crystals, which include graphene-based systems and the family of transition metal dichalcogenides, and provide a direct experimental demonstration using a 2H-MoSe$_{2}$ monolayer at room temperature. We deliberately use the simplest possible setup, where a single pulsed laser beam simultaneously controls the valley imbalance and tracks the SHG process. We further developed a model of the transient population dynamics which analytically describes the valley-induced SHG rotation in very good agreement with the experiment. In addition to providing the first experimental demonstration of the effect, this work establishes a conceptually simple, com-pact and transferable way of measuring instantaneous valley polarization, with direct applicability in the nascent field of valleytronics.


title: Sequential Relational Decomposition

abstract: The concept of decomposition in computer science and engineering is considered a fundamental component of computational thinking and is prevalent in design of algorithms, software construction, hardware design, and more. We propose a simple and natural formalization of sequential decomposition, in which a task is decomposed into two sequential sub-tasks, with the first sub-task to be executed before the second sub-task is executed. These tasks are specified by means of input/output relations. We define and study decomposition problems, which is to decide whether a given specification can be sequentially decomposed. Our main result is that decomposition itself is a difficult computational problem. More specifically, we study decomposition problems in three settings: where the input task is specified explicitly, by means of Boolean circuits, and by means of automatic relations. We show that in the first setting decomposition is NP-complete, in the second setting it is NEXPTIME-complete, and in the third setting there is evidence to suggest that it is undecidable. Our results indicate that the intuitive idea of decomposition as a system-design approach requires further investigation. In particular, we show that adding a human to the loop by asking for a decomposition hint lowers the complexity of decomposition problems considerably.


title: Eliminating Interference in LOS Massive Multi-User MIMO with a Few Transceivers

abstract: Wireless cellular communication networks are bandwidth and interference limited. An important means to overcome these resource limitations is the use of multiple antennas. Base stations equipped with a very large (massive) number of antennas have been the focus of recent research. A bottleneck in such systems is the limited number of transmit/receive chains. In this work, a line-of-sight (LOS) channel model is considered. It is shown that for a given number of interferers, it suffices that the number of transmit/receive chains exceeds the number of desired users by one, assuming a sufficiently large antenna array. From a theoretical point of view, this is the first result proving the near-optimal performance of antenna selection, even when the total number of signals (desired and interfering) is larger than the number of receive chains. Specifically, a single additional chain suffices to reduce the interference to any desired level. We prove that using the proposed selection, a simple linear receiver/transmitter for the uplink/downlink provides near-optimal rates. In particular, in the downlink direction, there is no need for complicated dirty paper coding; each user can use an optimal code for a single user interference-free channel. In the uplink direction, there is almost no gain in implementing joint decoding. The proposed approach is also a significant improvement both from system and computational perspectives. Simulation results demonstrating the performance of the proposed method are provided.


title: Identification and Validation of Virtual Battery Model for Heterogeneous Devices

abstract: The potential of distributed energy resources in providing grid services can be maximized with the recent advancements in demand side control. Effective utilization of this control strategy requires the knowledge of aggregate flexibility of the distributed energy resources (DERs). Recent works have shown that the aggregate flexibility of DERs can be modeled as a virtual battery (VB) whose state evolution is governed by a first order system including self-dissipation. The VB parameters (self-dissipation rate, energy capacity) are obtained by solving an optimization problem which minimizes the tracking performance of the ensemble and the proposed first-order model. For the identified first order model, time-varying power limits are calculated using binary search algorithms. Finally, this proposed framework is demonstrated for different homogeneous and heterogeneous ensembles consisting of air conditioners (ACs) and electric water heaters (EWHs).


title: Some corrections to an old paper

abstract: We give some corrections of our paper "Primes in arithmetic progressions to large moduli"[BFI]. The corrections do not affect the statements of any of the theorems in the paper. The contents of our two sequel papers [BFI2, BFI3] also remain unchanged.


title: Efficient Millimeter-Wave Infrastructure Placement for City-Scale ITS

abstract: Millimeter Waves (mmWaves) will play a pivotal role in the next-generation of Intelligent Transportation Systems (ITSs). However, in deep urban environments, sensitivity to blockages creates the need for more sophisticated network planning. In this paper, we present an agile strategy for deploying road-side nodes in a dense city scenario. In our system model, we consider strict Quality-of-Service (QoS) constraints (e.g. high throughput, low latency) that are typical of ITS applications. Our approach is scalable, insofar that takes into account the unique road and building shapes of each city, performing well for both regular and irregular city layouts. It allows us not only to achieve the required QoS constraints but it also provides up to $50\%$ reduction in the number of nodes required, compared to existing deployment solutions.


title: $\alpha$-Rank: Multi-Agent Evaluation by Evolution

abstract: We introduce $\alpha$-Rank, a principled evolutionary dynamics methodology for the evaluation and ranking of agents in large-scale multi-agent interactions, grounded in a novel dynamical game-theoretic solution concept called Markov-Conley chains (MCCs). The approach leverages continuous- and discrete-time evolutionary dynamical systems applied to empirical games, and scales tractably in the number of agents, the type of interactions, and the type of empirical games (symmetric and asymmetric). Current models are fundamentally limited in one or more of these dimensions and are not guaranteed to converge to the desired game-theoretic solution concept (typically the Nash equilibrium). $\alpha$-Rank provides a ranking over the set of agents under evaluation and provides insights into their strengths, weaknesses, and long-term dynamics. This is a consequence of the links we establish to the MCC solution concept when the underlying evolutionary model's ranking-intensity parameter, $\alpha$, is chosen to be large, which exactly forms the basis of $\alpha$-Rank. In contrast to the Nash equilibrium, which is a static concept based on fixed points, MCCs are a dynamical solution concept based on the Markov chain formalism, Conley's Fundamental Theorem of Dynamical Systems, and the core ingredients of dynamical systems: fixed points, recurrent sets, periodic orbits, and limit cycles. $\alpha$-Rank runs in polynomial time with respect to the total number of pure strategy profiles, whereas computing a Nash equilibrium for a general-sum game is known to be intractable. We introduce proofs that not only provide a unifying perspective of existing continuous- and discrete-time evolutionary evaluation models, but also reveal the formal underpinnings of the $\alpha$-Rank methodology. We empirically validate the method in several domains including AlphaGo, AlphaZero, MuJoCo Soccer, and Poker.