title: Combinatorial Algorithms for Multidimensional Necklaces

abstract: A necklace is an equivalence class of words of length $n$ over an alphabet under the cyclic shift (rotation) operation. As a classical object, there have been many algorithmic results for key operations on necklaces, including counting, generating, ranking, and unranking. This paper generalises the concept of necklaces to the multidimensional setting. We define multidimensional necklaces as an equivalence classes over multidimensional words under the multidimensional cyclic shift operation. Alongside this definition, we generalise several problems from the one dimensional setting to the multidimensional setting for multidimensional necklaces with size $(n_1,n_2,...,n_d)$ over an alphabet of size $q$ including: providing closed form equations for counting the number of necklaces; an $O(n_1 \cdot n_2 \cdot ... \cdot n_d)$ time algorithm for transforming some necklace $w$ to the next necklace in the ordering; an $O((n_1 \cdot n_2 \cdot ... \cdot n_d)^5)$ time algorithm to rank necklaces (determine the number of necklaces smaller than $w$ in the set of necklaces); an $O((n_1\cdot n_2 \cdot ... \cdot n_d)^{6(d + 1)} \cdot \log^d(q))$ time algorithm to unrank multidimensional necklace (determine the $i^{th}$ necklace in the set of necklaces). Our results on counting, ranking, and unranking are further extended to the fixed content setting, where every necklace has the same Parikh vector, in other words every necklace shares the same number of occurrences of each symbol. Finally, we study the $k$-centre problem for necklaces both in the single and multidimensional settings. We provide strong approximation algorithms for solving this problem in both the one dimensional and multidimensional settings.


title: Lung Sound Classification Using Co-tuning and Stochastic Normalization

abstract: In this paper, we use pre-trained ResNet models as backbone architectures for classification of adventitious lung sounds and respiratory diseases. The knowledge of the pre-trained model is transferred by using vanilla fine-tuning, co-tuning, stochastic normalization and the combination of the co-tuning and stochastic normalization techniques. Furthermore, data augmentation in both time domain and time-frequency domain is used to account for the class imbalance of the ICBHI and our multi-channel lung sound dataset. Additionally, we apply spectrum correction to consider the variations of the recording device properties on the ICBHI dataset. Empirically, our proposed systems mostly outperform all state-of-the-art lung sound classification systems for the adventitious lung sounds and respiratory diseases of both datasets.


title: Spatial Search on Johnson Graphs by Continuous-Time Quantum Walk

abstract: Spatial search on graphs is one of the most important algorithmic applications of quantum walks. To show that a quantum-walk-based search is more efficient than a random-walk-based search is a difficult problem, which has been addressed in several ways. Usually, graph symmetries aid in the calculation of the algorithm's computational complexity, and Johnson graphs are an interesting class regarding symmetries because they are regular, Hamilton-connected, vertex- and distance-transitive. In this work, we show that spatial search on Johnson graphs by continuous-time quantum walk achieves the Grover lower bound $\pi\sqrt{N}/2$ with success probability $1$ asymptotically for every fixed diameter, where $N$ is the number of vertices. The proof is mathematically rigorous and can be used for other graph classes.


title: Deformation and breakup dynamics of droplets within a tapered channel

abstract: In this paper we numerically investigate the breakup dynamics of droplets in an emulsion flowing in a tapered microchannel with a narrow constriction. The mesoscale approach for multicomponent fluids with near contact interactions is shown to capture the deformation and breakup dynamics of droplets interacting within the constriction, in agreement with experimental evidences. In addition, it permits to investigate in detail the hydrodynamic phenomena occurring during the breakup stages. Finally, a suitable deformation parameter is introduced and analyzed to characterize the state of deformation of the system by inspecting pairs of interacting droplets flowing in the narrow channel.


title: Staged trees and asymmetry-labeled DAGs

abstract: Bayesian networks are a widely-used class of probabilistic graphical models capable of representing symmetric conditional independence between variables of interest using the topology of the underlying graph. For categorical variables, they can be seen as a special case of the much more general class of models called staged trees, which can represent any type of non-symmetric conditional independence. Here we formalize the relationship between these two models and introduce a minimal Bayesian network representation of the staged tree, which can be used to read conditional independences in an intutitive way. A new labeled graph termed asymmetry-labeled directed acyclic graph is defined, whose edges are labeled to denote the type of dependence existing between any two random variables. We also present a novel algorithm to learn staged trees which only enforces a specific subset of non-symmetric independences. Various datasets are used to illustrate the methodology, highlighting the need to construct models which more flexibly encode and represent non-symmetric structures.


title: Robustness of convolutional neural networks to physiological ECG noise

abstract: The electrocardiogram (ECG) is one of the most widespread diagnostic tools in healthcare and supports the diagnosis of cardiovascular disorders. Deep learning methods are a successful and popular technique to detect indications of disorders from an ECG signal. However, there are open questions around the robustness of these methods to various factors, including physiological ECG noise. In this study we generate clean and noisy versions of an ECG dataset before applying Symmetric Projection Attractor Reconstruction (SPAR) and scalogram image transformations. A pretrained convolutional neural network is trained using transfer learning to classify these image transforms. For the clean ECG dataset, F1 scores for SPAR attractor and scalogram transforms were 0.70 and 0.79, respectively, and the scores decreased by less than 0.05 for the noisy ECG datasets. Notably, when the network trained on clean data was used to classify the noisy datasets, performance decreases of up to 0.18 in F1 scores were seen. However, when the network trained on the noisy data was used to classify the clean dataset, the performance decrease was less than 0.05. We conclude that physiological ECG noise impacts classification using deep learning methods and careful consideration should be given to the inclusion of noisy ECG signals in the training data when developing supervised networks for ECG classification.


title: Exact and Heuristic Approaches to Drone Delivery Problems

abstract: The Flying Sidekick Traveling Salesman Problem (FSTSP) considers a delivery system composed by a truck and a drone. The drone launches from the truck with a single package to deliver to a customer. Each drone must return to the truck to recharge batteries, pick up another package, and launch again to a new customer location. This work proposes a novel Mixed Integer Programming (MIP) formulation and a heuristic approach to address the problem. The proposedMIP formulation yields better linear relaxation bounds than previously proposed formulations for all instances, and was capable of optimally solving several unsolved instances from the literature. A hybrid heuristic based on the General Variable Neighborhood Search metaheuristic combining Tabu Search concepts is employed to obtain high-quality solutions for large-size instances. The efficiency of the algorithm was evaluated on 1415 benchmark instances from the literature, and over 80% of the best known solutions were improved.


title: DuCN: Dual-children Network for Medical Diagnosis and Similar Case Recommendation towards COVID-19

abstract: Early detection of the coronavirus disease 2019 (COVID-19) helps to treat patients timely and increase the cure rate, thus further suppressing the spread of the disease. In this study, we propose a novel deep learning based detection and similar case recommendation network to help control the epidemic. Our proposed network contains two stages: the first one is a lung region segmentation step and is used to exclude irrelevant factors, and the second is a detection and recommendation stage. Under this framework, in the second stage, we develop a dual-children network (DuCN) based on a pre-trained ResNet-18 to simultaneously realize the disease diagnosis and similar case recommendation. Besides, we employ triplet loss and intrapulmonary distance maps to assist the detection, which helps incorporate tiny differences between two images and is conducive to improving the diagnostic accuracy. For each confirmed COVID-19 case, we give similar cases to provide radiologists with diagnosis and treatment references. We conduct experiments on a large publicly available dataset (CC-CCII) and compare the proposed model with state-of-the-art COVID-19 detection methods. The results show that our proposed model achieves a promising clinical performance.


title: Adversarial Energy Disaggregation for Non-intrusive Load Monitoring

abstract: Energy disaggregation, also known as non-intrusive load monitoring (NILM), challenges the problem of separating the whole-home electricity usage into appliance-specific individual consumptions, which is a typical application of data analysis. {NILM aims to help households understand how the energy is used and consequently tell them how to effectively manage the energy, thus allowing energy efficiency which is considered as one of the twin pillars of sustainable energy policy (i.e., energy efficiency and renewable energy).} Although NILM is unidentifiable, it is widely believed that the NILM problem can be addressed by data science. Most of the existing approaches address the energy disaggregation problem by conventional techniques such as sparse coding, non-negative matrix factorization, and hidden Markov model. Recent advances reveal that deep neural networks (DNNs) can get favorable performance for NILM since DNNs can inherently learn the discriminative signatures of the different appliances. In this paper, we propose a novel method named adversarial energy disaggregation (AED) based on DNNs. We introduce the idea of adversarial learning into NILM, which is new for the energy disaggregation task. Our method trains a generator and multiple discriminators via an adversarial fashion. The proposed method not only learns shard representations for different appliances, but captures the specific multimode structures of each appliance. Extensive experiments on real-world datasets verify that our method can achieve new state-of-the-art performance.


title: On simulation of rough Volterra stochastic volatility models

abstract: Rough Volterra volatility models are a progressive and promising field of research in derivative pricing. Although rough fractional stochastic volatility models already proved to be superior in real market data fitting, techniques used in simulation of these models are still inefficient in terms of speed and accuracy. This paper aims to present accurate and efficient tools and techniques for Monte-Carlo simulations for a wide range of rough volatility models. In particular, we compare three commonly used simulation methods: the Cholesky method, the Hybrid scheme, and the rDonsker scheme. We also comment on the implementation of variance reduction techniques. In particular, we show the obstacles of the so-called turbocharging technique whose performance is sometimes counter-productive. To overcome these obstacles, we suggest several modifications.


title: Do What You Know: Coupling Knowledge with Action in Discrete-Event Systems

abstract: An epistemic model for decentralized discrete-event systems with non-binary control is presented. This framework combines existing work on conditional control decisions with existing work on formal reasoning about knowledge in discrete-event systems. The novelty in the model presented is that the necessary and sufficient conditions for problem solvability encapsulate the actions that supervisors must take. This direct coupling between knowledge and action -- in a formalism that mimics natural language -- makes it easier, when the problem conditions fail, to determine how the problem requirements should be revised.


title: Deep Neural Network Approach to Estimate Early Worst-Case Execution Time

abstract: Estimating Worst-Case Execution Time (WCET) is of utmost importance for developing Cyber-Physical and Safety-Critical Systems. The system's scheduler uses the estimated WCET to schedule each task of these systems, and failure may lead to catastrophic events. It is thus imperative to build provably reliable systems. WCET is available to us in the last stage of systems development when the hardware is available and the application code is compiled on it. Different methodologies measure the WCET, but none of them give early insights on WCET, which is crucial for system development. If the system designers overestimate WCET in the early stage, then it would lead to the overqualified system, which will increase the cost of the final product, and if they underestimate WCET in the early stage, then it would lead to financial loss as the system would not perform as expected. This paper estimates early WCET using Deep Neural Networks as an approximate predictor model for hardware architecture and compiler. This model predicts the WCET based on the source code without compiling and running on the hardware architecture. Our WCET prediction model is created using the Pytorch framework. The resulting WCET is too erroneous to be used as an upper bound on the WCET. However, getting these results in the early stages of system development is an essential prerequisite for the system's dimensioning and configuration of the hardware setup.


title: Online unsupervised Learning for domain shift in COVID-19 CT scan datasets

abstract: Neural networks often require large amounts of expert annotated data to train. When changes are made in the process of medical imaging, trained networks may not perform as well, and obtaining large amounts of expert annotations for each change in the imaging process can be time consuming and expensive. Online unsupervised learning is a method that has been proposed to deal with situations where there is a domain shift in incoming data, and a lack of annotations. The aim of this study is to see whether online unsupervised learning can help COVID-19 CT scan classification models adjust to slight domain shifts, when there are no annotations available for the new data. A total of six experiments are performed using three test datasets with differing amounts of domain shift. These experiments compare the performance of the online unsupervised learning strategy to a baseline, as well as comparing how the strategy performs on different domain shifts. Code for online unsupervised learning can be found at this link: https://github.com/Mewtwo/online-unsupervised-learning


title: Robust direct acoustic impedance control using two microphones for mixed feedforward-feedback controller

abstract: This paper presents an acoustic impedance control architecture for an electroacoustic absorber combining both a feedforward and a feedback microphone-based strategies on a current-driven loudspeaker. Feedforward systems enable good performance for direct impedance control. However, inaccuracies in the required actuator model can lead to a loss of passivity, which can cause unstable behaviors. The feedback contribution allows the absorber to better handle model errors and still achieve an accurate impedance, preserving passivity. Numerical and experimental studies were conducted to compare this new architecture against a state-of-the-art feedforward control method.


title: The maximum of the complementary of a semigroup with restricted conditions

abstract: We consider the set $$\mathcal{A} = \left\{10\cdot a + 11\cdot b \ | \gcd(a,b)=1, a\geq 1, b\geq 2a+1 \right\}.$$ We will prove that $\mathcal{A}$ is unbounded and that there exists a natural number $M\notin \mathcal{A}$ for which $$\left\{M+m:m\geq 1,m\in\mathbb N\right\}\subset \mathcal{A}.$$ Indeed, such number is $M = 1674$.


title: A purely data-driven framework for prediction, optimization, and control of networked processes: application to networked SIS epidemic model

abstract: Networks are landmarks of many complex phenomena where interweaving interactions between different agents transform simple local rule-sets into nonlinear emergent behaviors. While some recent studies unveil associations between the network structure and the underlying dynamical process, identifying stochastic nonlinear dynamical processes continues to be an outstanding problem. Here we develop a simple data-driven framework based on operator-theoretic techniques to identify and control stochastic nonlinear dynamics taking place over large-scale networks. The proposed approach requires no prior knowledge of the network structure and identifies the underlying dynamics solely using a collection of two-step snapshots of the states. This data-driven system identification is achieved by using the Koopman operator to find a low dimensional representation of the dynamical patterns that evolve linearly. Further, we use the global linear Koopman model to solve critical control problems by applying to model predictive control (MPC)--typically, a challenging proposition when applied to large networks. We show that our proposed approach tackles this by converting the original nonlinear programming into a more tractable optimization problem that is both convex and with far fewer variables.


title: On the Importance of Domain-specific Explanations in AI-based Cybersecurity Systems (Technical Report)

abstract: With the availability of large datasets and ever-increasing computing power, there has been a growing use of data-driven artificial intelligence systems, which have shown their potential for successful application in diverse areas. However, many of these systems are not able to provide information about the rationale behind their decisions to their users. Lack of understanding of such decisions can be a major drawback, especially in critical domains such as those related to cybersecurity. In light of this problem, in this paper we make three contributions: (i) proposal and discussion of desiderata for the explanation of outputs generated by AI-based cybersecurity systems; (ii) a comparative analysis of approaches in the literature on Explainable Artificial Intelligence (XAI) under the lens of both our desiderata and further dimensions that are typically used for examining XAI approaches; and (iii) a general architecture that can serve as a roadmap for guiding research efforts towards the development of explainable AI-based cybersecurity systems -- at its core, this roadmap proposes combinations of several research lines in a novel way towards tackling the unique challenges that arise in this context.


title: Estimation of road traffic state at a multi-lanes controlled junction

abstract: We present in this paper a method for the estimation of traffic state at road junctions controlled with traffic lights. We assume mixed traffic where a proportion of vehicles are equipped with communication resources. The estimation of road traffic state uses information given by communicating vehicles. The method we propose is built upon a previously published method which was applied to estimate the traffic in the case where roads are composed of two lanes. In this paper, we consider the case where roads are composed of three lanes and we show that this solution can address the general case, where roads are composed of any number of lanes. We assume the geometry of the road junction is known, as well as its connections between incoming and outgoing lanes and roads. Using the location data provided by the communicating vehicles, first, we estimate some primary parameters including the penetration ratio of the probe vehicles, as well as the arrival rates of vehicles (equipped and non-equipped) per lane by introducing the assignment onto the lanes. Second, we give estimations of the queue length of the 3-lanes road, without and with the additional information provided by the location of the communicating vehicles in the queue. We illustrate and discuss the proposed model with numerical simulations.


title: Personal Devices for Contact Tracing: Smartphones and Wearables to Fight Covid-19

abstract: Digital contact tracing has emerged as a viable tool supplementing manual contact tracing. To date, more than 100 contact tracing applications have been published to slow down the spread of highly contagious Covid-19. Despite subtle variabilities among these applications, all of them achieve contact tracing by manipulating the following three components: a) use a personal device to identify the user while designing a secure protocol to anonymize the user's identity; b) leverage networking technologies to analyze and store the data; c) exploit rich sensing features on the user device to detect the interaction among users and thus estimate the exposure risk. This paper reviews the current digital contact tracing based on these three components. We focus on two personal devices that are intimate to the user: smartphones and wearables. We discuss the centralized and decentralized networking approaches that use to facilitate the data flow. Lastly, we investigate the sensing feature available on smartphones and wearables to detect the proximity between any two users and present experiments comparing the proximity sensing performance between these two personal devices.


title: Isolation Intervals of the Real Roots of the Parametric Cubic Equation and Improved Complete Root Classification

abstract: The isolation intervals of the real roots of the real symbolic monic cubic polynomial $p(x) = x^3 + a x^2 + b x + c\,\,$ are found in terms of simple functions of the coefficients of the polynomial (such as: $-a$, $-a/3$, $-c/b$, $\pm \sqrt{-b}$, when $b$ is negative), and the roots of some auxiliary quadratic equations whose coefficients are also simple functions of the coefficients of the cubic. All possible cases are presented with clear and very detailed diagrams. It is very easy to identify which of these diagrams is the relevant one for any given cubic equation and to read from it the isolation intervals of the real roots of the equation. A much-improved complete root classification, addressing the signs (together with giving the isolation intervals) of the individual roots, is also presented. No numerical approximations or root finding techniques are used. Instead of considering the discriminant of the cubic, criterion for the existence of a single real root or three real roots is found as conditions on the coefficients of the cubic, resulting from the roots of the auxiliary quadratic equations. It is also shown that, if a cubic equation has three real roots, then these lie in an interval $I$ such that $\sqrt{3}\sqrt{a^2/3 - b} \le I \le 2 \sqrt{a^2/3 - b}$, independent of $c$. A detailed algorithm for applying the method for isolation of the roots of the cubic is also given and it is illustrated through examples, including the full mathematical analysis of the cubic equation associated with the Rayleigh elastic waves and finding the isolation intervals of its real roots.


title: On the Exploitability of Audio Machine Learning Pipelines to Surreptitious Adversarial Examples

abstract: Machine learning (ML) models are known to be vulnerable to adversarial examples. Applications of ML to voice biometrics authentication are no exception. Yet, the implications of audio adversarial examples on these real-world systems remain poorly understood given that most research targets limited defenders who can only listen to the audio samples. Conflating detectability of an attack with human perceptibility, research has focused on methods that aim to produce imperceptible adversarial examples which humans cannot distinguish from the corresponding benign samples. We argue that this perspective is coarse for two reasons: 1. Imperceptibility is impossible to verify; it would require an experimental process that encompasses variations in listener training, equipment, volume, ear sensitivity, types of background noise etc, and 2. It disregards pipeline-based detection clues that realistic defenders leverage. This results in adversarial examples that are ineffective in the presence of knowledgeable defenders. Thus, an adversary only needs an audio sample to be plausible to a human. We thus introduce surreptitious adversarial examples, a new class of attacks that evades both human and pipeline controls. In the white-box setting, we instantiate this class with a joint, multi-stage optimization attack. Using an Amazon Mechanical Turk user study, we show that this attack produces audio samples that are more surreptitious than previous attacks that aim solely for imperceptibility. Lastly we show that surreptitious adversarial examples are challenging to develop in the black-box setting.


title: High-performance Passive Eigen-model-based Detectors of Single Emitter Using Massive MIMO Receivers

abstract: For a passive direction of arrival (DoA) measurement system using massive multiple input multiple output (MIMO), it is mandatory to infer whether the emitter exists or not before performing DOA estimation operation. Inspired by the detection idea from radio detection and ranging (radar), three high-performance detectors are proposed to infer the existence of single passive emitter from the eigen-space of sample covariance matrix of receive signal vector. The test statistic (TS) of the first method is defined as the ratio of maximum eigen-value (Max-EV) to minimum eigen-value (R-MaxEV-MinEV) while that of the second one is defined as the ratio of Max-EV to noise variance (R-MaxEV-NV). The TS of the third method is the mean of maximum eigen-value (EV) and minimum EV(M-MaxEV-MinEV). Their closed-form expressions are presented and the corresponding detection performance is given. Simulation results show that the proposed M-MaxEV-MinEV and R-MaxEV-NV methods can approximately achieve the same detection performance that is better than the traditional generalized likelihood ratio test method with false alarm probability being less than 0.3.


title: Symplectic foliations induced by harmonic forms on 3-manifolds

abstract: This article details a construction of symplectic foliations on 3-dimensional orientable riemannian manifolds from harmonic forms; and how it suggests a topological approach to Poisson's equation and newtonian gravity.


title: Regularized Limit, analytic continuation and finite-part integration

abstract: Finite-part integration is a recent method of evaluating a convergent integral in terms of the finite-parts of divergent integrals deliberately induced from the convergent integral itself [E. A. Galapon, Proc. R. Soc., A 473, 20160567 (2017)]. Within the context of finite-part integration of the Stieltjes transform of functions with logarithmic growths at the origin, the relationship is established between the analytic continuation of the Mellin transform and the finite-part of the resulting divergent integral when the Mellin integral is extended beyond its strip of analyticity. It is settled that the analytic continuation and the finite-part integral coincide at the regular points of the analytic continuation. To establish the connection between the two at the isolated singularities of the analytic continuation, the concept of regularized limit is introduced to replace the usual concept of limit due to Cauchy when the later leads to a division by zero. It is then shown that the regularized limit of the analytic continuation at its isolated singularities equals the finite-part integrals at the singularities themselves. The treatment gives the exact evaluation of the Stieltjes transform in terms of finite-part integrals and yields the dominant asymptotic behavior of the transform for arbitrarily small values of the parameter in the presence of arbitrary logarithmic singularities at the origin.


title: A unified picture of Balance puzzles and Group testing: Some lessons from quantum mechanics for the pandemic

abstract: Balance (Counterfeit coin) puzzles have been part of recreational mathematics for a few decades. A particular type of Counterfeit coin puzzle is known in the literature as the "Beam balance puzzle". An abstract solution to it is provided by Iwama et.al as a modification of the Bernstein-Vazirani algorithm, making use of quantum parallelism and entanglement. Moreover, during this pandemic, group testing has proved to be an efficient algorithm to save time and cost of testing specimens for the presence of infection. In this article, we propose a "Binary Spring Balance" (BSB) puzzle, to facilitate a unified picture of the counterfeit coin problem and the testing for infection problem, as both aim to reduce the number of queries. We then showcase two solutions to the BSB problem, one using bits and other using classical-qubits ('cebits") for querying. Both solutions are demonstrated using circuits. In this pursuit, we develop a modified optical implementation of Bernstein-Vazirani algorithm using only polarizers (no need of beam splitters), which has surprisingly not yet been proposed earlier. Under the pretext of this demonstration we question why we have not yet developed testing mechanisms inspired by Bernstein-Vazirani algorithm for the pandemic, as they solve the problem in single query, they have no issues related to prevalence of infection in the population, nor are they plagued by the issue of dilution of samples due to pooling. The modified implementation of Bernstein-Vazirani algorithm using polarizers can also be a cost-effective demonstration in an undergraduate lab.


title: LMC N132D: A mature supernova remnant with a power-law gamma-ray spectrum extending beyond 8 TeV

abstract: We analyzed 252 hours of High Energy Stereoscopic System (H.E.S.S.) observations towards the supernova remnant (SNR) LMC N132D that were accumulated between December 2004 and March 2016 during a deep survey of the Large Magellanic Cloud, adding 104 hours of observations to the previously published data set to ensure a > 5 sigma detection. To broaden the gamma-ray spectral coverage required for modeling the spectral energy distribution, an analysis of Fermi-LAT Pass 8 data was also included. We unambiguously detect N132D at very high energies (VHE) with a significance of 5.7 sigma. We report the results of a detailed analysis of its spectrum and localization based on the extended H.E.S.S. data set. The joint analysis of the extended H.E.S.S and Fermi-LAT data results in a spectral energy distribution in the energy range from 1.7 GeV to 14.8 TeV, which suggests a high luminosity of N132D at GeV and TeV energies. We set a lower limit on a gamma-ray cutoff energy of 8 TeV with a confidence level of 95%. The new gamma-ray spectrum as well as multiwavelength observations of N132D when compared to physical models suggests a hadronic origin of the VHE gamma-ray emission. SNR N132D is a VHE gamma-ray source that shows a spectrum extending to the VHE domain without a spectral cutoff at a few TeV, unlike the younger oxygen-rich SNR Cassiopeia A. The gamma-ray properties of N132D may be affected by an interaction with a nearby molecular cloud that partially lies inside the 95% confidence region of the source position. [Abridged]


title: OncoNet: Weakly Supervised Siamese Network to automate cancer treatment response assessment between longitudinal FDG PET/CT examinations

abstract: FDG PET/CT imaging is a resource intensive examination critical for managing malignant disease and is particularly important for longitudinal assessment during therapy. Approaches to automate longtudinal analysis present many challenges including lack of available longitudinal datasets, managing complex large multimodal imaging examinations, and need for detailed annotations for traditional supervised machine learning. In this work we develop OncoNet, novel machine learning algorithm that assesses treatment response from a 1,954 pairs of sequential FDG PET/CT exams through weak supervision using the standard uptake values (SUVmax) in associated radiology reports. OncoNet demonstrates an AUROC of 0.86 and 0.84 on internal and external institution test sets respectively for determination of change between scans while also showing strong agreement to clinical scoring systems with a kappa score of 0.8. We also curated a dataset of 1,954 paired FDG PET/CT exams designed for response assessment for the broader machine learning in healthcare research community. Automated assessment of radiographic response from FDG PET/CT with OncoNet could provide clinicians with a valuable tool to rapidly and consistently interpret change over time in longitudinal multi-modal imaging exams.


title: Critical Droplets and sharp asymptotics for Kawasaki dynamics with weakly anisotropic interactions. Extended version

abstract: In this paper we analyze metastability and nucleation in the context of the Kawasaki dynamics for the two-dimensional Ising lattice gas at very low temperature with periodic boundary conditions. Let $\beta>0$ be the inverse temperature and let $\Lambda\subset\Lambda^\beta\subset\mathbb{Z}^2$ be two boxes. We consider the asymptotic regime corresponding to the limit as $\beta\rightarrow\infty$ for finite volume $\Lambda$ and $\lim_{\beta\rightarrow\infty}\frac{1}{\beta}\log|\Lambda^\beta|=\infty$. We study the simplified model, in which particles perform independent random walks on $\Lambda^\beta\setminus\Lambda$ and inside $\Lambda$ particles perform simple exclusion, but when they occupy neighboring sites they feel a binding energy $-U_1<0$ in the horizontal direction and $-U_2<0$ in the vertical one. Thus the Kawasaki dynamics is conservative inside the volume $\Lambda^\beta$. The initial configuration is chosen such that $\Lambda$ is empty and $\rho|\Lambda^\beta|$ particles are distributed randomly over $\Lambda^\beta\setminus\Lambda$. Our results will use a deep analysis of a local model, i.e., particles perform Kawasaki dynamics inside $\Lambda$ and along each bond touching the boundary of $\Lambda$ from the outside to the inside, particles are created with rate $\rho=e^{-\Delta\beta}$, while along each bond from the inside to the outside, particles are annihilated with rate $1$, where $\Delta>0$ is an activity parameter. Thus, in the local model the boundary of $\Lambda$ plays the role of an infinite gas reservoir with density $\rho$. We take $\Delta\in{(U_1,U_1+U_2)}$, so that the empty (respectively full) configuration is a metastable (respectively stable) configuration. We investigate how the transition from empty to full takes place in the local model with particular attention to the critical configurations that asymptotically have to be crossed with probability 1.


title: Backscatter Communication System With Dumb Diffusing Surface

abstract: Ambient backscatter communications have been identified for ultra-low energy wireless communications. Indeed, a tag can send a message to a reader without emitting any wave and without battery, simply by backscattering the waves generated by a source. In the simplest implementation of such a system, the tag sends a binary message by oscillating between two states and the reader detects the bits by comparing the two distinct received powers. In this paper, for the first time, we propose to study an ambient backscatter communication system, in the presence of a diffusing surface, a simple flat panel that diffuses in all directions. We establish the analytical closed form expression of the power contrast in the presence of the surface. We show that the diffusing surface improves the power contrast. Moreover our approach allows us to express the contrast to noise ratio, and therefore to establish the BER performance. Furthermore, we derive the optimum source transmit power for a given target power contrast. This makes it possible to quantify the amount of energy that can be saved at the source side, thanks to the diffusing surface.


title: Photo-Physical Characteristics of Boron Vacancy-Derived Defect Centers in Hexagonal Boron Nitride

abstract: Single photon emitter (SPE) sources are important building blocks for photonics-based quantum technologies. Recently, the highly bright and versatile SPEs from the two-dimensional insulator material hexagonal boron nitride (hBN) have attracted significant research interest. However, due to the variability of emitter species and properties, an exact correlation between the underlying atomistic structures and their photo-physical properties is still lacking. In this work, we study six boron vacancy-derived defect centers in hBN ($V_B^-$, $V_B$+$H$, $V_{B2}$, $V_BC_N$, $V_BC_N^-$, and $V_BC_NC_N$) using advanced first principles techniques, characterizing their quasiparticle defect levels, optical spectra and excitation energies, and magneto-resonance properties. These defects have been chosen because of their relatively low formation energies, and, importantly, because they are amenable to intentional creation under experimental conditions. We establish the correlation between the underlying defect atomic structure and their photo-physical properties, thus facilitating the identification of SPEs that have already been observed in experiments. Our results lead to clear insights into very recent experiments where hBN SPEs can be controlled intentionally. On the other hand, our results also serve as guidelines for the bottom-up design of defect emitter centers in hBN for target applications that require specific defect properties, such as emission in the telecom wavelength, optical addressability and high radiative decay rates. This work thus provides a comprehensive understanding of the photo-physical characteristics of $V_B$-derived defect emitting centers, aiding in their identification and manipulation for tailored applications.


title: First year of energetic particle measurements in the inner heliosphere with Solar Orbiter's Energetic Particle Detector

abstract: Solar Orbiter strives to unveil how the Sun controls and shapes the heliosphere and fills it with energetic particle radiation. To this end, its Energetic Particle Detector (EPD) has now been in operation, providing excellent data, for just over a year. EPD measures suprathermal and energetic particles in the energy range from a few keV up to (near-) relativistic energies (few MeV for electrons and about 500 MeV/nuc for ions). We present an overview of the initial results from the first year of operations and we provide a first assessment of issues and limitations. During this first year of operations of the Solar Orbiter mission, EPD has recorded several particle events at distances between 0.5 and 1 au from the Sun. We present dynamic and time-averaged energy spectra for ions that were measured with a combination of all four EPD sensors, namely: the SupraThermal Electron and Proton sensor (STEP), the Electron Proton Telescope (EPT), the Suprathermal Ion Spectrograph (SIS), and the High-Energy Telescope (HET) as well as the associated energy spectra for electrons measured with STEP and EPT. We illustrate the capabilities of the EPD suite using the 10-11 December 2020 solar particle event. This event showed an enrichment of heavy ions as well as $^3$He, for which we also present dynamic spectra measured with SIS. The high anisotropy of electrons at the onset of the event and its temporal evolution is also shown using data from these sensors. We discuss the ongoing in-flight calibration and a few open instrumental issues using data from the 21 July and the 10-11 December 2020 events and give guidelines and examples for the usage of the EPD data. We explain how spacecraft operations may affect EPD data and we present a list of such time periods in the appendix. A list of the most significant particle enhancements as observed by EPT during this first year is also provided.


title: Probabilistically nilpotent groups of class two

abstract: For $G$ a finite group, let $d_2(G)$ denote the proportion of triples $(x, y, z) \in G^3$ such that $[x, y, z] = 1$. We determine the structure of finite groups $G$ such that $d_2(G)$ is bounded away from zero: if $d_2(G) \geq \epsilon > 0$, $G$ has a class-4 nilpotent normal subgroup $H$ such that $[G : H] $ and $|\gamma_4(H)|$ are both bounded in terms of $\epsilon$. We also show that if $G$ is an infinite group whose commutators have boundedly many conjugates, or indeed if $G$ satisfies a certain more general commutator covering condition, then $G$ is finite-by-class-3-nilpotent-by-finite.


title: The Jacobian of Cyclic Voltage Covers of $K_n$

abstract: This paper proves results about the Jacobians of a certain family of covering graphs, $Y$, of a base graph $X$, that is constructed from an assignment of elements from a group $G$ to the edges of $X$ ($G$ is called the voltage group and $Y$ is called the derived graph). Of particular interest is when the voltage assignment is given by mapping a generator of the cyclic group of order $d$ to a single edge of $X$ (all other edges are assigned the identity), called a single voltage assignment. Both the order and abelian group structure of the Jacobian of single voltage assignment derived graphs are determined when the base graph $X$ is the complete graph on $n$ vertices, for every $n$ and $d$. Using zeta-functions, general product formulas that relate the order of the Jacobian of $Y$ to that of $X$ are developed; these formulas become very simple and explicit in the special case of single voltage covers of $X$.


title: DFSynthesizer: Dataflow-based Synthesis of Spiking Neural Networks to Neuromorphic Hardware

abstract: Spiking Neural Networks (SNN) are an emerging computation model, which uses event-driven activation and bio-inspired learning algorithms. SNN-based machine-learning programs are typically executed on tile- based neuromorphic hardware platforms, where each tile consists of a computation unit called crossbar, which maps neurons and synapses of the program. However, synthesizing such programs on an off-the-shelf neuromorphic hardware is challenging. This is because of the inherent resource and latency limitations of the hardware, which impact both model performance, e.g., accuracy, and hardware performance, e.g., throughput. We propose DFSynthesizer, an end-to-end framework for synthesizing SNN-based machine learning programs to neuromorphic hardware. The proposed framework works in four steps. First, it analyzes a machine-learning program and generates SNN workload using representative data. Second, it partitions the SNN workload and generates clusters that fit on crossbars of the target neuromorphic hardware. Third, it exploits the rich semantics of Synchronous Dataflow Graph (SDFG) to represent a clustered SNN program, allowing for performance analysis in terms of key hardware constraints such as number of crossbars, dimension of each crossbar, buffer space on tiles, and tile communication bandwidth. Finally, it uses a novel scheduling algorithm to execute clusters on crossbars of the hardware, guaranteeing hardware performance. We evaluate DFSynthesizer with 10 commonly used machine-learning programs. Our results demonstrate that DFSynthesizer provides much tighter performance guarantee compared to current mapping approaches.


title: Well-defined Erd\"{o}s-Straus equations and L.C.M

abstract: The Erd\"{o}s-Straus conjecture states that the equation $\frac{4}{n}=\frac{1}{x}+\frac{1}{y}+\frac{1}{z}$ has positive integer solutions $x, y, z$ for every postive integers $n\ge 2$. We generalize the Erd\"{o}s-Straus equation, state several methods for obtaining well-defined equations, and some features of well-defined equations. Using well-defined equations, we explain the inefficiency of the solutions that did not lead to a complete proof of the problem. Finally, we prove that this conjecture cannot be completely proven by the methods expressed in various articles. For this reason, we express conjectures equivalent to the Erd\"{o}s-Straus conjecture and make the concept of conjecture clearer. The well-defined equations of Erd\"{o}s-Straus have a lot to do with the concept of the least common multiple and the greatest common divisor, and in the last section, we will independently express the functions and features based on the subject of the least common multiple.


title: High-Performance Level-1 and Level-2 BLAS

abstract: The introduction of the Basic Linear Algebra Subroutine (BLAS) in the 1970s paved the way for different libraries to solve the same problem with an improved approach and hardware. The new BLAS implementation led to High-Performance Computing (HPC) innovation. All the love went to the level 3 BLAS due to its humongous application in different fields, not bounded by computer science. However, level 1 and level 2 got neglected; we tried to solve the problem by introducing the new algorithm for the Vector-Vector dot product, Vector-Vector outer product and Matrix-Vector product, which improves the performance of these operations in a significant way. We are not introducing any library but algorithms, which improves upon the current state of art algorithms. Also, we rely on the FMA instruction, OpenMP, and the compiler to optimize the code rather than implementing the algorithm in assembly. Therefore, our current implementation is machine oblivious and depends on the compilers ability to optimize the code.


title: Predictions for gamma-rays from clouds associated with supernova remnant PeVatrons

abstract: Interstellar clouds can act as target material for hadronic cosmic rays; gamma-rays produced through inelastic proton-proton collisions and spatially associated with the clouds provide a key indicator of efficient particle acceleration. However, even for PeVatron sources reaching PeV energies, the system of cloud and accelerator must fulfil several conditions in order to produce a detectable gamma-ray flux. In this contribution, we characterise the necessary properties of both cloud and accelerator. Using available Supernova Remnant (SNR) and interstellar cloud catalogues, and assuming particle acceleration to PeV energies in a nearby SNR, we produce a ranked shortlist of the most promising target systems; those for which a detectable gamma-ray flux is predicted. We discuss detection prospects for future facilities including CTA and SWGO; and compare our predictions with known gamma-ray sources, including the Ultra-High-Energy sources recently detected by LHAASO. A range of model scenarios are tested, including variation in the diffusion coefficient and particle spectrum, under which the best candidate clouds in our shortlist are consistently bright. On average, a detectable gamma-ray flux is more likely for more massive clouds; for systems with lower separation distance between the SNR and cloud; and for slightly older SNRs, due to the time required for particles to traverse the separation distance.


title: On singularity of energy measures for symmetric diffusions with full off-diagonal heat kernel estimates II: Some borderline examples

abstract: We present a concrete family of fractals, which we call the (two-dimensional) thin scale irregular Sierpi\'{n}ski gaskets and each of which is equipped with a canonical strongly local regular symmetric Dirichlet form. We prove that any fractal $K$ in this family satisfies the full off-diagonal heat kernel estimates with some space-time scale function $\Psi_{K}$ and the singularity of the associated energy measures with respect to the canonical volume measure (uniform distribution) on $K$, and also that the decay rate of $r^{-2}\Psi_{K}(r)$ to $0$ as $r\downarrow 0$ can be made arbitrarily slow by suitable choices of $K$. These results together support the energy measure singularity dichotomy conjecture [Ann. Probab. 48 (2020), no. 6, 2920--2951, Conjecture 2.15] stating that, if the full off-diagonal heat kernel estimates with space-time scale function $\Psi$ satisfying $\lim_{r\downarrow 0}r^{-2}\Psi(r)=0$ hold for a strongly local regular symmetric Dirichlet space with complete metric, then the associated energy measures are singular with respect to the reference measure of the Dirichlet space.


title: Incorporating Learnt Local and Global Embeddings into Monocular Visual SLAM

abstract: Traditional approaches for Visual Simultaneous Localization and Mapping (VSLAM) rely on low-level vision information for state estimation, such as handcrafted local features or the image gradient. While significant progress has been made through this track, under more challenging configuration for monocular VSLAM, e.g., varying illumination, the performance of state-of-the-art systems generally degrades. As a consequence, robustness and accuracy for monocular VSLAM are still widely concerned. This paper presents a monocular VSLAM system that fully exploits learnt features for better state estimation. The proposed system leverages both learnt local features and global embeddings at different modules of the system: direct camera pose estimation, inter-frame feature association, and loop closure detection. With a probabilistic explanation of keypoint prediction, we formulate the camera pose tracking in a direct manner and parameterize local features with uncertainty taken into account. To alleviate the quantization effect, we adapt the mapping module to generate 3D landmarks better to guarantee the system's robustness. Detecting temporal loop closure via deep global embeddings further improves the robustness and accuracy of the proposed system. The proposed system is extensively evaluated on public datasets (Tsukuba, EuRoC, and KITTI), and compared against the state-of-the-art methods. The competitive performance of camera pose estimation confirms the effectiveness of our method.


title: Signature Verification using Geometrical Features and Artificial Neural Network Classifier

abstract: Signature verification has been one of the major researched areas in the field of computer vision. Many financial and legal organizations use signature verification as access control and authentication. Signature images are not rich in texture; however, they have much vital geometrical information. Through this work, we have proposed a signature verification methodology that is simple yet effective. The technique presented in this paper harnesses the geometrical features of a signature image like center, isolated points, connected components, etc., and with the power of Artificial Neural Network (ANN) classifier, classifies the signature image based on their geometrical features. Publicly available dataset MCYT, BHSig260 (contains the image of two regional languages Bengali and Hindi) has been used in this paper to test the effectiveness of the proposed method. We have received a lower Equal Error Rate (EER) on MCYT 100 dataset and higher accuracy on the BHSig260 dataset.


title: Spin-conserving Boltzmann theory for carriers and excitons in organic semiconductors

abstract: The rise of organic electronics calls for versatile modeling tools. In this context, we develop a semiclassical Boltzmann theory that describes transport and excitonic processes in crystalline organic semiconductors on equal footing. The generation of singlet and triplet excitons out of the ground state, their formation from free electrons and holes, the reverse processes, as well as the fusion and fission of excitons are included. The corresponding scattering integrals respect spin conservation, which requires matrix-valued distribution functions. They also include fermionic and bosonic many-particle effects such as Pauli blocking. We employ a multipole expansion of the distribution functions, where quadrupolar terms turn out to be essential for the triplet excitons. This work provides a basis for the modeling of organic solar cells, in which excitonic processes are crucial for the performance. Moreover, the theory is of general interest for transport and transitions of multiple (quasi-) particle species carrying spin in nonequilibrium systems.


title: Crack and pull-off dynamics of adhesive, viscoelastic solids

abstract: When quickly detaching an elastomer from a counterface, viscoelasticity dramatically increases the perceived adhesion relative to its adiabatic or equilibrium value. Here, we report simulations on the sticking contact between a rigid cylinder and a viscoelastic half space revealing a maximum in the work of adhesion at intermediate pull-off velocities. Maximum tensile forces yet increase monotonically with the pull-off speed and the crack-tip speed in accordance with the Persson-Brener approach. As predicted theoretically, the fracture mode transitions from interfacial crack propagation to quasi-uniform bond breaking with increasing range of adhesion.


title: Multi-Label Gold Asymmetric Loss Correction with Single-Label Regulators

abstract: Multi-label learning is an emerging extension of the multi-class classification where an image contains multiple labels. Not only acquiring a clean and fully labeled dataset in multi-label learning is extremely expensive, but also many of the actual labels are corrupted or missing due to the automated or non-expert annotation techniques. Noisy label data decrease the prediction performance drastically. In this paper, we propose a novel Gold Asymmetric Loss Correction with Single-Label Regulators (GALC-SLR) that operates robust against noisy labels. GALC-SLR estimates the noise confusion matrix using single-label samples, then constructs an asymmetric loss correction via estimated confusion matrix to avoid overfitting to the noisy labels. Empirical results show that our method outperforms the state-of-the-art original asymmetric loss multi-label classifier under all corruption levels, showing mean average precision improvement up to 28.67% on a real world dataset of MS-COCO, yielding a better generalization of the unseen data and increased prediction performance.


title: Reversible $G^k$-Codes with Applications to DNA Codes

abstract: In this paper, we give a matrix construction method for designing DNA codes that come from group matrix rings. We show that with our construction one can obtain reversible $G^k$-codes of length $kn,$ where $k, n \in \mathbb{N},$ over the finite commutative Frobenius ring $R.$ We employ our construction method to obtain many DNA codes over $\mathbb{F}_4$ that satisfy the Hamming distance, reverse, reverse-complement and the fixed GC-content constraints. Moreover, we improve many lower bounds on the sizes of some known DNA codes and we also give new lower bounds on the sizes of some DNA codes of lengths $48, 56, 60, 64$ and $72$ for some fixed values of the Hamming distance $d.$


title: Dyn-ASR: Compact, Multilingual Speech Recognition via Spoken Language and Accent Identification

abstract: Running automatic speech recognition (ASR) on edge devices is non-trivial due to resource constraints, especially in scenarios that require supporting multiple languages. We propose a new approach to enable multilingual speech recognition on edge devices. This approach uses both language identification and accent identification to select one of multiple monolingual ASR models on-the-fly, each fine-tuned for a particular accent. Initial results for both recognition performance and resource usage are promising with our approach using less than 1/12th of the memory consumed by other solutions.


title: Knowledgeable Prompt-tuning: Incorporating Knowledge into Prompt Verbalizer for Text Classification

abstract: Tuning pre-trained language models (PLMs) with task-specific prompts has been a promising approach for text classification. Particularly, previous studies suggest that prompt-tuning has remarkable superiority in the low-data scenario over the generic fine-tuning methods with extra classifiers. The core idea of prompt-tuning is to insert text pieces, i.e., template, to the input and transform a classification problem into a masked language modeling problem, where a crucial step is to construct a projection, i.e., verbalizer, between a label space and a label word space. A verbalizer is usually handcrafted or searched by gradient descent, which may lack coverage and bring considerable bias and high variances to the results. In this work, we focus on incorporating external knowledge into the verbalizer, forming a knowledgeable prompt-tuning (KPT), to improve and stabilize prompt-tuning. Specifically, we expand the label word space of the verbalizer using external knowledge bases (KBs) and refine the expanded label word space with the PLM itself before predicting with the expanded label word space. Extensive experiments on zero and few-shot text classification tasks demonstrate the effectiveness of knowledgeable prompt-tuning.


title: An estimate of approximation of an analytic function of a matrix by a rational function

abstract: Let $A$ be a square complex matrix; $z_1$, ..., $z_{N}\in\mathbb C$ be arbitrary (possibly repetitive) points of interpolation; $f$ be an analytic function defined on a neighborhood of the convex hull of the union of the spectrum $\sigma(A)$ of the matrix $A$ and the points $z_1$, ..., $z_{N}$; and the rational function $r=\frac uv$ (with the degree of the numerator $u$ less than $N$) interpolates $f$ at these points (counted according to their multiplicities). Under these assumptions estimates of the kind $$ \bigl\Vert f(A)-r(A)\bigr\Vert\le \max_{t\in[0,1];\mu\in\text{convex hull}\{z_1,z_{2},\dots,z_{N}\}}\biggl\Vert\Omega(A)[v(A)]^{-1} \frac{\bigl(vf\bigr)^{{(N)}} \bigl((1-t)\mu\mathbf1+tA\bigr)}{N!}\biggr\Vert, $$ where $\Omega(z)=\prod_{k=1}^N(z-z_k)$, are proposed. As an example illustrating the accuracy of such estimates, an approximation of the impulse response of a dynamic system obtained using the reduced-order Arnoldi method is considered, the actual accuracy of the approximation is compared with the estimate based on this paper.


title: The MIT Supercloud Dataset

abstract: Artificial intelligence (AI) and Machine learning (ML) workloads are an increasingly larger share of the compute workloads in traditional High-Performance Computing (HPC) centers and commercial cloud systems. This has led to changes in deployment approaches of HPC clusters and the commercial cloud, as well as a new focus on approaches to optimized resource usage, allocations and deployment of new AI frame- works, and capabilities such as Jupyter notebooks to enable rapid prototyping and deployment. With these changes, there is a need to better understand cluster/datacenter operations with the goal of developing improved scheduling policies, identifying inefficiencies in resource utilization, energy/power consumption, failure prediction, and identifying policy violations. In this paper we introduce the MIT Supercloud Dataset which aims to foster innovative AI/ML approaches to the analysis of large scale HPC and datacenter/cloud operations. We provide detailed monitoring logs from the MIT Supercloud system, which include CPU and GPU usage by jobs, memory usage, file system logs, and physical monitoring data. This paper discusses the details of the dataset, collection methodology, data availability, and discusses potential challenge problems being developed using this data. Datasets and future challenge announcements will be available via https://dcc.mit.edu.


title: Direct detection of quantum non-Gaussian light from a dispersively coupled single atom

abstract: Many applications in quantum communication, sensing and computation need provably quantum non-Gaussian light. Recently such light, witnessed by a negative Wigner function, has been estimated using homodyne tomography from a single atom dispersively coupled to a high-finesse cavity. This opens an investigation of quantum non-Gaussian light for many experiments with atoms and solid-state emitters. However, at their early stage, an atom or emitter in a cavity system with different channels to the environment and additional noise are insufficient to produce negative Wigner functions. Moreover, homodyne detection is frequently challenging for such experiments. We analyse these issues and prove that such cavities can be used to emit quantum non-Gaussian light employing single-photon detection in the Hanbury Brown and Twiss configuration and quantum non-Gaussianity criteria suitable for this measurement. We investigate in detail cases of considerable cavity leakage when the negativity of the Wigner function disappears completely. Advantageously, quantum non-Gaussian light can be still conclusively proven for a large set of the cavity parameters at the cost of overall measurement time, even if noise is present.


title: Random Convolution Kernels with Multi-Scale Decomposition for Preterm EEG Inter-burst Detection

abstract: Linear classifiers with random convolution kernels are computationally efficient methods that need no design or domain knowledge. Unlike deep neural networks, there is no need to hand-craft a network architecture; the kernels are randomly generated and only the linear classifier needs training. A recently proposed method, RandOm Convolutional KErnel Transforms (ROCKETs), has shown high accuracy across a range of time-series data sets. Here we propose a multi-scale version of this method, using both high- and low-frequency components. We apply our methods to inter-burst detection in a cohort of preterm EEG recorded from 36 neonates <30 weeks gestational age. Two features from the convolution of 10,000 random kernels are combined using ridge regression. The proposed multi-scale ROCKET method out-performs the method without scale: median (interquartile range, IQR) Matthews correlation coefficient (MCC) of 0.859 (0.815 to 0.874) for multi-scale versus 0.841 (0.807 to 0.865) without scale, p<0.001. The proposed method lags behind an existing feature-based machine learning method developed with deep domain knowledge, but is fast to train and can quickly set an initial baseline threshold of performance for generic and biomedical time-series classification.


title: Convergence of gradient descent for learning linear neural networks

abstract: We study the convergence properties of gradient descent for training deep linear neural networks, i.e., deep matrix factorizations, by extending a previous analysis for the related gradient flow. We show that under suitable conditions on the step sizes gradient descent converges to a critical point of the loss function, i.e., the square loss in this article. Furthermore, we demonstrate that for almost all initializations gradient descent converges to a global minimum in the case of two layers. In the case of three or more layers we show that gradient descent converges to a global minimum on the manifold matrices of some fixed rank, where the rank cannot be determined a priori.


title: Node Connectivity Augmentation via Iterative Randomized Rounding

abstract: Many network design problems deal with the design of low-cost networks that are resilient to the failure of their elements, such as nodes or links. One such problem is Connectivity Augmentation, where the goal is to cheaply increase the connectivity of a network from a value k to k+1. The most studied setting focuses on edge-connectivity, which reduces to k=2, called Cactus Augmentation. From an approximation perspective, Byrka, Grandoni, and Jabal Ameli (2020) were the first to break the 2-approximation barrier for this problem, by exploiting a connection to the Steiner Tree problem, and by tailoring the analysis of the iterative randomized rounding technique for Steiner Tree to the specific instances arising from this connection. Recently, Nutov (2020) observed that a similar reduction to Steiner Tree holds for a node-connectivity problem called Block-Tree Augmentation, where the goal is to add edges to a given spanning tree so that the resulting graph 2-node-connected. Combining Nutov's result with the algorithm of Byrka, Grandoni, and Jabal Ameli yields a 1.91-approximation for Block-Tree Augmentation, that is the best bound known so far. In this work, we give a 1.892-approximation algorithm for the problem of augmenting the node-connectivity of any graph from 1 to 2. As a corollary, we improve upon the state-of-the-art approximation factor for Block-Tree Augmentation. Our result is obtained by developing a different and simpler analysis of the iterative randomized rounding technique. Our results also imply a 1.892-approximation algorithm for Cactus Augmentation. While this does not beat the best approximation by Cecchetto, Traub, and Zenklusen (2021) known for this problem, our analysis is quite simple compared to previous results in the literature. In addition, our work gives new insights on the iterative randomized rounding method, that might be of independent interest.


title: A generalization of the root function

abstract: We consider the interpretation and the numerical construction of the inverse branches of $n$ factor Blaschke-products on the disk and show that these provide a generalization of the $n$-th root function. The inverse branches can be defined on pairwise disjoint regions, whose union provides the disk. An explicit formula can be given for the $n$ factor Blaschke-products on the torus, which can be used to provide the inverses on the torus. The inverse branches can be thought of as the solutions $z=z_t(r) (0\le r\le 1)$ to the equation $B(z )=re^{it}$, where $B$ denotes an $n$ factor Blaschke-product. We show that starting from a known value $z_t(1)$, any $z_t(r)$ point of the solution trajectory can be reached in finite steps. The appropriate grouping of the trajectories leads to two natural interpretations of the inverse branches (see Figure 2). We introduce an algorithm which can be used to find the points of the trajectories.


title: The next generation: Impact of high-order analytical information on effective one body waveform models for noncircularized, spin-aligned black hole binaries

abstract: We explore the performance of an updated effective-one-body (EOB) model for spin-aligned coalescing black hole binaries designed to deal with any orbital configuration. The model stems from previous work involving the \TEOBResumS{} waveform model, but incorporates recently computed analytical information up to fifth post-Newtonian (PN) order in the EOB potentials. The dynamics is then informed by Numerical Relativity (NR) quasi-circular simulations (incorporating also recently computed 4PN spin-spin and, optionally, 4.5PN spin-orbit terms). The so-constructed model(s) are then compared to various kind of NR simulations, covering either quasi-circular inspirals, eccentric inspirals and scattering configurations. For quasi-circular (534 datasets) and eccentric (28 datasets) inspirals up to coalescence, the EOB/NR unfaithfulness is well below $1\%$ except for a few outliers in the high, positive, spin corner of the parameter space, where however it does not exceed the $3\%$ level. The EOB values of the scattering angle are found to agree ($\lesssim 1\%$) with the NR predictions for most configurations, with the largest disagreement of only $\sim 4\%$ for the most relativistic one. The inclusion of some high-order analytical information in the orbital sector is useful to improve the EOB/NR agreement with respect to previous work, although the use of NR-informed functions is still crucial to accurately describe the strong-field dynamics and waveform.


title: A Comparison of Different Source Code Representation Methods for Vulnerability Prediction in Python

abstract: In the age of big data and machine learning, at a time when the techniques and methods of software development are evolving rapidly, a problem has arisen: programmers can no longer detect all the security flaws and vulnerabilities in their code manually. To overcome this problem, developers can now rely on automatic techniques, like machine learning based prediction models, to detect such issues. An inherent property of such approaches is that they work with numeric vectors (i.e., feature vectors) as inputs. Therefore, one needs to transform the source code into such feature vectors, often referred to as code embedding. A popular approach for code embedding is to adapt natural language processing techniques, like text representation, to automatically derive the necessary features from the source code. However, the suitability and comparison of different text representation techniques for solving Software Engineering (SE) problems is rarely studied systematically. In this paper, we present a comparative study on three popular text representation methods, word2vec, fastText, and BERT applied to the SE task of detecting vulnerabilities in Python code. Using a data mining approach, we collected a large volume of Python source code in both vulnerable and fixed forms that we embedded with word2vec, fastText, and BERT to vectors and used a Long Short-Term Memory network to train on them. Using the same LSTM architecture, we could compare the efficiency of the different embeddings in deriving meaningful feature vectors. Our findings show that all the text representation methods are suitable for code representation in this particular task, but the BERT model is the most promising as it is the least time consuming and the LSTM model based on it achieved the best overall accuracy(93.8%) in predicting Python source code vulnerabilities.


title: Spin-scattering asymmetry at half-metallic ferromagnet/ferromagnet interface

abstract: We study spin-scattering asymmetry at the interface of two ferromagnets (FMs) based on a half-metallic Co$_{2}$Fe$_{0.4}$Mn$_{0.6}$Si (CFMS)/CoFe interface. First-principles ballistic transport calculations based on Landauer formula for (001)-CoFe/CFMS/CoFe indicate strong spin-dependent conductance at the CFMS/CoFe interface, suggesting a large interface spin-scattering asymmetry coefficient ($\gamma$). Fully epitaxial current-perpendicular-to-plane giant magnetoresistance (CPP-GMR) pseudo-spin-valve (PSV) devices involving CoFe/CFMS/Ag/CFMS/CoFe structures exhibit an enhancement in MR output owing to the formation of the CFMS/CoFe interface at room temperature (RT). This is well reproduced qualitatively by a simulation based on a generalized two-current series-resistor model with considering the presence of $\gamma$ at the CFMS/CoFe interface, half-metallicity of CFMS, and combinations of terminated atoms at the interfaces in the CPP-GMR PSV structure. We show direct evidence for a large $\gamma$ at a half-metallic FM/FM interface and its impact on CPP-GMR effect even at RT.


title: Rethinking real numbers as infinite decimals

abstract: We give a~detailed construction of the complete ordered field of real numbers by means of infinite decimal expansions. We prove that in the canonical encoding of decimals neither addition nor multiplication is {\em computable}, but that both operations are {\em weakly computable}; we introduce both kinds of computability in greater generality. We determine which additive and multiplicative shifts (restrictions of addition and multiplication to one variable) are computable, and prove that each of these shifts becomes computable after a~permutation of encoding. We ask if it is the case for the bivariate addition and multiplication.


title: Plasma-induced surface cooling

abstract: Here we show that, despite a massive incident flux of energetic species, plasmas can induce transient cooling of a material surface. Using time-resolved optical thermometry in-situ with this plasma excitation, we reveal the novel underlying physics that drive this `plasma cooling' that is driven by the diverse chemical and energetic species that comprise this fourth state of matter. We show that the photons and massive particles in the plasma impart energy to different chemical species on a surface, leading to local and temporally changing temperatures that lead to both increases and decreases in temperature on the surface of the sample, even though energy is being imparted to the material. This balance comes from the interplay between chemical reactions, momentum transfer, and energy exchange which occur on different time scales over the course of picoseconds to milliseconds. Thus, we show that through energetically exciting a material with a plasma, we can induce cooling, which can lead to revolutionary advances in refrigeration and thermal mitigation with this new process that is not inhibited by the same limitations in the current state-of-the-art systems.


title: Asymptotic analysis of Poisson shot noise processes, and applications

abstract: Poisson shot noise processes are natural generalizations of compound Poisson processes that have been widely applied in insurance, neuroscience, seismology, computer science and epidemiology. In this paper we study sharp deviations, fluctuations and the stable probability approximation of Poisson shot noise processes. Our achievements extend, improve and complement existing results in the literature. We apply the theoretical results to Poisson cluster point processes, including generalized linear Hawkes processes, and risk processes with delayed claims. Many examples are discussed in detail.


title: Anisotropic curvature measures and volume preserving flows

abstract: In the first part of this paper, we develop the theory of anisotropic curvature measures for convex bodies in the Euclidean space. It is proved that any convex body whose boundary anisotropic curvature measure equals a linear combination of other lower order anisotropic curvature measures with nonnegative coefficients is a scaled Wulff shape. This generalizes the classical results by Schneider [Comment. Math. Helv. \textbf{54} (1979), 42--60] and by Kohlmann [Arch. Math. (Basel) \textbf{70} (1998), 250--256] to the anisotropic setting. The main ingredients in the proof are the generalized anisotropic Minkowski formulas and an inequality of Heintze--Karcher type for convex bodies. In the second part, we consider the volume preserving flow of smooth closed convex hypersurfaces in the Euclidean space with speed given by a positive power $\alpha $ of the $k$th anisotropic mean curvature plus a global term chosen to preserve the enclosed volume of the evolving hypersurfaces. We prove that if the initial hypersurface is strictly convex, then the solution of the flow exists for all time and converges to the Wulff shape in the Hausdorff sense. The characterization theorem for Wulff shapes via the anisotropic curvature measures will be used crucially in the proof of the convergence result. Moreover, in the cases $k=1$, $n$ or $\alpha\geq k$, we can further improve the Hausdorff convergence to the smooth and exponential convergence.


title: ICECAP: Information Concentrated Entity-aware Image Captioning

abstract: Most current image captioning systems focus on describing general image content, and lack background knowledge to deeply understand the image, such as exact named entities or concrete events. In this work, we focus on the entity-aware news image captioning task which aims to generate informative captions by leveraging the associated news articles to provide background knowledge about the target image. However, due to the length of news articles, previous works only employ news articles at the coarse article or sentence level, which are not fine-grained enough to refine relevant events and choose named entities accurately. To overcome these limitations, we propose an Information Concentrated Entity-aware news image CAPtioning (ICECAP) model, which progressively concentrates on relevant textual information within the corresponding news article from the sentence level to the word level. Our model first creates coarse concentration on relevant sentences using a cross-modality retrieval model and then generates captions by further concentrating on relevant words within the sentences. Extensive experiments on both BreakingNews and GoodNews datasets demonstrate the effectiveness of our proposed method, which outperforms other state-of-the-arts. The code of ICECAP is publicly available at https://github.com/HAWLYQ/ICECAP.


title: Forced Imbibition in Stratified Porous Media: Fluid Dynamics and Breakthrough Saturation

abstract: Imbibition, the displacement of a nonwetting fluid by a wetting fluid, plays a central role in diverse energy, environmental, and industrial processes. While this process is typically studied in homogeneous porous media with uniform permeabilities, in many cases, the media have multiple parallel strata of different permeabilities. How such stratification impacts the fluid dynamics of imbibition, as well as the fluid saturation after the wetting fluid breaks through to the end of a given medium, is poorly understood. We address this gap in knowledge by developing an analytical model of imbibition in a porous medium with two parallel strata, combined with a pore network model that explicitly describes fluid crossflow between the strata. By numerically solving these models, we examine the fluid dynamics and fluid saturation left after breakthrough. We find that the breakthrough saturation of nonwetting fluid is minimized when the imposed capillary number Ca is tuned to a value Ca$^*$ that depends on both the structure of the medium and the viscosity ratio between the two fluids. Our results thus provide quantitative guidelines for predicting and controlling flow in stratified porous media, with implications for water remediation, oil/gas recovery, and applications requiring moisture management in diverse materials.


title: SIMPT: Process Improvement Using Interactive Simulation of Time-aware Process Trees

abstract: Process mining techniques including process discovery, conformance checking, and process enhancement provide extensive knowledge about processes. Discovering running processes and deviations as well as detecting performance problems and bottlenecks are well-supported by process mining tools. However, all the provided techniques represent the past/current state of the process. The improvement in a process requires insights into the future states of the process w.r.t. the possible actions/changes. In this paper, we present a new tool that enables process owners to extract all the process aspects from their historical event data automatically, change these aspects, and re-run the process automatically using an interface. The combination of process mining and simulation techniques provides new evidence-driven ways to explore "what-if" questions. Therefore, assessing the effects of changes in process improvement is also possible. Our Python-based web-application provides a complete interactive platform to improve the flow of activities, i.e., process tree, along with possible changes in all the derived activity, resource, and process parameters. These parameters are derived directly from an event log without user-background knowledge.


title: Integrable spin chains and cellular automata with medium range interaction

abstract: We study integrable spin chains and quantum and classical cellular automata with interaction range $\ell\ge 3$. This is a family of integrable models for which there was no general theory so far. We develop an algebraic framework for such models, generalizing known methods from nearest neighbor interacting chains. This leads to a new integrability condition for medium range Hamiltonians, which can be used to classify such models. A partial classification is performed in specific cases, including $U(1)$-symmetric three site interacting models, and Hamiltonians that are relevant for interaction-round-a-face models. We find a number of models which appear to be new. As an application we consider quantum brickwork circuits of various types, including those that can accommodate the classical elementary cellular automata on light cone lattices. In this family we find that the so-called Rule150 and Rule105 models are Yang-Baxter integrable with three site interactions. We present integrable quantum deformations of these models, and derive a set of local conserved charges for them. For the famous Rule54 model we find that it does not belong to the family of integrable three site models, but we can not exclude Yang-Baxter integrability with longer interaction ranges.


title: Partial Reuse AMG Setup Cost Amortization Strategy for the Solution of Non-Steady State Problems

abstract: The partial reuse algebraic multigrid (AMG) setup cost amortization strategy is presented for the solution of non-steady state problems. The transfer operators are reused from the previous time steps, and the system matrices and the smoother operators are rebuilt on each of the AMG hierarchy levels. It is shown on the example of modelling a two-fluid dam break scenario that the strategy may decrease the AMG preconditioner setup cost by 40% to 200%. The total compute time is decreased by up to 20%, but the specific outcome depends on the fraction of time that the setup step initially takes.


title: Recovery of Sobolev functions restricted to iid sampling

abstract: We study $L_q$-approximation and integration for functions from the Sobolev space $W^s_p(\Omega)$ and compare optimal randomized (Monte Carlo) algorithms with algorithms that can only use iid sample points, uniformly distributed on the domain. The main result is that we obtain the same optimal rate of convergence if we restrict to iid sampling, a common assumption in learning and uncertainty quantification. The only exception is when $p=q=\infty$, where a logarithmic loss cannot be avoided.


title: On the non-existence of stepped-pressure equilibria far from symmetry

abstract: The Stepped Pressure Equilibrium Code (SPEC) [Hudson et al., Phys. Plasmas 19, 112502 (2012)] has been successful in the construction of equilibria in 3D configurations that contain a mixture of flux surfaces, islands and chaotic magnetic field lines. In this model, the plasma is sliced into sub-volumes separated by ideal interfaces, and in each volume the magnetic field is a Beltrami field. In the cases where the system is far from possessing a continuous symmetry, such as in stellarators, the existence of solutions to a stepped-pressure equilibrium with given constraints, such as a multi-region relaxed MHD minimum energy state, is not guaranteed but is often taken for granted. Using SPEC, we have studied two different scenarios in which a solution fails to exist in a slab with analytic boundary perturbations. We found that with a large boundary perturbation, a certain interface becomes fractal, corresponding to the break up of a Kolmogorov-Arnold-Moser (KAM) surface. Moreover, an interface can only support a maximum pressure jump while a solution of the magnetic field consistent with the force balance condition can be found. An interface closer to break-up can support a smaller pressure jump. We discovered that the pressure jump can push the interface closer to being non-smooth through force balance, thus significantly decreasing the maximum pressure it can support. Our work shows that a convergence study must be performed on a SPEC equilibrium with interfaces close to break-up. These results may also provide insights into the choice of interfaces and have applications in finding out the maximum pressure a machine can support.


title: Manipulating polariton condensates by Rashba-Dresselhaus coupling at room temperature

abstract: The spin-orbit coupling plays an important role in the spin Hall effect and the topological insulators. In addition, the spin-orbit coupled Bose-Einstein condensates show remarkable quantum many-body phase transition. In this work we tune the exciton polariton condensate by virtue of the Rashba-Dresselhaus (RD) spin-orbit coupling in a liquid-crystal filled microcavity where perovskite CsPbBr3 microplates act as the gain material at room temperature. We realize an artificial gauge field on the CsPbBr3 exciton polariton condensate, which splits the condensates with opposite spins in both momentum and real spaces. Our work paves the way to manipulate the exciton polariton condensate with a synthetic gauge field based on the RD spin-orbit coupling at room temperature.


title: Laplace Transforms of Improper Random Variables

abstract: The probabilistic interpretation of Laplace transforms is used to help to describe the Laplace Transform $L(s)$ of improper random variables. In particular, busy periods in queueing models are examined. The value of $L(0)$ is explained in certain special cases.


abstract: For an image with multiple scene texts, different people may be interested in different text information. Current text-aware image captioning models are not able to generate distinctive captions according to various information needs. To explore how to generate personalized text-aware captions, we define a new challenging task, namely Question-controlled Text-aware Image Captioning (Qc-TextCap). With questions as control signals, this task requires models to understand questions, find related scene texts and describe them together with objects fluently in human language. Based on two existing text-aware captioning datasets, we automatically construct two datasets, ControlTextCaps and ControlVizWiz to support the task. We propose a novel Geometry and Question Aware Model (GQAM). GQAM first applies a Geometry-informed Visual Encoder to fuse region-level object features and region-level scene text features with considering spatial relationships. Then, we design a Question-guided Encoder to select the most relevant visual features for each question. Finally, GQAM generates a personalized text-aware caption with a Multimodal Decoder. Our model achieves better captioning performance and question answering ability than carefully designed baselines on both two datasets. With questions as control signals, our model generates more informative and diverse captions than the state-of-the-art text-aware captioning model. Our code and datasets are publicly available at https://github.com/HAWLYQ/Qc-TextCap.


title: CellProfiler Analyst 3.0: Accessible data exploration and machine learning for image analysis

abstract: Image-based experiments can yield many thousands of individual measurements describing each object of interest, such as cells in microscopy screens. CellProfiler Analyst is a free, open-source software package designed for the exploration of quantitative image-derived data and the training of machine learning classifiers with an intuitive user interface. We have now released CellProfiler Analyst 3.0, which in addition to enhanced performance adds support for neural network classifiers, identifying rare object subsets, and direct transfer of objects of interest from visualisation tools into the Classifier tool for use as training data. This release also increases interoperability with the recently released CellProfiler 4, making it easier for users to detect and measure particular classes of objects in their analyses.


title: gammaALPs: An open-source python package for computing photon-axion-like-particle oscillations in astrophysical environments

abstract: Axions and axion-like particles (ALPs) are hypothetical particles that occur in extensions of the Standard Model and are candidates for cold dark matter. They could be detected through their oscillations into photons in the presence of external electromagnetic fields. gammaALPs is an open-source python framework that computes the oscillation probability between photons and axions/ALPs. In addition to solving the photon-ALP equations of motion, gammaALPs includes models for magnetic fields in different astrophysical environments such as jets of active galactic nuclei, intra-cluster and intergalactic media, and the Milky Way. Users are also able to easily incorporate their own custom magnetic-field models. We review the basic functionality and features of gammaALPs, which is heavily based on other open-source scientific packages such as numpy and scipy. Although focused on gamma-ray energies, gammaALPs can be easily extended to arbitrary photon energies.


title: Morphodynamics of barchan-barchan interactions investigated at the grain scale

abstract: Corridors of size-selected crescent-shaped dunes, known as barchans, are commonly found in water, air, and other planetary environments. The growth of barchans results from the interplay between a fluid flow and a granular bed, but their size regulation involves intricate exchanges between different barchans within a field. One size-regulating mechanism is the binary interaction between nearby dunes, when two dunes exchange mass via the near flow field or by direct contact (collision). In a recent Letter (Assis and Franklin, Geophys. Res. Lett., 2020), we identified five different patterns arising from binary interactions of subaqueous barchans, and proposed classification maps. In this paper, we further inquire into binary exchanges by investigating the motion of individual grains while barchans interact with each other. The experiments were conducted in a water channel where the evolution of pairs of barchans in both aligned and off-centered configurations was recorded by conventional and high-speed cameras. Based on image processing, we obtained the morphology of dunes and motion of grains for all interaction patterns. We present the trajectories of individual grains, from which we show the origin and destination of moving grains, and their typical lengths and velocities. We also show that grains from the impacting dune spread with a diffusion-like component over the target barchan, and we propose a diffusion length. Our results provide new insights into the size-regulating mechanisms of barchans and barchanoid forms found on Earth and other planets.


title: Air leakage in seals with application to syringes

abstract: We study the leakage of air in syringes with Teflon coated rubber stopper and glass barrel. The leakrate depends on the interfacial surface roughness, the viscoelastic properties of the rubber and on the elastoplastic properties of the Teflon coating. The measured leakage rates are compared to the predictions of a simple theory for gas flow, which takes into account both the diffusive and ballistic air flow, and the elastoplastic multiscale contact mechanics which determines the probability distribution of interfacial separations. The theory shows that the interfacial air flow (leakage) channels are so narrow that the gas flow is mainly ballistic (the so called Knudsen limit). The implications for container closure integrity is discussed.


title: Analysis of an Incomplete Binary Outcome Dichotomized From an Underlying Continuous Variable in Clinical Trials

abstract: In many clinical trials, outcomes of interest include binary-valued endpoints. It is not uncommon that a binary-valued outcome is dichotomized from a continuous outcome at a threshold of clinical interest. To reach the objective, common approaches include (a) fitting the generalized linear mixed model (GLMM) to the dichotomized longitudinal binary outcome and (b) imputation method (MI): imputing the missing values in the continuous outcome, dichotomizing it into a binary outcome, and then fitting the generalized linear model for the "complete" data. We conducted comprehensive simulation studies to compare the performance of GLMM with MI for estimating risk difference and logarithm of odds ratio between two treatment arms at the end of study. In those simulation studies, we considered a range of multivariate distribution options for the continuous outcome (including a multivariate normal distribution, a multivariate t-distribution, a multivariate log-normal distribution, and the empirical distribution from a real clinical trial data) to evaluate the robustness of the estimators to various data-generating models. Simulation results demonstrate that both methods work well under those considered distribution options, but MI is more efficient with smaller mean squared errors compared to GLMM. We further applied both the GLMM and MI to 29 phase 3 diabetes clinical trials, and found that the MI method generally led to smaller variance estimates compared to GLMM.


title: Galaxy Zoo: 3D -- Crowd-sourced Bar, Spiral and Foreground Star Masks for MaNGA Target Galaxies

abstract: The challenge of consistent identification of internal structure in galaxies - in particular disc galaxy components like spiral arms, bars, and bulges - has hindered our ability to study the physical impact of such structure across large samples. In this paper we present Galaxy Zoo: 3D (GZ: 3D) a crowdsourcing project built on the Zooniverse platform which we used to create spatial pixel (spaxel) maps that identify galaxy centres, foreground stars, galactic bars and spiral arms for 29831 galaxies which were potential targets of the MaNGA survey (Mapping Nearby Galaxies at Apache Point Observatory, part of the fourth phase of the Sloan Digital Sky Surveys or SDSS-IV), including nearly all of the 10,010 galaxies ultimately observed. Our crowd-sourced visual identification of asymmetric, internal structures provides valuable insight on the evolutionary role of non-axisymmetric processes that is otherwise lost when MaNGA data cubes are azimuthally averaged. We present the publicly available GZ:3D catalog alongside validation tests and example use cases. These data may in the future provide a useful training set for automated identification of spiral arm features. As an illustration, we use the spiral masks in a sample of 825 galaxies to measure the enhancement of star formation spatially linked to spiral arms, which we measure to be a factor of three over the background disc, and how this enhancement increases with radius.


title: Fast frequency modulation is encoded according to the listener expectations in the human subcortical auditory pathway

abstract: Expectations aid and bias our perception. In speech, expected words are easier to recognise than unexpected words, particularly in noisy environments, and incorrect expectations can make us misunderstand our conversational partner. Expectations are combined with the output from the sensory pathways to form representations of speech in the cerebral cortex. However, it is unclear whether expectations are propagated further down to subcortical structures to aid the encoding of the basic dynamic constituent of speech: fast frequency-modulation (FM). Fast FM-sweeps are the basic invariant constituent of consonants, and their correct encoding is fundamental for speech recognition. Here we tested the hypothesis that subjective expectations drive the encoding of fast FM-sweeps characteristic of speech in the human subcortical auditory pathway. We used fMRI to measure neural responses in the human auditory midbrain (inferior colliculus) and thalamus (medial geniculate body). Participants listened to sequences of FM-sweeps for which they held different expectations based on the task instructions. We found robust evidence that the responses in auditory midbrain and thalamus encode the difference between the acoustic input and the subjective expectations of the listener. The results indicate that FM-sweeps are already encoded at the level of the human auditory midbrain and that encoding is mainly driven by subjective expectations. We conclude that the subcortical auditory pathway is integrated in the cortical network of predictive speech processing and that expectations are used to optimise the encoding of even the most basic acoustic constituents of speech.


title: Discovering outliers in the Mars Express thermal power consumption patterns

abstract: The Mars Express (MEX) spacecraft has been orbiting Mars since 2004. The operators need to constantly monitor its behavior and handle sporadic deviations (outliers) from the expected patterns of measurements of quantities that the satellite is sending to Earth. In this paper, we analyze the patterns of the electrical power consumption of MEX's thermal subsystem, that maintains the spacecraft's temperature at the desired level. The consumption is not constant, but should be roughly periodic in the short term, with the period that corresponds to one orbit around Mars. By using long short-term memory neural networks, we show that the consumption pattern is more irregular than expected, and successfully detect such irregularities, opening possibility for automatic outlier detection on MEX in the future.


title: Derivative corrections to the Heisenberg-Euler effective action

abstract: We show that the leading derivative corrections to the Heisenberg-Euler effective action can be determined efficiently from the vacuum polarization tensor evaluated in a homogeneous constant background field. After deriving the explicit parameter-integral representation for the leading derivative corrections in generic electromagnetic fields at one loop, we specialize to the cases of magnetic- and electric-like field configurations characterized by the vanishing of one of the secular invariants of the electromagnetic field. In these cases, closed-form results and the associated all-orders weak- and strong-field expansions can be worked out. One immediate application is the leading derivative correction to the renowned Schwinger-formula describing the decay of the quantum vacuum via electron-positron pair production in slowly-varying electric fields.


title: Extending the sample of core-collapse supernovae for searches of axion-like-particle induced gamma-ray bursts with the Fermi LAT

abstract: During a core-collapse supernova (SN), axion-like particles (ALPs) could be produced through the Primakoff process and subsequently convert into gamma rays in the magnetic field of the Milky Way. Using a sample of well studied extragalactic SNe at optical wavelengths, we estimate the time of the core collapse and search for a coincident gamma-ray burst with the Fermi Large Area Telescope (LAT). Under the assumption that at least one SN was contained within the LAT field of view, we exclude photon-ALP couplings within a factor of $\sim$5 of previous limits from SN1987A. With the increasing number of SNe observed with optical surveys, our results demonstrate the potential to probe ALP dark matter with combined optical and gamma-ray observations. We also provide preliminary results for the estimation of explosion times of 15 close-by SNe observed recently with ZTF. Our findings show that the explosion time can be estimated within one day (statistical uncertainty only) making them promising targets for a follow-up LAT analysis.