title: Transmission Game in MIMO Interference Channels With Radio-Frequency Energy Harvesting

abstract: For multi-user transmissions over MIMO interference channels, each user designs the transmit covariance matrix to maximize its information rate. When passive radio-frequency (RF) energy harvesters are present in the network, the transmissions are constrained by both the transmit power limits and the energy harvesting requirements. A passive RF energy harvester collects the radiated energy from nearby wireless information transmitters instead of using a dedicated wireless power source. It needs multiple transmitters to concentrate their RF radiation on it because typical electric field strengths are weak. In this paper, strategic games are proposed for the multi-user transmissions. First, in a non-cooperative game, each transmitter has a best-response strategy for the transmit covariance matrix that follows a multi-level water-filling solution. A pure-strategy Nash equilibrium exists. Secondly, in a cooperative game, there is no need to estimate the proportion of the harvested energy from each transmitter. Rather, the transmitters bargain over the unit-reward of the energy contribution. An approximation of the information rate is used in constructing the individual utility such that the problem of network utility maximization can be decomposed and the bargaining process can be implemented distributively. The bargaining solution gives a point of rates that is superior to the Nash equilibria and close to the Pareto front. Simulation results verify the algorithms that provide good communication performance while satisfying the RF energy-harvesting requirements.


title: Generalized Theory of Smallest Diameter of Metallic Nanorods

abstract: This Letter reports a generalized theory of the smallest diameter of metallic nanorods from physical vapor deposition. The generalization incorporates the effects of nanorod separation and those of van der Waals interactions on geometrical shadowing. In contrast, the previous theory for idealized geometrical shadowing [Phys. Rev. Lett. 110, 136102 (2013)] does not include any dependence on nanorod separation and it predicts the diameter to be about to of what the generalized theory does. As verification, numerical solutions and the generalized theory in closed-form agree in terms of effective deposition flux. As validation, experiments of physical vapor deposition and the generalized theory agree in terms of the diameter as a function of the separation of nanorods.


title: Gauged BPS baby Skyrmions with quantised magnetic flux

abstract: A new type of gauged BPS baby Skyrme model is presented, where the derivative term is just the Schroers current (i.e., gauge invariant and conserved version of the topological current) squared. This class of models has a topological bound saturated for solutions of the pertinent Bogomolnyi equations supplemented by a so-called superpotential equation. In contrast to the gauged BPS baby Skyrme models considered previously, the superpotential equation is linear and, hence, completely solvable. Furthermore, the magnetic flux is quantized in units of $2\pi$, which allows, in principle, to define this theory on a compact manifold without boundary, unlike all gauged baby Skyrme models considered so far.


title: The Sardinia Radio Telescope: From a Technological Project to a Radio Observatory

abstract: [Abridged] The Sardinia Radio Telescope (SRT) is the new 64-m dish operated by INAF (Italy). Its active surface will allow us to observe at frequencies of up to 116 GHz. At the moment, three receivers, one per focal position, have been installed and tested. The SRT was officially opened in October 2013, upon completion of its technical commissioning phase. In this paper, we provide an overview of the main science drivers for the SRT, describe the main outcomes from the scientific commissioning of the telescope, and discuss a set of observations demonstrating the SRT's scientific capabilities. One of the main objectives of scientific commissioning was the identification of deficiencies in the instrumentation and/or in the telescope sub-systems for further optimization. As a result, the overall telescope performance has been significantly improved. As part of the scientific commissioning activities, different observing modes were tested and validated, and first astronomical observations were carried out to demonstrate the science capabilities of the SRT. In addition, we developed astronomer-oriented software tools, to support future observers on-site. The astronomical validation activities were prioritized based on technical readiness and scientific impact. The highest priority was to make the SRT available for joint observations as part of European networks. As a result, the SRT started to participate (in shared-risk mode) in EVN (European VLBI Network) and LEAP (Large European Array for Pulsars) observing sessions in early 2014. The validation of single-dish operations for the suite of SRT first light receivers and backends continued in the following years, and was concluded with the first call for shared-risk/early-science observations issued at the end of 2015.


title: Global regularity and fast small scale formation for Euler patch equation in a disk

abstract: It is well known that the Euler vortex patch in $\mathbb{R}^{2}$ will remain regular if it is regular enough initially. In bounded domains, the regularity theory for patch solutions is less complete. We study here the Euler vortex patch in a disk. We prove global in time regularity by providing the upper bound of the growth of curvature of the patch boundary. For a special symmetric scenario, we construct an example of double exponential curvature growth, showing that such upper bound is qualitatively sharp.


title: Higgsino Dark Matter or Not: Role of Disappearing Track Searches at the LHC and Future Colliders

abstract: Higgsino in supersymmetric standard models is known to be a promising candidate for dark matter in the Universe. Its phenomenological property is strongly affected by the gaugino fraction in the Higgsino-like state. If this is sizable, in other words, if gaugino masses are less than ${\cal O}(10)$ TeV, we may probe the Higgsino dark matter in future non-accelerator experiments such as dark matter direct searches and measurements of electric dipole moments. On the other hand, if gauginos are much heavier, then it is hard to search for Higgsino in these experiments. In this case, due to a lack of gaugino components, the mass difference between the neutral and charged Higgsinos is uniquely determined by electroweak interactions to be around $350$ MeV, which makes the heavier charged state rather long-lived, with a decay length of about $1$ cm. In this letter, we argue that a charged particle with a flight length of ${\cal O}(1)$ cm can be probed in disappearing-track searches if we require only two hits in the pixel detector. Even in this case, we can reduce background events with the help of the displaced-vertex reconstruction technique. We study the prospects of this search strategy at the LHC and future colliders for the Higgsino dark matter scenario. It is found that an almost pure Higgsino is indeed within the reach of the future $33$ TeV collider experiments. We then discuss that the interplay among collider and non-accelerator experiments plays a crucial role in testing the Higgsino dark matter scenarios. Our strategy for disappearing-track searches can also enlarge the discovery potential of pure wino dark matter as well as other electroweak-charged dark matter candidates.


title: Bright Breathers in Nonlinear Left-Handed Metamaterial Lattices

abstract: In the present work, we examine a prototypical model for the formation of bright breathers in nonlinear left-handed metamaterial lattices. Utilizing the paradigm of nonlinear transmission lines, we build a relevant lattice and develop a quasi-continuum multiscale approximation that enables us to appreciate both the underlying linear dispersion relation and the potential for bifurcation of nonlinear states. We focus here, more specifically, on bright discrete breathers which bifurcate from the lower edge of the linear dispersion relation at wavenumber $k=\pi$. Guided by the multiscale analysis, we calculate numerically both the stable inter-site centered and the unstable site-centered members of the relevant family. We quantify the associated stability via Floquet analysis and the Peierls-Nabarro barrier of the energy difference between these branches. Finally, we explore the dynamical implications of these findings towards the potential mobility or lack thereof (pinning) of such breather solutions.


title: A multiplier algebra functional calculus

abstract: This paper generalizes the classical Sz.-Nagy--Foias $H^{\infty}(\mathbb{D})$ functional calculus for Hilbert space contractions. In particular, we replace the single contraction $T$ with a tuple $T=(T_1, \dots, T_d)$ of commuting bounded operators on a Hilbert space and replace $H^{\infty}(\mathbb{D})$ with a large class of multiplier algebras of Hilbert function spaces on the unit ball in $\mathbb C^d$.


title: Isolating Majorana fermions with finite Kitaev nanowires and temperature: the universality of the zero-bias conductance

abstract: The zero-bias peak (ZBP) is understood as the definite signature of a Majorana bound state (MBS) when attached to a semi-infinite Kitaev nanowire (KNW) nearby zero temperature. However, such characteristics concerning the realization of the KNW constitute a profound experimental challenge. We explore theoretically a QD connected to a topological KNW of finite size at non-zero temperatures and show that overlapped MBSs of the wire edges can become effectively decoupled from each other and the characteristic ZBP can be fully recovered if one tunes the system into the leaked Majorana fermion fixed point. At very low temperatures, the MBSs become strongly coupled similarly to what happens in the Kondo effect. We derive universal features of the conductance as a function of the temperature and the relevant crossover temperatures. Our findings offer additional guides to identify signatures of MBSs in solid state setups.


title: The bumpy light curve of supernova iPTF13z

abstract: A Type IIn supernova (SN) is dominated by the interaction of SN ejecta with the circumstellar medium (CSM). Some SNe IIn (e.g., SN 2006jd) have episodes of re-brightening ("bumps") in their light curves. We present iPTF13z, a SN IIn discovered by the intermediate Palomar Transient Factory (iPTF) and characterised by several bumps in its light curve. We analyse this peculiar behaviour trying to infer the properties of the CSM and of the SN explosion, as well as the nature of its progenitor star. We obtained multi-band optical photometry for over 1000 days after discovery with the P48 and P60 telescopes at Palomar Observatory. We obtained low-resolution optical spectra in the same period. We did an archival search for progenitor outbursts. We analyse our photometry and spectra, and compare iPTF13z to other SNe IIn. A simple analytical model is used to estimate properties of the CSM. iPTF13z was a SN IIn showing a light curve with five bumps during its decline phase. The bumps had amplitudes between 0.4 and 0.9 mag and durations between 20 and 120 days. The most prominent bumps appeared in all our different optical bands. The spectra showed typical SN IIn characteristics, with emission lines of H$\alpha$ (with broad component FWHM ~$10^{3}-10^{4} ~{\rm ~km ~s^{-1}}$ and narrow component FWHM ~$10^2 \rm ~km ~s^{-1}$) and He I, but also with Fe II, Ca II, Na I D and H$\beta$ P-Cygni profiles (with velocities of ~$10^{3}$ ${\rm ~km ~s^{-1}}$). A pre-explosion outburst was identified lasting $\gtrsim 50$ days, with $M_r \approx -15$ mag around 210 days before discovery. Large, variable progenitor mass-loss rates (~> 0.01 $M_{\odot} \rm ~yr^{-1}$) and CSM densities (~> 10$^{-16}$ g cm$^{-3}$) are derived. We suggest that the light curve bumps of iPTF13z arose from SN ejecta interacting with denser regions in the CSM, possibly produced by the eruptions of a luminous blue variable star.

categories: astro-ph.SR astro-ph.CO astro-ph.GA astro-ph.HE


title: Certifying numerical estimates of spectral gaps

abstract: We establish a lower bound on the spectral gap of the Laplace operator on special linear groups using conic optimisation. In particular, this provides a constructive (but computer assisted) proof that these groups have Kazhdan property (T). A software for such optimisation for other finitely presented groups is provided.


title: Paths to equilibrium in non-conformal collisions

abstract: We extend our previous analysis of holographic heavy ion collisions in non-conformal theories. We provide a detailed description of our numerical code. We study collisions at different energies in gauge theories with different degrees of non-conformality. We compare four relaxation times: the hydrodynamization time (when hydrodynamics becomes applicable), the EoSization time (when the average pressure approaches its equilibrium value), the isotropization time (when the longitudinal and transverse pressures approach each other) and the condensate relaxation time (when the expectation value of a scalar operator approaches its equilibrium value). We find that these processes can occur in several different orderings. In particular, the condensate can remain far from equilibrium even long after the plasma has hydrodynamized and EoSized. We also explore the rapidity distribution of the energy density at hydrodynamization. This is far from boost-invariant and its width decreases as the non-conformality increases. Nevertheless, the velocity field at hydrodynamization is almost exactly boost-invariant regardless of the non-conformality. This result may be used to constrain the initialization of hydrodynamic fields in heavy ion collisions.


title: On the Profile of Multiplicities of Complete Subgraphs

abstract: Let $G$ be a $2$-coloring of a complete graph on $n$ vertices, for sufficiently large $n$. We prove that $G$ contains at least $n^{(\frac{1}{4} - o(1))\log n}$ monochromatic complete subgraphs of size $r$, where \[ 0.3\log n < r < 0.7\log n. \] The previously known lower bound on the total number of monochromatic complete subgraphs, due to Sz\'{e}kely was $n^{0.1576\log n}$. We also prove that $G$ contains at least $n^{\frac{1}{7} \log n} $ monochromatic complete subgraphs of size $\frac{1}{2}\log n$. If furthermore one assumes that the largest monochromatic complete subgraph in $G$ is of size $(\frac{1}{2} + o(1))\log n$ (it is a well known open question whether such graphs exist), then for every constant $0 \le c \le \frac{1}{2}$ we determine (up to low order terms) the number of monochromatic complete subgraphs of size $c \log n$. We do so by proving a lower bound that matches (up to low order terms) a previous upper bound of Sz\'{e}kely. For example, the number of monochromatic complete subgraphs of size $\frac{1}{2} \log n$ is $n^{\frac{1}{8}(4 - \log e \pm o(1))\log n} \simeq n^{0.32 \log n}$.


title: Whitney theorem for complex polynomial mappings

abstract: We describe the topology of a general polynomial mapping $F=(f, g):X\to\Bbb C^2$, where $X$ is a complex plane or a complex sphere.


title: An Analysis of Visual Question Answering Algorithms

abstract: In visual question answering (VQA), an algorithm must answer text-based questions about images. While multiple datasets for VQA have been created since late 2014, they all have flaws in both their content and the way algorithms are evaluated on them. As a result, evaluation scores are inflated and predominantly determined by answering easier questions, making it difficult to compare different methods. In this paper, we analyze existing VQA algorithms using a new dataset. It contains over 1.6 million questions organized into 12 different categories. We also introduce questions that are meaningless for a given image to force a VQA system to reason about image content. We propose new evaluation schemes that compensate for over-represented question-types and make it easier to study the strengths and weaknesses of algorithms. We analyze the performance of both baseline and state-of-the-art VQA models, including multi-modal compact bilinear pooling (MCB), neural module networks, and recurrent answering units. Our experiments establish how attention helps certain categories more than others, determine which models work better than others, and explain how simple models (e.g. MLP) can surpass more complex models (MCB) by simply learning to answer large, easy question categories.


title: Renormalization of scalar field theories in rational spacetime dimensions

abstract: We renormalize various scalar field theories with a $\phi^n$ self interaction such as $n$ $=$ $5$, $7$ and $9$ in their respective critical dimensions which are non-integer. The renormalization group functions for the $O(N)$ symmetric extensions are also computed.


title: Nonminimal couplings, gravitational waves, and torsion in Horndeski's theory

abstract: The Horndeski Lagrangian brings together all possible interactions between gravity and a scalar field that yield second-order field equations in four-dimensional spacetime. As originally proposed, it only addresses phenomenology without torsion, which is a non-Riemannian feature of geometry. Since torsion can potentially affect interesting phenomena such as gravitational waves and early Universe inflation, in this paper we allow torsion to exist and propagate within the Horndeski framework. To achieve this goal, we cast the Horndeski Lagrangian in Cartan's first-order formalism, and introduce wave operators designed to act covariantly on p-form fields that carry Lorentz indices. We find that nonminimal couplings and second-order derivatives of the scalar field in the Lagrangian are indeed generic sources of torsion. Metric perturbations couple to the background torsion and new torsional modes appear. These may be detected via gravitational waves but not through Yang-Mills gauge bosons.


title: On multicolor Ramsey numbers for loose $k$-paths of length three

abstract: We show that there exists an absolute constant $A$ such that for each $k\ge2$ and every coloring of the edges of the complete $k$-uniform hypergraph on $ Ar$ vertices with $r$ colors, one of the color classes contains a loose path of length three.


title: Glassiness and Heterogeneous Dynamics in Dense Solutions of Ring Polymers

abstract: Understanding how topological constraints affect the dynamics of polymers in solution is at the basis of any polymer theory and it is particularly needed for melts of rings. These polymers fold as crumpled and space-filling objects and, yet, they display a large number of topological constraints. To understand their role, here we systematically probe the response of solutions of rings at various densities to "random pinning" perturbations. We show that these perturbations trigger non-Gaussian and heterogeneous dynamics, eventually leading to non-ergodic and glassy behaviours. We then derive universal scaling relations for the values of solution density and polymer length marking the onset of vitrification in unperturbed solutions. Finally, we directly connect the heterogeneous dynamics of the rings with their spatial organisation and mutual interpenetration. Our results suggest that deviations from the typical behaviours observed in systems of linear polymers may originate from architecture-specific (threading) topological constraints.


title: The world of long-range interactions: A bird's eye view

abstract: In recent years, studies of long-range interacting (LRI) systems have taken centre stage in the arena of statistical mechanics and dynamical system studies, due to new theoretical developments involving tools from as diverse a field as kinetic theory, non-equilibrium statistical mechanics, and large deviation theory, but also due to new and exciting experimental realizations of LRI systems. In this invited contribution, we discuss the general features of long-range interactions, emphasizing in particular the main physical phenomenon of non-additivity, which leads to a plethora of distinct effects, both thermodynamic and dynamic, that are not observed with short-range interactions: Ensemble inequivalence, slow relaxation, broken ergodicity. We also discuss several physical systems with long-range interactions: mean-field spin systems, self-gravitating systems, Euler equations in two dimensions, Coulomb systems, one-component electron plasma, dipolar systems, free-electron lasers, atoms trapped in optical cavities.


title: Efficient Two-Dimensional Sparse Coding Using Tensor-Linear Combination

abstract: Sparse coding (SC) is an automatic feature extraction and selection technique that is widely used in unsupervised learning. However, conventional SC vectorizes the input images, which breaks apart the local proximity of pixels and destructs the elementary object structures of images. In this paper, we propose a novel two-dimensional sparse coding (2DSC) scheme that represents the input images as the tensor-linear combinations under a novel algebraic framework. 2DSC learns much more concise dictionaries because it uses the circular convolution operator, since the shifted versions of atoms learned by conventional SC are treated as the same ones. We apply 2DSC to natural images and demonstrate that 2DSC returns meaningful dictionaries for large patches. Moreover, for mutli-spectral images denoising, the proposed 2DSC reduces computational costs with competitive performance in comparison with the state-of-the-art algorithms.


title: Self-organized mechano-chemical dynamics in amoeboid locomotion of Physarum fragments

abstract: The aim of this work is to quantify the spatio-temporal dynamics of flow-driven amoeboid locomotion in small ($\sim$100 $\mu$m) fragments of the true slime mold \phys {\it polycephalum}. In this model organism, cellular contraction drives intracellular flows, and these flows transport the chemical signals that regulate contraction in the first place. As a consequence of these non-linear interactions, a diversity of migratory behaviors can be observed in migrating \phys fragments. To study these dynamics, we measure the spatio-temporal distributions of the velocities of the endoplasm and ectoplasm of each migrating fragment, the traction stresses it generates on the substratum, and the concentration of free intracellular calcium. Using these unprecedented experimental data, we classify migrating \phys fragments according to their dynamics, finding that they often exhibit spontaneously coordinated waves of flow, contractility and chemical signaling. We show that \phys fragments exhibiting symmetric spatio-temporal patterns of endoplasmic flow migrate significantly slower than fragments with asymmetric patterns. In addition, our joint measurements of ectoplasm velocity and traction stress at the substratum suggest that forward motion of the ectoplasm is enabled by a succession of stick-slip transitions, which we conjecture are also organized in the form of waves. Combining our experiments with a simplified convection-diffusion model, we show that the convective transport of calcium ions may be key for establishing and maintaining the spatio-temporal patterns of calcium concentration that regulate the generation of contractile forces.


title: pySecDec: a toolbox for the numerical evaluation of multi-scale integrals

abstract: We present pySecDec, a new version of the program SecDec, which performs the factorisation of dimensionally regulated poles in parametric integrals, and the subsequent numerical evaluation of the finite coefficients. The algebraic part of the program is now written in the form of python modules, which allow a very flexible usage. The optimization of the C++ code, generated using FORM, is improved, leading to a faster numerical convergence. The new version also creates a library of the integrand functions, such that it can be linked to user-specific codes for the evaluation of matrix elements in a way similar to analytic integral libraries.


title: Fast link prediction for large networks using spectral embedding

abstract: Many link prediction algorithms require the computation of a similarity metric on each vertex pair, which is quadratic in the number of vertices and infeasible for large networks. We develop a class of link prediction algorithms based on a spectral embedding and the k closest pairs algorithm that are scalable to very large networks. We compare the prediction accuracy and runtime of these methods to existing algorithms on several large link prediction tasks. Our methods achieve comparable accuracy to standard algorithms but are significantly faster.


title: Strained graphene based highly efficient quantum heat engine operating at maximum power

abstract: A strained graphene monolayer is shown to operate as a highly efficient quantum heat engine delivering maximum power. The efficiency and power of the proposed device exceeds that of recent proposals. The reason for these excellent characteristics is that strain enables complete valley separation in transmittance through the device, implying that increasing strain leads to very high Seebeck coefficient as well as lower conductance. In addition, since time-reversal symmetry is unbroken in our system, the proposed strained graphene quantum heat engine can also act as a high performance refrigerator.


title: Semi and Weakly Supervised Semantic Segmentation Using Generative Adversarial Network

abstract: Semantic segmentation has been a long standing challenging task in computer vision. It aims at assigning a label to each image pixel and needs significant number of pixellevel annotated data, which is often unavailable. To address this lack, in this paper, we leverage, on one hand, massive amount of available unlabeled or weakly labeled data, and on the other hand, non-real images created through Generative Adversarial Networks. In particular, we propose a semi-supervised framework ,based on Generative Adversarial Networks (GANs), which consists of a generator network to provide extra training examples to a multi-class classifier, acting as discriminator in the GAN framework, that assigns sample a label y from the K possible classes or marks it as a fake sample (extra class). The underlying idea is that adding large fake visual data forces real samples to be close in the feature space, enabling a bottom-up clustering process, which, in turn, improves multiclass pixel classification. To ensure higher quality of generated images for GANs with consequent improved pixel classification, we extend the above framework by adding weakly annotated data, i.e., we provide class level information to the generator. We tested our approaches on several challenging benchmarking visual datasets, i.e. PASCAL, SiftFLow, Stanford and CamVid, achieving competitive performance also compared to state-of-the-art semantic segmentation method


title: Reversing the irreversible: from limit cycles to emergent time symmetry

abstract: In 1979 Penrose hypothesized that the arrows of time are explained by the hypothesis that the fundamental laws are time irreversible. That is, our reversible laws, such as the standard model and general relativity are effective, and emerge from an underlying fundamental theory which is time irreversible. In Cort\^{e}s and Smolin (2014a, 2014b, 2016) we put forward a research program aiming at realizing just this. The aim is to find a fundamental description of physics above the planck scale, based on irreversible laws, from which will emerge the apparently reversible dynamics we observe on intermediate scales. Here we continue that program and note that a class of discrete dynamical systems are known to exhibit this very property: they have an underlying discrete irreversible evolution, but in the long term exhibit the properties of a time reversible system, in the form of limit cycles. We connect this to our original model proposal in Cort\^{e}s and Smolin (2014a), and show that the behaviours obtained there can be explained in terms of the same phenomenon: the attraction of the system to a basin of limit cycles, where the dynamics appears to be time reversible. Further than that, we show that our original models exhibit the very same feature: the emergence of quasi-particle excitations obtained in the earlier work in the space-time description is an expression of the system's convergence to limit cycles when seen in the causal set description.

categories: gr-qc astro-ph.CO nlin.AO physics.class-ph


title: A Stabilised Nodal Spectral Element Method for Fully Nonlinear Water Waves, Part 2: Wave-body interaction

abstract: We present a new stabilised and efficient high-order nodal spectral element method based on the Mixed Eulerian Lagrangian (MEL) method for general-purpose simulation of fully nonlinear water waves and wave-body interactions. In this MEL formulation a standard Laplace formulation is used to handle arbitrary body shapes using unstructured - possibly hybrid - meshes consisting of high-order curvilinear iso-parametric quadrilateral/triangular elements to represent the body surfaces and for the evolving free surface. Importantly, our numerical analysis highlights that a single top layer of quadrilaterals elements resolves temporal instabilities in the numerical MEL scheme that are known to be associated with mesh topology containing asymmetric element orderings. The 'surface variable only' free surface formulation based on introducing a particle-following (Lagrangian) reference frame contains quartic nonlinear terms that require proper treatment by numerical discretisation due to the possibility of strong aliasing effects. We demonstrate how to stabilise this nonlinear MEL scheme using an efficient combination of (i) global L2 projection without quadrature errors, (ii) mild nonlinear spectral filtering and (iii) re-meshing techniques. Numerical experiments revisiting known benchmarks are presented, and highlights that modelling using a high-order spectral element method provides excellent accuracy in prediction of nonlinear and dispersive wave propagation, and of nonlinear wave-induced loads on fixed submerged and surface-piercing bodies.


title: On singular limit equations for incompressible fluids in moving thin domains

abstract: We consider the incompressible Euler and Navier-Stokes equations in a three-dimensional moving thin domain. Under the assumption that the moving thin domain degenerates into a two-dimensional moving closed surface as the width of the thin domain goes to zero, we give a heuristic derivation of singular limit equations on the degenerate moving surface of the Euler and Navier-Stokes equations in the moving thin domain and investigate relations between their energy structures. We also compare the limit equations with the Euler and Navier-Stokes equations on a stationary manifold, which are described in terms of the Levi-Civita connection.


title: Preface to the 19th workshop "What comes beyond the standard models", Bled July 11--19, 2016, and links to the talks in the proceedings

abstract: The contribution contains the preface to the Proceedings to the 19th Workshop "What Comes Beyond the Standard Models", Bled, July 11 - 19, 2016, published in Bled workshops in physics, Vol.17, No. 2, DMFA-Zaloznistvo, Ljubljana, Dec. 2016, links to (most of) the published contributions and section (by M.Yu. Khlopov) on VIA at Bled 2016.


title: Inverse Reinforcement Learning from Summary Data

abstract: Inverse reinforcement learning (IRL) aims to explain observed strategic behavior by fitting reinforcement learning models to behavioral data. However, traditional IRL methods are only applicable when the observations are in the form of state-action paths. This assumption may not hold in many real-world modeling settings, where only partial or summarized observations are available. In general, we may assume that there is a summarizing function $\sigma$, which acts as a filter between us and the true state-action paths that constitute the demonstration. Some initial approaches to extending IRL to such situations have been presented, but with very specific assumptions about the structure of $\sigma$, such as that only certain state observations are missing. This paper instead focuses on the most general case of the problem, where no assumptions are made about the summarizing function, except that it can be evaluated. We demonstrate that inference is still possible. The paper presents exact and approximate inference algorithms that allow full posterior inference, which is particularly important for assessing parameter uncertainty in this challenging inference situation. Empirical scalability is demonstrated to reasonably sized problems, and practical applicability is demonstrated by estimating the posterior for a cognitive science RL model based on an observed user's task completion time only.


title: Sampling Errors in Nested Sampling Parameter Estimation

abstract: Sampling errors in nested sampling parameter estimation differ from those in Bayesian evidence calculation, but have been little studied in the literature. This paper provides the first explanation of the two main sources of sampling errors in nested sampling parameter estimation, and presents a new diagrammatic representation for the process. We find no current method can accurately measure the parameter estimation errors of a single nested sampling run, and propose a method for doing so using a new algorithm for dividing nested sampling runs. We empirically verify our conclusions and the accuracy of our new method.


title: Factorization and Resummation for Massive Quark Effects in Exclusive Drell-Yan

abstract: Exclusive differential spectra in color-singlet processes at hadron colliders are benchmark observables that have been studied to high precision in theory and experiment. We present an effective-theory framework utilizing soft-collinear effective theory to incorporate massive (bottom) quark effects into resummed differential distributions, accounting for both heavy-quark initiated primary contributions to the hard scattering process as well as secondary effects from gluons splitting into heavy-quark pairs. To be specific, we focus on the Drell-Yan process and consider the vector-boson transverse momentum, $q_T$, and beam thrust, $\mathcal T$, as examples of exclusive observables. The theoretical description depends on the hierarchy between the hard, mass, and the $q_T$ (or $\mathcal T$) scales, ranging from the decoupling limit $q_T \ll m$ to the massless limit $m \ll q_T$. The phenomenologically relevant intermediate regime $m \sim q_T$ requires in particular quark-mass dependent beam and soft functions. We calculate all ingredients for the description of primary and secondary mass effects required at NNLL$'$ resummation order (combining NNLL evolution with NNLO boundary conditions) for $q_T$ and $\mathcal T$ in all relevant hierarchies. For the $q_T$ distribution the rapidity divergences are different from the massless case and we discuss features of the resulting rapidity evolution. Our results will allow for a detailed investigation of quark-mass effects in the ratio of $W$ and $Z$ boson spectra at small $q_T$, which is important for the precision measurement of the $W$-boson mass at the LHC.


title: On the Assembly Rate of Highly Eccentric Binary Black Hole Mergers

abstract: In this {\it Letter} we calculate the fraction of highly eccentric binary black hole (BBH) mergers resulting from binary-single interactions. Using an $N$-body code that includes post-Newtonian correction terms, we show that $\gtrsim 1\%$ of all BBH mergers resulting from this channel will have an eccentricity $e>0.1$ when coming into the LIGO frequency band. As the majority of BBH mergers forming in globular clusters are assembled through three-body encounters, we suggest that such interactions are likely to dominate the population of high eccentricity BBH mergers detectable by LIGO. The relative frequency of highly eccentric events could eventually help to constrain the astrophysical origin of BBH mergers.


title: Using HAWC to Discover Invisible Pulsars

abstract: Observations by HAWC and Milagro have detected bright and spatially extended TeV gamma-ray sources surrounding the Geminga and Monogem pulsars. We argue that these observations, along with a substantial population of other extended TeV sources coincident with pulsar wind nebulae, constitute a new morphological class of spatially extended TeV halos. We show that HAWCs wide field-of-view unlocks an expansive parameter space of TeV halos not observable by atmospheric Cherenkov telescopes. Under the assumption that Geminga and Monogem are typical middle-aged pulsars, we show that ten-year HAWC observations should eventually observe 37$^{+17}_{-13}$ middle-aged TeV halos that correspond to pulsars whose radio emission is not beamed towards Earth. Depending on the extrapolation of the TeV halo efficiency to young pulsars, HAWC could detect more than 100 TeV halos from mis-aligned pulsars. These pulsars have historically been difficult to detect with existing multiwavelength observations. TeV halos will constitute a significant fraction of all HAWC sources, allowing follow-up observations to efficiently find pulsar wind nebulae and thermal pulsar emission. The observation and subsequent multi-wavelength follow-up of TeV halos will have significant implications for our understanding of pulsar beam geometries, the evolution of PWN, the diffusion of cosmic-rays near energetic pulsars, and the contribution of pulsars to the cosmic-ray positron excess.


title: A Systematic Survey of the Effects of Wind Mass Loss Algorithms on the Evolution of Single Massive Stars

abstract: Mass loss is a key uncertainty in the evolution of massive stars. Stellar evolution calculations must employ parametric algorithms for mass loss, and usually only include stellar winds. We carry out a parameter study of the effects of wind mass loss on massive star evolution using the open-source stellar evolution code MESA. We provide a systematic comparison of wind mass loss algorithms for solar-metallicity, nonrotating, single stars in the initial mass range of $15-35\,M_\odot$. We consider combinations drawn from two hot phase algorithms, three cool phase algorithms, and two Wolf-Rayet algorithms. We consider linear wind efficiency scale factors of $1$, $0.33$, and $0.1$ to account for reductions in mass loss rates due to wind inhomogeneities. We find that the initial to final mass mapping for each zero-age main-sequence (ZAMS) mass has a $\sim 50\%$ uncertainty if all algorithm combinations and wind efficiencies are considered. The ad-hoc efficiency scale factor dominates this uncertainty. While the final total mass and internal structure of our models vary tremendously with mass loss treatment, final observable parameters are much less sensitive for ZAMS mass $\lesssim 30\,M_\odot$. This indicates that uncertainty in wind mass loss does not negatively affect estimates of the ZAMS mass of most single-star supernova progenitors from pre-explosion observations. Furthermore, we show that the internal structure of presupernova stars is sensitive to variations in both main sequence and post main-sequence mass loss. We find that the compactness parameter $\xi\propto M/R(M)$ varies by as much as $30\%$ for a given ZAMS mass evolved with different wind efficiencies and mass loss algorithm combinations. [abridged]


title: Dissolution of topological Fermi arcs in a dirty Weyl semimetal

abstract: Weyl semimetals (WSMs) have recently attracted a great deal of attention as they provide condensed matter realization of chiral anomaly, feature topologically protected Fermi arc surface states and sustain sharp chiral Weyl quasiparticles up to a critical disorder at which a continuous quantum phase transition (QPT) drives the system into a metallic phase. We here numerically demonstrate that with increasing strength of disorder the Fermi arc gradually looses its sharpness, and close to the WSM-metal QPT it completely dissolves into the metallic bath of the bulk. Predicted topological nature of the WSM-metal QPT and the resulting bulk-boundary correspondence across this transition can directly be observed in angle-resolved-photo-emmision-spectroscopy (ARPES) and Fourier transformed scanning-tunneling-microscopy (STM) measurements by following the continuous deformation of the Fermi arcs with increasing disorder in recently discovered Weyl materials.

categories: cond-mat.mes-hall cond-mat.dis-nn cond-mat.str-el hep-th


title: Accelerating gravitational microlensing simulations using the Xeon Phi coprocessor

abstract: Recently Graphics Processing Units (GPUs) have been used to speed up very CPU-intensive gravitational microlensing simulations. In this work, we use the Xeon Phi coprocessor to accelerate such simulations and compare its performance on a microlensing code with that of NVIDIA's GPUs. For the selected set of parameters evaluated in our experiment, we find that the speedup by Intel's Knights Corner coprocessor is comparable to that by NVIDIA's Fermi family of GPUs with compute capability 2.0, but less significant than GPUs with higher compute capabilities such as the Kepler. However, the very recently released second generation Xeon Phi, Knights Landing, is about 5.8 times faster than the Knights Corner, and about 2.9 times faster than the Kepler GPU used in our simulations. We conclude that the Xeon Phi is a very promising alternative to GPUs for modern high performance microlensing simulations.


title: Effects of disc midplane evolution on CO snowline location

abstract: Temperature changes in the planet forming disc midplanes carry important physico-chemical consequences, such as the effect on the locations of the condensation fronts of molecules - the snowlines. Snowlines impose major chemical gradients and possibly foster grain growth. The aim of this paper is to understand how disc midplane temperature changes with gas and dust evolution, and identify trends that may influence planet formation or allow to constrain disc evolution observationally. We calculate disc temperature, hydrostatic equilibrium and dust settling in a mutually consistent way from a grid of disc models at different stages of gas loss, grain growth and hole opening. We find that the CO snowline location depends very strongly on disc properties. The CO snowline location migrates closer to the star for increasing degrees of gas dispersal and dust growth. Around a typical A type star, the snowline can be anywhere between several tens and a few hundred au, depending on the disc properties such as gas mass and grain size. In fact, gas loss is as efficient as dust evolution in settling discs, and flat discs may be gas-poor counterparts of flared discs. Our results, in the context of different pre-main sequence evolution of the luminosity in low- and intermediate-mass stars suggests very different thermal (and hence chemical) histories in these two types of discs. Discs of T Tauri stars settle and cool down while discs of Herbig Ae stars may remain rather warm throughout the pre-main sequence.


title: On the effective turbulence driving mode of molecular clouds formed in disc galaxies

abstract: We determine the physical properties and turbulence driving mode of molecular clouds formed in numerical simulations of a Milky Way-type disc galaxy with parsec-scale resolution. The clouds form through gravitational fragmentation of the gas, leading to average values for mass, radii and velocity dispersion in good agreement with observations of Milky Way clouds. The driving parameter (b) for the turbulence within each cloud is characterised by the ratio of the density contrast (sigma_rho) to the average Mach number (Mach) within the cloud, b = sigma_rho/Mach. As shown in previous works, b ~ 1/3 indicates solenoidal (divergence-free) driving and b ~ 1 indicates compressive (curl-free) driving. We find that the average b value of all the clouds formed in the simulations has a lower limit of b > 0.2. Importantly, we find that b has a broad distribution, covering values from purely solenoidal to purely compressive driving. Tracking the evolution of individual clouds reveals that the b value for each cloud does not vary significantly over their lifetime. Finally, we perform a resolution study with minimum cell sizes of 8, 4, 2 and 1 pc and find that the average b value increases with increasing resolution. Therefore, we conclude that our measured b values are strictly lower limits and that a resolution better than 1 pc is required for convergence. However, regardless of the resolution, we find that b varies by factors of a few in all cases, which means that the effective driving mode alters significantly from cloud to cloud.


title: Fast and scalable Gaussian process modeling with applications to astronomical time series

abstract: The growing field of large-scale time domain astronomy requires methods for probabilistic data analysis that are computationally tractable, even with large datasets. Gaussian Processes are a popular class of models used for this purpose but, since the computational cost scales, in general, as the cube of the number of data points, their application has been limited to small datasets. In this paper, we present a novel method for Gaussian Process modeling in one-dimension where the computational requirements scale linearly with the size of the dataset. We demonstrate the method by applying it to simulated and real astronomical time series datasets. These demonstrations are examples of probabilistic inference of stellar rotation periods, asteroseismic oscillation spectra, and transiting planet parameters. The method exploits structure in the problem when the covariance function is expressed as a mixture of complex exponentials, without requiring evenly spaced observations or uniform noise. This form of covariance arises naturally when the process is a mixture of stochastically-driven damped harmonic oscillators -- providing a physical motivation for and interpretation of this choice -- but we also demonstrate that it can be a useful effective model in some other cases. We present a mathematical description of the method and compare it to existing scalable Gaussian Process methods. The method is fast and interpretable, with a range of potential applications within astronomical data analysis and beyond. We provide well-tested and documented open-source implementations of this method in C++, Python, and Julia.

categories: astro-ph.IM astro-ph.EP astro-ph.SR physics.data-an stat.AP


title: Evidence for Two Hot Jupiter Formation Paths

abstract: Disk migration and high-eccentricity migration are two well-studied theories to explain the formation of hot Jupiters. The former predicts that these planets can migrate up until the planet-star Roche separation ($a_{Roche}$) and the latter predicts they will tidally circularize at a minimum distance of 2$a_{Roche}$. Considering long-running radial velocity and transit surveys have identified a couple hundred hot Jupiters to date, we can revisit the classic question of hot Jupiter formation in a data-driven manner. We approach this problem using data from several exoplanet surveys (radial velocity, Kepler, HAT, and WASP) allowing for either a single population or a mixture of populations associated with these formation channels, and applying a hierarchical Bayesian mixture model of truncated power laws of the form $x^{\gamma-1}$ to constrain the population-level parameters of interest (e.g., location of inner edges, $\gamma$, mixture fractions). Within the limitations of our chosen models, we find the current radial velocity and Kepler sample of hot Jupiters can be well explained with a single truncated power law distribution with a lower cutoff near 2$a_{Roche}$, a result that still holds after a decade, and $\gamma=-0.51\pm^{0.19}_{0.20}$. However, the HAT and WASP data show evidence for multiple populations (Bayes factor $\approx 10^{21}$). We find that $15\pm^{9}_{6}\%$ reside in a component consistent with disk migration ($\gamma=-0.04\pm^{0.53}_{1.27}$) and $85\pm^{6}_{9}\%$ in one consistent with high-eccentricity migration ($\gamma=-1.38\pm^{0.32}_{0.47}$). We find no immediately strong connections with some observed host star properties and speculate on how future exoplanet surveys could improve upon hot Jupiter population inference.


title: The splashback radius of halos from particle dynamics. I. The SPARTA algorithm

abstract: Motivated by the recent proposal of the splashback radius as a physical boundary of dark matter halos, we present a parallel computer code for Subhalo and PARticle Trajectory Analysis (SPARTA). The code analyzes the orbits of all simulation particles in all host halos, billions of orbits in the case of typical cosmological N-body simulations. Within this general framework, we develop an algorithm that accurately extracts the location of the first apocenter of particles after infall into a halo, or splashback. We define the splashback radius of a halo as the smoothed average of the apocenter radii of individual particles. This definition allows us to reliably measure the splashback radii of 95% of host halos above a resolution limit of 1000 particles. We show that, on average, the splashback radius and mass are converged to better than 5% accuracy with respect to mass resolution, snapshot spacing, and all free parameters of the method.


title: The VLA-COSMOS 3 GHz Large Project: Continuum data and source catalog release

abstract: We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys.


title: A catalogue of white dwarf candidates in VST ATLAS

abstract: The Sloan Digital Sky Survey (SDSS) has created a knowledge gap between the Northern and the Southern hemispheres which is very marked for white dwarfs: only $\simeq 15$% of the known white dwarfs are south of the equator. Here we make use of the VST ATLAS survey, one of the first surveys obtaining deep, optical, multi-band photometry over a large area of the southern skies, to remedy this situation. Applying the colour and proper-motion selection developed in Gentile Fusillo et al. (2015a) to the most recent internal data release (2016 April 25) of VST ATLAS we created a catalogue of $\simeq 4200$ moderately bright ($g\leq19$), high-confidence southern white dwarf candidates, which can be followed up individually with both the large array of southern telescopes or in bulk with forthcoming multi-object spectrographs.


title: A Nonlocal Approach to the Cosmological Constant Problem

abstract: We construct a model in which the cosmological constant is canceled from the gravitational equations of motion. Our model relies on two key ingredients: a nonlocal constraint on the action, which forces the spacetime average of the Lagrangian density to vanish, and a dynamical way for this condition to be satisfied classically with arbitrary matter content. We implement the former condition with a spatially-constant Lagrange multiplier associated with the volume form and the latter by including a free four-form gauge field strength in the action. These two features are enough to remove the cosmological constant from the Einstein equation. The model is consistent with all cosmological and experimental bounds on modification of gravity and allows for both cosmic inflation and the present epoch of acceleration.


title: The splashback radius of halos from particle dynamics. II. Dependence on mass, accretion rate, redshift, and cosmology

abstract: The splashback radius $R_{\rm sp}$, the apocentric radius of particles on their first orbit after falling into a dark matter halo, has recently been suggested as a physically motivated halo boundary that separates accreting from orbiting material. Using the SPARTA code presented in Paper I, we analyze the orbits of billions of particles in cosmological simulations of structure formation and measure $R_{\rm sp}$ for a large sample of halos that span a mass range from dwarf galaxy to massive cluster halos, reach redshift 8, and include WMAP, Planck, and self-similar cosmologies. We analyze the dependence of $R_{\rm sp}/R_{\rm 200m}$ and $M_{\rm sp}/M_{\rm 200m}$ on the mass accretion rate $\Gamma$, halo mass, redshift, and cosmology. The scatter in these relations varies between 0.02 and 0.1 dex. While we confirm the known trend that $R_{\rm sp}/R_{\rm 200m}$ decreases with $\Gamma$, the relationships turn out to be more complex than previously thought, demonstrating that $R_{\rm sp}$ is an independent definition of the halo boundary that cannot trivially be reconstructed from spherical overdensity definitions. We present fitting functions for $R_{\rm sp}/R_{\rm 200m}$ and $M_{\rm sp}/M_{\rm 200m}$ as a function of accretion rate, peak height, and redshift, achieving an accuracy of 5% or better everywhere in the parameter space explored. We discuss the physical meaning of the distribution of particle apocenters and show that the previously proposed definition of $R_{\rm sp}$ as the radius of the steepest logarithmic density slope encloses roughly three-quarters of the apocenters. Finally, we conclude that no analytical model presented thus far can fully explain our results.


title: The Calculation of Sparticle and Higgs Decays in the Minimal and Next-to-Minimal Supersymmetric Standard Models: SOFTSUSY4.0

abstract: We describe a major extension of the SOFTSUSY spectrum calculator to include the calculation of the decays, branching ratios and lifetimes of sparticles into lighter sparticles, covering the next-to-minimal supersymmetric standard model (NMSSM) as well as the minimal supersymmetric standard model (MSSM). This document acts as a manual for the new version of SOFTSUSY, which includes the calculation of sparticle decays. We present a comprehensive collection of explicit expressions used by the program for the various partial widths of the different decay modes in the appendix.


title: Possible $D^{(*)}\bar{D}^{(*)}$ and $B^{(*)}\bar{B}^{(*)}$ molecular states in the extended constituent quark models

abstract: The possible neutral $D^{(*)}\bar{D}^{(*)}$ and $B^{(*)}\bar{B}^{(*)}$ molecular states are studied in the framework of the constituent quark models, which is extended by including the $s$-channel one gluon exchange. Using different types of quark-quark potentials, we solve the four-body Schr\"{o}dinger equation by means of the Gaussian expansion method. The bound states of $D^{(*)}\bar{D}^{(*)}$ with $J^{PC}=1^{++},2^{++}$ and $B^{(*)}\bar{B}^{(*)}$ with $J^{PC}=0^{++},1^{+-},1^{++},2^{++}$ are obtained. The molecular states $D^{*}\bar{D}$ with $J^{PC}=1^{++}$ and $B^{*}\bar{B}$ with $J^{PC}=1^{+-}$ are good candidates for the $X(3872)$ and $Z^0_b(10610)$, respectively.


title: The VLA-COSMOS 3 GHz Large Project: Multiwavelength counterparts and the composition of the faint radio population

abstract: (abridged) We study the composition of the faint radio population selected from the VLA-COSMOS 3GHz Large Project. The survey covers a 2.6sq.deg. area with a mean rms of ~2.3uJy/b, cataloging 10830 sources (>5sigma). Combining these radio data with optical, near-infrared (UltraVISTA), mid-infrared (Spitzer/IRAC) data, and X-ray data (Chandra), we find counterparts to radio sources for ~93% of the radio sample (in the areas of the COSMOS field not affected by saturated or bright sources in the optical to NIR bands), reaching out to z<6. We further classify the sources as star forming galaxies or AGN based on various criteria, such as X-ray luminosity, observed MIR color, UV-FIR spectral-energy distribution, rest-frame NUV-optical color corrected for dust extinction, and radio-excess relative to that expected from the hosts' star-formation rate. We separate the AGN into sub-samples dominated by low-to-moderate and moderate-to-high radiative luminosity AGN, candidates for high-redshift analogues to local low- and high-excitation emission line AGN, respectively. We study the fractional contributions of these sub-populations down to radio flux levels of ~11uJy at 3GHz (or ~20uJy at 1.4GHz assuming a spectral index of -0.7). We find that the dominant fraction at 1.4GHz flux densities above ~200uJy is constituted of low-to-moderate radiative luminosity AGN. Below densities of ~100uJy the fraction of star-forming galaxies increases to ~60%, followed by the moderate-to-high radiative luminosity AGN (~20%), and low-to-moderate radiative luminosity AGN (~20%). Based on this observational evidence, we extrapolate the fractions down to sensitivities of the SKA. Our estimates suggest that at the faint flux limits to be reached by the SKA1 surveys, a selection based only on radio flux limits can provide a simple tool to efficiently identify samples highly (>75%) dominated by star-forming galaxies.


title: The VLA-COSMOS 3~GHz Large Project: AGN and host-galaxy properties out to z$\lesssim$6

abstract: We explore the multiwavelength properties of AGN host galaxies for different classes of radio-selected AGN out to z$\lesssim$6 via a multiwavelength analysis of about 7700 radio sources in the COSMOS field. The sources were selected with the Very Large Array (VLA) at 3 GHz (10 cm) within the VLA-COSMOS 3 GHz Large Project, and cross-matched with multiwavelength ancillary data. This is the largest sample of high-redshift (z$\lesssim$6) radio sources with exquisite photometric coverage and redshift measurements available. We constructed a sample of moderate-to-high radiative luminosity AGN (HLAGN) via spectral energy distribution (SED) decomposition combined with standard X-ray and mid-infrared diagnostics. Within the remainder of the sample we further identified low-to-moderate radiative luminosity AGN (MLAGN) via excess in radio emission relative to the star formation rates in their host galaxies. We show that AGN power in HLAGN occurs predominantly in radiative form, while MLAGN display a substantial mechanical AGN luminosity component. We found significant differences in the host properties of the two AGN classes, as a function of redshift. At z$<$1.5, MLAGN appear to reside in significantly more massive and less star-forming galaxies compared to HLAGN. At z$>$1.5, we observed a reversal in the behaviour of the stellar mass distributions with the HLAGN populating the higher stellar mass tail. We interpret this finding as a possible hint of the downsizing of galaxies hosting HLAGN, with the most massive galaxies triggering AGN activity earlier than less massive galaxies, and then fading to MLAGN at lower redshifts. Our conclusion is that HLAGN and MLAGN samples trace two distinct galaxy and AGN populations in a wide range of redshifts, possibly resembling the radio AGN types often referred to as radiative- and jet-mode (or high- and low-excitation), respectively.


title: The Galactic Contribution to IceCube's Astrophysical Neutrino Flux

abstract: High energy neutrinos have been detected by IceCube, but their origin remains a mystery. Determining the sources of this flux is a crucial first step towards multi-messenger studies. In this work we systematically compare two classes of sources with the data: Galactic and extragalactic. We assume that the neutrino sources are distributed according to a class of Galactic models. We build a likelihood function on an event by event basis including energy, event topology, absorption, and direction information. We present the probability that each high energy event with deposited energy $E_{\rm dep}>60$ TeV in the HESE sample is Galactic, extragalactic, or background. For Galactic models considered the Galactic fraction of the astrophysical flux has a best fit value of $1.3\%$ and is $<9.5\%$ at 90\% CL. A zero Galactic flux is allowed at $<1\sigma$.


title: Science with the space-based interferometer LISA. V: Extreme mass-ratio inspirals

abstract: The space-based Laser Interferometer Space Antenna (LISA) will be able to observe the gravitational-wave signals from systems comprised of a massive black hole and a stellar-mass compact object. These systems are known as extreme-mass-ratio inspirals (EMRIs) and are expected to complete $\sim 10^4$-$10^5$ cycles in band, thus allowing exquisite measurements of their parameters. In this work, we attempt to quantify the astrophysical uncertainties affecting the predictions for the number of EMRIs detectable by LISA, and find that competing astrophysical assumptions produce a variance of about three orders of magnitude in the expected intrinsic EMRI rate. However, we find that irrespective of the astrophysical model, at least a few EMRIs per year should be detectable by the LISA mission, with up to a few thousands per year under the most optimistic astrophysical assumptions. We also investigate the precision with which LISA will be able to extract the parameters of these sources. We find that typical fractional statistical errors with which the intrinsic parameters (redshifted masses, massive black hole spin and orbital eccentricity) can be recovered are $\sim 10^{-6}$-$10^{-4}$. Luminosity distance (which is required to infer true masses) is inferred to about $10\%$ precision and sky position is localized to a few square degrees, while tests of the multipolar structure of the Kerr metric can be performed to percent-level precision or better.


title: VLA-COSMOS 3GHz Large Project: The infrared-radio correlation of star-forming galaxies and AGN to $z\lesssim6$

abstract: We examine the behaviour of the infrared-radio correlation (IRRC) over the range $0<z<6$ using new, highly sensitive 3GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared data from the Herschel Space Observatory in the 2deg$^{2}$ COSMOS field. We distinguish between objects where emission is believed to arise solely from star-formation, and those where an active galactic nucleus (AGN) is thought to be present. We account for non-detections in the radio or in the infrared using a doubly-censored survival analysis. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4GHz radio luminosity ratio ($q_{\rm TIR}$), decreases with increasing redshift: $q_{\rm TIR}(z)=(2.88\pm0.03)(1+z)^{-0.19\pm0.01}$. Moderate-to-high radiative luminosity AGN do not follow the same $q_{\rm TIR}$$(z)$ trend, having a lower normalisation and steeper decrease with redshift. We cannot rule out the possibility that unidentified AGN contributions only to the radio regime may be steepening the observed $q_{\rm TIR}(z)$ trend of the star-forming population. An increasing fractional contribution to the observed 3GHz flux by free-free emission of star-forming galaxies may also affect the derived evolution. However, we find that the standard (M82-based) assumption of the typical radio spectral energy distribution (SED) for star-forming galaxies is inconsistent with our results. This suggests a more complex shape of the typical radio SED for star-forming galaxies, and that imperfect $K$ corrections in the radio may govern the derived redshift trend of $q_{\rm TIR}$. Lastly, we present a redshift-dependent relation between rest-frame 1.4GHz radio luminosity and star formation rate taking the derived redshift trend into account.


title: The VLA-COSMOS 3 GHz Large Project: The cosmic star formation history since z~5

abstract: We make use of the deep Karl G. Jansky Very Large Array (VLA) COSMOS radio observations at 3 GHz to infer radio luminosity functions of star-forming galaxies up to redshifts of z~5 based on approximately 6000 detections with reliable optical counterparts. This is currently the largest radio-selected sample available out to z~5 across an area of 2 square degrees with a sensitivity of rms=2.3 ujy/beam. By fixing the faint and bright end shape of the radio luminosity function to the local values, we find a strong redshift trend that can be fitted with a pure luminosity evolution L~(1+z)^{(3.16 +- 0.2)-(0.32 +- 0.07) z}. We estimate star formation rates (SFRs) from our radio luminosities using an infrared (IR)-radio correlation that is redshift dependent. By integrating the parametric fits of the evolved luminosity function we calculate the cosmic SFR density (SFRD) history since z~5. Our data suggest that the SFRD history peaks between 2<z<3 and that the ultraluminous infrared galaxies (ULIRGs; 100 Msol/yr<SFR<1000 Msol/yr) contribute up to ~25% to the total SFRD in the same redshift range. Hyperluminous infrared galaxies (HyLIRGs; SFR>1000 Msol/yr) contribute an additional <2% in the entire observed redshift range. We find evidence of a potential underestimation of SFRD based on ultraviolet (UV) rest-frame observations of Lyman break galaxies (LBGs) at high redshifts (z>4) on the order of 15-20%, owing to appreciable star formation in highly dust-obscured galaxies, which might remain undetected in such UV observations.


title: An Epipolar Line from a Single Pixel

abstract: Computing the epipolar geometry from feature points between cameras with very different viewpoints is often error prone, as an object's appearance can vary greatly between images. For such cases, it has been shown that using motion extracted from video can achieve much better results than using a static image. This paper extends these earlier works based on the scene dynamics. In this paper we propose a new method to compute the epipolar geometry from a video stream, by exploiting the following observation: For a pixel p in Image A, all pixels corresponding to p in Image B are on the same epipolar line. Equivalently, the image of the line going through camera A's center and p is an epipolar line in B. Therefore, when cameras A and B are synchronized, the momentary images of two objects projecting to the same pixel, p, in camera A at times t1 and t2, lie on an epipolar line in camera B. Based on this observation we achieve fast and precise computation of epipolar lines. Calibrating cameras based on our method of finding epipolar lines is much faster and more robust than previous methods.


title: Ruling out FPT algorithms for Weighted Coloring on forests

abstract: Given a graph $G$, a proper $k$-coloring of $G$ is a partition $c = (S_i)_{i\in [1,k]}$ of $V(G)$ into $k$ stable sets $S_1,\ldots, S_{k}$. Given a weight function $w: V(G) \to \mathbb{R}^+$, the weight of a color $S_i$ is defined as $w(i) = \max_{v \in S_i} w(v)$ and the weight of a coloring $c$ as $w(c) = \sum_{i=1}^{k}w(i)$. Guan and Zhu [Inf. Process. Lett., 1997] defined the weighted chromatic number of a pair $(G,w)$, denoted by $\sigma(G,w)$, as the minimum weight of a proper coloring of $G$. For a positive integer $r$, they also defined $\sigma(G,w;r)$ as the minimum of $w(c)$ among all proper $r$-colorings $c$ of $G$. The complexity of determining $\sigma(G,w)$ when $G$ is a tree was open for almost 20 years, until Ara\'ujo et al. [SIAM J. Discrete Math., 2014] recently proved that the problem cannot be solved in time $n^{o(\log n)}$ on $n$-vertex trees unless the Exponential Time Hypothesis (ETH) fails. The objective of this article is to provide hardness results for computing $\sigma(G,w)$ and $\sigma(G,w;r)$ when $G$ is a tree or a forest, relying on complexity assumptions weaker than the ETH. Namely, we study the problem from the viewpoint of parameterized complexity, and we assume the weaker hypothesis $FPT \neq W[1]$. Building on the techniques of Ara\'ujo et al., we prove that when $G$ is a forest, computing $\sigma(G,w)$ is $W[1]$-hard parameterized by the size of a largest connected component of $G$, and that computing $\sigma(G,w;r)$ is $W[2]$-hard parameterized by $r$. Our results rule out the existence of $FPT$ algorithms for computing these invariants on trees or forests for many natural choices of the parameter.


title: A Numerical Approach to Virasoro Blocks and the Information Paradox

abstract: We chart the breakdown of semiclassical gravity by analyzing the Virasoro conformal blocks to high numerical precision, focusing on the heavy-light limit corresponding to a light probe propagating in a BTZ black hole background. In the Lorentzian regime, we find empirically that the initial exponential time-dependence of the blocks transitions to a universal $t^{-\frac{3}{2}}$ power-law decay. For the vacuum block the transition occurs at $t \approx \frac{\pi c}{6 h_L}$, confirming analytic predictions. In the Euclidean regime, due to Stokes phenomena the naive semiclassical approximation fails completely in a finite region enclosing the `forbidden singularities'. We emphasize that limitations on the reconstruction of a local bulk should ultimately stem from distinctions between semiclassical and exact correlators.


title: $B_s \to D_s \ell \nu$ Form Factors and the Fragmentation Fraction Ratio $f_s/f_d$

abstract: We present a lattice quantum chromodynamics determination of the scalar and vector form factors for the $B_s \rightarrow D_s \ell \nu$ decay over the full physical range of momentum transfer. In conjunction with future experimental data, our results will provide a new method to extract $|V_{cb}|$, which may elucidate the current tension between exclusive and inclusive determinations of this parameter. Combining the form factor results at non-zero recoil with recent HPQCD results for the $B \rightarrow D \ell \nu$ form factors, we determine the ratios $f^{B_s \rightarrow D_s}_0(M_\pi^2) / f^{B \rightarrow D}_0(M_K^2) = 1.000(62)$ and $f^{B_s \rightarrow D_s}_0(M_\pi^2) / f^{B \rightarrow D}_0(M_\pi^2) = 1.006(62)$. These results give the fragmentation fraction ratios $f_s/f_d = 0.310(30)_{\mathrm{stat.}}(21)_{\mathrm{syst.}}(6)_{\mathrm{theor.}}(38)_{\mathrm{latt.}} $ and $f_s/f_d = 0.307(16)_{\mathrm{stat.}}(21)_{\mathrm{syst.}}(23)_{\mathrm{theor.}}(44)_{\mathrm{latt.}}$, respectively. The fragmentation fraction ratio is an important ingredient in experimental determinations of $B_s$ meson branching fractions at hadron colliders, in particular for the rare decay ${\cal B}(B_s \rightarrow \mu^+ \mu^-)$. In addition to the form factor results, we make the first prediction of the branching fraction ratio $R(D_s) = {\cal B}(B_s\to D_s\tau\nu)/{\cal B}(B_s\to D_s\ell\nu) = 0.301(6)$, where $\ell$ is an electron or muon. Current experimental measurements of the corresponding ratio for the semileptonic decays of $B$ mesons disagree with Standard Model expectations at the level of nearly four standard deviations. Future experimental measurements of $R(D_s)$ may help understand this discrepancy.


title: Flux Flattening in Axion Monodromy Inflation

abstract: String theory models of axion monodromy inflation exhibit scalar potentials which are quadratic for small values of the inflaton field and evolve to a more complicated function for large field values. Oftentimes the large field behaviour is gentler than quadratic, lowering the tensor-to-scalar ratio. This effect, known as flattening, has been observed in the string theory context through the properties of the DBI+CS D-brane action. We revisit such flattening effects in type IIB flux compactifications with mobile D7-branes, with the inflaton identified with the D7-brane position. We observe that, with a generic choice of background fluxes, flattening effects are larger than previously observed, allowing to fit these models within current experimental bounds. In particular, we compute the cosmological observables in scenarios compatible with closed-string moduli stabilisation, finding tensor-to-scalar ratios as low as r ~ 0.04. These are models of single field inflation in which the inflaton is much lighter than the other scalars through a mild tuning of the compactification data.


title: The Discovery and Origin of A Very-High Velocity Cloud Toward M33

abstract: We report the detection of a largely ionized very-high velocity cloud (VHVC; $v_{\rm LSR}\sim-350$ km/s) toward M33 with the Hubble Space Telescope/Cosmic Origin Spectrograph. The VHVC is detected in OI, CII, SiII, and SiIII absorption along five sightlines separated by ~0.06-0.4 degree. On sub-degree scales, the velocities and ionic column densities of the VHVC remain relatively smooth with standard deviations of +/-14 km/s and +/-0.15 dex between the sightlines, respectively. The VHVC has a metallicity of [OI/HI]=-0.56+/-0.17 dex (Z=0.28+/-0.11 Z$_{\odot}$). Despite the position-velocity proximity of the VHVC to the ionized Magellanic Stream, the VHVC's higher metallicity makes it unlikely to be associated with the Stream, highlighting the complex velocity structure of this region of sky. We investigate the VHVC's possible origin by revisiting its surrounding HI environment. We find that the VHVC may be: (1) a MW CGM cloud, (2) related to a nearby HI VHVC -- Wright's Cloud, or (3) connected to M33's northern warp. Furthermore, the VHVC could be a bridge connecting Wright's Cloud and M33's northern warp, which would make it a Magellanic-like structure in the halo of M33.


title: Survival asymptotics for branching random walks in IID environments

abstract: We first study a model, introduced recently in \cite{ES}, of a critical branching random walk in an IID random environment on the $d$-dimensional integer lattice. The walker performs critical (0-2) branching at a lattice point if and only if there is no `obstacle' placed there. The obstacles appear at each site with probability $p\in [0,1)$ independently of each other. We also consider a similar model, where the offspring distribution is subcritical. Let $S_n$ be the event of survival up to time $n$. We show that on a set of full $\mathbb P_p$-measure, as $n\to\infty$, (i) Critical case: P^{\omega}(S_n)\sim\frac{2}{qn}; (ii) Subcritical case: P^{\omega}(S_n)= \exp\left[\left( -C_{d,q}\cdot \frac{n}{(\log n)^{2/d}} \right)(1+o(1))\right], where $C_{d,q}>0$ does not depend on the branching law. Hence, the model exhibits `self-averaging' in the critical case but not in the subcritical one. I.e., in (i) the asymptotic tail behavior is the same as in a "toy model" where space is removed, while in (ii) the spatial survival probability is larger than in the corresponding toy model, suggesting spatial strategies. We utilize a spine decomposition of the branching process as well as some known results on random walks.


title: The Colin de Verdi\`ere parameter, excluded minors, and the spectral radius

abstract: In this paper we characterize graphs which maximize the spectral radius of their adjacency matrix over all graphs of Colin de Verdi\`ere parameter at most $m$. We also characterize graphs of maximum spectral radius with no $H$ as a minor when $H$ is either $K_r$ or $K_{s,t}$. Interestingly, the extremal graphs match those which maximize the number of edges over all graphs with no $H$ as a minor when $r$ and $s$ are small, but not when they are larger.


title: Modeling the Vertical Structure of Nuclear Starburst Discs: A Possible Source of AGN Obscuration at $z\sim 1$

abstract: Nuclear starburst discs (NSDs) are star-forming discs that may be residing in the nuclear regions of active galaxies at intermediate redshifts. One dimensional (1D) analytical models developed by Thompson et al. (2005) show that these discs can possess an inflationary atmosphere when dust is sublimated on parsec scales. This make NSDs a viable source for AGN obscuration. We model the two dimensional (2D) structure of NSDs using an iterative method in order to compute the explicit vertical solutions for a given annulus. These solutions satisfy energy and hydrostatic balance, as well as the radiative transfer equation. In comparison to the 1D model, the 2D calculation predicts a less extensive expansion of the atmosphere by orders of magnitude at the parsec/sub-parsec scale, but the new scale-height $h$ may still exceed the radial distance $R$ for various physical conditions. A total of 192 NSD models are computed across the input parameter space in order to predict distributions of a line of sight column density $N_H$. Assuming a random distribution of input parameters, the statistics yield 56% of Type 1, 23% of Compton-thin Type 2s (CN), and 21% of Compton-thick (CK) AGNs. Depending on a viewing angle ($\theta$) of a particular NSD (fixed physical conditions), any central AGN can appear to be Type 1, CN, or CK which is consistent with the basic unification theory of AGNs. Our results show that $\log[N_H(\text{cm}^{-2})]\in$ [23,25.5] can be oriented at any $\theta$ from 0$^\circ$ to $\approx$80$^\circ$ due to the degeneracy in the input parameters.


title: Extending functions from nonisotropic Nikolskii-Besov spaces and approximating their derivatives (with supplement)

abstract: The article examines nonisotropic Nikolskii and Besov spaces with norms defined using $L_p$-averaged moduli of continuity of functions of appropriate orders along the coordinate directions, instead of moduli of continuity of known orders for derivative functions along the same directions. The author builds continuous linear mappings of such spaces of functions defined in domains of certain type to ordinary nonisotropic Nikolskii and Besov spaces in $\mathbb{R}^d $ that are function extension operators, thus incurring coincidence of both kinds of spaces in the said domains. The article also provides weak asymptotics of approximation characteristics related to the problem of derivative reconstruction from function values at a given number of points, the S.B.Stechkin's problem for differential operator, and the problem of width asymptotics for nonisotropic Nikolskii-Besov classes in those domains. The supplement provides a broad-sense definition of $\alpha$-type domain. It is demonstrated that all of the earlier statements established by the author with respect to the tasks of function extension, derivative reconstruction from function values, differential operator approximation with bounded operators, and width asymptotics for nonisotropic Nikolskii-Besov function classes defined in $\alpha$-type domains (in a narrow sense), remain valid in the same formulation even for an $\alpha$-type domain as construed in a broad sense.


title: Theory of electron cooling

abstract: We consider the process of cooling of a heavy particle beam in a co-moving electron beam of low temperature guided by a solenoidal magnetic field. This paper summarizes the main results of theoretical studies of this process conducted by the author during a period of several years. The main result of these studies is a conclusion that magnetization of the electron beam can provide the possibility of drastic enhancement of the cooling rate of a heavy particle beam with achieving equilibrium temperatures that are much lower than the transverse temperatures of the electron beam. Magnetized electron cooling is proposed for cooling of ion beams in high-luminosity colliders.


title: Graphene-based Josephson junction single photon detector

abstract: We propose to use graphene-based Josephson junctions (gJjs) to detect single photons in a wide electromagnetic spectrum from visible to radio frequencies. Our approach takes advantage of the exceptionally low electronic heat capacity of monolayer graphene and its constricted thermal conductance to its phonon degrees of freedom. Such a system could provide high sensitivity photon detection required for research areas including quantum information processing and radio-astronomy. As an example, we present our device concepts for gJj single photon detectors in both the microwave and infrared regimes. The dark count rate and intrinsic quantum efficiency are computed based on parameters from a measured gJj, demonstrating feasibility within existing technologies.


title: Inverting the Sweep Map on (kn,n)-Dyck Paths: A Simple Algorithm

abstract: We show that our algorithm for inverting the sweep map on (2n, n)-Dyck paths works for any (kn, n)-Dyck path, where k is an arbitrary positive integer.


title: Symbolic and Numerical Analysis in General Relativity with Open Source Computer Algebra Systems

abstract: We study three computer algebra systems, namely SageMath (with SageManifolds package), Maxima (with ctensor package) and Python language (with GraviPy module), which allow tensor manipulation for general relativity calculations along with general algebraic calculations. We present a benchmark of these systems using simple examples. After the general analysis, we focus on the SageMath and SageManifolds system to derive, analyze and visualize the solutions of the massless Klein-Gordon equation and geodesic motion with Hamilton-Jacobi formalism. We compare our numerical result of the Klein-Gordon equation with the asymptotic form of the analytical solution to see that they agree.


title: Quantum-Electron Back Action on Hybridization of Radiative and Evanescent Field Modes

abstract: A back action from Dirac electrons in graphene on the hybridization of radiative and evanescent fields is found as an analogy to Newton's third law. Here, the back action appears as a localized polarization field which greatly modifies an incident surface-plasmon-polariton (SPP) field. This yields a high sensitivity to local dielectric environments and provides a scrutiny tool for molecules or proteins selectively bounded with carbons. A scattering matrix is shown with varied frequencies nearby the surface-plasmon (SP) resonance for the increase, decrease and even a full suppression of the polarization field, which enables accurate effective-medium theories to be constructed for Maxwell-equation finite-difference time-domain methods. Moreover, double peaks in the absorption spectra for hybrid SP and graphene-plasmon modes are significant only with a large conductor plasma frequency, but are overshadowed by a round SPP peak at a small plasma frequency as the graphene is placed close to conductor surface. These resonant absorptions facilitate the polariton-only excitations, leading to polariton condensation for a threshold-free laser.


title: Quantum quench dynamics

abstract: Quench dynamics is an active area of study encompassing condensed matter physics and quantum information, with applications to cold-atomic gases and pump-probe spectroscopy of materials. Recent theoretical progress in studying quantum quenches is reviewed. Quenches in interacting one dimensional systems as well as systems in higher spatial dimensions are covered. The appearance of non-trivial steady states following a quench in exactly solvable models is discussed, and the stability of these states to perturbations is described. Proper conserving approximations needed to capture the onset of thermalization at long times are outlined. The appearance of universal scaling for quenches near critical points, and the role of the renormalization group in capturing the transient regime, are reviewed. Finally the effect of quenches near critical points on the dynamics of entanglement entropy and entanglement statistics is discussed. The extraction of critical exponents from the entanglement statistics is outlined.


title: The Origin of Planetary Ring Systems

abstract: This chapter of the book Planetary Ring Systems addresses the origin of planetary rings, one of the least understood processes related to planet formation and evolution. Whereas rings seem ubiquitous around giant planets, their great diversity of mass, structure and composition is a challenge for any formation scenario. Recent advances in our understanding of ring and satellite formation and destruction suggest that these processes are interconnected, so that rings and satellites may be two aspects of the same geological system. However, no single theory seems able to explain the origin of the different planetary rings known in our Solar System, and it now seems evident that rings may result from a variety of processes like giant collisions, tidal stripping of comets or satellites, as well as planet formation itself. In order to build any theory of ring formation it is important to specify physical processes that affect the long-term evolution of rings, as well as to describe the different observations that any ring formation model should explain. This is the topic of section 2. In section 3, we focus our attention on Saturn's rings and their main properties, and then discuss the pros and cons of a series of ring formation models. We also discuss the link between rings and satellites. In section 4, we extend the discussion to the other giant planets (Jupiter, Uranus, and Neptune). Section 5 is devoted to new types of rings -- the recent discovery of rings orbiting small outer Solar System bodies (Centaurs), and the possible rings around extrasolar planets. In section 6, we conclude and try to identify critical observations and theoretical advances needed to better understand the origin of rings and their significance in the global evolution of planets.


title: Isobaric Reconstruction of the Baryonic Acoustic Oscillation

abstract: In this paper, we report a significant recovery of the linear baryonic acoustic oscillation (BAO) signature by applying the isobaric reconstruction algorithm to the non-linear matter density field. Assuming only the longitudinal component of the displacement being cosmologically relevant, this algorithm iteratively solves the coordinate transform between the Lagrangian and Eulerian frames without requiring any specific knowledge of the dynamics. For dark matter field, it produces the non-linear displacement potential with very high fidelity. The reconstruction error at the pixel level is within a few percent, and is caused only by the emergence of the transverse component after the shell-crossing. As it circumvents the strongest non-linearity of the density evolution, the reconstructed field is well-described by linear theory and immune from the bulk-flow smearing of the BAO signature. Therefore this algorithm could significantly improve the measurement accuracy of the sound horizon scale. For a perfect large-scale structure survey at redshift zero without Poisson or instrumental noise, the fractional error is reduced by a factor of 2.7, very close to the ideal limit with linear power spectrum and Gaussian covariance matrix.


title: A holographic stress-energy tensor near the Cauchy horizon inside a rotating black hole

abstract: We investigate a stress-energy tensor for a CFT at strong coupling inside a small five-dimensional rotating Myers-Perry black hole with equal angular momenta by using the holographic method. As a gravitational dual, we perturbatively construct a black droplet solution by applying the "derivative expansion" method, generalizing the work of Haddad (arXiv:1207.2305), and analytically compute the holographic stress-energy tensor for our solution. We find that the stress-energy tensor is finite at both the future and past outer (event) horizons, and that the energy density is negative just outside the event horizons due to the Hawking effect. Furthermore, we apply the holographic method to the question of quantum instability of the Cauchy horizon since, by construction, our black droplet solution also admits a Cauchy horizon inside. We analytically show that the null-null component of the holographic stress-energy tensor negatively diverges at the Cauchy horizon, suggesting that a singularity appears there, in favor of strong cosmic censorship.


title: Feature Analysis and Selection for Training an End-to-End Autonomous Vehicle Controller Using the Deep Learning Approach

abstract: Deep learning-based approaches have been widely used for training controllers for autonomous vehicles due to their powerful ability to approximate nonlinear functions or policies. However, the training process usually requires large labeled data sets and takes a lot of time. In this paper, we analyze the influences of features on the performance of controllers trained using the convolutional neural networks (CNNs), which gives a guideline of feature selection to reduce computation cost. We collect a large set of data using The Open Racing Car Simulator (TORCS) and classify the image features into three categories (sky-related, roadside-related, and road-related features).We then design two experimental frameworks to investigate the importance of each single feature for training a CNN controller.The first framework uses the training data with all three features included to train a controller, which is then tested with data that has one feature removed to evaluate the feature's effects. The second framework is trained with the data that has one feature excluded, while all three features are included in the test data. Different driving scenarios are selected to test and analyze the trained controllers using the two experimental frameworks. The experiment results show that (1) the road-related features are indispensable for training the controller, (2) the roadside-related features are useful to improve the generalizability of the controller to scenarios with complicated roadside information, and (3) the sky-related features have limited contribution to train an end-to-end autonomous vehicle controller.


title: Profiling Users by Modeling Web Transactions

abstract: Users of electronic devices, e.g., laptop, smartphone, etc. have characteristic behaviors while surfing the Web. Profiling this behavior can help identify the person using a given device. In this paper, we introduce a technique to profile users based on their web transactions. We compute several features extracted from a sequence of web transactions and use them with one-class classification techniques to profile a user. We assess the efficacy and speed of our method at differentiating 25 users on a dataset representing 6 months of web traffic monitoring from a small company network.


title: Coordinating Filters for Faster Deep Neural Networks

abstract: Very large-scale Deep Neural Networks (DNNs) have achieved remarkable successes in a large variety of computer vision tasks. However, the high computation intensity of DNNs makes it challenging to deploy these models on resource-limited systems. Some studies used low-rank approaches that approximate the filters by low-rank basis to accelerate the testing. Those works directly decomposed the pre-trained DNNs by Low-Rank Approximations (LRA). How to train DNNs toward lower-rank space for more efficient DNNs, however, remains as an open area. To solve the issue, in this work, we propose Force Regularization, which uses attractive forces to enforce filters so as to coordinate more weight information into lower-rank space. We mathematically and empirically verify that after applying our technique, standard LRA methods can reconstruct filters using much lower basis and thus result in faster DNNs. The effectiveness of our approach is comprehensively evaluated in ResNets, AlexNet, and GoogLeNet. In AlexNet, for example, Force Regularization gains 2x speedup on modern GPU without accuracy loss and 4.05x speedup on CPU by paying small accuracy degradation. Moreover, Force Regularization better initializes the low-rank DNNs such that the fine-tuning can converge faster toward higher accuracy. The obtained lower-rank DNNs can be further sparsified, proving that Force Regularization can be integrated with state-of-the-art sparsity-based acceleration methods. Source code is available in https://github.com/wenwei202/caffe


title: Randomized Benchmarking, Correlated Noise, and Ising Models

abstract: We compute the expected randomized benchmarking sequence fidelity for a system subject to Gaussian time-correlated noise. For single qubit benchmarking we show that the expected sequence fidelity is given by the partition function of a long-range coupled spin-one Ising model, with each site in the Ising model corresponding to a free evolution interval. For d-state systems, the expected sequence fidelity is given by an Ising-like model partition function whose site variables are given by the weights of the adjoint representation of SU(d). A high effective temperature expansion for the partition function in the single qubit case shows decay of sequence fidelity varying from exponential for uncorrelated noise to a power law for quasistatic noise. Fitting an exponential to the sequence fidelity decay under correlated noise gives unreliable estimates of the average gate error rate.


title: Smallest order closed sublattices and option spanning

abstract: Let $Y$ be a sublattice of a vector lattice $X$. We consider the problem of identifying the smallest order closed sublattice of $X$ containing $Y$. It is known that the analogy with topological closure fails. Let $\overline{Y}^o$ be the order closure of $Y$ consisting of all order limits of nets of elements from $Y$. Then $\overline{Y}^o$ need not be order closed. We show that in many cases the smallest order closed sublattice containing $Y$ is in fact the second order closure $\overline{\overline{Y}^o}^o$. Moreover, if $X$ is a $\sigma$-order complete Banach lattice, then the condition that $\overline{Y}^o$ is order closed for every sublattice $Y$ characterizes order continuity of the norm of $X$. The present paper provides a general approach to a fundamental result in financial economics concerning the spanning power of options written on a financial asset.


title: Developpement de Methodes Automatiques pour la Reutilisation des Composants Logiciels

abstract: The large amount of information and the increasing complexity of applications constrain developers to have stand-alone and reusable components from libraries and component markets.Our approach consists in developing methods to evaluate the quality of the software component of these libraries, on the one hand and moreover to optimize the financial cost and the adaptation's time of these selected components. Our objective function defines a metric that maximizes the value of the software component quality by minimizing the financial cost and maintenance time. This model should make it possible to classify the components and order them in order to choose the most optimized. MOTS-CLES : d{\'e}veloppement de m{\'e}thode, r{\'e}utilisation, composants logiciels, qualit{\'e} de composant KEYWORDS:method development, reuse, software components, component quality .


title: Max Dehn, Axel Thue, and the Undecidable

abstract: This is a short essay on the roles of Max Dehn and Axel Thue in the formulation of the word problem for (semi)groups, and the story of the proofs showing that the word problem is undecidable.


title: Fourth order Superintegrable systems separating in Cartesian coordinates I. Exotic quantum potentials

abstract: A study is presented of two-dimensional superintegrable systems separating in Cartesian coordinates and allowing an integral of motion that is a fourth order polynomial in the momenta. All quantum mechanical potentials that do not satisfy any linear differential equation are found. They do however satisfy nonlinear ODEs. We show that these equations always have the Painlev\'e property and integrate them in terms of known Painlev\'e transcendents or elliptic functions.


title: Collective Anomaly Detection based on Long Short Term Memory Recurrent Neural Network

abstract: Intrusion detection for computer network systems becomes one of the most critical tasks for network administrators today. It has an important role for organizations, governments and our society due to its valuable resources on computer networks. Traditional misuse detection strategies are unable to detect new and unknown intrusion. Besides, anomaly detection in network security is aim to distinguish between illegal or malicious events and normal behavior of network systems. Anomaly detection can be considered as a classification problem where it builds models of normal network behavior, which it uses to detect new patterns that significantly deviate from the model. Most of the cur- rent research on anomaly detection is based on the learning of normally and anomaly behaviors. They do not take into account the previous, re- cent events to detect the new incoming one. In this paper, we propose a real time collective anomaly detection model based on neural network learning and feature operating. Normally a Long Short Term Memory Recurrent Neural Network (LSTM RNN) is trained only on normal data and it is capable of predicting several time steps ahead of an input. In our approach, a LSTM RNN is trained with normal time series data before performing a live prediction for each time step. Instead of considering each time step separately, the observation of prediction errors from a certain number of time steps is now proposed as a new idea for detecting collective anomalies. The prediction errors from a number of the latest time steps above a threshold will indicate a collective anomaly. The model is built on a time series version of the KDD 1999 dataset. The experiments demonstrate that it is possible to offer reliable and efficient for collective anomaly detection.


title: Self semi conjugations of Ulam's Tent-map

abstract: We study the self-semiconjugations of the Tent-map $f:\, x\mapsto 1-|2x-1|$ for $x\in [0,\, 1]$. We prove that each of these semi-conjugations $\xi$ is piecewise linear. For any $n\in \mathbb{N}$ we denote $A_n = f^{-n}(0)$ and describe the maps $\psi:\, A_n\rightarrow [0,\, 1]$ such that $\psi\circ f = f\circ \psi$. Also we describe all possible restrictions, of self-semiconjugations of the Tent-map onto $A_n$ and prove that for any $\alpha\in A_n\setminus A_{n-1}$ a restriction is completely determined by its value at $\alpha$.