title: Cosmological Tests With Strong Gravitational Lenses using Gaussian Processes

abstract: Strong gravitational lenses provide source/lens distance ratios D_obs useful in cosmological tests. Previously, a catalog of 69 such systems was used in a one-on-one comparison between the standard model, LCDM, and the R_h=ct universe, which has thus far been favored by the application of model selection tools to many other kinds of data. But in that work, the use of model parametric fits to the observations could not easily distinguish between these two cosmologies, in part due to the limited measurement precision. Here, we instead use recently developed methods based on Gaussian Processes (GP), in which D_obs may be reconstructed directly from the data without assuming any parametric form. This approach not only smooths out the reconstructed function representing the data, but also reduces the size of the 1-sigma confidence regions, thereby providing greater power to discern between different models. With the current sample size, we show that analyzing strong lenses with a GP approach can definitely improve the model comparisons, producing probability differences in the range ~10-30%. These results are still marginal, however, given the relatively small sample. Nonetheless, we conclude that the probability of R_h=ct being the correct cosmology is somewhat higher than that of LCDM, with a degree of significance that grows with the number of sources in the subsamples we consider. Future surveys will significantly grow the catalog of strong lenses and will therefore benefit considerably from the GP method we describe here. In addition, we point out that if the R_h=ct universe is eventually shown to be the correct cosmology, the lack of free parameters in the study of strong lenses should provide a remarkably powerful tool for uncovering the mass structure in lensing galaxies.

categories: astro-ph.CO astro-ph.GA astro-ph.HE gr-qc hep-ph


title: Confounder Detection in High Dimensional Linear Models using First Moments of Spectral Measures

abstract: In this paper, we study the confounder detection problem in the linear model, where the target variable $Y$ is predicted using its $n$ potential causes $X_n=(x_1,...,x_n)^T$. Based on an assumption of rotation invariant generating process of the model, recent study shows that the spectral measure induced by the regression coefficient vector with respect to the covariance matrix of $X_n$ is close to a uniform measure in purely causal cases, but it differs from a uniform measure characteristically in the presence of a scalar confounder. Then, analyzing spectral measure pattern could help to detect confounding. In this paper, we propose to use the first moment of the spectral measure for confounder detection. We calculate the first moment of the regression vector induced spectral measure, and compare it with the first moment of a uniform spectral measure, both defined with respect to the covariance matrix of $X_n$. The two moments coincide in non-confounding cases, and differ from each other in the presence of confounding. This statistical causal-confounding asymmetry can be used for confounder detection. Without the need of analyzing the spectral measure pattern, our method does avoid the difficulty of metric choice and multiple parameter optimization. Experiments on synthetic and real data show the performance of this method.


title: A turnover in the radio light curve of GW170817

abstract: We present 2-9 GHz radio observations of GW170817 covering the period 125-200 days post-merger, taken with the Australia Telescope Compact Array and the Karl G. Jansky Very Large Array. Our observations demonstrate that the radio afterglow peaked at $149\pm2$ days post-merger and is now declining in flux density. We see no evidence for evolution in the radio-only spectral index, which remains consistent with optically-thin synchrotron emission connecting the radio, optical, and X-ray regimes. The peak implies a total energy in the synchrotron-emitting component of a ${\rm few}\times 10^{50}\,$erg. The temporal decay rate is most consistent with mildly- or non-relativistic material and we do not see evidence for a very energetic off-axis jet, but we cannot distinguish between a lower-energy jet and more isotropic emission.


title: MONICA in Hamburg: Towards Large-Scale IoT Deployments in a Smart City

abstract: Modern cities and metropolitan areas all over the world face new management challenges in the 21st century primarily due to increasing demands on living standards by the urban population. These challenges range from climate change, pollution, transportation, and citizen engagement, to urban planning, and security threats. The primary goal of a Smart City is to counteract these problems and mitigate their effects by means of modern ICT to improve urban administration and infrastructure. Key ideas are to utilise network communication to inter-connect public authorities; but also to deploy and integrate numerous sensors and actuators throughout the city infrastructure - which is also widely known as the Internet of Things (IoT). Thus, IoT technologies will be an integral part and key enabler to achieve many objectives of the Smart City vision. The contributions of this paper are as follows. We first examine a number of IoT platforms, technologies and network standards that can help to foster a Smart City environment. Second, we introduce the EU project MONICA which aims for demonstration of large-scale IoT deployments at public, inner-city events and give an overview on its IoT platform architecture. And third, we provide a case-study report on SmartCity activities by the City of Hamburg and provide insights on recent (on-going) field tests of a vertically integrated, end-to-end IoT sensor application.


title: Einstein@Home discovers a radio-quiet gamma-ray millisecond pulsar

abstract: Millisecond pulsars (MSPs) are old neutron stars that spin hundreds of times per second and appear to pulsate as their emission beams cross our line of sight. To date, radio pulsations have been detected from all rotation-powered MSPs. In an attempt to discover radio-quiet gamma-ray MSPs, we used the aggregated power from the computers of tens of thousands of volunteers participating in the Einstein@Home distributed computing project to search for pulsations from unidentified gamma-ray sources in Fermi Large Area Telescope data. This survey discovered two isolated MSPs, one of which is the only known rotation-powered MSP to remain undetected in radio observations. These gamma-ray MSPs were discovered in completely blind searches without prior constraints from other observations, raising hopes for detecting MSPs from a predicted Galactic bulge population.


title: Ultrafast Photodynamics of Glucose

abstract: We have investigated the photodynamics of $\beta$-D-glucose employing our field-induced surface hopping method (FISH), which allows us to simulate the coupled electron-nuclear dynamics, including explicitly nonadiabatic effects and light-induced excitation. Our results reveal that from the initially populated S$_{1}$ and S$_{2}$ states, glucose returns nonradiatively to the ground state within about 200 fs. This takes place mainly via conical intersections (CIs) whose geometries in most cases involve the elongation of a single O-H bond, while in some instances ring-opening due to dissociation of a C-O bond is observed. Experimentally, excitation to a distinct excited electronic state is improbable due to the presence of a dense manifold of states bearing similar oscillator strengths. Our FISH simulations explicitly including a UV laser pulse of 6.43 eV photon energy reveals that after initial excitation the population is almost equally spread over several close-lying electronic states. This is followed by a fast nonradiative decay on the time scale of 100-200 fs, with the final return to the ground state proceeding via the S$_{1}$ state through the same types of CIs as observed in the field-free simulations.


title: Anomalous heat equation in a system connected to thermal reservoirs

abstract: We study anomalous transport in a one-dimensional system with two conserved quantities in presence of thermal baths. In this system we derive exact expressions of the temperature profile and the two point correlations in steady state as well as in the non-stationary state where the later describes the relaxation to the steady state. In contrast to the Fourier heat equation in the diffusive case, here we show that the evolution of the temperature profile is governed by a non-local anomalous heat equation. We provide numerical verifications of our results.

categories: cond-mat.stat-mech math-ph math.MP nlin.CD


title: An improved isomorphism test for bounded-tree-width graphs

abstract: We give a new fpt algorithm testing isomorphism of $n$-vertex graphs of tree width $k$ in time $2^{k\operatorname{polylog} (k)}\operatorname{poly} (n)$, improving the fpt algorithm due to Lokshtanov, Pilipczuk, Pilipczuk, and Saurabh (FOCS 2014), which runs in time $2^{\mathcal{O}(k^5\log k)}\operatorname{poly} (n)$. Based on an improved version of the isomorphism-invariant graph decomposition technique introduced by Lokshtanov et al., we prove restrictions on the structure of the automorphism groups of graphs of tree width $k$. Our algorithm then makes heavy use of the group theoretic techniques introduced by Luks (JCSS 1982) in his isomorphism test for bounded degree graphs and Babai (STOC 2016) in his quasipolynomial isomorphism test. In fact, we even use Babai's algorithm as a black box in one place. We also give a second algorithm which, at the price of a slightly worse running time $2^{\mathcal{O}(k^2 \log k)}\operatorname{poly} (n)$, avoids the use of Babai's algorithm and, more importantly, has the additional benefit that it can also used as a canonization algorithm.


title: Approximation of non-archimedean Lyapunov exponents and applications over global fields

abstract: Let $K$ be an algebraically closed field of characteristic 0 that is complete with respect to a non-archimedean absolute value. We establish a locally uniform approximation formula of the Lyapunov exponent of a rational map $f$ of $\mathbb{P}^1$ of degree $d>1$ over $K$, in terms of the multipliers of $n$-periodic points of $f$, with an explicit control in terms of $n$, $f$ and $K$. As an immediate consequence, we obtain an estimate for the blow-up of the Lyapunov exponent near a pole in one-dimensional families of rational maps over $K$. Combined with our former archimedean version, this non-archimedean quantitative approximation allows us to show: - a quantified version of Silverman's and Ingram's recent comparison between the critical height and any ample height on the moduli space $\mathcal{M}_d(\bar{\mathbb{Q}})$, - two improvements of McMullen's finiteness of the multiplier maps: reduction to multipliers of cycles of exact given period and an effective bound from below on the period, - a characterization of non-affine isotrivial rational maps defined over the function field $\mathbb{C}(X)$ of a normal projective variety $X$ in terms of the growth of the degree of the multipliers.


title: Combined readout of a triple-GEM detector

abstract: Optical readout of GEM based devices by means of high granularity and low noise CMOS sensors allows to obtain very interesting tracking performance. Space resolution of the order of tens of $\mu$m were measured on the GEM plane along with an energy resolution of 20%$\div$30%. The main limitation of CMOS sensors is represented by their poor information about time structure of the event. In this paper, the use of a concurrent light readout by means of a suitable photomultiplier and the acquisition of the electric signal induced on the GEM electrode are exploited to provide the necessary timing informations. The analysis of the PMT waveform allows a 3D reconstruction of each single clusters with a resolution on z of 100 $\mu$m. Moreover, from the PMT signals it is possible to obtain a fast reconstruction of the energy released within the detector with a resolution of the order of 25% even in the tens of keV range useful, for example, for triggering purpose.


title: Analytical and numerical study of the non-linear noisy voter model on complex networks

abstract: We study the noisy voter model using a specific non-linear dependence of the rates that takes into account collective interaction between individuals. The resulting model is solved exactly under the all-to-all coupling configuration and approximately in some random networks environments. In the all-to-all setup we find that the non-linear interactions induce "bona fide" phase transitions that, contrary to the linear version of the model, survive in the thermodynamic limit. The main effect of the complex network is to shift the transition lines and modify the finite-size dependence, a modification that can be captured with the introduction of an effective system size that decreases with the degree heterogeneity of the network. While a non-trivial finite-size dependence of the moments of the probability distribution is derived from our treatment, mean-field exponents are nevertheless obtained in the thermodynamic limit. These theoretical predictions are well confirmed by numerical simulations of the stochastic process.


title: Reconstruction of cosmic ray air showers with Tunka-Rex data using template fitting of radio pulses

abstract: We present an improved method for the precise reconstruction of cosmic ray air showers above $10^{17}$ eV with sparse radio arrays. The method is based on the comparison of predictions for radio pulse shapes by CoREAS simulations to measured pulses. We applied our method to the data of Tunka-Rex, a 1 km$^2$ radio array in Siberia operating in the frequency band of 30-80 MHz. Tunka-Rex is triggered by the air-Cherenkov detector Tunka-133 and by scintillators (Tunka-Grande). The instrument collects air-shower data since 2012. The present paper describes updated data and signal analyses of Tunka-Rex and details of a new method applied. After efficiency cuts, when Tunka-Rex reaches its full efficiency, the energy resolution of about 10% given by the new method has reached the limit of systematic uncertainties due to the calibration uncertainty and shower-to-shower fluctuations. At the same time the shower maximum reconstruction is significantly improved up to an accuracy of 35 g/cm$^2$ compared to the previous method based on the slope of the lateral distribution. We also define and now achieved conditions of the measurements, at which the shower maximum resolution of Tunka-Rex reaches a value of 25 g/cm$^2$ and becomes competitive to optical detectors. To check and validate our reconstruction and efficiency cuts we compare individual events to the reconstruction of Tunka-133. Furthermore, we compare the mean of shower maximum as a function of primary energy to the measurements of other experiments.


title: Approximation property on entropies for surface diffeomorphisms

abstract: In this paper, we prove that for any $C^1$ surface diffeomorphism $f$ with positive topological entropy, there exists a diffeomorphism $g$ arbitrarily close (in the $C^1$ topology) to $f$ exhibiting a horseshoe $\Lambda$, such that the topological entropy of $g$ restricted on $\Lambda$ can arbitrarily approximate the topological entropy of $f$. This extends the Theorem \cite[Theorem 1.1]{Gan} of Gan.


title: On the hierarchical structure of Pareto critical sets

abstract: In this article we show that the boundary of the Pareto critical set of an unconstrained multiobjective optimization problem (MOP) consists of Pareto critical points of subproblems considering subsets of the objective functions. If the Pareto critical set is completely described by its boundary (e.g. if we have more objective functions than dimensions in the parameter space), this can be used to solve the MOP by solving a number of MOPs with fewer objective functions. If this is not the case, the results can still give insight into the structure of the Pareto critical set. This technique is especially useful for efficiently solving many-objective optimization problems by breaking them down into MOPs with a reduced number of objective functions.


title: Spatio-temporal Poisson processes for visits to small sets

abstract: For many measure preserving dynamical systems $(\Omega,T,m)$ the successive hitting times to a small set is well approximated by a Poisson process on the real line. In this work we define a new process obtained from recording not only the successive times $n$ of visits to a set $A$, but also the position $T^n(x)$ in $A$ of the orbit, in the limit where $m(A)\to0$. We obtain a convergence of this process, suitably normalized, to a Poisson point process in time and space under some decorrelation condition. We present several new applications to hyperbolic maps and SRB measures, including the case of a neighborhood of a periodic point, and some billiards such as Sinai billiards, Bunimovich stadium and diamond billiard.


title: The planar 3-body problem II:reduction to pure shape and spherical geometry (2nd version)

abstract: Geometric reduction of the Newtonian planar three-body problem is investigated in the framework of equivariant Riemannian geometry, which reduces the study of trajectories of three-body motions to the study of their moduli curves, that is, curves which record the change of size and shape, in the moduli space of oriented mass-triangles. The latter space is a Riemannian cone over the shape 2-sphere, and the shape curve is the image curve on this sphere. It is shown that the time parametrized moduli curve is in general determined by the relative geometry of the shape curve and the shape potential function. This also entails the reconstruction of time, namely the geometric shape curve determines the time parametrization of the moduli curve, hence also the three-body motion itself, modulo a fixed rotation of the plane. The first version of this work is an (unpublished) paper from 2012, and the present version is an editorial revision of this.


title: Cloud Infrastructure Provenance Collection and Management to Reproduce Scientific Workflow Execution

abstract: The emergence of Cloud computing provides a new computing paradigm for scientific workflow execution. It provides dynamic, on-demand and scalable resources that enable the processing of complex workflow-based experiments. With the ever growing size of the experimental data and increasingly complex processing workflows, the need for reproducibility has also become essential. Provenance has been thought of a mechanism to verify a workflow and to provide workflow reproducibility. One of the obstacles in reproducing an experiment execution is the lack of information about the execution infrastructure in the collected provenance. This information becomes critical in the context of Cloud in which resources are provisioned on-demand and by specifying resource configurations. Therefore, a mechanism is required that enables capturing of infrastructure information along with the provenance of workflows executing on the Cloud to facilitate the re-creation of execution environment on the Cloud. This paper presents a framework, ReCAP, along with the proposed mapping approaches that aid in capturing the Cloud-aware provenance information and help in re-provisioning the execution resource on the Cloud with similar configurations. Experimental evaluation has shown the impact of different resource configurations on the workflow execution performance, therefore justifies the need for collecting such provenance information in the context of Cloud. The evaluation has also demonstrated that the proposed mapping approaches can capture Cloud information in various Cloud usage scenarios without causing performance overhead and can also enable the re-provisioning of resources on Cloud. Experiments were conducted using workflows from different scientific domains such as astronomy and neuroscience to demonstrate the applicability of this research for different workflows.


title: Generation of ten kilotesla longitudinal magnetic fields in ultraintense laser-solenoid target interactions

abstract: Production of the huge longitudinal magnetic fields by using an ultraintense laser pulse irradiating a solenoid target is considered. Through three-dimensional particle-in-cell simulations, it is shown that the longitudinal magnetic field up to ten kilotesla can be observed in the ultraintense laser-solenoid target interactions. The finding is associated with both fast and return electron currents in the solenoid target. The huge longitudinal magnetic field is of interest for a number of important applications, which include controlling the divergence of laser-driven energetic particles for medical treatment, fast-ignition in inertial fusion, etc., as an example, the well focused and confined directional electron beams are realized by using the solenoid target.


title: Realistic Floquet semimetal with exotic topological linkages between arbitrarily many nodal loops

abstract: Valence and conduction bands in nodal loop semimetals (NLSMs) touch along closed loops in momentum space. If such loops can proliferate and link intricately, NLSMs become exotic topological phases with unconventional topological characteristics and potentially peculiar transport properties. In conventional quantum materials or cold atom systems alike, such exotic phases necessarily require non-local hopping and are therefore intrinsically unrealistic. In this work, we show how this hurdle can be surmounted through an experimentally feasible periodic driving scheme. In particular, by tuning the period of a two-step periodic driving or some experimentally accessible parameters, we show how to generate arbitrarily many nodal loops that are linked with various levels of complexity. Furthermore, we propose to use both the Berry phase winding and the Alexander polynomial topological invariant to characterize the fascinating linkages among the nodal loops. This work thus presents a class of exotic Floquet topological phase that has hitherto not been proposed in any realistic setup.


title: An invitation to higher Teichm\"uller theory

abstract: The goal of this article is to invite the reader to get to know and to get involved into higher Teichm\"uller theory by describing some of its many facets.


title: Symbol-Level Precoding Design for Max-Min SINR in Multiuser MISO Broadcast Channels

abstract: In this paper, we address the symbol level precoding (SLP) design problem under max-min SINR criterion in the downlink of multiuser multiple-input single-output (MISO) channels. First, we show that the distance preserving constructive interference regions (DPCIR) are always polyhedral angles (shifted pointed cones) for any given constellation point with unbounded decision region. Then we prove that any signal in a given unbounded DPCIR has a norm larger than the norm of the corresponding vertex if and only if the convex hull of the constellation contains the origin. Using these properties, we show that the power of the noiseless received signal lying on an unbounded DPCIR is an strictly increasing function of two parameters. This allows us to reformulate the originally non-convex SLP max-min SINR as a convex optimization problem. We discuss the loss due to our proposed convex reformulation and provide some simulation results.


title: The group generated by Riordan involutions

abstract: We prove that any element in the group generated by the Riordan involutions is the product of at most four of them. We also give a description of this subgroup as a semidirect product of a special subgroup of the commutator subgroup and the Klein four-group.


title: The eBDIMS path-sampling server: generation, classification and interactive visualization of protein ensembles and transition pathways in 2D-motion space

abstract: The recent rise of cryo-EM and X-ray high-throughput techniques is providing a wealth of new structures trapped in different conformations. Understanding how proteins transition between different conformers, and how they relate to each other in terms of function is not straightforward, and highly depends on the choice of the right set of degrees of freedom. Here we present eBDIMS server, an online tool and software for automatic classification of structural ensembles and reconstruction of transition pathways using coarse-grained (CG) simulations. The server generates CG-pathways between two protein conformations along with a representation in a simplified 2D-motion landscape based on the Principal Components (PCs) from experimental structures. For a conformationally rich ensemble, the PCs provide powerful reaction coordinates for automatic structure classification, detection of on-pathway intermediates and validation of in silico pathways. When the number of available structures is low or sampling is limited, Normal Modes (NMs) provide alternative motion axes for trajectory analysis. The path-generation eBDIMS method is available at a user-friendly website: https://login.biophysics.kth.se/eBDIMS/ or as standalone software. The server incorporates a powerful interactive graphical interface for simultaneous visualization of transition pathways in 2D-motion space and 3D-molecular graphics, which greatly facilitates the exploration of the relationships between different conformations.


title: Multi-Codec DASH Dataset

abstract: The number of bandwidth-hungry applications and services is constantly growing. HTTP adaptive streaming of audio-visual content accounts for the majority of today's internet traffic. Although the internet bandwidth increases also constantly, audio-visual compression technology is inevitable and we are currently facing the challenge to be confronted with multiple video codecs. This paper proposes a multi-codec DASH dataset comprising AVC, HEVC, VP9, and AV1 in order to enable interoperability testing and streaming experiments for the efficient usage of these codecs under various conditions. We adopt state of the art encoding and packaging options and also provide basic quality metrics along with the DASH segments. Additionally, we briefly introduce a multi-codec DASH scheme and possible usage scenarios. Finally, we provide a preliminary evaluation of the encoding efficiency in the context of HTTP adaptive streaming services and applications.


title: On the streaming complexity of fundamental geometric problems

abstract: In this paper, we focus on lower bounds and algorithms for some basic geometric problems in the one-pass (insertion only) streaming model. The problems considered are grouped into three categories: (i) Klee's measure (ii) Convex body approximation, geometric query, and (iii) Discrepancy Klee's measure is the problem of finding the area of the union of hyperrectangles. Under convex body approximation, we consider the problems of convex hull, convex body approximation, linear programming in fixed dimensions. The results for convex body approximation implies a property testing type result to find if a query point lies inside a convex polyhedron. Under discrepancy, we consider both the geometric and combinatorial discrepancy. For all the problems considered, we present (randomized) lower bounds on space. Most of our lower bounds are in terms of approximating the solution with respect to an error parameter $\epsilon$. We provide approximation algorithms that closely match the lower bound on space for most of the problems.


title: Generalised Net Convergence Structures in Posets

abstract: In this paper, we introduce the notion of $\mathcal{M}$-convergence and $\mathcal{MN}$-convergence structures in posets, which, in some sense, generalise the well-known Scott-convergence and order-convergence structures. As results, we give a necessary and sufficient conditions for each generalised convergence structures being topological. These results then imply the following two well-established results: (1) The Scott-convergence structure in a poset $P$ is topological if and only if $P$ is continuous, and (2) The order-convergence structure in a poset $P$ is topological if and only if $P$ is $\mathcal{R}^*$-doubly continuous.


title: The mm-wave compact component of AGN

abstract: mm-wave emission from Active Galactic Nuclei (AGN) may hold the key to understanding the physical origin of their radio cores. The correlation between radio/mm and X-ray luminosity may suggest a similar physical origin of the two sources. Since synchrotron self absorption decreases with frequency, mm-waves probe smaller length scales than cm-waves. We report on 100 GHz (3 mm) observations with CARMA of 26 AGNs selected from the hard X-ray Swift/BAT survey. 20/26 targets were detected at 100 GHz down to the 1 mJy (3 $\sigma$) sensitivity, which corresponds to optically thick synchrotron source sizes of 10$^{-4}$ - 10$^{-3}$ pc). Most sources show a 100 GHz flux excess with respect to the spectral slope extrapolated from low frequencies. This mm spectral component likely originates from smaller scales than the few-GHz emission. The measured mm sources lie roughly around the L$_{mm}$ (100 GHz) $\sim$ 10$^{-1}$ L$_{X}$ (2-10 keV) relation, similar to a few previously published X-ray selected sources, and hinting perhaps at a common coronal origin.


title: Parameterized Complexity of Fair Vertex Evaluation Problems

abstract: A prototypical graph problem is centered around a graph-theoretic property for a set of vertices and a solution to it is a set of vertices for which the desired property holds. The task is to decide whether, in the given graph, there exists a solution of a certain quality, where we use size as a quality measure. In this work, we are changing the measure to the fair measure [Lin&Sahni: Fair edge deletion problems. IEEE Trans. Comput. 89]. The measure is k if the number of solution neighbors does not exceed k for any vertex in the graph. One possible way to study graph problems is by defining the property in a certain logic. For a given objective an evaluation problem is to find a set (of vertices) that simultaneously minimizes the assumed measure and satisfies an appropriate formula. In the presented paper we show that there is an FPT algorithm for the MSO Fair Vertex Evaluation problem for formulas with one free variable parameterized by the twin cover number of the input graph. Here, the free variable corresponds to the solution sought. One may define an extended variant of MSO Fair Vertex Evaluation for formulas with l free variables; here we measure a maximum number of neighbors in each of the l sets. However, such variant is W[1]-hard for parameter l even on graphs with twin cover one. Furthermore, we study the Fair Vertex Cover (Fair VC) problem. Fair VC is among the simplest problems with respect to the demanded property (i.e., the rest forms an edgeless graph). On the negative side, Fair VC is W[1]-hard when parameterized by both treedepth and feedback vertex set of the input graph. On the positive side, we provide an FPT algorithm for the parameter modular width.


title: The tree of numerical semigroups with low multiplicity

abstract: We show that the number of numerical semigroups with multiplicity three, four or five and fixed genus is increasing as a function in the genus. To this end we use the Kunz polytope for these multiplicities. Counting numerical semigroups with fixed multiplicity and genus is then an integer partition problem with some extra conditions (those of membership to the Kunz polytope). For the particular case of multiplicity four, we are able to prove that the number of numerical semigroups with multiplicity four and genus $g$ is the number of partitions $x+y+z=g+6$ with $0<x\le y\le z$, $x\neq 1$, $y\neq 2$ and $z\neq 3$.


title: Scaling of energy amplification in viscoelastic channel and Couette flow

abstract: The linear amplification of disturbances is critical in setting up transition scenarios in viscoelastic channel and Couette flow, and may also play an important role when such flows are fully turbulent. As such, it is of interest to assess how this amplification, defined as the steady-state variance maintained under Gaussian white noise forcing, scales with the main nondimensional parameters: the Reynolds ($Re$) and Weissenberg ($Wi$) numbers. This scaling is derived analytically in the two limits of strong and weak elasticity for when the forcing is streamwise-constant. The latter is the relevant forcing for capturing the overall behaviour because it was previously shown to have the dominant contribution to amplification. The final expressions show that for weak elasticity the scaling retains a form similar to the well-known O($Re^3$) relationship with an added elastic correction. For strong elasticity, however, the scaling is O($Wi^3$) with a viscous correction. The key factor leading to such a mirroring in the scaling is the introduction of forcing in the polymer stress. The results demonstrate that energy amplification in a viscoelastic flow can be very sensitive to the model parameters even at low $Re$. They also suggest that energy amplification can be significantly increased by forcing the polymer stress, thereby opening up possibilities such as flow control using systematically designed polymer stress perturbations.


title: Convex resource theory of non-Markovianity

abstract: We establish a convex resource theory of non-Markovianity under the constraint of small time intervals within the temporal evolution. We construct the free operations, free states and a generalized bona-fide measure of non-Markovianity. The framework satisfies the basic properties of a consistent resource theory. The proposed resource quantifier is lower bounded by the optimization free Rivas-Huelga-Plenio (RHP) measure of nonMarkovianity. We further define the robustness of non-Markovianity and show that it can directly be expressed as a function of the RHP measure of non-Markovianity. This enables a physical interpretation of the RHP measure.


title: The correct formulation of Gleason's theorem in quaternionic Hilbert spaces

abstract: From the viewpoint of the theory of orthomodular lattices of elementary propositions, Quantum Theories can be formulated in real, complex or quaternionic Hilbert spaces as established in Sol\'er's theorem. The said lattice eventually coincides with the lattice of all orthogonal projectors on a separable Hilbert space over R, C, or over the algebra of quaternions H. Quantum states are $\sigma$-additive probability measures on that non-Boolean lattice. Gleason's theorem proves that, if the Hilbert space is separable with dimension >2 and the Hilbert space is either real or complex, then states are one-to-one with standard density matrices (self-adjoint, positive, unit-trace, trace-class operators). The extension of this result to quaternionic Hilbert spaces was obtained by Varadarajan in 1968. Unfortunately, even if the hard part of the proof is correct, the formulation of this extension is mathematically incorrect. This is due to some peculiarities of the notion of trace in quaternionic Hilbert spaces, e.g., basis dependence, making the theory of trace-class operators in quaternionic Hilbert spaces different from the standard theory in real and complex Hilbert spaces. A minor issue also affects Varadarajan's statement for real Hilbert space formulation. This paper is mainly devoted to present Gleason-Varadarajan's theorem into a technically correct form valid for the three types of Hilbert spaces. After having develped part of the general mathematical technology of trace-class operators in (generally non-separable) quaternionic Hilbert spaces, we prove that only the {\em real part} of the trace enters the formalism of quantum theories (also dealing with unbounded observables and symmetries) and it can be safely used to formulate and prove a common statement of Gleason's theorem.

categories: math-ph hep-th math.CV math.MP math.OA quant-ph


title: Tetraquarks in the 1/N Expansion: a New Appraisal

abstract: We discuss the necessary, albeit not sufficient, conditions for tetraquark poles to occur in the 1/N expansion of QCD and find the minimum order at which such poles may appear. Assuming tetraquark poles, we find a new non-planar solution with the minimal number of topologies and tetraquark species. The solution implies narrow states. Mixing with quarkonium states is allowed so that P-wave tetraquarks with J^PC=1^-- would couple to e^+e^-.


title: Trade-off Between Work and Correlations in Quantum Thermodynamics

abstract: Quantum thermodynamics and quantum information are two frameworks for employing quantum mechanical systems for practical tasks, exploiting genuine quantum features to obtain advantages with respect to classical implementations. While appearing disconnected at first, the main resources of these frameworks, work and correlations, have a complicated yet interesting relationship that we examine here. We review the role of correlations in quantum thermodynamics, with a particular focus on the conversion of work into correlations. We provide new insights into the fundamental work cost of correlations and the existence of optimally correlating unitaries, and discuss relevant open problems.


title: Analysis of Angular Observables of $\Lambda_b \to \Lambda (\to p\pi)\mu^{+}\mu^{-}$ Decay in Standard and $Z^{\prime}$ Models

abstract: In 2015, the LHCb collaboration has measured $\frac{d{\mathcal{B}}}{dq^2}$, the lepton- and hadron-side forward-backward asymmetries, denoted by $A^\ell_{FB}$ and $A^{\Lambda}_{FB}$, respectively in the range $15 < q^2(=s) < 20$ GeV$^2$ with 3 fb$^{-1}$ of data. Motivated by these measurements, we perform an analysis of $q^2$ dependent $\Lambda_b \to \Lambda (\to p \pi ) \mu^+\mu^-$ angular observables at large- and low-recoil in the SM and in a family non-universal $Z^{\prime}$ model. In the present study we use the recently performed high-precision lattice QCD calculations of the form factors that have well controlled uncertainties especially in $15 < s < 20$ GeV$^2$ bin. Using the full four-folded angular distribution of $\Lambda_b \to \Lambda (\to p \pi ) \mu^+\mu^-$ decay, firstly we calculate the values of these measured quantitites in the SM and compare their numerical values with the measurements in appropriate bins of $s$. In case of the possible discrepancy between the SM prediction and measurements, we try to see if these can be accommodated though the extra neutral $Z^{\prime}$ boson. In addition, the fraction of longitudinal polarization of the dimuon $F_{L}$ is measured to be $0.61^{+0.11}_{-0.14}\pm 0.03$ in $15 < s < 20$ GeV$^2$ at the LHCb. We find that in this bin the value found in the $Z^{\prime}$ model is close to the observed values. After comparing the results of these observables, we have proposed a number of other observables whose values are calculated in different bins of $s$ in the SM and $Z^{\prime}$ model. We illustrate that the experimental observations of these observables in several bins of $s$ can help to test the predictions of the SM and unravel NP contributions arises due to $Z^{\prime}$ model in these decays.


title: Exchanging role of the phase space and symmetry group of integrable Hamiltonian systems related to Lie bialgebras of bi-symplectic types

abstract: We construct integrable Hamiltonian systems with Lie bialgebras $({\bf g} , {\bf \tilde{g}})$ of the bi-symplectic type for which the Poisson-Lie groups ${\bf G}$ play the role of the phase spaces, and their dual Lie groups ${\bf {\tilde {G}}}$ play the role of the symmetry groups of the systems. We give the new transformations to exchange the role of phase spaces and symmetry groups and obtain the relations between integrals of motions of these integrable systems. Finally, we give some examples of real four-dimensional Lie bialgebras of bi-symplectic type.


title: Lossless Analog Compression

abstract: We establish the fundamental limits of lossless analog compression by considering the recovery of arbitrary m-dimensional real random vectors x from the noiseless linear measurements y=Ax with n x m measurement matrix A. Our theory is inspired by the groundbreaking work of Wu and Verdu (2010) on almost lossless analog compression, but applies to the nonasymptotic, i.e., fixed-m case, and considers zero error probability. Specifically, our achievability result states that, for almost all A, the random vector x can be recovered with zero error probability provided that n > K(x), where K(x) is given by the infimum of the lower modified Minkowski dimension over all support sets U of x. We then particularize this achievability result to the class of s-rectifiable random vectors as introduced in Koliander et al. (2016); these are random vectors of absolutely continuous distribution---with respect to the s-dimensional Hausdorff measure---supported on countable unions of s-dimensional differentiable submanifolds of the m-dimensional real coordinate space. Countable unions of differentiable submanifolds include essentially all signal models used in the compressed sensing literature. Specifically, we prove that, for almost all A, s-rectifiable random vectors x can be recovered with zero error probability from n>s linear measurements. This threshold is, however, found not to be tight as exemplified by the construction of an s-rectifiable random vector that can be recovered with zero error probability from n<s linear measurements. This leads us to the introduction of the new class of s-analytic random vectors, which admit a strong converse in the sense of n greater than or equal to s being necessary for recovery with probability of error smaller than one. The central conceptual tools in the development of our theory are geometric measure theory and the theory of real analytic functions.


title: Gamma-Ray Observations of Nova Sgr 2015 No. 2 with INTEGRAL

abstract: INTEGRAL observed the nova V5668 Sgr around the time of its optical maximum on March 21, 2015. Studies at UV wavelengths showed spectral lines of freshly produced Be-7. This could be measurable also in gamma-rays at 478 keV from the decay to Li-7. Novae are also expected to synthesise Na-22 which decays to Ne-22, emitting a 1275 keV photon. About one week before the optical maximum, a strong gamma-ray flash on time-scales of hours is expected from short-lived radioactive nuclei, such as N-13 and F-18. These beta-plus-unstable nuclei should yield emission up to 511 keV, but which has never been observed. The spectrometer SPI aboard INTEGRAL pointed towards V5668 by chance. We use these observations to search for possible gamma-ray emission of decaying Be-7, and to directly measure the synthesised mass during explosive burning. We also aim to constrain possible burst-like emission days to weeks before the optical maximum using the SPI anticoincidence shield (ACS). We extract spectral and temporal information to determine the fluxes of gamma-ray lines at 478 keV, 511 keV, and 1275 keV. A measured flux value directly converts into abundances produced by the nova. The SPI-ACS rates are analysed for burst-like emission using a nova model light-curve. For the obtained nova flash candidate events, we discuss possible origins. No significant excess for the expected gamma-ray lines is found. Our upper limits on the synthesised Be-7 and Na-22 mass depend on the uncertainties of the distance to the nova: The Be-7 mass is constrained to less than $4.8\times10^{-9}\,(d/kpc)^2$, and Na-22 to less than $2.4\times10^{-8}\,(d/kpc)^2$ solar masses. For the Be-7 mass estimate from UV studies, the distance to V5668 Sgr must be larger than 1.2 kpc. During three weeks before the optical maximum, we find 23 burst-like events in the ACS rate, of which six could possibly be associated with V5668 Sgr.


title: Terahertz emission from laser-driven gas-plasmas: a plasmonic point of view

abstract: We disclose an unanticipated link between plasmonics and nonlinear frequency down-conversion in laser-induced gas-plasmas. For two-color femtosecond pump pulses, a plasmonic resonance is shown to broaden the terahertz emission spectra significantly. We identify the resonance as a leaky mode, which contributes to the emission spectra whenever electrons are excited along a direction where the plasma size is smaller than the plasma wavelength. As a direct consequence, such resonances can be controlled by changing the polarization properties of elliptically-shaped driving laser pulses. Both, experimental results and 3D Maxwell consistent simulations confirm that a significant terahertz pulse shortening and spectral broadening can be achieved by exploiting the transverse driving laser beam shape as an additional degree of freedom.


title: Frequency and mode identification of {\gamma} Doradus from photometric and spectroscopic observations

abstract: The prototype star for the {\gamma} Doradus class of pulsating variables was studied em- ploying photometric and spectroscopic observations to determine the frequencies and modes of pulsation. The four frequencies found were self-consistent between the obser- vation types and almost identical to those found in previous studies (1.3641 d-1 ,1.8783 d-1 , 1.4742 d-1 and 1.3209 d-1). Three of the frequencies are classified as l, m = (1, 1) pulsations and the other is ambiguous between l = 2 modes. Two frequencies are shown to be stable over twenty years since their first identification. The agreement in ground-based work makes this star an excellent calibrator for the upcoming TESS observations and a standard for continued asteroseismic modelling.


title: Correlations between X-ray properties and Black Hole Mass in AGN: towards a new method to estimate black hole mass from short exposure X-ray observations

abstract: Several investigations of the X-ray variability of active galactic nuclei (AGN) using the normalised excess variance (${\sigma^2_{\rm NXS}}$) parameter have shown that variability has a strong anti-correlation with black hole mass ($M_{\rm BH}$) and X-ray luminosity ($L_{\rm X}$). In this study we confirm these previous correlations and find no evidence of a redshift evolution. Using observations from XMM-Newton, we determine the ${\sigma^2_{\rm NXS}}$ and $L_{\rm X}$ for a sample of 1091 AGN drawn from the XMM-Newton Cluster Survey (XCS) - making this the largest study of X-ray spectral properties of AGNs. We created light-curves in three time-scales; 10 ks, 20 ks and 40 ks and used these to derive scaling relations between ${\sigma^2_{\rm NXS}}$, $L_{\rm X}$ (2.0-10 keV range) and literature estimates of $M_{\rm BH}$ from reverberation mapping. We confirm the anti-correlation between $M_{\rm BH}$ and ${\sigma^2_{\rm NXS}}$ and find a positive correlation between $M_{\rm BH}$ and $L_{\rm X}$. The use of ${\sigma^2_{\rm NXS}}$ is practical only for pointed observations where the observation time is tens of kiloseconds. For much shorter observations one cannot accurately quantify variability to estimate $M_{\rm BH}$. Here we describe a method to derive $L_{\rm X}$ from short duration observations and used these results as an estimate for $M_{\rm BH}$. We find that it is possible to estimate $L_{\rm X}$ from observations of just a few hundred seconds and that when correlated with $M_{\rm BH}$, the relation is statistically similar to the relation of $M_{\rm BH}$-$L_{\rm X}$ derived from a spectroscopic analysis of full XMM observations. This method may be particularly useful to the eROSITA mission, an all-sky survey, which will detect $>$10$^{6}$ AGN.


title: Dynamics and Stability of Meshed Multiterminal HVDC Networks

abstract: This paper investigates the existence of an equilibrium point in multiterminal HVDC (MT-HVDC) grids, assesses its uniqueness and defines conditions to ensure its stability. An offshore MT-HVDC system including two wind farms is selected as application test case. At first, a generalized dynamic model of the network is proposed, using hypergraph theory. Such model captures the frequency dependence of transmission lines and cables, it is non-linear due to the constant power behavior of the converter terminals using droop regulation, and presents a suitable degree of simplifications of the MMC converters, under given conditions, to allow system level studies over potentially large networks. Based on this model, the existence and uniqueness of the equilibrium point is demonstrated by returning the analysis to a load-flow problem and using the Banach fixed point theorem. Additionally, the stability of the equilibrium is analyzed by obtaining a Lyapunov function by the Krasovskii's theorem. Computational results obtained for the selected 4 terminals MT-HVDC grid corroborate the requirement for the existence and stability of the equilibrium point.


title: On reference solutions and the sensitivity of the 2D Kelvin-Helmholtz instability problem

abstract: Two-dimensional Kelvin-Helmholtz instability problems are popular examples for assessing discretizations for incompressible flows at high Reynolds number. Unfortunately, the results in the literature differ considerably. This paper presents computational studies of a Kelvin-Helmholtz instability problem with high order divergence-free finite element methods. Reference results in several quantities of interest are obtained for three different Reynolds numbers up to the beginning of the final vortex pairing. A mesh-independent prediction of the final pairing is not achieved due to the sensitivity of the considered problem with respect to small perturbations. A theoretical explanation of this sensitivity to small perturbations is provided based on the theory of self-organization of 2D turbulence. Possible sources of perturbations that arise in almost any numerical simulation are discussed.


title: Main features of detectors and isotopes to investigate double beta decay with increased sensitivity

abstract: The current situation in double beta decay experiments , the characteristics of modern detectors and the possibility of increasing the sensitivity to neutrino mass in future experiments are discussed. The issue of the production and use of enriched isotopes in double beta decay experiments is discussed in addition.


title: Global multiplicity bounds and Spectral Statistics Random Operators

abstract: In this paper, we consider Anderson type operators on a separable Hilbert space where the random perturbations are finite rank and the random variables have full support on $\mathbb{R}$. We show that spectral multiplicity has a uniform lower bound whenever the lower bound is given on a set of positive Lebesgue measure on the point spectrum away from the continuous one. We also show a deep connection between the multiplicity of pure point spectrum and local spectral statistics, in particular, we show that spectral multiplicity higher than one always gives non-Poisson local statistics in the framework of Minami theory. In particular, in higher rank Anderson models with pure-point spectrum, with the randomness having support equal to $\mathbb{R}$, there is a uniform lower bound on spectral multiplicity and in case this is larger than one the local statistics is not Poisson.


title: Prehistory of Transit Searches

abstract: Nowadays the more powerful method to detect extrasolar planets is the transit method. We review the planet transits which were anticipated, searched, and the first ones which were observed all through history. Indeed transits of planets in front of their star were first investigated and studied in the solar system. The first observations of sunspots were sometimes mistaken for transits of unknown planets. The first scientific observation and study of a transit in the solar system was the observation of Mercury transit by Pierre Gassendi in 1631. Because observations of Venus transits could give a way to determine the distance Sun-Earth, transits of Venus were overwhelmingly observed. Some objects which actually do not exist were searched by their hypothetical transits on the Sun, as some examples a Venus satellite and an infra-mercurial planet. We evoke the possibly first use of the hypothesis of an exoplanet transit to explain some periodic variations of the luminosity of a star, namely the star Algol, during the eighteen century. Then we review the predictions of detection of exoplanets by their transits, those predictions being sometimes ancient, and made by astronomers as well as popular science writers. However, these very interesting predictions were never published in peer-reviewed journals specialized in astronomical discoveries and results. A possible transit of the planet beta Pic b was observed in 1981. Shall we see another transit expected for the same planet during 2018? Today, some studies of transits which are connected to hypothetical extraterrestrial civilisations are published in astronomical refereed journals. Some studies which would be classified not long ago as science fiction are now considered as scientific ones.


title: Generic features of the neutron-proton interaction

abstract: We show that fully aligned neutron-proton pairs play a crucial role in the low-energy spectroscopy of nuclei. with valence nucleons in a high-j orbital. Their dominance is valid in nuclei with valence neutrons and protons in different high-j orbitals as well as in N = Z nuclei, where all nucleons occupy the same orbital. We demonstrate analytically this generic feature of the neutron-proton interaction for a variety of systems with four valence nucleons interacting through realistic, effective forces. The dominance of fully aligned neutron-proton pairs results from the combined effect of (i) angular momentum coupling and (ii) basic properties of the neutron-proton interaction.


title: A Mixture of Views Network with Applications to the Classification of Breast Microcalcifications

abstract: In this paper we examine data fusion methods for multi-view data classification. We present a decision concept which explicitly takes into account the input multi-view structure, where for each case there is a different subset of relevant views. This data fusion concept, which we dub Mixture of Views, is implemented by a special purpose neural network architecture. It is demonstrated on the task of classifying breast microcalcifications as benign or malignant based on CC and MLO mammography views. The single view decisions are combined by a data-driven decision, according to the relevance of each view in a given case, into a global decision. The method is evaluated on a large multi-view dataset extracted from the standardized digital database for screening mammography (DDSM). The experimental results show that our method outperforms previously suggested fusion methods.


title: Limit Theorems for Cylindrical Martingale Problems associated with L\'evy Generators

abstract: We prove limit theorems for cylindrical martingale problems associated to L\'evy generators. Furthermore, we give sufficient and necessary conditions for the Feller property of well-posed problems with continuous coefficients. We discuss two applications. First, we derive continuity and linear growth conditions for the existence of weak solutions to infinite-dimensional stochastic differential equations driven by L\'evy noise. Second, we derive continuity, local boundedness and linear growth conditions for limit theorems and the Feller property of weak solutions to stochastic partial differential equations driven by Wiener noise.


title: Analytical Impedance Models for Very Short Bunches

abstract: We discuss several analytical models for impedances of very short bunches. The approximate analytical models are compared with direct solution of Maxwells equations.


title: Cyclic Sieving and Cluster Duality of Grassmannian

abstract: We introduce a decorated configuration space $\mathscr{C}\!{\rm onf}_n^\times(a)$ with a potential function $\mathcal{W}$. We prove the cluster duality conjecture of Fock-Goncharov for Grassmannians, that is, the tropicalization of $\big(\mathscr{C}\!{\rm onf}_n^\times(a), \mathcal{W}\big)$ canonically parametrizes a linear basis of the homogeneous coordinate ring of the Grassmannian $\operatorname{Gr}_a(n)$ with respect to the Pl\"ucker embedding. We prove that $\big(\mathscr{C}\!{\rm onf}_n^\times(a), \mathcal{W}\big)$ is equivalent to the mirror Landau-Ginzburg model of the Grassmannian considered by Eguchi-Hori-Xiong, Marsh-Rietsch and Rietsch-Williams. As an application, we show a cyclic sieving phenomenon involving plane partitions under a sequence of piecewise-linear toggles.


title: Filters for anisotropic wavelet decompositions

abstract: Like the continous shearlet transform and their relatives, discrete transformations based on the interplay between several filterbanks with anisotropic dilations provide a high potential to recover directed features in two and more dimensions. Due to simplicity, most of the directional systems constructed so far were using prediction--correction methods based on interpolatory subdivision schemes. In this paper, we give a simple but effective construction for QMF (quadrature mirror filter) filterbanks which are the discrete object between orthogonal wavelet analysis. We also characterize when the filterbank gives rise to the existence of refinable functions and hence wavelets and give a generalized shearlet construction for arbitrary dimensions and arbitrary scalings for which the filterbank construction ensures the existence of an orthogonal wavelet analysis.


title: On class groups of random number fields

abstract: The main aim of the present paper is to disprove the Cohen--Lenstra--Martinet heuristics in two different ways and to offer possible corrections. We also recast the heuristics in terms of Arakelov class groups, giving an explanation for the probability weights appearing in the general form of the heuristics. We conclude by proposing a rigorously formulated Cohen--Lenstra--Martinet conjecture.


title: Aerial LaneNet: Lane Marking Semantic Segmentation in Aerial Imagery using Wavelet-Enhanced Cost-sensitive Symmetric Fully Convolutional Neural Networks

abstract: The knowledge about the placement and appearance of lane markings is a prerequisite for the creation of maps with high precision, necessary for autonomous driving, infrastructure monitoring, lane-wise traffic management, and urban planning. Lane markings are one of the important components of such maps. Lane markings convey the rules of roads to drivers. While these rules are learned by humans, an autonomous driving vehicle should be taught to learn them to localize itself. Therefore, accurate and reliable lane marking semantic segmentation in the imagery of roads and highways is needed to achieve such goals. We use airborne imagery which can capture a large area in a short period of time by introducing an aerial lane marking dataset. In this work, we propose a Symmetric Fully Convolutional Neural Network enhanced by Wavelet Transform in order to automatically carry out lane marking segmentation in aerial imagery. Due to a heavily unbalanced problem in terms of number of lane marking pixels compared with background pixels, we use a customized loss function as well as a new type of data augmentation step. We achieve a very high accuracy in pixel-wise localization of lane markings without using 3rd-party information. In this work, we introduce the first high-quality dataset used within our experiments which contains a broad range of situations and classes of lane markings representative of current transportation systems. This dataset will be publicly available and hence, it can be used as the benchmark dataset for future algorithms within this domain.


title: TBD: Benchmarking and Analyzing Deep Neural Network Training

abstract: The recent popularity of deep neural networks (DNNs) has generated a lot of research interest in performing DNN-related computation efficiently. However, the primary focus is usually very narrow and limited to (i) inference -- i.e. how to efficiently execute already trained models and (ii) image classification networks as the primary benchmark for evaluation. Our primary goal in this work is to break this myopic view by (i) proposing a new benchmark for DNN training, called TBD (TBD is short for Training Benchmark for DNNs), that uses a representative set of DNN models that cover a wide range of machine learning applications: image classification, machine translation, speech recognition, object detection, adversarial networks, reinforcement learning, and (ii) by performing an extensive performance analysis of training these different applications on three major deep learning frameworks (TensorFlow, MXNet, CNTK) across different hardware configurations (single-GPU, multi-GPU, and multi-machine). TBD currently covers six major application domains and eight different state-of-the-art models. We present a new toolchain for performance analysis for these models that combines the targeted usage of existing performance analysis tools, careful selection of new and existing metrics and methodologies to analyze the results, and utilization of domain specific characteristics of DNN training. We also build a new set of tools for memory profiling in all three major frameworks; much needed tools that can finally shed some light on precisely how much memory is consumed by different data structures (weights, activations, gradients, workspace) in DNN training. By using our tools and methodologies, we make several important observations and recommendations on where the future research and optimization of DNN training should be focused.


title: Understanding dynamics of looping of a long chain polymer in solution: Wilemski-Fixman Approach

abstract: We investigate theoretically the end-to-end looping time of a long chain polymer molecule immersed in a solvent. The dynamics of the end-to-end distance is governed by a Smoluchowski-like equation of a particle moving under the influence of a parabolic potential in presence of a Dirac delta sink of arbitrary strength and location. Using Wilemski-Fixman [ G. Wilemski and M. Fixman, J. Chem. Phys. {\bf 60}, {\it 866} (1974)] approach we calculate the looping time for a long chain molecule immersed in a solvent. We find that looping time varies with several parameters such as length of the polymer (N), bond length (b) and the relaxation time ${\tau_R}$.


title: Auxiliary information : the raking-ratio empirical process

abstract: We study the empirical measure associated to a sample of size $n$ and modified by $N$ iterations of the raking-ratio method. This empirical measure is adjusted to match the true probability of sets in a finite partition which changes each step. We establish asymptotic properties of the raking-ratio empirical process indexed by functions as $n\rightarrow +\infty$, for $N$ fixed. We study nonasymptotic properties by using a Gaussian approximation which yields uniform Berry-Esseen type bounds depending on $n, N$ and provides estimates of the uniform quadratic risk reduction. A closed-form expression of the limiting covariance matrices is derived as $N\rightarrow +\infty$. In the two-way contingency table case the limiting process has a simple explicit formula.


title: Table Based Detection of Degenerate Predicates in Free Space Construction

abstract: The key to a robust and efficient implementation of a computational geometry algorithm is an efficient algorithm for detecting degenerate predicates. We study degeneracy detection in constructing the free space of a polyhedron that rotates around a fixed axis and translates freely relative to another polyhedron. The structure of the free space is determined by the signs of univariate polynomials, called angle polynomials, whose coefficients are polynomials in the coordinates of the vertices of the polyhedra. Every predicate is expressible as the sign of an angle polynomial $f$ evaluated at a zero $t$ of an angle polynomial $g$. A predicate is degenerate (the sign is zero) when $t$ is a zero of a common factor of $f$ and $g$. We present an efficient degeneracy detection algorithm based on a one-time factoring of every possible angle polynomial. Our algorithm is 3500 times faster than the standard algorithm based on greatest common divisor computation. It reduces the share of degeneracy detection in our free space computations from 90% to 0.5% of the running time.


title: Row-finite systems of ordinary differential equations in a scale of Banach spaces

abstract: Motivated by the study of dynamics of interacting spins for infinite particle systems, we consider an infinite family of first order differential equations in a Euclidean space, parameterized by elements $x$ of a fixed countable set. We suppose that the system is row-finite, that is, the right-hand side of the $x$-equation depends on a finite but in general unbounded number $n_x$ of variables. Under certain dissipativity-type conditions on the right-hand side and a bound on the growth of $n_x$, we show the existence of the solutions with infinite life-time, and prove that they live in an increasing scale of Banach spaces. For this, we obtain uniform estimates for solutions to approximating finite systems using a version of Ovsyannikov's method for linear systems in a scale of Banach spaces. As a by-product, we develop an infinite-time generalization of the Ovsyannikov method.


title: Exploring neutrino mass and mass hierarchy in the scenario of vacuum energy interacting with cold dark matter

abstract: We investigate the constraints on total neutrino mass in the scenario of vacuum energy interacting with cold dark matter. We focus on two typical interaction forms, i.e., $Q=\beta H\rho_{\rm c}$ and $Q=\beta H\rho_{\Lambda}$. To avoid the occurrence of large-scale instability in interacting dark energy cosmology, we adopt the parameterized post-Friedmann approach to calculate the perturbation evolution of dark energy. We employ observational data, including the Planck cosmic microwave background temperature and polarization data, baryon acoustic oscillation data, a JLA sample of type Ia supernovae observation, direct measurement of the Hubble constant, and redshift space distortion data. We find that, compared with those in the $\Lambda$CDM model, much looser constraints on $\sum m_{\nu}$ are obtained in the $Q=\beta H\rho_{\rm c}$ model, whereas slightly tighter constraints are obtained in the $Q=\beta H\rho_{\Lambda}$ model. Consideration of the possible mass hierarchies of neutrinos reveals that the smallest upper limit of $\sum m_{\nu}$ appears in the degenerate hierarchy case. By comparing the values of $\chi^2_{\rm min}$, we find that the normal hierarchy case is favored over the inverted one. In particular, we find that the difference $\Delta \chi^2_{\rm min} \equiv \chi^2_{\rm IH; min}-\chi^2_{\rm NH; min}> 2$ in the $Q=\beta H\rho_{\rm c}$ model. In addition, we find that $\beta=0$ is consistent with the current observations in the $Q=\beta H\rho_{\rm c}$ model, and $\beta < 0$ is favored at more than the $1\sigma$ level in the $Q=\beta H\rho_{\Lambda}$ model.


title: Unsupervised Semantic Deep Hashing

abstract: In recent years, deep hashing methods have been proved to be efficient since it employs convolutional neural network to learn features and hashing codes simultaneously. However, these methods are mostly supervised. In real-world application, it is a time-consuming and overloaded task for annotating a large number of images. In this paper, we propose a novel unsupervised deep hashing method for large-scale image retrieval. Our method, namely unsupervised semantic deep hashing (\textbf{USDH}), uses semantic information preserved in the CNN feature layer to guide the training of network. We enforce four criteria on hashing codes learning based on VGG-19 model: 1) preserving relevant information of feature space in hashing space; 2) minimizing quantization loss between binary-like codes and hashing codes; 3) improving the usage of each bit in hashing codes by using maximum information entropy, and 4) invariant to image rotation. Extensive experiments on CIFAR-10, NUSWIDE have demonstrated that \textbf{USDH} outperforms several state-of-the-art unsupervised hashing methods for image retrieval. We also conduct experiments on Oxford 17 datasets for fine-grained classification to verify its efficiency for other computer vision tasks.


title: Non-linear double-peeling: experimental vs. theoretical predictions

abstract: The double peeling of detachment of non-linear adhesive tapes from a flat Poly(methylmethacrylate) (PMMA) surface has been investigated from both experimental and theoretical point of view. Double peeling tests show that, as the detachment process advances, the peeling angle stabilizes on a limiting value {\theta}lim corresponding to a critical pull-off force Fc above which the tape is completely detached from the substrate. This observed behavior is in good agreement with results obtained following the new theory of multiple peeling and taking into account the hardening-softening non-linear behavior of the experimentally tested adhesive tapes and clarifies some aspects of the experimental data. In particular, the theoretical model shows that the value of the limiting peeling angle depends on the geometry of the adhesive tape as well as on the stiffness properties and on the interfacial energy {\Delta}{\gamma}. Finally, theoretical predictions confirm that solutions with a peeling angle lower than {\theta}lim are unstable.


title: Newton: Gravitating Towards the Physical Limits of Crossbar Acceleration

abstract: Many recent works have designed accelerators for Convolutional Neural Networks (CNNs). While digital accelerators have relied on near data processing, analog accelerators have further reduced data movement by performing in-situ computation. Recent works take advantage of highly parallel analog in-situ computation in memristor crossbars to accelerate the many vector-matrix multiplication operations in CNNs. However, these in-situ accelerators have two significant short-comings that we address in this work. First, the ADCs account for a large fraction of chip power and area. Second, these accelerators adopt a homogeneous design where every resource is provisioned for the worst case. By addressing both problems, the new architecture, Newton, moves closer to achieving optimal energy-per-neuron for crossbar accelerators. We introduce multiple new techniques that apply at different levels of the tile hierarchy. Two of the techniques leverage heterogeneity: one adapts ADC precision based on the requirements of every sub-computation (with zero impact on accuracy), and the other designs tiles customized for convolutions or classifiers. Two other techniques rely on divide-and-conquer numeric algorithms to reduce computations and ADC pressure. Finally, we place constraints on how a workload is mapped to tiles, thus helping reduce resource provisioning in tiles. For a wide range of CNN dataflows and structures, Newton achieves a 77% decrease in power, 51% improvement in energy efficiency, and 2.2x higher throughput/area, relative to the state-of-the-art ISAAC accelerator.


title: Mixing Time of Markov chain of the Knapsack Problem

abstract: To find the number of assignments of zeros and ones satisfying a specific Knapsack Problem is $\#P$ hard, so only approximations are envisageable. A Markov chain allowing uniform sampling of all possible solutions is given by Luby, Randall and Sinclair. In 2005, Morris and Sinclair, by using a flow argument, have shown that the mixing time of this Markov chain is $\mathcal{O}(n^{9/2+\epsilon})$, for any $\epsilon > 0$. By using a canonical path argument on the distributive lattice structure of the set of solutions, we obtain an improved bound, the mixing time is given as $\tau_{_{x}}(\epsilon) \leq n^{3} \ln (16 \epsilon^{-1})$.


title: Exploiting symmetry in network analysis

abstract: Virtually all network analyses involve structural measures between pairs of vertices, or of the vertices themselves, and the large amount of symmetry present in real-world complex networks is inherited by such measures. This has practical consequences which have not yet been explored in full generality, nor systematically exploited by network practitioners. Here we study the effect of network symmetry on arbitrary network measures, and show how this can be exploited in practice in a number of ways, from redundancy compression, to computational reduction. We also uncover the spectral signatures of symmetry for an arbitrary network measure such as the graph Laplacian. Computing network symmetries is very efficient in practice, and we test real-world examples up to several million nodes. Since network models are ubiquitous in the Applied Sciences, and typically contain a large degree of structural redundancy, our results are not only significant, but widely applicable.

categories: math.CO cs.SI physics.data-an physics.soc-ph


title: Simulating the future urban growth in Xiongan New Area: a upcoming big city in China

abstract: China made the announement to create the Xiongan New Area in Hebei in April 1,2017. Thus a new magacity about 110km south west of Beijing will emerge. Xiongan New Area is of great practial significant and historical significant for transferring Beijing's non-capital function. Simulating the urban dynamics in Xiongan New Area can help planners to decide where to build the new urban and further manage the future urban growth. However, only a little research focus on the future urban development in Xiongan New Area. In addition, previous models are unable to simulate the urban dynamics in Xiongan New Area. Because there are no original high density urbna for these models to learn the transition rules.In this study, we proposed a C-FLUS model to solve such problems. This framework was implemented by coupling a modified Cellular automata(CA). An elaborately designed random planted seeds machanism based on local maximums is addressed in the CA model to better simulate the occurrence of the new urban. Through an analysis of the current driving forces, the C-FLUS can detect the potential start zone and simulate the urban development under different scenarios in Xiongan New Area. Our study shows that the new urban is most likely to occur in northwest of Xiongxian, and it will rapidly extend to Rongcheng and Anxin until almost cover the northern part of Xiongan New Area. Moreover, the method can help planners to evaluate the impact of urban expansion in Xiongan New Area.


title: Universal features of price formation in financial markets: perspectives from Deep Learning

abstract: Using a large-scale Deep Learning approach applied to a high-frequency database containing billions of electronic market quotes and transactions for US equities, we uncover nonparametric evidence for the existence of a universal and stationary price formation mechanism relating the dynamics of supply and demand for a stock, as revealed through the order book, to subsequent variations in its market price. We assess the model by testing its out-of-sample predictions for the direction of price moves given the history of price and order flow, across a wide range of stocks and time periods. The universal price formation model is shown to exhibit a remarkably stable out-of-sample prediction accuracy across time, for a wide range of stocks from different sectors. Interestingly, these results also hold for stocks which are not part of the training sample, showing that the relations captured by the model are universal and not asset-specific. The universal model --- trained on data from all stocks --- outperforms, in terms of out-of-sample prediction accuracy, asset-specific linear and nonlinear models trained on time series of any given stock, showing that the universal nature of price formation weighs in favour of pooling together financial data from various stocks, rather than designing asset- or sector-specific models as commonly done. Standard data normalizations based on volatility, price level or average spread, or partitioning the training data into sectors or categories such as large/small tick stocks, do not improve training results. On the other hand, inclusion of price and order flow history over many past observations is shown to improve forecasting performance, showing evidence of path-dependence in price dynamics.


title: Correcting Observation Model Error in Data Assimilation

abstract: Standard methods of data assimilation assume prior knowledge of a model that describes the system dynamics and an observation function that maps the model state to a predicted output. An accurate mapping from model state to observation space is crucial in filtering schemes when adjusting the estimate of the system state during the filter's analysis step. However, in many applications the true observation function may be unknown and the available observation model may have significant errors, resulting in a suboptimal state estimate. We propose a method for observation model error correction within the filtering framework. The procedure involves an alternating minimization algorithm used to iteratively update a given observation function to increase consistency with the model and prior observations, using ideas from attractor reconstruction. The method is demonstrated on the Lorenz 1963 and Lorenz 1996 models, and on a single-column radiative transfer model with multicloud parameterization.


title: Definition and Identification of Information Storage and Processing Capabilities as Possible Markers for Turing-universality in Cellular Automata

abstract: To identify potential universal cellular automata, a method is developed to measure information processing capacity of elementary cellular automata. We consider two features of cellular automata: Ability to store information, and ability to process information. We define local collections of cells as particles of cellular automata and consider information contained by particles. By using this method, information channels and channels' intersections can be shown. By observing these two features, potential universal cellular automata are classified into a certain class, and all elementary cellular automata can be classified into four groups, which correspond to S. Wolfram's four classes: 1) Homogeneous; 2) Regular; 3) Chaotic and 4) Complex. This result shows that using abilities of store and processing information to characterize complex systems is effective and succinct. And it is found that these abilities are capable of quantifying the complexity of systems.


title: On the Keldysh Problem of Flutter Suppression

abstract: This work is devoted to the Keldysh model of flutter suppression and rigorous approaches to its analysis. To solve the stabilization problem in the Keldysh model we use an analog of direct Lyapunov method for differential inclusions. The results obtained here are compared with the results of Keldysh obtained by the method of harmonic balance (describing function method), which is an approximate method for analyzing the existence of periodic solutions. The limitations of the use of describing function method for the study of systems with dry friction and stationary segment are demonstrated.


title: Approximating Flexibility in Distributed Energy Resources: A Geometric Approach

abstract: With increasing availability of communication and control infrastructure at the distribution systems, it is expected that the distributed energy resources (DERs) will take an active part in future power systems operations. One of the main challenges associated with integration of DERs in grid planning and control is in estimating the available flexibility in a collection of (heterogeneous) DERs, each of which may have local constraints that vary over time. In this work, we present a geometric approach for approximating the flexibility of a DER in modulating its active and reactive power consumption. The proposed method is agnostic about the type and model of the DERs, thereby facilitating a plug-and-play approach, and allows scalable aggregation of the flexibility of a collection of (heterogeneous) DERs at the distributed system level. Simulation results are presented to demonstrate the performance of the proposed method.


title: Approximation of Some Multivariate Risk Measures for Gaussian Risks

abstract: Gaussian random vectors exhibit the loss of dimension phenomena, which relate to their joint survival tail behaviour. Besides, the fact that the components of such vectors are light-tailed complicates the approximations of various multivariate risk measures significantly. In this contribution we derive precise approximations of marginal mean excess, marginal expected shortfall and multivariate conditional tail expectation of Gaussian random vectors and highlight links with conditional limit theorems. Our study indicates that similar results hold for elliptical and Gaussian like multivariate risks.


title: Spontaneous symmetry breaking and Higgs mode: comparing Gross-Pitaevskii and nonlinear Klein-Gordon equations

abstract: We discuss the mechanism of spontaneous symmetry breaking and the elementary excitations for a weakly-interacting Bose gas at finite temperature. We consider both the non-relativistic case, described by the Gross-Pitaevskii equation, and the relativistic one, described by the cubic nonlinear Klein-Gordon equation. We analyze similarities and differences in the two equations and, in particular, in the phase and amplitude modes (i.e. Goldstone and Higgs modes) of the bosonic matter field. We show that the coupling between phase and amplitude modes gives rise to a single gapless Bogoliubov spectrum in the non-relativistic case. Instead, in the relativistic case the spectrum has two branches: one is gapless and the other is gapped. In the non-relativistic limit we find that the relativistic spectrum reduces to the Bogoliubov one. Finally, as an application of the above analysis, we consider the Bose-Hubbard model close to the superfluid-Mott quantum phase transition and we investigate the elementary excitations of its effective action, which contains both non-relativistic and relativistic terms.


title: Cloud Workload Prediction based on Workflow Execution Time Discrepancies

abstract: Infrastructure as a service clouds hide the complexity of maintaining the physical infrastructure with a slight disadvantage: they also hide their internal working details. Should users need knowledge about these details e.g., to increase the reliability or performance of their applications, they would need solutions to detect behavioural changes in the underlying system. Existing runtime solutions for such purposes offer limited capabilities as they are mostly restricted to revealing weekly or yearly behavioural periodicity in the infrastructure. This article proposes a technique for predicting generic background workload by means of simulations that are capable of providing additional knowledge of the underlying private cloud systems in order to support activities like cloud orchestration or workflow enactment. Our technique uses long-running scientific workflows and their behaviour discrepancies and tries to replicate these in a simulated cloud with known (trace-based) workloads. We argue that the better we can mimic the current discrepancies the better we can tell expected workloads in the near future on the real life cloud. We evaluated the proposed prediction approach with a biochemical application on both real and simulated cloud infrastructures. The proposed algorithm has shown to produce significantly (~20%) better workload predictions for the future of simulated clouds than random workload selection.


title: (Parametrized) First Order Transport Equations: Realization of Optimally Stable Petrov-Galerkin Methods

abstract: We consider ultraweak variational formulations for (parametrized) linear first order transport equations in time and/or space. Computationally feasible pairs of optimally stable trial and test spaces are presented, starting with a suitable test space and defining an optimal trial space by the application of the adjoint operator. As a result, the inf-sup constant is one in the continuous as well as in the discrete case and the computational realization is therefore easy. In particular, regarding the latter, we avoid a stabilization loop within the greedy algorithm when constructing reduced models within the framework of reduced basis methods. Several numerical experiments demonstrate the good performance of the new method.


title: Supersymmetric Solutions of N = (1,1) General Massive Supergravity

abstract: We construct supersymmetric solutions of three dimensional N = (1,1) General Massive Supergravity (GMG). Solutions with a null Killing vector are in general pp-waves. We identify those that appear at critical points of the model some of which do not exist in N = (1,1) New Massive Supergravity (NMG). In the timelike case, we find that many solutions are common with NMG but there is a new class that is genuine to GMG, two members of which are stationary Lifshitz and timelike squashed AdS spacetimes. We also show that in addition to the fully supersymmetric AdS vacuum, there is a second AdS background with a non-zero vector field that preserves 1/4 supersymmetry.


title: Ultra Diffuse Galaxies are a Subset of Cluster Dwarf Elliptical/Spheroidal Galaxies

abstract: Since 2015 there has been a great deal of interest in a supposed new class of galaxy called Ultra Diffuse Galaxies (UDGs). These are large systems with sizes $> 1.5$ kpc and have surface brightness values which are $\mu > 25$ mag arcsec$^{-2}$. Because of their low-surface brightness they are proposed to be `failed' Milky Way type galaxies given their similar size, but much lower stellar masses. As such, these systems are considered by some as a new type of galaxy, yet we show that they are a subset of a well-established and well studied population of low-surface brightness galaxies found mostly in dense areas of the universe - clusters of galaxies. We argue based on previous literature that the most likely method for forming these galaxies is through cluster processes such as `Galaxy Harassment', where through multiple high speed encounters an infalling galaxy is gradually removed of its mass, until it resembles a dwarf elliptical. Future studies of UDGs should consider the above and their more general connection to previously studied populations.


title: Optimization Based Solutions for Control and State Estimation in Non-holonomic Mobile Robots: Stability, Distributed Control, and Relative Localization

abstract: Interest in designing, manufacturing, and using autonomous robots has been rapidly growing during the most recent decade. The main motivation for this interest is the wide range of potential applications these autonomous systems can serve in. The applications include, but are not limited to, area coverage, patrolling missions, perimeter surveillance, search and rescue missions, and situational awareness. In this thesis, the area of control and state estimation in non-holonomic mobile robots is tackled. Herein, optimization based solutions for control and state estimation are designed, analyzed, and implemented to such systems. One of the main motivations for considering such solutions is their ability of handling constrained and nonlinear systems such as non-holonomic mobile robots. Moreover, the recent developments in dynamic optimization algorithms as well as in computer processing facilitated the real-time implementation of such optimization based methods in embedded computer systems.


title: Stochastic filtering and optimal control of pure jump Markov processes with noise-free partial observation

abstract: We consider an infinite horizon optimal control problem for a pure jump Markov process $X$, taking values in a complete and separable metric space $I$, with noise-free partial observation. The observation process is defined as $Y_t = h(X_t)$, $t \geq 0$, where $h$ is a given map defined on $I$. The observation is noise-free in the sense that the only source of randomness is the process $X$ itself. The aim is to minimize a discounted cost functional. In the first part of the paper we write down an explicit filtering equation and characterize the filtering process as a Piecewise Deterministic Process. In the second part, after transforming the original control problem with partial observation into one with complete observation (the separated problem) using filtering equations, we prove the equivalence of the original and separated problems through an explicit formula linking their respective value functions. The value function of the separated problem is also characterized as the unique fixed point of a suitably defined contraction mapping.


title: Explicit formula for the density of local times of Markov Jump Processes

abstract: In this note we show a simple formula for the joint density of local times, last exit tree and cycling numbers of continuous-time Markov Chains on finite graphs, which involves the modified Bessel function of the first type.


title: Note on Calder\'on's inverse problem for measurable conductivities

abstract: The unique determination of a measurable conductivity from the Dirichlet-to-Neumann map of the equation $\mathrm{div} (\sigma \nabla u) = 0$ is the subject of this note. A new strategy, based on Clifford algebras and a higher dimensional analogue of the Beltrami equation, is here proposed. This represents a possible first step for a proof of uniqueness for the Calder\'on problem in three and higher dimensions in the $L^\infty$ case.


title: Orders of Tate-Shafarevich groups for the cubic twists of $X_0(27)$

abstract: This paper continues the authors previous investigations concerning orders of Tate-Shafarevich groups in quadratic twists of a given elliptic curve, and for the family of the Neumann-Setzer type elliptic curves. Here we present the results of our search for the (analytic) orders of Tate-Shafarevich groups for the cubic twists of $X_0(27)$. Our calculations extend those given by Zagier and Kramarz \cite{ZK} and by Watkins \cite{Wat}. Our main observations concern the asymptotic formula for the frequency of orders of Tate-Shafarevich groups. In the last section we propose a similar asymptotic formula for the class numbers of real quadratic fields.


title: The canonical projection associated to certain possibly infinite generalized iterated function system as a fixed point

abstract: In this paper, influenced by the ideas from A. Mihail, The canonical projection between the shift space of an IIFS and its attractor as a fixed point, Fixed Point Theory Appl., 2015, Paper No. 75, 15 p., we associate to every generalized iterated function system F (of order m) an operator H defined on C^m and taking values on C, where C stands for the space of continuous functions from the shift space on the metric space corresponding to the system. We provide sufficient conditions (on the constitutive functions of F) for the operator H to be continuous, contraction, phi-contraction, Meir-Keeler or contractive. We also give sufficient condition under which H has a unique fixed point. Moreover, we prove that, under these circumstances, the closer of the imagine of the fixed point is the attractor of F and that the fixed point is the canonical projection associated to F. In this way we give a partial answer to the open problem raised on the last paragraph of the above mentioned Mihail's paper.


title: PyGOM - A Python Package for Simplifying Modelling with Systems of Ordinary Differential Equations

abstract: Ordinary Differential Equations (ODE) are used throughout science where the capture of rates of change in states is sought. While both pieces of commercial and open software exist to study such systems, their efficient and accurate usage frequently requires deep understanding of mathematics and programming. The package we present here, PyGOM, seeks to remove these obstacles for models based on ODE systems. We provide a simple interface for the construction of such systems backed by a comprehensive and easy to use tool--box. This tool--box implements functions to easily perform common operations for ODE systems such as solving, parameter estimation, and stochastic simulation. The package source is freely available and organized in a way that permits easy extension. With both the algebraic and numeric calculations performed automatically (but still accessible), the end user is freed to focus on model development.


title: Mathematical Analysis of Anthropogenic Signatures: The Great Deceleration

abstract: Distributions of anthropogenic signatures (impacts and activities) are mathematically analysed. The aim is to understand the Anthropocene and to see whether anthropogenic signatures could be used to determine its beginning. A total of 23 signatures were analysed and results are presented in 31 diagrams. Some of these signatures contain undistinguishable natural components but most of them are of purely anthropogenic origin. Great care was taken to identify abrupt accelerations, which could be used to determine the beginning of the Anthropocene. Results of the analysis can be summarised in three conclusions. 1. Anthropogenic signatures cannot be used to determine the beginning of the Anthropocene. 2. There was no abrupt Great Acceleration around 1950 or around any other time. 3. Anthropogenic signatures are characterised by the Great Deceleration in the second half of the 20th century. The second half of the 20th century does not mark the beginning of the Anthropocene but most likely the beginning of the end of the strong anthropogenic impacts, maybe even the beginning of a transition to a sustainable future. The Anthropocene is a unique stage in human experience but it has no clearly marked beginning and it is probably not a new geological epoch.


title: Inverse Visual Question Answering: A New Benchmark and VQA Diagnosis Tool

abstract: In recent years, visual question answering (VQA) has become topical. The premise of VQA's significance as a benchmark in AI, is that both the image and textual question need to be well understood and mutually grounded in order to infer the correct answer. However, current VQA models perhaps `understand' less than initially hoped, and instead master the easier task of exploiting cues given away in the question and biases in the answer distribution. In this paper we propose the inverse problem of VQA (iVQA). The iVQA task is to generate a question that corresponds to a given image and answer pair. We propose a variational iVQA model that can generate diverse, grammatically correct and content correlated questions that match the given answer. Based on this model, we show that iVQA is an interesting benchmark for visuo-linguistic understanding, and a more challenging alternative to VQA because an iVQA model needs to understand the image better to be successful. As a second contribution, we show how to use iVQA in a novel reinforcement learning framework to diagnose any existing VQA model by way of exposing its belief set: the set of question-answer pairs that the VQA model would predict true for a given image. This provides a completely new window into what VQA models `believe' about images. We show that existing VQA models have more erroneous beliefs than previously thought, revealing their intrinsic weaknesses. Suggestions are then made on how to address these weaknesses going forward.


title: Pearson's correlation coefficient in the theory of networks: A comment

abstract: In statistics, the Pearson correlation coefficient $r_{x,y}$ determines the degree of linear correlation between two variables and it is known that $-1 \le r_{x,y} \le 1$. In the theory of networks, a curious expression proposed in [PRL {\bf 89} 208701 (2002)] for degree-degree correlation coefficient $r_{j_i,k_i}, i\in [1,M]$ has been in use. We realize that the suggested form is the conventional Pearson's coefficient for $\{(j_i,k_i), (k_i,j_i)\}$ for $2M$ data points and hence it is rightly dedicated to undirected networks.


title: Conformal amplitude hierarchy and the Poincare disk

abstract: The amplitude for the singlet channels in the 4-point function of the fundamental field in the conformal field theory of the 2d $O(n)$ model is studied as a function of $n$. For a generic value of $n$, the 4-point function has infinitely many amplitudes, whose landscape can be very spiky as the higher amplitude changes its sign many times at the simple poles, which generalize the unique pole of the energy operator amplitude at $n=0$. In the stadard parameterization of $n$ by angle in unit of $\pi$, we find that the zeros and poles happen at the rational angles, forming a hierarchical tree structure inherent in the Poincar\'{e} disk. Some relation between the amplitude and the Farey path, a piecewise geodesic that visits these zeros and poles, is suggested. In this hierarchy, the symmetry of the congruence subgroup $\Gamma(2)$ of $SL(2,\mathbb{Z})$ naturally arises from the two clearly distinct even/odd classes of the rational angles, in which one respectively gets the truncated operator algebras and the logarithmic 4-point functions.

categories: hep-th cond-mat.stat-mech math-ph math.MP


title: Motion of vortices in ferromagnetic spin-1 BEC

abstract: The paper investigates dynamics of nonsingular vortices in a ferromagnetic spin-1 BEC, where spin and mass superfluidity coexist in the presence of uniaxial anisotropy (linear and quadratic Zeeman effect). The analysis is based on hydrodynamics following from the Gross-Pitaevskii theory. Cores of nonsingular vortices are skyrmions with charge, which is tuned by uniaxial anisotropy and can have any fractal value between 0 and 1. There are circulations of mass and spin currents around these vortices. The results are compared with the equation of vortex motion derived earlier in the Landau-Lifshitz-Gilbert theory for magnetic vortices in easy-plane ferromagnetic insulators. In the both cases the transverse gyrotropic force (analog of the Magnus force in superfluid and classical hydrodynamics) is proportional to the charge of skyrmions in vortex cores.


title: The ECMWF Ensemble Prediction System: Looking Back (more than) 25 Years and Projecting Forward 25 Years

abstract: This paper has been written to mark 25 years of operational medium-range ensemble forecasting. The origins of the ECMWF Ensemble Prediction System are outlined, including the development of the precursor real-time Met Office monthly ensemble forecast system. In particular, the reasons for the development of singular vectors and stochastic physics - particular features of the ECMWF Ensemble Prediction System - are discussed. The author speculates about the development and use of ensemble prediction in the next 25 years.