title: Parallel and Distributed Algorithms for the housing allocation Problem

abstract: We give parallel and distributed algorithms for the housing allocation problem. In this problem, there is a set of agents and a set of houses. Each agent has a strict preference list for a subset of houses. We need to find a matching such that some criterion is optimized. One such criterion is Pareto Optimality. A matching is Pareto optimal if no coalition of agents can be strictly better off by exchanging houses among themselves. We also study the housing market problem, a variant of the housing allocation problem, where each agent initially owns a house. In addition to Pareto optimality, we are also interested in finding the core of a housing market. A matching is in the core if there is no coalition of agents that can be better off by breaking away from other agents and switching houses only among themselves. In the first part of this work, we show that computing a Pareto optimal matching of a house allocation is in {\bf CC} and computing the core of a housing market is {\bf CC}-hard. Given a matching, we also show that verifying whether it is in the core can be done in {\bf NC}. We then give an algorithm to show that computing a maximum Pareto optimal matching for the housing allocation problem is in {\bf RNC}^2 and quasi-{\bf NC}^2. In the second part of this work, we present a distributed version of the top trading cycle algorithm for finding the core of a housing market. To that end, we first present two algorithms for finding all the disjoint cycles in a functional graph: a Las Vegas algorithm which terminates in $O(\log l)$ rounds with high probability, where $l$ is the length of the longest cycle, and a deterministic algorithm which terminates in $O(\log^* n \log l)$ rounds, where $n$ is the number of nodes in the graph. Both algorithms work in the synchronous distributed model and use messages of size $O(\log n)$.


title: The prime index function

abstract: In this paper we introduce the prime index function \begin{align}\iota(n)=(-1)^{\pi(n)},\nonumber \end{align} where $\pi(n)$ is the prime counting function. We study some elementary properties and theories associated with the partial sums of this function given by\begin{align}\xi(x):=\sum \limits_{n\leq x}\iota(n).\nonumber \end{align}We show that a prime $p>2$ is a twin prime if and only if $\xi(p)=\xi(p+2)$. We also relate the prime index function to Cramer's conjecture by showing that \begin{align}|\xi(p_{n+1})-\xi(p_n)|+2=p_{n+1}-p_n.\nonumber \end{align}That is, Cramer's conjecture can be stated as \begin{align}\xi(p_{n+1})-\xi(p_n)\ll (\log p_n)^2.\nonumber \end{align}This reduces the problem to obtaining very good estimates of the second prime index function.


title: Locality-Sensitive Sketching for Resilient Network Flow Monitoring

abstract: Network monitoring is vital in modern clouds and data center networks for traffic engineering, network diagnosis, network intrusion detection, which need diverse traffic statistics ranging from flow size distributions to heavy hitters. To cope with increasing network rates and massive traffic volumes, sketch based approximate measurement has been extensively studied to trade the accuracy for memory and computation cost, which unfortunately, is sensitive to hash collisions. In addition, deploying the sketch involves fine-grained performance control and instrumentation. This paper presents a locality-sensitive sketch (LSS) to be resilient to hash collisions. LSS proactively minimizes the estimation error due to hash collisions with an autoencoder based optimization model, and reduces the estimation variance by keeping similar network flows to the same bucket array. To illustrate the feasibility of the sketch, we develop a disaggregated monitoring application that supports non-intrusive sketching deployment and native network-wide analysis. Testbed shows that the framework adapts to line rates and provides accurate query results. Real-world trace-driven simulations show that LSS remains stable performance under wide ranges of parameters and dramatically outperforms state-of-the-art sketching structures, with over $10^3$ to $10^5$ times reduction in relative errors for per-flow queries as the ratio of the number of buckets to the number of network flows reduces from 10\% to 0.1\%.


title: Transportless equilibration in isolated many-body quantum systems

abstract: A general analytical theory of temporal relaxation processes in isolated quantum systems with many degrees of freedom is elaborated, which unifies and substantially amends several previous approximations. Specifically, the Fourier transform of the initial energy distribution is found to play a key role, which is furthermore equivalent to the so-called survival probability in case of a pure initial state. The main prerequisite is the absence of any notable transport currents, caused for instance by some initially unbalanced local densities of particles, energy, and so on. In particular, such a transportless relaxation scenario naturally arises when both the system Hamiltonian and the initial non-equilibrium state do not exhibit any spatial inhomogeneities on macroscopic scales. A further requirement is that the relaxation must not be notably influenced by any approximate (but not exact) constant of motion or metastable state. The theoretical predictions are compared with various experimental and numerical results from the literature.


title: Nonlinear wave interactions as a model for naturally enhanced ion-acoustic lines in the ionosphere

abstract: Incoherent radar scatter from the ionosphere will, for equilibrium conditions, show two symmetric ion-acoustic lines, one for each direction of wave propagation. Many observation, from the EISCAT Svalbard Radar (ESR) for instance, demonstrate that the symmetry of this ion line can be broken, accompanied by an enhanced, nonthermal, level of fluctuations, i.e., Naturally Enhanced Ion-Acoustic Lines (NEIALs). Several models have been proposed for explaining these naturally enhanced lines. Here, we consider one of these, suggesting that decay of electron beam excited Langmuir waves gives rise to enhanced asymmetric ion lines in the ionosphere. We use a weak-turbulence approximation, and identify crucial parameters for Langmuir decay processes to be effective in generating the observed signals.


title: ABJM quantum spectral curve at twist 1: algorithmic perturbative solution

abstract: We present an algorithmic perturbative solution of ABJM quantum spectral curve at twist 1 in sl(2) sector for arbitrary spin values, which can be applied to, in principle, arbitrary order of perturbation theory. We determined the class of functions -- products of rational functions in spectral parameter with sums of Baxter polynomials and Hurwitz functions -- closed under elementary operations, such as shifts and partial fractions, as well as differentiation. It turns out, that this class of functions is also sufficient for finding solutions of inhomogeneous Baxter equations involved. For the latter purpose we present recursive construction of the dictionary for the solutions of Baxter equations for given inhomogeneous parts. As an application of the proposed method we present the computation of anomalous dimensions of twist 1 operators at six loop order. There is still a room for improvements of the proposed algorithm related to the simplifications of arising sums. The advanced techniques for their reduction to the basis of generalized harmonic sums will be the subject of subsequent paper. We expect this method to be generalizable to higher twists as well as to other theories, such as N=4 SYM.


title: Cycles and Patterns in the Sieve of Eratosthenes

abstract: We describe recurring patterns of numbers that survive each wave of the Sieve of Eratosthenes, including symmetries, uniform subdivisions, and quantifiable, predictive cycles that characterize their distribution across the number line. We generalize these results to numbers that are relatively prime to arbitrary sets of prime numbers and derive additional insights about the distribution of integers counted by Euler's phi-function.


title: Waveform systematics for binary neutron star gravitational wave signals: Effects of spin, precession, and the observation of electromagnetic counterparts

abstract: Extracting the properties of a binary system emitting gravitational waves relies on models describing the last stages of the compact binary coalescence. In this article, we study potential biases inherent to current tidal waveform approximants for spinning and precessing systems. We perform a Bayesian study to estimate intrinsic parameters of highly spinning binary neutron star systems. Our analysis shows that one has to include the quadrupolar deformation of the neutron stars due to their rotation once dimensionless spins above $\chi \sim 0.20$ are reached, otherwise the extracted intrinsic parameters are systematically biased. We find that at design sensitivity of Advanced LIGO and Virgo, it seems unlikely that for GW170817-like sources a clear imprint of precession will be visible in the analysis of the signal employing current waveform models. However, precession effects might be detectable for unequal mass configurations with spins larger than $\chi>0.2$. We finalize our study by investigating possible benefits of a combined gravitational wave and electromagnetic detection. The presence of electromagnetic counterparts help in reducing the dimensionality of the parameter space with constraints on the sky location, source distance, and inclination. However, we note that although a small improvement in the estimation of the tidal deformability parameter is seen in these cases, changes in the intrinsic parameters are overall very small.


title: Development and simulations of Enhanced Lateral Drift Sensors

abstract: We present the concept of a new type of silicon tracking sensor called Enhanced Lateral Drift (ELAD) sensor. In ELAD sensors the spatial resolution of the impact position of ionising particles is improved by a dedicated charge sharing mechanism, which is achieved by a non-homogeneous electric field in the lateral direction in the sensor bulk. The non-homogeneous electric field is created by buried doping implants with a higher concentration with respect to the background concentration of the bulk. The resulting position-dependent charge sharing allows for an improved interpolation of the impact position. TCAD-based electric field simulations for 2D and 3D geometries as well as transient simulations with a traversing particle for the 2D geometry have been carried out. The electric field profiles have further been optimised for position resolution. The simulations show a strong dependence of the charge sharing mechanism on the buried implant concentration. Optimal values for the buried implant concentration allow for nearly linear charge sharing between two readout electrodes as a function of the impact position. Additionally, the foreseen production technique combining silicon epitaxy and ion beam implantation is outlined.


title: Feed-down effect on $\Lambda$ spin polarization

abstract: We develop a theoretical framework to study the feed-down effect of higher-lying strange baryons on the spin polarization of the $\Lambda$ hyperon. In this framework, we consider two-body decays through strong, electromagnetic, and weak processes and derive general formulas for the angular distribution and spin polarization of the daughter particle by adopting the helicity formalism. Using the realistic experimental data as input, we explore the feed-down contribution to the global and the local $\Lambda$ polarizations and find that such a contribution suppresses the primordial $\Lambda$ polarization, which is not strong enough to resolve the discrepancy between the current theoretical and the experimental results on the azimuthal-angle dependence of $\Lambda$ polarization. Our paper may also be useful for the measurement of spin polarization of baryons heavier than $\Lambda$ (e.g., $\Xi^-$) in future experiments.


title: A First Course in Data Science

abstract: Data science is a discipline that provides principles, methodology and guidelines for the analysis of data for tools, values, or insights. Driven by a huge workforce demand, many academic institutions have started to offer degrees in data science, with many at the graduate, and a few at the undergraduate level. Curricula may differ at different institutions, because of varying levels of faculty expertise, and different disciplines (such as Math, computer science, and business etc) in developing the curriculum. The University of Massachusetts Dartmouth started offering degree programs in data science from Fall 2015, at both the undergraduate and the graduate level. Quite a few articles have been published that deal with graduate data science courses, much less so dealing with undergraduate ones. Our discussion will focus on undergraduate course structure and function, and specifically, a first course in data science. Our design of this course centers around a concept called the data science life cycle. That is, we view tasks or steps in the practice of data science as forming a process, consisting of states that indicate how it comes into life, how different tasks in data science depend on or interact with others until the birth of a data product or the reach of a conclusion. Naturally, different pieces of the data science life cycle then form individual parts of the course. Details of each piece are filled up by concepts, techniques, or skills that are popular in industry. Consequently, the design of our course is both "principled" and practical. A significant feature of our course philosophy is that, in line with activity theory, the course is based on the use of tools to transform real data in order to answer strongly motivated questions related to the data.


title: On a stochastic version of Lanchester's model of combat

abstract: Lanchester's model of combat has certain deficiencies in its standard form arising from the neglect of the influence of random fluctuations. Several approaches to rectify this have been proposed and various results are scattered throughout the literature. Here, a discrete-time stochastic version, which is amenable to exact solution, is revisited with the aim of deriving key results within one setting. The exposition simplifies and provides refinements to earlier derivations and analysis.


title: Polarization transfer in hyperon decays and its effect in relativistic nuclear collisions

abstract: We calculate the contribution to the polarization of $\Lambda$ hyperons in relativistic nuclear collisions at high energy from the decays of $\Sigma^*(1385)$ and $\Sigma^0$, which are the predominant sources of $\Lambda$ production besides the primary component, as a function of the $\Lambda$ momentum. Particularly, we estimate the longitudinal component of the mean spin vector as a function of the azimuthal angle in the transverse plane, assuming that primary $\Sigma^*$ and $\Sigma^0$ polarization follow the predictions of local thermodynamic equilibrium in a relativistic fluid. Provided that the rapidity dependence around midrapidity of polarization is negligible, we find that this component of the overall spin vector has a very similar pattern to the primary one. Therefore, we conclude that the secondary decays cannot account for the discrepancy in sign between experimental data and hydrodynamic model predictions of the longitudinal polarization of $\Lambda$ hyperons recently measured by the STAR experiment at RHIC.


title: Key-agreement based on automaton groups

abstract: We suggest several automaton groups as key-agreement platforms for Anshl-Anshel-Goldfeld metascheme, they include Grigorchuk and universal Grigorchuk groups, Hanoi 3-Towers group, Basilica group and a subgroup of the affine group with the unsolvable conjugacy problem


title: Batch-Size Independent Regret Bounds for the Combinatorial Multi-Armed Bandit Problem

abstract: We consider the combinatorial multi-armed bandit (CMAB) problem, where the reward function is nonlinear. In this setting, the agent chooses a batch of arms on each round and receives feedback from each arm of the batch. The reward that the agent aims to maximize is a function of the selected arms and their expectations. In many applications, the reward function is highly nonlinear, and the performance of existing algorithms relies on a global Lipschitz constant to encapsulate the function's nonlinearity. This may lead to loose regret bounds, since by itself, a large gradient does not necessarily cause a large regret, but only in regions where the uncertainty in the reward's parameters is high. To overcome this problem, we introduce a new smoothness criterion, which we term \emph{Gini-weighted smoothness}, that takes into account both the nonlinearity of the reward and concentration properties of the arms. We show that a linear dependence of the regret in the batch size in existing algorithms can be replaced by this smoothness parameter. This, in turn, leads to much tighter regret bounds when the smoothness parameter is batch-size independent. For example, in the probabilistic maximum coverage (PMC) problem, that has many applications, including influence maximization, diverse recommendations and more, we achieve dramatic improvements in the upper bounds. We also prove matching lower bounds for the PMC problem and show that our algorithm is tight, up to a logarithmic factor in the problem's parameters.


title: Dark Matter freezeout in modified cosmological scenarios

abstract: We study the effects of modifying the expansions history of the Universe on Dark Matter freezeout. We derived a modified Boltzmann equation for freeze-out for an arbitrary energy density in the early Universe and provide an analytic approach using some approximations. We then look at the required thermally averaged cross sections needed to obtain the correct relic density for the specific case where the energy density consists of radiation plus one extra component which cools faster. We compare our analytic approximation to a numerical solutions. We find that it gives reasonable results for most of the parameter space explored, being at most a factor of order one away from the measured value. We find that if the new contribution to the energy density is comparable to the radiation density, then a much smaller cross section for Dark Matter annihilation is required. This would lead to weak scale Dark Matter being much more difficult to detect and opens up the possibility that much heavier Dark Matter could undergo freezeout without violating perturbative unitarity.


title: A relativistic electron in an anisotropic conduction band

abstract: The Dirac equation is extended for a relativistic electron in an orthorhombically-anisotropic conduction band. Its covariance is established with general proper and improper Lorentz transformations. In the non-relativistic limit, the kinetic and Zeeman energy terms of the Hamiltonian are both determined by the same three effective masses, and the quantum spin Hall effect is derived. This has important consequences for magnetic measurements of many classes of clean anisotropic semiconductors, metals, and superconductors. The Zeeman energy is vanishingly small for magnetic fields parallel to clean monolayers and in all directions in quasi-one-dimensional materials.


title: Diffusive transport in a quasiperiodic Fibonacci chain: absence of many-body localization at small interactions

abstract: We study high-temperature magnetization transport in a many-body spin-1/2 chain with on-site quasiperiodic potential governed by the Fibonacci rule. In the absence of interactions it is known that the system is critical with the transport described by a continuously varying dynamical exponent (from ballistic to localized) as a function of the on-site potential strength. Upon introducing weak interactions, we find that an anomalous noninteracting dynamical exponent becomes diffusive for any potential strength. This is borne out by a boundary-driven Lindblad dynamics as well as unitary dynamics, with agreeing diffusion constants. This must be contrasted to random potential where transport is subdiffusive at such small interactions. Mean-field treatment of the dynamics for small U always slows down the non-interacting dynamics to subdiffusion, and is therefore unable to describe diffusion in an interacting quasiperiodic system. Finally, briefly exploring larger interactions we find a regime of interaction-induced subdiffusive dynamics, despite the on-site potential itself having no "rare-regions".


title: Extremely elastic soliton crystals generated in a passively mode-locked tunable high-repetition-rate fiber laser

abstract: We present the first direct observation of the bound state of multiple dissipative optical solitons in which bond length and bond strength can be individually controlled in a broad range in a regular manner. We have observed experimentally a new type of stable and extremely elastic soliton crystals that can be stretched and compressed many times conserving their structure by adjusting the bond properties in real time in a specially designed passively mode-locked fiber laser incorporating highly asymmetric tunable Mach-Zehnder interferometer. The temporal structure and dynamics of the generated soliton crystals have been studied using an asynchronous optical sampling system with picosecond resolution. We demonstrated that stable and robust soliton crystal can be formed by two types of primitive structures: single dissipative solitons, and(or) pairs of dissipative soliton and pulse with lower amplitude. Continuous stretching and compression of a soliton crystal with extraordinary high ratio of more than 30 has been demonstrated with a smallest recorded separation between pulses as low as 5 ps corresponding to an effective repetition frequency of 200 GHz. Collective pulse dynamics, including soliton crystal self-assembling, cracking and transformation of crystals comprising pulse pairs to the crystals of similar pulses has been observed experimentally.


title: Cooperative Distributed Robust Control of Modular Mobile Robots with Bounded Curvature and Velocity

abstract: A novel motion control system for Compliant Framed wheeled Modular Mobile Robots (CFMMR) is studied in this paper. This type of wheeled mobile robot uses rigid axles coupled by compliant frame modules to provide both full suspension and enhanced steering capability without additional hardware. The proposed control system is developed by combining a bounded curvature-based kinematic controller and a nonlinear damping dynamic controller. In particular, multiple forms of controller interaction are examined. A twoaxle scout CFMMR configuration is used to evaluate the different control structures. Experimental results verify efficient motion control of posture regulation.


title: Vision-based Unscented FastSLAM for Mobile Robot

abstract: This paper presents a vision-based Unscented FastSLAM (UFastSLAM) algorithm combing the Rao-Blackwellized particle filter and Unscented Kalman filte(UKF). The landmarks are detected by a binocular vision to integrate localization and mapping. Since such binocular vision system generally inherits larger measurement errors, it is suitable to adopt Unscented FastSLAM to improve the performance of localization and mapping. Unscented FastSLAM takes advantage of UKF instead of the linear approximations of the nonlinear function where the effective number of particles is used as the criteria to reduce the particle degeneration. Simulations and experiments are carried out to demonstrate that the Unscented FastSLAM algorithm can achieve much better performance in the vision-based system than FastSLAM2.0 algorithm on the accuracy and robustness.


title: Terrain Inclination Aided Three Dimensional Localization and Mapping for an Outdoor Mobile Robot

abstract: A new 3D localization and mapping techinque with terrain inclination assistance is proposed in this paper to allow a robot to identify its location and build a global map in an outdoor environment. The Iterative Closest Points (ICP) algorithm and terrain inclination-based localization are combined together to achieve accurate and fast localization and mapping. Inclinations of the terrains the robot navigates are used to achieve local localization during the interval between two laser scans. Using the results of the above localization as the initial condition, the ICP algorithm is then applied to align the overlapped laser scan maps to update the overhanging obstacles for building a global map of the surrounding area. Comprehensive experiments were carried out for the validation of the proposed 3D localization and mapping technique. The experimental results show that the proposed technique could reduce time consumption and improve the accuracy of the performance.


title: Driving Driven Lattice Gases to Identify Their Universality Classes

abstract: The critical behavior of driven lattice gas models has been studied for decades as a paradigm to explore nonequilibrium phase transitions and critical phenomena. However, there exists a long-standing controversy in the universality classes to which they belong. This is of paramount importance as it implies the question of whether or not a microscopic model and its mesoscopic field theory may possess different symmetries in nonequilibrium critical phenomena in contrast to their equilibrium counterparts. Here, we heat with finite rates two generic models of driven lattice gases across their respective nonequilibrium critical points into further nonequilibrium conditions. Employing the theory of finite-time scaling, we are able to unambiguously discriminate the universality classes between the two models. In particular, the infinitely driven lattice gas and the randomly driven lattice gas models belong to different universality classes. These results show that finite-time scaling is effective even in nonequilibrium phase transitions.


title: Finding cuts of bounded degree: complexity, FPT and exact algorithms, and kernelization

abstract: A matching cut is a partition of the vertex set of a graph into two sets $A$ and $B$ such that each vertex has at most one neighbor in the other side of the cut. The MATCHING CUT problem asks whether a graph has a matching cut, and has been intensively studied in the literature. Motivated by a question posed by Komusiewicz et al. [IPEC 2018], we introduce a natural generalization of this problem, which we call $d$-CUT: for a positive integer $d$, a $d$-cut is a bipartition of the vertex set of a graph into two sets $A$ and $B$ such that each vertex has at most $d$ neighbors across the cut. We generalize (and in some cases, improve) a number of results for the MATCHING CUT problem. Namely, we begin with an NP-hardness reduction for $d$-CUT on $(2d+2)$-regular graphs and a polynomial algorithm for graphs of maximum degree at most $d+2$. The degree bound in the hardness result is unlikely to be improved, as it would disprove a long-standing conjecture in the context of internal partitions. We then give FPT algorithms for several parameters: the maximum number of edges crossing the cut, treewidth, distance to cluster, and distance to co-cluster. In particular, the treewidth algorithm improves upon the running time of the best known algorithm for MATCHING CUT. Our main technical contribution, building on the techniques of Komusiewicz et al. [IPEC 2018], is a polynomial kernel for $d$-CUT for every positive integer $d$, parameterized by the distance to a cluster graph. We also rule out the existence of polynomial kernels when parameterizing simultaneously by the number of edges crossing the cut, the treewidth, and the maximum degree. Finally, we provide an exact exponential algorithm slightly faster than the naive brute force approach running in time $O^*(2^n)$.


title: Optimal Statistical Rates for Decentralised Non-Parametric Regression with Linear Speed-Up

abstract: We analyse the learning performance of Distributed Gradient Descent in the context of multi-agent decentralised non-parametric regression with the square loss function when i.i.d. samples are assigned to agents. We show that if agents hold sufficiently many samples with respect to the network size, then Distributed Gradient Descent achieves optimal statistical rates with a number of iterations that scales, up to a threshold, with the inverse of the spectral gap of the gossip matrix divided by the number of samples owned by each agent raised to a problem-dependent power. The presence of the threshold comes from statistics. It encodes the existence of a "big data" regime where the number of required iterations does not depend on the network topology. In this regime, Distributed Gradient Descent achieves optimal statistical rates with the same order of iterations as gradient descent run with all the samples in the network. Provided the communication delay is sufficiently small, the distributed protocol yields a linear speed-up in runtime compared to the single-machine protocol. This is in contrast to decentralised optimisation algorithms that do not exploit statistics and only yield a linear speed-up in graphs where the spectral gap is bounded away from zero. Our results exploit the statistical concentration of quantities held by agents and shed new light on the interplay between statistics and communication in decentralised methods. Bounds are given in the standard non-parametric setting with source/capacity assumptions.


title: Performance Engineering for Real and Complex Tall & Skinny Matrix Multiplication Kernels on GPUs

abstract: General matrix-matrix multiplications with double-precision real and complex entries (DGEMM and ZGEMM) in vendor-supplied BLAS libraries are best optimized for square matrices but often show bad performance for tall & skinny matrices, which are much taller than wide. NVIDIA's current CUBLAS implementation delivers only a fraction of the potential performance as indicated by the roofline model in this case. We describe the challenges and key characteristics of an implementation that can achieve close to optimal performance. We further evaluate different strategies of parallelization and thread distribution, and devise a flexible, configurable mapping scheme. To ensure flexibility and allow for highly tailored implementations we use code generation combined with autotuning. For a large range of matrix sizes in the domain of interest we achieve at least 2/3 of the roofline performance and often substantially outperform state-of-the art CUBLAS results on an NVIDIA Volta GPGPU.


title: The 'Delayed Choice Quantum Eraser' Neither Erases Nor Delays

abstract: It is demonstrated that 'quantum eraser' (QE) experiments do not erase any information. Nor do they demonstrate retrocausation or 'temporal nonlocality' in their 'delayed choice' form, beyond standard EPR correlations. It is shown that the erroneous erasure claims arise from assuming that the improper mixed state of the signal photon physically prefers either the 'which way' or 'both ways' basis, when no such preference is warranted. The latter point is illustrated through comparison of the QE spatial state space with the spin-1/2 space of particles in the EPR-spin experiment.


title: Somatic mutations render human exome and pathogen DNA more similar

abstract: Immunotherapy has recently shown important clinical successes in a substantial number of oncology indications. Additionally, the tumor somatic mutation load has been shown to associate with response to these therapeutic agents, and specific mutational signatures are hypothesized to improve this association, including signatures related to pathogen insults. We sought to study in silico the validity of these observations and how they relate to each other. We first addressed whether somatic mutations typically involved in cancer may increase, in a statistically meaningful manner, the similarity between common pathogens and the human exome. Our study shows that common mutagenic processes increase, in the upper range of biologically plausible frequencies, the similarity between cancer exomes and pathogen DNA at a scale of 12-16 nucleotide sequences and established that this increased similarity is due to the specific mutation distribution of the considered mutagenic processes. Next, we studied the impact of mutation rate and showed that increasing mutation rate generally results in an increased similarity between the cancer exome and pathogen DNA, at a scale of 4-5 amino acids. Finally, we investigated whether the considered mutational processes result in amino-acid changes with functional relevance that are more likely to be immunogenic. We showed that functional tolerance to mutagenic processes across species generally suggests more resilience to mutagenic processes that are due to exposure to elements of nature than to mutagenic processes that are due to exposure to cancer-causing artificial substances. These results support the idea that recognition of pathogen sequences as well as differential functional tolerance to mutagenic processes may play an important role in the immune recognition process involved in tumor infiltration by lymphocytes.

categories: q-bio.GN cs.IT math.IT math.PR stat.AP stat.CO


title: J-PLUS: Discovery and characterisation of ultracool dwarfs using Virtual Observatory tools

abstract: Ultracool dwarfs (UCDs) are objects with spectral types equal or later than M7. Most of them have been discovered using wide-field imaging surveys. The Virtual Observatory (VO) has proven to be of great utility to efficiently exploit these astronomical resources. We aim to validate a VO methodology designed to discover and characterize UCDs in the J-PLUS photometric survey. J-PLUS is a multiband survey carried out with the wide angle T80Cam optical camera mounted on the 0.83-m telescope JAST/T80 in the Observatorio Astrof\'isico de Javalambre. In this work we make use of the Internal Data Release (IDR) covering 528 deg$^2$. We complement J-PLUS photometry with other catalogues in the optical and IR using VOSA, a VO tool that estimates physical parameters from the spectral energy distribution fitting to collections of theoretical models. Objects identified as UCDs are distinguished from background M giants and highly reddened stars using parallaxes and proper motions from Gaia DR2. We identify 559 UCDs, ranging from i=16.2 to 22.4 mag, of which 187 are candidate UCDs not previously reported in the literature. This represents an increase in the number of known UCDs of about 50% in the studied region of the sky, particularly at the faint end of our sensitivity, which is interesting as reference for future wide and deep surveys such as Euclid. Three candidates constitute interesting targets for exoplanet surveys because of their proximity (<40 pc). We also analyze the kinematics of UCDs in our catalogue and find evidence that it is consistent with a Galactic thin-disk population, except for 6 objects that might be members of the thick disk. The results obtained validate the proposed methodology, which will be used in future J-PLUS and J-PAS releases. Considering the region of the sky covered by the IDR used, we foresee to discover 3,000-3,500 new UCDs at the end of the J-PLUS project.


title: K\"ahler packings of projective complex manifolds

abstract: In this note we show that the multipoint Seshadri constant determines the maximum possible radii of embeddings of K\"ahler balls and vice versa.


title: Interpolation by Linear Functions on an $n$-Dimensional Ball

abstract: By $B=B(x^{(0)};R)$ we denote the Euclidean ball in ${\mathbb R}^n$ given by the inequality $\|x-x^{(0)}\|\leq R$. Here $x^{(0)}\in{\mathbb R}^n, R>0$, $\|x\|:=\left(\sum_{i=1}^n x_i^2\right)^{1/2}$. We mean by $C(B)$ the space of continuous functions $f:B\to{\mathbb R}$ with the norm $\|f\|_{C(B)}:=\max_{x\in B}|f(x)|$ and by $\Pi_1\left({\mathbb R}^n\right)$ the set of polynomials in $n$ variables of degree $\leq 1$, i.e., linear functions on ${\mathbb R}^n$. Let $x^{(1)}, \ldots, x^{(n+1)}$ be the vertices of $n$-dimensional nondegenerate simplex $S\subset B$. The interpolation projector $P:C(B)\to \Pi_1({\mathbb R}^n)$ corresponding to $S$ is defined by the equalities $Pf\left(x^{(j)}\right)=f\left(x^{(j)}\right).$ We obtain the formula to compute the norm of $P$ as an operator from $C(B)$ into $C(B)$ via $x^{(0)}$, $R$ and coefficients of basic Lagrange polynomials of $S$. In more details we study the case when $S$ is a regular simplex inscribed into $B_n=B(0,1)$.


title: The kinetic exclusion process: a tale of two fields

abstract: We introduce a general class of stochastic lattice gas models, and derive their fluctuating hydrodynamics description in the large size limit under a local equilibrium hypothesis. The model consists in energetic particles on a lattice subject to exclusion interactions, which move and collide stochastically with energy-dependent rates. The resulting fluctuating hydrodynamics equations exhibit nonlinear coupled particle and energy transport, including particle currents due to temperature gradients (Soret effect) and energy flow due to concentration gradients (Dufour effect). The microscopic dynamical complexity is condensed in just two matrices of transport coefficients: the diffusivity matrix (or equivalently the Onsager matrix) generalizing Fick-Fourier's law, and the mobility matrix controlling current fluctuations. Both transport coefficients are coupled via a fluctuation-dissipation theorem, suggesting that the noise terms affecting the local currents have Gaussian properties. We further prove the positivity of entropy production in terms of the microscopic dynamics. The so-called kinetic exclusion process has as limiting cases two of the most paradigmatic models of nonequilibrium physics, namely the symmetric simple exclusion process of particle diffusion and the Kipnis-Marchioro-Presutti model of heat flow, making it the ideal testbed where to further develop modern theories of nonequilibrium behavior.

categories: cond-mat.stat-mech cond-mat.soft math-ph math.MP physics.flu-dyn


title: Wide-field SU(1,1) interferometer

abstract: An ${\rm SU(1,1)}$ interferometer uses a sequence of two optical parametric amplifiers for achieving sub-shot-noise sensitivity to a phase shift introduced in between. We present the first realization of a wide-field ${\rm SU(1,1)}$ interferometer, where the use of a focusing element enables spatially multimode operation within a broad angle. Over this angle, the interference phase is found to be flat. This property is important for the high sensitivity to the phase front disturbance. Further, $-4.3\pm0.7$ dB quadrature squeezing, an essential requirement to the high sensitivity, is experimentally demonstrated for plane-wave modes inside the interferometer. Such an interferometer is useful not only for quantum metrology, but also in remote sensing, enhanced sub-shot-noise imaging, and quantum information processing.


title: Blitz-starting QUIC Connections

abstract: In this paper, we revisit the idea to remove Slow Start from congestion control. To do so, we build upon the newly gained freedom of transport protocol extendability offered by QUIC to hint bandwidth estimates from a typical web client to a server. Using this bandwidth estimate, we bootstrap congestion windows of new connections to quickly utilize available bandwidth. This custom flow initialization removes the common early exit of Slow Start and thus fuels short flow fairness with long-running connections. Our results indicate that we can drastically reduce flow completion time accepting some losses and thereby an inflated transmission volume. For example, for a typical DSL client, loading a 2 MB YouTube video chunk is accelerated by nearly 2x. In the worst case, we find an inflation of the transfer volume by 12% due to losses.


title: The performance guarantee of randomized perfect voting trees

abstract: In this note we study randomized voting trees, previously introduced by Fisher, Procaccia and Samorodnitsky. They speculate that a non-trivial performance guarantee may be achievable using randomized, balanced trees whose height is carefully chosen. We explore some connections to the so-called Volterra quadratic stochastic operators, and show that uniformly random voting trees cannot provide a performance guarantee that is linear in the number of individuals.


title: Persistence of Technosignatures: A Comment on Lingam and Loeb

abstract: In a recent paper in this journal, Lingam and Loeb (2018) develop an excellent heuristic for searches for biosignatures vs. technosignatures. We consider two ways in which their approach could be extended and sharpened, with focus on durability of technosignatures. We also note an important consequence of the adopted heuristic which offers strong support to the ideas of the Dysonian SETI.


title: Plasmon-polariton from a helical state in Dirac magnet

abstract: Optical field interacting with a topologically protected one-dimensional helical state is shown to support a one-dimensional plasmon-polariton that is characterized by a non-linear dispersion. In a two-dimensional Dirac magnet these electro-optical excitations are confined to domain walls, thus, offering a possibility to manipulate quantum optical states by altering magnetic domain configurations. An exact spectral equation for such a topological plasmon-polariton is derived.


title: The asymptotic induced matching number of hypergraphs: balanced binary strings

abstract: We compute the asymptotic induced matching number of the $k$-partite $k$-uniform hypergraphs whose edges are the $k$-bit strings of Hamming weight $k/2$, for any large enough even number $k$. Our lower bound relies on the higher-order extension of the well-known Coppersmith-Winograd method from algebraic complexity theory, which was proven by Christandl, Vrana and Zuiddam. Our result is motivated by the study of the power of this method as well as of the power of the Strassen support functionals (which provide upper bounds on the asymptotic induced matching number), and the connections to questions in tensor theory, quantum information theory and theoretical computer science. Phrased in the language of tensors, as a direct consequence of our result, we determine the asymptotic subrank of any tensor with support given by the aforementioned hypergraphs. In the context of quantum information theory, our result amounts to an asymptotically optimal $k$-party stochastic local operations and classical communication (slocc) protocol for the problem of distilling GHZ-type entanglement from a subfamily of Dicke-type entanglement.


title: Updated orbital ephemeris of the ADC source X 1822-371: a stable orbital expansion over 40 years

abstract: The source X 1822-371 is an eclipsing compact binary system with a period close to 5.57 hr and an orbital period derivative $\dot{P}_{\rm orb}$ of 1.51(7)$\times 10^{-10}$ s s$^{-1}$. The very large value of $\dot{P}_{\rm orb}$ is compatible with a super-Eddington mass transfer rate from the companion star, as suggested by X-ray and optical data. The XMM-Newton observation taken in 2017 allows us to update the orbital ephemeris and verify whether the orbital period derivative has been stable over the last 40 yr. We added to the X-ray eclipse arrival times from 1977 to 2008 two new values obtained from the RXTE and XMM-Newton observations performed in 2011 and 2017, respectively. We estimated the number of orbital cycles and the delays of our eclipse arrival times spanning 40 yr using as reference time the eclipse arrival time obtained from the Rossi-XTE observation taken in 1996. Fitting the delays with a quadratic model, we found an orbital period $P_{\rm orb}=5.57062957(20)$ hr and a $\dot{P}_{\rm orb}$ value of $1.475(54) \times 10^{-10}$ s s$^{-1}$. The addition of a cubic term to the model does not significantly improve the quality of the fit. We also determined a spin-period value of $P_{\rm spin}=0.5915669(4)$ s and its first derivative $\dot{P}_{\rm spin}= -2.595(11) \times 10^{-12}$ s s$^{-1}$. The obtained results confirm the scenario of a super-Eddington mass transfer rate; we also exclude a gravitational coupling between the orbit and the change in the oblateness of the companion star triggered by the nuclear luminosity of the companion star.


title: Demonstration of Adiabatic Variational Quantum Computing with a Superconducting Quantum Coprocessor

abstract: Adiabatic quantum computing enables the preparation of many-body ground states. This is key for applications in chemistry, materials science, and beyond. Realisation poses major experimental challenges: Direct analog implementation requires complex Hamiltonian engineering, while the digitised version needs deep quantum gate circuits. To bypass these obstacles, we suggest an adiabatic variational hybrid algorithm, which employs short quantum circuits and provides a systematic quantum adiabatic optimisation of the circuit parameters. The quantum adiabatic theorem promises not only the ground state but also that the excited eigenstates can be found. We report the first experimental demonstration that many-body eigenstates can be efficiently prepared by an adiabatic variational algorithm assisted with a multi-qubit superconducting coprocessor. We track the real-time evolution of the ground and exited states of transverse-field Ising spins with a fidelity up that can reach about 99%.


title: Unrestricted Permutation forces Extrapolation: Variable Importance Requires at least One More Model, or There Is No Free Variable Importance

abstract: This paper reviews and advocates against the use of permute-and-predict (PaP) methods for interpreting black box functions. Methods such as the variable importance measures proposed for random forests, partial dependence plots, and individual conditional expectation plots remain popular because they are both model-agnostic and depend only on the pre-trained model output, making them computationally efficient and widely available in software. However, numerous studies have found that these tools can produce diagnostics that are highly misleading, particularly when there is strong dependence among features. The purpose of our work here is to (i) review this growing body of literature, (ii) provide further demonstrations of these drawbacks along with a detailed explanation as to why they occur, and (iii) advocate for alternative measures that involve additional modeling. In particular, we describe how breaking dependencies between features in hold-out data places undue emphasis on sparse regions of the feature space by forcing the original model to extrapolate to regions where there is little to no data. We explore these effects across various model setups and find support for previous claims in the literature that PaP metrics can vastly over-emphasize correlated features in both variable importance measures and partial dependence plots. As an alternative, we discuss and recommend more direct approaches that involve measuring the change in model performance after muting the effects of the features under investigation.


title: Groupoid symmetry, constraints and quantization of General Relativity

abstract: The purpose of the current paper is twofold: to provide a conceptual link between the quantization framework based on Lie integration of algebroids proposed by N.P. Landsman in the book "Mathematical Topics between Classical and Quantum Mechanics" (1998) and the dynamical formulation of the Einstein's equation, and to clarify how from the relevant groupoid the Poisson bracket between constraints of the hamiltonian formulation emerges as an algebroid bracket. To do so, we adapt the groupoid proposed by C. Blohmann, M. Fernandes and A. Weinstein in the paper "Groupoid symmetry and constraints in General Relativity" (2010) by using a different, in our view simpler and more natural, diffeological structure. As a result, we get a different algebroid associated to it, such that the algebroid bundle can be understood as a configuration space of the dynamical formulation of General Relativity, the bracket structure between its sections being of the desired form. We point out the conceptual advantages of this perspective, as well as some interesting issues that require further investigation. We also discuss some of the difficulties that need to be overcome to apply our approach to quantize the Einstein's theory in this sense. This is an extract from the Master's thesis written under the supervision of Prof. Klaas Landsman.


title: The Generalized Nonlinear Ohm's Law: How a Strong Electric Field Influences Non-ideal MHD Effects in Dusty Protoplanetary Disks

abstract: The magnetohydrodynamics (MHD) of protoplanetary disks are strongly subject to the non-ideal MHD effects arising from the low ionization fraction of the disk gas. A strong electric field induced by gas motions can heat ionized gas particles and can thereby affect the ionization balance in the disks. Our previous studies revealed that in dusty protoplanetary disks, the Ohmic conductivity decreases with increasing electric field strength until the electrical breakdown of the disk gas occurs. In this study, we extend our previous work to more general cases where both electric and magnetic fields affect the motion of plasma particles, allowing us to study the impacts of plasma heating on all non-ideal MHD effects: Ohmic, Hall, and ambipolar diffusion. We find that the upper limit on the electric current we previously derived applies even in the presence of magnetic fields. Although the Hall and ambipolar resistivities can either increase or decrease with electric field strength depending on the abundance of charged dust grains, the Ohmic resistivity always increases with electric field strength. An order-of-magnitude estimate suggests that a large-scale electric current generated by gas motions in the inner part of protoplanetary disks could exceed the upper limit. This implies that MHD motions of the inner disk, such as the motion driven by the Hall-shear instability, could either get suppressed or trigger electrical breakdown (lightning discharge). This may have important implications for gas accretion and chondrule formation in the inner part of protoplanetary disks.


title: On pure complex spectrum for truncations of random orthogonal matrices and Kac polynomials

abstract: Let $O(2n+\ell)$ be the group of orthogonal matrices of size $\left(2n+\ell\right)\times \left(2n+\ell\right)$ equipped with the probability distribution given by normalized Haar measure. We study the probability \begin{equation*} p_{2n}^{\left(\ell\right)} = \mathbb{P}\left[M_{2n} \, \mbox{has no real eigenvalues}\right], \end{equation*} where $M_{2n}$ is the $2n\times 2n$ left top minor of a $(2n+\ell)\times(2n+\ell)$ orthogonal matrix. We prove that this probability is given in terms of a determinant identity minus a weighted Hankel matrix of size $n\times n$ that depends on the truncation parameter $\ell$. For $\ell=1$ the matrix coincides with the Hilbert matrix and we prove \begin{equation*} p_{2n}^{\left(1\right)} \sim n^{-3/8}, \mbox{ when }n \to \infty. \end{equation*} We also discuss connections of the above to the persistence probability for random Kac polynomials.


title: Unidirectional Distributed Acoustic Reflection Transducers for Quantum Applications

abstract: Recent significant advances in coupling superconducting qubits to acoustic wave resonators has led to demonstrations of quantum control of surface and bulk acoustic resonant modes as well Wigner tomography of quantum states in these modes. These advances were achieved through the efficient coupling afforded by piezoelectric materials combined with GHz-frequency acoustic Fabry-Perot cavities. Quantum control of itinerant surface acoustic waves appears in reach, but is challenging due to the limitations of conventional transducers in the appropriate GHz-frequency band. In particular, GHz-frequency unidirectional transducers would provide an important addition to the desired quantum toolbox, promising unit efficiency with directional control over the surface acoustic wave emission pattern. Here we report the design, fabrication and experimental characterization of unidirectional distributed acoustic reflection transducers (DARTs) demonstrating a high transduction frequency of 4.8 GHz with a peak directivity larger than 25 dB and a directivity greater than 15 dB over a bandwidth of 17 MHz. A numerical model reproduces the main features of the transducer response quite well, with ten adjustable parameters (most of which are constrained by geometric and physical considerations). This represents a significant step towards quantum control of itinerant quantum acoustic waves.


title: Evaluating Cascading Impact of Attacks on Resilience of Industrial Control Systems: A Design-Centric Modeling Approach

abstract: A design-centric modeling approach was proposed to model the behaviour of the physical processes controlled by Industrial Control Systems (ICS) and study the cascading impact of data-oriented attacks. A threat model was used as input to guide the construction of the CPS model where control components which are within the adversary's intent and capabilities are extracted. The relevant control components are subsequently modeled together with their control dependencies and operational design specifications. The approach was demonstrated and validated on a water treatment testbed. Attacks were simulated on the testbed model where its resilience to attacks was evaluated using proposed metrics such as Impact Ratio and Time-to-Critical-State. From the analysis of the attacks, design strengths and weaknesses were identified and design improvements were recommended to increase the testbed's resilience to attacks.


title: Hypercyclic Algebras for convolution operators of unimodular constant term

abstract: We study the existence of hypercyclic algebras for convolution operators $\Phi(D)$ on the space of entire functions whose symbol $\Phi$ has unimodular constant term. In particular, we provide new eigenvalue criteria for the existence of densely strongly algebrable sets of hypercyclic vectors.


title: One-Dimensional Transport of Ultracold Bosons

abstract: Advances in cooling and trapping of atoms have enabled unprecedented experimental control of many-body quantum systems. This led to the observation of numerous quantum phenomena, important for fundamental science, indispensable for high-precision simulations of condensed-matter systems and promising for technological applications. However, transport measurements in neutral quantum gases are still in their infancy in contrast to the central role they play in electronics. In these lectures, after reviewing nascent experiments on quantum fermionic transport, I will focus on our theoretical prediction sand the possibility of experimental observations of qualitatively new phenomena in transport of ultracold bosons which do not have a direct counterpart in quantum electronic transport in condensed matter systems. The description of this transport is based on the Luttinger liquid (LL) theory. So in the first part of the lectures I will introduce main concepts of the LL based on the functional bosonisation approach.


title: A diagrammatic presentation and its characterization of non-split compact surfaces in the 3-sphere

abstract: We give a presentation for a non-split compact surface embedded in the 3-sphere by using diagrams of spatial trivalent graphs equipped with signs and we define Reidemeister moves for such signed diagrams. We show that two diagrams of embedded surfaces are related by Reidemeister moves if and only if the surfaces represented by the diagrams are ambient isotopic in the 3-sphere.


title: Decentralized Massive MIMO Processing Exploring Daisy-chain Architecture and Recursive Algorithms

abstract: Algorithms for Massive MIMO uplink detection and downlink precoding typically rely on a centralized approach, by which baseband data from all antenna modules are routed to a central node in order to be processed. In the case of Massive MIMO, where hundreds or thousands of antennas are expected in the base-station, said routing becomes a bottleneck since interconnection throughput is limited. This paper presents a fully decentralized architecture and an algorithm for Massive MIMO uplink detection and downlink precoding based on the Stochastic Gradient Descent (SGD) method, which does not require a central node for these tasks. Through a recursive approach and very low complexity operations, the proposed algorithm provides a good trade-off between performance, interconnection throughput and latency. Further, our proposed solution achieves significantly lower interconnection data-rate than other architectures, enabling future scalability.


title: Opening the Treasure Chest in Carina

abstract: We have mapped the G287.84-0.82 cometary globule (with the Treasure Chest cluster embedded in it) in the South Pillars region of Carina (i) in [CII], 63micron [OI], and CO(11-10) using upGREAT on SOFIA and (ii) in J=2-1 transitions of CO, 13CO, C18O and J=3-2 transitions of H2CO using the APEX telescope in Chile. We probe the morphology, kinematics, and physical conditions of the molecular gas and the photon dominated regions (PDRs) in G287.84-0.82. The [CII] and [OI] emission suggest that the overall structure of the pillar (with red-shifted photo evaporating tails) is consistent with the effect of FUV radiation and winds from eta-Car and O stars in Trumpler 16. The gas in the head of the pillar is strongly influenced by the embedded cluster, whose brightest member is an O9.5V star, CPD-59 2661. The emission of the [CII] and [OI] lines peak at a position close to the embedded star, while all other tracers peak at another position lying to the north-east consistent with gas being compressed by the expanding PDR created by the embedded cluster. The molecular gas inside the globule is probed with the J=2-1 transitions of CO and isotopologues as well as H2CO, and analyzed using a non-LTE model (escape-probability approach), while we use PDR models to derive the physical conditions of the PDR. We identify at least two PDR gas components; the diffuse part (~10^4 cm^-3) is traced by [CII], while the dense (n~ 2-8x10^5 cm^-3) part is traced by [CII], [OI], CO(11-10). Using the F=2-1 transition of [13CII] detected at 50 positions in the region, we derive optical depths (0.9-5), excitation temperatures of [CII] (80-255 K), and N(C+) of 0.3-1x10^19 cm^-2. The total mass of the globule is ~1000 Msun, about half of which is traced by [CII]. The dense PDR gas has a thermal pressure of 10^7-10^8 K cm^-3, which is similar to the values observed in other regions.


title: Entanglement, and unsorted database search in noise-based logic

abstract: We explore the collapse of "wavefunction" and the measurement of entanglement in the superpositions of hyperspace vectors in classical physical instantaneous-noise-based logic (INBL). We find both similarities with and major differences from the related properties of quantum systems. Two search algorithms utilizing the observed features are introduced. For the first one we assume an unsorted names database set up by Alice that is a superposition (unknown by Bob) of up to n=2^N strings; those we call names. Bob has access to the superposition wave and to the 2N reference noises of the INBL system of N noise bits. For Bob, to decide if a given name x is included in the superposition, once the search has begun, it takes N switching operations followed by a single measurement of the superposition wave. Thus the time and hardware complexity of the search algorithm is O[log(n)] which indicates an exponential speedup compared to Grover's quantum algorithm in a corresponding setting. An extra advantage is that the error probability of the search is zero. Moreover, the scheme can also check the existence of a fraction of a string, or several separate string fractions embedded in an arbitrarily long, arbitrary string. In the second algorithm, we expand the above scheme to a phonebook with n names and s phone numbers. When the names and numbers have the same bit resolution, once the search has begun, the time and hardware complexity of this search algorithm is O[log(n)]. In the case of one-to-one correspondence between names and phone numbers (n=s), the algorithm offers inverse phonebook search too. The error probability of this search algorithm is also zero.


title: Miniature Multi-Level Optical Memristive Switch Using Phase Change Material

abstract: The optical memristive switches are electrically activated optical switches that can memorize the current state. They can be used as optical latching switches in which the switching state is changed only by applying an electrical Write/Erase pulse and maintained without external power supply. We demonstrate an optical memristive switch based on a silicon MMI structure covered with nanoscale-size Ge2Sb2Te5 (GST) material on top. The phase change of GST is triggered by resistive heating of the silicon layer beneath GST with an electrical pulse. Experimental results reveal that the optical transmissivity can be tuned in a controllable and repeatable manner with the maximum transmission contrast exceeding 20 dB. Partial crystallization of GST is obtained by controlling the width and amplitude of the electric pulses. Crucially, we also demonstrate that both Erase and Write operations, to and from any intermediate level, are possible with accurate control of the electrical pulses. Our work marks a significant step forward towards realizing photonic memristive switches without static power consumption, which are highly demanded in high-density large-scale integrated photonics.


title: Stabilized energy factorization approach for Allen-Cahn equation with logarithmic Flory-Huggins potential

abstract: The Allen--Cahn equation is one of fundamental equations of phase-field models, while the logarithmic Flory--Huggins potential is one of the most useful energy potentials in various phase-field models. In this paper, we consider numerical schemes for solving the Allen--Cahn equation with logarithmic Flory--Huggins potential. The main challenge is how to design efficient numerical schemes that preserve the maximum principle and energy dissipation law due to the strong nonlinearity of the energy potential function. We propose a novel energy factorization approach with the stability technique, which is called stabilized energy factorization approach, to deal with the Flory--Huggins potential. One advantage of the proposed approach is that all nonlinear terms can be treated semi-implicitly and the resultant numerical scheme is purely linear and easy to implement. Moreover, the discrete maximum principle and unconditional energy stability of the proposed scheme are rigorously proved using the discrete variational principle. Numerical results are presented to demonstrate the stability and effectiveness of the proposed scheme.


title: Distance-balanced graphs and Travelling Salesman Problems

abstract: For every probability $p\in[0,1]$ we define a distance-based graph property, the $p$TS-distance-balancedness, that in the case $p=0$ coincides with the standard distance-balancedness, and in the case $p=1$ is related to the Hamiltonian-connectedness. In analogy with the classical case, where the distance-balancedness of a graph is equivalent to the property of being self-median, we characterize the class of $p$TS-distance-balanced graphs in terms of their equity with respect to certain probabilistic centrality measures, inspired by the Travelling Salesman Problem. We prove that it is possible to detect this property looking at the classical distance-balancedness (and therefore looking at the classical centrality problems) of a suitable graph composition, namely the wreath product of graphs. More precisely, we characterize the distance-balancedness of a wreath product of two graphs in terms of the $p$TS-distance-balancedness of the factors.


title: Spin-valley Hall transport induced by spontaneous symmetry breaking in half-filled zero Landau level of bilayer graphene

abstract: Intrinsic Hall conductivity, emerging when chiral symmetry is broken, is at the heart of future low energy consumption devices because it can generate non-dissipative charge neutral current. A symmetry breaking state is also induced by electronic correlation even for the centro-symmetric crystalline materials. However, generation of non-dissipative charge neutral current by intrinsic Hall conductivity induced by such spontaneous symmetry breaking is experimentally elusive. Here we report intrinsic Hall conductivity and generation of a non-dissipative charge neutral current in a spontaneous antiferromagnetic state of zero Landau level of bilayer graphene, where spin and valley contrasting Hall conductivity has been theoretically predicted. We performed nonlocal transport experiment and found cubic scaling relationship between the local and nonlocal resistance, as a striking evidence of the intrinsic Hall effect. Observation of such spontaneous Hall transport is a milestone toward understanding the electronic correlation effect on the non-dissipative transport. Our result also paves a way toward electrical generation of a spin current in non-magnetic graphene via coupling of spin and valley in this symmetry breaking state combined with the valley Hall effect.


title: Study of edge states and conductivity in spin-orbit coupled bilayer graphene

abstract: We present an elaborate and systematic study of the conductance properties of a zigzag bilayer graphene nanoribbon modeled by a Kane-Mele (KM) Hamiltonian. The interplay of the Rashba and the intrinsic spin-orbit couplings with the edge states, electronic band structures, charge and spin transport are explored in details. We have analytically derived the conditions for the edge states for a bilayer KM nanoribbon and show how these modes decay for lattice sites inside the bulk. It is particularly interesting to note that for a finite-size ribbon an even number of zigzag ribbon hosts a finite energy gap at the Dirac points, while the odd ones do not. This asymmetry is present both in presence and absence of a bias voltage that may exist between the layers. The interlayer Rashba spin-orbit coupling, along with the intralayer intrinsic spin-orbit and intralayer Rashba spin-orbit couplings seem to destroy the quantum spin Hall (QSH) phase where the QSH phase is identified by the presence of a conductance plateau (of magnitude 4e/h) in the vicinity of zero Fermi energy. The plateau is sensitive to the values of the spin-orbit coupling parameters. Further, the spin polarized conductance data reveal that a bilayer KM ribbon is found to be more efficient for spintronic applications compared to a monolayer graphene. Finally, a quick check with experiments is done via computing the effective mass of electrons.


title: Convolutional Neural Network for Intrusion Detection System In Cyber Physical Systems

abstract: The extensive use of Information and Communication Technology in critical infrastructures such as Industrial Control Systems make them vulnerable to cyber-attacks. One particular class of cyber-attacks is advanced persistent threats where highly skilled attackers can steal user authentication information's and move in the network from host to host until a valuable target is reached. The detection of the attacker should occur as soon as possible in order to take appropriate response, otherwise the attacker will have enough time to reach sensitive assets. When facing intelligent threats, intelligent solutions have to be designed. Therefore, in this paper, we take advantage of recent progress in deep learning to build a convolutional neural networks that can detect intrusions in cyber physical system. The Intrusion Detection System is applied on the NSL-KDD dataset and the performances of the proposed approach are presented and compared with the state of art. Results show the effectiveness of the techniques.


title: The contact structure induced by a line fibration of R^3 is standard

abstract: Building on the work of and answering a question by Michael Harrison, we show that any contact structure on Euclidean 3-space induced by a line fibration is diffeomorphic to the standard contact structure.


title: Surface braid groups, finite Heisenberg covers and double Kodaira fibrations

abstract: We exhibit new examples of double Kodaira fibrations by using finite Galois covers of a product $\Sigma_b \times \Sigma_b$, where $\Sigma_b$ is a smooth projective curve of genus $b \geq 2$. Each cover is obtained by providing an explicit group epimorphism from the pure braid group $\mathsf{P}_2(\Sigma_b)$ to some finite Heisenberg group. In this way, we are able to show that every curve of genus $b$ is the base of a double Kodaira fibration; moreover, the number of pairwise non-isomorphic Kodaira fibred surfaces fibering over a fixed curve $\Sigma_b$ is at least $\boldsymbol{\omega}(b+1)$, where $\boldsymbol{\omega} \colon \mathbb{N} \to \mathbb{N}$ stands for the arithmetic function counting the number of distinct prime factors of a positive integer. As a particular case of our general construction, we obtain a real $4$-manifold of signature $144$ that can be realized as a real surface bundle over a surface of genus $2$, with fibre genus $325$, in two different ways. This provides (to our knowledge) the first "double" solution to a problem from Kirby's list in low-dimensional topology.


title: Retrieval of the fluid Love number $k_2$ in exoplanetary transit curves

abstract: We are witness to a great and increasing interest in internal structure, composition and evolution of exoplanets. However, direct measurements of exoplanetary mass and radius cannot be uniquely interpreted in terms of interior structure, justifying the need for an additional observable. The second degree fluid Love number, $k_2$, is proportional to the concentration of mass towards the body's center, hence providing valuable additional information about the internal structure. When hydrostatic equilibrium is assumed for the planetary interior, $k_2$ is a direct function of the planetary shape. Previous attempts were made to link the observed tidally and rotationally induced planetary oblateness in photometric light curves to $k_2$ using ellipsoidal shape models. Here, we construct an analytical 3D shape model function of the true planetary mean radius, that properly accounts for tidal and rotational deformations. Measuring the true planetary mean radius is critical when one wishes to compare the measured $k_2$ to interior theoretical expectations. We illustrate the feasibility of our method and show, by applying a Differential Evolution Markov Chain to synthetic data of WASP-121b, that a precision $\leq$ 65 ppm/$\sqrt{2\,min}$ is required to reliably retrieve $k_2$ with present understanding of stellar limb darkening, therefore improving recent results based on ellipsoidal shape models. Any improvement on stellar limb darkening would increase such performance.


title: Parameters Calibration for Power Grid Stability Models using Deep Learning Methods

abstract: This paper presents a novel parameter calibration approach for power system stability models using automatic data generation and advanced deep learning technology. A PMU-measurement-based event playback approach is used to identify potential inaccurate parameters and automatically generate extensive simulation data, which are used for training a convolutional neural network (CNN). The accurate parameters will be predicted by the well-trained CNN model and validated by original PMU measurements. The accuracy and effectiveness of the proposed deep learning approach have been validated through extensive simulation and field data.


title: The Gerrymandering Jumble: Map Projections Permute Districts' Compactness Scores

abstract: In political redistricting, the compactness of a district is used as a quantitative proxy for its fairness. Several well-established, yet competing, notions of geographic compactness are commonly used to evaluate the shapes of regions, including the Polsby-Popper score, the convex hull score, and the Reock score, and these scores are used to compare two or more districts or plans. In this paper, we prove mathematically that any map projection from the sphere to the plane reverses the ordering of the scores of some pair of regions for all three of these scores. Empirically, we demonstrate that the effect of using the Cartesian latitude-longitude projection on the order of Reock scores is quite dramatic.


title: Index of minimal spheres and isoperimetric eigenvalue inequalities

abstract: In the present paper we use twistor theory in order to solve two problems related to harmonic maps from surfaces to Euclidean spheres $\mathbb{S}^n$. First, we propose a new approach to isoperimetric inequalities based on energy index. Using this approach we show that for any positive $k$, the $k$-th non-zero eigenvalue of the Laplacian on the real projective plane endowed with a metric of unit area, is maximized on the sequence of metrics converging to a union of $(k-1)$ identical copies of round sphere and a single round projective plane. This extends the results of P. Li and S.-T. Yau for $k=1$ (1982); N. Nadirashvili and A. Penskoi for $k=2$ (2018); and confirms the conjecture made in [KNPP]. Second, we improve the known upper bounds for the area index of minimal two-dimensional spheres and minimal projective planes in $\mathbb{S}^n$. In the course of the proof we establish a twistor correspondence for Jacobi fields, which could be of independent interest for the study of moduli space of harmonic maps.


title: A Hardware-Oriented and Memory-Efficient Method for CTC Decoding

abstract: The Connectionist Temporal Classification (CTC) has achieved great success in sequence to sequence analysis tasks such as automatic speech recognition (ASR) and scene text recognition (STR). These applications can use the CTC objective function to train the recurrent neural networks (RNNs), and decode the outputs of RNNs during inference. While hardware architectures for RNNs have been studied, hardware-based CTCdecoders are desired for high-speed CTC-based inference systems. This paper, for the first time, provides a low-complexity and memory-efficient approach to build a CTC-decoder based on the beam search decoding. Firstly, we improve the beam search decoding algorithm to save the storage space. Secondly, we compress a dictionary (reduced from 26.02MB to 1.12MB) and use it as the language model. Meanwhile searching this dictionary is trivial. Finally, a fixed-point CTC-decoder for an English ASR and an STR task using the proposed method is implemented with C++ language. It is shown that the proposed method has little precision loss compared with its floating-point counterpart. Our experiments demonstrate the compression ratio of the storage required by the proposed beam search decoding algorithm are 29.49 (ASR) and 17.95 (STR).


title: Multi-target Detection with an Arbitrary Spacing Distribution

abstract: Motivated by the structure reconstruction problem in single-particle cryo-electron microscopy, we consider the multi-target detection model, where multiple copies of a target signal occur at unknown locations in a long measurement, further corrupted by additive Gaussian noise. At low noise levels, one can easily detect the signal occurrences and estimate the signal by averaging. However, in the presence of high noise, which is the focus of this paper, detection is impossible. Here, we propose two approaches---autocorrelation analysis and an approximate expectation maximization algorithm---to reconstruct the signal without the need to detect signal occurrences in the measurement. In particular, our methods apply to an arbitrary spacing distribution of signal occurrences. We demonstrate reconstructions with synthetic data and empirically show that the sample complexity of both methods scales as 1/SNR^3 in the low SNR regime.


title: Does Data Augmentation Lead to Positive Margin?

abstract: Data augmentation (DA) is commonly used during model training, as it significantly improves test error and model robustness. DA artificially expands the training set by applying random noise, rotations, crops, or even adversarial perturbations to the input data. Although DA is widely used, its capacity to provably improve robustness is not fully understood. In this work, we analyze the robustness that DA begets by quantifying the margin that DA enforces on empirical risk minimizers. We first focus on linear separators, and then a class of nonlinear models whose labeling is constant within small convex hulls of data points. We present lower bounds on the number of augmented data points required for non-zero margin, and show that commonly used DA techniques may only introduce significant margin after adding exponentially many points to the data set.


title: Time-Frequency Analysis in the Unit Ball

abstract: As an appropriate analog of the Euclidean short-time Fourier transform, we study a windowed version of the Helgason-Fourier transform on the complex unit ball and translate the theory of modulation/coorbit spaces. As a result, atomic decompositions by means of Banach frames are obtained.


title: Anytime Multi-arm Task and Motion Planning for Pick-and-Place of Individual Objects via Handoffs

abstract: Automation applications are pushing the deployment of many high DoF manipulators in warehouse and manufacturing environments. This has motivated many efforts on optimizing manipulation tasks involving a single arm. Coordinating multiple arms for manipulation, however, introduces additional computational challenges arising from the increased DoFs, as well as the combinatorial increase in the available operations that many manipulators can perform, including handoffs between arms. The focus here is on the case of pick-and-place tasks, which require a sequence of handoffs to be executed, so as to achieve computational efficiency, asymptotic optimality and practical anytime performance. The paper leverages recent advances in multi-robot motion planning for high DoF systems to propose a novel multi-modal extension of the dRRT* algorithm. The key insight is that, instead of naively solving a sequence of motion planning problems, it is computationally advantageous to directly explore the composite space of the integrated multi-arm task and motion planning problem, given input sets of possible pick and handoff configurations. Asymptotic optimality guarantees are possible by sampling additional picks and handoffs over time. The evaluation shows that the approach finds initial solutions fast and improves their quality over time. It also succeeds in finding solutions to harder problem instances relative to alternatives and can scale effectively as the number of robots increases.


title: Molecules with dipoles in periodic boundary conditions in a tetragonal cell

abstract: When a system which contains a dipole, and whose dimensionality is less than three, is studied in a code which imposes periodic boundary conditions in all three dimensions, an artificial electric field arises which keeps the potential periodic. This has an impact on the total energy of the system, and on any other attribute which would respond to an electric field. Simple corrections are known for 0D systems embedded in a cubic geometry, and 2D slab systems. This paper shows how the 0D result can be extended to tetragonal geometries, and that for a particular $c/a$ ratio the correction is zero. It also considers an exponential error term absent from the usual consideration of 2D slab geometries, and discusses an empirical form for this.


title: Precision measurement of the proton charge radius in electron proton scattering

abstract: This report presents a project of an experiment for precision studies of the elastic electron proton scattering in the low momentum transfer region. The project is based on an innovative experimental method which allows for detection of recoil protons in a hydrogen "active target". The goal of the experiment is to measure the $ep$ elastic scattering differential cross section in the $Q^2$ range from 0.001 GeV$^2$ to 0.04 GeV$^2$ with 0.1% relative and 0.2% absolute precision and to determine the proton charge radius with a sub-percent precision. The important advantages of the proposed experiment are low radiative corrections, inherent to the recoil proton method, and absolute measurements of the differential cross sections. The experiment will be performed in the 720 MeV electron beam of the Mainz electron accelerator MAMI. This accelerator can provide an electron beam with optimal for this experiment parameters. The Proposal was approved by the MAMI Program Advisory Committee, and a special Agreement aimed at realization of this experiment was signed between Petersburg Nuclear Physics Institute and Institute of Nuclear Physics of the Mainz University. Also, scientists from the following institutes will participate in this project: GSI Centre for Heavy Ion Research, Germany; Joint Institute for Nuclear Research, Russia; William and Mary College, USA; Mount Allison University, Canada; University of Regina, Canada; Saint Mary's University, Canada. The first measurements should be started in 2020.


title: Death by starvation in May-Leonard models

abstract: We consider the dynamics of spatial stochastic May-Leonard models with mutual predation interactions of equal strength between any two individuals of different species. Using two-dimensional simulations, with two and three pecies, we investigate the dynamical impact of the death of individuals after a given threshold number of successive unsuccessful predation attempts. We find that the death of these individuals can have a strong impact on the dynamics of population networks and provide a crucial contribution to the preservation of coexistence.


title: Reconstructing Classes of Non-bandlimited Signals from Time Encoded Information

abstract: We investigate time encoding as an alternative method to classical sampling, and address the problem of reconstructing classes of non-bandlimited signals from time-based samples. We consider a sampling mechanism based on first filtering the input, before obtaining the timing information using a time encoding machine. Within this framework, we show that sampling by timing is equivalent to a non-uniform sampling problem, where the reconstruction of the input depends on the characteristics of the filter and on its non-uniform shifts. The classes of filters we focus on are exponential and polynomial splines, and we show that their fundamental properties are locally preserved in the context of non-uniform sampling. Leveraging these properties, we then derive sufficient conditions and propose novel algorithms for perfect reconstruction of classes of non-bandlimited signals such as: streams of Diracs, sequences of pulses and piecewise constant signals. Next, we extend these methods to operate with arbitrary filters, and also present simulation results on synthetic noisy data.


title: Implementing Efficient Message Logging Protocols as MPI Application Extensions

abstract: Message logging protocols are enablers of local rollback, a more efficient alternative to global rollback, for fault tolerant MPI applications. Until now, message logging MPI implementations have incurred the overheads of a redesign and redeployment of an MPI library, as well as continued performance penalties across various kernels. Successful research efforts for message logging implementations do exist, but not a single one of them can be easily deployed today by more than a few experts. In contrast, in this work we build efficient message logging capabilities on top of an MPI library with no message logging capabilities; we do so for two different HPC kernels, one with a global exchange pattern (CG), and one with a neighbourhood exchange pattern (LULESH). While our library of choice ULFM detects failure and recovers MPI communicators, we build on that to then restore the intra- and inter-process data consistency of both applications. This task turns out to be challenging, and we present the methodology for doing so in this work. In the end, we achieve message logging capabilities for each kernel, without the need for an actual message logging runtime underneath. On the performance side, we match state-of-the-art solutions and (a) eliminate event logging and the event logger component altogether, and (b) design a hybrid protocol, which gracefully shifts between global and local rollback, depending on the available payload logging memory. Such a hybrid protocol between local and global rollback has not been previously proposed to our knowledge. Our extensions span a few hundred lines of code for each kernel, are open-sourced, and enable local and global rollback after process failure.


title: Phase oscillator model for noisy oscillators

abstract: The Kuramoto model has become a paradigm to describe the dynamics of nonlinear oscillator under the influence of external perturbations, both deterministic and stochastic. It is based on the idea to describe the oscillator dynamics by a scalar differential equation, that defines the time evolution for the phase of the oscillator. Starting from a phase and amplitude description of noisy oscillators, we discuss the reduction to a phase oscillator model, analogous to the Kuramoto model. The model derived shows that the phase noise problem is a drift-diffusion process. Even in the case where the expected amplitude remains unchanged, the unavoidable amplitude fluctuations do change the expected frequency, and the frequency shift depends on the amplitude variance. We discuss different degrees of approximation, yielding increasingly accurate phase reduced descriptions of noisy oscillators.


title: A mathematical framework for amplitude and phase noise analysis of coupled oscillators

abstract: Synchronization of coupled oscillators is a paradigm for complexity in many areas of science and engineering. Any realistic network model should include noise effects. We present a description in terms of phase and amplitude deviation for nonlinear oscillators coupled together through noisy interactions. In particular, the coupling is assumed to be modulated by white Gaussian noise. The equations derived for the amplitude deviation and the phase are rigorous, and their validity is not limited to the weak noise limit. We show that using Floquet theory, a partial decoupling between the amplitude and the phase is obtained. The decoupling can be exploited to describe the oscillator's dynamics solely by the phase variable. We discuss to what extent the reduced model is appropriate and some implications on the role of noise on the frequency of the oscillators.


title: Path-following methods for calculating linear surface wave dispersion relations on vertical shear flows

abstract: The path-following scheme in [Loisel and Maxwell, SIAM J. Matrix Anal. Appl., 39-4 (2018), pp. 1726-1749] is adapted to efficiently calculate the dispersion relation curve for linear surface waves on an arbitrary vertical shear current. This is equivalent to solving the Rayleigh instability equation with linearised free-surface boundary condition for each sought point on the curve. Taking advantage of the analyticity of the dispersion relation, a path-following or continuation approach is adopted. The problem is discretized using a collocation scheme, parametrised along either a radial or angular path in the wave vector plane, and differentiated to yield a system of ODEs. After an initial eigenproblem solve using QZ decomposition, numerical integration proceeds along the curve using linear solves as the Runge--Kutta $F(\cdot)$ function; thus, many QZ decompositions on a $2N$ companion matrix are exchanged for one QZ decomposition and a small number of linear solves on a size $N$ matrix. A piecewise interpolant provides dense output. The integration represents a nominal setup cost afterwhich very many points can be computed at negligible cost whilst preserving high accuracy. Furthermore, a 2-dimensional interpolant suitable for scattered data query points in the wave vector plane is described. Finally, a comparison is made with existing numerical methods for this problem, revealing that the path-following scheme is asymptotically two orders of magnitude faster in number of query points.


title: Model-Based Reconstruction for Quantitative MRI using the Bloch Equations

abstract: In this work a generic model-based reconstruction for the quantification of relaxation parameters is developed. In contrast to previous approaches that rely on simplified models derived from the Bloch equations, this work includes the Bloch equations directly into the reconstruction. Therefore, non-linear calibrationless parallel imaging is combined with a generic Runge-Kutta 5(4) based forward operator to simulate spin dynamics described by arbitrary sequences. Gradients are determined by using a sensitivity analysis and solving the resulting ordinary differential equations parallel to the signal simulation. Based on this formulation an IRGNM-FISTA algorithm is used to estimate quantitative maps for $T_1$, $T_2$, $M_0$, and the coil profiles $c_N$ from fully-sampled multi-inversion and golden-angle single-shot inversion-recovery radial bSSFP measurement with a custom-built $T_1$-$T_2$ phantom.


title: Online reviews can predict long-term returns of individual stocks

abstract: Online reviews are feedback voluntarily posted by consumers about their consumption experiences. This feedback indicates customer attitudes such as affection, awareness and faith towards a brand or a firm and demonstrates inherent connections with a company's future sales, cash flow and stock pricing. However, the predicting power of online reviews for long-term returns on stocks, especially at the individual level, has received little research attention, making a comprehensive exploration necessary to resolve existing debates. In this paper, which is based exclusively on online reviews, a methodology framework for predicting long-term returns of individual stocks with competent performance is established. Specifically, 6,246 features of 13 categories inferred from more than 18 million product reviews are selected to build the prediction models. With the best classifier selected from cross-validation tests, a satisfactory increase in accuracy, 13.94%, was achieved compared to the cutting-edge solution with 10 technical indicators being features, representing an 18.28% improvement relative to the random value. The robustness of our model is further evaluated and testified in realistic scenarios. It is thus confirmed for the first time that long-term returns of individual stocks can be predicted by online reviews. This study provides new opportunities for investors with respect to long-term investments in individual stocks.


title: Convex choice, finite choice and sorting

abstract: We study the Weihrauch degrees of closed choice for finite sets, closed choice for convex sets and sorting infinite sequences over finite alphabets. Our main results are: One, that choice for finite sets of cardinality $i + 1$ is reducible to choice for convex sets in dimension $j$, which in turn is reducible to sorting infinite sequences over an alphabet of size $k + 1$, iff $i \leq j \leq k$. Two, that convex choice in dimension two is not reducible to the product of convex choice in dimension one with itself. Three, that sequential composition of one-dimensional convex choice is not reducible to convex choice in any dimension. The latter solves an open question raised at a Dagstuhl seminar on Weihrauch reducibility in 2015. Our proofs invoke Kleene's recursion theorem, and we describe in some detail how Kleene's recursion theorem gives rise to a technique for proving separations of Weihrauch degrees.


title: On the asymptotic Plateau problem for area minimizing surfaces in $\mathbb{E}(-1,\tau)$

abstract: We prove some existence and non-existence results for complete area minimizing surfaces in the homogeneous space $\mathbb{E}(-1,\tau)$. As one of our main results, we present sufficient conditions for a curve $\Gamma$ in $\partial_{\infty} \mathbb{E}(-1,\tau)$ to admit a solution to the asymptotic Plateau problem, in the sense that there exists a complete area minimizing surface in $\mathbb{E}(-1,\tau)$ having $\Gamma$ as its asymptotic boundary.


title: A hybrid recommendation algorithm based on weighted stochastic block model

abstract: Hybrid recommendation usually combines collaborative filtering with content-based filtering to exploit merits of both techniques. It is widely accepted that hybrid filtering outperforms the single algorithm, thus it has been the new trend in electronic commerce these years. In this paper, we propose a novel hybrid recommendation system based on weighted stochastic block model (WSBM). Our algorithm not only makes full use of content-based and collaborative filtering recommendation to solve the cold-start problem but also improves the accuracy of recommendation by selecting the nearest neighbor with WSBM. The experiment result shows that our proposed approach has better prediction and classification accuracy than traditional hybrid recommendation.


title: Blockchain-enabled Authentication Handover with Efficient Privacy Protection in SDN-based 5G Networks

abstract: 5G mobile networks provide additional benefits in terms of lower latency, higher data rates, and more coverage, in comparison to 4G networks, and they are also coming close to standardization. For example, 5G has a new level of data transfer and processing speed that assures users are not disconnected when they move from one cell to another; thus, supporting faster connection. However, it comes with its own technical challenges relating to resource management, authentication handover and user privacy protection. In 5G, the frequent displacement of the users among the cells as a result of repeated authentication handovers often lead to a delay, contradicting the 5G objectives. In this paper, we propose a new authentication approach that utilizes blockchain and software defined networking (SDN) techniques to remove the re-authentication in repeated handover among heterogeneous cells. The proposed approach is designed to assure the low delay, appropriate for the 5G network in which users can be replaced with the least delay among heterogeneous cells using their public and private keys provided by the devised blockchain component while protecting their privacy. In our comparison between Proof-of-Work (POW)-based and network-based models, the delay of our authentication handover was shown to be less than 1ms. Also, our approach demonstrated less signaling overhead and energy consumption compared to peer models.


title: Zeros and approximations of Holant polynomials on the complex plane

abstract: We present fully polynomial approximation schemes for a broad class of Holant problems with complex edge weights, which we call Holant polynomials. We transform these problems into partition functions of abstract combinatorial structures known as polymers in statistical physics. Our method involves establishing zero-free regions for the partition functions of polymer models and using the most significant terms of the cluster expansion to approximate them. Results of our technique include new approximation and sampling algorithms for a diverse class of Holant polynomials in the low-temperature regime and approximation algorithms for general Holant problems with small signature weights. Additionally, we give randomised approximation and sampling algorithms with faster running times for more restrictive classes. Finally, we improve the known zero-free regions for a perfect matching polynomial.


title: Four-point conformal blocks with three heavy background operators

abstract: We study CFT$_2$ Virasoro conformal blocks of the 4-point correlation function $\langle \mathcal{O}_L \mathcal{O}_H \mathcal{O}_H \mathcal{O}_H \rangle $ with three background operators $\mathcal{O}_H$ and one perturbative operator $\mathcal{O}_L$ of dimensions $\Delta_L/\Delta_H \ll1$. The conformal block function is calculated in the large central charge limit using the monodromy method. From the holographic perspective, the background operators create $AdS_3$ space with three conical singularities parameterized by dimensions $\Delta_H$, while the perturbative operator corresponds to the geodesic line stretched from the boundary to the bulk. The geodesic length calculates the perturbative conformal block. We propose how to address the block/length correspondence problem in the general case of higher-point correlation functions $\langle \mathcal{O}_L \cdots \mathcal{O}_L \mathcal{O}_H \cdots \mathcal{O}_H \rangle $ with arbitrary numbers of background and perturbative operators.


title: Discovering and Proving Invariants in Answer Set Programming and Planning

abstract: Answer set programming (ASP) and planning are two widely used paradigms for solving logic programs with declarative programming. In both cases, the quality of the input programs has a major influence on the quality and performance of the solving or planning process. Hence, programmers need to understand how to make their programs efficient and still correct. In my PhD studies, I explore how input programs can be improved and verified automatically as a means to support programmers. One of my research directions consists in discovering invariants in planning programs without human support, which I implemented in a system called ginkgo. Studying dynamic systems in greater depth, I then developed plasp 3 with members of my research group, which is a significant step forward in effective planning in ASP. As a second research direction, I am concerned with automating the verification of ASP programs against formal specifications. For this joint work with Lifschitz's group at the University of Texas, I developed a verification system called anthem. In my future PhD studies, I will extend my research concerning the discovery and verification of ASP and planning problems.