title: United States FDA drug approvals are persistent and polycyclic: Insights into economic cycles, innovation dynamics, and national policy

abstract: It is challenging to elucidate the effects of changes in external influences (such as economic or policy) on the rate of US drug approvals. Here, a novel approach, termed the Chronological Hurst Exponent (CHE), is proposed, which hypothesizes that changes in the long-range memory latent within the dynamics of time series data may be temporally associated with changes in such influences. Using the monthly number the FDA Center for Drug Evaluation and Research (CDER) approvals from 1939 to 2019 as the data source, it is demonstrated that the CHE has a distinct S-shaped structure demarcated by an 8-year (1939-1947) Stagnation Period, a 27-year (1947-1974) Emergent (time-varying Period, and a 45-year (1974-2019) Saturation Period. Further, dominant periodicities (resolved via wavelet analyses) are identified during the most recent 45-year CHE Saturation Period at 17, 8 and 4 years; thus, US drug approvals have been following a Juglar-Kuznet mid-term cycle with Kitchin-like bursts. As discussed, this work suggests that (1) changes in extrinsic factors (e.g., of economic and/or policy origin ) during the Emergent Period may have led to persistent growth in US drug approvals enjoyed since 1974, (2) the CHE may be a valued method to explore influences on time series data, and (3) innovation-related economic cycles exist (as viewed via the proxy metric of US drug approvals).


title: Dependence of the old star clusters' dynamical clock on the host galaxy gravitational field

abstract: I report outcomes of the analysis of the $A^+$ parameter, which measures the level of radial segregation of blue straggler stars in old star clusters, commonly known as the dynamical clock for the long-term internal dynamical evolution. I used $A^+$ values available in the literature for 48 Milky Way globular clusters. I found that the relationship of $A^+$ and the number of central relaxation times which have elapsed ($N_{relax}$) shows a non negligible dependence on the strength of the host galaxy gravitational potential, in addition to depending on the two-body relaxation mechanism. Indeed, a measured $A^+$ value corresponds to relatively smaller or larger $N_{relax}$ values for star clusters located farther or closer to the galaxy center. From an observational point of view, this finding reveals the possibility of disentangling for the first time the dynamical evolutionary stage due to two-body relaxation and tidal effect, that affect the whole star clusters' body concurrently.


title: Teaching quantum computing with an interactive textbook

abstract: Quantum computing is a technology that promises to offer significant advantages during the coming decades. Though the technology is still in a prototype stage, the last few years have seen many of these prototype devices become accessible to the public. This has been accompanied by the open-source development of the software required to use and test quantum hardware in increasingly sophisticated ways. Such tools provide new education opportunities, not just for quantum computing specifically, but also more broadly for quantum information science and even quantum physics as a whole. In this paper we present a case study of one education resource which aims to take advantage of the opportunities: the open-source online textbook `Learn Quantum Computation using Qiskit'. An overview of the topics covered is given, as well as an explanation of the approach taken for each.


title: Predictive K-means with local models

abstract: Supervised classification can be effective for prediction but sometimes weak on interpretability or explainability (XAI). Clustering, on the other hand, tends to isolate categories or profiles that can be meaningful but there is no guarantee that they are useful for labels prediction. Predictive clustering seeks to obtain the best of the two worlds. Starting from labeled data, it looks for clusters that are as pure as possible with regards to the class labels. One technique consists in tweaking a clustering algorithm so that data points sharing the same label tend to aggregate together. With distance-based algorithms, such as k-means, a solution is to modify the distance used by the algorithm so that it incorporates information about the labels of the data points. In this paper, we propose another method which relies on a change of representation guided by class densities and then carries out clustering in this new representation space. We present two new algorithms using this technique and show on a variety of data sets that they are competitive for prediction performance with pure supervised classifiers while offering interpretability of the clusters discovered.


title: Learning active learning at the crossroads? evaluation and discussion

abstract: Active learning aims to reduce annotation cost by predicting which samples are useful for a human expert to label. Although this field is quite old, several important challenges to using active learning in real-world settings still remain unsolved. In particular, most selection strategies are hand-designed, and it has become clear that there is no best active learning strategy that consistently outperforms all others in all applications. This has motivated research into meta-learning algorithms for "learning how to actively learn". In this paper, we compare this kind of approach with the association of a Random Forest with the margin sampling strategy, reported in recent comparative studies as a very competitive heuristic. To this end, we present the results of a benchmark performed on 20 datasets that compares a strategy learned using a recent meta-learning algorithm with margin sampling. We also present some lessons learned and open future perspectives.


title: From Weakly Supervised Learning to Biquality Learning: an Introduction

abstract: The field of Weakly Supervised Learning (WSL) has recently seen a surge of popularity, with numerous papers addressing different types of "supervision deficiencies". In WSL use cases, a variety of situations exists where the collected "information" is imperfect. The paradigm of WSL attempts to list and cover these problems with associated solutions. In this paper, we review the research progress on WSL with the aim to make it as a brief introduction to this field. We present the three axis of WSL cube and an overview of most of all the elements of their facets. We propose three measurable quantities that acts as coordinates in the previously defined cube namely: Quality, Adaptability and Quantity of information. Thus we suggest that Biquality Learning framework can be defined as a plan of the WSL cube and propose to re-discover previously unrelated patches in WSL literature as a unified Biquality Learning literature.


title: Visualization and Selection of Dynamic Mode Decomposition Components for Unsteady Flow

abstract: Dynamic Mode Decomposition (DMD) is a data-driven and model-free decomposition technique. It is suitable for revealing spatio-temporal features of both numerically and experimentally acquired data. Conceptually, DMD performs a low-dimensional spectral decomposition of the data into the following components: The modes, called DMD modes, encode the spatial contribution of the decomposition, whereas the DMD amplitudes specify their impact. Each associated eigenvalue, referred to as DMD eigenvalue, characterizes the frequency and growth rate of the DMD mode. In this paper, we demonstrate how the components of DMD can be utilized to obtain temporal and spatial information from time-dependent flow fields. We begin with the theoretical background of DMD and its application to unsteady flow. Next, we examine the conventional process with DMD mathematically and put it in relationship to the discrete Fourier transform. Our analysis shows that the current use of DMD components has several drawbacks. To resolve these problems we adjust the components and provide new and meaningful insights into the decomposition: We show that our improved components describe the flow more adequately. Moreover, we remove redundancies in the decomposition and clarify the interplay between components, allowing users to understand the impact of components. These new representations ,which respect the spatio-temporal character of DMD, enable two clustering methods that segment the flow into physically relevant sections and can therefore be used for the selection of DMD components. With a number of typical examples, we demonstrate that the combination of these techniques allow new insights with DMD for unsteady flow.


title: GN-z11-flash in the context of Gamma-Ray Burst Afterglows

abstract: The recently discovered rapid transient GN-z11-flash has been suggested to be the prompt-emission ultraviolet flash associated with a gamma-ray burst serendipitously exploding in the ultra-high-$z$ galaxy GN-z11. We here place the flash into the context of the early ultraviolet emission of gamma-ray bursts, and find it is in agreement with the luminosity distribution of these events.


title: Quantum d-separation and quantum belief propagation

abstract: The goal of this paper is to generalize classical d-separation and classical Belief Propagation (BP) to the quantum realm. Classical d-separation is an essential ingredient of most of Judea Pearl's work. It is crucial to all 3 rungs of what Pearl calls the 3 rungs of Causation. So having a quantum version of d-separation and BP probably implies that most of Pearl's Bayesian networks work, including his theory of causality, can be translated in a straightforward manner to the quantum realm.


title: XAI-P-T: A Brief Review of Explainable Artificial Intelligence from Practice to Theory

abstract: In this work, we report the practical and theoretical aspects of Explainable AI (XAI) identified in some fundamental literature. Although there is a vast body of work on representing the XAI backgrounds, most of the corpuses pinpoint a discrete direction of thoughts. Providing insights into literature in practice and theory concurrently is still a gap in this field. This is important as such connection facilitates a learning process for the early stage XAI researchers and give a bright stand for the experienced XAI scholars. Respectively, we first focus on the categories of black-box explanation and give a practical example. Later, we discuss how theoretically explanation has been grounded in the body of multidisciplinary fields. Finally, some directions of future works are presented.


title: Dynamical tides in neutron stars: The impact of the crust

abstract: We consider the dynamical tidal response of a neutron star in an inspiralling binary, focussing on the impact of the star's elastic crust. Within the context of Newtonian gravity, we add the elastic aspects to the theoretical formulation of the problem and quantify the dynamical excitation of different classes of oscillation modes. The results demonstrate the expectation that the fundamental mode dominates the tidal response and show how the usual tidal deformability (and the Love number) emerge in the static limit. In addition, we consider to what extent the different modes may be excited to a level where the breaking strain of the crust would be exceeded (locally). The results show that the fundamental mode may fracture the crust during the late stages of inspiral. This is also the case for the first gravity mode, which reaches the breaking threshold in strongly stratified stars. In our models with a fluid ocean, interface modes associated with the crust-ocean transition may also induce crust fracture. If this happens it does so earlier in the inspiral, at a lower orbital frequency.


title: Chaos Game Representation

abstract: The chaos game representation (CGR) is an interesting method to visualize one-dimensional sequences. In this paper, we show how to construct a chaos game representation. The applications mentioned here are biological, in which CGR was able to uncover patterns in DNA or proteins that were previously unknown. We also show how CGR might be introduced in the classroom, either in a modelling course or in a dynamical systems course. Some sequences that are tested are taken from the Online Encyclopedia of Integer Sequences, and others are taken from sequences that arose mainly from a course in experimental mathematics.


title: Preface to the 22nd workshop "What comes beyond the standard models", Bled July 06--July 14, 2019, and links to the talks appearing in the proceedings, and/or in arXiv

abstract: The contribution contains the preface to the Proceedings to the 22nd Workshop "What Comes Beyond the Standard Models", Bled, July 06 - July 14, 2019, published in Bled workshops in physics, Vol.20, No. 2, DMFA-Zalo\v{z}nistvo, Ljubljana, Dec. 2019, links to (most of) the published contributions and section (by M.Yu. Khlopov) on VIA at Bled 2019.


title: The achievements of the {\it spin-charge-family} theory so far

abstract: Fifty years ago the {\it standard model} offered an elegant new step towards understanding elementary fermion and boson fields, making several assumptions, suggested by experiments. The assumptions are still waiting for an explanation. There are many proposals in the literature for the next step. The {\it spin-charge-family} theory, proposing a simple starting action in $ d\ge (13+1)$-dimensional space with fermions interacting with the gravity only (the vielbeins and the two kinds of the spin connection fields), is offering the explanation for not only all by the {\it standard model} assumed properties of quarks and leptons and antiquarks and antileptons, with the families included, of the vectors gauge fields, of the Higgs's scalar and Yukawa couplings, of the appearance of the {\it dark matter}, of the {\it matter-antimatter asymmetry}, making several predictions, but explains as well the second quantization postulates for fermions and bosons by using the odd and the even Clifford algebra ''basis vectors'' to describe the internal space of fermions and bosons, respectively. Consequently the single fermion and single boson states already anticommute and commute, respectively. I present in this talk a very short overview of the achievements of the {\it spin-charge-family} theory so far, concluding with presenting not yet solved problems, for which the collaborators are very welcome.


title: Spatial-Temporal Fusion Graph Neural Networks for Traffic Flow Forecasting

abstract: Spatial-temporal data forecasting of traffic flow is a challenging task because of complicated spatial dependencies and dynamical trends of temporal pattern between different roads. Existing frameworks typically utilize given spatial adjacency graph and sophisticated mechanisms for modeling spatial and temporal correlations. However, limited representations of given spatial graph structure with incomplete adjacent connections may restrict effective spatial-temporal dependencies learning of those models. To overcome those limitations, our paper proposes Spatial-Temporal Fusion Graph Neural Networks (STFGNN) for traffic flow forecasting. SFTGNN could effectively learn hidden spatial-temporal dependencies by a novel fusion operation of various spatial and temporal graphs, which is generated by a data-driven method. Meanwhile, by integrating this fusion graph module and a novel gated convolution module into a unified layer, SFTGNN could handle long sequences. Experimental results on several public traffic datasets demonstrate that our method achieves state-of-the-art performance consistently than other baselines.


title: Limit of Weierstrass Measure on Stable Curves

abstract: The goal of the paper is to study the limiting behavior of the Weierstrass measures on a smooth curve of genus $g\geqslant 2$ as the curve approaches a certain nodal stable curve represented by a point in the Deligne-Mumford compactification $\bar{\mathcal M}_g$ of the moduli $\mathcal{M}_g$, including irreducible ones or those of compact type. As a consequence, the Weierstrass measures on a stable rational curve at the boundary of $\mathcal{M}_g$ are completely determined. In the process, the asymptotic behavior of the Bergman measure is also studied.


title: Automatic source localization and spectra generation from sparse beamforming maps

abstract: Beamforming is an imaging tool for the investigation of aeroacoustic phenomena and results in high dimensional data that is broken down to spectra by integrating spatial Regions Of Interest. This paper presents two methods that enable the automated identification of aeroacoustic sources in sparse beamforming maps and the extraction of their corresponding spectra to overcome the manual definition of Regions Of Interest. The methods are evaluated on two scaled airframe half-model wind-tunnel measurements and on a generic monopole source. The first relies on the spatial normal distribution of aeroacoustic broadband sources in sparse beamforming maps. The second uses hierarchical clustering methods. Both methods are robust to statistical noise and predict the existence, location, and spatial probability estimation for sources based on which Regions Of Interest are automatically determined.


title: Splitting of quaternions and octonions over purely inseparable extensions in characteristic 2

abstract: We give examples of quaternion and octonion division algebras over a field $F$ of characteristic $2$ that split over a purely inseparable extension $E$ of $F$ of degree $\geq 4$ but that do not split over any subextension of $F$ inside $E$ of lower exponent, or, in the case of octonions, over any simple subextension of $F$ inside $E$. Thus, we get a negative answer to a question posed by Bernhard M\"uhlherr and Richard Weiss. We study this question in terms of the isotropy behaviour of the associated norm forms.


title: On the Importance of Diversity in Re-Sampling for Imbalanced Data and Rare Events in Mortality Risk Models

abstract: Surgical risk increases significantly when patients present with comorbid conditions. This has resulted in the creation of numerous risk stratification tools with the objective of formulating associated surgical risk to assist both surgeons and patients in decision-making. The Surgical Outcome Risk Tool (SORT) is one of the tools developed to predict mortality risk throughout the entire perioperative period for major elective in-patient surgeries in the UK. In this study, we enhance the original SORT prediction model (UK SORT) by addressing the class imbalance within the dataset. Our proposed method investigates the application of diversity-based selection on top of common re-sampling techniques to enhance the classifier's capability in detecting minority (mortality) events. Diversity amongst training datasets is an essential factor in ensuring re-sampled data keeps an accurate depiction of the minority/majority class region, thereby solving the generalization problem of mainstream sampling approaches. We incorporate the use of the Solow-Polasky measure as a drop-in functionality to evaluate diversity, with the addition of greedy algorithms to identify and discard subsets that share the most similarity. Additionally, through empirical experiments, we prove that the performance of the classifier trained over diversity-based dataset outperforms the original classifier over ten external datasets. Our diversity-based re-sampling method elevates the performance of the UK SORT algorithm by 1.4$.


title: DAG-based Scheduling with Resource Sharing for Multi-task Applications in a Polyglot GPU Runtime

abstract: GPUs are readily available in cloud computing and personal devices, but their use for data processing acceleration has been slowed down by their limited integration with common programming languages such as Python or Java. Moreover, using GPUs to their full capabilities requires expert knowledge of asynchronous programming. In this work, we present a novel GPU run time scheduler for multi-task GPU computations that transparently provides asynchronous execution, space-sharing, and transfer-computation overlap without requiring in advance any information about the program dependency structure. We leverage the GrCUDA polyglot API to integrate our scheduler with multiple high-level languages and provide a platform for fast prototyping and easy GPU acceleration. We validate our work on 6 benchmarks created to evaluate task-parallelism and show an average of 44% speedup against synchronous execution, with no execution time slowdown compared to hand-optimized host code written using the C++ CUDA Graphs API.


title: Ultra-Fast, Low-Storage, Highly Effective Coarse-grained Selection in Retrieval-based Chatbot by Using Deep Semantic Hashing

abstract: We study the coarse-grained selection module in retrieval-based chatbot. Coarse-grained selection is a basic module in a retrieval-based chatbot, which constructs a rough candidate set from the whole database to speed up the interaction with customers. So far, there are two kinds of approaches for coarse-grained selection module: (1) sparse representation; (2) dense representation. To the best of our knowledge, there is no systematic comparison between these two approaches in retrieval-based chatbots, and which kind of method is better in real scenarios is still an open question. In this paper, we first systematically compare these two methods from four aspects: (1) effectiveness; (2) index stoarge; (3) search time cost; (4) human evaluation. Extensive experiment results demonstrate that dense representation method significantly outperforms the sparse representation, but costs more time and storage occupation. In order to overcome these fatal weaknesses of dense representation method, we propose an ultra-fast, low-storage, and highly effective Deep Semantic Hashing Coarse-grained selection method, called DSHC model. Specifically, in our proposed DSHC model, a hashing optimizing module that consists of two autoencoder models is stacked on a trained dense representation model, and three loss functions are designed to optimize it. The hash codes provided by hashing optimizing module effectively preserve the rich semantic and similarity information in dense vectors. Extensive experiment results prove that, our proposed DSHC model can achieve much faster speed and lower storage than sparse representation, with limited performance loss compared with dense representation. Besides, our source codes have been publicly released for future research.


title: Dynamic Reinsurance in Discrete Time Minimizing the Insurer's Cost of Capital

abstract: In the classical static optimal reinsurance problem, the cost of capital for the insurer's risk exposure determined by a monetary risk measure is minimized over the class of reinsurance treaties represented by increasing Lipschitz retained loss functions. In this paper, we consider a dynamic extension of this reinsurance problem in discrete time which can be viewed as a risk-sensitive Markov Decision Process. The model allows for both insurance claims and premium income to be stochastic and operates with general risk measures and premium principles. We derive the Bellman equation and show the existence of a Markovian optimal reinsurance policy. Under an infinite planning horizon, the model is shown to be contractive and the optimal reinsurance policy to be stationary. The results are illustrated with examples where the optimal policy can be determined explicitly.


title: Engineering Three Dimensional Moir\'e Flat Bands

abstract: We demonstrate that the concept of moir\'e flat bands can be generalized to achieve electronic band engineering in all three spatial dimensions. For many two dimensional van der Waals materials, twisting two adjacent layers with respect to each other leads to flat electronic bands in the two corresponding spatial directions -- a notion sometimes referred to as twistronics as it enables a wealth of physical phenomena. Within this two dimensional plane, large moir\'e patterns of nanometer size form. The basic concept we propose here is to stack multiple twisted layers on top of each other in a predefined pattern. If the pattern is chosen such that with respect to the stacking direction of layers, the large spatial moir\'e features are spatially shifted from one twisted layer to the next, the system exhibits twist angle controlled flat bands in all of the three spatial directions. With this, our proposal extends the use of twistronic to three dimensions. We exemplify the general concept by considering graphitic systems, boron nitride and WSe$_2$ as candidate materials, but the approach is applicable to any two-dimensional van der Waals material. For hexagonal boron nitride we develope an ab initio fitted tight binding model that captures the corresponding three dimensional low-energy electronic structure. We outline that interesting three dimensional correlated phases of matter can be induced and controlled following this route, including quantum magnets and unconventional superconducting states.

categories: cond-mat.mtrl-sci cond-mat.stat-mech cond-mat.str-el cond-mat.supr-con


title: A White Box Analysis of ColBERT

abstract: Transformer-based models are nowadays state-of-the-art in ad-hoc Information Retrieval, but their behavior is far from being understood. Recent work has claimed that BERT does not satisfy the classical IR axioms. However, we propose to dissect the matching process of ColBERT, through the analysis of term importance and exact/soft matching patterns. Even if the traditional axioms are not formally verified, our analysis reveals that ColBERT: (i) is able to capture a notion of term importance; (ii) relies on exact matches for important terms.


title: Uniform Convolution and Fourier Restriction estimates for complex polynomial curves in $\mathbb{C}^3$

abstract: We establish optimal $(p,q)$ ranges for two types of estimates associated to three dimensional complex polynomial curves. These are the estimates for the weighted restriction of the Fourier Transform to a complex polynomial curve, and the weighted Convolution Operator associated to a complex polynomial curve. Establishing these estimates comes down to establishing a lower bound for the Jacobian of a mapping associated to the complex curve in question.


title: Constructible sheaves and functions up to infinity

abstract: We introduce the category of b-analytic manifolds, a natural tool to define constructible sheaves and functions up to infinity. We study with some details the operations on these objects and also recall the Radon transform for constructible functions.


title: Bargmann Invariants, Geometric Phases and Recursive Parametrization with Majorana Fermions

abstract: A generalized connection between the quantum mechanical Bargmann invariants and the geometric phases was established for the Dirac fermions. We extend that formalism for the Majorana fermions by defining proper quantum mechanical ray and Hilbert spaces. We then relate both the Dirac and Majorana type Bargmann invariants to the rephasing invariant measures of CP violation with the Majorana neutrinos, assuming that the neutrinos have lepton number violating Majorana masses. We then generalize the recursive parametrization for studying any unitary matrices to include the Majorana fermions, which could be useful for studying the neutrino mixing matrix.


title: Detection and Prediction of Nutrient Deficiency Stress using Longitudinal Aerial Imagery

abstract: Early, precise detection of nutrient deficiency stress (NDS) has key economic as well as environmental impact; precision application of chemicals in place of blanket application reduces operational costs for the growers while reducing the amount of chemicals which may enter the environment unnecessarily. Furthermore, earlier treatment reduces the amount of loss and therefore boosts crop production during a given season. With this in mind, we collect sequences of high-resolution aerial imagery and construct semantic segmentation models to detect and predict NDS across the field. Our work sits at the intersection of agriculture, remote sensing, and modern computer vision and deep learning. First, we establish a baseline for full-field detection of NDS and quantify the impact of pretraining, backbone architecture, input representation, and sampling strategy. We then quantify the amount of information available at different points in the season by building a single-timestamp model based on a UNet. Next, we construct our proposed spatiotemporal architecture, which combines a UNet with a convolutional LSTM layer, to accurately detect regions of the field showing NDS; this approach has an impressive IOU score of 0.53. Finally, we show that this architecture can be trained to predict regions of the field which are expected to show NDS in a later flight -- potentially more than three weeks in the future -- maintaining an IOU score of 0.47-0.51 depending on how far in advance the prediction is made. We will also release a dataset which we believe will benefit the computer vision, remote sensing, as well as agriculture fields. This work contributes to the recent developments in deep learning for remote sensing and agriculture, while addressing a key social challenge with implications for economics and sustainability.


title: Analytical and fast Fiber Orientation Distribution reconstruction in 3D-Polarized Light Imaging

abstract: Three dimensional Polarized Light Imaging (3D-PLI) is an optical technique which allows mapping the spatial fiber architecture of fibrous postmortem tissues, at sub-millimeter resolutions. Here, we propose an analytical and fast approach to compute the fiber orientation distribution (FOD) from high-resolution vector data provided by 3D-PLI. The FOD is modeled as a sum of K orientations/Diracs on the unit sphere, described on a spherical harmonics basis and analytically computed using the spherical Fourier transform. Experiments are performed on rich synthetic data which simulate the geometry of the neuronal fibers and on human brain data. Results indicate the analytical FOD is computationally efficient and very fast, and has high angular precision and angular resolution. Furthermore, investigations on the right occipital lobe illustrate that our strategy of FOD computation enables the bridging of spatial scales from microscopic 3D-PLI information to macro- or mesoscopic dimensions of diffusion Magnetic Resonance Imaging (MRI), while being a means to evaluate prospective resolution limits for diffusion MRI to reconstruct region-specific white matter tracts. These results demonstrate the interest and great potential of our analytical approach.


title: Large Deviation in Continuous Time Random Walks

abstract: We discuss large deviation properties of continuous-time random walks (CTRW) and present a general expression for the large deviation rate in CTRW in terms of the corresponding rates for the distributions of steps' lengths and waiting times. In the case of Gaussian distribution of steps' lengths the general expression reduces to a sequence of two Legendre transformations applied to the cumulant generating function of waiting times. The discussion of several examples (Bernoulli and Gaussian random walks with exponentially distributed waiting times, Gaussian random walks with one-sided L\'evy and Pareto-distributed waiting times) reveals interesting general properties of such large deviations.


title: Formation of Singularities in Plasma Ion Dynamics

abstract: We study the formation of singularity for the Euler-Poisson system equipped with the Boltzmann relation, which describes the dynamics of ions in an electrostatic plasma. In general, it is known that smooth solutions to nonlinear hyperbolic equations fail to exist globally in time. We establish criteria for $C^1$ blow-up of the Euler-Poisson system, both for the isothermal and pressureless cases. In particular, our blow-up condition for the presureless model does not require that the gradient of velocity is negatively large. In fact, our result particularly implies that the smooth solutions can break down even if the gradient of initial velocity is trivial. For the isothermal case, we prove that smooth solutions leave $C^1$ class in a finite time when the gradients of the Riemann functions are initially large.


title: Metabolic response to point mutations reveals principles of modulation of in vivo enzyme activity and phenotype

abstract: The relationship between sequence variation and phenotype is poorly understood. Here we use metabolomic analysis to elucidate the molecular mechanism underlying the filamentous phenotype of E. coli strains that carry destabilizing mutations in the Dihydrofolate Reductase (DHFR). We find that partial loss of DHFR activity causes SOS response indicative of DNA damage and cell filamentation. This phenotype is triggered by an imbalance in deoxy nucleotide levels, most prominently a disproportionate drop in the intracellular dTTP. We show that a highly cooperative (Hill coefficient 2.5) in vivo activity of Thymidylate Kinase (Tmk), a downstream enzyme that phosphorylates dTMP to dTDP, is the cause of suboptimal dTTP levels. dTMP supplementation in the media rescues filamentation and restores in vivo Tmk kinetics to almost perfect Michaelis-Menten, like its kinetics in vitro. Overall, this study highlights the important role of cellular environment in sculpting enzymatic kinetics with system level implications for bacterial phenotype.


title: A convolution type model for the intensity of spatial point processes applied to eye-movement data

abstract: Estimating the first-order intensity function in point pattern analysis is an important problem, and it has been approached so far from different perspectives: parametrically, semiparametrically or nonparametrically. Our approach is close to a semiparametric one. Motivated by eye-movement data, we introduce a convolution type model where the log-intensity is modelled as the convolution of a function $\beta(\cdot)$, to be estimated, and a single spatial covariate (the image an individual is looking at for eye-movement data). Based on a Fourier series expansion, we show that the proposed model \rev{can be viewed as a} log-linear model with an infinite number of coefficients, which correspond to the spectral decomposition of $\beta(\cdot)$. After truncation, we estimate these coefficients through a penalized Poisson likelihood. We illustrate the efficiency of the proposed methodology on simulated data and on eye-movement data.


title: Fast whole-brain imaging of seizures in zebrafish larvae by two-photon light-sheet microscopy

abstract: Light-sheet fluorescence microscopy (LSFM) enables real-time whole-brain functional imaging in zebrafish larvae. Conventional one photon LSFM can however induce undesirable visual stimulation due to the use of visible excitation light. The use of two-photon (2P) excitation, employing near-infrared invisible light, provides unbiased investigation of neuronal circuit dynamics. However, due to the low efficiency of the 2P absorption process, the imaging speed of this technique is typically limited by the signal-to-noise-ratio. Here, we describe a 2P LSFM setup designed for non-invasive imaging that enables quintuplicating state-of-the-art volumetric acquisition rate of the larval zebrafish brain (5 Hz) while keeping low the laser intensity on the specimen. We applied our system to the study of pharmacologically-induced acute seizures, characterizing the spatial-temporal dynamics of pathological activity and describing for the first time the appearance of caudo-rostral ictal waves (CRIWs).


title: Geometric Brownian motion with affine drift and its time-integral

abstract: The joint distribution of a geometric Brownian motion and its time-integral was derived in a seminal paper by Yor (1992) using Lamperti's transformation, leading to explicit solutions in terms of modified Bessel functions. In this paper, we revisit this classic result using the simple Laplace transform approach in connection to the Heun differential equation. We extend the methodology to the geometric Brownian motion with affine drift and show that the joint distribution of this process and its time-integral can be determined by a doubly-confluent Heun equation. Furthermore, the joint Laplace transform of the process and its time-integral is derived from the asymptotics of the solutions. In addition, we provide an application by using the results for the asymptotics of the double-confluent Heun equation in pricing Asian options. Numerical results show the accuracy and efficiency of this new method.

categories: q-fin.MF math.CA math.PR q-fin.CP q-fin.PR


title: Firearm Detection via Convolutional Neural Networks: Comparing a Semantic Segmentation Model Against End-to-End Solutions

abstract: Threat detection of weapons and aggressive behavior from live video can be used for rapid detection and prevention of potentially deadly incidents such as terrorism, general criminal offences, or even domestic violence. One way for achieving this is through the use of artificial intelligence and, in particular, machine learning for image analysis. In this paper we conduct a comparison between a traditional monolithic end-to-end deep learning model and a previously proposed model based on an ensemble of simpler neural networks detecting fire-weapons via semantic segmentation. We evaluated both models from different points of view, including accuracy, computational and data complexity, flexibility and reliability. Our results show that a semantic segmentation model provides considerable amount of flexibility and resilience in the low data environment compared to classical deep model models, although its configuration and tuning presents a challenge in achieving the same levels of accuracy as an end-to-end model.


title: Architecture aware compilation of quantum circuits via lazy synthesis

abstract: Qubit routing is a key problematic related to quantum circuit compilation. It consists in rewriting a quantum circuit by adding the least possible number of instructions to make the circuit compliant with some architecture's connectivity constraints. Usually, this problem is tackled via either SWAP insertion techniques or re-synthesis of portions of the circuit using architecture aware synthesis algorithms. In this work, we propose a meta-heuristic that couples the iterative approach of SWAP insertion techniques with greedy architecture aware synthesis routines. We propose two new compilation algorithms based on this meta-heuristic and compare their performances to state-of-the-art quantum circuit compilation techniques for several standard classes of quantum circuits and show significant reduction in the entangling gate overhead due to compilation.


title: Multiperipheral final states in crowded twin-jet events at the LHC

abstract: The 13 TeV run of the LHC has provided a unique opportunity to explore multi-jet final states with unprecedented accuracy. An interesting region for study is that of events where one jet is tagged in the forward direction and another one in the backward direction and a plethora of low energy mini-jets populate the possibly large rapidity span in between them. Since the number of these events is very high, it is possible to introduce stringent constraints on the transverse momentum of the two leading jets which can be kept in small windows not very different from each other, defining what we call "twin jets". The associated "crowd" of mini-jets can also have a restricted span in transverse momentum. The study of these events for a fixed multiplicity is an ideal playground to investigate different models of multi-particle production in hadron-hadron collisions. We set up an exploratory analysis by using an ancient model of Chew and Pignotti to describe the gross features one can expect for the structure of single and double differential-in-rapidity cross sections and for particle-particle rapidity correlations when the longitudinal phase space completely decouples from the transverse degrees of freedom.


title: Large response of charge stripes to uniaxial stress in $\textrm{La}_{1.475}\textrm{Nd}_{0.4}\textrm{Sr}_{0.125}\textrm{Cu}\textrm{O}_{4}$

abstract: The La-based '214' cuprates host several symmetry breaking phases including superconductivity, charge and spin order in the form of stripes, and a structural othorhombic-to-tetragonal phase transition. Therefore, these materials are an ideal system to study the effects of uniaxial stress onto the various correlations that pervade the cuprate phase diagram. We report resonant x-ray scattering experiments on $\textrm{La}_{1.475}\textrm{Nd}_{0.4}\textrm{Sr}_{0.125}\textrm{Cu}\textrm{O}_{4}$ (LNSCO-125) that reveal a significant response of charge stripes to uniaxial tensile-stress of $\sim$ 0.1 GPa. These effects include a reduction of the onset temperature of stripes by $\sim$ 50 K, a 29 K reduction of the low-temperature orthorhombic-to-tetragonal transition, competition between charge order and superconductivity, and a preference for stripes to form along the direction of applied stress. Altogether, we observe a dramatic response of the electronic properties of LNSCO-125 to a modest amount of uniaxial stress.


title: A fully pipelined FPGA accelerator for scale invariant feature transform keypoint descriptor matching,

abstract: The scale invariant feature transform (SIFT) algorithm is considered a classical feature extraction algorithm within the field of computer vision. SIFT keypoint descriptor matching is a computationally intensive process due to the amount of data consumed. In this work, we designed a novel fully pipelined hardware accelerator architecture for SIFT keypoint descriptor matching. The accelerator core was implemented and tested on a field programmable gate array (FPGA). The proposed hardware architecture is able to properly handle the memory bandwidth necessary for a fully-pipelined implementation and hits the roofline performance model, achieving the potential maximum throughput. The fully pipelined matching architecture was designed based on the consine angle distance method. Our architecture was optimized for 16-bit fixed-point operations and implemented on hardware using a Xilinx Zynq-based FPGA development board. Our proposed architecture shows a noticeable reduction of area resources compared with its counterparts in literature, while maintaining high throughput by alleviating memory bandwidth restrictions. The results show a reduction in consumed device resources of up to 91 percent in LUTs and 79 percent of BRAMs. Our hardware implementation is 15.7 times faster than the comparable software approach.


title: Multi-Modal Depth Estimation Using Convolutional Neural Networks

abstract: This paper addresses the problem of dense depth predictions from sparse distance sensor data and a single camera image on challenging weather conditions. This work explores the significance of different sensor modalities such as camera, Radar, and Lidar for estimating depth by applying Deep Learning approaches. Although Lidar has higher depth-sensing abilities than Radar and has been integrated with camera images in lots of previous works, depth estimation using CNN's on the fusion of robust Radar distance data and camera images has not been explored much. In this work, a deep regression network is proposed utilizing a transfer learning approach consisting of an encoder where a high performing pre-trained model has been used to initialize it for extracting dense features and a decoder for upsampling and predicting desired depth. The results are demonstrated on Nuscenes, KITTI, and a Synthetic dataset which was created using the CARLA simulator. Also, top-view zoom-camera images captured from the crane on a construction site are evaluated to estimate the distance of the crane boom carrying heavy loads from the ground to show the usability in safety-critical applications.


title: Field-Tunable Topological Phase Transitions and Spin-Hall Effects in 2D Crystals

abstract: As recent additions to the catalog of 2D crystals, silicene and other silicene-class crystals have numerous unique properties currently being investigated and considered for use in novel device applications. In this paper, we investigate electronic and transport properties of silicene in a field effect transistor geometry. We find that the Berry curvature of silicene-class crystals can be continuously tuned by a perpendicular electric field. By direct calculation of the $Z_2$ invariant, we confirm that an electronic phase transition from a topological insulator to a band insulator occurs when the electric field passes a critical value. In a device setting with asymmetric gate voltages, this field-tunable Berry curvature generates a large spin current transverse to the charge current. When the electric field strength surpasses the critical value, the bulk spin current is found to change direction and greatly decrease in magnitude. This finding of a large magnitude, switchable spin current suggests that the silicene family of 2D crystals could be an attractive candidate for field-tunable charge-spin conversion. Such field-tunable phase transitions between topologically distinct phases could be useful for robust qubits as well.


title: New bi-harmonic superspace formulation of $4D, \mathcal{N}=4$ SYM theory

abstract: We develop a novel bi-harmonic $\mathcal{N}=4$ superspace formulation of the $\mathcal{N}=4$ supersymmetric Yang-Mills theory (SYM) in four dimensions. In this approach, the $\mathcal{N}=4$ SYM superfield constraints are solved in terms of on-shell $\mathcal {N}=2$ harmonic superfields. Such an approach provides a convenient tool of constructing the manifestly $\mathcal{N}=4$ supersymmetric invariants and further rewriting them in $\mathcal{N}= 2$ harmonic superspace. In particular, we present $\mathcal{N}=4$ superfield form of the leading term in the $\mathcal{N}=4$ SYM effective action which was known previously in $\mathcal{N}=2$ superspace formulation.


title: RainBench: Towards Global Precipitation Forecasting from Satellite Imagery

abstract: Extreme precipitation events, such as violent rainfall and hail storms, routinely ravage economies and livelihoods around the developing world. Climate change further aggravates this issue. Data-driven deep learning approaches could widen the access to accurate multi-day forecasts, to mitigate against such events. However, there is currently no benchmark dataset dedicated to the study of global precipitation forecasts. In this paper, we introduce \textbf{RainBench}, a new multi-modal benchmark dataset for data-driven precipitation forecasting. It includes simulated satellite data, a selection of relevant meteorological data from the ERA5 reanalysis product, and IMERG precipitation data. We also release \textbf{PyRain}, a library to process large precipitation datasets efficiently. We present an extensive analysis of our novel dataset and establish baseline results for two benchmark medium-range precipitation forecasting tasks. Finally, we discuss existing data-driven weather forecasting methodologies and suggest future research avenues.


title: Kerr-like nonlinearities in an optomechanical system with an asymmetric anharmonic mechanical resonator

abstract: In the framework of the nonsecular perturbation theory based on the Bogoliubov averaging method, an optomechanical system with an asymmetric anharmonic mechanical resonator is studied. The cross-Kerr interaction and the Kerr-like self-interaction of photons and vibration quanta arise in the Hamiltonian. These interactions are induced by both cubic and quartic nonlinearities of oscillations of the mechanical resonator and the cavity-resonator interaction that is linear in mechanical displacements. We demonstrate a bistable behavior of the number of vibration quanta and find that this behavior is controlled by the cross-Kerr interaction. It is shown that, without driving and dissipation, the constructed superposition Yurke-Stoler-like states of the cavity (or the mechanical resonator) disentangle at certain times the entangled modes of the system. The obtained results offer new possibilities for control of optomechanical systems with asymmetric mechanical oscillations.


title: Unitarity in KK-graviton production: A case study in warped extra-dimensions

abstract: The Kaluza-Klein (KK) decomposition of higher-dimensional gravity gives rise to a tower of KK-gravitons in the effective four-dimensional (4D) theory. Such massive spin-2 fields are known to be connected with unitarity issues and easily lead to a breakdown of the effective theory well below the naive scale of the interaction. However, the breakdown of the effective 4D theory is expected to be controlled by the parameters of the 5D theory. Working in a simplified Randall-Sundrum model we study the matrix elements for matter annihilations into massive gravitons. We find that truncating the KK-tower leads to an early breakdown of perturbative unitarity. However, by considering the full tower we obtain a set of sum rules for the couplings between the different KK-fields that restore unitarity up to the scale of the 5D theory. We prove analytically that these are fulfilled in the model under consideration and present numerical tests of their convergence. This work complements earlier studies that focused on graviton self-interactions and yields additional sum rules that are required if matter fields are incorporated into warped extra-dimensions.


title: Combating Mode Collapse in GAN training: An Empirical Analysis using Hessian Eigenvalues

abstract: Generative adversarial networks (GANs) provide state-of-the-art results in image generation. However, despite being so powerful, they still remain very challenging to train. This is in particular caused by their highly non-convex optimization space leading to a number of instabilities. Among them, mode collapse stands out as one of the most daunting ones. This undesirable event occurs when the model can only fit a few modes of the data distribution, while ignoring the majority of them. In this work, we combat mode collapse using second-order gradient information. To do so, we analyse the loss surface through its Hessian eigenvalues, and show that mode collapse is related to the convergence towards sharp minima. In particular, we observe how the eigenvalues of the $G$ are directly correlated with the occurrence of mode collapse. Finally, motivated by these findings, we design a new optimization algorithm called nudged-Adam (NuGAN) that uses spectral information to overcome mode collapse, leading to empirically more stable convergence properties.


title: Harish-Chandra Cuspidal Pairs

abstract: The irreducible characters of a finite reductive group are partitioned into Harish-Chandra series that are labelled by cuspidal pairs. In this note, we describe how one can algorithmically calculate those cuspidal pairs using results of Lusztig.


title: On a problem of Chen and Fang related to infinite additive complements

abstract: Two infinite sets $A$ and $B$ of nonnegative integers are called additive complements if their sumset contains every nonnegative integer. In 1964, Danzer constructed infinite additive complements $A$ and $B$ with $A(x)B(x) = (1 + o(1))x$ as $x \rightarrow \infty$, where $A(x)$ and $B(x)$ denote the counting function of the sets $A$ and $B$, respectively. In this paper we solve a problem of Chen and Fang by extending the construction of Danzer.


title: Transition between Instability and Seeded Self-Modulation of a Relativistic Particle Bunch in Plasma

abstract: We use a relativistic ionization front to provide various initial transverse wakefield amplitudes for the self-modulation of a long proton bunch in plasma. We show experimentally that, with sufficient initial amplitude ($\ge(4.1\pm0.4)$ MV/m), the phase of the modulation along the bunch is reproducible from event to event, with 3 to 7% (of 2$\pi$) rms variations all along the bunch. The phase is not reproducible for lower initial amplitudes. We observe the transition between these two regimes. Phase reproducibility is essential for deterministic external injection of particles to be accelerated.


title: Floquet Gauge Pumps as Sensors for Spectral Degeneracies Protected by Symmetry or Topology

abstract: We introduce the concept of a Floquet gauge pump whereby a dynamically engineered Floquet Hamiltonian is employed to reveal the inherent degeneracy of the ground state in interacting systems. We demonstrate this concept in a one-dimensional XY model with periodically driven couplings and transverse field. In the high-frequency limit, we obtain the Floquet Hamiltonian consisting of the static XY and dynamically generated Dzyaloshinsky-Moriya interaction (DMI) terms. The dynamically generated magnetization current depends on the phases of complex coupling terms, with the XY interaction as the real and DMI as the imaginary part. As these phases are cycled, the current reveals the ground-state degeneracies that distinguish the ordered and disordered phases. We discuss experimental requirements needed to realize the Floquet gauge pump in a synthetic quantum spin system of interacting trapped ions.


title: On the chromatic numbers of signed triangular and hexagonal grids

abstract: A signed graph is a simple graph with two types of edges. Switching a vertex $v$ of a signed graph corresponds to changing the type of each edge incident to $v$. A homomorphism from a signed graph $G$ to another signed graph $H$ is a mapping $\varphi: V(G) \rightarrow V(H)$ such that, after switching any number of the vertices of $G$, $\varphi$ maps every edge of $G$ to an edge of the same type in $H$. The chromatic number $\chi_s(G)$ of a signed graph $G$ is the order of a smallest signed graph $H$ such that there is a homomorphism from $G$ to $H$. We show that the chromatic number of signed triangular grids is at most 10 and the chromatic number of signed hexagonal grids is at most 4.


title: Polynomial-Time Algorithms for Counting and Sampling Markov Equivalent DAGs

abstract: Counting and uniform sampling of directed acyclic graphs (DAGs) from a Markov equivalence class are fundamental tasks in graphical causal analysis. In this paper, we show that these tasks can be performed in polynomial time, solving a long-standing open problem in this area. Our algorithms are effective and easily implementable. Experimental results show that the algorithms significantly outperform state-of-the-art methods.


title: Compiler Design for Distributed Quantum Computing

abstract: In distributed quantum computing architectures, with the network and communications functionalities provided by the Quantum Internet, remote quantum processing units (QPUs) can communicate and cooperate for executing computational tasks that single NISQ devices cannot handle by themselves. To this aim, distributed quantum computing requires a new generation of quantum compilers, for mapping any quantum algorithm to any distributed quantum computing architecture. With this perspective, in this paper, we first discuss the main challenges arising with compiler design for distributed quantum computing. Then, we analytically derive an upper bound of the overhead induced by quantum compilation for distributed quantum computing. The derived bound accounts for the overhead induced by the underlying computing architecture as well as the additional overhead induced by the sub-optimal quantum compiler -- expressly designed through the paper to achieve three key features, namely, general-purpose, efficient and effective. Finally, we validate the analytical results and we confirm the validity of the compiler design through an extensive performance analysis.


title: Computational Overhead of Locality Reduction in Binary Optimization Problems

abstract: Recently, there has been considerable interest in solving optimization problems by mapping these onto a binary representation, sparked mostly by the use of quantum annealing machines. Such binary representation is reminiscent of a discrete physical two-state system, such as the Ising model. As such, physics-inspired techniques -- commonly used in fundamental physics studies -- are ideally suited to solve optimization problems in a binary format. While binary representations can be often found for paradigmatic optimization problems, these typically result in k-local higher-order unconstrained binary optimization cost functions. In this work, we discuss the effects of locality reduction needed for the majority of the currently available quantum and quantum-inspired solvers that can only accommodate 2-local (quadratic) cost functions. General locality reduction approaches require the introduction of ancillary variables which cause an overhead over the native problem. Using a parallel tempering Monte Carlo solver on Microsoft Azure Quantum, as well as k-local binary problems with planted solutions, we show that post reduction to a corresponding 2-local representation the problems become considerably harder to solve. We further quantify the increase in computational hardness introduced by the reduction algorithm by measuring the variation of number of variables, statistics of the coefficient values, and the population annealing entropic family size. Our results demonstrate the importance of avoiding locality reduction when solving optimization problems.


title: Evidence for the coexistence of time-reversal symmetry breaking and Bardeen-Cooper-Schrieffer-like superconductivity in La$_{7}$Pd$_{3}$

abstract: Time-reversal symmetry breaking (TRSB) with a Bardeen-Cooper-Schrieffer (BCS) -like superconductivity occurs in a small, but growing number of noncentrosymmetric (NCS) materials, although the mechanism is poorly understood. We present heat capacity, magnetization, resistivity, and muon spin resonance/relaxation ($\mu$SR) measurements on polycrystalline samples of NCS La$_{7}$Pd$_{3}$. Transverse-field $\mu$SR and heat capacity data show La$_{7}$Pd$_{3}$ is a type-II superconductor with a BCS-like gap structure, while zero-field $\mu$SR results provide evidence of TRSB. We discuss the implications of these results for both the La$_{7}X_{3}$ (where $X=$ Ni, Pd, Rh, Ir) group of superconductors and other CS and NCS superconductors for which TRSB has been observed.


title: Detecting Reconnection Events in Kinetic Vlasov Hybrid Simulations Using Clustering Techniques

abstract: Kinetic turbulence in magnetized space plasmas has been extensively studied via in situ observations, numerical simulations and theoretical models. In this context, a key point concerns the formation of coherent current structures and their disruption through magnetic reconnection. We present automatic techniques aimed at detecting reconnection events in large data set of numerical simulations. We make use of clustering techniques known as K-means and DBscan (usually referred in literature as unsupervised machine learning approaches), and other methods based on thresholds of standard reconnection proxies. All our techniques use also a threshold on the aspect ratio of the regions selected. We test the performance of our algorithms. We propose an optimal aspect ratio to be used in the automated machine learning algorithm: AR=18. The performance of the unsupervised approach results to be strongly competitive with respect to those of other methods based on thresholds of standard reconnection proxies.


title: Topological Bulk Lasing Modes Using an Imaginary Gauge Field

abstract: Topological edge modes, which are robust against disorders, have been used to enhance the spatial stability of lasers. Recently, it was revealed that topological lasers can be further stabilized using a novel topological phase in non-Hermitian photonic topological insulators. Here we propose a procedure to realize topologically protected modes extended over a d-dimensional bulk by introducing an imaginary gauge field. This generalizes the idea of zero-energy extended modes in the one-dimensional Su-Schrieffer-Heeger lattice into higher dimensional lattices allowing a d-dimensional bulky mode that is topologically protected. Furthermore, we numerically demonstrate that the topological bulk lasing mode can achieve high temporal stability superior to topological edge mode lasers. In the exemplified topological extended mode in the kagome lattice, we show that large regions of stability exist in its parameter space.


title: Parabolic Harnack estimates for anisotropic slow diffusion

abstract: We prove a Harnack inequality for positive solutions of a parabolic equation with slow anisotropic spatial diffusion. After identifying its natural scalings, we reduce the problem to a Fokker-Planck equation and construct a self-similar Barenblatt solution. We exploit translation invariance to obtain positivity near the origin via a self-iteration method and deduce a sharp anisotropic expansion of positivity. This eventually yields a scale invariant Harnack inequality in an anisotropic geometry dictated by the speed of the diffusion coefficients. As a corollary, we infer H\"older continuity, an elliptic Harnack inequality and a Liouville theorem.


title: Hate Speech detection in the Bengali language: A dataset and its baseline evaluation

abstract: Social media sites such as YouTube and Facebook have become an integral part of everyone's life and in the last few years, hate speech in the social media comment section has increased rapidly. Detection of hate speech on social media websites faces a variety of challenges including small imbalanced data sets, the findings of an appropriate model and also the choice of feature analysis method. further more, this problem is more severe for the Bengali speaking community due to the lack of gold standard labelled datasets. This paper presents a new dataset of 30,000 user comments tagged by crowd sourcing and varified by experts. All the comments are collected from YouTube and Facebook comment section and classified into seven categories: sports, entertainment, religion, politics, crime, celebrity and TikTok & meme. A total of 50 annotators annotated each comment three times and the majority vote was taken as the final annotation. Nevertheless, we have conducted base line experiments and several deep learning models along with extensive pre-trained Bengali word embedding such as Word2Vec, FastText and BengFastText on this dataset to facilitate future research opportunities. The experiment illustrated that although all deep learning models performed well, SVM achieved the best result with 87.5% accuracy. Our core contribution is to make this benchmark dataset available and accessible to facilitate further research in the field of in the field of Bengali hate speech detection.


title: Height function delocalisation on cubic planar graphs

abstract: The interest is in models of integer-valued height functions on shift-invariant planar graphs whose maximum degree is three. We prove delocalisation for models induced by convex nearest-neighbour potentials, under the condition that each potential function is an excited potential, that is, a convex symmetric potential function $V$ with the property that $V(\pm1)\leq V(0)+\log2$. Examples of such models include the discrete Gaussian and solid-on-solid models at inverse temperature $\beta\leq\log2$, as well as the uniformly random $K$-Lipschitz function for fixed $K\in\mathbb N$. In fact, $\beta V$ is an excited potential for any convex symmetric potential function $V$ whenever $\beta$ is sufficiently small. To arrive at the result, we develop a new technique for symmetry breaking, and then study the geometric percolation properties of sets of the form $\{\varphi\geq a\}$ and $\{\varphi\leq a\}$, where $\varphi$ is the random height function and $a$ a constant. Along the same lines, we derive delocalisation for models induced by convex symmetric nearest-neighbour potentials which force the parity of the height of neighbouring vertices to be distinct. This includes models of uniformly random graph homomorphisms on the honeycomb lattice and the truncated square tiling, as well as on the same graphs with each edge replaced by $N$ edges linked in series. The latter resembles cable-graph constructions which appear in the analysis of the Gaussian free field.


title: PCT: Point cloud transformer

abstract: The irregular domain and lack of ordering make it challenging to design deep neural networks for point cloud processing. This paper presents a novel framework named Point Cloud Transformer(PCT) for point cloud learning. PCT is based on Transformer, which achieves huge success in natural language processing and displays great potential in image processing. It is inherently permutation invariant for processing a sequence of points, making it well-suited for point cloud learning. To better capture local context within the point cloud, we enhance input embedding with the support of farthest point sampling and nearest neighbor search. Extensive experiments demonstrate that the PCT achieves the state-of-the-art performance on shape classification, part segmentation and normal estimation tasks.


title: Trajectory Planning Under Stochastic and Bounded Sensing Uncertainties Using Reachability Analysis

abstract: Trajectory planning under uncertainty is an active research topic. Previous works predict state and state estimation uncertainties along trajectories to check for collision safety. They assume either stochastic or bounded sensing uncertainties. However, GNSS pseudoranges are typically modeled to contain stochastic uncertainties with additional biases in urban environments. Thus, given bounds for the bias, the planner needs to account for both stochastic and bounded sensing uncertainties. In our prior work we presented a reachability analysis to predict state and state estimation uncertainties under stochastic and bounded uncertainties. However, we ignored the correlation between these uncertainties, leading to an imperfect approximation of the state uncertainty. In this paper we improve our reachability analysis by predicting state uncertainty as a function of independent quantities. We design a metric for the predicted uncertainty to compare candidate trajectories during planning. Finally, we validate the planner for GNSS-based urban navigation of fixed-wing UAS.


title: Isochrone fitting in the Gaia era. III. Distances, ages and masses from UniDAM using Gaia eDR3 data

abstract: We present estimates of distances, ages and masses for over 6 million stars. These estimates are derived from the combination of spectrophotometric data and Gaia eDR3 parallaxes. For that, we used the previously published Unified tool to estimate Distances, Ages, and Masses (UniDAM).


title: Topological Superconductivity Induced by Magnetic Texture Crystals

abstract: We present a detailed investigation of the topological phases and Majorana fermion (MF) excitations that arise from the bulk interplay between (un)conventional one/two-band spin-singlet superconductivity and a number of magnetic texture crystals. The latter define inhomogeneous magnetization profiles which consist of a periodically-repeating primitive cell. Here we focus on magnetic texture crystals with a primitive cell of the helix, whirl, and skyrmion types, which feature distinct symmetry properties. We identify a multitude of accessible topological phases which harbor flat, uni- or bi-directional, (quasi-)helical, or chiral MF edge modes. This rich variety originates from the interplay between topological phases with gapped and nodal bulk energy spectra. The types of the emerging topological superconducting phases and the features of the arising MFs are solely determined by the properties and compatibility of the so-called magnetic and pairing point/space groups. Our analysis is general and does not rely on specific parameters of the models employed here to exemplify the topological scenarios which become accessible. Therefore, our results can be extended to systems with multiple bands, are relevant for a wide range of layered materials and hybrid devices, and provide predictions for strong, weak and crystalline topological phases.

categories: cond-mat.supr-con cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el


title: Five Psycholinguistic Characteristics for Better Interaction with Users

abstract: When two people pay attention to each other and are interested in what the other has to say or write, they almost instantly adapt their writing/speaking style to match the other. For a successful interaction with a user, chatbots and dialogue systems should be able to do the same. We propose a framework consisting of five psycholinguistic textual characteristics for better human-computer interaction. We describe the annotation processes used for collecting the data, and benchmark five binary classification tasks, experimenting with different training sizes and model architectures. The best architectures noticeably outperform several baselines and achieve macro-averaged F$_1$-scores between 72\% and 96\% depending on the language and the task. The proposed framework proved to be fairly easy to model for various languages even with small amount of manually annotated data if right architectures are used.


title: 6 GHz hyperfast rotation of an optically levitated nanoparticle in vacuum

abstract: We report an experimental observation of a record-breaking ultra-high rotation frequency about 6 GHz in an optically levitated nanoparticle system. We optically trap a nanoparticle in the gravity direction with a high numerical aperture (NA) objective lens, which shows significant advantages in compensating the influences of the scattering force and the photophoretic force on the trap, especially at intermediate pressure (about 100 Pa). This allows us to trap a nanoparticle from atmospheric to low pressure ($10^{-3}$ Pa) without using feedback cooling. We measure a highest rotation frequency about 4.3 GHz of the trapped nanoparticle without feedback cooling and a 6 GHz rotation with feedback cooling, which is the fastest mechanical rotation ever reported to date. Our work provides useful guides for efficiently observing hyperfast rotation in the optical levitation system, and may find various applications such as in ultrasensitive torque detection, probing vacuum friction, and testing unconventional decoherence theories.


title: Jacob's ladder as generator of new class of iterated $L_2$-orthogonal systems and their dependence on the Riemann's function

abstract: In this paper new classes of $L_2$-orthogonal functions are constructed as iterated $L_2$-orthogonal systems. In order to do this we use the theory of the Riemann's zeta-function as well as our theory of Jacob's ladders. The main result is new one in the theory of the Riemann's zeta-function and simultaneously in the theory of $L_2$-orthogonal systems.


title: Quantum statistics effects and fluctuations of particle numbers near the critical point of nuclear matter

abstract: Equation of state with quantum statistics corrections is derived for a multi-component gas of particles interacting through the repulsive and attractive van der Waals (vdW) forces up to first few orders over a small parameter $\delta \approx \hbar^3 n(mT)^{-3/2}[g(1- bn)]^{-1}$, where $n$ and $T$ are the particle number density and temperature, $m$ and $g$ the particle mass and degeneracy factor. The parameter $b$ corresponds to the vdW excluded volume. For interacting system of Fermi nucleon and Bose $\alpha$ particles, a small impurity of $\alpha$ particles to the nucleon system at leading first order in both $\alpha $ particle and nucleon small parameters $\delta $ does not change much the basic results for the symmetric nuclear matter. The particle number fluctuations $\omega$ determined by the isothermal in-compressibility $\mathcal{K}(n,T)$ can be obtained analytically at the same first order quantum-statistics approximation for symmetric nucleon matter. Our approximate analytical results appear to be in good agreement with the accurate numerical calculations.


title: Multi-FinGAN: Generative Coarse-To-Fine Sampling of Multi-Finger Grasps

abstract: While there exists many methods for manipulating rigid objects with parallel-jaw grippers, grasping with multi-finger robotic hands remains a quite unexplored research topic. Reasoning and planning collision-free trajectories on the additional degrees of freedom of several fingers represents an important challenge that, so far, involves computationally costly and slow processes. In this work, we present Multi-FinGAN, a fast generative multi-finger grasp sampling method that synthesizes high quality grasps directly from RGB-D images in about a second. We achieve this by training in an end-to-end fashion a coarse-to-fine model composed of a classification network that distinguishes grasp types according to a specific taxonomy and a refinement network that produces refined grasp poses and joint angles. We experimentally validate and benchmark our method against a standard grasp-sampling method on 790 grasps in simulation and 20 grasps on a real Franka Emika Panda. All experimental results using our method show consistent improvements both in terms of grasp quality metrics and grasp success rate. Remarkably, our approach is up to 20-30 times faster than the baseline, a significant improvement that opens the door to feedback-based grasp re-planning and task informative grasping. Code is available at https://irobotics.aalto.fi/multi-fingan/.


title: Some stability inequalities for hybrid inverse problems

abstract: We study some hybrid inverse problems associated to BVP's for Schr\"odinger and Helmholtz type equations. The inverse problems we consider consist in the determination of coefficients from the knowledge of internal energies. We establish local Lipschitz stability inequalities as well as conditional H\"older stability inequalities.


title: A Near-Conformal Composite Higgs Model

abstract: We analyze a composite Higgs model based on the confining $SU(3)$ gauge theory with $N_f = 8$ Dirac fermions in the fundamental representation. This gauge theory has been studied on the lattice and shown to be well described by a dilaton effective field theory (EFT). Here we modify the EFT by assigning standard-model quantum numbers such that four of the composite pseudo-Nambu-Goldstone boson (pNGB) fields form the standard-model Higgs doublet, by coupling it to the top quark, and by adding to the potential a term that triggers electroweak symmetry breaking. The model contains a pNGB Higgs boson, a set of heavier pNGBs, and an approximate dilaton in the same mass range. We study the phenomenology of the model, and discuss the amount of tuning required to insure consistency with current direct and indirect bounds on new physics, highlighting the role of the dilaton field.


title: A Generalization of Transformer Networks to Graphs

abstract: We propose a generalization of transformer neural network architecture for arbitrary graphs. The original transformer was designed for Natural Language Processing (NLP), which operates on fully connected graphs representing all connections between the words in a sequence. Such architecture does not leverage the graph connectivity inductive bias, and can perform poorly when the graph topology is important and has not been encoded into the node features. We introduce a graph transformer with four new properties compared to the standard model. First, the attention mechanism is a function of the neighborhood connectivity for each node in the graph. Second, the positional encoding is represented by the Laplacian eigenvectors, which naturally generalize the sinusoidal positional encodings often used in NLP. Third, the layer normalization is replaced by a batch normalization layer, which provides faster training and better generalization performance. Finally, the architecture is extended to edge feature representation, which can be critical to tasks s.a. chemistry (bond type) or link prediction (entity relationship in knowledge graphs). Numerical experiments on a graph benchmark demonstrate the performance of the proposed graph transformer architecture. This work closes the gap between the original transformer, which was designed for the limited case of line graphs, and graph neural networks, that can work with arbitrary graphs. As our architecture is simple and generic, we believe it can be used as a black box for future applications that wish to consider transformer and graphs.


title: RainNet: A Large-Scale Imagery Dataset and Benchmark for Spatial Precipitation Downscaling

abstract: AI-for-science approaches have been applied to solve scientific problems (e.g., nuclear fusion, ecology, genomics, meteorology) and have achieved highly promising results. Spatial precipitation downscaling is one of the most important meteorological problem and urgently requires the participation of AI. However, the lack of a well-organized and annotated large-scale dataset hinders the training and verification of more effective and advancing deep-learning models for precipitation downscaling. To alleviate these obstacles, we present the first large-scale spatial precipitation downscaling dataset named RainNet, which contains more than $62,400$ pairs of high-quality low/high-resolution precipitation maps for over $17$ years, ready to help the evolution of deep learning models in precipitation downscaling. Specifically, the precipitation maps carefully collected in RainNet cover various meteorological phenomena (e.g., hurricane, squall), which is of great help to improve the model generalization ability. In addition, the map pairs in RainNet are organized in the form of image sequences ($720$ maps per month or 1 map/hour), showing complex physical properties, e.g., temporal misalignment, temporal sparse, and fluid properties. Furthermore, two deep-learning-oriented metrics are specifically introduced to evaluate or verify the comprehensive performance of the trained model (e.g., prediction maps reconstruction accuracy). To illustrate the applications of RainNet, 14 state-of-the-art models, including deep models and traditional approaches, are evaluated. To fully explore potential downscaling solutions, we propose an implicit physical estimation benchmark framework to learn the above characteristics. Extensive experiments demonstrate the value of RainNet in training and evaluating downscaling models. Our dataset is available at https://neuralchen.github.io/RainNet/.


title: Multi-Wavelength, Optical (VI) and Near-Infrared (JHK) Calibration of the Tip of the Red Giant Branch Method based on Milky Way Globular Clusters

abstract: Using high precision ground-based photometry for 46 low-reddening Galactic globular clusters, in conjunction with Gaia DR2 proper motions for member star selection, we have calibrated the zero point of the tip of the red giant branch (TRGB) method at two optical ($VI$) and three near-infrared ($JHK$) wavelengths. In doing so, we utilized the sharply-defined zero-age horizontal branch (ZAHB) of these clusters to relatively calibrate our cluster sample into a composite color-magnitude diagram spanning a wide range of metallicities, before setting the absolute zero point of this composite using the geometric detached eclipsing binary distance to the cluster $\omega$ Centauri. The $I-$band zero point we measure [$M_I = -4.056 \pm 0.02 \text{ (stat}) \pm 0.10 \text{ (sys)} $] agrees to within one sigma of the two previously published independent calibrations, using TRGB stars in the LMC [$M_I = $ -4.047 mag; Freedman et al. 2019, 2020] and in the maser galaxy NGC 4258 [$M_{F814W} = $ -4.051 mag; Jang et al. 2020]. We also find close agreement for our $J,H,K$ zero points to several literature studies.


title: SRoll3: A neural network approach to reduce large-scale systematic effects in the Planck High Frequency Instrument maps

abstract: In the present work, we propose a neural network based data inversion approach to reduce structured contamination sources, with a particular focus on the mapmaking for Planck High Frequency Instrument (Planck-HFI) data and the removal of large-scale systematic effects within the produced sky maps. The removal of contamination sources is rendered possible by the structured nature of these sources, which is characterized by local spatiotemporal interactions producing couplings between different spatiotemporal scales. We focus on exploring neural networks as a means of exploiting these couplings to learn optimal low-dimensional representations, optimized with respect to the contamination source removal and mapmaking objectives, to achieve robust and effective data inversion. We develop multiple variants of the proposed approach, and consider the inclusion of physics informed constraints and transfer learning techniques. Additionally, we focus on exploiting data augmentation techniques to integrate expert knowledge into an otherwise unsupervised network training approach. We validate the proposed method on Planck-HFI 545 GHz Far Side Lobe simulation data, considering ideal and non-ideal cases involving partial, gap-filled and inconsistent datasets, and demonstrate the potential of the neural network based dimensionality reduction to accurately model and remove large-scale systematic effects. We also present an application to real Planck-HFI 857 GHz data, which illustrates the relevance of the proposed method to accurately model and capture structured contamination sources, with reported gains of up to one order of magnitude in terms of contamination removal performance. Importantly, the methods developed in this work are to be integrated in a new version of the SRoll algorithm (SRoll3), and we describe here SRoll3 857 GHz detector maps that will be released to the community.


title: Parallel WaveNet conditioned on VAE latent vectors

abstract: Recently the state-of-the-art text-to-speech synthesis systems have shifted to a two-model approach: a sequence-to-sequence model to predict a representation of speech (typically mel-spectrograms), followed by a 'neural vocoder' model which produces the time-domain speech waveform from this intermediate speech representation. This approach is capable of synthesizing speech that is confusable with natural speech recordings. However, the inference speed of neural vocoder approaches represents a major obstacle for deploying this technology for commercial applications. Parallel WaveNet is one approach which has been developed to address this issue, trading off some synthesis quality for significantly faster inference speed. In this paper we investigate the use of a sentence-level conditioning vector to improve the signal quality of a Parallel WaveNet neural vocoder. We condition the neural vocoder with the latent vector from a pre-trained VAE component of a Tacotron 2-style sequence-to-sequence model. With this, we are able to significantly improve the quality of vocoded speech.


title: Probing sub-eV Dark Matter decays with PTOLEMY

abstract: When the Dark Matter mass is below the eV-scale, its cosmological occupation number exceeds the ones of photons from the cosmic microwave background as well as of relic neutrinos. If such Dark Matter decays to pairs of neutrinos, it implies that experiments that seek the detection of the cosmic neutrino background may as well be sensitive to this additional form of "dark radiation". Here we study the prospects for detection taking into account various options for the forecasted performance of the future PTOLEMY experiment. From a detailed profile likelihood analysis we find that Dark Matter decays with lifetime as large as $10^4$ Gyr or a sub-% Dark Matter fraction decaying today can be discovered. The prospects are facilitated by the distinct spectral event shape that is introduced from galactic and cosmological neutrino dark radiation fluxes. In the process we also clarify the importance of Pauli-blocking in the Dark Matter decay. The scenarios presented in this work can be considered early physics targets in the development of these instruments with relaxed demands on performance and energy resolution.


title: Trellis Code Error Exponent From Results for Asynchronous Multiple Access Channels

abstract: An asynchronous multiple access error exponent result implicates a new result for time invariant trellis codes of memory 1.


title: Unified Emergence of Energy Scales and Cosmic Inflation

abstract: In the quest for unification of the Standard Model with gravity, classical scale invariance can be utilized to dynamically generate the Planck mass $M_\mathrm{Pl}$. Then, the relation of Planck scale physics to the scale of electroweak symmetry breaking $\mu_H$ requires further explanation. In this paper, we propose a model that uses the spontaneous breaking of scale invariance in the scalar sector as a unified origin for dynamical generation of both scales. Using the Gildener-Weinberg approximation, only one scalar acquires a vacuum expectation value of $v_S \sim (10^{16-17})\,\mathrm{GeV}$, thus radiatively generating $M_\mathrm{Pl} \approx \beta_S^{1/2} v_S$ and $\mu_H$ via the neutrino option with right handed neutrino masses $m_N = y_M v_S \sim 10^7 \,\mathrm{GeV}$. Consequently, active SM neutrinos are given a mass with the inclusion of a type-I seesaw mechanism. Furthermore, we adopt an unbroken $Z_2$ symmetry and a $Z_2$-odd set of right-handed Majorana neutrinos $\chi$ that do not take part in the neutrino option and are able to produce the correct dark matter relic abundance (dominantly) via inflaton decay. The model also describes cosmic inflation and the inflationary CMB observables are predicted to interpolate between those of $R^2$ and linear chaotic inflationary model and are thus well within the strongest experimental constraints.


title: Intrusion Detection and identification System Design and Performance Evaluation for Industrial SCADA Networks

abstract: In this paper, we present a study that proposes a three-stage classifier model which employs a machine learning algorithm to develop an intrusion detection and identification system for tens of different types of attacks against industrial SCADA networks. The machine learning classifier is trained and tested on the data generated using the laboratory prototype of a gas pipeline SCADA network. The dataset consists of three attack groups and seven different attack classes or categories. The same dataset further provides signatures of 35 different types of sub-attacks which are related to those seven attack classes. The study entailed the design of three-stage machine learning classifier as a misuse intrusion detection system to detect and identify specifically each of the 35 attack subclasses. The first stage of the classifier decides if a record is associated with normal operation or an attack signature. If the record is found to belong to an attack signature, then in the second stage, it is classified into one of seven attack classes. Based on the identified attack class as determined by the output from the second stage classifier, the attack record is provided for a third stage sub-attack classification, where seven different classifiers are employed. The output from the third stage classifier identifies the sub-attack type to which the record belongs. Simulation results indicate that designs exploring specialization to domains or executing the classification in multiple stages versus single-stage designs are promising for problems where there are tens of classes. Comparison with studies in the literature also indicated that the multi-stage classifier performed markedly better.


title: Efficient CNN-LSTM based Image Captioning using Neural Network Compression

abstract: Modern Neural Networks are eminent in achieving state of the art performance on tasks under Computer Vision, Natural Language Processing and related verticals. However, they are notorious for their voracious memory and compute appetite which further obstructs their deployment on resource limited edge devices. In order to achieve edge deployment, researchers have developed pruning and quantization algorithms to compress such networks without compromising their efficacy. Such compression algorithms are broadly experimented on standalone CNN and RNN architectures while in this work, we present an unconventional end to end compression pipeline of a CNN-LSTM based Image Captioning model. The model is trained using VGG16 or ResNet50 as an encoder and an LSTM decoder on the flickr8k dataset. We then examine the effects of different compression architectures on the model and design a compression architecture that achieves a 73.1% reduction in model size, 71.3% reduction in inference time and a 7.7% increase in BLEU score as compared to its uncompressed counterpart.


title: TESS observations of Cepheid stars: first light results

abstract: We present the first analysis of Cepheid stars observed by the TESS space mission in Sectors 1 to 5. Our sample consists of 25 pulsators: ten fundamental mode, three overtone and two double-mode classical Cepheids, plus three Type II and seven anomalous Cepheids. The targets were chosen from fields with different stellar densities, both from the Galactic field and from the Magellanic System. Three targets have 2-minute cadence light curves available by the TESS Science Processing Operations Center: for the rest, we prepared custom light curves from the full-frame images with our own differential photometric FITSH pipeline. Our main goal was to explore the potential and the limitations of TESS concerning the various subtypes of Cepheids. We detected many low amplitude features: weak modulation, period jitter, and timing variations due to light-time effect. We also report signs of non-radial modes and the first discovery of such a mode in an anomalous Cepheid, the overtone star XZ Cet, which we then confirmed with ground-based multicolor photometric measurements. We prepared a custom photometric solution to minimize saturation effects in the bright fundamental-mode classical Cepheid, $\beta$ Dor with the lightkurve software, and we revealed strong evidence of cycle-to-cycle variations in the star. In several cases, however, fluctuations in the pulsation could not be distinguished from instrumental effects, such as contamination from nearby sources which also varies between sectors. Finally, we discuss how precise light curve shapes will be crucial not only for classification purposes but also to determine physical properties of these stars.


title: Factorized lightcone expansion of conformal blocks

abstract: We present a factorized decomposition of 4-point scalar conformal blocks near the lightcone, which applies to arbitrary intermediate spin and general spacetime dimensions. Then we discuss the systematic expansion in large intermediate spin and the resummations of the large-spin tails of Regge trajectories. The basic integrals for the Lorentzian inversion are given by Wilson functions.


title: Kaluza-Klein from Colour-Kinematics Duality for Massive Fields

abstract: We consider a broad class of massive four dimensional effective theories describing an infinite tower of charged massive spin 1 states, interacting with massless spin 1 and spin 0. The spectrum is chosen to be the same as that appears in the Kaluza-Klein theory reduction of 5d Yang-Mills to ensure the absence of any spurious poles in a possible double copy formulation. The effective theories are characterized by multiple different couplings between different fields which are unconstrained by symmetries and low energy criteria. Remarkably, by demanding that the scattering amplitudes preserve colour-kinematics duality for different scattering processes, required for the existence of a massive double copy, we find that our 4d Lagrangian is fixed uniquely to the Kaluza-Klein compactification of 5d Yang-Mills theory together with its known double copy consistent higher derivative operators.


title: Deep Molecular Dreaming: Inverse machine learning for de-novo molecular design and interpretability with surjective representations

abstract: Computer-based de-novo design of functional molecules is one of the most prominent challenges in cheminformatics today. As a result, generative and evolutionary inverse designs from the field of artificial intelligence have emerged at a rapid pace, with aims to optimize molecules for a particular chemical property. These models 'indirectly' explore the chemical space; by learning latent spaces, policies, distributions or by applying mutations on populations of molecules. However, the recent development of the SELFIES string representation of molecules, a surjective alternative to SMILES, have made possible other potential techniques. Based on SELFIES, we therefore propose PASITHEA, a direct gradient-based molecule optimization that applies inceptionism techniques from computer vision. PASITHEA exploits the use of gradients by directly reversing the learning process of a neural network, which is trained to predict real-valued chemical properties. Effectively, this forms an inverse regression model, which is capable of generating molecular variants optimized for a certain property. Although our results are preliminary, we observe a shift in distribution of a chosen property during inverse-training, a clear indication of PASITHEA's viability. A striking property of inceptionism is that we can directly probe the model's understanding of the chemical space it was trained on. We expect that extending PASITHEA to larger datasets, molecules and more complex properties will lead to advances in the design of new functional molecules as well as the interpretation and explanation of machine learning models.


title: Digraphs and variable degeneracy

abstract: Let $D$ be a digraph, let $p \geq 1$ be an integer, and let $f: V(D) \to \mathbb{N}_0^p$ be a vector function with $f=(f_1,f_2,\ldots,f_p)$. We say that $D$ has an $f$-partition if there is a partition $(D_1,D_2,\ldots,D_p)$ into induced subdigraphs of $D$ such that for all $i \in [1,p]$, the digraph $D_i$ is weakly $f_i$-degenerate, that is, in every non-empty subdigraph $D'$ of $D_i$ there is a vertex $v$ such that $\min\{d_{D'}^+(v), d_{D'}^-(v)\} < f_i(v)$. In this paper, we prove that the condition $f_1(v) + f_2(v) + \ldots + f_p(v) \geq \max \{d_D^+(v),d_D^-(v)\}$ for all $v \in V(D)$ is almost sufficient for the existence of an $f$-partition and give a full characterization of the bad pairs $(D,f)$. Moreover, we describe a polynomial time algorithm that (under the previous conditions) either verifies that $(D,f)$ is a bad pair or finds an $f$-partition. Among other applications, this leads to a generalization of Brooks' Theorem as well as the list-version of Brooks' Theorem for digraphs, where a coloring of digraph is a partition of the digraph into acyclic induced subdigraphs. We furthermore obtain a result bounding the $s$-degenerate chromatic number of a digraph in terms of the maximum of maximum in-degree and maximum out-degree.


title: Attosecond Pulse-shaping using a seeded free-electron laser

abstract: Attosecond pulses are fundamental for the investigation of valence and core-electron dynamics on their natural timescale. At present the reproducible generation and characterisation of attosecond waveforms has been demonstrated only through the process of high-order harmonic generation. Several methods for the shaping of attosecond waveforms have been proposed, including metallic filters, multilayer mirrors and manipulation of the driving field. However, none of these approaches allow for the flexible manipulation of the temporal characteristics of the attosecond waveforms, and they suffer from the low conversion efficiency of the high-order harmonic generation process. Free Electron Lasers, on the contrary, deliver femtosecond, extreme ultraviolet and X-ray pulses with energies ranging from tens of $\mathrm{\mu}$J to a few mJ. Recent experiments have shown that they can generate sub-fs spikes, but with temporal characteristics that change shot-to-shot. Here we show the first demonstration of reproducible generation of high energy ($\mathrm{\mu}$J level) attosecond waveforms using a seeded Free Electron Laser. We demonstrate amplitude and phase manipulation of the harmonic components of an attosecond pulse train in combination with a novel approach for its temporal reconstruction. The results presented here open the way to perform attosecond time-resolved experiments with Free Electron Lasers.


title: Approximating inverse cumulative distribution functions to produce approximate random variables

abstract: For random variables produced through the inverse transform method, approximate random variables are introduced, which are produced by approximations to a distribution's inverse cumulative distribution function. These approximations are designed to be computationally inexpensive, and much cheaper than exact library functions, and thus highly suitable for use in Monte Carlo simulations. Two approximations are presented for the Gaussian distribution: a piecewise constant on equally spaced intervals, and a piecewise linear using geometrically decaying intervals. The error of the approximations are bounded and the convergence demonstrated, and the computational savings measured for C and C++ implementations. Implementations tailored for Intel and Arm hardwares are inspected, alongside hardware agnostic implementations built using OpenMP. The savings are incorporated into a nested multilevel Monte Carlo framework with the Euler-Maruyama scheme to exploit the speed ups without losing accuracy, offering speed ups by a factor of 5--7. These ideas are empirically extended to the Milstein scheme, and the Cox-Ingersoll-Ross process' non central chi-squared distribution, which offer speed ups by a factor of 250 or more.