title: Existence of cube terms in finite algebras

abstract: We study the problem of whether a given finite algebra with finitely many basic operations contains a cube term; we give both structural and algorithmic results. We show that if such an algebra has a cube term then it has a cube term of dimension at most $N$, where the number $N$ depends on the arities of basic operations of the algebra and the size of the basic set. For finite idempotent algebras we give a tight bound on $N$ that, in the special case of algebras with more than $\binom{|A|}2$ basic operations, improves an earlier result of K. Kearnes and A. Szendrei. On the algorithmic side, we show that deciding the existence of cube terms is in P for idempotent algebras and in EXPTIME in general. Since an algebra contains a $k$-ary near unanimity operation if and only if it contains a $k$-dimensional cube term and generates a congruence distributive variety, our algorithm also lets us decide whether a given finite algebra has a near unanimity operation.


title: AIDA: Associative DNN Inference Accelerator

abstract: We propose AIDA, an inference engine for accelerating fully-connected (FC) layers of Deep Neural Network (DNN). AIDA is an associative in-memory processor, where the bulk of data never leaves the confines of the memory arrays, and processing is performed in-situ. AIDA area and energy efficiency strongly benefit from sparsity and lower arithmetic precision. We show that AIDA outperforms the state of art inference accelerator, EIE, by 14.5x (peak performance) and 2.5x (throughput).


title: An open-source sensor platform for analysis of group dynamics

abstract: The collaboration of several people in groups is becoming more and more important nowadays. Teamwork is often used for decision-making processes and for solving complex problems. Research in this area focuses on the quantification and analysis of behavior within and between groups. By using wearable electronic devices, such as badges, this quantification can be performed automatically and with high accuracy. The goal of this work is the design and implementation of a new firmware for the badges of the Rhythm project - an open-source project of the Human Dynamics Group of the MIT Media Lab. The firmware is characterized by a modular and extensible architecture and combines different techniques. These include a filesystem, which efficiently stores sequentially generated data based on a virtual memory abstraction, a serialization library, which enables a platform-independent exchange of structured data, and a time synchronization technique, which compensates frequency deviations of the oscillator. In addition, an automated test environment was designed for the verification of individual functional software components of the application. As a result of detailed analysis and special measures, the power consumption could be significantly reduced compared to the previous implementation. Due to the hierarchical and modular structure and the high degree of abstraction, the developed techniques can also be integrated into other projects and platforms.


title: PI-Edge: A Low-Power Edge Computing System for Real-Time Autonomous Driving Services

abstract: To simultaneously enable multiple autonomous driving services on affordable embedded systems, we designed and implemented {\pi}-Edge, a complete edge computing framework for autonomous robots and vehicles. The contributions of this paper are three-folds: first, we developed a runtime layer to fully utilize the heterogeneous computing resources of low-power edge computing systems; second, we developed an extremely lightweight operating system to manage multiple autonomous driving services and their communications; third, we developed an edge-cloud coordinator to dynamically offload tasks to the cloud to optimize client system energy consumption. To the best of our knowledge, this is the first complete edge computing system of a production autonomous vehicle. In addition, we successfully implemented {\pi}-Edge on a Nvidia Jetson and demonstrated that we could successfully support multiple autonomous driving services with only 11 W of power consumption, and hence proving the effectiveness of the proposed {\pi}-Edge system.


title: $B_s\to D_s^{(*)}l\nu$ Form Factors with Heavy HISQ Quarks

abstract: We present progress on an ongoing calculation of the $B_s\to D_s^{(*)} l \nu$ form factors calculated on the $n_f=2+1+1$ MILC ensembles and using the Highly Improved Staggered Quark action for all valence quarks. We perform the calculation at a range of $b$ quark masses (and lattice spacings) so that we can extrapolate to the physical $b$-quark mass.


title: Maximum and shape of interfaces in 3D Ising crystals

abstract: Dobrushin (1972) showed that the interface of a 3D Ising model with minus boundary conditions above the $xy$-plane and plus below is rigid (has $O(1)$-fluctuations) at every sufficiently low temperature. Since then, basic features of this interface -- such as the asymptotics of its maximum -- were only identified in more tractable random surface models that approximate the Ising interface at low temperatures, e.g., for the (2+1)D Solid-On-Solid model. Here we study the large deviations of the interface of the 3D Ising model in a cube of side-length $n$ with Dobrushin's boundary conditions, and in particular obtain a law of large numbers for $M_n$, its maximum: if the inverse-temperature $\beta$ is large enough, then $M_n / \log n \to 2/\alpha_\beta$ as $n\to\infty$, in probability, where $\alpha_\beta$ is given by a large deviation rate in infinite volume. We further show that, on the large deviation event that the interface connects the origin to height $h$, it consists of a 1D spine that behaves like a random walk, in that it decomposes into a linear (in $h$) number of asymptotically-stationary weakly-dependent increments that have exponential tails. As the number $T$ of increments diverges, properties of the interface such as its surface area, volume, and the location of its tip, all obey CLTs with variances linear in $T$. These results generalize to every dimension $d\geq 3$.


title: Tree-decorated planar maps

abstract: We introduce the set of (non-spanning) tree-decorated planar maps, and show that they are in bijection with the Cartesian product between the set of trees and the set of maps with a simple boundary. As a consequence, we count the number of tree decorated triangulations and quadrangulations with a given amount of faces and for a given size of the tree. Finally, we generalise the bijection to study other types of decorated planar maps and obtain explicit counting formulas for them.


title: Mechanism to Mitigate AVX-Induced Frequency Reduction

abstract: Modern Intel CPUs reduce their frequency when executing wide vector operations (AVX2 and AVX-512 instructions), as these instructions increase power consumption. The frequency is only increased again two milliseconds after the last code section containing such instructions has been executed in order to prevent excessive numbers of frequency changes. Due to this delay, intermittent use of wide vector operations can slow down the rest of the system significantly. For example, previous work has shown the performance of web servers to be reduced by up to 10% if the SSL library uses AVX-512 vector instructions. These performance variations are hard to predict during software development as the performance impact of vectorization depends on the specific workload. We describe a mechanism to reduce the slowdown caused by wide vector instructions without requiring extensive changes to existing software. Our design allows the developer to mark problematic AVX code regions. The scheduler then restricts execution of this code to a subset of the cores so that only these cores' frequency is affected. Threads are automatically migrated to a suitable core whenever necessary. We identify a suitable load balancing policy to ensure good utilization of all available cores. Our approach is able to reduce the performance variability caused by AVX2 and AVX-512 instructions by over 70%.


title: Adaptive virtual organisms: A compositional model for complex hardware-software binding

abstract: The relation between a structure and the function running on that structure is of central interest in many fields, including computer science, biology (organ vs. function), psychology (body vs. mind), architecture (designs vs. functionality), etc. Our paper addresses this question with reference to computer science recent hardware and software advances, particularly in areas as robotics, AI-hardware, self-adaptive systems, IoT, CPS, etc. At the modelling, conceptual level, our main contribution is the introduction of the concept of "virtual organism" (VO), to populate the intermediary level between rigid, slightly reconfigurable, hardware agents and abstract, intelligent, adaptive software agents. A virtual organism has a structure, resembling the hardware capabilities, and it runs low-level functions, implementing the software requirements. The model is compositional in space (allowing the virtual organisms to aggregate into larger organisms) and in time (allowing the virtual organisms to get composed functionalities). Technically, the virtual organisms studied here are in 2D and their structures are described by regular 2D pattens; adding the time dimension, we conclude the VO model is in 3D. By reconfiguration, an organism may change its structure to another structure belonging to the same 2D pattern. We illustrate the VO concept with three increasingly more complex VOs: (1) a tree collector; (2) a feeding cell; and (3) a collection of connected feeding cells. We implemented a simulator for tree collector organisms and the quantitative results show reconfigurable structures are better suited than fixed structures in dynamically changing environments. We briefly show how Agapia - a structured parallel, interactive programming language where dataflow and control flow structures can be freely mixed - may be used for getting quick implementations for VO's simulation.


title: Thermalization and chaos in QED$_{3}$

abstract: We study the real time dynamics of $N_F$ flavors of fermions coupled to a $U(1)$ gauge field in $2+1$ dimensions to leading order in a $1/N_F$ expansion. For large enough $N_{F}$, this is an interacting conformal field theory and describes the low energy properties of the Dirac spin liquid. We focus on thermalization and the onset of many-body quantum chaos which can be diagnosed from the growth of initally anti-commuting fermion field operators. We compute such anti-commutators in this gauge theory to leading order in $1/N_F$. We find that the anti-commutator grows exponentially in time and compute the quantum Lyapunov exponent. We briefly comment on chaos, locality, and gauge invariance.

categories: cond-mat.str-el cond-mat.stat-mech hep-th quant-ph


title: NNStreamer: Stream Processing Paradigm for Neural Networks, Toward Efficient Development and Execution of On-Device AI Applications

abstract: We propose nnstreamer, a software system that handles neural networks as filters of stream pipelines, applying the stream processing paradigm to neural network applications. A new trend with the wide-spread of deep neural network applications is on-device AI; i.e., processing neural networks directly on mobile devices or edge/IoT devices instead of cloud servers. Emerging privacy issues, data transmission costs, and operational costs signifies the need for on-device AI especially when a huge number of devices with real-time data processing are deployed. Nnstreamer efficiently handles neural networks with complex data stream pipelines on devices, improving the overall performance significantly with minimal efforts. Besides, nnstreamer simplifies the neural network pipeline implementations and allows reusing off-shelf multimedia stream filters directly; thus it reduces the developmental costs significantly. Nnstreamer is already being deployed with a product releasing soon and is open source software applicable to a wide range of hardware architectures and software platforms.


title: Systimator: A Design Space Exploration Methodology for Systolic Array based CNNs Acceleration on the FPGA-based Edge Nodes

abstract: The evolution of IoT based smart applications demand porting of artificial intelligence algorithms to the edge computing devices. CNNs form a large part of these AI algorithms. Systolic array based CNN acceleration is being widely advocated due its ability to allow scalable architectures. However, CNNs are inherently memory and compute intensive algorithms, and hence pose significant challenges to be implemented on the resource-constrained edge computing devices. Memory-constrained low-cost FPGA based devices form a substantial fraction of these edge computing devices. Thus, when porting to such edge-computing devices, the designer is left unguided as to how to select a suitable systolic array configuration that could fit in the available hardware resources. In this paper we propose Systimator, a design space exploration based methodology that provides a set of design points that can be mapped within the memory bounds of the target FPGA device. The methodology is based upon an analytical model that is formulated to estimate the required resources for systolic arrays, assuming multiple data reuse patterns. The methodology further provides the performance estimates for each of the candidate design points. We show that Systimator provides an in-depth analysis of resource-requirement of systolic array based CNNs. We provide our resource estimation results for porting of convolutional layers of TINY YOLO, a CNN based object detector, on a Xilinx ARTIX 7 FPGA.


title: Tango: A Deep Neural Network Benchmark Suite for Various Accelerators

abstract: Deep neural networks (DNNs) have been proving the effectiveness in various computing fields. To provide more efficient computing platforms for DNN applications, it is essential to have evaluation environments that include assorted benchmark workloads. Though a few DNN benchmark suites have been recently released, most of them require to install proprietary DNN libraries or resource-intensive DNN frameworks, which are hard to run on resource-limited mobile platforms or architecture simulators. To provide a more scalable evaluation environment, we propose a new DNN benchmark suite that can run on any platform that supports CUDA and OpenCL. The proposed benchmark suite includes the most widely used five convolution neural networks and two recurrent neural networks. We provide in-depth architectural statistics of these networks while running them on an architecture simulator, a server- and a mobile-GPU, and a mobile FPGA.


title: A Survey of FPGA Based Deep Learning Accelerators: Challenges and Opportunities

abstract: With the rapid development of in-depth learning, neural network and deep learning algorithms have been widely used in various fields, e.g., image, video and voice processing. However, the neural network model is getting larger and larger, which is expressed in the calculation of model parameters. Although a wealth of existing efforts on GPU platforms currently used by researchers for improving computing performance, dedicated hardware solutions are essential and emerging to provide advantages over pure software solutions. In this paper, we systematically investigate the neural network accelerator based on FPGA. Specifically, we respectively review the accelerators designed for specific problems, specific algorithms, algorithm features, and general templates. We also compared the design and implementation of the accelerator based on FPGA under different devices and network models and compared it with the versions of CPU and GPU. Finally, we present to discuss the advantages and disadvantages of accelerators on FPGA platforms and to further explore the opportunities for future research.


title: Application-Specific System Processor for the SHA-1 Hash Algorithm

abstract: This work proposes an Application-Specific System Processor (ASSP) hardware for the Secure Hash Algorithm 1 (SHA-1) algorithm. The proposed hardware was implemented in a Field Programmable Gate Array (FPGA) Xilinx Virtex 6 xc6vlx240t-1ff1156. The throughput and the occupied area were analyzed for several implementations in parallel instances of the hash algorithm. The results showed that the hardware proposed for the SHA-1 achieved a throughput of 0.644 Gbps for a single instance and slightly more than 28 Gbps for 48 instances in a single FPGA. Various applications such as password recovery, password validation, and high volume data integrity checking can be performed efficiently and quickly with an ASSP for SHA1.


title: Design and Analysis of Microfluidic Cell Counter using Spice Simulation

abstract: Microfluidic cytometers based on coulter principle have recently shown a great potential for point of care biosensors for medical diagnostics. In this study, the design of coulter based microfluidic cytometer is investigated by using electrical circuit simulations. We explore the effects of physical dimensions of the microelectrodes, the measurement volume, size/morphology of the targeted cells, electrical properties of the reagents in the measurement volume, and, the impedance of external readout circuit, on the sensitivity of the sensor. We show that the effect of microelectrode's surface area and the dielectric properties of the suspension medium should be carefully considered when characterizing the output response of the sensor. In particular, the area of microelectrodes can have significant effect on cells electrical opacity( the ratio of the cell impedance at high to low frequency) which is commonly used to distinguish between sub-population of the target cells( e.g. lymphocytes vs monocytes when counting white blood cells).Moreover, we highlight that the opacity response vs frequency can significantly vary depending upon whether the absolute cell impedance or the differential output impedance is used in the calculation. These insights can provide valuable guidelines for the design and characterization of coulter based microfluidic sensors.


title: The stochastic extinction and stability conditions for a class of malaria epidemic models

abstract: The stochastic extinction and stability in the mean of a family of SEIRS malaria models with a general nonlinear incidence rate is presented. The dynamics is driven by independent white noise processes from the disease transmission and natural death rates. The basic reproduction number $R^{*}_{0}$, the expected survival probability of the plasmodium $E(e^{-(\mu_{v}T_{1}+\mu T_{2})})$, and other threshold values are calculated. A sample Lyapunov exponential analysis for the system is utilized to obtain extinction results. Moreover, the rate of extinction of malaria is estimated, and innovative local Martingale and Lyapunov functional techniques are applied to establish the strong persistence, and asymptotic stability in the mean of the malaria-free steady population. %The extinction of malaria depends on $R^{*}_{0}$, and $E(e^{-(\mu_{v}T_{1}+\mu T_{2})})$. Moreover, for either $R^{*}_{0}<1$, or $E(e^{-(\mu_{v}T_{1}+\mu T_{2})})<\frac{1}{R^{*}_{0}}$, whenever $R^{*}_{0}\geq 1$, respectively, extinction of malaria occurs. Furthermore, the robustness of these threshold conditions to the intensity of noise from the disease transmission rate is exhibited. Numerical simulation results are presented.


title: CFOF: A Concentration Free Measure for Anomaly Detection

abstract: We present a novel notion of outlier, called the Concentration Free Outlier Factor, or CFOF. As a main contribution, we formalize the notion of concentration of outlier scores and theoretically prove that CFOF does not concentrate in the Euclidean space for any arbitrary large dimensionality. To the best of our knowledge, there are no other proposals of data analysis measures related to the Euclidean distance for which it has been provided theoretical evidence that they are immune to the concentration effect. We determine the closed form of the distribution of CFOF scores in arbitrarily large dimensionalities and show that the CFOF score of a point depends on its squared norm standard score and on the kurtosis of the data distribution, thus providing a clear and statistically founded characterization of this notion. Moreover, we leverage this closed form to provide evidence that the definition does not suffer of the hubness problem affecting other measures. We prove that the number of CFOF outliers coming from each cluster is proportional to cluster size and kurtosis, a property that we call semi-locality. We determine that semi-locality characterizes existing reverse nearest neighbor-based outlier definitions, thus clarifying the exact nature of their observed local behavior. We also formally prove that classical distance-based and density-based outliers concentrate both for bounded and unbounded sample sizes and for fixed and variable values of the neighborhood parameter. We introduce the fast-CFOF algorithm for detecting outliers in large high-dimensional dataset. The algorithm has linear cost, supports multi-resolution analysis, and is embarrassingly parallel. Experiments highlight that the technique is able to efficiently process huge datasets and to deal even with large values of the neighborhood parameter, to avoid concentration, and to obtain excellent accuracy.


title: Large-Scale Joint Topic, Sentiment & User Preference Analysis for Online Reviews

abstract: This paper presents a non-trivial reconstruction of a previous joint topic-sentiment-preference review model TSPRA with stick-breaking representation under the framework of variational inference (VI) and stochastic variational inference (SVI). TSPRA is a Gibbs Sampling based model that solves topics, word sentiments and user preferences altogether and has been shown to achieve good performance, but for large data set it can only learn from a relatively small sample. We develop the variational models vTSPRA and svTSPRA to improve the time use, and our new approach is capable of processing millions of reviews. We rebuild the generative process, improve the rating regression, solve and present the coordinate-ascent updates of variational parameters, and show the time complexity of each iteration is theoretically linear to the corpus size, and the experiments on Amazon data sets show it converges faster than TSPRA and attains better results given the same amount of time. In addition, we tune svTSPRA into an online algorithm ovTSPRA that can monitor oscillations of sentiment and preference overtime. Some interesting fluctuations are captured and possible explanations are provided. The results give strong visual evidence that user preference is better treated as an independent factor from sentiment.


title: Structures of cobordism-like modules induced from generic maps of codimension -2

abstract: The Reeb space of a smooth map whose codimension is minus is the space defined as the space of all connected components of inverse images. For generic maps such as Morse functions and their higher dimensional versions, they are polyhedra whose dimensions are equal to those of the target manifolds and which have simplicial structures compatible with (the canonical) simplicial structures of the source and the target manifolds, and in considerable cases they inherit fundamental and important invariants of source manifolds. In fact, Reeb spaces are fundamentall tools in the algebraic and differential topological theory of generic maps or the global singularity theory. As one of studies of global topological properties of Reeb spaces, Hiratuka and Saeki showed in 2013 that for generic maps or more precisely, maps compatible with simplicial structures of the manifolds, inducing simplicial structures on the Reeb spaces and having connected components of inverse images of regular values being not null-cobordant, the top-dimensional homology groups with appropriate coefficient rings of the Reeb spaces do not vanish. Later the author extended this theorem: the author has introduced cobordism-like groups based on adjacent relations of connected components of inverse images of regular values and shown a similar theorem.


title: RPS(1) Preferences

abstract: We consider a model for decision making based on an adaptive, k-period, learning process where the priors are selected according to Von Neumann-Morgenstern expected utility principle. A preference relation between two prospects is introduced, defined by the condition which prospect is selected more often. We show that the new preferences have similarities with the preferences obtained by Kahneman and Tversky (1979) in the context of the prospect theory. Additionally, we establish that in the limit of large learning period, the new preferences coincide with the expected utility principle.


title: Contact discontinuities for 3-D axisymmetric inviscid compressible flows in infinitely long cylinders

abstract: We prove the existence of a subsonic axisymmetric weak solution $({\bf u},\rho,p)$ with ${\bf u}=u_x{\bf e}_x+u_r{\bf e}_r+u_\theta{\bf e}_{\theta}$ to steady Euler system in a three-dimensional infinitely long cylinder $\mathcal{N}$ when prescribing the values of the entropy $(=\frac{p}{\rho^{\gamma}})$ and angular momentum density $(=ru_{\theta})$ at the entrance by piecewise $C^2$ functions with a discontinuity on a curve on the entrance of $\mathcal{N}$. Due to the variable entropy and angular momentum density (=swirl) conditions with a discontinuity at the entrance, the corresponding solution has a nonzero vorticity, nonzero swirl, and contains a contact discontinuity $r=g_D(x)$. We construct such a solution via Helmholtz decomposition. The key step is to decompose the Rankine-Hugoniot conditions on the contact discontinuity via Helmholtz decomposition so that the compactness of approximated solutions can be achieved. Then we apply the method of iteration to obtain a piecewise smooth subsonic flow with a contact discontinuity, nonzero vorticity, and nonzero angular momentum density. We also analyze the asymptotic behavior of the solution at far field.


title: MAD-GAN: Multivariate Anomaly Detection for Time Series Data with Generative Adversarial Networks

abstract: The prevalence of networked sensors and actuators in many real-world systems such as smart buildings, factories, power plants, and data centers generate substantial amounts of multivariate time series data for these systems. The rich sensor data can be continuously monitored for intrusion events through anomaly detection. However, conventional threshold-based anomaly detection methods are inadequate due to the dynamic complexities of these systems, while supervised machine learning methods are unable to exploit the large amounts of data due to the lack of labeled data. On the other hand, current unsupervised machine learning approaches have not fully exploited the spatial-temporal correlation and other dependencies amongst the multiple variables (sensors/actuators) in the system for detecting anomalies. In this work, we propose an unsupervised multivariate anomaly detection method based on Generative Adversarial Networks (GANs). Instead of treating each data stream independently, our proposed MAD-GAN framework considers the entire variable set concurrently to capture the latent interactions amongst the variables. We also fully exploit both the generator and discriminator produced by the GAN, using a novel anomaly score called DR-score to detect anomalies by discrimination and reconstruction. We have tested our proposed MAD-GAN using two recent datasets collected from real-world CPS: the Secure Water Treatment (SWaT) and the Water Distribution (WADI) datasets. Our experimental results showed that the proposed MAD-GAN is effective in reporting anomalies caused by various cyber-intrusions compared in these complex real-world systems.


title: Violation of Universal Lower Bound for the Shear Viscosity to Entropy Density Ratio in Dark Energy Dominated Accretion

abstract: The universal lower bound of the ratio of shear viscosity to entropy density is suggested by the string theory and gauge duality for any matter. We examined the ratio of shear viscosity to entropy density for viscous accretion flow towards a central gravitating object in the presence of dark energy. The ratio appears close to the universal lower bound for certain optically thin, hot accretion flows as they are embedded by strong magnetic field. Dark energy is a kind of exotic matter which has negative pressure. So dark energy creates repulsive force between the accreting particles, which indicates that shear viscosity of the flow becomes very low. Dark energy as accreting fluid has very high entropy density. The ratio should reach near to the lowest value for dark energy accretion. We wish to study what happens to the shear viscosity to entropy density ratio for viscous dark energy accretion flow.


title: On Classical Solutions of Rayleigh--Taylor Instability in Inhomogeneous Incompressible Viscous Fluids in Bounded Domains

abstract: We study the existence of unstable classical solutions of the Rayleigh--Taylor instability problem (abbr. RT problem) of an inhomogeneous incompressible viscous fluid in a bounded domain. We find that, by using an existence theory of (steady) Stokes problem and an iterative technique, the initial data of classical solutions of the linearized RT problem can be modified to new initial data, which can generate local-in-time classical solutions of the RT problem, and are close to the original initial data. Thus, we can use a classical bootstrap instability method to further obtain classical solutions of (nonlinear) RT instability based on the ones of linear RT instability.


title: CKM 2018 Summary of Working Group 4: Mixing and mixing-related CP violation in the B system $\Delta M$, $\Delta \Gamma$, $\phi_s$, $\phi_1/\beta$, $\phi_2/\alpha$, $\phi_3/\gamma$

abstract: We summarise the excellent talks of Working group 4 at the 10th CKM Workshop in Heidelberg taking place from 17th till 21st September 2018.


title: Discrete Spectra of Convolutions on Disks using Sturm-Liouville Theory

abstract: This paper presents a systematic study for analytic aspects of discrete spectra methods for convolution of functions supported on disks, according to the Sturm-Liouville theory. We then investigate different aspects of the presented theory in the cases of zero-value boundary condition and derivative boundary condition.


title: Light-weighted Saliency Detection with Distinctively Lower Memory Cost and Model Size

abstract: Deep neural networks (DNNs) based saliency detection approaches have succeed in recent years, and improved the performance by a great margin via increasingly sophisticated network architecture. Despite the performance improvement, the computational cost is excessively high for such low level visual task. In this work, we propose a light-weighted saliency detection approach with distinctively lower runtime memory cost and model size. We evaluated the performance of our approach on multiple benchmark datasets, and achieved competitive results comparing with state-of-the-art methods on multiple metrics. We also evaluated the computational cost of our approach with multiple measurements. The runtime memory cost of our approach is 42 to 99 times fewer comparing with the previous DNNs based methods. The model size of our approach is 63 to 129 times smaller, and takes less than 1 Megabytes storage space with out any deep compression technique.


title: Quantum Teleportation-Inspired Algorithm for Sampling Large Random Quantum Circuits

abstract: We show that low-depth random quantum circuits can be efficiently simulated by a quantum teleportation-inspired algorithm. By using logical qubits to redirect and teleport the quantum information in quantum circuits, the original circuits can be renormalized to new circuits with a smaller number of logical qubits. We demonstrate the algorithm to simulate several random quantum circuits, including 1D-chain 1000-qubit 42-depth, 2D-grid 125*8-qubit 42-depth and 2D-Bristlecone 72-qubit 32-depth circuits. Our results present a memory-efficient method with a clear physical picture to simulate low-depth random quantum circuits.


title: On the Ubiquity and Stellar Luminosity Dependence of Exocometary CO Gas: Detection around M Dwarf TWA 7

abstract: Millimeter observations of CO gas in planetesimal belts show a high detection rate around A stars, but few detections for later type stars. We present the first CO detection in a planetesimal belt around an M star, TWA 7. The optically thin CO (J=3-2) emission is co-located with previously identified dust emission from the belt, and the emission velocity structure is consistent with Keplerian rotation around the central star. The detected CO is not well shielded against photodissociation, and must thus be continuously replenished by gas release from exocomets within the belt. We analyze in detail the process of exocometary gas release and destruction around young M dwarfs and how this process compares to earlier type stars. Taking these differences into account, we find that CO generation through exocometary gas release naturally explains the increasing CO detection rates with stellar luminosity, mostly because the CO production rate from the collisional cascade is directly proportional to stellar luminosity. More luminous stars will therefore on average host more massive (and hence more easily detectable) exocometary CO disks, leading to the higher detection rates observed. The current CO detection rates are consistent with a ubiquitous release of exocometary gas in planetesimal belts, independent of spectral type.


title: Phenotypic redshifts with self-organizing maps: A novel method to characterize redshift distributions of source galaxies for weak lensing

abstract: Wide-field imaging surveys such as the Dark Energy Survey (DES) rely on coarse measurements of spectral energy distributions in a few filters to estimate the redshift distribution of source galaxies. In this regime, sample variance, shot noise, and selection effects limit the attainable accuracy of redshift calibration and thus of cosmological constraints. We present a new method to combine wide-field, few-filter measurements with catalogs from deep fields with additional filters and sufficiently low photometric noise to break degeneracies in photometric redshifts. The multi-band deep field is used as an intermediary between wide-field observations and accurate redshifts, greatly reducing sample variance, shot noise, and selection effects. Our implementation of the method uses self-organizing maps to group galaxies into phenotypes based on their observed fluxes, and is tested using a mock DES catalog created from N-body simulations. It yields a typical uncertainty on the mean redshift in each of five tomographic bins for an idealized simulation of the DES Year 3 weak-lensing tomographic analysis of $\sigma_{\Delta z} = 0.007$, which is a 60% improvement compared to the Year 1 analysis. Although the implementation of the method is tailored to DES, its formalism can be applied to other large photometric surveys with a similar observing strategy.


title: Formation of Hot Jupiters through Secular Chaos and Dynamical Tides

abstract: The population of giant planets on short-period orbits can potentially be explained by some flavours of high-eccentricity migration. In this paper we investigate one such mechanism involving "secular chaos", in which secular interactions between at least three giant planets push the inner planet to a highly eccentric orbit, followed by tidal circularization and orbital decay. In addition to the equilibrium tidal friction, we incorporate dissipation due to dynamical tides that are excited inside the giant planet. Using the method of Gaussian rings to account for planet-planet interactions, we explore the conditions for extreme eccentricity excitation via secular chaos and the properties of hot Jupiters formed in this migration channel. Our calculations show that once the inner planet reaches a sufficiently large eccentricity, dynamical tides quickly dissipate the orbital energy, producing an eccentric warm Jupiter, which then decays in semi-major axis through equilibrium tides to become a hot Jupiter. Dynamical tides help the planet avoid tidal disruption, increasing the chance of forming a hot Jupiter, although not all planets survive the process. We find that the final orbital periods generally lie in the range of 2-3 days, somewhat shorter than those of the observed hot Jupiter population. We couple the planet migration to the stellar spin evolution to predict the final spin-orbit misalignments. The distribution of the misalignment angles we obtain shows a lack of retrograde orbits compared to observations. Our results suggest that high-eccentricity migration via secular chaos can only account for a fraction of the observed hot Jupiter population.


title: Does radiative feedback make faint z>6 galaxies look small?

abstract: Recent observations of lensed sources have shown that the faintest ($M_{\mathrm{UV}} \approx -15\,\mathrm{mag}$) galaxies observed at z=6-8 appear to be extremely compact. Some of them have inferred sizes of less than 40 pc for stellar masses between $10^6$ and $10^7\,\mathrm{M}_{\odot}$, comparable to individual super star clusters or star cluster complexes at low redshift. High-redshift, low-mass galaxies are expected to show a clumpy, irregular morphology and if star clusters form in each of these well-separated clumps, the observed galaxy size would be much larger than the size of an individual star forming region. As supernova explosions impact the galaxy with a minimum delay time that exceeds the time required to form a massive star cluster, other processes are required to explain the absence of additional massive star forming regions. In this work we investigate whether the radiation of a young massive star cluster can suppress the formation of other detectable clusters within the same galaxy already before supernova feedback can affect the galaxy. We find that in low-mass ($M_{200} \lesssim 10^{10}\,\mathrm{M}_{\odot}$) haloes, the radiation from a compact star forming region with an initial mass of $10^{7}\,\mathrm{M}_{\odot}$ can keep gas clumps with Jeans masses larger than $\approx 10^{7}\,\mathrm{M}_{\odot}$ warm and ionized throughout the galaxy. In this picture, the small intrinsic sizes measured in the faintest $z=6-8$ galaxies are a natural consequence of the strong radiation field that stabilises massive gas clumps. A prediction of this mechanism is that the escape fraction for ionizing radiation is high for the extremely compact, high-z sources.


title: Path integrals for the relativistic particle: Some conceptual and pedagogical comments

abstract: In my textbook on Quantum Field Theory \cite{tpqft} and in a recent paper \cite{tpejc2018}, I advocated a lattice regularization procedure for defining the path integral for the relativistic particle, using the non-quadratic action containing a square root. I also provided an interpretation of this result in terms of the Jacobi action principle. This note clarifies several conceptual and pedagogical issues related to this approach and highlights some interesting open questions which this result leads to.


title: Coherent-State-Based Twin-Field Quantum Key Distribution

abstract: Large-scale quantum communication networks are still a huge challenge due to the rate-distance limit of quantum key distribution (QKD). Recently, twin-field (TF) QKD has been proposed to overcome this limit. Here, we prove that coherent-state-based TF-QKD is a time-reversed entanglement protocol, where the entanglement generation is realized with entanglement swapping operation via an entangled coherent state measurement. We propose a coherent-state-based TF-QKD with optimal secret key rate under symmetric and asymmetric channels by using coherent state and cat state coding. Furthermore, we show that our protocol can be converted to all recent coherent-state-based TF-QKD protocols by using our security proof. By using the entanglement purification with two-way classical communication, we improve the transmission distance of all coherent-state-based TF-QKD protocols.


title: Predicting Growth Rate from Gene Expression

abstract: Growth rate is one of the most important and most complex phenotypic characteristics of unicellular microorganisms, which determines the genetic mutations that dominate at the population level, and ultimately whether the population will survive. Translating changes at the genetic level to their growth rate consequences remains a subject of intense interest, since such a mapping could rationally direct experiments to optimize antibiotic efficacy or bioreactor productivity. In this paper, we directly map transcriptional profiles to growth rates by gathering published gene-expression data from Escherichia coli and Saccharomyces cerevisiae with corresponding growth-rate measurements. Using a machine-learning technique called k-nearest-neighbors regression, we build a model which predicts growth rate from gene expression. By exploiting the correlated nature of gene expression and sparsifying the model, we capture 81% of the variance in growth rate of the E. coli dataset while reducing the number of features from over 4,000 to nine. In S. cerevisiae, we account for 89% of the variance in growth rate while reducing from over 5,500 dimensions to 18. Such a model provides a basis for selecting successful strategies from among the combinatorial number of experimental possibilities when attempting to optimize complex phenotypic traits like growth rate.


title: Detecting Ocean Glint on Exoplanets Using Multiphase Mapping

abstract: Rotational mapping and glint are two proposed methods to directly detect liquid water on the surface of habitable exoplanets. However, false positives for both methods may prevent the unambiguous detection of exoplanet oceans. We use simulations of Earth as an exoplanet to introduce a combination of multiwavelength, multiphase, time-series direct-imaging observations and accompanying analyses that may improve the robustness of exoplanet ocean detection by spatially mapping the ocean glint signal. As the planet rotates, the glint spot appears to "blink" as Lambertian scattering continents interrupt the specular reflection from the ocean. This manifests itself as a strong source of periodic variability in crescent-phase reflected light curves. We invert these light curves to constrain the longitudinal slice maps and apparent albedo of two surfaces at both quadrature and crescent phase. At crescent phase, the retrieved apparent albedo of ocean-bearing longitudinal slices is increased by a factor of 5, compared to the albedo at quadrature phase, due to the contribution from glint. The land-bearing slices exhibit no significant change in apparent albedo with phase. The presence of forward-scattering clouds in our simulated observation increases the overall reflectivity toward crescent, but clouds do not correlate with any specific surfaces, thereby allowing for the phase-dependent glint effect to be interpreted as distinct from cloud scattering. Retrieving the same longitudinal map at quadrature and crescent phases may be used to tie changes in the apparent albedo with phase back to specific geographic surfaces, although this requires ideal geometries. We estimate that crescent-phase time-dependent glint detections are feasible for between 1-10 habitable zone exoplanets orbiting the nearest G, K, and M dwarfs using a space-based, high-contrast, direct-imaging telescope with a diameter >6 m.


title: Ultraviolet HST Spectroscopy of Planck Cold Clumps

abstract: We report results of the first study utilizing the ultraviolet capabilities of the Hubble Space Telescope to investigate a sample of Planck Galactic Cold Clump (PGCC) sources. We have selected high-resolution spectra toward 25 stars that contain a multitude of interstellar absorption lines associated with the interstellar medium (ISM) gas within these PGCC sources, including carbon monoxide (CO), C I and O I. By building cloud-component models of the individual absorption components present in these spectra, we can identify and isolate components associated with the PGCC sources, allowing for a more accurate investigation of the ISM behavior within these sources. Despite probing a broad range of overall sightline properties, we detect CO along each sightline. Sightlines with CO column density N(CO)$~>~$10$^{15}~$cm$^{-2}$ exhibit spatial dependence in N(CO) and CO/C I, while sightlines with N(CO)$~<~$10$^{15}~$cm$^{-2}$ show no such spatial dependence. Differences between N(H$_2$) values derived from UV absorption and dust emission suggest structure in the spatial distribution of N(H$_2$), where "CO-bright" sightlines are associated with PGCC sources embedded within smooth translucent envelopes, and "CO-dark" sightlines are associated with PGCC sources embedded in patchier environments containing more diffuse gas.


title: Deciphering an evolutionary sequence of merger stages in infrared-luminous starburst galaxies at z ~ 0.7

abstract: Based on optical/near-IR Magellan FIRE spectra of 25 starburst galaxies at 0.5 < z < 0.9, Calabr\`o et al.(2018) showed that their attenuation properties can be explained by a single-parameter sequence of total obscurations ranging from A(V)=2 to A(V)=30 towards the starburst core centers in a mixed stars and dust configuration. We investigate here the origin of this sequence for the same sample. We show that total attenuations anti-correlate with the starburst sizes in radio (3 GHz) with a significance larger than 5sigma and a scatter of 0.26 dex. More obscured and compact starbursts also show enhanced N2 (=[NII]/Halpha) ratios and larger line velocity widths that we attribute to an increasing shock contribution toward later merger phases, driven by deeper gravitational potential wells at the coalescence. Additionally, the attenuation is also linked to the equivalent width (EW) of hydrogen recombination lines, which is sensitive to the luminosity weighted age of the relatively unobscured stellar populations. Overall, the correlations among A(V), radio size, line width, N2 and EW of Balmer/Paschen lines converge towards suggesting an evolutionary sequence of merger stages: all of these quantities are likely to be good time-tracers of the merger phenomenon, and their large spanned range appears to be characteristic of the different merger phases. Half of our sample at higher obscurations have radio sizes approximately 3 times smaller than early type galaxies at the same redshift, suggesting that, in analogy with local Ultraluminous Infrared galaxies (ULIRGs), these cores cannot be directly forming elliptical galaxies. Finally, we detect mid-IR AGN torus for half of our sample and additional X-ray emission for 6 starbursts; intriguingly, the latter have systematically more compact sizes, suggestive of emerging AGNs towards later merger stages, possibly precursors of a later QSO phase.


title: A mass-dependent slope of the galaxy size-mass relation out to z~3: further evidence for a direct relation between median galaxy size and median halo mass

abstract: We reassess the galaxy size-mass relation out to z~3 using a new definition of size and a sample of >29,000 galaxies from the 3D-HST, CANDELS, and COSMOS-DASH surveys. Instead of the half-light radius r_50 we use r_80, the radius containing 80% of the stellar light. We find that the r_80 -- M_star relation has the form of a broken power law, with a clear change of slope at a pivot mass M_p. Below the pivot mass the relation is shallow (r_80 \propto M_star^0.15) and above it it is steep (r_80\propto M_star^0.6). The pivot mass increases with redshift, from log(M_p/M_sun)~ 10.2 at z=0.4 to log(M_p/M_sun)~ 10.9 at z=1.7-3. We compare these r_80-M_star relations to the M_halo-M_star relations derived from galaxy-galaxy lensing, clustering analyses, and abundance matching techniques. Remarkably, the pivot stellar masses of both relations are consistent with each other at all redshifts, and the slopes are very similar both above and below the pivot when assuming M_halo \propto r_80^3. The implied scaling factor to relate galaxy size to halo size is r_80 / R_vir = 0.047, independent of stellar mass and redshift.From redshift 0 to 1.5, the pivot mass also coincides with the mass where the fraction of star-forming galaxies is 50%, suggesting that the pivot mass reflects a transition from dissipational to dissipationless galaxy growth. Finally, our results imply that the scatter in the stellar-to-halo mass ratio is relatively small for massive halos (~0.2 dex for M_halo>10^12.5 M_sun).


title: Speed of qubit states during thermalisation

abstract: Classifying quantum states usually demands to observe properties such as the amount of correlation at one point in time. Further insight may be gained by inspecting the dynamics in a given evolution scheme. Here we attempt such a classification looking at single-qubit and two-qubit states at the start of thermalisation with a heat bath. The speed with which the evolution starts is influenced by quantum aspects of the state, however, such signatures do not allow for a systematic classification.


title: A Low-Flux State in IRAS 00521-7054 seen with NuSTAR and XMM-Newton: Relativistic Reflection and an Ultrafast Outflow

abstract: We present results from a deep, coordinated $XMM$-$Newton$+$NuSTAR$ observation of the Seyfert 2 galaxy IRAS 00521-7054. The $NuSTAR$ data provide the first detection of this source in high-energy X-rays ($E > 10$ keV), and the broadband data show this to be a highly complex source which exhibits relativistic reflection from the inner accretion disc, further reprocessing by more distant material, neutral absorption, and evidence for ionised absorption in an extreme, ultrafast outflow ($v_{\rm{out}} \sim 0.4c$). Based on lamppost disc reflection models, we find evidence that the central supermassive black hole is rapidly rotating ($a > 0.77$), consistent with previous estimates from the profile of the relativistic iron line, and that the accretion disc is viewed at a fairly high inclination ($i \sim 59^{\circ}$). Based on extensive simulations, we find the ultrafast outflow is detected at $\sim$4$\sigma$ significance (or greater). We also estimate that the extreme outflow should be sufficient to power galaxy-scale feedback, and may even dominate the energetics of the total output from the system.


title: A New View of the Size-Mass Distribution of Galaxies: Using $r_{20}$ and $r_{80}$ instead of $r_{50}$

abstract: When investigating the sizes of galaxies it is standard practice to use the half-light radius, $r_{50}$. Here we explore the effects of the size definition on the distribution of galaxies in the size -- stellar mass plane. Specifically, we consider $r_{20}$ and $r_{80}$, the radii that contain 20% and 80% of a galaxy's total luminosity, as determined from a Sersic profile fit, for galaxies in the 3D-HST/CANDELS and COSMOS-DASH surveys. These radii are calculated from size catalogs based on a simple calculation assuming a Sersic profile. We find that the size-mass distributions for $r_{20}$ and $r_{80}$ are markedly different from each other and also from the canonical $r_{50}$ distribution. The most striking difference is in the relative sizes of star forming and quiescent galaxies at fixed stellar mass. Whereas quiescent galaxies are smaller than star forming galaxies in $r_{50}$, this difference nearly vanishes for $r_{80}$. By contrast, the distance between the two populations increases for $r_{20}$. Considering all galaxies in a given stellar mass and redshift bin we detect a significant bimodality in the distribution of $r_{20}$, with one peak corresponding to star forming galaxies and the other to quiescent galaxies. We suggest that different measures of the size are tracing different physical processes within galaxies; $r_{20}$ is closely related to processes controlling the star formation rate of galaxies and $r_{80}$ may be sensitive to accretion processes and the relation of galaxies with their halos.


title: A circumbinary protoplanetary disc in a polar configuration

abstract: Nearly all young stars are initially surrounded by `protoplanetary' discs of gas and dust, and in the case of single stars at least 30\% of these discs go on to form planets. The process of protoplanetary disc formation can result in initial misalignments, where the disc orbital plane is different to the stellar equator in single star systems, or to the binary orbital plane in systems with two stars. A quirk of the dynamics means that initially misaligned `circumbinary' discs -- those that surround two stars -- are predicted to evolve to one of two possible stable configurations, one where the disc and binary orbital planes are coplanar, and one where they are perpendicular (a `polar' configuration). Prior work has found coplanar circumbinary discs, but no polar examples were known until now. Here we report the first discovery of a protoplanetary circumbinary disc in the polar configuration, supporting the predictions that such discs should exist. The disc shows some characteristics that are similar to discs around single stars, and that are attributed to dust growth. Thus, the first stages of planet formation appear able to proceed in polar circumbinary discs.


title: Decorrelating the errors of the galaxy correlation function with compact transformation matrices

abstract: Covariance matrix estimation is a persistent challenge for cosmology, often requiring a large number of synthetic mock catalogues. The off-diagonal components of the covariance matrix also make it difficult to show representative error bars on the 2-point correlation function (2PCF), since errors computed from the diagonal values of the covariance matrix greatly underestimate the uncertainties. We develop a routine for decorrelating the projected and anisotropic 2PCF with simple and scale-compact transformations on the 2PCF. These transformation matrices are modeled after the Cholesky decomposition and the symmetric square root of the Fisher matrix. Using mock catalogues, we show that the transformed projected and anisotropic 2PCF recover the same structure as the original 2PCF, while producing largely decorrelated error bars. Specifically, we propose simple Cholesky based transformation matrices that suppress the off-diagonal covariances on the projected 2PCF by ~95% and that on the anisotropic 2PCF by ~87%. These transformations also serve as highly regularized models of the Fisher matrix, compressing the degrees of freedom so that one can fit for the Fisher matrix with a much smaller number of mocks.


title: The Galaxy Stellar Mass Function and Low Surface Brightness Galaxies from Core-Collapse Supernovae

abstract: We introduce a method for producing a galaxy sample unbiased by surface brightness and stellar mass, by selecting star-forming galaxies via the positions of core-collapse supernovae (CCSNe). Whilst matching $\sim$2400 supernovae from the SDSS-II Supernova Survey to their host galaxies using IAC Stripe 82 legacy coadded imaging, we find $\sim$150 previously unidentified low surface brightness galaxies (LSBGs). Using a sub-sample of $\sim$900 CCSNe, we infer CCSN-rate and star-formation rate densities as a function of galaxy stellar mass, and the star-forming galaxy stellar mass function. Resultant star-forming galaxy number densities are found to increase following a power-law down to our low mass limit of $\sim10^{6.4}$ M$_{\odot}$ by a single Schechter function with a faint-end slope of $\alpha = -1.41$. Number densities are consistent with those found by the EAGLE simulations invoking a $\Lambda$-CDM cosmology. Overcoming surface brightness and stellar mass biases is important for assessment of the sub-structure problem. In order to estimate galaxy stellar masses, a new code for the calculation of galaxy photometric redshifts, zMedIC, is also presented, and shown to be particularly useful for small samples of galaxies.


title: Revealing the dust grain size in the inner envelope of the Class I protostar Per-emb-50

abstract: A good constraint of when the growth of dust grains from sub-micrometer to millimeter sizes occurs, is crucial for planet formation models. This provides the first step towards the production of pebbles and planetesimals in protoplanetary disks. Currently, it is well established that Class II objects have large dust grains. However, it is not clear when in the star formation process this grain growth occurs. We use multi-wavelength millimeter observations of a Class I protostar to obtain the spectral index of the observed flux densities $\alpha_\mathrm{mm}$ of the unresolved disk and the surrounding envelope. Our goal is to compare our observational results with visibility modeling at both wavelengths simultaneously. We present data from NOEMA at 2.7 mm and SMA at 1.3 mm of the Class I protostar, Per-emb-50. We model the dust emission with a variety of parametric and radiative transfer models to deduce the grain size from the observed emission spectral index. We find a spectral index in the envelope of Per-emb-50 of $\alpha_{\rm env}$=$3.3\pm0.3$, similar to the typical ISM values. The radiative transfer modeling of the source confirms this value of $\alpha_{\rm env}$ with the presence of dust with a $a_\mathrm{max}$$\leq$100 $\mu$m. Additionally, we explore the backwarming effect, where we find that the envelope structure affects the millimeter emission of the disk. Our results reveal grains with a maximum size no larger than $100$ $\mu$m in the inner envelope of the Class I protostar Per-emb-50, providing an interesting case to test the universality of millimeter grain growth expected in these sources.


title: Vision Zero: on a Provable Method for Eliminating Roadway Accidents without Compromising Traffic Throughput

abstract: We propose an economical, viable, approach to eliminate almost all car accidents. Our method relies on a mathematical model of safety and can be applied to all modern cars at a mild cost.


title: $\lambda \phi^4$ Theory II: The Broken Phase Beyond NNNN(NNNN)LO

abstract: We extend the study of the two-dimensional euclidean $\phi^4$ theory initiated in ref. [1] to the $\mathbb Z_2$ broken phase. In particular, we compute in perturbation theory up to N$^4$LO in the quartic coupling the vacuum energy, the vacuum expectation value of $\phi$ and the mass gap of the theory. We determine the large order behavior of the perturbative series by finding the leading order finite action complex instanton configuration in the $\mathbb Z_2$ broken phase. Using an appropriate conformal mapping, we then Borel resum the perturbative series. Interestingly enough, the truncated perturbative series for the vacuum energy and the vacuum expectation value of the field is reliable up to the critical coupling where a second order phase transition occurs, and breaks down around the transition for the mass gap. We compute the vacuum energy using also an alternative perturbative series, dubbed exact perturbation theory, that allows us to effectively reach N$^8$LO in the quartic coupling. In this way we can access the strong coupling region of the $\mathbb Z_2$ broken phase and test Chang duality by comparing the vacuum energies computed in three different descriptions of the same physical system. This result can also be considered as a confirmation of the Borel summability of the theory. Our results are in very good agreement (and with comparable or better precision) with those obtained by Hamiltonian truncation methods. We also discuss some subtleties related to the physical interpretation of the mass gap and provide evidence that the kink mass can be obtained by analytic continuation from the unbroken to the broken phase.

categories: hep-th cond-mat.stat-mech hep-lat hep-ph


title: Econophysics of Asset Price, Return and Multiple Expectations

abstract: This paper describes asset price and return disturbances as result of relations between transactions and multiple kinds of expectations. We show that disturbances of expectations can cause fluctuations of trade volume, price and return. We model price disturbances for transactions made under all types of expectations as weighted sum of partial price and trade volume disturbances for transactions made under separate kinds of expectations. Relations on price allow present return as weighted sum of partial return and trade volume "return" for transactions made under separate expectations. Dependence of price disturbances on trade volume disturbances as well as dependence of return on trade volume "return" cause dependence of volatility and statistical distributions of price and return on statistical properties of trade volume disturbances and trade volume "return" respectively.


title: Effect of random anisotropy on NMR frequency shift in the polar phase of the superfluid 3He

abstract: The orbital anisotropy induced in the superfluid 3He by nematic aerogel is, generally speaking, spatially non-uniform. This anisotropy in its turn induces spatial fluctuations of the order parameter. It is shown here that for the polar phase these fluctuations decrease the overall amplitude of the order parameter and the value of the frequency shift of the transverse NMR. Different contributions to this effect are discussed and estimated. Their temperature dependencies are discssed as well.


title: A study of the persistence and extinction of malaria in a family of epidemic models

abstract: This paper investigates the deterministic extinction and permanence of a family of SEIRS malaria models with multiple random delays, and with a general nonlinear incidence rate. The conditions for the extinction and permanence of the disease are presented. For the persistence of disease, improved analytical techniques are employed to obtain eventual lower bounds for the states of the system, while extinction results are obtained via analyzing the Lyapunov exponent. The basic reproduction number $R^{*}_{0}$ is calculated, and for $R^{*}_{0}<1$, extinction of disease occurs, while for $R^{*}_{0}\geq 1$, the disease persists. But if $R^{*}_{0}\geq 1$, and the expected survival probability of the plasmodium $E(e^{-(\mu_{v}T_{1}+\mu T_{2})})<\frac{1}{R^{*}_{0}}$, then extinction occurs. Results are interpreted, and numerical simulation examples are presented to (a) compare and evaluate the impacts of the intensity of the incidence rate, and the intensity of response to malaria control, on the persistence of malaria in the human population, and (2) to justify the asymptotic stability of the disease-free equilibrium.


title: Diagonal Subalgebras of Residual Intersections

abstract: Let ${\sf k}$ be a field, $S$ be a bigraded ${\sf k}$-algebra, and $S_\Delta$ denote the diagonal subalgebra of $S$ corresponding to $\Delta = \{ (cs,es) \; | \; s \in \mathbb{Z} \}$. It is know that the $S_\Delta$ is Koszul for $c,e \gg 0$. In this article, we find bounds for $c,e$ for $S_\Delta$ to be Koszul, when $S$ is a geometric residual intersection. Furthermore, we also study the Cohen-Macaulay property of these algebras. Finally, as an application, we look at classes of linearly presented perfect ideals of height two in a polynomial ring, show that all their powers have a linear resolution, and study the Koszul, and Cohen-Macaulay property of the diagonal subalgebras of their Rees algebras.


title: The Bayesian Prophet: A Low-Regret Framework for Online Decision Making

abstract: We develop a new framework for designing online policies given access to an oracle providing statistical information about an offline benchmark. Having access to such prediction oracles enables simple and natural Bayesian selection policies, and raises the question as to how these policies perform in different settings. Our work makes two important contributions towards this question: First, we develop a general technique we call *compensated coupling* which can be used to derive bounds on the expected regret (i.e., additive loss with respect to a benchmark) for any online policy and offline benchmark. Second, using this technique, we show that a natural greedy policy, which we call *the Bayes Selector*, has constant expected regret (i.e., independent of the number of arrivals and resource levels) for a large class of problems we refer to as Online Allocation with finite types, which includes widely-studied Online Packing and Online Matching problems. Our results generalize and simplify several existing results for Online Packing and Online Matching, and suggest a promising pathway for obtaining oracle-driven policies for other online decision-making settings.


title: Resolved young binary systems and their disks

abstract: We have conducted a survey of young single and multiple systems in the Taurus-Auriga star-forming region with the Atacama Large Millimeter Array (ALMA), substantially improving both the spatial resolution and sensitivity with which individual protoplanetary disks in these systems have been observed. These ALMA observations can resolve binary separations as small as 25--30 AU and have an average 3$\sigma$ detection level of 0.35 mJy, equivalent to a disk mass of $4 \times 10^{-5}$ M$_{\odot}$ for an M3 star. Our sample was constructed from stars that have an infrared excess and/or signs of accretion and have been classified as Class II. For the binary and higher order multiple systems observed, we detect $\lambda = 1.3$ mm continuum emission from one or more stars in all of our target systems. Combined with previous surveys of Taurus, our 21 new detections increase the fraction of millimeter-detected disks to over 75% in all categories of stars (singles, primaries, and companions) earlier than spectral type M6 in the Class II sample. Given the wealth of other information available for these stars, this has allowed us to study the impact of multiplicity with a much larger sample. While millimeter flux and disk mass are related to stellar mass as seen in previous studies, we find that both primary and secondary stars in binary systems with separations of 30 to 4200 AU have lower values of millimeter flux as a function of stellar mass than single stars. We also find that for these systems, the circumstellar disk around the primary star does not dominate the total disk mass in the system and contains on average 62% of the total mass.


title: Achlys : Towards a framework for distributed storage and generic computing applications for wireless IoT edge networks with Lasp on GRiSP

abstract: Internet of Things (IoT) has gained substantial attention over the past years. And the main discussion has been how to process the amount of data that it generates which has lead to the edge computing paradigm. Wether it is called fog1, edge or mist, the principle remains that cloud services must become available closer to clients. This documents presents ongoing work on future edge systems that are built to provide steadfast IoT services to users by bringing storage and processing power closer to peripheral parts of networks. Designing such infrastructures is becoming much more challenging as the number of IoT devices keeps growing. Production grade deployments have to meet very high performance requirements, and end-to-end solutions involve significant investments. In this paper, we aim at providing a solution to extend the range of the edge model to the very farthest nodes in the network. Specifically, we focus on providing reliable storage and computation capabilities immediately on wireless IoT sensor nodes. This extended edge model will allow end users to manage their IoT ecosystem without forcibly relying on gateways or Internet provider solutions. In this document, we introduce Achlys, a prototype implementation of an edge node that is a concrete port of the Lasp programming library on the GRiSP Erlang embedded system. This way, we aim at addressing the need for a general purpose edge that is both resilient and consistent in terms of storage and network. Finally, we study example use cases that could take advantage of integrating the Achlys framework and discuss future work for the latter.


title: Analysis and algorithms for $\ell_p$-based semi-supervised learning on graphs

abstract: This paper addresses theory and applications of $\ell_p$-based Laplacian regularization in semi-supervised learning. The graph $p$-Laplacian for $p>2$ has been proposed recently as a replacement for the standard ($p=2$) graph Laplacian in semi-supervised learning problems with very few labels, where Laplacian learning is degenerate. In the first part of the paper we prove new discrete to continuum convergence results for $p$-Laplace problems on $k$-nearest neighbor ($k$-NN) graphs, which are more commonly used in practice than random geometric graphs. Our analysis shows that, on $k$-NN graphs, the $p$-Laplacian retains information about the data distribution as $p\to \infty$ and Lipschitz learning ($p=\infty$) is sensitive to the data distribution. This situation can be contrasted with random geometric graphs, where the $p$-Laplacian forgets the data distribution as $p\to \infty$. We also present a general framework for proving discrete to continuum convergence results in graph-based learning that only requires pointwise consistency and monotonicity. In the second part of the paper, we develop fast algorithms for solving the variational and game-theoretic $p$-Laplace equations on weighted graphs for $p>2$. We present several efficient and scalable algorithms for both formulations, and present numerical results on synthetic data indicating their convergence properties. Finally, we conduct extensive numerical experiments on the MNIST, FashionMNIST and EMNIST datasets that illustrate the effectiveness of the $p$-Laplacian formulation for semi-supervised learning with few labels. In particular, we find that Lipschitz learning ($p=\infty$) performs well with very few labels on $k$-NN graphs, which experimentally validates our theoretical findings that Lipschitz learning retains information about the data distribution (the unlabeled data) on $k$-NN graphs.


title: A dynamical model of remote-control model cars

abstract: Simple experiments for which differential equations cannot be solved analytically can be addressed using an effective model that satisfactorily reproduces the experimental data. In this work, the one-dimensional kinematics of a remote-control model (toy) car was studied experimentally and its dynamical equation modelled. In the experiment, maximum power was applied to the car, initially at rest, until it reached its terminal velocity. Digital video recording was used to obtain the relevant kinematic variables that enabled to plot trajectories in the phase space. A dynamical equation of motion was proposed in which the overall frictional force was modelled as an effective force proportional to the velocity raised to the power of a real number. Since such an equation could not be solved analytically, a dynamical model was developed and the system parameters were calculated by non-linear fitting. Finally, the resulting values were substituted in the motion equation and the numerical results thus obtained were compared with the experimental data, corroborating the accuracy of the model.


title: Linear stability analysis of transient electrodeposition in charged porous media: suppression of dendritic growth by surface conduction

abstract: We study the linear stability of transient electrodeposition in a charged random porous medium, whose pore surface charges can be of any sign, flanked by a pair of planar metal electrodes. Discretization of the linear stability problem results in a generalized eigenvalue problem for the dispersion relation that is solved numerically, which agrees well with the analytical approximation obtained from a boundary layer analysis valid at high wavenumbers. Under galvanostatic conditions in which an overlimiting current is applied, in the classical case of zero surface charges, the electric field at the cathode diverges at Sand's time due to electrolyte depletion. The same phenomenon happens for positive charges but earlier than Sand's time. However, negative charges allow the system to sustain an overlimiting current via surface conduction past Sand's time, keeping the electric field bounded. Therefore, at Sand's time, negative charges greatly reduce surface instabilities and suppress dendritic growth, while zero and positive charges magnify them. We compare theoretical predictions for overall surface stabilization with published experimental data for copper electrodeposition in cellulose nitrate membranes and demonstrate good agreement between theory and experiment. We also apply the stability analysis to how crystal grain size varies with duty cycle during pulse electroplating.


title: Top-antitop charge asymmetry measurements in the lepton+jets channel with the ATLAS detector

abstract: We report a measurement of the charge asymmetry in top quark pair production with the ATLAS experiment at the LHC using 20.3 fb$^{-1}$ of pp collision data at $\sqrt{s}$ = 8 TeV. The data are unfolded using a fully Bayesian unfolding method. Differential measurements of the charge asymmetry are performed.


title: An informal introduction to quantitative stochastic homogenization

abstract: Divergence-form operators with random coefficients homogenize over large scales. Over the last decade, an intensive research effort focused on turning this asymptotic statement into quantitative estimates. The goal of this note is to review one approach for doing so based on the idea of renormalization. The discussion is highly informal, with pointers to mathematically precise statements.


title: Decay of periodic entropy solutions to degenerate nonlinear parabolic equations

abstract: Under a precise nonlinearity-diffusivity condition we establish the decay of space-periodic entropy solutions of a multidimensional degenerate nonlinear parabolic equation.


title: The Finite Horizon impulse control Problem with arbitrary cost functions : the Viscosity Solution Approach

abstract: We consider stochastic impulse control problems when the impulses cost functions are arbitrary. We use the dynamic programming principle and viscosity solutions approach to show that the value function is a unique viscosity solution for the associated Hamilton-Jacobi-Bellman equation (HJB) partial differential equation (PDE) of stochastic impulse control problems


title: 3-D Radiation Mapping in Real-Time with the Localization and Mapping Platform LAMP from Unmanned Aerial Systems and Man-Portable Configurations

abstract: Real-time, meter-resolution gamma-ray mapping is relevant in the detection and mapping of radiological materials, and for applications ranging from nuclear decommissioning, waste management, and environmental remediation to homeland security, emergency response, and international safeguards. We present the Localization and Mapping Platform (LAMP) as a modular, contextual and radiation detector sensor suite, which performs gamma-ray mapping in three dimensions (3-D) and in real time, onboard an unmanned aerial vehicle (UAV) or in a man-portable configuration. The deployment of an unmanned aerial system (UAS) for gamma-ray mapping can be advantageous, as the UAS provides a means of measuring large areas efficiently and improving accessibility to some environments, such as multi-story structures. In addition, it is possible to increase measurement robustness through autonomous navigation, and to reduce radiation exposure to users as a result of the remote measurement. LAMP enables meter-resolution gamma-ray mapping through Scene Data Fusion (SDF) [1], a capability that fuses radiation and scene data via voxelized 3-D Maximum Likelihood Expectation Maximization (MLEM) to produce 3-D maps of radioactive source distributions in real-time. Results are computed onboard LAMP while it is flying on the UAV and streamed from the system to the user, who can view the 3-D map on a tablet in real-time. We present results that demonstrate the SDF concept, including a set of UAS flights where a 133Ba source is localized at a test site in Berkeley, CA and a handheld measurement in Fukushima Prefecture, Japan where the distribution of radiocesium(137,134Cs) released from the accident of the Fukushima Daiichi Nuclear Power Plant is mapped. The reconstruction parameters used for each measurement were identical, indicating that the same algorithm can be used for both point or distributed sources.


title: Local symmetry rank bound for positive intermediate Ricci curvatures

abstract: We use a local argument to prove if an $r$-dimensional torus acts isometrically and effectively on a connected $n$-dimensional manifold which has positive $k^\mathrm{th}$-intermediate Ricci curvature at some point, then $r \leq \lfloor \frac{n+k}{2} \rfloor$. This symmetry rank bound generalizes those established by Grove and Searle for positive sectional curvature and Wilking for quasipositive curvature. As a consequence, we show that the symmetry rank bound in the Maximal Symmetry Rank Conjecture for manifolds of non-negative sectional curvature holds for those which also have positive intermediate Ricci curvature at some point. In the process of proving our symmetry rank bound, we also obtain an optimal dimensional restriction on isometric immersions of manifolds with non-positive intermediate Ricci curvature into manifolds with positive intermediate Ricci curvature, generalizing a result by Otsuki.


title: Anomalous Hall Conductivity of a Non-Collinear Magnetic Antiperovskite

abstract: The anomalous Hall effect (AHE) is a well-known fundamental property of ferromagnetic metals, commonly associated with the presence of a net magnetization. Recently, an AHE has been discovered in non-collinear antiferromagnetic (AFM) metals. Driven by non-vanishing Berry curvature of AFM materials with certain magnetic space group symmetry, anomalous Hall conductivity (AHC) is very sensitive to the specific type of magnetic ordering. Here, we investigate the appearance of AHC in antiperovskite GaNMn$_{3}$ as a representative of broader materials family ANMn$_{3}$ (A is a main group element), where different types of non-collinear magnetic ordering can emerge. Using symmetry analyses and first-principles density-functional theory calculations, we show that with almost identical band structure, the nearly degenerate non-collinear AFM $\Gamma_{5g}$ and $\Gamma_{4g}$ phases of GaNMn$_{3}$ have zero and finite AHC, respectively. In a non-collinear ferrimagnetic $M$-1 phase, GaNMn$_{3}$ exhibits a large AHC due to the presence of a sizable net magnetic moment. In the non-collinear antiperovskite magnets, transitions between different magnetic phases, exhibiting different AHC states, can be produced by doping, strain, or spin transfer torque, which makes these materials promising for novel spintronic applications.


title: Answering Comparative Questions: Better than Ten-Blue-Links?

abstract: We present CAM (comparative argumentative machine), a novel open-domain IR system to argumentatively compare objects with respect to information extracted from the Common Crawl. In a user study, the participants obtained 15% more accurate answers using CAM compared to a "traditional" keyword-based search and were 20% faster in finding the answer to comparative questions.


title: Imaging with flat optics: metalenses or diffractive lenses?

abstract: Recently, there has been an explosion of interest in metalenses for imaging. The interest is primarily based on their sub-wavelength thicknesses. Diffractive lenses have been used as thin lenses since the late 19th century. Here, we show that multi-level diffractive lenses (MDLs), when designed properly can exceed the performance of metalenses. Furthermore, MDLs can be designed and fabricated with larger constituent features, making them accessible to low-cost, large area volume manufacturing, which is generally challenging for metalenses. The support substrate will dominate overall thickness for all flat optics. Therefore the advantage of a slight decrease in thickness (from ~2{\lambda} to ~{\lambda}/2) afforded by metalenses may not be useful. We further elaborate on the differences between these approaches and clarify that metalenses have unique advantages when manipulating the polarization states of light.


title: On complexity of branching droplets in electrical field

abstract: Decanol droplets in a thin layer of sodium decanoate with sodium chloride exhibit bifurcation branching growth due to interplay between osmotic pressure, diffusion and surface tension. We aimed to evaluate if morphology of the branching droplets changes when the droplets are subject to electrical potential difference. We analysed graph-theoretic structure of the droplets and applied several complexity measures. We found that, in overall, the current increases complexity of the branching droplets in terms of number of connected components and nodes in their graph presentations, morphological complexity and compressibility.


title: A linear programming approach to the tracking of partials

abstract: A new approach to the tracking of sinusoidal chirps using linear programming is proposed. It is demonstrated that the classical algorithm of McAulay and Quatieri is greedy and exhibits exponential complexity for long searches, while approaches based on the Viterbi algorithm exhibit factorial complexity. A linear programming (LP) formulation to find the best $L$ paths in a lattice is described and its complexity is shown to be less than previous approaches. Finally it is demonstrated that the new LP formulation outperforms the classical algorithm in the tracking of sinusoidal chirps in high levels of noise.


title: Solving inverse problems via auto-encoders

abstract: Compressed sensing (CS) is about recovering a structured signal from its under-determined linear measurements. Starting from sparsity, recovery methods have steadily moved towards more complex structures. Emerging machine learning tools such as generative functions that are based on neural networks are able to learn general complex structures from training data. This makes them potentially powerful tools for designing CS algorithms. Consider a desired class of signals $\cal Q$, ${\cal Q}\subset{R}^n$, and a corresponding generative function $g:{\cal U}^k\rightarrow {R}^n$, ${\cal U}\subset {R}$, such that $\sup_{{\bf x}\in {\cal Q}}\min_{{\bf u}\in{\cal U}^k}{1\over \sqrt{n}}\|g({\bf u})-{\bf x}\|\leq \delta$. A recovery method based on $g$ seeks $g({\bf u})$ with minimum measurement error. In this paper, the performance of such a recovery method is studied, under both noisy and noiseless measurements. In the noiseless case, roughly speaking, it is proven that, as $k$ and $n$ grow without bound and $\delta$ converges to zero, if the number of measurements ($m$) is larger than the input dimension of the generative model ($k$), then asymptotically, almost lossless recovery is possible. Furthermore, the performance of an efficient iterative algorithm based on projected gradient descent is studied. In this case, an auto-encoder is used to define and enforce the source structure at the projection step. The auto-encoder is defined by encoder and decoder (generative) functions $f:{R}^n\to{\cal U}^k$ and $g:{\cal U}^k\to{R}^n$, respectively. We theoretically prove that, roughly, given $m>40k\log{1\over \delta}$ measurements, such an algorithm converges to the vicinity of the desired result, even in the presence of additive white Gaussian noise. Numerical results exploring the effectiveness of the proposed method are presented.


title: Boundaries of relative factor graphs and subgroup classification for automorphisms of free products

abstract: Given a countable group $G$ splitting as a free product $G=G_1\ast\dots\ast G_k\ast F_N$, we establish classification results for subgroups of the group $Out(G,\mathcal{F})$ of all outer automorphisms of $G$ that preserve the conjugacy classes of each $G_i$. We show that every finitely generated subgroup $H\subseteq Out(G,\mathcal{F})$ either contains a relatively fully irreducible automorphism, or else it virtually preserves the conjugacy class of a proper free factor relative to the decomposition (the finite generation hypothesis on $H$ can be dropped for $G=F_N$, or more generally when $G$ is toral relatively hyperbolic). In the first case, either $H$ virtually preserves a nonperipheral conjugacy class in $G$, or else $H$ contains an atoroidal automorphism. The key geometric tool to obtain these classification results is a description of the Gromov boundaries of relative versions of the free factor graph $\mathrm{FF}$ and the $\mathcal{Z}$-factor graph $\mathcal{Z}\mathrm{F}$, as spaces of equivalence classes of arational trees (respectively relatively free arational trees). We also identify the loxodromic isometries of $\mathrm{FF}$ with the fully irreducible elements of $Out(G,\mathcal{F})$, and loxodromic isometries of $\mathcal{Z}\mathrm{F}$ with the fully irreducible atoroidal outer automorphisms.


title: Disorder-induced exceptional points and nodal lines in Dirac superconductors

abstract: We consider the effect of disorder on the spectrum of quasiparticles in the point-node and nodal-line superconductors. Due to the anisotropic dispersion of quasiparticles disorder scattering may render the Hamiltonian describing these excitations non-Hermitian. Depending on the dimensionality of the system, we show that the nodes in the spectrum are replaced by Fermi arcs or Fermi areas bounded by exceptional points or exceptional lines, respectively. These features are illustrated by first considering a model of a proximity-induced superconductor in an anisotropic two-dimensional (2D) Dirac semimetal, where a Fermi arc in the gap bounded by exceptional points can be realized. We next show that the interplay between disorder and supercurrents can give rise to a 2D Fermi surface bounded by exceptional lines in three-dimensional (3D) nodal superconductors.


title: Plurisuperharmonicity of reciprocal energy function on Teichmuller space and Weil-Petersson metrics

abstract: We consider harmonic maps$u(z): \mathcal{X}_z\to N$ in a fixed homotopy class from Riemann surfaces $\mathcal{X}_z$ of genus $g\geq 2$ varying in the Teichm\"u{}ller space $\mathcal T$ to a Riemannian manifold $N$ with non-positive Hermitian sectional curvature. The energy function $E(z)=E(u(z))$ can be viewed as a function on $\mathcal T$ and we study its first and the second variations. We prove that the reciprocal energy function $E(z)^{-1}$ is plurisuperharmonic on Teichm\"uller space. We also obtain the (strict) plurisubharmonicity of $\log E(z)$ and $E(z)$. As an application, we get the following relationship between the second variation of logarithmic energy function and the Weil-Petersson metric if the harmonic map $u(z)$ is holomorphic or anti-holomorphic and totally geodesic, i.e., $$ \sqrt{-1}\p\b{\p}\log E(z)=\frac{\omega_{WP}}{2\pi(g-1)}. $$ We consider also the energy function $E(z)$ associated to the harmonic maps from a fixed compact K\"ahler manifold $M$ to Riemann surfaces ${\mathcal{X}_z\}_{z\in\mathcal{T}}$ in a fixed homotopy class. If $u(z)$ is holomorphic or anti-holomorphic, then the above equation is also proved.


title: Bonseyes AI Pipeline -- bringing AI to you. End-to-end integration of data, algorithms and deployment tools

abstract: Next generation of embedded Information and Communication Technology (ICT) systems are collaborative systems able to perform autonomous tasks. The remarkable expansion of the embedded ICT market, together with the rise and breakthroughs of Artificial Intelligence (AI), have put the focus on the Edge as it stands as one of the keys for the next technological revolution: the seamless integration of AI in our daily life. However, training and deployment of custom AI solutions on embedded devices require a fine-grained integration of data, algorithms, and tools to achieve high accuracy. Such integration requires a high level of expertise that becomes a real bottleneck for small and medium enterprises wanting to deploy AI solutions on the Edge which, ultimately, slows down the adoption of AI on daily-life applications. In this work, we present a modular AI pipeline as an integrating framework to bring data, algorithms, and deployment tools together. By removing the integration barriers and lowering the required expertise, we can interconnect the different stages of tools and provide a modular end-to-end development of AI products for embedded devices. Our AI pipeline consists of four modular main steps: i) data ingestion, ii) model training, iii) deployment optimization and, iv) the IoT hub integration. To show the effectiveness of our pipeline, we provide examples of different AI applications during each of the steps. Besides, we integrate our deployment framework, LPDNN, into the AI pipeline and present its lightweight architecture and deployment capabilities for embedded devices. Finally, we demonstrate the results of the AI pipeline by showing the deployment of several AI applications such as keyword spotting, image classification and object detection on a set of well-known embedded platforms, where LPDNN consistently outperforms all other popular deployment frameworks.


title: The Use of Activity Trackers for Health Empowerment and Commitment: The Philippine Cycling Perspective

abstract: Activity tracking devices have found its way in the world of cycling. With its projected market demand and increasing popularity of cycling in the Philippines, cyclists are slowly adopting this technology in their daily cycling routines. Activity trackers demonstrate real-time data which allow cyclists to adjust physical efforts to achieve their personal goals. Six common features of activity trackers were formed as constructs to explore its influence on health empowerment in the context of cycling using Partial Least Squares Structural Equation Model. A total of 393 cyclists in the Philippines participated in the study. Some features demonstrated strong evidence of positive influence in achieving health empowerment and commitment. Implications for future design and development of this technology device are discussed.


title: Machine-learning a virus assembly fitness landscape

abstract: Realistic evolutionary fitness landscapes are notoriously difficult to construct. A recent cutting-edge model of virus assembly consists of a dodecahedral capsid with $12$ corresponding packaging signals in three affinity bands. This whole genome/phenotype space consisting of $3^{12}$ genomes has been explored via computationally expensive stochastic assembly models, giving a fitness landscape in terms of the assembly efficiency. Using latest machine-learning techniques by establishing a neural network, we show that the intensive computation can be short-circuited in a matter of minutes to astounding accuracy.


title: Classification of load forecasting studies by forecasting problem to select load forecasting techniques and methodologies

abstract: The key contribution of this paper is to propose a classification into two dimensions of the load forecasting studies to decide which forecasting tools to use in which case. This classification aims to provide a synthetic view of the relevant forecasting techniques and methodologies by forecasting problem. In addition, the key principles of the main techniques and methodologies used are summarized along with the reviews of these papers. The classification process relies on two couples of parameters that define a forecasting problem. Each article is classified with key information about the dataset used and the forecasting tools implemented: the forecasting techniques (probabilistic or deterministic) and methodologies, the data cleansing techniques, and the error metrics. The process to select the articles reviewed in this paper was conducted into two steps. First, a set of load forecasting studies was built based on relevant load forecasting reviews and forecasting competitions. The second step consisted in selecting the most relevant studies of this set based on the following criteria: the quality of the description of the forecasting techniques and methodologies implemented, the description of the results, and the contributions. This paper can be read in two passes. The first one by identifying the forecasting problem of interest to select the corresponding class into one of the four classification tables. Each one references all the articles classified across a forecasting horizon. They provide a synthetic view of the forecasting tools used by articles addressing similar forecasting problems. Then, a second level composed of four Tables summarizes key information about the forecasting tools and the results of these studies. The second pass consists in reading the key principles of the main techniques and methodologies of interest and the reviews of the articles.


title: An Agent-Based Model to Explain the Emergence of Stylised Facts in Log Returns

abstract: This paper outlines an agent-based model of a simple financial market in which a single asset is available for trade by three different types of traders. The model was first introduced in the PhD thesis of one of the authors, see reference [1]. The simulated log returns are examined for the presence of the stylised facts of financial data. The features of leptokurtosis, volatility clustering and aggregational Gaussianity are especially highlighted and studied in detail. The following ingredients are found to be essential for the production of these stylised facts: the memory of noise traders who make random trade decisions; the inclusion of technical traders that trade in line with trends in the price and the inclusion of fundamental traders who know the "fundamental value" of the stock and trade accordingly. When these three basic types of traders are included log returns are produced with a leptokurtic distribution and volatility clustering as well as some further statistical features of empirical data. This enhances and broadens our understanding of the fundamental processes involved in the production of empirical data by the market.


title: Solution of the equation $y^{\prime}=f(y)$ and Bell Polynomials

abstract: In this paper we use Faa di Bruno's formula to associate Bell polynomial values to differential equations of the form $y^{\prime}=f(y)$. That is, we use partial Bell polynomials to represent the solution of such an equation and use the solution to compute special values of partial Bell polynomials.


title: Sextic double solids with Artin-Mumford obstructions to rationality

abstract: We study a double solid X branched along a nodal sextic surface in a projective space and the 2-torsion subgroup in the third integer cohomology group of a resolution of singularities of X. This group can be considered as an obstruction to rationality of X. Studying this group we conclude that all sextic double solids admitting non-trivial obstructions to rationality are branched along determinantal surfaces of very specific type and we provide an explicit list of them.


title: A nonparametric super-efficient estimator of the average treatment effect

abstract: Doubly robust estimators of causal effects are a popular means of estimating causal effects. Such estimators combine an estimate of the conditional mean of the outcome given treatment and confounders (the so-called outcome regression) with an estimate of the conditional probability of treatment given confounders (the propensity score) to generate an estimate of the effect of interest. In addition to enjoying the double-robustness property, these estimators have additional benefits. First, flexible regression tools, such as those developed in the field of machine learning, can be utilized to estimate the relevant regressions, while the estimators of the treatment effects retain desirable statistical properties. Furthermore, these estimators are often statistically efficient, achieving the lower bound on the variance of regular, asymptotically linear estimators. However, in spite of their asymptotic optimality, in problems where causal estimands are weakly identifiable, these estimators may behave erratically. We propose two new estimation techniques for use in these challenging settings. Our estimators build on two existing frameworks for efficient estimation: targeted minimum loss estimation and one-step estimation. However, rather than using an estimate of the propensity score in their construction, we instead opt for an alternative regression quantity when building our estimators: the conditional probability of treatment given the conditional mean outcome. We discuss the theoretical implications and demonstrate the estimators' performance in simulated and real data.


title: Constraining Type Ia Supernova Asymmetry with the Gamma-Ray Escape Timescale

abstract: We calculate the effects of an asymmetric $^{56}\mathrm{Ni}$ distribution in Type Ia supernova (SN Ia) ejecta on the gamma-ray escape timescale ($t_0$) that characterizes the late light curve ($>40$ days after peak) and find the effect is modest compared to other possible variations in ejecta structure. We parameterize asymmetry in the $^{56}\mathrm{Ni}$ distribution and calculate $t_0$ for a grid of SN ejecta models spanning a large volume of the asymmetry parameter space. The models have spherical density profiles while the $^{56}\mathrm{Ni}$ distribution in them has various levels of asymmetry. By placing constraints based on the observational measurement of $t_0$ and other general properties of SN Ia ejecta, we find the range of allowed asymmetry in the $^{56}\mathrm{Ni}$ distribution. Considering these constraints we find that some level of asymmetry in the distribution is not ruled out. However, models with a single ejecta mass and varying $^{56}\mathrm{Ni}$ distributions cannot explain the full range of observed $t_0$ values. This strengthens the claim that both Chandrasekhar mass and sub-Chandrasekhar mass explosions are required to explain the diversity of SN Ia observations.


title: Going back to the roots: how to create stroboscopic photos from digital videos

abstract: We show here how to create from digital films using the well-known software Tracker stroboscopic photos in order to analyze different types of movements. The advantage of this procedure is that it is possible to analyze the printed photo or on a computer screen in an intuitive way for the students. After presenting a historical perspective of the use of stroboscopic photos in secondary education we discuss several examples: the movement of a remote control car, an elastic planar collision, the movement of a projectile and also an experiment in electromagnetism, more specifically, the discharge of a capacitor measured using an analog multimeter.