title: mmWave Channel Estimation via Approximate Message Passing with Side Information

abstract: This work considers millimeter-wave channel estimation in a setting where parameters of the underlying mmWave channels are varying dynamically over time and there is a single drifting path. In this setting, channel estimates at time block $k$ can be used as side information (SI) when estimating the channel at block $k+1$. To estimate channel parameters, we employ an SI-aided (complex) approximate message passing algorithm and compare its performance to a benchmark based on orthogonal matching pursuit.


title: Overcoming the inconsistences of the variance inflation factor: a redefined VIF and a test to detect statistical troubling multicollinearity

abstract: Multicollinearity is relevant to many different fields where linear regression models are applied, and its existence may affect the analysis of ordinary least squares (OLS) estimators from both the numerical and statistical points of views. Thus, multicollinearity can lead to incoherence in the statistical significance of the independent variables and the global significance of the model. The variance inflation factor (VIF) is traditionally applied to diagnose the possible existence of multicollinearity, but it is not always the case that detection by VIF of a troubling degree of multicollinearity corresponds to negative effects on the statistical analysis. The reason for the lack of specificity of VIF is that there are other factors, such as the size of the sample and the variance of the random disturbance, that can lead to high values of the VIF but not to problematic variance in the OLS estimators (see O'Brien 2007). This paper presents a new variance inflation factor (TVIF) that consider all these additional factors. Thresholds for this new measure and from the index provided by Stewart (1987) are also provided. These thresholds are reinterpreted and presented as a new statistical test to diagnose the existence of statistical troubling multicollinearity. The contributions of this paper are illustrated with two real data examples previously applied in the scientific literature.


title: Hodge symmetry for rigid varieties via log hard Lefschetz

abstract: Motivated by a question of Hansen and Li, we show that a smooth and proper rigid analytic space $X$ with projective reduction satisfies Hodge symmetry in the following situations: (1) the base non-archimedean field $K$ is of residue characteristic zero, (2) $K$ is $p$-adic and $X$ has good ordinary reduction, (3) $K$ is $p$-adic and $X$ has "combinatorial reduction."' We also reprove a version of their result, Hodge symmetry for $H^1$, without the use of moduli spaces of semistable sheaves. All of this relies on cases of Kato's log hard Lefschetz conjecture, which we prove for $H^1$ and for log schemes of "combinatorial type."


title: A Linear Algebra Approach to Linear Metatheory

abstract: Linear typed $\lambda$-calculi are more delicate than their simply typed siblings when it comes to metatheoretic results like preservation of typing under renaming and substitution. Tracking the usage of variables in contexts places more constraints on how variables may be renamed or substituted. We present a methodology based on linear algebra over semirings, extending McBride's kits and traversals approach for the metatheory of syntax with binding to linear usage-annotated terms. Our approach is readily formalisable, and we have done so in Agda.


title: Error bounds for dynamical spectral estimation

abstract: Dynamical spectral estimation is a well-established numerical approach for estimating eigenvalues and eigenfunctions of the Markov transition operator from trajectory data. Although the approach has been widely applied in biomolecular simulations, its error properties remain poorly understood. Here we analyze the error of a dynamical spectral estimation method called "the variational approach to conformational dynamics" (VAC). We bound the approximation error and estimation error for VAC estimates. Our analysis establishes VAC's convergence properties and suggests new strategies for tuning VAC to improve accuracy.


title: A robust algorithm for explaining unreliable machine learning survival models using the Kolmogorov-Smirnov bounds

abstract: A new robust algorithm based of the explanation method SurvLIME called SurvLIME-KS is proposed for explaining machine learning survival models. The algorithm is developed to ensure robustness to cases of a small amount of training data or outliers of survival data. The first idea behind SurvLIME-KS is to apply the Cox proportional hazards model to approximate the black-box survival model at the local area around a test example due to the linear relationship of covariates in the model. The second idea is to incorporate the well-known Kolmogorov-Smirnov bounds for constructing sets of predicted cumulative hazard functions. As a result, the robust maximin strategy is used, which aims to minimize the average distance between cumulative hazard functions of the explained black-box model and of the approximating Cox model, and to maximize the distance over all cumulative hazard functions in the interval produced by the Kolmogorov-Smirnov bounds. The maximin optimization problem is reduced to the quadratic program. Various numerical experiments with synthetic and real datasets demonstrate the SurvLIME-KS efficiency.


title: Homogeneous sets, clique-separators, critical graphs, and optimal $\chi$-binding functions

abstract: Given a set $\mathcal{H}$ of graphs, let $f_\mathcal{H}^\star\colon \mathbb{N}_{>0}\to \mathbb{N}_{>0}$ be the optimal $\chi$-binding function of the class of $\mathcal{H}$-free graphs, that is, $$f_\mathcal{H}^\star(\omega)=\max\{\chi(G): G\text{ is } \mathcal{H}\text{-free, } \omega(G)=\omega\}.$$ In this paper, we combine the two decomposition methods by homogeneous sets and clique-separators in order to determine optimal $\chi$-binding functions for subclasses of $P_5$-free graphs and of $(C_5,C_7,\ldots)$-free graphs. In particular, we prove the following for each $\omega\geq 1$: (i) $\ f_{\{P_5,banner\}}^\star(\omega)=f_{3K_1}^\star(\omega)\in \Theta(\omega^2/\log(\omega)),$ (ii) $\ f_{\{P_5,co-banner\}}^\star(\omega)=f^\star_{\{2K_2\}}(\omega)\in\mathcal{O}(\omega^2),$ (iii) $\ f_{\{C_5,C_7,\ldots,banner\}}^\star(\omega)=f^\star_{\{C_5,3K_1\}}(\omega)\notin \mathcal{O}(\omega),$ and (iv) $\ f_{\{P_5,C_4\}}^\star(\omega)=\lceil(5\omega-1)/4\rceil.$ We also characterise, for each of our considered graph classes, all graphs $G$ with $\chi(G)>\chi(G-u)$ for each $u\in V(G)$. From these structural results, we can prove Reed's conjecture -- relating chromatic number, clique number, and maximum degree of a graph -- for $(P_5,banner)$-free graphs.


title: Raccoons vs Demons: multiclass labeled P300 dataset

abstract: We publish dataset of visual P300 BCI performed in Virtual Reality (VR) game Raccoons versus Demons (RvD). Data contains reach labels incorporating information about stimulus chosen enabling us to estimate model's confidence at each stimulus prediction stage. Data and experiments code are available at https://gitlab.com/impulse-neiry_public/raccoons-vs-demons


title: The Role of the Effective Range in Resonantly Interacting Fermi Gases: How Breaking Scale Symmetry Affects the Bulk Viscosity

abstract: We investigate the role of the effective range on the bulk viscosity of $s$- and $p$-wave Fermi gases. At resonance, the presence of the effective range breaks the scale invariance of the system, and hence results in a non-zero bulk viscosity. However, we show that the effective range plays a very different role in the two cases. In the $s$-wave case, the role of the effective range is perturbative, and its contribution to the bulk viscosity vanishes in the limit of zero effective range. On the other hand, the effective range in $p$-wave Fermi gases leads to a non-zero bulk viscosity, even in the zero-range limit. We setup the general diagrammatic approach to compute the bulk viscosity spectral function that includes the effects of the effective range. We then compute the analytic expressions for the spectral function in the high temperature limit, at low- and high-frequencies. We also derive the sum rules for the bulk viscosity spectral function for both $s$- and $p$-wave gases.


title: Quasi-linear irreversible thermodynamics of a low-temperature-differential kinematic Stirling heat engine

abstract: Low-temperature-differential (LTD) Stirling heat engines are able to operate with a small temperature difference between low-temperature heat reservoirs that exist in our daily lives, and thus they are considered to be an important sustainable energy technology. The author recently proposed a nonlinear dynamics model of an LTD kinematic Stirling heat engine to study the rotational mechanism of the engine [Y. Izumida, EPL \textbf{121}, 50004 (2018)]. This paper presents our study of the nonequilibrium thermodynamics analysis of this engine model, where a load torque against which the engine does work is introduced. We demonstrate that the engine's rotational state is in a quasi-linear response regime where the thermodynamic fluxes show a linear dependency on the thermodynamic forces. Significantly, it is found that the response coefficients of the quasi-linear relations are symmetric, which is similar to Onsager symmetry in linear irreversible thermodynamics. Based on these relations, we formulate the maximum efficiency of the engine. We also elucidate that the symmetry of the quasi-linear response coefficients emerges by reflecting the (anti-)reciprocity of the Onsager kinetic coefficients identified for the relaxation dynamics of the engine in the vicinity of an equilibrium state. We expect that the present study paves the way for developing nonequilibrium thermodynamics of autonomous heat engines described as a nonlinear dynamical system.


title: Fluctuations of extreme eigenvalues of sparse Erd\H{o}s-R\'enyi graphs

abstract: We consider a class of sparse random matrices which includes the adjacency matrix of the Erd\H{o}s-R\'enyi graph $\mathcal{G}(N,p)$. We show that if $N^{\varepsilon} \leq Np \leq N^{1/3-\varepsilon}$ then all nontrivial eigenvalues away from 0 have asymptotically Gaussian fluctuations. These fluctuations are governed by a single random variable, which has the interpretation of the total degree of the graph. This extends the result [19] on the fluctuations of the extreme eigenvalues from $Np \geq N^{2/9 + \varepsilon}$ down to the optimal scale $Np \geq N^{\varepsilon}$. The main technical achievement of our proof is a rigidity bound of accuracy $N^{-1/2-\varepsilon} \, (Np)^{-1/2}$ for the extreme eigenvalues, which avoids the $(Np)^{-1}$-expansions from [9,19,24]. Our result is the last missing piece, added to [8, 12, 19, 24], of a complete description of the eigenvalue fluctuations of sparse random matrices for $Np \geq N^{\varepsilon}$.


title: Kernels of Perturbed Toeplitz Operators in vector-valued Hardy spaces

abstract: Recently, Liang and Partington \cite{YP} show that kernels of finite-rank perturbations of Toeplitz operators are nearly invariant with finite defect under the backward shift operator acting on the scalar-valued Hardy space. In this article we provide a vectorial generalization of a result of Liang and Partington. As an immediate application we identify the kernel of perturbed Toeplitz operator in terms of backward shift-invariant subspaces in various important cases by applying the recent theorem (\cite{CDP, OR}) in connection with nearly invariant subspaces of finite defect for the backward shift operator acting on the vector-valued Hardy space.


title: Regional Boundary Gradient Strategic Sensors Characterizations

abstract: The purpose of this paper is to characterize regional boundary gradient strategic sensors notions for different cases of regional boundary gradient observation to be achieved. Then, the characterizations based on how to a cross from internal gradient region to the boundary gradient region. Thus, the obtained results are applied in two dimensional linear infinite distributed systems in Hilbert space where the dynamics are governed by strongly continuous semi-group. Moreover, we give the relation between the regional gradient observability on a subregion small Gama of the spatial domain Omega with the regional boundary gradient observability on a subregion Gama of the boundary Omega of. Finally, sufficient conditions of regional boundary gradient strategic sensors notions are explored, analyzed and discussed in connection with the regional boundary gradient of exact (weak) observability, positive definite observability operator and rank conditions.


title: Flexoelectric radial polarization of single-walled nanotubes from first-principles

abstract: Flexoelectricity is the polar response of an insulator to strain gradients such as bending. While the size dependence of this effect makes it weak in bulk systems in comparison to piezoelectricity, it suggests that it plays a bigger role in nanoscale systems such as thin films and nanotubes (NTs). In this paper we demonstrate using first-principles calculations that the walls of carbon nanotubes (CNTs) and transition metal dichalcogenide nanotubes (TMD NTs) are polarized in the radial direction, the strength of the polarization increasing as the size of the NT decreases. For CNTs and TMD NTs with chiral indices $(n,m)$, the radial polarization of the walls $P_R$ starts to diverge below $ C(n,m)/a = \sqrt{n^2 +nm + m^2} \sim 10$, where $C(n,m)$ is the circumference of the NT and $a$ is the lattice constant of the 2D monolayer. For CNTs, $P_R$ drops to zero above this value but for TMD NTs there is a non-zero polarization above this value, which is an ionic rather than electronic effect. The size dependence of $P_R$ in the TMD NTs is interesting: it increases gradually and reaches a maximum of $P_R \sim 100 \ \text{C/cm}^2$ at $C(n,m)/a \sim 15$, then decreases until $C(n,m)/a \sim 10$ where it starts to diverge. Measurements of the radial strain on the chalcogen atoms with respect to the 2D monolayers shows that this polarization is the result of a significantly larger strain on the outer bonds than the inner bonds, but did not offer an explanation for the peculiar size dependence. These results suggest that while the walls of smaller CNTs and TMD NTs are polarized, the walls of larger TMD NTs are also polarized due to a difference in strain on the inner and outer atoms in the walls. This result may prove useful for the application of NTs for screening in liquid or biological systems.


title: AGE Challenge: Angle Closure Glaucoma Evaluation in Anterior Segment Optical Coherence Tomography

abstract: Angle closure glaucoma (ACG) is a more aggressive disease than open-angle glaucoma, where the abnormal anatomical structures of the anterior chamber angle (ACA) may cause an elevated intraocular pressure and gradually lead to glaucomatous optic neuropathy and eventually to visual impairment and blindness. Anterior Segment Optical Coherence Tomography (AS-OCT) imaging provides a fast and contactless way to discriminate angle closure from open angle. Although many medical image analysis algorithms have been developed for glaucoma diagnosis, only a few studies have focused on AS-OCT imaging. In particular, there is no public AS-OCT dataset available for evaluating the existing methods in a uniform way, which limits progress in the development of automated techniques for angle closure detection and assessment. To address this, we organized the Angle closure Glaucoma Evaluation challenge (AGE), held in conjunction with MICCAI 2019. The AGE challenge consisted of two tasks: scleral spur localization and angle closure classification. For this challenge, we released a large dataset of 4800 annotated AS-OCT images from 199 patients, and also proposed an evaluation framework to benchmark and compare different models. During the AGE challenge, over 200 teams registered online, and more than 1100 results were submitted for online evaluation. Finally, eight teams participated in the onsite challenge. In this paper, we summarize these eight onsite challenge methods and analyze their corresponding results for the two tasks. We further discuss limitations and future directions. In the AGE challenge, the top-performing approach had an average Euclidean Distance of 10 pixels (10um) in scleral spur localization, while in the task of angle closure classification, all the algorithms achieved satisfactory performances, with two best obtaining an accuracy rate of 100%.


title: Learning programs by learning from failures

abstract: We describe an inductive logic programming (ILP) approach called learning from failures. In this approach, an ILP system (the learner) decomposes the learning problem into three separate stages: generate, test, and constrain. In the generate stage, the learner generates a hypothesis (a logic program) that satisfies a set of hypothesis constraints (constraints on the syntactic form of hypotheses). In the test stage, the learner tests the hypothesis against training examples. A hypothesis fails when it does not entail all the positive examples or entails a negative example. If a hypothesis fails, then, in the constrain stage, the learner learns constraints from the failed hypothesis to prune the hypothesis space, i.e. to constrain subsequent hypothesis generation. For instance, if a hypothesis is too general (entails a negative example), the constraints prune generalisations of the hypothesis. If a hypothesis is too specific (does not entail all the positive examples), the constraints prune specialisations of the hypothesis. This loop repeats until either (i) the learner finds a hypothesis that entails all the positive and none of the negative examples, or (ii) there are no more hypotheses to test. We introduce Popper, an ILP system that implements this approach by combining answer set programming and Prolog. Popper supports infinite problem domains, reasoning about lists and numbers, learning textually minimal programs, and learning recursive programs. Our experimental results on three domains (toy game problems, robot strategies, and list transformations) show that (i) constraints drastically improve learning performance, and (ii) Popper can outperform existing ILP systems, both in terms of predictive accuracies and learning times.


title: Some observations on the properness of identity plus linear powers: part 2

abstract: This paper develops our previous work on properness of a class of maps related to the Jacobian conjecture. The paper has two main parts: - In part 1, we explore properties of the set of non-proper values $S_f$ (as introduced by Z. Jelonek) of these maps. In particular, using a general criterion for non-properness of these maps, we show that under a "generic condition" (to be precise later) $S_f$ contains $0$ if it is non-empty. This result is related to a conjecture in our previous paper. We obtain this by use of a "dual set" to $S_f$, particularly designed for the special class of maps. - In part 2, we use the non-properness criteria obtained in our work to construct a counter-example to the proposed proof in arXiv:2002.10249 of the Jacobian conjecture. In the conclusion, we present some comments pertaining the Jacobian conjecture and properness of polynomial maps in general.


title: Crucial Inflammatory Mediators and Efficacy of Drug Interventions in Pneumonia Inflated COVID-19: An Invivo Mathematical Modelling Study

abstract: The virus SARS-COV-2 caused disease COVID-19 has been declared a pandemic by WHO. Currently, over 210 countries and territories have been affected. Careful, well-designed drugs and vaccine for the total elimination of this virus seem to be the need of the hour. In this context, the invivo mathematical modelling studies can be extremely helpful in understanding the efficacy of the drug interventions. These studies can also help understand the role of the crucial inflammatory mediators and the behaviour of immune response towards this novel coronavirus. Motivated by these facts, in this paper, we study the invivo dynamics of Covid-19. The results obtained here are inline with some of the clinical findings for Covid-19. This invivo modelling study involving the crucial biomarkers of Covid-19 is the first of its kind and the results obtained from this can be helpful to researchers, epidemiologists, clinicians and doctors who are working in this field.


title: PolymoRF: Polymorphic Wireless Receivers Through Physical-Layer Deep Learning

abstract: Today's wireless technologies are largely based on inflexible designs, which makes them inefficient and prone to a variety of wireless attacks. To address this key issue, wireless receivers will need to (i) infer on-the-fly the physical-layer parameters currently used by transmitters; and if needed, (ii) change their hardware and software structures to demodulate the incoming waveform. In this paper, we introduce PolymoRF, a deep learning-based polymorphic receiver able to reconfigure itself in real time based on the inferred waveform parameters. Our key technical innovations are (i) a novel embedded deep learning architecture, called RFNet, which enables the solution of key waveform inference problems; (ii) a generalized hardware/software architecture that integrates RFNet with radio components and signal processing. We prototype PolymoRF on a custom software-defined radio platform, and show through extensive over-the-air experiments that (i) RFNet achieves similar accuracy to that of state-of-the-art yet with 52x and 8x latency and hardware reduction; (ii) PolymoRF achieves throughput within 87% of a perfect-knowledge Oracle system, thus demonstrating for the first time that polymorphic receivers are feasible and effective.


title: Trace of canonical modules, annihilator of Ext, and classes of rings close to being Gorenstein

abstract: In this note we study trace ideals of canonical modules. Characterizations of the trace ideals in terms of annihilators of certain Ext modules are given. We apply our results to study many classes of rings close to being Gorenstein that appear in recent literature. We discuss the behavior of nearly Gorensteinness, introduced by Herzog-Hibi-Stamate, under reductions by regular sequences. A nearly Gorenstein version of the Tachikawa conjecture is posed, and an affirmative answer given for type 2 rings. We also introduce the class of weakly almost Gorenstein rings using the canonical module, and show that they are almost Gorenstein rings in dimension one, and are closely related to Teter rings in the artinian case. We explain their basic properties and give some examples.


title: LandCover.ai: Dataset for Automatic Mapping of Buildings, Woodlands, Water and Roads from Aerial Imagery

abstract: Monitoring of land cover and land use is crucial in natural resources management. Automatic visual mapping can carry enormous economic value for agriculture, forestry, or public administration. Satellite or aerial images combined with computer vision and deep learning enable precise assessment and can significantly speed up change detection. Aerial imagery usually provides images with much higher pixel resolution than satellite data allowing more detailed mapping. However, there is still a lack of aerial datasets made for the segmentation, covering rural areas with a resolution of tens centimeters per pixel, manual fine labels, and highly publicly important environmental instances like buildings, woods, water, or roads. Here we introduce LandCover.ai (Land Cover from Aerial Imagery) dataset for semantic segmentation. We collected images of 216.27 sq. km rural areas across Poland, a country in Central Europe, 39.51 sq. km with resolution 50 cm per pixel and 176.76 sq. km with resolution 25 cm per pixel and manually fine annotated four following classes of objects: buildings, woodlands, water, and roads. Additionally, we report simple benchmark results, achieving 85.56% of mean intersection over union on the test set. It proves that the automatic mapping of land cover is possible with a relatively small, cost-efficient, RGB-only dataset. The dataset is publicly available at https://landcover.ai.linuxpolska.com/


title: Predicting atmospheric optical properties for radiative transfer computations using neural networks

abstract: The radiative transfer equations are well-known, but radiation parametrizations in atmospheric models are computationally expensive. A promising tool for accelerating parametrizations is the use of machine learning techniques. In this study, we develop a machine learning-based parametrization for the gaseous optical properties by training neural networks to emulate a modern radiation parameterization (RRTMGP). To minimize computational costs, we reduce the range of atmospheric conditions for which the neural networks are applicable and use machine-specific optimised BLAS functions to accelerate matrix computations. To generate training data, we use a set of randomly perturbed atmospheric profiles and calculate optical properties using RRTMGP. Predicted optical properties are highly accurate and the resulting radiative fluxes have average errors within \SI{0.5}{\flux} compared to RRTMGP. Our neural network-based gas optics parametrization is up to 4 times faster than RRTMGP, depending on the size of the neural networks. We further test the trade-off between speed and accuracy by training neural networks for the narrow range of atmospheric conditions of a single large-eddy simulation, so smaller and therefore faster networks can achieve a desired accuracy. We conclude that our machine learning-based parametrization can speed-up radiative transfer computations whilst retaining high accuracy.


title: Optimized FTO seeding enables the growth of high efficient Ta-doped TiO$_2$ nanorod photoanodes

abstract: Tantalum doped rutile nanorods were hydrothermally grown on FTO substrates using a new seeding approach. This approach allows the incorporation of high concentrations of up to 4.8 at% tantalum as active doping and results in a significant enhancement of photoelectrochemical water splitting rate (1.8 mA/cm2 at a potential of +1.5 V vs RHE) which corresponds to ca. 1% photocurrent conversion efficiency under AM 1.5, 100 mW/cm2 simulated sunlight irradiation.


title: Testing triplet fermions at the electron-positron and electron-proton colliders using fat jet signatures

abstract: The addition of $SU(2)_L$ triplet fermions of zero hypercharge with the Standard Model (SM) helps to explain the origin of the neutrino mass by the so-called seesaw mechanism. Such a scenario is commonly know as the type-III seesaw model. After the electroweak symmetry breaking the mixings between the light and heavy mass eigenstates of the neutral leptons are developed which play important roles in the study of the charged and neutral multiplets of the triplet fermions at the colliders. In this article we study such interactions to produce these multiplets of the triplet fermion at the electron-positron and electron-proton colliders at different center of mass energies. We focus on the heavy triplets, for example, having mass in the TeV scale so that their decay products including the SM the gauge bosons or Higgs boson can be sufficiently boosted, leading to a fat jet. Hence we probe the mixing between light-heavy mass eigenstates of the neutrinos and compare the results with the bounds obtained by the electroweak precision study.


title: Breaking RSA Security With A Low Noise D-Wave 2000Q Quantum Annealer: Computational Times, Limitations And Prospects

abstract: The RSA cryptosystem could be easily broken with large scale general purpose quantum computers running Shor's factorization algorithm. Being such devices still in their infancy, a quantum annealing approach to integer factorization has recently gained attention. In this work, we analyzed the most promising strategies for RSA hacking via quantum annealing with an extensive study of the low noise D-Wave 2000Q computational times, current hardware limitations and challenges for future developments.


title: Towards explainable classifiers using the counterfactual approach -- global explanations for discovering bias in data

abstract: The paper proposes summarized attribution-based post-hoc explanations for the detection and identification of bias in data. A global explanation is proposed, and a step-by-step framework on how to detect and test bias is introduced. Since removing unwanted bias is often a complicated and tremendous task, it is automatically inserted, instead. Then, the bias is evaluated with the proposed counterfactual approach. The obtained results are validated on a sample skin lesion dataset. Using the proposed method, a number of possible bias causing artifacts are successfully identified and confirmed in dermoscopy images. In particular, it is confirmed that black frames have a strong influence on Convolutional Neural Network's prediction: 22% of them changed the prediction from benign to malignant.


title: Hacking the Waveform: Generalized Wireless Adversarial Deep Learning

abstract: This paper advances the state of the art by proposing the first comprehensive analysis and experimental evaluation of adversarial learning attacks to wireless deep learning systems. We postulate a series of adversarial attacks, and formulate a Generalized Wireless Adversarial Machine Learning Problem (GWAP) where we analyze the combined effect of the wireless channel and the adversarial waveform on the efficacy of the attacks. We propose a new neural network architecture called FIRNet, which can be trained to "hack" a classifier based only on its output. We extensively evaluate the performance on (i) a 1,000-device radio fingerprinting dataset, and (ii) a 24-class modulation dataset. Results obtained with several channel conditions show that our algorithms can decrease the classifier accuracy up to 3x. We also experimentally evaluate FIRNet on a radio testbed, and show that our data-driven blackbox approach can confuse the classifier up to 97% while keeping the waveform distortion to a minimum.


title: Estimating the number of available states for normal and tumor tissues in gene expression space

abstract: The topology of gene expression space for a set of 12 cancer types is studied by means of an entropy-like magnitude, which allows the characterization of the regions occupied by tumor and normal samples. The comparison indicates that the number of available states in gene expression space is much greater for tumors than for normal tissues, suggesting the irreversibility of the progression to the tumor phase. The entropy is nearly constant for tumors, whereas exhibits a higher variability in normal tissues, probably due to tissue differentiation. In addition, we show an interesting correlation between the fraction of available states and the overlapping between the tumor and normal sample clouds, interpreted as a way of reducing the decay rate to the tumor phase in more ordered or structured tissues.


title: Modeling Urban Growth and Form with Spatial Entropy

abstract: Entropy is one of physical bases for fractal dimension definition, and the generalized fractal dimension was defined by Renyi entropy. Using fractal dimension, we can describe urban growth and form and characterize spatial complexity. A number of fractal models and measurements have been proposed for urban studies. However, the precondition for fractal dimension application is to find scaling relations in cities. In absence of scaling property, we can make use of entropy function and measurements. This paper is devoted to researching how to describe urban growth by using spatial entropy. By analogy with fractal dimension growth models of cities, a pair of entropy increase models can be derived and a set of entropy-based measurements can be constructed to describe urban growing process and patterns. First, logistic function and Boltzmann equation are utilized to model the entropy increase curves of urban growth. Second, a series of indexes based on spatial entropy are used to characterize urban form. Further, multifractal dimension spectrums are generalized to spatial entropy spectrums. Conclusions are drawn as follows. Entropy and fractal dimension have both intersection and different spheres of application to urban research. Thus, for a given spatial measurement scale, fractal dimension can often be replaced by spatial entropy for simplicity. The models and measurements presented in this work are significant for integrating entropy and fractal dimension into the same framework of urban spatial analysis and understanding spatial complexity of cities.


title: Systematic study on the quark-hadron mixed phase in compact stars

abstract: We investigate systematically the quark-hadron mixed phase in dense stellar matter, and its influence on compact star structures. The properties of quark matter and hadronic matter are fixed based on various model predictions. Beside adopting constant values, the surface tension $\Sigma$ for the quark-hadron interface is estimated with the multiple reflection expansion method and equivparticle model. To fix the structures of quark-hadron pasta phases, a continuous dimensionality of the structure is adopted as proposed by Ravenhall, Pethick, and Wilson. The corresponding properties of hybrid stars are then obtained and confronted with pulsar observations. It is found that the correlation between radius and tidal deformability in traditional neutron stars preserves in hybrid stars. For those permitted by pulsar observations, in almost all cases the quark phase persists inside the most massive compact stars. The quark-hadron interface plays an important role on hybrid star structures once quark matter emerges. The surface tension $\Sigma$ estimated with various methods increases with density, which predicts stiffer EOSs for the quark-hadron mixed phase and increases the maximum mass of hybrid stars. The EOSs of hybrid star matter are well constrained at densities $n\lesssim 0.8$ fm${}^{-3}$, while larger uncertainty is expected at higher densities.


title: Online Convex Optimization with Binary Constraints

abstract: We consider online optimization with binary decision variables and convex loss functions. We design a new algorithm, binary online gradient descent (bOGD) and bound its expected dynamic regret. We provide a regret bound that holds for any time horizon and a specialized bound for finite time horizons. First, we present the regret as the sum of the relaxed, continuous round optimum tracking error and the rounding error of our update in which the former asymptomatically decreases with time under certain conditions. Then, we derive a finite-time bound that is sublinear in time and linear in the cumulative variation of the relaxed, continuous round optima. We apply bOGD to demand response with thermostatically controlled loads, in which binary constraints model discrete on/off settings. We also model uncertainty and varying load availability, which depend on temperature deadbands, lockout of cooling units and manual overrides. We test the performance of bOGD in several simulations based on demand response. The simulations corroborate that the use of randomization in bOGD does not significantly degrade performance while making the problem more tractable.


title: Masur--Veech volumes of quadratic differentials and their asymptotics

abstract: Based on the Chen--M\"oller--Sauvaget formula, we apply the theory of integrable systems to derive three equations for the generating series of the Masur--Veech volumes ${\rm Vol} \, \mathcal{Q}_{g,n}$ associated with the principal strata of the moduli spaces of quadratic differentials, and propose refinements of the conjectural formulas given in [12,4] for the large genus asymptotics of ${\rm Vol} \, \mathcal{Q}_{g,n}$ and of the associated area Siegel--Veech constants.


title: On Absolute Continuity and Singularity of Multidimensional Diffusions

abstract: Consider two laws \(P\) and \(Q\) of multidimensional possibly explosive diffusions with common diffusion coefficient \(\mathfrak{a}\) and drift coefficients \(\mathfrak{b}\) and \(\mathfrak{b} + \mathfrak{a} \mathfrak{c}\), respectively, and the law \(P^\circ\) of an auxiliary diffusion with diffusion coefficient \(\langle \mathfrak{c},\mathfrak{a}\mathfrak{c}\rangle^{-1}\mathfrak{a}\) and drift coefficient \(\langle \mathfrak{c}, \mathfrak{a}\mathfrak{c}\rangle^{-1}\mathfrak{b}\). We show that \(P \ll Q\) if and only if the auxiliary diffusion \(P^\circ\) explodes almost surely and that \(P\perp Q\) if and only if the auxiliary diffusion \(P^\circ\) almost surely does not explode. As applications we derive a Khasminskii-type integral test for absolute continuity and singularity, an integral test for explosion of time-changed Brownian motion, and we discuss applications to mathematical finance.


title: Supercharacter theory for the Borel contraction of the group GL(n,F_q)

abstract: The notion of a supercharacter theory was proposed by P. Diaconis and I.M. Isaacs in 2008. A supercharacter theory for a given finite group is a pair of the system of certain complex characters and the partition of group into classes that have properties similar to the system of irreducible characters and the partition into conjugacy classes. In the present paper, we consider the group obtained by the Borel contraction from the general linear group over a finite field. For this group, we construct the supercharacter theory. In terms of rook placements, we classify supercharacters and superclasses, calculate values of supercharacters on superclasses


title: Measuring the expansion of the universe

abstract: We draw a figure from where it is possible to measure the number of e-folds of expansion of the universe with a ruler. We find model independent bounds for the number of e-folds during inflation, reheating and radiation. We also give a lower bound to the size of the universe at the beginning of observable inflation. Finally, we show that consistency with a relevant diagram requires the existence of a new form of energy to drive the present expansion of the universe.


title: Graph Homomorphisms Based On Particular Total Colorings of Graphs and Graphic Lattices

abstract: Lattice-based cryptography is not only for thwarting future quantum computers, and is also the basis of Fully Homomorphic Encryption. Motivated from the advantage of graph homomorphisms we combine graph homomorphisms with graph total colorings together for designing new types of graph homomorphisms: totally-colored graph homomorphisms, graphic-lattice homomorphisms from sets to sets, every-zero graphic group homomorphisms from sets to sets. Our graph-homomorphism lattices are made up by graph homomorphisms. These new homomorphisms induce some problems of graph theory, for example, Number String Decomposition and Graph Homomorphism Problem.


title: Coronavirus and sports leagues: obtaining a fair ranking when the season cannot resume

abstract: Many sports leagues are played in a tightly scheduled round-robin format, leaving a limited time window to postpone matches. If the season cannot resume due to an external shock such as the coronavirus pandemic in 2020, the ranking of the teams becomes non-trivial: it is necessary to account for schedule imbalances and possibly for the different number of matches played. First in the literature, we identify a set of desired axioms for ranking in these incomplete tournaments. It is verified that the generalised row sum, a parametric family of scoring rules, satisfies all of them. In particular, the well-established least squares method maximizes the influence of the strength of opponents on the ranking. Our approach is applied for six major premier European soccer competitions, where the rankings are found to be robust concerning the weight of adjustment for the strength of the opponents. Some methodologically simpler alternative policies are also discussed.


title: Synchronous oscillations and symmetry breaking in a model of two interacting ultrasound contrast agents

abstract: We study nonlinear dynamics in a system of two coupled oscillators, describing the motion of two interacting microbubble contrast agents. In the case of identical bubbles, the corresponding symmetry of the governing system of equations leads to the possibility of existence of asymptotically stable synchronous oscillations. However, it may be difficult to create absolutely identical bubbles and, moreover, one can observe in experiments regimes that are unstable with respect to perturbations of equilibrium radii of bubbles. Therefore, we investigate the stability of various synchronous and asynchronous dynamical regimes with respect to the breaking of this symmetry. We show that the main factors determining stability or instability of a synchronous attractor are the presence/absence and the type of an asynchronous attractor coexisting with the synchronous attractor. On the other hand, asynchronous hyperchaotic attractors are stable with respect to the symmetry breaking in all the situations we have studied. Therefore, they are likely to be observed in physically realistic scenarios and can be beneficial for suitable applications when chaotic behavior is desirable.


title: Bayesian longitudinal models for exploring European sardine fishing in the Mediterranean Sea

abstract: In the Mediterranean Sea, catches are dominated by small pelagic fish, representing nearly the 49\% of the total harvest. Among them, the European sardine (Sardina pilchardus) is one of the most commercially important species showing high over-exploitation rates in recent last years. In this study we analysed the European sardine landings in the Mediterranean Sea from 1970 to 2014. We made use of Bayesian longitudinal linear mixed models in order to assess differences in the temporal evolution of fishing between and within countries. Furthermore, we modelled the subsequent joint evolution of artisanal and industrial fisheries. Overall results confirmed that Mediterranean fishery time series were highly diverse along their dynamics and this heterogeneity was persistent throughout the time. In addition, results highlighted a positive relationship between the two types of fishing.


title: How to manage the post pandemic opening? A Pontryagin Maximum Principle approach

abstract: The COVID-19 pandemic has completely disrupted the operation of our societies. Its elusive transmission process, characterized by an unusually long incubation period, as well as a high contagion capacity, has forced many countries to take quarantine and social isolation measures that conspire against the performance of national economies. This situation confronts decision makers in different countries with the alternative of reopening the economies, thus facing the unpredictable cost of a rebound of the infection. This work tries to offer an initial theoretical framework to handle this alternative.


title: High-quality electrical transport using scalable CVD graphene

abstract: Producing and manipulating graphene on fab-compatible scale, while maintaining its remarkable carrier mobility, is key to finalize its technological application. We show that a large-scale approach (chemical vapor deposition on Cu followed by polymer-mediated semi-dry transfer) yields single-layer graphene crystals fully comparable, in terms of electronic transport, to micro-mechanically exfoliated flakes. hBN is used to encapsulate the graphene crystals $-$ without taking part to their detachment from the growth catalyst $-$ and study their intrinsic properties in field-effect devices. At room temperature, the electron-phonon coupling sets the mobility to $\sim1.3 \times10^5$ cm$^2$V$^{-1}$s$^{-1}$ at $\sim10^{11}$ cm$^{-2}$ concentration. At cryogenic temperatures, the mobility ($ > 6\times10^5$ cm$^2$V$^{-1}$s$^{-1}$ at $\sim10^{11}$ cm$^{-2}$) is limited by the devices' physical edges, and charge fluctuations $ < 7\times10^9$ cm$^{-2}$ are detected. Under perpendicular magnetic fields, we observe early onset of Landau quantization ($B\sim50$ mT) and signatures of electronic correlation, including the fractional quantum Hall effect.


title: Flattening of the tokamak current profile by a fast magnetic reconnection with implications for the solar corona

abstract: During tokamak disruptions the profile of the net parallel current is observed to flatten on a time scale that is so fast that it must be due to a fast magnetic reconnection. After a fast magnetic reconnection has broken magnetic surfaces, a single magnetic field line covers an entire volume and not just a magnetic surface. The current profile, given by $K\equiv\mu_0j_{||}/B$, relaxes to a constant within that volume by Alfv\'en waves propagating along the chaotic magnetic field lines. The time scale for this relaxation determines the commonly observed disruption phenomena of a current spike and a sudden drop in the plasma internal inductance. An efficient method for studying this relaxation is derived, which allows a better understanding of the information encoded in the current spike and the associated sudden drop in the plasma internal inductance. Implications for coronal heating are also discussed.


title: Reducing the impact of radioactivity on quantum circuits in a deep-underground facility

abstract: As quantum coherence times of superconducting circuits have increased from nanoseconds to hundreds of microseconds, they are currently one of the leading platforms for quantum information processing. However, coherence needs to further improve by orders of magnitude to reduce the prohibitive hardware overhead of current error correction schemes. Reaching this goal hinges on reducing the density of broken Cooper pairs, so-called quasiparticles. Here, we show that environmental radioactivity is a significant source of nonequilibrium quasiparticles. Moreover, ionizing radiation introduces time-correlated quasiparticle bursts in resonators on the same chip, further complicating quantum error correction. Operating in a deep-underground lead-shielded cryostat decreases the quasiparticle burst rate by a factor fifty and reduces dissipation up to a factor four, showcasing the importance of radiation abatement in future solid-state quantum hardware.


title: Convergence analysis of the scaled boundary finite element method for the Laplace equation

abstract: The scaled boundary finite element method (SBFEM) is a relatively recent boundary element method that allows the approximation of solutions to PDEs without the need of a fundamental solution. A theoretical framework for the convergence analysis of SBFEM is proposed here. This is achieved by defining a space of semi-discrete functions and constructing an interpolation operator onto this space. We prove error estimates for this interpolation operator and show that optimal convergence to the solution can be obtained in SBFEM. These theoretical results are backed by a numerical example.


title: Transverse-momentum-dependent gluon distribution functions in a spectator model

abstract: We present a model calculation of transverse-momentum-dependent distributions (TMDs) of gluons in the nucleon. The model is based on the assumption that a nucleon can emit a gluon, and what remains after the emission is treated as a single spectator particle. This spectator particle is considered to be on-shell, but its mass is allowed to take a continuous range of values, described by a spectral function. The nucleon-gluon-spectator coupling is described by an effective vertex containing two form factors. We fix the model parameters to obtain the best agreement with collinear gluon distributions extracted from global fits. We study the tomography in momentum space of gluons inside nucleons for various combinations of their polarizations. These can be used to make predictions of observables relevant for gluon TMD studies at current and future collider facilities.


title: In Silico Trial to test COVID-19 candidate vaccines: a case study with UISS platform

abstract: SARS-CoV-2 is a severe respiratory infection that infects humans. Its outburst entitled it as a pandemic emergence. To get a grip on this, outbreak specific preventive and therapeutic interventions are urgently needed. It must be said that, until now, there are no existing vaccines for coronaviruses. To promptly and rapidly respond to pandemic events, the application of in silico trials can be used for designing and testing medicines against SARS-CoV-2 and speed-up the vaccine discovery pipeline, predicting any therapeutic failure and minimizing undesired effects. Here, we present an in silico platform that showed to be in very good agreement with the latest literature in predicting SARS- CoV-2 dynamics and related immune system host response. Moreover, it has been used to predict the outcome of one of the latest suggested approach to design an effective vaccine, based on monoclonal antibody. UISS is then potentially ready to be used as an in silico trial platform to predict the outcome of vaccination strategy against SARS-CoV-2.


title: Convex Bodies with affinely equivalent projections and affine bodies of revolution

abstract: Abstract In this paper, we study affine bodies of revolution. This will allow us to prove that a convex body all whose orthogonal $n$-projections are affinely equivalent is an ellipsoid, provided $n\equiv 0,1, 2$ mod $4$, $n>1$ with the possible exemption of $n=133$. Our proof uses convex geometry and topology of compact Lie groups. AMS classification subject: 22E10, 52A05


title: NTIRE 2020 Challenge on Video Quality Mapping: Methods and Results

abstract: This paper reviews the NTIRE 2020 challenge on video quality mapping (VQM), which addresses the issues of quality mapping from source video domain to target video domain. The challenge includes both a supervised track (track 1) and a weakly-supervised track (track 2) for two benchmark datasets. In particular, track 1 offers a new Internet video benchmark, requiring algorithms to learn the map from more compressed videos to less compressed videos in a supervised training manner. In track 2, algorithms are required to learn the quality mapping from one device to another when their quality varies substantially and weakly-aligned video pairs are available. For track 1, in total 7 teams competed in the final test phase, demonstrating novel and effective solutions to the problem. For track 2, some existing methods are evaluated, showing promising solutions to the weakly-supervised video quality mapping problem.


title: Dark Matter with Light and Ultralight St\"uckelberg Axions

abstract: Scenarios with axion-like particles of variable masses, from the milli-eV range to ultralight, can be generated in a natural way in the context of effective field theory models with anomalous $U(1)$ gauge symmetries. They include dimension-5 operators which define a coupling of the axion to the gauge anomaly. They provide a realization, in the domain of ordinary gauge theories, of models for such particles, which evade the usual mass/coupling constraint of ordinary (Peccei-Quinn) axions and are natural dark matter candidates. As an illustration of these models, we present an overview of two of these scenarios. One of them is built around the next-to-minimal MSSM (NMSSM), a model called the USSM-A, which illustrates how the St\"uckelberg supermultiplet can be used to generate two dark matter candidates, a neutralino containing an axino component coming from the St\"uckelberg sector, plus the axion ($\textrm{Im}\, b$). The real component of the complex St\"uckelberg field carries dilaton-like ($\textrm{Re}\, b\, FF$) interactions. In a second model, non-supersymmetric, the St\"uckelberg scale is raised up to the GUT epoch. In this case the axion mass can be ultralight ($\sim 10^{-20}$ eV). The periodic potential generated at the GUT phase transition and the corresponding oscillations are related to a particle whose De Broglie wavelength can be sub-galactic. A similar analysis is also possible for the supersymmetric scenario.


title: Statistics of the Number of Records for Random Walks and L\'evy Flights on a ${1D}$ Lattice

abstract: We study the statistics of the number of records $R_n$ for a symmetric, $n$-step, discrete jump process on a $1D$ lattice. At a given step, the walker can jump by arbitrary lattice units drawn from a given symmetric probability distribution. This process includes, as a special case, the standard nearest neighbor lattice random walk. We derive explicitly the generating function of the distribution $P(R_n)$ of the number of records, valid for arbitrary discrete jump distributions. As a byproduct, we provide a relatively simple proof of the generalized Sparre Andersen theorem for the survival probability of a random walk on a line, with discrete or continuous jump distributions. For the discrete jump process, we then derive the asymptotic large $n$ behavior of $P(R_n)$ as well as of the average number of records $E(R_n)$. We show that unlike the case of random walks with symmetric and continuous jump distributions where the record statistics is strongly universal (i.e., independent of the jump distribution for all $n$), the record statistics for lattice walks depends on the jump distribution for any fixed $n$. However, in the large $n$ limit, we show that the distribution of the scaled record number $R_n/E(R_n)$ approaches a universal, half-Gaussian form for any discrete jump process. The dependence on the jump distribution enters only through the scale factor $E(R_n)$, which we also compute in the large $n$ limit for arbitrary jump distributions. We present explicit results for a few examples and provide numerical checks of our analytical predictions.

categories: cond-mat.stat-mech math-ph math.MP math.PR


title: Using posterior predictive distributions to analyse epidemic models: COVID-19 in Mexico City

abstract: Epidemiological models contain a set of parameters that must be adjusted based on available observations. Once a model has been calibrated, it can be used as a forecasting tool to make predictions and to evaluate contingency plans. It is customary to employ only point estimators for such predictions. However, some models may fit the same data reasonably well for a broad range of parameter values, and this flexibility means that predictions stemming from such models will vary widely, depending on the particular parameter values employed within the range that give a good fit. When data are poor or incomplete, model uncertainty widens further. A way to circumvent this problem is to use Bayesian statistics to incorporate observations and use the full range of parameter estimates contained in the posterior distribution to adjust for uncertainties in model predictions. Specifically, given the epidemiological model and a probability distribution for observations, we use the posterior distribution of model parameters to generate all possible epidemiological curves via the posterior predictive distribution. From the envelope of all curves one can extract the worst-case scenario and study the impact of implementing contingency plans according to this assessment. We apply this approach to the potential evolution of COVID-19 in Mexico City and assess whether contingency plans are being successful and whether the epidemiological curve has flattened.


title: Code-switching patterns can be an effective route to improve performance of downstream NLP applications: A case study of humour, sarcasm and hate speech detection

abstract: In this paper we demonstrate how code-switching patterns can be utilised to improve various downstream NLP applications. In particular, we encode different switching features to improve humour, sarcasm and hate speech detection tasks. We believe that this simple linguistic observation can also be potentially helpful in improving other similar NLP applications.


title: On building a cluster watch-list for identifying strongly lensed supernovae, gravitational waves and kilonovae

abstract: Motivated by discovering strongly-lensed supernovae, gravitational waves, and kilonovae in the 2020s, we investigate whether to build a watch-list of clusters based on observed cluster properties (i.e. lens-plane selection) or on the detectability of strongly-lensed background galaxies (i.e. source-plane selection). First, we estimate the fraction of high-redshift transient progenitors that reside in galaxies that are themselves too faint to be detected as being strongly-lensed. We find $\sim15-50$ per cent of transient progenitors reside in $z = 1-2$ galaxies too faint to be detected in surveys that reach ${\rm AB}\simeq23$, such as the Dark Energy Survey. This falls to $\ls10$ per cent at depths that will be probed by early data releases of LSST (${\rm AB}\simeq25$). Second, we estimate a conservative lower limit on the fraction of strong lensing clusters that will be missed by magnitude limited searches for multiply-imaged galaxies and giant arcs due to the faintness of such images. We find that DES-like surveys will miss $\sim75$ per cent of $10^{15}$M$_\odot$ strong lensing clusters, rising to $\sim100$ per cent of $10^{14}$M$_\odot$ clusters. Deeper surveys, such as LSST, will miss $\sim40$ per cent at $10^{15}$M$_\odot$, and $\sim95$ per cent at $10^{14}$M$_\odot$. Our results motivate building a cluster watch-list for strongly-lensed transients that includes those found by lens-plane selection.


title: NOMA Visible Light Communication System with Angle Diversity Receivers

abstract: In this paper, a non-orthogonal multiple access (NOMA) visible light communication (VLC) system is investigated. The system uses angle diversity receivers (ADRs) to provide high data rates. The ADR has 4 branches, each directed to a different direction. An 8m x 4m sized room is modelled to study the resource allocation to users according to their channel conditions to maximize the data rate. The results show that using ADRs improves the data rate by an average of 35% compared to a system using wide FOV receivers.


title: Reducing Uncertainty by Fusing Dynamic Occupancy Grid Maps in a Cloud-based Collective Environment Model

abstract: Accurate environment perception is essential for automated vehicles. Since occlusions and inaccuracies regularly occur, the exchange and combination of perception data of multiple vehicles seems promising. This paper describes a method to combine perception data of automated and connected vehicles in the form of evidential Dynamic Occupany Grid Maps (DOGMas) in a cloud-based system. This system is called the Collective Environment Model and is part of the cloud system developed in the project UNICARagil. The presented concept extends existing approaches that fuse evidential grid maps representing static environments of a single vehicle to evidential grid maps computed by multiple vehicles in dynamic environments. The developed fusion process additionally incorporates self-reported data provided by connected vehicles instead of only relying on perception data. We show that the uncertainty in a DOGMa described by Shannon entropy as well as the uncertainty described by a non-specificity measure can be reduced. This enables automated and connected vehicles to behave in ways not before possible due to unknown but relevant information about the environment.


title: Parameters Estimation from the 21 cm signal using Variational Inference

abstract: Upcoming experiments such as Hydrogen Epoch of Reionization Array (HERA) and Square Kilometre Array (SKA) are intended to measure the 21cm signal over a wide range of redshifts, representing an incredible opportunity in advancing our understanding about the nature of cosmic Reionization. At the same time these kind of experiments will present new challenges in processing the extensive amount of data generated, calling for the development of automated methods capable of precisely estimating physical parameters and their uncertainties. In this paper we employ Variational Inference, and in particular Bayesian Neural Networks, as an alternative to MCMC in 21 cm observations to report credible estimations for cosmological and astrophysical parameters and assess the correlations among them.


title: When Votes Change and Committees Should (Not)

abstract: Electing a single committee of a small size is a classical and well-understood voting situation. Being interested in a sequence of committees, we introduce and study two time-dependent multistage models based on simple Plurality voting. Therein, we are given a sequence of voting profiles (stages) over the same set of agents and candidates, and our task is to find a small committee for each stage of high score. In the conservative model we additionally require that any two consecutive committees have a small symmetric difference. Analogously, in the revolutionary model we require large symmetric differences. We prove both models to be NP-hard even for a constant number of agents, and, based on this, initiate a parameterized complexity analysis for the most natural parameters and combinations thereof. Among other results, we prove both models to be in XP yet W[1]-hard regarding the number of stages, and that being revolutionary seems to be "easier" than being conservative: If the (upper- resp. lower-) bound on the size of symmetric differences is constant, the conservative model remains NP-hard while the revolutionary model becomes polynomial-time solvable.


title: Regional Strategic Sensors Characterizations

abstract: In this paper, a linear infinite dimensional distributed systems in a Hilbert space has been discussed and analyzed where the dynamics of system is governed by strongly continuous semi-groups. The characterizations of regional strategic sensors have been given and tackled for different cases of regional observability. Furthermore, the results so obtained are applied to two-dimensional systems. Various cases of sensors are considered and analyzed. Also, the authors show that, the existent of a sensor for the diffusion system is not strategic in the usual sense, but it may be regionally strategic of this system.


title: Bayesian Inference for Johnson's SB and Weibull distributions

abstract: The four-parameter Johnson's SB (JSB) and three-parameter Weibull distributions have received much attention in the field of forestry for characterizing diameters at breast height (DBH). In this work, we suggest the Bayesian method for estimating parameters of the JBS distribution. The maximum likelihood approach uses iterative methods such as Newton-Raphson (NR) algorithm for maximizing the logarithm of the likelihood function. But there is no guarantee that the NR method converges. This fact that the NR method for estimating the parameters of the JSB distribution sometimes fails to converge was verified through simulation in this study. Further, it was shown that the Bayesian estimators presented in this work were robust with respect to the initial values and estimate the parameters of the JSB distribution efficiently. The performance of the JSB and three-parameter Weibull distributions was compared in a Bayesian paradigm when these models were fitted to DBH data of three plots that randomly selected from a study established in 107 plots of mixed-age ponderosa pine (Pinus ponderosa Dougl. ex Laws.) with scattered western junipers at the Malheur National Forest in south end of the Blue Mountains near Burns, Oregon, USA. Bayesian paradigm demonstrated that JBS was superior model than the three-parameter Weibull for characterizing the DBH distribution when these models were fitted to the DBH data of the three plots.


title: Effective temperatures in inhomogeneous passive and active bidimensional Brownian particle systems

abstract: We study the stationary dynamics of an active interacting Brownian particle system. We measure the violations of the fluctuation dissipation theorem, and the corresponding effective temperature, in a locally resolved way. Quite naturally, in the homogeneous phases the diffusive properties and effective temperature are also homogeneous. Instead, in the inhomogeneous phases (close to equilibrium and within the MIPS sector) the particles can be separated in two groups with different diffusion properties and effective temperatures. Notably, at fixed activity strength the effective temperatures in the two phases remain distinct and approximately constant within the MIPS region, with values corresponding to the ones of the whole system at the boundaries of this sector of the phase diagram. We complement the study of the globally averaged properties with the theoretical and numerical characterization of the fluctuation distributions of the single particle diffusion, linear response, and effective temperature in the homogeneous and inhomogeneous phases. We also distinguish the behavior of the (time-delayed) effective temperature from the (instantaneous) kinetic temperature, showing that the former is independent on the friction coefficient.


title: Resonating Experiences of Self and Others enabled by a Tangible Somaesthetic Design

abstract: Digitalization is penetrating every aspect of everyday life including a human's heart beating, which can easily be sensed by wearable sensors and displayed for others to see, feel, and potentially "bodily resonate" with. Previous work in studying human interactions and interaction designs with physiological data, such as a heart's pulse rate, have argued that feeding it back to the users may, for example support users' mindfulness and self-awareness during various everyday activities and ultimately support their wellbeing. Inspired by Somaesthetics as a discipline, which focuses on an appreciation of the living body's role in all our experiences, we designed and explored mobile tangible heart beat displays, which enable rich forms of bodily experiencing oneself and others in social proximity. In this paper, we first report on the design process of tangible heart displays and then present results of a field study with 30 pairs of participants. Participants were asked to use the tangible heart displays during watching movies together and report their experience in three different heart display conditions (i.e., displaying their own heart beat, their partner's heart beat, and watching a movie without a heart display). We found, for example that participants reported significant effects in experiencing sensory immersion when they felt their own heart beats compared to the condition without any heart beat display, and that feeling their partner's heart beats resulted in significant effects on social experience. We refer to resonance theory to discuss the results, highlighting the potential of how ubiquitous technology could utilize physiological data to provide resonance in a modern society facing social acceleration.


title: Generalized Planning With Deep Reinforcement Learning

abstract: A hallmark of intelligence is the ability to deduce general principles from examples, which are correct beyond the range of those observed. Generalized Planning deals with finding such principles for a class of planning problems, so that principles discovered using small instances of a domain can be used to solve much larger instances of the same domain. In this work we study the use of Deep Reinforcement Learning and Graph Neural Networks to learn such generalized policies and demonstrate that they can generalize to instances that are orders of magnitude larger than those they were trained on.


title: Tilting modules, dominant dimensions and Brauer-Schur-Weyl duality

abstract: Let $A$ be a standardly stratified algebra over a field $K$ and $T$ a tilting module over $A$. Let $\Lambda^+$ be an indexing set of all simple modules in $A\lmod$. We show that if there is an integer $r\in\N$ such that for any $\lambda\in\Lambda^+$, there is an embedding $\Delta(\lambda)\hookrightarrow T^{\oplus r}$ as well as an epimorphism $T^{\oplus r}\twoheadrightarrow\overline{\nabla}(\lambda)$ as $A$-modules, then $T$ is a faithful $A$-module and $A$ has the double centraliser property with respect to $T$. As applications, we prove that if $A$ is quasi-hereditary with a simple preserving duality and $T$ a given faithful tilting $A$-module, then $A$ has the double centralizer property with respect to $T$. This provides a simple and useful criterion which can be applied in many situations in algebraic Lie theory. We affirmatively answer a question of Mazorchuk and Stroppel by proving the existence of a unique minimal basic tilting module $T$ over $A$ for which $A=\End_{\End_A(T)}(T)$. We also establish a Schur-Weyl duality between the symplectic Schur algebra $S^{sy}(m,n)$ and $\bb_{n}/\mathfrak{B}_{n}^{(f)}$ on $V^{\otimes n}/V^{\otimes n}\mathfrak{B}_{n}^{(f)}$ when $\cha K>\min\{n-f+m,n\}$, where $V$ is a $2m$-dimensional symplectic space over $K$, $\mathfrak{B}_{n}^{(f)}$ is the two-sided ideal of the Brauer algebra $\bb_{n}(-2m)$ generated by $e_1e_3\cdots e_{2f-1}$ with $1\leq f\leq [\frac{n}{2}]$.


title: Modeling Self-Healing Behavior of Vitrimers using Molecular Dynamics with Dynamic Cross-Linking Capability

abstract: Vitrimers are a special class of polymers that undergo dynamic cross-linking under thermal stimuli. Their ability to exchange covalent bonds can be harnessed to mitigate damage in a composite or to achieve recyclable aerospace composites. This work addresses the primary challenge of modeling dynamic cross-linking reactions in vitrimers during thermomechanical loading. Dynamic bond exchange reaction probability change during heating and its effect on dilatometric and mechanical response are simulated for the first time in large scale molecular dynamics simulations. Healing of damage under thermal cycling is computed with mechanical properties predicted before and after self--healing.


title: Precoder Design and Statistical Power Allocation for MIMO-NOMA via User-Assisted Simultaneous Diagonalization

abstract: In this paper, we investigate the downlink precoder design for two-user power-domain multiple-input multiple-output (MIMO) non-orthogonal multiple access (NOMA). We propose a novel user-assisted (UA) simultaneous diagonalization (SD) based MIMO-NOMA scheme that achieves SD of the MIMO channels of both users through a combination of precoder design and low-complexity self-interference cancellation at the users, thereby considerably lowering the overall decoding complexity compared to joint decoding. The achievable ergodic user rates of the proposed scheme are analyzed for Rayleigh fading channels based on a finite-size random matrix theory framework, which is further exploited to develop a statistical power allocation algorithm. Simulation and numerical results show that the proposed UA-SD MIMO-NOMA scheme significantly outperforms orthogonal multiple access and a benchmark precoder design performing SD via generalized singular value decomposition in terms of the achievable ergodic rate region for most user rates. The ergodic rate region is further enhanced by a hybrid scheme which performs time sharing between the proposed UA-SD MIMO-NOMA scheme and single-user MIMO.


title: Properties of a recent quantum extension of the Kruskal geometry

abstract: Recently it was shown that, in an effective description motivated by loop quantum gravity, singularities of the Kruskal space-time are naturally resolved [1,2]. In this note we explore a few properties of this quantum corrected effective metric. In particular, we (i) calculate the Hawking temperature associated with the horizon of the effective geometry and show that the quantum correction to the temperature is completely negligible for macroscopic black holes, just as one would hope; (ii) discuss the subtleties associated with the asymptotic properties of the space-time metric, and show that the metric is asymptotically flat in a precise sense; (iii) analyze the asymptotic fall-off of curvature; and, (iv) show that the ADM energy is well-defined (and agrees with that determined by the horizon area), even though the curvature falls off less rapidly than in the standard asymptotically flat context.


title: Testing Compilers for Programmable Switches Through Switch Hardware Simulation

abstract: Programmable switches have emerged as powerful and flexible alternatives to fixed-function forwarding devices. But because of the unique hardware constraints of network switches, the design and implementation of compilers targeting these devices is tedious and error prone. Despite the important role that compilers play in software development, there is a dearth of tools for testing compilers for programmable network devices. We present Druzhba, a programmable switch simulator used for testing compilers targeting programmable packet-processing substrates. We show that we can model the low-level behavior of a switch's programmable hardware. We further show how our machine model can be used by compiler developers to target Druzhba as a compiler backend. Generated machine code programs are fed into Druzhba and tested using a fuzzing-based approach that allows compiler developers to test the correctness of their compilers. Using a program-synthesis-based compiler as a case study, we demonstrate how Druzhba has been successful in testing compiler-generated machine code for our simulated switch pipeline instruction set.


title: Solutions for nonlinear Fokker-Planck equations with measures as initial data and McKean-Vlasov equations

abstract: One proves the existence and uniqueness of a generalized (mild) solution for the nonlinear Fokker-Planck equation (FPE) \begin{align*} &u_t-\Delta (\beta(u))+{\mathrm{ div}}(D(x)b(u)u)=0, \quad t\geq0,\ x\in\mathbb{R}^d,\ d\ne2, \\ &u(0,\cdot)=u_0,\mbox{in }\mathbb{R}^d, \end{align*} where $u_0\in L^1(\mathbb{R}^d)$, $\beta\in C^2(\mathbb{R})$ is a nondecreasing function, $b\in C^1$, bounded, $b\ge0$, $D\in {L^\infty}(\mathbb{R}^d;\mathbb{R}^d)$, ${\rm div}\,D\in L^2(\mathbb{R}^d)+L^\infty(\mathbb{R}^d),$ with ${({\rm div}\, D)^-}\in L^\infty(\mathbb{R}^d)$, $\beta$ strictly increasing, if $b$ is not constant. Moreover, $t\to u(t,u_0)$ is a semigroup of contractions in $L^1(\mathbb{R}^d)$, which leaves invariant the set of probability density functions in $\mathbb{R}^d$. If ${\rm div}\,D\ge0$, $\beta'(r)\ge a|r|^{\alpha-1}$, and $|\beta(r)|\le C r^\alpha$, $\alpha\ge1,$ $d\ge3$, then $|u(t)|_{L^\infty}\le Ct^{-\frac d{d+(\alpha-1)d}}\ |u_0|^{\frac2{2+(m-1)d}},$ $t>0$, and, if $D\in L^2(\mathbb{R}^d;\mathbb{R}^d)$, the existence extends to initial data $u_0$ in the space $\mathcal{M}_b$ of bounded measures in $\mathbb{R}^d$. As a consequence for arbitrary initial laws, we obtain weak solutions to a class of McKean-Vlasov SDEs with coefficients which have singular dependence on the time marginal laws.


title: Thermodynamic uncertainty relation for time-dependent driving

abstract: Thermodynamic uncertainty relations yield a lower bound on entropy production in terms of the mean and fluctuations of a current. We derive their general form for systems under arbitrary time-dependent driving from arbitrary initial states and extend these relations beyond currents to state variables. The quality of the bound is discussed for various types of observables for an interacting pair of colloidal particles in a moving laser trap and for the dynamical unfolding of a small protein. Since the input for evaluating these bounds does not require specific knowledge of the system or its coupling to the time-dependent control, they should become widely applicable tools for thermodynamic inference in time-dependently driven systems.


title: Adversarial Training against Location-Optimized Adversarial Patches

abstract: Deep neural networks have been shown to be susceptible to adversarial examples -- small, imperceptible changes constructed to cause mis-classification in otherwise highly accurate image classifiers. As a practical alternative, recent work proposed so-called adversarial patches: clearly visible, but adversarially crafted rectangular patches in images. These patches can easily be printed and applied in the physical world. While defenses against imperceptible adversarial examples have been studied extensively, robustness against adversarial patches is poorly understood. In this work, we first devise a practical approach to obtain adversarial patches while actively optimizing their location within the image. Then, we apply adversarial training on these location-optimized adversarial patches and demonstrate significantly improved robustness on CIFAR10 and GTSRB. Additionally, in contrast to adversarial training on imperceptible adversarial examples, our adversarial patch training does not reduce accuracy.


title: Superconducting proximity effect in a transparent van der Waals superconductor-metal junction

abstract: We report on Andreev reflections at clean NbSe2-bilayer graphene junctions. The high transparency of the junction, which manifests as a large conductance enhancement of up to 1.8, enables us to see clear evidence of a proximity-induced superconducting gap in bilayer graphene and two Andreev reflections through a vertical NbSe2-graphene and a lateral graphene-graphene junction respectively. Quantum transport simulations capture the complexity of the experimental data and illuminate the impact of various microscopic parameters on the transmission of the junction. Our work establishes the practice and understanding of an all-van-der-Waals, high-performance superconducting junction. The realization of a highly transparent proximized graphene-graphene junction opens up possibilities to engineer emergent quantum phenomena.


title: Multi-interactive Dual-decoder for RGB-thermal Salient Object Detection

abstract: RGB-thermal salient object detection (SOD) aims to segment the common prominent regions of visible image and corresponding thermal infrared image that we call it RGBT SOD. Existing methods don't fully explore and exploit the potentials of complementarity of different modalities and multi-type cues of image contents, which play a vital role in achieving accurate results. In this paper, we propose a multi-interactive dual-decoder to mine and model the multi-type interactions for accurate RGBT SOD. In specific, we first encode two modalities into multi-level multi-modal feature representations. Then, we design a novel dual-decoder to conduct the interactions of multi-level features, two modalities and global contexts. With these interactions, our method works well in diversely challenging scenarios even in the presence of invalid modality. Finally, we carry out extensive experiments on public RGBT and RGBD SOD datasets, and the results show that the proposed method achieves the outstanding performance against state-of-the-art algorithms. The source code has been released at:https://github.com/lz118/Multi-interactive-Dual-decoder.


title: Laplacian Spectra of Comaximal Graph of $\mathbb{Z}_{n}$

abstract: This article focuses on finding the eigenvalues of the Laplacian matrix of the comaximal graph $\Gamma(\mathbb Z_n)$ of the ring $\mathbb Z_n$ for $n> 2$. We determine the eigenvalues of $\Gamma(\mathbb Z_n)$ for various $n$ and also provide a procedure to find the eigenvalues of $\Gamma(\mathbb Z_n)$ for any $n> 2$. We show that $\Gamma(\mathbb Z_n)$ is Laplacian Integral for $n=p^\alpha q^\beta$ where $p,q$ are primes and $\alpha, \beta$ are non-negative integers. The algebraic and vertex connectivity of $\Gamma(\mathbb Z_n)$ have been shown to be equal for all $n> 2$. An upper bound on the second largest eigenvalue of $\Gamma(\mathbb Z_n)$ has been obtained and a necessary and sufficient condition for its equality has also been determined. Finally we discuss the multiplicity of the spectral radius and the multiplicity of the algebraic connectivity of $\Gamma(\mathbb Z_n)$. Some problems have been discussed at the end of this article for further research.


title: Construction of Modern Robust Nodal Discontinuous Galerkin Spectral Element Methods for the Compressible Navier-Stokes Equations

abstract: Discontinuous Galerkin (DG) methods have a long history in computational physics and engineering to approximate solutions of partial differential equations due to their high-order accuracy and geometric flexibility. However, DG is not perfect and there remain some issues. Concerning robustness, DG has undergone an extensive transformation over the past seven years into its modern form that provides statements on solution boundedness for linear and nonlinear problems. This chapter takes a constructive approach to introduce a modern incarnation of the DG spectral element method for the compressible Navier-Stokes equations in a three-dimensional curvilinear context. The groundwork of the numerical scheme comes from classic principles of spectral methods including polynomial approximations and Gauss-type quadratures. We identify aliasing as one underlying cause of the robustness issues for classical DG spectral methods. Removing said aliasing errors requires a particular differentiation matrix and careful discretization of the advective flux terms in the governing equations.


title: A neural network model for solvency calculations in life insurance

abstract: Insurance companies make extensive use of Monte Carlo simulations in their capital and solvency models. To overcome the computational problems associated with Monte Carlo simulations, most large life insurance companies use proxy models such as replicating portfolios. In this paper, we present an example based on a variable annuity guarantee, showing the main challenges faced by practitioners in the construction of replicating portfolios: the feature engineering step and subsequent basis function selection problem. We describe how neural networks can be used as a proxy model and how to apply risk-neutral pricing on a neural network to integrate such a model into a market risk framework. The proposed model naturally solves the feature engineering and feature selection problems of replicating portfolios.


title: Tuning Rules for a Class of Port-Hamiltonian Mechanical Systems

abstract: In this extended abstract, we propose a tuning approach for nonlinear mechanical systems to modify the behavior of the closed-loop system, where we are particularly interested in attenuating oscillations from the transient response. Towards this end, we inject damping into the system, and we provide two tuning methods to select the gains that are appropriate for our purposes. Furthermore, we apply these tuning rules to a 2DoF planar manipulator and present its simulation results.


title: Model studies of fluctuations in the background for jets in heavy ion collisions

abstract: Jets produced in high energy heavy ion collisions are quenched by the quark gluon plasma. Measurements of these jets are influenced by the methods used to suppress and subtract the large, fluctuating background and the assumptions inherent in these methods. We compare the measurements of the background in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV by the ALICE Collaboration $\underline{(B. Abelev \textit{et al.}, J. High Energy Phys. \mathbf{2012}, 053)}$ to calculations in TennGen (a data-driven random background generator) and PYTHIA Angantyr. A detailed understanding of the width of these fluctuations is important for reducing uncertainties due to unfolding and extending measurements to lower momenta and larger resolution parameters. The standard deviation of the energy in random cones in TennGen is approximately in agreement with the form predicted in the ALICE paper, with deviations of 1--6%. The standard deviation of energy in random cones in Angantyr exceeds the same predictions by approximately 13%. Deviations in both models can be explained by the assumption that the single-particle $d^2N/dy$ $dp_T$ is a gamma distribution in the derivation of the prediction, whereas the model uses a different distribution. This indicates that model comparisons are potentially sensitive to the treatment of the background. We demonstrate that unfolding methods used to remove background fluctuations from jets can affect the comparisons between models and data, $\textit{even in the absence of detector effects}$. Our findings suggest the need to more carefully consider methods for comparing simulations and data.


title: Finite Sample Smeariness of Fr\'echet Means and Application to Climate

abstract: Fr\'echet means on non-Euclidean spaces may exhibit nonstandard asymptotic rates rendering quantile-based asymptotic inference inapplicable. We show here that this affects, among others, all circular distributions whose support exceeds a half circle. We exhaustively describe this phenomenon and introduce a new concept which we call finite samples smeariness (FSS). In the presence of FSS, it turns out that quantile-based tests for equality of Fr\'echet means systematically feature effective levels higher than their nominal level which perseveres asymptotically in case of Type I FSS. In contrast, suitable bootstrap-based tests correct for FSS and asymptotically attain the correct level. For illustration of the relevance of FSS in real data, we apply our method to directional wind data from two European cities. It turns out that quantile based tests, not correcting for FSS, find a multitude of significant wind changes. This multitude condenses to a few years featuring significant wind changes, when our bootstrap tests are applied, correcting for FSS.


title: Almost Kenmotsu metric as Ricci-Yamabe soliton

abstract: The object of the present paper is to characterize two classes of almost Kenmotsu manifolds admitting Ricci-Yamabe soliton. It is shown that a $(k,\mu)'$-almost Kenmotsu manifold admitting a Ricci-Yamabe soliton or gradient Ricci-Yamabe soliton is locally isometric to the Riemannian product $\mathbb{H}^{n+1}(-4) \times \mathbb{R}^n$. For the later case, the potential vector field is pointwise collinear with the Reeb vector field. Also, a $(k,\mu)$-almost Kenmotsu manifold admitting certain Ricci-Yamabe soliton with the curvature property $Q \cdot P = 0$ is locally isometric to the hyperbolic space $\mathbb{H}^{2n+1}(-1)$ and the non-existense of the curvature property $Q \cdot R = 0$ is proved.


title: Study of single crystalline SrAgSb and SrAuSb semimetals

abstract: Given renewed interest in the electronic properties of semimetallic compounds with varying degrees of spin orbit coupling we have grown single crystals of SrAgSb and SrAuSb, measured their temperature and field dependent electrical resistivity and magnetization and performed density functional theory (DFT) band structure calculations. Magnetization measurements are consistent with a diamagnetic host with a small amount of local moment bearing impurities. Although the residual resistivity ratio (RRR) for all samples studied was relatively low, ranging between 2.4 and 3.4, the compounds had non-saturating magnetoresistance (MR), reaching values of $\sim$ 17% and $\sim$ 70% at 4 K and 9 T for SrAgSb and SrAuSb respectively. Band structure calculations, using the experimentally determined Wyckoff positions for the Sr, Ag/Au, and Sb atoms, show that whereas SrAgSb is a topologically trivial, but compensated, semimetal; SrAuSb is a topologically non-trivial, Dirac semimetal.


title: Neural CRF Model for Sentence Alignment in Text Simplification

abstract: The success of a text simplification system heavily depends on the quality and quantity of complex-simple sentence pairs in the training corpus, which are extracted by aligning sentences between parallel articles. To evaluate and improve sentence alignment quality, we create two manually annotated sentence-aligned datasets from two commonly used text simplification corpora, Newsela and Wikipedia. We propose a novel neural CRF alignment model which not only leverages the sequential nature of sentences in parallel documents but also utilizes a neural sentence pair model to capture semantic similarity. Experiments demonstrate that our proposed approach outperforms all the previous work on monolingual sentence alignment task by more than 5 points in F1. We apply our CRF aligner to construct two new text simplification datasets, Newsela-Auto and Wiki-Auto, which are much larger and of better quality compared to the existing datasets. A Transformer-based seq2seq model trained on our datasets establishes a new state-of-the-art for text simplification in both automatic and human evaluation.


title: Digraph of Senegal s local languages: issues, challenges and prospects of their transliteration

abstract: The local languages in Senegal, like those of West African countries in general, are written based on two alphabets: supplemented Arabic alphabet (called Ajami) and Latin alphabet. Each writing has its own applications. Ajami writing is generally used by people educated in Koranic schools for communication, business, literature (religious texts, poetry, etc.), traditional religious medicine, etc. Writing with Latin characters is used for localization of ICT (Web, dictionaries, Windows and Google tools translated in Wolof, etc.), the translation of legal texts (commercial code and constitution translated in Wolof) and religious ones (Quran and Bible in Wolof), book edition, etc. To facilitate both populations general access to knowledge, it is useful to set up transliteration tools between these two scriptures. This work falls within the framework of the implementation of project for a collaborative online dictionary Wolof (Nguer E. M., Khoule M, Thiam M. N., Mbaye B. T., Thiare O., Cisse M. T., Mangeot M. 2014), which will involve people using Ajami writing. Our goal will consist, on the one hand in raising the issues related to the transliteration and the challenges that this will raise, and on the other one, presenting the perspectives.


title: Modeling the 10 September 2017 solar energetic particle event using the iPATH model

abstract: On September 10 2017, a fast coronal mass ejection (CME) erupted from the active region (AR) 12673, leading to a ground level enhancement (GLE) event at Earth. Using the 2D improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, we model the large solar energetic particle (SEP) event of 10 September 2017 observed at Earth, Mars and STEREO-A. Based on observational evidence,we assume that the CME-driven shock experienced a large lateral expansion shortly after the eruption, which is modelled by a double Gaussian velocity profile in this simulation. We use the in-situ shock arrival times and the observed CME speeds at multiple spacecraft near Earth and Mars as constraints to adjust the input model parameters. The modelled time intensity profiles and fluence for energetic protons are then compared with observations. Reasonable agreements with observations at Mars and STEREO-A are found. The simulated results at Earth differ from observations of GOES-15. Instead, the simulated results at a heliocentric longitude 20 degree west to Earth fit reasonably well with the GOES observation. This can be explained if the pre-event solar wind magnetic field at Earth is not described by a nominal Parker field. Our results suggest that a large lateral expansion of the CME-driven shock and a distorted interplanetary magnetic field due to previous events can be important in understanding this GLE event.


title: Primality of numbers of the form $ap^{k}+1$

abstract: In 1876, Edouard Lucas showed that if an integer $b$ exists such that $b^{n-1} \equiv 1 (\mathrm{mod} \ n)$ and $b^{(n-1)/p} \not\equiv 1( \mathrm{mod} \ n)$ for all prime divisors $p$ of $n-1$ , then $n$ is prime, a result known as Lucas's converse of Fermat's little theorem. This result was considerably improved by Henry Pocklington in 1914 when he showed that it's not necessary to know all the prime factors of $n-1$ in order to determine if $n$ is prime. In this paper we optimize Pocklington's primality test for integers of the form $ap^{k}+1$ where $p$ is prime, $a<p$, $k\ge 1$. An extension of Lucas's converse of Fermat's little theorem is given. We also prove a new general-purpose primality test that requires that only a single odd prime divisor of $n-1$ be found for the test to be implemented. Contrary to the well-known result: There are infinitely many Fermat pseudoprimes to any base $b$; In this paper we prove the finitude of Fermat pseudoprimes in some forms of integers.


title: Designing Accurate Emulators for Scientific Processes using Calibration-Driven Deep Models

abstract: Predictive models that accurately emulate complex scientific processes can achieve exponential speed-ups over numerical simulators or experiments, and at the same time provide surrogates for improving the subsequent analysis. Consequently, there is a recent surge in utilizing modern machine learning (ML) methods, such as deep neural networks, to build data-driven emulators. While the majority of existing efforts has focused on tailoring off-the-shelf ML solutions to better suit the scientific problem at hand, we study an often overlooked, yet important, problem of choosing loss functions to measure the discrepancy between observed data and the predictions from a model. Due to lack of better priors on the expected residual structure, in practice, simple choices such as the mean squared error and the mean absolute error are made. However, the inherent symmetric noise assumption made by these loss functions makes them inappropriate in cases where the data is heterogeneous or when the noise distribution is asymmetric. We propose Learn-by-Calibrating (LbC), a novel deep learning approach based on interval calibration for designing emulators in scientific applications, that are effective even with heterogeneous data and are robust to outliers. Using a large suite of use-cases, we show that LbC provides significant improvements in generalization error over widely-adopted loss function choices, achieves high-quality emulators even in small data regimes and more importantly, recovers the inherent noise structure without any explicit priors.


title: Many visits TSP revisited

abstract: We study the Many Visits TSP problem, where given a number $k(v)$ for each of $n$ cities and pairwise (possibly asymmetric) integer distances, one has to find an optimal tour that visits each city $v$ exactly $k(v)$ times. The currently fastest algorithm is due to Berger, Kozma, Mnich and Vincze [SODA 2019, TALG 2020] and runs in time and space $\mathcal{O}^*(5^n)$. They also show a polynomial space algorithm running in time $\mathcal{O}^*(16^{n+o(n)})$. In this work, we show three main results: (i) A randomized polynomial space algorithm in time $\mathcal{O}^*(2^nD)$, where $D$ is the maximum distance between two cities. By using standard methods, this results in $(1+\epsilon)$-approximation in time $\mathcal{O}^*(2^n\epsilon^{-1})$. Improving the constant $2$ in these results would be a major breakthrough, as it would result in improving the $\mathcal{O}^*(2^n)$-time algorithm for Directed Hamiltonian Cycle, which is a 50 years old open problem. (ii) A tight analysis of Berger et al.'s exponential space algorithm, resulting in $\mathcal{O}^*(4^n)$ running time bound. (iii) A new polynomial space algorithm, running in time $\mathcal{O}(7.88^n)$.


title: Secure Single-Server Nearly-Identical Image Deduplication

abstract: Cloud computing is often utilized for file storage. Clients of cloud storage services want to ensure the privacy of their data, and both clients and servers want to use as little storage as possible. Cross-user deduplication is one method to reduce the amount of storage a server uses. Deduplication and privacy are naturally conflicting goals, especially for nearly-identical (``fuzzy'') deduplication, as some information about the data must be used to perform deduplication. Prior solutions thus utilize multiple servers, or only function for exact deduplication. In this paper, we present a single-server protocol for cross-user nearly-identical deduplication based on secure locality-sensitive hashing (SLSH). We formally define our ideal security, and rigorously prove our protocol secure against fully malicious, colluding adversaries with a proof by simulation. We show experimentally that the individual parts of the protocol are computationally feasible, and further discuss practical issues of security and efficiency.


title: Revealing the modal content of obstructed beams

abstract: In this work, we propose a predictor/indicator of the self-healing ability of coherent structured light beams: the field's modal content. Specifically, the fidelity between the obstructed and unobstructed beams' modal spectrum serves as a useful measure of the degree to which the beam will likely self-heal after an arbitrary obstacle. Since any optical field can be decomposed in terms of any chosen orthonormal basis, this analysis is, therefore, less restrictive than other methods for determining self-healing ability. Furthermore, since modal content is propagation invariant, this allows beam self-reconstruction to be studied in this way at any convenient transverse plane. As a case study, we present convincing experimental evidence for the superiority of the self-healing properties of Laguerre-Gaussian over Bessel-Gaussian beams; analysis that is facilitated primarily by the proposed measure.