title: Analyzing gender inequality through large-scale Facebook advertising data

abstract: Online social media are information resources that can have a transformative power in society. While the Web was envisioned as an equalizing force that allows everyone to access information, the digital divide prevents large amounts of people from being present online. Online social media in particular are prone to gender inequality, an important issue given the link between social media use and employment. Understanding gender inequality in social media is a challenging task due to the necessity of data sources that can provide large-scale measurements across multiple countries. Here we show how the Facebook Gender Divide (FGD), a metric based on aggregated statistics of more than 1.4 Billion users in 217 countries, explains various aspects of worldwide gender inequality. Our analysis shows that the FGD encodes gender equality indices in education, health, and economic opportunity. We find gender differences in network externalities that suggest that using social media has an added value for women. Furthermore, we find that low values of the FGD are associated with increases in economic gender equality. Our results suggest that online social networks, while suffering evident gender imbalance, may lower the barriers that women have to access informational resources and help to narrow the economic gender gap.


title: Linear response for random dynamical systems

abstract: We study for the first time linear response for random compositions of maps, chosen independently according to a distribution $\PP$. We are interested in the following question: how does an absolutely continuous stationary measure (acsm) of a random system change when $\PP$ changes smoothly to $\PP_{\eps}$? For a wide class of one dimensional random maps, we prove differentiability of acsm with respect to $\eps$; moreover, we obtain a linear response formula. We apply our results to iid compositions, with respect to various distributions $\PP_{\eps}$, of uniformly expanding circle maps, Gauss-R\'enyi maps (random continued fractions) and Pomeau-Manneville maps. Our results yield an exact formula for the invariant density of random continued fractions; while for Pomeau-Manneville maps our results provide a precise relation between their linear response under certain random perturbations and their linear response under deterministic perturbations.


title: Lamplighters admit weakly aperiodic SFTs

abstract: Let $A$ be a finite set and $G$ a group. A closed subset $X$ of $A^G$ is called a subshift if the action of $G$ on $A^G$ preserves $X$. If $K$ is a closed subset of $A^G$ such that membership in $K$ is determined by looking at a fixed finite set of coordinates, and $X$ is the intersection of all translates of $K$ under the action of $G$, then $X$ is called a subshift of finite type (SFT). If an SFT is nonempty and contains no finite $G$-orbits, it is said to be weakly aperiodic. A virtually cyclic group has no weakly aperiodic SFT, and Carroll and Penland have conjectured that a group with no weakly aperiodic SFT must be virtually cyclic. Answering a question of Jeandel, we show that lamplighters always admit weakly aperiodic SFTs.


title: On the (In)stability of the Identity Map in Optimal Transportation

abstract: We collect some examples of optimal transports in order to explore the (in)stability of the identity map as an optimal transport. First, we consider density and domain perturbations near regular portions of domains. Second, we investigate density and domains deformations to non-regular parts of domains. Here, we restrict our attention to two dimensions and focus near 90 degree corners.


title: Strong electromagnetic pulses generated in high-intensity short-pulse laser interactions with thin foil targets

abstract: Measurements are reported of the target neutralization current, the target charge, and the tangential component of the magnetic field generated as a result of laser-target interaction by pulses with the energy in the range of 45 mJ to 92 mJ on target and the pulse duration from 39 fs to 1000 fs. The experiment was performed at the Eclipse facility in CELIA, Bordeaux. The aim of the experiment was to extend investigations performed for the thick (mm scale) targets to the case of thin (micrometer thickness) targets in a way that would allow for a straightforward comparison of the results. We found that thin foil targets tend to generate 20 to 50 percent higher neutralization current and the target charge than the thick targets. The measurement of the tangential component of the magnetic field had shown that the initial spike is dominated by the 1 ns pulse consistent with the 1 ns pulse of the neutralization current, but there are some differences between targets of different type on sub-ns scale, which is an effect going beyond a simple picture of the target acting as an antenna. The sub-ns structure appears to be reproducible to surprising degree. We found that there is in general a linear correlation between the maximum value of the magnetic field and the maximum neutralization current, which supports the target-antenna picture, except for pulses hundreds of fs long.


title: LaSalle Invariance Principle for Discrete-time Dynamical Systems: A Concise and Self-contained Tutorial

abstract: LaSalle invariance principle was originally proposed in the 1950's and has become a fundamental mathematical tool in the area of dynamical systems and control. In both theoretical research and engineering practice, discrete-time dynamical systems have been at least as extensively studied as continuous-time systems. For example, model predictive control is typically studied in discrete-time via Lyapunov methods. However, there is a peculiar absence in the standard literature of standard treatments of Lyapunov functions and LaSalle invariance principle for discrete-time nonlinear systems. Most of the textbooks on nonlinear dynamical systems focus only on continuous-time systems. In Chapter 1 of the book by LaSalle [11], the author establishes the LaSalle invariance principle for difference equation systems. However, all the useful lemmas in [11] are given in the form of exercises with no proof provided. In this document, we provide the proofs of all the lemmas proposed in [11] that are needed to derive the main theorem on the LaSalle invariance principle for discrete-time dynamical systems. We organize all the materials in a self-contained manner. We first introduce some basic concepts and definitions in Section 1, such as dynamical systems, invariant sets, and limit sets. In Section 2 we present and prove some useful lemmas on the properties of invariant sets and limit sets. Finally, we establish the original LaSalle invariance principle for discrete-time dynamical systems and a simple extension in Section~3. In Section 4, we provide some references on extensions of LaSalle invariance principles for further reading. This document is intended for educational and tutorial purposes and contains lemmas that might be useful as a reference for researchers.


title: On the Calculation of Lorenz Numbers for Complex Thermoelectric Materials

abstract: A first-principles informed approach to the calculation of Lorenz numbers for complex thermoelectric materials is presented and discussed. Example calculations illustrate the importance of using accurate band structures and energy-dependent scattering times. Results obtained by assuming that the scattering rate follows the density-of-states show that in the non-degenerate limit, Lorenz numbers below the commonly assumed lower limit of 2(kB/q)^2 can occur. The physical cause of low Lorenz numbers is explained by the shape of the transport distribution. The numerical and physical issues that need to be addressed in order to produce accurate calculations of the Lorenz number are identified. The results of this study provide a general method that should contribute to the interpretation of measurements of total thermal conductivity and to the search for materials with low Lorenz numbers, which may provide improved thermoelectric figures of merit, zT.


title: On the $\mathbf{\rm\Psi}-$fractional integral and applications

abstract: Motivated by the ${\rm \Psi}$-Riemann-Liouville $({\rm \Psi-RL})$ fractional derivative and by the ${\rm \Psi}$-Hilfer $({\rm \Psi-H})$ fractional derivative, we introduced a new fractional operator the so-called $\rm\Psi-$fractional integral. We present some important results by means of theorems and in particular, that the $\rm\Psi-$fractional integration operator is limited. In this sense, we discuss some examples, in particular, involving the Mittag-Leffler $({\rm M-L})$ function, of paramount importance in the solution of population growth problem, as approached. On the other hand, we realize a brief discussion on the uniqueness of nonlinear $\Psi$-fractional Volterra integral equation (${\rm VIE}$) using $\beta-$distance functions.


title: All Killing Superalgebras for Warped AdS Backgrounds

abstract: We present all the symmetry superalgebras $\mathfrak{g}$ of all warped AdS$_k\times_w M^{d-k}$, $k>2$, flux backgrounds in $d=10, 11$ dimensions preserving any number of supersymmetries. First we give the conditions for $\mathfrak{g}$ to decompose into a direct sum of the isometry algebra of AdS$_k$ and that of the internal space $M^{d-k}$. Assuming this decomposition, we identify all symmetry superalgebras of AdS$_3$ backgrounds by showing that the isometry groups of internal spaces act transitively on spheres. We demonstrate that in type II and $d=11$ theories the AdS$_3$ symmetry superalgebras may not be simple and also present all symmetry superalgebras of heterotic AdS$_3$ backgrounds. Furthermore, we explicitly give the symmetry superalgebras of AdS$_k$, $k>3$, backgrounds and prove that they are all classical.


title: A study of accurate exchange-correlation functionals through adiabatic connection

abstract: A systematic way of improving exchange-correlation energy functionals of density functional theory has been to make them satisfy more and more exact relations. Starting from the initial GGA functionals, this has culminated into the recently proposed SCAN(Strongly constrained and appropriately normed) functional that satisfies several known constraints and is appropriately normed. The ultimate test for the functionals developed is the accuracy of energy calculated by employing them. In this paper, we test these exchange-correlation functionals -the GGA hybrid functionals B3LYP and PBE0, and the meta-GGA functional SCAN- from a different perspective. We study how accurately these functionals reproduce the exchange-correlation energy when electron-electron interaction is scaled as scaling parameter times Vee with this parameter varying between 0 and 1. Our study reveals interesting comparison between these functionals and the associated difference Tc between the interacting and the non-interacting kinetic energy for the same density.


title: Salient Signature of van der Waals Interactions

abstract: van der Waals interactions govern the physics of a plethora of molecular structures. It is well known that the leading term in the distance-based London expansion of the van der Waals energy for atomic and molecular dimers decays as $1/R^6$, where $R$ is the dimer distance. Using perturbation theory, we find the leading term in the distance-based expansion of the intracule pair density at the interatomic distance. Our results unveil a universal $1/R^3$ decay, which is less prone to numerical errors than the $1/R^6$ dependency, and it is confirmed numerically in H$_2$ and He$_2$ molecules. This \textit{signature} of van der Waals interactions can be directly used in the construction of approximate pair density and energy functionals including vdW corrections.


title: Multiwavelength observations of nova SMCN 2016-10a --- one of the brightest novae ever observed

abstract: We report on multiwavelength observations of nova SMCN 2016-10a. The present observational set is one of the most comprehensive for any nova in the Small Magellanic Cloud, including: low, medium, and high resolution optical spectroscopy and spectropolarimetry from SALT, FLOYDS, and SOAR; long-term OGLE $V$- and $I$- bands photometry dating back to six years before eruption; SMARTS optical and near-IR photometry from $\sim$ 11 days until over 280 days post-eruption; $Swift$ satellite X-ray and ultraviolet observations from $\sim$ 6 days until 319 days post-eruption. The progenitor system contains a bright disk and a main sequence or a sub-giant secondary. The nova is very fast with $t_2 \simeq$ 4.0 $\pm$ 1.0 d and $t_3 \simeq$ 7.8 $\pm$ 2.0 d in the $V$-band. If the nova is in the SMC, at a distance of $\sim$ 61 $\pm$ 10 kpc, we derive $M_{V,\mathrm{max}} \simeq - 10.5$ $\pm$ 0.5, making it the brightest nova ever discovered in the SMC and one of the brightest on record. At day 5 post-eruption the spectral lines show a He/N spectroscopic class and a FWHM of $\sim$ 3500 kms$^{-1}$ indicating moderately high ejection velocities. The nova entered the nebular phase $\sim$ 20 days post-eruption, predicting the imminent super-soft source turn-on in the X-rays, which started $\sim$ 28 days post-eruption. The super-soft source properties indicate a white dwarf mass between 1.2 M$_{\odot}$ and 1.3 M$_{\odot}$ in good agreement with the optical conclusions.


title: The three-dimensional electronic structure of the nematic and antiferromagnetic phases of NaFeAs from detwinned ARPES measurements

abstract: We report a comprehensive ARPES study of NaFeAs, a prototypical parent compound of the Fe-based superconductors. By mechanically detwinning the samples, we show that in the nematic phase (below the structural transition at $T_s$ = 54 K but above the antiferromagnetic transition at $T_N$ = 43 K) spectral weight is detected on only the elliptical electron pocket along the longer $a_{orth}$ axis. This dramatic anisotropy is likely to arise as a result of coupling to a fluctuating antiferromagnetic order in the nematic phase. In the long-range ordered antiferromagnetic state below $T_N$, this single electron pocket is backfolded and hybridises with the hole bands, leading to the reconstructed Fermi surface. By careful analysis of the $k_z$ variation, we show that the backfolding of spectral weight in the magnetic phase has a wavector of ($\pi$,0,$\pi$), with the $c$-axis component being in agreement with the magnetic ordering in NaFeAs observed by neutron scattering. Our results clarify the origin of the tiny Fermi surfaces of NaFeAs at low temperatures and highlight the importance of the three-dimensional aspects of the electronic and magnetic properties of Fe-based superconductors.


title: Magnetar-like X-ray bursts suppress pulsar radio emission

abstract: Rotation-powered pulsars and magnetars are two different observational manifestations of neutron stars: rotation powered pulsars are rapidly spinning objects that are mostly observed as pulsating radio sources, while magnetars, neutron stars with the highest known magnetic fields, often emit short-duration X-ray bursts. Here we report simultaneous observations of the high-magnetic-field radio pulsar PSR J1119-6127 at X-ray, with XMM-Newton & NuSTAR, and at radio energies with Parkes radio telescope, during a period of magnetar-like bursts. The rotationally powered radio emission shuts off coincident with the occurrence of multiple X-ray bursts, and recovers on a time scale of ~70 seconds. These observations of related radio and X-ray phenomena further solidify the connection between radio pulsars and magnetars, and suggest that the pair plasma produced in bursts can disrupt the acceleration mechanism of radio emitting particles.


title: Spectral and spatial shaping of Smith Purcell Radiation

abstract: The Smith Purcell effect, observed when an electron beam passes in the vicinity of a periodic structure, is a promising platform for the generation of electromagnetic radiation in previously-unreachable spectral ranges. However, most of the studies of this radiation were performed on simple periodic gratings, whose radiation spectrum exhibits a single peak and its higher harmonics predicted by a well-established dispersion relation. Here, we propose a method to shape the spatial and spectral far-field distribution of the radiation using complex periodic and aperiodic gratings. We show, theoretically and experimentally, that engineering multiple peak spectra with controlled widths located at desired wavelengths is achievable using Smith-Purcell radiation. Our method opens the way to free-electron driven sources with tailored angular and spectral response, and gives rise to focusing functionality for spectral ranges where lenses are unavailable or inefficient.


title: Practical Integer Overflow Prevention

abstract: Integer overflows in commodity software are a main source for software bugs, which can result in exploitable memory corruption vulnerabilities and may eventually contribute to powerful software based exploits, i.e., code reuse attacks (CRAs). In this paper, we present IntGuard , a tool that can repair integer overflows with high-quality source code repairs. Specifically, given the source code of a program, IntGuard first discovers the location of an integer overflow error by using static source code analysis and satisfiability modulo theories (SMT) solving. IntGuard then generates integer multi-precision code repairs based on modular manipulation of SMT constraints as well as an extensible set of customizable code repair patterns. We have implemented and evaluated IntGuard with 2052 C programs (approx. 1 Mil. LOC) available in the currently largest open- source test suite for C/C++ programs and with a benchmark containing large and complex programs. The evaluation results show that IntGuard can precisely (i.e., no false positives are accidentally repaired), with low computational and runtime overhead repair programs with very small binary and source code blow-up. In a controlled experiment, we show that IntGuard is more time-effective and achieves a higher repair success rate than manually generated code repairs.


title: Radiative gravitational collapse to spherical, toroidal and higher genus black holes

abstract: We derive the matching conditions between FLRW and generalised Vaidya spacetimes with spherical, planar or hyperbolic symmetry, across timelike hypersurfaces. We then construct new models of gravitational collapse of FLRW spacetimes with a negative cosmological constant having electromagnetic radiation in the exterior. The final state of the collapse are asymptotically AdS black holes with spherical, toroidal or higher genus topologies. We analyse the collapse dynamics including trapped surface formation, for various examples.


title: First-Principles Supercell Calculations of Small Polarons with Proper Account for Long-Range Polarization Effects

abstract: We present a density functional theory (DFT) based supercell approach for modeling small polarons with proper account for the long-range elastic response of the material. Our analysis of the supercell dependence of the polaron properties (e.g., atomic structure, binding energy, and the polaron level) reveals long-range electrostatic effects and the electron-phonon interaction as the two main contributors. We develop a correction scheme for DFT polaron calculations that significantly reduces the dependence of polaron properties on the DFT exchange-correlation functional and the size of the supercell in the limit of strong electron-phonon coupling. Using our correction approach, we present accurate all-electron full-potential DFT results for small polarons in rocksalt MgO and rutile TiO$_2$.


title: A Convex Optimization Approach for Backstepping PDE Design: Volterra and Fredholm Operators

abstract: Backstepping design for boundary linear PDE is formulated as a convex optimization problem. Some classes of parabolic PDEs and a first-order hyperbolic PDE are studied, with particular attention to non-strict feedback structures. Based on the compactness of the Volterra and Fredholm type operators involved, their Kernels are approximated via polynomial functions. The resulting Kernel-PDEs are optimized using Sum-of-Squares(SOS) decomposition and solved via semidefinite programming, with sufficient precision to guarantee the stability of the system in the L2-norm. The effectiveness and limitations of the approach proposed are illustrated by numerical solutions of some Kernel-PDEs.


title: Life of \Pi

abstract: The heavy-quark contribution to the polarisation function $\Pi$ at higher perturbative orders is presently only known approximately. We scrutinise the accuracy of state-of-the-art approximations at three- and four-loop order. At three loops, we present for the first time a result with arbitrary numerical precision for general kinematics and compare to the best Pad\'e estimate. At four loops, we calculate the fourth (inverse) moment of the non-singlet heavy-quark vacuum polarisation in order to test the prediction for this moment based on Pad\'e approximation.


title: Astronomical bounds on a cosmological model allowing a general interaction in the dark sector

abstract: Non-gravitational interaction between two barotropic dark fluids, namely the pressureless dust and the dark energy in a spatially flat Friedmann-Lema\^{i}tre-Robertson-Walker model has been discussed. It is shown that for the interactions which are linear in terms the energy densities of the dark components and their first order derivatives, the net energy density is governed by a second order differential equation with constant coefficients. Taking a generalized interaction, which includes a number of already known interactions as special cases, the dynamics of the universe is described for three types of the dark energy equation of state, namely that of interacting quintessence, interacting vacuum energy density and interacting phantom. The models have been constrained using the standard cosmological probes, Supernovae type Ia data from joint light curve analysis and the observational Hubble parameter data. Two geometric tests, the cosmographic studies and the $Om$ diagnostic have been invoked so as to ascertain the behaviour of the present model vis-a-vis the $\Lambda$-cold dark matter model. We further discussed the interacting scenarios taking into account the thermodynamic considerations.


title: Route Optimization of Electric Vehicles based on Dynamic Wireless Charging

abstract: One of the barriers to adoption of Electric Vehicles (EVs) is the anxiety around the limited driving range. Recent proposals have explored charging EVs on the move, using dynamic wireless charging which enables power exchange between the vehicle and the grid while the vehicle is moving. In this article, we focus on the intelligent routing of EVs in need of charging so that they can make most efficient use of the so-called {\it Mobile Energy Disseminators} (MEDs) which operates as mobile charging stations. We present a method for routing EVs around MEDs on the road network, which is based on constraint logic programming and optimisation using a graph-based shortest path algorithm. The proposed method exploits Inter-Vehicle (IVC) communications in order to eco-route electric vehicles. We argue that combining modern communications between vehicles and state of the art technologies on energy transfer, the driving range of EVs can be extended without the need for larger batteries or overtly costly infrastructure. We present extensive simulations in city conditions that show the driving range and consequently the overall travel time of electric vehicles is improved with intelligent routing in the presence of MEDs.


title: Building topological quantum circuits: Majorana nanowire junctions

abstract: Topological quantum computation using non-Abelian Majorana zero modes localized in proximitized semiconductor nanowires requires careful electrostatic control of wire-junctions so as to manipulate and braid the zero modes enabling anyonic fault-tolerant gate operations. We theoretically investigate the topological superconducting properties of such elementary wire-junctions, the so-called T junctions, finding that the existence of the junction may nonperturbatively affect the Majorana behavior by introducing spurious non-topological subgap states mimicking zero-modes. We propose a possible solution to this potentially serious problem by showing that junctions made lithographically from two-dimensional (2D) electron gas systems may manifest robust subgap topological properties without any spurious zero modes. We propose a 2D structure that enables multiprobe tunneling experiments providing position-dependent spectroscopy, which can decisively settle outstanding open questions related to the origin of the zero-bias conductance peaks observed experimentally. We also find that junctions with trivial superconductors may result in local perturbations that induce extrinsic low-energy states similar to those associated with wire junctions.


title: Stable manifolds of two-dimensional biholomorphisms asymptotic to formal curves

abstract: Let $F\in\mathrm{Diff}(\mathbb{C}^2,0)$ be a germ of a holomorphic diffeomorphism and let $\Gamma$ be an invariant formal curve of $F$. Assume that the restricted diffeomorphism $F|_{\Gamma}$ is either hyperbolic attracting or rationally neutral non-periodic (these are the conditions that the diffeomorphism $F|_{\Gamma}$ should satisfy, if $\Gamma$ were convergent, in order to have orbits converging to the origin). Then we prove that $F$ has finitely many stable manifolds, either open domains or parabolic curves, consisting of and containing all converging orbits asymptotic to $\Gamma$. Our results generalize to the case where $\Gamma$ is a formal periodic curve of $F$.


title: Evidence of a topological Hall effect in Eu$_{1-x}$Sm$_x$TiO$_3$

abstract: We report on the observation of a possible topological Hall effect in thin films of the itinerant ferromagnet Eu1-xSmxTiO3. EuTiO3 and Eu0.955Sm0.045TiO3 films were grown by molecular beam epitaxy. The EuTiO3 film is insulating. The Hall resistivity of the Eu0.955Sm0.045TiO3 films exhibits the anomalous Hall effect below the Curie temperature of ~ 5 K as well as additional features that appear at 2 K. It is shown that these features are magnetic in origin and consistent with the topological Hall effect seen in materials systems with topologically nontrivial spin textures, such as skyrmions. The results open up interesting possibilities for epitaxial hybrid heterostructures that combine topological magnetic states, tunable carrier densities, and other phenomena.


title: The Horizontal Branch population of NGC 1851 as revealed by the Ultra-violet Imaging Telescope (UVIT)

abstract: We present the UV photometry of the globular cluster NGC 1851 using images acquired with the Ultra-violet Imaging Telescope (UVIT) onboard the ASTROSAT satellite. PSF-fitting photometric data derived from images in two far-UV (FUV) filters and one near-UV (NUV) filter are used to construct color-magnitude diagrams (CMD), in combination with HST and ground-based optical photometry. In the FUV, we detect only the bluest part of the cluster horizontal branch (HB); in the NUV, we detect the full extent of the HB, including the red HB, blue HB and a small number of RR Lyrae stars. UV variability was detected in 18 RR Lyrae stars, and 3 new variables were also detected in the central region. The UV/optical CMDs are then compared with isochrones of different age and metallicity (generated using Padova and BaSTI models) and synthetic HB (using helium enhanced $Y^2$ models). We are able to identify two populations among the HB stars, which are found to have either an age range of 10-12~Gyr, or a range in Y$_{ini}$ of 0.23 - 0.28, for a metallicity of [Fe/H] =$-$1.2 to $-$1.3. These estimations from the UV CMDs are consistent with those from optical studies. The almost complete sample of the HB stars tend to show a marginal difference in spatial/azimuthal distribution among the blue and red HB stars. This study thus show cases the capability of UVIT, with its excellent resolution and large field of view, to study the hot stellar population in Galactic globular clusters.


title: On weighted inductive limits of spaces of ultradifferentiable functions and their duals

abstract: In the first part of this paper we discuss the completeness of two general classes of weighted inductive limits of spaces of ultradifferentiable functions. In the second part we study their duals and characterize these spaces in terms of the growth of convolution averages of their elements. This characterization gives a canonical way to define a locally convex topology on these spaces and we give necessary and sufficient conditions for them to be ultrabornological. In particular, our results apply to spaces of convolutors for Gelfand-Shilov spaces.


title: RLT2-based Parallel Algorithms for Solving Large Quadratic Assignment Problems on Graphics Processing Unit Clusters

abstract: This paper discusses efficient parallel algorithms for obtaining strong lower bounds and exact solutions for large instances of the Quadratic Assignment Problem (QAP). Our parallel architecture is comprised of both multi-core processors and Compute Unified Device Architecture (CUDA) enabled NVIDIA Graphics Processing Units (GPUs) on the Blue Waters Supercomputing Facility at the University of Illinois at Urbana-Champaign. We propose novel parallelization of the Lagrangian Dual Ascent algorithm on the GPUs, which is used for solving a QAP formulation based on Level-2 Refactorization Linearization Technique (RLT2). The Linear Assignment sub-problems (LAPs) in this procedure are solved using our accelerated Hungarian algorithm [Date, Ketan, Rakesh Nagi. 2016. GPU-accelerated Hungarian algorithms for the Linear Assignment Problem. Parallel Computing 57 52-72]. We embed this accelerated dual ascent algorithm in a parallel branch-and-bound scheme and conduct extensive computational experiments on single and multiple GPUs, using problem instances with up to 42 facilities from the QAPLIB. The experiments suggest that our GPU-based approach is scalable and it can be used to obtain tight lower bounds on large QAP instances. Our accelerated branch-and-bound scheme is able to comfortably solve Nugent and Taillard instances (up to 30 facilities) from the QAPLIB, using modest number of GPUs.


title: The Tensor Networks Anthology: Simulation techniques for many-body quantum lattice systems

abstract: We present a compendium of numerical simulation techniques, based on tensor network methods, aiming to address problems of many-body quantum mechanics on a classical computer. The core setting of this anthology are lattice problems in low spatial dimension at finite size, a physical scenario where tensor network methods, both Density Matrix Renormalization Group and beyond, have long proven to be winning strategies. Here we explore in detail the numerical frameworks and methods employed to deal with low-dimension physical setups, from a computational physics perspective. We focus on symmetries and closed-system simulations in arbitrary boundary conditions, while discussing the numerical data structures and linear algebra manipulation routines involved, which form the core libraries of any tensor network code. At a higher level, we put the spotlight on loop-free network geometries, discussing their advantages, and presenting in detail algorithms to simulate low-energy equilibrium states. Accompanied by discussions of data structures, numerical techniques and performance, this anthology serves as a programmer's companion, as well as a self-contained introduction and review of the basic and selected advanced concepts in tensor networks, including examples of their applications.


title: Market impact with multi-timescale liquidity

abstract: We present an extended version of the recently proposed "LLOB" model for the dynamics of latent liquidity in financial markets. By allowing for finite cancellation and deposition rates within a continuous reaction-diffusion setup, we account for finite memory effects on the dynamics of the latent order book. We compute in particular the finite memory corrections to the square root impact law, as well as the impact decay and the permanent impact of a meta-order. The latter is found to be linear in the traded volume and independent of the trading rate, as dictated by no-arbitrage arguments. In addition, we consider the case of a spectrum of cancellation and deposition rates, which allows us to obtain a square root impact law for moderate participation rates, as observed empirically. Our multi-scale framework also provides an alternative solution to the so-called price diffusivity puzzle in the presence of a long-range correlated order flow.


title: Saturation of Berge Hypergraphs

abstract: Given a graph $F$, a hypergraph is a Berge-$F$ if it can be obtained by expanding each edge in $F$ to a hyperedge containing it. A hypergraph $H$ is Berge-$F$-saturated if $H$ does not contain a subgraph that is a Berge-$F$, but for any edge $e\in E(\overline{H})$, $H+e$ does. The $k$-uniform saturation number of Berge-$F$ is the minimum number of edges in a $k$-uniform Berge-$F$-saturated hypergraph on $n$ vertices. For $k=2$ this definition coincides with the classical definition of saturation for graphs. In this paper we study the saturation numbers for Berge triangles, paths, cycles, stars and matchings in $k$-uniform hypergraphs.


title: Geodesic behavior for Finsler metrics of constant positive flag curvature on $S^2$

abstract: We study non-reversible Finsler metrics with constant flag curvature 1 on S^2 and show that the geodesic flow of every such metric is conjugate to that of one of Katok's examples, which form a 1-parameter family. In particular, the length of the shortest closed geodesic is a complete invariant of the geodesic flow. We also show, in any dimension, that the geodesic flow of a Finsler metrics with constant positive flag curvature is completely integrable. Finally, we give an example of a Finsler metric on~$S^2$ with positive flag curvature such that no two closed geodesics intersect and show that this is not possible when the metric is reversible or have constant flag curvature


title: On annihilators of bounded $(\frak g, \frak k)$-modules

abstract: Let $\frak g$ be a semisimple Lie algebra and $\frak k\subset\frak g$ be a reductive subalgebra. We say that a $\frak g$-module $M$ is a bounded $(\frak g, \frak k)$-module if $M$ is a direct sum of simple finite-dimensional $\frak k$-modules and the multiplicities of all simple $\frak k$-modules in that direct sum are universally bounded. The goal of this article is to show that the "boundedness" property for a simple $(\frak g, \frak k)$-module $M$ is equivalent to a property of the associated variety of the annihilator of $M$ (this is the closure of a nilpotent coadjoint orbit inside $\frak g^*$) under the assumption that the main field is algebraically closed and of characteristic 0. In particular this implies that if $M_1, M_2$ are simple $(\frak g, \frak k)$-modules such that $M_1$ is bounded and the associated varieties of the annihilators of $M_1$ and $M_2$ coincide then $M_2$ is also bounded. This statement is a geometric analogue of a purely algebraic fact due to I. Penkov and V. Serganova and it was posed as a conjecture in my Ph.D. thesis.


title: Diffusivities bounds in the presence of Weyl corrections

abstract: In this paper, we investigate the behavior of the thermoelectric DC conductivities in the presence of Weyl corrections with momentum dissipation in the incoherent limit. Moreover, we compute the butterfly velocity and study the charge and energy diffusion with broken translational symmetry. Our results show that the Weyl coupling $\gamma$, violates the bounds on the charge and energy diffusivity. It is also shown that the Weyl corrections violate the bound on the DC electrical conductivity in the incoherent limit.


title: Attacks on the Search-RLWE problem with small errors

abstract: The Ring Learning-With-Errors (RLWE) problem shows great promise for post-quantum cryptography and homomorphic encryption. We describe a new attack on the non-dual search RLWE problem with small error widths, using ring homomorphisms to finite fields and the chi-squared statistical test. In particular, we identify a "subfield vulnerability" (Section 5.2) and give a new attack which finds this vulnerability by mapping to a finite field extension and detecting non-uniformity with respect to the number of elements in the subfield. We use this attack to give examples of vulnerable RLWE instances in Galois number fields. We also extend the well-known search-to-decision reduction result to Galois fields with any unramified prime modulus q, regardless of the residue degree f of q, and we use this in our attacks. The time complexity of our attack is O(nq2f), where n is the degree of K and f is the residue degree of q in K. We also show an attack on the non-dual (resp. dual) RLWE problem with narrow error distributions in prime cyclotomic rings when the modulus is a ramified prime (resp. any integer). We demonstrate the attacks in practice by finding many vulnerable instances and successfully attacking them. We include the code for all attacks.


title: Mixed Precision Training

abstract: Deep neural networks have enabled progress in a wide variety of applications. Growing the size of the neural network typically results in improved accuracy. As model sizes grow, the memory and compute requirements for training these models also increases. We introduce a technique to train deep neural networks using half precision floating point numbers. In our technique, weights, activations and gradients are stored in IEEE half-precision format. Half-precision floating numbers have limited numerical range compared to single-precision numbers. We propose two techniques to handle this loss of information. Firstly, we recommend maintaining a single-precision copy of the weights that accumulates the gradients after each optimizer step. This single-precision copy is rounded to half-precision format during training. Secondly, we propose scaling the loss appropriately to handle the loss of information with half-precision gradients. We demonstrate that this approach works for a wide variety of models including convolution neural networks, recurrent neural networks and generative adversarial networks. This technique works for large scale models with more than 100 million parameters trained on large datasets. Using this approach, we can reduce the memory consumption of deep learning models by nearly 2x. In future processors, we can also expect a significant computation speedup using half-precision hardware units.


title: Symplectic Boundary Conditions and Cohomology

abstract: We introduce new boundary conditions for differential forms on symplectic manifolds with boundary. These boundary conditions, dependent on the symplectic structure, allows us to write down elliptic boundary value problems for both second-order and fourth-order symplectic Laplacians and establish Hodge theories for the cohomologies of primitive forms on manifolds with boundary. We further use these boundary conditions to define a relative version of the primitive cohomologies and to relate primitive cohomologies with Lefschetz maps on manifolds with boundary. As we show, these cohomologies of primitive forms can distinguish certain K\"ahler structures of K\"ahler manifolds with boundary.


title: Feasibility of efficient room-temperature solid-state sources of indistinguishable single photons using ultrasmall mode volume cavities

abstract: Highly efficient sources of indistinguishable single photons that can operate at room temperature would be very beneficial for many applications in quantum technology. We show that the implementation of such sources is a realistic goal using solid-state emitters and ultrasmall mode volume cavities. We derive and analyze an expression for photon indistinguishability that accounts for relevant detrimental effects, such as plasmon-induced quenching and pure-dephasing. We then provide the general cavity and emitter conditions required to achieve efficient indistinguishable photon emission, and also discuss constraints due to phonon sideband emission. Using these conditions, we propose that a nanodiamond negatively charged silicon-vacancy center combined with a plasmonic-Fabry-Perot hybrid cavity is an excellent candidate system.


title: Confidence through Attention

abstract: Attention distributions of the generated translations are a useful bi-product of attention-based recurrent neural network translation models and can be treated as soft alignments between the input and output tokens. In this work, we use attention distributions as a confidence metric for output translations. We present two strategies of using the attention distributions: filtering out bad translations from a large back-translated corpus, and selecting the best translation in a hybrid setup of two different translation systems. While manual evaluation indicated only a weak correlation between our confidence score and human judgments, the use-cases showed improvements of up to 2.22 BLEU points for filtering and 0.99 points for hybrid translation, tested on English<->German and English<->Latvian translation.


title: Probing Leptogenesis at Future Colliders

abstract: We investigate the question whether leptogenesis, as a mechanism for explaining the baryon asymmetry of the universe, can be tested at future colliders. Focusing on the minimal scenario of two right-handed neutrinos, we identify the allowed parameter space for successful leptogenesis in the heavy neutrino mass range between $5$ and $50$ GeV. Our calculation includes the lepton flavour violating contribution from heavy neutrino oscillations as well as the lepton number violating contribution from Higgs decays to the baryon asymmetry of the universe. We confront this parameter space region with the discovery potential for heavy neutrinos at future lepton colliders, which can be very sensitive in this mass range via displaced vertex searches. Beyond the discovery of heavy neutrinos, we study the precision at which the flavour-dependent active-sterile mixing angles can be measured. The measurement of these mixing angles at future colliders can test whether a minimal type I seesaw mechanism is the origin of the light neutrino masses, and it can be a first step towards probing leptogenesis as the mechanism of baryogenesis. We discuss how a stronger test could be achieved with an additional measurement of the heavy neutrino mass difference.


title: Erdos-Hajnal conjecture for graphs with bounded VC-dimension

abstract: The Vapnik-Chervonenkis dimension (in short, VC-dimension) of a graph is defined as the VC-dimension of the set system induced by the neighborhoods of its vertices. We show that every $n$-vertex graph with bounded VC-dimension contains a clique or an independent set of size at least $e^{(\log n)^{1 - o(1)}}$. The dependence on the VC-dimension is hidden in the $o(1)$ term. This improves the general lower bound, $e^{c\sqrt{\log n}}$, due to Erdos and Hajnal, which is valid in the class of graphs satisfying any fixed nontrivial hereditary property. Our result is almost optimal and nearly matches the celebrated Erdos-Hajnal conjecture, according to which one can always find a clique or an independent set of size at least $e^{\Omega(\log n)}$. Our results partially explain why most geometric intersection graphs arising in discrete and computational geometry have exceptionally favorable Ramsey-type properties. Our main tool is a partitioning result found by Lov\'asz-Szegedy and Alon-Fischer-Newman, which is called the "ultra-strong regularity lemma" for graphs with bounded VC-dimension. We extend this lemma to $k$-uniform hypergraphs, and prove that the number of parts in the partition can be taken to be $(1/\varepsilon)^{O(d)}$, improving the original bound of $(1/\varepsilon)^{O(d^2)}$ in the graph setting. We show that this bound is tight up to an absolute constant factor in the exponent. Moreover, we give an $O(n^k)$-time algorithm for finding a partition meeting the requirements. Finally, we establish tight bounds on Ramsey-Tur\'an numbers for graphs with bounded VC-dimension.


title: Bayesian Attitude Estimation with the Matrix Fisher Distribution on SO(3)

abstract: This paper focuses on a stochastic formulation of Bayesian attitude estimation on the special orthogonal group. In particular, an exponential probability density model for random matrices, referred to as the matrix Fisher distribution is used to represent the uncertainties of attitude estimates and measurements in a global fashion. Various stochastic properties of the matrix Fisher distribution are derived on the special orthogonal group, and based on these, two types of intrinsic frameworks for Bayesian attitude estimation are constructed. These avoid complexities or singularities of the attitude estimators developed in terms of quaternions. The proposed approaches are particularly useful to deal with large estimation errors or large uncertainties for complex maneuvers to obtain accurate estimates of the attitude.


title: Measuring the Small-Scale Matter Power Spectrum with High-Resolution CMB Lensing

abstract: We present a method to measure the small-scale matter power spectrum using high-resolution measurements of the gravitational lensing of the Cosmic Microwave Background (CMB). To determine whether small-scale structure today is suppressed on scales below 10 kiloparsecs (corresponding to M < 10^9 M_sun), one needs to probe CMB-lensing modes out to L ~ 35,000, requiring a CMB experiment with about 20 arcsecond resolution or better. We show that a CMB survey covering 4,000 square degrees of sky, with an instrumental sensitivity of 0.5 uK-arcmin at 18 arcsecond resolution, could distinguish between cold dark matter and an alternative, such as 1 keV warm dark matter or 10^(-22) eV fuzzy dark matter with about 4-sigma significance. A survey of the same resolution with 0.1 uK-arcmin noise could distinguish between cold dark matter and these alternatives at better than 20-sigma significance; such high-significance measurements may also allow one to distinguish between a suppression of power due to either baryonic effects or the particle nature of dark matter, since each impacts the shape of the lensing power spectrum differently. CMB temperature maps yield higher signal-to-noise than polarization maps in this small-scale regime; thus, systematic effects, such as from extragalactic astrophysical foregrounds, need to be carefully considered. However, these systematic concerns can likely be mitigated with known techniques. Next-generation CMB lensing may thus provide a robust and powerful method of measuring the small-scale matter power spectrum.


title: Emergent Complexity via Multi-Agent Competition

abstract: Reinforcement learning algorithms can train agents that solve problems in complex, interesting environments. Normally, the complexity of the trained agent is closely related to the complexity of the environment. This suggests that a highly capable agent requires a complex environment for training. In this paper, we point out that a competitive multi-agent environment trained with self-play can produce behaviors that are far more complex than the environment itself. We also point out that such environments come with a natural curriculum, because for any skill level, an environment full of agents of this level will have the right level of difficulty. This work introduces several competitive multi-agent environments where agents compete in a 3D world with simulated physics. The trained agents learn a wide variety of complex and interesting skills, even though the environment themselves are relatively simple. The skills include behaviors such as running, blocking, ducking, tackling, fooling opponents, kicking, and defending using both arms and legs. A highlight of the learned behaviors can be found here: https://goo.gl/eR7fbX


title: Nijenhuis operators on pre-Lie algebras

abstract: First we use a new approach to give a graded Lie algebra whose Maurer-Cartan elements characterize pre-Lie algebra structures. Then using this graded Lie bracket we define the notion of a Nijenhuis operator on a pre-Lie algebra which generates a trivial deformation of this pre-Lie algebra. There are close relationships between O-operators, Rota-Baxter operators and Nijenhuis operators on a pre-Lie algebra. In particular, a Nijenhuis operator "connects" two O-operators on a pre-Lie algebra whose any linear combination is still an O-operator in certain sense and hence compatible L-dendriform algebras appear naturally as the induced algebraic structures. For the case of the dual representation of the regular representation of a pre-Lie algebra, there is a geometric interpretation by introducing the notion of a pseudo-Hessian-Nijenhuis structure which gives rise to a sequence of pseudo-Hessian and pseudo-Hessian-Nijenhuis structures. Another application of Nijenhuis operators on pre-Lie algebras in geometry is illustrated by introducing the notion of a para-complex structure on a pre-Lie algebra and then studying paracomplex quadratic pre-Lie algebras and paracomplex pseudo-Hessian pre-Lie algebras in detail. Finally, we give some examples of Nijenhuis operators on pre-Lie algebras.


title: Influence of classical anisotropy fields on the properties of Heisenberg antiferromagnets within unified molecular field theory

abstract: A comprehensive study of the influence of classical anisotropy fields on the magnetic properties of Heisenberg antiferromagnets within unified molecular field theory versus temperature T, magnetic field H, and anisotropy field parameter hA1 is presented for systems comprised of identical crystallographically-equivalent local moments. The anisotropy field for collinear z-axis antiferromagnetic (AFM) ordering is constructed so that it is aligned in the direction of each ordered and/or field-induced thermal-average moment with a magnitude proportional to the moment, whereas that for XY anisotropy is defined to be in the direction of the projection of the moment onto the xy plane, again with a magnitude proportional to the moment. Properties studied include the zero-field Neel temperature TN, ordered moment, heat capacity and anisotropic magnetic susceptibility of the AFM phase versus T with moments aligned either along the z axis or in the xy plane. Also determined are the high-field magnetization perpendicular to the axis or plane of collinear or planar noncollinear AFM ordering, the high-field magnetization along the z axis of a collinear z-axis AFM, spin-flop (SF), and paramagnetic (PM) phases, and the free energies of these phases versus T, H, and hA1. Phase diagrams at T = 0 in the Hz-hA1 plane and at T > 0 in the Hz-T plane are constructed for spins S = 1/2. For hA1 = 0 the SF phase is stable at low field and the PM phase at high field with no AFM phase present. As hA1 increases, the phase diagram contains the AFM, SF and PM phases. Further increases in hA1 lead to the disappearance of the SF phase and the appearance of a tricritical point on the AFM-PM transition curve. Applications of the theory to extract hA1 from experimental low-field magnetic susceptibility data and high-field magnetization versus field isotherms for single crystals of AFMs are discussed.


title: Bosonic Pairings

abstract: We extend the inner product from bosonic Fock space to a pairing between suitable antifunctionals on the symmetric algebra. Our account is illustrated by the (Gaussian) half-form pairing between positive polarizations in the form needed for geometric quantization.


title: Heavy Higgs Searches: Flavour Matters

abstract: We point out that the stringent lower bounds on the masses of additional electrically neutral and charged Higgs bosons crucially depend on the flavour structure of their Yukawa interactions. We show that these bounds can easily be evaded by the introduction of flavour-changing neutral currents in the Higgs sector. As an illustration, we study the phenomenology of a two Higgs doublet model with a Yukawa texture singling out the third family of quarks and leptons. We combine constraints from low-energy flavour physics measurements, LHC measurements of the 125 GeV Higgs boson rates, and LHC searches for new heavy Higgs bosons. We propose novel LHC searches that could be performed in the coming years to unravel the existence of these new Higgs bosons.


title: Optimizing Long Short-Term Memory Recurrent Neural Networks Using Ant Colony Optimization to Predict Turbine Engine Vibration

abstract: This article expands on research that has been done to develop a recurrent neural network (RNN) capable of predicting aircraft engine vibrations using long short-term memory (LSTM) neurons. LSTM RNNs can provide a more generalizable and robust method for prediction over analytical calculations of engine vibration, as analytical calculations must be solved iteratively based on specific empirical engine parameters, making this approach ungeneralizable across multiple engines. In initial work, multiple LSTM RNN architectures were proposed, evaluated and compared. This research improves the performance of the most effective LSTM network design proposed in the previous work by using a promising neuroevolution method based on ant colony optimization (ACO) to develop and enhance the LSTM cell structure of the network. A parallelized version of the ACO neuroevolution algorithm has been developed and the evolved LSTM RNNs were compared to the previously used fixed topology. The evolved networks were trained on a large database of flight data records obtained from an airline containing flights that suffered from excessive vibration. Results were obtained using MPI (Message Passing Interface) on a high performance computing (HPC) cluster, evolving 1000 different LSTM cell structures using 168 cores over 4 days. The new evolved LSTM cells showed an improvement of 1.35%, reducing prediction error from 5.51% to 4.17% when predicting excessive engine vibrations 10 seconds in the future, while at the same time dramatically reducing the number of weights from 21,170 to 11,810.


title: Solution by convex minimizationof the Cauchy problem forhyperbolic systems of conservation laws with convex entropy

abstract: We show that, for first-order systems of conservation laws with a strictly convex entropy,in particular for the very simple so-called "inviscid" Burgers equation,it is possible to address the Cauchy problem by a suitable convex minimizationproblem, quite similar to some problems arising in optimal transport or variational mean-field game theory.In the general case, we show that smooth, shock-free, solutions can be recoveredon some sufficiently small interval of time. In the special situation of the Burgers equation, we furthershow that every "entropy solution" (in the sense of Kruzhkov)including shocks, can be recovered, for arbitrarily long time intervals.


title: Towards a Practical Pedestrian Distraction Detection Framework using Wearables

abstract: Pedestrian safety continues to be a significant concern in urban communities and pedestrian distraction is emerging as one of the main causes of grave and fatal accidents involving pedestrians. The advent of sophisticated mobile and wearable devices, equipped with high-precision on-board sensors capable of measuring fine-grained user movements and context, provides a tremendous opportunity for designing effective pedestrian safety systems and applications. Accurate and efficient recognition of pedestrian distractions in real-time given the memory, computation and communication limitations of these devices, however, remains the key technical challenge in the design of such systems. Earlier research efforts in pedestrian distraction detection using data available from mobile and wearable devices have primarily focused only on achieving high detection accuracy, resulting in designs that are either resource intensive and unsuitable for implementation on mainstream mobile devices, or computationally slow and not useful for real-time pedestrian safety applications, or require specialized hardware and less likely to be adopted by most users. In the quest for a pedestrian safety system that achieves a favorable balance between computational efficiency, detection accuracy, and energy consumption, this paper makes the following main contributions: (i) design of a novel complex activity recognition framework which employs motion data available from users' mobile and wearable devices and a lightweight frequency matching approach to accurately and efficiently recognize complex distraction related activities, and (ii) a comprehensive comparative evaluation of the proposed framework with well-known complex activity recognition techniques in the literature with the help of data collected from human subject pedestrians and prototype implementations on commercially-available mobile and wearable devices.


title: Charge excitation dynamics in bosonic fractional Chern insulators

abstract: The experimental realization of the Harper-Hofstadter model in ultra-cold atomic gases has placed fractional states of matter in these systems within reach---a fractional Chern insulator state (FCI) is expected to emerge for sufficiently strong interactions when half-filling the lowest band. The experimental setups naturally allow to probe the dynamics of this topological state, yet little is known about its out-of-equilibrium properties. We explore, using density matrix renormalization group (DMRG) simulations, the response of the FCI state to spatially localized perturbations. After confirming the static properties of the phase we show that the characteristic, gapless features are clearly visible in the edge dynamics. We find that a local edge perturbation in this model propagates chirally independent of the perturbation strength. This contrasts the behavior of single particle models with counter-propagating edge states, such as the non-interacting Harper-Hofstadter model, where the chirality is manifest only for weak perturbations. Additionally, our simulations show that there is inevitable density leakage from the first row of sites into the bulk, preventing a naive chiral Luttinger theory interpretation of the dynamics.


title: Monte Carlo Tensor Network Renormalization

abstract: Techniques for approximately contracting tensor networks are limited in how efficiently they can make use of parallel computing resources. In this work we demonstrate and characterize a Monte Carlo approach to the tensor network renormalization group method which can be used straightforwardly on modern computing architectures. We demonstrate the efficiency of the technique and show that Monte Carlo tensor network renormalization provides an attractive path to improving the accuracy of a wide class of challenging computations while also providing useful estimates of uncertainty and a statistical guarantee of unbiased results.


title: Dark Matter Freeze-out During Matter Domination

abstract: We highlight the general scenario of dark matter freeze-out whilst the energy density of the universe is dominated by a decoupled non-relativistic species. Decoupling during matter domination changes the freeze-out dynamics, since the Hubble rate is parametrically different for matter and radiation domination. Furthermore, for successful Big Bang Nucleosynthesis the state dominating the early universe energy density must decay, this dilutes (or repopulates) the dark matter. As a result, the masses and couplings required to reproduce the observed dark matter relic density can differ significantly from radiation dominated freeze-out.


title: Hydrogen-rich Core Collapse Supernovae

abstract: Hydrogen-rich core collapse supernovae, known as "Type II" supernovae, are the most common type of stellar explosion realized in nature. They are defined by the presence of prominent hydrogen lines in their spectra. Type II supernovae are observed only in star-forming galaxies, and several events have been directly linked to massive star progenitors. Five main subclasses are identified: Type IIP (displaying a plateau in their light curve), Type IIL (displaying a light curve decline), Type IIn (displaying narrow emission lines), Type IIb (displaying increasingly strong He features with time) and 87A-likes (displaying long-rising light curves similar to that of SN 1987A). Type IIP supernovae have been robustly established as the explosions of red supergiants, while the progenitors of Type IIL's remain elusive. Type IIn's are likely linked to luminous blue variables, Type IIb progenitors may be interacting binary systems and the prototype of the 87A-like class was observed to be the explosion of a blue supergiant. The diversity in progenitor mass, metallicity, binarity and rotation is likely responsible for the diversity in observed explosion types, but the connection between progenitor parameters and supernova properties is not yet entirely understood theoretically nor fully mapped observationally. New observational methods for constraining this connection are currently being implemented, including the analyses of large samples of events, making use of very early data (obtained hours to days from explosion) and statistical studies of host-galaxy properties.


title: FlexibleSUSY 2.0: Extensions to investigate the phenomenology of SUSY and non-SUSY models

abstract: We document major new features and improvements of FlexibleSUSY, a Mathematica and C++ package with a dependency on the external package SARAH, that generates fast and precise spectrum generators. The extensions presented here significantly increase the generality and capabilities of the FlexibleSUSY package, which already works with a wide class of models, while maintaining an elegant structure and easy to use interfaces. The FlexibleBSM extension makes it possible to also create spectrum generators for non-supersymmetric extensions of the Standard Model. The FlexibleCPV extension adds the option of complex parameters to the spectrum generators, allowing the study of many interesting models with new sources of $CP$ violation. FlexibleMW computes the decay of the muon for the generated model and thereby allows FlexibleSUSY to predict the mass of the $W$ boson from the input parameters by using the more precise electroweak input of $\{ G_F, M_Z, \alpha_{\text{em}} \}$ instead of $\{ M_W, M_Z, \alpha_{\text{em}} \}$. The FlexibleAMU extension provides a calculator of the anomalous magnetic moment of the muon in any model FlexibleSUSY can generate a spectrum for. FlexibleSAS introduces a new solver for the boundary value problem which makes use of semi-analytic expressions for dimensionful parameters to find solutions in models where the classic two-scale solver will not work such as the constrained E$_6$SSM. FlexibleEFTHiggs is a hybrid calculation of the Higgs mass which combines the virtues of both effective field theory calculations and fixed-order calculations. All of these extensions are included in FlexibleSUSY 2.0, which is released simultaneously with this manual.


title: The JCMT Transient Survey: Identifying Submillimetre Continuum Variability over Several Year Timescales Using Archival JCMT Gould Belt Survey Observations

abstract: Investigating variability at the earliest stages of low-mass star formation is fundamental in understanding how a protostar assembles mass. While many simulations of protostellar disks predict non-steady accretion onto protostars, deeper investigation requires robust observational constraints on the frequency and amplitude of variability events characterised across the observable SED. In this study, we develop methods to robustly analyse repeated observations of an area of the sky for submillimetre variability in order to determine constraints on the magnitude and frequency of deeply embedded protostars. We compare \mbox{850 $\mu$m} JCMT Transient Survey data with archival JCMT Gould Belt Survey data to investigate variability over 2-4 year timescales. Out of 175 bright, independent emission sources identified in the overlapping fields, we find 7 variable candidates, 5 of which we classify as \textit{Strong} and the remaining 2 as \textit{Extended} to indicate the latter are associated with larger-scale structure. For the \textit{Strong} variable candidates, we find an average fractional peak brightness change per year of |4.0|\% yr$^{-1}$ with a standard deviation of $2.7\%\mathrm{\:yr}^{-1}$. In total, 7\% of the protostars associated with \mbox{850 $\mu$m} emission in our sample show signs of variability. Four of the five \textit{Strong} sources are associated with a known protostar. The remaining source is a good follow-up target for an object that is anticipated to contain an enshrouded, deeply embedded protostar. In addition, we estimate the \mbox{850 $\mu$m} periodicity of the submillimetre variable source, EC 53, to be \mbox{567 $\pm$ 32 days} based on the archival Gould Belt Survey data.


title: Evidence of Absence of Tidal Features in the Outskirts of Ultra Diffuse Galaxies in the Coma Cluster

abstract: We study the presence of tidal features associated with ultra diffuse galaxies (UDGs) in galaxy clusters. Specifically, we stack deep Subaru images of UDGs in the Coma cluster to determine whether they show position angle twists at large radii. Selecting galaxies with central surface brightness (g-band) >24 magarcsec^{-2} and projected half-light radius >1.5 kpc, we identify 287 UDGs in the Yagi et al. 2016 catalog of low surface brightness Coma objects. The UDGs have apparent spheroidal shapes with median Sersic index ~ 0.8 and median axis-ratio ~ 0.7. The images are processed by masking all background objects and rotating to align the major axis before stacking them in bins of properties such as axis ratio, angle of major axis with respect to the cluster center and separation from cluster center. Our image stacks reach further than 7 kpc (~4r_e). Analysis of the isophotes of the stacks reveal that the ellipticity remains constant up to the last measured point, which means that the individual galaxies have a non-varying position angle and axis ratio and show no evidence for tidal disruption out to ~ 4r_e. We demonstrate this explicitly by comparing our stacks with stacks of model UDGs with and without tidal features in their outskirts. We infer that the average tidal radius of the Coma UDGs is >7 kpc and estimate that the average dark matter fraction within the tidal radius of the UDGs inhabiting the innermost 0.5 Mpc of Coma is > 99%.


title: Is the Milky Way still breathing? RAVE-Gaia streaming motions

abstract: We use data from the Radial Velocity Experiment (RAVE) and the Tycho-Gaia astrometric solution catalogue (TGAS) to compute the velocity fields yielded by the radial (VR), azimuthal (Vphi) and vertical (Vz) components of associated Galactocentric velocity. We search in particular for variation in all three velocity components with distance above and below the disc midplane, as well as how each component of Vz (line-of-sight and tangential velocity projections) modifies the obtained vertical structure. To study the dependence of velocity on proper motion and distance we use two main samples: a RAVE sample including proper motions from the Tycho-2, PPMXL and UCAC4 catalogues, and a RAVE-TGAS sample with inferred distances and proper motions from the TGAS and UCAC5 catalogues. In both samples, we identify asymmetries in VR and Vz. Below the plane we find the largest radial gradient to be dVR / dR = -7.01+- 0.61 km\s kpc, in agreement with recent studies. Above the plane we find a similar gradient with dVR / dR= -9.42+- 1.77 km\s kpc. By comparing our results with previous studies, we find that the structure in Vz is strongly dependent on the adopted proper motions. Using the Galaxia Milky Way model, we demonstrate that distance uncertainties can create artificial wave-like patterns. In contrast to previous suggestions of a breathing mode seen in RAVE data, our results support a combination of bending and breathing modes, likely generated by a combination of external or internal and external mechanisms.


title: A viable QCD axion in the MeV mass range

abstract: The QCD axion is one of the most compelling solutions of the strong CP problem. There are major current efforts into searching for an ultralight, invisible axion, which is believed to be the only phenomenologically viable realization of the QCD axion. Visible axions with decay constants at or below the electroweak scale are believed to have been long excluded by laboratory searches. Considering the significance of the axion solution of the strong CP problem, we revisit experimental constraints on QCD axions in the O(10 MeV) mass window. In particular, we find a variant axion model that remains compatible with existing constraints. This model predicts new states at the GeV scale coupled hadronically, and a variety of low-energy axion signatures, such as rare meson decays, nuclear de-excitations via axion emission, production in $e^+e^-$ annihilation and fixed target experiments. This reopens the possibility of solving the strong CP problem at the GeV scale.


title: Gas flows in the circumgalactic medium around simulated high-redshift galaxies

abstract: We analyse the properties of circumgalactic gas around simulated galaxies in the redshift range z >= 3, utilising a new sample of cosmological zoom simulations. These simulations are intended to be representative of the observed samples of Lyman-alpha emitters recently obtained with the MUSE instrument (halo masses ~10^10-10^11 solar masses). We show that supernova feedback has a significant impact on both the inflowing and outflowing circumgalactic medium by driving outflows, reducing diffuse inflow rates, and by increasing the neutral fraction of inflowing gas. By temporally stacking simulation outputs we find that significant net mass exchange occurs between inflowing and outflowing phases: none of the phases are mass-conserving. In particular, we find that the mass in neutral outflowing hydrogen declines exponentially with radius as gas flows outwards from the halo centre. This is likely caused by a combination of both fountain-like cycling processes and gradual photo/collisional ionization of outflowing gas. Our simulations do not predict the presence of fast-moving neutral outflows in the CGM. Neutral outflows instead move with modest radial velocities (~ 50 kms^-1), and the majority of the kinetic energy is associated with tangential rather than radial motion.


title: Status of a Deep Learning Based Measurement of the Inclusive Muon Neutrino Charged-current Cross Section in the NOvA Near Detector

abstract: NOvA is a long-baseline neutrino oscillation experiment. It uses the NuMI beam from Fermilab and two sampling calorimeter detectors placed off-axis from the beam. The 293 ton Near Detector measures the unoscillated neutrino energy spectrum, which can be used to predict the neutrino energy spectrum observed at the 14 kton Far Detector. The Near Detector also provides an excellent opportunity to measure neutrino interaction cross sections with high statistics, which will benefit current and future long-baseline neutrino oscillation experiments. This analysis implements new algorithms to identify $\nu_{\mu}$ charge-current events by using visual deep learning tools such as convolutional neural networks. We present the status of a measurement of the inclusive $\nu_{\mu}$ CC cross section in the NOvA Near Detector.


title: Electromagnetic transitions in the algebraic cluster model

abstract: We study electromagnetic transition rates in the framework of the algebraic cluster model. The concept of shape-phase transitions is used to propose a mechanism that allows to have interband and intraband quadrupole transitions of comparable strength, as observed in 12C.


title: Global Anomalies, Discrete Symmetries, and Hydrodynamic Effective Actions

abstract: We derive effective actions for parity-violating fluids in both $(3+1)$ and $(2+1)$ dimensions, including those with anomalies. As a corollary we confirm the most general constitutive relations for such systems derived previously using other methods. We discuss in detail connections between parity-odd transport and underlying discrete symmetries. In (3+1) dimensions we elucidate connections between anomalous transport coefficients and global anomalies, and clarify a previous puzzle concerning transports and local gravitational anomalies.


title: Quantum jumps on Anderson attractors

abstract: In a closed single-particle quantum system, spatial disorder induces Anderson localization of eigenstates and halts wave propagation. The phenomenon is vulnerable to interaction with environment and decoherence, that is believed to restore normal diffusion. We demonstrate that for a class of experimentally feasible non-Hermitian dissipators, which admit signatures of localization in asymptotic states, quantum particle opts between diffusive and ballistic regimes, depending on the phase parameter of dissipators, with sticking about localization centers. In diffusive regime, statistics of quantum jumps is non-Poissonian and has a power-law interval, a footprint of intermittent locking in Anderson modes. Ballistic propagation reflects dispersion of an ordered lattice and introduces a new timescale for jumps with non-monotonous probability distribution. Hermitian dephasing dissipation makes localization features vanish, and Poissonian jump statistics along with normal diffusion are recovered.


title: Examining the time dependence of DAMA's modulation amplitude

abstract: If dark matter is composed of weakly interacting particles, Earth's orbital motion may induce a small annual variation in the rate at which these particles interact in a terrestrial detector. The DAMA collaboration has identified at a 9.3$\sigma$ confidence level such an annual modulation in their event rate over two detector iterations, DAMA/NaI and DAMA/LIBRA, each with $\sim7$ years of observations. We statistically examine the time dependence of the modulation amplitudes, which "by eye" appear to be decreasing with time in certain energy ranges. We perform a chi-squared goodness of fit test of the average modulation amplitudes measured\ by the two detector iterations which rejects the hypothesis of a consistent modulation amplitude at greater than 80\%, 96\%, and 99.6\% for the 2--4~keVee, 2--5~keVee and 2--6~keVee energy ranges, respectively. We also find that among the 14 annual cycles there are three $\gtrsim 3\sigma$ departures from the average in the 5-6~keVee energy range. In addition, we examined several phenomenological models for the time dependence of the modulation amplitude. Using a maximum likelihood test, we find that descriptions of the modulation amplitude as decreasing with time are preferred over a constant modulation amplitude at anywhere between 1$\sigma$ and 3$\sigma$, depending on the phenomenological model for the time dependence and the signal energy range considered. A time dependent modulation amplitude is not expected for a dark matter signal, at least for dark matter halo morphologies consistent with the DAMA signal. New data from DAMA/LIBRA--phase2 will certainly aid in determining whether any apparent time dependence is a real effect or a statistical fluctuation.


title: Fourier-Mukai transforms of slope stable torsion-free sheaves on Weierstrass elliptic threefolds

abstract: We focus on a class of Weierstrass elliptic threefolds that allows the base of the fibration to be a Fano surface or a numerically $K$-trivial surface. We define the notion of limit tilt stability, and show that the Fourier-Mukai transform of a slope stable torsion-free sheaf satisfying a vanishing condition in codimension 2 (e.g. a reflexive sheaf) is a limit tilt stable object. We also show that the inverse Fourier-Mukai transform of a limit tilt semistable object of nonzero fiber degree is a slope semistable torsion-free sheaf, up to modification in codimension 2.


title: Event Reconstruction in the NOvA Experiment

abstract: The NOvA experiment observes oscillations in two channels (electron-neutrino appearance and muon-neutrino disappearance) using a predominantly muon-neutrino NuMI beam. The Near Detector records multiple overlapping neutrino interactions in each event and the Far Detector has a large background of cosmic rays due to being located on the surface. The oscillation analyses rely on the accurate reconstruction of neutrino interactions in order to precisely measure the neutrino energy and identify the neutrino flavor and interaction mode. Similarly, measurements of neutrino cross sections using the Near Detector require accurate identification of the particle content of each interaction. A series of pattern recognition techniques have been developed to split event records into individual spatially and temporally separated interactions, to estimate the interaction vertex, and to isolate and classify individual particles within the event. This combination of methods to achieve full event reconstruction in the NOvA detectors has discussed.


title: The MUSE Hubble Ultra Deep Field Survey: II. Spectroscopic redshifts and comparisons to color selections of high-redshift galaxies

abstract: We have conducted a two-layered spectroscopic survey (1'x1' ultra deep and 3'x3' deep regions) in the Hubble Ultra Deep Field (HUDF) with the Multi Unit Spectroscopic Explorer (MUSE). The combination of a large field of view, high sensitivity, and wide wavelength coverage provides an order of magnitude improvement in spectroscopically confirmed redshifts in the HUDF; i.e., 1206 secure spectroscopic redshifts for HST continuum selected objects, which corresponds to 15% of the total (7904). The redshift distribution extends well beyond z>3 and to HST/F775W magnitudes as faint as ~30 mag (AB, 1-sigma). In addition, 132 secure redshifts were obtained for sources with no HST counterparts that were discovered in the MUSE data cubes by a blind search for emission-line features. In total, we present 1338 high quality redshifts, which is a factor of eight increase compared with the previously known spectroscopic redshifts in the same field. We assessed redshifts mainly with the spectral features [OII] at z<1.5 (473 objects) and Lya at 2.9<z<6.7 (692 objects). With respect to F775W magnitude, a 50% completeness is reached at 26.5 mag for ultra deep and 25.5 mag for deep fields, and the completeness remains >~20% up to 28-29 mag and ~27 mag, respectively. We used the determined redshifts to test continuum color selection (dropout) diagrams of high-z galaxies. The selection condition for F336W dropouts successfully captures ~80% of the targeted z~2.7 galaxies. However, for higher redshift selections (F435W, F606W, and F775W dropouts), the success rates decrease to ~20-40%. We empirically redefine the selection boundaries to make an attempt to improve them to ~60%. The revised boundaries allow bluer colors that capture Lya emitters with high Lya equivalent widths falling in the broadbands used for the color-color selection. Along with this paper, we release the redshift and line flux catalog.


title: AI Buzzwords Explained: Multi-Agent Path Finding (MAPF)

abstract: Explanation of the hot topic "multi-agent path finding".


title: SPECULOOS exoplanet search and its prototype on TRAPPIST

abstract: One of the most significant goals of modern science is establishing whether life exists around other suns. The most direct path towards its achievement is the detection and atmospheric characterization of terrestrial exoplanets with potentially habitable surface conditions. The nearest ultracool dwarfs (UCDs), i.e. very-low-mass stars and brown dwarfs with effective temperatures lower than 2700 K, represent a unique opportunity to reach this goal within the next decade. The potential of the transit method for detecting potentially habitable Earth-sized planets around these objects is drastically increased compared to Earth-Sun analogs. Furthermore, only a terrestrial planet transiting a nearby UCD would be amenable for a thorough atmospheric characterization, including the search for possible biosignatures, with near-future facilities such as the James Webb Space Telescope. In this chapter, we first describe the physical properties of UCDs as well as the unique potential they offer for the detection of potentially habitable Earth-sized planets suitable for atmospheric characterization. Then, we present the SPECULOOS ground-based transit survey, that will search for Earth-sized planets transiting the nearest UCDs, as well as its prototype survey on the TRAPPIST telescopes. We conclude by discussing the prospects offered by the recent detection by this prototype survey of a system of seven temperate Earth-sized planets transiting a nearby UCD, TRAPPIST-1.


title: The sounds of the Little and Big Bangs

abstract: Studies of heavy ion collisions have discovered that tiny fireballs of new phase of matter -- quark gluon plasma (QGP) -- undergoes explosion, called the Little Bang. In spite of its small size, it is not only well described by hydrodynamics, but even small perturbations on top of the explosion turned to be well described by hydrodynamical sound modes. The cosmological Big Bang also went through phase transitions, the QCD and electroweak ones, which are expected to produce sounds as well. We discuss their subsequent evolution and hypothetical inverse acoustic cascade, amplifying the amplitude. Ultimately, collision of two sound waves leads to formation of gravity waves, with the smallest wavelength. We briefly discuss how those can be detected.


title: On the velocity of sound in water: theoretical aspects of Colladon's nineteenth century experiments

abstract: In 1827, Colladon carried out a series of experiments in Lac Leman (Lake Geneva, Switzerland) to measure the speed of sound in water. The purpose of our contribution is to treat this measurement as an inverse problem, and show, by theory how to solve the latter. It is thus revealed under what circumstances it is legitimate to employ the time-of-flight scheme underlying the Colladon experiments and how to bypass this scheme in order to fully account for the temporal and geometric characteristics of the source (of sound), the temporal characteristics of the received signal and the error incurred by the finite distance between the source and receiver.


title: Joint Weakly and Semi-Supervised Deep Learning for Localization and Classification of Masses in Breast Ultrasound Images

abstract: We propose a framework for localization and classification of masses in breast ultrasound (BUS) images. We have experimentally found that training convolutional neural network based mass detectors with large, weakly annotated datasets presents a non-trivial problem, while overfitting may occur with those trained with small, strongly annotated datasets. To overcome these problems, we use a weakly annotated dataset together with a smaller strongly annotated dataset in a hybrid manner. We propose a systematic weakly and semi-supervised training scenario with appropriate training loss selection. Experimental results show that the proposed method can successfully localize and classify masses with less annotation effort. The results trained with only 10 strongly annotated images along with weakly annotated images were comparable to results trained from 800 strongly annotated images, with the 95% confidence interval of difference -3.00%--5.00%, in terms of the correct localization (CorLoc) measure, which is the ratio of images with intersection over union with ground truth higher than 0.5. With the same number of strongly annotated images, additional weakly annotated images can be incorporated to give a 4.5% point increase in CorLoc, from 80.00% to 84.50% (with 95% confidence intervals 76.00%--83.75% and 81.00%--88.00%). The effects of different algorithmic details and varied amount of data are presented through ablative analysis.


title: On the Weyl anomaly of 4D Conformal Higher Spins: a holographic approach

abstract: We present a first attempt to derive the full (type-A and type-B) Weyl anomaly of four dimensional conformal higher spin (CHS) fields in a holographic way. We obtain the type-A and type-B Weyl anomaly coefficients for the whole family of 4D CHS fields from the one-loop effective action for massless higher spin (MHS) Fronsdal fields evaluated on a 5D bulk Poincar\'e-Einstein metric with an Einstein metric on its conformal boundary. To gain access to the type-B anomaly coefficient we assume, for practical reasons, a Lichnerowicz-type coupling of the bulk Fronsdal fields with the bulk background Weyl tensor. Remarkably enough, our holographic findings under this simplifying assumption are certainly not unknown: they match the results previously found on the boundary counterpart under the assumption of factorization of the CHS higher-derivative kinetic operator into Laplacians of "partially massless" higher spins on Einstein backgrounds.


title: Approximation of Bergman kernels by rational functions with fixed poles

abstract: We solved the problem of the best rational approximation of the Bergman kernels on the unit circle of the complex plane in the quadratic and uniform metrics.


title: Active Disassembly and Reassembly of Actin Networks Induces Distinct Biphasic Mechanics

abstract: Actin is a key component of the cytoskeleton, which plays central roles in cell motility, division, growth, and tensile strength. To enable this wide range of transient mechanical processes and properties, networks of actin filaments continuously disassemble and reassemble via active de/re-polymerization. However, the question remains as to how de/re-polymerization kinetics of individual actin filaments translate to time-varying mechanics of dis/re-assembling networks. To address this question, to ultimately elucidate the molecular mechanisms that enable cells to exhibit a myriad of transitory mechanical properties, we couple time-resolved active microrheology with microfluidics to measure the time-varying viscoelastic moduli of entangled and crosslinked actin networks during chemically-triggered network disassembly and reassembly. We also develop a corresponding mathematical model that relates the time-evolution of filament length to the resulting time-dependent storage moduli of evolving networks. During disassembly, we find that the moduli exhibit two distinct exponential decays, with experimental time constants of $\sim$169 and $\sim$47 min, which we show arises from a phase transition from a rigid percolated network to a non-rigid regime occurring after $\sim$90 min. During reassembly, measured moduli exhibit linear increase with time for $\sim$90 min, after which steady-state values are achieved. Our theoretical model shows that reassembly mechanics are dominated by elongation kinetics, and that elongation rates are much slower than has been recently assumed. Our measurements shed much needed light onto how polymerization kinetics map to time-varying mechanics of actively evolving cytoskeleton networks, and provide a powerful platform for studying other active and dynamic systems currently under intense investigation.


title: Approximation theorems for parabolic equations and movement of local hot spots

abstract: We prove a global approximation theorem for a general parabolic operator $L$, which asserts that if $v$ satisfies the equation $Lv=0$ in a spacetime region $\Omega \subset \mathbb{R}^{n+1}$ satisfying certain necessary topological condition, then it can be approximated in a H\"older norm by a global solution $u$ to the equation. If $\Omega$ is compact and $L$ is the usual heat operator, one can instead approximate the local solution $v$ by the unique solution that falls off at infinity to the Cauchy problem with a suitably chosen smooth, compactly supported initial datum. These results are next applied to prove the existence of global solutions to the equation $Lu=0$ with a local hot spot that moves along a prescribed curve for all time, up to a uniformly small error. Global solutions that exhibit isothermic hypersurfaces of prescribed topologies for all times and applications to the heat equation on the flat torus are discussed too.


title: Why Black Hole Information Loss is Paradoxical

abstract: I distinguish between two versions of the black hole information-loss paradox. The first arises from apparent failure of unitarity on the spacetime of a completely evaporating black hole, which appears to be non-globally-hyperbolic; this is the most commonly discussed version of the paradox in the foundational and semipopular literature, and the case for calling it `paradoxical' is less than compelling. But the second arises from a clash between a fully-statistical-mechanical interpretation of black hole evaporation and the quantum-field-theoretic description used in derivations of the Hawking effect. This version of the paradox arises long before a black hole completely evaporates, seems to be the version that has played a central role in quantum gravity, and is genuinely paradoxical. After explicating the paradox, I discuss the implications of more recent work on AdS/CFT duality and on the `Firewall paradox', and conclude that the paradox is if anything now sharper. The article is written at a (relatively) introductory level and does not assume advanced knowledge of quantum gravity.


title: Monochromatic CT Image Reconstruction from Current-Integrating Data via Deep Learning

abstract: In clinical CT, the x-ray source emits polychromatic x-rays, which are detected in the current-integrating mode. This physical process is accurately described by an energy-dependent non-linear integral model on the basis of the Beer-Lambert law. However, the non-linear model is too complicated to be directly solved for the image reconstruction, and is often approximated to a linear integral model in the form of the Radon transform, basically ignoring energy-dependent information. This model approximation would generate inaccurate quantification of attenuation image and significant beam-hardening artifacts. In this paper, we develop a deep-learning-based CT image reconstruction method to address the mismatch of computing model to physical model. Our method learns a nonlinear transformation from big data to correct measured projection data to accurately match the linear integral model, realize monochromatic imaging and overcome beam hardening effectively. The deep-learning network is trained and tested using clinical dual-energy dataset to demonstrate the feasibility of the proposed methodology. Results show that the proposed method can achieve a high accuracy of the projection correction with a relative error of less than 0.2%.


title: Monomial ideals of weighted oriented graphs

abstract: Let I=I(D) be the edge ideal of a weighted oriented graph D. We determine the irredundant irreducible decomposition of I. Also, we characterize the associated primes and the unmixed property of I. Furthermore, we give a combinatorial characterization for the unmixed property of I, when D is bipartite, D is a whisker or D is a cycle. Finally, we study the Cohen-Macaulay property of I.


title: momentuHMM: R package for generalized hidden Markov models of animal movement

abstract: Discrete-time hidden Markov models (HMMs) have become an immensely popular tool for inferring latent animal behaviors from telemetry data. Here we introduce an open-source R package, momentuHMM, that addresses many of the deficiencies in existing HMM software. Features include: 1) data pre-processing and visualization; 2) user-specified probability distributions for an unlimited number of data streams and latent behavior states; 3) biased and correlated random walk movement models, including "activity centers" associated with attractive or repulsive forces; 4) user-specified design matrices and constraints for covariate modelling of parameters using formulas familiar to most R users; 5) multiple imputation methods that account for measurement error and temporally-irregular or missing data; 6) seamless integration of spatio-temporal covariate raster data; 7) cosinor and spline models for cyclical and other complicated patterns; 8) model checking and selection; and 9) simulation. momentuHMM considerably extends the capabilities of existing HMM software while accounting for common challenges associated with telemetery data. It therefore facilitates more realistic hypothesis-driven animal movement analyses that have hitherto been largely inaccessible to non-statisticians. While motivated by telemetry data, the package can be used for analyzing any type of data that is amenable to HMMs. Practitioners interested in additional features are encouraged to contact the authors.


title: Indoor Massive MIMO Deployments for Uniformly High Wireless Capacity

abstract: Providing consistently high wireless capacity is becoming increasingly important to support the applications required by future digital enterprises. In this paper, we propose Eigen-direction-aware ZF (EDA-ZF) with partial coordination among base stations (BSs) and distributed interference suppression as a practical approach to achieve this objective. We compare our solution with Zero Forcing (ZF), entailing neither BS coordination or inter-cell interference mitigation, and Network MIMO (NeMIMO), where full BS coordination enables centralized inter-cell interference management. We also evaluate the performance of said schemes for three sub-6 GHz deployments with varying BS densities -- sparse, intermediate, and dense -- all with fixed total number of antennas and radiated power. Extensive simulations show that: (i) indoor massive MIMO implementing the proposed EDA-ZF provides uniformly good rates for all users; (ii) indoor network densification is detrimental unless full coordination is implemented; (iii) deploying NeMIMO pays off under strong outdoor interference, especially for cell-edge users.


title: Link Quality Aware Channel Allocation for Multichannel Body Sensor Networks

abstract: Body Sensor Network (BSN) is a typical Internet-of-Things (IoT) application for personalized health care. It consists of economically powered, wireless and implanted medical monitoring sensor nodes, which are designed to continually collect the medical information of the target patients. Multichannel is often used in BSNs to reduce the spectrum competition of the tremendous sensor nodes and the problem of channel assignment has attracted much research attention. The health sensing data in BSNs is often required to be delivered to a sink node (or server) before a certain deadline for real time monitoring or health emergency alarm. Therefore, deadline is of significant importance for multichannel allocation and scheduling. The existing works, though designed to meet the deadline, often overlook the impact of the unreliable wireless links. As a result, the health sensing data can still be overdue because of the scheduled lossy links. Besides, potential collisions in the schedules also incur considerable delay in delivering the sensing data. In this paper, we propose a novel deadline- driven Link quality Aware Channel Assignment scheme (LACA), where link quality, deadlines and collisions are jointly considered. LACA prioritizes links with urgent deadlines and heavy collisions. Besides, LACA allows the exploition of the spare slots for retransmissions on lossy links, which can further reduce the retransmission delay. Extensive simulation experiments show that compared to the existing approaches, LACA can better utilize the wireless spectrum and achieve higher packet delivery ratio before the deadline.


title: The Case for a Single System Image for Personal Devices

abstract: Computing technology has gotten cheaper and more powerful, allowing users to have a growing number of personal computing devices at their disposal. While this trend is beneficial for the user, it also creates a growing management burden for the user. Each device must be managed independently and users must repeat the same management tasks on the each device, such as updating software, changing configurations, backup, and replicating data for availability. To prevent the management burden from increasing with the number of devices, we propose that all devices run a single system image called a personal computing image. Personal computing images export a device-specific user interface on each device, but provide a consistent view of application and operating state across all devices. As a result, management tasks can be performed once on any device and will be automatically propagated to all other devices belonging to the user. We discuss evolutionary steps that can be taken to achieve personal computing images for devices and elaborate on challenges that we believe building such systems will face.


title: Performance of a crystalline silicon photovoltaic power plant during sandstorms

abstract: Solar photovoltaic power generation has achieved grid parity in a number of countries during the last decade. This revolution has begun in countries where the solar radiation intensity is not very high. Solar resource maps created by a number of agencies worldwide indicate that areas of high intensity of solar radiation suitable solar for power generation are in deserts and semi-arid zones. These regions are dusty and are subjected to both anthropogenic pollution and natural sand storms. The impact of these events on the solar photovoltaic (PV) plant yield is significant. An increase in diffuse component of solar radiation is accompanied by a naturally occurring sand storm. Often the solar panels are tilted at an optimal angle and the global tilted radiation value falling on the panels changes significantly during and after the sand storms. It is noted that the decrease in the power output of a solar photovoltaic plant during such events is high enough to warrant further consideration while forecasting the plant output. Inverter performance during such sand storms is found to be satisfactory. Integrating dust storm monitoring with PV power plant yield prediction algorithms is necessary for an improved operation of PV plants.