title: Large-scale peculiar velocity fields: Newtonian vs relativistic treatment

abstract: We employ a perturbative analysis to study the evolution of large-scale peculiar velocity fields within the framework of Newtonian gravity and then compare our results to those of the corresponding relativistic treatment. In so doing, we use the same mathematical formalism and apply the same physical approach. This facilitates a direct and transparent comparison between the two treatments. Our study recovers and extends the familiar Newtonian results on the one hand, while on the other it shows that the Newtonian analysis leads to substantially weaker growth-rates for the peculiar velocity field, compared to the relativistic approach. This implies that, by using Newton's rather than Einstein's theory, one could seriously underestimate the overall kinematic evolution of cosmological peculiar motions. We are also in the position to identify the reason the two theories arrive at such considerably different results and conclusions.


title: MEGA: Merger graphs of structure formation

abstract: When following the growth of structure in the Universe, we propose replacing merger trees with merger graphs, in which haloes can both merge and split into separate pieces. We show that this leads to smoother mass growth and eliminates catastrophic failures in which massive haloes have no progenitors or descendants. For those who prefer to stick with merger trees, we find that trees derived from our merger graphs have similar mass growth properties to previous methods, but again without catastrophic failures. For future galaxy formation modelling, two different density thresholds can be used to distinguish host haloes (extended galactic haloes, groups and clusters) from higher-density subhaloes: sites of galaxy formation.


title: Quasinormal modes, stability and shadows of a black hole in the 4D Einstein-Gauss-Bonnet gravity

abstract: Recently a $D$-dimensional regularization approach leading to the non-trivial $(3+1)$-dimensional Einstein-Gauss-Bonnet (EGB) effective description of gravity was formulated which was claimed to bypass the Lovelock's theorem and avoid Ostrogradsky instability. Later it was shown that the regularization is possible only for some broad, but limited, class of metrics and Aoki, Gorji and Mukohyama [arXiv:2005.03859] formulated a well-defined four-dimensional EGB theory, which breaks the Lorentz invariance in a theoretically consistent and observationally viable way. The black-hole solution of the first naive approach proved out to be also the exact solution of the well-defined theory. Here we calculate quasinormal modes of scalar, electromagnetic and gravitational perturbations and find the radius of shadow for spherically symmetric and asymptotically flat black holes with Gauss-Bonnet corrections. We show that the black hole is gravitationally stable when ($-16 M^2<\alpha \lessapprox 0.6 M^2$). The instability in the outer range is the eikonal one and it develops at high multipole numbers. The radius of the shadow $R_{Sh}$ obeys the linear law with a remarkable accuracy.


title: A Szemer\'{e}di-type theorem for subsets of the unit cube

abstract: We investigate gaps of $n$-term arithmetic progressions $x, x+y, \ldots, x+(n-1)y$ inside a positive measure subset $A$ of the unit cube $[0,1]^d$. If lengths of their gaps $y$ are evaluated in the $\ell^p$-norm for any $p$ other than $1, 2, \ldots, n-1$, and $\infty$, and if the dimension $d$ is large enough, then we show that the numbers $\|y\|_{\ell^p}$ attain all values from an interval, the length of which depends only on $n$, $p$, $d$, and the measure of $A$. Known counterexamples prevent generalizations of this result to the remaining values of the exponent $p$. We also give an explicit bound for the length of the aforementioned interval. The proof makes the bound depend on the currently available bounds in Szemer\'{e}di's theorem on the integers, which are used as a black box. A key ingredient of the proof are power-type cancellation estimates for operators resembling the multilinear Hilbert transforms. As a byproduct of the approach we obtain a quantitative improvement of the corresponding (previously known) result for side lengths of $n$-dimensional cubes with vertices lying in a positive measure subset of $([0,1]^2)^n$.


title: Adversarial Generation of Informative Trajectories for Dynamics System Identification

abstract: Dynamic System Identification approaches usually heavily rely on the evolutionary and gradient-based optimisation techniques to produce optimal excitation trajectories for determining the physical parameters of robot platforms. Current optimisation techniques tend to generate single trajectories. This is expensive, and intractable for longer trajectories, thus limiting their efficacy for system identification. We propose to tackle this issue by using multiple shorter cyclic trajectories, which can be generated in parallel, and subsequently combined together to achieve the same effect as a longer trajectory. Crucially, we show how to scale this approach even further by increasing the generation speed and quality of the dataset through the use of generative adversarial network (GAN) based architectures to produce a large databases of valid and diverse excitation trajectories. To the best of our knowledge, this is the first robotics work to explore system identification with multiple cyclic trajectories and to develop GAN-based techniques for scaleably producing excitation trajectories that are diverse in both control parameter and inertial parameter spaces. We show that our approach dramatically accelerates trajectory optimisation, while simultaneously providing more accurate system identification than the conventional approach.


title: On the Temperature Dependence of Grain Boundary Mobility

abstract: The grain boundary (GB) mobility relates GB velocity to the thermodynamic driving forces and is central to our understanding of microstructure evolution in polycrystals. Recent molecular dynamics (MD) and experimental studies have shown that the temperature-dependence of the GB mobility is much more varied than is commonly thought. GB mobility may increase, decrease, remain constant or show multiple peaks with increasing temperature. We propose a mechanistic model for GB migration, based on the formation and migration of line defects (disconnection) within the GB. We implement this model in a kinetic Monte Carlo and statistical mechanics framework; the results capture all of these observed temperature dependencies and are shown to be in quantitative agreement with each other and direct MD simulations of GB migration for a set of specific grain boundaries. Examination of the dependence of GB mobility on disconnection mode and temperature provides new insight into how GBs migrate in polycrystalline materials.


title: Persistence probabilities of weighted sums of stationary Gaussian sequences

abstract: With $\{\xi_i\}_{i\ge 0}$ being a centered stationary Gaussian sequence with non-negative correlation function $\rho(i):=\mathbb{E}[ \xi_0\xi_i]$ and $\{\sigma(i)\}_{i\ge 1}$ a sequence of positive reals, we study the asymptotics of the persistence probability of the weighted sum $\sum_{i=1}^\ell \sigma(i) \xi_i$, $\ell\ge 1$. For summable correlations $\rho$, we show that the persistence exponent is universal. On the contrary, for non-summable $\rho$, even for polynomial weight functions $\sigma(i)\sim i^p$ the persistence exponent depends on the rate of decay of the correlations (encoded by a parameter $H$) and on the polynomial rate $p$ of $\sigma$. In this case, we show existence of the persistence exponent $\theta(H,p)$ and study its properties as a function of $(p,H)$. During the course of our proofs, we develop several tools for dealing with exit problems for Gaussian processes with non-negative correlations -- e.g.\ a continuity result for persistence exponents and a necessary and sufficient criterion for the persistence exponent to be zero -- that might be of independent interest.


title: Approximate modularity: Kalton's constant is not smaller than 3

abstract: Kalton and Roberts [Trans. Amer. Math. Soc., 278 (1983), 803--816] proved that there exists a universal constant $K\leqslant 44.5$ such that for every set algebra $\mathcal{F}$ and every 1-additive function $f\colon \mathcal{F}\to \mathbb R$ there exists a finitely-additive signed measure $\mu$ defined on $\mathcal{F}$ such that $|f(A)-\mu(A)|\leqslant K$ for any $A\in \mathcal{F}$. The only known lower bound for the optimal value of $K$ was found by Pawlik [Colloq. Math., 54 (1987), 163--164], who proved that this constant is not smaller than $1.5$; we improve this bound to $3$ already on a non-negative 1-additive function.


title: A soft robotic tongue to develop solutions to manage swallowing disorders

abstract: The development of novel texture adaptations for the management of swallowing disorders could be accelerated if reliable in vitro tests were made available. This study addresses some of the limitations of swallowing in vitro models, by introducing a soft robotic actuator inspired by the tongue. The wettability of the soft-robotic actuator was engineered to achieve physiologically relevant contact angles to allow comparing dry and lubricated conditions. The actuator design and the control algorithm are designed to offer flexibility in the swallowing patterns to consider different scenarios, including poor lingual coordination. In vitro swallowing tests with shear thinning liquids were performed as a proof of concept and showed physiologically relevant oral transit time, bolus velocity and palatal pressure. Both the bolus rheology and the coordination of the peristaltic contractions influenced the temporal evolution of the bolus velocity and the bolus transit time. This novel soft robotic tongue, its integration in an in vitro oral cavity and its flexible control can therefore contribute to developing better solutions to manage swallowing disorders.


title: Homotopy lifting property for actions of finite abelian groups on Hausdorff topological spaces

abstract: In this paper we prove the homotopy lifting property for actions of finite abelian groups on Hausdorff topological spaces.


title: DEEVA: A Deep Learning and IoT Based Computer Vision System to Address Safety and Security of Production Sites in Energy Industry

abstract: When it comes to addressing the safety/security related needs at different production/construction sites, accurate detection of the presence of workers, vehicles, equipment important and formed an integral part of computer vision-based surveillance systems (CVSS). Traditional CVSS systems focus on the use of different computer vision and pattern recognition algorithms overly reliant on manual extraction of features and small datasets, limiting their usage because of low accuracy, need for expert knowledge and high computational costs. The main objective of this paper is to provide decision makers at sites with a practical yet comprehensive deep learning and IoT based solution to tackle various computer vision related problems such as scene classification, object detection in scenes, semantic segmentation, scene captioning etc. Our overarching goal is to address the central question of What is happening at this site and where is it happening in an automated fashion minimizing the need for human resources dedicated to surveillance. We developed Deep ExxonMobil Eye for Video Analysis (DEEVA) package to handle scene classification, object detection, semantic segmentation and captioning of scenes in a hierarchical approach. The results reveal that transfer learning with the RetinaNet object detector is able to detect the presence of workers, different types of vehicles/construction equipment, safety related objects at a high level of accuracy (above 90%). With the help of deep learning to automatically extract features and IoT technology to automatic capture, transfer and process vast amount of realtime images, this framework is an important step towards the development of intelligent surveillance systems aimed at addressing myriads of open ended problems in the realm of security/safety monitoring, productivity assessments and future decision making.


title: Learning to Collide: An Adaptive Safety-Critical Scenarios Generating Method

abstract: Long-tail and rare event problems become crucial when autonomous driving algorithms are applied in the real world. For the purpose of evaluating systems in challenging settings, we propose a generative framework to create safety-critical scenarios for evaluating specific task algorithms. We first represent the traffic scenarios with a series of autoregressive building blocks and generate diverse scenarios by sampling from the joint distribution of these blocks. We then train the generative model as an agent (or a generator) to investigate the risky distribution parameters for a given driving algorithm being evaluated. We regard the task algorithm as an environment (or a discriminator) that returns a reward to the agent when a risky scenario is generated. Through the experiments conducted on several scenarios in the simulation, we demonstrate that the proposed framework generates safety-critical scenarios more efficiently than grid search or human design methods. Another advantage of this method is its adaptiveness to the routes and parameters.


title: Archaeoastronomical study of the historic churches of La Gomera

abstract: In this paper we discuss the importance of studying the orientation of ancient Christian churches, as a complement to the historical and cultural research of the temples. We present preliminary results of the analysis of the precise spatial orientation of 38 colonial churches located on the Canary Island of La Gomera (Spain). The sample suggests that, although several churches have a canonical orientation within the solar range, a large proportion of them follow a pattern of orientations that is compatible with the steep orography of the island and, therefore, contrasts with the prescriptions contained in the texts of the early Christian writers.


title: Structural and Magnetic Properties of Molecular Beam Epitaxy Grown Chromium Selenide Thin Films

abstract: Chromium selenide thin films were grown epitaxially on Al${_2}$O${_3}$(0001) and Si(111)-(7${\times}$7) substrates using molecular beam epitaxy (MBE). Sharp streaks in reflection high-energy electron diffraction and triangular structures in scanning tunneling microscopy indicate a flat smooth film growth along the c-axis, and is very similar to that from a hexagonal surface. X-ray diffraction pattern confirms the growth along the c-axis with c-axis lattice constant of 17.39 {\AA}. The grown film is semiconducting, having a small band gap of about 0.034 eV, as calculated from the temperature dependent resistivity. Antiferromagnetic nature of the film with a N\'eel temperature of about 40 K is estimated from the magnetic exchange bias measurements. A larger out-of-plane exchange bias, along with a smaller in-plane exchange bias is observed below 40 K. Exchange bias training effects are analyzed based on different models and are observed to be following a modified power-law decay behavior.


title: Natural Language Processing Advancements By Deep Learning: A Survey

abstract: Natural Language Processing (NLP) helps empower intelligent machines by enhancing a better understanding of the human language for linguistic-based human-computer communication. Recent developments in computational power and the advent of large amounts of linguistic data have heightened the need and demand for automating semantic analysis using data-driven approaches. The utilization of data-driven strategies is pervasive now due to the significant improvements demonstrated through the usage of deep learning methods in areas such as Computer Vision, Automatic Speech Recognition, and in particular, NLP. This survey categorizes and addresses the different aspects and applications of NLP that have benefited from deep learning. It covers core NLP tasks and applications and describes how deep learning methods and models advance these areas. We further analyze and compare different approaches and state-of-the-art models.


title: Spectral Bayesian Estimation for General Stochastic Hybrid Systems

abstract: General Stochastic Hybrid Systems (GSHS) have been formulated to represent various types of uncertainties in hybrid dynamical systems. In this paper, we propose computational techniques for Bayesian estimation of GSHS. In particular, the Fokker-Planck equation that describes the evolution of uncertainty distributions along GSHS is solved by spectral techniques, where an arbitrary form of probability density of the hybrid state is represented by a mixture of Fourier series. The method is based on splitting the Fokker-Planck equation represented by an integro-partial differential equation into the partial differentiation part for continuous diffusion and the integral part for discrete transition, and integrating the solution of each part. The propagated density function is used with a likelihood function in the Bayes' formula to estimate the hybrid state for given sensor measurements. The unique feature of the proposed technique is that the probability density describing complete stochastic properties of the hybrid state is constructed without relying on the common Gaussian assumption, in contrast to other methods that compute certain expected values such as mean or variance. We apply the proposed method to the estimation of two examples: the bouncing ball model and the Dubins vehicle model. We show that the proposed technique yields propagated densities consistent with Monte Carlo simulations, more accurate estimates over the Gaussian based approach and some computational benefits over a particle filter.


title: Toward evading the strong coupling problem in Horndeski genesis

abstract: It is of interest to understand whether or not one can construct a classical field theory description of early cosmology which would be free of the initial singularity and stable throughout the whole evolution. One of the known possibilities is genesis within the Horndeski theory, which is thought to be an alternative to or a possible completion of the inflationary scenario. In this model, the strong coupling energy scale tends to zero in the asymptotic past, $t \to - \infty$, making the model potentially intractable. We point out that despite the latter property, the classical setup may be trustworthy since the energy scale of the classical evolution (the inverse of its timescale) also vanishes as $t \to -\infty$. In the framework of a concrete model belonging to the Horndeski class, we show that the strong coupling energy scale of the cubic interactions vastly exceeds the classical energy scale in a certain range of parameters, indicating that the classical description is possible.


title: Concurrent Disjoint Set Union

abstract: We develop and analyze concurrent algorithms for the disjoint set union (union-find) problem in the shared memory, asynchronous multiprocessor model of computation, with CAS (compare and swap) or DCAS (double compare and swap) as the synchronization primitive. We give a deterministic bounded wait-free algorithm that uses DCAS and has a total work bound of $O(m \cdot (\log(np/m + 1) + \alpha(n, m/(np)))$ for a problem with $n$ elements and $m$ operations solved by $p$ processes, where $\alpha$ is a functional inverse of Ackermann's function. We give two randomized algorithms that use only CAS and have the same work bound in expectation. The analysis of the second randomized algorithm is valid even if the scheduler is adversarial. Our DCAS and randomized algorithms take $O(\log n)$ steps per operation, worst-case for the DCAS algorithm, high-probability for the randomized algorithms. Our work and step bounds grow only logarithmically with $p$, making our algorithms truly scalable. We prove that for a class of symmetric algorithms that includes ours, no better step or work bound is possible.


title: Energy-efficient and Robust Cumulative Training with Net2Net Transformation

abstract: Deep learning has achieved state-of-the-art accuracies on several computer vision tasks. However, the computational and energy requirements associated with training such deep neural networks can be quite high. In this paper, we propose a cumulative training strategy with Net2Net transformation that achieves training computational efficiency without incurring large accuracy loss, in comparison to a model trained from scratch. We achieve this by first training a small network (with lesser parameters) on a small subset of the original dataset, and then gradually expanding the network using Net2Net transformation to train incrementally on larger subsets of the dataset. This incremental training strategy with Net2Net utilizes function-preserving transformations that transfers knowledge from each previous small network to the next larger network, thereby, reducing the overall training complexity. Our experiments demonstrate that compared with training from scratch, cumulative training yields ~2x reduction in computational complexity for training TinyImageNet using VGG19 at iso-accuracy. Besides training efficiency, a key advantage of our cumulative training strategy is that we can perform pruning during Net2Net expansion to obtain a final network with optimal configuration (~0.4x lower inference compute complexity) compared to conventional training from scratch. We also demonstrate that the final network obtained from cumulative training yields better generalization performance and noise robustness. Further, we show that mutual inference from all the networks created with cumulative Net2Net expansion enables improved adversarial input detection.


title: The Solubility Graph Associated With a Finite Group

abstract: Let $G$ be a finite group. The solubility graph associated with the finite group $G$, denoted by $\Gamma_{\cal S}(G)$, is a simple graph whose vertices are the non-trivial elements of $G$, and there is an edge between two distinct elements $x$ and $y$ if and only if $\langle x, y\rangle$ is a soluble subgroup of $G$. In this paper, we examine some properties of solubility graphs.


title: A Scalar Parameterized Mechanism for Two-Sided Markets

abstract: We consider a market in which both suppliers and consumers compete for a product via scalar-parameterized supply offers and demand bids. Scalar-parameterized offers/bids are appealing due to their modeling simplicity and desirable mathematical properties with the most prominent being bounded efficiency loss and price markup under strategic interactions. Our model incorporates production capacity constraints and minimum inelastic demand requirements. Under perfect competition, the market mechanism yields allocations that maximize social welfare. When market participants are price-anticipating, we show that there exists a unique Nash equilibrium, and provide an efficient way to compute the resulting market allocation. Moreover, we explicitly characterize the bounds on the welfare loss and prices observed at the Nash equilibrium.


title: BARD: A structured technique for group elicitation of Bayesian networks to support analytic reasoning

abstract: In many complex, real-world situations, problem solving and decision making require effective reasoning about causation and uncertainty. However, human reasoning in these cases is prone to confusion and error. Bayesian networks (BNs) are an artificial intelligence technology that models uncertain situations, supporting probabilistic and causal reasoning and decision making. However, to date, BN methodologies and software require significant upfront training, do not provide much guidance on the model building process, and do not support collaboratively building BNs. BARD (Bayesian ARgumentation via Delphi) is both a methodology and an expert system that utilises (1) BNs as the underlying structured representations for better argument analysis, (2) a multi-user web-based software platform and Delphi-style social processes to assist with collaboration, and (3) short, high-quality e-courses on demand, a highly structured process to guide BN construction, and a variety of helpful tools to assist in building and reasoning with BNs, including an automated explanation tool to assist effective report writing. The result is an end-to-end online platform, with associated online training, for groups without prior BN expertise to understand and analyse a problem, build a model of its underlying probabilistic causal structure, validate and reason with the causal model, and use it to produce a written analytic report. Initial experimental results demonstrate that BARD aids in problem solving, reasoning and collaboration.


title: Unveiling of the mechanisms of acoustic streaming induced by sharp edges

abstract: Acoustic waves can generate steady streaming within a fluid owing to the generation of viscous boundary layers near walls, of typical thickness $\delta$. In microchannels, the acoustic wavelength $\lambda$ is adjusted to twice the channel width $w$ to ensure a resonance condition, which implies the use of MHz transducers. Recently though, intense acoustic streaming was generated by acoustic waves of a few kHz (hence with $\lambda \gg w$), owing to the presence of sharp-tipped structures of curvature radius at the tip $r_c$ smaller than $\delta$. The present study quantitatively investigates this sharp-edge acoustic streaming via the direct resolution of the full Navier-Stokes equation, using Finite Element Method. The influence of $\delta$, $r_c$ and viscosity $\nu$ on the acoustic streaming performance are quantified. Our results suggest choices of operating conditions and geometrical parameters, via dimensionless quantities $r_c/\delta$ and $\delta/w$ and provide guidelines on how to obtain strong, optimal sharp-edge acoustic streaming.


title: Log orthogonal functions: approximation properties and applications

abstract: We present two new classes of orthogonal functions, log orthogonal functions (LOFs) and generalized log orthogonal functions (GLOFs), which are constructed by applying a $\log$ mapping to Laguerre polynomials. We develop basic approximation theory for these new orthogonal functions and apply them to solve several typical fractional differential equations whose solutions exhibit weak singularities. Our error analysis and numerical results show that our methods based on the new orthogonal functions are particularly suitable for functions which have weak singularities at one endpoint, and can lead to exponential convergence rate, as opposed to low algebraic rates if usual orthogonal polynomials are used.


title: $p$-harmonic functions by way of intrinsic mean value properties

abstract: Let $\Omega\subset\mathbb{R}^n$ be a bounded domain satisfying the uniform exterior cone condition. We establish existence and uniqueness of continuous solutions of the Dirichlet Problem associated to certain intrinsic nonlinear mean value properties in $\Omega $. Furthermore we show that, when properly normalized, such functions converge to the $p$-harmonic solution of the Dirichlet problem in $\Omega$, for $p\in[2,\infty)$. The proof of existence is constructive and the methods are entirely analytic, a fundamental tool being the construction of explicit, $p$-independent barrier functions in $\Omega$.


title: A bijective proof of Kohnert's rule for Schubert polynomials

abstract: Kohnert proposed the first monomial positive formula for Schubert polynomials as the generating polynomial for certain unit cell diagrams obtained from the Rothe diagram of a permutation. Billey, Jockusch and Stanley gave the first proven formula for Schubert polynomials as the generating polynomial for compatible sequences of reduced words of a permutation. In this paper, we give an explicit bijection between these two models, thereby definitively proving Kohnert's rule for Schubert polynomials.


title: SrCu$_2$(BO$_3$)$_2$ under pressure: a first-principle study

abstract: Using density-functional theory (DFT) band-structure calculations, we study the crystal structure, the lattice dynamics, and the magnetic interactions in the Shastry-Sutherland magnet SrCu$_2$(BO$_3$)$_2$ under pressure, concentrating on experimentally relevant pressures up to 4 GPa. All tendencies inside the candidate intermediate phases are thoroughly worked out, including specific predictions for some Raman active phonon modes that could be used to pin down the nature of the intermediate phase.


title: Extrinsic Calibration of a 3D-LIDAR and a Camera

abstract: This work presents an extrinsic parameter estimation algorithm between a 3D LIDAR and a Projective Camera using a marker-less planar target, by exploiting Planar Surface Point to Plane and Planar Edge Point to back-projected Plane geometric constraints. The proposed method uses the data collected by placing the planar board at different poses in the common field of view of the LIDAR and the Camera. The steps include, detection of the target and the edges of the target in LIDAR and Camera frames, matching the detected planes and lines across both the sensing modalities and finally solving a cost function formed by the aforementioned geometric constraints that link the features detected in both the LIDAR and the Camera using non-linear least squares. We have extensively validated our algorithm using two Basler Cameras, Velodyne VLP-32 and Ouster OS1 LIDARs.


title: Inferring high-resolution human mixing patterns for disease modeling

abstract: Mathematical and computational modeling approaches are increasingly used as quantitative tools in the analysis and forecasting of infectious disease epidemics. The growing need for realism in addressing complex public health questions is however calling for accurate models of the human contact patterns that govern the disease transmission processes. Here we present a data-driven approach to generate effective descriptions of population-level contact patterns by using highly detailed macro (census) and micro (survey) data on key socio-demographic features. We produce age-stratified contact matrices for 277 sub-national administrative regions of countries covering approximately 3.5 billion people and reflecting the high degree of cultural and societal diversity of the focus countries. We use the derived contact matrices to model the spread of airborne infectious diseases and show that sub-national heterogeneities in human mixing patterns have a marked impact on epidemic indicators such as the reproduction number and overall attack rate of epidemics of the same etiology. The contact patterns derived here are made publicly available as a modeling tool to study the impact of socio-economic differences and demographic heterogeneities across populations on the epidemiology of infectious diseases.


title: Lagrangian-Eulerian Multi-Density Topology Optimization with the Material Point Method

abstract: In this paper, a hybrid Lagrangian-Eulerian topology optimization (LETO) method is proposed to solve the elastic force equilibrium with the Material Point Method (MPM). LETO transfers density information from freely movable Lagrangian carrier particles to a fixed set of Eulerian quadrature points. This transfer is based on a smooth radial kernel involved in the compliance objective to avoid the artificial checkerboard pattern. The quadrature points act as MPM particles embedded in a lower-resolution grid and enable a sub-cell multi-density resolution of intricate structures with a reduced computational cost. A quadrature-level connectivity graph-based method is adopted to avoid the artificial checkerboard issues commonly existing in multi-resolution topology optimization methods. Numerical experiments are provided to demonstrate the efficacy of the proposed approach.


title: High Performance Parallel Sort for Shared and Distributed Memory MIMD

abstract: We present four high performance hybrid sorting methods developed for various parallel platforms: shared memory multiprocessors, distributed multiprocessors, and clusters taking advantage of existence of both shared and distributed memory. Merge sort, known for its stability, is used to design several of our algorithms. We improve its parallel performance by combining it with Quicksort. We present two models designed for shared memory MIMD (OpenMP): (a) a non-recursive Merge sort and (b) a hybrid Quicksort and Merge sort. The third model presented is designed for distributed memory MIMD (MPI) using a hybrid Quicksort and Merge sort. Our fourth model is designed to take advantage of the shared memory within individual nodes of todays cluster systems, and to eliminate all internal data transfers between different nodes, Our model implements a one-step MSD-Radix to distribute data in ten packets (MPI) while parallel cores of each node use Quicksort to sort their data partitions sequentially then merge and sort them in parallel employing the OpenMp. The performances of all developed models outperform the baseline performance. Hybrid Quicksort and Merge sort outperformed Hybrid Memory Parallel Merge Sort using Hybrid MSD-Radix and Quicksort in Cluster Platforms when sorting small size data, but with larger data the speedup of Hybrid Memory Parallel Merge Sort Using Hybrid MSD-Radix and Quicksort in Cluster Platforms becomes bigger and it keeps improving. The speedup of Distributed Memory Parallel Hybrid Quicksort and Merge Sort is the best.


title: MRI Super-Resolution with GAN and 3D Multi-Level DenseNet: Smaller, Faster, and Better

abstract: High-resolution (HR) magnetic resonance imaging (MRI) provides detailed anatomical information that is critical for diagnosis in the clinical application. However, HR MRI typically comes at the cost of long scan time, small spatial coverage, and low signal-to-noise ratio (SNR). Recent studies showed that with a deep convolutional neural network (CNN), HR generic images could be recovered from low-resolution (LR) inputs via single image super-resolution (SISR) approaches. Additionally, previous works have shown that a deep 3D CNN can generate high-quality SR MRIs by using learned image priors. However, 3D CNN with deep structures, have a large number of parameters and are computationally expensive. In this paper, we propose a novel 3D CNN architecture, namely a multi-level densely connected super-resolution network (mDCSRN), which is light-weight, fast and accurate. We also show that with the generative adversarial network (GAN)-guided training, the mDCSRN-GAN provides appealing sharp SR images with rich texture details that are highly comparable with the referenced HR images. Our results from experiments on a large public dataset with 1,113 subjects showed that this new architecture outperformed other popular deep learning methods in recovering 4x resolution-downgraded images in both quality and speed.


title: A First Step Towards Understanding Real-world Attacks on IoT Devices

abstract: With the rapid growth of Internet of Things (IoT) devices, it is imperative to proactively understand the real-world cybersecurity threats posed to them. This paper describes our initial efforts towards building a honeypot ecosystem as a means to gathering and analyzing real attack data against IoT devices. A primary condition for a honeypot to yield useful insights is to let attackers believe they are real systems used by humans and organizations. IoT devices pose unique challenges in this respect, due to the large variety of device types and the physical-connectedness nature. We thus create a multiphased approach in building a honeypot ecosystem, where researchers can gradually increase a low-interaction honeypot's sophistication in emulating an IoT device by observing real-world attackers' behaviors. We deployed honeypots both on-premise and in the cloud, with associated analysis and vetting infrastructures to ensure these honeypots cannot be easily identified as such and appear to be real systems. In doing so we were able to attract increasingly sophisticated attack data. We present the design of this honeypot ecosystem and our observation on the attack data so far. Our data shows that real-world attackers are explicitly going after IoT devices, and some captured activities seem to involve direct human interaction (as opposed to scripted automatic activities). We also build a low interaction honeypot for IoT cameras, called Honeycamera, that present to attackers seemingly real videos. This is our first step towards building a more comprehensive honeypot ecosystem that will allow researchers to gain concrete understanding of what attackers are going after on IoT devices, so as to more proactively protect them.


title: Exactly Computing the Local Lipschitz Constant of ReLU Networks

abstract: The local Lipschitz constant of a neural network is a useful metric with applications in robustness, generalization, and fairness evaluation. We provide novel analytic results relating the local Lipschitz constant of nonsmooth vector-valued functions to a maximization over the norm of the generalized Jacobian. We present a sufficient condition for which backpropagation always returns an element of the generalized Jacobian, and reframe the problem over this broad class of functions. We show strong inapproximability results for estimating Lipschitz constants of ReLU networks, and then formulate an algorithm to compute these quantities exactly. We leverage this algorithm to evaluate the tightness of competing Lipschitz estimators and the effects of regularized training on the Lipschitz constant.


title: Semi-supervised learning of glottal pulse positions in a neural analysis-synthesis framework

abstract: This article investigates into recently emerging approaches that use deep neural networks for the estimation of glottal closure instants (GCI). We build upon our previous approach that used synthetic speech exclusively to create perfectly annotated training data and that had been shown to compare favourably with other training approaches using electroglottograph (EGG) signals. Here we introduce a semi-supervised training strategy that allows refining the estimator by means of an analysis-synthesis setup using real speech signals, for which GCI ground truth does not exist. Evaluation of the analyser is performed by means of comparing the GCI extracted from the glottal flow signal generated by the analyser with the GCI extracted from EGG on the CMU arctic dataset, where EGG signals were recorded in addition to speech. We observe that (1.) the artificial increase of the diversity of pulse shapes that has been used in our previous construction of the synthetic database is beneficial, (2.) training the GCI network in the analysis-synthesis setup allows achieving a very significant improvement of the GCI analyser, (3.) additional regularisation strategies allow improving the final analysis network when trained in the analysis-synthesis setup.


title: Graph covers with two new eigenvalues

abstract: A certain signed adjacency matrix of the hypercube, which Hao Huang used last year to resolve the sensitivity conjecture, is closely related to the unique, 4-cycle free, 2-fold cover of the hypercube. We develop a framework in which this connection is a natural first example of the relationship between group labeled adjacency matrices with few eigenvalues, and combinatorially interesting covering graphs. In particular, we define a two-eigenvalue cover to be a covering graph whose adjacency spectra differs (as a multiset) from that of the graph it covers by exactly two eigenvalues. We show that walk regularity of a graph implies walk regularity of any abelian two-eigenvalue cover. We also give a spectral characterization for when a cyclic two-eigenvalue cover of a strongly-regular graph is distance regular.


title: The Koala: A Fast Blue Optical Transient with Luminous Radio Emission from a Starburst Dwarf Galaxy at $z=0.27$

abstract: We present ZTF18abvkwla (the "Koala"), a fast blue optical transient discovered in the Zwicky Transient Facility (ZTF) One-Day Cadence (1DC) Survey. ZTF18abvkwla has a number of features in common with the groundbreaking transient AT2018cow: blue colors at peak ($g-r\approx-0.5$ mag), a short rise time from half-max of under two days, a decay time to half-max of only three days, a high optical luminosity ($M_{g,\mathrm{peak}}\approx-20.6$mag), a hot ($\gtrsim 40,000$K) featureless spectrum at peak light, and a luminous radio counterpart. At late times ($\Delta t>80$d) the radio luminosity of ZTF18abvkwla ($\nu L_\nu \gtrsim 10^{40}$erg/s at 10 GHz, observer-frame) is most similar to that of long-duration gamma-ray bursts (GRBs). The host galaxy is a dwarf starburst galaxy ($M\approx5\times10^{8}M_\odot$, $\mathrm{SFR}\approx7 M_\odot$/yr) that is moderately metal-enriched ($\log\mathrm{[O/H]} \approx 8.5$), similar to the hosts of GRBs and superluminous supernovae. As in AT2018cow, the radio and optical emission in ZTF18abvkwla likely arise from two separate components: the radio from fast-moving ejecta ($\Gamma \beta c >0.38c$) and the optical from shock-interaction with confined dense material ($<0.07M_\odot$ in $\sim 10^{15}$cm). Compiling transients in the literature with $t_\mathrm{rise} <5$d and $M_\mathrm{peak}<-20$mag, we find that a significant number are engine-powered, and suggest that the high peak optical luminosity is directly related to the presence of this engine. From 18 months of the 1DC survey, we find that transients in this rise-luminosity phase space are at least two to three orders of magnitude less common than CC SNe. Finally, we discuss strategies for identifying such events with future facilities like the Large Synoptic Survey Telescope, and prospects for detecting accompanying X-ray and radio emission.


title: A Deep learning Approach to Generate Contrast-Enhanced Computerised Tomography Angiography without the Use of Intravenous Contrast Agents

abstract: Contrast-enhanced computed tomography angiograms (CTAs) are widely used in cardiovascular imaging to obtain a non-invasive view of arterial structures. However, contrast agents are associated with complications at the injection site as well as renal toxicity leading to contrast-induced nephropathy (CIN) and renal failure. We hypothesised that the raw data acquired from a non-contrast CT contains sufficient information to differentiate blood and other soft tissue components. We utilised deep learning methods to define the subtleties between soft tissue components in order to simulate contrast enhanced CTAs without contrast agents. Twenty-six patients with paired non-contrast and CTA images were randomly selected from an approved clinical study. Non-contrast axial slices within the AAA from 10 patients (n = 100) were sampled for the underlying Hounsfield unit (HU) distribution at the lumen, intra-luminal thrombus and interface locations. Sampling of HUs in these regions revealed significant differences between all regions (p<0.001 for all comparisons), confirming the intrinsic differences in the radiomic signatures between these regions. To generate a large training dataset, paired axial slices from the training set (n=13) were augmented to produce a total of 23,551 2-D images. We trained a 2-D Cycle Generative Adversarial Network (cycleGAN) for this non-contrast to contrast (NC2C) transformation task. The accuracy of the cycleGAN output was assessed by comparison to the contrast image. This pipeline is able to differentiate between visually incoherent soft tissue regions in non-contrast CT images. The CTAs generated from the non-contrast images bear strong resemblance to the ground truth. Here we describe a novel application of Generative Adversarial Network for CT image processing. This is poised to disrupt clinical pathways requiring contrast enhanced CT imaging.


title: Photometry of Kuiper Belt Object (486958) Arrokoth from New Horizons LORRI

abstract: On January 1st 2019, the New Horizons spacecraft flew by the classical Kuiper belt object (486958) Arrokoth (provisionally designated 2014 MU69), possibly the most primitive object ever explored by a spacecraft. The I/F of Arrokoth is analyzed and fit with a photometric function that is a linear combination of the Lommel-Seeliger (lunar) and Lambert photometric functions. Arrokoth has a geometric albedo of p_V = 0.21_(-0.04)^(+0.05) at a wavelength of 550 nm and ~0.24 at 610 nm. Arrokoth's geometric albedo is greater than the median but consistent with a distribution of cold classical Kuiper belt objects whose geometric albedos were determined by fitting a thermal model to radiometric observations. Thus, Arrokoth's geometric albedo adds to the orbital and spectral evidence that it is a cold classical Kuiper belt object. Maps of the normal reflectance and hemispherical albedo of Arrokoth are presented. The normal reflectance of Arrokoth's surface varies with location, ranging from ~0.10-0.40 at 610 nm with an approximately Gaussian distribution. Both Arrokoth's extrema dark and extrema bright surfaces are correlated to topographic depressions. Arrokoth has a bilobate shape and the two lobes have similar normal reflectance distributions: both are approximately Gaussian, peak at ~0.25 at 610 nm, and range from ~0.10-0.40, which is consistent with co-formation and co-evolution of the two lobes. The hemispherical albedo of Arrokoth varies substantially with both incidence angle and location, the average hemispherical albedo at 610 nm is 0.063 +/- 0.015. The Bond albedo of Arrokoth at 610 nm is 0.062 +/- 0.015.


title: Fourier temporal ghost imaging

abstract: Ghost imaging is a fascinating framework which constructs the image of an object by correlating measurements between received beams and reference beams, none of which carries the structure information of the object independently. Recently, by taking into account space-time duality in optics, computational temporal ghost imaging has attracted attentions. Here, we propose a novel Fourier temporal ghost imaging (FTGI) scheme to achieve single-shot non-reproducible temporal signals. By sinusoidal coded modulation, ghost images are obtained and recovered by applying Fourier transformation. For demonstration, non-repeating events are detected with single-shot exposure architecture. It's shown in results that the peak signal-to-noise ratio (PSNR) of FTGI is significantly better (13dB increase) than traditional temporal ghost imaging in the same condition. In addition, by using the obvious physical meaning of Fourier spectrum, we show some potential applications of FTGI, such as frequency division multiplexing demodulation in the visible light communications.


title: Reachability Analysis for Feed-Forward Neural Networks using Face Lattices

abstract: Deep neural networks have been widely applied as an effective approach to handle complex and practical problems. However, one of the most fundamental open problems is the lack of formal methods to analyze the safety of their behaviors. To address this challenge, we propose a parallelizable technique to compute exact reachable sets of a neural network to an input set. Our method currently focuses on feed-forward neural networks with ReLU activation functions. One of the primary challenges for polytope-based approaches is identifying the intersection between intermediate polytopes and hyperplanes from neurons. In this regard, we present a new approach to construct the polytopes with the face lattice, a complete combinatorial structure. The correctness and performance of our methodology are evaluated by verifying the safety of ACAS Xu networks and other benchmarks. Compared to state-of-the-art methods such as Reluplex, Marabou, and NNV, our approach exhibits a significantly higher efficiency. Additionally, our approach is capable of constructing the complete input set given an output set, so that any input that leads to safety violation can be tracked.


title: Fast Predictive Uncertainty for Classification with Bayesian Deep Networks

abstract: In Bayesian Deep Learning, distributions over the output of classification neural networks are often approximated by first constructing a Gaussian distribution over the weights, then sampling from it to receive a distribution over the softmax outputs. This is costly. We reconsider old work (Laplace Bridge) to construct a Dirichlet approximation of this softmax output distribution, which yields an analytic map between Gaussian distributions in logit space and Dirichlet distributions (the conjugate prior to the Categorical distribution) in the output space. Importantly, the vanilla Laplace Bridge comes with certain limitations. We analyze those and suggest a simple solution that compares favorably to other commonly used estimates of the softmax-Gaussian integral. We demonstrate that the resulting Dirichlet distribution has multiple advantages, in particular, more efficient computation of the uncertainty estimate and scaling to large datasets and networks like ImageNet and DenseNet. We further demonstrate the usefulness of this Dirichlet approximation by using it to construct a lightweight uncertainty-aware output ranking for ImageNet.


title: The Use of Implicit Human Motor Behaviour in the Online Personalisation of Prosthetic Interfaces

abstract: In previous work, the authors proposed a data-driven optimisation algorithm for the personalisation of human-prosthetic interfaces, demonstrating the possibility of adapting prosthesis behaviour to its user while the user performs tasks with it. This method requires that the human and the prosthesis personalisation algorithm have same pre-defined objective function. This was previously ensured by providing the human with explicit feedback on what the objective function is. However, constantly displaying this information to the prosthesis user is impractical. Moreover, the method utilised task information in the objective function which may not be available from the wearable sensors typically used in prosthetic applications. In this work, the previous approach is extended to use a prosthesis objective function based on implicit human motor behaviour, which represents able-bodied human motor control and is measureable using wearable sensors. The approach is tested in a hardware implementation of the personalisation algorithm on a prosthetic elbow, where the prosthetic objective function is a function of upper-body compensation, and is measured using wearable IMUs. Experimental results on able-bodied subjects using a supernumerary prosthetic elbow mounted on an elbow orthosis suggest that it is possible to use a prosthesis objective function which is implicit in human behaviour to achieve collaboration without providing explicit feedback to the human, motivating further studies.


title: Self-driven oscillation in Coulomb blockaded suspended carbon nanotubes

abstract: Suspended carbon nanotubes are known to support self-driven oscillations due to electromechanical feedback under certain conditions, including low temperatures and high mechanical quality factors. Prior reports identified signatures of such oscillations in Kondo or high-bias transport regimes. Here, we observe self-driven oscillations that give rise to significant conduction in normally Coulomb-blockaded low-bias transport. Using a master equation model, the self-driving is shown to result from strongly energy-dependent electron tunneling, and the dependencies of transport features on bias, gate voltage, and temperature are well reproduced.


title: Laboratory investigations of Lunar ice imaging in permanently shadowed regions using reflected starlight

abstract: A proof of concept for a frost detection imager using reflected starlight is presented; the limitations of this technique are explored experimentally. An ice-covered lunar surface is simulated inside a vacuum chamber, which is then illuminated with a lamp containing UV and visible output to simulate the wavelengths of the background starfield. The simulated lunar surface is imaged with a camera utilizing a UV and visible filter pairing. At Lyman-alpha wavelengths, ice has low reflectivity, and on average appears darker than the regolith in the UV image. In visible wavelengths, this behaviour is reversed, with ice appearing brighter than regolith. UV/VIS image ratioing is subsequently performed in order to discern frost from the lunar regolith simulant in order to demonstrate the capability of this technology for locating the presence of ice on the lunar surface. When the two images are ratioed, the signal to noise ratio to distinguish ice from regolith improves by 36%. In cases where the presence of shadows and specular reflection make distinguishing ice from regolith in either a single UV or visible image difficult, ratioing the images makes the distinction clear.


title: Impurity induced topological phase transitions in Cd$_3$As$_2$ and Na$_3$Bi Dirac semimetals

abstract: Using first-principles density functional theory calculations, combined with a topological analysis, we have investigated the electronic properties of $Cd_3As_2$ and $Na_3Bi$ Dirac topological semimetals doped with non-magnetic and magnetic impurities. Our systematic analysis shows that the selective breaking of the inversion, rotational and time-reversal symmetry, controlled by specific choices of the impurity doping, induces phase transitions from the original Dirac semimetal to a variety of topological phases such as, topological insulator, trivial semimetal, non-magnetic and magnetic Weyl semimetal, and Chern insulator. The Dirac semimetal phase can exist only if the rotational symmetry $C_n$ with $n > 2$ is maintained. One particularly interesting phase emerging in doped $Cd_3As_2$ is a coexisting Dirac-Weyl phase, which occurs when only inversion symmetry is broken while time-reversal symmetry and rotational symmetry are both preserved. To further characterize the low-energy excitations of this phase, we have complemented our density functional results with a continuum four-band $k\cdot p$ model, which indeed displays nodal points of both Dirac and Weyl type. The coexisting phase appears as a transition point between two topologically distinct Dirac phases, but may also survive in a small region of parameter space controlled by external strain.


title: The orthogonal Lie algebra of operators: ideals and derivations

abstract: We study in this paper the infinite-dimensional orthogonal Lie algebra $\mathcal{O}_C$ which consists of all bounded linear operators $T$ on a separable, infinite-dimensional, complex Hilbert space $\mathcal{H}$ satisfying $CTC=-T^*$, where $C$ is a conjugation on $\mathcal{H}$. By employing results from the theory of complex symmetric operators and skew-symmetric operators, we determine the Lie ideals of $\mathcal{O}_C$ and their dual spaces. We study derivations of $\mathcal{O}_C$ and determine their spectra. These results complete some results of P. de la Harpe and provide interesting contrasts between $\mathcal{O}_C$ and the algebra $\mathcal{B(H)}$ of all bounded linear operators on $\mathcal{H}$.


title: An autopilot for energy models -- automatic generation of renewable supply curves, hourly capacity factors and hourly synthetic electricity demand for arbitrary world regions

abstract: Energy system models are increasingly being used to explore scenarios with large shares of variable renewables. This requires input data of high spatial and temporal resolution and places a considerable preprocessing burden on the modeling team. Here we present a new code set with an open source license for automatic generation of input data for large-scale energy system models for arbitrary regions of the world, including sub-national regions, along with an associated generic capacity expansion model of the electricity system. We use ECMWF ERA5 global reanalysis data along with other public geospatial datasets to generate detailed supply curves and hourly capacity factors for solar photovoltaic power, concentrated solar power, onshore and offshore wind power, and existing and future hydropower. Further, we use a machine learning approach to generate synthetic hourly electricity demand series that describe current demand, which we extend to future years using regional SSP scenarios. Finally, our code set automatically generates costs and losses for HVDC interconnections between neighboring regions. The usefulness of our approach is demonstrated by several different case studies based on input data generated by our code. We show that our model runs of a future European electricity system with high share of renewables are in line with results from more detailed models, despite our use of global datasets and synthetic demand.


title: MVC-Net: A Convolutional Neural Network Architecture for Manifold-Valued Images With Applications

abstract: Geometric deep learning has attracted significant attention in recent years, in part due to the availability of exotic data types for which traditional neural network architectures are not well suited. Our goal in this paper is to generalize convolutional neural networks (CNN) to the manifold-valued image case which arises commonly in medical imaging and computer vision applications. Explicitly, the input data to the network is an image where each pixel value is a sample from a Riemannian manifold. To achieve this goal, we must generalize the basic building block of traditional CNN architectures, namely, the weighted combinations operation. To this end, we develop a tangent space combination operation which is used to define a convolution operation on manifold-valued images that we call, the Manifold-Valued Convolution (MVC). We prove theoretical properties of the MVC operation, including equivariance to the action of the isometry group admitted by the manifold and characterizing when compositions of MVC layers collapse to a single layer. We present a detailed description of how to use MVC layers to build full, multi-layer neural networks that operate on manifold-valued images, which we call the MVC-net. Further, we empirically demonstrate superior performance of the MVC-nets in medical imaging and computer vision tasks.


title: Measuring the length of stretched wires with a avalanche oscillator

abstract: We measure the length of stretched wire pairs by using them as the delay line pulse shaping element in an avalanche oscillator. The circuitry and method are extremely simple, and insensitive to oscillator supply voltage, to the particular transistor (2N3904) used and to wire tension. Rudimentary tests with simulated broken wires show that these can be detected easily. It remains to be seen whether the technique will work at scale with realistic wire planes, but it should be at least as good as reflecting short pulses.


title: The integration of photonic crystal waveguides with atom arrays in optical tweezers

abstract: Integrating nanophotonics and cold atoms has drawn increasing interest in recent years due to diverse applications in quantum information science and the exploration of quantum many-body physics. For example, dispersion-engineered photonic crystal waveguides (PCWs) permit not only stable trapping and probing of ultracold neutral atoms via interactions with guided-mode light, but also the possibility to explore the physics of strong, photon-mediated interactions between atoms, as well as atom-mediated interactions between photons. While diverse theoretical opportunities involving atoms and photons in 1-D and 2-D nanophotonic lattices have been analyzed, a grand challenge remains the experimental integration of PCWs with ultracold atoms. Here we describe an advanced apparatus that overcomes several significant barriers to current experimental progress with the goal of achieving strong quantum interactions of light and matter by way of single-atom tweezer arrays strongly coupled to photons in 1-D and 2-D PCWs. Principal technical advances relate to efficient free-space coupling of light to and from guided modes of PCWs, silicate bonding of silicon chips within small glass vacuum cells, and deterministic, mechanical delivery of single-atom tweezer arrays to the near fields of photonic crystal waveguides.


title: A remark on a conjecture of Erd\H{o}s and Straus

abstract: The aim of this note is to show that given a positive integer $n \geq 5$, the positive integral solutions of the diophantine equation $4/n = 1/x + 1/y+1/z$ cannot have solution such that $x$ and $y$ are coprime with $xy < \sqrt{z/2}$. The proof uses the continued fraction expansion of $4/n$.


title: Robot Mindreading and the Problem of Trust

abstract: This paper raises three questions regarding the attribution of beliefs, desires, and intentions to robots. The first one is whether humans in fact engage in robot mindreading. If they do, this raises a second question: does robot mindreading foster trust towards robots? Both of these questions are empirical, and I show that the available evidence is insufficient to answer them. Now, if we assume that the answer to both questions is affirmative, a third and more important question arises: should developers and engineers promote robot mindreading in view of their stated goal of enhancing transparency? My worry here is that by attempting to make robots more mind-readable, they are abandoning the project of understanding automatic decision processes. Features that enhance mind-readability are prone to make the factors that determine automatic decisions even more opaque than they already are. And current strategies to eliminate opacity do not enhance mind-readability. The last part of the paper discusses different ways to analyze this apparent trade-off and suggests that a possible solution must adopt tolerable degrees of opacity that depend on pragmatic factors connected to the level of trust required for the intended uses of the robot.


title: Rapidly Adaptable Legged Robots via Evolutionary Meta-Learning

abstract: Learning adaptable policies is crucial for robots to operate autonomously in our complex and quickly changing world. In this work, we present a new meta-learning method that allows robots to quickly adapt to changes in dynamics. In contrast to gradient-based meta-learning algorithms that rely on second-order gradient estimation, we introduce a more noise-tolerant Batch Hill-Climbing adaptation operator and combine it with meta-learning based on evolutionary strategies. Our method significantly improves adaptation to changes in dynamics in high noise settings, which are common in robotics applications. We validate our approach on a quadruped robot that learns to walk while subject to changes in dynamics. We observe that our method significantly outperforms prior gradient-based approaches, enabling the robot to adapt its policy to changes based on less than 3 minutes of real data.


title: A State of Art Review on Wireless Power Transmission Approaches for Implantable Medical Devices

abstract: Wireless power transmission (WPT) is a critical technology that provides a secure alternative mechanism for wireless power and communication with implantable medical devices. WPT approaches for implantable medical devices have been utilized based on applications. For instance, the inductive coupling tactic is mostly employed for transmission of energy to neuro-stimulators, and the ultrasonic method is used for deep-seated implants. This article provides a study concentrating on popular WPT techniques for implantable medical devices (IMDs) including inductive coupling, microwave, ultrasound, and hybrid WPT systems consisting of two approaches combined. Moreover, an overview of the major works is analyzed with a comparison of their major design elements, operating frequency, distance, efficiency, and harvested power.


title: Deep modeling of quasar variability

abstract: Quasars have long been known as intrinsically variable sources, but the physical mechanism underlying the temporal optical/UV variability is still not well understood. We propose a novel nonparametric method for modeling and forecasting the optical variability of quasars utilizing an autoencoder neural network to gain insight into the underlying processes. The autoencoder is trained with ~15,000 decade-long quasar light curves obtained by the Catalina Real-time Transient Survey selected with negligible flux contamination from the host galaxy. The autoencoder's performance in forecasting the temporal flux variation of quasars is superior to that of the damped random walk process. We find a temporal asymmetry in the optical variability and a novel relation - the amplitude of the variability asymmetry decreases as luminosity and/or black hole mass increases - is suggested with the help of autoencoded features. The characteristics of the variability asymmetry are in agreement with those from the self-organized disk instability model, which predicts that the magnitude of the variability asymmetry decreases as the ratio of the diffusion mass to inflow mass in the accretion disk increases.


title: Improving trial generalizability using observational studies

abstract: Complementary features of randomized controlled trials (RCTs) and observational studies (OSs) can be used jointly to estimate the average treatment effect of a target population. We propose a calibration weighting estimator that enforces the covariate balance between the RCT and OS, therefore improving the trial-based estimator's generalizability. Exploiting semiparametric efficiency theory, we propose a doubly robust augmented calibration weighting estimator that achieves the efficiency bound derived under the identification assumptions. A nonparametric sieve method is provided as an alternative to the parametric approach, which enables the robust approximation of the nuisance functions and data-adaptive selection of outcome predictors for calibration. We establish asymptotic results and confirm the finite sample performances of the proposed estimators by simulation experiments and an application on the estimation of the treatment effect of adjuvant chemotherapy for early-stage non-small cell lung patients after surgery.


title: A Multi-Core Fibre Photonic Lantern Based Spectrograph for Raman Spectroscopy

abstract: We report on the development of a compact (volume $\approx$ 100\:cm$^3$), multimode diffraction-limited Raman spectrograph and probe designed to be compact as possible. The spectrograph uses `off the shelf' optics, a custom 3D-printed two-part housing and harnesses a multi-core fibre (MCF) photonic lantern (multimode to few-mode converter), which slices a large 40~\textmu m multimode input into a near-diffraction-limited 6~\textmu m aperture. Our unique design utilises the hexagonal geometry of our MCF, permitting high multimode collection efficiency with near-diffraction-limited performance in a compact design. Our approach does not require a complex reformatter or mask and thus preserves spectral information and throughput when forming the entrance slit of the spectrograph. We demonstrate the technology over the interval 800~nm to 940~nm (200~cm$^{-1}$ to 2000~cm$^{-1}$) with a resolution of 0.3\:nm (4\:cm$^{-1}$), but other spectral regions and resolutions from the UV to the near infrared are also possible. We demonstrate the performance of our system by recording the Raman spectra of several compounds, including the pharmaceuticals paracetamol and ibuprofen.


title: Universal quivers

abstract: We show that for any positive integer $n$, there exists a quiver $Q$ with $O(n^2)$ vertices and $O(n^2)$ edges such that any quiver on $n$ vertices is a full subquiver of a quiver mutation equivalent to $Q$. We generalize this statement to skew-symmetrizable matrices and obtain other related results. In particular, we show that any quiver is a full subquiver of a quiver mutation equivalent to a quiver of a plabic graph.


title: Rectangular Pyramid Partitioning using Integrated Depth Sensors (RAPPIDS): A Fast Planner for Multicopter Navigation

abstract: We present RAPPIDS: a novel collision checking and planning algorithm for multicopters that is capable of quickly finding local collision-free trajectories given a single depth image from an onboard camera. The primary contribution of this work is a new pyramid-based spatial partitioning method that enables rapid collision detection between candidate trajectories and the environment. By leveraging the efficiency of our collision checking method, we shown how a local planning algorithm can be run at high rates on computationally constrained hardware, evaluating thousands of candidate trajectories in milliseconds. The performance of the algorithm is compared to existing collision checking methods in simulation, showing our method to be capable of evaluating orders of magnitude more trajectories per second. Experimental results are presented showing a quadcopter quickly navigating a previously unseen cluttered environment by running the algorithm on an ODROID-XU4 at 30 Hz.


title: An integral method for the calculation of the reduction in interfacial free energy due to interfacial segregation

abstract: A method based on the Gibbs' adsorption isotherm is developed to calculate the decrease in interfacial free energy resulting from solute segregation at an internal interface, built on measured concentration profiles. Utilizing atom-probe tomography (APT), we first measure a concentration profile of the relative interfacial excess of solute atoms across an interface. To accomplish this we utilize a new method based on J. W. Cahn's formalism for the calculation of the Gibbs interfacial excess. We also introduce a method to calculate the decrease in interfacial free energy that is caused by the segregating solute atoms. This method yields a discrete profile of the decrease in interfacial free energies, which takes into account the measured concentration profile and calculated Gibbsian excess profile. We demonstrate that this method can be used for both homo- and hetero-phase interfaces and takes into account the actual distribution of solute atoms across an interface as determined by APT. It is applied to the case of the semiconducting system PbTe-PbS 12 mol.%-Na 1 mol.%, where Na segregation at the PbS/PbTe interface is anticipated to reduce the interfacial free energy of the {100} facets. We also consider the case of the nickel-based Alloy 600, where B and Si segregation are suspected to impede inter-granular stress corrosion cracking (IGSCC) at homo- (GB) and hetero-phase metal carbide (M7C3) interfaces. The concentration profiles associated with internal interfaces are measured by APT using an ultraviolet (wavelength = 355 nm) laser to dissect nanotips on an atom-by-atom and atomic plane-by-plane basis.


title: Iterative Averaging in the Quest for Best Test Error

abstract: We analyse and explain the increased generalisation performance of iterate averaging using a Gaussian process perturbation model between the true and batch risk surface on the high dimensional quadratic. We derive three phenomena \latestEdits{from our theoretical results:} (1) The importance of combining iterate averaging (IA) with large learning rates and regularisation for improved regularisation. (2) Justification for less frequent averaging. (3) That we expect adaptive gradient methods to work equally well, or better, with iterate averaging than their non-adaptive counterparts. Inspired by these results\latestEdits{, together with} empirical investigations of the importance of appropriate regularisation for the solution diversity of the iterates, we propose two adaptive algorithms with iterate averaging. These give significantly better results compared to stochastic gradient descent (SGD), require less tuning and do not require early stopping or validation set monitoring. We showcase the efficacy of our approach on the CIFAR-10/100, ImageNet and Penn Treebank datasets on a variety of modern and classical network architectures.


title: Night-to-Night Variability of Sleep Electroencephalography-Based Brain Age Measurements

abstract: Objective Brain Age Index (BAI), calculated from sleep electroencephalography (EEG), has been proposed as a biomarker of brain health. This study quantifies night-to-night variability of BAI and establishes probability thresholds for inferring underlying brain pathology based on a patient's BAI. Methods 86 patients with multiple nights of consecutive EEG recordings were selected from Epilepsy Monitoring Unit patients whose EEGs reported as being within normal limits. BAI was calculated for each 12-hour segment of patient data using a previously described algorithm, and night-to-night variability in BAI was measured. Results The within-patient night-to-night standard deviation in BAI was 7.5 years. Estimates of BAI derived by averaging over 2, 3, and 4 nights had standard deviations of 4.7, 3.7, and 3.0 years, respectively. Conclusions Averaging BAI over n nights reduces night-to-night variability of BAI by a factor of the square root of n, rendering BAI more suitable as a biomarker of brain health at the individual level. Significance With increasing ease of EEG acquisition including wearable technology, BAI has the potential to track brain health and detect deviations from normal physiologic function. In a clinical setting, BAI could be used to identify patients who should undergo further investigation or monitoring.


title: Hidden Cost of Randomized Smoothing

abstract: The fragility of modern machine learning models has drawn a considerable amount of attention from both academia and the public. While immense interests were in either crafting adversarial attacks as a way to measure the robustness of neural networks or devising worst-case analytical robustness verification with guarantees, few methods could enjoy both scalability and robustness guarantees at the same time. As an alternative to these attempts, randomized smoothing adopts a different prediction rule that enables statistical robustness arguments which easily scale to large networks. However, in this paper, we point out the side effects of current randomized smoothing workflows. Specifically, we articulate and prove two major points: 1) the decision boundaries of smoothed classifiers will shrink, resulting in disparity in class-wise accuracy; 2) applying noise augmentation in the training process does not necessarily resolve the shrinking issue due to the inconsistent learning objectives.


title: Explicitly Trained Spiking Sparsity in Spiking Neural Networks with Backpropagation

abstract: Spiking Neural Networks (SNNs) are being explored for their potential energy efficiency resulting from sparse, event-driven computations. Many recent works have demonstrated effective backpropagation for deep Spiking Neural Networks (SNNs) by approximating gradients over discontinuous neuron spikes or firing events. A beneficial side-effect of these surrogate gradient spiking backpropagation algorithms is that the spikes, which trigger additional computations, may now themselves be directly considered in the gradient calculations. We propose an explicit inclusion of spike counts in the loss function, along with a traditional error loss, causing the backpropagation learning algorithms to optimize weight parameters for both accuracy and spiking sparsity. As supported by existing theory of over-parameterized neural networks, there are many solution states with effectively equivalent accuracy. As such, appropriate weighting of the two loss goals during training in this multi-objective optimization process can yield an improvement in spiking sparsity without a significant loss of accuracy. We additionally explore a simulated annealing-inspired loss weighting technique to increase the weighting for sparsity as training time increases. Our preliminary results on the Cifar-10 dataset show up to 70.1% reduction in spiking activity with iso-accuracy compared to an equivalent SNN trained only for accuracy and up to 73.3% reduction in spiking activity if allowed a trade-off of 1% reduction in classification accuracy.


title: Point-GNN: Graph Neural Network for 3D Object Detection in a Point Cloud

abstract: In this paper, we propose a graph neural network to detect objects from a LiDAR point cloud. Towards this end, we encode the point cloud efficiently in a fixed radius near-neighbors graph. We design a graph neural network, named Point-GNN, to predict the category and shape of the object that each vertex in the graph belongs to. In Point-GNN, we propose an auto-registration mechanism to reduce translation variance, and also design a box merging and scoring operation to combine detections from multiple vertices accurately. Our experiments on the KITTI benchmark show the proposed approach achieves leading accuracy using the point cloud alone and can even surpass fusion-based algorithms. Our results demonstrate the potential of using the graph neural network as a new approach for 3D object detection. The code is available https://github.com/WeijingShi/Point-GNN.


title: Quasi-periodic oscillations from the accretion disk around distorted black holes

abstract: We study the quasi-periodic oscillations from the accretion disk around the distorted Schwarzschild black hole in the framework of the resonant models. We confine ourselves to the case of a quadrupole distortion which can be caused for example by the accreting matter flow in the vicinity of the compact object. For the purpose we examine the linear stability of the circular geodesic orbits in the equatorial plane and derive analytical expressions for the radial and vertical epicyclic frequencies. We investigate their properties in comparison with the isolated Schwarzschild black hole. Due to the influence of the external matter the vertical epicyclic frequency is not always positive anymore, and the stability of the circular orbits is determined by the interplay between both of the frequencies. As a result, the stable circular orbits do not extend to infinity, but are confined to a finite annular region between an inner and an outer marginally stable orbit. In addition, the degeneracy between the vertical epicyclic and the orbital frequency, which is characteristic for the Schwarzschild solution, is broken, and there are regions in the parametric space where the radial epicyclic frequency is larger than the vertical one. All these properties allow for much more diverse types of non-linear resonances to be excited than for the isolated Schwarzschild black hole, which can provide an explanation for the observed 3:2 ratio between the twin-peak frequencies of the quasi-periodic oscillations from the accretion disk.


title: The absolute discriminant of the endomorphism ring of most reductions of a non-CM elliptic curve is close to maximal

abstract: Let $E/\mathbb{Q}$ be a non-CM elliptic curve. Assuming GRH, we prove that, for a set of primes $p$ of density $1$, the absolute discriminant of the $\mathbb{F}_p$-endomorphism ring of the reduction of $E$ modulo $p$ is close to maximal.


title: Massively Parallel Algorithms for Distance Approximation and Spanners

abstract: Over the past decade, there has been increasing interest in distributed/parallel algorithms for processing large-scale graphs. By now, we have quite fast algorithms -- usually sublogarithmic-time and often $poly(\log\log n)$-time, or even faster -- for a number of fundamental graph problems in the massively parallel computation (MPC) model. This model is a widely-adopted theoretical abstraction of MapReduce style settings, where a number of machines communicate in an all-to-all manner to process large-scale data. Contributing to this line of work on MPC graph algorithms, we present $poly(\log k) \in poly(\log\log n)$ round MPC algorithms for computing $O(k^{1+{o(1)}})$-spanners in the strongly sublinear regime of local memory. To the best of our knowledge, these are the first sublogarithmic-time MPC algorithms for spanner construction. As primary applications of our spanners, we get two important implications, as follows: -For the MPC setting, we get an $O(\log^2\log n)$-round algorithm for $O(\log^{1+o(1)} n)$ approximation of all pairs shortest paths (APSP) in the near-linear regime of local memory. To the best of our knowledge, this is the first sublogarithmic-time MPC algorithm for distance approximations. -Our result above also extends to the Congested Clique model of distributed computing, with the same round complexity and approximation guarantee. This gives the first sub-logarithmic algorithm for approximating APSP in weighted graphs in the Congested Clique model.


title: Height Gap Conjectures, $D$-Finiteness, and Weak Dynamical Mordell-Lang

abstract: In previous work, the first author, Ghioca, and the third author introduced a broad dynamical framework giving rise to many classical sequences from number theory and algebraic combinatorics. Specifically, these are sequences of the form $f(\Phi^n(x))$, where $\Phi\colon X\to X$ and $f\colon X\to\mathbb{P}^1$ are rational maps defined over $\overline{\mathbb{Q}}$ and $x\in X(\overline{\mathbb{Q}})$ is a point whose forward orbit avoids the indeterminacy loci of $\Phi$ and $f$. They conjectured that if the sequence is infinite, then $\limsup \frac{h(f(\Phi^n(x)))}{\log n} > 0$. They also made a corresponding conjecture for $\liminf$ and showed that it implies the Dynamical Mordell-Lang Conjecture. In this paper, we prove the $\limsup$ conjecture as well as the $\liminf$ conjecture away from a set of density $0$. As applications, we prove results concerning the growth rate of coefficients of $D$-finite power series as well as the Dynamical Mordell-Lang Conjecture up to a set of density $0$.


title: Augmented Reality on the Large Scene Based on a Markerless Registration Framework

abstract: In this paper, a mobile camera positioning method based on forward and inverse kinematics of robot is proposed, which can realize far point positioning of imaging position and attitude tracking in large scene enhancement. Orbit precision motion through the framework overhead cameras and combining with the ground system of sensor array object such as mobile robot platform of various sensors, realize the good 3 d image registration, solve any artifacts that is mobile robot in the large space position initialization problem, effectively implement the large space no marks augmented reality, human-computer interaction, and information summary. Finally, the feasibility and effectiveness of the method are verified by experiments.


title: Orders of Growth and Generalized Entropy

abstract: We construct the complete set of orders of growth and we define on it the generalized entropy of a dynamical systems. With this object we provide a framework where we can study the separation of orbits of a map beyond the scope of exponential growth. We are going to show that this construction is particularly useful to study families of dynamical systems with vanishing entropy. Moreover, we are going to see that the space of orders of growth in which orbits are separated is wilder than expected. This is going to be achieved with different types of examples.


title: A note on the refined Strichartz estimates and maximal extension operator

abstract: There are two parts for this paper. In the first part, we extend some results in a recent paper by Du, Guth, Li and Zhang to a more general class of phase functions. The main methods are Bourgain-Demeter's $l^2$ decoupling theorem and induction on scales. In the second part, we prove some positive results for the maximal extension operator for hypersurfaces with positive principal curvatures. The main methods are sharp $L^2$ estimates by Du and Zhang, and the bilinear method by Wolff and Tao.


title: Discontinuity gravity modes in hybrid stars: assessing the role of rapid and slow phase conversions

abstract: Discontinuity gravity modes may arise in perturbed quark-hadron hybrid stars when a sharp density jump exists in the stellar interior and are a potential fingerprint to infer the existence of quark matter cores in compact objects. When a hybrid star is perturbed, conversion reactions may occur at the quark-hadron interface and may have a key role in global stellar properties such as the dynamic stability and the quasi-normal mode spectrum. In this work we study the role of the conversion rate at the interface. To this end, we first derive the junction conditions that hold at the sharp interface of a non-radially perturbed hybrid star in the case of slow and rapid conversions. Then, we analyse the discontinuity $g$-mode in both cases. For rapid conversions, the discontinuity $g$-mode has zero frequency because a displaced fluid element near the phase splitting surface adjusts almost immediately its composition to its surroundings and gravity cannot provide a buoyancy force. For slow conversions, a $g$-mode exists and its properties are analysed here using modern hadronic and quark equations of state. Moreover, it has been shown recently that in the case of slow conversions an extended branch of stable hybrid configurations arises for which $\partial M/ \partial \epsilon_c <0$. We show that $g$-modes of the standard branch (that is, the one with $\partial M/ \partial \epsilon_c > 0$) have frequencies and damping times in agreement with previous results in the literature. However, $g$-modes of the extended branch have significantly larger frequencies (in the range $1-2 \, \mathrm{kHz}$) and much shorter damping times (few seconds in some cases). We discuss the detectability of $g$-mode GWs with present and planned GW observatories.


title: A Fixed Point Framework for Recovering Signals from Nonlinear Transformations

abstract: We consider the problem of recovering a signal from nonlinear transformations, under convex constraints modeling a priori information. Standard feasibility and optimization methods are ill-suited to tackle this problem due to the nonlinearities. We show that, in many common applications, the transformation model can be associated with fixed point equations involving firmly nonexpansive operators. In turn, the recovery problem is reduced to a tractable common fixed point formulation, which is solved efficiently by a provably convergent, block-iterative algorithm. Applications to signal and image recovery are demonstrated. Inconsistent problems are also addressed.


title: Adversarial Network Traffic: Towards Evaluating the Robustness of Deep Learning-Based Network Traffic Classification

abstract: Network traffic classification is used in various applications such as network traffic management, policy enforcement, and intrusion detection systems. Although most applications encrypt their network traffic and some of them dynamically change their port numbers, Machine Learning (ML) and especially Deep Learning (DL)-based classifiers have shown impressive performance in network traffic classification. In this paper, we evaluate the robustness of DL-based network traffic classifiers against Adversarial Network Traffic (ANT). ANT causes DL-based network traffic classifiers to predict incorrectly using Universal Adversarial Perturbation (UAP) generating methods. Since there is no need to buffer network traffic before sending ANT, it is generated live. We partition the input space of the DL-based network traffic classification into three categories: packet classification, flow content classification, and flow time series classification. To generate ANT, we propose three new attacks injecting UAP into network traffic. AdvPad attack injects a UAP into the content of packets to evaluate the robustness of packet classifiers. AdvPay attack injects a UAP into the payload of a dummy packet to evaluate the robustness of flow content classifiers. AdvBurst attack injects a specific number of dummy packets with crafted statistical features based on a UAP into a selected burst of a flow to evaluate the robustness of flow time series classifiers. The results indicate injecting a little UAP into network traffic, highly decreases the performance of DL-based network traffic classifiers in all categories.


title: Selectivity considered harmful: evaluating the causal impact of class selectivity in DNNs

abstract: The properties of individual neurons are often analyzed in order to understand the biological and artificial neural networks in which they're embedded. Class selectivity-typically defined as how different a neuron's responses are across different classes of stimuli or data samples-is commonly used for this purpose. However, it remains an open question whether it is necessary and/or sufficient for deep neural networks (DNNs) to learn class selectivity in individual units. We investigated the causal impact of class selectivity on network function by directly regularizing for or against class selectivity. Using this regularizer to reduce class selectivity across units in convolutional neural networks increased test accuracy by over 2% for ResNet18 trained on Tiny ImageNet. For ResNet20 trained on CIFAR10 we could reduce class selectivity by a factor of 2.5 with no impact on test accuracy, and reduce it nearly to zero with only a small ($\sim$2%) drop in test accuracy. In contrast, regularizing to increase class selectivity significantly decreased test accuracy across all models and datasets. These results indicate that class selectivity in individual units is neither sufficient nor strictly necessary, and can even impair DNN performance. They also encourage caution when focusing on the properties of single units as representative of the mechanisms by which DNNs function.


title: Approximation of noisy data using multivariate splines and finite element methods

abstract: We compare a recently proposed multivariate spline based on mixed partial derivatives with two other standard splines for the scattered data smoothing problem. The splines are defined as the minimiser of a penalised least squares functional. The penalties are based on partial differentiation operators, and are integrated using the finite element method. We compare three methods to two problems: to remove the mixture of Gaussian and impulsive noise from an image, and to recover a continuous function from a set of noisy observations.


title: Quantal-classical fluctuation relation and the second law of thermodynamics: The quantum linear oscillator

abstract: In this work, we study the fluctuation relation and the second law of thermodynamics within a quantum linear oscillator externally driven over the period of time t = tau. To go beyond the standard approach (the two-point projective measurement one) to this subject and also render it discussed in both quantum and classical domains on the single footing, we recast this standard approach in terms of the Wigner function and its propagator in the phase space (x,p). With the help of the canonical transformation from (x,p) to the angle-action coordinates (\phi,I), we can then derive a measurement-free (classical-like) form of the Crooks fluctuation relation in the Wigner representation. This enables us to introduce the work W_{I_0,I_{tau}} associated with a single run from (I_0) to (I_{tau}) over the period tau, which is a quantum generalization of the thermodynamic work with its roots in the classical thermodynamics. This quantum work differs from the energy difference e_{I_0,I_{tau}} = e(I_{tau}) - e(I_0) unless beta, hbar --> 0. Consequently, we will obtain the quantum second-law inequality Delta F_{beta} \leq <W>_{P} \leq <e>_{P} = Delta U, where P, Delta F_{beta}, and <W>_P denote the work (quasi)-probability distribution, the free energy difference, and the average work distinguished from the internal energy difference Delta U, respectively, while <W>_P --> Delta U in the limit of beta, hbar --> 0 only. Therefore, we can also introduce the quantum heat Q_q = Delta U - W even for a thermally isolated system, resulting from the quantum fluctuation therein. This is a more fine-grained result than <W>_P = Delta U obtained from the standard approach. Owing to the measurement-free nature of the thermodynamic work W_{I_0,I_{tau}}, our result can also apply to the (non-thermal) initial states rho_0 = (1-gamma) rho_{beta} + gamma sigma with sigma \ne rho_{beta}.


title: Towards Global Optimal Control via Koopman Lifts

abstract: This paper introduces a framework for solving time-autonomous nonlinear infinite horizon optimal control problems, under the assumption that all minimizers satisfy Pontryagin's necessary optimality conditions. In detail, we use methods from the field of symplectic geometry to analyze the eigenvalues of a Koopman operator that lifts Pontryagin's differential equation into a suitably defined infinite dimensional symplectic space. This has the advantage that methods from the field of spectral analysis can be used to characterize globally optimal control laws. A numerical method for constructing optimal feedback laws for nonlinear systems is then obtained by computing the eigenvalues and eigenvectors of a matrix that is obtained by projecting the Pontryagin-Koopman operator onto a finite dimensional space. We illustrate the effectiveness of this approach by computing accurate approximations of the optimal nonlinear feedback law for a Van der Pol control system, which cannot be stabilized by a linear control law.


title: Chemical freeze-out parameters of net-kaons in heavy-ion collisions

abstract: We study chemical freeze-out parameters for heavy-ion collisions by performing two different thermal analyses. We analyze results from thermal fits for particle yields, as well as, net-charge fluctuations in order to characterize the chemical freeze-out. The Hadron Resonance Gas (HRG) model is employed for both methods. By separating the light hadrons from the strange hadrons in thermal fits, we study the proposed flavor hierarchy. For the net-charge fluctuations, we calculate the mean-over-variance ratio of the net-kaon fluctuations in the HRG model at the five highest energies of the RHIC Beam Energy Scan (BES) for different particle data lists. We compare these results with recent experimental data from the STAR collaboration in order to extract sets of chemical freeze-out parameters for each list. We focused on particle lists which differ largely in the number of resonant states. By doing so, our analysis determines the effect of the amount of resonances included in the HRG model on the freeze-out conditions. Our findings have potential impact on various other models in the field of relativistic heavy-ion collisions.


title: Single-Shot Pose Estimation of Surgical Robot Instruments' Shafts from Monocular Endoscopic Images

abstract: Surgical robots are used to perform minimally invasive surgery and alleviate much of the burden imposed on surgeons. Our group has developed a surgical robot to aid in the removal of tumors at the base of the skull via access through the nostrils. To avoid injuring the patients, a collision-avoidance algorithm that depends on having an accurate model for the poses of the instruments' shafts is used. Given that the model's parameters can change over time owing to interactions between instruments and other disturbances, the online estimation of the poses of the instrument's shaft is essential. In this work, we propose a new method to estimate the pose of the surgical instruments' shafts using a monocular endoscope. Our method is based on the use of an automatically annotated training dataset and an improved pose-estimation deep-learning architecture. In preliminary experiments, we show that our method can surpass state of the art vision-based marker-less pose estimation techniques (providing an error decrease of 55% in position estimation, 64% in pitch, and 69% in yaw) by using artificial images.


title: Impact of temporal correlations on high risk outbreaks of independent and cooperative SIR dynamics

abstract: We first propose a quantitative approach to detect high risk outbreaks of independent and coinfective SIR dynamics on three empirical networks: a school, a conference and a hospital contact network. This measurement is based on the k-means clustering method and identifies proper samples for calculating the mean outbreak size and the outbreak probability. Then we systematically study the impact of different temporal correlations on high risk outbreaks over the original and differently shuffled counterparts of each network. We observe that, on the one hand, in the coinfection process, randomization of the sequence of the events increases the mean outbreak size of high risk cases. On the other hand, these correlations don't have a consistent effect on the independent infection dynamics, and can either decrease or increase this mean. While randomization of the daily pattern correlations has no significant effect on the size of outbreak in either of the coinfection or independent spreading cases. We also observer that an increase in the mean outbreak size doesn't always coincide with an increase in the outbreak probability; therefore we argue that merely considering the mean outbreak size of all realizations may lead us into misestimating the outbreak risks. Our results suggest that some sort of randomizing contacts in organization level of schools, events or hospitals might help to suppress the spreading dynamics while the risk of an outbreak is high.


title: Families of vector fields with many numerical invariants

abstract: We study bifurcations in finite-parameter families of vector fields on $S^2$. Recent papers by Yu. Ilyashenko, N. Goncharuk, Yu. Kudryashov, I. Schurov, and N. Solodovnikov provide examples of (locally generic) structurally unstable families of 3-parameter vector fields: generic close 3-parameter families experience different bifurcations. In this paper, we use these results to construct new examples of few-parameter generic families of planar vector fields such that their classification has many invariants. In particular, we construct (a) 3-parameter families with infinitely many numerical invariants; (b) 4-parameter families with arbitrarily many "robust" numerical invariants; (c) 5-parameter families with functional invariants.


title: Influence Of Climate Change On The Corn Yield In Ontario And Its Impact On Corn Farms Income At The 2068 Horizon

abstract: Our study aims at quantifying the impact of climate change on corn farming in Ontario under several warming scenarios at the 2068 horizon. It is articulated around a discrete-time dynamic model of corn farm income with an annual time-step, corresponding to one agricultural cycle from planting to harvest. At each period, we compute the income given the corn yield, which is highly dependent on weather variables. We also provide a reproducible forecast of the yearly distribution of corn yield for 10 cities in Ontario. The price of corn futures at harvest time is taken into account and we fit our model by using 49 years of historical data. We then conduct out-of-sample Monte-Carlo simulations to obtain the farm income forecasts under a given climate change scenario.