title: Quantifying the quantum heat contribution from a driven superconducting circuit

abstract: Unveiling the impact in thermodynamics of the phenomena specific to quantum mechanics is a crucial step to identify fundamental costs for quantum operations and quantum advantages in heat engines. We propose a two-reservoir setup to detect the quantum component in the heat flow exchanged by a coherently driven atom with its thermal environment. Tuning the driving parameters switches on and off the quantum and classical contributions to the heat flows, enabling their independent characterization. We demonstrate the feasibility of the measurement in a circuit-QED setup. Our results pave the road towards the first experimental verification of this quantum thermodynamic signature ubiquitous in quantum technologies.


title: Extreme quasars as distance indicators in cosmology

abstract: Quasars accreting matter at very high rates (known as extreme Population A [xA] quasars, possibly associated with super-Eddington accreting massive black holes) may provide a new class of distance indicators covering cosmic epochs from present day up to less than 1 Gyr from the Big Bang. At a more fundamental level, xA quasars are of special interest in studies of the physics of AGNs and host galaxy evolution. However, their observational properties are largely unknown. xA quasars can be identified in relatively large numbers from major optical surveys over a broad range of redshifts, thanks to selection criteria defined from the systematic changes along the quasars main sequence. It has been possible to build a sample of about 250 quasars at low and intermediate redshift, and larger samples can be easily selected from the latest data releases of the Sloan Digital Sky Survey. A large sample can clarify the main properties of xA quasars which are expected - unlike the general population of quasars - to radiate at an extreme, well defined Eddington ratio with small scatter. As a result of the small scatter in Eddington ratio shown by xA quasars, we propose a method to derive the main cosmological parameters based on redshift-independent "virial luminosity" estimates from measurements of emission line widths, roughly equivalent to the luminosity estimates based from line width in early and late type galaxies. A major issue related to the cosmological application of the xA quasar luminosity estimates from line widths is the identification of proper emission lines whose broadening is predominantly virial over a wide range of redshift and luminosity. We report on preliminary developments using the AlIII 1860 intermediate ionization line and the Hydrogen Balmer line H-beta as virial broadening estimators, and we briefly discuss the perspective of the method based on xA quasars.


title: Distinguished Limits and Vibrogenic Force revealed by Newton's Equation with Oscillating Force

abstract: In this paper, we analyse the basic ideas of Vibrodynamics and the two-timing method. To make our analysis most instructive, we have chosen the Newton's equation with a general oscillating force. We deal with its asymptotic solutions in the high frequency limit. Our treatment is simple but general. The targets of our study are \emph{the distinguished limits} and \emph{the universal vibrogenic force}. The aim of \emph{the distinguished limit procedure} is to identify how the small parameter can appear in an equation. The proper appearance of a small parameter leads to \emph{valid successive approximations}, and, in particular, to closed systems of averaged equations. We show, that there are only two distinguished limits. This means that Newton's equation, with high-frequency forcing, has two types of interesting asymptotic solutions. The key item in the averaged equations for all distinguished limits is \emph{the unique vibrogenic force}. The current state-of-the-art in this area is: a large number of particular examples are well-known, effective and advanced general methods (like the Krylov-Bogolyubov approach) are well developed. However, the presented general and simple analysis of distinguished limits and the vibrogenic force, formulated as a compact practical guide, is novel. An advantage of our treatment is the possibility of its straightforward use for various ODEs and PDEs with oscillating coefficients.


title: Hyperbolic Model Reduction for Kinetic Equations

abstract: We make a brief historical review to the moment model reduction to the kinetic equations, particularly the Grad's moment method for Boltzmann equation. The focus is on the hyperbolicity of the reduced model, which is essential to the existence of its classical solution as a Cauchy problem. The theory of the framework we developed in last years is then introduced, which may preserve the hyperbolic nature of the kinetic equations with high universality. Some lastest progress on the comparison between models with/without hyperbolicity is presented to validate the hyperbolic moment models for rarefied gases.


title: Improving operational flexibility of integrated energy system with uncertain renewable generations considering thermal inertia of buildings

abstract: Insufficient flexibility in system operation caused by traditional "heat-set" operating modes of combined heat and power (CHP) units in winter heating periods is a key issue that limits renewable energy consumption. In order to reduce the curtailment of renewable energy resources through improving the operational flexibility, a novel optimal scheduling model based on chance-constrained programming (CCP), aiming at minimizing the lowest generation cost, is proposed for a small-scale integrated energy system (IES) with CHP units, thermal power units, renewable generations and representative auxiliary equipments. In this model, due to the uncertainties of renewable generations including wind turbines and photovoltaic units, the probabilistic spinning reserves are supplied in the form of chance-constrained; from the perspective of user experience, a heating load model is built with consideration of heat comfort and inertia in buildings. To solve the model, a solution approach based on sequence operation theory (SOT) is developed, where the original CCP-based scheduling model is tackled into a solvable mixed-integer linear programming (MILP) formulation by converting a chance constraint into its deterministic equivalence class, and thereby is solved via the CPLEX solver. The simulation results on the modified IEEE 30-bus system demonstrate that the presented method manages to improve operational flexibility of the IES with uncertain renewable generations by comprehensively leveraging thermal inertia of buildings and different kinds of auxiliary equipments, which provides a fundamental way for promoting renewable energy consumption.


title: Efficient Optimization Method for Finding Minimum Energy Paths of Magnetic Transitions

abstract: Efficient algorithms for the calculation of minimum energy paths of magnetic transitions are implemented within the geodesic nudged elastic band (GNEB) approach. While an objective function is not available for GNEB and a traditional line search can, therefore, not be performed, the use of limited memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) and conjugate gradient algorithms in conjunction with orthogonal spin optimization (OSO) approach is shown to greatly outperform the previously used velocity projection and dissipative Landau-Lifschitz dynamics optimization methods. The implementation makes use of energy weighted springs for the distribution of the discretization points along the path and this is found to improve performance significantly. The various methods are applied to several test problems using a Heisenberg-type Hamiltonian, extended in some cases to include Dzyaloshinskii-Moriya and exchange interactions beyond nearest neighbours. Minimum energy paths are found for magnetization reversals in a nano-island, collapse of skyrmions in two-dimensional layers and annihilation of a chiral bobber near the surface of a three-dimensional magnet. The LBFGS-OSO method is found to outperform the dynamics based approaches by up to a factor of 8 in some cases.


title: Stochastic homogenization on randomly perforated domains

abstract: We study the existence of uniformly bounded extension and trace operators for $W^{1,p}$-functions on randomly perforated domains, where the geometry is assumed to be stationary ergodic. Such extension and trace operators are important for compactness in stochastic homogenization. In contrast to former approaches and results, we use very weak assumptions on the geometry which we call local $(\delta,M)$-regularity, isotropic cone mixing and bounded average connectivity. The first concept measures local Lipschitz regularity of the domain while the second measures the mesoscopic distribution of void space. The third is the most tricky part and measures the "mesoscopic" connectivity of the geometry. In contrast to former approaches we do not require a minimal distance between the inclusions and we allow for globally unbounded Lipschitz constants and percolating holes. We will illustrate our method by applying it to the Boolean model based on a Poisson point process and to a Delaunay pipe process. We finally introduce suitable Sobolev spaces on $\mathbb{R}^d$ and $\Omega$ in order to construct a stochastic two-scale convergence method and apply the resulting theory to the homogenization of a $p$-Laplace problem on a randomly perforated domain.


title: The Enron Corpus: Where the Email Bodies are Buried?

abstract: To probe the largest public-domain email database for indicators of fraud, we apply machine learning and accomplish four investigative tasks. First, we identify persons of interest (POI), using financial records and email, and report a peak accuracy of 95.7%. Secondly, we find any publicly exposed personally identifiable information (PII) and discover 50,000 previously unreported instances. Thirdly, we automatically flag legally responsive emails as scored by human experts in the California electricity blackout lawsuit, and find a peak 99% accuracy. Finally, we track three years of primary topics and sentiment across over 10,000 unique people before, during and after the onset of the corporate crisis. Where possible, we compare accuracy against execution times for 51 algorithms and report human-interpretable business rules that can scale to vast datasets.


title: Assigning probabilities to non-Lipschitz mechanical systems

abstract: We present a method for assigning probabilities to the solutions of initial value problems that have a Lipschitz singularity. To illustrate the method, we focus on the following toy example: $\frac{d^2r(t)}{dt^2} = r^a$, $r(t=0) =0$, and $\frac{dr(t)}{dt}\mid_{r(t=0)} =0$, with $a \in ]0,1[$. This example has a physical interpretation as a mass in a uniform gravitational field on a frictionless, rigid dome of a particular shape; the case with $a=1/2$ is known as Norton's dome. Our approach is based on (1) finite difference equations, which are deterministic; (2) elementary techniques from alpha-theory, a simplified framework for non-standard analysis that allows us to study infinitesimal perturbations; and (3) a uniform prior on the canonical phase space. Our deterministic, hyperfinite grid model allows us to assign probabilities to the solutions of the initial value problem in the original, indeterministic model.


title: Advaita: Bug Duplicity Detection System

abstract: Bugs are prevalent in software development. To improve software quality, bugs are filed using a bug tracking system. Properties of a reported bug would consist of a headline, description, project, product, component that is affected by the bug and the severity of the bug. Duplicate bugs rate (% of duplicate bugs) are in the range from single digit (1 to 9%) to double digits (40%) based on the product maturity , size of the code and number of engineers working on the project. Duplicate bugs range are between 9% to 39% in some of the open source projects like Eclipse, Firefox etc. Detection of duplicity deals with identifying whether any two bugs convey the same meaning. This detection of duplicates helps in de-duplication. Detecting duplicate bugs help reduce triaging efforts and saves time for developers in fixing the issues. Traditional natural language processing techniques are less accurate in identifying similarity between sentences. Using the bug data present in a bug tracking system, various approaches were explored including several machine learning algorithms, to obtain a viable approach that can identify duplicate bugs, given a pair of sentences(i.e. the respective bug descriptions). This approach considers multiple sets of features viz. basic text statistical features, semantic features and contextual features. These features are extracted from the headline, description and component and are subsequently used to train a classification algorithm.


title: Saddlepoint approximations for spatial panel data models

abstract: We develop new higher-order asymptotic techniques for the Gaussian maximum likelihood estimator in a spatial panel data model, with fixed effects, time-varying covariates, and spatially correlated errors. Our saddlepoint density and tail area approximation feature relative error of order $O(1/(n(T-1)))$ with $n$ being the cross-sectional dimension and $T$ the time-series dimension. The main theoretical tool is the tilted-Edgeworth technique in a non-identically distributed setting. The density approximation is always non-negative, does not need resampling, and is accurate in the tails. Monte Carlo experiments on density approximation and testing in the presence of nuisance parameters illustrate the good performance of our approximation over first-order asymptotics and Edgeworth expansions. An empirical application to the investment-saving relationship in OECD (Organisation for Economic Co-operation and Development) countries shows disagreement between testing results based on first-order asymptotics and saddlepoint techniques.


title: MetaSelector: Meta-Learning for Recommendation with User-Level Adaptive Model Selection

abstract: Recommender systems often face heterogeneous datasets containing highly personalized historical data of users, where no single model could give the best recommendation for every user. We observe this ubiquitous phenomenon on both public and private datasets and address the model selection problem in pursuit of optimizing the quality of recommendation for each user. We propose a meta-learning framework to facilitate user-level adaptive model selection in recommender systems. In this framework, a collection of recommenders is trained with data from all users, on top of which a model selector is trained via meta-learning to select the best single model for each user with the user-specific historical data. We conduct extensive experiments on two public datasets and a real-world production dataset, demonstrating that our proposed framework achieves improvements over single model baselines and sample-level model selector in terms of AUC and LogLoss. In particular, the improvements may lead to huge profit gain when deployed in online recommender systems.


title: SANST: A Self-Attentive Network for Next Point-of-Interest Recommendation

abstract: Next point-of-interest (POI) recommendation aims to offer suggestions on which POI to visit next, given a user's POI visit history. This problem has a wide application in the tourism industry, and it is gaining an increasing interest as more POI check-in data become available. The problem is often modeled as a sequential recommendation problem to take advantage of the sequential patterns of user check-ins, e.g., people tend to visit Central Park after The Metropolitan Museum of Art in New York City. Recently, self-attentive networks have been shown to be both effective and efficient in general sequential recommendation problems, e.g., to recommend products, video games, or movies. Directly adopting self-attentive networks for next POI recommendation, however, may produce sub-optimal recommendations. This is because vanilla self-attentive networks do not consider the spatial and temporal patterns of user check-ins, which are two critical features in next POI recommendation. To address this limitation, in this paper, we propose a model named SANST that incorporates spatio-temporal patterns of user check-ins into self-attentive networks. To incorporate the spatial patterns, we encode the relative positions of POIs into their embeddings before feeding the embeddings into the self-attentive network. To incorporate the temporal patterns, we discretize the time of POI check-ins and model the temporal relationship between POI check-ins by a relation-aware self-attention module. We evaluate the performance of our SANST model with three real-world datasets. The results show that SANST consistently outperforms the state-of-theart models, and the advantage in nDCG@10 is up to 13.65%.


title: Investigating Classification Techniques with Feature Selection For Intention Mining From Twitter Feed

abstract: In the last decade, social networks became most popular medium for communication and interaction. As an example, micro-blogging service Twitter has more than 200 million registered users who exchange more than 65 million posts per day. Users express their thoughts, ideas, and even their intentions through these tweets. Most of the tweets are written informally and often in slang language, that contains misspelt and abbreviated words. This paper investigates the problem of selecting features that affect extracting user's intention from Twitter feeds based on text mining techniques. It starts by presenting the method we used to construct our own dataset from extracted Twitter feeds. Following that, we present two techniques of feature selection followed by classification. In the first technique, we use Information Gain as a one-phase feature selection, followed by supervised classification algorithms. In the second technique, we use a hybrid approach based on forward feature selection algorithm in which two feature selection techniques employed followed by classification algorithms. We examine these two techniques with four classification algorithms. We evaluate them using our own dataset, and we critically review the results.


title: Detecting State Transitions of a Markov Source: Sampling Frequency and Age Trade-off

abstract: We consider a finite-state Discrete-Time Markov Chain (DTMC) source that can be sampled for detecting the events when the DTMC transits to a new state. Our goal is to study the trade-off between sampling frequency and staleness in detecting the events. We argue that, for the problem at hand, using Age of Information (AoI) for quantifying the staleness of a sample is conservative and therefore, introduce \textit{age penalty} for this purpose. We study two optimization problems: minimize average age penalty subject to an average sampling frequency constraint, and minimize average sampling frequency subject to an average age penalty constraint; both are Constrained Markov Decision Problems. We solve them using linear programming approach and compute Markov policies that are optimal among all causal policies. Our numerical results demonstrate that the computed Markov policies not only outperform optimal periodic sampling policies, but also achieve sampling frequencies close to or lower than that of an optimal clairvoyant (non-causal) sampling policy, if a small age penalty is allowed.


title: Selective Weak Supervision for Neural Information Retrieval

abstract: This paper democratizes neural information retrieval to scenarios where large scale relevance training signals are not available. We revisit the classic IR intuition that anchor-document relations approximate query-document relevance and propose a reinforcement weak supervision selection method, ReInfoSelect, which learns to select anchor-document pairs that best weakly supervise the neural ranker (action), using the ranking performance on a handful of relevance labels as the reward. Iteratively, for a batch of anchor-document pairs, ReInfoSelect back propagates the gradients through the neural ranker, gathers its NDCG reward, and optimizes the data selection network using policy gradients, until the neural ranker's performance peaks on target relevance metrics (convergence). In our experiments on three TREC benchmarks, neural rankers trained by ReInfoSelect, with only publicly available anchor data, significantly outperform feature-based learning to rank methods and match the effectiveness of neural rankers trained with private commercial search logs. Our analyses show that ReInfoSelect effectively selects weak supervision signals based on the stage of the neural ranker training, and intuitively picks anchor-document pairs similar to query-document pairs.


title: Studies of Strange and Non-Strange Beauty Productions in PbPb Collisions with the CMS Detector

abstract: Beauty quarks are considered as one of the best probes of the strongly interacting medium created in relativistic heavy-ion collisions because they are predominantly produced via initial hard scatterings. Our measurements of B mesons productions provide information on the diffusion of beauty quarks and the flavor dependence of in-medium energy loss. In these studies, clarifying the hadronization mechanism is crucial for understanding the transport properties of beauty quarks. We propose to show our experimental results on three sets of particles: nonprompt $J/\psi$ and $D^0$, $B^0_s$ and $B^+$, and prompt $D^+_s$ and $D^0$. Our measurements of $B^0_s$ production can shed light on the mechanisms of beauty recombination in the medium. In addition, the measurements of mesons productions containing both strange and beauty quarks can provide more information about strangeness enhancement in the quark-gluon plasma. The ratio of $B^0_s$ over $B^+$ nuclear modification factors at a nucleon-nucleon center-of-mass energy of 5.02 TeV, using fully reconstructed B mesons with the CMS detector, are presented. We will also show the nonprompt $J/\psi$ and $D^0$ from the partial B mesons decay chains as well as our precise measurement the prompt $D^+_s$ to $D^0$ ratio in pp and PbPb collisions.


title: Change of measure under the hard-to-borrow model

abstract: As the Securities and Exchange Commission(SEC) has implemented a new regulation on short-sellings, short-sellers are required to repurchase stocks once the clearing risk rises to a certain level. Avellaneda and Lipkin proposed a fully coupled SDE system to describe the mechanism which is referred as Hard-To-Borrow(HTB) models. Guiyuan Ma obtained the PDE system for both American and European options. There is a technical error in Guiyuan Ma where two correlated Brownian motion should be converted before change of measure. In this paper, I will provide supplement conditions.


title: Impact of dijet and D-meson data from 5.02 TeV p+Pb collisions on nuclear PDFs

abstract: We discuss the new constraints on gluon parton distribution function (PDF) in lead nucleus, derivable with the Hessian PDF reweighting method from the 5.02 TeV p+Pb measurements of dijet (CMS) and $D^0$-meson (LHCb) nuclear modification ratios. The impact is found to be significant, placing stringent constraints in the mid- and previously unconstrained small-$x$ regions. The CMS dijet data confirm the existence of gluon anti-shadowing and the onset of small-$x$ shadowing, as well as reduce the gluon PDF uncertainties in the larger-$x$ region. The gluon constraints from the LHCb $D^0$ data, reaching down to $x \sim 10^{-5}$ and derived in a NLO perturbative QCD approach, provide a remarkable reduction of the small-$x$ uncertainties with a strong direct evidence of gluon shadowing. Furthermore, we find a good description of the data even down to zero $D^0$-meson transverse momentum within a purely DGLAP-based approach without a need for imposing any non-linear effects. Importantly, the constraints obtained from the dijet and $D^0$ data are mutually fully consistent, supporting the universality of nuclear PDFs in hard-scattering processes.


title: Taking Recoveries to Task: Recovery-Driven Development for Recipe-based Robot Tasks

abstract: Robot task execution when situated in real-world environments is fragile. As such, robot architectures must rely on robust error recovery, adding non-trivial complexity to highly-complex robot systems. To handle this complexity in development, we introduce Recovery-Driven Development (RDD), an iterative task scripting process that facilitates rapid task and recovery development by leveraging hierarchical specification, separation of nominal task and recovery development, and situated testing. We validate our approach with our challenge-winning mobile manipulator software architecture developed using RDD for the FetchIt! Challenge at the IEEE 2019 International Conference on Robotics and Automation. We attribute the success of our system to the level of robustness achieved using RDD, and conclude with lessons learned for developing such systems.


title: An operational information decomposition via synergistic disclosure

abstract: Multivariate information decompositions hold promise to yield insight into complex systems, and stand out for their ability to identify synergistic phenomena. However, the adoption of these approaches has been hindered by there being multiple possible decompositions, and no precise guidance for preferring one over the others. At the heart of this disagreement lies the absence of a clear operational interpretation of what synergistic information is. Here we fill this gap by proposing a new information decomposition based on a novel operationalisation of informational synergy, which leverages recent developments in the literature of data privacy. Our decomposition is defined for any number of information sources, and its atoms can be calculated using elementary optimisation techniques. The decomposition provides a natural coarse-graining that scales gracefully with the system's size, and is applicable in a wide range of scenarios of practical interest.


title: CSNNs: Unsupervised, Backpropagation-free Convolutional Neural Networks for Representation Learning

abstract: This work combines Convolutional Neural Networks (CNNs), clustering via Self-Organizing Maps (SOMs) and Hebbian Learning to propose the building blocks of Convolutional Self-Organizing Neural Networks (CSNNs), which learn representations in an unsupervised and Backpropagation-free manner. Our approach replaces the learning of traditional convolutional layers from CNNs with the competitive learning procedure of SOMs and simultaneously learns local masks between those layers with separate Hebbian-like learning rules to overcome the problem of disentangling factors of variation when filters are learned through clustering. We investigate the learned representation by designing two simple models with our building blocks, achieving comparable performance to many methods which use Backpropagation, while we reach comparable performance on Cifar10 and give baseline performances on Cifar100, Tiny ImageNet and a small subset of ImageNet for Backpropagation-free methods.


title: Eigen-Stratified Models

abstract: Stratified models depend in an arbitrary way on a selected categorical feature that takes $K$ values, and depend linearly on the other $n$ features. Laplacian regularization with respect to a graph on the feature values can greatly improve the performance of a stratified model, especially in the low-data regime. A significant issue with Laplacian-regularized stratified models is that the model is $K$ times the size of the base model, which can be quite large. We address this issue by formulating eigen-stratifed models, which are stratified models with an additional constraint that the model parameters are linear combinations of some modest number $m$ of bottom eigenvectors of the graph Laplacian, i.e., those associated with the $m$ smallest eigenvalues. With eigen-stratified models, we only need to store the $m$ bottom eigenvectors and the corresponding coefficients as the stratified model parameters. This leads to a reduction, sometimes large, of model size when $m \leq n$ and $m \ll K$. In some cases, the additional regularization implicit in eigen-stratified models can improve out-of-sample performance over standard Laplacian regularized stratified models.


title: Dise\~no de un controlador de \'angulo en un balanc\'in

abstract: This document describes the design of a PID controller for a rotation pitch plant of a degree freedom. In the controller design, the tuning methods of Astr\"on H\"agglund (AH), Kaiser Chaira (KC) and Kaiser Rajka (KR) will be used, verifying the performance in simulations and in the plant. Finally, the development for the implementation of an analog PID controller through circuits with operational amplifiers is described.


title: Low-rank matrix denoising for count data using unbiased Kullback-Leibler risk estimation

abstract: Many statistical studies are concerned with the analysis of observations organized in a matrix form whose elements are count data. When these observations are assumed to follow a Poisson or a multinomial distribution, it is of interest to focus on the estimation of either the intensity matrix (Poisson case) or the compositional matrix (multinomial case) when it is assumed to have a low rank structure. In this setting, it is proposed to construct an estimator minimizing the regularized negative log-likelihood by a nuclear norm penalty. Such an approach easily yields a low-rank matrix-valued estimator with positive entries which belongs to the set of row-stochastic matrices in the multinomial case. Then, as a main contribution, a data-driven procedure is constructed to select the regularization parameter in the construction of such estimators by minimizing (approximately) unbiased estimates of the Kullback-Leibler (KL) risk in such models, which generalize Stein's unbiased risk estimation originally proposed for Gaussian data. The evaluation of these quantities is a delicate problem, and novel methods are introduced to obtain accurate numerical approximation of such unbiased estimates. Simulated data are used to validate this way of selecting regularizing parameters for low-rank matrix estimation from count data. For data following a multinomial distribution, the performances of this approach are also compared to $K$-fold cross-validation. Examples from a survey study and metagenomics also illustrate the benefits of this methodology for real data analysis.


title: Commutators and images of noncommutative polynomials

abstract: Let $A$ be an algebra and let $f$ be a nonconstant noncommutative polynomial. In the first part of the paper, we consider the relationship between $[A,A]$, the linear span of commutators in $A$, and span$f(A)$, the linear span of the image of $f$ in $A$. In particular, we show that $[A,A]=A$ implies span$f(A)=A$. In the second part, we establish some Waring type results for images of polynomials. For example, we show that if $C$ is a commutative unital algebra over a field $F$ of characteristic $0$, $A$ is the matrix algebra $M_n(C)$, and the polynomial $f$ is neither an identity nor a central polynomial of $M_n(F)$, then every commutator in $A$ can be written as a difference of two elements, each of which is a sum of $7788$ elements from $f(A)$ (if $C=F$ is an algebraically closed field, then $4$ elements suffice). Similar results are obtained for some other algebras, in particular for the algebra $B(H)$ of all bounded linear operators on a Hilbert space $H$.


title: A random forest based approach for predicting spreads in the primary catastrophe bond market

abstract: We introduce a random forest approach to enable spreads' prediction in the primary catastrophe bond market. We investigate whether all information provided to investors in the offering circular prior to a new issuance is equally important in predicting its spread. The whole population of non-life catastrophe bonds issued from December 2009 to May 2018 is used. The random forest shows an impressive predictive power on unseen primary catastrophe bond data explaining 93% of the total variability. For comparison, linear regression, our benchmark model, has inferior predictive performance explaining only 47% of the total variability. All details provided in the offering circular are predictive of spread but in a varying degree. The stability of the results is studied. The usage of random forest can speed up investment decisions in the catastrophe bond industry.


title: Graph Neighborhood Attentive Pooling

abstract: Network representation learning (NRL) is a powerful technique for learning low-dimensional vector representation of high-dimensional and sparse graphs. Most studies explore the structure and metadata associated with the graph using random walks and employ an unsupervised or semi-supervised learning schemes. Learning in these methods is context-free, because only a single representation per node is learned. Recently studies have argued on the sufficiency of a single representation and proposed a context-sensitive approach that proved to be highly effective in applications such as link prediction and ranking. However, most of these methods rely on additional textual features that require RNNs or CNNs to capture high-level features or rely on a community detection algorithm to identify multiple contexts of a node. In this study, without requiring additional features nor a community detection algorithm, we propose a novel context-sensitive algorithm called GAP that learns to attend on different parts of a node's neighborhood using attentive pooling networks. We show the efficacy of GAP using three real-world datasets on link prediction and node clustering tasks and compare it against 10 popular and state-of-the-art (SOTA) baselines. GAP consistently outperforms them and achieves up to ~9% and ~20% gain over the best performing methods on link prediction and clustering tasks, respectively.


title: Direct observation of photonic Landau levels and helical edge states in strained honeycomb lattices

abstract: We report the realization of a synthetic magnetic field for photons and polaritons in a honeycomb lattice of coupled semiconductor micropillars. A strong synthetic field is induced in both the s and p orbital bands by engineering a uniaxial hopping gradient in the lattice, giving rise to the formation of Landau levels at the Dirac points. We provide direct evidence of the sublattice symmetry breaking of the lowest-order Landau level wavefunction, a distinctive feature of synthetic magnetic fields. Our realization implements helical edge states in the gap between n=0 and n=1 Landau levels, experimentally demonstrating a novel way of engineering propagating edge states in photonic lattices. In light of recent advances in the enhancement of polariton-polariton nonlinearities, the Landau levels reported here are promising for the study of the interplay between pseudomagnetism and interactions in a photonic system.


title: Bandit optimisation of functions in the Mat\'ern kernel RKHS

abstract: We consider the problem of optimising functions in the reproducing kernel Hilbert space (RKHS) of a Mat\'ern kernel with smoothness parameter $\nu$ over the domain $[0,1]^d$ under noisy bandit feedback. Our contribution, the $\pi$-GP-UCB algorithm, is the first practical approach with guaranteed sublinear regret for all $\nu>1$ and $d \geq 1$. Empirical validation suggests better performance and drastically improved computational scalablity compared with its predecessor, Improved GP-UCB.


title: A partially overdetermined problem in domains with partial umbilical boundary in space forms

abstract: In the first part of this paper, we consider a partially overdetermined mixed boundary value problem in space forms and generalize the main result in \cite{GX} into the case of general domains with partial umbilical boundary in space forms. Precisely, we prove that a partially overdetermined problem in a domain with partial umbilical boundary admits a solution if and only if the rest part of the boundary is also part of an umbilical hypersurface. In the second part of this paper, we prove a Heintze-Karcher-Ros type inequality for embedded hypersurfaces with free boundary lying on a horosphere or an equidistant hypersurface in the hyperbolic space. As an application, we show Alexandrov type theorem for constant mean curvature hypersurfaces with free boundary in these settings.


title: A Kernel Mean Embedding Approach to Reducing Conservativeness in Stochastic Programming and Control

abstract: We apply kernel mean embedding methods to sample-based stochastic optimization and control. Specifically, we use the reduced-set expansion method as a way to discard sampled scenarios. The effect of such constraint removal is improved optimality and decreased conservativeness. This is achieved by solving a distributional-distance-regularized optimization problem. We demonstrated this optimization formulation is well-motivated in theory, computationally tractable and effective in numerical algorithms.


title: QActor: On-line Active Learning for Noisy Labeled Stream Data

abstract: Noisy labeled data is more a norm than a rarity for self-generated content that is continuously published on the web and social media. Due to privacy concerns and governmental regulations, such a data stream can only be stored and used for learning purposes in a limited duration. To overcome the noise in this on-line scenario we propose QActor which novel combines: the selection of supposedly clean samples via quality models and actively querying an oracle for the most informative true labels. While the former can suffer from low data volumes of on-line scenarios, the latter is constrained by the availability and costs of human experts. QActor swiftly combines the merits of quality models for data filtering and oracle queries for cleaning the most informative data. The objective of QActor is to leverage the stringent oracle budget to robustly maximize the learning accuracy. QActor explores various strategies combining different query allocations and uncertainty measures. A central feature of QActor is to dynamically adjust the query limit according to the learning loss for each data batch. We extensively evaluate different image datasets fed into the classifier that can be standard machine learning (ML) models or deep neural networks (DNN) with noise label ratios ranging between 30% and 80%. Our results show that QActor can nearly match the optimal accuracy achieved using only clean data at the cost of at most an additional 6% of ground truth data from the oracle.


title: On the history of negative and complex numbers interpretation

abstract: The history of the development of the concept of complex numbers from the 16th to 19th centuries. The origin and refinement of the geometric and physical meaning of complex numbers, the emergence of vectoral analysis.


title: Black hole gluing in de Sitter space

abstract: We construct dynamical many-black-hole spacetimes with well-controlled asymptotic behavior as solutions of the Einstein vacuum equation with positive cosmological constant. We accomplish this by gluing Schwarzschild-de Sitter or Kerr-de Sitter black hole metrics into neighborhoods of points on the future conformal boundary of de Sitter space, under certain balance conditions on the black hole parameters. We give a self-contained treatment of solving the Einstein equation directly for the metric, given the scattering data we encounter at the future conformal boundary. The main step in the construction is the solution of a linear divergence equation for trace-free symmetric 2-tensors; this is closely related to Friedrich's analysis of scattering problems for the Einstein equation on asymptotically simple spacetimes.


title: Convergence of Update Aware Device Scheduling for Federated Learning at the Wireless Edge

abstract: We study federated learning (FL) at the wireless edge, where power-limited devices with local datasets collaboratively train a joint model with the help of a remote parameter server (PS). We assume that the devices are connected to the PS through a bandwidth-limited shared wireless channel. At each iteration of FL, a subset of the devices are scheduled to transmit their local model updates to the PS over orthogonal channel resources, while each participating device must compress its model update to accommodate to its link capacity. We design novel scheduling and resource allocation policies that decide on the subset of the devices to transmit at each round, and how the resources should be allocated among the participating devices, not only based on their channel conditions, but also on the significance of their local model updates. We then establish convergence of a wireless FL algorithm with device scheduling, where devices have limited capacity to convey their messages. The results of numerical experiments show that the proposed scheduling policy, based on both the channel conditions and the significance of the local model updates, provides a better long-term performance than scheduling policies based only on either of the two metrics individually. Furthermore, we observe that when the data is independent and identically distributed (i.i.d.) across devices, selecting a single device at each round provides the best performance, while when the data distribution is non-i.i.d., scheduling multiple devices at each round improves the performance. This observation is verified by the convergence result, which shows that the number of scheduled devices should increase for a less diverse and more biased data distribution.


title: Improper Gaussian Signaling for the $K$-user MIMO Interference Channels with Hardware Impairments

abstract: This paper investigates the performance of improper Gaussian signaling (IGS) for the $K$-user multiple-input, multiple-output (MIMO) interference channel (IC) with hardware impairments (HWI). HWI may arise due to imperfections in the devices like I/Q imbalance, phase noise, etc. With I/Q imbalance, the received signal is a widely linear transformation of the transmitted signal and noise. Thus, the effective noise at the receivers becomes improper, which means that its real and imaginary parts are correlated and/or have unequal powers. IGS can improve system performance with improper noise and/or improper interference. In this paper, we study the benefits of IGS for this scenario in terms of two performance metrics: achievable rate and energy efficiency (EE). We consider the rate region, the sum-rate, the EE region and the global EE optimization problems to fully evaluate the IGS performance. To solve these non-convex problems, we employ an optimization framework based on majorization-minimization algorithms, which allow us to obtain a stationary point of any optimization problem in which either the objective function and/or constraints are linear functions of rates. Our numerical results show that IGS can significantly improve the performance of the $K$-user MIMO IC with HWI and I/Q imbalance, where its benefits increase with the number of users, $K$, and the imbalance level, and decrease with the number of antennas.


title: The Impact of Metallicity on the Evolution of Rotation and Magnetic Activity of Sun-Like Stars

abstract: The rotation rates and magnetic activity of Sun-like and low-mass (< 1.4 Msun) main-sequence stars are known to decline with time, and there now exist several models for the evolution of rotation and activity. However, the role that chemical composition plays during stellar spin-down has not yet been explored. In this work, we use a structural evolution code to compute the rotational evolution of stars with three different masses (0.7, 1.0, and 1.3Msun and six different metallicities, ranging from [Fe/H]=-1.0 to [Fe/H]=+0.5. We also implement three different wind-braking formulations from the literature (two modern and one classical) and compare their predictions for rotational evolution. The effect that metallicity has on stellar structural properties, and in particular the convective turnover timescale, leads the two modern wind-braking formulations to predict a strong dependence of the torque on metallicity. Consequently, they predict that metal rich stars spin-down more effectively at late ages (> 1 Gyr) than metal poor stars, and the effect is large enough to be detectable with current observing facilities. For example, the formulations predict that a Sun-like (solar-mass and solar-aged) star with [Fe/H]=-0.3 will have a rotation period of less than 20 days. Even though old, metal poor stars are predicted to rotate more rapidly at a given age, they have larger Rossby numbers and are thus expected to have lower magnetic activity levels. Finally, the different wind-braking formulations predict quantitative differences in the metallicity-dependence of stellar rotation, which may be used to test them.


title: Language-Based Web Session Integrity

abstract: Session management is a fundamental component of web applications: despite the apparent simplicity, correctly implementing web sessions is extremely tricky, as witnessed by the large number of existing attacks. This motivated the design of formal methods to rigorously reason about web session security which, however, are not supported at present by suitable automated verification techniques. In this paper we introduce the first security type system that enforces session security on a core model of web applications, focusing in particular on server-side code. We showcase the expressiveness of our type system by analyzing the session management logic of HotCRP, Moodle, and phpMyAdmin, unveiling novel security flaws that have been acknowledged by software developers.


title: Splitting methods and short time existence for the master equations in mean field games

abstract: We develop a splitting method to prove the well-posedness, in short time, of solutions for two master equations in mean field game (MFG) theory: the second order master equation, describing MFGs with a common noise, and the system of master equations associated with MFGs with a major player. Both problems are infinite dimensional equations stated in the space of probability measures. Our new approach simplifies, shortens and generalizes previous existence results for second order master equations and provides the first existence result for systems associated with MFG problems with a major player.


title: On Polynomials in Primes, Ergodic Averages and Monothetic Groups

abstract: Let $G$ denote a compact monothetic group, and let $$\rho (x) = \alpha_k x^k + \ldots + \alpha_1 x + \alpha_0,$$ where $\alpha_0, \ldots , \alpha_k$ are elements of $G$ one of which is a generator of $G$. Let $(p_n)_{n\geq 1}$ denote the sequence of rational prime numbers. Suppose $f \in L^{p}(G)$ for $p> 1$. It is known that if $$A_{N}f(x) := {1 \over N} \sum_{n=1}^{N} f(x + \rho (p_n)) \qquad (N=1,2, \ldots ),$$ then the limit $\lim _{n\to \infty} A_Nf(x)$ exists for almost all $x$ with respect Haar measure. We show that if $G$ is connected then the limit is $\int_{G} f d\lambda$. In the case where $G$ is the $a$-adic integers, which is a totally disconnected group, the limit is described in terms of Fourier multipliers which are generalizations of Gauss sums.


title: Second cohomology group of the finite-dimensional simple Jordan superalgebra $\mathcal{D}_{t}$, $t\neq 0$

abstract: The second cohomology group (SCG) of the Jordan superalgebra $\mathcal{D}_{t}$, $t\neq 0$, is calculated by using the coefficients which appear in the regular superbimodule $\mathrm{Reg}\mathcal{D}_t$. Contrary to the case of algebras, this group is nontrivial thanks to the non-splitting caused by the Wedderburn Decomposition Theorem \cite{Faber1}. First, to calculate the SCG of a Jordan superalgebra we use split-null extension of the Jordan superalgebra and the Jordan superalgebra representation. We prove conditions that satisfy the bilinear forms $h$ that determine the SCG in Jordan superalgebras. We use these to calculate the SCG for the Jordan superalgebra $\mathcal{D}_{t}$ , $t\neq 0$. Finally, we prove that $\mathcal{H}^2(\mathcal{D}_{t}, \textrm{Reg}\mathcal{D}_{t})=0\oplus\mathbb{F}^2$, $t\neq 0$.


title: Magnetic oscillation modes in square lattice artificial spin ice

abstract: Small amplitude dipolar oscillations are considered in artificial spin ice on a square lattice in two dimensions. The net magnetic moment of each elongated magnetic island in the spin ice is assumed to have Heisenberg-like dynamics. Each island's magnetic moment is assumed to be influenced by shape anisotropies and by the dipolar interactions with its nearest neighbors. The magnetic dynamics is linearized around one of the ground states, leading to an $8\times 8$ matrix to be diagonalized for the magnetic spin wave modes. Analytic solutions are found and classified as antisymmetric and symmetric with regard to their in-plane dynamic fluctuations. Although only the leading dipolar interactions are included, modes similar to these may be observable experimentally.


title: A Self-Matched Leaky-Wave Antenna for Ultrahigh-Field MRI with Low SAR

abstract: The technology of magnetic resonance imaging is developing towards higher magnetic fields to improve resolution and contrast. However, whole-body imaging at 7 T or even higher fields remains challenging due to wave interference, tissue inhomogneities and high RF power deposition. Nowadays, proper RF excitation of a human body in prostate and cardiac MRI is only possible to achieve by using phased arrays of antennas attached to the body (so-called surface coils). Due to safety concerns, the design of such coils aims to minimize the local specific absorption rate (SAR) keeping the highest possible RF signal in the region of interest. All previously demonstrated approaches were based on resonant structures such as e. g.dipoles, capacitively-loaded loops, TEM-line sections. In this study, we show that there is a better compromise between the transmit signal and the local SAR using non-resonant surface coils due to weaker RF near fields in the close proximity of their conductors. With this aim, we propose and experimentally demonstrate a first leaky-wave surface coil implemented as a periodically-slotted microstrip transmission line. Due to its non-resonant radiation, the proposed coil induces only slightly over half the peak local SAR compared to a state-of-the-art dipole coil, but has the same transmit efficiency in prostate imaging at 7 T. Unlike other coils, the leaky-wave coil intrinsically matches its input impedance to the averaged wave impedance of body tissues in a broad frequency range, which makes it very attractive for future clinical applications of 7 T MRI.


title: Global Sensitivity Analysis on Numerical Solver Parameters of Particle-In-Cell Models in Particle Accelerator Systems

abstract: Every computer model depends on numerical input parameters that are chosen according to mostly conservative but rigorous numerical or empirical estimates. These parameters could for example be the step size for time integrators, a seed for pseudo-random number generators, a threshold or the number of grid points to discretize a computational domain. In case a numerical model is enhanced with new algorithms and modelling techniques the numerical influence on the quantities of interest, the running time as well as the accuracy is often initially unknown. Usually parameters are chosen on a trial-and-error basis neglecting the computational cost versus accuracy aspects. As a consequence the cost per simulation might be unnecessarily high which wastes computing resources. Hence, it is essential to identify the most critical numerical parameters and to analyze systematically their effect on the result in order to minimize the time-to-solution without losing significantly on accuracy. Relevant parameters are identified by global sensitivity studies where Sobol' indices are common measures. These sensitivities are obtained by surrogate models based on polynomial chaos expansion. In this paper, we first introduce the general methods for uncertainty quantification. We then demonstrate their use on numerical solver parameters to reduce the computational costs and discuss further model improvements based on the sensitivity analysis. The sensitivities are evaluated for neighbouring bunch simulations of the existing PSI Injector II and PSI Ring as well as the proposed Daedalus Injector cyclotron and simulations of the rf electron gun of the Argonne Wakefield Accelerator.


title: Local Existence of Analytic Sharp Fronts for Singular SQG

abstract: In this paper, we prove local existence and uniqueness of analytic sharp-front solutions to a generalised SQG equation by the use of an abstract Cauchy--Kowalevskaya theorem. Here, the velocity is determined by $u = |\nabla|^{-2\beta}\nabla^\perp\theta $ which (for $1<\beta\leq 2$) is more singular than in SQG. This is achieved despite the appearance of pseudodifferential operators of order higher than one in our equation, by recasting our equation in a suitable integral form. We also provide a full proof of the abstract version of the Cauchy--Kowalevskaya theorem we use.


title: On the density of sumsets

abstract: Recently introduced by the authors in [Proc. Edinb. Math. Soc. 60 (2020), 139-167], quasi-densities form a large family of real-valued functions partially defined on the power set of the integers that serve as a unifying framework for the study of many known densities (including the asymptotic density, the Banach density, the logarithmic density, the analytic density, and the P\'olya density). We further contribute to this line of research by proving that (i) for each $n \in \mathbf N^+$ and $\alpha \in [0,1]$, there is $A \subseteq \mathbf{N}$ with $kA \in \text{dom}(\mu)$ and $\mu(kA) = \alpha k/n$ for every quasi-density $\mu$ and every $k=1,\ldots, n$, where $kA:=A+\cdots+A$ is the $k$-fold sumset of $A$ and $\text{dom}(\mu)$ denotes the domain of definition of $\mu$; (ii) for each $\alpha \in [0,1]$ and every non-empty finite $B\subseteq \mathbf{N}$, there is $A \subseteq \mathbf{N}$ with $A+B \in \mathrm{dom}(\mu)$ and $\mu(A+B)=\alpha$ for every quasi-density $\mu$; (iii) for each $\alpha \in [0,1]$, there exists $A\subseteq \mathbf{N}$ with $2A = \mathbf{N}$ such that $A \in \text{dom}(\mu)$ and $\mu(A) = \alpha$ for every quasi-density $\mu$. Proofs rely on the properties of a little known density first considered by R.C. Buck and the "structure" of the set of all quasi-densities; in particular, they are rather different than previously known proofs of special cases of the same results.


title: Predicting the Impact of Electric Field Stimulation in a Detailed Computational Model of Cortical Tissue

abstract: Neurostimulation using weak electric fields has generated excitement in recent years due to its potential as a medical intervention. However, study of this stimulation modality has been hampered by inconsistent results and large variability within and between studies. In order to begin addressing this variability we need to properly characterise the impact of the current on the underlying neuron populations. To develop and test a computational model capable of capturing the impact of electric field stimulation on networks of neurons. We construct a cortical tissue model with distinct layers and explicit neuron morphologies. We then apply a model of electrical stimulation and carry out multiple test case simulations. The cortical slice model is compared to experimental literature and shown to capture the main features of the electrophysiological response to stimulation. Namely, the model showed 1) a similar level of depolarisation in individual pyramidal neurons, 2) acceleration of intrinsic oscillations, and 3) retention of the spatial profile of oscillations in different layers. We then apply alternative electric fields to demonstrate how the model can capture differences in neuronal responses to the electric field. We demonstrate that the tissue response is dependent on layer depth, the angle of the apical dendrite relative to the field, and stimulation strength. We present publicly available computational modelling software that predicts the neuron network population response to electric field stimulation.


title: On variation functions and their moduli of continuity

abstract: We study the moduli of continuity of functions of bounded variation and of their variation functions. It is easy to see that the modulus of continuity of a function of bounded variation is always smaller or equal to the modulus of continuity of its variation function. We show that we cannot make any reasonable conclusion on the modulus of continuity of the variation function if we only know the modulus of continuity of the parent function itself. In particular, given two moduli of continuity, the first being weaker than Lipschitz continuity, we show that there exists a function of bounded variation with minimal modulus of continuity less than the first modulus of continuity, but with a variation function with minimal modulus of continuity greater than the second modulus of continuity. In particular, this negatively resolves the open problem whether the variation function of an $\alpha$-H\"older continuous function is $\alpha$-H\"older continuous.


title: Science Under Attack

abstract: Some attack scientific rationality, others defend it, but both miss the point. What both parties take to be scientific rationality is actually a species of irrationality masquerading as scientific rationality. The current orthodox conception of science, taken for granted by scientists and non-scientists alike, is irrational because it suppresses problematic assumptions, inherent in the aims of science, having to do with metaphysics, values, and political and social issues. We urgently need a more rigorous conception of science to be adopted and implemented that honestly acknowledges the problematic character of the aims of science, and seeks to improve them as science proceeds.


title: Power spectrum response of large-scale structure in 1D and in 3D: tests of prescriptions for post-collapse dynamics

abstract: The power spectrum response function of the large-scale structure of the Universe describes how the evolved power spectrum is modified by a small change in initial power through non-linear mode coupling of gravitational evolution. It was previously found that the response function for the coupling from small to large scales is strongly suppressed in amplitude, especially at late times, compared to predictions from perturbation theory (PT) based on the single-stream approximation. One obvious explanation for this is that PT fails to describe the dynamics beyond shell-crossing. We test this idea by comparing measurements in $N$-body simulations to prescriptions based on PT but augmented with adaptive smoothing to account for the formation of non-linear structures of various sizes in the multi-stream regime. We first start with one-dimensional (1D) cosmology, where the Zel'dovich approximation provides the exact solution in the single stream regime. Similarly to the 3D case, the response function of the large-scale modes exhibits a strong suppression in amplitude at small scales which cannot be explained by the Zel'dovich solution alone. However, by performing adaptive smoothing of initial conditions to identify haloes of different sizes and solving approximately post-collapse dynamics in the 3-streams regime, agreement between theory and simulations drastically improves. We extend our analyses to the 3D case using PINOCCHIO algorithm, in which similar adaptive smoothing is implemented on the Lagrangian PT fields to identify haloes and is combined with a spherical halo prescription to account for post-collapse dynamics. Again, a suppression is found in the coupling between small- and large-scale modes and the agreement with simulations is improved.


title: Counter-propagating photon pair generation in a nonlinear waveguide

abstract: Counter-propagating parametric conversion processes in non-linear bulk crystals have been shown to feature unique properties for efficient narrowband frequency conversion. In quantum optics, the generation of photon pairs with a counter-propagating parametric down-conversion process (PDC) in a waveguide, where signal and idler photons propagate in opposite directions, offers unique material-independent engineering capabilities. However, realizing counter-propagating PDC necessitates quasi-phase-matching (QPM) with extremely short poling periods. Here, we report on the generation of counter-propagating single-photon pairs in a self-made periodically poled lithium niobate waveguide with a poling period on the same order of magnitude as the generated wavelength. The single photons of the biphoton state bridge GHz and THz bandwidths with a separable joint temporal-spectral behavior. Furthermore, they allow the direct observation of the temporal envelope of heralded single photons with state-of-the art photon counters.


title: Structure theory of p.p. rings and their generalizations

abstract: In this paper, new and significant advances on the understanding the structure of p.p. rings and their generalizations have been made. Especially among them, it is proved that a commutative ring $R$ is a generalized p.p. ring if and only if $R$ is a generalized p.f. ring and its minimal spectrum is Zariski compact, or equivalently, $R/\mathfrak{N}$ is a p.p. ring and $R_{\mathfrak{m}}$ is a primary ring for all $\mathfrak{m}\in\rm{Max}(R)$. Some of the major results of the literature either are improved or are proven by new methods. In particular, we give a new and quite elementary proof to the fact that a commutative ring $R$ is a p.p. ring if and only if $R[x]$ is a p.p. ring.


title: OPFython: A Python-Inspired Optimum-Path Forest Classifier

abstract: Machine learning techniques have been paramount throughout the last years, being applied in a wide range of tasks, such as classification, object recognition, person identification, and image segmentation. Nevertheless, conventional classification algorithms, e.g., Logistic Regression, Decision Trees, and Bayesian classifiers, might lack complexity and diversity, not suitable when dealing with real-world data. A recent graph-inspired classifier, known as the Optimum-Path Forest, has proven to be a state-of-the-art technique, comparable to Support Vector Machines and even surpassing it in some tasks. This paper proposes a Python-based Optimum-Path Forest framework, denoted as OPFython, where all of its functions and classes are based upon the original C language implementation. Additionally, as OPFython is a Python-based library, it provides a more friendly environment and a faster prototyping workspace than the C language.


title: Nonlocal PageRank

abstract: In this work we introduce and study a nonlocal version of the PageRank. In our approach, the random walker explores the graph using longer excursions than just moving between neighboring nodes. As a result, the corresponding ranking of the nodes, which takes into account a \textit{long-range interaction} between them, does not exhibit concentration phenomena typical of spectral rankings which take into account just local interactions. We show that the predictive value of the rankings obtained using our proposals is considerably improved on different real world problems.


title: NAS-Bench-1Shot1: Benchmarking and Dissecting One-shot Neural Architecture Search

abstract: One-shot neural architecture search (NAS) has played a crucial role in making NAS methods computationally feasible in practice. Nevertheless, there is still a lack of understanding on how these weight-sharing algorithms exactly work due to the many factors controlling the dynamics of the process. In order to allow a scientific study of these components, we introduce a general framework for one-shot NAS that can be instantiated to many recently-introduced variants and introduce a general benchmarking framework that draws on the recent large-scale tabular benchmark NAS-Bench-101 for cheap anytime evaluations of one-shot NAS methods. To showcase the framework, we compare several state-of-the-art one-shot NAS methods, examine how sensitive they are to their hyperparameters and how they can be improved by tuning their hyperparameters, and compare their performance to that of blackbox optimizers for NAS-Bench-101.


title: Rate of Estimation for the Stationary Distribution of Stochastic Damping Hamiltonian Systems with Continuous Observations

abstract: We study the problem of the non-parametric estimation for the density $\pi$ of the stationary distribution of a stochastic two-dimensional damping Hamiltonian system $(Z_t)_{t\in[0,T]}=(X_t,Y_t)_{t \in [0,T]}$. From the continuous observation of the sampling path on $[0,T]$, we study the rate of estimation for $\pi(x_0,y_0)$ as $T \to \infty$. We show that kernel based estimators can achieve the rate $T^{-v}$ for some explicit exponent $v \in (0,1/2)$. One finding is that the rate of estimation depends on the smoothness of $\pi$ and is completely different with the rate appearing in the standard i.i.d.\ setting or in the case of two-dimensional non degenerate diffusion processes. Especially, this rate depends also on $y_0$. Moreover, we obtain a minimax lower bound on the $L^2$-risk for pointwise estimation, with the same rate $T^{-v}$, up to $\log(T)$ terms.


title: Parallel solution of saddle point systems with nested iterative solvers based on the Golub-Kahan Bidiagonalization

abstract: We present a scalability study of Golub-Kahan bidiagonalization for the parallel iterative solution of symmetric indefinite linear systems with a 2x2 block structure. The algorithms have been implemented within the parallel numerical library PETSc. Since a nested inner-outer iteration strategy may be necessary, we investigate different choices for the inner solvers, including parallel sparse direct and multigrid accelerated iterative methods. We show the strong and weak scalability of the Golub-Kahan bidiagonalization based iterative method when applied to a two-dimensional Poiseuille flow and to two- and three-dimensional Stokes test problems.


title: Purity in chromatically localized algebraic $K$-theory

abstract: We prove a purity property in telescopically localized algebraic $K$-theory of ring spectra: For $n\geq 1$, the $T(n)$-localization of $K(R)$ only depends on the $T(0)\oplus \dots \oplus T(n)$-localization of $R$. This complements a classical result of Waldhausen in rational $K$-theory. Combining our result with work of Clausen--Mathew--Naumann--Noel, one finds that $L_{T(n)}K(R)$ in fact only depends on the $T(n-1)\oplus T(n)$-localization of $R$, again for $n \geq 1$. As consequences, we deduce several vanishing results for telescopically localized $K$-theory, as well as an equivalence between $K(R)$ and $\mathrm{TC}(\tau_{\geq 0} R)$ after $T(n)$-localization for $n\geq 2$.


title: Multi-trek separation in Linear Structural Equation Models

abstract: Building on the theory of causal discovery from observational data, we study interactions between multiple (sets of) random variables in a linear structural equation model with non-Gaussian error terms. We give a correspondence between structure in the higher order cumulants and combinatorial structure in the causal graph. It has previously been shown that low rank of the covariance matrix corresponds to trek separation in the graph. Generalizing this criterion to multiple sets of vertices, we characterize when determinants of subtensors of the higher order cumulant tensors vanish. This criterion applies when hidden variables are present as well. For instance, it allows us to identify the presence of a hidden common cause of k of the observed variables.


title: 30$^\circ$-twisted bilayer graphene quasicrystals from chemical vapor deposition

abstract: The artificial stacking of atomically thin crystals suffers from intrinsic limitations in terms of control and reproducibility of the relative orientation of exfoliated flakes. This drawback is particularly severe when the properties of the system critically depend on the twist angle, as in the case of the dodecagonal quasicrystal formed by two graphene layers rotated by 30$^\circ$. Here we show that large-area 30$^\circ$-rotated bilayer graphene can be grown deterministically by chemical vapor deposition on Cu, eliminating the need of artificial assembly. The quasicrystals are easily transferred to arbitrary substrates and integrated in high-quality hBN-encapsulated heterostructures, which we process into dual-gated devices exhibiting carrier mobility up to $10^5$ cm$^2$/Vs. From low-temperature magnetotransport, we find that the graphene quasicrystals effectively behave as uncoupled graphene layers, showing 8-fold degenerate quantum Hall states: this result indicates that the Dirac cones replica detected by previous photo-emission experiments do not contribute to the electrical transport.


title: Sum rule improved double parton distributions in position space

abstract: Models for double parton distributions that are realistic and consistent with theoretical constraints are crucial for a reliable description of double parton scattering. We show how an ansatz that has the correct behaviour in the limit of small transverse distance between the partons can be improved step by step, such as to fulfil the sum rules for double parton distributions with an accuracy around 10%.


title: Singularity-Free Inverse Dynamics for Underactuated Systems with a Rotating Mass

abstract: Motion control of underactuated systems through the inverse dynamics contains configuration singularities. These limitations in configuration space mainly stem from the inertial coupling that passive joints/bodies create. In this study, we present a model that is free from singularity while the trajectory of the rotating mass has a small-amplitude sine wave around its circle. First, we derive the modified non-linear dynamics for a rolling system. Also, the singularity regions for this underactuated system is demonstrated. Then, the wave parameters are designed under certain conditions to remove the coupling singularities. We obtain these conditions from the positive definiteness of the inertia matrix in the inverse dynamics. Finally, the simulation results are confirmed by using a prescribed Beta function on the specified states of the rolling carrier. Because our algebraic method is integrated into the non-linear dynamics, the proposed solution has a great potential to be extended to the Lagrangian mechanics with multiple degrees-of-freedom.


title: Electrically reconfigurable silicon photonic filter enabled by embedded phase change material in microring resonator

abstract: We report a tunable optical filter based on phase change material $Ge_{2}Sb_{2}Te_{5}$ embedded in a silicon microring resonator. The high thermo-optic coefficient of $Ge_{2}Sb_{2}Te_{5}$ in amorphous phase enables tuning of resonance wavelength in broad range with a very small active volume. The low-loss indium-tin-oxide electrodes are employed to induce Joule heating in $Ge_{2}Sb_{2}Te_{5}$-Si active waveguide region. The electrically induced heating in the active region alters the effective refractive index of hybrid microring resulting in a wavelength tuning of 1.04 nm for an applied voltage of only 3V. The device exhibits high extinction ratios in the range of 20-41 dB and a compact active footprint 0.96 $\mu m^{2}$ only and is suitable for the large scale reconfigurable integrated photonic circuits.


title: Weight truncation for wall-crossings in birational cobordisms

abstract: We develop the technique of weight truncation in the context of wall-crossings in birational cobordisms, parallel to that in [HL15, BFK19]. More precisely, for each such wall-crossing, we embed the bounded above derived category of coherent sheaves of the semistable part as a semi-orthogonal summand of that of the stack in question. Our construction does not require any smoothness assumptions, and exhibits a strong symmetry across the two sides of the wall-crossing. As an application, we show that for wall-crossings satisfying suitable regularity conditions, a certain duality of local cohomology complexes implies the existence of a fully faithful functor/equivalence between the derived categories under wall-crossings.


title: Investment behavior and firms' financial performance: A comparative analysis using firm-level data from the wine industry

abstract: This paper assesses the role of financial performance in explaining firms' investment dynamics in the wine industry from the three European Union (EU) largest producers. The wine sector deserves special attention to investigate firms' investment behavior given the high competition imposed by the latecomers. More precisely, we investigate how the capitalization, liquidity and profitability influence the investment dynamics using firm-level data from the wine industry from France (331 firms), Italy (335) firms and Spain (442) firms. We use data from 2007 to 2014, drawing a comparison between these countries, and relying on difference-and system-GMM estimators. Specifically, the impact of profitability is positive and significant, while the capitalization has a significant and negative impact on the investment dynamics only in France and Spain. The influence of the liquidity ratio is negative and significant only in the case of Spain. Therefore, we notice different investment strategies for wine companies located in the largest producer countries. It appears that these findings are in general robust to different specifications of liquidity and profitability ratios, and to the different estimators we use.


title: Many-body perturbation theories for finite nuclei

abstract: In recent years many-body perturbation theory encountered a renaissance in the field of ab initio nuclear structure theory. In various applications it was shown that perturbation theory, including novel flavors of it, constitutes a useful tool to describe atomic nuclei, either as a full-fledged many-body approach or as an auxiliary method to support more sophisticated non-perturbative many-body schemes. In this work the current status of many-body perturbation theory in the field of nuclear structure is discussed and novel results are provided that highlight its power as a efficient and yet accurate (pre-processing) approach to systematically investigate medium-mass nuclei. Eventually a new generation of chiral nuclear Hamiltonians is benchmarked using several state-of-the-art flavours of many-body perturbation theory.


title: Equation of State and Progenitor Dependence of Stellar-Mass Black-Hole Formation

abstract: The core collapse of a massive star results in the formation of a proto-neutron star (PNS). If enough material is accreted onto a PNS it will become gravitationally unstable and further collapse into a black-hole (BH). We perform a systematic study of failing core-collapse supernovae in spherical symmetry for a wide range of presupernova progenitor stars and equations of state (EOSs) of nuclear matter. We analyze how variations in progenitor structure and the EOS of dense matter above nuclear saturation density affect the PNS evolution and subsequent BH formation. Comparisons of core-collapse for a given progenitor star and different EOSs show that the path traced by the PNS in mass-entropy phase space $M_{\mathrm{grav}}^{\mathrm{PNS}}-\tilde{s}$ is well correlated with the progenitor compactness and almost EOS independent, apart from the final endpoint. Furthermore, BH formation occurs, to a very good approximation, soon after the PNS overcomes the maximum \textit{gravitational} mass supported by a hot NS with constant entropy equal to $\tilde{s}$. These results show a path to constraining the temperature dependence of the EOS through the detection of neutrinos from a failed galactic supernova.


title: Bases of Quantum Group Algebras in Terms of Lyndon Words

abstract: We have reviewed some results on quantized shuffling, and in particular, the grading and structure of this algebra. In parallel, we have summarized certain details about classical shuffle algebras, including Lyndon words (primes) and the construction of bases of classical shuffle algebras in terms of Lyndon words. We have explained how to adapt this theory to the construction of bases of quantum group algebras in terms of Lyndon words. This method has a limited application to the specific case of the quantum group parameter being a root of unity, with the requirement that specialization to the root of unity is non-restricted. As an additional, applied part of this work, we have implemented a Wolfram Mathematica package with functions for quantum shuffle multiplication and constructions of bases in terms of Lyndon words.


title: Characterisation of the pressure term in the incompressible Navier-Stokes equations on the whole space

abstract: We characterise the pressure term in the incompressible 2D and 3D Navier-Stokes equations for solutions defined on the whole space.


title: Differentiation of neural-type cells on multi-scale ordered collagen-silica bionanocomposites

abstract: Cells respond to biophysical and biochemical signals. We developed a composite filament from collagen and silica particles modified to interact with collagen and/or present a laminin epitope (IKVAV) crucial for cell-matrix adhesion and signal transduction. This combines scaffolding and signaling and shows that local tuning of collagen organization enhances cell differentiation.


title: Two-dimensional patterns in dip coating -- first steps on the continuation path

abstract: We present a brief comparative investigation of the bifurcation structure related to the formation of two-dimensional deposition patterns as described by continuum models of Cahn-Hilliard type. These are, on the one hand a driven Cahn-Hilliard model for Langmuir-Blodgett transfer of a surfactant layer from the surface of a bath onto a moving plate and on the other hand a driven thin-film equation modelling the surface acoustic wave-driven coating of a plate by a simple liquid. In both cases, we present selected two-dimensional steady states corresponding to deposition patterns and discuss the main structure of the corresponding bifurcation diagrams.


title: How Many Simulations Do We Exist In? A Practical Mathematical Solution to the Simulation Argument

abstract: The Simulation Argument has gained significant traction in the public arena. It has offered a hypothesis based on probabilistic analysis of its assumptions that we are likely to exist within a computer simulation. This has been derived from factors including the prediction of computing power, human existence, extinction and population dynamics, and suggests a very large value for the number of possible simulations within which we may exist. On evaluating this argument through the application of tangible real-world evidence and projections, it is possible to calculate real numerical solutions for the Simulation Argument. This reveals a much smaller number of possible simulations within which we may exist, and offers a novel practicable approach in which to appraise the variety and multitude of conjectures and theories associated with the Simulation Hypothesis.


title: Analyzing Factors Associated with Fatal Road Crashes: A Machine Learning Approach

abstract: Road traffic injury accounts for a substantial human and economic burden globally. Understanding risk factors contributing to fatal injuries is of paramount importance. In this study, we proposed a model that adopts a hybrid ensemble machine learning classifier structured from sequential minimal optimization and decision trees to identify risk factors contributing to fatal road injuries. The model was constructed, trained, tested, and validated using the Lebanese Road Accidents Platform (LRAP) database of 8482 road crash incidents, with fatality occurrence as the outcome variable. A sensitivity analysis was conducted to examine the influence of multiple factors on fatality occurrence. Seven out of the nine selected independent variables were significantly associated with fatality occurrence, namely, crash type, injury severity, spatial cluster-ID, and crash time (hour). Evidence gained from the model data analysis will be adopted by policymakers and key stakeholders to gain insights into major contributing factors associated with fatal road crashes and to translate knowledge into safety programs and enhanced road policies.


title: Orthant-Strictly Monotonic Norms, Generalized Top-k and k-Support Norms and the L0 Pseudonorm

abstract: The so-called l0 pseudonorm on the Euclidean space Rd counts the number of nonzero components of a vector. We say that a sequence of norms is strictly increasingly graded (with respect to the l0 pseudonorm) if it is nondecreasing and that the sequence of norms of a vector~x becomes stationary exactly at the index l0(x). In this paper, with any (source) norm, we associate sequences of generalized top-k and k-support norms, and we also introduce the new class of orthant-strictly monotonic norms (that encompasses the lp norms, but for the extreme ones). Then, we show that an orthant-strictly monotonic source norm generates a sequence of generalized top-k norms which is strictly increasingly graded. With this, we provide a systematic way to generate sequences of norms with which the level sets of the l0 pseudonorm are expressed by means of the difference of two norms. Our results rely on the study of orthant-strictly monotonic norms.


title: Four Points and a Quadric

abstract: This paper provides a new simple proof of Hesse's theorem in projective geometry for any dimension.


title: Why the same key can open different padlocks?

abstract: We accidentally observed that one key can open two padlocks from different brands. Based on this observation, here we derived equations to estimate the number of different padlocks that can be produced. We focus on simple padlocks that are usually used in most households. It is made of a rectangular block of metal and crafted at any position. Different keys have different sizes and positions of the crafting (hole) and un-crafting (hill) portions. We limit the study to only padlocks having a maximum of four holes. This is the typical architecture of household padlocks. We observed that the number of different padlocks depends on the minimum size of the holes and hills. We also observe a scaling relationship between the number of different padlocks that can be produced with the size of holes or hills.


title: Description of the set of strictly regular quadratic bistochastic operators and examples

abstract: The present paper focuses on the dynamical systems of the quadratic bistochastic operators (QBO) on the standard simplex. In the paper, we show the character of connection of the dynamical systems of a bistochastic operator with the dynamical systems of the extreme bistochastic operators. In addition, we prove that almost all quadratic bistochastic operators are strictly regular and give a description of the strictly regular quadratic bistochastic operators in the convex polytope QBO. Furthermore, the density of the set of strictly regular QBO in the set of QBO is proved and nontrivial examples of strictly regular bistochastic operators are given.


title: Prime Role of the Curvature of Potentials in Physics; Application to Inertia

abstract: The curvature of the inertial or gravitational potentials defined as a Hodge-Helmholtz decomposition of acceleration into an irrotational and a solenoidal components, enable to federate certain domains of macroscopic physics. After two verifications in physics, one on the calculation of the curvatures for capillary effects and the second on the deflection of light by a gravitational effect, the concept of curvature of the potential is applied to the inertia. The physical analysis of each of the contributions of the decomposition is carried out on a classic example of fluid mechanics, the backward-facing step flow where inertia plays a preponderant role on the recirculation length.


title: Studies on tetrafluoropropene-based gas mixtures with low environmental impact for Resistive Plate Chambers

abstract: Gaseous detectors are widely used in high-energy physics experiments, and in particular at the CERN Large Hadron Collider (LHC), to provide tracking and triggering over large volumes. It has been recently estimated that Resistive Plate Chambers (RPC), used for muon detection, have the highest contribution on the overall greenhouse gas (GHG) emissions at the LHC experiments. Gas mixtures for RPCs are mainly made of C$_{2}$H$_{2}$F$_{4}$, which is a greenhouse gas with a high environmental impact in the atmosphere. C$_{2}$H$_{2}$F$_{4}$ is already phasing out of production, due to recent European Union (EU) regulations, and its cost is expected to increase in the near future. Therefore, finding alternative gas mixtures made of gas components with a low Global Warming Potential (GWP) has become extremely important for limiting the GHG emissions as well as for economic reasons. The novel hydrofluoroolefins are likely appropriate candidates to replace C$_{2}$H$_{2}$F$_{4}$ due to their similar chemical structures. This study is focused on the characterization of innovative gas mixtures with tetrafluoropropene HFO1234ze(E) (C$_{3}$H$_{2}$F$_{4}$) that is one of the most eco-friendly hydrofluoroolefins, thanks to its very low GWP. HFO1234ze(E)-based gas mixtures with the addition of Ar, N$_{2}$, O$_{2}$ and CO$_{2}$ are extensively discussed in this paper as well as the role of i-C$_{4}$H$_{10}$ and SF$_{6}$ in such mixtures.


title: Nonexistence of subcritical solitary waves

abstract: We prove the nonexistence of two-dimensional solitary gravity water waves with subcritical wave speeds and an arbitrary distribution of vorticity. This is a longstanding open problem, and even in the irrotational case there are only partial results relying on sign conditions or smallness assumptions. As a corollary, we obtain a relatively complete classification of solitary waves: they must be supercritical, symmetric, and monotonically decreasing on either side of a central crest. The proof introduces a new function which is related to the so-called flow force and has several surprising properties. In addition to solitary waves, our nonexistence result applies to "half-solitary" waves (e.g. bores) which decay in only one direction.


title: Active Brownian Heat Engines

abstract: When do non-equilibrium forms of disordered energy qualify as heat? \textcolor{blue}{We address this question in the context of cyclically operating heat engines in contact with a non-equilibrium energy reservoir that defies the zeroth law of thermodynamics. To consistently address the latter as a heat bath requires the existence of a precise mapping to an equivalent cycle with an equilibrium bath at a time-dependent effective temperature. We identify the most general setup for which this can generically be ascertained and thoroughly discuss an analytically tractable, experimentally relevant scenario}: a Brownian particle confined in a \textcolor{blue}{periodically} modulated harmonic potential and coupled to some non-equilibrium bath of variable activity. We deduce formal limitations for its thermodynamic performance, including maximum efficiency, efficiency at maximum power, and maximum efficiency at fixed power. They can guide the design of new micro-machines and clarify how much these can outperform passive-bath designs, which has been a debated issue for recent experimental realizations. To illustrate the general principles for practical quasi-static and finite-rate protocols, we further analyze a specific realization of such an active heat engine based on the paradigmatic Active Brownian Particle (ABP) model. This reveals some non-intuitive features of the explicitly computed dynamical effective temperature, illustrates various conceptual and practical limitations of the effective-equilibrium mapping, and clarifies the operational relevance of various coarse-grained measures of dissipation.


title: Range-Max Enhanced Ultra-Wideband Micro-Doppler Signatures of Behind Wall Indoor Human Activities

abstract: Penetrating detection and recognition of behind wall indoor human activities has drawn great attentions from social security and emergency service department in recent years since intelligent surveillance aforehand could avail the proper decision making before operations being carried out. However, due to the influence of the wall effects, the obtained micro-Doppler signatures would be severely degenerated by strong near zero-frequency DC components, which would inevitably smear the detailed characteristic features of different behind wall motions in time-frequency (TF) map and further hinder the motion recognition and classification. In this paper, an ultra-wideband (UWB) radar system is first employed to probe through the opaque wall to detect the behind wall motions, which often span a certain number of range bin cells. By employing such a system, a high resolution range map can be obtained, in which the embedded rich range information is expected to be fully exploited to improve the subsequent recognition and classification performance. Secondly, a high-pass filter is applied to remove the effect of the wall in the raw range map. Then, with the aim of enhancing the characteristic features of different behind wall motions in TF maps, a novel range-max enhancement strategy is proposed to extract the most significant micro-Doppler feature of each TF cell along all range bins for a specific motion. Lastly, the effectiveness of the proposed micro-Doppler signature enhancement strategy is investigated by means of onsite experiments and comparative classification. Both the feature enhanced TF maps and classification results show that the proposed approach outperforms other state-of-art Short-Time Fourier Transform (STFT) based TF feature extraction methods.


title: Observation of a smooth polaron-molecule transition in a degenerate Fermi gas

abstract: Understanding the behavior of an impurity strongly interacting with a Fermi sea is a long-standing challenge in many-body physics. When the interactions are short-ranged, two vastly different ground states exist: a polaron quasiparticle and a molecule dressed by the majority atoms. In the single-impurity limit, it is predicted that at a critical interaction strength, a first-order transition occurs between these two states. Experiments, however, are always conducted in the finite temperature and impurity density regime. The fate of the polaron-to-molecule transition under these conditions, where the statistics of quantum impurities and thermal effects become relevant, is still unknown. Here, we address this question experimentally and theoretically. Our experiments are performed with a spin-imbalanced ultracold Fermi gas with tunable interactions. Utilizing a novel Raman spectroscopy combined with a high-sensitivity fluorescence detection technique, we isolate the quasiparticle contribution and extract the polaron energy, spectral weight, and the contact parameter. As the interaction strength is increased, we observe a continuous variation of all observables, in particular a smooth reduction of the quasiparticle weight as it goes to zero beyond the transition point. Our observation is in good agreement with a theoretical model where polaron and molecule quasiparticle states are thermally occupied according to their quantum statistics. At the experimental conditions, polaron states are hence populated even at interactions where the molecule is the ground state and vice versa. The emerging physical picture is thus that of a smooth transition between polarons and molecules and a coexistence of both in the region around the expected transition.


title: PF: A C++ Library for Fast Particle Filtering

abstract: Particle filters are a class of algorithms that are used for "tracking" or "filtering" in real-time for a wide array of time series models. Despite their comprehensive applicability, particle filters are not always the tool of choice for many practitioners, due to how difficult they are to implement. This short article presents PF, a C++ header-only template library that provides fast implementations of many different particle filters. A tutorial along with an extensive fully-worked example is provided.


title: Radiation detection and energy conversion in nuclear reactor environments by hybrid photovoltaic perovskites

abstract: Detection and direct power conversion of high energy and high intensity ionizing radiation could be a key element in next generation nuclear reactor safety systems and space-born devices. For example, the Fukushima catastrophe in 2011 could have been largely prevented if 1\% of the reactor's remnant radiation ($\gamma$-rays of the nuclear fission) were directly converted within the reactor to electricity to power the water cooling circuit. It is reported here that the hybrid halide perovskite methylammonium lead triiodide could perfectly play the role of a converter. Single crystals were irradiated by a typical shut-down $\gamma$-spectrum of a nuclear reactor with 7.61E14 Bq activity exhibit a high-efficiency of $\gamma$-ray to free charge carrier conversion with radiation hardening. The power density of 0.3 mW/kg of methylammonium lead triiodide at 50 Sv/h means a four times higher efficiency than that for silicon-based cells. The material was stable to the limits of the experiment without changing its performance up to 100 Sv/h dose rate and 57 Sv H*(10) ambient total $\gamma$-dose. Moreover, the {\gamma}-shielding performance of methylammonium lead triiodide was found to be superior to both ordinary and barite concrete.


title: Limit theorems for a stable sausage

abstract: In this article, we study fluctuations of the volume of a stable sausage defined via a $d$-dimensional rotationally invariant $\alpha$-stable process. As the main results, we establish a functional central limit theorem (in the case when $d/\alpha>3 /2$) with a standard one-dimensional Brownian motion in the limit, and Khintchine's and Chung's laws of the iterated logarithm (in the case when $d/\alpha>9 /5$).


title: Stimulated-emission based model of fast radio bursts

abstract: Fast radio bursts (FRBs) are bright, short-duration radio transients with very high brightness temperatures implying highly coherent emission. We suggest that the FRBs are caused by the self-focusing of an electron beam interacting with an ambient plasma right beyond the light cylinder radius of a neutron star. The magnetic field at the light cylinder radius is relatively high which can accommodate both young Crab-like systems and old millisecond pulsars addressing the diverse environments of FRBs. At the first stage, the intense pulsed-beam passing through the background plasma causes instabilities such that the trapped particles in local Buneman-type cavitons saturate the local field. The beam is then radially self-focused due to the circular electric field developed by the two-stream instability which leads to Weibel instability in the transverse direction. Finally, the non-linear saturation of the Weibel instability results in the self-modulational formation of solitons due to plasmoid instability. The resonant solitary waves are the breather-type solitons hosting relativistic particles with self-excited oscillations. The analytical solutions obtained for non-linear dispersion and solitons suggest that, near the current sheets, the relativistic bunches are accelerated/amplified by klystron-like structures due to self-excited oscillations by the induced local electric field. Boosted coherent radio emission propagates through a narrow cone with strong focusing due to radial electric field and magnetic pinching. The non-linear evolution of solitons and the stimulated emission are associated with the Buneman instability and the possibility of the presence of nanosecond shots in FRBs are investigated.


title: Robust edge modes in dislocated systems of subwavelength resonators

abstract: Robustly manipulating waves on subwavelength scales can be achieved by, firstly, designing a structure with a subwavelength band gap and, secondly, introducing a defect so that eigenfrequencies fall within the band gap. Such frequencies are well known to correspond to localized modes. We study a one-dimensional array of subwavelength resonators, proving that there is a subwavelength band gap, and showing that by introducing a dislocation we can place localized modes at any point within the band gap. We complement this analysis by studying the stability properties of the corresponding finite array of resonators, demonstrating the value of being able to customize the position of eigenvalues within the band gap.


title: A Cancellation Nanoflare Model for Solar Chromospheric and Coronal Heating III. 3D Simulations and Atmospheric Response

abstract: Inspired by recent observations suggesting that photospheric magnetic flux cancellation occurs much more frequently than previously thought, we analytically estimated the energy released from reconnection driven by photospheric flux cancellation, and proposed that it can act as a mechanism for chromospheric and coronal heating (Priest et al., 2018). Using two-dimensional simulations we validated the analytical estimates and studied the resulting atmospheric response (Syntelis et al. 2019). In the present work, we set up three-dimensional resistive MHD simulations of two cancelling polarities in a stratified atmosphere with a horizontal external field to further validate and improve upon the analytical estimates. The computational evaluation of the parameters associated with the energy release are in good qualitative agreement with the analytical estimates. The computational Poynting energy flux into the current sheet is in good qualitative agreement with the analytical estimates, after correcting the analytical expression to better account for the horizontal extent of the current sheet. The atmospheric response to the cancellation is the formation of hot ejections, cool ejections, or a combination of both hot and cool ejections, which can appear with a time difference and/or be spatially offset, depending on the properties of the cancelling region and the resulting height of the reconnection. Therefore, during the cancellation, a wide spectrum of ejections can be formed, which can account for the variety of multi-thermal ejections associated with Ellerman bombs, UV bursts and IRIS bombs, and also other ejections associated with small-scale cancelling regions and spicules.