title: Interplay of Game Incentives, Player Profiles and Task Difficulty in Games with a Purpose

abstract: How to take multiple factors into account when evaluating a Game with a Purpose? How is player behaviour or participation influenced by different incentives? How does player engagement impact their accuracy in solving tasks? In this paper, we present a detailed investigation of multiple factors affecting the evaluation of a GWAP and we show how they impact on the achieved results. We inform our study with the experimental assessment of a GWAP designed to solve a multinomial classification task.


title: A novel gravitational lensing feature by wormholes

abstract: Horizonless compact objects with light rings (or photon spheres) are becoming increasingly popular in recent years for several reasons. In this paper, we show that a horizonless object such as a wormhole of Morris-Thorne type can have two photon spheres. In particular, we show that, in addition to the one present outside a wormhole throat, the throat can itself act as an effective photon sphere. Such wormholes exhibit two sets of relativistic Einstein ring systems formed due to strong gravitational lensing. We consider a previously obtained wormhole solution as a specific example. If such type of wormhole casts a shadow at all, then the inner set of the relativistic Einstein rings will form the outer bright edge of the shadow. Such a novel lensing feature might serve as a distinguishing feature between wormholes and other celestial objects as far as gravitational lensing is concerned.


title: On the Riemann hypothesis for self-dual weight enumerators of genus three

abstract: In this note, we give an equivalent condition for a self-dual weight enumerator of genus three to satisfy the Riemann hypothesis. We also observe the truth and falsehood of the Riemann hypothesis for some families of invariant polynomials.


title: Analysis of a moving mask approximation for martensitic transformations

abstract: In this work we introduce a moving mask approximation to describe the dynamics of austenite to martensite phase transitions at a continuum level. In this framework, we prove a new type of Hadamard jump condition, from which we deduce that the deformation gradient must be of the form $\mathsf{1} +\mathbf{a}\otimes \mathbf{n}$ a.e. in the martensite phase. This is useful to better understand the complex microstructures and the formation of curved interfaces between phases in new ultra-low hysteresis alloys such as Zn45Au30Cu25, and provides a selection mechanism for physically-relevant energy-minimising microstructures. In particular, we use the new type of Hadamard jump condition to deduce a rigidity theorem for the two well problem. The latter provides more insight on the cofactor conditions, particular conditions of supercompatibility between phases believed to influence reversibility of martensitic transformations.


title: On the Asymptotics of a Prime Spin Relation II

abstract: For K a cyclic cubic number field with odd class number containing a unit w such that Norm(w)=Norm(1-w)=-1, we prove that the density of rational primes p that satisfy the given spin relation is equal to 1/2. Furthermore, we prove that this density restricted to rational primes that are 1 mod 4 is 1/4 and this density restricted to rational primes that are -1 mod 4 is 3/4.


title: Revised estimates of the frequency of Earth-like planets in the Kepler field

abstract: The search for Earth-like planets around Sun-like stars and the evaluation of their occurrence rate is a major topic of research for the exoplanetary community. Two key characteristics in defining a planet as 'Earth-like' are having a radius between 1 and 1.75 times the Earth's radius and orbiting inside the host star's habitable zone; the measurement of the planet's radius and related error is however possible only via transit observations and is highly dependent on the precision of the host star's radius. A major improvement in the determination of stellar radius is represented by the unprecedented precision on parallax measurements provided by the Gaia astrometry satellite. We present a new estimate of the frequency of Earth-sized planets orbiting inside the host stars's habitable zones, obtained using Gaia measurements of parallax for solar-type stars hosting validated planets in the Kepler field as input for reassessing the values of planetary radius and incident stellar flux. This updated occurrence rate can usefully inform future observational efforts searching for Earth-like system in the Sun backyard using a variety of techniques such as the spectrograph ESPRESSO, the space observatory PLATO and the proposed astrometric satellite Theia.


title: Subleading power corrections to $B\to \gamma l\nu $ decay in PQCD approach

abstract: The leptonic radiative decay $B \to \gamma l\nu$ is of great importance in the determination of $B$ meson wave functions, and evaluating the form factors $ F_{V,A}$ are the essential problem on the study of this channel. We computed next-to-leading power corrections to the form factors within the framework of PQCD approach, including the power suppressed hard kernel, the contribution from a complete set of three-particle $B$ meson wave functions up to twist-4 and two-particle off light-cone wave functions, the $1/m_b$ corrections in heavy quark effective theory, and the contribution from hadronic structure of photon. In spite of large theoretical uncertainties, the overall power suppressed contributions decreases about $50\%$ of the leading power result. The $\lambda_B$ dependence of the integrated branching ratio is reduced after including the subleading power contributions, thus the power corrections lead to more ambiguity in the determination of $\lambda_B$ from $B \to \gamma l\nu$ decay.


title: The maximal discrete extension of $SL_2(\mathcal{\scriptstyle{O}}_K)$ for an imaginary-quadratic number field $K$

abstract: Let $\mathcal{\scriptstyle{O}}_K$ be the ring of integers of an imaginary quadratic number field $K$. In this paper we give a new description of the maximal discrete extension of the group $SL_2(\mathcal{\scriptstyle{O}}_K)$ inside $SL_2(\mathbb{C})$, which uses generalized Atkin-Lehner involutions. Moreover we find a natural characterization of this group in $SO(1,3)$.


title: Deep Unfolded Robust PCA with Application to Clutter Suppression in Ultrasound

abstract: Contrast enhanced ultrasound is a radiation-free imaging modality which uses encapsulated gas microbubbles for improved visualization of the vascular bed deep within the tissue. It has recently been used to enable imaging with unprecedented subwavelength spatial resolution by relying on super-resolution techniques. A typical preprocessing step in super-resolution ultrasound is to separate the microbubble signal from the cluttering tissue signal. This step has a crucial impact on the final image quality. Here, we propose a new approach to clutter removal based on robust principle component analysis (PCA) and deep learning. We begin by modeling the acquired contrast enhanced ultrasound signal as a combination of a low rank and sparse components. This model is used in robust PCA and was previously suggested in the context of ultrasound Doppler processing and dynamic magnetic resonance imaging. We then illustrate that an iterative algorithm based on this model exhibits improved separation of microbubble signal from the tissue signal over commonly practiced methods. Next, we apply the concept of deep unfolding to suggest a deep network architecture tailored to our clutter filtering problem which exhibits improved convergence speed and accuracy with respect to its iterative counterpart. We compare the performance of the suggested deep network on both simulations and in-vivo rat brain scans, with a commonly practiced deep-network architecture and the fast iterative shrinkage algorithm, and show that our architecture exhibits better image quality and contrast.


title: Multithermal-multibaric molecular simulations from a variational principle

abstract: We present a method for performing multithermal-multibaric molecular dynamics simulations that sample entire regions of the temperature-pressure (TP) phase diagram. The method uses a variational principle [Valsson and Parrinello, Phys. Rev. Lett. 113, 090601 (2014)] in order to construct a bias that leads to a uniform sampling in energy and volume. The intervals of temperature and pressure are taken as inputs and the relevant energy and volume regions are determined on the fly. In this way the method guarantees adequate statistics for the chosen TP region. We show that our multithermal-multibaric simulations can be used to calculate all static physical quantities for all temperatures and pressures in the targeted region of the TP plane. We illustrate our approach by studying the density anomaly of TIP4P/Ice water.


title: Geometric polarized log Hodge structures with a base of log rank one

abstract: We prove that a projective vertical exact log smooth morphism of fs log analytic spaces with a base of log rank one yields polarized log Hodge structures in the canonical way.


title: An updated review of (sub-)optimal diversification models

abstract: In the past decade many researchers have proposed new optimal portfolio selection strategies to show that sophisticated diversification can outperform the na\"ive 1/N strategy in out-of-sample benchmarks. Providing an updated review of these models since DeMiguel et al. (2009b), I test sixteen strategies across six empirical datasets to see if indeed progress has been made. However, I find that none of the recently suggested strategies consistently outperforms the 1/N or minimum-variance approach in terms of Sharpe ratio, certainty-equivalent return or turnover. This suggests that simple diversification rules are not in fact inefficient, and gains promised by optimal portfolio choice remain unattainable out-of-sample due to large estimation errors in expected returns. Therefore, further research effort should be devoted to both improving estimation of expected returns, and possibly exploring diversification rules that do not require the estimation of expected returns directly, but also use other available information about the stock characteristics.


title: Bubbling the NHEK

abstract: We build the first family of smooth bubbling microstate geometries that are asymptotic to the near-horizon region of extremal five-dimensional Kerr black holes (NHEK). These black holes arise as extremal non-supersymmetric highly-rotating D1-D5-P solutions in type IIB string theory on T$^4\times$S$^1$. Our solutions are asymptotically NHEK in the UV and end in the IR with a smooth cap. In the context of the Kerr/CFT correspondence, these bubbling geometries are dual to pure states of the 1+1 dimensional chiral conformal field theory dual to NHEK. Since our solutions have a bubbling structure in the IR, they correspond to an IR phase of broken conformal symmetry, and their existence supports the possibility that all the pure states whose counting gives the Kerr black hole entropy correspond to horizonless bulk configurations.


title: FALCON: A Fourier Transform Based Approach for Fast and Secure Convolutional Neural Network Predictions

abstract: Machine learning as a service has been widely deployed to utilize deep neural network models to provide prediction services. However, this raises privacy concerns since clients need to send sensitive information to servers. In this paper, we focus on the scenario where clients want to classify private images with a convolutional neural network model hosted in the server, while both parties keep their data private. We present FALCON, a fast and secure approach for CNN predictions based on Fourier Transform. Our solution enables linear layers of a CNN model to be evaluated simply and efficiently with fully homomorphic encryption. We also introduce the first efficient and privacy-preserving protocol for softmax function, which is an indispensable component in CNNs and has not yet been evaluated in previous works due to its high complexity. We implemented the FALCON and evaluated the performance on real-world CNN models. The experimental results show that FALCON outperforms the best known works in both computation and communication cost.


title: Modeling the interaction of an ultra-high intensity laser pulse with nano-layered flat-top cone targets for ion acceleration

abstract: We propose new nanotargets as nano-layered flat-top cone targets to be used for future laser-ion acceleration experiments at ELI-NP. We study their behaviour at the interaction with an ultra-high intensity laser pulse by performing Particle-In-Cell simulations. We analyze spatio-temporal the electromagnetic field based on the finite-difference time-domain method for a complementary description. We find the optimum diameter of the nanospheres and the proper nano-flat-top foil thickness for which one can obtain monoenergetic beams of very energetic ions with low angular divergence.


title: Extracting the Spectrum by Spatial Filtering

abstract: We show that the spectrum of a flow field can be extracted within a local region by straightforward filtering in physical space. We find that for a flow with a certain level of regularity, the filtering kernel must have a sufficient number of vanishing moments in order for the "filtering spectrum" to be meaningful. Our derivation follows a similar analysis by Perrier et al. 1995 for the wavelet spectrum, where we show that the filtering kernel has to have at least $p$ vanishing moments in order to correctly extract a spectrum $k^{-\alpha}$ with $\alpha < p+2$. For example, any flow with a spectrum shallower than $k^{-3}$ can be extracted by a straightforward average on grid-cells of a stencil. We construct two new "simple stencil" kernels, ${\mathcal M}^{I}$ and ${\mathcal M}^{II}$, with only two and three fixed stencil weight coefficients, respectively, and that have sufficient vanishing moments to allow for extracting spectra steeper than $k^{-3}$. We demonstrate our results using synthetic fields, 2D turbulence from a Direct Numerical Simulation, and 3D turbulence from the JHU Database. Our method guarantees energy conservation and can extract spectra of non-quadratic quantities self-consistently, such as kinetic energy in variable density flows, which the wavelet spectrum cannot. The method can be useful in both simulations and experiments when a straightforward Fourier analysis is not justified, such as within coherent flow structures covering non-rectangular regions, in multi-phase flows, or in geophysical flows on Earth's curved surface.


title: Results from the Pierre Auger Observatory

abstract: I describe some of the results on ultrahigh-energy cosmic rays that have been obtained with the Pierre Auger Observatory. These include measurements of the spectrum, composition and anisotropies. Possible astrophysical scenarios that account for these results are discussed.


title: On double sum generating functions in connection with some classical partition theorems

abstract: We focus on writing closed forms of generating functions for the number of partitions with gap conditions as double sums starting from a combinatorial construction. Some examples of the sets of partitions with gap conditions to be discussed here are the set of Rogers--Ramanujan, G\"ollnitz--Gordon, and little G\"ollnitz partitions. This work also includes finding the finite analogs of the related generating functions and the discussion of some related series and polynomial identities. Additionally, we present a different construction and a double sum representation for the products similar to the ones that appear in the Rogers--Ramanujan identities.


title: Giant dielectric difference in chiral asymmetric bilayers

abstract: Twistronics rooted in the twist operation towards bilayer van der Waals crystals is of both theoretical and technological importance. The realize of the correlated electronic behaviors under this operation encourages enormous effort to the research on magic-angle systems which possess sensitive response to the external field. Here, a giant dielectric difference between 30 plus or minus degree twist case is observed in a typical magnetic system 2H-VSe2 bilayer. It is shown that due to the structural inversion asymmetry in its monolayer, the different stacking of the two cases corresponds to the two kind of valley polarized states: interlayer ferrovalley and interlayer antiferrovalley. Further investigations reveal that such different dielectric response between the two states stems from the different Fermi wave vectors coupled to the electric field. More interestingly, we even obtain the selective circularly polarized optical absorption by tuning the interlayer twist. These findings open an appealing route toward functional 2D materials design for electric and optical devices.


title: A Deep Dive into Blockchain Selfish Mining

abstract: This paper studies a fundamental problem regarding the security of blockchain on how the existence of multiple misbehaving pools influences the profitability of selfish mining. Each selfish miner maintains a private chain and makes it public opportunistically for the purpose of acquiring more rewards incommensurate to his Hashrate. We establish a novel Markov chain model to characterize all the state transitions of public and private chains. The minimum requirement of Hashrate together with the minimum delay of being profitable is derived in close-form. The former reduces to 21.48% with the symmetric selfish miners, while their competition with asymmetric Hashrates puts forward a higher requirement of the profitable threshold. The profitable delay increases with the decrease of the Hashrate of selfish miners, making the mining pools more cautious on performing selfish mining.


title: Transferable Interactiveness Knowledge for Human-Object Interaction Detection

abstract: Human-Object Interaction (HOI) Detection is an important problem to understand how humans interact with objects. In this paper, we explore Interactiveness Knowledge which indicates whether human and object interact with each other or not. We found that interactiveness knowledge can be learned across HOI datasets, regardless of HOI category settings. Our core idea is to exploit an Interactiveness Network to learn the general interactiveness knowledge from multiple HOI datasets and perform Non-Interaction Suppression before HOI classification in inference. On account of the generalization of interactiveness, interactiveness network is a transferable knowledge learner and can be cooperated with any HOI detection models to achieve desirable results. We extensively evaluate the proposed method on HICO-DET and V-COCO datasets. Our framework outperforms state-of-the-art HOI detection results by a great margin, verifying its efficacy and flexibility. Code is available at https://github.com/DirtyHarryLYL/Transferable-Interactiveness-Network.


title: Average optimal cost for the Euclidean TSP in one dimension

abstract: The traveling-salesman problem is one of the most studied combinatorial optimization problems, because of the simplicity in its statement and the difficulty in its solution. We study the traveling salesman problem when the positions of the cities are chosen at random in the unit interval and the cost associated to the travel between two cities is their distance elevated to an arbitrary power $p\in \mathbb{R}$. We characterize the optimal cycle and compute the average optimal cost for every number of cities when the measure used to choose the position of the cities is flat and asymptotically for large number of cities in the other cases. We also show that the optimal cost is a self-averaging quantity, and we test our analytical results with extensive simulations.


title: Asymptotic velocity for four celestial bodies

abstract: Asymptotic velocity is defined as the Ces\`aro limit of velocity. As such, its existence has been proven for bounded interaction potentials. This is known to be wrong in celestial mechanics with four or more bodies. Here we show for a class of pair potentials including the homogeneous ones of degree -a for 0<a<2, that asymptotic velocities exist for up to four bodies, dimension three or larger, for any energy and almost all initial conditions on the energy surface.


title: Better than a lens -- Increasing the signal-to-noise ratio through pupil splitting

abstract: Lenses are designed to fulfill Fermats principle such that all light interferes constructively in its focus, guaranteeing its maximum concentration. It can be shown that imaging via an unmodified full pupil yields the maximum transfer strength for all spatial frequencies transferable by the system. Seemingly also the signal-to-noise ratio (SNR) is optimal. The achievable SNR at a given photon budget is critical especially if that budget is strictly limited as in the case of fluorescence microscopy. In this work we propose a general method which achieves a better SNR for high spatial frequency information of an optical imaging system, without the need to capture more photons. This is achieved by splitting the pupil of an incoherent imaging system such that two (or more) sub-images are simultaneously acquired and computationally recombined. We compare the theoretical performance of split pupil imaging to the non-split scenario and implement the splitting using a tilted elliptical mirror placed at the back-focal-plane (BFP) of a fluorescence widefield microscope.


title: Nilpotent Lie groups and hyperbolic automorphisms

abstract: A connected Lie group admitting an expansive automorphism is known to be nilpotent, but all nilpotent Lie groups do not admit expansive automorphism. In this article, we find sufficient conditions for a class of nilpotent Lie groups to admit expansive automorphism.


title: Human Intention Estimation based on Hidden Markov Model Motion Validation for Safe Flexible Robotized Warehouses

abstract: With the substantial growth of logistics businesses the need for larger warehouses and their automation arises, thus using robots as assistants to human workers is becoming a priority. In order to operate efficiently and safely, robot assistants or the supervising system should recognize human intentions in real-time. Theory of mind (ToM) is an intuitive human conception of other humans' mental state, i.e., beliefs and desires, and how they cause behavior. In this paper we propose a ToM based human intention estimation algorithm for flexible robotized warehouses. We observe human's, i.e., worker's motion and validate it with respect to the goal locations using generalized Voronoi diagram based path planning. These observations are then processed by the proposed hidden Markov model framework which estimates worker intentions in an online manner, capable of handling changing environments. To test the proposed intention estimation we ran experiments in a real-world laboratory warehouse with a worker wearing Microsoft Hololens augmented reality glasses. Furthermore, in order to demonstrate the scalability of the approach to larger warehouses, we propose to use virtual reality digital warehouse twins in order to realistically simulate worker behavior. We conducted intention estimation experiments in the larger warehouse digital twin with up to 24 running robots. We demonstrate that the proposed framework estimates warehouse worker intentions precisely and in the end we discuss the experimental results.


title: Graph Convolutional Neural Networks via Motif-based Attention

abstract: Many real-world problems can be represented as graph-based learning problems. In this paper, we propose a novel framework for learning spatial and attentional convolution neural networks on arbitrary graphs. Different from previous convolutional neural networks on graphs, we first design a motif-matching guided subgraph normalization method to capture neighborhood information. Then we implement subgraph-level self-attentional layers to learn different importances from different subgraphs to solve graph classification problems. Analogous to image-based attentional convolution networks that operate on locally connected and weighted regions of the input, we also extend graph normalization from one-dimensional node sequence to two-dimensional node grid by leveraging motif-matching, and design self-attentional layers without requiring any kinds of cost depending on prior knowledge of the graph structure. Our results on both bioinformatics and social network datasets show that we can significantly improve graph classification benchmarks over traditional graph kernel and existing deep models.


title: A Low-latency Secure Data Outsourcing Scheme for Cloud-WSN

abstract: With the support of cloud computing, large quantities of data collected from various WSN applications can be managed efficiently. However, maintaining data security and efficiency of data processing in cloud-WSN (C-WSN) are important and challenging issues. In this paper, we present an efficient data outsourcing scheme based on CP-ABE, which can not only guarantee secure data access, but also reduce overall data processing time. In our proposed scheme, a large file is divided into several data blocks by data owner (DO) firstly. Then, the data blocks are encrypted and transferred to the cloud server in parallel. For data receiver (DR), data decryption and data transmission is also processed in parallel. In addition, data integrity can be checked by DR without any master key components. The security analysis shows that the proposed scheme can meet the security requirement of C-WSN. By performance evaluation, it shows that our scheme can dramatically improve data processing efficiency compared to the traditional CP-ABE method.


title: Science Pipelines for the Square Kilometre Array

abstract: The Square Kilometre Array (SKA) will be both the largest radio telescope ever constructed and the largest Big Data project in the known Universe. The first phase of the project will generate on the order of 5 zettabytes of data per year. A critical task for the SKA will be its ability to process data for science, which will need to be conducted by science pipelines. Together with polarization data from the LOFAR Multifrequency Snapshot Sky Survey (MSSS), we have been developing a realistic SKA-like science pipeline that can handle the large data volumes generated by LOFAR at 150 MHz. The pipeline uses task-based parallelism to image, detect sources, and perform Faraday Tomography across the entire LOFAR sky. The project thereby provides a unique opportunity to contribute to the technological development of the SKA telescope, while simultaneously enabling cutting-edge scientific results. In this paper, we provide an update on current efforts to develop a science pipeline that can enable tight constraints on the magnetised large-scale structure of the Universe.


title: Security management for backhaul-aware 5G-V2X

abstract: Security is a primary concern for the networks aiming at the utilization of Cellular (C) services for connecting Vehicles to Everything (V2X). At present, C-V2X is observing a paradigm shift from Long Term Evolution (LTE) - Evolved Universal Terrestrial Radio Access Network (E-UTRAN) to Fifth Generation (5G) based functional architecture. However, security and credential management are still concerns to be resolved under 5G-V2X. A sizably voluminous number of key updates and non-availability of sub-functions at the edge cause adscititious overheads and decrement the performance while alarming the possibilities of variants of cyber attacks. In this paper, security management is studied as a principle of sustainability and its tradeoff is evaluated with the number of key-updates required to maintain an authenticated connection of a vehicle to the 5G-terminals keeping intact the security functions at the backhaul. A numerical study is presented to determine the claims and understand the proposed tradeoff.


title: Fullness of crossed products of factors by discrete groups

abstract: Let $M$ be an arbitrary factor and $\sigma : \Gamma \curvearrowright M$ an action of a discrete group. In this paper, we study the fullness of the crossed product $M \rtimes_\sigma \Gamma$. When $\Gamma$ is amenable, we obtain a complete characterization: the crossed product factor $M \rtimes_\sigma \Gamma$ is full if and only if $M$ is full and the quotient map $\overline{\sigma} : \Gamma \rightarrow \mathrm{Out}(M)$ has finite kernel and discrete image. This answers a question of Jones from 1981. When $M$ is full and $\Gamma$ is arbitrary, we give a sufficient condition for $M \rtimes_\sigma \Gamma$ to be full which generalizes both Jones' criterion and Choda's criterion. In particular, we show that if $M$ is any full factor (possibly of type $\mathrm{III}$) and $\Gamma$ is a non-inner amenable group, then the crossed product $M \rtimes_\sigma \Gamma$ is full.


title: Autonomous Extraction of a Hierarchical Structure of Tasks in Reinforcement Learning, A Sequential Associate Rule Mining Approach

abstract: Reinforcement learning (RL) techniques, while often powerful, can suffer from slow learning speeds, particularly in high dimensional spaces. Decomposition of tasks into a hierarchical structure holds the potential to significantly speed up learning, generalization, and transfer learning. However, the current task decomposition techniques often rely on high-level knowledge provided by an expert (e.g. using dynamic Bayesian networks) to extract a hierarchical task structure; which is not necessarily available in autonomous systems. In this paper, we propose a novel method based on Sequential Association Rule Mining that can extract Hierarchical Structure of Tasks in Reinforcement Learning (SARM-HSTRL) in an autonomous manner for both Markov decision processes (MDPs) and factored MDPs. The proposed method leverages association rule mining to discover the causal and temporal relationships among states in different trajectories, and extracts a task hierarchy that captures these relationships among sub-goals as termination conditions of different sub-tasks. We prove that the extracted hierarchical policy offers a hierarchically optimal policy in MDPs and factored MDPs. It should be noted that SARM-HSTRL extracts this hierarchical optimal policy without having dynamic Bayesian networks in scenarios with a single task trajectory and also with multiple tasks' trajectories. Furthermore, it has been theoretically and empirically shown that the extracted hierarchical task structure is consistent with trajectories and provides the most efficient, reliable, and compact structure under appropriate assumptions. The numerical results compare the performance of the proposed SARM-HSTRL method with conventional HRL algorithms in terms of the accuracy in detecting the sub-goals, the validity of the extracted hierarchies, and the speed of learning in several testbeds.


title: Band crossings in $^{168}$Ta: a particle-number conserving analysis

abstract: The structures of two observed high-spin rotational bands in the doubly-odd nucleus ${}^{168}$Ta are investigated using the cranked shell model with pairing correlations treated by a particle-number conserving method, in which the blocking effects are taken into account exactly. The experimental moments of inertia and alignments are reproduced very well by the calculations, which confirms the configuration assignments for these two bands in previous works. The backbending and upbending mechanisms in these two bands are analyzed in detail by calculating the occupation probabilities of each orbital close to the Fermi surface and the contributions of each orbital to the total angular momentum alignments. The investigation shows that the level crossings in the first backbending is the neutron $i_{13/2}$ crossing and the in second upbending is the proton $h_{11/2}$ crossing.


title: Influence of proton bunch parameters on a proton-driven plasma wakefield acceleration experiment

abstract: We use particle-in-cell (PIC) simulations to study the effects of variations of the incoming 400 GeV proton bunch parameters on the amplitude and phase of the wakefields resulting from a seeded self-modulation (SSM) process. We find that these effects are largest during the growth of the SSM, i.e. over the first five to six meters of plasma with an electron density of $7 \times10^{14}$ cm$^{-3}$. However, for variations of any single parameter by $\pm$5%, effects after the SSM saturation point are small. In particular, the phase variations correspond to much less than a quarter wakefield period, making deterministic injection of electrons (or positrons) into the accelerating and focusing phase of the wakefields in principle possible. We use the wakefields from the simulations and a simple test electron model to estimate the same effects on the maximum final energies of electrons injected along the plasma, which are found to be below the initial variations of $\pm$5%. This analysis includes the dephasing of the electrons with respect to the wakefields that is expected during the growth of the SSM. Based on a PIC simulation, we also determine the injection position along the bunch and along the plasma leading to the largest energy gain. For the parameters taken here (ratio of peak beam density to plasma density $n_{b0}/n_0 \approx 0.003$), we find that the optimum position along the proton bunch is at $\xi \approx -1.5 \; \sigma_{zb}$, and that the optimal range for injection along the plasma (for a highest final energy of $\sim$1.6 GeV after 10 m) is 5-6 m.


title: Assessing four Neural Networks on Handwritten Digit Recognition Dataset (MNIST)

abstract: Although the image recognition has been a research topic for many years, many researchers still have a keen interest in it[1]. In some papers[2][3][4], however, there is a tendency to compare models only on one or two datasets, either because of time restraints or because the model is tailored to a specific task. Accordingly, it is hard to understand how well a certain model generalizes across image recognition field[6]. In this paper, we compare four neural networks on MNIST dataset[5] with different division. Among them, three are Convolutional Neural Networks (CNN)[7], Deep Residual Network (ResNet)[2] and Dense Convolutional Network (DenseNet)[3] respectively, and the other is our improvement on CNN baseline through introducing Capsule Network (CapsNet)[1] to image recognition area. We show that the previous models despite do a quite good job in this area, our retrofitting can be applied to get a better performance. The result obtained by CapsNet is an accuracy rate of 99.75\%, and it is the best result published so far. Another inspiring result is that CapsNet only needs a small amount of data to get excellent performance. Finally, we will apply CapsNet's ability to generalize in other image recognition field in the future.


title: Dynamical system analysis of Quintom Dark Energy Model

abstract: The present work deals with dynamical system analysis of a Quintom Model of Dark Energy. By suitable transformation of variables the Einstein field equations are converted to an autonomous system. The critical points are determined and stability of hyperbolic critical points are determined by Hartman-Grobman theorem. To analyze non-hyperbolic critical points different tools (notably Center Manifold theory) are used. Possible bifurcation scenarios have also been explained.


title: Use of Enumerative Combinatorics for proving the applicability of an asymptotic stability result on discrete-time SIS epidemics in complex networks

abstract: In this paper, we justify by the use of Enumerative Combinatorics, that the results obtained in \cite{Alarcon1}, where is analysed the complex dynamics of an epidemic model to identify the nodes that contribute the most to the propagation process and because of that are good candidates to be controlled in the network in order to stabilize the network to reach the extinction state, is applicable in almost all the cases. The model analysed was proposed in \cite{Gomez1} %et al. [Phys.Rev.E 84, 036105(2011)] and results obtained in \cite{Alarcon1} implies that it is not necessary to control all nodes, but only a minimal set of nodes if the topology of the network is not regular. This result could be important in the spirit of considering policies of isolation or quarantine of those nodes to be controlled. Simulation results were presented in \cite{Alarcon1} for large free-scale and regular networks, that corroborate the theoretical findings. In this article we justify the applicability of the controllability result obtained in \cite{Alarcon1} in almost all the cases by means of the use of Combinatorics. {\em Mathematics Subjects Classification}: 05A16,34H20,58E25 {\em Keywords}: Asymptotic Graph Enumeration Problems; Network control; virus spreading.


title: Modeling Tripartite Entanglement in Quantum Protocols using Evolving Entangled Hypergraphs

abstract: The notion of entanglement is the most well-known nonclassical correlation in quantum mechanics and a fundamental resource in quantum information and computation. This correlation, which is displayed by certain classes of quantum states, is of utmost importance when dealing with protocols such as quantum teleportation, cryptography, and quantum key distribution. In this paper, we exploit a classification of tripartite entanglement by introducing the concepts of entangled hypergraph and evolving entangled hypergraph as data structures suitable to model quantum protocols which use entanglement. Finally, we present a few examples to provide applications of this model.


title: Applying the swept rule for solving explicit partial differential equations on heterogeneous computing systems

abstract: Applications that exploit the architectural details of high-performance computing (HPC) systems have become increasingly invaluable in academia and industry over the past two decades. The most important hardware development of the last decade in HPC has been the General Purpose Graphics Processing Unit (GPGPU), a class of massively parallel devices that now contributes the majority of computational power in the top 500 supercomputers. As these systems grow, small costs such as latency---due to the fixed cost of memory accesses and communication---accumulate in a large simulation and become a significant barrier to performance. The swept time-space decomposition rule is a communication-avoiding technique for time-stepping stencil update formulas that attempts to reduce latency costs. This work extends the swept rule by targeting heterogeneous, CPU/GPU architectures representing current and future HPC systems. We compare our approach to a naive decomposition scheme with two test equations using an MPI+CUDA pattern on 40 processes over two nodes containing one GPU. The swept rule produces a factor of 1.9 to 23 speedup for the heat equation and a factor of 1.1 to 2.0 speedup for the Euler equations, using the same processors and work distribution, and with the best possible configurations. These results show the potential effectiveness of the swept rule for different equations and numerical schemes on massively parallel computing systems that incur substantial latency costs.


title: CheMixNet: Mixed DNN Architectures for Predicting Chemical Properties using Multiple Molecular Representations

abstract: SMILES is a linear representation of chemical structures which encodes the connection table, and the stereochemistry of a molecule as a line of text with a grammar structure denoting atoms, bonds, rings and chains, and this information can be used to predict chemical properties. Molecular fingerprints are representations of chemical structures, successfully used in similarity search, clustering, classification, drug discovery, and virtual screening and are a standard and computationally efficient abstract representation where structural features are represented as a bit string. Both SMILES and molecular fingerprints are different representations for describing the structure of a molecule. There exist several predictive models for learning chemical properties based on either SMILES or molecular fingerprints. Here, our goal is to build predictive models that can leverage both these molecular representations. In this work, we present CheMixNet -- a set of neural networks for predicting chemical properties from a mixture of features learned from the two molecular representations -- SMILES as sequences and molecular fingerprints as vector inputs. We demonstrate the efficacy of CheMixNet architectures by evaluating on six different datasets. The proposed CheMixNet models not only outperforms the candidate neural architectures such as contemporary fully connected networks that uses molecular fingerprints and 1-D CNN and RNN models trained SMILES sequences, but also other state-of-the-art architectures such as Chemception and Molecular Graph Convolutions.


title: Feature exploration for almost zero-resource ASR-free keyword spotting using a multilingual bottleneck extractor and correspondence autoencoders

abstract: We compare features for dynamic time warping (DTW) when used to bootstrap keyword spotting (KWS) in an almost zero-resource setting. Such quickly-deployable systems aim to support United Nations (UN) humanitarian relief efforts in parts of Africa with severely under-resourced languages. Our objective is to identify acoustic features that provide acceptable KWS performance in such environments. As supervised resource, we restrict ourselves to a small, easily acquired and independently compiled set of isolated keywords. For feature extraction, a multilingual bottleneck feature (BNF) extractor, trained on well-resourced out-of-domain languages, is integrated with a correspondence autoencoder (CAE) trained on extremely sparse in-domain data. On their own, BNFs and CAE features are shown to achieve a more than 2% absolute performance improvement over baseline MFCCs. However, by using BNFs as input to the CAE, even better performance is achieved, with a more than 11% absolute improvement in ROC AUC over MFCCs and more than twice as many top-10 retrievals for two evaluated languages, English and Luganda. We conclude that integrating BNFs with the CAE allows both large out-of-domain and sparse in-domain resources to be exploited for improved ASR-free keyword spotting.


title: On Conformal Spectral Gap Estimates of the Dirichlet-Laplacian

abstract: We study spectral stability estimates of the Dirichlet eigenvalues of the Laplacian in non-convex domains $\Omega\subset\mathbb R^2$. With the help of these estimates we obtain asymptotically sharp inequalities of ratios of eigenvalues in the frameworks of the Payne-P\'olya-Weinberger inequalities. These estimates are equivalent to spectral gap estimates of the Dirichlet eigenvalues of the Laplacian in non-convex domains in terms of conformal (hyperbolic) geometry.


title: Accelerating the Evolution of Convolutional Neural Networks with Node-Level Mutations and Epigenetic Weight Initialization

abstract: This paper examines three generic strategies for improving the performance of neuro-evolution techniques aimed at evolving convolutional neural networks (CNNs). These were implemented as part of the Evolutionary eXploration of Augmenting Convolutional Topologies (EXACT) algorithm. EXACT evolves arbitrary convolutional neural networks (CNNs) with goals of better discovering and understanding new effective architectures of CNNs for machine learning tasks and to potentially automate the process of network design and selection. The strategies examined are node-level mutation operations, epigenetic weight initialization and pooling connections. Results were gathered over the period of a month using a volunteer computing project, where over 225,000 CNNs were trained and evaluated across 16 different EXACT searches. The node mutation operations where shown to dramatically improve evolution rates over traditional edge mutation operations (as used by the NEAT algorithm), and epigenetic weight initialization was shown to further increase the accuracy and generalizability of the trained CNNs. As a negative but interesting result, allowing for pooling connections was shown to degrade the evolution progress. The best trained CNNs reached 99.46% accuracy on the MNIST test data in under 13,500 CNN evaluations -- accuracy comparable with some of the best human designed CNNs.


title: On-sky performance of the CLASS Q-band telescope

abstract: The Cosmology Large Angular Scale Surveyor (CLASS) is mapping the polarization of the Cosmic Microwave Background (CMB) at large angular scales ($2<\ell\lesssim200$) in search of a primordial gravitational wave B-mode signal down to a tensor-to-scalar ratio of $r \approx 0.01$. The same data set will provide a near sample-variance-limited measurement of the optical depth to reionization. Between June 2016 and March 2018, CLASS completed the largest ground-based Q-band CMB survey to date, covering over 31 000~square-degrees (75% of the sky), with an instantaneous array noise-equivalent temperature (NET) sensitivity of $32~\mu \mbox{K}_{cmb}\sqrt{\mbox{s}}$. We demonstrate that the detector optical loading ($1.6~\mbox{pW}$) and noise-equivalent power ($19~\mbox{aW}\sqrt{\mbox{s}}$) match the expected noise model dominated by photon bunching noise. We derive a $13.1\pm0.3~\mbox{K/pW}$ calibration to antenna temperature based on Moon observations, which translates to an optical efficiency of $0.48\pm0.04$ and a $27~\mbox{K}$ system noise temperature. Finally, we report a Tau A flux density of $308\pm11~\mbox{Jy}$ at $38.4\pm0.2~\mbox{GHz}$, consistent with the WMAP Tau A time-dependent spectral flux density model.


title: The Gamma filtrations for the Spin groups

abstract: Let $G$ be a compact Lie group and $T$ its maximal torus. In this paper, we try to compute $gr_{\gamma}^*(G/T)$ the graded ring associated with the gamma filtration of the complex $K$-theory $K^0(G/T)$ for $G=Spin(n)$. In particular, we give a counterexample for a conjecture by Karpenko when $G=Spin(17)$. The arguments for $E_7$ in $\S 11$ of the old version were not correct, and they are deleted in this version.


title: A Lagrangian Model to Predict Microscallop Motion in non Newtonian Fluids

abstract: The need to develop models to predict the motion of microrobots, or robots of a much smaller scale, moving in fluids in a low Reynolds number regime, and in particular, in non Newtonian fluids, cannot be understated. The article develops a Lagrangian based model for one such mechanism - a two-link mechanism termed a microscallop, moving in a low Reynolds number environment in a non Newtonian fluid. The modelling proceeds through the conventional Lagrangian construction for a two-link mechanism and then goes on to model the external fluid forces using empirically based models for viscosity to complete the dynamic model. The derived model is then simulated for different initial conditions and key parameters of the non Newtonian fluid, and the results are corroborated with a few existing experimental results on a similar mechanism under identical conditions. Lastly, with a view to implementing control algorithms we explore accessibility of the system at certain configurations.


title: An Efficient Optical Flow Based Motion Detection Method for Non-stationary Scenes

abstract: Real-time motion detection in non-stationary scenes is a difficult task due to dynamic background, changing foreground appearance and limited computational resource. These challenges degrade the performance of the existing methods in practical applications. In this paper, an optical flow based framework is proposed to address this problem. By applying a novel strategy to utilize optical flow, we enable our method being free of model constructing, training or updating and can be performed efficiently. Besides, a dual judgment mechanism with adaptive intervals and adaptive thresholds is designed to heighten the system's adaptation to different situations. In experiment part, we quantitatively and qualitatively validate the effectiveness and feasibility of our method with videos in various scene conditions. The experimental results show that our method adapts itself to different situations and outperforms the state-of-the-art real-time methods, indicating the advantages of our optical flow based method.


title: A Two Phase Investment Game for Competitive Opinion Dynamics in Social Networks

abstract: We propose a setting for two-phase opinion dynamics in social networks, where a node's final opinion in the first phase acts as its initial biased opinion in the second phase. In this setting, we study the problem of two camps aiming to maximize adoption of their respective opinions, by strategically investing on nodes in the two phases. A node's initial opinion in the second phase naturally plays a key role in determining the final opinion of that node, and hence also of other nodes in the network due to its influence on them. More importantly, this bias also determines the effectiveness of a camp's investment on that node in the second phase. To formalize this two-phase investment setting, we propose an extension of Friedkin-Johnsen model, and hence formulate the utility functions of the camps. There is a tradeoff while splitting the budget between the two phases. A lower investment in the first phase results in worse initial biases for the second phase, while a higher investment spares a lower available budget for the second phase. We first analyze the non-competitive case where only one camp invests, for which we present a polynomial time algorithm for determining an optimal way to split the camp's budget between the two phases. We then analyze the case of competing camps, where we show the existence of Nash equilibrium and that it can be computed in polynomial time under reasonable assumptions. We conclude our study with simulations on real-world network datasets, in order to quantify the effects of the initial biases and the weightage attributed by nodes to their initial biases, as well as that of a camp deviating from its equilibrium strategy. Our main conclusion is that, if nodes attribute high weightage to their initial biases, it is advantageous to have a high investment in the first phase, so as to effectively influence the biases to be harnessed in the second phase.


title: Freeze-in leptogenesis with 3 right-handed neutrinos

abstract: We provide the first systematic study of the viable parameter space for leptogenesis in the type-I seesaw model with three right-handed neutrinos whose Majorana masses lie below the electroweak scale. We highlight the very rich phenomenology of this scenario and discuss several mechanisms that can help to enhance the baryon asymmetry of the Universe. This allows for much larger heavy neutrino mixing angles than in the minimal scenario with two heavy neutrinos, resulting in solutions that can be probed with the LHC and other experiments. The light neutrino masses can be protected from radiative corrections by an approximate symmetry related to generalised lepton number.


title: Pressure function and limit theorems for almost Anosov flows

abstract: We obtain limit theorems (Stable Laws and Central Limit Theorems, both Gaussian and non-Gaussian) and thermodynamic properties for a class of non-uniformly hyperbolic flows: almost Anosov flows, constructed here. The proofs of the limit theorems for these flows are applications of corresponding theorems in an abstract functional analytic framework, which may be applicable to other classes of non-uniformly hyperbolic flows.


title: Minimal surfaces with mixed three-form flux

abstract: We study minimal area world sheets ending on two concentric circumferences on the boundary of Euclidean $AdS_{3}$ with mixed R-R and NS-NS three-form fluxes. We solve the problem by reducing the system to a one-dimensional integrable model. We find that the NS-NS flux term either brings the surface near to the boundary or separates the circumferences. In the limit of pure NS-NS flux the solution adheres to the boundary in the former case and the outer radius diverges in the latter. We further construct the underlying elliptic spectral curve, which allows us to analyze the deformation of other related minimal surfaces. We show that in the regime of pure NS-NS flux the elliptic curve degenerates.


title: T-CGAN: Conditional Generative Adversarial Network for Data Augmentation in Noisy Time Series with Irregular Sampling

abstract: In this paper we propose a data augmentation method for time series with irregular sampling, Time-Conditional Generative Adversarial Network (T-CGAN). Our approach is based on Conditional Generative Adversarial Networks (CGAN), where the generative step is implemented by a deconvolutional NN and the discriminative step by a convolutional NN. Both the generator and the discriminator are conditioned on the sampling timestamps, to learn the hidden relationship between data and timestamps, and consequently to generate new time series. We evaluate our model with synthetic and real-world datasets. For the synthetic data, we compare the performance of a classifier trained with T-CGAN-generated data, against the performance of the same classifier trained on the original data. Results show that classifiers trained on T-CGAN-generated data perform the same as classifiers trained on real data, even with very short time series and small training sets. For the real world datasets, we compare our method with other techniques of data augmentation for time series, such as time slicing and time warping, over a classification problem with unbalanced datasets. Results show that our method always outperforms the other approaches, both in case of regularly sampled and irregularly sampled time series. We achieve particularly good performance in case with a small training set and short, noisy, irregularly-sampled time series.


title: A passive Stokes flow rectifier for Newtonian fluids

abstract: Non-linear effects of the Navier-Stokes equations disappear under the Stokes regime of Newtonian fluid flows disallowing the fluid flow rectification. Here we show mathematically and experimentally that passive flow rectification of Newtonian fluids is obtainable under the Stokes regime of both compressible and incompressible flows by introducing nonlinearity into the otherwise linear Stokes equations. Asymmetric flow resistances arise in shallow nozzle/diffuser microchannels with deformable ceiling, in which the fluid flow is governed by a non-linear coupled fluid-solid mechanics equation. Fluid flow rectification has been demonstrated for low-Reynolds-number flows (Re ~ O(0.001)-O(1)) of common Newtonian fluids such as air, water, and alcohol. This mechanism can pave the way for regulating the low-Reynolds-number fluid flows with potential applications in precise low-flow-rate micropumps, drug delivery systems, etc.


title: Finite Mixture Model of Nonparametric Density Estimation using Sampling Importance Resampling for Persistence Landscape

abstract: Considering the creation of persistence landscape on a parametrized curve and structure of sampling, there exists a random process for which a finite mixture model of persistence landscape (FMMPL) can provide a better description for a given dataset. In this paper, a nonparametric approach for computing integrated mean of square error (IMSE) in persistence landscape has been presented. As a result, FMMPL is more accurate than the another way. Also, the sampling importance resampling (SIR) has been presented a better description of important landmark from parametrized curve. The result, provides more accuracy and less space complexity than the landmarks selected with simple sampling.


title: A novel derivation of the Marchenko-Pastur law through analog bipartite spin-glasses

abstract: In this work we consider the {\em analog bipartite spin-glass} (or {\em real-valued restricted Boltzmann machine} in a neural network jargon), whose variables (those quenched as well as those dynamical) share standard Gaussian distributions. First, via Guerra's interpolation technique, we express its quenched free energy in terms of the natural order parameters of the theory (namely the self- and two-replica overlaps), then, we re-obtain the same result by using the replica-trick: a mandatory tribute, given the special occasion. Next, we show that the quenched free energy of this model is the functional generator of the moments of the correlation matrix among the weights connecting the two layers of the spin-glass (i.e., the Wishart matrix in random matrix theory or the Hebbian coupling in neural networks): as weights are quenched stochastic variables, this plays as a novel tool to inspect random matrices. In particular, we find that the Stieltjes transform of the spectral density of the correlation matrix is determined by the (replica-symmetric) quenched free energy of the bipartite spin-glass model. In this setup, we re-obtain the Marchenko-Pastur law in a very simple way.


title: Exact rank-reduction of network models

abstract: With the advent of the big data era, generative models of complex networks are becoming elusive from direct computational simulation. We present an exact, linear-algebraic reduction scheme of generative models of networks. By exploiting the bilinear structure of the matrix representation of the generative model, we separate its null eigenspace, and reduce the exact description of the generative model to a smaller vector space. After reduction, we group generative models in universality classes according to their rank and metric signature, and work out, in a computationally affordable way, their relevant properties (e.g., spectrum). The reduction also provides the environment for a simplified computation of their properties. The proposed scheme works for any generative model admitting a matrix representation, and will be very useful in the study of dynamical processes on networks, as well as in the understanding of generative models to come, according to the provided classification.


title: Theoretical Constraints on Supersymmetric Models: Perturbative Unitarity vs. Vacuum Stability

abstract: There are nowadays strong experimental constraints on supersymmetric theories from the Higgs measurements as well as from the null results in Sparticle searches. However, even the parameter spaces which are in agreement with experimental data can be further constrained by using theoretical considerations. Here, we discuss for the MSSM and NMSSM the impact of perturbative unitarity as well as of the stability of the one-loop effective potential. We find in the case of the MSSM, that vacuum stability is always the stronger constraint. On the other side, the situation is more diverse in the NMSSM and one should always check both kind of constraints.


title: Fingerprints of tilted Dirac cones on the RKKY exchange interaction in 8-Pmmn Borophene

abstract: We theoretically investigate the indirect signatures of the tilted anisotropic Dirac cones on Ruderman-Kittel-Kasuya-Yosida (RKKY) exchange interaction in a two dimensional polymorph of boron atoms. Unlike the case of isotropic non-tilted Dirac material-graphene, here we observe that the tilting of the Dirac cones exhibits a significant impact on the RKKY exchange interaction in terms of the suppression of oscillation frequency. The reason can be attributed to the behavior of the Fermi level and the corresponding density of states with respect to the tilting parameter. The direct measurement of the period of the RKKY interaction can thus be a possible probe of the tilt parameter associated with the tilted Dirac cones. We also obtain the direction dependent analytical expressions of the RKKY exchange interaction, in terms of Meijer G-function. However, the effects due to tilting of the Dirac cones on the RKKY interaction depend on the spatial alignments of the two magnetic impurities with respect to the direction of tilting.


title: How Quantum is the Speedup in Adiabatic Unstructured Search?

abstract: In classical computing, analog approaches have sometimes appeared to be more powerful than they really are. This occurs when resources, particularly precision, are not appropriately taken into account. While the same should also hold for analog quantum computing, precision issues are often neglected from the analysis. In this work we present a classical analog algorithm for unstructured search that can be viewed as analogous to the quantum adiabatic unstructured search algorithm devised by Roland and Cerf [Phys. Rev. A 65, 042308 (2002)]. We show that similarly to its quantum counterpart, the classical construction may also provide a quadratic speedup over standard digital unstructured search. We discuss the meaning and the possible implications of this result in the context of adiabatic quantum computing.


title: Abelian Decomposition and Weyl Symmetric Effective Action of SU(3) QCD

abstract: We show how to calculate the effective potential of SU(3) QCD which tells that the true minimum is given by the monopole condensation. To do this we make the gauge independent Weyl symmetric Abelian decomposition of the SU(3) QCD which decomposes the gluons to the color neutral neurons and the colored chromons. In the perturbative regime this decomposes the Feynman diagram in such a way that the conservation of color is explicit. Moreover, this shows the existence of two gluon jets, the neuron jet and chromon jet, which can be verified by experiment. In the non-perturbative regime, the decomposition puts QCD to the background field formalism and reduces the non-Abelian gauge symmetry to a discrete color reflection symmetry, and provides us an ideal platform to calculate the one-loop effective action of QCD. Integrating out the chromons from the Weyl symmetric Abelian decomposition of QCD gauge invariantly imposing the color reflection invariance, we obtain the SU(3) QCD effective potential which generates the stable monopole condensation and the mass gap. We discuss the physical implications of our result, in particular the possible existence of the vacuum fluctuation mode of the monopole condensation in QCD.


title: Search for $CP$ violation through an amplitude analysis of $D^0 \rightarrow K^+ K^- \pi^+ \pi^-$ decays

abstract: A search for $CP$ violation in the Cabibbo-suppressed $D^0 \rightarrow K^+ K^- \pi^+ \pi^-$ decay mode is performed using an amplitude analysis. The measurement uses a sample of $pp$ collisions recorded by the LHCb experiment during 2011 and 2012, corresponding to an integrated luminosity of 3.0 fb$^{-1}$. The $D^0$ mesons are reconstructed from semileptonic $b$-hadron decays into $D^0\mu^- X$ final states. The selected sample contains more than 160000 signal decays, allowing the most precise amplitude modelling of this $D^0$ decay to date. The obtained amplitude model is used to perform the search for $CP$ violation. The result is compatible with $CP$ symmetry, with a sensitivity ranging from 1% to 15% depending on the amplitude considered.


title: IVD-Net: Intervertebral disc localization and segmentation in MRI with a multi-modal UNet

abstract: Accurate localization and segmentation of intervertebral disc (IVD) is crucial for the assessment of spine disease diagnosis. Despite the technological advances in medical imaging, IVD localization and segmentation are still manually performed, which is time-consuming and prone to errors. If, in addition, multi-modal imaging is considered, the burden imposed on disease assessments increases substantially. In this paper, we propose an architecture for IVD localization and segmentation in multi-modal MRI, which extends the well-known UNet. Compared to single images, multi-modal data brings complementary information, contributing to better data representation and discriminative power. Our contributions are three-fold. First, how to effectively integrate and fully leverage multi-modal data remains almost unexplored. In this work, each MRI modality is processed in a different path to better exploit their unique information. Second, inspired by HyperDenseNet, the network is densely-connected both within each path and across different paths, granting the model the freedom to learn where and how the different modalities should be processed and combined. Third, we improved standard U-Net modules by extending inception modules with two dilated convolutions blocks of different scale, which helps handling multi-scale context. We report experiments over the data set of the public MICCAI 2018 Challenge on Automatic Intervertebral Disc Localization and Segmentation, with 13 multi-modal MRI images used for training and 3 for validation. We trained IVD-Net on an NVidia TITAN XP GPU with 16 GBs RAM, using ADAM as optimizer and a learning rate of 10e-5 during 200 epochs. Training took about 5 hours, and segmentation of a whole volume about 2-3 seconds, on average. Several baselines, with different multi-modal fusion strategies, were used to demonstrate the effectiveness of the proposed architecture.


title: Study of Multi-Step Knowledge-Aided Iterative Nested MUSIC for Direction Finding

abstract: In this work, we propose a subspace-based algorithm for direction-of-arrival (DOA) estimation applied to the signals impinging on a two-level nested array, referred to as multi-step knowledge-aided iterative nested MUSIC method (MS-KAI-Nested-MUSIC), which significantly improves the accuracy of the original Nested-MUSIC. Differently from existing knowledge-aided methods applied to uniform linear arrays (ULAs), which make use of available known DOAs to improve the estimation of the covariance matrix of the input data, the proposed Multi-Step KAI-Nested-MU employs knowledge of the structure of the augmented sample covariance matrix, which is obtained by exploiting the difference co-array structure covariance matrix, and its perturbation terms and the gradual incorporation of prior knowledge, which is obtained on line. The effectiveness of the proposed technique can be noticed by simulations focusing on uncorrelated closely-spaced sources.


title: A Criterion for the Existence of Relaxation Oscillations with Applications to Predator-Prey Systems and an Epidemic Model

abstract: We derive characteristic functions to determine the number and stability of relaxation oscillations for a class of planar systems. Applying our criterion, we give conditions under which the chemostat predator-prey system has a globally orbitally asymptotically stable limit cycle. Also we demonstrate that a prescribed number of relaxation oscillations can be constructed by varying the perturbation for an epidemic model studied by Li et al. [SIAM J. Appl. Math, 2016].


title: Economics of disagreement -- financial intuition for the R\'enyi divergence

abstract: Disagreement is an essential element of science and life in general. The language of probabilities and statistics is often used to describe disagreements quantitatively. In practice, however, we want much more than that. We want disagreements to be resolved. This leaves us with a substantial knowledge gap which is often perceived as a lack of practical intuition regarding probabilistic and statistical concepts. Take for instance the R\'enyi divergence which is a well-known statistical quantity specifically designed as a measure of disagreement between probabilistic models. Despite its widespread use in science and engineering, the R\'enyi divergence remains a highly abstract axiomatically-motivated measure. Certainly, it offers no practical insight as to how disagreements can be resolved. Here we propose to address disagreements using the methods of financial economics. In particular, we show how a large class of disagreements can be transformed into investment opportunities. The expected financial performance of such investments quantifies the amount of disagreement in a tangible way. This provides intuition for statistical concepts such as the R\'enyi divergence which becomes connected to the financial performance of optimized investments. Investment optimization takes into account individual opinions as well as attitudes towards risk. The result is a market-like social mechanism by which funds flow naturally to support a more accurate view. Such social mechanisms can help us with difficult disagreements (e.g., financial arguments concerning the future climate). In terms of scientific validation, we used the findings of independent neurophysiological experiments as well as our own research on the equity premium.

categories: q-fin.GN cs.GT cs.IT math.IT math.ST stat.TH


title: Modeling Deep Learning Accelerator Enabled GPUs

abstract: The efficacy of deep learning has resulted in its use in a growing number of applications. The Volta graphics processor unit (GPU) architecture from NVIDIA introduced a specialized functional unit, the "tensor core", that helps meet the growing demand for higher performance for deep learning. In this paper we study the design of the tensor cores in NVIDIA's Volta and Turing architectures. We further propose an architectural model for the tensor cores in Volta. When implemented a GPU simulator, GPGPU-Sim, our tensor core model achieves 99.6\% correlation versus an NVIDIA Titan~V GPU in terms of average instructions per cycle when running tensor core enabled GEMM workloads. We also describe support added to enable GPGPU-Sim to run CUTLASS, an open-source CUDA C++ template library providing customizable GEMM templates that utilize tensor cores.


title: Selective chaos of travelling waves in feedforward chains of bistable maps

abstract: We study the chaos of travelling waves (TW) in unidirectional chains of bistable maps. Previous numerical results suggested that this property is selective, {\sl viz.}\ given the parameters, there is at most a single (non-trivial) velocity for which the corresponding set of wave profiles has positive topological entropy. However, mathematical proofs have remained elusive, in particular because the related symbolic dynamics involves entire past sequences. Here, we consider instead inite (short) rank approximations for which the symbolic dynamics has finite memory. For every possible velocity, we compute the existence domains of all possible finite type subshifts of TW with positive entropy. In all examples, chaos of TW turns out to be selective, indeed.


title: PHENIX results on collision energy dependent Levy HBT correlations from $\sqrt{s_{NN}}$ = 15 to 200 GeV

abstract: Different regions on the QCD phase diagram can be investigated by varying the collision energy and the centrality in heavy-ion collisions. In our latest measurements at the PHENIX experiment at RHIC, we utilize L\'evy-type sources to describe the measured HBT correlation functions. In this paper we report the current status of the analysis of the centrality and beam energy dependence of the L\'evy source parameters in Au+Au collisions from $\sqrt{s_{NN}}$ = 15 GeV to $\sqrt{s_{NN}}$ = 200 GeV.


title: Review of Calculators for BSM Higgs bosons

abstract: We have reached a new era of particle physics in which the properties of the Higgs boson, in particular its mass, turned into precision observables. Therefore, it is necessary to have accurate predictions of these properties in models for new physics. I give an overview of available tools for supersymmetric and non-supersymmetric models which are publically available. Afterwards, the main idea behind generic tools to study a large variety of models is summarised. Also some remarks about the validity of checks for perturbative unitarity, which are often applied for non-supersymmetric models, are given.


title: Two-temperatures overlap distribution for the 2D discrete Gaussian free field

abstract: In this paper, we prove absence of temperature chaos for the two-dimensional discrete Gaussian free field using the convergence of the full extremal process, which has been obtained recently by Biskup and Louidor. This means that the overlap of two points chosen under Gibbs measures at different temperatures has a nontrivial distribution. Whereas this distribution is the same as for the random energy model when the two points are sampled at the same temperature, we point out here that they are different when temperatures are distinct: more precisely, we prove that the mean overlap of two points chosen under Gibbs measures at different temperatures for the DGFF is strictly smaller than the REM's one. Therefore, although neither of these models exhibits temperature chaos, one could say that the DGFF is more chaotic in temperature than the REM.


title: Counting Words Avoiding a Short Increasing Pattern and the Pattern 1k...2

abstract: We find finite-state recurrences to enumerate the words on the alphabet $[n]^r$ which avoid the patterns 123 and $1k(k-1)\dots2$, and, separately, the words which avoid the patterns 1234 and $1k(k-1)\dots2$.


title: On the Isochronous analytic motions and the quantum spectrum

abstract: The problem of the characterization of all analytic potentials which give rise to isochronous oscillatory motions still open. However, there are several approaches to highlight motions with period $T(E) \equiv T_0$ independent on the energy. It is proposed in this paper to give necessary and sufficient conditions for a center to be isochronous. The proofs produced here are self contained. As corollaries we find again lot of known conditions previously established by different authors. This allows us to produce several class of isochronous analytic motions. We are also interested in the quantum spectrum. We then use the perturbation method WKB to derive an expression for the corrections to the equally spaced valid for analytic isochronous potentials.


title: Constructing high-order discontinuity-capturing schemes with linear-weight polynomials and boundary variation diminishing algorithm

abstract: In this study, a new framework of constructing very high order discontinuity-capturing schemes is proposed for finite volume method. These schemes, so-called $\mathrm{P}_{n}\mathrm{T}_{m}-\mathrm{BVD}$ (polynomial of $n$-degree and THINC function of $m$-level reconstruction based on BVD algorithm), are designed by employing high-order linear-weight polynomials and THINC (Tangent of Hyperbola for INterface Capturing) functions with adaptive steepness as the reconstruction candidates. The final reconstruction function in each cell is determined with a multi-stage BVD (Boundary Variation Diminishing) algorithm so as to effectively control numerical oscillation and dissipation. We devise the new schemes up to eleventh order in an efficient way by directly increasing the order of the underlying upwind scheme using linear-weight polynomials. The analysis of the spectral property and accuracy tests show that the new reconstruction strategy well preserves the low-dissipation property of the underlying upwind schemes with high-order linear-weight polynomials for smooth solution over all wave numbers and realizes $n+1$ order convergence rate. The performance of new schemes is examined through widely used benchmark tests, which demonstrate that the proposed schemes are capable of simultaneously resolving small-scale flow features with high resolution and capturing discontinuities with low dissipation. With outperforming results and simplicity in algorithm, the new reconstruction strategy shows great potential as an alternative numerical framework for computing nonlinear hyperbolic conservation laws that have discontinuous and smooth solutions of different scales.


title: A power law decay evolution scenario for polluted single white dwarfs

abstract: Planetary systems can survive the stellar evolution, as evidenced by the atmospheric metal pollution and dusty disks of single white dwarfs. Recent observations show that 1 to 4 percent of single white dwarfs are accompanied by dusty disks, while the occurrence rate of metal pollution is about 25 to 50 percent . The dusty disks and metal pollution have been associated with accretion of remnant planetary systems around white dwarfs, yet the relation between these two phenomena is still unclear. Here we suggest an evolutionary scenario to link the two observational phenomena. By analyzing a sample of metal polluted white dwarfs, we find that the mass accretion rate onto the white dwarf generally follows a broken power law decay, which matches well with the theoretical prediction, if assuming dust accretion is primarily driven by Poynting-Robertson drag and the dust source is primarily delivered via dynamically falling asteroids perturbed by a Jovian planet. The presence of disks is mainly at the early stage (about 0.1 to 0.7 billion years) of the whole process of metal pollution, which is detectable until about 8 billion years, naturally explaining the fraction (about 2-16 percent) of metal-polluted white dwarfs having dusty disks. The success of this scenario also implies that the configuration of an asteroid belt with an outer gas giant might be common around stars of several solar masses.


title: Self-Organizing Maps for Storage and Transfer of Knowledge in Reinforcement Learning

abstract: The idea of reusing or transferring information from previously learned tasks (source tasks) for the learning of new tasks (target tasks) has the potential to significantly improve the sample efficiency of a reinforcement learning agent. In this work, we describe a novel approach for reusing previously acquired knowledge by using it to guide the exploration of an agent while it learns new tasks. In order to do so, we employ a variant of the growing self-organizing map algorithm, which is trained using a measure of similarity that is defined directly in the space of the vectorized representations of the value functions. In addition to enabling transfer across tasks, the resulting map is simultaneously used to enable the efficient storage of previously acquired task knowledge in an adaptive and scalable manner. We empirically validate our approach in a simulated navigation environment, and also demonstrate its utility through simple experiments using a mobile micro-robotics platform. In addition, we demonstrate the scalability of this approach, and analytically examine its relation to the proposed network growth mechanism. Further, we briefly discuss some of the possible improvements and extensions to this approach, as well as its relevance to real world scenarios in the context of continual learning.


title: Advances in Reduced Order Methods for Parametric Industrial Problems in Computational Fluid Dynamics

abstract: Reduced order modeling has gained considerable attention in recent decades owing to the advantages offered in reduced computational times and multiple solutions for parametric problems. The focus of this manuscript is the application of model order reduction techniques in various engineering and scientific applications including but not limited to mechanical, naval and aeronautical engineering. The focus here is kept limited to computational fluid mechanics and related applications. The advances in the reduced order modeling with proper orthogonal decomposition and reduced basis method are presented as well as a brief discussion of dynamic mode decomposition and also some present advances in the parameter space reduction. Here, an overview of the challenges faced and possible solutions are presented with examples from various problems.


abstract: Hadronic light-by-light scattering in the anomalous magnetic moment of the muon $a_\mu$ is one of two hadronic effects limiting the precision of the Standard Model prediction for this precision observable, and hence the new-physics discovery potential of direct experimental determinations of $a_\mu$. In this contribution, we report on recent progress in the calculation of this effect achieved both via dispersive and lattice QCD methods.


title: Stability Based Filter Pruning for Accelerating Deep CNNs

abstract: Convolutional neural networks (CNN) have achieved impressive performance on the wide variety of tasks (classification, detection, etc.) across multiple domains at the cost of high computational and memory requirements. Thus, leveraging CNNs for real-time applications necessitates model compression approaches that not only reduce the total number of parameters but reduce the overall computation as well. In this work, we present a stability-based approach for filter-level pruning of CNNs. We evaluate our proposed approach on different architectures (LeNet, VGG-16, ResNet, and Faster RCNN) and datasets and demonstrate its generalizability through extensive experiments. Moreover, our compressed models can be used at run-time without requiring any special libraries or hardware. Our model compression method reduces the number of FLOPS by an impressive factor of 6.03X and GPU memory footprint by more than 17X, significantly outperforming other state-of-the-art filter pruning methods.


title: The generalized distance matrix of digraphs

abstract: Let $D(G)$ and $D^Q(G)= Diag(Tr) + D(G)$ be the distance matrix and distance signless Laplacian matrix of a simple strongly connected digraph $G$, respectively, where $Diag(Tr)=\textrm{diag}(D_1,D_2,$ $\ldots,D_n)$ be the diagonal matrix with vertex transmissions of the digraph $G$. To track the gradual change of $D(G)$ into $D^Q(G)$, in this paper, we propose to study the convex combinations of $D(G)$ and $Diag(Tr)$ defined by $$D_\alpha(G)=\alpha Diag(Tr)+(1-\alpha)D(G), \ \ 0\leq \alpha\leq1.$$ This study reduces to merging the distance spectral and distance signless Laplacian spectral theories. The eigenvalue with the largest modulus of $D_\alpha(G)$ is called the $D_\alpha$ spectral radius of $G$, denoted by $\mu_\alpha(G)$. We determine the digraph which attains the maximum (or minimum) $D_\alpha$ spectral radius among all strongly connected digraphs. Moreover, we also determine the digraphs which attain the minimum $D_\alpha$ spectral radius among all strongly connected digraphs with given parameters such as dichromatic number, vertex connectivity or arc connectivity.


title: On the Alexander Theorem for the oriented Thompson group $\vec{F}$

abstract: In [Jo14] and [Jo18] Vaughan Jones introduced a construction which yields oriented knots and links from elements of the oriented Thompson group $\vec{F}$. In this paper we prove, by analogy with Alexander's classical theorem establishing that every knot or link can be represented as a closed braid, that given an oriented knot/link $\vec{L}$, there exists an element $g$ in $\vec{F}$ whose closure $\vec{\mathcal{L}}(g)$ is $\vec L$.


title: Well-posedness for a system of quadratic derivative nonlinear Schr\"odinger equations with radial initial data

abstract: In the present paper, we consider the Cauchy problem of the system of quadratic derivative nonlinear Schr\"odinger equations. This system was introduced by M. Colin and T. Colin (2004). The first and second authors obtained some well-posedness results in the Sobolev space $H^{s}(\mathbb{R}^d)$. We improve these results for conditional radial initial data by rewriting the system radial form.


title: Geometrical properties of the space of idempotent probability measures

abstract: We establish some geometrical properties of the space of idempotent probability measures. In particular, for a compact $X$ it is established that if the space $I_{3}(X)\backslash X$ is hereditary normally, then $X$ is metrizable; some subsets allocate of the space of idempotent probability measures which are, respectively, $Z$-sets, $\max$-$\mbox{plus}$-convex subsets, $G_\delta$-sets.


title: Photonic Jackiw-Rebbi states in all-dielectric structures controlled by bianisotropy

abstract: Electric and magnetic resonances of dielectric particles have recently uncovered a range of exciting applications in steering of light at the nanoscale. Breaking of particle inversion symmetry further modifies its electromagnetic response giving rise to bianisotropy known also as magneto-electric coupling. Recent studies suggest the crucial role of magneto-electric coupling in realization of photonic topological metamaterials. To further unmask this fundamental link, we design and test experimentally one-dimensional array composed of dielectric particles with overlapping electric and magnetic resonances and broken mirror symmetry. Flipping over half of the meta-atoms in the array, we observe the emergence of interface states providing photonic realization of the celebrated Jackiw-Rebbi model. We trace the origin of these states to the fact that local modification of particle bianisotropic response affects its effective coupling with the neighboring meta-atoms which provides a promising avenue to engineer topological states of light.


title: Solution Formulas for Differential Sylvester and Lyapunov Equations

abstract: The differential Sylvester equation and its symmetric version, the differential Lyapunov equation, appear in different fields of applied mathematics like control theory, system theory, and model order reduction. The few available straight-forward numerical approaches if applied to large-scale systems come with prohibitively large storage requirements. This shortage motivates us to summarize and explore existing solution formulas for these equations. We develop a unifying approach based on the spectral theorem for normal operators like the Sylvester operator $\mathcal S(X)=AX+XB$ and derive a formula for its norm using an induced operator norm based on the spectrum of $A$ and $B$. In view of numerical approximations, we propose an algorithm that identifies a suitable Krylov subspace using Taylor series and use a projection to approximate the solution. Numerical results for large-scale differential Lyapunov equations are presented in the last sections.


title: Sensor Adaptation for Improved Semantic Segmentation of Overhead Imagery

abstract: Semantic segmentation is a powerful method to facilitate visual scene understanding. Each pixel is assigned a label according to a pre-defined list of object classes and semantic entities. This becomes very useful as a means to summarize large scale overhead imagery. In this paper we present our work on semantic segmentation with applications to overhead imagery. We propose an algorithm that builds and extends upon the DeepLab framework to be able to refine and resolve small objects (relative to the image size) such as vehicles. We have also investigated sensor adaptation as a means to augment available training data to be able to reduce some of the shortcomings of neural networks when deployed in new environments and to new sensors. We report results on several datasets and compare performance with other state-of-the-art architectures.


title: Unfolding overlapped slitless imaging spectrometer data for extended sources

abstract: Slitless spectrometers can provide simultaneous imaging and spectral data over an extended field of view, thereby allowing rapid data acquisition for extended sources. In some instances, when the object is greatly extended or the spectral dispersion is too small, there may be locations in the focal plane where contributions from emission lines at different wavelengths contribute. It is then desirable to unfold the overlapped regions in order to isolate the contributions from the individual wavelengths. In this paper, we describe a method for such an unfolding, using an inversion technique developed for an extreme ultraviolet imaging spectrometer and coronagraph named the COronal Spectroscopic Imager in the EUV (COSIE). The COSIE spectrometer wavelength range (18.6 - 20.5 nm) contains a number of strong coronal emission lines and several density sensitive lines. We focus on optimizing the unfolding process to retrieve emission measure maps at constant temperature, maps of spectrally pure intensity in the Fe XII and Fe XIII lines and density maps based on both Fe XII and Fe XIII diagnostics.


title: Automatic Test Improvement with DSpot: a Study with Ten Mature Open-Source Projects

abstract: In the literature, there is a rather clear segregation between manually written tests by developers and automatically generated ones. In this paper, we explore a third solution: to automatically improve existing test cases written by developers. We present the concept, design, and implementation of a system called \dspot, that takes developer-written test cases as input (junit tests in Java) and synthesizes improved versions of them as output. Those test improvements are given back to developers as patches or pull requests, that can be directly integrated in the main branch of the test code base. We have evaluated DSpot in a deep, systematic manner over 40 real-world unit test classes from 10 notable and open-source software projects. We have amplified all test methods from those 40 unit test classes. In 26/40 cases, DSpot is able to automatically improve the test under study, by triggering new behaviors and adding new valuable assertions. Next, for ten projects under consideration, we have proposed a test improvement automatically synthesized by \dspot to the lead developers. In total, 13/19 proposed test improvements were accepted by the developers and merged into the main code base. This shows that DSpot is capable of automatically improving unit-tests in real-world, large-scale Java software.


title: Updated Core Libraries of the ALPS Project

abstract: The open source ALPS (Algorithms and Libraries for Physics Simulations) project provides a collection of physics libraries and applications, with a focus on simulations of lattice models and strongly correlated electron systems. The libraries provide a convenient set of well-documented and reusable components for developing condensed matter physics simulation codes, and the applications strive to make commonly used and proven computational algorithms available to a non-expert community. In this paper we present an update of the core ALPS libraries. We present in particular new Monte Carlo libraries and new Green's function libraries.


title: muCool: A next step towards efficient muon beam compression

abstract: A novel device to compress the phase space of a muon beam by a factor of $10^{10}$ with a $10^{-3}$ efficiency is under development. A surface muon beam is stopped in a helium gas target consisting of several compression stages, wherein strong electric and magnetic fields are applied. The spatial extent of the stopped muon swarm is decreased by means of these fields until muons with eV energy are extracted into vacuum through a small orifice. It was measured that a 20 cm long muon stop distribution can be compressed in longitudinal direction to sub-mm extent within 2 ${\rm \mu s}$. Additionally, a drift perpendicular to the magnetic field of the compressed low-energy muon swarm was successfully demonstrated, paving the way towards the extraction from the gas and re-acceleration of the muons.


title: The construction of $\Psi^\star$-algebra by commutator method containing the Bergman projection on the unit ball

abstract: The manuscript is devoted to the construction of $\Psi^\star$- algebras containing the Bergman projection on the unit ball $B_n$ of $\mathbb C^n$. We consider the $C^\star$-algebra $\mathcal{L}(L^2(B_n))$ of bounded operator acting on the Hilbert space of square integrable functions on $B_n$ with respect to the standard probability measure. We search for spectral invariant Fr\'{e}chet subalgebras of $\mathcal{L}(L^2(B_n))$ containing the Bergman projection $P$ which is defined on $L^2(B_n)$ with values in the space of holomorphic functions. We use the commutator method as introduced by Gramsch. This method generalizes the work of Beals for the spectral invariance of pseudodifferential operators. One need not to work on the H\"{o}rmander classes but can describe the algebra of pseudodifferential operators as continuous commutator between scales of Sobolev spaces. To connect the analytical properties to the geometry of the unit ball, we search for linear tangent vector fields $X$ on the the unit sphere $\partial B_n$ so that the commutator $[X,P]$ has a continuous extension to $L^2(B_n)$. In comparison to the work of Bauer in for the case of the Fischer-Fock space it turns out that the results are completely different. Our first contribution states that the set of vector fields which provides the $\Psi^\star$-algebra for the Fock space fails even to provide a continuity with the Bergman projection in the case of the unit ball. Our second contribution, is the construction of vector fields being tangent to the unit sphere which provides a new $\Psi^\star$-algebra containing the Bergman projection of the unit ball. We prove that every linear vector field on the unit sphere commutes with the Bergman projection on the unit ball.


title: An asymptotic resolution of a problem of Plesn\'{i}k

abstract: Fix $d \ge 3$. We show the existence of a constant $c>0$ such that any graph of diameter at most $d$ has average distance at most $d-c \frac{d^{3/2}}{\sqrt n}$, where $n$ is the number of vertices. Moreover, we exhibit graphs certifying sharpness of this bound up to the choice of $c$. This constitutes an asymptotic solution to a longstanding open problem of Plesn\'{i}k. Furthermore we solve the problem exactly for digraphs if the order is large compared with the diameter.


title: Renormalization in Some 2D su(2) Coset Models

abstract: We consider a RG flow in certain 2D coset models perturbed by the least relevant field. In the case of the symmetric su(2) coset model we show, up to second order of the perturbation theory, that there exists a nontrivial IR fixed point.We obtain the structure constants and the four-point functions of certain fields by deriving specific recursive relations. This allows us to compute the anomalous dimensions and the mixing coefficients of these fields in the UV and IR theories. In the case of another su(2) coset model, describing the N=2 superconformal theories, we show that there does not exists a nontrivial IR fixed point up to second order.


title: Identifiability from a few species for a class of biochemical reaction networks

abstract: Under mass-action kinetics, biochemical reaction networks give rise to polynomial autonomous dynamical systems whose parameters are often difficult to estimate. We deal in this paper with the problem of identifying the kinetic parameters of a class of biochemical networks which are abundant, such as multisite phosphorylation systems and phosphorylation cascades (for example, MAPK cascades). For any system of this class we explicitly exhibit a single species for each connected component of the associated digraph such that the successive total derivatives of its concentration allow us to identify all the parameters occurring in the component. The number of derivatives needed is bounded essentially by the length of the corresponding connected component of the digraph. Moreover, in the particular case of the cascades, we show that the parameters can be identified from a bounded number of successive derivatives of the last product of the last layer. This theoretical result induces also a heuristic interpolation-based identifiability procedure to recover the values of the rate constants from exact measurements.


title: Black-Box Autoregressive Density Estimation for State-Space Models

abstract: State-space models (SSMs) provide a flexible framework for modelling time-series data. Consequently, SSMs are ubiquitously applied in areas such as engineering, econometrics and epidemiology. In this paper we provide a fast approach for approximate Bayesian inference in SSMs using the tools of deep learning and variational inference.


title: Causal Inference by String Diagram Surgery

abstract: Extracting causal relationships from observed correlations is a growing area in probabilistic reasoning, originating with the seminal work of Pearl and others from the early 1990s. This paper develops a new, categorically oriented view based on a clear distinction between syntax (string diagrams) and semantics (stochastic matrices), connected via interpretations as structure-preserving functors. A key notion in the identification of causal effects is that of an intervention, whereby a variable is forcefully set to a particular value independent of any prior propensities. We represent the effect of such an intervention as an endofunctor which performs `string diagram surgery' within the syntactic category of string diagrams. This diagram surgery in turn yields a new, interventional distribution via the interpretation functor. While in general there is no way to compute interventional distributions purely from observed data, we show that this is possible in certain special cases using a calculational tool called comb disintegration. We demonstrate the use of this technique on a well-known toy example, where we predict the causal effect of smoking on cancer in the presence of a confounding common cause. After developing this specific example, we show this technique provides simple sufficient conditions for computing interventions which apply to a wide variety of situations considered in the causal inference literature.


title: Theoretically motivated search and detection of non-thermal pulsations from PSRs J1747-2958, J2021+3651, and J1826-1256

abstract: Based on a theoretical selection of pulsars as candidates for detection at X-ray energies, we present an analysis of archival X-ray observations done with Chandra and XMM-Newton of PSR J1747-2958 (the pulsar in the "Mouse" nebula), PSR J2021+3651 (the pulsar in the "Dragonfly" nebula), and PSR J1826-1256. X-ray pulsations from PSR J1747-2958 and PSR J1826-1256 are detected for the first time, and a previously reported hint of an X-ray pulsation from PSR J2021+3651 is confirmed with a higher significance. We analyze these pulsars' spectra in regards to the theoretically predicted energy distribution, finding a remarkable agreement, and provide here a refined calculation of the model parameters taking into account the newly derived X-ray spectral data.


title: On the eigenvalues of truncations of random unitary matrices

abstract: We consider the empirical eigenvalue distribution of an $m\times m$ principle submatrix of an $n\times n$ random unitary matrix distributed according to Haar measure. Earlier work of Petz and R\'effy identified the limiting spectral measure if $\frac{m}{n}\to\alpha$, as $n\to\infty$; under suitable scaling, the family $\{\mu_\alpha\}_{\alpha\in(0,1)}$ of limiting measures interpolates between uniform measure on the unit disc (for small $\alpha$) and uniform measure on the unit circle (as $\alpha\to1$). In this note, we prove an explicit concentration inequality which shows that for fixed $n$ and $m$, the bounded-Lipschitz distance between the empirical spectral measure and the corresponding $\mu_\alpha$ is typically of order $\sqrt{\frac{\log(m)}{m}}$ or smaller. The approach is via the theory of two-dimensional Coulomb gases and makes use of a new "Coulomb transport inequality" due to Chafa\"i, Hardy, and Ma\"ida.